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.


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.


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.


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, chief geophysicist (reservoir), at Arcis Seismic Solutions, Calgary, Canada, 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.


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 4067 Book

See Also: Bulletin Article

Understanding the distribution and geometry of reservoir geobodies is crucial for net-to-gross estimates and to model subsurface flow. This article focuses on the process of dolomitization and resulting geometry of diagenetic geobodies in an outcrop of Jurassic host rocks from northern Oman. Field and petrographic data show that a first phase of stratabound dolomite is crosscut by a second phase of fault-related dolomite. The stratabound dolomite geobodies are laterally continuous for at least several hundreds of meters (sim1000 ft) and probably regionally and are one-half meter (1.6 ft) thick. Based on petrography and geochemistry, a process of seepage reflux of mesosaline or hypersaline fluids during the early stages of burial diagenesis is proposed for the formation of the stratabound dolomite. In contrast, the fault-related dolomite geobodies are trending along a fault that can be followed for at least 100 m (328 ft) and vary in width from a few tens of centimeters to as much as 10 m (sim1–33 ft). Petrography, geochemistry, and high homogenization temperature of fluid inclusions all point to the formation of the dolomite along a normal fault under deep burial conditions during the Middle to Late Cretaceous. The high 87Sr/86Sr ratio in the dolomite and the high salinity measured in fluid inclusions indicate that the dolomitizing fluids are deep basinal brines that interacted with crystalline basement. The dolomitization styles have an impact on the dimension, texture, and geochemistry of the different dolomite geobodies, and a modified classification scheme (compared to the one from Jung and Aigner, 2012) is proposed to incorporate diagenetic geobodies in future reservoir modeling.
Desktop /Portals/0/PackFlashItemImages/WebReady/linking-process-dimension-texture-and-geochemistry.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3772 Bulletin Article

See Also: Learn! Blog

AAPG offers two short courses in conjunction with this year’s Unconventional Resources Technology Conference (URTeC). A wealth of information in a short period of time, theses short courses are an effective and efficient way to learn about the industry.

Desktop /Portals/0/PackFlashItemImages/WebReady/AAPG-learn-blog-200x200.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 11160 Learn! Blog

Looking at a formation as a source rock, then turning around and considering it a viable reservoir requires you to be able to shift your thinking and to analyze a great deal of data in a new way. If you don’t, you risk not understanding the nature of “sweet spots” and how to accurately complete or use reservoir characterization studies.

Desktop /Portals/0/PackFlashItemImages/WebReady/ec-fec-my-source-rock-is-now-my-reservoir.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 12765 Learn! Blog

See Also: Short Course

This course is a non-numerical introduction to the use of geochemistry and BPSM to better understanding unconventional resources. This course is designed to provide participants with new information on unconventional and sweet spot identification that is not normally available in routine service company courses.

Desktop /Portals/0/PackFlashItemImages/WebReady/hero-geochemistry-and-modeling-unconventional-petroleum-systems.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 20059 Short Course