scholarly journals Large-offset P-wave traveltime in layered transversely isotropic media

Geophysics ◽  
2021 ◽  
pp. 1-49
Author(s):  
Mohammad Mahdi Abedi ◽  
David Pardo

Large-offset seismic data processing, imaging, and velocity estimation require an accurate traveltime approximation over a wide range of offsets. In layered transversely isotropic media with vertical symmetry axis (VTI), the accuracy of traditional traveltime approximations is limited to near offsets. Herein, we propose a new traveltime approximation that maintains the accuracy of the classical equations around zero offset, and exhibits the correct curvilinear asymptote at infinitely large offsets. Our approximation is based on the conventional acoustic assumption. Its equation incorporates six parameters. To define them, we use the Taylor series expansion of the exact traveltime around zero offset, and a new asymptotic series for infinite offset. Our asymptotic equation shows that the traveltime behavior at infinitely large offsets is dominated by the properties of the layer with the maximum horizontal velocity in the sequence. The parameters of our approximation depend on: the effective zero offset traveltime, the normal moveout velocity, the anellipticity, a new large-offset heterogeneity parameter, and the properties of the layer with the maximum horizontal velocity in the sequence. We apply our traveltime approximation: (1) to directly calculate traveltime and ray parameter at given offsets, as analytical forward modeling; and (2) to estimate the first four of the aforementioned parameters for the layers beneath a known high-velocity layer. Our large-offset heterogeneity parameter includes the layering effect on the reflections traveltime.

Geophysics ◽  
1994 ◽  
Vol 59 (4) ◽  
pp. 591-596 ◽  
Author(s):  
Suhas Phadke ◽  
S. Kapotas ◽  
N. Dai ◽  
Ernest R. Kanasewich

Wave propagation in transversely isotropic media is governed by the horizontal and vertical wave velocities. The quasi‐P(qP) wavefront is not an ellipse; therefore, the propagation cannot be described by the wave equation appropriate for elliptically anisotropic media. However, for a limited range of angles from the vertical, the dispersion relation for qP‐waves can be approximated by an ellipse. The horizontal velocity necessary for this approximation is different from the true horizontal velocity and depends upon the physical properties of the media. In the method described here, seismic data is migrated using a 45-degree wave equation for elliptically anisotropic media with the horizontal velocity determined by comparing the 45-degree elliptical dispersion relation and the quasi‐P‐dispersion relation. The method is demonstrated for some synthetic data sets.


Geophysics ◽  
1989 ◽  
Vol 54 (12) ◽  
pp. 1564-1574 ◽  
Author(s):  
B. S. Byun ◽  
D. Corrigan ◽  
J. E. Gaiser

A new velocity analysis technique is presented for analyzing moveout of signals on multichannel surface seismic or VSP data. An approximate, skewed hyperbolic moveout formula is derived for horizontally layered, transversely isotropic media. This formula involves three measurement parameters: the average vertical velocity and horizontal and skew moveout velocities. By extending Dix‐type hyperbolic moveout analysis, we obtain improved coherence over large source‐geophone offsets for more accurate moveout correction. Compared with the stacking velocity obtained by simple hyperbolic analysis methods, the three velocity parameters estimated by this technique contain more physically meaningful geologic information regarding the anisotropy and/or velocity heterogeneity of the subsurface. Synthetic P‐wave model experiments demonstrate that the skewed hyperbolic moveout formula yields an excellent fit to time‐distance curves over a wide range of ray angles. Consequently, the measurement parameters are shown to reflect adequately the characteristics of velocity dependence on ray angle, i.e., velocity anisotropy. The technique is then applied to two field offset VSP data sets to measure and analyze the velocity parameters. The results show that the apparent anisotropy, defined as the ratio between the horizontal moveout velocity and average vertical velocity, correlates reasonably well with lithology. Highly anisotropic shale and chalk exhibit higher horizontal‐to‐vertical velocity ratios and sandstones show lower ratios.


