TITLE(2)="AIR WATER MODEL TEST 3: 3D, 2 Boundary conditions"

DESCRIPTION()=
"LENGTH (FT) : 15"
"WIDTH (FT) : 15"
"THICKNESS (FT) : 15"
"GRID BLOCKS : 50X11X5   1000 GRID ELEMENTS"
"DATE : 9/1/99"

AIR_WATER_MODEL 

TIMEEND = 3.

MAXNEWT = 200
CVTOL= 1.0E-07

$ if use GMRES
$LSOL_TOL = 1.E-5
LSOL_TOL = 1.E-8
LSOL_ITMAX = 1000
GMRES_PREC = 16
$ the parameter below is for GMRES : controls the number of Gauss-Seidel
$ block pressure solves
N_GS_STEP = 8

$ I/O OPTIONS

OUTLEVEL = 2  SPLINEOUT  GEOMOUT  PROCOUT
WELLOUTKEY = 3   WELLFILE = "air.wel"

$ FAULT BLOCK AND MESH DATA

DOWN() = 0 0 -1
NX(1) = 50  NY(1) = 11   NZ(1) = 5

DX()  = 3.36 3.06 2.78 2.53 2.30 2.09 1.90 1.73 1.57 1.43
1.30 1.18 1.07 0.97 0.89 0.81 0.73 0.67 0.61 0.55 10*0.50  
0.55 0.61 0.67 0.73 0.81 0.89 0.97 1.07 1.18 1.30
1.43 1.57 1.73 1.90 2.09 2.30 2.53 2.78 3.06 3.36

DY() = 3.60 2.78 2.30 1.90 1.50 0.50 1.50 1.90 2.30 2.78 3.60  

DZ() = 2.00 1.80 1.75 1.50 1.00

XYZ111(,1) = 3*0


ISTEP() = 1
JFIRST() =1 
JLAST() = 6
JSTEP() = 1
KSTEP() = 1

$ Boundary Conditions: geometry
$NBND_REG = 2
$NBND_SUBREG = 2

$             nreg nblk  x1 y1 z1      x2 y2 z2
$BND_VOL(,1) = 1     1,   0. 0. 0.      0. 15. 15.      $ x = 0, z-
$BND_VOL(,2) = 2     1,   15. 0. 0.     15. 15. 15.     $ x = 15, z+
$BND_VOL(,1) = 1     1,   0. 0. 0.      15. 15. 0.      $ z = 0, z-
$BND_VOL(,2) = 2     1,   0. 0. 15.     15. 15. 15.     $ z = 15, z+

$ boundary conditions: type and conditions

$BOUND_TYPE(1) = 1
$ABOUND1(1) Block $ PWAT
$  Interpolation Linear
$  Extrapolation Constant
$  Data  0.  12.5
$EndBlock

$ABOUND2(1) Block $ SWAT
$  Interpolation Linear
$  Extrapolation Constant
$  Data  0.  0.15
$EndBlock

$BOUND_TYPE(2) = 1
$ABOUND1(2) Block $ PWAT
$  Interpolation Linear
$  Extrapolation Constant
$  Data  0.  16.
$EndBlock

$ABOUND2(2) Block $ SWAT
$  Interpolation Linear
$  Extrapolation Constant
$  Data  0.  1.
$EndBlock

$ FLUID PROPERTIES

WATST = 62.42
STDENW = 62.42
AIRST = 0.288
STDENG = 0.288
ROCKCMP = 0.0
RGAS  = 10.732
MOLWT = 28.97     
TEMPERATURE = 60.

WCOMP =6.8E-5 
WATVIS = 1.0[cp]

$ INITIAL CONDITIONS

$DINIT = 0.0
$PWINIT = 14.269443
$SWINIT = .15

AREADIN
$WATERPRES1() = 14.7
WATERPRES1(,,1) = 17.536
WATERPRES1(,,2) = 16.713
WATERPRES1(,,3) = 15.945
WATERPRES1(,,4) = 15.241
WATERPRES1(,,5) = 14.7
WATERPSAT1() = 0.2

$ POROSITY

POROSITY1() = .4

$ PERMEABILITIES

XPERM1() = 18000.
YPERM1() = 18000.
ZPERM1() = 18000.

$ TECPLOT VISUALIZATION

$$$$$$ request to visualize table and lookup spline data : choose one
$VIS_TABOUTTYPE = 0         $ default : no vis output is produced
VIS_TABOUTTYPE = 1         $ suggested: output for all tables is produced 
$VIS_TABOUTTYPE = 2         $ ask for specific tables and stop simulation

$ TABLES

