魏風(fēng)清，　程沛棟，　張鈞祥，　杜曉奇

(1.河南理工大學(xué) 安全科學(xué)與工程學(xué)院， 河南 焦作　454003；2.河南理工大學(xué) 能源科學(xué)與工程學(xué)院， 河南 焦作　454003)

?

煤層鉆孔有效抽采半徑數(shù)值模擬研究

魏風(fēng)清1，程沛棟1，張鈞祥2，杜曉奇1

(1.河南理工大學(xué) 安全科學(xué)與工程學(xué)院， 河南 焦作454003；2.河南理工大學(xué) 能源科學(xué)與工程學(xué)院， 河南 焦作454003)

摘要：在煤層瓦斯流動理論和彈性力學(xué)基礎(chǔ)上，考慮煤體滲透率變化對瓦斯流動的影響，建立了含瓦斯煤巖滲流-應(yīng)力耦合數(shù)學(xué)模型，根據(jù)該模型對COMSOL Multiphysics數(shù)值模擬軟件進(jìn)行二次開發(fā)，以某突出礦井為例建立了煤層鉆孔瓦斯抽采模型，并進(jìn)行求解計算，結(jié)果表明煤層瓦斯壓力、煤體滲透率均隨抽采時間的增加而逐漸降低，鉆孔有效抽采半徑隨抽采時間的增加而逐漸擴(kuò)大，但均具有一定的時效性；采用壓力法進(jìn)行現(xiàn)場測試，實測結(jié)果與模擬結(jié)果基本一致，證明了含瓦斯煤巖滲流-應(yīng)力耦合數(shù)學(xué)模型的正確性及采用數(shù)值模擬方法計算煤層鉆孔有效抽采半徑的可行性。

關(guān)鍵詞：煤層鉆孔； 瓦斯抽采； 有效抽采半徑； 含瓦斯煤； 滲流-應(yīng)力耦合； 數(shù)值模擬

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/32.1627.TP.20160705.1458.007.html

0引言

煤層瓦斯抽采是煤礦防治瓦斯事故的根本性措施，在瓦斯災(zāi)害防治方面可起到較好的防突和綜合治理效果[1]。鉆孔的有效抽采半徑是煤層瓦斯抽采的一個重要參數(shù)，直接關(guān)系到鉆孔的設(shè)計、布置和抽采時間。鉆孔煤層透氣性、煤體孔隙壓力等都會對鉆孔有效抽采半徑及抽采效果產(chǎn)生影響，因此，準(zhǔn)確測定鉆孔有效抽采半徑對煤層消突效果具有重要的現(xiàn)實意義[2]。

目前國內(nèi)外學(xué)者不斷完善了煤層瓦斯流動理論[3]研究，并在此基礎(chǔ)上，將基于計算機(jī)數(shù)值模擬的鉆孔有效抽采半徑確定方法應(yīng)用到實際中，如采用瓦斯壓力[4-5]、瓦斯流量衰減[6]、瓦斯含量[7]等指標(biāo)現(xiàn)場測定鉆孔有效抽采半徑；劉軍等[8]通過數(shù)值模擬方法計算鉆孔有效抽采半徑，但未對鉆孔有效抽采半徑的變化規(guī)律進(jìn)行研究。本文以某突出礦井為例，建立考慮煤體孔隙壓力、吸附膨脹應(yīng)力的滲流-應(yīng)力耦合模型，利用COMSOL Multiphysics多物理場軟件模擬建立鉆孔模型，并對其求解計算，模擬研究抽采過程中鉆孔周圍瓦斯流動規(guī)律，確定鉆孔有效抽采半徑并分析其變化規(guī)律，最后采用瓦斯含量指標(biāo)法進(jìn)行現(xiàn)場驗證。

1瓦斯流動模型

1.1瓦斯?jié)B流場方程

由于氣體在多孔介質(zhì)內(nèi)的流動符合氣相質(zhì)量守恒定律，所以在單位時間內(nèi)，控制體流入、流出的質(zhì)量差等于控制體的質(zhì)量變化，則瓦斯流動連續(xù)性方程為

(1)

式中：Q為單位體積含瓦斯煤巖中的瓦斯質(zhì)量，kg/m3；t為時間，s；ρ為煤層瓦斯密度，kg/m3；V為瓦斯?jié)B流速度，m/s。

煤層內(nèi)瓦斯以游離態(tài)和吸附態(tài)2種方式存在。游離態(tài)瓦斯符合理想氣體狀態(tài)方程，吸附態(tài)瓦斯?jié)M足朗格繆爾等溫吸附方程。根據(jù)參考文獻(xiàn)[9]的煤層瓦斯含量方程和參考文獻(xiàn)[10]的有效應(yīng)力方程可得煤層瓦斯流動方程：

(2)

式中：β為煤層瓦斯的壓縮因子，kg/(m3·Pa)；p為瓦斯壓力，Pa；φ為煤體孔隙率；a為單位質(zhì)量煤的極限吸附量，m3/kg；b為煤的吸附常數(shù)，MPa-1；c為單位體積煤中可燃物的質(zhì)量，kg/m3；pn為標(biāo)準(zhǔn)狀況下瓦斯壓力，Pa；k為煤體滲透率，m2；μ為瓦斯的動力黏度，Pa·s；m為Klinkenberg系數(shù)，Pa。

1.2煤體滲透率模型

滲透率是煤層內(nèi)瓦斯流動難易程度的重要指標(biāo)，也是瓦斯抽采的一項重要參數(shù)。含瓦斯煤體是一種多孔介質(zhì)，在瓦斯抽采過程中，煤體孔隙壓力降低、煤基質(zhì)解吸效應(yīng)及有效應(yīng)力的改變會影響滲透率的變化[11]?？紤]滲流場內(nèi)瓦斯壓力、瓦斯解吸和應(yīng)力場煤體體積應(yīng)變的影響，可得到參考文獻(xiàn)[12]中的孔隙率方程。根據(jù)Kozeny-Carman方程可得到孔隙率與滲透率的關(guān)系，進(jìn)而得到煤體滲透率方程：

(3)

式中：k0為煤體初始滲透率,m2；εv為煤體體積應(yīng)變；φ0為煤體初始孔隙率；Δp為瓦斯壓力變化，Pa；ks為煤體體積模量，Pa；ρs為煤體視密度，kg/m3；R為普氏氣體常數(shù)，R=8.314 5 J/(kg·K)；T為煤體溫度，K；Vm為氣體摩爾體積，m3/mol;p0為初始瓦斯壓力,Pa。

1.3應(yīng)力場方程

在煤層瓦斯抽采過程中，煤體多屬于壓縮變形。若忽略溫度影響，煤體變形可認(rèn)為是由地應(yīng)力、瓦斯壓力及瓦斯吸附膨脹產(chǎn)生的膨脹應(yīng)力共同作用影響[13]。由含瓦斯煤體的有效應(yīng)力平衡方程、煤體變形在線彈性變形階段遵循的廣義胡克定律、參考文獻(xiàn)[14]中的煤體膨脹應(yīng)力的計算公式、參考文獻(xiàn)[12]中的煤體體積應(yīng)變張量表示形式，可得出應(yīng)力場煤體變形控制方程：

(4)

式中：G為剪切模量,GPa；ui,uj分別為水平方向和垂直方向的位移分量，m；υ為含瓦斯煤巖泊松比；α為Biot系數(shù)；p′為瓦斯壓力分量,Pa；F為體應(yīng)力，N/m3。

式(3)、式(4)共同構(gòu)成了含瓦斯煤巖滲流-應(yīng)力耦合模型。

2計算機(jī)數(shù)值模擬

2.1幾何模型建立

計算機(jī)數(shù)值模擬以理論數(shù)學(xué)模型為基礎(chǔ)，根據(jù)建立的含瓦斯煤巖滲流-應(yīng)力耦合模型，通過對COMSOL Multiphysics軟件進(jìn)行二次開發(fā)，根據(jù)某突出礦井實際情況建立抽采鉆孔模型，如圖1所示。煤層邊界為100 m×80 m，煤層中布置孔徑為94 mm的抽采鉆孔，孔深70 m，抽采負(fù)壓為13 kPa。

