Homotopy solutions of the acoustic eikonal equation for strongly attenuating transversely isotropic media with a vertical symmetry axis

Geophysics ◽  
2021 ◽  
pp. 1-33
Author(s):  
Xingguo Huang ◽  
Stewart Greenhalgh ◽  
Yun Long ◽  
Lin Jun ◽  
Xu Liu
Keyword(s):  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Taylor Series ◽  
Eikonal Equation ◽  
Symmetry Axis ◽  
Taylor Series Expansion ◽  
Transversely Isotropic ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Homotopy Analysis

Many applications in seismology involve the modeling of seismic wave traveltimes in anisotropic media. We present homotopy solutions of the acoustic eikonal equation for P-waves traveltimes in attenuating transversely isotropic media with a vertical symmetry axis. Instead of the commonly used perturbation theory, we use the homotopy analysis method to express the traveltimes by a Taylor series expansion over an embedding parameter. For the derivation, we first perform homotopy analysis of the eikonal equation and derive the linearized ordinary differential equations for the coefficients of the Taylor series expansion. Then, we obtain the homotopy solutions for the traveltimes by solving the linearized ordinary differential equations. Results on approximate formulae investigations demonstrate that the analytical expressions are efficient methods for the computation of traveltimes from the eikonal equation. In addition, these formulas are also effective methods for benchmarking approximated solutions in strongly attenuating anisotropic media.

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.

Transverse isotropy of the upper mantle in the vicinity of Pacific fracture zones

1969 ◽  
Vol 59 (1) ◽  
pp. 59-72
Author(s):  
Robert S. Crosson ◽  
Nikolas I. Christensen
Keyword(s):  
Upper Mantle ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Seismic Wave Propagation ◽  
Fracture Zones ◽  
Mantle Material ◽  
Transversely Isotropic Media ◽  
Horizontal Symmetry ◽  
Isotropic Media

Abstract Several recent investigations suggest that portions of the Earth's upper mantle behave anisotropically to seismic wave propagation. Since several types of anisotropy can produce azimuthal variations in Pn velocities, it is of particular geophysical interest to provide a framework for the recognition of the form or forms of anisotropy most likely to be manifest in the upper mantle. In this paper upper mantle material is assumed to possess the elastic properties of transversely isotropic media. Equations are presented which relate azimuthal variations in Pn velocities to the direction and angle of tilt of the symmetry axis of a transversely isotropic upper mantle. It is shown that the velocity data of Raitt and Shor taken near the Mendocino and Molokai fracture zones can be adequately explained by the assumption of transverse isotropy with a nearly horizontal symmetry axis.

scholarly journals The offset-midpoint traveltime pyramid in 3D transversely isotropic media with a horizontal symmetry axis

Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. T51-T62 ◽  
Author(s):  
Qi Hao ◽  
Alexey Stovas ◽  
Tariq Alkhalifah
Keyword(s):  
Symmetry Axis ◽  
Transversely Isotropic ◽  
P Wave ◽  
Analytic Representation ◽  
Velocity Inversion ◽  
Transversely Isotropic Media ◽  
Horizontal Symmetry ◽  
Time Domains ◽  
Isotropic Media ◽  
Hti Media

Analytic representation of the offset-midpoint traveltime equation for anisotropy is very important for prestack Kirchhoff migration and velocity inversion in anisotropic media. For transversely isotropic media with a vertical symmetry axis, the offset-midpoint traveltime resembles the shape of a Cheops’ pyramid. This is also valid for homogeneous 3D transversely isotropic media with a horizontal symmetry axis (HTI). We extended the offset-midpoint traveltime pyramid to the case of homogeneous 3D HTI. Under the assumption of weak anellipticity of HTI media, we derived an analytic representation of the P-wave traveltime equation and used Shanks transformation to improve the accuracy of horizontal and vertical slownesses. The traveltime pyramid was derived in the depth and time domains. Numerical examples confirmed the accuracy of the proposed approximation for the traveltime function in 3D HTI media.

