李 非

(遼寧省水文局，遼寧 沈陽 110003)

流域防洪安全是對流域防洪體系中可能存在的危險和不安全程度進行評估，流域防洪安全的準確評估對于區(qū)域防洪安全至關重要，為此國內許多學者開展流域防洪安全評估的研究[1- 5]，但當前對防洪安全大都采用層次分析方法進行評估研究，而流域防洪安全影響要素較多，防洪安全評估存在模糊和隨機性的特點，傳統(tǒng)層次分析方法在區(qū)域防洪安全評估存在一定的局限性。近些年來，云模型在數(shù)據(jù)挖掘的基礎上結合數(shù)據(jù)模糊和隨機性的特點，降低了評估指標的不確定性，在許多領域的指標評估方面得到具體應用[6- 9]，考慮到不同指標的權重，為對指標權重進行客觀評價，本文將云模型和熵權方法進行組合，并將組合模型應用到中小河流防洪安全的評估研究，從而為中小河流防洪體系的架構提高科學的參考依據(jù)。

1 熵權云模型組合原理

(1)

式中，μ(x)—模型論域U上的正態(tài)分布云函數(shù)；EX—云模型的期望值；En—云模型的熵權值。

云模型在具體計算時，首先生成具有正態(tài)分布的隨機數(shù)，計算方程為：

En′～N(En，He2)

(2)

式中，En—隨機數(shù)組的期望；He2—隨機數(shù)組的方程值，并在此基礎上進一步在計算正態(tài)隨機數(shù)組，計算方程為：

x～N(Ex，En2)

(3)

式中，Ex—正態(tài)隨機數(shù)組序列的期望，En—正態(tài)隨機數(shù)組序列的方差。

在兩次正態(tài)隨機數(shù)組計算的基礎上，利用特征值x和期望值Ex代入到式(1)計算云模型的指標隸屬度。為客觀定量確定評價指標的權重，本文結合熵權方法對各個評價指標進行權重的確定，熵權方法的計算步驟為：

(1)假定m個評價對象，每個對象包含n個評價指標，并構建評價判定矩陣，該判定矩陣的表達式為：

R=(rij)m×n

(4)

式中，R—判定矩陣；r—判定指數(shù)；i—評價對象；j—評價指標。在判定矩陣確定的基礎上，確定各個指標的熵值，計算方程為：

(5)

式中，Hj—各評價指標的熵權值；fij—指標的熵值參數(shù)，其計算方程為：

(6)

式中，rij、m的變量含義同式(4)中變量含義。在各指標熵值確定的基礎上，需要對各評價指標的熵權值進行確定，計算方程為：

(7)

式中，Hj—各評價指標的熵值。

2 熵權云模型組合計算步驟

熵權云組合模型將熵權和云模型進行組合，模型的主要建立步驟為：

(1)建立評估對象各個指標的集合矩陣u={u1，u2，…un}，評價矩陣為v={v1，v2，…vm}；

(2)結合熵權方法對各評價指標權重進行計算，確定權重矩陣W={w1，w2，…wn}；

(3)結合指標和評價權重進行單要素評價，建立各指標模糊矩陣R，矩陣R對應的不同評價指標的不同等級下的隸屬度，本文采用云模型計算不同指標的隸屬度，其中期望值Ex和En的計算表達式為：

(8)

(9)

結合各要素指標評價值，確定不同評價指標的不同等級的云模型隸屬度矩陣Z=(zij)n×m，采用云模型進行正向運算，計算得到不同等級隸屬度下的指標評估綜合值Zij，計算方程為：

(10)

式中，N—運算次數(shù)；K—評價指標數(shù)。

3 實例應用

3.1 研究區(qū)域概況

本文以遼寧中部某中小河流為研究實例，該河流的防洪標準為20年一遇，為對該河流的防洪安全進行定量評價，選取12個指標進行評價指標體系的構建，各指標的屬性值見表1。此外本文為對區(qū)域不同分區(qū)的防洪安全進行定量評價，將區(qū)域劃分成6個分區(qū)，結合熵權云模型組合方法對各分區(qū)的防洪安全進行定量評估。

表1 區(qū)域各指標的屬性值

3.2 防洪安全評估指標熵權值

結合熵權分析方法對區(qū)域不同頻率洪水下各防洪安全評價指標進行熵權值的計算，計算結果見表2、3。

表2 防洪安全評定指標評估熵權值計算結果(10年一遇洪水)

表3 防洪安全評定指標評估熵權值計算(20年一遇洪水)

從表2可以看出，在10年一遇洪水下，人均GDP不同頻率洪水下的熵權值在0.11～0.29之間，防汛指揮體系下不同頻率洪水的熵權值在0.42～0.59之間，可見人均GDP對區(qū)域防洪安全影響較小，而完善的防汛指揮體系對區(qū)域中小河流的防洪安全至關重要。熵權云模型可對各個指標進行熵權計算，實現(xiàn)對各指標熵權的客觀賦值。

