CMG组分模块GEM教程

1GEM概述

在提市采收率项目，包括济气或济溶解剂中，蓁过程可以混相或者非混相，它取决于入流和油臧压略和温度。例如，富气驱、高气驱、CO2驱和凝析气臧的循环注气。其模拟需根热力学和流体流动的特殊处理。

GEM就是一有效的多维状态方程组分模拟器全部相注气的重根机理，即油的蒸发和膨胀、气凝析、粘度和界面张力降低，通过多次接触形成混相溶解剂段塞。 自适应隐式公式、全隐式自应隐式三种运行。在大多数情形下，只有很少量的网格需根全隐式求解。大多数网格可以显式求解。自适应隐式方法在计算中动态选择网格的隐式度，它对于井筒附近的锥进，或非常薄层的状油臧非常用。 使用自适应隐式方式，可以比隐式节约三分之一到一半的时间，而时间不步可以和全部隐式方法一样长。

用户选择全隐式格块，GEM然后进一步动态择隐式格块。 双孔隙度/双渗透率

双孔模型将裂缝油藏的孔隙空间分为两部分：主孔隙和次孔隙。主孔隙（基岩）表示岩块基质中小的粒间孔隙，次孔隙（裂缝）、结点、孔洞组成。双孔模型以一个油藏体积两种介质表示为特征。具有较大存贮性的裂缝，是流体流动的主要通道，而基质则为存贮空间。

在GEM中，可指定单的双孔模型。每个单元都分一基质和缝形状因子。在这种情况下，基质和裂缝之间的交换项处于半稳态流动。 除双孔模型以外，还可指定双渗模型，该模型假定相邻基岩块之一间存在流体流动。在基质、基质之间质量交换不能忽略的情况下有用，例如气油重力驱占主要地位的情况。 状态方程

GEM使用Peng—Robinson或Soave—Redlich—Kwong状态方程预测油相和气相的平衡组成和密度，应用Jossi和Thodos方程预测油和气的粘度。 闪蒸计算

拟牛顿逐次换法QNSS用于求解闪蒸计算的非线方程，以Gibbs能量分析为基础的稳定性测试监测单相状态。 复杂油藏

油藏定义关键字用来描述油藏，网格可以是可变厚度一可变深度型，也可以是角点类型，断层也是可以定义的，可支持笛卡尔坐标和径向网格。 拟组分选择

拟组分（流动组分）可有效模拟一次接触混相过程和其它单相流体系统。 全隐式井

井可以完善的方式求解，井底流压和完井段的格块变量以全隐式方式求解。如果完井段不止一个层，其井底流压以完全对耦的方法求职解，即考虑所有完井段。这个可以解决典型层状油藏中多完井段的收敛问题。而且在GEM中，复杂的井控制性能是千变万化的。可输入各种约束条件（最大井底或井口压力、最小井底和井口压力、最大产量、最大气油比、最小气油比等）。当违反其中一条约束条件时，用户可指定另一新的约束条件。气循环选择允许对组分进行剥离，可附加补充气到循环气流中。 矩阵求解法

GEM使用AIMSOL，它是一优秀的线性求解程序，基于不完全高斯消去法。AIMSOL尤其针对自适隐式雅可比矩阵而开发。

油藏初始比

初使油藏条件可给定气油和油水界面深度。流体组成可随深度而改变，油藏温度也可随深度呈线性变化。 水区模型

水层可以另增加仅有水的边界单元，或使用Carter—Tracy不层模型模拟。前者用于水体大小和位置已知，而附加的水体网格又相对小的情况。后者一般用于不体较大或未知，希望近似计算水侵量，而附加水体网格不可行的情况。 单位：

SI或矿场单位 模拟结果文件

GEM产生SRF格式，第二代模拟结果文件*IRF和*MRF文件。结果文件可用于二维和三维可视化软件，也可用于动态曲线图输出。

如果希望重启，需要其它几个存在的文件，同时产生另外三个文件，如下： 关键字输入部分数据组 a）七个不同组

b）数据组应遵循确定的顺序 I/O控制 油藏描述 组分特性

岩石一流体数据 初始条件控制

