的实验涉及测定氯残留在不同时间和不同初始氯剂量下使用,水样品进行了常规治疗包括快速砂过滤和消毒之前收集形成Matsapha水处理厂,随后被用于实验的试验。收集的样本被带到Eswatini大学的化学实验室。不同剂量的氯被添加到多个样品,和氯残留测定在预先确定的时间间隔的时间。方法用于测定余氯碘量滴定的滴定在考试标准方法如上所述的水和废水 26]。方法是基于当碘化钾的原则是添加到样本含有余氯的水的pH值8或更少,余氯释放碘碘化钾和解放碘与硫代硫酸钠滴定(Na₂年代₂O₃滴定剂)。

为了使滴定化学计量,样本的pH值降低,通过添加醋酸在3和4之间。预计方法检出限为0.04 mg / L如果0.001 N硫代硫酸钠滴定液使用。已经使用的所有化学试剂级质量和新鲜的解决方案准备每次实验重复了在不同时期。

根据预期的浓度,合适的样本体积被滴定剂,硫代硫酸钠(Na₂年代₂O₃)体积不超过20毫升和淀粉体积指标也高于0.2毫升。样品体积,5毫升醋酸添加之后,1克碘化钾粉衡量重量的平衡。在这之后,样本与硫代硫酸钠滴定(Na₂年代₂O₃)适当的常态。滴定剂对于预期的氯浓度较低的样品,准备是正常的低。滴定是直到黄色几乎消失。然后,1毫升的淀粉溶液添加到蓝色消失了。为了弥补方法错误,空白滴定也使用distilled-deionised水,通过执行相同的程序用于实际样品。至少三个重复执行每个决心。滴定后,mg / L氯残留在示例使用以下公式确定: (1) 毫克 l 氯 = 一个 ± B ∗ N ∗ 35450年 毫升的样品 , 在哪里 一个滴定剂的体积用于样本, B是用于空白滴定剂的体积,然后呢 NNa的正常吗₂年代₂O₃滴定液。

随时间变化的反应速率常数方程所建议的菲利普et al。 5)用于数学公式。根据这个方程,concentration-averaged反应速率随时间变化是由以下表达式: (2) d k t d t = C t k t 2 − ∑ 我 = 1 N k 我 2 X 我 t X T 。

最初的二阶变异氯衰减率与摩尔浓度的氯和反应物和方程( 2)是由派生而来 (3) d C t d t = − k t C t X t = d X T d t 。 在方程( 2)和( 3), k _t摩尔重反应速率, C _t是氯残留测量时间吗 t,X _它反应物的摩尔浓度吗 我测量时间 t与氯反应, X _T是反应物施加的总物质的量浓度氯需求,然后呢 N反应物氯消耗的总数目前在水里。

与此同时,molar-averaged反应速率 k _t被定义为 (4) k t = ∑ 我 = 1 N k 我 X 我 t ∑ 我 = 1 N X 我 t = ∑ 我 = 1 N k 我 X 我 t X T , 在哪里 k _我反应物的反应速率 我在摩尔浓度 x _我并施加氯需求。其他符号定义方程( 2)。方程( 2)表示的二阶特性的变化molar-averaged反应速率常数 k _t随着时间的推移。很明显从方程( 2),其他因素保持不变,反应速率常数 k _t变化是一个一阶反应对氯残留。的聚合效应的存在氯残留这二阶变化的速率常数,因为随着氯浓度的增加速率常数趋于减少。类似的观察也指出一些研究人员通过实证实验测定反应速率常数之间的关系和最初的实验中使用氯 13, 18, 22, 23]。增加的速率常数似乎减少反比氯剂量用于实验。菲利普et al。 5)还认为,速率常数随时间变化总是负在方程( 2这个方程),括号中的表达式可以证明只有负面或零。换句话说,它不能采取积极的价值观。

