CMG中IMEX中文操作手册(一)

1、 目 录 imex 的介绍1指导段5关键字输入系统中的数据段5如何建立数据文件的文档方式6如何执行重新启动运行7控制打印输出文件的内容8图形文件(sr2)的内容控制9网格系统描述10局部加密网格描述11双孔/双渗模型的使用12死结点的输入17水区选项的使用18拟混相选项的使用19注聚合物模型的使用20程序执行时间长或时间步长太小问题21单相油藏的模拟23水平井25垂直平衡计算25定义多个 pvt 区29井的定义30井的类型定义31如何关井以及重新开井32亏空填充(voidage replacement)34井筒流动模型的使用36操作及监测限制37井指数的输入39中止模拟运行41在井的列表中使用

2、通配符41关于井数据段设置的指导42在循环数据段内可用的其他段关键字44关键字数据输入系统46关键字系统介绍46注释行 (可选择)52空行52包含文件53控制数据文件列表54基岩网格性质的输入57裂缝网格性质的输入57加密网格性质的输入58由 i 方向确定 j 和 k 方向的数据59常数值数组59以 ijk 方式输入数组60数组输入值沿 i 方向变化61数组输入值沿 j 方向变化62数组输入值沿 k 方向变化62大多数或所有网格的值都是不同的63输入/输出控制段67基本网格定义99k 坐标方向105i 方向的网格尺寸107j 方向的网格尺寸108k 方向的网格尺寸110网格的顶部深度112加密

3、网格的位置123双孔介质minc方法126死结点标识符128孔隙度130岩石压缩系数131渗透率134压力影响函数147流体组分性质数据段153流体模型153油藏温度154油和气的 pvt 表155作为压力函数的油压缩系数168密度173水相粘度179原始气油比180聚合物粘度混合（条件）183岩石流体数据段186相对渗透率表186油水相对渗透率表187气液相对渗透率表194评价三相共存时油相相对渗透率的方法（任选）203初始条件数据段209初始条件标识209油藏初始油相压力212初始泡点压力212初始泡点压力与深度关系213初始含油饱和度215初始含水饱和度216初始聚合物浓度216参考深度

4、和参考压力217油水界面深度218数值计算方法数据控制段225井数据段241井的改变日期241井的改变时间242设定井底流压初始化的频率248附录32065 imex 的介绍简介 imex 是一个考虑重力及毛细管力的三相黑油模拟软件,网络系统可采用直角坐标,径向坐标,变深度/变厚度坐标,在任何网络系统中.都可建立两维或三维模型.在处理气相的出现及消失情况时,程序采用了变量替换方法.imex 的一些特征和功能为：自适应隐式方法 imex可以在显示,全隐式以及自适应隐式三种方式下运行.在大多数情况下,只有很少一部分网格需要采用全隐式求解,而大部分网格都可采用显式方法求解.自适应隐式方法正是适合于这

5、种情况的解法,并且在井附近以及层状油藏的薄层中,开采时会产生高速流动的锥进问题,采用自适应隐式处理这类问题是很有效的. 采用自适应隐式选项可节省三分之一到一半的运行时间。计算时可采用和全隐式方法同样大的时间步长.用户可以指定采用全隐式方法计算的网格，可根据用户确定的界限或矩阵转换临界值，动态地选择采用全隐式计算的网络网格。双孔/双渗 双孔隙度选项允许采用两种方法对基岩模型进行离散化处理，其中一种为嵌套格式，成为“多重内部作用连续域”(minc)方法，另一种为层状格式，称作“子区域”方法。双孔隙模型对裂缝油藏进行了理想化的近似处理，认为裂隙油藏由两部分组成：主要孔隙度和次要孔隙度， 主要孔隙度(

6、基岩)代表岩块中的微小粒间孔隙，次要孔隙度(裂缝)由裂缝，通道和溶洞组成。双孔隙模型将油藏分为两个连续域，裂隙是流体流动的主要通道，只具有很小的储集性能；而基岩具有较低的流体传导能力，但具有较大的储存能力。使用单个基岩岩块/一个裂缝系统就可确定一个简单的双孔隙模型。形状因子的计算基于 warren 和 root 或 gilman 和 kszwmi 的工作。 在这种情况下，假设基岩和裂缝是在半稳定流状态下进行传输的。为了恰当地表示基岩与裂缝间的传导机理，有必要将基岩进一步划分为较小的体积单元。minc 法通过将基岩划分成一系列嵌套的体积单元来达到这种目的，离散化时基于这样一种假设，即体积单元表面

