
- •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.

An alternative is therefore to model the aquifer by some artificial means. There are several ways to do this, but we will concentrate on the so-called numerical aquifers, which also can be modelled in Eclipse. The reasons for using a special model for the aquifer, are
•Practical: it saves memory space and computer time
•No loss of accuracy. Since aquifers are large one-phase volumes, the computations are simple, and as accurate on large grid blocks as on small.
•Detail water flow pattern is not needed in the water zone, only the energy (pressure) support from aquifer to reservoir, and water flux between the two are needed.
Actually, the flux and energy support can be computed with sufficient accuracy in a one-dimensional abstraction of the aquifer, which is how it is implemented. With this knowledge we see that a sufficient description of an aquifer doesn’t require shape, only volumes. However, the interface between the aquifer and the reservoir requires a shape definition, which is only logical.
Hence there are two aspects of setting up an aquifer, the aquifer definition (volume, properties), and how it interfaces to the reservoir.
Aquifer definition
The way Eclipse does it, ordinary grid cells in the reservoir are used to define aquifers. The cells will have all their properties redefined, but apart from that they are still part of the grid, and cannot be regarded as having been removed from the grid. E.g. connectivity between grid cells is not changed just because a cell is defined as an aquifer cell. This means that aquifer cells should be completely isolated from the rest of the reservoir, i.e. completely surrounded by inactive cells.
Most aquifers can be satisfactory defined by one grid cell only. However, if the aquifer properties and/or depth have large variation, it can be a good idea to try to capture this variation by using more cells to define the aquifer. More than three or four cells should seldom or never be necessary.
To verify this it is a nice exercise to build a simple grid including a large water zone which is explicitly modelled. Then model the water zone as a numerical aquifer using one, two, four and ten cells. Most probably, the results from these cases will be as good as equal.
Figure 20 shows some examples of ways to use grid cells to define aquifers, and the connectivity to expect.
A
B1 B2 B3
C1C2
Active cells
Inactive cells
D
Figure 20. Examples of use of inactive cells to define aquifers (grid viewed from above)
Cell A is a typical example of using an isolated inactive cell to define an aquifer. The cell has no connection to the surrounding world at all, so the necessary flow connections will be defined by the user.
75
Cells B1, B2, B3 are connected to each other, but isolated from the rest of the grid. The three cells can either be used to define a single aquifer, which obviously will have internal communication, or the cells could define two or three different aquifers, such that the aquifers communicate with each other. Cells C1 and C2 are isolated from each other, and from the rest of the world. They can be used to define two different aquifers. Using such a configuration for a single aquifer would be meaningless, as there would be no connectivity between the two different parts of the aquifer.
Cell D is not suited for use as an aquifer cell, since it is in direct contact with an active cell in the grid. In addition to the connectivity defined by the user, fluid in cell D will flow directly into the reservoir, which is probably not intended.
Having decided which (inactive) cells from the grid to use as aquifer cells, the keyword to define the aquifer is AQUNUM (numerical aquifer). The keyword data is comprised of arbitrary many lines, each line containing definition for exactly one aquifer cell on the original grid. The syntax for each line is,
AQ-id AI AJ AK X-sect_area Length Poro Perm Depth P_init PVT-table SatNum
AQ-id
an integer identifying the aquifer. If several lines are used to define the same aquifer, only the cell defined in the first line of definition will be in contact with the reservoir, as defined by the AQUCON keyword below.
AI AJ AK
The (I, J, K)-index for the grid cell to use as aquifer cell
X-sect_area Length
The gross volume of the aquifer cell is the product of these two numbers. Recall that the aquifer is modelled as one-dimensional – flow is assumed to be in the direction defined by Length, and flux out of the aquifer is through the cross sectional area defined by X-sect_area.
Poro Perm Depth
Normally we will define porosity, permeability and a typical depth for the aquifer, taken from known properties of the water zone. If these items are defaulted, the values will be taken from the grid cell as it is defined in the grid.
P_init
Initial pressure in the aquifer cell. It is recommended to default this value, whereby it will be computed from the fluid gradients and the depth defined above.
PVT-table SatNum
These items are best defaulted. (Ref Eclipse documentation if alternative is needed.)
Example.
We will define four different aquifers, corresponding to examples A, B, and C in Figure 20. Cell A has aquifer-id 1, aquifer-id 2 is used for cells B, and aquifer-id 3 and 4 for cells C. The view in Figure 20 is assumed to be the top layer (K=1), at the edge of the reservoir (I=1 at left end, J=1 at top).
AQUNUM |
AI AJ AK |
X-sect |
Length |
Poro |
Perm |
Depth P_init |
|
|||
-- AQ-ID |
/ |
|||||||||
1 |
3 |
1 |
1 |
2E7 |
3000 |
0.3 |
800 |
3000 |
1* |
|
2 |
3 |
3 |
1 |
6E6 |
2500 |
0.28 |
1000 |
2500 |
1* |
/ |
2 |
2 |
3 |
1 |
6E6 |
2500 |
0.28 |
800 |
2800 |
1* |
/ |
2 |
1 |
3 |
1 |
6E6 |
2500 |
0.28 |
500 |
3100 |
1* |
/ |
3 |
1 |
5 |
1 |
5E7 |
6000 |
0.25 |
300 |
3500 |
1* |
/ |
4 |
3 |
5 |
1 |
4E6 |
4000 |
0.24 |
100 |
2700 |
1* |
/ |
/ |
|
|
|
|
|
|
|
|
|
|
Note that the dimensions, petrophysics, and depth of the aquifer cells bear no relation to what was defined in the grid originally at all.
Flow between the reservoir and aquifer 2 will only be to/from cell (3, 3, 1), since this cell is the first to be defined in the keyword. Also note the reason for using three cells for aquifer 2, to capture the depth increase (higher fluid pressure, larger compressive energy) and decrease in permeability with depth.
76