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, 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.


See Also: Bulletin Article

Thirty-seven mudstone samples were collected from the uppermost Lower Mudstone Member of the Potrerillos Formation in El Gordo minibasin within La Popa Basin, Mexico. The unit is exposed in a circular pattern at the earth's surface and is intersected by El Gordo diapir in the northeast part of the minibasin. Vitrinite reflectance (Ro) results show that samples along the eastern side of the minibasin (i.e., south of the diapir) are mostly thermally immature to low maturity (Ro ranges from 0.53% to 0.64%). Vitrinite values along the southern, western, and northwestern part of the minibasin range between 0.67% and 0.85%. Values of Ro immediately northwest of the diapir are the highest, reaching a maximum of 1.44%. The results are consistent with two different possibilities: (1) that the diapir plunges to the northwest, or (2) that a focused high-temperature heat flow existed along just the northwest margin of the diapir. If the plunging diapir interpretation is correct, then the thermally immature area south of the diapir was in a subsalt position, and the high-maturity area northwest of the diapir was in a suprasalt position prior to Tertiary uplift and erosion. If a presumed salt source at depth to the northwest of El Gordo also fed El Papalote diapir, which is located just to the north of El Gordo diapir, then the tabular halokinetic sequences that are found only along the east side of El Papalote may be subsalt features. However, if the diapir is subvertical and the high-maturity values northwest of the diapir are caused by prolonged, high-temperature fluid flow along just the northwestern margin of the diapir, then both of these scenarios are in disagreement with previously published numerical models. This disagreement arises because the models predict that thermal anomalies will extend outward from a diapir a distance roughly 1.5 times the radius of the diapir, but the results reported here show that the anomalous values on one side of the diapir are about two times the radius, whereas they are as much as five times the radius on the other side of the diapir. The results indicate that strata adjacent to salt margins may experience significantly different heat histories adjacent to different margins of diapirs that result in strikingly different diagenetic histories, even at the same depth.
Desktop /Portals/0/PackFlashItemImages/WebReady/a-surprising-asymmetric-paleothermal-anomaly-around.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 5769 Bulletin Article

See Also: DL Abstract

The Gulf of Mexico (GOM) is the 9th largest body of water on earth, covering an area of approximately 1.6 million km2 with water depths reaching 4,400 m (14,300’). The basin formed as a result of crustal extension during the early Mesozoic breakup of Pangaea. Rifting occurred from the Late Triassic to early Middle Jurassic. Continued extension through the Middle Jurassic combined with counter-clockwise rotation of crustal blocks away from North America produced highly extended continental crust in the subsiding basin center. Subsidence eventually allowed oceanic water to enter from the west leading to thick, widespread, evaporite deposition. Seafloor spreading initiated in the Late Jurassic eventually splitting the evaporite deposits into northern (USA) and southern (Mexican) basins. Recent work suggests that this may have been accomplished by asymmetric extension, crustal delamination, and exposure of the lower crust or upper mantle rather than true sea floor spreading (or it could be some combination of the two).
Desktop /Portals/0/images/_site/AAPG-newlogo-vertical-morepadding.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 844 DL Abstract

See Also: Energy Policy Blog

Comment period is open for the EPA proposed rule that would require CCS, carbon capture and storage. The Environmental Protection Agency (EPA) proposed rule, "Standards of Performance for Greenhouse Gas Emissions for New Stationary Sources: Electric Utility Generating Units," is open for comment until March 10.

Desktop /Portals/0/PackFlashItemImages/WebReady/epa-power-Plant-regs-depend-on-enhanced-oil-recovery-2014-02feb-17.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 12375 Energy Policy Blog

The Center for Strategic and International Studies (CSIS) recently published a series of reports assessing how rapid growth in U.S., and possibly global, oil and gas production from shales may impact various net-energy exporting or importing countries.

Desktop /Portals/0/PackFlashItemImages/WebReady/Hero-Global-Impacts-of-US-Shale-Production.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 11996 Energy Policy Blog

See Also: 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