Approximate ray tracing for qP-waves in inhomogeneous layered media with weak structural anisotropy

Geophysics ◽  
2007 ◽  
Vol 72 (5) ◽  
pp. SM47-SM60 ◽  
Author(s):  
Kaveh Dehghan ◽  
Véronique Farra ◽  
Laurence Nicolétis
Keyword(s):  
Ray Tracing ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Structural Anisotropy ◽  
First Order ◽  
Wave Ray ◽  
Order Quantities ◽  
Paraxial Ray ◽  
Approximate Equations

We present approximate equations for qP-wave ray tracing and paraxial ray tracing in inhomogeneous layered weakly transversely isotropic (TI) media. Inside layers, the symmetry axis direction of the TI medium is allowed to vary continuously. Approximate equations are based on perturbation theory in which deviations of anisotropy from isotropy are considered to be first-order quantities. For qP-waves propagating in a TI medium, the approximate ray-tracing and paraxial ray-tracing equations depend on three parameters and two angles defining the direction of the symmetry axis. We also present the boundary conditions at interfaces for first-order rays and paraxial rays and compute reflection/transmission coefficients to the first order. All the quantities required for evaluation of the Green’s function are calculated to the first order, except the traveltime that is calculated to the second order. The accuracy of the presented algorithm is verified on simple models with respect to exact ray-tracing results. For anisotropy of about [Formula: see text], considered in the examples presented, the relative errors in traveltime and amplitude are less than [Formula: see text] and [Formula: see text], respectively. We then consider media where the symmetry axis is conformable with the structure, named structural transverse isotropy (STI). The method is applied to simple STI models and to a realistic overthrust model, showing significant differences with vertical transverse isotropy (VTI) models not only for traveltimes, but also for amplitudes.

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.

Transverse isotropy of the upper mantle in the vicinity of Pacific fracture zones

1969 ◽  
Vol 59 (1) ◽  
pp. 59-72
Author(s):  
Robert S. Crosson ◽  
Nikolas I. Christensen
Keyword(s):  
Upper Mantle ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Seismic Wave Propagation ◽  
Fracture Zones ◽  
Mantle Material ◽  
Transversely Isotropic Media ◽  
Horizontal Symmetry ◽  
Isotropic Media

Abstract Several recent investigations suggest that portions of the Earth's upper mantle behave anisotropically to seismic wave propagation. Since several types of anisotropy can produce azimuthal variations in Pn velocities, it is of particular geophysical interest to provide a framework for the recognition of the form or forms of anisotropy most likely to be manifest in the upper mantle. In this paper upper mantle material is assumed to possess the elastic properties of transversely isotropic media. Equations are presented which relate azimuthal variations in Pn velocities to the direction and angle of tilt of the symmetry axis of a transversely isotropic upper mantle. It is shown that the velocity data of Raitt and Shor taken near the Mendocino and Molokai fracture zones can be adequately explained by the assumption of transverse isotropy with a nearly horizontal symmetry axis.

Normal moveout from dipping reflectors in anisotropic media

Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 268-284 ◽  
Author(s):  
Ilya Tsvankin
Keyword(s):  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Weak Anisotropy ◽  
Anisotropic Models ◽  
Transversely Isotropic Media ◽  
Wide Range ◽  
Isotropic Media ◽  
The Difference

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.

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.

Inversion of azimuthally dependent NMO velocity in transversely isotropic media with a tilted axis of symmetry

Geophysics ◽  
2000 ◽  
Vol 65 (1) ◽  
pp. 232-246 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Wave Velocity ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Model Parameters ◽  
Symmetry Direction ◽  
Time Migration ◽  
Inversion Procedure ◽  
Tti Media

