李小明，姜 騫，3，袁 強，任奕瑋，3

(1.重慶大學 資源與安全學院，重慶 400044；2.重慶大學 煤礦災害動力學與控制國家重點實驗室，重慶 400044；3.中國煤炭科工集團重慶研究院有限公司，重慶 400037)

我國煤礦煤與瓦斯突出(簡稱“突出”)事故較為嚴重，常常造成非常嚴重的破壞和次生災害，造成生命和財產(chǎn)損失[1，2]。大家普遍接受的一個觀點認為突出的發(fā)生是煤的結(jié)構(gòu)力學特性、地應力和瓦斯壓力相互作用的結(jié)果[3]。《防治煤與瓦斯突出規(guī)定》[4]中，用于預測煤層突出危險性的有甲烷初始擴散速率ΔP、煤的破壞方式、恒定系數(shù)f和煤層瓦斯壓力P；用于預測揭煤面危險性的有瓦斯地質(zhì)信息、綜合指數(shù)(D，K)和鉆屑的解吸指數(shù)(h2或K1)；用于預測煤巷工作面危險的有鉆孔瓦斯初始釋放速率(Qm)和R值指數(shù)法。這些不連續(xù)且基于經(jīng)驗抽樣的傳統(tǒng)煤與瓦斯突出預測方法，不能完全滿足高產(chǎn)高效礦井安全生產(chǎn)的需求[5，6]。鑒于此，人們不斷深入研究以尋求更好的預測方法，隨著研究手段、技術(shù)的不斷豐富，提出了在時間和空間上連續(xù)動態(tài)的預警方法，即：電磁波預警[7，8]、聲發(fā)射預警[9，10]、微震測量預警[11，12]、煤體溫度輻射預警[13]。這些方法只需較少工作量，而且不會占用工作時間而影響生產(chǎn)。通過對以上預測方法的比較分析，表明傳統(tǒng)預測方法能夠較好地判斷突出危險性，但存在與礦井安全生產(chǎn)效率相沖突的問題，新近研發(fā)的部分監(jiān)測預警系統(tǒng)能減輕工作過程對生產(chǎn)的直接影響，但存在覆蓋面窄和對象單一等不足。因此，考慮到安全監(jiān)控系統(tǒng)在我國瓦斯礦井的普及，并且已有學者和研究機構(gòu)利用瓦斯數(shù)據(jù)結(jié)合相關(guān)資料進行瓦斯涌出分析[14]，進而預測煤與瓦斯突出危險性，本文擬構(gòu)建瓦斯量指標、瓦斯解吸指標和瓦斯波動指標相結(jié)合的煤與瓦斯突出辨識綜合性指標體系[15]，并給出指標及臨界值的優(yōu)選方法以提高辨識效率，最后在平舒礦區(qū)進行應用以檢驗其可行性。

1 基于指標考察的突出預警理論

1.1 預警指標體系

1.1.1 瓦斯量指標模型

1)瓦斯涌出動態(tài)表征指標A模型。建立指標A的目的是為了能夠?qū)⒕蜻M面前方煤體可解吸瓦斯量的發(fā)展狀態(tài)或發(fā)展趨勢數(shù)值化，指標值越大，越有可能發(fā)生突出，指標A計算模型如下：

2)T時間內(nèi)瓦斯涌出量QT模型。建立QT指標是為了能夠衡量落煤瓦斯涌出量的大小，其值越大越可能發(fā)生突出危險。掘進一段時間內(nèi)瓦斯涌出量QT指標的計算公式如下：

1.1.2 瓦斯解吸指標模型

1)瓦斯涌出解吸特征指標B模型。建立指標B的依據(jù)是煤體存在不同的瓦斯解吸速率，故其能實現(xiàn)量化掘進落煤瓦斯解吸特征，進而反映落煤物理力學性質(zhì)的目的。B值越大表明落煤解吸越快，越有可能發(fā)生突出危險。參考井下采煤實際情況提出瓦斯涌出解吸特征指標B計算公式：

