Analysis of stress field and determination of safe mud window in borehole drilling (case study: SW Iran)

Abstract

Having information on the magnitude and direction of stress in high depth is of critical importance in geology and engineering sciences. Of uppermost applicants of stress data in oil industry, determination of borehole stability can be mentioned. To achieve scientific and engineering solutions for problems such as optimal mud weight, stable routes, casing, sand production, etc., precise information on stress position in depth is required. So, the present study is aimed to determine magnitude and directions of in situ stresses in a borehole located in south western Iran. Minimum horizontal stress is observed along the direction of borehole breakout and maximum horizontal one, along the direction of derived from drilling. Direction of these two horizontal stresses is determined by means of image log. Position of in situ stresses in the under-study field is of normal stress regime. Having determined direction and magnitude of in situ stresses, we find safe mud window to design a stable borehole.

J Petrol Explor Prod Technol (2013) 3:105110 DOI 10.1007/s13202-013-0053-2

ORIGINAL PAPER - EXPLORATION ENGINEERING

Analysis of stress eld and determination of safe mud window in borehole drilling (case study: SW Iran)

Mohammad Abdideh Mohammad Reza Fathabadi

Received: 15 November 2011 / Accepted: 3 March 2013 / Published online: 29 March 2013 The Author(s) 2013. This article is published with open access at Springerlink.com

Abstract Having information on the magnitude and direction of stress in high depth is of critical importance in geology and engineering sciences. Of uppermost applicants of stress data in oil industry, determination of borehole stability can be mentioned. To achieve scientic and engineering solutions for problems such as optimal mud weight, stable routes, casing, sand production, etc., precise information on stress position in depth is required. So, the present study is aimed to determine magnitude and directions of in situ stresses in a borehole located in south western Iran. Minimum horizontal stress is observed along the direction of borehole breakout and maximum horizontal one, along the direction of derived from drilling. Direction of these two horizontal stresses is determined by means of image log. Position of in situ stresses in the under-study eld is of normal stress regime. Having determined direction and magnitude of in situ stresses, we nd safe mud window to design a stable borehole.

Keywords Reservoir rock Geomechanics Stress eld

Well bore stability Safe mud window Petrophysical data

Introduction

Hydrocarbon geomechanic reservoirs play daily-increasing important role in evaluation and development of oil and gas elds. Geomechanics is, indeed, a eld of science exploring and analyzing earth behavior against stresses. Stresses are of two types: natural stress inside the earth and the

other, those induced by men during various operations like drilling. So, as implied by the term Geomechanics, keys to physical and mechanical problems resulted from materials operating in various depth of the earth and from the existing stresses can be found in domain of this science. However, chemical disturbance and problems related with it, should be studied and investigated by other sciences especially geo mechanics (Amadie and Stephansson 1997). For the purpose of analyzing borehole instability, direction and magnitude of in situ is of great criticality.

In situ stresses are usually expressed in terms of a tensor of six components. However, in oil eld studies, the stresses are determined only for three principal components of the stress: a vertical component (rv), and two minimum and maximum horizontal stress components (rH and rh). Accordingly, a number of direct and indirect methods (empirical equations) have been presented to determine the magnitude of stresses in oil reservoirs in great depth. However, the principal horizontal stresses, especially rH,

are still major challenge to be determine in geomechanical oil eld studies (Zobac et al. 2003).

During drilling, we may encounter two main problems in borehole stability: breakout (BO) and drilling induced tensile fracture (DITF) that may lead to shing, stuck pipe, sidetracking, reaming and uid loss (Nguyen et al. 2004), all of which can be resolved through determination of safe mud window for borehole drilling. Safe mud window consists of a pressure with the magnitude between pore pressure and minimum horizontal stress. In situ stress eld should be analyzed so as to determine safe mud window. Perfect illustration of stress eld in the earth only requires magnitude and direction of main stresses to be determined.

Minimum horizontal stress is observed along the direction of borehole breakout and maximum horizontal one, along the direction of DITFs derived from drilling.

M. Abdideh (&) M. R. Fathabadi

Department of Petroleum Engineering, Omidiyeh Branch, Islamic Azad University, Omidiyeh, Irane-mail: [email protected]

123

106 J Petrol Explor Prod Technol (2013) 3:105110

