UTCHEM Users Guide PDF

Volume I:
User's Guide
for
UTCHEM-9.0
A Three-Dimensional
Chemical Flood Simulator
Prepared by
Reservoir Engineering Research Program

Center for Petroleum and Geosystems Engineering
The University of Texas at Austin
Austin, Texas 78712
July, 2000
UTCHEM User's Guide
TABLE OF CONTENTS
Page
1.0 INTRODUCTION 1
2.0 OPERATION OF THE SIMULATOR 3
2.1 Input and Output Files 3
2.2 Description of HEAD file 6
2.3 Compilation and Execution on Workstations 7
2.4 Compilation and Execution on a WINDOWS-based PC 9
3.0 INPUT DATA DESCRIPTION 11
3.1 Title and Reservoir Description Data 13
3.2 Output Option Data 23
3.3 Reservoir Properties 27
3.4 General Physical Property Data 43
3.5 Physical Property Data for Geochemical Options 89
3.6 Data for Biodegradation Option 97
3.7 Recurrent Injection/Production Data Set 107
APPENDICES
A OUTPUT FILES 125
A.1 Data Written to Profile Data File 125
A.2 Restart Run Procedure 126
A.3 Data Written to Stored Restart Run Data File 127
A.4 Data Written to Well History Plotting Data File(s) 129
A.5 Data Written to Overall History Plotting Data File 130
A.6 Data Written to Tracer Concentration Observation Point Data
File(s) 133
A.7 Output Files Created When Running GMS 132
B GEOCHEMISTRY OPTION (IREACT>1) 133
C MAIN PROGRAM FLOW OUTLINE 137
D PHASES AND SPECIES IN UTCHEM 141
E TIME-STEP SELECTION 143
E.1 Method of Relative Changes for First Three Components
(IMES=2) 143
E.2 Method of Relative Changes for All Components (IMES=3) 145
E.3 Method of Relative Changes Using Dimensionless
Concentration for All Components (IMES=4) 146
F CHANGES TO USER'S GUIDE 147
F.1 New Variables 147
F.2 Modified Variables 147
F.3 Other Changes 148
INPUT VARIABLE INDEX 149

UTCHEM User's Guide
INTRODUCTION
1.0 INTRODUCTION
UTCHEM is a three-dimensional chemical flooding simulator. The solution
scheme is analogous to IMPES, where pressure is solved for implicitly, but
concentrations rather than saturations are then solved for explicitly. Phase
saturations and concentrations are then solved in a flash routine. An energy balance
equation is solved explicitly for reservoir temperature. The energy balance equation
includes heat flow between the reservoir and the over- and under-burden rocks.
The major physical phenomena modeled in the simulator are:
• dispersion • organic biodegradation capability

• diffusion • multiple organic species
• dilution effects • equilibrium and nonequilibrium
• adsorption for oil, surfactant and organic dissolution in aqueous
polymer phase
• interfacial tension • dual porosity option for simple
• relative permeability phase tracer flow
• capillary pressure • polymer properties: shear thinning
• hysteresis in relative permeability viscosity, inaccessible pore
and capillary pressure volume, permeability reduction,
• capillary trapping adsorption
• cation exchange • gel properties: viscosity,
• phase density permeability reduction, adsorption
• compositional phase viscosity • tracer properties: partitioning,
• phase behavior (pseudoquaternary) adsorption, radioactive decay,
• aqueous reactions reaction (ester hydrolization),
• partitioning of chemical species dead-end pore (capacitance)
between oil and water • temperature dependent properties:
• dissolution/precipitation viscosity, tracer reaction, gel
• cation exchange reactions involving reactions, Surfactant phase
more than two cations behavior
• in-situ generation of surfactant from • gas mobility reduction due to foam
acidic crude oil • mixed-wet oil/water capillary
• pH dependent surfactant adsorption pressure and relative permeability
See Section 2 of the UTCHEM Technical Documentation for the general

formulation of the simulator.
1
UTCHEM User's Guide
INTRODUCTION
User notes:
2
UTCHEM User's Guide
OPERATION OF THE SIMULATOR
2.0 OPERATION OF THE SIMULATOR

The UTCHEM simulator is run on a CRAY SV1 at the Texas Advanced
Computing Center affiliated with The University of Texas at Austin (UNICOS
operating system), a DEC Alpha system (DEC/Compaq XP1000) at the Center for
Petroleum and Geosystems Engineering (OSF/1 operating system), and several
WINDOWS-based PC workstations (Windows operating system, running at 450
MHz with 256MB RAM; depending on the size of the problem, the code will also
run on slower PCs with less RAM).
2.1 Input and Output Files

UTCHEM requires two input files (HEAD and INPUT) for non-restart
runs. For restart runs, a binary restart file (INPUT2) is required in addition to the
original input data files used for the previous run. The HEAD input file gives
UTCHEM information regarding the name and size of the problem to be run; it is
described in Section 2.2. A detailed description of the data in the second input file
is given in Section 3.0 and the data in the restart data file is documented in
Appendix A.3. A number of UTCHEM input example files demonstrating a variety
of petroleum oil-field and groundwater applications are available to UTCHEM
users. The oil-field applications include water, single-well tracer, interwell tracer,
polymer, profile control using gel, surfactant/polymer, and high pH
alkaline/surfactant/polymer flooding. The groundwater applications include
contaminant infiltration, water flushing, partitioning interwell tracer, surfactant
enhanced aquifer remediation, and bioremediation.
We provide all users with two sample input files for testing purposes (see
the distribution package for copies of the files). The first sample input file (EX01)
is for a 3D surfactant/polymer flood. In that file, the surfactant properties are for
petroleum sulfonate and the polymer properties are based on the typical data for
xanthan gum. CPU usage for the EX01 example run is about 55 seconds on the
DEC/Compaq XP1000 and about 286 seconds on the 450MHz Windows PC. The
second sample input file (EX21) is for a 2-D contamination event in the saturated
zone of an aquifer. CPU usage for the EX21 example run is about 14 seconds on
the DEC/Compaq XP1000 and about 90 seconds on the 450MHz Windows PC. If
3
UTCHEM User's Guide
you would like to receive additional sample input files, please contact Joanna L.
Castillo ([email protected] or 512-471-3229) for details.
The input files are summarized in the following table.
Input File Name Contents

HEAD Run number and array dimensioning information
INPUT Input data; described in detail in Section 3.0
INPUT2 Restart run data (binary input file RESTAR created by an
earlier run)
The number of output files generated by UTCHEM varies depending upon
several control flags set by the user in the input file. See Section 3.0 for more
information on how to set those flags. Output file names are created by appending
a suffix to the run number specified in the HEAD input file. For example, if the run
number is TEST01, the echo print output file generated by UTCHEM would be
TEST01.ECHO. All output files are ASCII files which can be edited with a
standard text editor unless otherwise specified. The output files are summarized in
the following table.
Output file Contents

extension
ALKP Alkaline option related profile plotting data, created if
IREACT>1 and IPALK=1
BIOD Biodegradation option plotting data, created if IPBIO=1 and
IBIO=1
CAPP Capacitance property profile plotting data, created if IPCAP=1
and ICAP=1
COMP_AQ Aqueous phase composition plotting data, created if ICKL=1
COMP_ME Microemulsion phase composition plotting data, created if
surfactant is present and ICKL=1
COMP_OIL Oleic phase composition plotting data, created if ICKL=1
CONCP Component concentration profile plotting data, created if
IPCTOT=1
SALT Effective salinity plotting data, created if surfactant is present and
ICSE=1
DISS Mass transfer option plotting data, created if INONEQ=1 and
IMASS=2 and ICOR>0
ECHO Echo print of the input file information
FOAM Foam option plotting data, created if IGAS=2 and IFOAMP=2
FRAC Dual porosity plotting data, created if IPCAP=1 and ICAP=2
GFILEP Gel property profile plotting data, created if IREACT=1 or 4 and
IPGEL=1
4
UTCHEM User's Guide
HIST01 Well history plotting data for well #1; described in Appendix A
HIST02 Well history plotting data for well #2
… …
HISTxx Well history plotting data for last well
HYST Hysteresis option plotting data, created if (IHYST=1 or
IPERM=3) and IHYSTP=1
MESH Number of gridblocks and distances to center of gridblocks in
each direction
OVERAL History of overall properties; described in Appendix A
PERM Permeability and porosity data
PRESP Phase pressure profile plotting data, created if IPPRES=1
PROF Formatted profile data; described in Appendix A
RESTAR Stored restart run data; described in Appendix A (binary file)
RPERM Phase relative permeability plotting data, creating if IPER=1
SATP Phase saturation profile plotting data, created if IPSAT=1
TEMPP Temperature profile, created if IENG=1 and IPTEMP=1
TRAC01 Aqueous (or gas) phase tracer concentration for the 1st tracer at
observation points, created if IPOBS>0; described in Appendix
A
TRAC02 Aqueous (or gas) phase tracer concentration for 2nd tracer at
observation points
… …
TRACxx Aqueous (or gas) phase tracer concentration for the last tracer at
observation points
TRAP Interfacial tension and phase trapping plotting data, created if
ICNM=1 and ITRAP>0
TTABLE Table of time steps and Courant numbers
VISC Phase viscosity plotting data, created if IVIS=1
WARN Warning messages
5
UTCHEM User's Guide
2.2 Description of HEAD file
The variables listed in the table below give UTCHEM information regarding
the name and size of the problem to be run. All values should be equal to those
specified in the main input file (INPUT). The values can be larger than those
specified in the main input file but it is recommended they not be in order to
maintain as small an executable file as possible. Note that the values must not be
smaller than those in the main input file or the program will not execute properly.
Example HEAD file:
utex21
NX NY NZ N NWELL
49 1 24 8 3
NTW NTA
0 0
NO NPHASE
0 4
NSUB MSUB
0 0
The first line of the HEAD file contains the run number which is a 6-character string
used as a prefix for output file names generated by UTCHEM (as described in
Section 2.1). We strongly recommend this be the same as RUNNO in the main
input file so as to avoid confusion. The following table contains definitions of the
variables used in the HEAD file.
Variable Definition
NX Number of gridblocks in X-direction
NY Number of gridblocks in Y-direction
NZ Number of gridblocks in Z-direction
N Number of components (cannot be less than 8)
NWELL Number of wells
NTW Number of water tracers
NTA Number of air tracers
NO Number of organic components
NPHASE Number of phases (set to 3 when there is no gas phase or 4 if gas is
present)
NSUB Number of subgrids in lateral direction (check for ICAP=2)
MSUB Number of subgrids in vertical direction (check for ICAP=2)
6
UTCHEM User's Guide
2.3 Compilation and Execution on Workstations
The UTCHEM distribution package for workstations contains the
FORTRAN-90 source files (all ending with a ".f" file extension), two files used
for compiling the code (Commodule and Makev9), input files for the ex01
sample run (ex01.data and ex01.head), and input files for the ex21 sample
run (ex21.data and ex21.head). Make sure you place all files in the same
directory on your workstation. Then, follow these steps in order to compile and
execute the code:
1. Issue the command:
Commodule
This will compile the modules before building the executable file.
2. To build an executable file called utchem9.exe, issue the command:
make -f Makev9 FC=f90
Note that the example makefile assumes the use of the f90 FORTRAN
compiler.
3. Run UTCHEM in the background using a command file (in the example
below, the file is called work.job). Note that the work.job file needs
to be executable. Use any text editor to create your job file and then issue
the command:
chmod +x work.job
at the UNIX prompt to make sure the work.job file is executable. Then,
submit the job file by issuing the following command:
work.job &
at the UNIX prompt.
Example work.job file:
rm -r ex01.dir
mkdir ex01.dir
cd EX01.dir
ln -s ../ex01.head HEAD
ln -s ../ex01.data INPUT
time ../utchem9.exe
gzip *
7
UTCHEM User's Guide
The work.job file needs to be modified to reflect the directory structure you
create for running UTCHEM jobs. See the table on the following page for a
description of the contents of the work.job file.
Command line Description

rm -r ex01.dir Removes the ex01.dir directory if it already
exists. Make sure you've copied files from
previous runs (if you want to save them) to
another location before executing the
w o r k . j o b file because all files in the
ex01.dir directory will be deleted when the
directory is deleted.
mkdir ex01.dir Create a subdirectory in which to place the new
simulation results.
cd ex01.dir Make e x 0 1 . d i r the current working
directory.
ln -s ../ex01.head HEAD Create symbolic links to the ex01.head and
ln -s ../ex01.data INPUT ex01.data input files (which are located up
one level in the directory structure). When the
program looks for the files HEAD and INPUT,
it will automatically be pointed to the
e x 0 1 . h e a d and e x 0 1 . d a t a files,
respectively.
time ../utchem9.exe Run the program (which is located up one level
in the directory structure).
gzip * Add this line to the work.job file to compress
all files in order to save space (if necessary).
8
UTCHEM User's Guide
2.4 Compilation and Execution on a WINDOWS-based PC
The UTCHEM distribution package for PCs contains the FORTRAN-90
source files (all ending with a ".f" file extension), input files for the ex01 sample
run (ex01.data and ex01.head), and input files for the ex21 sample run
(ex21.data and ex21.head). At The University of Texas at Austin, we
compile UTCHEM on PCs using the Digital FOFTRAN compiler. UTCHEM will
probably compile successfully using other compilers with little or no modification if
those compilers have been updated within the last year and contain FORTRAN-90
features. Copy the source files (all files ending with the ".f" file extension in the
distribution package) to a directory on your hard disk and then follow the steps
below to compile UTCHEM using the Digital FORTRAN compiler:
1. Launch Visual Studio.
2. Select "New" from the File menu. In the dialog box that appears,
click on the "Projects" tab and select Win32 Console Application,
provide a project name ("utchem9" for this example), and select the
location on your hard disk in which you'd like the executable code
to be placed. Click on OK.
3. Select "Add to Project…" and then "File…" from the Project menu.
Navigate to the location of the source files. Select "All files (*.*)"
so all the files can be seen. Control-click on the three module source
files (module1.f, module11.f, and module12.f). Click on
OK.
4. Select "Build utchem9.exe" from the Build menu. This will

generate a couple of error messages regarding the Library module
file and error executing df.exe. Disregard these message.
5. Select "Add to Project…" and then "File…" from the Project menu.
Navigate to the location of the source files. Select "All files (*.*)"
so all the files can be seen. Select all the source files. Click on OK.
6. Select "Build utchem9.exe" from the Build menu. You will get a
warning that the module files have already been compiled. Click on
OK.
7. Wait patiently while the code compiles. Upon completion, exit

Visual Studio. An executable version of the code called
utchem9.exe will be placed in the Debug folder within the
utchem9 project folder.
9
UTCHEM User's Guide
8. Copy the input and head files to the directory containing

utchem9.exe and name them INPUT and HEAD. Make sure the
files have no file extensions (like ".txt" or ".dat").
9. Double click on utchem9.exe. A screen will appear and will

remain visible until execution of the job has completed.
Periodically, the simulation time, the number of time steps, time step
in days, and Courant number will be written to the screen. Your
output files will appear in the same directory from which you
launched utchem9.exe.
10
UTCHEM User's Guide
INPUT DATA DESCRIPTION
3.0 INPUT DATA DESCRIPTION
The main UTCHEM input file (called INPUT) consists of comment lines
and data lines. All comment lines are ignored by the UTCHEM simulator. It is
important to note, however, that the number of comment lines between data lines is
fixed. The first twenty-two lines of the input file are reserved for comment lines
used to briefly describe the input file. Each data line is preceded by three comment
lines (except for the input described in Section 3.5). The input file is basically
divided into seven sections and each of those input sections (except Section 3.5) is
preceded by an additional seven comment lines. The user should update the
comment lines as the input file is modified in order to make using the simulator
easier.
All non-character data is free-formatted. This means that for each read
statement, it is only necessary to leave a blank space between data elements. Note
that the first data element for a given read statement must be on a new line in the
input file. Subsequent data elements for that read statement can span as many lines
as necessary. Implicit type matching is used; that is, all REAL variables begin with
the letters A-H or O-Z and all integer variables begin with the letters I-N.
The following is a list of variables as they are read by UTCHEM. The
variable names appear in all-caps on a single line in the order they are read by the
program (variables that are new to the latest version of UTCHEM are printed in
italicized boldface as well). Every list of variables is followed by a description
of each variable and corresponding units or possible values if applicable. All of the
variables listed in the input description will be read by the program unless otherwise
noted; therefore, a dummy value will be read by the program for variables not
pertinent to the problem being run.
11
UTCHEM User's Guide
User notes:
12
UTCHEM User's Guide
3.1. Title and Reservoir Description Data
The first input section consists of the title and reservoir description data.
Please remember that there are 22 comment lines at the beginning of this section and
that each data line is preceded by three comment lines.
3.1.1 RUNNO
RUNNO - Run number.

Note: The run number can consist of any combination of six non-blank
alphanumeric characters. This information will be printed as the
first line of every output file and used as a prefix for file names of
output files. We strongly recommend that this be identical to the
first line of the HEAD file so as to avoid confusion. If you are
using the UTCHEM-GUI program, this will be done automatically.
3.1.2 TITLE
TITLE - Title and run description.

Note: The title can consist of any combination of alphanumeric characters
spanning three lines in the input file (not to exceed 80 characters per
line). The title must span three lines and any of those lines can be
blank.
3.1.3 IMODE, IMES, IDISPC, ICWM, ICAP, IREACT, IBIO, ICOORD, ITREAC,
ITC, IGAS, IENG
IMODE - Flag indicating if the problem to be run is a first run or a restart

simulation.
Possible Values:
1 - First simulation run
2 - Restart simulation
Note: See Appendix A for more details on how to run restart simulations.
IMES - Flag indicating if a constant or automatic time-step is to be used.

Possible Values:
1 - Constant time-step size is used
2 - Automatic time-step size selector based on method of relative
changes for the first three components is used
3 - Automatic time-step size selector based on method of relative
changes for all the components is used
4 - Automatic time-step size selector based on changes in
dimensionless concentration for all the components is used
Note: The automatic time-step selector is recommended. See input lines
3.7.9 through 3.7.13 and Appendix E for more details on the above
options.
IDISPC - Flag indicating which type of numerical dispersion control is used.

Possible Values:
0 - Single point upstream method is used
2 - Two point upstream method is used
3 - Improved total variation diminishing third order method is used
13
UTCHEM User's Guide
Note: These methods are applied to both concentration and relative

permeability.
ICWM - Flag indicating if the concentration well model is used or not.

Possible Values:
0 - Concentration well model is not used
1 - Concentration well model is used
Note: The concentration well model (ICWM=1) can only be used with
vertical wells (IDIR(M)=3).
ICAP - Flag indicating if the capacitance model is used or not.

Possible Values:
0 - Capacitance model or dual porosity option is not used
1 - Capacitance model is used
2 - Dual porosity option is used for single phase tracer flow
Note: The dual porosity option (ICAP=2) is available only if IMODE=1,
IUNIT=0, and ICOORD=1.
IREACT - Flag indicating if gel reactions or alkaline options are used or not.
Possible Values:
0 - Gel reactions are not used
1 - Gel reactions are used
2 - Geochemistry option with no acidic crude is used
3 - Geochemistry option with acidic crude is used
4 - IREACT=2 and gel reactions are used
IBIO - Flag indicating whether or not biodegradation reactions occur

Possible Values:
0 - No biodegradation reactions
1 - Biodegradation reactions occur
ICOORD - Flag indicating which coordinate system is used.

Possible Values:
1 - Cartesian coordinate system is used
2 - Radial coordinate system is used
3 - Cartesian coordinate system with variable-width gridblock is
used (2-D cross section only)
4 - Curvilinear grid definition of the X-Z cross section is used (2-D
or 3-D)
Note: For ICOORD=4, the 3-D grid consists of the 2-D cross sectional
grid repeated at specified intervals (uniform or non-uniform) in the
Y direction, according to the definition of DY1. The curvilinear grid
option is not available with the temperature equation option.
ITREAC - Flag indicating if a tracer reaction is used or not.

Possible Values:
0 - Tracer reactions are not used
1 - Tracer reactions are used
ITC - Flag indicating if second-order time approximation is used or not.

Possible Values:
0 - Second-order time approximation is not used
1 - Second-order time approximation is used
14
UTCHEM User's Guide
Note: We recommend that second-order time approximation (ITC=1) only

be used with higher-order dispersion methods (IDISPC>1).
IGAS - Flag indicating if gas phase is considered or not.

Possible Values:
0 - Gas is not present
1 - Gas is present
2 - Gas is present and foam option is used
IENG - Flag indicating if temperature variation is considered or not.

Possible Values:
0 - Isothermal simulation
1 - Temperature equation is solved
3.1.4 NX, NY, NZ, IDXYZ, IUNIT
NX - Number of gridblocks along X-axis (ICOORD=1, 3, or 4) or number of

gridblocks in radial direction (ICOORD=2).
Note: This value should be equal to or smaller than the NNX parameter in
UTCHEM.
NY - Number of gridblocks along Y-axis.

Note: This value should be equal to or smaller than the NNY parameter in
UTCHEM. It should be set equal to 1 if the user is running a 1-D
problem or a 2-D cross sectional problem. If ICOORD=2, this
value is automatically set equal to 1.
NZ - Number of gridblocks along Z-axis.

Note: This value should be equal to or smaller than the NNZ parameter in
UTCHEM. It should be set equal to 1 if the user is running a 1-D
problem or a 2-D areal problem.
IDXYZ - Flag indicating constant or variable grid size.

Possible Values:
0 - Constant grid size
1 - Variable grid size on a regional basis
2 - Variable grid size
Note: IDXYZ must be set equal to 2 if ICOORD=3.
IUNIT - Flag indicating English or Metric units.

Possible Values:
0 - English unit
1 - Metric unit
Note: These NX, NY, and NZ input values must be less than or equal to the NX,
NY, and NZ values listed in the HEAD input file.
15
UTCHEM User's Guide
Simulation Gridblock Sizes (input lines 3.1.5-3.1.23)
3.1.5 XCORD(I), ZCORD(I), for I=1, (NX+1)×(NZ+1) (This line is read only if
ICOORD=4)
Note: One line is required for each coordinate pair.
XCORD - Gridblock coordinate of Ith corner point in X-direction.

Units: feet (IUNIT=0) or m (IUNIT=1)
ZCORD - Gridblock coordinate of Ith corner point in Z-direction.

Note: The coordinates of the corners (or vertices) of the 2-D X-Z cross section
gridblocks are input in pairs as follows:
XCORD(1), ZCORD(1)
. .
. .
. .
XCORD(nodes), ZCORD(nodes)
where nodes = (NX+1) × (NZ+1) and is the total number of corner points
defining the X-Z cross section and Z is positive downward. The following
figure illustrates the input order for an example X-Z cross section grid:
Top (surface) of reservoir
XCORD(1), ZCORD(1) XCORD(2), ZCORD(2)
1 2 3
4
5
7 6
8
X
Z 9
XCORD(9), ZCORD(9)
The number of gridblocks is equal to NX × NZ and the number of

coordinate pairs (or nodes) is equal to (NX+1) × (NZ+1).
Cautionary warning: The X-Z cross section grid should be constructed by
the user such that the curvilinear coordinate system is at least quasi-
orthogonal. Departure from orthogonality will lead to numerical
errors in the solution.
16
UTCHEM User's Guide
Refer to the following table to help determine which input lines should be used to specify
the gridblock size input values for different options:
ICOORD IDXYZ Input line(s)
1 (Cartesian coordinate system) 0 (constant grid size) 3.1.6

1 (regional variable grid) 3.1.9–3.1.11
2 (variable grid) 3.1.16, 3.1.17,
3.1.19
2 (Radial coordinate system) 0 (constant grid size) 3.1.7

1 (regional variable grid) 3.1.12–3.1.14
2 (variable grid) 3.1.20–3.1.22
3 (Cartesian coordinate system 0 (constant grid size) (not applicable)

with variable-width gridblock) 1 (regional variable grid) (not applicable)
2 (variable grid) 3.1.16, 3.1.18,
3.1.19
4 (Curvilinear grid) 0 (constant grid size) 3.1.8

1 (regional variable grid) 3.1.15
2 (variable grid) 3.1.17
3.1.6 DX1, DY1, DZ1 (This line is read only if IDXYZ=0 and ICOORD=1)
DX1 - Gridblock size in X direction.

DY1 - Gridblock size in Y direction.

DZ1 - Gridblock size in Z direction.

3.1.7 R(1), DX1, DZ1 (This line is read only if IDXYZ=0 and ICOORD=2)
R(1) - Wellbore radius.

DX1 - Distance between nodes in radial direction.


3.1.8 DY1 (This line is read only if IDXYZ=0 and ICOORD=4)
17
UTCHEM User's Guide
3.1.9 II1, II2, DX1 (This line is read only if IDXYZ=1 and ICOORD=1)
II1 - First index for gridblocks with same size in X direction.
II2 - Last index for gridblocks with same size in X direction.
DX1 - Gridblock size in X direction

Note: This line is repeated until sizes for each of the NX gridblocks in the X
direction have been specified. The first line in the set must have II1=1 and
the last line must have II2=NX.
Example: If NX=11 and the first three gridblocks in the X direction are 3 feet in
size, the fourth through ninth gridblocks in the X direction are 2 feet in size,
and the last two gridblocks in the X direction are 2.5 feet in size, this line
would need to be repeated three times to fully describe the X direction
gridblocks as follows:
1 3 3.0
4 9 2.0
10 11 2.5
3.1.10 JJ1, JJ2, DY1 (This line is read only if IDXYZ=1 and ICOORD=1)
JJ1 - First index for gridblocks with same size in Y direction.
JJ2 - Last index for gridblocks with same size in Y direction.

Note: This line is repeated until sizes for each of the NY gridblocks in the Y
direction have been specified. The first line in the set must have JJ1=1 and
the last line must have JJ2=NY. See the example for input line 3.1.9.
3.1.11 KK1, KK2, DZ1 (This line is read only if IDXYZ=1 and ICOORD=1)
KK1 - First index for gridblocks with same size in Z direction.
KK2 - Last index for gridblocks with same size in Z direction.

Note: This line is repeated until sizes for each of the NZ gridblocks in the Z
direction have been specified. The first line in the set must have KK1=1
and the last line must have KK2=NZ. See the example for input line
3.1.9.
18
UTCHEM User's Guide
3.1.12 R(1) (This line is read only if IDXYZ=1 and ICOORD=2)

3.1.13 II1, II2, DX1 (This line is read only if IDXYZ=1 and ICOORD=2)
II1 - First index for radial node distances of the same size.
II2 - Last index for radial node distances of the same size.
DX1 - Distance between nodes in radial direction.

Note: This line is repeated until the NX-1 distances between the NX nodes in the
radial direction have been specified. The first line in the set must have
II1=1 and the last line must have II2=NX-1.
Example: If NX=35 and the first ten gridblocks in the X direction are 1 foot in
size and the rest are 2 feet in size, this line would need to be repeated twice
to fully describe the radial direction nodes as follows:
1 10 1.0
11 34 2.0
3.1.14 KK1, KK2, DZ1 (This line is read only if IDXYZ=1 and ICOORD=2)
KK1 - First index for gridblocks with same size in Z direction.
KK2 - Last index for gridblocks with same size in Z direction.

Note: This line is repeated until sizes for each of the NZ gridblocks in the Z
direction have been specified. The first line in the set must have KK1=1
and the last line must have KK2=NZ. See the example for input line
3.1.9.
3.1.15 JJ1, JJ2, DY1 (This line is read only if IDXYZ=1 and ICOORD=4)
JJ1 - First index for gridblocks with same size in Y direction.
JJ2 - Last index for gridblocks with same size in Y direction.

Note: This line is repeated until sizes for each of the NY gridblocks in the Y
direction have been specified. The first line in the set must have JJ1=1 and
the last line must have JJ2=NY. See the example for input line 3.1.9.
19
UTCHEM User's Guide
3.1.16 DX(I), for I=1, NX (This line is read only if IDXYZ=2 and ICOORD=1 or 3)
DX(I) - Grid size of Ith block in X direction.

