黃 輝, 周衛(wèi)華, 程志鵬, 李 倩, 費相琴, 徐繼明
(1.淮陰師范學院 化學化工學院 江蘇省生物質(zhì)能與酶技術(shù)重點實驗室, 江蘇 淮安 223300;2.淮安市環(huán)境監(jiān)測中心站, 江蘇 淮安 223001; 3.淮安自來水有限公司, 江蘇 淮安 223002)
常規(guī)氯化消毒工藝中典型消毒副產(chǎn)物生成的數(shù)學模型
黃 輝1, 周衛(wèi)華2, 程志鵬1, 李 倩1, 費相琴3, 徐繼明1
(1.淮陰師范學院 化學化工學院 江蘇省生物質(zhì)能與酶技術(shù)重點實驗室, 江蘇 淮安 223300;2.淮安市環(huán)境監(jiān)測中心站, 江蘇 淮安 223001; 3.淮安自來水有限公司, 江蘇 淮安 223002)
快速有效監(jiān)控飲用水中消毒副產(chǎn)物的生成,是飲用水安全管理中的重要基礎(chǔ)性工作.基于供水企業(yè)常規(guī)氯化消毒工藝過程,采用均勻試驗設(shè)計,研究了水源地源水UV254值、氨氮含量、加氯量和pH值等關(guān)鍵因素對典型消毒副產(chǎn)物三鹵甲烷(THMs)和鹵乙酸(HAAs)生成的影響,并以SPSS統(tǒng)計軟件對數(shù)據(jù)進行了線性和非線性擬合.結(jié)果表明:線性模型能較好地描述THMs和HAAs的生成(相關(guān)度均大于0.9),而非線性模型中二者的相關(guān)度分別為0.924和0.883;考慮各因素交互作用下,THMs和HAAs生成的線性模型的相關(guān)度均為1,進一步改善了擬合效果.建立的快速預測上述消毒副產(chǎn)物生成的新方法,可為完善檢測條件制約下飲用水安全風險管理提供理論依據(jù).
氯化消毒; 消毒副產(chǎn)物; 數(shù)學模型; 均勻設(shè)計
當前我國大多數(shù)自來水廠仍采用以去除源水中懸浮物質(zhì)和病源微生物為主要目的的混凝、沉淀、過濾、液氯消毒四段式傳統(tǒng)凈水工藝.液氯消毒是一種比較安全、簡便而又廉價的消毒方法,它能有效殺滅水中的微生物病原體,大大降低了人們因飲水而感染痢疾、霍亂等水傳播疾病致死的幾率,因而液氯消毒一直是國內(nèi)外城市給水處理中廣泛采用的消毒手段,在過去很長一段時間內(nèi),它在保障人體健康方面起了巨大作用[1,2].然而,出廠水在漫長的輸送過程中,其中所含的有機污染物與管網(wǎng)中余氯產(chǎn)生取代反應(yīng)生成多種具有“三致”作用的消毒副產(chǎn)物,致使自來水的致突變性比源水更強[2,3].另一方面,隨著城鎮(zhèn)化的快速推進,我國飲用水源水的污染態(tài)勢依然嚴峻.飲用水安全已成為當前政府、社會、公眾日益關(guān)注的焦點問題.
近年來,作為飲用水安全評價的重要指標——“三致”消毒副產(chǎn)物受到人們廣泛關(guān)注..其中,三鹵甲烷(THMs)和鹵乙酸(HAAs)是最典型也是研究最多的鹵代有機消毒副產(chǎn)物[2-5].我國現(xiàn)施行的《生活飲用水衛(wèi)生標準》(GB5749-2006)中明確規(guī)定了6種上述消毒副產(chǎn)物的監(jiān)測限值和要求——全部4種THMs(即CHCl3、CHBrCl2、CHBr2Cl和CHBr3)和HAAs中的二氯乙酸(DCAA)、三氯乙酸(TCAA),其中僅CHCl3是常規(guī)指標,其余均為非常規(guī)指標.受現(xiàn)階段我國大部分地區(qū)飲用水鹵代有機消毒副產(chǎn)物檢測條件的制約,其實施的有效性有待進一步提高,飲用水安全的風險存在被低估的可能.如何實現(xiàn)對飲用水消毒副產(chǎn)物的有效監(jiān)控,降低此類“隱形殺手”對人們健康的威脅,是當前及今后一段時期重要而迫切的課題.在此背景下,基于供水企業(yè)常規(guī)氯化消毒工藝過程,研究水源地源水生成上述兩種典型消毒副產(chǎn)物的數(shù)學模型,可為科學預測不同工況條件下消毒副產(chǎn)物的生成量,實現(xiàn)對飲用水中消毒副產(chǎn)物的快速監(jiān)控、評價和飲用水安全管理提供方法支撐和理論依據(jù).可為現(xiàn)有工藝提標改造或新工藝實施下飲用水中消毒副產(chǎn)物的快速監(jiān)控提供參考和借鑒.
在常規(guī)氯化消毒工藝中,影響THMs和HAAs生成的因素主要有源水水質(zhì)(水溫、TOC、UV254、NH3-N、藻類等)和工藝參數(shù)(氯的投加量、反應(yīng)時間、pH值、共存金屬離子等)[3-5],由于水溫在特定的供水期間變化不大,而藻類和共存離子在經(jīng)混凝、沉淀和過濾(及活性炭吸附)工藝處理后殘余濃度極低[2,6-8],同時考慮到水廠對水樣測試的簡便化、快速化要求,以UV254表征源水有機物濃度.因此在本研究中,僅考慮水樣UV254值、氨氮含量、加氯量和pH值四個關(guān)鍵因素,源水經(jīng)充分氯化反應(yīng)后,監(jiān)測THMs和HAAs的變化情況.為保證各因素各水平分布均勻且能減少試驗次數(shù),本研究選用均勻設(shè)計法開展實驗.
均勻設(shè)計(Uniform Design)[9]是中國數(shù)學家方開泰和王元于1978年首先提出來的,它是一種只考慮實驗點在實驗范圍內(nèi)均勻散布的實驗設(shè)計方法.該設(shè)計方法的優(yōu)勢在于設(shè)計較少次數(shù)的實驗就可取得較充分的信息,尤其適用于多因素、多水平的實驗設(shè)計,近年來獲得廣泛應(yīng)用[10].
由于均勻設(shè)計因素水平的選取不具有整齊可比性,無法通過指標觀測值直接得到指標與因素之間的關(guān)系,因此需要對數(shù)據(jù)進行回歸擬合.通過對擬合公式的系統(tǒng)分析,不僅可以定性地得到各因素對指標的主效應(yīng)以及各因素之間的交互效應(yīng),而且可以定量地求得達到最優(yōu)響應(yīng)值時的實驗條件.
2.1 均勻?qū)嶒炘O(shè)計
表1 均勻試驗工況設(shè)計表
注: 表格中括號前面數(shù)值表示該因素對應(yīng)的水平,括號內(nèi)表示其具體數(shù)值.
2.2 消毒副產(chǎn)物檢測方法
表2 均勻試驗工況條件下消毒副產(chǎn)物的生成量
THMs采用頂空氣相色譜法(APHA 5710 B)[11]測定,HAAs采用本研究組建立的液-液萃取、衍生化氣相色譜法[12]測定.色譜儀為惠普公司HP-6890型.
按照表1所設(shè)計的均勻試驗工況開展實驗,結(jié)果見表2.根據(jù)所得試驗數(shù)據(jù),采用SPSS18.0進行回歸分析.
3.1 線性擬合結(jié)果
假設(shè)消毒副產(chǎn)物生成量與UV254值、氨氮含量(以NH3-N表示)、加氯量(以Cl2表示)和pH值符合以下模型:
CDBP=k0+k1×[UV254]+k2×[NH3-N]+k3×[Cl2]+k4×[pH]
(1)
經(jīng)SPSS18.0線性擬合后,置信度在95%以上時模型中各參數(shù)值見表3.SPSS在進行回歸時,首先對所有的自變量進行了檢驗,判斷是否對因變量有影響且影響是否足夠大,之后才進行參數(shù)的確定.這樣影響小的變量就被剔除了,表現(xiàn)為模型中變量的系數(shù)為0(如表3中的k1).這說明,UV254值對消毒副產(chǎn)物生成的影響很小,這可能是由于在均勻?qū)嶒炛?UV254數(shù)值的變化幅度過小(原水UV254值較低)所致.從其余參數(shù)來看,氨氮含量和pH值與消毒副產(chǎn)物呈負相關(guān),且對THMs的影響高于對HAAs的影響;加氯量與消毒副產(chǎn)物呈正相關(guān).
表3 消毒副產(chǎn)物生成的線性模型參數(shù)
由表3可知,THMs和HAAs為因變量時線性模型的相關(guān)度分別為0.911和0.982,表明線性模型能較好地描述消毒副產(chǎn)物的生成.將表3中各參數(shù)代入式(1)可得:
CTHMs=47.189-3.545×[NH3-N]+3.861×[Cl2]-5.307×[pH]
(2)
CHAAs=13.980-1.179×[NH3-N]+1.217×[Cl2]-1.325×[pH]
(3)
3.2 非線性擬合結(jié)果
假設(shè)消毒副產(chǎn)物生成量與UV254、加氯量(以Cl2表示)和pH值符合以下模型:
CDBP=k0×[UV254]k1×[NH3-N]k2×[Cl2]k3×[pH]k4
(4)
兩邊取對數(shù)
lgCDBP=lgk0+k1×lg[UV254]+k2×lg[NH3-N]+k3×lg[Cl2]+k4×lg[pH]
(5)
以lgCDBP為因變量,以lg[UV254]、lg[NH3-N]、lg[Cl2]、lg[pH]為因變量,經(jīng)SPSS18.0擬合回歸后,置信度在95%以上時模型中各參數(shù)值見表4.
表4 消毒副產(chǎn)物生成的非線性模型參數(shù)
由表4可知, THMs和HAAs生成的非線性模型的相關(guān)度分別為0.924和0.883, 總體上低于線性模型.將表4中相關(guān)參數(shù)代入式(4)可得消毒副產(chǎn)物生成的非線性模型為:
CTHMs=101.735×[UV254]0.293×[NH3-N]-0.087×[Cl2]0.416×[pH]-1.106
(6)
CHAAs=100.048×[UV254]1.053×[NH3-N]0.001×[Cl2]0.377×[pH]0.077
(7)
3.3 包含各因素交互作用的線性擬合結(jié)果
進一步考慮各因素間存在交互作用下的線性擬合結(jié)果,假設(shè)消毒副產(chǎn)物生成量符合以下模型:
CDBP=k0+k1×[UV254]+k2×[NH3-N]+k3×[Cl2]+k4×[pH]+k5×[UV254]·[NH3-N]+
k6×[UV254]·[Cl2]+k7×[UV254]·[pH]+k8×[NH3-N]·[Cl2]+k9×[NH3-N]·[pH]+
k10×[Cl2]·[pH]+k11×[UV254]·[NH3-N]·[Cl2]+k12×[UV254]·[NH3-N]·[pH]+
k13×[UV254]·[Cl2]·[pH]+k14×[NH3-N]·[Cl2]·[pH]+
k15×[UV254]·[NH3-N]·[Cl2]·[pH]
(8)
經(jīng)SPSS18.0線性擬合后,置信度在95%以上時模型中各參數(shù)值見表5.
表5包含各因素交互作用的線性模型參數(shù)
由表5可知,擬合過程中除變量[NH3-N]、[Cl2]、[pH]、[UV254]·[pH]和[UV254]·[NH3-N]·[Cl2]·[pH]以外,其余變量均被剔除,表明其對消毒副產(chǎn)物生成的影響很小.考慮各因素交互作用下,THMs和HAAs生成的線性模型的相關(guān)度均為1,取得了極佳的擬合效果.影響THMs生成的因素排序為:氨氮含量>pH值(負相關(guān))>加氯量>UV254和pH的交互作用(負相關(guān))>四因素的交互作用(負相關(guān));影響HAAs生成的因素排序為:氨氮含量(負相關(guān))>pH值(負相關(guān))>加氯量>四因素的交互作用>UV254和pH的交互作用(負相關(guān)).
將表5中相關(guān)參數(shù)代入式(8)可得包含各因素交互作用的消毒副產(chǎn)物生成的線性模型為:
CTHMs=73.405+13.109×[NH3-N]+4.808×[Cl2]-5.439×[pH]-
1.177×[UV254]·[pH]-0.189×[UV254]·[NH3-N]·[Cl2]·[pH]
(9)
CHAAs=17.305-13.222×[NH3-N]+0.532×[Cl2]-1.401×[pH]-
0.006×[UV254]·[pH]+0.137×[UV254]·[NH3-N]·[Cl2]·[pH]
(10)
采用均勻試驗設(shè)計,討論了水樣UV254值、氨氮含量、加氯量和pH值等關(guān)鍵因素對消毒副產(chǎn)物THMs和HAAs生成的影響,并以SPSS統(tǒng)計軟件進行了線性及非線性擬合,結(jié)果表明:
1) 線性模型能較好地描述THMs和HAAs的生成(相關(guān)度均大于0.9),具體關(guān)系式為:
CTHMs=47.189-3.545×[NH3-N]+3.861×[Cl2]-5.307×[pH]
CHAAs=13.980-1.179×[NH3-N]+1.217×[Cl2]-1.325×[pH]
2) THMs和HAAs生成的非線性模型的相關(guān)度分別為0.924和0.883,總體上低于線性模型.具體關(guān)系式為:
CTHMs=101.735×[UV254]0.293×[NH3-N]-0.087×[Cl2]0.416×[pH]-1.106
CHAAs=100.048×[UV254]1.053×[NH3-N]0.001×[Cl2]0.377×[pH]0.077
3) 考慮各因素交互作用下,THMs和HAAs生成的線性模型的相關(guān)度均為1,進一步改善了擬合效果.具體關(guān)系式為:
CTHMs=73.405+13.109×[NH3-N]+4.808×[Cl2]-5.439×[pH]-
1.177×[UV254]·[pH]-0.189×[UV254]·[NH3-N]·[Cl2]·[pH]
CHAAs=17.305-13.222×[NH3-N]+0.532×[Cl2]-1.401×[pH]-
0.006×[UV254]·[pH]+0.137×[UV254]·[NH3-N]·[Cl2]·[pH]
由于試驗數(shù)據(jù)有限,上述模型還需進行進一步實踐驗證以及修正.
[1] 鄂學禮,王麗,邢方瀟.飲水消毒副產(chǎn)物及其標準研究進展[J].環(huán)境與健康雜志,2010, 27(1):2-4.
[2] Banna M H,Imran S,Francisque A, et al. Online drinking water quality monitoring: review on available and emerging technologies[J].Critical Reviews in Environmental Science and Technology,2013,106:389-395.
[3] Farre M J, Day S, Neale P A, et al. Bioanalytical and chemical assessment of the disinfection by-product formation potential: Role of organic matter[J].Water Research,2013,47:5409-5421.
[4] Jeyong Yoon,Youshik Choi,Soonhang Cho,et al. Low trihalomethane formation in Korean drinking water[J].Science of the Total Environment,2003,302:157-166.
[5] Villanueva C M,Kogevinasa M, Grimaltb J O. Haloacetic acids and trihalomethanes in finished drinking waters from heterogeneous sources[J]. Water research,2003,37:953-958.
[6] Richard A L,Ariene L R. Chloroform and chlorophenol production by decarboxylation of natural acids during aqueous chlorination[J].Environmental Science and Technology,1979,13(3):325-329.
[7] Boyce S D,Horning J F. Reaction pathways of trihalomethan formation from the halogenation of dihydroxyaromatic model compounds for humic acid[J]. Environmental Science andTechnology,1983,17(4):202-211.
[8] White D M,Garland D S,Narr J,et al.Natural organic mater and DBPs formation potential in Alaskan water supplies[J].Water research,2003,37(4):939-947.
[9] 方開泰.均勻設(shè)計與均勻設(shè)計表[M].北京:科學出版社,1994.
[10] 石瑞,史垠垠,曹福亮.均勻設(shè)計優(yōu)化渾濁型白果飲料中復配穩(wěn)定劑[J].食品工業(yè)科技,2012,33(19):222-223.
[11] American Public Health Association. Standard methods forthe examination of water andwastewater[M].NewyoK:19th Washington D C,1995.
[12] 黃輝,林芳,程志鵬,等.飲用水中鹵乙酸衍生-氣相色譜測定法改進[J].環(huán)境與健康雜志,2013,30(2):150-153.
MathematicalModelfortheFormationofTypicalDisinfectionby-ProductsduringConventionalChlorinationProcess
HUANG Hui1, ZHOU Wei-hua2, CHENG Zhi-peng1, LI Qian1, FEI Xiang-qin3, XU Ji-ming1
(1.School of Chemistry and Chemical Engineering of Huaiyin Normal University, Jiangsu Key Laboratory for Biomass based Energy and Enzyme Technology, Huaian Jiangsu 223300, China)(2.Huaian Environmental Monitoring Central Station, Huaian Jiangsu 223001, China)(3.Huaian Water Co., Ltd., Huaian Jiangsu 223002, China)
Rapid and effective monitoring of the formation of disinfection by-products is critically fundamental forsecurity management of drinking water. On the basis of conventional chlorination process in waterworks, influences of four key factors, i.e. value of UV254and ammonia nitrogen of the raw water, chlorine dosage andpHvalue, on the formation of typical disinfection by-products-trihalomethanes (THMs) and haloacetic acids(HAAs), were investigated by employing uniform design method, and linear and nonlinear fitting for the data were achieved by utilization of the statistical software-SPSS. Results showed that linear modelgot a good fitting for THMs and HAAs generation with the correlation coefficients both greater than 0.9, while correlation coefficientsof nonlinear model for the above twowere 0.924 and 0.883,respectively.Moreover, correlation coefficientsof linear model for THMs and HAAs generation were both 1 in consideration of the interactions between various factors,which further improved the fitting effects. The new approach forrapid prediction of THMs and HAAsgeneration proposed in this study, may provides a theoretical basis for improvingrisk management of drinking watersecurity under detection-restricted conditions.
chlorination disinfection; disinfection by-products; mathematical model; uniformdesign
2013-05-15
淮安市科技支撐計劃項目(HAS2010003); 2011年江蘇省大學生實踐創(chuàng)新訓練計劃項目
黃輝(1983-), 男, 江西九江人, 講師, 博士研究生, 研究方向為水污染防治與資源化利用.
X703
A
1671-6876(2013)03-0243-05
[責任編輯蔣海龍]