圖1　抽采鉆孔模型

煤層內(nèi)初始瓦斯壓力為1.5 MPa，應(yīng)力場初始位移為0。假設(shè)瓦斯僅在煤層中流動，煤層頂?shù)装鍨椴煌笟鈳r層，煤層四周為自由邊界。模型主要計算參數(shù)見表1。

2.2瓦斯壓力模擬結(jié)果分析

圖2為仿真得到的瓦斯壓力分布曲線?？煽闯鲭S著抽采時間的增加，煤層瓦斯壓力不斷降低，且與鉆孔的距離越近，煤層瓦斯壓力下降速率越大；隨著與鉆孔距離的增加，瓦斯壓力逐漸增大，至煤層深部恢復(fù)至初始瓦斯壓力1.5 MPa。這是由于煤層鉆孔施工破壞了原始煤體的應(yīng)力狀態(tài)，鉆孔周圍形成一個卸壓區(qū)域范圍，該范圍內(nèi)煤層透氣性系數(shù)增大，在抽采初期高瓦斯壓力的影響下，抽采效果顯著；隨著時間的推移，鉆孔周圍的應(yīng)力集中帶逐漸向深部轉(zhuǎn)移，煤體中孔隙被壓縮，煤體滲透率減小，抽采難度加大，導(dǎo)致瓦斯抽采量降低。

表1　模型主要計算參數(shù)

圖2　瓦斯壓力分布曲線

2.3煤體滲透率模擬結(jié)果分析

圖3為距鉆孔3 m處一點煤體滲透率變化曲線?？煽闯鲭S著抽采時間的增加，煤體滲透率不斷降低，且降低速率也逐漸減小。這是因為隨著抽采工作的進(jìn)行，煤層內(nèi)瓦斯壓力持續(xù)降低，煤體受到的有效應(yīng)力增大，導(dǎo)致煤體受壓，煤體內(nèi)孔隙逐漸閉合?？芍趯嶋H瓦斯抽采工作中，隨著抽采時間的增加，雖然煤層中瓦斯壓力繼續(xù)降低，但由于煤體滲透率降低,導(dǎo)致瓦斯壓力變化量逐漸減小，鉆孔瓦斯抽采量也趨于穩(wěn)定。

圖3　煤體滲透率變化曲線

2.4有效抽采半徑的確定

在瓦斯抽采工作中需要確定鉆孔間距，就必須確定煤層鉆孔的有效抽采半徑。依據(jù)《煤礦瓦斯抽采達(dá)標(biāo)暫行規(guī)定》，預(yù)期的防突效果達(dá)標(biāo)瓦斯壓力或瓦斯含量按煤層始突深度處的瓦斯壓力或瓦斯含量取值；沒有考察出煤層始突深度處的煤層瓦斯壓力或瓦斯含量時，分別按照0.74 MPa、8 m3/t取值[15]。因此，本文將煤層內(nèi)距鉆孔最遠(yuǎn)處殘余瓦斯壓力小于0.74 MPa的范圍確定為鉆孔有效抽采半徑。鉆孔有效抽采半徑模擬結(jié)果如圖4所示。

(a) 抽采30 d

(b) 抽采60 d

(c) 抽采90 d

(d) 抽采120 d

從圖4可看出，鉆孔有效抽采半徑隨抽采時間的增加而增大，在抽采30，60，90，120 d條件下，鉆孔有效抽采半徑分別為0.52，0.96，1.27，1.41 m。為了得到有效抽采半徑與抽采時間的變化規(guī)律，對模擬結(jié)果進(jìn)行曲線擬合，如圖5所示?？芍谝欢ǔ椴蓵r間內(nèi)，有效抽采半徑與抽采時間呈正比，之后隨著抽采時間的增加，鉆孔有效抽采半徑增幅逐漸降低。結(jié)合圖3可知，由于煤體滲透率隨抽采時間增加而降低，鉆孔抽采量不斷降低，鉆孔抽采瓦斯的效果會逐漸減弱，每個鉆孔的控制范圍也趨于穩(wěn)定。