井数据和循环数据 怎样建立数据文件注释 a）*TITLE1 b）*TITLE2 c）*TITLE3 d）CASEID

这四个关键字都必须在输入/输出控制部分。 这四个关键字是可选的，并可以去掉。但是，它们对于辩别不同数据文件很有用。至少应使用一行标题。全部标题和标识必须包含在单引号之内。

TITLE1和CASEID均在模拟结果文件中使用，该文件用来产生模拟图形。*TITLE1可以有40个字符，*TITLE2和*TITLE3每个可有80个字符。*CASEID最多少个字符。

也可使用两个关键字标识符即‘**’插入注释，注释可以出现在数据文件任何地方。

例：*TITLE1

‘Simulation Run #1—1989—01—23’ *TITLE2

‘Dual porosity problem using MINC option’ *TITLE3

‘This is a 12×12×10Cattesian grid system' *CASEID 'Runl' **注释部分 怎样重启

什么是重启文件?

重启文件是一个二进制文件,初始数据和主要变量按用户指明频率写入,写重启文件是可选择的. 为什么需要重启？ 有以下几点理由：

a）做敏感性分析或历史模拟 b）修改井定义

c）在运行一大的长作业之前，做一短的模拟运行，看结果是否满意。

d）随意运行后的运行中节约运行时间。例如，当完成一模拟运行且初步结果盾起来不错，则需做预测运行。

因为已经由前次运行创建一重启文件，则可以选取一中间时间步重启运行。模拟器不需从初始日期启运行，可以在选择的时间步继续运行。 怎样做一重启？

重启记录是可选的。如果需要做重启运行，则要在前次运行时产生。

使用*WRST划*RESTART创建重启文件，它们必须位于输入/输出控制部分。然而，当井工作制度改变时*WRST可出现在井数据部分，*WRST指示写重启记录的频率。*RESTART表示当前值，否则，输入时间步数。 例：*RESTART 30 当做重启时：

a）不改变原始油藏数据，但在输入/输出控制部分加一*RESTART关键字。 b）增加最大时间步数，如果必要删除*MAXSTEPS。

c）需要附加输入文件：输入重启文件和输入模拟结果文件。这些文件在前次运行期间产生

d）需要新的一组输出文件名。 控制输出打印文件内容

控制输出打印文件内容，使用： a）*WPRN b）*OUTPRN

这些关键字可以出现在输入/输出控制部分，其参数也可在随后井数据部分修改。 *WPRN表示写网格块数据、井数据和数值方法控制数据的频率。 如果不希望输出网格或井数据，则频率置零。 例如：

*WPRN *WELL 0 *WPRN *GRID 0

如果将这两个关键字删除，其缺省选择为每一次井变化时打印。这种打印控制会产生一个非常大的输出打印文件，会很快占满计算机的可用空间。 *OUTPRN限制打印什么样的油藏信息和井信息。

详细的网格信息和流体性质数据可使用*OUTPRN *RES 产生 控制模拟结果文件内容

模拟结果文件（SRF）是模拟运行期间产生的数据文件。SRF文件由图形和控制读入，产生绘图输出。为控制是SRF的内容，使用： a）*WSRF b）*OUTSRF

这两个关键字一般首次出现在输入/输出控制部分。

*SWRF表示写网格数据和井数据的频率。

如果不希望网格数据或井数据输出，则输出频率置零。 网格系统描述

描述网格系统，使用： a）*GRID b）*KDIR

描述网格长度，使用 c）*DJ d）*DJ e）*DK

描述油藏深度和倾角，使用 f）*DEPTH g）*DIP 或

h）*DTOP

在以上关键字中，只有*KDIR和*DIP完全可选，并可能以从数据集中删除。 上述列出的关键字必须在油藏描述部分。*GRID必须是这一部分的第一个关键字。

*GRID表示网格系统类型，有四种选择：直角坐标、变深度/变厚度、径向井筒和角点网格。其中每个均要求I，J，K方向网格块数目： 例如：

*GRID *CART 10 10 6 *GRID*VARI 10 10 6 *GRID *RADIAL 10 1 15 *GRID *CORNER 10 10 6

