趙紅澤，張超力

(中國礦業(yè)大學(xué)(北京) 能源與礦業(yè)學(xué)院，北京 100083)

我國近年來煤炭百萬噸死亡率逐年下降，煤礦安全生產(chǎn)形勢實現(xiàn)根本性好轉(zhuǎn)，但與發(fā)達(dá)國家相比，我國的煤礦安全生產(chǎn)水平依然有很大差距。由于煤礦生產(chǎn)環(huán)境的復(fù)雜性，特別是近年來煤礦開采深度不斷加大，煤礦安全生產(chǎn)依然面臨嚴(yán)峻的挑戰(zhàn)[1]。從2014年的馬航事件，2015年的“東方之星”客輪翻沉事件，天津港爆炸事件，再到今年的新冠疫情，引發(fā)了社會對我國公共安全應(yīng)急救援能力建設(shè)問題的思考。著眼于構(gòu)建統(tǒng)一領(lǐng)導(dǎo)、權(quán)責(zé)一致、權(quán)威高效的國家應(yīng)急能力體系，在2018年十三屆全國人大一次會議上，國務(wù)院組建應(yīng)急管理部[2]。新部門的成立在很大程度上有利于提升應(yīng)對復(fù)雜性突發(fā)事件的能力，是我國加強現(xiàn)代應(yīng)急管理的一個重要舉措[3]。因此，如何在事故發(fā)生初期結(jié)合有限的信息，對應(yīng)急救援過程中所需要的應(yīng)急物資需求進(jìn)行預(yù)測，保證應(yīng)急救援物資準(zhǔn)確高效地調(diào)配至事故地點，最大程度上減輕事故造成的損失，挽救更多的生命，對于提高煤礦企業(yè)應(yīng)急救援能力具有重大意義[4-6]。因此，研究煤礦應(yīng)急救援物資需求的問題不僅具有重要的理論意義，而且具有重大的社會意義[7]。

1 案例推理理論

案例推理技術(shù)是人工智能發(fā)展成熟的一個分支，CBR基于這樣一個事實：人們在解決新的問題時總是先試圖從以往的經(jīng)驗中找到案例從而得到答案或者啟發(fā)。例如對于煤礦事故應(yīng)急救援，對于相同種類的事故，在事故救援中的一些救援措施具有很強的可參考性。

CBR正是對新問題求解的抽象解釋：把當(dāng)前問題與過往曾經(jīng)成功解決的問題進(jìn)行匹配，從中得到答案或者啟發(fā)，在案例匹配的過程中，可能需要對案例庫中檢索到匹配案例的解決方案進(jìn)行修改，使其更加符合當(dāng)前問題的具體情況。

案例推理可以看成一個4R的循環(huán)過程，即案例檢索、案例重用、案例修改調(diào)整、案例學(xué)習(xí)四個步驟的循環(huán)，其流程如圖1所示[8]。

圖1 案例推理流程

2 基于案例推理的煤礦應(yīng)急物資需求預(yù)測模型構(gòu)建

2.1 需求預(yù)測方法對比確定

煤礦應(yīng)急救援物資需求預(yù)測是一個多輸入單輸出或者多輸入多輸出的一種預(yù)測問題，并且輸入和輸出的結(jié)果具有邏輯上的因果聯(lián)系。經(jīng)過對預(yù)測方法的研究，對于多個因變量的預(yù)測方法主要有多元線性回歸預(yù)測、案例推理預(yù)測、神經(jīng)網(wǎng)絡(luò)預(yù)測等方法。各預(yù)測方法優(yōu)缺點見表1。

表1 預(yù)測方法對比表

由于多元線性回歸預(yù)測對于自變量和預(yù)測量之間的線性關(guān)系要求高，樣本數(shù)據(jù)的質(zhì)量和數(shù)量要求較高；神經(jīng)網(wǎng)絡(luò)預(yù)測方法的結(jié)果不穩(wěn)定，具有波動性大等限制和不足，樣本數(shù)量要求大；案例推理是一種定性、定量結(jié)合的方法，學(xué)習(xí)能力和適應(yīng)性較強，求解效率高，在煤礦事故應(yīng)急救援的過程中，不同事故的危害程度不同，也決定了應(yīng)急資源需求的多少。所以不同煤礦事故的屬性相似性使得采用案例推理理論在煤礦事故應(yīng)急救援資源需求預(yù)測領(lǐng)域具有了理論基礎(chǔ)，因此，根據(jù)煤礦事故應(yīng)急救援的特點，本文采用案例推理對煤礦應(yīng)急救援物資需求進(jìn)行預(yù)測。

2.2 煤礦事故應(yīng)急救援案例組成與表示

煤礦事故應(yīng)急救援同地震、水災(zāi)等自然災(zāi)害事故的救援不同，煤礦事故種類繁多，按照事故性質(zhì)的不同可分為頂板、瓦斯、機電、運輸、放炮、火災(zāi)、水害和其他八類事故，對于不同類型的事故，也因為發(fā)生原因、環(huán)境、地理位置、地質(zhì)條件的不同導(dǎo)致處理方式的不同，例如，對于瓦斯事故，根據(jù)事故原因的不同又分為瓦斯爆炸、煤塵爆炸、煤與瓦斯突出等類型。同時，對于不同類型的煤礦事故，應(yīng)急救援所需的救援物資也不同，例如，對于火災(zāi)事故，不僅需要氧氣瓶、呼吸機等生命救援設(shè)備，也需要針對火災(zāi)事故的CO2滅火器。對于透水事故，需要水泥進(jìn)行井下筑壩和大功率水泵進(jìn)行抽水。針對煤礦應(yīng)急救援案例的復(fù)雜性，本文采用公共屬性與私有屬性結(jié)合的案例表示方式，公共類型屬性為事故類型、年產(chǎn)量、傷亡人數(shù)、經(jīng)濟損失等，私有屬性為不同類型煤礦事故的細(xì)分屬性。

煤礦應(yīng)急救援案例庫包括案例屬性和應(yīng)急資源數(shù)量，可以將煤礦應(yīng)急救援案例模型化[9]，如式(1)所示。

case(F，W，R)

(1)

式中，F(xiàn)為煤礦應(yīng)急救援案例屬性集；W為屬性權(quán)重集；R為案例應(yīng)急資源需求集合。

2.3 屬性權(quán)重的確定

在建立煤礦應(yīng)急救援的案例庫之后，當(dāng)新的事故發(fā)生時，就可以在案例庫中檢索匹配最佳案例，在案例檢索之前還要確定各屬性值的權(quán)重，目前屬性權(quán)重的確定主要有主觀賦權(quán)法和客觀賦權(quán)法兩類，由于上述兩種方法各有利弊，為了能夠全面反映屬性的重要性，本文采用二者相結(jié)合的方法，首先采用AHP方法計算各屬性的權(quán)重，然后利用熵技術(shù)對得到的權(quán)重系數(shù)進(jìn)行修正，把修正后的系數(shù)作為最終的權(quán)重系數(shù)。

2.3.1 層次分析法確定權(quán)重

第一步：構(gòu)建層次結(jié)構(gòu)模型。根據(jù)所描述的問題建立目標(biāo)層、準(zhǔn)則層、方案層三層層次結(jié)構(gòu)模型。

第二步：構(gòu)建判斷矩陣C。

第三步：進(jìn)行層次單排列計算。根據(jù)式(2)計算矩陣每行標(biāo)度乘積的n次方根，通過式(3)進(jìn)行歸一化處理，得到指標(biāo)權(quán)重。

第四步：通過一致性檢驗，以確定其是否有滿意的一致性，如果不滿足一致性要求，則要調(diào)整判斷矩陣的數(shù)據(jù)，直到滿足一致性要求為止。

2.3.2 基于熵權(quán)法修正權(quán)系數(shù)

