Velocity analysis using nonhyperbolic moveout in transversely isotropic media

Geophysics ◽  
1997 ◽  
Vol 62 (6) ◽  
pp. 1839-1854 ◽  
Author(s):  
Tariq Alkhalifah

P‐wave reflections from horizontal interfaces in transversely isotropic (TI) media have nonhyperbolic moveout. It has been shown that such moveout as well as all time‐related processing in TI media with a vertical symmetry axis (VTI media) depends on only two parameters, [Formula: see text] and η. These two parameters can be estimated from the dip‐moveout behavior of P‐wave surface seismic data. Alternatively, one could use the nonhyperbolic moveout for parameter estimation. The quality of resulting estimates depends largely on the departure of the moveout from hyperbolic and its sensitivity to the estimated parameters. The size of the nonhyperbolic moveout in TI media is dependent primarily on the anisotropy parameter η. An “effective” version of this parameter provides a useful measure of the nonhyperbolic moveout even in v(z) isotropic media. Moreover, effective η, [Formula: see text], is used to show that the nonhyperbolic moveout associated with typical TI media (e.g., shales, with η ≃ 0.1) is larger than that associated with typical v(z) isotropic media. The departure of the moveout from hyperbolic is increased when typical anisotropy is combined with vertical heterogeneity. Larger offset‐to‐depth ratios (X/D) provide more nonhyperbolic information and, therefore, increased stability and resolution in the inversion for [Formula: see text]. The X/D values (e.g., X/D > 1.5) needed for obtaining stability and resolution are within conventional acquisition limits, especially for shallow targets. Although estimation of η using nonhyperbolic moveouts is not as stable as using the dip‐moveout method of Alkhalifah and Tsvankin, particularly in the absence of large offsets, it does offer some flexibility. It can be applied in the absence of dipping reflectors and also may be used to estimate lateral η variations. Application of the nonhyperbolic inversion to data from offshore Africa demonstrates its usefulness, especially in estimating lateral and vertical variations in η.

Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 623-631 ◽  
Author(s):  
Tariq Alkhalifah

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.


Geophysics ◽  
2000 ◽  
Vol 65 (4) ◽  
pp. 1316-1325 ◽  
Author(s):  
Tariq Alkhalifah

Prestack Kirchhoff time migration for transversely isotropic media with a vertical symmetry axis (VTI media) is implemented using an offset‐midpoint traveltime equation, Cheop’s pyramid equivalent equation for VTI media. The derivation of such an equation for VTI media requires approximations that pertain to high frequency and weak anisotropy. Yet the resultant offset‐midpoint traveltime equation for VTI media is highly accurate for even strong anisotropy. It is also strictly dependent on two parameters: NMO velocity and the anisotropy parameter, η. It reduces to the exact offset‐midpoint traveltime equation for isotropic media when η = 0. In vertically inhomogeneous media, the NMO velocity and η parameters in the offset‐midpoint traveltime equation are replaced by their effective values: the velocity is replaced by the rms velocity and η is given by a more complicated equation that includes summation of the fourth power of velocity.


Geophysics ◽  
1997 ◽  
Vol 62 (6) ◽  
pp. 1855-1866 ◽  
Author(s):  
Jack K. Cohen

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.


Geophysics ◽  
2008 ◽  
Vol 73 (4) ◽  
pp. D53-D62 ◽  
Author(s):  
Huub Douma ◽  
Mirko van der Baan

The [Formula: see text] domain is the natural domain for anisotropy parameter estimation in horizontally layered media. The need to transform the data to the [Formula: see text] domain or to pick traveltimes in the [Formula: see text] domain is, however, a practical disadvantage. To overcome this, we combine [Formula: see text]-derived traveltimes and offsets in horizontally layered transversely isotropic media with a vertical symmetry axis (VTI) with a rational interpolation procedure applied in the [Formula: see text] domain. This combination results in an accurate and efficient [Formula: see text]-based semblance analysis for anisotropy parameter estimation from the moveout of qP-waves in horizontally layered VTI media. The semblance analysis is applied to the moveout to search directly for the interval values of the relevant parameters. To achieve this, the method is applied in a layer-stripping fashion. We demonstrate the method using synthetic data examples and show that it is robust in the presence of random noise and moderate statics.


