李世民
(山西省臨汾市蒲縣黑龍關(guān)鎮(zhèn)豹子溝煤業(yè), 山西 臨汾 041000)
煤與瓦斯突出是煤礦的主要災害,煤與瓦斯突出的預測技術(shù)對保證煤礦安全生產(chǎn)具有重要的意義。多位學者針對煤與瓦斯突出的預測方法進行了研究,劉曉光等[1]提出運用PCA-聚類分析預測煤與瓦斯突出;龔星宇等[2]開展了獨立成分分析在瓦斯?jié)舛阮A測中的應用研究;馬晟翔等[3]將因子分析與BP神經(jīng)網(wǎng)絡結(jié)合,提出一種改進BP神經(jīng)網(wǎng)絡預測的方法作為煤與瓦斯突出結(jié)果的預測;溫廷新等[4]開展了基于PSOBP-AdaBoost模型的瓦斯涌出量分源預測研究;董曉雷等[5]基于SVM耦合遺傳算法,作為回采工作面瓦斯涌出量結(jié)果的預測;金洪偉等[6]運用主成分分析法對煤礦瓦斯涌出量進行了預測。以上方法在預測煤與瓦斯突出的結(jié)果方面誤差偏大,主要是因為影響瓦斯突出的因素復雜多樣,各因素之間存在非線性關(guān)系。且影響因素有主要因素和次要因素,需要運用數(shù)學方法處理這些影響結(jié)果權(quán)重的問題。熵權(quán)法和灰色關(guān)聯(lián)度數(shù)學方法的結(jié)合能夠簡化預測指標間的相互聯(lián)系,獲得主要影響因素與突出結(jié)果的關(guān)聯(lián)度。為此,提出將灰色熵權(quán)法和BP神經(jīng)網(wǎng)絡法相結(jié)合,通過對預測指標加權(quán)法處理得到影響因素與結(jié)果的關(guān)聯(lián)度,將主要影響因素作為BP神經(jīng)網(wǎng)絡輸入層參數(shù)的個數(shù),簡化BP神經(jīng)網(wǎng)絡結(jié)構(gòu),提高迭代與運算效率。以某礦典型的煤與瓦斯突出為數(shù)據(jù)樣本,運用灰色熵權(quán)關(guān)聯(lián)度與BP神經(jīng)網(wǎng)絡法的煤與瓦斯突出預測,驗證所提預測方法的可行性和準確性。
灰色關(guān)聯(lián)分析法是對于一個目標系統(tǒng)進行定量的直觀描述和比較各因素之間的關(guān)聯(lián)性,基本理論方法為選出目標序列作為母序列,各影響目標因素作為子序列,運用數(shù)學的方法來定量和比較各影響因素對于結(jié)果的緊密程度。進行灰色關(guān)聯(lián)度分析時,根據(jù)以往的數(shù)據(jù)處理方法,要進行數(shù)據(jù)的無量綱化處理,便于得到合理的、科學的結(jié)論[7].
考慮整個體系中各指標所占權(quán)重并非完全相等的特點,在使用客觀的事實數(shù)據(jù)中,依靠數(shù)學理論分析的基礎,在灰色關(guān)聯(lián)度數(shù)學方法的基礎上,再進行熵權(quán)計算,得到的權(quán)重關(guān)聯(lián)度客觀地分析了各屬性對于結(jié)果的影響,有利于提高分析的計算精度。
(1)
式中:γi為關(guān)聯(lián)度;n為元素的個數(shù);k為目數(shù)據(jù)序列,取1,2,3…n;ωj為關(guān)聯(lián)系數(shù);ξi為關(guān)聯(lián)函數(shù)。
BP具有包含多個隱含層的網(wǎng)絡特點,具備處理線性不可分問題的能力。應用最廣的是多層前饋網(wǎng)絡及誤差反向傳播學習算法,網(wǎng)絡采用誤差反向傳播算法進行學習。BP神經(jīng)網(wǎng)絡中的內(nèi)部計算,數(shù)據(jù)由輸入層經(jīng)過隱含層逐層逐步向后傳播。網(wǎng)絡權(quán)值由輸出層逐漸向前修正網(wǎng)絡的權(quán)值大小,網(wǎng)絡層層相連,實現(xiàn)彼此的映射關(guān)系。BP網(wǎng)絡是前向神經(jīng)網(wǎng)絡的核心,在人工智能等領域廣泛解決分類識別的問題。
選用某礦現(xiàn)場實際煤與瓦斯突出的具體數(shù)據(jù)作為研究煤與瓦斯突出的影響因素數(shù)據(jù),見表1.
表1 影響因素數(shù)據(jù)表
2.2.1 熵權(quán)法
應用Matlab2018a數(shù)學計算軟件,重點分析數(shù)據(jù)主要影響因素的權(quán)重大小。由于煤與瓦斯突出各影響因素之間的數(shù)量級和綱量不同,數(shù)據(jù)之間具有比較大的差異性,建立了19×8原始數(shù)據(jù),運用灰色熵權(quán)法分析之前,需對原始數(shù)據(jù)進行歸一化處理,見表2.
對歸一化后的數(shù)據(jù)運用公式得到各指標的權(quán)重,其中煤層最大地應力0.38 N、瓦斯壓力0.18 MPa、瓦斯含量0.06 m3/t、頂?shù)装鍘r性0.41%、與斷裂距離0.11 m、煤層厚度0.12 m、開采深度0.04 m、絕對瓦斯涌出量0.05 m3/d. 在影響煤與瓦斯突出因素中,頂?shù)装鍘r性0.41%在整個體系中所占權(quán)重值最大,這也說明它對于煤與瓦斯突出危險性預測起到的作用最大。
2.2.2 灰色關(guān)聯(lián)熵法
應用Matlab2018a數(shù)學計算軟件,對影響因素的數(shù)據(jù)進行預處理,重點分析各影響因素對結(jié)果的影響程度。
1) 原始數(shù)據(jù)均值化。
由于煤與瓦斯突出各影響因素之間的數(shù)量級和綱量不同,數(shù)據(jù)之間具有比較大的差異性。建立了19×9原始數(shù)據(jù),運用灰色熵權(quán)法分析之前,有必要進行原始數(shù)據(jù)均值化。
2) 灰色關(guān)聯(lián)熵指標系數(shù)。
應用均值化后的數(shù)據(jù)進行灰色關(guān)聯(lián)熵分析,確定反應瓦斯突出特征的參考數(shù)列與影響結(jié)果的比較數(shù)列,計算灰色關(guān)聯(lián)熵指標系數(shù),如表3.
表3 指標系數(shù)表
由表3可知,各影響因素對煤與瓦斯突出危險性影響的關(guān)聯(lián)度,其中瓦斯壓力和頂?shù)装鍘r性的加權(quán)灰色關(guān)聯(lián)熵值分別為0.126 7和0.240 6,說明兩者對煤與瓦斯突出的程度影響最大。而灰色關(guān)聯(lián)度值中,瓦斯含量為0.803 7 m3/t,表明它對煤與瓦斯突出的程度影響大。因此,在進行煤與瓦斯突出的預測時,可重點對頂?shù)装鍘r性、瓦斯壓力、瓦斯含量作為BP神經(jīng)網(wǎng)絡的輸入層,這樣能夠使預測效果更加準確。
選取煤與瓦斯突出的3個主要影響因素(瓦斯壓力、瓦斯含量、頂?shù)装鍘r性)作為BP神經(jīng)網(wǎng)絡輸入層參數(shù)。將發(fā)生煤與瓦斯突出設置成1,將不發(fā)生煤與瓦斯突出設置成0,作為BP神經(jīng)網(wǎng)絡輸出層。在BP網(wǎng)絡中,如果預測值≥0.5,輸出值為1,預測值<0.5時,輸出結(jié)果為0. 使用feedforward網(wǎng)絡作為輸入和輸出的映射函數(shù),其中trainafcn屬于訓練網(wǎng)絡性能的函數(shù),運用訓練函數(shù)trainbfg-BFGS算法(擬牛頓反向傳播算法)作為訓練函數(shù)。編輯劃分數(shù)據(jù)程序設計,分別任意選取突出10個與非突出5個作為訓練樣本,建立BP神經(jīng)網(wǎng)絡,最后用檢驗樣本檢驗訓練好的網(wǎng)絡作為預測。經(jīng)過多次訓練,隱含層包括10個隱含層效果最好。即模型結(jié)構(gòu)為3-10-1. 用最后4個樣本檢驗改進后的BP神經(jīng)網(wǎng)絡,結(jié)果見表4.
表4 煤與瓦斯突出預測表
由表4可見,采用改進后的BP神經(jīng)網(wǎng)絡預測瓦斯發(fā)生突出的情況下,預測相對誤差最大為5.11%,最小為0.36%.預測瓦斯不突出相對誤差為2.33%. 預測結(jié)果的相對誤差均小于6%,表明該設計的灰色熵權(quán)BP網(wǎng)絡模型預測煤與瓦斯突出是可行的,并且結(jié)果符合現(xiàn)場實際情況。
1) 設計瓦斯突出預測的新模型,引入了灰色熵權(quán)法來分析各影響因素的大小,并對各影響因素客觀賦予權(quán)重進行了排序,避免了主觀性判斷的影響。建立的灰色關(guān)聯(lián)熵分析模型能很好對煤與瓦斯突出起到降維的作用,且具有較高的分辨率,并提高了BP神經(jīng)網(wǎng)絡的運算速度和準確率。
2) 在得到影響因素與突出的關(guān)聯(lián)度大小的基礎上,選取主要的影響因素數(shù)據(jù)作為BP模型的輸入層。由訓練及測試結(jié)果可知,基于灰色關(guān)聯(lián)熵的BP神經(jīng)網(wǎng)絡模型適用于煤與瓦斯突出預測,且有強的預測能力及更高的運算效率。
3) 采用改進后的BP神經(jīng)網(wǎng)絡預測方法對某礦進行煤與瓦斯突出預測,4個預測樣本的相對誤差分別為5.11%,0.36%,4.64%,2.33%,預測結(jié)果的相對誤差均小于6%,表明該設計的灰色熵權(quán)BP網(wǎng)絡模型預測煤與瓦斯突出是可行的,并且結(jié)果符合現(xiàn)場實際情況。