add a simplified version of the SPE1 dataset

the only remarkable aspect of this deck is that the current version of
the CPR preconditioner in opm-autodiff breaks down:

```
and@heuristix:~/src/opm-autodiff|master > ./bin/sim_fibo_ad deck_filename=../opm-data/simplespe1/SIMPLESPE1.DATA

================    Test program for fully implicit three-phase black-oil flow     ===============

---------------    Reading parameters     ---------------
deck_filename found at /, value is ../opm-data/simplespe1/SIMPLESPE1.DATA
strict_parsing not found. Using default value 'true'.
deck_filename found at /, value is ../opm-data/simplespe1/SIMPLESPE1.DATA
output not found. Using default value 'true'.
output_interval not found. Using default value '1'.
output_dir not found. Using default value '.'.
pvt_tab_size not found. Using default value '-1'.
sat_tab_size not found. Using default value '-1'.
threephase_model not found. Using default value 'gwseg'.
use_cpr not found. Using default value 'true'.
output not found. Using default value 'true'.
output_dir not found. Using default value 'output'.
output not found. Using default value 'true'.
output_vtk not found. Using default value 'true'.
output_dir not found. Using default value 'output'.
output_interval not found. Using default value '1'.

================ Starting main simulation loop ===============

---------------    Simulation step number 0    ---------------
      Current time (days)     0
      Current stepsize (days) 1
      Total time (days)       1216

dp_max_rel not found. Using default value '1e+09'.
ds_max not found. Using default value '0.2'.
drs_max_rel not found. Using default value '1e+09'.
relax_max not found. Using default value '0.5'.
max_iter not found. Using default value '15'.
relax_type not found. Using default value 'dampen'.

Iteration         Residual
        0        1.29402173
terminate called after throwing an instance of 'Dune::ISTLError'
Aborted
```

on the other hand, everything works fine if the CPR preconditioner is
not used:

```
and@heuristix:~/src/opm-autodiff|master > ./bin/sim_fibo_ad deck_filename=../opm-data/simplespe1/SIMPLESPE1.DATA use_cpr=false

================    Test program for fully implicit three-phase black-oil flow     ===============

---------------    Reading parameters     ---------------
deck_filename found at /, value is ../opm-data/simplespe1/SIMPLESPE1.DATA
strict_parsing not found. Using default value 'true'.
deck_filename found at /, value is ../opm-data/simplespe1/SIMPLESPE1.DATA
output not found. Using default value 'true'.
output_interval not found. Using default value '1'.
output_dir not found. Using default value '.'.
pvt_tab_size not found. Using default value '-1'.
sat_tab_size not found. Using default value '-1'.
threephase_model not found. Using default value 'gwseg'.
use_cpr found at /, value is 'false'.
linsolver not found. Using default value 'umfpack'.
output not found. Using default value 'true'.
output_dir not found. Using default value 'output'.
output not found. Using default value 'true'.
output_vtk not found. Using default value 'true'.
output_dir not found. Using default value 'output'.
output_interval not found. Using default value '1'.

================ Starting main simulation loop ===============

---------------    Simulation step number 0    ---------------
      Current time (days)     0
      Current stepsize (days) 1
      Total time (days)       1216

dp_max_rel not found. Using default value '1e+09'.
ds_max not found. Using default value '0.2'.
drs_max_rel not found. Using default value '1e+09'.
relax_max not found. Using default value '0.5'.
max_iter not found. Using default value '15'.
relax_type not found. Using default value 'dampen'.

Iteration         Residual
        0        1.29402173
        1      0.0614937239
        2    1.45760381e-05
Fully implicit solver took: 0.097702 seconds.

---------------    Simulation step number 1    ---------------
      Current time (days)     1
      Current stepsize (days) 14
      Total time (days)       1216

...
```
This commit is contained in:
Andreas Lauser 2014-09-23 12:15:55 +02:00
parent 66bce3d00d
commit 0ccc2e2608
2 changed files with 499 additions and 0 deletions

20
simplespe1/README Normal file
View File

