淺埋煤層走向煤柱對支承壓力影響的數(shù)值模擬研究

劉文濤,許家林,朱衛(wèi)兵,王曉振 邊樹興

(1.中國礦業(yè)大學(xué)煤炭資源與安全開采國家重點實驗室,江蘇省徐州市,221008;2.中國礦業(yè)大學(xué)礦業(yè)工程學(xué)院,江蘇省徐州市,221116;3.中國人民解放軍73151部隊,福建省廈門市,361100)

目前煤柱下開采多集中于對礦壓顯現(xiàn)規(guī)律的研究,對工作面超前支承壓力的研究較少,且對傾向煤柱研究的廣泛、深入程度遠(yuǎn)大于走向煤柱,因此,對走向煤柱下超前支承壓力的研究具有很大的必要性。本文結(jié)合典型淺埋煤層條件——神東礦區(qū)大柳塔煤礦活雞兔井21306工作面的實際開采情況,運(yùn)用FLAC3D軟件對工作面超前支承壓力進(jìn)行模擬研究,掌握超前支承壓力在煤柱下的變化規(guī)律,評價其對工作面開采的影響程度,為今后類似條件下工作面的開采提供參考依據(jù)。

1 地質(zhì)條件

神東煤炭分公司大柳塔煤礦活雞兔井三盤區(qū)為上煤層1-2上煤已采,正在開采下煤層1-2煤的主力采區(qū),目前正在開采21306工作面,平均采高4.3 m,傾向長度255.7 m,共安裝151架支架,走向長度2699.3 m,其中前1050 m范圍內(nèi)地表溝谷縱橫,1-2煤層底板標(biāo)高1100～1106 m,起伏不大,屬近水平煤層,采深75～155 m。地面標(biāo)高1173～1246 m,落差最大達(dá)60 m,其余部分地表較為平緩。距21306工作面機(jī)頭173～193 m的位置為98#～110#支架,上部為已采21上305和21上306工作面間保護(hù)煤柱,煤柱高3.2 m,寬20 m,與工作面走向一致,下煤層1-2煤與上煤層1-2上平均間距為30 m,屬近距離煤層(見圖1所示),相互間影響較為明顯。

圖1 21306工作面的開采條件

2 模型的建立

根據(jù)21306工作面的實際地質(zhì)條件,按照工作面的實際寬度,對工作面距切眼200～550 m范圍內(nèi)的開采區(qū)域進(jìn)行建模,模型在XYZ 3個方向劃分138600個單元,產(chǎn)生147112個節(jié)點(見圖2所示)。

圖2 模型塊體劃分

模型塊體根據(jù)現(xiàn)場柱狀進(jìn)行劃分,但由于塊體劃分時受軟件自身的影響,先劃分網(wǎng)格再生成幾何模型,為避免網(wǎng)格的過度集中,保持網(wǎng)格的均一性,對模擬柱狀做了一定的簡化,同時在工作面推進(jìn)方向研究范圍之外各加設(shè)100 m,以抵消模型計算過程中的邊界效應(yīng),實際模型左右邊界分別對應(yīng)距切眼100 m及650 m位置。

計算過程為:生成模型→原巖應(yīng)力計算→開采1-2上煤→平衡→開采1-2煤→結(jié)果輸出。

模型采用摩爾庫侖屈服準(zhǔn)則:

式中:σ1——最大主應(yīng)力;

σ3——最小主應(yīng)力;

c——表土的粘結(jié)力;

φ——表土的摩擦角。

當(dāng)fs<0時,將發(fā)生剪切破壞。建模參數(shù)如表1所示。

表1 物理力學(xué)參數(shù)

3 模擬結(jié)果分析

取6張垂直應(yīng)力分布圖,見圖3所示。首先由圖3(d)、3(e)和3(f)可以看出工作面前端和兩翼為應(yīng)力集中顯著區(qū)段,前部支承壓力由于煤柱的存在呈現(xiàn)一系列的波狀起伏。輸出圖3(f)中工作面前方集中應(yīng)力峰值繪制的曲線如圖4所示。曲線整體上由3個波峰和2個波谷組成:第一個波峰位于機(jī)頭(Y=170 m)與煤柱之間,中部壓力較大但變化幅度較小,形成一個鍋底狀的壓力升高區(qū)段;第二個波峰位于煤柱正下方,在煤柱集中壓力的作用下,峰值遠(yuǎn)大于其他位置,壓力表現(xiàn)最為明顯;第三個波峰位于煤柱與機(jī)尾(Y=430 m)之間略靠近機(jī)尾側(cè),由于走向煤柱產(chǎn)生的集中應(yīng)力影響范圍有限,在煤柱側(cè)向10～30 m之內(nèi),其作用效果就會降低至最小值,并且煤柱影響范圍的近機(jī)尾側(cè)邊界距側(cè)向煤壁較近,使得側(cè)向煤壁承受了很大一部分載荷,工作面前方煤壁承受的載荷大為減小,因此,峰值與其余兩個相比較小,在此不著重研究。

