Nonsplit complex-frequency-shifted perfectly matched layer combined with symplectic methods for solving second-order seismic wave equations — Part 2: Wavefield simulations

Abstract

The perfectly matched layer (PML) is an efficient artificial boundary condition that has been routinely implemented in seismic wave modeling. However, the effective combination of PML with symplectic numerical schemes for solving seismic wave equations has rarely been studied. In a companion paper, we have developed a complex-frequency-shifted convolutional PML (CPML) with a nonconstant compression grid parameter for solving the time-domain second-order seismic wave equation. Subsequently, we combine this CPML with two classes of symplectic methods to formulate symplectic partitioned Runge-Kutta (SPRK) + CPML and nearly analytic SPRK (NSPRK) + CPML, both of which are properly synchronized. To further investigate their validity, the two algorithms are then applied to acoustic and elastic wave simulations in typical geologic models, including a heterogeneous acoustic model, several isotropic and orthotropic elastic models, and an isotropic elastic model with a free-surface boundary. Relevant numerical results demonstrate the effectiveness of our CPML and combination algorithms. Specifically, the numerical accuracy and stability of the CPML that we develop are greatly improved compared with the classic split-field PML. Moreover, the final model with the free-surface boundary condition indicates that the nonconstant grid-compression parameter can eliminate the unstable modes at the free surface in the PML domain. The (N)SPRK + CPML that we propose is prospective for future application in other complex models and wave-equation-based migration and inversion.