TY - JOUR
T1 - Failures of the perfectly-matched layer method in frequency-domain seismic wave modelling in elastic anisotropic media
AU - Zhou, Bing
AU - Greenhalgh, Stewart
AU - Liu, Xu
AU - Bouzidi, Youcef
AU - Riahi, Mohamed Kamel
AU - Al-Khaleel, Mohammad
N1 - Funding Information:
This work is supported by the Khalifa University of Science and Technology under Award No. CIRA-2018-48.
Publisher Copyright:
© 2019 SEG
PY - 2019/8/10
Y1 - 2019/8/10
N2 - The perfectly-matched layer (PML) technique is popular approach to remove the artificial edge effects in numerical wave modelling. However, there are stability conditions of the PML which all wavefronts must satisfy, such as the spherical wavefronts of P and S waves in elastic isotropic media. But the qSV wavefronts in anisotropic media may present severe instability of the PML because of triplications or cusps. In this paper, we give examples showing the failures of the PML in frequency-domain seismic wave modelling and how serious the problem can be. Our results demonstrate that the frequency-domain wave solutions may be destroyed by intersecting qSV wavefronts in singular directions. To overcome this issue, we introduce a generalized stiffness reduction method (GSRM) that aims at stabilizing the absorbing layer. We show that the GSRM successfully dissipates the wave energy to such an extent at the edge of the grid so that the spurious boundary reflections are minimized. Our numerical examples prove that the GSRM can replace the PML in frequency-domain seismic wave modelling in arbitrary elastic anisotropic media.
AB - The perfectly-matched layer (PML) technique is popular approach to remove the artificial edge effects in numerical wave modelling. However, there are stability conditions of the PML which all wavefronts must satisfy, such as the spherical wavefronts of P and S waves in elastic isotropic media. But the qSV wavefronts in anisotropic media may present severe instability of the PML because of triplications or cusps. In this paper, we give examples showing the failures of the PML in frequency-domain seismic wave modelling and how serious the problem can be. Our results demonstrate that the frequency-domain wave solutions may be destroyed by intersecting qSV wavefronts in singular directions. To overcome this issue, we introduce a generalized stiffness reduction method (GSRM) that aims at stabilizing the absorbing layer. We show that the GSRM successfully dissipates the wave energy to such an extent at the edge of the grid so that the spurious boundary reflections are minimized. Our numerical examples prove that the GSRM can replace the PML in frequency-domain seismic wave modelling in arbitrary elastic anisotropic media.
UR - http://www.scopus.com/inward/record.url?scp=85121864523&partnerID=8YFLogxK
U2 - 10.1190/segam2019-3209086.1
DO - 10.1190/segam2019-3209086.1
M3 - Conference article
AN - SCOPUS:85121864523
SN - 1052-3812
SP - 3750
EP - 3754
JO - SEG Technical Program Expanded Abstracts
JF - SEG Technical Program Expanded Abstracts
T2 - Society of Exploration Geophysicists International Exposition and 89th Annual Meeting, SEG 2019
Y2 - 15 September 2019 through 20 September 2019
ER -