Modeling sensitivity of 3D, 9-C wide azimuth data to changes in fluid content and crack density in cracked reservoirs

Geophysics ◽  
2010 ◽  
Vol 75 (5) ◽  
pp. T155-T165 ◽  
Author(s):  
Herurisa Rusmanugroho ◽  
George A. McMechan
Keyword(s):  
Practical Importance ◽  
Crack Density ◽  
Maximum Amplitude ◽  
P Wave ◽  
Staggered Grid ◽  
Wave Reflections ◽  
S Wave ◽  
Reservoir Properties ◽  
Fluid Content ◽  
Isotropic Media

The volume density of cracks and the fluids contained in them are salient aspects of characterization of cracked reservoirs. Thus, it is of practical importance to investigate whether variations in these reservoir properties are detectable in seismic observations. Eighth-order staggered-grid, 3D finite-difference simulations generate nine-component amplitude variations with offset and azimuth (AVOAZ) for reflections from the top of a vertically cracked zone embedded in an isotropic host. The T-matrix method is used to calculate elastic stiffness tensors. Responses for various crack densities and fluid contents show sensitivity to the spatial orientation of, and variation in, anisotropy. In isotropic media, when source and receiver components have the same orientation (such as XX and YY), reflection amplitude contours align approximately perpendicular to the particle motion. Mixed components (such as XY and YX) have amplitude patterns thatare symmetrical pairs of the same, or opposite, polarity on either side of the diagonal of the 9-C response matrix. In anisotropic media, AVOAZ data show the same basic patterns and symmetries as for isotropic media but with a superimposed tendency for alignment parallel to the strike of the vertical cracks. The data contain combined effects related to the source, receiver, and crack orientations. The sensitivity of data to changes in fluid content is quantified by comparing the differences between responses to various fluid conditions, to the maximum amplitude of oil-filled crack responses. For a crack density of 0.1, amplitude differences are [Formula: see text] for oil-dry and [Formula: see text] for oil-brine. The corresponding values for S-wave reflections are [Formula: see text] for oil-dry and [Formula: see text] for oil-brine. Amplitude changes caused by changing the oil-filled crack density from 0.1 to 0.2 are [Formula: see text] for P-wave reflections and [Formula: see text] for S-wave reflections. These differences are visible in AVOAZ data and are potentially diagnostic for reservoir characterization.

Petrophysical properties variation of bitumen-bearing carbonates at increasing temperatures from laboratory to model

Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. MR297-MR308
Author(s):  
Roberta Ruggieri ◽  
Fabio Trippetta
Keyword(s):  
Central Italy ◽  
Carbonate Reservoir ◽  
P Wave ◽  
Petrophysical Properties ◽  
S Wave ◽  
Reservoir Properties ◽  
Velocity Change ◽  
Seismic Properties ◽  
Wave Velocities ◽  
Velocity Changes

Variations in reservoir seismic properties can be correlated to changes in saturated-fluid properties. Thus, the determination of variation in petrophysical properties of carbonate-bearing rocks is of interest to the oil exploration industry because unconventional oils, such as bitumen (HHC), are emerging as an alternative hydrocarbon reserve. We have investigated the temperature effects on laboratory seismic wave velocities of HHC-bearing carbonate rocks belonging to the Bolognano Formation (Majella Mountain, central Italy), which can be defined as a natural laboratory to study carbonate reservoir properties. We conduct an initial characterization in terms of porosity and density for the carbonate-bearing samples and then density and viscosity measurements for the residual HHC, extracted by HCl dissolution of the hosting rock. Acoustic wave velocities are recorded from ambient temperature to 90°C. Our acoustic velocity data point out an inverse relationship with temperature, and compressional (P) and shear (S) wave velocities show a distinct trend with increasing temperature depending on the amount of HHC content. Indeed, samples with the highest HHC content show a larger gradient of velocity changes in the temperature range of approximately 50°C–60°C, suggesting that the bitumen can be in a fluid state. Conversely, below approximately 50°C, the velocity gradient is lower because, at this temperature, bitumen can change its phase in a solid state. We also propose a theoretical model to predict the P-wave velocity change at different initial porosities for HHC-saturated samples suggesting that the velocity change mainly is related to the absolute volume of HHC.

Seismic reflectivity of hydraulic fractures approximated by thin fluid layers

Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. T79-T87 ◽  
Author(s):  
A. Oelke ◽  
D. Alexandrov ◽  
I. Abakumov ◽  
S. Glubokovskikh ◽  
R. Shigapov ◽  
...  
Keyword(s):  
Analytical Solutions ◽  
Fluid Layer ◽  
Synthetic Data ◽  
P Wave ◽  
Hydraulic Fractures ◽  
Reflection Coefficients ◽  
S Wave ◽  
Reservoir Properties ◽  
Thin Fluid ◽  
Theoretical Results

We have analyzed the angle-dependent reflectivity of microseismic wavefields at a hydraulic fracture, which we modeled as an ideal thin fluid layer embedded in an elastic, isotropic solid rock. We derived full analytical solutions for the reflections of an incident P-wave, the P-P and P-S reflection coefficients, as well as for an incident S-wave, and the S-S and S-P reflection coefficients. The rather complex analytical solutions were then approximated and we found that these zero-thickness limit approximations are in good agreement with the linear slip model, representing a fracture at slip contact. We compared the analytical solutions for the P-P reflections with synthetic data that were derived using finite-difference modeling and found that the modeling confirmed our theoretical results. For typical parameters of microseismic monitoring by hydraulic fracturing, e.g., a layer thickness of [Formula: see text] and frequencies of [Formula: see text], the reflection coefficients depend on the Poisson’s ratio. Furthermore, the reflection coefficients of an incident S-wave are remarkably high. Theoretical results suggested that it is feasible to image hydraulic fractures using microseismic events as a source and to solve the inverse problem, that is, to interpret reflection coefficients extracted from microseismic data in terms of reservoir properties.

P-wave and S-wave azimuthal anisotropy at a naturally fractured gas reservoir, Bluebell‐Altamont Field, Utah

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1312-1328 ◽  
Author(s):  
Heloise B. Lynn ◽  
Wallace E. Beckham ◽  
K. Michele Simon ◽  
C. Richard Bates ◽  
M. Layman ◽  
...  
Keyword(s):  
Azimuthal Anisotropy ◽  
P Wave ◽  
Fracture Density ◽  
Wave Reflections ◽  
S Wave ◽  
Gas Reservoir ◽  
Green River ◽  
Wave Data ◽  
Reflection Data ◽  
Naturally Fractured

Reflection P- and S-wave data were used in an investigation to determine the relative merits and strengths of these two data sets to characterize a naturally fractured gas reservoir in the Tertiary Upper Green River formation. The objective is to evaluate the viability of P-wave seismic to detect the presence of gas‐filled fractures, estimate fracture density and orientation, and compare the results with estimates obtained from the S-wave data. The P-wave response to vertical fractures must be evaluated at different source‐receiver azimuths (travelpaths) relative to fracture strike. Two perpendicular lines of multicomponent reflection data were acquired approximately parallel and normal to the dominant strike of Upper Green River fractures as obtained from outcrop, core analysis, and borehole image logs. The P-wave amplitude response is extracted from prestack amplitude variation with offset (AVO) analysis, which is compared to isotropic‐model AVO responses of gas sand versus brine sand in the Upper Green River. A nine‐component vertical seismic profile (VSP) was also obtained for calibration of S-wave reflections with P-wave reflections, and support of reflection S-wave results. The direction of the fast (S1) shear‐wave component from the reflection data and the VSP coincides with the northwest orientation of Upper Green River fractures, and the direction of maximum horizontal in‐situ stress as determined from borehole ellipticity logs. Significant differences were observed in the P-wave AVO gradient measured parallel and perpendicular to the orientation of Upper Green River fractures. Positive AVO gradients were associated with gas‐producing fractured intervals for propagation normal to fractures. AVO gradients measured normal to fractures at known waterwet zones were near zero or negative. A proportional relationship was observed between the azimuthal variation of the P-wave AVO gradient as measured at the tops of fractured intervals, and the fractional difference between the vertical traveltimes of split S-waves (the “S-wave anisotropy”) of the intervals.

scholarly journals Seismic Attenuation Tomography From 2018 Lombok Earthquakes, Indonesia

