趙國貞 馬占國 龔 鵬 王 坤 蔣眾喜

（1.中國礦業(yè)大學(xué)深部巖土力學(xué)與地下工程國家重點(diǎn)實(shí)驗(yàn)室，江蘇省徐州市，221116；2.中國礦業(yè)大學(xué)力學(xué)與建筑工程學(xué)院，江蘇省徐州市，221116）

表土層厚度對地面瓦斯鉆孔穩(wěn)定性影響研究＊

趙國貞1，2馬占國1，2龔 鵬1，2王 坤2蔣眾喜2

（1.中國礦業(yè)大學(xué)深部巖土力學(xué)與地下工程國家重點(diǎn)實(shí)驗(yàn)室，江蘇省徐州市，221116；2.中國礦業(yè)大學(xué)力學(xué)與建筑工程學(xué)院，江蘇省徐州市，221116）

針對淮南礦區(qū)煤礦11-2煤層覆巖物理力學(xué)特性，利用ANSYS軟件中死活單元技術(shù)與材料非線性分析方法，模擬在不同表土層厚度下覆巖的水平變形情況，研究表土層厚度對瓦斯抽采鉆孔穩(wěn)定性的影響。研究結(jié)果表明，不同表土層厚度下，覆巖水平變形 “O”形圈現(xiàn)象明顯，并表現(xiàn)出明顯的對稱性，Ⅰ測線和Ⅲ測線變形錯動明顯，呈現(xiàn)趨勢對稱，梯度差值不對稱；Ⅱ測線表土層變形平緩；隨表土層厚度增加，覆巖水平位移最大變形位置明顯由地表、表土層底部向表土層下部基巖轉(zhuǎn)移，基巖梯度變化區(qū)間個數(shù)明顯減少，梯度差值增大；工作面中部Ⅱ測線位置不容易剪斷，但中部的破碎巖體易被壓實(shí)，后期抽氣效率不高。

瓦斯抽采 表土層厚度 覆巖 水平位移 數(shù)值模擬 死活單元

煤和瓦斯共采問題是采礦工程中普遍關(guān)注的熱點(diǎn)問題，如何在保障煤礦安全開采的前提下，對瓦斯進(jìn)行有效穩(wěn)定的抽采成為研究的重點(diǎn)。在地面打設(shè)瓦斯抽放井是瓦斯抽采的一種行之有效的途徑，與此同時，煤礦采動影響下，地面瓦斯鉆孔的穩(wěn)定性直接影響瓦斯抽放的安全性和可靠性，從而影響到礦井的正常生產(chǎn)。煤層群開采條件下，當(dāng)保護(hù)層開采后，被保護(hù)層得到卸壓，圍巖應(yīng)力狀態(tài)、圍巖結(jié)構(gòu)和瓦斯動力參數(shù)發(fā)生顯著變化，而影響各參數(shù)變化和瓦斯抽采鉆孔穩(wěn)定性的主要因素之一便是表土層厚度。表土層厚度的改變，直接影響著圍巖應(yīng)力場和位移場，因此，伴隨采深逐漸增加，諸多礦區(qū)針對自身地質(zhì)特征展開了關(guān)于影響瓦斯抽采鉆孔穩(wěn)定性的研究。

本文通過ANSYS數(shù)值模擬的方法研究了不同表土層厚度下，覆巖水平變形對瓦斯抽采鉆孔穩(wěn)定性的影響，總結(jié)表土層厚度變化和瓦斯抽采鉆孔變形規(guī)律。研究不僅對鉆孔打設(shè)位置及抽采鉆孔參數(shù)的確定具有指導(dǎo)作用，而且對煤與瓦斯安全共采具有重要意義。

1 模型建立

研究以淮南礦區(qū)煤礦11-2煤層為原型。井田煤系地層被第四系表土層覆蓋，表土層厚230～500 m。工作面走向長2801 m，傾斜長250 m，該階段回風(fēng)巷標(biāo)高-790～-680 m，地面標(biāo)高＋23～＋24 m，煤層厚度1.97～2.99 m，平均2.53 m。結(jié)合覆巖的物理力學(xué)特性，以煤層綜合柱狀圖為依據(jù)，在對網(wǎng)格單元優(yōu)化的基礎(chǔ)上建立力學(xué)分析模型，模型中測線布置如圖1所示。

