萬(wàn)　蕾，張曼玉，陸　晟

(徐州工程學(xué)院，江蘇 徐州　221018)

?

基于多元統(tǒng)計(jì)分析的奎河(徐州段)水質(zhì)污染現(xiàn)狀

萬(wàn)蕾，張曼玉，陸晟

(徐州工程學(xué)院，江蘇 徐州221018)

摘要：奎河是徐州境內(nèi)一條重要的河流.因近年生活污水、工業(yè)廢水的隨意排放，奎河水體的污染問(wèn)題愈發(fā)嚴(yán)重.文章分析了奎河(徐州段)7個(gè)監(jiān)測(cè)斷面的水質(zhì)基本情況，并采用聚類分析和主成分分析對(duì)監(jiān)測(cè)斷面和水質(zhì)指標(biāo)進(jìn)行了統(tǒng)計(jì)分析.結(jié)果表明：1)七個(gè)監(jiān)測(cè)斷面的濁度值的變化范圍1.05到1.99 NTU，pH值的變化范圍為7.23~7.83，溶解氧的變化范圍在0.48~2.98 mg/L；2)DO、TP和硝酸鹽氮含量較高，7個(gè)斷面總體屬于V類~劣V類水；3)聚類分析表明7個(gè)斷面可以分為城區(qū)-城鄉(xiāng)結(jié)合部以及農(nóng)村段兩大類；4)主成分分析表明，奎河富營(yíng)養(yǎng)化問(wèn)題嚴(yán)重，磷是首要污染因子，有機(jī)物和氮污染也較為嚴(yán)重.

關(guān)鍵詞：奎河；水質(zhì)；污染；多元分析

奎河開挖于明朝，是一條人工開鑿的泄洪河道，距今已有400多年的歷史.奎河全長(zhǎng)76 km，徐州市境內(nèi)長(zhǎng)25.75 km，流域總面積達(dá)1231 km2.它發(fā)源于云龍湖，流過(guò)奎山，貫穿銅山縣和三堡鎮(zhèn)，在黃橋進(jìn)入安徽，然后匯入淮河.奎河原本是一條十分清澈的小河，自從云龍湖被辟為風(fēng)景游覽區(qū)，便切斷了云龍湖與下游河流的聯(lián)系，奎河便成了無(wú)源之水.徐州市區(qū)黃河故道以北的許多工廠，像造紙、污水處理廠等污染較重的企業(yè)，徐州有關(guān)方面鋪設(shè)起了地下管道，從故黃河的下面穿過(guò)，調(diào)頭朝南，也全都被引入奎河.這樣，奎河就已經(jīng)不是實(shí)際意義上的河流，而成了徐州市地地道道的一條總的下水道.近年來(lái)，對(duì)奎河的整頓一直不斷，但因仍有污水的不斷排入，效果微弱，總體達(dá)不到地表水Ⅴ類水質(zhì)要求[1-2].本文采用多元統(tǒng)計(jì)方法對(duì)奎河(徐州段)的水質(zhì)狀況進(jìn)行分析，以期掌握奎河的污染現(xiàn)狀和污染規(guī)律，為奎河水質(zhì)改善和環(huán)境保護(hù)措施的制定提供參考.

關(guān)于多元統(tǒng)計(jì)的方法用于水質(zhì)分析，前人已經(jīng)有所研究和應(yīng)用.例如，富天乙等[3]通過(guò)層次聚類分析將遼陽(yáng)太子河2009年至2012年3個(gè)斷面的144個(gè)樣本點(diǎn)分為11組，以各組的樣本均值為基礎(chǔ)，采用水質(zhì)標(biāo)識(shí)指數(shù)方法得到各組的水質(zhì)評(píng)價(jià)結(jié)果.范良千等[4]以主成分分析探討主要污染來(lái)源；以聚類分析劃分監(jiān)測(cè)點(diǎn)類別并識(shí)別其空間相似性，結(jié)果顯示梅梁灣水質(zhì)主要受農(nóng)業(yè)非點(diǎn)源、外源輸入的有機(jī)懸浮物等5方面影響，9個(gè)監(jiān)測(cè)點(diǎn)位劃歸為4類.劉海平等[5]根據(jù)Pearson相關(guān)性矩陣，將尼洋河26項(xiàng)理化指標(biāo)分為3大類；析因分析顯示，可將26項(xiàng)理化指標(biāo)分為4 大類；主成分分析發(fā)現(xiàn)，平水期、豐水期以及枯水期水環(huán)境特征明顯不同，以上3期各自聚為一類.多元統(tǒng)計(jì)分析可以較好的識(shí)別水質(zhì)數(shù)據(jù)的時(shí)空特征，識(shí)別監(jiān)測(cè)樣本的內(nèi)在差異性和相似性，結(jié)果客觀可信.

