Adding Texture Attributes to the 3-D Mix

In a previous Geophysical Corner (“A New Approach to Stratigraphic Interpretation,” September EXPLORER), we introduced a new set of seismic attributes that play an important role in extracting detailed stratigraphic information from seismic data.

The attributes in question were derived from a HorizonCube, an interpretation technique that provides fully interpreted seismic volumes where horizons are automatically tracked between a given set of framework horizons and faults.

This month, we go further and examine one other set of attributes – specifically, texture attributes, and how they can combine with HorizonCube attributes for 3-D segmentation.

Neural network-based waveform segmentation workflows have proved to be a highly valuable instrument for seismic interpreters to quickly visualize seismic patterns in a relatively thin interval of interest (see figure 1). These networks compare the input vector (the waveform) with a set of pre-calculated vectors that represent segment (cluster) centers.

The resulting segmentation maps show for each position the winning segment center – and such maps often reveal patterns that can be interpreted in terms of geological features.

Neural network-based segmentation also can be performed in 3-D. However, this cannot be done by feeding the neural network with waveforms, as is done in the horizon-guided approach. This is because waveforms change along the application window while segment centers are fixed.

The solution to this problem is to feed the neural work with phase-independent seismic attributes, such as energy, or amplitude spectra derivatives.

There are not that many uncorrelated attributes that can be used effectively for 3-D segmentation, which may explain why 3-D segmentation is not as successful as horizon-guided waveform segmentation.

There are two particular sets of attributes, however, that can combine together to address the challenges of 3-D segmentation:

  • Attributes derived from a HorizonCube and that were described in the previous article.
  • Texture attributes.

Texture attributes are popular in image processing and are used in seismic interpretation to carry out seismic facies analysis and highlight geo-morphological features. They originate from image analysis and are based on the Grey Level Co-occurrence Matrix (GLCM), which describes the relationship between pixels and was developed to capture the roughness or smoothness of an image.

A GLCM is a 2-D matrix of N x N dimensions representing the amplitude values of the reference pixel versus the amplitudes of the neighboring pixel. The matrix is filled by comparing each amplitude in the input area (volume) with its direct neighbor and increasing the occurrence of the corresponding matrix cell.

This is repeated for all amplitude pairs in the input cube, which then are converted into probabilities.

The GLCM thus captures how probable it is to find pairs of neighbouring amplitudes in the area (volume) around the evaluation point.

Texture attributes are computed in two steps:

  • First, the GLCM is computed for an area (volume) around the evaluation point.
  • Second, a statistical property from the GLCM is returned.

The GLCM input volume can be “dip-steered,” meaning that the input follows the stratigraphic layering, which results in sharper attribute responses for dipping strata.

Three groups of texture attributes are computed from the GLCM:

  • Contrast (contrast, dissimilarity and homogeneity), where measurements are calculated through the use of weights related to the distance from the GLCM diagonal.
  • Orderliness (ASM, energy and entropy), where interpreters measure how regular the pixel values are within the window.
  • Statistics (mean, variance and correlation), which are derived from the GLCM.

In each group, the attributes are highly correlated.

Through the use of both texture and the already described HorizonCube attributes, seismic interpreters now have two new groups of attributes that can be used as inputs to UVQ (Unsupervised Vector Quantization) networks for 3-D seismic facies analysis.

Figure 1 shows a horizon-guided UVQ waveform segmentation map that captured the seismic response below a maximum flooding surface in a wave-dominated Pliocene deltaic setting, offshore Netherlands. The UVQ network was trained to segment the interval of interest into 10 segments.

The UVQ segmentation map reveals several key geomorphological patterns that help to understand the depositional environment and the influence of salt tectonics.

The NW-SE oriented dark brown-red features on the right are sand ridges of 10-20 meters in height, developed parallel to the coast. These features are analogous to present day deposits observed along the Dutch coast.

Furthermore, NE-SW oriented deepwater channel systems are recognized (purple-red, on the left). These narrow channel-levee systems are developed as a result of halo kinetic movement of Zechstein salt in the northeast (upper right) corner of the image. Up-dip these channels cross-cut the sand ridges while down-dip they meander and bifurcate into the basin, where turbiditic deposits could be developed.

Figure 2 shows a seismic display for UVQ segmentation using both texture and HorizonCube attributes with energy, variance and correlation for the texture attributes; density and thickness for the HorizonCube attributes; and dip-steered energy as the input attribute for the neural network.

There are six classes in the illustration with it possible to see different behaviors in different areas.

If interpreters have carried out a systems tracts interpretation, they can use this interpretation to more or less force the segmentation to adhere to one’s own interpretation through a biased 3-D UVQ as seen in figure 3.

In addition to texture and energy attributes, figure 3 illustrates how a systems tract’s interpretation is used as bias for UVQ segmentation with, in this case, HorizonCube attributes not used as input.

