Inversion of ultrasonic data for transversely isotropic media

Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. C1-C7 ◽  
Author(s):  
Yevhen Kovalyshen ◽  
Joel Sarout ◽  
Jeremie Dautriat
Keyword(s):  
Numerical Algorithm ◽  
Seismic Data ◽  
Symmetry Axis ◽  
Rock Sample ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Wave Velocities ◽  
Transversely Isotropic Media ◽  
Isotropic Media

We have developed a new numerical algorithm for inversion of ultrasonic data in transversely isotropic media. This algorithm is able to determine from the measured P-wave velocities the orientation of the symmetry axis of a rock sample and the Thomsen’s parameters, only assuming transverse isotropy. The inversion of ultrasonic data acquired on natural and potentially heterogeneous shale samples produced reasonable results. In addition, the algorithm was successfully tested on ultrasonic data acquired on synthetic samples with predefined orientations of the symmetry axis. An additional outcome of the algorithm is a simple approximation of Thomsen’s formulation, which can be effectively used for interpretation of seismic data in transversely isotropic media.

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.

Analytic study of the effective parameters for determination of the NMO velocity function in transversely isotropic media

Geophysics ◽  
1997 ◽  
Vol 62 (6) ◽  
pp. 1855-1866 ◽  
Author(s):  
Jack K. Cohen
Keyword(s):  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Velocity Function ◽  
Transversely Isotropic Media ◽  
Wave Phase Velocity ◽  
Isotropic Media ◽  
Anisotropy Parameters ◽  
Vti Media ◽  
Ray Parameter

In their studies of transversely isotropic media with a vertical symmetry axis (VTI media), Alkhalifah and Tsvankin observed that, to a high numerical accuracy, the normal moveout (NMO) velocity for dipping reflectors as a function of ray parameter p depends mainly on just two parameters, each of which can be determined from surface P‐wave observations. They substantiated this result by using the weak‐anisotropy approximation and exploited it to develop a time‐domain processing sequence that takes into account vertical transverse isotropy. In this study, the two‐parameter Alkhalifah‐Tsvankin result was further examined analytically. It was found that although there is (as these authors already observed) some dependence on the remaining parameters of the problem, this dependence is weak, especially in the practically important regimes of weak to moderately strong transverse isotropy and small ray parameter. In each of these regimes, an analytic solution is derived for the anisotropy parameter η required for time‐domain P‐wave imaging in VTI media. In the case of elliptical anisotropy (η = 0), NMO velocity expressed through p is fully controlled just by the zero‐dip NMO velocity—one of the Alkhalifah‐ Tsvankin parameters. The two‐parameter representation of NMO velocity also was shown to be exact in another limit—that of the zero shear‐wave vertical velociy. The analytic results derived here are based on new representations for both the P‐wave phase velocity and normal moveout velocity in terms of the ray parameter, with explicit expressions given for the cases of vanishing onaxis shear speed, weak to moderate transverse isotropy, and small to moderate ray parameter. Using these formulas, I have rederived and, in some cases, extended in a uniform manner various results of Tsvankin, Alkhalifah, and others. Examples include second‐order expansions in the anisotropy parameters for both the P‐wave phase‐velocity function and NMO‐velocity function, as well as expansions in powers of the ray parameter for both of these functions. I have checked these expansions against the corresponding exact functions for several choices of the anisotropy parameters.

The effects of transverse isotropy on reflection amplitude versus offset

Geophysics ◽  
1987 ◽  
Vol 52 (4) ◽  
pp. 564-567 ◽  
Author(s):  
J. Wright
Keyword(s):  
Wave Velocity ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
S Wave ◽  
Reflection Amplitude ◽  
Amplitude Versus Offset ◽  
Transversely Isotropic Media ◽  
P Wave Velocity ◽  
Isotropic Media

