金菊良，湯 睿，周戎星，王友貞，夏小林，張浩宇

(1.合肥工業(yè)大學土木與水利工程學院,安徽 合肥 230009;2.合肥工業(yè)大學水資源與環(huán)境系統(tǒng)工程研究所,安徽 合肥 230009;3.安徽省·水利部淮河水利委員會水利科學研究院水利水資源安徽省重點實驗室,安徽 蚌埠 233000)

水生態(tài)文明建設是人類生產(chǎn)生活、社會經(jīng)濟發(fā)展的重要保證和基礎[1-3]。城市水生態(tài)文明建設評價和診斷研究的現(xiàn)有基礎較為薄弱，尚缺乏統(tǒng)一的理論框架和有效的定量方法[4-5]，其主要問題有：一是如何合理構(gòu)建適用于不同城市水生態(tài)文明建設評價的指標體系和定量評價模型；二是如何診斷識別城市水生態(tài)文明建設中的脆弱性指標和判別建設狀態(tài)的發(fā)展趨勢。目前城市水生態(tài)文明建設評價方法和模型主要有主成分分析法[6]、模糊物元模型[7]、隨機森林回歸算法[8]、模糊綜合評價法[9]、投影尋蹤法[10]、壓力-狀態(tài)-響應概念模型[11]、物元可拓模型[12]等，這些方法和模型在評價指標值與評價等級之間的不確定性處理方面仍有改進空間，有關(guān)城市水生態(tài)文明建設的定量評價研究成果還很少。加速遺傳算法與模糊層次分析法相結(jié)合的遺傳模糊層次分析法(accelerating genetic algorithm based fuzzy analytic hierarchy process， AGA-FAHP)是由金菊良等[13]提出的一種改進的層次分析法，有效地克服了傳統(tǒng)層次分析法中構(gòu)造、檢驗和修正判斷矩陣的一致性問題，計算結(jié)果較為客觀穩(wěn)定，已在資源、環(huán)境[14]等領(lǐng)域取得較好的應用效果。集對分析是由趙克勤[15]基于事物的普遍聯(lián)系原理和對立統(tǒng)一的關(guān)系結(jié)構(gòu)提出的一種不確定性定量分析理論，通過構(gòu)造聯(lián)系數(shù)可以從同、異、反三方面定量地分析系統(tǒng)確定性與不確定性特征關(guān)系，使評價結(jié)果更為全面精確，已經(jīng)在社會、經(jīng)濟、信息管理、資源與環(huán)境科學[15-18]等領(lǐng)域廣泛應用。本文在實際調(diào)研、深入分析城市區(qū)域水資源-生態(tài)環(huán)境-經(jīng)濟社會復雜系統(tǒng)并結(jié)合專家咨詢意見的基礎上，采用AGA-FAHP篩選城市水生態(tài)文明建設評價指標，采用基于聯(lián)系數(shù)的方法對城市水生態(tài)文明建設進行定量評價，并通過安徽省4城市評價驗證了方法的實用性。

1 評價指標篩選

采用AGA-FAHP篩選評價指標的步驟如下：

步驟1：通過咨詢專家、比較各要素(子系統(tǒng)或指標)兩兩之間的重要性，得到模糊互補判斷矩陣A=(aij)n×n，其中n為要素的個數(shù)，0≤aij≤1，aij+aji=1，aij為要素i的重要性高于要素j的程度，aji為要素j的重要性高于要素i的程度[13,19-21]。

步驟2：檢驗模糊互補判斷矩陣的一致性，計算各要素的權(quán)重。若模糊互補判斷矩陣A滿足加性傳遞性：

(aik-0.5)+(akj-0.5)=aij-0.5
(i,j,k=1,2,…,n)

(1)

則根據(jù)式(1)可構(gòu)造模糊互補判斷矩陣A的一致性指標ρ：

(2)

若ρ≤0.2，則判斷矩陣A具有滿意的一致性[19-20]。

(3)

當式(3)左端的值小于或等于0.1時，則判斷矩陣A具有滿意的一致性[13]。當矩陣A的一致性不滿足要求時，就需要修正A，假設修正后的矩陣為B=(bij)n×n，矩陣B的各要素權(quán)重仍表示為{wk|k=1，2，…，n}。于是可用下式得到A的最優(yōu)模糊一致性判斷矩陣B和指標權(quán)重wk[13]：

(4)

式中：目標函數(shù)f(n)稱為一致性指標系數(shù)(consistency index coefficient， CIC)；d為非負參數(shù)，經(jīng)驗取值范圍為[0,0.5]。式(4)用加速遺傳算法(accelerated genetic algorithm， AGA)[22]求解。

步驟3：根據(jù)權(quán)重wk的大小，篩選掉權(quán)重相對較小的要素，所剩下的要素組成城市水生態(tài)文明建設評價指標體系。

2 評價方法與步驟

為對不同城市水生態(tài)文明建設進行定量評價，針對評價指標隸屬于評價等級的不確定性問題(如原有評價方法缺陷等)，在篩選得到的評價指標體系的基礎上，先計算不同城市水生態(tài)文明建設的聯(lián)系數(shù)，再綜合級別特征值法、屬性識別法、加權(quán)綜合得分法和減法集對勢法的評價結(jié)果得到城市水生態(tài)文明建設評價結(jié)果。該方法具體評價步驟如下：

