admin 管理员组

文章数量: 887021


2023年12月18日发(作者:spring请求流程)

Eclipse 100 油藏数模软件使用手册

二OO四年十月

目 录

1 Eclipse 油藏模拟软件特点……………………………………… 1

1.1 Eclipse软件91年A版本的新进展概况………………………… 1

1.2 Eclipse100软件特点……………………………………………… 1

2 数据文件综述……………………………………………………… 12

2.1 RUNSPEC部分………………………………………………… 15

2.2 GRID部分………………………………………………………19

2.3 EDIT部分………………………………………………………24

2.4 PROPS部分…………………………………………………… 25

2.5 REGIONS部分………………………………………………… 31

2.6 SOLUTION部分…………………………………………………32

2.7 SUMMARY(汇总)部分………………………………………… 35

2.8 SCHEDULE部分…………………………………………………42

3 关键字描述(按字母顺序排列)…………………………………… 47

ACTNUM 活节点的识别………………………………………………47

ADD 在当前BOX中指定的数组加一个常数……………………48

ADDREG 给某一流动区域内指定的数组加一个常数………………49

ADDZCORN 给角点深度数组加一个常数………………………………49

APIGROUP 给API追踪中的油PVT表分组……………………………51

APIVD API追踪平衡的深度与原油比重(API)的关系……………51

AQANTRC 指定分析水层的示踪剂浓度………………………………51

AQUANCON 定义分析水层的相关数据…………………………………52

AQUCON 数值化水层与油藏的连接…………………………………53

AQUCT 说明Carter—Tracy水层的特征数据……………………54

AQUFET Tetkovich水层说明数据…………………………………55

AQUFETP 说明Fetkovich水层的特征数据…………………………56

AQUNUM 给一个网格块赋值一个数值化水层………………………57

AQUTAB Carter—Tracy水层的影响函数表………………………58

BDENSITY 盐水地面密度………………………………………………59

BOUNDARY 定义在打印网格表中显示的网格范围………………… 59

BOX 重新定义当前输入的BOX…………………………………60

CECON 生产井射开节点的经济极限……………………………… 61

COLLAPSE 识别在压缩VE选择中可压塌的单元………………………62

COLUMNS 设置输入数据文件的左右范围…………………………… 62

COMPDAT 井完井段说明数据………………………………………… 63

1

COMPFLSH 井射孔段的闪蒸转化比………………………………………65

COMPIMB 井射开网格的渗吸表号………………………………………67

COMPINJK 用户定义的注入井相对渗透率………………………………68

COMPLUMP 为自动修井而将射开网格归在一起…………………………69

COMPRP 重新标定井射开节点的饱和度数据……………………… 70

COMPVE 垂直平衡(V.E.)运行时,井射孔深度的重设定 …………72

COORD 坐标线……………………………………………………… 75

COORDSYS

COPY

COPYBOX

CRITPERM

DATE

DATES

DATUM

BEBUG

DENSITY

DEPTH

DIFFC

DIFFDP

DIFFMMF

DIFFMR

DIFFMTHT

DIFFMX

DIFFMY

DIFFMZ ZDIFFR

DIFFTHT

DIFFX

DIFFY YDIFFZ

DPGRID

DR RDRSDT

DRV RDRVDT

DTHETA

DTHETAV

坐标系统信息…………………………………………………76

从一个数组拷贝数据到另一数组 ………………………… 77

从一个BOX向另外一个拷贝一组网格数据 ……………… 77

对VE节点压缩的渗透率标准 ………………………………78

输出日期到汇总文件…………………………………………79

模拟者事先指定报告日期……………………………………79

基准面深度,用于深度校正压力的输出 ………………… 80

控制检测输出 ………………………………………………80

地面条件下流体密度……………………………………… 81

网块中心深度……………………………………………… 82

每一个PVT区域的分子扩散数据……………………… 82

在双重介质运行中,限制分子扩散……………………… 83

基质一裂缝的扩散乘子…………………………………… 83

方向的扩散乘子……………………………………………83

θ方向扩散系数乘子…………………………………………84

方向的扩散乘子……………………………………………84

方向的扩散乘子……………………………………………85

方向的扩散乘子 ………………………………………… 85

方向的扩散系数……………………………………………86

θ方向的扩散系数………………………………………… 86

方向扩散系数………………………………………………87

方向扩散系数………………………………………………87

方向扩散系数 …………………………………………… 88

对裂缝单元使用基质单元的网格数据…………………… 88

方向网格的大小……………………………………………88

溶解GOR的增加的最大速度…………………………………89

方向网格大小(矢量) …………………………………… 89

挥发油的OGR的增加的最大速度…………………………… 90

θ方向的网格大小…………………………………………… 90

网格的角度大小(向量) …………………………………… 91

2

R X Y R X Z

DX X方向的网格大小………………………………………… 91

DXV X方向网格大小(向量)………………………………………91

DY Y方向网格大小………………………………………………92

DYV Y方向网格大小(向量)………………………………………92

DZ Z方向网格大小………………………………………………92

DZMTRX 基质块的垂直尺寸……………………………………………93

DZMTRXV 基质岩体块的垂直尺寸(向量)…………………………… 93

DZNET

ECHO

EDITNNC

EHYSTR

END

ENDBOX

ENDNUM

ENKRVD

ENPTVD

EQLNUM

EQUALS

EQUIL

EXTRAPMS

FIPNUM

GCONINJE

GCONPRI

GCONPROD

GCONSALE

GCONSUMP

GCONTOL

GECON

GLIFTLIM

GRAVITY

GRIDFILE

GRUPRIG

GRUPTREE

GSEPCOND

IMBNUM

IMBNUMMF

IMPES

净厚度……………………………………………………… 93

接通重复输出开关……………………………………………94

改变非相邻连接………………………………………………94

滞后作用参数和模型选择………………… ……………… 95

标志SCHEDULE部分的结束 ………………………………95

将BOX恢复到包含全部网格…………………………………95

端点标定与深度区域号 ……………………………………95

相对渗透率端点与深度关系表 …………………………… 96

饱和度端点与深度关系表……………………………………97

平衡区号数 ……………………………………………… 98

在目前的BOX中设置数组为常数………………………… 99

平衡数据详述…………………………………………… 99

对表的外插请求预告信息 ……………………………… 101

流体储量区域号…………………………………………… 102

对井组井/油田注入率的控制/限制 ………………………102

为“优先”而设的井组或油田产量限制………………………104

井组或油田的产率控制或限制…………………………… 104

井组或油田的售气控制产率……………………………… 107

井组的气消耗率和引进率 ……………………………… 109

井组控制目标(产率)允许差额 ……………………………110

井组或油田的经济极限数据……………………………… 111

最大井组人工举升能力…………………………………… 112

地面条件下的流体密度 ……………………………………113

控制几何文件网格的容量 ……………………………… 113

给井组配置修井设备……………………………………… 113

建立多级井组控制的树状结构 ……………………………114

井组设置分离器…………………………………………… 115

渗吸饱和度函数据区域号………………………………… 115

基质—裂缝渗吸区域号…………………………………… 116

建立IMPES求解过程 …………………………………… 117

3

IMPLICIT 重建全隐式求解 …… …………………………………… 117

INCLUDE 包含数据文件名 ……………………………………………117

INIT 要求输出初始文件 ……………………………………118

INRAD 径向模型的内径 ………………………………………… 118

KRG 标定气相对渗透率的端点 ………………………………… 118

KRNUM 方向性相对渗透率表格数 ……………………………… 119

KRNUMMF 基岩—裂缝流动饱和度表号 ……………………………… 120

KRO

KRW

LOAD

MESSAGES

MINPV

MINPVV

MISCNUM

MONITOR

MULTIPLY

MULTR RMULTTHT THETAMULTX XMULTY

MULTZ ZNEWTON

NEWTRAN

NEXTSTEP

NNC

NOECHO

NOGGF

NODPPM

NOWARN

NTG

OILAPI

OLDTRAN

OLDTRANR

OPTIONS

OUTRAD

OVERBURD

PERMR R

标定油相对渗透率端点 …………………………………… 120

标定水相对渗透率端点 …………………………………… 121

调入一个SAVE文件以便执行一个快速重起动…………… 122

重设置打印和停止限定的信息 …………………………… 123

设置活动网格的最小孔隙体积………………………………124

建立一个有效网格的最小孔隙空间 ……………………… 124

混合区数目 ……………………………………………… 125

请求实时显示输出 ……………………………………… 125

当前定义区中的数组……………………………………… 126

方向传导率乘子……………………………………… 126

方向传导率乘子………… …………………………127

方向传导率乘子………………………………………… 127

方向传导率乘子………………………………………… 127

方向传导率乘子 …………………………………………128

输出迭代计数到汇总文件 …………………………………128

标定使用块拐角传导率 ……………………………………128

建立下一时间步最大值 ……………………………………129

非相邻连接的直接输入 ……………………………………129

关闭输出的响应 ………………………………………… 130

压缩网格几何模型文件 ……………………………………130

非双孔的渗透率乘子 ………………………………………130

压制ECLIPSE警报信息…………………………………… 130

厚度净毛比………………………………………………… 130

初始原油API值,以便API示踪选择………………………131

标定块中心传导率……………………………………………131

标定任意一块中心传导率………………………………… 131

开启特别程序选择………………………………………… 132

径向模型外半径 …………………………………………… 134

岩石负载压力表……… ……………………………………135

方向绝对渗透率………………………………………… 135

4

Y

PERMTHT θ方向绝对渗透率………………………………………… 136

PERMX X方向绝对渗透率…………………………………………136

PERMY Y方向绝对渗透率 …… ……………………………… 136

PERMZ Z方向绝对渗透率…………………………………………137

PINCH 建立尖灭层上下的连接 ………………………………137

PINCHOUT 建立尖灭层上下的连接…………………………………138

PMAX

PMISC

PORO

PORV

PRESSURE

PRIORITY

PRVD

PSEUDOS

PVCO

PVDG

PVDO

PVTG

PVTNUM

PVTO

PVTW

PVTWSALT

QDRILL

RESTART

RESVNUM

ROCK

ROCKNUM

ROCKTAB

ROCKTABH

RPTGRID

RPTONLY

RPTPROPS

RPTREGS

RPTRST

RPTRUNSP

RPTSCHED

RPTSMRY

模拟中的最大压力………………………………………138

与压力有关的可混性表…………………………………138

网格孔隙度………………………………………………139

网格孔隙体积……………………………………………140

初始压力……………………………………………… 140

为井的优先级选项设置系数……………………………140

原始压力与深度关系表…………………………………142

为PSEUDO包要求输出的数据………………………… 142

含气原油PVT性质………………………………………142

干气的PVT性质(无挥发油) ………………………… 144

死油的PVT性质(无挥发气) ………………………… 145

湿气的PVT性质(有挥发油) ………………………… 145

PVT区数目………………………………………………146

活性油的PVT^性质(有溶解气) ………………………147

水PVT性质………………………………………………148

含盐的水PVT函数………………………………………149

在钻井队列中安置井……………………………………150

设置重启动………………………………………………151

对一给定油藏输入角点坐标数据………………………153

岩石压缩系数……………………………………………153

岩石压实表格区数………………………………………154

岩石压实数据表…………………………………………154

滞后岩石压实数据表……………………………………155

从GRID部分输出控制………………………………… 156

摘要输出的常规限制……………………………………158

控制PROPS部分的输出…………………………………158

控制REGIONS部分的输出………………………………159

输到RESTART文件的控制………………………………159

控制RUNSPEC部分的数据输出…………………………160

控制SCHEDULE部分的输出…………………………… 160

控制SUMARY部分的输出……………………………… 163

5

RPTSOL

RS

RSCONST

RSCONSTT

RSVD

RUNSUM

RV

RVCONST

RVCONSTT

RVVD

SALT

SALTVD

SAVE

SCALELIM

SDENSITY

SEPVALS

SGAS

SGCR

SGFN

SGL

SGOF

SGU

SIGMA

SIGMAV

SLGOF

SOF2

SOF3

SOGCR

SOMGAS

SOMWAT

SORWMIS

SOWCR

SPECGRID

STOG

STONE1

STONE2

STOW

控制SOLUTION部分的输出………………………………163

初始溶解气油比……………………………………………165

为死油设置的一个常数Rs值…………………………… 165

为每一个死油PVT表设置的一个常数Rs值…………… 166

用于平衡选择的RWJ深度关系表…………………………166

所需的SUMMARY数据的制表输出…………………………167

初始挥发油气比……………………………………………167

为干气设置的一个常数Rv值…………………………… 167

为每个干气PVT表设置一个常数Rv值……………………168

用于平衡选择的Rv与深度关系表…………………………168

初始盐浓度 ………………………………………………169

用于平衡的盐浓度与深度关系……………………………169

用于快速重启文件而需输出的SAVE文件……………… 170

设置饱和度表的标度限制…………………………………170

在地面条件的混相气密度…………………………………170

分离测试的Bo和Rs值……………………………………171

初始气饱和度………………………………………………173

临界气饱和度的标度………………………………………173

气体饱和度函数……………………………………………174

原生气饱和度的标度…………………………………… 175

气/油饱和度函数与气饱和度…………………………… 176

最大气饱和度的饱和度表的标度…………………………177

双重孔隙基岩—裂缝的连结………………………………178

