某城市排污河水質狀況分析評價

郝 園，孫彥敏，仲曉倩，裴 琨，梁立業

(石家莊市環境監測中心，河北石家莊 050022)

某城市排污河水質狀況分析評價

郝 園，孫彥敏，仲曉倩，裴 琨，梁立業

(石家莊市環境監測中心，河北石家莊 050022)

以城市排污河道為研究對象，測定了表層沉積物中的重金屬總量、水質常規指標及DBP和DEHP，采用改進的主成分分析法和等標有機污染指數分析對排污河污染狀況進行分析評價。運用改進的主成分分析法找出了主要污染源和主要污染物。結果表明：不同排污口的主要污染物不同，B化工廠CODCr及pH值為主要污染物；C泵站污染嚴重的為CODCr，pH值和TP；D排污口主要污染物為CODCr，NH3-N，Cr和pH值；E化工廠的主要污染物為CODCr，NH3-N和Cd。D排污口是第一主成分的主要污染源點，B化工廠無論是DBP等標污染指數還是DEHP等標污染指數均明顯高于其他排污口，是A排污河PAEs污染的主要污染源。

主成分分析法；DBP；DEHP；等標污染指數分析

A排污河數十年來一直是市政污水的主要受納河體。統計、監測資料顯示，該市的市區及郊縣的化工廠、煤氣廠、造紙廠、印染廠、制革廠、染化廠、電池廠等均直接向A排污河排放工業廢水，其運行情況直接影響著工農業生產和人民生活。A排污河為地表水，屬于地表水環境質量標準(GB 3838—2002)中的V類水體。水體環境受到諸多指標的影響，加大了水質評價的難度，客觀、合理評價水體綜合水質尤顯重要[1-2]。

主成分分析法是一種成熟的數據降維或特征提取的方法，屬于數理統計的應用范疇[3]。主成分分析法能夠較為合理地對評價因子進行賦值，在水質評價因子賦值中具有廣闊的應用前景。李經偉等運用改進的主成分分析法評價了白洋淀水質狀況[4]。地表水環境質量評價應根據水域功能類別，選取相應類別標準，進行改進的主成分分析評價。近年來發現酞酸酯可對環境、生物與食品等形成污染[5]。鄰苯二甲酸酯急性毒性不大，對Ames試驗呈陰性反應，但在大劑量情況下，對動物有致畸和致突變的作用，顯示較強的內分泌干擾性[6]。

本研究選取A排污河具有代表性的污染源點，對其進行了多種水質指標以及有機污染物的監測。運用數學統計方法找出主要污染因子和重要污染源點，初步探討了水質的影響因素。其結果不僅可為該市排污河道的管理和治理提供指導，而且對于重點污染源的治理具有重要的參考價值。

1 材料與方法

1.1檢測指標

選取某年5月到10月期間采集的24個水樣進行分析，按照《水和廢水監測分析方法》[7]對水樣進行測定。常規檢測指標：pH值，SS，好氧有機物(CODCr類)，生源要素(TN，TP，NH3-N)，重金屬元素(Cu，Fe，Zn，Mn，Pb，Cd，Ni，Cr)；有機物檢測指標：鄰苯二甲酸二丁酯(DBP)和鄰苯二甲酸異辛酯(DEHP)。

1.2評價方法

目前，水環境質量的綜合評價有很多方法，本研究采用改進的主成分分析法對排污河水質狀況進行評價。采用SPSS軟件進行計算得到各指標的特征值、主成分的貢獻率和累計貢獻率。累計貢獻率說明主成分所包含全部指標信息的比例。前n個主成分的累計方差貢獻率達到85%以上，就可以充分反映原始數據的主要信息。因此，可以利用前n個主成分來對排污河的水質污染狀況進行可比性研究。

設pH值為x1，CODCr為x2，SS為x3，總氮為x4，氨氮為x5，總磷為x6，銅為x7，鋅為x8，鉻為x9，鎘為x10，鉛為x11，鐵為x12，錳為x13，鎳為x14。

主成分分析方法是一種將多維因子納入同一系統中進行定量化研究、理論比較完善的多元統計分析方法。改進的主成分分析法是對傳統的主成分分析法的改進，適用于原始變量之間呈現出非線性關系的情況。它是一種新的評價模型，使得研究有一定的理論及實踐價值[8]。改進的主成分分析可有效地弱化異常值的影響，增強分析結果的穩定性，提高分析的可信度，具有良好的降維效果[9]。

改進的主成分分析法的主要步驟如下：

第一，構造水質數據矩陣X=[x1,x2,…,xp]；

第三，計算矩陣s的特征值和特征向量，并按累計方差貢獻率大于85%來確定一個m值；

