-
地表水中有害藻类爆发(harmful algae bloom, HAB)引发的水体污染已经引起了全球公众的关注[1]。微囊藻是淡水生态系统中分布最普遍的藻类之一,其代谢产物——微囊藻毒素-LR(Microcystin-LR, MC-LR)会随着藻体的死亡和裂解而大量释放到水中[2]。MC-LR是强致癌物质,能够通过影响蛋白磷酸酶活性等途径破坏肝脏细胞,引发癌变。因此,MC-LR的存在会影响饮用水的水质安全,不利于渔业和城市居民的身心健康[3]。
去除饮用水中MC-LR的一种潜在解决方案是高级氧化工艺(AOP)。常规高级氧化技术包括紫外/过氧化氢(UV/H2O2)、紫外/次氯酸(UV/HClO)、H2O2/O3和基于Fenton的铁/H2O2体系可有效去除MC-LR[4]。相比于UV/H2O2体系,UV/HClO体系在pH低于7时的降解效率更高,因此,部分水处理厂已采用UV/HClO替代UV/H2O2[5]。此外,由于UV/HClO体系的余氯可以在饮用水处理中提供多种消毒屏障,从而简化工艺[6]。UV/HClO体系对MC-LR的降解途径包括:与HClO/ClO−的直接反应;通过紫外线辐射的直接光解;活性自由基(氯自由基(Cl·)和羟基自由基(HO·))等介导的间接光解。HO·能够以接近扩散的速率与污染物反应,达到降解的目的。水中的其他化合物(溶解的有机物(DOM)和碳酸氢盐)可以与污染物竞争HO·,但由于其具有快速反应活性,因此,依然可以用于多种类型难降解化合物的降解[7]。
Kintecus模型能够模拟高级氧化工艺反应,当前主要用于化学、生物、核和大气化学动力学及平衡过程的反应进行建模和回归/拟合/优化。一些研究[8-9]用该模型对污染物降解过程进行模拟,模拟得到的数据不仅可以用来支持和验证实验结果,也可以通过模拟实验来预测污染物的降解程度,论述建立的降解模型对实际水处理情况的适用性,并对运行参数的优化提供一定的参考。虽然利用kintecus模型模拟自由基稳态浓度的研究较多,但已有研究仅模拟了特定情境下UV/HClO体系降解污染物时自由基的变化规律,用以验证实验结果。然而,对于UV/HClO体系降解MC-LR过程中主要操作条件的模拟优化却鲜有报导。因此,本研究通过改变pH、次氯酸投加量和紫外光波长来模拟UV/HClO体系降解MC-LR的过程。利用已报导的实验数据来验证本研究建立动力学模型的准确性,并预测不同条件下UV/HClO体系对MC-LR的降解过程,优化UV/HClO体系参数。
全文HTML
-
在UV/HClO体系中,当HClO吸收光子<511 nm时,会产生羟基自由基(HO·)和氯自由基(Cl·)(式(1))[10],形成的Cl·会与HClO/ClO−溶液中的氯离子反应形成
Cl−2⋅ (式(2))[11-12]。HO·、Cl·和
Cl−2⋅ 均是强氧化剂,Cl·和Cl−2⋅ 标准还原电位分别为2.4 V和2.0 V,与HO·相当[13]。HO·、Cl·和Cl−2⋅ 通过单电子氧化,氢-吸附和不饱和C-C键的加成与有机物发生反应。Cl·是一种选择性氧化剂,在某些情况下其反应性可能高于HO·,但Cl−2⋅ 的反应性通常比HO·或Cl·低得多[5]。此外,自由基清除作用将会进一步使UV/HClO体系复杂化。例如HClO和ClO−会对Cl·和HO·具有不同程度的清除作用,Cl·会与OH−反应生成HClO−·,HO·会与Cl−反应也会形成HClO−·。其中HO·和Cl−之间的化学反应以很高的速率常数发生,但该反应是可逆的,并且逆反应的速率常数稍高(式(3))[12]。尽管Cl−会大量清除HO·,但HClO−·分解为HO·,会使HO·清除的影响极小,但Cl−的存在会引起体系内其他自由基稳态浓度的变化[14]。
碱度的增加会导致MC-LR降解速率的降低。含有碱度的MC-LR溶液在254 nm的吸光度显著降低[15]。在水溶液中存在
CO2−3 和HCO−3 ,当pH为6~8时,HCO−3 的比例更大[16]。式(4)和式(5)表明确实发生了HO·清除[17],但实际对于较低含量的自由氯体系,HCO−3 对体系中MC-LR的降解没有明显影响[16]。HO·与Cl·可能是反应体系中主要的污染物清除剂,这两者的生成速率与稳态浓度将成为UV/HClO体系去除污染物过程效率的限制性因素。因此,本次模型的建立主要关注HO·与Cl·在UV/HClO体系中的反应与转换。
-
Kintecus模型主要由模型方程和反应速率常数构成。模型方程由UV/HClO体系基本自由基转化方程[18]和MC-LR主要降解方程组成[9, 19],而光解速率常数由实验条件决定,在UV/HClO体系中,参与紫外光解的主要物质为HClO和ClO−,光解速率常数由式(6)计算[20]。
式中:
rHClO/ClO− 为UV光解过程中HClO或ClO−衰减的速率;Φ为HClO/ClO−光解的表观量子产率,mol·E−1;Ep为单位面积的光通量,mE·(s·cm2)−1;d是有效光程长度,cm;ε为HClO/ClO−的摩尔吸收系数,L·(mol·cm)−1;C为HClO或ClO−的初始浓度,mol·L−1。游离氯光解速率的差异可以用氯的形态和其光子反应性的差异(即摩尔吸收系数和量子产率)来解释。氯的形态与pH的大小有关,可以根据pKa进行计算。对于光子反应性,HClO(εHClO,254 nm= 58 L·(mol·cm)−1)和ClO−(
εClO−,254nm = 62 L·(mol·cm)−1)的摩尔吸收系数相似[21],另外也与其他报导测出的实验结果基本一致(εHClO,254 nm = 60 L·(mol·cm)−1;εClO−,254nm = 58 L·(mol·cm)−1)[22]。本研究中HClO与ClO−摩尔吸收系数均取60 L·(mol·cm)−1。尽管有许多研究在不同的HClO投加量下测出了不同的HClO/ClO−量子产率(表1),但HClO投加量浓度小于71 mg·L−1(1 mmol·L−1),HClO的量子产率就基本恒定在1.0±0.1。ClO−浓度在3.5~640 mg·L−1(0.05~9.0 mmol·L−1),ClO−的量子产率就基本恒定在0.9±0.1[23],本研究的HClO投加量与ClO−变化均在上述范围内,因此本研究选取量子产率ΦHClO=1,
ΦClO− =0.9。此外,根据文献的实验条件[27],确定其他实验参数:反应皿直径为5 cm;溶液体积为30 mL;溶液深度(有效光程长度d)为4.8 cm;辐照面积为6.25 cm2;光源强度(光通量)为9.31×10−9 E·s−1;HClO投加量为42 μmol·L−1;MC-LR浓度为1 μmol·L−1;实验温度为25 ℃。
1.1. UV/HClO体系对有机污染物的降解机理
1.2. kintecus模型的建立
-
pH通过影响HClO/ClO−的形态来决定氯降解的吸光度和量子产率,进而决定自由基的生成量[18, 24]。从理论上讲,次氯酸在pH为4.0~6.0时主要以次氯酸的形式存在;碱性条件(pH> 7.0)会导致次氯酸根离子占主导地位;而pH为1的溶液中次氯酸将转化为Cl2[10]。调节pH显著影响以HClO形式存在的氯的比例[28],以不同的试剂去调节pH会在不同程度上改变HClO与ClO−的比例,但在弱酸性条件下,HClO均为氯的主要存在形式[23]。本研究设置了一系列pH,并利用式(7)~式(8)计算不同pH下HClO和ClO−的浓度(pKa=7.5),并对HO·与Cl·的稳态浓度进行了计算,结果见表2、表3和图1。
式中:Ka为水解速率常数,取值为2.95×10−8。
根据表3及图1中模拟结果与原文献实验值的对比,在pH为6、7和9的条件下,HO·稳态浓度和Cl·稳态浓度误差均控制在1.3倍左右。在pH为8时,HO·稳态浓度的预测偏差略大,预测值为实验值的1.5倍左右;Cl·稳态浓度的实验值为预测值的1.5倍左右。模型预测值和实验值的误差主要来自于对参与该过程反应的不完整描述与光解速率常数计算参数的选择。在紫外辐照下MC-LR分子会裂解产生NH2·[29],NH2·经过复杂的反应会生成一系列含氮化合物(NH2·,NO·,
NO−2⋅ 等)[30],HO·会与含氮化合物发生反应导致HO·含量降低,但该类反应的研究目前尚处于起步阶段,因此,对HO·预测精准性尚有待深入研究[31]。此外,HClO/ClO−浓度是根据理想公式计算的,HClO/ClO−计算值会略大于实际含量,因此,在酸性条件下,HO·与Cl·的稳态浓度预测值高于实验值。而随着pH的增加,与实验值相比,模型预测的Cl·稳态浓度减小幅度更大,这可能是由于Cl−对Cl·清除作用(式(2))。在UV/HClO体系中Cl−主要由ClO−光解产生(式(9)),由于模型预测的ClO−的含量可能较高,且在中性和碱性条件下OCl−的含量高于酸性条件,因此,容易高估Cl−的含量,从而促进Cl·的转化,因此预测的Cl·的稳态浓度在pH为7和8时低于实验值(图1(b))。不同pH对MC-LR降解效果影响明显。如图2(a)所示,pH由6.0增加到9.0时,MC-LR降解的一级反应速率常数大大降低,该结果与原文献的实验结果一致[27]。说明当HClO为主要反应性物质时,反应体系能够生成足量的HO·和Cl·等活性物质,而ClO−的反应性可忽略不计,该模拟结果与先前的研究报导一致[19]。且在酸性溶液中反应速率的增加表明该反应易受反应体系的pH影响。因此,本研究模拟了pH为4和pH为11条件下的HO·和Cl·稳态浓度,以研究低pH与高pH条件对MC-LR降解过程的影响,发现HO·稳态浓度与Cl·稳态浓度均在pH为4时达到最高,在pH 为11时最低,且呈单调递减趋势(图1);当pH为4时,HO·和Cl·的稳态浓度与pH为6时模拟结果基本相同,且当pH由6降到4时,其表观降解速率常数仅由5.26×10−3 s−1提升到5.30×10−3 s−1,仅提升了0.76%,可以忽略不计。从降解效率和绿色环保的角度可以推测比较合适的pH为6。
-
HClO与MC-LR的摩尔比(nHClO∶nMC-LR)为12∶1时,MC-LR的氯化反应符合伪一阶反应[32]。在pH为7.78的条件下,将HClO与MC-LR的摩尔比为1∶1、2∶1和5∶1的样品避光氯化180 min后体系中的自由氯可被消耗殆尽,说明HClO投加量会对MC-LR的降解产生显著的影响[33]。在pH为7.78时,计算得到HClO和ClO−的含量分别为2.69×10−5 mol·L−1和1.51×10−5 mol·L−1,其他参数与上述讨论保持一致,计算光解速率常数后将参数代入模型计算得出自由基稳态浓度(表4)。
当MC-LR初始浓度始终为1 μmol·L−1时,随着HClO的增加,HClO的光解速率常数、自由基稳态浓度和降解速率均逐步升高,且HO·浓度的升高幅度大于Cl·(图3(a)、图3(b)和图4(a)),这更有利于污染物的降解。当HClO投加量低于12 μmol·L−1时,在进行60 min紫外辐照后MC-LR仍无法完全降解;但随着HClO投加量提升,MC-LR的降解率明显提高(图5(b)),在氯添加量为20 μmol·L−1时,在经过60 min降解后能够将MC-LR完全降解 (图4(b))。
MC-LR降解的表观速率常数与添加的HClO量呈正比的关系(图5(a))。在低HClO浓度下(<20 μmol·L−1),氯投加量增加能显著提升降解速率(图5(b)),氯添加量由12 μmol·L−1增加到20 μmol·L−1时,降解率达到90%所需的时间分别由42.9 min减少到25.2 min。在HClO投加量为30 μmol·L−1时,MC-LR降解率达到90%的时间为16.4 min,HClO投加量由42 μmol·L−1提升到60 μmol·L−1时,MC-LR降解率达到90%的时间仅由11.5 min变为8.5 min。HClO投加量大于40 μmol·L−1之后降解率趋向于稳定。考虑到HClO投加量与降解所需时间,根据模型结果推测HClO适宜投加量应为30~42 μmol·L−1。
-
光源的性质通过辐照光波长和光源强度来表现,不同性质的光源可以影响MC-LR降解效果。不同的波长会影响MC-LR对光的吸收程度以及体系内氧化剂的量子产率,为了准确计算光解速率常数,必须保证光源强度一致的前提下计算HClO/ClO−在对应辐照波长下的摩尔吸收系数以及量子产率。本研究列举了HClO/ClO−在波长257.7、268、282.3和301.2 nm下的摩尔吸收系数和量子产率(表5)。并通过固定光源强度[27]模拟不同紫外光波长对污染物降解的影响。
在弱酸性、中性和弱碱性3个状态下(pH=6、7和8)对4个不同波长(257.7、268、282.3和301.2 nm)光源条件下对初始浓度为1 μmol·L−1的MC-LR进行了降解模拟。结果表明,在pH为6时MC-LR具有最佳的降解效果,且pH的增加使MC-LR的表观降解速率常数降低(图6(a)),增加波长可以使MC-LR的表观降解速率常数下降(图6(b)),这可能是由于波长增加导致量子产率减小的缘故。当pH由6提升至8时,MC-LR的表观降解速率由5.07×10−3 s−1下降到3.84×10−3 s−1,减小了24%,当波长由257.7 nm提升到301.2 nm时,表观降解速率由5.07×10−3 s−1下降到4.69×10−3 s−1,下降了7.5%,pH变化引起的MC-LR的表观降解速率提升比波长变化的影响更加显著,可以推测pH才是MC-LR实验中最应优先考虑的操作因素。而波长变化引起的反应速率常数的变化不明显,可能是由于本次模拟所选取的波长变化范围较窄仅为43.5 nm所造成的。
当波长由257.7 nm增加到301.2 nm时,HClO的摩尔吸收系数由48.35 L·(mol·cm)−1下降到25.2 L·(mol·cm)−1为原来0.5倍,ClO−的摩尔系数从83.32 L·(mol·cm)−1增加到316.15 L·(mol·cm)−1提升了3.8倍(表5),增加波长ClO−的摩尔吸收系数增加的幅度显著大于HClO,且增加波长能大幅度提升ClO−光解速率(图7(a)),但MC-LR在碱性条件下的降解效率低于酸性条件,这归因于尽管在不同的波长下或在不同的pH下HO·均是主要的反应性物质(图7(b)),这与这与原文献实验研究结果一致[27]。在pH为6的条件下,在波长由257.7 nm增加到301.2 nm的过程中,HO·与Cl·的稳态浓度均随波长增加而降低,说明低波长下能产生更多的HO·与Cl·。
2.1. HO·与Cl·在不同pH下的稳态浓度
2.2. 次氯酸投加量对自由基稳态浓度和MC-LR降解率的影响
2.3. 紫外光波长对自由基稳态浓度和MC-LR降解率的影响
-
1) pH的模拟结果与实验结果相近,误差基本控制在1.5倍以内。当辐照波长为254 nm时,当体系为弱酸性环境(pH=6)时,对MC-LR的降解效果较好,继续降低pH,MC-LR的降解效果提升不明显。
2)在254 nm的辐照波长下,当MC-LR为1 μmol·L−1时,HClO投加量超过40 μmol·L−1之后,降解效率的提升趋势较为平缓,结果表明,HClO投加量在30~42 μmol·L−1时UV/HClO体系降解MC-LR效果最好。
3)在相同辐照功率下,在低波长(257.7~301.2 nm)的范围内,低波长具有更好的降解效果,但相对pH和HClO投加量,波长对MC-LR的降解效果影响并不明显。