SPE/ISRM 78163 Mud Chemistry Effects On Time-Delayed Borehole Stability Problems in Shales
SPE/ISRM 78163 Mud Chemistry Effects On Time-Delayed Borehole Stability Problems in Shales
SPE/ISRM 78163 Mud Chemistry Effects On Time-Delayed Borehole Stability Problems in Shales
Pore pressure. If the fluid pressure outside the shale is initial effect of the chemical activity is controlled by pore
different from the initial pore pressure, there will be a pressure diffusion, while the time scale of the relaxation
disturbance in pore pressure that propagates by diffusion from process is controlled by ion diffusion, which is usually a
the surface into the rock. This may give rise to a delayed slower process. Thus, even if the membrane efficiency is zero,
instability in a drilling situation: If the well pressure is higher there will be a temporary effect of the chemical activity.
than the pore pressure, the drilling fluid provides a support for Traditionally3,4, the effects of mud chemistry have only
the borehole wall, and may prevent failure. As the pore been assigned to the chemical activity. It has been shown
pressure inside the rock gradually increases towards the well experimentally however5, that some ions typically cations
pressure (provided there is no isolating membrane or like K+ and Ca2+ may be exchanged with Na+ in the clay
mudcake at the borehole wall), the support of the well mineral lattice structure, resulting in a change in the size of the
pressure will fade away, and the borehole wall may become structure. Thus the rock may expand or shrink as these alien
mechanically unstable. ions move in. The consequences for a borehole are similar to
Temperature. If the temperature of the fluid outside the the effect of a difference in temperature, giving a shift in the
rock is different from the initial temperature of the rock, the axial stress and the hoop stress at the borehole wall:
step in temperature will propagate by diffusion from the
surface into the rock. Rocks, like most other materials, expand E
or shrink when the temperature is changing. In a borehole z = = ch (3)
geometry, thermal expansion is restricted (except for the radial 1
direction). A difference between the mud temperature (Tm) and
the formation temperature (Tf) will therefore result in changes ch is the chemically induced strain (positive in contraction).
in the axial stress (z) and the hoop stress () at the borehole It is assumed to be a function of the cation exchange capacity
wall, given as2 and the relative amount of swelling clay minerals (e.g.
smectite) in the rock. The degree to which the ion exchange
process is expected to take place is assumed to depend on the
E
z = = Th (Tm T f ) (1) level of concentration of alien ions in the pore fluid of the
1 shale, up to a certain limit where the lattice is saturated. The
ion exchange process is considered to be much faster than the
E is the Youngs modulus, the Poissons ratio, and Th the diffusion process that brings the ions into the shale; hence the
coefficient of thermal expansion. It is apparent from Eq. (1) latter process will effectively control the time scale for the
that a cooling of the drilling fluid will tend to reduce the axial total process.
stress and the hoop stress, and thus enable the hole to Diffusion. Each of the processes described above implies
withstand lower well pressures. transport of some quantity from the surface into the rock, and
Chemical effects. The pore fluid in the shale has a the transport mechanism is diffusion. Thus, the distribution in
chemical activity (ash) given by the type and amount of time and space of the processes are controlled by the diffusion
dissolved ions. If the surface of the rock is exposed to a fluid equation
with a different chemical activity (adf), water molecules will
move into or out of the shale to compensate for this Q
disturbance, thus creating a disturbance (p) in pore pressure = C D 2Q (4)
at the surface3:
t
where Q is the actual quantity (ion concentrations, pore
adf pressure or temperature) and CD is the diffusion constant.
p = coT f ln (2)
ash The consequences of these processes will therefore be
delayed relative to the time of exposure of the free surface,
and the delay is given by the diffusion constants and the size
This disturbance (commonly called osmotic pressure) will and geometry of the rock involved.
propagate into the rock like any other pore pressure
disturbance, with the same consequences. p may be positive Laboratory tests on small samples
or negative, depending on whether the activity of the fluid is Diffusion in shales is a rather slow process, due to the low
larger or smaller than the activity of the shale. If adf < ash, we permeability. Typically, the diffusion constant is in the order
have that p < 0, which is beneficial for borehole stability of 10-8 10-9 m2/s for pore pressure diffusion, while for ion
since that increases the effective support of the drilling fluid diffusion the value is in the order of 10-10 10-11 m2/s. The
on the borehole wall. However, if the rock surface does not act typical time needed to propagate a distance L by diffusion is
as a perfect membrane (i.e. a membrane that allows only water given as
molecules to pass) ions will also move through the surface and
gradually change the activity of the shale towards that of the L2
= (5)
drilling fluid. After a while, the pore pressure shift p is CD
reduced to p, where is a number (called the membrane
efficiency) between zero and one. The time scale for the
SPE/ISRM 78163 MUD CHEMISTRY EFFECTS ON TIME-DELAYED BOREHOLE STABILITY PROBLEMS IN SHALES 3
Thus, the time needed for ion diffusion over a distance of concentrations of KCl make the shale shrink faster, however
19 mm (= radius of a standard 1.5 x 3 core plug) is typically the final amount of shrinkage is not correspondingly higher.
10 100 hrs for pore pressure diffusion, and 42 420 days for Thus the effect seems to reach a saturation level even at
ion diffusion. For a sample of radius 5 mm (which is the relatively low concentrations of the salt. The effect of
typical sample size used in the tests reported here) the reducing the salt concentration at a later stage is seen to be
corresponding diffusion time is reduced by a factor of 14. relatively low.
Notice that the ions will follow the water molecules if there is Quantitative analysis of these tests is not straightforward,
a pore pressure gradient causing the pore water to flow. This because the measured strain represents the integrated response
effect (called advection) will slow down or speed up the over the entire sample, while the various parts of the sample
overall movement of the ions, depending on the direction of respond differently according to the distance to the exposed
fluid flow versus the direction of ion diffusion. surface. As the deformation of the sample is actually resticted
Clearly, the aspect of time consumption favours the use of by the pistons, the sample must deform uniformly. This
smaller samples for laboratory studies of diffusion processes implies that parts of the sample will experience restricted
in shales. Furthermore, limited availability of core material deformation, such that the stress will vary over the sample,
and difficulties in sample preparation also encourage the use however in such a way that the total stress remains constant.
of smaller samples. The following test procedure has therefore Denoting as loc(r,t) the strain that would have been induced at
been established for fluid exposure test on shales6: the distance r from the centre of the sample at time t, if the
The shale sample is shaped like a small cylinder, typically deformation were not restricted, we find that the measured
10 mm in diameter and 2 mm high. The sample is positioned strain of the sample at time t, (t), is given as
between two pistons and subjected to a predefined, constant
uniaxial load, while the free surface of the sample is exposed R
2
to various fluids. The induced axial deformation of the sample ( t ) = 2 r loc ( r , t ) dr (6)
is measured as displacement of the pistons. The test geometry R 0
is shown in Fig. 1. Dehydration of the samples is carefully
avoided during the mounting procedure, to prevent capillarity where R is the sample radius. This implies that the measured
effects. response is dominated by the events close to the surface of the
In a typical test, the sample is first exposed to a 3.5 w% sample, since the part of the sample outside the radius
NaCl solution or sea water (which is assumed to be fairly
close to the salinity of the pore water) and left to stabilise for 1 RM = R / 2 will have the same impact as the part inside RM.
- 2 days. Then the sample is exposed to the solution of Furthermore, when the solution of the exposing fluid is
interest, and the response of the shale in terms of induced changed, several processes are triggered simultaneously. The
deformation is monitored. observed strain is thus also a sum of several effects. In order to
For the ion exchange process, the deformation of the facilitate quantitative analysis of the tests, a mathematical
sample is directly related to the chemically induced strain. For model, based on Eq. (4) for the actual geometry, has therefore
the effect induced by a contrast in chemical activity, the been established. The model estimates loc(r,t) due to pore
deformation of the sample is an elastic response to the change pressure diffusion (induced by chemical activity) and ion
in effective stress caused by the shift in pore pressure. Thus, exchange processes. By fitting the model to the measured
both effects have an impact on the sample that is detectable by strain data, the key parameters for the various processes can be
the measurement of sample deformation. extracted.
Results from tests on a field core material are shown in The modelled strain is shown in Fig. 2 as solid lines. The
Fig. 2. The material is a Tertiary shale from the North Sea, figure shows that the match between model and experiment is
with about 28% porosity7. The total clay content is about 44%, reasonably good. The key parameters extracted from these
and about half of this is smectite. Prior to the tests the samples KCl exposure tests were:
were stabilised in a 3.5 w% NaCl solution. In the first test
shown in Fig. 2, the sample was exposed to 5 w% KCl. It is CD (ion diffusion) = 4.210-11 m2/s
seen that the sample shrinks (the strain is positive) due to this Total strain amplitude ch = 0.0102
exposure. As the chemical activity of 5 w% KCl is almost Saturation level for ion exchange = 6.6 w%
exactly the same as for 3.5 w% NaCl, we can deduce that
chemical activity is not responsible for the observed effect, The theoretical value8 for the pore pressure diffusion
and we assume that it is caused by an ion exchange process5,7. constant is about 710-9 m2/s for a material with the actual
In the second test the sample was first exposed to 10 w% KCl. permeability and stiffness. For diffusion over a distance of 1.5
After about 24 hours the salinity was reduced to 5 w% KCl. In mm ( R - RM) the characteristic time is then about 5 minutes.
the third test the sample was first exposed to 20 w% KCl, Thus, on the time scale seen here, the chemical activity effect
followed by 5 w% KCl after about 24 hours. For the 10 w% occurs almost immediately. For the modeling, the theoretical
and 20 w% cases the activity of the water is lower than for the value was therefore used for the pore pressure diffusion
3.5 w% case (0.956 and 0.892 versus 0.982, respectively). constant.
Thus for these cases we expect to see effects of both chemical The membrane efficiency is expected to vary with the
activity and ion exchange. It is seen that the higher salinity of the exposing fluid. This has been accounted for, in
4 E. FJR, R.M. HOLT, O.-M. NES, E.F. SNSTEB SPE/ISRM 78163
accordance with the model presented by Lomba et al.9. The by Eq. (5)). The corresponding time for ion diffusion is about
values used here were = 0.17, 0.12 and 0.10 for KCl 60 hours.
concentrations of 5 w%, 10 w% and 20 w%, respectively. Fig. 6 shows the effect of using a 5 w% KCl solution for
The effect of fluid exposure clearly depends on the type of the drilling. This solution has the same chemical activity as the
ions involved. Fig. 3 shows the response to a CaCl2 solution in solution shown in Fig. 4; hence the differences between Fig. 4
comparison to a KCl solution for the same shale (but different and Fig. 6 are due to the ion exchange process only. It is seen
batches) as shown in Fig. 2. The samples shown in Fig. 3 had that the shrinkage induced by the ion exchange process relaxes
been exposed to sea water during the stablisation period prior the axial and tangential stresses in the vicinity of the borehole.
to the tests. Notice that the sample exposed to CaCl2 expands, This significantly improves the stability of the borehole with
while the sample exposed to KCl contracts. The total strain respect to possible shear failure, as it reduces the shear stress
amplitude estmated for the CaCl2 exposure was ch = -0.0012. (Fig. 7). The effect appears gradually, at about the same time
Other tests indicate that the expansion induced by CaCl2 as the positive effect of the chemical activity is declining (Fig.
exposure depends on the concentration of the solution. An 5). For practical applications, one should try to obtain the
apparent difference in saturation level for K+ ion exchange combined benefits of both effects, such that the chemical
between Fig. 2 and Fig. 3 may possibly be due to the activity ensures short term stability, while the ion exchange
difference in chemical content of the stabilizing fluid (3.5 w% takes over and ensures stability when the chemical activity
NaCl versus sea water). effect is declining.
The laboratory tests on small samples represent a fast and The borehole model has the ability to simulate the stress
inexpensive way to obtain information about how shales state in the vicinity of a borehole in a fairly general way, as it
respond to various mud compositions. In particular, they can also includes the general expressions (based on linear
tell us the magnitude of the response, and the characteristic elasticity) for stresses in the vicinity of a wellbore10.
time scale for the response process. This is essential Anisotropic horizontal stresses and deviating boreholes are
information for modelling the effect of mud composition on thus accounted for. Tests of the effective stresses versus
borehole stability. criteria for shear and tensile failure are also included, and the
model thus represents a powerful tool for analysis and
Numerical modeling of borehole stability prediction of borehole stability.
The effects of fluid exposure observed in the laboratory tests The model uses analytical solutions extensively, in order to
are likely to occur also in a borehole. To illustrate the reduce the computing time. Some of these solutions are exact
consequences for wellbore stability, a model based on the only at the borehole wall, and have reduced precision deeper
same equations and principles as the one used for the into the formation. This is not considered to be a significant
laboratory tests has been established for borehole geometry. problem for the model as a tool for borehole stability analysis,
Figs. 4 - 6 show examples of calculations using this model. as the critical events with respect to stability occurs in the
The data for the shale are representative for the material vicinity of the borehole wall.
described in the previous section. A simple geometry and
stress state has been chosen here in order to illustrate the Discussion
effects of mud chemistry: The well is vertical, the horizontal The behaviour of shales in the vicinity of a borehole is fairly
stresses are equal, and the well pressure is equal to the in situ complex, and several aspects have not been included in the
pore pressure. The shale properties are the same as derived modeling, or have only been treated approximately. Some of
from the laboratory tests described in the previous section. these are:
Data for the simulation are given in Table 1. - Anisotropy of chemical strain. Shales are highly anisotropic
Fig. 4 shows the situation when the drilling fluid is in on a microscopic scale, and to some extent the anisotropy
perfect balance with the pore fluid. There is an induced pore prevails also on a macroscopic scale. This may a.o. imply that
pressure effect, and no ion exchange effect. The stresses and the strain observed in the laboratory tests (Fig. 2) is
the fluid pressure remain constant in time. representative only for that specific orientation of the material.
Fig. 5 shows the effect of using a 10 w% NaCl solution for To compensate for this possibility, the chemically induced
the drilling fluid. This composition gives a chemical activity shrinkage has been scaled by a factor 1/3 in the borehole
(osmotic) effect, but no ion exchange effect. It is seen that the model. This would be the correct value if there is no strain in
pore pressure in the vicinity of the borehole is affected, and the lateral directions in the laboratory tests (not measured),
hence the effective stresses are also affected. The effect is and the shale downhole consist of a randomly oriented
positive for borehole stability, since it implies that the assembly of particles, each of which behaving like the sample
borehole wall gets some support from the fluid pressure in the observed in the laboratory.
well. This shifts the stress state away from the failure - Alteration of strength. It is possible (in fact evident in the
envelope, as illustrated in a Mohr-circle plot (Fig. 7). The data shown by Horsrud et al.5) that the shift in ion
positive effect declines gradually with time (as ions from the composition of the pore fluid, and the exchange of ions in the
drilling fluid move into the rock and reduces the contrast in shale microstructure, also have an impact on the strength of
chemical activity). With the actual data set, the characteristic the shale. Such an effect would of course have a significant
time for fluid pressure diffusion from the borehole wall to a impact on the stability of the borehole, as it would shift the
point 3 mm into the formation is about 20 minutes (as given position of the failure envelope (Fig. 7).
SPE/ISRM 78163 MUD CHEMISTRY EFFECTS ON TIME-DELAYED BOREHOLE STABILITY PROBLEMS IN SHALES 5
Constant load 6
Exposed 5 w% CaCl2
4
surface
Strain (mStr)
2
Piston
0
Measured
Sample -2 0 10 20 30 40 50
deformation
-4
Piston
-6
Time (h)
12
10
Strain (mStr)
5 w% KCl
8
Fig. 1 Fluid exposure test for small samples.
6
12 4
10 2
Strain (mStr)
8 5 w% KCl 0
6 0 10 20 30 40 50
4 Time (h)
2 Fig. 3 Effect of fluid exposure on a shale material, for two
0 different types of solutions. Solid line: mathematical model
0 20 40 60 80 including chemical activity effects and ion exchange.
Time (h)
3.5 w% NaCl
12
40
10 10 w% KCl
Strain (mStr)
8
6 35
Axial
Principal stresses [MPa]
4
Tangential
2 5 w% KCl 30
0
0 20 40 60 80
25
Time (h) Radial
12 20 w% KCl 20 Fluid
10 pressure
Strain (mStr)
8
15
6
0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7
4
Distance from borehole centre (relative to hole radius)
2 5 w% KCl
0 Fig. 4 Principal stresses and fluid pressure in the vicinity of the
borehole wall, as predicted by the borehole model, for a drilling
0 20 40 60 80 fluid containing 3.5 w% NaCl. Formation properties and in situ
Time (h) stresses are given in Tab. 1. Solid lines: 0.1 days after drilling;
Dashed lines: 1 day after drilling; Dotted lines: 10 days after
drilling. (The lines fall on top of each other.)
Fig. 2 Effect of fluid exposure on a shale material, for various
concentrations of KCl. Solid line: mathematical model including
chemical activity effects and ion exchange.
SPE/ISRM 78163 MUD CHEMISTRY EFFECTS ON TIME-DELAYED BOREHOLE STABILITY PROBLEMS IN SHALES 7
15 w% NaCl
40 Failure envelope
35
Shear stress
Axial
Principal stresses [MPa]
Initial state
Tangential
30
25
Radial
Pore pressure
reduction
20 Fluid
Shrinkage
pressure
Effective normal stress
15
0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 Fig. 7 Mohr-circle plot illustrating schematically the initial stress
state at the borehole wall, and the effects of pore pressure
Distance from borehole centre (relative to hole radius) reduction (induced by a contrast in chemical activity) and
shrinkage (induced by ion exchange).
Fig. 5 Principal stresses and fluid pressure in the vicinity of the
borehole wall, as predicted by the borehole model, for a drilling
fluid containing 15 w% NaCl. Formation properties and in situ
stresses are given in Tab. 1. Solid lines: 0.1 days after drilling;
Dashed lines: 1 day after drilling; Dotted lines: 10 days after
drilling.
5 w% KCl
40
35
Axial
Principal stresses [MPa]
Tangential
30
25
Radial
20 Fluid
pressure
15
0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7
Distance from borehole centre (relative to hole radius)