CMG组分模块GEM教程.docx

《CMG组分模块GEM教程.docx》由会员分享，可在线阅读，更多相关《CMG组分模块GEM教程.docx（198页珍藏版）》请在课桌文档上搜索。

1、IGEM概述在提市采收率项目，包括济气或济溶解剂中，藁过程可以混相或者非混 相，它取决于入流和油臧压略和温度。例如，富气驱、高气驱、C02驱和凝析 气臧的循环注气。其模拟需根热力学和流体流动的特殊处理。GEM就是一有效的多维状态方程组分模拟器全部相注气的重根机理，即油 的蒸发和膨胀、气凝析、粘度和界面张力降低，通过多次接触形成混相溶解剂 段塞。自适应隐式公式、全隐式自应隐式三种运行。在大多数情形下，只有很少 量的网格需根全隐式求解。大多数网格可以显式求解。自适应隐式方法在计算 中动态选择网格的隐式度，它对于井筒附近的锥进，或非常薄层的状油臧非常 用。使用自适应隐式方式，可以比隐式节约三分之一到

2、一半的时间，而时间不 步可以和全部隐式方法一样长。用户选择全隐式格块，GEM然后进一步动态择隐式格块。双孔隙度/双渗透率双孔模型将裂缝油藏的孔隙空间分为两部分：主孔隙和次孔隙。主孔隙（基岩）表示岩块基质中小的粒间孔隙，次孔隙 （裂缝）、结点、孔洞组成。双孔模型以一个油藏体积两种介质表示为特征。 具有较大存贮性的裂缝，是流体流动的主要通道，而基质则为存贮空间。在GEM中，可指定单的双孔模型。每个单元都分一基质和缝形状因子。在 这种情况下，基质和裂缝之间的交换项处于半稳态流动。除双孔模型以外，还可指定双渗模型，该模型假定相邻基岩块之一间存在 流体流动。在基质、基质之间质量交换不能忽略的情况下有用，

3、例如气油重力 驱占主要地位的情况。状态方程GEM使用Peng-Robinson或Soave-Redlich-Kwong状态方程预测油相和 气相的平衡组成和密度，应用Jossi和Thodos方程预测油和气的粘度。闪蒸计算拟牛顿逐次换法QNSS用于求解闪蒸计算的非线方程，以Gibbs能量分析为 基础的稳定性测试监测单相状态。复杂油藏油藏定义关键字用来描述油藏，网格可以是可变厚度一可变深度型，也可 以是角点类型，断层也是可以定义的，可支持笛卡尔坐标和径向网格。拟组分选择拟组分（流动组分）可有效模拟一次接触混相过程和其它单相流体系统。全隐式井井可以完善的方式求解，井底流压和完井段的格块变量以全隐式方式

4、求 解。如果完井段不止一个层，其井底流压以完全对耦的方法求职解，即考虑所 有完井段。这个可以解决典型层状油藏中多完井段的收敛问题。而且在GEM中，复杂 的井控制性能是千变万化的。可输入各种约束条件（最大井底或井口压力、最 小井底和井口压力、最大产量、最大气油比、最小气油比等）。当违反其中一 条约束条件时，用户可指定另一新的约束条件。气循环选择允许对组分进行剥 离，可附加补充气到循环气流中。矩阵求解法GEM使用AlMSOL,它是一优秀的线性求解程序，基于不完全高斯消去法。AIMSoL尤其针对自适隐式雅可比矩阵而开发。油藏初始比初使油藏条件可给定气油和油水界面深度。流体组成可随深度而改变，油 藏温

5、度也可随深度呈线性变化。水区模型水层可以另增加仅有水的边界单元，或使用Carter-Tracy不层模型模拟。 前者用于水体大小和位置已知，而附加的水体网格又相对小的情况。后者一般 用于不体较大或未知，希望近似计算水侵量，而附加水体网格不可行的情况。单位：Sl或矿场单位模拟结果文件GEM产生SRF格式，第二代模拟结果文件*IRF和*MRF文件。结果文件可用 于二维和三维可视化软件，也可用于动态曲线图输出。如果希望重启，需要其它几个存在的文件，同时产生另外三个文件，如 下：关键字输入部分数据组a）七个不同组b）数据组应遵循确定的顺序I/O控制油藏描述组分特性岩石一流体数据初始条件控制井数据和循环数

