计算物理 ›› 2013, Vol. 30 ›› Issue (6): 843-854.

• 论文 • 上一篇    下一篇

双重弹性波波场分离数值模拟及相关理论证明

陈可洋1, 吴清岭1, 范兴才1, 陈树民1, 李来林1, 刘振宽2, 王建民1, 关昕1   

  1. 1. 中国石油大庆油田有限责任公司勘探开发研究院, 黑龙 江大庆 163712;
    2. 中国石油大庆油田有限责任公司勘探事业部, 黑龙江 大庆 163453
  • 收稿日期:2012-12-31 修回日期:2013-05-13 出版日期:2013-11-25 发布日期:2013-11-25
  • 作者简介:陈可洋(1983-),男,汉族,硕士,工程师,主要从事高精度弹性波传播模拟与逆时成像方法和地震资料数字处理方法研究,E-mail:keyangchen@163.com
  • 基金资助:
    国家重点基础研究发展计划(973)(2009CB219307)资助项目

Dual Elastic Wave Wavefield Separating Simulation Method and Related Theory Derivation

CHEN Keyang1, WU Qingling1, FAN Xingcai1, CHEN Shumin1, LI Lailin1, LIU Zhenkuan2, WANG Jianmin1, GUAN Xin1   

  1. 1. Exploration and Development Research Institute of Daqing Oilfield Company Limited of CNPC, Daqing 163712, China;
    2. Exploration Company of Daqing Oil Field Co. Ltd of CNPC, Daqing 163453, China
  • Received:2012-12-31 Revised:2013-05-13 Online:2013-11-25 Published:2013-11-25

摘要: 提出一种等效的双重弹性波波场分离数值模拟方法,用于模拟纯纵波和纯横波分离模式的质点振动速度、位移以及散度场和旋度场,并将该方法应用于全弹性波波动方程数值模拟中.同时,详细推导双重弹性波波场分离波动方程的高阶交错网格有限差分数值计算公式及其稳定性条件、数值频散关系和完全匹配层(PML)吸收边界条件.理论分析和数值计算均表明,该方法可以实现高精度双重弹性波波场分离数值模拟,且纯纵波和纯横波得到完全分离,边界吸收效果较好.与前人工作相比,存储量和计算时间均得到有效改善,数值计算结果进一步验证了该方法的优越性.

关键词: 双重波场分离, 弹性波, 稳定性条件, PML吸收边界, 高阶交错网格, 时域有限差分

Abstract: We present an equivalent dual elastic wave separation equation, which simulates particle-velocity, pressure, divergence and curl fields in pure P- and S- modes. The method is used in full elastic wave numerical simulations. We give complete derivations of explicit high-order staggered-grid finite difference discrete equations, together with stability condition, dispersion relation and perfectly matched layer (PML) absorbing boundary condition. Theoretical analysis and numerical simulations show that pare P-waves and S-waves in final numerical results are completely separated in the method. Effect of absorbing boundary is perfect. Storage and computing time requirements are greatly reduced compared with previous works.

Key words: dual wavefield separation, elastic wave, stability condition, PML absorbing boundary condition, high-order staggered-grid, time-domain finite-difference

中图分类号: