本發(fā)明涉及煤礦掘進(jìn)工作面瓦斯災(zāi)害的預(yù)警領(lǐng)域，特別涉及一種基于正態(tài)分布檢驗(yàn)辨識瓦斯涌出異常的瓦斯災(zāi)害預(yù)警方法。

背景技術(shù)：

《煤礦安全規(guī)程》規(guī)定，高瓦斯礦井、煤(巖)與瓦斯突出礦井，都應(yīng)裝備煤礦安全生產(chǎn)監(jiān)測監(jiān)控系統(tǒng)。但目前的瓦斯災(zāi)害預(yù)警方法側(cè)重于實(shí)時(shí)的超限判斷和指標(biāo)型預(yù)警，屬于事后評估和預(yù)警，當(dāng)被測地點(diǎn)瓦斯?jié)舛瘸^規(guī)定濃度時(shí)，進(jìn)行斷電、報(bào)警，而缺乏對預(yù)測前后瓦斯?jié)舛葦?shù)據(jù)的有效分析。就已發(fā)生瓦斯災(zāi)害的事故來看，有些事故發(fā)生前，瓦斯?jié)舛纫恢痹谝?guī)定的濃度范圍內(nèi)。因此，研究瓦斯監(jiān)測數(shù)據(jù)本身的特征，辨識瓦斯?jié)舛犬惓｜c(diǎn)(瓦斯?jié)舛冗^高或者過低的點(diǎn))，進(jìn)而確定工作面前方應(yīng)力、瓦斯壓力和強(qiáng)度等一些影響工作面瓦斯?jié)舛鹊囊蛩匕l(fā)生著變化，以此發(fā)出預(yù)警信息，勢必能為礦井瓦斯災(zāi)害的風(fēng)險(xiǎn)預(yù)警提供有效的輔助手段。

技術(shù)實(shí)現(xiàn)要素：

針對現(xiàn)有技術(shù)側(cè)重于實(shí)時(shí)的超限判斷和指標(biāo)型預(yù)警的不足，本發(fā)明提出了基于正態(tài)分布檢驗(yàn)辨識瓦斯涌出異常的瓦斯災(zāi)害預(yù)警方法，可廣泛應(yīng)用于煤礦瓦斯災(zāi)害的預(yù)測預(yù)警中。

為了解決上述技術(shù)問題，本發(fā)明是通過以下技術(shù)方案實(shí)現(xiàn)的：

基于正態(tài)分布檢驗(yàn)辨識瓦斯涌出異常的瓦斯災(zāi)害預(yù)警方法，包括：

步驟1、從煤礦安全生產(chǎn)監(jiān)測監(jiān)控系統(tǒng)中獲得煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列X_(t){x_t，t＝1，2，…，n}，其時(shí)間序列長度為T₀＝n，t表示采樣時(shí)刻，x_t表示t時(shí)刻的瓦斯?jié)舛龋?/p>

步驟2、對煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列X_(t)進(jìn)行Lilliefors正態(tài)分布檢驗(yàn)；

步驟3、對煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列X_(t)進(jìn)行Shapiro-Wilk正態(tài)分布檢驗(yàn)；

步驟4、經(jīng)過Lilliefors正態(tài)分布檢驗(yàn)和Shapiro-Wilk正態(tài)分布檢驗(yàn)，判斷瓦斯?jié)舛鹊臅r(shí)間序列X_(t)是否服從正態(tài)分布：是，則當(dāng)前t＝s時(shí)刻煤礦掘進(jìn)工作面瓦斯涌出處于正常階段，取下一時(shí)刻更新后瓦斯?jié)舛鹊臅r(shí)間序列X_(s+1)，執(zhí)行步驟2至步驟4；否，則當(dāng)前t＝s時(shí)刻煤礦掘進(jìn)工作面瓦斯涌出處于異常階段，執(zhí)行步驟5；

