When Thin is In, Enhancement Helps

Expanding the frequency bandwidth of surface seismic data is an unending quest for geophysicists, because increased seismic resolution is essential for extracting stratigraphic detail from seismic images.

While both vertical resolution and horizontal resolution are important for interpreting small geologic features on seismic data, we focus our attention here on vertical resolution – recognizing that migration procedures usually enhance horizontal resolution.

If the frequency spectrum of a seismic wavelet is centered around 30 Hz, which is usually achievable, and the seismic interval velocity is greater than 3000 m/s, reservoirs having a thickness less than 25 meters may not be resolved. “Not resolved” means there is no distinct reflection peak or trough centered on the top and bottom interfaces of the reservoir unit.

This interval thickness, where seismic data can no longer position a distinct reflection peak or trough at the top and base of the interval, is called “tuning thickness.”

Because numerous stratigraphic targets have thicknesses of 10 meters or less – which is thinner than tuning thickness for most seismic profiles – frequency enhancement procedures need to be applied to seismic data to study reservoir targets in this “thinner than tuning thickness” domain.

One post-stack spectral inversion method that resolves thin layers having a thickness less than tuning thickness was published in 2005 (contact the authors for the reference). This method is driven by geological principles rather than by mathematical assumptions and uses spectral decomposition to enhance the frequency spectrum local to a thin-bed unit.

This spectral, or thin-bed, reflectivity inversion outputs a reflectivity series, and the apparent resolution of the inversion product is superior to the resolution of the input seismic data used to generate the reflectivity response.

Applications of this method in deconvolving complex seismic interference patterns are changing the mindset of many seismic interpreters, because the technique shows stratigraphic patterns with such remarkable detail.

The method consists of the following steps:

  1. A set of time-varying and space-varying wavelets is estimated from the seismic data.

For this purpose, it is good to have well control data to aid in selecting optimal space and time dependencies that should be expressed by these wavelets. In the absence of well control, a statistical method of wavelet estimation can be adopted.

  1. The wavelets estimated in step 1 are removed from the seismic data using an inversion procedure in which spectral constraints are derived on the basis of spectral decomposition procedures.

It is important to note that no Earth model or any assumption about stratigraphic layering is used in this inversion procedure – the trace-by-trace inversion procedure requires no starting geologic model and has no lateral continuity constraints.

Figure 1 shows a comparison of a segment of a seismic section from Alberta, Canada, before and after reflectivity inversion. After reflectivity inversion, more reflection detail can be seen, and faults are shown with improved clarity.

Once thin-bed reflectivity is derived from an input seismic volume – using, for example, a wavelet derived from an existing well – an interpreter can determine the amount of uncertainty involved in the inversion process by using a blind-well test.

Our experience with such exercises suggests that thin-bed spectral inversion creates data that tie favorably with other wells positioned in the same 3-D seismic volume.

Figure 2 shows a comparison between a segment of an input seismic section (figure 2a) and an equivalent segment of thin-bed reflectivity that has been convolved with a bandpass wavelet that extends the high end of the frequency spectrum to 120 Hz (figure 2b). Enhanced resolution of the reflectivity section is indicated by the extra reflection cycles.

More individual reflection cycles can now be tracked, leading to more detailed interpretation of the data.

Seismic attributes are a great help in extracting geologic information and are widely used to map geologic features at many scales. Geologic information not revealed by conventional displays of seismic data can often be seen on displays of one or more attributes derived from the data.

As a result, there has been an explosive growth in the development and application of seismic attributes.

Attribute computation done on data with enhanced resolution proves to be particularly useful for mapping onlap and offlap patterns or other stratigraphic features, which facilitates the mapping of parasequences and the direction of sediment transport.

Figure 3 shows a comparison of a stratal slice through a coherence-attribute volume generated for both input seismic data and for enhanced-resolution data.

Notice the significant impact that enhanced resolution has on the coherence attribute, as evidenced by the increased lateral resolution of the channel system and by the improved faulting picture seen in figure 3b.