1材料與方法

1.1監(jiān)測(cè)方案

為了全面了解奎河的水質(zhì)污染狀況，2015年4月設(shè)定了奎河的源頭即云龍湖、積水壩、十里閘、黃河路橋、珠江路(奎河銅山段)、楊山頭閘及黃橋斷面7個(gè)采樣點(diǎn)，每個(gè)采樣點(diǎn)采3組平行樣，取其平均值進(jìn)行分析.監(jiān)測(cè)了水的pH值、溶解氧、濁度、有機(jī)物、總磷、硝酸鹽氮6個(gè)指標(biāo).

1.2水質(zhì)指標(biāo)分析方法

水質(zhì)指標(biāo)的測(cè)定按照《水和廢水監(jiān)測(cè)分析方法(第四版)》進(jìn)行[6].pH值、濁度和DO的測(cè)定采用便攜式儀器現(xiàn)場(chǎng)測(cè)定.

水中總磷(TP)的測(cè)定采用鉬銻抗分光光度法，取25 mL樣品于比色管中，加4 mL過(guò)硫酸鉀溶液，密封，高壓消解30 min，加入抗壞血酸和鉬酸鹽溶液，于700 nm波長(zhǎng)下進(jìn)行比色.

水中有機(jī)質(zhì)(COD)含量通過(guò)重鉻酸鉀氧化滴定法測(cè)定，取20.00 mL混合均勻的水樣置于250 mL磨口的回流錐形瓶中，加入0.4 g硫酸汞，搖動(dòng)使溶解，準(zhǔn)確加入10.00 mL重鉻酸鉀標(biāo)準(zhǔn)溶液及數(shù)粒小玻璃珠或沸石，連接磨口回流冷凝管，從冷凝管上口慢慢加入30 mL硫酸-硫酸銀溶液，加熱回流2 h.溶液冷卻后，滴入指示劑，用硫酸亞鐵銨標(biāo)準(zhǔn)溶液滴定，記錄用量.

硝酸鹽氮的測(cè)定方法為酚二磺酸分光光度法，取50 mL水樣于蒸發(fā)皿中，蒸干后加入酚二磺酸試劑，充分研磨，加入氨水，轉(zhuǎn)移到比色管中，于410 nm下進(jìn)行比色.

1.3統(tǒng)計(jì)方法

1.3.1聚類分析

聚類分析(cluster analysis)是一組將研究對(duì)象分為相對(duì)同質(zhì)的群組(clusters)的統(tǒng)計(jì)分析技術(shù).它依據(jù)所研究的樣品或變量之間存在不同程度的相似性，以一些能夠度量樣品或指標(biāo)之間相似程度的統(tǒng)計(jì)量作為劃分類型的依據(jù)，把一些相似程度較高的樣品(或指標(biāo))聚合為一類，其余一些彼此之間相似程度較大的樣品(或指標(biāo))又聚合為另一類，直到把所有的樣品(或指標(biāo))聚合完畢.通常根據(jù)研究對(duì)象的不同，將其分為Q型聚類分析和R型聚類分析，其中Q型聚類分析是對(duì)樣本進(jìn)行分類處理，R型聚類分析則對(duì)變量進(jìn)行分類處理.本文利用SPSS(19.0)統(tǒng)計(jì)分析軟件，選用離差平方和算法和歐式距離，對(duì)研究區(qū)共7個(gè)監(jiān)測(cè)點(diǎn)位進(jìn)行Q型聚類分析.聚類分析結(jié)果最終用聚類樹狀圖直觀地表現(xiàn)出來(lái)[7].

1.3.2主成分分析

主成分分析(principal component analysis，PCA)，將多個(gè)變量通過(guò)線性變換以選出較少個(gè)數(shù)重要變量的一種多元統(tǒng)計(jì)分析方法.近年來(lái)，主成分分析被廣泛的應(yīng)用在分析水環(huán)境污染特征、水環(huán)境管理等方面.主成分分析計(jì)算步驟：1)對(duì)監(jiān)測(cè)數(shù)據(jù)矩陣X進(jìn)行數(shù)據(jù)標(biāo)準(zhǔn)化處理；2)計(jì)算出標(biāo)準(zhǔn)化數(shù)據(jù)的相關(guān)系數(shù)矩陣，并求出其特征根；3)根據(jù)累積方差貢獻(xiàn)率或特征值大于1原則確定主成分個(gè)數(shù)；4)對(duì)主成分進(jìn)行分析解釋[8].

