李 爽，李丁煒，猶夢潔

（1.中國礦業(yè)大學(xué) 管理學(xué)院，江蘇 徐州221008；2.中國礦業(yè)大學(xué) 安全科學(xué)與應(yīng)急管理研究中心，江蘇 徐州 221008）

根據(jù)國家煤礦安監(jiān)局網(wǎng)站信息顯示，2015 年到2018 年我國共發(fā)生煤礦生產(chǎn)重大事故20 起，共造成288 人死亡，煤礦安全生產(chǎn)形勢依然嚴(yán)峻。國內(nèi)外學(xué)者對(duì)煤礦安全生產(chǎn)預(yù)控工作開展了大量研究，構(gòu)建了應(yīng)用人工智能技術(shù)的煤礦瓦斯突出、沖擊地壓、地板水災(zāi)等災(zāi)害的預(yù)測系統(tǒng)[1-3]，并將卷積神經(jīng)網(wǎng)絡(luò)和自動(dòng)編碼器等深度學(xué)習(xí)模型應(yīng)用于健康預(yù)警、設(shè)備預(yù)警和環(huán)境預(yù)警[4-6]。但目前的煤礦安全生產(chǎn)預(yù)控系統(tǒng)仍存在缺陷：風(fēng)險(xiǎn)預(yù)測預(yù)警方法相對(duì)分散和獨(dú)立，處理的風(fēng)險(xiǎn)與數(shù)據(jù)類型較少，無法保證預(yù)警的準(zhǔn)確度；在風(fēng)險(xiǎn)處于產(chǎn)生的臨界點(diǎn)或已經(jīng)產(chǎn)生時(shí)進(jìn)行報(bào)警，無法真正實(shí)現(xiàn)風(fēng)險(xiǎn)預(yù)控；全面系統(tǒng)的安全生產(chǎn)風(fēng)險(xiǎn)智能預(yù)警平臺(tái)尚未建立等。為了實(shí)現(xiàn)煤礦風(fēng)險(xiǎn)的預(yù)控，需要全面掌握煤礦在某一時(shí)間各項(xiàng)風(fēng)險(xiǎn)發(fā)生的分布情況及其對(duì)煤礦整體安全的影響程度。安全態(tài)勢感知（security situation awareness，SSA）是通過智能技術(shù)從時(shí)間和空間維度來提取所有與安全相關(guān)的因素，對(duì)其進(jìn)行分析理解來感知整體的安全狀態(tài)并對(duì)未來的安全狀態(tài)進(jìn)行預(yù)測[7],在航空和軍事、電網(wǎng)、計(jì)算機(jī)網(wǎng)絡(luò)、醫(yī)學(xué)研究[8-10]等領(lǐng)域有較多的研究。將SSA 的概念引入到煤礦安全生產(chǎn)中，從而構(gòu)建智能感知、風(fēng)險(xiǎn)動(dòng)態(tài)評(píng)估與智能預(yù)警一體化的礦山安全態(tài)勢感知系統(tǒng)。

1 煤礦安全態(tài)勢感知系統(tǒng)

煤礦安全態(tài)勢的評(píng)估因素著眼于某一時(shí)間、空間內(nèi)各種風(fēng)險(xiǎn)產(chǎn)生的可能性。其基本流程為：通過對(duì)煤礦事故因素動(dòng)態(tài)變化、內(nèi)外部因素相互作用的分析，總結(jié)對(duì)煤礦存在潛在威脅的內(nèi)容，提取煤礦安全態(tài)勢要素；探究多因素作用下的煤礦事故誘發(fā)機(jī)理，分析識(shí)別煤礦風(fēng)險(xiǎn)源頭及其影響，并建立態(tài)勢分析模型；通過態(tài)勢分析模型識(shí)別已經(jīng)獲取的煤礦數(shù)據(jù)中所隱含的風(fēng)險(xiǎn)源對(duì)煤礦安全態(tài)勢的影響，評(píng)估煤礦系統(tǒng)中存在的潛在威脅，以從全局視角提升煤礦對(duì)安全威脅的識(shí)別、分析、響應(yīng)能力。

