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.

Higher-order 3D Parabolic Equation for Radio Wave Propagation Modeling in a Street Canyon

Application of a 3D Parabolic Equation Model with Higher-Order Operator Approximation

Numerical Simulation and Analytical Solution for Plane Waves Focused by a Parabolic Reflector

Modeling Processes of Radio Wave Propagation in Urban Areas

Modeling Of Rayleigh-Wave Propagation In Three-Dimensional Media With Topography Using Staggered-Grid Finite Difference Scheme

A 3D Parabolic Equation Model with Higher-Order Operator Approximation

Numerical Investigation of Rayleigh-wave Propagation on Topography Surface

Radio Wave Propagation Prediction in Urban Based on Parabolic Equation Methods

A Radio Wave Propagation Modeling Method Based on High-Precision 3-D Mapping in Urban Scenarios

Finite-Difference Solution of the Parabolic Equation Under Horizontal Polar Coordinates

Application of a spectral scheme to simulate horizontally slowly varying three-dimensional ocean acoustic propagation

Secure Counting: counting members of a subset without revealing their identities

Stability And Accuracy Of The Cubic Interpolated Propagation Scheme

A versatile wave propagation model for the VHF/UHF range considering three-dimensional terrain

A graph theoretic method for modeling urban scenarios radio wave propagation

A vector Helmholtz electromagnetic wave propagator for inhomogeneous media

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

Concepts and results for 3D digital terrain-based wave propagation models: an overview

A Ray-Tracing Method Based on the Triangular Grid Approach and Application to Propagation Prediction in Urban Environments

Three-dimensional boundary fitted parabolic-equation model of underwater sound propagation

A New Series Solution Method for Two-Dimensional Elastic Wave Scattering along A Canyon in Half-Space