Module 8 Relative Permeability
Module 8 Relative Permeability
Module 8 Relative Permeability
Synopsis
What is water-oil relative permeability and why does it matter?
endpoints and curves, fractional flow, what curve shapes mean
How do we measure it (in the lab)? How do we quality control and refine data?
Page 2
Applications
To predict movement of fluid in the reservoir
e.g velocity of water and oil fronts
To predict and bound ultimate recovery factor Application depends on reservoir type
gas-oil water-oil gas-water
Page 3
Definitions
Absolute Permeability
permeability at 100% saturation of single fluid
e.g. brine permeability, gas permeability
Effective Permeability
permeability to one phase when 2 or more phases present
e.g. ko(eff) at Swi
Relative Permeability
ratio of effective permeability to a base (often absolute) permeability
e.g. ko/ka or ko/ko at Swi
Page 4
Requirements
Gas-Oil Relative Permeability (kg-ko)
solution gas drive gas cap drive
Jargon Buster!
Relative permeability curves are known as rel perms Endpoints are the (4) points at the ends of the curves The displacing phase is always first, i.e.:
kw-ko is water(w) displacing oil (o) kg-ko is gas (g) displacing oil (o) kg-kw is gas displacing water
Page 6
Sro = 0.25
Kw = 24 mD
Endpoints
1.0 0.9 0.8 0.7 0.6
Endpoint- oil kro = ko/ko @ Swir = 80/80 =1 Swir = 0.20 Sro = 0.25
0.0 0.0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7 0.8 0.9 1.0
Page 8
Endpoints
1.0 0.9 0.8 0.7 0.6
Swir = 0.20
0.5 0.4 0.3 0.2 0.1 0.0 0.0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7
Sro = 0.25
0.8
0.9
1.0
Page 9
Curves - 1
1.0 0.9 0.8 0.7 0.6
Swir = 0.20
0.5 0.4 0.3 0.2 0.1 0.0 0.0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7
Sro = 0.25
0.8
0.9
1.0
Page 10
Curves - 2
1.0 0.9 0.8 0.7 0.6
Swir = 0.20
0.5 0.4 0.3 0.2 0.1 0.0 0.0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7
Sro = 0.25
0.8
0.9
1.0
Page 11
Curves - 3
1.0 0.9 0.8 0.7 0.6
Swir = 0.20
0.5 0.4 0.3 0.2 0.1 0.0 0.0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7
Sro = 0.25
0.8
0.9
1.0
Page 12
Relative Permeability
1
0.9
0.8
0.7
0.6
0.5
kro krw
0.4
viscous forces
Darcys Law
0.3
0.2
0.1
capillary forces
low flood rates
Page 13
Waterflood Interpretation
Welge
fw=1
fw only after BT
Average Saturation behind flood front
Swf , fw | S
wf
fw
Sw at BT
fw =
1 +
k ro . k rw
w o
Page 15
Swc
Sw
1-Sor
fw =
1 +
k ro . k rw
k rw o . M= k ro w
M< 1: piston-like M > 1: unstable
w o
Page 16
dS wa = fo2 dW
d( 1 ) f WI r = o2 1 k ro 2 d( ) W
pt =0 Ir = pt =i
Page 17
Page 18
Leverett
capillary boundary effects on short cores boundary effects negligible in reservoir
Page 19
End Effect
Pressure Trace for Flood zero p (no injection) start of injection water nears exit
p increases abruptly until Sw(exit) = 1-Sro and Pc nears zero suppresses krw Sw(exit) = 1-Sro, Pc ~0 rate of p increase reduces as krw increases
BT
Page 20
After BT
Scaling Coefficient
Breakthrough Recovery (Rappaport & Leas) Affected by Pc end effects At lengths > 25 cm Little effect on BT recovery (LVw > 1) Hence composite samples or high rates
Page 21
Overcome by:
flooding at high rate
300 ml/hour +
Page 23
S wn
Page 24
S w S wi = 1 S wi S ro
Normalisation
Swn = 1
1 0.9 0.8
0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7 0.8 0.9 1 krw at Sro krwn = 1
Sample 1 Sample 2
krwn = 1
Page 25
Corey Exponents
Depend on wettability
Wettability Water-Wet Intermediate Wet Oil-Wet No (kro) 2 to 4 3 to 6 6 to 8 Nw (krw) 5 to 8 3 to 5 2 to 3
Uses:
interpolate & extrapolate data lab data quality control
Page 26
Gas is non wetting takes easiest flow path kro drops rapidly as Sg increases krg higher than krw Srog > Srow in lab tests
end effects
Sgc ~ 2% - 6%
1-(Srog+Swi)
kro krg
Sgc
Page 28
0.1
1-(Srog+Swi)
kro krg
0.01
0.001 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 Gas Saturation (fractional)
Page 29
Gas-Oil Curves
Most lab data are artefacts
due to capillary end effects
Tests should be carried out on long cores
krg = 1 at Srog = 0
Page 30
kro = Son No
Son = 1 Sg Swir Srog 1 Swir Srog
krg = Sgn Ng
Sgn = Sg Sgc 1 Swir Srog Sgc
0.1
0.01
0.001
Sgc = 0.03
0.0001
0.00001 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 Gas Saturation (-)
Page 32
0.01
0.001
Composite Gas-Oil Curves Ng : No : Sgc: Srog: krg' : 2.3 4.0 0.03 0.10 1.0
Ng = 2.3; Swir = 0.15 Ng = 2.3; Swir = 0.20 11a-5 # 4 11a-5 # 31 11a-5 # 34 11a-5 #39 11a-7 BEA5 11a-7 BEA7 11a-7 BEB5 11a-7 BEC5
0.0001
0.00001 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 Swi+Sg (fraction)
Page 33
Laboratory Methods
Core Selection
all significant reservoir flow units often constrained by preserved core availability core CT scanning to select plugs
Core Size
at least 25 cm long to overcome end effects butt samples (but several end effects?) flood at high rate to overcome end effects?
Page 34
Test States
Fresh or Preserved State
tested as is (no cleaning) probably too oil wet (e.g OBM, long term storage) Native state term also used (defines bland mud) Some labs fresh state is other labs restored state
Cleaned State
Cleaned (soxhlet or miscible flush) water-wet by definition (but could be oil-wet!!!!!!)
Test State
Fresh-State Tests
too oil wet data unreliable data unreliable
Cleaned-State Tests
too water wet (or oil-wet)
Restored-State Tests
Page 36
native wettability restored data reliable (?) if GOR low can use dead crude ageing (cheaper) if GOR high must use live crude ageing (expensive) if wettability restored - use synthetic fluids at ambient ensure cores water-wet prior to restoration
Centrifuge
faster than others Swir can be non-uniform
Porous Plate
slow, grain loss, loss of capillary contact
Page 37
Swir uniform
25
20 Swi (%)
180 psi ???
15
10
200 psi
Page 38
Centrifuge Tests
Displaced phase relative permeability only
oil-displacing-brine : krw drainage brine-displacing-oil : kro imbibition assume no hysteresis for krw imbibition
oil-wet or neutral wet rocks?
1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 0.0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7 0.8 0.9 1.0
Good for low kro data (near Sro) Computer simulation used Problems
uncontrolled imbibition at Swirr mobilisation of trapped oil sample fracturing
Page 39
Relative Permeability (-)
Test Conditions
fresh state cleaned state restored state ambient or reservoir conditions
Page 40
Unsteady-State Waterflood
Saturate in brine Desaturate to Swirr Oil permeability at Swirr (Darcy analysis) Waterflood (matched viscosity)
o w o = w res lab
Incremental oil recovery measured kw at Srow (Darcy analysis) Relative permeability (JBN Analysis)
Page 42
Unsteady-State Procedures
Water Oil Only oil produced Measure oil volume
Page 43
Unsteady-State
Rel perm calculations require
fractional flow data at core outlet (JBN) pressure data versus water injected
o/w = 30:1 Unstable flood front Early BT Prolonged 2 phase flow Oil recovery lower o/w = 3:1 Stable flood front BT delayed Suppressed 2 phase flow Oil recovery higher
0.2 0.1 0.0 0.0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7 0.8 0.9 1.0
Page 45
Unsteady-State Tests
Only post BT data are used for rel perm calculations
Sw range restricted if matched viscosities
Advantages
appropriate Buckley-Leverett shock-front reservoir flow rates possible fast and low throughput (fines)
Disadvantages
inlet and outlet boundary effects at lower rates complex interpretation
Page 46
Steady-State Tests
Intermediate relative permeability curves
Saturate in brine Desaturate to Swir Oil permeability at Swir (Darcy analysis) Inject oil and water simultaneously in steps Determine So and Sw at steady state conditions kw at Srow (Darcy analysis) Relative Permeability (Darcy Analysis)
Page 47
Oil in
Water in
Mixing Sections
Coreholder
Page 48
Steady-State Procedures
Summary
100% Oil:
ko & kw at Sw(n)
Page 49
Easier analysis
Darcy vs JBN
Slower
days versus hours
Laboratory Tests
You can choose from:
matched or high oil-water viscosity ratio cleaned state, fresh state, restored-state tests ambient or reservoir condition high rate or low rate USS versus SS
Oil Recovery
70 Fixed - 120 ml/hour 60 Oil Recovery (% OIIP) Preferred 360
50
40
120
30
120
Page 52
Steady-State
kg-ko, kg-kw and kw-kg saturation determination difficult much slower
Page 54
Imbibition Tests
Waterflood
low rate waterflood from Swi to Sgr obtain krg and krw on imbibition Sgr too low (viscous force dominates)
129.90 g
Page 56
Page 57
Page 58
Imbibition Kw-Kg
1
Drainage
kr
Imbibition
krw 0
Page 59
Sw
Page 60
Wettability
Page 61
Wettability
Page 62
Wettability
Waterflood of Water-Wet Rock
Page 63
front moves at uniform rate oil displaced into larger pores and produced water moves along pore walls oil trapped at centre of large pores - snap-off BT delayed oil production essentially complete at BT water invades smaller pores earlier BT oil remains continuous oil produced at low rate after BT krw higher - fewer water channels blocked by oil
Effects of Wettability
Water-Wet
better kro lower krw krw = kro > 50% better flood performance poorer kro higher krw kro = krw < 50% poorer flood performance
Oil-Wet
Page 64
Preserved Core Neutral to oil-wet low kro - high krw Extracted Core Water wet high kro - low krw
Page 65
Intermediate Wet
No = 4 Nw = 4 Swir = 0.15 Sro = 0.25, krw = 0.5, ultimate recovery = 0.706 OIIP
Oil Wet
No = 8 Nw = 2 Swir = 0.10 Sro = 0.20, krw = 0.75, ultimate recovery = 0.778 OIIP
Page 66
o/w = 3:1
Page 67
Page 68
Page 69
Page 70
WW fw IW fw
Page 71
Page 72
Rock Texture
Page 74
Viscosity Ratio
krw and kro - no effect ? End-Points - viscosity dependent Hence: use high viscosity ratio for curves use matched for end-points
Page 75
Saturation History
100 %
Primary Drainage
Primary Imbibition
NW
kr
Swi W
kr
Sro
W 0% 0%
Page 76
0%
Sw
100 %
0%
Sw
100 %
Flow Rate
Reservoir Frontal Advance Rate
about 1 ft/day
Flow Parameters
End Effect Capillary Number Flood Capillary Number
Nc end
k o vL
Ncend 2.3 0.07 0.02 0.02 0
Nc =
Rate (ml/h) 4 120 360 400 Reservoir
Nc
For reservoir-appropriate data Nclab ~ Ncreservoir If Ncend > 0.1 kro and krw decrease as Ncend increases
Page 78
Bump Flood
1.0 0.9 0.8 0.7 0.6 0.5 High Rate krw ??? 0.4 0.3 0.2 0.1 0.0 0.0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7 0.8 0.9 1.0 Low Rate krw' Bump Flood krw'
Page 79
Considerations
ensure Swi is representative low rate floods for Sro: bump for krw steady-state tests
Page 80
Considerations
high rate floods (minimum Dp = 50 psid) to minimise end effects steady-state tests with ISSM low rates with ISSM and simulation
Page 81
Recommendations
data acquired at representative rates (e.g. near wellbore, grid block rates)
Page 82
JBN Validity
High Viscosity Ratio
viscous fingering invalidates 1D flow assumption
Low Rate
end effects invalidate JBN
Test Recommendations
Wettability Conditioning
flood rate selected on basis of wettability Amott and USBM tests required Wettability pre-study
reservoir wettability? fresh-state, cleaned-state, restored-state wettabilities
beware fresh-state tests (often waste of time) reservoir condition tests most representative
but expensive and difficult
Page 84
Wettability Restoration
Hot soxhlet does not make cores water wet! Restored-state cores too oil wet Lose 10% OIIP potential recovery
USBM
0.0 1.0 STRONGLY WATER-WET
Amott
Page 85
Viscosity Ratio
matched viscosity ratio for end-points investigate viscosity dependency for rel perms normalise then denormalise to matched end-points
Page 86
Saturation Determination
conventional
grain loss, flow processes unknown
NISM
can reveal heterogeneity, end effects, etc
Page 87
Use of NISM
Examples from North Sea Core Laboratories SMAX System
low rate waterflood followed by bump flood X-ray scanning along length of core end-points some plugs scanned during waterflood
Fresh-State Tests
core drilled with oil-based mud
Page 88
X-Ray Scanner
Coreholder (invisible to Xrays)
X-rays detected X-rays emitted
Scanning Bed
X-ray adsorption
0 %
Sw(NaI)
100%
Page 90
bump flood produces oil from body of core neutral wet data reliable
Page 92
oil-wet end effect bump produces incremental oil from body of core but does not remove end effect neutral to oil-wet
Page 93
data unreliable
Page 94
Page 95
Simulation required
e.g. SENDRA, SCORES
Page 96
Fluid properties
viscosity, IFT, density
Imbibition Pc curve (option) ISSM or NISM Scans (option) Beware several non-unique solutions possible
Page 97
History Matching
Pressure and production
1.66 cc/min
800 Differential Pressure (kPa) 700 600 500 400 300 200 100 0 0,1 1,0 10,0 100,0 Time (min) 1000,0 Measured differential pressure Simulated differential pressure Measured oil production Simulated oil production 1,0 2,0 4,0 Oil Production (cc) 6,0
5,0
3,0
0,0 10000,0
Page 98
History Matching
Saturation profiles
0.8 0.7
Water Saturation
0.6 0.5 0.4 0.3 0.2 0.0 0.2 0.4 0.6 0.8 1.0
Page 99
Relative Permeability
0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Krw Kro low rate end point high rate end point
Water saturation
Page 100
Relative Permeability
0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Krw Kro low rate end point high rate end point Krw Simulation Kro Simulation
Water saturation
Page 101
Quality Control
Most abused measurement in core analysis Wide and unacceptable laboratory variation Quality Control essential
test design detailed test specifications and milestones contractor supervision modify test programme if required
Benefits
better data more cost effective
Page 102
Curve Shapes
1 0.9 0.8 0.7 0.6
Kr
1
0 0.2 0.4 Sw 0.6 0.8 1
0.1
Kr
0.01
Kro Krw
0.001
Page 104
Sro Determination
Compute Son high, medium and low Sro
1
curves down Sro too high curves up Sro just right straight line
Page 105
1.000 0.100
Son = (1-Sw-Sor)/(1-Swi-Sor)
0.0001 0.010
Refine krw
Refined krw
Use refined Sro Plot krw versus Swn Fit line to last few points
Krw 0.1 1
Page 106
No Nw
Plot vs Sw Take No and Nw from flat sections Least influenced by end effects
0.2
0.4 Sw
0.6
0.8
Page 107
No
Nw
Normalisation Equations
Water-Oil Data
k ro n =
k ro k ro end
k ro n =
k ro k ro end
krgn =
krg krgend
Page 109
1 0.9 0.8
0.7 0.6 0.5 Swirr Swn = 0 0.4 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7 0.8 0.9 1 Sw = 1-Sro Swn = 1
Sample 1 Sample 2
Page 110
1 0.9 0.8
0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 Water Saturation (-) 0.6 0.7 0.8 0.9 1 krw at Sro krwn = 1
Sample 1 Sample 2
Page 111
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Page 112
Page 113
Denormalisation
Group data by zone, HU, lithology etc Determine Swir (e.g. logs, saturation-height model) Determine ultimate Sro
e.g. from centrifuge core tests
Denormalisation Equations
Water Oil
Denormalised Endpoints Water-Oil Swi kro (@Swi) krw (@1-Srow) From correlations & average data