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.

Anisotropic pseudo-acoustic wave pre-stack reverse time migration and the application of the imaging condition

Isotropic Elastic Reverse Time Migration

Elastic Reverse Time Migration Without Wavefield Decomposition

Reflected wave least squares reverse time migration with angle illumination compensation

An efficient least-squares reverse time migration in image domain

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

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

Anisotropic Wave Separation Elastic Reverse Time Migration Based on the Pseudo-Decoupled Wave Equations in VTI Media

Wavefield simulation and reverse time migration for acoustic TTI media

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

Wave Equation Migration for 3-D VSP Data

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

Target-oriented Q-compensated reverse-time migration by using optimized pure-mode wave equation in anisotropic media

Least-squares RTM Based on Nonstationary Preconditioner for Viscoacoustic Media

Reverse Time Migration Based on Normalized Wavefield Decomposition Imaging Condition

Anisotropic least squares reverse time migration of target region based on Marchenko double-focusing

Acoustic prestack reverse time migration based on the symplectic method

Consistency of the inverse scattering imaging condition, the energy norm imaging condition and the impedance kernel in acoustic and elastic reverse-time migration

Frequency-Domain Reverse-Time Migration With Attenuation Compensation

Attenuation-compensated Reverse Time Migration Based on the Stereo-Modeling Operator