張嘉勇， 陳 建，周鳳增，張新雙，張金海，李鳳志

(1.華北理工大學(xué) 礦業(yè)工程學(xué)院，河北 唐山 063009；2.開(kāi)灤(集團(tuán))有限責(zé)任公司，河北 唐山 063000；3.開(kāi)灤(集團(tuán))有限責(zé)任公司 錢(qián)家營(yíng)礦業(yè)分公司，河北 唐山 063000)

基于數(shù)值模擬的煤層瓦斯抽采半徑合理化試驗(yàn)研究

張嘉勇1， 陳 建1，周鳳增2，張新雙2，張金海3，李鳳志3

(1.華北理工大學(xué) 礦業(yè)工程學(xué)院，河北 唐山 063009；2.開(kāi)灤(集團(tuán))有限責(zé)任公司，河北 唐山 063000；3.開(kāi)灤(集團(tuán))有限責(zé)任公司 錢(qián)家營(yíng)礦業(yè)分公司，河北 唐山 063000)

針對(duì)井下煤與瓦斯的突出問(wèn)題，在開(kāi)灤錢(qián)家營(yíng)礦進(jìn)行了煤層瓦斯抽放半徑測(cè)定的實(shí)驗(yàn)。運(yùn)用壓力指標(biāo)法對(duì)5煤層的瓦斯抽放半徑進(jìn)行測(cè)定，結(jié)果表明5煤層的抽放半徑大于1m；采用改進(jìn)含量指標(biāo)法現(xiàn)場(chǎng)測(cè)定了5煤層瓦斯的有效抽放半徑，在3～4m范圍內(nèi)；采用RFPA2D-GasFlow模擬軟件對(duì)瓦斯在抽放過(guò)程中的滲流情況進(jìn)行模擬，最終確定模擬有效半徑為3.4m；通過(guò)數(shù)學(xué)分析最終確定抽放鉆孔布置最佳間距為4.8m，測(cè)定結(jié)果為錢(qián)家營(yíng)礦防突措施的開(kāi)展提供關(guān)鍵數(shù)據(jù)和依據(jù)。

含量指標(biāo)法；礦井防突；抽放半徑；數(shù)值分析

隨著礦井開(kāi)采年限增加，可采煤層的賦存深度隨之加深[1]。煤層深部瓦斯壓力高[2]、地應(yīng)力增大，煤層發(fā)生突出可能性大，礦井的安全生產(chǎn)受到威脅，急需采取安全有效的方法治理瓦斯、預(yù)防突出。煤層瓦斯抽放是治理礦井瓦斯的重要措施之一[3]。煤層瓦斯抽放措施的關(guān)鍵在于確定煤層瓦斯抽放半徑[4]。錢(qián)家營(yíng)礦-850m水平5煤層瓦斯涌出量較大，煤的硬度小，瓦斯壓力大，煤層滲水且水量不穩(wěn)定[5]，煤層突出危險(xiǎn)性較高。因此，在5煤層采用抽放瓦斯的方法降低煤層瓦斯壓力。瓦斯抽放過(guò)程中，確定煤層瓦斯抽放半徑是關(guān)鍵，煤層瓦斯抽放的效果取決于瓦斯抽采半徑[6]。因此，展開(kāi)了對(duì)錢(qián)家營(yíng)礦5煤層瓦斯抽放半徑的試驗(yàn)研究。

1 測(cè)定方法選擇及現(xiàn)場(chǎng)測(cè)定

煤層瓦斯抽放半徑的測(cè)定方法一般分為3種：壓力指標(biāo)法、含量指標(biāo)法、相對(duì)壓力指標(biāo)法[7]。首先選擇方便快捷的壓力指標(biāo)法對(duì)錢(qián)家營(yíng)礦5煤層瓦斯抽放半徑進(jìn)行測(cè)定。

1.1 布置鉆孔

