Open access peer-reviewed chapter

# Hydraulic and Sleeve Fracturing Laboratory Experiments on 6 Rock Types

By Sebastian Brenne, Michael Molenda, Ferdinand Stöckhert and Michael Alber

Submitted: August 1st 2012Reviewed: March 4th 2013Published: May 17th 2013

DOI: 10.5772/56301

## 1. Introduction

Hydraulic tensile strength is a crucial value for planning reservoir stimulation and stress measurements. It is used in the classical breakdown pressure (Pb) relation by Hubbert & Willis [1], where Pbis a function of major and minor principal horizontal stresses SHand Sh, hydraulic tensile strength σTand pore pressure P0:

Pb=3ShSH+σTP0E1

For hydraulic fracturing laboratory experiments (MiniFrac – MF) under isostatic confining pressure Pmthis might be reduced to:

Pb=cPm+σTP0E2

The coefficient cshould be equal to two when porepressure is neglected. However, many laboratory experiments [2,3] resulted in values of about 1 for c, which might be explained by poroelastic effects.

Thus, when poroelasticity is excluded in the experiments by taking dry samples and sealing off the central borehole by an impermeable membrane (like a polymer tube), one would expect that c equals two and σTwill be in the range of the tensile strength as determined by other tensile strength tests.

However, experiments with jacketed boreholes (sleeve MiniFrac – SMF) yield remarkable high values for c(about 6 to 8) and also for σT(about 3 to 5 times the tensile strength of the material) [4]. As a consequence we use a linear elastic fracture mechanics approach to evaluate our experiments.

### 1.1. Theory of hydraulic and sleeve fracturing on hollow cylinders

Fracture mechanics deal with stress concentrations around fractures and the definition of propagation criteria for fractures. The theory is essentially based on the works of Griffith [5] and Irwin [6], which led to the introduction of the stress intensity factor K.

K=σπaE3

Krepresents the magnitude of the elastic stress singularity at the tip of a fracture of the length 2asubjected to a uniform stress σ. With this concept, it is possible to formulate a simple fracture propagation criterion K=KC. The fracture propagates when Kreaches a critical value KC(fracture toughness) with the fracture toughness assumed to be a property of the rock.

Mode I stress intensity factors (KI)for arbitrary tractions (σ(x))applied to the surface of a fracture of the length 2amay be computed by following formula [7,8]:

KI=1πaaaσ(x)(a+xax)12dxE4

The direction of propagation is the x-axis and the stresses are applied perpendicular to the fracture. As can be seen from equation (4), KIincreases with growing fracture length. A simple, 2-dimensional model was assumed for determination of stress intensity factors at the crack tips of the hydraulically induced fractures in MF and SMF tests.

Two fractures of length aare radially emanating from a circular hole of radius rin an infinite plate subjected to a compressive far field stress of the magnitude  Pm. A fluid pressure Pinjis acting on the borehole wall and the pressure inside the fractures is either zero (SMF) or equal to the pressure in the borehole (MF:  Pfrac = Pinj). Stress intensities on the fracture tips can be determined by superposition of stress intensity factors resulting from each loading type [2,3]:

KIMF=KI(Pm)+KI(Pinj)+KI(Pfrac)E5
KISMF=KI(Pm)+KI(Pinj)E6

KIMF/SMFare not only dependent on the fracture length a(cf. Equations (3) and (4)) but also on the borehole radius r(see Appendix).

KIMF(full pressure in the fracture) gives an upper bound for stress intensities in this geometry (actual KIMF might be lower due to a negative pressure gradient inside the fracture), while KISMFis only induced by the pressure in the borehole and far-field stresses and is therefore substantially lower than KISMF(Figure 1).

