B5 · 出版卷 10

水溶液地球化学

Speciation、complexing、pH、Eh、solubility

学习目标

完成本课后，你应能够区分分析总量与溶解物种、浓度与活度、严格定义的 pH 与随意的浓度标签，以及实测氧化还原指标与假定的氧化态图解。你应能计算离子电荷平衡和饱和指数，描述络合及酸碱形态，并建立明确说明组分、相、温度和标准态的反应模型。

本课把水溶液计算视为有条件的热力学和动力学模型。一份溶液分析不会自动给出唯一的自由离子清单、矿物反应或水—岩历史；这些内容需要通过实测总量、温度、必要时的压力、活度模型、平衡常数和约束来推断。只有输入与假设保持可见，模型才可审核。

组分、物种与活度

从水样到组分、物种、活度和矿物饱和度的水溶液地球化学证据链
从水样到组分、物种、活度和矿物饱和度的水溶液地球化学证据链

分析结果常报告 组分 或总量，例如总溶解钙或总硫。化学模型再把这些组分分配到 物种：自由离子、质子化形式、水络合物、离子对、气体和矿物。以溶解无机碳为例，它可存在为溶解二氧化碳、重碳酸根、碳酸根及金属络合物，其比例取决于 pH、温度、离子强度和所选平衡常数。

热力学平衡使用活度。对物种 i，

$a_i=\gamma_i m_i,$

其中 m_i 为质量摩尔浓度，\gamma_i 为声明约定下的活度系数。低离子强度时活度与浓度数值可能接近；盐度较高时差异可能决定结论。活度模型有自己的标定范围，选择模型是科学决策，不是格式选项。

以质量摩尔浓度为基准的离子强度为

$I=\frac{1}{2}\sum_i m_i z_i^2,$

其中 z_i 是电荷。多价离子因电荷平方而贡献很大。只有总量、没有物种和电荷的表格不能直接给出离子强度；需要把质量平衡和平衡关系联立迭代求解。

酸碱平衡、pH 与碱度

pH 由氢离子活度定义：

$\mathrm{pH}=-\log_{10}a_{H^+}.$

它不是某个质量浓度的负对数。测量尺度、缓冲液标定、温度、液接电位和基体都会影响结果。样品若长时间接触空气，测得 pH 可能已不代表原位二氧化碳条件。

碱度 是在明确终点和方法下测得的酸中和容量，以当量每体积或质量表示。它并不等同于重碳酸根，尽管中性附近常由重碳酸根主导。碳酸根、氢氧根、硼酸盐、磷酸盐、有机碱和金属络合物都可能贡献。碱度和酸度约束电荷与质子平衡，常比由模型直接计算的重碳酸根更稳健。

主要离子电荷平衡误差可写为

$\mathrm{CBE}(\%)=100\frac{\sum \mathrm{Cat}-\sum \mathrm{An}}{\sum \mathrm{Cat}+\sum \mathrm{An}},$

两端都必须使用电荷当量，并说明正负号约定。平衡较小不能证明全部分析正确，因为错误可能相互抵消；平衡较大则提示缺失物种、单位、价态、基准或分析问题，必须在解释前查清。

氧化还原、络合与矿物饱和

氧化还原状态由具体电子转移对描述，不是一个万能电极数字。使用 Eh 测量时，应记录电极类型、标定、温度、稳定过程和样品处理。天然水可能没有达到内部氧化还原平衡：铁、硫、碳、氮和砷等氧化还原对会受不同动力学与微环境控制。把 pH 和电位换成单一 pe 并不会迫使所有体系一致。

络合会改变有效迁移能力。若形成稳定配体络合物，即使自由金属离子活度很低，总溶解金属仍可较高。配体可以是氯、碳酸根、硫酸根和氢氧根等无机组分，也可以是有机组分。因此，总浓度趋势可能反映配体供应，而不是物源强度。拟议络合反应所需组分必须实测或合理约束。

矿物溶解反应的离子活度积 IAP 与平衡常数 K 比较：

$SI=\log_{10}\left(\frac{IAP}{K}\right).$

SI<0 表示模型下不饱和，SI=0 表示平衡，SI>0 表示过饱和。过饱和只说明热力学上有沉淀潜力，不证明已经沉淀；成核障碍、抑制剂、表面积、运移和数据库不确定性都可能控制过程。接近平衡也不能唯一识别控制矿物。

