CMG组分模块GEM教程

1GEM概述在提市釆收率项目，包括济气或济溶解剂中，蕖过程可以混相或者非混相，它取决于入流和油臧压略和温度。例如，富气驱、高气驱、CO2驱和凝析气臧的循环注气。其模拟需根热力学和流体流动的特殊处理。GEM就是一有效的多维状态方程组分模拟器全部相注气的重根机理，即油的蒸发和膨胀、气凝析、粘度和界面张力降低，通过多次接触形成混相溶解剂段塞。自适应隐式公式、全隐式自应隐式三种运行。在大多数情形下，只有很少量的网格需根全隐式求解。大多数网格可以显式求解。自适应隐式方法在计算中动态选择网格的隐式度，它对于井筒附近的锥进，或非常薄层的状油臧非常用。使用自适应隐式方式，可以比隐式节约三分之一到一半的时间，而时间不步可以和全部隐式方法一样长。用户选择全隐式格块，GEM然后进一步动态择隐式格块。双孔隙度/双渗透率双孔模型将裂缝油藏的孔隙空间分为两部分：主孔隙和次孔隙。主孔隙（基岩）表示岩块基质中小的粒间孔隙，次孔隙（裂缝）、结点、孔洞组成。双孔模型以一个油藏体积两种介质表示为特征。具有较大存贮性的裂缝，是流体流动的主要通道，而基质则为存贮空间。在GEM中，可指定单的双孔模型。每个单元都分一基质和缝形状因子。在这种情况下，基质和裂缝之间的交换项处于半稳态流动。除双孔模型以外，还可指定双渗模型，该模型假定相邻基岩块之一间存在流体流动。在基质、基质之间质量交换不能忽略的情况下有用，例如气油重力驱占主要地位的情况。状态方程GEM使用Peng—Robmson或Soave—Redlich—Kwong状态方程预测油相和气相的平衡组成和密度，应用Jossi和Thodos方程预测油和气的粘度。闪蒸计算拟牛顿逐次换法QNSS用于求解闪蒸计算的非线方程，以Gibbs能量分析为基础的稳定性测试监测单相状态。复杂油藏油藏定义关键字用來描述油藏，网格可以是可变厚度一可变深度型，也可以是角点类型，断层也是可以定义的，可支持笛卡尔坐标和径向网格。拟组分选择拟组分（流动组分）可有效模拟一次接触混相过程和其它单相流体系统。全隐式井井可以完善的方式求解，井底流压和完井段的格块变量以全隐式方式求解。如果完井段不止一个层，其井底流压以完全对耦的方法求职解，即考虑所有完井段。这个可以解决典型层状油藏中多完井段的收敛问题。而且在GEM中，复杂的井控制性能是「变7J化的。可输入各种约束条件（最大井底或井口压力、最小井底和井口压力、最大产量、最大气油比、最小气油比等）。当违反其中一条约束条件时，用户可指定另一新的约束条件。气循环选择允许对组分进行剥离，可附加补充气到循环气流中。矩阵求解法GEM使用AIMSOL,它是一优秀的线性求解程序，基于不完全高斯消去法。AIMSOL尤其针对自适隐式雅可比矩阵而开发。油藏初始比初使油藏条件可给定气油和油水界面深度。流体组成可随深度而改变，油藏温度也可随深度呈线性变化。水区模型水层可以另增加仅有水的边界单元，或使用Cartel—Tiacy不层模型模拟。前者用于水体大小和位置己知，而附加的水体网格乂相对小的情况。后者一般用于不体较大或未知，希望近似计算水侵量，而附加水体网格不可行的情况。单位：SI或矿场单位模拟结果文件GEM产生SRF格式，第二代模拟结果文件*吨和*“町文件。结果文件可用于二维和三维可视化软件，也可用于动态曲线图输出。如果希望重启，需要其它几个存在的文件，同时产生另外三个文件，如下：关键字输入部分数据组七个不同组数据组应遵循确定的顺序I/O控制油藏描述组分特性岩石一流体数据初始条件控制井数据和循环数据怎样建立数据文件注释*TITLE1*TITLE2*TITLE3CASEID这四个关键字都必须在输入/输出控制部分。这四个关键字是可选的，并可以去掉。但是，它们对于辩别不同数据文件很有用。至少应使用一行标题。全部标题和标识必须包含在单引号之内。TITLE1和CASEID均在模拟结果文件中使用，该文件用來产生模拟图形。*TITLE1可以有40个字符，*TITLE2和*TITLE3每个可有80个字符。*CASELD最多少个字符。也可使用两个关键字标识符即’**'插入注释，注释可以出现在数据文件任何地方。例：*TITLE1'SimulationRun#1—1989—01—23'*TITLE2'DualporosityproblemusingMINCoption'*TITLE3'Thisisa12X12XlOCattesiangridsystem'*CASEID'RunT**注释部分怎样重启什么是重启文件？重启文件是一个二进制文件，初始数据和主要变量按用户指明频率写入，写重启文件是可选择的.为什么需要重启？有以下几点理由：a） 做敏感性分析或历史模拟b） 修改井定义c） 在运行一大的长作业之前，做一短的模拟运行，看结果是否满意。d） 随意运行后的运行中节约运行时间。例如，当完成一模拟运行且初步结果盾起来不错，则需做预测运行。因为己经由前次运行创建一重启文件，则可以选取一中间时间步重启运行。模拟器不需从初始日期启运行，可以在选择的时间步继续运行。怎样做一重启？重启记录是可选的。如果需要做重启运行，则要在前次运行时产生。使用*WRST划*RESTART创建重启文件，它们必须位于输入/输出控制部分。然而，当井工作制度改变时*WRST可出现在井数据部分，*WRST指示写重启记录的频率。*RESTART表示当前值，否则，输入时间步数。例：*RESTART30当做重启时：a） 不改变原始油藏数据，但在输入/输出控制部分加一*RESTART关键字。b） 增加最大时间步数，如果必要删除*MAXSTEPSoc） 需要附加输入文件：输入重启文件和输入模拟结果文件。这些文件在前次运行期间产生d） 需要新的一组输出文件名。控制输出打印文件内容控制输出打印文件内容，使用：a） *WPRNb） *OUTPRN这些关键字可以出现在输入/输出控制部分，其参数也可在随后井数据部分修改。*WPRN表示写网格块数据、井数据和数值方法控制数据的频率。如果不希望输出网格或井数据，则频率置零。例如：*WPRN*WELL0*WPRN*GRLD0如果将这两个关键字删除，其缺省选择为每一次井变化时打印。这种打印控制会产生一个非常大的输出打印文件，会很快占满计算机的可用空间。*OUTPRN限制打印什么样的油藏信息和井信息。详细的网格信息和流体性质数据可使用*OUTPRN*RES产生控制模拟结果文件内容模拟结果文件（SRF）是模拟运行期间产生的数据文件。SRF文件由图形和控制读入，产生绘图输出。为控制是SRF的内容，使用：a） *WSRFb） *OUTSRF这两个关键字一般首次出现在输入/输出控制部分。*SWRF表示写网格数据和井数据的频率。如果不希望网格数据或井数据输出，则输出频率置零。网格系统描述描述网格系统，使用：*GRID*KDIR描述网格长度，使用TOC\o"1-5"\h\z*DJ*DJ*DK描述油藏深度和倾角，使用*DEPTH*DIP或*DTOP在以上关键字中，只有*KDIR和*DIP完全可选，并可能以从数据集中删除。上述列出的关键字必须在油藏描述部分。*GRID必须是这一部分的第一个关键字。*GRID表示网格系统类型，有四种选择：直角坐标、变深度/变厚度、径向井筒和角点网格。其中每个均要求I,J,K方向网格块数目：例如：*GRID*CART10106*GRID*VARI10106*GRID*RADIAL10115*GRID*CORNER10106第一个为标准直角人坐标10X10X6,第二个为变深度/变厚度网格10X10X6,第三个为径向网格，最后一个为角点网格。KDIR定义Z轴的方向，。缺省值为UP,第一层为最底层。死结点定义死结点有两种方法*NULL*POR两种关键字都必须在油藏描述部分，在数据集中必须有*POR,但*NULL可选。当用*NULL关键字表示无效网格时，0表示无效网格，1表示有效网可知，例如:*NULL*IJK1：101：10 1：311：41：3 1：30第二行覆盖第一行。*NULL是可选的，如果不出现，则假定全部网格均不有效网格。NULL关键字覆盖POR关键字。即使给一网格赋予非零孔隙度，而在*NULL关键字中赋给零值，则该网格在流动计算中也不予考虑。上述例子可以下列形式表示：例如：*POR*IJK1：10 1：10 1：3 031:4 1:3 1:3 0.0双孔隙度/双渗透率调用孔/双渗选择，可以使用(只能选一个);*DUALPOR*DUALPERM除裂缝介质标准网格信息外，还需输入：*DIFRAC*DJFRAC*DKFRAC*TRASFER在组分性质部分，还可以指定：*DIFFUSION孔隙度值需要对基质和裂缝输入例:*POR*MATRIX*IJK1：101：10 1：310.31:41:3 1:30.0*POR*FRACTURE*UK1:101:10 1:30.087:9 1:20.4*MOD8 7:8 1=0.45水层选择：水层模型如下：*AQUIFER*AQPROPAQFUNC步长过小或运行时间过长问题收敛失败可能由于：内迭代收敛失败牛顿收敛失败导致时步截短物质平衡误差如果在输出文件中发现“迭代不收敛”频繁出现，那么可以：1减小时间步长。将*DTMAX改小，或用*NORM*PRESS和*”0阳*SAUR减小第个时间步的改变量。2用关键子*ITERMAX增加迭代次数3使用*SDEGREE增加因子分解度。注意：这个方法增加了存贮需求。如果问题原因是最大改变量，则若出夙不频繁，问题不大。如果时间步大量重复出现，那么，应釆取以下措施：4检查岩石和PVT曲线是否出现非线性，该曲线应当平滑。5检查网格和其他性质6检查井的限制。一个有效的办法是对每一注入井总是指定一最大井底压力，对每一生产井总是指定一最小井底压力。7如果不收敛由于超过最大牛顿迭代次数引起，使用关键字*NEWTONCYC,增加牛顿循环次数。如果由于最大改变量振荡引起，达到最大迭代次数，那么，使蜚线性段平滑，或减小时间步长是最好的解决办法。8油藏是否在迭代之间出现在单相油和单相气之间浮动。如果出现浮动，设置适当的*PHASELD或*REFDENo9置油藏的某些区域或整个油藏为隐式。缺省开关符，*AIM*STAB检查网格的显式隐式开关，仅当该网格是一隐式格块的相邻格块。如果某些油藏区域变量发生剧烈变化，并且不是井的毗邻区域，或流体前沿运动太快，而在一个时间步扫过两个格块，那么置这些区域为隐式a） 当不使用垂向平衡初始化时。在某些情况下，这可能引起初始压力饱和度发生大的改变，即使全部井关井。当该情况发生时使用全隐式。b） 当出现气顶时。如果出现较强的指进，置气顶的底层为隐式，至少在指进出现的区域。c） 在具有极端高渗透率的格块，压力的微小改变可能引起饱和度很大变化。在这些区域，置格为全隐式。物质平衡误差可能由收敛误差超过设定的改变量而引起。在大多数情况下，*NORM和Converge的缺省值足够。然而，当模拟裂缝油藏或使用裂缝代表水平井时，推荐使用较小的值。对锥进问题，推荐使用小值。模拟一单相油藏有两种主要方法模拟气水两相的干气，单相气藏或一次接触混相溶解剂注入。A） 使用饱和压力使用关键字*PSAT置饱和压力为一很小的值例如，一1。在这种情况下，所有油藏流体是欠饱和的，并跳过每个可知块的闪蒸计算。所有流体性质仍然使用状态方程计算。*PSAT不影响井流动计算。B） 使用拟合组分选择1使用*PSEUDOo一般来说，置油藏油为一拟组分，每个注入流体作为不同的单个拟组分。2用*PHASEID将单相标为'油'或'气'。一般建议用'油'，从而流体相对渗透率依据油水相对渗透率曲线计算。该方法也不做网格闪蒸计算。另外，GEM还将单个组分流动方程分组为较小数目的拟组分流动方程。这个可以大大减少每个网格坟程的数目。水平井水平井可作为一线源（注入井）或汇（生产井）。GEM使用AXIX等流动图计算体滞留量，同时也计算井筒磨擦压力降。注意，当前方法并不能解决回流问题。垂赂平衡*VERTICAL控制垂向平衡计算。处理垂向平衡计算时，使用下列儿种方法：A） *VERTICAL*ON同时使用下列关键字*DWOC,*REFDEPTH,*REFPRES,*XGLOBALand*SWOCB） ^VERTICAL*GASCAP及*DWOC*DGOC*REFDEPTH*REFPRES*ZOIL*ZGAS*SWOC*VERTICAL*COMP及*DWOC*REFDEPTH*REFPRES*SWOC*CDEPTH*ZDEPTH这些关键字必须位于初始条件部分。VERTICALON主要用于欠饱和油藏，而VERTICAL*GASCAP可用于饱和油藏，也可用于欠饱和油藏。对于VERTICAL*ON,用重力一毛管压力平衡计算全部格压力和水饱和度。当网格包含油相和气相时，不做特殊处理。因此，在第一个时间步可能出夙剧烈的流体变化。*VERTICAL*GASCAP是更一般的选择，假定气顶油区处于平衡状态，用穗力一毛管压力平衡计算全部网格压力和油、气、水饱和度。闪蒸计算用于确油和气相组成。油相和气相依饱和度混合，然后计算网格总的组成。因此，油藏中总的组成可能随深度变化。*VERTIACL*COMP与用户指定的随度变化的组成一起，用于欠饱和油藏或者和油藏的重力初始化。这时，要指定水油界面深度正气参考压力。用户必须输入总的组成与深度表。用户还必须输入临界深度。在临界深度以下，单炷相网格为油网格，而在临界深度以一，单烧相网格为气相网格。如果不做垂赂平衡计算，输入VERTICAL*OFF及*PRES*SW*ZGLOBAL下述第一个例子表示对饱和油藏做垂计算的初始条件数据：**INITIALVERTICAL*GASCAP*REFDEPTH900.0*REFPRES3000.0*SWOC1250.0*DGOC78&0*ZOIL0.300.10 0.200.300.120.08*ZGAS0.780.190.030.00.00.0*SWOC1.0第二个例子表示不使用垂赂平衡计算的数据：*INITIALVERTICAL*OFF*PRES*CON3400.0*SW*CON0.2*GLOBAL*CON0.300.10 0.20 0.300.120.08井的定义*井使用以下关键字定义，注意严格按照关键字排序*WELL井号 ，井名，井类型定义^PRODUCER井号或者*INJECTOR井号或者*CYCLPROD井号在井类型以下，输入：*PWELLBORE（如果需要计算井筒压力，需要井筒数据）或*WELLBOREINCOMP（如果定义了注入井时需要，紧接着INJECTOR）*OPERATE（至少定义一个工作限制）*ONITOR（监视限制任）定义井的位置和射孔，使用：*GEOMETRY（可选的，在完井卡之前）*PERF（必需）这些关键字必须位于井数据部分。怎样关井？怎样重新打开井？在关井之前：（1）该井必须用以下方式完全定义：（a） 作约束和任意监视约束（b） 完井关键字在完全定义一口井以后，该井可以在使用TIME或DATE关键字的任一时间关闭。在关井之后任一时间可打开一口闭的井。WELL1'watermjector'WELL2r producer*WELL3'Solventnijector,**全部井具有相同的几何因了******radgefacwfiacskin*GEOMETRY*K0.250.34 1.00.0两口源程序入井初始关井，只有生产井开井。这咱工作制度持续两年。*INJECTOR1*INCOMP*WATER*OPERATE*MAX*STW12000.0*OPERATE*MAX*BHP10000.0*PERF*GEO1**ifjfkfff1 1 3 1.0**关1号井SHUTIN1*PRODUVER2*OPERATE*MAX*STO12000.0*OPERATE*MIN*BHP1000.0*MIONTOR*GOR10000.0*STOP*MONITOR*WCUT0.8330002*STOP*PERF*GEO2**ifjfKfff7711.0*INJECTOR3*INCOMP*SOLVENT0.780.190.030,00.00.0*OPERAE*MAX*STG1.2E7*OPERATE*MAX*BHP10000.0*PERF*GEO3**ifjfkfff]1 1 3 1.0-号井和1号井在相同网格射孔**关3号井*SHUTIN3*TIME730.000*OPEN1 **打开注水井*TIME1095.0糾在这个方案中，注水和注溶解剂井在同一网格定义；在打开第二个以前，关闭第一个。*SHUTIM1 **关注水井*OPEN3和开注溶解剂井*TIME同14600*SHUTIN3 **关注溶解剂井*OPEN1 **打开注水井*TIME18250SHUTIN1 **关注水井*OPEN3和开注溶解剂井井筒模型说明井筒模型时，需要：PWELLBORE用于生产井，IWELLBORE用于注入井，位于井类型卡之后。当要求井口压力时需要使用井筒模型，用于：a）仅是一种信息；b）当井口压力作为一个约束条件时。例如：对2号注入井建立井筒模型*INJECTOR2*IWELLBORE**wepthlengthroughwhtempbhtempwi1100.1100. 1365200.3490.033工作和监视约束OPERATE和MONITOR说明一给定井的约朿条件。至少要求一个工作约束，监视约束是可选的。每一口井引入一新的未知变量pbh（井底流压），要求一约束方程确定该变量。在工作和监视约束表中第一个工作约束条件为主工作约束。模拟器首先在主工作约束下运行，同时监视其它约束条件。当监视约束之一违反，而且使用CONT关键字，则该约束变为工作约束。如果不止一个工作约束违反，那么使用第一个工作约束。生产井对一生产井，应当：a） 在产量约束下生产（主工作约束）b） 在最小井底流压约束下生产。如果生产井为产油井，则选一油产量约束。如果生产井为产气井，则选一气产量约束。对一生产井使用的约束可以是一最小产量。例如：^PRODUCER1*OPERATE*MAC*STO120000.0*CONT*OPERATE*MIN*BHP1500.0*CONT这个例子表明：a） 使用油产量作为油井的主约束b） 使用井底流压作为从约束如果有一条违反，运行继续，工作约束改变为刚刚被违反的约朿。CONT是缺省值不需输入。注入井对一注入井，应选：a） 最大注入时作为主工作约束b） 最大井底流压约束如果为注水井，选水产量作为约束。如果为注气井，则选气井，则选产量。例如：*INJECTOR2*OPERATE*MAX*STW1000.0*STOP*OPERATE*MAX*BHP2250.0*STOP该例表示：a） 注水井的注水量是主约束b） 同时监视爐底流压，作为从约束只要有一个违反，则模拟停止.监视约束监视约束格式与Operate格式相同。而且，如果同时有多个约朿违反，则采取约束表中第一个约束的动作。极力推荐监控生产井的GOR和含水率：这个可以避免作业运行中的某些问题。例如：*PRODUCER1*OPERATE*MAX*OIL1200.0*CONT*OPERATE*MIN*BHP2500.0*CONT*MONITOR*MAX*GOR15000.0*STOP*MONITOR*MIN*WCUT60.0*SHUTIN如果在一时间步，GOR和含水率均违反，模拟停止。输入井指数井指数可由下述两种方法获得：a） 从井模型的井几何因子。使用下述关键字：*GEOMETRY*PERF*GEOb） 使用关键字直接指定：*PERF这些关键字位于数据部分。*GEMTRY指定内部计算井指数的必要参数。完井关键字PERF包括完井段的位置。如果*GEO子关键字不在*PERF中，则应指定井指数。*PERF可用作水平井、斜井和垂直井。例如：*WELL1'12-09-18-56'*PERF1**ifjfkfwisetn112：41.241定义井指数为1.24或例如：*WELL1'12-09-18-56'**radgefacwfiacskin*GEOMETRY*K37524881.00.0**完井关键字必须紧接着儿何关键字*PERF*GEO1**lfjfkfffsetn1 12：411停止模拟运行可使用关键字*STOP中止模拟运行。另一方法是，在一井改变之后，由*。丁\¥£!±指定的第一个时间步长为0。例：*DATE19980908*DTWELL0.0怎样建立井数据当首次输入井数据时，数据必须依据以下顺序：1要法语*TIME或*DATE卡2定义一*DTWEL值，该数据是井定义之后，立刻使用的第一个时间步长。3使用*WELL标识所有新井每组井定义由如下部分组成：4a）使用*producer和*INJECTOR定义一口新井或主工作制度改变。b） 如果选择*WHP作为约束之一，必须指明井筒模型要求。c） 定义注入流体组d） 定义该井的工作或监控约束5使用*GEOMETRY定义井位置、儿何因子或井指数，使用*PERF定义完井段。6在第1步到第5步后可使用*SHUTN7使用*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卡。假定10X10X12直角坐标系统，井数据的构成如下：糾井数据*RUN*DATE199011*DTWELL0.025*WELLrOILPRODUCER5*PRODUVER1**吸有生产井的工作和监控约束*OPERATE*MAX*STO1000.00*OPERATE*MIN*BHP3000.00*MONITOR*MAC*WCUT0.40**radgefacwfiacskin*GEOMETRY*K0.250.341.00.0*PERF*GEO1**ifjfkfff1131*TIME10.0**在10天时，主约束值由1000改为100*ALTER1100.0**在50天时，主约束值再改为50.0*TIME50.0ALTER11000.0和定义一注入井，而井对生产井再次修改主约束*TIME720*INJECTOR1*INCOMP*WATER*OPERATE*MAC*STW12000.0*OPERATE*MAX*BHP10000.0**radgeofacwfiacskin*GEOMETRY*K0.300.341.00.0*PERF*GEO1**ifjfkfff773 1ALTER2100.0**在第10年停止*DATE19901231*STOP其它部分循环数据其它部分也有一些关键字可在模拟运行期间改变，这些关键字必须位于*DATE或*TIME卡之后，并且后一个值覆盖前一个值。在输入/输出控制部分出现的这些关键字可以再赋值：*WRST*WPRN*ORTPRN*WSRF*OUTSRF*GRID*OUTSRFGRID可以改变，但OUTSRF*WELL不能改变。数值方法控制部分，出现的关键字可以再赋值的有：*DTMAX*DTMIN关键字输入系统关键字系统介始关键字标识符关键字前可带星号“和”或不带星号如：poi或POR或*POR两个星号表示注释关键字顺序所有关键字分成关键字组。不同组关键字必须严格按手册中顺序排列。若无特殊说明同一组内关键字内关键字可以改变顺序。*LIST,*NOLIST,*INCLUDE,*COMMENTfl*TRANSLATE可以出现在数据文件的任何地方。有些关键字可以出现在关键字组内，出可在循环数据中。每一关键字的描述说明该关键字是可选的或必须的。有些关键字与其它关键字一起使用时是可选的或必须的。如果可选的关键字在数据文件中不出现，则使用其缺省值。串关键字若一关键字不是另一关键字或另一组关键字，不能与主关键字在同一行。每一行是一个关键字加上它的子关键字，除非另有说明。例如，水层选项需要3个关键字。*AQUIFER具有一系列关键字，并需要定义水层的位置。该选择还需要水层性质（*AQPROP）和水层影响函数（*AQFUNC）o三个关键字是独立的实体，不能放在同一行。不正确的关键字表示：*AQUIFER*BOUNDARY*AQPROP1000.4.002230*AQFUNC0.010.11•••••••行长度一行最多130个字符。界定符空格，逗号或新行可作为关键字，数字或字符串的分隔符。字符串字符串用单引号标出，例如'5-35-48-W5',内嵌的单引号不能使用。可内嵌双引号。4TinsIsTHE"nght“Way'更名关键字：对主关键字可使用TRANSLATE更名注释（可选）格式桝注释文本空行（可选）分隔不同组关键字。检查错误扫描（可选）仅用于数据格式检查，不做运行。格式：*CHECKONLY该关键字必须放在I/O控制部分。包含文件（可选）格式：*INCLUDE'文件名'控制数据文件列表（可选）格式*LIST*NOLIST缺省值：*LIST改变注释符：（可选）格式：*COMMENT，ab'缺省值：*COMMENT''命名用更名规则改变关键字（可选）格式：TRANSLATE，用户关键字'，CMG关键字'一个关键字字可有几个别名：例如：TRANSLATE'KX''PERMI'TRANSLATE4X—permerbility,'PERMI'TRANSLATE'PERMX''PERMI'网格性质数组输入数组读选择数组赋值由五部分组成，其中两个个任选。句法：网格数组（数组限制符）读选择数据（数组修改符）定义：网格数级要赋值的性制，例哪：*POR数组限制符对网格网格的不同元素赋值（例如，基岩和裂缝），数组限制符可选。选择是:*MATRIX*FRACTURE*EQUALSI如果没有数组限制符，则假定％*MATRIX读选择：*CON*IVAR*JVAR*KVAR*ALL*UK数据实际网格数组值数组修改符*MOD输入基岩网格性质关键字：*MATRIX输入裂缝网格性质关键字：*FRACTUREJ和K方赂数据与I方向数据相同关键字：*EQUALSI例1：单孔隙系统*PERMI*CON100.0*PERMJ*EQUALSI*PERMK*EQUALSI*2例2：天然裂缝*PERMI*MATRIX*CON100.0*PERMI*FRACTURE*CON10000*PERMJ*EQUALSI*PERMK*EQUALSI*2常值数组*CONIJK标识的数组输入关键字：*IJKII：12JI：J2K1：K2值数组输入值随I方赂变化关键字：*IVAR值(1),……,值(m)数组输入值随J方向变化关键字L*JVAR值(1)……，值(刃)数组输入值随K方向变化关键字：*KVAR值(1),……值(nk)值随全部网可知改变关键字：*ALL值(1),……值(mnjnk)数值I方向下标增加最快，J方向次之，k方赂最慢。修改数组数据格式：*MOD11：12jl：j2kl：k2(+)value(-)(*)(/)(=)或MOD(*)value(-)(十)(/)(=)输入输出控制输入输出控制说明GEM可以生成两种不同格式输出显式文件，SRF和SR2,缺省输出文件仅有SRF文件。使关键字*RESULTFILE可使GEM仅生成SR2或生成SRF和SR2文件。输入输出文件名格式：*FILENAMES文件类型名称选择或：*FILENAME文件类型名称选择其中文件类型是下列之一或更多：*OUTPUT*RESTARTIN*RESTARTOUT*SRFIN*SRFOUT*INDEX—OUT*MAIN—RESULTS—OUT*INDEX—IN*MAIN—RESULTS—IN*SR2ASCOUT对文件类型为*INDEX—OUT,*MAIN—RESULTS—OUT,*MAIN~RESULTS-IN,可以指定下列名称选择之一。文件名'*PROMPT*SCREEN对文件类型*INDEX—IN和*RESTARTIN,可以指定下列名称选择之一：文件名*PROMPT格式：FILENAMES(OUTPUT)(，')('文件名')(*PROMPT)(*SCREEN)(*RESTARTIN)(文件名)(*PROMPT)('')如果指定*INDEXIN(*RESTARTOUT)(，')('文件名')(^PROMPT)(*SRFIN)('')('文件名')(*PROMPT)(*NONE)(*SRFOUT)(，')('文件名')(*PROMPT)(INDEX—OUT)('')('文件名')(*PROMPT)(*MAIN-RESULTS-OUT)('')('文件名')(*PROMPT)(*INDEX-IN)('文件名')(*PROMPT)(*MAIN-RESULTS-IN)('')('文件名f')(*PROMPT)(*SR2ASCOUT)('')(文件名f')(*PROMPT)定义：*FILENAMES文件名关键字*OUTPUT模拟结果主输出文件RESTARTIN如果选RESTART,该关键字扌旨明强制读入的重启动数据文件。RESTARTOUT如果选WRST,为强制定入的重启动数据文件。SRFIN如果选RESTART,为读入的SRF数据文件。如果没有该数据文件，使*NONESRFOUT如果已选SRF写选择WSRF,表示要写入的模拟结果文件（SRF）INDEX—OUT如果指定RESUKTFILE*SR2,定入的SR2ASCII下标文件。*MAIN—RESULTS—OUT定入的SR2主二进制结果文件。*INDEX—IN重启运行读入的SR2ASCII下标文件*MAIN-RESULTS-IN重启运行读入的SR2主二进制结果文件。*SR2ASCOUTSR2主二进制结果的SACII文件。'空串，表示使用缺省文件名。'文件名'文件名。PROMPT提示用户从键般输入文件名。如果数据文件中不使用FILENAMES,这是缺省值。*SCREEN该文件在屏幕上显示（标准输出设备）。*NONE表示不提供输入文件。当输入文件丢失，或用户希望不考虑以前历史数据时，该选择是有用的。缺省文件名仅当名称选择使用''时，用缺省文件名缺省文件名具有形式'xxxxxxx.yyy',其中'XXXXXXXX'为根，’YYY'为后缀，后缀依赖于有关的文件。例如：’OUT'为主输出文件，’RST为重启，'srf为SRF文件。根据的取法规则如下：（a）取掉包括'。'分界符的后缀；（b）除最右边［了符串和'一'号以外，取掉所有字符；（c）取掉80列以后字符，例如：DOS路径，\CMGQATA\TEST_3.DAT'具有根'TEST_3'下列表总结了缺省文件名的根和后缀：文件根基二后缀*OUTPUT数据文件名Out*INDEX—OUT*OUTPUTIlf*RESTARTOUT*OUTPUTRst*MAIN—RESULTS—OUT*OUTPUTMif*INDEX—IN（无缺省值）*RESTARTIN*INDEX—INRst*MAIN—RESULTS—OUT*INDEX—INMif*SRFOUT*OUTPUTSrf*MAIN*RESTARTINSff使用这套缺省系统，用户对每一次运行，仅修改*OUTPUT和*INDEX-IN或*RESTARTIN文件名，就可以执行一系列'自举'重启运行例1：数据文件是'Cycle。Dat\因此*OUTPUT缺省是'Cycle.out如果使用*OUTPUT缺省值，那么*RESTARTOUT缺省值是'Tycle.rsf,*SRFOUT缺省值是''Cycle.sif.例2：数据'Cycle.dat'第一次运行：**FILEMAMES*OUTPUTtuall*RESTARTOUT缺省值是<tnall.ist,*SRFOUT缺省值是'心all.sif'第二次运行：在'Cycle.dat，中增加*RESTART关键字，使用*FILEMAMES^OUTPUT，'tnal2.out',*FILENAMES*RESTARTIN' 'triall.ist'。不要更名或删除文件。*RESTARTOUT缺省值是：，＜trial2.rst,*SRFOUT缺省值是''Mall.sif'*SRFIN缺省值是''Mall.srf'例3：数据文件是：’Cycle.daf,因此，*OUTPUT缺省是Tycle.oufo如果使用*OUTPUT缺省值，那么*INDEX—OUT缺省值'Cycle.uf,MAIN—RESULTS—OUT缺省值是，Cvcle.mif条件：*FILENAMES必须是数据文件的第一个关键字。否则，将不提示用户输入文件名。文件名不能在I/O控制部分的任一位置指定，它们必须是该部分第一个关键字。所有输入出文件打开状态为’unknown 因此，可写入任何信息。所有输入文件打开状态为'OLD',因此，模拟运行时，必须存在。项目主标题（可选）格式：*TITLE1stung定义：suing字符串，最多40个字符。条件：该关键字必须位于输入/输出控制部分，在数据文件的开始。项目次标题：（可选）格式：*TITLE2stung定义：suing字串符（最多80个字符）项目第三条标题（可选）格式：*TITLE3定义：同*TITLE2方案标识符目的CASEID标识物定的方案运行，它还用于模拟结果文件标识绘图的数据曲线。格式：*CASEIDsuing定义：suing字符串，最多8个字符例：*CASEED'case21c'CHECKONLY（可选）目的CJECKONLY表示对数据做格式检查，不做运行。格式：*CHECKONLY条件：该关键字位于I/O控制部。在数据文件中越早越好，但它不能放在重启运行灵敏据中。日期格式：（可选）目的：*DIARY指明在模拟运行期间，输入出现到幕输出文件的时间步信息类型。格式：*DIARY（*CHANGES）（*WELL—INFO）定义：*CHANGES该选择表示在每一个时间步后，佃出到终端屏幕和输出文件中的信息为该时间步压力、饱和度和组分组成最大变化量的大小和位置。*WELL—INFO该选择表示在每一个时间步后，输出到终端屏幕和输出文件中的信息为当前油田油、气、水产量，全油田含水率和气油比、油田气、水注入量和平均油藏压力。缺省值：^CHANGES输入数据单位(可选)格式：*INUNIT(*SI)(*FIFLD)缺省值:*SI单位制表*SI0*FIELD时间天天液体体积M3Ft3或bbl气体体积M3Ft3压力KpaPsi绝对渗透率MdMd粘度CpCp质量密度Kg/M3Lbni/ft3摩尔密度Gmol/M3Gmol/ft3长度、距离MFt温度°C°F临界压力atmAtm临界温度临界体积M3/(k-mol)M3/(k-gmol)数据范围检査格式：*RANGECHECK(*0N)(*FF)定义*0N打开范围检查*OFF关闭范围检查缺省值：*0N说明：检查大部分输入数据以确定它是否在允许范围之内。指定*RANGECHECK*OFF将取消非临界数据范围检査。*RANGECHECK*OFF还将取消所有警告信息的打印。最大错误信息数(可选)格式：*MACERROR数字缺省值：20条件：数字范围为1到100重启时间步(可选)格式：*RESTART数字定义：数字一时间步数。缺省值：如果在输入数据文件中，没有该关键字，则不读重启记录，模拟从零开始。如果使用*RESTART,但其后没有时间步，则重启在输入重启文件的最后一个时间步开始。重启记录写频率(可选)格式：*WRST(fieq)(*TIME)*REWIND(num)定义：fieq频率一正整数或零。说明每几个时间步写一次重启记录。如果频率为零，不写重启记录，不打开输出重启文件。重启频率不能在循环数据中重启意义。如果频率是正数，打开输出重启文件，在循环数据中可重定义频率，重启记录在初始化(零时间步)，模拟终止的最后一个时间步，及中间按定义的频率写入。TIME按循环数据TIMEn£*DATE卡指定的时间写重启文件。Num数字一反绕重启文件的频率。该数字是重启文件反绕前累计的重启记录最大数。该数字必须大于零。如果该数字是1,则只有最后一个重启记录可用。缺省值：如果不使用*WRST,则在仃吨或*。^^指定的时间步写重启记录；如果不使用*REWIND,重启文件不反绕；如果缺num,假定num=l条件：该关键字可以在I/O控制部分，也可做为循环数据的一部分。因此，只要开始不出现*WRSTO重启文件记录频率可在模拟期间改变。如果指定*WRSTO,不打开任何重启文件，在循环数据部分也不能对*WRST再定义。如果开始不需要重启记录，但以后需要，则使用：*WRST100000并在适当的时间对WRST再定义。输出打印频率(可选)格式：*RPRN(*WELL)(fieq)(*TINE)(*GRLD)(freq)(*TIME)(*ITER)(*ALL)(*MATRIX)(*NEWTON)(*BRIEF)(*NONE)定义：*WELLfieq每fleq个时间步写井结果到输出文件，其中fieq为整数；如果fieq为零，不写井结果。WELLTIME按循环数据中TIME卡或DATE卡指定时间写井结果到输出文件。*GRIDfieq每fieq个时间步写网格结果到输出文件，其中fieq为整数。如果fieq为零，不写网格结果。GRIDTIME按循环数据中TIME或DATE卡扌旨定时间写网格结果到输出文件。*ITER*ALL提供每个时间步矩阵解法和牛顿迭代的全部细节。*ITER*MATRIX仅提供每个时间步矩阵解法细节*ITER*BRIEF提供时间步收敛状态汇总*ITER*NONE不提供时间步状态信息。缺省值：每个*TIME或*DATE卡写网格结果。每当写网格结果时，写井结果。*ITER的缺省值是*BRIEFo条件：打印文件中细节的多少可在模拟期间改变。输出打印文件项(可选)格式：*OUTPRN(*WELL)(*ALL)(*BRIEF)严GRID)(*NONE)(gild—list)(*RES)(*ALL)(*NONE)定义：*WELL该键字表示井结果将写到输出打印文件*GRID该关键字表示网格结果将写到输出打印文件*ES该关键字表示在模拟运行开始时打凶输入的油藏性质。ALL写所有可能的变量到输出文件。不能与GRID一起使用。*BRIEF写井变量汇总。*NONE没有变量写到输出文件。grid-list一系列字符串标识哪个变量写到输出文件*GRIDRgrid—list是：*PES油藏压力*DROP从时间零开始的压力降*SO油饱和度*SG气饱和度*SW水饱和度*TSO圈闭油饱和度*RHOO油的摩尔密度*RHOG气的摩尔密度*FRG油/气系统的气相摩尔分数*PCW水的毛宇航局丈夫力*PCG气的毛管压力*KRO油的相对深透率*KRG气的相对深透率*KRW水的相对浓透率*VISO油的粘度*VISG气的粘度*MWO油的分子量*MWG气的分子量*DENO油质量密度*DENG气质量密度*DENW水质量密度*SIG油气之间界面张力吃m组分M总摩尔分数*Xin油相中组分m的摩尔分数*Ym气相中组分m的摩尔分数*Wm水相中组分m的摩尔分数*Kin组分m的气/油K值条件：模拟期间可改变输出量的多少。结果文件格式选择(可选)目的：*RESULTFILE表示写入的图形显式文件格式——SRF或SR2或两者兼有之。格式：RESULTFILE(SRF)(SR2)(*BOTH)定义：*SRF仅生产sif文件SR2产生iff和mf文件*BOTH产生三个文件缺省值：*RESULTFILE的缺省形式是*RESULTFILESR2。如果RESULTFILE后无关键字，则等价于*RESULTFILE*SR2附加的普通文本SR2主文件(可选)目的：*SR2ASCII表示若图形显示文件以SR2格式定入，也可产生一二进制主结果文件的ASCII副本，扩展名为asc格式：*SR2ASCD模拟结果定频率(可选)格式：*WSRF(*WELL)fieq(*TINE)(*GRLD)fieq(*TIME)缺省值：如果数据文件中缺省*WSRF,则在每个时间步写井数据，在所有井修改的进间写明网格结果到SRF或SR2文件中。缺省值与下列数据行等价。*WSRF*WELL1WSRF*GRID*TIME条件：该关健字也可以放在循环数据中。说明：如果指定如下形式：*WSRF*WELL0*WSRF*GRLD0则不打开任何SRF或SR2文件，在循环数据部分也不能对，WSRF重定义。如果开始不需要结果信息，而在随后的模拟运行中需要，可对freq输入一大数，例如：100000,然后在循环数据的适当位置覆盖掉。模拟结果文件项(可选)格式：*XDR(*ON/*OFF)OUTSRF(WELL)(welljist)严GRID)*NONE(纠d」ist)(*RES)(*ALL)(*NONE)定义：*XDR二井制文件可以所在平台二进制格式，也可以外部数据格式(XDR)写入。使用XDR可以使SR2二进制文件在一个平台上写入，而在另一平台上读出。例如，在UNIX工作站服务器上产生SR2文件，可用PC上的Results或RepotWntei访问。如果SR2以XDR格式定入，则关键字XDR位于RRF文件的顶部*WELL该关键字表示将井或定变量写到SRF或SR2文件中。至少要对每口井生成下列井变量：井底流压地面油产量地面气产量地面水产量累计油产量/注入量累计气产量/注入量累计水产量/注入量气油比水油比*GRID该关键字表示将网格结果写入模拟结果文件SRF或SR2RES该关键字单独或与ALL关键字一起表示将某些油藏始条件写入SR2文件的零时记录，用于Results图形显示，比如，渗透率、参考孔隙度、网格相对渗透率端点。*NONE无变量写入ALL初始RES性质的所有变量定入SR2文件中。Welljist一系列关键字表示哪一口井或哪一个变量写入SRF或SR2文件中。这些值按时间顺序显示。尽管出现在井变量列表中，许多变量与指定网格有关。*WELL中的Welljist的有效关键字如下。I、J、K指格的I、J、K值。*DELT时间步长度*PAVG平均油藏压力*PRES_I_J_K网格I、J、K的压力*SO_I_J_K网格I、J、K的油饱和度*SG_I_J_K网格I、J、K的气饱和度*SW_I_J_K网格I、J、K的水饱和度*RHOO_I_J_K网格I、J、K的油摩密度*RHOG_I_J_K网格I、J、K气摩尔密度*DENO_I_J_K网格I、J、K的油质量密度*DENG_I_J_K网格I、J、K的气质量密度*KRO_I_J_K网格I、J、K油相对渗透率*KRG_I_J_K网格I、J、K气相对渗透率*KRW_I_J_K水相对渗透率*PCG_I_J_K网格I、J、K气毛管压力*PCW_I_J_K网格I、J、K水毛管压力*SIG_I_J_K网格I、J、K油气之间界面张力*TSO_I_J_K网格I、J、K圈闭油饱和度*FRG_I_J_K网格I、J、K的气相摩尔分数*VISO_I_J_K网格I、J、K油粘度*VISG_I_J_K网格I、J、K气粘度*MWO_I_J_K网格I、J、K的油分子量*MWG_I_J_K水粘度BHP_IW_IK井IW在IL层的井筒压力*TROW_rW_IL井IW在IL层的油传导率*TRGW_IAV_IL在IW在IL层的气传导率*TRWW_rW_IL井在IW在IL层的水传导率*WSPR_IW井IW的井定义*WKRO_r\V_IL井IW在IL层的油相对渗透率*WKRG_r\V_IL井IW在IL层的气相对渗透率*WKRW_IW_IL井IW在IL层的水相对渗透率网格I、J、K组分M的气/油K值网格I、J、K油相中组分M组成网格I、J、K气相中组分M组成吃网格I、J、K组分M的总组成*RMOLM_I_J_K网格I、J、K单位网格体积组分M的摩尔数*DROP从时刻零开始的压力降*ZWELM_IW井IW组分M的总组成*FOIP裂琢中油储量*FGIP裂缝中气储量*FWIP裂缝中水储量*MOIP基岩中油储量*MGIP基岩中气储量*MWIP基岩中水储量*TOIP油藏中油的总储量*TGIP油藏中气的总储量*TWIP油藏中水的总储量*RECO总的油釆收率*RECG总的气釆收率*gnd_list一系列关键字标识哪一个网格变量写到SRF或SR2文件中，对油藏中每个网格写入一个值。*GRID的gndjist为:*PRES油藏压另*DATUMPRES校正到基准深度的网格压力*DROP从时间零开始的压力降*SO油饱和度*SG气饱和度*SW水饱和度*TSO圈闭油饱和度*RHOO油摩尔密度*RHOG气摩尔密度*FRG油气系统中气相摩尔分数*PCW水的毛管压力*PCG气的毛管压力*KRO油的相对渗透率*KRG气的相对渗透率*KRW水的相对渗透率*VISO油的粘度*VISG气的粘度*MWO油的子分量*MWG气的分子量*DENO油质量密度*DENG气质量密度*DENW水质量密度*SIG油气界面张力*Zm组分m总摩尔分数*Xin油相中级分M的摩尔分数*Ym气相中组分M的摩洋分数*Wm水相中组分M的摩尔分数*Kin组分M的气/油K值缺省值：如果缺省*OUTSRFWELL关键字，在上述WELL下指定的井性质按*WSRFWELL中指定的频率定入SRF文件中。如果缺省*OUTSRF*RES,则不将初始信息定入SR2文件中。而*OUTSRFRES等价于OUTSRF*RES*ALL条件：该关键字也可在模拟期间修改。井、水层和区域汇总表(可选)格式：*SUMMARY油藏描述网格卡(必须)格式：*GRID(*CART)ni刃ilk(*VARI)mnj11k(*CORNER)niq11k(*RADIAL)mnj11k*RWrw定义：*CART标准矩形网可知*\4RI矩形网格，深度和厚度变化*CORNER角点网格*RADIAL径向网格*RW井筒半径关键字rw井筒半径，非零值。缺省值：如果RW缺省或赋零值，则使用0.0762m(0.25ft)条件：该关键字应是油藏描述关键字组第一个关键字。说明：在*VARI关键字中，同一层的相邻网格间的全部流动保持不变，即使它们处于不同深度。半径“BV”所表示的井筒体积从径赂网格的所有层的最内环取掉，因而最内环有一与井筒相连的圆柱形内边界。网格类型IJK层顺序*CARTXYZ向上（缺省）*VARIXYZ向上（缺省）*CORNERXYZ总是向下*RADIALR0Z向上（缺省）K下标的方向（可选）格式：*KDIR（*UP）（*DOWN）I方向网格长度（条件）格式：*DIJ方向网格长度（条件）格式：*DJK方向网格长度（条件）格式：*DK深度（条件）目的：*DEPTH后输入单个网格的油藏深度，该深度通常以网格中心为准，除非*TOP出现。格式：*DEPTH（*TOP）Ijkdepth（*CENTER）定义：*TOP指到参考网格顶部顶面中心深度*CENTER指到参考网格中心深度。Ijk参考网格I、j、k方向下标。Depth深度，值可正可负，取决于参考面的位置。到网格顶部深度（条件）格式：*DTOP条件：必须输入（mnj）个值。到有效百姓工中心深度（条件）格式：*PAYDEPTH网格倾角（条件）格式：*DIPldipjdip（对*GRID*CART）DIPldip（对GRID*RADIAL）定义：idipX倾角ljdipY轴倾角说明：idip和jdip的范围为一90°到+90°角点网格的角点深度（条件）格式：*ZCORN角点网格的层面角点位置（条件）格式：*XCORN或*丫（30阳基于线的角点位置（条件）格式：*COORD角点网格完全角点位置格式：*CORNERS双孔隙度（可选）格式：*DUALPOR双渗透率（可选）格式：*DUALPORM基岩一裂缝交换项计算（条件）目的：TRANSFER表示双孔或双渗油藏中使用的基岩一裂缝模型。格式：TRANSFERntran定义：ntian基岩一裂缝模型号：=0；无特殊处理。该选择不能精确模拟重力驱过程。一般不推荐该选择。=1;完全相分离模型=2；拟毛管压力模型=3；具有相界面校正的拟毛管压力模型缺省值：3条件：仅对于由Dualpor或DUALPERM关键字指定的天然裂缝油藏要求。裂缝间距（条件）格式：*DIFRAC*DJFRAC*DKFRAC条件：若裂缝间距为零，表示无裂缝面与该轴垂直，即裂缝间距无穷大。说明：若某一网格DIFRAC=DJFRAC=DKFRAC=0,则假定无裂缝孔隙度。死结点标识符网格（可选）格式：*NULL条件：紧接着数组值，对死结点是零，对有效网格为1。说明：其它性质也可使网格无效。如果一个网格标为无效网格，设一非无效性质,包括一非零孔隙度，也不能使网格成为有效网格，所赋性质不使用。如果使用一双孔隙模型，用*NULL可选择使用两个孔隙度。使用*NULL,无*MATRIX之一无效，然后再使另一个无效具有与上述方法相同的效果。使用*NULLMATRIX并且给DIFRAC,*DJFRAC赋零值，其结果亦相同。使用*NULLMATRIX或NULL*F[actuie仅使用得一种孔隙无效，流动仍可在另一种孔隙中发生。需要注意的是，岩石尖灭覆盖*NULL设置。这就意味着，如果一网格出夙在*PINCHOUT中，或具有零厚度，流体在其中会发生流动，而不管该网网格己被定义成无效格块（或给定零孔隙度）。孔隙度（必须）格式：*POR岩石压缩系数（必须）格式：*PRPOR（*MATRIXief_piessure（*FRACTURE）CPOR（MATRIX）cpoi（*FRACTURE）DCPOR（MATRIX）dcopoi（*FRACTURE）定义：ref_pressure参考压力（Kpa|psi）cpoi依赖于压力的地层孔隙度(岩石压缩系数)(l/Kpa|l/psi)dcpo［依赖于压力的岩石压缩系数(1Kpa**2|l/|psi**2)缺小值：DCPOR的缺省值是0。说明：在压力为P时，包含流体的地层孔隙度po［的计算使用下列公式：por(P)=poi*［l+cpoi*(p-ief_piessure)十dcpor*(P—ref_pressui)**2其中“poi”是*POR关键字输入的原始孔隙度，P是油相压力。使用cpor和ref_piessui-计算所有网格的孔隙度，除非在岩石流体性质中由*CROCK定义了压缩系数。在双孔模型中，应输入基岩参考孔隙度、裂缝参考孔隙度和岩石压缩系数。孔隙体积修改(可选)格式：*VOLMOD说明：在双孔模型中，或者同时使用VOLMODMATRIX和VOLMOD*FRACTURE或者都不用。当油藏边界穿过一网网格时，使得网格的某一部位于油藏之外，或者模拟一井网驱动时，遇到流动单元边界，需要调整孔隙体积。渗透率（必须）格式：*PERMI*PERNJ*PERMK传导率乘子（可选）格式：*TRANSI*TRANSJ*TRANSK缺省值：1.0条件：也可位于循环数据中。说明：因为网格间流动计算涉及到两个网格，因此需要确定哪一个网格对传导率因子做贡献。如果考虑两个网格间的流动，则选择在I方向为最小的I下标的网格，J方向为最小数点的J下标的网格，K方向为最小的K下标的网格。即这些具有方向性的传导率乘子应用于该方向具有较大下标邻近网格的界面。同理，当使用*KDIR中*!/（缺省值）时，*TRANSK适用于K方向增加的网格表面：如果选*KDIR*DOWN,则为K减少的网格表面。双孔隙度对双孔模型，考虑非空基岩块，同时在某一方向具有非零裂缝间隔，那么*MATRIX指定的传导率用于基岩裂缝间流动。双渗透率基岩传导率乘子用于基岩与基岩间流动。典型的情况是，I、J方向传导率乘子是零，而K方向非零，因为相密度不同使得最重要的双渗透率效果体现在垂向上。有效厚度（可选）格式：*NETPAY厚度因子（可选）目的：*NETROSS标识有效厚度与总厚度之比。格式：*NETGROSS条件：如果使用*NETPAY,则不应使用该关键字。尖灭（可选）格式：*PINCHOUT订：12jl：j2kl：k2条件：零厚度网格（DK为零）自动尖灭。说明：*PINCHOUT模拟尖灭层。这些网网格不参与任何流动计算，即：定为无效网格，但是，允许垂向流动。尖灭网格用于从某些区域中抽掉一些层，它们不在该区域而有基它区域。注意：尖灭说明覆盖*NUL说明。如果一网格为尖灭或零厚度，则流动会穿过该网格，尽管该网格己用*NULL定为无效网格。断层（可选）格式：*FAULTtluowII:12Jl：J2定义：throw从地质上來说，“山row”是地质断层块与邻近油藏岩石之间的深度差。说明：当一部分油藏与另一部分油藏错位时，就形成了地质断层。这些错位部分形成地质层块。断层网格之间的平面流动不遵循通常的地质层位。水层（可选）格式：*AQUIFER(*BOUNDARY)(*BOTTOM)(*REGION)il：i2jl：j2kl：k2(*I)(*J)(*K)(*U)(*JK)(*IK)(*UK)*AQPROP厚度孔隙度渗透率半径角度定义：*BOUNDARY水层网格均为边缘网格，而非顶部或底部。(*BOTTOM)水层网格为油藏度部(*REGION水层网格为一网格区域*1水层与油藏边界面在I方向*J水层与油藏边界面在J方向*K水层与油藏边界面在K方向*U水层与油藏边界面在I和J方向*JK水层与油藏边界面在J和K方向*IK水层与油藏边界面在I和K方向*UK水层与油藏边界面在I、J、K方向厚度水层厚度，缺省值为平均油藏厚度。孔隙度水层孔隙度，无缺省值渗透率水层渗透率，无缺省值半径有效油藏半径，缺省值依据油藏儿何因子计算角度水层影响角度，表示为一圆周的分数缺省值选择关键字。*AQUIFER无缺省值。*AQPROP的缺省值：厚度：油藏厚度孔隙度：无缺省值渗透率：无缺省值半径：对*BOUNDARY,为油藏半径。对*BOTTOM,为油藏高度的一半。角度：以*BOUNDARY而言，为与水层相边的油藏半径分数。对*BOTTOM而言，为界面长度/(pio油藏高度)。说明：在*REGION之后，可以输入多个行，定义复杂的水层。这些独立的行定义同一个水层。例如：*AQUIFER*REGION1:347：9*124：67：9*J4：777：9*KK*AQPROP0.0.3 150.00.0.5不能对不同的区域定义不同的水层性质。

