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

更新时间:2024-07-01 21:38:01 阅读量: 综合文库 文档下载

说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。

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

COMPFLSH 井射孔段的闪蒸转化比???????????????65 COMPIMB 井射开网格的渗吸表号???????????????67 COMPINJK 用户定义的注入井相对渗透率????????????68 COMPLUMP 为自动修井而将射开网格归在一起??????????69 COMPRP 重新标定井射开节点的饱和度数据????????? 70 COMPVE 垂直平衡(V.E.)运行时,井射孔深度的重设定 ????72 COORD 坐标线????????????????????? 75

COORDSYS 坐标系统信息???????????????????76

COPY 从一个数组拷贝数据到另一数组 ?????????? 77 COPYBOX 从一个BOX向另外一个拷贝一组网格数据 ?????? 77 CRITPERM 对VE节点压缩的渗透率标准 ????????????78 DATE 输出日期到汇总文件????????????????79 DATES 模拟者事先指定报告日期??????????????79 DATUM 基准面深度,用于深度校正压力的输出 ??????? 80 BEBUG 控制检测输出 ??????????????????80 DENSITY 地面条件下流体密度??????????????? 81 DEPTH DIFFC DIFFDP DIFFMMF DIFFMR DIFFMTHT DIFFMX DIFFMY DIFFMZ ZDIFFR DIFFTHT DIFFX DIFFY YDIFFZ DPGRID DR RDRSDT DRV RDRVDT DTHETA DTHETAV DX XDXV XDY YDYV YDZ ZDZMTRX DZMTRXV DZNET ECHO

EDITNNC EHYSTR END ENDBOX ENDNUM 网块中心深度?????????????????? 82 每一个PVT区域的分子扩散数据????????? 82 在双重介质运行中,限制分子扩散????????? 83 基质一裂缝的扩散乘子?????????????? 83 方向的扩散乘子?????????????????83 θ方向扩散系数乘子????????????????84 方向的扩散乘子?????????????????84 方向的扩散乘子?????????????????85 方向的扩散乘子 ???????????????? 85 方向的扩散系数?????????????????86 θ方向的扩散系数???????????????? 86 方向扩散系数??????????????????87 方向扩散系数??????????????????87 方向扩散系数 ????????????????? 88 对裂缝单元使用基质单元的网格数据???????? 88 方向网格的大小?????????????????88 溶解GOR的增加的最大速度?????????????89 方向网格大小(矢量) ?????????????? 89 挥发油的OGR的增加的最大速度??????????? 90 θ方向的网格大小????????????????? 90 网格的角度大小(向量) ?????????????? 91 方向的网格大小???????????????? 91 方向网格大小(向量)???????????????91 方向网格大小??????????????????92 方向网格大小(向量)???????????????92 方向网格大小??????????????????92 基质块的垂直尺寸?????????????????93 基质岩体块的垂直尺寸(向量)??????????? 93 净厚度????????????????????? 93

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

改变非相邻连接??????????????????94 滞后作用参数和模型选择??????? ?????? 95 标志SCHEDULE部分的结束 ????????????95 将BOX恢复到包含全部网格?????????????95 端点标定与深度区域号 ??????????????95

R X Y R X Z ENKRVD 相对渗透率端点与深度关系表 ??????????? 96 ENPTVD 饱和度端点与深度关系表??????????????97 EQLNUM 平衡区号数 ?????????????????? 98 EQUALS 在目前的BOX中设置数组为常数?????????? 99 EQUIL 平衡数据详述????????????????? 99 EXTRAPMS 对表的外插请求预告信息 ???????????? 101 FIPNUM 流体储量区域号????????????????? 102 GCONINJE 对井组井/油田注入率的控制/限制 ?????????102 GCONPRI 为“优先”而设的井组或油田产量限制?????????104 GCONPROD GCONSALE

GCONSUMP GCONTOL GECON GLIFTLIM GRAVITY GRIDFILE GRUPRIG GRUPTREE GSEPCOND IMBNUM IMBNUMMF IMPES IMPLICIT INCLUDE INIT INRAD KRG KRNUM KRNUMMF KRO KRW LOAD MESSAGES MINPV MINPVV MISCNUM MONITOR MULTIPLY MULTR RMULTTHT THETAMULTX XMULTY MULTZ Z 井组或油田的产率控制或限制??????????? 104 井组或油田的售气控制产率???????????? 107