圖1 模型測線布置圖

圖1中工作面長度為240 m，模擬工作面正上方表土層厚度為300～500 m，并在模擬過程中設(shè)3條測線，分別為距離工作面端頭10 m、100 m、198 m處。利用ANSYS軟件中死活單元技術(shù)、材料非線性分析方法，模擬煤礦開采工作面圍巖的應(yīng)力與變形，分析在不同表土層厚度 （300 m、400 m、500 m）下，覆巖的應(yīng)力分布特征、破壞場特征及巖層移動規(guī)律，研究表土層厚度變化對瓦斯抽采鉆孔穩(wěn)定性的影響。

2 表土層厚度對鉆孔穩(wěn)定性影響

表土層為地表土層，其厚度受地質(zhì)構(gòu)造及地形影響嚴(yán)重，內(nèi)部結(jié)構(gòu)復(fù)雜，厚度參差不齊，隨著煤炭開采深度增加，表土層厚度對煤礦開采的直接影響為工作面載荷負(fù)擔(dān)，同時其多變的地質(zhì)構(gòu)造對鉆孔的穩(wěn)定性也有非常大的影響。

2.1 300 m厚表土層覆巖變形規(guī)律

當(dāng)表土層厚300 m，瓦斯抽放鉆孔深度為550 m時，在采動影響下，根據(jù)工作面不同推進(jìn)距離，利用ANSYS軟件中死活單元技術(shù)研究覆巖水平變形規(guī)律。

（1）當(dāng)采場推進(jìn)至60 m時，工作面采空區(qū)前、后覆巖內(nèi)部都出現(xiàn)很小的水平變形，滯后影響角約為46°。

（2）當(dāng)采場推進(jìn)至100 m時，采空區(qū)前、后方覆巖的水平變形范圍快速增大；受超前壓力影響，工作面前方巖體的水平變形范圍較大，但變形量較小。

（3）當(dāng)采場推進(jìn)至140 m時，采空區(qū)后方覆巖上部表土層變形范圍局部快速增大，在切眼偏后約60 m處正上方的表土層內(nèi)產(chǎn)生了一個變形核區(qū)；超前工作面10 m范圍內(nèi)表土層受超前壓力影響出現(xiàn)不對稱的水平位移。

（4）當(dāng)采場推進(jìn)至180 m時，采空區(qū)前后方上部基巖部分呈現(xiàn)基本對稱的水平變形，邊界移動影響角約為45°，采空區(qū)中偏后方的表土層出現(xiàn)梯度水平位移。

（5）當(dāng)采場推進(jìn)至220 m時，采空區(qū)中后方上部表土層內(nèi)水平位移明顯大于工作面上方的表土層，而采空區(qū)中后方基巖部分的水平位移卻比工作面上方的小，邊界移動影響角約為44°。

表土層厚300 m時，3條測線測得覆巖水平位移如圖2所示。

圖2 表土層厚300 m時覆巖水平位移

由圖2可知，當(dāng)表土層厚300 m時，隨著工作面的推進(jìn)，由于受到剪切力的作用，各測線都產(chǎn)生了一定的水平位移。

2.2 400 m厚表土層覆巖變形規(guī)律

當(dāng)表土層厚400 m，瓦斯抽放鉆孔深度為650 m時，在采動影響下，根據(jù)工作面不同推進(jìn)距離，利用ANSYS軟件中死活單元技術(shù)研究覆巖水平變形規(guī)律。

（1）當(dāng)采場推進(jìn)至60 m、100 m時，覆巖表現(xiàn)出與300 m表土層相似的特征。

（2）當(dāng)采場推進(jìn)至140 m時，切眼偏后方70 m處產(chǎn)生變形核區(qū)，但位置滯后300 m厚表土層10 m。

（3）當(dāng)采場推進(jìn)至180 m時，采空區(qū)前后方上部基巖同樣部分呈現(xiàn)出基本對稱的水平變形，邊界移動影響角稍有增加，約為46°，呈梯度的水平位移仍出現(xiàn)在采空區(qū)中后方。

