馬忠　石必明　穆朝民

摘要：瓦斯抽采對于立井揭突出煤層起到重要的作用，準(zhǔn)確的確定鉆孔瓦斯有效抽采半徑和合理的在待抽煤層中布置抽采鉆場對煤層消突具有關(guān)鍵性作用?；诙嗫捉橘|(zhì)中流體流動達(dá)西定律理論，采用COMSOL Multiphysics軟件對該煤層瓦斯抽采進(jìn)行了模擬。模擬結(jié)果表明，此煤層的瓦斯有效抽采半徑為3 m，隨著抽采時間的增加，煤層瓦斯壓力逐漸的降低，但降低的速率會逐漸的減小。瓦斯抽采30天后，其殘余的瓦斯壓力為0.18 MPa，這與現(xiàn)場實測的最大殘余瓦斯壓力0.2 MPa相接近，這說明了模型的可信性，其模擬結(jié)果可為瓦斯抽采設(shè)計提供參考。

關(guān)鍵詞：立井揭煤；抽采半徑；瓦斯壓力；COMSOL Multiphysics

中圖分類號：TD713.37 文獻(xiàn)標(biāo)志碼：A

文章編號：1672-1098（2014）01-0083-04

我國煤礦煤與瓦斯突出災(zāi)害比較嚴(yán)重，特別是在石門或井筒揭煤的過程中[1]。對于立井揭突出煤層，需要在揭穿煤層前采取相應(yīng)的防治煤與瓦斯突出措施[2]，而瓦斯抽采便是其中一項重要工作。瓦斯抽采設(shè)計需考慮待抽煤層有效抽采半徑和抽采鉆場中瓦斯壓力的動態(tài)演化規(guī)律[3]。本文借助COMSOL Multiphysics軟件，結(jié)合現(xiàn)場實測的數(shù)據(jù)，通過計算數(shù)學(xué)模型，模擬了瓦斯在單個抽采鉆孔和抽采鉆場條件下的流動情況，得出了該煤層的有效抽采半徑，并且對多鉆孔抽采條件下瓦斯的賦存和運(yùn)移進(jìn)行了分析，這對立井揭煤過程中瓦斯抽采設(shè)計具有一定指導(dǎo)意義。

1幾何模型

以潘三礦立井揭13-1煤為研究對象，根據(jù)煤層特性參數(shù)計算此煤層的有效抽采半徑，結(jié)合現(xiàn)場施工情況，設(shè)定計算區(qū)域如圖1所示。煤層長寬各為50 m，鉆孔半徑為005 m，各鉆孔間距離為6 m。

2瓦斯流動模型

21基本假設(shè)

為了模擬出鉆孔抽采煤層瓦斯，需建立描述瓦斯運(yùn)移的方程，瓦斯流動方程是煤層瓦斯流動模型的核心，假設(shè)如下：

1） 認(rèn)為鉆孔穿透煤層，鉆孔抽采形成的流場為徑向流場。

2） 認(rèn)為吸附瓦斯含量符合朗格繆爾方程。

3） 認(rèn)為瓦斯在煤層中的流動符合達(dá)西定律。

4） 認(rèn)為瓦斯為理想氣體，滲流過程按等溫過程處理。

5） 認(rèn)為煤層各向同性，煤層滲透率為定值。

6） 認(rèn)為煤層透氣性系數(shù)和孔隙率不受煤層壓力的變化。

7） 認(rèn)為煤層瓦斯流動不受水的影響。

22建立方程

在煤層中，滲流的為游離瓦斯。吸附瓦斯對滲流場來說是質(zhì)量源，吸附瓦斯含量由朗格繆爾方程計算得到。

達(dá)西定律：

3數(shù)值模擬及結(jié)果分析

本次模擬使用COMSOL中的多孔介質(zhì)滲流模塊，結(jié)合潘三礦立井揭煤實測的煤層及瓦斯參數(shù)完成鉆孔抽采瓦斯模擬。

31模擬參數(shù)

模擬應(yīng)用潘三礦揭煤中所取煤樣測得的煤層特性參數(shù)（見表1）。

32單鉆孔抽采瓦斯流場

煤體距離鉆孔越遠(yuǎn)，所含的瓦斯壓力越大，直至接近原始瓦斯壓力（見圖2）。從鉆孔中心到模型邊緣這條直線上，瓦斯壓力隨著時間的變化，隨著抽采時間的增加，煤層瓦斯壓力逐漸的降低，鉆孔抽采瓦斯所影響的范圍逐漸增大。距離鉆孔1 m處瓦斯壓力隨時間的變化如圖3所示，瓦斯抽采的最初階段瓦斯壓力下降明顯， 隨著時間的推移， 瓦斯壓力變化逐漸放緩， 直到保持穩(wěn)定。 根據(jù)文獻(xiàn)[4]的規(guī)定， 在鉆孔有效抽采半徑內(nèi)， 需將煤層瓦斯壓力降至074 MPa以下。此煤層在抽采30 d后可知，距鉆孔3 m處的煤體的瓦斯壓力約為074 MPa，如僅依據(jù)抽采鉆孔周圍殘余瓦斯壓力值判定有效抽采半徑，那么可認(rèn)為此煤層的有效抽采半徑為3 m。

33多鉆孔抽采瓦斯流場

