等价各向异性介质一阶拟声波方程正演模拟及其效率分析.pdfVIP

第 37卷 第 4期 地 震 学 报 Vo1．37，NO．4 2015年 7月 (648—660) ACTA SEISMOLOGICA SINICA Ju1．，2O15 杨富森，李振春 ，王小丹 ，张凯，刘成斋．2015．等价各向异性介质一阶拟声波方程正演模拟及其效率分析．地震学报 ， 37(4)：648—660．doi：10．11939／jass．2015．04．011． YangF S，LiZC，WangXD ，ZhangK，LiuC Z．2015．Equivalentanisotropicfirst～orderpseudo—acousticwaveequa— tionsmodelinganditsefficiencyanalysis．ActaSeismologicaSinica，37(4)：648—660．doi：10．11939／)ass．2015．04．OlI． 等价各 向异性介质一阶拟声波 方程正演模拟及其效率分析 杨富森”· 李振春D 王小丹” 张 凯 刘成斋 1)中国山东青岛266580中国石油大学(华东)地球科学与技术学院 2)中国山东东营257022中国石化胜利油 田物探研究院 摘要 为克服各 向异性弹性波动方程正演模拟的局限，本文研究了各 向异性介质拟声波方程 的交错网格有限差分数值解法．首先，从VTI介质胡克定律和 qP-qSV波频散关系两种思路 出发，通过声假设近似 ，给出了两种不同形式的VTI介质一阶拟声波方程 ，并通过引入波场 的伪速度分量，推导了一种新的VTI介质一阶应力一速度方程，并通过旋转坐标系将其推广 到TTI介质中；其次，构造了一阶拟声波方程的交错网格高阶有限差分格式，并推导了相应 的PML边界条件 ；最后 ，对本文方法中固有 的qSV人为干扰波的产生机制和压制方法进行 了简单讨论．数值结果表明：3种一阶拟声波方程在运动学和动力学上是等价的，相对于各 向异性弹性波正演模拟，其节省了内存，提高了计算效率；各向异性因素会影响反射波旅行 时和振幅等波场特征，在后续的处理、反演和解释 中不可忽略；VTI介质 HESS模型的逆时 偏移结果也验证 了本文方法 的合理性． 关键词 各向异性介质 一阶拟声波 PML边界条件 交错 网格 正演模拟 qSV人为 干扰波 doi：10．11939／jass．2015．04．011 中图分类号：P315．3 1 文献标志码 ：A Equivalentanisotropicfirst—-orderpseudo-a·cousticwave equationsmodelingand itsefficiencyanalysis YangFusen LiZhenchun W angXiaodan ZhangKai LiuChengzhai 1)SchoolofGeosciences，ChinaUniversityofPetroleum，ShandongQingdao266580，China 2)GeophysicalResearchInstituteofSINOPECShengliOilfield，ShandongDongying 257022，China Abstract：Inordertoovercomethelimitationsofanisotropicelasticwaveequa— tionsmodeling，thispaperstudiedanisotropicpseudo～acousticwavemodelingby usingstaggered—gridfinite—differencemethod．Startingfrom twodifferentideas， i．e．，Hooke’Slaw andqP—qSV dispersionrelation forVTImedia，twofirst— orderpseudo—acousticwaveequationsin two formsforVTImediaareg