油藏数值模拟基础培训(第四讲)ppt课件.ppt

《油藏数值模拟基础培训(第四讲)ppt课件.ppt》由会员分享，可在线阅读，更多相关《油藏数值模拟基础培训(第四讲)ppt课件.ppt（82页珍藏版）》请在三一办公上搜索。

1、油藏数值模拟基础培训（第四讲）,张世明，2009年4月，职工大学,油藏数值模拟基础培训,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,RUNSPEC/GRID要点分析,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,主要作用（分配存储空间并启动应用程序选项）,在主存储空间中顺序为模拟的各个组成部分分配内存；模拟的组成部分包括井、表格数据、模拟网格及解算器存储栈；一些关键字是必选的，大部分关键字是可选的；忽略的关键字等价于将这些关键字设为默认值。,模型大小：DIMENS 模型流体：OIL,GAS

2、,WATER,DISGAS,VAPOIL,根据油藏实际情况选用不同组合。 模型单位：METRIC, FIELD, LAB 模型维数：WELLDIMS(井，射孔和井组维数），VFPIDIMS(注水井VFP表），VFPPDIMES(生产井VFP表），EQLDIMS(平衡区），TABDIMS(流体及岩石维数） 模型输出格式控制：UNIFOUT,FMTOUT,UNIFIN,FMTIN 特殊功能：API,BRINE,DUALPERM,DUALPORO,MISCIBLE,TRACES 标题： TITLES 模拟开始时间：START,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVA

3、TION,TRANSCENDENCY,PVT和相渗曲线数据占用的存储空间相对较小；油藏模拟网格数据及解算器存储栈占用存储空间最多；在每一个事件步，eclipse就会求解一个关于每个网格的压力、饱和度、溶解气油比（或油气比）的方程组。连续的模拟迭代步的解必须是正交的，所以Eclipse需要快速访问前几个时间步的解，这就是解算器的存储栈。默认的是快速访问前10个时间步的解，也可以用NSTACK来修改这个值。,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,RUNSPEC部分常用关键字：,AQUDIM 指定解析水体的个数和大小。DIME

4、NS 指定模型的大小。这个关键字是必须的。DISGAS，GAS，OIL，VAPOIL，WATER 指定模型中的相。至少要有一相。DUALPORO 启动双孔选项。ENDSCALE 启动端点平衡选项。EQLOPTS 设置一系列选项，用来定义初始压力和饱和度。FIELD，METRIC或LAB 指定使用的单位制。一个模型的所有数据使用一个相同的单位制。FMTIN，FMTOUT 分别指定格式输入和格式输出。,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,RUNSPEC部分常用关键字：,NOSIM 关闭SCHEDULE部分的模拟。这对于在模

5、拟运算前检查关键字的错误是非常有用的。NSTACK 指定存储前面时间步结果的堆栈的长度。困难的问题需要多于10步的值，但是设置的值不要大于最大线形迭代数减一。NUPCOL 定义更新井目标的每一个时间步的非线性迭代的次数。有时设置的值需要大于默认值3。RADIAL 设定径向几何模型。RPTRUNSP 在PRT文件中创建一个关于RUNSPEC选项的列表。START 设置模拟的起始时间。这个关键字是必须的。TABDIMS 设置输入的PVT和饱和度表的大小。,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,RUNSPEC部分常用关键字：,

6、TITLE 给出运算的名字。这个关键字是必须的。UNIFIN，UNIFOUT 分别设置独立的输入和输出文件。VE 启动垂向平衡选项。VFPIDIMS，VFPPDIMS 分别定义注入井和生产井的VFP表的个数。WELLDIMS 指定模型中井和井组的个数。这个关键字是必须的。,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,一般性的两相PVT相图,A区：代表我们通常所说的黑油。线AA：代表泡点以上的油 （死油）。线BB：代表初始时高于泡点的黑油（活油）。 线CC：代表油气两相混合状态（活油）。 线DD：代表一种初始状态接近于临界点的流

7、体 线EE：代表初始状态是单相的气体（干气）。 线FF：也是初始是油气两相混合状态（湿气）。线GG：代表初始状态是高于露点的挥发油（湿气）。,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY, 死油油藏：AA油藏在整个开发阶段压力一直保持在泡点压力以上，在油藏中不会发生脱气。生产油气比是常数（脱气发生在井筒或地面）。 活油油藏：BB油藏在开发阶段压力会降到泡点压力以下，在油藏中发生脱气，溶解油气比降低，生产油气比增加（自由气会生产出来）。 干气藏：EE气藏压力在生产过程中不会穿过露点线，气藏中没有挥发油产生。 湿气藏：FF气藏压力在

8、生产过程中会穿过露点线，气藏中产生挥发油。,一般性的两相PVT相图,黑油模型可以模拟的流体类型,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,流体类型的关键字选择,黑油模拟器，如ECLIPSE，是不能显式模拟组分的变化的。如果需要模拟气体的分离或油滴的析出过程，则要间接的通过改变溶解气的GOR（Rs）和挥发油的OGR（Rv）来实现。具有这种特征的流体分别称为live oil和湿气。,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,文件输出格式汇总,RUNSPEC部分

9、要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,ECLIPSE输出文件格式主要是以文件管理形式，并可定义为有格式和无格式，也可定义为单个文件或合成文件（特别是与时间有关的结果文件）。 无格式：在缺省状态条件下，文件存储和输出为不可读，便于程序读取，但转换有困难。 有格式：可读，便于编辑、交换，但文件较大。使用定义关键字：FMTOUT、FMTIN。 单个文件：在缺省状态条件下，用户可根据需要自行定义某些文件，对中途终止运行的作业，不会影响已做完的计算步，单个文件文件总数可达999个。 合成文件：为与时间有关的文件。如：重启文件、计算总结文件。使用定义

10、的关键字：UNIFOUT、UNIFIN。已合成一个文件，对时间步无限定，但对中途终止运算的作业没有用于GRAF的文件。,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,ECLIPSE输出文件,RUNSPEC部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,ECLIPSE输出文件,.GRID或.FGRID: 网格文件（NOSIM控制输出） .EGRID: 网格文件,与GRID格式不同，文件要小的多。（用关键字GRIDFILE来控制输出类型） .INIT或.FINIT: 属性文件。（用关

11、键字INIT来控制输出） .PRT: 报告输出。 文件很大，模型处理及计算结果详细报告。（RPTGRID,RPTPROP,RPTSOL,RPTSCHED控制输出） .LOG: 后台作业时的输出报告，文件比PRT要小很多。可用于错误检查。 .DBG: Debug文件，一般不用。可用于检查ECLIPSE如何处理输入参数。 .SAVE: 用于快速重启。（用关键字SAVE来控制输出） .RFT: RFT计算结果。（用关键字WRFTPLT来控制输出) .FLUX: 流动边界。 (用关键字DUMPFLUX来控制输出) .Snnnn或.UNSMRY: 图形文件输出（在SUMMARY部分定义） .Xnnnn或

12、.UNRST: 重启文件输出（用RPTRST，RPTSOL或RPTSCHED来控制输出）,RUNSPEC部分作业练习,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,1、根据自己实际区块的流体性质，选择合适的关键字；2、改变RUSPEC中文件输出类型关键字FMTOUT和UNIFOUT，看输出文件格式的改变；3、运行概念模型后分析产生的文件类型，找出对应的控制关键字。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,主要作用,油藏的几何模型和基本的岩石属性；网格块的大小尺寸和深度，DX、DY、DZ、TO

13、PS或COORD和ZCORN；建立网格间的流体流动连接关系；每个网格的孔隙度、净毛比、渗透率（PERMX/Y/Z或PERMR、PERMTHT、PERMZ）；传导系数修改（MULT）；创建各网格孔隙体积、深度和传导系数数组；数值水体和网格水体的定义；有效网格定义，尖灭控制和最小有效网格孔隙体积定义；,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,最简化数据需求,网格大小和深度 DX或DXV，DY或DYV，DZ，TOPS 或COORD，ZCORN孔隙度 PORO渗透率 PERMX，PERMY，PERMZ 或PERMR，PERMTHT，PE

14、RMZ净毛比或净厚度 NTG或DZNET（默认为1）,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,关于净毛比NTG,储层总厚: H=h1+h2油层总厚：h1水层总厚：h2假设油层净毛比：ntg1假设水层净毛比：ntg2则储层净毛比：NTG=(h1Xntg1+h2Xntg2)/ (h1+h2)地质概念净毛比：NTG=ntg1模拟概念净毛比：NTGNTG,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,主要作用,达西方程：,适于地面条件下水平单相流动，其中KA/L项视为表征流体在

15、模拟网格间流动能力的量，即传导系数。,考虑净毛比后的传导系数定义：,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,主要作用,在离散网格中，X方向传导系数：,传导系数不能在GRID部分显示设定，但可以通过其它方式修改。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,数据排列,读取网格数据时是X方向优先，然后是Y方向，再是Z方向,网格数据读取规则,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,数据排列,笛卡儿网格：起始点

16、在左上角。在页面上，数值从左至右沿着X方向排列，从上到下沿着Y方向排列。对于3D网格，网格块和节点是按照从左到右，从后到前，从上到下的顺序排列的。网格的起始点并不是一定要从网格块（1，1，1）开始。X，Y，Z轴并不是一定要与I，J，K方向平行。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,数据排列,径向网格数据读取规则,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,数据排列,径向网格：起始点是网格模型的中心点。在R方向上，最内部的一圈是第1圈，然后圈数从内到外依次增加。角是

17、沿着顺时针方向度量的；Z方向保持不变。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,数据排列规则,对于每一个网格都有的数据，应用下列读写规则：第一个读写的网格编号为（1，1，1）。对于笛卡儿网格的显示，总是把网格块（1，1，1）放在上、后、左的位置。径向网格的网格块（1，1，1）在最接近网格模型中心的位置。在读取笛卡儿网格数据时，总是X方向循环最快，其次是Y方向，最后是Z方向。在读取径向网格数据时，总是R方向循环最快，然后是THEAT A方向，最后是Z方向。所有的GeoQuest模拟软件都使用这种数据读取规则。,GRID部分要点分析

18、,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,网格几何形状,块中心和角点网格,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,ECLIPSE支持多种网格类型。块中心（BC）网格。角点（CP）网格。每一种网格类型都各有优劣。BC或CP的选择取决于模型的类型。在同一个数据文件中不能够同时使用BC和CP网格。,网格几何形状,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,网格几何形状,块中心网格：块中心网格需要知道各网格的顶深及X、Y、Z方向上的

19、网格尺寸。各网格的上、下底面都是水平的，各侧面都是平的且垂直的。所有的网格块都是矩形的。要创建笛卡儿BC网格需要指定关键字TOPS，DX（或DXV），DY（或DYV）和DZ。要创建径向BC网格需要指定关键字DR（或DRV），DTHETA（或DTHETAV）和DZ。以V结尾的关键字指的是向量格式的关键字，是其他关键字的替代。由于每个网格块都只用四个实数来定义，所以BC网格比它响应的CP网格的数据量要少一些。用BC网格建立的简单模型甚至可以不用GRID之类的前处理软件就可以创建。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,网格几何形状

20、,块中心网格表示的含有断层的倾斜构造,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,网格几何形状,考虑相邻的BC网格块。对于倾斜结构，各网格的顶深TOPS是不同的。那么断层两边的网格的顶深TOPS也是不同的。BC网格有效地忽略了倾斜。面和断层的差别。例如，在没有构造图指导的情况下，图中的横断面可以看作是倾斜面，也可以看作是断层。BC网格没有足够的信息来计算相邻网格间的交迭面积，这是因为不知道网格块拐点处的深度。流体通常在相邻网格间流动，所以需要建立相邻网格之间的连接关系。ECLIPSE假设有相邻标号的网格（如图中箭头所连接的网格）之间

21、是有连接的，即使它们之间实际可能并没有，如图中所示的有交叠面的网格间就没有连接。这就是因为在BC网格中不区分断层和倾斜面。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,网格几何形状,角点网格表示的含有断层的倾斜构造,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,网格几何形状,角点网格：角点（CP）网格是建立在坐标线和拐点深度的概念基础上的。一个坐标线定义一个网格柱的一个边。坐标线总是直线，但不一定是垂直的。给出直线上分别在网格上下的两个点的X、Y、Z坐标就可以定义一条坐标线

22、。然后，只要确定其各拐点在各坐标线上的位置就可以确定网格块了。这种定义网格的方式使得网格块可以是任何可能的形状：倾斜表面、断层面、尖灭和剥蚀面都可以正确的表示出来。因为每一个网格块都是用4条坐标线和8各拐点坐标定义的，所以CP网格比BC网格有更大的数据量，而且几乎只能借助于如FloGrid/GRID之类的前处理软件才可以创建。角点网格包含足够的信息来计算相邻网格块间的交迭区域，因为知道各网格块的拐点深度。这意味着，在图中，只有可以看到真正有公共交接面的网格之间才有流体的交换。因此，箭头所指的网格间的流体不可能穿过断层流动。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVAT

23、ION,TRANSCENDENCY,角点网格与块中心网格,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,网格模型格块属性,定义网格模型格块属性,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,描述各网格格块的大小尺寸和深度的参数称为几何参数。用来指定几何参数的关键字是：对于笛卡儿网格是TOPS，DX（或DXV），DY（或DYV）；对于角点网格是TOPS，DR（或DRV），DTHETA（或DTHETAV）和DZ；对于角点网格是COORD和ZCORN。描述孔隙度、渗透率的参数称为属

24、性参数。描述属性的参数是：PORO（），PERMX（Kx），PERMY（Ky），PERMZ（Kz）。净厚比可以用关键字NTG（净厚比率）或DZNET（净厚度）来定义。用来显式定义一个网格是否是有效网格的关键字是ACTNUM。对于无效网格，此值为0；对于有效网格则有一个统一的值。,网格模型格块属性,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,网格模型格块属性,属性关键字是一个网格一个值。有许多种方法，可以用不同的ECLIPSE关键字的组合来来给这些关键字赋值。各属性值认为是网格块中心处的值，而且是网格块内各点处此属性的平均值。这种求平

25、均值的运算是在数据输入ECLIPSE之前做的，而且求平均值的方法是由客户自己决定的。同时注意，虽然有些网格会被排除在模拟运算之外（无效网格），但是仍然要提供计算其孔隙体积和传导系数的数据，这是因为ECLIPSE的内部工具是根据其最小孔隙体积来设定一个网格是无效网格的。ECLIPSE有一个规则，那就是所有提供的数据都必须是显式的，因为ECLIPSE没有提供输入函数形式数据的工具。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,如何指定网格格块属性,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSC

26、ENDENCY,如何指定网格格块属性,典型格式为：KEYWORDValue1 value2 value3.value(NX*NY*NZ) /,如何为每个网格格块设定一个属性值,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,如何指定网格格块属性,如何用Boxes设定网格格块属性值,输入的box是用I、J、和K的范围来定义的，是一个由多个网格格块组成的矩形块。Box可以是网格的一个区域、一层、一列或一行。可以为Box内的网格块设置所有属性值，应该为Box内的每一个网格块（包括区域端点的网格块）都指定一个参数值。一个Box会一直保持有效，直

27、到又读入了一个新的Box或读入了ENDBOX关键字。这个新的box会关闭前一个box，并打开另一个box。默认的box是模型中所有网格格块的集合。,BOX-I1I2J1J2K1K2 12011110/PERMX200*100/ENDBOX,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,如何指定网格格块属性,如何用EQUALS设定网格格块属性值,EQUALS关键字是用数组的形式对网格的属性值进行操作，可以用作关键字BOX的替代品。,EQUILS同样可以用来对第3层的网格格块设置其PERMX值：BOX-I1I2J1J2K1K2 12011

28、33 /EQUALSPERMX 100 /ENDBOX不用BOX关键字，EQUILS也可以显式定义box。例如：EQUALS-ArrayvalueI1I2J1J2K1K2PERMX2000/默认是当前打开的box- i.e. entire reservoirPERMX101201511/,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,如何指定网格格块属性,Box作为EQUALS关键字的一部分被打开，然后被退出EQUALS关键字的第二个正斜杠关闭。为了保证为每一个网格格块都指定了参数值，建议现给全体网格格块赋值，然后再用BOX和EQUI

29、LS选择性的覆盖其中某些值。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,如何指定网格格块属性,如何拷贝网格格块属性值,COPY关键字用来在当前的box内从一个数组向另一个数组拷贝数据。如果油田的渗透率是统一的，并且当前的box是整个油田时，可以这样设置其渗透率值：COPYPERMXPERMY /PERMXPERMZ /如果只有一个特定区域的关键字是各向同性的，则可以把BOX关键字和COPY关键字结合起来使用。例如：BOX-I1I2J1J2K1K2 1201511 /COPYPERMXPERMY /PERMXPERMZ /ENDBO

30、X在COPY关键字内部，对数组拷贝的个数没有限制。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,如何指定网格格块属性,如何加、减、乘、除网格格块属性值,KEYWORDValueI1I2J1J2K1K2 /KEYWORDValueI1I2J1J2K1K2 /KEYWORDValueI1I2J1J2K1K2 /或BOXI1I2J1J2K1K2 /KeywordValue /KeywordValue /KeywordValue /ENDBOX,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCEND

31、ENCY,如何指定网格格块属性,如何用MULTPV对网格格块孔隙体积做乘法,对于网格块的孔隙体积可以用MULTPV做乘法操作，例如：MULTPV200*1.01 / 能够增加前200个网格格块1的孔隙体积。这种乘法操作是很必要的，例如当要确保拟合的地下流体（FIP）与来自于其它来源的估计相一致时。然而，要注意的是，孔隙体积越大，则这个网格所提供的压力支持度就越高。因此，孔隙体积的修改不应该过大，而且对油藏全区压力的拟合应该在调整孔隙体积之后进行。在历史拟合过程中，有时需要调整孔隙体积来达到拟合一个或多个目标的目的。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,

32、TRANSCENDENCY,如何指定网格格块属性,如何用COPYBOX把一部分网格数据拷贝给另一部分网格,可以用关键字COPYBOX把一部分网格数据拷贝给另一部分网格。源box和目标box必须在I、J、K方向上有相同的大小尺寸。一般格式如下：COPYBOXKeywordSource boxDestination box /,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,如何指定网格格块属性,如何用INCLUDE读取另一个文件的数据,外部文件数据可以用INCLUDE关键字直接调入ECLIPSE数据文件。这个外部数据文件可以存在于文件系统

33、的任何地方，可以用绝对路径指定，也可以用相对路径指定。如果这个外部数据文件是在读取ECLIPSE数据文件的同一个路径内，则不需要指定其路径。 INCLUDE可以用于数据文件的任何位置，一般用于调用较大的数据体。如果路径名中含有一些会引起ECLIPSE歧义的字符，如前斜杠或句号则需要给路径加上引号。在UNIX系统中，文件名是大小写敏感的。用INCLUDE调入的文件可以是嵌套的，也就是说此文件内部还可以包含INCLUDE关键字。但是关键字INCLUDE不能够用于另一个关键字内部。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,如何指定网格

34、格块属性,如何用ACTNUM设定网格格块为无效网格,任何网格格块都可以用关键字ACTNUM显式设定为无效网格。ACTNUM是一个网格格块属性数组，因此可以用于前面的所有例子。ACTNUM只能有两个值：0表示无效网格，1表示有效网格。一般用FloGrid/GRID等前处理软件来定义ACTNUM。无效网格会从模拟中删除出去，因为ECLIPSE不需要计算这些网格块的流动，但是不会改变图形输出，因为网格的位置是保持不变的。ECLIPSE必须有足够的数据才可以计算出无效网格的孔隙体积，深度和传导系数。虽然它们是无效网格，但是仍要计算它们的PORO，PERMX等参数。无效网格在输出时给以指定，也不会为无效

35、网格生成很多相关量，而且网格是用连字符标出来。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,页岩模型,页岩的表示,每一个岩层都可以模拟为一个单独的网格层。页岩层可以合并到砂岩层中。页岩可以模拟为网格层间的间隙。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,页岩模型,将页岩显式模拟为网格层,这种方法的优点是可以精确的控制层的属性。通过改变渗透率来影响传导系数，可以控制穿过页岩的流体流动。NTG不会出现在垂向传导系数的计算中，但是会用于水平传导系数的计算，而页岩隔层的NTG几

36、乎是0，所以在页岩内部是没有水平流动的。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,页岩模型,将页岩显式模拟为网格层,页岩一般相对较薄，而且与其相邻网格相比，其孔隙体积通常很小，这会带来一些困难：,如果模型运用PINCH和（或）MINPV关键字来设定一些网格块为无效网格，或建立跨尖灭网格间的连接，则很有可能把页岩网格从模拟中除去。用MINPV产生的无效网格成了流动的阻碍。比尖灭厚度门限值薄的网格会从模拟中除去，而此网格两边的网格会重新建立相互的流体流动，就像页岩不存在一样。从效果上来说，前者把页岩转化为流动障碍，而后者则把页岩从模

37、型中除去。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,页岩模型,将页岩显式模拟为网格层,一个大孔隙体积的网格块与一个小孔隙体积的网格块相邻常常会导致流通量相关的收敛性问题。考虑两个孔隙体积分别为1Rb和1000Rb的相邻网格块。ECLIPSE求解的油相饱和度具有精度。如果已经计算出油相饱和度为0.5，则各网格的地下原油量分别为Rb和Rb。那么其中一个网格块的地下原油含量就比另一个其相邻网格块内地下原油含量的误差还要小，这就使得大网格会决定小网格的流通量，其孔隙体积是小网格孔隙体积的许多倍。这就很可能会导致收敛性问题，我们应该避免这

38、种情况的发生。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,页岩模型,将页岩显式模拟为网格层,当有很多个大范围的页岩层时，如果显式建立它们的模型就会给整个模型引来大量没有意义的网格，这将会极大的影响运算时间。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,通过将页岩合并到大的砂岩网格块中来模拟页岩,页岩模型,页岩可以合并到砂岩层中去，要实现这种方法必须要用关键字DZNET或NTG提供砂岩的净厚度，因为DZNET和NTG表达的是相同的量，所以这两个关键字不能同时使用。网格块的

39、孔隙体积应该与纯砂岩部分的孔隙度相同，而不是整个网格块孔隙度的平均值。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,这种方法为了增加效率而把页岩部分从整个网格块中除去，因此会导致不正确的网格块孔隙体积。同时，地质页岩层与其相邻砂岩层间的垂向传导系数也不再是0，页岩所在位置的水平传导系数也不再是0，而用的是砂岩区域的水平传导系数。因此，所有这些传导系数，尤其是垂向传导系数都需要用一个或多个“MULT”关键字来进行修正。这种做法的优点是消除了流通量相关的收敛性问题，并减少了模拟中网格块的数目。而且事先不知道用MINPV和MINPVV会把

40、哪些网格变为无效网格所带来的困难也消除了。,通过将页岩合并到大的砂岩网格块中来模拟页岩,页岩模型,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,页岩模型,把页岩模拟为砂岩层间的间隙,使用一些建立网格的前处理软件，如GRID，是可以把一个层与另一个层分离开的。层K=1与K2之间的空间并没有被网格占据，而是一个无效空间。然而K=1与K2之间的垂向渗透率并不是0，而是用通常的计算方法来计算的，层间的间隙对它并没有影响。从效果上来说，页岩层已经不存在了，必须用MULTZ和（或）MULTZ-调整层间的垂向渗透率，以减小或消除垂向流动。对于层1、

41、2的水平传导系数和孔隙体积分布来说，同显式把页岩模拟为网格的处理方法一样，保持不变。这种模拟页岩的方法的优点是可以减少模拟的网格数目，消除流通量相关的收敛性问题，并且可以不依赖于页岩孔隙度和渗透率而对层间的垂向流动进行调整。但这种方法的主要缺点是当页岩层的含油量较高，并随着时间不断被提取出来时，这部分油不会被包含在模拟过程当中。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,非相邻连接NNC,非相邻连接的来源,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,非相邻连接NNC,N

42、NC允许没有相邻IJK标号的网格间建立相互流动关系。当不相邻的两个网格之间有流体的流动时，就要用到NNC。默认的是产生NNC，在RUNSPEC部分用NONNC关键字可以关闭NNC的生成。在下面一些情况中会自动产生NNC和其传导系数：断层处，当使用了关键字NEWTRAN时。尖灭和侵蚀面，当使用了PINCH或PINCHOUT关键字时。局部网格加密和粗化时。双孔模型。当在径向模型中“完成整个圆圈”时，COORDSYS关键字会激活NNC的生成，ECLIPSE会计算它们的传导系数值。水体的NNC必须显式创建，ECLIPSE会计算它们的传导系数。关键字NNC可以为模型中的任意一对网格块显式建立非相邻连接，

43、但是这些非相邻连接网格间的传导系数也需要显式设定。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,非相邻连接NNC,我们通常认为相邻网格之间存在流体流动，而且会将其模拟为有相邻（I，J，K）标号的网格块之间的直接连接。然而，在对许多构造的精确定义中，常常需要把没有相邻（I，J，K）标号的网格块放在相邻的位置上，而且实际相互之间是有流体流动的。穿越大落差断层的流动就是一个典型的例子。其它需要用到NNC的情况就是由于ECLIPSE所用几何描述方法的本质特点了，如局部网格加密（LGRs）和水体的连接。,GRID部分要点分析,求实 创新 超越

44、,REALITY，INNOVATION,TRANSCENDENCY,非相邻连接NNC,断层处NNC的生成,断层处的非相邻连接,断层NNC会自动生成，除非在RUNSPEC部分用了NONNC关键字。如果使用NEWTRAN，则会自动生成NNC。自动计算传导系数。,GRID部分要点分析,求实 创新 超越,REALITY，INNOVATION,TRANSCENDENCY,非相邻连接NNC,断层处NNC的生成,NNC的使用是默认的，会在断层两边必要的地方自动创建。对于角点网格，默认使用的是NEWTRAN传导系数计算方法，NNC的生成不需要用户的干预。对于角点网格，我们不推荐使用OLDTRAN方法，因为它会

45、禁止在断层处生成NNC。块中心网格默认使用的是OLDTRAN传导系数计算方法，由于ECLIPSE缺少判断几何信息来判断断层两边的网格是否真的相邻，所以不会在断层处生成NNC。虽然对于块中心网格可以使用NEWTRAN，但是这样做并不会提供额外的网格几何信息。,非相邻连接NNC,尖灭处NNC的生成,/ 孔隙体积5000的网格会被设为无效,GRID部分要点分析,非相邻连接NNC,ECLIPSE允许层厚度为0，可以用来模拟尖灭和侵蚀面。如果没有用关键字NONNC，关键字PINCH会为DZ小于设定门限值的层两边的网格块建立非相邻连接。如果网格块已经被MINPV设定为无效网格，则PINCH中仍有选项可以把

46、DZ大于门限值的层两边设定为非相邻连接。可以用PINCHREG和PINCHNUM基于一个区域定义尖灭。,尖灭处NNC的生成,GRID部分要点分析,非相邻连接NNC,角点网格允许网格具有非矩形的形状，这对于精确的模拟尖灭和侵蚀面是一个非常有效的工具。虽然尖灭和侵蚀面的地质构造差别很大，但是从模拟网格的角度来说，它们可以用相同的方法来定义。PINCH关键字会在厚度小于指定值的无效网格两边建立非相邻连接NNC。例如，PINCH/会在任何DZ小于0.01英尺（或米）的网格的上下两个相邻网格之间建立NNC。从98a版本之后，就可以基于一个区域定义其中的尖灭了。关键字PINCHNUM就是用来标识由PINC

47、HREG定义的尖灭区域的。,尖灭处NNC的生成,GRID部分要点分析,非相邻连接NNC,尖灭处NNC的生成,然而，除了这些用户用关键字自己定义无效网格的方法以外，ECLIPSE还会自动的把很多网格设置为无效网格。为了避免流通量相关的收敛性问题，ECLIPSE会把所有孔隙体积小于10-6 m3、Rb 或 cc的网格都设为无效。可以用关键字MINPV增大此限制值，但是不能再减小了。例如，MINPV5000 /会把一个公制单位制的模型中所有孔隙体积小于5000Rb的网格块都设为无效网格。这些网格不会参与模拟运算，也不会求解它们的流动方程。默认的这些网格内是没有流动的，是流体流动的障碍，但是并不能保证

48、所有这些网格都会比PINCH所设定的无效网格厚度门限值要薄。如果PINCH的第二项设为“GAP”，则在被MINPV设定的无效网格上下两边也建立非相邻连接。若是这样，则对于模型中那些孔隙体积大于MINPV，但是厚度小于尖灭门限值的网格块来说，它们也会被PINCH设为无效网格。,GRID部分要点分析,非相邻连接NNC,局部网格加密（LGR）处NNC的生成,LGR NNC生成,对于一个LGR，ECLIPSE自动在生成一个全局主网格块与加密网格块之间的NNC表。,GRID部分要点分析,非相邻连接NNC,局部网格加密（LGR）处NNC的生成,每个局部网格加密内的网格块用的都是的此LGR内的局部I，J，K

49、标号。LGR网格块很少会和与其相邻的全局主网格有相邻的网格标号。对于包含LGR的模拟，ECLIPSE会报告LGR网格块与其相邻全局网格块间的NNC。从严格意义上来说，这并不是真正意义上的NNC，而只是为了输出报告的方便，因为LGR网格与全局网格之间是分开求解的，在它们之间并没有直接的流动。对LGR的处理与ECLIPSE 200 FLUX BOUNDARY OPTION的处理是非常类似的，这两方面内容在GeoQuest提供的其它教程中都有详细的分析。,GRID部分要点分析,非相邻连接NNC,双孔模型中的NNC,双孔模型中NNC的生成,在双孔模型中，ECLIPSE会把层数增加一倍。上半部分是基质网

50、格块，下半部分是裂缝网格块。基质网格块和其对应的裂缝网格块之间会自动的用NNC建立连接。,GRID部分要点分析,非相邻连接NNC,对于双孔双渗油藏，ECLIPSE在模拟的时候是把基质和裂缝分开的，分别将其模拟为单独的网格组，而它们之间是通过自动生成的NNC建立连接的。用户可以像往常一样建立油藏的网格模型，但是在RUNSPEC部分需要分配两倍的层数。而后面关键字DPGRID的使用则会指导ECLIPSE把模型网格的层数增加一倍，然后把现有层的属性拷贝到新增加的层上去。上面一层被处理为基质网格，而下面一层则被处理为裂缝。用户必须提供裂缝网格块的属性、基质裂缝耦合系数以及基质和裂缝分别的饱和度函数。在