曹 毅

(宿遷市環(huán)境監(jiān)測中心站，江蘇 宿遷 223800)

·解析評價·

物元分析法在水質(zhì)監(jiān)測優(yōu)化布點中的應(yīng)用

曹 毅

(宿遷市環(huán)境監(jiān)測中心站，江蘇 宿遷 223800)

物元分析法是處理不相容問題的一種有效方法。以竺山灣2010年7月水質(zhì)監(jiān)測數(shù)據(jù)為例，利用物元分析法分析太湖竺山灣水質(zhì)監(jiān)測點位分布情況，并將分析結(jié)果與系統(tǒng)聚類法的分析結(jié)果相比較，發(fā)現(xiàn)物元分析法的分析結(jié)果是準(zhǔn)確、可信的。綜合使用物元分析法和系統(tǒng)聚類法，最終確定水質(zhì)監(jiān)測點位可由14個優(yōu)化至11個。

水質(zhì)監(jiān)測;物元分析法;優(yōu)化布點

環(huán)境監(jiān)測布點通常涉及多項污染指標(biāo)，且各單項指標(biāo)優(yōu)選出的點位往往是不相容的。物元分析法是處理不相容問題的一種有效方法［1，2］，可用于解決環(huán)境監(jiān)測中多指標(biāo)優(yōu)化布點的問題［3-7］。合理設(shè)置的水環(huán)境監(jiān)測點位不但能夠節(jié)約大量的人力、物力、財力，而且能夠使監(jiān)測數(shù)據(jù)更為準(zhǔn)確、有效。

1 材料與方法

1.1 樣品采集及處理

竺山灣位于太湖常州水域(120°00′E—120°05′E，31°20′E—31°28′N)，面積約 88 km2，約占整個太湖總面積的4%，平均水深約2 m，多條入湖河道由此入湖，是太湖暴發(fā)藍藻水華頻率最大的湖灣之一［8］。竺山灣水環(huán)境質(zhì)量狀況直接影響整個太湖的水環(huán)境質(zhì)量。

7月是太湖藍藻水華最嚴(yán)重的時期，選擇7月水質(zhì)監(jiān)測數(shù)據(jù)來優(yōu)化竺山灣水質(zhì)監(jiān)測點位，共設(shè)置監(jiān)測點位14個，按照其相對位置，可分為5個區(qū)域(表1)。每周一現(xiàn)場采樣并帶回實驗室分析，現(xiàn)場操作嚴(yán)格按照《水質(zhì)采樣技術(shù)指導(dǎo)》完成［9］。

1.2 分析項目和方法

水樣的高錳酸鹽指數(shù)(CODMn)采用酸性法測定，總磷(TP)采用鉬酸銨分光光度法，總氮(TN)采用堿性過硫酸鉀消解紫外分光光度法，氨氮(NH3-N)采用納氏試劑比色法，14個點位共取得4批次224個水質(zhì)監(jiān)測數(shù)據(jù)，各點位監(jiān)測數(shù)據(jù)月平均值如表1所示。

表1 2010年7月竺山灣水質(zhì)監(jiān)測平均值統(tǒng)計結(jié)果mg/L

續(xù)表

1.3 數(shù)據(jù)統(tǒng)計

1.3.1 物元分析法

通過建立各點位物元分析模型和節(jié)域分析模型，計算各點位的線性關(guān)聯(lián)函數(shù)值和綜合關(guān)聯(lián)函數(shù)值，根據(jù)計算結(jié)果繪制點聚圖，根據(jù)圖中點位分布和實際情況，確定優(yōu)化的點位［2］。

1.3.1.1 建立各點位的物元分析模型

比較全部監(jiān)測點位的污染指標(biāo)，選出每項指標(biāo)的最佳值A(chǔ)、最劣值B及期望值C，4項指標(biāo)的最大值(最劣值)、最小值(最佳值)及平均值(期望值)，用每項指標(biāo)的最佳值A(chǔ)、最劣值B及期望值C構(gòu)造2個標(biāo)準(zhǔn)物元矩陣:

由2個標(biāo)準(zhǔn)物元矩陣構(gòu)成節(jié)域物元矩陣;每個測點作為一個物元，由全部污染指標(biāo)測值構(gòu)成各點位物元矩陣:

式中:a1，a2，…am——各指標(biāo)的最佳值;b1，b2，…bm——各指標(biāo)的最劣值;c1，c2，…cm——各指標(biāo)的期望值;Q1，Q2…Qm——各項指標(biāo);Xi1，Xi2…Xim——i測點 m 指標(biāo)的監(jiān)測值;i=1，2，3，…n;j=1，2，3，…m。

1.3.1.2 計算各點位綜合關(guān)聯(lián)函數(shù)

將表1中數(shù)據(jù)代入式(5)、(6)，分別計算各點位污染指標(biāo)的對標(biāo)準(zhǔn)物元A、B的線性關(guān)聯(lián)函數(shù):

式中:Xij——第 i個測點第 j項指標(biāo)測值;Cj——第 j項指標(biāo)期望值;Aj——第 j項指標(biāo)最小值(最佳值);Bj——第j項指標(biāo)最大值(最劣值)。

關(guān)聯(lián)函數(shù)值的大小表示測點符合標(biāo)準(zhǔn)物元的程度，意義見表2。

表2 關(guān)聯(lián)函數(shù)的意義

由以上結(jié)果計算出全部污染指標(biāo)對A、B的綜合關(guān)聯(lián)函數(shù):

式中:Wj——j污染指標(biāo)的權(quán)值。

采用地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)的指數(shù)超標(biāo)法計算［10］，得到各項污染指標(biāo)歸一化權(quán)值如表3所示。

表3 地表水污染指數(shù)分級標(biāo)準(zhǔn)及歸一化權(quán)值

其中，Sj為5級標(biāo)準(zhǔn)的均值，Xj/Sj為指數(shù)超標(biāo)法計算的權(quán)值，Xj為第j種污染物所有測點的均值，Wj為第j種污染物的歸一化權(quán)值。

將Wj值代入(7)、(8)式，計算出各采樣點的綜合關(guān)聯(lián)函數(shù)的值見表4。條件最差的點。

表4 綜合關(guān)聯(lián)函數(shù)的值

圖1 14個點的Ka和Kb

1.3.2 系統(tǒng)聚類法

利用SPSS 17.0系統(tǒng)聚類法對表1中數(shù)據(jù)進行聚類分析，利用組間連接的方法，以歐式距離為判斷點位間關(guān)系的依據(jù)［11］，從距離上看，可將點位分為9類，見圖2。

圖2 聚類分析樹狀圖

1.3.1.3 繪制點聚圖

根據(jù)表中的數(shù)據(jù)，結(jié)合關(guān)聯(lián)函數(shù)的含義，以Ka(Xi)為橫坐標(biāo)，Kb(Xi)為縱坐標(biāo)，畫出點聚圖(圖1)。

根據(jù)綜合關(guān)聯(lián)函數(shù)的含義，可以發(fā)現(xiàn)14個點位都分布在第二、第四兩個象限內(nèi)。在第二象限中的點Ka＞0，符合最佳點條件，其中，以11號點符合程度最好;靠近坐標(biāo)中心的點在一定條件下可以向第三、第四象限轉(zhuǎn)化，當(dāng)條件變化時，是最易發(fā)生變化的點;在第四象限中的點Kb＞0，符合最劣點條件，以4號點和2號點最為符合，是點位中相對

2 結(jié)果與討論

物元分析法和聚類分析法結(jié)果比較見表5。

表5 物元分析法和聚類分析法結(jié)果比較

物元分析法將各點位進行優(yōu)選至10個，與聚類分析法優(yōu)選至9個點位進行比較，見表5。從物元分析法與聚類分析法優(yōu)選出的結(jié)果可以看出:第7、8、13、14 號點是一致的獨立點，第 2、4 號點，第3、9號點，第11、12號點，是3組可以優(yōu)化合并的點，這3組點分別在地理上也是鄰近的點位，其他點位2種分類法有著不同的區(qū)分，如第10、14號點，第1、6號點，這2組點在地理上是隔開的，因為水質(zhì)監(jiān)測的結(jié)果相似而歸化為一類，不能作為點位優(yōu)化的依據(jù)。

