煤層注水促抽瓦斯及其影響因素的數值模擬

陳紹杰， 陳舉師， 汲銀鳳， 李改革， 周逸飛

(1. 北京科技大學 土木與資源工程學院，北京 100083； 2. 華北科技學院 安全工程學院，北京 101601)

煤層注水促抽瓦斯及其影響因素的數值模擬

陳紹杰1,2， 陳舉師1， 汲銀鳳1， 李改革2， 周逸飛2

(1. 北京科技大學 土木與資源工程學院，北京 100083； 2. 華北科技學院 安全工程學院，北京 101601)

為有效預防煤礦瓦斯災害，獲取煤層注水促抽瓦斯的合理參數，以常村煤礦2103工作面為例，依據多相滲流理論，采用Fluent軟件的VOF模型及多孔介質模型耦合求解，對煤層注水促抽瓦斯技術及其影響因素進行數值模擬，并將模擬結果應用于現場，對比分析數值模擬與現場實測數據，二者基本吻合. 研究結果表明：煤層瓦斯含量以注水孔為中心徑向逐步降低，以抽采孔為中心徑向逐步升高；注水前抽采階段，隨著抽采時間的增加，抽采范圍逐漸增大，抽采孔瓦斯流量先快速下降，后逐步緩慢降低；注水促抽階段，隨著注水時間的增加，注水范圍逐漸增大，注水流量逐步降低，煤層瓦斯含量緩慢升高，抽采孔瓦斯流量逐漸增加；注水后抽采階段，隨著抽采時間的增加，壓力水覆蓋范圍持續增大，煤層瓦斯含量逐漸降低，抽采孔瓦斯流量逐漸減小. 注水時機、注水時間、注水壓力、注水方式、布置方式及鉆孔間距是影響煤層注水促抽瓦斯效果的6個主要因素. 瓦斯正常抽采20 d后，按照一注一抽方式及5 m間距布置注抽鉆孔，在8 MPa煤層注水壓力下間歇注水10 d，煤層注水促抽瓦斯效果較好.

煤層注水；多相滲流；瓦斯抽采；影響因素；數值模擬

煤礦開采過程中，瓦斯時刻威脅著井下工人的生命安全，我國煤礦的重大災害中約70%都是瓦斯事故[1]. 近年來，隨著我國煤礦開采深度逐漸向下延伸，煤層瓦斯含量和涌出量均呈現出大幅增加的趨勢，亟待采取有效措施防患瓦斯災害于未然[2]. 煤層瓦斯抽采是一條防治瓦斯災害的積極有效措施，不僅可以降低煤層瓦斯含量，還可以消除和預防煤與瓦斯突出[3].

一般來說，煤層瓦斯抽采效率主要取決于煤層瓦斯壓力及煤層滲透性. 煤層滲透性越差，瓦斯在煤層中滲流難度越大，抽采效率越低[4]. 我國煤層賦存條件極其復雜，部分松軟低透氣煤層由于原始應力區透氣性極低，需采取相關輔助措施方能實現有效抽采[5]. 目前提高本煤層瓦斯抽采效率的途徑主要有兩種：一是采取人為方法提高煤層的透氣性，二是合理布置鉆孔和改變鉆孔參數[6]. 而采取高壓注水驅替瓦斯具有一定的實際應用效果，通過將高壓水注入煤層驅替孔隙及裂隙中的瓦斯，迫使瓦斯沿煤層孔隙及裂隙通道往抽采孔滲透，從而實現瓦斯的有效抽采[7]. 但由于該技術在工藝參數選擇上尚未形成統一認識，嚴重制約該技術的發展[8].

因此，本文通過對煤層注水促抽瓦斯效果及其影響因素進行數值模擬研究，掌握煤層注水條件下瓦斯抽采流量的變化規律，獲取最佳的工藝參數，可為煤層注水促抽瓦斯現場應用提供技術支持，對于防治煤礦瓦斯災害具有十分重大的指導意義.

1 煤層注水驅替瓦斯機理

煤層可看做是由孔隙及裂隙結構組成的多孔介質. 煤層注水時，一方面水在高壓作用下進入煤層裂隙，并驅逐裂隙通道中的瓦斯向低壓區滲透；另一方面，裂隙中水在孔隙毛管力的吸引下，會沿著基質煤塊固體顆粒表面進入孔隙，并將瓦斯從孔隙通道中驅逐而出，進入裂隙通道. 隨著注水過程的推進，煤層孔隙和裂隙通道將會被水大量占據，只殘留少量低于束縛飽和度的瓦斯. 孔隙通道中瓦斯進入裂隙通道后，在裂隙內即存在水與瓦斯兩相流體，形成高壓作用下沿裂隙通道的兩相滲流. 煤層內水-瓦斯兩相滲流與瓦斯單相滲流的交界區域內，所有瓦斯飽和度為零的質點組成了水與瓦斯突變界面，而研究煤層注水驅替瓦斯問題的關鍵就是要追蹤該界面的運動情況，以此來確定煤層中水與瓦斯含量的分布及變化規律.

2 控制方程

煤層注水促抽瓦斯模擬采用VOF模型追蹤水與瓦斯突變的自由界面，確定該自由界面在煤層中的運動規律，通過求解其中一相體積分數方程和多孔介質的動量方程，掌握水與瓦斯兩相流體在煤層中的滲流特征.

2.1體積分數方程

對于煤層注水驅替瓦斯來說，假設瓦斯和水所占體積分數分別為α1和α2，則可得出瓦斯體積分數守恒方程為

(1)

由于煤層中瓦斯和水體積分數之和等于1，則煤層中混合流體的密度ρ為

ρ=α1ρ1+(1-α1)ρ2.

(2)

式中：ui為i方向上的速度，xi為i方向上的坐標，t為時間,ρ為混合流體密度,ρ1為瓦斯密度,ρ2為水密度.

煤層注水驅替瓦斯過程其他參數如粘性系數等均可按上述方法計算. 在不同時刻下，對瓦斯流場的體積分數方程進行求解，可得出水與瓦斯體積分數分布情況，再運用相關方法和手段重構運動界面，即可追蹤兩相自由界面在煤層中隨時間的變化.

2.2動量守恒方程

通過求解煤層中混合流體的動量守恒方程，可得出混合流體的速度場，結合瓦斯與水所占體積分數，可求出水和瓦斯的速度分布. 動量守恒方程為

(3)

煤層可看做由孔隙和裂隙結構組成的均勻多孔介質，且各向同性，則其具有附加的動量源項：

(4)

由于煤層中流體流動表現為層流，屬于線性變化的達西滲流，粘性阻力起主要作用，內部阻力可忽略不計. 根據達西公式，則可進一步將煤層的附加動量源項簡化為

(5)

根據Blake-Kozeny半經驗公式：

(6)

則可得出煤層滲透性系數為

(7)

式中：uj為j方向上的速度，xj為j方向上的坐標，p為壓力，μ為動力粘性系數，gi為i方向上的重力加速度，K為煤層滲透性系數，C2為內部阻力系數，Dp為顆粒平均直徑，n為孔隙率.

3 幾何模型的建立及求解

3.1工程概況

常村煤礦為高瓦斯礦井，2103工作面位于+470水平21采區，開采3#煤層，埋藏深度為+462.5～+504.1 m，主應力方向北東向，平均N40°E，最大垂直應力11 MPa，最大水平應力13 MPa. 工作面呈矩形，走向長度724 m，傾斜長度294 m，斜面積187 046 m2. 煤層平均煤厚6.4 m，采高3.3±0.1 m，采放比1.1∶1，煤體容重1.4 t/m3. 工作面老頂為粉砂巖，直接頂為細粒砂巖，直接底為中粒砂巖，老底為細粒砂巖. 2103工作面所采3#煤層本質上屬于松軟低滲透性煤層，瓦斯含量為7.6 m3/t，可解吸量為5.7 m3/t，煤塵具有爆炸性，無自燃現象，地溫小于26 ℃，屬溫度正常區，無熱害威脅.

3.2幾何模型建立及網格劃分

煤層注水促抽瓦斯屬于三維非穩態問題，對于計算機性能及計算時長要求均比較高，綜合考慮計算機配置及數值模擬效率，采用二維模型進行煤層注水促抽瓦斯模擬，并選取鉆孔徑向平面開展研究. 但由于煤層本身賦存條件極其復雜，注水促抽工藝影響因素較多，現場情況復雜多變，要完全復制現場細節難度較大. 因此，需結合煤層注水促抽瓦斯工藝參數設計及現場布置情況，對本次數值模擬做出如下簡化及假設：

