稳定的二维TTI介质一阶qP波方程RSGFD数值模拟  

RSGFD numerical simulation of stable 2D TTI first-order qP-wave equations

作　　者：王小丹 印兴耀[1] 杨富森 

机构地区：[1]中国石油大学(华东)地球科学与技术学院,山东青岛266580

出　　处：《物探与化探》2016年第1期135-142,共8页Geophysical and Geochemical Exploration

基　　金：国家自然科学基金(41204085);国家重点基础研究发展计划(“973计划”)项目(2013CB228604)共同资助

摘　　要：相比于VTI介质q P波数值模拟方法,考虑倾角因素的TTI介质q P波数值模拟方法能够更加准确地描述各向异性介质中波场的传播特征。常规的声学近似方法往往会造成TTI介质中倾角急剧变化的区域出现数值不稳定,笔者首先引入一个各向异性控制参数,推导了稳定形式的TTI介质二阶q P波方程,保证q P波场的稳定传播;其次,通过引入波场的伪速度分量,推导了其等价的一阶应力—速度形式及相应的PML边界条件,并应用优化的旋转交错网格有限差分(RSGFD)方法实现了精确的数值模拟。数值结果表明:TTI介质一阶q P波方程能够稳定、有效地模拟q P波的运动学特征,利用优化的RSGFD方法可以得到精确的合成地震记录,同时可以相对地提高计算效率。Compared with VTI qP-wave modeling, TFI qP-wave numerical simulation method which considers dip angle factor could be more accurate in describing the characteristics of wave-field propagation in real anisotropic media.Because of the rapid change of dip an- gle in TYI media,using the conventional acoustic approximated method may cause numerical instability.In this paper, the authors first derived the stable TI＇I second-order qP-wave equations by introducing an anisotropic control parameter to ensure the stability of qP-wave propagation.Then ,the equivalent first-order stress-velocity forms were deduced through introducing the pseudo-velocity components of the wave-fields, and the corresponding PML boundary condition was given.Finally, the authors implemented accurate numerical simula- tion using the optimal rotated stagger-grid finite-difference （RSGFD） method. Numerical results demonstrate that TTI first-order qP- wave equations can describe the kinematics features of qP-wave effectively and stably, and the application of optimal RSGFD method can acquire accurate synthetic seismic recordings and enhance computational efficiency relatively.

关 键 词：TTI介质 一阶应力-速度 数值模拟 旋转交错网格有限差分 控制参数 稳定性 PML边界条件 

分 类 号：P631.4[天文地球—地质矿产勘探]