7、为等势面，通过一个裂缝网格时，主要变量的变化很小。 minc 模型对于通过基岩一裂缝的非稳定流动是一种很好的表示方法，但是当基岩的高度很大时，不能适当的处理重力影响，是模型的一个严重缺陷。一种可以表示重力影响和流体相态分异现象的改进模型是子区域方法，在这种模型中，对基岩中的流体和压力分布进行计算，但假设重力分异现象只存在于裂缝之中。 除了 minc 模型和子区域模型外，双渗模型允许基岩之间的流体传导，当裂缝的连通具有方向性时，或者当基岩之间的连通性显的重要时，例如气一油的重力驱替过程，采用双渗模型进行计算是十分有用的。拟混相驱选项 拟混相驱模型可模拟一次接触及多次接触混相过程。并能模拟水相中的

8、溶剂及驱动气。聚合物模型选项 聚合物模型选项可模拟岩石吸附，聚合物不可进入孔隙体积，阻力系数，混合粘度的扩散情况。断层油藏选项 可模拟具有一条或多条断层的层状油藏，可精确地计算错断层间的流体流动过程，模拟执行过程中不带有任何不真实的性质平均现象。不平行于坐标轴的下滑断层也能进行模拟计算。井的全隐式处理 采用了一种非常有效的方式处理井的问题，井底压力以及井网格的变量使用全隐式方法求解。假如一口井完井时射开一个以上的层，在求解时将采用充分偶合方式计算井底压力，这样便消除了层状油藏射开多层完井情况下的不收敛问题。同时 imex 软件还提供了完整有效的井控制功能，可对井底和井口压力，产量、油气比等进行

9、限制。当达到某一限制时，可根据用户要求选择新的限制条件。矩阵求解方法 imex 软件使用一种基于不完全高斯消去法的先进求解过程，作为正交加速的预处理部骤，并且专门开发了用于自适应隐式贾可比矩阵的线性方程解法。如需要得到进一步的信息请查阅 aimsol 技术手册。 对于大多数模拟应用，imex 选择的缺省控制值选项都能使计算能有效地执行，所以用户并不需要具有矩阵求解方法方面的详细知识。局部网格加密 imex 中包括局部网格加密功能，用户可以确定油藏中的一个区域进行细分，程序将自动计算出内部网格的连通情况及传导率，这种功能可在油田规模的模拟中，用于研究井周围地区的各种影响，而不需对井进行拟函数处理

10、。使用这种技术也可有效地模拟静态裂缝情况。加密所产生的额外项都可在矩阵求解程序中得到正确处理。尖灭层 使用尖灭选项可以有效地模拟尖灭情况，程序将尖灭网格从有效网格表中移走，自动地将尖灭层的上下层偶合在一起。油藏初始化 可利用给定油气和油水界面的方式建立油藏的初始情况。灵活的网格系统 在imex中可使用几种网格选项，笛卡尔坐标，圆柱坐标和变厚度/变深度网格，在以上任何一种坐标中，都可使用二维或三维系统。可变泡点使用了一种严密的采用变量替换方法的在可变泡点公式，在低饱和区内，对使用不同 pvt 数据油的混合情况能进行适当的模拟。水区模型可采用两种方式对水区进行模拟，可采用在边界增加一定的含水网格的

11、方法，或者使用 carter 和 tracy 提出的水区解析模型。当了解水体大小和位置，并且对于油藏可用包含相对少量的附加网格的方式就可进行模拟的情况下，前一种方法是很有用的；对于非常大甚至无限大的水体，后一种方法更加有用，在这里需要对进入油藏的水侵量做近似计算，而通过附加油藏边界网格的描述方法是不可行的。输入/输出单位 国际标准单位(si)，矿场单位(英制)和试验室(修改的si)三种单位制可供选择。可移植性 imex 的程序编写采用了标准的 fortran77 语言，曾在各种各样的硬件平台上运行过，这些包括：ibm大型机，cdc，cray，honeywell，dec，prime， gould