步驟5、取下一時(shí)刻的瓦斯?jié)舛鹊臅r(shí)間序列X_(s+1)，執(zhí)行步驟2至步驟4，若該瓦斯?jié)舛鹊臅r(shí)間序列X_(s+1)服從正態(tài)分布，則繼續(xù)取下一時(shí)刻瓦斯?jié)舛鹊臅r(shí)間序列X_(s+2)執(zhí)行步驟2至步驟4；若該瓦斯?jié)舛鹊臅r(shí)間序列X_(s+1)仍不服從正態(tài)分布，則判定煤礦掘進(jìn)工作面瓦斯涌出已進(jìn)入異常階段，即瓦斯涌出異常時(shí)刻處于t＝s-1時(shí)刻和t＝s時(shí)刻之間，在t＝s時(shí)刻掘進(jìn)工作面已進(jìn)入瓦斯涌出異常階段，即t＝s時(shí)刻是瓦斯涌出異常開始點(diǎn)。

所述步驟2包括：

步驟2-1、將煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列進(jìn)行劃分：隨著從煤礦安全生產(chǎn)監(jiān)測監(jiān)控系統(tǒng)中獲得的瓦斯?jié)舛缺O(jiān)測數(shù)據(jù)的動態(tài)更新，分別生成新的瓦斯?jié)舛鹊臅r(shí)間序列X_(t+1)，X_(t+2)，X_(t+3)，……，長度依次為T₁＝n+1，T₂＝n+2，T₃＝n+3，……；

各瓦斯?jié)舛葧r(shí)間序列的相關(guān)信息詳見表1：

表1 Lilliefors檢驗(yàn)各時(shí)間序列信息表

步驟2-2、對煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列X_(t)采用置信水平1-α(這里，α為顯著性水平)進(jìn)行Lilliefors正態(tài)分布檢驗(yàn)。

所述步驟2-2包括：

步驟2-2-1、對煤礦掘進(jìn)工作面瓦斯?jié)舛葧r(shí)間序列X_(t)所包含的瓦斯?jié)舛葦?shù)據(jù)x₁，x₂，…，x_n進(jìn)行標(biāo)準(zhǔn)化變換：令瓦斯?jié)舛葦?shù)據(jù)的期望瓦斯?jié)舛葦?shù)據(jù)的標(biāo)準(zhǔn)差標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)為假設(shè)經(jīng)過標(biāo)準(zhǔn)化變換后的瓦斯?jié)舛葦?shù)據(jù)z＝{z_t，t＝1，2，3…，n}服從標(biāo)準(zhǔn)正態(tài)分布，其數(shù)學(xué)語言表達(dá)為：H₀：z～N(0，1)；

步驟2-2-2、由標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)z₁，z₂，…，z_n的順序統(tǒng)計(jì)量z₍₁₎，z₍₂₎，…，z_(n)，求標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)的經(jīng)驗(yàn)分布函數(shù)F_n(z)：

式中，f_t為標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)z∈[z_(j)，z_(j+1))的頻數(shù)，這里j表示標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)z₁，z₂，…，z_n的順序統(tǒng)計(jì)量z₍₁₎，z₍₂₎，…，z_(n)的下標(biāo)，其中m為標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)z∈[z₍₁₎，z_(j))的頻數(shù)；

步驟2-2-3、標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)的總體分布函數(shù)F₀(z)為：

步驟2-2-4、根據(jù)標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)計(jì)算檢驗(yàn)統(tǒng)計(jì)量D_n的值：

這里，j＝1，2，…，n，約定F_n(z_(n+1))＝1；

步驟2-2-5、作出檢驗(yàn)決策：

在置信水平1-α下，檢驗(yàn)的決策法則如下：

當(dāng)D_n≥D_n，1-α時(shí)，拒絕H₀，即認(rèn)為F_n(z)≠F₀(z)；

當(dāng)D_n＜D_{n，1-α}時(shí)，接受H₀，即認(rèn)為F_n(z)＝F₀(z)；

其中，Lilliefors分布的上側(cè)α分位數(shù)D_n，1-α的近似計(jì)算公式為：

所述步驟3包括：

步驟3-1、將煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列進(jìn)行劃分：隨著從煤礦安全生產(chǎn)監(jiān)測監(jiān)控系統(tǒng)中獲得的瓦斯?jié)舛缺O(jiān)測數(shù)據(jù)的動態(tài)更新，分別生成新的瓦斯?jié)舛葧r(shí)間序列X_(t+1)，X_(t+2)，X_(t+3)，……，其時(shí)間序列的長度均為n，在該步驟中n取8≤n≤50，各時(shí)間序列的信息詳見表2。