1)煤層頂底板不具備滲透性且不含瓦斯；瓦斯在煤層內流動狀態為層流滲流，服從達西定律；煤層中瓦斯吸附特性符合朗格繆爾方程且解吸在瞬間完成；

2)煤層沿鉆孔徑向各向滲透性保持不變；沿鉆孔軸向滲透性可變，可通過設置不同滲透性系數模擬沿鉆孔軸向各平面煤層滲流情況；

3)煤層瓦斯視為可壓縮氣體，密度可隨壓力變化；水不可壓縮，密度保持不變；忽略滲流過程水與瓦斯溫度變化，按等溫過程處理.

基于上述簡化及假設，運用Gambit建立計算域為15 m×6.4 m的煤層注水促抽瓦斯二維幾何模型. 模型建立3個鉆孔，鉆孔直徑均為0.094 m，鉆孔間距5 m，其中，2#鉆孔兼做抽采孔及注水孔，坐標原點位于模型中心. 采用由線至面的方式對該幾何模型進行網格劃分，并對鉆孔附近區域網格進行局部加密. 煤層注水促抽瓦斯模型及網格劃分見圖1.

圖1 煤層注水促抽瓦斯幾何模型

Fig.1 2D geometric model of coal seam water infusion promoting methane drainage

3.3模擬參數的設定及求解

通過廣泛查閱煤層注水促抽瓦斯相關文獻及技術資料，結合現場實際測定數據，遵循Fluent軟件中多孔介質模型及VOF模型的應用條件，確定出煤層注水促抽瓦斯模擬所必須的參數并進行設定，最終求解出煤層中水與瓦斯的滲流規律，從而獲得煤層注水促抽瓦斯的合理工藝參數. 數值模擬參數設定見表1[9-15].

表1 計算模型參數設定

4 數值模擬結果及分析

4.1注水前煤層瓦斯抽采效果及分析

為探索注水前煤層瓦斯抽采效果，保持1#、2#及3#抽采孔正常抽采20 d，在Fluent中設置抽采負壓8 000 Pa，抽采孔徑0.094 m，鉆孔間距5 m，煤層滲透性系數根據常村煤礦2103工作面實驗數據取0.185×10-15m2，最終求解出煤層瓦斯隨時間滲流規律. 圖2為煤層注水前不同抽采時間下瓦斯含量分布，圖3為鉆孔中心線瓦斯含量分布. 由圖2、圖3及圖8可知：

1)抽采條件下，煤層瓦斯含量以抽采孔為中心徑向逐步升高；隨著抽采時間的增加，抽采范圍逐漸增大，煤層瓦斯含量及抽采孔瓦斯流量均呈現出先快速下降，后逐步緩慢降低的趨勢.

2)當抽采時間t抽1由1 h增加至20 d時，煤層瓦斯含量由1.25～5.5 m3/t逐步下降至0.65～1.14 m3/t，鉆孔瓦斯抽采流量由3.75×10-5m3/s逐步降低至3.50×10-7m3/s.

3)抽采過程中，煤層瓦斯在原始瓦斯壓力及抽采孔負壓共同作用下，沿徑向持續往抽采孔滲透；隨著抽采過程的推進，瓦斯抽采總量不斷累積，煤層瓦斯總量不斷減少，導致瓦斯壓力不斷降低，供給瓦斯滲流的壓差逐漸減小，抽采孔瓦斯流量逐步降低.

(a) t抽1=1 h～2 d煤層瓦斯含量分布

(b) t抽1=5 d～20 d煤層瓦斯含量分布

Fig.2 Distribution of methane content in coal seam before water infusion

圖3 注水前鉆孔中心線瓦斯含量分布

Fig.3 Distribution of methane content in the center line of the hole before water infusion

4.2注水時煤層瓦斯抽采效果及分析

為研究煤層注水促進瓦斯抽采效果，在煤層正常抽采20 d的基礎上，按照一注兩抽的方式間隔布置注水孔和抽采孔(2#作為注水孔，1#、3#繼續作為抽采孔)，連續注水10 d，注水壓力8 MPa，其余參數與注水前保持一致. 圖4、5分別為注水時煤層及鉆孔中心線上瓦斯含量分布，由圖4、5、9可知：

1)注水促抽瓦斯條件下，煤層瓦斯含量以注水孔為中心徑向逐步降低，以抽采孔為中心徑向逐步升高；隨著注水時間的增加，注水范圍逐漸增大，注水流量逐步降低，煤層瓦斯含量緩慢升高，抽采孔瓦斯流量逐漸增加.