双重孔隙度基岩—裂缝的连结(向量) ………………… 178

气/油饱和度函数与液体饱和度………………………… 179

油饱和度函数(2相) ………………………………………180

油饱和度参数(3相) ………………………………………181

临界的气中含油饱和度的标度……………………………182

STONE1模型中含油饱和度最小值……………………… 183

STONE1模型中最小油饱和度值……………………………184

混相残余油饱和度数表……………………………………185

标度临界水中含油的饱和度值……………………………186

网格特性的详细说明………………………………………187

油气表面张力与压力………………………………………187

三相油相对渗透率模型……………………………………188

三相油相对渗透率模型……………………………………188

油水表面张力与对应压力…………………………………188

6

SWAT

SWATINIT

SWCR

SWFN

SWL

SWLPC

SWOF

SWU

TBLK

THPRES

TLMIXPAR

TNUM

TOPS

TRACER

TRACTVD

TRANR

TRANTHT

TRANX

TRANY

TRANZ

TSTEP

TUNING

TVDP

TZONE

VAPPARS

VEDEBUG

VEFRAC

VEFRACP

VEFRACPV

VFPINJ

VFPPROD

WBOREVOL

WCONHIST

WCONINJ

WCONINJE

WCONPROD

WCUTBACK

初始水饱和度 ………………………………………189

标定毛管压力的初始水饱和度 …………………190

临界水饱和度的标度………………………………190

水饱和度函数………………………………………191

原生水饱和度的标定………………………………192

仅对毛管压力曲线标定原生水饱和度……………193

水/油饱和度函数和对应的水饱和度…………… 193

饱和度数表中最大的含水饱和度的标定…………195

示踪剂的初始浓度…………………………………196

门限压力……………………………………………196

Todd-Longstaff混合参数…………………………197

示踪剂浓度区 ……………………………………198

每个网格的顶面深度………………………………198

被动的示踪剂名……………………………………199

为示踪剂要求“流率极限传输” ………………199

R方向的传导率……………………………………199

θ方向的传导率……………………………………200

X方向的传导率……………………………………200

Y方向的传导率……………………………………201

Z方向的传导率……………………………………201

把模拟器推向新的报告时间………………………202

设置模拟器控制参数………………………………202

初始示踪浓度与深度表……………………………204

过度带控制选择……………………………………205

油挥发控制 …………………………………… 205

对垂向平衡和压缩垂向平衡选择控制调整………205

垂向平衡曲线系数的应用…………………………206

垂向平衡拟毛管压力系数的使用…………………207

垂向平衡拟毛管压力系数的使用……………… 207

对注水井输入V.F.P表………………………… 208

对生产井输入V.F.P表………………………… 209

对井筒贮存设置体积………………………… 212

历史拟合井观测产量…………………………… 213

设有组控制的注入井的控制数据……………… 215

对注入井控制数据…………………………… 217

对生产井控制数据……………………………… 218

井减少限制…………………………………… 220

7

WCYCLE

WDRILRES

WDRILTIM

WECON

WEFAC

WELDEBUG

WELDRAW

WELOPEN

WELPI

WELPRI

WELSOMIN

WELSPECS

WELTARG

WGASPROD

WGRUPCON

WHISTCTL

WLIFT

WLIMTOL

WORKLIM

WPIMULT

WPLUG

WSALT

WTEST

WTRACER

ZCORN

井自动循环开与关……………………………… 222

防止在同一网格中同时开两口井…………………222

新井自动开钻的控制条件……………………… 223

生产井的经济极限数据………………………… 224

设置井的效率系数(为停工期)……………… 226

个别井的跟踪输出控制………………………… 226

设置生产井的最大允许压差…………………… 227

关闭或重开井或井的射开层…………………… 228

设置井的生产/注入指数值………………………229

设置井的优先数………………………………… 229

自动开井的最小含油饱和度…………………… 230

井的综合说明数据……………………………… 230

重新设置井的操作目标或限制………………… 232

为控制销气而设置的特别产气井…………… 233

为井组控制而给井设置指导产率……………… 234

给历史拟合井设置覆盖控制…………………… 235

自动换管串和升举的开关数据……………… 235

经济和其它限制的容差分数…………………… 236

每次自动修井所花的时间……………………… 237

用给定值乘以井射开层地地层系数…………… 237

设置井的回堵长度……………………………… 238

设置注入井的盐浓度…………………………… 238

命令对已关着的井进行周期性测试………… 239

给注水井设置示踪剂浓度……………………… 240

网格块角点的深度……………………………… 241

8

1 Eclipse 油藏模拟软件特点

1.1 Eclipse软件91年A版本的新进展概况

详细说明见附录B1

1.新功能

(1)提供了可供选择的通用的油PVT数据和饱和度数据的输入关键词;

(2)对每一个PVT区设计了恒量Rs或Rv值;

(3)分子扩散选择能模拟气的扩散和油的组份;

(4)盐水选择能模拟不同矿化度盐水的流动。

(5)独立的饱和度表能用于基质—裂缝介质中的流动计算;

(6) 设置了先用驱替曲线平衡计算,后用吸吮曲线模拟的选择;

(7)端点刻度选择可以用来刻度相对渗透率端点;

(8)在模拟时间一节中,饱和度和PVT区可以改变;

(9)可以设置有效网格节点最小孔隙体积的下限值;

(10) 若网格扩展,对O基面深度可以用负值;

(11) 用块中心模型,设置了垂向平衡计算的块内流动新选择;

(12) 当分离器条件改变时,井/井组流动可以变换;

(13) 井压差限制可用于气井;

(14) 对注水井可以设置恒量相对渗透率;

(15) 用一个给定系数,刻度井的连通系数;

(16) 一个关键字能控制单井历史拟合的全过程;

(17) 在汇总数据文件中输出油田和区块的油采收率;

(18) 在汇总数据文件中输出每个时间步长非线性重复数;

(19) 在汇总数据文件中输出监测天然气流动(GRAF)运行时间的附

加资料;

(20) 在重新起动文件中还能输出水、气相的压力值;

(21) 网格文件中还包括了无效节点的资料。

2.对原有程序的改进

(1)在没有变化WCONHIST关键词中,垂向平衡压力(VEP)表号和ALQ值的缺省;

(2)一定方向的相对渗透率表能用于断层产生的不相邻连接;

(3)能减少模拟软盘的读/写时间;

(4)用缺省值代替示踪剂追踪的寄存选择;

(5)当选用等Rs/Rv值时,在模拟运行停止前就打印出网格内压力、饱和度交会表;

(6)在向量计算机上,提高了油藏特性表的查阅速度。

9

手册和关键词描述部分都是新的,或者说,从上一版手册以来,对这二部分进行了修改,这个修改可在题目一行开始处加一个星号来识别,即:*

1.2 *Eclipse软件特点

概 论

Eclipse100是一个全隐式的,三维、三相、还包括天然气、凝析油选择功能综合的黑油模拟软件。该软件用Fortran77语言编写,无论是虚拟储存,还是足够的实际储存,都能在用ANSI码标准Fortran77编译程序的任何计算机上运算,能运算该软件的计算机有:

microvax,vax,sun,Apollo,DG,convex,IBM,Cray等机。包括在Eclipse100软件包中的辅助程序有:

GRAF:一个独立图形后处理程序;

PSEUDO:生成三维拟函数程序;

Fill:角点模型前处理程序;

VEP:井筒水动力计算前处理程序;

EDIT:特别为准备Eclipse模拟数据而设计的屏幕编辑程序。

自由格式输入

Eclipse的输入资料是用关键字系统自由格式输入的。任何标准的编辑程序常来编辑输入文件。EDIT是一个可供选择的专门用于屏幕编辑的Eclipse编辑程序。当数据输入后,EDIT程序能检查输入数据,Edit的辅助程序包括大部分Eclipse参考手册。

相选择

Eclipse100油藏模拟软件常用于1、2或3相系统的模拟。当油藏为二组份(油水,油气,气水)系统时,选用双相模拟,它既能节省计算机的储存,又能节省计算时间。除了模拟气溶解于油中以外(可变的泡点压力或可变的气/油比),Eclipse还常用于油在天然气中挥发的模拟(可变的露点压力或可变油/气比)。

图形选择

Eclipse软件提供的角点模型和常规块中心模型是很有用的。在1、2、3维模拟中,经向和笛卡尔块中心点模型选择是很有用的。三维经向模型能模拟0—360°界面上的园形流动(关键词COORDSYS)。

Eclipse角点模型是唯一的能对正确代表油藏最复杂的几何地质图形进行模拟。前处理程序Fill和GRID常用来准备Eclipse油藏模拟的角点数据。独立图形的处理GRAF和GRID程序能用多种方法显示网格。例如,在进行大型的三维模拟时,用户可要求同时显示在XZ方向的多条横剖面。网格的平面透明覆盖图是十分有用的,它能正确地对地质平面图进行检查。

10

对多断层的油藏,角点模型特别有用。网格能沿着断层线平面上变形,垂向上移动,对复杂性的剪状断层进行模拟。使用前处理程序(Fill),甚至对倾斜断层亦能容易地和精确地描述出来。Eclipse软件能自动的计算不相邻网格的传导率。Eclipse的计算程序能有效的对穿过平移断层的流动进行模拟。

角点模型别的功能还包括了管流模拟和局部网格加密模拟。

全隐技术

Eclipse软件利用全隐法来确保超长时间步长模拟的稳定性。务必确保用减少所有残差到最小的允许误差对非线性全隐方程进行精确解。物质平衡误差(残差和)是非常小的。牛顿法常用来解非线性方程。用所有变量彻底扩展雅各宾矩阵来保证二次方程的(快速)收敛。各种各样的专门方法常用来加速非线性模型的收敛。同时,用正交极小法来加速插入因式分解,来解每次牛顿迭代产生的线性方程。

惯常,全隐技术用于小范围的锥形研究,在单个时间步长内,研究区内许多节点孔隙体积内的流体能穿过靠近井筒的小网格。隐压显饱和半隐式技术绝不能用来解这样的问题,除非它们的时间步长被减到不切实际的极小值时。当全隐技术用于消元式联立线性方程时,只能用联立方法求解,不能用简单的序列法来解。这类小作用直接方法,如D4高斯消除法求解。大型模拟需用迭代法,如强隐函数程序方程(SIP)和线松驰法(ISOR),肯定是不能收敛的。这样,大部分的程序不能用全隐法来解大型的作业。在Eclipse软件中,同插入因式分解计算能取消这些限制,能有效和可靠的进行大型作业的模拟。

隐压显饱(IMPES)

虽然,全隐技术是Eclipse软件使用的标准解法,但有的时候用隐压显饱(IMPES)技术亦十分有利的,在Eclipse中这种选择是很有用的。隐压显饱技术可能是不稳定的,它只用于正规网格系统(不是小网格)和小时间步长的模拟计算,如历史拟合。进行Eclipse模拟时,在同一个运行期间,属不同次的模拟,同时用IMPES和全隐技术都是可能的。

插入因式分解计算程序

用正交极小法加入插入因式分解法解每次牛顿迭代产生的线性方程。这是最快的迭代法,并且还用来解大型多个稀疏性方程。在每次迭代时,插入因式分解技术准确地保存了资料。因此,物质平衡法产生的误差仅与非线性方程解的结果有关。因不相邻连接面形成的较远带状基质单元被包含在因式分解程序中,所以它运行时十分有效。这种计算与用直接法(N**

3)计算相比,计算费用的增加略快于油藏中有效节点数的增加(N* *5/4)。所以插入因式分解计算程序能特别好的适合于大型作业模拟。同时,也能解二相和三相的作业。附录中线性方程解部分详细说明了这种方法。

11

在每次迭代时,用插入因式分解计算程序计算一个新的搜索方向(近似解)。能使残差(误差)平方和最小的正交极小法要求每一个新的搜索方向垂直于以前的搜索方向。这样,前一个搜索方向必须存储起来,因而,正交极小法将用完有用的计算机记忆。幸运的是,计算机通常足以存储最后几个搜索方向。若干个搜索方向称为NSTACK,并且用户可在运行定义卡数据那部分中设置它。设置NSTACK=10,为缺省值。对于有限记忆计算机进行大型作业模拟,这种记忆管理的选择是非常有效的。设置NSTACK=0,进行记忆管理。于是,不受长度限制的堆栈在磁盘上储存起来,一当需要的时候,解算器仍能回放(回读)出来。由于内存数据存盘转换速度比较慢而增加了计算机CPU的钟时,故内存管理选择是用增加执行钟时来节省内存的。令NSTACK为负值,则可将完全机内叠加记忆和完全磁盘储存之间进行折衷。例如,若NSTACK为-9,于是,第一个8搜索方向保持记忆,而其余的搜索方向则储存在磁盘中。因为在每次运行开始,Eclipse软件能计算所需要的记忆量,并且用户能用试验的方法对NSTACK的特殊,找到最佳值。

由于使用了有效节点址,解算器完全是矢量性的。例如,所有的正交极小法的堆栈叠加计算全是矢量计算。矩阵相乘亦是使用矢量运行技术进行矢量计算的。然而,插入因式分解计算核心的三角转换是用回归法,不能用矢量计算。

不相邻连接

在正常的网格系统里,在每个方向上,每个节点仅有一个相邻的节点(维上一维的每个节点有二个相邻节点)。这样在二维网格里,每个节点有四个相邻节点,在三维网格里,每个节点可以有六个相邻的节点。流体流动仅发生在相邻节点之间。

Eclipse软件中选择性连接成对不相邻节点是有可能的,处于这种连接时,该软件能模拟流体从一个节点直接流动到另一对节点。不相邻节点连接的典型应用有:1.能进行3D经向园形模拟;

2.允许穿过平移断层在不同层之间进行流动模拟;

3.允许在井筒附近进行局部网格加密,它还包括了在一个笛卡尔座标系统内对锥形网格模拟的功能。

若前处理程序(FILL)常用来形成网格,那么Eclipse就能自动的计算上面第1、2二个典型例子中不相邻网格间的传导率。用一般的不相邻连接设置,能进行局部的网格加密。如何进行局部网格加密详见后面的手册。

不相邻网格连接能产生雅各宾矩阵中边缘带状单元,这个方法包括在插入因式分解程序中,并且该方法常用来进行线性方程解。

运行次数的标定

Eclipse所有的内部的数组都标有运行的次数,目的是减少计算机内

12

存的使用。例如在二相运算中,雅各宾矩阵中每个单元的运算次数为22,*而在三相运算中则为33。一个二维的运算作业就要使1个雅各宾矩阵产生5个带状单元,而三维运算作业中就会使一个雅各宾矩阵产生7个带状单元。显然,要包含最不利的可能情况这种内部数组运算次数的标定是低效率的。在运行定义卡数据准备部分中提供了ECLIPSE计算内部数组的大小所需的大部分数据。像自动形成不相邻连接,死节点,垂向平衡选择等另外所需的资料则被减少了,它们在网格部分中提供。不存储死节点不必要数据,节省了总的计算空间。在模拟开始前,每次运算所需要的总内存可打印出来。

*垂向平衡

在三维模型中,为了模拟分层流动,Eclipse软件还包含了唯一的垂向平衡选择。垂向平衡程序既可用于块中心模型,亦能用于角点模型。用角点模型,垂向平衡程序考虑每个变形网格的形态和方向。一个混合参数能代表用户划分的流体的分隔程度。一种情况,假设总流体是分层的,另一种极端情况假设流过每一个网格流体完全是分散的。例如,垂向平衡程序考虑了,当油侵入水区和被驱替后,留下的临界的残余油饱和度引起的滞后效应。经压缩的垂向平衡选择程序能将3D模型当成2D平面模型来运算。

双孔隙度

Eclipse双孔隙度/渗透率选择程序是适用于多裂缝油藏的模拟。用一个基质网格和一个裂缝网格代表每一个双重介质网格。在双孔隙度选择程序中,油藏内的流动只能通过裂网格才能发生。同样在每个基质网格和与它相对应的裂缝网格之间亦产生流动。在双孔单渗的情况下,在相邻的基质网格之间决不会发生流体的流动,在双孔隙度/渗透率选择程序中,在相邻基质网格之间才可以产生流动。

在基质网格和裂缝网格之间流体的重力吸吮/驱替和分子扩散等现象都能解释。

原油的高压物性(PVT)数据和岩石数据。

Eclipse软件可靠地保证用户所确定的压力和饱和度函数的数据。它决不会出现工业上常见的用内插法修改相同的饱和度间隔值园滑数据。

若干不同数据表可用于油藏的不同部分。这样,例如,一些单独的饱和度表可以输入每种岩石类型的模型中。

相对渗透率的方向性

相对渗透率曲线可以分别给定为水平方向流动和垂直方向流动的二个相对渗透率曲线。对每个饱和度函数区,用户可以有给定3条(X、Y、Z)或6条(+X,-X,+Y,-Y,+Z,-Z)相对渗透率曲线的选择,而不是常用一条曲线。这种新增的自由给定相对渗透率曲线的条数,油藏工程师能用

13

比较简单的方法提高油藏物理性质的模拟研究。例如,大部分油藏节点在垂向上分布少,在平面上分布多,如果水平相对渗透率小于巴克里-莱弗里特饱和度,那么水平流动模拟将会能更真实的模拟油藏实际情况。在Eclipse拟函数选择程序中,同样也需要相对渗透率的方向性,它可以用于三维粗网格的模拟中。

饱和度数据表的刻度

相对渗透率和毛细管压力曲线的束缚、临界和最大饱和度端点值可分别输入油藏内的每一个网格中或者选用这些端点值随深度变化表分别输入网格内的每一个模拟区内。饱和度刻定选择程序能让用户给定随标准化饱和度而变化的相对渗透率数据和毛细管压力数据。此外,该程序还能用初始临界饱和度或束缚流体饱和度随深度变化进行油藏模拟。用户还能给定饱和度表各向是同性的,但在X、Y、Z座标方向上有区别,或者在6个座标(+X,-X,+Y,-Y,+Z,-Z)方向上有区别。

饱和度刻度程序常可用于一个油藏内,原始地层可动流体是校正模拟,该校正是借助于流体接触界面处的一个特殊校正。

滞后效应

相对渗透率和毛细管压力的二个滞后效应都可模拟。对非润湿相,相对渗透率的滞后效应有二个模型可供选择:即Carlson模型和Killough模型。在二相模型的运算中,也可以选用Killough润湿相滞后模型。

岩石压实

当流体压力减少时,岩石压实选择程序能模拟孔隙通道的减小。这个过程可以是可逆的,不可逆的,或者需要时还可滞后。这个程序在双重介质的运算中特别有用:

示踪剂追踪

在模拟运算的期间内,示踪剂跟踪程序能确定流体组分的运动“标志”。例如,它可常用于辨别不同注水井的水运动情况和原始地层水的移动情况,或者预测地层水矿化度和其他种类化学物质浓度的变化。

示踪剂还可定义为存在于一个烃类相中,并且它可用来模拟在游离和溶解状态之间流体成分上的交换,无论什么时候,这种交换总是发生在伴生的烃类相中。所以,在预测原始游离气或溶解气的运动中是十分有用的。原始含硫量是可以用移动油的跟踪与深度的函数来定义。

在示踪剂界面上,控制数值扩散效应选择程序是非常有用的。

原油比重(API)跟踪

不像上面所说的被动的示踪剂跟踪程序那样,原油比重跟踪是主动的。它充分考虑了完全不同(PVT)性质原油的混合。与常规的假设相对照,原油比重跟踪程序完美的假设了在网格内流动原油的地层体积系数,原油粘度,原油密度和泡点压力。原始的原油比重在每一个平衡区内把它作为与

14

深度呈函数关系,而对每个网格来说,它是独立的。

海水盐度跟踪

在水相中,盐水跟踪程序是一个主动的示踪剂跟踪程序。表列了水的性质与盐的浓度有关。因此,若原始地层水与注入水有不同的矿化度,那么,就可说明不同水的密度和粘度来自何处了。

混相驱动能力

Eclipse软件具有二组份混相驱替模型。该程序用于混相注气过程的模拟。模型假设油藏流体由三个组分组成:油藏油(储藏油+溶解气),注入气(溶剂)和水。假设油藏油和溶剂气的组份按比例溶混。用Todd-Longstaff技术论述混相组份的物理弥散现象。Todd-Longstaff技术是借助于一个地区给定一个基本的混合参数来论述的。溶剂气的比重与溶解气的比重是不同的,该模型还允许油/溶剂在高含水饱和度研究区内的混驱结果在屏幕显示出来。一个专门数字扩散控制程序对混相驱替选择是很有用的,尤其在大的混合参数研究中,该程序特别有用。

单井控制程序

Eclipse软件有一套综合的单井控制选择程序。生产井能在给定的油量、水量、液量、气量,油藏流体废弃产量、井底压力和井口压力等数值下进行运算。油藏工程师应对上述数值中一项定出一个目标值,而对其余数值的各项定出限制值。生产井只要不超过限值,在规定目标值下工作。若一个限制值将要超过它的限制界限,那么生产井将自动的改变它的控制模式,使井保持在允许的限制条件下工作。效率系数还可用于考虑正常的停产时间。该程序将能计算处于正常流动条件下的流量和流动压力,但累计流量按照有效系数进行扣除。注水井亦有一套类似的控制程序。它能提供地面条件下的注入量,油藏条件下的注入量,井底压力,和井口压力的目标值和限制值。

有助于历史拟合的全过程,生产井可处于另处一种控制类型。输入这些井实际生产的油、水和气量,使井处于相同的产液量或油藏流体体积量的条件下生产。这样,甚至当含水率和气/油比不完全拟合时,压力下降的速度大致是正确的。实际生产的和计算的这二个产量都可编入汇总文件,以便进行图形对比。

生产井同样还可受到另一种的经济约束。如果一口生产井的产油量或产气量低于经济界限,那么该井将自动关井或重新完井。如果生产井的含水率、气/油比或水/气比超过规定的上限,那么生产井就要修井或关井。修井速度取决于占机提供修井的能力。生产井的含水率、气/油比和水/气比的生产界限同样可以用于同生产井连通的其它井,若这些界限超过规定,也应关井。此外,每当生产井的超过界限时,我们可以逐渐地减少生产井的产量来限制生产井含水量,气/油比或水/气比的上升。生产井还受到最

15

小井口压力的限制,一旦生产井的产油量或产液量低于规定的界限,可以用自动开关或人工举升方法来提高它们的生产能力。对于人工举升的生产井还受到井组或油田产能的限制。一些生产井由于某种理由关井要设置一个使模拟自动停止的标记,而其它井要有自动开井的标记。若油田上所有的生产井,或一个或更多的井组关井,那么这个运算同样可以停止。关闭的井可以进行定期测试,了解它们是否具有再次生产的经济效益。

井组、油田的生产控制

整个说来,模拟程序包含了井组、油田动态资料采集的控制程序。用优先或目标产量控制法,使一个井组或几个井组总的产油量、产水量、产气量或产液量都能达到规定的目标值。

用优先产量控制,将井按优先的次序排列,按优先的递减次序生产直到一个井组的产量超过界限为止。超过井组界限产量的井应取消它的生产量以适应井组的界限。井组中其于低优先级的井应关井,直到需要时间开井。井的优先率既可以由用户自己设置,亦可以用用户定义系数通式,在规定地井段中进行计算。

用指导产量控制,井组的目标产量是按各个生产井的可能生产量或规定的指导产量按比例的分配到各个井去,但必须确保没有一口井会超过它的单井产量或压力界限。当井组的产能达不到井组的目标产量时,井组产量就要递减。然而,一当需要补充产量时,该程序可以按顺序自动打开备用新井来延迟产量的递减,同样该程序也受到钻速和每个井组最多允许开井数的限制。

井组的上限值同样可以用于井组的产油量、产水量、产气量和产液量。一当这些界限值的其中之一超过时,那么,为了保持井组的生产量在界限以上时,从最差井的修井工作量到井组控制的应用范围必须采取选择。此外,上面所说的一套单井经济界限和相应的描述能应用于井组的所有动态分析中。

上述同样一套井组控制和界限同样可以应用到全油田的生产动态分析中。油田处于指导产量的控制下,油田的目标产量或是按规定的井组指导产量按比例的分配到各个井组,或是按各井可能生产量或指导产量按比例直接分配到各生产井。油田的生产量同时可用一个井组或它的更多子井组的生产来控制。控制流体相态(油、水、气或液)的一致性,在每一种情况下不需要是相同的。例如,油田可按一定的需要量生产油,而其中的一个井组正在按其最大的产气能力生产。这样,油田将按照它的目标产量生产,每个井组将严格地受它目标产量的约束,而每一口单井则严格地受它的单井压力和单井产量限制。

多级开采系统

用表明某些井组是高一级井组的子井组的办法,可将井组按任一所需

16

的层次排列成一个多级树式结构。所有井组的控制和限制条件都可应用任何一级的井组(包括油田亦可看作成树状结构顶部的0级井组)。

注水井组的控制

井组和油田可以给定目标注入量和限制注入量来控制它们的注水量。油藏工程师可以设置下列各项参数目标值和限制值:

(1)单相流体的地面注入量;

(2)多相流体注入油藏的总体积量;

(3)单相流体的回注比例;

(4)多相注入流体总的注采比。

所以,例如一个井组可以按它产气量的一定比例进行回注,并且从注入水量得到一个油藏总的注入目标或总注采比。

一个井组的目标注入量或是按它们的注入能力或是按它们额定的注入量按比例分配到各注水井中,一个井组的目标注入量还受到各注水的分层流动和压力的限制。一个高一级井组的目标注入量(或油田注入量)可按它们井组的额定注入量有比例的分配到子井组中去,井组的额定注入量可以直接设置,或者是它们生产的亏空量,或者是需要的排替亏空量。若一个注水井组(或油田)注入能力不足以满足于它的目标注入量时,那么,为了打开一口合适的注入井,应当仔细研究预备的注入井钻井顺序。

销售气生产控制

销售气生产控制程序提供了一个从油田伴生产输出的控制方法,这个方法同时可对油的生产量进行控制。一个井组或油田的商品气的生产量可看作为总的产气量减去它总的注气量,再减去井组和子井组气的消耗量,加上输入到井组或子井组的气量(若有的话)。对每个井组可单独设置气的消耗量和输入量。

销售气的生产量还受到不需要销售过剩的回注气量控制。销售气控制程序能独立的应用于油田或井组另外的生产控制,同时,该程序还能提供有足够的注入能力来注入过剩气、井组或油田自动的处于伴生气回注控制之下,它的目标回注比例是按照每个时间步长内注入过剩气量的动态来确定。

