蔣艷君,謝悅波,黃　旻

(1.河海大學水文水資源學院,江蘇 南京　210098; 2.金華市浦江縣環(huán)境監(jiān)測站,浙江 金華　322200)

?

改進的物元分析法在水質監(jiān)測斷面布設優(yōu)化中的應用

蔣艷君1,謝悅波1,黃旻2

(1.河海大學水文水資源學院,江蘇 南京210098; 2.金華市浦江縣環(huán)境監(jiān)測站,浙江 金華322200)

摘要：依據(jù)傳統(tǒng)物元分析法的原理,使用層次分析法結合熵權法的組合權重法對物元分析法的權重計算進行改進,對關聯(lián)度進行聚類分析,使傳統(tǒng)的物元分析法得到改進與拓展。以浦陽江干流9個水質監(jiān)測斷面,主要監(jiān)測指標CODMn、COD、TP、NH3-N及氟化物的2012年監(jiān)測數(shù)據(jù)為基礎,采用改進的物元分析法對水質監(jiān)測斷面進行優(yōu)化,將原有9個斷面優(yōu)化為7個斷面,并對優(yōu)化結果進行一致性檢驗,結果顯示,改進的物元分析法優(yōu)化結果準確、有效。

關鍵詞：層次分析法;熵權法;組合權重法;斷面優(yōu)化;改進的物元分析法;聚類分析

水質監(jiān)測斷面的布設方案直接影響水質監(jiān)測結果,斷面優(yōu)化布設體現(xiàn)監(jiān)測工作的科學性、合理性,水質監(jiān)測斷面布設的總原則為以最少的斷面、測點數(shù),取得最全面、科學、合理的水質監(jiān)測信息[1]。傳統(tǒng)的水質監(jiān)測斷面常依據(jù)水資源保護管理的需要以及相關條例的規(guī)定而布設[2],具有較強的主觀性,并不能完全符合水質監(jiān)測斷面布設的總原則。可利用數(shù)理統(tǒng)計的方法對水質監(jiān)測斷面布設進行優(yōu)化,常用的數(shù)理統(tǒng)計方法有物元分析法[3-4]、歷史數(shù)據(jù)相關分析法[5]、模糊聚類分析法[6]、人工神經(jīng)網(wǎng)絡模型[7]等。

物元分析法能有效地解決多項水質指標優(yōu)選監(jiān)測斷面存在的不相容問題[8],筆者將層次分析法與熵權法結合,對物元分析法的權重計算進行改進,對關聯(lián)度進行聚類分析,提出利用改進的物元分析法對水質監(jiān)測斷面布設進行優(yōu)化。

1　改進的物元分析法

基于監(jiān)測斷面多項水質監(jiān)測指標的監(jiān)測值和將其標準化處理后得到的標準值,利用物元分析法建立系列物元矩陣進行對比分析,計算多項指標的綜合關聯(lián)函數(shù),分析其關聯(lián)度,劃分斷面的親近關系。

1.1構造物元矩陣和節(jié)域物元矩陣

根據(jù)各監(jiān)測斷面指標監(jiān)測數(shù)據(jù),計算各項監(jiān)測指標的平均值,構造平均值矩陣C,對于值越大表示水質越差的指標,分別取其平均值的最小值和最大值,構成最優(yōu)值矩陣A和最劣值矩陣B;對于值越小表示水質越差的指標,分別取其平均值的最小值和最大值的倒數(shù),構成相應矩陣。監(jiān)測指標的平均值、最優(yōu)值、最劣值構成的標準物元矩陣[8]表達式為

式中:RAC為最優(yōu)值標準物元矩陣;RCB為最劣值標準物元矩陣;MAC、MCB為標準事物;n為監(jiān)測指標數(shù);Qj為各項監(jiān)測指標;aj、bj、cj分別為各項監(jiān)測指標年均值的最小值、最大值、平均值。

由監(jiān)測指標的最大值與最小值構成的節(jié)域物元矩陣為

(3)

再將每一個采樣點的水質監(jiān)測數(shù)據(jù)作為一個物元,構成物元矩陣

(4)

式中:RAB為節(jié)域物元矩陣;Ri為物元矩陣;Mi為對象;xij為第i個監(jiān)測斷面各監(jiān)測指標的平均值。

1.2改進權重計算

筆者將計算主觀權重的層次分析法與計算客觀權重的熵權法結合,構成組合權重法對物元分析法的權重計算進行改進,組合權重法能很好地表達監(jiān)測指標對優(yōu)選結果的作用,同時反映出各監(jiān)測指標之間的相互作用,形成互補,使權重系數(shù)更具合理性。

1.2.1層次分析法

根據(jù)水質監(jiān)測指標與地表水環(huán)境質量標準的5級水質類別[9],構建水質的遞階層次結構[10],見圖1。