現(xiàn)有的煤礦安全監(jiān)控和報(bào)警系統(tǒng)核心思想在于當(dāng)風(fēng)險(xiǎn)處于產(chǎn)生的臨界點(diǎn)或已經(jīng)產(chǎn)生時(shí)，通過實(shí)時(shí)的數(shù)據(jù)監(jiān)控和報(bào)警功能及時(shí)發(fā)現(xiàn)并控制風(fēng)險(xiǎn)。區(qū)別于此，煤礦安全態(tài)勢預(yù)測的目的在于，在煤礦進(jìn)入風(fēng)險(xiǎn)狀態(tài)之前，對(duì)影響煤礦風(fēng)險(xiǎn)的主要因素的變化趨勢進(jìn)行分析預(yù)測，在對(duì)未來節(jié)點(diǎn)這些數(shù)據(jù)的預(yù)測的基礎(chǔ)上評(píng)估該未來節(jié)點(diǎn)風(fēng)險(xiǎn)產(chǎn)生的可能性?；谶@種目的，需要構(gòu)建適應(yīng)于度量煤礦安全風(fēng)險(xiǎn)態(tài)勢的預(yù)測模型，分析預(yù)測煤礦內(nèi)風(fēng)險(xiǎn)產(chǎn)生可能性的變化趨勢，并動(dòng)態(tài)生成煤礦安全風(fēng)險(xiǎn)態(tài)勢圖，以支持系統(tǒng)對(duì)于風(fēng)險(xiǎn)的提前識(shí)別和預(yù)警。

1.1 系統(tǒng)整體架構(gòu)

煤礦安全態(tài)勢感知系統(tǒng)主要由物聯(lián)網(wǎng)信息收集云平臺(tái)數(shù)據(jù)管理、煤礦安全態(tài)勢預(yù)測模型、信息可視化平臺(tái)、預(yù)警通知4 個(gè)部分組成，煤礦安全態(tài)勢感知系統(tǒng)技術(shù)架構(gòu)圖如圖1。系統(tǒng)利用基于物聯(lián)網(wǎng)的信息收集云平臺(tái)獲取煤礦實(shí)時(shí)監(jiān)測數(shù)據(jù)和非實(shí)時(shí)隱患排查數(shù)據(jù)，依據(jù)安全態(tài)勢預(yù)測指標(biāo)體系構(gòu)建安全態(tài)勢預(yù)測模型，對(duì)監(jiān)測數(shù)據(jù)進(jìn)行挖掘和分析，進(jìn)而得出安全態(tài)勢預(yù)測結(jié)論；通過可視化平臺(tái)將態(tài)勢預(yù)測結(jié)果可視化展現(xiàn)，并經(jīng)由預(yù)警通知模塊通知相關(guān)的風(fēng)險(xiǎn)負(fù)責(zé)人。

圖1 煤礦安全態(tài)勢感知系統(tǒng)技術(shù)架構(gòu)圖Fig.1 Technical architecture diagram of coal mine safety situational awareness system

1.2 物聯(lián)網(wǎng)信息收集云平臺(tái)數(shù)據(jù)管理

物聯(lián)網(wǎng)是互聯(lián)網(wǎng)、傳統(tǒng)電信網(wǎng)等信息承載體，讓所有能行使獨(dú)立功能的普通物體實(shí)現(xiàn)互聯(lián)互通的網(wǎng)絡(luò)。物聯(lián)網(wǎng)及其信息處理、存儲(chǔ)云平臺(tái)是基于大數(shù)據(jù)的煤礦安全態(tài)勢預(yù)測系統(tǒng)的起點(diǎn)。為實(shí)現(xiàn)對(duì)煤礦整體安全態(tài)勢的評(píng)估預(yù)測和及時(shí)預(yù)警，需要對(duì)煤礦實(shí)時(shí)與非實(shí)時(shí)安全監(jiān)測數(shù)據(jù)進(jìn)行收集和分析。基于物聯(lián)網(wǎng)技術(shù)，通過煤礦監(jiān)測傳感器收集海量的煤礦實(shí)時(shí)與非實(shí)時(shí)監(jiān)測數(shù)據(jù)并將其上傳至計(jì)算云平臺(tái)。計(jì)算云平臺(tái)通過對(duì)大數(shù)據(jù)的挖掘處理，提取并為煤礦安全態(tài)勢感知系統(tǒng)提供其所需指標(biāo)的數(shù)據(jù)，以進(jìn)行煤礦的安全態(tài)勢分析和風(fēng)險(xiǎn)預(yù)測。