2結(jié)果與討論

2.1奎河水質(zhì)污染基本情況

奎河水質(zhì)污染基本情況如表1所示.

表1　奎河監(jiān)測(cè)斷面水質(zhì)基本情況

從表1中可以看出，7個(gè)監(jiān)測(cè)斷面的濁度值的變化范圍1.05~1.99 NTU，積水壩的濁度值最小，為1.05 NTU.黃河橋路斷面的濁度值也較小，為1.11 NTU，黃河橋監(jiān)測(cè)斷面處的奎河河道寬，且水流流速要大于其他監(jiān)測(cè)點(diǎn).其他5個(gè)監(jiān)測(cè)斷面的濁度值都處于較高水平.由于條件的限制，云龍湖斷面的采樣點(diǎn)位于湖岸不遠(yuǎn)處，這可能會(huì)與云龍湖水質(zhì)的真實(shí)情況有所出入.7個(gè)監(jiān)測(cè)斷面的pH值的變化范圍為7.23~7.83，云龍湖斷面的pH值最大，為7.83，水呈弱堿性水，云龍湖斷面采樣點(diǎn)為第一采樣點(diǎn)，采樣時(shí)間在上午9點(diǎn)左右.pH值受很多因素的影響，溫度、時(shí)間、水中溶解的CO2、水中的酸性或堿性的溶解性物質(zhì)等，如溫度升高，水中的溶解度會(huì)下降，CO2會(huì)以其體型是釋放到大氣中，這樣會(huì)導(dǎo)致水體的pH值明顯的增大.黃河橋的pH值最小，其他五個(gè)監(jiān)測(cè)點(diǎn)的pH值處在差不多的水平上.水體溶解氧的變化范圍在0.48~2.98 mg/L，十里閘監(jiān)測(cè)斷面的DO最低，為0.48 mg/L，十里閘監(jiān)測(cè)斷面附近有一定量的生態(tài)浮床，種植大量的植物，主要有千屈菜和黃菖蒲.監(jiān)測(cè)時(shí)間為上午九點(diǎn)半作用，植物才開始新一天的光合作用，夜間植物沒(méi)有光合作用，僅有呼吸作用，加上水生動(dòng)物、微生物等的呼吸作用，水總?cè)芙庋鯐?huì)大幅度下降，這就是十里閘DO值較低的的原因，DO最高的是云龍湖，云龍湖水面開闊，與空氣接觸面大，所以DO值較高，又因?yàn)闆](méi)有大面積水生植物生長(zhǎng)，所以DO值最大.其他五個(gè)監(jiān)測(cè)斷面的溶解氧值在相似水平上.

從表1還可以看出，7個(gè)水樣TP的變化范圍在0.27~0.43 mg/L，參考地表水水水質(zhì)標(biāo)準(zhǔn)可知，Ⅳ類水的限值為0.3 mg/L，所有監(jiān)測(cè)斷面僅黃河橋和珠江路為Ⅳ類水，云龍湖和黃橋監(jiān)測(cè)斷面超過(guò)V類水限值的0.4 mg/L，為劣V類水，其他四個(gè)監(jiān)測(cè)斷面水質(zhì)均為V類水.7個(gè)監(jiān)測(cè)斷面的COD值中，黃橋監(jiān)測(cè)面的COD值最高，為37.11 mg/L超過(guò)了地表水Ⅳ類水，為V類水.最低的是十里閘監(jiān)測(cè)斷面18.95 mg/L，屬Ⅲ類水，這個(gè)值要低于其他監(jiān)測(cè)斷面的COD值，可能與十里閘附近水面種植了大量的浮床植物有關(guān).其他監(jiān)測(cè)斷面的COD值在相似水平，為Ⅲ類水.黃橋監(jiān)測(cè)斷面的水樣硝酸鹽氮值處于最高水平，為1.17 mg/L，黃橋的污染較重，農(nóng)村的衛(wèi)生條件較差，生活垃圾多，生活污水大量排放，也是導(dǎo)致硝酸鹽氮值過(guò)高的原因.黃河橋路的硝酸鹽氮的含量也較高.積水壩的硝酸鹽氮處于最低水平，為0.78 mg/L.其它斷面的硝酸鹽氮的含量水平較為接近.