第一个为标准直角人坐标10×10×6,第二个为变深度/变厚度网格10×10×6,第三个为径向网格,最后一个为角点网格。 KDIR定义Z轴的方向，。缺省值为UP，第一层为最底层。 死结点

定义死结点有两种方法 a）*NULL b）*POR

两种关键字都必须在油藏描述部分，在数据集中必须有*POR，但*NULL可选。 当用*NULL关键字表示无效网格时，0表示无效网格，1表示有效网可知，例如： *NULL *IJK

1：10 1：10 1：3 1 1：4 1：3 1：3 0

第二行覆盖第一行。*NULL是可选的，如果不出现，则假定全部网格均不有效网格。

NULL关键字覆盖POR关键字。即使给一网格赋予非零孔隙度，而在 *NULL关键字中赋给零值，则该网格在流动计算中也不予考虑。 上述例子可以下列形式表示： 例如：*POR *IJK

1：10 1：10 1：3 0.3

1:4 1:3 1:3 0.0 双孔隙度/双渗透率

调用孔/双渗选择，可以使用（只能选一个）； a）*DUALPOR b）*DUALPERM

除裂缝介质标准网格信息外，还需输入： c）*DIFRAC d）*DJFRAC e）*DKFRAC f）*TRASFER

在组分性质部分，还可以指定： g）*DIFFUSION

孔隙度值需要对基质和裂缝输入 例：*POR *MATRIX *IJK

1：10 1：10 1：3 10.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 水层选择：

水层模型如下： a）*AQUIFER b）*AQPROP c）AQFUNC

步长过小或运行时间过长问题 收敛失败可能由于： a）内迭代收敛失败

b）牛顿收敛失败导致时步截短 c）物质平衡误差

如果在输出文件中发现“迭代不收敛”频繁出现，那么可以：

1减小时间步长。将*DTMAX改小，或用*NORM*PRESS和*NORM *SAUR减小第个时间步的改变量。

2用关键子 *ITERMAX增加迭代次数

3使用*SDEGREE增加因子分解度。注意：这个方法增加了存贮需求。如果问题原因是最大改变量，则若出夙不频繁，问题不大。如果时间步大量重复出现，那么，应采取以下措施：

4检查岩石和PVT曲线是否出现非线性，该曲线应当平滑。 5检查网格和其他性质

6检查井的限制。一个有效的办法是对每一注入井总是指定一最大井底压力，对每一生产井总是指定一最小井底压力。

7如果不收敛由于超过最大牛顿迭代次数引起，使用关键字*NEWTONCYC，增加牛顿循环次数。如果由于最大改变量振荡引起，达到最大迭代次数，那么，使

蜚线性段平滑，或减小时间步长是最好的解决办法。

8油藏是否在迭代之间出现在单相油和单相气之间浮动。如果出现浮动，设置适当的*PHASEID或*REFDEN。

9置油藏的某些区域或整个油藏为隐式。缺省开关符，*AIM *STAB检查网格的显式隐式开关，仅当该网格是一隐式格块的相邻格块。如果某些油藏区域变量发生剧烈变化，并且不是井的毗邻区域，或流体前沿运动太快，而在一个时间步扫过两个格块，那么置这些区域为隐式

a）当不使用垂向平衡初始化时。在某些情况下，这可能引起初始压力饱和度发生大的改变，即使全部井关井。当该情况发生时使用全隐式。

b）当出现气顶时。如果出现较强的指进，置气顶的底层为隐式，至少在指进出现的区域。

c）在具有极端高渗透率的格块，压力的微小改变可能引起饱和度很大变化。在这些区域，置格为全隐式。

物质平衡误差可能由收敛误差超过设定的改变量而引起。在大多数情况下，*NORM和*Converge的缺省值足够。然而，当模拟裂缝油藏或使用裂缝代表水平井时，推荐使用较小的值。对锥进问题，推荐使用小值。 模拟一单相油藏

有两种主要方法模拟气水两相的干气，单相气藏或一次接触混相溶解剂注入。 A） 使用饱和压力

使用关键字*PSAT置饱和压力为一很小的值

