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

Geophysics ◽  
2002 ◽  
Vol 67 (3) ◽  
pp. 711-726 ◽  
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.

QVOA analysis: P-wave attenuation anisotropy for fracture characterization

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. C37-C48 ◽  
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.

Model-based amplitude versus offset and azimuth inversion for estimating fracture parameters and fluid content

Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. M1-M17 ◽  
Author(s):  
Jiao Xue ◽  
Hanming Gu ◽  
Chengguo Cai
Keyword(s):  
Seismic Data ◽  
Shear Fracture ◽  
Influencing Factor ◽  
Fractured Reservoirs ◽  
Fracture Density ◽  
Inversion Algorithm ◽  
Fluid Content ◽  
Fracture Parameters ◽  
Amplitude Versus Offset ◽  
Amplitude Versus

The normal-to-shear fracture compliance ratio is commonly used as a fluid indicator. In the seismic frequency range, the fluid indicator lies between the values for isolated fluid-filled fractures and dry fractures, and it is not easy to discriminate the fluid content. Assuming that the fracture surfaces are smooth, we use [Formula: see text], with [Formula: see text] and [Formula: see text] representing the normal fracture weakness of the saturated and dry rock, to indicate fluid types, and to define a fluid influencing factor. The fluid influencing factor is sensitive to the fluid properties, the aspect ratio of the fractures, and the frequency. Conventionally, the amplitude versus offset and azimuth (AVOA) inversion is formulated in terms of the contrasts of the fracture weaknesses across the interface, assuming that the fractures are vertical with the same symmetry axis. We consider fractures with arbitrary azimuths, and develop a method to estimate fracture parameters from wide-azimuth seismic data. The proposed AVOA inversion algorithm is tested on real 3D prestack seismic data from the Tarim Basin, China, and the inverted fracture density show good agreement with well log data, except that there are some discrepancies for one of the fractured reservoir sections. The discrepancies can be ascribed to neglect of the dip angle for the tilted fractures and the conjugate fracture sets, and to the validity of the linear-slip model. The fractured reservoirs are expected to be liquid saturated, under the assumption of smooth fractures. Overall, the inverted fracture density and fluid influencing factor can be potentially used for better well planning in fractured reservoirs and quantitatively estimating the fluid effects.

Use of AVOA data to estimate fluid indicator in a vertically fractured medium

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. C15-C24 ◽  
Author(s):  
Ranjit K. Shaw ◽  
Mrinal K. Sen
Keyword(s):  
Reflection Coefficient ◽  
Random Noise ◽  
Crack Density ◽  
Synthetic Data ◽  
Quantitative Measure ◽  
P Wave ◽  
Fluid Content ◽  
Data Inversion ◽  
Microstructural Property ◽  
Fractured Medium

Microstructural attributes of cracks and fractures, such as crack density, aspect ratio, and fluid infill, determine the elastic properties of a medium containing a set of parallel, vertical fractures. Although the tangential weakness [Formula: see text] of the fractures does not vary with the fluid content, the normal weakness [Formula: see text] exhibits significant dependence on fluid infill. Based on linear-slip theory, we used the ratio [Formula: see text] — termed the fluid indicator — as a quantitative measure of the fluid content in the fractures, with g representing the square of the ratio of S- and P-wave velocity in the unfractured medium. We used a Born formalism to derive the sensitivity to fracture weakness of PP- and PS-reflection coefficients for an interface separating an unfractured medium from a vertically fractured medium. Our formulae reveal that the PP-reflection coefficient does not depend on the 2D microcorrugation/surface roughness with ridges and valleys parallel to the fracture strike, whereas the PS-reflection coefficient is sensitive to this microstructural property of the fractures. Based on this formulation, we developed a method to compute the fluid indicator from wide-azimuth PP-AVOA data. Inversion of synthetic data corrupted with 10% random noise reliably estimates the normal and tangential fracture weaknesses and hence the fluid indicator can be determined accurately when the fractures are liquid-filled or partially saturated. As the gas saturation in the fractures increases, the quality of inversion becomes poorer. Errors of 15%–20% in g do not affect the estimation of fluid indicator significantly in case of liquid infill or partial saturation. However, for gas-saturated fractures, incorrect values of g may have a significant effect on fluid-indicator estimates.

Amplitude variation with offset and azimuth inversion for fluid indicator and fracture weaknesses in an oil-bearing fractured reservoir

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. N41-N53 ◽  
Author(s):  
Xinpeng Pan ◽  
Guangzhi Zhang ◽  
Xingyao Yin
Keyword(s):  
Scattering Theory ◽  
Fractured Reservoirs ◽  
Fractured Reservoir ◽  
Amplitude Variation ◽  
Data Set ◽  
Amplitude Variation With Offset ◽  
Fluid Identification ◽  
Background Medium ◽  
Fracture Parameters ◽  
Fractured Medium

