! SPE10 benchmark, from:
! 'Tenth SPE Comparative Solution Project: A Comparison of Upscaling Techniques'
! M.A. Christie and M.J. Blunt, SPE 72469
! Note: this models the full-size version of the problem with no upscaling

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

SUBSURFACE

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

GRID
  TYPE grdecl spe10.grdecl
END

#=========================== solver options ===================================

NEWTON_SOLVER FLOW
  ATOL 1D-8
  RTOL 1D-8
  STOL 1D-50
  MAXIT 40
END

LINEAR_SOLVER FLOW
  
  CPR_OPTIONS
    CPR_HYPRE_OPTIONS
      BOOMERAMG_AGG_NL 1
    END
  END

  GMRES_RESTART 200
  RTOL 1e-5

END

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

# Actual final time is 2000 days

TIME
  FINAL_TIME 2000 d
  INITIAL_TIMESTEP_SIZE 1.0D-1  d 
  MINIMUM_TIMESTEP_SIZE 1.0D-10 d
  MAXIMUM_TIMESTEP_SIZE 100.0 d at 0. d 
END

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

OUTPUT
  MASS_BALANCE_FILE
   PERIODIC TIMESTEP 1  
  END
  ECLIPSE_FILE
    PERIOD_SUM TIMESTEP 1
    PERIOD_RST TIMESTEP 10
  END
  LINEREPT
END

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

MATERIAL_PROPERTY formation
  ID 1
  ROCK_COMPRESSIBILITY    1.450377E-5 1/Bar  # Use 1.0E-6 1/psi from website
  ROCK_REFERENCE_PRESSURE 344.73786   Bar    # Reference to 5000 psi
  CHARACTERISTIC_CURVES ch1
/

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

CHARACTERISTIC_CURVES ch1
  CAP_PRESSURE_FUNCTION_OW VAN_GENUCHTEN
   WATER_CONNATE_SATURATION 0.2
   WATER_RESIDUAL_SATURATION 0.20  
   M 0.5d0
   ALPHA 5.792e-05
   MAX_CAPILLARY_PRESSURE 0.d0
  /
  PERMEABILITY_FUNCTION_WAT MOD_BROOKS_COREY
   M 2.0d0
   WATER_CONNATE_SATURATION 0.2
   WATER_RESIDUAL_SATURATION 0.20
   MAX_REL_PERM 1.0d0 !default is 1 anyway
  END
  PERMEABILITY_FUNCTION_OW MOD_BROOKS_COREY
   M 2.0d0
   OIL_RESIDUAL_SATURATION 0.2d0
   MAX_REL_PERM 1.0d0 !default is 1 anyway
  END  
/
#=========================== EOSs =============================================

# Water data is Bw=1.01, Cw=3.0E-6, viscw=0.3, surface density 64 lbs/cuft
# Converted to reference 5000 psi, Bw=1.01, and density=64/1.01 lbs/cuft

EOS WATER
  SURFACE_DENSITY 1025.18
  DENSITY EXPONENTIAL 1015.0313 34473786 4.351132E-10
  VISCOSITY CONSTANT  0.0003   ! 0.3 cp
END

# Oil data is specified by a table, and a surface density of 53 lbs/cuft
# Converted to compressibility form at reference 5000 psi

EOS OIL
  SURFACE_DENSITY 848.979 kg/m^3
  FORMULA_WEIGHT 1000.0d0
  DENSITY INVERSE_LINEAR ! den = den0 / [1- C(p-p1)] [1 + E(t-t1)]
    REFERENCE_VALUE 837.1193678         !den0 [kg/m3]  
    PRES_REF_VALUE 34473786             !p1 [Pa]
    TEMP_REF_VALUE 15.0d0               !t1 [C]
    COMPRESS_COEFF 1.986274E-10         !C [1/Pa]
    THERMAL_EXPANSION_COEFF 0.0         !E [1/C]
  END
  ENTHALPY QUADRATIC_TEMP ! h = c1*(t-t1) + 0.5*c2*(t-t2) # Not used
    TEMP_REF_VALUES   15.6d0  15.6d0    !t1, t2 [C]
    TEMP_COEFFICIENTS  0.0     0.0d0    !c1, c2 [J/kg/C]
  END
  VISCOSITY CONSTANT 0.002994  ! # Use a constant viscosity of visco(4000)=2.994 cp
END

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

WELL_DATA INJW
  RADIUS    0.1524 m
  WELL_TYPE WATER_INJECTOR
  BHPL      1000    Bar
  TARG_WSV  794.9365703125 m^3/day
  CIJK_D 30 110 1 85  
END

WELL_DATA prod11
  RADIUS         0.1524 m
  THETA_FRACTION 0.25
  SKIN_FACTOR    0.693147 
  WELL_TYPE      PRODUCER
  BHPL           275.7903 Bar
  TARG_OSV       1000     m^3/day
  CIJK_D           1 1 1 85
END

WELL_DATA prod21
  RADIUS         0.1524 m
  THETA_FRACTION 0.25
  SKIN_FACTOR    0.693147 
  WELL_TYPE      PRODUCER
  BHPL           275.7903 Bar
  TARG_OSV       1000     m^3/day
  CIJK_D           60 1 1 85
END

WELL_DATA prod12
  RADIUS         0.1524 m
  THETA_FRACTION 0.25
  SKIN_FACTOR    0.693147 
  WELL_TYPE      PRODUCER
  BHPL           275.7903 Bar
  TARG_OSV       1000     m^3/day
  CIJK_D           1 220 1 85
END

WELL_DATA prod22
  RADIUS         0.1524 m
  THETA_FRACTION 0.25
  SKIN_FACTOR    0.693147 
  WELL_TYPE      PRODUCER
  BHPL           275.7903 Bar
  TARG_OSV       1000     m^3/day
  CIJK_D           60 220 1 85
END

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

EQUILIBRATION

  PRESSURE 413.6854 Bar
  DATUM_D 3657.6 m
  OWC_D   4000.0d0 m
  PCOW_OWC 0.0    Bar
  RTEMP 15.0d0 C
/

END_SUBSURFACE