Nonhyperbolic reflection moveout in anisotropic media

Geophysics ◽  
1994 ◽  
Vol 59 (8) ◽  
pp. 1290-1304 ◽  
Author(s):  
Ilya Tsvankin ◽  
Leon Thomsen
Keyword(s):  
Fourth Order ◽  
Taylor Series ◽  
Transverse Isotropy ◽  
Taylor Series Expansion ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Analytical Approximation ◽  
Qualitative Description ◽  
Weak Anisotropy ◽  
Anisotropic Models

The standard hyperbolic approximation for reflection moveouts in layered media is accurate only for relatively short spreads, even if the layers are isotropic. Velocity anisotropy may significantly enhance deviations from hyperbolic moveout. Nonhyperbolic analysis in anisotropic media is also important because conventional hyperbolic moveout processing on short spreads is insufficient to recover the true vertical velocity (hence the depth). We present analytic and numerical analysis of the combined influence of vertical transverse isotropy and layering on long‐spread reflection moveouts. Qualitative description of nonhyperbolic moveout on “intermediate” spreads (offset‐to‐depth ratio x/z  < 1.7–2) is given in terms of the exact fourth‐order Taylor series expansion for P, SV, and P‐SV traveltime curves, valid for multilayered transversely isotropic media with arbitrary strength of anisotropy. We use this expansion to provide an analytic explanation for deviations from hyperbolic moveout, such as the strongly nonhyperbolic SV‐moveout observed numerically in the case where δ < ε. With this expansion, we also show that the weak anisotropy approximation becomes inadequate (to describe nonhyperbolic moveout) for surprisingly small values of the anisotropies δ and ε. However, the fourth‐order Taylor series rapidly loses numerical accuracy with increasing offset. We suggest a new, more general analytical approximation, and test it against several transversely isotropic models. For P‐waves, this moveout equation remains numerically accurate even for substantial anisotropy and large offsets. This approximation provides a fast and effective way to estimate the behavior of long‐spread moveouts for layered anisotropic models.

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.

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.

scholarly journals Dip‐moveout processing by Fourier transform in anisotropic media

Geophysics ◽  
1997 ◽  
Vol 62 (4) ◽  
pp. 1260-1269 ◽  
Author(s):  
John E. Anderson ◽  
Ilya Tsvankin
Keyword(s):  
Transverse Isotropy ◽  
Computing Time ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
P Wave ◽  
Surface Reflection ◽  
Anisotropic Models ◽  
Reflection Data ◽  
Isotropic Media ◽  
Vti Media

Conventional dip‐moveout (DMO) processing is designed for isotropic media and cannot handle angle‐dependent velocity. We show that Hale's isotropic DMO algorithm remains valid for elliptical anisotropy but may lead to serious errors for nonelliptical models, even if velocity anisotropy is moderate. Here, Hale's constant‐velocity DMO method is extended to anisotropic media. The DMO operator, to be applied to common‐offset data corrected for normal moveout (NMO), is based on the analytic expression for dip‐dependent NMO velocity given by Tsvankin. Since DMO correction in anisotropic media requires knowledge of the velocity field, it should be preceded by an inversion procedure designed to obtain the normal‐moveout velocity as a function of ray parameter. For transversely isotropic models with a vertical symmetry axis (VTI media), P‐wave NMO velocity depends on a single anisotropic coefficient (η) that can be determined from surface reflection data. Impulse responses and synthetic examples for typical VTI media demonstrate the accuracy and efficiency of this DMO technique. Once the inversion step has been completed, the NMO-DMO sequence does not take any more computing time than the genetic Hale method in isotropic media. Our DMO operator is not limited to vertical transverse isotropy as it can be applied in the same fashion in symmetry planes of more complicated anisotropic models such as orthorhombic.

Geometrical spreading in a layered transversely isotropic medium with vertical symmetry axis

Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 2082-2091 ◽  
Author(s):  
Bjørn Ursin ◽  
Ketil Hokstad
Keyword(s):  
Ray Tracing ◽  
Seismic Data ◽  
Symmetry Axis ◽  
Transversely Isotropic ◽  
Analytical Approximation ◽  
Geometrical Spreading ◽  
S Wave ◽  
Weak Anisotropy ◽  
P Waves ◽  
Amplitude Versus