2)當注水時間t注由1 h增加至10 d時，注水流量由3.96×10-6m3/s逐步下降至1.85×10-6m3/s；煤層瓦斯含量由0.65～51.82 m3/t逐步升高至0.67～53.89 m3/t，瓦斯流量由3.50×10-7m3/s逐步增大至1.50×10-6m3/s.

3)注水促抽瓦斯過程中，壓力水自注水孔進入煤層后，幾乎占據了煤層中所有的主要通道，只留下少量低于束縛飽和度的瓦斯分布在孔隙和裂隙的邊、角及基質煤塊表面；隨著注水進程的推移，煤層內水滲流阻力不斷增大，導致注水流量逐漸降低；在注水范圍內，處于分散相的瓦斯被高壓水不斷壓縮，導致該區域瓦斯含量逐步升高；在注水范圍外，大量瓦斯被壓力水驅逐進入，導致該區域瓦斯總量增加，瓦斯含量不斷升高.

圖4 注水時煤層瓦斯含量分布

Fig.4 Distribution of methane content in coal seam when water infusion

圖5 注水時鉆孔中心線瓦斯含量分布

Fig.5 Distribution of methane content in the center line of the hole when water infusion

4.3注水后煤層瓦斯抽采效果及分析

為掌握注水結束后煤層瓦斯抽采效果，在上一階段持續注水10 d的基礎上，關閉注水閥門，停止注水，僅留下1#、3#鉆孔繼續抽采. 圖6、圖7分別為注水后煤層及鉆孔中心線上瓦斯含量分布，由圖6、圖7及圖9中可知：

1)注水后抽采條件下，煤層瓦斯含量分布規律和注水時基本保持一致，但數值上有較大程度的降低；隨著抽采時間的增加，壓力水覆蓋范圍持續增大，煤層瓦斯含量逐漸降低，抽采孔瓦斯流量逐漸減小.

2)當抽采時間t抽2由1 d增加至10 d時，煤層瓦斯含量由0.66～2.27 m3/t逐步降低至0.64～0.95 m3/t，瓦斯流量由1.20×10-6m3/s逐步減小至2.35×10-7m3/s.

3)注水后抽采過程中，壓力水繼續在煤層內往抽采孔方向滲透，注水范圍內壓力不斷衰減，水的飽和度不斷降低，瓦斯飽和度不斷增加，覆蓋范圍逐漸增大，導致該區域瓦斯含量逐步降低；在注水范圍外，雖然壓力水的持續驅逐作用會在一定程度上增加該區域瓦斯總量，但其增長速率遠遠低于瓦斯抽采速率，因此該區域內瓦斯含量也逐步降低；而抽采孔瓦斯流量也由于水壓力衰減、瓦斯抽采總量增加等原因逐步降低.

圖6 注水后煤層瓦斯含量分布

Fig.6 Distribution of methane content in coal seam after water infusion

圖7 注水后鉆孔中心線瓦斯含量分布

Fig.7 Distribution of methane content in the center line of the hole after water infusion

4.4煤層注水促抽瓦斯影響因素及分析

煤層注水促抽瓦斯工藝參數復雜，影響因素較多，目前對于如何合理、有效地選取煤層注水促抽瓦斯的技術參數尚未形成統一的認識. 為此，通過廣泛查閱相關文獻資料，結合該項技術的實際特點，選取了8個較為常見的技術參數進行模擬，以期獲得煤層注水促抽瓦斯的主要影響因素，并確定一個比較合理的工況范圍. 由于影響因素較多，開展正交模擬難度較大，本次模擬選擇單因素分析法對注水時機、注水時間、注水壓力、注水方式、鉆孔布置、鉆孔間距、鉆孔直徑及煤層滲透率等8個因素逐一進行對比.

4.4.1 注水時機

保持其余參數不變，分別對抽采t0=0 d、10 d、20 d及30 d后開始注水的促抽瓦斯效果進行模擬，不同注水時機下抽采孔瓦斯流量隨時間變化見圖8，可知：

不同注水時機下煤層注水對于瓦斯抽采均有一定促進作用，且注水時機越早，抽采孔內瓦斯流量增加值越大，促抽效果越好；在實際工程應用中，開始階段瓦斯抽采效率均比較高，從充分發揮煤層注水工藝優勢的角度出發，可考慮在抽采一段時間后，待抽采流量降低至一個較小值時，再介入煤層注水. 當抽采20 d后，抽采孔瓦斯流量下降至約3.50×10-7m3/s，瓦斯抽采效率較低，此時介入煤層注水較為合理.

