基于交错网格的纵横波波场分离数值模拟方法-工业技术创新.PDF

E X P L O R A T I O N T E C H N I Q U E S 勘探技术 DOI ：10.3969/j .issn.1672-7703.2010.03.008 基于交错网格的纵横波波场分离数值模拟方法 陈可洋 （中国石油大庆油田公司勘探开发研究院，黑龙江省大庆市 163712） 摘 要：为了更精确地研究混合波场中纯纵波和纯横波波场的传播规律，常采用算子分裂方法构造等效的一阶双 曲型弹性波动方程，实现弹性波波场分离数值模拟，这不但增加了计算量和额外的计算时间，而且还增加了计算的复杂 度。提出了基于交错网格的二阶弹性波动方程波场分离数值模拟的快速实现方法，对层状介质模型和固液界面模型进行 模拟试验，并根据交错网格的特点计算出相应的纯纵波 （散度场）和纯横波 （旋度场），从而克服了上述困难。数值计 算结果表明，该方法能够从混合波场中准确地分离出纯纵波和纯横波波场，计算精度和效率得到有效提高，同时指出波 场分离数值模拟方法不适用于固液界面等情况。 关键词：高阶交错网格；波场分离；弹性波动方程；散度和旋度；VSP 中图分类号：P631.41文献标识码：A 弹性波数值模拟技术一直是勘探地球物理界最为 本文提出了基于交错网格的二阶弹性波动方程波 活跃的研究内容之一。利用精确的弹性波动方程数值 场分离数值模拟的快速实现方法，不必进行算子分裂 解模拟复杂的地下弹性波场，为研究弹性波传播机理、 或化为一阶双曲型方程，仅需在每一个时间步中，根据 多波多分量资料的特殊处理方法，以及复杂地层的岩 交错网格来求取散度场（即纯纵波场）与旋度场（即纯 性、构造解释等许多方面提供更为科学的数学物理依 横波场），它们的计算过程又独立于二阶弹性波动方程 据。因此弹性波数值模拟方法已成为认识弹性波传播 高阶交错网格数值模拟方法的整个计算过程，因此该 规律和指导野外多波勘探的一种重要手段，其中研究 方法对整个弹性波数值模拟的精度和稳定性没有任何 纯纵波和纯横波波场的弹性波波场分离数值模拟方法 影响，同时通过应用实例研究了本文方法的适用条件。 正逐步兴起，开始作为弹性数值模拟理论和反演问题 该方法以期作为研究弹性波场传播规律的一种快速有 [1] 新的研究方向。马德堂等 首次提出可实现纯纵波和纯 效方法。 横波波场分离的等价方程，采用虚谱法求解该等价的 1 波场分离数值模拟原理 二阶弹性波动方程，实现了波场分离数值模拟，获得较 好的数值模拟结果，但是伪谱法边界处理困难，计算量 1.1 二阶弹性波动方程及其高阶交错网格有限差分格 [2,3] [4] 大，使得该方法的推广应用受到限制 。李振春等 和 式的构造 [5] Zhang Jianlei 等 提出基于一阶速度—应力弹性波波 场分离等价的波动方程，并且采用高阶交错网格有限 这里仅考虑二维情况，二阶各向同性介质弹性波 [10, 11] 差分法实现波场分离数值模拟，其计算结果精度较高。 动方程的形式为 : 但由于上述这些算法需要进行算子分裂，这不但增加 [6～9] (1) 了计算过程的复杂性，同时也增加了许多计算量 。 作者简介：陈可洋，男，硕士，2009年毕业于大庆石油学院地球探测与信息技术专业，主要从事地震资料数字处理方法研究。 收稿日期：2009-06-17；修改日期：2010-02-24