1998 ◽  
Vol 29 (3-4) ◽  
pp. 550-553
Author(s):  
Patrick N. Okoye ◽  
Norman F. Uren ◽  
John A. McDonald

Geophysics ◽  
2002 ◽  
Vol 67 (3) ◽  
pp. 939-951 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin

Reflection moveout of pure modes recorded on conventional‐length spreads is described by a normal‐moveout (NMO) velocity that depends on the orientation of the common‐midpoint (CMP) line. Here, we introduce the concept of NMO‐velocity surfaces, which are obtained by plotting the NMO velocity as the radius‐vector along all possible directions in 3‐D space, and use it to develop Dix‐type averaging and differentiation algorithms in anisotropic heterogeneous media. The intersection of the NMO‐velocity surface with the horizontal plane represents the NMO ellipse that can be estimated from wide‐azimuth reflection data. We demonstrate that the NMO ellipse and conventional‐spread moveout as a whole can be modeled by Dix‐type averaging of specifically oriented cross‐sections of the NMO‐velocity surfaces along the zero‐offset reflection raypath. This formalism is particularly simple to implement for a stack of homogeneous anisotropic layers separated by plane dipping boundaries. Since our method involves computing just a single (zero‐offset) ray for a given reflection event, it can be efficiently used in anisotropic stacking‐velocity tomography. Application of the Dix‐type averaging to layered transversely isotropic media with a vertical symmetry axis (VTI) shows that the presence of dipping interfaces above the reflector makes the P‐wave NMO ellipse dependent on the vertical velocity and anisotropic coefficients ε and δ. In contrast, P‐wave moveout in VTI models with a horizontally layered overburden is fully controlled by the NMO velocity of horizontal events and the Alkhalifah‐Tsvankin coefficient η ≈ ε − δ. Hence, in some laterally heterogeneous, layered VTI models P‐wave reflection data may provide enough information for anisotropic depth processing.


Geophysics ◽  
1996 ◽  
Vol 61 (3) ◽  
pp. 835-845 ◽  
Author(s):  
John Anderson ◽  
Tariq Alkhalifah ◽  
Ilya Tsvankin

The main advantage of Fowler’s dip‐moveout (DMO) method is the ability to perform velocity analysis along with the DMO removal. This feature of Fowler DMO becomes even more attractive in anisotropic media, where imaging methods are hampered by the difficulty in reconstructing the velocity field from surface data. We have devised a Fowler‐type DMO algorithm for transversely isotropic media using the analytic expression for normal‐moveout velocity. The parameter‐estimation procedure is based on the results of Alkhalifah and Tsvankin showing that in transversely isotropic media with a vertical axis of symmetry (VTI) the P‐wave normal‐moveout (NMO) velocity as a function of ray parameter can be described fully by just two coefficients: the zero‐dip NMO velocity [Formula: see text] and the anisotropic parameter η (η reduces to the difference between Thomsen parameters ε and δ in the limit of weak anisotropy). In this extension of Fowler DMO, resampling in the frequency‐wavenumber domain makes it possible to obtain the values of [Formula: see text] and η by inspecting zero‐offset (stacked) panels for different pairs of the two parameters. Since most of the computing time is spent on generating constant‐velocity stacks, the added computational effort caused by the presence of anisotropy is relatively minor. Synthetic and field‐data examples demonstrate that the isotropic Fowler DMO technique fails to generate an accurate zero‐offset section and to obtain the zero‐dip NMO velocity for nonelliptical VTI models. In contrast, this anisotropic algorithm allows one to find the values of the parameters [Formula: see text] and η (sufficient to perform time migration as well) and to correct for the influence of transverse isotropy in the DMO processing. When combined with poststack F-K Stolt migration, this method represents a complete inversion‐processing sequence capable of recovering the effective parameters of transversely isotropic media and producing migrated images for the best‐fit homogeneous anisotropic model.


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.


Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 268-284 ◽  
Author(s):  
Ilya Tsvankin

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.


Sign in / Sign up

Export Citation Format

Share Document