羅日洪，黃錦林

(1.廣東省水利水電科學(xué)研究院，廣東 廣州 510635;2. 廣東省粵港澳大灣區(qū)水安全保障工程技術(shù)研究中心，廣東 廣州 510635)

1 概述

截至2018年底，我國已建成各類水庫98 822座，水庫總庫容8 953億m3[1]。這些水庫的擋水建筑物絕大部分為土石壩[2]。土石壩在建設(shè)和運用過程中由于各類不利因素組合，容易導(dǎo)致工程質(zhì)量變差，進而演變?yōu)椴‰U工程。因此，為掌握土石壩安全狀態(tài)和加強管理而進行定期安全評價具有重要意義。

土石壩的安全狀態(tài)是一個受多種因素影響的復(fù)雜系統(tǒng)。傳統(tǒng)的土石壩安全評價大多以定性的方式(規(guī)范法)，從偏安全的角度考慮，一般以評為最差的單項評價作為最終評價結(jié)果，理論簡單、直接，但缺乏綜合考量和定量結(jié)果,決策時各部分的重要性沒有得到很好的突出。目前，國內(nèi)外很多學(xué)者針對堤壩工程安全評價進行了研究。Chauha和Bowles[3]提出了大壩安全風(fēng)險評價中的不確定性分析框架，用與滿足可容忍風(fēng)險準(zhǔn)則相關(guān)的置信度表示風(fēng)險分析結(jié)果中的不確定性；Vuillet等[4]基于多準(zhǔn)則聚合的確定性模型，采用蒙特卡羅模擬建立每個堤段的概率分布；Park等[5]提出了一種能同時考慮堤壩漫水、入滲、侵蝕等多種災(zāi)害因素的洪水風(fēng)險指數(shù)，并對堤防安全進行了可靠的評價；Coelho等[6]基于材料點法為重新定義了堤壩抗剪破壞安全系數(shù)，并對堤防的漸進性破壞及其行為進行了分析；胡建平等[7]根據(jù)突變理論，計算了土石壩安全級別評價指標(biāo)值；蔡新等[8]引進灰色理論構(gòu)建堤壩安全風(fēng)險評價指標(biāo)體系及堤壩安全評價數(shù)學(xué)模型。以上研究中，有些僅是較為單一因素的研究，未綜合考慮各種因素的影響，有些對原始數(shù)據(jù)要求較高，并且在擬合度方面存在局限。因此，本文提出一種堤防安全綜合評價方法，引入改進AHP法、熵權(quán)法進行組合權(quán)重優(yōu)化計算，基于熵權(quán)的專家賦分群組決策方法，對土石壩安全進行綜合評價，并結(jié)合具體工程，確定出較為合理的評價結(jié)果，為后續(xù)采取措施提供科學(xué)依據(jù)。

2 基本原理

2.1 土石壩安全評價體系建立

在分析土石壩結(jié)構(gòu)特征和破壞機理的基礎(chǔ)上，構(gòu)建土石壩安全評價體系。為了使評價方法具有可操作性性，依據(jù)SL 258—2017《水庫大壩安全評價導(dǎo)則》對土石壩安全評價的規(guī)定[9]，并結(jié)合實際工程出險分析[10]，構(gòu)建由目標(biāo)層、準(zhǔn)則層、指標(biāo)層組成的安全評價指標(biāo)體系，指標(biāo)層共計26個評價指標(biāo)(見圖1)。

圖1 土石壩安全綜合評價指標(biāo)體系示意

2.2 改進AHP-熵權(quán)法組合優(yōu)化賦權(quán)

AHP法是一種系統(tǒng)性分析方法，在各種決策問題中具有廣泛的應(yīng)用，具有簡單明了、邏輯清晰等特點，在各種系統(tǒng)決策問題中應(yīng)用廣泛。但用AHP方法決策易受個人偏好影響，主管性較強，甚至出現(xiàn)判斷失真或誤判。為了解決AHP法權(quán)重的不確定性的問題，通過組建不同專業(yè)的專家組成群組決策，構(gòu)建專家自身客觀權(quán)重熵模型，用熵表示專家和各專家與理想專家的水平差異，反映了各專家權(quán)重的貢獻度，得到專家自身客觀權(quán)重。利用專家自身客觀權(quán)重對AHP法權(quán)重進行修正，目的是提高權(quán)重的合理性和精確性。

2.2.1基于改進AHP法的主觀權(quán)重

