李曉英,曾宇軒
(1.河海大學(xué) 水利水電學(xué)院,江蘇 南京 210098;2.河海大學(xué) 水利工程實(shí)驗(yàn)教學(xué)中心,江蘇 南京 210098)
安全是可持續(xù)發(fā)展的基礎(chǔ)。水利工程作為涉及能源、環(huán)境、社會(huì)穩(wěn)定等諸多方面的復(fù)雜系統(tǒng)工程,其投資大、不確定因素多,一旦出現(xiàn)安全事故將造成巨大的經(jīng)濟(jì)損失,甚至?xí)?duì)社會(huì)公共安全產(chǎn)生巨大的威脅[1];同時(shí),水利行業(yè)屬于資金密集型產(chǎn)業(yè),如何在其建設(shè)和運(yùn)營階段進(jìn)行有效地管理對(duì)水利事業(yè)的發(fā)展具有重要參考價(jià)值[2]。雖然水利工程具有高風(fēng)險(xiǎn)性,但我國水電企業(yè)的經(jīng)營者和管理者風(fēng)險(xiǎn)意識(shí)普遍不強(qiáng),缺乏系統(tǒng)和明確的風(fēng)險(xiǎn)管理目標(biāo)[3]。因此,構(gòu)建現(xiàn)代化水電企業(yè)風(fēng)險(xiǎn)管理體系具有重要的意義。
風(fēng)險(xiǎn)評(píng)估能夠用具體的數(shù)值直觀反映風(fēng)險(xiǎn)的大小,是構(gòu)建系統(tǒng)風(fēng)險(xiǎn)管理體系的基礎(chǔ),相關(guān)方法被廣泛應(yīng)用于安全管理領(lǐng)域[4-5]。風(fēng)險(xiǎn)指標(biāo)的權(quán)重計(jì)算是風(fēng)險(xiǎn)評(píng)估的核心所在,其計(jì)算方法主要可分為主觀賦權(quán)法和客觀賦權(quán)法兩類,在主觀賦權(quán)法方面,現(xiàn)階段已有不少相關(guān)研究,李創(chuàng)[6]運(yùn)用風(fēng)險(xiǎn)評(píng)估中傳統(tǒng)的AHP-模糊綜合評(píng)價(jià)法對(duì)納米比亞歐曼德爾管道工程項(xiàng)目施工中的風(fēng)險(xiǎn)因素進(jìn)行了綜合評(píng)價(jià),由于只考慮風(fēng)險(xiǎn)嚴(yán)重程度一個(gè)輸入變量,其主觀性和隨機(jī)性較強(qiáng);在此基礎(chǔ)上,林倩[7]利用風(fēng)險(xiǎn)矩陣綜合考慮風(fēng)險(xiǎn)發(fā)生概率和后果嚴(yán)重程度兩個(gè)輸入變量來劃分風(fēng)險(xiǎn)等級(jí),并引入?yún)^(qū)間數(shù)削弱主觀賦權(quán)的隨機(jī)性,但其通過專家賦分選取區(qū)間數(shù),仍有較強(qiáng)的主觀性。近年來,部分學(xué)者將智能算法引入主觀賦權(quán),取得了較好的效果。江新等[8]利用集對(duì)理論和云模型削弱傳統(tǒng)AHP(analytic hierarchy process)法的模糊性和不確定性,建立基于SPA(set pair analysis)-云模型的工程運(yùn)營風(fēng)險(xiǎn)評(píng)估模型,該模型能有效處理專家判斷的同一性和差異性,減少主觀誤差;程江洲等[9]結(jié)合貝葉斯網(wǎng)絡(luò)和Dirichlet模型分解水力發(fā)電系統(tǒng)以提高主觀賦權(quán)的精度,并在此基礎(chǔ)上運(yùn)用權(quán)威度賦權(quán)和模糊隸屬函數(shù)進(jìn)一步削弱主觀賦權(quán)的隨機(jī)性;此外,傅鶴林等[10]引入語言評(píng)價(jià)值來反映指標(biāo)的重要程度,該方法能較好地反映系統(tǒng)的模糊性、不確定性和復(fù)雜性,進(jìn)而削弱賦分的隨機(jī)性。以上研究嘗試多種方法削弱主觀賦權(quán)法的主觀性與隨機(jī)性,取得了一定的成效,但水利工程風(fēng)險(xiǎn)復(fù)雜性高,專家難以給出確切評(píng)分,主觀賦權(quán)的隨機(jī)性大且很難規(guī)避。在客觀賦權(quán)法方面,沈菊琴等[11]采用熵權(quán)法計(jì)算指標(biāo)權(quán)重以削弱主觀性,形成基于熵權(quán)可拓物元模型的水利工程穩(wěn)定風(fēng)險(xiǎn)評(píng)估方法;盧丹等[12]運(yùn)用反熵權(quán)法對(duì)電網(wǎng)風(fēng)險(xiǎn)進(jìn)行客觀賦權(quán),并運(yùn)用博弈理論結(jié)合相關(guān)主觀賦權(quán)結(jié)果綜合兩種賦權(quán)方法求解最終權(quán)重。但目前多數(shù)客觀賦權(quán)法的基礎(chǔ)數(shù)據(jù)仍源于專家評(píng)分或問卷調(diào)查,不能直接接觸實(shí)際的客觀數(shù)據(jù)。
由于傳統(tǒng)風(fēng)險(xiǎn)評(píng)估中賦權(quán)的主觀性和隨機(jī)性強(qiáng),且多為基于文字和經(jīng)驗(yàn)意義上的風(fēng)險(xiǎn)描述,不能直接接觸實(shí)際數(shù)據(jù)[13-15],需尋找一種基于實(shí)際客觀數(shù)據(jù)的評(píng)分法以代替主觀評(píng)分??紤]到通過保險(xiǎn)理賠數(shù)據(jù)能較精確地識(shí)別和定量風(fēng)險(xiǎn)[16],同時(shí)近年來國內(nèi)水電企業(yè)實(shí)現(xiàn)全國運(yùn)營項(xiàng)目保險(xiǎn)統(tǒng)一安排,保險(xiǎn)的標(biāo)準(zhǔn)化、規(guī)范化情況較好[17],且出險(xiǎn)資料實(shí)時(shí)更新,具有實(shí)用價(jià)值,故其可以較好地被選作代替主觀評(píng)分的客觀數(shù)據(jù)。本研究以水電企業(yè)的保險(xiǎn)理賠數(shù)據(jù)為切入點(diǎn)展開分析,提出了一種基于保險(xiǎn)理賠數(shù)據(jù)的水電企業(yè)風(fēng)險(xiǎn)評(píng)估方法。
本文研究對(duì)象為某水電企業(yè),收集到的數(shù)據(jù)資料包括其運(yùn)營的分布于新疆、湖南、云南、寧夏、四川、河南、江西7個(gè)省份的24座水電站2016-2019年的所有保險(xiǎn)理賠數(shù)據(jù),數(shù)據(jù)質(zhì)量良好。
2.2.1 風(fēng)險(xiǎn)因素指標(biāo)體系構(gòu)建 通過分析某水電企業(yè)的保險(xiǎn)理賠數(shù)據(jù)并結(jié)合相關(guān)專著[1],總結(jié)近年來水電企業(yè)存在的主要運(yùn)營風(fēng)險(xiǎn),按成因?qū)⑵鋭澐譃槿鐖D1所示的3個(gè)類別和14個(gè)因素。
圖1 水電企業(yè)運(yùn)營風(fēng)險(xiǎn)因素指標(biāo)體系
2.2.2 風(fēng)險(xiǎn)損失區(qū)間及其標(biāo)度 目前水電企業(yè)的保險(xiǎn)理賠數(shù)據(jù)較少,難以精確定量風(fēng)險(xiǎn)損失,故通過區(qū)間數(shù)來提高風(fēng)險(xiǎn)損失估計(jì)的準(zhǔn)確性。
(1)
K=
(2)
同一準(zhǔn)則下可能包含多項(xiàng)指標(biāo),為全面體現(xiàn)各項(xiàng)指標(biāo)避免遺漏,定義一種最優(yōu)風(fēng)險(xiǎn)損失區(qū)間(x-,x+)。設(shè)損失處于該區(qū)間內(nèi)的概率為P,則該區(qū)間需滿足下列條件:
①P≥95%
(n=1,2,…,∞)
(3)
(4)
(2)指標(biāo)層風(fēng)險(xiǎn)損失區(qū)間。指標(biāo)層通常涉及發(fā)生概率很低的風(fēng)險(xiǎn),其出險(xiǎn)數(shù)據(jù)資料少,不足以完成擬合。為保證指標(biāo)層風(fēng)險(xiǎn)損失區(qū)間在構(gòu)造上的一致性,采用百分位數(shù)求解指標(biāo)層風(fēng)險(xiǎn)損失區(qū)間。
設(shè)某理賠數(shù)據(jù)樣本的偏度系數(shù)為SK,第1四分位數(shù)為Q1,中位數(shù)為Q2,第3四分位數(shù)為Q3,則其風(fēng)險(xiǎn)損失區(qū)間可表示為:
(5)
(6)
其標(biāo)度C可表示為:
(7)
(8)
2.2.4 AHP法計(jì)算指標(biāo)權(quán)重 AHP法是將決策問題按目標(biāo)、準(zhǔn)則、指標(biāo)的順序分解為不同的層次結(jié)構(gòu),然后求解判斷矩陣特征向量進(jìn)而計(jì)算各指標(biāo)對(duì)總目標(biāo)的權(quán)重的方法。其具體步驟可歸納為建立層次結(jié)構(gòu)模型、構(gòu)造判斷矩陣、層次單排序及其一致性檢驗(yàn)、層次總排序及其一致性檢驗(yàn)。
(1)判斷矩陣。風(fēng)險(xiǎn)指數(shù)是風(fēng)險(xiǎn)的后果嚴(yán)重程度與發(fā)生概率的乘積,用以劃分風(fēng)險(xiǎn)等級(jí)[20],可利用風(fēng)險(xiǎn)指數(shù)作為判斷矩陣的計(jì)算依據(jù)。設(shè)某風(fēng)險(xiǎn)的發(fā)生概率為Pi,后果嚴(yán)重程度為Ci,風(fēng)險(xiǎn)指數(shù)為Ri,則Ri可按以下公式計(jì)算:
Ri=Pi·Ci
(9)
(10)
(2)權(quán)重向量。判斷矩陣A=(αij)n×m如滿足αij=αik·αkj,則稱其為一致性矩陣。一致性矩陣A存在唯一的非零特征值λ=n,其規(guī)范化特征向量w=(w1,w2,…,wn)T稱為權(quán)重向量,且αij=wi/wj。
(3)一致性檢驗(yàn)?,F(xiàn)實(shí)中判斷矩陣的賦值具有非等比性,階數(shù)大于2時(shí),很難滿足一致性要求,這時(shí)可通過λmax求權(quán)重向量,但矩陣偏離一致性的程度必須在可接受范圍內(nèi)。
設(shè)λmax為矩陣A=(αij)n×m的最大特征值,λmax≥n,等號(hào)當(dāng)且僅當(dāng)A為一致性矩陣時(shí)成立。因此,可根據(jù)λmax與n的偏差來檢驗(yàn)矩陣的一致性,λmax與n相差越多,A的一致性程度就越差。由此定義一致性指標(biāo)CI和平均隨機(jī)一致性指標(biāo)RI。
(11)
RI的值可參照平均隨機(jī)一致性指標(biāo)標(biāo)準(zhǔn)值表選取。取CI對(duì)RI的比值作為一致性檢驗(yàn)判別式,稱其為一致性比率CR,即:
(12)
如CR<0.1,則認(rèn)為該判斷矩陣通過一致性檢驗(yàn)。
3.1.1 保險(xiǎn)理賠數(shù)據(jù)統(tǒng)計(jì)計(jì)算 利用圖1的指標(biāo)體系,對(duì)該水電企業(yè)的保險(xiǎn)理賠數(shù)據(jù)進(jìn)行統(tǒng)計(jì)計(jì)算,以尋找其分布規(guī)律。
(1)準(zhǔn)則層風(fēng)險(xiǎn)。以設(shè)備風(fēng)險(xiǎn)為例,利用SPSS統(tǒng)計(jì)分析得到設(shè)備風(fēng)險(xiǎn)損失數(shù)據(jù)的概括性數(shù)據(jù)度量指標(biāo)見表1。
表1 設(shè)備風(fēng)險(xiǎn)損失數(shù)據(jù)度量指標(biāo)
由表1可看出,設(shè)備風(fēng)險(xiǎn)出險(xiǎn)次數(shù)多,損失數(shù)據(jù)為左偏、尖峰分布。選用半正態(tài)分布、對(duì)數(shù)正態(tài)分布、威布爾分布對(duì)其進(jìn)行擬合,對(duì)數(shù)正態(tài)分布擬合設(shè)備風(fēng)險(xiǎn)損失數(shù)據(jù)的P-P圖如圖2所示。由圖2可看出,對(duì)數(shù)正態(tài)分布的擬合效果較好,故采用對(duì)數(shù)正態(tài)分布對(duì)設(shè)備風(fēng)險(xiǎn)損失數(shù)據(jù)進(jìn)行擬合。
圖2 對(duì)數(shù)正態(tài)分布擬合設(shè)備風(fēng)險(xiǎn)損失數(shù)據(jù)的P-P圖
取理賠數(shù)據(jù)的自然對(duì)數(shù),可得其μ=1.464,σ=1.520,由此可得相應(yīng)的概率分布函數(shù)。劃分風(fēng)險(xiǎn)損失計(jì)算區(qū)間并利用Excel完成各區(qū)間分布累計(jì)概率與理論概率的計(jì)算,利用公式(3)計(jì)算各區(qū)間的概率貢獻(xiàn)率,將計(jì)算結(jié)果列于表2。
表2 風(fēng)險(xiǎn)損失區(qū)間計(jì)算表
按區(qū)間概率貢獻(xiàn)率將損失區(qū)間降序排序,利用公式(4)可得到理論概率為95.82%的風(fēng)險(xiǎn)損失區(qū)間為(0,60],進(jìn)一步劃分(30,60]并重復(fù)上述計(jì)算過程,可得到設(shè)備風(fēng)險(xiǎn)損失發(fā)生概率為95%的風(fēng)險(xiǎn)損失區(qū)間為(0,52.68]。由于風(fēng)險(xiǎn)損失不可能為0,取理賠數(shù)據(jù)中的最小值對(duì)風(fēng)險(xiǎn)損失區(qū)間進(jìn)行修正,將風(fēng)險(xiǎn)損失區(qū)間調(diào)整為[0.08,52.68]。
同理可求得環(huán)境風(fēng)險(xiǎn)和其他風(fēng)險(xiǎn)的風(fēng)險(xiǎn)損失區(qū)間,將計(jì)算結(jié)果列入表3。
表3 準(zhǔn)則層風(fēng)險(xiǎn)損失區(qū)間 104元
(2)指標(biāo)層風(fēng)險(xiǎn)。通過公式(1)、(5)可求出指標(biāo)層各風(fēng)險(xiǎn)因素的損失區(qū)間如表4。
表4 指標(biāo)層各風(fēng)險(xiǎn)因素的損失區(qū)間 104元
結(jié)合表5,通過公式(10)可得到準(zhǔn)則層判斷矩陣A以及指標(biāo)層判斷矩陣B1、B2、B3。
表5 準(zhǔn)則層和指標(biāo)層各風(fēng)險(xiǎn)因素風(fēng)險(xiǎn)后果嚴(yán)重程度及相對(duì)發(fā)生概率
結(jié)合2.2.4節(jié)中的理論可分別求解準(zhǔn)則層和指標(biāo)層的權(quán)重,將兩者相乘即可得到該水電企業(yè)的風(fēng)險(xiǎn)評(píng)估結(jié)果,如表6所示。
表6 某水電企業(yè)的風(fēng)險(xiǎn)評(píng)估結(jié)果
由表4、5可看出,2016-2019年該水電企業(yè)運(yùn)營過程中環(huán)境類風(fēng)險(xiǎn)和其他類風(fēng)險(xiǎn)損失較大,但發(fā)生概率較低,設(shè)備類風(fēng)險(xiǎn)損失較小,但發(fā)生概率較高,這符合水電站的風(fēng)險(xiǎn)統(tǒng)計(jì)規(guī)律[1]。
進(jìn)一步對(duì)照表6中的風(fēng)險(xiǎn)權(quán)重,權(quán)重最高的3個(gè)風(fēng)險(xiǎn)因素為設(shè)備故障、暴雨、洪水,其后果嚴(yán)重程度的評(píng)估值分別為4.00、5.21、4.00,相對(duì)發(fā)生概率的評(píng)估值分別為0.48、0.33、0.42,其中暴雨因素的相對(duì)發(fā)生概率低于洪水因素,但其后果嚴(yán)重程度高于后者,經(jīng)綜合計(jì)算分析可得其權(quán)重略高于后者,由此可見風(fēng)險(xiǎn)權(quán)重能綜合反映水電企業(yè)事故的嚴(yán)重程度與發(fā)生概率。水電企業(yè)不同風(fēng)險(xiǎn)的權(quán)重存在明顯差異,可為不同風(fēng)險(xiǎn)劃分不同的防范等級(jí),進(jìn)而分級(jí)制定具體的防控措施。
工程風(fēng)險(xiǎn)有一定的規(guī)律性,建立合理的風(fēng)險(xiǎn)防控體系可有效預(yù)防風(fēng)險(xiǎn),降低財(cái)產(chǎn)損失。本研究結(jié)合企業(yè)保險(xiǎn)理賠數(shù)據(jù),建立風(fēng)險(xiǎn)因素指標(biāo)體系,并進(jìn)一步分析得到各風(fēng)險(xiǎn)因素權(quán)重,這與前人通過文獻(xiàn)檢索、現(xiàn)場(chǎng)走訪等方式歸納風(fēng)險(xiǎn)因素指標(biāo)體系并結(jié)合專家打分計(jì)算得到風(fēng)險(xiǎn)權(quán)重[6-8,11,13]的方法不同,能更好地結(jié)合水電企業(yè)的特征,得到更有針對(duì)性的風(fēng)險(xiǎn)評(píng)估數(shù)據(jù)以供水電企業(yè)風(fēng)險(xiǎn)防控決策參考。但目前水電企業(yè)可用保險(xiǎn)理賠數(shù)據(jù)較少,基于保險(xiǎn)理賠數(shù)據(jù)的客觀評(píng)分法精度可能不足,本文引入風(fēng)險(xiǎn)相對(duì)發(fā)生概率與風(fēng)險(xiǎn)損失區(qū)間及其標(biāo)度確定方法以提高風(fēng)險(xiǎn)評(píng)估的準(zhǔn)確性。有學(xué)者在客觀權(quán)重的基礎(chǔ)上通過博弈模型[12]耦合主觀權(quán)重,該方法亦能較好地提高評(píng)估精度。
本研究得出了某水電企業(yè)各風(fēng)險(xiǎn)因素的權(quán)重與排序,但并未對(duì)風(fēng)險(xiǎn)成因及其防控措施作深入研究,水電企業(yè)長(zhǎng)周期安全穩(wěn)定運(yùn)行的技術(shù)及管理模式還需進(jìn)一步探討。環(huán)境風(fēng)險(xiǎn)和設(shè)備風(fēng)險(xiǎn)是眾多水電企業(yè)長(zhǎng)周期安全穩(wěn)定運(yùn)行中亟需改進(jìn)的突出問題,環(huán)境風(fēng)險(xiǎn)損失大、機(jī)理復(fù)雜、難以避免,有待進(jìn)一步研究其事故機(jī)理及預(yù)警方法以減小事故損失,設(shè)備風(fēng)險(xiǎn)出險(xiǎn)頻率高、類別多,需深入研究運(yùn)營管理模式以減小設(shè)備事故概率。
(1)提出了一種基于保險(xiǎn)理賠數(shù)據(jù)的水電企業(yè)風(fēng)險(xiǎn)評(píng)估方法,該方法通過對(duì)保險(xiǎn)理賠數(shù)據(jù)的統(tǒng)計(jì)分析直接得到各風(fēng)險(xiǎn)因素的權(quán)重,能較好地規(guī)避傳統(tǒng)風(fēng)險(xiǎn)評(píng)估方法的主觀性和隨機(jī)性,且保險(xiǎn)理賠數(shù)據(jù)可與風(fēng)險(xiǎn)評(píng)估數(shù)據(jù)實(shí)現(xiàn)同步實(shí)時(shí)更新,評(píng)估精度將隨時(shí)間的積累而不斷提高,具有實(shí)用意義。
(2)引入了風(fēng)險(xiǎn)相對(duì)發(fā)生概率與風(fēng)險(xiǎn)損失區(qū)間及其標(biāo)度的確定方法,適用于水電企業(yè)出險(xiǎn)數(shù)據(jù)少的情況,且其求解方法簡(jiǎn)單,有較強(qiáng)的可操作性,可推廣至水電企業(yè)的基層管理人員。
(3)利用基于保險(xiǎn)理賠數(shù)據(jù)的水電企業(yè)風(fēng)險(xiǎn)評(píng)估方法對(duì)某水電企業(yè)進(jìn)行了運(yùn)營風(fēng)險(xiǎn)評(píng)估,得到了該企業(yè)的風(fēng)險(xiǎn)因素權(quán)重及排序,可為該企業(yè)的風(fēng)險(xiǎn)防控工作提供參考。結(jié)果表明:該水電企業(yè)應(yīng)加大在設(shè)備故障、暴雨、洪水等方面的風(fēng)險(xiǎn)防控投入,適當(dāng)減緩在電壓波動(dòng)、絕緣擊穿、過電壓等方面的風(fēng)險(xiǎn)防控投入??筛鶕?jù)風(fēng)險(xiǎn)權(quán)重排序,有區(qū)別地制定不同等級(jí)的風(fēng)險(xiǎn)防控措施。