Embedded boundary methods for modeling 3D finite-difference Laplace-Fourier domain acoustic-wave equation with free-surface topography

Author(s): Alsalem, H; Petrov, P; Newman, G; Rector, J | Abstract: © 2018 Society of Exploration Geophysicists. We have developed embedded boundary methods to handle arbitrarily shaped topography to accurately simulate acoustic seismic wave propagation in the Laplace-Fourier domain. The purpose is to use this method to enhance accurate wave simulation near the surface. Unlike most existing methods such as the ones using curvilinear grids to fit irregular surface topography, we use a regular Cartesian grid system without suffering from the staircasing error that occurs in conventional implementations. In this improved embedded-boundary method, we use the method of images, by imposing ghost nodes above the surface and approximating their acoustic pressures using linear extrapolation, quadratic interpolation, or cubic interpolation, to account for an arbitrarily curved surface. Implementing this method instead of using curvilinear grids near the boundaries greatly reduces the complexity of preprocessing procedures and the computational cost. Furthermore, using numerical examples, we found the accuracy gain and performance of our embedded-boundary methods in comparison with conventional finite-difference implementation of the problem.

[1]  P. Morse,et al.  Methods of theoretical physics , 1955 .

[2]  K. Marfurt.  Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations , 1984 .

[3]  Martin B. van Gijzen,et al.  IDR(s): A Family of Simple and Fast Algorithms for Solving Large Nonsymmetric Systems of Linear Equations , 2008, SIAM J. Sci. Comput..

[4]  S. Operto,et al.  Mixed‐grid and staggered‐grid finite‐difference methods for frequency‐domain acoustic wave modelling , 2004 .

[5]  Alfred Behle,et al.  Elastic wave propagation simulation in the presence of surface topography , 1992 .

[6]  Maik Moeller,et al.  Introduction to Electrodynamics , 2017 .

[7]  Effect of a simple mountain range on underground seismic motion , 1984 .

[8]  Wei Zhang,et al.  Traction image method for irregular free surface boundaries in finite difference seismic wave simulation , 2006 .

[9]  M. Spivak.  A comprehensive introduction to differential geometry , 1979 .

[10]  Zhang Jianfeng,et al.  P–SV-wave propagation in heterogeneous media: grid method , 1999 .

[11]  Peter Moczo,et al.  Testing four elastic finite-difference schemes for behavior at discontinuities , 1993 .

[12]  Ekkehart Tessmer,et al.  3-D Elastic Modeling With Surface Topography By a Chebychev Spectral Method , 1994 .

[13]  S. Operto,et al.  3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study , 2007 .

[14]  Stig Hestholm,et al.  Three-dimensional finite difference viscoelastic wave modelling including surface topography , 1999 .

[15]  R. Pratt.  Inverse theory applied to multisource cross-hole tomography, Part2 : Elastic wave-equation method , 1990 .

[16]  John Lysmer,et al.  A Finite Element Method for Seismology , 1972 .

[17]  C. Shin,et al.  An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator , 1996 .

[18]  Erik H. Saenger,et al.  Accuracy of heterogeneous staggered-grid finite-difference modeling of Rayleigh waves , 2006 .

[19]  John B. Schneider,et al.  Application of the perfectly matched layer (PML) absorbing boundary condition to elastic wave propagation , 1996 .

[20]  I. Štekl,et al.  Accurate viscoelastic modeling by frequency‐domain finite differences using rotated operators , 1998 .

[21]  A. Fichtner.  Full Seismic Waveform Modelling and Inversion , 2011 .

[22]  Robert W. Graves,et al.  Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences , 1996, Bulletin of the Seismological Society of America.

[23]  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.

[24]  Yang Zhang,et al.  Frequency-Domain Finite-Difference Acoustic Modeling With Free Surface Topography Using Embedded Boundary Method , 2010 .

[25]  Jose Pujol,et al.  Elastic Wave Propagation and Generation in Seismology: Anelastic attenuation , 2003 .

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

[27]  Werner C. Rheinboldt,et al.  Methods for solving systems of nonlinear equations , 1987 .

[28]  Jean Virieux,et al.  Free and smooth boundaries in 2-D finite-difference schemes for transient elastic waves , 2007, 0706.3825.

[29]  D. A. Dunnett.  Classical Electrodynamics , 2020, Nature.

[30]  Gregory A. Newman,et al.  3D finite-difference modeling of elastic wave propagation in the Laplace-Fourier domain , 2012 .

[31]  Wei Zhang,et al.  Three-dimensional elastic wave numerical modelling in the presence of surface topography by a collocated-grid finite-difference method on curvilinear grids , 2012 .

[32]  Heinz-Otto Kreiss,et al.  An Embedded Boundary Method for the Wave Equation with Discontinuous Coefficients , 2005, SIAM J. Sci. Comput..

[33]  Werner C. Rheinboldt,et al.  Methods for Solving Systems of Nonlinear Equations: Second Edition , 1998 .

[34]  Bent O. Ruud,et al.  2D finite‐difference viscoelastic wave modelling including surface topography , 2000 .

[35]  Ekkehart Tessmer,et al.  3-D elastic modeling with surface topography by a Chebychev spectral method , 1994 .

[36]  Joseph E. Pasciak,et al.  Analysis of a Cartesian PML approximation to acoustic scattering problems in R2 , 2010 .

[37]  Ellen Gottschämmer,et al.  Accuracy of the Explicit Planar Free-Surface Boundary Condition Implemented in a Fourth-Order Staggered-Grid Velocity-Stress Finite-Difference Scheme , 2001 .