徐邦耘，李希建

?

層次分析法?突變理論在煤礦煤與瓦斯突出危險性預(yù)測中的應(yīng)用研究*

徐邦耘1, 2, 3，李希建1, 2, 3

(1.貴州大學(xué)礦業(yè)學(xué)院，貴州 貴陽 550025；2.復(fù)雜地質(zhì)礦山開采安全技術(shù)工程中心，貴州 貴陽 550025；3.貴州大學(xué)瓦斯災(zāi)害防治與煤層氣開發(fā)研究所，貴州 貴陽 550025)

煤與瓦斯突出事故是嚴(yán)重影響煤礦安全生產(chǎn)的事故之一。為了提高對煤礦煤與瓦斯突出危險性預(yù)測的準(zhǔn)確度，將層次分析法?突變理論引入煤礦的煤與瓦斯突出危險性預(yù)測中，通過層次分析法優(yōu)化突變理論中各層指標(biāo)的排序，以達(dá)到提高預(yù)測模型精確度的目的。通過對某煤礦3個煤層的實例計算，由其突變級數(shù)預(yù)測該煤層突出危險性的大小，表明基于層次分析法?突變理論預(yù)測模型在煤礦煤與瓦斯突出危險性預(yù)測中是有效的。

煤與瓦斯突出；預(yù)測模型；突變理論；層次分析法

0 引 言

目前，對于煤礦煤與瓦斯突出預(yù)測的方法很多，常用的有鉆屑指標(biāo)法與綜合指標(biāo)法[1]。傳統(tǒng)的預(yù)測方法在時間成本與經(jīng)濟成本上的消耗都比較大，所以越來越多的學(xué)者開始選擇用數(shù)學(xué)建模的方式來進行煤礦煤與瓦斯突出危險性的預(yù)測[2?4]。但是單一的預(yù)測模型用于實際生產(chǎn)中存在諸多問題，目前的預(yù)測模型大多都是基于數(shù)學(xué)方法優(yōu)化后的預(yù)測模型。溫廷新等[5]通過基于灰色關(guān)聯(lián)熵的神經(jīng)網(wǎng)絡(luò)預(yù)測模型進行了煤層與瓦斯預(yù)測，通過灰色關(guān)聯(lián)熵優(yōu)化了神經(jīng)網(wǎng)絡(luò)預(yù)測模型，使得預(yù)測精確度更好；匡亮等[6]運用BP神經(jīng)網(wǎng)絡(luò)預(yù)測法對成貴高鐵7座瓦斯隧道進行了煤與瓦斯突出預(yù)測；郭德勇等[7]在進行礦井瓦斯預(yù)測時，引入了突變理論的方法，通過研究，認(rèn)為突變理論在煤與瓦斯突出預(yù)測中是合理且可行的。

突變理論是根據(jù)現(xiàn)有的模型勢函數(shù)以及分歧方程推導(dǎo)歸一公式，通過歸一公式對底層指標(biāo)按照重要度排序，從而得到頂層目標(biāo)的突變級數(shù)的預(yù)測模型。但是，在實際預(yù)測中，底層指標(biāo)的重要度受主觀影響過大，所以直接應(yīng)用突變理論進行預(yù)測，預(yù)測結(jié)果精確度不夠高，本文通過在突變理論的基礎(chǔ)上引入層次分析法進一步優(yōu)化突變理論，使得預(yù)測結(jié)果準(zhǔn)確度更高。

1 基于AHP的隧道煤與瓦斯突出突變理論預(yù)測模型

1.1 模型構(gòu)建步驟

(1) 構(gòu)建煤與瓦斯突出危險性預(yù)測指標(biāo)體系。根據(jù)煤與瓦斯突出的機理，對影響煤與瓦斯突出的指標(biāo)進行劃分，分解為多個層級，由多個指標(biāo)構(gòu)成多層預(yù)測系統(tǒng)。

(2) 對各層指標(biāo)分別進行重要度排序。根據(jù)所構(gòu)建的多層預(yù)測系統(tǒng)，利用層次分析法進行重要度計算排序。首先需要在每個層級中對其指標(biāo)進行重要度比較，按1-9的賦值構(gòu)造判斷矩陣，其最大特征值max要與對應(yīng)的特征向量滿足以下等式：

所得的max就是指標(biāo)重要度值，從大到小進行排序；然后需要進行一致性檢驗，判斷矩陣是否具有滿意的一致性，否則就需對判斷矩陣進行調(diào)整。

(3) 突出危險性預(yù)測。根據(jù)排序過后的多層預(yù)測系統(tǒng)進行突變理論預(yù)測。

表1 常見的突變數(shù)學(xué)模型勢函數(shù)及其分歧方程

