Nonhyperbolic reflection moveout for horizontal transverse isotropy

Geophysics ◽  
1998 ◽  
Vol 63 (5) ◽  
pp. 1738-1753 ◽  
Author(s):  
AbdulFattah Al‐Dajani ◽  
Ilya Tsvankin
Keyword(s):  
Transverse Isotropy ◽  
Single Layer ◽  
Azimuthal Velocity ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Analytic Description ◽  
Pure Mode ◽  
Quartic Term ◽  
Hti Media ◽  
Vertical Transverse Isotropy

The transversely isotropic model with a horizontal axis of symmetry (HTI) has been used extensively in studies of shear‐wave splitting to describe fractured formations with a single system of parallel vertical penny‐shaped cracks. Here, we present an analytic description of longspread reflection moveout in horizontally layered HTI media with arbitrary strength of anisotropy. The hyperbolic moveout equation parameterized by the exact normal‐moveout (NMO) velocity is sufficiently accurate for P-waves on conventional‐length spreads (close to the reflector depth), although the NMO velocity is not, in general, usable for converting time to depth. However, the influence of anisotropy leads to the deviation of the moveout curve from a hyperbola with increasing spread length, even in a single‐layer model. To account for nonhyperbolic moveout, we have derived an exact expression for the azimuthally dependent quartic term of the Taylor series traveltime expansion [t2(x2)] valid for any pure mode in an HTI layer. The quartic moveout coefficient and the NMO velocity are then substituted into the nonhyperbolic moveout equation of Tsvankin and Thomsen, originally designed for vertical transverse isotropy (VTI). Numerical examples for media with both moderate and uncommonly strong nonhyperbolic moveout show that this equation accurately describes azimuthally dependent P-wave reflection traveltimes in an HTI layer, even for spread lengths twice as large as the reflector depth. In multilayered HTI media, the NMO velocity and the quartic moveout coefficient reflect the influence of layering as well as azimuthal anisotropy. We show that the conventional Dix equation for NMO velocity remains entirely valid for any azimuth in HTI media if the group‐velocity vectors (rays) for data in a common‐midpoint (CMP) gather do not deviate from the vertical incidence plane. Although this condition is not exactly satisfied in the presence of azimuthal velocity variations, rms averaging of the interval NMO velocities represents a good approximation for models with moderate azimuthal anisotropy. Furthermore, the quartic moveout coefficient for multilayered HTI media can also be calculated with acceptable accuracy using the known averaging equations for vertical transverse isotropy. This allows us to extend the nonhyperbolic moveout equation to horizontally stratified media composed of any combination of isotropic, VTI, and HTI layers. In addition to providing analytic insight into the behavior of nonhyperbolic moveout, these results can be used in modeling and inversion of reflection traveltimes in azimuthally anisotropic media.

Feasibility of nonhyperbolic moveout inversion in transversely isotropic media

Geophysics ◽  
1998 ◽  
Vol 63 (3) ◽  
pp. 957-969 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Vertical Velocity ◽  
Transverse Isotropy ◽  
Horizontal Velocity ◽  
P Wave ◽  
Vacuum Field ◽  
Percentage Error ◽  
Correlated Errors ◽  
Data Set ◽  
Time Migration ◽  
Quartic Term