圖5　鉆孔有效抽采半徑曲線擬合結(jié)果

3現(xiàn)場實測考察

3.1測試方案

測試地點位于豫南某突出礦井回采工作面。該礦僅有二1煤層全區(qū)可采，無保護(hù)層可采。煤層瓦斯含量隨埋深遞變規(guī)律明顯，由于煤體為強(qiáng)破壞—全粉煤破壞類型，本煤層局部地區(qū)有可能形成高壓瓦斯聚集區(qū)域。在回采工作面掘進(jìn)過程中，實測煤層初始瓦斯含量最大值為7.88 m3/t，煤層原始瓦斯壓力最大值為1.5 MPa。

為確保回采工作的安全進(jìn)行，采用本煤層順層鉆孔抽采煤層瓦斯的區(qū)域防突措施。采用鉆孔流量法對鉆孔有效抽采半徑進(jìn)行測試，測試鉆孔布置如圖6所示。共布置3個測試鉆孔，其中2號鉆孔為抽采孔，1號、3號鉆孔為觀測孔，封孔深度為18 m，煤層傾角為16°，鉆孔施工深度為20 m。鉆孔施工完成后，每天記錄數(shù)據(jù)。

圖6　測試鉆孔布置

3.2數(shù)值模擬與實測結(jié)果

采用鉆孔流量法對鉆孔有效抽采半徑進(jìn)行考察，實測結(jié)果見表2。

抽采半徑為94 mm、抽采負(fù)壓為13 kPa時，鉆孔瓦斯流量的模擬值與實測值對比如圖7所示。可看出模擬值與實測值具有很好的一致性，驗證了本文建立的含瓦斯煤巖滲流-應(yīng)力耦合模型的正確性及采用COMSOL Multiphysics軟件模擬計算煤層鉆孔有效抽采半徑的可行性。

表2　1號鉆孔瓦斯流量測試數(shù)據(jù)

圖7　鉆孔瓦斯流量的模擬值與實測值對比

4結(jié)語

(1) 在瓦斯流動理論及彈性力學(xué)的基礎(chǔ)上，考慮煤體滲透率的動態(tài)變化，建立了含瓦斯煤巖滲流-應(yīng) 力耦合模型，根據(jù)該模型對COMSOL Multiphysics計算機(jī)模擬軟件進(jìn)行二次開發(fā)，以某突出礦井為例，建立了煤層鉆孔瓦斯抽采模型。

(2) 數(shù)值模擬結(jié)果表明，煤層瓦斯壓力、煤體滲透率均隨抽采時間的增加而逐漸降低，鉆孔有效抽采半徑隨抽采時間的增加而逐漸擴(kuò)大，但存在一定的時效性。根據(jù)模擬結(jié)果，確定抽采120 d后鉆孔的抽采范圍趨于穩(wěn)定，得出鉆孔有效抽采半徑為1.41 m。

(3) 基于瓦斯流動理論的數(shù)值模擬方法與壓力法實測結(jié)果基本吻合，說明了數(shù)值模擬的準(zhǔn)確性。因此，可通過掌握煤層瓦斯相關(guān)參數(shù)及煤層概況，方便、準(zhǔn)確地確定煤層鉆孔有效抽采半徑。

參考文獻(xiàn)：

[1]李艷增,王耀鋒,高中寧,等.水力割縫(壓裂)綜合增透技術(shù)在丁集煤礦的應(yīng)用[J].煤礦安全,2011,42(9):108-110.

[2]李建兵,陳向軍.瓦斯抽采鉆孔有效影響半徑測定[J].能源技術(shù)與管理,2012(3):10-11.

[3]韓穎,張飛燕,余偉凡,等.鉆孔瓦斯動態(tài)涌出規(guī)律的實驗研究[J].煤炭學(xué)報,2011,36(11):1874-1878.

[4]孫培德.徑向流場瓦斯壓力分布[J].山東科技大學(xué)學(xué)報(自然科學(xué)版),1989(4):33-39.

[5]馬宏宇,王兆豐,楊宏民,等.寺家莊礦下向穿層鉆孔有效抽放半徑研究[J].煤炭工程,2011(5):96-98.

