P-Wave Azimuthal AVO in a Carbonate Reservoir: An Integrated Seismic Anisotropy Study
- Mary Johns (ExxonMobil Upstream Research Co.) | David Y. Wang (ExxonMobil) | Shiyu Xu (ExxonMobil Upstream Research Co.) | Chih-Ping Lu (ExxonMobil) | Sam Sun (ExxonMobil) | Da Zhou (Exxon Mobil)
- Document ID
- Society of Petroleum Engineers
- SPE Reservoir Evaluation & Engineering
- Publication Date
- August 2008
- Document Type
- Journal Paper
- 719 - 725
- 2008. Society of Petroleum Engineers
- 4.3.4 Scale, 1.6.9 Coring, Fishing, 1.14 Casing and Cementing, 1.2.3 Rock properties, 5.6.1 Open hole/cased hole log analysis, 5.1.7 Seismic Processing and Interpretation, 3 Production and Well Operations, 7.2.2 Risk Management Systems, 3.2.3 Hydraulic Fracturing Design, Implementation and Optimisation, 5.1.6 Near-Well and Vertical Seismic Profiles, 5.1.8 Seismic Modelling, 5.8.7 Carbonate Reservoir, 5.1.2 Faults and Fracture Characterisation
- 2 in the last 30 days
- 527 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 5.00|
|SPE Non-Member Price:||USD 35.00|
Seismic anisotropy is sensitive to the intensity and orientation of fractures aligned by tectonic or local stress fields. Azimuthal anisotropy of amplitude vs. offset (AzAVO) is one common measure of seismic anisotropy for fracture detection. For this study, we calculated P-wave AzAVO attributes for a Lower Cretaceous limestone reservoir in east Texas. Our study tests AzAVO inversion for the relatively low-fold, fair data quality typical of some land surveys. We developed a workflow to integrate the AzAVO inversion with rock physics, forward seismic modeling, seismic-scale fault interpretation, image logs, and cores.
Our AzAVO inversion method was based on Rüger's equation, which describes the observed reflectivity (R) as a function of incidence angle (?i ) and acquisition azimuth (Øi) to equal a function of normal incidence reflectivity (A 0), amplitude vs. offset (AVO) slope (B 0), anisotropy magnitude (B 1), orientation of symmetry of anisotropy (Ø 0), and noise (ni ):
R(?i , Øi ) ˜ A 0 + [B 0 + B 1 cos2 (?i - Ø 0 )] sin2 ?i + ni
Inversion of Rüger's equation yielded four attribute volumes (corresponding to A 0 , B 0 , B 1, and Ø 0 ). Additionally, we proposed an additional data-quality parameter that describes the relative amount of the total reflection energy predicted by the model. Preinversion processing carefully preserved signal via random noise attenuation, superbinning, bandwidth balancing, and phase alignment.
In general, AzAVO inversion displayed geologically interpretable high-anisotropy anomalies. Noise, however, limited our confidence for quantitative fracture prediction. On the flanks of the anticline, data quality was high (greater than 50% fit of observed to predicted amplitude). Here the magnitude of the anisotropy (B 1) volume showed strong lineations parallel to northeast-southwest trending faults. Over the crest of the anticline, the inversion quality was degraded by an overburden effect, coincident with poor data quality (less than 50% fit). AzAVO orientation maps were noisy and in some places obscured by an acquisition footprint. Locally, however, AzAVO orientations were parallel to east-northeast fracture orientations from cores, maximum horizontal stress from image logs, and fast-velocity direction from a dipole-sonic log.
A network of open, connected fractures increases the permeability of many reservoirs, especially carbonate reservoirs, worldwide. Understanding the distribution of fractures helps us to predict the spatial variations in production rates and flow anisotropy. Geologic fracture characterization tends to be limited by the number of cores or image logs, as well as sampling biases inherent to measuring vertical fractures in vertical wells. Seismic anisotropy is a fieldwide fracture detection technique most sensitive to the open fractures that are important for flow.
We define seismic anisotropy as the directional response of any seismic attribute, such as travel time or amplitude. Theoretical models show that both travel time and amplitude vary with azimuth in rocks with vertically aligned fractures (Rüger 1996, 1998; Tsvankin 2001; Grechka et al. 1999). Azimuthal normal moveout (AzNMO) detects anisotropy by use of the normal moveout travel time attribute of a prestack common midpoint (CMP) gather. We focus on AzAVO, which detects anisotropy from the amplitude variation of a particular event, such as Primary-wave (P-wave) or converted shear-wave (S-wave). Previous seismic models have predicted that shear-wave data may have a clearer anisotropic response than P-wave data (Lynn and Thomsen 1990; Li 1997). This study, however, considered only P-wave anisotropy. Because of the wide availability of P-wave seismic data, we expect the development of P-wave anisotropic techniques to have a broad business impact.
There are several known limitations of azimuthal anisotropy for fracture detection. Seismic waves respond to open fractures, so they are relatively insensitive to fractures closed by stress or cementation. In a similar manner, if fracture density is low, fracture-induced anisotropy is too weak to detect. Fracture sets in multiple orientations may diminish the magnitude of the azimuthal anisotropy. Azimuthal anisotropy is enhanced for vertical fractures and reduced for more shallow dipping fractures. Steeply dipping beds may induce an azimuthal anisotropy. Anisotropy at a reservoir level may be influenced by the anisotropy of rocks that overlie it. Finally, like any seismic attribute, azimuthal attributes are sensitive to seismic data quality. Interpreting azimuthal attributes requires understanding the geologic sources of anisotropy, the effects of azimuthal and offset distribution, the assumptions of the data processing, and the effects of noise.
For this study, our objectives were the following:
- Develop AzAVO algorithms and workflows for fracture detection.
- Apply the workflows for a fractured carbonate reservoir.
- Test the workflows on a seismic survey with fair land P-wave seismic acquisition and data quality.
- Integrate the seismic anisotropy results with mapped structures, stacked data attribute volumes, dipole sonic logs, fractures in cores and image logs, and production data.
Insight gained from this study may apply to other fractured carbonate reservoirs worldwide.
|File Size||2 MB||Number of Pages||7|
Caine, J.S., Evans, J.P., and Forster, C.B. 1996. Fault Zone Architecture and Permeability Structure. Geology24 (11): 1025-1028.DOI:10.1130/0091-7613(1996)024<1025:FZAAPS>2.3.CO;2.
Gray, F.D., Roberts, G., and Head, K.J. 2002. Recent Advances in Determination ofFracture Strike and Crack Density From P-Wave Seismic Data. TheLeading Edge 21 (3): 280-285. DOI:10.1190/1.1463778.
Grechka, V., and Kachanov, M. 1999. Seismic Characterization of MultipleFracture Sets: Does Orthotropy Suffice? Geophysics 71 (3):D93-D105. DOI:10.1190/1.2196872.
Grechka, V., Tsvankin, I., and Cohen, J.K. 1999. Generalized DixEquation and Analytic Treatment of Normal-Moveout Velocity for AnisotropicMedia. Geophysical Prospecting 47 (2): 117-148.DOI:10.1046/j.1365-2478.1999.00120.x.
Kim, Y.-S., Peacock, D.C.P., and Sanderson, D.J. 2004. Fault Damage Zones.Journal of Structural Geology 26 (3): 503-517.DOI:10.1016/j.jsg.2003.08.002.
Li, X.-Y. 1997. Fractured ReservoirDelineation Using Multicomponent Seismic Data. GeophysicalProspecting 45 (1): 39-64.DOI:10.1046/j.1365-2478.1997.3200262.x.
Lynn, H.B. and Thomsen, L.A. 1990. Reflection Shear-Wave Data CollectedNear Principal Axes of Azimuthal Anisotropy. Geophysics 55(2): 147-156. DOI:10.1190/1.1442821.
Rüger, A. 1996. Reflection Coefficients and Azimuthal AVO Analysis inAnisotropic Media. PhD thesis, Colorado School of Mines, Leadville,Colorado.
Rüger, A. 1998. Variation ofP-Wave Reflectivity With Offset and Azimuth in Anisotropic Media.Geophysics 63 (3): 935-947. DOI:10.1190/1.1444405.
Tsvankin, I. 2001. Seismic Signatures and Analysis of Reflection Data inAnisotropic Media. New York: Pergamon Press.
Varela, I., Liu, E., and Chapman, M. 2006. Azimuthal Multiattribute Analysisfor Fracture Properties in Time-Lapse VSP Data. In SEG InternationalExposition and 76th Annual Meeting, Expanded Abstracts, 3436-3440. Tulsa:Society of Exploration Geophysicists.
Weeks, A.W. 1945. Balcones, Luling, and Mexia Fault Zones in Texas. AAPGBulletin 29: 1733-1737.
Zheng, Y., Todorovic-Marinic, D., and Larson, G. 2004. Seismic FractureDetection: Ambiguity and Practical Solution. SEG 74th Annual InternationalMeeting, Expanded Abstracts, 1575-1578. Tulsa: Society of ExplorationGeophysicists.