急傾斜煤層開采沖擊危險性數值模擬研究★

于麗艷 蓋芳芳 于月民

(黑龍江科技大學理學院，黑龍江 哈爾濱 150022)

·巖土工程·地基基礎·

急傾斜煤層開采沖擊危險性數值模擬研究★

于麗艷 蓋芳芳 于月民

(黑龍江科技大學理學院，黑龍江 哈爾濱 150022)

對北京昊華能源股份責任公司的長溝峪煤礦的急傾斜煤層進行三維數值模擬，利用ANSYS有限元軟件單元“處死”的功能真實模擬4槽工作面的開挖過程，對沿巷道走向和沿煤層走向推采不同距離時的塑性應變云圖和應力云圖進行分析，得出數值模擬結果與現場實際開采情況基本吻合，說明ANSYS有限元數值模擬方法可為沖擊地壓的防治工作和煤礦的安全開采提供參考。

急傾斜煤層，沖擊危險性，三維數值模擬

0 引言

沖擊地壓[1-4]是全球范圍內巖石工程及煤礦中遇到的最為嚴重的自然災害之一，沖擊地壓發生時煤體坍塌、崩落、巷道底板發生底鼓、頂板下沉、巷道支護物破壞等現象，并造成人員傷亡，嚴重威脅井下安全生產，是煤礦井下開采過程中嚴重的災害之一，是目前國際采礦工程和巖石力學界迫切需要研究解決的科學難題。我國的急傾斜煤層分布較廣泛，在南部地區80%的礦山賦存急傾斜煤層。在地方和鄉鎮煤礦中急傾斜煤層的礦井數約占37.58%，產煤量占7.0%，國有重點煤礦開采急傾斜煤層的礦井數占17%，產煤量占3.88%。隨著采深的增加，由于急傾斜煤層地質條件比較復雜，開采時易發生沖擊地壓現象，這已成為制約我國煤礦安全生產的重大災害之一。鑒于此在開采前結合礦區地質條件用數值模擬方法判斷出煤層的沖擊危險性是很有必要的，這樣可為煤礦的安全順利開采提供保障，為預測和防治沖擊地壓提供依據[5-7]。

1 沖擊危險性的判定原則

沖擊地壓自身規律顯示，沖擊地壓通常發生在應力值比較高且集中的煤巖體上。因此要判斷煤巖體哪個區域的沖擊危險性大，首先要判定煤巖體哪個區域的應力值比較大，煤巖體高應力區域確定了，易于發生沖擊地壓的危險區域也就確定了。當然在應力比較高的區域，并非每處都要發生沖擊地壓，只是此處沖擊危險性較大，至于開采時是否發生沖擊地壓還要根據煤巖體本身的物理力學性質、所處的地質構造特點、開采時的生產條件以及為降低煤巖體應力所采取的措施等有關。

依據生產的技術條件、經驗及地質條件綜合判斷礦井的高應力區[8-11]主要集中在：

1)采空區周圍，工作面前方回采巷道或其他巷道。

2)本層或鄰層的開采邊界或遺留煤柱影響區。

3)地質構造比較復雜處如向斜軸部、斷層、褶曲、煤層突變處、火成巖侵入體兩側。

4)煤礦開采形成的垂直應力集中處，開采支撐壓力增高處，由于應力發生變化，受外界微小干擾而發生失穩破壞，此處為沖擊地壓危險區域。

2 三維數值模擬

2.1 地質概況

采區位于長溝峪向斜南翼,為單斜地層,該區地質構造屬于較簡單類型，煤層的傾角為45°～80°，平均傾角為55°，在向斜軸部附近傾角比較小，為45°左右。采區4槽煤層最小厚度為0.05 m,煤層最大厚度為7 m,平均煤厚3.3 m,平均傾角為55°, 煤層厚度較為穩定,大部分可采，基本為單層結構,局部為雙層結構,地質儲量29.4萬t,可采儲量18.6萬t。底板為粉砂巖,厚度10.7 m；老頂為砂巖,厚度22.0 m; 直接頂為粉砂巖,厚度15.0 m;偽頂為粉砂巖,厚度0.5 m。

2.2 三維數值計算有限元模型

據回采巷道地質、煤巖等條件在回采巷道推進方向建立煤層的三維數值計算有限元模型。為了減少模型邊界效應的影響，在模型垂直方向上取200 m，模型走向方向上取200 m,厚度取100 m。模型左右端面和背面施加平面約束，模型底部端面固定，模型頂部和正面為自由面。取模型上邊界距地表575 m,在模型上部施加均布載荷14.37 MPa，在模型左右端面施加載荷16.78 MPa,并且考慮自重[12-14]。

建立三維數值計算有限元模型，如圖1所示。

2.3 三維數值計算所取參數

數值計算所用的煤巖層的物理、力學參數以力學實驗結果和現場地質調查給定,具體見表1。

表1 三維數值計算參數表

2.4 三維數值計算結果與分析

由圖2可以看出應力集中區都出現在開采工作面的上下端。隨著沿巷道走向逐漸推進，應力集中區域也隨之向里擴展，即工作面上下端頭具有一定的沖擊危險性。

通過圖3可以看出塑性變形出現在煤層與底板交界處和工作面上、下端頭，隨著開采推進，塑性應變出現的范圍在逐漸擴大，最大塑性變形都發生在工作面上下端頭。

