Processing‐induced anisotropy

Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1920-1928 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Seismic Data ◽  
Symmetry Axis ◽  
Homogeneous Medium ◽  
Transversely Isotropic ◽  
P Wave ◽  
Well Logs ◽  
Velocity Analysis ◽  
Induced Anisotropy ◽  
Vertical Heterogeneity ◽  
Vti Media

Processing of seismic data is often performed under the assumption that the velocity distribution in the subsurface can be approximated by a macromodel composed of isotropic homogeneous layers or blocks. Despite being physically unrealistic, such models are believed to be sufficient for describing the kinematics of reflection arrivals. In this paper, we examine the distortions in normal‐moveout (NMO) velocities caused by the intralayer vertical heterogeneity unaccounted for in velocity analysis. To match P‐wave moveout measurements from a horizontal or a dipping reflector overlaid by a vertically heterogeneous isotropic medium, the effective homogeneous overburden has to be anisotropic. This apparent anisotropy is caused not only by velocity monotonically increasing with depth, but also by random velocity variations similar to those routinely observed in well logs. Assuming that the effective homogeneous medium is transversely isotropic with a vertical symmetry axis (VTI), we express the VTI parameters through the actual depth‐dependent isotropic velocity function. If the reflector is horizontal, combining the NMO and vertical velocities always results in nonnegative values of Thomsen's coefficient δ. For a dipping reflector, the inversion of the P‐wave NMO ellipse yields a nonnegative Alkhalifah‐Tsvankin coefficient η that increases with dip. The values of η obtained by two other methods (2‐D dip‐moveout inversion and nonhyperbolic moveout analysis) are also nonnegative but generally differ from that needed to fit the NMO ellipse. For truly anisotropic (VTI) media, the influence of vertical heterogeneity above the reflector can lead to a bias toward positive δ and η estimates in velocity analysis.

Common-reflection-point migration velocity analysis of 2D P-wave data from TTI media

Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. C65-C79 ◽  
Author(s):  
Ernesto V. Oropeza ◽  
George A. McMechan
Keyword(s):  
Symmetry Axis ◽  
Transversely Isotropic ◽  
Migration Velocity ◽  
P Wave ◽  
Velocity Analysis ◽  
Reflection Point ◽  
Axis Orientation ◽  
Velocity Models ◽  
Migration Velocity Analysis ◽  
Free Data

We have developed a common-reflection-point (CRP)-based kinematic migration velocity analysis for 2D P-wave reflection data to estimate the four transversely isotropic (TI) parameters [Formula: see text], [Formula: see text], and [Formula: see text], and the tilt angle [Formula: see text] of the symmetry axis in a TI medium. In each iteration, the tomographic parameter was updated alternately with prestack anisotropic ray-based migration. Iterations initially used layer stripping to reduce the number of degrees of freedom; after convergence was reached, a couple of more iterations over all parameters and all CRPs ensured global interlayer coupling and parameter interaction. The TI symmetry axis orientation was constrained to be locally perpendicular to the reflectors. The [Formula: see text] dominated the inversion, and so it was weighted less than [Formula: see text] and [Formula: see text] in the parameter updates. Estimates of [Formula: see text] and [Formula: see text] were influenced if the error in [Formula: see text] was [Formula: see text]; estimates of [Formula: see text] were also influenced if the error in [Formula: see text] was [Formula: see text]. Examples included data for a simple model with a homogeneous TI layer whose dips allowed recovery of all anisotropy parameters from noise-free data, and a more realistic model (the BP tilted transversely isotropic (TTI) model) for which only [Formula: see text], [Formula: see text], and [Formula: see text] were recoverable. The adequacy of the traveltimes predicted by the inverted anisotropic models was tested by comparing migrated images and common image gathers, with those produced using the known velocity models.