采样、保存与模型输入

明确结果属于现场、过滤后、溶解态、可回收总量还是全量，并记录滤膜材质与孔径、容器、保存、保存时间、温度和顶空。过滤可能去除胶体，也可能让更小胶体通过；酸化可能溶解颗粒。“溶解态”是操作定义，不保证所有物质处于分子溶液。

温度、pH、电导率、溶解氧和氧化还原指标等不稳定参数，应按方法要求尽量接近原位测量。记录流量、抽排或稳定过程以及气体交换是否受控。保留碱度滴定原始数据，而不是只存最终值。主要离子、无机碳以及密度或盐度可能是可靠电荷和平衡模型的必要输入。

建模前把单位转换为同一质量或体积基准上的摩尔量或当量，正确处理删失值并核查分析平衡。选择适合温度、压力和组成的热力学数据库与活度模型，记录版本。分析平衡不良时，不能通过调整不确定矿物相来“修好”模型。

合成算例

一份合成水样含下列毫摩尔每升总量：Ca 2.00、Mg 1.00、Na 3.00、K 0.20、重碳酸根 4.00、硫酸根 1.50、氯离子 2.00、硝酸根 0.10。按相应价态，阳离子当量为

$2(2.00)+2(1.00)+1(3.00)+1(0.20)=9.20\ \mathrm{meq/L},$

阴离子当量为

$1(4.00)+2(1.50)+1(2.00)+1(0.10)=9.10\ \mathrm{meq/L}.$

按上述正负号约定，\mathrm{CBE}=100(9.20-9.10)/(9.20+9.10)=0.55\%。这是一项有用的一致性核查，但不证明形态分配或每个测量都正确。

假设 pH 为 7.00，模型给出 \gamma_{H^+}=0.80，则 a_{H^+}=10^{-7}，相应质量摩尔浓度为 1.25\times10^{-7} mol/kg。若某候选矿物反应在模型温度下 IAP=2.0\times10^{-9}、K=1.0\times10^{-8}，

$SI=\log_{10}(0.20)=-0.70,$

因此在所选数据库和活度模型下为不饱和。温度、络合分配或分析输入变化都可能改变结果；饱和指数是计算输出，不是直接观测属性。

解释流程

  1. 定义水体、采样时间、流动状态和操作性组分。
  2. 保存现场条件、制备和不稳定参数测量。
  3. 把总量转换到一致摩尔或当量单位并核查电荷平衡。
  4. 声明组分、相、气体、温压、活度模型和数据库版本。
  5. 在解释自由离子活度和饱和度前先求形态分配。
  6. 比较多个氧化还原对和矿物学观测，不强迫单一平衡状态。
  7. 通过现场梯度和质量平衡检验反应路径、混合与气体交换。
  8. 改变不确定分析输入和热力学选择，评价敏感性。
  9. 分开报告实测量、计算量和地质解释。

常见失败与不确定性

典型错误包括把 mg/L 当作 mol/L，用质量而不是当量核查电荷，把“过滤后”视为真正溶解态，在没有活度约定时把 pH 当作氢离子浓度，以及不核查数据库与温度便接受饱和指数。另一个错误是在多个矿物共享组分时，由接近平衡推断唯一控制矿物。

不确定性来自水体瞬变、气体交换、过滤效应、滴定和标定、删失数据、活度系数、平衡常数以及缺失的动力学约束。多模型集合通常比单一确定性物种表更诚实。

练习与复核

  1. 把一份合成主要离子表从 mg/L 转换为 mmol/L 和 meq/L，并记录式量与价态。
  2. 省略一个二价阴离子后重算电荷平衡，解释诊断特征。
  3. 定性绘制碳酸体系物种随 pH 的变化，并列出为何该图不能直接定量。
  4. 设计区分矿物溶解、蒸发和混合三种总溶解固体上升原因的现场观测。
  5. 比较两种活度模型，说明在哪种组成范围内差异会变得重要。

完成任务要求建立一份合成水—岩反应账本，包含实测总量、现场条件、电荷平衡、候选矿物反应以及低值—中值—高值输入情景。至少报告一个失败模型，并说明失败为何具有诊断意义。

来源与延伸阅读