第四，寫出前m個主成分計算公式，則第i個主成分為

A=(aij)n×p為單位正交化的特征向量矩陣(e1,e2,e3,…,ep)中的各項，又稱主成分載荷矩陣；

第五，對所選主成分作物理解釋；

第六，求主成分值及綜合主成分值[10]。

2 結果與分析

A排污河屬于地表水，加之其排放污水的性質，按照GB 3838—2002《地表水環境質量標準》[11]中的V類水為主要評價依據，對各個污染源進行水質評價。

2.1改進的主成分分析法評價結果

2.1.1 各斷面間的相對污染程度

據4月到5月對A排污河4個污染源點采樣和分析得到的數據，以污染源點和監測指標為基礎，建立水質數據矩陣。

利用SPSS軟件計算得到協方差矩陣的特征值、特征向量，以及主成分方差貢獻率、累計方差貢獻率。

方差分解主成分提取分析見表1，前m個主成分的載荷矩陣見表2，A排污河污染源點分析結果見表3—表6。

表1 方差分解主成分提取分析

表2 前m個主成分的載荷矩陣

表3 5月主成分得分及污染源排名

表4 7月主成分得分及污染源排名

表5 8月主成分得分及污染源排名

表6 10月主成分得分及污染源排名

改進的主成分分析中，由于同趨勢化后所有指標為正指標，所以只需通過計算綜合得分進行排序，無需考慮正負值關系，即可得到評價結果。首先選取A河5月、7月、8月、10月份監測數據進行主成分分析。由表1可知，5月份第一主成分和第二主成分累計貢獻率達到91.23%，并且10月份累計貢獻率為100%。因此可以用其代表A河重點污染源的綜合污染情況。通過表2的5月份的載荷矩陣可知第一主成分為pH值及CODCr，SS，TP，Cr，Cd，Pb，Fe，Mn；第二主成分為TN，Cu，Zn，Ni。由表1—表3可知，B化工廠在第一主成分的污染上最為嚴重，D排污口在第二主成分的污染最為嚴重。由表3可知，5月份C泵站總排名第2，且其在第一主成分上排名第2，說明5月份C泵站主要污染現象與第一主成分相關性較大。E化工廠在5月份的排名均在第4，說明污染相對較輕。

A河7月份數據的第一主成分和第二主成分累計貢獻率已經超過85%，因此可以用其表征7月份主要污染。由表2可知：第一主成分為Mn，Cu，Ni，Zn，Fe，Cr，pH值，CODCr，TP；第二主成分為SS，NH3-N。可以看出，重金屬污染和pH值指標有很大的關系。同時由表4可知，整個7月份D排污口的第一主成分和綜合排名均為第1，說明其是7月份的主要污染源點，B化工廠為第4污染源點，和5月份相比B化工廠污染暫時減輕，說明污染情況相對于其他點稍微減輕。而D由5月份的綜合第3上升為第1，成為最嚴重的污染源。E化工廠排名呈上升趨勢。8月份前2個主成分累積貢獻率為92.6%。第一主成分包括SS，TN，NH3-N，TP，pH值，Cu，Zn，Cr，Fe，Mn；第二主成分包括Pb，Cd，CODCr。由表5可知，8月份B化工廠污染有反復現象，綜合污染排名仍為D排污口第1。表6中10月份B排污河的數據變化起伏較大。由表1、表2可知：10月份分為3個主成分，其特征值分別為7.88，4.48和1.64，第一主成分為pH值，TP，Cu，Zn，Cr，Fe，Mn；第二主成分為CODCr，SS，Cd，Ni；第三主成分為NH3-N和TN。第一主成分可以歸納為主要重金屬污染，第二主成分主要為無機鹽類污染，第三主成分為營養鹽類污染。

由表1—表6可知，D排污口為重金屬污染和pH值污染相當嚴重的源點，這與該點pH值呈強酸性吻合，并且氮污染也較為嚴重。B化工廠為有機類污染，其中CODCr的污染載荷最高，也與實際監測值吻合。C泵站為營養鹽污染最嚴重的點。從綜合排名可以看出，D排污口和C泵站排污口在后幾個月的污染呈逐漸加重的趨勢，B化工廠有緩和的趨勢。但是通過幾個月的主成分分析可以認為，D排污口是第一主成分的主要污染源點。E化工廠污染的一個曲折性可能與E化工廠的生產狀況相關。

2.1.2 各污染源點污染物的主成分分析

