Seismic Aids Heterogeneity Hunt

Contributors: Satinder Chopra, Yong Xu

Canada’s Athabasca oil sands represent the biggest petroleum accumulation in the world and presently produce more than one million barrels of oil per day. These are the Lower Cretaceous bitumen sand reservoirs comprising the McMurray Formation that varies in thickness from 10 to 90 meters and occurs at depth of 0 to 400 meters.

While the shallow oil sands are exploited by open-pit mining, the deeper reservoirs are produced through some type of insitu bitumen production like SAGD (Steam-Assisted Gravity Drainage). SAGD operations require extensive, bitumen-saturated, homogeneous formations for optimum production.

The McMurray formation, however, is heterogeneous in terms of reservoir continuity, mineralogy, sedimentary facies and water-saturation, and is too complex to be understood from the sparse available core database.

Surface seismic data are one option for characterizing this reservoir heterogeneity, with a common approach being to use neural networks or statistical analysis at well locations to deduce relationships between seismic attributes and lithology. These relationships are then used to determine lateral lithology variations between wells.

We describe here a two-step approach to understand the heterogeneity of Athabasca oil-sand reservoirs.

  • The first step involves a rock physics study to understand relationships between lithology and petrophysical parameters. From this effort, lithology-sensitive rock parameters are selected that can be detected seismically.
  • The second step is to derive these lithology parameters from seismic data.

The first step – rock physics analysis – is carried out for various rock physics parameters across the zone of interest. Parameters that exhibit the best sensitivity to lithology are selected. For example, on figure 1 (below) we show P-impedance, Vp/Vs velocity ratio and density cross-plotted against gamma ray using log samples from the study area.

For these study wells, P-impedance shows a limited ability to distinguish lithology; Vp/Vs ratio shows a gentle variation with lithology; and density appears to be the best indicator of clay content.

Now that the desired rock parameter – density – has been determined, step two is to do AVO analyses of prestack seismic data to estimate rock density along seismic profiles.

Normally, density determination is done using a three-term AVO analysis that requires prestack seismic data with long offsets. In our study we improved this conventional approach by adopting innovative ideas like using a windowed approach instead of a sample-by-sample computation for deriving AVO attributes, reducing distortion due to NMO stretch and offset-dependent tuning, using error-based weights and accounting for the strong reflections from the McMurray Devonian interface.

On figure 2 we show a density section derived from seismic data. The lateral variation in seismic-based density shows the richest sand areas (in green color) within the mid-McMurray are around wells 5 and 6, with good shaley cap rocks in the upper-McMurray.

These predictions are verified by gamma ray logs acquired in both wells.

Recently drilled wells 3 and 7 served as blind well tests. Well 3 found mainly a shaley facies within the McMurray, and the seismic inversion for density agrees with these log calibration data. Well 7 was drilled at the edge of the richest sand zone, and its reservoir also matches the seismic-based density results.

In addition, the sandy cap rock within the upper-McMurray in well 7 is predicted by the seismic inversion.

This study demonstrates the application of rock physics analysis and the determination of rock density from seismic data can be used to characterize heterogeneity with the McMurray formation portion of the Athabasca oil sands. While other rock physics parameters can be used, density seems to be a good indicator of lithology at this study site.

This same methodology can be applied to other areas where the objective is to determine heterogeneity within any formation of interest.

Comments (0)


Geophysical Corner - Rongfeng Zhang

Rongfeng Zhang is a senior geoscientist with Geomodeling Technology Corp.

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.

Arthur Barnes, an AAPG member, is with Landmark Graphics Corp., Highlands Ranch, Colo. He can be contacted at Landmark .

Yong Xu are with Arcis Corp., Calgary, Canada

Geophysical Corner - Alistair Brown

Alistair Brown, a consultant from Allen, Texas, is a former editor of the EXPLORER’s Geophysical Corner and in 2009 received an AAPG Distinguished Service award.

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