12、，unisys，apollo，sun，hp，fps 和 ibmpc 386 和 486 的兼容机等。绘图系统cmg 的绘图系统 results，将 sr2 文件系统用于模拟输出的后处理。results 也可用于包括网格设计在内的输入数据准备工作。 指导段 _ 介绍 _ 这一指导段是用于对缺乏经验的用户在关键字输入系统方面的指导，它并不能替代这个文件中的用户参考手册。在这一指导段中只讨论了一些特殊的关键字和问题，虽然当你在建立数据文件时可能会遇到指导段中“如何做”的问题，但用户手册中包含有每个关键字的详细描述。imex 使用你初始建立的数据文件，然后生成三个或四个其他文件。imex 的每次运行产

13、生一个输出文件 (out)，一个结果索引文件 (irf) 和一个主结果文件 (mrf)。此外根据用户的选择可以产生或不产生一个可回绕结果文件 (rrf)。 +-+ | | - .out (输出) | | | | - .irf (索引) 数据文件 - | imex | | | - .mrf (主) | | | | - .rrf (可回绕) +-+ (可选择) 假如要求进行重新启动运行，那么将需要几个现存文件并产生另外三个文件，说明见下图： +-+ 数据文件 - | | - .out (输出) | | input.irf - | | - .irf (索引) | imex | input.mrf -

14、 | | - .mrf (主) | | input.rrf - | | - .rrf (可回绕) (可选择) +-+ (可选择)关键字输入系统中的数据段 _当你使用关键字输入系统建立数据文件时，有几点应当记住：a)在关键字输入系统中，有7个不同的数据组。b)这7个数据组必须遵循下列顺序输入: 输入/输出控制 油藏描述 组份性质 岩石，流体数据 初始情况 数值计算方法控制 井的数据和循环数据。c)属于各组的关键字不能出现在其他组内，除非是在某些特殊情况下。通常是在循环数据组中某些别的数据需要改变。d)同样，也应注意关键字在组内的输入顺序。 如何建立数据文件的文档方式_使用这些关键字建立你的数据文

15、件文档方式： a) *title1, b) *title2, c) *title3 d) *caseid. 它们必须位于输入/输出控制段。这些关键字是可选择的，可以不包括在数据文件中，然而它们对于文件的文档化以及区别相似的数据文件是十分有用的，至少应使用一个标题，所有的标题和情况标识必须包括在单引号之间。 *title1 和 *caseid 均用于 sr2 文件系统，这一系统用于产生模拟结果的图形，*title1 最长可以是40个字符，而 *title2 和 *title3 每个最多允许为80个字符。你也可以使用两个关键字标识符即“*”在你的数据文件中插入注释，注释可以出现在数据文件中的任何地

16、方。 e.g. *title1 simulation run #1 - 1989-01-23 *title2 dual porosity problem using the minc option *title3 this is a 12 x 12 x 10 cartesian grid system *caseid run1 * 如果这些标题行没有足够的空间用于这个数据文件 * 的文档化，你也可以使用注释来描述你输入的数据。 * 这样你就可以在这里或其它任何地方加入附加信息。如何执行重新启动运行_ 什么是重新启动文件？重新启动文件是一套sr2文件，它至少包括一个结果索引文件(irf)和一个主

17、结果文件，这套文件可以包括也可以不包括一个可回绕重启文件，这取决于用户使用的选项，关于模拟变量的重新启动信息（通常称作重新启动记录）根据用户选择的频率写入这些文件中，只能在写有重新启动记录的时间步长进行重新启动运行，是否写重新启动记录是可选择的。为什么你需要做重新启动运行？由于下述原因你可能想重新启动：a) 做敏感性研究或历史拟合，b) 改变井的工作制度，c) 在做一个较大的，时间较长的作业之前，执行一个短的模拟运行以观察结果是否满意。d) 为在一系列运行中节省执行时间，例如：你已经完成了一个模拟运行并且初步结果看上去不错，现在你想进行预测运行。因为在初始运行中你已经生成了重新启动记录，你可以

18、从你的运行中选择一个时间步并且重新启动模拟，模拟软件并不需要在开始时启动，而是从你选择的时间步开始继续执行。如何进行重新启动运行？重新启动记录是可选择的并不是必须要写的，但是如果你计划做重新启动，你就需要在初始运行中生成重新启动记录。使用 *wrst 和 restart 关键字就可生成重新启动文件， 它们必须位于输入/输出控制段中，然而当井的情况改变时，*wrst 也可出现在井的数据段中。 *wrst 用于说明写重新启动记录的频率。*restart 说明当前的模拟运行是一个重新启动运行。假如你希望在最后一个时间步开始重启，那么就可在 *restart 这一行的关键字后其余部分保持空缺，这与缺省

