! SPE5 benchmark, from:
! 'Fifth Comparative Solution Project: Evaluation of Miscible Flood Simulators'
! J.E. Killough and C.A. Kossack, SPE 16000
! Note this models the 4-component Todd-Longstaff version of the problem

SIMULATION
  SIMULATION_TYPE SUBSURFACE
  PROCESS_MODELS
    SUBSURFACE_FLOW Flow
      MODE SOLVENT_TL
      OPTIONS
       ISOTHERMAL
       TL_OMEGA 0.6
       RESERVOIR_DEFAULTS
      /
    / ! end of subsurface_flow
  / ! end of process models
END  !! end simulation block

SUBSURFACE

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

GRID
  TYPE grdecl spe5.grdecl
END

#=========================== times ============================================
TIME
  FINAL_TIME 20 y
  INITIAL_TIMESTEP_SIZE 1.0d-5 d
  MINIMUM_TIMESTEP_SIZE 1.0D-10 d
  MAXIMUM_TIMESTEP_SIZE 5 d at 0. d
END

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

OUTPUT
  MASS_BALANCE_FILE
    PERIODIC TIMESTEP 1
  END
  ECLIPSE_FILE
    PERIOD_SUM TIMESTEP 5 d
    PERIOD_RST TIME 200 d
  END
  LINEREPT
END

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

MATERIAL_PROPERTY formation
  ID 1
  ROCK_COMPRESSIBILITY    7.25d-5 1/Bar
  ROCK_REFERENCE_PRESSURE 1.0     Bar
  CHARACTERISTIC_CURVES ch1
/

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

CHARACTERISTIC_CURVES ch1

 TABLE swfn_table
   PRESSURE_UNITS Pa
   SWFN
     0.2    0      3.1026
     0.2809 0.0022 1.3120
     0.3778 0.0180 0.6943
     0.4667 0.0607 0.3378
     0.5556 0.1438 0.1241
     0.6444 0.2809 0.0344 
     0.7000 0.4089 0.0034
     0.7333 0.4855 0.0007
     0.8222 0.7709 0.0
     0.9111 1.0    0.0
     1.0    1.0    0.0
   /
 END

 TABLE sgfn_table
  PRESSURE_UNITS Pa
    SGFN
     0.0     0     0
     0.05    0     0
     0.0889  0.001 0
     0.1778  0.01  0
     0.2667  0.03  0.000069
     0.3556  0.05  0.000069
     0.4445  0.10  0.002068
     0.5333  0.20  0.055158
     0.6222  0.35  0.206843
     0.65    0.39  0.275790
     0.7111  0.56  0.551581
     0.8     1.0   2.068427
    /
 END

 TABLE sof3_table
   SOF3
     0.0    0.0     0.0
     0.0889 0.0     0.0
     0.1778 0.0     0.0
     0.2667 0.0     0.0110
     0.3000 0.0     0.0370
     0.3556 0.0123  0.0878 
     0.4444 0.0835  0.1715
     0.5333 0.2178  0.2963 
     0.6221 0.4153  0.4705
     0.7191 0.6769  0.7023
     0.8    1.0     1.0
   /
 END
 !TEST 
/

#=========================== EOSs =============================================

EOS WATER
  SURFACE_DENSITY 995.5
  DENSITY CONSTANT 995.5
  VISCOSITY CONSTANT 0.0007 ! 0.7 cp
END

EOS OIL
  FORMULA_WEIGHT 800.0d0
  SURFACE_DENSITY 617.19139 kg/m^3
  PVCO
    DATA_UNITS ! Metric in the Eclipse sense
      PRESSURE Bar           !Bar is default
      RS  m^3/m^3            !m^3/m^3 is default
      FVF m^3/m^3            !m^3/m^3 is default
      VISCOSITY cP           !cP default
      COMPRESSIBILITY 1/Bar  !1/Bar is default
      VISCOSIBILITY   1/Bar  !1/Bar is default
    END
    DATA
      TEMPERATURE 15.0
          1.0135     0	  1.034	 0.310	2.92780E-04  0.00211302
         34.473	  20.945  1.1017 0.295	2.92780E-04  0.00211302
         68.947	  39.646  1.1478 0.274	2.92780E-04  0.00211302
         82.737	  47.679  1.1677 0.264	2.92780E-04  0.00211302
        103.421	  60.805  1.1997 0.249	2.92780E-04  0.00211302
        124.106	  75.072  1.235	 0.234	2.92780E-04  0.00211302
        137.895	  85.313  1.26	 0.224	2.92780E-04  0.00211302
        158.717	 102.020  1.301	 0.208	2.92780E-04  0.00211302
        172.369	 112.938  1.327	 0.200	2.92780E-04  0.00211302
        206.843	 140.580  1.395	 0.187	2.92780E-04  0.00211302
        241.317	 168.204  1.463	 0.175	2.92780E-04  0.00211302
        275.790	 195.829  1.5312 0.167	2.92780E-04  0.00211302
        310.264	 223.471  1.5991 0.159	2.92780E-04  0.00211302
        330.948	 240.053  1.6398 0.155	2.92780E-04  0.00211302
      END !end TEMP block
    END !endDATA
  END !end PVDO
  ENTHALPY QUADRATIC_TEMP ! h = c1*(t-t1) + 0.5*c2*(t-t2)
    TEMP_REF_VALUES   15.6d0  15.6d0    !t1, t2 [°C]
    TEMP_COEFFICIENTS  0.0     0.0d0    !c1, c2 [J/kg/°C]
  END !end Enthlapy
