1. 引言

中国共产党第二十次全国代表大会所作报告中指出：“大自然是人类赖以生存发展的基本条件。尊重自然、顺应自然、保护自然，是全面建设社会主义现代化国家的内在要求。必须牢固树立和践行绿水青山就是金山银山的理念，站在人与自然和谐共生的高度谋划发展。”近年来，随着我国的城市化和工业化水平有了进一步增长，人口密度也不断增加，经济发展的速度不断加快，同时也给水生生态环境造成了较大压力。受化学工业发展的影响，一些污水在排放时会携带不同含量的Fe²⁺、Fe³⁺、Zn²⁺、Cu²⁺、Ca²⁺、Hg²⁺、Ag⁺、Mg²⁺等重金属离子。这些离子掺杂在环境水体中，经由水循环过程后很可能对人体健康带来严重危害。为了进一步加强对环境水质的分析和保护，对水体中的重金属离子进行检测和相应防治之举尤显必要。本文基于此背景下设定一具体且较为复杂的滴定分析体系，旨在探讨如何利用配位滴定原理借助去公式化解析策略，对复杂体系中目标离子含量实现精确分析，并总结这一解析策略在配位滴定体系中的一般性应用方法以及相关延伸问题。

2. 复杂滴定体系的构建与解析

2.1. 复杂定量解析体系的设定构想

基于环境背景，并结合分析化学课程所学内容，在设定一具体的定量解析体系时需要考虑到以下几个层面的问题和要求：

1) 被污染水体中往往存在多种重金属离子，因难以迅速分离去除，这就不可避免地对滴定体系引入了多种副反应，增大了对该体系定量解析的复杂度。考虑到工业生产中这种情况较为常见，故所设定的定量解析体系理应包含多种常见的重金属离子；

2) 为最大程度地实现利用EDTA配位剂准确滴定某一种目标重金属离子的目的，在实际操作中还需排除其他离子的干扰，特别是与目标离子的滴定突跃小区间非常接近的离子。为此，可向体系中引入有效的沉淀剂(或离子)，通过沉淀反应尽力缓解上述干扰的不良影响；

3) 正如2)中所言，实际操作中往往需引入其他试剂来实现准确滴定的目的，但这很可能也同时给滴定体系增添了新的副反应，其原因在于所引入试剂或离子有可能与原体系中存在的金属离子发生配位络合、酸碱交换等反应，因此在设定滴定情景时需充分考虑诸如此类的可能性。

依据上述分析，结合课本附表所列物理化学常数，现对该定量解析体系大致设定如下：

从工业废水中抽样并作相关处理后得到含有Zn²⁺、Ag⁺、Hg²⁺三者的混合溶液，以一定浓度的EDTA (为简化计算，这里不妨认为EDTA和目标离子浓度相同)滴定其中的Zn²⁺ (选择Zn²⁺为目标离子是由于与其他两种离子相比，Zn²⁺的定量解析方法更为有限)，并加入NH₄I作掩蔽剂，在合理误差范围内可认为I⁻与Hg²⁺发生配位络合，与Ag⁺发生沉淀反应，且对目标离子Zn²⁺无直接影响。

2.2. 复杂定量解析体系的分析与突破路径

与大多数分析化学课本所涉及的简单配位滴定体系相比，本文所设定的体系在解析过程中有以下几个难点：

1) 涉及三种金属离子，且从理论角度考虑均可与配位剂EDTA发生络合反应，也即EDTA的副反应较多；

2) 体系中所引入的掩蔽剂NH₄I (主要是I⁻)既可与Hg²⁺发生配位反应，又可与Ag⁺发生沉淀反应，两种干扰金属离子的副反应较多；

3) 值得关注的是，若再考虑Zn²⁺与H₂O的配位反应及 NH 4 + 水解对pH等的影响，将会进一步加大体系的解析难度，同时也会降低本文所探讨原理和方法的实际应用价值，因此在下述解析中忽略Zn²⁺与H₂O的配位反应和 NH 4 + 水解。