3.1.17 DY(J), for J=1, NY (This line is read only if IDXYZ=2 and ICOORD=1 or 4)
DY(J) - Grid size of Jth block in Y direction.

3.1.18 DY(I), for I=1, NX (This line is read only if IDXYZ=2 and ICOORD=3)
DY(I) - Thickness of Ith block.

3.1.19 DZ(K), for K=1, NZ (This line is read only if IDXYZ=2 and ICOORD=1 or 3)
DZ(K) - Grid size of Kth block in Z direction.

3.1.20 R(1) (This line is read only if IDXYZ=2 and ICOORD=2)

3.1.21 DX(I), for I=1, NX-1 (This line is read only if IDXYZ=2 and ICOORD=2)
DX(I) - Distance between the Ith node and the I+1th node in the radial direction.
3.1.22 DZ(K), for K=1, NZ (This line is read only if IDXYZ=2 and ICOORD=2)
DZ(K) - Grid size of Kth block in Z direction.

3.1.23 N, NO, NTW, NTA, NGC, NG, NOTH
N - Total number of components in the run (including tracers and reactive

components).
Value must be set equal to: N=8+NO+NTW+NTA+NGC+NG+NOTH
NO - Total number of NAPL phase organic components in the run.

Note: If IBIO=1, set NO=0 if no NAPL phase is present and all
biodegrading species are present only in the aqueous phase;
otherwise, set NO to the number of organic species.
NTW - Number of water/oil tracers.
NTA - Number of oil/gas tracers.
NGC - Number of components for geochemistry option.
NG - Number of gel components.
20
UTCHEM User's Guide
NOTH - Total number of other chemical and biological species that are
considered in biodegradation reactions, including products generated by
biodegradation reactions, nutrients required for biological growth, electron
acceptors, and biological species.
Note: See Appendix D for the component numbering scheme used in UTCHEM
and Section 9 of the UTCHEM Technical Documentation for more details
on the microbiological population model options.
3.1.24 SPNAME(I), for I=1, N
SPNAME(I) - Name of Ith species.

Note: The name of each component may not exceed 8 characters and each
name must be on a separate line of the input file.
3.1.25 ITRU(I), for I=1, NTW (This line is read only if NTW>0 and ITREAC=1)
ITRU(I) - Flag indicating the units of the Ith water tracer.

Possible Values:
1 - Ith tracer units are in volume %
2 - Ith tracer units are in weight %
3.1.26 ICF(KC), for KC=1, N
ICF(KC) - Flag indicating if KCth component is included in the calculations or

not.
Possible Values:
0 - The KCth component is not included in the calculations
1 - The KCth component is included in the calculations
Example: If 11 components are considered but Alcohol 2 is not present, this line
would appear as follows:
1 1 1 1 1 1 1 0 1 1 1
21
UTCHEM User's Guide
User notes:
22
UTCHEM User's Guide
3.2 Output Option Data
The second input section consists of output options. Please remember that
there are seven comment lines at the beginning of this section and that each data line
is preceded by three comment lines.
3.2.1 ICUMTM, ISTOP, IOUTGMS
ICUMTM - Flag indicating if the output intervals indicated by the CUMPR1,

CUMHI1, WRHPV, WRPRF and RSTC variables are specified in pore
volumes or days.
Possible Values:
0 - Data will be written in day intervals
1 - Data will be written in pore volume intervals
Note: The day interval output option (ICUMTM=0) is particularly useful if
there is a shut in period.
ISTOP - Flag indicating if the maximum and injection times (variables TMAX and
TINJ) are specified in pore volumes or days.
Possible Values:
0 - TMAX and TINJ are specified in days
1 - TMAX and TINJ are specified in pore volumes
IOUTGMS - Flag indicating whether or not Ground Water Modeling System

(GMS) files will be created.
Possible Values:
0 - GMS files will not be created
1 - GMS files will be created
Note: IOUTGMS must be set to 0 unless you are running GMS.
3.2.2 IPRFLG(KC), for KC=1, N
IPRFLG(KC) - Flag indicating if profile of KCth component should be written to

CONCP, COMP_AQ, COMP_OIL, COMP_ME, and GMS concentration
output files.
Possible Values:
0 - Profile of KCth component will not be written
1 - Profile of KCth component will be written
Note: If IPCTOT=0, none of the component profiles will be written.
Example: If 11 components are present and only profiles for the oil, surfactant, and
polymer components are desired, this line would appear as follows:
0 1 1 1 0 0 0 0 0 0 0
3.2.3 IPPRES, IPSAT, IPCTOT, IPBIO, IPCAP, IPGEL, IPALK, IPTEMP, IPOBS
IPPRES - Flag indicating if profile of phase pressures should be written to

PRESP output file.
Possible Values:
0 - Profile of phase pressures will not be written
23
UTCHEM User's Guide
1 - Profile of phase pressures will be written
IPSAT - Flag indicating if profile of phase saturations should be written to SATP

output file.
Possible Values:
0 - Profile of phase saturations will not be written
1 - Profile of phase saturations will be written
IPCTOT - Flag indicating if profile of component concentrations should be written

to CONCP output file.
Possible Values:
0 - Profile of component concentrations will not be written
1 - Profile of component concentrations will be written
IPBIO - Flag indicating if chemical and biological species concentrations in the

aqueous phase and within attached biomass should be written to BIO output
file. IPBIO also controls whether aqueous phase concentrations of
components 1 through 8 plus oil components should be written to the
HISTxx output files.
Possible Values:
0 - Aqueous and intra-biomass concentrations will not be written
1 - Aqueous and intra-biomass concentrations will be written
IPCAP - Flag indicating if profile of capacitance (if ICAP=1) or dual porosity (if
ICAP=2) properties should be written to CAPP (capacitance, ICAP=1) or
FRAC (dual porosity, ICAP=2) output files.
Possible Values:
0 - Profile of capacitance or dual porosity properties will not be
written
1 - Profile of capacitance or dual porosity properties will be written
IPGEL - Flag indicating if profile of gel properties should be written to GFILEP

output file.
Possible Values:
0 - Profile of gel properties will not be written
1 - Profile of gel properties will be written
IPALK - Flag indicating if profile of properties related to the alkaline option

should be written to ALKP output file.
Possible Values:
0 - Profile of properties related to the alkaline option will not be
written
1 - Profile of properties related to the alkaline option will be written
IPTEMP - Flag indicating if profile of reservoir temperature should be written to

TEMPP output file.
Possible Values:
0 - Profile of temperature will not be written
1 - Profile of temperature will be written
IPOBS - Flag indicating if aqueous (or gas) phase tracer concentration at

observation points should be written to the TRACxx output files.
Possible Values:
24
UTCHEM User's Guide
0 - Phase tracer concentrations at observation points will not be

written
1 - Tracer concentrations at observation points will be written
3.2.4 ICKL, IVIS, IPER, ICNM, ICSE, IHYSTP, IFOAMP, INONEQ
ICKL - Flag indicating if component concentration data in each phase should be

written to COMP_AQ, COMP_OIL, and COMP_ME output files.
Possible Values:
0 - Component concentration data in each phase will not be written
1 - Component concentration data in each phase will be written
IVIS - Flag indicating if phase viscosities should be written to VIS output file.
Possible Values:
0 - Phase viscosities will not be written
1 - Phase viscosities will be written
IPER - Flag indicating if relative permeabilities should be written to RPERM

output file.
Possible Values:
0 - Relative permeabilities will not be written
1 - Relative permeabilities will be written
ICNM - Flag indicating if phase capillary numbers and interfacial tensions should
be written to TRAP output file.
Possible Values:
0 - Capillary numbers, residual saturations, and interfacial tensions
will not be written
1 - Capillary numbers, residual saturation, and interfacial tensions
will be written
ICSE - Flag indicating if effective salinity should be written to CSE output file.
Possible Values:
0 - Effective salinity information will not be written
1 - Effective salinity information will be written
IHYSTP - Flag indicating if properties related to the hysteresis option should be

written to HYST output file.
Possible Values:
0 - Hysteresis data will not be written
1 - Hysteresis data will be written
IFOAMP - Flag indicating if properties related to the foam option should be

written to FOAM output file.
Possible Values:
0 - Foam data will not be written
1 - Foam data will be written
INONEQ - Flag indicating if properties related to the rate-limited mass transfer

option should be written to DISS output file.
Possible Values:
0 - Rate-limited mass transfer data will not be written
1 - Rate-limited mass transfer data will be written
25
UTCHEM User's Guide
3.2.5 IADS, IVEL, IRKF, IPHSE
IADS - Flag indicating if surfactant, polymer, calcium, gel, chromium, hydrogen,

and sodium adsorption data should be written to PROF output file.
Possible Values:
0 - Adsorption data will not be written
1 - Adsorption data will be written
IVEL - Flag indicating if X, Y, and Z direction phase fluxes should be written to

PROF output file.
Possible Values:
0 - X, Y, and Z direction phase fluxes will not be written
1 - X, Y, and Z direction phase fluxes will be written
IRKF - Flag indicating if permeability reduction factors should be written to

PROF output file.
Possible Values:
0 - Permeability reduction factors, polymer viscosities, and
equivalent shear rate will not be written
1 - Permeability reduction factors, polymer viscosities, and
equivalent shear rate will be written
IPHSE - Flag indicating if phase environment indexing should be written to PROF

output file.
Possible Values:
0 - Phase environment indexing will not be written
1 - Phase environment indexing will be written
Note: The indices for the phase environment are as follows:
1 - single phase
2 - two phase oil/water or oil/microemulsion or
water/microemulsion
3 - three phase oil/microemulsion/water
4 - lobe II(+) of type III
5 - lobe II(-) of type III
3.2.6 NOBS (This line is read only if IPOBS=1)
NOBS - Number of tracer concentration observation points.
3.2.7 IOBS(I), JOBS(I), KOBS(I), for I=1, NOBS (This line is read only if IPOBS=1
and NOBS>0)
IOBS(I) - Index of Ith observation point in X direction.
JOBS(I) - Index of Ith observation point in Y direction.
KOBS(I) - Index of Ith observation point in Z direction.

Note: See the note for input line 3.3.6 for a description of how the gridblocks are
ordered in UTCHEM.
26
UTCHEM User's Guide
3.3. Reservoir Properties
The third input section consists of the reservoir properties. Please
remember that there are seven comment lines at the beginning of this section and
that each data line is preceded by three comment lines.
3.3.1 TMAX
TMAX - Total injection period (maximum simulation time).

Units: days or pore volumes (dependent on value of ISTOP flag)
3.3.2 COMPR, PSTAND
COMPR - Rock compressibility.

Units: 1/psi (IUNIT=0) or 1/kPa (IUNIT=1)
PSTAND - Reference pressure at which pore volume and fluid compressibilities

are specified.
Units: psi (IUNIT=0) or kPa (IUNIT=1)
Reservoir/Aquifer Properties (input lines 3.1.3-3.1.17)
3.3.3 IPOR1, IPERMX, IPERMY, IPERMZ, IMOD
IPOR1 - Flag indicating constant or variable porosity for reservoir.

Possible Values:
0 - Constant porosity for whole reservoir
1 - Constant porosity for each layer
2 - Variable porosity over reservoir
IPERMX - Flag indicating constant or variable X direction permeability

(ICOORD=1 or 3) or radial direction permeability (ICOORD=2) for
reservoir.
Possible Values:
0 - Constant permeability for whole reservoir
1 - Constant permeability for each layer in the X direction
(ICOORD=1 or 3) or radial direction (ICOORD=2)
2 - Variable permeability over reservoir
IPERMY - Flag indicating constant or variable Y direction permeability for

reservoir.
Possible Values:
1 - Constant permeability for each layer in the Y direction
3 - Y direction permeability is dependent on X direction
permeability
IPERMZ - Flag indicating constant or variable Z direction permeability for

reservoir.
Possible Values:
27
UTCHEM User's Guide

1 - Constant permeability for each layer in the Z direction
3 - Z direction permeability is dependent on X direction
permeability
IMOD - Flag indicating whether the reservoir properties are modified or not.
Possible Values:
0 - No property is modified
1 - Allow for property modification
the porosity and permeability values for different options:
IPOR1 (porosity flag) 0 (constant porosity for reservoir) 3.3.4

1 (constant porosity for each layer) 3.3.5
2 (variable porosity over reservoir) 3.3.6
IPERMX (X or radial 0 (constant permeability for reservoir) 3.3.7

direction permeability 1 (constant permeability for each layer) 3.3.8
flag) 2 (variable permeability over reservoir) 3.3.9
IPERMY (Y direction 0 (constant permeability for reservoir) 3.3.10

permeability flag; 1 (constant permeability for each layer) 3.3.11
only if ICOORD≠2) 2 (variable permeability over reservoir) 3.3.12
3 (permeability dependent on X direction) 3.3.13
IPERMZ (Z direction 0 (constant permeability for reservoir) 3.3.14

permeability flag) 1 (constant permeability for each layer) 3.3.15
2 (variable permeability over reservoir) 3.3.16
3 (permeability dependent on X direction) 3.3.17
3.3.4 PORC1 (This line is read only if IPOR1=0)
PORC1 - Reservoir porosity.

Units: fraction
Note: All elements of the POR array will be set equal to PORC1.
3.3.5 POR(K), for K=1, NZ (This line is read only if IPOR1=1)
POR(K) - Porosity of Kth layer.

Units: fraction
Note: NZ values are actually read into a workspace array (WKSP1) and
then the first set of NX × NY elements (corresponding to layer 1) of the
POR array are set equal to WKSP1(1), the second set of NX × NY
elements (corresponding to layer 2) of the POR array are set equal to
WKSP1(2), etc.
28
UTCHEM User's Guide
3.3.6 POR(I), for I=1, NX × NY × NZ (This line is read only if IPOR1=2)
POR(I) - Porosity of Ith gridblock

Units: fraction
Note: The three-dimensional grid system is being read into a one-
dimensional array. The first index (column) of the three-
dimensional system varies fastest, the second index (row) varies
next fastest, and the third index (layer) varies slowest.
Example: If you had a 4 × 3 × 2 system (4 columns—NX=4, 3
rows—NY=3, and 2 layers—NZ=2), the values would be read in
the following order:
1,1,1 2,1,1 3,1,1 4,1,1
1,2,1 2,2,1 3,2,1 4,2,1
1,3,1 2,3,1 3,3,1 4,3,1
1,1,2 2,1,2 3,1,2 4,1,2
1,2,2 2,2,2 3,2,2 4,2,2
1,3,2 2,3,2 3,3,2 4,3,2
UTCHEM sets the transmissibilities to zero for gridblocks with

porosity values less than or equal to 0.01 to specify inactive cells
(ICOORD=1).
3.3.7 PERMXC (This line is read only if IPERMX=0)
PERMXC - Permeability of the reservoir in the X direction or in the radial

direction (ICOORD=2).
Units: millidarcies = 10-3 µm2
Note: All elements of the PERMX array will be set equal to PERMXC.
Note: In order to specify clay gridblocks, the user should set the permeability for
those gridblocks to a value less than or equal to 10-5 md. UTCHEM then
sets the porosity to 0.01 and water saturation to 1.0 for those gridblocks.
3.3.8 PERMX(K), for K=1, NZ (This line is read only if IPERMX=1)
PERMX(K) - Permeability of the Kth layer in the X direction or in the radial

Note: See the note for input line 3.3.5 for information on how the data is
read. See the note for input line 3.3.7 for information on specifying
clay gridblocks.
3.3.9 PERMX(I), for I=1, NX × NY × NZ (This line is read only if IPERMX=2)
PERMX(I) - Permeability of the Ith gridblock in the X direction or in the radial

29
UTCHEM User's Guide
Note: See the note and example for input line 3.3.6 for the order of the
permeability values. See the note for input line 3.3.7 for
information on specifying clay gridblocks.
3.3.10 PERMYC (This line is read only if IPERMY=0 and ICOORD≠2)
PERMYC - Permeability of the reservoir in the Y direction.

Note: All elements of the PERMY array will be set equal to PERMYC.
3.3.11 PERMY(K), for K=1, NZ (This line is read only if IPERMY=1 and ICOORD≠2)
PERMY(K) - Permeability of the Kth layer in the Y direction.

read.
3.3.12 PERMY(I), for I=1, NX × NY × NZ (This line is read only if IPERMY=2 and
ICOORD≠2)
PERMY(I) - Permeability of the Ith gridblock.

permeability values.
3.3.13 FACTY (This line is read only if IPERMY=3 and ICOORD≠2)
FACTY - Constant permeability multiplier for Y direction permeability.

Units: dimensionless
Note: The X direction permeabilities are multiplied by FACTY to obtain
the Y direction permeabilities.
3.3.14 PERMZC (This line is read only if IPERMZ=0)
PERMZC - Permeability of the reservoir in the Z direction.

Note: All elements of the PERMZ array will be set equal to PERMZC.
3.3.15 PERMZ(K), for K=1, NZ (This line is read only if IPERMZ=1)
PERMZ(K) - Permeability of the Kth layer in the Z direction.

read.
3.3.16 PERMZ(I), for I=1, NX × NY × NZ (This line is read only if IPERMZ=2)
PERMZ(I) - Permeability of the Ith gridblock.

Units: millidarcies (10-3 µm2)
30
UTCHEM User's Guide
permeability values.
3.3.17 FACTZ (This line is read only if IPERMZ=3)
FACTZ - Constant permeability multiplier for Z direction permeability.

Note: The X direction permeabilities are multiplied by FACTZ to obtain
the Z direction permeabilities.
Initial Reservoir/Aquifer Data (input lines 3.3.18-3.3.35)
3.3.18 IDEPTH, IPRESS, ISWI, ICWI
IDEPTH - Flag indicating type of depth measurement of the top layer.

Possible Values:
0 - Single value for depth of the top layer is specified
1 - Depth of top gridblock (1,1,1) and the reservoir dip angles are
specified
2 - Depth of each gridblock in the top layer is specified
Note: If ICOORD=2, this value is automatically set equal to 0. The depth
is specified at the middle of a gridblock.
IPRESS - Flag indicating type of reservoir initial pressure measurement.

Possible Values:
0 - Single value for reservoir initial pressure is used for all
gridblocks
1 - Initial pressure for a point at a specified depth is specified by
user
2 - Initial pressure for each gridblock is specified by user
ISWI - Flag indicating type of initial water saturation measurement.

Possible Values:
0 - Single value for initial water saturation is used for all
gridblocks
1 - Constant value for water saturation for each layer is specified
by user
2 - Initial water saturation for each gridblock is specified by user
ICWI - Flag indicating type of initial aqueous phase compositions.

Possible Values:
-1 - Backward compatibility with previous versions of UTCHEM
(only species in aqueous phase are anion and cation)
0 - Constant initial aqueous phase concentration for whole
reservoir
1 - Constant initial aqueous phase concentration for each layer
2 - Initial aqueous phase concentration specified for each
gridblock
31
UTCHEM User's Guide
the initial properties such as depth, pressure, and initial water saturations:
IDEPTH (depth flag) 0 (single value for top layer) 3.3.19

1 (top gridblock and dip angles) 3.3.20
2 (each gridblock in top layer) 3.3.21
IPRESS (initial pressure 0 (single value for reservoir) 3.3.22

flag) 1 (one point at a specified depth) 3.3.23
2 (each gridblock) 3.3.24
ISWI (water saturation 0 (single value for reservoir) 3.3.25

flag) 1 (constant value for each layer) 3.3.26
ICWI (aqueous phase 0 (single value for reservoir) 3.3.32

composition flag) 1 (constant value for each layer) 3.3.33
3.3.19 D111 (This line is read only if IDEPTH=0)
D111 - Depth of the top layer of the reservoir measured from the surface
(reference plane), positive downward.
Note: If IDEPTH=0 and ICOORD=4, D111 is the reference depth of the first
gridblock.
3.3.20 D111, THETAX, THETAY (This line is read only if IDEPTH=1)
D111 - Depth of the first gridblock (1,1,1).

THETAX - Reservoir dip angle in X direction, positive downward.

Units: radians
THETAY - Reservoir dip angle in Y direction, positive downward.

Units: radians
Note: If ICOORD=4, set THETAY equal to 0 (dip angle in X-Z plane).
3.3.21 EL(I), for I=1, NX × NY (This line is read only if IDEPTH=2)
EL(I) - Depth of Ith gridblock in the top layer (K=1).

gridblock depths.
32
UTCHEM User's Guide
3.3.22 PRESS1 (This line is read only if IPRESS=0)
PRESS1 - Initial reservoir pressure.

3.3.23 PINIT, HINIT (This line is read only if IPRESS=1)
PINIT - Initial reservoir pressure at HINIT depth.

Units: psia (IUNIT=0) or kPa (IUNIT=1)
HINIT - Depth of the point where PINIT initial pressure is specified.

Note: Initial pressure is assumed to be the aqueous phase pressure.
3.3.24 P(I), for I=1, NX × NY × NZ (This line is read only if IPRESS=2)
P(I) - Initial pressure of each gridblock in the reservoir.

initial pressure values. This is assumed to be the aqueous phase
pressure.
3.3.25 SWI (This line is read only if ISWI=0)
SWI - Initial water saturation for all gridblocks of the reservoir.

Units: fraction of pore volume
3.3.26 S(K,1), for K=1, NZ (This line is read only if ISWI=1)
S(K,1) - Initial water saturation for Kth layer.

read.
3.3.27 S(I,1), for I=1, NX × NY × NZ (This line is read only if ISWI=2)
S(I,1) - Initial water saturation for Ith block.

initial water saturation values.
Initial Gas Saturations (input lines 3.3.28-3.3.31) — This section is required

only if gas is present (IGAS≥1).
3.3.28 ISGI (This line is read only if IGAS≥1)
ISGI - Flag indicating type of initial gas saturation.

Possible Values:
0 - Constant initial gas saturation for whole reservoir
1 - Constant initial gas saturation for each layer
2 - Initial gas saturation specified for each gridblock
33
UTCHEM User's Guide
3.3.29 SGI (This line is read only if IGAS≥1 and ISGI=0)
SGI - Initial gas saturation for all gridblocks of the reservoir.

3.3.30 S(K,4), for K=1, NZ (This line is read only if IGAS≥1 and ISGI=1)
S(K,4) - Initial gas saturation for Kth layer.

read.
3.3.31 S(I,4), for I=1, NX × NY × NZ (This line is read only if IGAS≥1 and ISGI=2)
S(I,4) - Initial gas saturation for Ith gridblock.

initial gas saturation values.
Initial Aqueous Phase Compositions (input lines 3.3.32-3.3.34)
Note: Please see the example following input line 3.3.34 for more information on how
UTCHEM uses the ICF flag to control exactly how many values are read on input
lines 3.3.32-3.3.34.
3.3.32 CWI(KW), for KW=1, N << or >>

CWI(KW), for KW=1, 8+NO (This line is read only if ICWI=0)
CWI(KW) - Initial concentration of species KW in the aqueous phase.

Units: volume fraction
Note: If N≤8, at most N values are required by the program. If N>8, at most
8+NO values are required by the program. Please see example following
input line 3.3.34.
3.3.33 CWI(K,KW), for K=1, NZ, for KW=1, N << or >>

CWI(K,KW), for K=1, NZ, for KW=1, 8+NO (This line is read only if
ICWI=1)
CWI(K, KW) - Initial concentration of species KW in the aqueous phase at Kth

layer.
Note: If N≤8, at most N×NZ values are required by the program. If N>8, at
most (8+NO)×NZ values are required by the program. Please see example
following input line 3.3.34.
34
UTCHEM User's Guide
3.3.34 CWI(I,KW), for I=1, NX × NY × NZ, for KW=1, N << or >>

CWI(I,KW), for I=1, NX × NY × NZ, for KW=1, 8+NO (This line is read
only if ICWI=2)
CWI(I,KW) - Initial concentration of species KW in the aqueous phase at Ith

gridblock.
Note: If N≤8, at most N×NX×NY×NZ values are required by the program. If

N>8, at most (8+NO) ×NX×NY×NZ values are required by the program.
Please see the example below.
Example:
Values are only read for the 8 "basic" components (water, oil, surfactant,
polymer, anion, cation, alcohol 1, alcohol 2) and any multiple organic
components present. See Appendix D for more information on UTCHEM
components.
In addition, values are only read for a component if the ICF flag for that
component has been set to 1.
In this example, we are providing initial concentrations for the water, oil,
and anion components for a case with 5 components and 5×3×2 gridblocks
where the ICF array has been initialized as:
1 1 0 0 1
For ICWI=0:
CC
CC
*-- Initial water conc. in aq. phase
0.999
CC
CC
*-- Initial oil conc. in aq. phase
0.001
CC
CC
*-- Initial chloride conc. in aq. phase
0.4
For ICWI=1:
CC
CC
*-- Initial water conc. in aq. phase in each layer
0.999 0.999
CC
CC
*-- Initial oil conc. in aq. phase in each layer
0.001 0.001
CC
35
UTCHEM User's Guide
CC
*-- Initial chloride conc. in aq. phase in each layer
0.4 0.4
For ICWI=2:
CC
CC
*-- Initial water conc. in aq. phase in each gridblock
30*0.999
CC
CC
*-- Initial oil conc. in aq. phase in each gridblock
30*0.001
CC
CC
*-- Initial chloride conc. in aq. phase in each gridblock
15*0.4 15*0.3
Initial Organic Concentrations (input lines 3.3.35-3.3.38) — This section is

required only if the multiple organic option is used (NO>1).
3.3.35 ICOI (This line is read only if NO>1)
ICOI - Flag indicating type of initial oil phase compositions.

Possible Values:
0 - Constant initial oil phase concentration for whole reservoir
1 - Constant initial oil phase concentration for each layer
2 - Initial oil phase concentration specified for each gridblock
3.3.36 COI(KO), for KO=1, NO (This line is read only if NO>1 and ICOI=0)
COI(KO) - Initial oil phase concentration for oil component KO for the reservoir.
3.3.37 COI(K,KO), for K=1, NZ, for KO=1, NO (This line is read only if NO>1 and
ICOI=1)
COI(K, KO) - Initial oil phase concentration for oil component KO at Kth layer.
3.3.38 COI(I,KO), for I=1, NX × NY × NZ, for KO=1,NO (This line is read only if
NO>1 and ICOI=2)
COI(I,KO) - Initial oil phase concentration for oil component KO at Ith gridblock.
36
UTCHEM User's Guide
Reservoir Property Modification Data (input lines 3.3.39-3.3.49) — This

section is required only if IMOD=1.
3.3.39 IMPOR, IMKX, IMKY, IMKZ, IMSW (This line is read only if IMOD=1)
IMPOR - Flag indicating whether the porosity is modified or not.

Possible Values:
0 - No modification is considered in porosity values
1 - Allow modification in porosity
IMKX- Flag indicating whether the permeability in the X direction is modified or

not.
Possible Values:
0 - No modification is considered in X permeability
1 - Allow modification in X permeability
IMKY- Flag indicating whether the permeability in the Y direction is modified or

not.
Possible Values:
0 - No modification is considered in Y permeability
1 - Allow modification in Y permeability
IMKZ- Flag indicating whether the permeability in the Z direction is modified or

not.
Possible Values:
0 - No modification is considered in Z permeability
1 - Allow modification in Z permeability
IMSW- Flag indicating whether the initial water saturation is modified or not.
Possible Values:
0 - No modification is considered in initial water saturation values
1 - Allow modification in initial water saturation
3.3.40 NMOD0 (This line is read only if IMOD=1 and IMPOR=1)
NMOD0 - Number of regions with modified porosity.
3.3.41 IMIN, IMAX, JMIN, JMAX, KMIN, KMAX, IFACT, FACTX (This line is
read only if IMOD=1 and NMOD>0)
IMIN - The first index in X direction.
IMAX - The last index in X direction.
JMIN - The first index in Y direction.
JMAX - The last index in Y direction.
KMIN - The first index in Z direction.
KMAX - The last index in Z direction.
IFACT - Flag indicating how porosity is modified.

Possible Values:
37
UTCHEM User's Guide
1 - Replace porosity with FACTX

