TABLE OF CONTENTS

1. INTRODUCTION............................................................................................................................... 1

1.1 INSTALLING AND RUNNING S3D................................................................................................... 1

1.2 PROGRAM UPDATES ..................................................................................................................... 2

2. INPUT AND OUTPUT FILES .......................................................................................................... 2

2.1 FILE NAME DEFINITION................................................................................................................ 2

2.2 PARAMETERS ............................................................................................................................... 4
2.3 DETAILS OF INPUT DATA.............................................................................................................. 4
2.3.1 Well and Control Parameters (first input file, data01)........................................................... 4
2.3.2 Gridblock and Rock Properties Data (second input file, data02)......................................... 10
2.3.3 Streamline Data (third input file, data03) ............................................................................ 12
2.3.4 Fluid Property Data (fourth input file, data04).................................................................... 14
2.3.5 Infill drilling data (fifth input file, data05) ........................................................................... 16
2.4 DETAILS OF OUTPUT DATA AND FILES .............................................................................................. 17

3. TECHNICAL REFERENCES......................................................................................................... 20

4. DISCLAIMER AND COPYRIGHT INFORMATION................................................................. 20

 1

1. Introduction
This is a user manual for S3D Ver 3.2, a three-dimensional streamline simulator. It can
be used to model two-phase waterflooding and also single/two-phase tracer transport
under incompressibility assumptions. The source code (not supplied) is written in
FORTRAN 77 / FORTRAN 90 and runs under any operating system. However, the
executables supplied with the GUI is compiled for Windows. The streamline approach
comprises of two steps: generating streamlines in 3D space and then solving the 1D
transport equations analytically or numerically along the streamlines. The following are
the important features of S3D.

1. Waterflood/tracer simulation (Maximum grid size = 100x100x30, Maximum

number of producers = 35, Maximum number of injectors = 35).
2. Pressure updating.
3. Analytic/Numerical solution for 1D solution.
4. Infill drilling /changing well configuration.

This simulator is meant for educational/research purposes only and is not intended for
commercial use.

1.1 Installing and running S3D

Follow the steps outlined below for installation of S3D and getting started.

• Double click on Set-up.exe.

• Follow the steps to install the program.
• It is recommended to install the program with a simple path name, for example,
C:/S3D.
• The installation process with generate an ‘icon’ for the program.
• Double click on the ‘icon’ (or use the Start menu) to start the program.
 2

• After starting the program, follow these steps: File -> Existing Project. Now
browse to the directory ‘dataFiles’ and inside the directory you will see several
directories related to the computer exercises.
• For testing purposes, enter the directory ‘Homogeneous-5Spot’. Select the file
with extension, *.s3d and click open.
• Now you should be able to view all the data related to this case under the menu
‘Input Data’.
• You can now run the program (Run ->Start) and view the results in Plots.

1.2 Program Updates

This program is under continuous development at Texas A&M University. Updates to the
program, if any, will be notified at www.pe.tamu.edu/datta-gupta.

2. Input and Output Files

2.1 File Name Definition

The names of the input and output files should be specified in the file, “stream.dat”.
Sample input and output files are supplied. There are five input files and five (or more)
output files. A sample “stream.dat” file may look like the following.

data01.test
data02.test
data03.test
data04.test
data05.test
welltracer.out
gridpress.out
intermed.out
ratelayer.out
n_of_failsline.out

The first five lines in “stream.dat” are the names of the input files. The first file (e.g.
“data01.test” as above) corresponds to unit 01, second file to unit 02 and so forth for all
the five input files. The next five files are the various output files. A description of the
files are given in the following table.
 3

Note : the file extension can be other than ‘.test’ also, e.g. data01.dat or data01.txt.

File name Contents of files

stream.dat File name definition
data01.test Well data/solution control parameters (input file)
data02.test Gridblock/rock properties data (input file)
data03.test Streamline/tracer injection data (input file)
data04.test Fluid property data (input file)
data05.test Infill drilling/recurrent data (input file)
welltracer.out Well wise tracer response output (output file)
gridpress.out Pressure solution (output file)
intermed.out Intermediate calculation results when run is terminated (output file)
ratelayer.out Production and injection rate information for each layer (output file)
n_of_failsline.out Streamline tracing failure information (output file)