物聯(lián)網(wǎng)通過監(jiān)測傳感器收集煤礦監(jiān)測數(shù)據(jù)，上傳至云平臺(tái)分析處理。煤礦安全態(tài)勢感知系統(tǒng)通過云平臺(tái)數(shù)據(jù)管理模塊借助云平臺(tái)提供的數(shù)據(jù)接口實(shí)現(xiàn)云平臺(tái)與系統(tǒng)的數(shù)據(jù)交流，以實(shí)現(xiàn)應(yīng)用于安全態(tài)勢預(yù)測的煤礦監(jiān)測數(shù)據(jù)的下載，以及煤礦預(yù)測預(yù)警結(jié)果的上傳。實(shí)現(xiàn)煤礦安全數(shù)據(jù)集約化、網(wǎng)絡(luò)化遠(yuǎn)程管理，充分發(fā)揮物聯(lián)網(wǎng)及云技術(shù)在煤礦安全態(tài)勢預(yù)測中的作用。物聯(lián)網(wǎng)信息收集云平臺(tái)數(shù)據(jù)管理主要包括物聯(lián)網(wǎng)數(shù)據(jù)獲取以及本地?cái)?shù)據(jù)上傳2 部分。

1.3 煤礦安全態(tài)勢預(yù)測模型

煤礦安全態(tài)勢預(yù)測模型為系統(tǒng)的核心部分，在構(gòu)建安全態(tài)勢預(yù)測模型的基礎(chǔ)上對(duì)煤礦整體安全態(tài)勢進(jìn)行分析和預(yù)測，以實(shí)現(xiàn)對(duì)煤礦可能風(fēng)險(xiǎn)的預(yù)先知悉、評(píng)估和診斷，為煤礦安全建設(shè)提供依據(jù)。并將過往、當(dāng)前以及未來時(shí)間節(jié)點(diǎn)和安全態(tài)勢評(píng)估綜合，形成一個(gè)體現(xiàn)對(duì)煤礦風(fēng)險(xiǎn)的總體評(píng)估的安全態(tài)勢走向圖，以及各項(xiàng)風(fēng)險(xiǎn)產(chǎn)生可能性變化趨勢圖，通過信息可視化平臺(tái)展現(xiàn)給使用者。在安全態(tài)勢預(yù)測模型的評(píng)估結(jié)果中，如果未來時(shí)間節(jié)點(diǎn)某一區(qū)塊存在某種風(fēng)險(xiǎn)，系統(tǒng)將對(duì)風(fēng)險(xiǎn)進(jìn)行預(yù)警，并通過通知模塊告知相關(guān)負(fù)責(zé)人。

1.4 信息可視化平臺(tái)

信息可視化平臺(tái)是系統(tǒng)與使用者交互的接口，將管理者在煤礦生產(chǎn)管理中所要的各種信息（煤礦安全態(tài)勢預(yù)測結(jié)果、煤礦風(fēng)險(xiǎn)預(yù)測結(jié)果、煤礦作業(yè)信息、系統(tǒng)管理信息）等內(nèi)容進(jìn)行綜合可視化展現(xiàn)。

1）煤礦作業(yè)信息展現(xiàn)。模塊將煤礦作業(yè)信息以圖表的形式展示給礦區(qū)負(fù)責(zé)人，并支持其對(duì)作業(yè)信息的新增和更改。礦區(qū)需要進(jìn)行作業(yè)時(shí)，負(fù)責(zé)人創(chuàng)建作業(yè)信息，系統(tǒng)通過信息通知模塊通知相應(yīng)礦區(qū)的負(fù)責(zé)人以協(xié)調(diào)工作。同時(shí)模塊還會(huì)協(xié)調(diào)安全態(tài)勢預(yù)測模塊，在作業(yè)期間停止對(duì)作業(yè)所影響的風(fēng)險(xiǎn)指標(biāo)進(jìn)行預(yù)測預(yù)警。

