Gaussian-type weighted hybrid absorbing boundary for elastic wave simulation and its acceleration on GPU
WANG Shaowen1, SONG Peng1,2,3, TAN Jun1,2,3, XIE Chuang1, MAO Shibo1, WANG Qianqian1
1. College of Marine Geo-sciences, Ocean University of China, Qingdao, Shandong 266100, China; 2. Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong 266100, China; 3. Key Lab of Submarine Geosciences and Prospecting Techniques, MOE, Ocean University of China, Qingdao, Shandong 266100, China
Abstract:Usually, hybrid absorbing boundary conditions (HABC) in simulating elastic wave based on elastic wave equation use linear or exponential weighted coefficient. It is hard to get excellent absorbing effect in both inner and outer boundaries, and the boundary difference scheme is not suitable for efficient acceleration of GPU. This study proposes a HABC method based on a Gaussian-type weighted coefficient to achieve better absorption, and a finite-difference scheme of HABC based on the first-order Higdon absorbing boundary conditions to realize GPU acceleration. Numerical simulation has demonstrated that the HABC based on a Gaussian-type hybrid weighted coefficient can achieve better absorbing effect in both inner and outer boundaries than the HABC based on a linear and exponential hybrid weighted coefficient. In addition, the application of the finite-difference scheme can improve GPU acceleration efficiency, making more suitable for large-scale numerical simulation on elastic wave field based on a complex model.
王绍文, 宋鹏, 谭军, 解闯, 毛士博, 王倩倩. 弹性波数值模拟中的高斯型混合吸收边界条件及其GPU并行[J]. 石油地球物理勘探, 2021, 56(3): 485-495.
WANG Shaowen, SONG Peng, TAN Jun, XIE Chuang, MAO Shibo, WANG Qianqian. Gaussian-type weighted hybrid absorbing boundary for elastic wave simulation and its acceleration on GPU. Oil Geophysical Prospecting, 2021, 56(3): 485-495.
Israeli M,Orszag,S A. Approximation of radiation boundary conditions[J].Journal of Computational Physics,1981,41(1):115-135.
[2]
Cerjan C,Kosloff R,Reshef M. A nonreflecting boundary condition for discrete acoustic and elastic wave equations[J].Geophysics,1985,50(4):705-708.
[3]
Kosloff R,Kosloff D. Absorbing boundaries for wave propagation problems[J].Journal of Computational Physics,1986,63(2):363-376.
[4]
Berenger J P. A perfectly matched layer for the absorption of electromagnetic waves[J].Journal of Computational Physics,1994,114(2):185-200.
[5]
王守东. 声波方程完全匹配层吸收边界[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.
[6]
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.
[7]
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.
[8]
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.
[9]
Martin R,Komatitsch D. An unsplit convolutional perfectly matched layer technique improved at grazing incidence for the viscoelastic wave equation[J].Geophysical Journal International,2009,179(1):334-344.
罗玉钦,刘财. 多轴复频移近似完全匹配层弹性波模拟[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.
[12]
杨茜娜,张伟. 复频移完全匹配层在复杂速度模型中的应用[J].石油地球物理勘探,2020,55(6):1271-1281.YANG Xina,ZHANG Wei. Application of CFS-PML in complex velocity models[J].Oil Geophysical Prospecting,2020,55(6):1271-1281.
[13]
陈可洋. 完全匹配层吸收边界条件研究[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.
[14]
Wang H X,Wang H X. A study of damping factors in perfectly matched layers for the numerical simulation of seismic waves[J].Applied Geophysics,2013,10(1):63-77.
[15]
Enquist B,Madja A. Absorbing boundary conditions for the numerical simulation of waves[J].Applied Mathematical Sciences,1977,74(5):1765-1766.
[16]
Clayton R,Engquist B. Absorbing boundary conditions for acoustic and elastic wave equations[J].Bulletin of the Seismological Society of America,1977,67(6):1529-1540.
[17]
Higdon R L. Absorbing boundary conditions for diffe-rence approximations to the multidimensional wave equation[J].Mathematics of Computation,1986,47(176):437-459.
[18]
Higdon R L. Numerical absorbing boundary conditions for the wave equation[J].Mathematics of Computation,1987,49(179):65-90.
[19]
Higdon R L. Absorbing boundary conditions for elastic waves[J].Geophysics,1991,56(2):231-241.
[20]
宋鹏,王修田. 优化系数的四阶吸收边界条件[J].中国海洋大学学报(自然科学版),2008,38(2):251-258.SONG Peng,WANG Xiutian. The fourth-order absorbing boundary conditions with optimal parameters[J].Periodical of Ocean University of China(Natural Science Edition),2008,38(2):251-258.
[21]
Song P,Liu Z L,Zhang X B,et al. The fourth-order absorbing boundary condition with optimized coefficients for the simulation of the acoustic equation[J].Journal of Geophysics and Engineering,2015,12(6):996-1007.
[22]
Baffet D,Hagstrom T,Givoli D. Double absorbing boundary formulations for acoustics and elastodyna-mics[J].Society for Industrial and Applied Mathema-tics,2014,36(3):A1277-A1312.
[23]
Hagstrom T,Givoli D,Rabinovich D,et al. The dou-ble absorbing boundary method[J].Journal of Computational Physics,2014,259:220-241.
[24]
Potter T,Shragge J,Lumley D. Performance and stability of the double absorbing boundary method for acoustic-wave propagation[J].Geophysics,2019,84(2):T59-T72.
[25]
Liu Y,Sen M K. A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation[J].Geophysics,2010,75(2):A1-A6.
[26]
Liu Y,Sen M K. 3D acoustic-wave modeling with time-space domain dispersion-relation-based finite-difference schemes and hybrid absorbing boundary conditions[J].Exploration Geophysics,2011,42(3):176-189.
[27]
Liu Y,Sen M K. A hybrid absorbing boundary condition for elastic staggered-grid modeling[J].Geophysical Prospecting,2012,60(6):1114-1132.
[28]
任志明,刘洋. 一阶弹性波方程数值模拟中的混合吸收边界条件[J].地球物理学报,2014,57(2):595-606.REN Zhiming,LIU Yang. Numerical modeling of the first-order elastic equations with the hybrid absorbing boundary condition[J].Chinese Journal of Geophy-sics,2014,57(2):595-606.
[29]
Liu X,Liu Y,Ren Z M,et al. Hybrid absorbing boundary condition for three-dimensional elastic wave modeling[J].Applied Geophysics,2017,14(2):270-278.
[30]
Ren Z M,Liu Y. A hybrid absorbing boundary condition for frequency-domain finite-difference modelling[J].Journal of Geophysics and Engineering,2013,10(5):054003.
[31]
刘学义,何兵寿. VTI介质中弹性波方程正演的一阶混合吸收边界[J].中国煤炭地质,2015,27(3):59-63.LIU Xueyi,HE Bingshou. First-order hybrid absor-bing boundary condition in VTI media elastic wave equation for modeling[J].Coal Geology of China,2015,27(3):59-63.
[32]
Wang E J,Liu Y. The hybrid absorbing boundary condition for one-step extrapolation and its application in wavefield decomposition-based reverse time migration[J].Journal of Geophysics and Engineering,2017,14(5):1177-1188.
[33]
薛浩,刘洋,杨宗青. 基于优化时空域频散关系的声波方程有限差分最小二乘逆时偏移[J].石油地球物理勘探,2018,53(4):745-753.XUE Hao,LIU Yang,YANG Zongqing. Least-square reverse time migration of finite-difference acoustic wave equation based on an optimal time-space dispersion relation[J].Oil Geophysical Prospecting,2018,53(4):745-753.
[34]
Liu Y,Sen M K. An improved hybrid absorbing boundary condition for wave equation modeling[J].Journal of Geophysics and Engineering,2018,15(6):2602-2613.
[35]
Levander A R. Fourth-order finite difference P-SV seismograms[J].Geophysics,1988,53(11):1425-1436.
[36]
董良国,马在田,曹景忠. 一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报,2000,43(6):856-864.DONG Liangguo,MA Zaitian,CAO Jingzhong. A study on stability of the staggered-grid high-order difference method of first-order elastic wave equation[J].Chinese Journal of Geophysics,2000,43(6):856-864.
[37]
董良国,马在田,曹景忠,等. 一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419.DONG Liangguo,MA Zaitian,CAO Jingzhong,et al. A staggered-grid high-order difference method of one-order elastic wave equation[J].Chinese Journal of Geophysics,2000,43(3):411-419.
[38]
杨庆节,刘财,耿美霞,等. 交错网格任意阶导数有限差分格式及差分系数推导[J].吉林大学学报(地球科学版),2014,44(1):375-385.YANG Qingjie,LIU Cai,GENG Meixia,et al. Staggered grid finite difference scheme and coefficients deduction of any number of derivatives[J].Journal of Jilin University(Earth Science Edition),2014,44(1):375-385.