黃繼廣，馬漢鵬，2，范春姣，姚絨絨

(1.華北科技學(xué)院，北京 065201；2.中國煤礦安全技術(shù)培訓(xùn)中心，北京 065201)

0 引言

現(xiàn)今全球煤炭資源儲(chǔ)量合計(jì)1.14萬億t，其中美國以2 516億t居首位，占比22%；而中國以2 440億t居第二位，占全球21.4%。近些年來，隨著我國煤炭產(chǎn)量的不斷增加，煤炭產(chǎn)量從2009年29.1億t增加到2018年36.8億t，增長(zhǎng)20.9%。我國現(xiàn)有煤礦的開采方式多為井工開采，但是井下的條件比較特殊，生產(chǎn)的過程繁瑣、復(fù)雜，另外還要受到各類自然災(zāi)害的威脅，因此，煤礦事故發(fā)生的概率比較大[1]。通過持續(xù)多年的綜合整治，煤礦安全生產(chǎn)形勢(shì)持續(xù)穩(wěn)定向好，重特大事故多發(fā)勢(shì)頭得到有效遏制。雖然安全形勢(shì)逐年好轉(zhuǎn)，但是與美國等發(fā)達(dá)國家相比仍有很大的差距，2004年美國的百萬噸死亡率為0.027，而我國2017年煤礦百萬噸死亡率為0.106。因此，我國的煤礦安全形勢(shì)依然嚴(yán)峻。

近年來，國內(nèi)外很多學(xué)者都對(duì)煤礦生產(chǎn)事故展開了統(tǒng)計(jì)研究。1978年，歐洲共同體所屬礦業(yè)安全與衛(wèi)生委員會(huì)第一次發(fā)表了關(guān)于煤礦事故統(tǒng)計(jì)的報(bào)告，統(tǒng)計(jì)了各成員國的事故傷亡數(shù)字，分析了事故傷亡原因，并對(duì)改進(jìn)安全工作提出了相關(guān)建議[3]；程磊[2]通過對(duì)煤礦百萬噸死亡率、事故區(qū)域和類型、月份等6個(gè)維度對(duì)煤礦事故發(fā)生規(guī)律進(jìn)行分析研究，并提出相關(guān)對(duì)策措施；朱云飛等[4]收集2000—2016年全國煤礦重特大事故數(shù)據(jù)資料，從時(shí)間、地域、經(jīng)濟(jì)、事故類別四個(gè)維度統(tǒng)計(jì)分析，得出預(yù)防重特大事故困難的原因；諸利一等[5]對(duì)近年全國煤礦事故統(tǒng)計(jì)數(shù)據(jù)，從事故級(jí)別、類型、地點(diǎn)3個(gè)維度開展研究，得出貴州煤礦形勢(shì)嚴(yán)峻及頂板事故發(fā)生率最高等規(guī)律，并構(gòu)建了重大、特重大事故與瓦斯、火災(zāi)的關(guān)系模型。通過以上文獻(xiàn)，可以從中了解到當(dāng)前大多學(xué)者主要研究方向基本只針對(duì)煤礦事故統(tǒng)計(jì)分析，并提出相關(guān)措施和建議。為此，依據(jù)2009—2018年我國煤礦生產(chǎn)安全事故統(tǒng)計(jì)資料，分析我國近10年煤礦事故發(fā)生的規(guī)律；并且通過運(yùn)用數(shù)據(jù)分析軟件SPSS 25.0構(gòu)建Brown線性趨勢(shì)模型，對(duì)2019—2021年我國百萬噸死亡率進(jìn)行預(yù)測(cè)分析，以此為政府和企業(yè)對(duì)煤礦安全管理及相關(guān)政策制定提供理論依據(jù)。

1 2009—2018年煤礦生產(chǎn)事故數(shù)據(jù)統(tǒng)計(jì)

1.1 事故統(tǒng)計(jì)

為保證分析的真實(shí)準(zhǔn)確，數(shù)據(jù)收集于國家或地方煤礦安全監(jiān)察局網(wǎng)站、公開發(fā)表的文獻(xiàn)資料[5-7]。得到2009—2018年我國煤礦事故死亡人數(shù)、事故起數(shù)和百萬噸死亡率統(tǒng)計(jì)表，見表1。