Bm=Xm，max-Xm，30min

(3)

式中，Bm為第m班次瓦斯涌出解吸指標B；Xm，max為第m班次瓦斯?jié)舛茸畲笾担?%；Xm，30min為第m班次30min內(nèi)瓦斯?jié)舛纫苿幼钚≈担?%。

2)瓦斯解吸量特征指標i模型。煤礦安監(jiān)系統(tǒng)數(shù)據(jù)可以運用到瓦斯解吸指標計算模型建立當中。特征指標i是煤體瓦斯自然解吸曲線的一個擬合數(shù)值，特征值增大那么解吸曲線將會更加陡峭，其具體的數(shù)學計算模型為：

式中，Q為t時間內(nèi)瓦斯解吸量，m3；α為瓦斯自然解吸擬合曲線參數(shù)；i為瓦斯解吸量特征指標參數(shù)i；t為解吸時間，min。

1.1.3 瓦斯波動指標模型

1)瓦斯分源表征指標C模型。瓦斯分源表征是指瓦斯在落煤期間還是非落煤期間涌出的傾向性，而這取決于煤壁瓦斯流通通道是否順暢，因此指標C的建立能夠反映掘進面煤壁瓦斯涌出的能力，進而間接表征煤體前方的地應力。依據(jù)井下采煤的現(xiàn)實情況提出表征指標C計算公式：

2)變異系數(shù)指標V模型的建立。瓦斯涌出波動最大幅度實際上就是瓦斯最大值和最小值之間的差值，但其監(jiān)控結(jié)果往往會由于系統(tǒng)誤差的存在而出現(xiàn)錯誤數(shù)據(jù)。這類數(shù)據(jù)可由系統(tǒng)的濾噪功能實現(xiàn)自動濾除，其余瞬時數(shù)據(jù)造成的變異系數(shù)指標V的誤差需要由模型減少誤差。變異系數(shù)指標V的計算模型為：

1.2 指標辨識效果考察方法

參照突出危險指標評價方法，提出利用常規(guī)預測信息與辨識結(jié)果之間的對比關(guān)系確定指標的辨識“四率”(辨識率η1、辨識準確率η2、辨識誤報率η3、辨識漏報率η4)，從而量化辨識效果。

設ai(ti，si)為井下發(fā)生突出危險的時間位置信息，由于ai(ti，si)的出現(xiàn)往往是不連續(xù)的，并且巷道的掘進或者采煤的進尺是由工作排班表決定的，因此位置s和時間t是存在函數(shù)關(guān)系s=Y(t)，可令ai(ti-3σ，ti+σ)為經(jīng)查驗存有突出風險的區(qū)域，Pj(tj)為辨識指標發(fā)布威脅或者危險信息。

在給定臨界值XLk下根據(jù)Pj與ai在時間軸上的聯(lián)系，可以得出第L個辨識指標XL的辨識“四率”。

式中，n0為辨識總次數(shù)，次；nLk1為預警次數(shù)，次；nLk2為發(fā)布預警時間處于ai區(qū)間的次數(shù)，次；nLk3為發(fā)布預警但不屬于ai區(qū)間的次數(shù)，次；nLk0為獨立的ai區(qū)間個數(shù)，個；nLk4為獨立的ai區(qū)間內(nèi)不存在預警信息的區(qū)間個數(shù)，個。

2 指標優(yōu)選模型和實時預警

2.1 辨識指標及臨界值優(yōu)選

為使辨識效果最佳，辨識“四率”應滿足線性方程R=a2η2-a1η1-a3η3-a4η4，R越大辨識效果越佳，同時η4恒等于0才能安全生產(chǎn)。臨界值不同所對應的辨識指標“四率”也是不同的，故“四率”應是指標XL的函數(shù)，即等價于求R(XL)=(a2+a3)·H(XL)-a3-a1·F(XL)-a4·G(XL)的最大值，當R(XL)滿足下列條件，則可認為XL存在局部極值XLm：