2 - Multiply porosity by FACTX
3 - Add FACTX to porosity
FACTX - The constant used to modify the porosity value.
ordered in UTCHEM. This line is repeated NMOD0 times.
3.3.42 NMOD1 (This line is read only if IMOD=1 and IMKX=1)
NMOD1 - Number of regions with modified X permeability.
read only if IMOD=1 and NMOD1>0)
IFACT - Flag indicating how X permeability is modified.

Possible Values:
1 - Replace X permeability with FACTX
2 - Multiply X permeability by FACTX
3 - Add FACTX to X permeability
FACTX - The constant used to modify the X permeability value.
3.3.44 NMOD2 (This line is read only if IMOD=1 and IMKY=1)
NMOD2 - Number of regions with modified Y permeability.
IMIN - The first index in direction.
38
UTCHEM User's Guide
IFACT - Flag indicating how Y permeability is modified.

Possible Values:
1 - Replace Y permeability with FACTX
2 - Multiply Y permeability by FACTX
3 - Add FACTX to Y permeability
FACTX - The constant used to modify the Y permeability value.
3.3.46 NMOD3 (This line is read only if IMOD=1 and IMKZ=1)
NMOD3 - Number of regions with modified Z permeability.
IFACT - Flag indicating how Z permeability is modified.

Possible Values:
1 - Replace Z permeability with FACTX
2 - Multiply Z permeability by FACTX
3 - Add FACTX to Z permeability
FACTX - The constant used to modify the Z permeability value.
3.3.48 NMOD4 (This line is read only if IMOD=1 and IMSW=1)
NMOD4 - number of regions with modified initial water saturation.
39
UTCHEM User's Guide
IFACT - Flag indicating how initial water saturation is modified.

Possible Values:
1 - Replace initial water saturation with FACTX
2 - Multiply initial water saturation by FACTX
3 - Add FACTX to initial water saturation
FACTX - The constant used to modify the initial water saturation value.
3.3.50 NMOD5 (This line is read only if IMOD=1, IMSW=1, and IGAS≥1)
NMOD5 - number of regions with modified initial gas saturation.
IFACT - Flag indicating how initial gas saturation is modified.

Possible Values:
1 - Replace initial gas saturation with FACTX
2 - Multiply initial gas saturation by FACTX
3 - Add FACTX to initial gas saturation
FACTX - The constant used to modify the initial gas saturation value.
3.3.52 C50, C60 (This line is read only if ICWI=-1)
C50 - Constant initial brine salinity.
40
UTCHEM User's Guide
Units: meq/ml of water

Note: This is assumed to be all the anions (in equivalents).
C60 - Initial divalent cation concentration of brine.

Units: meq/ml of water
Note: C50 and C60 are replaced by the input values of C5I and C6I when
IREACT>1.
41
UTCHEM User's Guide
User notes:
42
UTCHEM User's Guide
3.4 General Physical Property Data
The fourth input section consists of the general physical property data.
Please remember that there are seven comment lines at the beginning of this section
and that each data line is preceded by three comment lines.
Surfactant/Cosolvent Phase Behavior Data (input lines 3.4.1-3.4.21)
3.4.1 C2PLC, C2PRC, EPSME, IHAND
C2PLC - Oil concentration at plait point in type II(+) region.

C2PRC - Oil concentration at plait point in type II(-) region.

EPSME - Critical micelle concentration (CMC)—minimum surfactant

concentration for the formation of micelles.
IHAND - Flag to specify whether modified Hand's rule is considered or not.

Possible Values:
0 - Original Hand's rule is considered for phase behavior (default)
1 - Modified Hand's rule is considered for Phase behavior
Note: The option of IHAND=1 is available only for oil/microemulsion,
Type II(-) phase behavior, and IMASS=1.
3.4.2 IFGHBN
IFGHBN - Flag indicating type of phase behavior parameters.

Possible Values:
0 - Input height of binodal curve (default)
1 - Input solubilization ratio (new option)
Note: The input solubilization parameter (IFGHBN=1) is currently only available

for the multiple organic option (NO>1). The effect of temperature or
alcohol on phase behavior is not currently modeled for IFGHBN=1. See
Sections 2 & 11 of the UTCHEM Technical Documentation for more details
on the input height of binodal curve option and Section 7 of the UTCHEM
Technical Documentation for more details on the input solubilization ratio
option.
Binodal Curve Input Option (input lines 3.4.3-3.4.11) — This section is

required only if IFGHBN=0.
3.4.3 HBNS70, HBNC70, HBNS71, HBNC71, HBNS72, HBNC72 (This line is read
only if IFGHBN=0)
HBNS70 - Slope for maximum height of binodal curve vs. fraction of Alcohol 1
associated with surfactant at zero salinity.
43
UTCHEM User's Guide
HBNC70 - Intercept of maximum height of binodal curve at zero fraction of

Alcohol 1 (associated with surfactant at zero salinity.
associated with surfactant at optimal salinity.

Alcohol associated with surfactant at optimal salinity.
associated with surfactant at twice optimal salinity.

Alcohol 1 associated with surfactant at twice optimal salinity.
Note: If alcohol is not present, the maximum height of binodal curves at
three different salinities are the only parameters used in the phase
behavior calculations.
3.4.4 HBNT0, HBNT1, HBNT2, CSET (This line is read only if IFGHBN=0 and
IENG=1)
HBNT0 - Slope of height of binodal curve versus temperature at zero salinity

Units: volume fraction/(˚F) (IUNIT=0) or volume fraction/˚C (IUNIT=1)
HBNT1 - Slope of height of binodal curve versus temperature at optimal salinity

Units: volume fraction /(˚F) (IUNIT=0) or volume fraction/˚C (IUNIT=1)
HBNT2 - Slope of height of binodal curve versus temperature at twice optimal

salinity
Units: volume fraction /(˚F) (IUNIT=0) or volume fraction/˚C (IUNIT=1)
CSET - The Slope parameter, βT, for temperature dependency of the three-phase
window
Units: (˚F)-1 (IUNIT=0) or (˚C)-1 (IUNIT=1)
3.4.5 HBNS80, HBNC80, HBNS81, HBNC81, HBNS82, HBNC82 (This line is read
only if IFGHBN=0)
associated with surfactant at zero salinity.

Alcohol 2 associated with surfactant at zero salinity.
44
UTCHEM User's Guide
HBNS81 - Slope of maximum height of binodal curve vs. fraction of Alcohol 2

associated with surfactant at optimal salinity.

Alcohol 2 associated with surfactant at optimal salinity.
associated with surfactant at twice optimal salinity.

Alcohol 2 associated with surfactant at twice optimal salinity.
3.4.6 CSEL7, CSEU7, CSEL8, CSEU8 (This line is read only if IFGHBN=0)
CSEL7 - Lower effective salinity limit for type III phase region determined when
Alcohol 1 and calcium approach zero.
Units: meq/ml water
CSEU7 - Upper effective salinity limit for type III phase region determined when
Units: meq/ml water
CSEL8 - Lower effective salinity limit for type III phase region determined when
Units: meq/ml water
CSEU8 - Upper effective salinity limit for type III phase region determined when
Units: meq/ml water
3.4.7 BETA6, BETA7, BETA8 (This line is read only if IFGHBN=0)
BETA6 - The effective salinity slope parameter for calcium.

BETA7 - The effective salinity slope parameter for Alcohol 1.

BETA8 - The effective salinity slope parameter for Alcohol 2.

3.4.8 IALC, OPSK7O, OPSK7S, OPSK8O, OPSK8S (This line is read only if
IFGHBN=0)
IALC - Flag indicating choice of alcohol partition model to use.

Possible Values:
0 - Hirasaki's model will be used
1 - Prouvost's model will be used
OPSK7O - Alcohol partition coefficient (oil/water) for Alcohol 1.
45
UTCHEM User's Guide
OPSK7S - Alcohol partition coefficient (surfactant/water) for Alcohol 1.

OPSK8O - Alcohol partition coefficient (oil/water) for Alcohol 2.

OPSK8S - Alcohol partition coefficient (surfactant/water) for Alcohol 2.

Note: If IALC=0 then OPSK7O, OPSK7S, OPSK8O, and OPSK8S remain

fixed. If OPSK7O, OPSK7S, OPSK8O, and OPSK8S are equal to zero
and IALC=0, then alcohol is lumped with surfactant as a single component
(total chemical). OPSK7O, OPSK7S, OPSK8O, and OPSK8S are only
used when Hirasaki's model is chosen. See Section 11 of the UTCHEM
Technical Documentation for more details on the alcohol partition models.
3.4.9 NALMAX, EPSALC (This line is read only if IFGHBN=0)
NALMAX - Maximum number of iterations for alcohol partitioning for two

alcohol system.
Note: The suggested value is 20 and a value of zero would result in no
iterations.
EPSALC - Tolerance for convergence of iterations for two alcohol system.

Note: Suggested values are 10-3 and 10-4.
3.4.10 AKWC7, AKWS7, AKM7, AK7, PT7 (This line is read only if IFGHBN=0)
AKWC7, AKWS7 - Parameters used to determine partition coefficient of

monomeric Alcohol 1 between aqueous and oleic pseudophases.
AKM7 - Partition coefficient of monomeric Alcohol 1 between surfactant and

oleic pseudophases.
AK7 - Self-association constant of Alcohol 1 in oleic pseudophase.

PT7 - Ratio of molar volume of Alcohol 1 to equivalent molar volume of

surfactant.
3.4.11 AKWC8, AKWS8, AKM8, AK8, PT8 (This line is read only if IFGHBN=0)
AKWC8, AKWS8 - Parameters used to determine partition coefficient of

monomeric Alcohol 2 between aqueous and oleic pseudophases.
AKM8 - Partition coefficient of monomeric Alcohol 2 between surfactant and

oleic pseudophases.
46
UTCHEM User's Guide
AK8 - Self-association constant of Alcohol 2 in oleic pseudophase.

PT8 - Ratio of molar volume of Alcohol 2 to equivalent molar volume of

surfactant.
Solubilization Ratio Input Option (input lines 3.4.12-3.4.21) — These lines

are required only if the phase behavior calculation is based on the solubilization ratio
(IFGHBN=1).
3.4.12 IOD (This line is read only if IFGHBN=1 and NO>1)
IOD - Flag indicating whether phase behavior depends on organic composition

Possible Values:
0 - Phase behavior and properties depend on organic composition
1 - Phase behavior and properties are independent of organic
composition (default)
3.4.13 NCOMP (This line is read only if IFGHBN=1 and IOD=0)
NCOMP - Number of organic components in the first solubility measurement.
47
UTCHEM User's Guide
3.4.14 ONAME(I) for I=1, NCOMP (This line is read only if IFGHBN=1 and IOD=0)
ONAME(I) - Name of organic species in the first solubility measurement.

Currently the values in the following table can be specified.
ONAME Formula Name Molecular Equivalent

Value Weight Alkane
Carbon No.
(EACN)
DECANE C10H22 Decane 142 10
OCTANE C8H18 Octane 114 8
HEXANE C6H14 Hexane 86 6
PCE C2Cl4 Tetrachloro- 165.8 2.9
ethylene
PXYLEN P-xylene 106 2
TOLUEN Toluene 92 1
CCL4 CCl4 Carbon 153.8 -0.06
tetrachloride
TCE C2HCl3 Trichloro- 131.4 -3.81
ethylene
DCB 1,2-C6H4Cl2 1,2-dichloro- 146.9 -4.89
benzene
TCA CH3CCl3 1,1,1-trichloro- 133.35 -2.5
ethene
DCE 1,2-C2H4Cl2 98.9 -12.1
CHCL3 CHCl3 119.4 -13.67
CH2CL2 CH2Cl2 84.9 -13.79
C2CL4 1,1,2,2- 167.8 -22.15
C2H2Cl4
Note: Each organic species name must appear on a separate line of the input file.
3.4.15 OCOMP(I), for I=1, NCOMP (This line is read only if IFGHBN=1 and IOD=0)
OCOMP(I) - Concentration of Ith organic component in the first solubility

measurement.
Units: mole fraction
3.4.16 CS0, SCS0, CS1, SCS1, CS2, SCS2, DCS20 (This line is read only if
IFGHBN=1)
CS0 - Effective salinity which is between the lower and optimal effective salinity
limits for type III phase region in the first solubility measurement;
CSEL<CS0<CSEOP.
Units: meq/ml
SCS0 - Solubilization parameter at CS0 in the first solubility measurement.

48
UTCHEM User's Guide
CS1 - Optimal effective salinity in the first solubility measurement (CSEOP).

Units: meq/ml
SCS1 - Solubilization parameter at CS1.

CS2 -Effective salinity which is between the optimal and the upper effective
salinity limits for type III phase region in the first solubility measurement;
CSEOP<CS2<CSEU.
Units: meq/ml
SCS2 - Solubilization parameter at CS2 in the first solubility measurement.

DCS20 - The difference of the upper and the lower effective salinity limits for
type III phase region in the first solubility measurement; CSEU - CSEL.
3.4.17 NCOMP (This line is read only if IFGHBN=1 and IOD=0)
NCOMP - Number of organic components in the second solubility measurement.
3.4.18 ONAME(I), for I=1, NCOMP (This line is read only if IFGHBN=1 and IOD=0)
ONAME(I) - Name of organic species in the second solubility measurement.

Note: See input line 3.4.14 for a list of available species names. Each
organic species name must appear on a separate line of the input
file.
3.4.19 OCOMP(I), for I=1, NCOMP (This line is read only if IFGHBN=1 and IOD=0)
OCOMP(I) - Concentration of Ith organic component in the second solubility

measurement.
Units: mole fraction
3.4.20 CS0, SCS0, CS1, SCS1, CS2, SCS2, DCS20 (This line is read only if
IFGHBN=1 and IOD=0)
CS0 - Effective salinity which is between the lower and optimal effective salinity
limits for type III phase region in the second solubility measurement;
CSEL<CS0<CSEOP.
Units: meq/ml
SCS0 - Solubilization parameter at CS0 in the second solubility measurement.

CS1 - Optimal effective salinity in the second solubility measurement.

Units: meq/ml

49
UTCHEM User's Guide
CS2 -Effective salinity which is between the optimal and the upper effective
salinity limits for type III phase region in the second solubility
measurement; CSEOP<CS2<CSEU.
Units: meq/ml

DCS20 - The difference of the upper and the lower effective salinity limits for
type III phase region in the second solubility measurement; CSEU -
C SEL .
3.4.21 BETA6 (This line is read only if IFGHBN=1)
BETA6 - The effective salinity slope parameter for calcium.

Interfacial Tension Data (input lines 3.4.22-3.4.25)
3.4.22 IFT
IFT - This flag indicates which interfacial tension correlation is used.

Possible Values:
0 - Healy and Reed's correlation will be used
1 - Huh's correlation will be used
Note: See Section 2 of the UTCHEM Technical Documentation for more details
on the interfacial tension model options.
3.4.23 G11, G12, G13, G21, G22, G23 (This line is read only if IFT=0)
G11, G12, G13 - Interfacial tension parameters for water-microemulsion

interface.
G21, G22, G23 - Interfacial tension parameters for oil-microemulsion interface.

Units: Dimensionless
3.4.24 CHUH, AHUH (This line is read only if IFT=1)
CHUH - Constant in modified Huh's interfacial tension correlation.

Typical Values: 0.1 - 0.35
AHUH - Constant in modified Huh's interfacial tension correlation

Typical values: 5 - 20
3.4.25 XIFTW
XIFTW - log10 σwo where σwo is the interfacial tension of the water-oil interface.
Units: dynes/cm = mN/m
50
UTCHEM User's Guide
Organic Mass Transfer Data (input lines 3.4.26-3.4.30)
3.4.26 IMASS, ICOR
IMASS - Flag indicating the choice of oil solubility in water.

Possible Values
0 - No solubility of oil in water in the absence of surfactant
1 - Allow for solubility of oil in water in the absence of surfactant
2- Allow for non-equilibrium solubility of oil in water
ICOR - Flag indicating whether mass transfer coefficient is constant or calculated.

Possible Values
0 - Mass transfer coefficient in constant
1 - Mass transfer coefficient in calculated using correlation
on the organic dissolution model used in UTCHEM.
3.4.27 WSOL (This line is read only if NO≤1 and IMASS>0)
WSOL - Equilibrium concentration of oil in water in the absence of surfactant.

3.4.28 WSOLK(K), for K=1, NO (This line is read only if NO>1 and IMASS>1)
WSOLK(K) - Water/oil equilibrium partition coefficient for oil component K in

the absence of surfactant
on the rate limited multiple organic dissolution model.
3.4.29 ISOL, CNEM (This line is read only if IMASS=2 and ICOR=0)
ISOL - Flag indicating the solution scheme for the non-equilibrium mass transfer
calculations
Possible Values
0 - Implicit method is used
1 - Explicit method is used
Note: The explicit method (ISOL=1) is the only option available when gas
is present (IGAS≥1).
CNEM - Coefficient of non-equilibrium mass transfer of oil in aqueous phase

with or without surfactant.
Units: vol. of water/(bulk vol.-day)
Note: The non-equilibrium mass transfer calculation is valid for type II(-)
with the plait point in the corner (C2PLC=0) and in the absence of
gas phase (IGAS=0).
3.4.30 ISOL, EQB0, EQB1, EQB2, EQB3 (This line is read only if IMASS=2 and
ICOR=1)
ISOL - Flag indicating the solution scheme for the non-equilibrium mass transfer
calculations
51
UTCHEM User's Guide
Possible Values
0 - Implicit method is used
1 - Explicit method is used
Note: The explicit method (ISOL=1) is the only option available when gas
is present (IGAS≥1).
EQB0 - Sherwood number constant.
EQB1 - Reynold's number exponent.
EQB2 - Oil content exponent.
EQB3 - Schmidt number exponent.
3.4.31 ITRAP, T11, T22, T33
ITRAP - Flag indicating whether residual saturations and relative permeabilities

are dependent on capillary number or not.
Possible Values:
0 - Residual saturations are not dependent on capillary number;
endpoint and exponent of relative permeability curves are
constant
1 - Residual saturations and relative permeabilities are dependent
on capillary number
2 - Residual saturations and relative permeabilities are dependent
on trapping number
Note: ITRAP=2 is currently not available with the curvilinear grid option
(ICOORD=4) or when gas is present (IGAS≥1)
T11 - Capillary desaturation curve parameter for aqueous phase.
T22 - Capillary desaturation curve parameter for oleic phase.
T33 - Capillary desaturation curve parameter for microemulsion phase.
Note: Options ITRAP=1 and ITRAP=2 are identical for 1-d displacement in the
vertical direction with zero capillary pressure. See Section 2 of the
UTCHEM Technical Documentation for more information on the capillary
and trapping number options.
Relative Permeability Data (input lines 3.4.32-3.4.60)
3.4.32 IPERM
IPERM - Flag indicating which relative permeability and capillary pressure model
is used.
Possible Values:
0 - Imbibition Corey
1 - First drainage Corey (only for two phase water/oil flow)
2- Parker and Lenhard's model for water-wet porous media with
hysteresis option
3 - Lenhard's model for mixed-wet porous media including
hysteresis
52
UTCHEM User's Guide
Note: Lenhard's model (IPERM=3) is for two phase oil/water only. See
Section 14 of the UTCHEM Technical Documentation for more
details on this option.
3.4.33 IHYST (This line is read only if IPERM=2)
IHYST - Flag indicating whether the hysteresis is used with Parker and
Lenhard's model (IPERM=2)
Possible Values:
0 - Hysteretic model is not used
1 - Hysteretic model is used
Note: See Section 3 of the UTCHEM Technical Documentation for more
details on this option.
3.4.34 IPARK (This line is read only if IPERM=2 and IHYST=1)
IPARK - Flag indicating the model used to calculate the oil trapping for hysteretic
model
Possible Values:
0 - Kalurachchi and Parker's model is used
1 - Parker and Lenhard's model is used
Low Capillary Number Data (input lines 3.4.35-3.4.57)
3.4.35 ISRW, IPRW, IEW (This line is read only for IPERM< 2)
ISRW - Flag indicating type of residual saturation.

Possible Values:
0 - Constant residual saturation for entire reservoir
1 - Constant residual saturation for each layer
2 - Residual saturation specified for each gridblock
IPRW - Flag indicating type of endpoint relative permeability.

Possible Values:
0 - Constant endpoint relative permeability for entire reservoir
1 - Constant endpoint relative permeability for each layer
2 - Endpoint relative permeability specified for each gridblock
IEW - Flag indicating type of relative permeability exponent.

Possible Values:
0 - Constant relative permeability exponent for entire reservoir
1 - Constant relative permeability exponent for each layer
2 - Relative permeability exponent specified for each gridblock
3.4.36 ISRW (This line is read only if IPERM=2 or 3)
ISRW - Flag indicating type of residual saturation.

Possible Values:
0 - Constant residual saturation for entire reservoir
1 - Constant residual saturation for each layer
2 - Residual saturation specified for each gridblock
53
UTCHEM User's Guide
the relative permeability parameters for different options:
IPERM = 0 or 1 ISRW [residual saturation]

(imbibition Corey or first 0 (constant for reservoir) 3.4.37
drainage Corey) 1 (constant for each layer) 3.4.38-3.4.40
2 (each gridblock) 3.4.41-3.4.43
IPRW [endpoint rel. perm.]

0 (constant for reservoir) 3.4.44
1 (constant for each layer) 3.4.45-3.4.47
IEW [rel. perm. exponent]

IPERM = 2 or 3 ISRW [residual saturation]

(Parker and Lenhard for 0 (constant for reservoir) 3.4.37
water-wet or Lenhard for 1 (constant for each layer) 3.4.38-3.4.40
mixed-wet) 2 (each gridblock) 3.4.41-3.4.43
3.4.37 S1RWC, S2RWC, S3RWC (This line is read only if ISRW=0)
S1RWC - Residual saturation of aqueous phase displaced by oil at low capillary

number for entire reservoir.
Units: fraction
S2RWC - Residual saturation of oleic phase displaced by water at low capillary

Units: fraction
S3RWC - Residual saturation of microemulsion phase displaced by water at low

capillary number for entire reservoir.
Units: fraction
3.4.38 S1RW(K), for K=1, NZ (This line is read only if ISRW=1)
S1RW(K) - Residual saturation of aqueous phase displaced by oil or gas at low

capillary number for Kth layer.
Units: fraction
Note: S1RW(K) must begin a separate line in the input file for each layer. See
the note for input line 3.3.5 for information on how the data is read.
54
UTCHEM User's Guide
S2RW(K) - Residual saturation of oleic phase displaced by water at low capillary

number for Kth layer.
Units: fraction
S3RW(K) - Residual saturation of microemulsion phase displaced by water or oil

at low capillary number for Kth layer.
Units: fraction
3.4.41 S1RW(I), for I=1, NX × NY × NZ (This line is read only if ISRW=2)
S1RW(I) - Residual saturation of aqueous phase displaced by oil or gas at low

capillary number for Ith gridblock.
Units: fraction
Note: See the note for input line 3.3.6 for information on how this data should be
ordered.
S2RWC(K) - Residual saturation of oleic phase displaced by water at low

capillary number for Ith gridblock.
Units: fraction
ordered.
S3RW(I) - Residual saturation of microemulsion phase displaced by water or oil

at low capillary number for Ith gridblock.
Units: fraction
ordered.
55
UTCHEM User's Guide
Relative Permeability Endpoint and Exponent Data (input lines 3.4.44-

3.4.57) — These lines are required only if Corey function (IPERM<2) is used.
3.4.44 P1RWC, P2RWC, P3RWC (This line is read only if IPERM<2 and IPRW=0)
P1RWC - End point relative permeability of water at low capillary number for
entire reservoir.
P2RWC - End point relative permeability of oil at low capillary number for entire
reservoir.
P3RWC - End point relative permeability of microemulsion at low capillary

3.4.45 P1RW(K), for K=1, NZ (This line is read only if IPERM<2 and IPRW=1)
P1RW(K) - Constant endpoint relative permeability of water at low capillary

Note: See the note for input line 3.3.5 for information on how the data is read.
P2RW(K) - Constant endpoint relative permeability of oil at low capillary number

for Kth layer.
P3RW(K) - Constant endpoint relative permeability of microemulsion at low

3.4.48 P1RW(I), for I=1, NX × NY × NZ (This line is read only if IPERM<2 and
IPRW=2)
P1RW(I) - Endpoint relative permeability of water at low capillary number for Ith
gridblock.
ordered.
56
UTCHEM User's Guide
IPRW=2)
P2RW(I) - Endpoint relative permeability of oil at low capillary number for Ith
gridblock.
ordered.
IPRW=2)
P3RW(I) - Endpoint relative permeability of microemulsion at low capillary

number for Ith gridblock.
ordered.
3.4.51 E1WC, E2WC, E3WC (This line is read only if IPERM<2 and IEW=0)
E1WC - Phase relative permeability exponent for aqueous phase at low capillary
E2WC - Phase relative permeability exponent for oleic phase at low capillary
E3WC - Phase relative permeability exponent for microemulsion phase at low

capillary number system for entire reservoir.
3.4.52 E1W(K), for K=1, NZ (This line is read only if IPERM<2 and IEW=1)
E1W(K) - Relative permeability exponent of aqueous phase at low capillary

E2W(K) - Relative permeability exponent of oleic phase at low capillary number

for Kth layer.
57
UTCHEM User's Guide
E3W(K) - Relative permeability exponent of microemulsion phase at low

3.4.55 E1W(I), for I=1, NX × NY × NZ (This line is read only if IPERM<2 and
IEW=2)
E1W(I) - Relative permeability exponent of aqueous phase at low capillary

ordered.
IEW=2)
E2W(I) - Relative permeability exponent of oleic phase at low capillary number

for Ith gridblock.
ordered.
IEW=2)
E3W(I) - Relative permeability exponent of microemulsion phase at low capillary

ordered.
High Capillary Number Data (input lines 3.4.58-3.4.60) — These lines are
required only if ITRAP=1 or ITRAP=2.
3.4.58 S1RC, S2RC, S3RC (This line is read only if ITRAP=1 or 2)
S1RC - Residual saturation of aqueous phase at high capillary number.

Units: fraction
S2RC - Residual saturation of oleic phase at high capillary number.

Units: fraction
S3RC - Residual saturation of microemulsion phase at high capillary number.

Units: fraction
58
UTCHEM User's Guide
Note: The residual saturations at high capillary number can not be set equal to
those at low capillary number.
3.4.59 P1RC, P2RC, P3RC (This line is read only for ITRAP=1 or 2 and IPERM< 2)
P1RC - End point relative permeability of aqueous phase at high capillary number
condition.
P2RC - End point relative permeability of oleic phase at high capillary number
condition.
P3RC - End point relative permeability of microemulsion phase at high capillary

number condition.
3.4.60 E13C, E23C, E31C (This line is read only for ITRAP=1 or 2 and IPERM< 2)
E13C, E23C, E31C - Parameters used for calculating exponents for relative
permeability calculations at high capillary number.
Viscosity Data (input lines 3.4.61-3.4.66)
3.4.61 VIS1, VIS2, TSTAND
VIS1 - Water viscosity at reference temperature.

Units: cp = mPa.s
VIS2 - Oil viscosity at reference temperature.

Units: cp = mPa.s
TSTAND - Reference temperature

Units: ˚F (IUNIT=0) or ˚C (IUNIT=1)
3.4.62 IOVIS (This line is read only if NO>1)
IOVIS - Flag indicating whether the viscosity is a function of organic

composition.
Possible Values:
0 - Viscosity does not depend on the organic species
concentration.
1 - Viscosity depends on the organic species concentration.
information on the IOVIS=0 option and Section 7 of the UTCHEM
Technical Documentation for more information on the IOVIS=1
option.
59
UTCHEM User's Guide
3.4.64 OVIS(K), for K=1, NO (This line is read only if NO>1 and IOVIS=1)
OVIS(K) - Viscosity for organic component K at reference temperature.

Units: cp = mPa.s
3.4.64 VIS4, VSLOPG (This line is read only if IGAS≥1)
VIS4 - Gas viscosity at reference temperature and reference pressure.

