化学驱数值模仿软件
CDNS（4.0）
利用阐明
戴 家 林
2008年1月
目 录
一、组分设计 1
二、化学离子设计 2
三、CDNS软件运转流程 3
四、数据流的填写办法 4
五、要害字细致形貌 7
1、TTT 7
2、DATE 7
3、TMAX 7
4、UNIT 8
5、INJECT 8
6、GRIDMD 9
7、GRIDVL 15
8、WELLVL 19
9、PARA 21
10、A－DATA 28
11、S-DATA 30
12、P－DATA 35
13、T－DATA 41
14、SLUG 44
15、DT 46
16、ENDFW 47
17、PERFOR 48
18、TIME 48
19、STOP 49
六、数据流的例子 49
七、后果输入文件 57
1、数据反省信息文件 57
2、目标剖析信息文件 58
3、曲线绘制数据文件 58
4、等值图和三维可视图绘制数据文件 59
5、流线图绘制数据文件 59
化学驱数值模仿软件CDNS利用阐明
一、组分设计
CDNS可以模仿碱、表活剂、聚合物、示踪剂等化学剂参加注入液的驱油历程。软件设计了下列物质组分：
1、水
2、油
3、表活剂
4、聚合物
5、阴离子
6、二阶阳离子
7、钠
8、氢
9、碳酸根（拟组分）
10、煤油酸（拟组分）
11、第1种示踪剂
12、第2种示踪剂
13、第3种示踪剂
依据注入化学剂的状况和原油的性子来确定组分数。
1、假如注入液中没有任何化学剂，即为纯水驱，组分数为2（1－2）。
2、假如注入液中没有碱，也没有示踪剂，组分数为6（1－6）。
3、假如注入液中没有碱，有示踪剂，组分数为6加示踪剂的个数。
4、假如注入碱，原油中没有煤油酸，组分数为9加示踪剂的个数。
5、假如注入碱且原油中有煤油酸，组分数为10加示踪剂个数。
二、化学离子设计
当注入液中参加碱当前，油藏流体中将会有一系列的化学反响产生，这些反响包罗水相中的化学反响，固体的沉淀消融反响、岩石外表的离子互换反响和表活剂上的离子互换反响。为了形貌这些化学反响，CDNS设计了下列化学离子成份。
1、水相中的化学离子
⑴ H+
⑵ Na+
⑶ Ca2+ (Mg2+)
⑷ CO32-
⑸ HAo（油里的煤油酸）
⑹ H2O
（以上6种是独立的离子）
⑺ HCO3-
⑻ OH-
⑼ A-
⑽ HAw（水里的煤油酸）
2、固体化学离子
⑴ CaCO3
⑵ Ca (OH)2
3、在岩石外表上吸附的离子
⑴
⑵()
三、CDNS软件运转流程
1、CDNS共有两个步伐：
①CDNS1
②CDNS2
CDNS1是初始化步伐，CDNS2是主模子模仿盘算步伐。
CDNS1必要一个用户预备的输出数据文件（数据流）。CDNS2必要的输出数据从CDNS1运转后的输入文件失掉。
2、CDNS运转流程
四、数据流的填写办法
1、要害字
CDNS的所有输出数据部由要害字举行引导。共有18个一级要害字，它们是：
（1）TTT
（2）DATE
（3）TMAX
（4）UNIT
（5）INJECT
（6）GRIDMD
（7）GRIDVL
（8）WELLVL
（9）PARA
（10）A－DATA
（11）S－DATA
（12）P－DATA
（13）T－DATA
（14）SLUG
（15）DT
（16）ENDFW
（17）PERFOR
（18）TIME
（19）STOP
2、要害字阐明
（1）要害字引导输出的各种数据均接纳自在F款式。二数之间至多相距一个空格。
（2）若统一要害字在数据流文件呈现屡次，新的将取代老的。
（3）当某些要害字在数据流文件中缺失时，步伐主动对响应的数据置缺省值。
3、要害字的序次
CDNS对要害字之间的序次没有限定，但从模仿盘算的必要动身，一样平常序次该当是：
TTT（＊）
DATE
TMAX
UNIT（＊）
INJECT
GRIDMD
GRIDVL或WELLVL
PARA
A－DATA（＊）
S－DATA（＊）
P－DATA（＊）
T－DATA（＊）
SLUG
DT
ENDFW (*)
PERFOR（＊）
TIME
STOP
此中有（＊）的要害字是可选择的。
五、要害字细致形貌
1、TTT
款式：
TTT
ttt
阐明：
ttt是题目内容，它是一个不凌驾50个字符的字符串。
例：
TTT
Gudong 52＋3
2、DATE
款式：
DATE nn yy dd
阐明：
开端模仿的日期或油田投产日期
例：
DATE 1998 07 11
3、TMAX
款式：
TMAX tmax ooip
阐明：
tmax－这次模仿的最大天数
ooip－原始地质储量
单元：104m3
例：
TMAX 5000 1538
4、UNIT
当接纳英制单元时，才利用此要害字，缺省单元是公制单元
款式：
UNIT k
阐明：
k――k＝0 公制单元
k＝1 英制单元
例：
UNIT 0
5、INJECT
款式：
INJECT ka ks kp kgel kfoam kt kgas kacid
阐明
ka——ka＝0 不注碱
ka＝1 注碱,
ks——ks＝0 不注表活剂
ks＝1 注表活剂
kp——kp＝0 不注聚合物
kp＝1 注聚合物
kgel——kgel＝0 不注GEL组分
kgel＝1 注GEL组分
kfoam——kfoam＝0 不注FOAM组分
kfoam＝1 注FOAM组分
kt——kt＝0 不注示踪剂
kt＝1 注示踪剂
kgas——kgas＝0 无气相
kgas＝1 有气相
kacid——kacid＝0 无煤油酸
kacid＝1 有煤油酸
例：
INJECT 1 1 1 0 0 0 0 1
6、GRIDMD
款式：
GRIDMD k1 k2
阐明：
k1――坐标体系标识
k1＝1 直角坐标体系
k1＝2 柱坐标体系
k1＝3 变宽度剖面坐标体系
k2――分别网格办法标识
k2＝1 人工分别网格
k2＝2 步伐主动分别网格
例：
GRIDMD 1 1
6.1假如k2＝1，在GRIDMD行前面有下列数据行：
（1）数据行1
款式：
nx ny nz nw nwa
阐明：
nx――x偏向的网格块数
ny――y偏向的网格块数
nx――z偏向的网格块数
nw――消费井数（包罗注入井和采出井）
nwa――材料井数。假如nwa＝0则所有消费井也是材料井。
例：
10 10 3 2 0
（2）数据行2
款式：
dxi （i=1,nx）
阐明：
dxi－x偏向第i个网格块的步长，
单元：m或feet。
例：(nx=10)
10＊50.0
(3)数据行3
款式
dyi (i＝1，ny)
阐明：
dyi－y偏向第i个网格块的步长，
单元：m或feet。
例：(ny=10)
3＊50 20 30 5＊40.0
(4)消费井信息行，每口井一行，共nw行。
款式：
wname iw jw flag rw swell
idir ifirst ilast pwfmin pwfmax qmin qmax
阐明：
wname－井的名字，不凌驾10个字符
iw－井位在x偏向上的网格号
jw－井位在y偏向上的网格号
iflag－井别标识
iflag＝1，定产注入井
iflag＝2，定压消费井
iflag＝3，定压注入井
iflag＝4，定发生产井
rw－井筒的半径
单元：m或feet。
swell－井的表皮因子
idir－井筒偏向标识
idir=1，沿x偏向的程度井
idir=2，沿y偏向的程度井
idir=3，垂直井
ifirst－井射开层段的第一个网格块号
filast－井射开层段的最初一个网格块号
pwfmin－最小井底压力限定
单元：Kpa 或pai
pwfmax－最大井底压力限定
单元：Kpa 或pai
qwin－最小速率限定
单元：m3/day或ft3/day
qmax－最大速率限定
单元：m3/day或ft3/day
例：
A－15 5 3 1 0.15 0 3 1 3 5000 15000 0 500
(5)材料井信息行。假如nwa＞0，则有nwa个材料井信息行，每口材料井一行，共nwa行。
款式：
wname iw jw
阐明
wname－井名，不凌驾10个字符
iw－材料井在x偏向的网格块号
jw－材料井在y偏向的网格块号
例：B－301 8 17
6.2假如k2=2，在GRIDMD行前面有下列数据行：
（1）数据行1
款式
xoo yoo nz nw nwa nbou
阐明：
xoo yoo－坐标系原点的坐标
单元：m或feet。
nz－Z偏向上的网格数（层数）
nw－消费井数（包罗注入井和采出井）
nwa－材料井数。若nwa为0，则所有消费井是材料井
nbou－界限线的点数
例：
0 0 5 11 0 20
(2)消费井信息行，每口井一行，共NW行。
款式：
wname xw yw iflag rw swell idir dfirst dlast pwfmin pwfmax qmin qmax
阐明：
wname－井名，不凌驾10个字符
xw,yw－井在x－y立体上的坐标
单元：m或feet。
dfirst－井射开段顶部的深度
单元：m或feet。
dlast－井射开段底部的深度
单元：m或feet。
别的与5.1中消费井信息行的阐明相反。
例：
A－487 70 25 4 0.15 0 3 1250 1380 8000 16000 0 400
(3)材料井信息行
假如nwa＞0，每口材料井一行，共nwa行。
款式：
wname xw yw
阐明：
wname－井名
xw,yw－井的坐标
单元：m或feet。
例：
GD－33 84 35
(4)界限点信息行
若nbou＞0，则有nbou个界限点信息行，共nbou行
款式：
xb yb
阐明
xb,yb－界限点的坐标
单元：m或feet。
界限点按逆时针偏向排序，第一个界限点与第nbou个界限点重合。
例：
15 3
7、GRIDVL
这个要害字用于输出地质属性在网格块上的散布值。
GRIDVL用双要害字办法输出数据。双要害字的第一个要害字是地质属性的名字，它可以是下列名字中的任何一个：
DEPTH（深度）
DZ (网格厚度)
DZNET(净厚度)
POR（孔隙度）
PERMX（x偏向上的浸透率）
PERMY（y偏向上的浸透率）
PERMZ（z偏向上的浸透率）
PW（初始水相压力）
SW（初始水饱和度）
双要害字的第二个要害字是赋值办法，它可以是下列要害字之一：
CC XX YY ZZ VV L1 DIP MLT MOD
阐明:
（1）CC
标明后面的数组网格值是一个常数。
例：
POR CC 0.33
（2）xx
标明数组值只在x偏向上变革，应输出nx个值。
例：
nx＝10
PERMX xx
1500 1400 1300 1260 1250 1200 1000 950 900 850
（3）yy
标明数组值只在y偏向上变革，应输出ny个值。
例：（ny=5）
PERMX YY
1400 1800 1900 2000 2100
（4）ZZ
标明数组值只在Z偏向上变革，应输出nz个值。
例：（nz=3）
DZNET ZZ
4.0 5.5 3.2
(5)VV
标明输出所有网格块上的值
例：（nx=10,ny=5,nz=2）
DEPTH VV
1500 1510 1520 1530 1540 1550 1560 1570 1580 1590
1510 1520 1530 1540 1550 1560 1570 1580 1590 1600
1520 1530 1540 1550 1560 1570 1580 1590 1600 1610
1530 1540 1550 1560 1570 1580 1590 1600 1610 1620
1540 1550 1560 1570 1580 1590 1600 1610 1620 1630
1505 1515 1525 1535 1545 1555 1565 1575 1585 1595
1515 1525 1535 1545 1555 1565 1575 1585 1595 1605
1525 1535 1545 1555 1565 1575 1585 1595 1605 1615
1535 1545 1555 1565 1575 1585 1595 1605 1615 1625
1545 1555 1565 1575 1585 1595 1605 1615 1625 1635
(6)L1
标明只输出第一层的网格值，别的层由步伐主动盘算。L1只实用于第一个要害字是DEPTH（深度）。
例：（nx=10 ny=5 nz=3）
DEPTH L1
50*1230
(7)DIP
DIP只用于DEPTH（深度），标明只输第一个网格块的值，别的网格块的深度由步伐主动盘算。
款式：
DEPTH DIP d111 xdd ydd
阐明：
d111－第一个网格块（1,1,1）的深度
单元：m或feet
xdd－x偏向上的倾角
单元：度
ydd－y偏向上的倾角
单元：度
例：（nx＝10，ny＝5，nz=3）
DEPTH DIP 1230 10 7.5
(8)MLT
标明该数组的值是另一个数组的倍数
款式：
数组名1 MLT fact 数组名2
阐明：
fact－倍数因子
例：
PERMY MLT 1.0 PERMX
（9）MOD
MOD用于对数组一局部网格上的值举行修正，别的值不动。
款式：
数组名 MOD i1 i2 j1 j2 k1 k2 # bb
阐明：
i1－x偏向肇始网格号
i2－x偏向完毕网格号
j1－y偏向肇始网格号
j2－y偏向完毕网格号
k1－z偏向肇始网格号
k2－z偏向完毕网格号
＃－算术运算符，它是“＋”、“－”、“＊”、“/”、“＝”之一。
bb－数值
例：PERMX MOD 2 3 3 4 1 1 ＝ 1500
8、WELLVL
外地质属性参数只在井点有值时，必要用WELLVL来输出井点的地质参数，然后步伐主动插值盘算网格块上的值。
款式：
WELLVL intpl
阐明：
intpl－插值办法标识
intpl＝1 改良的间隔加权插值法
intpl＝2 克里金插值法
intpl＝3 分形插值法
例:
WELLVL 1
8.1材料井数据信息行
要害字WELLVL行的前面跟有nwa口材料井的数据信息，每口井有nz行，共有nwa＊nz行。
款式：
jh kz wdepth wdz wpor wkx wky wkz wpw wsw
阐明：
jh－材料井号，不凌驾10个字符
kz－层号，1≤kz≤nz
wdepth－该井层顶深
单元：m或feet
wdz－该井层的厚度
单元：m或feet
wpor－该井层的孔隙度
单元：分数
wkx－该井层x偏向上的浸透率
单元：10-3um2
wky－该井层y偏向上的浸透率
单元：10-3um2
wkz－该井层z偏向上的浸透率
单元：10-3um2
wpw－该井层的初始水相压力
单元：kpa或ps;
wsw－该井层的初始水饱和度
单元：分数
例：
A8 1 1200 3.5 0.35 1300 1300 130 12000 0.45
9、PARA
要害字PARA引导用户输出模仿必要的一些关于油藏的常数，包罗岩石、水、油的物感性质、地层水和注入水的化学成份、绝对浸透率曲线和毛管压力曲线等。
常数用下列要害字输出，这些要害字是可选择的，假如没有输出，响应的常数将取缺省值。
ROCK
WATER
OIL
KRHI
KRLO
PC
FWCHEM
IWCHE1
IWCHE2
9.1 ROCK
用于输出岩石的紧缩系数和油藏的尺度条件
款式：
ROCK compcr pstand tstand
阐明：
compcr－岩石的紧缩系数
单元：1/kpa 或1/psi
pstand－尺度条件下的参照压力
单元：Kpa 或psi；
tstand－尺度条件下的参照温度
单元：℃或℉
例：
ROCK 0.00001 12000 40.0
9.2 WATER
输出水的物感性质
款式：
WATER denw visw conpcw xiftw
阐明：
denw－水的密度
单元：g/cm3或psi/ft
visw－水的粘度
单元：mPa.s
compcw－水的紧缩系数
单元：1/kpa 或1/psi
xiftw――油水界面张力的对数值
例：
WATER 0.99 0.45 0.00005 1.3
9.3 OIL
输出油的物感性质
款式：
OIL deno viso conpco acidv
阐明：
deno－油的密度
单元：g/cm3或psi/ft
viso－水的粘度
单元：mPa.s
compco－油的紧缩系数
单元：1/kpa 或1/psi
acidv－煤油酸的酸值
单元：无量纲
例：
OIL 0.92 70.0 0.0001 1.7
9.4 KRHI
输出高界面张力时的绝对浸透率曲线
款式：
KRHI S1r S2r P1r P2r e1 e2
阐明：
S1r－水相的约束水饱和度
单元：分数
S2r－油相的剩余油饱和度
单元：分数
P1－水相绝对浸透率的最大油
单元：分数
P2－油相绝对浸透率的最大值
单元：分数
e1－水相绝对浸透率曲线的指数
单元：无量纲
e2－油相绝对浸透率曲线的指数
单元：无量纲
例：
KRHI 0.30 0.20 0.2 0.95 1.8 2.3
9.5 KRLO
输出低界面张力时的绝对浸透率曲线。
款式：
KRLO s1rc s2rc p1rc p2rc e1c e2c
阐明：
s1rc——水相的约束水饱和度
单元：分数
s2rc——油相的剩余油饱和度
单元：分数
p1rc——水相绝对浸透率的最大值
单元：分数
p2rc——油相绝对浸透率的最大值
单元：分数
e1c——水相绝对浸透率曲线的指数，无量纲
e2c——油相绝对浸透率曲线的指数，无量纲
例：
KRLO 0. 0. 0.8 1.0 1.1 0.9
9.6 PC
输出毛管压力曲线数据
款式：
PC iow cpc epc
阐明：
iow－岩石的润湿性标识
iow=0 亲水
iow=1 亲油
cpc－毛压曲线端点值
单元：或
epc－毛压曲线指数
单元：无量纲
例：
PC 0 62 2.1
9.7 FWCHEM
输出地层水的化学构成
款式：
FWCHEM C－CL C－H C－Na C－Ca C－Mg C－CO3 C－HCO3
C－H2CO3 C－OH
阐明
C－CL－地层水cl－的浓度
C－H－地层水H＋的浓度
C－Na－地层水Na+的浓度
C－Ca－地层水Ca2+的浓度
C－Mg－地层水Mg2＋的浓度
C－CO3－地层水CO32－的浓度
C－HCO3－地层水HCO3－的浓度
C－H2CO3－地层水H2CO3的浓度
C－OH－地层水OH－的浓度
单元：PPm （即mg/L）
例：
FWCHEM 5600 0.000006 400 80 23 37 220 110 13
9.8 IWCHE1
输出第一种注入水的化学构成
款式：
IWCHE1 C－CL C－H C－Na C－Ca C－Mg C－CO3 C－HCO3
C－H2CO3 C－OH
阐明：
C－CL－注入水cl－的浓度
C－H－注入水H＋的浓度
C－Na－注入水Na+的浓度
C－Ca－注入水Ca2+的浓度
C－Mg－注入水Mg2＋的浓度
C－CO3－注入水CO32－的浓度
C－HCO3－注入水HCO3－的浓度
C－H2CO3－注入水H2CO3的浓度
C－OH－注入水OH－的浓度
单元：ppm（即mg/L）
例：
IWCHE1 4300 0.0000001 230 120 55 41 60 95 18
9.9 IWCHE2
输出第一种注入水的化学构成
款式：
IWCHE2 C－CL C－H C－Na C－Ca C－Mg C－CO3 C－HCO3
C－H2CO3 C－OH
阐明：
C－CL－注入水cl－的浓度
C－H－注入水H＋的浓度
C－Na－注入水Na+的浓度
C－Ca－注入水Ca2+的浓度
C－Mg－注入水Mg2＋的浓度
C－CO3－注入水CO32－的浓度
C－HCO3－注入水HCO3－的浓度
C－H2CO3－注入水H2CO3的浓度
C－OH－注入水OH－的浓度
单元：ppm（即mg/L）
例：
IWCHE2 4300 0.0000001 230 120 55 41 60 95 18
10、A－DATA
要害字A－DATA引导用户输出关于碱的参数。
碱的数据用下列要害字输出：
A－TYPE
A－PH
A－ACID
10.1 A－TYPE
款式：
A－TYPE itype
阐明：
itype——碱的范例标记
itype=1 NaOH
itype=2 NaCO3
10.2 A－PH
款式：
A－PH phc pht
阐明：
phc——PH的临界值，大于该值后表活剂吸附依赖于PH值
pht——PH的临界值，大于该值后表活剂不产生吸附
例：
A－PH 7.0 11.2
10.3 A－ACID
款式：
A－ACID cselp cseup cac2i acidis eqwps
阐明：
cselp——可以发生表活剂的含盐量低限
单元：meq/ml
cseup——可以发生表活剂的含盐量高限
单元：meq/ml
cac2i——油里的初始煤油酸浓度
单元：moles/l 油
acidis——煤油酸的消融系数
eqwps——煤油酸确当量
例：
A－ACID 0.02 0.25 0.003 0.006 500
11、S-DATA
要害字S-DATA引导用户输出关于表活剂的参数。
表活剂的数据用下列要害字输出：
IFTTAB
S－DEN
S－ADSO
S－EXCH
S－CSEL
11.1 IFTTAB
款式：
IFTTAB nift cseift
阐明：
nift——实行室界面张力曲线的数据点数
cseift——测试时含盐量
单元：mg/l
假如nift>0,则前面跟有nift个数据行。
款式：
xift yift zift
阐明：
xift――碱的浓度
单元：wt%
yift――表活剂浓度
单元：wt%
zift――界面张力
单元：mN/m
例：
IFTTAB 6
0 0 20
0 0.11 5
0 0.1 0.5
0 0.2 0.01
0 0.3 0.002
0 0.4 0.001
11.2 S－DEN
款式：
S－DEN den3 compc3
阐明：
den3——表活剂的密度
单元：g/m3或psi/ft
compc3——表活剂的紧缩系数
单元：1/kPa或1/psi
例：
S－DEN 0.99 0.0
11.3 S－ADSO
输出吸附参数
此要害字有两种款式任选其一 。
款式1：
S－ADSO 1 ephi3 nsads csesad
阐明：
ephi3－外表活性剂的可及体积分数
nsads－吸附曲线的数据点数量
csesad－测吸附曲线的含盐量 单元：mg/L
假如nsads＞0则前面跟有nsads个数据行，
款式：
xsads ysads
阐明：
xsads－外表活性剂浓度 单元mg/L
ysads－外表活性剂的吸附量 单元mg/g
例：S－ADSO 1 1.0 15 5000
0.0 0.0
500 0.15
10000.21
15000.25
2000 0.31
3000 0.4
4000 0.52
5000 0.64
60000.75
7000 0.86
8000 0.98
9000 1.02
10000 1.02
11000 1.02
12000 1.02
款式2：
S－ADSO 2 ephi3 ad31 ad32 b3d
阐明：
ephi3 －表活剂的可及体积分数
ad31——表活剂的吸附参数，a31
单元：体积比
ad32——表活剂的吸附参数，a32
单元：ml/meq
b3d——表活剂吸附参数，b3
单元：1/表活剂体积
例：
S－ADSO 2 1.0 0.6 0.25 1010
11.4 S－EXCH
输出离子互换参数
款式：
S－EXCH qv xkc xks eqw
阐明：
qv——粘土的离子互换才能
单元：(meq/ml)/pv
xkc——粘土的阴离子互换常数，c
单元：(meq/ml)-1
xks——表活剂的阴离子互换常数，s
单元：(meq/ml)-1
eqw——表活剂确当量
例：
S－EXCH 0.05 0.25 0.2 419
11.5 S－CSEL
输出含盐量的参数
款式：
S－CSEL csel cseu 6
阐明：
csel——进入Ⅲ相的最低的含盐量
单元：meq/ml
cseu——进入Ⅲ相的最高的含盐量
单元：meq/ml
6——关于钙的含盐量的斜率参数，无量纲
例：
S－CSEL 0.012 0.54 0.8
12、P－DATA
要害字P－DATA引导用户输出有关聚合物的参数。在P－DATA行前面用下列要害字输出聚合物的详细参数：
VISTAB
P－GAMM
P－REDU
P－ADSO
12.1 VISTAB
要害字P－VIST用于输出一种聚合物的粘浓曲线表。
款式:
VISTAB nvis csevis
阐明:
nvis－聚合物粘浓曲线的数据点数
csevis－测试时的含盐量
单元：mg/l
假如nvis＞0，则前面跟有nvis个数据行。
款式：
xvis yvis zvis
阐明：
xvis－碱的浓度
单元：分量百分数
yvis－聚合物的浓度
单元：ppm(即mg/l)
zvis－聚合物溶液的粘度
单元：mpa.s
例：
VISTAB 6 5000
0 0 0.45
0 100 3.0
0 500 8.4
0 1000 12.5
0 1500 16.8
0 2000 25.1
12.2 P－GAMM
要害字 P－GAMM用于输出聚合物的剪切参数。
款式：
P－GAMM gammac gamhf pown
阐明：
gammc－剪切方程中的系数
单元：
gamhf－零剪切速率时的粘度半值时的剪切速率
单元：sec－1
Pown－聚合物粘度随剪切速率变革的指数
单元：无量纲
例：
P－GAMM 3.2 20 1.1
12.3 P－REDU
要害字P－REDU用于输出浸透率低落的参数。
此要害字有两种款式任选其一。
款式1：P－RERU 1 npredu csered
阐明：
npredu－浸透率降落曲线数据点数量
csered－测试浸透率降落曲线的含盐量 单元：mg/L
假如npredu＞0，前面跟有npredu个数据行。
款式：
xpredu ， rk
阐明：
xpredu－聚合物浓度 单元：wt%
rk－浸透率降落系数 单元：无量纲
例：
P－REDU 1 6 5000
0.01.0
0.05 3.0
0.10 5.0
0.15 6.1
0.20 7.0
0.25 7.9
款式2：
P－REDU 2 brk crk
阐明：
brk－盘算浸透率低落因子RK的参数
单元：体积分数/分量百分数
crk－盘算Rk的参数
单元：
例：
P－REDU 2 1000 0.0186
12.4 P－ADSO
要害字P－ADSO用于输出聚合物的可及体积和吸附方面的参数。
此要害字有两种款式任选其一。
款式1：
P－ADSO 1 ephi4 npads csepad
阐明：
ephi4 －聚合物的可及体积分数
npads－吸附曲线的数据点数
csepad－侧吸附曲线的含盐量 单元：mg/L
假如npads＞0，则前面跟有npads个数据行，
款式：
xpads ypads
阐明：
xpads－聚合物浓度 单元：mg/L
ypads－聚合物吸附量 单元：mg/g
例：
P－ADSO 1 0.8 8 5000
0.0 0.0
2500.02
500 0.04
7500.065
10000.07
12500.07
15000.071
1750 0.071
款式2：
P－ADSO 2 ephi4 ad41 ad42 b4d
阐明：
ephi4 －聚合物的可及体积分数
单元：分数
ad41－聚合物的吸附参数
单元：无量纲
ad42－聚合物的吸附参数
单元：ml/meq
b4d－聚合物的吸附参数
单元：体积分数/分量百分数
例：
P－ADSO 0.85 0.7 0 100
13、T－DATA
要害字T－DATA引导用户输出关于示踪剂的参数。在T－DATA行前面用下列要害字输出具休数据：
T－NT
T－NAME
T－TK
T－TKS
T－RDC
T－RET
13.1 T－NT
用于输出示踪剂的
T—NT 用于输出示踪剂的数量.
款式:
T－NT nt
阐明:
nt ----- 示踪剂数量
nt = 1,2,3
13.2 T－NAME
T－NAME用于输出示踪剂的名字
款式:
T－NAME tname i (i=1,NT)
阐明：
tanme i－第i个示踪剂的名字不凌驾10个字符。
例：（NT＝3）
T－NAME TRACER1 TRACER2 TRACER3
13.3 T－TK
T－TK用于输出示踪剂在油水中的分派系数
款式:
T－TK tki (i=1,NT)
阐明:
tki－第i种示踪剂在初始条件下在水中的分派系数，
假如tki＝0，表现是水示踪剂
假如tki＝－1，表现是油示踪剂
例：(NT=1)
T－TK 0
13.4 T－TKS
T－TKS用于输出盘算示踪剂分派系数的参数，此参数作为含盐质变化的系数来盘算分派系数的变革。
款式：
T－TKS tksi (i=1,NT)
阐明：
tksi－第i个示踪剂的参数
单元：（meq/ml）－1
例：(NT=1)
T－TKS 0.003
13.5 T－RDC
T－RDC用于输出示踪剂的放射性衰减系数。
款式：
T－RDC rdci （i=1,NT）
阐明：
rdci－第i个示踪剂的放射性衰减系数
单元：1/day
例：(NT=1)
T－RDC 0.03
13.6 T－RET
T－RET用于输出示踪剂的滞后因子。
款式:
T－RET reti (i＝1，NT)
阐明：
reti－第i个示踪剂的滞后因子
单元：无量纲
例：(NT=1)
T－RET 0.061
14、SLUG
要害字SLUG引导用户输出静态井数据（事情制度）。当事情制度改动或注入浓度改动时就要用SLUG输出一次。因而SLUG即将在数据流中屡次呈现。
SLUG行后紧跟nw个数据行，每口井占一行。数据行关于差别的井别有差别的款式，共有4种款式：
款式1：（用于定产注入井，iflag=1）
iwd q cca ccs ccp iiw
GELCC cgel1 cgel2 cgel3 cgel4(*)
FOAMCC cfoam1 cfoam2 cfoam3 (*)
TRACERCC cti (I=1,NT) (*)
阐明：
iwd－井的序号
q－井的注入速率
单元：m3/day或ft3/day
cca－碱的注入浓度
单元：wt%
ccs－表活剂的注入浓度
单元：wt%
ccp－聚合物的注入浓度
单元：wt%
iiw－注入水号
iiw=1 或 2
cgell1,cgel2,cgel3,cgel4 是GEL各组分的浓度,如没有GEL注入,则没有GELCC这一行.
Cfoam1,cfoam2,cfoam3 是FOAM各组分的浓度, 如没有FOAM注入, 则没有FOAMCC这一行.
cti－第i种示踪剂的注入浓度
单元：wt%
如没有示踪剂,即NT=0, 则没有TRACERCC这一行.
例：
1 200 1.2 0.3 0.15
款式2：（用于定压注入井，iflag＝3）
iwd pwf cca ccs ccp iiw
GELCC cgel1 cgel2 cgel3 cgel4(*)
FOAMCC cfoam1 cfoam2 cfoam3 (*)
TRACERCC cti (I=1,NT) (*)
阐明:
iwd－井的序号
pwf－井底压力
单元：kpa 或psi
cca ccs ccp 等,同款式1
例:
2 12000 1.2 0.3 0.18
款式3：（用于定产采出井，iflag=4）
iwd q
阐明:
iwd－井的序号
q－井的产液速率
单元：m3/day或ft3/day
例：
5 106
款式4：（用于定压采出井，iflag=2）
iwd Pwf
阐明:
iwd－井的序号
pwf－井底压力
单元：kpa或psi
例：
8 10100
15、DT
要害字DT用于输出工夫步长控制，它在数据流可呈现屡次。
款式：
DT kdt dto dclim dtmin dtmax
阐明：
kdt－工夫步长控制方法标识
kdt=1 定步长
kdt=2 变步长
dto－当kdt=1时，表现定长长值
当kdt=2时，表现初始步长值
单元：day
delim－每工夫步浓度变革限定
单元：分数
dtmin－最小的工夫步长限定
（kdt＝2）
单元：day
dtmax－最大的工夫步长限定
（kdt＝2）
单元：day
例:
DT 1 0.5
DT 2 0.05 0.005 0.01 5
16、ENDFW
要害字ENDFW输出最大含水值,用以控制模仿完毕工夫.。
款式：
ENDFW endfw
阐明：
endfw－含水值.
例：
ENDFW 0.98
17、PERFOR
要害字PERFOR输出井的射孔信息，它可以屡次呈现。假如没有PERFOR行，标明井对一切层射孔。
一个PERFOR行输出一口井的射孔信息。这口井改动射孔时再用PERFOR行。
款式：
PERFOR wname Ki (i=L1,L2)
阐明：
wname－井的名字
Ki－第i段的射孔标记
Ki=0不射孔
Ki=1射孔
L1－井的肇始段
L2－井的完毕段
例：（NZ=4）
PERFOR A-71 1 1 0 1
18、TIME
要害字TIME输出用户必要输入模仿盘算后果的工夫，它该当呈现很多次。
款式：
TIME t
阐明：
t－输入模仿盘算后果的工夫
单元：day
例：
TIME 30
TIME 1500
19、STOP
STOP是数据流完毕的要害字
款式：
STOP
例:
STOP
六、数据流的例子
C --- TITLE
TTT
GUDAO A-S-P FLOOD TEST
C DATE YERS MONTH DAYS
DATE 1997 10 1
C
TMAX 15000 52.5
C
C --- METRIC
METRIC
C
C --- INJEC k-a k-s k-p k-t k-g k-acid
INJECT 1 1 1 0 0 0 0 1
C
C --- GRIDMD kcoord kgrid
GRIDMD 1 1
C
C *** nx ny nz nw
10 10 2 2 0 0
C *** dx(i) i=1 nx
10*50.
C *** dy(i) i=1 ny
10*50.
C *** wname(m) iw(m) type(m) rw(m) swell(m) dir(m)
C ifirst(m) ilast(m) prf(m) pwfmin(m) pwfmax(m)
C qtmin(m) qtmax(m) m=1 nw
C
AA1 2 2 1 .15 0. 3 1 2 0.0 125000. .0 500.
AA2 9 9 4 .15 0. 3 1 2 0.0 125000. .0 500.
C --- GRIDVL
GRIDVL
C
C
POR VV
200*0.3
C
DZNET VV
200*5.
C
DEPTH L1
100*1200.
C
PERMX VV
100*5000. 100*500.
C
C
PERMY MLT 1.00 PERMX
C
C
PERMZ MLT 0.01 PERMX
C
C
PW VV
100*12000. 100*12050.
C
C
SW VV
200*0.45
C
C
PARA
ROCK 0.000001 12000. 0.
OIL 0.96 70.0 0.000001 1.7
WATER 1.00 0.46 0.000001 1.3
KRHI .30 .25 .25 .95 1.45 2.7
KRLO .0 .0 0.7 1. 1.95 0.9
PC 0. 62.0689 2.0
C
C --- 1 2 3 4 5 6 7 8 9
C --- Cl- H+ Na+ Ca++ Mg++ CO3-- HCO3- H2CO3 OH-
C
FWCHEM 6000. 0.000002 100. 99. 20. 20. 700. 0. 0.
C
C --WCI1 CL- H+ Na+ Ca++ Mg++ CO3--HCO3-H2CO3 OH-(ppm)
IWCHE1 6000. 0.000002 100. 99. 20. 20. 700. 0. 0.
C --WCI2 CL- H+ Na+ Ca++ Mg++ CO3--HCO3-H2CO3 OH-(ppm)
IWCHE2 6000. 0.000002 100. 99. 20. 20. 700. 0. 0.
C
C
C
C
C --- SURFAC
S-DATA
IFTTAB 36
0. 0. 20.
0. 0.1 0.5
0. 0.2 0.08
0. 0.3 0.006
0. 0.4 0.0009
0. 0.5 0.0008
0.5 0. 5.
0.5 0.1 0.2
0.5 0.2 0.05
0.5 0.3 0.004
0.5 0.4 0.0007
0.5 0.5 0.0006
1. 0. 3.
1. 0.1 0.015
1. 0.2 0.02
1. 0.3 0.002
1. 0.4 0.0005
1. 0.5 0.0004
1.5 0. 1.
1.5 0.1 0.10
1.5 0.2 0.009
1.5 0.3 0.0009
1.5 0.4 0.0003
1.5 0.5 0.0002
2.0 0. 0.5
2.0 0.1 0.05
2.0 0.2 0.006
2.0 0.3 0.0006
2.0 0.4 0.0001
2.0 0.5 0.00009
3.0 0. 0.1
3.0 0.1 0.01
3.0 0.2 0.003
3.0 0.3 0.0003
3.0 0.4 0.0001
3.0 0.5 0.0001
C
C
C S-CMC 0. 1. 0.00008
C S-IFT 13. -16.5 .007 13. -16.1 .010
C S-TRAP 1865. 59074 364.2
S-DEN 0.98 0.
C S-ALPH 0.5 0.5 0. 0.9 0.7
C S-PC 0. 62.0689 2.0
C
S-ADSO 1 1.0 15 5000.
0. 0.
500. 0.15
1000. 0.21
1500. 0.25
2000. 0.31
3000. 0.4
4000. 0.52
5000. 0.64
6000. 0.75
7000. 0.86
8000. 0.98
9000. 1.02
10000. 1.02
11000. 1.02
12000. 1.02
C
S-EXCH 0.05 .25 .2 419.
S-CSEL 0.012 0.544 0.8
C
C
C
C
C --- POLYM
P-DATA
C
VISTAB 18
1000. 0. 0.5
1000. 200. 4.5
1000. 600. 18.
1000. 1000. 35.
1000. 1500. 44.
1000. 2000. 50.
4000. 0. 0.5
4000. 200. 4.5
4000. 600. 17.
4000. 1000. 34.
4000. 1500. 43.
4000. 2000. 49.
7000. 0. 0.5
7000. 200. 4.4
7000. 600. 16.
7000. 1000. 33.
7000. 1500. 42.
7000. 2000. 48.
C
C
C P-AP 60. 1600. 0. 2. .01 -0.175
P-GAMM 3.219 20. 1.1
C P-REDU 1000. 0.0186
P-REDU 1 5
0.0 1.0
1000. 1.2
1500. 1.5
2000. 2.0
3000. 2.5
C
P-ADSO 1 0.8 8 5000.
0.0 0.0
250. 0.02
500. 0.04
750. 0.065
1000. 0.07
1250. 0.07
1500. 0.071
1750. 0.071
C
C
C
C --- ALKAL
A-DATA
A-TYPE 2
C A-PH 7.22 11.25 10.06 0.
A-PH 7.22 11.25
A-ACID 0.30 0.41 1.7 0.00000001 500.
C
C
C
C --- DT kdt dt dclim dtmin dtmax
DT 2 0.05 0.005 0.005 5.0
C
C --- ENDFW endfw
ENDFW 0.98
C
SLUG
C *** id(m) qt(m) c(1 m) c(2 m) .... c(n m) for flag(m)=1
C *** id(m) pwf(m) for flag(m)=2
C *** id(m) pwf(m) c(1 m) c(2 m) .... c(n m) for flag(m)=3
C *** id(m) qt(m) for flag(m)=4
C *** m=1 nw
1 100. 0. 0. 0. 1
2 100.
C
C
C --- TIME time (days)
TIME 1.
TIME 2.
TIME 5.
TIME 30.
C
C
C --- SLUG
SLUG
C *** id(m) qt(m) c(1 m) c(2 m) .... c(n m) for flag(m)=1
C *** id(m) pwf(m) for flag(m)=2
C *** id(m) pwf(m) c(1 m) c(2 m) .... c(n m) for flag(m)=3
C *** id(m) qt(m) for flag(m)=4
C *** m=1 nw
1 100. 0. 0. .05 1
2 100.
C
C
C --- TIME time (days)
TIME 60.
TIME 90.
TIME 120.
TIME 150.
TIME 200.
TIME 250.
TIME 300.
C
C
C --- SLUG
SLUG
C *** id(m) qt(m) c(1 m) c(2 m) .... c(n m) for flag(m)=1
C *** id(m) pwf(m) for flag(m)=2
C *** id(m) pwf(m) c(1 m) c(2 m) .... c(n m) for flag(m)=3
C *** id(m) qt(m) for flag(m)=4
C *** m=1 nw
1 100. 1.2 0.3 0.15 1
2 100.
C
C
C
TIME 350.
TIME 400.
TIME 450.
TIME 500.
TIME 550.
TIME 600.
TIME 650.
TIME 700.
TIME 750.
TIME 800.
TIME 850.
TIME 900.
TIME 950.
TIME 1000.
TIME 1050.
TIME 1100.
TIME 1200.
C
C --- SLUG
SLUG
C *** id(m) qt(m) c(1 m) c(2 m) .... c(n m) for flag(m)=1
C *** id(m) pwf(m) for flag(m)=2
C *** id(m) pwf(m) c(1 m) c(2 m) .... c(n m) for flag(m)=3
C *** id(m) qt(m) for flag(m)=4
C *** m=1 nw
1 100. 0. 0. 0. 1
2 100.
C
C
TIME 1300.
TIME 1400.
TIME 1500.
TIME 1600.
TIME 1700.
TIME 1800.
TIME 1900.
TIME 2000.
C STOP
C
C DT 2 0.05 0.006 0.005 5.
C
TIME 2200.
TIME 2600.
TIME 3000.
TIME 3600.
ENDFW 0.98
TIME 4000.
TIME 4500.
TIME 5000.
TIME 5500.
TIME 6000.
TIME 7000.
TIME 8000.
TIME 9000.
TIME 10000.
TIME 11000.
TIME 12000.
TIME 13000.
TIME 14000.
TIME 15000.
C
C
STOP
END
七、后果输入文件
用户在用CDNS软件模仿盘算一个方案后,可以失掉5类后果数据文件..
这些数据文件有下列用处:
1.反省运转历程的准确性和错误;
2.油藏工程技能职员剖析目标;
3.绘制种种曲线;
4.绘制种种等值图和三维可视图;
5.绘制流线图.
1、数据反省信息文件
CDNS输入2个数据反省信息文件：
$1.check1
$1.check2
此中, $1 表现模仿盘算方案数据流文件的名字.
$1.check1是初始化数据反省信息文件.
$1.check2是模子盘算数据反省信息文件.
2、目标剖析信息文件
CDNS输入3个目标剖析信息文件：
$1.report0
$1.report1
$1.report2
此中, $1 表现模仿盘算方案数据流文件的名字.
$1.report1是油田目标数据文件.
$1.report2是井目标数据文件.
$1.report0 是从$1.report1中提出的14项目标.
3、曲线绘制数据文件
CDNS输入2个曲线绘制数据文件：
$1.field
$1.wells
此中, $1 表现模仿盘算方案数据流文件的名字.
$1. field是针对油田目标绘制曲线的数据文件.
$1. wells是针对井目标绘制曲线的数据文件.
4、等值图和三维可视图绘制数据文件
CDNS输入4个等值图和三维图绘制的数据文件：
$1.stat
$1.psvi
$1.conc
$1.chem
此中, $1 表现模仿盘算方案数据流文件的名字.
$1.stat是静态数据场,包罗深度,厚度,孔隙度,浸透率,初始水相压力和初始水饱和度.
$1. psvi是静态数据场,包罗相压力,相饱和度,相粘度和界面张力.
$1.conc 是静态数据场,包罗各组分的浓度.
$1.chem. 是静态数据场,包罗各化学离子的浓度.
5、流线图绘制数据文件
CDNS输入1个流线图绘制的数据文件：
$1.lx000
此中, $1 表现模仿盘算方案数据流文件的名字.
$1.lx000纪录了各工夫阶段各层流线图的轨迹,可间接用于绘制流线图.
