Abstract:There are two kinds of problems in perfectly matched layers (PML) introduced in the seismic wave field simulation:one is how to reduce the extra storage capacity in the practical boundary application and enhance the calculation speed; the other is how to absorb incident wave more quickly and make the absorption better.The first problem is well solved with a few advanced non-split perfectly matched layer approaches.For this issue,we adopt a newly-developed non-split perfectly matched layer approach,which transforms directly the wavefield,and will not change wavefield form and use convolution.For the second issue,there is not much progress in boundary absorption enhancement due to the discrepancy complexity caused by discretization.Currently,we decrease the discretization discrepancy by adding the number of layers.In this paper,we analyze the property of damping term and design a new attenuation function which can improve the absorption effect.The function could weaken the reflection amplitude value by 20%~60% when the number of layers is 5~10 compared with conventional attenuation function.
Berenger J P.A perfectly matched layer for the absorption of electromagnetic waves.Journal of Computational Physics,1994,114(2):185-200.
[2]
Chew W C,Weedon W H.A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates.Microwave and Optical Technology Letters,1994,7(13):599-604.
[3]
Collino F,Monk P B.Optimizing the perfectly matched layer.Computer Methods in Applied Mechanics & Engineering,1998,164(12):157-171.
[4]
Chew W C,Liu Q H.Perfectly matched layers for ela-stodynamics:A new absorbing boundary condition.Journal of Computational Acoustics,1996,4(4):341-359.
[5]
Hastings F D,Schneider J B,Broschat S L.Application of the perfectly matched layer(PML)absorbing boundary condition to elastic wave propagation.Journal of the Acoustical Society of America,1996,100(5):3061-3069.
[6]
Komatitsch D,Tromp J.A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation.Geophysics,2003,72(5):SM155-SM167.
[7]
Drossaert F,Giannopoulos A.A nonsplit complex frequency-shifted PML based on recursive integration for FDTD modeling of elastic waves.Geophysics,2007,72(2):9-17.
[8]
Cummer S A.A simple nearly perfectly matched layer for general electromagnetic media.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.Geophysics,2007,72(5):SM169-SM175.
[10]
Chen J Y.Application of the nearly perfectly matched layer for seismic wave propagation in 2D homogeneous isotropic media.Geophysical Prospecting,2011,59(4):662-672.
[11]
Chen J Y.Nearly perfectly matched layer method for seismic wave propagation in poroelastic media.Canadian Journal of Exploration Geophysics,2012,37(1):22-27.
[12]
Kuzuoglu M,Mittra R.Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers.IEEE Microwave and Guided Wave Letters,1996,6(12):447-449.
[13]
Roden J A,Gedney S D.Convolution PML (C-PML):An efficient FDTD implementation of the CFS-PML for arbitrary media.Microwave and Optical Technology Letters,2000,27(5):334-339.
[14]
Komatitsch D,Martin R.An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation.Geophysics,2007,72(5):SM155-SM167.
高正辉,孙建国,孙章庆等.基于完全匹配层构建新方法的2.5维声波近似方程数值模拟.石油地球物理勘探,2016,51(6):1128-1133.Gao Zhenghui,Sun Jianguo,Sun Zhangqing et al.2.5D acoustic approximate equation numerical simulation with a new construction method of the perfectly matched layer.OGP,2016,51(6):1128-1133
[17]
刘洋.波动方程时空域有限差分数值解及吸收边界条件研究进展.石油地球物理勘探,2014,49(1):35-46.Liu Yang.Research progress on time-space domain finite difference numerical solution and absorption boundary conditions of wave equation.OGP,2014,49(1):35-46.
[18]
李佩笑,林伟军,张秀梅等.常规分裂和非分裂完全匹配层吸收边界比较研究.声学学报,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.Acta Acustica,2015,40(1):44-53
[19]
张衡,刘洪,李博等.TTI介质声波方程分裂式PML吸收边界条件研究.石油物探,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.GPP,2017,56(3):349-361
[20]
张显文,韩立国,黄玲等.基于递归积分的复频移PML弹性波方程交错网格高阶差分.地球物理学报,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.Chinese Journal of Geophysics,2009,52(7):1800-1807.
[21]
黄建平,黄金强,李振春等.TTI介质一阶qP波方程交错网格伪谱法正演模拟.石油地球物理勘探,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.OGP,2016,51(3):487-496.
[22]
廉西猛,单联瑜,隋志强.地震正演模拟完全匹配层吸收边界条件研究综述.地球物理学进展,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 numerical si-mulation.Progress in Geophysics,2015,30(4):1725-1733.
[23]
王守东.声波方程完全匹配层吸收边界.石油地球物理勘探,2003,38(1):31-34.Wang Shoudong.Absorbing boundary condition for acoustic wave equation by perfectly matched layer.OGP,2003,38(1):31-34.
[24]
秦臻,任培罡.弹性波正演模拟中PML吸收边界条件的改进.中国地质大学学报,2009,34(4):658-664.Qin Zhen,Ren Peigang.Improvement of PML absor-bing boundary conditions in elastic wave forward mo-deling.Journal of China University of Geosciences,2009,34(4):658-664.
Collino F,Tsogka C.Application of the PML absor-bing layer model to the linear elastodynamic problem in anisotropic heterogeneous media.Geophysics,2001,66(1):294-307.
[27]
Groby J,Tsogka C.A time domain method for mode-ling viscoacoustic wave propagation.Journal of Computational Acoustics,2006,14(2):201-236.
[28]
陈可洋.完全匹配层吸收边界条件研究.石油物探,2010,49(5):472-477.Chen Keyang.Study on perfectly matched layer absorbing boundary condition.GPP,2010,49(5):472-477.