表1 2009—2018年我國煤礦事故死亡人數(shù)、事故起數(shù)和百萬噸死亡率統(tǒng)計(jì)表

根據(jù)表1相關(guān)數(shù)據(jù)，得知2009—2018年我國煤礦生產(chǎn)安全事故起數(shù)共發(fā)生7 107起，死亡人數(shù)共12 205人；平均每年約發(fā)生711起，死亡人數(shù)約1 221人。2009年事故起數(shù)和死亡人數(shù)最多，分別為1 616起、2 631人，2009—2018年間，事故起數(shù)和死亡人數(shù)都有大幅下降。煤礦百萬噸死亡率從2009年的0.892下降到2018年的0.093，同比下降了89.6%，并且其首次降到0.1以下。

1.2 事故發(fā)展趨勢(shì)

根據(jù)表1相關(guān)數(shù)據(jù)，繪制2009—2018年我國煤礦生產(chǎn)安全事故發(fā)展趨勢(shì)圖，如圖1所示。從圖1可以直觀看出2009—2018年我國煤礦事故起數(shù)和死亡人數(shù)一直下降，我國煤礦生產(chǎn)安全形勢(shì)持續(xù)向好。由此可見，近年來我國相關(guān)部門和煤礦企業(yè)采取的一系列預(yù)防煤礦事故發(fā)生的舉措已得到實(shí)施，并取得了很好的效果。

圖1 2009—2018年我國煤礦生產(chǎn)安全事故發(fā)展趨勢(shì)圖

2 2009—2018年我國煤礦生產(chǎn)安全事故

從煤礦事故發(fā)生級(jí)別、事故類型2個(gè)角度對(duì)我國近10年煤礦生產(chǎn)安全事故統(tǒng)計(jì)分析，并對(duì)其事故發(fā)生規(guī)律進(jìn)行研究。受事故資料詳細(xì)度限制，部分內(nèi)容屬不完全統(tǒng)計(jì)，但由于事故總體數(shù)量龐大和抽樣的隨機(jī)性，并不影響統(tǒng)計(jì)結(jié)果的客觀性[4]。

2.1 事故發(fā)生級(jí)別

我國2007年頒布的《生產(chǎn)安全事故報(bào)告和調(diào)查處理?xiàng)l例》規(guī)定，按事故造成的人員傷亡情況將事故劃分為一般事故(3人以下)、較大事故(3～9人)、重大事故(10～29人)以及特大事故(30人及以上)4個(gè)等級(jí)[8]。2009—2018年我國煤礦發(fā)生的各大事故不完全統(tǒng)計(jì)概況見表2，2009—2018年我國煤礦不同級(jí)別安全事故占比概況如圖2所示。

表2 2009—2018年我國煤礦發(fā)生的各大事故起數(shù)和死亡人數(shù)(不完全統(tǒng)計(jì))

圖2 2009—2018年我國煤礦各大生產(chǎn)安全事故起數(shù)餅狀圖

由此看出，2009—2018年間，我國煤礦生產(chǎn)安全事故發(fā)生次數(shù)最多的是一般事故(占事故總起數(shù)的93.35%)，其次是較大事故(占事故總起數(shù)的5.33%)，再次是重大事故(占事故總起數(shù)的1.17%)。從表2可以看出，2009—2018年較大事故發(fā)生379起，其中2016—2017年出現(xiàn)增加趨勢(shì)；重大事故發(fā)生83起，其中2009—2011年、2012—2013年、2016—2017年都有不同程度的增加；特別重大事故發(fā)生12起，其中2014年、2015年、2017年及2018年都不曾發(fā)生特別重大事故，而2016年發(fā)生2起特大事故，造成65人死亡的悲劇。從中看來，2016年煤礦安全監(jiān)管部門在安全監(jiān)管監(jiān)察力度上有所松懈，思想上準(zhǔn)備不充分；現(xiàn)階段，我國發(fā)展對(duì)煤礦的依賴度較高，許多礦井都提高了煤炭資源采出率；我國有相當(dāng)數(shù)量的礦井的開采條件復(fù)雜，故完全消除重特大事故十分困難[4]。重大事故起數(shù)和死亡人數(shù)總體呈現(xiàn)下降趨勢(shì)，2009年和2018年重大事故起數(shù)和死亡人數(shù)相比，分別下降到77.8%和70.9%。在煤礦生產(chǎn)安全事故中，重大事故、較大事故和一般事故發(fā)生率均在逐年降低，煤礦安全生產(chǎn)形勢(shì)得到進(jìn)一步好轉(zhuǎn)，但是特大事故時(shí)有發(fā)生，因此我國的煤礦安全生產(chǎn)形勢(shì)依然嚴(yán)峻。所以在今后的煤礦事故預(yù)防工作中，對(duì)重大事故的預(yù)防要給予高度重視，對(duì)構(gòu)建大型礦井事故的應(yīng)急救援體系也要高度重視，并盡可能減少事故發(fā)生時(shí)造成的人員傷亡及財(cái)產(chǎn)損失[8]。

2.2 事故發(fā)生類型

事故起數(shù)和死亡人數(shù)統(tǒng)計(jì)：通過查閱國家煤監(jiān)局網(wǎng)站[9]和文獻(xiàn)資料[7，10]，不完全統(tǒng)計(jì)了2009—2018年我國煤礦各類死亡事故和事故起數(shù)，見表3。

表3 2009—2018年我國煤礦各類事故起數(shù)和死亡人數(shù)不完全統(tǒng)計(jì)表

事故起數(shù)占比分析：從表3和圖3可以看出，2009—2018年全國發(fā)生的各類型死亡事故起數(shù)中，頂板事故、運(yùn)輸事故、瓦斯事故和其他事故發(fā)生比較多，分別是3 255起(占事故總起數(shù)的45.8%)、1 301起(占事故總起數(shù)18.3%)、722起(占事故總起數(shù)10.2%)、904起(占事故總起數(shù)12.7%)，而機(jī)電事故、水害事故、火災(zāi)事故及爆破事故分別是448起、234起、36起和189起。

圖3 2009—2018年煤礦各類事故起數(shù)所占比例

死亡人數(shù)占比分析：從表3和圖4可以看出，2009—2018年全國發(fā)生的各類型事故死亡人數(shù)中，頂板、瓦斯、運(yùn)輸與其他事故類型死亡人數(shù)相對(duì)較多，分別是3 952人(占總死亡人數(shù)的32.4%)、3 433人(占總死亡人數(shù)的28.1%)、1 530人(占總死亡人數(shù)的12.5%)、1 176人(占總死亡人數(shù)的 9.6%)。

圖4 2009—2018年煤礦各類事故死亡人數(shù)所占比例

事故原因分析：由此可以得出，在2009—2018年發(fā)生的8類煤礦事故類型中，頂板事故發(fā)生率最高，累計(jì)造成死亡的人數(shù)也是最多，原因大都是由于工作面支護(hù)強(qiáng)度或者支護(hù)質(zhì)量差，初撐力不符合要求，造成頂板離層，在控頂區(qū)內(nèi)發(fā)生頂板垮塌和冒落或是支護(hù)不及時(shí)，造成頂板懸空時(shí)間較長(zhǎng)，引起頂板垮落；在維檢修巷道時(shí)安全退路不暢，安全確認(rèn)不到位，冒險(xiǎn)作業(yè)等造成頂板事故的發(fā)生[11]。其次是運(yùn)輸事故、其他事故和瓦斯事故的發(fā)生起數(shù)占總事故起數(shù)比例較大，但在這三種事故類型中，瓦斯事故造成的傷害程度更高，死亡人數(shù)更多。其中造成瓦斯事故發(fā)生的原因是存在個(gè)別礦井抽采不達(dá)標(biāo)，2個(gè)“四位一體”措施落實(shí)不到位；通風(fēng)系統(tǒng)設(shè)計(jì)不合理，局部通風(fēng)管理紊亂，甚至存在長(zhǎng)期無風(fēng)作業(yè)；有的低瓦斯礦井忽視安全管理，存在監(jiān)測(cè)監(jiān)控系統(tǒng)運(yùn)行不正常等問題。自2017年國家煤監(jiān)局開展煤礦安全生產(chǎn)標(biāo)準(zhǔn)化以來，各煤礦企業(yè)也成立相關(guān)部門進(jìn)行專項(xiàng)管理，煤礦事故發(fā)生起數(shù)和死亡人數(shù)逐年減少，主要發(fā)生的事故類型得到有效控制，尤其是瓦斯事故起數(shù)從2009年發(fā)生157起下降到2018年的11起，減少了93%。因此，從總體上看瓦斯事故起數(shù)呈現(xiàn)持續(xù)降低的趨勢(shì)，這也足以說明自國家提出的“通風(fēng)可靠，抽采達(dá)標(biāo)，監(jiān)控有效，管理到位”的瓦斯綜合治理工作體系同國家先進(jìn)的科學(xué)技術(shù)水平相融合，基本控制了煤礦瓦斯事故的發(fā)生。