徐州的云龍湖水質(zhì)相對(duì)較好，云龍湖在2005年4月—5月期間，TP的質(zhì)量濃度為0.035~0.061mg/L，DO濃度為6.86~8.17 mg/L，高錳酸鹽指數(shù)6.8 mg/L左右[9].2009年8月到2010年3月期間，TP的質(zhì)量濃度為0.088~0.136 mg/L，DO質(zhì)量濃度為6.67~7.55 mg/L，高錳酸鹽指數(shù)質(zhì)量濃度為4.54~6.27 mg/L[10].可見(jiàn)，奎河水體DO含量低，TP和有機(jī)物污染嚴(yán)重.

2.2奎河水質(zhì)的多元統(tǒng)計(jì)分析

2.2.1聚類分析

奎河7個(gè)點(diǎn)位聚類分析見(jiàn)圖1.從圖中可以看出，7個(gè)點(diǎn)位可以分為4組.其中1(云龍湖)和3(十里閘)可以分為一組，2(積水壩)和4(黃河橋路)為一組，5(珠江路橋)和7(黃橋斷面)為一組，6(楊山頭閘)單獨(dú)分為一組.但總體來(lái)說(shuō)，可以分為兩大類，1(云龍湖)、2(積水壩)、3(十里閘)、4(黃河橋路)監(jiān)測(cè)點(diǎn)為第一大類，5(珠江路橋)、6(楊山頭閘)、7(黃橋斷面)監(jiān)測(cè)點(diǎn)為第二大類.第一大類的4個(gè)斷面位于城區(qū)—城鄉(xiāng)結(jié)合部，第二大類的3個(gè)斷面位于農(nóng)村.說(shuō)明從珠江路以南河段主要受農(nóng)業(yè)非點(diǎn)源污染的影響，而從云龍湖出口到黃河橋路河段受到點(diǎn)源和城市非點(diǎn)源及農(nóng)業(yè)非點(diǎn)源污染的共同影響.

圖1　監(jiān)測(cè)點(diǎn)聚類分析圖

2.2.2主成分分析

奎河7個(gè)監(jiān)測(cè)點(diǎn)的水質(zhì)指標(biāo)的相關(guān)系數(shù)見(jiàn)表2.從表2中可以看出，pH值和DO及TP的相關(guān)系數(shù)較高，分別為0.572和0.602.濁度和TP的相關(guān)系數(shù)也較高，為0.544，說(shuō)明降雨導(dǎo)致奎河濁度升高，同時(shí)TP濃度增加，非點(diǎn)源污染也是奎河流域的一大污染源.奎河水質(zhì)主成分矩陣見(jiàn)表3.從表3中可以看出，前3個(gè)主成分累計(jì)反映了原始數(shù)據(jù)總信息量的81.03%，起到了既能精簡(jiǎn)數(shù)據(jù)又能較完整地保持信息量的作用，故抽取前3個(gè)因子作為分析對(duì)象.從因子載荷大小看，第一主成分方差貢獻(xiàn)率為37.58%.與第一主成分相關(guān)系數(shù)較高的水質(zhì)參數(shù)是TP和pH值.其中TP的因子荷載最高為0.864，說(shuō)明奎河富營(yíng)養(yǎng)化問(wèn)題嚴(yán)重，磷是首要污染因子.第二主成分因子載荷占23.92%，與第二主成分相關(guān)系數(shù)較高的是COD、硝態(tài)氮和DO，可見(jiàn)水體有機(jī)質(zhì)含量較高，有機(jī)物和氮污染也較為嚴(yán)重.第三主成分因子載荷占19.53%，與第三主成分相關(guān)系數(shù)較高的主要是濁度.

表2　水質(zhì)指標(biāo)相關(guān)系數(shù)矩陣

表3　奎河水質(zhì)主成分矩陣

3結(jié)論

通過(guò)對(duì)奎河(徐州段)7個(gè)監(jiān)測(cè)斷面的水質(zhì)監(jiān)測(cè)，可以得出以下結(jié)論：

1)7個(gè)監(jiān)測(cè)斷面的濁度值的變化范圍1.05~1.99 NTU，pH值的變化范圍為7.23~7.83，溶解氧的變化范圍在0.48~2.98 mg/L；COD、TP和硝酸鹽氮含量較高，7個(gè)斷面總體屬于V類~劣V類水，水體DO含量低，TP和有機(jī)物污染嚴(yán)重.

