Near Well-Bore Condition (Skin) : Chapter 5 (Petroleum Production Systems by Economides)
Near Well-Bore Condition (Skin) : Chapter 5 (Petroleum Production Systems by Economides)
Near Well-Bore Condition (Skin) : Chapter 5 (Petroleum Production Systems by Economides)
qm Skin factor
Dps = s (dimensionless)
2p kh
Note the pressure drop changes with the flow rate and the viscosity for
the same perm alteration (skin factor). k is the undamaged permeability. 4
Skin Added to the Radial Flow Equation
• Skin or skin factor (s) is dimensionless and is
analogous to film coefficient in heat transfer
• Non dimensionalizing is done similar to PD function
done in well testing
• It can then be added to the radial flow equation as
shown below
8
Effective Wellbore Radius of Well with Skin
• With simple rearrangement of s to ln(es)
12
Bottom-hole pressures (with and w/o skin)
Reduced perm zone Unaltered perm zone
Dps
Dps = pwf ,ideal - pwf ,real
13
Hawkins’ (1956) Formula
• Calculates skin for a given depth of damaged zone
with reduced perm
• For pwf,ideal
• For pwf,real
14
Hawkins’ (1956) Formula
• Difference pwf,ideal – pwf,real = Dps, which is equal to
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 N bundle of capillaries of equal length but
different x-sectional area in unit bulk volume
• Suppose there are n, no. of capillaries with area
between A & A+dA
• The number of capillaries, n, with area between A and
A+dA is calculated as hdA where h is known as
distribution or density function.
22
Capillary Model for Porous Medium (cont.)
• If l is the pore length, lA will be the pore volume of a
capillary with area A
• Pore volume of n capillaries with area between A and
A+dA equals lA*n = lA*hdA
• Pore volume of all capillaries -> summing pore volume
of capillaries of area between 0 and dx, then dx and
2dx & so on until infinity
Hagen-Poiseuille
equation for flow in a
capillary/tube/pipe
24
Capillary Model for Porous Medium (cont.)
• Flow rate from tubes with area between A & A+dA =
c1vAhdA = c2A2hdA
• Total flow rate from all the capillaries is
¥
q = c2 ò A h dA
2
0
• Applying Darcy’s law to this bundle of capillary model
(Ab is the bulk area = some constant as Ab*l=1)
kAb
q= DP = ck
mL
25
Capillary Model for Porous Medium (cont.)
• Equating the two models and calculating k, we get
¥
k = c3 ò A h dA
2
26
27
28
Conservation Equation for Particle Transport in Porous Media
Continuity equation which states
Derived from continuity equation applicable to any conserved quantity (mass, energy, momentum, charge)
The above Convective-Diffusive equation is similar to Boltzmann transport (statistical behavior of thermodynamic
systems not in equilibrium) and Navier-Stokes equation (motion of viscous fluids)
Continuity equation
Quantity can be concentration, density, temperature
u = Darcy velocity
quantity/volume Flux = qt./area/time
c = Concentration of particles/solids (solid volume / fluid volume)
= Porosity
= Concentration of deposited particles (volume of deposited particles / bulk volume) Generation of
qt./vol/time
D = Dispersion Coefficient
29
Reasonable Approximations
• Incompressible flow (both fluids and solids)
• Dispersion is negligible
• Concentration of solids is assumed to be low
• Deposition follows a relation as proposed by Iwasaki (1937) ds
= luc
dt
30
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 )
31
Example Particle Concentration & Permeability
Profile Around the Well
32
Permeability Reduction Model
Pang and Sharma
2
é ù
ê 1+ s / (1- f ) ú
æ f 3 (1- fo )2 ö kds = ê o ú æ 1 ö
kdp = ç 3 ê s ú kdt = ç
2÷
è fo (1- f ) ø ê 1+
f
(d g
/ d )
p ú è 1+ bs ÷ø
ë (1- o
) û
Perm reduction Perm reduction due to Perm reduction due to
due to porosity increase in surface area increase in tortuosity
reduction
33
Filtration Coefficient () Using Trajectory Analysis (Rajagopalan and Tien,
1976)
3(1- f )
l= h Rajagopalan & Tien
(b,S)
Liquid 2d g (1976)
Shell (r,)
Grain
2rP 2 ( r P - r f )g dp
Pore
dg NG = NR =
throat
diameter b Limiting trajectory
9um dg
(dth)
(ap+ac , ) H ud g
N LO = N PE =
9pmrP 2u DBM
h is the efficiency of the grain, cin and cout are
35
the conc. of particles at inlet and outlet
Filtration Coefficient Model
(Compared to 106 Experiments)
13 Researchers
10000
Sakthivadivel Roque et al
Iwasaki (1995)
(1966) Exp
(1937) Exp:92-93
Exp:24-25
Exp:1-10 RT Model
Eliassen
1000 (1941 )
Exp:11 Fox, Cleasby Gruesbeck & Kau et al
(1966) Collins (1995)
Mackrle
Exp:26-27 (1982) Exp:94-97
(1960)
Exp:12-13 Ives Exp:88-91
100 (1967) Ison
Maroudas (1967)
Filtration Coefficient (1/m)
Exp:28-40
(1961) Exp:41-64
Exp:14-23 Fitzpatrick
(1973) Close et 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
36
Experiment No.
Formation of Internal & External Filter Cake
(A Multi-Component Filtration Model, Suri et al., SPEJ, March 2004)
Mud particles
Dh1, k1,dg1
Dh2, k2,dg2
Depth of
damage
Formation
grains
37
Simulator: UTDAMAGE (Suri et al.)
38
UTDAMAGE: Output Windows (Suri et al.)
39
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
40
Electrical Double Layer
41
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
42
43
Critical Salt Concentrations in Filtrate for
Minimum Damage (Schechter)
44
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
45
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
46
Inorganic Precipitation Example
47
Organic Precipitation Example
• Organic precipitation are waxes and asphaltenes
• Wax (paraffins) precipitates when temp. is reduced or
when oil composition changes due to gas liberation
• Asphaltenes (aromatic and napthenes) are colloidally
dispersed in the crude and are stable due to presence
of resins
• If resin conc. is reduced, asphaltenes can flocculate
causing damage
48
Fluid Damage: Emulsions, Rel Perm, Wettability
• Changes in the formation oil can cause
• Increase in apparent viscosity (emulsification)
• Decrease in effective perm (water block)
• Water-in-oil emulsions have
• Order of magnitude higher than oil viscosity
• Yield stress to be overcome to flow
• Typically caused by mechanical mixing dispersing one
phase into another but in formation
• Likely formed by surfactants or fines stabilizing small
droplets
49
Water Damage
• Increase in water saturation around the wellbore (also
called as water block) results in reduction in effective
perm to oil
50
Wettability Damage
• Certain chemicals can alter the wettability from
water-wet to oil-wet
• This would reduce the oil perm greatly
51
Perforation Damage (Kruger, 1986)
• Damage near perforations is unavoidable
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
52
Mechanical Damage
• Collapse of weak formation around wellbore
• In friable formations
• In formations weakened by acidizing
53
Biological Damage
• Particularly water injection wells are susceptible to
bacteria (particularly anaerobic that grows rapidly)
• Using these bacteria to plug the thief zones is studied
for EOR
• Bactericides should be used for prevention
54
Damage during Well Operations
• Drilling Damage
• Completion Damage
• Production Damage
• Injection Damage
55
Drilling Damage
• Invasion of drilling fluid particles and filtrate
• Particle damage is more severe typically
• Depth of invasion is typically an inch to a max of 1 ft
• Correctly sized particles invaded less deep
• Perm recovery is around 70-80% for correctly sized
drilling fluids
• Rule of thumb, 5 vol% particles should have diameter
> 1/3rd of mean pore size
56
Drilling Damage (cont.)
• Perforations or acids typically overcome drilling damage
• Filter cake grows after sudden sport loss
• Filter cake reduces the rate of filtrate loss, common filtrate
depths are 1-6 ft
• Filter cake is eroded by the shear force of drilling fluid
58
59
60
61
Drilling Damage (cont.)
• Filtrate damage can be reduced by tailoring ionic
composition compatible with formation water
• If water blocking or clay swelling can be a serious
problem water-based muds have to avoided
62
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
63
Underbalance needed to minimize perforation damage in
gas zones based on its perm (King et al., 1985)
64
Underbalance needed to minimize perforation damage in
oil zones based on its perm (King et al., 1985)
65
Stimulation Damage
• Stimulation fluids
• Similar solids and filtrate invasion as drilling and
completion fluids causing solids plugging and
precipitation
66
Production Damage
• Fines are mobilized due to
• high velocity (critical velocity, Schechter, 1992)
• If they are water-wet and when water production starts
67
Production Damage
• Precipitation of solids (similar to before)
• Inorganic minerals from brine
• Organic solids from crude
• Occasional stimulation can be done
• acids for carbonate precipitates
• solvents for waxes
• Chemicals can be squeezed for prevention
68
Injection Damage
• Suspended solids in injection fluids
• Typically suspended solids size is kept < 2 mm
• Precipitation due to incompatibility issues
• Injection of sulfate or carbonate when Ca2+, Mg2+, or Ba2+ is
present
• Cation exchange with clays can release divalent cations and
cause precipitation even when injection water is
compatible
• Growth of bacteria
• Bactericide should be added
69
Skin Components
• sd = damage skin
• sc+ = skin due to partial completion and slant
• sp = perforation skin
• spseudo = Rate dependent & phase skins
70
Rate Dependent Skin
• Well tests done at different rates can isolate non-rate dependent skin
71
Phase Skin
• Producing below bubble point point leads to gas evolution &
reduction in effective perm to oil
• When the well is shut-in the gas re-enters the solution
• Similarly gas retrograde condensate will lead to liquid drop around
the gas well causing effective perm to gas to reduce
• When the well is shut-in, or the pressure is increased, much of the liquid will
not re-enter the gas (Fussell, 1973; Cvetkovic et al., 1990)
• Natural gas needs to be injected that may dissolve the condensate (Huff &
Puff operation can be repeated periodically)
• Perforating, acidizing and fracking will alter most pre-treatment skins
72
Skin from Partial Completion (sc)
• Frequently, wells are partially completed
(open well height < total reservoir height)
• If completed interval > 75% of total, sc is
negligible
• Also referred to as partial penetration
• Can be a result of a bad perf job
• Can be deliberate to retard or avoid coning
• Generates early time spherical flow which
allows calculating vertical permeability
73
Skin from Slant (s)
• While partial completion reduces the area of
contact (hence +skin), slant well, increases the
area of contact (hence –skin)
• Combined skin is sc+
• Cinco-Ley et al. (1975) calculated sc+ semi-
analytically & presented tables for various
combinations: h/rw, zw/h, hw/h,
• hD = h/rw is dimensionless reservoir thickness
• Elevation ratio (zw/h), zw is elevation of mid of
perfs & completion ratio (hw/h)
74
75
small
76
Higher
than for
hD=100
(sc~21)
77
78
Perforating
• Perforating guns are attached to wireline, coiled
tubing, or tubing
• Shaped charges arranged helically for good perf
density and small angle between adjoining perfs.
(phasing)
• Cable head connects string to wireline/tubing and
provides a weak point to break if problem arises
• Corelation device ensures right depth by locating
casing collars and matching with previous logs
• Positioning device orients shots for optimum
geometry
79
Perforating (cont.)
• Shaped charges consist of a case, explosive and a
liner
• Electric current initiates the detonation
• Perforations with diameter 0.25 – 0.4 in. and a
tunnel between 6-12 inches typically created
• Usually done underbalance for immediate
flowback carrying debris resulting in a cleaner perf
• Dimensions, number, and phasing have a
controlling role in well performance
80
Perforation Skin (sp)
• Karakas & Tariq (1988) calculated
perforation skin semi-analytically
• Divided the skin into components
• Constant a1, a2, b1, b2 are also functions of phasing as given in the
table, sV is the potentially the largest contributor, for ex. for low spf
84
Wellbore Effect (Calculation of swb)
85
86
87
88
89
90
Near Well Damage and Perforation Skin
Combined
• Karakas and Tariq (1988) also showed that near well damage & perf.
can be combined
91
Near Well Damage and Perforation Skin
Combined (cont.)
• If perforations terminate outside the damage zone
92
Horizontal Well Production Model (Joshi)
Excellent producers for reservoirs with
thickness < 50 ft or with good vertical
permeability
Also needed in tight reservoirs and shale
reservoirs for larger access & fracking
94
Horizontal Well Damage Shape (Frick &
Economides, 1991)
• Described the shape of damage along and normal to a horizontal well
Damage is highest near the
aH,max is the maximum damage depth along
vertical section due to higher
the horizontal near the vertical section cone
residence time and almost
zero at the tip
98
99
100
101
102
103
104
105
106
107
108
109