碳酸平衡

2.3 碳酸平衡

地下水中的碳酸盐组分，大气中的CO2，岩石中的碳酸盐共同组成了一个完整的碳酸平衡体系，它们之间的化学反应对地下水化学成分的形成与演化起着重要的控制作用，是理解地下水系统中许多地球化学过程和现象的基础。因此在天然水化学或水文地球化学研究中，都离不开对碳酸平衡深入讨论和研究。

2.3.1 CO2－H2O体系

1. CO2分压（Pco2）已知

CO2在水中的溶解可用下述的化学反应来表达：

CO2 + H2O = H2CO3 （2-3-1）

溶液中溶解CO2的浓度取决于与其平衡的大气CO2的分压Pco2，可根据Henry定律由下式计算：

[H2CO3] = KH×Pco2 （2-3-2）

式中，KH为Henry常数，它是温度的函数，随着温度的升高，KH值减小，表2-3-1给出了不同温度下的KH值，25℃时KH =10

-1.47

。在表2-3-1中，平衡常数K的数值是以pK的形式给出

的，与pH值的定义类似，pK= -lgK。

除了CO2的溶解反应外，在CO2－H2O体系中的重要反应还包括H2CO3的一级和二级电离反应以及水的离解反应：

+

-

H2CO3 = H + HCO3 （2-3-3）

-+

2-

HCO3 = H + CO3 （2-3-4）

+

-

H2O = H+ OH （2-3-5）

这些反应的平衡常数可依次表示为：

（25℃） （2-3-6）

（25℃） （2-3-7）

（25℃） （2-3-8）

表2-3-1也给出了不同温度下的Ka1、Ka1、KW的值，同时给出的还有方解石的溶度积Ks0。

由式（2-3-2）、（2-3-6~8）有：

lg [H2CO3] =lg KH + lg Pco2

= -1.47 +lg Pco2

lg [HCO3] =lg Ka1 +lg [H2CO3] –lg [H]

=lg Ka1 +lg KH +lg Pco2 + pH = -7.82 +lg Pco2 + pH

lg [CO3] =lg Ka2 +lg [HCO3] –lg [H]

=lg Ka2 +lg Ka1 +lg KH +lg Pco2 + 2pH = -18.15 +lg Pco2 +2pH

+2--+

-+

因为 lg [H] = -pH

故由式（2-3-7）有：

-

lg [ OH ] = -14 + pH

表2-3-1 CO2－H2O体系中的主要化学反应的平衡常数随温度的变化

t（℃）

0 5 10 15 20 25 30 35 40 45 50 100 150 200

pKH 1.11 1.19 1.27 1.33 1.41 1.47 1.53 1.59 1.64 1.68 1.72 1.99 2.07 2.05

pKa1 6.579 6.517 6.464 6.419 6.381 6.352 6.327 6.309 6.298 6.290 6.285 6.45 6.73 7.08

pKa2 10.625 10.557 10490 10.430 10.377 10.329 10.290 10.250 10.220 10.195 10.172 10.16 10.33 10.71

pKW 14.955 14.734 14.534 14.337 14.161 13.999 13.833 13.676 13.533 13.394 13.263 12.27 11.64 11.28

pKcal 8.03 8.09 8.15 8.22 8.28 8.34 8.40 8.46 8.51 8.56 8.62 9.62 10.54 11.62

（据Butler, James N, 1982）

一般情况下，大气中的CO2分压（Pco2）为一常数，等于10组分浓度随pH值的变化关系曲线（图2-1-1）。

-3.5

atm，据此可绘制上述各

在该系统中，水溶液的电荷平衡方程为：

+

--2-

[H] = [ OH ] + [HCO3] + 2[CO3] （2-3-8）

-

在所限定的条件下，水溶液的pH值不会超过7，水中带负电荷的主要组分显然为HCO3，与[HCO3]相比，[OH]与 2[CO3]被忽略，这一点也可以由图2-1-1看出来。故有：

+

---2-