Some more output files are listed below.

File name # Contents of files assigned to the unit

S3dblackw.log This log file contains useful diagnostic information about the status of the run. It
is a good idea to check this file after every run. It contains the contents of the
DOS window that appears during execution of the program.
Production.his Water-cut and oil-rate vs. time for all wells
Sat_Num.dat Water saturation distribution at each gridblock for each pressure update and infill
drilling. Can be directly input to Tecplot (Ver. 7 or higher). Plot by “zones”. Each
zone represents each pressure update or infill drilling.
Tau.dat Time of flight distribution for streamlines at each gridblock. Can be directly input
to Tecplot ( Ver 7 or higher)
Plot_stream.dat Configuration of three-dimensional streamlines. Gives the 3-D streamlines co-
ordinates. Can be directly input to Tecplot.
Vsweep.out Pattern and layer wise volumetric sweep calculations from the time of flight
connectivity and also from method of moments using tracer response at
producers.
Stream.dbg Streamline information for debugging purposes
Ratealloc.out Pattern wise rate allocation of injected water at each pressure update and at the
end of the simulation run
Prod.out Writes out the total flow rate after summation over the streamlines and supplied
or calculated flow rate for each production well. The content of the file is mainly
used to compare the calculated/supplied flow rate with rate obtained from
summing up the streamlines for a well.
PoreVolume.out Writes out (for each pressure update) the percentage of fluid injected and
produced as a function of total pore volume of the system.
Recovery.out Gives the fractional recovery at each time steps
 4

2.2 Parameters
The parameters specifying the array size are declared in the parameter statements in the
source code (not supplied) and in the common header file "flocom". Prior to the
compilation and execution of the program, these parameters should be correctly defined
corresponding to the expected array size. This version of the program is set for the
following maximum size restrictions.

Maximum grid size = 100x100x30

Maximum number of producers = 35

Maximum number of injectors = 35

Maximum number of streamlines = 40000

2.3 Details of Input Data

In S3D, all input data should be entered as a combination of a "label or keyword", a
“descriptor” as shown in the sample input data set and then the “data values”. This
combination is called a data group. An example of a data group is given below.
index label or keyword
idxa descriptor
0 data value
The order of input for each data group is flexible. Blank lines can be input between data
groups but not within a group unless specified.

2.3.1 Well and Control Parameters (first input file, data01)

The following data groups describe the input data required in the FORTRAN unit 01
input file (e.g. data01.test).

Label Input variables

index idxa

where:
 5

idxa Control parameter for the use of the analytical solution for tracer transport
with longitudinal dispersivity
=0 analytical solution without longitudinal dispersivity
=1 analytical solution with longitudinal dispersivity. This option will work
only if no pressure update is selected (see below).

Label Input variables

nxyz nx, ny, nz, nphase

where:
nx: Number of gridblocks in the x direction in Cartesian coordinate
ny: Number of gridblocks in the y direction in Cartesian coordinate
nz: Number of gridblocks in the z direction in Cartesian coordinate
nphase: Number of phases — note: set nphase = 1 for single phase simulation (e.
g. single phase tracer case); for partioning tracer case and two phase
simulation ( e.g. water flood) set nphase=2.
MaxP: Total number of producers during the entire simulation time, includes
any new producer resulting from infill drilling, e.g. 4 producer initially
and then 1 additional producer after infill drilling, then Maxp=5.

Label Input variables

wprod nwelp
iwp(i), jwp(i), kwpt(i), kwpb(i), rwp(i), rop(i), qp(i), iflagp, pwfp_obs,
NumWell(i) for i=1, nwelp

where:

nwelp: Number of production wells

iwp(i): i th location of the gridblock in the x direction indicating the location of
(i)-th production well
jwp(i): j the location of gridblock in the y direction indicating the location of (i)-
th production well
 6

kwpt(i): k th location of the gridblock in the z direction indicating the top

