SBP-SAT finite difference discretization of acoustic wave equations on staggered block-wise uniform grids

We consider the numerical simulation of the acoustic wave equations arising from seismic applications, for which staggered grid finite difference methods are popular choices due to their simplicity and efficiency. We relax the uniform grid restriction on finite difference methods and allow the grids to be block-wise uniform with nonconforming interfaces. In doing so, variations in the wave speeds of the subterranean media can be accounted for more efficiently. Staggered grid finite difference operators satisfying the summation-by-parts (SBP) property are devised to approximate the spatial derivatives appearing in the acoustic wave equation. These operators are applied within each block independently. The coupling between blocks is achieved through simultaneous approximation terms (SATs), which impose the interface condition weakly, i.e., by penalty. Ratio of the grid spacing of neighboring blocks is allowed to be rational number, for which specially designed interpolation formulas are presented. These interpolation formulas constitute key pieces of the simultaneous approximation terms. The overall discretization is shown to be energy-conserving and examined on test cases of both theoretical and practical interests, delivering accurate and stable simulation results.

[1]  Ken Mattsson,et al.  Boundary Procedures for Summation-by-Parts Operators , 2003, J. Sci. Comput..

[2]  Wei Zhang,et al.  Stable discontinuous grid implementation for collocated-grid finite-difference seismic wave modelling , 2013 .

[3]  Martin Galis,et al.  Stable discontinuous staggered grid in the finite-difference modelling of seismic motion , 2010 .

[4]  D. Komatitsch,et al.  The spectral element method: An efficient tool to simulate the seismic response of 2D and 3D geological structures , 1998, Bulletin of the Seismological Society of America.

[5]  Roelof Versteeg,et al.  The Marmousi experience; velocity model determination on a synthetic complex data set , 1994 .

[6]  Gerard T. Schuster,et al.  Wave-equation traveltime inversion , 1991 .

[7]  Siyang Wang,et al.  High Order Finite Difference Methods for the Wave Equation with Non-conforming Grid Interfaces , 2015, Journal of Scientific Computing.

[8]  Michael Dumbser,et al.  An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes - IV. Anisotropy , 2007 .

[9]  C. Loan.  The ubiquitous Kronecker product , 2000 .

[10]  N. Anders Petersson,et al.  Treatment of the polar coordinate singularity in axisymmetric wave propagation using high-order summation-by-parts operators on a staggered grid , 2017 .

[11]  D. Komatitsch,et al.  Introduction to the spectral element method for three-dimensional seismic wave propagation , 1999 .

[12]  Magnus Svärd,et al.  Review of summation-by-parts schemes for initial-boundary-value problems , 2013, J. Comput. Phys..

[13]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[14]  Bertil Gustafsson,et al.  The convergence rate for difference approximations to general mixed initial boundary value problems , 1981 .

[15]  Jeremy E. Kozdon,et al.  Boundary conditions and stability of a perfectly matched layer for the elastic wave equation in first order form , 2015, J. Comput. Phys..

[16]  Jean Virieux,et al.  Using time filtering to control the long-time instability in seismic wave simulation , 2016 .

[17]  M. Nafi Toksöz,et al.  Discontinuous-Grid Finite-Difference Seismic Modeling Including Surface Topography , 2001 .

[18]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[19]  William W. Symes,et al.  Migration velocity analysis and waveform inversion , 2008 .

[20]  M. Dumbser,et al.  An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes - I. The two-dimensional isotropic case with external source terms , 2006 .

[21]  David E. Keyes,et al.  On long-time instabilities in staggered finite difference simulations of the seismic acoustic wave equations on discontinuous grids , 2018 .

[22]  D. Komatitsch,et al.  An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation , 2007 .

[23]  Magnus Svärd,et al.  Steady-State Computations Using Summation-by-Parts Operators , 2005, J. Sci. Comput..

[24]  Ken Mattsson,et al.  Boundary optimized diagonal-norm SBP operators , 2018, J. Comput. Phys..

[25]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .

[26]  B. Gustafsson.  The convergence rate for difference approximations to mixed initial boundary value problems , 1975 .

[27]  H. Kreiss,et al.  Finite Element and Finite Difference Methods for Hyperbolic Partial Differential Equations , 1974 .

[28]  Jason E. Hicken,et al.  Summation-by-parts operators and high-order quadrature , 2011, J. Comput. Appl. Math..

[29]  Mark H. Carpenter,et al.  Stable and Accurate Interpolation Operators for High-Order Multiblock Finite Difference Methods , 2009, SIAM J. Sci. Comput..

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

[31]  Jean Virieux,et al.  SH-wave propagation in heterogeneous media: velocity-stress finite-difference method , 1984 .

[32]  A. Tarantola.  Inversion of seismic reflection data in the acoustic approximation , 1984 .

[33]  Erik Schnetter,et al.  Optimized High-Order Derivative and Dissipation Operators Satisfying Summation by Parts, and Applications in Three-dimensional Multi-block Evolutions , 2005, J. Sci. Comput..

[34]  Jan Nordström,et al.  Energy stable and high-order-accurate finite difference methods on staggered grids , 2017, J. Comput. Phys..

[35]  Kenneth Duru,et al.  The Role of Numerical Boundary Procedures in the Stability of Perfectly Matched Layers , 2014, SIAM J. Sci. Comput..

[36]  Jean Virieux,et al.  An overview of full-waveform inversion in exploration geophysics , 2009 .

[37]  Daniel J. Bodony,et al.  Accuracy of the Simultaneous-Approximation-Term Boundary Condition for Time-Dependent Problems , 2010, J. Sci. Comput..

[38]  David C. Del Rey Fernández,et al.  Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations , 2014 .

[39]  N. Anders Petersson,et al.  Stable Grid Refinement and Singular Source Discretization for Seismic Wave Simulations , 2009 .

[40]  Ken Mattsson,et al.  Compatible diagonal-norm staggered and upwind SBP operators , 2018, J. Comput. Phys..

[41]  A. Levander.  Fourth-order finite-difference P-SV seismograms , 1988 .

[42]  Willi-Hans Steeb,et al.  Matrix Calculus and the Kronecker Product with Applications and C++ Programs , 1997 .