第十章-地球化学热力学数据库
第十章地球化学热力学数据库
(殷辉安，1988；Nordstrom and Munoz, 1994; Gottschalk, 1997)热力学数据内洽性、热力学数据的测量、相平衡数据处理、数据库实例
为了对矿物岩石体系平衡关系作出预测，需要有各地球化学相 (纯矿物、熔体组分、流体)的内洽的热力学数据库。地球化学热力学数据库在过去的几十年间得到了发展。1995年以来，又有6套相对高精度的地球化学热力学数据库(Robie and Hemingway, 1995;
Berman and Aranovich, 1996; Gottschalk, 1997; Chatterjee et al., 1998; Gerya et al., 1998; Holland and Powell, 1998)问世，给地质学研究带来了不少方便。
10.1 热力学数据的内洽性
所谓热力学数据的内洽性(internal consistency),指的是数据之间不存在相互矛盾、热力学数据与已有的实验之间也不存在矛盾的现象。只有内洽的热力学数据，才能客观地反映实际,所以,在研究中,应该从内洽的热力学数据库获取数据。
10.2 热力学数据的测量完整的地球化学热力学数据库,应该包括每一地球化学相(端元矿物、熔体组分、流体等)的标准摩尔生成焓(A f H0)、标准摩尔熵(S0)、标准摩尔体积(V0)、摩尔恒压热容(C p)、等压热膨胀系数(a V)、等温压缩系数(B V)以及电解质的溶解度、电导率、电动势等参数。这些参数的测试方法不完全相同,有量热法、光谱法、电化学测定法等化学上传统的方法,以及新方法例如高温量热、实验相平衡等方法。其中,量热测量(calorimetry)可以获得恒压摩尔热容，相应地可以计算得到标准摩尔生成焓和熵。对于矿物,通常的标准状态选取1 bar、298.15 K。
实验相平衡所获得的矿物热力学数据在目前可能是最可靠的。因为不仅其热力学内洽性可以得到保证,而且其精度可以得到改善。可以说,除量热法之外,实验相平衡研究是减小矿物热力学数据不确定度的主要方法。
对于固体矿物，通过X光衍射，可以获得矿物的晶胞参数，即单位晶胞的轴长a b、c(nm)及轴间角a、B、丫(度)。晶体的单位晶胞体积V(nm3)表达为Vc abc 1 cos2a cos2p cos2丫2cos a cos B cos 丫
矿物的密度P (g/cm3)定义为单位晶胞的质量与体积之比，即
P = MZ/(N A Vc) = MZ/(602.252Vc)
式中M为摩尔质量(g),Z为单位晶胞中的矿物分子数，N A为Avogadro常数(=6.022 m
X1023/mol)，其中M 斗a j n j。m为组分数，n i为矿物化学式中元素i的系数，a i为元
素i的原子量。摩尔体积为V(cm3/mol) = M/ p=VcN A/Z O
10.3相平衡实验数据的处理
由实验相平衡研究获得矿物热力学数据的条件是：(1)原始实验物质的成分和性质已知；(2)反应的平衡关系采用“逆转”实验(reverse experiment)方法圈定；(3)求得之反应方向发生逆转的P-T区间范围(bracket)较窄。
关于实验相平衡数据处理，通常我们面临三个基本问题：(1)确定其热力学内洽性；
(2)求出理想的、符合热力学原理的平衡P-T曲线；(3)导出“未知相”的热力学性
质。
曲线拟合法、线性规划法是目前文献中处理相平衡实验数据的几种常用方法。此处
以线性规划法为例，说明相平衡实验数据的处理。
10.3.1线性规划确定热力学数据的内洽集
(殷辉安，1988)
处理实验相平衡数据的一种较严格的数学解析方法是所谓的线性规划(li near programmi ng)方法。线性规戈U是数学规戈U( mathematical programmi ng)的一种。
Greenwood (1972) Gorden (1973)强调指出，应该从彼此内洽的、代表平衡真实位置极限的实验数据提取有关的热力学参数。这样，采用线性规划法，可以同时达到检验实验数据的热力学内洽性、提取相应热力学参数的目的。
若将矿物化学反应之△ H、△ S作为未知数，则在平衡条件下(A GrO)，可有△H(298.15K) = T A S(298.15K)- A G^
在A H(298.15K)- A S(298.15K)坐标系中，上式是一条以T为斜率、-A G'为截距的直
根据逆转实验(reverse experiment)的特点，对应于每一个实验点(每一个P-T半区间)，要么A H(298.15K) < T A S(298.15K)-A G' 要么A H(298.15K) > T A S(298.15K)-A G'这样，测定热力学常数的问题就变成了求解上述两式定义的一组线性不等式。根据线性规划分析要求，此二式应该写成w或》的形式。对于不等式，其未知数不可能有确定解，但是通常有一组可行解。
A H(298.15K)- A S(298.15K)图中，满足A H(298.15K) < T A S(298.15K) - A G'的可行解在理论上就应当是下图箭头所指示的阴影半区间内的任一点。
满足A H(298.15K) < T A S(298.15K) - A G'的可行解在理论上就应当是图中箭头所指
示的阴影半区间内的任一点。而满足A H(298.15K) > T A S(298.15K)- A G'的可行解就应当
是另一半区间内的任一点
这样，对应于逆转实验的所有P-T 半区间，只要它们在热力学上是彼此内洽的，贝U 在
A H(298.15K) -A S(298.15K)图上就可以得到一个满足所有不等式的凸多边形，见下图 所示。
或者说，只有这一凸多边形区域才是关于这些不等式的唯一的数学上的集。虽然， 从实验角度无法知道平衡的A H(298.15K)、A S(298.15K)值的确切位置，但是它们必然位 于该区域内。倘若在实验上确定的逆转实验点
(P-T 半区间)已经不止一对，却不能得到
这样一个区域，则说明各实验数据彼此不洽合，即从热力学原理看，它们彼此是矛盾的， 其中某些实验数据必然有错误。这就为我们确认实验数据的可靠性提供了依据。
上述凸多边形内洽区可借助线性规划法求得。
所谓线性规划，简言之，就是把一个
称为目标函数的线性函数加上一组约束条件，使之成为极大或极小。这些约束条件可以 包括等式或不等式。计算出目标函数在可行的凸多边形之某一顶点之值，然后沿着此多 边形之棱边依次地求出目标函数在其他顶点处之值， 即可得到同时满足所有不等式之凸
多边形区域。
如下图所示，设A H(298.15K) < T A S(298.15K) -4G'所表达的两个实验点分别为l 、J 。 一般
K)
酒足所有实验虑的二 维凸翕边
潜内梏区城￡图中阴
2 3
番).-巧
而言，在△ H(298.15K)- A S(298.15K)坐标中便有交于0点的两条直线I、J。显然，在0点有：
A S(298.15K) = [ A G' (JG' (I)]/[T(J-)T(I)]
A H(298.15K) = T(l) A S-A G' (I)或A H(298.15K) = T(J) A S-A G' (J)
一般均有上述结果。将那些同所有实验点都洽合的交点都保留下来，即得上述凸多边形内洽区域。
当然，矿物相的任何性质或其结合均可选做目标函数。例如，Gorden (1973,1977) 采用A G和A S做目标函数，而Yin and Greenwood (1983, 1984则取A H、A S做目标函数。
由此可见，线性规划法可以提供一组保证同所有实验数据都洽合的结果。不仅如此，线性规划法可以将量热测量与相平衡实验数据结合一起处理，以导出所研究矿物相的热
力学数据内洽集。但是，这种结果是一个可行解范围，即实质上为无穷多组解，而不是唯一的一组解。因此，不可能以通常概念论及这种解的误差大小或不确定度范围。
已知A H、A S数据后，根据相应矿物的热力学性质便很容易求出未知相的标准生成
焓A f H(298.15K)、标准生成熵A f S(298.15K)。
K)
1032二分逼近法确定平衡P-T曲线
(殷辉安，1988)
应用二分逼近法(interval halving iteration)，可以求出较为理想的、符合热力学原理的平衡P-T曲线。
实例：透闪石热力学分解平衡实验研究
透闪石角闪石族里的一种端元矿物，其低压和高压热力学稳定性曾先后由Boyd (1959)、Gilbert and Troll (1974)、Veblen (1982)、Yin and Greenwood (1983, 1984做过实验研究。
Yin and Greenwood (1984以氧化物混合物、凝胶混合物作为原始试料，合成透闪石。实验研究的反应是
Ca2Mg5Si8O22(OH)2 = 3 MgSiO3 + 2CaMgSi2O6 + SiO2 + H2O
tremolite enstatite diopside quartz vapor
透闪石顽辉石透辉石石英水实验采用采用冷封式和内热式高压弹做逆转反应实验，获得如下表的数据
注：括弧内数字为2 c值，且指最末位数字；De代表组合3En + 2 Di + Qtz + H20;
***以Cal2（0.08 M）、Mgl2（0.2M）溶液代替水注入样品囊管；****在Tr-D-2-9的基础上再反应14天