completion layer of (i)-th production well
kwpb(i): k th location of the gridblock in the z direction indicating the bottom
completion layer of (i)-th production well
rwp(i): Wellbore radius of (i)-th production well (ft)
rop(i): Equivalent radius of (i)-th production well (ft) — note: for example set
rop(i) = 0.208dx for square gridblocks (see details in Peaceman 1978,
1983)
qp(i): Production rate (total) of (i)-th production well (STB/day)
iflagp(i): Control parameter for bottom hole pressure vs. total rate constraint.
= 0 rate constrained.
= 1 pressure constrained. Note: The production rate will be ignored.
pwfp_obs(i): Bottom hole pressure for pressure constrained well.
NumWellP(i): Well number for each producer

Label Input variables

winj Nweli
iwi(i), jwi(i), kwit(i), kwib(i), rwi(i), roi(i) qi(i), iflagi(i), pwfp_obs(i),
NumWell(i) for i=1, nweli

where:
nweli: number of injection wells
iwi(i) : i th location of the gridblock in x direction indicating the location of (i)-
th injection well
jwi(i): j th location of the gridblock in y direction indicating the location of (i)-
th injection well
kwit(i): k th location of the gridblock in z direction indicating the top completion
layer of (i)-th injection well
kwib(i): k th location gridblock in z direction indicating the bottom completion
layer of (i)-th injection well
rwi(i): wellbore radius of (i)-th injection well (ft)
roi(i): equivalent radius of (i)-th injection well (ft) — note: set roi(i) = 0.208dx
for square gridblocks (see details in Peaceman 1978, 1983)
 7

qi(i): injection rate of (i)-th injection well (STB/day)

iflagi(i): Control parameters for bottom hole pressure vs. total rate constraint.
= 0 rate constrained.
= 1 pressure constrained. Note: The injection rate will be ignored.
pwfi_obs(i): Pressure for pressure constrained well.
NumWellI(i): Well number for each injector

Important Note: For incompressible flow as assumed here, the total production
must equal total injection for closed systems. Hence, at least one of the
injection or production wells must be pressure constrained if no flow
boundaries are specified on all sides. Otherwise, the pressure matrix will
become singular (see Aziz and Settari, 1979).

Label Input variables

update i_update
i_pup_sol
update_time_total, n_pup, nfactor, dtaufrac, time_write

where:
i_update: Control parameter for pressure update.
= 0 No pressure update is required.
= 1 pressure update is required.
Note: If you want to use analytic solutions along streamlines, then specify
i_update = 0 and infill = 0. The ‘infill’ keyword is discussed later.
Production history output is controlled with time_write parameter
discussed below.

i_pup_sol: Control parameter for solution scheme for 1-D waterflood solution
along streamlines.
= 0 Analytical solution.
= 1 Numerical solution.
 8

update_time_total: Total time the program going to simulate before starting

infill drilling, if any. Effective for both i_update = 0 and 1.
n_pup: no. of pressure updates in between the starting time and
update_time_total. e.g. if the update_time_total = 400 days and n_pup =
4, then the pressures will be updated four times at intervals of 100 days.
nfactor: It is a multiplier to determine the number of the nodes along streamline
for 1-D numerical solution. For example, if a streamline passes through
20 grid blocks and nfactor = 2, then the discretization along the
streamline for numerical solution uses 40 nodes and saturation
propagations are calculated on these nodes. The nfactor should be
minimum 1 and integer values.
dtaufrac: It is a multiplier that specifies the time step for 1-D numerical solution
along streamlines. In fact it is the same as the CFL number here and
dtaufrac should be < 1 for the stability of the numerical solution.
time_write: The time interval (days) at which the output water-cut data will be
written.

Label Input variables

Infill n_recurrent_total

where:
infill: Specifies infill drilling, new well rates or changes in well-configuration.
n_recurrent_total: Number of infill drilling or rate changes. If no infill drilling or
rate changes are desired, then set infill = 0. If infill =1, then we must
specify i_update=1 as we require pressure updates to account for the
changes.

Label Input variables

tracer Itracer
Ntrc
(part(i, j), i=1, 2), (diffus(i, j), i=1, 3), (disp(i, j), i=1, 2), for j=1, ntrc

where:
 9

itracer: Control parameter for tracer option.

= 0 tracer option off.
= 1 tracer option on.
ntrc: number of tracers
part(1, j): partition coefficient of (j)-th tracer to oil phase (option disabled: keep
default values)
part(2, j): partition coefficient of (j)-th tracer to gas phase (option disabled: keep
default values)
diffus(1, j): molecular diffusion coefficient of (j)-th tracer in water phase (ft2/day)
(option disabled: keep default values)
diffus(2, j): molecular diffusion coefficient of (j)-th tracer in oil phase (ft2/day)
(option disabled: keep default values)
diffus(3, j): molecular diffusion coefficient of (j)th tracer in gas phase (ft2/day)
(option disabled: keep default values)
disp(1, j): longitudinal dispersivity of (j)-th tracer (ft) (you can use this option
provided no pressure update has been selected)
disp(2, j): transverse dispersivity of (j)-th tracer (ft) (option disabled: keep default
values)

Label Input variables

outfmt Idbg, i_plot, i_zone

where:
idbg: Control parameters for debugging option.
= 0 debugging option off.
= 1 debugging option on.
i_plot: Control parameters for streamline file generation.
= 0 streamline file generation off.
= 1 streamline file generation on.
i_zone: Control parameters for writing out saturation map at each pressure
update
=0 only at the last update time
 10

=1 write out at each pr. update

2.3.2 Gridblock and Rock Properties Data (second input file, data02)
The following data groups describe the input data required in the FORTRAN unit 02
input file.

Label Input variables

dx dx(i), for i=1, nx

where:
dx(i): length of (i)-th gridblock in x direction (ft)

Label Input variables

dy dy(i), for i=1, ny

where:
dy(i): length of (i)-th gridblock in y direction (ft)

Label Input variables

dz (((dz(i, j, k), for i=1, nx), for j=1, ny), for k=1, nz)

where:
dz(i, j, k) : thickness of (i, j, k)-th gridblock in z direction (ft)

Label Input variables

depth (((depth(i, j, k), for i=1, nx), for j=1, ny), for k=1, nz)

where:
depth(i, j, k) : top depth of (i, j, k)-th gridblock (ft)
 11

Label Input variables

por (((por(i, j, k), for i=1, nx), for j=1, ny), for k=1, nz)

where:
por(i, j, k): porosity at (i, j, k)-th gridblock

Label Input variables

permx (((permx(i, j, k), for i=1, nx), for j=1, ny), for k=1, nz)

where:
permx(i, j, k): permeability in x direction at (i, j, k)-th gridblock (md)

Label Input variables

permy (((permy(i, j, k), for i=1, nx), for j=1, ny), for k=1, nz)

where:
permy(i, j, k) : permeability in y direction at (i, j, k)-th gridblock (md)
Note: permy can also be represented as a fraction of permx. Using ky=kx*1.0
will result in ky=kx.
Label Input variables
permz (((permz(i, j, k), for i=1, nx), for j=1, ny), for k=1, nz)

where:
permz(i, j, k): permeability in z direction at (i, j, k)-th gridblock (md)
Note: permz can also be represented as a fraction of permx. Using kz=kx*1.0
will result in kz=kx.

All the data for porosity and permeability can also be input from an external file
and the user has to specify the file name after each label index.
 12

Label Input variables

ngr (((rng(i, j, k), for i=1, nx), for j=1, ny), for k=1, nz)

where:
rng(i, j, k): net to gross ratio at (i, j, k)-th gridblock.

Label Input variables

rock (((irock(i, j, k), for i=1, nx), for j=1, ny), for k=1, nz)

where:
irock(i, j, k): rock type for (i, j, k)-th gridblock. If more than one rock type is
used, different relative permeability curves could also be used (in data04).

All the data for ngr and rock type can also be input from an external file and the
user has to specify the file name after each label index.

Label Input variables