[H] ≈ [HCO3] （2-3-9）

即

[H] = Ka1×KH×Pco2 =10

+

-5.65

+2-11.3

[H] =10

pH = 5.65

这就是雨水的pH通常低于6.5的原因所在。

图2-1-1 Pco2恒定时CO2－H2O体系中各组分含量随pH值的变化

2. 水中含碳组分总量CT已知

水中碳酸盐组分的总量可表示为：

CT = [H2CO3] + [HCO3-] + [CO32-] （2-3-10）

由式（2-3-6）、（2-3-7）有：

（2-3-11）

（2-3-12）

以式（2-3-11）、（2-3-12）代入式（2-3-10）得：

（2-3-13）

由式（2-3-13）有：

（2-3-14）

代（2-3-14）入式（2-3-11） 、（2-3-12）得：

（2-3-15）

（2-3-16）

若以α0、α1、α2分别表示[H2CO3] 、[HCO3] 、[CO3]在碳酸盐组分总量CT中所占的比例，即：

-2-

（2-3-17）

（2-3-18）

（2-3-19）

则有：

（2-3-20）

（2-3-21）

（2-3-22）

根据式（2-3-20~22）所绘制的水溶液中各碳酸盐组分在CT中所占比例随pH值的变化关系如图2-3-2。

图2-3-2 α0、α1和α2随pH值的变化关系

-2-

同时，根据式（2-3-20~22）也可对H2CO3 、HCO3 、CO3含量随pH值的变化进行定性分析。

-

当[H2CO3]=[HCO3]时，α0=α1，故有：

[H] = Ka1 = 10

+

-6.35

pH=6.35

-

考虑到水溶液中H2CO3与HCO3之间的化学平衡关系（式2-3-3），可知当pH<6.35时，[H2CO3]>[HCO3]；当pH>6.35时，[H2CO3]<[HCO3]。

-2---

当[HCO3]=[CO3]时，α1=α2，故：

[H] = Ka2 = 10

+

-10.33

pH = 10.33

-2--2-

考虑到HCO3与CO3之间的化学平衡关系（式2-3-4），当pH<10.33时，[HCO3]>[ CO3]；当pH>10.33时，[HCO3]<[ CO3]。

2--2-

当[H2CO3]=[CO3]时，α0=α2，故有：

[H] = Ka2Ka1 = 10

+2

-16.68

pH = 8.34

2-2-

显见，当pH<8.34时，[H2CO3]>[CO3]；当pH>8.34时，[H2CO3]<[CO3]。

综合上述分析可知，当pH<6.35时，H2CO3在各碳酸盐组分中的含量最大；当pH>10.33时，CO3在各碳酸盐组分中的含量最大；而当6.35

--2--

即： -[H]+ Ka2Ka1 = 0

[H] = Ka2Ka1 = 10

+2

-16.68

+2

pH = 8.34

-2-

故pH = 8.34时，水溶液中的HCO3浓度达到了最大，同时在该pH点上，[H2CO3]=[CO3]。

2.3.2 碳酸盐的溶解平衡

在地球壳层可溶于地下水的矿物中，碳酸盐可能是丰度最大且最重要的矿物。其中最常见的就是CaCO3，它是地球上石灰岩沉积物及岩石的主要成分。自然界中 CaCO3有两种矿

物，一种是方解石，另一种是文石。在海相沉积岩中，白云石（CaMg（CO3）2）及镁质方解石（

<1）也是常见的碳酸盐矿物。

25℃时，方解石溶解反应的平衡常数为：

Kcal=[Ca2+][CO32-]=10-8.37 （2-3-23）

文石的平衡常数与此类似，25℃时Karag=10

-8.16

。但由于内部结构不同，两者还是有一定程度

的差别。Karag> Kcal表明方解石比文石更稳定，并且会首先从溶液中沉淀出来。

白云石的溶解反应可表示如下：

2+

2+

2-

CaMg（CO3）2 = Ca＋ Mg ＋ 2 CO3 （2-3-24）