步驟1：評價指標評分。已知城市水生態(tài)文明建設評價指標體系為{Cj|j=1,2,…,m1}，其中Cj為指標體系中第j個評價指標，m1為評價指標數(shù)目，根據(jù)各指標的變化范圍及與評價目標的作用關(guān)系可直接對評價指標在0～100分之間評分，得到的評價指標評分樣本數(shù)據(jù)集合為{xij|i=1,2,…,m2;j=1,2,…,m1}，其中m2為評價樣本數(shù)目。

步驟2：構(gòu)建城市水生態(tài)文明建設評價等級標準。根據(jù)DB 34/T 3321—2019《水生態(tài)文明城市評價準則》，本文采用5個評價等級(l=5)[23]：得分(90,100]為1級，(80,90]為2級，(60,80]為3級，(40,60]為4級，(0,40]為5級，等級越高則越不和諧。

步驟3：應用AGA-FAHP確定評價指標的權(quán)重{wj|j=1,2,…,m1}。

步驟4：計算城市水生態(tài)文明建設評價樣本的聯(lián)系數(shù)[15-17]：

ui=vi1+vi2I1+vi3I2+vi4I3+vi5J

(5)

(6)

步驟5：計算級別特征值法[25]等級值：

(7)

式中：hi為用級別特征值法計算的樣本i的評價等級值。根據(jù)均分原則，1.0≤hi<1.8為1級，1.8≤hi<2.6為2級，2.6≤hi<3.4為3級，3.4≤hi<4.2為4級，4.2≤hi≤5.0為5級，等級越高代表越不和諧。

步驟6：計算屬性識別法[26]等級值：

(8)

式中:gi為用屬性識別法計算的樣本i的等級值；λ為置信度，一般可在[0.50,0.70]內(nèi)取值，λ越大說明置信水平越高[26]；l為評價標準等級數(shù)目。gi代表城市水生態(tài)文明建設屬性識別法中樣本i的等級值。

步驟7：計算加權(quán)綜合得分法[27]評分值：

(9)

式中：yi為綜合評分值；wj為指標j的權(quán)重；xij為樣本i指標j的評分值，根據(jù)步驟2的分級標準可判斷樣本i所屬的等級，得分越高越和諧。

步驟8：用減法集對勢法[17]識別、判斷城市水生態(tài)文明建設脆弱性指標和建設狀態(tài)的發(fā)展趨勢[28-30]。五元聯(lián)系數(shù)ui的減法集對勢可表達為[17,31]

sf(ui)=(vi1-vi5)(1+vi2+vi3+vi4)+

0.5(vi2-vi4)(vi2+vi3+vi4)

(10)

基于最大信息熵的均分原則可把減法集對勢的變化范圍均等劃分為5個等級[17,31]：-1≤sf(u)<0.6為5級，-0.6≤sf(u)<-0.2為4級，-0.2≤sf(u)≤0.2是3級，0.2

步驟9：將上述4種方法計算出的結(jié)果進行對比分析，定性綜合得出當前城市水生態(tài)文明建設的評價結(jié)果。

3 實例驗證

選取馬鞍山、淮北、六安和池州4座具有代表性的安徽省城市運用上述方法進行實例分析驗證。馬鞍山市地處安徽東部，交通便捷，是一座工業(yè)發(fā)達的沿江城市；淮北市是皖北礦產(chǎn)資源豐富的能源城市；六安市位于皖西，處在江淮之間，生態(tài)環(huán)境良好，是安徽省文明城市；池州市位于安徽南部，森林覆蓋率高，自然資源豐富，旅游業(yè)發(fā)達。4市在地理分布、自然條件和經(jīng)濟發(fā)展上均具有一定的代表性，本文以4市2018年的水生態(tài)文明建設評價為例進行分析。

3.1 評價結(jié)果

根據(jù)指標體系構(gòu)建原則[2-4]，從城市水生態(tài)文明建設評價因子系統(tǒng)與研究區(qū)域社會經(jīng)濟系統(tǒng)交互作用角度，基于所選4市的水生態(tài)文明建設實際情況，構(gòu)建由水安全、水生態(tài)、水管理、水景觀、水文化5個子系統(tǒng)、各子系統(tǒng)下若干二級評價指標組成的安徽省城市水生態(tài)文明建設初步評價指標體系。用AGA-FAHP對初步評價指標體系進行篩選，得到安徽省城市水生態(tài)文明建設評價指標體系如圖1所示(括號中數(shù)值為指標權(quán)重)。

4市2018年評價指標評分結(jié)果見表1[4,10]，評分細則參考安徽省水資源公報等。由式(6)計算得水生態(tài)文明建設評價樣本聯(lián)系數(shù)分量，再將評價年的聯(lián)系數(shù)分量分別代入式(7)～(10)，根據(jù)表1可求得馬鞍山、淮北、六安和池州4市2018年水生態(tài)文明建設評價的級別特征值、屬性識別值(置信度取0.70)[26,32]、加權(quán)綜合得分法評分值以及相應的減法集對勢值如表2所示。

