摘 要：為了確定合適的瓦斯抽采參數(shù)，本文采用COMSOL Multiphysics多物理場耦合模擬軟件建立了煤巖-瓦斯固氣耦合數(shù)值分析模型，對瓦斯抽采鉆孔周圍的瓦斯含量變化規(guī)律進行了模擬。模擬結果與山西麥捷煤業(yè)有限公司15號煤層150205工作面傾向鉆孔瓦斯抽采現(xiàn)場實測數(shù)據(jù)進行比較，結果趨勢基本一致，進一步驗證了該模型的可行性，為推廣該模型更大范圍的應用于現(xiàn)場提供了可靠的實踐支持。

關鍵詞：鉆孔、瓦斯抽采、數(shù)值模擬

中圖分類號：TD712 文獻標志碼：A

目前，在煤炭開采深度不斷增加的過程中，瓦斯的安全問題也越來越嚴重，瓦斯含量、瓦斯涌出量突然增大和瓦斯突出成為最為棘手的問題。公認的最有效的治理和利用瓦斯的方法和途徑是鉆孔抽采瓦斯技術。

在抽放工程的施工過程中，特別是抽放鉆孔的施工過程中，抽放鉆孔的布置方式、孔間距、抽放時間等參數(shù)往往根據(jù)實踐經(jīng)驗進行設計，存在較大的隨意性和盲目性，造成抽采效果不理想和人力物力的浪費等現(xiàn)象，如何在保證安全的前提下，準確測定抽放參數(shù)，合理布置抽放工程，將抽放工程的投入和取得的經(jīng)濟效益有效結合是一個急需研究的問題。

一、數(shù)值計算軟件簡介

COMSOL Multiphysics固氣耦合模擬軟件，是一個基于偏微分方程對科學和工程問題進行建模和仿真計算的大型高級數(shù)值仿真軟件，由瑞典的COMSOL公司開發(fā)，廣泛應用于各個領域的科學研究以及工程計算，被當今世界科學家稱為“第一款真正的任意多物理場直接耦合分析軟件”，適用于模擬科學和工程領域的各種物理過程。

二、模型建立

1、幾何模型設定

幾何模型以麥捷煤業(yè)15號煤層瓦斯抽采工程為原型，同時為了保證運算速度以及模擬文件大小不超出模擬計算機實際內(nèi)存范圍，本報告將幾何模型進行如下設定：

（1）15號煤層為近水平煤層，忽略煤層傾角的影響，將幾何模型設定為水平模型；

（2）15號煤層平均厚度為3.5m。

（3）煤體范圍遠大于研究范圍，為了將無限流場模型轉化為有限流場模型，同時考慮到邊界效應，將模型設置為長為30m高為3.5m的矩形；

（4）鉆孔直徑為113mm；

（5）鉆孔破壞半徑為0.48m；

（6）封孔長度為10m。

2、模型的主要參數(shù)

3、抽放指標

依據(jù)《煤礦瓦斯抽采達標暫行規(guī)定》和《煤礦安全規(guī)程》的要求，考慮一定的安全系數(shù)，確定麥捷煤業(yè)150205工作面15號煤層瓦斯有效抽放半徑的指標為：抽采一定時間后，鉆孔周圍一定范圍內(nèi)的煤層瓦斯含量降至4m3/t以下。

4、建立數(shù)值模型

考慮到計算機運算的特點、鉆孔瓦斯抽放半徑和順層鉆孔的施工條件、鉆孔的分布對稱性，若進行二維數(shù)值模型建立的，首先必須對物理環(huán)境進行合理的簡化：

（1）假設煤層是各向同性并且連續(xù)的，而且是水平分布，對于煤體中存在的斷層和裂隙發(fā)育暫不考慮；

（2）只將截取的部分鉆孔進行模擬。

模擬區(qū)域取工作面前方受采動影響較小區(qū)域，煤層瓦斯分布均勻，鉆孔取在煤層中部，模擬鉆孔間距為2m條件下的瓦斯抽放效果。模擬計算過程按平面問題求解，模型尺寸取30m×3.5m（長×寬），模擬過程取4個鉆孔，瓦斯?jié)B透只需考慮在煤層中的滲透，處理邊界條件如下：