19、情况相同，否则输入一个时间步数。 例： *restart 30 *wrst 10 为了执行一个重启运行：a) 除了将 *restart 加在数据文件中的输入/输出控制段之外，不需要对原来的油藏数据做任何修改。b) 如果必要的话，增加时间步的最大值，或者不与 *maxsteps 一起考虑。c) 你需要输入文件(在初始运行期间生成的)，包括一个 irf 输入文件，一个 mrf 输入文件，并且可能还有一个 rrf 输入文件。d) 你也需要为新的一组输出文件命名，确保所有的输入文件都是由相同的模拟运行所产生的。输出文件的控制内容控制打印输出文件的内容_控制输出文件内容的打印时，需使用关键字： a) *

20、wprn b) *outprn. 这些关键字必须出现在输入/输出控制段中，也可能在以后参数改变时，出现在井的数据段中。 *wprn 控制打印频率，打印内容包括网格数据，段和租赁面数据，井的数据和计算方法控制数据，例如牛顿迭代和时间步的收敛情况。如果不要求在输出时打印网格及井的数据，则将打印频率设置为零。 例： *wprn *well 0 *wprn *grid 0 *wprn *sector 0 假如它们之中的某些不出现在数据文件中，缺省情况是在循环数据段中每个 *time 或 *date 关键字所表示的时间阶段打印 *grid 和 *sector 信息。对于 *well 缺省是在每个时间步打

21、印信息，这种缺省将会产生非常大的输出文件，能很快地将计算机设备的可用空间添满。*outprn 用于限制打印什么样的井数据，网格数据，油藏数据，以及多少性质表。实际上你可以列出你想打印的网格数据类型。井数据的处理差别相当大，你可以打印出关于井的全部数据或仅仅打印井的总结，使用 *outprn *well *layer 可打印出所有井的各层信息，使用 *outprn *well *reservoir 打印油藏条件的生产情况，使用 *outprn *well *brief 可对每口井印出一行总结，这也是缺省条件。 图形文件(sr2)的内容控制 _ 为控制 sr2 文件系统的内容，使用： a) *ws

22、rf b) *outsrf. 这些关键字必须出现在输入/输出控制段中，也可能在以后参数改变时，出现在井的数据段中。 *wsrf 控制网格数据，段数据和井数据的输出频率，这些信息可以以一定的频率输出。如果不要求将网格及井的数据输出到 sr2 文件，则将频率设置为零。这可以用于减小文件的输出量，然而你也可以在随后的井改变时对输出频率进行修改。如果 *wsrf 不存在于数据段内，那么缺省条件为在每个时间步打印井的信息，当在井和循环数据段遇到 *time 关键字时则输出 *grid 和 *sector 信息。可通过使用 *outsrf 关键字控制输出信息量。关键字组合 *wsrf *grid 0 确保

23、了在模拟结果文件中不写任何网格信息，在以后井的数据改变时可对其进行修改。*outsrf 用于限制输出那种网格数据，你可以要求在给定的网格区域输出确定的变量，对于网格信息和井的信息可使用分别的变量表。 网格系统描述_对网格系统进行描述时，你需要使用： a) *grid, b) *kdir, c) *di, d) *dj, e) *dk, 使用 f) *depth 或 g) *paydepth 或是 h) *dtop, 和 i) *dip.在上述关键字中，只有 *kdir 和 *dip 是可选择项，可以不出现在数据段内。上述列出的关键字必须出现在油藏描述数据段中，并且位于 *null 和 *por

24、 关键字之前。*grid 描述使用的网格系统类型，有三种选择：规则直角坐标，变深度/变厚度坐标，径向圆柱坐标。每种选择都要求输入各方向的网格数。 例： *grid *cart 10 10 6 例： *grid *vari 10 10 6 例： *grid *radial 10 1 151)描述一个直角坐标，网格为10x10x62)描述一个变深度/变厚度坐标，网格与1)相同。3)径向坐标系统用于锥进研究，网格为10x1x15。*di, *dj 和 *dk 是要求的关键字，用于输入网格尺寸，在输入时必须使用数组读人选项。 例： *grid *cart 10 10 12 *di *con 100.0

