Assessment of intervention measures on trihalomethane in finished water by interrupted time series analysis
-
摘要:背景
青草沙水库面临着水藻繁殖、富营养化的问题,由此引起水中的氯化消毒副产物含量增加一直是人们关注的重点。
目的评估“青草沙水库藻类监测与应对方案”(后称“方案”)的实施对常规处理出厂水中三卤甲烷(THMs)的治理效果。
方法从2011—2019年,每个季节采集陆家嘴水厂1次出厂水水样,每次采集1个样品,检测出厂水中四种THMs(三氯甲烷、二氯一溴甲烷、一氯二溴甲烷、三溴甲烷)的质量浓度(后称浓度)。以2014年“方案”的实施作为干预分界点,将整个研究阶段分为干预措施实施前(2011—2013年)和实施后(2014—2019年)两个阶段。应用中断时间序列分析方法进行分段线性回归,评估干预措施实施前后出厂水中THMs浓度和变化趋势。
结果在“方案”实施前,时间每延长一个季节,出厂水中总THMs浓度上升1.561 µg·L−1(
P =0.010),最高浓度达25.97 µg·L−1。干预措施实施后,出厂水中总THMs的浓度变化未见统计学意义(P =0.198);但出厂水中总THMs浓度回归线的斜率为-0.626(P =0.001),即THMs浓度随时间变化呈下降趋势。自2017年起至2019年,陆家嘴水厂出厂水中THMs的平均浓度保持在10 µg·L−1以下。结论青草沙水库藻类和富营养化治理措施达到了较好的效果,将出厂水中的THMs控制在较低水平,THMs的变化趋势由干预前的逐年上升,变为干预后的逐年降低。
Abstract:BackgroundThe Qingcaosha Reservoir is facing issues of algal blooms and eutrophication, and the resulting increase in the level of chlorination disinfection by-products in the water has been a major concern.
ObjectiveTo evaluate the impact of "Algae Monitoring and Control Program in Qingcaosha Reservoir" (hereinafter referred to as the program) on the control of trihalomethanes (THMs) in conventional finished water.
MethodsFrom 2011 to 2019, water samples were collected from the Lujiazui Water Plant once per season, one sample each time, and the concentrations of four THMs (trichloromethane, dichlorobromomethane, monochlorodibromomethane, and tribromomethane) were measured in the samples. Using 2014 when the program was implemented as a cut-off point, the entire study period was divided into two phases: pre-implementation (2011–2013) and post-implementation(2014–2019). Segmented linear regression with interrupted time series analysis was applied to assess the concentrations and trends of THMs in the finished water before and after the program launch.
ResultsThe concentration of total THMs in finished water increased by 1.561 µg·L−1 (
P =0.010) for each season of time extension before launching the program. The change in the concentration of total THMs in finished water was not statistically significant after the program launch, but the THMs concentration showed a decreasing trend as the slope was −0.626 (P =0.001). From 2017 until the end of 2019, the average concentration of THMs in finished water of Lujiazui Water Plant dropped to 10 μg·L−1 or less.ConclusionsThe algae and eutrophication control measures in Qingcaosha Reservoir have achieved good results, controlling THMs in finished water at a low level, and the trend of THMs has changed from a yearly increase pattern before the program to a yearly decrease pattern after the program.
-
消毒是饮用水水处理过程中不可或缺的工序,可以杀灭病原微生物,防止介水传染病的发生,保护公众健康,然而消毒过程中会不可避免地产生各种消毒副产物。常见的消毒副产物包括受控的三卤甲烷(trihalomethanes,THMs)、卤乙酸及非受控的卤乙腈、卤酮等。研究显示这些物质具有致癌性,美国环保局已将THMs中的三氯甲烷定义为B1组(人类证据有限的致癌物),将二氯一溴甲烷和三溴甲烷定义为B2组(动物实验证据充足的可能人类致癌物)[1]。
位于长江口的青草沙水库于2011年1月正式运营,作为水源地向上海市的陆家嘴自来水厂供应原水。由于该水库一直面临着水体富营养化、藻类繁殖的问题,导致出厂水中消毒副产物的含量增加[2]。为确保供水安全,控制自来水厂出厂水中以三卤甲烷为代表的消毒副产物的含量,上海市水务部门2013年制定了高峰供水期“青草沙水库藻类监测及应对方案”(后简称:“方案”),通过加强水质监控,优化水源地环境等措施对原水的富营养化问题进行治理[3]。中断时间序列分析(interrupted time series, ITS)经常用于评估在明确规定的时期内人群水平上采取的卫生干预措施的有效性[4–5]。在国外,ITS已被广泛用于评估公共卫生干预措施,包括新疫苗、自行车头盔立法、对乙酰氨基酚(扑热息痛)包装的变化、交通区间测速、预防鼻腔感染的措施,以及评估诸如全球金融危机等意外事件对健康的影响[6–11]。但ITS分析方法目前在国内公共卫生干预措施效果评价中的应用并不多见。
本研究旨在观察青草沙水库供水时期,出厂水中THMs的质量浓度(后称浓度),将2014年作为实施“方案”的干预分界点,分为“方案”实施前(2011—2013年)和实施后(2014—2019年)两个阶段,采用ITS方法分析“方案”对出厂水中THMs浓度的控制效果。
1. 材料与方法
1.1 样品采集
从2011—2019年,每季度采集1次陆家嘴水厂出厂水的水样,每年采集时间分别为2月、5月、8月、11月,研究期间共采集36个样品,每个样品采集60 mL。用60 mL的棕色螺纹口样品瓶(德国CNW)采集水样至满瓶(不应有气泡),样品在4 ℃冷藏保存。采样后加盐酸至样品PH小于2,样本尽快送至实验室,在24 h内完成测定。
1.2 样品处理和分析
采用USEPA 524.2吹脱捕集气质联用方法对样品进行检测[12]。将样品平衡至室温后,沿瓶壁倒满40 mL样品瓶,拧紧顶空瓶盖。利用吹扫捕集自动进样器精确移取5.0 mL水样至吹扫管,利用吹扫捕集装置,用氦气将水样中的挥发性有机物吹脱、捕集在捕集管内。吹脱程序完成后,将捕集管瞬间加热并以氦气反吹,将所吸附的目标化合物解吸进入毛细管气相色谱质谱联用仪进行分离测定,根据保留时间及特征离子定性,内标法定量。具体检测条件如下。吹扫捕集条件:捕集阱材料2,6-二苯基呋喃多孔聚合物树脂(Tenax)、硅胶(silicol)、椰壳碳(charcol)混合吸附剂(美国Tekmar公司),吹扫进样量5 mL,吹扫温度25 ℃,吹扫时间11 min,解吸温度220 ℃,解吸时间2 min,烘烤温度260 ℃,烘烤时间2 min。色谱条件:色谱柱HP-5MS UI (30 m × 0.25 mm,0.25 μm);进样口温度220 ℃;分流比为20∶1;程序升温为35 ℃保持5 min,然后以6 ℃·min−1升至100 ℃,再以20 ℃·min−1升至210 ℃,保持1 min。质谱条件:离子源类型EI,离子源温度230 ℃,传输线温度260 ℃,溶剂延迟1.6 min,全扫描方式,扫描范围为35~450 u。
为保证分析的准确性,每批样品测试空白样品,并且每20个样品运行一个质控样,控制相对偏差小于10%。本检测方法的检出限:0.50 µg·L−1;回收率:90%~110%;相对标准偏差:<5%;方法线性范围:0.50~25 µg·L−1。将未检出浓度设为检出限的1/2,即0.25 µg·L−1。
1.3 ITS回归模型构建
将水源地“方案”的实施作为一个干预事件,对中断时间序列数据进行分段线性回归,估计与干预措施实施效果相关的2个参数:(I)干预前后THMs的浓度变化;(II)干预引起的趋势或斜率的变化。干预前的时间段为2011—2013年,干预后的时间段为2014—2019年。拟合的水平和斜率改变回归模型如下:
$$ {Y}_{t}={\beta }_{0}+{\beta }_{1}T+{\beta }_{2}{X}_{t}+{\beta }_{3}\left(T-{T}_{0}\right)\cdot {X}_{t} $$ (1) 式(1)中:应变量Yt为每个季节出厂水中THMs的浓度;β0是时间序列开始前THMs的基线水平;β1是干预前THMs浓度随时间单位的变化趋势;β2是干预措施引起的浓度改变;β3是干预措施引起斜率的变化;(β1 + β3)是干预后THMs浓度随时间的变化趋势(干预后回归线的斜率);T代表时间单位,在本研究中以一个季度为一个时间单位;$ {T}_{0} $为干预措施开始的时间,在本研究中$ {T}_{0} $=12;Xt代表干预措施的不同阶段,本研究分干预前(Xt=0)和干预后(Xt=1)两个阶段。本研究基于12个干预前数据点(2011—2013年)和24个干预后的数据点(2014—2019年)进行ITS分析。利用Durbin-Watson(DW)法检验序列是否存在1阶自相关。利用增广Dickey-Fuller法检验干预前后序列的平稳性,发现数据达到开展时间序列分析的弱平稳性要求。利用赤池信息准则(Akaike information criterion, AIC)最小原则选择最佳回归模型,以确定纳入模型中的回归项,提高统计效力。通过计算模型预测值和真实值的差,绘制残差图,并计算序列数据的偏自相关系数,对模型的拟合效果进行检验。
1.4 统计学分析
采用Excel2000软件建立数据库,数据保存为.csv格式。采用R 3.5.2软件和Graphpad 5.0软件对数据进行统计分析和绘图。采用"lmtest","Epi","tsModel","vcd"等R包,对“方案”实施前后的监测数据进行ITS分析[13]。本研究建立ITS模型的R语言编码见链接https://blog.csdn.net/franksui/article/details/129819879。在ITS分析中,对回归模型的检验采用双侧检验,检验水准为α=0.05。
2. 结果
2.1 饮用水中三卤甲烷的浓度
2011—2019年陆家嘴水厂出厂水中THMs的平均浓度为14.31 µg·L−1。图1展示了2011—2019年逐年出厂水中THMs的平均浓度,可以发现2013年前THMs的浓度呈上升趋势,最高达25.97 µg·L−1;而自2014年起,浓度逐年降低,在2017年达到最低,为8.38 µg·L−1,至2019年底一直保持在10 µg·L−1以下。
干预措施实施前后出厂水中THMs的浓度见表1。“方案”实施后THMs的平均浓度降至12.51 µg·L−1,与“方案”实施前的浓度(17.92 µg·L−1)相比,下降了30.2%。
表 1 干预前后出厂水中THMs浓度的变化Table 1. Variation of THMs concentrations in finished water before and after interventionTHMs 平均浓度/(µg·L−1) 范围/(µg·L−1) 干预前 干预后 干预前 干预后 总THMs 17.92 12.51 4.37~31.70 3.98~34.36 三氯甲烷 7.09 3.95 0.25~17.70 0.25~15.20 一氯二溴甲烷 3.26 2.67 0.25~11.10 0.56~9.80 二氯一溴甲烷 6.43 4.84 0.25~11.00 1.50~11.80 三溴甲烷 1.17 0.80 0.25~6.67 0.25~4.11 2011—2019年每一季度出厂水中THM 浓度和建立模型所需变量的设置(Yt ,T,Xt)的数据(以2013、2014年数据为例)格式示意,见表2。
表 2 陆家嘴水厂出厂水中THMs浓度变化趋势中断时间序列模型构建的数据格式(以2013年、2014年数据为例)Table 2. Data format for interrupted time series modeling of the trend of THMs concentrations in finished water from Lujiazui Waterworks (taking 2013 and 2014 as examples)年份 THMs (Yt) 季节 时间单位(T) 干预阶段 (Xt) 2013 24.82 1 9 0 2013 28.23 2 10 0 2013 19.02 3 11 0 2013 31.70 4 12 0 2014 21.65 1 13 1 2014 34.36 2 14 1 2014 22.10 3 15 1 2014 15.65 4 16 1 [注] 季节:每年的季节编号从春夏秋冬依序为1~4;时间单位从2011年第一季度开始顺序编号为1~36;干预阶段:方案实施前设为0,方案实施后设为1。 2.2 “方案”实施后THMs的浓度变化趋势
在图2中用阴影的深浅来表示干预的不同时期。干预之前,在2011—2013年,出厂水中THMs的浓度逐年增加;对水库中的藻类开始干预和控制后,斜率由正值变为负值,出厂水中THMs的浓度呈现降低趋势。施行了水源藻类控制“方案”后,THMs浓度(深灰色区域内的蓝色实线)要低于未干预状态的预测值(深灰色区域内的红色虚线)。
ITS回归模型的各项统计学参数见表3。结果显示,在干预措施实施之前的阶段,出厂水中THMs的斜率为1.561,呈上升趋势,即随着时间因素每延长一个季节,出厂水中THMs浓度上升1.561 µg·L−1,有统计学意义(b1=1.561,P=0.010)。干预措施实施后,出厂水中THMs水平的改变没有统计学意义(b2=−6.180,P=0.198)。干预措施导致的出厂水中THMs斜率的改变(b3=−2.187),具有统计学意义(P=0.001)。干预措施实施以后的阶段,THMs的斜率=b1+ b3,即1.561 + (−2.187)=−0.626,即THMs浓度随时间变化呈下降趋势。
表 3 ITS回归模型参数Table 3. Parameters of ITS regression model变量 回归系数(b) t P 常数项 7.736 1.846 0.074 时间($ T $) 1.561 2.728 0.010 干预($ {X}_{t} $) −6.180 −1.314 0.198 斜率改变($ \left(T-{T}_{0}\right)\cdot {X}_{t} $) −2.187 −3.604 0.001 3. 讨论
出厂水中的THMs是饮水消毒过程中产生的主要氯化消毒副产物,包括三氯甲烷(CHCl3)、一氯二溴甲烷(CHBr2Cl)、二氯一溴甲烷(CHBrCl2)、三溴甲烷(CHBr3)。Richardson等[14]的研究表明,终生暴露于饮用水消毒副产物会增加癌症风险,尤其是膀胱癌和结肠癌。另外,THMs对人体的肝脏、肾脏、神经系统、生殖系统等器官具有非致癌风险[15–17]。陆家嘴水厂采用的是“絮凝-沉淀-过滤-消毒”的常规水处理工艺,加氯消毒过程主要应用氯和氯胺作为消毒剂。本研究对2011—2019年陆家嘴水厂的出厂水中的THMs的浓度进行了连续监测,过去9年间,THMs的平均浓度为14.31 µg·L−1,低于伊朗Bushehr省(92.9 µg·L−1)[18]、中国北方的山东省(35 µg·L−1)[19]出厂水中的检出浓度,与中国南方的珠江三角洲出厂水中的水平(17.7µg·L−1)接近[20]。在本研究中,该水厂所有出厂水样品中THMs的总浓度远低于WHO和欧盟的限值(100 µg·L−1)以及US EPA的限值(80 µg·L−1)[21]。
当青草沙水库自2011年建成开始供水后,它供应了上海市所有自来水厂原水总用量的60%,覆盖人口约1000万。青草沙水库作为上海市饮用水的新水源地,具有水量充足,水质优良稳定,水源防护与运行成本低等优势。然而,有学者对青草沙水库的水质进行了调查,发现水库内氮/磷比值较高(16~29 : 1),而且夏季水库蓝藻、微囊藻数量明显增加,这表明青草沙水库面临着水体富营养化和蓝藻繁殖的问题[22]。水体中的藻类繁殖会为THMs的形成提供了更多的有机前体物质(腐殖酸、富里酸、氨基酸),饮用水加氯消毒后,出厂水中的THMs含量增加[15,23]。
在上海市水务管理部门制定的青草沙高峰供水时期的《藻类监测与应对方案》中,主要采取对水体中氮/磷营养盐、藻类进行环境监测,利用水动力调控、物理拦截等控藻技术和生态工程联合治理技术,优化水源地周边环境等水体富营养化综合治理措施[24–25]。陆家嘴水厂自2010年底开始由青草沙水库供水,截至2019年底,出厂水中THMs的平均浓度降至10 µg·L−1以下,与高峰年份的25.97 µg·L−1相比,降低了67.73%。在2014年的第二季度出厂水中THMs的浓度达到最大值(34.36 µg·L−1),随后浓度开始逐年下降。由ITS分析的结果可见,干预措施不仅引起了出水中THMs浓度的降低,而且扭转了THMs浓度的变化趋势,即由干预前的逐年上升趋势转变为干预后的逐年下降趋势。这一趋势的变化表明,干预措施有效地遏制了原水中藻类迅速增长的态势,从源头上降低了形成消毒副产物所需要的前体物质,有效降低了水中消毒副产物的含量。今后,饮用水监管部门应继续加强水源地的保护和监管,减少水源地上游区域含磷、含氮元素废水的排放量,研发富营养化原水深度处理新工艺,从而进一步降低饮用水中消毒副产物含量,提高饮用水的卫生安全性。
本研究具有一定的局限性。首先,采样点位于浦东新区的一家水厂,一个采样点的监测数据无法代表整个上海市的消毒副产物的水平。其次是监测频率低,每3个月采样监测一次,无法反映水中污染物逐月的动态变化,季节性的变化也不明显,后续的研究需要提高监测频率,如果能做到每月采样检测一次,将会有效地弥补这一缺陷。另外,适当增加样本量、平行样也是必要的,以减少抽样误差。
综上所述,青草沙水库供水高峰期《藻类监测及应对方案》实施后,陆家嘴水厂的出厂水中THMs的浓度自2014年起明显下降,THMs浓度的发展趋势由干预前的逐年上升,变为干预后的逐年下降,说明上海市水务管理部门采取的藻类监测与原水富营养化治理措施取得了较好的效果。
-
表 1 干预前后出厂水中THMs浓度的变化
Table 1 Variation of THMs concentrations in finished water before and after intervention
THMs 平均浓度/(µg·L−1) 范围/(µg·L−1) 干预前 干预后 干预前 干预后 总THMs 17.92 12.51 4.37~31.70 3.98~34.36 三氯甲烷 7.09 3.95 0.25~17.70 0.25~15.20 一氯二溴甲烷 3.26 2.67 0.25~11.10 0.56~9.80 二氯一溴甲烷 6.43 4.84 0.25~11.00 1.50~11.80 三溴甲烷 1.17 0.80 0.25~6.67 0.25~4.11 表 2 陆家嘴水厂出厂水中THMs浓度变化趋势中断时间序列模型构建的数据格式(以2013年、2014年数据为例)
Table 2 Data format for interrupted time series modeling of the trend of THMs concentrations in finished water from Lujiazui Waterworks (taking 2013 and 2014 as examples)
年份 THMs (Yt) 季节 时间单位(T) 干预阶段 (Xt) 2013 24.82 1 9 0 2013 28.23 2 10 0 2013 19.02 3 11 0 2013 31.70 4 12 0 2014 21.65 1 13 1 2014 34.36 2 14 1 2014 22.10 3 15 1 2014 15.65 4 16 1 [注] 季节:每年的季节编号从春夏秋冬依序为1~4;时间单位从2011年第一季度开始顺序编号为1~36;干预阶段:方案实施前设为0,方案实施后设为1。 表 3 ITS回归模型参数
Table 3 Parameters of ITS regression model
变量 回归系数(b) t P 常数项 7.736 1.846 0.074 时间($ T $) 1.561 2.728 0.010 干预($ {X}_{t} $) −6.180 −1.314 0.198 斜率改变($ \left(T-{T}_{0}\right)\cdot {X}_{t} $) −2.187 −3.604 0.001 -
[1] 李永珍, 何更生, 詹铭, 等. 上海市水源水及出厂水中卤甲烷、卤乙酸、卤乙腈类消毒副产物含量及健康风险评估[J]. 环境与职业医学, 2021, 38(5): 460-466. LI Y Z, HE G S, ZHAN M, et al. Levels and health risk assessments of halomethanes, haloacetic acids, and haloacetonitriles disinfection by-products in source and finished water in Shanghai[J]. J Environ Occup Med, 2021, 38(5): 460-466.
[2] OU H S, WEI C H, DENG Y, et al. Principal component analysis to assess the composition and fate of impurities in a large river-embedded reservoir: Qingcaosha reservoir[J]. Environ Sci Process Impacts, 2013, 15(8): 1613-1621. doi: 10.1039/c3em00154g
[3] 上海水资源保护基金会. 2013年高峰供水青草沙藻类监测工作启动[EB/OL]. [2023-06-06]. http://www.swrf.org.cn/news-detail.asp?nid=4444. Shanghai Water Resources Foundation. In 2013, the algae monitoring work of Qingcaosha in peak water supply was started[EB/OL]. [2023-06-06]. http://www.swrf.org.cn/news-detail.asp?nid=4444.
[4] SOUMERAI S B, STARR D, MAJUMDAR S R. How do you know which health care effectiveness research you can trust? A guide to study design for the perplexed[J]. Prev Chronic Dis, 2015, 12: E101.
[5] WAGNER A K, SOUMERAI S B, ZHANG F, et al. Segmented regression analysis of interrupted time series studies in medication use research[J]. J Clin Pharm Ther, 2002, 27(4): 299-309. doi: 10.1046/j.1365-2710.2002.00430.x
[6] LAU W C Y, MURRAY M, EL-TURKI A, et al. Impact of pneumococcal conjugate vaccines on childhood otitis media in the United Kingdom[J]. Vaccine, 2015, 33(39): 5072-5079. doi: 10.1016/j.vaccine.2015.08.022
[7] DENNIS J, RAMSAY T, TURGEON A F, et al. Helmet legislation and admissions to hospital for cycling related head injuries in Canadian provinces and territories: interrupted time series analysis[J]. BMJ, 2013, 346: f2674. doi: 10.1136/bmj.f2674
[8] HAWTON K, BERGEN H, SIMKIN S, et al. Long term effect of reduced pack sizes of paracetamol on poisoning deaths and liver transplant activity in England and Wales: interrupted time series analyses[J]. BMJ, 2013, 346: f403. doi: 10.1136/bmj.f403
[9] GRUNDY C, STEINBACH R, EDWARDS P, et al. Effect of 20 mph traffic speed zones on road injuries in London, 1986-2006: controlled interrupted time series analysis[J]. BMJ, 2009, 339: b4469. doi: 10.1136/bmj.b4469
[10] DERDE L P G, COOPER B S, GOOSSENS H, et al. Interventions to reduce colonisation and transmission of antimicrobial-resistant bacteria in intensive care units: an interrupted time series study and cluster randomised trial[J]. Lancet Infect Dis, 2014, 14(1): 31-39. doi: 10.1016/S1473-3099(13)70295-0
[11] LOPEZ BERNAL J A, GASPARRINI A, ARTUNDO C M, et al. The effect of the late 2000s financial crisis on suicides in Spain: an interrupted time-series analysis[J]. Eur J Public Health, 2013, 23(5): 732-736. doi: 10.1093/eurpub/ckt083
[12] ALFORD-STEVENS A, EICHELBERGER J W, BUDDE W L. Measurement of purgeable organic compounds in water by capillary column gas chromatography/mass spectrometry: EPA Method 524.2[R]. Cincinnati: U. S. Environmental Protection Agency, 1992.
[13] BERNAL J L, CUMMINS S, GASPARRINI A. Corrigendum to: interrupted time series regression for the evaluation of public health interventions: a tutorial[J]. Int J Epidemiol, 2020, 49(4): 1414.
[14] RICHARDSON S D, PLEWA M J, WAGNER E D, et al. Occurrence, genotoxicity, and carcinogenicity of regulated and emerging disinfection by-products in drinking water: a review and roadmap for research[J]. Mutat Res, 2007, 636(1/3): 178-242.
[15] RISTOIU D, VON GUNTEN U, MOCAN A, et al. Trihalomethane formation during water disinfection in four water supplies in the Somes river basin in Romania[J]. Environ Sci Pollut Res Int, 2009, 16 Suppl 1: S55-S65.
[16] SADIQ R, RODRIGUEZ M J. Disinfection by-products (DBPs) in drinking water and predictive models for their occurrence: a review[J]. Sci Total Environ, 2004, 321(1/3): 21-46.
[17] WOO Y T, LAI D, MCLAIN J L, et al. Use of mechanism-based structure-activity relationships analysis in carcinogenic potential ranking for drinking water disinfection by-products[J]. Environ Health Perspect, 2002, 110 Suppl 1(Suppl 1): 75-87.
[18] DOBARADARAN S, SHABANKAREH FARD E, TEKLE-RÖTTERING A, et al. Age-sex specific and cause-specific health risk and burden of disease induced by exposure to trihalomethanes (THMs) and haloacetic acids (HAAs) from drinking water: an assessment in four urban communities of Bushehr province, Iran, 2017[J]. Environ Res, 2020, 182: 109062. doi: 10.1016/j.envres.2019.109062
[19] YAO Z, SUN S, WANG M, et al. The occurrence of THMs and AOX in drinking water of Shandong province, China[J]. Environ Sci Pollut Res Int, 2019, 26(18): 18583-18592. doi: 10.1007/s11356-019-05094-1
[20] GAN W, GUO W, MO J, et al. The occurrence of disinfection by-products in municipal drinking water in China's Pearl River Delta and a multipathway cancer risk assessment[J]. Sci Total Environ, 2013, 447: 108-115. doi: 10.1016/j.scitotenv.2012.12.091
[21] CHOWDHURY S, RODRIGUEZ M J, SADIQ R. Disinfection byproducts in Canadian provinces: associated cancer risks and medical expenses[J]. J Hazard Mater, 2011, 187(1/3): 574-584.
[22] HUANG Z, XIE B, YUAN Q, et al. Microbial community study in newly established Qingcaosha Reservoir of Shanghai, China[J]. Appl Microbiol Biotechnol, 2014, 98(23): 9849-9858. doi: 10.1007/s00253-014-5928-8
[23] GALLARD H, VON GUNTEN U. Chlorination of natural organic matter: kinetics of chlorination and of THM formation[J]. Water Res, 2002, 36(1): 65-74. doi: 10.1016/S0043-1354(01)00187-7
[24] 余梦. 青草沙监控水质严防藻类暴发[EB/OL]. [2023-06-06]. https://news.sina.com.cn/o/2013-07-02/100027553516.shtml. YU M. Qingcaosha monitors water quality to prevent algae outbreaks[EB/OL]. [2023-06-06]. https://news.sina.com.cn/o/2013-07-02/100027553516.shtml.
[25] 申一尘, 王绍祥, 朱宜平, 等. 青草沙水库藻类防控与多水源调配示范研究成果[J]. 净水技术, 2017, 36(6): 1-5. SHEN Y C, WANG S X, ZHU Y P, et al. Research achievements on algae control of Qingcaosha Reservoir and demonstration of multi-water sources distribution[J]. Water Purif Technol, 2017, 36(6): 1-5.