因此其溶解反应的平衡常数可表示为：

Kcal=[ Mg2+][Ca2+][CO32-]2=10-16.79 （25℃） （2-3-25）

在地下水中，白云石的溶解速度通常要比方解石慢，明确这一点，对于解释地下水中的一些溶解沉淀作用是非常重要的。

在纯水中，碳酸盐矿物是微溶的，但在CO2的影响下，碳酸盐矿物的溶解度可大大增加。由于系统相对CO2的性质不同，水溶液中CO2的溶解量差别很大，因此系统的性质是决定碳酸盐矿物溶解度的另一个重要因素。在封闭系统中，CO2在反应过程中会被不断地消耗，直至消耗殆尽，因此相对于能够获得CO2持续补充的开放系统来说，碳酸盐在其中的溶解度要小一些。下面方解石为例，分别对开放和封闭系统中碳酸盐的溶解平衡予以讨论。

1. CaCO3－H2O系统

+

--

在该系统中，需要对其浓度进行计算的组分有6种，分别为：H、OH、H2CO3、HCO3、CO3、Ca，因此求解该问题需要有6个方程，分别为：

2-2+

（1） 化学计量关系方程：

[Ca] = [H2CO3] + [HCO3] + [CO3] （2-3-26）

2+

2+-2-

这是由系统中的Ca和碳酸盐组分均来自于方解石的溶解所决定的。

（2） 电中性方程：

+

2+

--2-

[H] +2[Ca] = [OH] + [HCO3] + 2[CO3] （2-3-27）

（3） H2CO3的一级电离方程（式2-3-6）；

（3） H2CO3的二级电离方程（式2-3-7）；

（5） 水的电离平衡方程（式2-3-8）；

（6） 方解石的溶解平衡方程（式2-3-23）。

考虑到该系统中水溶液的pH>7，舍去式（2-3-26） 和（2-3-27）中含量较小的各项有：

2+

-2-

[Ca] = [HCO3] + [CO3] （2-3-28）

2+

--2-

2[Ca] = [ OH ] + [HCO3] + 2[CO3] （2-3-29）

（2-3-28）×2-（2-3-29）得：

--

[HCO3] = [ OH ] （2-3-30）

代（2-3-30）入式（2-3-28）得：

2+

-2-

[Ca] = [OH] + [CO3] （2-3-31）

由式（2-3-8）有：

（2-3-32）

由式（2-3-8）有：

（2-3-33）

以式（2-3-30）、（2-3-32）代入（2-3-33）得：

（2-3-34）

由式（2-3-23）并结合式（2-3-34）有：

（2-3-35）

把式（2-3-35）、（2-3-32）、（2-3-34）代入式（2-3-31）并整理得：

（2-3-36）

可使用迭代法或图解法求解方程（2-3-36）。使用迭代法求解时，给定氢离子浓度的初值为10，通过五次迭代（10解为[H]=10[OH]=10

-+

-9.937

-9

-9.980

、10

-9.933

、10

-9.938

、10

-9.937

、10

-9.937

）便可求得方程（2-3-36）的

2--4.456

。这样水溶液中其他组分的浓度依次为（单位为mol/L）：[CO3]=10

--4.063

，

-4.063

，[HCO3]=10，[H2CO3]= 10

-7.650

，[Ca]=10

2+-3.914

。据此求得方解石在纯水中的

溶解度为12.19 mg/L。这与文献（袁道先等，1988；中国科学院地质研究所岩溶研究组，1979）中报道的14 mg/L的方解石的溶解度相比有一定程度偏小，主要是由于未考虑各组分活度系数的影响而造成的。

2. CaCO3－H2O－CO2系统

在这种情况下所研究的问题实际上就是在一定PCO2下，方解石在纯水中的溶解问题。即CaCO3－H2O－CO2系统中水－岩－气的平衡问题。在该系统中，需要对其浓度进行计算的组分有9种，分别为：H、OH、H2CO3、HCO3-、CO3、Ca、CaHCO3、CaOH、CaCO3。因此求解该问题需要有9个方程，其中8个方程为质量作用方程，它们依次是：