As an analytical solution for KI(Pm)and KI(Pinj)for the ring geometry (corresponding to the hollow cylinder) is quite complex, we used the simpler solutions for a circular hole in an infinite plate as described by Rummel and Winter [2,3] (cf. Appendix). We compared the results of numerical simulations for the ring geometry with analytical solutions for the infinite plate. These results indicate that the simplification might be valid for fracture lengths smaller then aRr10with R=10r(Ris the outer radius of the ring geometry (cf. Figure 1).

Solving KIMFand KISMFfor Pinjand setting KIMF=KISMF=KIC(mode I fracture toughness) yields a critical injection pressure (PC(a))for each crack length a. If Pinjreaches PC (a), the fracture will propagate. From Figure 2 it can be seen, that PC (a)is very large for very small crack lengths. In consequence, the presence of microcracks is required for the formation of macroscopic fractures.

MF-equation (Equation 13) with full injection pressure in the fracture yields unstable fracture propagation at constant injection pressures as soon as microcracks start to propagate. On the other hand, the SMF-equations (Equation 14) show a minimum. Thus, after a fracture reaches the crack length corresponding to the minimum critical injection pressure, stable fracture propagation (i.e. to propagate the fracture, the injection pressure has to be increased) could be expected.

To calculate the coefficient cfrom Equation 2, we assume the presence of microcracks of a fixed length a0in the sample. The corresponding PC(a0)versus Pmfor the MF case (pressure in fracture = injection pressure) yields a coefficient c=1, which is independent of a0. PC(a0)while for SMF the c value depends strongly on the assumed microcrack length a0and gives c>2(increasing a0yield higher c).

## 2. Sample preparation and rock testing

The core specimens are drilled either with 40 mm or 62 mm water cooled diamond core drills. Core end planes are cut with a water flushed diamond saw blade and ground coplanar to a maximum deviation of ± 0.02 mm. The length and diameter ratio is chosen between 1.5:1 and 2.25:1. After sample preparation core specimens were dried for two days at a temperature of 105°C. For calculations of porosity Φ, measurements of bulk density ρdand of grain density ρsvia pycnometer were done. Static geomechanical parameters were determined by uniaxial and triaxial compressive as well as Brazilian disc tensile strength test series according to ISRM and DGGT suggested methods [9,10]. Mode I fracture toughness was determined using the Chevron notched three-point bending test according to [9]. Furthermore, a dynamic rock parameter, the compressional ultra-sonic wave velocity (vp) was measured. For MF/SMF specimens a central axial borehole was drilled into cores, using a water flushed diamond hollow drill with an outer diameter of 4 mm or 6 mm.

 rock type era & period quarry localization Microstructure marble TriassicUpper CarraraItaly coarse monocrystalline polygonal fabric limestone Jurassic upper Malm TreuchtlingenSouth Germany micritic limestone with abundant fossils and stylolites sandstone CarboniferousMississippian Dortmund/HagenWest Germany fine-grained arcose andesite D PermianRotliegend DoenstedtN German Basin porphyric fine-grained partly altered and pre-fractured rhyolite PermianRotliegend FlechtingenN German Basin porphyric fine-grained partly pre-fractured and sealed joints andesite R PermianRotliegend Thuringian ForestRotkopf porphyric coarse-grained and pre-fractured

### Table 1.

Rock types used in our experiments.

### 2.1. Stress field and injection

Figure 3 shows schematically the components of the MF and SMF experimental set-up. The stress field is induced by a hydraulic ram (capacity 4500 kN) through a servo controlled MTS Test Star II system with a Hoek triaxial cell which is pressurized using a hand pump to achieve simultaneous pressure increase of confining pressure and axial load. In all tests axial stress is set to be 2.5 MPa higher than Pmto prevent leakage. Distilled water is pumped into borehole as the injection fluid (MF) or into a polymer tube inside the borehole (SMF). A servo controlled pressure intensifier with a maximum injection pressure of 105 MPa was used to perform a constant pumping rate of 0.1 ml/s. With this apparatus also steady-state flow tests were conducted to obtain rock permeability values (according to the procedure described in [11]).

### 2.2. Acoustic emission monitoring

Acoustic Emission (AE) signals are acquired with an AMSY5 Acoustic Emission Measurement System (Vallen Systeme GmbH, Germany) equipped with 5 Sensors of type VS150-M. The VS150-M Sensors operate over a frequency range of 100-450 kHz with a resonance frequency at 150 kHz. Due to machine noise in the range below 100 kHz incoming signals are filtered by a digital bandpass-filter that passes a frequency range of 95-850 kHz. AE data are sampled with a sampling rate of 10 MHz. The sensors are fixed using hot-melt adhesive to ensure best coupling characteristics. Pencil-break tests (Hsu-Nielsen source [12]) and sensor pulsing runs (active acoustic emission by one sensor) are used to test the actual sensor coupling on the sample.

## 3. Results

### 3.1. Petrophysical and mechanical parameters

An overview of the rock properties is given in Table 2. A wide range of low porosity/permeability rocks with KICfrom 1 to 2 MPamwere tested.

 rock type ρd [g/cm3] Φ[-] k[m²] vp [m/s] KIC[MPa⋅m] C/φ[MPa]/[°] Estat [GPa] σT[MPa] marble 2.71±0.002 0.40±0.08 1E-19 5.67±0.06 1.57±0.11 (N=3) 29/22 36.0±1.0 6.4± 1.5 limestone 2.56±0.008 5.64±0.04 1E-18 5.59±0.05 1.19±0.14 (N=8) 27/53 32.2±1.6 8.2± 2.2 sandstone 2.57±0.006 4.39±0.06 8E-18 4.61±0.13 1.54±0.13 (N=4) 36/50 29.4±1.6 13.2± 2.1 rhyolite 2.63±0.015 1.02±0.12 9E-19 5.39±0.34 2.16±0.10 (N=4) 20…36/55 30.2±1.9 15.8± 3.2 andesite D 2.72±0.023 0.51±0.09 6E-19 5.26±0.28 1.90±0.08 (N=2) 20…41/50 28.7±3.1 14.6± 4.5 andesite R 2.60±0.013 1.70±0.08 4E-20 4.35±0.27 1.63±0.24 (N=4) 31/46 21.3±0.9 11.4± 2.8

### Table 2.

Mean values and standard deviations of petrophysical and mechanical parameters of tested rocks: dry bulk density ρd, porosity Φ, permeability k, compressive wave velocity vp, fracture toughness from Chevron notched three-point bending tests KIC, cohesion Cand friction angle φfrom a Mohr-Coulomb fit, Young’s modulus Estat, σTas determined by Brazilian disc tensile strength tests.

### 3.2. MF and SMF experiments

A schematic example of typical experiment data for MF and SMF tests is shown in Figure 4. Acoustic emission recordings are used to identify fracture processes in the test specimens. AE counts (threshold crossings per time interval – corresponding to AE activity) can directly be linked to localized fracture propagation [4]. The pressure at which the AE count rate raises rapidly is defined as PAE, which is further used as initial fracture propagation pressure. PAEis picked where the AE count rate permanently exceeds 1/10of the test’s average (see Figure 4).

In MF experiments, there is almost no AE activity prior to failure. Failure occurs in a very short time span just before sample breakdown (which occurs at maximum injection pressure Pinj max= Pb), therefore in MF experiments PAE Pb. In contrast, SMF experiments show an exponential increase in AE activity at injection pressures that are substantially lower than the actual breakdown pressure (PAE< Pb), but much higher than PAEin MF experiments. Therefore, it is possible to interrupt the experiment after AE activity started but before sample breakdown. The latter occurs in SMF experiments when the sample is completely splitted into two parts, which results in a tube breakdown and therefore in an injection pressure drop. Thin sections of specimens, where the experiment was interrupted, show macroscopic fractures emanating several millimeters into the sample but without any connection to the outer surface.

Noteworthy is the discrepancy between the MF and SMF initial fracture propagation pressures PAE at zero confining pressure. This result would imply different hydraulic tensile strength values for the same rock type when using equation (2). Furthermore there is a significant difference between the values of coefficient ccalculated for MF and SMF experiments. This can be seen clearly in Figure 5. Scale effects in PAE(Figure 2) with borehole radius are not evident for the 2 and 3 mm borehole radius samples due to data scattering. One single SMF test of a sandstone with a 6.35 mm borehole radius showed a significantly lower PAEas can be seen in Figure 5.

 MF SMF rock type Boreholediam. N PAE0[MPa] c N PAE0[MPa] c marble 4 mm 8 7.7 1.03 6 31.7 6.97 6 mm 8 9.4 0.96 4 19.6 8.54 limestone 4 mm 9 10.3 1.00 6 26.7 6.06 6 mm 8 8.2 1.01 7 29.1 5.79 sandstone 4 mm 8 18.2 1.13 5 41.7 6.29 6 mm 8 18.5 1.14 4 40.5 7.26 rhyolite 4 mm 11 18.2 0.89 4 51.6 6.04 6 mm 8 16.0 0.85 5 50.9 5.88 andesite D 4 mm 9 16.1 1.00 3 64.2 4.17 6 mm 6 10.9 0.87 4 48.1 4.83 andesite R 4 mm 10 10.0 1.17 4 47.4 6.26 6 mm 6 8.2 1.17 5 29.7 7.44 ∑ 93 - Ø 1.02 ∑ 57 - Ø 6.33

### Table 3.

Results of all MF and SMF rock type test series in form of PAE0and coefficient c (see equation (2)). N gives the number of tested samples per lithology and borehole diameter.

## 4. Conclusion

With SMF tests, stable fracture propagation was achieved over a wide range of injection pressure. Fracture initiation can be confidently linked to the AE count rates. This can be concluded from experiments that were interrupted after PAEbut below breakdown pressure. Physical examination revealed the presence of distinct fractures in these specimens (see Figure 6).

Due to high data scatter, the theoretical scale effect (critical injection pressure Pcis higher for smaller borehole radii) cannot be resolved by our data. However, tests with a larger (r=6.35mm) borehole give some support to the notion.

The simple fracture mechanics model is able to explain the higher PAEin SMF experiments. Equations 5 and 6 include the influence of fractures (with or without pressure inside), which is omitted in the classical approach (Equation 1). The high coefficient c in SMF test can only be explained by assuming high microcrack lengths (a06 mm).

We excluded poroelastic effects in our analysis due to the use of initially dry rocks with low permeabilities.

## Appendix

Superposition of stress intensity factors for two radial cracks of length aemanating from an internally pressurized (Pinj- injection pressure in the borehole, Pfrac- pressure inside the fracture) circular hole of radius rin an infinite plate subjected to an isostatic far-field stress Pmas described by [2] and [3] :

KI(Pm)=Pmr*fPm(a,r)E7
KI(Pinj)=Pinjr*fPinj(a,r)E8
KI(Pfrac)=Pfracr*fPfrac(a,r)E9
fPm(a,r)=2(1+ar)2((1+ar)21(π(1+ar))7)12+(π(1+ar))12(12πsin1(11+ar))E10
fPinj(a,r)=(1.3ar1+(1+ar)32+7.8(sin(2ar))2(1+ar)521.7) E11
fPfrac(a,r)=(π(1+ar))12(12πsin1(11+ar) )E12

Note: In equations 10 and 12 the borehole was excluded from the integration of stresses (cf. equation 4). The critical fracture propagation pressure at a given fracture length a, borehole radius rand mode I fracture toughness KICfor the unjacketed (PcMF)and the jacketed (PcSMF)case:

PcMF=1fPinj*fPfrac(KICr+PmfPm)E13
PcSMF=1fPinj(KICr+PmfPm)E14

## Acknowledgments

The authors wish to thank the German Federal Ministry for the Environment, Nature Conservation and Nuclear Safety for financing our project (FKZ 0325279B). Many core specimens were prepared and analyzed by our student staff: T. Hoferichter, J. Braun, S. Hönig, K. Bartmann and A. Kraft. A great praise to the precision mechanics workshop guys for the construction of the fine working pressure intensifier system. We appreciate fruitful discussions with geomecon GmbH, Potsdam.

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2013 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Sebastian Brenne, Michael Molenda, Ferdinand Stöckhert and Michael Alber (May 17th 2013). Hydraulic and Sleeve Fracturing Laboratory Experiments on 6 Rock Types, Effective and Sustainable Hydraulic Fracturing, Andrew P. Bunger, John McLennan and Rob Jeffrey, IntechOpen, DOI: 10.5772/56301. Available from:

### chapter statistics

3Crossref citations

### Related Content

#### Effective and Sustainable Hydraulic Fracturing

Edited by Rob Jeffrey

Next chapter

#### Microseismic Monitoring Developments in Hydraulic Fracture Stimulation

By Mirko van der Baan, David Eaton and Maurice Dusseault

First chapter

#### Introduction to Infrared Spectroscopy

By Theophile Theophanides

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.