在工作面布置5個(gè)測(cè)壓鉆孔，分別標(biāo)號(hào)1，2，3，4，5號(hào)。5個(gè)鉆孔承立體式布置，減小測(cè)壓鉆孔對(duì)煤層瓦斯抽放的影響。3號(hào)孔為測(cè)試孔，2，4，1，5號(hào)孔分別距3號(hào)孔1，2，3，4m，并分別布置在3號(hào)鉆孔的下、右、左、上方向。具體分布如下圖1所示。

圖1 鉆孔分布示意

1，2，3，4，5號(hào)鉆孔現(xiàn)場(chǎng)布置數(shù)據(jù)如表1所示。

表1 5煤層瓦斯參數(shù)測(cè)定鉆孔布置參數(shù)

具體施工過(guò)程：在-850m水平副石門(mén)處，按照鉆孔布置圖所示從巷道左幫向5煤層打測(cè)壓鉆孔，分別為1，2，3，4，5號(hào)鉆孔。成孔后迅速清洗鉆孔，安裝鉆孔套管進(jìn)行封孔。具體操作步驟如下所述。

(1)第1鉆孔成孔后，放入8～10m的直徑為108mm的套管，進(jìn)行封孔。封孔過(guò)程中，要保證有3m以上的距離為聚氨酯乳液和煤孔壁直接接觸，并在管口處用速凝水泥封孔。

(2)鉆孔口處，用法蘭盤(pán)與套管對(duì)接，法蘭盤(pán)上連接有壓力表和抽放閥門(mén)。

(3)在距離第1鉆孔1m處打第2鉆孔，其他鉆孔施工方式同第1鉆孔施工過(guò)程類(lèi)似。

(4)將所有鉆孔封閉，進(jìn)行氣密性試驗(yàn)，確保鉆孔密封效果在誤差范圍內(nèi)。

1.2 壓力指標(biāo)法現(xiàn)場(chǎng)試驗(yàn)

分別向1，2，4，5號(hào)鉆孔內(nèi)注入2MPa的壓縮空氣，關(guān)閉閥門(mén)，并在3號(hào)鉆孔進(jìn)行抽放瓦斯，若某鉆孔內(nèi)氣體壓力逐步減小，可認(rèn)定其鉆孔位置在3號(hào)鉆孔的抽放半徑內(nèi)。

現(xiàn)場(chǎng)記錄各個(gè)壓力表數(shù)據(jù)，如表2所示。

表2 5煤層瓦斯測(cè)壓鉆孔壓力數(shù)據(jù)

由表2數(shù)據(jù)分析可知，3號(hào)鉆孔進(jìn)行抽放工作時(shí)，1，2，4號(hào)鉆孔內(nèi)氣體壓力小幅度降低，后逐漸增高。氣體充入鉆孔后，會(huì)逐漸滲透到孔壁周?chē)牧严秲?nèi)，導(dǎo)致鉆孔內(nèi)氣體壓力降低。氣體壓力逐步增高是由煤層滲水引起的，煤層水緩慢滲入鉆孔中導(dǎo)致鉆孔內(nèi)部氣體被壓縮，氣體壓力升高。由于抽放作用于裂隙之間，滲流速度緩慢，無(wú)法排除滲水對(duì)實(shí)驗(yàn)的干擾。因此，需要采用新方法對(duì)3號(hào)鉆孔抽放半徑進(jìn)行測(cè)定。

1.3 含量指標(biāo)法測(cè)定有效抽放半徑

由于壓力指標(biāo)法所測(cè)數(shù)據(jù)無(wú)法排除誤差原因，因此，采用含量指標(biāo)法測(cè)定有效抽放半徑。