AHP法由Saaty提出[11]，主要應(yīng)用于系統(tǒng)問題決策。AHP法采用的是1～9的標(biāo)度，有研究認(rèn)為 1～9標(biāo)度中各等級之間關(guān)系存在不合理性[12-13]。鑒于此，有學(xué)者提出指數(shù)an,a8標(biāo)度[14]、9n/9標(biāo)度[13]、0.1～0.9標(biāo)度[15]?；谧匀恢笖?shù)e0/5～e8/5的標(biāo)度在保序性、一致性、標(biāo)度均勻性、標(biāo)度可記憶性、標(biāo)度可感知性、標(biāo)度權(quán)重擬合性具有優(yōu)勢，建議用于精度要求較高的多準(zhǔn)則排序問題(見表1)[16]。

表1 自然指數(shù)標(biāo)度和1～9標(biāo)度對比

本文采用自然指數(shù)標(biāo)度計算權(quán)重，步驟如下：

1)假設(shè)有m個專家，根據(jù)專家k意見，將評價對象同層元素X1,X2,…,Xn兩兩比較，采用e0/5～e8/5標(biāo)度進行賦值，構(gòu)造出符合下列要求的判斷矩陣A=(aijk)n×n：

(1)

式中aijk為第i個因素相對于第j個因素的重要程度的兩兩比較標(biāo)度。

2)首先歸一化處理判斷矩陣，將得到的矩陣按行相加得特征向量。對特征向量歸一化處理，即得專家k的各個指標(biāo)主觀權(quán)重：

(2)

其中：

式中Sijk為元素aijk歸一化結(jié)果；wik為歸一化后的元素Sijk按行相加的結(jié)果；Wik為歸一化后的主觀權(quán)重，則專家k的主觀權(quán)重向量為Wk=(W1k,W2k,…,Wnk)T。

3)分別對專家組中m個專家意見確定指標(biāo)的主觀權(quán)重向量，形成主觀權(quán)重矩陣如下：

(3)

2.2.2基于熵權(quán)法的專家自身客觀權(quán)重

建立熵模型步驟如下：

2)確定專家主觀權(quán)重水平向量。各專家主觀權(quán)重結(jié)果與最優(yōu)專家S*主觀權(quán)重結(jié)果進行比較，以差異的大小來衡量主觀權(quán)重質(zhì)量的好壞。第k個專家的主觀權(quán)重水平向量Ek為：

Ek=(ek1,ek2,…,ekn)k=1,2,…,m

(4)

其中：

3)建立專家自身客觀權(quán)重的熵模型：

(5)

其中：

式中Hk為專家k主觀權(quán)重的不確定度；hik為專家k對第i個指標(biāo)的主觀權(quán)重熵值。

4)確定專家自身客觀權(quán)重。專家評價結(jié)果的可信度用不確定度Hk表示，Hk越大說明可信度越低。專家客觀權(quán)重Ck可表示為：

(6)

式中Ck越大，表示專家k的意見在評價結(jié)果中比重越大。

2.2.3改進AHP-熵權(quán)法組合優(yōu)化賦權(quán)

采用改進AHP法計算的主觀權(quán)重矩陣為W，采用熵權(quán)法計算的專家權(quán)重向量為C=(C1,C2,…,Cm)T，基于加權(quán)融合獲得的組合優(yōu)化權(quán)重為：

(7)

2.3 土石壩安全綜合評價

2.3.1安全分類和評價標(biāo)準(zhǔn)

根據(jù)SL 258—2017，大壩安全分類劃分為3類。大壩安全的綜合評價是在現(xiàn)狀安全檢查和監(jiān)測資料的基礎(chǔ)上，根據(jù)防洪能力、滲流安全、結(jié)構(gòu)安全、抗震安全、金屬結(jié)構(gòu)安全等專項復(fù)核成果，并參考工程質(zhì)量與大壩運行管理評價結(jié)論，對大壩安全進行綜合評價，評定其類別。鑒于SL 258—2017的方法應(yīng)用廣泛、方法成熟、可操作性強，本文在借鑒SL 258—2017對大壩安全分類的基礎(chǔ)上，制定評分標(biāo)準(zhǔn)(見表2)。底層指標(biāo)的安全級別也與表2相對應(yīng)，劃分為{A，B，C}3級，分別對應(yīng){安全，基本安全，不安全}。

2.3.2專家評分

