A comparative study of Mur second-order absorbing boundary condition and unsplit recursive convolutional perfectly matched layer method under multi-source concurrency
CUI Fan1,2(), CHEN Yi1(), XUE Han-Peng1, PENG Su-Ping1,2, DU Yun-Fei1
1. School of Geosciences and Surveying Engineering,China University of Mining and Technology(Beijing),Beijing 100083,China 2. State Key Laboratory of Coal Resources and Safe Mining,China University of Mining and Technology(Beijing),Beijing 100083,China
Plane beam signals form when multiple excitation sources simultaneously emit pulses with the same center frequency (multiple-source concurrency),thus enhancing the quality of data records.This paper compares and analyzes the electromagnetic wave absorption effects of unsplit recursive convolutional perfectly matched layer (PML) as the absorbing boundary condition and Mur second-order absorbing boundary condition under multi-source concurrency through numerical simulation.According to study results,the traditional Mur second-order absorbing boundary condition did not perform well in absorbing electromagnetic waves under the conditions of multi-source concurrency and multi-angle grazing,and it will cause waveform distortion and spurious reflections in the case of large offsets.For the unsplit recursive convolutional perfectly matched layer as the absorbing boundary condition under multi-source concurrency,coordinate scale factors were introduced into the finite-difference time-domain (FDTD) algorithm.Then,the PML equation for coordinate stretching was transformed from frequency domain into time domain through the inverse Fourier transform.Finally,the electric and magnetic field values were solved using the recursive convolution method in the discrete state,thus avoiding the complicated calculation involved in directly determining the numerical solution of convolution.This allows less memory space and high calculation efficiency while ensuring accuracy.Therefore,the unsplit recursive convolutional perfectly matched layer method improves the electromagnetic wave absorption effect at the positions where the grid terminate without inducing wave-field splitting.
Bayliss A, Turkel E. Radiation boundary conditions for wave-like equations[J]. Communications on Pure and Applied Mathematics, 1980, 33(6):707-725.
Engquist B, Majda A. Absorbing boundary conditions for numerical simulation of waves[J]. Proceedings of the National Academy of Sciences of the United States of America, 1977, 74(5):1765-1766.
Mur G. Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic-field equations[J]. IEEE Transactions on Electromagnetic Compatibility, 1981, 23(4):377-382.
Berenger J P. Three-dimensional perfectly matched layer for the absorption of electromagnetic waves[J]. Journal of Computational Physics, 1996, 127(2):363-379.
Chen Y H, Weng C C, Oristagilo M L. Application of perfectly matched layers to the transient modeling of subsurface EM problems[J]. Geophysics, 1997, 97(6):1730-1736.
Festa G. Interaction between surface waves and absorbing boundaries for wave propagation in geological basins:2D numerical simulations[J]. Geophysical Research Letters, 2005, 32:L20306.
Qin Z, Lu M, Zheng X, et al. The implementation of an improved NPML absorbing boundary condition in elastic wave modeling[J]. Applied Geophysics, 2009, 6(2):113-121.
Higdon R L. Absorbing boundary conditions for elastic waves[J]. Siam Journal on Numerical Analysis, 2012, 31(2):64-100.
Liao Z P, Wong H L, Yang B P, et al. A transmitting boundary for transient wave analysis[J]. Science in China Series A-Mathematics Physics Astronomy & Technological Science, 1984, 27(10):1063-1076.
Marfurt K J. Accuracy of finite difference and finite element modeling of the scalar and elastic wave equation[J]. Geophysics, 1984, 49(4):533-549.
Feng D S, Xie Y. Full Wave field forward modeling of 3D GPR based on boundary conditions of uniaxial anisotropic completely matched layer[J]. Journal of Central South University:Science and Technology, 2011, 42(8):202-210.
Feng D S, Chen C S, Dai Q W. Numerical simulation of GPR full wave field with alternate direction implicit finite difference method based on UPML boundary condition[J]. Chinese Journal of Geophysics, 2010, 53(10):2484-2496.
Feng D S, Wang X Y. Elastic wave propagation simulation in anisotropic media and random media using heig-order difference method of rotation staggered grids based on convolutional perfectly matched layer[J]. Geophysical and Geochemical Exploration, 2018, 42(4):766-776.