Just as the transversely isotropic model with a vertical symmetry axis (VTI media) is typical for describing horizontally layered sediments, transverse isotropy with a tilted symmetry axis (TTI) describes dipping TI layers (such as tilted shale beds near salt domes) or crack systems. P-wave kinematic signatures in TTI media are controlled by the velocity [Formula: see text] in the symmetry direction, Thomsen’s anisotropic coefficients ε and δ, and the orientation (tilt ν and azimuth β) of the symmetry axis. Here, we show that all five parameters can be obtained from azimuthally varying P-wave NMO velocities measured for two reflectors with different dips and/or azimuths (one of the reflectors can be horizontal). The shear‐wave velocity [Formula: see text] in the symmetry direction, which has negligible influence on P-wave kinematic signatures, can be found only from the moveout of shear waves. Using the exact NMO equation, we examine the propagation of errors in observed moveout velocities into estimated values of the anisotropic parameters and establish the necessary conditions for a stable inversion procedure. Since the azimuthal variation of the NMO velocity is elliptical, each reflection event provides us with up to three constraints on the model parameters. Generally, the five parameters responsible for P-wave velocity can be obtained from two P-wave NMO ellipses, but the feasibility of the moveout inversion strongly depends on the tilt ν. If the symmetry axis is close to vertical (small ν), the P-wave NMO ellipse is largely governed by the NMO velocity from a horizontal reflector Vnmo(0) and the anellipticity coefficient η. Although for mild tilts the medium parameters cannot be determined separately, the NMO-velocity inversion provides enough information for building TTI models suitable for time processing (NMO, DMO, time migration). If the tilt of the symmetry axis exceeds 30°–40° (e.g., the symmetry axis can be horizontal), it is possible to find all P-wave kinematic parameters and construct the anisotropic model in depth. Another condition required for a stable parameter estimate is that the medium be sufficiently different from elliptical (i.e., ε cannot be close to δ). This limitation, however, can be overcome by including the SV-wave NMO ellipse from a horizontal reflector in the inversion procedure. While most of the analysis is carried out for a single layer, we also extend the inversion algorithm to vertically heterogeneous TTI media above a dipping reflector using the generalized Dix equation. A synthetic example for a strongly anisotropic, stratified TTI medium demonstrates a high accuracy of the inversion (subject to the above limitations).

First-order ray tracing for qP waves in inhomogeneous, weakly anisotropic media

Geophysics ◽  
2005 ◽  
Vol 70 (6) ◽  
pp. D65-D75 ◽  
Author(s):  
Ivan Pšenčík ◽  
Véronique Farra
Keyword(s):  
Ray Tracing ◽  
Anisotropic Media ◽  
Order Accuracy ◽  
Weak Anisotropy ◽  
Higher Symmetry ◽  
First Order ◽  
Isotropic Media ◽  
Anisotropy Parameters ◽  
Order Quantities ◽  
Relative Errors

We propose approximate ray-tracing equations for qP-waves propagating in smooth, inhomogeneous, weakly anisotropic media. For their derivation, we use perturbation theory, in which deviations of anisotropy from isotropy are considered to be the first-order quantities. The proposed ray-tracing equations and corresponding traveltimes are of the first order. Accuracy of the traveltimes can be increased by calculating a secondorder correction along first-order rays. The first-order ray-tracing equations for qP-waves propagating in a general weakly anisotropic medium depend on only 15 weak-anisotropy parameters (generalization of Thomsen’s parameters). The equations are thus considerably simpler than the exact ray-tracing equations. For higher-symmetry anisotropic media the equations differ only slightly from equations for isotropic media. They can thus substitute for the traditional isotropic ray tracers used in seismic processing. For vanishing anisotropy, the first-order ray-tracing equations reduce to standard, exact ray-tracing equations for isotropic media. Numerical tests for configuration and models used in seismic prospecting indicate negligible dependence of accuracy of calculated traveltimes on inhomogeneity of the medium. For anisotropy of about 8%, considered in the examples presented, the relative errors of the traveltimes, including the second-order correction, are well under 0.05%; for anisotropy of about 20%, they do not exceed 0.3%.