例如，—1。在这种情况下，所有油藏流体是欠饱和的，并跳过每个可知块的闪蒸计算。所有流体性质仍然使用状态方程计算。*PSAT不影响井流动计算。 B） 使用拟合组分选择

1使用*PSEUDO。一般来说，置油藏油为一拟组分，每个注入流体作为不同的单个拟组分。

2用*PHASEID将单相标为‘油’或‘气’。一般建议用‘油’，从而流体相对渗透率依据油水相对渗透率曲线计算。

该方法也不做网格闪蒸计算。另外，GEM还将单个组分流动方程分组为较小数目的拟组分流动方程。这个可以大大减少每个网格坟程的数目。 水平井

水平井可作为一线源（注入井）或汇（生产井）。

GEM使用AXIX等流动图计算体滞留量，同时也计算井筒磨擦压力降。注意，当前方法并不能解决回流问题。 垂赂平衡

*VERTICAL控制垂向平衡计算。

处理垂向平衡计算时，使用下列几种方法： A）*VERTICAL *ON 同时使用下列关键字 *DWOC，

*REFDEPTH， *REFPRES， *XGLOBAL and *SWOC

B）*VERTICAL *GASCAP及

*DWOC *DGOC

*REFDEPTH *REFPRES *ZOIL *ZGAS *SWOC

C）*VERTICAL *COMP及 *DWOC *REFDEPTH *REFPRES *SWOC *CDEPTH *ZDEPTH

这些关键字必须位于初始条件部分。

*VERTICAL ON主要用于欠饱和油藏，而VERTICAL * GASCAP可用于饱和油藏，也可用于欠饱和油藏。 对于*VERTICAL *ON ，用重力一毛管压力平衡计算全部格压力和水饱和度。当网格包含油相和气相时，不做特殊处理。因此，在第一个时间步可能出夙剧烈的流体变化。

*VERTICAL *GASCAP是更一般的选择，假定气顶油区处于平衡状态，用穗 力一毛管压力平衡计算全部网格压力和油、气、水饱和度。闪蒸计算用于确油和气相组成。油相和气相依饱和度混合，然后计算网格总的组成。因此，油藏中总的组成可能随深度变化。

*VERTIACL *COMP与用户指定的随度变化的组成一起，用于欠饱和油藏或者和油藏的重力初始化。这时，要指定水油界面深度正气参考压力。用户必须输入总的组成与深度表。用户还必须输入临界深度。在临界深度以下，单烃相网格为油网格，而在临界深度以一，单烃相网格为气相网格。 如果不做垂赂平衡计算，输入 *VERTICAL *OFF 及 *PRES *SW

*ZGLOBAL

下述第一个例子表示对饱和油藏做垂计算的初始条件数据： **INITIAL

*VERTICAL *GASCAP *REFDEPTH 900.0 *REFPRES 3000.0 *SWOC 1250.0 *DGOC 788.0

*ZOIL 0.30 0.10 0.20 0.30 0.12 0.08 *ZGAS 0.78 0.19 0.03 0.0 0.0 0.0 *SWOC 1.0

第二个例子表示不使用垂赂平衡计算的数据：

*INITIAL

*VERTICAL *OFF *PRES *CON 3400.0 *SW *CON 0.2 *GLOBAL *CON

0.30 0.10 0.20 0.30 0.12 0.08 井的定义

*井使用以下关键字定义，注意严格按照关键字排序 *WELL 井号 ‘井名’ 井类型定义

*PRODUCER 井号 或者

*INJECTOR 井号 或者

*CYCLPROD 井号 在井类型以下，输入：

*PWELLBORE（如果需要计算井筒压力，需要井筒数据） 或*IWELLBORE

INCOMP（如果定义了注入井时需要，紧接着INJECTOR） *OPERATE（至少定义一个工作限制） *ONITOR（监视限制任） 定义井的位置和射孔，使用：

*GEOMETRY（可选的，在完井卡之前） *PERF（必需）

这些关键字必须位于井数据部分。 怎样关井？怎样重新打开井？ 在关井之前：

（1） 该井必须用以下方式完全定义： （a）作约束和任意监视约束 （b）完井关键字

