*AC ac *PCRIT pc *VCRIT vc *TCRIT tc *ZCRIT Zc
*PCHOR parach *HCFLAG hcflag 定义：
sg,tb ,amolw,ac,pc,vc,tc,zc,parach 和bcflag全部由nuser个实数组成。 置HCFLAG
1如果用户组分为烃组分 0如果用户组分为非烃组分
说明：这些关键字用于输入用户组分性质。如果其中某一个未知数，输入零值。如果对所有用户组分该性质未知，则跳过该关键字。并非所有性质必须输入。 如果使用Jossi—Stiel—Thodos相关式，计算粘度时需要vc参数。Zc可代表vc因为vc可由zc直接计算而得。如果*ZCRIT和*VCRIT都缺省，则ZC由偏心因子计算而得。
状态方程流体性质计算需要所有组分的tc。Pc，ac和amolw。如果其中某一性质未知，则输零值或缺省，它们由其它可用数据依据经验公式内部计算而得。使用*CRIT关键字选择所需的经验公式。
Sg和tb仅用于计算所需的临界性质和amolw。如果所有必须的tc。Pc，ac和amolw已积压，则不需要*SG和*TB。
如果要求sg或tb，并且它们未知，由使用与sg，tb和amolw有并的经验公式内部计算而得。因此，对sg，tb，pc，tc，ac，和amolw来说，至少需要sg，tb或amolw之一。
一般来说，给出的参数越多，相态预测越准确。
等张经容用于计算界面张力。如果未知，则由分子量计算。 Hcflag指明用哀悼组分是否为烃组分。烃组分间二元相互影响系数由内部公式计算（参见*PVC3和*PV3S）；然而，对烃组分与非烃组分之间，非烃组分与非烃组分之间的相互影响系数使用BIN和BINS输入。因此，在GEM中，一组分分为烃组分，可使模拟器或者由输入值，或由公式得到二元相互影响系数。如果用户希望GEM使有特定二元相互影响系数值，则其相应组分必须标定为非烃组分。否则，二元相互影响系数值依计算而得，不使用读入值。 状态方程OMEGA参数（可选） 目的：*OMEGA状态方程Ωa 状态方程Ωb
格式*OMEGA omega *OMEGB omegb 定义：
omega nc个实数，表示流体Ω a参数： omegb nc个实数，表示流体的Ωb参数 缺省值：对Peng—Robinson状态方程为： Ω a=4.57235530e-01 Ωb=8.6640350e-01
说明：*OMEGAS和*OMEGBS用于计算分离器和地面条件流体性质。缺省值为*OMEGA和*OMEGB。
用户组分二元相互影响系数（条件） 格式：*BIN delij 定义：
delij nuser*（2*nc—nuser—1）/2个实数 缺省值：零
条件：仅当有用户组分时需要
说明：*BINS用于计算分离器和地面条件的流体性质。 缺省为*BIN
烃类相互影响参数（选择） 格式：*PVC3、PVC3 定义：
PVC3定义用于计算烃组分间二元相互影响系数的实数 Delij定义如下：
|2*sqrt（vc3（I）*vc3（j）|**PVC3 Delij=1-|-----------------------------------|
| vc3（I）+vc3（j） |
其中：VC3I=（组分I的临界体积）1/3
pvc3=0，表示零二元相互影响系数。PVC3典型值为1.2 缺省值为0。
说明：*PVC3S用于计算分离器和地面条件流体性质。缺省为*PVC3。 用户组分临界性质关联（条件） 格式：*CRIT （*LEEKESLER） （*RIAZI） （*Twu） 定义：
*LEEKESLER Lee—kesler公式 *RIAZI Riazi—Daubert公式 *Twu Twu化工
缺省值：* LEEKESLER
说明：仅当用户组分存在，而部分临界性质未知时使用。 体积平移参数（选择）
目的：*VSHIFT用于指定体积平移量，计算状态方程摩尔体积。 格式：*VSHIFT Shv 缺省值：零。
说明：在GEM中，使用peneloux等的体积平移方法，改善密度预测。计算的相压缩因子（z因子）和摩尔体积V，修正如下： Z=Zo-P（V~）/（RT） V=Vo-V~
V~=Sun[x（I）*Ti（I）] T（I）=shv（I）*bc（I）
Bc（I）=Omegb（I）*R*Tc（I）/pc（I） 其中：Sun=所有nc个组分求和
zo=状态方程直接计算的压缩因子
Vo=状态方程直接计算的摩尔体积，m3/kmol Nc=总组分数
X（I）=组分I的摩尔分数
Omegb（I）=状态方程Ωb参数 对PR EOS，Ωb=0.07780 对SRK EOS，Ωb=0.08664
对小于C6的组分，shv〈0；对大于C6的组分，shv〉0
VSHIFTS指定体积平移量，用于计算分离器和地面条件流体性质。如果VSHIFTIS缺省，则使用*VSHIFT值。 地面条件状态方程参数（选择） 格式：*BINS delijs *PVC3S pvc3s *PCRITS pcs *OMEGAS omegas *OMEGBS omegbs *TCRITS tcs *VSHIFTS shvs
说明：当用少量组分代表油藏流体时，想用一组EOS参数同时精确到地面条件和油藏条件相态是因难的。用两组参数可减少组分数量。 油藏温度（必须） 格式：*TRES tres 或
*TRES depthl tresl depth2 tres2 定义：
tres 油藏温度 depthl深度1
tresl在深度1处的温度 depth2深度2
tres2在深度2处的温度
条件：如果TRES后跟一单实数，则油藏中所有位置温度定为该值。如果TRES后跟 4个实数，则油藏中网格的温度取决于网格中心的深度。
T（网格）= tresl+[（tres2—tresl）/（depth2—depthl）*[depth（网格）—depthl] 在任一给定深度，温度在整个模拟中取常数。 平均饱和压力（选择） 格式：*PSAT psat 定义：
psat 原始油藏流体饱和压力 缺省值：1.0e+20
说明：当网格压力大于psat的绝对值时，流体处于单相状态。这可大大降低闪蒸计算的次数。为避免不一致性，psat必须大于或等于模拟器中使用EOS计算的饱和压力。该值可使用CMG相态软件包CMGProp计算而得。
如果psat为正数，网格流体组成原油藏流体组成相同，则跳过闪蒸计算。当网格压力大于psat的绝对值时，用户希望跳过闪蒸计算，而不管网格流体组成，
则将psat输入为负数。
因此，如果PSAT=-1.0，则整个油藏始终保持单相状态。这用于一次接触混相过程。
单相流体的识别（可选） 格式：*PHASEID（*OIL） （*GAS） （*CRIT） （*DEN） 或
*REFDEN dvref 定义：
*OIL表示模拟中所有单相网格为油相网格 *GAS表示模拟中所有单相网格为气相网格 *CRIT单相识别临依据下列准则确定： a）在超临界条件，称为气。
b）在次临界条件，当摩尔体符号小于临界摩尔体积时，称为气
为避免复杂的临界点计算，假定依据混合规则计算的混合物EOS临界性质为实际临界值。
*DEN依据质量密度更接近于参考气相密度或参考液相密度分类单相。参考密度自动计算，必要时使用稳定性测试中的不稳定相，不必用户输入。 Dvref参考质量密度（kg/m3|lb/ft3）
缺省值：*REFDEN不能和*PHASEID同时输入。 说明：有四种方法定义可知块中单相流体：
（1） 直接由*PHASEID OIL或PHASEID *GAS定义。
（2） 将其质量密度与*PEFDEN输入的参考质量密度较，当质量密度
〈dvref-delt，则假设流体为气相；当质量密度〉dvref+delta时，则假设流体为油相，其中delta=3.12psi(50kpa)
当dvref-delta〈质量密度〈dvref+delta时，如果格块中气相烃摩尔分数在前一时间步结束时（对第一个时间步而言，则为初始化结果时）超过0.5，则流体为气相。在初始化时，假定油区和水区气相烃摩尔分数为零（水油界面和气油界面相等的气藏除外），气顶中的气相烃摩尔分数由给定的烃组成闪蒸确定。
（3） Gosset方法：用混合物近似临界性质确定相一致性。
（4） 质量密度方法：内部计算的参考气相和参考油相质量密度判定。
当同时使用*PHASEID和*PEFDEN时，密度准则覆盖相识别设置*PHASEID；当对两系统，低质量密度相总是假定为气。 粘度关联定义（可选） 格式：*VISCOR（*MIX） （*GZYT） （*LE） 定义：
**MIX依据用户输入的油藏条件组分粘度，直接使用混合规则计算油藏条件烃相粘度
*HZYT YOON—Thodos公式：
mu（low，I）=4,610*tR0.618-2.040*exp（-0.449*Tr） +1.94*exp(--4.058*Tr)+0.1)*10**(-4)/mup 其中Tr=Tabs/Tc
331/6
Mup=[Tc/（MP4）]Tc—单位为K： Pc—为大气压： M为g/g mol；
Mu（low，I）—单位为厘泊。 Hernirg—Zipperer公式：
Sum mu（low，I）*x（I）*sqrt（M（I） I=1，nc
Mu（mix）=------------------------------------------------------- Sum mux（I）*sqrt（M（I） I=1，nc
Mu（mix）=（Σi=1ncmu（low，I）*Xi*sqrt（Mi））/（Σi=1nc Xi*sqrt（Mi）） 依据Hernirg-Zipperer公式，与Yoon_thodos公式计算的低压粘度一起，计算烃混合物的低压粘度。油藏条件混合物粘度使用Jossi，Stiel thosos公式由低压粘度计算而得。
{[ Mu（mix）—Mu（low，mix）]。Mup（mix）+10-4}1/4 =0.10230+0.23364.rr+0.58533.rr2-0.040758.rr3+0.0093324.rr4 其中：
rr=临界摩尔体积/摩尔体积
*LE混合物低粘度是绝对温度和平均分子量的函数。 Lee-Eakin公式如下：
Mu（mix）=[10-4*（17.94+0.321M）*T3/2]/(1.8T+75.4+13.9M) T是绝对温度，M为混合物平均分子量g/gmol。
缺省值：如果缺省VISCOR，则假定*VISCOR HCYT。如果在数据中说明VISCOSITY，则缺省值改变为MIX。如果在*TRES关键字下输入油藏非零温度梯度，缺省值改为*LE，对整个油藏计算与温度相关的粘度。
说明：由于计算油藏中不同温度下*HZYT公式，需要对油藏中每个组分及每个网格块计算几次指数函数，当油藏温度梯度非零时，如果说明*HZYT，则不计算与温度有关的烃粘度。计算油藏中部深度的组分粘度，并且在模拟期间不改变。 当有一非零温度梯度时，粘度公式缺省为*LE。 计算烃流体粘度参数（选择） 目的：
*VISCOSITY表示用直接混合规则计算烃粘度时使用的组分粘度。 MIXINGRULE表示直接混合规则计算烃粘度时的指数。
*VISC表示Jossi，stiel Tjodos关联式中用于计算混合物临界体积的组分拟临界体积。
*MIXVC表示Jossi，stiel Tjodos关联式中用于计算混合物临界体积的指数。 *MVISCOEFF表示Jossi，stiel Tjodos关联式系数。 格式：*VISCOSITY amus MIXINGRU*LE emixmu 或