Velocity analysis for tilted transversely isotropic media: A physical modeling example

Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 904-910 ◽  
Author(s):  
Vladimir Grechka ◽  
Andres Pech ◽  
Ilya Tsvankin ◽  
Baoniu Han
Keyword(s):  
Physical Modeling ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Velocity Analysis ◽  
Isotropic Layer ◽  
Anisotropic Parameter ◽  
Data Set ◽  
Symmetry Direction

Transverse isotropy with a tilted symmetry axis (TTI media) has been recognized as a common feature of shale formations in overthrust areas, such as the Canadian Foothills. Since TTI layers cause serious problems in conventional imaging, it is important to be able to reconstruct the velocity model suitable for anisotropic depth migration. Here, we discuss the results of anisotropic parameter estimation on a physical‐modeling data set. The model represents a simplified version of a typical overthrust section from the Alberta Foothills, with a horizontal reflector overlaid by a bending transversely isotropic layer. Assuming that the TTI layer is homogeneous and the symmetry axis stays perpendicular to its boundaries, we invert P-wave normal‐moveout (NMO) velocities and zero‐offset traveltimes for the symmetry‐direction velocity V0 and the anisotropic parameters ε and δ. The coefficient ε is obtained using the traveltimes of a wave that crosses a dipping TTI block and reflects from the bottom of the model. The inversion for ε is based on analytic expressions for NMO velocity in media with intermediate dipping interfaces. Our estimates of both anisotropic coefficients are close to their actual values. The errors in the inversion, which are associated primarily with the uncertainties in picking the NMO velocities and traveltimes, can be reduced by a straighforward modification of the acquisition geometry. It should be emphasized that the moveout inversion also gives an accurate estimate of the thickness of the TTI layer, thus reconstructing the correct depth scale of the section. Although the physical model used here was relatively simple, our results demonstrate the principal feasibility of anisotropic velocity analysis and imaging in overthrust areas. The main problems in anisotropic processing for TTI models are likely to be caused by the lateral variation of the velocity field and overall structural complexity.

Ray-based gridded tomography for tilted transversely isotropic media

Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. C11-C23 ◽  
Author(s):  
Xiaoxiang Wang ◽  
Ilya Tsvankin
Keyword(s):  
Symmetry Axis ◽  
A Priori ◽  
Transversely Isotropic ◽  
P Wave ◽  
Model Parameters ◽  
Velocity Analysis ◽  
Symmetry Direction ◽  
Reflection Data ◽  
Transversely Isotropic Media ◽  
Isotropic Media

Reflection tomography in the migrated domain can help reconstruct heterogeneous, anisotropic velocity fields needed for accurate depth imaging of complex geologic structures. The presence of anisotropy, however, increases the uncertainty in velocity analysis and typically requires a priori constraints on the model parameters. Here, we develop a 2D P-wave tomographic algorithm for heterogeneous transversely isotropic media with a tilted symmetry axis (TTI) and investigate the conditions necessary for stable estimation of the symmetry-direction velocity [Formula: see text] and the anisotropy parameters [Formula: see text] and [Formula: see text]. The model is divided into rectangular cells, and the parameters [Formula: see text], [Formula: see text], [Formula: see text], and the tilt [Formula: see text] of the symmetry axis are defined at the grid points. To increase the stability of the inversion, the symmetry axis is set orthogonal to the imaged reflectors, with the tilt interpolated inside each layer. The iterative migration velocity analysis involves efficient linearized parameter updating designed to minimize the residual moveout in image gathers for all available reflection events. The moveout equation in the depth-migrated domain includes a nonhyperbolic term that describes long-offset data, which are particularly sensitive to [Formula: see text]. Synthetic tests for models with a “quasi-factorized” TTI syncline (i.e., [Formula: see text] and [Formula: see text] are constant inside the anisotropic layer) and a TTI thrust sheet demonstrate that stable parameter estimation requires either strong smoothness constraints or additional information from walkaway VSP (vertical seismic profiling) traveltimes. If the model is quasi-factorized with a linear spatial variation of [Formula: see text], it may be possible to obtain the interval TTI parameters just from long-spread reflection data.