Units: cp = mPa.s
VSLOPG - Slope of gas viscosity.

Units: (psi)-1 (IUNIT=0) or (kPa)-1 (IUNIT=1)
3.4.65 BVI(1), BVI(2) (This line is read only if IENG=1)
BVI(1) - Parameter for calculating water viscosity as a function of reservoir

temperature.
Units: (˚K)-1
BVI(2) - Parameter for calculating oil viscosity as a function of reservoir

temperature.
Units: (˚K)-1
3.4.66 BVI(4) (This line is read only if IGAS≥1 and IENG=1)
BVI(4) - Parameter for calculating gas viscosity as a function of reservoir

temperature.
Units: (˚K)-1
Corey Option Gas Relative Permeability Data (input lines 3.4.67-3.4.78)

— These lines are required only if IGAS≥1 and IPERM=0.
the gas relative permeability parameters for different options:
ISRW (residual saturation 0 (constant for reservoir) 3.4.67

flag) 1 (constant for each layer) 3.4.68-3.4.69
2 (variable over reservoir) 3.4.70-3.4.71
IPRW (endpoint relative 0 (constant for reservoir) 3.4.72

permeability flag) 1 (constant for each layer) 3.4.73
2 (variable over reservoir) 3.4.74
IEW (relative permeability 0 (constant for reservoir) 3.4.75

exponent flag) 1 (constant for each layer) 3.4.76
60
UTCHEM User's Guide
3.4.67 S2RWC4, S4RWC (This line is read only if ISRW=0)
S2RWC4 - Constant residual oil saturation to displacing gas phase for entire
reservoir.
Units: fraction
S4RWC - Constant residual gas saturation for entire reservoir.

Units: fraction
3.4.68 S2RW4(K), for K=1, NZ (This line is read only if ISRW=1)
S2RW4(K) - Constant residual oil saturation to displacing gas phase for Kth
layer.
Units: fraction
S4RW(K) - Constant residual gas saturation for Kth layer.

Units: fraction
3.4.70 S2RW4(I), for I=1, NX × NY × NZ (This line is read only if ISRW=2)
S2RW4(I) - Constant residual oil saturation to displacing gas phase for Ith
gridblock.
Units: fraction
ordered.
S4RW(I) - Residual gas saturation for Ith gridblock.

Units: fraction
ordered.
3.4.72 P4RWC (This line is read only if IPRW=0)
P4RWC - Constant gas endpoint relative permeability for entire reservoir.

3.4.73 P4RW(K), for K=1, NZ (This line is read only if IPRW=1)
P4RW(K) - Constant gas endpoint relative permeability for Kth layer.

61
UTCHEM User's Guide
3.4.74 P4RW(I), for I=1, NX × NY × NZ (This line is read only if IPRW=2)
P4RW(I) - Constant gas endpoint relative permeability for Ith gridblock.

ordered.
3.4.75 E4WC (This line is read only if IEW=0)
E4WC - Constant gas relative permeability exponent for entire reservoir.

3.4.76 E4W(K), for K=1, NZ (This line is read only if IEW=1)
E4W(K) - Constant gas relative permeability exponent for Kth layer.

3.4.77 E4W(I), for I=1, NX × NY × NZ (This line is read only if IEW=2)
E4W(I) - Constant gas relative permeability exponent for Ith gridblock.

ordered.
3.4.78 S4RC, P4RC, E4C, T44 (This line is read only if ITRAP=1)
S4RC - Residual gas saturations at high capillary number.

Units: fraction
P4RC - Gas endpoint relative permeability at high capillary number.

E4C - Gas relative permeability exponent at high capillary number.

T44 - Gas phase trapping parameter.

3.4.79 XIFTG, XIFTGW (This line is read only if IGAS≥1 or IPERM=2)
XIFTG - Log of interfacial tension between gas and oil.

Units: dyne/cm = mN/m
XIFTGW - Log of interfacial tension between gas and water.

Units: dyne/cm = mN/m
62
UTCHEM User's Guide
Microemulsion Viscosity Data (input line 3.4.80)
3.4.80 ALPHAV(I), for I=1, 5
ALPHAV(I) - Compositional phase viscosity parameters.

Note: All five viscosity parameters must be positive values.

information on the compositional viscosity model.
Polymer Property Data (input lines 3.4.81-3.4.84) — See Section 2 of the

UTCHEM Technical Documentation for information on the polymer property models.
3.4.81 AP1, AP2, AP3
AP1, AP2, AP3 - Parameters used for calculating polymer viscosity at zero shear
rate as a function of polymer and electrolyte concentrations.
Units: (wt. %)-1, (wt. %)-2, (wt. %)-3
3.4.82 BETAP, CSE1, SSLOPE
BETAP - Parameter for calculating the effective divalent salinity used to calculate
polymer viscosity.
CSE1 - Value below which the polymer viscosity is considered to be independent

of salinity.
Units: meq/ml
SSLOPE - Slope of viscosity vs. effective salinity on a log-log plot—assumed to

be constant.
Note: This value is usually large and negative for hydrolyzed
polyacrylamides and small and positive for polysaccharides.
3.4.83 GAMMAC, GAMHF, POWN
GAMMAC - Coefficient in shear rate equation below.
( )
12
day(darcy)
12 day µm 2
Units: (IUNIT=0) or (IUNIT=1)
ft − sec m − sec
GAMHF - Shear rate at which polymer viscosity is one half polymer viscosity at
zero shear rate.
Units: sec-1
POWN - Exponent for calculating shear rate dependence of polymer viscosity.
63
UTCHEM User's Guide
3.4.84 IPOLYM, EPHI3, EPHI4, BRK, CRK
IPOLYM - Flag indicating type of polymer partitioning.

Possible Values:
0 - All polymer exists in aqueous phase if aqueous phase exists;
otherwise, it exists completely in microemulsion phase
1 - Partitioning of polymer to water component is constant
EPHI3 - Effective porosity for surfactant—ratio of apparent porosity for

surfactant to actual porosity.
EPHI4 - Effective porosity for polymer—ratio of apparent porosity for polymer

to actual porosity.
BRK - Parameter for calculating permeability reduction factor.

volume of polymer − rich phase
Units:
weight% polymer
CRK - Parameter for calculating permeability reduction factor.

Units: (darcy)1/2 (100 g/g)-1/3 = (µm2)1/2 (100 g/g)-1/3)
Note: EPHI3 and EPHI4 are used to account for inaccessible pore volume in the
case of surfactant and polymer.
φsurfactant = φ × EPHI3
φpolymer = φ × EHPI4
The effect of permeability reduction or residual resistance is to reduce the

mobility of the polymer rich phase. This is accounted for by multiplying
the viscosity of the phase by BRK.
Component Density Data (input lines 3.4.85-3.4.92)
3.4.85 DEN1, DEN2, DEN23, DEN3, DEN7, DEN8, IDEN (This line is read only if
NO≤1)
DEN1 - Specific weight or density of water (Component 1).

Units: psi/ft (IUNIT=0) or g/cm3 (IUNIT=1)
DEN2 - Specific weight or density of oil (Component 2).
DEN23 - Coefficient of oil in microemulsion phase density calculations.
DEN3 - Specific weight or density of surfactant (Component 3).
64
UTCHEM User's Guide
DEN7 - Specific weight or density of Alcohol 1 (Component 7).

DEN8 - Specific weight or density of Alcohol 2 (when IGAS=0) or gas (when
IGAS≥1) (Component 7).
IDEN - Flag indicating if gravity effect should be considered.
Possible Values:
1 - Do not consider gravity effect
2 - Consider gravity effect
3.4.86 DEN1, DEN2, DEN3, DEN7, DEN8, IDEN, IODEN (This line is read only if
NO>1)
DEN1 - Specific weight or density of water (Component 1).

DEN2 - Specific weight or density of oil (Component 2).
DEN3 - Specific weight or density of surfactant (Component 3).
DEN7 - Specific weight or density of Alcohol 1 (Component 7).
DEN8 - Specific weight or density of Alcohol 2 (when IGAS=0) or gas (when
IDEN - Flag indicating if gravity effect should be considered.
Possible Values:
1 - Do not consider gravity effect
2 - Consider gravity effect
IODEN - Flag indicating if specific weight/density is a function of organic species

concentration.
Possible Values:
0 - Does not depend on organic species concentration
1 - Depends on organic species concentration
Note: See Section 2 of the UTCHEM Technical Documentation for
information on the IODEN=0 option or Section 7 of the UTCHEM
Technical Documentation for information on the IODEN=1 option.
Note: Specific weight for pure water is 0.433 psi/ft (density of 1 g/cm3 ).
IODEN must be set to 1 if any non-aqueous phase species (those with
indices < (8+NO)) participate in biodegradation equations.
65
UTCHEM User's Guide
Multiple Organic Density Data (input lines 3.4.87-3.4.88) — These lines are
required only if IODEN=1 and NO>1.
3.4.87 DNOILC(K), for K=1, NO (This line is read only if IODEN=1 and NO>1)
DNOILC(K) - Specific weight or density of organic component K for oleic

phase.
3.4.88 DNOME(K), for K=1, NO (This line is read only if IODEN=1 and NO>1)
DNOME(K) - Specific weight or density of organic component K for

microemulsion phase.
Multiple Organic Data (input lines 3.4.89-3.4.92) — These lines are required
only if NO>1. Furthermore, if (IFGHBN=1 and IOD=0) or IOVIS=1.
3.4.89 INAME (This line is read only if NO>1 and ((IFGHBN=1 and IOD=0) or
IOVIS=1)))
INAME - Flag indicating whether name of the organic components will be

