Abstract:Numerical simulation of wave propagation in heterogeneous medium provides a fundamental tool in many inverse scattering problems. Although real data are acquired in three-dimensional (3D) world, we can still apply some approximation to simplify the model and to accelerate the computation. For example, when the data are acquired along a straight line, they are mainly sensitive to the formation in the vertical plane of the acquisition line. Therefore, we can use a two-dimensional (2D) model as an approximation while maintaining the 3D acquisition layout. This defines the settings of the two-and-a-half dimensional (2.5D) modeling schemes. Here we study a two-and-a-half dimensional (2.5D) finite-difference algorithm to model elastic wave propagation in heterogeneous media. In 2.5D problems, it is assumed that the elastic properties of models are invariant along a certain direction. Therefore, we can convert the three-dimensional (3D) problem into a set of two-dimensional (2D) problems in the spectral domain. The 3D solutions are then obtained by applying a numerical integration in the spectral domain. Usually the quadrature points used in this spectral-domain integration are sampled from the real axis in the spectral domain. However, because of the oscillating nature of the integrand, the convergence of this quadrature can be very slow especially at high frequencies. In this work, we apply the near-optimal quadrature scheme to the spectral integration. This is equivalent to transform the contour of integration from the real axis into a path in the complex plane and near-optimal quadrature is computed along the path numerically. Numerical studies show more than ten times of reduction in the number of quadrature points compared with sampling along the real axis in the spectral domain. This scheme alleviates the bottleneck of computing speed in the 2.5D elastic wave modeling. Furthermore, it can improve the computational efficiency of 2.5D elastic full waveform inversion algorithms. This technique can also be applied to solve 2.5D electromagnetic problems.

DPSM Modeling of Wave Propagation in Anisotropic Half Space

NG-DPSM: A neural green-distributed point source method for modelling ultrasonic field emission near fluid-solid interface using physics informed neural network

Coupled Finite Particle Method for Simulations of Wave and Structure Interaction

Active Sourced Wavefield Modeling for Layered Half-Space

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

A Finite Volume Method for Elastic Waves in Heterogeneous, Anisotropic and Fractured Media

Introduction to the smoothed particle hydrodynamics method in electromagnetics

Preliminary Understanding of Near-Field Effect of CSAMT in Electric Azimuthally Anisotropic Half-Space

An Accurate and Efficient Time-domain Model for Simulating Wave Propagation in Half-space

Application of the Near-Optimal Quadrature to 2.5d Wave Modeling

Efficient Calculation of Electromagnetic Characteristics of 2-D Periodic Dielectric Objects Above Half-Space

A finite element method for modelling electromechanical wave propagation in anisotropic piezoelectric media

1D Time-Domain Method for 2D Wave Motion in Elastic Half-Space by Antiplane Wave Oblique Incidence

Generalized Transfer Matrix Method for Propagation of Surface Waves in Layered Azimuthally Anisotropic Half-Space

Multi-physics and Multi-scale Electromagnetic Modeling and High Performance Algorithms

Symplectic Stereomodelling Method for Solving Elastic Wave Equations in Porous Media

Fully Coupled Electromechanical Elastodynamic Model for Guided Wave Propagation Analysis

The Meshfree Radial Point Interpolation Method (RPIM) for Wave Propagation Dynamics in Non-Homogeneous Media

Time Domain Green Functions Technique for a Point Source over a Dissipative Stratified Half-Space with a Phase Velocity Mismatch at the Surface

Modeling 3-D Elastic Wave Propagation in TI Media Using Discontinuous Galerkin Method on Tetrahedral Meshes

Time domain Green function technique for a point source over a dissipative stratified half-space