劉 銘，朱長軍，頓珠加措，尼瑪次仁

（1.河北工程大學(xué)城建學(xué)院，河北邯鄲 057150;2.西藏自治區(qū)水文水資源勘測局日喀則水文水資源分局，西藏日喀則 857000;3.西藏自治區(qū)水利電力規(guī)劃勘測設(shè)計研究院西藏拉薩 850009）

水質(zhì)評價是指按照評價目標，選擇相應(yīng)的水質(zhì)參數(shù)，采用相應(yīng)的水質(zhì)標準和評價方法，對水體的質(zhì)量利用價值及水的處理要求做出評定，使用廣泛的常規(guī)環(huán)境質(zhì)量評價方法大致分為綜合指數(shù)法、數(shù)理統(tǒng)計法、分級評價法和其他方法四大類［1］。在復(fù)雜的環(huán)境系統(tǒng)中存在很多模糊性，例如環(huán)境影響的“大”與“小”、污染物濃度的“高”與“低”、方案的“優(yōu)”與“劣”等均無法用確定的界限加以劃分［2］。在對區(qū)域水環(huán)境質(zhì)量進行水質(zhì)評價時，僅將綜合指數(shù)評價法應(yīng)用于環(huán)境質(zhì)量綜合評價，其理論方法和實際結(jié)果存在著一定的缺陷，往往難以實現(xiàn)對區(qū)域水環(huán)境的總體評價。目前，聚類分析法已被廣泛應(yīng)用于環(huán)境質(zhì)量評價中，聚類分析方法能夠較好評價生物指標，綜合考慮生化指標，有助于合理評價水質(zhì)分布與特性［3］;通過聚類分析能確定海域主要污染指標和污染程度，能夠較合理地反映海域的環(huán)境質(zhì)量特征［4］;許多土壤研究者與環(huán)境工作者運用聚類分析在土壤質(zhì)量評價中開展了卓有成效的工作［5］。大量的研究結(jié)果表明，多指標存在時應(yīng)用聚類分析方法能夠?qū)Νh(huán)境質(zhì)量做出更為客觀、準確的綜合評價。本次水質(zhì)分析采用聚類分析的統(tǒng)計方法，通過對水質(zhì)監(jiān)測監(jiān)測數(shù)據(jù)的統(tǒng)計分析，研究比較各監(jiān)測站點水質(zhì)指標數(shù)據(jù)的相關(guān)關(guān)系，對具有相近數(shù)據(jù)關(guān)系的監(jiān)測站點進行歸類，另外結(jié)合監(jiān)測站點水質(zhì)綜合評價結(jié)果和區(qū)域內(nèi)地理位置及當?shù)厮廴緦嶋H情況，從而對區(qū)域水環(huán)境質(zhì)量做出總體初步評價。

1 監(jiān)測站點資料

濟寧市位于山東省西南部，京杭大運河穿過濟寧市，素有“運河之都”的美譽，我國北方最大的淡水湖南四湖也位于其中。隨著城市的發(fā)展，濟寧市水污染情況越來越嚴重，同時南水北調(diào)工程的啟動也對濟寧市的水質(zhì)提出嚴格要求。濟寧市內(nèi)河流縱橫分布，各段污染狀況不同，水質(zhì)監(jiān)測站點的布點也因此密布整個濟寧市而且數(shù)量較多，這最終導(dǎo)致水質(zhì)評價時會出現(xiàn)繁雜等困難。本次分析選取濟寧市水質(zhì)監(jiān)測數(shù)據(jù)，共選取包括25個地下水監(jiān)測站點的水質(zhì)監(jiān)測結(jié)果，選取C1、SO4、總硬度、BOD5、CODMn、氨氮、亞硝酸鹽氮、硝酸鹽氮、酚共9項水質(zhì)監(jiān)測指標，濟寧市監(jiān)測站點區(qū)域分布見下圖1，各個監(jiān)測站點水質(zhì)監(jiān)測數(shù)據(jù)見表1:

表1 監(jiān)測站點水質(zhì)監(jiān)測數(shù)據(jù)Tab.1 The water quality monitoring data from monitoring sites

