趙莊礦瓦斯涌出趨勢(shì)分形預(yù)測(cè)研究

郭雪勇

（山西晉城無(wú)煙煤礦業(yè)集團(tuán)有限責(zé)任公司 趙莊煤業(yè)有限責(zé)任公司，山西 晉城 046605）

近年來(lái)，對(duì)于瓦斯涌出研究的方法關(guān)注神經(jīng)網(wǎng)絡(luò)、灰色系統(tǒng)等理論，但其主觀性過(guò)強(qiáng)，無(wú)法完整表達(dá)瓦斯涌出的客觀性能，實(shí)際效果不佳。由于分形理論和方法可以彌補(bǔ)這些不足，本文以分形理論中的R/S分析法為基礎(chǔ)分析工具，針對(duì)趙莊礦瓦斯涌出量的預(yù)測(cè)進(jìn)行研究。

1 分形預(yù)測(cè)法的應(yīng)用概述

R/S分析法是屬于分形理論中的時(shí)間序列分析法，本文利用已知的時(shí)間序列，計(jì)算趙莊礦瓦斯涌出的分形維數(shù)，并考察該系統(tǒng)的持久性和反持久性，從而對(duì)該時(shí)間序列所表達(dá)的內(nèi)容做出預(yù)測(cè)。礦井瓦斯涌出量的大小受多方面因素影響、錯(cuò)綜復(fù)雜，主要表現(xiàn)為礦井瓦斯涌出量既非確定的、也非完全隨機(jī)的，而是一種模糊系統(tǒng)狀態(tài)。由于瓦斯涌出量的測(cè)定是等時(shí)間間隔的，這正與分形理論中的時(shí)間序列交相呼應(yīng)。所以根據(jù)赫斯特指數(shù)的大小，對(duì)礦井瓦斯涌出的增量趨勢(shì)做出基本預(yù)測(cè)是可行的。

2 R/S分析法對(duì)礦井瓦斯的分形預(yù)測(cè)

首先建立一個(gè)統(tǒng)計(jì)量：赫斯特指數(shù)（H）。瓦斯的涌出是圍繞其平均水平漲落的。漲落的極差是變化的，它依賴(lài)于用于度量的時(shí)間長(zhǎng)度。趨勢(shì)的強(qiáng)度和噪聲的水平可按重標(biāo)極差隨時(shí)間的變化情況度量，即看H比0.5大多少。采用回歸的方法求其斜率，赫斯特指數(shù)、瓦斯涌出量時(shí)間序列的分形維數(shù)以及平均循環(huán)長(zhǎng)度，具體算法如下：

3 實(shí)例分析

3.1 礦井概述

礦井通風(fēng)系統(tǒng)基本情況。趙莊礦井采用中央分區(qū)式通風(fēng)方式，通風(fēng)方法為機(jī)械抽出式。全礦6個(gè)井筒，其中4個(gè)進(jìn)風(fēng)井，2個(gè)回風(fēng)井。張店1號(hào)回風(fēng)井主扇45°運(yùn)行，總回風(fēng)量18 957 m3/min，主扇負(fù)壓2 396 Pa。張店2號(hào)回風(fēng)井主扇現(xiàn)-2.5°運(yùn)行，總進(jìn)風(fēng)量27 161 m3/min，主扇負(fù)壓2 677 Pa。這4臺(tái)主要通風(fēng)機(jī)均按規(guī)定進(jìn)行了性能測(cè)定。目前有3個(gè)大采高綜采工作面，2個(gè)在回采、1個(gè)準(zhǔn)備，回采工作面采用三進(jìn)兩回和一進(jìn)兩回的通風(fēng)方式，掘進(jìn)工作面使用局部通風(fēng)機(jī)進(jìn)行通風(fēng)。為了研究瓦斯的涌出規(guī)律，在趙莊礦進(jìn)行了實(shí)地的調(diào)研和勘查工作，如表1所示。瓦斯涌出的趨勢(shì)預(yù)測(cè)，以趙莊礦瓦斯涌出量為研究對(duì)象，采用R/S分析法預(yù)測(cè)瓦斯涌出的未來(lái)趨勢(shì)。時(shí)間序列是從2012年4月至2012年6月止，根據(jù)表1資料，求出該時(shí)間序列的分形維數(shù)。

3.2 瓦斯涌出量的赫斯特指數(shù)

由表2和表3，可得到圖1。根據(jù)式1～式6計(jì)算

表1 趙莊礦1號(hào)總回風(fēng)風(fēng)流瓦斯涌出情況統(tǒng)計(jì)表 （m3/min）

如下：

由上結(jié)果可知，礦井瓦斯涌出量的分形維數(shù)D為1.8，并且其平均循環(huán)長(zhǎng)度為2d。

圖1 lnR/S隨時(shí)間變化的趨勢(shì)圖

3.3 瓦斯涌出量的H討論