Direction of these two horizontal stresses is determined by means of image log. Magnitude of main stress is obtained through integration of density log. Due to the absence of leak-off test data, magnitude of minimum and maximum horizontal stresses is estimated via equations available.

Discussion

The eld under study consists of three reservoirs of Asmari, Bangestan and Khami, of which the rst two are oil reservoirs and the latter, a gas reservoir. The big Asmari reservoir is complicated and heterogeneous in terms of reservoir rock features. Lithologic status of the considered borehole in Asmari structure is determined applying cross plots of (photoelectric and Neutron) density logs. Both cross plots appropriately reveal that the dominant lithology is carbonate (dolomite enjoys higher percentage than that of limestone); furthermore, amount of sandstone found is insignicant when compared to the abundance of carbonate. Figure 1 shows the above-mentioned cross plots.

Direction of main stresses

Breakouts and DITFs of borehole walls appropriately reveal stress direction. In almost vertical boreholes, BO axis is straightly along minimum horizontal stress and DITF, along maximum horizontal stress (Tingay 2008). Therefore, these fractures can be applied to determine stress direction. Fracture direction can also be determined by means of geophysical logs such as televiewer ultrasonic borehole (BHTV), formation micro imager (FMI), and Caliper Logs (Fig. 2).

Magnitude of in situ stresses

Assuming verticality of one component of the main stress, vertical boreholes are parallel to vertical main stress (Sv).

Therefore, SHmax and Shmin will be two other main stresses.

Vertical stress is obtained through integration of rock density from surface to the considered depth:

sv Z

z

qzgdz qgz: 1

where qz is rock density which itself is a function of

depth; g represents gravity constant and q stands for the average density. In the under study borehole, q has been assumed equal to 2.56 g

cm3.

Magnitude of shmin can be measured via different

methods like hydraulic fracture, leak-off tests, micro-fracture test and mini-fracture test. Willis and Hubbert (1957) presented a thorough discussion on that hydraulic fracture, in which a fracture developed, extended vertical to the direction of minimum horizontal stress. The reason they offered was that the work done to open a crack of a given magnitude was proportionate to the stress vertically imposed to crack plane against opening magnitude.

Another method used to calculate the horizontal stresses is based on the poroelastic theory. In a tectonically active basin, tectonic stresses and strains arise from tectonic plate movement. If tectonic strains are applied to rock formations, these strains add a stress component in an elastic rock. The poroelastic horizontal strain model takes tectonic strains into account, and therefore accommodates anisotropic horizontal stresses (Blanton and Olson 1999).

rh

0

m1 m

rv

m1 m

app app

E1 m2

ex

mE1 m2

ey

2

rH

m1 m

rv

m1 m

app app

E1 m2

ey

mE1 m2

ex

3

where rh is minimum horizontal stress, rH represents maximum horizontal stress and m stands for Poissons ratio, rv is vertical stress, a indicates Biot coefcient, Pr stands for pore pressure and E for youngs modulus, and nally ex and ey present strain toward minimum and maximum horizontal stresses, respectively.

Elasticity coefcients can be obtained applying empirical relations and also measurement well logs. Gradient of the average pore pressure in the under examination bore-hole equals 0.365 psi/ft; Biot coefcient amounts to 1 and ex and ey are considered 1.5 and 1, respectively.

mdyn

1=2 Dts

Dtc

2

Dts Dtc

2

1

1 : 4

Fig. 1 Variation of lithology in oil eld

123

J Petrol Explor Prod Technol (2013) 3:105110 107

Fig. 2 Identication of DITF and BO using FMI log

Edyn

qb3 4

2

Dt2s Dt2c

5

Dtc

Dts

where Dts is the slowness of shear wave s

km

, Dtc rep-

resents slowness of compression wave s

km

, qb stands for

density g

cm3

, mdyn is dynamic Poissons ratio, and Edyn indicates dynamic youngs modulus (Gpa).

Dynamic data cannot directly be utilized to develop mechanical models. So, they should be rst converted into static data through some calculation changes made and then used in geomechanical model. Poissons ratio and static youngs modulus are both calculated via the following relations in SW Iran. The results show good conformity with laboratorial data:

msat

0:7mdyn 6 Esta 0:414 Edyn 1:0593 7

