QVOA analysis: P-wave attenuation anisotropy for fracture characterization

This paper investigates [Formula: see text]-anisotropy for characterizing fractured reservoirs — specifically, the variation of the seismic quality factor [Formula: see text] versus offset and azimuth (QVOA). We derive an analytical expression for P-wave attenuation in a transversely isotropic medium with horizontal symmetry axis (HTI) and provide a method (QVOA) for estimating fracture direction from azimuthally varying [Formula: see text] in PP-wave reflection data. The QVOA formula is similar to Rüger’s approximation for PP-wave reflection coefficients, the theoretical basis for amplitude variation with angle offset (AVOA) analysis. The technique for QVOA analysis is similar to azimuthal AVO analysis. We introduce two new seismic attributes: [Formula: see text] versus offset (QVO) gradient and intercept. QVO gradient inversion not only indicates fracture orientation but also characterizes [Formula: see text]-anisotropy. We relate the [Formula: see text]-anisotropy parameter [Formula: see text] to fractured-medium parameters and invert the QVO gradient to estimate [Formula: see text]. The attenuation parameter [Formula: see text] and Thomsen-style anisotropy parameter [Formula: see text] are found to be interdependent. The attenuation anisotropy magnitude strongly depends on the host rock’s [Formula: see text] parameter, whereas the dependence on fracture parameters is weak. This complicates the QVO gradient inversion for the fracture parameters. This result is independent of the attenuation mechanism. To illustrate the QVOA method in synthetic data, we use Hudson’s first-order effective-medium model of a dissipative fractured reservoir with fluid flow between aligned cracks and random pores as a possible mechanism for P-wave attenuation.

Download Full-text

P-wave attenuation anisotropy in TI media and its application in fracture parameters inversion

Applied Geophysics ◽

10.1007/s11770-016-0592-7 ◽

2016 ◽

Vol 13 (4) ◽

pp. 649-657 ◽

Cited By ~ 1

Author(s):

Yi-Yuan He ◽

Tian-Yue Hu ◽

Chuan He ◽

Yu-Yang Tan

Keyword(s):

Wave Attenuation ◽

P Wave ◽

Fracture Parameters ◽

Parameters Inversion ◽

Attenuation Anisotropy ◽

Ti Media

Download Full-text

Plane-wave attenuation anisotropy in orthorhombic media

Geophysics ◽

10.1190/1.2387137 ◽

2007 ◽

Vol 72 (1) ◽

pp. D9-D19 ◽

Cited By ~ 44

Author(s):

Yaping Zhu ◽

Ilya Tsvankin

Keyword(s):

Attenuation Coefficient ◽

Wave Attenuation ◽

Fractured Reservoirs ◽

The Body ◽

P Wave ◽

Attenuation Coefficients ◽

Velocity Anisotropy ◽

Homogeneous Wave ◽

Anisotropy Parameters ◽

Attenuation Anisotropy

Orthorhombic models are often used in the interpretation of azimuthally varying seismic signatures recorded over fractured reservoirs. Here, we develop an analytic framework for describing the attenuation coefficients in orthorhombic media with orthorhombic attenuation (i.e., the symmetry of both the real and imaginary parts of the stiffness tensor is identical) under the assumption of homogeneous wave propagation. The analogous form of the Christoffel equation in the symmetry planes of orthorhombic and VTI (transversely isotropic with a vertical symmetry axis) media helps to obtain the symmetry-plane attenuation coefficients by adapting the existing VTI equations. To take full advantage of this equivalence with transverse isotropy, we introduce a parameter set similar to the VTI attenuation-anisotropy parameters [Formula: see text], [Formula: see text], and [Formula: see text]. This notation, based on the same principle as Tsvankin’s velocity-anisotropy parameters for orthorhombic media, leads to concise linearized equations for thesymmetry-plane attenuation coefficients of all three modes (P, [Formula: see text], and [Formula: see text]).The attenuation-anisotropy parameters also allow us to simplify the P-wave attenuation coefficient [Formula: see text] outside the symmetry planes under the assumptions of small attenuation and weak velocity and attenuation anisotropy. The approximate coefficient [Formula: see text] has the same form as the linearized P-wave phase-velocity function, with the velocity parameters [Formula: see text] and [Formula: see text] replaced by the attenuation parameters [Formula: see text] and [Formula: see text]. The exact attenuation coefficient, however, also depends on the velocity-anisotropy parameters, while the body-wave velocities are almost unperturbed by the presence of attenuation. The reduction in the number of parameters responsible for the P-wave attenuation and the simple approximation for the coefficient [Formula: see text] provide a basis for inverting P-wave attenuation measurements from orthorhombic media. The attenuation processing must be preceded by anisotropic velocity analysis that can be performed (in the absence of pronounced velocity dispersion) using existing algorithms for nonattenuative media.