表2 Shapiro-Wilk檢驗(yàn)各時(shí)間序列信息表

步驟3-2、對煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列X_(t)采用置信水平1-α(這里，α為顯著性水平)進(jìn)行Shapiro-Wilk正態(tài)分布檢驗(yàn)。

所述步驟3-2，包括：

步驟3-2-1、由瓦斯?jié)舛鹊臅r(shí)間序列X_(t)所包含瓦斯?jié)舛葦?shù)據(jù)x₁，x₂，…，x_n的順序統(tǒng)計(jì)量x₍₁₎，x₍₂₎，…，x_(n)計(jì)算檢驗(yàn)統(tǒng)計(jì)量W：

式中，j表示瓦斯?jié)舛葦?shù)據(jù)x₁，x₂，…，x_n的順序統(tǒng)計(jì)量x₍₁₎，x₍₂₎，…，x_(n)的下標(biāo)，n為時(shí)間序列長度，L為Shapiro-Wilk檢驗(yàn)的順序統(tǒng)計(jì)量的回歸直線斜率的最佳無偏估計(jì)值與一個(gè)使線性系數(shù)標(biāo)準(zhǔn)化的常數(shù)的積，L的計(jì)算公式如下：

式中，a_j為與j相關(guān)的系數(shù)，其值由表3給出；

表3 Shapiro-Wilk檢驗(yàn)所用的系數(shù)a_j

步驟3-2-2、作出檢驗(yàn)決策，在置信水平1-α下，檢驗(yàn)的決策法則如下：

當(dāng)W≤W_α，則瓦斯?jié)舛鹊臅r(shí)間序列X_(t)拒絕正態(tài)性假設(shè)；

當(dāng)W_α＜W≤1，則瓦斯?jié)舛鹊臅r(shí)間序列X_(t)接受正態(tài)性假設(shè)；

其中，W_α是樣本容量為n時(shí)W分布的α分位數(shù)，其取值由表4給出。

表4 Shapiro-Wilk檢驗(yàn)統(tǒng)計(jì)量W的α分位數(shù)表

所述步驟4中，瓦斯?jié)舛鹊臅r(shí)間序列X_(t)是否服從正態(tài)分布的判斷標(biāo)準(zhǔn)是：

若瓦斯?jié)舛葧r(shí)間序列X_(t)所包含數(shù)據(jù)的正態(tài)分布檢驗(yàn)結(jié)果為D_n＜D_n，1-α且W_α＜W≤1，則該瓦斯?jié)舛鹊臅r(shí)間序列X_(t)服從正態(tài)分布，否則，該瓦斯?jié)舛葧r(shí)間序列X_(t)不服從正態(tài)分布，其中，D_n為標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)的Lilliefors檢驗(yàn)統(tǒng)計(jì)量，D_n，1-α為Lilliefors分布的上側(cè)α分位數(shù)，W為瓦斯?jié)舛葦?shù)據(jù)的Shapiro-Wilk檢驗(yàn)統(tǒng)計(jì)量，W_α為Shapiro-Wilk分布的α分位數(shù)。

有益效果：

本發(fā)明基于正態(tài)分布檢驗(yàn)辨識瓦斯涌出異常的瓦斯災(zāi)害預(yù)警方法，具有以下方面的優(yōu)點(diǎn)：對事前危險(xiǎn)狀態(tài)進(jìn)行實(shí)時(shí)監(jiān)控，通過發(fā)現(xiàn)工作面瓦斯涌出的異常點(diǎn)，進(jìn)而確定工作面前方應(yīng)力、瓦斯壓力和強(qiáng)度等一些影響工作面瓦斯?jié)舛鹊囊蛩匕l(fā)生著變化，以此發(fā)出預(yù)警信息，使得瓦斯災(zāi)害預(yù)警從事后向事前轉(zhuǎn)移，以主動式預(yù)防取代傳統(tǒng)被動式的事故防治，有利于瓦斯事故預(yù)防和控制。

附圖說明

圖1是本發(fā)明一實(shí)施例的煤礦掘進(jìn)工作面瓦斯?jié)舛葧r(shí)間曲線圖；

