Let’s Get Small: Benefits of Reduced Bin Size

Contributors: Kurt Marfurt

Ideally, seismic data should be acquired at high spatial and temporal sampling, so that the small subsurface features of interest can be clearly seen on the seismic display.

Such interpretation is easiest when the geological features of interest are uniformly illuminated, which in turn is possible by recording the scattered seismic wavefield on a uniform surface grid. The “nominal grid” is defined by the sourceto- source spacing within a shot line and the geophone-group-to-geophone-group spacing within a receiver line.

One also can increase the trace density by reducing the spacing between shot lines and between receiver lines.

Once acquired, data processing workflows are designed to retain the highest possible lateral and vertical resolution of the geologic target. Because of the different ray paths, dense acquisition (closer source and receiver lines) provides greater leverage against backscattered ground roll and interbed multiples, as well as decreased migration artifacts.

Uniform acquisition results in increased lateral continuity of amplitudes that otherwise may be contaminated by acquisition footprint.

In reality, economic, permitting and physical access constraints result in 3-D seismic data that are not uniformly sampled in all directions.

Vibrators require access along roads or open country, while geophones easily can be placed in forest or cultivated farmland. For this reason, “shot lines” may be more coarsely spaced than “receiver lines,” resulting in rectangular rather than square bins. Obstacles such as ponds, road, buildings and archaeological sites give rise to holes in an otherwise uniform acquisition grid.

Such non-uniformity in offsets and azimuths introduce additional artifacts in the final image. This non-uniformity also affects the performance of the processing algorithms – and so could lead to suboptimally processed data, affecting subsequent data interpretation.

In principal, any processing algorithm can be modified to handle sparse data. In practice, it is much easier to write an accurate 5-D interpolation algorithm, and thereby precondition the seismic data such that well-calibrated (and perhaps more complicated) algorithms such as prestack migration and prestack inversion work well.

While there is no substitute for acquiring good quality seismic data that has the above-mentioned qualities, it is possible to mimic or address some of the problems that crop up due to the sub-optimum parameterization used in the acquisition, during processing of the data.

In the June 2013 EXPLORER, we described the advantages of regularizing the offsets and azimuths of the input seismic data during processing by way of 5-D interpolation, which then aids the computation of seismic attributes.

Regularization of seismic data has proven to be a successful method – not only for superior imaging of post-stack data but for prestack analysis as well, such as AVO (amplitude versus offset) or AVAz (amplitude versus azimuth). The advantages accrue from the enhanced spatial and azimuthal sampling of the 3-D seismic data.

After doing the trace edits, amplitude recovery, refraction statics, preliminary velocity analysis and trim statics, one can run 5-D interpolation to “regularize” the data to have a uniform coverage of offsets and azimuths, thereby conditioning them for more detailed velocity analysis, noise rejection, prestack time migration and prestack impedance inversion.

In our June 2013 article we showed how such 5-D interpolation reduced artifacts in amplitude as well as in coherence and curvature attributes, but somewhat reduced the lateral resolution.

In this article we demonstrate the results of reducing the bin size of the seismic data as part of the 5-D interpolation process.

The source and receiver spacings are both 40 meters, giving rise to a nominal bin size of 20 meters by 20 meters. The primary use of 5-D interpolation is to fill in missing shots, receivers, offsets and azimuths corresponding to the nominal grid.

However, there is nothing preventing us in postulating a denser 10-meter by 10-meter grid and interpolating the corresponding unmigrated surface data. This smaller bin size is the interpolated equivalent of a (four times) more densely acquired survey.

We then use the same processing and prestack migration procedure applied to the interpolated data on the 20-meter by 20-meter grid.

The results we show are from a land seismic data volume from western Canada. We used 5-D interpolation to build missing traces for both 20 by 20 and 10 by 10 bin sizes.

In figure 1 we show a comparison of coherence horizon slices generated from data with 5-D interpolation at the nominal 20-meter by 20-meter bin size and at the “dense survey” 10-meter by 10-meter bin size.

Notice the enhanced resolution of the faults and the suppression of the NE-SW trending acquisition footprint.

A comparison of the seismic amplitude data before and after bin size reduction is shown in figure 2.

In the dashed box we notice somewhat clearer seismic signatures corresponding to the channel features that can be seen on the coherence phantom horizon slices comparison shown in figures 3.

Although the distributary channel system seen on the coherence slices is well imaged at the nominal grid size, the interpolated surface data provide much sharper individual channel limbs. Such enhanced quality imaging of data in terms of suitable seismic attributes helps squeeze out more information from the seismic data – and contributes in a generous way to the overall interpretation of the data as well.

We thank Arcis Seismic Solutions and 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

Reservoir properties of Upper Triassic–Middle Jurassic sandstones, Spitsbergen, are studied as part of a CO2 storage pilot project in Longyearbyen. The reservoir formations show large contrasts in sandstone compositions, with unexpected low permeability despite moderate porosity values. Petrographic analyses were performed to investigate the influence and distribution of diagenesis. It is concluded that, because of various compaction, cementation, and dissolution processes, the sandstone porosity is mainly isolated molds and micropores and associated with fibrous illite and chamosite, explaining the low permeability. Diagenesis and the distribution of quartz cement is influenced by lithofacies and detrital compositions. Mineralogically immature sandstones (De Geerdalen Formation) show a homogeneous distribution of quartz cement overgrowths on quartz grains, distributed interstitial to labile grains and other cements (e.g., late calcite). The main silica source was from the dissolution of adjacent feldspar and labile grains as part of the chemical compaction. In contrast, quartz-dominated sandstones (Knorringfjellet Formation) show a heterogeneous patchy distribution of quartz cement influenced by the sedimentary bioturbation pattern, with silica sourced also from dissolution at clay-rich microstylolites. Phosphatic beds at the base and top of the formation are strongly influenced by marine eogenesis and reworking processes and associated with concentration of iron-rich authigenic minerals. The highest porosity appears in sand-supported conglomerate where moldic clay-mineral ooids contributed to reduce quartz cementation. The stratigraphic change from mineralogical immature (Triassic) to mature (uppermost Triassic–Jurassic) sandstone compositions is detected in wide areas of the Barents Shelf and has considerable implications for the distribution of sandstone reservoir properties.
Desktop /Portals/0/PackFlashItemImages/WebReady/diagenesis-and-quartz-cement-distribution.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3713 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 4564 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 4560 CD-DVD

See Also: DL Abstract

Hydrothermal alteration occurs when relatively high-pressure, high-temperature fluids flow up active faults and into permeable formations that underlie sealing shales, or other low permeability strata. This process can and commonly does occur at relatively shallow burial depths of less than a kilometer and in many cases less than 500 meters.

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

See Also: Online e Symposium

The presentation will focus on hydraulic fracture geometry in shales, the materials used in the fracturing process, and treatment monitoring via microseismic.

Desktop /Portals/0/PackFlashItemImages/WebReady/oc-es-hydraulic-fracturing-of-shale-gas.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 1463 Online e-Symposium