2)聚類分析表明7個(gè)斷面可以分為城區(qū)—城鄉(xiāng)結(jié)合部以及農(nóng)村段兩大類，從珠江路以南河段主要受農(nóng)業(yè)非點(diǎn)源污染的影響，而從云龍湖出口到黃河橋路河段受到點(diǎn)源和城市非點(diǎn)源及農(nóng)業(yè)非點(diǎn)源污染的共同影響.

3)主成分分析表明，奎河富營(yíng)養(yǎng)化問(wèn)題嚴(yán)重，磷是首要污染因子，有機(jī)物和氮污染也較為嚴(yán)重.

參考文獻(xiàn)：

[1] 魏燦,解恒參,朱亦仁.奎河水污染現(xiàn)狀及對(duì)策[J].徐州師范大學(xué)學(xué)報(bào)(自然科學(xué)版)，2005,23(2):63-67.

[2] 李炬.略談奎河的水環(huán)境保護(hù)與水資源利用[J].江蘇環(huán)境科技,2001,14(3):28-30.

[3] 富天乙,鄒志紅,王曉靜.基于多元統(tǒng)計(jì)和水質(zhì)標(biāo)識(shí)指數(shù)的遼陽(yáng)太子河水質(zhì)評(píng)價(jià)研究[J].環(huán)境科學(xué)學(xué)報(bào),2014,34(2):473-480.

[4] 范良千,吳祖成,張清宇,等.多元統(tǒng)計(jì)方法用于太湖梅梁灣水質(zhì)特征識(shí)別[J].浙江大學(xué)學(xué)報(bào)(理學(xué)版),2013,40(3):308-313.

[5] 劉海平,鐘國(guó)輝,葉少文,等.西藏尼洋河水環(huán)境特征多元統(tǒng)計(jì)分析[J].湖泊科學(xué),2015，27(6):1187-1196.

[6] 國(guó)家環(huán)保局.水和廢水監(jiān)測(cè)分析方法[M].4版,北京：中國(guó)環(huán)境科學(xué)出版社,2002.

[7] 沈楊,何江濤,王俊杰,等.基于多元統(tǒng)計(jì)方法的地下水水化學(xué)特征分析：以沈陽(yáng)市李官堡傍河水源地為例[J].現(xiàn)代地質(zhì),2013,27(2):440-447.

[8] 范良千,吳祖成,張清宇,等.多元統(tǒng)計(jì)方法用于太湖梅梁灣水質(zhì)特征識(shí)別[J].浙江大學(xué)學(xué)報(bào)(理學(xué)版),2013,40(3):308-313.

[9] 蔣滔,李多松,張凱.徐州市云龍湖水質(zhì)的模糊數(shù)學(xué)綜合評(píng)價(jià)[J].環(huán)境研究與監(jiān)測(cè),2007,20(3): 16-19.

[10] 呂順,余莉琳,周猛,等.徐州市云龍湖水質(zhì)監(jiān)測(cè)與評(píng)價(jià)[J].徐州工程學(xué)院學(xué)報(bào)(自然科學(xué)版),2010,25(4):42-46.

(編輯崔思榮)

The Current Situation of Water Pollution of Kui River(Xuzhou Section) Based on Multivariate Statistical Analysis

WAN Lei,ZHANG Manyu,LU Sheng

(College of Environmental Engineering Xuzhou Institute of Technology, Xuzhou 221018, China)

Abstract:Kui River is an important river in Xuzhou. In recent years, the pollution of the river is more serious because of the free discharge of sewage and industrial wastewater.In this paper,the water quality of the 7 monitoring sections of the Kui River (Xuzhou section) was analyzed,and the statistical analysis was carried out using the cluster analysis and principal component analysis.The results showed that:(1)the turbidity value of seven monitoring sections changed from1.05 to 1.99NTU,pH 7.23~7.83, dissolved oxygen 0.48~2.98mg/L;(2)higher content of COD,TP and nitrate nitrogen,seven section overall belonged to class V~inferior V class of water;(3)clustering analysis showed that seven section could be divided for urban and rural combination section and rural section;(4)principal component analysis showed that eutrophication problem of Kui River was serious and the phosphorus was the main pollution factor, organic matter and nitrogen pollution was also serious.

Key words:Kui River; water quality; pollution; multivariate statistical analysis

中圖分類號(hào)：X171

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1674-358X(2016)01-0053-05

作者簡(jiǎn)介：萬(wàn)蕾(1981-)，女，副教授，博士，主要從事水污染防治研究.

基金項(xiàng)目：住房和城鄉(xiāng)建設(shè)部科學(xué)技術(shù)項(xiàng)目(2013-K7-14)

收稿日期：2016-01-03