由改進的主成分分析法，分析各個排污口重點污染物。污染物排名見表7。D排污口的第一主成分的貢獻率為99.829%，因此可以用第一主成分分析綜合污染。D排污口排名前3的是CODCr，NH3-N，Cr，反映了該點有機污染嚴重。從整個污染物來看，D排污口是最為主要的污染源且CODCr是其主要污染現象。綜合主成分分析及單因子分析可知，D排污口主要污染物為CODCr，NH3-N，Cr和pH值。C泵站前2個主成分累計貢獻率約為85%。可以用前2個主成分表示污染物的綜合排名，污染嚴重的為CODCr，pH值和TP。而由因子分析法可知，C泵站5月份的主要污染現象為氮污染，這與主成分分析稍有出入，可能是數據量不足的緣故。B化工廠采用前2個主成分累計貢獻率達到87.216%，因此用其表示綜合污染的CODCr，pH值為主要污染物。由于E化廠工缺少一定數據，故不適合進行主成分分析。

表7 排污河污染物綜合排名

2.2排污口有機物污染狀況的等標污染指數分析

中國大部分主要河流水體中近年也多檢出PAEs的存在。調查結果同時顯示水體中PAEs的主要成分為DBP和DEHP[12]。基于水樣全掃描數據、研究報告及環境中普遍存在的有毒有害物質及其對人體的危害，本實驗以鄰苯二甲酸二丁酯(DBP)、鄰苯二甲酸異辛酯(DEHP)作為各排污口有毒有害物質的監測項目。

2.2.1 污染源點污染物分析

圖1和圖2分別為A排污河各排污口PAEs等標污染指數。可以看出A排污河各污染源點在6月、7月、10月份的等標污染指數較高，8月、9月份指數較低。B化工廠無論是DBP等標污染指數還是DEHP等標污染指數均明顯高于其他排污口，是A排污河PAEs污染的主要污染源。D排污口和E化工廠次之。

圖1 A排污河DBP等標污染指數Fig.1 A sewage river DBP equivalent standard pollution index

圖2 A排污河DEHP等標污染指數Fig.2 A sewage river DEHP equivalent standard pollution index

2.2.2 B化工廠主要污染源分析

由于B化工廠為有機污染物的主要污染源，因此對其進行詳細分析。B化工排污口水質全掃描圖譜見圖3和圖4。其中圖3為非極性物質全掃描譜圖，圖4為極性物質全掃描譜圖。

圖3 B化工廠全掃描非極性物質Fig.3 B Chemical factory full scan non-polar substances

圖4 B化工廠全掃描極性物質Fig.4 B Chemical factory full scan polar substances

圖5 B化工廠PAEs等標污染指數圖Fig.5 B PAEs and other standard chemical pollution index chart

圖3中，11.547 min出峰物質為八甲基環四硅氧烷，與C泵站檢出物相同。圖4中，20.259，25.334 min出峰物質均含有羥基與苯環，26.481 min出峰物質為2,4-二叔丁基苯酚。由此可見，B化工排污口中易揮發性有機污染物種類大部分帶苯環，對環境危害較大。B化工排污口水質PAEs分析見圖5。由圖5可見，B化工排污口DEHP污染程度嚴重于DBP。DBP呈現先降低、后升高的趨勢，6月份DBP等污染指數最高，且已超標將近1.4倍，而8月份其指數為0，7月份等標污染指數與10月份類似。DEHP等標污染指數變化規律與DBP類似，6月份最高，超標近2.4倍，7月份和8月份達標排放，等標污染指數將近0.3，9月份含量顯著增多，10月份指數有所下降但仍比較大。在監測時段內，6月份DBP超標率均為20%。

3 結 論

通過改進的主成分分析法和有機污染物的等標污染指數分析，對A排污河水質狀況進行了研究，得到如下結論。

B化工廠CODCr和pH值為主要污染物；C泵站污染嚴重的為CODCr，pH值和TP；D排污口主要污染物為CODCr，NH3-N，Cr和pH值；E化工廠的主要污染物為CODCr，NH3-N和Cd。

從綜合排名可以看出，D排污口和C泵站排污口在后幾個月的污染有逐漸加重的趨勢，B化工廠有緩和的趨勢。但是通過幾個月的主成分分析可以認為，D排污口是第一主成分的主要污染源點。

無論是DBP等標污染指數還是DEHP等標污染指數，B化工廠均明顯高于其他排污口，是A排污河PAEs污染的主要污染源。

/

[1] 汪天祥，許士國，韓成偉. 改進主成分分析法在南淝河水質評價中的應用[J].水電能源科學,2012,30(10)：33-36.

WANG Tianxiang, XU Shiguo,HAN Chengwei.Application of improved principal component analysis method in Nansi river water quality evaluation[J]. Water Resources and Power,2012,30(10)：33-36.

[2] 楊會龍，張志軍，張雅楠.污水廠改良氧化溝處理工藝的優化控制[J].河北科技大學學報，2009，26(3)：257-261.