熵最早由信息論之父香農(nóng)引入信息論中，信息熵越大，說明該屬性提供的信息量越大，在案例中所起的作用越大，權(quán)重也應(yīng)該越大。因此信息熵可作為一種客觀賦權(quán)法，結(jié)合層次分析法對指標(biāo)權(quán)重進(jìn)行修正，最終得到客觀的指標(biāo)權(quán)重。

第一步：計算屬性的熵φj和差異系數(shù)ηj。首先根據(jù)式(4)對判斷矩陣做歸一化處理，第j項的熵φj的計算公式如(5)所示。其中，k=1/lnm，m為案例庫中指標(biāo)屬性的個數(shù)，從而得到第j項屬性的差異系數(shù)為：ηj=1-φj。

第二步：求修正系數(shù)θj。

第三步：根據(jù)修正系數(shù)θj得到通過熵技術(shù)修正后的權(quán)重向量δj，如式(7)所示。

通過上述計算完成了基于熵技術(shù)的屬性權(quán)值修正，可以進(jìn)行下一步的案例相似度計算。

2.4 改進(jìn)的案例檢索方法

目前常用的案例檢索算法有最近鄰法、神經(jīng)網(wǎng)絡(luò)法、知識引導(dǎo)法、歸納索引法。由于煤礦應(yīng)急救援案例庫的復(fù)雜性，屬性的多樣性，不僅包含數(shù)值型屬性，而且包含文本型屬性，難以在多屬性情況下得出一個綜合相似度度量值。本文針對數(shù)值型屬性，采用基于歐氏距離的最近鄰算法計算相似度，針對文本型屬性，以模計算其相似度，最后根據(jù)數(shù)值型屬性和文本型屬性的權(quán)重計算案例的全局相似度。

第一步，對目標(biāo)案例及源案例中的數(shù)值型屬性根據(jù)式(8)進(jìn)行歸一化處理。

第二步，計算局部相似度。對于數(shù)值型屬性，采用最近鄰算法計算其局部相似度，可以用歐氏距離計算兩個案例特征屬性之間的距離，計算出它們之間的相似度，其相似度公式如(9)所示。對于文本型屬性，一種是屬性取值僅有兩個狀態(tài)，一種是有多于兩個的狀態(tài)值，用模計算其相似度，其相似度定義如式(10)所示。

第三步，計算全局相似度。全局相似度為數(shù)值型屬性局部相似度和文本型屬性局部相似度的加權(quán)之和。

最終計算出的全局相似度越大，源案例與目標(biāo)案例越匹配，最后找出相似度最大的源案例即為最相似案例，則檢索到的最相似案例的物資調(diào)配情況可以作為目標(biāo)案例所需資源量進(jìn)行應(yīng)急救援物資調(diào)配。

3 實例驗證

3.1 案例概述

2014年4月27日內(nèi)蒙古興通煤礦11B07綜采工作面發(fā)生瓦斯燃燒事故，由于火勢發(fā)展迅速，近距離無法封閉，被迫進(jìn)行全礦封閉，事故造成經(jīng)濟損失兩千萬元。事故發(fā)生后，平莊煤業(yè)公司救護(hù)大隊出動液態(tài)二氧化碳滅火裝備、氣體化驗車等救援裝備，對已封閉的礦井利用灌注液態(tài)二氧化碳的方法進(jìn)行惰化，合計注入151.5t液態(tài)二氧化碳，既達(dá)到了抑制瓦斯煤塵爆炸、保護(hù)救災(zāi)人員安全的目的，又收到了快速撲滅災(zāi)區(qū)火源、降低災(zāi)區(qū)溫度的顯著效果，啟封過程中氧氣濃度始終在10% 以下，溫度低于18℃，為救護(hù)隊員啟封火區(qū)創(chuàng)造了相對有利的條件。本文以煤礦火災(zāi)事故應(yīng)急救援物資需求量的預(yù)測為實例，對構(gòu)建的基于案例推理的煤礦應(yīng)急救援物資需求預(yù)測模型進(jìn)行驗證。

3.2 層次分析法確定權(quán)重

這里選擇我國過去礦井火災(zāi)事故案例作為案例庫[10-12]，對新案例進(jìn)行預(yù)測分析。具體的火災(zāi)事故案例見表2。