Fluid identification and fracture discrimination play an important role in the exploration and development of an oil-bearing fractured reservoir. The most common fluid indicator in fractured reservoirs, the normal-to-shear fracture compliance ratio, is influenced by the fluid content and the fracture intensity. To reduce the ambiguities in the discrimination of the fluid and fracture parameters, we have aimed to extend the scattering theory to implement the fluid identification and fracture detection by incorporating the azimuthal data in an oil-bearing fractured reservoir via the proposed Bayesian amplitude variation with offset and azimuth (AVOAz) inversion approach. The background medium is, as far as the scattering theory is concerned, an isotropic medium without fractures, and the fractured medium is corresponding to a perturbed medium. The elastic parameters of a saturated anisotropic medium can be parameterized as a perturbation over a homogeneous isotropic background medium. We used the scattering theory to derive a generalized AVOAz approximation that provided the iterative estimates of hydrocarbon fluid indicator, shear modulus, density, and fracture weaknesses in a Bayesian scheme. The inversion algorithm is based on a convolutional model and a weak-contrast and small-weakness PP-wave reflection coefficient. The approach is applied to an oil-bearing field data set from a fractured marlstone reservoir. We observe that reasonable estimates of fluid indicator and fracture weaknesses are inverted, which can be used to perform the discrimination of fluid and fracture parameters. We conclude that the proposed approach provides us a potentially powerful tool to estimate the reservoir fluid and fracture properties in a more straightforward and efficient manner than those previous methods.

Fracture parameter estimation from well-log data

Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. D129-D134 ◽  
Author(s):  
Samik Sil
Keyword(s):  
Wave Velocity ◽  
P Wave ◽  
Well Log ◽  
Fracture Parameter ◽  
Fracture Density ◽  
Log Data ◽  
Vertical Well ◽  
Fracture Parameters ◽  
P Wave Velocity ◽  
Fractured Medium

We evaluated a method of deriving seismic fracture parameters from vertical-well-log data with the assumption that the fractured medium is transversely isotropic with a horizontal axis of symmetry (HTI). One approximation we used is that the observed vertical P-wave velocity is the same as the background isotropic P-wave velocity of the HTI medium. Another assumption was that the fractures and cracks are noninteractive and penny shaped. Using these approximations, we generated the fracture compliance matrix for each layer. Fracture parameters were then derived by constructing the HTI stiffness matrix for those layers. We tested our method using vertical-well-log data from a tight sand reservoir in Colorado, USA. “Thomsen-style” parameters were derived, and gas-filled fractures were identified on this log. The identified gas-filled fractures were compared to the production log data. The fracture density was also obtained at the well location within the depth of interest. We also found some problems and limitations caused by approximating vertical P-wave velocity the same as the background isotropic P-wave velocity.

Multiscale fracture density analysis at Stromboli Volcano, Italy: implications to flank stability

2021 ◽  
Author(s):  
Thomas Alcock ◽  
Sergio Vinciguerra ◽  
Phillip Benson ◽  
Federico Vagnon
Keyword(s):  
Electrical Resistivity ◽  
Wave Velocity ◽  
Rift Zone ◽  
Pulse Generator ◽  
Average Length ◽  
Crack Density ◽  
P Wave ◽  
Fracture Density ◽  
Stromboli Volcano ◽  
P Wave Velocity

<p>Stromboli volcano has experienced four sector collapses over the past 13 thousand years, resulting in the formation of the Sciara del Fuoco (SDF) horseshoe-shaped depression and an inferred NE / SW striking rift zone across the SDF and the western sector of the island. These events have resulted in the formation of steep depressions on the slopes on the volcano where episodes of instability are continuously being observed and recorded. This study aims to quantify the fracture density inside and outside the rift zone to identify potential damaged zones that could reduce the edifice strength and promote fracturing. In order to do so we have carried out a multiscale analysis, by integrating satellite observations, field work and seismic and electrical resistivity analyses on cm scales blocks belonging to 11 lava units from the main volcanic cycles that have built the volcano edifice, ie. Paleostromboli, Nestromboli and Vancori. 0.5 m resolution Pleiades satellite data has been first used to highlight 23635 distinct linear features across the island. Fracture density has been calculated using Fracpaq based on the Mauldon et al (2001) method to determine the average fracture density of a given area on the basis of the average length of drawn segments within a predetermined circular area. 41.8 % of total fracture density is found around intrusions and fissures, with the summit area and the slopes of SDF having the highest average fracture density of 5.279  . Density, porosity, P- wave velocity in dry and wet conditions and electrical resistivity (in wet conditions) were measured  via an ultrasonic pulse generator and acquisition system (Pundit) and an on purpose built measuring quadrupole on cm scale blocks of lavas collected from both within and outside the proposed rift zone to assess the physical state and the crack damage of the different lava units.  Preliminary results show that P-wave velocity between ~ 2.25 km/s < Vp < 5km/s decreases with porosity while there is high variability electrical resistivity with 21.7 < ρ < 590 Ohm * m. This is presumably due to the lavas texture and the variable content of bubble/vesicles porosity and crack damage, that is reflected by an effective overall porosity between 0 and 9 %. Higher porosity is generally mirrored by lower p-wave velocity values. Neostromboli blocks show the most variability in both P-wave velocity and electrical resistivity. Further work will assess crack density throughout optical analyses and systematically investigate the UCS and elastic moduli. This integrated approach is expected to provide a multiscale fracture density and allow to develop further laboratory testing on how slip surfaces can evolve to a flank collapse at Stromboli.</p>

