魯婷婷,俞 淞,王紅瑞,丁建新

(北京師范大學(xué)水科學(xué)研究院,北京 100875)

地下水是水資源的重要組成部分,對(duì)經(jīng)濟(jì)和社會(huì)的可持續(xù)發(fā)展具有不可替代的作用,特別是在地表水缺少或污染嚴(yán)重的地區(qū),地下水常常被用作人們的飲用水源。地下水污染趨勢(shì)日益嚴(yán)重,正面臨著由點(diǎn)到面、由淺到深、由城市到農(nóng)村不斷擴(kuò)展的趨勢(shì)[1]。2011年,中國(guó)工程院與中國(guó)環(huán)保部聯(lián)合發(fā)布的中國(guó)環(huán)境宏觀戰(zhàn)略研究報(bào)告指出,中國(guó)地下水污染嚴(yán)重,其中118個(gè)大中城市地下水已普遍遭到污染,嚴(yán)重污染的城市占64%,平原區(qū)約有54%的地下水不符合生活飲用水水質(zhì)標(biāo)準(zhǔn)。

近幾年,隨著科學(xué)技術(shù)的不斷發(fā)展,地下水質(zhì)量綜合評(píng)價(jià)方法也越來越多,主要包括單因子評(píng)價(jià)法[2]、灰色聚類法[3]、神經(jīng)網(wǎng)絡(luò)法[4]和支持向量機(jī)[5]等,這些方法都取得了良好的效果。在1994年實(shí)施的GB/T 14848—93《地下水質(zhì)量標(biāo)準(zhǔn)》中明確給出地下水質(zhì)量分類指標(biāo),但由于影響地下水質(zhì)量的因素眾多,評(píng)價(jià)指標(biāo)與評(píng)價(jià)因子之間仍存在模糊、非線性的復(fù)雜關(guān)系,使得各種方法在應(yīng)用中有其優(yōu)點(diǎn)和弊端。例如單因子評(píng)價(jià)法,評(píng)價(jià)方法簡(jiǎn)單、明了,能清晰的判斷出主要污染因子及其分布區(qū)域,但不能全面反應(yīng)水質(zhì)的整體狀況。模糊綜合評(píng)價(jià)法能充分考慮邊界的模糊性,但是它的計(jì)算比較復(fù)雜,在評(píng)價(jià)面積較大、評(píng)價(jià)對(duì)象數(shù)量較多的情況下,應(yīng)用比較困難,需要借助計(jì)算機(jī)輔助編程完成[6]。支持向量機(jī)雖然評(píng)價(jià)精度高,能簡(jiǎn)化評(píng)價(jià)過程,但在相關(guān)參數(shù)選擇方面有待完善[7]。為解決水質(zhì)評(píng)價(jià)中評(píng)價(jià)指標(biāo)與評(píng)價(jià)因子之間模糊、非線性的復(fù)雜關(guān)系,將趙克勤[8]提出的集對(duì)分析理論應(yīng)用于地下水水質(zhì)評(píng)價(jià)中。

1 基于集對(duì)分析的地下水質(zhì)量評(píng)價(jià)模型的建立

1.1 集對(duì)分析基本原理

集對(duì)分析是趙克勤于1989年在包頭召開的全國(guó)系統(tǒng)理論會(huì)上提出的一種新型的處理模糊和不確定知識(shí)的系統(tǒng)分析方法,目前已廣泛應(yīng)用于地表水環(huán)境質(zhì)量評(píng)價(jià)中[9],其核心思想是對(duì)不確定性系統(tǒng)中有一定聯(lián)系的2個(gè)集合構(gòu)建集對(duì),再對(duì)集對(duì)的特性做同一性、差異性、對(duì)立性分析,并加以度量刻畫,建立集對(duì)的同異反聯(lián)系度表達(dá)式。

設(shè)有聯(lián)系的集合X和Y,X有n項(xiàng)表征其特征,即X=(x1,x2,…，xn),Y亦有n項(xiàng)表征其特征,即Y=(y1,y2,…，yn)。X和Y構(gòu)成集對(duì)H(X,Y)。描述H(X,Y)間關(guān)系的聯(lián)系度定義為

(1)

