Abstract: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. Accurate wavefield simulation is paramount in subsurface imaging applications such as reverse-time migration (RTM) and full-waveform inversion (FWI), requiring suitable topography handling. Methods currently detailed in the literature are generally limited in scope application-wise or non-trivial to apply to real-world geometries. 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 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. 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 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. Numerical examples and convergence tests demonstrate the accuracy of boundary treatments devised by the prescribed approach, their suitability to practical applications, and the feasibility of automatically generating treatments to suit each case.

RBF-FD analysis of 2D time-domain acoustic wave propagation in heterogeneous media

2D Frequency-Domain Finite-Difference Acoustic Wave Modeling Using Optimized Perfectly Matched Layers

Application of selective filtering non-staggered finite difference method to seismic wave simulation

Oscillatory Responses of a Spatially Random Seabed in the Vicinity of Structures: Validation and Application of the Integrated CFD-SFEM

Distributional finite-difference modeling of seismic waves

A stabilized radial basis-finite difference (RBF-FD) method with hybrid kernels

Computational Electromagnetics with the RBF-FD Method

An unfitted RBF-FD method in a least-squares setting for elliptic PDEs on complex geometries

Application of an RBF-FD solver for the Helmholtz equation to full-waveform inversion

Radial basis function collocation method for decoupled fractional Laplacian wave equations

SIMULATION OF ACOUSTIC WAVE PROPAGATION IN COMPLEX MEDIA USING MRFD METHOD

A space-time domain RBF method for 2D wave equations

An Efficient Immersed Free Surface Boundary Method for 3-D Scalar Seismic Waves Finite-Difference Modeling in Presence of Topography

A Compact High-Order Finite-Difference Method with Optimized Coefficients for 2D Acoustic Wave Equation

2-D seismic wave propagation using the distributional finite-difference method: further developments and potential for global seismology

FDwave3D: a MATLAB solver for the 3D anisotropic wave equation using the finite-difference method

Simulate the elastic wavefields in media with an irregular surface topography based on staggered grid finite difference

A superconvergence result in the RBF-FD method

A Nearly Analytic Discrete Method for Acoustic and Elastic Wave Equations in Anisotropic Media

A Runge-Kutta method using high-order interpolation approximation for solving 2D acoustic and elastic wave equations

A novel immersed boundary approach for irregular topography with acoustic wave equations