（4）當(dāng)采場推進(jìn)至220 m時，采空區(qū)中后方60 m處上部表土層的水平位移比工作面上方的大，基巖部分的水平位移呈現(xiàn)出采空區(qū)中后方比工作面上方的小的特征，邊界移動影響角穩(wěn)定在46°。

表土層厚400 m時，3條測線測得覆巖水平位移如圖3所示。

由圖3可知，表土層厚400 m時3條測線整體表現(xiàn)出與表土層厚300 m時不同的變化趨勢。

（1）Ⅰ測線水平位移區(qū)段錯動明顯，同樣表現(xiàn)出明顯的梯度變化，但最大位移出現(xiàn)在表土層下部318 m處，值為37.0 mm，地表處的位移值為27.1 mm，表土層區(qū)間內(nèi)，同樣伴隨表土層深度增加水平位移呈現(xiàn)低降高升的曲線變形；在表土層下部基巖段和圍巖段表現(xiàn)出較表土層厚300 m時更明顯的小范圍區(qū)間性梯度變化，其區(qū)間水平位移峰值最大值分別為 （445 m，36.1 mm）、 （508 m，3.0 mm）、（572 m，19.5 mm），谷值最大值分別為 （447 m，-11.4 mm）、（540 m，-3.8 mm），峰值與谷值最大梯度差值為47.5 mm，峰值與谷值較表土層厚300 m時波動明顯，增大16.4%。

（2）Ⅱ測線水平位移整體變形平穩(wěn)，梯度變化較小，表土層內(nèi)水平位移最大變形出現(xiàn)在地表8.3 mm；基巖內(nèi)部變形平穩(wěn)，僅有小范圍波動，其區(qū)間水平位移峰值最大值分別為 （477 m，1.5 mm）、（541 m，0.7 mm），谷值最大值分別為（445 m，-6.6 mm）、（508 m，-1.8 mm）、（572 m，-8.6 mm），區(qū)間內(nèi)峰值和谷值的最大梯度差值為10.1 mm，較Ⅰ測線減少78.7%。

圖3 表土層400 m時覆巖水平位移

（3）Ⅲ測線表現(xiàn)出與Ⅰ測線相對稱的梯度變化規(guī)律，其最大位移出現(xiàn)在地表下部318 m處，值為37.6 mm，高于Ⅰ測線；在表土層下部基巖段和圍巖段表現(xiàn)出小范圍的區(qū)間性梯度變化，其區(qū)間水平位移峰值最大值分別為 （477 m，10.1 mm）、（541 m，4.0 mm），谷值最大值分別為 （445 m，-36.7 mm）、（509 m，-4.5 mm）、（572 m，-19.1 mm），峰值與谷值最大梯度差值為46.8 mm，整條測線的最大差值為47.6 mm。

（4）對比表土層厚300 m時的3條測線，表土層厚400 m時，Ⅰ測線梯度差值明顯增大，對瓦斯鉆孔的影響更加明顯；Ⅱ測線與Ⅲ測線變化并不明顯，但表土層厚度的增加導(dǎo)致梯度變化區(qū)間減少，小范圍內(nèi)位移變化劇烈，由此可以推斷為表土層載荷增大導(dǎo)致瓦斯鉆孔在小范圍內(nèi)將承受劇烈沖擊，不利于鉆孔穩(wěn)定。

2.3 500 m厚表土層覆巖變形規(guī)律

當(dāng)表土層厚500 m，瓦斯抽放鉆孔深度為750 m時，在采動影響下，根據(jù)工作面不同推進(jìn)距離，利用ANSYS軟件中死活單元技術(shù)研究覆巖水平變形規(guī)律。

（1）當(dāng)采場推進(jìn)至60 m時，滯后影響角約為40°。

（2）當(dāng)采場推進(jìn)至100 m時，切眼偏后約60 m位置便產(chǎn)生了一個變形核區(qū)，早于表土層厚300 m、400 m時約40 m。

（3）當(dāng)采場推進(jìn)至180 m時，采空區(qū)前后方上部基巖部分不再呈現(xiàn)對稱的水平變形，邊界移動影響角繼續(xù)增大為55°，采空區(qū)中后方的厚表土層出現(xiàn)梯度水平位移。