圖2是本發(fā)明一實(shí)施例的基于正態(tài)分布檢驗(yàn)辨識瓦斯涌出異常的瓦斯災(zāi)害預(yù)警方法流程圖；

圖3是本發(fā)明一實(shí)施例的瓦斯?jié)舛葧r(shí)間序列X₍₅₀₎的經(jīng)驗(yàn)分布函數(shù)和總體分布函數(shù)圖；

圖4是本發(fā)明一實(shí)施例的瓦斯?jié)舛葧r(shí)間序列X₍₅₀₎的回歸直線圖。

具體實(shí)施方式

下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例?；诒景l(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒有作出創(chuàng)造性勞動前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。

根據(jù)文獻(xiàn)資料，2005年9月23日研究區(qū)淮南潘三礦東四8號煤運(yùn)輸上山綜掘時(shí)瓦斯突出前后一段時(shí)間內(nèi)的瓦斯?jié)舛葦?shù)據(jù)(采樣周期為一分鐘，共150個(gè)瓦斯?jié)舛葦?shù)據(jù))，繪制瓦斯?jié)舛葧r(shí)間曲線如圖1所示。以研究區(qū)淮南潘三礦東四8號煤運(yùn)輸上山掘進(jìn)工作面瓦斯災(zāi)害孕育時(shí)期為例，基于正態(tài)分布檢驗(yàn)辨識瓦斯涌出異常的瓦斯災(zāi)害預(yù)警方法的流程如圖2所示，包括如下步驟：

步驟1、從潘三煤礦安全生產(chǎn)監(jiān)測監(jiān)控系統(tǒng)中獲得煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列X₍₅₀₎{x_t，t＝1，2，…，50}，其時(shí)間序列長度為T₀＝50，t表示采樣時(shí)刻，x_t表示t時(shí)刻的瓦斯?jié)舛取?/p>

步驟2、對煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列X₍₅₀₎進(jìn)行Lilliefors正態(tài)分布檢驗(yàn)。

步驟2-1、將煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列進(jìn)行劃分：隨著從煤礦安全生產(chǎn)監(jiān)測監(jiān)控系統(tǒng)中獲得的監(jiān)測數(shù)據(jù)動態(tài)更新，分別生成新的瓦斯?jié)舛鹊臅r(shí)間序列X₍₅₁₎，X₍₅₂₎，X₅₃₎，……，其時(shí)間序列的長度依次為T₁＝51，T₂＝52，T₃＝53，……各瓦斯?jié)舛葧r(shí)間序列的相關(guān)信息詳見表5。

表5 Lilliefors正態(tài)分布檢驗(yàn)各時(shí)間序列信息表

步驟2-2、對煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列X₍₅₀₎采用1-α＝0.95的置信水平進(jìn)行Lilliefors正態(tài)分布檢驗(yàn)，具體檢驗(yàn)過程如下所述：

步驟2-2-1、對煤礦掘進(jìn)工作面瓦斯?jié)舛葧r(shí)間序列X₍₅₀₎所包含的瓦斯?jié)舛葦?shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化變換：令瓦斯?jié)舛葦?shù)據(jù)的期望瓦斯?jié)舛葦?shù)據(jù)的標(biāo)準(zhǔn)差標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)為假設(shè)經(jīng)過標(biāo)準(zhǔn)化變換后的瓦斯?jié)舛葦?shù)據(jù)z＝{z_t，t＝1，2，…，50}服從標(biāo)準(zhǔn)正態(tài)分布，其數(shù)學(xué)語言表達(dá)為：H₀：z～N(0，1)，進(jìn)而進(jìn)行Lilliefors正態(tài)分布檢驗(yàn)，并查看檢驗(yàn)結(jié)論是否接受檢驗(yàn)假設(shè)H₀。

步驟2-2-2、由標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)z₁，z₂，…，z₅₀的順序統(tǒng)計(jì)量z₍₁₎，z₍₂₎，…，z₍₅₀₎，求標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)的經(jīng)驗(yàn)分布函數(shù)F₅₀(z)：

式中，f_t為標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)z∈[z_(j)，z_(j+1))的頻數(shù)，這里j表示標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)z₁，z₂，…，z₅₀的順序統(tǒng)計(jì)量z₍₁₎，z₍₂₎，…，z₍₅₀₎的下標(biāo)，其中m為標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)z∈[z₍₁₎，z_(j))的頻數(shù)。

