A Novel Immersed Boundary Approach for Irregular Topography with Acoustic Wave Equations

Irregular terrain has a pronounced effect on the propagation of seismic and acoustic wavefields but is not straightforwardly reconciled with structured finite-difference (FD) methods used to model such phenomena. Methods currently detailed in the literature are generally limited in scope application-wise or non-trivial to apply to real-world geometries. With this in mind, a general immersed boundary treatment capable of imposing a range of boundary conditions in a relatively equation-agnostic manner has been developed, alongside a framework implementing this approach, intending to complement emerging code-generation paradigms. The approach is distinguished by the use of N-dimensional Taylor-series extrapolants constrained by boundary conditions imposed at some suitably-distributed set of surface points. The extrapolation process is encapsulated in modified derivative stencils applied in the vicinity of the boundary, utilizing hyperspherical support regions. This method ensures boundary representation is consistent with the FD discretization: both must be considered in tandem. Furthermore, high-dimensional and vector boundary conditions can be applied without approximation prior to discretization. A consistent methodology can thus be applied across free and rigid surfaces with the first and second-order acoustic wave equation formulations. Application to both equations is demonstrated, and numerical examples based on analytic and real-world topography implementing free and rigid surfaces in 2D and 3D are presented.

[1]  R. Gioria,et al.  spyro: a Firedrake-based wave propagation and full-waveform-inversion finite-element solver , 2022, Geoscientific Model Development.

[2]  J. Lyons,et al.  Local Explosion Detection and Infrasound Localization by Reverse Time Migration Using 3-D Finite-Difference Wave Propagation , 2021, Frontiers in Earth Science.

[3]  William J. Pringle,et al.  SeismicMesh: Triangular meshing for seismology , 2021, J. Open Source Softw..

[4]  Chengyu Sun,et al.  Mesh-free radial-basis-function-generated finite differences and their application in reverse time migration , 2020 .

[5]  G. Yao,et al.  An immersed boundary method with iterative symmetric interpolation for irregular surface topography in seismic wavefield modelling , 2020, Journal of Geophysics and Engineering.

[6]  Parashkev Nachev,et al.  Full-waveform inversion imaging of the human brain , 2019, bioRxiv.

[7]  Felix J. Herrmann,et al.  Devito: an embedded domain-specific language for finite differences and geophysical exploration , 2018, Geoscientific Model Development.

[8]  Philipp A. Witte,et al.  Architecture and Performance of Devito, a System for Automated Stencil Computation , 2018, ACM Trans. Math. Softw..

[9]  Michael Dumbser,et al.  A simple diffuse interface approach on adaptive Cartesian grids for the linear elastic wave equations with complex topography , 2018, J. Comput. Phys..

[10]  G. Yao,et al.  Accurate seabed modeling using finite difference methods , 2018, Computational Geosciences.

[11]  W. Mulder,et al.  A simple finite-difference scheme for handling topography with the first-order wave equation , 2017 .

[12]  W. Mulder A simple finite-difference scheme for handling topography with the second-order wave equation , 2017 .

[13]  Wenyi Hu,et al.  An improved immersed boundary finite-difference method for seismic wave propagation modeling with arbitrary surface topography , 2016 .

[14]  H. Maurer,et al.  Ground topography effects on near-surface elastic full waveform inversion , 2016 .

[15]  N. Anders Petersson,et al.  Wave propagation in anisotropic elastic materials and curvilinear coordinates using a summation-by-parts finite difference method , 2015, J. Comput. Phys..

[16]  Romain Brossier,et al.  An immersed free-surface boundary treatment for seismic wave simulation , 2015 .

[17]  Hermann F. Fasel,et al.  A locally stabilized immersed boundary method for the compressible Navier-Stokes equations , 2015, J. Comput. Phys..

[18]  Yi Wang,et al.  Three-dimensional full waveform inversion of short-period teleseismic wavefields based upon the SEM–DSM hybrid method , 2015 .

[19]  Bengt Fornberg,et al.  Seismic modeling with radial-basis-function-generated finite differences , 2015 .

[20]  M. Nafi Toksöz,et al.  Finite difference elastic wave modeling with an irregular free surface using ADER scheme , 2015 .

[21]  Junichi Takekawa,et al.  A mesh-free method with arbitrary-order accuracy for acoustic wave propagation , 2015, Comput. Geosci..

[22]  Takashi Furumura,et al.  Scattering of high-frequency seismic waves caused by irregular surface topography and small-scale velocity inhomogeneity , 2015 .

[23]  Andrew T. T. McRae,et al.  Firedrake , 2015, ACM Trans. Math. Softw..

[24]  Jeffrey Shragge,et al.  Solving the 3D acoustic wave equation on generalized structured meshes: A finite-difference time-domain approach , 2014 .

[25]  J. Lees,et al.  Local Volcano Infrasound and Source Localization Investigated by 3D Simulation , 2014 .

[26]  D. Fee,et al.  Network‐Based Evaluation of the Infrasonic Source Location at Sakurajima Volcano, Japan , 2014 .

[27]  Shaolin Liu,et al.  A mixed-grid finite element method with PML absorbing boundary conditions for seismic wave modelling , 2014 .

[28]  Wim A. Mulder,et al.  A comparison of continuous mass‐lumped finite elements with finite differences for 3‐D wave propagation , 2014 .