Download Full-text

Effects of fractures on NMO velocities and P‐wave azimuthal AVO response

Geophysics ◽

10.1190/1.1484514 ◽

2002 ◽

Vol 67 (3) ◽

pp. 711-726 ◽

Cited By ~ 16

Author(s):

Feng Shen ◽

Xiang Zhu ◽

M. Nafi Toksöz

Keyword(s):

Aspect Ratio ◽

Fractured Reservoirs ◽

Crack Density ◽

P Wave ◽

Fractured Reservoir ◽

Fracture Density ◽

Fracture Detection ◽

Fluid Content ◽

Fracture Parameters ◽

Fractured Medium

This paper attempts to explain the relationships between fractured medium properties and seismic signatures and distortions induced by geology‐related influences on azimuthal AVO responses. In the presence of vertically aligned fractures, the relationships between fracture parameters (fracture density, fracture aspect ratio, and saturated fluid content) and their seismic signatures are linked with rock physics models of fractured media. The P‐wave seismic signatures studied in this paper include anisotropic parameters (δ(v), (v), and γ(v)), NMO velocities, and azimuthal AVO responses, where δ(v) is responsible for near‐vertical P‐wave velocity variations, (v) defines P‐wave anisotropy, and γ(v) governs the degree of shearwave splitting. The results show that in gas‐saturated fractures, anisotropic parameters δ(v) and (v) vary with fracture density alone. However, in water‐saturated fractures δ(v) and (v) depend on fracture density and crack aspect ratio and are also related to Vp/VS and Vp of background rocks, respectively. Differing from δ(v) and (v), γ(v) is the parameter most related to crack density. It is insensitive to the saturated fluid content and crack aspect ratio. The P‐wave NMO velocities in horizontally layered media are a function of δ(v), and their properties are comparable with those of δ(v). Results from 3‐D finite‐difference modeling show that P‐wave azimuthal AVO variations do not necessarily correlate with the magnitude of fracture density. Our studies reveal that, in addition to Poisson's ratio, other elastic properties of background rocks have an effect on P‐wave azimuthal AVO variations. Varying the saturated fluid content of fractures can lead to azimuthal AVO variations and may greatly change azimuthal AVO responses. For a thin fractured reservoir, a tuning effect related to seismic wavelength and reservoir thickness can result in variations in AVO gradients and in azimuthal AVO variations. Results from instantaneous frequency and instantaneous bandwidth indicate that tuning can also lead to azimuthal variations in the rates of changes of the phase and amplitude of seismic waves. For very thin fractured reservoirs, the effect of tuning could become dominant. Our numerical results show that AVO gradients may be significantly distorted in the presence of overburden anisotropy, which suggests that the inversion of fracture parameters based on an individual AVO response would be biased unless this influence were corrected. Though P‐wave azimuthal AVO variations could be useful for fracture detection, the combination of other types of data is more beneficial than using P‐wave amplitude signatures alone, especially for the quantitative characterization of a fractured reservoir.

Download Full-text

P-wave attenuation anisotropy in fractured media: A seismic physical modelling study

Geophysical Prospecting ◽

10.1111/j.1365-2478.2012.01127.x ◽

2012 ◽

Vol 61 ◽

pp. 420-433 ◽

Cited By ~ 32

Author(s):

A.M. Ekanem ◽

J. Wei ◽

X.-Y. Li ◽

M. Chapman ◽

I.G. Main

Keyword(s):

Wave Attenuation ◽

Physical Modelling ◽

Fractured Media ◽

P Wave ◽

Attenuation Anisotropy ◽

Modelling Study

Download Full-text