KASW(1) Block        $ AIR RELATIVE PERMEABILITY VS Sw - ROCK TYPE 1
  Interpolation Spline2
  Extrapolation Constant
  Constraint  1.0 At 0.2
  Constraint  0.001 At 0.95
  Derivative  0.0 At 0.95  
  Data 
  0.20 1.0000, 0.22 0.9538, 0.23 0.9089, 0.25 0.8653, 0.26 0.8229, 0.28 0.7818, 0.29 0.7420,
  0.31 0.7035, 0.32 0.6661, 0.34 0.6300, 0.35 0.5951, 0.37 0.5613, 0.38 0.5288, 0.40 0.4974,
  0.41 0.4671, 0.43 0.4380, 0.44 0.4100, 0.46 0.3831, 0.47 0.3572, 0.49 0.3325, 0.50 0.3088,
  0.52 0.2862, 0.53 0.2646, 0.55 0.2440, 0.56 0.2243, 0.58 0.2057, 0.59 0.1880, 0.61 0.1713,
  0.62 0.1555, 0.64 0.1406, 0.65 0.1266, 0.67 0.1135, 0.68 0.1012, 0.70 0.0897, 0.71 0.0791,
  0.73 0.0693, 0.74 0.0602, 0.76 0.0519, 0.77 0.0443, 0.79 0.0374, 0.80 0.0312, 0.82 0.0257,
  0.83 0.0208, 0.85 0.0165, 0.86 0.0128, 0.88 0.0097, 0.89 0.0070, 0.91 0.0049, 0.92 0.0032,
  0.94 0.0019, 0.95 0.0010
EndBlock

KWSW(1) Block        $ WATER RELATIVE PERMEABILITY VS Sw - ROCK TYPE 1
  Interpolation Spline2
  Extrapolation Constant
  Constraint  0.0 At .2
  Constraint  .5446 At .95
  Derivative  0.0 At .2
  Data  
  0.20 0.0000, 0.22 0.0015, 0.23 0.0044, 0.25 0.0080, 0.26 0.0123, 0.28 0.0172, 0.29 0.0226,
  0.31 0.0285, 0.32 0.0349, 0.34 0.0416, 0.35 0.0487, 0.37 0.0562, 0.38 0.0640, 0.40 0.0722,
  0.41 0.0807, 0.43 0.0895, 0.44 0.0986, 0.46 0.1080, 0.47 0.1176, 0.49 0.1276, 0.50 0.1378,
  0.52 0.1482, 0.53 0.1590, 0.55 0.1699, 0.56 0.1811, 0.58 0.1926, 0.59 0.2042, 0.61 0.2161,
  0.62 0.2282, 0.64 0.2406, 0.65 0.2531, 0.67 0.2659, 0.68 0.2789, 0.70 0.2920, 0.71 0.3054,
  0.73 0.3190, 0.74 0.3327, 0.76 0.3467, 0.77 0.3609, 0.79 0.3752, 0.80 0.3897, 0.82 0.4044,
  0.83 0.4193, 0.85 0.4344, 0.86 0.4496, 0.88 0.4650, 0.89 0.4806, 0.91 0.4964, 0.92 0.5123,
  0.94 0.5284, 0.95 0.5446
EndBlock

$ AIR-WATER CAPILLARY PRESSURE - ROCK TYPE 1
PCAW(1) Block        $ AIR-WATER CAPILLARY PRESSURE - ROCK TYPE 1
  Interpolation Linear
  Extrapolation Constant
  Data
  0.10 0.29814240, 0.125 0.28580647, 0.15 0.27364553, 0.175 0.26166215, 0.20 0.24985898,
  0.225 0.23823882, 0.25 0.22680461, 0.275 0.21555941, 0.30 0.20450646, 0.325 0.19364917,
  0.35 0.18299112, 0.375 0.17253611, 0.40 0.16228817, 0.425 0.15225156, 0.45 0.14243083,
  0.475 0.13283084, 0.50 0.12345679, 0.525 0.11431425, 0.55 0.10540926, 0.575 0.09674831,
  0.60 0.08833849, 0.625 0.08018754, 0.65 0.07230395, 0.675 0.06469713, 0.70 0.05737753,
  0.725 0.0503569, 0.75 0.04364857, 0.775 0.03726780, 0.80 0.03123237, 0.825 0.02556331,
  0.85 0.02028602, 0.875 0.0154321, 0.895 0.01188075, 0.90 0.01104231, 0.95 0.00390405,
  0.97 0.00181444, 1.0 0. 
EndBlock

ZFAC(1) Block    
  Interpolation Linear
  Extrapolation Constant
  Data
  147. 1.0 , 368. 1.001 , 735. 1.005 , 1470. 1.018
EndBlock

PHI(1) Block  
  Interpolation Linear
  Extrapolation Constant
  Data
  10.   0.4 , 4000.   0.4
EndBlock

AIRVIS(1) Block   
  Interpolation Linear
  Extrapolation Constant
  Data
  147. 0.015 , 735. .015 , 1470. .015 , 2940. .015
EndBlock


NUMWELL=10