Inversion of reflection traveltimes in anisotropic media can provide estimates of anisotropic coefficients required for seismic processing and lithology discrimination. Nonhyperbolic P-wave moveout for transverse isotropy with a vertical symmetry axis (VTI media) is controlled by the parameter η (or, alternatively, by the horizontal velocity Vhor), which is also responsible for the influence of anisotropy on all time‐processing steps, including dip‐moveout (DMO) correction and time migration. Here, we recast the nonhyperbolic moveout equation, originally developed by Tsvankin and Thomsen, as a function of Vhor and normal‐moveout (NMO) velocity Vnmo and introduce a correction factor in the denominator that increases the accuracy at intermediate offsets. Then we apply this equation to obtain Vhor and η from nonhyperbolic semblance analysis on long common midpoint (CMP) spreads and study the accuracy and stability of the inversion procedure. Our error analysis shows that the horizontal velocity becomes relatively well‐constrained by reflection traveltimes if the spreadlength exceeds twice the reflector depth. There is, however, a certain degree of tradeoff between Vhor and Vnmo caused by the interplay between the quadratic and quartic term of the moveout series. Since the errors in Vhor and Vnmo have opposite signs, the absolute error in the parameter η (which depends on the ratio Vhor/Vnmo) turns out to be at least two times bigger than the percentage error in Vhor. Therefore, the inverted value of η is highly sensitive to small correlated errors in reflection traveltimes, with moveout distortions on the order of 3–4 ms leading to errors in η up to ±0.1—even in the simplest model of a single VTI layer. Similar conclusions apply to vertically inhomogeneous media, in which the interval horizontal velocity can be obtained with an accuracy often comparable to that of the NMO velocity, while the interval values of η are distorted by the tradeoff between Vhor and Vnmo that gets amplified by the Dix‐type differentiation procedure. We applied nonhyperbolic semblance analysis to a walkaway VSP data set acquired at Vacuum field, New Mexico, and obtained a significant value of η = 0.19 indicative of nonnegligible anisotropy in this area. Then we combined moveout inversion results with the known vertical velocity to resolve the anisotropic coefficients ε and δ. However, in agreement with our modeling results, η estimation was significantly compounded by the scatter in the measured traveltimes. Certain instability in η inversion has no influence on the results of anisotropic poststack time migration because all kinematically equivalent models obtained from nonhyperbolic moveout give an adequate description of long‐spread reflection traveltimes. Also, inversion of nonhyperbolic moveout provides a relatively accurate horizontal‐velocity function that can be combined with the vertical velocity (if it is available) to estimate the anisotropic coefficient ε. However, η represents a valuable lithology indicator that can be obtained from surface P-wave data. Therefore, for purposes of lithology discrimination, it is preferable to find η by means of the more stable DMO method of Alkhalifah and Tsvankin.

Inversion of reflection traveltimes for transverse isotropy

Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 1095-1107 ◽  
Author(s):  
Ilya Tsvankin ◽  
Leon Thomsen
Keyword(s):  
Wave Velocity ◽  
Vertical Velocity ◽  
Joint Inversion ◽  
Transverse Isotropy ◽  
Taylor Series Expansion ◽  
High Sensitivity ◽  
Transversely Isotropic ◽  
P Wave ◽  
S Wave ◽  
Quartic Term

In anisotropic media, the short‐spread stacking velocity is generally different from the root‐mean‐square vertical velocity. The influence of anisotropy makes it impossible to recover the vertical velocity (or the reflector depth) using hyperbolic moveout analysis on short‐spread, common‐midpoint (CMP) gathers, even if both P‐ and S‐waves are recorded. Hence, we examine the feasibility of inverting long‐spread (nonhyperbolic) reflection moveouts for parameters of transversely isotropic media with a vertical symmetry axis. One possible solution is to recover the quartic term of the Taylor series expansion for [Formula: see text] curves for P‐ and SV‐waves, and to use it to determine the anisotropy. However, this procedure turns out to be unstable because of the ambiguity in the joint inversion of intermediate‐spread (i.e., spreads of about 1.5 times the reflector depth) P and SV moveouts. The nonuniqueness cannot be overcome by using long spreads (twice as large as the reflector depth) if only P‐wave data are included. A general analysis of the P‐wave inverse problem proves the existence of a broad set of models with different vertical velocities, all of which provide a satisfactory fit to the exact traveltimes. This strong ambiguity is explained by a trade‐off between vertical velocity and the parameters of anisotropy on gathers with a limited angle coverage. The accuracy of the inversion procedure may be significantly increased by combining both long‐spread P and SV moveouts. The high sensitivity of the long‐spread SV moveout to the reflector depth permits a less ambiguous inversion. In some cases, the SV moveout alone may be used to recover the vertical S‐wave velocity, and hence the depth. Success of this inversion depends on the spreadlength and degree of SV‐wave velocity anisotropy, as well as on the constraints on the P‐wave vertical velocity.

scholarly journals Prior Analysis in Determination of TI Anisotropy Type through Well-Based Modelling: A Case Study on the Deep Water Environment

2021 ◽  
Vol 873 (1) ◽  
pp. 012102
Author(s):  
Madaniya Oktariena ◽  
Wahyu Triyoso ◽  
Fatkhan Fatkhan ◽  
Sigit Sukmono ◽  
Erlangga Septama ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Deep Water ◽  
Initial Conditions ◽  
Transverse Isotropy ◽  
Water Environment ◽  
Horizontal Stress ◽  
P Wave ◽  
S Wave ◽  
Geomechanical Analysis ◽  
Vertical Transverse Isotropy