Geophysics ◽  
1987 ◽  
Vol 52 (12) ◽  
pp. 1654-1664 ◽  
Author(s):  
N. C. Banik

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.


Geophysics ◽  
2014 ◽  
Vol 79 (1) ◽  
pp. C19-C26 ◽  
Author(s):  
Alexey Stovas ◽  
Tariq Alkhalifah

Moveout approximations play a very important role in seismic modeling, inversion, and scanning for parameters in complex media. We developed a scheme to map one-way moveout approximations for transversely isotropic media with a vertical axis of symmetry (VTI), which is widely available, to the tilted case (TTI) by introducing the effective tilt angle. As a result, we obtained highly accurate TTI moveout equations analogous with their VTI counterparts. Our analysis showed that the most accurate approximation is obtained from the mapping of generalized approximation. The new moveout approximations allow for, as the examples demonstrate, accurate description of moveout in the TTI case even for vertical heterogeneity. The proposed moveout approximations can be easily used for inversion in a layered TTI medium because the parameters of these approximations explicitly depend on corresponding effective parameters in a layered VTI medium.


2004 ◽  
Vol 41 (3) ◽  
pp. 349-354 ◽  
Author(s):  
P F Daley ◽  
L R Lines

This paper advocates a modified parameterization for transversely isotropic (TI) media in seismology. Transversely isotropic media may be parameterized by the reference quasi-compressional (qP) and quasi-shear (qSV) velocities, α and β, together with the two parameters, ε and δ. These last two variables account for the deviation of the coupled qP and qSV modes of wave propagation from the isotropic case. The dimensionless quantity ε is a measure of the ellipticity of the qP wavefront. The "strange" parameter δ has been employed as a measure of deviation of the qP wavefront or slowness surface from an ellipsoid of revolution and also of the qSV wavefront or slowness surface from a sphere. As the parameter δ has been described as "conceptually inaccessible"; it is logical to determine an alternative characterization in physically realizable quantities. As in the earlier literature, the defining parameter for the degree of ellipticity of the qP slowness or energy surface in a TI medium is chosen here to be ε. In addition, the deviation of both the qP and qSV surfaces from the degenerate ellipsoidal case will be specified here as σ. An examination of the linearized phase (wavefront normal) velocities, and the linearized PP and SVSV reflection coefficients at an interface between two TI media is undertaken to emphasize the merits of this modified parameterization. All formulae used here have appeared in one form or another in the cited literature. In this paper, a reformulation of existing equations as functions of the two independent variables, (ε, σ), is presented. This results in a more physically realizable context for the linearized formulae.


Geophysics ◽  
2013 ◽  
Vol 78 (5) ◽  
pp. WC103-WC111 ◽  
Author(s):  
Umair bin Waheed ◽  
Tariq Alkhalifah ◽  
Alexey Stovas

Diffractions in seismic data contain valuable information that can help improve our modeling capability for better imaging of the subsurface. They are especially useful for anisotropic media because they inherently possess a wide range of dips necessary to resolve the angular dependence of velocity. We develop a scheme for diffraction traveltime computations based on perturbation of the anellipticity anisotropy parameter for transversely isotropic media with tilted axis of symmetry (TTI). The expansion, therefore, uses an elliptically anisotropic medium with tilt as the background model. This formulation has advantages on two fronts: first, it alleviates the computational complexity associated with solving the TTI eikonal equation, and second, it provides a mechanism to scan for the best-fitting anellipticity parameter [Formula: see text] without the need for repetitive modeling of traveltimes, because the traveltime coefficients of the expansion are independent of the perturbed parameter [Formula: see text]. The accuracy of such an expansion is further enhanced by the use of Shanks transform. We established the effectiveness of the proposed formulation with tests on a homogeneous TTI model and complex media such as the Marmousi and BP models.


Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. T51-T62 ◽  
Author(s):  
Qi Hao ◽  
Alexey Stovas ◽  
Tariq Alkhalifah

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.


Sign in / Sign up

Export Citation Format

Share Document