Borehole stability

Borehole instability factors are generally of three types: mechanical factors, chemical factors or a combination of these two. Mechanical factors mostly depend on inappropriate mud weight (too light or too heavy) and improper drilling methods (drilling extent, pipes moment and vibration, and lifting pipes up and down), whereas chemical factors are dramatically inuenced by drilling mud, that is to say, improper mud and insufcient inhibitors. Of course, in most cases borehole instability is the simultaneous result of chemical and mechanical factors (Hawkes and Mclellan 2002).

Above all, borehole instability relies on stress distribution and centralization around the drilled borehole.

In case rock resistance is more than stresses induced, borehole will be a stable one; otherwise, there will be a fracture. Furthermore, borehole instability lengthens drilling period and consequently increases borehole drilling costs.

For a safe mud window to be obtained, mud weight should be something between pore pressure and minimum horizontal stress. Mud weight less than structure pressure results in uid eruption and uid ows from reservoir into the well. In the event mud pressure exceeds minimum horizontal stress, DITFs develop and we will have mud loss. Range of stable mud window is something between least drilling mud weight (to prevent breakout) and minimum horizontal stress. From geomechanical point of view, stable mud window keeps borehole safe against DITFs and/ or stuck pipes which happen due to high mud weight as well as BOs which result from low mud weight (Al-Ajmi and Zimmerman 2006). Safe and stable mud weight windows have been shown in Fig. 3.

Fig. 3 Safe and stable mud weight windows

123

108 J Petrol Explor Prod Technol (2013) 3:105110

Stress around a vertical borehole

Underground structures are always under pressure because of vertical and tectonic stresses. When a borehole is drilled inside a structure, some parts of (rock) materials are missed. Boreholes side walls are maintained only by uid pressure. Since this uid pressure is not consistent with in situ structure stresses, new stress redistribution occurs around borehole that may lead to rock fracture. Therefore, identication of stresses existent around a borehole is of great importance to examine borehole problems. There are two common types of borehole fracture: rst one is drilling induced tension fracture (DITF) and the second one, with 90 of difference, is called breakout. In a substance with linear elastic behavior, centralization of bigger stress occurs on borehole wall, so borehole fracture is expected to start from this area.

Drilling operation disturbs the equilibrium condition of in situ stresses. Therefore, new (induced) stresses are provided around the wellbore. If the wellbores are affected by induced stresses, they have some failures. Thus, after the estimation of the principal stresses, they are converted into the induced stresses using existing relations. Stress concentration around a vertical wellbore, which has been drilled in parallel with the principal vertical stress and is in isotropic conditions in elastic environment, was described by Kirsch equations (Kirsch 1898).

rrr

Fig. 4 Variation of effective principal stresses around a vertical wellbore as a function of azimuth

Table 1 Multiple modes of wellbore failure

Multiple modes of tensile failure

Modes