式中：μX～Y為集合X和Y的聯(lián)系度,S為集對(duì)H(X,Y)中特性同一的個(gè)數(shù);F為特性差異的個(gè)數(shù);P為特性對(duì)立的個(gè)數(shù);S+F+P=n。I為差異度系數(shù),取值[-1,1]。當(dāng)I取1時(shí),差異性轉(zhuǎn)化為同一性；當(dāng)I取-1時(shí),差異性轉(zhuǎn)化為對(duì)立性；當(dāng)I在(-1,1)取值時(shí),反映了同一性與對(duì)立性所占的比例,I越趨于0,不確定性越明顯。J為對(duì)立度系數(shù),規(guī)定其恒取值為-1。

μX～Y=a+bI+cJ

(2)

式中:a,b,c滿足a+b+c=1。a,b,c分別為集對(duì)H(X,Y)的同一度、差異度和對(duì)立度。當(dāng)a越接近于1時(shí),表示這2個(gè)集合的特性越傾向于同一;當(dāng)c越接近于1時(shí),表示這2個(gè)集合的特性越傾向于對(duì)立;當(dāng)b越接近于1時(shí),表示這2個(gè)集合的特性越傾向于差異(既不同一也不對(duì)立)。

在地下水質(zhì)量評(píng)價(jià)中,地下水質(zhì)量分類及分類指標(biāo)是確定的,實(shí)際采用的評(píng)價(jià)因子是不確定的,基于集對(duì)分析的原理,將實(shí)際采用的評(píng)價(jià)因子與地下水評(píng)價(jià)標(biāo)準(zhǔn)之間構(gòu)成一個(gè)集對(duì),通過計(jì)算2個(gè)集合的聯(lián)系度,評(píng)判地下水的環(huán)境質(zhì)量。

1.2 聯(lián)系度的確定

聯(lián)系度是對(duì)不確定系統(tǒng)中有一定聯(lián)系的2個(gè)集合關(guān)系的量化,能清晰定量的顯示2個(gè)集合之間的復(fù)雜關(guān)系。在確定聯(lián)系度之前,首先要進(jìn)行單因子評(píng)價(jià),將樣本的實(shí)測(cè)值與地下水質(zhì)量標(biāo)準(zhǔn)進(jìn)行對(duì)比,分別得出每個(gè)指標(biāo)所屬的水質(zhì)級(jí)別,然后再計(jì)算單個(gè)指標(biāo)所處的水質(zhì)級(jí)別與其他水質(zhì)級(jí)別的聯(lián)系度。確定準(zhǔn)確的聯(lián)系度是進(jìn)行集對(duì)分析的關(guān)鍵,具體確定方法是:將單因子評(píng)價(jià)的結(jié)果與5類水質(zhì)標(biāo)準(zhǔn)進(jìn)行對(duì)比,若單因子評(píng)價(jià)結(jié)果與該評(píng)價(jià)級(jí)別相同,聯(lián)系度為1;當(dāng)單因子評(píng)價(jià)結(jié)果與該評(píng)價(jià)級(jí)別相隔,聯(lián)系度為-1;當(dāng)單因子評(píng)價(jià)結(jié)果與該評(píng)價(jià)級(jí)別相鄰,根據(jù)評(píng)價(jià)指標(biāo)與評(píng)價(jià)標(biāo)準(zhǔn)距離遠(yuǎn)近確定聯(lián)系度。具體量化公式為

a.樣本i第j項(xiàng)評(píng)價(jià)指標(biāo)與Ⅰ級(jí)水質(zhì)聯(lián)系度

(3)

式中：xij為樣本i第j項(xiàng)評(píng)價(jià)指標(biāo)實(shí)測(cè)值;s1j,s2j為各評(píng)價(jià)標(biāo)準(zhǔn)中分類指標(biāo)的限值。

b.樣本i第j項(xiàng)評(píng)價(jià)指標(biāo)與Ⅱ級(jí)水質(zhì)聯(lián)系度

(4)

式中：s3j為評(píng)價(jià)標(biāo)準(zhǔn)中分類指標(biāo)的限值。

c.樣本i第j項(xiàng)評(píng)價(jià)指標(biāo)與Ⅲ級(jí)水質(zhì)聯(lián)系度

(5)

式中：s4j為評(píng)價(jià)標(biāo)準(zhǔn)中分類指標(biāo)的限值。

