Traveltime approximations for a layered transversely isotropic medium

Geophysics ◽  
2006 ◽  
Vol 71 (2) ◽  
pp. D23-D33 ◽  
Author(s):  
Bjørn Ursin ◽  
Alexey Stovas
Keyword(s):  
Taylor Series ◽  
Isotropic Medium ◽  
Continued Fraction ◽  
Transversely Isotropic ◽  
Horizontal Distance ◽  
Weak Anisotropy ◽  
Sh Waves ◽  
Numerical Studies ◽  
Anisotropic Layers ◽  
Zero Offset

We consider multiple transmitted, reflected, and converted qP-qSV-waves or multiple transmitted and reflected SH-waves in a horizontally layered medium that is transversely isotropic with a vertical symmetry axis (VTI). Traveltime and offset (horizontal distance) between a source and receiver, not necessarily in the same layer, are expressed as functions of horizontal slowness. These functions are given in terms of a Taylor series in slowness in exactly the same form as for a layered isotropic medium. The coefficients depend on the parameters of the anisotropic layers through which the wave has passed, and there is no weak anisotropy assumption. Using classical formulas, the traveltime or traveltime squared can then be expressed as a Taylor series in even powers of offset. These Taylor series give rise to a shifted hyperbola traveltime approximation and a new continued-fraction approximation, described by four parameters that match the Taylor series up to the sixth power in offset. Further approximations give several simplified continued-fraction approximations, all of which depend on three parameters: zero-offset traveltime, NMO velocity, and a heterogeneity coefficient. The approximations break down when there is a cusp in the group velocity for the qSV-wave. Numerical studies indicate that approximations of traveltime squared are generally better than those for traveltime. A new continued-fraction approximation that depends on three parameters is more accurate than the commonly used continued-fraction approximation and the shifted hyperbola.

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.

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.

Imaging reflections in elliptically anisotropic media

Geophysics ◽  
1988 ◽  
Vol 53 (12) ◽  
pp. 1616-1618 ◽  
Author(s):  
Joe Dellinger ◽  
Francis Muir
Keyword(s):  
Point Source ◽  
Isotropic Medium ◽  
Short Note ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Virtual Image ◽  
Transversely Isotropic Medium ◽  
Sh Waves ◽  
Special Case

In an isotropic medium, waves reflected from a mirror form a virtual image of their source. This property of planar reflectors is generally not true in the presence of anisotropy. In their short note, Blair and Korringa (1987) show that for the special case of SH waves from a point source in a transversely isotropic medium, an aberration‐free image is formed for any orientation of the mirror. While their proof is mathematical, we show the same result in an intuitive, pictorial fashion and in the process discover that although the image is indeed aberration free, it is still distorted.

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.

Upscaling in orthorhombic media: Behavior of elastic parameters in heterogeneous fractured earth

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. C113-C126 ◽  
Author(s):  
Yuriy Ivanov ◽  
Alexey Stovas
Keyword(s):  
Seismic Waves ◽  
Homogeneous Medium ◽  
Transversely Isotropic ◽  
Processing Parameters ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
Weak Anisotropy ◽  
Fluid Saturation ◽  
Fractured Medium ◽  
Anisotropic Layers

A stack of horizontal homogeneous elastic arbitrary anisotropic layers in welded contact in the long-wavelength limit is equivalent to an elastic anisotropic homogeneous medium. Such a medium is characterized by an effective average description adhering to previously derived closed-form formalism. We have used this formalism to study three different inhomogeneous orthorhombic (ORT) models that could represent real geologic scenarios. We have determined that a stack of thin orthorhombic layers with arbitrary azimuths of vertical symmetry planes can be approximated by an effective orthorhombic medium. The most suitable approach for this is to minimize the misfit between the effective anisotropic medium, monoclinic in that case, and the desirable orthorhombic medium. The second model is an interbedding of VTI (transversely isotropic with a vertical symmetry axis) layers with the same layers containing vertical fractures (shales are intrinsically anisotropic and often fractured). We have derived a weak-anisotropy approximation for important P-wave processing parameters as a function of the relative amount of the fractured lithology. To accurately characterize fractures, inversion for the fracture parameters should use a priori information on the relative amount of a fractured medium. However, we have determined that the cracks’ fluid saturation can be estimated without prior knowledge of the relative amount of the fractured layer. We have used field well-log data to demonstrate how fractures can be included in the interval of interest during upscaling. Finally, the third model that we have considered is a useful representation of tilted orthorhombic medium in the case of two-way propagation of seismic waves through it. We have derived a weak anisotropy approximation for traveltime parameters of the reflected P-wave that propagates through a stack of thin beds of tilted orthorhombic symmetry. The tilt of symmetry planes in an orthorhombic medium significantly affects the kinematics of the reflected P-wave and should be properly accounted for to avoid mispositioning of geologic structures in seismic imaging.

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.

Sign in / Sign up

Export Citation Format

Share Document