Two-step joint PP- and PS-wave three-term amplitude-variation with offset inversion

Interpretation ◽

10.1190/int-2016-0056.1 ◽

2016 ◽

Vol 4 (4) ◽

pp. T613-T625 ◽

Cited By ~ 1

Author(s):

Qizhen Du ◽

Bo Zhang ◽

Xianjun Meng ◽

Chengfeng Guo ◽

Gang Chen ◽

...

Keyword(s):

Reflection Coefficient ◽

Wave Reflection ◽

Coefficient Matrix ◽

P Wave ◽

Inversion Method ◽

S Wave ◽

Amplitude Variation ◽

Amplitude Variation With Offset ◽

Avo Inversion ◽

Pp Wave

Three-term amplitude-variation with offset (AVO) inversion generally suffers from instability when there is limited prior geologic or petrophysical constraints. Two-term AVO inversion shows higher instability compared with three-term AVO inversion. However, density, which is important in the fluid-type estimation, cannot be recovered from two-term AVO inversion. To reliably predict the P- and S-waves and density, we have developed a robust two-step joint PP- and PS-wave three-term AVO-inversion method. Our inversion workflow consists of two steps. The first step is to estimate the P- and S-wave reﬂectivities using Stewart’s joint two-term PP- and PS-AVO inversion. The second step is to treat the P-wave reflectivity obtained from the first step as the prior constraint to remove the P-wave velocity related-term from the three-term Aki-Richards PP-wave approximated reflection coefficient equation, and then the reduced PP-wave reflection coefficient equation is combined with the PS-wave reflection coefficient equation to estimate the S-wave and density reﬂectivities. We determined the effectiveness of our method by first applying it to synthetic models and then to field data. We also analyzed the condition number of the coefficient matrix to illustrate the stability of the proposed method. The estimated results using proposed method are superior to those obtained from three-term AVO inversion.

Download Full-text

Seismic inversion for geologic fractures and fractured media

Geophysics ◽

10.1190/geo2016-0123.1 ◽

2017 ◽

Vol 82 (5) ◽

pp. C145-C161 ◽

Cited By ~ 3

Author(s):

Xiaoqin Cui ◽

Edward S. Krebes ◽

Laurence R. Lines

Keyword(s):

Seismic Data ◽

Wave Reflection ◽

Rock Properties ◽

Reflection And Transmission ◽

Transmission Coefficients ◽

Avo Inversion ◽

Impedance Contrast ◽

Fractured Medium ◽

Pp Wave ◽

Contact Part

Amplitude variation with offset (AVO) inversion attempts to use the available surface seismic data to estimate the density, P-wave velocity, and S-wave velocity of the earth model. Under linear slip interface theory, synthetic seismograms for models with fractures prove that fractures are also reflection generators. Consequently, observed reflections are not necessarily due to lithologic variations only, but they could be due in part to the effect of fractures. To obtain approximate equations for AVO inversion for fractured media, denoted by AVO with fracture (AVOF), we derived new equations for PP-wave reflection and transmission coefficients that are based on nonwelded contact boundary conditions. In particular, along with the fracture compliances, azimuth has also been taken into account in the equations because the fractures can have any orientation. The new approximate AVOF equations for a horizontally fractured medium with impedance contrast are developed by simplifying the equations for the new PP-wave reflection and transmission coefficients. In the new approximate AVOF equations, the reflection coefficients are divided into a welded contact part (a conventional impedance contrast part) and a nonwelded contact part (a fracture part). This makes the equations flexible enough to separately invert for the rock properties of the fracture and the background medium in the case of a fractured medium with impedance contrast. The new approximate AVOF equations state that fractures could cause the seismic reflectivity to be frequency dependent, and that the fractures not only influence the wave amplitude but also change the wave phase. The linear least-squares and nonlinear conjugate gradient inversion algorithms are applied to estimate the elastic reflectivity using the new approximate AVOF equations. The inverted results for seismic data for a horizontally fractured medium with impedance contrast are evaluated to find a more accurate delineation of the subsurface rock properties.

Download Full-text

Reflection coefficients in attenuative anisotropic media

Geophysics ◽

10.1190/1.3142874 ◽

2009 ◽

Vol 74 (5) ◽