d.樣本i第j項(xiàng)評(píng)價(jià)指標(biāo)與Ⅳ級(jí)水質(zhì)聯(lián)系度

(6)

e.樣本i第j項(xiàng)評(píng)價(jià)指標(biāo)與Ⅴ級(jí)水質(zhì)聯(lián)系度

(7)

式中：s5j為各評(píng)價(jià)標(biāo)準(zhǔn)中分類指標(biāo)的限值，是一個(gè)不確定的值,為了便于計(jì)算樣本單因子位于Ⅴ級(jí)時(shí)與Ⅳ級(jí)評(píng)價(jià)標(biāo)準(zhǔn)之間的聯(lián)系度,將s5j賦值為2s4j-s3j。

1.3 確定指標(biāo)權(quán)重的熵權(quán)法

在信息論中,熵是對(duì)不確定性的一種度量,信息量越大,不確定性就越小,熵也就越小;信息量越小,不確定性越大,熵也就越大。根據(jù)熵的特性,可以根據(jù)評(píng)價(jià)指標(biāo)值構(gòu)成的矩陣判斷各指標(biāo)的權(quán)重,避免人為給予權(quán)重的主觀性。主要的計(jì)算步驟如下:

a.假定有m個(gè)評(píng)價(jià)對(duì)象,每個(gè)評(píng)價(jià)對(duì)象有n個(gè)評(píng)價(jià)指標(biāo),構(gòu)成矩陣R

R=(Cij)m×n

i=(1,2,…,m),j=(1,2,…,n)

(8)

式中:Cij是第i個(gè)評(píng)價(jià)對(duì)象的第j個(gè)評(píng)價(jià)指標(biāo)的實(shí)測(cè)值。

b.將判斷矩陣R歸一化,得到歸一化矩陣B,歸一化公式:

(9)

式中:bij為矩陣B中的元素；Cmax，Cmin分別為同指標(biāo)下,不同方案中最滿意者和最不滿意者的值(本文評(píng)價(jià)為越小越滿意)。

c.根據(jù)傳統(tǒng)的熵的概念可以定義各評(píng)價(jià)指標(biāo)的熵為

(i=1,2,…,m;j=1,2,…，n)

(10)

d.計(jì)算各指標(biāo)的熵權(quán)

W=(wj)1×n

(11)

式中:wj為第j個(gè)評(píng)價(jià)指標(biāo)的權(quán)重;wt為第t個(gè)評(píng)價(jià)指標(biāo)的權(quán)重。

1.4 地下水環(huán)境質(zhì)量等級(jí)的確定

根據(jù)式(3)～(7),計(jì)算樣本與各評(píng)價(jià)等級(jí)之間的聯(lián)系度,并加權(quán)處理:

(12)

式中:uik為樣本i對(duì)評(píng)價(jià)等級(jí)k級(jí)之間的聯(lián)系度;uijk為樣本i中評(píng)價(jià)指標(biāo)j對(duì)k級(jí)評(píng)價(jià)標(biāo)準(zhǔn)的聯(lián)系度。

依據(jù)集對(duì)分析基本原理,選取聯(lián)系度最大的評(píng)價(jià)等級(jí)為該樣本地下水環(huán)境質(zhì)量等級(jí),即

(13)

判定地下水水質(zhì)級(jí)別為p級(jí)。

2 實(shí)例應(yīng)用

為了驗(yàn)證集對(duì)分析方法的合理性,本文采用山東省聊城市東昌府區(qū)5個(gè)辦事處地下水監(jiān)測(cè)數(shù)據(jù)[10],以鐵、錳、總硬度、氯化物、高錳酸鹽指數(shù)和硫酸鹽為評(píng)價(jià)指標(biāo),以GB/T14848—93《地下水水質(zhì)標(biāo)準(zhǔn)》為水質(zhì)分級(jí)原則,進(jìn)行地下水水質(zhì)評(píng)價(jià)。地下水環(huán)境質(zhì)量標(biāo)準(zhǔn)見表1。

表1 地下水環(huán)境質(zhì)量標(biāo)準(zhǔn) mg/L