傳統(tǒng)含量指標(biāo)法測(cè)定有效抽放半徑，主要通過(guò)測(cè)定游離瓦斯和吸附瓦斯之和，但整個(gè)實(shí)驗(yàn)過(guò)程較復(fù)雜緩慢，所以對(duì)此方法進(jìn)行改進(jìn)，主要改進(jìn)過(guò)程為：實(shí)驗(yàn)前抽測(cè)到鉆孔內(nèi)氣體成分，然后進(jìn)行封孔，依次從5，1，4，2號(hào)鉆孔進(jìn)行注氮，對(duì)3號(hào)鉆孔抽放瓦斯，監(jiān)測(cè)該鉆孔內(nèi)氣體成分變化。一段時(shí)間后，確定首先使3號(hào)鉆孔內(nèi)氮?dú)庠隽砍^(guò)誤差范圍的注氮鉆孔，此鉆孔距3號(hào)鉆孔的距離為有效抽放半徑。此方法為氮?dú)鈽?biāo)定法。

1.4 氮?dú)鈽?biāo)定法現(xiàn)場(chǎng)試驗(yàn)

抽放前工作面氣體濃度色譜分析結(jié)果：氮?dú)?8.53%，氧氣20.78%，甲烷0.41%，其他氣體0.28%。

1.4.1 5號(hào)鉆孔注氮

5號(hào)鉆孔持續(xù)注入氮?dú)?，不關(guān)閉注氮閥門(mén)，保證整個(gè)過(guò)程中鉆孔內(nèi)氣體壓力為0.04MPa。對(duì)3號(hào)鉆孔進(jìn)行抽放，監(jiān)測(cè)3號(hào)鉆孔內(nèi)氣體濃度變化。監(jiān)測(cè)數(shù)據(jù)如表3所示。

表3 5號(hào)鉆孔注氮后3號(hào)鉆孔氣體濃度

由表3可以分析出，3號(hào)鉆孔內(nèi)的甲烷濃度先升高后降低，并降到抽放前工作面測(cè)定濃度；氮?dú)鉂舛认冉档秃笊撸⒒謴?fù)到抽放前工作面測(cè)定濃度。以上分析表明，由于3號(hào)鉆孔的抽放作用，煤層賦存瓦斯涌入鉆孔，導(dǎo)致甲烷氣體濃度增加，氮?dú)鉂舛冉档?。隨著抽放過(guò)程的進(jìn)行，15d后瓦斯?jié)舛然痉€(wěn)定，并與抽放前工作面測(cè)定濃度一致，此時(shí)氮?dú)鉂舛任蠢^續(xù)升高，因此可以判定5號(hào)鉆孔距3號(hào)鉆孔的距離大于5煤層的抽放半徑。

1.4.2 鉆孔注氮

1號(hào)鉆孔持續(xù)注入氮?dú)?，不關(guān)閉注氮閥門(mén)，保證整個(gè)過(guò)程中鉆孔內(nèi)氣體壓力為0.04MPa。對(duì)3號(hào)鉆孔進(jìn)行抽放，監(jiān)測(cè)3號(hào)鉆孔內(nèi)氣體濃度變化。監(jiān)測(cè)數(shù)據(jù)如表4所示。

表4 1號(hào)鉆孔注氮后3號(hào)鉆孔氣體濃度

分析表4可知，由于向5號(hào)鉆孔持續(xù)注入氮?dú)鈱?shí)驗(yàn)中3號(hào)鉆孔的瓦斯抽放工作已經(jīng)完成，所以數(shù)據(jù)中瓦斯?jié)舛茸兓^小，氮?dú)鉂舛戎鸩缴?，并趨于穩(wěn)定。分析表明，1號(hào)鉆孔內(nèi)的氮?dú)庠诔榉抛饔孟逻M(jìn)入3號(hào)鉆孔，1號(hào)鉆孔到3號(hào)鉆孔的距離小或者等于5煤層的抽放半徑。所以5煤層的瓦斯抽放半徑范圍為3～4m。

2 數(shù)值模擬

2.1 抽放參數(shù)確定

