穿層抽采鉆孔抽采效果的數(shù)值分析

鄧鵬江

摘 要：文章對穿層抽采鉆孔條件下瓦斯運(yùn)移規(guī)律進(jìn)行了探討，以某礦工程實例為背景，運(yùn)用數(shù)值模擬軟件COMSOL Multiphysics模擬了鉆孔孔徑分別為94mm、180mm、300mm的瓦斯流場壓力、壓力梯度的分布，結(jié)合模擬結(jié)果分析了大孔徑鉆孔對穿層抽采鉆孔抽采效果，結(jié)果表明：鉆孔半徑越大，瓦斯的抽采效果越好，但綜合大鉆孔作業(yè)的成本與其對瓦斯抽放效果的貢獻(xiàn)，大鉆孔對某些礦井的瓦斯抽放來說意義不大。

關(guān)鍵詞：穿層抽采鉆孔 瓦斯抽采 抽采效果 COMSOL Multiphysics

中圖分類號：TD713 文獻(xiàn)標(biāo)識碼：A 文章編號：1674-098X（2017）07（c）-0099-03

煤層是非勻質(zhì)且各向異性介質(zhì)，是瓦斯賦存和運(yùn)移的場所。瓦斯在煤層中的運(yùn)移規(guī)律屬于滲流力學(xué)知識體系，穿層抽采鉆孔條件下的煤層瓦斯流動規(guī)律屬于徑向不穩(wěn)定流動，遵循質(zhì)量守恒定律、理想氣體狀態(tài)方程與運(yùn)動方程。

目前，煤礦抽采半徑的確定方法主要分為井下實測法與理論計算法，本文通過數(shù)值模擬軟件COMSOL Multiphysics對提出的變系數(shù)非線性滲流方程進(jìn)行模擬計算，模擬了穿層抽采鉆孔模式下，不同鉆孔半徑對抽采效果的影響，為穿層抽采鉆孔模式下抽采鉆孔的合理布置提供進(jìn)一步的理論計算依據(jù)。

1 煤層瓦斯流動理論方程的建立

為便于對穿層抽采鉆孔條件下瓦斯運(yùn)移規(guī)律研究，特做以下假設(shè)：（1）煤層為均勻連續(xù)性介質(zhì)；（2）煤層頂?shù)装鍑鷰r不透氣且不含瓦斯；（3）煤層的透氣性和孔隙率不受瓦斯壓力變化影響；（4）瓦斯是理想氣體且在煤層中的流動屬于層流滲透，服從達(dá)西定律；（5）瓦斯在煤層中的流動按等溫過程處理；（6）煤層中的瓦斯含量可由煤層瓦斯含量方程近似表示，如式（1）所示：

穿層抽采鉆孔條件下瓦斯的運(yùn)移屬于徑向不穩(wěn)定流動，描述該流動的偏微分方程為：

式（2）中：t為煤層瓦斯向鉆孔流動的時間，d；r為煤體內(nèi)的點距鉆孔中心的距離，m；λ為煤層透氣性系數(shù)，m2/（MPa2·d）；α為煤層瓦斯含量系數(shù)，m3/（m3·MPa1/2）；p0為煤層原始瓦斯壓力，MPa；p1為鉆孔中的瓦斯壓力，MPa。

2 COMSOL Multiphysics的數(shù)值計算及分析

2.1 應(yīng)用實例

某礦實測數(shù)據(jù)：煤厚m；煤層透氣性系數(shù)m2/（MPa2·d）；煤層原始瓦斯壓力p0=2.6MPa；吸附常數(shù)a=32.27；吸附常數(shù)b=0.89；灰分A=18.34%；水分M=0.59%；孔隙率=0.9；煤層視密度m3/t；經(jīng)擬合后求得的瓦斯含量系數(shù)15.48m3/（m3·MPa1/2）。

2.2 建立模型

COMSOL Multiphysics是一款基于有限元分析的數(shù)值仿真軟件，能夠通過求解偏微分方程或偏微分方程組來實現(xiàn)對單場或多場的物理現(xiàn)象進(jìn)行仿真模擬。

穿層鉆孔條件下，選取與鉆孔平行的矩形平面為研究對象，通過COMSOL Multiphysics中二維軸對稱模塊對該矩形進(jìn)行旋轉(zhuǎn)得到模擬所需的幾何模型，且為更逼近實際情況與計算收斂性更好，定義模型右邊界極距為dGeomCharm、物力寬度為1E3*dGeomCharm的無限元域，即該無限元域代表了長為1000m的煤層，如圖1所示。

式（2）是二階變系數(shù)偏微分方程，調(diào)用COMSOL Multi

