Abstract:Seismic imaging necessitates precisely separating P and S waves to mitigate undesirable crosstalk between them. Failure to properly handle this crosstalk can lead to distorted or misinterpreted images. In an elastic anisotropic medium, the polarization of a P wave is not necessarily parallel to its propagating direction, making the separation of P and S waves more difficult than in the case of an elastic isotropic medium. A pure acoustic transversely isotropic (TI) wave equation involves a mixed-domain problem. This usually means the space dependencies and wavenumbers are coupled in the wavefield extrapolator. An efficient solution to the wavefield requires a separation between the space dependencies and wavenumbers in the extrapolator. There are two primary approaches: (1) an analytical separation of the mixed-domain term; and (2) a numerical solution to the mixed-domain term, such as low-rank and Poynting vector methods. In this study, the acoustic approximation is derived by setting the S-wave phase velocity to zero in all phase directions. An analytical separation is developed in this study by employing a Taylor series to linearize the anellipticity term, thereby decoupling the space-dependent and wavenumber-dependent variables within the wavefield extrapolator. The kinematics analysis demonstrates that our method has competitive phase velocity accuracy compared with existing methods. Moreover, the corresponding TI acoustic wave equation can be efficiently solved by the hybrid finite-difference/pseudospectral method. In a 2D case, it only requires one fast Fourier transform, two inverse Fourier transforms, and a few additional spatial differentiations at each time step. Remarkably, our wavefield solution remains independent of the wavefield heterogeneity (as observed in the low-rank approach) and is immune to wavefield interference (as encountered in the Poynting vector method). The developed acoustic approximation ensures the absence of any S-wave artifacts with relatively cheap computation. Numerical examples indicate our acoustic tilted transversely isotropic wave equation is reliable for accurate wavefield modeling and reverse time migration.

An anti-dispersion wave equation for modeling and reverse-time migration

An Optimized Wave Equation for Seismic Modeling and Reverse Time Migration

Elastic Reverse Time Migration Without Wavefield Decomposition

Wave Equation Migration for 3-D VSP Data

Isotropic Elastic Reverse Time Migration

Dispersion-Relation Preserving Stereo-Modeling Method Beyond Nyquist Frequency for Acoustic Wave Equation

An anti-dispersion reverse-time migration method with local nearly-analytic operators and its application

Multidimensional Q‐compensated reverse time migration using a high‐efficient decoupled viscoacoustic wave equation

A New High-Order Scheme Based on Numerical Dispersion Analysis of the Wave Phase Velocity for Semidiscrete Wave Equations

2-D wave equation modeling and migration by a new finite difference scheme based on the Galerkin method

Simulation of Acoustic Wavefields in Heterogeneous Media: A Robust Method for Automatic Suppression of Numerical Dispersion

Wavefield simulation and reverse time migration for acoustic TTI media

Analysis of Numerical Dispersion in Wave-Field Simulation

Dispersive Problem in Seismic Wave Propagation Numerical Modeling

Time-domain Explicit Finite-Difference Method Based on the Mixed-Domain Function Approximation for Acoustic Wave Equation

Efficient and Stable Finite Difference Modelling of Acoustic Wave Propagation in Variable-density Media

A Low-Dispersive Symplectic Partitioned Runge-Kutta Method for Solving Seismic-Wave Equations: I. Scheme and Theoretical Analysis

3D Weak-Dispersion Reverse Time Migration with a Stereomodeling Method

An novel finite difference dispersion error elimination mechanism in the Lax–Wendroff high‐order time discretization

Complex-valued Adaptive-Coefficient Finite-Difference Frequency-Domain Method for Wavefield Modeling Based on the Diffusive-Viscous Wave Equation

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