Geophysical Corner

Extracting Information From Texture Attributes

Figure 1 – Comparison of strat slices through (a) input seismic amplitude, (b) coherence, (c) GLCM energy, (d) GLCM entropy and (e) GLCM homogeneity attributes. On 1(f), a chair display with the vertical seismic amplitude correlated with a horizontal strat-slice through the most-positive curvature (long-wavelength) volume. The channel features are seen better defined on the texture attribute displays than the seismic display.
Figure 1 – Comparison of strat slices through (a) input seismic amplitude, (b) coherence, (c) GLCM energy, (d) GLCM entropy and (e) GLCM homogeneity attributes. On 1(f), a chair display with the vertical seismic amplitude correlated with a horizontal strat-slice through the most-positive curvature (long-wavelength) volume. The channel features are seen better defined on the texture attribute displays than the seismic display.

There are a number of seismic attributes that are derived from seismic amplitudes to facilitate the interpretation of geologic structure, stratigraphy and rock/pore fluid properties.

The earliest attributes were extracted by treating seismic amplitudes as analytic signals for aiding feature identification and interpretation.

As the computation of these attributes is carried out at each sample of the seismic trace, they are referred to as instantaneous attributes.

This development was followed by attributes that are derived by transforming seismic amplitudes into impedance or velocity.

Also called seismic impedance inversion attributes, these attributes yield lithology or fluid information that can be calibrated with well logs.

A third class of attributes quantifies the lateral changes in waveform using an ensemble of windowed traces in the inline and crossline directions.

Such geometric attributes include dip, coherence and curvature, and are routinely used to accelerate and quantify the interpretation of faults, fractures and folds from 3-D seismic data.

While texture attributes are less familiar to seismic interpreters, seismic texture forms the basis of seismic stratigraphy, giving rise to descriptions of “concordant,” “blocky,” “hummocky” and “chaotic” pictures.

Quantitative texture analysis is one of the primary tools in remote sensing of forestry, agriculture and urban planning.

The classic definition of texture defines a window, such as the human thumb, sampling subtle changes in elevation. Rubbing your thumb across nearby surfaces may give rise to textures you may describe as smooth, rough, silky, corrugated, wavy or chaotic.

Most people can easily recognize pine, oak, maple, ash, mahogany, teak and many other woods from their grain, but may have difficulty explaining how they are able to distinguish them.

For this reason, it is difficult to teach a computer to recognize such patterns.

Most remote sensing and industrial applications use statistical measures of the gray-level co-occurrence matrix, or GLCM, which measures the repetition of a pattern from point-to-point.

Thus a “brick pattern” in North America would have mortar every 12 inches horizontally and four inches vertically.

GLCM seismic analysis might search for vertical patterns such as onlap, frequency and parallelism.

In this article we search for lateral patterns in the seismic data along structural dip.

We find three texture attributes to be the most useful in extracting lateral changes in reflectivity. They are:

  • GLCM energy.
  • GLCM entropy.
  • GLCM homogeneity.

Somewhat confusingly, the GLCM energy is a measure of the energy of the GLCM matrix and not of the seismic data itself.

For this reason, a checkerboard pattern, which has many adjacent red and black pixels, will have high GLCM energy, high homogeneity and low entropy. A smooth pattern will have high homogeneity, moderate energy and low entropy.

We illustrate the application of these texture attributes and their usefulness on an area in south central Alberta, Canada.

In figure 1a we see a strat slice through a seismic volume showing some Mannville channels. Not all these channels are incised, as the main channel on the left (blue arrow) is seen to have a signature somewhat different from the channel seen to the right and indicated with a green arrow. This is because of the greater measure of differential compaction noticed on the curvature strat slice (shown in figure 1f), and described in the July 2012 Geophysical Corner.

This main channel is seen to have a definite outline in blue on the seismic display, and at the location of the pink arrow it merges with the vertical channel to the right (green arrow), which appears to have undergone lesser differential compaction.

A thin vertical channel seen on the seismic amplitude display in figure 1a (yellow arrow) is seen with a better definition on the coherence.

While coherence shows the edges of the channel, it gives little indication of the heterogeneity or uniformity of the channel fill. Notice the clear definition of this channel on the three texture attributes shown in figures 1c-e, especially the complete thin high entropy, low homogeneity N-S running channel seen in figures 1d and e.

We interpret a similar high entropy, low homogeneity feature in figures 1d and e to be a point bar in the middle of the incised valley (green arrows). This internal architecture was not delineated by coherence.


Unlike geometric attributes, which are clearly linked to faults, folds and fractures, texture attributes are more difficult to interpret.

In remote sensing of forestry and agriculture, calibration is obtained by control sites, with a human being visiting a given location and literally providing ground truth.

In seismic texture the ground truth is provided by well control, interpreter experience and an understanding of geologic processes.

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.


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 4461 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 4046 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 8146 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: Short Course

This course is designed for geologists who wish to interpret and explore lacustrine microbialites and associated carbonate facies for hydrocarbons, or just want to familiarize themselves with microbialites and lacustrine systems.

Desktop /Portals/0/PackFlashItemImages/WebReady/ace2015-sc14-Microbialites-hero1.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 14608 Short Course