physics系數(shù)型偏微分方程模塊，選擇在瞬態(tài)條件下求解的方式，模型的假設(shè)方程如式（3）所示：

由式（2）知，式（3）中：擴(kuò)散系數(shù)c=1且各向同性；阻尼或質(zhì)量系數(shù)；對流系數(shù)β在r方向為-1/r，z方向為0；其余各項系數(shù)均為0。

設(shè)置如圖（1）所示瓦斯流場的初始瓦斯壓力平方P=p，在該模型左側(cè)邊界處新增狄氏邊界條件P=r=P12；在該模型右側(cè)邊界處新增狄氏邊界條件P=r=P02。

網(wǎng)格劃分合理與否與方程解算結(jié)果收斂性直接相關(guān)，為提高該模型的解算精度，除選擇特別細(xì)化的自由剖分三角形網(wǎng)格之外，在網(wǎng)格模塊下，調(diào)整模型左側(cè)邊界的邊界層屬性，設(shè)置邊界層數(shù)為8，邊界層拉伸因子為1.2，模型的網(wǎng)格剖分圖如圖（2）所示。

在瞬態(tài)求解器中采用向后拆分公式并將求解器采用的步長改為精確，時間步長設(shè)為0.1d；根據(jù)實際需要調(diào)整計算的總天數(shù)與鉆孔半徑。

2.3 結(jié)果分析

圖（3）為鉆孔半徑R=94mm、抽放時間t=100d后煤層瓦斯壓力分布圖。

由圖（3）可知，穿層抽采鉆孔條件下，鉆孔周圍瓦斯壓力呈同心圓狀分布，靠近鉆孔周圍的瓦斯壓力變化幅度大且距鉆孔距離越近瓦斯壓力越低。由達(dá)西定律可知，煤層中的瓦斯流速與瓦斯壓力梯度成正比，在煤層滲透率一定的條件下，瓦斯流速將隨距鉆孔距離的增加而減小，當(dāng)鉆孔瓦斯壓力梯度為零時（瓦斯壓力不隨距鉆孔距離發(fā)生改變）瓦斯流速亦為零。

分別以R=94mm、R=180mm、R=300mm的鉆孔半徑為模擬條件，繪制t=60d的鉆孔周圍瓦斯壓力梯度分布圖，如圖4所示。

由圖4可知，鉆孔周圍瓦斯壓力梯度在煤壁處取得最大值，且隨鉆孔半徑不斷增大而減小，這與圖3表現(xiàn)一致；對不同半徑的鉆孔而言，瓦斯壓力梯度的最大值隨鉆孔半徑的不斷增大而減小，而抽采影響半徑隨鉆孔半徑的不斷增大而擴(kuò)大。

3 結(jié)語

文章對穿層抽采鉆孔條件下的煤層瓦斯徑向不穩(wěn)定流動規(guī)律進(jìn)行了探討，通過借用數(shù)值模擬軟件COMSOL Multiphysics，以某礦穿層鉆孔實例為研究背景，模擬了不同鉆孔半徑條件下，煤層瓦斯流場的壓力與壓力梯度分布。

由模擬結(jié)果可知：鉆孔周圍瓦斯壓力呈同心圓狀分布且隨距鉆孔距離的增加而增大，而瓦斯壓力梯度隨距鉆孔距離的增加而減??；不同鉆孔半徑條件下，抽采影響半徑隨鉆孔半徑的增加而擴(kuò)大，但擴(kuò)大并不明顯。理論上為擴(kuò)大瓦斯的抽采效果應(yīng)選取鉆孔半徑較大的鉆孔，但綜合大鉆孔作業(yè)的成本與其對瓦斯抽放效果的貢獻(xiàn)，大鉆孔對該礦瓦斯抽放來說意義不大。

參考文獻(xiàn)

[1] 曹新奇，辛海會，徐立華，等.瓦斯抽放鉆孔有效抽放半徑的測定[J].煤炭工程，2009（9）：88-90.

[2] 姚偉，金龍哲，張君.采空區(qū)高位鉆孔瓦斯抽放的數(shù)值模擬[J].北京科技大學(xué)學(xué)報，2010（12）：1521-1525.

[3] 王兆豐，周少華，李志強(qiáng).瓦斯抽采鉆孔有效抽采半徑的數(shù)值計算方法[J].煤炭工程，2011（6）：82-84.

[4] 侯振海，趙耀江，韓升，等.基于COMSOL Multiphysics的鉆孔周圍瓦斯流動規(guī)律數(shù)值模擬[J].煤礦安全，2016 （2）：14-17.

[5] 周福寶，王鑫鑫，夏同強(qiáng).瓦斯安全抽采及其建模[J].煤炭學(xué)報，2014（8）：1659-1666.endprint