add a 300k cells version of the SPE1 problem

i.e., it provides a simple deck which is suitable for performance
comparisons.
This commit is contained in:
Andreas Lauser 2015-03-02 19:33:34 +01:00
parent 2c8e9ae71b
commit b472253edd
2 changed files with 482 additions and 0 deletions

464
spe1_large/LARGE_SPE1.DATA Normal file
View File

@ -0,0 +1,464 @@
-- =====================================================================
-- 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
100 100 30 /
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
100*100.0
/
DYV
100*100.0
/
DZV
10*2.0 10*3.0 10*5.0
/
PORO
300000*0.3
/
PERMX
100000*500.0
100000*50.0
100000*200.0
/
PERMY
100000*500.0
100000*50.0
100000*200.0
/
-- Note: ignoring MULTZ!
-- layer 1-2 'MULTZ' 0.64 /
-- layer 2-3 'MULTZ' 0.265625 /
-- Reducing PERMZ a little instead
PERMZ
100000*300.0
100000*30.0
100000*50.0
/
TOPS
10000*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
-- 100000*3.297832774859256e7
-- 100000*3.302313357125603e7
-- 100000*3.309483500720813e7
-- /
PRESSURE
100000*4783.10205078125
100000*4789.60058593750
100000*4800.00000000000
/
SWAT
300000*0.12
/
SGAS
300000*0.0
/
-- RS
-- 300000*226.1966570852417
-- /
RS
300000*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
95 95 25 /
5 5 5 /
/
BOSAT
95 95 25 /
5 5 5 /
/
-- PRESSURE IN INJECTION AND PRODUCTION CELL
BPR
95 95 25 /
5 5 5 /
/
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' 95 95 8400 'OIL' /
'INJECTOR' 'G' 5 5 8335 'GAS' /
/
-- COMPLETION SPECIFICATION DATA
--
-- WELL -LOCATION- OPEN/ SAT CONN WELL
-- NAME I J K1 K2 SHUT TAB FACT DIAM
COMPDAT
'PRODUCER' 95 95 21 30 'OPEN' 0 -1 0.5 /
'INJECTOR' 5 5 1 10 'OPEN' 1 -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
-- 'INJECTOR' 'GAS' 'OPEN' 'RATE' 100000 /
WCONINJE
'INJECTOR' 'GAS' 'OPEN' 'RATE' 100000 100000 50000/
/
--------------
-- YEAR 1
--------------
TSTEP
-- 1 3 7 12 19 23 27
2*0.05 4*0.1 5*0.2 7*0.5 4*1.0 4*2 2.5
/
-- t=19.5
-- at this point in time, gas around the producer well starts to
-- become significant, so convergence rates are reduced for a while...
TSTEP
-- 28 33 34 35 37
5*0.5 0.15 0.32 2*0.39 12*0.45
/
-- t=28.65
TSTEP
-- 49 60 74 76 78 132
11*1.25 14*2.5 2*5.0 2*2.5 54*5.0 2.6
/
-- YEAR 2
TSTEP
73*5.0
/
-- YEAR 3
TSTEP
146*2.5
/
END ================================================================

18
spe1_large/README Normal file
View File

@ -0,0 +1,18 @@
This directory contains the input data file for the first comparison
project of the society of petroleum engineers [1] with each gridcell
refinded by a factor of 10 into each direction. The dataset consists
of a single file called 'LARGE_SPE1.DATA' and uses the ECL format
which is used by many standard commercial simulators for oil
reservoirs. Its purpose is to allow performance comparisons for a
much simpler deck compared to the "Norne" data set.
The ownership of the copyright of this dataset 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