井组的气消耗率和引进率 ???????????? 109 井组控制目标(产率)允许差额 ???????????110 井组或油田的经济极限数据???????????? 111 最大井组人工举升能力?????????????? 112 地面条件下的流体密度 ??????????????113 控制几何文件网格的容量 ???????????? 113 给井组配置修井设备??????????????? 113 建立多级井组控制的树状结构 ???????????114 井组设置分离器????????????????? 115 渗吸饱和度函数据区域号????????????? 115 基质—裂缝渗吸区域号?????????????? 116 建立IMPES求解过程 ?????????????? 117 重建全隐式求解 ?? ?????????????? 117 包含数据文件名 ?????????????????117 要求输出初始文件 ??????????????118 径向模型的内径 ???????????????? 118 标定气相对渗透率的端点 ????????????? 118 方向性相对渗透率表格数 ???????????? 119 基岩—裂缝流动饱和度表号 ???????????? 120 标定油相对渗透率端点 ?????????????? 120 标定水相对渗透率端点 ?????????????? 121 调入一个SAVE文件以便执行一个快速重起动????? 122 重设置打印和停止限定的信息 ??????????? 123 设置活动网格的最小孔隙体积????????????124 建立一个有效网格的最小孔隙空间 ????????? 124 混合区数目 ?????????????????? 125 请求实时显示输出 ??????????????? 125 当前定义区中的数组??????????????? 126 方向传导率乘子??????????????? 126 方向传导率乘子???? ??????????127 方向传导率乘子???????????????? 127 方向传导率乘子???????????????? 127 方向传导率乘子 ????????????????128

YNEWTON 输出迭代计数到汇总文件 ?????????????128 NEWTRAN 标定使用块拐角传导率 ??????????????128 NEXTSTEP 建立下一时间步最大值 ??????????????129 NNC 非相邻连接的直接输入 ??????????????129 NOECHO 关闭输出的响应 ???????????????? 130 NOGGF 压缩网格几何模型文件 ??????????????130 NODPPM 非双孔的渗透率乘子 ???????????????130 NOWARN 压制ECLIPSE警报信息?????????????? 130 NTG 厚度净毛比??????????????????? 130 OILAPI OLDTRAN OLDTRANR

OPTIONS OUTRAD OVERBURD PERMR RPERMTHT PERMX PERMY PERMZ PINCH

PINCHOUT PMAX PMISC PORO PORV PRESSURE PRIORITY PRVD PSEUDOS PVCO PVDG PVDO PVTG

PVTNUM PVTO PVTW

PVTWSALT QDRILL RESTART RESVNUM ROCK ROCKNUM ROCKTAB

初始原油API值,以便API示踪选择?????????131 标定块中心传导率?????????????????131 标定任意一块中心传导率????????????? 131

开启特别程序选择???????????????? 132 径向模型外半径 ????????????????? 134 岩石负载压力表??? ??????????????135 方向绝对渗透率???????????????? 135 θ方向绝对渗透率???????????????? 136 X方向绝对渗透率????????????????136 Y方向绝对渗透率 ?? ???????????? 136 Z方向绝对渗透率????????????????137 建立尖灭层上下的连接 ????????????137

建立尖灭层上下的连接?????????????138 模拟中的最大压力???????????????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

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

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

不相邻连接

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

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

2.允许穿过平移断层在不同层之间进行流动模拟; 3.允许在井筒附近进行局部网格加密,它还包括了在一个笛卡尔座标系统内对锥形网格模拟的功能。

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

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

运行次数的标定

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

垂向平衡

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

的垂向平衡选择程序能将3D模型当成2D平面模型来运算。

双孔隙度

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

在基质网格和裂缝网格之间流体的重力吸吮/驱替和分子扩散等现象都能解释。 原油的高压物性(PVT)数据和岩石数据。

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

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

相对渗透率的方向性

相对渗透率曲线可以分别给定为水平方向流动和垂直方向流动的二个相对渗透率曲线。对每个饱和度函数区,用户可以有给定3条(X、Y、Z)或6条(+X,-X,+Y,-Y,+Z,-Z)相对渗透率曲线的选择,而不是常用一条曲线。这种新增的自由给定相对渗透率曲线的条数,油藏工程师能用比较简单的方法提高油藏物理性质的模拟研究。例如,大部分油藏节点在垂向上分布少,在平面上分布多,如果水平相对渗透率小于巴克里-莱弗里特饱和度,那么水平流动模拟将会能更真实的模拟油藏实际情况。在Eclipse拟函数选择程序中,同样也需要相对渗透率的方向性,它可以用于三维粗网格的模拟中。

饱和度数据表的刻度

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

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

滞后效应

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

岩石压实

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

示踪剂追踪

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

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

义。

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

原油比重(API)跟踪

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

海水盐度跟踪

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

混相驱动能力

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

单井控制程序

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

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

生产井同样还可受到另一种的经济约束。如果一口生产井的产油量或产气量低于经济界限,那么该井将自动关井或重新完井。如果生产井的含水率、气/油比或水/气比超过规定的上限,那么生产井就要修井或关井。修井速度取决于占机提供修井的能力。生产井的含水率、气/油比和水/气比的生产界限同样可以用于同生产井连通的其它井,若这些界限超过规定,也应关井。此外,每当生产井的超过界限时,我们可以逐渐地减少生产井的产量来限制生产井含水量,气/油比或水/气比的上升。生产井还受到最小井口压力的限制,一旦生产井的产油量或产液量低于规定的界限,可以用自动开关或人工举升方法来提高它们的生产能力。对于人工举升的生产井还受到井组或油田产能的限制。一些生产井由于某种理由关井要设置一个使模拟自动停止的标记,而其它井要有自动开井的标记。若油田上所有的生产井,或一个或更多的井组关井,那么这个运算同样可以停止。关闭的井可以进行定期测试,了解它们是否具有再次生产的经济效益。

