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.

A hybrid eikonal solver for accurate first-arrival traveltime computation in anisotropic media with strong contrasts

An Iterative Fast Sweeping Method for the Eikonal Equation in 2D Anisotropic Media on Unstructured Triangular Meshes

An Adaptive Finite-Difference Method for Seismic Traveltime Modeling Based on 3D Eikonal Equation

Calculating Traveltimes in 2D General Tilted Transversely Isotropic Media Using Fast Sweeping Method

On the Solution of the Complex Eikonal Equation in Acoustic VTI Media: a Perturbation Plus Optimization Scheme

Homotopy Solutions of the Acoustic Eikonal Equation for Strongly Attenuating Transversely Isotropic Media with a Vertical Symmetry Axis

Traveltime calculations for qP, qSV, and qSH waves in 3D tilted transversely isotropic media using fast sweeping method with a Newton method local solver

Adjoint-state Traveltime Tomography for Azimuthally Anisotropic Media in Spherical Coordinates

Traveltime Computation and Imaging from Rugged Topography in 3D TTI Media

Wavefield simulation and reverse time migration for acoustic TTI media

A complex point-source solution of the acoustic eikonal equation for Gaussian beams in transversely isotropic media

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

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

An in-silico study of travel time estimation using Eikonal equation with finite element method in ultrasound imaging

2D High-Resolution Crosswell Seismic Traveltime Tomography

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

Kirchhoff migration using eikonal-based computation of traveltime source derivatives

3D Traveltime Computation in Orthorhombic Media Using Layer-by-layer Wavefront Marching

3D Traveltime Computation for Quasi-P-wave in Orthorhombic Media Using Dynamic Programming

First-arrival Traveltime Computation for Quasi-P Waves in 2D Transversely Isotropic Media Using Fermat's Principle-Based Fast Marching