（4）當(dāng)采場推進(jìn)至220 m時，表土層的水平位移表現(xiàn)出采空區(qū)中后方40 m處比工作面上方的大，比400 m時提前了20 m；基巖部分的水平位移呈現(xiàn)采空區(qū)中偏后方比工作面上方的小，邊界移動影響角仍穩(wěn)定為46°。

表土層厚500 m時，3條測線測得覆巖水平位移如圖4所示。

圖4 表土層厚500 m時覆巖水平位移

由圖4可知，表土層水平位移最大變形位置明顯向表土層下部基巖推移，由地表、表土層底部轉(zhuǎn)移到表土層下部基巖，同時，基巖梯度變化區(qū)間個數(shù)明顯減少，梯度差值變化更加劇烈。

（1）Ⅰ測線水平位移區(qū)段仍舊表現(xiàn)出明顯的梯度變化，但最大位移出現(xiàn)在表土層下部662 m處，位移值為53.4 mm，地表處的位移值為14.8 mm，相比表土層厚300 m時，表土層底部覆巖水平位移明顯增大，地表變形明顯減小，表土層區(qū)間內(nèi)，同樣伴隨表土層深度增加水平位移呈現(xiàn)低降高升的曲線變形；在表土層下部基巖段和圍巖段表現(xiàn)出較表土層厚400 m時更明顯的小范圍區(qū)間性梯度變化，梯度區(qū)間個數(shù)明顯減少，其區(qū)間水平位移峰值最大值為 （662 m，53.4 mm），谷值最大值為（589 m，-30.3 mm），峰值與谷值最大梯度差值為83.7 mm，峰值與谷值較表土層厚300 m時波動明顯，增大105.1%。

（2）Ⅱ測線水平位移整體變形平穩(wěn)，梯度變化較小，表土層內(nèi)水平位移最大變形出現(xiàn)在表土層底部478.4 m處，值為7.5 mm；基巖內(nèi)部變形平穩(wěn)，僅有小范圍波動，其區(qū)間水平位移峰值最大值為 （589 m，3.3 mm），谷值最大值為 （662 m，-26.2 mm），區(qū)間內(nèi)峰值和谷值的最大梯度差值為29.5 mm，較Ⅰ測線減少64.8%。

（3）Ⅲ測線表現(xiàn)出與Ⅰ測線相對稱的梯度變化規(guī)律，其最大位移出現(xiàn)在表土層下部基巖662 m處，值為53.0 mm；在表土層下部基巖段和圍巖段表現(xiàn)出小范圍的區(qū)間性梯度變化，其區(qū)間水平位移峰值最大值為 （589 m，29.6 mm），谷值最大值為 （662 m，-53.0 mm），峰值與谷值最大梯度差值為82.6 mm，與Ⅰ測線相近。

（4）對比表土層厚300 m、400 m時的3條測線，表土層厚500 m時，受載荷影響，Ⅰ測線梯度差值明顯增大，區(qū)間個數(shù)明顯減小，最大水平位移明顯向下部基巖轉(zhuǎn)移，對瓦斯鉆孔的影響也更加明顯。

3 結(jié)論

通過模擬不同表土層厚度下覆巖的水平變形情況，得到以下結(jié)論：

（1）在不同的表土層厚度下，覆巖水平變形“O”形圈現(xiàn)象明顯，并表現(xiàn)出明顯的對稱性。Ⅰ測線和Ⅲ測線變形錯動明顯，關(guān)于工作面中軸線呈現(xiàn)趨勢對稱，梯度差值不對稱；Ⅱ測線表土層變形平緩，基巖層變化劇烈。

（2）受表土層厚度影響，表土層水平位移最大變形位置明顯由地表、表土層底部向表土層下部基巖轉(zhuǎn)移，同時，隨表土層厚度增加，基巖梯度變化區(qū)間個數(shù)明顯減少，梯度差值變化更加劇烈。