圖8 不同注水時機下抽采孔瓦斯流量隨時間變化

Fig.8 Methane flow of drainage hole variations with time under different water injection timings

4.4.2 注水時間

待抽采20 d后，分別對注水時間t=0 d、5 d、10 d及15 d時煤層注水促抽瓦斯效果進行模擬，不同注水時間下抽采孔瓦斯流量隨時間變化見圖9，可知：

隨著注水時間的增加，抽采孔瓦斯流量亦不斷增加，但當注水時間達到10 d后，瓦斯流量開始迅速降低，且注水達到13 d后，抽采孔內瓦斯流量開始上下劇烈波動，這是由于壓力水滲流至抽采孔內導致的. 因此，從煤層注水促抽瓦斯效果及工藝安全方面考慮，在實際工程應用中，注水時間選取10 d比較合適.

圖9 不同注水時間下抽采孔瓦斯流量隨時間變化

Fig.9 Methane flow of drainage hole variations with time under different water injection time

4.4.3 注水壓力

在抽采20 d后注水10 d的工藝基礎上，分別對注水壓力P=4、6、8、10、12 MPa下煤層注水促抽瓦斯效果進行模擬，不同注水壓力下抽采孔瓦斯流量隨時間變化見圖10，可知：

在本次模擬所選擇幾何模型及參數設置條件下，隨著注水壓力的升高，抽采孔瓦斯流量不斷升高，但當注水壓力為10 MPa及12 MPa時，瓦斯流量分別在注水時間8 d及7 d時開始迅速降低，且在注水10 d及9 d后開始劇烈波動. 因此，從煤層注水促抽瓦斯效果及工藝安全方面考慮，在實際工程應用中，注水壓力選取8 MPa較為合適.

圖10 不同注水壓力下抽采孔瓦斯流量隨時間變化

Fig.10 Methane flow of drainage hole variations with time under different water injection pressure

4.4.4 注水方式

在以往煤層注水的經驗中，常采用連續及間歇兩種方式將壓力水注入煤層. 本文在保證注水總時長為10 d的基礎上，分別對連續注水及間歇注水(注8 h，停8 h)兩種注水方式促抽瓦斯效果進行模擬，不同注水方式下抽采孔瓦斯流量隨時間變化見圖11，可知：

間歇注水方式下抽采孔瓦斯流量增加值較之連續注水同期要高出將近1倍，而注水總時長卻僅為連續注水的一半. 因此，從煤層注水促抽瓦斯效果和節約能耗方面考慮，在實際工程應用中，采用間歇注水更為恰當.

圖11 不同注水方式下抽采孔瓦斯流量隨時間變化

Fig.11 Methane flow of drainage hole variations with time under different water injection modes

4.4.5 布置方式

本文主要采用注抽鉆孔間隔布置的方式開展煤層注水促抽瓦斯效果研究，分別對兩注一抽、一注一抽及一注兩抽三種注抽鉆孔間隔布置方式進行模擬，不同布置方式下抽采孔瓦斯流量隨時間變化見圖12，可知：

三種布置方式均能取得一定的促抽效果，且兩注一抽、一注一抽及一注兩抽布置方式分別能將單孔抽采效率提高18.6、11.4及4.3倍. 其中，兩注一抽布置方式下容易引發壓力水滲流至抽采孔內，安全性較差. 本文結合常村煤礦2103工作面實際情況，并充分考慮節約能耗、提高整體抽采效率等因素，選擇一注一抽方式布置注抽鉆孔開展現場工業性試驗.

圖12 不同布置方式下抽采孔瓦斯流量隨時間變化

Fig.12 Methane flow of drainage hole variations with time under different arrangement form

4.4.6 鉆孔間距

為了解鉆孔間距對煤層注水促抽瓦斯效果的影響，分別對鉆孔間距d=3、4、5、6及7 m條件下注水促抽瓦斯進行模擬，不同鉆孔間距下抽采孔瓦斯流量隨時間變化見圖13，可知：

鉆孔間距越大，抽采階段瓦斯抽采流量越高，注水階段瓦斯抽采流量升高幅度越小，注水促抽瓦斯效果越差. 當鉆孔間距為3 m及4 m時，注水過程中壓力水會滲流至抽采孔內，極不安全. 因此，推薦選擇5 m間距進行鉆孔布置.

