Near Well-Bore Condition (Skin) : Dr. Ajay Suri Associate Professor Dept. of Petroleum Engineering IIT (ISM) Dhanbad
Near Well-Bore Condition (Skin) : Dr. Ajay Suri Associate Professor Dept. of Petroleum Engineering IIT (ISM) Dhanbad
Near Well-Bore Condition (Skin) : Dr. Ajay Suri Associate Professor Dept. of Petroleum Engineering IIT (ISM) Dhanbad
8
Effective Wellbore Radius of Well with Skin
12
Bottom-hole pressures (with and w/o skin)
Reduced perm zone Unaltered perm zone
Dps
D ps =pwf ,ideal - pwf ,real
13
Hawkins’ (1956) Formula
• For pwf,real
14
Hawkins’ (1956) Formula
17
Plugging due to Solids
• Most common
• Sources
– Solids from the wellbore fluids (mostly during
drilling) and some could come during completion,
although completion fluids are typically solids
free)
– Dispersion of clays present in the rock
– Precipitation of minerals in solution
– Growth of bacteria aggregates
18
SEM of a Fontainebleau sandstone rock sample. Institut
Francais du Petrole. Grain 200 mm, surface roughness 1 mm
19
20
Particle Entrapment &
Perm Reduction
(Schechter, 1992)
• Surface deposition of
particles
• Reduced porosity
• Increased surface area
• Increases tortuosity
• Internal Pore blockage
• External filter cake
21
Capillary Model for Porous Medium
• Imagine bundle of capillaries of equal length
but different x-sectional area
• Probability of capillaries is the relative
abundance or no. of capillaries with area
between A & A+dA to the total no. of capillaries
• Probability of capillaries in unit volume (bulk)
having area between A and A+dA is hdA where
h is known as probability distribution function
or pore density function
22
Capillary Model for Porous Medium (cont.)
• If l is the pore length, lA is the pore volume of
a single capillary with area A
• Note the probability of pores with area exactly
equal to A is zero
• Area of pores with area between A and A+dA
in unit volume equals A*no. of pores = A*hdA
¥
f =l òAh dA
0
• For k, we know DP is same across all tubes
23
Capillary Model for Porous Medium (cont.)
• Velocity in a tube is proportional to it’s area
• Flow rate in a tube is proportional to it’s area
times velocity = vA
• No. of tubes in unit volume with area between
A & A+dA is hdA
• Flow rate from tubes with area between A &
A+dA = vAhdA = cA2hdA
• Total flow rate per unit volume and k will be
¥ ¥
proportional
q =c1 òA2h dA k =c2 òA2h dA
24
0 0
DP same in all capillaries
v a A in a capillary
Q = Av of all capillaries
25
26
Steady State Conservation Equation for Particle
Transport in Porous Media
u = Darcy velocity
c = Concentration of solids (solid vol. / fluid vol.)
= Porosity
= Fraction of deposited particles (vol. of deposited particles / bulk vol.)
D = Dispersion Coefficient
27
Reasonable Approximations
• Incompressible flow (both fluids and solids)
• Dispersion is negligible
• Concentration of solids is assumed to be low
• Deposition follows an empirical relation as proposed by
Iwasaki (1937)
ds
=l uc
dt
28
Filtration Coefficient () Using Trajectory Analysis
(Rajagopalan and Tien, 1976)
29
Happel’s Sphere-in-cell Model (1958)
1.2
h =(1- f )2/3 AS N LO1/8 N R15/8 + 3.375 ´ 10- 3 (1- f)2/3 AS N G N R - 0.4 + 4 AS 1/3 N PE - 2/3
Flow
3(1- f ) Rajagopalan & Tien (1976)
l= h
2d g
(b,S)
Liquid
Shell (r,)
Grain
2rP 2 ( r P - r f )g dp
Pore
dg NG = NR =
throat
diameter b Limiting trajectory 9um dg
(dth)
H ud g
(ap+ac , ) N LO = N PE =
9pm rP 2u DBM
30
Filtration Coefficient Model (Compared to 106 Experiments)
13 Researchers
10000
S akthivadive l Ro que e t al
Iwas aki
(1966) Exp (1995)
(1937) Exp:92-93
Exp:24-25
Exp:1-10 RT Model
Elias s e n
1000 (1941 )
Exp:11 Fo x, Cle as by Grue s be c k & Kau e t al
(1966) Co llins (1995)
Mac krle
Exp:26-27 (1982) Exp:94-97
(1960)
Exp:12-13 Ive s Exp:88-91
100 (1967) Is o n
Maro udas (1967)
Filtratio n Co e ffic ie nt (1/m)
Exp:28-40
(1961) Exp:41-64
Exp:14-23 Fitzpatric k
(1973) Clo s e e t al
10 Exp:65-87 (2005)
Exp:98-106
0.1
0.01
0.001
0 10 20 30 40 50 60 70 80 90 100 110 31
Expe rime nt No .
Concentration of Suspended Particles in the
Formation (Linear Geometry)
fx
c(x,t) =0 t<
u
fx
c(x,t) =cin exp(- l x) t>
u
( x , t ) 0 ( x, t )
32
Example Particle Concentration & Permeability
Profile Around the Well
33
Permeability Reduction Model
Pang and Sharma
k / ko =kdp k ds kdt é ù
2
ê 1+ s / (1- f ) ú
æf 3 (1- fo )2 ö k ds =ê o ú æ 1 ö
kdp =ç 3 ê s ú kdt =ç
2÷ 1+ (d / d ) è 1+ bs ÷
f
è o (1- f ) ø ê (1- f ) g
ë o
p ú
û
ø
Perm reduction due to Perm reduction due to increase Perm reduction due to
porosity reduction in surface area increase in tortuosity
34
Formation of Internal & External Filter Cake
(A Multi-Component Filtration Model, Suri et al., SPEJ, March 2004)
Mud particles
h1, k1,dg1
h2, k2,dg2
Depthof
damage
Formation
grains
35
Simulator: UTDAMAGE (Suri et al.)
UTDAMAGE: Output Windows (Suri et al.)
Fines Mobilization
• Due to change in chemical composition of water
(reduced salinity or ions between the formation
water and filtrate)
• Due to shear forces applied by moving fluid
• Sources of filtrate
– Drilling fluids
– Completion fluids
– Stimulation fluids
– Injection fluids
• Essential to check filtrate compatibility with the
formation 38
Dispersion of Clay Particles
• Sudden decrease in salinity in sandstone may
cause dispersion of clay particles (called water
sensitivity)
• Depends on cation type, pH, rate of salinity
change
39
Critical Salt Concentrations in Filtrate for
Minimum Damage (Schechter)
40
Prevent of Clay Dispersion
• A critical salt concentration is typically
recommended
• A minimal conc. of a monovalent ions for ex.
as given for Berea SS and a sufficient fraction
of divalent ions should be present
• A common criteria is to have 2 wt % of KCl and
at least 1/10th of salt is should be divalent
cations
41
Chemical Precipitation
• Precipitation of solids from brine or crude can
cause severe plugging
• Triggers are change in temp., pressure, or
composition alteration of the phases
• Precipitates can be organic or inorganic
• Inorganic precipitates are usually divalent ions
such as Ca2+, Ba2+ combined with carbonate or
sulphate ions
42
Inorganic Precipitation Example
46
Wettability Damage
47
Perforation Damage (Kruger, 1986)
Lab testing of
perforating into
sandstone cores
showed
damaged zone to
be 1/4-1/2 inch
thick with perm
of 7-20% of
undamaged
perm
48
Mechanical Damage
49
Biological Damage
50
Damage during Well Operations
• Drilling Damage
• Completion Damage
• Production Damage
• Injection Damage
51
Drilling Damage
54
Completion Damage
• Completion fluids
– Cements
– Perforating fluids
– Stimulation fluids
• Similar solids and filtrate invasion as drilling
fluid
• Typically solids should be < 2 ppm with size < 2
mm
• Cements have high Ca2+, so potential for
precipitation 55
Underbalance needed to minimize perforation damage in
gas zones based on its perm (King et al., 1985)
56
Underbalance needed to minimize perforation damage in oil
zones based on its perm (King et al., 1985)
57
Stimulation Damage
• Stimulation fluids
• Similar solids and filtrate invasion as drilling
and completion fluids causing solids plugging
and precipitation
58
Production Damage
59
Production Damage
• Precipitation of solids
– high velocity (critical velocity, Schechter, 1992)
– If they are water-wet and when water production
starts
60
Production Damage
61
Injection Damage
• sd = damage skin
• sc+q = skin due to partial completion and slant
• sp = perforation skin
• spseudo = Rate dependent & phase skins
63
Rate Dependent Skin
• Well tests done at different rates can isolate
non-rate dependent skin
At high rates
rate dependent
skin Dq can be
much larger
than s
64
Phase Skin
• Producing below bubble point point leads to
gas evolution & reduction in effective perm to
oil
• Similarly gas retrograde condensate will lead
to liquid drop around the gas well causing
effective perm to gas to reduce
65
66
67
68
69
70
71