scholarly journals Analysis of fluid substitution in a porous and fractured medium

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. WA157-WA166 ◽  
Author(s):  
Samik Sil ◽  
Mrinal K. Sen ◽  
Boris Gurevich
Keyword(s):  
Water Saturation ◽  
Transversely Isotropic ◽  
P Wave ◽  
Reflection Coefficients ◽  
Fluid Substitution ◽  
Fracture Density ◽  
Maximum Increase ◽  
Amplitude Variation With Offset ◽  
Fractured Medium ◽  
Medium Change

To improve quantitative interpretation of seismic data, we analyze the effect of fluid substitution in a porous and fractured medium on elastic properties and reflection coefficients. This analysis uses closed-form expressions suitable for fluid substitution in transversely isotropic media with a horizontal symmetry axis (HTI). For the HTI medium, the effect of changing porosity and water saturation on (1) P-wave moduli, (2) horizontal and vertical velocities, (3) anisotropic parameters, and (4) reflection coefficients are examined. The effects of fracture density on these four parameters are also studied. For the model used in this study, a 35% increase in porosity lowers the value of P-wave moduli by maximum of 45%. Consistent with the reduction in P-wave moduli, P-wave velocities also decrease by maximum of 17% with a similar increment in porosity. The reduction is always larger for the horizontal P-wave modulus than for the vertical one and is nearly independent of fracture density. The magnitude of the anisotropic parameters of the fractured medium also changes with increased porosity depending on the changes in the value of P-wave moduli. The reflection coefficients at an interface of the fractured medium with an isotropic medium change in accordance with the above observations and lead to an increase in anisotropic amplitude variation with offset (AVO) gradient with porosity. Additionally, we observe a maximum increase in P-wave modulus and velocity by 30% and 8%, respectively, with a 100% increase in water saturation. Water saturation also changes the anisotropic parameters and reflection coefficients. Increase in water saturation considerably increases the magnitude of the anisotropic AVO gradient irrespective of fracture density. From this study, we conclude that porosity and water saturation have a significant impact on the four studied parameters and the impacts are seismically detectable.

Analysis of multiazimuthal VSP data for anisotropy and AVO

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1172-1180 ◽  
Author(s):  
W. Scott Leaney ◽  
Colin M. Sayers ◽  
Douglas E. Miller
Keyword(s):  
Fractured Reservoirs ◽  
Vertical Axis ◽  
Horizontal Stress ◽  
Azimuthal Anisotropy ◽  
Carbonate Reservoir ◽  
P Wave ◽  
Fractured Reservoir ◽  
Data Set ◽  
Vsp Data ◽  
The Impact

Multioffset vertical seismic profile (VSP) experiments, commonly referred to as walkaways, enable anisotropy to be measured reliably in the field. The results can be fed into modeling programs to study the impact of anisotropy on velocity analysis, migration, and amplitude versus offset (AVO). Properly designed multioffset VSPs can also provide the target AVO response measured under optimum conditions, since the wavelet is recorded just above the reflectors of interest with minimal reflection point dispersal. In this paper, the multioffset VSP technique is extended to include multioffset azimuths, and a multiazimuthal multiple VSP data set acquired over a carbonate reservoir is analyzed for P-wave anisotropy and AVO. Direct arrival times down to the overlying shale and reflection times and amplitudes from the carbonate are analyzed. Data analysis involves a three‐term fit to account for nonhyperbolic moveout, dip, and azimuthal anisotropy. Results indicate that the overlying shale is transversely isotropic with a vertical axis of symmetry (VTI), while the carbonate shows 4–5% azimuthal anisotropy in traveltimes. The fast direction is consistent with the maximum horizontal stress orientation determined from break‐out logs and is also consistent with the strike of major faults. AVO analysis of the reflection from the top of the carbonate layer shows a critical angle reduction in the fast direction and maximum gradient in the slow direction. This agrees with modeling and indicates a greater amplitude sensitivity in the slow direction—the direction perpendicular to fracture strike. In principle, 3-D surveys should have wide azimuthal coverage to characterize fractured reservoirs. If this is not possible, it is important to have azimuthal line coverage in the minimum horizontal stress direction to optimize the use of AVO for fractured reservoir characterization. This direction can be obtained from multiazimuthal walkaways using the azimuthal P-wave analysis techniques presented.

