scholarly journals An acoustic eikonal equation for attenuating transversely isotropic media with a vertical symmetry axis

Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. C9-C20 ◽  
Author(s):  
Qi Hao ◽  
Tariq Alkhalifah
Keyword(s):  
Eikonal Equation ◽  
Wave Attenuation ◽  
Transversely Isotropic ◽  
S Wave ◽  
P Waves ◽  
The Earth ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Vti Media ◽  
Complex Valued

Seismic-wave attenuation is an important component of describing wave propagation. Certain regions, such as gas clouds inside the earth, exert highly localized attenuation. In fact, the anisotropic nature of the earth induces anisotropic attenuation because the quasi P-wave dispersion effect should be profound along the symmetry direction. We have developed a 2D acoustic eikonal equation governing the complex-valued traveltime of quasi P-waves in attenuating, transversely isotropic media with a vertical-symmetry axis (VTI). This equation is derived under the assumption that the complex-valued traveltime of quasi P-waves in attenuating VTI media are independent of the S-wave velocity parameter [Formula: see text] in Thomsen’s notation and the S-wave attenuation coefficient [Formula: see text] in Zhu and Tsvankin’s notation. We combine perturbation theory and Shanks transform to develop practical approximations to the acoustic attenuating eikonal equation, capable of admitting an analytical description of the attenuation in homogeneous media. For a horizontal-attenuating VTI layer, we also derive the nonhyperbolic approximations for the real and imaginary parts of the complex-valued reflection traveltime. These equations reveal that (1) the quasi SV-wave velocity and the corresponding quasi SV-wave attenuation coefficient given as part of Thomsen-type notation barely affect the ray velocity and ray attenuation of quasi P-waves in attenuating VTI media; (2) combining the perturbation method and Shanks transform provides an accurate analytic eikonal solution for homogeneous attenuating VTI media; (3) for a horizontal attenuating VTI layer with weak attenuation, the real part of the complex-valued reflection traveltime may still be described by the existing nonhyperbolic approximations developed for nonattenuating VTI media, and the imaginary part of the complex-valued reflection traveltime still has the shape of nonhyperbolic curves. In addition, we have evaluated the possible extension of the proposed eikonal equation to realistic attenuating media, an alternative perturbation solution to the proposed eikonal equation, and the feasibility of applying the proposed nonhyperbolic equation for the imaginary part of the complex-valued traveltime to invert for interval attenuation parameters.

Three-parameter normal moveout correction in layered anisotropic media: A stretch-free approach

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. C129-C142 ◽  
Author(s):  
Mohammad Mahdi Abedi ◽  
Mohammad Ali Riahi ◽  
Alexey Stovas
Keyword(s):  
Transversely Isotropic ◽  
Real Data ◽  
Data Sets ◽  
P Waves ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Nmo Correction ◽  
Recorded Data ◽  
Normal Moveout Correction ◽  
Vti Media

In conventional normal moveout (NMO) correction, some parts of the recorded data at larger offsets are discarded because of NMO distortions. Deviation from the true traveltime of reflections due to the anisotropy and heterogeneity of the earth, and wavelet stretching are two reasons of these distortions. The magnitudes of both problems increase with increasing the offset to depth ratio. Therefore, to be able to keep larger offsets of shallower reflections, both problems should be obviated. Accordingly, first, we have studied different traveltime approximations being in use, alongside new parameterizations for two classical functional equations, to select suitable equations for NMO correction. We numerically quantify the fitting accuracy and uncertainty of known nonhyperbolic traveltime approximations for P-waves in transversely isotropic media with vertical symmetry axis (VTI). We select three suitable three-parameter approximations for NMO in layered VTI media as the VTI generalized moveout approximation, a double-square-root approximation, and a perturbation-based approximation. Second, we have developed an extension of the earlier proposed stretch-free NMO method, using the selected moveout approximations. This method involves an automatic modification of the input parameters in anisotropic NMO correction, for selected reflections. Our anisotropic stretch-free NMO method is tested on synthetic and three real data sets from Gulf of Mexico and Iranian oil fields. The results verify the success of the method in extending the usable offsets, by generating flat and stretch-free NMO corrected reflections.

Computed waveforms in transversely isotropic media

Geophysics ◽  
1982 ◽  
Vol 47 (5) ◽  
pp. 771-783 ◽  
Author(s):  
J. E. White
Keyword(s):  
Asymptotic Theory ◽  
Transversely Isotropic ◽  
Inversion Method ◽  
S Wave ◽  
Fourier Inversion ◽  
P Waves ◽  
Arrival Times ◽  
S Waves ◽  
Transversely Isotropic Media ◽  
Isotropic Media