在完全定义一口井以后，该井可以在使用TIME或DATE关键字的任一时间关闭。 在关井之后任一时间可打开一口闭的井。 *WELL1’water injector’

*WELL2’ producer’ *WELL3’Solvent injector’ **全部井具有相同的几何因了 ** **

** rad gefac wfrac skin

*GEOMETRY *K 0.25 0.34 1.0 0.0

两口源程序入井初始关井，只有生产井开井。 这咱工作制度持续两年。 *INJECTOR 1

*INCOMP *WATER

*OPERATE *MAX *STW 12000.0 *OPERATE *MAX *BHP 10000.0 *PERF *GEO 1 **if jf kf ff 1 1 3 1.0 **关1号井 *SHUTIN 1 *PRODUVER 2

*OPERATE *MAX *STO 12000.0 *OPERATE *MIN *BHP 1000.0

*MIONTOR *GOR 10000.0 *STOP *MONITOR *WCUT 0.8330002 *STOP *PERF *GEO 2 **if jf Kf ff 7 7 1 1.0

*INJECTOR 3

*INCOMP *SOLVENT

0.78 0.19 0.03 0,0 0.0 0.0 *OPERAE *MAX *STG 1.2E7 *OPERATE *MAX *BHP 10000.0 *PERF *GEO 3 **if jf kf ff ]1 1 3 1.0-

**3号井和1号井在相同网格射孔 **关3号井 *SHUTIN 3 *TIME730.000

*OPEN 1 **打开注水井 *TIME 1095.0

**在这个方案中，注水和注溶解剂井在同一网格定义；在打开第二个以前，关闭第一个。

*SHUTIM 1 **关注水井 *OPEN 3 **开注溶解剂井 *TIME 同1460 0

*SHUTIN 3 **关注溶解剂井 *OPEN 1 **打开注水井 *TIME 1825 0

*SHUTIN 1 **关注水井 *OPEN 3 **开注溶解剂井 井筒模型

说明井筒模型时，需要：

PWELLBORE用于生产井，IWELLBORE用于注入井，位于井类型卡之后。 当要求井口压力时需要使用井筒模型，用于： a）仅是一种信息；

b）当井口压力作为一个约束条件时。 例如：对2号注入井建立井筒模型 *INJECTOR 2 *IWELLBORE

**wepth length rough whtemp bhtemp wr 1100. 1100. 1365 200.3 490.0 33 工作和监视约束

OPERATE和MONITOR说明一给定井的约束条件。至少要求一个工作约束，监视约束是可选的。

每一口井引入一新的未知变量pbh（井底流压），要求一约束方程确定该变量。 在工作和监视约束表中第一个工作约束条件为主工作约束。模拟器首先在主工作约束下运行，同时监视其它约束条件。当监视约束之一违反，而且使用CONT关键字，则该约束变为工作约束。

如果不止一个工作约束违反，那么使用第一个工作约束。 生产井

对一生产井，应当：

a）在产量约束下生产（主工作约束） b）在最小井底流压约束下生产。

如果生产井为产油井，则选一油产量约束。如果生产井为产气井，则选一气产量约束。

对一生产井使用的约束可以是一最小产量。 例如：

*PRODUCER 1

*OPERATE *MAC *STO 120000.0 *CONT *OPERATE *MIN *BHP 1500.0 *CONT 这个例子表明：

a）使用油产量作为油井的主约束 b）使用井底流压作为从约束

如果有一条违反，运行继续，工作约束改变为刚刚被违反的约束。CONT是缺省值不需输入。 注入井

对一注入井，应选：

a）最大注入时作为主工作约束 b）最大井底流压约束

如果为注水井，选水产量作为约束。如果为注气井，则选气井，则选产量。 例如：

*INJECTOR 2

*OPERATE *MAX *STW 1000.0 *STOP *OPERATE *MAX *BHP 2250.0 *STOP 该例表示：

a）注水井的注水量是主约束

b）同时监视蟛底流压,作为从约束 只要有一个违反,则模拟停止. 监视约束

监视约束格式与*Operate格式相同。而且，如果同时有多个约束违反，则采取约束表中第一个约束的动作。

极力推荐监控生产井的GOR和含水率：这个可以避免作业运行中的某些问题。 例如：