由于上述难点及多种干扰因素的局限，在解析该体系时难以使用传统解析策略得到较为精确的结果。为此可采用《分析化学 第二版》(邵利民编著)中所述的去公式化解析策略，具体求解过程如下所述。

2.3. 滴定体系的构建与解析

基于2.1中设定的框架，查阅相关物理化学常数后对该问题涉及到的相关数据作以细化如下：

现从工业废水中抽样并作相关处理后得到Zn²⁺、Ag⁺、Hg²⁺三者的混合溶液，三者的分析浓度均为0.020 mol∙L⁻¹。以同浓度的EDTA溶液滴定其中的Zn²⁺，控制pH = 5.50，并加入NH₄I使得[I⁻] = 0.10 mol∙L⁻¹。相关参数如下：

K Z n Y = 3.16 × 10 16 ； K A g Y = 2.09 × 10 7 ； K H g Y = 5.01 × 10 21 ； α Y ( H ) | p H = 5.50 = 3.24 × 10 5 ；

Hg²⁺-I⁻配离子： β 1 = 10 12.87 ； β 2 = 10 23.82 ； β 3 = 10 27.60 ； β 4 = 10 29.83 ； K s p , A g I = 9.3 × 10 − 7 。

依据上述设定，对该滴定体系的解析较为复杂，主要涉及到三种金属离子和EDTA各自的一个或多个反应。为得到精确解析的结果，最关键的部分则是滴定过程等量关系式的建立，据此可绘制出滴定曲线并计算滴定突跃。下基于去公式化解析策略对此复杂滴定体系作简明解析： [1]

分别用V和V_sp表示某一滴定时刻和化学计量点时加入EDTA溶液的体积，并令： R = V e p V s p

根据反应物浓度和滴定反应的化学计量关系可知：被测物Zn²⁺溶液的体积为V_sp；

当加入V mL滴定剂时，溶液的总体积为(V + V_sp)，各物质的分析浓度为：

c Z n = c A g = c H g = 0.020 V s p V + V s p = 0.020 R + 1

c E D T A = 0.020 V V + V s p = 0.020 R R + 1

基于上述条件，可列出Zn、Ag、Hg以及EDTA的MBE式：

[ Z n Y ] + [ Z n ] = c Z n = 0.020 R + 1 (1)

[ A g Y ] + [ A g ] + [ A g I ( s ) ] = c A g = 0.020 R + 1 (2)

[ H g Y ] + [ H g ] + [ H g I ] + [ H g I 2 ] + [ H g I 3 ] + [ H g I 4 ] = c H g = 0.020 R + 1 (3)

[ Z n Y ] + [ A g Y ] + [ H g Y ] + [ Y ″ ] = c E D T A = 0.020 R R + 1 (4)

其中， [ A g I ( s ) ] 表示沉淀AgI的假想浓度，仅用于建立等量关系；以 [ Y ″ ] 代替 [ Y ] 是为了减小数值运算误差，考虑到：

[ Y ″ ] = [ Y ] ⋅ α Y ( H ) (5)

可知 [ Y ″ ] (加双撇号是为了区别于真实表观浓度)数值远大于 [ Y ] ，故可以实现减小误差的预设目的。

为消去(4)式中金属配离子项，现用分布分数来替换[ZnY]和[HgY]，应有：

[ Z n Y ] = K ″ Z n Y [ Y ″ ] K ″ Z n Y [ Y ″ ] + 1 ⋅ 0.020 R + 1 (6)

[ H g Y ] = K ″ H g Y [ Y ″ ] K ″ H g Y [ Y ″ ] + α H g ⋅ 0.020 R + 1 (7)

上式中加双撇号是为了区别于真实条件稳定常数，且有：

K ″ Z n Y = K Z n Y α Y ( H ) (8)

K ″ H g Y = K H g Y α Y ( H ) (9)

α H g = 1 + β 1 [ I − ] + β 2 [ I − ] 2 + β 3 [ I − ] 3 + β 4 [ I − ] 4 (10)

再考虑通过建立平衡关系替换(4)式中[AgY]。首先利用 K s p , A g I 建立等量关系，即有：

K s p , A g I = [ A g ] [ I − ] (11)

