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

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

QVOA analysis: P-wave attenuation anisotropy for fracture characterization

Geophysics ◽

10.1190/1.2194531 ◽

2006 ◽

Vol 71 (3) ◽

pp. C37-C48 ◽

Cited By ~ 48

Author(s):

Tatiana Chichinina ◽

Vladimir Sabinin ◽

Gerardo Ronquillo-Jarillo

Keyword(s):

Wave Reflection ◽

Wave Attenuation ◽

Fractured Reservoirs ◽

Anisotropy Parameter ◽

P Wave ◽

Fracture Parameters ◽

Fractured Medium ◽

Seismic Quality Factor ◽

Pp Wave ◽

Attenuation Anisotropy

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

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