Dip moveout in anisotropic media

Geophysics ◽  
1990 ◽  
Vol 55 (7) ◽  
pp. 863-867 ◽  
Author(s):  
N. F. Uren ◽  
G. H. F. Gardner ◽  
J. A. McDonald
Keyword(s):  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Point Scatterer ◽  
Velocity Function ◽  
Processing Step ◽  
Axis Of Symmetry ◽  
Zero Offset ◽  
Seismic Lines ◽  
Seismic Measurements ◽  
Simple Numerical Model

When a common‐midpoint gather is collected above a dipping reflector, the point at which reflection occurs moves updip as the source‐receiver offset increases. Stacking velocity is constant, but it is a function of dip. Hence a stacked trace is not equivalent to a zero‐offset recorded trace. Dip moveout (DMO) is a processing step which converts traces to the equivalent of true zero‐offset records, making migration after stack (MAS) equivalent to migration before stack (MBS). The theory of velocity‐independent Gardner DMO is extended to triaxially elliptically anisotropic media in this paper. It is shown that the transformation is exact for homogeneous elliptically anisotropic media and that, after DMO, the stacking velocity is the horizontal component of the elliptically anisotropic velocity function. By taking three distinct seismic lines, three of the six constants of triaxial elliptical anisotropy may be determined. The remaining three cannot be obtained from surface seismic measurements. A simple numerical model of a point scatterer in a transversely isotropic medium with a tilted axis of symmetry is used to generate examples. The DMO process works when the velocity function is elliptical, but is not exact when the velocity function is nonelliptical.

scholarly journals Transformation to zero offset in transversely isotropic media

Geophysics ◽  
1996 ◽  
Vol 61 (4) ◽  
pp. 947-963 ◽  
Author(s):  
Tariq Alkhalifah
Keyword(s):  
Ray Tracing ◽  
Impulse Response ◽  
Homogeneous Medium ◽  
Anisotropy Parameter ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Amplitude Distribution ◽  
Isotropic Media ◽  
Zero Offset ◽  
Vertical Inhomogeneity

Nearly all dip‐moveout correction (DMO) implementations to date assume isotropic homogeneous media. Usually, this has been acceptable considering the tremendous cost savings of homogeneous isotropic DMO and considering the difficulty of obtaining the anisotropy parameters required for effective implementation. In the presence of typical anisotropy, however, ignoring the anisotropy can yield inadequate results. Since anisotropy may introduce large deviations from hyperbolic moveout, accurate transformation to zero‐offset in anisotropic media should address such nonhyperbolic moveout behavior of reflections. Artley and Hale’s v(z) ray‐tracing‐based DMO, developed for isotropic media, provides an attractive approach to treating such problems. By using a ray‐tracing procedure crafted for anisotropic media, I modify some aspects of their DMO so that it can work for v(z) anisotropic media. DMO impulse responses in typical transversely isotropic (TI) models (such as those associated with shales) deviate substantially from the familiar elliptical shape associated with responses in homogeneous isotropic media (to the extent that triplications arise even where the medium is homogeneous). Such deviations can exceed those caused by vertical inhomogeneity, thus emphasizing the importance of taking anisotropy into account in DMO processing. For isotropic or elliptically anisotropic media, the impulse response is an ellipse; but as the key anisotropy parameter η varies, the shape of the response differs substantially from elliptical. For typical η > 0, the impulse response in TI media tends to broaden compared to the response in an isotropic homogeneous medium, a behavior opposite to that encountered in typical v(z) isotropic media, where the response tends to be squeezed. Furthermore, the amplitude distribution along the DMO operator differs significantly from that for isotropic media. Application of this anisotropic DMO to data from offshore Africa resulted in a considerably better alignment of reflections from horizontal and dipping reflectors in common‐midpoint gather than that obtained using an isotropic DMO. Even the presence of vertical inhomogeneity in this medium could not eliminate the importance of considering the shale‐induced anisotropy.

Theory of traveltime shifts around compacting reservoirs: 3D solutions for heterogeneous anisotropic media