2021 ◽  
Vol 9 ◽  
Author(s):  
Awali Priyono ◽  
Andri Dian Nugraha ◽  
Muzli Muzli ◽  
Ardianto Ardianto ◽  
Atin Nur Aulia ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Source Region ◽  
Bouguer Anomaly ◽  
High Ratio ◽  
Seismic Attenuation ◽  
P Wave ◽  
S Wave ◽  
Fluid Content ◽  
P Wave Velocity ◽  
Lombok Island

Local earthquake data was used to determine a three-dimensional (3D) seismic attenuation structure around the aftershock source region of the 2018 Lombok earthquake in Indonesia. The aftershocks were recorded by 13 seismic stations from August 4 to September 9, 2018. The selected data consist of 6,281 P-wave t∗ values from 914 events, which had good t∗ quality in at least four stations. Our results show that the two aftershock clusters northwest and northeast of Lombok Island have different attenuation characteristics. A low P-wave quality factor (low-Qp), low P-wave velocity (Vp), and high ratio of P-wave velocity and S-wave velocity (Vp/Vs), which coincide with a shallower earthquake (<20 km) northwest of Lombok Island, might be associated with a brittle area of basal and imbricated faults influenced by high fluid content. At the same time, the high-Qp, low Vp, and low Vp/Vs, which coincide with a deeper earthquake (>20 km) northeast of Lombok Island, might be associated with an area that lacks fluid content. The difference in fluid content between the northwest and northeast regions might be the cause of the early generation of aftershocks in the northwest area. The significant earthquake that happened on August 5, 2018, took place in a region with moderate Qp, close to the contrast of high and low-Qp and high Vp, which suggests that the earthquake started in a strong material before triggering the shallower aftershocks occurring in an area affected by fluid content. We also identified an old intrusive body on the northeast flank of the Rinjani volcano, which was characterized by a high-Qp, high-velocity, and a high Bouguer anomaly.

A rock-physics model to determine the pore microstructure of cracked porous rocks

2020 ◽  
Vol 223 (1) ◽  
pp. 622-631
Author(s):  
Lin Zhang ◽  
Jing Ba ◽  
José M Carcione
Keyword(s):  
Crack Density ◽  
Rock Physics ◽  
Differential Pressure ◽  
P Wave ◽  
Porous Rocks ◽  
S Wave ◽  
Wave Velocities ◽  
Pore Microstructure ◽  
Physics Model ◽  
Rock Physics Model

SUMMARY Determining rock microstructure remains challenging, since a proper rock-physics model is needed to establish the relation between pore microstructure and elastic and transport properties. We present a model to estimate pore microstructure based on porosity, ultrasonic velocities and permeability, assuming that the microstructure consists on randomly oriented stiff equant pores and penny-shaped cracks. The stiff pore and crack porosity varying with differential pressure is estimated from the measured total porosity on the basis of a dual porosity model. The aspect ratio of pores and cracks and the crack density as a function of differential pressure are obtained from dry-rock P- and S-wave velocities, by using a differential effective medium model. These results are used to invert the pore radius from the matrix permeability by using a circular pore model. Above a crack density of 0.13, the crack radius can be estimated from permeability, and below that threshold, the radius is estimated from P-wave velocities, taking into account the wave dispersion induced by local fluid flow between pores and cracks. The approach is applied to experimental data for dry and saturated Fontainebleau sandstone and Chelmsford Granite.

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.

The effects of transverse isotropy on reflection amplitude versus offset

Geophysics ◽  
1987 ◽  
Vol 52 (4) ◽  
pp. 564-567 ◽  
Author(s):  
J. Wright
Keyword(s):  
Wave Velocity ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
S Wave ◽  
Reflection Amplitude ◽  
Amplitude Versus Offset ◽  
Transversely Isotropic Media ◽  
P Wave Velocity ◽  
Isotropic Media

Studies have shown that elastic properties of materials such as shale and chalk are anisotropic. With the increasing emphasis on extraction of lithology and fluid content from changes in reflection amplitude with shot‐to‐group offset, one needs to know the effects of anisotropy on reflectivity. Since anisotropy means that velocity depends upon the direction of propagation, this angular dependence of velocity is expected to influence reflectivity changes with offset. These effects might be particularly evident in deltaic sand‐shale sequences since measurements have shown that the P-wave velocity of shales in the horizontal direction can be 20 percent higher than the vertical P-wave velocity. To investigate this behavior, a computer program was written to find the P- and S-wave reflectivities at an interface between two transversely isotropic media with the axis of symmetry perpendicular to the interface. Models for shale‐chalk and shale‐sand P-wave reflectivities were analyzed.