2 聚類分析與綜合評價

2.1 監(jiān)測站點的聚類分析

聚類是按照某個特定標準把一個數(shù)據(jù)分割成不同的類或簇，使得類內(nèi)相似性盡可能的大，同時類間的差異性也盡可能的大［6］。聚類分析方法的基本思想是研究對象之間存在著不同程度的相似性，通過分析比較各事物之間的性質(zhì)，將性質(zhì)相近的歸為一類，將性質(zhì)差別較大的歸入不同的類，聚類分析中常見的有系統(tǒng)聚類法和K均值聚類法。本次數(shù)據(jù)分析選用系統(tǒng)聚類法，這是目前實踐中使用最多的，該方法的分析步驟是:先將所有樣本看成一類，選用一定的方法計算樣本與樣本之間的距離和類與類之間的距離;開始，每個樣本自成一類，類類之間和樣本之間的距離相同;然后，在所有的類中選擇距離最小的合并為一個新類，并計算出其與其它類的距離;接著再將距離接近的類合并，直至將所有樣本合并為一類為止。由聚類系譜圖可以清楚的看出全部樣本的聚集過程，進而做出樣本的分類。

由于水質(zhì)各項指標監(jiān)測數(shù)據(jù)存在計量量綱的不同，數(shù)值大小差異較大，因此需要首先對數(shù)據(jù)進行標準化處理，本次選用極差標準化法處理數(shù)據(jù)。假設(shè)有m個聚類對象，每個對象有n個要素構(gòu)成，其所對應(yīng)要素數(shù)據(jù)如表2所示。

表2 聚類對象與要素數(shù)據(jù)Tab.2 The clustering objects and data elements

極差標準化公式如下:

式中Xij為監(jiān)測指標原值，其中1≤i≤m，1≤j≤n為標準化后數(shù)據(jù);為監(jiān)測指標的最小值;為檢測指標的最大值。

經(jīng)過標準化處理所得的數(shù)據(jù)，極大值為1，極小值為0，其余的數(shù)值均在0與1之間。

系統(tǒng)聚類分析的依據(jù)和基礎(chǔ)是數(shù)據(jù)之間的差異性，差異性的測度即為距離的計算，距離較大則樣本相似性越小。距離的計算有多種選擇方法，本次選用Pearson相關(guān)性來計算各樣本之間距離。

系統(tǒng)聚類的過程中類與類之間距離計算方法可以分為單連接法、完全連接法、平均連接法、組平均連接法和離差平方和法等，本次計算類與類之間距離選用組平均連接法，該方法采用的距離定義為兩類之間的平均平方距離，如下:

根據(jù)選取的濟寧市各個監(jiān)測站點得到水質(zhì)監(jiān)測數(shù)據(jù)，利用統(tǒng)計分析軟件SPSS進行聚類分析計算，計算過程如下:

（1）用極差標準化方法對各水質(zhì)指標的原始數(shù)據(jù)進行處理;

（2）采用Pearson相關(guān)性測度監(jiān)測站點之間的樣本間距離;

（3）選用組平均法計算類間的距離，并對樣本進行歸類。

經(jīng)過上述聚類個步驟方法的選擇得到聚類的結(jié)果如圖2示:

SPSS軟件聚類分析所得的樹狀圖能夠清楚顯示出這些監(jiān)測站點的聚類過程，同時根據(jù)聚類分析樹狀圖的結(jié)果，能夠?qū)⒈O(jiān)測站點分為六類，具體監(jiān)測站點分類如下，第一類監(jiān)測站點編號:7、13、1、3、12、4、24、21、8、10、2、6;第二類:5、25、14、15、18、23、17、16;第三類:9;第四類 19;第五類:20;第六類:11、22。

2.2 監(jiān)測站點的綜合評價

采用內(nèi)梅羅綜合指數(shù)評價法，計算各監(jiān)測站點的綜合指數(shù)。首先進行單組分評價，劃分組分所屬質(zhì)量類別，對各類別按表3分別確定單項組分評價分值Fi，然后代入內(nèi)梅羅綜合指數(shù)評價公式計算綜合評價分值F，根據(jù)F值，按照表4劃分地下水質(zhì)量級別。

