必须建立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方向上扩散率
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,免费范文网,提供经典小说综合文库Eclipse 100 油藏数模软件使用手册(6)在线全文阅读。
相关推荐: