侯冠慧，張英梅，魏紅霞

（太原理工大學(xué)a.工程訓(xùn)練中心；b.電力工程學(xué)院；c.測控技術(shù)研究所，太原 030024)

礦井突水是煤礦最頻繁發(fā)生的災(zāi)害。尋求煤礦井下突水發(fā)生原理，尋找主動防治技術(shù)是井下工作的一項重要內(nèi)容。造成礦井突水的原因眾多，且一般均為各種因素共同影響的結(jié)果。然而，這種共同作用很難用精確的數(shù)學(xué)語言來描述［1］。

目前，地質(zhì)環(huán)境評價分析、數(shù)學(xué)與計算機(jī)模型評價分析、物探勘探評價分析、水化學(xué)離子及同位素評價分析是國內(nèi)外常用的四種分析方法。但礦井由于評價體系方法使用單一，或設(shè)備不夠先進(jìn)等多方面原因，使得挖掘過程中缺乏有效的動態(tài)監(jiān)測手段。筆者提出結(jié)合地質(zhì)勘探評價分析、自動化檢測以及多傳感器數(shù)據(jù)融合三種技術(shù)融合用于井下水災(zāi)害防治方面，重點(diǎn)突出了在生產(chǎn)過程中的動態(tài)監(jiān)測［2-3］。

1 工作面突水動態(tài)評價模型的設(shè)計

礦井突水是指在特定的水文地質(zhì)環(huán)境條件下，受水壓作用的隔水巖層若同時受到強(qiáng)大外力影響，導(dǎo)致其發(fā)生嚴(yán)重的破壞，從而使得大量承壓水涌入礦井。同時，煤礦中不同的地質(zhì)環(huán)境、巖石性質(zhì)、地應(yīng)力的大小與分布和不同的斷裂構(gòu)造以及破碎帶的發(fā)育情況等都有可能誘發(fā)不同機(jī)理的礦井突水。因此，礦井突水是由多種因素共同作用造成的，這種相互作用并不是直接的線性關(guān)系。礦井突水的評價分析客觀上要求用一種整合多種因素的動態(tài)的非線性處理方法。

人工神經(jīng)網(wǎng)絡(luò)具有非線性自適應(yīng)性信息處理能力，其主要優(yōu)點(diǎn)是能夠自適應(yīng)樣本數(shù)據(jù)，當(dāng)數(shù)據(jù)中有噪聲、形變和非線性時，它能夠正常地工作，很容易繼承現(xiàn)有的領(lǐng)域知識，能夠處理來自多個資源和決策系統(tǒng)的數(shù)據(jù)；提供簡單工具進(jìn)行自動特征選取，產(chǎn)生有用的數(shù)據(jù)表示，可作為專家系統(tǒng)的預(yù)處理器?？梢酝ㄟ^把已獲取的突水因素與預(yù)報結(jié)果等數(shù)據(jù)按照數(shù)學(xué)映射理論，看作輸入到輸出的映射，從而分析它們之間的規(guī)律。最終根據(jù)得到映射關(guān)系用新的輸入數(shù)據(jù)來推算輸出結(jié)果［4］。人工神經(jīng)網(wǎng)絡(luò)技術(shù)在非線性問題或非結(jié)構(gòu)性問題解決方面的獨(dú)特優(yōu)勢在突水預(yù)報中得到極大地應(yīng)用。只要存在足夠的可供學(xué)習(xí)的突水?dāng)?shù)據(jù)，那么神經(jīng)網(wǎng)絡(luò)便可建立起它們之間的映射關(guān)系，進(jìn)而用于礦井突水預(yù)報［5］。

本文進(jìn)一步提出將D-S證據(jù)理論與BP神經(jīng)網(wǎng)絡(luò)相結(jié)合的融合方法應(yīng)用于井下突水評價模型中。一方面利用BP神經(jīng)網(wǎng)絡(luò)為D-S證據(jù)理論分配基本可信度，另一方面在一定程度上改善了單一的神經(jīng)網(wǎng)絡(luò)輸出結(jié)果不穩(wěn)定的現(xiàn)象。信息融合處理與系統(tǒng)評價框圖如圖1所示［6］。

1.1 影響工作面突水的主要因素

要建立突水系統(tǒng)評價模型，需全面分析與礦井工作面突水有關(guān)的水文、地質(zhì)環(huán)境、采礦影響等諸多因素。經(jīng)與礦區(qū)礦井水文地質(zhì)專家討論，提出造成礦井工作面水害的主要因素：頂板和底板的壓力、底板含水層的水位、巷道水的水位和掘進(jìn)工作面前方的水量。

圖1 礦井工作面動態(tài)評價系統(tǒng)框圖

1.1.1 傳感器的布置及參數(shù)測量

提出了根據(jù)現(xiàn)場地質(zhì)結(jié)構(gòu)條件，沿工作面優(yōu)化布置傳感器，從而進(jìn)行水情動態(tài)多參數(shù)測量。其中，頂板、底板布置壓力傳感器，巷道內(nèi)布置液位傳感器。采礦工作面周圍空間含水構(gòu)造的超前探測，采用目前比較先進(jìn)的瞬變電磁法。礦井瞬變電磁法是目前最先進(jìn)的探水方法之一，適用于斷層及裂隙帶富水性評價陷落柱探測、煤層頂?shù)装搴羲畬觿澐值?。在進(jìn)行煤礦工作面開采之前用此方法進(jìn)行探測，可以判斷地下地質(zhì)體的電性、規(guī)模、產(chǎn)狀等。

上文所提到的測量因素便成了一個5維模型（本模型選定突水影響因子5個，即N＝5)的輸入矢量，而模型的輸出則為突水的安全等級S＝｛S1，S2，S3，S4｝。從左到右分別代表四種安全等級，很安全、較安全、一般、危險。神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)如圖2，隱含層層數(shù)、隱含層的神經(jīng)元數(shù)量可依照具體實際需要選取。

圖2 BP網(wǎng)絡(luò)模型

1.2 樣本庫的建立

根據(jù)突水影響因素建立模型樣本數(shù)據(jù)庫。如含水層水位（鉆孔得到)、頂板壓力、底板壓力等煤礦已有數(shù)據(jù)，可以直接導(dǎo)入數(shù)據(jù)庫；而巷道水及掘進(jìn)面前方水量，要不斷測量，積累數(shù)據(jù)，逐步建立完善的數(shù)據(jù)樣本及數(shù)據(jù)庫。影響因數(shù)便是數(shù)據(jù)庫中字段，每一個安全生產(chǎn)實例或礦井發(fā)生突水案例、都是數(shù)據(jù)庫中的一條記錄［8］。

1.3 BP網(wǎng)絡(luò)的建立

1.3.1 生成網(wǎng)絡(luò)

利用函數(shù)newff（)可建立一個BP神經(jīng)網(wǎng)絡(luò)。其調(diào)用格式為：

式中：Xr是由R維的輸入樣本最大最小值構(gòu)成的R×2維矩陣；［S1S2…SN1］為各層的神經(jīng)元個數(shù)；｛TF1TF2…TFN1｝為各層神經(jīng)元的傳遞函數(shù)；BLF表示反向權(quán)值學(xué)習(xí)函數(shù)；PF表示性能數(shù)；BTF表示網(wǎng)絡(luò)的反向訓(xùn)練用函數(shù)的名稱［11］。

1.3.2 權(quán)值初始化

BP神經(jīng)網(wǎng)絡(luò)在訓(xùn)練之前必須對權(quán)值和偏值進(jìn)行初始化，在已知網(wǎng)絡(luò)輸入向量范圍，各層神經(jīng)元個數(shù)及傳遞函數(shù)，就可以用init（)函數(shù)進(jìn)行初始化，得到每層網(wǎng)絡(luò)的權(quán)值和偏值。

1.3.3 仿真網(wǎng)絡(luò)

給定BP網(wǎng)絡(luò)結(jié)構(gòu)和輸入變量，每一層都從前一層得到輸入數(shù)據(jù)，就可以應(yīng)用simuff（)函數(shù)仿真相應(yīng)的網(wǎng)絡(luò)輸出。

1.4 D-S證據(jù)理論

D-S算法是由數(shù)學(xué)家A.P.Dempster提出，并由他的學(xué)生G.Shafer做進(jìn)一步發(fā)展，形成的基于“證據(jù)”和“組合”來處理不確定性推理問題的數(shù)學(xué)方法。作為一 種不精確推理方法，證據(jù)理論的主要特點(diǎn)是：滿足比貝葉斯概率論更弱的條件；具有直接表達(dá) “不確定”和 “不知道”的能力。

1.4.1 應(yīng)用D-S證據(jù)理論的依據(jù)

將D-S證據(jù)理論應(yīng)用到水害評價模型中，D-S證據(jù)理論是將每一傳感器對突水狀態(tài)的影響轉(zhuǎn)化為概率值來進(jìn)行融合的，D-S證據(jù)理論n個證據(jù)源的合成公式如下［10］：

式中，K表示各個證據(jù)沖突程度。系數(shù)1／（1-k)稱為歸一化因子，其作用是為了避免在合成時將非零的概率賦給空集。

式（1)中，若k＝1，分母為零，此時改式無法進(jìn)行合成；當(dāng)k→1時，即證據(jù)高度沖突時，公式可能導(dǎo)致違反常理的結(jié)果。

為克服這一不足，Yager對D-S合成公式做了改進(jìn)，公式如下（以2個證據(jù)源為例)［9］。

2 實例分析

國家煤礦防治水規(guī)定提出了“預(yù)測預(yù)報、有疑必探、先探后掘、先治后采”的原則。在采礦工作前采用鉆探、物探等方法，結(jié)合在掘進(jìn)面多點(diǎn)傳感器得到的數(shù)據(jù)，再結(jié)合BP神經(jīng)網(wǎng)絡(luò)和D-S證據(jù)理論進(jìn)行預(yù)測，從而提出預(yù)防措施。

根據(jù)水文地質(zhì)環(huán)境分析，山西某礦3＃煤層突水因素判定和指標(biāo)分級標(biāo)準(zhǔn)，確定各級評價指標(biāo)值。突水的安全等級定位4級，分別為：

S1級：0.1，0.1，0.1，0.9；

S2級：0.1，0.1，0.9，0.1；

S3級：0.1，0.9，0.1，0.1；

S4級：0.9，0.1，0.1，0.1；

從該礦開采、探測資料中提取6組數(shù)據(jù)作為構(gòu)建預(yù)測模型的樣本，樣本如表1所示。

表1 系統(tǒng)訓(xùn)練樣本表

采用MATLAB編程建立BP神經(jīng)網(wǎng)絡(luò)，實驗證明單個BP網(wǎng)絡(luò)輸出結(jié)果不穩(wěn)定，所以采用2個BP網(wǎng)絡(luò)作為證據(jù)輸入，通過D-S證據(jù)理論的組合規(guī)則進(jìn)行計算。各數(shù)據(jù)之間的相關(guān)性很小，可以認(rèn)為輸出結(jié)果之間是相互獨(dú)立的。由于神經(jīng)網(wǎng)絡(luò)和測量數(shù)據(jù)都存在誤差，用神經(jīng)網(wǎng)絡(luò)的平方和誤差函數(shù)（SSE)的值來表示不確定性證據(jù)。BP1訓(xùn)練中應(yīng)用的函數(shù)名稱為trainbp，BP2網(wǎng)絡(luò)訓(xùn)練中應(yīng)用的函數(shù)名稱為trainlm，最大循環(huán)次數(shù)為300（epochs)、訓(xùn)練目標(biāo)誤差為10-5（goal)、顯示訓(xùn)練迭代次數(shù)為10（show)。并經(jīng)sim（)函數(shù)仿真驗證。圖3和圖4分別為BP1和BP2的訓(xùn)練結(jié)果。