Studies have shown that elastic properties of materials such as shale and chalk are anisotropic. With the increasing emphasis on extraction of lithology and fluid content from changes in reflection amplitude with shot‐to‐group offset, one needs to know the effects of anisotropy on reflectivity. Since anisotropy means that velocity depends upon the direction of propagation, this angular dependence of velocity is expected to influence reflectivity changes with offset. These effects might be particularly evident in deltaic sand‐shale sequences since measurements have shown that the P-wave velocity of shales in the horizontal direction can be 20 percent higher than the vertical P-wave velocity. To investigate this behavior, a computer program was written to find the P- and S-wave reflectivities at an interface between two transversely isotropic media with the axis of symmetry perpendicular to the interface. Models for shale‐chalk and shale‐sand P-wave reflectivities were analyzed.

An effective anisotropy parameter in transversely isotropic media

Geophysics ◽  
1987 ◽  
Vol 52 (12) ◽  
pp. 1654-1664 ◽  
Author(s):  
N. C. Banik
Keyword(s):  
Transverse Isotropy ◽  
Anisotropy Parameter ◽  
Transversely Isotropic ◽  
P Wave ◽  
Weak Anisotropy ◽  
Reflection Seismic ◽  
The North ◽  
Reflection Amplitude ◽  
Transversely Isotropic Media ◽  
Isotropic Media

An interesting physical meaning is presented for the anisotropy parameter δ, previously introduced by Thomsen to describe weak anisotropy in transversely isotropic media. Roughly, δ is the difference between the P-wave and SV-wave anisotropies of the medium. The observed systematic depth errors in the North Sea are reexamined in view of the new interpretation of the moveout velocity through δ. The changes in δ at an interface adequately describe the effects of transverse isotropy on the P-wave reflection amplitude, The reflection coefficient expression is linearized in terms of changes in elastic parameters. The linearized expression clearly shows that it is the variation of δ at the interface that gives the anisotropic effects at small incidence angles. Thus, δ effectively describes both the moveout velocity and the reflection amplitude variation, two very important pieces of information in reflection seismic prospecting, in the presence of transverse isotropy.

Approximate dispersion relations for qP-qSV-waves in transversely isotropic media

Geophysics ◽  
2000 ◽  
Vol 65 (3) ◽  
pp. 919-933 ◽  
Author(s):  
Michael A. Schoenberg ◽  
Maarten V. de Hoop
Keyword(s):  
Wave Propagation ◽  
Order Approximation ◽  
Transverse Isotropy ◽  
Slight Modification ◽  
Transversely Isotropic ◽  
P Wave ◽  
Slowness Surface ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
First Order Approximation

To decouple qP and qSV sheets of the slowness surface of a transversely isotropic (TI) medium, a sequence of rational approximations to the solution of the dispersion relation of a TI medium is introduced. Originally conceived to allow isotropic P-wave processing schemes to be generalized to encompass the case of qP-waves in transverse isotropy, the sequence of approximations was found to be applicable to qSV-wave processing as well, although a higher order of approximation is necessary for qSV-waves than for qP-waves to yield the same accuracy. The zeroth‐order approximation, about which all other approximations are taken, is that of elliptical TI, which contains the correct values of slowness and its derivative along and perpendicular to the medium’s axis of symmetry. Successive orders of approximation yield the correct values of successive orders of derivatives in these directions, thereby forcing the approximation into increasingly better fit at the intervening oblique angles. Practically, the first‐order approximation for qP-wave propagation and the second‐order approximation for qSV-wave propagation yield sufficiently accurate results for the typical transverse isotropy found in geological settings. After only slight modification to existing programs, the rational approximation allows for ray tracing, (f-k) domain migration, and split‐step Fourier migration in TI media—with little more difficulty than that encountered presently with such algorithms in isotropic media.

P‐wave signatures and notation for transversely isotropic media: An overview

Geophysics ◽  
1996 ◽  
Vol 61 (2) ◽  
pp. 467-483 ◽  
Author(s):  
Ilya Tsvankin
Keyword(s):  
Wave Propagation ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Seismic Inversion ◽  
P Wave ◽  
Symmetry Direction ◽  
Time Migration ◽  
Transversely Isotropic Media ◽  
Isotropic Media