2.2 預警系統(tǒng)構(gòu)建

2.2.1 系統(tǒng)硬件架構(gòu)

KJA突出預警系統(tǒng)運行流程如圖1所示，借助礦井瓦斯監(jiān)控系統(tǒng)來獲取井下瓦斯?jié)舛葦?shù)據(jù)。系統(tǒng)能夠在不同的巷道預設參數(shù)下自動計算目前采掘面基于瓦斯涌出動態(tài)表征的突出危險性指標，并進行全方位的突出危險性分析、評價，最后通過客戶端向工作人員報警，同時向班組長及礦山領導發(fā)送預警短信且能將預警信息發(fā)布到礦山及集團的局域網(wǎng)內(nèi)。

圖1 KJA突出預警系統(tǒng)架構(gòu)

井下瓦斯監(jiān)控系統(tǒng)、監(jiān)控系統(tǒng)服務器、KJA服務器、KJA管理員控制終端、預警信息發(fā)布終端和預警信息短信發(fā)布終端共同組成突出危險快速辨識系統(tǒng)。監(jiān)控系統(tǒng)服務器是礦井監(jiān)控系統(tǒng)的服務器，其作用是儲存井下環(huán)境參數(shù)監(jiān)控數(shù)據(jù)并實時向KJA服務器傳輸采掘工作面的瓦斯、風速等監(jiān)控數(shù)據(jù)。KJA服務器安裝在礦井安全監(jiān)控中心，按設計要求接收、處理數(shù)據(jù)，并將預警信息發(fā)送給用戶終端。KJA管理員控制終端的主要功能是對系統(tǒng)進行維護，如根據(jù)井下巷道施工情況輸入合適的巷道特征參數(shù)。該系統(tǒng)為C/S、B/S雙構(gòu)架，既可以通過C/S模式設置辨識參數(shù)，查詢、訪問辨識結(jié)果、打印辨識表單，也可以通過基于B/S模型的系統(tǒng)網(wǎng)站，在局域網(wǎng)或互聯(lián)網(wǎng)上及時查詢、訪問辨識結(jié)果。

2.2.2 系統(tǒng)軟件設計

系統(tǒng)軟件接收來自監(jiān)控系統(tǒng)服務器的瓦斯信息后，通過預先設計好的算法進行處理后，判斷工作面是否有突出危險，如果有則發(fā)出預警信號。通過專用網(wǎng)絡和短信平臺可以將預警信息傳送至相關(guān)人員或機構(gòu)。

軟件程序基于C#語言和Microsoft SQL Server 2005 數(shù)據(jù)庫開發(fā)，操作系統(tǒng)平臺為Microsoft Windows XP Professional，能實現(xiàn)預警信息在網(wǎng)絡上便捷地發(fā)布接收，并支持遠程網(wǎng)絡在線訪問Web頁面獲取預警信息。該系統(tǒng)具有實時數(shù)據(jù)和各種統(tǒng)計數(shù)據(jù)存儲能力，具有數(shù)據(jù)濾噪功能、實時辨識功能、辨識結(jié)果查詢功能、聯(lián)動辨識功能、辨識報表功能和系統(tǒng)自動診斷功能。

3 現(xiàn)場應用

3.1 工作面實際危險性傳統(tǒng)判識

對平舒礦81111尾巷掘進工作面進尺165m期間進行統(tǒng)計得知，掘進循環(huán)共27個，其中有6個循環(huán)存在異常，見表1。

表1 尾巷異常信息

(9)

式中，B為預測指標的權(quán)系數(shù)矩陣；C為預測指標賦值矩陣；Ci，j表示第i個指標在第j個循環(huán)內(nèi)的值。