步驟2-2-3、標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)的總體分布函數(shù)F₀(z)為：

步驟2-2-4、根據(jù)標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)計(jì)算檢驗(yàn)統(tǒng)計(jì)量D₅₀的值：

式中，j＝1，2，…，50，約定F₅₀(z₍₅₀₊₁₎)＝1。

標(biāo)準(zhǔn)化后的瓦斯?jié)舛葦?shù)據(jù)的經(jīng)驗(yàn)分布函數(shù)和總體分布函數(shù)如圖3所示。

步驟2-2-5、作出檢驗(yàn)決策：

在置信水平1-α＝0.95下，檢驗(yàn)的決策法則如下：

當(dāng)D₅₀≥D_{50，0.95}時(shí)，拒絕H₀，即認(rèn)為F₅₀(z)≠F₀(z)；

當(dāng)D₅₀＜D_{50，0.95}時(shí)，接受H₀，即認(rèn)為F₅₀(z)＝F₀(z)；

其中，D_50，0.95為樣本容量為50時(shí)，Lilliefors分布的上側(cè)α＝0.05分位數(shù)的值，其近似計(jì)算公式為：

瓦斯?jié)舛葧r(shí)間序列X₍₅₀₎的Lilliefors檢驗(yàn)結(jié)果為D₅₀＝0.1110＜D_{50，0.95}。

步驟3、對煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列X₍₅₀₎進(jìn)行Shapiro-Wilk正態(tài)分布檢驗(yàn)。

步驟3-1、將煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列進(jìn)行劃分：隨著從煤礦安全生產(chǎn)監(jiān)測監(jiān)控系統(tǒng)中獲得的監(jiān)測數(shù)據(jù)動態(tài)更新，分別生成新的瓦斯?jié)舛葧r(shí)間序列X₍₅₁₎，X₍₅₂₎，X₍₅₃₎，……，其時(shí)間序列的長度均為T＝50，各瓦斯?jié)舛葧r(shí)間序列的相關(guān)信息詳見表6。

表6 Shapiro-Wilk正態(tài)分布檢驗(yàn)各時(shí)間序列信息表

步驟3-2、對煤礦掘進(jìn)工作面瓦斯?jié)舛鹊臅r(shí)間序列X₍₅₀₎采用1-α＝0.95的置信水平進(jìn)行Shapiro-Wilk正態(tài)分布檢驗(yàn)，具體檢驗(yàn)過程如下所述：

步驟3-2-1、由瓦斯?jié)舛葧r(shí)間序列X₍₅₀₎所包含瓦斯?jié)舛葦?shù)據(jù)x₁，x₂，…，x₅₀的順序統(tǒng)計(jì)量x₍₁₎，x₍₂₎，…，x₍₅₀₎計(jì)算檢驗(yàn)統(tǒng)計(jì)量W：

式中，j表示瓦斯?jié)舛葦?shù)據(jù)x₁，x₂，…，x₅₀的順序統(tǒng)計(jì)量x₍₁₎，x₍₂₎，…，x₍₅₀₎的下標(biāo)，時(shí)間序列長度為50，L為Shapiro-Wilk檢驗(yàn)的順序統(tǒng)計(jì)量的回歸直線斜率的最佳無偏估計(jì)值與一個(gè)使線性系數(shù)標(biāo)準(zhǔn)化的常數(shù)的積，L的計(jì)算公式如下：

式中，α_j為與順序統(tǒng)計(jì)量x₍₁₎，x₍₂₎，…，x₍₅₀₎相關(guān)的系數(shù)，其值由表3給出。

時(shí)間序列X₍₅₀₎所包含數(shù)據(jù)x₁，x₂，…，x₅₀的回歸直線如圖4所示。

步驟3-2-2、作出檢驗(yàn)決策，在置信水平1-α＝0.95下，檢驗(yàn)的決策法則如下：

當(dāng)W≤W_0.05，則瓦斯?jié)舛鹊臅r(shí)間序列X_(t)拒絕正態(tài)性假設(shè)；

當(dāng)W_0.05＜W≤1，則瓦斯?jié)舛鹊臅r(shí)間序列X_(t)接受正態(tài)性假設(shè)；

