Abstract:Seismic anisotropy has been observed in many exploration areas, especially in sedimentary rocks. The most commonly used anisotropic model is the transverse isotropy (TI) with a vertical (VTI) or a titled (TTI) symmetric axis. For practical application, transverse isotropy is commonly considered under the acoustic approximation. In general, such acoustic approximation can be divided into two categories: coupled acoustic anisotropic wave equations and decoupled pure P-wave equations. One of the problems of the coupled acoustic equations is that they are not really free of shear waves. The artificial SV component is usually considered as numerical artifacts and may cause numerical instabilities in TTI media. However, decoupled pure P-wave equations are completely free from shear wave. In this paper, we firstly evaluate some commonly used phase velocity approximations and conclude that the Harlan's approximation, widely used to derive pure P-wave equations, has low accuracy in a medium with large epsilon. Then we present a method to improve its accuracy by modifying the anelliptic (fractional) term and obtain three modifications. Based on the most accurate one, we derive an improved TTI pure P-wave equation. The resulting wave equation can be formulated either as a second-order equation or as a system of first-order equations. The proposed TTI pure P-wave equation involves pseudo-differential operators which are difficult to be calculated by the finite difference method. We use a hybrid pseudo-spectral/finite-difference (PS/FD) scheme to solve the TTI pure P-wave equation. The pseudo-differential operators are calculated by the PS method and the differential operators are calculated by the high-order rotated-staggered-grid FD method (RSGFD). We test the hybrid PS/RSGFD scheme with various kinds of TTI models. In the first numerical example, we compare the first-order TTI pure P-wave equation with the counterpart elastic equation. We conclude that pure P-wave equation can provide good kinematic and dynamic approximations to the elastic wave equation for homogeneous media with weak-to-moderate anisotropy. The second test model is a heterogeneous wedge TTI model which contains sharp contrasts. This numerical example demonstrates such a hybrid PS/RSGFD scheme can provide a stable P-wave propagator. To further verify the proposed algorithm, we implement an anisotropic reverse-time migration (RTM) example on a small region of BP2007 TTI model. These numerical examples demonstrate the potential of the proposed method.

Modeling of Pure Visco-Qp-wave Propagation in Attenuating Tilted Transversely Isotropic Media Based on Decoupled Fractional Laplacians

Accurate Simulations of Pure-Viscoacoustic Wave Propagation in Tilted Transversely Isotropic Media

Modeling of pure qP- and qSV-waves in tilted transversely isotropic media with the optimal quadratic approximation

A New Highly Accurate and Efficient Pure Visco‐acoustic Wave Equation for Tilted Transversely Isotropic Attenuating Media

Modeling Pure Qp-Wave in Attenuative TI Media Using Complex-Valued Fractional Viscoacoustic Wave Equation

Viscoacoustic Wave Equation for Qp-Wave in Transversely Isotropic Media

Efficient Pure Qp-Wave Modeling and Reverse Time Migration in Tilted Transversely Isotropic Media Calculated by a Finite-Difference Approach

Efficient Pure Qp-Wave Simulation and Reverse Time Migration Imaging for Vertical Transverse Isotropic (VTI) Media

Numerical simulation and reverse time migration using an improved pure P-wave equation in tilted transversely isotropic media

Anisotropic Attenuation Compensated Reverse Time Migration of Pure Qp-Wave in Transversely Isotropic Attenuating Media

Modeling Viscoacoustic Wave Propagation Using a New Spatial Variable-Order Fractional Laplacian Wave Equation

Modeling and Reverse Time Migration of Pure Quasi-P-waves in Complex TI Media with a Low-Rank Decomposition

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

Modified Pure-Viscoacoustic Wave Propagation and Compensated Reverse-Time Migration in Transversely Isotropic Media

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

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

Viscoelastic Wave Propagation Simulation Using New Spatial Variable-Order Fractional Laplacians

Attenuation and Compensation Rtm of Anisotropic Pure Quasi P-Wave Based on Nearly Constant Q Model

Pure quasi-P-wave modeling and imaging using an approximated space-domain pseudo-differential operator in the VTI media

Theory and modelling of constant-Q viscoelastic anisotropic media using fractional derivative

Reverse Time Migration Ofviscoacoustic Anisotropic Pure Qp-Wave Based on Nearly Constant Q Model