圖1　層次分析法結構

層次分析法步驟[11]如下:

a. 篩選水質監(jiān)測指標。筆者選取CODMn、COD、TP、NH3-N及氟化物作為主要監(jiān)測指標。

b. 構造判斷矩陣。以各準則層相對重要度構建判斷矩陣D=(dij)n×n,見表1,Bk為準則層中的指標準則(k=1,2,…,5),dCODMn、dCOD、dNH3-N、dTP、dF分別為CODMn、COD、TP、NH3-N及氟化物的監(jiān)測值。

表1　判斷矩陣

c. 求解最大特征值。采用特征向量法,使用Matlab軟件編程,對判斷矩陣的最大特征值λmax進行求解。

d. 一致性檢驗。采用判斷矩陣D除去最大特征值λmax的其他特征值的負均值,作為衡量判斷矩陣是否一致的數(shù)值指標,即一致性檢驗指標CI:

(5)

其中,CI值越小,表示一致性越好;反之越差。

1.2.2熵權法

熵權法具體步驟[12]如下:

c. 信息熵值計算。第j個監(jiān)測指標的信息熵值計算表達式為

(6)

(7)

d. 熵權計算。熵權計算表達式為

(8)

熵權向量為ω2=(ω21,ω22,…,ω2j,…,ω2n),其中ω2j為權重向量的分量。

1.2.3組合權重

筆者將層次分析法與熵權法結合,構成組合權重法對物元分析法賦權,計算式[13]為

(9)

式中:ω1i、ω2i分別是層次分析法與熵權法計算得到的權重向量的分量,i=1,2,…,5。

1.3關聯(lián)函數(shù)

關聯(lián)函數(shù)K(xij)用來描述監(jiān)測斷面物元與標準物元的符合程度[14],其含義為:①當K(xij)≥0時,表示比較對象符合標準對象要求,其值越大,符合程度越高;②當-1≤K(xij)<0時,表示比較對象不符合標準對象要求,但具備轉化為標準對象的條件,其絕對值越大,越容易轉化;③當K(xij)<-1時,表示比較對象不符合標準對象要求,且不具備轉化為標準對象的條件。

分別建立每一監(jiān)測斷面各項監(jiān)測指標對A、B的線性關聯(lián)函數(shù)[15]:

(10)

(11)

所有監(jiān)測指標對A、B的綜合關聯(lián)函數(shù)為

(12)

(13)

式中:KA(xij)、KB(xij)分別為第i個監(jiān)測斷面第j個監(jiān)測指標對A、B的線性關聯(lián)函數(shù);KA(xi)、KB(xi)分別為第i個監(jiān)測斷面對A、B的綜合關聯(lián)函數(shù);xij為斷面各項指標的監(jiān)測值;aj、bj、cj分別為監(jiān)測指標年均值的最小值、最大值與平均值;ωj為各監(jiān)測指標的權重系數(shù);xi為第i個監(jiān)測斷面監(jiān)測指標的最小值和最大值的平均值。以最優(yōu)、最劣綜合關聯(lián)函數(shù)值為基礎,采用SPSS軟件實現(xiàn)系統(tǒng)聚類分析過程[16]。

2　實例應用

2.1研究區(qū)域簡介

以浦陽江干流9個水質監(jiān)測斷面,2012年主要監(jiān)測指標CODMn、COD、TP、NH3-N及氟化物的年均值數(shù)據(jù)為基礎,上游至下游監(jiān)測斷面依次編號1～9號,位置見圖2。應用改進的物元分析法對9個斷面進行優(yōu)化分析,主要監(jiān)測指標的質量濃度年均值見表2。

圖2　浦陽江水系及斷面位置分布

2.2物元矩陣及權重確定

根據(jù)式(1)～(2),計算所得標準物元矩陣及節(jié)域物元矩陣如下,Q1為CODMn;Q2為COD;Q3為TP;Q4為NH3-N;Q5為氟化物。

(14)

(15)

(16)

根據(jù)式(16)計算得到9個監(jiān)測斷面的物元矩陣Ri(i=1,2,…,9),給出1號斷面的物元矩陣R1作為示例:

(17)

運用層次分析法及熵權法計算所得權重結果分別為:ω1=(0.202,0.717,0.049,0.008,0.024),ω2=(0.232,0.252,0.134,0.204,0.181)。組合權重為:ω=(0.261,0.513,0.098,0.049,0.079)。

2.3關聯(lián)函數(shù)值確定及聚類分析

根據(jù)式(14)～(15)計算9個斷面的線性關聯(lián)函數(shù)結果,見表3。

根據(jù)式(6)～(7),計算最優(yōu)、最劣綜合關聯(lián)函數(shù),結果見表4。

表4　最優(yōu)、最劣綜合關聯(lián)函數(shù)值

