Abstract:The perfectly matched layer (PML) is a very efficient method and has become widely used in the seismic forward modeling.The development from split methods to various non-split methods greatly improves the computational efficiency of PML and reduces the required memory.The nearly perfectly matched layer (NPML) belongs to non-split PML,and it transforms directly the wavefield.Therefore,it does not change the form of the governing equations and is easy to be implemented.However,the NPML suffers from large spurious reflections at grazing incidence and evanescent wave.So,the convolution perfect matching layer (CPML) is still usually used.We implement the complex frequency shifted-NPML (CFS-NPML) to improve the capacity of absorbing grazing-incidence wave.We introduce different damping profiles that are proportional to each other in the orthogonal direction to make CFS-NPML stable called multi-axial CFS-NPML (MCFS-NPML).Meanwhile we design a new attenuation function which can improve the NPML absorption further.Eventually,the MCFS-NPML is able to absorb any incident wave more stably and efficiently.
Berenger J P.A perfectly matched layer for the ab-sorption of electromagnetics waves[J].Journal of Computational Physics,1994,114(2):185-200.
[2]
Berenger J P.On the reflection from Cummer's nearly perfectly matched layer[J].IEEE Microwave and Wireless Components Letters,2004,14(7):334-336.
[3]
Chew W C,Weedon W H.A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates[J].Microwave and Optical Technology Letters,1994,7(13):599-604.
[4]
Collino F,Monk P B.Optimizing the perfectly matched layer[J].Computer Methods in Applied Mechanics & Engineering,1998,164(1):157-171.
[5]
Kuzuoglu M,Mittra R.Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers[J].IEEE Microwave and Guided Wave Letters,1996,6(12):447-449.
[6]
Roden J A,Gedney S D.Convolution PML (C-PML):An efficient FDTD implementation of the CFS-PML for arbitrary media[J].Microwave and Optical Technology Letters,2000,27(5):334-339.
[7]
Komatitsch D,Martin R.An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation[J].Geophysics,2007,72(5):SM155-SM167.
[8]
Cummer S A.A simple nearly perfectly matched layer for general electromagnetic media[J].IEEE Microwave and Wireless Components Letters,2003,13(3):128-130.
[9]
Hu W Y,Abubakar A,Habashy T M.Application of the nearly perfectly matched layer in acoustic wave modeling[J].Geophysics,2007,72(5):SM169-SM175.
[10]
Chen J Y.Application of the nearly perfectly matched layer for seismic wave propagation in 2D homoge-neous isotropic media[J].Geophysical Prospecting,2011,59(4):662-672.
[11]
Chen J Y.Nearly perfectly matched layer method for seismic wave propagation in poroelastic media[J].Canadian Journal of Exploration Geophysics,2012,37(1):22-27.
[12]
Festa G,Delavaud E,Vilotte J P.Interaction between surface waves and absorbing boundaries for wave propagation in geological basins:2D numerical simulations[J].Geophysical Research Letters,2005,32(20):55-102.
[13]
Mezafajardo K C,Papageorgiou A S.A nonconvolutional,split-field,perfectly matched layer for wave propagation in isotropic and anisotropic elastic media:stability analysis[J].Bulletin of the Seismological Society of America,2008,98(4):1811-1836.
[14]
Collino F,Tsogka C.Application of the PML abso-rbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics,2001,66(1):294-307.
[15]
Groby J,Tsogka C.A time domain method for mode-ling viscoacoustic wave propagation[J].Journal of Computational Acoustics,2006,14(2):201-236.
[16]
王守东.声波方程完全匹配层吸收边界[J].石油地球物理勘探,2003,38(1):31-34.WANG Shoudong.Absorbing boundary condition for acoustic wave equation by perfectly matched layer[J].Oil Geophysical Prospecting,2003,38(1):31-34.
[17]
秦臻,任培罡,姚姚,等.弹性波正演模拟中PML吸收边界条件的改进[J].中国地质大学学报,2009,34(4):658-664.QIN Zhen,REN Peigang,YAO Yao,et al.Improvement of PML absorbing boundary conditions in elastic wave forward modeling[J].Journal of China University of Geosciences,2009,34(4):658-664.
[18]
张显文,韩立国,黄玲,等.基于递归积分的复频移PML弹性波方程交错网格高阶差分[J].地球物理报,2009,52(7):1800-1807.ZHANG Xianwen,HAN Liguo,HUANG Ling,et al.A staggered-grid high-order difference method of complex frequency-shifted PML based on recursive intergration for elastic wave equation[J].Chinese Journal of Geophysics,2009,52(7):1800-1807.
刘洋.波动方程时空域有限差分数值解及吸收边界条件研究进展[J].石油地球物理勘探,2014,49(1):35-46.LIU Yang.Research progress on time-space domain finite difference numerical solution and absorption boundary conditions of wave equation[J].Oil Geophysical Prospecting,2014,49(1):35-46.
[21]
李佩笑,林伟军,张秀梅,等.常规分裂和非分裂完全匹配层吸收边界比较研究[J].声学学报,2015,40(1):44-53.LI Peixiao,LIN Weijun,ZHANG Xiumei,et al.Comparisons for regular splitting and non-splitting perfectly matched layer absorbing boundary conditions[J].Acta Acustica,2015,40(1):44-53.
[22]
廉西猛,单联瑜,隋志强.地震正演模拟完全匹配层吸收边界条件研究综述[J].地球物理学进展,2015,30(4):1725-1733.LIAN Ximeng,SHAN Lianyu,SUI Zhiqiang.An overview of research on perfectly matched layers absorbing boundary condition of seismic forward numeri-cal simulation[J].Progress in Geophysics,2015,30(4):1725-1733.
[23]
黄建平,黄金强,李振春,等.TTI介质一阶qP波方程交错网格伪谱法正演模拟[J].石油地球物理勘探, 2016,51(3):487-496.HUANG Jianping,HUANG Jinqiang,LI Zhenchun,et al.Staggered grid pseudo-spectral numerical simulation of first-order quasi P-wave equation for TTI media[J].Oil Geophysical Prospecting,2016,51(3):487-496.
[24]
张衡,刘洪,李博,等.TTI介质声波方程分裂式PML吸收边界条件研究[J].石油物探,2017,56(3):349-361.ZHANG Heng,LIU Hong,LI Bo,et al.The research on split PML absorbing boundary conditions of acoustic equation for TTI media[J].Geophysical Prospecting for Petroleum,2017,56(3):349-361.
[25]
刘东洋,彭苏萍,师素珍,等.基于Lebedev网格的TTI介质二维三分量正演模拟[J].石油地球物理勘探,2018,53(2):288-296.LIU Dongyang,PENG Suping,SHI Suzhen,et al.2D three-component seismic forward modeling in TTI media based on the Lebedev grid[J].Oil Geophysical Prospecting,2018,53(2):288-296.
[26]
成景旺,毛宁波,吕晓春,等.非分裂完全匹配层边界存储时间域全波形反演[J].石油地球物理勘探,2018,53(4):754-764.CHENG Jingwang,MAO Ningbo,LYU Xiaochun,et al.Time-domain full waveform inversion with CFS-NPML boundary storage[J].Oil Geophysical Prospecting,2018,53(4):754-764.
[27]
胡建林,宋维琪,张建坤,等.交错网格有限差分正演模拟的联合吸收边界[J].石油地球物理勘探,2018,53(5):914-920.HU Jianlin,SONG Weiqi,ZHANG Jiankun,et al.Joint absorbing boundary in the staggered-grid finite difference forward modeling simulation[J].Oil Geophysical Prospecting,2018,53(5):914-920.
[28]
汪勇,徐佑德,高刚,等.二维黏滞声波方程的优化组合型紧致有限差分数值模拟[J].石油地球物理勘探,2018,53(6):1152-1164.WANG Yong,XU Youde,GAO Gang,et al.Numerical simulation of 2D visco-acoustic wave equation with an optimized combined compact difference scheme[J].Oil Geophysical Prospecting,2018,53(6):1152-1164.