Velocity analysis for tilted TI media: A physical‐modeling example

1999 ◽  
Author(s):  
Vladimir Grechka ◽  
Andres Pech ◽  
Ilya Tsvankin ◽  
Baoniu Han
Geophysics ◽  
2006 ◽  
Vol 71 (4) ◽  
pp. D135-D143 ◽  
Author(s):  
Pawan Dewangan ◽  
Ilya Tsvankin ◽  
Mike Batzle ◽  
Kasper van Wijk ◽  
Matthew Haney

Mode-converted PS-waves can provide critically important information for velocity analysis in transversely isotropic (TI) media. We demonstrate, with physical-modeling data, that the combination of long-spread reflection traveltimes of PP- and PS-waves can be inverted for the parameters of a horizontal TI layer with a tilted symmetry axis. The 2D multicomponent reflection data are acquired over a phenolic sample manufactured to simulate the effective medium formed by steeply dipping fracture sets or shale layers. The reflection moveout of PS-waves in this model is asymmetric with respect to the source and receiver positions, and the moveout-asymmetry attributes play a crucial role in constraining the TI parameters. Applying the modified [Formula: see text] method to the PP and PS traveltimes recorded in the symmetry-axis plane, we compute the time and offset asymmetry attributes of the PS-waves along with the traveltimes of the pure SS reflections. The algorithm of Dewangan and Tsvankin is then used to invert the combination of the moveout attributes of PP-, SS-, and PS-waves for the medium parameters and the thickness of the sample. It should be emphasized that the pure-mode (PP and SS) traveltimes alone are insufficient for the inversion, even if 3D wide-azimuth data are available. Our estimates of the symmetry axis tilt and layer thickness almost coincide with the actual values. The inverted model is also validated by reproducing the results of transmission experiments with both P- and S-wave sources. The transmitted SV wavefield exhibits a prominent cusp (triplication) accurately predicted by the parameter-estimation results.


Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. U87-U96 ◽  
Author(s):  
Mamoru Takanashi ◽  
Ilya Tsvankin

One of the most serious problems in anisotropic velocity analysis is the trade-off between anisotropy and lateral heterogeneity, especially if velocity varies on a scale smaller than the maximum offset. We have developed a P-wave MVA (migration velocity analysis) algorithm for transversely isotropic (TI) models that include layers with small-scale lateral heterogeneity. Each layer is described by constant Thomsen parameters [Formula: see text] and [Formula: see text] and the symmetry-direction velocity [Formula: see text] that varies as a quadratic function of the distance along the layer boundaries. For tilted TI media (TTI), the symmetry axis is taken orthogonal to the reflectors. We analyzed the influence of lateral heterogeneity on image gathers obtained after prestack depth migration and found that quadratic lateral velocity variation in the overburden can significantly distort the moveout of the target reflection. Consequently, medium parameters beneath the heterogeneous layer(s) are estimated with substantial error, even when borehole information (e.g., check shots or sonic logs) is available. Because residual moveout in the image gathers is highly sensitive to lateral heterogeneity in the overburden, our algorithm simultaneously inverts for the interval parameters of all layers. Synthetic tests for models with a gently dipping overburden demonstrate that if the vertical profile of the symmetry-direction velocity [Formula: see text] is known at one location, the algorithm can reconstruct the other relevant parameters of TI models. The proposed approach helps increase the robustness of anisotropic velocity model-building and enhance image quality in the presence of small-scale lateral heterogeneity in the overburden.


2005 ◽  
Author(s):  
Pawan Dewangan ◽  
Ilya Tsvankin ◽  
Mike Batzle ◽  
Kasper van Wijk ◽  
Matt Haney

Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1550-1566 ◽  
Author(s):  
Tariq Alkhalifah ◽  
Ilya Tsvankin

The main difficulty in extending seismic processing to anisotropic media is the recovery of anisotropic velocity fields from surface reflection data. We suggest carrying out velocity analysis for transversely isotropic (TI) media by inverting the dependence of P‐wave moveout velocities on the ray parameter. The inversion technique is based on the exact analytic equation for the normal‐moveout (NMO) velocity for dipping reflectors in anisotropic media. We show that P‐wave NMO velocity for dipping reflectors in homogeneous TI media with a vertical symmetry axis depends just on the zero‐dip value [Formula: see text] and a new effective parameter η that reduces to the difference between Thomsen parameters ε and δ in the limit of weak anisotropy. Our inversion procedure makes it possible to obtain η and reconstruct the NMO velocity as a function of ray parameter using moveout velocities for two different dips. Moreover, [Formula: see text] and η determine not only the NMO velocity, but also long‐spread (nonhyperbolic) P‐wave moveout for horizontal reflectors and the time‐migration impulse response. This means that inversion of dip‐moveout information allows one to perform all time‐processing steps in TI media using only surface P‐wave data. For elliptical anisotropy (ε = δ), isotropic time‐processing methods remain entirely valid. We show the performance of our velocity‐analysis method not only on synthetic, but also on field data from offshore Africa. Accurate time‐to‐depth conversion, however, requires that the vertical velocity [Formula: see text] be resolved independently. Unfortunately, it cannot be done using P‐wave surface moveout data alone, no matter how many dips are available. In some cases [Formula: see text] is known (e.g., from check shots or well logs); then the anisotropy parameters ε and δ can be found by inverting two P‐wave NMO velocities corresponding to a horizontal and a dipping reflector. If no well information is available, all three parameters ([Formula: see text], ε, and δ) can be obtained by combining our inversion results with shear‐wave information, such as the P‐SV or SV‐SV wave NMO velocities for a horizontal reflector. Generalization of the single‐layer NMO equation to layered anisotropic media with a dipping reflector provides a basis for extending anisotropic velocity analysis to vertically inhomogeneous media. We demonstrate how the influence of a stratified anisotropic overburden on moveout velocity can be stripped through a Dix‐type differentiation procedure.


Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 904-910 ◽  
Author(s):  
Vladimir Grechka ◽  
Andres Pech ◽  
Ilya Tsvankin ◽  
Baoniu Han

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.


Sign in / Sign up

Export Citation Format

Share Document