表2 煤礦火災(zāi)事故應(yīng)急救援案例

通過向煤礦安全技術(shù)人員、行業(yè)專家發(fā)放調(diào)查問卷，收集評判數(shù)據(jù)得到了指標(biāo)判斷矩陣如下：

計算判斷矩陣一致性指標(biāo)：

隨機一致性比率CR小于0.1，表明判斷矩陣符合一致性檢驗指標(biāo)，基本符合現(xiàn)場實際情況，一致性檢驗較好，得到煤礦火災(zāi)事故案例的五個屬性權(quán)重為：

ω1=0.0629，ω2=0.1150，ω3=0.2494，

ω4=0.4063，ω5=0.1664

3.3 基于信息熵技術(shù)的案例權(quán)重修正

根據(jù)上文方法對判斷矩陣進(jìn)行歸一化處理，得到屬性權(quán)重的歸一化判斷矩陣：

將歸一化后的判斷矩陣代入式(6)—式(7)中，依次得到信息熵φj，差異系數(shù)η，修正系數(shù)θj，使用熵信息技術(shù)修正后得權(quán)重δj，各參數(shù)數(shù)據(jù)見表3。

表3 權(quán)重修正參數(shù)表

最終得到經(jīng)過熵信息技術(shù)修正后的煤礦火災(zāi)事故應(yīng)急救援屬性權(quán)重為：年產(chǎn)量δ1=0.0390，傷亡人數(shù)δ2=0.0850，可燃物種類δ3=0.3156，火災(zāi)類型δ4=0.3288，持續(xù)時間δ5=0.2316。

接下來對目標(biāo)案例和源案例中數(shù)值型屬性進(jìn)行歸一化處理，然后計算相似度，按照式(8)—式(11)分別計算數(shù)值型屬性和文本型屬性的局部相似度以及全局相似度，計算結(jié)果見表4。

表4 案例相似度計算表

經(jīng)過計算，案例D為目標(biāo)案例的最佳匹配案例，案例D在應(yīng)急救援過程中調(diào)配液態(tài)二氧化碳罐車10輛，即190t，目標(biāo)案例興通煤礦“4.27”事故救援過程中，實際使用的液態(tài)二氧化碳罐車8輛，151.5t。經(jīng)過案例分析，可見構(gòu)建的方法符合實際情況。如果在事故發(fā)生初期，按照該方法對應(yīng)急資源進(jìn)行預(yù)測，將大大增加救援的科學(xué)性和有效性。

4 結(jié) 論

1)針對煤礦事故應(yīng)急救援的特點，建立了公共屬性結(jié)合私有屬性的煤礦應(yīng)急救援物資需求案例庫，為煤礦應(yīng)急救援物資需求案例庫的構(gòu)建提供了框架。

2)由于傳統(tǒng)的案例推理檢索方法無法滿足煤礦應(yīng)急救援物資需求案例檢索的要求，本文構(gòu)建了針對文本屬性以模求其相似度，針對數(shù)值型屬性以最近鄰法求其相似度的改進(jìn)混合案例檢索方法。

3)結(jié)合實際煤礦應(yīng)急救援案例，建立煤礦火災(zāi)事故應(yīng)急救援物資需求案例庫，結(jié)合本文的方法進(jìn)行計算，結(jié)果顯示所需物資預(yù)測量符合實際情況，驗證了模型的有效性。

4)基于以上研究成果，開發(fā)煤礦應(yīng)急物資儲備及調(diào)度虛擬演練系統(tǒng)，可有效提高煤礦應(yīng)急物資儲備和調(diào)度的可靠性及時效性，結(jié)合GIS系統(tǒng)和三維可視化系統(tǒng)應(yīng)用于教學(xué)培訓(xùn)中，可直觀展示應(yīng)急物資的滿足程度、儲存狀態(tài)、調(diào)配流程和欠缺情況，為礦山應(yīng)急救援等課程的實驗教學(xué)提供支撐，豐富教學(xué)形式與內(nèi)容。