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.

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 Phased array compaction cell for measurement of the transversely isotropic elastic properties of compacting sediments

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. WA113-WA123 ◽  
Author(s):  
Kurt T. Nihei ◽  
Seiji Nakagawa ◽  
Frederic Reverdy ◽  
Larry R. Myer ◽  
Luca Duranti ◽  
...  
Keyword(s):  
Elastic Properties ◽  
Phased Array ◽  
Transversely Isotropic ◽  
P Wave ◽  
S Wave ◽  
P Waves ◽  
Array Measurements ◽  
Experimental Apparatus ◽  
Marine Sediment Sample ◽  
Plane P Waves

Sediments undergoing compaction typically exhibit transversely isotropic (TI) elastic properties. We present a new experimental apparatus, the phased array compaction cell, for measuring the TI elastic properties of clay-rich sediments during compaction. This apparatus uses matched sets of P- and S-wave ultrasonic transducers located along the sides of the sample and an ultrasonic P-wave phased array source, together with a miniature P-wave receiver on the top and bottom ends of the sample. The phased array measurements are used to form plane P-waves that provide estimates of the phase velocities over a range of angles. From these measurements, the five TI elastic constants can be recovered as the sediment is compacted, without the need for sample unloading, recoring, or reorienting. This paper provides descriptions of the apparatus, the data processing, and an application demonstrating recovery of the evolving TI properties of a compacting marine sediment sample.

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.

scholarly journals An acoustic eikonal equation for attenuating transversely isotropic media with a vertical symmetry axis

Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. C9-C20 ◽  
Author(s):  
Qi Hao ◽  
Tariq Alkhalifah
Keyword(s):  
Eikonal Equation ◽  
Wave Attenuation ◽  
Transversely Isotropic ◽  
S Wave ◽  
P Waves ◽  
The Earth ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Vti Media ◽  
Complex Valued

Seismic-wave attenuation is an important component of describing wave propagation. Certain regions, such as gas clouds inside the earth, exert highly localized attenuation. In fact, the anisotropic nature of the earth induces anisotropic attenuation because the quasi P-wave dispersion effect should be profound along the symmetry direction. We have developed a 2D acoustic eikonal equation governing the complex-valued traveltime of quasi P-waves in attenuating, transversely isotropic media with a vertical-symmetry axis (VTI). This equation is derived under the assumption that the complex-valued traveltime of quasi P-waves in attenuating VTI media are independent of the S-wave velocity parameter [Formula: see text] in Thomsen’s notation and the S-wave attenuation coefficient [Formula: see text] in Zhu and Tsvankin’s notation. We combine perturbation theory and Shanks transform to develop practical approximations to the acoustic attenuating eikonal equation, capable of admitting an analytical description of the attenuation in homogeneous media. For a horizontal-attenuating VTI layer, we also derive the nonhyperbolic approximations for the real and imaginary parts of the complex-valued reflection traveltime. These equations reveal that (1) the quasi SV-wave velocity and the corresponding quasi SV-wave attenuation coefficient given as part of Thomsen-type notation barely affect the ray velocity and ray attenuation of quasi P-waves in attenuating VTI media; (2) combining the perturbation method and Shanks transform provides an accurate analytic eikonal solution for homogeneous attenuating VTI media; (3) for a horizontal attenuating VTI layer with weak attenuation, the real part of the complex-valued reflection traveltime may still be described by the existing nonhyperbolic approximations developed for nonattenuating VTI media, and the imaginary part of the complex-valued reflection traveltime still has the shape of nonhyperbolic curves. In addition, we have evaluated the possible extension of the proposed eikonal equation to realistic attenuating media, an alternative perturbation solution to the proposed eikonal equation, and the feasibility of applying the proposed nonhyperbolic equation for the imaginary part of the complex-valued traveltime to invert for interval attenuation parameters.

Finite-difference frequency-domain modeling of viscoacoustic wave propagation in 2D tilted transversely isotropic (TTI) media

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. T75-T95 ◽  
Author(s):  
Stéphane Operto ◽  
Jean Virieux ◽  
A. Ribodetti ◽  
J. E. Anderson
Keyword(s):  
Finite Difference ◽  
Frequency Domain ◽  
Symmetry Axis ◽  
Transversely Isotropic ◽  
Wave Mode ◽  
Absorbing Boundary Conditions ◽  
P Wave ◽  
Domain Modeling ◽  
S Wave ◽  
Absorbing Boundary