井的层间交互和混流

当一口井完井于几个网格区内,那么每一个网格的流量与产生的下列三个参数成正比:

(1)网格和井筒之间的传导率;

(2)网格内射孔段相的流度;

(3)网格和井筒之间的压力降。

当在不同射开网格间分配井的目标产量时,Eclipse软件充分考虑了上述的三个参数。这对于在有二个或二个以上的,垂向连通性比较差的油

17

层,每个油层的压力降又是各不相同的。完井特别重要。

在某种情况下,一个层的压力降可能与另外几个油层压力降不一样。那么,该程序将允许通过井筒在油层之间产生交互窜。在井口能观察到物质平衡的约束,以至于流出井筒的混合相相当于进入井筒的平均混合相。Eclipse软件能进行任何控制模式下生产井的任何类型的交互窜模拟,还包括井内上部地层被封堵,而其余部分的油层生产情况的模拟。如果油藏工程师不愿意有交互窜产生,则可以关闭交互窜流模拟的装置。

大斜度井和水平井

每口井完井的位置上是没有限制的。一口大斜度井可以射开不属于同一垂直列的几个网格。实际上,一口井的同一个层,可以射开几个网格,这种情况,这个程序能模拟射开几个相邻网格的水平井。这个功能同样在3D径向锥进模型研究中十分有用,3D径向锥进模型研究地区是指:与井连接最近的地层,与井连通的楔状砂体形态以及与井连通的几个网格区内的已知油层。Eclipse200井筒摩擦选择程序常用于研究射孔井段上的摩擦损失效应,摩擦损失效应对水平井是十分重要的。

专门的气井程序

Eclipse软件有二个提高气体流动模拟精度的程序。第一个是用与表皮系数相关的流动来模拟因井筒附近紊流而形成的压力降非达西指数流动。程序用保证层段上机械表皮系数是互相一致的方法计算每个完井层段的紊流表皮系数和气体流入量。为了在超长的时间步长内保持稳定,用全隐技术进行计算。当因紊流引起的压力降与总压降成正比时,这些事先警告是十分重要的。

第二个或是用地层游离气正常流入动态关系拟压力方程选择程序,或是用Russell Goodrich方程选择程序。这二个方程都能模拟气的粘度和密度变化以及网格和井筒之间压力降变化。对于Russell Goodrich方程,是用井筒和网格之间的平均压力代替网格压力值,来确定上述参数。对于拟压力方程,是用网格和井筒之间压力的积分平均值来确定上述参数。

细网格平衡程序

油藏网格区内的原始压力和原始饱和度可用平衡程序来计算。已知油藏内每一个相的静水柱压力梯度,该梯度与局部的流体密度有关。用毛细管压力表,根据各相之间局部的压力差确定原油饱和度。

通常,使每一个网格区内的流体饱和度等于网格区内中心的那一个点的饱和度。然而,当一个大的网格区含有流体接触面时或具有流体的过渡带部分时,这种方法在流体储量的估算能产生明显的误差。在这种情况下,Eclipse软件具有提高地层流体储量计算精度的功能。在网格区内,可以确定若干面上流体的饱和度,将这些值平均就可以得到网格区内总的流体饱和度。这种方法还考虑了倾斜网格的形状,用网格区内每一个面上的流

18

体饱和度与网格区内每一个水平横切面的面积进行权衡得到倾斜网格的流体平均饱和度。

水层模拟

水层模型可用解析模型和模拟模型任选一个来代表。解析模型是Fetkovich水层模型和Carter-Tracy水层模型。模拟模型是由在一个方向上按顺序连接在一起的水层网格所组成。模拟水层序列的最里面的网格能与油藏的若干个网格相连接。每一个水层网格的深度、方向(维)、孔隙度和渗透率等参数可以分别给定,这就给用户在定义每个水层时,给予充分灵活性。

辅助程序

辅助程序是Eclipse软件系统的一部分,现描述如下:

Fill——图形角点前处理程序

Fill程序根据图形角点产生特殊的网格系统。为使网格沿断层线分布,可将网格变形。用这个程序,都能很容易地产生倾斜的平移断层,复杂的剪状断层,变形楔状等特殊网格。正如程序定义所指,该程序专门设计用于内插空白数据,这些数据仅在网格的几个角点上有。为了产生完整的Eclipse网格数据文件应填入缺失的值(数据)。例如,岩石物性中的孔隙度仅需要在每个网格每个层的几个点上给定,Fill程序就能提供其余部分的数据。

VEP——井筒水动力计算前处理程序

VEP程序能计算井和管线中的压力降,并产生垂向流动动态表。这个动态表是井底压力与流量、井口压力、含水和含气比例有关的多维数据表。当模拟有井口压力限制时,须输入这些数据表。将这些表进行内插,Eclipse程序就能进行包含井口压力在内的计算。这种方法比Eclipse实际压力降计算更为效,这种方法在每个时间步长内必须做几次迭代。表的查阅方法也能使得程序能更容易的用全隐的多个变量进行井口压力的计算。

VEP程序选用一个标准的多相流动对比,在给定的综合流动条件下,对压力进行全面的计算。这些计算是指气举效应,井下泵效应,气压缩机效应和地面油咀等。

有五个多相流动对比的选择

(1)Aziz,Govier和Fogarasi;

(2)Orkiszewski;

(3)Hagedorn和Brown;

(4)Baggs和Brill;

(5)Mukherjee和Brill。

不同的多相流动对比可用于不同部分的油管计算,即不同的对比能使

19

VEP程序对水平管线,波状管线和直井油管及主管等进行计算。

EDIT——编辑程序

EDIT是一个人机对话全屏幕的编辑程序,它专用于简化Eclipse输入数据文件的编辑。在EDIT程序一开始,屏幕上显示输入文件的综述,还显示若干关键词,但这些不是它们想要的数据。然后,才有可能将注意力集中选择所要的关键词,并且EDIT将会显示合适的数据表格。例如,关键词SWEN(含水饱和度函数)表由若干的数据表组成,每个表由3个柱列数据组成(含水饱和度,水相对渗透率,水—油毛细管压力)。不同数据表用于不同的油藏网格区,倘若已使用了相对渗透率方向选择,那么不同的数据表还可用于辩别流动的方向性。当数据输入表内,EDIT程序会进行简单的数据检查,还提醒用户是否有输入的误差。在屏幕顶部的帮助行列出了每项数据有用资料。在帮助行上,输入Help(帮助)指令,就能进入Eclipse手册中,结果显示在屏幕的顶部或底部,在那儿,几个另外菜单的选择同样也显示了。在编辑关键词SWEN时,若按Help键时,则描述了SWEN关键词的内容就会显示在屏幕上。

GRAF—独立图形后处理程序

GRAF是一个容易使用,具菜单式操作性能的绘图程序。GRAF程序典型的用途有:

1.线图形显示,例如,油田的压力与时间关系图;

2.包括断层在内的油藏网格图形显示;

3.动态图形显示,如饱和度或压力随时间变化的彩色动画片显示。

在每一个报告时间,Eclipse产生的汇总文件可绘制线状图形。每一个汇总文件包含在中间时间步长中记录的与油田,井组,井,连接区和网格数目有关的数据。这样汇总文件(像restart—重新启动文件)可常于监督作业的进行。

在模拟开始时,由Eclipse产生的GRID文件能产生网格图形。GRID(网格)文件包含了油藏内每个网格角点的位置,并且使用图形软件包,以便用各种方法显示网格。

在事先规定的时间,Eclipse产生的Restart文件可以产生动态图形。Restart文件包含了每个网格详细的资料,例如饱和度和压力的资料。包含这些数据主要目的是在事先规定好的时间里,使Eclipse模拟重新启动(restart)然而,Restart文件同样可以与GRID文件连用从而产生饱和度和压力的彩色图形。图形显示软件包能让用户用一系列的彩色动画片前前后后的显示水驱、气锥等现象。

PSEUDO—三维拟函数生成程序

PSEUDO程序是Eclipse软件为了进行低费用敏感性研究而产生粗网格拟函数功能的软件包、为了能产生粗网格模型的拟函数,PSEUDO程序使用

20

了在细网格模拟期间产生的Eclipse restart文件。细网格模型可以是一条多层的横剖面,模型或者是一个大型的3D模型。粗网格模型可以是1D、2D或3D模型,它可以把相邻细网格模型合并成粗网格。因为PSEUDO程序根据Eclipse单独运算,故重新定义粗网格并产生新的经济拟函数是有可能的。由PSEUDO程序所产生的定向拟相对渗透率和相应粗网格的每个面,可以输到Eclipse粗网格模型中去,从而重新得到细网格模型的结果,或者进行低费用的敏感性分析。

口令

有些公司购买或租用有限制的已经修改过的各种程序的Eclipse软件版本。有限制的Eclipse100版本只能用Ecl提供的特别的口令来运算。Eclipse程序的限制程度现已发展到用户可要求Eclipse公司从Ecl中发行一个新的运算口令。

附加的程序

用户可以购买专门新增到基本的Eclipse100软件上的其它程序,这些程序描述如下:

Eclipse200程序

Eclipse200由一系列专门对Eclipse100扩展的程序所组成,这些扩展的程序增加了模拟的能力。它们是:

(ⅰ)聚合物驱替模型;

(ii)地面网络程序;

(iii)灵活边界条件程序;

(iv)局部网格加密和网格粗化程度;

(v)井筒摩擦选择程序;

(vi)气田作业模型;

(vii)溶剂模型;

(viii)GI拟组分模型;

(ix)表面活性剂模型。

这些专门扩展的程序按照需要可以分别购买。

Polymer Flood Model—聚合物驱替模型

聚合物选择程序是使用全隐式5组份模型(油、水、气、聚合物、盐水)来模拟研究聚合物驱替过程的详细机理。该模型解释了由于聚合物的溶解作用的增加使水相的粘度增加情况,并且还解释了由于高速非牛顿剪切效应使粘度降低的现象。聚合物与岩石表面的吸附和吸附后导致水相相对渗透率减少的现象都可以用模型来处理。有关聚合溶解流动不能进入孔隙体积的效应也包含在程序中。油藏工程师可以选择性的使用模型来研究因盐水与聚合物相互作用而使聚合物溶解粘度降低的原因。一个全隐公式用来解5相组分的联立方程,以确保在所有的时间步长中聚合物完全稳定

21

溶解。

Surface Networks—地面网络程序

地面网络选择程序是用于提供井组可变的井口压力界限,按照一套管线压力损失的关系,井口压力界限随井的产量而变化。一个简单的应用是模拟水下完井,在水下,几口生产井生产的油进入一个单管汇,该管汇用流动管线与分离器连接。管汇处在压力则取决于井的总产量,气油比和含水量。为了更复杂的应用,可以建立管线的多级分层网络。每条管线已知有一个垂向流动动态数据表(该表可用VEP程序产生),用这个表可以确定沿管线的压力损失,用同法计算井筒的压力损失。在用户定义的井段上,修改网络每个节点压力,可以进行网络的平衡计算。同样也可分别建注水管线网络,为了模拟进入注水管汇的泵动态,注水网络是十分有用的。

The Flux Boundary Option—灵活边界选择程序

灵活边界选择程序只在油田的小范围内进行模拟。研究全油田的基本运算和穿过油田小部分边界的流动时间都被编入flux文件中去。读flux文件,建立与全油田一致的小部分边界条件,接着进行油田小部分的模拟运算。这个选择程序对只有几口剩余井的历史拟合是十分有用,当仅仅是油田剩余小范围需要重新拟合时,这个程序可以减少重复的运算费用。该程序可与泄油区的程序相结合(见关键词OPTIONS,键Switch25)以减少敏感性分析的计算费用。

The Local Grid Refinement Option—局部网格加密选择程序

该程序能完善井周围网格的定义。局部网格模型可以是2D径向,3D径向,或3D笛卡儿座标系统。与整体模型相比,局部模型可以是多层的。在局部模型和整体模型之间的传导率可以自动的计算。分别对每个局部模型解能取得运算的效率。总的来说,局部模型与整体模型相比,它的时间步长更短,因而对全油田中某些井出现的局部维形效应可以重新得到解答,而不需要停止模拟的运算。粗化网格程序同样亦常于油田模拟时按区合并网格,这些合并的地区,模拟的精度显得十分重要。

The Wellbore Friction Optionmodel—井筒摩擦选择程序

井筒摩擦选择程序能模拟由于井内射孔管壁和射孔与井底参考点之间引起的压力损失效应。该程序主要用于水平井的模拟,因为在水平井的水平部分,其摩擦力的压力损失是明显的。当决策一口水平井的最佳水平段和直径时,考虑摩擦力压力损失是很重要的。为了保持数值的稳定,摩擦力的计算用全隐式进行。

The Gas Field Model—气田作业模型

气田作业模型是用于模拟气田生产作业和限制的综合性程序。该程序包括下列五个子程序:

(ⅰ)用月产量乘平均日产量或日合同是(OCQ)来设置月生产的目标

22

产量;

(ii)调正每一合同年的日合同量,使其适应调峰系数的需要;

(iii)控制和输出商品气产量;

(iv)考虑到穿过网络增加的压力损失,估算每个时间步长的输气能力;