Geophysics ◽  
2009 ◽  
Vol 74 (1) ◽  
pp. D25-D36 ◽  
Author(s):  
Rodrigo Felício Fuck ◽  
Andrey Bakulin ◽  
Ilya Tsvankin
Keyword(s):  
Isotropic Medium ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Time Lapse ◽  
Theory Of Elasticity ◽  
Closed Form Expression ◽  
Reservoir Compaction ◽  
Velocity Changes ◽  
Zero Offset ◽  
Induced Velocity

Time-lapse traveltime shifts of reflection events recorded above hydrocarbon reservoirs can be used to monitor production-related compaction and pore-pressure changes. Existing methodology, however, is limited to zero-offset rays and cannot be applied to traveltime shifts measured on prestack seismic data. We give an analytic 3D description of stress-related traveltime shifts for rays propagating along arbitrary trajectories in heterogeneous anisotropic media. The nonlinear theory of elasticity helps to express the velocity changes in and around the reservoir through the excess stresses associated with reservoir compaction. Because this stress-induced velocity field is both heterogeneous and anisotropic, it should be studied using prestack traveltimes or amplitudes. Then we obtain the traveltime shifts by first-order perturbation of traveltimes that accounts not only for the velocity changes but also for 3D deformation of reflectors. The resulting closed-form expression can be used efficiently for numerical modeling of traveltime shifts and, ultimately, for reconstructing the stress distribution around compacting reservoirs. The analytic results are applied to a 2D model of a compacting rectangular reservoir embedded in an initially homogeneous and isotropic medium. The computed velocity changes around the reservoir are caused primarily by deviatoric stresses and produce a transversely isotropic medium with a variable orientation of the symmetry axis and substantial values of the Thomsen parameters [Formula: see text] and [Formula: see text]. The offset dependence of the traveltime shifts should play a crucial role in estimating the anisotropy parameters and compaction-related deviatoric stress components.

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).

The migrator’s equation for anisotropic media

Geophysics ◽  
1990 ◽  
Vol 55 (11) ◽  
pp. 1429-1434 ◽  
Author(s):  
N. F. Uren ◽  
G. H. F. Gardner ◽  
J. A. McDonald
Keyword(s):  
Physical Model ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Vertical Axis ◽  
Original Structure ◽  
Model Data ◽  
P Waves ◽  
Reflection Seismic ◽  
Isotropic Media ◽  
Zero Offset

The migrator’s equation, which gives the relationship between real and apparent dips on a reflector in zero‐offset reflection seismic sections, may be readily implemented in one step with a frequency‐domain migration algorithm for homogeneous media. Huygens’ principle is used to derive a similar relationship for anisotropic media where velocities are directionally dependent. The anisotropic form of the migrator’s equation is applicable to both elliptically and nonelliptically anisotropic media. Transversely isotropic media are used to demonstrate the performance of an f-k implementation of the migrator’s equation for anisotropic media. In such a medium SH-waves are elliptically anisotropic, while P-waves are nonelliptically anisotropic. Numerical model data and physical model data demonstrate the performance of the algorithm, in each case recovering the original structure. Isotropic and anisotropic migration of anisotropic physical model data are compared experimentally, where the anisotropic velocity function of the medium has a vertical axis of symmetry. Only when anisotropic migration is used is the original structure recovered.

Polarization, phase velocity, and NMO velocity of qP-waves in arbitrary weakly anisotropic media

Geophysics ◽  
1998 ◽  
Vol 63 (5) ◽  
pp. 1754-1766 ◽  
Author(s):  
Ivan Pšenčk ◽  
Dirk Gajewski
Keyword(s):  
Phase Velocity ◽  
Polarization Vector ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Vertical Axis ◽  
Natural Generalization ◽  
Weak Anisotropy ◽  
Isotropic Media ◽  
Axis Of Symmetry ◽  
Approximate Formulas