井组、油田的生产控制

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

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

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

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

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

多级开采系统

用表明某些井组是高一级井组的子井组的办法,可将井组按任一所需的层次排列成一个多级树式结构。所有井组的控制和限制条件都可应用任何一级的井组(包括油田亦可看作成树状结构顶部的0级井组)。

注水井组的控制

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

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

(2)多相流体注入油藏的总体积量; (3)单相流体的回注比例;

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

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

销售气生产控制

销售气生产控制程序提供了一个从油田伴生产输出的控制方法,这个方法同时可对油的生产量进行控制。一个井组或油田的商品气的生产量可看作为总的产气量减去它总的注气

量,再减去井组和子井组气的消耗量,加上输入到井组或子井组的气量(若有的话)。对每个井组可单独设置气的消耗量和输入量。

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

井的层间交互和混流

当一口井完井于几个网格区内,那么每一个网格的流量与产生的下列三个参数成正比: (1)网格和井筒之间的传导率; (2)网格内射孔段相的流度; (3)网格和井筒之间的压力降。

当在不同射开网格间分配井的目标产量时,Eclipse软件充分考虑了上述的三个参数。这对于在有二个或二个以上的,垂向连通性比较差的油层,每个油层的压力降又是各不相同的。完井特别重要。

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

大斜度井和水平井

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

专门的气井程序

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

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

细网格平衡程序

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

通常,使每一个网格区内的流体饱和度等于网格区内中心的那一个点的饱和度。然而,当一个大的网格区含有流体接触面时或具有流体的过渡带部分时,这种方法在流体储量的估算能产生明显的误差。在这种情况下,Eclipse软件具有提高地层流体储量计算精度的

功能。在网格区内,可以确定若干面上流体的饱和度,将这些值平均就可以得到网格区内总的流体饱和度。这种方法还考虑了倾斜网格的形状,用网格区内每一个面上的流体饱和度与网格区内每一个水平横切面的面积进行权衡得到倾斜网格的流体平均饱和度。

水层模拟

水层模型可用解析模型和模拟模型任选一个来代表。解析模型是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。

不同的多相流动对比可用于不同部分的油管计算,即不同的对比能使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程序使用了在细网格模拟期间产生的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相组分的联立方程,以确保在所有的时间步长中聚合物完全稳定溶解。

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)来设置月生产的目标产量;

(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用的多级分类管理系统,但另外还带有生产井向分离器和天然气工厂供气的复杂部分。根据流体的地面体积,油藏体积,井底压力等可以强有力的对单井、井组、油田各级的生产量进行控制。对特定的井或井组的一部分流体实行回注是有可能的。在组分模型里,注入量可以专门规定,或减少其它井或井组的注入量。

GRID—网格程序

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

ECLMAP—Eclipse绘图软件包

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

PVT—状态方程软件包

PVT是一个人机对话状态方程(EOS)软件包,为了确定油藏流体的相态,该软件包用

于实验室的资料分析。用物理平衡检验,可测试实验室数据测量的质量。利用立体状态方程的变化来模拟实验室的分析数据,并且用回归一个或多个EOS的参数使测量值和计算值之间的误差减到最小。于是EOS模型常可用于产生适用于Eclipse100,200或300所需的数据。

WELLTEST—试井程序

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

CLAN—测井分析软件包

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

2

综 述

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

数据文件综述

关 键 字

内 容(简述)

RUNSPEC(必须) GRID(必须) EDIT(选择) PROPS(必须) REGIONS(选择)

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

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

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

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

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

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

对油藏原始条件的确定:a)用规定的流体接触面深度到可能的流体高度(平衡)计算;b)从上一次运算建立的重启文件中读出;或c)用户自己确定每一个网格区的原始条件(一般不用此项)。 在每一个时间步长后,确定的数据被编入汇总文件中。如果必要的话在每次运算结束后输出某种类型成果的图形(如,含水量与时间关系曲线)。若这部分省略,则不产生汇总文件。

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

SOLUTION(必须)

SUMMARY(选择)

SCHEDULE(必须)

和油管模拟参数。

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

注意每部分文头的关键字(像所有其他的文头关键字一样)必须从第1列开始。所有的

字符在1—8格内有效。该行上的其余字符(从第9列开始)当作为注解。

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

3*450 等于 450

450

450

400

350

350

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

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

*

*

400 2*350

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

关 键 字

COLUMS DEBUG ECHO EXTRAPMS INCLUDE MESSAGES NOECHO NOWARN OPTIONS

简 述

为了能读出数据输入文件,对左右二边重新定位键。 输出控制键调试设置。