現(xiàn)以東昌府區(qū)新區(qū)辦事處為例,根據(jù)水質(zhì)監(jiān)測(cè)數(shù)據(jù)和地下水環(huán)境質(zhì)量標(biāo)準(zhǔn),確定單因子評(píng)價(jià)結(jié)果和聯(lián)系度,見表2和表3。

表2 新區(qū)辦事處單因子評(píng)價(jià)結(jié)果

表3 各評(píng)價(jià)因子與評(píng)價(jià)等級(jí)之間的聯(lián)系度

根據(jù)熵權(quán)計(jì)算公式,得出各指標(biāo)在地下水環(huán)境質(zhì)量評(píng)價(jià)中的權(quán)重分別為(0.20,0.18,0.13,0.16,0.15,0.18)。由式(12)計(jì)算出新區(qū)辦事處地下水水質(zhì)與Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ 5類評(píng)價(jià)等級(jí)之間的綜合聯(lián)系度分別為(-0.55,-0.14,-0.08,-0.06,-0.37),由式(13)可判定該地區(qū)地下水水質(zhì)屬于Ⅳ類。同理可得北城辦事處、柳園辦事處、鳳凰辦事處、湖西辦事處地下水水質(zhì)與各評(píng)價(jià)等級(jí)之間的綜合聯(lián)系度分別為(-0.61,-0.11,-0.05,-0.60,-0.34)、(-0.47,-0.21,0.10,-0.06,-0.63)、(-0.34,-0.17,-0.03,0.05,-0.63)、(0.31,0.37,-0.31,-0.85,-1)。相應(yīng)的地下水水質(zhì)等級(jí)分別為:Ⅲ類,Ⅲ類,Ⅳ類,Ⅱ類。

3 評(píng)價(jià)結(jié)果的比較與分析

根據(jù)上述方法,將所得的結(jié)果與文獻(xiàn)[11]中模糊綜合評(píng)價(jià)結(jié)果進(jìn)行比較,評(píng)價(jià)結(jié)果見表4.

表4 集對(duì)分析方法與模糊綜合評(píng)價(jià)法評(píng)價(jià)結(jié)果比較

從表4可以看出,新區(qū)辦事處、柳園辦事處和湖西辦事處2種評(píng)價(jià)方法的評(píng)價(jià)結(jié)果相同,鳳凰辦事處水質(zhì)評(píng)價(jià)結(jié)果相差一級(jí),北城辦事處評(píng)價(jià)結(jié)果相差較大。產(chǎn)生這種差異的原因可能有2個(gè)方面：第一方面是2種評(píng)價(jià)方法之間的差異。2種方法都是通過建立函數(shù)進(jìn)行質(zhì)量評(píng)價(jià),但模糊綜合評(píng)價(jià)法是通過計(jì)算隸屬度確定評(píng)價(jià)等級(jí),隸屬度在[0,1]取值,沒有對(duì)水質(zhì)等級(jí)之間的過渡性進(jìn)行詳細(xì)描述,而集對(duì)分析方法是通過計(jì)算聯(lián)系度確定評(píng)價(jià)等級(jí),不僅將跨越區(qū)間擴(kuò)大至[-1,1],而且從同一性、對(duì)立性、差異性3個(gè)方面進(jìn)行分析,特別是將差異度引入函數(shù)中,描述水質(zhì)等級(jí)之間的模糊過渡性,使評(píng)價(jià)結(jié)果更加精細(xì)、合理。另一方面的原因是評(píng)價(jià)過程中權(quán)重確定的方法不同。在模糊綜合評(píng)價(jià)過程中,權(quán)重的求法是根據(jù)該地區(qū)指標(biāo)的實(shí)測(cè)值與該指標(biāo)標(biāo)準(zhǔn)限值的算術(shù)平均值的比值求得。在計(jì)算過程中,權(quán)重的大小只與該指標(biāo)實(shí)測(cè)值有關(guān),這就會(huì)出現(xiàn)因?yàn)槟骋粋€(gè)指標(biāo)實(shí)測(cè)值偏高而導(dǎo)致評(píng)價(jià)結(jié)果偏高的現(xiàn)象,如北城辦事處鐵值明顯高于其他地區(qū),最后的結(jié)果也使模糊綜合評(píng)價(jià)法的結(jié)果高于集對(duì)分析評(píng)價(jià)方法結(jié)果。而且模糊綜合評(píng)價(jià)法只是根據(jù)一個(gè)單一的數(shù)值確定指標(biāo)的權(quán)重,誤差較大。本文在確定權(quán)重的過程中引進(jìn)信息論中熵的概念,根據(jù)評(píng)價(jià)指標(biāo)值構(gòu)成的矩陣來判斷各指標(biāo)的權(quán)重,權(quán)重的大小與每個(gè)地區(qū)的實(shí)測(cè)值都存在關(guān)系,既減少了單一數(shù)值確定權(quán)重產(chǎn)生的誤差,又能夠反應(yīng)各指標(biāo)的集中趨勢(shì)。

