艾 莉 程加堂
(紅河學(xué)院工學(xué)院,云南 蒙自 661199)
如何可靠、準(zhǔn)確地預(yù)測(cè)煤與瓦斯突出強(qiáng)度,是煤礦安全技術(shù)研究中的重要課題。國(guó)內(nèi)開(kāi)采實(shí)踐表明,煤與瓦斯突出的發(fā)生呈現(xiàn)區(qū)域性分布,而災(zāi)害發(fā)生區(qū)域只占開(kāi)采區(qū)的8%~20%。因此,首先要在宏觀上研究煤與瓦斯突出的影響因素,對(duì)突出危險(xiǎn)進(jìn)行區(qū)域預(yù)測(cè)。由于引起煤與瓦斯突出的因素很多,而且影響突出的因素多數(shù)是定性的或是模糊相似的,各影響因素間又呈現(xiàn)出復(fù)雜的非線性關(guān)系,從而難以用經(jīng)典的理論建立精確的數(shù)學(xué)預(yù)測(cè)模型[1]。
人工神經(jīng)網(wǎng)絡(luò)是能表現(xiàn)極為復(fù)雜動(dòng)態(tài)行為的網(wǎng)絡(luò),具有高度的非線性和容錯(cuò)性,在現(xiàn)代高科技發(fā)展過(guò)程中,應(yīng)用越來(lái)越廣泛,它已是控制工程、信號(hào)處理等領(lǐng)域不可缺少的工具。BP神經(jīng)網(wǎng)絡(luò)是目前應(yīng)用比較普遍的一種神經(jīng)網(wǎng)絡(luò),但存在收斂速度慢和易陷入局部極值點(diǎn)等缺陷[2]。PSO利用種群中個(gè)體對(duì)信息的共享,自適應(yīng)的調(diào)整搜索方向,具有內(nèi)在的隱并行性和很好的全局尋優(yōu)能力,并且流程簡(jiǎn)單易于實(shí)現(xiàn),算法參數(shù)簡(jiǎn)潔,無(wú)需復(fù)雜的調(diào)整。本文將此法應(yīng)用于煤與瓦斯突出的預(yù)測(cè),即:利用PSO訓(xùn)練BP網(wǎng)絡(luò)的權(quán)值和閾值,綜合神經(jīng)網(wǎng)絡(luò)能自動(dòng)在海量數(shù)據(jù)中提取隱藏規(guī)律和PSO全局搜索、收斂速度快的PSO—BP模型,本文收集了國(guó)內(nèi)具有突出代表性的8個(gè)礦井作為學(xué)習(xí)訓(xùn)練樣本,對(duì)煤與瓦斯突出強(qiáng)度進(jìn)行了預(yù)測(cè),通過(guò)對(duì)預(yù)測(cè)精度和誤差方面的分析,驗(yàn)證了該方法的準(zhǔn)確性和適用性。
在煤與瓦斯突出預(yù)測(cè)中,煤與瓦斯突出影響因素的合理選取是預(yù)測(cè)的前提,它直接影響著預(yù)測(cè)結(jié)果的準(zhǔn)確性。中華人民共和國(guó)煤炭工業(yè)部頒布的煤與瓦斯突出礦井鑒定規(guī)范中規(guī)定,判定煤層是否具有煤與瓦斯突出危險(xiǎn)的指標(biāo)可用煤的破壞類型、瓦斯放散初速度指標(biāo)、煤的堅(jiān)固性系數(shù)和煤層瓦斯壓力。實(shí)際上煤與瓦斯突出是地應(yīng)力、瓦斯和煤結(jié)構(gòu)的物理學(xué)性質(zhì)等多個(gè)因素綜合作用的結(jié)果。影響煤與瓦斯突出強(qiáng)度的大小與眾多因素相關(guān),結(jié)合我國(guó)煤礦突出的特點(diǎn),本文選取煤層瓦斯含量(m3/t)、瓦斯壓力(MPa)、瓦斯放散初速度Δp (m/s)、煤的堅(jiān)固性系數(shù)f、軟分層煤體厚度(m)、煤體破壞類型6項(xiàng)突出指標(biāo)作為神經(jīng)網(wǎng)絡(luò)的輸入向量。
在d維搜索空間中,有m個(gè)粒子組成一個(gè)群體,對(duì)每個(gè)粒子i(i=1,2,…,m)包含一個(gè)速度向量Vi=(vi1,vi2,…,vid)和位置向量Xi=(xi1,xi2,…,xid)。該粒子經(jīng)歷的最優(yōu)位置記為pi=(pi1,pi2,…,pid),整個(gè)粒子群所經(jīng)歷的最優(yōu)位置記為pg= (pg1,pg2,…,pgd)。在每次迭代中,粒子i通過(guò)跟蹤這兩個(gè)最優(yōu)值,按照下式來(lái)更新速度和位置。
其中,w為慣性權(quán)因子,它使粒子保持運(yùn)動(dòng)慣性;c1、c2為學(xué)習(xí)因子,一般c1=c2;k為迭代次數(shù);r1、r2為兩隨機(jī)數(shù),取值區(qū)間為[0,1]。粒子在調(diào)整自己的位置時(shí),還受到最大速度vmax的限制。如果對(duì)粒子的加速導(dǎo)致其在某維的速度超過(guò)在該維的最大速度,則該維的速度就被限制在最大速度上[3]。
為了平衡粒子群算法的全局和局部搜索能力,慣性因子常采用(3)式進(jìn)行自適應(yīng)調(diào)整。
式中,wmax、wmin分別表示w的最大值和最小值,k表示當(dāng)前迭代步數(shù),kmax表示最大迭代步數(shù),通常取wmax=0.9,wmin=0.4。
盡管標(biāo)準(zhǔn)粒子群算法收斂速度快,具有很強(qiáng)的通用性,但其存在早熟收斂問(wèn)題。借鑒遺傳算法的變異思想,在標(biāo)準(zhǔn)粒子群算法中引入變異操作,即對(duì)某些變量以一定的概率重新初始化。變異操作拓展了在迭代中不斷縮小的種群搜索空間,使粒子能夠跳出先前搜索到的最優(yōu)值位置,在更大的空間中繼續(xù)搜索。這種方法既保持了種群的多樣性,又可提高搜索出最優(yōu)值的可能性[4]。
粒子群神經(jīng)網(wǎng)絡(luò)算法是通過(guò)粒子群各粒子在其搜索空間內(nèi)搜索BP網(wǎng)絡(luò)各層間的權(quán)值及閾值。首先,將BP網(wǎng)絡(luò)的權(quán)、閾值賦予初始值。其次,根據(jù)所確定BP網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu),計(jì)算出網(wǎng)絡(luò)的輸出值。再次,用粒子群算法搜索出最優(yōu)位置,使得適應(yīng)度函數(shù)達(dá)到最小。在這里,對(duì)給定的8組數(shù)據(jù),分別選用標(biāo)準(zhǔn)粒子群算法和帶變異操作的粒子群算法進(jìn)行BP網(wǎng)絡(luò)權(quán)、閾值尋優(yōu)。最后,采用BP網(wǎng)絡(luò)對(duì)所獲得的最優(yōu)值進(jìn)行二次學(xué)習(xí)訓(xùn)練,得到最終的煤與瓦斯突出強(qiáng)度的預(yù)測(cè)模型。
在仿真實(shí)驗(yàn)中,本文所采用的粒子群神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)為6-10-1,訓(xùn)練誤差0.00001。粒子群優(yōu)化算法的參數(shù)設(shè)置為:種群規(guī)模25,最大迭代次數(shù)50,學(xué)習(xí)因子c1=c2=2,慣性權(quán)重w隨迭代次數(shù)由0.9線性地減小為0.3。帶變異操作的粒子群算法的變異概率p=0.2,MATLAB代碼為:
利用表1中煤與瓦斯突出影響因素原始數(shù)據(jù)[5]對(duì)粒子群網(wǎng)絡(luò)進(jìn)行訓(xùn)練,在訓(xùn)練過(guò)程中當(dāng)滿足目標(biāo)精度要求或達(dá)到最大迭代次數(shù)時(shí),自動(dòng)停止訓(xùn)練。粒子群神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型在經(jīng)過(guò)27次迭代后達(dá)到預(yù)定誤差精度0.00001,煤與瓦斯突出強(qiáng)度的預(yù)測(cè)結(jié)果見(jiàn)表2。
表1 煤與瓦斯突出影響因素原始數(shù)據(jù)Table1 The original data of coal and gas outburst influence factors
表2 煤與瓦斯突出強(qiáng)度預(yù)測(cè)結(jié)果表Table2 The forecast results of the strength of coal and gas outburst
由表2可知,本文所建立的粒子群神經(jīng)網(wǎng)絡(luò)煤與瓦斯突出強(qiáng)度的預(yù)測(cè)模型效果較好,相對(duì)誤差的范圍在-0.04%~1.52%之間,說(shuō)明本文所采用的方法輸出數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)的吻合程度較高,能夠滿足實(shí)際應(yīng)用的需要。
BP神經(jīng)網(wǎng)絡(luò)的高度非映射性和自我學(xué)習(xí)能力,很好地解決了各影響因素之間與煤與瓦斯突出強(qiáng)度的非線性關(guān)系,準(zhǔn)確地把握了各影響因素對(duì)煤與瓦斯突出規(guī)律的影響程度。粒子群優(yōu)化算法提高了BP神經(jīng)網(wǎng)絡(luò)的收斂速度,改善其泛化能力。實(shí)例表明,本文提出的預(yù)測(cè)方法是有效的,對(duì)于準(zhǔn)確預(yù)報(bào)煤與瓦斯突出危險(xiǎn)狀況,采取針對(duì)性的防突和安全措施具有較好的指導(dǎo)作用。
[1]張子戌.基于模糊聚類分析和模糊模式識(shí)別的煤與瓦斯突出預(yù)測(cè)[J].煤田地質(zhì)與勘探,2007,35(3):22-24.
[2]楊艷春,趙瑋燁.BP神經(jīng)網(wǎng)絡(luò)在煤與瓦斯突出預(yù)測(cè)的研究[J].蘭州交通大學(xué)學(xué)報(bào),2009,28(6):26-28.
[3]陸寧,周建中,何耀耀.粒子群優(yōu)化的神經(jīng)網(wǎng)絡(luò)模型在短期負(fù)荷預(yù)測(cè)中的應(yīng)用[J].電力系統(tǒng)保護(hù)與控制,2010,38 (12):65-68.
[4]MATLAB中文論壇.MATLAB神經(jīng)網(wǎng)絡(luò)30個(gè)案例分析[M].北京:北京航空航天大學(xué)出版社,2010.
[5]李春輝,陳日輝,蘇恒瑜.BP神經(jīng)網(wǎng)絡(luò)在煤與瓦斯突出預(yù)測(cè)的應(yīng)用[J].礦冶,2010,19(3):21-23.