25、 *dj *con 100.0 *dk *kvar 25.0 2*50.0 3*40.0 75.0 3*40 2*50例中采用的是规则的笛卡尔网格系统，i 和 j 方向的每个网格都是100米宽，而 k 方向的每个层都具有相同的厚度，而各层的厚度互不相同。应当注意到，当使用 *kdir *up 关键字时，你的数据是从最底层开始输入的。 局部加密网格描述_ 使用 *refine 关键字描述加密网格的位置，也可和 *hybrid 一起使用描述杂交加密网格，杂交网格是在笛卡尔网格系统中用径向网格细分。关键字 *refine 必须出现在油藏描述段中，并位于 *null 和 *por 关键字之前。*ref

26、ine 要求给出加密网格在三个坐标方向的细分网格数。使用 *hybrid 关键字，允许你将一个直角坐标下的井网格，使用径向坐标进行细分。在网格加密时，必须使用 *range 关键字指出加密网格的范围。例如：在一个 10x13x3 的规则笛卡尔坐标系统中，要将网格(1, 1, 3 )进行加密，在三个坐标方向的加密网格数分别为2,3,2.其表示方式如下: *refine 2 3 2 *range 1:1 1:1 3:3对一个基本网格在任一坐标方向最多可划分为4个加密网格，如果你对不同范围内的网格采取不同的细分形式，可顺序输入 *refine 关键字，并指明基本网格的正确位置。对于使用 *hybri

27、d 关键字进行杂交网格加密，如果在不同范围采用不同的细分形式，两个细分区域的相邻关系将受到某种限制(详见 *refine 关键字)。在某些情况下，细分程度不同的区域之间至少要有一个不细分的基本网格将两个区域分开。同时也要注意到，在双孔隙度条件下，不允许进行网格局部加密。在变深度，变厚度坐标下可以进行局部网格加密，然而，对于每个网格来说，加密网格的划分厚度对于每个基础网格是相同的。如果加密网格的性质与对应的基本网格性质不一致，使用关键字 *rg 描述加密网格的性质，否则，认为两者具有相同的性质。双孔/双渗模型的使用_可以通过使用双孔隙度模型进行天然裂缝性油藏的模拟，在这些模型中，油藏被离散化为两

28、个搭配在一起的连续域（位于同一空间的两套网格，其中一个称作基岩，另一个叫裂缝。假设基岩连续域是由在空间被裂缝隔离的基岩组成，这些基岩的尺寸在整个油藏中是变化的，并且是裂缝间隔，方位和宽度的函数。在计算的网格中，基岩与裂缝之间的传导由一个单独的流动项表示，通过油藏的流体流动由裂缝系统产生，这是流体流动的主要通道，在一个有裂缝的网格，井始终是与裂缝连接的。imex有四种不同类型的双孔隙度模型可以使用。当使用任何一种双孔隙度选项时，基岩和裂缝的性质都必须确定，这通过在主关键字后使用 *matrix 或 *fracture 修正因子进行。例如：确定基岩和裂缝的孔隙度和渗透率，基岩孔隙度是 0.10 而

29、裂缝孔隙度为 0.05，基岩渗透率为 50.0 md，kv/kh = 0.10；裂缝渗透率为 500 md，kv/kh = 0.10。 在油藏描述段输入如下: *por *matrix *con 0.10 *por *fracture *con 0.05 *permi *matrix *con 50.0 *permj *matrix *con 50.0 *permj *matrix *con 5.0 *permi *fracture *con 500.0 *permj *fracture *con 500.0 *permj *fracture *con 50.0 以相似的方式，可以完成对于基岩

30、和裂缝网格设定相对渗透率表或采用impes/全隐式。例：在开始将基岩网格设定为impes方式而裂缝网格采用全隐式方式，在井和循环数据段第一个 *date 关键字之后，输入如下： *aimset *matrix *con 0 *aimset *fracture *con 1 此外在确定使用的双孔隙度选项类型时，也必须确定计算形状因子的方法，imex 有两种有效方法用于计算形状因子，它们分别基于 warren 和 root 模型或 gilman-kazemi 模型，请参考用户手册和技术手册的有关段落近一步了解形状因子的两种计算方法。在数据文件的油藏描述段中确定双孔隙度选项类型之后设定 *shape