對表1中的突變數(shù)學(xué)模型的分歧點集方程進行變換和推導(dǎo)，得出尖點突變和燕尾突變兩種突變模型的歸一公式如下：

通過歸一公式的計算求出最頂層變量值，該值就為該變量的突變級數(shù)，越靠近1危險性越大。

1.2 指標(biāo)的確定與模型建立

根據(jù)相關(guān)學(xué)者[8-10]的研究，本文將煤礦煤與瓦斯突出的主要因素歸結(jié)于地應(yīng)力、瓦斯參數(shù)、煤的物理特性3類，在這3類中根據(jù)隧道可測量的數(shù)據(jù)分別選取以下7個指標(biāo)來描述上述的3個因素，其中影響地應(yīng)力的指標(biāo)為：煤層埋深、瓦斯壓力以及煤的孔隙率；瓦斯參數(shù)包括：瓦斯含量、瓦斯放散初速度以及瓦斯壓力；煤的物理特性包括：煤的堅固性系數(shù)、煤的孔隙率和煤的焦渣特征。

由于在突變理論的運用中要求每個指標(biāo)有排序，因此通過層次分析法來將每個指標(biāo)進行量化，對比其重要程度，再運用突變理論來預(yù)測穿煤層隧道的煤與瓦斯突出危險性。

在地應(yīng)力下的指標(biāo)排序確定中，根據(jù)文獻總結(jié)與專家打分的方法，確定判斷矩陣見表2。

表2 地應(yīng)力判斷矩陣

注：其中C1代表煤的埋深，C2代表煤的孔隙率，C3代表瓦斯壓力。

經(jīng)計算max=3.053，對應(yīng)的特征向量為[0.825 0.218 0.521]，一致性檢驗符合一致性。

同理在瓦斯參數(shù)的指標(biāo)中，確定的判斷矩陣見 表3。

表3 瓦斯參數(shù)判斷矩陣

注：其中C4代表瓦斯含量，C5代表瓦斯放散初速度。

經(jīng)計算max=3.018，對應(yīng)的特征向量為[0.852 0.488 0.186]，一致性檢驗符合一致性。

同理在煤的物理特性的指標(biāo)中，確定的判斷矩陣見表4。

表4 煤的物理特性判斷矩陣

注：C6代表煤的堅固系數(shù)，C7代表煤的焦渣特征。

經(jīng)計算max=3.073，對應(yīng)的特征向量為[0.394 0.902 0.172]，一致性檢驗符合一致性。

通過計算，準(zhǔn)則層通過一致性檢驗的最大特征值對應(yīng)的特征向量為[0.491 0.707 0.509]。通過計算，總排序為：

C3=0.853>C6=0.459>C1=0.419>C4=0.345>C2= 0.307>C5=0.132>C7=0.088

由此構(gòu)造突變級數(shù)評價指標(biāo)，如圖1所示。

圖1 瓦斯突出危險性預(yù)測評價指標(biāo)體系

注：各層指標(biāo)從左自右進行重要性排序。

2 實例計算

通過上文所介紹的模型，選取某煤礦的3個煤層，分別為A1煤層、A3煤層以及A5煤層，進行實例計算，根據(jù)該礦所提供的相關(guān)參數(shù)，整理出表5所示的7個煤與瓦斯突出危險性預(yù)測指標(biāo)的數(shù)據(jù)，對于指標(biāo)體系中數(shù)值不在0～1之間的指標(biāo)，不能直接進行數(shù)值計算，需進行無量綱轉(zhuǎn)換，將其取值轉(zhuǎn)換為0～1之間的數(shù)量值才能進行突變級數(shù)計算，具體數(shù)據(jù)結(jié)果如表6所示。

表5 桃子埡煤與瓦斯突出危險性預(yù)測指標(biāo)原始數(shù)據(jù)

表6 規(guī)范化后的煤與瓦斯危險性預(yù)測指標(biāo)及突變級數(shù)數(shù)據(jù)

以A1煤層為例進行突變級數(shù)的演算。由圖1可知，地應(yīng)力下的C1、C2、C3指標(biāo)構(gòu)成燕尾突變模型，C1、C2、C3的重要程度排序為C3>C1>C2，采用燕尾突變歸一公式有：

由于C1、C2、C3對煤與瓦斯突出的危險性具有互補性，共同對準(zhǔn)則層評價指標(biāo)起作用，因此，按照“互補”原則，求取平均值有：

同理，在瓦斯參數(shù)準(zhǔn)則層中，C3、C4、C5構(gòu)成燕尾突變模型，其重要度排序為C3>C4>C5。由燕尾歸一公式有，x3=0.959, x4=0.913，x5=0.966，按照互補原則，求其平均值，得2(瓦斯參數(shù))=0.946。