*PRODUCER 1

*OPERATE *MAX *OIL 1200.0 *CONT *OPERATE *MIN *BHP 2500.0 *CONT *MONITOR *MAX *GOR 15000.0 *STOP *MONITOR *MIN *WCUT 60.0 *SHUTIN

如果在一时间步，GOR和含水率均违反，模拟停止。 输入井指数

井指数可由下述两种方法获得： a） 从井模型的井几何因子。使用下述关键字： *GEOMETRY *PERF *GEO

b）使用关键字直接指定： *PERF

这些关键字位于数据部分。

*GEMTRY指定内部计算井指数的必要参数。完井关键字PERF包括完井段的位置。

如果*GEO子关键字不在*PERF中，则应指定井指数。 *PERF可用作水平井、斜井和垂直井。 例如：

*WELL 1’12-09-18-56’ *PERF 1

**if jf kf wi setn 1 1 2：4 1.24 1 定义井指数为1.24 或 例如：

*WELL 1’12-09-18-56’ **rad gefac wfrac skin

*GEOMETRY *K 375 2488 1.0 0.0 **完井关键字必须紧接着几何关键字 *PERF *GEO 1 ** if jf kf ff setn 1 1 2：4 1 1 停止模拟运行

可使用关键字*STOP中止模拟运行。

另一方法是，在一井改变之后，由*DTWELL指定的第一个时间步长为0。 例：*DATE 1998 09 08 *DTWELL 0.0 怎样建立井数据

当首次输入井数据时，数据必须依据以下顺序：

1要法语*TIME或*DATE卡

2定义一*DTWEL值，该数据是井定义之后，立刻使用的第一个时间步长。 3使用*WELL标识所有新井 每组井定义由如下部分组成：

4a）使用*producer和*INJECTOR定义一口新井或主工作制度改变。 b）如果选择*WHP作为约束之一，必须指明井筒模型要求。 c）定义注入流体组

d）定义该井的工作或监控约束

5使用*GEOMETRY定义井位置、几何因子或井指数，使用*PERF定义完井段。 6在第1步到第5步后可使用*SHUTN 7使用*OPEN重新找开一关闭井

8第9步到第14步可在不同时间或日期对井进行修改，注意依据使用什么选择要求不同的关键字。

9定义一口新井，做第1步到第6步，应在调整现有井参数之前。

10使用*ALTER可改变任一口井的主工作约束，同时使用*TEME或*DATE 下述关键字中，只有一种关键字可与*TIME或*DATE卡一起使用。换句话说，在一个*TIME或*DATE下不能同时使用两组*AIMSET，两组*AIMWELL等。 11可以按需要调整隐式度/显式度（*AIMSET，*AIMWELL）。 12可以按需要调整输入/输出控制

13关键字*DTWELL，*DTMAX，*DTMIN也可以出现在随后的井变化卡中。 14终止模拟运行总是用*STOP或零*DTWELL卡。 假定10×10×12直角坐标系统，井数据的构成如下： **井数据 *RUN

*DATE 1990 1 1 *DTWELL 0. 025

*WELL 1’OIL PRODUCER’ *PRODUVER 1

**吸有生产井的工作和监控约束 *OPERATE *MAX *STO 1000.00 *OPERATE *MIN *BHP 3000.00 *MONITOR *MAC *WCUT 0.40 ** rad gefac wfrac skin

*GEOMETRY *K 0.25 0.34 1.0 0.0 *PERF *GEO 1 **if jf kf ff 1 1 3 1

*TIME 10.0

**在10天时，主约束值由1000改为100 *ALTER 1 100.0

**在50天时，主约束值再改为50.0 *TIME 50.0

*ALTER 1 1000.0

**定义一注入井，而井对生产井再次修改主约束 *TIME 720 *INJECTOR 1

*INCOMP *WATER

*OPERATE *MAC *STW 12000.0 *OPERATE *MAX *BHP 10000.0 **rad geofac wfrac skin

*GEOMETRY *K 0.30 0.34 1.0 0.0 *PERF *GEO 1 **if jf kf ff 7 7 3 1 *ALTER 2 100.0