由圖4可知,未開采時垂直應力為壓應力，大小分布的總體趨勢均勻，在煤層傾角較大處容易出現應力升高區和應力降低區，如模型傾角大的上下底板處，急傾斜煤層的存在導致應力集中的產生。整個模型應力值都不大，應力最小值出現在底板和煤層邊界處，最小應力值5.73 MPa，相對大的應力值出現在老頂處，值為54.3 MPa，且范圍很小。

由圖5可以看出：

1)雖然沿煤層走向推采距離不同，但應力集中的區域出現在同處，最大的拉應力區出現在工作面下端處，最大壓應力區出現在工作面上端處。

2)推采15 m時壓應力最大值為55.8 MPa，拉應力最大值為10.5 MPa；推采25 m和35 m時應力值減小，推采45 m時應力值又增大為53.5 MPa，8.9 MPa。在推采不同距離時應力值變化不大，原因是在數值計算過程中沒有忽略支護的作用。

3)工作面下巷臨近工作面65 m范圍內和工作面上端頭25 m的范圍內區域應力值會比較大，為高應力區域，此區域沖擊危險性大。

3 結語

1)通過對長溝峪煤礦向斜南翼4槽煤層沿巷道走向和沿煤層走向推采不同距離時的三維Y向塑性應變云圖和三維垂直應力云圖分析可知：在工作面下巷臨近工作面65 m范圍內和工作面上端頭25 m范圍內為高應力集中區，此處沖擊危險性大，是發生沖擊地壓的危險區域。

2)運用ANSYS軟件對長溝峪煤礦急傾斜煤層開采進行了三維數值模擬計算，找到了開采時沖擊危險性大的區域，這一模擬結果與開采實際基本吻合,說明數值模擬可為煤礦的安全開采和沖擊地壓的防治提供借鑒。

3)急傾斜煤層大多數具有沖擊傾向性,為了達到預防沖擊地壓事故發生的目的,在開采前必須進行相關實驗研究,確定煤層沖擊性的強度，根據實驗結果并結合數值模擬結果確定相對危險的區域，開采時采取相應的措施，確保急傾斜煤層的安全開采。

[1] 潘一山.沖擊地壓發生和破壞過程研究[D].北京:清華大學，1999.

[2] 潘一山，徐秉業.沖擊地壓定量預測的研究[J].煤礦開采，1998(3):35-38.

[3] 李忠華，潘一山.采煤工作面沖擊地壓的解析分析[J].遼寧工程技術大學學報，2002，21(1):40-42.

[4] Brauner G.Rock Bursts in Coal Mines and Their Prevention[M].Rotterdam：A.A.Balkema，1994:2-64.

[5] 王世娟,潘一山，李忠華，等.通化八寶礦煤層沖擊傾向性實驗與數值模擬[J].中國地質災害與防治學報,2010,21(3):105-108.

[6] 鄧利民,張 芳，于 華.傾斜煤層沖擊地壓危險狀況的數值模擬研究[J].遼寧工程技術大學學報，2001，20(4)：455- 456.

[7] 關 杰，孫可明，朱月明.急傾斜煤層開采解放層防治沖擊地壓數值模擬[J].遼寧工程技術大學學報，2002，21(4):411-413.

[8] 王永秀，齊慶新.煤柱應力分布規律的數值模擬分析[J].煤炭科學技術，2004，3(10):59-62.

[9] 毛德兵.沖擊礦壓發生危險性評價方法[J].煤礦開采,2000(4)：52-53.

[10] 李英杰，潘一山，唐巨鵬，等.五龍礦沖擊地壓危險區劃分研究[J].礦山壓力與頂板管理，2005(1):94-96.

[11] 于麗艷,潘一山，李忠華，等.吉林省道清礦北斜井南平峒煤層沖擊危險性數值模擬[J].中國地質災害與防治學報, 2011,22(1)：94-98.

[12] 朱伯芳.有限元法原理及應用[M].北京：水電出版社,1979.

[13] 唐巨鵬，潘一山.ANSYS在煤礦開采數值模擬中應用研究[J].巖土力學，2004，25(sup):329-342.

[14] 于麗艷,閆龍海.急傾斜煤層的沖擊地壓危險區域劃分研究[J].山西建筑，2012,38(24)：84-85.

Research of dimensional numerical simulationon rockburst risk of mining steeply inclined coal seam★

YU Li-yan GAI Fang-fang YU Yue-min

(College of Science, Heilongjiang Institute of Science & Technology, Harbin 150022, China)

In the paper, there is a analysis of three dimensional numerical simulation on the steeply inclined coal seam in the Beijing Haohua energy shares the responsibility of the company’s Changgouyu coal mine, based on the function of element is “killed” of ANSYS finite element analysis software to really simulate mining process of 4 slot mining face, there are anlysises on the cloud picture of vertical stress and plastic strain under the different pushed distance along the roadway and the seam of pushing. Results obtained by the numerical simulation basically consistent with the actual situation. So the method of numerical simulation provides the reference for the prevention of rockburst and safely mining of the coal.

steeply inclined coal seam, rockburst risk, three dimensional numerical simulation

1009-6825(2014)31-0070-03

2014-08-30★:國家自然科學基金項目資助(項目編號：11202070)

于麗艷(1979- )，女，碩士，副教授； 蓋芳芳(1981- )，女，博士，講師； 于月民(1972- )，女，碩士，教授

TD324

A