GHG选项 – 水相模型（可选）
Aqueous Phase Model (Optional)
目的
使用内部创建的包括矿物成分的湮灭矩阵来计算减少组分的物质平衡误差。
格式
| *ACTIVITY-MODEL | *IDEAL | *DEBYE-HUCKEL | *B-DOT *PITZER | *PITZER2 |
| *ACTIVCOEF-SOL | *ON | *OFF |
| *ACTIVCOEF | |
| *SALINITY-CALC | *ON | *OFF |
| *AQFILL | *ON | *OFF |
定义
*ACTIVITY-MODEL
水相的活性模型。其后必须跟随IDEAL、DEBYE-HUCKEL、B-DOT、PITZER或PITZER2之一。对于IDEAL，组分活性系数等于1。
对于其他情况，活性系数通过Debye-Hückel模型、B-dot模型或Pitzer模型计算得出。PITZER是PITZER2的简化版本，运行速度更快。PITZER2是更严格、完整的Pitzer版本。
默认情况下使用IDEAL。有关Pitzer模型详细信息，请参见解释部分。
*ACTIVCOEF-SOL
在2013.12版之前的GEM版本中，可溶性组分（如CO2(aq)）的活性系数设为1。因此，可溶性组分的活性等于它们的摩尔法。而在2013.12版中，使用ACTIVCOEF-SOL ON将可溶性组分的活性系数设为与盐分系数相等。ACTIVCOEF-SOL的默认值为OFF，与GEM的之前版本保持一致。
*ACTIVCOEF
此关键词允许输入水相中所有组分、nc油气组分、H2O和nc_aqu水相组分的活性系数。因此需要nc + 1 + nc_aqu个实值。对于不可溶于水相的组分，请输入1。这些活性系数将代替Debye-Hückel模型或B-dot模型。请注意，除H2O外的所有组分的活性等于活性系数与摩尔法的乘积。H2O的活性等于活性系数与H2O的摩尔分数的乘积。
*SALINITY-CALC
盐度的计算。其后必须跟随ON或OFF。如果输入的是OFF，则假定水相的盐度恒定，并可使用关键词SALINITY指定。如果输入的是ON，则默认从水相中的Cl-浓度计算盐度。如果未定义Cl-，则使用Na+。如果使用COMPNAME-SAL指定了盐度组分，则盐度计算中将使用盐度组分的浓度。
*AQFILL
此关键词指定含水层完全充满水。其后必须跟随ON或OFF。如果未使用此关键词或输入的是OFF，则含水层中始终存在少量烃相痕迹。使用ON，水相的饱和度等于1。如果使用此关键词，溶于水相的气组分的最大数量必须为1。
解释
Pitzer方程是一组用于建模电解质活性的方程。Pitzer活性模型基于过量吉布斯能的维里尔展开（Pitzer，1973年和Pitzer，1974年）。Pitzer方程的应用范围包括高离子强度溶液，特别是高NaCl浓度的溶液。由于Pitzer活性模型的计算比其他活性模型更费力，所以计算成本更高。在大规模计算中选择适当的活性模型时要小心。
Pitzer方程是用来描述电解质溶液中活度系数的一种模型。这个方程特别适用于高离子强度条件下的溶液，例如盐湖或者海水。
假设我们有一个含有钠离子（Na⁺）和氯离子（Cl⁻）的盐水溶液。在这个溶液中，离子之间会发生相互作用，导致它们的活度系数与理想行为有所不同。
使用Pitzer方程，我们可以更准确地描述这种相互作用。例如，我们可以通过Pitzer方程计算钠离子和氯离子的活度系数，以更好地预测它们在这种高离子强度环境中的行为。
总的来说，Pitzer方程提供了一种对非理想电解质溶液中活度系数进行模拟的方法，通过这种方式可以更精确地描述溶液的化学特性。
