!CO2 Storage - symmetrical cold injection

SIMULATION
  SIMULATION_TYPE SUBSURFACE
  PROCESS_MODELS
    SUBSURFACE_FLOW Flow
      MODE GAS_WATER
      OPTIONS
        RESERVOIR_DEFAULTS
      /
    / ! end of subsurface_flow
  / ! end of process models
END  !! end simulation block

SUBSURFACE

#=========================== discretization ===================================

GRID
  TYPE grdecl inj_th.grdecl
END


#=========================== times ============================================

TIME
  !uses default start date JAN 2000
  FINAL_TIME 7 y
  INITIAL_TIMESTEP_SIZE 0.1 d
  MAXIMUM_TIMESTEP_SIZE 100 d at 0 d
  MAXIMUM_TIMESTEP_SIZE 365 d at 3000 d
  MAXIMUM_TIMESTEP_SIZE 1000 d at 18250 d
END

#=========================== output options ===================================

OUTPUT
  MASS_BALANCE_FILE
   PERIODIC TIMESTEP 1  
  END
  ECLIPSE_FILE
    TIMES y 1 2 3 4 5 6 7
    PERIOD_SUM TIMESTEP 1
    OUTFILE
    WRITE_DENSITY
    WRITE_VISCOSITY
    SUMMARY_D
      BPR   1  1  10
      BGSAT 1  1  10
      BVGAS 1  1  10
      BVWAT 1  1  10
      BDENG 1  1  10
      BDENW 1  1  10
      BTEMP 1  1  10
      BGSAT 7  7   4
      BPR   7  7   4
      BTEMP 7  7   4
    END_SUMMARY_D
  END
  LINEREPT
END

#=========================== material properties ==============================

MATERIAL_PROPERTY formation
  ID 1
  ROCK_DENSITY 2.350d3          kg/m^3
  SPECIFIC_HEAT 1.0d3           J/kg-C
  THERMAL_CONDUCTIVITY_DRY 1.6  W/m-C
  THERMAL_CONDUCTIVITY_WET 4.3  W/m-C
  ROCK_COMPRESSIBILITY    4.35d-5 1/Bar
  ROCK_REFERENCE_PRESSURE 1.0     Bar
  CHARACTERISTIC_CURVES ch1
/

#=========================== saturation functions =============================

CHARACTERISTIC_CURVES ch1
 
 TABLE swfn_table
   PRESSURE_UNITS Bar
   external_file SWFN_drain.dat
 END

 TABLE sgfn_table
  PRESSURE_UNITS Bar
  external_file SGFN_drain.dat
 END

/

#=========================== fluid properties =================================

FLUID_PROPERTY
  PHASE LIQUID
  DIFFUSION_COEFFICIENT 2.0d-9
/
FLUID_PROPERTY
  PHASE GAS
  DIFFUSION_COEFFICIENT 2.0d-5
/

BRINE 1.0 MOLAL

EOS WATER
 SURFACE_DENSITY 1000.0 kg/m^3
END

EOS GAS
  SURFACE_DENSITY 1.8 kg/m^3
  CO2_DATABASE co2_dbase.dat
END

#=========================== Equilibration ==================================

EQUILIBRATION

  PRESSURE 120.0 Bar
  DATUM_D  1200 m

  GAS_IN_LIQUID_MOLE_FRACTION 0.0

  TEMPERATURE_TABLE
    D_UNITS m
    TEMPERATURE_UNITS C !cannot be otherwise
    RTEMPVD
     1150 50
     1300 54.5
    /
  END
 
/

#=========================== thermal aquifers =================================

TBC_DATA tbctop
  TEMPVD m C
     1150 50
     1300 54.5
  END TEMPVD
  CONN_D 1 20 1 20 1 1 Z-     !top faces 
END

TBC_DATA tbcbot
  TEMPVD m C
     1150 50
     1300 54.5
  END TEMPVD
  CONN_D 1 20 1 20 16 16 Z+ ! bottom faces
END

AQUIFER_DATA north
  BACKFLOW
  DEPTH     1225.0 m
  THICKNESS  50 m
  RADIUS    3000.0 m
  PERM      300  mD
  TEMPERATURE 53.0 C
  COMPRESSIBILITY 2.0E-9 1/Pa
  POROSITY  0.3
  VISCOSITY 1 cP
  CONN_D 1  20 20 20 4 13 Y+
END

AQUIFER_DATA east
  BACKFLOW
  DEPTH     1225.0 m
  THICKNESS  50 m
  RADIUS    3000.0 m
  PERM      300  mD
  TEMPERATURE 53.0 C
  COMPRESSIBILITY 2.0E-9 1/Pa
  POROSITY  0.3
  VISCOSITY 1 cP
  CONN_D 20 20  1 20 4 13 X+
END



#=========================== Wells ==================================

WELL_DATA injg
  WELL_TYPE GAS_INJECTOR
  INJECTION_ENTHALPY_P  138 Bar
  INJECTION_ENTHALPY_T  15 C 
  THETA_FRACTION 0.25
  BHPL     1000    Bar
  TARG_GM  0.25 Mt/y
  CIJK_D   1 1 9 13 
  TIME 7 y
  SHUT
END

#================================================================================

END_SUBSURFACE