又有：

K ″ A g Y = [ A g Y ] [ A g ] [ Y ″ ] (12)

联立(11) (12)可表示出[AgY]，得到：

[ A g Y ] = K ″ A g Y [ Y ″ ] ⋅ K s p , A g I [ I − ] (13)

类似(8) (9)式，有：

K ″ A g Y = K A g Y α Y ( H ) (14)

现已将(4)式中各金属配离子项替换，将(6) (7) (13)式代入(4)式中得到：

K ″ Z n Y [ Y ″ ] K ″ Z n Y [ Y ″ ] + 1 ⋅ 0.020 R + 1 + K ″ A g Y [ Y ″ ] ⋅ K s p , A g I [ I − ] + K ″ H g Y [ Y ″ ] K ″ H g Y [ Y ″ ] + α H g ⋅ 0.020 R + 1 + [ Y ″ ] = c E D T A = 0.020 R R + 1 (15)

考虑R的定义式，对上式作相应替换，有：

K ″ Z n Y [ Y ″ ] K ″ Z n Y [ Y ″ ] + 1 ⋅ 0.020 R + 1 + K ″ A g Y [ Y ″ ] ⋅ K s p , A g I [ I − ] + K ″ H g Y [ Y ″ ] K ″ H g Y [ Y ″ ] + α H g ⋅ 0.020 R + 1 + [ Y ″ ] = c E D T A = 0.020 V V + V s p (16)

2.4. 滴定曲线的绘制与准确滴定判别

由上述分析可知，V (某一滴定时刻加入EDTA溶液的体积)与 [ Y ″ ] 之间的函数关系已经建立，再考虑建立 [ Y ″ ] 与[Zn]之间的等量关系。为此，可借助Zn的分布分数：

[ Z n ] = 1 K ″ Z n Y [ Y ″ ] + 1 ⋅ 0.020 V s p V + V s p (17)

至此，绘制滴定曲线所需的等量关系式已全部得到。为同时得到可视化兼具高精度的结果，可借助MATLAB程序实现滴定曲线的绘制(取V_sp = 20 mL)，只需将(17)式代入(16)式中以消去中间变量 [ Y ″ ] ，得到包含[Zn]与V的等式即可。相关代码参见附录。

得到的滴定曲线如图1所示。

为对该滴定体系进行准确滴定判别，需获取R与[Zn]的等量关系式，可间接得到此关系。由上述分析可知，R与 [ Y ″ ] 之间的函数关系已经建立，再考虑建立 [ Y ″ ] 与[Zn]之间的等量关系。为此，可借助Zn的分布分数： [1]

[ Z n ] = 1 K ″ Z n Y [ Y ″ ] + 1 ⋅ 0.020 R + 1 (18)

至此，计算滴定突跃所需的等量关系式已经全部得到。仍借助MATLAB程序辅助分析，只需将(18)式代入(15)式中以消去中间变量 [ Y ″ ] ，得到包含[Zn]与R的等式，然后求解 R = 0.999 和 R = 1.001 对应的[Zn]，从而获得滴定突跃。值得一提的是，使用MATLAB程序求解时，不必写出这个方程，以 [ Y ″ ] 为中间变量，即可快速完成输入。所用MATLAB代码参见附录。

代入相关数据可求得：

p Z n R = 0.999 = 4.99

p Z n R = 1.001 = 6.81

然后计算出滴定突跃：

| p Z n R = 0.999 − p Z n R = 1.001 | = 1.82

滴定突跃大于0.4，因此可以实现准确滴定。

2.5. 传统解析方法的求解思路

本问题也可使用传统分析化学方法中的林邦公式近似求解，但是需作一定的近似处理。具体求解思路如下： [1]

林邦公式的原始形式为：

E t = 10 Δ p M ′ − 10 − Δ p M ′ c M , s p K ′ M Y × 100% (19)

c Z n , s p ⋅ K ′ Z n Y ≥ ( 10 0.2 − 10 − 0.2 0.001 ) 2 (20)

c Z n , s p = 0 .010 mol / L (21)

