AQUNUM
Recognised by the OPM Flow parserGRID
The AQUNUM keyword defines the properties of numerical aquifers, including which grid blocks in the model should be utilized as part of the numerical aquifer. Each row entry in the AQUNUM keyword defines one numerical aquifer. Note that a numerical aquifer may consists of more than one grid cell, in order to better describe the water influx from the aquifer to the grid.
Notes:
Where NX, NY and NZ are the dimensions of the model as defined on the DIMENS keyword in the RUNSPEC section.
The keyword is followed by up to MXNAQN records as defined on the AQUDIMS keyword in the RUNSPEC section
Each record is terminated by a “/” and the keyword should be terminated by a “/”.
Numerical aquifers are modeled as one-dimensional, with aquifer flow assumed to be in the direction defined by LENGTH, and flux out of the aquifer is through the cross sectional area defined by AREA. Note that, all the aquifer cells must be isolated from the reservoir cells, with only the AQUCON connections connecting to the actual reservoir cells.
The values entered on the AQUNUM keyword are used to calculate the aquifer's pore volume and the transmissibility between the aquifer and the connected cell faces defined on the AQUCON keyword. Thus:
The aquifer’s pore volume is always calculated from the data entered on the AQUNUM keyword using[Pore~ Volume = AREA times LENGTH times PORO]and any modifications to the host cell values performed in either the and any modifications to the host cell values performed in either the GRID or EDIT sections are always ignored for cells declared as numerical aquifers cells.
For the transmissibility calculation either the cross-sectional area (AREA) defined on the AQUNUM keyword may be used or the connecting cell cross-sectional area by setting the AQUOPT1 variable on the AQUCON keyword.
In order to fully define a numerical aquifer one has to define the aquifer properties via the AQUNUM keyword, and how the aquifer is connected to the reservoir using the AQUCON keyword in the GRID or SOLUTION sections.
Note
If the AQUCON keyword has been utilized in the run deck then OPM Flow will write the AQUIFERN array to the *.INIT file in order to visualize the aquifer connections in OPM ResInsight.
This is accomplished by setting the AQUIFERN value to 2^(AQUID-1) for cells connected to aquifer AQUID. If a cell is connected to multiple numerical aquifers then AQUIFERN is summed for all aquifers connected to a cell. Note that connecting cells to multiple aquifers is best avoided.
Finally for cells representing the numerical aquifers themselves, AQUIFERN is set to minus AQUID.
Using one aquifer cell should generally be sufficient, provided the aquifer properties are constant. However, employing multiple cells may be appropriate if the aquifer properties vary with distance or depth, which is not uncommon. Secondly, having multiple cells may help to minimize throughput-related convergence problems by increasing the pore volumes exponentially away from the reservoir cells. Typically, three to five aquifer cells are employed under theses circumstances.
The aquifer pore volume can be used to define the amount of pressure support; whereas the aquifer transmissibility will influence the responsiveness of the aquifer.
See also the AQUDIMS keyword in the RUNSPEC section that defines the numerical aquifer dimensions.
- Records
- One or more records, each terminated by
/ - Items per record
- 12
Parameters
| No. | Name | Description | Type | Dimension | Field | Metric | Lab | Default |
|---|---|---|---|---|---|---|---|---|
| 1 | AQUIFER_manual: AQUID | AQUID is a positive integer greater than or equal to one and less than or equal to the maximum number of numerical aquifers as defined by the MXNAQN variable on the AQUDIMS keyword in the RUNSPEC section, that defines the aquifer to be connected to the grid. | INT | None | ||||
| 2 | I | A positive integer that defines the cell in the I-direction that represents the AQUID aquifer, and which must be greater than or equal to one and less than or equal to NX. | INT | None | ||||
| 3 | J | A positive integer that defines the cell in the J-direction that represents the AQUID aquifer, and which must be greater than or equal to one and less than or equal to NY. | INT | None | ||||
| 4 | K | A positive integer that defines the cell in the K-direction that represents the AQUID aquifer, and which must be greater than or equal to one and less than or equal to NZ. | INT | None | ||||
| 5 | CROSS_manual: AREA | AREA is a real positive value that defines the cross-sectional area of the aquifer used in calculating the aquifer connection transmissibility. The actual transmissibility between the numerical aquifer and the connected grid cell is the harmonic average of the aquifer connection transmissibility and the calculated connected cell transmissibility. The value entered for AREA does not effect the visualization of the cell in OPM ResInsight, as it is only used in calculating the transmissibility. Note that the AQUOPT1 variable on the AQUCON keyword in the GRID section allows one to use the value entered for AREA or to use the grid cell cross-sectional area instead for the transmissibility calculation. | DOUBLE | Length* | ft2 | m2 | cm2 | None |
| 6 | LENGTH | LENGTH is a real positive value that defines the length of the numerical aquifer. Similar to the AREA variable, LENGTH is not constrained by the original host cell size and the value entered does not effect the visualization of the cell in OPM ResInsight. | DOUBLE | Length | feet | m | cm | None |
| 7 | PORO | PORO is a real positive number greater than zero and less than or equal to one that assigns the porosity to the numerical aquifer. | DOUBLE | 1 | dimensionless | dimensionless | dimensionless | None |
| 8 | PERM | PERM is a real positive number that assigns the permeability to the numerical aquifer. | DOUBLE | Permeability | mD | mD | mD | None |
| 9 | DEPTHmanual: DATUM | DATUM is a real positive number that sets the reference datum depth of the numerical aquifer. Similar to the AREA variable, DATUM is not constrained by the original host cell depth and the value entered does not effect the visualization of the cell in OPM ResInsight. If defaulted then the depth of cell (I, J, K) will be used | DOUBLE | Length | feet | m | cm | Cell Depth |
| 10 | INITIAL_manual: PRESS | PRESS is a single positive value that defines the numerical aquifer pressure at DATUM. If PRESS is defaulted then the simulator will set the aquifer’s initial reservoir pressure to be in equilibrium with the cells the aquifer is contacted to. This is the preferred manner to initialize the numerical aquifer. | DOUBLE | Pressure | psia | barsa | atma | 1* |
| 11 | PVT_manual: PVTNUM | PVTNUM is positive integer greater than zero and less than the NTPVT variable on the TABDIMS keyword in the RUNSPEC section, that defines the PVTW table allocated to the numerical aquifer. If defaulted then the PVT tables allocated to cell (I, J, K) will be used. | INT | Cell PVTNUM | ||||
| 12 | SAT_manual: SATNUM | SATNUM is positive integer greater than zero and less than the NTSFUN variable on the TABDIMS keyword in the RUNSPEC section, that defines the saturation tables allocated to the numerical aquifer. If defaulted then the saturation tables allocated to cell (I, J, K) will be used. | INT | Cell SATNUM |
Example
Given the following grid and aquifer dimensions in the RUNSPEC section:
--
-- MAX MAX MAX
-- NDIVIX NDIVIY NDIVIZ
DIMENS
96 58 28 /
-- AQF AQF AQF AQF AQF AQF AQF AQF
-- MXAQN MXNAQC NIFTBL NRIFTB NANAQU NCAMAX MXNALI MXAAQL
AQUDIMS
3 3 1* 1* 1* 1* 1* 1* /
The following numerical aquifer definition:
--
-- NUMERICAL AQUIFER DESCRIPTION
--
-- AQUID - LOCATION - AQF AQF AQF AQF AQF AQF PVT SATNUM
-- NUMBER I J K AREA LENGTH PORO PERM DATUM PRES TAB TAB
--
AQUNUM
1 1 3 28 5.20E03 2.0E3 0.05 0.3 1* 1* 1* 1* /
1 1 2 28 5.20E06 2.0E6 0.05 0.3 1* 1* 1* 1* /
1 1 1 28 5.20E09 2.0E9 0.05 0.3 1* 1* 1* 1* /
/
defines one numerical aquifer consisting of three cells. The connection to the grid cells would take the form of:
--
-- NUMERICAL AQUIFER CONNECTIONS
--
-- ID ---------- BOX --------- CONNECT TRANS TRANS ADJOIN
-- NUMBER I1 I2 J1 J2 K1 K2 FACE MULT OPTN CELLS
AQUCON
1 1 1 4 58 28 28 'K+' 1.3 1* 1* /
1 2 96 1 58 28 28 'K+' 1.3 1* 1* /
/
that creates a basal aquifer115
Basal Aquifer: An aquifer located at the bottom of a geologic unit..
Basal Aquifer: An aquifer located at the bottom of a geologic unit.
| Note If the AQUCON keyword has been utilized in the run deck then OPM Flow will write the AQUIFERN array to the *.INIT file in order to visualize the aquifer connections in OPM ResInsight. This is accomplished by setting the AQUIFERN value to 2^(AQUID-1) for cells connected to aquifer AQUID. If a cell is connected to multiple numerical aquifers then AQUIFERN is summed for all aquifers connected to a cell. Note that connecting cells to multiple aquifers is best avoided. Finally for cells representing the numerical aquifers themselves, AQUIFERN is set to minus AQUID. |
|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|Manual source: parts/chapters/subsections/6.3/AQUNUM.fodt