Compensation for geometrical spreading is important in prestack Kirchhoff migration and in amplitude versus offset/amplitude versus angle (AVO/AVA) analysis of seismic data. We present equations for the relative geometrical spreading of reflected and transmitted P‐ and S‐wave in horizontally layered transversely isotropic media with vertical symmetry axis (VTI). We show that relatively simple expressions are obtained when the geometrical spreading is expressed in terms of group velocities. In weakly anisotropic media, we obtain simple expressions also in terms of phase velocities. Also, we derive analytical equations for geometrical spreading based on the nonhyperbolic traveltime formula of Tsvankin and Thomsen, such that the geometrical spreading can be expressed in terms of the parameters used in time processing of seismic data. Comparison with numerical ray tracing demonstrates that the weak anisotropy approximation to geometrical spreading is accurate for P‐waves. It is less accurate for SV‐waves, but has qualitatively the correct form. For P waves, the nonhyperbolic equation for geometrical spreading compares favorably with ray‐tracing results for offset‐depth ratios less than five. For SV‐waves, the analytical approximation is accurate only at small offsets, and breaks down at offset‐depth ratios less than unity. The numerical results are in agreement with the range of validity for the nonhyperbolic traveltime equations.

Reflection moveout approximations for P-waves in a moderately anisotropic homogeneous tilted transverse isotropy layer

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. C175-C185 ◽  
Author(s):  
Ivan Pšenčík ◽  
Véronique Farra
Keyword(s):  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Weak Anisotropy ◽  
P Waves ◽  
3D Space ◽  
Axis Of Symmetry ◽  
Relative Errors ◽  
Symmetry Tests ◽  
Ti Media

We have developed approximate nonhyperbolic P-wave moveout formulas applicable to weakly or moderately anisotropic media of arbitrary anisotropy symmetry and orientation. Instead of the commonly used Taylor expansion of the square of the reflection traveltime in terms of the square of the offset, we expand the square of the reflection traveltime in terms of weak-anisotropy (WA) parameters. No acoustic approximation is used. We specify the formulas designed for anisotropy of arbitrary symmetry for the transversely isotropic (TI) media with the axis of symmetry oriented arbitrarily in the 3D space. Resulting formulas depend on three P-wave WA parameters specifying the TI symmetry and two angles specifying the orientation of the axis of symmetry. Tests of the accuracy of the more accurate of the approximate formulas indicate that maximum relative errors do not exceed 0.3% or 2.5% for weak or moderate P-wave anisotropy, respectively.

The effect of transverse isotropy on isotropic traveltime inversion of vertical seismic profile data—A modeling study

Geophysics ◽  
1990 ◽  
Vol 55 (9) ◽  
pp. 1235-1241 ◽  
Author(s):  
Jan Douma
Keyword(s):  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Vertical Axis ◽  
Seismic Profile ◽  
Isotropic Layer ◽  
Low Velocity ◽  
Traveltime Inversion ◽  
Axis Of Symmetry ◽  
Transversely Isotropic Layer

Traveltime inversion of multioffset VSP data can be used to determine the depths of the interfaces in layered media. Many inversion schemes, however, assume isotropy and consequently may introduce erroneous structures for anisotropic media. Synthetic traveltime data are computed for layered anisotropic media and inverted assuming isotropic layers. Only the interfaces between these layers are inverted. For a medium consisting of a horizontal isotropic low‐velocity layer on top of a transversely isotropic layer with a horizontal axis of symmetry (e.g., anisotropy due to aligned vertical cracks), 2-D isotropic inversion results in an anticline. For a given axis of symmetry the form of this anticline depends on the azimuth of the source‐borehole direction. The inversion result is a syncline (in 3-D a “bowl” structure), regardless of the azimuth of the source‐borehole direction for a vertical axis of symmetry of the transversely isotropic layer (e.g., anisotropy due to horizontal bedding).

Analytic study of the geometrical spreading of P-waves in a layered transversely isotropic medium with a vertical symmetry axis

Geophysics ◽  
2000 ◽  
Vol 65 (4) ◽  
pp. 1305-1315 ◽  
Author(s):  
Hongbo Zhou ◽  
George A. McMechan
Keyword(s):  
Isotropic Medium ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Analytical Formula ◽  
Transversely Isotropic ◽  
P Wave ◽  
Geometrical Spreading ◽  
Transversely Isotropic Medium ◽  
Weak Anisotropy ◽  
Special Cases