3.3 評價指標的正態(tài)云標準矩陣

結合各防洪安全評價指標熵權值和屬性值，結合式(9)和(10)計算了不同頻率洪水下的正態(tài)云模型標準值，建立各評價指標的正態(tài)云標準矩陣，計算結果見表4、5。

表4 各評價指標正態(tài)云計算標準值計算結果(10年一遇)

表5 各評價指標正態(tài)云計算標準值結果(20年一遇)

3.4 各評價指標的隸屬度矩陣

結合式(1)和各指標正態(tài)云標準值，假定N=1000，結合對應的指標正態(tài)云模型標準值。結合正向標準算法計算了不同指標的隸屬度矩陣，各頻率洪水的隸屬度矩陣計算結果見表6、7。

表6 各評價指標正態(tài)云隸屬度計算結果(10年一遇)

表7 各評價指標正態(tài)云隸屬度計算結果(20年一遇)

從表中可以看出，各個評價指標的正態(tài)云隸屬度函數(shù)呈現(xiàn)標準的正向正態(tài)分布，這主要是因為熵權云模型結合正向標準正態(tài)函數(shù)對各個指標的隸屬度進行計算，綜合考慮了各評價指標的隨機性和模糊性的特點，并可對數(shù)據(jù)系列進行多層次的數(shù)據(jù)挖掘，實現(xiàn)不同層次指標的綜合定量評價。

3.5 基于熵權云模型的組合的防洪安全評價結果

在各指標隸屬度確定的基礎上，結合熵權云模型組合方法對各分區(qū)不同頻率下的洪水進行防洪安全綜合評估，綜合評估結果見表8、9。

表8 區(qū)域防洪安全綜合評價結果(10年一遇)

表9 區(qū)域防洪安全綜合評價結果(20年一遇)

從表8可以看出，當發(fā)生10年一遇洪水時，各分區(qū)綜合評價值在0.1622～0.2921，經綜合評價區(qū)域防洪形勢總體達到較安全程度，這和區(qū)域實際情況較為吻合，區(qū)域中小河流發(fā)生10年一遇標準洪水時，由于區(qū)域具有一定的防洪體系和防洪工程措施，因此在發(fā)生10年一遇標準洪水時，區(qū)域防洪還較為安全，而當發(fā)生20年一遇標準洪水時，各分區(qū)防洪安全綜合評估至在0.0979～0.2407之間，區(qū)域防洪總體為較不安全的程度，這主要是因為區(qū)域中小河流的防洪體系可抵御的洪水標準還未達到20年一遇，區(qū)域的防洪體系、防洪工程、防洪非工程措施力度還需要進一步加大。

4 結語

本文將熵權與云模型進行組合，并將組合模型應用于中小河流防洪安全的定量綜合評估研究中，研究取得下述結論。

(1)熵權云模型組合方法可綜合考慮評價指標的隨機性和模糊性，并可對數(shù)據(jù)系列進行多層次的數(shù)據(jù)挖掘，實現(xiàn)不同層次指標的綜合定量評價，可客觀合理對各評價對象進行綜合評估，該方法可為中下河流和山洪災害防洪安全綜合評估提供一種新的思路和方法。

(2)該方法可以分區(qū)為計算單元，對各分區(qū)的防洪安全進行綜合評估，可為區(qū)域防洪措施規(guī)劃提供決策依據(jù)，對于綜合評估較不安全的區(qū)域，應加強防洪工程和非工程措施力度。

[1] 王欣, 高青峰. 運糧河流域防洪安全等級評價研究[J]. 水利科技與經濟, 2015(12): 90- 91.

[2] 趙洪杰, 唐德善. 流域防洪體系防洪安全評價研究[J]. 災害學, 2006(04): 31- 35.

[3] 余江順. 大連市防洪安全策略研究[J]. 水利規(guī)劃與設計, 2015(09): 26- 27+123.

[4] 周健. 中小河流治理規(guī)劃設計基本思路與要點[J]. 水利規(guī)劃與設計, 2012(05): 6- 7+52.

[5] 吳新新. 希尼爾水庫防洪安全復核分析[J]. 水利技術監(jiān)督, 2016(03): 61- 63.

[6] 楊理智, 張韌, 洪梅, 等. 基于云模型的我國西南邊境水資源安全風險評估[J]. 長江流域資源與環(huán)境, 2014(S1): 1- 5.

[7] 江新, 朱沛文, 沈力. 基于ANP和云模型的水電項目群資源沖突風險評估[J]. 中國安全科學學報, 2014(11): 152- 158.

[8] 宮凡荔. 基于界殼論和云模型的區(qū)域水資源可持續(xù)利用方案分析與風險評估[D]. 東北農業(yè)大學, 2014.

[9] 江迎. 基于云模型和GIS/RS的壩堤潰決風險分析及災害損失評估研究[D]. 華中科技大學, 2012.