通過對單個鉆孔的模擬，得出了瓦斯有效抽采半徑，但是揭煤過程中，是通過多鉆孔來抽取瓦斯的，那么各鉆孔之間必然會相互影響，可以推測抽采一段時間后，鉆場中某鉆孔周圍瓦斯壓力分布必與單個鉆孔周圍瓦斯壓力分布不同。

在9個鉆孔所形成的正方形區(qū)域內(nèi)瓦斯壓力明顯低于此區(qū)域外的瓦斯壓力，說明在由鉆孔所構(gòu)成的鉆場內(nèi)，瓦斯壓力下降明顯，其抽采效果好。鉆場中瓦斯壓力隨著抽采時間的增加而逐漸減?。ㄒ妶D4）。當(dāng)煤層瓦斯抽采到25 000 s這個時刻時，從圖4中可以看出鉆場周圍瓦斯壓力梯度大，說明此時瓦斯流動迅速，大量的瓦斯流向鉆孔，而在抽采末期（259×106 s），瓦斯壓力梯度小，此時瓦斯緩慢的流向鉆孔，如此時繼續(xù)增加抽采時間則無任何意義。

L/m

1-0 s；2-2.5×104 s；3-864×105 s；4-1728×106 s；5-2592×106 s

圖4不同抽采時間瓦斯壓力分布

距鉆孔3 m處煤體瓦斯壓力隨時間變化從28 MPa下降到02 MPa以下（見圖5），比單一鉆孔抽采條件下的瓦斯壓力下降的幅度明顯增大，這說明各鉆孔在抽采過程中產(chǎn)生了相互產(chǎn)生了影響，對瓦斯抽采起到了強(qiáng)化作用。

T/Ms

圖5瓦斯壓力隨時間變化

4模擬結(jié)果與現(xiàn)場實測對比分析

根據(jù)測壓鉆孔設(shè)計參數(shù)（見表2），在距13-1煤頂板最小法距10 m前施工6個測壓孔（2-1#、2-2#、2-3#、2-4#、2-5#、2-6#），測定13-1煤原始瓦斯壓力，鉆孔穿透13-1煤進(jìn)入13-1煤底板05 m，其中2-4#孔控制到井筒輪廓線15 m外。

5結(jié)論

1） 煤層內(nèi)的瓦斯分為吸附瓦斯和游離瓦斯兩種狀態(tài)，吸附瓦斯隨游離瓦斯壓力的下降而解吸出來，進(jìn)入裂隙系統(tǒng)，參與瓦斯的流動。模擬結(jié)果與現(xiàn)場實測相近，驗證了吸附瓦斯量可通過朗格繆爾方程求得。

2） 通過建立鉆孔瓦斯流動模型，結(jié)合現(xiàn)場實測煤層相關(guān)參數(shù)，模擬得出了瓦斯抽采半徑為3 m，為煤層抽采瓦斯鉆場提供參考。

3） 模擬了多鉆孔抽采煤層瓦斯過程中，瓦斯壓力在空間和時間上的變化規(guī)律，并且這一變化規(guī)律可量化研究，對預(yù)測抽采瓦斯量具有重要意義。

4） 瓦斯抽采過程中，各鉆孔之間會產(chǎn)生相互影響，可起到強(qiáng)化瓦斯抽采的作用。

參考文獻(xiàn)：

[1]雷文杰.千米埋深特厚煤層井筒揭煤綜合防突技術(shù)[J].煤炭科學(xué)技術(shù)，2011，39（8）：50-53.

[2]俞啟香.礦井瓦斯防治[M].徐州：中國礦業(yè)大學(xué)出版社，1992：66.

[3]唐兵.鉆孔瓦斯抽采半徑的確定方法及實踐[J].礦業(yè)安全與環(huán)保，2012，39（4）：44-48.

[4]國家安全生產(chǎn)監(jiān)督管理總局.AQ1026-2006煤礦瓦斯抽采基本指標(biāo)[S].北京：煤炭工業(yè)出版社，2007.

[5]陸占金.潘三礦下保護(hù)層開采條件下地面鉆孔瓦斯抽采技術(shù)及效果[J].煤礦開采，2012，17（6）：89-91.

[6]吳厚華.“漸進(jìn)式”揭煤法在立井近水平突出煤層的實踐[J].煤炭工程，2012，12（3）：45-47.

[7]韓穎.緩傾斜突出煤層石門揭煤瓦斯災(zāi)害治理模式[J].煤礦安全，2012，43（10）：128-130.

[8]周國平.考慮基質(zhì)收縮效應(yīng)的煤層氣應(yīng)力場-滲流場耦合作用分析[J].巖土力學(xué)，2010，31（7）：2 317-2 323.

[9]王兆豐.基于COMSOL的順層鉆孔有效抽采半徑的數(shù)值模擬[J].煤礦安全，2012，43（10）：4-6.

[10]韓云龍.順層鉆孔抽放瓦斯數(shù)值模擬與應(yīng)用研究[D].淮南：安徽理工大學(xué)，2003.

[11]楊勝來， 崔飛飛， 楊思松.煤層氣滲流特征實驗研究[J].中國煤層氣， 2005，2（1）：36-39.

（責(zé)任編輯：何學(xué)華，吳曉紅）