煤層群條件下高抽巷合理層位布置的數(shù)值模擬研究

李 忠，張慧杰，賈恒義，郝秀明，孟 濤

(1.晉能控股煤業(yè)集團(tuán) 興峪煤業(yè)有限責(zé)任公司， 山西 陽泉 045100； 2.煤炭科學(xué)技術(shù)研究院有限公司 安全分院 北京 100013； 3.煤炭資源高效開采與潔凈利用國家重點(diǎn)實(shí)驗(yàn)室(煤炭科學(xué)研究總院) 北京 100013)

高瓦斯煤層群賦存條件的煤層回采時(shí)，鄰近層卸壓瓦斯易涌入采空區(qū)，造成回采工作面瓦斯超限，威脅煤礦的安全生產(chǎn)[1]. 現(xiàn)有采空區(qū)瓦斯抽采方法較多，如采空區(qū)埋管抽采瓦斯、高位鉆孔抽采瓦斯、留尾巷抽采瓦斯、高抽巷抽采瓦斯等，其中高抽巷具有抽采時(shí)間長(zhǎng)、輻射半徑大、抽采效果好等優(yōu)點(diǎn)，在采空區(qū)瓦斯抽采中應(yīng)用廣泛[2-3]. 高抽巷抽采瓦斯技術(shù)是將高抽巷布置在頂板裂隙帶內(nèi)，鄰近層瓦斯在負(fù)壓作用下通過采動(dòng)裂隙進(jìn)入高抽巷并被抽出。顯然，合理確定高抽巷的位置對(duì)于提高瓦斯抽采效果、延長(zhǎng)高抽巷使用時(shí)間及減少高抽巷工程量都具有重要意義。

學(xué)者們針對(duì)高抽巷的布置形式及層位開展了大量的研究。周華東等[4]研究了綜采工作面初采期局部高抽巷的瓦斯治理效果，提出局部高抽巷分為初采傾向高抽巷段、走向高抽巷段和輔助傾斜高抽巷段3段。王紅勝等[5]研究了外錯(cuò)高抽巷布置方式，實(shí)現(xiàn)一條高抽巷服務(wù)于相鄰兩工作面、提高高抽巷利用效率的目的。肖峻峰等[6]采用非線性大變形程序數(shù)值模擬了采空區(qū)頂板覆巖應(yīng)力分布及裂隙演化規(guī)律，提出走向高抽巷應(yīng)布置在采空區(qū)頂板裂隙帶的中上部。陳遠(yuǎn)峰等[7]利用FLAC3D研究了高抽巷頂?shù)装濉蓭偷南鄬?duì)移近量，獲得了高抽巷的變形破壞規(guī)律，在此基礎(chǔ)上優(yōu)化了高抽巷的空間分布。林海飛等[8]構(gòu)建了走向高抽巷條件下的采空區(qū)瓦斯運(yùn)移模型，模擬了高抽巷與回風(fēng)巷不同平距、與煤層頂板不同垂距條件下的采空區(qū)瓦斯抽采效果。肖峻峰等[9]研究了近距離高瓦斯煤層群傾向高抽巷的優(yōu)化布置問題，得出傾向高抽巷最佳布置位置在豎直方向上15倍采高、水平方向深入工作面約35 m處。張慧杰等[10]提出了基于巖層控制關(guān)鍵層理論確定走向高抽巷合理層位的方法，并提出由沿傾向的煤層卸壓角和采空區(qū)上方“O”形圈的寬度可確定走向高抽巷的水平位置。

本文基于興峪煤礦煤巖賦存及開采條件，利用FLAC3D模擬軟件建立數(shù)值模型，研究了15號(hào)煤層回采后上覆巖層應(yīng)力、位移及損傷破壞特征，探討了高抽巷的合理層位，以期為高抽巷布置參數(shù)的確定提供借鑒。

1 礦井概況

興峪煤礦位于山西省盂縣城南7.8 km處，石坡峪村東北，行政區(qū)劃屬盂縣路家村鎮(zhèn)。井田內(nèi)含煤地層為石炭系上統(tǒng)太原組及二疊系下統(tǒng)山西組，所含煤層由淺至深依次為1、2、3、4、6、8、9、9下、11、12、13和15號(hào)煤層。其中8、9號(hào)煤層為賦存區(qū)大部可采穩(wěn)定煤層，9下號(hào)煤層為不穩(wěn)定局部可采煤層，15號(hào)煤層為全區(qū)穩(wěn)定可采煤層，其余煤層均為不穩(wěn)定的不可采煤層或薄煤線。