圖13 不同鉆孔間距下抽采孔瓦斯流量隨時間變化

Fig.13 Methane flow of drainage hole variations with time under different distances between the drill holes

4.4.7 鉆孔直徑

為研究鉆孔直徑對煤層注水促抽瓦斯效果的影響，分別對鉆孔直徑D=0.075、0.094、0.113、0.133及0.153 m條件下注水促抽瓦斯進行模擬，不同鉆孔直徑下抽采孔瓦斯流量隨時間變化見圖14，可知：

隨著鉆孔直徑的增大，抽采孔瓦斯流量略微有所增加，但彼此之間差距并不大. 由此可見鉆孔直徑對于煤層注水促抽瓦斯影響較小，在實際工程應用中，可以不考慮鉆孔直徑因素的影響.

圖14 不同鉆孔直徑下抽采孔瓦斯流量隨時間變化

Fig.14 Methane flow of drainage hole variations with time under different hole diameters

4.4.8 煤層滲透性

煤層滲透率決定著流體在煤層內的滲流情況，本文分別選擇煤層滲透率K=0.046×10-15、0.104×10-15、0.185×10-15、0.289×10-15及0.417×10-15m2條件下注水促抽瓦斯進行模擬，不同煤層滲透率下抽采孔瓦斯流量隨時間變化見圖15，可知：

在抽采階段，煤層滲透率越大，抽采孔瓦斯流量越大；在注水階段，煤層滲透率越大，注水促抽瓦斯效果越好，但也越容易導致壓力水滲流至抽采孔內，安全性能越差. 因煤層滲透率是煤層固有特性參數，不能列為煤層注水促抽瓦斯工藝的主要影響因素，在實際工程應用中，可直接根據煤層相關特性實驗確定.

圖15 不同煤層滲透性下抽采孔瓦斯流量隨時間變化

Fig.15 Methane flow of drainage hole variations with time under different coal seam permeabilities

5 煤層注水促抽瓦斯的現場試驗

5.1現場試驗方案

本次試驗在常村煤礦+470水平21采區2103皮帶順槽展開，選擇煤層瓦斯賦存穩定、無地質構造、頂底板透氣性較差，且未抽采區域施工3個抽采鉆孔(1#、2#和3#)，其中2#鉆孔兼做抽采孔和注水孔，注抽鉆孔布置方式按照一注一抽進行.

5.2工藝參數選擇

根據前期數值模擬結果，結合煤層注水工程實踐，施工3個間距為5 m的鉆孔，孔徑0.094 m，設計鉆孔深度80 m，封孔長度8 m，封孔后連接抽采系統進行抽采，抽采負壓8 000 Pa. 當抽采20 d后，將2#鉆孔作為注水孔開始注水，采用間歇注水方式(注8 h，停8 h)，注水壓力8 MPa. 當注水10 d后，關閉注水閥門，停止注水，1#、3#鉆孔繼續抽采10 d.

5.3數據測定及分析

整個試驗期間，每天按時測定1#和3#抽采孔混合瓦斯流量、瓦斯濃度及鉆孔負壓. 為了對比煤層注水促抽瓦斯數值模擬結果的準確性，取相同參數設置及邊界條件下模擬結果與現場試驗結果進行對比分析. 但由于本次模擬采用二維模型，抽采孔瓦斯流量為單位厚度煤層流量(鉆孔軸向長度為1 m)，需乘以鉆孔有效抽采長度72 m(鉆孔深度與封孔長度之差)后才具有可比性. 由于1#鉆孔及3#鉆孔各方面參數均相同，故只選擇1#鉆孔瓦斯抽采流量數據進行對比分析，煤層注水促抽瓦斯現場試驗及數值模擬對比見圖16,可知：

煤層注水促抽瓦斯模擬結果及實測數據在規律上基本保持一致，在抽采階段(0～20 d)，抽采孔瓦斯流量均隨時間推移逐步降低；在注水階段(20～30 d)，隨注水時間的增加，瓦斯抽采流量也迅速增加，注水促抽瓦斯效應顯著；注水后抽采階段(30～40 d)，瓦斯流量均隨時間推移逐步降低. 經仔細對比，模擬結果與實測數據在具體數值上還是略微有所偏差. 據分析，造成此類偏差的主要原因是由于現場封孔不嚴、數據測定誤差較大及模擬采用二維代替三維等因素造成的. 通過本文的研究，能看出煤層注水對于瓦斯抽采具有明顯促進作用，模擬所得出的結果對于現場煤層注水促抽瓦斯試驗具有較強的指導意義.

圖16 模擬結果與試驗數據對比

Fig.16 Comparison between simulated model and field measurement

6 結 論

1)注水前抽采階段，煤層瓦斯含量及壓力均以抽采孔為中心徑向逐步升高；隨著抽采時間的增加，抽采范圍逐漸增大，煤層瓦斯含量及抽采孔瓦斯流量均呈現出先快速下降，后逐步緩慢降低的趨勢.