**在第10年停止 *DATE 1990 12 31 *STOP

其它部分循环数据

其它部分也有一些关键字可在模拟运行期间改变，这些关键字必须位于*DATE或*TIME卡之后，并且后一个值覆盖前一个值。

在输入/输出控制部分出现的这些关键字可以再赋值： a）*WRST b）*WPRN c）*ORTPRN d）*WSRF

e）*OUTSRF *GRID

*OUTSRF GRID可以改变，但OUTSRF *WELL不能改变。 数值方法控制部分，出现的关键字可以再赋值的有： a）*DTMAX b）*DTMIN 关键字输入系统 关键字系统介始 关键字标识符

关键字前可带星号“**”或不带星号 如：por或POR或*POR 两个星号表示注释 关键字顺序

所有关键字分成关键字组。

不同组关键字必须严格按手册中顺序排列。若无特殊说明同一组内关键字内关键字可以改变顺序。*LIST，*NOLIST，*INCLUDE，*COMMENT和*TRANSLATE可以出现在数据文件的任何地方。有些关键字可以出现在关键字组内，出可在循环数据中。每一关键字的描述说明该关键字是可选的或必须的。有些关键字与其它关键字一起使用时是可选的或必须的。如果可选的关键字在数据文件中不出

现，则使用其缺省值。 串关键字

若一关键字不是另一关键字或另一组关键字，不能与主关键字在同一行。 每一行是一个关键字加上它的子关键字，除非另有说明。

例如，水层选项需要3个关键字。*AQUIFER具有一系列关键字，并需要定义水层的位置。该选择还需要水层性质（*AQPROP）和水层影响函数（*AQFUNC）。三个关键字是独立的实体，不能放在同一行。 不正确的关键字表示：

*AQUIFER *BOUNDARY *AQPROP 1000.4.002 230 *AQFUNC 0.11 … …. 行长度

一行最多130个字符。 界定符

空格，逗号或新行可作为关键字，数字或字符串的分隔符。 字符串

字符串用单引号标出，例如‘5-35-48-W5’，内嵌的单引号不能使用。可内嵌双引号。

‘This Is THE “right“ Way’ 更名关键字：

对主关键字可使用 TRANSLATE更名 注释（可选） 格式 **注释文本 空行（可选）

分隔不同组关键字。 检查错误扫描（可选）

仅用于数据格式检查，不做运行。 格式：*CHECKONLY

该关键字必须放在I/O控制部分。 包含文件（可选）

格式：*INCLUDE ‘文件名’ 控制数据文件列表（可选） 格式*LIST *NOLIST

缺省值：*LIST 改变注释符：（可选）

格式：*COMMENT‘ab’ 缺省值：*COMMENT‘’

命名用更名规则改变关键字（可选） 格式：*TRANSLATE ‘用户关键字’‘CMG关键字’ 一个关键字字可有几个别名： 例如：

*TRANSLATE ‘KX’‘PERMI’

*TRANSLATE‘X—permerbility’‘PERMI’ *TRANSLATE ‘PERMX’‘PERMI’ 网格性质数组输入 数组读选择

数组赋值由五部分组成，其中两个个任选。 句法：

网格数组（数组限制符）读选择 数据 （数组修改符） 定义： 网格数级

要赋值的性制，例哪：*POR 数组限制符

对网格网格的不同元素赋值（例如，基岩和裂缝），数组限制符可选。选择是： *MATRIX *FRACTURE *EQUALSI

如果没有数组限制符，则假定为*MATRIX 读选择： *CON *IVAR *JVAR *KVAR *ALL *IJK 数据

实际网格数组值 数组修改符 *MOD

输入基岩网格性质 关键字：*MATRIX 输入裂缝网格性质 关键字：*FRACTURE

J和K方赂数据与I方向数据相同 关键字：*EQUALSI 例1：单孔隙系统 *PERMI *CON 100.0 *PERMJ *EQUALSI *PERMK *EQUALSI *2 例2：天然裂缝

*PERMI *MATRIX *CON 100.0 *PERMI *FRACTURE *CON 10000 *PERMJ *EQUALSI *PERMK *EQUALSI *2 常值数组

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