Normal moveout from dipping reflectors in anisotropic media

Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 268-284 ◽  
Author(s):  
Ilya Tsvankin
Keyword(s):  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Weak Anisotropy ◽  
Anisotropic Models ◽  
Transversely Isotropic Media ◽  
Wide Range ◽  
Isotropic Media ◽  
The Difference

Description of reflection moveout from dipping interfaces is important in developing seismic processing methods for anisotropic media, as well as in the inversion of reflection data. Here, I present a concise analytic expression for normal‐moveout (NMO) velocities valid for a wide range of homogeneous anisotropic models including transverse isotropy with a tilted in‐plane symmetry axis and symmetry planes in orthorhombic media. In transversely isotropic media, NMO velocity for quasi‐P‐waves may deviate substantially from the isotropic cosine‐of‐dip dependence used in conventional constant‐velocity dip‐moveout (DMO) algorithms. However, numerical studies of NMO velocities have revealed no apparent correlation between the conventional measures of anisotropy and errors in the cosine‐of‐dip DMO correction (“DMO errors”). The analytic treatment developed here shows that for transverse isotropy with a vertical symmetry axis, the magnitude of DMO errors is dependent primarily on the difference between Thomsen parameters ε and δ. For the most common case, ε − δ > 0, the cosine‐of‐dip–corrected moveout velocity remains significantly larger than the moveout velocity for a horizontal reflector. DMO errors at a dip of 45 degrees may exceed 20–25 percent, even for weak anisotropy. By comparing analytically derived NMO velocities with moveout velocities calculated on finite spreads, I analyze anisotropy‐induced deviations from hyperbolic moveout for dipping reflectors. For transversely isotropic media with a vertical velocity gradient and typical (positive) values of the difference ε − δ, inhomogeneity tends to reduce (sometimes significantly) the influence of anisotropy on the dip dependence of moveout velocity.

DERIVATION OF 2-POINT ZERO STABLENUMERICAL ALGORITHM OF BLOCK BACKWARD DIFFERENTIATION FORMULA FOR SOLVING FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS

2021 ◽  
Vol 5 (2) ◽  
pp. 579-583
Author(s):  
Muhammad Abdullahi ◽  
Bashir Sule ◽  
Mustapha Isyaku
Keyword(s):  
Differential Equation ◽  
Ordinary Differential Equation ◽  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Taylor Series ◽  
Taylor Series Expansion ◽  
Order Ordinary Differential Equation ◽  
Differentiation Formula ◽  
First Order ◽  
Backward Differentiation Formula

This paper is aimed at deriving a 2-point zero stable numerical algorithm of block backward differentiation formula using Taylor series expansion, for solving first order ordinary differential equation. The order and zero stability of the method are investigated and the derived method is found to be zero stable and of order 3. Hence, the method is suitable for solving first order ordinary differential equation. Implementation of the method has been considered

Reflection and transmission responses in a layered transversely isotropic medium with horizontal symmetry axis

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. C143-C157 ◽  
Author(s):  
Song Jin ◽  
Alexey Stovas
Keyword(s):  
Good Accuracy ◽  
Symmetry Axis ◽  
Local Property ◽  
Transversely Isotropic ◽  
Weak Anisotropy ◽  
Reflection And Transmission ◽  
Numerical Tests ◽  
Transversely Isotropic Media ◽  
Horizontal Symmetry ◽  
Isotropic Media

Seismic wave reflection and transmission (R/T) responses characterize the subsurface local property, and the widely spread anisotropy has considerable influences even at small incident angles. We have considered layered transversely isotropic media with horizontal symmetry axes (HTI), and the symmetry axes were not restricted to be aligned. With the assumption of weak contrast across the interface, linear approximations for R/T coefficients normalized by vertical energy flux are derived based on a simple layered HTI model. We also obtain the approximation with the isotropic background medium under an additional weak anisotropy assumption. Numerical tests illustrate the good accuracy of the approximations compared with the exact results.

Sign in / Sign up

Export Citation Format

Share Document