基于PCA-T-S模糊神經(jīng)網(wǎng)絡(luò)底板破壞深度預(yù)測

高 兵,岳 沖,趙 偉

(棗礦(集團)付村煤業(yè)有限公司,山東 濟寧 277605)

隨著煤層淺部資源的減少,煤層開采深度加深,煤層所處地質(zhì)構(gòu)造更加復(fù)雜,各種地下水威脅越發(fā)嚴(yán)重,礦井水突涌時常發(fā)生,影響了礦井的安全生產(chǎn)。采動下的底板破壞深度對底板突水有著重要影響,對工作面安全開采評價有著重大參考價值。我國諸多學(xué)者們從不同方向、采取不同方法對煤層底板破壞深度進行了研究,現(xiàn)行的確定底板破壞深度的方法主要有:《建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采規(guī)程》[1](以下簡稱 “三下”開采規(guī)程)收錄的底板破壞深度計算公式、經(jīng)驗公式[2-3]、現(xiàn)場觀測、數(shù)值模擬[4-6]、神經(jīng)網(wǎng)絡(luò)[7-8]等方法,這些方法各具優(yōu)點,也存在著不足。如經(jīng)驗公式計算雖然簡單,但底板破壞深度受影響因素較多,因此預(yù)測結(jié)果精度低;現(xiàn)場觀測可靠性高但需消耗大量時間、物力及財力,實測過程復(fù)雜;數(shù)值模擬簡單、經(jīng)濟,但受限于參數(shù)的選取,精度較低;神經(jīng)網(wǎng)絡(luò)具有較強的學(xué)習(xí)和預(yù)測能力,在底板破壞深度預(yù)測應(yīng)用廣泛,但未考慮各影響因素之間的相關(guān)性,當(dāng)計算數(shù)據(jù)較多時導(dǎo)致建模效率低。而本文則利用主成分分析與T-S模糊神經(jīng)網(wǎng)絡(luò)組合的方式來預(yù)測煤層底板破壞深度,建立PCA-T-S模糊神經(jīng)網(wǎng)絡(luò)組合預(yù)測模型。

由于底板破壞深度影響因素較多,會導(dǎo)致模型建立過程中維數(shù)過大,較多因素帶有太多重疊信息,導(dǎo)致模型計算效率低,通過采用主成分分析法對輸入數(shù)據(jù)進行預(yù)處理,消除各個因素之間的相關(guān)性,同時降低輸入T-S神經(jīng)網(wǎng)絡(luò)的輸入維數(shù),極大地縮短了模型學(xué)習(xí)時間,再利用T-S神經(jīng)網(wǎng)絡(luò)的模糊性、灰色性的優(yōu)點來預(yù)測底板破壞深度,建立一種預(yù)測精度較高的模型。

1 理論分析

1.1 主成分分析法

主成分分析法是一種數(shù)據(jù)集簡化技術(shù),通過對原始數(shù)據(jù)指標(biāo)信息進行降維,提取出主要指標(biāo)信息,同時減少信息損失[9-10]。經(jīng)分析,所提取的指標(biāo)信息即主成分之間互不相關(guān),每個主成分都是原始變量的線性組合。主成分方差貢獻率越大,說明占比越大。主成分分析步驟如下:

1) 對原始數(shù)據(jù)A進行標(biāo)準(zhǔn)化處理;

2) 計算標(biāo)準(zhǔn)化數(shù)據(jù)的協(xié)方差矩陣Σ的特征向量為λ1≥λ2≥λn,相應(yīng)的特征向量為T1,T2,…Tn.其轉(zhuǎn)換矩陣為A=T′,即A的第i行就是Σ的第i大特征根對應(yīng)的單位特征向量Ti.且第i個主成分Yi的方差就等于Σ的第i大特征根λi;

4) 計算方差貢獻率后,選取的主成分個數(shù)取決于累計方差貢獻率。一般取m個主成分使得方差累積貢獻率達85%以上,說明前m個主成分的樣本信息量包含了n個原始變量所提供的絕大部分信息,即Z1、Z2、Z3、…、Zm.

5) 把標(biāo)準(zhǔn)化后的數(shù)據(jù)帶入主成分表達式,得出所需的新數(shù)據(jù)。

1.2 T-S模糊神經(jīng)網(wǎng)絡(luò)

