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.

Prestack Reverse-Time Migration of Acoustic Wave Equations in VTI Media

Isotropic Elastic Reverse Time Migration

Elastic Reverse Time Migration Without Wavefield Decomposition

Wave Equation Migration for 3-D VSP Data

Prestack Depth Migration Method Based on Quasi-Linear Born Approximation

Generalized Screen Method for Prestack Time Migration

A Method of Plane Wave Depth Migration for Pre-Stack Seismic Data

2-D PP and PS Imaging for VTI with Weak Equation-Converted S-waves

Reverse-time Migration in Acoustic VTI Media Using a High-Order Stereo Operator

Steeply-dipping structures migration based on the wavefield vertical-and horizontal-extrapolation

An efficient least-squares reverse time migration in image domain

Reflected wave least squares reverse time migration with angle illumination compensation

Wavefield simulation and reverse time migration for acoustic TTI media

Viscoacoustic VTI media reverse time migration method in the complex-domain based on the staining algorithm

Reverse-Time Migration for Pure qP-Wave Based on Elliptical Decomposition Vector Equation

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

An Efficient Symplectic Reverse Time Migration Method Using A Local Nearly Analytic Discrete Operator in Acoustic Transversely Isotropic Media with A Vertical Symmetry Axis

Visco‐acoustic Offset Plane‐wave Finite‐difference Prestack Time Migration

Viscoelastic true-amplitude prestack reverse-time depth migration

Frequency-Domain Reverse-Time Migration With Attenuation Compensation

The Optimization and Application of The Propagated Riemannian Wavefield Extrapolator in VTI Media- Pseudo-Depth Domain Least-Squares Reverse-Time Migration