A 2D finite-difference, frequency-domain method was developed for modeling viscoacoustic seismic waves in transversely isotropic media with a tilted symmetry axis. The medium is parameterized by the P-wave velocity on the symmetry axis, the density, the attenuation factor, Thomsen’s anisotropic parameters [Formula: see text] and [Formula: see text], and the tilt angle. The finite-difference discretization relies on a parsimonious mixed-grid approach that designs accurate yet spatially compact stencils. The system of linear equations resulting from discretizing the time-harmonic wave equation is solved with a parallel direct solver that computes monochromatic wavefields efficiently for many sources. Dispersion analysis shows that four grid points per P-wavelength provide sufficiently accurate solutions in homogeneous media. The absorbing boundary conditions are perfectly matched layers (PMLs). The kinematic and dynamic accuracy of the method wasassessed with several synthetic examples which illustrate the propagation of S-waves excited at the source or at seismic discontinuities when [Formula: see text]. In frequency-domain modeling with absorbing boundary conditions, the unstable S-wave mode is not excited when [Formula: see text], allowing stable simulations of the P-wave mode for such anisotropic media. Some S-wave instabilities are seen in the PMLs when the symmetry axis is tilted and [Formula: see text]. These instabilities are consistent with previous theoretical analyses of PMLs in anisotropic media; they are removed if the grid interval is matched to the P-wavelength that leads to dispersive S-waves. Comparisons between seismograms computed with the frequency-domain acoustic TTI method and a finite-difference, time-domain method for the vertical transversely isotropic elastic equation show good agreement for weak to moderate anisotropy. This suggests the method can be used as a forward problem for viscoacoustic anisotropic full-waveform inversion.

Seismic data processing in vertically inhomogeneous TI media

Geophysics ◽  
1997 ◽  
Vol 62 (2) ◽  
pp. 662-675 ◽  
Author(s):  
Tariq Alkhalifah
Keyword(s):  
Seismic Data ◽  
Transversely Isotropic ◽  
Inhomogeneous Media ◽  
P Wave ◽  
Seismic Data Processing ◽  
Time Migration ◽  
Processing Data ◽  
Two Parameters ◽  
Ti Media ◽  
Ray Parameter

The first and most important step in processing data in transversely isotropic (TI) media for which velocities vary with depth is parameter estimation. The multilayer normal‐moveout (NMO) equation for a dipping reflector provides the basis for extending the TI velocity analysis of Alkhalifah and Tsvankin to vertically inhomogeneous media. This NMO equation is based on a root‐mean‐square (rms) average of interval NMO velocities that correspond to a single ray parameter, that of the dipping event. Therefore, interval NMO velocities [including the normal‐moveout velocity for horizontal events, [Formula: see text]] can be extracted from the stacking velocities using a Dix‐type differentiation procedure. On the other hand, η, which is a key combination of Thomsen's parameters that time‐related processing relies on, is extracted from the interval NMO velocities using a homogeneous inversion within each layer. Time migration, like dip moveout, depends on the same two parameters in vertically inhomogeneous media, namely [Formula: see text] and η, both of which can vary with depth. Therefore, [Formula: see text] and ε estimated using the dip dependency of P‐wave moveout velocity can be used for TI time migration. An application of anisotropic processing to seismic data from offshore Africa demonstrates the importance of considering anisotropy, especially as it pertains to focusing and imaging of dipping events.

Seismic vertical transversely isotropic parameter inversion from P- and S-wave cross-borehole measurements in an aquifer environment

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. D101-D116
Author(s):  
Julius K. von Ketelhodt ◽  
Musa S. D. Manzi ◽  
Raymond J. Durrheim ◽  
Thomas Fechner
Keyword(s):  
Transversely Isotropic ◽  
P Wave ◽  
Perturbation Equation ◽  
S Wave ◽  
P Waves ◽  
Data Set ◽  
Near Surface ◽  
S Waves ◽  
Wave Measurements ◽  
Wave Splitting

