
- •Basics of Reservoir Simulation
- •with the
- •Eclipse Reservoir Simulator
- •Lecture Notes
- •Øystein Pettersen
- •Introduction
- •Prerequisites
- •1. Overview – minimum required data input
- •1.1 The grid
- •1.2 Petrophysics
- •1.3 Fluid properties
- •1.4 Soil properties
- •1.5 Equilibration
- •1.6 Well specification
- •1.7 Dynamics
- •1.8 Output
- •1.9 Simple Eclipse data file contents
- •A. Syntax
- •B. Data file (“BASIC data input example”)
- •A note on units
- •2. The RUNSPEC section
- •Grid dimension – keyword DIMENS
- •Phases
- •Unit system
- •Start date
- •Unified / Non-unified files (chapter 11)
- •Data checking only
- •Table dimensions
- •EQLDIMS
- •NSTACK (chapters 17-18)
- •Aquifer specifications – AQUDIMS (chapter 14)
- •Grid options (chapter 13)
- •Rock compressibility options (chapter 6)
- •Local Grid Refinement (chapter 15)
- •3. Structured grids (Corner point grids) (GRID section)
- •The Corner Point Grid
- •Defining a corner point grid in Eclipse
- •Moderately complex grids – FILL
- •4. Petrophysics (GRID section)
- •Average permeability
- •Transmissibility
- •Inactive cells
- •5. Fluid properties (PROPS section)
- •Tables in Eclipse
- •Relative permeability and Capillary Pressure
- •Two-phase curves (water – oil)
- •Three-phase relative permeabilities
- •PVT data
- •Water
- •Dead Oil
- •Live Oil
- •6. Soil compressibility (PROPS section)
- •7. Initialisation (SOLUTION section)
- •Datum depth
- •Contacts
- •Equilibrium – discussion – advanced issues
- •8. Time dependent input data (SCHEDULE section)
- •8.1 Well definitions and control
- •Well Specification (WELSPECS keyword)
- •Well Completions (COMPDAT keyword)
- •Production / Injection data (Keywords WCONPROD / WCONINJE)
- •Economic well constraints (keywords WECON, WECONINJ)
- •Other often used Well control keywords
- •8.2 Time stepping
- •Order of actions
- •8.3 Convergence Control I (keyword TUNING)
- •9. Regions
- •10. Simplified input and modification of Eclipse arrays
- •EQUALS
- •ADD, MULTIPLY
- •COPY
- •COPYBOX
- •11. Eclipse output, formats and files
- •File names
- •Textual output
- •The RPTXXX keywords
- •Time dependent vectors – SUMMARY data
- •Restart data and restart files
- •12. Restarting a simulation
- •The SKIPREST keyword
- •13. Fault modelling – Non-neighbour connections
- •The 7-point stencil
- •The fault layout – non-neighbour connections
- •Fault transmissibility multipliers
- •Defining a fault manually – the ADDZCORN keyword
- •14. Aquifer Modelling (GRID section)
- •Aquifer definition
- •Aquifer connection to reservoir
- •15. Local Grid Refinement
- •15.2 LGR on an irregular volume – Amalgamation
- •15.3 Wells on local grids – Horizontal wells
- •15.4 Horizontal wells and friction
- •16. Numerical Solution of the Flow Equations
- •The IMPES method
- •Solution of Non-linear Equations – the Newton-Raphson method
- •17. Iteration methods for linear systems
- •Direct, simple approach
- •The Gauss-Seidel method
- •Accelerators – the point SOR method
- •Conjugate Gradients – ORTHOMIN
- •Preconditioning
- •Preconditioning and Orthomin
- •Determining a preconditioner – Nested Factorisation
- •18. Convergence Control II – TUNING parameters
- •TUNING keyword summarized
- •19. Non-neighbour Connections and System Structure
- •A. GRF files in GRAF
- •A simple straightforward GRF file
- •Advanced GRF file
- •B. Some Considerations Regarding Grid Consistency
- •Grids planned for use in rock mechanics simulations
- •Embedding
- •Non-vertical coordinate lines
- •Honouring material properties of non-reservoir rock.
B. Data file (“BASIC data input example”)
Following is an example of minimum required data input to Eclipse. The keywords are commented using lines starting with two dashes. Full understanding of the file is not to be expected at this stage. Note: Section headers are in bold for clarity only. An eclipse data file must be in text-only format, not. e.g. saved as a Word-document.
RUNSPEC ==============================================================
--Title is used as header on output, to identify run TITLE
Example simple Eclipse file, lecture notes
--Specify dimension of model, NX, NY, NZ
DIMENS |
NY |
NZ |
|
-- NX |
/ |
||
10 |
3 |
3 |
--Phases included (oil and water, i.e. 2-phase run)
OIL WATER
--Units to use, alternatives are METRIC, FIELD or LAB FIELD
--Specify maximum values for well data (# means “max number of”) WELLDIMS
-- |
#wells #cell-connections |
Next three values are defaulted |
|
|
4 |
5 |
3* / |
-- |
Start date, |
i.e. “Day 0” in simulation |
|
START |
/ |
|
|
|
1 'JAN' 2004 |
|
GRID ===========================================================
--Turn off writing of data to report file NOECHO
--Definition of grid cells. From DIMENS keyword above, we have:
--N = NX*NY*NZ = 10*3*3 = 90.
--In each cell DX = DY = 1000 ft and DZ = 50 ft
--(Dimension measured in feet since units were defined as FIELD above)
--The grid cells are then defined by:
DX
90*1000 /
DY
90*1000 /
DZ
90*50 /
--So far the grid shape has been defined, but we also need to define at
--which depth (measured from sea level) the cell tops are. We assume a
--horizontal reservoir, and input depths which are consistent with the
--DZ-values we defined above.
TOPS
30*5000 30*5050 30*5100 /
--Permeabilities are constant, so define 90 equal values. PERMX
90*200 / PERMY
90*200 / PERMZ
90*200 /
--Porosity is constant in each layer, so define 30 values in layer 1,
--30 in layer 2, and 30 in layer 3. Each layer contains 30 cells. PORO
30*0.3
30*0.23
30*0.18 /
11
--Request to write an INIT (initial) file. Contains all data used to
--initialize the model
INIT
-- Turn report writing back on ECHO
PROPS ===============================================================
--Relative permeability for water and oil,
--and water-oil capillary pressure, as a function of water saturation SWOF
-- Sw |
Krw |
Krow |
Pcow |
|
0.220 |
0.0000 |
1.0000 |
0 |
|
0.300 |
0.0700 |
0.4000 |
0 |
|
0.400 |
0.1500 |
0.1250 |
0 |
|
0.500 |
0.2400 |
0.0649 |
0 |
|
0.600 |
0.3300 |
0.0048 |
0 |
|
0.800 |
0.6500 |
0.0 |
0 |
|
0.900 |
0.8300 |
0.0 |
0 |
/ |
1.000 |
1.0000 |
0.0 |
0 |
--PVT properties for water.
--(Pref: Reference pressure for rest of data (psi)
--Bw: Volume formation factor for water
--Cw: Water compressibiblity
--ViscW: Water viscosity )
PVTW |
Bw |
Cw |
|
ViscW |
|
-- Pref |
|
0 / |
|||
4014.7 |
1.029 |
3.13D-6 |
0.31 |
||
-- PVT properties for oil |
|
|
|||
PVDO |
Bo |
viscO |
|
|
|
-- P |
|
|
|
||
3337 |
1.2600 |
1.042 |
|
|
|
3725 |
1.2555 |
1.072 |
|
|
|
4139.5 |
1.2507 |
1.096 |
|
|
|
4573.2 |
1.2463 |
1.118 |
|
|
|
5053.9 |
1.24173 |
1.151 |
|
|
|
5487.5 |
1.2377 |
1.174 |
/ |
|
|
5813.9 |
1.2356 |
1.2 |
|
|
--Dead oil: Rs (Gas resolution factor) is constant RSCONST
--Rs Bubble-point-pressure
0.43337.0 /
--Specify constant rock compressibility.
ROCK |
Cr |
|
-- Pref |
/ |
|
14.7 |
3.0D-6 |
-- |
Fluid densities at |
surface conditions |
|
DENSITY |
Water Gas |
|
|
-- |
Oil |
|
49.164.79 0.06054 /
SOLUTION =============================================================
EQUIL |
|
|
|
|
|
|
-- |
DD = Datum depth, the depth to which all reports will be referenced. |
|||||
-- |
DD |
Pressure@DD |
OWC |
Pcow(OWC) |
Default rest of data items |
|
|
5000 |
4800 |
6000 |
0.0 |
6* |
/ |
12
SUMMARY ==============================================================
--List data vectors which we want stored for graphics post-processing
--Field Oil Production Rate
FOPR
--Field Oil Production Total FOPT
--Field Water Cut
FWCT
--Field Pressure (averaged reservoir pressure)
FPR
--Field Oil In Place
FOIP
--Field Water Production Rate FWPR
--Field Water Injection Rate FWIR
--Well Water Cut for all wells WWCT
/
SCHEDULE ==============================================================
--Specify output of graphics result files for cell data, and times which
--to write these. (Details to come later)
RPTRST
BASIC=5 NORST=1 FREQ=6 /
-- Well specification: Give names, positions (i, j) and main phase of wells
WELSPECS |
Z(bhp) |
prefPhase |
--wname group i j |
||
'WP1' 'G' 10 2 |
1* |
'OIL' / |
/ |
|
|
--(Note two slashes, one terminates each well, one terminates the keyword)
--Completion data, the well is open to the reservoir in cells in layers
--from k_hi to k_lo.
COMPDAT |
ic |
jc |
k_hi k_lo open/shut 2*Don’tCare well_diam |
Default... |
||||
--wname |
||||||||
'WP1' |
10 |
2 |
1 |
2 |
'OPEN' |
2* |
0.5 4* |
/ |
/ |
|
|
|
|
|
|
|
|
--Production data, specify producing rate and constraints
--Well WP1 will produce with an oil rate of 2000 STB/day constrained by
--that well’s flowing bottomhole pressure is minimum 3350 psi.
--(Details to come later)
WCONPROD |
|
|
|
orat 4*Default bhpmin |
Rest default... |
||
--wname open/shut ctrlmode |
|||||||
'WP1' |
'OPEN' |
'ORAT' |
2000 |
4* |
3350 |
/ |
|
/ |
|
|
|
|
|
|
|
-- Milestone dates, e.g. to write results |
|
|
|||||
DATES |
|
/ |
|
|
|
|
|
1 'JUL' 2004 |
|
|
|
|
|
||
1 'JAN' 2005 |
/ |
|
|
|
|
|
|
/ |
|
|
|
|
|
|
|
-- Change production rate to 4000 STB/day |
|
|
|||||
WCONPROD |
|
|
|
orat 4*Default bhpmin |
Rest default... |
||
--wname open/shut ctrlmode |
|||||||
'WP1' |
'OPEN' |
'ORAT' |
4000 |
4* |
3350 |
/ |
|
/ |
|
|
|
|
|
|
|
DATES
1 'JUL' 2005 /
/
END
13
A note on units
Eclipse accepts data in metric units, as well as field or lab units. Specification of which unit convention to use is done in the RUNSPEC section. Since metric units are widely used by Norwegian oil companies and international companies based in Norway, unless otherwise specified all the examples in these notes will assume metric units are used.
14

Figure 1. Regular Cartesian Grid
Figure 2. Regular XY-Cartesian Grid
15

Figure 3. Irregular structured grid
Figure 4. Unstructured grid
16