This article addresses the controls exerted by sedimentologic and diagenetic factors on the preservation and modification of pore-network characteristics (porosity, pore types, sizes, shapes, and distribution) of carbonates belonging to the Bolognano Formation. This formation, exposed at the Majella Mountain, Italy, is composed of Oligocene–Miocene carbonates deposited in middle- to outer-ramp settings. The carbonates consist of (1) grainstones predominantly composed of either larger benthic foraminifera, especially Lepidocyclina, or bryozoans; (2) grainstones to packstones with abundant echinoid plates and spines; and (3) marly wackestones to mudstones with planktonic foraminifera.

The results of this field- and laboratory-based study are consistent with skeletal grain assemblages, grain sizes, sorting, and shapes, all representing the sedimentologic factors responsible for high values of connected primary macroporosity in grainstones deposited on the high-energy, middle to proximal outer ramp. Cementation, responsible for porosity reduction and overall macropore shape and distribution in grainstones to packstones deposited on the intermediate outer ramp, was mainly dependent on the following factors: (1) amount of echinoid plates and spines, (2) grain size, (3) grain sorting and shapes, and (4) clay amount. Differently, in the wackestones to mudstones, laid down on the low-energy, distal outer ramp, matrix is the key sedimentologic factor responsible for low values of scattered macroporosity and dominance of microporosity. The aforementioned results may be useful to improve the prediction of reservoir quality by means of mapping, simulating, and assessing individual carbonate facies with peculiar pore-network characteristics.

Desktop /Portals/0/PackFlashItemImages/WebReady/sedimentologic-and-diagenetic-controls-italy.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3706 Bulletin Article

A series of short and steep unidirectionally migrating deep-water channels, which are typically without levees and migrate progressively northeastward, are identified in the Baiyun depression, Pearl River Mouth Basin. Using three-dimensional seismic and well data, the current study documents their morphology, internal architecture, and depositional history, and discusses the distribution and depositional controls on the bottom current–reworked sands within these channels.

Unidirectionally migrating deep-water channels consist of different channel-complex sets (CCSs) that are, overall, short and steep, and their northeastern walls are, overall, steeper than their southwestern counterparts. Within each CCS, bottom current–reworked sands in the lower part grade upward into muddy slumps and debris-flow deposits and, finally, into shale drapes.

Three stages of CCSs development are recognized: (1) the early lowstand incision stage, during which intense gravity and/or turbidity flows versus relatively weak along-slope bottom currents of the North Pacific intermediate water (NPIW-BCs) resulted in basal erosional bounding surfaces and limited bottom current–reworked sands; (2) the late lowstand lateral-migration and active-fill stage, with gradual CCS widening and progressively northeastward migration, characterized by reworking of gravity- and/or turbidity-flow deposits by vigorous NPIW-BCs and the CCSs being mainly filled by bottom current–reworked sands and limited slumps and debris-flow deposits; and (3) the transgression abandonment stage, characterized by the termination of the gravity and/or turbidity flows and the CCSs being widely draped by marine shales. These three stages repeated through time, leading to the generation of unidirectionally migrating deep-water channels.

The distribution of the bottom current–reworked sands varies both spatially and temporally. Spatially, these sands mainly accumulate along the axis of the unidirectionally migrating deep-water channels and are preferentially deposited to the side toward which the channels migrated. Temporally, these sands mainly accumulated during the late lowstand lateral-migration and active-fill stage.

The bottom current–reworked sands developed under the combined action of gravity and/or turbidity flows and along-slope bottom currents of NPIW-BCs. Other factors, including relative sea level fluctuations, sediment supply, and slope configurations, also affected the formation and distribution of these sands. The proposed distribution pattern of the bottom current–reworked sands has practical implications for predicting reservoir occurrence and distribution in bottom current–related channels.

Desktop /Portals/0/PackFlashItemImages/WebReady/upper-miocene-to-quaternary.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3665 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 4317 CD-DVD

See Also: Map

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

See Also: Online e Symposium

Desktop /Portals/0/PackFlashItemImages/WebReady/oc-es-eagle-ford-shale-prospecting-with-3d-seismic-data.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 1482 Online e-Symposium