方程( 2),由菲利普et al。 5摘要]扩展提供一个公式,与时间变化的反应速率 k _t用氯残留现在和反应物的摩尔浓度。处理不确定的词的难度在支架的右边方程( 2反应物)和涉及个人可以通过数学来解决产品的平均个人的反应率和反应物的摩尔浓度。这个平均后验证了理论公式的误差项和蒙特卡罗模拟的平均误差在不同价值观这些产品条款,即。反应物的摩尔浓度和个人的反应速率。

从括号中的表达式方程( 2),这样调整 (5) ∑ 我 = 1 N k 我 2 X 我 t = ∑ 我 = 1 N k 我 2 X 我 t 2 X 我 t = ∑ 我 = 1 N k 我 X 我 t 2 X 我 t 。

重新排列的右边上面的表达式,我们得到的 (6) ∑ 我 = 1 N k 我 X 我 t 2 X 我 t = ∑ k 我 X 我 t 。 k 我 X 我 t X 我 t 。

现在定义的平均值 k _我 X _它这样 (7) k 我 X 我 t av = ∑ k 我 X 我 t N 。

因为它已经定义在方程( 4), (8) k t = ∑ k 我 X 我 t X T 。

方程( 4)和( 7)相结合 (9) k 我 X 我 t av = ∑ k 我 X 我 t N = k t X T N 。