END !end EOS OIL

EOS GAS
  FORMULA_WEIGHT 16.04d0
  SURFACE_DENSITY 1.0995 kg/m^3
  PVDG
    DATA_UNITS ! Metric in the Eclipse sense
      PRESSURE  Bar
      FVF       m^3/m^3
      VISCOSITY cP
    END
    DATA
      TEMPERATURE 15.0
         1.0135  1.186979 0.0107
        34.4738  0.033261 0.0127
        68.9476  0.016001 0.0134
        82.7371  0.013160 0.0138
        103.421  0.010358 0.0145
        124.105  0.008535 0.0153
        137.895  0.007636 0.0159
        158.717  0.006597 0.0170
        172.369  0.006190 0.0177
        206.843  0.005531 0.0195
        241.317  0.005118 0.0214
        275.790  0.004840 0.0232
        310.264  0.004617 0.0250
        330.948  0.004509 0.0261
      END !end TEMP block
    END !endDATA
    !specify temperature dependency for viscosity
  END !end PVDG
END

EOS SOLVENT
  FORMULA_WEIGHT 23.76d0
  SURFACE_DENSITY 1 kg/m^3
  PVDS
    DATA_UNITS ! Metric in the Eclipse sense
      PRESSURE  Bar
      FVF       m^3/m^3
      VISCOSITY cP
    END
    DATA
      TEMPERATURE 15.0
          1.0135  1.25323  0.0110
         34.4738  0.03145  0.0120
         68.9476  0.01421  0.0130
         82.7371  0.01142  0.0140
        103.4214  0.00875  0.0160
        124.1057  0.00710  0.0180
        137.8952  0.00633  0.0190
        158.7173  0.00550  0.0220
        172.3690  0.00510  0.0230
        206.8428  0.00438  0.0270
        241.3166  0.00392  0.0310
        275.7904  0.00361  0.0340
        310.2642  0.00337  0.0370
        330.9485  0.00326  0.0380
      END !end TEMP block
    END !endDATA
    !specify temperature dependency for viscosity
  END !end PVDS
END

#=========================== wells ==================================

WELL_DATA injg

  RADIUS 0.0762 m
  CIJK_D 1 1 1 1  
  WELL_TYPE SOLVENT_INJECTOR
  BHPL     689.4757    Bar
  SHUT

  TIME 2 y
  OPEN
  WELL_TYPE WATer_INJECTOR
  TARG_WSV 1907 m^3/day

  TIME 3 y
  WELL_TYPE SOLVENT_INJECTOR
  TARG_SSV 339802 m^3/day

  TIME 4 y
  WELL_TYPE WATER_INJECTOR
  TARG_WSV 1907 m^3/day

  TIME 5 y
  WELL_TYPE SOLVENT_INJECTOR
  TARG_SSV 339802 m^3/day

  TIME 6 y
  WELL_TYPE WATER_INJECTOR
  TARG_WSV 1907 m^3/day

  TIME 7 y
  WELL_TYPE SOLVENT_INJECTOR
  TARG_SSV 339802 m^3/day

  TIME 8 y
  WELL_TYPE WATER_INJECTOR
  TARG_WSV 1907 m^3/day 

  TIME 9 y
  WELL_TYPE SOLVENT_INJECTOR
  TARG_SSV 339802 m^3/day

  TIME 10 y
  WELL_TYPE WATER_INJECTOR
  TARG_WSV 1907 m^3/day

  TIME 11 y
  WELL_TYPE SOLVENT_INJECTOR
  TARG_SSV 339802 m^3/day  

  TIME 12 y
  WELL_TYPE WATer_INJECTOR
  TARG_WSV 1907 m^3/day

  TIME 13 y
  WELL_TYPE SOLVENT_INJECTOR
  TARG_SSV 339802 m^3/day

  TIME 14 y
  WELL_TYPE WATER_INJECTOR
  TARG_WSV 1907 m^3/day

  TIME 15 y
  WELL_TYPE SOLVENT_INJECTOR
  TARG_SSV 339802 m^3/day

  TIME 16 y
  WELL_TYPE WATER_INJECTOR
  TARG_WSV 1907 m^3/day

  TIME 17 y
  WELL_TYPE SOLVENT_INJECTOR
  TARG_SSV 339802 m^3/day

  TIME 18 y
  WELL_TYPE WATER_INJECTOR
  TARG_WSV 1907 m^3/day

  TIME 19 y
  WELL_TYPE SOLVENT_INJECTOR
  TARG_SSV 339802 m^3/day

  TIME 20 y
  WELL_TYPE WATER_INJECTOR
  TARG_WSV 1907 m^3/day
END

WELL_DATA prod

  RADIUS 0.0762 m
  CIJK_D  7 7 3 3
  WELL_TYPE PRODUCER
  BHPL      68.9476 Bar
  TARG_OSV  1907 m^3/day

END

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

EQUILIBRATION

  PRESSURE 275.8 Bar
  DATUM_D  2560.32 m
  OGC_D    2500.0 m
  OWC_D    2600.0 m
  PCOG_OGC 0.0    Bar
  PCOW_OWC 0.0    Bar

  BUBBLE_POINT_TABLE
   D_UNITS m  
   PRESSURE_UNITS Bar
   PBVD
     1000.0 158.8
     3000.0 158.8
   /
  END

  RTEMP 15.0 C
/

END_SUBSURFACE