################################################################################ 7TH SPE COMPARATIVE STUDY (CASE 4B) The problem is from SPE-21221 paper (case 4b): Nghiem, L. et al. 1991. Seventh SPE Comparative Solution Project: Modelling of Horizontal Wells in Reservoir Simulation. SPE Symposium on Reservoir Simulation, 17-20 February, Anaheim, California. DOI: 10.2118/21221-MS. See description at: http://dx.doi.org/10.2118/21221-MS -------------------------------------------------------------------------------- Last update on 10.01.2017 Program Version: 2016.C01 ################################################################################ Note: We type all comments in lowercase letters while all keywords and mnemonics must be in uppercase letters. -------------------------------------------------------------------------------- Note: -- any line beginning with '!' or '--' is a comment line; any line not beginning with an uppercase letter outside keyword ! instruction is also a comment line. -------------------------------------------------------------------------------- Note: Any data line must be terminated by the slash '/' sign. -------------------------------------------------------------------------------- Note: We denote by the repeated symbols the following structural elements of this RUN-file: '##############' - delimits the sections of the RUN-file. '<<<<<<<<<<<<<<' - designates an opening bracket. '>>>>>>>>>>>>>>' - designates a closing bracket. '^^^^^^^^^^^^^^' - designates the keyword, which affects the order of the data input for all the following keywords (both in the current and in the following sections). '==============' - merge the keywords and the comments in paragraphs. -------------------------------------------------------------------------------- RUNSPEC ################### RUNSPEC section begins here ###################### BLACKOIL ^^^^^^ We switch the program to the BLACKOIL EOS module ^^^^^^^^^^^^^^ We enable: WATER - water phase; OIL - oil phase; GAS - gas phase; DISGAS - gas dissolution in oil phase. FIELD ^^^^^^^^^^^^ We use FIELD units ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ Note, this keywords affects the units in many keywords below. We enable: ROCKCOMP - rock compressibility; CAPPRES - capillary pressure; MUSEGWEL - multisegmented wells. OPTIONS Piecewise linear interpolation 53* 1 / of relative permeabilities. Use this option if you want a better fit to the SPE-21221 paper. FAST We use FAST option. GRID ##################### GRID section begins here ####################### The grid is specified within brackets MAKE-ENDMAKE MAKE <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< -- cartesian We select Cartesian gridding -- grid nx ny nz option and specify the number of CART 9 9 6 / grid blocks along every axis. XYZBOUND -- xmin-xmax ymin-ymax zmin-zmax We specify the domain extent in feet. 0 1 0 1 3590 3591 / The top of reservoir is at depth 3590 feet. DXV The grid blocks length along X axis 9*300 / is 300 feet. DYV Here, we specify the grid blocks 620 400 200 100 60 100 200 400 620 / length along Y axis. DZV Here, we specify the thickness of the 4*20 30 50 / layers. WELSPECS -- -------- -- ---- ---- -------- -- --- -- ------ -- well i- j- datum mod flux -- name ind ind depth el num -- -------- -- ---- ---- -------- -- --- -- ------ Here, we specify INJECTOR 1* 1 5 1* 6* 1* 4* 222 / both wells PRODUCER 1* 8 5 1* 6* SEG 4* 111 / -- -------- -- ---- ---- -------- -- --- -- ------ / The AVG well model (default option) is used for the INJECTOR. and the SEG model is used for the PRODUCER. The pipe segment associated with the INJECTOR can be farther referred to by using FLUXNUM=222, and the pipe segments asscoiated with the PRODUCER can be reffered to by FLUXNUM=111. COMPDAT -- -------- ---- ---- ---- ---- ---- ------ -------- -------- -- --- -- ----- -- well i- j- k- k- mode satnum well well dir rough -- name ind ind min max index diameter ness -- -------- ---- ---- ---- ---- ---- ------ -------- -------- -- --- -- ----- INJECTOR 1 5 6 6 OPEN 1* 243.457 / 243.457 cP*bbl/day/psi INJECTOR 2 5 6 6 OPEN 1* 243.457 / =2.16e5 ft*mD INJECTOR 3 5 6 6 OPEN 1* 243.457 / INJECTOR 4 5 6 6 OPEN 1* 243.457 / INJECTOR 5 5 6 6 OPEN 1* 243.457 / INJECTOR 6 5 6 6 OPEN 1* 243.457 / INJECTOR 7 5 6 6 OPEN 1* 243.457 / INJECTOR 8 5 6 6 OPEN 1* 243.457 / INJECTOR 9 5 6 6 OPEN 1* 243.457 / -- -------- ---- ---- ---- ---- ---- ------ -------- -------- -- --- -- ----- PRODUCER 8 5 1 1 OPEN 1* 1* 0.375 3* 'X' 1* 0.001 / PRODUCER 7 5 1 1 OPEN 1* 1* 0.375 3* 'X' 1* 0.001 / PRODUCER 6 5 1 1 OPEN 1* 1* 0.375 3* 'X' 1* 0.001 / PRODUCER 5 5 1 1 OPEN 1* 1* 0.375 3* 'X' 1* 0.001 / PRODUCER 4 5 1 1 OPEN 1* 1* 0.375 3* 'X' 1* 0.001 / PRODUCER 3 5 1 1 OPEN 1* 1* 0.375 3* 'X' 1* 0.001 / PRODUCER 2 5 1 1 OPEN 1* 1* 0.375 3* 'X' 1* 0.001 / -- -------- ---- ---- ---- ---- ---- ------ -------- -------- -- --- -- ----- / Here, we specify the INJECTOR completions in the bottom layer. The well index is specified explicitly. And we specify the PRODUCER completions in the top layer. DZWELL This is for visualization. 300 30 / ENDMAKE >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> EQUALS -- -------- -------- We equal PORO 0.2 / - porosity to 0.2; PERMX 3000 / - X-axis permeability to 300; PERMY 3000 / - Y-axis permeability to 300; PERMZ 300 / - Z-axis permeability to 30. -- -------- -------- / RPTGRID This properties are reported from the DEPTH PORO PERMX PERMY PERMZ / section GRID. PROPS ####################### PROPS section begins here #################### Rock properties are specified within brackets ROCK-ENDROCK ROCK <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< / ROCKECL We specify the reference pressure (psia) -- -refpres- --compr-- and the rock compressibility (1/psi). 3600 4.0D-6 / ENDROCK >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The relative permeabilities are specified within brackets SAT-ENDSAT SAT <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< / STONE2 We use the 2nd Stone formula for the oil relative permeability in the three-phase-flow regions. SWOF -- swat krw krow pcow 0.22 0 1 6.30 / Here, we specify the water relative 0.3 0.07 0.4 3.60 / permeability (krw), oil relative 0.4 0.15 0.125 2.70 / permeability for the oil-water flows 0.5 0.24 0.0649 2.25 / (krow) and oil-water capillary 0.6 0.33 0.0048 1.80 / pressure as the functions of water 0.8 0.65 0 0.90 / saturation (swat). 0.9 0.83 0 0.45 / 1 1 0 0 / / SGOF -- sgas krg krog pcog 0 0 1 0 / Here, we specify the gas relative 0.04 0 0.6 0.2 / permeability (krg), oil relative 0.1 0.022 0.33 0.5 / permeability for the oil-gas flows 0.2 0.1 0.1 1.0 / (krg) and gas-oil capillary pressure 0.3 0.24 0.02 1.5 / as the functions of gas saturation 0.4 0.34 0 2.0 / (sgas). 0.5 0.42 0 2.5 / 0.6 0.5 0 3.0 / 0.7 0.8125 0 3.5 / 0.78 1 0 3.9 / / ENDSAT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The fluid properties are specified within brackets EOS-ENDEOS EOS <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< / DENSITY We specify the phases densities at -- oil water gas surface conditions (lb/cft). 45 62.14 0.0702 / PVTW We specify water properties: reference -- refpres bwat comprw viswat pressure (refpres), form. vol. factor at 14.7 1.0142 3.0D-6 0.96 / refpres (bwats), viscosity (viswat) and water compressibility (comprw). PVDG -- pgas bgas visgas We specify the properies of dry gas, 400 5.9 0.013 / i.e. the formation volume factor (bgas), 800 2.95 0.0135 / and viscosity (visgas) on pressure 1200 1.96 0.014 / (pgas). 1600 1.47 0.0145 / 2000 1.18 0.015 / 2400 0.98 0.0155 / 2800 0.84 0.016 / 3200 0.74 0.0165 / 3600 0.65 0.017 / 4000 0.59 0.0175 / 4400 0.54 0.018 / 4800 0.49 0.0185 / 5200 0.45 0.019 / 5600 0.42 0.0195 / / PVTO -- rs poil boil visoil We specify the properties of live oil, 0.165 400 1.012 1.17 / i.e. the oil pressure (poil), 0.335 800 1.0255 1.14 / formation volume factor (boil), 0.500 1200 1.038 1.11 / viscosity (visoil) on gas-oil ratio 0.665 1600 1.051 1.08 / of the solution (rs). 0.828 2000 1.063 1.06 / 0.985 2400 1.075 1.03 / 1.130 2800 1.087 1.00 / 1.270 3200 1.0985 0.98 / 1.390 3600 1.11 0.95 / 1.500 4000 1.12 0.94 / 1.600 4400 1.13 0.92 / 1.676 4800 1.14 0.91 / 1.750 5200 1.148 0.9 / 1.810 5600 1.155 0.89 6000 1.1504 0.89 6400 1.1458 0.89 6800 1.1412 0.89 7200 1.1367 0.89 / / ENDEOS >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> INIT ####################### INIT section begins here ##################### We specify pressure and water saturation EQUALS for every layer: PRES 3600 4* 2*1 / - first (top) layer; SWAT 0.289 / - first (top) layer; PRES 3608 4* 2*2 / - second layer; SWAT 0.348 / - second layer; PRES 3616 4* 2*3 / etc. SWAT 0.473 / PRES 3623 4* 2*4 / SWAT 0.649 / PRES 3633 4* 2*5 / SWAT 0.869 / PRES 3650 4* 2*6 / SWAT 1.0 / / COPY We copy pressure array to the buble PRES PBUB / pressure array. / REGALL This keyword affects the following keyword EQUALREG. We specify initial parameters of the EQUALREG fluid in the boreholes. PRES 3600 FLUXNUM 111 / -- in the PRODUCER. SWAT 0.289 / PBUB 3600 / PRES 3650 FLUXNUM 222 / -- in the INJECTOR. PBUB 3650 / SWAT 1.0 / / RPTSUM This parameters are saved at every PRES SWAT SOIL SGAS PBUB RS / report time. RPTWELL WBHP WWIR WOPR WWPR WGPR WLPR WOPT WWPT WGPT WLPT WWOR WGOR/ SCHEDULE ################### SCHEDULE section begins here #################### REPORTS No reports in the LOG-file. NOTHING / WEEKTOL We use the WEEKTOL option. WELLPROD -- -------- ---- ---- -------- -------- -------- -------- -------- -------- -- name mode ctrl volume mass bhp max pump tanktype vol rate -- rate rate rate dimens -- -------- ---- ---- -------- -------- -------- -------- -------- -------- PRODUCER OPEN LRAT 9000 1* 1500 1E5 / -- -------- ---- ---- -------- -------- -------- -------- -------- -------- / We produce 9000 bbl/day of liquid (oil and water) at the minimum bottom-hole pressure 1500 psi. WELLINJE -- -------- ---- ---- -------- -------- -------- -------- -------- -------- -- name mode ctrl volume mass bhp max pump tanktype vol rate -- rate rate rate dimens -- -------- ---- ---- -------- -------- -------- -------- -------- -------- INJECTOR OPEN WRAT 6000. 1* 4000 1E5 WATER / -- -------- ---- ---- -------- -------- -------- -------- -------- -------- / We inject 6000 bbl/day of water at the maximum bottom-hole pressure 4000 psi. VARS Maximimum pressure alteration at PRES DMAX 500 / a time step is 500 psi. / TUNING The next timestep is limited to 0.1 25 / 0.1 days. The maximum timestep is 25 days. TSTEP We simulate the production over 15*100 / 1500 days reporting every 100 days. POST ####################### POST section begins here ##################### RPTPOST We save time series report for TIME WBHP WWIR / the INJECTOR. POSTWELL INJECTOR / / RPTPOST NOTHING TIME WBHP WOPR WOPT WWPT WWOR WLPR WGPT WWPR WGOR WGPR / We save time series report for the POSTWELL PRODUCER. PRODUCER / / RPTPOST We report pressure in the nodes NOTHING TIME PRES / 489 and 495. These nodes are the first POSTCELL and the last nodes of the PRODUCER. This 489 / can be seen through formatted output of 495 / any summary. The difference between / these pressures is the pressure drop along in the PRODUCER. CONVERT We convert the output to ParaView compatible format. END #####################################################################