1. School of Earth and Space Sciences, University of Science and Technology of China, Hefei, Anhui 230026, China; 2. Department of Earth and Space Sciences, Sou-thern University of Science and Technology, Shen-zhen, Guangdong 518055, China; 3. Shenzhen Key Laboratory of Deep offshore Oil and Gas Exploration Technology, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
Abstract：In numerical simulation of seismic wave,the limited space size of discrete grids usually leads to artificial reflections when calculating the outermost layer of the model.The perfectly matched layer (PML) method has achieved good results in absorbing artificial reflections,of which the complex frequency-shifted perfectly matched layer (CFS-PML) is particularly effective.Many studies that discussed the optimal settings of the CFS-PML parameters only consider homogeneous velocity models.However,there are usually complex velocity models with strong heterogeneous in field applications.Therefore,the CFS-PML parameters setting method of complex velocity models is explored based on the setting method of homogeneous velocity models.It is considered that the parameter setting is not only related to the source and grid properties,but also to the velocity of the absorption layer and the wavelength,that is,highly dependent on the velocity model.For a complex velocity model,the key to parameter setting is how to select a reference velocity,and three practical selection methods of reference velocity are given.By testing and comparing the absorption effects of three complex models,i.e.double-layer,multi-layer and Marmousi models, the best reference velocity selection method is defined as the median velocity of boundary medium,and then the setting method of CFS-PML optimization parameters for complex velocity models is proposed.
Berenger J P.A perfectly matched layer for the absorption of electromagnetic waves[J].Journal of Computational Physics,1994,114(2):185-200.
[2]
Chew W C,Liu Q H.Perfectly matched layers for elastodynamics:A new absorbing boundary condition[J].Journal of Computational Acoustics,1996,4(4):341-359.
[3]
Collino F,Tsogka C.Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics,2001,66(1):294-307.
[4]
Komatitsch D,Tromp J.A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation[J].Geophysical Journal International,2003,154(1):146-153.
[5]
Festa G,Vilotte J P.The Newmark scheme as velocity- stress time-staggering:an efficient PML implementation for spectral element simulations of elastodyna-mics[J].Geophysical Journal International,2005,161(3):789-812.
[6]
郭念民,吴国忱.基于PML边界的变网格高阶有限差分声波方程逆时偏移[J].石油地球物理勘探,2012,47(2):256-265.GUO Nianmin,WU Guochen.High-order finite difference method in reverse-time migration with variable grids based on PML boundary condition[J].Oil Geophysical Prospecting,2012,47(2):256-265.
[7]
高正辉,孙建国,孙章庆,等.基于完全匹配层构建新方法的2.5维声波近似方程数值模拟[J].石油地球物理勘探,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[J].Oil Geophysical Prospecting,2016,51(6):1128-1133.
[8]
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.
[9]
Roden J A,Gedney S D.Convolution PML (CPML):An efficient FDTD implementation of the CFS-PML for arbitrary media[J].Microwave and Optical Technology Letters,2000,27(5):334-339.
[10]
李佩笑,林伟军,张秀梅,等.常规分裂和非分裂完全匹配层吸收边界比较研究[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.
[11]
廉西猛,单联瑜,隋志强.地震正演数值模拟完全匹配层吸收边界条件研究综述[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 nume-rical simulation[J].Progress in Geophysics,2015,30(4):1725-1733.
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.
[14]
Drossaert F H,Giannopoulos A.A nonsplit complex frequency-shifted PML based on recursive integration for FDTD modeling of elastic waves[J].Geophysics,2007,72(2):T9-T17.
[15]
Li J,Dai J.Z-transform implementations of the CFS-PML[J].IEEE Antennas and Wireless Propagation Letters,2006,5(1):410-413.
[16]
Zhang W,Shen Y.Unsplit complex frequency-shifted PML implementation using auxiliary differential equations for seismic wave modeling[J].Geophysics,2010,75(4):T141-T154.
Correia D,Jin J M.On the development of a higher-order PML[J].IEEE Transactions on Antennas and Propagation, 2005,53(12):4157-4163.
[19]
冯海珂.复频移及双极点复频移完全匹配层在地震波高效数值模拟中的应用研究[D].安徽合肥:中国科学技术大学,2017.FENG Haike.Applications of the CFS-PML and Double-pole CFS-PML for Efficient Seismic Wave Numerical Simulation[D].University of Science and Technology of China,Hefei,Anhui,2017.
[20]
Gedney S D.The perfectly matched layer absorbing medium[J].Advances in Computational Electrodynamic,1998,263-344.
[21]
Marcinkovich C,Olsen K.On the implementation of perfectly matched layers in a three-dimensional fourth-order velocity-stress finite difference scheme[J].Journal of Geophysical Research,2003,108(B5):2276-2292.
[22]
Appelo D,Kreiss G.A new absorbing layer for elastic waves[J].Journal of Computatioal Physics,2006,215(2):642-660.
[23]
陈可洋.声波完全匹配层吸收边界条件的改进算法[J].石油物探,2009,48(1):76-79.CHEN Keyang.Improved algorithm for absorbing boundary condition of acoustic perfectly matched layer[J].Geophysical Prospecting for Petroleum,2009, 48(1):76-79.
[24]
陈可洋.完全匹配层吸收边界条件研究[J].石油物探, 2010,49(5):472-477.CHEN Keyang.Study on perfectly matched layer absorbing boundary condition[J].Geophysical Prospecting for Petroleum,2010,49(5):472-477.
[25]
罗玉钦,刘财.多轴复频移近似完全匹配层弹性波模拟[J].石油地球物理勘探,2019,54(5):1024-1033.LUO Yuqin,LIU Cai.Multi-axial complex-frequency shifting nearly perfectly matched layer for seismic forward modeling in elastic media[J].Oil Geophysical Prospecting,2019,54(5):1024-1033.
[26]
罗玉钦,刘财.近似完全匹配层边界条件吸收效果分析及衰减函数的改进[J].石油地球物理勘探,2018,53(5):903-913.LUO Yuqin,LIU Cai.Absorption effects in nearly perfectly matched layers and damping factor improvement[J].Oil Geophysical Prospecting,2018,53(5):903-913.
[27]
Bérenger J P.Numerical reflection from FDTD-PMLs:a comparison of the split PML with the unsplit and CFS PMLs[J].IEEE Transactions on Antennas and Pro-pagation,2002,50(3):258-265.
[28]
Bérenger J P.Application of the CFS PML to the absorption of evanescent waves in waveguides[J].IEEE Microwave and Wireless Components Letters,2002,12(6):218-220.
[29]
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):L20306.
