王福生,張志明,董憲偉
(1.華北理工大學(xué) 礦業(yè)工程學(xué)院,河北 唐山 063210;2.河北省礦業(yè)開發(fā)與安全技術(shù)實驗室,河北 唐山 063210)
礦井災(zāi)害事故對人民的生命財產(chǎn)安全造成了巨大威脅,事故起源以煤的自燃發(fā)火居多。煤的自燃傾向性是煤的固有屬性。煤本身的復(fù)雜性致使科研難度較大、未知領(lǐng)域較多[1],科研人員對此研究不斷深入:鄧軍等[2]建立了礦井溫度預(yù)測的PSO-SVR模型用于判別煤自燃溫度;王福生等[3]通過進行復(fù)采與原生煤樣的實驗發(fā)現(xiàn)官能團差異導(dǎo)致的煤自燃能力的不同;王海燕[4]等發(fā)現(xiàn)煤在氧化過程存在“自限制”特征;唐一博[5]從活性基團角度出發(fā)模擬煤自燃過程,發(fā)現(xiàn)植酸等抗氧化劑對煤的低溫氧化有良好的抑制效果。當(dāng)前對煤自燃的研究多基于單一因素,而煤自燃影響因素眾多,因此,考慮多因素影響煤自燃傾向的研究顯得非常重要。本文基于煤的多項基礎(chǔ)實驗數(shù)據(jù)建立BP神經(jīng)網(wǎng)絡(luò)模型,將訓(xùn)練成功的網(wǎng)絡(luò)模型用于預(yù)測煤的自燃傾向性,以便豐富礦井煤自燃預(yù)測的技術(shù)手段。
影響煤自燃的原因可分為外部原因和內(nèi)在原因。外部原因包含礦井壓力、水文地質(zhì)等不可控因素,在此不做贅述。內(nèi)在原因可細分為煤的組成和結(jié)構(gòu)兩方面:煤的組成包括元素成分、煤巖組成成分等;結(jié)構(gòu)包括煤的孔隙結(jié)構(gòu)、分子結(jié)構(gòu)等。從煤的組成與結(jié)構(gòu)方面入手能準確、全面地探究煤自燃傾向性的差異性原因[6]。
從內(nèi)在原因說,影響煤的自燃傾向性的主要因素有:①碳含量。碳為煤的骨架元素,也是煤燃燒放熱中產(chǎn)生熱能的主導(dǎo)元素,化學(xué)性質(zhì)穩(wěn)定的芳香結(jié)構(gòu)中也是以碳元素為主。煤中碳含量增加,煤的自燃傾向性減小。②鏡質(zhì)組含量。鏡質(zhì)組是指植物殘體轉(zhuǎn)變的有機成分。文獻表明[6-7],對于中低階煤,鏡質(zhì)組含量越高,自燃能力越大;對于中高階煤,鏡質(zhì)組含量越高,自燃能力越低。③固定碳含量。煤經(jīng)過熱解過程殘留下來的有機質(zhì)為固定碳,是高分子化合物的混合物,其含量越高,煤階越高,煤自燃傾向性越低。④比表面積。比表面積值越大,煤與氧氣接觸的面積越大,越易發(fā)生化學(xué)反應(yīng)致使煤熱量集聚引發(fā)自燃。⑤微孔占比。微孔相比較于中大孔,內(nèi)徑小,空氣阻力大,分子結(jié)構(gòu)緊密,吸附能力強,氣體流通性低。微孔占比大,在煤自燃初期會造成反應(yīng)遲緩,煤自燃能力下降[8]。⑥羥基含量。通過傅里葉紅外光譜實驗發(fā)現(xiàn)煤中存在大量羥基結(jié)構(gòu),其結(jié)構(gòu)穩(wěn)定性差,氧化性能活躍,通過參與各項基元反應(yīng)影響煤自燃傾向性,即羥基含量越高,煤自燃傾向性越大。
BP神經(jīng)網(wǎng)絡(luò)原理表明,樣本選取數(shù)量盡可能大,才能做量化處理[9]。結(jié)合實際情況與數(shù)據(jù)完整性,本文選取內(nèi)蒙古、新疆、河北、陜西等地的12個礦區(qū)新鮮煤樣,遵照《煤層煤樣采取方法》(GB/T482-2008)分別進行煤的氧化動力學(xué)實驗、工業(yè)分析實驗、元素分析實驗、煤巖組成實驗、傅里葉紅外光譜實驗及低溫氮氣吸附實驗,將取得的煤樣基礎(chǔ)實驗數(shù)據(jù)作為測試樣本(見表1)。特別說明的是:羥基數(shù)據(jù)為傅里葉紅外光譜實驗結(jié)果中的晰光度數(shù)據(jù)經(jīng)ONMIC軟件的“自動基線校正”后運用Peakfit軟件進行分峰處理的3 685~3 610 cm-1處的含量;根據(jù)《煤自燃傾向性的氧化動力學(xué)測定方法》(AQ/T1068-2008)進行氧化動力學(xué)實驗,測定煤的自燃傾向性綜合判定指數(shù)可信度高,多為業(yè)界采用。構(gòu)建BP神經(jīng)網(wǎng)絡(luò),利用Matlab R2014b編寫代碼進行網(wǎng)絡(luò)訓(xùn)練,表1為模型的訓(xùn)練樣本數(shù)據(jù),其中1-11號煤樣為學(xué)習(xí)樣本,12號煤樣用于檢驗網(wǎng)絡(luò)可靠性。BP神經(jīng)網(wǎng)絡(luò)中,P為輸入矩陣,T為輸出矩陣,P(6×11)為煤自燃傾向性影響因素,T(1×11)為煤自燃傾向性綜合判定指數(shù)。
表1 模型的訓(xùn)練樣本數(shù)據(jù)
樣本輸入及輸出數(shù)據(jù)利用歸一化公式(式1)與反歸一化公式(式2)進行處理,使網(wǎng)絡(luò)訓(xùn)練樣本數(shù)據(jù)處于區(qū)間[-1,1]中[10]。
(1)
(2)
式中,Xn為各參數(shù)第n個歸一化處理數(shù)據(jù);X1為各參數(shù)第1個數(shù)據(jù);X2為反歸一化輸出數(shù)據(jù);Xmax,Xmin為各類參數(shù)中最大值與最小值。
選用3層改進的BP神經(jīng)網(wǎng)絡(luò),基本結(jié)構(gòu)為:輸入層6個節(jié)點,分別代表6個影響因素;輸出層1個節(jié)點,代表煤的自燃傾向性綜合判定指數(shù);依據(jù)經(jīng)驗公式及模擬訓(xùn)練效果將隱含層節(jié)點數(shù)選為5個。網(wǎng)絡(luò)框架如圖1所示。
圖1 煤自燃傾向性預(yù)測的神經(jīng)網(wǎng)絡(luò)框架
利用Matlab R2014b編寫改進的BP神經(jīng)網(wǎng)絡(luò)代碼,激活軟件自帶的工具箱(Neural Network Training),采用L-M優(yōu)化算法,訓(xùn)練函數(shù)選用trainlm,傳遞函數(shù)選用tansig,目標誤差設(shè)為0.001,詳細代碼見圖2。
圖2 BP神經(jīng)網(wǎng)絡(luò)的Matlab代碼
運行編寫后的代碼,結(jié)果如表2,圖3,圖4所示。根據(jù)運行結(jié)果可知,模型經(jīng)75次迭代后達到最佳效果,相關(guān)系數(shù)R為0.999 94,呈高度相關(guān)。經(jīng)計算分析發(fā)現(xiàn):12號煤樣預(yù)測值與計算值相差極小,相對誤差為1.01%,遠低于10%,滿足工程預(yù)測精度要求。由此可見,該網(wǎng)絡(luò)預(yù)測準確度高,預(yù)測效果顯著,模型通過了檢驗。
表2 實際值與預(yù)測值對比
圖3 迭代次數(shù)效果圖
圖4 回過效果圖
為說明基于BP神經(jīng)網(wǎng)絡(luò)的煤自燃傾向性預(yù)測系統(tǒng)的實用性,選擇唐山礦9煤層、11煤層和荊各莊礦的新鮮煤樣進行各項實驗,實驗結(jié)果如表3所示。
表3 唐山礦9煤層、11煤層和荊各莊礦煤樣實驗數(shù)據(jù)
分別將各參數(shù)值按上述次序輸入神經(jīng)網(wǎng)絡(luò)代碼中進行網(wǎng)絡(luò)仿真預(yù)測,并利用相對誤差與絕對誤差法對氧化動力學(xué)煤自燃傾向性計算值與預(yù)測值進行對比分析,結(jié)果見表4。
表4 唐山礦9煤層、11煤層和荊各莊礦煤自燃傾向性綜合判定指數(shù)計算值與預(yù)測值對比
經(jīng)BP神經(jīng)網(wǎng)絡(luò)仿真預(yù)測,得到唐山礦9煤層、11煤層與荊各莊礦的煤自燃傾向性綜合判定指數(shù)的預(yù)測值,與理論計算值的相對誤差分別為4.08%,3.34%,3.19%,誤差值較小,說明了該網(wǎng)絡(luò)模型的可行性與準確性。
(1)影響煤自燃傾向性的因素有多種,量綱不一,且部分因素?zé)o法量化,據(jù)此本文選取碳含量、鏡質(zhì)組含量、固定碳含量、比表面積、微孔占比與羥基含量6項煤的組成與結(jié)構(gòu)方面的基礎(chǔ)數(shù)據(jù)進行建模,因此,此模型具有全面性、適應(yīng)性特點。
(2)基于BP神經(jīng)網(wǎng)絡(luò)的煤自燃傾向性預(yù)測模型的實際運行結(jié)果與行業(yè)標準誤差均小于5%,說明該模型可信度高。
(3)神經(jīng)網(wǎng)絡(luò)模型需要大量樣本數(shù)據(jù)才能達到高度準確,但由于目前受實際煤樣采集及實驗難度所限,完整數(shù)據(jù)僅限于此,相對誤差較大。筆者在后續(xù)研究中會盡可能多地收集數(shù)據(jù)來完善神經(jīng)網(wǎng)絡(luò)模型,提高計算精度,而且在大樣本條件下亦考慮增添煤自燃傾向性的其他影響因素。
(4)基于BP神經(jīng)網(wǎng)絡(luò)的煤自燃傾向性預(yù)測模型豐富了礦井煤自燃預(yù)測的技術(shù)手段,且模型可依據(jù)所選的關(guān)鍵指標進行更改,操作便捷,有較大的實用價值。