+

-2-2+

+

+

0

（1） 水的电离平衡方程（式2-3-8）；

（2） CO2的溶解平衡方程（式2-3-2）；

（3） H2CO3的一级平衡方程（式2-3-6）；

（4） H2CO3的二级平衡方程（式2-3-7）；

（5） CaCO3的溶解平衡方程（式2-3-23）；

（6） 形成CaHCO+

3的络合平衡方程：

Ca2+

＋ HCO-+

3 =CaHCO3

（25℃） （7） 形成CaOH+

的络合平衡方程：

Ca2+

＋ OH- =CaOH+

（25℃） （8） 形成CaCO0

3的络合平衡方程：

Ca2+

＋ CO2-0

3 =CaCO3

K2+2-3=[Ca][CO3]=103.1 （25℃） 另一个方程就是下述的电荷平衡方程：

[H+

]+2[Ca2+

]+[CaOH+

]+[ CaHCO+

--] + 2[CO2-

3] = [ OH ] + [HCO33] 这样，系统中各组分的浓度可用氢离子浓度（pH）和Pco2表示如下：

lg [H2CO3] = -1.47 +lg Pco2

（2-3-37）（2-3-38）（2-3-38）2-3-39）

（

lg [HCO3] = -7.82 +lg Pco2 + pH lg [CO3] = -18.15 +lg Pco2 + 2pH

lg [Ca]=-8.37-lg [CO3]=9.78-lg Pco2-2pH

+

-2+

2+

2-2-

-

lg[CaHCO3]=1.22+lg[HCO3]+lg [Ca]=3.18-pH

0

2-2+

lg[CaCO3]=3.1+lg[CO3]+lg [Ca]=-5.27

lg[CaOH]=1.3+lg[ OH ]+lg [Ca]=-2.92-lg Pco2-pH

+

--2+

lg [H] = -pH

-

lg [ OH ] = -14 + pH

根据上述公式可绘制出不同Pco2条件下各组分的浓度随水溶液pH值的变化关系，图2-3-3和图2-3-4分别为Pco2=10

-3.5

atm和Pco2=1 atm时的情况。

图2-3-3 Pco2=10

-3.5

atm 时CaCO3－H2O－CO2体系中各组分浓度随pH值的变化

由图2-3-3可见，当Pco2=10

-3.5

atm时，在水溶液中的所有阳离子中，Ca的含量明显

2+

最大；而在所有的阴离子中，HCO3的含量显著大于其他离子。因此在这种情况下水溶液的电荷平衡方程（2-3-39）可简化为：

2+

-

-

2[Ca] ≈ [HCO3] （2-3-40）

即

[H]=10

+3-7.82

×10

-0.3

×10

-9.78

×Pco2=10

2-17.9

×Pco2

2

（2-3-41）

与Pco2=10

-3.5

atm相比，Pco2=1 atm 时（图2-3-4）溶液中CaHCO3的含量明显增大，

+

这时水溶液的电荷平衡方程（2-3-39）可简化为：

2+

+

-

2[Ca]+[CaHCO3] ≈ [HCO3] （2-3-42）

即

[H]=10

+3-7.82

×10

-0.3

×10

-9.78

×Pco2=10

2-17.9

×Pco2

2

（2-3-43）

图2-3-4 Pco2=1 atm时CaCO3－H2O－CO2体系中各组分浓度随pH值的变化 令Pco2=1 atm，使用迭代法或图解法求解方程（2-3-43）可得[H]=10（2-3-41）可求得[H]=10

+

-5.97

++

-5.98

；直接使用式

。可见，尽管Pco2=1 atm 时溶液中CaHCO3的含量有明显增大，