An analytical formula for geometrical spreading is derived for a horizontally layered transversely isotropic medium with a vertical symmetry axis (VTI). With this expression, geometrical spreading can be determined using only the anisotropy parameters in the first layer, the traveltime derivatives, and the source‐receiver offset. Explicit, numerically feasible expressions for geometrical spreading are obtained for special cases of transverse isotropy (weak anisotropy and elliptic anisotropy). Geometrical spreading can be calculated for transversly isotropic (TI) media by using picked traveltimes of primary nonhyperbolic P-wave reflections without having to know the actual parameters in the deeper subsurface; no ray tracing is needed. Synthetic examples verify the algorithm and show that it is numerically feasible for calculation of geometrical spreading. For media with a few (4–5) layers, relative errors in the computed geometrical spreading remain less than 0.5% for offset/depth ratios less than 1.0. Errors that change with offset are attributed to inaccuracy in the expression used for nonhyberbolic moveout. Geometrical spreading is most sensitive to errors in NMO velocity, followed by errors in zero‐offset reflection time, followed by errors in anisotropy of the surface layer. New relations between group and phase velocities and between group and phase angles are shown in appendices.

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.

Building tilted transversely isotropic depth models using localized anisotropic tomography with well information

Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. D27-D36 ◽  
Author(s):  
Andrey Bakulin ◽  
Marta Woodward ◽  
Dave Nichols ◽  
Konstantin Osypov ◽  
Olga Zdraveva
Keyword(s):  
Seismic Data ◽  
Model Building ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Synthetic Data ◽  
Rock Physics ◽  
Transversely Isotropic ◽  
Data Set ◽  
Anisotropic Models ◽  
Well Data

Tilted transverse isotropy (TTI) is increasingly recognized as a more geologically plausible description of anisotropy in sedimentary formations than vertical transverse isotropy (VTI). Although model-building approaches for VTI media are well understood, similar approaches for TTI media are in their infancy, even when the symmetry-axis direction is assumed known. We describe a tomographic approach that builds localized anisotropic models by jointly inverting surface-seismic and well data. We present a synthetic data example of anisotropic tomography applied to a layered TTI model with a symmetry-axis tilt of 45 degrees. We demonstrate three scenarios for constraining the solution. In the first scenario, velocity along the symmetry axis is known and tomography inverts for Thomsen’s [Formula: see text] and [Formula: see text] parame-ters. In the second scenario, tomography inverts for [Formula: see text], [Formula: see text], and velocity, using surface-seismic data and vertical check-shot traveltimes. In contrast to the VTI case, both these inversions are nonunique. To combat nonuniqueness, in the third scenario, we supplement check-shot and seismic data with the [Formula: see text] profile from an offset well. This allows recovery of the correct profiles for velocity along the symmetry axis and [Formula: see text]. We conclude that TTI is more ambiguous than VTI for model building. Additional well data or rock-physics assumptions may be required to constrain the tomography and arrive at geologically plausible TTI models. Furthermore, we demonstrate that VTI models with atypical Thomsen parameters can also fit the same joint seismic and check-shot data set. In this case, although imaging with VTI models can focus the TTI data and match vertical event depths, it leads to substantial lateral mispositioning of the reflections.

Fluid dependence of anisotropy parameters in weakly anisotropic porous media

Geophysics ◽  
2013 ◽  
Vol 78 (5) ◽  
pp. WC137-WC145 ◽  
Author(s):  
Olivia Collet ◽  
Boris Gurevich
Keyword(s):  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Threshold Value ◽  
Seismic Velocities ◽  
Anisotropic Porous Media ◽  
Weak Anisotropy ◽  
Shear Moduli ◽  
Saturated Media ◽  
Anisotropy Parameters ◽  
Saturated Medium

Predicting seismic velocities in isotropic fluid-saturated rocks is commonly done using the isotropic Gassmann theory. For anisotropic media, the solution is expressed in terms of stiffness or compliance, which does not provide an intuitive understanding on how the fluid affects wave propagation in anisotropic media. Assuming weak anisotropy, we expressed the anisotropy parameters of transversely isotropic saturated media as a function of the anisotropy parameters in the dry medium, the bulk and shear moduli of the saturated and dry media, the grain and fluid bulk moduli, and the porosity. By deriving an approximation of the anellipticity parameter [Formula: see text], we discovered that if the dry medium was elliptical, the saturated medium was also elliptical but only if the porosity exceeded a certain threshold value. This result can provide a way of differentiating between stress- and fracture-induced anisotropy.

Sign in / Sign up

Export Citation Format

Share Document