Confident displays

Finding the Path Past 'Attribute Footprints'

Contributors: Kurt Marfurt

Three-D seismic surveys usually are designed in a way that the subsurface features are regularly sampled in different dimensions, comprising the spatial coordinates, offsets and azimuths.

Many processing algorithms require this regularity for their optimum performance. For example:

Marine acquisition suffers from platforms, shallow shoals, and tides and currents that give rise to feathering, all of which result in irregularity in spatial sampling of the data.

For older marine surveys, inlines are well sampled while crosslines are more coarsely spaced.

Land acquisition encounters a different suite of obstacles, such as highways, buildings and lakes.

Such obstacles, coupled with limited recording capacity and greater cost, results in missing data or “holes” in seismic data coverage.

  Finally, recording equipment malfunctions and noise bursts during acquisition may add more missing traces to the usable recordable data.

Sparse or missing data create problems while processing, as the different algorithms applied pre-stack or post-stack demand regularity in the offset and azimuth dimensions for optimum performance.

Non-uniformity in offsets and azimuths leads to inconsistencies in fold that follow a regular pattern we refer to as “acquisition footprint.” This imprint is an undesirable artifact that masks geologic features or amplitude variations seen on time slices from the seismic data, especially at shallow times.

Besides, the seismic data-derived attribute volumes also show acquisition footprint and other artifacts.

Obviously, the ideal way to fill in the missing data gaps would be to reshoot the data in those areas – although such infill acquisition would be extremely expensive per data point, if the equipment could be made available for such a small time in the field.

Such problems have been addressed at the processing stage since the advent of digital processing. The most common preconditioning of seismic data improves the signal-to-noise ratio of the seismic data by removing spatial noise or enhancing the coherency and alignment of the reflection events, without unnecessary smoothing or smearing of the discontinuities.

Although we usually think of removing unwanted features, we also can improve the signal-to-noise ratio by predicting unmeasured signal, such as dead traces and lower-fold areas corresponding to unrecorded offsets and azimuths in the gathers.

Prediction or population of missing traces in seismic data is referred to as interpolation. Initial interpolation methods used localized information, such as creating a missing trace by taking the lateral average of measured adjacent traces.

These methods worked well on stacked data volumes. A geologic feature seen on seismic data has three dimensions – time, an easting (x) and a northing (y). However, modern seismic processes such as prestack impedance inversion uses a fourth dimension of source-receiver offset (h), while azimuthal anisotropy uses a fifth dimension of source-receiver azimuth.

The gap between the measured traces of a given offset-azimuth volume can be quite large, such that simple interpolation of spatial neighbors no longer works.

During the last decade or so, more sophisticated methods for data interpolation have evolved that interpolate the missing traces using not only in with neighboring samples in t, x and y, but also in offset and azimuth. Such “5-D” interpolators operate simultaneously in all dimensions, and are able to predict the missing data with more accurate amplitude and phase variations.

As expected, these methods are compute intensive and have longer run-times than the simplistic interpolation methods.

We demonstrate here the application of one such method of 5-D interpolation on seismic data and show how it aids some of the seismic attributes derived from them.

In figure 1a we show a representative vertical slice through a merged 3-D amplitude volume that has many dead traces. Such dead traces are seen on other inlines as well.

The location of this inline is shown in figure 2a, where we show a horizon slice through the corresponding coherence volume.

The dead traces result in the speckled pattern indicated with yellow ellipses. To regularize the data, 5-D interpolation was run on the seismic data prior to migration with the equivalent displays shown in figures 1b and 2b, respectively.

Notice in figure 1b that not only are the missing traces interpolated, but the overall signal-to-noise ratio and reflector continuity is improved.

Similarly, note the absence of the speckles associated with the missing traces and the greater continuity of the channel and other discontinuity features as indicated by the red arrows.

The inference we draw from this example is that regularization by 5-D interpretation yields better-focused images. Interpretation carried out on such attributes will definitely be more accurate than the one carried out on data without regularization.