（3）工作面中部Ⅱ測線位置不容易剪斷，但中部的破碎巖體易被壓實(shí)，后期抽氣效率不高；Ⅰ測線和Ⅲ測線水平位移梯度變化明顯，覆巖劇烈變動，瓦斯抽放鉆孔易被剪斷或拉斷，Ⅰ測線位于“O”形圈的外邊緣，相對穩(wěn)定、變形梯度較Ⅲ測線??；從長遠(yuǎn)來看，鉆孔位置適宜打在工作面邊緣Ⅰ測線位置。

［1］汪理全等.煤層 （群）上行開采技術(shù) ［M］.北京：煤炭工業(yè)出版社.1995

［2］周世寧，鮮學(xué)福，朱旺喜.煤礦瓦斯災(zāi)害防治理論戰(zhàn)略研討 ［M］.徐州：中國礦業(yè)大學(xué)出版社，2001

［3］馬占國，涂敏，馬繼剛等.遠(yuǎn)距離下保護(hù)層開采煤巖體變形特征 ［J］.采礦與安全工程學(xué)報，2008（3）

［4］袁亮.瓦斯治理理念和煤與瓦斯共采技術(shù) ［J］.中國煤炭，2010（6）

［5］謝文強(qiáng)，李煥.煤層鉆孔瓦斯抽放參數(shù)模擬與應(yīng)用［J］.中國煤炭，2012（4）

［6］李霄尖，姚精明，劉會田等.基于UDEC的高位鉆孔抽放瓦斯數(shù)值模擬研究 ［J］.中國煤炭，2008（8）

［7］羅志中.淮南礦區(qū)高瓦斯煤層群深井開拓實(shí)踐［J］.中國煤炭，2010（1）

［8］葛春貴.開采下保護(hù)層時瓦斯抽采方案的優(yōu)化［J］.煤炭科技，2009（2）

［9］程遠(yuǎn)平，周德永，俞啟香等.保護(hù)層卸壓瓦斯抽采及涌出規(guī)律研究 ［J］.采礦與安全工程學(xué)報，2006（1）

［10］涂敏，黃乃斌，劉寶安等.遠(yuǎn)距離下保護(hù)層開采上覆煤巖體卸壓效應(yīng)研究 ［J］.采礦與安全工程學(xué)報，2007（4）

On impact of thickness of topsoil on stability of ground gas drilling

Zhao Guozhen1，2，Ma Zhanguo1，2，Gong Peng1，2，Wang Kun2，Jiang Zhongxi2
（1.State Key Laboratory for Geomechanics ＆ Deep Underground Engineering，China University of Mining ＆Technology，Xuzhou，Jiangsu 221116，China；2.School of Mechanics ＆Civil Engineering，China University of Mining＆Technology，Xuzhou，Jiangsu 221116，China）

Based on the physical and mechanical properties of overlying strata of 11-2 coal seam located in Huainan mining area，the paper researches the impact of the thickness of topsoil on stability of gas drainage boreholes.It uses the birth-death element technology of ANSYS software and the material nonlinear analysis methods to simulate the horizontal deformation of overlying strata with different thickness of topsoil.The results show that O-ring phenomenon of overlying strata is obvious and the horizontal deformation presents obvious symmetry with different thickness of topsoil；ⅠandⅢsurvey lines show the tendency symmetry and gradient difference asymmetry with evident dislocation deformation；however，deformation of topsoil ofⅡsurvey line is gentle；with the increases of topsoil thickness，the maximum horizontal displacement deformation position of overlying strata transfers to the bedrock from the surface and the bottom of topsoil，and the number of gradient interval significantly decreases and gradient difference increases；the position ofⅡsurvey line in middle part of working face is not easy to cut，but the broken rock is easily compacted，so the gas post-extraction is not efficient.

gas drainage，thickness of topsoil，overlying strata，horizontal displacement，numerical simulation，birth-death element

TD325 TD712

A

國家自然科學(xué)基金項(xiàng)目 （51074163），國家自然科學(xué)基金重點(diǎn)項(xiàng)目 （50834005），江蘇省自然科學(xué)基金 （BK2009092），中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助（2012LWB66），江蘇省研究生培養(yǎng)創(chuàng)新工程 （CXZZ12＿0923）

趙國貞 （1985-），男，山東泰安人，博士，從事采動巖體力學(xué)與圍巖控制理論方面的科研工作。

（責(zé)任編輯 張毅玲）