根據(jù)煤層地質(zhì)條件和方案設(shè)計(jì)要求等因素，將錢(qián)家營(yíng)礦5煤瓦斯抽放數(shù)值模擬主要參數(shù)設(shè)定為：煤層瓦斯壓力分別為0.89MPa (模擬時(shí)確定為0.9，0.6MPa)；煤層埋藏深度約為850m，其產(chǎn)生的地應(yīng)力約為19MPa(模擬選取20MPa)；根據(jù)礦井的瓦斯抽放實(shí)際情況，一般選擇抽放鉆孔直徑約為94～108mm(模擬時(shí)選取100mm)；鉆孔抽放負(fù)壓選取為0.04MPa；巖層透氣性非常小，模擬時(shí)將參數(shù)設(shè)置為零[8]；煤層瓦斯抽放效果主要受煤層的透氣性影響，由于井下煤層賦存環(huán)境的差異性，導(dǎo)致了煤層透氣性變化較大。因此，對(duì)煤層抽放半徑的研究可轉(zhuǎn)化為對(duì)煤層透氣性的研究。

2.2 數(shù)值模擬軟件及模型

RFPA2D-GasFlow是以彈性力學(xué)、彈性損傷理論為基礎(chǔ)理論，耦合的模擬巖石破壞過(guò)程的2維系統(tǒng)。

2.2.1 數(shù)值模型建立

根據(jù)現(xiàn)場(chǎng)工作面實(shí)際情況，建立數(shù)值模型，如圖2所示。

圖2 有效抽放半徑數(shù)值模型

圖2(a)模擬模塊選擇尺寸為11m×11m，將整個(gè)部分按110×110的形式劃分單元，最外側(cè)1m為巖層，中間10m×10m為煤層，其中黑體部分為取芯鉆孔(2m×0.1m)；在瓦斯抽放鉆孔的兩側(cè)設(shè)置的邊界條件是應(yīng)力邊界條件，其值固定為20MPa；兩外側(cè)限制位移。

數(shù)值模型圖2(b)尺寸為4m×4m，劃分為200×200個(gè)單元。模型四周為固定壓力邊界20MPa，且無(wú)瓦斯流動(dòng)補(bǔ)給，鉆孔就抽放此區(qū)域內(nèi)的瓦斯，鉆孔直徑100mm，無(wú)抽放負(fù)壓自然排放。

2.2.2 數(shù)值模擬結(jié)果

確定抽放煤層相關(guān)參數(shù)、抽放時(shí)間和抽放負(fù)壓。根據(jù)錢(qián)家營(yíng)礦實(shí)際情況，進(jìn)行數(shù)值模擬，確定抽放時(shí)間為30d(時(shí)間跨度為1d，模擬抽放時(shí)間為30d)，抽放負(fù)壓確定為0.04MPa，排放負(fù)壓為0。

(1)抽放鉆孔模擬結(jié)果 5煤層的透氣性為0.0669m2/(MPa2·d)，其模擬此透氣性系數(shù)時(shí)的鉆孔瓦斯流動(dòng)圖和鉆孔瓦斯流量圖如圖3(a)、(b)所示。

圖3 5煤鉆孔瓦斯抽放數(shù)值模擬

由鉆孔瓦斯流動(dòng)圖3(a)和鉆孔瓦斯流量圖3(b)可知，當(dāng)抽放孔孔徑為100mm，抽放負(fù)壓為0.04MPa時(shí)，抽放鉆孔的有效影響半徑分別約為3.5，3.4m，由圖3(b)可知臨界流量和鉆孔周?chē)咚沽髁肯嘟挥诰嚯x鉆孔中心3.4m的位置，可得鉆孔有效抽放范圍為3.4m。

(2)抽放半徑模擬結(jié)果分析 RFPA2D-GasFlow數(shù)值模擬生成有2組圖片，分別是抽放鉆孔周?chē)后w的鉆孔瓦斯流動(dòng)圖和流量曲線(xiàn)圖。當(dāng)抽放負(fù)壓和抽放孔徑一定時(shí)，抽放鉆孔有效影響半徑大小和煤層透氣性系數(shù)成正比。煤層透氣性越小，鉆孔瓦斯流動(dòng)阻力越大，抽放鉆孔有效影響半徑就越小，鉆孔流量就越小(鉆孔流量曲線(xiàn)與臨界流量曲線(xiàn)閉合區(qū)域小)。

