徐啸,李刚,郭英杰,王震,李鹏程,郭锦桦,吕晓磊

• 1:国网能源哈密煤电有限公司
• 2:中煤科工西安研究院(集团)有限公司
• 3:西安煤科透明地质科技有限公司

摘要(Abstract)：

地球内部可视为由含流体孔隙介质构成，即双相介质模型。为了进一步认识双相介质中波场传播规律和特征，基于双相各向异性介质，以BIOT理论为基础推导了双相各向异性介质二维三分量一阶速度—应力弹性波方程及其差分格式。传统交错网格有限差分技术在双相各向异性介质波场模拟中经常出现数值频散现象，为了解决这一难题，利用传统交错网格有限差分技术，对均匀和非均匀双相介质进行了波场模拟，改进了传统交错网格点处应力的定义方式，并将改进前后的波场快照进行了对比。对比结果表明：改进后的网格定义方式能够较好的模拟弹性波在均匀和非均匀双相介质中的传播；改进后的网格定义方式在计算中无需应力的插值，减少了计算误差；改进后的正演方法在非均匀双相介质波场模拟中无波场异常现象，有效降低了数值频散,提高了模拟精度，是一种有效且实用的交错网格有限差分方法。

关键词(KeyWords)： 双相各向异性介质;交错网格;有限差分;网格定义;数值频散

Abstract:

Keywords:

基金项目(Foundation): 中国煤炭科工集团重点项目“透明地质综合保障系统集成与应用示范”（2024-TD-ZD001-02）;中国煤炭科工集团有限公司科技创新创业资金专项重点项目“基于动态更新地质模型的智能采掘交互技术研究”（2024-TDZD001-01）

作者(Author): 徐啸,李刚,郭英杰,王震,李鹏程,郭锦桦,吕晓磊

参考文献(References)：

• [1]李振春，卫二祥，黄建平，等.基于一阶速度-应力方程的起伏地表最小二乘逆时偏移[J].中国石油大学学报（自然科学版），2025,49(2):108-117.
• [2]孙辉，岳玉波，李猛.基于Born散射理论的二维黏声介质高斯波束正演[J].地球物理学报，2021, 64(2):637-644.
• [3]BIOT M A. Theory of propagation of elastic waves in a fluidsaturated porous solid, I:Low-frequency range[J]. The Journal of the Acoustic Society of America, 1956, 28:168-178.
• [4]BIOT M A. Theory of propagation of elastic waves in a fluidsaturated porous solid, II:High-frequency range[J]. The Journal of the Acoustic Society of America, 1956, 28:179-191.
• [5]BIOT M A. Mechanics of deformation and acoustic propagation in porous media[J]. Journal of Applied Physics, 1962, 33:1482-1489.
• [6]ZHU X, MCMECHAN G A. Numerical simulation of seismic responses of poroelastic reservoirs using Biot theory[J]. Geophysics,1991, 56(3):328-339.
• [7]杨顶辉.双相各向异性介质中弹性波方程的有限元解法及波场模拟[J].地球物理学报，2002, 45(4):575-583.
• [8]王秀明，张海澜，王东.利用高阶交错网格有限差分法模拟地震波在非均匀孔隙介质中的传播[J].地球物理学报，2003, 46(6):842-849.
• [9]裴正林.双相各向异性介质弹性波传播交错网格高阶有限差分法模拟[J].石油地球物理勘探，2006,41(2):137-143.
• [10]李卫志.双相介质地震波场数值模拟研究[D].四川成都：成都理工大学，2007.
• [11]张文忠. Biot介质的交错网格差分法波场模拟研究[D].北京：中国地质大学（北京），2007.
• [12]苏云.各向异性及双相介质弹性波正演模拟与波场特征分析[D].四川成都：成都理工大学，2010.
• [13]邓万九.双相煤岩介质中的地震波传播特性研究[D].河南焦作：河南理工大学，2021.
• [14]王家麟，孙成，李昊.高精度双相各向异性弹性介质地震波正演与波场特征分析[C]//中国地球物理学会油气地球物理专业委员会.第五届油气地球物理学术年会论文集.中国矿业大学（北京）地球科学与测绘工程学院；2023:678-683.
• [15]BERENGER J P. A perfectly matched layer for the absorption of electromagnetic waves[J]. Journal of computational physics, 1994, 114(2):185-200.
• [16]COLLINO F, TSOGKA C. Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J]. Geophysics, 2001, 66(1):294-307.
• [17]孙卫涛，杨慧珠.双相各向异性介质弹性波场有限差分正演模拟[J].固体力学学报，2004, 25(1):21-28.
• [18]孙瑞艳. TTI介质旋转交错网格有限差分及其组合边界条件[D].山东青岛：中国石油大学，2010.