传统与旋转交错网格有限差分在双相介质中模拟对比.docxVIP

传统与旋转交错网格有限差分在双相介质中模拟对比林朋，卢勇旭（中国矿业大学（北京）煤炭资源与安全开采国家重点实验室，北京 100083）摘要：以Biot双相介质模型为背景，推导了双相各向同性介质二维三分量一阶速度—应力弹性波方程方程，建立了各向同性双相介质波动方程的二维三分量有限差分格式。分别采用传统交错网格有限差分技术与旋转交错网格有限差分技术对均匀与非均匀双相各向同性介质进行了波场模拟。结果表明：旋转交错网格有限差分技术能够有效模拟双相各向同性介质中弹性波的传播情况；通过传统与旋转交错网格有限差分技术的对比，说明了旋转交错网格有限差分算法在稳定性更强，避免了插值带来的误差，是一种有效的地震波场模拟方法。关键词：双相介质；传统交错网格；旋转交错网格；波场模拟随着人们对地球内部认识的逐步加深与地球物理学的快速发展，含饱与流体的孔隙介质被认为是与含油气储层最为接近的介质模型，其是由固体骨架颗粒与孔隙中的流体（如气、水等）组成双相或多相介质。双相介质中弹性波传播规律的正确认识，对于油气的开采具有重要意义。地震波数值模拟是人们用来描述与认识地震波传播规律的有效途径。交错网格有限差分技术作为地震波场模拟的常用方法之一，在地球物理中一直被广范应用。Biot理论描述了饱与流体空隙介质中地震波的传播，是双相介质波动理论的基础[1-3]；Schmitt讨论了柱坐标系下地震波在横向各向异性介质中的传播问题[4]；Crampin通过一系列实验研究发现，双相各向异性介质中存在横波分裂现象[5]；王尚旭研究了双相介质地震波传播规律，并利用有限元法实现了双相介质地震波场模拟[6]；牟永光应用有限差分技术对孔隙各向同性介质进行了波场分析[7]；刘洋等通过虚谱法对双相各向异性介质中弹性波的传播特征[8]；王秀明等使用高阶交错网格有限差分技术实现了非均匀孔隙介质的正演模拟[9]；裴正林通过交错网格有限差分法实现了双相各向异性介质与三维横向各向同性介质弹性波的高阶波场模拟[10-11]。本文使用传统与旋转交错网格两种有限差分技术对双相介质进行了波场模拟，通过对比，说明了旋转交错网格有限差分技术的有效性与优越性。1 方法原理双相各向异性介质二维三分量弹性波方程根据双相介质地震波传播理论[2]，可得到饱与流体孔隙介质运动方程[11]（1）（2）式中，i、j表示x,y,z三个不同方向分量，bij表示耗散系数，为固相应力分量，s为作用在流体上的有效应力，为介质固体骨架速度矢量,为介质流体速度矢量。令为孔隙度，与分别表示单位体积内固体骨架与流体部分的有效质量，为流体相对固体骨架运动时的视质量[12]。三者与固相密度与流相密度之间满足[2]（3）（4）由饱与流体孔隙介质传播理论，易得双相各向同性介质二维三分量一阶速度—应力表达式。令Qi表示固体骨架与流体空隙之间体积变化的耦合参数，R表示描述流体的弹性参数，即将一定体积流体注入孔隙介质体积元，为保持总体积不变时在流体上施加的一种力的度量[12]。对于双相各向同性介质，固相部分速度与应力之间的关系为（5）由广义达西定律可得Biot介质流相部分速度与应力关系式（6）交错网格有限差分技术交错网格有限差分技术是一种有效的地震波正演模拟方法。传统交错网格是将网格剖分成整网格点与半网格点，在相邻的两个时间层上的网格点处分别定义速度与应力分量，并且在空间分布上相邻的两个时间层上的物理量恰好交错半个网格，导数值在半网格点处计算，以实现时间与空间的交错[13]。由于速度与应力的相对关系，对于不同的介质模型，在计算过程中需要对部分场量与模型参数进行插值，增加了计算误差，降低了计算精度。网格定义、波场分量及弹性参数位置如图1（a）所示。旋转交错网格由传统交错网格发展而来，通过旋转对网格进行了重新划分与定义。与传统交错网格相比，不同之处在于，在同一网格点处仅定义同一物理量（速度、应力），通过计算沿网格对角线物理量的差分来计算微分，由于速度分量与应力分量各定义整网格点与半网格点上，避免了部分场量与模型参数的插值[14]，降低了计算误差，提高了计算精度[15-16]。网格定义、波场分量及弹性参数位置如图1（b）所示[17]。（a）传统交错网格（b）旋转交错网格图1 交错网格示意图稳定性分析稳定性问题是地球物理学科研究中所必须解决的难题，直接关系到数值模拟方法的成败。对于传统交错网格有限差分技术，在步长相等的情况下，时间域二阶、空间2M阶的稳定性条件如下：（7）式中，是时间步长，D是空间维数，为空间步长，Vmax为最大相速度，Ck时是空间差分系数。对于旋转交错网格有限差分技术，Saenger[14]在Neumann稳定性条件下，给出了在空间步长相等时，时间域二阶、空间2M阶的稳定性条件为：（8）式中参数意义同（7）式。通过对比两种交错网格的稳定性条件可