模糊系統(tǒng)在模糊建模的過程中常存在學(xué)習(xí)能力缺乏、辨識過程復(fù)雜、模型參數(shù)優(yōu)化困難等問題,而人工神經(jīng)網(wǎng)絡(luò)具有自學(xué)習(xí)、自組織和自適應(yīng)的能力,具有強大的非線性處理能力[11]。二者的結(jié)合構(gòu)成的模糊神經(jīng)網(wǎng)絡(luò),可以有效地發(fā)揮模糊邏輯與神經(jīng)網(wǎng)絡(luò)的各自優(yōu)勢,彌補各自的不足。

T-S模糊系統(tǒng)由Takagi和Sugeno在1985年提出[12],是模糊神經(jīng)網(wǎng)絡(luò)形式之一,具有能夠自動更新、不斷修正模糊子集的隸屬函數(shù)的特點,并且擬合性能強、預(yù)測準(zhǔn)確度高。T-S模糊系統(tǒng)用“if-then”規(guī)則形式來定義,在規(guī)則Ri的情況下,推理如下[13]:

(1)

輸入部分為模糊的,輸出部分是確定的,該推理表示輸出是輸入的線性組合。

假設(shè)對于輸入量x=[x1,x2,…,xk],首先根據(jù)模糊規(guī)則計算各輸入變量xj的隸屬度。

(2)

將各隸屬度進行模糊計算,采用模糊算子為連乘算子:

(3)

式中:i=1,2, … ,n;j=1,2,…,k.

根據(jù)模糊計算結(jié)果計算模糊模型的輸出值。

(4)

T-S模糊神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)為輸入層、模糊化層、模糊規(guī)則計算層和輸出層,如圖1所示。輸入層與輸入向量x=[x1,x2, … ,xk]連接,節(jié)點數(shù)與輸入向量的維數(shù)相同。模糊化層采用隸屬度函數(shù)(2)對輸入值進行模糊化得到模糊隸屬度值u,模糊規(guī)則計算采用模糊計算公式(3)計算得到w,輸出層采用公式(4)計算模糊神經(jīng)網(wǎng)絡(luò)輸出。

圖1 T-S模糊神經(jīng)網(wǎng)絡(luò)基本結(jié)構(gòu)圖

2 底板破壞深度影響因素選取

本文參考相關(guān)文獻,依據(jù)研究成果[7,14-15]選取煤層埋藏深度(X1)、煤層傾角(X2)、煤層開采厚度(X3)、工作面斜長(X4)、底板抗破壞能力(X5)以及斷層的存在性(X6)作為底板破壞深度的影響因素,Y為破壞深度,其中底板抗破壞能力以公式(5)進行計算,斷層存在性以“0”代表不存在,“1”為存在,選取25個工作面數(shù)據(jù)作為訓(xùn)練樣本。

(5)

式中:D為底板抗破壞能力;Rc為巖石抗壓強度;C1為節(jié)理裂隙影響系數(shù);C2為分層厚度影響系數(shù)。

3 底板破壞深度影響因素選取

3.1 主成分分析

通過使用SPSS19.0軟件對表1中的樣本數(shù)據(jù)進行標(biāo)準(zhǔn)化處理,然后進行因子降維分析,計算得到的各成分特征值、方差貢獻率及累計貢獻率如表2所示。

表1 原始數(shù)據(jù)[7]

表2 主成分及其貢獻率

對表2中的數(shù)據(jù)分析,累計方差貢獻率大于85%的主成分有4個,其特征值為2.192、1.401、1.072、0.726,第四個主成分累計方差貢獻率達89.852%,能夠很好地解釋原始變量88.815%的信息,說明前4個主成分的提取能夠很好地表達原始變量的總體信息。根據(jù)所求出的成分矩陣載荷進行計算得到相應(yīng)主成分系數(shù),如表3所示。

表3 主成分線性組合系數(shù)

根據(jù)主成分線性組合系數(shù)得到主成分F1,F2,F3,F4的函數(shù)表達式,即:

F1=0.443 8Z(X1)+0.497 1Z(X2)-

0.392 4Z(X3)+ 0.589Z(X4)+0.233Z(X5)-

0.007 4Z(X6)

F2=0.294 9Z(X1)-0.379 4Z(X2)+