因此问题的关键转化为求解上述体系下的 K ′ Z n Y 。由定义：

K ′ Z n Y = K Z n Y α Z n α Y = K Z n Y 1 ⋅ ( α Y ( H ) + α Y ( A g ) + α Y ( H g ) − 2 ) (22)

上式中各系数可由以下各式求解：

α Y ( A g ) = 1 + K A g Y [ A g + ] = 1 + K A g Y ⋅ K s p , A g I [ I − ] (23)

α Y ( H g ) = 1 + K H g Y [ H g 2 + ] (24)

[ H g Y ] + [ H g ] + [ H g I ] + [ H g I 2 ] + [ H g I 3 ] + [ H g I 4 ] = c H g = 0.020 R + 1 (3)

c H g ≈ c H g , s p = 0.010 mol / L (25)

据此，联立(3) (21) (22) (23) (24) (25)并代入(20)式中验证发现：该不等式在题述体系下成立。

也即：由传统求解模式同样可以得到能够准确滴定Zn²⁺的结论。

3. 结果与讨论

3.1. 实例分析结论与结果

3.1.1. 分析结论

上文选取一含有三种金属离子、涉及到配位平衡和沉淀溶解平衡等诸多反应的复杂体系为例进行解析，以Zn²⁺作为目标金属离子，对滴定分析其含量和能否准确滴定判别分别进行计算讨论，所得结论如下：

1) 在此体系下，可以通过配位滴定原理利用EDTA对目标离子(Zn²⁺)的含量进行分析测算；

2) 因对Zn²⁺的配位滴定曲线中滴定突跃大于0.4，故在本框架条件下可以实现对Zn²⁺的准确滴定。

3.1.2. 成果与不足

本文以废水中重金属离子含量的分析为例，基于传统分析化学原理，结合《分析化学 第二版》(邵利民编著)中所述的去公式化解析策略，对上文所述一较为复杂的配位滴定定量体系作了相对精准的解析，主要取得了以下成果：

1) 摆脱了传统近似计算策略的桎梏，从原理和思路上对该定量滴定体系作以精准的解析，更符合分析化学以“分析”为主的核心和本质；

2) 借助MATLAB程序较为精确地绘制了配位滴定曲线，由曲线可以直观地判断能否准确滴定目标离子，并由去公式化计算策略从理论上证明之；

3) 尽可能减小了解析过程中由于近似处理而产生的误差，较复杂的计算均借助MATLAB代码实现；

在应用层面更为广阔、应用方式更为多样的具体生产生活中，本文所述的定量解析方法具有更强的应用价值；在响应我国加快生态文明建设进程的号召之下，上述环境分析原理兼具实用性与现实意义。然而，与许多理论分析方法一样，文中所述方法原理在实际生产应用中存在着诸多的局限性，具体而言，可分为以下几个层面：

1) 成本问题：经济方面的局限性在众多实际应用体系下成本问题都是一个不可忽略且难以规避的层面。具体而言，实际工业或农业等含重金属离子的废水量往往是一个庞大的数字，相应所需要的试剂用量也远高于实验室分析用量。以本文所述体系为例，配位滴定体系中除过最常用的EDTA用量较大以外，还可能需要多种掩蔽剂、沉淀剂等预处理试剂，这就可能带来了数额不菲的成本预算，为体系的解析造成了根本性的阻碍。

2) 复杂度问题：本文所述体系相比各高等院校所用分析化学课本提及的配位滴定体系，具有相当高的解析难度，然而与工农业生产中的实际废水样相比，仍已做了许多简化处理。具体而言，上文体系中仅涉及目标金属离子与干扰金属离子和EDTA的配位平衡、金属离子和干扰配体的络合平衡、沉淀的溶解平衡等，实际体系中可能还存在着不溶性物质、特殊金属离子、酸碱度等多种因素的干扰，故从此角度上考虑，上述解析策略的应用性被大大削弱；