Velocity analysis using nonhyperbolic moveout in transversely isotropic media

Geophysics ◽  
1997 ◽  
Vol 62 (6) ◽  
pp. 1839-1854 ◽  
Author(s):  
Tariq Alkhalifah
Keyword(s):  
Anisotropy Parameter ◽  
Transversely Isotropic ◽  
P Wave ◽  
Transversely Isotropic Media ◽  
Vertical Heterogeneity ◽  
Isotropic Media ◽  
D Values ◽  
Two Parameters ◽  
Vti Media ◽  
Ti Media

P‐wave reflections from horizontal interfaces in transversely isotropic (TI) media have nonhyperbolic moveout. It has been shown that such moveout as well as all time‐related processing in TI media with a vertical symmetry axis (VTI media) depends on only two parameters, [Formula: see text] and η. These two parameters can be estimated from the dip‐moveout behavior of P‐wave surface seismic data. Alternatively, one could use the nonhyperbolic moveout for parameter estimation. The quality of resulting estimates depends largely on the departure of the moveout from hyperbolic and its sensitivity to the estimated parameters. The size of the nonhyperbolic moveout in TI media is dependent primarily on the anisotropy parameter η. An “effective” version of this parameter provides a useful measure of the nonhyperbolic moveout even in v(z) isotropic media. Moreover, effective η, [Formula: see text], is used to show that the nonhyperbolic moveout associated with typical TI media (e.g., shales, with η ≃ 0.1) is larger than that associated with typical v(z) isotropic media. The departure of the moveout from hyperbolic is increased when typical anisotropy is combined with vertical heterogeneity. Larger offset‐to‐depth ratios (X/D) provide more nonhyperbolic information and, therefore, increased stability and resolution in the inversion for [Formula: see text]. The X/D values (e.g., X/D > 1.5) needed for obtaining stability and resolution are within conventional acquisition limits, especially for shallow targets. Although estimation of η using nonhyperbolic moveouts is not as stable as using the dip‐moveout method of Alkhalifah and Tsvankin, particularly in the absence of large offsets, it does offer some flexibility. It can be applied in the absence of dipping reflectors and also may be used to estimate lateral η variations. Application of the nonhyperbolic inversion to data from offshore Africa demonstrates its usefulness, especially in estimating lateral and vertical variations in η.

Inversion of ultrasonic data for transversely isotropic media

Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. C1-C7 ◽  
Author(s):  
Yevhen Kovalyshen ◽  
Joel Sarout ◽  
Jeremie Dautriat
Keyword(s):  
Numerical Algorithm ◽  
Seismic Data ◽  
Symmetry Axis ◽  
Rock Sample ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Wave Velocities ◽  
Transversely Isotropic Media ◽  
Isotropic Media

We have developed a new numerical algorithm for inversion of ultrasonic data in transversely isotropic media. This algorithm is able to determine from the measured P-wave velocities the orientation of the symmetry axis of a rock sample and the Thomsen’s parameters, only assuming transverse isotropy. The inversion of ultrasonic data acquired on natural and potentially heterogeneous shale samples produced reasonable results. In addition, the algorithm was successfully tested on ultrasonic data acquired on synthetic samples with predefined orientations of the symmetry axis. An additional outcome of the algorithm is a simple approximation of Thomsen’s formulation, which can be effectively used for interpretation of seismic data in transversely isotropic media.

scholarly journals New acoustic approximation for transversely isotropic media with a vertical symmetry axis

Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. C1-C12 ◽  
Author(s):  
Shibo Xu ◽  
Alexey Stovas ◽  
Tariq Alkhalifah ◽  
Hitoshi Mikada
Keyword(s):  
Wave Propagation ◽  
Seismic Data ◽  
Eikonal Equation ◽  
Symmetry Axis ◽  
Transversely Isotropic ◽  
P Wave ◽  
Seismic Data Processing ◽  
S Wave ◽  
P Waves ◽  
Acoustic Approximation