Quartic moveout coefficient for a dipping azimuthally anisotropic layer

Geophysics ◽  
2004 ◽  
Vol 69 (3) ◽  
pp. 699-707 ◽  
Author(s):  
Andrés Pech ◽  
Ilya Tsvankin
Keyword(s):  
Reservoir Characterization ◽  
Fractured Reservoirs ◽  
Symmetry Plane ◽  
Exact Expression ◽  
Naturally Fractured Reservoirs ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
Fracture Density ◽  
Weak Anisotropy ◽  
Velocity Anisotropy

Interpretation and inversion of azimuthally varying nonhyperbolic reflection moveout requires accounting for both velocity anisotropy and subsurface structure. Here, our previously derived exact expression for the quartic moveout coefficient A4 is applied to P‐wave reflections from a dipping interface overlaid by a medium of orthorhombic symmetry. The weak‐anisotropy approximaton for the coefficient A4 in a homogeneous orthorhombic layer is controlled by the anellipticity parameters η(1), η(2), and η(3), which are responsible for time processing of P‐wave data. If the dip plane of the reflector coincides with the vertical symmetry plane [x1, x3], A4 on the dip line is proportional to the in‐plane anellipticity parameter η(2) and always changes sign for a dip of 30○. The quartic coefficient on the strike line is a function of all three η–parameters, but for mild dips it is mostly governed by η(1)—the parameter defined in the incidence plane [x2, x3]. Whereas the magnitude of the dip line A4 typically becomes small for dips exceeding 45○, the nonhyperbolic moveout on the strike line may remain significant even for subvertical reflectors. The character of the azimuthal variation of A4 depends on reflector dip and is quite sensitive to the signs and relative magnitudes of η(1), η(2), and η(3). The analytic results and numerical modeling show that the azimuthal pattern of the quartic coefficient can contain multiple lobes, with one or two azimuths of vanishing A4 between the dip and strike directions. The strong influence of the anellipticity parameters on the azimuthally varying coefficient A4 suggests that nonhyperbolic moveout recorded in wide‐azimuth surveys can help to constrain the anisotropic velocity field. Since for typical orthorhombic models that describe naturally fractured reservoirs the parameters η(1,2,3) are closely related to the fracture density and infill, the results of azimuthal nonhyperbolic moveout analysis can also be used in reservoir characterization.

scholarly journals Physical and numerical modeling of tuning and diffraction in azimuthally anisotropic media

Geophysics ◽  
2000 ◽  
Vol 65 (5) ◽  
pp. 1613-1621 ◽  
Author(s):  
Richard L. Gibson ◽  
Stephen Theophanis ◽  
M. Nafi Toksöz
Keyword(s):  
Numerical Modeling ◽  
Fractured Reservoirs ◽  
Homogeneous Model ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
Fractured Reservoir ◽  
Sh Wave ◽  
Wave Data ◽  
Ray Methods

Fractured reservoirs are an important target for exploration and production geophysics, and the azimuthal anisotropy often associated with these reservoirs can strongly influence seismic wave propagation. We created a physical model of a fractured reservoir to simulate some of these propagation effects. The reservoir is represented by a phenolite disk that is thin with respect to the elastic wavelengths in the experiment, creating model dimensions that are representative of realistic reservoirs. Phenolite is strongly anisotropic with orthorhombic symmetry, which suggests that azimuthal amplitude versus offset (AVO) effects should be obvious in data. We acquired both SH- and P-wave data in common‐offset gathers with a near offset and a far offset and found that although the SH-wave data show clear azimuthal variations in AVO, the P-wave signals show no apparent changes with azimuth. We then applied numerical modeling to analyze the data. Because ray methods cannot model diffractions from the disk edge, we first used a ray‐Born technique to simulate variations in waveforms associated with such scattering. The synthetic seismograms reproduced variations in the SH-wave waveforms accurately, though the amplitude contrast between acquisition azimuths was overestimated. Assuming a laterally homogeneous model, we then applied ray methods to simulate tuning effects in SH- and P-wave data and confirmed that in spite of the large contrasts in elastic properties, the tuning of the P-wave reflections from the thin disk changed so there was negligible contrast in AVO with azimuth. Models of field scale reservoirs showed that the same effects could be expected for field applications.

Sign in / Sign up

Export Citation Format

Share Document