3) 预处理及后处置等问题：正如2)中所言，欲将去公式化方法策略应用于实际体系中，需考虑可能存在的各种干扰因素，并对其进行预处理以保证配位滴定顺利进行，这就对仪器、操作、成本等多方面施以更高的要求。除此之外，在某些特殊体系中，出于节约资源、保护环境等目的，在完成主要滴定部分后还需对母液等进行相应的后处置，这在无形中增加了实际应用的难度与挑战性。

3.2. 去公式化解析策略在配位滴定体系中的理论计算方法总结

3.2.1. 相关名词定义

分析化学对配位滴定体系的传统解析中需引入表观浓度、条件稳定常数等定义，若利用去公式化策略进行解析，为叙述方便，需在此基础上定义几个新的概念：

1) 精确条件稳定常数：作如下定义

K ″ M Y = [ M Y ] [ M ] [ Y ″ ] = K M Y α Y ( H )

以代替条件稳定常数，加双撇号是为了区别于一般的条件稳定常数；

2) 精确表观浓度：作如下定义：

[ Y ″ ] = [ Y ] ⋅ α Y ( H )

以代替表观浓度[Y’]，减小数值运算误差；加双撇号是为了区别于一般的表观浓度。

由此，在解析过程中便不再出现EDTA的一般表观浓度和配合物的条件稳定常数。但对于逐级稳定常数β_n、金属离子和EDTA相应副反应系数α等广泛应用的物理量，在去公式化解析策略中仍发挥着重要的建立等量关系作代换的作用。

3.2.2. 涉及多种金属离子配位滴定体系的去公式化解析策略

首先需根据已有信息确定出相关物质或离子的分析浓度和平衡浓度，然后列出体系中各物质和离子(包括目标金属离子、干扰金属离子、EDTA、干扰配体等)的MBE式，接着对式中量通过相应代换建立新的更为简明的等量关系式，以减少未知量的个数，最后分析求解即可。在上述分析步骤中，作代换是最为核心的一步，下列举出部分作代换的思考角度：

1) 对金属与EDTA形成配合物的平衡浓度[MY]用经由精确条件稳定常数 K ″ M Y 确定的分布分数作代换，如：

[ M Y ] = K ″ M Y [ Y ″ ] K ″ M Y [ Y ″ ] + α M ⋅ c M

当金属离子在体系中不存在副反应时，上述代换式中 α M = 1 ；

2) 对EDTA的平衡浓度[Y]利用精确表观浓度 [ Y ″ ] 的定义式作代换，即：

[ Y ″ ] = [ Y ] ⋅ α Y ( H )

上式中，EDTA的酸效应系数 α Y ( H ) 在一定pH条件下为一已知量；

3) 对金属的平衡浓度[M]用经由精确条件稳定常数 K ″ M Y 确定的分布分数作代换，如：

[ M ] = 1 K ″ M Y [ Y ″ ] + α M ⋅ c M

当金属离子在体系中不存在副反应时，上述代换式中 α M = 1 ；

4) 存在干扰配体时，对金属离子与干扰配体L形成配合物的平衡浓度[ML_n]利用逐级稳定常数β_n作代换，也即：

[ M L n ] = β n [ M ] [ L ] n

上式中，逐级稳定常数β_n在一定配体条件下为一已知量；

5) 存在干扰配体时，对金属离子的副反应系数 α M ( L ) 利用逐级稳定常数β作代换，也即：

α M ( L ) = 1 + β 1 [ L ] + β 2 [ L ] 2 + ⋅ ⋅ ⋅ + β n [ L ] n

上式中，逐级稳定常数β_n在一定配体条件下为一已知量；

6) 存在干扰金属离子N时，对EDTA的副反应系数 α Y ( N ) 利用其稳定常数K_NY作代换，也即：

α Y ( N ) = 1 + K N 1 Y [ N 1 ] + K N 2 Y [ N 2 ] + ⋅ ⋅ ⋅ + K N n Y [ N n ]

上式中，稳定常数K_NY在一定条件下为一已知量；

7) 涉及酸碱平衡时，可对对弱酸及所有酸根离子的加和式用代求酸根离子和其分布系数 δ 的关系式作代换，例如：