4 結(jié) 語(yǔ)

本文通過構(gòu)建評(píng)價(jià)因子的實(shí)測(cè)值與水質(zhì)評(píng)價(jià)標(biāo)準(zhǔn)之間的集對(duì)關(guān)系計(jì)算聯(lián)系度來確定地下水的質(zhì)量等級(jí)。在計(jì)算過程中,引用信息論中的熵值理論確定各個(gè)評(píng)價(jià)指標(biāo)的權(quán)重,避免人為賦予權(quán)重的主觀性,使評(píng)價(jià)結(jié)果更趨合理性。集對(duì)分析方法的主要優(yōu)點(diǎn)在于不僅能夠通過構(gòu)建同、異、反聯(lián)系度將評(píng)價(jià)過程量化,而且評(píng)價(jià)的原理易懂、計(jì)算簡(jiǎn)便、結(jié)果直觀,適合大部分地區(qū)的地下水的質(zhì)量評(píng)價(jià)。但是集對(duì)分析方法在評(píng)價(jià)過程中也有不成熟的地方,例如水質(zhì)分級(jí)時(shí)s5的確定,本文雖然根據(jù)水質(zhì)分級(jí)的一般規(guī)律將s5賦值為2s4j-s3j,但水質(zhì)分級(jí)的指標(biāo)限值并不是標(biāo)準(zhǔn)線性關(guān)系,在確定過程中存在一定的主觀性。這也需要對(duì)集對(duì)分析方法不斷改進(jìn)優(yōu)化,使其能客觀有效地對(duì)地下水質(zhì)量進(jìn)行評(píng)價(jià)。

[1] 羅蘭.我國(guó)地下水污染現(xiàn)狀與防治對(duì)策研究[J].中國(guó)地質(zhì)大學(xué)學(xué)報(bào):社會(huì)科學(xué)版,2008,8(2):72-75.

[2] 胡浩云,王海峰,劉新霞.邯鄲市地下水水質(zhì)現(xiàn)狀分析[J].人民長(zhǎng)江,2005,36(10):31-32.

[3] 張祖亮.灰色聚類法在地下水質(zhì)量評(píng)價(jià)中的應(yīng)用[J].云南環(huán)境科學(xué),2004,23(1):60-62.

[4] 羅定貴,王學(xué)軍,郭青.基于MATLAB實(shí)現(xiàn)ANN方法在地下水質(zhì)評(píng)價(jià)中的應(yīng)用[J].北京大學(xué)學(xué)報(bào):自然科學(xué)版,2004,40(2):296-302.

[5] 徐勁力.支持向量機(jī)在水質(zhì)評(píng)價(jià)中的應(yīng)用[J].中國(guó)農(nóng)村水利水電,2007(3):7-9.

[6] 張新鈺,辛寶東,劉文臣,等.三種地下水水質(zhì)評(píng)價(jià)方法的應(yīng)用對(duì)比分析[J].城市地質(zhì),2011,6(1):30-36.

[7] CHAPELLE O, VAPNIK V.Choosing multiple parameters for support vector machines [J].Machine Learning, 2002, 46(1):131-159.

[8] 趙克勤.集對(duì)分析及其初步應(yīng)用[J].大自然探索,1994,13(47):67-72.

[9] 高軍省.湖泊富營(yíng)養(yǎng)化綜合評(píng)價(jià)的五元聯(lián)系數(shù)法[J].人民長(zhǎng)江,2010,41(21):81-84.

[10] 彭小金,張艷紅,李輝輝.模糊綜合評(píng)價(jià)在地下水質(zhì)評(píng)價(jià)中的應(yīng)用[J].河南水利與南水北調(diào),2009(1):34-35,37.