sat (((sat(l, i, j, k), for l=1, nphase-1), for i=1, nx), for j=1, ny), for k=1, nz

where:
sat(l, i, j, k): saturation of (l)-th phase at (i, j, k)th gridblock
note: l=1 : water, l=2 : oil
Note: For oil-water, keep nphase=2 for two-phase case. Here the user
can input either the uniform or non-uniform initial saturation. Non-
uniform initial saturation can be input from an external file and the
file name can be mentioned after the label index.

2.3.3 Streamline Data (third input file, data03)

The following data groups describe the input data required in the FORTRAN unit 03
input file.
 13

Label Input variables

title Title

where:
title : title on the 3D streamline model run (up to 80 characters)

Label Input variables

nline nlinei

where:
nlinei: expected total number of streamlines to be traced — note: the number of
streamlines actually traced by the model will be close to nlinei, but not
necessarily equal to nlinei, because streamlines are drawn in equal areal
increments at the faces of gridblocks of production wells

Label Input variables

tcutoff tcutoff

where:
tcutoff: the threshold time-flight-cutoff value for calculating volumetric sweep
efficiency

Label Input variables

Cinj Nstepc
(((cinj(i, j, k), for i=1, nstepc), j=1, ntrc), for k=1, nweli)

where:
nstepc : number of time steps for input of injection tracer concentration
cinj(i, j, k): injection concentration of (j)-th tracer at (i)th time step in (k)th
injection well — note: if nstepc is less than nstep, cinj(nstepc, j, k) is
 14

regarded as injection concentration after the (nstepc)-th time step;

hence in the case of constant concentration injection, only one time
step of cinj(i, j, k) may be input

2.3.4 Fluid Property Data (fourth input file, data04)

The following data groups describe the input data required in the FORTRAN unit 04
input file.

Label Input variables

Fluid vis,fvf,fvfw,viso

where:
vis : viscosity of water (cp)
fvf : formation volume factor of oil (bbl/STB)
fvfw : formation volume factor of water (bbl/STB)
viso: viscosity of oil (cp)

Label Input variables

Press datum, Pdatum, woc

where:
datum : datum depth (ft)
pdatum: average reservoir pressure (psia).
woc : depth of oil water contact (ft)
Note: The pressure of the reservoir is adjusted based on the pressure constrained
well(s). Recall that at least one production well has to be pressure
constrained to solve the pressure solution for closed systems.

Label Input variables

 15

Sides ibl, pbl

ibr, pbr
ibn, pbn
ibs, pbs
ibt, pbt
ibb, pbb

where:
sides: Boundary condition of the reservoir. Six sides of the reservoir can be
“closed” or “open”.
ibl: Control parameter for the boundary condition of the left side of the
reservoir. “l” stands for left.
=0 closed boundary (left side).
=1 open boundary (left side).
pbl: Value of pressure at the open boundary. Only considered if left boundary
is open i.e. ibl = 1
Note: Similarly for other sides. The abbreviations are as follows:
r: right
n: north
s: south
t: top
b: bottom

Label Input variables

Twoph itwo_ph itable
end_kro, end_krw, sor, swc, real_no, real_nw
or
itwo_ph itable
ntab, jval
TSw, TKrw, Tkro

itwo_ph: Control parameter for two-phase oil-water flow

=1 oil-water flow.
 16

=0 only tracer flow.

itable flag to select the input of relative permeability from table or end
points
=0 From end points
=1 From table
end_kro relative permeability of oil at connate water saturation (swc)
end_krw relative permeability of water at residual oil saturation (sor)
sor residual oil saturation
swc connate water saturation
real_no saturation exponent for oil.
real_nw saturation exponent for water.
ntab Number of saturation table
jval Number of saturation values in the table
TSw,TKrw,Tkro Water Saturation, Relative Perm to water, Relative Perm to
oil, repeat this table ntab times

Label Input variables

density denw, deno

denw : density of water, lbm/ft^3.

deno : density of oil, lbm/ft^3.

2.3.5 Infill drilling data (fifth input file, data05)