[29]  Mauricio Hanzich,et al.  Mimetic seismic wave modeling including topography on deformed staggered grids , 2014 .

[30]  William Gropp,et al.  CFD Vision 2030 Study: A Path to Revolutionary Computational Aerosciences , 2014 .

[31]  Youngseok Kim,et al.  Laplace-domain full waveform inversion using irregular finite elements for complex foothill environments , 2013 .

[32]  Hermann F. Fasel,et al.  A novel concept for the design of immersed interface methods , 2013, J. Comput. Phys..

[33]  M. Warner,et al.  Anisotropic 3D full-waveform inversion , 2013 .

[34]  Hiroshi Takenaka,et al.  FDM simulation of seismic-wave propagation for an aftershock of the 2009 Suruga bay earthquake: Effects of ocean-bottom topography and seawater layer , 2012 .

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

[36]  Richard D. Miller,et al.  An improved vacuum formulation for 2D finite-difference modeling of Rayleigh waves including surface topography and internal discontinuities , 2012 .

[37]  Romain Brossier,et al.  Shallow-structure characterization by 2D elastic full-waveform inversion , 2011 .

[38]  Jung Hee Seo,et al.  A high-order immersed boundary method for acoustic wave scattering and low-Mach number flow-induced sound in complex geometries , 2011, J. Comput. Phys..

[39]  Jean Roman,et al.  High-performance finite-element simulations of seismic wave propagation in three-dimensional nonlinear inelastic geological media , 2010, Parallel Comput..

[40]  G. Lauder,et al.  Computational modelling and analysis of the hydrodynamics of a highly deformable fish pectoral fin , 2010, Journal of Fluid Mechanics.

[41]  Yang Liu,et al.  A new time-space domain high-order finite-difference method for the acoustic wave equation , 2009, J. Comput. Phys..

[42]  Stéphane Rondenay,et al.  Effects of surface scattering in full-waveform inversion , 2009 .

[43]  William W. Symes,et al.  Interface error analysis for numerical wave propagation , 2009 .

[44]  Martin Galis,et al.  A 3-D hybrid finite-difference—finite-element viscoelastic modelling of seismic wave motion , 2008 .

[45]  Haibo Dong,et al.  A computational study of the aerodynamic performance of a dragonfly wing section in gliding flight , 2008, Bioinspiration & biomimetics.

[46]  Rajat Mittal,et al.  A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries , 2008, J. Comput. Phys..

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

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

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

[50]  H. Igel,et al.  Seismic wave simulation in the presence of real volcano topography , 2003 .

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

[52]  S. Hestholm,et al.  2D surface topography boundary conditions in seismic wave modelling , 2001 .

[53]  M. Lai,et al.  An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced Numerical Viscosity , 2000 .

[54]  E. Husebye,et al.  3-D versus 2-D finite-difference seismic synthetics including real surface topography , 1999 .

[55]  I. Opršal,et al.  Elastic finite-difference method for irregular grids , 1999 .

[56]  Bent O. Ruud,et al.  3-D finite-difference elastic wave modeling including surface topography , 1998 .

[57]  Craig A. Schultz,et al.  A density‐tapering approach for modeling the seismic response of free‐surface topography , 1997 .

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

[59]  Johan O. A. Robertsson,et al.  A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography , 1996 .

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

[61]  C. Tam,et al.  Dispersion-relation-preserving finite difference schemes for computational acoustics , 1993 .

[62]  L. Sirovich,et al.  Modeling a no-slip flow boundary with an external force field , 1993 .

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

[64]  F. Muir,et al.  Modeling elastic fields across irregular boundaries , 1992 .

[65]  G. Bollinger,et al.  The effect of Appalachian Mountain topography on seismic waves , 1979, Bulletin of the Seismological Society of America.

[66]  David M. Boore,et al.  The effect of simple topography on seismic waves: Implications for the accelerations recorded at Pacoima Dam, San Fernando Valley, California , 1973, Bulletin of the Seismological Society of America.

[67]  D. Boore A note on the effect of simple topography on seismic SH waves , 1972, Bulletin of the Seismological Society of America.

[68]  G. Yao,et al.  Waveform inversion of seismic first arrivals acquired on irregular surface , 2022 .

[69]  J. Tromp,et al.  3D elastic full-waveform inversion of surface waves in the presence of irregular topography using an envelope-based misfit function , 2018 .

[70]  R. Shamasundar,et al.  Delft University of Technology Performance of continuous mass-lumped tetrahedral elements for elastic wave propagation with and without global assembly , 2016 .

[71]  Dimitri Komatitsch,et al.  A hybrid method to compute short-period synthetic seismograms of teleseismic body waves in a 3-D regional model , 2013 .

[72]  G. Schuster,et al.  Multi-source Least-squares Reverse Time Migration with Topography , 2013 .

[73]  D. H. Hien,et al.  Reverse time migration for tilted transversely isotropic (TTI) media , 2011 .

[74]  L. Bartel,et al.  Graded Boundary Simulation of Air/earth Interfaces In Finite-difference Elastic Wave Modeling , 2000 .

[75]  Eduardo Reinoso,et al.  Three-dimensional scattering of seismic waves from topographical structures , 1997 .