浸液流体饱与孔隙柱腔中的震电耦合波.pdfVIP

中国地球物理2011 ·655· (18)储层地球物理 浸液流体饱和孔隙柱腔中的震电耦合波 wave in in Seismoelectric shellfluidmedia cylindrical propagationporous 高文阳+张言崔志文王克协 Gao YanCuiZhiwen Kexie Wenyang+Zhang Wang 吉林大学物理学院长春130012 当弹性波通过流体饱和孔隙介质时，引起固体骨架与孔隙中的电解质溶液的相对运动，就会诱导伴随 的电磁场，孔隙介质中的这种弹性波与电磁场的耦合称为声电效应。在地震频率范围内称为震电效应。近 些年来，震电效应倍受关注，理论和实验研究取得了长足的进展。本文考虑流体饱和孔隙柱腔浸在流体内， 研究震电耦合波在流体饱和孔隙柱腔中的传播问题，对震电耦合波在流体孔隙柱腔中的传播进行数值模拟 分析。以柱腔的中心轴为Z轴建立柱坐标系，声源换能器放在Z轴上。从描述流体饱和孔隙介质中弹性波 和电磁场耦合现象的Pride方程组出发，通过引入势函数采用亥姆赫兹分解，推导震电耦合波在孔隙柱腔 中耦合场的理论公式，得出孔隙腔体中势函数表示的声场和电磁场表达式。同时也采用似稳法来解决该问 题。对腔体内外声、电磁场、声电耦合波场进行数值模拟，并比较两种方法所得的结果相互验证。 理Pride方程组可得以位移U、渗流位移W、电场E为基本场量的孔隙介质声电耦合方程组。对这三个矢量 进行亥姆霍兹分解，在轴对称条件下，可得亥姆霍兹方程势函数的通解，进而求得位移“、渗流位移w、 电场E。通过Pride方程组中其它场量与这三个基本场量的关系，即可求得声压、应力、磁场的计算公式。 在单相介质中声场和电场是同时存在，但二者是脱耦的。同样，对场量进行亥姆霍兹分解，若加上位于井 轴上的点声源的辐射场，可得在腔内和腔外的流体中声场和电场表达式。对于孔内的声场为有源场，其余 为无源场。将场量的表达式带入流体和孔隙柱腔的内外边界条件，可得出场量中的未知系数。利用这些计 算公式可以模拟孔内外声场和电磁场。为使问题简化，采用胡恒山等(2003)提出的似稳法考虑此问题。作 如下假设：第～，忽略电场对声场得影响，使声场不再受电场的影响，可独立求解。第二，把电磁场看成 似稳场，则电磁场的时间导数项可以忽略。在似稳条件下电场的散度为零，则电场可以表示为电势的负梯 导率，妒为电势，L(∞)为动电耦合系数，P为流体压强，∥为流体密度，U为固相位移。方程的左边是关于 电势的，右边是关于声场的，显然是声场作为了电场的源来诱导电场。这样把求解电场的问题转化为求满 足腔中电势方程和一定边界条件的电势妒。由腔中电势的方程可得电势的解，包括两项：一项是齐次方程 在柱坐标系下的通解，另一项是非齐次方程的特解。通过电场与电势，电流与电势的关系可得到电场和电 流的表达式。同样，在似稳条件下可得腔内和腔外的电势，通过电场和电流与电势的关系可得腔内和腔外 的电场和电流的表达式。通过边界上的电势和径向电流连续，确定腔内流体、腔外流体、腔体中的电场。 通过计算考察了不同腔厚(柱腔厚度为孔隙柱腔内外半径的差)时腔内流体中的声场和电磁场，数值模 拟得到了腔内流体中转换声场和电场、磁场的全波波形。考查了壳厚比、孔隙度和渗透率对转换电磁场的 影响。对通常的模型参数，基于pride方程组所得腔内流体电场与似稳法所得到腔内流体的电场进行数值 模拟对比表明，基本看不出两者的所得结果的差异。从而验证了似稳法的正确性，也说明可以利用似稳法 来解决流体饱和孔隙柱壳问题以提高计算效率。与此同时还详细地考察了波场的频率波数域响应特点。本 文的研究为理论和实验研究震电耦合波在浸液流体饱和孔隙柱腔中传播特性提供了参考依据。 本研究由国家自然科学基金和声场声信息国家重点实验室开放课题研究基金资助。