但其对溶液pH值的计算结果并无太大影响。在Pco2≤1 atm的情况下，完全可以使用式（2-3-41）计算CaCO3－H2O－CO2系统中水溶液的pH值。由式（2-3-41）所反映的水溶液pH值随CO2分压的变化关系见图2-3-5。

这样，根据所求得的水溶液的pH值以及所给定的CO2的分压，便可方便地求出水溶液中各组分的浓度。例如当Pco2=1 atm时，水溶液中各组分的浓度分别为（单位为mol/L）：[H2CO3] = 10[CaCO3]=10

0

-1.47

，[HCO3] =10

-

--1.84

，[CO3] = 10

+

-5.98

2--6.19

，[Ca]=10

-8.02

2+-2.18

，[CaHCO3]=10

+-2.8

，

-5.27

，[CaOH]=10

-8.9

，[H] =10，[OH] =10

2+

-

。

在CaCO3－H2O－CO2系统中，水溶液中的平衡Ca浓度反映了水中方解石的溶解度。根据式 （2-3-6）、（2-3-7）有：

由式（2-3-2）、（2-3-23）有：

即 （2-3-44）

以式（2-3-40）代入式（2-3-44）得：

（2-3-45）

把25℃时的平衡常数值代入式（2-3-45）有：

（2-3-46）

根据式（2-3-46）作图，可得如图2-3-6所示的曲线。

图2-3-6 混合溶蚀作用图解

在图2-3-6中，位于曲线之上的部分为方解石的过饱和区，位于曲线之下的部分为方解石的未饱和区，而处于曲线之上的点则正好与方解石处于平衡状态。根据该曲线可以发现一个有趣的现象，这就是当两种与方解石处于平衡状态的水混合时，混合水相对方解石是处于未饱和状态的。例如当图中处于A点和B点的水发生混合时，尽管A、B两点的水与方解石正好处于平衡状态，但位于A、B之间直线上的混合水则相对方解石是处于未饱和状态的，

这就是混合溶蚀作用依据所在。关于混合溶蚀作用更详细的讨论，读者可参阅有关的文献（Bogli，1964；钱会，1996）。

2.3.3 其他离子的影响

-2+

当通过在水溶液中加入相关易溶盐类而使溶液中的HCO3或Ca浓度增加时，溶液中与方解石平衡的Ca浓度也会相应改变，这种现象被称为同离子效应（Common Ion Effect）。例如，当在水溶液中加入NaHCO3时，电荷平衡方程（2-3-40）变为：

+

2+

-2+

[Na] + 2[Ca] ≈ [HCO3] （2-3-47）

-

把式（2-3-47）中的[HCO3]代入式（2-3-44）得：

（2-3-48）

据此可求得加入不同量的NaHCO3时（由溶液中的Na浓度所反映），溶液中与方解石平衡的Ca浓度随Pco2的变化关系曲线（图2-3-7）。由图可见，随着NaHCO3加入量的增加，与方解石平衡的Ca浓度不断减小。因此，NaHCO3的加入使得方解石在水溶液中的溶解度不断降低。

2+

2+

+

图2-3-7 NaHCO3的加入对方解石溶解度的影响

2+

2-

-

在与方解石平衡的地下水中加入石膏时的情况与此类似。随着石膏的溶解，Ca与SO4

进入了水中，与上述情况不同的是，这时的同离子不是HCO3而是Ca。在这种情况下，Ca

--2+

2+

的加入使方解石过饱和并发生沉淀，从而使水溶液中的HCO3含量减少。根据Appelo和Postma（1993）的研究，因同离子效应而发生沉淀的前提条件是，两种矿物的溶解度相差在3个数量级以上。对于石膏和方解石，这一条件只有在温度大于25℃时才满足。

除了上述的同离子效应外，另外一种常见的情况就是向水溶液中加入不含同离子的其他易溶盐类，例如NaCl。在这种情况下，与方解石平衡的Ca含量的变化可用下述的溶度积公式来解释：

2+

（2-3-49）

式中，

2+

、

2-

分别表示Ca和CO3的活度系数。NaCl的加入使水溶液的离子强度增大，

