Analyzing PML absorbing impact factors with wave equation simulation
LIU Rui-He1, YIN Xing-Yao1, PU Yi-Tao2
1. School of Geoscience, China University of Petroleum, Qingdao 266580, China;
2. Bureau of Geophysical Prospecting INC., China National Petroleum Corporation, Korla 841001, China
In seismic numerical simulation, the perfectly matched layer (PML) boundary can attenuate the seismic wave effectively, and the snapshots and record have no boundary reflection interferes.The authors further analyzed the attenuation mechanism of split perfectly matched layer (SPML).By wave equation simulation, the authors analyzed the problem as to how the parameters affect absorption and attenuation characteristics of SPML boundary, such as source dominant frequency, spatial grid spacing and medium speed.As a result, the authors obtained the results of PML attenuation based on different conditions.Through analyzing and comparing the results, the conclusions have been reached that source dominant frequency has no direct influence on PML attenuation, spatial grid spacing is in inverse proportion to the PML attenuation, and high-speed medium slows down PML attenuation effect.By Marmousi2 test, the authors confirmed the validity of the conclusions.The selection method of SPML parameter in complex model or different scales of numerical simulation was given, which lays a foundation for quantitative selection of SPML boundary parameters in numerical simulation.
[1] Cerjan C, Kosloff D, Kosloff R, et al. A nonreflecting boundary condition for discrete acoustic and elastic wave equations[J].Geophysics,1985,50(4):705-708.
[2] Berenger J. A perfectly matched layer for the absorption of electromagnetic waves[J].Journal of Computational Physics,1994,114(2):185-200.
[3] Clayton R, Engquist B. Absorbing boundary conditions for wave-equation migration[J].Geophysics,1980, 45(5):895-904.
[4] 郭念民, 吴国忱.基于PML边界的变网格高阶有限差分声波方程逆时偏移[J].石油地球物理勘探,2012,47(2):256-265.
[5] 孙林洁, 印兴耀.基于PML边界条件的高倍可变网格有限差分数值模拟方法[J].地球物理学报,2011,54(6):1614-1623.
[6] 张岩, 吴国忱.TTI介质qP波逆时偏移中伪横波噪声压制方法[J].地球物理学报,2013,56(6):2065-2076.
[7] 印兴耀,浦义涛,梁锴,等.VTI介质准P波旋转交错有限差分数值模拟[J].CT理论与应用研究,2014,23(5):771-783.
[8] Hastings F D,Schneider J B,Broschat S L.Application of the perfectly matched layer (PML) absorbing boundary condition to elastic wave propagation[J]. The Journal of the Acoustical Society of America,1996,100(5):3061-3069.
[9] 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.
[10] Liu Q, Tao J. The perfectly matched layer for acoustic waves in absorptive media[J]. The Journal of the Acoustical Society of America,1997,102(4):2072-2082.
[11] Zeng Y, He J, Liu Q. The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media[J].Geophysics,2001,66(4):1258-1266.
[12] 张鲁新,符力耘,裴正林. 不分裂卷积完全匹配层与旋转交错网格有限差分在孔隙弹性介质模拟中的应用[J]. 地球物理学报,2010,53(10):2470-2483.
[13] Pu Y T, Yin X Y, Cao D.The PML boundary of small-scale seismic modeling[J]. SEG Technical Program Expanded Abstracts,2014:3493-3497.
[14] 王守东. 声波方程完全匹配层吸收边界[J]. 石油地球物理勘探, 2003,38(1):31-34.
[15] 张剑锋,徐义. 地震波数值模拟的非规则网格PML吸收边界[J]. 地球物理学报, 2008, 51(5):1520-1526.