Radiation of elastic waves from a point force or from a localized torque into a transversely isotropic medium has been formulated in terms of displacement potentials, and transient waveforms have been computed by numerical Fourier inversion. For isotropic sandstone, this procedure yields P‐ and S‐wave pulses whose arrival times and magnitudes agree with theory. For a range of anisotropic rocks, arrival times of quasi‐P‐waves and quasi‐S‐waves agree with asymptotic theory. For extreme anisotropy, some quasi‐S‐wave pulses arrive at times which are not predicted by asymptotic theory. Magnitudes have not been compared with results of asymptotic theory, but decrease with distance appears to be in agreement. This Fourier inversion method gives near‐source changes in waveform which are not obtainable from the asymptotic theory.

A convenient expression for the NMO velocity function in terms of ray parameter

Geophysics ◽  
1998 ◽  
Vol 63 (1) ◽  
pp. 275-278 ◽  
Author(s):  
Jack K. Cohen
Keyword(s):  
Phase Velocity ◽  
Phase Angle ◽  
Transversely Isotropic ◽  
P Waves ◽  
Velocity Function ◽  
Analysis Methodology ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Vti Media ◽  
Ray Parameter

Tsvankin (1995) derived an expression for the normal‐moveout (NMO) velocity from horizontal or dipping reflectors that is valid in symmetry planes of anisotropic media. This equation is written in terms of the phase velocity and its derivatives with respect to the phase angle with vertical. The phase angle, however, is not a seismic observable, while the ray parameter p is. Indeed, Alkhalifah and Tsvankin (1995) pointed out the importance of representing the NMO velocity function as a function of ray parameter and obtained a number of significant results via an indirect use of this variable. Subsequently, Cohen (1997) directly transformed the Tsvankin expression for the moveout velocity obtaining an expression involving the phase velocity expressed as a function of ray parameter. Applying this result to the normal moveout of P-waves in transversely isotropic media with a vertical symmetry axis (VTI media), he was able to get various approximations for the ray‐parameter form of the moveout function and thus give additional analytic insight into the success of the velocity‐analysis methodology Alkhalifah and Tsvankin.

Acoustic approximations for processing in transversely isotropic media

Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 623-631 ◽  
Author(s):  
Tariq Alkhalifah
Keyword(s):  
Anisotropy Parameter ◽  
Transversely Isotropic ◽  
P Wave ◽  
Elastic Model ◽  
S Wave ◽  
Time Migration ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Two Parameters ◽  
Vti Media

When transversely isotropic (VTI) media with vertical symmetry axes are characterized using the zero‐dip normal moveout (NMO) velocity [[Formula: see text]] and the anisotropy parameter ηinstead of Thomsen’s parameters, time‐related processing [moveout correction, dip moveout (DMO), and time migration] become nearly independent of the vertical P- and S-wave velocities ([Formula: see text] and [Formula: see text], respectively). The independence on [Formula: see text] and [Formula: see text] is well within the limits of seismic accuracy, even for relatively strong anisotropy. The dependency on [Formula: see text] and [Formula: see text] reduces even further as the ratio [Formula: see text] decreases. In fact, for [Formula: see text], all time‐related processing depends exactly on only [Formula: see text] and η. This fortunate dependence on two parameters is demonstrated here through analytical derivations of time‐related processing equations in terms of [Formula: see text] and η. The time‐migration dispersion relation, the NMO velocity for dipping events, and the ray‐tracing equations extracted by setting [Formula: see text] (i.e., by considering VTI as acoustic) not only depend solely on [Formula: see text] and η but are much simpler than the counterpart expressions for elastic media. Errors attributed to this use of the acoustic assumption are small and may be neglected. Therefore, as in isotropic media, the acoustic model arising from setting [Formula: see text], although not exactly true for VTI media, can serve as a useful approximation to the elastic model for the kinematics of P-wave data. This approximation can boost the efficiency of imaging and DMO programs for VTI media as well as simplify their description.

scholarly journals An efficient Helmholtz solver for acoustic transversely isotropic media

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. C75-C83 ◽  
Author(s):  
Zedong Wu ◽  
Tariq Alkhalifah
Keyword(s):  
Wave Equation ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
P Wave ◽  
Seismic Exploration ◽  
S Wave ◽  
P Waves ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Separable Approximation

