Estimation of anisotropy parameters in VTI media from surfaceP‐wave seismic data and VSP in Tarim Basin, China

Author(s):  
Xiaoming Luo ◽  
Yun Lin ◽  
Lin Wu ◽  
Fengchang Yao
Minerals ◽  
2021 ◽  
Vol 11 (3) ◽  
pp. 293
Author(s):  
Wei Tian ◽  
Xiaomin Li ◽  
Lei Wang

Disparities between fold amplitude (A) and intrusion thickness (Hsill) are critical in identifying elastic or inelastic deformation in a forced fold. However, accurate measurements of these two parameters are challenging because of the limit in separability and detectability of the seismic data. We combined wireline data and 3-D seismic data from the TZ-47 exploring area in the Tarim Basin, Northwest China, to accurately constrain the fold amplitude and total thickness of sills that induced roof uplift in the terrain. Results from the measurement show that the forced fold amplitude is 155.0 m. After decompaction, the original forced fold amplitude in the area penetrated by the well T47 ranged from 159.9 to 225.8 m, which overlaps the total thickness of the stack of sills recovered by seismic method (171.4 m) and well log method (181.0 m). Therefore, the fold amplitude at T47 area is likely to be elastic. In contrast, the outer area of the TZ-47 forced fold is characterized by shear-style deformation, indicating inelastic deformation at the marginal area. It is suggested that interbedded limestone layers would play an important role in strengthening the roof layers, preventing inelastic deformation during the emplacement of intrusive magma.


Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. C219-C227 ◽  
Author(s):  
Hanjie Song ◽  
Yingjie Gao ◽  
Jinhai Zhang ◽  
Zhenxing Yao

The approximation of normal moveout is essential for estimating the anisotropy parameters of the transversally isotropic media with vertical symmetry axis (VTI). We have approximated the long-offset moveout using the Padé approximation based on the higher order Taylor series coefficients for VTI media. For a given anellipticity parameter, we have the best accuracy when the numerator is one order higher than the denominator (i.e., [[Formula: see text]]); thus, we suggest using [4/3] and [7/6] orders for practical applications. A [7/6] Padé approximation can handle a much larger offset and stronger anellipticity parameter. We have further compared the relative traveltime errors between the Padé approximation and several approximations. Our method shows great superiority to most existing methods over a wide range of offset (normalized offset up to 2 or offset-to-depth ratio up to 4) and anellipticity parameter (0–0.5). The Padé approximation provides us with an attractive high-accuracy scheme with an error that is negligible within its convergence domain. This is important for reducing the error accumulation especially for deeper substructures.


Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. C229-C237 ◽  
Author(s):  
Shibo Xu ◽  
Alexey Stovas

The moveout approximations are commonly used in seismic data processing such as velocity analysis, modeling, and time migration. The anisotropic effect is very obvious for a converted wave when estimating the physical and processing parameters from the real data. To approximate the traveltime in an elastic orthorhombic (ORT) medium, we defined an explicit rational-form approximation for the traveltime of the converted [Formula: see text]-, [Formula: see text]-, and [Formula: see text]-waves. To obtain the expression of the coefficients, the Taylor-series approximation is applied in the corresponding vertical slowness for three pure-wave modes. By using the effective model parameters for [Formula: see text]-, [Formula: see text]-, and [Formula: see text]-waves, the coefficients in the converted-wave traveltime approximation can be represented by the anisotropy parameters defined in the elastic ORT model. The accuracy in the converted-wave traveltime for three ORT models is illustrated in numerical examples. One can see from the results that, for converted [Formula: see text]- and [Formula: see text]-waves, our rational-form approximation is very accurate regardless of the tested ORT model. For a converted [Formula: see text]-wave, due to the existence of cusps, triplications, and shear singularities, the error is relatively larger compared with PS-waves.


Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. B15-B32 ◽  
Author(s):  
Shaun Hadden ◽  
R. Gerhard Pratt ◽  
Brendan Smithyman

Anisotropic waveform tomography (AWT) uses anisotropic traveltime tomography followed by anisotropic full-waveform inversion (FWI). Such an approach is required for FWI in cases in which the geology is likely to exhibit anisotropy. An important anisotropy class is that of transverse isotropy (TI), and the special case of TI media with a vertical symmetry axis (VTI) media is often used to represent elasticity in undeformed sedimentary layering. We have developed an approach for AWT that uses an acoustic approximation to simulate waves in VTI media, and we apply this approach to crosshole data. In our approach, the best-fitting models of seismic velocity and Thomsen VTI anisotropy parameters are initially obtained using anisotropic traveltime tomography, and they are then used as the starting models for VTI FWI within the acoustic approximation. One common problem with the acoustic approach to TI media is the generation of late-arriving (spurious) S-waves as a by-product of the equation system. We used a Laplace-Fourier approach that effectively damps the spurious S-waves to suppress artifacts that might otherwise corrupt the final inversion results. The results of applying AWT to synthetic data illustrate the trade-offs in resolution between the two parameter classes of velocity and anisotropy, and they also verify anisotropic traveltime tomography as a valid method for generating starting models for FWI. The synthetic study further indicates the importance of smoothing the anisotropy parameters before proceeding to FWI inversions of the velocity parameter. The AWT technique is applied to real crosshole field gathers from a sedimentary environment in Western Canada, and the results are compared with the results from a simpler (elliptical) anisotropy model. The transversely isotropic approach yields an FWI image of the vertical velocity that (1) exhibits a superior resolution and (2) better predicts the field data than does the elliptical approach.


Geophysics ◽  
1997 ◽  
Vol 62 (6) ◽  
pp. 1855-1866 ◽  
Author(s):  
Jack K. Cohen

In their studies of transversely isotropic media with a vertical symmetry axis (VTI media), Alkhalifah and Tsvankin observed that, to a high numerical accuracy, the normal moveout (NMO) velocity for dipping reflectors as a function of ray parameter p depends mainly on just two parameters, each of which can be determined from surface P‐wave observations. They substantiated this result by using the weak‐anisotropy approximation and exploited it to develop a time‐domain processing sequence that takes into account vertical transverse isotropy. In this study, the two‐parameter Alkhalifah‐Tsvankin result was further examined analytically. It was found that although there is (as these authors already observed) some dependence on the remaining parameters of the problem, this dependence is weak, especially in the practically important regimes of weak to moderately strong transverse isotropy and small ray parameter. In each of these regimes, an analytic solution is derived for the anisotropy parameter η required for time‐domain P‐wave imaging in VTI media. In the case of elliptical anisotropy (η = 0), NMO velocity expressed through p is fully controlled just by the zero‐dip NMO velocity—one of the Alkhalifah‐ Tsvankin parameters. The two‐parameter representation of NMO velocity also was shown to be exact in another limit—that of the zero shear‐wave vertical velociy. The analytic results derived here are based on new representations for both the P‐wave phase velocity and normal moveout velocity in terms of the ray parameter, with explicit expressions given for the cases of vanishing onaxis shear speed, weak to moderate transverse isotropy, and small to moderate ray parameter. Using these formulas, I have rederived and, in some cases, extended in a uniform manner various results of Tsvankin, Alkhalifah, and others. Examples include second‐order expansions in the anisotropy parameters for both the P‐wave phase‐velocity function and NMO‐velocity function, as well as expansions in powers of the ray parameter for both of these functions. I have checked these expansions against the corresponding exact functions for several choices of the anisotropy parameters.


Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1920-1928 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin

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.


Sign in / Sign up

Export Citation Format

Share Document