以最優(yōu)綜合關聯(lián)函數(shù)KA為縱坐標,以最劣綜合關聯(lián)函數(shù)KB為橫坐標,繪制每個監(jiān)測斷面關聯(lián)函數(shù)點聚圖,見圖3。

圖3　最優(yōu)、最劣綜合關聯(lián)函數(shù)點聚圖

將每個監(jiān)測斷面的最優(yōu)、最劣綜合關聯(lián)函數(shù)值作為監(jiān)測斷面系統(tǒng)聚類分析的基礎數(shù)據(jù),使用SPSS軟件對各監(jiān)測斷面進行系統(tǒng)聚類分析,結果見圖4。

圖4　系統(tǒng)聚類分析樹狀圖

根據(jù)水質監(jiān)測斷面的實際位置及檢測功能作用,本實例將9個水質監(jiān)測斷面分為3大類,聚類優(yōu)化結果見表5。

表5　聚類優(yōu)化結果

對于第一類的4個監(jiān)測斷面6、7、1、2,根據(jù)斷面實際地理位置及監(jiān)測功能分析,6、7號監(jiān)測斷面相距較近且監(jiān)測功能一致,1號監(jiān)測斷面位于流域入口且為監(jiān)測飲用水水源地之用,2號監(jiān)測斷面與6、7號監(jiān)測斷面相距較遠,故保留6、1、2號監(jiān)測斷面;對于第二類的4個監(jiān)測斷面5、9、4、3,5、4號監(jiān)測斷面相距較近檢測功能一致,9號監(jiān)測斷面位于支流與干流匯流處,故保留5、9號監(jiān)測斷面，3號斷面監(jiān)測功能與5、9號斷面不同，故作為保留；第三類的8號斷面作為保留。

將優(yōu)化前水質監(jiān)測數(shù)據(jù)與優(yōu)化后水質監(jiān)測數(shù)據(jù)作為兩個樣本,采用F檢驗與t檢驗對兩個樣本的均值做一致性檢驗,檢驗結果見表6。

表6　優(yōu)化前后斷面水質監(jiān)測數(shù)據(jù)及相關性檢驗結果

綜上所述,通過改進的物元分析法對浦江縣9個水質監(jiān)測斷面進行優(yōu)化分析,將9個監(jiān)測斷面最終優(yōu)化為7個監(jiān)測斷面1、2、3、5、6、8、9,刪減2個監(jiān)測斷面節(jié)約了資源與成本,實現(xiàn)斷面優(yōu)化布設的科學性。

3　結　語

筆者對傳統(tǒng)物元分析法的權重計算方法及關聯(lián)函數(shù)方面進行了改進,使傳統(tǒng)物元分析法得到拓展與改進,分析結果更為科學合理,適用于各種監(jiān)測斷面的優(yōu)化選擇。利用改進的物元分析法對浦江縣9個水質監(jiān)測斷面進行優(yōu)化選擇,結合主觀權重與客觀權重對各個斷面主要監(jiān)測指標進行賦權,使權重計算更為科學合理,使用SPSS軟件對綜合關聯(lián)函數(shù)值進行聚類分析,相比目測聚類更為科學。得到的優(yōu)化結果通過了一致性檢驗,改進的物元分析法為水質監(jiān)測斷面優(yōu)化提供了一種科學、實用、有效的方法。

參考文獻：

[1] 謝悅波.水信息技術[M].北京:中國水利水電出版社,2009.

[2] 樊引琴,李婳,劉婷婷,等.物元分析法在水質監(jiān)測斷面優(yōu)化中的應用[J].人民黃河,2012,34(11):82-84.(FAN Yinqin,LI Hua,LIU Tingting,et al.Application of matter element analysis method in optimized sections of water quality monitoring[J].Yellow River,2012,34(11):82-84.(in Chinese))

[3] 蔡文.物元分析[M].廣州:廣東高等教育出版社,1987.

[4] 朱慧君,于永斌.物元分析法在大氣監(jiān)測優(yōu)化布點中的應用[J].云南環(huán)境科學,1987，17(4):51-53.(ZHU Huijun,YU Yongbin.Application of matter element analysis to optimized sites of atmospheric monitoring[J].Yunnan Environmental Science,1987,17(4):51-53.(in Chinese))

[5] 趙吉國.東江流域地表水監(jiān)測省控斷面優(yōu)化布點[J].廣東水利水電,2004(6):27-29.(ZHAO Jiguo.Optimization of province controlled surface water monitoring section in the Dongjiang distribution[J].Guangdong Water Resources and Hydropower,2004(6):27-29.(in Chinese))

[6] 王晶,孫世群,邵超,等.模糊聚類分析法在水質監(jiān)測斷面優(yōu)化中的應用[J].廣州化工,2012,40(7):153-154,160.(WANG Jing,SUN Shiqun,SHAO Chao,et al.Fuzzy cluster analysis in the optimization of water quality monitoring sections[J].Guangdong Chemical Industry.2012,40(7):153-154,160.(in Chinese))