2）系統(tǒng)管理信息展現(xiàn)。模塊將系統(tǒng)用戶、用戶權(quán)限、系統(tǒng)參數(shù)設(shè)定、相關(guān)負(fù)責(zé)人設(shè)定、煤礦區(qū)塊劃分、預(yù)測預(yù)警指標(biāo)體系、安全態(tài)勢模型參數(shù)等系統(tǒng)管理信息展示給系統(tǒng)管理者，并持支其對(duì)各項(xiàng)參數(shù)的修改。同時(shí)對(duì)后臺(tái)管理模塊的一切操作進(jìn)行記錄，以便后續(xù)的審計(jì)工作開展。

1.5 預(yù)警通知

當(dāng)作業(yè)人員發(fā)布新的作業(yè)后，模塊根據(jù)作業(yè)信息中的作業(yè)位置，及時(shí)告知后臺(tái)管理模塊負(fù)責(zé)人管理中所設(shè)定的改作業(yè)區(qū)塊的負(fù)責(zé)人，以便區(qū)塊負(fù)責(zé)人協(xié)調(diào)工作。當(dāng)煤礦處于風(fēng)險(xiǎn)狀態(tài)并進(jìn)行風(fēng)險(xiǎn)預(yù)警后，模塊立即通知存在風(fēng)險(xiǎn)隱患的區(qū)塊的負(fù)責(zé)人，同時(shí)通知諸如煤礦礦長等管理人員。信息通知方式采用短信通知、郵件通知、系統(tǒng)手機(jī)客戶端多種方式并行通知。

2 煤礦安全態(tài)勢預(yù)測模型方案

煤礦安全態(tài)勢預(yù)測模型的建立包括2 個(gè)部分，安全態(tài)勢預(yù)測指標(biāo)體系和安全態(tài)勢預(yù)測模型。指標(biāo)體系的建立基于對(duì)煤礦事故記錄的分析以提取煤礦風(fēng)險(xiǎn)影響因子，并通過貝葉斯網(wǎng)絡(luò)（Bayesian Network，BN）提取事故致因鏈以區(qū)分多級(jí)指標(biāo)。預(yù)測模型的建立基于粗糙集理論（Rough Set, RS）和支持向量機(jī)（Support Vector Machine, SVM）的結(jié)合應(yīng)用，通過機(jī)器學(xué)習(xí)的方式對(duì)各項(xiàng)風(fēng)險(xiǎn)直接影響因素進(jìn)行一級(jí)預(yù)測，并根據(jù)貝葉斯網(wǎng)絡(luò)模型對(duì)風(fēng)險(xiǎn)發(fā)生的可能性進(jìn)行二級(jí)預(yù)測，以實(shí)現(xiàn)煤礦安全態(tài)勢的預(yù)測和風(fēng)險(xiǎn)的預(yù)警。

2.1 預(yù)處理階段

預(yù)處理階段的目的是明確煤礦存在的風(fēng)險(xiǎn)，確定進(jìn)行煤礦安全態(tài)勢預(yù)測所需要監(jiān)測獲得的數(shù)據(jù)。系統(tǒng)依據(jù)全國煤礦事故數(shù)據(jù)提取出5 個(gè)主要煤礦事故：瓦斯、頂板、運(yùn)輸、水災(zāi)以及機(jī)電事故,以及會(huì)導(dǎo)致各項(xiàng)煤礦事故的危險(xiǎn)源，并結(jié)合專家經(jīng)驗(yàn)習(xí)構(gòu)建貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)。通過貝葉斯網(wǎng)絡(luò)圖可以得到各項(xiàng)事故的致因鏈以及各個(gè)節(jié)點(diǎn)的概率分布。依據(jù)事故致因鏈，可以進(jìn)一步區(qū)分導(dǎo)致事故的直接因素也即一級(jí)指標(biāo)，以及間接因素也即多級(jí)指標(biāo)。具體流程如下：①根據(jù)全國煤礦安全事故記錄或地區(qū)煤礦安全事故提取出煤礦主要風(fēng)險(xiǎn)及危險(xiǎn)源，并對(duì)事故樣本進(jìn)行參數(shù)學(xué)習(xí)和結(jié)構(gòu)學(xué)習(xí)，生成貝葉斯網(wǎng)絡(luò)；②地區(qū)煤礦專家根據(jù)經(jīng)驗(yàn)調(diào)整概率網(wǎng)絡(luò)圖，分析獲得各風(fēng)險(xiǎn)的致因鏈以及各節(jié)點(diǎn)概率；③根據(jù)貝葉斯網(wǎng)絡(luò)概率圖和風(fēng)險(xiǎn)致因鏈，區(qū)分影響煤礦風(fēng)險(xiǎn)的一級(jí)指標(biāo)和多級(jí)指標(biāo)，并得到風(fēng)險(xiǎn)二級(jí)預(yù)測模型，用于后續(xù)的安全態(tài)勢預(yù)測；④根據(jù)物聯(lián)網(wǎng)收集的監(jiān)測數(shù)據(jù)動(dòng)態(tài)調(diào)整概率網(wǎng)絡(luò)，逐漸得到最適應(yīng)地區(qū)煤礦的指標(biāo)體系及概率分布。以煤礦瓦斯事故為例實(shí)施預(yù)處理階段，以2000 年至2018 年間全國瓦斯事故為樣本，總結(jié)并歸納事故原因，構(gòu)建瓦斯事故的貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)。根據(jù)貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)提取以下信息。

2.1.1 多級(jí)指標(biāo)體系

由貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)分析得知，會(huì)導(dǎo)致瓦斯事故風(fēng)險(xiǎn)的直接因素有明火、電器失爆、瓦斯積聚3 項(xiàng)，稱一級(jí)指標(biāo)，記為FWi0對(duì)瓦斯事故各條一級(jí)指標(biāo)致因鏈進(jìn)行分析，確認(rèn)對(duì)一級(jí)指標(biāo)存在直接或間接影響的因素，稱多級(jí)指標(biāo)，記為MWi0一級(jí)指標(biāo)明火FW1，對(duì)應(yīng)的多級(jí)指標(biāo)為：安全培訓(xùn)不到位MW11、安全執(zhí)行能力弱MW12、技術(shù)管理不到位MW13、礦工安全意識(shí)低MW14、摩擦撞擊MW15、“三違”行為MW16共6 項(xiàng)。一級(jí)指標(biāo)電器失爆F(xiàn)W2，對(duì)應(yīng)的多級(jí)指標(biāo)為MW11、MW12、MW13、MW14、MW16、設(shè)備老化MW17等6 項(xiàng)。一級(jí)指標(biāo)瓦斯積聚FW3，對(duì)應(yīng)的多級(jí)指標(biāo)為MW13、通風(fēng)系統(tǒng)混亂MW18、瓦斯監(jiān)控系統(tǒng)故障MW19、瓦斯涌出MW110等4 項(xiàng)。

2.1.2 風(fēng)險(xiǎn)二級(jí)預(yù)測模型

對(duì)于t 時(shí)刻瓦斯事故風(fēng)險(xiǎn)，該時(shí)刻此風(fēng)險(xiǎn)對(duì)應(yīng)的一級(jí)指標(biāo)樣本取值集合為FW＝f＝｛FW1＝f1，F(xiàn)W2＝f2，F(xiàn)W3＝f3，預(yù)測在FW＝f 條件下該時(shí)刻瓦斯事故風(fēng)險(xiǎn)的后驗(yàn)概率，以此評(píng)估該時(shí)刻PW＝Y(jié) 的可能性, 其中PW＝Y(jié) 表示瓦斯風(fēng)險(xiǎn)PW發(fā)生。概率模型計(jì)算公式為：

一級(jí)指標(biāo)與其對(duì)應(yīng)多級(jí)指標(biāo)處于同一致因鏈，兩者之間存在一定程度的聯(lián)系。但由于該預(yù)處理階段通過全國歷史事故數(shù)據(jù)學(xué)習(xí)與專家經(jīng)驗(yàn)結(jié)合構(gòu)建貝葉斯網(wǎng)絡(luò)，據(jù)此考察實(shí)際煤礦的兩者之間的映射關(guān)系存在局限性，因此系統(tǒng)將煤礦風(fēng)險(xiǎn)預(yù)測分為2級(jí)，在預(yù)處理階段并不進(jìn)一步探究兩者間的聯(lián)系。

