Euler Curvature Can Be a Calculated Success

Contributors: Satinder Chopra, Kurt Marfurt

Several (12, we think) types of seismic-based curvature attributes have been introduced the last few years – and of these, the most-positive and the most-negative curvatures described in last month’s article are the most popular.

Most-positive and most-negative curvatures provide more continuous maps of faults and flexures than do maximum and minimum curvatures, because the latter tend to rapidly change algebraic sign at fault and flexure intersections.

Other attributes, such as mean curvature, Gaussian curvature and shape index, also have been used by a few practitioners.

We describe here a technique called Euler curvature, which has valuable applications.

An attraction of Euler curvature is that it can be calculated in any desired azimuth across a 3-D volume to enhance the definition of specific lineaments. When this apparent curvature (the Euler curvature) is computed in several specific azimuths, the results are quite useful for interpreting azimuth-dependent structure.

The flow diagram in figure 1 explains the method for computing azimuth-dependent Euler curvature.

Applications

Mapping the intensities of fracture sets has been a major objective of reflection seismologists. Curvature, acoustic impedance and reflection coherence currently are the most effective attributes used to predict fractures in post-stack seismic data.

We describe here the application of Euler curvature to a 3-D seismic volume from northeast British Columbia, Canada. We use an interactive workflow to utilize Euler curvature much as we do in generating a suite of shaded relief maps.

Figure 2 shows 3-D chair displays through volumes of Euler curvature calculated at azimuths of 0, 45, 90 and 135 degrees from north. The left column shows long-wavelength curvature calculations, and the right column displays short-wavelength calculations.

Notice how east-west lineaments stand out when north-south curvature is estimated (azimuth = 0):

  • When curvature is estimated in an azimuth of 45 degrees, northwest-southeast lineaments are pronounced.
  • When east-west curvature is calculated (azimuth = 90 degrees), north-south features events are emphasized.
  • When northwest-southeast curvature is estimated (azimuth = 135 degrees), events slightly inclined away from north-south are better defined.

The analysis area shown in these figures spans approximately 100 square kilometers.

As emphasized in last month’s article, short-wavelength displays show more lineament detail and resolution than do long-wavelength displays. That principle is illustrated again by the displays in figure 2.

The important concept presented here is that there are obvious advantages in calculating Euler curvature on post-stack seismic volumes, because azimuth directions of curvature can be chosen to highlight lineaments in preferred directions.

Conclusions

Euler curvatures calculated in desired azimuthal directions produce better definitions of targeted lineaments.

Depending on the desired level of detail, either long- wavelength or short-wavelength estimates can be calculated. Short-wavelength Euler curvature would be more beneficial for observing fracture lineaments.

This work is in progress, and we hope to calibrate seismic-based lineaments determined with this technology with lineaments interpreted from image logs.

Acknowledgments

We thank Arcis Corporation for permission to show the data examples, as well as for the permission to publish this work.

Comments (0)

 

Geophysical Corner

Geophysical Corner - Satinder Chopra
Satinder Chopra, award-winning chief geophysicist (reservoir), at Arcis Seismic Solutions, Calgary, Canada, and a past AAPG-SEG Joint Distinguished Lecturer began serving as the editor of the Geophysical Corner column in 2012.

Geophysical Corner - Kurt Marfurt
AAPG member Kurt J. Marfurt is with the University of Oklahoma, Norman, Okla.

Geophysical Corner

The Geophysical Corner is a regular column in the EXPLORER that features geophysical case studies, techniques and application to the petroleum industry.

VIEW COLUMN ARCHIVES

Image Gallery

Part 2 of 3

This month's column, part two of a three-part series comparing structural and amplitude curvatures, deals with observing fault and fracture lineaments. See: Part 1

See Also: Book

Desktop /Portals/0/images/_site/AAPG-newlogo-vertical-morepadding.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 4588 Book

See Also: Bulletin Article

Prolific hydrocarbon discoveries in the subsalt, commonly known as the “presalt,” section of Brazil and the conjugate African margin have created a business imperative to predict reservoir quality in lacustrine carbonates. Geothermal convection is a style of groundwater flow known to occur in rift settings, which is capable of diagenetic modification of reservoir quality. We simulated variable density groundwater flow coupled with chemical reactions to evaluate the potential for diagenesis driven by convection in subsalt carbonates.

Rates of calcite diagenesis are critically controlled by temperature gradient and fluid flux following the principles of retrograde solubility. Simulations predict that convection could operate in rift carbonates prior to salt deposition, but with rates of dissolution in the reservoir interval only on the order of 0.01 vol. %/m.y., which is too low to significantly modify reservoir quality. The exception is around permeable fault zones and/or unconformities where flow is focused and dissolution rates are amplified to 1 to 10 vol. %/m.y. and could locally modify reservoir quality. After salt deposition, simulations also predict convection with a critical function for salt rugosity. The greatest potential for dissolution at rates of 0.1 to 1 vol. %/m.y. occurs where salt welds, overlying permeable carbonates thin to 500 m (1640 ft) or less. With tens of million years residence times feasible, convection under these conditions could locally result in reservoir sweet spots with porosity modification of 1% to 10% and potentially an order of magnitude or more in reservoir permeability. Integrating quantitative model–derived predictive diagenetic concepts with traditional subsurface data sets refines exploration to production scale risking of carbonate reservoir presence and quality.

Desktop /Portals/0/PackFlashItemImages/WebReady/geothermal-convection-in-south-atlantic-subsalt.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 5682 Bulletin Article

See Also: CD DVD

Desktop /Portals/0/images/_site/AAPG-newlogo-vertical-morepadding.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 4088 CD-DVD
Desktop /Portals/0/images/_site/AAPG-newlogo-vertical-morepadding.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 4576 CD-DVD

See Also: Online e Symposium

The goal of this e-symposium is to provide an overview of the latest trends and technologies for water management for oil and gas drilling, completions, and production.

Desktop /Portals/0/PackFlashItemImages/WebReady/oc-es-successful-oilfield-water-management.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 1496 Online e-Symposium