(v)一当需要完成生产目标时,在网络上的压缩机则自动投入作业。

The Solvent Model—溶剂模型

溶剂模型是Todd和Longstaff关于四组分混相驱替经验模型的扩展模型。Eclipse软件用全隐技术解四组分(油、水、气、溶剂)的四个联立方程。若油、水、气三相系统缺少溶剂组分时,它们可减化为传统的黑油模型。溶剂相进入模型后对油藏中的烃类流体能进行混相驱动模拟。用户可用Todd—Longstaff参数来控制混相能力的大小,以及用输入混相能力与溶剂饱和度关系表来控制混相能力的变化。大多数可能的混相流体仅在高压下才是混相的。用定义一个混相能力与压力关系表Eclipse软件中就能模拟压力的作用。

The GI Pseudo—Compcsitional Model—GI拟组分模型

GI组分模型是一个四方程模型,该模型改变全组分的描述就能提高挥发油和凝析油系统的模拟。GI参数是一个与原始烃类接触的干气的测量值。在饱和压力下的Rs、Rv、Bo和Bg值都应输入Eclipse程序中,这四个值是压力和GI参数的函数。例如,气的饱和值Rv随着GI变小而减少;这样,气能包含较少量的蒸发油。这个作用就解决了黑油模型中产生的一个问题,即解决了在研究循环注气时最佳采收率的预测问题。用PVT程序根据描述油藏流体的状态方程可以产生Eclipse程序所需的数据。

The Surfactant Model—表面活性剂模型

表面活性剂模型能提高水驱的模拟效果。在水相和油相之间因增加了表面活性剂而减少了毛细管压力。毛细管压力减少,就增大了从孔隙中排油的比例,从而就降低了残余油饱和度。表面活性剂驱替选择程序模拟了残余油饱和度的减少,并且还修改有关相对渗透率与毛细管压力的关系函数(即,毛细管力与粘滞力的比值)。该模型能研究表面活性剂浓度对水粘度的作用,和表面活性剂在岩表面的吸附作用。

Eclipse300模拟程序

Eclipse300是一个具有立体状态方程,能进行与K值有关的压力处理和黑油处理的组分模拟程序。状态方程有:Redlich—Kwong,Soave—Radlich—Kwong,Peng—Robinson,和Zudkevitch—Joffe四个状态方程。Eclipse300程序与GRID,PVT,VEP和GRAF兼容。半隐式和全隐式技术都是有用的。单井的模型仍是Eclipse100用的多级分类管理系统,但另外还带有生产井向分离器和天然气工厂供气的复杂部分。根据流体的地面体积,

23

油藏体积,井底压力等可以强有力的对单井、井组、油田各级的生产量进行控制。对特定的井或井组的一部分流体实行回注是有可能的。在组分模型里,注入量可以专门规定,或减少其它井或井组的注入量。

GRID—网格程序

GRID是一个常用于Eclipse油藏建模,设计模拟网格以及产生输入数据人机交互的模拟程序。基本的地质模型由等值线,井位,图形特征等部分所组成,地质模型可由数字化输入,或直接根据工业标准图形文件输入,地质模型可有助于当作建立模拟网格的背景显示。GRID是一个代替Fill程序人机对话的绘图程序。图形特征设置包含:彩色图形,角点,正常图形,倾斜或垂直断层,面积,横剖面或3D立体显示,为与原始输入数据进行对比对网格特征轮廓进行回放,灵活使用编辑和综合的“帮助”设置。局部网格加密和网格粗化在Eclipse300上都能准备和显示,为了使用这二个设置可用Eclipse200局部网格加密选择程序。

ECLMAP—Eclipse绘图软件包

ECLMAP是一个与GRID程序相结合的人机对话绘图软件包。该程序提供了数字化输入等值线、断层、图形特征等,网格化数据抽稀,以及在网格图上进行操作等设置。它还提供了图形投影,地震炮点,测线编入,体积计算和地震反射层的深度转换或网格图的深度转换等设置。

PVT—状态方程软件包

PVT是一个人机对话状态方程(EOS)软件包,为了确定油藏流体的相态,该软件包用于实验室的资料分析。用物理平衡检验,可测试实验室数据测量的质量。利用立体状态方程的变化来模拟实验室的分析数据,并且用回归一个或多个EOS的参数使测量值和计算值之间的误差减到最小。于是EOS模型常可用于产生适用于Eclipse100,200或300所需的数据。

WELLTEST—试井程序

WELLTEST是一个人机对话的试井程序。它提供了多方面的特征图,这些图形有助于油藏工程师选择油藏建模时所需的最合适数据。为了确定模型参数的第一个计算值,通过优化来交互的拟合直线。于是,优化能确定油藏的参数值,用这些参数值,在选择的模型和实际观察的数据之间能观察的数据之间能给予最好的拟合。该程序同样还包含了典型曲线分析,气井试井分析,和试井报告产生等设置。

CLAN—测井分析软件包

CLAN是一个人机对话的测井分析软件包,该程序在确定解释模型和使用资料上给予最大的灵活性。该软件具有输入编辑,标准化,电测井分析,泥浆录井和取心数据分析等功能。

24

2

综 述

数据文件综述

Eclipse数据输入文件可分成几个部分,每一个部分由一个关键字词引导。各个部分文头上的关键字列表如下,每个部分与关键字一起,带有一个简要的内容描述。每部分更详细的分类内容详见每部分的综述。

关 键 字

内 容(简述)

RUNSPEC(必须)

GRID(必须)

题目,作业维数,运算键,目前的相态等。

计算网格的图形规定(网格角点位置)和每个网格岩石物性的规定(孔隙度,绝对渗透率等)。

EDIT(选择) 对计算的孔隙体积,网格中心点深度和传导率进行修改。

PROPS(必须) 油层岩石和流体性质(密度、粘度、相对渗透率、毛管压力等)随压力、饱和度、组份变化表。

REGIONS(选择) 为计算流体物性(PVT,即流体的密度,和粘度)饱和度特性(相对渗透率和毛细管压力)原始条件(平衡)压力和平衡饱和度),流体储量(流体储量和区内流动)所需对计算网格分区。

如果这部分省略,所有的网格区块都放到第1区。

SOLUTION(必须) 对油藏原始条件的确定:a)用规定的流体接触面深度到可能的流体高度(平衡)计算;b)从上一次运算建立的重启文件中读出;或c)用户自己确定每一个网格区的原始条件(一般不用此项)。

SUMMARY(选择) 在每一个时间步长后,确定的数据被编入汇总文件中。如果必要的话在每次运算结束后输出某种类型成果的图形(如,含水量与时间关系曲线)。若这部分省略,则不产生汇总文件。

SCHEDULE(必须) 确定模拟的作业(产量、注水量控制和限制)和给定需要输出模拟结果的时间。在Schedule部分中同样还能确定垂向流动的动态曲线和油管模拟参数。

这几个部分的关键字必须按上面显示的次序排列。这样的排列则有助于对这几部分的记忆(如,Ruth Greets Edgr Properly at REGgies Social

Summer School)。

注意每部分文头的关键字(像所有其他的文头关键字一样)必须从第1列开始。所有的字符在1—8格内有效。该行上的其余字符(从第9列开

25

始)当作为注解。

所有行开始用二个字符“— —”表示,作为注解处理,而在Eclipse模拟中是不参加运算的。注解行(和空白行同样)可插入数据文件中的任一个地方。注解还可以加到用二个字符“— —”开始的数据行的末端,但是,在这种情况下,注解不能含有任何引号。在没有二个字符“— —”的情况下,注解也可以放在同一行的斜线符号(/)之后面,斜线符号(/)常用来表示结束数据的记录。(然而,用EDIT程序准备数据文件时,注解仅能立即插入每个关键字的前面。在关键字部分的数据里,EDIT决不能插入注解)。

在下列一个关键字的数据里,星号(*)可用来表示“重复计数”。在一个数前加一需要次数和*号,就能表示这一个数需要重复计数的次数。例如:

**3450 400 2350

等于

450 450 450 400 350 350

在*号的二边不能插入空的空间。

对具有内装缺省值功能的软件来说,有些项目的数据可以缺省。当要使用缺省值时,关键字的描述会提示。有二种方法来设置它们的缺省值。第一个方法是用一个斜线(/)过早的结束一个数据的记录,其余没有被给*定的数将作为它们的缺省值。第二个方法是在斜线(/)之前,加上n来表*示缺省值,这儿的n是表示连续缺省的值。例如,在关键字数据中,3使其后面的三个数是它们的缺省值。在连续数和*号之间没有空白的插入空*间。如果只有一项缺省一次,那么必须加入1。只写星号是不行的。

在数据部分出现的一些关键字列表如下,与这些关键字一起还有它们功能的简单描述。更详细的论述可查阅手册中关键字部分。

关 键 字

COLUMS

DEBUG

ECHO

EXTRAPMS

INCLUDE

MESSAGES

NOECHO

NOWARN

OPTIONS

简 述

为了能读出数据输入文件,对左右二边重新定位键。

输出控制键调试设置。

每一个运算开始时数据重新打印键。

要求外推PVT和VEP数据表信息预告。

插入规定文件的目录。

所有正规信息打印和停止限止重放键。

每一个运算开始的不进行数据重新打印。

Eclipse信息预告抑止键。

在特别的程序中设置特殊活动选择键。

26

建议将不经常变动各关键字部分的正文变成个别文件,并有关键字INCLUDE将这些文件输入数据文件中。

例:

RUNSPEC

GRID

RPTGRID

*60 11*

网格是不经常改变的,所以可将它放入INCLUDE数据文件中。各个分别文件插入Include文件中的具体格式见2.4页(Eclipse100 91A Reference

Manual)。

数据文件实例:

RUNSPEC

.

.

.

GRID

*===============================================================

===============================================================

*RPTGRID 60 111 /

— — GRID data is not often changed - so put it in an INCLUDE file

INCLUDE

‘’ /

PROPS

RPTPROPS

201

INCLUDE

‘’ /

REGIONS

RPTREGS

200

INCLUDE

‘’

SOLUTION

RPTSOL

201

INCLUDE

***===============================================================

/

=============================================================

/

/

============================================================

/

27

‘’

SUMMARY

INCLUDE

/

============================================================

‘SYSTEM_DEFAULT:’

SCHEDULE

RPTSCHED

0 1 0 1 1 0

/

============================================================

32

*/

INCLUDE

‘’ /

END

============================================================

2.1 RUNSPEC部分

综 述

SUNSPEC部分是Eclipse数据输入文件的第一个部分。文件中包括了运算题目,开始日期,单位,各种作业的维数(网格区数、井数、表数等)目前相标志和选择键。只有注解和关键字INCLUDE ECHO,NOECHO,COLUMNS

OPTIONS和DEBUG放在运行定义卡的前面。除非从装有Runspec数据文件的存盘中重新运算,运行定义卡这部分始终是存在的。

运行定义卡部分有一个固定的结构。它必须占30道记录或者是30行。第一行是Runspec关键字的名称,它开始于第1列。所有下面的奇数行是指导用户的注释,在Eclipse软件中是不参加计算的。偶数行是混合数和局部的程序键。每个数据行必须用斜线字符(/)结束。除了题目,网格区数和单位之外,所有的数据项已设置了缺省值(缺省值用角形括号表示,下面将会细述)。如果一个数据行过早的用(/)符号结束,那么其余的数据行则设置为缺省值。对大多数计算来说,主要的数据行都可留下它们的缺省值。

运行定义卡中所必须的最少数据行是:

Record 2 题目

Record 4 在X、Y、Z方向上的网格区块数。

Record 6 模拟的相态,即各相的饱和度,Rv或Rs变化值。

Record 8 单位

Record 18 最大生产井数、井组数、井组内生产井数以及单井射孔井段。

Record 28 模拟的开始日期。

Record 30 模拟运算(缺省)或数据检验。

28

还可以缺省的其他项目是:笛卡儿座标图形,不相邻连接的分散流动模型,一套PVT,饱和度和平衡数据表,报告区块和各个格式不正规重新启动和绘图文件。

在表数/网格数和最大的生产井数/节点数之间进行区别是十分重要的,前者指在数据文件中必须提供正确的项目数,后者则仅指能够给定项目数的上限。

至少,一开始,用户可以用EDIT编辑程序找到取简单建立Runspec部分的数据文件。数据文件描述的记忆方法和注释,以及下面的例子,这些都可以在EDIT编辑程序中产生。

Runspec的详细描述

Record 1 在1-8列,应输入关键字Runspec 。

Record⨪2⨪⨪应包含运算题目(限制为72个字符)。它不需加号。该行也需要用(/)斜线结束。

Record 3 本行以及下面的各奇数行(5、7、9…29)在Eclipse软件是不参加计算的,它们仅作为注释。但要注意,在这些行的第一个字符不应当是一个字母,而用EDIT文件说明时每行开头用一个关键字的字母。

Record 4 包括网格的尺寸,径向/笛卡儿开关键,不相邻连接键,水层数值信息,双孔隙度运算键。

〈error〉 NDIVIX—X方向上的网格数(或R—径向)

〈error〉 NDIVIY—Y方向上的网格数(或θ角)

〈error〉 NDIVIZ—Z方向上的网格数

选择:径向,NNC,水层数值分析,双孔隙度。

〈F〉QRDIAL—T为径向图形,F为笛卡儿座标图形(无引号)。