The acoustic approximation, even for anisotropic media, is widely used in current industry imaging and inversion algorithms mainly because P-waves constitute most of the energy recorded in seismic exploration. The resulting acoustic formulas tend to be simpler, resulting in more efficient implementations, and they depend on fewer medium parameters. However, conventional solutions of the acoustic-wave equation with higher-order derivatives suffer from S-wave artifacts. Thus, we separate the quasi-P-wave propagation in anisotropic media into the elliptic anisotropic operator (free of the artifacts) and the nonelliptic anisotropic components, which form a pseudodifferential operator. We then develop a separable approximation of the dispersion relation of nonelliptic-anisotropic components, specifically for transversely isotropic media. Finally, we iteratively solve the simpler lower-order elliptical wave equation for a modified source function that includes the nonelliptical terms represented in the Fourier domain. A frequency-domain Helmholtz formulation of the approach renders the iterative implementation efficient because the cost is dominated by the lower-upper decomposition of the impedance matrix for the simpler elliptical anisotropic model. In addition, the resulting wavefield is free of S-wave artifacts and has a balanced amplitude. Numerical examples indicate that the method is reasonably accurate and efficient.

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

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. T191-T207
Author(s):  
Xingguo Huang ◽  
Hui Sun ◽  
Zhangqing Sun ◽  
Nuno Vieira da Silva
Keyword(s):  
Point Source ◽  
Eikonal Equation ◽  
Taylor Series Expansion ◽  
Transversely Isotropic ◽  
Complex Point ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Anisotropy Parameters ◽  
Complex Point Source ◽  
Complex Eikonal

The complex traveltime solutions of the complex eikonal equation are the basis of inhomogeneous plane-wave seismic imaging methods, such as Gaussian beam migration and tomography. We have developed analytic approximations for the complex traveltime in transversely isotropic media with a titled symmetry axis, which is defined by a Taylor series expansion over the anisotropy parameters. The formulation for the complex traveltime is developed using perturbation theory and the complex point-source method. The real part of the complex traveltime describes the wavefront, and the imaginary part of the complex traveltime describes the decay of the amplitude of waves away from the central ray. We derive the linearized ordinary differential equations for the coefficients of the Taylor-series expansion using perturbation theory. The analytical solutions for the complex traveltimes are determined by applying the complex point-source method to the background traveltime formula and subsequently obtaining the coefficients from the linearized ordinary differential equations. We investigate the influence of the anisotropy parameters and of the initial width of the ray tube on the accuracy of the computed traveltimes. The analytical formulas, as outlined, are efficient methods for the computation of complex traveltimes from the complex eikonal equation. In addition, those formulas are also effective methods for benchmarking approximated solutions.

A fast-marching eikonal solver for tilted transversely isotropic media

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. S385-S393
Author(s):  
Umair bin Waheed
Keyword(s):  
Transversely Isotropic ◽  
Fast Marching Method ◽  
Fast Marching ◽  
P Waves ◽  
Eikonal Equations ◽  
Transversely Isotropic Media ◽  
Fast Sweeping ◽  
Fast Marching Algorithm ◽  
Isotropic Media ◽  
And Migration

Fast and accurate traveltime computation for quasi-P waves in anisotropic media is an essential ingredient of many seismic processing and interpretation applications such as Kirchhoff modeling and migration, microseismic source localization, and traveltime tomography. Fast-sweeping methods are widely used for solving the anisotropic eikonal equation due to their flexibility in solving general equations compared to the fast-marching method. However, it has been observed that fast sweeping can be much less efficient than fast marching for models with curved characteristics and practical grid sizes. By representing a tilted transversely isotropic (TTI) equation as a sequence of elliptically isotropic (EI) eikonal equations, we determine that the fast-marching algorithm can be used to compute fast and accurate traveltimes for TTI media. The tilt angle is absorbed into the description of the effective EI model; therefore, the adopted approach does not compromise on the solution accuracy. Through tests on benchmark synthetic models, we test our fast-marching algorithm and discover considerable improvement in accuracy by using factorization and a second-order finite-difference stencil. The adopted methodology opens the door to the possibility of using the fast-marching algorithm for a wider class of anisotropic eikonal equations.

Kirchhoff time migration for transversely isotropic media: An application to Trinidad data

Geophysics ◽  
2006 ◽  
Vol 71 (1) ◽  
pp. S29-S35 ◽  
Author(s):  
Tariq Alkhalifah
Keyword(s):  
Symmetry Axis ◽  
Transversely Isotropic ◽  
Data Set ◽  
Time Migration ◽  
Transversely Isotropic Media ◽  
Higher Order Terms ◽  
Isotropic Media ◽  
Migration Method ◽  
Vti Media ◽  
Vertical Inhomogeneity

Using a newly developed nonhyperbolic offset-mid-point traveltime equation for prestack Kirchhoff time migration, instead of the conventional double-square-root (DSR) equation, results in overall better images from anisotropic data. Specifically, prestack Kirchhoff time migration for transversely isotropic media with a vertical symmetry axis (VTI media) is implemented using an analytical offset-midpoint traveltime equation that represents the equivalent of Cheop's pyramid for VTI media. It includes higher-order terms necessary to better handle anisotropy as well as vertical inhomogeneity. Application of this enhanced Kirchhoff time-migration method to the anisotropic Marmousi data set demonstrates the effectiveness of the approach. Further application of the method to field data from Trinidad results in sharper reflectivity images of the subsurface, with the faults better focused and positioned than with images obtained using isotropic methods. The superiority of the anisotropic time migration is evident in the flatness of the image gathers.

Sign in / Sign up

Export Citation Format

Share Document