Abstract The existence of anisotropy phenomena in the subsurface will affect the image quality of seismic data. Hence a prior knowledge of the type of anisotropy is quite essential, especially when dealing with deep water targets. The preliminary result of the anisotropy of the well-based modelling in deep water exploration and development is discussed in this study. Anisotropy types are modelled for Vertical Transverse Isotropy (VTI) and Horizontal Transverse Isotropy (HTI) based on Thomsen Parameters of ε and γ. The parameters are obtained from DSI Logging paired with reference δ value for modelling. Three initial conditions are then analysed. The first assumption is isotropic, in which the P-Wave Velocity, S-Wave Velocity, and Density Log modelled at their in-situ condition. The second and third assumptions are anisotropy models that are VTI and HTI. In terms of HTI, the result shows that the model of CDP Gather in the offset domain has a weak distortion in Amplitude Variation with Azimuth (AVAz). However, another finding shows a relatively strong hockey effect in far offset, which indicates that the target level is a VTI dominated type. It is supported by the geomechanical analysis result in which vertical stress acts as the maximum principal axis while horizontal stress is close to isotropic one. To sum up, this prior anisotropy knowledge obtained based on this study could guide the efficiency guidance in exploring the deep water environment.

Seismic critical-angle reflectometry: A method to characterize azimuthal anisotropy?

Geophysics ◽  
2007 ◽  
Vol 72 (3) ◽  
pp. D41-D50 ◽  
Author(s):  
Martin Landrø ◽  
Ilya Tsvankin
Keyword(s):  
Critical Angle ◽  
Transverse Isotropy ◽  
Symmetry Plane ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Head Wave ◽  
Orthorhombic Symmetry ◽  
Amplitude Curve ◽  
Diagnostic Features ◽  
Azimuthal Variation

Existing anisotropic parameter-estimation algorithms that operate with long-offset data are based on the inversion of either nonhyperbolic moveout or wide-angle amplitude-variation-with-offset (AVO) response. We show that valuable information about anisotropic reservoirs can also be obtained from the critical angle of reflected waves. To explain the behavior of the critical angle, we develop weak-anisotropy approximations for vertical transverse isotropy and then use Tsvankin’s notation to extend them to azimuthally anisotropic models of orthorhombic symmetry. The P-wave critical-angle reflection in orthorhombic media is particularly sensitive to the parameters [Formula: see text] and [Formula: see text] responsible for the symmetry-plane horizontal velocity in the high-velocity layer. The azimuthal variation of the critical angle for typical orthorhombic models can reach [Formula: see text], which translates into substantial changes in the critical offset of the reflected P-wave. The main diagnostic features of the critical-angle reflection employed in our method include the rapid amplitude increase at the critical angle and the subsequent separation of the head wave. Analysis of exact synthetic seismograms, generated with the reflectivity method, confirms that the azimuthal variation of the critical offset is detectable on wide-azimuth, long-spread data and can be qualitatively described by our linearized equations. Estimation of the critical offset from the amplitude curve of the reflected wave, however, is not straightforward. Additional complications may be caused by the overburden noise train and by the influence of errors in the overburden velocity model on the computation of the critical angle. Still, critical-angle reflectometry should help to constrain the dominant fracture directions and can be combined with other methods to reduce the uncertainty in the estimated anisotropy parameters.

Moveout inversion of P-wave data for horizontal transverse isotropy

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1219-1229 ◽  
Author(s):  
Pedro Contreras ◽  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Vertical Velocity ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
Anisotropic Parameter ◽  
Data Set ◽  
Velocity Models ◽  
Wave Data ◽  
Hti Media

The transversely isotropic model with a horizontal symmetry axis (HTI media) has been extensively used in seismological studies of fractured reservoirs. In this paper, a parameter‐estimation technique originally developed by Grechka and Tsvankin for the more general orthorhombic media is applied to horizontal transverse isotropy. Our methodology is based on the inversion of azimuthally‐dependent P-wave normal‐moveout (NMO) velocities from horizontal and dipping reflectors. If the NMO velocity of a given reflection event is plotted in each azimuthal direction, it forms an ellipse determined by three combinations of medium parameters. The NMO ellipse from a horizontal reflector in HTI media can be inverted for the azimuth β of the symmetry axis, the vertical velocity [Formula: see text], and the Thomsen‐type anisotropic parameter δ(V). We describe a technique for obtaining the remaining (for P-waves) anisotropic parameter η(V) (or ε(V)) from the NMO ellipse corresponding to a dipping reflector of arbitrary azimuth. The interval parameters of vertically inhomogeneous HTI media are recovered using the generalized Dix equation that operates with NMO ellipses for horizontal and dipping events. High accuracy of our method is confirmed by inverting a synthetic multiazimuth P-wave data set generated by ray tracing for a layered HTI medium with depth‐varying orientation of the symmetry axis. Although estimation of η(V) can be carried out by the algorithm developed for orthorhombic media, for more stable results the HTI model has to be used from the outset of the inversion procedure. It should be emphasized that P-wave conventional‐spread moveout data provide enough information to distinguish between HTI and lower‐symmetry models. We show that if the medium has the orthorhombic symmetry and is sufficiently different from HTI, the best‐fit HTI model cannot match the NMO ellipses for both a horizontal and a dipping event. The anisotropic coefficients responsible for P-wave moveout can be combined to estimate the crack density and predict whether the cracks are fluid‐filled or dry. A unique feature of the HTI model that distinguishes it from both vertical transverse isotropy and orthorhombic media is that moveout inversion provides not just zero‐dip NMO velocities and anisotropic coefficients, but also the true vertical velocity. As a result, reflection P-wave data acquired over HTI formations can be used to build velocity models in depth and perform anisotropic depth processing.

