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.

Frequency-domain weighted-averaging finite-difference numerical simulation of qP wave propagation in TTI media

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

Efficient least-squares reverse time migration in TTI media using a finite-difference solvable pure qP-wave equation

Stable Q-compensated reverse time migration in TTI media based on a modified fractional Laplacian pure-viscoacoustic wave equation

Power-Law Frequency-Dependent Q Simulation and Reverse-Time Migration in Transversely Isotropic Acoustic Media

Broadband finite-difference <i>Q</i>-compensated reverse time migration algorithm for tilted transverse isotropic media

An improved pure quasi-P-wave equation for complex anisotropic media

Frequency-domain elastic-wave modeling for VTI media with an average-derivative optimal method

Wavefield simulation of acoustic VTI wave equation based on adaptive-coefficient finite-difference frequency-domain method

Wavefield simulation and reverse time migration for acoustic TTI media

A Nonbalanced Staggered-Grid FDTD Scheme for the First-Order Elastic-Wave Extrapolation and Reverse-Time Migration

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

Numerical simulation of elastic wave based on the staggered grid finite difference method

A time-space domain stereo finite difference method for 3D scalar wave propagation

Frequency-dependent Q simulation and viscoacoustic reverse time migration based on the fractional Zener model

Power-law frequency-dependent Q simulations in viscoacoustic media using decoupled fractional Laplacians

A stable Q-compensated reverse-time migration method using a modified fractional viscoacoustic wave equation

Time-Fractional Viscoacoustic Wave Equation-Based Frequency-Domain Stable Q-RTM

Average‐derivative optimized 21‐point and improved 25‐point forward modelling and full waveform inversion in frequency domain

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

An Optimized Time-Space-domain Finite Difference Method with Piecewise Constant Interpolation Coefficients for Scalar Wave Propagation