3 我國煤礦百萬噸死亡率預(yù)測(cè)分析

3.1 指數(shù)平滑預(yù)測(cè)法

指數(shù)平滑法介紹：指數(shù)平滑法又稱指數(shù)加權(quán)平均法，由美國學(xué)者BROWN R G提出的，常用的指數(shù)平滑法包括一次指數(shù)平滑、二次指數(shù)平滑和三次指數(shù)平滑[12]。指數(shù)平滑法是在移動(dòng)平均法基礎(chǔ)上發(fā)展起來的一種時(shí)間序列分析預(yù)測(cè)法，它是通過計(jì)算指數(shù)平滑值，配合一定的時(shí)間序列預(yù)測(cè)模型對(duì)現(xiàn)象的未來進(jìn)行預(yù)測(cè)?；舅枷胧菍?duì)原有的數(shù)據(jù)進(jìn)行處理，處理過后的數(shù)據(jù)稱為“平滑值”，然后根據(jù)平滑值構(gòu)造出預(yù)測(cè)模型，用來預(yù)測(cè)未來的數(shù)據(jù)。原理是對(duì)時(shí)間序列進(jìn)行修勻，不過不是求算數(shù)平均值，而是注重時(shí)間序列的長(zhǎng)期數(shù)值對(duì)未來預(yù)測(cè)值的影響，即對(duì)時(shí)間序列的各個(gè)數(shù)據(jù)進(jìn)行加權(quán)平均，時(shí)間越近的數(shù)據(jù)，其權(quán)值越大。

指數(shù)平滑法應(yīng)用：指數(shù)平滑法的基本思想是跟移動(dòng)加權(quán)平均相關(guān)聯(lián)的，而且平滑方法的選擇關(guān)鍵在于平滑常數(shù)α如何確定。平滑常數(shù)α是指移動(dòng)加權(quán)平均的權(quán)數(shù)，體現(xiàn)不同時(shí)期數(shù)據(jù)與現(xiàn)實(shí)聯(lián)系的緊密程度。因此，為達(dá)到最好的預(yù)測(cè)效果，必須合理確定α值。

③模型檢驗(yàn)。常用統(tǒng)計(jì)量指標(biāo)有可決系數(shù)R2和平均絕對(duì)百分誤差MAPE。其中

SPSS軟件的優(yōu)勢(shì)：SPSS 25.0軟件中有利用Brown線性趨勢(shì)指數(shù)平滑模型擬合數(shù)據(jù)的功能，為使用該方法計(jì)算預(yù)測(cè)結(jié)果提供了方便。

3.2 Brown線性趨勢(shì)模型的構(gòu)建及應(yīng)用

最佳擬合模型確定：軟件SPSS 25.0的時(shí)間序列建模器中的專家建模器，是可以自動(dòng)查找時(shí)間序列的最佳擬合模型。布朗線性趨勢(shì)模型，該模型適用于具有線性趨勢(shì)且沒有季節(jié)性的序列，其平滑參數(shù)是水平和趨勢(shì)，并假定二者等同。根據(jù)圖5及SPSS 25.0軟件分析，煤礦百萬噸死亡率的最佳擬合模型為Brown線性趨勢(shì)模型。

圖5 2009—2018年煤礦百萬噸死亡率趨勢(shì)圖