在鉆孔瓦斯流動(dòng)圖中，鉆孔周?chē)幱皡^(qū)域表示鉆孔瓦斯流動(dòng)，陰影未覆蓋的區(qū)域，因其瓦斯流量較小，軟件生成分辨率有限未顯示。因此，可以認(rèn)為鉆孔周?chē)幱皡^(qū)域即為瓦斯抽放鉆孔有效抽放范圍，按照模型尺寸(10m×10m)對(duì)比鉆孔周?chē)幱皡^(qū)域的范圍，可測(cè)定，上述透氣性條件下的瓦斯抽放鉆孔有效抽放半徑3.5m。

在流量曲線(xiàn)圖中，流量的正負(fù)表示流動(dòng)方向相反，離鉆孔越遠(yuǎn)，其流量越小。當(dāng)流量降低到一定數(shù)值后，該位置到鉆孔的距離為鉆孔抽放瓦斯的有效范圍，此數(shù)值也為臨界流量。由計(jì)算分析可知，5煤鉆孔周?chē)罕趩挝幻娣e瓦斯涌出臨界流量約為8.0×10-4m3/(min·m2)，根據(jù)鉆孔瓦斯流量曲線(xiàn)和臨界流量直線(xiàn)可確定瓦斯抽放鉆孔有效抽放范圍為3.4m。

對(duì)以上2種方式得出不同透氣性條件下瓦斯抽放鉆孔的有效抽放半徑進(jìn)行統(tǒng)計(jì)，根據(jù)2種數(shù)據(jù)結(jié)果對(duì)比可知，鉆孔有效抽放范圍為3.4m。

根據(jù)數(shù)值模擬結(jié)果，結(jié)合錢(qián)家營(yíng)礦5煤層實(shí)際情況，綜合考量后確定瓦斯抽放鉆孔有效抽放半徑約為3.4m。

有效抽放半徑是指鉆孔抽放一定時(shí)間后能消除突出的范圍[9]，這個(gè)范圍用以鉆孔為中心的半徑來(lái)表示。假定煤層是均質(zhì)煤體[10]，抽放鉆孔在一定抽放時(shí)間下的影響范圍就是以鉆孔為圓心，以該時(shí)間下的抽放影響半徑為半徑的一個(gè)圓[11]，見(jiàn)圖4。

圖4 合理抽放鉆孔間距

從圖4可知，A為DE的中點(diǎn)，B，C為抽放鉆孔的中心，抽放鉆孔在一定抽放時(shí)間下的影響半徑為R，AB的長(zhǎng)度為L(zhǎng)，陰影部分為抽放空白帶。由勾股定理可得：

式中，L為AB長(zhǎng)度，m；R為抽放鉆孔有效影響半徑，m。

經(jīng)計(jì)算錢(qián)家營(yíng)礦-850m水平5煤抽放鉆孔合理間距應(yīng)為4.8m。

3 結(jié) 論

(1)通過(guò)改進(jìn)的含量指標(biāo)法確定了錢(qián)家營(yíng)礦5煤層的有效抽放半徑在3～4m之間。

(2)運(yùn)用RFPA2D-GasFlow模擬軟件確定了錢(qián)家營(yíng)礦5煤層的有效抽放半徑為3.4m。

(3)通過(guò)理論分析和計(jì)算，確定錢(qián)家營(yíng)礦5煤層的最佳鉆孔間距為4.8m，為礦井瓦斯抽采鉆孔施工提供基礎(chǔ)數(shù)據(jù)。

[1]袁 亮.低透高瓦斯煤層群安全開(kāi)采關(guān)鍵技術(shù)研究[J].巖石力學(xué)與工程學(xué)報(bào)，2008，27(7):1370-1379.