這里，W_0.05為n＝50時(shí)W分布的α＝0.05分位數(shù)，其值由表4給出，在該步驟中W_0.05＝0.947。

時(shí)間序列X₍₅₀₎的Shapiro-Wilk檢驗(yàn)結(jié)果為W_0.05＜W＝0.9620≤1。

步驟4、經(jīng)過Lilliefors正態(tài)分布檢驗(yàn)和Shapiro-Wilk正態(tài)分布檢驗(yàn)，瓦斯?jié)舛葧r(shí)間序列X₍₅₀₎的檢驗(yàn)結(jié)果為D₅₀＜D_50，0.95且W_0.05＜W≤1，則該瓦斯?jié)舛葧r(shí)間序列服從正態(tài)分布，煤礦掘進(jìn)工作面瓦斯涌出處于正常階段，取下一時(shí)刻更新后的瓦斯?jié)舛葧r(shí)間序列X₍₅₁₎，執(zhí)行步驟2至步驟4；如果瓦斯?jié)舛葧r(shí)間序列X₍₅₁₎的檢驗(yàn)結(jié)果為不服從正態(tài)分布，則執(zhí)行步驟5。

取下一時(shí)刻更新后的瓦斯?jié)舛葧r(shí)間序列X₍₅₁₎，進(jìn)行如步驟2至步驟4，經(jīng)檢驗(yàn)，該瓦斯?jié)舛葧r(shí)間序列服從正態(tài)分布，繼續(xù)取再下一時(shí)刻更新后的瓦斯?jié)舛葧r(shí)間序列X₍₅₂₎，并進(jìn)行如步驟2至步驟4所述的正態(tài)分布檢驗(yàn)，依次類推，直至在t＝79時(shí)刻獲得更新后的瓦斯?jié)舛葧r(shí)間序列X₍₇₉₎不服從正態(tài)分布，瓦斯?jié)舛葧r(shí)序序列X₍₅₀₎至X₍₇₉₎的檢驗(yàn)結(jié)果如表7所示：

表7時(shí)序序列X₍₅₀₎至X₍₇₉₎的檢驗(yàn)結(jié)果

步驟5、取下一時(shí)刻的瓦斯?jié)舛葧r(shí)間序列X₍₈₀₎，對X₍₈₀₎執(zhí)行步驟2至步驟4，若該瓦斯?jié)舛葧r(shí)間序列服從正態(tài)分布，則繼續(xù)取下一時(shí)刻的瓦斯?jié)舛葧r(shí)間序列X₍₈₁₎執(zhí)行步驟2至步驟4；若該瓦斯?jié)舛葧r(shí)間序列仍不服從正態(tài)分布，則判定煤礦掘進(jìn)工作面瓦斯涌出已進(jìn)入異常階段，即異常時(shí)刻處于t＝78時(shí)刻和t＝79時(shí)刻之間，在t＝79時(shí)刻掘進(jìn)工作面已進(jìn)入瓦斯涌出異常階段，即t＝79時(shí)刻是瓦斯異常開始點(diǎn)，此時(shí)工作人員要提高警惕并排查安全隱患。

為了排除瓦斯?jié)舛葧r(shí)間序列X₍₇₉₎不服從正態(tài)分布的情況是由瓦斯?jié)舛葦?shù)據(jù)異常造成的，繼續(xù)再取t＝80時(shí)刻獲得的瓦斯?jié)舛葧r(shí)序序列X₍₈₀₎，對X₍₈₀₎進(jìn)行如步驟2至步驟4所述的正態(tài)分布檢驗(yàn)，其檢驗(yàn)結(jié)果為D₈₀＝0.0966＜D_{80，0.95}，W＝0.9349＜W_0.05，所以該瓦斯?jié)舛葧r(shí)間序列仍不服從正態(tài)分布，則判定掘進(jìn)工作面瓦斯涌出已進(jìn)入異常階段，即異常時(shí)刻處于t＝78和t＝79之間，認(rèn)為在t＝79時(shí)刻掘進(jìn)工作面已進(jìn)入瓦斯涌出異常階段，即t＝79時(shí)刻是瓦斯異常開始點(diǎn)，此時(shí)工作人員要提高警惕并排查安全隱患。