A primary factor in controlled lab work

Testing Curvature’s Impact on Fractures

Contributors: Kurt Marfurt, Evan Staples

Fracture zones can be critical to improving or creating sufficient porosity and permeability in hydrocarbon reservoirs – with strain, along with lithology and thickness being the major controls.

In this article we will use layer curvature as a proxy for strain induced by layer bending.

We test our hypothesis using small-scale clay models, which provide a time history of folding and fracturing, and are commonly considered analogs for large-scale field structures. The experimental apparatus consisted of a horizontal table with one moveable sidewall, one stationary sidewall and a deforming base.

We present four experiments – one extensional and three compressional.

In the extensional experiment, the clay cake was placed on top of two rigid, thin metal plates that were moved away from each other.

In the compressional experiments, the clay cake was placed on two metal wedges with inclinations 45 degrees, 30 degrees and 15 degrees, that moved toward each other to generate reverse basement faulting.

We use wet clay with dimensions of 20 centimeters long, 15 centimeters wide and five centimeters thick. A laser scanner positioned above the clay cakes captured 3-D surface images.

A typical experiment lasted approximately 30 minutes with a laser scan occurring every two minutes; this short duration eliminated clay drying as a variable.

The curvature of the clay surface was calculated from the laser scans using commercial software, and the fractures were mapped on digital photographs of the clay surface. Curvature over the deformation area was calculated in each stage by placing three polygons at fixed locations and averaging the curvature within each polygon.

The fracture intensity (FI) was calculated for these polygons by dividing total measured fracture length in each polygon by its area.

In the extensional experiment, a basin formed above the moving plates with a synthetic normal fault on one side, and a system of antithetic-normal faults, fractures and a flexure on the other side (top of figure 1). Curvature and fracture intensity measurements were conducted on the later side of the basin.

The bottom of figure 1 displays a suite of positive curvature images from initial through final stages.

Fractures were first visible on the clay surface with a measured curvature of 2.53 x 10-3 cm-1. We observed that curvature and fracture intensity increased with time, but was dependent on fault movement. Normal faults started as individual segments and eventually joined up to form faults that spanned the entire clay model.

At a critical point, fractures no longer accommodated the displacement; rather, normal faults accommodated displacement in the model.

At this point, we ended the experiment.

In the three compressional experiments, an anticline developed with fractures on the crest sub-parallel to the axial plane (top of figure 2). Fractures were first visible on the clay surface with curvature ranging between 1.40 x 10-2 cm-1 to 2.13 x 10-2 cm-1.

The bottom of figure 2 shows a suite of positive curvature images initial to final stages for the 30-degree ramp.

During the experiments we observed that curvature and fracture intensity qualitatively increased with time as initially hypothesized. Fractures initially accommodated the extensional displacement on the top of the anticline in individual segments; as the experiments progressed, the fracture segments began to connect and eventually generated faults to accommodate the increased displacement.

At this point, fractures and curvature were no longer the principle features that accommodated displacement, and we stopped the experiments.

Curvature increased over time with systematic deformation occurring in the extensional, 15-degree compressional ramp and 30-degree compressional ramp experiments. Deformation in the 45-degree compressional ramp was non-systematic, but followed a similar trend as that seen in the other two ramp experiments.

Maximum curvature values varied based on the experimental setting, with the extensional value (7.5 x 10-3 cm-1) being one order of magnitude lower than the compressional settings (2.28 x 10-2 cm-2 to 1.78 x 10-2 cm-1).

Calculated FI increased with curvature and correlations show a strong-positive-linear relationship (figure 3) – however, fracturing did not occur at the same curvature value in each experiment. The extensional experiments showed fracturing initiation at a significantly lower curvature value (2.53 x 10-3 cm-1) than compressional experiments, where fracturing initiated at values one order of magnitude higher than extensional experiments. (2.13 x 10-2 cm-1, 1.62 x 10-2 cm-1, and 1.40 x 10-2 cm-1 for 45-degree, 30-degree and 15-degree ramps, respectively).

We interpret differences in horizontal strain as the main reason there is such a disparity in fracturing and curvature values at fracture initiation. In the extensional experiment, we assume that horizontal strain occurred throughout the experiment, since the basement plates were constantly moving apart – though we only measured vertical displacement.

In the compressional settings, however, an anticline had to develop and grow before horizontal strain was great enough on its crest to induce fracturing.

