Dense sets of seismic horizons

A New Approach to Stratigraphic Interpretation

Contributors: Paul DeGroot, Nanne Hemstra, Farrukh Qayyum

This article introduces a new set of seismic attributes that play an important role in extracting detailed stratigraphic information from seismic data.

The attributes in question are derived from a HorizonCube, one of the emerging interpretation techniques that provide fully interpreted seismic volumes (figure 1).

In this particular technique, horizons are automatically tracked between a given set of framework horizons and faults. The tracking is done using a seismic dip volume.

Compared with conventional amplitude tracking, this algorithm is more robust in areas with low signal-to-noise ratio, where diachronous events can be tracked as well as events that are phase inconsistent.

As these horizons are guided by a continuous dip-field, they may converge and diverge according to the dip of a seismic reflector. In this way, the key geologic features such as unconformities, pinchouts and condensed sections are highlighted.

Horizons can be tracked in two different modes (figure 1):

As continuous horizons that stay together when they converge and never cross each other.

Such horizons help identify unconformities and condensed sections – particularly useful in 3-D attributes visualization and reservoir modelling.

As truncated horizons that may stop when they meet each other in space based on a user-defined threshold.

This helps to identify stratigraphic lapouts (onlaps, downlaps and top laps). These also are useful in Wheeler diagram creation and sequence stratigraphic interpretation.

New Family of Stratigraphic Attributes

For the continuous set of seismic horizons, a new family of attributes can be computed that visualize geologic features previously hidden.

Attributes include:

Isochron thickness – This attribute highlights not only sedimentary bodies but also picks up local pinch-outs, condensed intervals and local unconformities.

Curvature – This curvature attribute tends to be smoother than the conventional volume curvature attribute, computed using the seismic dip volume without mapping.

HorizonCube density – Events can be counted within user-defined time windows, with high-density values corresponding to horizon convergence and low density values corresponding to the horizon divergence.

The “HorizonCube density attribute” helps define the zones of pinch-outs, condensed sections and unconformities.

Arbitrary layers –
A dense set of seismic horizons can be separated into a set of arbitrary (stratigraphically uninterpreted) layers that divide the mapped seismic horizons into fix layers with a unique ID for each layer.

This plays an important role in quality control where visualizing hundreds of horizons can often be slow and laborious.

Derivatives – Finally, first or second derivatives on a set of horizon events can be extended to perform a derivative on an isochron between a given ranges of events. These attributes measure subtle geometrical changes and discontinuities.

From a stratigraphic interpretation standpoint, the mapped seismic horizons in truncated form can be sub-divided into sequence stratigraphic units through the co-visualization of a structural domain and a Wheeler domain in conjunction with well data.

In both domains various data, such as wells and seismic, can be combined to interpret a set of sequence stratigraphic units, such as systems tracts and sequences.

Typical attributes that are extracted using this sequence stratigraphic framework include stratigraphic unit IDs (identification number) and isochrons. The IDs of interpreted stratigraphic units can either be unique (each systems tract unit is assigned a unique identification number) or common (each identical systems tract is assigned a common number).

An isochron attribute calculated the thickness of a stratigraphic unit. The unit of this attribute depends on the seismic survey type (TWT or in depth). This is a key attribute to understand how sedimentation filled a sedimentary basin as a function of geologic time.

The relative rate of preservation – the ratio between an isochron volume and a known geologic time-span for a particular unit – also can be calculated with the results being relative measurements of the rate of preservation per geological time unit.

Finally, an attribute can be generated that defines the difference between two isochron grids as a volume at a certain trace location.

Extracting Deepwater Channel Stages

During the Neogene period, Carnarvon Basin (Australian North West Shelf) comprised several deepwater meandering channelized systems. One of those systems is studied using the new stratigraphic unit ID attribute.

To our knowledge, this is the first study in which this attribute was used to define an unconventional set of stratigraphic units.

Convention nomenclatures are not adopted, as the interval may fall in one systems tract. Therefore, a further sub-division of channel system requires a new sequence model.