最后對BP1、BP2的輸出進(jìn)行歸一化處理，即形成證據(jù)理論的證據(jù)。按照D-S證據(jù)組合的規(guī)則進(jìn)行融合，其輸出結(jié)果如表2所示。

圖3 BP1網(wǎng)絡(luò)訓(xùn)練性能圖

圖4 BP2網(wǎng)絡(luò)訓(xùn)練性能圖

表2 融合結(jié)果

3 結(jié)束語

以樣本1為例，BP1、BP2的不確定度分別位0.0205、0.0176，融合后的不確定度降低為 0.0004。實驗結(jié)果表明，通過D-S證據(jù)理論融合后對比單一使用BP神經(jīng)網(wǎng)絡(luò)的輸出數(shù)據(jù)極大地降低了不確定度；顯著提高了井下安全狀態(tài)判斷的準(zhǔn)確度，且與實際數(shù)據(jù)相對較符合。

本文所涉及的煤礦工作面突水評價模型的完成是建立在諸多影響因素基礎(chǔ)之上的，所以用戶提供資料的準(zhǔn)確程度及全面程度直接影響著對系統(tǒng)得出的結(jié)論的可信程度。礦區(qū)提供資料越豐富全面，越能得到較為精確的結(jié)論。即本突水評價模型適合有一定突水資料的煤礦。而對于那些暫時沒有足夠多資料的煤礦，則需根據(jù)實際情況先安置傳感器，然后對所測數(shù)據(jù)進(jìn)行日、月、年積累，建立完善的數(shù)據(jù)樣本及數(shù)據(jù)庫。才能利用本評價模型，進(jìn)行突水狀態(tài)的評價預(yù)測。

［1］張英梅，湯紅穎.基于LabView的煤礦井下排水監(jiān)控系統(tǒng)的設(shè)計［J］.煤炭科學(xué)技術(shù)，2008，36（1)：68-71.

［2］毛紅川，劉志，鄧春濤，王檔良.煤礦突水評價方法探討［J］.河北工業(yè)科技2008，25（4)：197-199.

［3］Dubois H.representation and combination of uncertainty with belief functions and possibility measures［J］.Computational Intelligence，1998（1)：88-101.

［4］姜成志，張紹兵.建立在神經(jīng)網(wǎng)絡(luò)基礎(chǔ)上的煤礦突水評價模型［J］.黑龍江科技學(xué)院學(xué)報，2006，116（1)：9-11.

［5］張永懷，劉君華.采用神經(jīng)網(wǎng)絡(luò)及其改進(jìn)算法改善傳感器特性［J］.傳感器技術(shù)學(xué)報，2002，15（3)：186-188.

［6］Zhang Ying-Mei，Cheng Zhen-Zhen.The Coal Mine Flood Prediction Research Based On Neural Network and D-S Theory of Evidence［J］.《IEEE IITA 2008》，2008，103（12)：311-315.

［7］吳景龍，王巨，王春山.人工神經(jīng)網(wǎng)絡(luò)系統(tǒng)在預(yù)防煤礦突水中應(yīng)用［J］.煤礦安全.2008，39（4)：55-57.

［8］董書寧，靳德武，馮安.煤礦防治水實用技術(shù)及裝備［J］.煤炭科學(xué)技術(shù)，2008，36（3)：8-11.

［9］Yager R R.On the Dempster-Shafer Framework and New Combination Rules［J］.Information Science，1998，41（2)：93-137.

［10］張英梅，程珍珍.D-S證據(jù)理論在煤礦水害預(yù)測中的應(yīng)用［J］.太原理工大學(xué)學(xué)報2008，39（6)：589-591.

［11］董長虹.matlab神經(jīng)網(wǎng)絡(luò)與應(yīng)用［M］.北京：國防工業(yè)出版社，2005：65-85.