〈I〉 NOMRES—研究的油层数,除了用几套坐标线的角点图形之外,该部分应设置(1见GRID网格部分中的COORD和ZCORN关键字)。

〈F〉 QNNCON—若是不相邻连接则设置T。常数于水层数值,断层的角点图形,双孔隙度,完成径向流动研究的园周,或者直接用NNC关键字来设置。(GRID部分)。

〈O〉 MXNAQN—最多的模拟水层数(见GRID部分中的AQUCON和AQUNUM)。

〈O〉 MXNAQC—在关键字AQUCON中,水层最多的连接线(见GRID部分的AQUCON和AQUNUM)。

〈F〉 QDPORO—T为双孔隙度运算。在这种情况下,NDIVIZ必须是偶数,第一个NDIVIZ/2层将当作为骨架网格,其余的则作为裂缝网格。裂缝网格渗透率和裂缝网格孔隙度相

29

乘,骨架和裂缝总成对存在,并且用SIGMA关键字来决定裂缝网格的渗透率。若QDPORO是T,则QNNCON同样应为T。

〈F〉 QDPERM—在双孔隙度运算中,T为骨架与骨架之间的传导率,不等于0(双孔隙度/双渗透率)。F为骨架与骨架之间的传导率,等于0。若QDPORO为F,则QDPERM则没有作用。

Record 6

包括相存在的标定(相存在为T,相不存在则为F),以及流体示踪选择键。

相存在:

〈F〉 OIL T或F

〈F〉 WATER T或F

〈F〉 GAS T或F

〈F〉 活油中的溶解气 T或F

(若油和气都存在,则只能为T)

〈F〉 湿气中的蒸发油 T或F

(若油和气都存在,则只能为T)

若油有一个恒定的,均匀的溶解气浓度,并且油层压力不再降到泡点以下,就用F设置气和溶解气的标志,则能更有效的进行模拟,把油当成死油看待,并且用关键字RSCONST或RSCONSTT设置气油比Rs值为常数。运算的作业当成没有气的死油模拟,但是在计算时则考虑了Rs值(溶解气浓度为常数)。

若气有一个恒定的,均匀的挥发油浓度,其气层的压力不同志降到露点压力之下,就用F设置油和蒸发油的标志,这时把气当成干气看待,并且用关键词RVCONST或RVCONSTT设置油气比Rv值为常数。运算作业当成没有油的干气模拟,但是,在计算时则考虑了Rv值(蒸发油浓度为常数)。

选择:油比重(API),盐水,示踪剂,混相的(200:聚合物,GI溶剂表面活性剂),

〈F〉QAPITR—T为需要油比重(API)的跟踪选择(见附录中的API追踪)。这个选择可用在活油或死油系统中,但它不能用在湿气系统中。RSCONST(T)关键词绝不能用于API追踪作业中。

〈F〉QPOLYM—T为需要进行聚合物驱替选择。这是Eclipse200中的一个选择程序,用户可参考Eclipse200参考手册。

〈F〉QBRINE—T为需要进行盐水选择(见附录中盐水追踪)。对应用聚合物驱替选择(QPOLYM=T)的运算,代替了QBRINE

30

控制聚合物—盐的敏感性分析选择(见Eclipse200参考手册)。

〈O〉NOTRAC—被注入示踪油的数量(见附录中的示踪剂追踪部分)。

〈O〉NWTRAC—被注入示踪水的数量。

〈O〉NGTRAC—被注入示踪气的数量。

〈F〉QMISCL—T为需要混相气驱替选择。若QMISCL=T,则不应使用垂向平衡选择(见第12记录道QVEOPT)。

〈F〉QGGICO—若需要模拟气藏凝析油/挥发油系统的GI模型,则设置T。这是Eclipse中选择程序;用所参考Eclipse200参考手册。

〈F〉 QSOLVT—若需要4组分溶剂模型,则设置T。这是一个Eclipse200的选择程序;用户可参考Eclipse200参考手册。

〈F〉QSACT—如果需要表面活性剂模型,则设置T。这是一个Eclipse200的选择程序;用户可参考Eclipse200参考手册。

Record 8 包括运算中所有数据的单位换算—除非你用EDIT程序,否则你就辅以引号(‘METRC’,‘Field’或‘LAB’)。

Record 10 包括压力表的维数。注意把活油(油中含有溶解气)的PVT性质当成2D数据表输入—即在不同Rs下Bo、Vo与压力关系数据表。同样,湿气的(气中含有挥发油)PVT性质当成2D数据表输入即在不同压力下Bg、rg与Rv的关系数据表。

〈20〉NRPVT—在活油PVT表中Rs值的最多节点数(PVTO,PVCO)或在湿气PVT表中,Rv值的节点数(PVTG)。如果没有活油或湿气,就设置1。

〈20〉NPPVT—在任何PVT表中压力值的最多节点数(PVDG,PVDO,PVTO,或PVTG)以及岩石压实性表(ROCKTAB)。

〈1〉 NTPVT—PVT表数(不同的PVT表可用于油藏的不同部分—见模拟区部分中的PVTNUM和DENSITY,PVDG,PVDO,PVTG,PVTO,PVCO,及油藏特性部分中的ROCK)。

选择:岩石压缩性,失压垂向平衡,毛细管压力/表面张力函数,分子扩散。

〈I〉NTROCC—岩石压缩性表数(不同岩石压缩性表可用于油藏的不同部分—见模拟区部分中的ROCKNUM和油藏特性部分中的ROCKTAB)。若岩石压缩性选择键没有打开,那么NTROCC不参与运算(见下面的QROCKC)。

〈F〉QROCKC—若使用岩石压缩选择程序,则设置T(用油藏特性部分

31

中的关键词ROCKTAB输入传导率与同压力有关的孔隙体积的乘积)。

〈T〉QRCREV—岩石压缩可逆键。若随着压力的增加,其压缩是可逆的,则设置T。F则为不可逆的,即当压力增加时,孔隙空间决不会膨胀。若QRCREY为F,则QRCREY不参加工作。

〈F〉QROCKH—T为使用岩石压缩滞后选择程序。与滞后压力有关的孔隙体积和传导率的乘积用油藏特性部分中的ROCKTABH关键词输入运算定义卡中。若QROCKH为T,则可逆键不参加作业。

〈F〉QVEOZD—T为使用垂向的压缩平衡选择程序。为了减少计算机的费用,将3D垂向平衡当成2D平面运算。仅使用于QVEOPT为T(Record12)和穿过油田垂向平衡是有效的—例如没有泥岩隔层引起垂向平衡的破坏(无效)。

〈F〉QSURFT—若油—气或油—水毛细管压力随表面张力而变化,则设置T。用关键词STOG和STOW输入表面张力表。

〈F〉QMDIFF—T为需要用分子扩散选择。根据网格数据来计算扩散率。这些计算值可以用DIFFMX,DIFFMY,DIFFMZ关键词来修改,或者用上述三键将它们输入EDIT文件中。用油藏特性部分中的DIFFC关键词输入扩散系数。(见附录分子扩散部分)。

Record 12 包括了饱和度表的维数,方向性相对渗透率的标志,垂向平衡,滞后效应,饱和度表端点刻度和混相驱选择程序。

〈20〉NSSFUN—在任何饱和度表上饱和度值的最多节点数(见油藏特性部分中的SGFN,SGOF,SLGOF,SOF2,SOF3,SWFN和SWOF)。

〈1〉NTSFUN—用于油藏不同部分的饱和度表数(见模拟区部分中的

SATNUM和油藏特性部分中的SGFN等)。

选择:方向相对渗透率(Kr),VE(垂向平衡),滞后效应,端点刻度,混相驱。

〈F〉QDIRKR—如果使用方向相对渗透率数据表设置T,要不然使用F(无引号)—(见模拟区部分中的KRNUMX,KRNUMY和KRNUMZ)。

〈T〉QREVER—如果方向相对渗透率是可逆的设置T(即,使用了相同的饱和度表流体是否从I→I+1或从I→I-1流动)—要不然使用F(无引号)。若QREVER是F,那末用KRNUMX,KRNUMX-,KRNUMY,KRNUMY,KRNUMZ,KRNUME键必须建

32

立6个方向的相对渗透率区(见模拟区部分)。若QREVER是T,QDIRKR是T,那末需设置3个方向相对渗透率区(可用KRNUMX,KRNUMY,KRNUMZ键)。注意,若QDIRKR是F,那么QREVER必须是T。

〈F〉QVEOPT—T为需要垂向平衡选择程序。

〈F〉QHYSTR—T为需要滞后效应选择程序(单个饱和度函数表可用于驱替和吸吮)—否则设置F。如果设T,则用模拟区部分中IMBNUM键设置分区吸吮表。若QDIRKR也设置T,那么用IMBNUX,IMBNUY和IMBNUZ来设置不同方向的吸吮表(若QREVKR是F,则还加上IMBNUMX-,IMBNUMY-,IMBNUMZ-)

〈F〉QSCALE—T为使用饱和度的端点刻度选择程序—否则设F。用PROPS(油藏特性)部分中的SWL,SWCR,SWU,SGL,SGCR,SGU,SOWCR,SOGCR中的一个或是用PROPS部分中的ENPTVD和模拟区部分中的ENDNUM可输入端点饱和度值。用KRW,KRG,KRO或ENKRVD关键字同样亦能刻度端点相对渗透率。

〈F〉QSCDIR—T为饱和度表的端点刻度选择是方向的(即,不同饱和度表的端点值可用于流入X、Y、Z方向的饱和度表)。每当QSCALE是T,则QSCDIR仅能是T。端点饱和度即可用 PROPS部分中的SWL,SWCR,SWU,SGL,SGCR,SGU,SOWCR,SOGCR关键字或用PROPS部分中的ENPTVDX,ENPTVDY,ENPTVDZ,和模拟区部分中的ENDNUM关键字输入合适的方向。若端点的相对渗透率也被刻度,则必须输入KRW,KRG,KRO或ENKRVD关键字的方向形式。

〈T〉QSCREV—T为饱和度表的端点刻度是可逆的(即,使用的饱和度端点的流动是否从I→I+1或从I→I-1)—否则设F。每当QSCDIR为T,则QSCREV仅能是F。无论是在合适的方向(+X,-X,+Y,-Y,+Z,-Z)F,用油藏特性部分中的SWL,SWCR,SWU,SGL,SGCR,SGU,SOWCR,SOGCR关键字还是油藏特性部分中的关键字ENPTVDX,ENPTVDX-,ENPTVDY,ENPTVDY-,ENPTVDZ,ENPTVDZ-和模拟区部分中的ENDNUM都可以输入不可逆的端点饱和度。若QSCREV=F,而端点相对渗透率同样也刻度了,那就必须输入关键字KRW,KRG,KRC,或ENKRVD+Ve和-Ve的方向形式。

〈20〉NSENDP—与深度表对应的任何饱和度表的端点的最多节点数(见

33

PROPS部分中的关键字ENPTVD,ENKRVD)。

〈1〉NTENDP—与深度表对应的饱和度端点的最多表数(见PROPS部分中的关键字ENPTVD,ENKRVD)。

〈F〉QTWOPT—当使用相驱替选择时,(QMISCL=T),则上游二点算法常可用于计算烃类组分的相对渗透率,这时应设置T。更详细的见附件中的混相驱替模拟。

〈20〉NSMISC—与含水饱和度对应的任何混相驱残余油饱和度表内的饱和度值的最多节点数(见PROPS部分中的关键字SORWMIS)。这个参数仅用于QMISCL=T。

〈1〉NTMISC—与含水饱和度对应的混相驱残余油饱和度的最多表数和最多混合驱参数的区块数。(见PROPS部分中的TLMIXPAR和SORWMIS)。这个参数仅用于QMISCL=T。

Record 14 包插控制平衡计算的参数和维数(见模拟解部分)。

〈20〉NDRXVD—在任何RVVD或RSVD表内深度的最多节点数(随深度变化的原始Rs、Rv)。

〈1〉NTEQUL—设置的平衡数据数(不同的平衡数据可用于油藏的不同部分—见Reglons部分中的EQLNUM和SOLUTION部分中的RSVD,RVVD,EQUIL)。

〈100〉NDPRVD—由平衡码内部建立的与深度对应的相压力深度表内深度的最多节点数。

选择:静止,门槛压力,可动油,示踪剂。

〈F〉QUIESC—如果为了达到原始静止状态需要对压力进行修正,设置T。(即,产生一个真正的稳定状态解)。一个块中心点平衡(关键字EQUIL中第9项设置0)将终是静止的,但是更精确的选择将不是绝对静止的。若QUIESC设为T,为了得到初始静止状态解,将对相压力进行修正。在运算期间,继续进行相压力的修改。若使用垂向平衡选择程序,(第12记录道,QVEOPT)则决不能用QUIESC选择程序。

〈F〉QTHPRS—T为能进行门槛压力选择。门槛力阻挡了在不同平衡区之间发生的流动一直到可能的压力差超过了门槛压力值。可用关键字THPRES给定门槛压力值,或者从初始解中进行缺省计算。

〈T〉QREVTH—若假设门槛压力在每个方向上的流动是相等的,则设置T。若设置F,则在平衡区之间流入各个方向上的门槛压力是不同的。若QTHPRS设为T,则QREVTH也仅相应设置为T。

34