Reflection moveout inversion for horizontal transverse isotropy: Accuracy, limitation, and acquisition

Geophysics ◽  
2000 ◽  
Vol 65 (1) ◽  
pp. 222-231 ◽  
Author(s):  
Abdulfattah Al‐Dajani ◽  
Tariq Alkhalifah
Keyword(s):  
Parameter Estimation ◽  
Vertical Velocity ◽  
Transverse Isotropy ◽  
Angular Separation ◽  
P Wave ◽  
Velocity Anisotropy ◽  
Hti Media ◽  
Inversion Procedure ◽  
Anisotropy Strength

Horizontal transverse isotropy (HTI) is the simplest azimuthally anisotropic model used to describe vertical fracturing in an isotropic matrix. Assuming that the subsurface is laterally homogeneous, and using the elliptical variation of P-wave NMO velocity with azimuth measured in at least three different source‐to‐receiver orientations, we can estimate three key parameters of HTI media: the vertical velocity, anisotropy, and the azimuth of the symmetry axis. Such parameter estimation is sensitive to the angular separation between the survey lines in 2-D acquisition or, equivalently, to source‐to‐receiver azimuths in 3-D acquisition and the set of azimuths used in the inversion procedure. The accuracy in estimating the azimuth, in particular, is also sensitive to the strength of anisotropy, while the accuracy in resolving vertical velocity and anisotropy is about the same for any strength of anisotropy. To maximize the accuracy and stability in parameter estimation, it is best to have the azimuths for the source‐to‐receiver directions 60° apart when only three directions are used. This requirement is feasible in land seismic data acquisition where wide azimuthal coverage can be designed. In marine streamer acquisition, however, the azimuthal data coverage is limited. Multiple survey directions are necessary to achieve such wide azimuthal coverage in streamer surveys. To perform the inversion using three azimuth directions, 60° apart, an HTI layer overlain by an azimuthally isotropic overburden should have a time thickness, relative to the total time, of at least the ratio of the error in the NMO (stacking) velocity to the interval anisotropy strength of the HTI layer. Having more than three source‐to‐receiver azimuths (e.g., full azimuthal coverage), however, provides a useful data redundancy that enhances the quality of the estimates, thus allowing acceptable parameter estimation at smaller relative thicknesses.

scholarly journals Reducing Residual Moveout for Long Offset Data in VTI Media Using Padé Approximation

2019 ◽  
Vol 125 ◽  
pp. 15005
Author(s):  
Indah Nur Pratiwi ◽  
Mohammad Syamsu Rosid ◽  
Humbang Purba
Keyword(s):  
Travel Time ◽  
Transverse Isotropy ◽  
Padé Approximation ◽  
Three Dimensional ◽  
Synthetic Data ◽  
Single Layer ◽  
Pade Approximation ◽  
Depth Ratio ◽  
Long Offset ◽  
Vertical Transverse Isotropy

Modification of the hyperbolic travel time equation into non-hyperbolic travel time equation is important to increase the reduction residual moveout for long offset data. Some researchers have modified hyperbolic travel time equation into a non-hyperbolic travel time equation to obtain a more accurate value NMO velocity and parameter an-ellipticity or etha on the large offset to depth ratio (ODR) so that the residual moveout value is smaller mainly in large offset to depth ratio. The aims of research is to increase the reduction value of error residue at long offset data using Padé approximation then compare with several approximations. The method used in this study is to conduct forward modeling of the subsurface coating structure. The results of the three-dimensional analysis show that the Padé approximation has the best accuracy compared to the other travel time equations for ODR value up to 4 with an-ellipticity parameter is varying from 0 to 0.5. Testing of synthetic data for single layer on vertical transverse isotropy (VTI) medium obtained the maximum residual error value produced by the Padé approximation is 0.25% in ODR=4. Therefore, Padé approximation is better than other methods for reducing residual moveout.