现在回归方程( 5再一次,我们得到的 (10) ∑ 我 = 1 N k 我 X 我 t 2 X 我 t = ∑ k 我 X 我 t · k 我 X 我 t X 我 t = k 我 X 我 t av · ∑ k 我 X 我 t X 我 t 。 在上面的表达式中,( k _我 X _它/ X _它)作为考虑因素。这也类似于一个连续函数的重,也就是说, (11) ∫ x 。 f x d x = x ¯ ∫ f x d x 。

此外,通过使用类似的平均方法,我们得到的 (12) ∑ k 我 X 我 t X 我 t = k 我 X 我 t av · ∑ 1 X 我 t 。

在这种情况下,1 / X _它作为考虑因素。

总的来说, (13) ∑ 我 = 1 N k 我 X 我 t 2 X 我 t = k 我 X 我 t av 2 。 ∑ 1 X 我 t = k t 2 X T 2 N 2 。 ∑ 1 X 我 t 。

使用加权平均方程( 13离散值)可以验证理论以及实践使用蒙特卡罗模拟使用随机变化的个人反应物的速率常数和反应物的摩尔浓度如下所述。首先,理论基础是解释如下: (14) ∑ 我 = 1 N k x 2 ≈ k x av · ∑ 我 = 1 N X = 1 N ∑ 我 = 1 N k x ∑ 我 = 1 N X 。

工作有点右边表达式,得到 (15) 1 N ∑ 我 = 1 N k x ∑ 我 = 1 N X = 1 N ∑ 我 = 1 N k x 2 + ∑ 我 = 1 N ∑ j = 1 , j ≠ 我 N k 我 x 我 x j 。

现在定义术语的区别在摩尔反应物Δ的表情 x _ij,我们得到 (16) Δ x 我 j = x 我 − x j , 这 x j = x 我 − Δ x 我 j 。

用上述的表达式 x _j项,我们得到 (17) 1 N ∑ 我 = 1 N k x ∑ 我 = 1 N X = 1 N ∑ 我 = 1 N k x 2 + ∑ 我 = 1 N ∑ j = 1 , j ≠ 我 N k 我 x 我 x 我 − Δ x 我 j 。

右边上面的表达变化 (18) 1 N ∑ 我 = 1 N k x ∑ 我 = 1 N X = 1 N ∑ 我 = 1 N k x 2 + N − 1 ∑ 我 = 1 N k x 2 − ∑ 我 = 1 N ∑ j = 1 , j ≠ 我 N k 我 x 我 Δ x 我 j , 1 N ∑ 我 = 1 N k x ∑ 我 = 1 N X = ∑ 我 = 1 N k x 2 − 1 N ∑ 我 = 1 N ∑ j = 1 , j ≠ 我 N k 我 x 我 Δ x 我 j 。

因为对称,很容易显示 (19) Δ x 我 j = − Δ x j 我 。

使用这个对称,表达式是进一步减少 (20) 1 N ∑ 我 = 1 N k x ∑ 我 = 1 N X = ∑ 我 = 1 N k x 2 − 1 N ∑ 我 = 1 N ∑ j = 我 + 1 N k 我 x 我 − k j x j Δ x 我 j 。

定义术语Δ的区别 kx这样 (21) Δ k x 我 j = k 我 x 我 − k j x j , 和使用这个术语在上面的方程,我们得到 (22) 1 N ∑ 我 = 1 N k x ∑ 我 = 1 N X = ∑ 我 = 1 N k x 2 − 1 N ∑ 我 = 1 N ∑ j = 我 + 1 N Δ k x 我 j Δ x 我 j 。

因此, (23) ∑ 我 = 1 N k x 2 = 1 N ∑ 我 = 1 N k x ∑ 我 = 1 N X + E , 在误差项 E被定义为 (24) E = 1 N ∑ 我 = 1 N ∑ j = 我 + 1 N Δ k x 我 j Δ x 我 j 。

上述误差项总和(即二阶差异。,Δ kx次Δ x与其他相比),往往是小条款包含的方程 E:

第二个实践验证是使用蒙特卡罗模拟的误差项,发现平均以及托架的变异误差项的范围 E。为此,假设有五氯要求与个人反应物反应率 k _我不同的最低10倍,最大1000倍和反应物的摩尔浓度 x _我也不同的最大最小10倍和1000倍。蒙特卡罗模拟是为每个范围的变化做了10000次,用下面的结果见表 1。

如表所示 1,平均误差 E约12%与95%信心躺在12.11和12.88之间范围非常狭窄的范围。因此误差项相对较小,使用速率方程的推导过程的考虑因素是有道理的。也是如此的几何和谐波误差项 E在8%和3%的顺序,分别如表所示 1。这些数据支持的启发式使用平均用于推导方程( 13)的范围动态变化的速率和反应物摩尔浓度。

现在回归方程( 13)和使用算术平均( X _一个)和调和平均数( X _H反应物的摩尔浓度的施加氯需求, (25) X ¯ 一个 = X T N , X ¯ H = N ∑ 我 = 1 N 1 / X 我 t 。

方程的表达式( 13)现在可以写成 (26) ∑ 我 = 1 N k 我 X 我 t 2 X 我 t = k t 2 X T 2 N 2 · ∑ 1 X 我 t = k t 2 X T X T N · 1 N · ∑ 1 X 我 t = k t 2 X T X ¯ 一个 N N X ¯ H 。

进一步简化,我们得到 (27) ∑ 我 = 1 N k 我 2 X 我 t = ∑ 我 = 1 N k 我 X 我 t 2 X 我 t = k t 2 X T X ¯ 一个 X ¯ H 。

现在,原来的表达式给出了方程( 2)是 (28) d k t d t = C t k t 2 − ∑ 我 = 1 N k 我 2 X 我 t X T , 可以写成 (29) d k t d t = C t k t 2 − ∑ 我 = 1 N k 我 2 X 我 t X T = C t k t 2 − k t 2 X T X ¯ 一个 / X ¯ H X T ,

这是一样的吗 (30) d k t d t = = C t k t 2 1 − X ¯ 一个 X ¯ H 。

在括号改变表达的顺序,我们获得 (31) d k t d t = − C t k t 2 X ¯ 一个 X ¯ H − 1 。

上面的表达式在方程( 31日)表明,反应速率与对二阶 k _t和一阶 C _t。反应速率随时间总是消极自括号表达式总是积极(算术平均值总是大于调和平均数)。

方程( 31日)因此取代了个人反应速率由concentration-weighed聚合反应速率 k _t因为它是很容易测量。concentration-weighed聚合反应速率 k _t影响算法的比例调和平均数的反应物。如果反应物的摩尔浓度接近级, (32) X ¯ 一个 X ¯ H ≈ 1 , X ¯ 一个 X ¯ H − 1 ⟶ 0。

率 k _t成为常数( dk _t/ dt在方程(= 0 31日))和氯衰减模型方法一阶速度反应。相比之下,如果有重要反应物的摩尔浓度的变化,比例 (33) X ¯ 一个 X ¯ H − 1 ≫ 0 , 适用。氯衰变遵循二级反应,深受反应物摩尔浓度。因为表达的有效性 (34) X ¯ 一个 X ¯ H − 1 ≥ 0 , 的算术平均值总是大于谐波的意思是,时间的衰减速度concentration-weighed反应速率总是负的,随时间减少。这也是通过经验观察到的文献中逆关系氯衰减率和初始氯试验得到的剂量。很明显从方程( 31日), k _t减少速度与初始氯剂量的增加。言外之意是,氯衰减曲线变得平坦(二阶)更高的初始氯剂量。换句话说,在配电系统中,需要长时间的氯衰减比一阶速率衰减和残差往往更大。

为了确定氯衰减率,有必要解决涉及方程(微分方程 31日)和方程( 36),即 (35) d k t d t = − C t k t 2 X ¯ 一个 X ¯ H − 1 , (36) d C t d t = − k t C t X t = d X T d t 。

确定衰变率从以上两个方程,它是必要的,以确定 X _t和 k _t值,这两个方程之间的交替。换句话说, X _t和 C _t价值观是决定从方程( 36)和 k _t从方程( 31日)。然而,首先,初始利率 k₀和 X₀以及反应物的比例调和平均数的计算需要确定。这就是最初的提议 k₀在方程(值的表达式 31日),也就是说, (37) X ¯ 一个 X ¯ H − 1 , 可以由校准。这意味着假设不同的值 k₀和 X ¯ 一个 / X ¯ H 和选择的值之间的误差平方和最小观察和计算氯残留浓度在不同的时间测量。可能的组合 k₀和 X ¯ 一个 / X ¯ H 可以确定使用进化算法或蒙特卡罗方法。

也可以联系 X ¯ 一个 / X ¯ H 涉及价值反应物化学物种参数,例如,TOC和紫外吸收值。一种可能性是使用TOC /紫外线比有关 X ¯ 一个 / X ¯ H 并利用这种关系模型变量衰变率 k _t由方程( 31日)。

也提出了初始反应物的浓度 X₀的摩尔当量总有机碳(TOC)的浓度。的实际价值 X₀也可以通过稀释实验通过运行一个氯衰减实验确定连续稀释和推断相对应的一级反应速率常数未稀释的样本。的 X₀摩尔单位价值大于或等于最大氯需求产生的样本。

程序求解方程( 31日)和( 36)同时遵循使用向前欧拉方法。前两个步骤的过程说明如下。让最初的假设值 k在 t= 0 X ¯ 一个 / X ¯ H 是 k₀和 X ¯ 一个 / X ¯ H = 一个 ,分别。此外,让初始反应物的浓度 X₀从样品的TOC值估计,实际的摩尔的价值 X₀通过氯衰减实验进行一系列的稀释的原始样本,和 X₀值外推到未稀释的样本。使用一个时间步Δ t和向前欧拉方法求解一阶微分方程,在下一个时间步的值 t₁= t₀+Δ t确定如下。

K₁决定使用方程( 31日),也就是说, (38) k 1 − k 0 Δ t = − C 0 k 0 2 一个 − 1 。

下一个 C₁和 X₁决心从方程( 36),也就是说, (39) C 1 − C 0 Δ t = − k 一个 v C 0 X 0 = X 1 − X 0 Δ t 。

平均衰减率 k_av在上面的方程确定从以下公式: (40) k av = k 0 + k 1 2 。

继续下一个时间步的过程 t₂= t₁+Δ t使用的新计算的值 X₁, C₁和 k₁计算的值 X₂, C₂和 k₂相应的时间步 t₂。这个过程继续进一步直到反应物和/或氯耗尽可用。正如上面所讨论的,的值 k₀和 一个 = X ¯ 一个 / X ¯ H 取决于校准的值给实验确定之间的最小误差平方和余氯浓度和计算氯残留根据上述步骤:

氯残留衰变率可以表示为一个一阶反应与整体反应速率 K,即 (41) K = k t X T 。

衰减率由以下表达式给出(一阶): (42) d C d t = d X T d t = − k t X T C = − K C 。

使用时间导数,得到 (43) d k t d t = d K / X T d t = d K / d T X T − d X T / d t K X T 2 , d k t d t = d K / X T d t = d K / d T X T − − K C K X T 2 。

进一步重新排列,得到 (44) d K d t = d k t d t X T − K 2 C X T 。

使用方程( 31日) dk _t/ dt,我们得到 (45) d K d t = − C k t 2 X ¯ 一个 X ¯ H − 1 X T − K 2 C X T , = − C K 2 X T 2 X ¯ 一个 X ¯ H − 1 X T − K 2 C X T , = − C K 2 X T X ¯ 一个 X ¯ H + K 2 C X T − K 2 C X T 。

最终, (46) d K d t = − C K 2 X T X ¯ 一个 X ¯ H 。

总体率包含类似的因素如方程( 31日)的初始剂量氯和算术的比例调和平均数的反应物的摩尔浓度。额外的因素是反应物的总摩尔浓度 X _T。更高的摩尔浓度降低了总体率导致陡峭的曲线接近一阶反应。降低总体摩尔浓度 X _T增加了衰变率接近二阶最终导致整体平坦衰减曲线的典型二阶衰减率。总体而言,以低剂量氯氯衰减率陡和更高的摩尔浓度的反应物和奉承相反在更高的初始氯剂量和降低反应物的摩尔浓度。因为总体率 K平均摩尔,这些值可以确定最初的初始氯剂量。

指的一阶方程提供的氯的衰减速度公式( 42),我们得到 (47) d C d t = K C , 在哪里 C的浓度 t, t天的时间通常转化为单位,然后呢 K总反应速率常数。

从方程( 46),我们得到 (48) d K d t = − C K 2 X T X ¯ 一个 X ¯ H 。

自从molar-averaged整体反应速率取决于反应物的初始摩尔浓度这意味着只考虑整体速度的偏导数与氯浓度 C, (49) ∂ K ∂ t = − α C K 2 , 在哪里 (50) α = X ¯ 一个 X ¯ H 1 X T = 常数 。

方程( 49)也可以写成 (51) ∂ 1 / K ∂ t = α C 。

以一阶导数对氯浓度,我们得到 (52) ∂ ∂ C ∂ 1 / K ∂ t = α 。

逆转分化的顺序,我们得到 (53) ∂ ∂ t ∂ 1 / K ∂ C = α 。

在任何给定的时间 t集成后,上面的表达式 (54) ∂ 1 / K ∂ C = α t = 常数 = K 0 , 或 (55) 1 K 2 ∂ K ∂ C = − X ¯ 一个 X ¯ H t X T = − K 0 , 自 (56) α = X ¯ 一个 X ¯ H 1 X T = 常数 。

最后, (57) ∂ K ∂ C = − K 0 K 2 。

初始速度之间的整合 C₀= 0,在任何给定的初始浓度 C₀给了 (58) ∫ β K d K K 2 = ∫ 0 C 0 − K 0 d C 0 , 在哪里 β是初始反应速率常数的初始浓度氯趋于零。

集成后,表情变得 (59) − 1 K − 1 β = − K 0 C 0 。

最后, (60) 1 K = 1 β + K 0 C 0 。

回归建模基于方程( 59)是由线性回归(1 / K)对初始浓度的氯衰减测试在不同初始浓度的氯。回归参数 β和 K₀决定从这一步。

;反应速率常数的表达式的回归参数确定后就变成了 (61) K = β 1 + β K 0 C 0 。

整个衰变率建模然后通过结合传统的第一顺序;还包括衰变速率与反应速率常数: (62) d C d t = − K C , C t = C 0 e − K t 。

用基于浓度的表达式 K价值在上面的方程,我们得到 (63) C t = C 0 e − β / 1 + β K 0 C 0 t 。

方程( 63年)可用于开发批量分销系统中的氯的衰减。水质建模项目比如用于EPANET可以适应总体率的变化对初始氯剂量考虑使用方程而不是在程序中使用的常量值。

瓶子测试实验大部分衰变氯残留在不同初始氯剂量下运行。测量的结果不同初始氯剂量下的氯在数据绘制 1- - - - - - 4。在每一批实验,衰变速率通常遵循一阶衰减率除了很高的初始氯(18.38 mg / L的初始剂量氯)显示早期的偏离线性衰减。一阶模型参数确定通过线性回归确定系数的指数衰减变化在0.89和0.95之间。图 5显示了实验运行的摘要氯残留在不同初始氯剂量使用。

反应速率常数的情节来自每个批量衰变实验策划针对相应的初始氯剂量使用,如图 6。从曲线图可以看出,之间存在反比关系大部分衰变反应速率常数和初始氯剂量用于实验。这种关系也类似于二阶衰变率数学模型,本文中开发和演示,最终显示在方程( 61年),描绘了逆氯衰变反应速率常数之间的关系 K和初始氯剂量使用。

曲线在图 6还表明,反应速率 K迅速减少,当初始氯剂量很低。这个地区氯剂量率限制氯氯衰减步反应是有限的。曲线的陡峭的斜坡在这个区域是一个指示的事实衰变氯反应是有限的。这个地区是一阶速率地区因为反应物浓度的衰减变化不明显反应进展有限的可用的氯。另一方面,作为初始氯剂量逐渐增加,衰减率曲线往往会奉承,证明弱初始氯剂和反应速率常数之间的关系。在这种情况下,反应反应物是有限的。这个区域是二阶速率区域因为反应物浓度变化明显随着衰变反应的进行。

二阶衰变率模型参数,即 β和 K₀所示的方程( 61年),是由回归方程的表达式(后 61年)是通过逆阴谋使直线化,即1 / K与初始氯剂量 C。结果数据图所示 7绘制一条直线与确定系数( R²)值为0.99。

为了能够使用数学模型的结果如上所述在方程( 31日)和不同反应条件下测试其有效性,它有助于运行一系列的实验在不同样品的水和确定实验常数。在这一过程中,进一步的必要与这些实验常数如总有机碳浓度和反应物特征紫外吸收,这样这些参数可以帮助设置在数学建模实验常数的值。

实验结果显示在本研究进行了一个相同的样本收集的水Matsapha水处理。因此,反应物的反应性物种的影响的模型参数方程( 63年),即 β和 K₀,不能完全确定,除非实验进行了广泛的样本。这是因为的价值 K₀取决于最易反应的反应物的反应在水里 β取决于如何使用快,和其他反应物的氯。为了能够使用方程( 63年)在水质模型如EPANET广泛的水质条件下,有必要再次关联的值 K₀和 β与反应物的化学物种参数如TOC和紫外线吸收,这样适当的值可以选择这些参数来建模的目的。能够做到这一点,建议基于初始氯试验进行了不同类型的水样在相同类型的反应物物种样本被测量适当的反应物的特征参数如TOC和紫外吸收。