表1 安徽省4市水生態(tài)文明建設評價指標評分結(jié)果

表2 安徽省4市水生態(tài)文明建設評價的聯(lián)系數(shù)分量和評價結(jié)果

3.2 分析與討論

a.從級別特征值法的評價結(jié)果看，馬鞍山、淮北、六安、池州市2018年水生態(tài)文明建設分別處于2級、1級、1級和1級等級狀態(tài)；從屬性識別法的評價結(jié)果看，馬鞍山、淮北、六安、池州市2018年水生態(tài)文明建設分別處于2級、1級、1級、1級等級狀態(tài)；從4市2018年水生態(tài)文明建設的加權(quán)綜合得分法評分結(jié)果來看，馬鞍山處于3級，淮北、六安和池州均屬于1級；由減法集對勢法的評價結(jié)果可知4市2018年水生態(tài)文明建設的發(fā)展趨勢依次處于2級、1級、1級、1級狀態(tài)。這4種方法得出的結(jié)果具有較高的契合度。

b.評價年的聯(lián)系數(shù)分量評價結(jié)果刻畫了4市2018年水生態(tài)文明建設處于5個評價等級的分布信息結(jié)構(gòu)，它所提供的各評價等級信息結(jié)構(gòu)明顯比只有一維實數(shù)值的級別特征值法、屬性識別法、加權(quán)綜合得分法、減法集對勢的評價結(jié)果豐富；級別特征值法、加權(quán)綜合得分法的評價結(jié)果比屬性識別法的評價結(jié)果更精細，屬性識別法的評價結(jié)果含有置信概率信息，它的評價結(jié)論比較可靠、穩(wěn)妥。綜合考慮這4種方法得出的評價結(jié)果更為穩(wěn)妥，4市2018年水生態(tài)文明建設的綜合評價結(jié)果分別處于2級、1級、1級、1級狀態(tài)。

c.由表1中的評分可知，馬鞍山市需改進的指標有規(guī)劃編報情況、城市排澇達標率、水土流失治理率、城市污水處理回用率和城市水功能區(qū)水質(zhì)達標率；對淮北市需改進的指標有水利工程管理到位率、城市防洪達標率、自然水景觀、生態(tài)需水維持度、水域及周邊景點觀賞性及水文化特色；六安市在萬元GDP用水量上需要改善；池州市城市污水處理回用率和跡地恢復治理率需要提高。

d.盡管池州市的級別特征值法、加權(quán)綜合得分法的評價結(jié)果優(yōu)于淮北市，但池州市的減法集對勢值劣于淮北市，這是由于池州市的跡地恢復治理率、城市污水處理回用率明顯處于不和諧狀態(tài)，存在明顯的短板指標，這是池州市亟需改進之處。4市中六安市水生態(tài)文明建設的發(fā)展趨勢相對最佳。

e.由表2評價結(jié)果與安徽省水資源公報、安徽省城市水生態(tài)相關(guān)資料對比分析可知，本文評價結(jié)果與4市水生態(tài)文明建設總體現(xiàn)狀有較高的吻合度，表明了本文方法的合理性。4市水生態(tài)文明建設現(xiàn)狀評價結(jié)果區(qū)分度明顯，評價過程穩(wěn)健，診斷分析合理，適用性廣，表明本文方法在其他城市水生態(tài)文明建設動態(tài)評價和診斷分析中具有推廣應用價值。

f.圖1評價指標體系與文獻[4]一致，但在確定指標權(quán)重時，本文采用的AGA-FAHP能有效克服文獻[4]中傳統(tǒng)層次分析法構(gòu)造、檢驗和修正判斷矩陣的一致性問題，計算結(jié)果更穩(wěn)定，通用性也更高；在進行水生態(tài)文明建設評價時，文獻[9]分別采用打分法和歐式距離法對10座城市進行評價，但兩種方法結(jié)果相差較大，而本文采用的聯(lián)系數(shù)方法得出的結(jié)果精確穩(wěn)定，說明本文采用的評價方法更具穩(wěn)定性。

4 結(jié) 語

本文將城市水生態(tài)文明建設評價與所在城市經(jīng)濟社會發(fā)展、日常生產(chǎn)生活結(jié)合，為了能篩選出適用當前城市水生態(tài)文明建設的評價指標，定量準確評價城市水生態(tài)文明建設現(xiàn)狀和發(fā)展趨勢，在運用AGA-FAHP篩選評價指標的基礎上，提出了基于聯(lián)系數(shù)的城市水生態(tài)文明建設的評價方法，并以安徽省4個代表城市為例進行了驗證，結(jié)果顯示，該方法合理可靠，評價結(jié)果與實際情況契合度高，較好地解決了評價指標與評價等級之間不確定性方面的問題。本文提出的方法評價過程穩(wěn)健，診斷分析合理，適用性廣，在水生態(tài)相關(guān)問題的動態(tài)評價和診斷分析中具有推廣應用價值。