邀請水工、地質(zhì)、水文、金屬結(jié)構(gòu)等專業(yè)的專家組成專家組(專家組專家數(shù)量應(yīng)為單數(shù))到現(xiàn)場進行檢查，并在查閱大壩現(xiàn)狀安全檢查報告和監(jiān)測資料分析報告，以及各分項安全復(fù)核評價報告的基礎(chǔ)上，獨立對2.1節(jié)所列的大壩評價指標(biāo)的安全狀況進行評分。在綜合各個專家評分的方法上可以選擇簡單的算術(shù)平均法。算術(shù)平均法默認(rèn)所有專家評分的權(quán)重都相等，鑒于專家不可能對所有專業(yè)領(lǐng)域都能熟練掌握，不同專家對同一個評價對象的評價結(jié)果是存在重要性差異的。所以仍然采用前述的熵權(quán)法，假定存在一個理想的最優(yōu)專家，對專家評分結(jié)果進行權(quán)重結(jié)算，在與對應(yīng)的評分相乘，得到綜合評分，使評價結(jié)果更為合理。

2.3.3安全類別評估計算

現(xiàn)行規(guī)范SL 258—2017偏安全地以防洪能力、結(jié)構(gòu)安全、滲流安全、抗震安全、金屬結(jié)構(gòu)安全等單項評價的最低等級來確定大壩的安全類別，某一單項評價對最終結(jié)果具有決定權(quán)，造成評定結(jié)果與實際情況存在較大出入。而根據(jù)各類不同因素對大壩安全影響程度的不同，賦予不同權(quán)重，通過一定手段把不同因素的多源信息進行融合，對大壩進行綜合評價的方法，更加科學(xué)、合理，得到的結(jié)果更加符合大壩的實際安全狀態(tài)。

根據(jù)2.3.2節(jié)的方法，底層指標(biāo)綜合得分與相應(yīng)權(quán)重相乘，得到準(zhǔn)則層的得分：

(8)

式中D為準(zhǔn)則層的得分值；Wi為第i個底層指標(biāo)的權(quán)重；Ci為第i個底層指標(biāo)的得分值。

目標(biāo)層評分的計算則采用準(zhǔn)則層各指標(biāo)加權(quán)得到。

(9)

式中F為目標(biāo)層評分值；Wj為第j個準(zhǔn)則層指標(biāo)權(quán)重。

表2 大壩安全評定分類

3 工程應(yīng)用

3.1 工程概況

銀盞水庫位于清遠市清城區(qū)東南20 km，位于北江二級支流銀盞河上游，為中型水利工程。1959年11月動工興建，分四期續(xù)建，1987年完成工程，2000年進行除險加固，水庫集雨面積為28.4 km2，總庫容為 2 547.3萬m3。正常水位為58.00 m，相應(yīng)庫容為 2 333萬m3。工程等別為Ⅲ等，防洪標(biāo)準(zhǔn)按50年一遇洪水設(shè)計，2000年一遇洪水校核。銀盞水庫工程由大壩、溢洪道、輸水涵管、壩后電站等建筑物組成。水庫大壩為均質(zhì)土壩，大壩長為280 m，壩頂寬為 6.5 m，壩頂高程為64.67 m，最大壩高為41.67 m。上游坡為干砌石護坡；下游坡為草皮護坡。溢洪道位于大壩右端，由兩扇弧形鋼閘門控制。輸水涵管為壩內(nèi)填埋式鋼筋混凝土壓力管，管道內(nèi)徑為1.2 m。經(jīng)過多年運行，工程存在的主要問題為壩體填土滲透系數(shù)較大，不滿足規(guī)范要求；輸水涵管下穿壩體，混凝土老化嚴(yán)重，并存在向壩體滲漏的現(xiàn)象。這些問題已嚴(yán)重影響大壩的安全運行。

3.2 指標(biāo)權(quán)重

以2.1節(jié)土石壩評價層次體系為基礎(chǔ)，采用 2.2節(jié)改進AHP和熵權(quán)組合優(yōu)化權(quán)重計算方法，計算各指標(biāo)層和準(zhǔn)則層各要素的權(quán)重。根據(jù)專家組對指標(biāo)采用改進AHP法兩兩比較結(jié)果，按式(1)建立判斷矩陣，然后按(2)式進行和法近似計算得到主觀權(quán)重，經(jīng)計算，指標(biāo)層C對準(zhǔn)則層B、準(zhǔn)則層B對目標(biāo)層A均滿足一致性檢驗要求；根據(jù)得到的主觀權(quán)重，按(4)式計算主觀權(quán)重水平向量，進而按(5)式計算得到主觀權(quán)重的不確定度，根據(jù)(6)式確定專家自身客觀權(quán)重。最后根據(jù)式(7)將AHP得到的專家主觀權(quán)重和熵權(quán)法計算的專家意見客觀權(quán)重融合計算，得到優(yōu)化后的權(quán)重結(jié)果見表3。

在土石壩安全評價中，準(zhǔn)則層指標(biāo)權(quán)重大小依次是結(jié)構(gòu)安全、防洪能力、滲流安全、金屬結(jié)構(gòu)安全、工程質(zhì)量、運行管理。