WELLNAME(1) = "FIRST WELL"
KINDWELL(1) = 2
DIAMETER(1,1) = 0.1
WELLTOP(1 TO 3,1,1) = 33.78 12.33 8.05
WELLBOTTOM(1 TO 3,1,1) = 33.78 12.33 7.05
WELLPQ(1) Block
  Interpolation Linear
  Extrapolation Constant
  Data  0.  0.5343
EndBlock

WELLNAME(2) = "SECOND WELL"
KINDWELL(2) = 2
DIAMETER(1,2) = 0.1
WELLTOP(1 TO 3,1,2) = 34.28 12.33 8.05
WELLBOTTOM(1 TO 3,1,2) = 34.28 12.33 7.05
WELLPQ(2) Block
  Interpolation Linear
  Extrapolation Constant
  Data  0.  0.5343
EndBlock

WELLNAME(3) = "Third WELL"
KINDWELL(3) = 3
DIAMETER(1,3) = 0.2
WELLTOP(1 TO 3,1,3) = 1.68 4.99 8.05
WELLBOTTOM(1 TO 3,1,3) = 1.68 4.99 0.0
WELLPQ(3) Block
  Interpolation Linear
  Extrapolation Constant
  Data  0.  14.7
EndBlock

WELLNAME(4) = "Fourth WELL"
KINDWELL(4) = 3
DIAMETER(1,4) = 0.2
WELLTOP(1 TO 3,1,4) = 1.68 11.33 8.05
WELLBOTTOM(1 TO 3,1,4) = 1.68 11.33 0.0
WELLPQ(4) Block
  Interpolation Linear
  Extrapolation Constant
  Data  0.  14.7
EndBlock

WELLNAME(5) = "Fifth WELL"
KINDWELL(5) = 3
DIAMETER(1,5) = 0.2
WELLTOP(1 TO 3,1,5) = 1.68 15.03 8.05
WELLBOTTOM(1 TO 3,1,5) = 1.68 15.03 0.0
WELLPQ(5) Block
  Interpolation Linear
  Extrapolation Constant
  Data  0.  14.7
EndBlock

WELLNAME(6) = "Sixth WELL"
KINDWELL(6) = 3
DIAMETER(1,6) = 0.2
WELLTOP(1 TO 3,1,6) = 1.68 19.67 8.05
WELLBOTTOM(1 TO 3,1,6) = 1.68 19.67 0.0
WELLPQ(6) Block
  Interpolation Linear
  Extrapolation Constant
  Data  0.  14.7
EndBlock

WELLNAME(7) = "Seventh WELL"
KINDWELL(7) = 3
DIAMETER(1,7) = 0.2
WELLTOP(1 TO 3,1,7) = 66.38 4.99 8.05
WELLBOTTOM(1 TO 3,1,7) = 66.38 4.99 0.0
WELLPQ(7) Block
  Interpolation Linear
  Extrapolation Constant
  Data  0.  14.7
EndBlock

WELLNAME(8) = "Eighth WELL"
KINDWELL(8) = 3
DIAMETER(1,8) = 0.2
WELLTOP(1 TO 3,1,8) = 66.38 11.33 8.05
WELLBOTTOM(1 TO 3,1,8) = 66.38 11.33 0.0
WELLPQ(8) Block
  Interpolation Linear
  Extrapolation Constant
  Data  0.  14.7
EndBlock


WELLNAME(9) = "Ninth WELL"
KINDWELL(9) = 3
DIAMETER(1,9) = 0.2
WELLTOP(1 TO 3,1,9) = 66.38 15.03 8.05
WELLBOTTOM(1 TO 3,1,9) = 66.38 15.03 0.0
WELLPQ(9) Block
  Interpolation Linear
  Extrapolation Constant
  Data  0.  14.7
EndBlock

WELLNAME(10) = "Tenth WELL"
KINDWELL(10) = 3
DIAMETER(1,10) = 0.2
WELLTOP(1 TO 3,1,10) = 66.38 19.67 8.05
WELLBOTTOM(1 TO 3,1,10) = 66.38 19.67 0.0
WELLPQ(10) Block
  Interpolation Linear
  Extrapolation Constant
  Data  0.  14.7
EndBlock

EndInitial

$ TRANSIENT DATA INPUT BLOCKS

BeginTime    0.
$DELTIM = .0001  DTIMMUL = 1.05  TIMOUT = .1  
DELTIM = .01  DTIMMUL = 1.05   TIMOUT = .0  
DTIMOUT = .1  
DTIMMAX = 1.0   DTIMMIN =0.0001
VISOUT = .1    DVISOUT = .1
VIS_SCL = 4
VIS_SCL_NAMES(1) = "PWAT"
VIS_SCL_NAMES(2) = "SWAT"
VIS_SCL_NAMES(3) = "PAIR"
VIS_SCL_NAMES(4) = "SAIR"
VISFLAG = 2
VIS_FNAME = "AIRWATER"
EndTime

BeginTime  1.
  KINDWELL(1) = 0
  KINDWELL(2) = 0
EndTime

$BeginTime    2000.
$DTIMMAX = 10.
$EndTime