Building tilted transversely isotropic depth models using localized anisotropic tomography with well information

Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. D27-D36 ◽  
Author(s):  
Andrey Bakulin ◽  
Marta Woodward ◽  
Dave Nichols ◽  
Konstantin Osypov ◽  
Olga Zdraveva
Keyword(s):  
Seismic Data ◽  
Model Building ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Synthetic Data ◽  
Rock Physics ◽  
Transversely Isotropic ◽  
Data Set ◽  
Anisotropic Models ◽  
Well Data

Tilted transverse isotropy (TTI) is increasingly recognized as a more geologically plausible description of anisotropy in sedimentary formations than vertical transverse isotropy (VTI). Although model-building approaches for VTI media are well understood, similar approaches for TTI media are in their infancy, even when the symmetry-axis direction is assumed known. We describe a tomographic approach that builds localized anisotropic models by jointly inverting surface-seismic and well data. We present a synthetic data example of anisotropic tomography applied to a layered TTI model with a symmetry-axis tilt of 45 degrees. We demonstrate three scenarios for constraining the solution. In the first scenario, velocity along the symmetry axis is known and tomography inverts for Thomsen’s [Formula: see text] and [Formula: see text] parame-ters. In the second scenario, tomography inverts for [Formula: see text], [Formula: see text], and velocity, using surface-seismic data and vertical check-shot traveltimes. In contrast to the VTI case, both these inversions are nonunique. To combat nonuniqueness, in the third scenario, we supplement check-shot and seismic data with the [Formula: see text] profile from an offset well. This allows recovery of the correct profiles for velocity along the symmetry axis and [Formula: see text]. We conclude that TTI is more ambiguous than VTI for model building. Additional well data or rock-physics assumptions may be required to constrain the tomography and arrive at geologically plausible TTI models. Furthermore, we demonstrate that VTI models with atypical Thomsen parameters can also fit the same joint seismic and check-shot data set. In this case, although imaging with VTI models can focus the TTI data and match vertical event depths, it leads to substantial lateral mispositioning of the reflections.

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.

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.

Azimuthal variations in scattering amplitudes induced by transverse isotropy

Geophysics ◽  
1991 ◽  
Vol 56 (10) ◽  
pp. 1584-1595 ◽  
Author(s):  
M. A. Pelissier ◽  
Anna Thomas‐Betts ◽  
Peter D. Vestergaard
Keyword(s):  
Seismic Waves ◽  
Symmetry Axis ◽  
Mode Conversion ◽  
Transverse Isotropy ◽  
Full Range ◽  
Transversely Isotropic ◽  
Anisotropy Axis ◽  
Scattering Coefficients ◽  
Transversely Isotropic Media ◽  
Isotropic Media

The study of amplitude variations of reflected and transmitted seismic waves due to anistropy has received considerable attention in recent years, but most investigations have concentrated on the effect of transverse isotropy with the symmetry axis either vertical or horizontal. The published results on the whole tend to exclude mode conversions. Amplitudes of all reflected and transmitted wave modes are addressed for [Formula: see text]-waves incident on boundaries between isotropic and transversely isotropic media, the symmetry axis of which is oriented at 45 degrees to the interface. The results cover the full range of incidence angles and all “acquisition azimuth” in the plane of the interface. When the anistropy axis is not normal to the interface, the scattering coefficients are shown to be highly dependent on the azimuth. The pattern of azimuthal variation is especially complicated in the case of mode conversion, and scattering coefficient profiles that are 180 degrees apart are not the same. This has the implication that source‐receiver interchangeability does not hold and could have serious consequences to amplitude studies in split‐spread surveys. Both the [Formula: see text] and [Formula: see text] reflections show strong azimuthal variations, dependent on both the dip and the strike of the anisotropy axis. It may be possible to recover shown that the scattered amplitude patterns are dominantly controlled by the value of the elastic modulus [Formula: see text].

Sign in / Sign up

Export Citation Format

Share Document