Progress in seismic inversion and processing in anisotropic media depends on our ability to relate different seismic signatures to the anisotropic parameters. While the conventional notation (stiffness coefficients) is suitable for forward modeling, it is inconvenient in developing analytic insight into the influence of anisotropy on wave propagation. Here, a consistent description of P‐wave signatures in transversely isotropic (TI) media with arbitrary strength of the anisotropy is given in terms of Thomsen notation. The influence of transverse isotropy on P‐wave propagation is shown to be practically independent of the vertical S‐wave velocity [Formula: see text], even in models with strong velocity variations. Therefore, the contribution of transverse isotropy to P‐wave kinematic and dynamic signatures is controlled by just two anisotropic parameters, ε and δ, with the vertical velocity [Formula: see text] being a scaling coefficient in homogeneous models. The distortions of reflection moveouts and amplitudes are not necessarily correlated with the magnitude of velocity anisotropy. The influence of transverse isotropy on P‐wave normal‐moveout (NMO) velocity in a horizontally layered medium, on small‐angle reflection coefficient, and on point‐force radiation in the symmetry direction is entirely determined by the parameter δ. Another group of signatures of interest in reflection seisimology—the dip‐dependence of NMO velocity, magnitude of nonhyperbolic moveout, time‐migration impulse response, and the radiation pattern near vertical—is dependent on both anisotropic parameters (ε and δ) and is primarily governed by the difference between ε and δ. Since P‐wave signatures are so sensitive to the value of ε − δ, application of the elliptical‐anisotropy approximation (ε = δ) in P‐wave processing may lead to significant errors. Many analytic expressions given in the paper remain valid in transversely isotropic media with a tilted symmetry axis. Moreover, the equation for NMO velocity from dipping reflectors, as well as the nonhyperbolic moveout equation, can be used in symmetry planes of any anisotropic media (e.g., orthorhombic).

Ray-based gridded tomography for tilted transversely isotropic media

Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. C11-C23 ◽  
Author(s):  
Xiaoxiang Wang ◽  
Ilya Tsvankin
Keyword(s):  
Symmetry Axis ◽  
A Priori ◽  
Transversely Isotropic ◽  
P Wave ◽  
Model Parameters ◽  
Velocity Analysis ◽  
Symmetry Direction ◽  
Reflection Data ◽  
Transversely Isotropic Media ◽  
Isotropic Media

Reflection tomography in the migrated domain can help reconstruct heterogeneous, anisotropic velocity fields needed for accurate depth imaging of complex geologic structures. The presence of anisotropy, however, increases the uncertainty in velocity analysis and typically requires a priori constraints on the model parameters. Here, we develop a 2D P-wave tomographic algorithm for heterogeneous transversely isotropic media with a tilted symmetry axis (TTI) and investigate the conditions necessary for stable estimation of the symmetry-direction velocity [Formula: see text] and the anisotropy parameters [Formula: see text] and [Formula: see text]. The model is divided into rectangular cells, and the parameters [Formula: see text], [Formula: see text], [Formula: see text], and the tilt [Formula: see text] of the symmetry axis are defined at the grid points. To increase the stability of the inversion, the symmetry axis is set orthogonal to the imaged reflectors, with the tilt interpolated inside each layer. The iterative migration velocity analysis involves efficient linearized parameter updating designed to minimize the residual moveout in image gathers for all available reflection events. The moveout equation in the depth-migrated domain includes a nonhyperbolic term that describes long-offset data, which are particularly sensitive to [Formula: see text]. Synthetic tests for models with a “quasi-factorized” TTI syncline (i.e., [Formula: see text] and [Formula: see text] are constant inside the anisotropic layer) and a TTI thrust sheet demonstrate that stable parameter estimation requires either strong smoothness constraints or additional information from walkaway VSP (vertical seismic profiling) traveltimes. If the model is quasi-factorized with a linear spatial variation of [Formula: see text], it may be possible to obtain the interval TTI parameters just from long-spread reflection data.

Sign in / Sign up

Export Citation Format

Share Document