每一个运算开始时数据重新打印键。 要求外推PVT和VEP数据表信息预告。 插入规定文件的目录。 所有正规信息打印和停止限止重放键。 每一个运算开始的不进行数据重新打印。 Eclipse信息预告抑止键。

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

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

RUNSPEC ? GRID RPTGRID 6*0

11*

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

数据文件实例:

RUNSPEC . . . GRID

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

RPTGRID 6*0 11*1 /

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

‘GRID.SECT’ / PROPS RPTPROPS 20*1 INCLUDE

‘PROPS.SECT’ / REGIONS RPTREGS 200 INCLUDE

‘REGIONS.SECT’ SOLUTION RPTSOL 20*1 INCLUDE

‘SOLUTION.SECT’ SUMMARY INCLUDE

‘SYSTEM_DEFAULT:SUMMARY.SECT’ SCHEDULE RPTSCHED 0

1

0

1

1

0

3*2

/

/

/

/

/

*

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

/

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

/

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

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

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

INCLUDE

‘SCHEDULE.SECT’ END

/

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

2.1 RUNSPEC部分

综 述

SUNSPEC部分是Eclipse数据输入文件的第一个部分。文件中包括了运算题目,开始日期,单位,各种作业的维数(网格区数、井数、表数等)目前相标志和选择键。只有注解和关键字INCLUDE ECHO,NOECHO,COLUMNS OPTIONS和DEBUG放在运行定义卡的前面。除非从装有Runspec数据文件的存盘中重新运算,运行定义卡这部分始终是存在的。

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

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

Record 2 Record 4 Record 6 Record 8 Record 18 Record 28 Record 30

题目

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

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

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

模拟的开始日期。

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

还可以缺省的其他项目是:笛卡儿座标图形,不相邻连接的分散流动模型,一套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

层将当作为骨架网格,其余的则作为裂缝网格。裂缝网格渗透率和裂缝

网格孔隙度相乘,骨架和裂缝总成对存在,并且用SIGMA关键字来决定

裂缝网格的渗透率。若QDPORO是T,则QNNCON同样应为T。

〈F〉 QDPERM—在双孔隙度运算中,T为骨架与骨架之间的传导率,不等于0(双孔隙

度/双渗透率)。F为骨架与骨架之间的传导率,等于0。若QDPORO为F,则QDPERM则没有作用。 Record 6

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

〈F〉 OIL 〈F〉 WATER 〈F〉 GAS

T或F

T或F T或F T或F

〈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控制聚合物—盐的敏感性分析选择

