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.

Hybrid Pseudospectral/Finite‐Difference Numerical Simulation and Stability Analysis of Pure P‐Wave Propagation in Vti Media

A Novel Pseudospectral Time Domain Method for Complex EM Simulations

Description of qP-wave propagation in anisotropic media, Part I: Pseudo-pure-mode wave equations

Numerical Solutions of Quasi-P Wave Equations in VTI Media

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

Alternative Stable Qp Wave Equations in TTI Media with Their Applications for Reverse Time Migration

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

An Approximate Acoustic Wave Equation for VTI Media

Decoupled Wave Equations for P and SV Waves in an Acoustic VTI Media

Quasi-P wave forward modeling in viscoelastic VTI media in frequency-space domain

Wavefield simulation and reverse time migration for acoustic TTI media

A Numerical Study of 1-D Nonlinear P-wave Propagation in Solid

A Two‐dimensional Hybrid Method for Modeling Seismic Wave Propagation in Anisotropic Media

Modeling of wave equations in anisotropic media using finite element method and its stability condition

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

Multidomain Finite-Difference and Chebyshev Pseudospectral Hybrid Method for Elastic Wave Simulation: Two-Dimensional Case

Pseudo-Helmholtz decomposition for elastic VTI wavefield based on wavefront phase direction#xD;#xD;

3D Staggered-Grid High-Order Finite Difference Numerical Simulation of Elastic Wave Equation in Porous Isotropic Media Based on BISQ Model

Helmholtz decomposition with a scalar Poisson equation in elastic anisotropic media

Novel hybrid numerical simulation of the wave equation by combining physical and numerical representation theorems and a review of hybrid methodologies

Simulation and Characteristics Analysis of a Wavefield in a Thermoelastic Medium Adopting the Rotated Staggered-Grid Pseudo-Spectral Method and L-S Theory