周 偉

(霍州煤電集團(tuán)有限責(zé)任公司團(tuán)柏煤礦，山西 霍州 031400)

1 前言

工作面回采破壞了煤層頂板原巖應(yīng)力平衡狀態(tài)，改變了巖體內(nèi)部的應(yīng)力場環(huán)境，使巖體內(nèi)應(yīng)力重新分布，最終會形成一個新的應(yīng)力平衡狀態(tài)[1]。頂板巖層結(jié)構(gòu)運(yùn)動變化是影響采場礦壓顯現(xiàn)的主要原因，科學(xué)掌握工作面回采過程中頂板運(yùn)動變化特征，掌握礦壓顯現(xiàn)規(guī)律，對于維護(hù)采場穩(wěn)定、保障工作面安全生產(chǎn)具有重要的作用。

眾多專家學(xué)者從采場圍巖控制以及頂板管理的角度出發(fā)，基本掌握了采場頂板破斷演化規(guī)律，比較成熟的有壓力拱假說、懸臂梁假說、鉸接巖塊假說、傳遞巖梁理論、砌體梁理論以及薄板理論等[2～3]。這些理論在不同程度上都揭示了工作面回采過程中頂板破壞演化特征。團(tuán)柏煤礦11- 101工作面采用綜采走向長壁采煤法。由于工作面還未開采，為給實際生產(chǎn)提供技術(shù)指導(dǎo)，本文采用F- RFPA2D數(shù)值模擬工作面回采過程中頂板破壞演化規(guī)律。

2 工程概況

團(tuán)柏煤礦井田西部位于呂梁山南端的東麓，隸屬于山西焦煤霍州煤電集團(tuán)，井田批準(zhǔn)開采9號、10號、11號煤。其中11號煤埋深約304～381m，煤層厚度為3～3.5m，平均厚度約為3.3m，傾角1°～9°，平均4°。煤層直接頂厚約5.4m，以粉砂巖為主，層面可見泥質(zhì)成份，具裂隙，裂隙呈半充填；基本頂厚約9.7m，以灰?guī)r為主，中部夾有0.4m泥巖；直接底厚約1.5m，以鋁質(zhì)泥巖為主，層間有少量巖屑。本文以11- 101工作面為研究對象，該工作面位于11號煤首采區(qū)右翼，走向長度約為850m，傾向長度約為135m。

3 數(shù)值分析

巖石破裂全過程分析系統(tǒng)F- RFPA2D以有限元理論為基礎(chǔ)，根據(jù)計算力學(xué)原理，通過弱化單元模擬材料的非線性變形和巖石破裂過程的非連續(xù)行為，采用計算機(jī)可視化技術(shù)展現(xiàn)材料的破裂過程。作為彈性力學(xué)分析工具，根據(jù)彈性損傷理論及修正的Mohr- Coulomb破壞準(zhǔn)則研究巖石材料從細(xì)觀損傷到宏觀破壞的全部過程，對解決實際工程問題具有參考意義[4～5]。

3.1 模型建立

根據(jù)11- 101工作面的地質(zhì)條件，建立數(shù)值模型，模型大小為250m×200m。根據(jù)表1巖石力學(xué)性質(zhì)參數(shù)實驗數(shù)據(jù)，對相關(guān)巖層賦值。考慮到計算精度，將模型劃分250mm×200mm共50 000個單元。為了降低計算量，模型未建到地表，按海姆假說，在模型上部加200m的上覆巖層，模擬上覆巖層載荷大小為5MPa。處理模型邊界時，模型兩側(cè)和底部邊界均為滑動支座，兩側(cè)限制水平和垂直方向位移，底部限制垂直方向位移；底部邊界與兩側(cè)邊界的角點位置采用固定支座，限制模型水平方向與垂直方向位移，如圖1所示。

表1 巖石力學(xué)性質(zhì)參數(shù)

圖1 數(shù)值模擬力學(xué)模型示意圖

3.2 模擬過程

11號煤層采動之前，巖體處于原巖應(yīng)力，所以數(shù)值計算時，首先要初始化模型應(yīng)力。煤層采動后，覆巖遭到破壞，原巖應(yīng)力受到擾動，應(yīng)力重新分布，采用該軟件中修正的帶拉伸截斷的庫侖準(zhǔn)則，根據(jù)11- 101工作面的推進(jìn)方向和推進(jìn)速度開挖11號煤層，每次開挖步距為6m，開挖20步，模擬11號煤層采動，分析頂巖層在工作面推進(jìn)過程中的破壞特征。

3.3 模擬結(jié)果

數(shù)值模擬的重點是分頂板的初次來壓步距、周期來壓步距和在不同來壓步距中頂板巖層的破壞特征。圖片中的灰度值反映圍巖應(yīng)力大小，灰度越亮，應(yīng)力越大，反之，越來越小。

(1)直接頂垮落。如圖2所示，當(dāng)工作面推進(jìn)到12m位置時，原巖應(yīng)力受到擾動較大，在切眼和煤壁前方位置出現(xiàn)應(yīng)力集中，上覆巖層位置出現(xiàn)一個“壓力拱”，前拱腳位于工作面前方實體煤，后拱腳位于切眼后方實體煤。由于直接頂巖性較弱，當(dāng)推進(jìn)到一定距離時，由于懸露面積增加，隨著計算步數(shù)的增加，由于沒有支撐，受到重力作用，直接頂發(fā)生彎曲，然后離層，最終發(fā)生垮落，其垮落高度為2m。