[2]謝廣祥,胡祖祥,王 磊.深部高瓦斯工作面煤體采動(dòng)擴(kuò)容特性研究[J].煤炭學(xué)報(bào),2014,29(1):91-96.

[3]張學(xué)博,姬忠超,高建良.順層鉆孔有效抽放半徑影響因素的數(shù)值分析[J].安全與環(huán)境學(xué)報(bào),2013,13(2):157-160.

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

[5]關(guān)維娟,張國(guó)樞,趙志根,等.煤與瓦斯突出多指標(biāo)綜合辨識(shí)與實(shí)時(shí)預(yù)警研究[J].采礦與安全工程學(xué)報(bào),2013,30(6):922-929.

[6]余 陶,盧 平,朱貴旺,等.穿層鉆孔水力壓裂強(qiáng)化抽采瓦斯消突技術(shù)應(yīng)用研究[J].安全與環(huán)境學(xué)報(bào),2010,10(6):172-175.

[7]鐘后選.有效壓力指標(biāo)法測(cè)定鉆孔有效抽放半徑[J].煤礦安全,2013,44(8):166-168.

[8]張 超，齊黎明，張志剛.低透氣性煤層鉆孔瓦斯運(yùn)移及抽放參數(shù)研究[J].煤炭技術(shù)，2014,33(9):5-8.

[9]馬永慶.松軟低透大采深突出煤層瓦斯的區(qū)域治理技術(shù)[J].煤炭科學(xué)技術(shù),2010,38(6):57-60.

[10]張嘉勇，郭立穩(wěn)，羅新榮.高壓水射流掏槽防突技術(shù)參數(shù)數(shù)值模擬與試驗(yàn)研究[J].采礦與安全工程學(xué)報(bào)，2013，30(5):785-790.

[11]劉 軍，孟兆豐，李學(xué)臣，等.消除礦井瓦斯抽采空白帶方法的研究[J].煤炭科學(xué)技術(shù)，2012，40(12)：59-61，87.

[責(zé)任編輯：李 青]

Experimental Study of Coal Seam Gas Drainage Radius Based on Numerical Simulation

ZHANG Jia-yong1，CHEN Jian1，ZHOU Feng-zeng2,ZHANG Xin-shuang2，ZHANG Jin-hai3,LI Feng-zhi3

(1.School of Mine Engineering，North China University of science and technology，Tangshan 063009，China；2.Kailuan(Corporation) Co.，Ltd.，Tangshan 063000，China；3.Qianjiaying Mine Branch Company，Kailuan(Corporation) Co.，Ltd.，Tangshan 063000，China)

To the problem of rockburst in deep，and experimental detection of coal seam gas drainage radius experimental was done in Qianjiaying mine of Kailuan，and gas drainage radius of 5# coal seam was detected with pressure index method，the results showed that the gas drainage radius of 5# coal seam more than 1m，and then an effective radius of 5# coal seam was field test with improved content index method，it is between 3～4m，then seepage situation during drainage process was simulated with software RFPA2D-GasFlow，an effective radius was determined as 3.4m，and the optimum spacing of drainage hole distribution was 4.8m with mathematics analysis，it references for outburst prevention development of Qianjiaying coal mine.

content index method；mine outburst prevention；drainage radius，numerical analysis

2016-11-12

10.13532/j.cnki.cn11-3677/td.2017.03.026

國(guó)家自然基金資助項(xiàng)目(51374089)

張嘉勇(1977-)，男，河南南陽(yáng)人，副教授，主要從事礦井瓦斯防治方面的研究工作。

張嘉勇， 陳 建，周鳳增，等.基于數(shù)值模擬的煤層瓦斯抽采半徑合理化試驗(yàn)研究[J].煤礦開(kāi)采，2017，22(3)：89-92.

TD712.6

A

1006-6225(2017)03-0089-04