2)注水促抽階段，煤層瓦斯含量以注水孔為中心徑向逐步降低，以抽采孔為中心徑向逐步升高；隨著注水時間的增加，注水范圍逐漸增大，注水流量逐步降低，煤層瓦斯含量緩慢升高，抽采孔瓦斯流量逐漸增加.

3)注水后抽采階段，煤層瓦斯含量分布規律和注水時基本保持一致，但數值上有較大程度的降低；隨著抽采時間的增加，壓力水覆蓋范圍持續增大，煤層瓦斯含量逐漸降低，抽采孔瓦斯流量逐漸減小.

4)注水時機、注水時間、注水壓力、注水方式、布置方式及鉆孔間距是影響煤層注水促抽瓦斯效果的6個主要因素；抽采20 d后，按照一注一抽方式及5 m間距布置注抽鉆孔，在8 MPa注水壓力下間歇注水10 d，煤層注水促抽瓦斯效果較好.

[1] 孟艷軍，湯達禎，許浩，等．煤巖氣水相對滲透率研究進展與展望[J]．煤炭科學技術，2014，42(8)：51-55．

MENG Yanjun, TANG Dazhen, XU Hao, et al. Progress and prospect of gas-water relative permeability of coal and rock[J]. Coal Science and Technology, 2014, 42(8): 51-55.

[2] JOHN H, STUART J, JOHN S, et al. Methane capacities of Bowen Basin coals related to coal properties [J]. Fuel, 1997, 76(9): 813-819.

[3] 姜楠楠，金龍哲，陳紹杰．煤層注水驅替瓦斯前緣的滲流規律研究[J]．煤礦安全，2015，46(4)：16-19．

JIANG Nannan, JIN Longzhe, CHEN Shaojie. Study on seepage law for the front of gas displacement by coal seam water injection[J]. Safety in Coal Mines, 2015, 46(4): 16-19.

[4] 宋維源，李大廣，章夢濤，等．煤層注水的水氣驅替理論研究[J]．中國地質災害與防治學報，2006，17(2)：147-150．

SONG Weiyuan, LI Daguang, ZHANG Mengtao, et al. The study on the theory of water and gas replacing during water injection into coal seams[J]. The Chinese Journal of Geological Hazard and Control, 2006, 17(2): 147-150.

[5] WACLAW D, ANDRZEJ K. Mathematical model of methane emission caused by a collapse of rock mass erump[J]. Archives of Mining Sciences, 2001, 46(4): 433-449.

[6] 程慶迎．低透煤層水力致裂增透與驅趕瓦斯效應研究[D]．徐州：中國礦業大學，2012．

CHENG Qingying. Research on permeability improvement and methane driven effect of hydraulic fracturing for low permeability coal seam[D]. Xuzhou: China University of Mining and Technology, 2012.

[7] 郜國肖，白云龍，侯濤，等．煤層非穩態法氣和水相對滲透率實驗研究[J]．煤礦安全，2015，46(5)：12-15．

GAO Guoxiao, BAI Yunlong, HOU Tao, et al. Experiment on relative permeability of gas and water based on coal unsteady-state method[J]. Safety in Coal Mines, 2015, 46(5): 12-15.

[8] 呂祥鋒，潘一山，劉建軍，等．煤層氣-水兩相流滲透率測定實驗研究[J]．水資源與水工程學報，2010，21(2)：29-32．

LV Xiangfeng, PAN Yishan, LIU Jianjun, et al. Experiment on the permeability rate in the gas-water double state flow of coal deposit[J]. Journal of Water Resources and Water Engineering, 2010, 21(2): 29-32.

[9] CLARKSON C, BUSTIN R. Binary gas adsorption/ desorption isotherms: effect of moisture and coal composition upon carbon dioxide selectivity over methane [J]. International Journal of Coal Geology, 2000,42(4): 241-272.