〈F〉QMOBIL—若设T,则要进行原始可动流体端点临界饱和度校正。这个选择适用于小规模平衡(见关键字EQUIL)计算。它不能用于VE运算(见附录中饱和度表的刻度部分)。

〈1〉NTTRVD—与深度对应的原始示踪剂浓度表数(见Solution部分TVDP关键字)。

〈20〉NSTRVD—在与深度对应的任何原始示踪剂浓度表内深度最多节点数(见Solution部分TVDP关键字)。

2.2 GRID部分

综 述

本部分定义了计算网格的图形形状和每个网格内各种岩石性质(孔隙度、绝对渗透率、净毛比)。根据这些资料,ECLIPSE软件就可以计算网格的孔隙体积、中部深度和网格内部传导率。

无论用笛卡儿坐标X、Y、Z、还是用径向坐标R、θ、Z油藏的图形形状都可以在RUNSPEC(运行定义)记录卡4中以两种方法都可设置:

(a)块中心图形网格:网格是水平的,8个角都是直角,每个网格用它的3个面大小和顶面深度来定义。

(b)角点图形网格:网格8个角的位置是独立设置的,不要求是直角。

ECLIPSE软件用关键字COORD和ZCORN来表示角点网格图形。其他规定网格尺寸的关键词是指块中的图形网格。通常,角点图形网格输入关键字NEWTRAN,就能在记录4中把不相邻网格连结起来。

模型中不同方向(NX、NY、NZ)的网格数在RUNSPEC记录4中设置。注意:X和R方向用I表示,Y和θ方向用J表示,Z方向用K表示。

所有的深度和厚度沿着垂直方向的Z轴量度,用一个较大值表示一个较深的深度。笛卡儿坐标的原点是左后角的顶点。X轴上的坐标值是从左向右增加的,Y轴上的坐标值是从后向前增加的。径向坐标的原点通常是1口井模型的中心。油藏的内半径必须用关键字INRAD设置。R与X一致,θ与Y一致。当θ=0时,R向右增加;θ值向顺时针方向增加。

输入资料必须按I、J、K顺序进行。以便从一页的左顶角沿着行读I,另一行读J、另一页读另一个层K的顺序读资料。

每个网格的形状、大小、总体积和位置都是由块中心网格和角点网格确定的。网格的其它特性是不考虑网格几何形态的平均值。当每个方向上多于一个网格时,但每个网格在不同方向上须有一个孔隙度和渗透率值。

资料可以用每个网格的NX、NY、NZ值按数组形式输入,或对部分网格区的所有网格用关键字BOX定义。在一个新的BOX被 定义之前,一个

35

BOX一直存在,或由关键字ENDBOX来读得。缺省BOX是整体模型。面积的大小、X、Y、R、θ值可以按照矢量输入,例如DXV(I),这里仅用I变化DX。

对于常数,可以用关键字EQUALS设置,或用关键字ADD、MULTIPLY、COPY处理。两者任选一种。

在块中心图形网格中,每个网格的长度用下表关键字说明。在径向图形网格中,可以用DR代替关键字OUTRAD,在径向上,ECLIPSE软件将自动产生对数空间网格。每个网格的顶部深度用关键字TOPS说明。

基 本 关 键 字 表

块中心图形网格

笛卡儿 径向 角点图形网格

INRAD

DXV或DX DRV或DR或OUTRAD COORD

DYV或DY DTHETAV或DTHETA ZCORN

DZ DZ

TOPS TOPS

PORO PORO PORO

PERMX PERMR PERMX/R

PERMY PERMTHT PERMY/THT

PERMZ PERMZ PERMZ

ECLIPSE软件定义了传导率,TX(I)表示(网格I)和(I+1)之间的传导率。传导率可以用关键MULTX、MULTY、MULTZ、MULTR或MULTTHT(传导率乘数)来修饰。

如果关键字NEWTRAN被设置了,ECLIPSE软件就能用角点数据计算横穿断层的传导率,并在需要处自动产生不相邻连结。

传导率进一步详细说明见传导率计算附录。

角点图形网格数据通常是大量的,推荐用FILL和GRID进行前处理构成它,注意:块中心图形网格和角点图形网格不能混用。

所有关键字必须从第1列开始,所有字符1到8列有效,这行剩余部分(从第9列起)由编辑程序使用。

网格部分中可能出现的所有关键字与其功能的简要描述列表如下,更详细的描述可在关键字手册中查到。

块中心图形网格

TOPS

(必须用某种方法设置)因为

Way是单数,若多数则:Seller-

Al ways几种方法!

36

目前输入盒子网格的

顶面深度。

DZ (必须用某种方法设置) 目前输入盒子Z方向

(垂直方向)的网格尺寸。

笛卡儿块中心图形网格

DX

DXV

DY

DYV

(必须用某种方法设置)

(必须用某种方法设置)

(DY的简化版本)

目前输入盒子X方向的网

格尺寸。

或X方向网格尺寸的矢量。

目前输入盒子Y方向的网

格尺寸。

或Y方向网格尺寸的矢量。

径向块中心图形网格

INRAD (径向块中心图形网格要求) 设置径向网格的内半径。

DR (必须用某种方法设置) 目前输入盒子R方向的网

格尺寸。

DRV (DR的简化版本) R方向上网格尺寸的矢量。

OUTRAD (径向网格任意选择) 设置径向网格外半径。

DTHETA (必须用某种方法设置) 目前输入盒子θ方向的网

格尺寸。

DTHETAV (Dθ的简化版本) θ方向网格尺寸的矢量。

角点图形网格

COORD (角点图形网格要求) 定义网格范围内每个油层

和包括所有网格角点(I、J)的范围。

网格角点的深度。

有关网格范围内每个油层

坐标系统的资料。指定了

角点图形网格和块中心图

形网格以循环方式完成。

ZCORN

COORDSYS

PORO

DZNET

(角点图形网格要求)

(任意选择)

所有图形网格的参数

(必须用某种方法设置) 目前输入盒子内网格的孔隙度。

(有效厚度或净毛比任意选择) 网格的有效厚度。

37

NTG 目前输入盒内网格的净毛比。

ACTNUM (任意选择) 识别活动网格。

MINPV (任意选择) 每个活动网格的最小孔隙体积。

MINPVV (任意选择) 目前盒内网格的最小孔隙体积。

PERMR (如果径向NR>1,必须 目前输入盒内R方向上渗透率值。

用某种方法设置)

PERMTHT (如果径向Nθ>1,必须 目前输入盒内θ方向上渗透率值。

用某种方法设置)

PERMX (如果笛卡儿座标NX>1, 目前输入盒内X方向上渗透率值。

必须用某种方法设置)

PERMY (如果笛卡儿座标NY>1, 目前输入盒内Y方向上渗透率值。

必须用某种方法设置)

PERMZ (任何系统NZ>1,必须 目前输入盒内z方向上渗透率值。

用某种方法设置)

传导率关键字

NNC

NEWTRAN

OLDTRAN

OLDTRANR

(任意选择)

(任意选择)

(任意选择)

(任意选择)

不相邻连接显式通道。

说明计算传导率的角点网格数据。

说明计算传导率的块中心网格数据。

选择网格特征计算传导率。

传导率修饰

MULTR (如果径向NR>1,任意选择)

MULTTHT (如果径向Nθ>1,任意选择)

MULTX (如果笛卡儿座标NX>1,任

意选择)

MULTY (如果笛卡儿座标NY>1,任

意选择)

MULTZ (如果NZ>1,任意选择)

目前输入盒内R方向上传

导率乘数。

目前输入盒内θ方向上传

导率乘数。

目前输入盒内X方向上传

导率乘数。

目前输入盒内Y方向上传

导率乘数。

目前输入盒内Z方向上传

导率乘数。

确定含水层范围

AQUCON

AQUNUM

(定义含水层要求)

(定义含水层要求)

定义含水层与油层连通。

定义含水层范围。

38

双孔隙度

用户仅能输入骨架网格数据。

全部网格中,典型骨架网格的

垂向长度。

DZMTRXV (仅选用于双孔隙度运行) 目前盒内典型骨架网格垂向长

度。

NODPPM (仅选用于双孔隙度运行) 规定裂缝网格内渗透率不能乘

裂缝孔隙度。

SIGMA 全部网格中,骨架与裂缝的连

结系数。⨪⨪⨪⨪⨪⨪⨪⨪⨪

SIGMAV⨪(仅选用于双孔隙度运行)目前盒子内,骨架与裂缝的连结系数。

DPGRID (仅选用于双孔隙度运行)

DZMTRX (仅选用于双孔隙度运行)

扩散率分布(仅用分子扩散选择)

DIFFMR (如果径向NR>1,任意选择)目前输入盒内R方向上扩散率乘

数。

DIFFMTHT(如果径向Nθ>1,任意选择)目前输入盒内θ方向上扩散率乘

DIFF (如果笛卡儿NX>1,任意选择) 目前输入盒内X方向上扩

散率乘数。

散率乘数。

DIFFMY (如果笛卡儿NY>1,任意选择) 目前输入盒内Y方向上扩

散率乘数。

DIFFMZ (如果NZ>1,任意选择) 目前输入盒内Z方向上扩散率

乘数。

DIFFMMF (仅选用于双孔隙度运行) 目前输入盒内,骨架与裂缝扩

散率的乘数。

操作关键字

ADD (任意选择)

BOUNDARY (任意选择)

BOX (任意选择)

COPY (任意选择)

ENDBOX (任意选择)

EQUALS (任意选择)

INIT (任意选择)

在目前输入盒内,对指定的数组加个指定常

数。

定义输出图形需要的网格范围。

重新定义目前输入盒子。

在目前盒子中,把资料从一个数组拷贝到另

一个数组。

重新定义含有整体网格的目前输入盒子。

在目前盒子中对指定的常数设置指定的数组

产生GRAF和PSEUDO原始文件。

39

GRIDFILE (任意选择) 控制网格图形文件目录。(为了图形输出)

MULTIPLY (任意选择) 在目前盒内,用指定常数乘以指定数组。

NOGGF (任意选择) 说明产生没有网格图形文件(为了图形输出)

PSEUDO (任意选择) 打开拟函数软件包数据的输出键。

RPTGRID (任意选择) 设置网格数据的输出格式。

SPECGRID (任意选择) 说明网格尺寸。

关 键 字

ACTNUM

ADD

AQUCON

AQUNUM

BOUNDARY

BOX

COORD

COORDSYS

COPY

DIFFMMF

DIFFMR

DIFFMTHT

DIFFMX

DIFFMY

DIFFMZ

DPGRID

DR

DRV

DIHETA

DTHETAV

DX

DXV

DY

DYV

DZ

关 键 字 按 字 母 顺 序 列 表

简 单 描 述

识别活动网格。

在目前盒子内,给指定数组加一个指定常数。

定义含水层和油层间的连通。

定义含水层范围。

定义输出图形需要的网格范围。

重新定义目前输入盒子。

定义网格范围内,每个油层和包括所有网格角点

(I、J)的范围。

指定角点和块中心两种图形网格按循环方式完成

。如果网格内有多个油层,对网格进行分层。

在目前盒子内,把资料从一个数组拷贝到另一个数组。

目前输入盒内骨架与裂缝扩散率的乘数。

目前输入盒内R方向上扩散率乘数。

目前输入盒内θ方向上扩散率乘数。

目前输入盒内X方向上扩散率乘数。

目前输入盒内Y方向上扩散率乘数。

目前输入盒内Z方向上扩散率乘数。

用户仅能输入骨架网格数据。

目前盒子内R方向上的网格尺寸。

R方向上网格尺寸的矢量。

目前盒子θ方向上的网格尺寸。

θ方向上网格尺寸的矢量。

目前盒子X方向上的网格尺寸。

X方向上网格尺寸的矢量。

目前盒子内Y方向上的网格尺寸。

Y方向上网格尺寸的矢量。

目前盒子内Z方向上网格尺寸。

40

DZMTRX

DZMTRXV

DZNET

ENDBOX

EQUALS

GRIDFILS

INIT

INRAD

MINPV

MINPVV

MULTIPLY

MULTR

MULTTHT

MULTX

MULTY

MULTZ

NEWTRAN

NNC

NOGGF

NODPPM

NTG

OLDTRAN

OLDTRANR

OUTRAD

PERMR

PERMTHT

PERMX

PERMY

PERMZ

PORO

PSEUDOS

RPTGRID

SAVE

SIGMA

SIGMAV

整体网格中,典型基质网格的垂向长度。

目前盒子内,典型骨架网格的垂向长度。

目前盒子内,网格的有效厚度。

重新定义含有全部网格的目前盒子。

在目前盒子内,给指定常数设置指定数组。

控制网格图形的文件目录(为了图形输出)。

产生GRAF和PSEUDO原始文件。

设置径向网格内半径。

设置每个有效网格的最小孔隙体积。

在目前盒子内,设置每个有效网格的最小孔隙体积的

矢量。

在目前盒子中,用指定常数乘以指定数组。

目前盒子内,R方向上传导率的乘数。

目前盒子内,θ方向上传导率的乘数。

目前盒子内,X方向上传导率的乘数。

目前盒子内,Y方向上传导率的乘数。

目前盒子内,Z方向上传导率的乘数。

说明用网格间接触面计算的传导率。角点图形网格缺

省。不用于块中心图形网格。

不相邻连结显式通道。

停止产生网格图形文件(为了图形输出)。