@ -0,0 +1,20 @@
This directory contains an simple input data file which is based on
the first comparison project of the society of petroleum engineers
[1]. The dataset consists of a single file called 'SIMPLESPE1.DATA'
and uses the ECL format which is used by many standard commercial
simulators for oil reservoirs.
The simplifications are that the grid consists only of a single layer
of cells, and that there is no injector well specified. (i.e., only a
producer is present.)
Just like for the SPE1 dataset, the ownership of the copyright of this
file is unclear: Some parts seem to be written by Schlumberger [2],
which in turn seem to be modified by Jon Kleppe of NTNU [3] and some
(mainly saturation function related) data seems to be calculated by
SINTEF using MRST. Finally some modifications where done by Ove
Savareid and Andreas Lauser in the context of the OPM project.
[1] http://www.spe.org/web/csp/datasets/set01.htm
[2] http://slb.com
[3] http://www.ipt.ntnu.no/~kleppe/TPG4160/ex3/ODEH.DATA

479
simplespe1/SIMPLESPE1.DATA Normal file
View File

@ -0,0 +1,479 @@
-- =====================================================================
-- THIS IS THE FIRST SPE COMPARISON PROBLEM,"COMPARISON OF SOLUTIONS TO A
-- THREE-DIMENSIONAL BLACK-OIL RESERVOIR SIMULATION PROBLEM", REPORTED
-- BY AZIS AND ODEH AT THE SPE SYMPOSIUM ON RESERVOIR SIMULATION ,
-- JANUARY 1981. IT IS A NON SWELLING AND SWELLING STUDY. A REGULAR
-- GRID WITH TWO WELLS (INJECTOR AND PRODUCER)AND A IMPES SOLUTION METHOD
-- IS USED FOR THIS SIMULATION.THE PRODUCTION IS CONTROLLED BY FLOW RATE
-- AND MIN. BHP. OIL RATE, GOR, PRESSURE AND GAS SATURATION ARE TO BE REPORTED.
-- =====================================================================
RUNSPEC
-- TITLE
-- ODEH PROBLEM - IMPLICIT OPTION - 1200 DAYS /
DIMENS
10 10 1 /
NONNC
OIL
WATER
GAS
DISGAS
FIELD
EQLDIMS
1 100 10 1 1 /
TABDIMS
1 1 17 12 1 12 /
WELLDIMS
2 1 1 2 /
NUPCOL
4 /
START
19 'OCT' 1982 /
NSTACK
24 /
--FMTOUT
--FMTIN
UNIFOUT
UNIFIN
--NOSIM
--IMPES
GRID ================================================================
-------- IN THIS SECTION , THE GEOMETRY OF THE SIMULATION GRID AND THE
-------- ROCK PERMEABILITIES AND POROSITIES ARE DEFINED.
------------------------------------------------------------------------
-- THE X AND Y DIRECTION CELL SIZES ( DX, DY ) AND THE POROSITIES ARE
-- CONSTANT THROUGHOUT THE GRID. THESE ARE SET IN THE FIRST 3 LINES
-- AFTER THE EQUALS KEYWORD. THE CELL THICKNESSES ( DZ ) AND
-- PERMEABILITES ARE THEN SET FOR EACH LAYER. THE CELL TOP DEPTHS
-- ( TOPS ) ARE NEEDED ONLY IN THE TOP LAYER ( THOUGH THEY COULD BE.
-- SET THROUGHOUT THE GRID ). THE SPECIFIED MULTZ VALUES ACT AS
-- MULTIPLIERS ON THE TRANSMISSIBILITIES BETWEEN THE CURRENT LAYER
-- AND THE LAYER BELOW.
DXV
10*1000.0
/
DYV
10*1000.0
/
DZV
100.0
/
PORO
100*0.3
/
PERMX
100*500.0
/
PERMY
100*500.0
/
-- Note: ignoring MULTZ!
-- layer 1-2 'MULTZ' 0.64 /
-- layer 2-3 'MULTZ' 0.265625 /
-- Reducing PERMZ a little instead
PERMZ
100*500.0
/
TOPS
100*8325.0
/
INIT
PROPS ===============================================================
-------- THE PROPS SECTION DEFINES THE REL. PERMEABILITIES, CAPILLARY
-------- PRESSURES, AND THE PVT PROPERTIES OF THE RESERVOIR FLUIDS
----------------------------------------------------------------------
-- WATER RELATIVE PERMEABILITY AND CAPILLARY PRESSURE ARE TABULATED AS
-- A FUNCTION OF WATER SATURATION.
--
-- Generated with MRST's family_1() function from the original deck.
-- SWAT KRW KRO PCOW
SWOF
0.120000000000000 0 1.000000000000000 0
0.121000000000000 0.000000011363636 1.000000000000000 0
0.140000000000000 0.000000227272727 0.997000000000000 0
0.170000000000000 0.000000568181818 0.980000000000000 0
0.240000000000000 0.000001363636364 0.700000000000000 0
0.319999999000000 0.000002272727261 0.350000004375000 0
0.370000000000000 0.000002840909091 0.200000000000000 0
0.420000000000000 0.000003409090909 0.090000000000000 0
0.520000000000000 0.000004545454545 0.021000000000000 0
0.570000000000000 0.000005113636364 0.010000000000000 0
0.620000000000000 0.000005681818182 0.001000000000000 0
0.720000000000000 0.000006818181818 0.000100000000000 0
0.820000000000000 0.000007954545455 0.000000000000000 0
1.000000000000000 0.000010000000000 0 0
/
SGOF
0 0 1.000000000000000 0
0.001000000000000 0 1.000000000000000 0
0.020000000000000 0 0.997000000000000 0
0.050000000000000 0.005000000000000 0.980000000000000 0
0.120000000000000 0.025000000000000 0.700000000000000 0
0.199999999000000 0.074999999375000 0.350000004375000 0
0.200001100000000 0.075000000000000 0.350000000000000 0
0.250000000000000 0.125000000000000 0.200000000000000 0
0.300000000000000 0.190000000000000 0.090000000000000 0
0.400000000000000 0.410000000000000 0.021000000000000 0
0.449999999000000 0.599999996200000 0.010000000220000 0
0.450001100000000 0.600000000000000 0.010000000000000 0
0.500000000000000 0.720000000000000 0.001000000000000 0
0.600000000000000 0.870000000000000 0.000100000000000 0
0.700000000000000 0.940000000000000 0 0
0.850000000000000 0.980000000000000 0 0
0.880000000000000 0.984000000000000 0 0
/
-- PVT PROPERTIES OF WATER
--
-- REF. PRES. REF. FVF COMPRESSIBILITY REF VISCOSITY VISCOSIBILITY
PVTW
4014.7 1.029 3.13E-6 0.31 0 /
-- ROCK COMPRESSIBILITY
--
-- REF. PRES COMPRESSIBILITY
ROCK
14.7 3.0E-6 /
-- SURFACE DENSITIES OF RESERVOIR FLUIDS
--
-- OIL WATER GAS
DENSITY
49.1 64.79 0.06054 /
-- PVT PROPERTIES OF DRY GAS (NO VAPOURISED OIL)
-- WE WOULD USE PVTG TO SPECIFY THE PROPERTIES OF WET GAS
--
-- PGAS BGAS VISGAS
PVDG
14.7 166.666 0.008
264.7 12.093 0.0096
514.7 6.274 0.0112
1014.7 3.197 0.014
2014.7 1.614 0.0189
2514.7 1.294 0.0208
3014.7 1.080 0.0228
4014.7 0.811 0.0268
5014.7 0.649 0.0309
9014.7 0.386 0.047 /
-- PVT PROPERTIES OF LIVE OIL (WITH DISSOLVED GAS)
-- WE WOULD USE PVDO TO SPECIFY THE PROPERTIES OF DEAD OIL
--
-- FOR EACH VALUE OF RS THE SATURATION PRESSURE, FVF AND VISCOSITY
-- ARE SPECIFIED. FOR RS=1.27 AND 1.618, THE FVF AND VISCOSITY OF
-- UNDERSATURATED OIL ARE DEFINED AS A FUNCTION OF PRESSURE. DATA
-- FOR UNDERSATURATED OIL MAY BE SUPPLIED FOR ANY RS, BUT MUST BE
-- SUPPLIED FOR THE HIGHEST RS (1.618).
--
-- RS POIL FVFO VISO
-- PVTO
-- 0.001 14.7 1.062 1.04 /
-- 0.0905 264.7 1.15 0.975 /
-- 0.18 514.7 1.207 0.91 /
-- 0.371 1014.7 1.295 0.83 /
-- 0.636 2014.7 1.435 0.695 /
-- 0.775 2514.7 1.5 0.641 /
-- 0.93 3014.7 1.565 0.594 /
-- 1.270 4014.7 1.695 0.51
-- 5014.7 1.671 0.549
-- 9014.7 1.579 0.74 /
-- 1.618 5014.7 1.827 0.449
-- 9014.7 1.726 0.605 /
-- /
PVTO
1.0000000000e-03
1.4700000000e+01 1.0620000000e+00 1.0400000000e+00
1.0147000000e+03 1.0469628319e+00 1.1195294118e+00
5.0147000000e+03 9.8932035398e-01 1.5090196078e+00
/
9.0500000000e-02
2.6470000000e+02 1.1500000000e+00 9.7500000000e-01
1.2647000000e+03 1.1337168142e+00 1.0495588235e+00
5.2647000000e+03 1.0712979351e+00 1.4147058824e+00
/
1.8000000000e-01
5.1470000000e+02 1.2070000000e+00 9.1000000000e-01
1.5147000000e+03 1.1899097345e+00 9.7958823529e-01
5.5147000000e+03 1.1243970501e+00 1.3203921569e+00
/
3.7100000000e-01
1.0147000000e+03 1.2950000000e+00 8.3000000000e-01
2.0147000000e+03 1.2766637168e+00 8.9347058824e-01
6.0147000000e+03 1.2063746313e+00 1.2043137255e+00
/
6.3600000000e-01
2.0147000000e+03 1.4350000000e+00 6.9500000000e-01
3.0147000000e+03 1.4146814159e+00 7.4814705882e-01
7.0147000000e+03 1.3367935103e+00 1.0084313725e+00
/
7.7500000000e-01
2.5147000000e+03 1.5000000000e+00 6.4100000000e-01
3.5147000000e+03 1.4787610619e+00 6.9001764706e-01
7.5147000000e+03 1.3973451327e+00 9.3007843137e-01
/
9.3000000000e-01
3.0147000000e+03 1.5650000000e+00 5.9400000000e-01
4.0147000000e+03 1.5428407080e+00 6.3942352941e-01
8.0147000000e+03 1.4578967552e+00 8.6188235294e-01
/
1.2700000000e+00
4.0147000000e+03 1.6950000000e+00 5.1000000000e-01
5.0147000000e+03 1.6710000000e+00 5.4900000000e-01
9.0147000000e+03 1.5790000000e+00 7.4000000000e-01
/
1.6180000000e+00
5.0147000000e+03 1.8270000000e+00 4.4900000000e-01
9.0147000000e+03 1.7260000000e+00 6.0500000000e-01
/
/
-- OUTPUT CONTROLS FOR PROPS DATA
-- ACTIVATED FOR SOF3, SWFN, SGFN, PVTW, PVDG, DENSITY AND ROCK KEYWORDS
RPTPROPS
1 1 1 0 1 1 1 1 /
SOLUTION ===============================================================
-------- THE SOLUTION SECTION DEFINES THE INITIAL STATE OF THE SOLUTION
-------- VARIABLES (PHASE PRESSURES, SATURATIONS AND GAS-OIL RATIOS)
------------------------------------------------------------------------
-- DATA FOR INITIALISING FLUIDS TO POTENTIAL EQUILIBRIUM
--
-- DATUM DATUM OWC OWC GOC GOC RSVD RVVD SOLN
-- DEPTH PRESS DEPTH PCOW DEPTH PCOG TABLE TABLE METH
-- EQUIL
-- 8400 4800 8500 0 8200 0 1 0 0 /
-- VARIATION OF INITIAL RS WITH DEPTH
--
-- DEPTH RS
-- RSVD
-- 8200 1.270
-- 8500 1.270 /
-- PRESSURE
-- 100*3.297832774859256e7
-- 100*3.302313357125603e7
-- 100*3.309483500720813e7
-- /
PRESSURE
100*4800.00000000000
/
SWAT
100*0.12
/
SGAS
100*0.0
/
-- RS
-- 300*226.1966570852417
-- /
RS
100*1.27
/
-- OUTPUT CONTROLS (SWITCH ON OUTPUT OF INITIAL GRID BLOCK PRESSURES)
RPTSOL
RESTART=2 PRES DENW DENO DENG POIL PGAS PWAT VOIL VGAS VWAT KRO KRG KRW RSSAT /
-- RESTART=2 DENSO DENSG DENSW POIL PGAS PWAR VOIL VGAS VWAT KRO KRG KRW RSSAT /
SUMMARY ===============================================================
-------- THIS SECTION SPECIFIES DATA TO BE WRITTEN TO THE SUMMARY FILES
-------- AND WHICH MAY LATER BE USED WITH THE ECLIPSE GRAPHICS PACKAGE
------------------------------------------------------------------------
EXCEL
SEPARATE
--REQUEST PRINTED OUTPUT OF SUMMARY FILE DATA
RUNSUM
-- FIELD OIL PRODUCTION
FOPR
-- WELL GAS-OIL RATIO FOR PRODUCER
WGOR
-- 'PRODUCER'
/
-- WELL BOTTOM-HOLE PRESSURE
WBHP
-- 'PRODUCER'
/
-- GAS AND OIL SATURATIONS IN INJECTION AND PRODUCTION CELL
BGSAT
10 10 3 /
1 1 1 /
/
BOSAT
10 10 3 /
1 1 1 /
/
-- PRESSURE IN INJECTION AND PRODUCTION CELL
BPR
10 10 3 /
1 1 1 /
/
SCHEDULE ===============================================================
-------- THE SCHEDULE SECTION DEFINES THE OPERATIONS TO BE SIMULATED
------------------------------------------------------------------------
-- CONTROLS ON OUTPUT AT EACH REPORT TIME
RPTSCHED
1 0 1 1 0 0 4 2 2 0 0 2 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 /
--IMPES
-- 1.0 1.0 10000.0 /
-- SET 'NO RESOLUTION' OPTION
--DRSDT
-- 0 /
-- SET INITIAL TIME STEP TO 1 DAY AND MAXIMUM TO 6 MONTHS
TUNING
1 182.5 /
1.0 0.5 1.0E-6 /
/
-- WELL SPECIFICATION DATA
--
-- WELL GROUP LOCATION BHP PI
-- NAME NAME I J DEPTH DEFN
WELSPECS
'PRODUCER' 'G' 10 10 8400 'OIL' /
/
-- COMPLETION SPECIFICATION DATA
--
-- WELL -LOCATION- OPEN/ SAT CONN WELL
-- NAME I J K1 K2 SHUT TAB FACT DIAM
COMPDAT
'PRODUCER' 10 10 1 1 'OPEN' 0 -1 0.5 /
/
-- PRODUCTION WELL CONTROLS
--
-- WELL OPEN/ CNTL OIL WATER GAS LIQU RES BHP
-- NAME SHUT MODE RATE RATE RATE RATE RATE
WCONPROD
'PRODUCER' 'OPEN' 'ORAT' 20000 4* 1000 /
/
-- WCONPROD
-- 'PRODUCER' 'OPEN' 'BHP' 5* 1000 /
-- /
-- INJECTION WELL CONTROLS
--
-- WELL INJ OPEN/ CNTL FLOW
-- NAME TYPE SHUT MODE RATE
-- WCONINJ
-- YEAR 1
TSTEP
--0.2343 0.1393 0.1840 0.2189 0.2235
--1.0 4.0 2*5.0 6*5.0 10*25.0
1.0 14.0 13*25.0
/
TSTEP
25.0
/
-- YEAR 2
TSTEP
13*20.0 7*13.0
/
TSTEP
14.0
/
-- YEAR 3
TSTEP
17*10.0
/
TSTEP
12.5
/
-- 912.50 --> 1000.0
TSTEP
8.5 16*5.0
/
TSTEP
5.0
/
-- 1000.0 --> 1100.0
TSTEP
19*5.0
/
TSTEP
5.0
/
-- 1100.0 --> 1200.0
TSTEP
19*5.0
/
TSTEP
5.0
/
TSTEP
10.0 /
END ================================================================