由式(9)求得Y值，再根據(jù)式(10)定義的三類范圍并設σ=24h，可得該尾巷危險度劃分結(jié)果，如圖2所示。該尾巷工作面在累計27個循環(huán)中，屬于危險區(qū)間的循環(huán)有1個，用三角形標注；屬于威脅區(qū)間的循環(huán)有5個，用正方形標注；其余正常無危險區(qū)間則用圓形標注。

圖2 尾巷實際危險度

3.2 辨識指標及臨界值確定

由于風流及瓦斯擴散的作用，井巷各處瓦斯?jié)舛扔休^大差異，為能準確反映工作面瓦斯涌出特征需要開展現(xiàn)場試驗確定傳感器的最佳位置。選取一新掘煤巷工作面作為試驗場所，三個傳感器在掘進初期均位于掘進面開口，其中，T1傳感器和T3傳感器隨著掘進面推進而移動并分別維持距掘進面5m和30m的間隔，T2傳感器始終固定在掘進面開口。試驗結(jié)果表明，T1傳感器和T3傳感器的數(shù)據(jù)一致性極高，平均絕對誤差僅為0.02%。據(jù)此可以認為T1傳感器能夠較理想的反應掘進面瓦斯流動。同時，煤礦監(jiān)控系統(tǒng)現(xiàn)有瓦斯傳感器能夠滿足這一位置要求，可以作為預警系統(tǒng)的辨識傳感器。

由實時預警系統(tǒng)KJA從監(jiān)控系統(tǒng)服務器中讀取上文所提尾巷的瓦斯數(shù)據(jù)并通過預置算法得出4個優(yōu)化指標及其臨界值，并對各辨識指標發(fā)布111次辨識信息，其中各指標最優(yōu)辨識臨界值XLk，辨識“四率”以及辨識效果R值選取結(jié)果見表2。

表2 尾巷辨識最優(yōu)值

當A=1.1時，辨識效果R值最大，其辨識結(jié)果如圖3所示，與圖2比較可知，該系統(tǒng)能夠提前近10h發(fā)出預警。將A指標及其臨界值為1.1推廣應用到平舒礦其他防突工作面，該指標僅以10%的辨識率η1，實現(xiàn)高達83.3%的辨識準確率η2，以及16.7%的較低辨識誤報率η3，并且辨識漏報率η4為0，辨識效果R值達到56.6%，取得了良好的應用效果。

圖3 指標A的辨識結(jié)果

研究表明，本文提出的KJA實時預警系統(tǒng)能夠?qū)Ａ康牡V山瓦斯監(jiān)控數(shù)據(jù)進行預警指標優(yōu)選并在線發(fā)布即時預警信息，在平舒礦已取得較好的應用效果。該預警系統(tǒng)的建立能夠為智能礦山建設中煤礦災害透明化、信息化奠定理論與技術(shù)基礎。

4 結(jié) 論

1)從突出機理以及突出事故發(fā)展的角度，建立了反映煤體瓦斯含量的瓦斯量指標、反映煤體解吸速度的解吸類指標和反映瓦斯紊亂程度的波動類指標的綜合預警指標體系，揭示煤與瓦斯突出機理及其早期孕育演化狀態(tài)，可實現(xiàn)煤與瓦斯突出的準確、實時預警。

2)依據(jù)筆者提出的指標效果考察方法，建立以服務端平臺為基礎，客戶端為依托，集報表、短信、網(wǎng)站信息發(fā)布為一體的自動化突出危險性快速辨識預警系統(tǒng)。實現(xiàn)了預警信息分級定向發(fā)送，為智能礦山建設和煤礦災害控制提供信息化基礎。

3)該預警系統(tǒng)在平舒礦上的成功應用，表明其能夠在不干擾正常生產(chǎn)的前提下提前近10h發(fā)布預警信息，綜合預警準確率達83.3%，預警效果出色，實現(xiàn)了局部突出防治及超前預警的研究目標，為突出危險區(qū)域采取區(qū)域綜合防突措施補充了一定的理論依據(jù)。