规定裂缝网格的渗透率不能乘裂缝孔隙度。

目前盒子内网格的净毛比。

规定块中心图形网格数据常用于计算传导率。

选用块中心资料计算传导率。

设置径向网格的外半径。

目前盒子内,R方向上的渗透率。

目前盒子内,θ方向上的渗透率。

目前盒子内,X方向上的渗透率。

目前盒子内,Y方向上的渗透率。

目前盒子内,Z方向上的渗透率。

目前盒子内,网格的孔隙度。

打开拟函数软件包数据的输出键。

设置网格数据的输出格式。

为了容易重新启动需要产生一个SAVE文件。

整体网格中,骨架~裂缝的连结系数。

目前盒子内骨架~裂缝的连结系数。

41

SPECGRID 网格大小和字符的规定。

TOPS 目前盒子内,网格的顶面深度。

ZCORN 网格的角点深度。

数据设置实例:

2.3 EDIT 部分

综 述

本部分包括了修改孔隙体积、块中心图形网格的深度、传导率、扩散率(分子扩散选择)和用ECLIPSE软件,用输入网格部分的资料计算不相邻连结(NNC),这是整体选择。

网格输出数组关键字:

PORV DEPTH TRANX TRANR DIFFX DIFFR

TRANY TRANTHT DIFFY DIFFTHT

TRANZ DIFFZ

无论对整个油层,还是用关键字BOX定义的一组网格都可以用EDIT程序来编辑资料。象网格部分一样,可以选择控制关键字准确修改数组。

BOX EQUALS ADD

ENDBOX COPY MULTIPLY

修饰不相邻连结

关键字EDITNNC可直接用于输入不相邻连接的修饰或在网格部分自动产生的不相邻连结。

注意:关键字TRANX等是不能改变不相邻连结的传导率。但在网格部分用MULTX等关键字就能改变相邻网格间以及沿着网格正向自动产生的不相邻连结的传导率。同样地,当使用分子扩散关键字DIFFX等时,不能改变不相邻连结的扩散率,而在网格部分使用DIFFMX等就能改变沿网格正向自动产生的不相邻连结的扩散率。

孔隙体积、块中心深度、传导率、扩散率和不相邻连结显式通道(NNCS)的输出由网格部分中的关键字RPTGRID开关控制。不过,每次任何有效改变的打印数值都是在编辑部分进行的。

本节中可能出现的所有关键字及其功能简要描述列表如下,更详细说明可在关键字手册中查阅。

所有的关键字必须从一行的第1列开始,所有字符从1到第8列是有效的,这行余下部分(从第9列起)由编辑程序使用。

关 键 字 简 要 描 述

ADD 在目前输入盒子中,给指定数组加指定常数。

BOUNDARY 定义输出图形需要的网格范围。

42

BOX

COPY

DEPTH

DIFFR

DIFFTHT

DIFFX

DIFFY

DIFFZ

EDITNNC

ENDBOX

EQUALS

MULTIPLY

PORV

TRANR

TRANTHT

TRANX

TRANY

TRANZ

重新定义目前输入盒子。

在目前盒子中,把资料从一个指定数组拷贝到另

一个数组。

目前输入盒子内的网格中心深度。

目前输入盒子中X方向上的扩散率(仅用于径向

分子扩散选择程序)。

目前输入盒子中,网格θ方向上的扩散率(仅用

于径向分子扩散选择程序)。

目前输入盒子中,网格X方向上的扩散率(仅用

于笛卡儿图形网格分子扩散选择程序)。

目前输入盒子中,网格Y方向上的扩散率(仅用

于笛卡儿图形网格的分子扩散选择程序)。

目前输入盒子中,网格Z方向上的扩散率(用于

径向或笛卡儿分子扩散选择)。

修饰不相邻连结。

重新定义整体网格的目前输入盒子。

在目前盒子中,给指定常数设置指定数组。

在目前盒子中,用指定常数乘以指定数组。

目前输入盒子中,网格的孔隙体积。

目前输入盒子中,网格R方向上的传导率(仅用

于径向图形)。

目前输入盒子中,网格θ方向上的传导率(仅用

于径向图形)。

目前输入盒子中,网格X方向上的传导率(仅用

于笛卡儿图形)。

目前输入盒子中,网格Y方向上的传导率(仅用

于笛卡儿图形)。

目前输入盒子中,网格Z方向上的传导率(用于

径向或笛卡儿图形)。

2.4 PROPS 部分

综 述

PROPS输入的数据包括与油藏流体和岩石性质有关的油层压力和饱和度。数据只能用一个能被接受的多表格式的关键字输入。每个类型表格的数量在RUNSPEC中说明。必须提供表格的准确数量。RUNSPEC中还规定了

43

每个表格的最大尺寸,当在一个关键字后面输入多表格时,每个表格用斜线号(/)终止。如:

SWFN

data for table 1 (表1数据)

/

data for table 2 (表2数据)

/

data for table 3 (表3数据)

/

需要的关键字由RUNSPEC记录6相位选择和记录6、10、12的特殊选择确定。

数据必须包括:流体的地面密度、计算的油藏流体密度和粘度、岩石压缩系数以及与饱和度有关的相对渗透率和毛管压力。

可动油的PVT性质可以用两种格式中的任何一种输入,(1)PVCO和PMAX;或(2)PVTO。第一种格式假定一个特定Rs值的 饱和油有一个压力无关的压缩系数,还假定 饱和油的粘度有一个与压力无关的变化率。第二种格式要复杂些,但在定义 饱和性质中允许有更大的灵活性。

在2组关键字中选1组输入与饱和度有关的相对渗透率和毛管压力。第一组:SWOF和(SGOF或SLGOF),允许用户在与气、水相对渗透率相同的表中输入油的相对渗透率。第二组:SWFN、SGFN和(SOF3或SOF2)需要在一个与油饱和度有关的独立表中输入油的相对渗透率。除了第二组必须用在气水二相以及混相驱或混相溶解选择的情况下,用户可根据原始数据选择一组最合适的关键字。两组关键字不能混用。设置油藏特性参数所需的最少的关键字是:

GRAVITY or DENSITY 用于地面油比重/密度。

ROCK 岩石压缩系数(除非使用岩石压缩

选择程序时设置)。

PVTO or (PVCO and PMAX) 如果是可动油。

PVTDO 如果是死油。

RSCONST or RSCONSTT 如果死油是用来模拟溶解气浓度恒

定时。

PVTG 如果是湿气。

PVTDG 如果是干气。

PVCONST or RVCONSTT 如果干气是用来模拟挥发油浓度恒

定时。

PVTW 如果是水(用PVTWSALT代替海水选

择)(并且可选用PVTWSALT替代海

44

水)。

第一组饱和度功能关键字:

SWOF 如果是油、水两相。

SGOF or SLGOF 如果是气、油两相。

第二组饱和度功能关键字:

SOF3 用于3相中的油相。

SOF2 用于2相中的油相。

SGFN 如果是气。

SWFN 如果是水。

有些关键字是相互独立的,不能混用。禁止混用的例子有:

PVTO and RSCONST (T) RS值恒定设置,仅用于死油。

PVTG and RVCONST (T) RV值恒定设置,仅用于干气。

SWOF and SGFN 两组饱和度关键字不能混用。

本部分的关键字可以用任何命令说明。所有关键字必须从一行的第1列开始,所有字符从1到第8列有效。这行余下部分(从第9列起)由编辑程序使用。

下表为不同相态下关键字例子综合表

相 态

+

+

+

+

+

关 键 字

PVTG

PVTG

PVDG

PVDG

SGFN

PVTG

PVTG

SGOF或SLGOF

SGFN

RSCONST(T) PVTW

SWOF

SWFN

PVTW

SWFN

PVTW

PVTW

SWEN

PVTW

PVTW

湿气

挥发油

干气

+

水 PVCO PMAX

或PVTO

或SOF3

PVCO PMAX

或PVTO

或SOF3

PVCO PMAX

或PVTO

或SOF2

PVDO

溶解气

溶解气 SGOF或SLGOF SWOF

湿气

挥发油 溶解气

饱和油 水

定量溶解气 或SOF2

干气

干气

+

PVDG

PVDG

SGFN

水 RVCONST(T)

定量挥发油

本部分可能出现的所有关键字与其功能的简要描述列表如下。最后的是

45

按字母顺序列表。更详细的描述可在关键字手册中查到。

关键字

输出

RPTPROPS

状态

简要描述

(任意选择) 设置油藏性质参数输出格式。

压力相关资料

DENSITY

GRAVITY

(需要的密度或比重) 储罐流体密度。

(需要的密度或比重) 储罐流体比重。

PMAX

PVCO

PVTO

PVDO

RSCONST

RSCONSTT

PVTG

PVDG

RVCONST

RVCONSTT

PVTW

PVTWSALT

ROCK

ROCKTAB

ROCKTABH

(使用PVCO时要求)(可动油系统要求用PVCO或PVTO)

(可动油系统要求用PVCO或PVTO)

(死油或 饱和油系

统要求)

(仅适用于死油的任意选择程序)

(仅适用于死油的任意选择程序)

(湿气系统要求)

(干气或 饱和气系

统要求)

(仅用于干气的任

意选择程序)

(仅用于干气的任

意选择程序)

(水系统要求,但

海水不能选择)

(在RUNSPEC中,

若QBRINE=T时要

求)

(若QROCKC=F时

要求)

(若QROCKC=T时

要求)

(如果QROCKH=T

时要求)

46

在模拟中预测的最大压力。

随泡点压力变化的可动油的体

积系数和粘度。

随压力和RS的可动油的体积系数和粘度。

随压力变化的死油的体积系数

和粘度。

全油田死油的RS值为定值。

对不同PVT区的死油设置RS值为恒定。

随压力和RV变化的气体的体积系数和粘度。

随压力变化的干气的体积系数

和粘度。

全气田干气的RV为定值。

对不同PVT区的干气,其RV为定值。

地层水的体积系数、压缩系数

和粘度。

随含盐量变化的水的体积系数

、压缩系数和粘度。

岩石的压缩系数。

岩石压实表。

岩石滞后压实表。

OVERBURD (选择用QROCKC 上覆岩层压力表。

或QROCKH设置T)

饱和度相关资料

SGOF (在气、油二相系统中,要用第

1组的关键字SGOF或SLGOF)

SLGOF (在气、油二相系统中,要用第

1组关键字SGOF或SLGOF)

SWOF (在油、水二相系统中,要用第

1组关键字)

SGFN (在2相或有气存在的3相系统

中,要用第2组关键字)

SOF2 (有油存在的2相系统中,要

用第2组关键字)

SOF3 (3相系统中,要用第2组关

键字)

SWFN (在2相或有水存在的3相系

统中,要用第2组关键字)

STONE1 (任意选择)

STONE2 (任意选择)

STOG (在RUNSPEC中,如果QSURFT

=T,任意选择)

STOW (在RUNSPEC中,如果QSURFT

47

在束缚水条件下,随

气体饱和度变化的气

的相对渗透率,毛管

压力和油对气的相对

渗透率。

在束缚水条件下,随

液体饱和度变化的气

的相对渗透率,毛管

压力和油对气的相对

渗透率。

随水饱和度变化的水

的相对渗透率,毛管

压力和油对水的相对

渗透率。

随气饱和度变化的气

相对渗透率和毛管压

力。

随油饱和度变化的油

的相对渗透率。

随油饱和度变化的油

的相对渗透率。

随水饱和度变化的水

的相对渗透率和毛管

压力。

用Stone修改的公式1

计算3相油的相对渗

透率。

用Stone修改的公式2

计算3相油的相对渗

透率。

油气表面张力表。

油水表面张力表。

=T,任意选择)

端点标定

ENKRVD (仅选用于饱和度表端点标定) 相对渗透率的端点值随深

度而变。

ENPTVD (仅选用于饱和度端点标定) 饱和度的端点值随深而

变。

KRG (仅选用于饱和度端点标定)

KRO (任意选择,仅用于饱和度表

终点刻度选择)

KRW (仅选用于饱和度表端点标定)

SCALELIM (仅选用于饱和度表端点标定)

SGCR (仅选用于饱和度表端点标定)

SGL (仅选用于饱和度表端点标定)

SGU (仅选用于饱和度表端点标定)

SOGCR (仅选用于饱和度表端点标定)

SOWCR (仅选用于饱和度表端点标定)

SWCR (仅选用于饱和度表端点标定)

SWL (仅选用于饱和度表端点标定)

SWU (仅选用于饱和度表端点标定)

TZONE (选用于可动流体校正情况)

48

在目前输入盒子中,网格

节点气的相对渗透率的端

点值。

在目前输入盒子中,网格

节点油的相对渗透率的端

点值。

在目前输入盒子中,网格

节点水的相对渗透率端点

值。

对每个端点的刻度和深度

范围设置饱和度标定界限

值。

目前输入盒子中,网格临

界气饱和度。

目前输入盒子中,网格原

始气饱和度。

目前输入盒子中,网格最

大气饱和度。

目前输入盒子中,网格残

留气中的油饱和度。

目前输入盒子中,网格残

留水中的油饱和度。

目前输入盒子中,网格临

界水饱和度。

目前输入盒子中,网格的

束缚水饱和度。

目前输入盒子中,网格最

大含水饱和度。

控制临界端点与深度的变

化。


本文标签: 网格 程序 压力