A simple implementation of PML for second-order elastic wave equations

Abstract When modeling time-domain elastic wave propagation in an unbound space, the standard perfectly matched layer (PML) is straightforward for the first-order partial differential equations (PDEs); by contrast, the PML requires tremendous re-constructions of the governing equations in the second-order PDE form, which is however preferable, because of much less memory and time consumption. Therefore, it is imperative to explore a simple implementation of PML for the second-order system. In this work, we first systematically extend the first-order Nearly PML (NPML) technique into second-order systems, implemented by the spectral element and finite difference time-domain algorithms. It merits the following advantages: the simplicity in implementation, by keeping the second-order PDE-based governing equations exactly the same; the efficiency in computation, by introducing a set of auxiliary ordinary differential equations (ODEs). Mathematically, this PML technique effectively hybridizes the second-order PDEs and first-order ODEs, and locally attenuates outgoing waves, thus efficiently avoid either spatial or temporal global convolutions. Numerical experiments demonstrate that the NPML for the second-order PDE has an excellent absorbing performance for elastic, anelastic and anisotropic media in terms of the absorption accuracy, implementation complexity, and computation efficiency.

[1]  A. Chopra,et al.  Perfectly matched layers for transient elastodynamics of unbounded domains , 2004 .

[2]  Jingyi Chen,et al.  Application of the Nearly Perfectly Matched Layer to Seismic-Wave Propagation Modeling in Elastic Anisotropic Media , 2011 .

[3]  T. Mukerji,et al.  The Rock Physics Handbook , 1998 .

[4]  Qing Huo Liu,et al.  Isotropic Riemann Solver for a Nonconformal Discontinuous Galerkin Pseudospectral Time-Domain Algorithm , 2017, IEEE Transactions on Geoscience and Remote Sensing.

[5]  Barbara Kaltenbacher,et al.  A modified and stable version of a perfectly matched layer technique for the 3-d second order wave equation in time domain with an application to aeroacoustics , 2013, J. Comput. Phys..

[6]  D. Komatitsch,et al.  Simulation of anisotropic wave propagation based upon a spectral element method , 2000 .

[7]  Bernard A. Chouet,et al.  A free-surface boundary condition for including 3D topography in the finite-difference method , 1997, Bulletin of the Seismological Society of America.

[8]  Jeroen Tromp,et al.  A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation , 2003 .

[9]  Qing Huo Liu,et al.  The rotated Cartesian coordinate method to remove the axial singularity of cylindrical coordinates in finite‐difference schemes for elastic and viscoelastic waves , 2018 .

[10]  Roland Martin,et al.  Improved forward wave propagation and adjoint-based sensitivity kernel calculations using a numerically stable finite-element PML , 2014 .

[11]  Yu Zhang,et al.  A multiaxial perfectly matched layer (M-PML) for the long-time simulation of elastic wave propagation in the second-order equations , 2014 .

[12]  José M. Carcione,et al.  Seismic modeling in viscoelastic media , 1993 .

[13]  J. Vilotte,et al.  The Newmark scheme as velocity–stress time-staggering: an efficient PML implementation for spectral element simulations of elastodynamics , 2005 .

[14]  Carlos Torres-Verdín,et al.  Back-propagating modes in elastic logging-while-drilling collars and their effect on PML stability , 2013, Comput. Math. Appl..

[15]  M. Korn,et al.  Incorporation of attenuation into time-domain computations of seismic wave fields , 1987 .

[16]  Qing Huo Liu,et al.  Perfectly matched layers for elastic waves in cylindrical and spherical coordinates , 1999 .

[17]  Guo Tao,et al.  Stability of finite difference numerical simulations of acoustic logging-while-drilling with different perfectly matched layer schemes , 2013, Applied Geophysics.

[18]  Jingyi Chen.  Application of the nearly perfectly matched layer for seismic wave propagation in 2D homogeneous isotropic media , 2011 .

[19]  Ladislav Halada,et al.  3D Fourth-Order Staggered-Grid Finite-Difference Schemes: Stability and Grid Dispersion , 2000 .

[20]  Ushnish Basu,et al.  Explicit finite element perfectly matched layer for transient three‐dimensional elastic waves , 2009 .

[21]  Don L. Anderson,et al.  Velocity dispersion due to anelasticity; implications for seismology and mantle composition , 1976 .

[22]  Erik H. Saenger,et al.  Finite-difference modeling of viscoelastic and anisotropic wave propagation using the rotated staggered grid , 2004 .

[23]  Qing Huo Liu,et al.  A three‐dimensional finite difference simulation of sonic logging , 1996 .

[24]  Kristel C. Meza-Fajardo,et al.  A Nonconvolutional, Split-Field, Perfectly Matched Layer for Wave Propagation in Isotropic and Anisotropic Elastic Media: Stability Analysis , 2008 .

[25]  Roland Martin,et al.  A perfectly matched layer for fluid-solid problems: Application to ocean-acoustics simulations with solid ocean bottoms. , 2016, The Journal of the Acoustical Society of America.

[26]  Qing Huo Liu,et al.  Incorporating Full Attenuation Mechanisms of Poroelastic Media for Realistic Subsurface Sensing , 2019, IEEE Transactions on Geoscience and Remote Sensing.

[27]  René Matzen.  An efficient finite element time‐domain formulation for the elastic second‐order wave equation: A non‐split complex frequency shifted convolutional PML , 2011 .