震电效应下平面波在流体和孔隙介质界面上的折反射

【摘要】：对于由弹性波场激发所引起的电磁场变化效应震电效应的研究,正成为勘探地球物理领域,特别是油气勘探开发业界的关注热点。震电效应应用于油气勘探开发,特别是在寻找剩余油方面前景十分诱人,其产生的电磁波场对地下油气藏分布也具有直接探测性。由此,本文旨在研究震电效应下,平面声波在流体与孔隙介质分界面上的折反射,围绕这个主题开展了研究。 本文对研究震电所需的理论基础——Biot理论进行了简单介绍,Biot建立了流体饱和孔隙介质弹性波动力学理论,该理论对研究的孔隙介质作了一些假设后,得到了应力应变关系和运动方程,给出了均匀孔隙介质中的平面波解,并且对孔隙介质中是否存在慢纵波进行了预测。 本文简单介绍了震电效应产生的机理,弹性波在孔隙介质中传播时,填充在孔隙中的流体会相对于固相骨架运动,形成波动渗流。若孔隙中的流体含净剩电荷,这种波动渗流则会引起电磁场。相反,电磁场可引起渗流和固相骨架运动,形成弹性波。这种叫做动电耦合现象,在地震频率范围称为震电效应。 Pride建立了声场和电磁场的耦合方程组,本文通过对电场矢量和位移矢量进行亥姆霍兹分解,把用电场矢量和位移矢量表示的Pride声电耦合方程组转化为关于势函数的波动方程组,在平面直角坐标中导出在流体饱和空隙介质中声电耦合场的形式解,通过声电耦合问题的边界条件,独立编制了计算程序,计算了能流折反射系数,计算表明边界上法向能流连续,证明程序的正确性。通过数值模拟,可以获得折反射系数随不同入射角和不同频率的变化规律,计算结果表明电磁波能量反射和透射系数都是在某一频率值达到最大后,随频率的增大而逐渐减小;流体与孔隙介质界面不同的边界状况对能量折反射系数也有明显的影响。折反射系数不仅与入射角度和频率有关,也和地层的其他参数有关,本文模型和编制的程序为进一步详细深入研究它们的关系提供理论基础和计算软件。 此外,为了提高计算效率,我们先独立求解了声场,分别采用非耦合和似稳法进行研究。非耦合法则是在忽略声波诱导电磁场对声场影响的前提下,将Pride声电耦合方程组简化为求解具有波动电流源的麦克斯韦方程组;然后引入赫兹矢量,把求解具有波动电流源的麦克斯韦方程组问题转化为求解关于赫兹矢量的非齐次矢量亥姆霍兹方程组;运用标量亥姆霍兹方程的解,找出非齐次矢量亥姆霍兹方程组的一组特解,与齐次亥姆霍兹方程组的通解叠加,得出电场、磁场的表达式,由于非耦合法计算上也相对复杂,本文只对其进行了理论研究。似稳法同样是忽略声波诱导电磁场对声场影响的前提下,把低频电磁场看成似稳场,引入电势后,利用了边界条件求出电场,再利用电场进一步确定电流密度和磁场强度。通过独立编程初步试算表明,简化的算法编制的程序相对简单,在低于1千赫兹频率以下与严格解求得的电场反射系数比相对误差在10-5,说明在地震频率范围内完全可以采用似稳法提高计算效率。 本文的研究结果表明:平面分层介质声电耦合波折反射系数是与入射角度、频率及地层参数有关的;采用完成的Pride动电耦合方程得到的折反射系数公式体系和编制的程序为其它简化方法提供严格比较标准;似稳法简单而且计算效率高在地震和测井频率范围值得采用。本文的理论研究结果、编制的程序对水平夹层介质中声电耦合波的特性等其他声电耦合波场的课题提供了参考价值。