2.2 數(shù)據(jù)處理階段

數(shù)據(jù)處理階段將來自云平臺(tái)的實(shí)時(shí)數(shù)據(jù)和非實(shí)時(shí)數(shù)據(jù)作為輸入，并進(jìn)行過濾、融合及精簡處理。

1）數(shù)據(jù)過濾的目的是剔除不符合約束要求的數(shù)據(jù)，數(shù)據(jù)過濾將剔除監(jiān)測數(shù)據(jù)中因各種原因?qū)е碌娜缑黠@數(shù)據(jù)異常、錯(cuò)誤格式數(shù)據(jù)等錯(cuò)誤數(shù)據(jù)，并通過RS 理論對(duì)缺失數(shù)據(jù)進(jìn)行補(bǔ)完。

2）數(shù)據(jù)精簡的目的是通過RS 理論在數(shù)據(jù)分析的基礎(chǔ)上剔除對(duì)事故風(fēng)險(xiǎn)影響可以忽略的指標(biāo)，從而減少SVM 模型的訓(xùn)練時(shí)間，降低其復(fù)雜程度。

3）數(shù)據(jù)合并的目的是將同一監(jiān)測區(qū)域的多個(gè)同類監(jiān)測數(shù)據(jù)進(jìn)行合并處理，以精確定位風(fēng)險(xiǎn)區(qū)域并減少風(fēng)險(xiǎn)分析中的重復(fù)計(jì)算。合并策略是風(fēng)險(xiǎn)最大估計(jì)策略，即在一個(gè)小范圍的檢測區(qū)域內(nèi)，同類型的監(jiān)測數(shù)據(jù)取其中最壞的結(jié)果，以避免風(fēng)險(xiǎn)漏報(bào)的可能性。具體流程為：取最差值：一個(gè)煤礦區(qū)塊同一項(xiàng)監(jiān)測指標(biāo)的全部傳感器取值集合為Vw={V1，V2，V3，…，Vw}，對(duì)這組數(shù)據(jù)正向化處理并取最小值min（Vw）。歸一化處理：對(duì)全部監(jiān)測數(shù)據(jù)歸一化，采用最大最小標(biāo)準(zhǔn)化法，以進(jìn)行一級(jí)指標(biāo)預(yù)測階段SVM 模型的訓(xùn)練。數(shù)據(jù)合并同時(shí)對(duì)非實(shí)時(shí)數(shù)據(jù)進(jìn)行實(shí)時(shí)歸類，依據(jù)專家意見為非實(shí)時(shí)數(shù)據(jù)設(shè)定的效應(yīng)時(shí)間t，當(dāng)一個(gè)非實(shí)時(shí)數(shù)據(jù)上傳至系統(tǒng)后，在其發(fā)生時(shí)間后的t 時(shí)間內(nèi)，將其視為實(shí)時(shí)數(shù)據(jù)用于態(tài)勢預(yù)測。

2.3 風(fēng)險(xiǎn)預(yù)測階段

一級(jí)指標(biāo)預(yù)測階段實(shí)現(xiàn)對(duì)煤礦風(fēng)險(xiǎn)可能性的預(yù)測。實(shí)際場合中，地區(qū)煤礦的事故數(shù)據(jù)較難獲得，將全國歷史數(shù)據(jù)用于地方煤礦建模適用程度有限；而人為設(shè)定的煤礦風(fēng)險(xiǎn)評(píng)估值較主觀，缺乏公信力。為此，系統(tǒng)將風(fēng)險(xiǎn)預(yù)測中的預(yù)測階段分為2 級(jí)。