scholarly journals Anisotropic P-wave travel-time tomography implementing Thomsen's weak approximation in TOMO3D

Solid Earth ◽  
2019 ◽  
Vol 10 (6) ◽  
pp. 1857-1876
Author(s):  
Adrià Meléndez ◽  
Clara Estela Jiménez ◽  
Valentí Sallarès ◽  
César R. Ranero
Keyword(s):  
Travel Time ◽  
Transverse Isotropy ◽  
P Wave ◽  
Weak Approximation ◽  
Weak Anisotropy ◽  
Simultaneous Inversion ◽  
Advantages And Disadvantages ◽  
Parameter Combination ◽  
Rough Approximation ◽  
Vertical Transverse Isotropy

Abstract. We present the implementation of Thomsen's weak anisotropy approximation for vertical transverse isotropy (VTI) media within TOMO3D, our code for 2-D and 3-D joint refraction and reflection travel-time tomographic inversion. In addition to the inversion of seismic P-wave velocity and reflector depth, the code can now retrieve models of Thomsen's parameters (δ and ε). Here, we test this new implementation following four different strategies on a canonical synthetic experiment in ideal conditions with the purpose of estimating the maximum capabilities and potential weak points of our modeling tool and strategies. First, we study the sensitivity of travel times to the presence of a 25 % anomaly in each of the parameters. Next, we invert for two combinations of parameters (v, δ, ε and v, δ, v⊥), following two inversion strategies, simultaneous and sequential, and compare the results to study their performance and discuss their advantages and disadvantages. Simultaneous inversion is the preferred strategy and the parameter combination (v, δ, ε) produces the best overall results. The only advantage of the parameter combination (v, δ, v⊥) is a better recovery of the magnitude of v. In each case, we derive the fourth parameter from the equation relating ε, v⊥ and v. Recovery of v, ε and v⊥ is satisfactory, whereas δ proves to be impossible to recover even in the most favorable scenario. However, this does not hinder the recovery of the other parameters, and we show that it is still possible to obtain a rough approximation of the δ distribution in the medium by sampling a reasonable range of homogeneous initial δ models and averaging the final δ models that are satisfactory in terms of data fit.

Geometrical spreading of P-waves in horizontally layered, azimuthally anisotropic media

Geophysics ◽  
2005 ◽  
Vol 70 (5) ◽  
pp. D43-D53 ◽  
Author(s):  
Xiaoxia Xu ◽  
Ilya Tsvankin ◽  
Andrés Pech
Keyword(s):  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Azimuthal Velocity ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Geometrical Spreading ◽  
Weak Anisotropy ◽  
P Waves ◽  
Anisotropic Models ◽  
Wide Range

For processing and inverting reflection data, it is convenient to represent geometrical spreading through the reflection traveltime measured at the earth's surface. Such expressions are particularly important for azimuthally anisotropic models in which variations of geometrical spreading with both offset and azimuth can significantly distort the results of wide-azimuth amplitude-variation-with-offset (AVO) analysis. Here, we present an equation for relative geometrical spreading in laterally homogeneous, arbitrarily anisotropic media as a simple function of the spatial derivatives of reflection traveltimes. By employing the Tsvankin-Thomsen nonhyperbolic moveout equation, the spreading is represented through the moveout coefficients, which can be estimated from surface seismic data. This formulation is then applied to P-wave reflections in an orthorhombic layer to evaluate the distortions of the geometrical spreading caused by both polar and azimuthal anisotropy. The relative geometrical spreading of P-waves in homogeneous orthorhombic media is controlled by five parameters that are also responsible for time processing. The weak-anisotropy approximation, verified by numerical tests, shows that azimuthal velocity variations contribute significantly to geometrical spreading, and the existing equations for transversely isotropic media with a vertical symmetry axis (VTI) cannot be applied even in the vertical symmetry planes. The shape of the azimuthally varying spreading factor is close to an ellipse for offsets smaller than the reflector depth but becomes more complicated for larger offset-to-depth ratios. The overall magnitude of the azimuthal variation of the geometrical spreading for the moderately anisotropic model used in the tests exceeds 25% for a wide range of offsets. While the methodology developed here is helpful in modeling and analyzing anisotropic geometrical spreading, its main practical application is in correcting the wide-azimuth AVO signature for the influence of the anisotropic overburden.

Sign in / Sign up

Export Citation Format

Share Document