淺埋煤層賦存條件下煤層內(nèi)部原巖應(yīng)力較低,結(jié)合模擬結(jié)果可以知道應(yīng)力為1.5～3.0 MPa,見圖3(a)。上煤層開采后,工作面間20 m保護(hù)煤柱承受較高的集中應(yīng)力,達(dá)到10.6～13.2 MPa(見圖5中1曲線),類似于400～500 m采深時煤巖內(nèi)部的原巖應(yīng)力,受開采地形影響中部略有降低,應(yīng)力集中系數(shù)為4.8～5.3。當(dāng)集中應(yīng)力傳遞至下煤層底板時減至8.5～9.6 MPa(見圖5中2曲線),為上煤層底板的73%～80%,受地形的影響變小,峰值變化較之前有所降低。下煤層開采過程中,對比圖3(d)、圖3(e)與圖3(f)可以看出工作面推進(jìn)至不同的開采位置時,超前支承壓力十分相近,為突出其微弱量變,將各推進(jìn)位置煤柱下、機(jī)頭與煤柱間超前支承壓力峰值統(tǒng)計如圖5所示。

圖3 垂直應(yīng)力分布圖

可以看到,下煤層開采過程中,工作面位于煤柱下的超前支承壓力顯著提高,達(dá)到16.0～16.7 MPa,相當(dāng)于610～640 m采深時煤巖內(nèi)部的原巖應(yīng)力,與下煤層開采前相比,增幅73%～88%,應(yīng)力集中系數(shù)達(dá)到6.0～8.0;與此同時,工作面機(jī)頭與煤柱間煤壁中的集中應(yīng)力為9.3～10.0 MPa,類似于360～380 m采深時煤巖內(nèi)部的原巖應(yīng)力,應(yīng)力集中系數(shù)為3.6～4.3。對比二者,煤柱下應(yīng)力峰值較機(jī)頭與煤柱間增加了6.6～7.2 MPa,增幅66%～77%。受煤柱下高支承壓力的影響,工作面開采此位置時的危險程度較高,更易產(chǎn)生片幫、漏頂現(xiàn)象。

4 試測驗證

圖6記錄了工作面推進(jìn)至距切眼250～550 m位置中班生產(chǎn)過程中的片幫情況。從圖中可以看出,片幫主要分布于煤柱下部及機(jī)頭與煤柱中間部位,且煤柱下片幫頻率更高。據(jù)不完全統(tǒng)計,煤柱下片幫達(dá)到1 m的位置對應(yīng)距切眼距離為404.4 m、461.2 m、467.2 m、475 m及523.2 m,均位于煤柱下部位置;而機(jī)頭與煤柱之間卻沒有發(fā)生嚴(yán)重的片幫情況,片幫深度多為0.2～0.7 m。 通過以上分析可以知道,由于工作面上部煤柱的存在,超前支承壓力會因煤柱下應(yīng)力集中而變得很大,由于煤層埋藏較淺,其作用結(jié)果最大僅相當(dāng)于不足650 m采深時的原巖應(yīng)力狀態(tài),對生產(chǎn)影響較小,僅煤柱下片幫較為嚴(yán)重,但如果煤層間距變小或是采深加大,這種作用的結(jié)果勢必加劇。

圖7 煤壁片幫位置實測記錄圖

5 結(jié)論

(1)由于煤柱的影響,超前支承壓力表現(xiàn)為一系列波狀彎曲,整體上由3個背斜2個向斜組成,背斜峰值大小依次為煤柱下、機(jī)頭與煤柱間、煤柱與機(jī)尾間。

(2)下煤層開采過程中,煤柱下超前支承壓力峰值較其他位置峰值最小增幅為66%～77%;與下煤層開采前相比,煤柱下超前支承壓力峰值增加了73%～88%,表現(xiàn)出上部遺留煤柱對超前支承壓力的顯著影響。

(3)由于煤層埋深較淺,應(yīng)力集中的效果僅相當(dāng)于不足650 m采深時的原巖應(yīng)力狀態(tài),對生產(chǎn)影響較小,如果煤層間距變小或是采深加大,作用的效果勢必加劇。

[1] 牛偉,董鳳寶,孫孝林.大采高工作面煤柱下開采與圍巖控制[J].山東煤炭科技,2003(3)

[2] 佟文旗,魏勝利.近距離煤層煤柱下的開采方法[J].河北能源職業(yè)技術(shù)學(xué)院學(xué)報,2009(2)

[3] 楊光玉,賀興元.局部煤柱下安全采煤技術(shù)[J].煤炭科學(xué)技術(shù),2001(10)

[4] 王曉振,許家林,朱衛(wèi)兵等.走向煤柱對近距離煤層大采高綜采面礦壓影響[J].煤炭科學(xué)技術(shù),2009(2)

[5] 鄭百生,謝文兵,竇林名等.不規(guī)則煤柱作用下工作面開采的三維數(shù)值模擬[J].煤炭學(xué)報,2006(2)

[6] 許清起.上方煤柱下綜采面巷道沖擊礦壓的防治[J].煤礦開采,2004(3)

[7] 肖雪峰.近距離煤層群開采采場支承壓力分布規(guī)律研究[J].華北科技學(xué)院學(xué)報,2009(1)

[8] 陳軼平.綜采工作面超前支護(hù)支承壓力觀測與分析[J].山西大同大學(xué)學(xué)報(自然科學(xué)版),2009(3)

[9] 魯巖,樊勝強(qiáng),鄒喜正等.工作面超前支承壓力分布規(guī)律[J].遼寧工程技術(shù)大學(xué)學(xué)報(自然科學(xué)版),2008(2)

[10] 侯忠杰,張杰.砂土基型淺埋煤層保水煤柱穩(wěn)定性數(shù)值模擬[J].巖石力學(xué)與工程學(xué)報,2005(13)