(见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(用油藏特性部分中的关键词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键

必须建立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—与深度表对应的任何饱和度表的端点的最多节点数(见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。

〈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被 定义之前,一个BOX一直存在,或由关键字ENDBOX来读得。缺省BOX是整体模型。面积的大小、X、Y、R、θ值可以按照矢量输入,例如DXV(I),这里仅用I变化DX。 对于常数,可以用关键字EQUALS设置,或用关键字ADD、MULTIPLY、COPY处理。两者任选一种。

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

基 本 关 键 字 表

块中心图形网格 笛卡儿

DXV或DX DYV或DY DZ TOPS PORO PERMX PERMY

径向 INRAD

角点图形网格 COORD ZCORN

PORO

DRV或DR或OUTRAD

DTHETAV或DTHETA DZ TOPS PORO

PERMR PERMTHT PERMX/R

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

DZ

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

目前输入盒子网格的

Way是单数,若多数则:Seller- Al ways几种方法!

(必须用某种方法设置)

顶面深度。

目前输入盒子Z方向 (垂直方向)的网格尺寸。

笛卡儿块中心图形网格

DX DXV

(必须用某种方法设置)

目前输入盒子X方向的网

格尺寸。

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

DY 或 DYV INRAD DR

或 DRV

(必须用某种方法设置)

目前输入盒子Y方向的网 格尺寸。

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

(DY的简化版本)

径向块中心图形网格

(径向块中心图形网格要求) 设置径向网格的内半径。 (必须用某种方法设置)

(DR的简化版本)

目前输入盒子R方向的网

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

设置径向网格外半径。 目前输入盒子θ方向的网 格尺寸。

OUTRAD (径向网格任意选择)

DTHETA (必须用某种方法设置)

DTHETAV (Dθ的简化版本)

θ方向网格尺寸的矢量。

角点图形网格

COORD ZCORN

(角点图形网格要求) (角点图形网格要求)

(任意选择)

定义网格范围内每个油层

和包括所有网格角点(I、J)的范围。 网格角点的深度。 有关网格范围内每个油层 坐标系统的资料。指定了 角点图形网格和块中心图 形网格以循环方式完成。 目前输入盒子内网格的孔隙度。

网格的有效厚度。

COORDSYS PORO DZNET 或 NTG

所有图形网格的参数

(必须用某种方法设置) (有效厚度或净毛比任意选择)

目前输入盒内网格的净毛比。 识别活动网格。

每个活动网格的最小孔隙体积。

ACTNUM (任意选择) MINPV (任意选择)

MINPVV (任意选择) 目前盒内网格的最小孔隙体积。 PERMR (如果径向NR>1,必须 目前输入盒内R方向上渗透率值。 用某种方法设置)

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

用某种方法设置)

(如果笛卡儿座标NX>1, 目前输入盒内X方向上渗透率值。 必须用某种方法设置)

(如果笛卡儿座标NY>1, 目前输入盒内Y方向上渗透率值。 必须用某种方法设置)

(任何系统NZ>1,必须 目前输入盒内z方向上渗透率值。 用某种方法设置)

PERMX PERMY PERMZ

传导率关键字

NNC NEWTRAN OLDTRAN OLDTRANR

(任意选择) (任意选择) (任意选择) (任意选择)

不相邻连接显式通道。

说明计算传导率的角点网格数据。 说明计算传导率的块中心网格数据。 选择网格特征计算传导率。

传导率修饰

MULTR

MULTX MULTY MULTZ

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

目前输入盒内R方向上传

导率乘数。 目前输入盒内θ方向上传 导率乘数。

目前输入盒内X方向上传 导率乘数。

目前输入盒内Y方向上传 导率乘数。

目前输入盒内Z方向上传 导率乘数。

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

(如果笛卡儿座标NX>1,任 意选择)

(如果笛卡儿座标NY>1,任 意选择) (如果NZ>1,任意选择)

确定含水层范围

AQUCON AQUNUM

(定义含水层要求) (定义含水层要求)

定义含水层与油层连通。 定义含水层范围。

双孔隙度

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

用户仅能输入骨架网格数据。 全部网格中,典型骨架网格的 目前盒内典型骨架网格垂向长

垂向长度。

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

度。

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

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

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

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

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

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

散率乘数。 散率乘数。

DIFFMY (如果笛卡儿NY>1,任意选择) 目前输入盒内Y方向上扩 散率乘数。 DIFFMZ (如果NZ>1,任意选择) 目前输入盒内Z方向上扩散率

乘数。

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

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

操作关键字

ADD

(任意选择)

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

数。

定义输出图形需要的网格范围。 重新定义目前输入盒子。

在目前盒子中,把资料从一个数组拷贝到另 重新定义含有整体网格的目前输入盒子。

BOUNDARY (任意选择) BOX (任意选择) COPY

(任意选择)

一个数组。

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

产生GRAF和PSEUDO原始文件。 控制网格图形文件目录。(为了图形输出)

ENDBOX (任意选择) EQUALS (任意选择)

INIT (任意选择) GRIDFILE (任意选择) MULTIPLY (任意选择) NOGGF (任意选择) PSEUDO (任意选择) RPTGRID (任意选择) SPECGRID

关 键 字 ACTNUM ADD AQUCON AQUNUM BOUNDARY BOX COORD

在目前盒内,用指定常数乘以指定数组。 说明产生没有网格图形文件(为了图形输出) 。

打开拟函数软件包数据的输出键。 设置网格数据的输出格式。

说明网格尺寸。

(任意选择)

关 键 字 按 字 母 顺 序 列 表 简 单 描 述

识别活动网格。

在目前盒子内,给指定数组加一个指定常数。 定义含水层和油层间的连通。 定义含水层范围。 定义输出图形需要的网格范围。

重新定义目前输入盒子。

定义网格范围内,每个油层和包括所有网格角点 (I、J)的范围。 指定角点和块中心两种图形网格按循环方式完成 。如果网格内有多个油层,对网格进行分层。

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

COORDSYS COPY

DIFFMMF DIFFMR DIFFMTHT DIFFMX DIFFMY DIFFMZ DPGRID DR

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

目前输入盒内θ方向上扩散率乘数。 目前输入盒内X方向上扩散率乘数。 目前输入盒内Y方向上扩散率乘数。 目前输入盒内Z方向上扩散率乘数。 用户仅能输入骨架网格数据。

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

DRV DIHETA DTHETAV DX DXV DY DYV

DZ 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

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

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

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

整体网格中,典型基质网格的垂向长度。 目前盒子内,典型骨架网格的垂向长度。 目前盒子内,网格的有效厚度。 重新定义含有全部网格的目前盒子。

在目前盒子内,给指定常数设置指定数组。 控制网格图形的文件目录(为了图形输出)。 产生GRAF和PSEUDO原始文件。 设置径向网格内半径。

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

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

在目前盒子中,用指定常数乘以指定数组。 目前盒子内,R方向上传导率的乘数。 目前盒子内,θ方向上传导率的乘数。 目前盒子内,X方向上传导率的乘数。 目前盒子内,Y方向上传导率的乘数。 目前盒子内,Z方向上传导率的乘数。

说明用网格间接触面计算的传导率。角点图形网格缺省。不用于块中心图形网格。

不相邻连结显式通道。

停止产生网格图形文件(为了图形输出)。 规定裂缝网格的渗透率不能乘裂缝孔隙度。 目前盒子内网格的净毛比。

规定块中心图形网格数据常用于计算传导率。 选用块中心资料计算传导率。 设置径向网格的外半径。

目前盒子内,R方向上的渗透率。 目前盒子内,θ方向上的渗透率。 目前盒子内,X方向上的渗透率。 目前盒子内,Y方向上的渗透率。 目前盒子内,Z方向上的渗透率。 目前盒子内,网格的孔隙度。 打开拟函数软件包数据的输出键。

设置网格数据的输出格式。 为了容易重新启动需要产生一个SAVE文件。 整体网格中,骨架~裂缝的连结系数。

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

SPECGRID 网格大小和字符的规定。 TOPS 目前盒子内,网格的顶面深度。 ZCORN 网格的角点深度。 数据设置实例:

2.3 EDIT 部分

综 述

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

网格输出数组关键字: PORV

DEPTH

TRANX

TRANY

TRANR DIFFX TRANTHT DIFFY

DIFFR DIFFTHT

TRANZ DIFFZ

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

BOX EQUALS ADD ENDBOX COPY MULTIPLY

修饰不相邻连结

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

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

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

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

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

ADD BOUNDARY BOX COPY DEPTH DIFFR

在目前输入盒子中,给指定数组加指定常数。 定义输出图形需要的网格范围。 重新定义目前输入盒子。

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

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

目前输入盒子中X方向上的扩散率(仅用于径向 分子扩散选择程序)。

目前输入盒子中,网格θ方向上的扩散率(仅用 于径向分子扩散选择程序)。

DIFFTHT

DIFFX DIFFY DIFFZ

目前输入盒子中,网格X方向上的扩散率(仅用 于笛卡儿图形网格分子扩散选择程序)。 目前输入盒子中,网格Y方向上的扩散率(仅用 于笛卡儿图形网格的分子扩散选择程序)。 目前输入盒子中,网格Z方向上的扩散率(用于 径向或笛卡儿分子扩散选择)。 修饰不相邻连结。

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

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

在目前盒子中,用指定常数乘以指定数组。 目前输入盒子中,网格的孔隙体积。

EDITNNC ENDBOX EQUALS MULTIPLY PORV TRANR

目前输入盒子中,网格R方向上的传导率(仅用 于径向图形)。

目前输入盒子中,网格θ方向上的传导率(仅用 于径向图形)。

目前输入盒子中,网格X方向上的传导率(仅用 于笛卡儿图形)。

目前输入盒子中,网格Y方向上的传导率(仅用 于笛卡儿图形)。 目前输入盒子中,网格Z方向上的传导率(用于 径向或笛卡儿图形)。

TRANTHT TRANX TRANY TRANZ

2.4 PROPS 部分

综 述

PROPS输入的数据包括与油藏流体和岩石性质有关的油层压力和饱和度。数据只能用一个能被接受的多表格式的关键字输入。每个类型表格的数量在RUNSPEC中说明。必须提供表格的准确数量。RUNSPEC中还规定了每个表格的最大尺寸,当在一个关键字后面输入多表格时,每个表格用斜线号(/)终止。如: SWFN data /

data /

data /

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

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

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

for for

table

2

1

(表1数据)

table

(表2数据) 3

(表3数据)

for table

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

GRAVITY or DENSITY ROCK PVTDO PVTG

用于地面油比重/密度。

岩石压缩系数(除非使用岩石压缩 选择程序时设置)。

如果是可动油。 如果是死油。

PVTO or (PVCO and PMAX) RSCONST or RSCONSTT

如果死油是用来模拟溶解气浓度恒 定时。

如果是湿气。

如果是干气。

如果干气是用来模拟挥发油浓度恒 定时。

如果是水(用PVTWSALT代替海水选 择)(并且可选用PVTWSALT替代海 水)。

PVTDG

PVCONST or RVCONSTT PVTW

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

SWOF 如果是油、水两相。 SGOF or SLGOF 如果是气、油两相。 第二组饱和度功能关键字: SOF3 用于3相中的油相。 SOF2 SGFN

用于2相中的油相。 如果是气。

SWFN 如果是水。

有些关键字是相互独立的,不能混用。禁止混用的例子有: PVTO and RSCONST (T) RS值恒定设置,仅用于死油。 PVTG and RVCONST (T) SWOF and SGFN

RV值恒定设置,仅用于干气。

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

本部分的关键字可以用任何命令说明。所有关键字必须从一行的第1列开始,所有字符

从1到第8列有效。这行余下部分(从第9列起)由编辑程序使用。 下表为不同相态下关键字例子综合表

相 态 油 + 油 +

关 键 字

PVTG PVTG PVDG PVDG

PVTW PVTW

PVTW PVTW

湿气 +

挥发油 干气

PVCO PMAX 或PVTO 或SOF3 PVCO PMAX 或PVTO

溶解气

溶解气 油

SGOF或SLGOF SWOF SGFN PVTG

SWEN

或SOF3 PVCO PMAX

湿气

+ +

或PVTO

PVTG

SGOF或SLGOF SGFN

RSCONST(T) PVTW

SWOF SWFN

溶解气 +

饱和油

挥发油

或SOF2 PVDO

定量溶解气

或SOF2

干气 干气 +

PVDG PVDG SGFN

PVTW SWFN

水 RVCONST(T)

定量挥发油

本部分可能出现的所有关键字与其功能的简要描述列表如下。最后的是按字母顺序列表。更详细的描述可在关键字手册中查到。

关键字 输出

RPTPROPS

状态

简要描述

设置油藏性质参数输出格式。

(任意选择)

压力相关资料

DENSITY GRAVITY PMAX PVCO PVTO PVDO

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

(使用PVCO时要求) 在模拟中预测的最大压力。 (可动油系统要求用 随泡点压力变化的可动油的体 PVCO或PVTO) 积系数和粘度。 (可动油系统要求用 PVCO或PVTO) (死油或 饱和油系 统要求)

随压力和RS的可动油的体积系

数和粘度。

随压力变化的死油的体积系数 和粘度。

RSCONST

RSCONSTT PVTG PVDG

(仅适用于死油的任 全油田死油的RS值为定值。 意选择程序)

(仅适用于死油的任 对不同PVT区的死油设置RS值 意选择程序) 为恒定。 (湿气系统要求) 随压力和RV变化的气体的体积 (干气或 饱和气系 统要求) (仅用于干气的任 意选择程序)

(仅用于干气的任 意选择程序)

(水系统要求,但 海水不能选择) (在RUNSPEC中, 若QBRINE=T时要 求)

(若QROCKC=F时 要求)

系数和粘度。

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

和粘度。

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

RVCONST RVCONSTT

PVTW

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

和粘度。

随含盐量变化的水的体积系数 、压缩系数和粘度。 岩石的压缩系数。

PVTWSALT ROCK

ROCKTAB ROCKTABH OVERBURD

(若QROCKC=T时 要求)

(如果QROCKH=T 时要求)

(选择用QROCKC 或QROCKH设置T)

岩石压实表。 岩石滞后压实表。 上覆岩层压力表。

饱和度相关资料

SGOF

SLGOF

SWOF

SGFN SOF2 SOF3 SWFN

STONE1

STONE2

STOG STOW

端点标定 (在气、油二相系统中,要用第

在束缚水条件下,随

1组的关键字SGOF或SLGOF)

气体饱和度变化的气 的相对渗透率,毛管

压力和油对气的相对 渗透率。 (在气、油二相系统中,要用第 在束缚水条件下,随 1组关键字SGOF或SLGOF) 液体饱和度变化的气

的相对渗透率,毛管

压力和油对气的相对 渗透率。

(在油、水二相系统中,要用第 随水饱和度变化的水

1组关键字)

的相对渗透率,毛管

压力和油对水的相对

渗透率。 (在2相或有气存在的3相系统 随气饱和度变化的气 中,要用第2组关键字) 相对渗透率和毛管压 力。 (有油存在的2相系统中,要 随油饱和度变化的油

用第2组关键字)

的相对渗透率。

(3相系统中,要用第2组关 随油饱和度变化的油 键字) 的相对渗透率。 (在2相或有水存在的3相系 随水饱和度变化的水 统中,要用第2组关键字) 的相对渗透率和毛管 压力。

(任意选择) 用Stone修改的公式1 计算3相油的相对渗 透率。 (任意选择) 用Stone修改的公式2

计算3相油的相对渗

透率。

(在RUNSPEC中,如果QSURFT 油气表面张力表。

=T,任意选择) (在RUNSPEC中,如果QSURFT 油水表面张力表。

=T,任意选择)

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

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

KRO (任意选择,仅用于饱和度表 终点刻度选择) KRW (仅选用于饱和度表端点标定)

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

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

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

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

SWCR (仅选用于饱和度表端点标定) SWL (仅选用于饱和度表端点标定) SWU (仅选用于饱和度表端点标定) TZONE (选用于可动流体校正情况)

混相选择

SDENSITY (仅用于混相驱选择程序)

SORWMIS (仅选用于混相驱选择程序) TLMIXPAR (仅用于混相驱选择程序)

分子扩散选择

DIFFC

(仅用于分子扩散选择程序)在目前输入盒子中,网格 节点气的相对渗透率的端 点值。

在目前输入盒子中,网格 节点油的相对渗透率的端 点值。

在目前输入盒子中,网格 节点水的相对渗透率端点 值。 对每个端点的刻度和深度 范围设置饱和度标定界限 值。

目前输入盒子中,网格临 界气饱和度。 目前输入盒子中,网格原 始气饱和度。 目前输入盒子中,网格最 大气饱和度。 目前输入盒子中,网格残 留气中的油饱和度。 目前输入盒子中,网格残 留水中的油饱和度。 目前输入盒子中,网格临 界水饱和度。 目前输入盒子中,网格的 束缚水饱和度。

目前输入盒子中,网格最 大含水饱和度。

控制临界端点与深度的变 化。

输入溶剂的地面密度。

为混相驱替输入残余油饱

和度。 输入Todd-Longstaff混合

参数。

每个PVT区的分子扩散系

数。

在双孔隙度运算时,基质 与裂缝间流动进行有限分 子扩散计算。

DIFFDP

(仅选用于与双孔隙度有关 的分子扩散选择)

其他选择

VEFRAC

(仅选用于垂向平衡选择程序)

整体网格模拟时,用于地 下岩石曲线的VE函数(小 数。

目前盒子中,用于地下岩 石曲线的VE函数(小数)。 曲率参数和滞后的模型选择。

VEFRACV EHYSTR AQUTAB

TRACER

AQUTAB DENSITY DIFFC DIFFDP

EHYSTR ENKRVD ENPTVD

(仅选用于垂向平衡选择程序) (滞后选择要求) (任意选择)

Carter-Tracy含水层使用的影

响函数。

(用于惰性示踪剂跟踪) 惰性示踪剂的名称和相态。

关键字按字母顺序列表

Carter-Tracy含水层影响函数。 储罐流体密度。

每个PVT区分子的扩散系数。

双孔隙度运算时,对基质与裂缝间流动进行有限分子扩散 计算。

曲率参数和滞后的模型选择。

相对渗透率的端点值与深度的变化表。 饱和度的端点值值与深度的变化表。

储罐流体的比重。

目前输入盒子中,网格气的相对渗透率曲线端点值。 目前输入盒子中,网格油的相对渗透率曲线端点值。 目前输入盒子中,网格水的相对渗透率曲线端点值。 上覆岩层压力表。

模拟预测的最大压力。

随泡点压力变化的可动油的体积系数和粘度。 随压力变化的干气的体积系数和粘度。 随压力变化的死油的体积系数和粘度。 随压力和Rv变化的湿气的体积系数和粘度。 随压力和Rs变化的可动油的体积系数和粘度。 水的体积系数、压缩系数和粘度。

随含盐量变化的水的体积系数、压缩系数和粘度。 岩石压缩系数。 岩石压实表。

岩石滞后的压实表。

油藏性质参数的输出格式。 对全油田令死油的Rs为定值。

对不同PVT区令死油的Rs值为定值。

GRAVITY KRG KRO KRW OVERBURD PMAX PVCO PVDG PVDO PVTG PVTO PVTW

PVTWSALT ROCK ROCKTAB ROCKTABH RPTPROPS RSCONST RSCONSTT

RVCONST RVCONSTT SCALELIM SDENSITY SGCR SGFN SGL SGOF SGU SLGOF SOF2 SOF3

令全油田干气的Rv为定值。

令不同PVT区干气的Rv为定值。

对每个端点的刻度与深度范围设置饱和度标定界限值。 输入溶剂的地面密度。

目前输入盒子中,网格气体的临界饱和度。 随Sg变化的气的相对渗透率和毛管压力。 目前输入网格原始气饱和度。

随Sg变化的气的相对渗透率和毛管压力和油对气的相对渗 透率。

目前输入盒子中,网格最大气饱和度。

随Sg变化的气的相对渗透率和毛管压力以及油对气的相对 渗透率。

在两相系统中,随So变化的油的相对渗透率。 在三相系统中,随So变化的油的相对渗透率。 目前盒子中,网格残留在气中的油饱和度。 对混相驱替,输入残余油饱和度。

目前盒子中,网格残留在水中的油饱和度。 油气表面张力表。

用Stone的公式1计算三相油的相对渗透率。 用Stone的公式2计算三相油的相对渗透率。 油水表面张力表。

目前输入盒子中,网格临界水饱和度。 随Sw变化的水的相对渗透率和毛管压力。 目前输入网格束缚水饱和度。

随Sw变化的水的相对渗透率和毛管压力以及油对水的相对 渗透率。

目前输入盒子中,网格最大含水饱和度。 输入Todd-Longstaff混和参数。 惰性示踪剂的名称和相态。 控制临界端点与深度的变化。

为整体网格模拟,用于地下岩石曲线的VE函数(小数)。 在目前盒子中,用于地下岩石曲线的VE函数(小数)。

SOGCR

SORWMIS SOWCR STOG STONE1 STONE2 STOW SWCR SWFN SWL SWOF SWU

TLMIXPAR TRACER TZONE VEFRAC VEFRACV

2.5 REGIONS部分

综 述

REGIONS把计算网格分成各区进行下列计算:

(a)计算饱和度函数(相对渗透率和毛管压力);

(b)计算PVT特性(流体密度、地层流体体积系数、粘度); (c)平衡计算(设置原始压力和饱和度); (d)原始流体和区间流体流动,结果输出;

(e)计算不同方向相对渗透率;

(f)计算吸吮饱和度函数(滞后选择);

(g)岩石压实选择,计算ROCRTAB(岩石特性表)特性;

本文来源:https://www.bwwdw.com/article/sm9.html

Top