Detection of near-surface hydrocarbon seeps using P- and S-wave reflections

2016 ◽  
Vol 4 (3) ◽  
pp. SH21-SH37 ◽  
Author(s):  
Mathieu J. Duchesne ◽  
André J.-M. Pugin ◽  
Gabriel Fabien-Ouellet ◽  
Mathieu Sauvageau
Keyword(s):  
P Wave ◽  
Unconsolidated Sediments ◽  
Fluid Migration ◽  
Wave Reflections ◽  
S Wave ◽  
Abrupt Decrease ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
Wave Data ◽  
Hydrocarbon Seeps

The combined use of P- and S-wave seismic reflection data is appealing for providing insights into active petroleum systems because P-waves are sensitive to fluids and S-waves are not. The method presented herein relies on the simultaneous acquisition of P- and S-wave data using a vibratory source operated in the inline horizontal mode. The combined analysis of P- and S-wave reflections is tested on two potential hydrocarbon seeps located in a prospective area of the St. Lawrence Lowlands in Eastern Canada. For both sites, P-wave data indicate local changes in the reflection amplitude and slow velocities, whereas S-wave data present an anomalous amplitude at one site. Differences between P- and S-wave reflection morphology and amplitude and the abrupt decrease in P-velocity are indirect lines of evidence for hydrocarbon migration toward the surface through unconsolidated sediments. Surface-gas analysis made on samples taken at one potential seeping site reveals the occurrence of thermogenic gas that presumably vents from the underlying fractured Utica Shale forming the top of the bedrock. The 3C shear data suggest that fluid migration locally disturbs the elastic properties of the matrix. The comparative analysis of P- and S-wave data along with 3C recordings makes this method not only attractive for the remote detection of shallow hydrocarbons but also for the exploration of how fluid migration impacts unconsolidated geologic media.

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.

Separation of S‐wave and P‐wave reflections offshore western Florida

Geophysics ◽  
1984 ◽  
Vol 49 (5) ◽  
pp. 493-508 ◽  
Author(s):  
Robert H. Tatham ◽  
Donald V. Goolsbee
Keyword(s):  
Mode Conversion ◽  
Hard Water ◽  
P Wave ◽  
Angle Of Incidence ◽  
Wave Reflections ◽  
S Wave ◽  
Seismic Reflection Data ◽  
Reflection Time ◽  
Processing Techniques ◽  
Independent Confirmation

Hard water‐bottom marine environments, such as offshore western Florida, have presented particular problems in the acquisition and processing of seismic reflection data. One problem has been the limited angle of incidence (less than critical) available to P‐wave penetration into the subsurface. Mode conversion from P‐wave to S‐waves (SV), however, is quite efficient over a broad range of angles of incidence. After the success of a previously reported physical model experiment, an experimental line was acquired offshore western Florida. The 19 mile line, located approximately 100 miles west of Key West, Florida, was shot and processed. Three key factors have contributed to the successful recording of mode‐converted S‐wave reflections: (1) recognition of the effect of the group length on attenuation of energy arriving at large angles of incidence; (2) tau‐p processing techniques that allow separation of energy by angle of incidence; and (3) velocity filtering over a range of hyperbolic normal‐moveout (NMO) velocities as part of the forward tau‐p transform. These three factors, two of them data processing techniques, have allowed separation of P‐ and S‐wave energy in the marine environment. Overall, S‐wave reflections have been unambiguously identified to a reflection time of 2 sec and may be interpreted to a reflection time of 2 sec. Integrating an S‐wave section with P‐wave interpretations of offshore Florida data allows an independent confirmation of structural events. This independent confirmation may be more significant than improvements in the P‐wave data quality alone. Lateraly stable [Formula: see text] values are computed in intervals 1500 to 5000 ft thick and to S‐wave reflection times as great as 3 sec. The opportunity of [Formula: see text], interpretations for lithologic identification, gas thickness estimates, and general stratigraphic trap exploration makes mode‐converted shear waves a new tool in this area.

Sign in / Sign up

Export Citation Format

Share Document