2-Dynamic Simulation Application On Wells Sept 07-2
2-Dynamic Simulation Application On Wells Sept 07-2
1
Dynamic Simulation
Separate
Closure laws
continuity
are semi-
equations for
the Gas, Liquid
mechanistic
(O+W) bulk and
Boundary Initial and required
Liquid droplets Conditions Conditions experimental
coupled verification
through
interphasial A mixture energy
mass transfer. conservation
equation is applied. Dynamic
Numerical solution scheme: semi-
Simulator implicit integration method –
Two momentum equations are used:
allows for relatively long time
1) a combined one for the Gas and
possible Liquid droplets,
steps with efficient run times
2) a separate one for the Liquid film
2
OLGA - Dynamic 3-phase Flow Simulator
the measurements 2o
Complex
sets
Research Oil & Gas
Scandpower
• 15 lower GOR pipeline data sets Institutes Industry
3
1 D - Well Dynamic Simulation
0
-500
-1000
-1500
-2000
Conduction
2
TFluid 4
3
1
Te
4
OLGA® for Wells Thermal Model
Warm-up of Well and Formation
Temperature
t0 t1 t2 t3 tss
Fo
rm
a
Depth
Wellbore
tio
nT
em
p.
Temperature
155
150
C
145
140
135
130
0 500 1000 1500 2000 2500 3000 3500 4000
Length [m]
10
5
OLGA® for Wells Thermal Model
Warm-up of Well and Formation
11000
10000
9000
8000
7000
sia
6000
p
5000
4000
3000
2000
1000
0
0 5 10 15 20 25
Time [h]
12
6
Transient vs. Steady State Solutions
13
SEPARATED
DISTRIBUTED
7
Potential Problems for Stable Multiphase Flow
BUBBLE BUBBLE
BUBBLE
SLUG FLOW
45 bar
20 bar ANNULAR
Down
Horiz. STRATIFIED
Up STRATIFIED
16
8
Multiphase Flow Modeling
17
9
Multiphase Flow Modeling
SPE24789
142
140
138
136
Temperature (deg C)
130
128
126
124
0 50 100 150 200 250 300 350 400
Gas Rate (MMscfd)
10
Well Integrity – WHT-WHP Calculations
141
140
Wellhead Temperature (deg C)
139
138
137
136
135
134
133
132
3500 3700 3900 4100 4300 4500 4700 4900 5100 5300 5500
Wellhead Pressure (psia)
140
135
130
6500 Psia
Temperature (deg C)
6000 Psia
125 5500 Psia
5000 Psia
4500 Psia
120 4000 Psia
3500 Psia
3000 Psia
115
110
105
100
0 50 100 150 200 250 300 350 400
Gas Rate (MMscfd)
11
Transient vs. Steady State Solutions
23
OIL
Poettman and Carpenter 1952
Gilbert 1954
Duns & Ros 1961 Higher rate wells
1962
Winkler & Smith
1965 Slugging wells
Hagedorn & Brown 1967
Orkiszewski 1969
Wallis 1972
Aziz & Govier 1973 0
Flowlines (wells > 70 )
Beggs & BrilI 1974
Chierici et al MMS (Moreland 1978 All wells
Mobil/Shell)
(A modified Duns & Ros Correlation)
DRY GAS 3 6 3
(LGR < 115 m /10 m ) 20 bbl/MMscf
Poettman & Carpenter 1952
Cullender and Smith Equation 1956 No liquid hold-up
1964
Dukler (Non-Slip) No liquid hold-up
Lacroix Equation
(Average temp. and Comp. Method)
Eaton and Dukler 1969
Aziz and Govier 1972
Gray (API 14B) 1974 All wells
Ros/Gray 1978 All Wells
WET GAS 3 6 3
(LGR > (115 m /10 m ) 20 bbl / MMscf
3 6 3
(LGR > (550 m /10 m (100 b/MMScf])
Modified Cullender and Smith Method
1952
Hagedorn and Brown 1965 Slugging wells
Aziz and Govier 1972
Beggs and Brill 1973 Flowlines
Govier and Fogarasi 1975
Gray (API 14B) 1974 All wells
Ros/Gray 1978 All wells
LGR: Liquid Gas Ratio
12
Transient vs. Steady State Solutions
0.70
0.60
PREDICTED LIQUID HOLDUP
0.50
OLGAS
0.40 Beggs & Brill
0.30
0.20
0.10
0.00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80
MEASURED LIQUID HOLDUP
Use of Correlations
25
26
13
Well Flow Dynamics
Production / Injection Wells
• Flow Stability / Minimum stable flow rates / Slug Mitigation
– Tubing sizing / Production Optimisation
• Water accumulation / Corrosion rates
– Material selection / contact time of water as film
• Flow assurance, Wax , Hydrates
– Location of SCSSV
– MeOH / Glycol requirements
14
Well Dynamic Simulation
Flow Stability Solutions
29
15
Well Dynamic Simulation
Water Acummulation / Corrosion Rates Solutions
Corrosion Modelling
32
16
Well Dynamic Simulation
Water Acummulation / Corrosion Rates Solutions
33
17
Well Dynamic Simulation
Flow Assurance / Hydrates, Wax Solutions
assurance problems? 10
DTHYD (°C)
0.5 hr
1 hr
-20
– Remove supply of water -30
2 hr
3 hr
5 hr
7 hr
9 hr
Hydrate Temperature
• Hydrate dissociation Temperature
10
0
0 hr
5 min
-10 15 min
DTHYD (°C)
0.5 hr
1 hr
-20 2 hr
3 hr
-30 5 hr
7 hr
9 hr
-40 12 hr
15 hr
30 hr
-50
-60
0 5000 10000 15000 20000 25000
Distance from wellheads (m)
Hydrate formation when re-starting a cooled down line
36
18
Hydrate control
• Insulation (passive thermal control)
– Used for tie-ins and short to medium pipelines
– Not usual for long gas-condensate lines
• Bundles (active thermal control)
– Complex bundles used for deep off-shore.
– Generally used for risers, tie-ins and short to
medium pipelines
• Inhibition
– Typical for gas-condensate systems
– Used for oil systems at critical points
(e.g well-heads) and during critical operational
phases shut-in-cool-down-start-up).
37
Hydrate curves
200
180
160
140
120
P (bara)
100 20 % MeOH
80
60
20 % MEG
40
20
0
0 5 10 15 20 25 30
T (C)
38
19
MEG Tracking Modelling
– MEG concentration
200
along well / flowline
– Time to reach
P (bara) 150
a desired 100 0 MEG
MEG concentration 50
10 % MEG
30 % MEG
40 % MEG
0
0 10 20 30
39
T (C)
40
20
Wax - Definitions
• Slightly soluble
• Pour Point
– The lowest temperature at which oil can be handled
without excessive amounts of wax crystals forming out of
solution.
41
• Wax Precipitation
• Molecular Diffusion and Shear Transport
• Wax Melting
42
21
What can be done with the OLGA Wax Module?
43
WAT
Dissolved wax Pour Point
44
22
Wax Modeling Concepts: WAT
WAT
Pour Point
•Wax crystals start appearing
•Uncertainty of measurements of WAT (Wax
Appearance Temperature) is within +/- 5 F
45
WAT
Pour Point
46
23
Wax Modeling Concepts: Molecular Diffusion
Tbulk
Deposition
Twall
Twall < Tbulk and Twall < WAT
Dissolved wax is transported to the wall
and precipitates instantaneously
47
Deposition
48
24
Wax Modeling Concepts: Shear Deposition
49
Melting/Dissolution T
WDT
WAT
• If the Wax Dissolution Temperature
Pour Point
(WDT) is reached, the wax layer will
melt
• The melting rate is governed by
diffusion
• Can be limited by a user-given
kinetic parameter
50
25
Wax Deposition
Wax deposition rate for two crudes
140
D if f u s io n d e p o s it io n ra t e
120 KNK crude
( k g /m 3 E - 0 6 )
80
60
40
20
0
0 10 20 30 40 50
Temperature (C)
• Two crudes differences below 30 Deg. C
Viscosity
Temperature and shear rate
450
400 New tonian
350 30 s-1
Visco sity (CP)
300
250 100 s-1
200 300 s-1
150 500 s-1
100
50
0
0 20 40 60 80 100
Temperature (C)
26
Wax Deposition
• Inhibitors
– Slow down deposition by up to 5x
– Wax deposits are softer
– Pigging frequency is reduced
– Does not prevent deposition
53
Model Experience
• Tends to under predict wax deposition for single-
1. RRR (Rygg, phase systems
Rydahl & • Good agreement for multiphase oil systems
Rønningsen) • Somewhat conservative for gas condensate
systems
54
27
Well Flow Dynamics
Production / Injection Wells
• Flow Stability / Minimum stable flow rates / Slug Mitigation
– Tubing sizing / Production Optimisation
• Water accumulation / Corrosion rates
– Material selection / contact time of water as film
• Flow assurance, Wax , Hydrates
– Location of SCSSV
– MeOH / Glycol requirements
• Is GL injection stable?
– Annular heading? GLV choke size? Compressors?
– Density wave instability? Reservoir? Horizontal well?
– Multi GLV design? Single injection point?
– Deepwater GL injection line? Liquid injection?
– Riser instability? Riser GL?
– Optimum amount of gas lift gas?
– GL Well unloading? Compositional tracking?
• Which are the flow conditions at ESP depth?
• Slug Size (free gas volume)? Slug frequency?
• Optimum setting depth to improve pump Ev?
• Optimum pump size to improve pump Ev?
• Liner-Tubing-Casing optimum ID?
• Optimum ESP design & operating conditions? 56
28
Well Dynamic Simulation
Artificial Lift Design and Production Optimisation Solutions
• Is GL injection stable?
– Annular heading? GLV choke size? Compressors?
– Density wave instability? Reservoir? Horizontal well?
– Multi GLV design? Single injection point?
– Deepwater GL injection line? Liquid injection?
– Riser instability? Riser GL?
– Optimum amount of gas lift gas?
– GL Well unloading? Compositional tracking?
• Which are the flow conditions at ESP depth?
• Slug Size (free gas volume)? Slug frequency?
• Optimum setting depth to improve pump Ev?
• Optimum pump size to improve pump Ev?
• Liner-Tubing-Casing optimum ID?
• Optimum ESP design & operating conditions? 57
58
29
Gas Lift Well Stability
59
60
30
Interaction Between Downhole & Surface Orifice
61
0,80
0,75
0,70 ٛ Increasing reservoir pressure and gas
0,65
0,60
injection rate increases stability.
0,55
0,50
ٛ Increasing well depth, tubing diameter,
0,45 PI and system pressure decreases
0,40
0,35
stability
0,30 ٛ Instability occurs only when
0,25
0,20 PR − Psep
0,15 <1
0,10 ρ l gL
0,05
0,00
30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280 290 300 310
PR-Psep (bar)
SPE 84917
Two-phase vertical flow under gravity domination often is unstable, particularly in gas lift wells.
62
31
Typical Gas Lift Well Optimisation
Modelling concerns:
c) Heat Transfer
d) Non-constant Composition in
Tubing above Injection Point
Casing
e) Injecting dry gas?
63
• Is GL injection stable?
– Annular heading? GLV choke size? Compressors?
– Density wave instability? Reservoir? Horizontal well?
– Multi GLV design? Single injection point?
– Deepwater GL injection line? Liquid injection?
– Riser instability? Riser GL?
– Optimum amount of gas lift gas?
– GL Well unloading? Compositional tracking?
• Which are the flow conditions at ESP depth?
• Slug Size (free gas volume)? Slug frequency?
• Optimum setting depth to improve pump Ev?
• Optimum pump size to improve pump Ev?
• Liner-Tubing-Casing optimum ID?
• Optimum ESP design & operating conditions? 64
32
ESP – Design/Performance Optimisation
Pump Efficiency – Gas Locking
Using Olga
• Using Olga to model transient flow.
• Can also see slugs on amp charts.
• Have determined “classes” of problems:
o Slug dominated production
o Free gas dominated production.
Operational Strategies
• Using VSD’s to slow pump when gas slug hits
• VSD varies frequency to stabilize amp chart
• Using combinations of gas separators/gas handlers, which they refer to as “fluid conditioners.”
• Concerning well testing, need to measure the actual production rates vs. time, not just the
average daily production rate, to evaluate the slugging problems.
• Able to influence the drilling process to produce well bores that are more “friendly” for
production.
Reduced downtime by about 50%.
Able to control the VSD based on the amps being drawn by the motor.
Using an “auto orienting” pump intake, to place the intake on one side of the wellbore.
33
Well Flow Dynamics
Production / Injection Wells
• Flow Stability / Minimum stable flow rates / Slug Mitigation
– Tubing sizing / Production Optimisation
• Water accumulation / Corrosion rates
– Material selection / contact time of water as film
• Flow assurance, Wax , Hydrates
– Location of SCSSV
– MeOH / Glycol requirements
70
38
40 36
10
0
0
34
Well Flow Dynamics
Production / Injection Wells
• Flow Stability / Minimum stable flow rates / Slug Mitigation
– Tubing sizing / Production Optimisation
• Water accumulation / Corrosion rates
– Material selection / contact time of water as film
• Flow assurance, Wax , Hydrates
– Location of SCSSV
– MeOH / Glycol requirements
35
Well Dynamic Simulation
Horizontal Wells / Smart Wells Solutions
0.9
0.8
0.5
•
0.4
0.2
0.1
36
Well Dynamic Simulation
WAG Injection Solutions
37
Well Dynamic Simulation
Well Clean-up and Kick-off Solutions
5 MMscfd – Purple
10 MMscfd – Light Blue
15 MMscfd – Green
20 MMscfd – Black
25 MMscfd – Blue
5 MMscfd – Red
Dowhnole shut-in at 60 min
Surface shut-in at 65 minutes
38
Well Dynamic Simulation
Well Clean-up and Kick-off Solutions
5 MMscfd – Purple
10 MMscfd – Light Blue
15 MMscfd – Green
20 MMscfd – Black
25 MMscfd – Blue
5 MMscfd – Red
Dowhnole shut-in at 60 min
Surface shut-in at 65 minutes
7" x Tree ⊕⊕
W/Hd &
SCSSV
7
CRA 30"
O
20"
9-5/8
CRA
7-5/8
CRA
13-3.8"
9-5/8" CRA
7" CRA
78
39
Well Dynamic Simulation
Well Testing Solutions
80
40
Well Dynamic Simulation
Well Testing Solutions
81
When using gas lift (300 scf/min of N2 during the first 4 hours), it took approximately 4.3 hours
to clean-up the well.
82
41
WHP Comparison
3500
1
0.9
0.7
Pressure (psia)
2000 0.5
•
Pressure (Measured)
0.3
1000
500
0.1
0 0
0.0 100.0 200.0 300.0 400.0 500.0 600.0 700.0 800.0
The interaction between the nearwellbore reservoir and the well can
play a dominant role in the description of the dynamic behaviour of
the complete system
OLGA Multiphase Flow Transient OLGA Multiphase Flow Transient OLGA Multiphase Flow Transient
Simulator Simulator Simulator
84
42
The Virtual Gauge - Well Testing
OLGA Quasi-dynamic Reservoir Input
85
Pwh Ps Ts
Twg
Surface Choke
• Improved DST measurements Test Manifold
Closed
–
Pg Tg Bottomhole P-T measurements
Gas expansion effect
– After-flow effect
Gauge P correction
ΔP=ρgh
in-situ density Flowing bottomhole P-T
Formation Top
Pwf Pr Tr
Twf Datum
Formation Bottom
Rathole
86
43
The Virtual Gauge - Well Testing
DST Simulation – Build-up Test (Surface Shut-in)
Closed Wellhead
• Dowhhole Shut-in
Rising Bubbles Expanding
Phase Segregation
Qgf
87
Pwh Ps Ts
Twg
Surface Choke
Test Manifold
Diesel Closed
Completion Mud
Packer
Tester Valve open Downhole Choke effects
Opening - Closing Test Valve
Gauge
Pg Tg Bottomhole P-T measurements
Gauge P correction
ΔP=ρgh
in-situ density Flowing bottomhole P-T
Formation Top
Pwf Pr Tr
Twf
Datum
Formation Bottom
Rathole
88
44
The Virtual Gauge - Well Testing
DST Simulation - Well Test Data
Wellhead
Gauge
Qg (t) Pg (t)
nflow Point
Qfbh(t) Pfbh(t)
89
45
The Virtual Gauge - Well Testing
DST Simulation - Well Test Data
6000
5950
Pressure (psia)
5900
5850
5800
800 1300 1800 2300 2800 3300 3800 4300
Time (minutes)
46
Well Dynamic Simulation
Well Control and Workovers Solutions
94
47
What We Need to Know About
Dynamic Simulation
ID=8-in,
Gas Lift Depth=120 m
ID=8-in, Length=4.6 km
Annulus
ID=0.2159 m
PC
Depth=2840 m
Gas
Outlet
PCV
ID=2 m, Length=6 m
W1 W2 W3 W4 NLL=0.842 m LC
HHLL=1.687
Tubing ID=0.1143 m LLLL=0.315
Production
Separator Liquid
Riser Outlet
• Quasi-dynamic Reservoir: LCV
96
48
Integrated Modelling Application
Gas Lift Example
Reservoir-Well-GL-Flowline-Riser-Separator-Facilities
Gas rate
97
49
e-DPM
Real-Time Solutions (using APIS)
1) Minimise downtime
4) Enhance Operability
5) Inhibitor Management
99
DELIVERABLES / VALUE
Feasibility
Design / Engineering
Max Operational window
Operating Phil.
Procedures
Operator Training
Production Optimization
Monitoring
Training
Modifications
MODELING TOOLS
OLGA Engineering tool
OLGA models
I
50
Economic Drivers
Huge savings and risk minimisation
Economic Drivers
Huge savings and risk minimisation
102
51
Conclusions
• Technical, operational and HSE
integrity during design &
operation of production system
• Overall CAPEX and OPEX
effectiveness + Production and
Profit Optimisation
be dynamic
52