目前,該礦主采15號(hào)煤層，該煤層厚6.33～8.90 m，平均7.06 m，上距9下號(hào)煤層78.65～95.50 m，平均87.92 m. 2020年度礦井瓦斯等級(jí)測(cè)定結(jié)果表明，礦井絕對(duì)瓦斯涌出量為98.62 m3/min，礦井相對(duì)瓦斯涌出量為55.23 m3/t. 在回采15號(hào)煤層的過程中上覆煤層卸壓瓦斯大量涌入采空區(qū)，礦井鄰近層瓦斯涌出量大。為此，該礦采用走向高抽巷方式治理鄰近層瓦斯，高抽巷沿11號(hào)煤層頂板布置，與15號(hào)煤層頂板的垂直距離在40～60 m.

2 模型構(gòu)建

根據(jù)興峪煤礦綜合柱狀圖和實(shí)際開采情況建立FLAC3D數(shù)值模型，模擬工作面回采后上覆巖層的應(yīng)力及位移變化。模型的長(zhǎng)、寬和高(x方向、y方向和z方向)分別為320 m、10 m和160 m，工作面回采寬度為220 m，上覆巖層厚度為100 m，下伏巖層厚度為60 m(含15號(hào)煤層厚度)。模型四周及底部采用位移邊界條件約束，頂部為自由邊界，施加未出現(xiàn)在模型中的上覆巖層的重力載荷Q. 上覆巖層平均容重γ取25 kN/m3，模型頂部距地面距離H′取300 m，計(jì)算得到重力載荷Q為7.5 MPa. 模型單元網(wǎng)格劃分見圖1.

圖1 覆巖運(yùn)移分析模型單元網(wǎng)格圖

模型定義材料為Mohr-Coulomb，遵循彈塑性屈服準(zhǔn)則，計(jì)算過程中考慮大變形，模擬計(jì)算所采用的巖體力學(xué)參數(shù)根據(jù)該礦地質(zhì)資料和實(shí)驗(yàn)室測(cè)試力學(xué)試驗(yàn)結(jié)果確定，具體參數(shù)見表1.

3 結(jié)果及分析

3.1 采動(dòng)覆巖應(yīng)力分布特征

圖2和圖3分別為采動(dòng)覆巖水平應(yīng)力和垂向應(yīng)力分布云圖。從圖2中可以看出，水平應(yīng)力分布特征存在明顯的“分層”現(xiàn)象，即不同距離上覆巖層水平應(yīng)力呈現(xiàn)出交替升降的特征，這可能是由于巖體層位的不同導(dǎo)致的采動(dòng)應(yīng)力分布差異。從圖3中可以看出，在兩側(cè)煤柱位置出現(xiàn)了應(yīng)力集中區(qū)，在采空區(qū)內(nèi)部出現(xiàn)了卸壓區(qū)。采空區(qū)中心區(qū)域的卸壓程度較四周小，這是由于上覆巖層垮落壓實(shí)采空區(qū)中心區(qū)域所致。

圖2 采動(dòng)覆巖水平應(yīng)力分布云圖

圖3 采動(dòng)覆巖垂向應(yīng)力分布云圖

為了研究不同距離上覆巖層采動(dòng)應(yīng)力分布的差異性，提取了0、20 m、40 m、50 m、60 m、70 m、80 m和90 m處上覆巖層的水平應(yīng)力、垂直應(yīng)力以及主應(yīng)力差，見圖4，5，6. 從圖中可以看出，水平應(yīng)力、垂直應(yīng)力及主應(yīng)力差均受采動(dòng)影響明顯，在頂板上覆巖層水平區(qū)域內(nèi)存在應(yīng)力集中區(qū)和卸壓區(qū)。隨著上覆巖層與15號(hào)煤層垂直距離的增加，應(yīng)力集中區(qū)的極值應(yīng)力不斷減小，并且極值點(diǎn)自兩側(cè)煤柱位置不斷向采空區(qū)內(nèi)部收緊，這種上覆巖層梯形臺(tái)狀應(yīng)力分布特征將影響裂隙分布。在垂直方向上高抽巷應(yīng)布置在梯形臺(tái)的“腰線”位置，且高抽巷布置層位越高，在水平位置上越應(yīng)該貼近回采工作面中心區(qū)域。卸壓區(qū)主要位于采空區(qū)內(nèi)部核心區(qū)域，且隨著上覆巖層與15號(hào)煤層垂直距離的增加卸壓區(qū)范圍逐漸減小，這也是導(dǎo)致隨距離增加保護(hù)層開采有效保護(hù)范圍不斷減小的原因。

表1 煤巖層力學(xué)參數(shù)表