1）一級(jí)預(yù)測。通過構(gòu)建一級(jí)指標(biāo)與其對(duì)應(yīng)的多級(jí)指標(biāo)間的SVM 一級(jí)預(yù)測模型對(duì)一級(jí)指標(biāo)進(jìn)行預(yù)測。相比于難以獲得的煤礦事故數(shù)據(jù)，煤礦內(nèi)一級(jí)指標(biāo)和多級(jí)指標(biāo)數(shù)據(jù)在監(jiān)測物聯(lián)網(wǎng)的支持下是可以精確獲取的，以多級(jí)指標(biāo)為輸入自變量，對(duì)應(yīng)下一時(shí)間節(jié)點(diǎn)的一級(jí)指標(biāo)為輸入因變量構(gòu)建一級(jí)指標(biāo)的預(yù)測模型。系統(tǒng)采用RBF 核函數(shù)作為SVM 模型的特征函數(shù)，核函數(shù)參數(shù)g 和錯(cuò)誤懲罰因子C 根據(jù)煤礦具體情況和專家意見確定。

2）二級(jí)預(yù)測。根據(jù)各項(xiàng)風(fēng)險(xiǎn)對(duì)應(yīng)的一級(jí)指標(biāo)預(yù)測結(jié)果通過貝葉斯網(wǎng)絡(luò)風(fēng)險(xiǎn)二級(jí)預(yù)測模型評(píng)估此未來時(shí)間節(jié)點(diǎn)風(fēng)險(xiǎn)產(chǎn)生的可能性。風(fēng)險(xiǎn)可能性預(yù)測后，據(jù)此對(duì)風(fēng)險(xiǎn)進(jìn)行分級(jí)，以評(píng)估其對(duì)煤礦整體安全態(tài)勢的影響，風(fēng)險(xiǎn)分級(jí)見表1。同時(shí)，為了滿足煤礦具體至某一區(qū)塊的報(bào)警需求，需要在風(fēng)險(xiǎn)預(yù)測階段識(shí)別風(fēng)險(xiǎn)發(fā)生區(qū)域，此階段對(duì)煤礦中所有區(qū)塊進(jìn)行單獨(dú)的風(fēng)險(xiǎn)預(yù)測。

表1 安全態(tài)勢影響的風(fēng)險(xiǎn)等級(jí)Table 1 Risk level of security situation impact

2.4 態(tài)勢評(píng)估階段

煤礦安全態(tài)勢就是煤礦在某一時(shí)間節(jié)點(diǎn)內(nèi)各項(xiàng)風(fēng)險(xiǎn)發(fā)生的分布情況及其對(duì)煤礦整體安全的影響程度,與時(shí)間變化、空間分布存在關(guān)系，并且表現(xiàn)為風(fēng)險(xiǎn)的嚴(yán)重程度和影響范圍。對(duì)于某一時(shí)刻煤礦整體的安全態(tài)勢，需要對(duì)此時(shí)刻所有風(fēng)險(xiǎn)的影響進(jìn)行計(jì)算和累加得到反應(yīng)煤礦安全狀態(tài)的態(tài)勢值。

t 時(shí)刻風(fēng)險(xiǎn)N 對(duì)于煤礦整體安全的影響以安全態(tài)勢評(píng)估值SN（t）表示，計(jì)算公式為：

式中：P 為風(fēng)險(xiǎn)可能性；L 為風(fēng)險(xiǎn)等級(jí)；T 為風(fēng)險(xiǎn)告警時(shí)間；S 為風(fēng)險(xiǎn)影響范圍；M 為該風(fēng)險(xiǎn)所有可能存在的區(qū)域的數(shù)量。

T 默認(rèn)值為1，當(dāng)某一區(qū)塊預(yù)測風(fēng)險(xiǎn)存在并預(yù)警后，未及時(shí)解除風(fēng)險(xiǎn)狀態(tài)，至下一時(shí)刻該區(qū)塊該風(fēng)險(xiǎn)仍然存在，則T 值加1。

對(duì)于煤礦整體安全態(tài)勢，需要考慮所有評(píng)估維度的影響。煤礦整體安全態(tài)勢評(píng)估值S（t）計(jì)算為：

式中：N 為主要風(fēng)險(xiǎn)的數(shù)量，N=5；Wi為根據(jù)煤礦實(shí)際情況參考專家經(jīng)驗(yàn)為各項(xiàng)風(fēng)險(xiǎn)設(shè)定的權(quán)重。

