基于散度和旋度的弹性波波场分离数值模拟方法-物探与化探.pdfVIP

　第３４卷第１期 物　探　与　化　探 Ｖｏｌ．３４，Ｎｏ．１　 　２０１０年２月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ Ｆｅｂ．，２０１０　 基于散度和旋度的弹性波波场分离数值模拟方法 １ １ ２ １ １ 陈可洋 ，吴清岭 ，杨微 ，李来林 ，范兴才 （１．中国石油大庆油田有限责任公司勘探开发研究院，黑龙江 大庆　１６３７１２；２．中国石油大庆油 田有限责任公司采油六厂，黑龙江 大庆　１６３７１２） 摘 要：弹性波波场分离数值模拟方法是研究混合波场中纯纵波和纯横波波场传播规律的一种重要手段，通常采用 算子分裂来实现，这不但增加了计算量，而且还浪费了大量计算时间。因此提出了基于散度和旋度的弹性波波场 分离数值模拟新方法，仅依赖于交错网格来计算出散度场（纯纵波）和旋度场（纯横波），从而提高计算效率。在人 工截断边界处采用完全匹配层吸收边界条件（ＰＭＬ）。数值模拟结果表明，该方法能够从混合波场中准确地分离出 纯纵横波波场，且精度高，边界吸收效果好，计算速度快，计算过程稳定，可用于指导实际多波多分量地震资料处 理，同时指出本文算法对相互耦合的波型（如面波）无法实现准确分离。 关键词：高阶交错网格；波场分离；弹性波动方程；散度和旋度 中图分类号：Ｐ６３１．４　　　文献标识码：Ａ　　　文章编号：１０００－８９１８（２０１０）０１－０１０３－０５ 　　目前，地震波数值模拟方法已成为准确认识弹 １　波场分离数值模拟原理 性波传播规律的一种重要手段，其中研究纯纵波和 纯横波波场的弹性波波场分离数值模拟方法逐步兴 １．１　一阶双曲型弹性波动方程及差分格式的构成 起，开始作为数值模拟理论和反演问题新的研究方 二维一阶双曲型各向同性介质弹性波动方程的 ［１］ ［６－７］ 向。马德堂等 首次提出满足纵波为无旋场、横波 形式如下 为无散场的等价方程思路，采用虚谱法求解等价的 ｖ ｖ τ    ｘｘ ２ ２ ｘ ２ ｚ ＝（ －２ ） ＋ ， ［２－３］  ｔ α β ｘ α ｚ 二阶弹性波动方程 ，实现了波场分离数值模拟，     ｖ ｖ 获得较好的数值模拟结果，但是伪谱法边界处理困 τ    ｚｚ ２ ２ ｚ ２ ｘ ＝（ －２ ） ＋ ， α β α ｔ