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

See Also: Bulletin Article

West Edmond field, located in central Oklahoma, is one of the largest oil accumulations in the Silurian–Devonian Hunton Group in this part of the Anadarko Basin. Production from all stratigraphic units in the field exceeds 170 million barrels of oil (MMBO) and 400 billion cubic feet of gas (BCFG), of which approximately 60 MMBO and 100 BCFG have been produced from the Hunton Group. Oil and gas are stratigraphically trapped to the east against the Nemaha uplift, to the north by a regional wedge-out of Hunton strata, and by intraformational diagenetic traps. Hunton Group reservoirs are the Bois d'Arc and Frisco Limestones, with lesser production from the Chimneyhill subgroup, Haragan Shale, and Henryhouse Formation.

Hunton Group cores from three wells that were examined petrographically indicate that complex diagenetic relations influence permeability and reservoir quality. Greatest porosity and permeability are associated with secondary dissolution in packstones and grainstones, forming hydrocarbon reservoirs. The overlying Devonian–Mississippian Woodford Shale is the major petroleum source rock for the Hunton Group in the field, based on one-dimensional and four-dimensional petroleum system models that were calibrated to well temperature and Woodford Shale vitrinite reflectance data. The source rock is marginally mature to mature for oil generation in the area of the West Edmond field, and migration of Woodford oil and gas from deeper parts of the basin also contributed to hydrocarbon accumulation.

Desktop /Portals/0/PackFlashItemImages/WebReady/Bulletin-hero-2013-07jul.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3770 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 4316 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 4448 CD-DVD

See Also: Short Course

This course will focus on practical techniques to investigate and optimize fracture treatments. Participants in this course will have access to more than 200 published field studies in which the productivity and profitability of fields have been improved by altering the treatment design.
Desktop /Portals/0/PackFlashItemImages/WebReady/ace2015-sc16-frac-engineer-hero.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 14631 Short Course