化学驱数值模拟软件
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记载了各时间阶段各层流线图的轨迹,可直接用于绘制流线图.