当使用双孔隙模型*DUALPOR时，水层与网格的裂缝孔隙度相连。如果裂缝孔隙不存在，则与基岩相连。子关键字*REGION表示水层边界网网格的直接输入。子关键字*1、*J、*K、*U、*IK.*JK和*IJK表示油藏的那一个面与水层相连。这些方向关键字可以累加。例如：*AQUIFER*REGION111 *1111*J等价于：*AQUIFER*REGIONill/*IJ并且：*AQUIFER*REGION111/*ijk///*I等价于:*AQUIFER*REGION111*ijk注意：*AQUIFER*REGION下指定的网网格或者是油藏边界，或者与水层网格中空网格毗连。活动油藏网格间水的流动在模拟中自动予以考虑。水层仅考虑从外部到油藏的水的流动。*AQPROP指定水层性质，用于计算水侵量。水层压力影响函数(条件)目的：*AQFUNC表示水层累计压力降函数表：AQFUNCdimtimedimpdiop定义：dimtime无量纲时间dimpdrop无量纲水层累计压力降缺省值：如果使用AQUIFER无*AQFUMC,那么使用缺省函数(看如下描述)说明：GEM使用Carter和Tiacy方法计算水侵量：水侵量=a(t)+b(t)压力降函数a和b要求无量纲累计压力降为无量的纲时间的函数。如果水层网格压力低于原始压力，则流量置零。因此，油藏水返流支水层的情况下不予模拟。0.050.100.150.050.100.150.200.300.500.70无量纲时间0.01累计无量纲压力降函数0.1120.22903150.3760.4240.5030.6160.7201.000.8021.50 0.9272.001.0203.00 1.169TOC\o"1-5"\h\z5.00 1.3627.00 1.5001.61510.001.61515.00 1.82920.00 1.96030.00 2.1472.28240.002.28250.00 2.38860.00 2.47670.00 2.55080.00 2.61590.00 2.672100.00200.00300.00区段（选择）格式：SECTOR2.7233.0643.263区名订：12jl：j2kl：k2缺省值：全油田区。如果用双孔模型，则有*MATRIX和*FRACTURE区。说明：一给定区可由几个*SECTOR定义，区的说明是累加性的，一网格可属于不同区。模拟运行结束时每个区的储罐油、气、水体积时在汇总表中找印。它们是通过初始条件部分定义的标准分离器，将区流体闪蒸而得。每个区执行三种流体闪蒸计算，一个为油相，一个为气相，一个为水相。三种结果累加得到区中油、气、水的地面体积。孔隙体积截止阀（选择）（目的：*PVCUTOFF控制的孔隙体积阀值，当网格的孔隙体积小于此阀值，网格就设成死结点。格式：*PVCUTOFFpvcut定义：pvcut孔隙体积（网格体积乘以孔隙度）低于此值，可以认为该网格为无效网格。组分模型GEM组分性质组分性质说明组分类型在模拟中使用桂组分（不包含水）数目由关键字NC定义。一般说來，模拟中可使用两种类型组分：（1） 程序组分：其性质存在程序内部数据库中。GEM内部存贮了油藏模拟中常见的60个组分的全部性质。（2） 用户组分：其性质不存在程序内部数据库中，由用户输入。所有程序组分必须先于用户组分出现。用户组分性质由以下方式输入：SG、*TB、*MW、*AC、*PCRIT.*VCRIT、*TCRIT、*ZCRIT、*PCHOR和*HCFLAG*HCFLAG全部炷的相性质使用状态方程（EOS）计算。状态方程的基本要求是每个组分的临界性质和组分间的相互影响系数（DIJ）对用户组分而言，其临界性质（Pc和Tc）一般是未知的。在这种情况下，其临界性质通过物于是性质（Sg,Tb和Mw）的相关式内部计算得出。相关式由关键字*CRIT定义。相互影响系数用户组分的相互影响系数由用户指定。对炷类组分而言，DU依据其临界体积和PVC3的半经验相关式计算。对所有包含用户组分的其它烧类一非烧类组分的DIJ必须由用户使用*BIN输入。分离器条件第二组状态方程参数为减少相物性拟合的难度，可输入两组状态方程参数。一组用于油藏条件闪蒸计算，另一组用于分离器及标准条件下闪蒸计算。使用下述关键字：BINS,PVC3S,*PCRITS,*OMEGAS,*OMEGBS,*TCRITS和*\公111尸丁$输入标准条件（分离器）计算的状态方程参数。单相一识别当状态方程确定炷流体作为一单相存在时，如果不经过复杂的临界点计算，状态方程并不能确定该相称为'油'或'气然而，相识别在油藏模拟中是重要的,因为油相和气相使用不同的相对渗透率和不同的毛管压力数据。在GEM模拟中,单相的识别由下列4种方法之一确定：（1） 由用户*PHASELDOIL或PHASELD*GAS确定（2） 将*EFDEN所确定的参考密度比较而得：（3） 依据*PHASEUD*CRIT指定的混合物的临界性质近似值确定：（4） 根据相质量密度接近参考气密度或参考油密度來划分相。如果必要，参考密度取自这些关键字只在模拟开始后使用。在油藏初始化中使用的关键字，参见初始条件部分。格式：*MODEL（*PR）（*SRK）定义：*PR用Peng—Robinson状态方程模拟油气相流体性质*SRK使用Soave—Redlich—Kwong状态方程模拟油气相流体性质条件：该关键字必须是组分性质关键字组的第一个关键字。组分数目（必须）格式：*NCnc（nuser）定义：nc烧组分数目（不包括水）nc包括程序组分和用户组分。Nuse［用户组分数目条件：该关键字是组分性质第二个关键字组分名（必须）目的：定义油藏流体组分格式：*COMPNAMEcname定义：cnamenc个组分名，每个最多8个字符。说明：当按顺序排列组分时，所有程序级分排在前面，用户组分排在后面，不允许两者交义。程序组分由其名称莱辨别。如果输入的名称与下述60个组分名称任一个相同，则该组分为程序组分；否则为用户组分。CH4(Cl)FC6FC26C2H6(C2)FC7FC27C3H8(C3)FCSFC28IC4FC9FC29NC4FC1OFC30IC5FC11FC31NC5FC12FC32NC6FC13FC33NC7FC14FC34NC8FC15FC35NC9FC16FC36NC10FC17FC37NC16FC18FC38N2FC19FC39CO2FC20FC40H2SFC21FC41H2OFC22FC42TOLUENEFC23FC43BENZENEFC24FC44CYCLO-C6FC25FC45例如，假设油藏流体由8个组分组成：Cl,C2,C3,NC4,NC5,FRACTION_1,FRACTION_2和FRACTION—3,因此，可将*COMPNAME输为：*NC83*COMPNAME'CH4' 'C2H6''C3H8''NC4''NC5'TRACIJ'FRAC2'<FRAC3,其中前五个组分是程序组分，后三个组分为用户组分。用户只需输入这三个组分性质。在上述例子中，如果希望给出甲烷性质，则组分名可输入如下：*NC84*COMPNAME'C2H6''C3H8''NC4''NC5''METHANE''FRACL''FRAC・''FRAC3'用户给分性质（条件）目的：*SG比重（无量纲）*TB平均常丈夫沸点（degC/degF）*MW分子量（gm/mol）*AC偏心因子（无量纲）*PCRIT临界压力（atm）*VCRIT临界体积（nf/K—MOL）*TCRIT临界温度（degK）*ZCRIT临界压缩因子*PCHOR等张比容（无量纲）*HCFLAG表示用户组分是不否为烧组（整数，0或1）格式：*SGsg*TBtb*MWamolw*ACac*PCRITpc*VCRITvc*TCRITtc*ZCRITZc*PCHORpaiach*HCFLAGhcflag定义：sgjb,amolw,ac,pc,v