表4 Brown線性趨勢(shì)模型的擬合結(jié)果

表5 2009—2018年煤礦百萬噸死亡率實(shí)際值與預(yù)測(cè)值的對(duì)比表

圖6 序列圖

2019—2021年煤礦百萬噸死亡率預(yù)測(cè)：運(yùn)用Brown線性趨勢(shì)模型對(duì)2019—2021年煤礦百萬噸死亡率進(jìn)行預(yù)測(cè)，得到預(yù)測(cè)值見表6。

表6 2019—2021年煤礦百萬噸死亡率預(yù)測(cè)值

4 結(jié)論

(1)2009—2018年我國煤礦發(fā)生的特大事故、重大事故、較大事故及一般事故的事故起數(shù)和死亡人數(shù)，總體上是持續(xù)下降的。但較大事故和特大事故分別在2014—2015年、2015—2016年，這兩個(gè)階段出現(xiàn)死亡人數(shù)上升的過程，尤其是2016年發(fā)生的2起特大事故造成了65人死亡，給國家和人民造成了巨大的損失，這也是政府和企業(yè)在現(xiàn)階段需要重視，采取什么樣的思路制定完善的安全管理制度，防止重特大事故的發(fā)生。

(2)2009—2018年我國煤礦在事故類型上造成的死亡人數(shù)相對(duì)較多的主要集中在頂板和瓦斯事故類型上，分別達(dá)到了32.4%和28.1%；在發(fā)生的事故起數(shù)上，頂板和運(yùn)輸事故類型發(fā)生的相對(duì)較多，分別達(dá)到了45.8%和18.3%，而瓦斯事故類型僅僅占總事故起數(shù)的10.2%。因此企業(yè)可以通過加大對(duì)煤礦從業(yè)人員的安全教育培訓(xùn)力度，防治“三違”的發(fā)生；嚴(yán)格根據(jù)規(guī)程要求進(jìn)行防治瓦斯工作，對(duì)于瓦斯事故頻發(fā)地區(qū)需高度重視，構(gòu)建“四位一體，五步配套”防突體系。政府可以加強(qiáng)監(jiān)管力度，嚴(yán)格高瓦斯礦井的準(zhǔn)入和退出，強(qiáng)化瓦斯抽采達(dá)標(biāo)檢查。

(3)基于SPSS 25.0軟件對(duì)2009—2018年我國煤礦百萬噸死亡率統(tǒng)計(jì)數(shù)據(jù)建模預(yù)測(cè)分析，構(gòu)建煤礦百萬噸死亡率的最佳擬合模型-Brown線性趨勢(shì)模型，并對(duì)2019—2021年我國煤礦百萬噸死亡率進(jìn)行預(yù)測(cè)，得到我國煤礦百萬噸死亡率未來三年呈現(xiàn)下降趨勢(shì)，說明我國煤礦安全系數(shù)越來越高，比較符合現(xiàn)在整個(gè)煤礦行業(yè)的安全形勢(shì)。通過2009—2018年我國煤礦百萬噸死亡率實(shí)際值和預(yù)測(cè)值的對(duì)比驗(yàn)證來看，平均相對(duì)誤差達(dá)到7.88%，這是由于預(yù)測(cè)值只是根據(jù)不完全數(shù)據(jù)統(tǒng)計(jì)預(yù)測(cè)的，因此實(shí)際值與預(yù)測(cè)值可能會(huì)有偏差。

(4)因?yàn)槊旱V安全事故具有一定的隨機(jī)性，文中只是單一的對(duì)近10年的煤礦百萬噸死亡率進(jìn)行二次指數(shù)平滑預(yù)測(cè)，得到2019—2021年的預(yù)測(cè)值，并不一定能達(dá)到理想效果。所以可以通過增加樣本數(shù)量和更新樣本數(shù)據(jù)，并結(jié)合多種預(yù)測(cè)方法對(duì)其綜合分析運(yùn)用，來建立新的預(yù)測(cè)模型，以保證未來煤礦百萬噸死亡率的預(yù)測(cè)值更科學(xué)可靠，進(jìn)而可以為煤礦企業(yè)的安全管理工作提供有力的理論依據(jù)。