[10]丁厚成，蔣仲安，韓云龍．順煤層鉆孔抽放瓦斯數值模擬與應用[J]．北京科技大學學報，2008，30(11)：1205-1210．

DING Houcheng, JIANG Zhongan, HAN Yunlong. Numerical simulation and application of boreholes along coal seam for methane drainage[J]. Journal of University of Science and Technology Beijing, 2008, 30(11): 1205-1210.

[11]韋善陽，陳學習，董利輝．順層鉆孔抽采煤層瓦斯有效影響半徑數值模擬及應用[J]．中國礦業，2015，24(8)：101-105．

WEI Shanyang, CHEN Xuexi, DONG Lihui. Numerical simulation and test of effective influence radius of gas drainage boreholes[J]. China Mining Magazine, 2015, 24(8): 101-105.

[12]RUPPEL T C. Adsorption of methane/ethane mixtures on dry coal at elevated pressures[J]. Fuel, 1972,51(10):297-303.

[13]徐明智，李希建．煤層瓦斯抽放半徑及其影響因素的數值模擬[J]．工業安全與環保，2012，38(12)：28-30．

XU Mingzhi, LI Xijian. Numerical simulation of gas drainage radius and its influencing factors in coal seam[J]. Industrial Safety and Environmental Protection, 2012, 38(12): 28-30.

[14]MAZUMDER S, HEMERT P, BRUINING J, et al. In-situ CO2-coal reaction in view of carbon dioxide storage in deep unminable coal seams[J]. Fuel, 2006,85(9):1904-1912.

[15]楊新樂，任常在，張永利．低滲透煤層氣注熱開采熱-流-固耦合數學模型及數值模擬[J]．煤炭學報，2013，38(6)：1044-1049．

YANG Xinle, REN Changzai, ZHANG Yongli. Numerical simulation of the coupled thermal-fluid-solid mathematical models during extracting methane in low-permeability coal bed by heat injection[J]. Journal of China Coal Society, 2013, 38(6): 1044-1049.

Numericalsimulationofcoalseamwaterinfusionpromotingmethanedrainageanditsinfluencefactors

CHEN Shaojie1,2, CHEN Jushi1, JI Yinfeng1, LI Gaige2, ZHOU Yifei2

(1. School of Civil & Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China;2.School of Safety Engineering, North China Institute of Science and Technology, Beijing 101601, China)

In order to prevent the methane disaster in coal mine effectively and obtain the reasonable parameters of promoting methane extraction by coal seam water infusion, the 2103 working face of Changcun coal mine was investigated as an example. Using the coupling of VOF model of Fluent software and porous medium model, promoting methane drainage by coal seam water infusion and its influence factors was simulated based on the multi-phase porous flow theory. The comparison shows that the simulation results were basically consistent with the measured data. The results are detailed as follow. The gas content in coal seam gradually decreased with the increase of the radial distance from centerline of water injection hole and it increased with the increase of the radial distance from centerline of gas drainage hole. During drainage period before water injection, the scope of drainage increased gradually and gas flow in drainage hole decreased quickly at first, and then decreased slowly. During the stage of water injection, the scope of drainage increased and water injection flow decreased gradually, while gas flow in coal seam and gas flow in drainage hole also increased gradually over time. And after water injection, the coverage area of pressure water continued to increase, while gas content in coal seam and gas flow in drainage hole both declined little by little over time. The opportunity, time, pressure, layout, pressure, space between boreholes and the way of water injection are the 6 main factors which influence the effect of promoting methane extraction by coal seam water infusion. It can achieve the best effect when intermittent water injection was carried out for 10 days under the pressure of 8 MPa water injection according to the way of one injection and one pumping and arranging the drilling holes spacing at 5 m after 20 days of the methane drainage.

coal seam water infusion; multi-phase porous flow; gas drainage; influence factors; numerical simulation

10.11918/j.issn.0367-6234.201703016

TD713

A

0367-6234(2017)11-0087-08

2017-03-04

中央高校基本科研業務費項目(3142015105); 國家自然科學基金(51074015); 國家自然科學基金青年科學基金項目(51604018); 北京市自然科學基金青年科學基金項目(8164060)

陳紹杰(1981—)，男，博士研究生，副教授

陳舉師，chenjushi@163.com

(編輯苗秀芝)