[6]余陶,盧平,孫金華,等.基于鉆孔瓦斯流量和壓力測定有效抽采半徑[J].采礦與安全工程學(xué)報,2012,29(4):596-600.

[7]劉三鈞,馬耕,盧杰,等.基于瓦斯含量的相對壓力測定有效半徑技術(shù)[J].煤炭學(xué)報,2011,36(10):1715-1719.

[8]劉軍,孫東玲,孫海濤,等.含瓦斯煤固氣耦合動力學(xué)模型及其應(yīng)用研究[J].中國礦業(yè),2013(11):126-130.

[9]張力,何學(xué)秋,李侯全.煤層氣滲流方程及數(shù)值模擬[J].天然氣工業(yè),2002,22(1):23-26.

[10]HU Guozhong, WANG Hongtu, TAN Haixing, et al. Gas seepage equation of deep mined coal seams and its application[J]. Journal of China University of Mining and Technology, 2008, 18(4):483-487.

[11]魏建平,秦恒潔,王登科,等.含瓦斯煤滲透率動態(tài)演化模型[J].煤炭學(xué)報,2015,40(7):1555-1561.

[12]郭平,曹樹剛,張遵國,等.含瓦斯煤體固氣耦合數(shù)學(xué)模型及數(shù)值模擬[J].煤炭學(xué)報,2012,37(增刊2):330-335.

[13]李傳亮,孔祥言,徐獻(xiàn)芝,等.多孔介質(zhì)的雙重有效應(yīng)力[J].自然雜志,1999,21(5):288-292.

[14]吳世躍.煤層氣與煤層耦合運(yùn)動理論及其應(yīng)用的研究[D].沈陽:東北大學(xué),2006.

[15]國家安全生產(chǎn)監(jiān)督管理總局,國家煤礦安全監(jiān)察局.防治煤與瓦斯突出規(guī)定[M].北京：煤炭工業(yè)出版社,2009.

Study on numerical simulation of effective drainage radius of drilling hole in coal seam

WEI Fengqing1,CHENG Peidong1,ZHANG Junxiang2,DU Xiaoqi1

(1.School of Safety Science and Engineering, Henan Polytechnic University, Jiaozuo 454003, China;2.School of Energy Science and Engineering, Henan Polytechnic University, Jiaozuo 454003, China)

Abstract:A seepage-stress coupling model of coal containing gas was established considering effect of coal permeability changing on gas flow based on coal seam gas flow theory and elastic mechanic. The second development of COMSOL Multiphysics numerical simulation software was taken according to the model. A model of gas drainage model of drilling hole in coal seam was established taking a gas outburst mine as an example and resolved with numerical solution. The results show that coal seam gas pressure and coal permeability decrease gradurally with time extending, meanwhile effective radius of drilling hole increases gradually, and both have certain timeliness. The field test is taken by use of a pressure testing method with a result in accordance with the simulated one basically, which verifies correctness of the seepage-stress coupling model of coal containing gas and feasibility of calculating effective drainage radius of drilling hole by numerical simulation method.

Key words:drilling hole in coal seam; gas drainage; effective drainage radius; coal containing gas; seepage-stress coupling; numerical simulation

文章編號：1671-251X(2016)07-0025-05

DOI:10.13272/j.issn.1671-251x.2016.07.007

收稿日期：2015-12-30；修回日期：2016-05-12；責(zé)任編輯：李明。

基金項目：長江學(xué)者和創(chuàng)新團(tuán)隊發(fā)展計劃資助項目(IRT1235)。

作者簡介：魏風(fēng)清(1966-)，男，河南新鄉(xiāng)人，教授級高級工程師，博士，主要從事瓦斯防治技術(shù)方面的研究及教學(xué)工作,E-mail:wfq@hpu.edu.cn。

中圖分類號：TD713

文獻(xiàn)標(biāo)志碼：A網(wǎng)絡(luò)出版時間：2016-07-05 14:58

魏風(fēng)清，程沛棟，張鈞祥，等.煤層鉆孔有效抽采半徑數(shù)值模擬研究[J].工礦自動化，2016,42(7)：25-29.