[7] 何璠.基于BP人工神經(jīng)網(wǎng)絡的環(huán)境質量評價模型研究[D].成都：四川大學,2006.

[8] 曹毅.物元分析法在水質監(jiān)測優(yōu)化布點中的應用[J].環(huán)境監(jiān)控與預警,2012,4(1):43-45,52.(CAO Yi.Application of metter-element analysis in optimized points selection of water quality monitoring[J].Environmental Monitoring and Forewarning,2012,4(1):43-45,52.(in Chinese))

[9] 中華人民共和國國家環(huán)境保護總局.GB3838—2002地表環(huán)境質量標準[S].北京: 中國環(huán)境科學出版社,2002.

[10] 鄧雪,李家銘,曾浩健,等.層次分析法權重計算方法分析及其應用研究[J].數(shù)學的實踐與認識,2012,42(7): 93-100.(DENG Xue,LI Jiaming,ZENG Haojian,et al.Research on computation methods of AHP wight vector and its applications[J].Mathematics in Practice and Theory,2012,42(7): 93-100.(in Chinese))

[11] 盧文喜,李迪,張蕾,等.基于層次分析法的模糊綜合評價在水質評價中的應用[J].節(jié)水灌溉,2011(3):43-46.(LU Wenxi,LI Di,ZHANG Lei,et al.Application of fuzzy comprehensive evaluation based on AHP in water quality evaluation[J].Water Saving and Irrigation,2011(3):43-46.(in Chinese))

[12] 馮斌.模糊綜合評價的熵權法在地下水水質評價中的應用[J].中國水運,2012,12(7): 186-187.(FENG Bin.Fuzzy comprehensive evaluation of entropy method in groundwater quality evaluation[J].China Water Transport,2012,12(7): 186-187.(in Chinese))

[13] 喬雨,梁秀娟,王宇博,等.組合權重模糊數(shù)學法在水質評價中的應用[J].人民黃河,2015,37(5):77-79.(QIAO Yu,LIANG Xiujuan,WANG Yubo,et al.Evaluation of groundwater quality using combined weight fuzzy mathematical method[J].Yellow River,2015,37(5):77-79.(in Chinese))

[14] 孫學穎,唐德善,張新嬌.改進物元分析模型在北之江流域水質評價中的應用[J].南水北調與水利科技,2014,12(3):55-58.(SUN Xueying,TANG Deshan,ZHANG Xinjiao.Application of improved metter-element model in water quality assessment of Beizhijiang Drainage Basin[J],South-to-North Water Transfers and Water Science & Technology,2014,12(3):55-58.(in Chinese))

[15] 馬飛,蔣莉.河流水質監(jiān)測斷面優(yōu)化設置研究: 以南運河為例[J].環(huán)境科學與管理,2006,31(8):171-172.(MA Fei,JIANG Li.Research of optimized setting of water quality monitoring sections in rivers: take South Canal as an example[J].Environmental Science and Management,2006,31(8):171-172.(in Chinese))

[16] 楊靜.改進的模糊綜合評價法在水質評價中的應用[D].重慶:重慶大學,2014.

DOI：10.3880/j.issn.1004-6933.2016.04.023

作者簡介:蔣艷君(1992—),女,碩士研究生,研究方向為水生態(tài)修復。E-mail:640483090@qq.com 通信作者:謝悅波,教授。E-mail:xyb@hhu.edu.cn

中圖分類號：TV213.4;X832

文獻標志碼：A

文章編號：1004-6933(2016)04-0136-06

(收稿日期：2015-09-13編輯：王芳)

Application of improved matter element analysis method to optimization of water quality monitoring sections

JIANG Yanjun1, XIE Yuebo1, Huang Min2

(1. College of Hydrology and Water Resources, Hohai University, Nanjing 210098, China;2.EnvironmentalMonitoringStationofPujiangCounty,JinhuaCity,Jinhua322200,China)

Abstract：According to the principles of traditional matter element analysis methods, the combined weight method, which is the combination of the AHP method and the entropy method, was used to modify the weights in the matter element analysis method. Then, the relevance was studied through cluster analysis, and thus the traditional matter element analysis method was improved and expanded. Based on monitoring data involving CODMn, COD, TP, NH3-N, and fluoride, at nine water quality monitoring sections on the main stream of the Puyang River in 2012, the improved matter element analysis method was used to optimize these water quality monitoring sections. As a result, the original nine sections were decreased to seven sections after optimization. The results were verified by the consistency test, indicating that the improved matter element analysis method can provide accurate and effective optimization results.

Key words：AHP; entropy method; combined weight method; section optimization; improved matter element analysis method; cluster analysis