In this case, a “channel stages” sequence model is defined. It contains a set of N stratigraphic units defined as channel stages. Each unit is defined by its own unique ID.

Upper and lower boundaries of the channel system were manually mapped to form an initial framework to process the densely mapped set of seismic horizons. Within such a framework, a data-driven collection of densely mapped seismic horizons was created.

These were then sub-divided into a channel stages model that contained a user-defined set of color-coded stages. The criteria of sub-dividing a channel system into its stages are based on observations of features such as crosscutting relationships, timing and geomorphologic patterns.

Five channel stages can be seen in figure 2. The same interval also can be studied using the density attribute (figure 3), which clearly shows channel base, levees and banks of channel systems.

Workflows such as these support building better stratigraphic models to predict reservoir quality sands and to perform reservoir characterization.

Dense sets of seismic horizons help in defining new stratigraphic attributes – and these attributes help in understanding stratigraphy and geomorphology, and can be used as a means of performing stratigraphic interpretation on seismic data.

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

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


Image Gallery

Meet the Authors

Farrukh Qayyum
Farrukh Qayyum
Nanne Hemstra
Nanne Hemstra
Paul de Groot
Paul de Groot

AAPG members Farrukh Qayyum, Nanne Hemstra and Paul de Groot all are with dGB Earth Sciences, Enschede, Netherlands

See Also: Book

Desktop /Portals/0/PackFlashItemImages/WebReady/book-m95-Lacustrine-Sandstone-Reservoirs-and-Hydrocarbon-Systems.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3991 Book

See Also: Bulletin Article

The Marcellus Shale is considered to be the largest unconventional shale-gas resource in the United States. Two critical factors for unconventional shale reservoirs are the response of a unit to hydraulic fracture stimulation and gas content. The fracture attributes reflect the geomechanical properties of the rocks, which are partly related to rock mineralogy. The natural gas content of a shale reservoir rock is strongly linked to organic matter content, measured by total organic carbon (TOC). A mudstone lithofacies is a vertically and laterally continuous zone with similar mineral composition, rock geomechanical properties, and TOC content. Core, log, and seismic data were used to build a three-dimensional (3-D) mudrock lithofacies model from core to wells and, finally, to regional scale. An artificial neural network was used for lithofacies prediction. Eight petrophysical parameters derived from conventional logs were determined as critical inputs. Advanced logs, such as pulsed neutron spectroscopy, with log-determined mineral composition and TOC data were used to improve and confirm the quantitative relationship between conventional logs and lithofacies. Sequential indicator simulation performed well for 3-D modeling of Marcellus Shale lithofacies. The interplay of dilution by terrigenous detritus, organic matter productivity, and organic matter preservation and decomposition affected the distribution of Marcellus Shale lithofacies distribution, which may be attributed to water depth and the distance to shoreline. The trend of normalized average gas production rate from horizontal wells supported our approach to modeling Marcellus Shale lithofacies. The proposed 3-D modeling approach may be helpful for optimizing the design of horizontal well trajectories and hydraulic fracture stimulation strategies.

Desktop /Portals/0/PackFlashItemImages/WebReady/organic-rich-marcellus-shale-lithofacies-modeling.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 5725 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 3932 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 4561 CD-DVD

See Also: Energy Policy Blog

The U.S. Energy Information Administration reports that Europe consumed 18.7 trillion cubic feet (Tcf) of natural gas in 2013; Russia supplied 30% (5.7 Tcf) of this volume, and 16% (3.0 Tcf) of the total natural gas consumed in Europe passed through Ukraine’s pipeline network; Slovakia gets 80 percent of its gas, and Bulgaria gets 90 percent of its gas from Russia via pipeline through Ukraine.

Desktop /Portals/0/PackFlashItemImages/WebReady/the-ukraine-crisis-and-european-natural-gas-supplies-2014-07jul-09-hero.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 12905 Energy Policy Blog