################################################################################ TEST FOR LONGITUDIAL MECHANICAL DISPERSION OF SOLUTES The problem is from the paper Afanasyev, A.A. 2018. Numerical modelling of solute flow dispersion in porous media using simulator MUFITS // IOP Conf. Series: Journal of Physics: Conf. Series. 1129: 012002 DOI:10.1088/1742-6596/1129/1/012002. See description at: http://doi.org/10.1088/1742-6596/1129/1/012002. -------------------------------------------------------------------------------- Last update on 17.07.2019 Program Version: 2019.B01 ################################################################################ RUNSPEC ################### RUNSPEC section begins here ###################### GASSTORE ^^^^^^ We switch the program to the GASSTORE EOS module ^^^^^^^^^^^^^^ The fluid consists of WATER - water SALT - salt (solute) ISOTHERM We switch the program to the isothermal modelling option. METRIC ^^^^^^^^^^^^^^^^^^^^^ We use METRIC units ^^^^^^^^^^^^^^^^^^^^^^^^^^^^ MECHDISP We enable the mechanical dispersion modelling option. FAST We use FAST option. AUXARRAY We create auxiliary array #XCOORD which #XCOORD / will be used for reporting x coordinate / from the SCHEDULE section. 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 12000 1 1 / along every axis. XYZBOUND -- xmin-xmax ymin-ymax zmin-zmax We specify the domain extent in metres. -20 100 0 1 0 1 / The domain length is 120 metres. BOUNDARY -- ---- ----------- ------------ -------- -------- 102 6* 'I+' 5* INFTHIN 4* 1 2 / -- ---- ----------- ------------ -------- -------- / We create infinetely thin gost grid block at the right boundary. The parameters in this block are fixed (ACTNUM=2), will be used as open boundary. Further we can refer to this block with FLUXNUM=102. SRCSPECG -- name i j k MYSRC 1 1 1 / We create a point source with the name / MYSRC and place it in the grid block i=1, j=1 and k=1. ENDMAKE >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> EQUALS We equate PORO 0.1 / - porosity to 0.1; PERMX 100 / - X permeability to 100 mD; PERMY 100 / - Y permeability to 100 mD. / RPTGRID We reporrt these properties from the DEPTH PORO PERMX PERMZ / GRID section. PROPS ####################### PROPS section begins here #################### Rock properties are specified within brackets ROCK-ENDROCK ROCK <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< / ROCKDH We specify that rock density and heat 2896.5 1.0 / capacity are 2896.5 kg/m3 and 1 kJ/kg/K. These data are irrelevant with the ISOTHERM simulation but must be specified according to input data rules. MDISPL This keyword specifies the longitudial -- phase -- sat -- dispersivity for the LIQ phase. The LIQ 0 1 / dispersivity is equal to the true fluid -- vel ------------- velocity and it does not depend on the 0 0 0 / phase saturation. 1 1 1 / / ENDROCK >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The relative permeabilities are specified within brackets SAT-ENDSAT SAT <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< / SSFN -- ssol krf We define the reduction factor for 0.0 1.0 / the rock permeability (krf). These data 1.0 0.0 / does not influence the simulation / results because the halite does not precipitate. ENDSAT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> INIT ####################### INIT section begins here ##################### REGALL ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ The following OPERAREG keyword will be applied to all grid blocks. OPERAREG -- -------- -------- -------- -------- -------- -------- -------- -------- PRES 1* ROCKNUM 1 EQUALS 50 / TEMPC 1* ROCKNUM 1 EQUALS 50 / XSM 1* FLUXNUM 1 EQUALS 0.0000000 / XSM XCOORD FLUXNUM 1 SETINT 0.0001 -3000 0 / XSM 1* FLUXNUM 102 EQUALS 0.0000000 / -- -------- -------- -------- -------- -------- -------- -------- -------- / With the above keyword we specfy that at t=0: - pressure = 50 bar in the whole domain - temperature = 50 C in the whole domain - salt concentration is equal to 0 in all internal grid blocks - salt concentration is equal to 0.0001 at x<0 - salt concentration is equal to 0 at the right boundary. EQUALNAM The fluid injected throuh the point source is PRES 50 MYSRC / as follows: T=50 C at P=50bar, its salinity TEMPC 50 / is 0.0001. XSM 0.0001 / / COPY We copy x coordinate to the user created array XCOORD #XCOORD / #XCOORD. / RPTSUM PRES TEMPC XSM PHST DENT SSOL #XCOORD / We report these parameters at every report time. SCHEDULE #################### SCHEDULE section begins here ################### REPORTS NOTHING / No reports to LOG-file. Linear solver setup: ILUTFILL - Ilut-fill parameter is 5; 5 / ILUTDROP - Ilut-drop parameter is 1e-4. 1E-4 / TUNING The maximum time step is 7.5 days. 1* 7.5 / SRCINJE Through the point source we inject MYSRC RATEIN 1* 100 1E-4 / fluid at a rate of 1e-4 m3/day. / TSTEP We advance simulation to 50000 days 100*500 / reporting results every 500 days. POST ######################## POST section begins here #################### CONVERT We convert the output to ParaView compatible format. MNEMPARA The change the output format of the XSM XSM 1PE12.5 / mnemonic. / RPTPOST The report list consits of two columns, #XCOORD XSM / x xoordinate and salinity. POSTSPEC With these kewords we report the 1* TIME25.CSV 25 / salinity distribution at the 25th report POSTBLK time, i.e. at t=12500 days. 1750 1 1 / 2000 1 1 / 2250 1 1 / 2500 1 1 / 2750 1 1 / 3000 1 1 / 3250 1 1 / 3500 1 1 / 3750 1 1 / 4000 1 1 / 4250 1 1 / 4500 1 1 / 4750 1 1 / 5000 1 1 / 5250 1 1 / 5500 1 1 / 5750 1 1 / 6000 1 1 / 6250 1 1 / / POSTSPEC With these kewords we report the 1* TIME50.CSV 50 / salinity distribution at the 50th report POSTBLK time, i.e. at t=25000 days. 2750 1 1 / 3000 1 1 / 3250 1 1 / 3500 1 1 / 3750 1 1 / 4000 1 1 / 4250 1 1 / 4500 1 1 / 4750 1 1 / 5000 1 1 / 5250 1 1 / 5500 1 1 / 5750 1 1 / 6000 1 1 / 6250 1 1 / 6500 1 1 / 6750 1 1 / 7000 1 1 / 7250 1 1 / / POSTSPEC With these kewords we report the 1* TIME75.CSV 75 / salinity distribution at the 75th report POSTBLK time, i.e. at t=37500 days. 3750 1 1 / 4000 1 1 / 4250 1 1 / 4500 1 1 / 4750 1 1 / 5000 1 1 / 5250 1 1 / 5500 1 1 / 5750 1 1 / 6000 1 1 / 6250 1 1 / 6500 1 1 / 6750 1 1 / 7000 1 1 / 7250 1 1 / 7500 1 1 / 7750 1 1 / 8000 1 1 / 8250 1 1 / / POSTSPEC With these kewords we report the 1* TIME100.CSV 100 / salinity distribution at the 100th POSTBLK report time, i.e. at t=50000 days. 4750 1 1 / 5000 1 1 / 5250 1 1 / 5500 1 1 / 5750 1 1 / 6000 1 1 / 6250 1 1 / 6500 1 1 / 6750 1 1 / 7000 1 1 / 7250 1 1 / 7500 1 1 / 7750 1 1 / 8000 1 1 / 8250 1 1 / 8500 1 1 / 8750 1 1 / 9000 1 1 / 9250 1 1 / / END #######################################################################