[ H 3 A ] + [ H 2 A − ] + [ H A 2 − ] + [ A 3 − ] = [ H 3 A ] δ H 3 A = [ H 2 A − ] δ H 2 A − = [ H A 2 − ] δ H A 2 − = [ A 3 − ] δ A 3 −

8) 涉及沉淀溶解平衡时，常利用沉淀的溶度积常数对其中的离子作代换，这里不作赘述。

在实际情景下，根据已知条件和信息可灵活处理，这里仅举几例说明去公式化解析策略的核心步骤。

3.2.3. 利用去公式化解析策略进行配位平衡体系准确滴定判别

1) 先做预处理：分别用V和V_sp表示滴定终点和化学计量点时加入EDTA的总体积，并令 R = V V s p ；根据反应物浓度和滴定反应的计量关系，容易知道被测物溶液的体积为V_sp；终点时溶液的总体积等于(V + V_sp)；

2) 类似3.2.3方式作代换，消去其中的[Y’]_ep或[Y”]_ep，得到包含[M’]_ep或[M]_ep和R的方程；

3) 分别代入R₁ = 0.999和R₂ = 1.001，得到相应[M’]₁和[M’]₂，进而计算 Δ p M ′ (若大于0.4则可以准确滴定) [1] 。

3.3. 去公式化解析策略在生产生活中应用举例

结合本文所述原理并查阅相关资料可知，该方法策略可用于众多体系中并取得较精准的结果，具体有如下几类：

1) 配位滴定原理应用于多体系下金 [ Y ″ ] 属离子含量检测：配位滴定法作为传统化学分析中四大滴定最重要的一类之一，主要用于样品溶液中金属离子含量的测定。由于EDTA的结构优势(EDTA可与几乎所有金属离子发生配位反应生成非常稳定的螯合物结构)，大部分的配位滴定体系均选用EDTA作为螯合剂或滴定剂。基于此，可将配位滴定原理应用于生产生活中的多个领域，例如：自来水硬度的测定、土壤样品中重金属离子含量的测定、化学需氧量的分析等；

2) 去公式化解析策略应用于多环境下化学分析与测量：从上文所述分析思路不难看出去公式化解析策略的优势与便捷。因此，只需把握住各个化学环境中的守恒思想，便能将该解析策略推广延伸，例如：酸碱平衡体系下利用电荷守恒(CBE)式结合分布分数、氧化还原平衡体系下结合能斯特方程、沉淀平衡体系下结合溶度积常数及物料守恒(MBE)式等即可实现去公式化解析的目的；

3) MATLAB等数据分析软件应用于多领域计算化学：除过去公式化解析策略，本文所述分析方法中另一大亮点是MATLAB程序的有效应用。在包含化学与材料科学的众多研究领域中，MATLAB等数据分析软件可高效应用于数据处理、数值计算和模拟实验等方向，可提供丰富的工具和相关函数以及各种数值计算方法和算法，非常便捷地实现对科学计算中复杂问题的求解，可作为计算化学中的一种方便高效的辅助工具。

3.4. 废水样处理方法的拓展延伸

3.4.1. 分光光度法测废水样中全铁含量

1) 方法原理：废水样中铁含量超标将会严重影响水的色、嗅、味等感官性状，还会对废水生化处理构筑物中的微生物有较强的致死作用，对人体的健康造成直接或间接的危害。分光光度法是测算废水样中全铁含量的一种精确度较高的仪器分析方法，常选用邻二氮菲作为测定过程中铁的显色剂，它是测定Fe(II)的高灵敏度高选择性试剂。具体而言，在pH约2~9的范围内，溶液中的Fe²⁺可与邻二氮菲反应生成一种极稳定的橘红色配合物，如图2所示：

其 lg K 稳 = 21.3 ( 20 ˚ C ) ，在510 nm波长下有最大吸收峰，相应摩尔吸光系数为 ε = 1.1 × 10 4 。