由以上計(jì)算結(jié)果可推得，趙莊礦瓦斯涌出量的H=0.2（0≤H≤0.5），更接近于 0，說(shuō)明趙莊礦的系統(tǒng)為反持久性或遍歷性的時(shí)間序列，在首個(gè)平均生產(chǎn)周期內(nèi)序列不斷增強(qiáng)，下一個(gè)平均周期序列將減弱，就趙莊礦實(shí)際而言，整個(gè)計(jì)算結(jié)果表現(xiàn)為瓦斯涌出量的突變性或易變性較強(qiáng)，相對(duì)比較容易發(fā)生逆轉(zhuǎn)。

3.4 相關(guān)系數(shù)r的討論

根據(jù)趙莊礦基礎(chǔ)瓦斯數(shù)據(jù)處理結(jié)果可知，相關(guān)系數(shù)r=0.73，說(shuō)明數(shù)據(jù)的ln-ln值的線性關(guān)系不是很好，即趙莊礦瓦斯涌出的時(shí)間序列無(wú)長(zhǎng)期記憶性，該瓦斯涌出量表現(xiàn)為異常復(fù)雜的非線性，這種特征存在的原因是多因素融合。

3.5 瓦斯涌出量分形行為的探討

趙莊礦瓦斯涌出的分形維數(shù)D=1.8，其特征表現(xiàn)出一定的不確定性，但分形規(guī)律性也存在，見(jiàn)圖1。lnR/S的值始終在0.2～1.2之間，但平均循環(huán)周期很短。這種狀態(tài)說(shuō)明：趙莊礦瓦斯涌出的影響因素，包括自身因素和外力因素。自身因素有煤本身的結(jié)構(gòu)、賦存等；外力因素有開(kāi)拓開(kāi)采及是否存在抽放等。這些因素相互聯(lián)系、相互制約、錯(cuò)綜復(fù)雜。趙莊礦瓦斯涌出的“諾亞效應(yīng)”特征，使其異常狀態(tài)的研究有了一個(gè)定量化的可能。

表2 趙莊礦1號(hào)總回風(fēng)風(fēng)流瓦斯涌出量的R/S值

表3 趙莊礦1號(hào)總回風(fēng)風(fēng)流瓦斯涌出量的(R/S)ln-ln值

4 結(jié)束語(yǔ)

趙莊礦瓦斯涌出的分形維數(shù)D=1.8，其時(shí)間序列突變性相對(duì)較強(qiáng)。趙莊礦瓦斯涌出預(yù)測(cè)期望值較低，前一個(gè)生產(chǎn)周期的數(shù)據(jù)無(wú)法完整預(yù)測(cè)下一個(gè)生產(chǎn)周期的瓦斯實(shí)際狀態(tài)。趙莊礦瓦斯涌出具有的“諾亞效應(yīng)”特征，為深入研究定量預(yù)測(cè)其瓦斯涌出狀態(tài)提供一個(gè)依據(jù)。

[1]呂伏，梁冰，孫維吉，等.基于主成分回歸分析法的回采工作面瓦斯涌出量預(yù)測(cè)[J].煤炭學(xué)報(bào)，2012，37(1)：113－116．

[2]付華，姜偉，單欣欣.基于耦合算法的煤礦瓦斯涌出量預(yù)測(cè)模型研究[J].煤炭學(xué)報(bào)，2012，37(4)：654－658．

[3]黃文標(biāo)，施式亮.基于改進(jìn)Lyapunov指數(shù)的瓦斯涌出時(shí)間序列預(yù)測(cè)[J].煤炭學(xué)報(bào)，2009，34(12)：1665－1668．

[4]徐精彩，趙慶賢，鄧軍，等.礦井瓦斯涌出量時(shí)間序列的分形特性分析[J].遼寧工程技術(shù)大學(xué)學(xué)報(bào)，2004（1）：1-4．

[5]劉新喜，王勇，趙云勝.回采工作面瓦斯涌出特征及灰色預(yù)測(cè)模型[J].中國(guó)安全科學(xué)學(xué)報(bào)，2001(2)：11-16．

[6]金芳勇.基于分形—混沌理論的煤與瓦斯突出預(yù)測(cè)研究(碩士論文)[D].淮南：安徽理工大學(xué)，2006．

[7]謝生榮.綜采工作面的瓦斯涌出規(guī)律及瓦斯涌出量的預(yù)測(cè)(碩士論文)[D].太原：太原理工大學(xué)，2005．

[8]聶百勝,等.煤與瓦斯突出預(yù)測(cè)技術(shù)研究現(xiàn)狀及發(fā)展趨勢(shì)[J].中國(guó)安全學(xué)報(bào)，2003(6)：40～43．

[9]辛厚文.分形理論及其應(yīng)用[M].合肥:中國(guó)科學(xué)技術(shù)大學(xué)出版社，1993．