YANG Huilong,ZHANG Zhijun,ZHANG Ya′nan.Optimized control of improved oxidation ditch craft for sewage treatment plant[J].Journal of Hebei University of Science and Technology,2009,26(3):257-261.

[3] 王 維，紀 枚，蘇亞楠.水質評價研究進展及水質評價方法綜述[J]. 科技情報開發與經濟，2012，22(30)：129-136.

WANG Wei, JI Mei，SU Ya′nan.Reviews of the progress in the research of water quality evaluation and the methods for water quality evaluation[J]. Sci-Tech Information Development & Economy, 2012,22(30):129-136.

[4] 李經偉，楊路華，梁寶成,等.改進的主成分分析法在白洋淀水質評價中的應用[J].海河水利，2007(3)：40-43.

LI Jingwei, YANG Luhua, LIANG Baocheng et al. The application of improved principal component analysis method to the water quality appraisement in Baiyangdian Lake[J].Haihe River Water Resources,2007(3):40-43.

[5] 趙文紅.酞酸酯類增塑劑毒理研究進展[J].環境與職業醫學,2003,20(2):135-138.

ZHAO Wenhong.Advances in toxicological studies on plasticizers phthalate esters[J].Journal of Labour Medicine, 2003,20(2):135-138.

[6] 胡雄星,韓中豪,劉必寅,等.鄰苯二甲酸酯的毒性及其在環境中的分布[J].環境科學與管理,2007,32(1):37-40.

HU Xiongxing, HAN Zhonghao, LIU Biyin, et al. Distribution on phthalic acid esters in environment and its toxicity[J].Environmental Science and Management,2007,32(1):37-40.

[7] 國家環境保護總局.水和廢水監測分析方法[M].北京：中國環境科學出版社，2002.

State Environmental Protection Administration.Water and Wastewater Monitoring Analysis Method[M]. Beijing: China Environmental Science Press,2002.

[8] 馬志明.基于改進主成分分析法的河流水質評價[J].才智，2012(19)：298-299.

MA Zhiming. Evaluation of river water quality improvement based on principal component analysis[J].Intelligence，2012(19)：298-299.

[9] 高 艷,于飛西.一種用于綜合評價的主成分分析改進方法[J].西安文理學院學報(自然科學版)，2011,14(1)：105-107

GAO Yan, YU Feixi. A modified principal component analysis algorithm for comprehensive evaluation[J].Journal of Xi′an University of Arts & Science(Natural Science Edition)，2011,14(1)：105-107.

[10] 姚煥玫，黃仁濤，甘復興，等.用改進的主成分分析法對東湖的水質污染進行評價[J].武漢大學學報(信息科學版)，2005，30(8)：732-735.

YAO Huanmei, HUANG Rentao, GAN Fuxing, et al. Principal component analysis of the water quality evaluation in East Lake[J].Geomatics and Information Science of Wuhan University, 2005,30(8):732-735.

[11] GB 3838—2002,地表水環境質量標準[S].

GB 3838—2002, Environmental Quality Standard for Surface Water[S].

[12] SHA Yujuan,XIA Xinghui,YANG Zhifeng,et al.Distribution of PAEs in the middle and lower eaches ofthe Yellow River of China[J].Environ Monit Assess,2007,124:277-287.

Water quality assessment of some city drainage river

HAO Yuan，SUN Yanmin，ZHONG Xiaoqian，PEI Kun，LIANG Liye

(Shijiazhuang Municipal Environmental Monitoring Center, Shijiazhuang Hebei 050022, China)

The A drainage river in Tianjin is studied. The pollutants in the surface layer of the river sediment is measured, including heavy metals, water quality indexes, DBP and DEHP, by improved principal component analysis method and organic pollution index method. The experimental results show that:(1) The main pollutants of each outfall is different, for example, the B chemical plant' main pollutants are CODCrand pH; the C pump station's main pollutants are TP, CODCrand pH; the D outlet's main pollutants are CODcr, NH3-N, Cr and pH; the E chemical plant's main pollutants are CODCr, NH3-N and Cd. (2)D outfall is the main source of the first component.(3) the DBP and DEHP of B chemical plant are significantly higher than that of other outfalls, which are the main source of PAEs of A drainage river.

principal component analysis method; DBP; DEHP; equal-standard organic pollutant index analysis

1008-1534(2013)05-0391-07

X824

A

10.7535/hbgykj.2013yx0519

2013-04-30;

2013-06-06

責任編輯：陳書欣

郝 園(1986-)，男，河北正定人，碩士，主要從事環境監測方面的研究。

E-mail:523583582@qq.com