0.356 5Z(X3)+ 0.082Z(X4)+0.625 2Z(X5)-

0.495 1Z(X6)

F3=0.514 8Z(X1)+0.179 6Z(X2)+

0.561 1Z(X3)- 0.174 8Z(X4)+0.041 5Z(X5)+

0.596 9Z(X6)

F4=-0.306 3Z(X1)-0.288 7Z(X2)-

0.296 9Z(X3)+0.042 2Z(X4)+0.609 1Z(X5)+

0.600 9Z(X6)

把標(biāo)準(zhǔn)化原始數(shù)據(jù)代入主成分線性組合公式中,得到新的主成分矩陣。

3.2 T-S模糊神經(jīng)網(wǎng)絡(luò)數(shù)據(jù)預(yù)測

將經(jīng)主成分分析后的新數(shù)據(jù)作為自變量進行底板破壞深度預(yù)測。運用MATLAB-R2014a 對T-S模糊神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練,把F1、F2、F3、F4作為影響因子,輸入數(shù)據(jù)維數(shù)為4,底板破壞深度為預(yù)測因子,輸出數(shù)據(jù)維數(shù)為1,每個輸入層神經(jīng)元語言變量的數(shù)目確定為2,故模糊化層節(jié)點數(shù)為8,模糊規(guī)則數(shù)為16,模糊隸屬度函數(shù)中心和寬度c、b進行隨機初始化。選取前24組數(shù)據(jù)進行訓(xùn)練,后8組作為預(yù)測數(shù)據(jù),并與實際數(shù)據(jù)進行對比,用于檢驗網(wǎng)絡(luò)的預(yù)測精度和預(yù)測誤差。

3.3 預(yù)測模型檢驗

用樣本后8組數(shù)據(jù)作為測試樣本,對訓(xùn)練好的預(yù)測模型進行計算和檢驗,并與傳統(tǒng)T-S模糊神經(jīng)網(wǎng)絡(luò)及 “三下”指南收錄的計算公式(6)的預(yù)測結(jié)果進行對比,由圖表可知:

h1=0.008 5H+0.166 5α+0.107 9L-4.357 9

(6)

式中:H為采深;α為煤層傾角;L為工作面斜長。

由圖2分析可知,傳統(tǒng)T-S模糊神經(jīng)網(wǎng)絡(luò)的預(yù)測最大相對誤差為19.98%,這是因數(shù)據(jù)之間相關(guān)性所致;“三下”開采規(guī)范收錄的公式預(yù)測最大相對誤差為20.46%;PCA-T-S模糊神經(jīng)網(wǎng)絡(luò)預(yù)測效果較好,最大相對誤差為15.25%.綜上可知,PCA-T-S模糊神經(jīng)網(wǎng)絡(luò)為底板破壞深度研究提供了一種較高精度的預(yù)測方法。

圖2 預(yù)測結(jié)果對比圖

4 結(jié) 語

1) 選取了埋藏深度、煤層傾角、煤層開采厚度、工作面斜長、底板抗破壞能力以及斷層的存在性等6個影響因素作為預(yù)測模型的輸入指標(biāo),運用主成分分析方法降維原理,在保證數(shù)據(jù)信息完整性的條件下,將底板破壞深度6個影響因素降維為4個主成分,消除各因素之間的相關(guān)性。

2) T-S模糊神經(jīng)網(wǎng)絡(luò)擁有擬合性能強、運算速度快、模擬精度高等特點,通過MALAB軟件把經(jīng)主成分分析后的數(shù)據(jù)輸入到T-S模糊神經(jīng)網(wǎng)絡(luò)中,提高了模型預(yù)測結(jié)果的準(zhǔn)確性。

3) 運用PCA-T-S模糊神經(jīng)網(wǎng)絡(luò)對24組數(shù)據(jù)進行訓(xùn)練,然后對8組數(shù)據(jù)進行預(yù)測,通過與傳統(tǒng)T-S模糊神經(jīng)網(wǎng)絡(luò)及“三下”指南收錄的公式進行對比,PCA-T-S模糊神經(jīng)網(wǎng)絡(luò)組合模型預(yù)測準(zhǔn)確性較高,最大預(yù)測相對誤差僅為15.25%,這為底板破壞深度預(yù)測提供了一種較為準(zhǔn)確的預(yù)測方法。