導(dǎo)語：煤礦人工神經(jīng)網(wǎng)絡(luò)預(yù)測管理論文一文來源于網(wǎng)友上傳，不代表本站觀點(diǎn)，若需要原創(chuàng)文章可咨詢客服老師，歡迎參考。

摘要：應(yīng)用人工神經(jīng)網(wǎng)絡(luò)的基本原理，建立了一個(gè)基于神經(jīng)網(wǎng)絡(luò)的煤礦立井井筒非采動(dòng)破裂的預(yù)測系統(tǒng)，實(shí)現(xiàn)了立井井筒破裂預(yù)測的智能化。最后將神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果與數(shù)值計(jì)算結(jié)果對比，認(rèn)為應(yīng)用人工神經(jīng)網(wǎng)絡(luò)對立井井筒破裂時(shí)間的預(yù)測比較準(zhǔn)確，實(shí)用。

關(guān)鍵詞：立井井筒非采動(dòng)破裂反向傳播網(wǎng)絡(luò)神經(jīng)網(wǎng)絡(luò)預(yù)測數(shù)值模擬

一、煤礦立井發(fā)生破壞問題的提出

徐淮地區(qū)是中國東部主要的煤炭開采基地，其煤炭的生產(chǎn)直接影響著我國煤炭的總產(chǎn)量，在國民經(jīng)濟(jì)建設(shè)中占有重要的地位。然而自20世紀(jì)80年代以來，在我國徐淮地區(qū)（徐州、淮北、淮南）地區(qū)，出現(xiàn)了一種新的礦井破裂災(zāi)害——井筒的非采動(dòng)破裂，即煤礦立井在不受地下采動(dòng)影響的條件下（井筒及其附近的工業(yè)廣場都留有足夠的保護(hù)煤柱），井壁發(fā)生嚴(yán)重變形和破裂，致使提升運(yùn)輸困難。90年代以來，在我國特大型煤炭企業(yè)'''');">企業(yè)兗州礦業(yè)集團(tuán)的9對井筒也先后發(fā)生了破壞(有關(guān)兗州礦區(qū)井筒破裂的基本情況如表1所示)。煤礦豎井是礦山生產(chǎn)運(yùn)輸?shù)难屎硪溃虼素Q井的破壞嚴(yán)重影響了礦山生產(chǎn)，給各煤礦造成了巨大的經(jīng)濟(jì)損失。

為了盡量減少立井的非采動(dòng)破裂所造成的經(jīng)濟(jì)損失，現(xiàn)在各大礦山都對井筒的變形進(jìn)行了預(yù)報(bào)和治理，到目前對井筒破裂的預(yù)報(bào)方法主要有兩種：

1）第一種方法是通過加強(qiáng)對井筒變形的監(jiān)測，以監(jiān)測數(shù)據(jù)為依據(jù)，對井筒的變形進(jìn)行分析，對其破裂進(jìn)行預(yù)報(bào)。

2）第二種方法是新近發(fā)展起來的預(yù)報(bào)方法［1，2，3］，即系統(tǒng)科學(xué)、智能技術(shù)方法，特別是研究非線性復(fù)雜系統(tǒng)的一些方法。其主要內(nèi)容使用系統(tǒng)科學(xué)原理或智能技術(shù)來建立模型的框架，用觀測的實(shí)測資料填充以實(shí)現(xiàn)建模。

本文采用實(shí)際與智能技術(shù)相結(jié)合的預(yù)報(bào)方法，建立井筒破裂的人工神經(jīng)網(wǎng)絡(luò)模型，模仿人腦的運(yùn)行機(jī)制，通過對井筒破壞規(guī)律的學(xué)習(xí)，使網(wǎng)絡(luò)具有根據(jù)特征值對井筒破壞進(jìn)行預(yù)報(bào)的能力，并據(jù)此來推測相關(guān)煤礦的井筒破壞規(guī)律。

表1兗州礦區(qū)井筒破裂的基本情況

Table1BasicsituationofshaftliningfractureinYanzhouCoalMine

礦區(qū)

序號

井筒名稱

竣工時(shí)間

破裂時(shí)間

凈徑/m

外徑/m

施工方法

井壁類型

表土厚度/m

破裂深度/m

破裂情況

兗州

1

鮑店副井

1979.11.26

1995.6.5

8.0

10.2

凍結(jié)法

雙層井壁

148.6

126.9

罐道縫壓實(shí),罐道,管路壓縮彎曲,混凝土表層剝落出現(xiàn)水平裂縫,豎筋彎曲外露

2

鮑店主井

1979.5.14

1995.7.12

6.5

8.5

凍結(jié)法

雙層井壁

148.69

136—144

3

鮑店北風(fēng)井

1979.10.21

1996.8.2

5.0

6.6

凍結(jié)法

雙層井壁

202.56

168.4,180,204

4

鮑店南風(fēng)井

1979.8.1

1996.8.9

凍結(jié)法

雙層井壁

157.92

158.1—159.3

5

興隆莊西風(fēng)井

1976.8

1995.10

5.5

7.4

凍結(jié)法

雙層井壁

183.9

165.5—171.6

6

興隆莊東風(fēng)井

1977.5.31

1997.6.7

5.0

6.4

凍結(jié)法

雙層井壁

176.45

157—180

7

興隆莊主井

1977.8.13

1997.6.23

凍結(jié)法

雙層井壁

189.31

150,184

在未出現(xiàn)嚴(yán)重破裂時(shí)進(jìn)行了治理

8

興隆莊副井

1978.9

1997.6.26

凍結(jié)法

雙層井壁

190.41

154,200

罐道縫壓實(shí),罐道,管路壓縮彎曲,混凝土表層剝落出現(xiàn)水平裂縫,豎筋彎曲外露

9

楊村主井

1984.12

1997.2.29

5.0

6.6

凍結(jié)法

雙層井壁

185.42

176.5,196

10

楊村副井

1985.1.23

1997.12.2

凍結(jié)法

雙層井壁

184.45

160,176,212

在未出現(xiàn)嚴(yán)重破裂時(shí)進(jìn)行了治理

11

楊村北風(fēng)井

1984.10.31

1997.2.4

4.5

5.9

凍結(jié)法

雙層井壁

173.40

179.6,150,156.6

罐道縫壓實(shí),罐道,管路壓縮彎曲,混凝土表層剝落出現(xiàn)水平裂縫,豎筋彎曲外露

二、人工神經(jīng)網(wǎng)絡(luò)的預(yù)報(bào)原理

人工神經(jīng)網(wǎng)絡(luò)實(shí)現(xiàn)井筒破裂的智能預(yù)報(bào)，是通過機(jī)器學(xué)習(xí)的方法［4］對破裂井筒的特征值進(jìn)行抽取，并對已知的井筒破裂規(guī)律進(jìn)行學(xué)習(xí)掌握規(guī)律性，然后運(yùn)用訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)對其他井筒進(jìn)行推理預(yù)測，并據(jù)此對其他井筒的破裂進(jìn)行預(yù)報(bào)（其流程見圖1）。

三、立井井筒破裂影響因素的選取

經(jīng)調(diào)查表明立井井壁破裂的主要原因?yàn)椋涸诿旱V開采過程中新生界底部第四系含水層（底含）的水頭的大幅疏降，使該含水層及上覆土層產(chǎn)生壓縮和變形，且引起地表沉降，在地層發(fā)生變形的過程中對井壁產(chǎn)生垂直向下的附加力，使得立井井壁發(fā)生破裂。

立井井筒破裂礦區(qū)的水文地質(zhì)與工程地質(zhì)條件都具有如下的特點(diǎn)：井筒都穿過第四系深厚表土，其厚度大都在200m左右。土層結(jié)構(gòu)復(fù)雜，但大體上都可分為四個(gè)含水層和三個(gè)隔水層共七個(gè)工程巖組，即由上至下常簡稱為一含、一隔、二含、二隔、三含、三隔和底含（四含）。

通過對立井井筒非采動(dòng)破裂機(jī)理及破裂礦區(qū)的水文地質(zhì)與工程地質(zhì)特點(diǎn)分析選取以下幾個(gè)因素作為影響立井井筒破裂的特征因素：

1、表土層厚度

由于立井井筒非采動(dòng)破裂只發(fā)生在厚沖積層中建成的立井井筒，因此表土層厚度是立井井筒非采動(dòng)破裂現(xiàn)象發(fā)生的必要因素。表土層厚度越大，土層對立井井筒的側(cè)壓力越大，且土層與井筒的相互作用的面積增大，底含沉降時(shí)產(chǎn)生的立井井筒附加力加大，立井井筒發(fā)生破裂的可能性越大。

2、底含厚度

底含厚度決定了立井井筒周圍土層的變形量，且土層變形量直接關(guān)系到立井井筒附加應(yīng)力的大小，因此底含厚度越大，井筒破裂的可能性增大，所以確定底含厚度為立井井筒破裂的主要影響因素。

3、底含水位降速

底含水位降速?zèng)Q定了立井井筒周圍土層變形的速率，從而決定了立井井壁破裂的時(shí)間。底含水頭降速直接決定了立井井筒破裂時(shí)間的大小。

4、井筒外徑

由于在確定的工程地質(zhì)條件下立井井筒外表面積與立井井筒附加力的大小成正比，則確定立井井筒外徑大小為立井井筒破裂的主要因素。

5、井壁厚度

井壁厚度越大，立井井筒的凈截面積越大，立井井壁內(nèi)壁應(yīng)力降低，有利于立井井筒的穩(wěn)定。

四、神經(jīng)網(wǎng)絡(luò)的設(shè)計(jì)與實(shí)現(xiàn)

根據(jù)以上對影響井筒變形的特征因素的選取，選擇反向傳播（backpropagation,BP）神經(jīng)網(wǎng)絡(luò)算法對井筒的破裂規(guī)律進(jìn)行訓(xùn)練，其網(wǎng)絡(luò)為包含兩層隱含層的神經(jīng)網(wǎng)絡(luò)，輸入層、隱含層、輸出層的神經(jīng)元的個(gè)數(shù)分別為5、20、10、1，

表2神經(jīng)網(wǎng)絡(luò)的輸入矢量p及目標(biāo)矢量t

Table2Inputvectorpandtargetvectortoftheneuralnetwork

輸入矢量p

輸出矢量t

表土層厚度(m)

井筒外徑(m)

井壁厚度(m)

底含水位降速（米/年）

底含厚度(m)

井壁破裂時(shí)間（月）

189.31

8.92

1.21

3.764

34.1

192

190.41

10.1

1.3

3.212

30

225

190.41

6.4

0.7

2.988

32.85

241

189.5

7.4

0.95

3.652

29.9

230

148.69

8.5

1

5.196

56.29

194

148.6

10

1

5.262

55.0

187

202.56

6.6

0.8

5.053

59.0

190

185.5

6.4

0.7

7.192

57.72

146

其訓(xùn)練函數(shù)分別采用雙曲正切函數(shù)tansig及線性激活函數(shù)purelin［5］，網(wǎng)絡(luò)學(xué)習(xí)采用的輸入矢量及目標(biāo)矢量如表2所示。

網(wǎng)絡(luò)訓(xùn)練后，其輸出值與期望值之及目標(biāo)誤差如表3所示。

表3BP網(wǎng)絡(luò)對井筒破壞規(guī)律的學(xué)習(xí)

Table3BPnetworklearningofshaftliningfractureregularity

輸入矢量p

目標(biāo)矢量t

目標(biāo)誤差

輸出矢量

a

189.31

8.92

1.21

3.764

34.1

192

0.0001

199.5

190.41

10.1

1.3

3.212

30

225

0.0001

222.3

190.41

6.4

0.7

2.988

32.85

241

0.0001

243.8

189.5

7.4

0.95

3.652

29.9

230

0.0001

221.3

148.69

8.5

1

5.196

56.29

194

0.0001

193.06

148.6

10

1

5.262

55.0

187

0.0001

185.6

202.56

6.6

0.8

5.053

59.0

190

0.0001

188.9

185.5

6.4

0.7

7.192

57.72

146

0.0001

147.23

五、實(shí)例應(yīng)用及與數(shù)值模擬結(jié)果的比較

兗州礦區(qū)楊村煤礦北風(fēng)井井筒表土段厚173.4m，采用凍結(jié)法施工，于1984年竣工，在1997年的檢查中發(fā)現(xiàn)井壁已發(fā)生了破裂，現(xiàn)在用學(xué)習(xí)后的神經(jīng)網(wǎng)絡(luò)對楊村北風(fēng)井的井筒破裂時(shí)間進(jìn)行預(yù)測，預(yù)測結(jié)果如所表4示。

據(jù)上表可以得出有神經(jīng)網(wǎng)絡(luò)預(yù)報(bào)得出的預(yù)測值與實(shí)際的目標(biāo)矢量之間的誤差僅為0.015，該誤差在現(xiàn)場的實(shí)際工作中是可以接受的，這說明由BP網(wǎng)來預(yù)測豎井的非采動(dòng)破壞在實(shí)際工作中是可行的，且行之有效。

表4神經(jīng)網(wǎng)絡(luò)對楊村北風(fēng)井破裂的預(yù)測

Table4NeuralnetworkforecasttothefracturetimeofYangcunnorthshaft

輸入矢量p

目標(biāo)矢量t

輸出矢量

a

誤差

表土層厚度(m)

井筒直徑(m)

井壁厚度(m)

底含水位降速（米/年）

底含厚度(m)

井壁破裂時(shí)間（月）

173.4

5.9

0.7

7.5

65.3

136

138

0.015

根據(jù)兗州礦區(qū)的工程地質(zhì)資料和及楊村立井井筒施工資料建立了立井井筒破裂的幾何計(jì)算模型，采用Flac3D進(jìn)行數(shù)值模擬計(jì)算，模型共19008個(gè)六面體、21600個(gè)結(jié)點(diǎn)。數(shù)值模擬計(jì)算后立井井壁最大z方向應(yīng)力隨底含水頭降變化如圖2所示，底含失水沉降情況下立井井壁發(fā)生破裂時(shí)的底含水頭降為0.8MPa左右，換算為水頭高度等于80m，此時(shí)井壁內(nèi)部的最大應(yīng)力為30MPa，達(dá)到了立井井筒的破裂強(qiáng)度。則立井井筒的破裂時(shí)間T為：

T＝底含水頭高度損失量/底含水位降速

＝(80÷7.5)×12

＝128月

根據(jù)底含水頭降速可得立井井筒破裂的時(shí)間為128月，與神經(jīng)網(wǎng)絡(luò)預(yù)測值相比，相差10月左右，因此可以認(rèn)為神經(jīng)網(wǎng)絡(luò)預(yù)測基本可以用于立井井筒破裂時(shí)間的預(yù)測。

圖2立井井壁最大z方向應(yīng)力隨底含水頭降變化

Fig.2Waterheadvariationinbottomaquifervs.maximumz-directionalstressinshaftlining

六、結(jié)論

采用人工神經(jīng)網(wǎng)絡(luò)方法來預(yù)測井筒的非采動(dòng)破裂，是基于神經(jīng)網(wǎng)絡(luò)所具有的學(xué)習(xí)、和容錯(cuò)功能進(jìn)行的非精確推理，他在一定程度上模擬了人腦的人工智能，對于煤礦豎井破裂的這種受到許多不確定因素的影響的工程現(xiàn)象往往能夠給出比較準(zhǔn)確的預(yù)測。所以應(yīng)用人工神經(jīng)網(wǎng)絡(luò)對井筒破裂進(jìn)行預(yù)測，可有效的預(yù)測井壁的破裂時(shí)間，且實(shí)用性強(qiáng)。但是由于人工神經(jīng)網(wǎng)絡(luò)是在對已發(fā)生破裂的立井井筒進(jìn)行學(xué)習(xí)后進(jìn)行預(yù)測，所以在學(xué)習(xí)范圍外的立井井筒破裂預(yù)測可能會產(chǎn)生比較大的誤差，因此對于超出訓(xùn)練數(shù)據(jù)外的立井井筒破裂時(shí)間預(yù)測應(yīng)結(jié)合其他方法進(jìn)行比較驗(yàn)證。