################################################################################ RADIAL FLOW FROM A CO2 INJECTION WELL The problem is from the paper: Pruess, K. et al. 2004. Code Intercomparison Builds Confidence in Numerical Simulation Models for Geologic Disposal of CO2. Energy, 29(9-10): 1431-1444. DOI:10.1016/j.energy.2004.03.077. Pruess, K. et al. 2002. Multiphase flow dynamics during CO2 injection into saline aquifers// Envirom. Geol. 42, 282-295. DOI:10.1007/s00254-001-0498-3 See description at: http://dx.doi.org/10.1016/j.energy.2004.03.077; http://dx.doi.org/10.1007/s00254-001-0498-3 http://esd1.lbl.gov/files/research/projects/tough/ documentation/TOUGH2-ECO2N_V2.0_Users_Guide.pdf -------------------------------------------------------------------------------- Last update on 12.11.2015 Program Version: 2015.G01 RUNSPEC ################### RUNSPEC section begins here ###################### METRIC ^^^^^^^^^^^^ We use METRIC units ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ Note, this keywords affects the units in many keywords below. GASSTORE ^^^^^^ We switch the program to the GASSTORE EOS module ^^^^^^^^^^^^^^ We enable NACL - salt component; CO2-SEQ - CO2 component; WATER - water component; DISGAS - CO2 dissolution in brine; VAPWAT - water vapour in the gas phase. CAPPRES We enable capillary pressure. ROCKCOMP We enable rock compressibility. ISOTHERM We use the Isothermal modelling option. FAST We use FAST option. AUXARRAY We introduce an auxiliary array #SIMVAR / #SIMVAR, where we will store values of / the similarity variable. GRID ##################### GRID section begins here ####################### The grid is specified within brackets MAKE-ENDMAKE MAKE <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< -- cartesian We select Radial gridding -- grid nr ntheta nz option and specify the number of RADIAL 1270 1 1 / grid blocks along every direction. RTZBOUND 0 6000 -0.0314159 0.0314159 0 100 / We specify the domain extent in meters along R and Z axes, and in radians along Theta direction. Thus, we consider the 0.01 fraction of the full circle. DRV 1100*1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 20 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 200 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 2000 2000 3000 4000 5000 6000 7000 8000 9000 10000 10*10000 10*20000 10*30000 10*40000 10*50000 10*60000 10*70000 10*80000 10*90000 10*100000 / Using the DRV keyword we specify the grid blocks extensions along axis R. We have a uniform grid for the first 1100 meters, and farther the extension of the grid blocks increases. SRCSPECG We introduce a point source INJECTOR 'INJECTOR' 1 1 1 / in the grid block i=1, j=1, k=1. / ENDMAKE >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> EQUALS We equal PORO 0.12 / porosity to 0.12; PERMI 100 / permeablity to 100 mD. / MULTIPLY Because we specified grid only for CELLVOL 100.0 / 0.01 fraction of the full circle we AREA 100.0 / multiply cell volumes and interface / areas (both in R and Z directions) by the factor 100. PROPS ###################################################################### INCLUDE We load the properties of CO2 from 'CO2TAB.METRIC.INC' / include file. Rock properties are specified within brackets ROCK-ENDROCK ROCK <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< / ROCKECL -- -refpres- --compr-- We specify the reference pressure (bar) 120 4.5E-5 / and the rock compressibility (1/bar). ENDROCK >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The relative permeabilities and capillary pressure are specified within brackets SAT-ENDSAT SAT <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< / SGLF -- sgas krgas krliq pcgl 0.00 0.0 1.0 0.0 / We define the relative 0.01 0.0 0.6284139343 0.0248270428 / permeability of the 0.025 0.0 0.4690103837 0.0413776570 / gas phase (krgas), the 0.05 0.0 0.3243551786 0.0616685866 / relative permeability 0.1 0.0008753194 0.1754738386 0.0942244008 / of the liquid (krliq) 0.2 0.0217429361 0.0560863397 0.1525231868 / and gas-liquid 0.3 0.0919085466 0.0163414095 0.2147871648 / capillary pressure 0.4 0.2281782851 0.0036709332 0.2901918260 / (pcgl) as the function 0.5 0.4339133784 0.0004812265 0.3906042421 / of the gas saturation 0.6 0.6990301461 0.0000159286 0.5385130956 / (sgas). 0.7 1.0 0.0 0.7873987963 / / SSFN -- ssol krf We define the reduction factor for 0.0 1.0 / the rock permeability (krf) as the 0.01 0.9 / the faction of the solid phase 0.05 0.7 / saturation (ssol). 0.1 0.5 / 0.2 0.0 / / ENDSAT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> INIT ####################### INIT section begins here ##################### EQUALS We equal in the whole domain: PRES 120 / - pressure to 120 bar; TEMPC 45 / - temperature to 45 C; SLIQ 1.0 / - saturation of the liquid phase to 1; XSM 0.15 / - salt mass fraction in brine to 0.15; XGM 0.0 / - gas mass fraction in brine to 0. / EQUALS We fix parameters in the grid block ACTNUM 2 1270 / i=1270, j=1, k=1. (The farthest block / from the well. We specify parameters in the stock tank EQUALNAM associated with the INJECTOR source: PRES 120 'INJECTOR' / - pressure to 120 bar; TEMPC 45 / - temperature to 45 C; SGAS 1.0 / - saturation of the gas phase to 1; YWM 0.0 / - water vapour mass fraction in the / gas phase to 0. Here, we calculate the values of the similarity variable corresponding to OPERATE the time 1000 days (8.64e7 seconds): #SIMVAR XCOORD COPY / - #SIMVAR:=XCOORD; #SIMVAR #SIMVAR MULTARR / - #SIMVAR:=#SIMVAR*#SIMVAR; #SIMVAR 1* DIVIDE 8.64E7 / - #SIMVAR:=#SIMVAR/8.64E7. / RPTSUM #SIMVAR PRES TEMPC XGM YWM SLIQ SGAS XSMR SSOL ASCII/ We report these properties in the grid blocks; RPTSRC We report these properties for the SMIR#1 SMIR#2 SMIT#1 SMIT#2 / INJECTOR. The 1st component is WATER, the 2nd component is CO2. SCHEDULE #################### SCHEDULE section begins here ################### REPORTS No reports in the LOG-file. NOTHING / VARS The maximal alteration of pressure PRES DMAX 10 / in a Newton iteration is 10 bars. / SRCINJE We specify that the injecton rate is 'INJECTOR' MASS 1* 400 1* 8640 / 8640 ton/day=100 kg/sec, and the / maximum pressure in the grid block i=1, j=1, k=1 is 400 bars. TUNING The maximun and next time step is 4*0.005 / 0.005 days. TIME We advance simulation to t=0.3 days. 0.3 / TUNING We increase the maximum time step to 1* 0.05 / 0.05 days. TIME We advance simulation to t=3 days. 1 3 / TUNING We increase the maximum time step to 1* 0.5 / 0.5 days. TIME We advanxe simulation to t=20 days. 10 20 / TUNING We increase the maximum time step to 1* 2 / 2 days. TIME We advance simulation to 300 days. 50 100 300 / TUNING We increase the maximum time step to 1* 5 / 5 days. TIME We advance simulation to 1000 days. 500 1000 / POST ####################### POST section begins here ##################### CONVERT We convert the output to ParaView compatible format. RPTPOST We define the properties output from the TIME SMIR#1 SMIR#2 SMIT#1 SMIT#2 / following POSTSRC keyword. POSTSRC We save consolidated time series INJECTOR / report for the INJECTOR. / END ######################################################################