綜合考慮物元分析法和聚類分析法的優(yōu)化結(jié)果，得到最優(yōu)化的水質(zhì)監(jiān)測斷面。將整個竺山灣點位優(yōu)化為11個點位，即雅浦港、竺山湖中、田雞山、大樁頭、竺山湖南、舊瀆東、田雞山西、沙塘港、殷村港、百瀆港、百瀆港南，既可以反映湖區(qū)的水環(huán)境質(zhì)量狀況，又可以提高監(jiān)測數(shù)據(jù)的有效性和代表性，并避免了重復(fù)設(shè)置帶來的經(jīng)濟投入，減少了相應(yīng)的工作量和經(jīng)費。

3 結(jié)語

根據(jù)物元分析法和聚類分析法分析的結(jié)果可以看出，物元分析法使用簡便，可直接編制程序上機運行，評價結(jié)果直觀、準(zhǔn)確、可靠，其結(jié)果與聚類分析法分析的結(jié)果近似，可用于水質(zhì)監(jiān)測優(yōu)化布點。但兩種方法都有一定的局限性，即兩種方法都以計算結(jié)果為分類依據(jù)(物元分析法以綜合關(guān)聯(lián)函數(shù)為依據(jù)，聚類分析法以歐式距離為依據(jù))，不能反映點位的空間關(guān)系與實際情況，在優(yōu)化點位時還需要根據(jù)實際情況具體分析，綜合考慮之后最終確定點位的優(yōu)化方案。

［1］ 蔡文.物元分析［M］.廣州:廣東高等教育出版社，1987.

［2］ 蔡文.物元模型及其應(yīng)用［M］.北京:科學(xué)技術(shù)文獻出版社，1994.

［3］ 孫中黨，畢雪玲，康莉萍，等.物元分析法在鄭州市大氣監(jiān)測優(yōu)化布點中的應(yīng)用［J］.河南科學(xué)，1999，11(2):194-199.

［4］ 弓曉峰，陳春麗，趙晉，等.鄱陽湖樂安河流域水質(zhì)監(jiān)測優(yōu)化布點［J］.湖泊科學(xué)，2006，18(5):545-549.

［5］ 高明惠.用物元分析進行水質(zhì)環(huán)境監(jiān)測優(yōu)化布點的研究［J］.環(huán)境科學(xué)進展，1997，5(3):77-81.

［6］ 朱慧君，于永斌.物元分析法在大氣監(jiān)測優(yōu)化布點中的應(yīng)用［J］.云南環(huán)境科學(xué)，1998，17(4):51-53.

［7］ 高明惠.物元分析在水質(zhì)營養(yǎng)化評價中的應(yīng)用［J］.環(huán)境科學(xué)進展，1995，3(5):43-49.

［8］ 翁建中，李繼影，梁柱，等.太湖藍藻水華時空分布與預(yù)警監(jiān)測響應(yīng)的分析［J］.環(huán)境監(jiān)控與預(yù)警，2010，2(3):1-4.

［9］ 國家環(huán)境保護局.HJ 494—2009 水質(zhì)采樣技術(shù)指導(dǎo)［S］.

［10］ 國家環(huán)境保護局.GB 3838—2002 地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)［S］.

［11］ 米紅，張文璋.實用現(xiàn)代統(tǒng)計分析方法及SPSS應(yīng)用［M］.北京:當(dāng)代中國出版社，2004.

Application of Matter-element Analysis in Optimized Points Selection of Water-quality Monitoring

CAO Yi
(Suqian Environmental Monitoring Central Station，Suqian，Jiangsu 223800，China)

Matter-element analysis is an effective method to manage incompatible problems.This article took the water-quality monitoring data of Zhushan Bay in July 2010 as example and analyzed the distribution of water quality monitoring point in Zhushan Bay with matter-element analysis and system cluster method.By comparing the results of two method，matter-element analysis method was more accurate and more credible.Combined the matter-element analysis with system cluster method，the distribution of water quality monitoring point was optimized，and the number of water quality monitoring point was reduced from 14 to 11.

water-quality monitoring;matter-element analysis;optimized points selection

X830.1

A

1674-6732(2012)-01-0043-04

10.3969/j.issn.1674-6732.2012.01.011

2010-11-20

曹毅(1985—)，男，助理工程師，本科，從事藍藻預(yù)警監(jiān)測工作。

(本欄目編輯 陸敏)