Joint P- and S-wave measurements for tomographic cross-borehole analysis can offer more reliable interpretational insight concerning lithologic and geotechnical parameter variations compared with P-wave measurements on their own. However, anisotropy can have a large influence on S-wave measurements, with the S-wave splitting into two modes. We have developed an inversion for parameters of transversely isotropic with a vertical symmetry axis (VTI) media. Our inversion is based on the traveltime perturbation equation, using cross-gradient constraints to ensure structural similarity for the resulting VTI parameters. We first determine the inversion on a synthetic data set consisting of P-waves and vertically and horizontally polarized S-waves. Subsequently, we evaluate inversion results for a data set comprising jointly measured P-waves and vertically and horizontally polarized S-waves that were acquired in a near-surface ([Formula: see text]) aquifer environment (the Safira research site, Germany). The inverted models indicate that the anisotropy parameters [Formula: see text] and [Formula: see text] are close to zero, with no P-wave anisotropy present. A high [Formula: see text] ratio of up to nine causes considerable SV-wave anisotropy despite the low magnitudes for [Formula: see text] and [Formula: see text]. The SH-wave anisotropy parameter [Formula: see text] is estimated to be between 0.05 and 0.15 in the clay and lignite seams. The S-wave splitting is confirmed by polarization analysis prior to the inversion. The results suggest that S-wave anisotropy may be more severe than P-wave anisotropy in near-surface environments and should be taken into account when interpreting cross-borehole S-wave data.

scholarly journals An acoustic eikonal equation for attenuating orthorhombic media

Geophysics ◽  
2017 ◽  
Vol 82 (4) ◽  
pp. WA67-WA81 ◽  
Author(s):  
Qi Hao ◽  
Tariq Alkhalifah
Keyword(s):  
Eikonal Equation ◽  
P Wave ◽  
S Wave ◽  
Seismic Wave Velocity ◽  
P Waves ◽  
Azimuthal Variation ◽  
Wave Parameters ◽  
Complex Valued ◽  
Complex Eikonal ◽  
Wave Complex

Attenuating orthorhombic models are often used to describe the azimuthal variation of the seismic wave velocity and attenuation in finely layered hydrocarbon reservoirs with vertical fractures. In addition to the P-wave related medium parameters, S-wave parameters are also present in the complex eikonal equation needed to describe the P-wave complex-valued traveltime in an attenuating orthorhombic medium, which increases the complexity of using the P-wave traveltime to invert for the medium parameters in practice. We have used the acoustic assumption to derive an acoustic eikonal equation that approximately governs the complex-valued traveltime of P-waves in an attenuating orthorhombic medium. For a homogeneous attenuating orthorhombic media, we solve the eikonal equation using a combination of the perturbation method and Shanks transform. For a horizontal attenuating orthorhombic layer, the real and imaginary parts of the complex-valued reflection traveltime have nonhyperbolic behaviors in terms of the source-receiver offset. Similar to the roles of normal moveout (NMO) velocity and anellipticity, the attenuation NMO velocity and the attenuation anellipticity characterize the variation of the imaginary part of the complex-valued reflection traveltime around zero source-receiver offset.

The joint nonhyperbolic moveout inversion of PP and PS data in VTI media

Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1929-1932 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Vertical Velocity ◽  
Single Layer ◽  
Transversely Isotropic ◽  
Accurate Estimate ◽  
P Wave ◽  
S Wave ◽  
P Waves ◽  
Layer Cake ◽  
Two Parameters ◽  
Vti Media

Nonhyperbolic moveout of P‐waves in horizontally layered transversely isotropic media with a vertical symmetry axis (VTI) can be used to estimate the anellipticity coefficient η in addition to the NMO velocity Vnmo,P. Those two parameters are sufficient for time processing of P‐wave data (despite a certain instability in the inversion for η), but they do not constrain the vertical velocity VP0 and the depth scale of the model. It has been suggested in the literature that this ambiguity in the depth‐domain velocity analysis for layer‐cake VTI media can be resolved by combining long‐spread reflection traveltimes of P‐waves and mode‐converted PSV‐waves. Here, we show that reflection traveltimes of horizontal PSV events help to determine the ratio of the P‐ and S‐wave vertical velocities and the NMO velocity of SV‐waves, and they give a more accurate estimate of η. However, nonhyperbolic moveout of PSV‐waves turns out to be mostly controlled by wide‐angle P‐wave traveltimes and does not provide independent information for the inversion. As a result, even for a single‐layer model and uncommonly large offsets, traveltimes of P‐ and PSV‐waves cannot be inverted for the vertical velocity and anisotropic parameters ε and δ. To reconstruct the horizontally layered VTI model from surface data, it is necessary to combine long‐spread traveltimes of pure P and SV reflections.

Sign in / Sign up

Export Citation Format

Share Document