1）頂部巖層看作為隔絕體無滲透，左右邊界均為自由界面；

2）鉆孔的抽放負壓為20KPa。

三、數(shù)值模擬結果分析

1、數(shù)值模擬分析

本文以順層鉆孔間距為2m為例進行分析。

順層鉆孔間距為2m，在抽采20天和40天的條件下， 15號煤層鉆孔周圍瓦斯含量分布云圖如圖4-1和圖4-2所示。

從圖中可以看出，鉆孔周圍瓦斯含量分布曲線近似漏斗狀，越靠近鉆孔，瓦斯含量越低，越遠離鉆孔，瓦斯含量越高。鉆孔周圍瓦斯含量隨抽采時間的增大而有明顯降低。

鉆孔間距為2m時，鉆孔控制范圍的瓦斯含量，在抽采20天后，水平方向上降至4.5m3/t以下，垂直方向上降至4.6m3/t以下；抽采40天后，水平方向上降至4m3/t以下，垂直方向上降至4m3/t以下；

2、模擬結果

研究結論如下：

抽采40天，順層鉆孔有效抽采半徑為1.0m，鉆孔間距2.0m。

4 150205工作面瓦斯綜合治理現(xiàn)場試驗研究

150205工作面采用采前預抽的方式，在150205軌道順槽巷道內(nèi)施工本煤層鉆孔進行瓦斯抽采。在2013年11月完成全部鉆孔的施工。鉆孔間距2m，孔徑113mm，孔深130m，共布置鉆孔522個。

抽采40天后，150205工作面實測瓦斯含量如表5-1

經(jīng)分析，150205工作面實測的瓦斯含量基本與模擬結果一致，進一步驗證了該模型的可行性。

四、結論

本文根據(jù)煤巖體變形與瓦斯?jié)B流的基本理論，建立了煤巖-瓦斯固氣耦合的數(shù)學模型，利用COMSOL Multiphysics多物理場耦合軟件模擬了麥捷煤業(yè)15號煤層150205工作面瓦斯抽采參數(shù)。結果表明，抽采40天，順層鉆孔有效抽采半徑為1.0m，鉆孔間距2.0m。

根據(jù)上述理論分析與數(shù)值模擬結果，在麥捷煤業(yè)150205工作面進行了現(xiàn)場試驗研究，抽采后工作面實測瓦斯含量結果表明，數(shù)值模擬結果基本正確，符合實踐要求，該模型具有推廣可行性。

參考文獻：

[1] 周世寧，林柏泉.《煤層瓦斯賦存與流動理論》[M]，北京：煤炭工業(yè)出版社，1999.2.

[2] 郭培紅，李海霞，朱建安.《煤層鉆孔瓦斯抽放數(shù)值模擬》[J]，遼寧工程技術大學學報（自然科學版），2009，28：260-262.

[3] 馬耕，蘇現(xiàn)波，魏慶喜.《基于瓦斯流態(tài)的抽放半徑確定方法》[J]，煤炭學報，2009，34：501-504.

[4] 林柏泉，張建國.《礦井瓦斯抽放理論與技術》[M]，徐州：中國礦業(yè)大學出版社，1996.8.165

[5] 俞啟香.《礦井瓦斯防治》[M]，徐州：中國礦業(yè)大學出版社，1992.2.

[6] 吳世躍.《煤層中的耦合運動理論及其應用》[M]，北京：科學出版社，2009.3.

[7] 楊宏民，魏晨慧，王兆豐.《基于多物理場耦合的井下注氣驅(qū)替煤層甲烷的數(shù)值模擬》[J]，煤炭學報，2010，35：109-114.

[8] 劉鶴年.《流體力學》[M]，北京：中國建筑工業(yè)出版社，2001.6.昂溫 G，昂溫 P S. 陳生錚，譯. 外出版史. 北京：中國書籍出版社，1988： 12-15