rzz = r3 D = rzz - T [ 0

Pm [ - [rV ? rh

(2t3) ? rH (2t3)]

horizontal (HOR)

rhh = r3 D = rhh - T [ 0

Pm \ 3rh - rH

-PP - T

Vertical (VER):

1

2 rH rh 2Pp 1

R2 r2

1

2 rH rh

Cylindrical (CYL)

Multiple modes of compressive failure

r1 r2 r3 Modes

Rhh rzz rrr Wide breakout (WBO):

conventional breakout

rhh rrr rzz Low angle echelon (LAE):

requires high mud weights.

Failed rock will not fall into the wellbore as rrr : r2

rzz rrr rhh High angle echelon (HAE):

forms on opposite side of well as a conventional breakout but the failed rock will not fall into the wellbore as rrr : r2

rzz rhh rrr Shallow knockout (SKO): results in failure all the way around the wellbore

rrr rzz rhh Narrow breakout (NBO): requires unreasonably high mud weights.

rrr rhh rzz Deep knockout (SDKO):

requires unreasonably high mud weights

r1, r2 and r3 maximum, mean and minimum stresses, rH and rh maximum and minimum horizontal stresses, h angle measured from the azimuth of rH, PP and Pm pore and mud pressure, rhh, rrr and rzz tangential, redial and vertical effective stresses

rrr=r3 D = rrr - T [ 0

Pm [ PP ?T

1 4

R2r2 3

R4 r4

cos 2h

DPR2r2 8

rhh

1

2 rH rh 2Pp 1

R2 r2

1

2 rH rh

1 3

R2 r2

cos 2h

DPR2r2 9

rzz

rV 2mrH rh cos 2h PP 10 where h: angle measured from the azimuth of rH, rhh:

effective tangential stress, rrr: effective redial stress and rzz: effective vertical stress.

rDT represents thermal stresses arising from the difference between the mud temperature and formation temperature. This will be ignored for the moment.

Mathematically, effective stresses around a vertical borehole of the radius R are dened in cylinder coordinate system of axis as follows:

rhh max 3SH

max

Sh

min

Pw Pp 11

rhh min 3Sh

min

SH

max

Pw Pp 12

rzz max Sv 2mSH max

Sh

min

Pp 13

123

J Petrol Explor Prod Technol (2013) 3:105110 109

rzzmin Sv 2mSH max

Sh

Pp 14

Considering fracture mechanics science, DITFs vertical to the main stress are minimum, and maximum along main stress (Fig. 4). If tangential stress is less than rock tension resistance, we will have a DITF and if it is greater than rock pressure resistance, a breakout will occur.

Depending on the status of the three components of stresses, 6 modes of shear failure, and 3 modes of tensile failure may occur around the wellbore (Table 1) (Fjaer et al. 2008). Because of the typical depth of oil reservoirs, the most common types of failure are of breakout and vertical tensile failure.

Conclusion

In this paper, the stresses were determined using the poroelastic method and based on petrophysical data. The result of this method shows relative compliance with the stress that was calculated on the basis of the width of the breakout and rock strength. The wellbores affected by these stresses are mostly stable, with no severe and extensive failures. The safe mud weight window was determined based on the estimated stress proles. This mud window showed that in some depth, the used mud pressure is not appropriate and cause the limited failures.

It should be noted that in the studied eld, there are other important factors in causing the failure such as collisions of the drill string with the wellbore, sudden decrease of drilling mud pressure and the presence of fractures.

Stress regime existing in under examination borehole is of normal type. Stress position and pore pressure, as shown in Figs. 5 and 6, are determined to the extent of safe mud window located somewhere between pore pressure and minimum horizontal stress. Considering type of stress eld, in case of horizontal drilling it is recommended that optimal drilling route is along minimum horizontal stress.

Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.

min

Fig. 5 Stress eld and determination of safe mud window in borehole

References

Al-Ajmi AM, Zimmerman RW (2006) Stability analysis of vertical boreholes using the Mogi-Coulomb failure criterion. Int J Rock Mech Mini Sci 43:12001211

Amadie B, Stephansson O (1997) Rock stress and its measurement,

1st edn. Chapman and Hall, LondonBlanton TL, Olson JE (1999) Stress magnitudes from logs-effects of tectonic strains and temperature. SPE Res Eval Eng 2(1):6268 Fjaer E, Holt RM, Horsrud P, Raaen AM, Risnes R (2008) Petroleum related rock mechanics, 2nd edn. Elsevier, Amsterdam

Fig. 6 Induced stress around borehole

123

110 J Petrol Explor Prod Technol (2013) 3:105110

Hawkes CD, Mclellan PJ (2002) Coupled modeling of borehole stability and multiphase ow for under balanced drilling, SPE pp 74447

Hubbert MK, Willis DG (1957) Mechanics of hydraulic fracturing.

Pet Trans AIME 210:5363Kirsch G (1898) Die theorie der elastizitat und die bedurfnisse der festigkeitslehre. Z Ver Dtsch Ing 42(29):797807Nguyen Binh T, Tomochica T, Akihiko O (2004) In situ stress and pore pressure elds in the North Cuu Long Basin, Offshore Vietnam, SPE pp 87055

Tingay M (2008) Borehole breakout and drilling-induced fracture analysis from image logs. World stress map project. Guidelines: image logs

Zoback MD, Barton CA, Brudy M, Castillo DA, Finkbeiner T, Grollimund BR, Moos DB, Peska P, Ward CD, Wiprut DJ (2003) Determination of stress orientation and magnitude in deep wells. Int J Rock Mech Mini Sci 40:10491076

123