The clay model experiment results suggest that curvature was the primary factor in fracture generation in controlled, laboratory settings.

Assuming results from clay models show fracture patterns similar to real rocks, then clay-model results imply that curvature and fracture intensity follow a 3-phase curve of elastic deformation (no fractures), through linear correlation (Figure 3), to fracture saturation (movement along faults) when subsurface deformation involves basement reverse faulting or layer extension.

In a future article we will show the correlation of curvature computed from 3-D seismic data to fractures seen in horizontal image logs.

Acknowledgements: We would like to thank AAPG member Ze’ev Reches from the University of Oklahoma for help with this project and guidance in clay model setup and analysis. Additionally, we would like to thank Schlumberger for providing software for this research at the University of Oklahoma.

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.

AAPG Member Evan Staples is with ConocoPhillips in Houston.

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: Book

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

See Also: Bulletin Article

Estimation of the dimensions of fluvial geobodies from core data is a notoriously difficult problem in reservoir modeling. To try and improve such estimates and, hence, reduce uncertainty in geomodels, data on dunes, unit bars, cross-bar channels, and compound bars and their associated deposits are presented herein from the sand-bed braided South Saskatchewan River, Canada. These data are used to test models that relate the scale of the formative bed forms to the dimensions of the preserved deposits and, therefore, provide an insight as to how such deposits may be preserved over geologic time. The preservation of bed-form geometry is quantified by comparing the alluvial architecture above and below the maximum erosion depth of the modern channel deposits. This comparison shows that there is no significant difference in the mean set thickness of dune cross-strata above and below the basal erosion surface of the contemporary channel, thus suggesting that dimensional relationships between dune deposits and the formative bed-form dimensions are likely to be valid from both recent and older deposits.

The data show that estimates of mean bankfull flow depth derived from dune, unit bar, and cross-bar channel deposits are all very similar. Thus, the use of all these metrics together can provide a useful check that all components and scales of the alluvial architecture have been identified correctly when building reservoir models. The data also highlight several practical issues with identifying and applying data relating to cross-strata. For example, the deposits of unit bars were found to be severely truncated in length and width, with only approximately 10% of the mean bar-form length remaining, and thus making identification in section difficult. For similar reasons, the deposits of compound bars were found to be especially difficult to recognize, and hence, estimates of channel depth based on this method may be problematic. Where only core data are available (i.e., no outcrop data exist), formative flow depths are suggested to be best reconstructed using cross-strata formed by dunes. However, theoretical relationships between the distribution of set thicknesses and formative dune height are found to result in slight overestimates of the latter and, hence, mean bankfull flow depths derived from these measurements.

This article illustrates that the preservation of fluvial cross-strata and, thus, the paleohydraulic inferences that can be drawn from them, are a function of the ratio of the size and migration rate of bed forms and the time scale of aggradation and channel migration. These factors must thus be considered when deciding on appropriate length:thickness ratios for the purposes of object-based modeling in reservoir characterization.

Desktop /Portals/0/PackFlashItemImages/WebReady/deposits-of-the-sandy-braided-saskatchewan-river.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3712 Bulletin Article

See Also: GIS Open File

Courtesy of  Dr. Ricardo J. Padilla y Sánchez for allowing the GIS project of the Tectonic Map of Mexico 2013, to be posted as a GIS Open file. This compilation shows key geologic information about Mexico.

Desktop /portals/0/images/gis-open-files/Tectonic-Map-of-Mexico.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 10137 GIS Open File

See Also: Learn! Blog

Designed to expose geoscientists to the language and methodologies commonly employed by reservoir engineers to estimate oil and gas reserves, this course will include a brief review of reservoir rock properties, introduction to fluid properties, reservoir drive mechanisms, and commonly used techniques for estimating reserves including analogy, volumetric analysis, material balance, decline curve analysis, and computer simulation.

Desktop /Portals/0/PackFlashItemImages/WebReady/Fifth-Annual-AAPG-SPE-Deepwater-Reservoirs-Geosciences-Technology-Workshop.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 13612 Learn! Blog

See Also: Online e Symposium

The geochemistry of formation fluids (water and hydrocarbon gases) in the Uinta Basin, Utah, is evaluated at the regional scale based on fluid sampling and compilation of past records.

Desktop /Portals/0/PackFlashItemImages/WebReady/oc-es-green-river-shales-geochemical-basin-study.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 1436 Online e-Symposium