We present approximate formulas for the qP-wave phase velocity, polarization vector, and normal moveout velocity in an arbitrary weakly anisotropic medium obtained with first‐order perturbation theory. All these quantities are expressed in terms of weak anisotropy (WA) parameters, which represent a natural generalization of parameters introduced by Thomsen. The formulas presented and the WA parameters have properties of Thomsen’s formulas and parameters: (1) the approximate equations are considerably simpler than exact equations for qP-waves, (2) the WA parameters are nondimensional quantities, and (3) in isotropic media, the WA parameters are zero and the corresponding equations reduce to equations for isotropic media. In contrast to Thomsen’s parameters, the WA parameters are related linearly to the density normalized elastic parameters. For the transversely isotropic media with vertical axis of symmetry, the equations presented and the WA parameters reduce to the equations and linearized parameters of Thomsen. The accuracy of the formulas presented is tested on two examples of anisotropic media with relatively strong anisotropy: on a transversely isotropic medium with the horizontal axis of symmetry and on a medium with triclinic anisotropy. Although anisotropy is rather strong, the approximate formulas presented yield satisfactory results.

scholarly journals Theory of interval traveltime parameter estimation in layered anisotropic media

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. C253-C263 ◽  
Author(s):  
Yanadet Sripanich ◽  
Sergey Fomel
Keyword(s):  
Parameter Estimation ◽  
Fourth Order ◽  
Taylor Expansion ◽  
Symmetry Plane ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Interval Parameter ◽  
Effective Parameters ◽  
Zero Offset ◽  
Vti Media

Moveout approximations for reflection traveltimes are typically based on a truncated Taylor expansion of traveltime squared around the zero offset. The fourth-order Taylor expansion involves normal moveout velocities and quartic coefficients. We have derived general expressions for layer-stripping second- and fourth-order parameters in horizontally layered anisotropic strata and specified them for two important cases: horizontally stacked aligned orthorhombic layers and azimuthally rotated orthorhombic layers. In the first of these cases, the formula involving the out-of-symmetry-plane quartic coefficients has a simple functional form and possesses some similarity to the previously known formulas corresponding to the 2D in-symmetry-plane counterparts in vertically transversely isotropic (VTI) media. The error of approximating effective parameters by using approximate VTI formulas can be significant in comparison with the exact formulas that we have derived. We have proposed a framework for deriving Dix-type inversion formulas for interval parameter estimation from traveltime expansion coefficients in the general case and in the specific case of aligned orthorhombic layers. The averaging formulas for calculation of effective parameters and the layer-stripping formulas for interval parameter estimation are readily applicable to 3D seismic reflection processing in layered anisotropic media.

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.

Stacking-velocity tomography for tilted orthorhombic media

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. C171-C180 ◽  
Author(s):  
Qifan Liu ◽  
Ilya Tsvankin
Keyword(s):  
Parameter Estimation ◽  
Subsurface Layer ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
P Waves ◽  
Additional Information ◽  
Reflection Data ◽  
Anisotropic Layers ◽  
Zero Offset ◽  
Velocity Tomography

Tilted orthorhombic (TOR) models are typical for dipping anisotropic layers, such as fractured shales, and can also be due to nonhydrostatic stress fields. Velocity analysis for TOR media, however, is complicated by the large number of independent parameters. Using multicomponent wide-azimuth reflection data, we develop stacking-velocity tomography to estimate the interval parameters of TOR media composed of homogeneous layers separated by plane dipping interfaces. The normal-moveout (NMO) ellipses, zero-offset traveltimes, and reflection time slopes of P-waves and split S-waves ([Formula: see text] and [Formula: see text]) are used to invert for the interval TOR parameters including the orientation of the symmetry planes. We show that the inversion can be facilitated by assuming that the reflector coincides with one of the symmetry planes, which is a common geologic constraint often employed for tilted transversely isotropic media. This constraint makes the inversion for a single TOR layer feasible even when the initial model is purely isotropic. If the dip plane is also aligned with one of the symmetry planes, we show that the inverse problem for [Formula: see text]-, [Formula: see text]-, and [Formula: see text]-waves can be solved analytically. When only [Formula: see text]-wave data are available, parameter estimation requires combining NMO ellipses from a horizontal and dipping interface. Because of the increase in the number of independent measurements for layered TOR media, constraining the reflector orientation is required only for the subsurface layer. However, the inversion results generally deteriorate with depth because of error accumulation. Using tests on synthetic data, we demonstrate that additional information such as knowledge of the vertical velocities (which may be available from check shots or well logs) and the constraint on the reflector orientation can significantly improve the accuracy and stability of interval parameter estimation.

Sign in / Sign up

Export Citation Format

Share Document