對(duì)于單項(xiàng)風(fēng)險(xiǎn)或是煤礦整體的安全態(tài)勢評(píng)估值，當(dāng)評(píng)估值為0 時(shí)，煤礦所有區(qū)塊都不存在風(fēng)險(xiǎn)。當(dāng)評(píng)估值不為0 時(shí)，表明煤礦內(nèi)存在風(fēng)險(xiǎn)，并且值越大，風(fēng)險(xiǎn)可能性越高、風(fēng)險(xiǎn)影響時(shí)間與范圍越廣或是存在的風(fēng)險(xiǎn)越多。

3 實(shí)例分析

安全態(tài)勢預(yù)測建模數(shù)據(jù)來源于山西長治霍爾辛赫試點(diǎn)煤礦物聯(lián)網(wǎng)建設(shè)，鑒于該項(xiàng)目建設(shè)尚未全部完成，無法提供完整的監(jiān)測數(shù)據(jù)，僅以瓦斯風(fēng)險(xiǎn)為例，從瓦斯風(fēng)險(xiǎn)這1 單一維度對(duì)該礦進(jìn)行安全態(tài)勢預(yù)測。樣本數(shù)據(jù)為30 min 上傳1 次的監(jiān)測數(shù)據(jù)，共計(jì)48 組樣本。使用36 組數(shù)據(jù)作為訓(xùn)練集構(gòu)建預(yù)測模型，12 組數(shù)據(jù)作為測試集對(duì)試點(diǎn)煤礦進(jìn)行瓦斯風(fēng)險(xiǎn)預(yù)測以及態(tài)勢評(píng)估結(jié)果見表2。由結(jié)果可知，在監(jiān)測時(shí)間內(nèi)，雖然由于短暫出現(xiàn)火源和瓦斯?jié)舛忍嵘龑?dǎo)致瓦斯風(fēng)險(xiǎn)可能性上升，但提升幅度較小，仍處于無風(fēng)險(xiǎn)范圍內(nèi)。由于監(jiān)測時(shí)間內(nèi)樣本區(qū)塊始終處于無風(fēng)險(xiǎn)狀態(tài)，瓦斯風(fēng)險(xiǎn)安全態(tài)勢評(píng)估值始終為0，證明在監(jiān)測時(shí)間段內(nèi)樣本區(qū)塊瓦斯單項(xiàng)風(fēng)險(xiǎn)的安全態(tài)勢處于良好的狀態(tài)，這與樣本區(qū)塊在該時(shí)間段未產(chǎn)生瓦斯風(fēng)險(xiǎn)這一事實(shí)相符。這一階段內(nèi)區(qū)塊瓦斯風(fēng)險(xiǎn)態(tài)勢評(píng)估值變化曲線是一條取值始終為0 的線段。

表2 瓦斯風(fēng)險(xiǎn)預(yù)測Table 2 Gas risk forecast

4 結(jié) 語

將安全態(tài)勢感知的概念引入煤礦安全領(lǐng)域之中，評(píng)估各時(shí)間節(jié)點(diǎn)內(nèi)煤礦風(fēng)險(xiǎn)發(fā)生的分布情況及其對(duì)煤礦整體安全的影響程度，以實(shí)現(xiàn)在風(fēng)險(xiǎn)產(chǎn)生前預(yù)測風(fēng)險(xiǎn)的到來并提前識(shí)別風(fēng)險(xiǎn)、消除風(fēng)險(xiǎn)的目的。同時(shí)提出了一種基于煤礦風(fēng)險(xiǎn)維度的煤礦整體安全態(tài)勢評(píng)價(jià)方法，以煤礦安全態(tài)勢走向圖的形式為煤礦管理者提供一種極其直觀簡明了解煤礦整體安全情況的方法，并在此基礎(chǔ)上構(gòu)建了集智能感知、風(fēng)險(xiǎn)動(dòng)態(tài)評(píng)估與智能預(yù)警于一體化的礦山安全態(tài)勢感知系統(tǒng)，以實(shí)現(xiàn)煤礦安全管理的定量化、自動(dòng)化、智能化，為管理者進(jìn)行管理決策提供相當(dāng)程度的支持。