As 5-D interpolation discussed above regularizes the geometry of the seismic data, it addresses the root cause of the footprint arising due to the acquisition irregularities.

In figure 3a we show time slices at 158 ms, where the acquisition footprint appears prominently on the coherence attribute as striations in the NE-SW direction masking the reflection detail behind them. Figure 3b shows the equivalent coherence slice after 5-D regularization exhibiting considerable improvement in data quality.

Seismic data usually have geometry regularization issues that give rise to artifacts on geometric attribute displays. Five-D interpolation methods adopted during processing help address issues such as missing data pockets and acquisition footprint striations.

Coherence and curvature attributes computed on regularized seismic data yield displays clear of these artifacts, and so lead to more confident displays.

We thank Arcis Seismic Solutions, TGS, for encouraging this work and for permission to present these results.

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.


Image Gallery

See Also: Bulletin Article

The Upper Jurassic Arab Formation in the Arabian Peninsula, the most prolific oil-bearing interval of the world, is a succession of interbedded thick carbonates and evaporites that are defined stratigraphically upsection as the Arab-D, Arab-C, Arab-B, and Arab-A. The Arab-D reservoir is the main reservoir in Khurais field, one of the largest onshore oil fields of the Kingdom of Saudi Arabia.

In Khurais field, the Arab-D reservoir is composed of the overlying evaporitic Arab-D Member of the Arab Formation and the underlying upper part of the Jubaila Formation. It contains 11 lithofacies, listed from deepest to shallowest: (1) hardground-capped skeletal wackestone and lime mudstone; (2) intraclast floatstone and rudstone; (3) pelletal wackestone and packstone; (4) stromatoporoid wackestone, packstone, and floatstone; (5) Cladocoropsis wackestone, packstone, and floatstone; (6) Clypeina and Thaumatoporella wackestone and packstone; (7) peloidal packstone and grainstone; (8) ooid grainstone; (9) crypt-microbial laminites; (10) evaporites; and (11) stratigraphically reoccurring dolomite.

The Arab-D reservoir lithofacies succession represents shallowing-upward deposition, which, from deepest to shallowest, reflects the following depositional environments: offshore submarine turbidity fans (lithofacies 1 and 2); lower shoreface settings (lithofacies 3); stromatoporoid reef (lithofacies 4); lagoon (lithofacies 5 and 6); shallow subtidal settings (lithofacies 7 and 8); peritidal settings (lithofacies 9); and sabkhas and salinas (lithofacies 10). The depositional succession of the reservoir represents a prograding, shallow-marine, reef-rimmed carbonate shelf that was subjected to common storm abrasion, which triggered turbidites.

Desktop /Portals/0/PackFlashItemImages/WebReady/arabian-carbonate-reservoirs-a-depositional-model.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3771 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 4578 CD-DVD

See Also: Energy Policy Blog

China plans to significantly increase its natural gas consumption to help cut its appalling air pollution. But natural gas is still a small part of its energy mix. In addition, and to confound environmentalists, a significant part of China’s gas supply comes from Coal-to-Gas technology, which generates large volumes of greenhouse gas and other pollutants, but does allow China to deliver clean-burning gas to locations with severe air pollution.

Desktop /Portals/0/PackFlashItemImages/WebReady/waiting-for-chinas-natural-gas-revolution-2014-08aug-01-hero.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 12908 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

See Also: Short Course

This course is oriented towards the recognition and characterization of uncertainty in unconventional reservoirs. Starting with resource/reservoir assessment methods, it moves through the full unconventional value-chain. This two-day exercise and example filled workshop provides participants with the techniques and reasoning needed to validly assess the merits of the search for, and development of, unconventional resource plays.
Desktop /Portals/0/PackFlashItemImages/WebReady/hero-assessment-forecasting-and-decision-making-in-unconventional-resource-plays.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 13769 Short Course