表3 指標(biāo)權(quán)重結(jié)果

3.3 綜合評價

專家組在充分了解工程狀態(tài)的基礎(chǔ)上，結(jié)合現(xiàn)場檢查、質(zhì)量檢測和單項安全復(fù)核結(jié)果[17]，以及自身豐富的工程評價經(jīng)驗，對底層指標(biāo)進行賦分，以群決策方法計算綜合評分，更能代表工程安全的實際狀態(tài)。基于熵權(quán)法和算術(shù)平均法計算的綜合評分結(jié)果見表4。

以表4得到的各個指標(biāo)的綜合評分與表3對應(yīng)的權(quán)重，按照式(8)可計算得到準(zhǔn)則層綜合評分，結(jié)果見表4和圖2所示，鑒于該工程檢測結(jié)果中壩體填土滲透系數(shù)大于規(guī)范規(guī)定，輸水涵管直穿壩體，老化嚴(yán)重、存在滲漏現(xiàn)象且?guī)Р」ぷ?，根?jù)專家群組決策，準(zhǔn)則層工程質(zhì)量評價B1綜合得分按熵權(quán)法計算結(jié)果為59.69，為C級，質(zhì)量不合格；若按算術(shù)平均法計算結(jié)果為60.17，為B級，則為合格；滲流安全評價B5得分按熵權(quán)法計算結(jié)果為59.94，屬于C級；若按算術(shù)平均法計算結(jié)果為61.24，則為B級。

表4 各指標(biāo)賦分

圖2 準(zhǔn)則層綜合得分情況示意

在得到準(zhǔn)則層綜合評分的基礎(chǔ)上按式(9)計算得到總目標(biāo)土石壩安全狀態(tài)的總得分，為71.64分。但根據(jù)表2的評價標(biāo)準(zhǔn)，按熵權(quán)法土石壩工程質(zhì)量不合格，有一項評為C，則土石壩綜合評為3類壩，若按算術(shù)平均法，則評價結(jié)果為2類壩。根據(jù)參考文獻[17]可知，該工程安全鑒定結(jié)果也為3類壩，熵權(quán)法評價結(jié)果和文獻[17]結(jié)果對比見表5。

表5 土石壩安全綜合評價結(jié)果對比

按照規(guī)范方法進行的安全評價，由于是定性評價，無法確定工程及各單項所處的類別在該評級中的水平。按照表5本文對各準(zhǔn)則層的評價結(jié)果，結(jié)合表4中的各自的綜合評分和工程實際安全狀態(tài)，工程質(zhì)量、滲流安全均為B級，但這兩項的評分與B差距較小，主要是因為壩體填土滲流系數(shù)達不到規(guī)范要求，以及壩體中輸水涵管的老化滲漏。后續(xù)應(yīng)進行除險加固，通過壩體灌漿、改建輸水管道等工程措施，消除隱患，使該項提升到B級標(biāo)準(zhǔn)。因此，本例中通過群組決策得到專家對工程各單項評價結(jié)果的綜合評分，能更準(zhǔn)確地評判土石壩的安全狀態(tài)，對工程后續(xù)運行管理和風(fēng)險防范措施制定的指導(dǎo)性更強。

4 結(jié)語

本文建立了與現(xiàn)行土石壩安全評價規(guī)范相一致的土石壩安全狀態(tài)綜合評價體系，綜合考慮了影響土石壩安全的要素，融合得到工程質(zhì)量、運行管理、防洪能力、結(jié)構(gòu)安全、滲流安全、抗震安全、金屬結(jié)構(gòu)安全七大方面因素的影響。在此基礎(chǔ)上，引入自然指數(shù)標(biāo)度對傳統(tǒng)AHP方法進行改進，得到評價指標(biāo)的主觀權(quán)重，然后利用熵權(quán)法計算指標(biāo)客觀權(quán)重，并進行組合優(yōu)化，使得各指標(biāo)的權(quán)重分配更加科學(xué)、合理；最后以現(xiàn)場檢查、檢測報告和各項安全評價復(fù)核結(jié)果為基礎(chǔ)，進行專家打分。同時采用算術(shù)平均法和熵權(quán)法獲取專家組群決策的綜合評分，參考現(xiàn)有土石壩安全評價規(guī)范的評級標(biāo)準(zhǔn)，得到定量化安全狀態(tài)結(jié)果，確定了土石壩的安全類別，可以直觀地反映出土石壩的實際安全狀況，為土石壩采取必要的風(fēng)險防范措施提供科學(xué)支撐。