2+2-

而Ca和CO3的活度系数则随着溶液离子强度的增大而减小。因此，在Kcal为常数的情况下，NaCl的加入将使与方解石平衡的Ca含量增大，从而使方解石的溶解度升高。一般情况下，矿物在浓溶液要比稀溶液中更易溶，这种效应称为离子强度效应（Ionic Strength Effect）。

2+

2.3.4 方解石和白云石的非全等溶解

当地下水在其流动过程中同时或在不同时间遇到方解石和白云石时，可发生一种矿物溶解、同时又有另一种矿物沉淀的情况，这就是所谓的非全等溶解。表2-3-2给出了不同温度下方解石和白云石的溶度积，据此可以对方解石和白云石的非全等溶解作用进行分析。

表2-3-2 不同温度下方解石和白云石溶解反应的平衡常数

t（℃）

0 5 10 15

pKcal 8.340 8.345 8.355 8.370

pKdol 16.56 16.63 16.71 16.79

pKdol 8.280 8.315 8.355 8.395

1/2

20 25 30

8.385 8.400 8.510

16.89 17.00 17.90

8.445 8.500 8.950

这里首先提出的一个问题是，在10℃的条件下与白云石平衡的地下水接触到了方解石后，将会发生什么样的情况？该问题可通过把白云石的溶度积写成下述的形式来说明：

（2-3-50）

由于通过白云石溶解在水中形成的Ca和Mg的活度相等，故式（2-3-50）中圆括弧内的两项是相等的，因此({Ca}{ CO3})= Kdol。由表2-3-2可见，10℃时Kdol=10

2+

2-1/2

1/2

-8.355

2+2+

，这刚好

等于该温度下方解石的溶度积Kcal，因此这时既不会发生方解石的溶解，也不会发生方解石的沉淀。

当温度大于10℃时（例如30℃），由于({Ca}{ CO3})= Kdol=10

2+

2-2+

2-1/2

-8.95

-8.51

，因此

这时方解石将发生溶解。方解石的溶解把Ca和CO3带入了水中，这便使得水溶液相对白云石处于过饱和状态，从而引起白云石的沉淀。我们把同时发生了方解石溶解和白云石沉淀的作用称为方解石的非全等溶解。

当温度小于10℃时（例如0℃），由于({Ca}{ CO3})= Kdol=10

2+

2+

2-1/2

-8.28

>Kcal=10

-8.34

，因此

这种情况下将发生方解石的沉淀。方解石的沉淀从水中带走了Ca和CO3，其结果必然使水中白云石的活度积小于Kdol。这样，白云石的非全等溶解作用即不可避免地发生了。

2-

白云石的非全等溶解作用也可以由方解石的溶解速度比白云石的溶解速度快而引起。在一个同时含有方解石和白云石的含水层中，如果方解石由于其较快的溶解速度首先达到了平衡状态，而白云石则沿着流动途径继续缓慢地溶解。这种情况下白云石的溶解必然使水溶液相对方解石处于过饱和状态，进而可能使其发生沉淀。这是白云石非全等溶解另一个例子。

参考文献

1. Appelo C A J and Postma D .Geochemistry, Groundwater and pollution. Balkema A

A, Rotterdam, 1993。

2. Bogli Alfred .Karst hydrology and physical speleology. springer-Verlag，1980。

3. Butler , James N .Carbon dioxide equilibrium and their applications.

Addison-wesley publishing company, 1982。

4. Kehew A E .Applied Chemical Hydrogeology. Prentice Hall, 2001。

5. 袁道先，蔡桂鸿 .岩溶环境学. 重庆：重庆出版社，1988。

6. 中国科学院地质研究所岩溶研究组著 .中国岩溶研究. 北京：科学出版社，1979。

7. 钱会 .Bogli混合溶蚀理论及其在实际应用中所存在的问题. 中国岩溶，1996（4）。

本文来源：https://www.bwwdw.com/article/8wmo.html