圖4 頂板上覆巖層不同水平處的水平應(yīng)力分布圖

圖5 頂板上覆巖層不同水平處的垂向應(yīng)力分布圖

3.2 采動(dòng)覆巖位移規(guī)律分析

圖7和圖8分別為采動(dòng)覆巖水平位移云圖和垂直位移云圖。由圖可以看出，不同距離上覆巖層水平位移也呈“分層”分布的現(xiàn)象，與水平應(yīng)力分布特征一致，可以看出應(yīng)力變化是導(dǎo)致位移發(fā)生變化的誘因。垂直位移主要分布在采空區(qū)正上方區(qū)域，且隨著上覆巖層距離的增加位移量不斷減小。

圖8 采動(dòng)覆巖垂向位移分布云圖

為了研究不同距離上覆巖層位移分布的差異性，提取了0 m、20 m、40 m、50 m、60 m、70 m、80 m和90 m處上覆巖層的水平位移和垂直位移，見圖9，圖10. 由圖可以看出，上覆巖層破壞“豎三帶”特征明顯，距離煤層頂板0～20 m處于垮落帶，20～60 m處于裂隙帶，60 m以上為彎曲下沉帶。另外，遠(yuǎn)距離上覆巖層的水平位移大于近距離上覆巖層，而垂向位移又小于近距離上覆巖層。表明遠(yuǎn)距離上覆巖層以水平位移為主，近距離上覆巖層以垂向位移為主，這種位移規(guī)律造成了近距離上覆巖層內(nèi)的穿層裂隙發(fā)育和遠(yuǎn)距離上覆巖層內(nèi)的順層裂隙發(fā)育現(xiàn)象。順層裂隙發(fā)育區(qū)域的瓦斯流動(dòng)狀態(tài)主要為層內(nèi)流動(dòng)，穿層裂隙發(fā)育區(qū)域的瓦斯流動(dòng)狀態(tài)包含煤巖層間瓦斯越流，因而，為實(shí)現(xiàn)采空區(qū)卸壓瓦斯抽采，高抽巷布設(shè)區(qū)域宜處于穿層裂隙發(fā)育區(qū)。

圖10 頂板上覆巖層不同水平處的垂向位移分布圖

3.3 采動(dòng)覆巖塑性區(qū)分析

圖11為采動(dòng)覆巖塑性區(qū)分布云圖。由圖可知，采空區(qū)上方的巖層主要表現(xiàn)為拉伸破壞，而兩側(cè)煤柱位置呈現(xiàn)出剪切破壞特征，這與采場(chǎng)應(yīng)力分布特征密切相關(guān)。煤層回采后在采空區(qū)正上方形成了卸壓區(qū)，各項(xiàng)應(yīng)力均大幅度減小，使得上覆巖層出現(xiàn)拉伸破壞。兩側(cè)煤柱位置為應(yīng)力集中區(qū)，顯著增大的差應(yīng)力使得該區(qū)域巖層發(fā)生剪切破壞。可根據(jù)巖層的破壞特征對(duì)裂隙發(fā)育區(qū)進(jìn)行預(yù)測(cè)，高抽巷應(yīng)布設(shè)于應(yīng)力集中區(qū)內(nèi)側(cè)的拉伸破壞區(qū)域。

圖11 采動(dòng)覆巖塑性區(qū)分布云圖

4 結(jié) 論

1) 采空區(qū)頂板上覆巖層水平區(qū)域內(nèi)存在應(yīng)力集中區(qū)和卸壓區(qū)，在垂直方向上覆巖層應(yīng)力呈梯形臺(tái)狀分布，高抽巷應(yīng)布置在梯形臺(tái)的“腰線”位置，且高抽巷布置層位越高，在水平位置上越應(yīng)該貼近回采工作面中心區(qū)域。

2) 遠(yuǎn)距離上覆巖層以水平位移為主，近距離上覆巖層以垂向位移為主，導(dǎo)致近距離上覆巖層內(nèi)的穿層裂隙發(fā)育和遠(yuǎn)距離上覆巖層內(nèi)的順層裂隙發(fā)育現(xiàn)象，為實(shí)現(xiàn)采空區(qū)卸壓瓦斯抽采，高抽巷宜布設(shè)于穿層裂隙發(fā)育區(qū)。

3) 采空區(qū)上方的巖層主要表現(xiàn)為拉伸破壞，而兩側(cè)煤柱位置呈現(xiàn)出剪切破壞特征，高抽巷應(yīng)布設(shè)于應(yīng)力集中區(qū)內(nèi)側(cè)的拉伸破壞區(qū)域。