同理，煤的物理特性準(zhǔn)則層也滿足燕尾突變模型，C2、C6、C7的重要度排序為C6>C2>C7。

通過計算得，x6=1.000，x2=0.990，x7=0.905，3(煤的物理特性)=0.965。

對于準(zhǔn)則層目標(biāo)B1、B2、B3構(gòu)成燕尾突變，其重要性排序為B2>B3>B1根據(jù)燕尾突變模型歸一公式有：x2=0.973，x3=0.988，x1=0.996，然后根據(jù)互補原則，取平均值，得到A1=0.986。即A1煤層的突變級數(shù)值為0.986。同理可計算其他兩個煤層的突變級數(shù)值。

根據(jù)3個煤層的突變級數(shù)計算結(jié)果顯示，其中A5煤層的突變級數(shù)值最大，其突出危險性最大，A3煤層的突出危險性最小，但是3個煤層的突變級數(shù)值都比較接近1，所以計算結(jié)果顯示3個煤層都具有突出危險性。根據(jù)該礦的開采事實，A5煤層的突出風(fēng)險為最高，這與計算結(jié)果相印證。

3 結(jié) 論

(1) 在綜合分析煤與瓦斯突出影響因素的基礎(chǔ)上，基于瓦斯突出綜合作用假說機理，本文提出了突變理論預(yù)測隧道煤與瓦斯突出危險性的方法，該方法綜合反映了眾多復(fù)雜的影響因素，具有更高效的判別效率。

(2) 本文構(gòu)建了基于層次分析法?突變理論預(yù)測模型，根據(jù)對某煤礦的實例預(yù)測計算，計算結(jié)果與實際情況相符合，驗證了基于層次分析法的突變理論模型應(yīng)用到煤礦煤與瓦斯突出危險性預(yù)測是可行的。

(3) 通過基于層次分析的突變理論模型對某煤礦進行煤與瓦斯突出危險性預(yù)判，并與實際情況相對比，結(jié)果顯示，該煤礦的3個煤層都具有突出危險性，其中A5煤層的突出危險性最高。

[1] 國家安全生產(chǎn)監(jiān)督管理局,國家煤礦安全監(jiān)察局.煤礦安全規(guī)程.(2014 版)[M].北京:煤炭工業(yè)出版社,2014.

[2] 李 苗,張紅軍,陳唐軍,劉夏臨,王紀(jì)鵬.灰色－突變理論在尾礦壩失穩(wěn)預(yù)測中的應(yīng)用[J].礦業(yè)研究與開發(fā),2014,34(5):66?68.

[3] 于海云,楊 力.基于灰色關(guān)聯(lián)與SVM 模型的煤礦通風(fēng)系統(tǒng)評價方法研究[J].煤礦安全,2013,44(1):181?184.

[4] 夏 龍,郝兵元,錢翰飛.模糊數(shù)學(xué)綜合評價在矸石山邊坡穩(wěn)定性評價中的應(yīng)用[J].礦業(yè)研究與開發(fā),2017,37(10):33?37.

[5] 溫廷新,于鳳娥,邵良杉.基于灰色關(guān)聯(lián)熵的煤與瓦斯突出概率神經(jīng)網(wǎng)絡(luò)預(yù)測模型[J].計算機應(yīng)用研究,2018,35(11):96?100.

[6] 匡 亮,趙萬強,喻 渝.BP神經(jīng)網(wǎng)絡(luò)法預(yù)測隧道瓦斯突出的模型與實例[J].鐵道工程學(xué)報,2018,233(2):56?61.

[7] 郭德勇,李佳乃,王彥凱.基于黏滑失穩(wěn)與突變理論的煤與瓦斯突出預(yù)警模型[J].北京科技大學(xué)學(xué)報,2013(11):1407?1412.

[8] 吳應(yīng)明.華鎣山隧道有害氣體監(jiān)測與綜合治理技術(shù)[J].現(xiàn)代隧道技術(shù),2003,40(4):68?73.

[9] 袁森林,劉 旭,周建中.瓦斯隧道工程勘查技術(shù)探討[J].四川地質(zhì)學(xué)報,2009,29(1):50?51.

[10] 袁真秀,張羽軍,王崇良.新寨隧道煤與瓦斯突出危險性評價及揭煤施工設(shè)計[J].四川工業(yè)學(xué)院學(xué)報,2003,22(3):88?91.

國家自然科學(xué)基金面上項目(51874107).

(2018?10?22)

徐邦耘(1995－)，男，貴州六盤水人，在讀碩士生，主要從事安全科學(xué)技術(shù)相關(guān)研究，Email：790606745@ qq.com。