6、据怎样建立数据文件注释a) *TITLE1b) *TITLE2c) *TITLE3d) CASEID这四个关键字都必须在输入/输出控制部分。这四个关键字是可选的，并可以去掉。但是，它们对于辩别不同数据文件 很有用。至少应使用一行标题。全部标题和标识必须包含在单引号之内。TITLEl和CAS日D均在模拟结果文件中使用，该文件用来产生模拟图形。*TITLE1可以有40个字符，*TITLE2和*TITLE3每个可有80个字符。*CASEID 最多少个字符。也可使用两个关键字标识符即*插入注释，注释可以出现在数据文件任何 地方。例：*TITLE1Simulation Run #1-1989-Ol-23

7、*TITLE2Dual porosity problem using MINC option/*TITLE3This is a 121210Cattesian grid system*CASEID,Run*注释部分怎样重启什么是重启文件？重启文件是一个二进制文件,初始数据和主要变量按用户指明频率写入,写重 启文件是可选择的.为什么需要重启？有以下几点理由：a）做敏感性分析或历史模拟b）修改井定义c）在运行一大的长作业之前，做一短的模拟运行，看结果是否满意。d）随意运行后的运行中节约运行时间。例如，当完成一模拟运行且初步结 果盾起来不错，则需做预测运行。因为已经由前次运行创建一重启文件，则可以选

8、取一中间时间步重启运 行。模拟器不需从初始日期启运行，可以在选择的时间步继续运行。怎样做一重启？重启记录是可选的。如果需要做重启运行，则要在前次运行时产生。使用*WRST划*RESTART创建重启文件，它们必须位于输入/输出控制部 分。然而，当井工作制度改变时*WRST可出现在井数据部分，*WRST指示写重 启记录的频率。*RESTART表示当前值，否则，输入时间步数。例：RESTART 30当做重启时：a）不改变原始油藏数据，但在输入/输出控制部分加一RESTART关键字。b）增加最大时间步数，如果必要删除*MAXSTEPS0C）需要附加输入文件：输入重启文件和输入模拟结果文件。这些文件在前

9、次运行期间产生d）需要新的一组输出文件名。控制输出打印文件内容控制输出打印文件内容，使用：a） *WPRNb） *OUTPRN这些关键字可以出现在输入/输出控制部分，其参数也可在随后井数据部分 修改。*WPRN表示写网格块数据、井数据和数值方法控制数据的频率。如果不希望输出网格或井数据，则频率置零。例如：*WPRN*WELLO*WPRN*GRIDO如果将这两个关键字删除，其缺省选择为每一次井变化时打印。这种打印 控制会产生一个非常大的输出打印文件，会很快占满计算机的可用空间。*OTPRN限制打印什么样的油藏信息和井信息。详细的网格信息和流体性质数据可使用*OUTPRN *RES产生控制模拟结果

10、文件内容模拟结果文件（SRF）是模拟运行期间产生的数据文件。SRF文件由图形和 控制读入，产生绘图输出。为控制是SRF的内容，使用：b) *OTSRF这两个关键字一般首次出现在输入/输出控制部分。*SWRF表示写网格数据 和井数据的频率。如果不希望网格数据或井数据输出，则输出频率置零。网格系统描述描述网格系统，使用：a) *GRIDb) *KDIR描述网格XX,使用C)*DJd) *DJe) *DK描述油藏XX和倾角，使用f) *DEPTHg) *DIP 或 h) *DTOP在以上关键字中，只有*KDIR和*DIP完全可选，并可能以从数据集中删 除。上述列出的关键字必须在油藏描述部分。*GRI

11、D必须是这一部分的第一个关 键字。*GRID表示网格系统类型，有四种选择：直角坐标、变深度/变厚度、径向井筒和角点网格。其中每个均要求I, J, K 方向网格块数目：例如：* GRID *CART 10 10 6* GRID*VARI 10 10 6* GRID *RADIAL 10 1 15* GRID*RNER10 10 6第一个为标准直角人坐标IOXIoX6,第二个为变深度/变厚度网格IOXlOX6,第 三个为径向网格,最后一个为角点网格。KDlR定义Z轴的方向，。缺省值为UP,第一层为最底层。死结点定义死结点有两种方法a) *NULLb) *POR两种关键字都必须在油臧描述部分，在数据

12、集中必须有*POR,但*NULL可 选。当用*NULL关键字表示无效网格时，0表示无效网格，1表示有效网可知， 例如：*NLL*IJK1： 101： 101： 311: 41： 31： 30第二行覆盖第一行。*NULL是可选的，如果不出现，则假定全部网格均不 有效网格。NULL关键字覆盖PoR关键字。即使给一网格赋予非零孔隙度，而在*NULL 关键字中赋给零值，则该网格在流动计算中也不予考虑。上述例子可以下列形式表示：例如：*POR*IJK1： 101： 101： 30.31:41:31:30.0双孔隙度/双渗透率调用孔/双渗选择，可以使用(只能选一个)；a) *DALPORb) *DALPE

13、RM除裂缝介质标准网格信息外，还需输入：c) *DIFRACd) *DJFRACe) *DKFRACf) *TRASFER在组分性质部分，还可以指定：g) *DIFFSION孔隙度值需要对基质和裂缝输入例：*POR*MATRIX*IJK1： 101： 101： 310.31:41:31:30,0*POR*FRACTRE*IJK1:101:101:30.087:91:20.4*M0D87:81=0.45水层选择：水层模型如下：a) *AQIFERb) *AQPROPc) AQFUNC步长过小或运行时间过长问题收敛失败可能由于：a)内迭代收敛失败b) XX收敛失败导致时步截短c)物质平衡误差如果在

14、输出文件中发现迭代不收敛频繁出现，那么可以：1减小时间步长。将*DTMAX改小，或用*NORM*PRESS和*NORM *SAUR减 小第个时间步的改变量。2用关键子*1TERMAX增加迭代次数3使用*SDEGREE增加因子分解度。注意：这个方法增加了存贮需求。如果问题原因是最大改变量，则若出夙不频 繁，问题不大。如果时间步大量重复出现，那么，应采取以下措施：4检查岩石和PVT曲线是否出现非线性，该曲线应当平滑。5检查网格和其他性质6检查井的限制。一个有效的办法是对每一注入井总是指定一最大井底压 力，对每一生产井总是指定一最小井底压力。7如果不收敛由于超过最大牛顿迭代次数引起，使用关键字*NE

15、WTONCYC, 增加牛顿循环次数。如果由于最大改变量振荡引起，达到最大迭代次数，那 么，使蜚线性段平滑，或减小时间步长是最好的解决办法。8油藏是否在迭代之间出现在单相油和单相气之间浮动。如果出现浮动，设 置适当的*PHASEID 或*REFDEN9置油藏的某些区域或整个油藏为隐式。缺省开关符，*AIM*STAB检查网 格的显式隐式开关，仅当该网格是一隐式格块的相邻格块。如果某些油藏区域 变量发生剧烈变化，并且不是井的毗邻区域，或流体前沿运动太快，而在一个 时间步扫过两个格块，那么置这些区域为隐式a）当不使用垂向平衡初始化时。在某些情况下，这可能引起初始压力饱和 度发生大的改变，即使全部井关井

16、。当该情况发生时使用全隐式。b）当出现气顶时。如果出现较强的指进，置气顶的底层为隐式，至少在指 进出现的区域。c）在具有极端高渗透率的格块，压力的微小改变可能引起饱和度很大变 化。在这些区域，置格为全隐式。物质平衡误差可能由收敛误差超过设定的改变量而引起。在大多数情况 下，*NORM和Converge的缺省值足够。然而，当模拟裂缝油藏或使用裂缝代 表水平井时，推荐使用较小的值。对锥进问题，推荐使用小值。模拟一单相油藏有两种主要方法模拟气水两相的干气，单相气藏或一次接触混相溶解剂注AoA）使用饱和压力使用关键字*PSAT置饱和压力为一很小的值例如，-Io在这种情况下，所有油藏流体是欠饱和的，并跳

17、过每个可知块 的闪蒸计算。所有流体性质仍然使用状态方程计算。*PSAT不影响井流动计算。B）使用拟合组分选择1使用*PSEUD00 一般来说，置油藏油为一拟组分，每个注入流体作为不同 的单个拟组分。2用*PHAS曰D将单相标为驰，或勺一般建议用由、从而流体相对渗透率 依据油水相对渗透率曲线计算。该方法也不做网格闪蒸计算。另外，GEM还将单个组分流动方程分组为较 小数目的拟组分流动方程。这个可以大大减少每个网格坟程的数目。水平井水平井可作为一线源（注入井）或汇（生产井）。GEM使用AXlX等流动图计算体滞留量，同时也计算井筒磨擦压力降。注 意，当前方法并不能解决回流问题。垂赂平衡VERTICAL

18、控制垂向平衡计算。处理垂向平衡计算时，使用下列几种方法：A) *VERTICAL *0N同时使用下列关键字*DWOC,*REFDEPTH,*REFPRES,*XGLOBALand*SWOCB) VERTICAL *GASCAP 及*DWOC*DGOC*REFDEPTH*REFPRES*ZOIL*ZGAS*SWOCC) *VERTICAL*COMP 及*DWOC*REFDEPTH*REFPRES*SWOC*CDEPTH*ZDEPTH这些关键字必须位于初始条件部分。VERTICAL ON主要用于欠饱和油藏，而VERTlCAL *GASCAP可用于饱和油藏，也可用于欠饱和油藏。对于VERTICAL

19、*0N,用重力一毛管压力平衡计算全部格压力和水饱和度。 当网格包含油相和气相时，不做特殊处理。因此，在第一个时间步可能出夙剧 烈的流体变化。*VERTICAL*GASCAP是更一般的选择，假定气顶油区处于平衡状态，用穗力 一毛管压力平衡计算全部网格压力和油、气、水饱和度。闪蒸计算用于确油和 气相组成。油相和气相依饱和度混合，然后计算网格总的组成。因此，油藏中 总的组成可能随深度变化。*VERTIACL*COMP与用户指定的随度变化的组成一起，用于欠饱和油藏或者 和油藏的重力初始化。这时，要指定水油界面深度正气参考压力。用户必须输 入总的组成与深度表。用户还必须输入临界深度。在临界深度以下，单煌

20、相网 格为油网格，而在临界深度以一，单燃相网格为气相网格。如果不做垂赂平衡计算，输入*VERTICAL*OFF 及*PRES*SW*ZGLOBAL下述第一个例子表示对饱和油藏做垂计算的初始条件数据：*INITIALVERTICAL *GASCAP*REFDEPTH 900.0*REFPRES 3000.0*SWOC1250.0*DGOC788.0*ZOIL0.300.100.200.300.120.08*ZGAS0.780.190.030.00.00.0*SWOC1.0第二个例子表示不使用垂赂平衡计算的数据：INITIALVERTICAL *OFF*PRES *CON 3400.0*SW*CO

21、N 0.2*GLOBAL*CON0.300.100.200.300.120.08井的定义*井使用以下关键字定义，注意严格按照关键字排序*WELL 井号xx井类型定义PRODUCER 井号或者*1N正CToR井号或者*CYCLPROD 井号在井类型以下，输入：*PWELLBORE （如果需要计算井筒压力，需要井筒数据）或 *1WELLBoREINCOMP （如果定义了注入井时需要，紧接着IN正CTOR）OPERATE （至少定义一个工作限制）*ONITOR （监视限制任）定义井的位置和射孔，使用：*GEOMETRY （可选的，在完井卡之前）*PERF （必需）这些关键字必须位于井数据部分。怎样关

22、井？怎样重新打开井？在关井之前：（1）该井必须用以下方式完全定义：（a）作约束和任意监视约束（b）完井关键字在完全定义一口井以后，该井可以在使用TIME或DATE关键字的任一时间 关闭。在关井之后任一时间可打开一口闭的井。*WELLl,water injector,*WELL2,producer,*WELL3,Solvent injector,*全部井具有相同的几何因了*rad gefac wfrac skin*GEOMETRY *K 0.25 0.341.00.0两口源程序入井初始关井，只有生产井开井。这咱工作制度持续两年。*INJECTOR 1*INCOMP *WATEROPERATE *

23、MAX *STW 12000.0OPERATE *MAX *BHP 100.0*PERF*GEO 1*ifjf kfff1131.0*关1号井*SHUTIN 1*PRODUVER2OPERATE *MAX *STO 12000.0OPERATE *MIN*BHP 1000.0*MIONTOR*GOR 100.0*STOP*MONITOR*WCUT 0.*STOP*PERF*GE02*ifjf Kfff* 1MlECTOR 3* INC0MP *SOLVENT0.780.190.030,0 0.0 0.0* OPERAE *MAX *STG 1.2E7* OPERATE *MAX *BHP 10

24、0.0* PERF *GEO3* *ifjf kfff1131.0-* *3号井和1号井在相同网格射孔* *关3号井*SHUTIN3*TIME730.000* OPEN 1*打开注水井*TIME1095.0* *在这个方案中，注水和注溶解剂井在同一网格定义；在打开第二个以 前，关闭第一个。* SHUTIM 1*关注水井* 0pen3*开注溶解剂井* TIME 同 14600*SHUTIN3*关注溶解剂井*OPEN 1*打开注水井*TIME18250*SHUTIN1*关注水井*OPEN3*开注溶解剂井井筒模型说明井筒模型时，需要：PWELLBoRE用于生产井，IWELLBORE用于注入井，位于井

25、类型卡之后。当要求井口压力时需要使用井筒模型，用于：a）仅是一种信息；b）当XX压力作为一个约束条件时。例如:对2号注入井建立井筒模型*INJECT0R2*IWELLBORE*wepthlength rough whtemp bhtemp wr1100.1100.3490.0 33工作和监视约束OPERATE和MONITOR说明一给定井的约束条件。至少要求一个工作约束, 监视约束是可选的。每一口井引入一新的未知变量Pbh （井底流压），要求一约束方程确定该 变量。在工作和监视约束表中第一个工作约束条件为主工作约束。模拟器首先在 主工作约束下运行，同时监视其它约束条件。当监视约束之一违反，而且使

26、用 CONT关键字，则该约束变为工作约束。如果不止一个工作约束违反，那么使用第一个工作约束。生产井对一生产井，应当：a）在产量约束下生产（主工作约束）b）在最小井底流压约束下生产。如果生产井为产油井，则选一油产量约束。如果生产井为产气井，则选一 气产量约束。对一生产井使用的约束可以是一最小产量。例如：PRODUCER 1OPERATE *MAC *ST01200.0*NTOPERATE *MIN *BHP 1500.0*NT这个例子表明：a）使用油产量作为油井的主约束b）使用井底流压作为从约束如果有一条违反，运行继续，工作约束改变为刚刚被违反的约束。CoNT是 缺省值不需输入。注入井对一注入井

27、，应选:a）最大注入时作为主工作约束b）最大井底流压约束如果为注水井，选水产量作为约束。如果为注气井，则选气井，则选产 量。例如：*INJECT0R2*OPERATE *MAX *STWIOOO.O*STOPOPERATE *MAX *BHP 2250.0*STOP该例表示：a）注水井的注水量是主约束b）同时监视蟆底流压,作为从约束只要有一个违反,则模拟停止.监视约束监视约束格式与*0Perate格式相同。而且，如果同时有多个约束违反，则 采取约束表中第一个约束的动作。极力推荐监控生产井的GOR和含水率：这个可以避免作业运行中的某些问题。例如：PRODUCER 1OPERATE *MAX *0

28、IL 1200.0 *CONT*OPERATE *MIN *BHP 2500.0 *CONT*MONITOR *MAX *GOR 15000.0*STQP*MONITOR *MIN *WCT 60.0 *SHTIN如果在一时间步，GoR和含水率均违反，模拟停止。输入井指数井指数可由下述两种方法获得：a）从井模型的井几何因子。使用下述关键字：*GEOMETRY*PERF *GEOb）使用关键字直接指定：*PERF这些关键字位于数据部分。* GEMTRY指定内部计算井指数的必要参数。完井关键字PERF包括完井段的 位置。如果*GEO子关键字不在*PERF中，则应指定井指数。* PERF可用作水平井

29、、斜井和垂直井。例如：*WELL l,12-09-18-56,* PERF 1* *if jf kf wi setn112： 4 1.24 1定义井指数为124或例如：*WELL l,12-09-18-56,* *rad gefac wfrac skin* GEOMETRY *K 375 2488 1.00.0* *完井关键字必须紧接着几何关键字* PERF *GEO 1* if jf kf ff setn112： 411停止模拟运行可使用关键字*STOP中止模拟运行。另一方法是，在一井改变之后，由*DTWELL指定的第一个时间步长为0。例：*DATE 1998 09 08*DTWELL0.0

30、怎样建立井数据当首次输入井数据时，数据必须依据以下顺序：1 要法语*TIME 或*DATExx2定义一*DTWEL值，该数据是井定义之后，立刻使用的第一个时间步长。3使用*WELL标识所有XX每组井定义由如下部分组成：4a）使用producer和*1NJECTOR定义一口新井或主工作制度改变。b）如果选择*WHP作为约束之一，必须指明井筒模型要求。c）定义注入流体组d)定义该井的工作或监控约束5使用*GEOMETRY定义井位置、几何因子或井指数，使用*PERF定义完井 段。6在第1步到第5步后可使用*SHUTN7使用*0PEN重新找开一关闭井8第9步到第14步可在不同时间或日期对井进行修改，注

31、意依据使用什么 选择要求不同的关键字。9定义一口新井，做第1步到第6步，应在调整现有井参数之前。10使用*ALTER可改变任一口井的主工作约束，同时使用*TEME或*DATE下 述关键字中，只有一种关键字可与*TIME或*DATE卡一起使用。换句话说，在一 个*TIME或*DATE下不能同时使用两组*AIMSET,两组*AIMWELL等。11可以按需要调整隐式度/显式度(*AIMSET, *AIMWELL) o12可以按需要调整输入/输出控制13关键字*DTWELL, *DTMAX, *DTMIN也可以出现在随后的井变化卡中。14终止模拟运行总是用*STOP或零*DTWELL卡。假定IOXIO

32、X12直角坐标系统，井数据的构成如下：*井数据*RUN*DATE 1990 1 1*DTWELL0.025*WELL 1,OILPRODUCER,*PRODUVER1*吸有生产井的工作和监控约束OPERATE *MAX *STO 1000.00*OPERATE *MIN *BHP 3000.00MONITOR *MAC *WCT 0.40* * rad gefac wfrac skin* GEOMETRY *K 0.25 0.34 1.0 0.0* PERF *GEO 1* *ifjf kfff1131* TIME 10.0* *在10天时，主约束值由1000改为100* ALTER 1100

33、.0* *在50天时，主约束值再改为50.0*TIME50.0*ALTER 11000.0定义一注入井，而井对生产井再次修改主约束*TIME720*INJECTOR 1*INCOMP *WATEROPERATE *MAC *STW 12000.0* OPERATE *MAX *BHP 100.0* *rad geofacwfrac skin* GEOMETRY *K 0.30 0.34 1.00.0* PERF*GE0 1* *ifjf kfff7731* ALTER 2100.0* *在第10年停止* DATE 1990 12 31*STOP其它部分循环数据其它部分也有一些关键字可在模拟运行

34、期间改变，这些关键字必须位于 *DATE或*TIME卡之后，并且后一个值覆盖前一个值。在输入/输出控制部分出现的这些关键字可以再赋值：a) *WRSTb) *WPRNc) *ORTPRNd) *WSRFe) *OTSRF *GRID*OTSRF GRID可以改变，但C)UTSRF *WELL不能改变。数值方法控制部分，出现的关键字可以再赋值的有：a) *DTMAXb) *DTMIN关键字输入系统关键字系统介始关键字标识符关键字前可带星号* 或不带星号如：Por 或 PoR 或*POR两个星号表示注释关键字顺序所有关键字分成关键字组。不同组关键字必须严格按手册中顺序排列。若无特殊说明同一组内关键

35、字 可以改变顺序。*LIST, *NOLIST, *INCLUDE, *C0MMENT 和TRANSLATE 可以出 现在数据文件的任何地方。有些关键字可以出现在关键字组内，出可在循环数 据中。每一关键字的描述说明该关键字是可选的或必须的。有些关键字与其它 关键字一起使用时是可选的或必须的。如果可选的关键字在数据文件中不出 现，则使用其缺省值。串关键字若一关键字不是另一关键字或另一组关键字，不能与主关键字在同一行。每一行是一个关键字加上它的子关键字，除非另有说明。例如，水层选项需要3个关键字。*AQUIFER具有一系列关键字，并需要定 义水层的位置。该选择还需要水层性质（*AQPROP）和水层

36、影响函数（*AQFNC） o三个关键字是独立的实体，不能放在同一行。不正确的关键字表示：*AQUIFER *BOUNDARY *AQPROP 1000.4.002 230*AQFNC0.11行XX一行最多130个字符。界定符空格，逗号或新行可作为关键字，数字或字符串的分隔符。字符串字符串用单引号标出，例如535-48-W5内嵌的单引号不能使用。可内嵌 双引号。This ISTHE “rightWay更名关键字：对主关键字可使用TRANSLATE更名注释（可选）格式*注释文本空行（可选）分隔不同组关键字。检查错误扫描（可选）仅用于数据格式检查，不做运行。格式：*CHECKONLY该关键字必须放在

37、I/O控制部分。包含文件（可选）格式：*1NCLUDE，文件名，控制数据文件列表（可选）格式*LIST*NOLIST缺省值：*LIST改变注释符：（可选）格式：*COMMENTab缺省值：*COMMEN,命名用更名规则改变关键字（可选）格式：TRANSLATE，用户关键字“CMG关键字，一个关键字字可有几个别名：例如：*TRANSLATEzKX,zPERMI,TRANSLATEX-PermerbiIWPERM1*TRANSLATEzPERMX,zPERMI,网格性质数组输入数组读选择数组赋值由五部分组成，其中两个个任选。句法：网格数组（数组限制符）读选择数据（数组修改符）定义：网格数级要赋值的

38、性制，例哪：*POR数组限制符对网格的不同元素赋值（例如，基岩和裂缝），数组限制符可选。选择 是:MATRIX*FRACTURE*EQUALSI如果没有数组限制符，则假定为*MATRIX读选择：*CON*IVAR*JVAR*KVAR*ALL*IJK数据实际网格数组值数组修改符*MOD输入基岩网格性质关键字：*MATRIX输入裂缝网格性质关键字：*FRACTUREJ和K方赂数据与I方向数据相同关键字：*EQUALSI例1：单孔隙系统*PERMI *CON 100.0*PERMJ *EQALSI*PERMK*EQALSI *2例2：天然裂缝*PERMI*MATRIX *CON 100.0*PERM

39、I *FRACTURE *CON 100*PERMJ *EQUALSI*PERMK*EQUALSI *2常值数组*CONIJK标识的数组输入关键字：*IJKII： 12 JI： J2 KI： K2 值数组输入值随I方赂变化关键字：*IVAR 值(1),值(ni)数组输入值随J方向变化关键字L*JVAR值(1)，值(nj)数组输入值随K方向变化关键字：*KVAR 值，值(nk)值随全部网可知改变关键字：*ALL 值(1), 值(ni nj nk)数值I方向下标增加最快，j方向次之，k方赂最慢。修改数组数据格式：*MODil：i2jl:j2 kl：k2 ( + ) value(-)(*)(/)(=

40、)或 MOD(*)ValUe(-)()(/)(=)输入输出控制输入输出控制说明GEM可以生成两种不同格式输出显式文件，SRF和SR2,缺省输出文件仅有 SRF文件。使关键字*RESULTFILE可使GEM仅生成SR2或生成SRF和SR2文件。输入输出文件名格式：FILENAMES文件类型名称选择成FILENAME文件类型名称选择其中文件类型是下列之一或更多：*0UTPT*RESTARTIN*RESTARTOT*SRFIN*SRFOT*INDEX-OUT*MAIN-RESULTS-OUT*INDEX-IN*MAIN-RESULTS-IN*SR2ASCOUT对文件类型为*NDEX-OUT , *M

41、AIN-RESULTS-OT, *MAIN-RESULTS- IN,可以指定下列名称选择之一。文件名PROMPT*SCREEN对文件类型*NDEX-IN和RESTARTIN,可以指定下列名称选择之一：文件名*PROMPT格式：FILENAMES (OUTPUT)()(，文件名，)(PROMPT)(SCREEN)(*RESTARTIN)(文件名)(*PROMPT)()如果指定*1NDEXIN(*RESTARTOUT)()(文件名)(*PROMPT)(*SRFIN)()(文件名)(*PROMPT) (*NONE)(*SRFOT)()(文件名)(*PROMPT)(INDEX-OUT)()。文件名，)

42、(*PROMPT)(*MAIN-RESULTS-OUT)( )。文件名，)(*PROMPT)(*INDEX-IN) (1 文件名，)(*PROMPT)(*MAIN-RESLTS-IN)()(，文件名f)(*SR2ASC0UT)()（文件名F）（PROMPT）定义：FILENAMES文件名关键字* OUTPT模拟结果主输出文件restartin如果选restart,该关键字指明强制读入的重启动数据文件。RESTARToUT如果选WRST,为强制定入的重启动数据文件。SRFIN如果选 RESTART,为读入的SRF数据文件。如果没有该数据文件，使*NONESRFOUT如果已选SRF写选择WSRF,

43、表示要写入的模拟结果文件（SRF） INDEX-OUT如果指定RESUKTFILE*SR2,定入的SR2 ASCII下标文件。* MAIN-RESULTS-OUT定入的SR2主二进制结果文件。* 1NDEX-IN重启运行读入的SR2 ASCII下标文件* MAIN-RESLTS-IN重启运行读入的SR2主二进制结果文件。* SR2ASCOUTSR2主二进制结果的SACII文件。空串，表示使用缺省文件名Q文件名，文件名。PRoMPT提示用户从键般输入文件名。如果数据文件中不使用FILENAMES, 这是缺省值。* SCREEN该文件在屏幕上显示（标准输出设备）。*NONE表示不提供输入文件。当输

44、入文件丢失，或用户希望不考虑以前历 史数据时，该选择是有用的。缺省文件名仅当名称选择使用时，用缺省文件名缺省文件名具有形式yyy,其中XXXX为根，YYY为后缀，后缀依赖于 有关的文件。例如：9UT，为主输出文件，RST为重启，srf为SRF文件。根据的取法规则如下：（a）取掉包括，。分界符的后缀；（b）除最右边了符串和一，号以外，取 掉所有字符；（C）取掉80列以后字符，例如：DOS 路径，CMGDATATEST_3.DAT，具有根TEST_3，下列表总结了缺省文件名的根和后缀：文件根基二后缀*OTPUTOut数据文件名*OUTPTlrf*INDEX-OUT*RESTARTOT*OUTPTRstMrf*MAIN-RESULTS-*OTPUTOUT*INDEX-IN （无缺省值）*RESTARTIN*INDEX-IN Rst*MAIN-RESULTS*1NDEX-IN MrfOUT*SRFOT*OUTPUTSrf*MAIN*RESTARTINsrf使用这套缺省系统，用户对每一次运行，仅修改*0UTPUT和*1NDEX-IN或*RestartiN