31、 *wr 以激活 warren 和 root 方法；相似地使用 *shape *gk 激活 gilman-kazemi 方法。为了计算形状因子，imex 需要裂缝的几何形状信息，这通过确定裂缝在三个主要方向(i,j,k)的间隔得到，裂缝间隔在整个油藏的空间内可以是变化的，油藏中的某些部分可以是双孔隙度而其他部分可通过以设置裂缝孔隙度为零的方式变为单孔隙度。在数据文件的油藏描述段中使用关键字 *difrac, *djfrac 和 *dkfrac 确定三个主要方向上的裂缝间隔，注意如果确定一个非常大的裂缝间隔则说明形状因子，sigma 较小，因而基岩裂缝的传导将会降低；确定小的裂缝间隔值表明基岩含

32、有许多裂缝，因而基岩裂缝的传导将会增加。确定裂缝间隔为零说明没有裂缝面垂直于那些坐标轴。注意裂缝间隔仅影响基岩裂缝的传导项，而不影响基岩与裂缝的连接数。为了将油藏的某一区域设定为单孔隙度，必须将网格的裂缝孔隙度置为零。注意为了将基岩或裂缝设为死结点，两者孔隙度都必须设为零，下面的用法是不可能的： *null *matrix *con 0 或者 *null *fracture *con 0这是因为 *null 应用于网格空间(在这种情况下为基岩和裂缝两者)而不是分别对基岩和裂缝网格(有效网格)。 例：一个由 10x10x5 个网格组成的油藏模型。所有层都有裂缝，裂缝间隔在 i，j 方向为 50

33、英尺，而在 k 方向为 10 英尺。在油藏描述段中输入： *difrac *con 50.0 *djfrac *con 50.0 *dkfrac *con 10.0层 3 和 4 不包含裂缝(可用单孔隙度系统代表)，相应地对数据进行修改，假设裂缝孔隙度为 0.05。 *por fracture *ijk 1:10 1:10 1:5 0.05 1:10 1:10 3:4 0.0四种不同的双孔隙度选项描述如下。在数据文件中只能使用一种类型的双孔隙度选项。 a) 标准双孔隙度 -油藏的流体流动通过裂缝网络产生，基岩网格的基本作用是作为源汇项，在油藏描述段通过使用关键字 *dualpor 激活这一模型

34、。例如：使用标准双孔隙度模型并用 gilman kazemi 方法计算形状因子。在油藏描述数据段中输入如下： *dualpor *shape *gk b) 双渗透率 - 这个模型与标准的双孔隙度模型的不同在于每个基岩网格与裂缝网格和周围基岩网格两者都连接，那么流体流动既通过裂缝网络也通过基岩网格。双渗透率选项存在连续的毛细管现象的情况可能是重要的。在这些情况下，双渗透率对于存在自由气，油并且沿深度有很大变化的油藏经常是重要的，通常只有在垂直 k 方向，基岩对基岩的传导才是重要的。通过在油藏描述段使用 *dualperm 激活双渗透率选项，通过使用基岩传导率乘子可减小基岩对基岩的流动或将其设为零

35、。例如：使用双渗透率模型以及 gilman kazemi 方法计算形状因子，仅在垂直方向设置基岩对基岩的传导。 在油藏描述数据段输入如下数据： *dualperm *shape *gk *transi *matrix *con 0.0 *transj *matrix *con 0.0 *transk *matrix *con 1.0 c) 子区域 - 这个模型是标准双孔隙度选项的一变形。它允许用户在垂直方向对基岩网格进行细分，其目的是为了更精确地表示由基岩到裂缝的重力排驱过程，这在垂直方向是很重要的。每个网格以确定数目进行细分，以此更精确地表示基岩网格内的流体压力和饱和度。对于裂缝并不需要更高

36、的分辨率，因此仅对基岩网格进行细分。在油藏描述段通过 *subdomain ndiv 激活子区域选项，这里的 ndiv 是每个基岩网格的细分层数。注意 ndiv 在这个油藏内是常数。例如：将上面双孔隙度的例子修改为子区域方式，假设垂直细分数为 3。 在油藏描述数据段： *subdomain 3 *shape *wr 假设油藏模型由 10x10x5 个网格组成，所有层都有裂缝，裂缝间隔在 i，j 方向为 50 英尺，而在 k 方向为 10 英尺。将下述数据加入油藏描述段： *difrac *ijk 1:10 1:10 1:5 50.0 *djfrac *ijk 1:10 1:10 1:5 50.