The thin-bed spectral inversion method discussed here is a novel way of removing wavelet effects from seismic data to create a pure reflectivity sequence. For data with a high signal-to-noise ratio, units with thicknesses less than the tuning thickness of the input data can be resolved.

The improved-resolution seismic data retrieved in the form of reflectivity data are not only important for more accurate geologic interpretations but prove to be advantageous for:

  • Convolving the extracted reflectivity with a wider bandpass wavelet (say 5-120 Hz) to provide a high-frequency section.
  • Providing high-frequency attributes that enhance lateral resolution of geologic features.

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.

Yong Xu are with Arcis Corp., Calgary, Canada

John P. Castagna is with the University of Houston/Fusion Geo Inc., Houston.

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

This article describes a 250-m (820-ft)-thick upper Eocene deep-water clastic succession. This succession is divided into two reservoir zones: the lower sandstone zone (LSZ) and the upper sandstone zone, separated by a package of pelitic rocks with variable thickness on the order of tens of meters. The application of sequence-stratigraphic methodology allowed the subdivision of this stratigraphic section into third-order systems tracts.

The LSZ is characterized by blocky and fining-upward beds on well logs, and includes interbedded shale layers of as much as 10 m (33 ft) thick. This zone reaches a maximum thickness of 150 m (492 ft) and fills a trough at least 4 km (2 mi) wide, underlain by an erosional surface. The lower part of this zone consists of coarse- to medium-grained sandstones with good vertical pressure communication. We interpret this unit as vertically and laterally amalgamated channel-fill deposits of high-density turbidity flows accumulated during late forced regression. The sandstones in the upper part of this trough are dominantly medium to fine grained and display an overall fining-upward trend. We interpret them as laterally amalgamated channel-fill deposits of lower density turbidity flows, relative to the ones in the lower part of the LSZ, accumulated during lowstand to early transgression.

The pelitic rocks that separate the two sandstone zones display variable thickness, from 35 to more than 100 m (115–>328 ft), indistinct seismic facies, and no internal markers on well logs, and consist of muddy diamictites with contorted shale rip-up clasts. This section is interpreted as cohesive debris flows and/or mass-transported slumps accumulated during late transgression.

The upper sandstone zone displays a weakly defined blocky well-log signature, where the proportion of sand is higher than 80%, and a jagged well-log signature, where the sand proportion is lower than 60%. The high proportions of sand are associated with a channelized geometry that is well delineated on seismic amplitude maps. Several depositional elements are identified within this zone, including leveed channels, crevasse channels, and splays associated with turbidity flows. This package is interpreted as the product of increased terrigenous sediment supply during highstand normal regression.

Desktop /Portals/0/PackFlashItemImages/WebReady/sequence-stratigraphy-of-the-eocene-turbidite.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 5773 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 4502 CD-DVD

See Also: DL Abstract

Seismic correlations and well data confirm that deep-water carbonate beds of Mesozoic age have been found above the shallow allochthonous salt canopy in the northern Gulf of Mexico. These rafts of carbonate strata often overlie equivalent age Mesozoic carbonates in their correct stratigraphic position below the salt canopy.

The origin of keel structures is presently not well understood. As deformation occurs after shallow canopy emplacement, the keels are fairly recent developments geologically. Volumetrically few but intriguing observations suggest possible basement involvement in keel formation.

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

See Also: Learn! Blog

We are just a few months out from the 2014 AAPG Fall Education Conference (FEC). This years FEC will be bringing together new technologies, techniques, and lessons learned to offer four great days of geosciencetraining to enhance and advance your career.

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

See Also: Short Course

Black Belt Ethics examines the various tenants that define the martial artist’s code of honor. The course reviews each of these tenants and discusses how they can be applied in our personal and professional lives.

Desktop /Portals/0/PackFlashItemImages/WebReady/ace2015-sc05-black-belt-hero.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 14562 Short Course