內(nèi)梅羅綜合指數(shù)評價法公式:

2.3 監(jiān)測站點分類

經(jīng)過聚類分析，濟寧市的監(jiān)測站點一共分為六大類，用綜合指數(shù)法計算各監(jiān)測站點的污染級別，最終將二者數(shù)據(jù)相結(jié)合，各監(jiān)測站點所屬類別及分布區(qū)域和污染級別如下表5所示:

3 結(jié)論

（1）第一類包括12個監(jiān)測站點，24號及21號分布于中山區(qū)外，其余均分布于任城區(qū)，水質(zhì)綜合評價級別除2號為IV類外其余均為II類，對2號監(jiān)測站點水質(zhì)數(shù)據(jù)分析知，其酚類污染物屬國家IV類標準，導(dǎo)致其水質(zhì)級別較差;

（2）第二類包括8個監(jiān)測站點，其中3個分布于任城區(qū)，5個分布于市中區(qū)，水質(zhì)綜合評價級別均分布于III－V類，由監(jiān)測站點數(shù)據(jù)分析知，25號、17號及5號、15號、16號監(jiān)測點均是由于總硬度超出IV、V類標準，因而水質(zhì)級別較差和極差;

（3）9號、19號各為一類，分別位于任城區(qū)和市中區(qū)，雖然綜合評價級別均為良好，但由監(jiān)測數(shù)據(jù)可知其主要污染類別不同，19號其氯化物以及硝酸鹽氮量明顯高于9號;

（4）20號自成一類，位于市中區(qū)，其主要污染物是亞硝酸鹽氮，其含量超出國家V類標準;

（5）11號和22號分別位于任城區(qū)和中山區(qū)，均系酚類超出國家V類標準，因而分為一類，而22號由于硝酸鹽氮屬V類標準，導(dǎo)致水質(zhì)較差。

通過聚類分析法對監(jiān)測站點的分類結(jié)果看出各類間的主要區(qū)別在于水質(zhì)污染級別的不同或是主要污染因素的不同，結(jié)合各個監(jiān)測站點以及污染源分布圖可以得出:市中區(qū)總體水質(zhì)質(zhì)量較差，其主要原因是市區(qū)內(nèi)人口以及工業(yè)企業(yè)密集分布，污染排放量大;任城區(qū)一些高污染排放企業(yè)周圍水質(zhì)較差外，其他地區(qū)水質(zhì)良好，這與任城區(qū)河流分布較多是分不開的，污染物排放后會隨河流流入湖泊，從而減輕了其對地下水污染的影響。

表3 單項組分評價分值Tab.3 Single component evaluation

表4 地下水質(zhì)量劃分級別Tab.4 Groundwater quality grade

表5 綜合評價與聚類結(jié)果Tab.5 Results of comprehensive evaluation and clustering

［1］熊德琪.環(huán)境系統(tǒng)模糊集分析理論與應(yīng)用［M］.大連:大連海事大學(xué)出版社，2001.

［2］馬風(fēng)才.項目環(huán)境影響模糊評價理論與應(yīng)用研究［D］.北京:北京航空航天大學(xué)，2001.

［3］曾淦寧，吳國權(quán)，徐曉群.多元聚類分析方法在杭州灣水質(zhì)分析上的應(yīng)用［J］.浙江工業(yè)大學(xué)學(xué)報，2009（1）:14－19.

［4］羅薇，邵秘華，周立新.聚類分析功能在大連港水域環(huán)境質(zhì)量評價中的應(yīng)用［J］.大連海事大學(xué)學(xué)報，2004（4）:51－55.

［5］蔡瞳，徐惠，吳群.土壤質(zhì)量聚類分析—以封丘縣為例［J］.安徽農(nóng)業(yè)科學(xué)，2008，36（25）:10998－10999.

［6］顧洪博，張繼杯.不確定數(shù)據(jù)的聚類分析研究及應(yīng)用［J］.河北工程大學(xué)學(xué)報:自然科學(xué)版，2012，29（1）:109－112.