Seismic data processing in the elastic anisotropic model is complicated due to multiparameter dependency. Approximations to the P-wave kinematics are necessary for practical purposes. The acoustic approximation for P-waves in a transversely isotropic medium with a vertical symmetry axis (VTI) simplifies the description of wave propagation in elastic media, and as a result, it is widely adopted in seismic data processing and analysis. However, finite-difference implementations of that approximation are plagued with S-wave artifacts. Specifically, the resulting wavefield also includes artificial diamond-shaped S-waves resulting in a redundant signal for many applications that require pure P-wave data. To derive a totally S-wave-free acoustic approximation, we have developed a new acoustic approximation for pure P-waves that is totally free of S-wave artifacts in the homogeneous VTI model. To keep the S-wave velocity equal to zero, we formulate the vertical S-wave velocity to be a function of the model parameters, rather than setting it to zero. Then, the corresponding P-wave phase and group velocities for the new acoustic approximation are derived. For this new acoustic approximation, the kinematics is described by a new eikonal equation for pure P-wave propagation, which defines the new vertical slowness for the P-waves. The corresponding perturbation-based approximation for our new eikonal equation is used to compare the new equation with the original acoustic eikonal. The accuracy of our new P-wave acoustic approximation is tested on numerical examples for homogeneous and multilayered VTI models. We find that the accuracy of our new acoustic approximation is as good as the original one for the phase velocity, group velocity, and the kinematic parameters such as vertical slowness, traveltime, and relative geometric spreading. Therefore, the S-wave-free acoustic approximation could be further applied in seismic processing that requires pure P-wave data.

scholarly journals Analysis of RTM extended images for VTI media

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. S139-S150 ◽  
Author(s):  
Vladimir Li ◽  
Ilya Tsvankin ◽  
Tariq Alkhalifah
Keyword(s):  
Symmetry Axis ◽  
Transversely Isotropic ◽  
P Wave ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Anisotropy Parameters ◽  
And Migration ◽  
Vti Media

Extended images obtained from reverse time migration (RTM) contain information about the accuracy of the velocity field and subsurface illumination at different incidence angles. Here, we evaluate the influence of errors in the anisotropy parameters on the shape of the residual moveout (RMO) in P-wave RTM extended images for VTI (transversely isotropic with a vertical symmetry axis) media. Using the actual spatial distribution of the zero-dip NMO velocity ([Formula: see text]), which could be approximately estimated by conventional techniques, we analyze the extended images obtained with distorted fields of the parameters [Formula: see text] and [Formula: see text]. Differential semblance optimization (DSO) and stack-power estimates are employed to study the sensitivity of focusing to the anisotropy parameters. We also build angle gathers to facilitate interpretation of the shape of RMO in the extended images. The results show that the signature of [Formula: see text] is dip-dependent, whereas errors in [Formula: see text] cause defocusing only if that parameter is laterally varying. Hence, earlier results regarding the influence of [Formula: see text] and [Formula: see text] on reflection moveout and migration velocity analysis remain generally valid in the extended image space for complex media. The dependence of RMO on errors in the anisotropy parameters provides essential insights for anisotropic wavefield tomography using extended images.

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.

scholarly journals The offset-midpoint traveltime pyramid in 3D transversely isotropic media with a horizontal symmetry axis

Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. T51-T62 ◽  
Author(s):  
Qi Hao ◽  
Alexey Stovas ◽  
Tariq Alkhalifah
Keyword(s):  
Symmetry Axis ◽  
Transversely Isotropic ◽  
P Wave ◽  
Analytic Representation ◽  
Velocity Inversion ◽  
Transversely Isotropic Media ◽  
Horizontal Symmetry ◽  
Time Domains ◽  
Isotropic Media ◽  
Hti Media

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.

Sign in / Sign up

Export Citation Format

Share Document