格式：*AQUIFER（*BOUNDARY） （*BOTTOM）
（*REGION）i1：i2 j1：j2 k1：k2（*I）
（*J） （*K） （*IJ） （*JK） （*IK） （*IJK）
*AQPROP
厚度 孔隙度 渗透率 半径 角度 定义：
*BOUNDARY水层网格均为边缘网格，而非顶部或底部。 （*BOTTOM）水层网格为油藏度部 （*REGION水层网格为一网格区域 *I水层与油藏边界面在I方向 *J水层与油藏边界面在J方向 *K水层与油藏边界面在K方向 *IJ水层与油藏边界面在I和J方向 *JK水层与油藏边界面在J和K方向 *IK水层与油藏边界面在I和K方向 *IJK水层与油藏边界面在I、J、K方向 厚度 水层厚度，缺省值为平均油藏厚度。 孔隙度 水层孔隙度，无缺省值 渗透率 水层渗透率，无缺省值
半径 有效油藏半径，缺省值依据油藏几何因子计算 角度 水层影响角度，表示为一圆周的分数 缺省值 选择关键字。*AQUIFER无缺省值。 *AQPROP的缺省值： 厚度：油藏厚度 孔隙度：无缺省值 渗透率：无缺省值
半径：对*BOUNDARY，为油藏半径。对*BOTTOM，为油藏高度的一半。 角度：以*BOUNDARY而言，为与水层相边的油藏半径分数。 对*BOTTOM而言，为界面长度/（pi。油藏高度）。
说明：在*REGION之后，可以输入多个行，定义复杂的水层。这些独立的行定义同一个水层。 例如：
*AQUIFER *REGION 1：3 4 7：9 *I 2 4：6 7：9 *J 4：7 7 7：9 *KK *AQPROP
0. 0.3 150.0 0.0 .5
不能对不同的区域定义不同的水层性质。
当使用双孔隙模型*DUALPOR时，水层与网格的裂缝孔隙度相连。如果裂缝孔隙不存在，则与基岩相连。
子关键字*REGION表示水层边界网网格的直接输入。子关键字*I、*J、*K、*IJ、*IK、*JK和*IJK表示油藏的那一个面与水层相连。这些方向关键字可以累加。例如：
*AQUIFER *REGION 111 *I
111 *J
等价于：
*AQUIFER *REGION 111/ *IJ
并且：*AQUIFER *REGION 111/ *ijk ///*I
等价于：
*AQUIFER *REGION 111 *ijk
注意：*AQUIFER *REGION 下指定的网网格或者是油藏边界，或者与水层网格中空网格毗连。活动油藏网格间水的流动在模拟中自动予以考虑。水层仅考虑从外部到油藏的水的流动。
*AQPROP指定水层性质，用于计算水侵量。 水层压力影响函数（条件）
目的：*AQFUNC表示水层累计压力降函数 表：AQFUNC
dimtime dimpdrop 定义：
dimtime 无量纲时间
dimpdrop无量纲水层累计压力降
缺省值：如果使用AQUIFER无*AQFUMC，那么使用缺省函数（看如下描述） 说明：GEM使用Carter和Tracy方法计算水侵量： 水侵量=a（t）+b（t）压力降
函数a和b要求无量纲累计压力降为无量的纲时间的函数。
如果水层网格压力低于原始压力，则流量置零。因此，油藏水返流支水层的情况下不予模拟。
GEM中，缺省的无量纲累计压力降函数如下： 无量纲时间 累计无量纲压力降函数 0.112
0.05 0.229 0.10 0.315 0 .15 0.376 0.20 0.424 0.30 0.503 0.50 0.616 0.70 0.720 1.00 0.802 1.50 0.927 2.00 1.020 3.00 1.169
5.00 1.362 1.500
10.00 1.615 15.00 1.829 20.00 1. 960 30.00 2.147 40.00 2.282 50.00 2.388 60.00 2 .476 70.00 2.550 80.00 2.615 90.00 2.672 100.00 2.723 200.00 3.064 300.00 3.263 区段（选择）
格式：SECTOR 区名 i1：i2 j1：j2 k1：k2
缺省值：全油田区。如果用双孔模型，则有*MATRIX和*FRACTURE区。 说明：一给定区可由几个*SECTOR定义，区的说明是累加性的，一网格可属于不同区。
模拟运行结束时每个区的储罐油、气、水体积时在汇总表中找印。它们是通过初始条件部分定义的标准分离器，将区流体闪蒸而得。
每个区执行三种流体闪蒸计算，一个为油相，一个为气相，一个为水相。三种结果累加得到区中油、气、水的地面体积。 孔隙体积截止阀（选择）（
目的：*PVCUTOFF控制的孔隙体积阀值，当网格的孔隙体积小于此阀值，网格就设成死结点。
格式：*PVCUTOFF pvcut 定义：
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定义。 相互影响系数
用户组分的相互影响系数由用户指定。对烃类组分而言，DIJ依据其临界体积和PVC3的半经验相关式计算。
对所有包含用户组分的其它烃类一非烃类组分的DIJ必须由用户使用*BIN输入。 分离器条件第二组状态方程参数 为减少相物性拟合的难度，可输入两组状态方程参数。一组用于油藏条件闪蒸计算，另一组用于分离器及标准条件下闪蒸计算。 使用下述关键字：
BINS，PVC3S，*PCRITS，*OMEGAS，*OMEGBS，*TCRITS和*VSHIFTS 输入标准条件（分离器）计算的状态方程参数。 单相一识别
当状态方程确定烃流体作为一单相存在时，如果不经过复杂的临界点计算，状态方程并不能确定该相称为‘油’或‘气’。然而，相识别在油藏模拟中是重要的，因为油相和气相使用不同的相对渗透率和不同的毛管压力数据。在GEM模拟中，单相的识别由下列4种方法之一确定：
（1） 由用户*PHASEID OIL或PHASEID *GAS确定 （2） 将*EFDEN所确定的参考密度比较而得：
（3） 依据*PHASEUD *CRIT指定的混合物的临界性质近似值确定：
（4） 根据相质量密度接近参考气密度或参考油密度来划分相。如果必要，参考
密度取自这些关键字只在模拟开始后使用。在油藏初始化中使用的关键字，参见初始条件部分。 格式：*MODEL（*PR） （*SRK） 定义：
*PR用Peng—Robinson状态方程模拟油气相流体性质
*SRK使用Soave—Redlich—Kwong状态方程模拟油气相流体性质 条件：该关键字必须是组分性质关键字组的第一个关键字。 组分数目（必须）
格式：*NC nc （nuser） 定义：
nc烃组分数目（不包括水）nc包括程序组分和用户组分。 Nuser用户组分数目
条件：该关键字是组分性质第二个关键字 组分名（必须）
目的：定义油藏流体组分 格式：*COMPNAME cname 定义：
cname nc个组分名，每个最多8个字符。
说明：当按顺序排列组分时，所有程序级分排在前面，用户组分排在后面，不允许两者交叉。程序组分由其名称莱辨别。如果输入的名称与下述60个组分名称
任一个相 同，则该组分为程序组分；否则为用户组分。 CH4（C1） FC6 FC26 C2H6（C2） FC7 FC27 C3H8（C3） FC8 FC28 IC4 FC9 FC29 NC4 FC10 FC30 IC5 FC11 FC31 NC5 FC12 FC32 NC6 FC13 FC33 NC7 FC14 FC34 NC8 FC15 FC35 NC9 FC16 FC36 NC10 FC17 FC37 NC16 FC18 FC38 N2 FC19 FC39 CO2 FC20 FC40 H2S FC21 FC41 H2O FC22 FC42 TOLUENE FC23 FC43 BENZENE FC24 FC44 CYCLO-C6 FC25 FC45 例如，假设油藏流体由8个组分组成：C1，C2，C3，NC4，NC5，FRACTION_1，FRACTION_2和FRACTION—3，因此，可将*COMPNAME输为： *NC 8 3
*COMPNAME ‘CH4’ ‘C2H6’‘C3H8’‘NC4’‘NC5’‘FRAC1’‘FRAC2’‘FRAC3’
其中前五个组分是程序组分，后三个组分为用户组分。用户只需输入这三个组分性质。在上述例子中，如果希望给出甲烷性质，则组分名可输入如下： *NC 8 4
*COMPNAME ‘C2H6’‘C3H8’‘NC4’‘NC5’‘METHANE’‘FRACL’‘FRAC·’‘FRAC3’ 用户给分性质（条件）
目的：*SG比重（无量纲）
*TB平均常丈夫沸点（deg C/deg F） *MW分子量（gm/mol） *AC偏心因子（无量纲） *PCRIT临界压力（atm）
*VCRIT临界体积（m3/K—MOL） *TCRIT临界温度（deg K） *ZCRIT临界压缩因子
*PCHOR等张比容（无量纲）
*HCFLAG表示用户组分是不否为烃组（整数，0或1） 格式： *SG sg *TB tb
*MW amolw