2) 简要步骤与测算方法：为得到更为精确的测量结果，实际操作中往往采用条件实验并控制单因素变量。为此可先用盐酸羟胺将试样溶液中的Fe³⁺全部还原为Fe²⁺ (前者也可与邻二氮菲发生配位络合反应，但生成的配合物稳定性及相应吸光度曲线精度、测算灵敏度等均不如后者；除此之外，由于Fe³⁺完全沉淀的pH较低，不利于对溶液pH的控制)，后待二者充分反应后通过分光光度计测算溶液吸光度A，并与事先测绘好的吸光度–铁含量标准曲线相比较从而确定试样中全铁含量 [2] 。

3.4.2. 重铬酸钾法测废水样中化学需氧量COD

1) 方法原理：化学需氧量(Chemical Oxygen Demand，简称COD)，是指水体中易被强氧化剂氧化的还原性物质所清耗的氧化剂的量，结果折算成氧的量(以mg/L计)，其反映了水中受还原性物质污染的程度。高锰酸钾、重铬酸钾等强氧化剂常用于测定废水样的COD，由于这类水中含有许多结构复杂的有机物，用高锰酸钾难以充分氧化且不易控制反应条件，因此往往选用氧化效果更好的重铬酸钾作氧化剂，并以硫酸亚铁铵作标准溶液通过滴定的分析方法测定COD。

2) 简要步骤与测算方法：对废水样作预处理(主要是去除不溶性杂质、酸化、排除可能具有氯化物的干扰等)后加入一定量重铬酸钾溶液，在Ag₂SO₄催化下加热回流2 h，充分反应并冷却后加入2~3滴亚铁灵指示剂，用硫酸亚铁铵标准溶液滴定至溶液由黄色到绿蓝色最终变为红蓝色为终点，记录消耗标准液的用量为V₁；做空白实验(即以等量蒸馏水代替废水样重复上述步骤)，记录消耗标准液的用量为V₀。根据化学反应计量关系，可根据下式计算COD：

C O D | C r = ( V 0 − V 1 ) ⋅ c M O 2 ⋅ 100 V 2 ( mg / L )

上式中，c为硫酸亚铁铵标准溶液的浓度(mol/L)；V₀为空白消耗标准液的体积(mL)；V₁为废水样消耗标准液的体积(mL)；V₂为水样体积 [3] [4] 。

附录

1. 滴定曲线绘制MATLAB原代码

K1=3.16e11/3.24;

K2=2.09e2/3.24;

K3=5.01e16/3.24;

Ksp=9.3e-17;

a=7.16e25;

%各金属离子在pH=5.5的条件稳定常数

VLowerlimit=0;

VUpperlimit=40;

Y=[linspace(0,1e-10,10000) linspace(1e-10,1e-8,10000) linspace(1e-8,1e-6,10000) l inspace(1e-6,1e-3,10000) linspace(1e-3,1e-2,10000)];

aux1=K1*Y./(K1*Y+1);

aux2=K2*Y./(K2*Y+a);

aux3=K3*Y*Ksp./0.1;

V=(0.4*aux1+0.4*aux2+20*aux3+20*Y)./(0.02-aux3-Y);

pZn=-log10(0.4./(V+20)./(K1*Y+1));

Filter=find((V>VLowerlimit)&(V

figure;plot(V(Filter),pZn(Filter));

2. 滴定突跃计算MATLAB原代码

K1=3.16e11/3.24;

K2=2.09e2/3.24;

K3=5.01e16/3.24;

Ksp=9.3e-17;

a=7.16e25;

%各金属离子在pH=5.5的条件稳定常数

VLowerlimit=0;

VUpperlimit=40;

Y=[linspace(0,1e-10,10000) linspace(1e-10,1e-8,10000) linspace(1e-8,1e-6,10000) linspace(1e-6,1e-3,10000) linspace(1e-3,1e-2,10000)];

aux1=K1*Y./(K1*Y+1);

aux2=K2*Y./(K2*Y+a);

aux3=K3*Y*Ksp./0.1;

V=(0.4*aux1+0.4*aux2+20*aux3+20*Y)./(0.02-aux3-Y);

pZn=-log10(0.4./(V+20)./(K1*Y+1));

Filter=find((V>VLowerlimit)&(V

figure;plot(V(Filter),pZn(Filter));

NOTES

^*通讯作者。

参考文献