37、0 *dkfrac *ijk 1:10 1:10 1:5 10.0 d) minc 方法 - 对于标准的双孔隙度概念的另一种扩展是由 pruess 和 narasimhan 提出的多重作用连续域(minc)方法。minc 方法的主要优点是使用基岩网格的嵌套式离散化进行基岩与裂缝的传导计算。这可以十分有效地描述流体的非稳定状况，而这在标准的双孔隙度模型中是经常被忽略的。这种嵌套式离散方法是一维的，可描述压力，粘滞力和毛细管力，然而在这种一维的基岩嵌套式细分方法中不考虑重力。通过在油藏描述段设定 *minc ndiv 激活 minc 选项，这里的 ndiv 是每个基岩网格的嵌套层数。注意 ndiv

38、 的值在整个油藏中是常数。例如：将上面双孔隙度的例子修改为 minc 方式，假设嵌套细分数为 3。在油藏描述段中加入： *minc 3 *shape *wr 假设油藏模型由 10x10x5 个网格组成，所有层都有裂缝，裂缝间隔在 i，j 方向为 50 英尺，而在 k 方向为 10 英尺。 在油藏描述段中加入下述数据： *difrac *ijk 1:10 1:10 1:5 50.0 *djfrac *ijk 1:10 1:10 1:5 50.0 *dkfrac *ijk 1:10 1:10 1:5 10.0 总结：为了激活双孔/双渗选项，你可以使用： a) *dualpor b) *minc,

39、c) *subdomain, d) *dualperm, e) *shape, f) *difrac, g) *djfrac, and h) *dkfrac 这些关键字必须出现在油藏描述段，这些不同的选项都是有效的，但在任何给定的数据文件中只能使用其中一种。 如果使用了其中任何一种关键字，应将它们放在 *null 和 *por 关键字之前。 死结点的输入_在给定的网格系统中，有两种方式可指出存在的死结点:a) *null b) *por它们必须出现在油藏描述数据段中， null 是可选择项，而 por 无论如何必须出现在数据段中。*por 是用于确定孔隙度的关键字。使用 *null 时，用 0

40、 表示死结点，用 1 表示有效结点。在下面的例子中，表示 i 方向的网格 1 到 4 ，j 方向的网格 1 到 3，k 方向的网格1 到 3 范围内是死结点，其余为有效结点。对于这个例子可以使用 *ijk 数组读人选项。 例: *null *ijk 1:10 1:10 1:3 1 1:4 1:3 1:3 0可以看出上述例子中的第二行覆盖了第一行。 *null 是可选择项，如果不输入，则假设所有结点为有效结点。*null 的设定始终覆盖 *por 的设定，假如对于一个被 *null 设定的死结点，即使用 *por 对此结点设置一个非零的孔隙度，也将被程序忽略。当有效结点设为非零值时，也可通过 *

41、por 用零孔隙度指出死结点。前面的例子同样可以使用 *por 代替: 例: *por *ijk 1:10 1:10 1:3 0.3 1:4 1:3 1:3 0.0在双孔/双渗模型情况下，需要对基岩和裂缝分别输入孔隙度值，必须先输入基岩的孔隙度值，然后输入裂缝系统的孔隙度，对于其他数据的输入过程也采用同样的顺序。 例: *por *matrix *ijk 1:10 1:10 1:3 0.3 1:4 1:3 1:3 0.0 *por *fracture *ijk 1:10 1:10 1:3 0.0 8 7:9 1:2 0.4 *mod 8 7:8 1 = 0.45这个例子同时也说明了 *mod

42、关键字的使用，它用来修改某些网格的性质，将孔隙度由0.40 改变为 0.45。在双孔隙系统中，只有当基岩和裂缝的孔隙度都为零时，才将它作为死结点处理，在一般情况下，基岩和裂缝的孔隙度可以一个设为零值，而另一个不为零。水区选项的使用_在 imex 中模拟水区的公式来自 carter 和 tracy 的“一种计算水侵量的改进方法”。这些公式是对 everdingen 和 hurst 公式的一种近似，请参考 van everdingen 和 w. hurst 的“对于油藏流动问题的拉普拉斯变换应用”和 l. p. dake 的“油藏工程基础”第九章。 建立水区模型时，可以使用下述关键字: a) *aquifer, b) *aqprop, c) *aqleak d) *aqfunc. 这些关键字放在油藏描述数据段中，它们完全是可选择的，如果使用这些选项，必