Note: Give new production well configuration only if “infill” is = 1. Specify all
the wells that will be active. Refer to the first input file for explanation of the
variables/keywords.

Label Input variables

wprod nwelp
iwp(i), jwp(i), kwpt(i), kwpb(i), rwp(i), rop(i) qp(i), iflagp(i),
pwfp_obs(I), NumWellP(i) for i=1, nwelp
 17

Label Input variables

winj Nweli
iwi(i), jwi(I), kwit(i), kwib(i), rwi(i), roi(i) qi(i), iflagi(i), pwfi_obs(i),
NumWellI(i) for i=1, nweli

Label Input variables

update i_update
i_pup_sol
update_time_total, n_pup, nfactor, dtaufrac, time_write
i_plot_stream

Note: See first input file for explanation of the keywords. The
“update_time_total” is the cumulative time elapsed from the beginning of
the simulation. The “update_time_total” has to be higher than the
“update_time_total” given in first input file.

Note: Repeat the above four labels for the number of infill drilling. For example,
if infill=2 is specified in first input file, then repeat the four labels twice,
one for each new well configurations.

2.4 Details of output data and files

Various file names and definitions for the output files are listed in the section 2.1. In
this section a brief overview of these output files and ways to visualize them are
discussed.

1. Production.his: This output file gives the output of water cut and oil
production data as a function of time for each well. First column is the time
step, second column is the water cut data and third column is the oil
production data for the producer number 1, then fourth and fifth columns give
 18

the water cut and oil production data respectively for the producer number 2
and so on for the other producers. This file can be directly taken to the
EXCEL SPREAD SHEET to make a plot of the water cut and oil rate against
time in days for all producing wells.
2. Tau.dat: This output file gives the time of flight at each grid block. This file
can be taken directly to TECPLOT to make a plot of the time-of-flight
contours. An example time-of-flight contour plot has been shown in the figure
below.

Figure 1-Time-of-Flight distribution for a heterogeneous Quarter 5 spot

3. Plot_stream.dat: This is the output file for plotting the streamlines in 2D or

3D with color codes showing the injector-producer relations and injection rate
allocations (see ratealloc.out for numerical values). This output file can be
directly taken to TECPLOT for plotting the streamlines in 3D or 2D. An
example of the streamline plot for a field case with arbitrary well pattern has
been shown in the next figure.
 19

I2 P7
I13 I14
I4
P23 P19 P24
P20 P15 P3 P12 P6
I11 I12 I9 I7
P17 P22 P18 P14 P2 P5
P21 P11
I5 I10 I8
P27 P9 P10 P1 P16 P4 P13
I15 I6
I1 I3
P25 P26 P8

Figure 2- Streamline plot(X-Y view) for a 3D field case with arbitrary well pattern
4. Saturation.dat: This output file gives the water saturation distribution for
each grid blocks at the final time step for the semi-analytical option. This
output can directly be taken to the TECPLOT for plotting the saturation
distribution.

Figure 3- Saturation distribution for a 3-D field case at certain time

5. Sat_Num.dat: This output file also gives the saturation distribution when
pressure updates are selected and numerical calculations are performed along
streamlines. For other detail see section 2.1. This output can be directly taken
to TECPLOT for plotting the saturation distribution as shown in the above
figure.
 20

6. Ratealloc.out: This output file gives the rate allocation for each injector to
each producer. The output is given at each pressure update time and also for
each infill drilling.
7. Vsweep.out: This output file gives the volumetric sweep calculations for each
layer for a particular threshold value of time-of-flight.
8. Recovery.out: This output file gives the recovery (fraction of original oil in
place recovered) at each time step.

3. Technical References
1. Datta-Gupta, A. and King, M.J.: “A Semianalytical Approach to Tracer Flow
Modeling in Heterogeneous Permeable Media,” Advances in Water Resources, 18(1),
9, pp189-24, (1995).

2. King, M. J. and Datta-Gupta, A., "Streamline Simulation: A Current Perspective," In

Situ, 22 (1), (1998).

3. Datta-Gupta, A., “Streamline Simulation: A Technology Update,” SPE Distinguished

Author Series, Journal of Petroleum Technology, 68-73 (December 2000).

4. Cheng, H., Osako, I., Datta-Gupta, A. and King, M. J. “A Rigorous Compressible

Streamline Formulation for Two and Three Phase Black Oil Simulation,” SPE
Journal, 11(4), December 2006, p407-417.

4. Disclaimer and Copyright information

The program is written for research purposes and is distributed WITHOUT ANY
WARRANTY. The authors do not warrant that the operation of the program will be
uninterrupted or error free. The authors are not liable for any cost, damages or any other
consequences arising from the use of this program or guarantee that the program will
serve the intended purpose. No portion of the information and documents provided may
be reproduced in any form without prior permission. This program is protected by
international copyright laws and treaties.
 21

Appendix : Sample Datafiles

Data01

index
idxa
0

nxyz
nx ny nz nphase
41 41 1 2

MaxP
Maximum producing wells(Vertical including all infills)
8

wprod
nwelp
4
iwp jwp kwptop kwpbot rwp rop qp i_constraint pressure
wellID
1 21 1 1 0.25 1.04 3.56 1 2500. 1
21 1 1 1 0.25 1.04 3.56 1 2500. 2
41 21 1 1 0.25 1.04 3.56 1 2500. 3
21 41 1 1 0.25 1.04 3.56 1 2500. 4

winj
nweli
1
iwi jwi kwitop kwibot rwi roi qi i_constraint Pressure wellID
21 21 1 1 0.25 1.04 800.00 0 0. 1

update
i_update
1
i_pup_sol
1
update_time_total n_pup n_factor dtaufrac time_write
500.0 5 2 0.25 20.0

infill
N_recurrent_total
1

tracer
itracer
0
ntrc
1
part diffus disp
0.0 0.0 0.0 0.0 0.0 0.0 0.0

outfmt
dbg
1 1 1
 22

Data02

dx
41*52.86

dy
41*52.86

dz
1681*37.

depth
1681*4000.0

por
1681*0.03

permx
perm_high.out

permy
ky=kx*1.0

permz
kz=kx*1.0

ngr
1681*1.0

rock
1
1681*1

sat
1681*0.2
 23

Data03

title
test run for a 2d reservoir with infill drilling

nline
number of streamlines
200

tcutoff
Swept Volume Calculation Time
500

cinj
number of time steps for input of injection tracer
2
1*1.0 0.0
 24

Data04

fluid
Visw Bo Bw Viso
1.0 1.2 1.0 0.80

press
datum dt.Pressure woc
4000. 4100. 4100.

sides
index pressure at boundary
left
0 4000.
right
0 4000.
north
0 4000.
south
0 4000.
top
0 0.
bottom
0 0.

twoph
itwo_phase itable
1 0
end_kro end_krw sor swc real_no real_nw
0.5 0.8 0.25 0.20 2.0 2.0

density
denw deno
62.5 50.0
 25

Data05

wprod
nwelp
8
iwp jwp kwptop kwpbot rwp rop qp i_constraint pressure
WellID
1 21 1 1 0.25 1.04 3.56 1 2500. 1
21 1 1 1 0.25 1.04 3.56 1 2500. 2
41 21 1 1 0.25 1.04 3.56 1 2500. 3
21 41 1 1 0.25 1.04 3.56 1 2500. 4
1 1 1 1 0.25 1.04 3.56 1 2500. 5
41 1 1 1 0.25 1.04 3.56 1 2500. 6
41 41 1 1 0.25 1.04 3.56 1 2500. 7
1 41 1 1 0.25 1.04 3.56 1 2500. 8

winj
nweli
1
iwi jwi kwitop kwibot rwi roi qi i_constraint Pressure WellID
21 21 1 1 0.25 1.04 600.00 0 0. 1

update
i_update
1 /0->no pr_update 1->pr_update
i_pup_sol
1 /0->analytical 1->numerical
update_time_total n_pup n_factor dtaufrac time_write
1000.0 5 2 0.25 20.0
i_plot_stream
1