The illustration shows the impact of using systems tracts as an input and the final segmentation. These classes (such as the yellow class) can then be extracted and seen in 3-D, allowing the interpreters to better define geobodies.

HorizonCube attributes are particularly useful for identification of stratigraphic features – such as pinch-outs, clino-forms, unconformities and condensed sections – whereas texture attributes play an important role in seismic facies analysis and in highlighting geo-morphological features.

When combined together, HorizonCube and texture attributes play an important role in 3-D segmentation and in identifying and mapping potential reservoir lithology in seismic data.

Comments (0)


Geophysical Corner

Geophysical Corner - Nanne Hemstra
Nanne Hemstra is with dGB Earth Sciences, Enschede, Netherlands.

Geophysical Corner

Geophysical Corner - Farrukh Qayyum
Farrukh Qayyum is with dGB Earth Sciences, Enschede, Netherlands.

Geophysical Corner

Geophysical Corner - Paul de Groot
Paul de Groot is with dGB Earth Sciences, Enschede, Netherlands.

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

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


Image Gallery

See Also: Book

Alternative Resources, Structure, Geochemistry and Basin Modeling, Sedimentology and Stratigraphy, Geophysics, Business and Economics, Engineering, Petrophysics and Well Logs, Environmental, Geomechanics and Fracture Analysis, Compressional Systems, Salt Tectonics, Tectonics (General), Extensional Systems, Fold and Thrust Belts, Structural Analysis (Other), Basin Modeling, Source Rock, Migration, Petroleum Systems, Thermal History, Oil Seeps, Oil and Gas Analysis, Maturation, Sequence Stratigraphy, Clastics, Carbonates, Evaporites, Seismic, Gravity, Magnetic, Direct Hydrocarbon Indicators, Resource Estimates, Reserve Estimation, Risk Analysis, Economics, Reservoir Characterization, Development and Operations, Production, Structural Traps, Oil Sands, Oil Shale, Shale Gas, Coalbed Methane, Deep Basin Gas, Diagenetic Traps, Fractured Carbonate Reservoirs, Stratigraphic Traps, Subsalt Traps, Tight Gas Sands, Gas Hydrates, Coal, Uranium (Nuclear), Geothermal, Renewable Energy, Eolian Sandstones, Sheet Sand Deposits, Estuarine Deposits, Fluvial Deltaic Systems, Deep Sea / Deepwater, Lacustrine Deposits, Marine, Regressive Deposits, Transgressive Deposits, Shelf Sand Deposits, Slope, High Stand Deposits, Incised Valley Deposits, Low Stand Deposits, Conventional Sandstones, Deepwater Turbidites, Dolostones, Carbonate Reefs, (Carbonate) Shelf Sand Deposits, Carbonate Platforms, Sebkha, Lacustrine Deposits, Salt, Conventional Drilling, Directional Drilling, Infill Drilling, Coring, Hydraulic Fracturing, Primary Recovery, Secondary Recovery, Water Flooding, Gas Injection, Tertiary Recovery, Chemical Flooding Processes, Thermal Recovery Processes, Miscible Recovery, Microbial Recovery, Drive Mechanisms, Depletion Drive, Water Drive, Ground Water, Hydrology, Reclamation, Remediation, Remote Sensing, Water Resources, Monitoring, Pollution, Natural Resources, Wind Energy, Solar Energy, Hydroelectric Energy, Bioenergy, Hydrogen Energy
Desktop /Portals/0/PackFlashItemImages/WebReady/book-s56-atlas-deep-water-outcrops-hero.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 7811 Book

See Also: Bulletin Article

We present a method of using fault displacement-distance profiles to distinguish fault-bend, shear fault-bend, and fault-propagation folds, and use these insights to guide balanced and retrodeformable interpretations of these structures. We first describe the displacement profiles associated with different end-member fault-related folding models, then provide examples of structures that are consistent with these model-based predictions. Natural examples are imaged in high-resolution two- and three dimensional seismic reflection data sets from the Niger Delta, Sichuan Basin, Sierras Pampeanas, and Cascadia to record variations in displacement with distance updip along faults (termed displacement-distance profiles). Fault-bend folds exhibit constant displacement along fault segments and changes in displacement associated with bends in faults, shear fault-bend folds demonstrate an increase in displacement through the shearing interval, and fault-propagation folds exhibit decreasing displacement toward the fault tip. More complex structures are then investigated using this method, demonstrating that displacement-distance profiles can be used to provide insight into structures that involve multiple fault-related folding processes or have changed kinematic behavior over time. These interpretations are supported by comparison with the kinematics inferred from the geometry of growth strata overlying these structures. Collectively, these analyses illustrate that the displacement-distance approach can provide valuable insights into the styles of fault-related folding.

Desktop /Portals/0/PackFlashItemImages/WebReady/fault-displacement-distance-relationships-as-indicators.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 5770 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 3929 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 4514 CD-DVD

See Also: Digital Download

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