(2)基本頂初次來壓。如圖3所示，隨著工作面不斷向前推進(jìn)，直接頂不斷垮落，垮落高度保持在2m位置。當(dāng)工作面推進(jìn)到18m位置時，采空區(qū)上方“壓力拱”的拱頂向上部延伸，工作面前方煤壁的應(yīng)力不斷增大，應(yīng)力集中系數(shù)約為2.5。

圖2 直接頂垮落模擬結(jié)果(推進(jìn)12m)

圖3 基本頂初次來壓模擬結(jié)果

當(dāng)工作面推進(jìn)到24m時，從模擬結(jié)果中可以看出，位于工作面切眼一側(cè)的基本頂巖層發(fā)生嚴(yán)重破壞。由于巖層自然平衡的結(jié)果，在發(fā)生塑性變形、破壞的巖體周邊，形成了一個“壓力拱”結(jié)構(gòu)，在前后拱腳之間，由于巖層結(jié)構(gòu)破壞，應(yīng)力得到釋放，形成了減壓區(qū)。應(yīng)力場發(fā)生變化的范圍隨著工作面的向前推進(jìn)逐漸增大。

當(dāng)工作面推進(jìn)到36m位置時，懸露的基本頂承受的載荷超過巖層抗拉強(qiáng)度，基本頂發(fā)生變形斷裂、回轉(zhuǎn)，導(dǎo)致工作面前方煤壁壓力突然增大，基本頂發(fā)生初次來壓。基本頂破斷后，積聚的能量得到釋放，煤壁前方承受的壓應(yīng)力變小。所以，根據(jù)數(shù)值模擬結(jié)果，11- 101綜采工作面的基本初次來壓步距約為36m。

(3)基本頂周期來壓。如圖4所示，在基本頂發(fā)生初次來壓后，上覆巖層趨于穩(wěn)定，工作面向前推進(jìn)，基本頂會再次發(fā)生斷裂，結(jié)構(gòu)失穩(wěn)。當(dāng)工作面推進(jìn)到54m位置時，巖層結(jié)構(gòu)的失穩(wěn)導(dǎo)致了工作面頂板來壓，這種來壓將隨工作面的推進(jìn)而呈周期性出現(xiàn)，即所謂的工作面頂板周期性來壓。根據(jù)數(shù)值模擬結(jié)果，基本頂?shù)牡谝淮沃芷趤韷翰骄嗉s為18m。

圖4 基本頂周期來壓模擬結(jié)果

當(dāng)工作面推進(jìn)到78m位置時，上覆巖層結(jié)構(gòu)發(fā)生了很大的變化，已經(jīng)發(fā)生斷裂、離層的巖塊在自重和開采擾動等因素的綜合影響下，部分發(fā)生閉合。由于模型未建立到地表，覆巖破壞已經(jīng)波及模型頂部邊界，基本頂發(fā)生第二次明顯的周期來壓。所以，根據(jù)數(shù)值模擬結(jié)果，第二次周期來壓步距約為24m。綜合以上結(jié)果，11- 101工作面基本頂周期來壓步距約為18～24m。

4 結(jié)論

根據(jù)團(tuán)柏井田11- 101綜采工作面的實際地質(zhì)條件，采用F- RFPA2D模擬軟件，對工作面回采過程中頂板破壞演化規(guī)律展開數(shù)值模擬研究，主要得到以下結(jié)論。

(1)隨著工作面的向前推進(jìn)，該工作面頂板覆巖結(jié)構(gòu)變化為“壓力拱”，前拱腳位于工作面前方煤體，后拱腳位于開切眼后方煤體，拱頂隨著工作面推進(jìn)，逐漸延伸到模型上部邊界。

(2)根據(jù)數(shù)值模擬結(jié)果，11- 101工作面直接頂初次垮落步距約為12m，垮落高度為2m；基本頂初次來壓步距約為36m；基本頂周期來壓步距約為18～24m。數(shù)值模擬結(jié)果對于現(xiàn)場預(yù)測預(yù)報頂板來壓與選擇合理有效支護(hù)方式具有重要的指導(dǎo)作用。

[1] 錢鳴高，繆協(xié)興，許家林，等．巖層控制的關(guān)鍵層理論[M]．徐州：中國礦業(yè)大學(xué)出版社，2003．

[2] 宋占國.煤礦采場礦壓的一般規(guī)律理論研究[J].煤炭技術(shù),2009,28(2):160-163.

[3] 王世炫,羅 勇,沈順平.大采高綜采工作面礦壓觀測及顯現(xiàn)特征研究[J].煤炭技術(shù),2014,33(10):177-180.

[4] 林崇德.層狀巖石頂板破壞機(jī)理數(shù)值模擬過程分析[J].巖石力學(xué)與工程學(xué)報,1999,(4):392-396.

[5] 謝文兵，陳曉翔，鄭百生．采礦工程問題數(shù)值模擬研究與分析[M]．徐州：中國礦業(yè)大學(xué)出版社，2005.