pp. WB193-WB202 ◽

Cited By ~ 19

Author(s):

Jyoti Behura ◽

Ilya Tsvankin

Keyword(s):

Plane Wave ◽

Wave Reflection ◽

Symmetry Plane ◽

Anisotropic Media ◽

P Wave ◽

Substantial Contribution ◽

Reflection Coefficients ◽

Slowness Vector ◽

Anisotropy Parameters ◽

Attenuation Anisotropy

Such reservoir rocks as tar sands are characterized by significant attenuation and, in some cases, attenuation anisotropy. Most existing attenuation studies are focused on plane-wave attenuation coefficients, which determine the amplitude decay along the raypath of seismic waves. Here we study the influence of attenuation on PP- and PS-wave reflection coefficients for anisotropic media with the main emphasis on transversely isotropic models with a vertical symmetry axis (VTI). Concise analytic solutions obtained by linearizing the exact plane-wave reflection coefficients are verified by numerical modeling. To make a substantial contribution to reflection coefficients, attenuation must be strong, with the quality factor [Formula: see text] not exceeding 10. For such highly attenuative media, it is also necessary to take attenuation anisotropy into account if the magnitude of the Thomsen-styleattenuation-anisotropy parameters is relatively large. In general, the linearized reflection coefficients in attenuative media include velocity-anisotropy parameters but have almost “isotropic” dependence on attenuation. Our formalism also helps evaluate the influence of the inhomogeneity angle (the angle between the real and imaginary parts of the slowness vector) on the reflection coefficients. A nonzero inhomogeneity angle of the incident wave introduces additional terms into the PP- and PS-wave reflection coefficients, which makes conventional amplitude-variation-with-offset (AVO) analysis inadequate for strongly attenuative media. For instance, an incident P-wave with a nonzero inhomogeneity angle generates a mode-converted PS-wave at normal incidence, even if both half-spaces have a horizontal symmetry plane. The developed linearized solutions can be used in AVO inversion for highly attenuative (e.g., gas-sand and heavy-oil) reservoirs.

Download Full-text

Physical modeling and analysis of P-wave attenuation anisotropy in transversely isotropic media

Geophysics ◽

10.1190/1.2374797 ◽

2007 ◽

Vol 72 (1) ◽

pp. D1-D7 ◽

Cited By ~ 50

Author(s):

Yaping Zhu ◽

Ilya Tsvankin ◽

Pawan Dewangan ◽

Kasper van Wijk

Keyword(s):

Attenuation Coefficient ◽

Symmetry Axis ◽

Wave Attenuation ◽

Transversely Isotropic ◽

P Wave ◽

Velocity Variation ◽

Ratio Method ◽

Fracture Detection ◽

Wide Range ◽

Attenuation Anisotropy

Anisotropic attenuation can provide sensitive attributes for fracture detection and lithology discrimination. This paper analyzes measurements of the P-wave attenuation coefficient in a transversely isotropic sample made of phenolic material. Using the spectral-ratio method, we estimate the group (effective) attenuation coefficient of P-waves transmitted through the sample for a wide range of propagation angles (from [Formula: see text] to [Formula: see text]) with the symmetry axis. Correction for the difference between the group and phase angles and for the angular velocity variation help us to obtain the normalized phase attenuation coefficient [Formula: see text] governed by the Thomsen-style attenuation-anisotropy parameters [Formula: see text] and [Formula: see text]. Whereas the symmetry axis of the angle-dependent coefficient [Formula: see text] practically coincides with that of the velocity function, the magnitude of the attenuation anisotropy far exceeds that of the velocity anisotropy. The quality factor [Formula: see text] increases more than tenfold from the symmetry axis (slow direction) to the isotropy plane (fast direction). Inversion of the coefficient [Formula: see text] using the Christoffel equation yields large negative values of the parameters [Formula: see text] and [Formula: see text]. The robustness of our results critically depends on several factors, such as the availability of an accurate anisotropic velocity model and adequacy of the homogeneous concept of wave propagation, as well as the choice of the frequency band. The methodology discussed here can be extended to field measurements of anisotropic attenuation needed for AVO (amplitude-variation-with-offset) analysis, amplitude-preserving migration, and seismic fracture detection.

Download Full-text