TTI介质纯准P波一阶压力-速度方程及求解方法TTI Medium Pure Quasi-P Wave First-Order Pressure-Velocity Equation and Solution Method
李佳珂,张会星,白冰,张建敏
摘要(Abstract):
以具有倾斜对称轴的横向各向同性(Titled Transverse Isotropic,TTI)介质中纯准P波二阶方程为基础,通过引入辅助波场实现了方程的降阶,推导出了TTI介质纯准P波一阶压力-速度方程,并利用旋转交错网格高阶有限差分方法对方程进行数值求解,给出了完全匹配层(Perfectly Matched Layers,PML)吸收边界条件在该方程中应用的方法。理论分析及数值模拟结果表明:相对于二阶形式的准P波方程,该TTI介质纯准P波一阶压力-速度方程在计算精度相同的情况下,具有更高的计算效率;同时能更准确地描述对称轴参数变化剧烈的TTI介质中准P波的传播规律。
关键词(KeyWords): TTI介质;纯准P波;旋转交错网格;PML边界;Wedge模型
基金项目(Foundation): 中央高校基本科研业务费专项(201822011);; 国家自然科学基金(41674118);; 国家科技重大专项(2016ZX05027-002)联合资助
作者(Author): 李佳珂,张会星,白冰,张建敏
参考文献(References):
- [1]王璐琛,常旭,王一博.TTI介质的交错网格伪P波正演方法[J].地球物理学报,2016,59(3):1046-1057.
- [2]王娟,李振春,孙小东,陶丽.TTI介质逆时偏移成像[J].石油地球物理勘探,2012,47(4):573-577.
- [3]Chu C,Macy B K,Anno P D.An accurate and stable wave equation for pure acoustic TTI modeling[C]//81th Annual International Meeting,Seg Technical Program Expanded Abstracts,2010:179-184.
- [4]Zhang Y,Zhang H.A stable TTI reverse time migration and its implementation[J].Geophysics,2009,76(3):2794-2799.
- [5]Zhou H,Zhang G,Bloor R.An anistropic acoustic wave equation for modeling and migration in 2D TTI media[C]//76th Annual International Meeting,SEG,Espanded Abstracts,2006.
- [6]Fletcher R,Du X,Fowler P J.A new pseudo-acoustic wave equation for TI media[J].SEG Technical Program Expanded Abstracts,2008,27(1):2082-2086.
- [7]Fletcher R P,Du X,Fowler P J.Reverse time migration in tilted transversely isotropic(TTI)media[J].Geophysics,2009,4(6):179-187.
- [8]Li B,Li M,Liu H W.Stability of reverse time migration in TTI media[J].Geophys(in Chinese),2012,55(4):1366-1375.
- [9]Zhan G,Pestana R C,Stoffa P L.An efficient hybrid pseudospectral/finite-difference scheme for solving the TTI pure P-wave equation[J].Journal of Geophysics and Engineering,2013,10(2):25004-25011.
- [10]Zhang J M,He B S,Tang H G.Pure quasi-P wave equation and numerical solution in 3D TTI media[J].Applied Geophysics,2017,14(1):125-132.
- [11]Alkhalifah T.An acoustic wave equation for anisotropic media[J].Geophysics,1998,65(4):1239-1250.
- [12]Alkhalifah T.The offset-midpoint traveltime pyramid in transversely isotropic media[J].Geophysics,2000,65(4):1316-1325.
- [13]吴国忱,梁锴.VTI介质频率-空间域准P波正演模拟[J].石油地球物理勘探,2006,40(5):535-545.
- [14]Du X,Fletcher R,Fowler P J.A new pseud-cracoustic wave equation for VTI media[C]//Presented at the 70th EAGE Conference&Exhibition,2008.
- [15]Duvenck E,P Milik,P M Bakker.Acoustic VTI wave equations and their application for anistropic reverse-time migration[J].SEG Technical Program Expanded Abstracts,2008,27(1):2186-2190.
- [16]Hestholm S.Acoustic VTI modeling using high-order finite differences[J].Geophysics,2009,74(5):67-72.
- [17]韩令贺,何兵寿.VTI介质一阶准P波方程正演模拟及边界条件[J].石油地球物理勘探,2010,45(6):819-825.
- [18]Fowler P J,Du X,Fletcher R P.Coupled equations for reverse time migration in transversely isotropic media[J].Geophysics,2010,75(1):S11-S22.
- [19]Xu S,Zhou H.Accurate simulation of pure quasi-P-waves in complex anisotropic media[J].Geophysics,2014,79(6):341-348.
- [20]Xu S,Tang B,Mu J,Zhou H.Quasi-P wave propagation with an elliptic differential operator[C]//85th New Orleans Annual Meeting,SEG,Expanded Abstracts,2015.
- [21]Thomsen L.Weak elastic anisotropy[J].Geophysics,1986,51(10):1954-1996.
- [22]严红勇,刘洋.黏弹TTI介质中旋转交错网格高阶有限差分数值模拟[J].地球物理学报,2012,55(4):1354-1365.
- [23]张晶,何兵寿,张会星.含直立裂隙介质的弹性波动方程正演模拟[J].中国煤炭地质,2008,20(1):50-54.
- [24]常海明,何兵寿,杨佳佳,陈婷.逆时偏移中的归一化加权互相关成像条件及应用[J].中国煤炭地质,2014,26(1):51-55.
- [25]田雪丰,张奋轩,曹新领.弹性波逆时偏移及其成像条件[J].中国煤炭地质(原中国煤田地质),2002,14(2):58-61.
- [26]Liu Y,Li C C and Mou Y G.Finite-difference numerical modeling of any even-order accuracy[J].Chinese J Oil Geophysics Prospecting(in Chinese),1998,33(1):1-10.
- [27]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419.
- [28]王芝尧,韩国猛,付东立,牟智全,刘志英.二维随机介质模型正演指导复杂地层地震解释——以薄层砂泥岩交互地层为例[J].中国煤炭地质,2012,24(1):45-48.
- [29]Berenger J P.A perfectly matched layer for the absorption of electromagnetic waves[J].Journal of Computational Physics,1994,114(2):185-200.
- [30]刘学义,何兵寿.VTI介质中弹性波方程正演的一阶混合吸收边界[J].中国煤炭地质,2015,27(3):59-63.