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

Understanding the distribution and geometry of reservoir geobodies is crucial for net-to-gross estimates and to model subsurface flow. This article focuses on the process of dolomitization and resulting geometry of diagenetic geobodies in an outcrop of Jurassic host rocks from northern Oman. Field and petrographic data show that a first phase of stratabound dolomite is crosscut by a second phase of fault-related dolomite. The stratabound dolomite geobodies are laterally continuous for at least several hundreds of meters (sim1000 ft) and probably regionally and are one-half meter (1.6 ft) thick. Based on petrography and geochemistry, a process of seepage reflux of mesosaline or hypersaline fluids during the early stages of burial diagenesis is proposed for the formation of the stratabound dolomite. In contrast, the fault-related dolomite geobodies are trending along a fault that can be followed for at least 100 m (328 ft) and vary in width from a few tens of centimeters to as much as 10 m (sim1–33 ft). Petrography, geochemistry, and high homogenization temperature of fluid inclusions all point to the formation of the dolomite along a normal fault under deep burial conditions during the Middle to Late Cretaceous. The high 87Sr/86Sr ratio in the dolomite and the high salinity measured in fluid inclusions indicate that the dolomitizing fluids are deep basinal brines that interacted with crystalline basement. The dolomitization styles have an impact on the dimension, texture, and geochemistry of the different dolomite geobodies, and a modified classification scheme (compared to the one from Jung and Aigner, 2012) is proposed to incorporate diagenetic geobodies in future reservoir modeling.
Desktop /Portals/0/PackFlashItemImages/WebReady/linking-process-dimension-texture-and-geochemistry.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3772 Bulletin Article

Umiat field in northern Alaska is a shallow, light-oil accumulation with an estimated original oil in place of more than 1.5 billion bbl and 99 bcf associated gas. The field, discovered in 1946, was never considered viable because it is shallow, in permafrost, and far from any infrastructure. Modern drilling and production techniques now make Umiat a more attractive target if the behavior of a rock, ice, and light oil system at low pressure can be understood and simulated.

The Umiat reservoir consists of shoreface and deltaic sandstones of the Cretaceous Nanushuk Formation deformed by a thrust-related anticline. Depositional environment imparts a strong vertical and horizontal permeability anisotropy to the reservoir that may be further complicated by diagenesis and open natural fractures.

Experimental and theoretical studies indicate that there is a significant reduction in the relative permeability of oil in the presence of ice, with a maximum reduction when connate water is fresh and less reduction when water is saline. A representative Umiat oil sample was reconstituted by comparing the composition of a severely weathered Umiat fluid to a theoretical Umiat fluid composition derived using the Pedersen method. This sample was then used to determine fluid properties at reservoir conditions such as bubble point pressure, viscosity, and density.

These geologic and engineering data were integrated into a simulation model that indicate recoveries of 12%–15% can be achieved over a 50-yr production period using cold gas injection from five well pads with a wagon-wheel configuration of multilateral wells.

Desktop /Portals/0/PackFlashItemImages/WebReady/Integrated-reservoir-characterization-and-simulation-of.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 7968 Bulletin Article
A new hierarchical architectural classification for clastic marginal-marine depositional systems is presented and illustrated with examples. In ancient rocks, the architectural scheme effectively integrates the scales of sedimentology (core, outcrop) and sequence stratigraphy (wireline-log correlation, reflection seismic). The classification also applies to modern sediments, which allows for direct comparison of architectural units between modern and ancient settings. In marginal-marine systems, the parasequence typically defines reservoir flow units. This classification addresses subparasequence scales of stratigraphy that commonly control fluid flow in these reservoirs. The scheme consists of seven types of architectural units that are placed on five architectural hierarchy levels: hierarchy level I: element (E) and element set (ES); hierarchy level II: element complex (EC) and element complex set (ECS); hierarchy level III: element complex assemblage (ECA); hierarchy level IV: element complex assemblage set (ECAS); and hierarchy level V: transgressive-regressive sequence (T-R sequence). Architectural units in levels I to III are further classified relative to dominant depositional processes (wave, tide, and fluvial) acting at the time of deposition. All architectural units are three-dimensional and can also be expressed in terms of plan-view and cross-sectional geometries. Architectural units can be linked using tree data structures by a set of familial relationships (parent-child, siblings, and cousins), which provides a novel mechanism for managing uncertainty in marginal-marine systems. Using a hierarchical scheme permits classification of different data types at the most appropriate architectural scale. The use of the classification is illustrated in ancient settings by an outcrop and subsurface example from the Campanian Bearpaw–Horseshoe Canyon Formations transition, Alberta, Canada, and in modern settings, by the Mitchell River Delta, northern Australia. The case studies illustrate how the new classification can be used across both modern and ancient systems, in complicated, mixed-process depositional environments.
Desktop /Portals/0/PackFlashItemImages/WebReady/a-hierarchical-approach-to-architectural-classification.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3769 Bulletin Article

See Also: DL Abstract

Production from the Marcellus gas shale generated international interest when methane accumulated in the surface housing of a water well pump and exploded.

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

Hydrocarbon exploration beneath the shallow allochthonous salt canopy of the ultra-deepwater central Gulf of Mexico has encountered three thick, sand-rich, submarine fan successions that punctuate an otherwise relatively condensed and fine-grained basin center stratigraphy. These sand-rich fans are Late Paleocene, Early Miocene, and Middle Miocene in age and each coincide with periods of very high sediment flux and basin margin instability. They are the primary exploration targets in most ultra-deepwater fields, recent discoveries, and failed exploration tests.

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