provided by user.
Possible Values:
0 - Name of the organic components will be provided; the
molecular weight and equivalent alkane carbon number will
be obtained from the built-in database
1 - Molecular weight and equivalent alkane carbon number for
each organic components will be provided.
3.4.90 ONAME(K), for K=1, NO (This line is read only if NO>1 and ((IFGHBN=1 and
IOD=0) or IOVIS=1)) and INAME=0)
ONAME(K) - Name of organic component K. See input line 3.4.14 for a list of
valid component names.
3.4.91 OMWT(K), for K=1, NO (This line is read only if NO>1 and ((IFGHBN=1 and
OMWT(K) - Molecular weight for organic component K.
3.4.92 OEACN(K), for K=1, NO (This line is read only if NO>1 and ((IFGHBN=1 and
OEACN(K) - equivalent alkane carbon number for organic component K.
Note: See Section 7 of the UTCHEM Technical Documentation for information

on the equivalent alkane carbon number.
66
UTCHEM User's Guide
3.4.93 ISTB
ISTB - Flag indicating the units to be used when printing injection and production
rates.
Possible Values:
0 - Rates printed at bottomhole condition in ft3 or m3
1 - Rates printed at surface condition in bbls
3.4.94 FVF(L), for L=1, MXP (This line is read only if ISTB=1 and IUNIT=0)
FVF(L) - Formation volume factor for Lth phase.

Units: SCF/ft3
Note: MXP=3 when IGAS=0 and MXP=4 when IGAS≥1.
Fluid Compressibility Data (input lines 3.4.95-3.4.97)
3.4.95 COMPC(1), COMPC(2), COMPC(3), COMPC(7), COMPC(8)
COMPC(1) - Compressibility of brine (Component 1).

COMPC(2) - Compressibility of oil (Component 2).

COMPC(3) - Compressibility of surfactant (Component 3).

COMPC(7) - Compressibility of Alcohol 1 (Component 7).

COMPC(8) - Compressibility of Alcohol 2 (when IGAS=0) or gas (when

Note: For incompressible fluids, values of zero should be used for the COMPC
values listed above.
3.4.96 ICOMPO (This line is read only if NO>1)
ICOMPO - Flag indicating whether each organic component has different

compressibility.
Possible Values:
0 - All organic components have the same compressibility as
COMPC(2)
1 - Each organic component has different compressibility
3.4.97 COMPC(KC), for KC=1, NO (This line is read only if NO>1 and ICOMPO=1)
COMPC(KC) - Compressibility of organic component KC.

67
UTCHEM User's Guide
Capillary Pressure Data (input lines 3.4.98-3.4.116) — See Sections 2, 3, and

14 of the UTCHEM Technical Documentation for capillary pressure information.
3.4.98 ENTP, CLAMDA, ALPDS, ALPWS, VNDS, VNWS, PNEG (This line is read
only if IPERM=3)
ENTP - Brooks-Corey entry pressure.

CLAMDA - Exponent of Brooks-Corey capillary pressure function.
ALPDS - Alpha parameter in the van Genuchten capillary pressure function for
drying path (water saturation is decreasing).
ALPWS - Alpha parameter in the van Genuchten capillary pressure function for
wetting path (water saturation is increasing).
VNDS - Exponent n in the van Genuchten capillary pressure function for drying
path.
VNWS - Exponent n in the van Genuchten capillary pressure function for wetting
path.
PNEG - Maximum negative capillary pressure.

Note: The value for PNEG must be a negative number.
3.4.99 ICPC, IEPC, IOW (This line is read only if IPERM≤2)
ICPC - Flag indicating type of capillary pressure endpoint.

Possible Values:
0 - Constant capillary pressure endpoint for entire reservoir
1 - Constant capillary pressure endpoint for each layer
2 - Capillary pressure endpoint specified for each gridblock
IEPC - Flag indicating type of capillary pressure exponent.

Possible Values:
0 - Constant capillary pressure exponent for entire reservoir
1 - Constant capillary pressure exponent for each layer
2 - Capillary pressure exponent specified for each gridblock
IOW - Flag indicating the wettability for two-phase oil/water capillary pressure
calculations using imbibition Corey function (IPERM=0).
Possible Values:
0 - The capillary pressure curve is for strongly water-wet rock
(default)
1 - The capillary pressure curve is for strongly oil-wet rock
2 - The capillary pressure curve is for mixed-wet rocks
Note: IOW=1 and 2 are only available with the imbibition Corey function
(IPERM=0).
68
UTCHEM User's Guide
the capillary pressure data for different values of IPERM:
If IPERM = 0
IOW=0 or 1 (strongly ICPC [capillary pressure endpoint]

water- or oil-wet) 0 (constant for reservoir) 3.4.100
1 (constant for each layer) 3.4.101
IEPC [capillary pressure exponent]

IOW = 2 (mixed-wet) ICPC [capillary pressure endpoint]

2 (variable over reservoir) 3.4.112-3.4.116
If IPERM = 1 or 2
IOW=0 (strongly ICPC [capillary pressure endpoint]

water-wet) 0 (constant for reservoir) 3.4.100
IEPC [capillary pressure exponent]

If IPERM=3 go to input line 3.4.117
Capillary Pressure Data for Strongly Water- or Oil-Wet Rocks (input lines
3.4.100-3.4.105) — These lines are required only if IPERM≤2 and IOW<2.
3.4.100 CPC0 (This line is read only if IPERM≤2 and IOW<2 and ICPC=0)
CPC0 - Capillary pressure endpoint for entire reservoir (if IPERM<2) or van
Genuchten capillary pressure parameter, α, for entire reservoir (if
IPERM=2).
69
UTCHEM User's Guide
Units: if IUNIT=0 if IUNIT=1

if IPERM<2: psi darcies kPa µm 2
if IPERM=2: ( psi darcies )−1 µm 2 / kPa
3.4.101 CPC(K), for K=1, NZ (This line is read only if IPERM≤2 and IOW<2 and
ICPC=1)
CPC(K) - Capillary pressure endpoint for Kth layer (if IPERM<2) or van
Genuchten capillary pressure parameter, α, for Kth layer (if IPERM=2).
3.4.102 CPC(I), for I=1, NX × NY × NZ (This line is read only if IPERM≤2 and
IOW<2 and ICPC=2)
CPC(I) - Capillary pressure endpoint for Ith gridblock (if IPERM<2) or van
Genuchten capillary pressure parameter, α , for Ith gridblock (if
IPERM=2).
3.4.103 EPC0 (This line is read only if IPERM≤2 and IOW<2 and IEPC=0)
EPC0 - Capillary pressure exponent for entire reservoir (if IPERM<2) or van
Genuchten capillary pressure parameter, n, for entire reservoir (if
IPERM=2).
3.4.104 EPC(K), for K=1, NZ (This line is read only if IPERM≤2 and IOW<2 and
IEPC=1)
EPC(K) - Capillary pressure exponent for Kth layer (if IPERM<2) or van
Genuchten capillary pressure parameter, n, for Kth layer (if IPERM=2).
3.4.105 EPC(I), for I=1, NX × NY × NZ (This line is read only if IPERM≤2 and
IOW<2 and IEPC=2)
EPC(I) - Capillary pressure exponent for Ith gridblock (if IPERM<2) or van
Genuchten capillary pressure parameter, n, for Ith gridblock (if
IPERM=2).
70
UTCHEM User's Guide
Capillary Pressure Data for Mixed-Wet Rocks using Corey's option (input
lines 3.4.106-3.4.116) — These lines are required only if I P E R M = 0 and
IOW=2.
3.4.106 CPCW, EPCW, CPCO, EPCO, SSTAR (This line is read only if IPERM=0 and
IOW=2 and ICPC=0)
CPCW - Capillary pressure endpoint for entire reservoir for positive branch.
Units: psi darcies (IUNIT=0) or kPa µm 2 (IUNIT=1)
EPCW - Capillary pressure exponent for entire reservoir for positive branch.
CPCO - Capillary pressure endpoint for entire reservoir for negative branch.
EPCO - Capillary pressure exponent for entire reservoir for negative branch.
SSTAR - Water saturation where the capillary pressure is zero.

3.4.107 CPC(K,1), for K=1, NZ (This line is read only if IPERM=0 and IOW=2 and
ICPC=1)
CPC(K,1) - Capillary pressure endpoint for Kth layer for positive branch.
3.4.108 EPC(K,1), for K=1, NZ (This line is read only if IPERM=0 and IOW=2 and
ICPC=1)
EPC(K,1) - Capillary pressure exponent for Kth layer for positive branch.
3.4.109 CPC(K,2), for K=1, NZ (This line is read only if IPERM=0 and IOW=2 and
ICPC=1)
CPC(K,2) - Capillary pressure endpoint for Kth layer for negative branch.
71
UTCHEM User's Guide
3.4.110 EPC(K,2), for K=1, NZ (This line is read only if IPERM=0 and IOW=2 and
ICPC=1)
EPC(K,2) - Capillary pressure exponent for Kth layer for negative branch.
3.4.111 SSTAR(K), for K=1, NZ (This line is read only if IPERM=0 and IOW=2 and
ICPC=1)
SSTAR(K) - Water saturation for Kth layer where the capillary pressure is zero.
3.4.112 CPC(I,1), for I=1, NX × NY × NZ (This line is read only if IPERM=0 and
IOW=2 and ICPC=2)
CPC(I,1) - Capillary pressure endpoint for Ith gridblock for positive branch.
Note: See the note for input line 3.3.6 for information on how this data should
be ordered.
3.4.113 EPC(I,1), for I=1, NX × NY × NZ (This line is read only if IPERM=0 and
IOW=2 and ICPC=2)
EPC(I,1) - Capillary pressure exponent for Ith gridblock for positive branch.
be ordered.
3.4.114 CPC(I,2), for I=1, NX × NY × NZ (This line is read only if IPERM=0 and
IOW=2 and ICPC=2)
CPC(I,2) - Capillary pressure endpoint for Ith gridblock for negative branch.
be ordered.
3.4.115 EPC(I,2), for I=1, NX × NY × NZ (This line is read only if IPERM=0 and
IOW=2 and ICPC=2)
EPC(I,2) - Capillary pressure exponent for Ith gridblock for negative branch.
72
UTCHEM User's Guide
be ordered.
3.4.116 SSTAR(I), for I= 1, NX × NY × NZ (This line is read only if IPERM=0 and

IOW=2 and ICPC=2)
SSTAR(I) - Water saturation in Ith gridblock where the capillary pressure is

zero.
be ordered.
Diffusion and Dispersion Data (input lines 3.4.117-3.4.124)
3.4.117 D(KC,1), for KC=1, N
D(KC,1) - Molecular diffusion coefficient of KCth component in aqueous

phase.
Units: ft2/day (IUNIT=0) or m2/day (IUNIT=1)
Note: The input diffusion coefficient should be divided by tortuosity
(D/τ) where the value of tortuosity is greater than one.
3.4.118 D(KC,2), for KC=1, N
D(KC,2) - Molecular diffusion coefficient of KCth component in oleic phase.

3.4.119 D(KC,3), for KC=1, N
D(KC,3) - Molecular diffusion coefficient of KCth component in microemulsion

phase.
3.4.120 D(KC,4), for KC=1, N (This line is read only if IGAS≥1)
D(KC,4) - Molecular diffusion coefficient of KCth component in gas phase.

3.4.121 ALPHAL(1), ALPHAT(1)
ALPHAL(1) - Longitudinal dispersivity of aqueous phase.

73
UTCHEM User's Guide
ALPHAT(1) - Transverse dispersivity of aqueous phase.

ALPHAL(2) - Longitudinal dispersivity of oleic phase.

ALPHAT(2) - Transverse dispersivity of oleic phase.

ALPHAL(3) - Longitudinal dispersivity of microemulsion phase.

ALPHAT(3) - Transverse dispersivity of microemulsion phase.

3.4.124 ALPHAL(4), ALPHAT(4) (This line is read only if IGAS≥1)
ALPHAL(4) - Longitudinal dispersivity of gas phase.

ALPHAT(4) - Transverse dispersivity of gas phase.

Adsorption Data (input lines 3.4.125-3.4.130)
3.4.125 IADSO
IADSO - Organic adsorption calculation flag.

Possible Values
0 - Organic adsorption is not considered
1 - Organic adsorption calculated using linear adsorption model
2 - Organic adsorption calculated using Langmuir isotherm
3 - Organic adsorption calculated using Freundlich isotherm

on the organic adsorption isotherm models.
3.4.126 IREV, DENS, IKADS (This line is read only if IADSO>0)
IREV - Flag indicating if organic adsorption is reversible or not.

Possible Values:
0 - Irreversible organic adsorption
1 - Reversible organic adsorption
DENS - Grain soil density

Units: lb/ft3 (IUNIT=0), g/cc (IUNIT=1)
74
UTCHEM User's Guide
IKADS - Flag indicating constant or variable organic adsorption parameter for

reservoir.
Possible Values:
0 - Constant organic adsorption parameter for entire reservoir
1 - Constant organic adsorption parameter for each layer
2 - Organic adsorption parameter specified for each gridblock
3.4.127 AD21(1,KO), AD22(1,KO), for KO=1, NO (This line is read only if IADSO>0
and IKADS=0)
AD21(1,KO) - Adsorption parameter for KOth organic component.
AD22(1,KO) - Adsorption parameter for KOth organic component.

Note: This line is read at least once even if NO=0.
All elements of the AD21(n,KO) array will be set equal to AD21(1,KO)

and all elements of the AD22(n,KO) array will be set equal to
AD22(1,KO).
Example: If there are two organic components (NO=2) with adsorption

parameters of 0.6 & 1.0 for the first component and 0.3 & 1.0 for the
second component, this line will be repeated twice (once for each
component) as follows:
0.6 1.0
0.3 1.0
3.4.128 AD21(K,KO), AD22(K,KO), for K=1, NZ, for KO=1, NO (This line is read
only if IADSO>0 and IKADS=1)
AD21(K,KO) - Adsorption parameter for KOth organic component in the Kth

layer.
AD22(K,KO) - Adsorption parameter for KOth organic component in the Kth

layer.
Note: This line is read at least once even if NO=0. See the note for input line
3.3.5 for information on how the data is read.
The user should provide NO sets (1 per organic component) of NZ lines

of data (1 per layer in the Z direction) of adsorption parameter pairs. Each
of these adsorption parameter sets is preceded by 3 comment lines.
Example: If there are two organic components (NO=2), and three layers (NZ=3),
the input data would appear as follows:
CC
CC
*----AD21 AD22 (1st organic component)
0.6 1.0
0.6 1.0
0.6 1.0
CC
75
UTCHEM User's Guide
CC
*----AD21 AD22 (2nd organic component)
0.3 1.0
0.3 1.0
0.3 1.0
3.4.129 The data on input lines 3.4.129.a and 3.4.129.b are repeated for KO=1 to NO
times. See note following input line 3.4.129.b for more information.
3.4.129.a AD21(I,KO), for I=1, NX × NY × NZ (This line is read only if IADSO>0 and
IKADS=2)
AD21(I,KO) - Adsorption parameter for KOth organic component in Ith

gridblock.
3.4.129.b AD22(I,KO), for I=1, NX × NY × NZ (This line is read only if IADSO>0 and
IKADS=2)
AD22(I,KO) - Adsorption parameter for KOth organic component in Ith

gridblock.
Note: This line is read at least once even if NO=0. See the note for input line
3.3.6 for information on how this data should be ordered.
The user should provide NO sets (1 per organic component) of

NX × NY × NZ adsorption parameters (1 per gridblock). Each set
consists of a separate input line for the AD21 and AD22 parameters.
Example: If there are two organic components (NO=2), NX=2, NY=2, and
NZ=1, the input data would appear as follows:
CC
CC
*----AD21 (1st organic component)
0.6 0.6 0.6 0.6
CC
CC
*----AD22 (1st organic component)
1.0 1.0 1.0 1.0
CC
CC
*----AD21 (2nd organic component)
4*0.3
CC
CC
*----AD22 (2nd organic component)
4*1.0
Note: This example also shows an input convention allowing the user
to specify multiple occurances of the same value. Using
4*0.3
is equivalent to using
76
UTCHEM User's Guide
0.3 0.3 0.3 0.3
3.4.130 AD31, AD32, B3D, AD41, AD42, B4D, IADK, IADS1, FADS, REFK
AD31 - Surfactant adsorption parameter.

AD32 - Surfactant adsorption parameter.

Units: ml/meq
B3D - Surfactant adsorption parameter.

volume of water
Units:
volume of surfactant
AD41 - Polymer adsorption parameter.

AD42 - Polymer adsorption parameter.

Units: ml/meq
B4D - Polymer adsorption parameter.

volume of water
Units:
weight% polymer
IADK - Permeability dependent adsorption flag

Possible Values:
0 - Surfactant and polymer adsorption is independent of
permeability
1 - Surfactant and polymer adsorption is dependent on
permeability
IADS1 - Flag to specify the competitive surfactant adsorption in the presence of

polymer
Possible Values:
0 - Surfactant and polymer adsorption are independent
1 - Competitive surfactant and polymer adsorption is considered
FADS - Parameter to adjust the competitive adsorption calculation
REFK - The reference permeability at which the surfactant and polymer

adsorption input parameters are measured.
Units: md

on the surfactant and polymer adsorption model.
3.4.131 QV, XKC, XKS, EQW
QV - Cation exchange capacity of clays.

Units: meq/ml of pore volume
77
UTCHEM User's Guide
XKC - Cation exchange constant for clays.

Units: (meq/ml)-1
XKS - Cation exchange constant for surfactant.
Units: (meq/ml)-1
EQW - Equivalent weight of surfactant.

on the cation exchange model.
Tracer Data (input lines 3.4.132-3.4.144) — These lines are required only if
NTW+NTA>0. See Section 4 of the UTCHEM Technical Documentation for more
details on tracer modeling in UTCHEM.
3.4.132 TK(I), for I=1, NTW+NTA (This line is read only if NTW+NTA>0)
TK(I) - Tracer partition coefficient for Ith water/oil tracer at initial chloride
concentration and reference temperature. A value of 0.0 indicates a water
or gas nonpartitioning tracer and a value of -1.0 indicates a
nonpartitioning oil tracer.
Units: fraction
3.4.133 TKS(I), for I=1, NTW, C5INI (This line is read only if NTW>0)
TKS(I) - Parameter for calculating water/oil tracer partition coefficient for Ith
tracer as a function of salinity.
Units: (meq/ml)-1
C5INI - Reference salinity at which the partition coefficients are measured.
Units: meq/ml water
3.4.134 TKT(I), for I=1, NTW+NTA (This line is read only if NTW+NTA>0 and
IENG=1)
TKT(I) - Parameter for calculating tracer partitioning coefficient for Ith tracer as
a function of reservoir temperature.
Units: (˚F)-1 (IUNIT=0) or (˚C)-1 (IUNIT=1)
3.4.135 RDC(I), for I=1, NTW+NTA (This line is read only if NTW+NTA>0)
RDC(I) - Radioactive decay coefficient for Ith tracer. A value of 0.0 indicates a
non–radioactive tracer.
Units: 1/days
3.4.136 RET(I), for I=1, NTW+NTA (This line is read only if NTW+NTA>0)
RET(I) - Tracer adsorption parameter (adsorbed concentration/flowing

concentration). A value of 0.0 indicates no retardation.
78
UTCHEM User's Guide
Dead-end Pore Model Data (input lines 3.4.137-3.4.140) — These lines are
required only if NTW+NTA>0 and ICAP=1.
3.4.137 FFL(1), FFH(1), CM(I,1), for I=1, NTW+NTA (This line is read only if
NTW+NTA>0 and ICAP=1)
FFL(1) - Water phase flowing fraction at fractional flow = 0.0.

FFH(1) - Water phase flowing fraction at fractional flow = 1.0.

CM(I,1) - Mass transfer coefficients for Ith tracer in water phase.

Units: 1/sec
FFL(2) - Oil phase flowing fraction at fractional flow = 0.0.

FFH(2) - Oil phase flowing fraction at fractional flow = 1.0.

CM(I,2) - Mass transfer coefficients for Ith tracer in oil phase.

Units: 1/sec
FFL(3) - Microemulsion phase flowing fraction at fractional flow = 0.0.

FFH(3) - Microemulsion phase flowing fraction for at fractional flow = 1.0.

CM(I,3) - Mass transfer coefficients for Ith tracer in microemulsion phase.

Units: 1/sec
NTW+NTA>0 and ICAP=1 and IGAS≥1)
FFL(4) - Value of flowing fraction for gas phase when fractional flow = 0.0.
FFH(4) - Value of flowing fraction for gas phase when fractional flow = 1.0.
CM(I,4) - Mass transfer coefficients for Ith tracer in gas phase.

Units: 1/sec
79
UTCHEM User's Guide
Reacting Tracer Data (input lines 3.4.141-3.4.144) — The following lines are
required only if NTW>0 and ITREAC=1.
3.4.141 NRT, TAK(I), for I=1, NRT (This line is read only if NTW>0 and ITREAC=1)
NRT - Number of reacting tracers.

Possible Values: 1 or 2
TAK(I) - Rate constant for a first-order aqueous phase reaction at reference

temperature for reacting tracer I
Units: days-1
Note: First reacting tracer is tracer 2 hydrolyzes to form tracer 3. The second
reacting tracer if present is tracer 4 hydrolyzes to form tracer 5.
3.4.142 TMW(I), for I=1, NTW (This line is read only if NTW>0 and ITREAC=1)
TMW(I) - Molecular weight of the Ith tracer.

Units: The user can specify the molecular weight in any unit as long as
the units are the same for all the tracers. It is assumed that the
reaction of 1 mole of primary tracer produces 1 mole of secondary
tracer. If not, use "equivalent" molecular weights.
3.4.143 TDEN(I), for I=1, NTW (This line is read only if NTW>0 and ITREAC=1)
TDEN(I) - Density of the Ith tracer.

Units: g/cm3
3.4.144 TAKT(I), for I=1, NRT (This line is read only if NTW>0 and ITREAC=1 and
IENG=1)
TAKT(I) - Parameter for calculating rate constant for a first-order aqueous phase
reaction as a function of reservoir temperature for reacting tracer I.
Units: (˚K)-1
Dual Porosity Data (input lines 3.4.145-3.4.163) — This section is required

for dual porosity option (ICAP=2) only. This option works only with English unit
(IUNIT=0) and the Cartesian coordinate (ICOORD=1). See Section 5 of the UTCHEM
Technical Documentation for more details on this option. The capability of dual porosity
option at this time is limited to single phase water and one tracer component. Restart
capability (IMODE=2) is not currently available for the dual porosity option.
3.4.145 NSUB, MSUB, ISUB (This line is read only if ICAP=2)
NSUB - Number of subgrids in lateral direction.
MSUB - Number of subgrids in vertical direction.
ISUB - Mode of subgridding

Possible Values:
0 - Uniform matrix block dimension
80
UTCHEM User's Guide
1 - Variable matrix block size dimension in each direction

2 - Variable matrix block dimension in whole reservoir
3.4.146 DXL(1), DYL(1), DZL(1) (This line is read only if ICAP=2 and ISUB=0)
DXL(1), DYL(1), DZL(1) - Uniform matrix block sizes in x y, and z directions.

Units: feet
Note: All values of the DXL, DYL, and DZL arrays are set to DXL(1), DYL(1),
and DZL(1), respectively.
3.4.147 II1, II2, XL1 (This line is read only if ICAP=2 and ISUB=1)
II1, II2 - First and last index for gridblocks with same size in X direction.
XL1− Matrix block size in X direction.

Units: feet
3.4.148 JJ1, JJ2, YL1 (This line is read only if ICAP=2 and ISUB=1)
JJ1, JJ2 - First and last index for gridblocks with same size in Y direction.
YL1 - Matrix block size in Y direction.

Units: feet
3.4.149 KK1, KK2, ZL1 (This line is read only if ICAP=2 and ISUB=1)
KK1, KK2 - First and last index for gridblocks with same size in Z direction
ZL1 - Matrix block size in Z direction.
Units: feet
3.4.150 DXL(I), for I=1, NX × NY × NZ (This line is read only if ICAP=2 and
ISUB=2)
DXL(I) - Matrix block size in X direction for Ith reservoir gridblock.

Units: feet
3.4.151 DYL(I), for I=1, NX × NY × NZ (This line is read only if ICAP=2 and
ISUB=2)
DYL(I) - Matrix block size in Y direction for Ith reservoir gridblock.

Units: feet
3.4.152 DZL(I), for I=1, NX × NY × NZ (This line is read only if ICAP=2 and
ISUB=2)
DZL(I), - Matrix block size in Z direction for Ith reservoir gridblock.

Units: feet
81
UTCHEM User's Guide
3.4.153 VFRACM(I), for I=1, NSUB (This line is read only if ICAP=2)
VFRACM(I) - Volume fraction of Ith subgrid. These values are used to generate
subgrids in the lateral direction.
3.4.154 ZLSUB(1,I), for I=1, MSUB (This line is read only if ICAP=2)
ZLSUB(1,I) - Thickness fraction of Ith vertical subgrid. These values are used
to generate subgrids in the vertical direction.
3.4.155 KPH, KKX, KKZ (This line is read only if ICAP=2)
KPH - Flag for matrix porosity distribution.

Possible Values:
0 - Uniform matrix porosity
1 - Variable matrix porosity
KKX - Flag for diffusion coefficient distribution in lateral direction in matrix.

Possible Values:
0 - Uniform matrix diffusion coefficient
1 - Uniform matrix diffusion coefficient in each reservoir layer
2 - Variable matrix diffusion coefficient at each reservoir node
KKZ - Flag for diffusion coefficient distribution in vertical direction in matrix.

Possible Values:
0 - Uniform matrix diffusion coefficient
1 - Uniform matrix diffusion coefficient in each reservoir layer
2 - Variable matrix diffusion coefficient at each reservoir node
3.4.156 PHIC (This line is read only if ICAP=2 and KPH=0)
PHIC - Matrix porosity.
3.4.157 PORCM(I,1,1), for I=1, NX × NY × NZ (This line is read only if ICAP=2 and
KPH=1)
PORCM(I,1,1) - Matrix porosity of the Ith reservoir node.
3.4.158 DMC (This line is read only if ICAP=2 and KKX=0)
DMC - Matrix diffusion coefficient in lateral direction.

Units: ft2/day
3.4.159 DMX(K), for K=1, NZ (This line is read only if ICAP=2 and KKX=1)
DMX(K) - Matrix diffusion coefficient in Kth reservoir layer in lateral direction.

Units: ft2/day
3.4.160 DMX(I), for I=1, NX × NY × NZ (This line is read only if ICAP=2 and
KKX=2)
DMX(I) - Matrix diffusion coefficient at Ith reservoir node in lateral direction.
82
UTCHEM User's Guide
Units: ft2/day
3.4.161 DMC (This line is read only if ICAP=2 and KKZ=0)
DMC - Matrix diffusion coefficient in vertical direction.

Units: ft2/day
3.4.162 DMZ(K), for K=1, NZ (This line is read only if ICAP=2 and KKZ=1)
DMZ(K) - Matrix diffusion coefficient in Kth reservoir layer in vertical direction.

Units: ft2/day
3.4.163 DMZ(I), for I=1, NX × NY × NZ (This line is read only if ICAP=2 and
KKZ=2)
DMZ(I) - Matrix diffusion coefficient at Ith reservoir node in vertical direction.

Units: ft2/day
Gel Reaction Data (input lines 3.4.164-3.4.168) — These lines are required
only if IREACT=1 or IREACT=4 and NG>0. Refer to Section 6 of the UTCHEM
Technical Documentation for more details on this option.
3.4.164 KGOPT, AK1, AK2, SCR, X4, X13, X14, X16, WM4 (This line is read only
if IREACT=1 or 4 and NG>0)
KGOPT - Flag to specify the gelation type used.

Possible Values:
1 - Polymer/chromium chloride gel
2 - Polymer/chromium malonate gel
3 - Silicate gel
AK1 - Kinetic rate coefficient for NG1 and NG2 at reference temperature
(KGOPT=1).
Units: ppm-1 days-1
AK2 - Kinetic rate coefficient for gel at reference temperature .

Units: (mole/liter)1-X4-X14+X16 days-1 for KGOPT=1
(mole/liter)1-X4-X13+X16 days-1 for KGOPT=2
(mole/liter)1-X4+X14 days-1 for KGOPT=3
SCR - Stoichiometric ratio in mass between Cr3+ and polymer.

X4 - Exponent to be used for polymer or silicate in gelation reaction.

X13 - Exponent to be used for component (NG2) in gelation reaction for

KGOPT=2.
83
UTCHEM User's Guide
X14 - Exponent to be used for the third gel option component (NG3) in gelation
reaction.
X16 - Exponent to be used for hydrogen ion component (NG5) of gelation

reaction.
WM4 - Molecular weight of polymer (KGOPT=1 or 2) or silicate (KGOPT=3).

Units: g/mole
Note: The following tables define different gelation options and corresponding
gel components.
For IREACT=1:
Component No. KGOPT=1 KGOPT=2 KGOPT=3
4 Polymer Polymer Silicate
NG1 Na2Cr2O7 - -
NG2 CSN2H4 Malonate ion -
NG3 Cr(III) Cr(III) OH-
NG4 Gel Gel Gel
NG5 Hydrogen Hydrogen -
For IREACT=4:
Component NO. KGOPT=1 KGOPT=2
4 Polymer Polymer
NGC1 Sodium Sodium
NGC2 Hydrogen Hydrogen
NGC3 Magnesium Magnesium
NGC4 Carbonate Carbonate
NGC5 Chromium* Chromium*
NGC6 Silica Silica
NG1 Na2Cr2O7 -
NG2 CSN2H4 Malonate ion
NG3 Cr(III)** Cr(III)**
NG4 Gel Gel
Where NG and NGC are the gel option and geochemistry option species.
3.4.165 AK1T, AK2T (This line is read only if IREACT=1 or 4 and NG>0 and
IENG=1)
AK1T - Parameter for calculating Kinetic rate coefficient for Cr3+ as a function
of reservoir temperature.
Units: (˚K)-1
AK2T - Parameter for calculating Kinetic rate coefficient for gel as a function of
reservoir temperature.
Units: (˚K)-1
84
UTCHEM User's Guide
3.4.166 AG1, AG2, CRG, AGK, BGK (This line is read only if IREACT=1 or 4 and
NG>0)
AG1 - Flory-Huggins parameter for gel viscosity.

Units: cp ppm-1 = m = Pa.s ppm-1
AG2 - Flory-Huggins parameter for gel viscosity.
Units: cp ppm-2 = mPa.s ppm-2
CRG - Constant in the dimensionless pore radius reduction group. This
constant depends on the gel type.
Units: darcy (wt%)1 3 = µm 2 (wt%)1 3
AGK, BGK - Permeability reduction parameters for Langmuir correlation with

gel concentration.
3.4.167 A15D, B15D, ICREX, A14D, B14D, CRNAK, HNAK, C160 (This line is read
only if IREACT=1 or 4 and NG>0)
A15D, B15D - Gel adsorption parameters.

Units: vol. of water / ppm gel
ICREX - Flag indicating if Cr3+ will be allowed to exchange with clays.

Possible Values:
0 - Cr3+ exchange with clays is not allowed
1 - Cr3+ exchange with clays is allowed
A14D, B14D - Chromium adsorption parameters.
Units: vol. of water / ppm chromium
CRNAK - Chromium-sodium exchange reaction equilibrium constant.
HNAK - Hydrogen-sodium exchange reaction equilibrium constant.
C160 - Initial hydrogen ion concentration.

Units: meq/ml
Note: The input values of CRNAK, HNAK, and C160 are ignored for
IREACT=4
3.4.168 IP1, IP2 (This line is read only if IREACT=1 or 4 and NG>0 and NY=1 and
NZ=1)
IP1, IP2 - Gridblock locations where calculated pressure values should be

printed to UNIT 19.
Note: These values are intended to be used for comparison with pressure tab
data of 1-D experiments.
85
UTCHEM User's Guide
Temperature Data (input lines 3.4.169-3.4.172) — These lines are required

only if temperature variation is considered in the simulation for IENG=1.
3.4.169 TEMPI (This line is read only if IENG=1)
TEMPI- Constant initial reservoir temperature.

3.4.170 DENS, CRTC, CVSPR, (CVSPL(L), for L=1, MXP) (This line is read only if
IENG=1)
DENS - Reservoir rock density.

Units: lb/ft3 (IUNIT=0) or g/cm3 (IUNIT=1)
CRTC - Reservoir thermal conductivity.
Units: Btu (day-ft-˚F)-1 (IUNIT=0) or kJ (day-m-˚K)-1 (IUNIT=1)
CVSPR - Reservoir rock heat capacity.

Units: Btu (lb-˚F)-1 (IUNIT=0) or kJ (kg-˚K)-1 (IUNIT=1)
CVSPL(L) - Phase L heat capacity (MXP is equal to 3 (IGAS=0) or 4

(IGAS≥1)).
3.4.171 IHLOS, IANAL (This line is read only if IENG=1)
IHLOS - Flag indicating if the heatloss calculation to overburden and

underburden rock is considered or not.
Possible Values:
0 - Heatloss is not considered
1 - Heatloss is considered
IANAL - Flag indicating if the temperature profile is calculated from analytical

solution (only 1-D).
Possible Values:
0 - Analytical solution is not considered
1 - Analytical solution is considered
3.4.172 TCONO, DENO, CVSPO, TCONU, DENU, CVSPU (This line is read only if
IHLOS=1 and if IENG=1)
TCONO - Thermal conductivity of overburden rock.

DENO - Density of overburden rock.

CVSPO - Heat capacity of overburden rock.

TCONU - Thermal conductivity of underburden rock.
86
UTCHEM User's Guide
DENU - Density of underburden rock.

CVSPU - Heat capacity of underburden rock.
Foam Model Data (input lines 3.4.173 and 3.4.174) — These lines are
required only if the foam option is considered (IGAS=2).
3.4.173 RFMAX, SOSTAR, CSTAR, EPXLO, SHRTN, VELGR (This line is read only
if IGAS=2)
RFMAX- Maximum foam "R" parameter.

SOSTAR - Critical oil saturation above which foam is not generated.

CSTAR - Critical surfactant concentration below which foam is not generated.

EPXLO - Water saturation tolerance parameter in foam model.

SHRTN - Gas shear thinning exponent.

VELGR - Reference gas velocity.

Units: ft/day (IUNIT=0) or m/day (IUNIT=1)
3.4.174 SWSTAR(I), for I=1, NX × NY × NZ (This line is read only if IGAS=2)
SWSTAR(I)- Water saturation at critical capillary pressure for Ith gridblock.

Note: SWSTAR(I) is assumed to be corrected for the permeability used in the

simulation.
87
UTCHEM User's Guide
User notes:
88
UTCHEM User's Guide
3.5 Physical Property Data for Geochemical Options
The fifth input section consists of physical property data for geochemistry
option and it is required only if IREACT>1. The data for this section is
generated by a preprocessor program (EQBATCH) and does not have the same
format as the rest of the input data for UTCHEM. This input section is read by a
separate routine called GEOREAD not preceded by the usual seven comment lines
and individual data lines are not preceded by three comment lines. Appendix B
gives an example for the list of elements, fluid species, solid species, and adsorbed
species for geochemical options. See Section 8 of the UTCHEM Technical
Documentation for information on the geochemical model.
3.5.1 IRSPS, IPHAD
IRSPS - Flag indicating if the reactive species concentrations should be printed.

Possible Values:
0 - Reactive species concentrations will not be printed
1 - Independent aqueous reactive species, solid species, and sorbed
species concentrations will be printed
2 - All aqueous species, solid species, and sorbed species
concentrations will be printed
IPHAD - Flag indicating whether surfactant adsorption is pH dependent or not.

Possible Values:
0 - Surfactant adsorption is not pH dependent
1 - Surfactant adsorption is pH dependent
3.5.2 PHC, PHT, PHT1, HPHAD (This line is read only if IPHAD>0)
PHC - Critical pH above which surfactant adsorption is pH dependent.
PHT - Extrapolated pH value at zero surfactant adsorption.
PHT1 - pH value above which surfactant adsorption is constant.
HPHAD - Fraction of the low-pH adsorption plateau retained at a pH above

PHT1.
3.5.3 CSELP, CSEUP (This line is read only if IREACT=3)
CSELP - Lower optimum salinity limit for generated surfactant.

Units: meq/ml water
CSEUP - Upper optimum salinity limit for generated surfactant.

Units: meq/ml water
89
UTCHEM User's Guide
3.5.4 NELET, NFLD, NSLD, NSORB, NACAT, ICHRGE
NELET - Total number of elements less nonreacting element.
NFLD - Total number of fluid species.
NSLD - Total number of solid species.
NSORB - Total number of sorbed species.
NACAT - Total number of surfactant associated cations.
ICHRGE - Flag indicating whether an oxygen balance or a charge balance will be

used.
Possible Values:
0 - Oxygen balance used
1 - Charge balance in solution used
Note: If solid SiO2 is considered, the oxygen balance must be used.
3.5.5 NIAQ, NEX, NSLEL, NSURF1
NIAQ - Total number of independent fluid species.
NEX - Total number of insoluble exchangers.
NSLEL - Total number of elements comprising the solid species.
NSURF1 - Position number corresponding to the insitu generated surfactant anion

in the fluid species array FLDSPS.
Note: NSURF1 is automatically set to 0 by the program if IREACT=2.
3.5.6 NH, NNA, NCA, NMG, NCARB
NH - Position number corresponding to the hydrogen element in the element array

ELEMNT.
NNA - Position number corresponding to the sodium element in the element array
ELEMNT.
NCA - Position number corresponding to the calcium element in the element array
ELEMNT.
NMG - Position number corresponding to the magnesium element in the element

array ELEMNT.
NCARB - Position number corresponding to the carbonate pseudo-element in the

element array ELEMNT.
Note: A value of zero is required if the element is not considered.
3.5.7 NALU, NSILI, NOXY
NALU - Position number corresponding to the aluminum element in the element

array ELEMNT.
90
UTCHEM User's Guide
NSILI - Position number corresponding to the silicon element in the element array
ELEMNT.
NOXY - Position number corresponding to the oxygen element in the element

array ELEMNT.
Note: A value of zero is required if the element is not considered.
3.5.8 NACD (This line is read only if IREACT=3)
NACD - Position number corresponding to the petroleum acid pseudo-element in

the element array ELEMNT.
3.5.9 NCR, NHFD, NCRFD (This line is read only if IREACT=4)
NCR - Position number corresponding to the chromium element in the element

array ELEMNT.
NHFD - Position number corresponding to the hydrogen ion element in the fluid
species array FLDSPS.
NCRFD - Position number corresponding to the CR(III) ion in the fluid species
array FLDSPS.
3.5.10 ELEMNT(I), ELCRG(I) for I=1, NELET
ELEMNT(I) - Name of the Ith element.
ELCRG(I) - Charge for Ith element

Note: The name of each element may not exceed 32 characters. Each name and
its charge must be on a separate line of the input file. The order in which
these elements must be listed corresponds to the order in which the
component names are specified on input line 3.1.24 with the exceptions of
calcium and chloride (if they exist) since Components 5 and 6 are reserved
for these elements.
3.5.11 FLDSPS(I), for I=1, NFLD
FLDSPS(I) - Name of the Ith fluid species.

Note: The name of each fluid species may not exceed 32 characters and
each name must be on a separate line of the input file. If
IREACT=3, the last fluid species must be HAw (petroleum acid in
water).
3.5.12 SLDSPS(I), for I=1, NSLD (This line is read only if NSLD>0)
SLDSPS(I) - Name of the Ith solid species.

Note: The name of each solid may not exceed 32 characters and each
name must be on a separate line of the input file.
3.5.13 SORBSP(I), for I=1, NSORB (This line is read only if NSORB>0)
SORBSP(I) - Name of the Ith adsorbed cation.
91
UTCHEM User's Guide
Note: The name of each adsorbed cation may not exceed 32 characters
and each name must be on a separate line of the input file.
3.5.14 ACATSP(I), for I=1, NACAT (This line is read only if NACAT>0)
ACATSP(I) - Name of the Ith surfactant adsorbed cation.

Note: The name of each surfactant adsorbed cation may not exceed 32
characters and each name must be on a separate line of the input
file.
3.5.15 NSORBX(I), for I=1, NEX (This line is read only if NSORB>0)
NSORBX(I) - Number of cations for Ith exchanger.

3.5.16 AR(I,J), for J=1, NFLD, for I=1, NELET or >>
AR(I,J), for J=1, NFLD, for I=1, NELET-1
AR(I,J) - Stoichiometric coefficient of Ith element in Jth fluid species.
Note: If ICHRGE=0, then NFLD × NELET values are required by the program.
If ICHRGE=1, then NFLD × (NELET-1) values are required by the
program.
3.5.17 BR(I,J), for J=1, NSLD, for I=1, NELET << or >>
BR(I,J), for J=1, NSLD, for I=1, NELET-1 (This line is read only if NSLD>0)
BR(I,J) - Stoichiometric coefficient of Ith element in Jth solid species.
Note: If ICHRGE=0, then NSLD × NELET values are required by the program.
If ICHRGE=1, then NSLD × (NELET-1) values are required by the
program.
3.5.18 DR(I,J), for J=1, NSORB, for I=1, NELET << or >>
DR(I,J), for J=1, NSORB, for I=1, NELET-1 (This line is read only if
NSORB>0)
DR(I,J) - Stoichiometric coefficient of Ith element in Jth sorbed species.
Note: If ICHRGE=0, then NSORB × NELET values are required by the

program. If ICHRGE=1, then NSORB × (NELET-1) values are required
by the program.
3.5.19 ER(I,J), for J=1, NACAT, for I=1, NELET << or >>
ER(I,J), for J=1, NACAT, for I=1, NELET-1 (This line is read only if
NACAT>0)
ER(I,J) - Stoichiometric coefficient of Ith element in Jth surfactant associated

cation.
92
UTCHEM User's Guide
Note: If ICHRGE=0, then NACAT × NELET values are required by the

program. If ICHRGE=1, then NACAT × (NELET-1) values are required
by the program.
3.5.20 BB(I,J), for J=1, NIAQ+NSORB+NACAT, for I=1, NFLD+NSORB+NACAT
BB(I,J) - Exponent of the Jth independent fluid species concentration when the
I th fluid species is expressed in terms of independent species
concentrations.
3.5.21 EXSLD(I,J), for J=1, NIAQ, for I=1, NSLD (This line is read only if NSLD>0)
EXSLD(I,J) - Exponent of the Jth independent fluid species concentration in the

solubility product definition of the Ith solid.
3.5.22 CHARGE(I), for I=1, NFLD
CHARGE(I) - Charge of the Ith fluid species.

3.5.23 SCHARG(I,J), for J=1, NSORBX(I), for I=1, NEX (This line is read only if
NSORB>0)
SCHARG(I,J) - Charge of the Jth sorbed species on the Ith exchanger.

3.5.24 EQK(I), for I=1, NFLD
EQK(I) - Equilibrium constant for I th fluid species when expressed in

independent species concentrations only.
3.5.25 EXK(I,J), for J=1, NSORBX(I)-1, for I=1, NEX (This line is read only if
NEX>0)
EXK(I,J) - Exchange equilibrium constant for Jth exchange equilibrium of the Ith
insoluble exchanger.
3.5.26 EXEX(I,J,K), for K=1, NIAQ+NSORB+NACAT, for J=1, NSORBX(I)-1, for

I=1, NEX (This line is read only if NEX>0)
EXEX(I,J,K) - Exponent of Kth independent species in Jth equilibrium relation

of the Ith exchanger
3.5.27 REDUC(I,J), for J=1, NSORBX(I)-1, for I=1, NEX (This line is read only if
NEX>0)
REDUC(I,J) - Valence difference of the two cations involved in the exchange

reaction J on exchanger I.
Note: This value is positive if the higher valence cation bulk concentration has a
positive exponent in EXEX(I,J) definition and is negative otherwise.
93
UTCHEM User's Guide
3.5.28 EXCAI(I), for I=1, NEX (This line is read only if NEX>0)
EXCAI(I) - Exchange capacity of Ith insoluble exchanger.

Units: meq/ml pore volume
3.5.29 SPK(I), for I=1, NSLD (This line is read only if NSLD≥1)
SPK(I) - Solubility product of Ith solid defined in terms of independent fluid

species concentrations only.
3.5.30 CHACAT(I), for I=1, NACAT (This line is read only if NACAT≥1)
CHACAT(I) - Charge of Ith surfactant associated cation.

3.5.31 ACATK(I), for I=1, NACAT-1 (This line is read only if NACAT≥1)
ACATK(I) - Equilibrium constant for Ith exchange equilibrium for cation

exchanges on surfactant.
3.5.32 EXACAT(I,J) for J=1, NIAQ+NSORB+NACAT, for I=1, NACAT-1 (This line
is read only if NACAT≥1)
EXACAT(I,J) - Exponent of Jth independent species in Ith equilibrium for cation

exchange on surfactant.
3.5.33 CI(J), for J=1, NACAT (This line is read only if NACAT≥1)
CI(J) - Initial concentration of Jth surfactant associated cation.

Units: moles/liter pore volume
3.5.34 C5I, C6I
C5I - Initial concentration of nonreacting anions.

Units: equivalents/liter
C6I - Initial concentration of calcium in aqueous phase.

3.5.35 CELAQI(J), for J=1, NGC
CELAQI(J) - Initial concentrations of Jth geochemistry component.

3.5.36 CAC2I (This line is read only if IREACT=3)
CAC2I - Initial concentration of acid in oil.

Units: moles/liter oil
3.5.37 CAQI(J), for J=1, NIAQ
CAQI(J) - Initial guesses for Jth independent species concentration.

Units: moles/liter water
94
UTCHEM User's Guide
3.5.38 CSLDI(I), for I=1, NSLD (This line is read only NSLD>0)
CSLDI(I) - Initial concentration of Ith solid.

3.5.39 CSORBI(I), for I=1, NSORB (This line is read only if NSORB>0)
CSORBI(I) - Initial concentration of Ith adsorbed cation.

3.5.40 C1I, C2I (This line is read only if IREACT=3)
C1I - Initial concentration of water in aqueous phase.

C2I - Initial concentration of oil in oleic phase.

3.5.41 ACIDIS, EQWPS (This line is read only if IREACT=3)
ACIDIS - Dissociation constant of the petroleum acid.
EQWPS - Equivalent weight of petroleum acid.
95
UTCHEM User's Guide
User notes:
96
UTCHEM User's Guide
3.6 Data for Biodegradation Option
The sixth input section consists of physical property data that is required
only if IBIO=1. This section includes the biodegradation and mass transfer
parameters required to model the biodegradation of chemical species. This section
is read only if IBIO=1. The data is read by a separate subroutine called BIOREAD,
and is input in the standard UTCHEM format. Section 9 of the UTCHEM
Technical Documentation gives more details on this option.
3.6.1 DENBLK
DENBLK - Bulk density of the porous medium (mass of porous medium per unit
total volume).
Units: g/cm3
3.6.2 CMIN, EPSBIO, IBTIM
CMIN - Minimum concentration of substrate and electron acceptor that is of

interest. This parameter is used for two purposes. First, if concentrations
of all substrates and electron acceptors in a gridblock are below CMIN, then
biodegradation reactions are assumed negligible at the gridblock and are not
modeled. Second, when the concentration of all substrates and electron
acceptors fall below CMIN during solution of the biodegradation reaction
expressions, further biodegradation reactions are assumed to be negligible
and program execution returns to the main program.
Units: mg/L
EPSBIO - Convergence tolerance for solution of the biodegradation equations.

Note: Values of 10-4 to 10-6 are recommended, although larger values can
also result in accurate simulations. Small values ensure accurate
solutions of the biodegradation equations but increase run times,
while larger values decrease run times at the expense of some
accuracy.
IBTIM - Flag indicating type of time step control for solution of biodegradation
equations.
Possible Values:
0 - No time step control; biodegradation equations are solved at
every transport time step
1 - Manual time step control; biodegradation time step is specified
by user
2 - Automatic time step control; biodegradation time step is
controlled by UTCHEM based on an acceptable error
specified by the user
97
UTCHEM User's Guide
3.6.3 BTMIN, BTMAX (This line is read only if IBTIM=1)
BTMIN - Minimum time step for solution of biodegradation equations.
BTMAX - Maximum potential biodegradation time step.
Note: If the transport time step is smaller than BTMIN, the biodegradation
equations will be solved at intervals of BTMIN. If the transport time step
falls between BTMIN and BTMAX, the biodegradation equations will be
solved at every transport time step. BTMAX should be set to a value
greater than BTMIN and is only used to initialize variables in this option.
3.6.4 BIORME, BTMAX (This line is read only if IBTIM=2)
BIORME - The maximum acceptable error in terms of the mass of substrate in the
system resulting from a biodegradation time step being larger than the
transport time step.
BTMAX - Maximum biodegradation time step allowed.
Note: A warning message is presented if the biodegradation time step computed

from the automatic time step controller is smaller than the transport time
step. See Section 9 of the UTCHEM Technical Documentation for more
information.
3.6.5 NBC, NMET, IBKIN, IBPP
NBC - Total number of chemical and biological species that are considered in
biodegradation reactions, including oil components, surfactants, products
generated by abiotic and biodegradation reactions, nutrients required for
biological growth, electron acceptors, and biological species.
NMET - Number of substrate/electron acceptor/biological species metabolic

combinations. Include combinations of biodegrading products/electron
acceptor/biological species for each product that also biodegrades.
IBKIN - Flag specifying the type of biodegradation kinetics.

Possible Values:
0 - No reaction; biodegradation parameters are read, but
biodegradation equations are never solved (useful for restart
runs)
1 - Monod kinetics with external mass transfer resistances
2 - Monod kinetics with no external mass transfer resistances
3 - Instantaneous kinetics (stoichiometric reactions)
4 - Monod kinetics with automatic control of external mass transfer
resistances
Note: First order kinetics with biomass growth can be also be modeled by
adjusting the values of the Monod parameters. See input line 3.6.12
below.
IBPP - Flag indicating whether or not porosity and permeability are affected by
biological growth.
Possible Values:
0 - Porosity and permeability are not affected by biomass growth
98
UTCHEM User's Guide
1 - Porosity and permeability are affected by biomass growth
3.6.6 IMTVAR (This line is read only if IBKIN=4)
IMTVAR - Flag indicating type of mass transfer control.

Possible Values:
1 - Mass transfer is considered in each gridblock at each time step
only when the Damkohler number is greater than the user-
specified value
2 - Mass transfer is considered in each gridblock at each time step
only when the effectiveness factor is less than the user-
specified value
3.6.7 DAMX (This line is read only if IBKIN=4 and IMTVAR=1)
DAMX - Value of Damkohler number used to control mass transfer in

biodegradation calculations.
Recommended value: 0.1
3.6.8 EFMIN (This line is read only if IBKIN=4 and IMTVAR=2)
EFMIN - Value of effectiveness factor used to control mass transfer in

biodegradation calculations.
Recommended value: greater than 0.95
3.6.9 KC(I), ITYPE(I), CINIT(I), RABIO(I), NPABIO(I), for I=1, NBC
Note: One line is required for each chemical and biological species that participates
in biodegradation reactions.
KC(I) - Index of the Ith chemical or biological species.
ITYPE(I) - Flag indicating whether the I'th component is a chemical or biological

species.
Possible values:
1 - The Ith component is a chemical species
2 - The Ith component is a biological species
CINIT(I) - Initial concentration of chemical or free-floating (unattached) biological
species I in the aqueous phase.
Units: mg/l
Note: Although a value of CINIT must be entered for organic species that
participate in biodegradation reactions, these values are ignored by
the program. Initial concentrations of these components are input on
input lines 3.3.32 through 3.3.35.
RABIO(I) - First-order abiotic reaction rate constant.

Units: 1/days
Note: Although a value of RABIO can be specified for biological as well
as chemical species, biomass decay should not be controlled with
RABIO. Instead, use the parameters ENDOG and ENDOGB to
control endogenous decay of unattached and attached biomass,
99
UTCHEM User's Guide
respectively. RABIO should normally be set to 0.0 for biological

species.
NPABIO(I) - Number of products generated by a first-order abiotic reaction of

chemical species I.
Note: A value must be entered for biological species as well, although the
value is ignored by the program because generation of products
from decay of biomass is not allowed.
3.6.10 KC(I), DENBIO(I), RCOL(I), TCOL(I), COLNUM(I), ENDOG(I),

ENDOGB(I), CBI(I), CBIOMN(I), ADSBIO(I), for I=1, NBS
Note: One line is required for each biological species.
KC(I) - Index of the biological species.
DENBIO(I) - density of attached biological species I (biofilm density).

Units: g cells / cm3 biomass
RCOL(I) - radius of an attached microcolony of biological species I.
Units: cm
Note: The parameter RCOL is used to calculate the surface area of a single
attached microcolony. Although microcolonies are assumed to be
disk-shaped by the model, the user may specify any desired surface
area per microcolony using RCOL. TCOL can then be adjusted to
obtain the desired volume of the microcolony since the thickness of
the microcolony does not affect external mass transfer into the
attached biomass.
TCOL(I) - Thickness of a single attached microcolony of biological species I.

Units: cm
COLNUM(I) - Number of bacterial cells per microcolony of biological species I.

Units: cells/colony
ENDOG(I) - endogenous decay coefficient of unattached cells of biological

species I.
Units: 1/days
ENDOGB(I) - endogenous decay coefficient of attached cells of biological

species I.
Units: 1/days
CBI(I) - Number of attached bacterial cells of biological species I per gram of dry
soil.
Units: cells/gram of solid
CBIOMN(I) - Lower limit of number of attached bacterial cells of biological

species I.
Units: cells/gram of solid
Note: A population of attached microorganisms, sustained by naturally
occurring organic matter, is assumed to exist in the porous medium
regardless of the concentration of other chemical species. This
100
UTCHEM User's Guide
concentration is CBIOMN, and the concentration of biomass is not

allowed to fall below this value.
ADSBIO(I) - Biomass partitioning coefficient.

mass of attached microorganisms
Units:
mass of unattached microorganisms
3.6.11 ISUB(I), IEA(I), IBS(I), BRMAX(I), BRMAXB(I), YXS(I), AKS(I), AKA(I),

FEA(I), for I=1, NMET
Note: One line is read for each metabolic combination. Values of all parameters
on this line must be supplied by the user even if instantaneous kinetics or
abiotic reactions are used to describe loss of substrate.
ISUB(I) - Substrate index for metabolic combination I.
IEA(I) - Electron acceptor index for metabolic combination I.
IBS(I) - Biological species index for metabolic combination I.
BRMAX(I) - maximum specific growth rate of unattached microorganisms for

metabolic combination I.
Units: 1/days
BRMAXB(I) - Maximum specific growth rate of attached microorganisms for

Units: 1/days
YXS(I) - Yield coefficient for metabolic combination I.

Units: biomass produced per mass of substrate biodegraded
AKS(I) - Substrate half-saturation coefficient for metabolic combination I.

Units: mg/l
AKA(I) - Electron acceptor half-saturation coefficient for metabolic combination

I.
Units: mg/l
FEA(I) - Electron acceptor utilization coefficient.

Units: mass of electron acceptor consumed per mass of substrate
biodegraded
3.6.12 ISUB(I), IEA(I), IBS(I), NCOMPS(I), NIHB(I), NPROD(I), NNUT(I),

ICOMET(I), for I=1, NMET
Note: One line is read for each metabolic combination.
101
UTCHEM User's Guide
NCOMPS(I) - Number of other substrates competing with substrate ISUB in

NIHB(I) - Number of other chemical species that inhibit metabolic combination I.
NPROD(I) - Number of products generated from metabolic combination I.
NNUT(I) - Number of nutrients that limit the biodegradation rate through Monod
terms for metabolic combination I.
ICOMET(I) - Flag indicating whether or not the substrate in metabolic

combination I is biodegraded through aerobic cometabolism.
Possible Values:
0 - Substrate ISUB(I) serves as a primary substrate
1 - Substrate ISUB(I) is biodegraded through aerobic
cometabolism
Note: Users may specify multiple competing substrates, nutrients, and inhibiting
constituents for each metabolic combination. Substrate competition,
inhibition and cometabolism cannot be modeled if instantaneous kinetics
are selected. However, values for the biodegradation rate parameters must
be specified even if instantaneous kinetics are specified.
First-order biodegradation kinetics can be modeled by using a very large

value of KS, and adjusting the ratio of µmax/KS to be equal to the desired
first-order biodegradation rate coefficient.
3.6.13 ISUB(I), IEA(I), IBS(I), (ICSUB(I,J), for J=1, Number of competing

substrates), for I=1, NMET (This line is read only if there are competing
substrates)
Note: One line is read for each metabolic combination for which there is substrate
competition between two or more substrates.
ICSUB(J) - Indices of other substrates that compete with substrate ISUB(I) in

Note: The number of input lines required must equal the number of species that
are competing, since complementary lines are required to fully describe the
competition. For example, if substrate 12 in metabolic combination 12-15-
16 must compete with substrates 13 and 14 that are also biodegraded by
biological species 16 using electron acceptor 15, then the required input
lines are:
12 15 16 13 14
13 15 16 12 14
14 15 16 12 13
102
UTCHEM User's Guide
3.6.14 ISUB(I), IEA(I), IBS(I), IHB(I), BSIHB(I), for I=1, Number of metabolic
combination and inhibiting compound associations (This line is only read if there
are metabolic combination and inhibiting compound associations)
Note: One input line is read for each association of metabolic combination and
inhibiting compound.
IHB(I) - Index of chemical species that inhibits metabolic combination I.
BSIHB(I) - Inhibition constant for metabolic combination I.

Units: mg/l
Note: The total number of lines are

NMET
I= ∑ NIHB(J) .
J=1
For example, suppose there are two substrate-electron acceptor-biological
species metabolic combinations: 9-10-12 and 9-11-12. Metabolic
combination 9-10-12 is inhibited by the substrate itself (9) and electron
acceptor 11, while metabolic combination 9-11-12 is inhibited by only the
substrate. Assuming that the imbibition constant for each combination is
0.001, the input lines for this scenario would be:
9 10 12 9 0.001
9 10 12 11 0.001
9 11 12 9 0.001
3.6.15 ISUB(I), IEA(I), IBS(I), IPR(I), FPR(I), for I=1, Number of biodegradation
product formation and metabolic combination associations (This line is read only
if there are biodegradation product formation described with Monod kinetics and
metabolic combination associations)
Note: One line is read for each association of product formation and metabolic
combination.
IPR(I) - Index of product generated by metabolic combination I.
FPR(I) - Product generation coefficient (stoichiometric ratio - mass of product

generated per mass of substrate biodegraded).
103
UTCHEM User's Guide
Note: Number of lines are

NMET
I= ∑ NPRODT(J) .
J=1
For example, if metabolic combination 9-14-15 generates products 10 and
11, and metabolic combination 12-14-15 generates product 13, then the
input lines would be (assuming a values of 0.5 for FPR(I)):
9 14 15 10 0.5
9 14 15 11 0.5
12 14 15 13 0.5
Parameters for generation of products through first-order reactions are

described on input line 3.6.16.
3.6.16 ISUB(I), IPR(I), FPR(I), for I=1, Number of products generated by first-order
reactions of the biodegradation species. (This line is read only if there are products
of first-order reactions of biological species)
ISUB(I) - Index of chemical species that reacts abiotically to generate a product.
IPR(I) - Index of product generated by abiotic reaction of ISUB(I).
FPR(I) - Product generation coefficient (stoichiometric ratio - mass of product

generated per mass of reactant reacted).
Note: Number of lines are

NBC
I= ∑ NPABIO(J) .
J=1
For example, if the abiotic products 11, 12 and 13 were generated from the
first-order reaction of biodegradation species 9, the input lines would be:
9 11 1.0
9 12 1.0
9 13 1.0
3.6.17 ISUB(I), IEA(I), IBS(I), INUT(I), AKN(I), FN(I), for I=1, Number of metabolic
combination and limiting nutrient associations (This line is read only if there are
metabolic combination and limiting nutrient associations)
Note: One line is read for each association of metabolic combination and limiting
nutrient.
INUT(I) - Index of nutrient limiting the rate of biodegradation through a Monod

term in metabolic combination I.
104
UTCHEM User's Guide
AKN(I) - Nutrient half-saturation coefficient for metabolic combination I.

Units: mg/l
FN(I) - Nutrient utilization coefficient (mass of nutrient consumed per mass of

substrates biodegraded) for metabolic combination I.
NMET
Note: Number of lines are I = ∑ NNUTT(J) .
J=1
3.6.18 ISUB(I), IEA(I), IBS(I), TC(I), IRLIM(I), for I=1, Number of cometabolic
combinations for which aerobic cometabolism exists (This line is read only if there
are cometabolic combinations for which aerobic cometabolism exists for at least
one metabolic combination specified on input line 3.6.12)
Note: One input line is required for each cometabolic combination for which
aerobic cometabolism exists.
TC(I) - Transformation capacity for cometabolism of substrate ISUB(I) (mass of

substrate utilized per mass of biomass destroyed).
IRLIM(I) - Flag indicating whether reducing power limitations are considered for
cometabolic combination I.
Possible Values:
0 - No reducing power limitations are considered
1 - Cometabolic reaction consumes reducing power
Note: Reducing power limits the biodegradation rate through Monod
terms in the manner of Chang and Alvarez-Cohen, 1995. The loss
of a biological species' reducing power reduces its activity toward
all substrates, not just the cometabolite.
NMET
Note: The total number of lines are I = ∑ COMET(J) .
J=1
3.6.19 ISUB(I), IEA(I), IBS(I), IGROW(I), REDI(I), AKR(I), FRP(I), FRC(I), for
I=1, Number of cometabolic biodegradation reactions in which reducing power
limitations are considered (This line is read only if there are cometabolic
biodegradation reactions in which reducing power limitations are considered and
IRLIM(I)>0 for at least one metabolic combination specified on input line 3.6.17)
Note: One line is required for each cometabolic biodegradation reaction in which
reducing power limitations are considered.
105
UTCHEM User's Guide
IGROW(I) - Index of growth substrate for cometabolism of substrate ISUB(I).
REDI(I) - Initial intracellular reducing power (NADH) concentration.

Units: mmol/mg of biomass
AKR(I) - Reducing power half-saturation coefficient.

Units: mmol/mg of biomass
FRP(I) - Reducing power generation coefficient for metabolic combination

(IGROW(I), IEA(I), IBS(I) (mmol reducing power generated per mg of
growth substrate consumed).
FRC(I) - Reducing power consumption coefficient for cometabolic combination

(ISUB(I), IEA(I), IBS(I)) (mmol reducing power consumed per mg of
substrate consumed by cometabolism).
Note: The total number of lines are

K
I= ∑ IRLIM(J)
J=1
where
NMET
K= ∑ COMET(J) .
J=1
106
UTCHEM User's Guide
3.7 Recurrent Injection/Production Data Set
The sixth input section consists of the recurrent injection/production well
data. Please remember that there are seven comment lines at the beginning of this
section and that each line is preceded by three comment lines.
3.7.1 IBOUND, IZONE
IBOUND - Flag indicating if constant potential boundaries at the top, left, and
right sides of the simulation model are specified.
Possible Values:
0 - No boundary is specified
1 - Boundary is specified
IZONE - Flag indicating the aquifer zone to be modeled.

Possible Values:
0 - Saturated zone is modeled
1 - Vadose zone is modeled
2 - Both saturated and vadose zones are modeled
Note: See Section2 of the UTCHEM Technical Documentation for more

information.
3.7.2 IBL, IBR, IBT, WTABL, WTABR (This line is read only if IBOUND=1)
IBL - Flag indicating if the left hand side constant potential boundary is specified.
Possible Values:
IBR - Flag indicating if the right hand side constant potential boundary is
specified.
Possible Values:
IBT - Flag indicating if the top constant pressure boundary of the vadose zone is
specified.
Possible Values:
0 - Top boundary is closed
1 - Top boundary is open
WTABL - The depth to the water table at the left boundary.

WTABR - The depth to the water table at the right boundary.

Note: The depth to the water table is only used if IZONE=2.
107
UTCHEM User's Guide
3.7.3 PBL, C1BL, C5BL, C6BL (This line is read only if IBOUND=1 and IZONE≠1
and IBL=1)
PBL- Pressure at the center of the top layer at the left boundary.
C1BL- Concentration of water in aqueous phase at the left boundary.

C5BL - Concentration of chloride in aqueous phase at the left boundary.

Units: meq/ml water
C6BL - Concentration of calcium in aqueous phase at the left boundary.

Units: meq/ml water
3.7.4 PBR, C1BR, C5BR, C6BR (This line is read only if IBOUND=1 and IZONE≠1
and IBR=1)
PBR - Pressure at the center of the top layer at the right boundary.
C1BR - Concentration of water in aqueous phase at the right boundary.

C5BR - Concentration of chloride in aqueous phase at the right boundary.

Units: meq/ml water
C6BR - Concentration of calcium in aqueous phase at the right boundary.

Units: meq/ml water
Note: For the biodegradation option (IBIO=1), the concentrations of all species
considered at the boundary are set to the initial concentrations. If the
variation in temperature is considered (IENG=1), the temperatures at the
boundaries are set to the initial temperature.
3.7.5 NWELL, IRO, ITSTEP, NWREL
NWELL - Number of wells used for the simulation including the pseudowells to
mimic an open boundary.
Note: If ICOORD=2, NWELL must be equal to 1 and the NWELL input
file in the HEAD file must be set equal to 2.
IRO - Flag indicating the equivalent well radius model to be used.

Possible Values:
1 - Babu and Odeh model is used
2 - Peaceman model is used
Note: The Babu and Odeh model (IRO=1) does not work for ICOORD=4.
ITSTEP - Flag indicating the units to be used when specifying the minimum and
maximum time step.
Possible Values:
0 - Minimum and maximum time steps are input in days
108
UTCHEM User's Guide
1 - Minimum and maximum time steps are input as Courant

numbers
Note: This option is only used if IMES>1 and is not a shut-in period. If
IMES=1, this flag is ignored. For a shut-in period you need to use
ITSTEP=0
NWREL - Number of actual wells used for the simulation excluding the
pseudowells.
Note: The history data are written only for NWREL wells.
Note: See Appendix E for more details on the Courant number and time step
selection options.
The following values for minimum and maximum Courant numbers are
recommended for different simulations as follows:
Process Min. Courant # Max. Courant #

Waterflood/tracer 0.04 0.4
Polymerflood 0.02 0.2
Surfactant/polymerflood 0.01 0.1
Geochemical process 0.01 0.1
See Section 10 of the UTCHEM Technical Documentation for well model

information.
3.7.6 The data on input lines 3.7.6.a through 3.7.6.d are repeated for M=1 to NWELL
times.
Important note: Input the actual wells first (NWREL) and then the
pseudowells.
3.7.6.a IDW(M), IW(M), JW(M), IFLAG(M), RW(M), SWELL(M), IDIR(M),

IFIRST(M), ILAST(M), IPRF(M)
IDW(M) - Well I.D. number for the Mth well.

Possible Values: Must be between 1 and MXW (the source code
parameter indicating the maximum number of wells)
Note: This number is used by UTCHEM to keep track of which well is
being described in the recurrent injection/production well section.
IW(M) - First index of the reservoir gridblock containing the Mth well.
Possible Values: Between 1 and the number of gridblocks in the pertinent
direction, inclusive
Note: If the Mth well is completed parallel to the X-axis, IW(M) is the Y
direction index—if the well is completed parallel to the Y- or Z-
axis, IW(M) is the X direction index. See example below.
If ICOORD=2, IW(1)=JW(1)=1.
JW(M) - Second index of the reservoir gridblock containing the Mth well.
109
UTCHEM User's Guide
Note: If the Mth well is completed parallel to the X- or Y-axis, JW(M) is

the Z direction index—if the well is completed parallel to the Z-
axis, JW(M) is the Y direction index. See example below.
If ICOORD=2, IW(M)=JW(M)=1.
IFLAG(M) - Flag indicating type of well constraint specification for Mth well.
Possible Values:
1 - Rate constrained injection well
2 - Pressure constrained production well (This option is
available only if ICOORD=1 or 3)
3 - Pressure constrained injection well (This option is available
only if ICOORD=1 or 3)
4 - Rate constrained production well
RW(M) - Radius of Mth well.

SWELL(M) - Skin factor for Mth well.

IDIR(M) - Flag indicating the direction in which the Mth well is completed.
Possible Values:
1 - Well completed parallel to the X-axis
2 - Well completed parallel to the Y-axis
3 - Well completed parallel to the Z-axis
Note: If ICOORD=2, IDIR(M) must be equal to 3.
IFIRST(M) - Index of the first block in which the Mth well is completed.
ILAST(M) - Index of the last block in which the Mth well is completed.
Possible Values: Between IFIRST(M) and the number of gridblocks in
the pertinent direction, inclusive
IPRF(M) - Flag indicating if partial completion of the well is considered.

Possible Values:
0 - The well is fully completed
1 - The well is partially completed
Example: For a vertical well (completed through all the layers) as illustrated in the
4 × 4 × 3 example below, note the values of IDIR(M), IW(M), JW(M),
IFIRST(M), and ILAST(M):
110
UTCHEM User's Guide
X
Z
Y
IDIR(M) = 3
IW(M) = 1
JW(M) = 1
IFIRST(M) = 1
ILAST(M) = 3
For a horizontal well (completed from the first to last gridblock in the X
direction and parallel to the X-axis) as illustrated in the 4 × 4 × 3 example
below, note the values of IDIR(M), IW(M), JW(M), IFIRST(M), and
ILAST(M):
X
Z
Y
IDIR(M) = 1
IW(M) = 2
JW(M) = 1
IFIRST(M) = 1
ILAST(M) = 4
Note: Horizontal wells can be used for 2-D X-Y or 3-D simulations.
3.7.6.b KPRF(M,IWB), for IWB=1, NWBC (This line is read only if IPRF=1)
KPRF(M,IWB) - Flag indicating if the IWBth well block of the Mth well is
perforated or not.
Possible Values:
0 - The well block is not perforated
1 - The well block is perforated
3.7.6.c WELNAM(M)
WELNAM(M) - Name of the Mth well.

Note: The name can consist of any combination of up to 8 alphanumeric
characters. This information will be printed—along with the well
I.D. number, IDW(M)—at the beginning of the history output
files.
111
UTCHEM User's Guide
3.7.6.d ICHEK(M), PWFMIN(M), PWFMAX(M), QTMIN(M), QTMAX(M)
ICHEK(M) - The flag to specify whether to check the rate or pressure caps for
the Mth well.
Possible Values:
0 - There will be no check on the rate or pressure limits and no
automatic shut in for the pressure constraint injector
1 - There will be no automatic shut in for the pressure constraint
injector but the pressure or rate limits are checked
2 - There will be both the automatic shut in and the check on the
pressure or rate limits
PWFMIN(M) - Minimum flowing bottom hole pressure (specified at the top

layer) for the Mth well.
PWFMAX(M) - Maximum flowing bottom hole pressure (specified at the top

layer) for the Mth well.
QTMIN(M) - Minimum total flow rate (specified at the top layer) for the Mth
well.
QTMAX(M) - Maximum total flow rate (specified at the top layer) for the Mth
well.
Note: QTMIN(M) and QTMAX(M) need to be input as negative values for
production wells (IFLAG=2 or 4).
PWFMIN(M) and PWFMAX(M) are the pressure caps for a rate

constraint injector or producer well. QTMIN(M) and QTMAX(M) are the
total rate caps for a pressure constraint injector or producer well. If the
Mth pressure constraint injector or producer produces at total rate less than
QTMIN(M), the Mth well will be switched to a rate constraint well with
total rate of QTMIN(M) for the rest of the injector or production period.
On the other hand, if the total rate is greater than the QTMAX(M), the Mth
well then will be switched to a rate constraint well with the total rate of
QTMAX(M). The similar concept is applied to a rate constraint injector or
producer.
The user can skip the well control calculation by specifying very small
values for QTMIN(M) and PWFMIN(M) and very large values for
QTMAX(M) and PWFMAX(M).
The code still has the automatic option for shut in of a pressure constraint
injector injecting at a rate of less than QTMIN(M).
112
UTCHEM User's Guide
3.7.7 The data on input lines 3.7.7.a, 3.7.7.b, 3.7.7.c, and 3.7.7.d are repeated for
M=1 to NWELL times.
Note: For injection wells that are on rate constraint only injection rates and
concentrations for each phase are listed. For injection wells that are on
pressure constraint the injection pressure is also specified. In this case the
injection rates are treated as phase cuts in the injected fluid. For producer
pressure constraint only the bottom hole pressure is specified. For
producer rate constraint only the total production rate is specified.
The user can shut in a pressure constraint well by specifying a negative

bottom hole pressure or a rate constraint well by specifying a value of zero
for rate (QI).
3.7.7.a ID(M), (QI(M,L), (C(M,KC,L), for KC=1, N), for L=1, MXP) (This set of
data is read only if IFLAG(M)=1 or 3)
ID(M) - Well I.D. number for the Mth well.

Note: See note for IDW(M) on input line 3.7.6.a.
QI(M,L) - Injection rate of Lth phase in Mth well (see note below).
C(M,KC,L) - Concentration of KCth component in Lth phase in Mth well.

Units: vary according to component (see note below)
Note: See Appendix D for component and phase numbering scheme and the
concentration units for each species.
The KC index changes the fastest, the L index changes the next fastest,
and the M index changes the slowest. A separate data line should be in
the input file for each phase—that is, M × L lines will be read by the
program. MXP is equal to 3 (IGAS=0) or 4 (IGAS≥1).
3.7.7.b ID(M), PWF(M) (This line is read only if IFLAG(M)=2 or 3)

Note: See note for IDW(M) on input line 3.7.6.a. For IFLAG(M)=3,
the rates (QI(M,L)) are used an injected volume fraction for each
phase.
PWF(M) - Flowing bottom hole pressure for the Mth well.

3.7.7.c ID(M), TEMINJ(M) (This line is read only if IENG=1 and IFLAG(M)=1 or 3)
113
UTCHEM User's Guide

TEMINJ(M) - Injection temperature for Mth well.

3.7.7.d ID(M), QI(M,1) (This line is read only if IFLAG(M)=4)

QI(L) - Total production rate for Mth well.

Note: This value needs to be input as a negative number.
3.7.8 TINJ, CUMPR1, CUMHI1, WRHPV, WRPRF, RSTC
TINJ - Cumulative injection time.

CUMPR1 - Indicates interval at which profiles should be written to the PROF

output file.
Units: pore volumes or days (dependent on value of ICUMTM flag)
CUMHI1 - Indicates interval at which production data should be written to the

PROF output file.
WRHPV - Indicates interval at which production histories should be written to

output file(s) for history plotting.
Note: If WRHPV > total pore volume injected or maximum simulation
time, the data will not be printed.
WRPRF - Indicates interval at which information should be written to the output

files controlled by the flags listed on input lines 3.2.3 and 3.2.4.
Note: If WRPRF > total pore volume injected or maximum simulation
time, the data will not be written.
RSTC - Indicates the interval at which restart data should be written to the
RESTAR output file.
114
UTCHEM User's Guide
Time Step Selection Data (input lines 3.7.9-3.7.13)
See Appendix E for more details on the time step selection options.
3.7.9 DT (This line is read only if IMES=1 and ITSTEP=0)
DT - Time step size for constant time step option.

Units: days
3.7.10 DT, DCLIM, DTMAX, DTMIN (This line is read only if IMES=2 and
ITSTEP=0)
DT - Initial time step size.

Units: days
DCLIM - Tolerance for concentration change for the first three components.
DTMAX - Maximum time step size.

Units: days
DTMIN - Minimum time step size.

Units: days
3.7.11 DT, DCLIM, CNMAX, CNMIN (This line is read only if IMES=2, ITSTEP=1,
and at least one well is not shut-in.)

Units: days
DCLIM - Tolerance for concentration changes for the first three components.
CNMAX - Maximum Courant number.

CNMIN - Minimum Courant number.

3.7.12 DT, (DELC(KC), for KC=1, N), DTMAX, DTMIN (This line is read only if
IMES=3 or 4 and ITSTEP=0)
DT - Initial time-step size.

Units: days
DELC(KC) - Tolerance for concentration change of KCth component (IMES=3)

or relative tolerance for concentration change of KCth component
(IMES=4).
Units: IMES=3: volume fraction, weight percent, meq/ml, or ppm
(depending on which component the tolerance is for—see note for
input line 3.7.7.a)
IMES=4: dimensionless
115
UTCHEM User's Guide
Note: DELC(KC) is the dimensionless relative change in concentration.

For example: DELC(3)=0.1 indicates a 10% change in
concentration of component 3.

Units: days

Units: days
3.7.13 DT, (DELC(KC), for KC=1, N), CNMAX, CNMIN (This line is read only if
IMES=3 or 4, ITSTEP=1, and all the wells are not shut-in)

Units: days

(IMES=4).
input line 3.7.7.a)


IMPORTANT NOTE: The data on input lines 3.7.14 through 3.7.29 describe the
changes in boundary conditions and are repeated until the injected time (TINJ) is greater
than or equal to the maximum simulation time (TMAX).
3.7.14 IBMOD
IBMOD - Flag to indicate if there are changes to the boundary conditions in this
set of recurrent data.
Possible Values:
0 - There are no changes to the boundary conditions
1 - There are changes to the boundary conditions
3.7.15 IBOUND (This line is read only if IBMOD=1)
IBOUND - Flag indicating if constant potential boundaries at the top, left, and
right sides of the simulation model are specified.
Possible Values:
116
UTCHEM User's Guide
3.7.16 IBL, IBR, IBT, WTABL, WTABR (This line is read only if IBMOD=1 and
IBOUND=1)
IBL - Flag indicating if the left hand side constant potential boundary is specified.
Possible Values:
IBR - Flag indicating if the right hand side constant potential boundary is
specified.
Possible Values:
IBT - Flag indicating if the top constant pressure boundary of the vadose zone is
specified.
Possible Values:
0 - Top boundary is closed
1 - Top boundary is open
WTABL - The depth to the water table at the left boundary.

WTABR - The depth to the water table at the right boundary.

3.7.17 PBL, C1BL, C5BL, C6BL (This line is read only if IBMOD=1 and IBOUND=1
and IZONE≠1 and IBL=1)
PBL- Pressure at the center of the top layer at the left boundary.
C1BL- Concentration of water in aqueous phase at the left boundary.

C5BL - Concentration of chloride in aqueous phase at the left boundary.

Units: meq/ml water
C6BL - Concentration of calcium in aqueous phase at the left boundary.

Units: meq/ml water
3.7.18 PBR, C1BR, C5BR, C6BR (This line is read only if IBMOD=1 and IBOUND=1
and IZONE≠1 and IBR=1)
PBR - Pressure at the center of the top layer at the right boundary.
C1BR - Concentration of water in aqueous phase at the right boundary.

C5BR - Concentration of chloride in aqueous phase at the right boundary.

Units: meq/ml water
C6BR - Concentration of calcium in aqueous phase at the right boundary.
117
UTCHEM User's Guide
Units: meq/ml water
Note: For the biodegradation option (IBIO=1), the concentrations of all species
considered at the boundary are set to the initial concentrations. If the
variation in temperature is considered (IENG=1), the temperatures at the
boundaries are set to the initial temperature.
3.7.19 IRO, ITSTEP, (IFLAG(M), for M=1, NWELL)
IRO - Flag indicating the equivalent well radius model to be used.

Possible Values:
1 - Babu and Odeh model is used
2 - Peaceman model is used
Note: The Babu and Odeh model (IRO=1) does not work for
ICOORD=4.
ITSTEP - Flag indicating the units to be used when specifying the minimum and
maximum time step.
Possible Values:
0 - Minimum and maximum time steps are input in days
1 - Minimum and maximum time steps are input as Courant
numbers
Note: This option is only used if IMES>1 and it is not a shut-in period.
If IMES=1, this flag is ignored. For a shut-in period of time, set
ITSTEP=0.
IFLAG(M) - Flag indicating type of well constraint specification for Mth well.
Possible Values:
1 - Rate constrained injection well
2 - Pressure constrained production well (This option is available
3 - Pressure constrained injection well (This option is available
4 - Rate constrained production well
3.7.20 NWEL1
NWEL1 - Number of wells with changes in location (IW(M), JW(M)), skin,

direction, perforation, name, or minimum and maximum bottomhole
pressure or minimum or maximum rate.
3.7.21 The data on input lines 3.7.21.a through 3.7.21.d are repeated for M=1 to
NWEL1 times.
3.7.21.a ID, IW(ID), JW(ID), RW(ID), SWELL(ID), IDIR(ID), IFIRST(ID),

ILAST(ID), IPRF(ID)
ID - Well ID number with changes from the previous slug injection period.
IW(ID) - First index of the reservoir gridblock containing the IDth well.
Note: See note for input line 3.7.6.a.
JW(ID) - Second index of the reservoir gridblock containing the IDth well.
118
UTCHEM User's Guide
Note: See note for input line 3.7.6.a.
RW(ID) - Radius of IDth well.

SWELL(ID) - Skin factor for IDth well.

IDIR(ID) - Flag indicating the direction in which the IDth well is completed.
Possible Values:
1 - Well completed parallel to the X-axis
2 - Well completed parallel to the Y-axis
3 - Well completed parallel to the Z-axis
Note: If ICOORD=2, IDIR(1) must be equal to 3.
IFIRST(ID) - Index of the first block in which the IDth well is completed.
ILAST(ID) - Index of the last block in which the IDth well is completed.
IPRF(ID) - Flag indicating if partial completion of the well is considered.
Possible Values:
0 - The well is fully completed
1 - The well is partially completed
3.7.21.b KPRF(ID,IWB), for IWB=1, NWBC (This line is read only if IPRF=1)
KPRF(ID,IWB) - Flag indicating if the IWBth well block of the IDth well is
perforated or not.
Possible Values:
0 - The well block is not perforated
1 - The well block is perforated
3.7.21.c WELNAM(ID)
WELNAM(ID) - Name of the IDth well.

3.7.21.d ICHEK(ID), PWFMIN(ID), PWFMAX(ID), QTMIN(ID), QTMAX(ID)
ICHEK(ID) - The flag to specify whether to check the rate or pressure caps for
the IDth well.
Possible Values:
0 - There will be no check on the rate or pressure limits and no
automatic shut in for the pressure constraint injector
1 - There will be no automatic shut in for the pressure constraint
injector but the user specified pressure or rate limits are
checked
2 - There will be both the automatic shut in and the check on the
user specified pressure or rate limits
PWFMIN(ID) - Minimum flowing bottom hole pressure (specified at the top

layer) for the IDth well.
119
UTCHEM User's Guide
PWFMAX(ID) - Maximum flowing bottom hole pressure (specified at the top

layer) for the IDth well.
QTMIN(ID) - Minimum total flow rate (specified at the top layer) for the IDth
well.
QTMAX(ID) - Maximum total flow rate (specified at the top layer) for the IDth
well.
Note: See the note for input line 3.7.6.d.
3.7.22 NWEL2, (IDWW(J), for J=1, NWEL2)
NWEL2 - Number of wells with changes in rate, concentration or bottomhole

pressure.
IDWW(J) - ID number for Jth well with changes.

3.7.23 The data on input lines 3.7.23.a through 3.7.23.d are repeated for M=1 to
NWEL2 times.
3.7.23.a ID, QI(ID,L), (C(ID,KC,L), for KC=1,N), for L=1, MXP (This set of data is
read only if IFLAG(ID)=1 or 3)
QI(ID,L) - Injection rate of Lth phase in IDth well (see note for input line
3.7.7.a).
C(ID,KC,L) - Concentration of KCth component in Lth phase for IDth well.

Units: vary according to component (see note for input line 3.7.7.a)
Note: If IGAS=0, then MXP=3. If IGAS≥1, then MXP=4.
3.7.23.b ID, PWF(ID) (This line is read only if IFLAG(ID)=2 or 3)
PWF(ID) - Flowing bottom hole pressure for the IDth well.

3.7.23.c ID, TEMINJ(ID) (This line is read only if IENG=1 and IFLAG(ID)=1 or 3)
TEMINJ(ID) - Injection temperature for the IDth well.

120
UTCHEM User's Guide
3.7.23.d ID, QI(ID,1) (This line is read only if IFLAG(ID)=4)
QI(ID,1) - Total production rate for IDth well.

Note: This value needs to be input as a negative number.
3.7.24 TINJ, CUMPR1, CUMHI1, WRHPV, WRPRF, RSTC
TINJ - Cumulative injection time.

CUMPR1 - Indicates interval at which profiles should be written to the PROF

output file.
CUMHI1 - Indicates interval at which production data should be written to the

PROF output file.
WRHPV - Indicates interval at which production histories should be written to

output file(s) for history plotting.
Note: If WRHPV > total pore volume injected or maximum simulation
time, the data will not be printed.
WRPRF - Indicates interval at which information should be written to the output

files controlled by the flags listed on input lines 3.2.3 and 3.2.4.
Note: If WRPRF > total pore volume injected or maximum simulation
time, the data will not be written.
RSTC - Indicates the interval at which restart data should be written to the
RESTAR output file.
3.7.25 DT (This line is read only if IMES=1 and ITSTEP=0)
DT - Time step size for constant time step option.

Units: days
3.7.26 DT, DCLIM, DTMAX, DTMIN (This line is read only if IMES=2 and
ITSTEP=0)

Units: days
DCLIM - Tolerance for concentration change for the first three components.
121
UTCHEM User's Guide
Units: days

Units: days
3.7.27 DT, DCLIM, CNMAX, CNMIN (This line is read only if IMES=2, ITSTEP=1,
and at least one well is not shut-in)

Units: days
DCLIM - Tolerance for concentration changes for the first three components.


3.7.28 DT, (DELC(KC), for KC=1, N), DTMAX, DTMIN (This line is read only if
IMES=3 or 4 and ITSTEP=0)
DT - Initial time-step size.

Units: days
DELC(KC) - Tolerance for concentration change, ∆Clim,κ, of KCth component

(IMES=3) or relative tolerance for concentration change of KCt h
component (IMES=4).
input line 3.7.7.a)

Units: days

Units: days
3.7.29 DT, (DELC(KC), for KC=1, N), CNMAX, CNMIN (This line is read only if
IMES=3 or 4, ITSTEP=1, and reservoir is not shut-in)

Units: days

(IMES=4).
input line 3.7.7.a)
122
UTCHEM User's Guide


123
UTCHEM User's Guide
User notes:
124
UTCHEM User's Guide
APPENDIX A: OUTPUT FILES
A OUTPUT FILES
The following sections describe: (A.1) data that is written to the profile data
file, (A.2) restart run procedure, (A.3) data written to stored restart data file, (A.4)
data written to history files for each well, (A.5) data written to history of reservoir
properties and overall injection and production rates from all the wells, (A.6) data
written to aqueous phase tracer concentration data files, and (A.7) output files
created when running GMS.
A.1 Data Written to Profile Data File
The information in the following lists is written to the profile data file
(PROF).
Printed for each well at each CUMHI1 interval:

Time, number of time steps
Time step size
Courant number
Cumulative pore volume injected
Original in place for each component
Cumulative injection for each component
Cumulative production for each component
Amount retained for each component
Relative error for each component
Fraction of oil recovered
If IREACT ≥ 2:
Average number of iterations, computation time
For each well:
Position of the well, first and last well block completed
Cumulative injection/production
Bottomhole pressure for each well block
All well related information (such as pressure for each phase, phase
concentration, phase cut, etc.)
Producer wellbore temperature and phase cut and concentration
Printed for each gridblock at each CUMPR1 interval:

Phase flux in each direction if IVEL=1
Permeability reduction factor if IRKF=1
Surfactant, polymer, or organic adsorbed concentration if IADS=1
Index of surfactant phase environment if IPHSE=1
125
UTCHEM User's Guide
A.2 Restart Run Procedure
The restart procedure is available with UTCHEM. This enables a user to
continue a run past the initial time period or to break a large run up into smaller
segments. Each time you run UTCHEM, a file called RESTAR is created. This
file (described in section 6.3) contains all the information necessary to continue the
run at a later time. In order to do so, the user needs to:
1.- Rename the output file RESTAR from the previous run to INPUT2
2.- Set the variable IMODE equal to 2 on input line 3.1.3 of file INPUT
3.- Change the value of TMAX on input line 3.3.1 of file INPUT to the new
injection period being simulated in the restart run
4.- Change the value of TINJ on input line 3.7.8 of file INPUT if appropriate
5.- Add additional information for input lines 3.7.14 through 3.7.29 of file
INPUT if the well conditions are different for the new injection period
Note: Make sure the source code you run the restart problem (IMODE=2) has the
same values for the array sizes in the parameter statement as the one used in
original run (IMODE=1).
126
UTCHEM User's Guide
A.3 Data Written to Stored Restart Run Data File
The information in the following list is always written to the stored restart
data file (RESTAR). If the user is running a RESTART run, this data file needs to
be renamed to correspond to the INPUT2 input file. The values in parentheses are
the FORTRAN variable names as they appear in the code.
Printed at the end of each run:

Values in the PARAM.INC file (see the table in Section2.2 for a
comprehensive list)
Time (T), injection time (TINJ), time step size (DT), number of time steps
(ICNT)
New slug injection or restart flag (IINJ), number of time step reduction
(INEC), cumulative pore volume injection (CUMPV), number of
blocks in X-direction minus 1 (NXM1)
Cumulative injection (CUMI), cumulative production (CUMP), original in
place (OIP) for each component
Cumulative injection/production (CUMQI and CUMQP) for each well
Phase concentration (C), phase saturation (S), effective salinity (CSE),
overall concentration (CTOT), number of phases (NPHASE)
If ICOORD=2:
Boundary concentration (CE), boundary pressure (PE)
Viscosity (VIS), relative permeability (RPERM), injection rate (QI), total
rate for each well (QT), phase rate (Q), bottomhole pressure (PWF)
Pressure (P)
If IADSO>1:
Organic adsorption (C2ADSS)
If LMO=1:
Multiple organic adsorption (C2ADSK)
Surfactant adsorption (C3ADSS), surfactant adsorption parameter (A3DS),
polymer adsorption (C4ADSS)
Permeability reduction factor (RKF), calcium concentration (C6JO),
calcium adsorbed by clay (C6ADSS), calcium adsorbed by
surfactant (C6HATS)
Alcohol 7 partitioning coefficient (X7OLD), alcohol 8 partitioning
coefficient (X8OLD)
Oil breakthrough (BTO), tracer breakthrough (TBT), tracer injection
concentration (CINJT), tracer retardation factor (TRD)
Lower effective salinity (CSEL), upper effective salinity (CSEU)
Density (DEN), capillary pressure (PRC)
Total surfactant (TSURF)
If IPERM=2 and IHYST=1:
Minimum water saturation (SWMIN)
If IPERM=3:
Minimum water saturation (SWMIN), last time step water saturation
(SWOLD), water saturation at reversal point (SWREV), last
time step capillary pressure (PCOLD), capillary pressure at
reversal point (PCREV)
Flag to define the scanning path (IPATH), saturation reversal flag
(ISWREV)
If ICAP=1:
127
UTCHEM User's Guide
Dendiritic concentration (CD), dendiritic saturation (SD), flowing

saturation (SF), total flowing concentration (CTF)
If NG≠0:
Chromium adsorption (C14ADS), gel adsorption (C15ADS), cation
exchange capacity of clay (QW)
If IENG=1:
Cum. heat inj. (CUMHI), cum. heat prod. (CUMHP), temperature
(TEM), total volumetric heat capacity (TVHC)
If IENG=1 and IHLOS=1:
Cum. heat loss (TQLOS), integral for overburden and underburden
heatloss calculations (RINO, RINU), time of change of overburden
temp. from the reservoir block (TTCHG), overburden temperature
(TEMPOB), underburden temperature (TEMPUB)
If IENG =1 and ICOORD =2:
Boundary enthalpy (ENTHE)
If IREACT >1:
Solid concentration (CSLDT), adsorbed concentration (CSORBT),
species concentration (CAQSP), surf. associated cation
concentration (CACATT), cation concentration (CACAT)
Cumulative no. of iteration for geochem option (ITCUM)
If IBIO =1:
Concentrations of attached biomass and intra-biomass
concentrations of chemical species (CB)
Concentrations of reducing power in aqueous phase biomass (RED)
and attached biomass (REDB)
128
UTCHEM User's Guide
A.4 Data Written to Well History Plotting Data File(s)
The information in the following list is always written to the well history
plotting data files (HIST01-HISTxx) for each production well.
Printed at each WRHPV interval:

Cumulative pore volume, time [days], cumulative production [ft3, m3, or
STB], water oil ratio, cumulative oil recovery, total production rate
[ft3/day, m3/day, or STB/day]
If IGAS =1: Water cut, oil cut, microemulsion cut, gas cut
Wellbore pressure for each well block [psi or kPa]
If IENG =1: Wellbore temperature [˚F or ˚C]
For I=1, N:
If ICF(I) =1: phase concentration for component N (C(I,L),
L=1,MXP), total concentration of component N (CTOT(I))
If IREACT>1 or ICF(3)=1: Lower effective salinity, upper effective
salinity, effective salinity
If IBIO =1 and IBPR =1:
Aqueous phase concentrations of components 1 through 8 and oil
components in units of mg/L
If IBIO=1, IBPR=1, and there are cometabolic biodegradation reactions in
which reducing power limitations are considered:
Concentrations of reducing power (NADH) in aqueous phase and
attached biomass
If IREACT >1:
Independent species concentration (CAQSP(KK), KK=1, NIAQ)
[mole/liter of water]
If IRSPS >0: Dependent species concentration (CAQSP(KK),
KK=NIAQ+1,NFLD) [mole/liter of water]
If IREACT =3: Phase concentration of injected + generated
surfactant (PSURF(I,L), L=1,3), total concentration of
injected + generated surfactant (TSURF)
If NSLD >0: Concentration of solid components (CSLDT(KK),
KK=1, NSLD) [mole/liter of pore volume]
If ICNM >0: Log 1 0 of interfacial tension between
water/microemulsion and oil/microemulsion (XIFT1,
XIFT2) [dyne/cm]
The information in the following list is always written to the well history plotting
data files for each injection well.

Cumulative pore volume, time in days, cumulative injection [ft3,
m3, or STB], injection rate [ft3/day, m3/day, or STB/day]
Wellbore pressure for each well block [psi or kPa]
Pressure drop between the wells (for the specific case of one injector
and one producer only) or pressure drop between the
pressure tabs (when NG>0, NY=1, NZ=1, see input line
3.4.168) [psi or kPa]
129
UTCHEM User's Guide
A.5 Data Written to Overall History Plotting Data File
The information in the following list is always written to the overall history
plotting data file (OVERAL).

Cumulative pore volume, time in days, volumetric averaged reservoir
pressure (psi or kPa), cumulative oil produced (%OOIP),
cumulative oil produced (bbls or m3), volumetric averaged reservoir
temperature (˚F or ˚C) (only if IENG=1)
Total injection rate (B/D or m3/day), total production rate (B/D or m3/day),
total fluid injected (1000 bbls or m3), total fluid produced (1000
bbls or m3 ) [Note: The fluid injected and produced values are
calculated for the last time step before the print interval.]
Overall production rate for each phase (QBAR(L) for L=1, MXP where
MXP=3 if IGAS=0 and MXP=4 if IGAS =1) (B/D or m3/day)
Average cut for each phase (FBAR(L) for L=1, MXP where MXP=3 if
IGAS=0 and MXP=4 if IGAS =1)
Average saturation for each phase (SBAR(L) for L=1, MXP where MXP=3
if IGAS=0 and MXP=4 if IGAS =1)
If ICF(3)=1: Cumulative surfactant injected (bbls or m3), Cumulative
surfactant produced (bbls or m3), adsorbed surfactant (bbls or m3),
retained surfactant (bbls or m3), adsorbed surfactant (ml/ml PV)
If ICF(4)=1: Cumulative polymer injected (wt%), Cumulative polymer
produced (wt%), adsorbed polymer (wt%), retained polymer
(wt%), adsorbed polymer (wt. % / PV)
130
UTCHEM User's Guide
A.6 Data Written to Tracer Concentration Observation Point Data File(s)
The information in the following list is written to the tracer observation
history plotting data files (TRAC01-TRACn) for each tracer (if IPOBS=0).

Time in days, cumulative pore volume
For NTW tracers:
Aqueous phase tracer concentration at NOBS observation locations
For NTA tracers:
Gas phase tracer concentration at NOBS observation locations
131
UTCHEM User's Guide
A.7 Output Files Created When Running GMS
The following binary output files are created when running GMS. The
"file" prefix for the file names comes from the prefix of the input file name. For
example, if the input file is named run1.dat, all GMS output file names will start
with "run1" in place of "file" in the table below. For the component concentration
files, spname(#) refers to the value of SPNAME on input line 3.1.24.
Output file name Contents

file_spname(1).dat Total concentration of component 1, created if IPCTOT=1
and ICF(1)≠0 and IPRFLG(1)≠1
ile_spname(2).dat Total concentration of component 2, created if IPCTOT=1
and ICF(2)≠0 and IPRFLG(2)≠1
… …
ile_spname(xx).dat Total concentration of last component, created if IPCTOT=1
and ICF(xx)≠0 and IPRFLG(xx)≠1
file_p001.dat Water pressure data, created if IPRESS=1
file_p002.dat Oil pressure data, created if IPRESS=1
file_p003.dat Microemulsion pressure data, created if IPRESS=1
file_p004.dat Gas pressure data, created if IPRESS=1 and IGAS>0
file_s001.dat Water saturation data, created if IPSAT=1
file_s002.dat Oil saturation data, created if IPSAT=1
file_s003.dat Microemulsion saturation data, created if IPSAT=1
file_s004.dat Gas saturation data, created if IPSAT=1 and IGAS>0
file_t001.dat Temperature data, created if IENG=1 and IPTEMP=1
132
UTCHEM User's Guide
APPENDIX B: GEOCHEMISTRY OPTION (IREACT>1)
B GEOCHEMISTRY OPTION (IREACT>1)
This section gives an example list of elements and reactive species for the
geochemistry options of IREACT=2 or IREACT=4.
Elements or pseudo-element: Hydrogen (reactive), Sodium, Calcium,

Magnesium, Carbonate, Aluminum,
Silicon, Oxygen, Chlorine, S (Injected
surfactant)
Independent aqueous or oleic species: H+, Na+,Ca2+, Mg2+, Al3+, CO2- , Cl-,
3
S-, H4SiO4, H2O
Dependent aqueous or oleic species: Ca(OH)+, Mg(OH)+, Al(OH)2-,Al(OH)2-,
Ca(HCO3)+, Mg(HCO3)+, OH-, HCO-3 ,
H3SiO4-, H2SiO42-, HSi2O63-, Si2O52-,
Al(OH)4-, H2CO3, CaCO3o , MgCO3o
Solid species: CaCO3 (Calcite), Al2Si2O5(OH)4
(Kaolinite), MgCO3 (Magnesite),
NaAlSi2O6.H2O (Analcite), SiO2 (Silica),
Mg(OH)2 (Magnesium Hydroxide)
__ __ __ __
Adsorbed cations on rock surface: H + , Na + , Ca 2+, Mg 2+
= = =
Adsorbed cations on micelles: Na +, Ca 2+, Mg 2+
Aqueous reactions Equilibrium constant

Keq + -
eq
K = [H+] [OH-]
H 2O → 1
← H + OH
1
eq
+ 2- K 2
H + CO 3 → ← HCO3
- eq
K2 =
[ HCO3]
-
[ H+] [ CO2-3]
Keq eq [Ca(OH)+] [ H+]
Ca2+ + H2O → 3 + + K3 =
← Ca(OH) + H [Ca2+]
2+ Keq + + eq [ Mg(OH)+] [ H+]
Mg + H2 O → 4
← Mg(OH) + H K4 =
[Mg2+]
Keq eq [Al(OH)2+][H+]
Al3+ + H2O → 5 2+
← Al(OH) + H
+ K5 =
[Al3+]
Keq eq [Al(OH)2+][H+]2
Al3+ + 2H2O → 6 +
← Al(OH)2 + 2H
+ K6 =
[Al3+]
133
UTCHEM User's Guide
Aqueous reactions (cont.) Equilibrium constant (cont.)

Keq eq [Al(OH)4-][H+]4
Al + 4H2O →
3+ 7 -
← Al(OH)4 + 4H
+ K7 =
[Al3+]
Keq + - eq [H+][H3SiO4-]
H4SiO4 → 8
← H + H3SiO4
K8 =
[H4SiO4]
Keq + 2- eq [ H +]2[H2SiO42-]
H4SiO4 → 9
← 2H + H2SiO4
K9 =
[H4SiO4]
eq
Ca2+ + H+ + CO2-
K 10
→
3 ← Ca(HCO3)
+ eq
K10 =
[ Ca(HCO3)+]
[ Ca2+] [ CO2-3][H+]
eq
Mg2+ + 2- K 11
+ H + CO 3 →
← Mg(HCO3)
+ eq
K11 =
[ Mg(HCO3)+]
[ Mg2+] [ CO23-][H+]
eq
K 12 eq [H2CO3]
2H+ + CO2- →
3 ← H2CO3 K12 =
[ CO2-3][H+]2
eq
K 13 + 3- eq [H+]3[HSi2O63-]
2H 4SiO4 →
← 2H2O + 3H + HSi2O6
K13 =
[H4SiO4]2
eq
K 14 + 2- eq [H+]2[Si2O52-]
2H 4SiO4 →
← 2H + 3H2O + Si2O5
K14 =
[H4SiO4]2
Solid species Solubility product

CaCO3
Ksp 2+
[ ]
2-
1 = [Ca ] CO 3
MgCO3
2 = [Mg ] [CO 3 ]
Ksp 2+ 2-
SiO2 Ksp
3 = [H4SiO4]
Al2Si2O5(OH)4 Ksp + -6 3+ 2
4 = [H ] [Al ] [H4SiO4]
2
NaAlSi2O6.H2O sp + -4 + 3+
K 5 = [ H ] [ Na ][Al ][H4SiO4]
2
Mg(OH)2 Ksp 2+ + -2
6 = [Mg ] [H ]
134
UTCHEM User's Guide
Exchange reactions (on matrix) Exchange equilibrium constant

ex __ 2+
__ +
2Na + Ca 2+ K 1
→
← 2Na + + __
Ca 2+ ex
K1 =
[ ]
Ca [ Na+]2
__ + 2
[ ]
2+
[Ca ] Na
ex __
__ +
2Na + Mg 2+ K 2
→
← 2Na
__ 2+
+ + Mg ex
K2 =
[ ]
Mg2+ [ Na+]2
__ + 2
[ ]
[Mg2+] Na
Kex _+
_ +
H + Na+ + OH- →
←
__ +
3 Na + H O
2
ex [ [ ]
Na+] H
K 3 = __ +
[ ]
Na [ H+]
Exchange reactions (on micelle) Exchange equilibrium constant

K exm
= +
2Na + Ca2+ →
←
1 + = 2+
2Na + Ca exm
K 1 =
[ ]
= 2+
Ca [ Na+]2
[ ]
=+ 2
Na [ Ca2+]
exm exm
where K 1 = β 1 [S-]
K exm
= +
2Na + Mg2+ →
←
2
2Na
+
+
= 2+
Mg exm
K 2 =
[ ]
= 2+
Mg [ Na+]2
[ ]
=+ 2
Na [ Mg2+]
exm exm
where K 2 = β 2 [S-]
135
UTCHEM User's Guide
User notes:
136
UTCHEM User's Guide
APPENDIX C: MAIN PROGRAM FLOW OUTLINE
C MAIN PROGRAM FLOW OUTLINE
The following outline represents the basic flow through the main program
of UTCHEM (that is, the order in which major subroutines are called).
Ι. INOUT
A. FILE1
B. PRINTS
C. PRINTI
D. RESREAD
E. INNAME (called only if NO>1 and IOD is not equal to 1)
F . MOPPST (called only if NO>1)
1. SLV2EQ
G . PRINTO (called only if NO>1)
H . WMEACN (called only if NO>1 and IOD is not equal to 1)
I. DPREAD (called only if ICAP=2 and NTW>0)
J . GEOREAD (called only if IREACT>1)
K . BIOREAD (called only if IBIO=1)
L. WELREAD
M. METRIC (called only if IUNIT=1)
N . INIT
1. WELLIX
a. RADIUS
O. FILE3
P . FILE2
Q. GRDFAC
R. GMSINI (called only if IOUTGMS=1)
ΙΙ. RSTART (called only if IMODE=2)

A. NSLUG
B. WELLIX
1. RADIUS
III. TIME0
A. OMOFR (called only if NO>1)
B. DENSTY
C. RPERM0
D. ADS0 (called only if ICWI≥0)
E. VIS0
IV. TRAN1
V. ASIGN1
VI. TRANS (Transmissibilities)
VII. SOLMAT (Pressure Eq.)

A. WELL
B. BUNDRY (called only if IBOUND=1)
C. JCG
VIII. QRATE
137
UTCHEM User's Guide
IX. CONEQ (Conservation Eq.)

A. GEL
X. ADOIL (called only if IADSO>0)
XI. ADSORB
XII. REACTR (called only if IREACT>1)

A. RENAM1
B. TOTALS
C. MANIPL
D. JACUP
1. GAUSS
2. SOLVE
E. RENAM2
XIII. CSECAL
A. ALCPTN (called only if IALC=1)
1. TWOALC
2. CUBIC
B. CSEOD (called only if NO>1 and IOD is not equal to 1)
C. IONCNG
D. PHASC (called only if surfactant is not present and IREACT<1)
1. MAYER (called only if IMASS=2 and ICOR=1)
a. FLUX (called only if IMASS=2 and ICOR=1)
2. NONEQ
3. ODISWO (called only if NO>1)
a. NONEQK
E. PHCOMP (called only if surfactant is present or IREACT>1)
1. MAYER (called only if IMASS=2 and ICOR=1)
a. FLUX (called only if IMASS=2 and ICOR=1)
2. NONEQ
3. ODISWO
4. TIELIN
a. TRY
5. REVISE
6. VGAMMA
7. SINGLE
8. ODISTM (called only if IHAND=0)
9. ODISTM1 (called only if IHAND=1)
XIV. BIOSOLVE (called only if IBIO=1)

A. FLUX
B. PHABIO (called only if IBPP=1)
C. THIRDD
D. INST (called only if IBKIN =3)
E. SDRIV2 (called only if IBKIN = 1 or 2 and running Cray version
of code)
1. F
2. G
F . DDRIV2 (called only if IBKIN = 1 or 2 and running double
precision version of code)
1. F
138
UTCHEM User's Guide
2. G
XV. OMOFR (called only if NO>1)
XVI. TCAP (called only if ICAP=1)
XVII. TDIFFU (called only if ICAP=2)
XVIII. DENSTY
XIX. ASIGN2
XX. ENGBAL (called only if IENG=1)
XXI. LAUWER (called only if IENG=1 and IANAL=1)
XXII. CAPNUM (called only if ITRAP=1)
XXIII. TRAPNO (called only if ITRAP=2 and IGAS=0)
XXIV. TRAP (called only if IPERM<2 and IGAS=0)
XXV. TRAPG (called only if IPERM<2 and IGAS≥1)
XXVI. HYSWWG (called only if IPERM=2 and ICF(3)=0)
XXVII. HYSWWS (called only if IPERM=2 and ICF(3)=1)
XXVIII. HYSTMW (called only if IPERM=3)

A. CUBIC2
XXIX. UTFOAM (called only if IGAS=2)
XXX. VISCOS
XXXI. WELLCK (called only if ICHEK>0)
XXXII. OUTDT1
XXXIII. OUTDT3
XXXIV. OUTDT2
A. PRINTI
B. PRINTS
XXXV. NSLUG
A. WELLIX
1. RADIUS
XXXVI. RSTART
XXXVII. Go to step V (ASIGN1) if not done
XXXVIII. SUMTAB
139
UTCHEM User's Guide
User notes:
140
UTCHEM User's Guide
APPENDIX D: PHASES AND SPECIES IN UTCHEM
D PHASES AND SPECIES IN UTCHEM
This section gives the component numbering scheme in UTCHEM and the
unit for each component.
The following values for L correspond to the indicated phase:
L Phase
1 Aqueous phase
2 Oleic phase
3 Microemulsion phase
4 Gas phase
The following indices correspond to the indicated components [corresponding

concentration units are listed in square brackets]:
For all values of IREACT:
Index Component [conc. units]

1 Water [volume fraction]
2 Oil [volume fraction]
3 Surfactant [volume fraction]
4 Polymer or silicate (KGOPT=3) [weight percent]
5 Total nonsorbing anions concentration, assumed to all be
chloride anions [meq/ml water]
6 Divalent cations, assumed to all be calcium for IREACT<2
[meq/ml water]
7 Alcohol 1 [volume fraction]
8 Alcohol 2 or Gas [volume fraction]
Organic species (IBIO=0 and NO>1):

9 First organic species [volume fraction]
8+NO Last organic species [volume fraction]
Tracers (NT>0):

9+NO First tracer [dependent on user input]
8+NO+NTW+NTA Last tracer [dependent on user input]
Geochemistry option species (IREACT=2 or IRECAT=3):

9+NO+NTW+NTA First geochemistry component
[meq/ml]
8+NO+NTW+NTA+NGC Last geochemistry component
[meq/ml]
Gel model species (IREACT=1 and NG>0):
141
UTCHEM User's Guide
APPENDIX D: PHASES AND SPECIES IN UTCHEM
KGOPT=1 KGOPT=2 KGOPT=3

Index Component Component Component
[conc. units] [conc. units] [conc. units]
NG1* Na2Cr2O7 [ppm] — —
NG2 CSN2H4 [ppm] Malonate ion [ppm] —
NG3 Cr3+ [ppm] Cr3+ [ppm] OH- [ppm]
NG4 Gel [ppm] Gel [ppm] Gel [ppm]
NG5 Hydrogen [meq/ml] Hydrogen [meq/ml] —
*where NG1 = 9+NO+NTW+NTA
Geochemistry/Gel option species (IREACT=4 and NG>0):

9+NO+NTW+NTA First geochemistry component
[meq/ml]
8+NO+NTW+NTA+NGC Last geochemistry component
[meq/ml]
KGOPT=1 KGOPT=2
Index Component [conc. units] Component [conc. units]
NG1* Na2Cr2O7 [ppm] —
NG2 CSN2H4 [ppm] Malonate ion [ppm]
NG3 Cr3+ [ppm] Cr3+ [ppm]
NG4 Gel [ppm] Gel [ppm]
*where NG1 = 9+NO+NTW+NTA+NGC
Biological model species (IBIO=1):
Index Component
[conc. units]
9+NO+NTW+NTA+NGC+NG First biological species
[mg/l]
8+NO+NTW+NTA+NGC+NG+NOTH Last biological species
[mg/l]
142
UTCHEM User's Guide
APPENDIX E: TIME-STEP SELECTION
E TIME-STEP SELECTION
This section discusses the automatic time-step selection options available in
UTCHEM: (E.1) selector based on method of relative changes for the first three
components, (E.2) selector based on method of relative changes for all the
components, and (E.3) selector based on changes in dimensionless concentration
for all the components.
The Courant number, C, is defined as:
Q ∆t
C =
∆x ∆y ∆z φ
where Q is maximum injection/production per wellblock.
E.1 Method of Relative Changes for First Three Components (IMES=2)
Minimum and maximum time step in days (option ITSTEP=0):
The time step selection is based on the method of relative changes for the
first three components (water, oil, and surfactant) as:
 
∆C
∆t n+1
= ∆t min NBL lim
n  κ = 1, 2, 3 (E.1)
 
 max ∆C i,κ 
 i=1 
where ∆tn+1 is limited to: ∆tmin≤∆tn+1≤∆tmax. ∆tmin and ∆tmax are the input time
step values and ∆Clim is the input value for concentration tolerance (DCLIM).
Minimum and maximum time step in days (option ITSTEP=1):
The minimum and maximum time steps in days are computed based on the
minimum and maximum Courant number as:
CNMIN
∆t min =
nwell nwbc Qi 
min  max 
M=1  i=1 ∆x i ∆y i ∆z i φ i 
and
143
UTCHEM User's Guide
CNMAX
∆t max =
nwell nwbc Qi 
min  max 
M=1  i=1 ∆x i ∆y i ∆z i φ i 
where Courant number (C) is limited to: CNMIN≤C≤CNMAX. CNMIN and
CNMAX are the input Courant number values in input line 3.7.11. The time step is
then calculated from Eq. E.1.
144
UTCHEM User's Guide
E.2 Method of Relative Changes for All Components (IMES=3)
For IMES=3, the method of relative changes is applied to all the
components in the simulation run:
 
∆C
∆t n+1
= ∆t min NBL lim,κ
n 
 κ = 1, ..., nc
 max ∆C i,κ 
 i=1 
where ∆tn+1 is limited to: ∆tmin≤∆tn+1≤∆tmax. ∆Clim,κ is a fraction of the initial or
injected concentration (whichever is larger) of the KCth component, DELC(KC) in

input lines 3.7.12 or 3.7.13. For example: ∆Clim,3 = 0.1 × C3 where C3 is the
total concentration of component 3. If ∆Clim,κ of the KCth component is entered as
zero, that component is not considered in the time-step size selection.
145
UTCHEM User's Guide
E.3 Method of Relative Changes Using Dimensionless Concentration for

All Components (IMES=4)
For IMES=4, the method of relative changes is applied to all the
components in the simulation run:
 
 
R lim,κ
∆t n+1
= ∆t min
n  κ = 1, ..., nc
 NBL ∆C 
 max i,κ

 i=1 C 
 i,κ 
where ∆tn+1 is limited to: ∆tmin≤∆tn+1≤∆tmax. Rlim,κ is the input dimensionless

relative change in concentration of component κ (DELC(KC)) on input lines 3.7.12
or 3.7.13. For example Rlim,3=0.1 indicates a 10% change in concentration of
Component 3.
146
UTCHEM User's Guide
APPENDIX F: CHANGES TO USER'S GUIDE
F CHANGES TO USER'S GUIDE
This section is designed to point out the major modifications to the
UTCHEM user's guide (as opposed to the code) that occurred between versions
7.0 and 9.0. This list is not comprehensive.
F.1 New Variables
This section lists new input variables introduced to UTCHEM with version
7.0. The variables are also flagged within the body of the user's guide by being
printed in italicized boldface.
Variable(s) Line number(s)

IOUTGMS 3.2.1
IPBIO 3.2.3
IHYSTP, IFOAMP, INONEQ 3.2.4
ICWI 3.3.18
CWI 3.3.32-3.3.34
REFK 3.4.130
C5INI 3.4.133
F.2 Modified Variables
This section lists deleted input variables and variables whose definitions
were modified or expanded in version 9.0.
Variable(s) Line number(s) Modification

RUNNO 3.1.1 Restricted size and
expanded usage
IPTRAC, IPBR 3.2.3 (in 7.0 guide) Variable deleted
IPHP 3.2.4 (in 7.0 guide) Variable deleted
IADS, IVEL, IRFK, IPHSE 3.2.5 Moved to new input line
147
UTCHEM User's Guide
APPENDIX F: CHANGES TO USER'S GUIDE
F.3 Other Changes
This section lists additional changes made to the user's guide for version
9.0.
Change Section number (s)

Updated Input and Output file discussion to 2.1
include HEAD file
Changed "Source Code Array Dimensions" 2.2
section to "Description of HEAD File" - major
overhaul of this section
Updated compilation instructions 2.3 and 2.4
Removed Cray compilation information Previously Section 2.5
Updated information A.1
Updated main program flow outline C
Moved work.job description 2.4, removed old Appendix F
148
UTCHEM User's Guide
INPUT VARIABLE INDEX
A14D 85 BTMIN 98 CSLDI 95 EL 32

A15D 85 BVI 60 CSORBI 95 ELCRG 91
ACATK 94 C 113, 120 CSTAR 87 ELEMNT 91
ACATSP 92 C160 85 CUMHI1 114, 121 ENDOG 100
ACIDIS 95 C1BL 108, 117 CUMPR1 114, 121 ENDOGB 100
AD21 75, 76 C1BR 108, 117 CVSPL 86 ENTP 68
AD22 75, 76 C1I 95 CVSPO 86 EPC 70, 71, 72
AD31 77 C2I 95 CVSPR 86 EPC0 70
AD32 77 C2PLC 43 CVSPU 86 EPCO 71
AD41 77 C2PRC 43 CWI 34, 35 EPCW 71
AD42 77 C50 40 D 73 EPHI3 64
ADSBIO 100 C5BL 108, 117 D111 32 EPHI4 64
AG1 85 C5BR 108, 117 DAMX 99 EPSALC 46
AG2 85 C5I 94 DCLIM 115, 121, 122 EPSBIO 97
AGK 85 C5INI 78 DCS20 48, 49 EPSME 43
AHUH 50 C60 40 DELC 115, 116, 122 EPXLO 87
AK1 83 C6BL 108, 117 DEN1 64, 65 EQB0 51
AK1T 84 C6BR 108, 117 DEN2 64, 65 EQB1 51
AK2 83 C6I 94 DEN23 64 EQB2 51
AK2T 84 CAC2I 94 DEN3 64, 65 EQB3 51
AK7 46 CAQI 94 DEN7 64, 65 EQK 93
AK8 46 CBI 100 DEN8 64, 65 EQW 77
AKA 101 CBIOMN 100 DENBIO 100 EQWPS 95
AKM7 46 CELAQI 94 DENBLK 97 ER 92
AKM8 46 CHACAT 94 DENO 86 EXACAT 94
AKN 104 CHARGE 93 DENS 74, 86 EXCAI 94
AKR 105 CHUH 50 DENU 86 EXEX 93
AKS 101 CI 94 DMC 82, 83 EXK 93
AKWC7 46 CINIT 99 DMX 82 EXSLD 93
AKWC8 46 CLAMDA 68 DMZ 83 FACTX 37, 38, 39, 40
AKWS7 46 CM 79 DNOILC 66 FACTY 30
AKWS8 46 CMIN 97 DNOME 66 FACTZ 31
ALPDS 68 CNEM 51 DR 92 FADS 77
ALPHAL 73, 74 CNMAX 115, 116, 122 DT 115, 116, 121, 122 FEA 101
ALPHAT 73, 74 CNMIN 115, 116, 122 DTMAX 115, 121, 122 FFH 79
ALPHAV 63 COI 36 DTMIN 115, 121, 122 FFL 79
ALPWS 68 COLNUM 100 DX 20 FLDSPS 91
AP1 63 COMPC 67 DX1 17, 18, 19 FN 104
AP2 63 COMPR 27 DXL 81 FPR 103, 104
AP3 63 CPC 70, 71, 72 DY 20 FRC 105
AR 92 CPC0 69 DY1 17, 18, 19 FRP 105
B14D 85 CPCO 71 DYL 81 FVF 67
B15D 85 CPCW 71 DZ 20 G11 50
B3D 77 CRG 85 DZ1 17, 18, 19 G12 50
B4D 77 CRK 64 DZL 81 G13 50
BB 93 CRNAK 85 E13C 59 G21 50
BETA6 45, 50 CRTC 86 E1W 57, 58 G22 50
BETA7 45 CS0 48, 49 E1WC 57 G23 50
BETA8 45 CS1 48, 49 E23C 59 GAMHF 63
BETAP 63 CS2 48, 49 E2W 57, 58 GAMMAC 63
BGK 85 CSE1 63 E2WC 57 HBNC70 43
BIORME 98 CSEL7 45 E31C 59 HBNC71 43
BR 92 CSEL8 45 E3W 58 HBNC72 43
BRK 64 CSELP 89 E3WC 57 HBNC80 44
BRMAX 101 CSET 44 E4C 62 HBNC81 44
BRMAXB 101 CSEU7 45 E4W 62 HBNC82 44
BSIHB 103 CSEU8 45 E4WC 62 HBNS70 43
BTMAX 98 CSEUP 89 EFMIN 99 HBNS71 43
149
UTCHEM User's Guide
HBNS72 43 IEW 53 IPR 103, 104 NCOMP 47, 49

HBNS80 44 IFACT 37, 38, 39, 40 IPRESS 31 NCOMPS 101
HBNS81 44 IFGHBN 43 IPRF 109, 118 NCR 91
HBNS82 44 IFIRST 109, 118 IPRFLG 23 NCRFD 91
HBNT0 44 IFLAG 109, 118 IPRW 53 NELET 90
HBNT1 44 IFOAMP 25 IPSAT 23 NEX 90
HBNT2 44 IFT 50 IPTEMP 23 NFLD 90
HINIT 33 IGAS 13 IREACT 13 NG 20
HNAK 85 IGROW 105 IREV 74 NGC 20
HPHAD 89 IHAND 43 IRKF 26 NH 90
IADK 77 IHB 103 IRLIM 105 NHFD 91
IADS 26 IHLOS 86 IRO 108, 118 NIAQ 90
IADS1 77 IHYST 53 IRSPS 89 NIHB 101
IADSO 74 IHYSTP 25 ISGI 33 NMET 98
IALC 45 II1 18, 19, 81 ISOL 51 NMG 90
IANAL 86 II2 18, 19, 81 ISRW 53 NMOD0 37
IBIO 13 IKADS 74 ISTB 67 NMOD1 38
IBKIN 98 ILAST 109, 118 ISTOP 23 NMOD2 38
IBL 107, 117 IMASS 51 ISUB 80, 101, 102, 103, NMOD3 39
IBMOD 116 IMAX 37, 38, 39, 40 104, 105 NMOD4 39
IBOUND 107, 116 IMES 13 ISWI 31 NMOD5 40
IBPP 98 IMIN 37, 38, 39, 40 ITC 13 NNA 90
IBR 107, 117 IMKX 37 ITRAP 52 NNUT 101
IBS 101, 102, 103, 104, IMKY 37 ITREAC 13 NO 20
105 IMKZ 37 ITRU 21 NOBS 26
IBT 107, 117 IMOD 27 ITSTEP 108, 118 NOTH 20
IBTIM 97 IMODE 13 ITYPE 99 NOXY 90
ICAP 13 IMPOR 37 IUNIT 15 NPABIO 99
ICF 21 IMSW 37 IVEL 26 NPROD 101
ICHEK 112, 119 IMTVAR 99 IVIS 25 NRT 80
ICHRGE 90 INAME 66 IW 109, 118 NSILI 90
ICKL 25 INONEQ 25 IZONE 107 NSLD 90
ICNM 25 INUT 104 JJ1 18, 19, 81 NSLEL 90
ICOI 36 IOBS 26 JJ2 18, 19, 81 NSORB 90
ICOMET 101 IOD 47 JMAX 37, 38, 39, 40 NSORBX 92
ICOMPO 67 IODEN 65 JMIN 37, 38, 39, 40 NSUB 80
ICOORD 13 IOUTGMS 23 JOBS 26 NSURF1 90
ICOR 51 IOVIS 59 JW 109, 118 NTA 20
ICPC 68 IOW 68 KC 99, 100 NTW 20
ICREX 85 IP1 85 KGOPT 83 NWEL1 118
ICSE 25 IP2 85 KK1 18, 19, 81 NWEL2 120
ICSUB 102 IPALK 23 KK2 18, 19, 81 NWELL 108
ICUMTM 23 IPARK 53 KKX 82 NWREL 108
ICWI 31 IPBIO 23 KKZ 82 NX 15
ICWM 13 IPCAP 23 KMAX 37, 38, 39, 40 NY 15
ID 113, 114, 118, 120, IPCTOT 23 KMIN 37, 38, 39, 40 NZ 15
121 IPER 25 KOBS 26 OCOMP 48, 49
IDEN 64, 65 IPERM 52 KPH 82 OEACN 66
IDEPTH 31 IPERMX 27 KPRF 111, 119 OMWT 66
IDIR 109, 118 IPERMY 27 MSUB 80 ONAME 48, 49, 66
IDISPC 13 IPERMZ 27 N 20 OPSK7O 45
IDW 109 IPGEL 23 NACAT 90 OPSK7S 45
IDWW 120 IPHAD 89 NACD 91 OPSK8O 45
IDXYZ 15 IPHSE 26 NALMAX 46 OPSK8S 45
IEA 101, 102, 103, 104, IPOBS 23 NALU 90 OVIS 60
105 IPOLYM 64 NBC 98 P 33
IENG 13 IPOR1 27 NCA 90 P1RC 59
IEPC 68 IPPRES 23 NCARB 90 P1RW 56
150
UTCHEM User's Guide
P1RWC 56 S2RWC 54 WRHPV 114, 121

P2RC 59 S2RWC4 61 WRPRF 114, 121
P2RW 56, 57 S3RC 58 WSOL 51
P2RWC 56 S3RW 55 WSOLK 51
P3RC 59 S3RWC 54 WTABL 107, 117
P3RW 56, 57 S4RC 62 WTABR 107, 117
P3RWC 56 S4RW 61 X13 83
P4RC 62 S4RWC 61 X14 83
P4RW 61, 62 SCHARG 93 X16 83
P4RWC 61 SCR 83 X4 83
PBL 108, 117 SCS0 48, 49 XCORD 16
PBR 108, 117 SCS1 48, 49 XIFTG 62
PERMX 29 SCS2 48, 49 XIFTGW 62
PERMXC 29 SGI 34 XIFTW 50
PERMY 30 SHRTN 87 XKC 77
PERMYC 30 SLDSPS 91 XKS 77
PERMZ 30 SORBSP 91 XL1 81
PERMZC 30 SOSTAR 87 YL1 81
PHC 89 SPK 94 YXS 101
PHIC 82 SPNAME 21 ZCORD 16
PHT 89 SSLOPE 63 ZL1 81
PHT1 89 SSTAR 71, 72, 73 ZLSUB 82
PINIT 33 SWELL 109, 118
PNEG 68 SWI 33
POR 28, 29 SWSTAR 87
PORC1 28 T11 52
PORCM 82 T22 52
POWN 63 T33 52
PRESS1 33 T44 62
PSTAND 27 TAK 80
PT7 46 TAKT 80
PT8 46 TC 105
PWF 113, 120 TCOL 100
PWFMAX 112, 119 TCONO 86
PWFMIN 112, 119 TCONU 86
QI 113, 114, 120, 121 TDEN 80
QTMAX 112, 119 TEMINJ 113, 120
QTMIN 112, 119 TEMPI 86
QV 77 THETAX 32
R 17, 19, 20 THETAY 32
RABIO 99 TINJ 114, 121
RCOL 100 TITLE 13
RDC 78 TK 78
REDI 105 TKS 78
REDUC 93 TKT 78
REFK 77 TMAX 27
RET 78 TMW 80
RFMAX 87 TSTAND 59
RSTC 114, 121 VELGR 87
RUNNO 13 VFRACM 82
RW 109, 118 VIS1 59
S 33, 34 VIS2 59
S1RC 58 VIS4 60
S1RW 54, 55 VNDS 68
S1RWC 54 VNWS 68
S2RC 58 VSLOPG 60
S2RW 55 WELNAM 111, 119
S2RW4 61 WM4 83
151

UTCHEM Users Guide PDF

Uploaded by

Copyright:

Available Formats

UTCHEM Users Guide PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

UTCHEM Users Guide PDF

Uploaded by

Copyright:

Available Formats

What is UTCHEM used for?

What is UTCHEM used for?

What physical phenomena does UTCHEM model?

What physical phenomena does UTCHEM model?

Volume I:

Reservoir Engineering Research Program

INPUT VARIABLE INDEX 149

• dispersion • organic biodegradation capability

See Section 2 of the UTCHEM Technical Documentation for the general

2.0 OPERATION OF THE SIMULATOR

2.1 Input and Output Files

Input File Name Contents

Output file Contents

2.2 Description of HEAD file

Example HEAD file:

variables used in the HEAD file.

2.3 Compilation and Execution on Workstations

The UTCHEM distribution package for workstations contains the

execute the code:

1. Issue the command:

2. To build an executable file called utchem9.exe, issue the command:

make -f Makev9 FC=f90

Example work.job file:

Command line Description

2.4 Compilation and Execution on a WINDOWS-based PC

The UTCHEM distribution package for PCs contains the FORTRAN-90

(ex21.data and ex21.head). At The University of Texas at Austin, we

probably compile successfully using other compilers with little or no modification if

below to compile UTCHEM using the Digital FORTRAN compiler:

1. Launch Visual Studio.

4. Select "Build utchem9.exe" from the Build menu. This will

7. Wait patiently while the code compiles. Upon completion, exit

8. Copy the input and head files to the directory containing

9. Double click on utchem9.exe. A screen will appear and will

3.0 INPUT DATA DESCRIPTION

The following is a list of variables as they are read by UTCHEM. The

italicized boldface as well). Every list of variables is followed by a description

pertinent to the problem being run.

3.1. Title and Reservoir Description Data

that each data line is preceded by three comment lines.

RUNNO - Run number.

TITLE - Title and run description.

IMODE - Flag indicating if the problem to be run is a first run or a restart

IMES - Flag indicating if a constant or automatic time-step is to be used.

IDISPC - Flag indicating which type of numerical dispersion control is used.

Note: These methods are applied to both concentration and relative

ICWM - Flag indicating if the concentration well model is used or not.

ICAP - Flag indicating if the capacitance model is used or not.

IBIO - Flag indicating whether or not biodegradation reactions occur

ICOORD - Flag indicating which coordinate system is used.

ITREAC - Flag indicating if a tracer reaction is used or not.

ITC - Flag indicating if second-order time approximation is used or not.

Note: We recommend that second-order time approximation (ITC=1) only

IGAS - Flag indicating if gas phase is considered or not.

IENG - Flag indicating if temperature variation is considered or not.

3.1.4 NX, NY, NZ, IDXYZ, IUNIT

NX - Number of gridblocks along X-axis (ICOORD=1, 3, or 4) or number of

NY - Number of gridblocks along Y-axis.

NZ - Number of gridblocks along Z-axis.