基于COMSOL的采空區(qū)埋管抽采參數(shù)優(yōu)化數(shù)值研究*

楊前意，羅伙根，石必明，張雷林，張鴻智

(1. 安徽理工大學 能源與安全學院，安徽 淮南 232001； 2. 中國神華神東煤炭集團保德煤礦，山西 忻州 036600)

0 引言

瓦斯所導致的煤礦安全事故頻繁發(fā)生，對礦山安全生產(chǎn)帶來巨大威脅，針對煤礦開采中局部地區(qū)瓦斯積聚，特別是工作面上隅角瓦斯超限所造成的各種問題，采用埋管抽采瓦斯的治理方法，通過采空區(qū)埋管抽采，使得工作面上隅角瓦斯?jié)舛认陆担_到國家安全規(guī)程標準。

國內(nèi)外學者對采空區(qū)瓦斯運移及埋管抽采參數(shù)進行了大量研究。Matuszewski采用一維模型對采空區(qū)瓦斯?jié)舛确植技捌錅y試開始對采空區(qū)的氣流場進行研究。劉劍[1]采用三維有限單元法來解算采場瓦斯?jié)舛确植迹瑫r介紹了計算機運行程序；丁廣驤等[2]通過線性滲流理論以及管路漏流流動規(guī)律，建立了采場風流參數(shù)計算的二維數(shù)學模型；張東明等[3]研究了采空區(qū)瓦斯流動分布規(guī)律，并運用有限元方法建立了非線性滲流數(shù)學模型；周愛桃等[4]設計了采空區(qū)埋管抽放的技術(shù)方案，建立了采空區(qū)瓦斯運移的數(shù)值計算模型，并通過FLUENT軟件對采空區(qū)埋管抽放效果進行了模擬；趙宇新[5]研究了綜采工作面及采空區(qū)瓦斯分布情況和運移規(guī)律，并分別闡明了煤層賦存條件和開采技術(shù)工藝對瓦斯涌出量產(chǎn)生的各種影響，建立三維數(shù)學模型，運用FLUENT數(shù)值模擬軟件對常易造成積聚超限的上隅角瓦斯運移情況進行數(shù)值模擬。

但是抽采條件下，采空區(qū)內(nèi)瓦斯運移復雜，目前使用的FLUENT軟件無法達到多物理場耦合模擬，導致相關(guān)抽采參數(shù)設置不合理[6]，本文在前人的基礎上，以保德煤礦81307工作面采空區(qū)抽采現(xiàn)狀為例，利用COMSOL 多物理場耦合模擬軟件建立采空區(qū)抽采三維模型，采用多孔介質(zhì)稀物質(zhì)流動和達西定律耦合，模擬研究布置間距、抽采負壓和抽采流量對治理上隅角瓦斯超限效果的影響，對抽采參數(shù)進行優(yōu)化，從而為采空區(qū)埋管抽采工作提供可靠的技術(shù)指導。

1 工作面概況

保德煤礦81307綜放工作面設計長2 489.5 m，工作面傾向長240 m，屬厚煤層，平均煤厚為7.3 m，設計工作面回采率為93%，工作面回采煤量為518.7萬t。工作面于2018年4月開始回采，通風方式為 “偏Y”型，即81307一號回風巷和81307膠帶巷進風、靠近切眼處聯(lián)巷回風，具體通風布置方式如圖1所示。

圖1 81307綜放工作面布置示意Fig.1 Schematic diagram of 81307 fully mechanized caving face

2 模擬方案的確定

2.1 模擬方案

為了研究抽采參數(shù)與上隅角瓦斯?jié)舛戎g的關(guān)系，確定最佳的埋管抽采參數(shù)。結(jié)合81307綜放工作面實際情況，采用COMSOL Multiphysics多物理場耦合軟件對以下4種情況進行數(shù)值模擬：①無埋管抽采條件下上隅角瓦斯分布規(guī)律；②埋管布置間距與上隅角瓦斯抽采效果間的規(guī)律；③抽采負壓與上隅角瓦斯抽采效果間的規(guī)律；④抽采流量與上隅角瓦斯抽采效果間的規(guī)律。分析結(jié)果最終得出最佳采空區(qū)埋管抽采參數(shù)。

2.2 幾何模型及網(wǎng)格劃分

由于采煤工作面存在各種設備，按實際尺寸建立幾何模型過程十分復雜，在設計模型時理想化參數(shù)條件，忽略采煤工作面各種設備的影響。將采煤工作面及其進回風巷都按矩形斷面處理。最終設計簡化模型偏Y型工作面的尺寸為：巷道寬5.2 m，高3.5 m，工作面長240 m，采空區(qū)長500 m，高30 m，同時對模型進行網(wǎng)格劃分并局部加密，如圖2～3所示。

圖2 幾何模型Fig.2 Geometric model

圖3 工作面及采空區(qū)模型網(wǎng)格劃分視Fig.3 Work area and goaf model meshing view

工作面模擬時以2個進風巷位置為入口邊界，忽略邊界湍流效應，設計入口處風流風速在進風巷斷面上大小一致，且入口邊界設定為法向速度場風流速度，入口初始風量2 200 m3/min。回風聯(lián)巷設置為出口邊界，埋管抽采采用低負壓抽采簡化為出口邊界。其余三維固體邊界設置為無通量邊界。開采層和鄰近層絕對瓦斯涌出量分別為7，10.15 m3/min；采空區(qū)及巷道內(nèi)初始壓力設為大氣壓；設定好模型和參數(shù)后，設定多孔介質(zhì)稀物質(zhì)流動和達西定律多物理場耦合求解，直至模型計算收斂為止。

3 模擬結(jié)果及分析

3.1 無埋管抽采條件下

設定采空區(qū)無埋管抽采，使用COMSOL軟件進行求解，最終得到的無抽采下采空區(qū)三維瓦斯分布云圖，如圖4所示。

圖4 無埋管條件下瓦斯分布Fig.4 Gas distribution diagram without buried pipe

無埋管抽采條件下，分析瓦斯?jié)舛仍茍D可以看出：采空區(qū)內(nèi)瓦斯?jié)舛入S著距離工作面的距離加大而不斷增加。鄰近工作面一側(cè)瓦斯?jié)舛容^低，隨著距離的不斷加大，瓦斯?jié)舛仍絹碓酱蟆６谧匀欢逊e區(qū)向壓實穩(wěn)定區(qū)的過渡區(qū)域，由于滲透率和孔隙率的改變，該區(qū)域瓦斯?jié)舛妊杆僭龃螅纬闪瞬煽諈^(qū)最高的瓦斯?jié)舛忍荻萚7-9]。

風流入口處瓦斯?jié)舛茸畹停拷ぷ髅嫖恢玫牟煽諈^(qū)的自然堆積區(qū)，由于碎脹性系數(shù)大和漏風量大導致其濃度較低。采空區(qū)無埋管抽采條件下，漏風流從工作面進風側(cè)漏向采空區(qū)，從回風側(cè)漏回到工作面，此時的漏風流會攜帶大量的瓦斯。在自然堆積區(qū)內(nèi)，沿工作面風流方向，風流隨自然堆積區(qū)的漏風擾動作用越來越小，從圖4中明顯看出自然堆積區(qū)內(nèi)靠近進風巷的地方被風流沖洗作用強，瓦斯?jié)舛鹊停拷仫L巷和上隅角處區(qū)域受風流擾動作用小，瓦斯大量積聚。此時，工作面上隅角瓦斯體積分數(shù)高達4.3%，回風巷瓦斯體積分數(shù)為2.2%，>1.0%(國家安全標準)，上隅角和回風巷瓦斯?jié)舛葒乐爻蓿仨氝M行瓦斯治理。

3.2 布置間距

研究抽采管布置間距與上隅角瓦斯?jié)舛鹊年P(guān)系，數(shù)值模擬布置間距為80，100，120 m時上隅角瓦斯?jié)舛确植家?guī)律。得到采空區(qū)三維瓦斯分布云圖如圖5～7所示。同時，將瓦斯分布云圖中數(shù)據(jù)進行量化，濃度數(shù)據(jù)如表1所示。

表1 不同布置間距下瓦斯體積分數(shù)Table 1 Gas concentration under different arrangement spacing

圖5 布置間距為80 m時的瓦斯分布Fig.5 Gas distribution when the arrangement spacing is 80 m

圖6 布置間距為100 m時的瓦斯分布Fig.6 Gas distribution when the arrangement spacing is 100 m

圖7 布置間距為120 m時的瓦斯分布Fig.7 Gas distribution when the arrangement spacing is 120 m

由圖5～7和表1可以發(fā)現(xiàn)，上隅角等地的瓦斯體積分數(shù)總體上隨著抽采管間的布置間距增大而增大，呈現(xiàn)正相關(guān)。當布置間距為120 m時，上隅角瓦斯體積分數(shù)為1.85%，回風流瓦斯體積分數(shù)為1.4%，上隅角和回風巷瓦斯體積分數(shù)嚴重超限。這是由于當布置間距較大時，抽采管離上隅角較遠，埋管抽采對采空區(qū)內(nèi)風流影響較小，此時抽采管很難及時的抽采上隅角處積聚的瓦斯[10]，治理效果差。當布置間距等于100 m時，上隅角瓦斯體積分數(shù)為0.50%，回風流瓦斯體積分數(shù)為0.35%，抽采管離上隅角距離近，能及時抽采上隅角處瓦斯，從而達到治理上隅角瓦斯體積分數(shù)超限目的。同時綜合考慮經(jīng)濟效益，埋管間距設置為100 m最佳。

3.3 抽采負壓

研究抽采負壓與上隅角瓦斯體積分數(shù)的關(guān)系，數(shù)值模擬抽采負壓為22，32，42 kPa時上隅角瓦斯體積分數(shù)分布規(guī)律，結(jié)果如圖8～10所示。

圖8 抽采負壓為22 kPa時的瓦斯分布Fig.8 Gas distribution when the pumping negative pressure is 22 kPa

圖9 抽采負壓為32 kPa時的瓦斯分布Fig.9 Gas distribution when the pumping negative pressure is 32 kPa

圖10 抽采負壓為42 kPa時的瓦斯分布Fig.10 Gas distribution when the pumping negative pressure is 42 kPa

分析抽采負壓為22，32，42 kPa的數(shù)值模擬結(jié)果，發(fā)現(xiàn)抽采負壓對上隅角瓦斯體積分數(shù)產(chǎn)生一定的影響，隨著抽采負壓的增大，上隅角瓦斯體積分數(shù)減小。將瓦斯分布云圖中數(shù)據(jù)進行量化，具體關(guān)系如表2所示。

表2 不同抽采負壓下瓦斯體積分數(shù)Table 2 Gas concentration under different suction negative pressure

通過表2可以看出，當抽采負壓為22 kPa時，上隅角瓦斯體積分數(shù)為3.5%，回風流瓦斯體積分數(shù)為2.8%。此時上隅角和回風巷瓦斯體積分數(shù)嚴重超限，雖然相比較無抽采條件下，采空區(qū)內(nèi)瓦斯流向采空區(qū)深部和抽采管處，但是抽采負壓較低，抽采瓦斯有限，沒有完全解決上隅角瓦斯超限問題。當抽采負壓為32 kPa時，上隅角體積分數(shù)為0.5%，達到國家規(guī)定的安全標準。從整體上來看，上隅角瓦斯體積分數(shù)隨著抽采負壓的增大而減小，呈現(xiàn)出減小的趨勢。但是隨著負壓的增大，有可能會造成采空區(qū)漏風嚴重，造成煤炭自燃[11-13]。同時考慮到礦井實際生產(chǎn)中加大負壓，經(jīng)濟投入加大，因此最佳抽采負壓可為32 kPa。

3.4 抽采流量

研究抽采流量與上隅角瓦斯體積分數(shù)的關(guān)系，數(shù)值模擬抽采流量為500，600，700 m3/min時上隅角瓦斯分布規(guī)律，結(jié)果如圖11～13所示。同時將瓦斯體積分數(shù)進行量化，得到數(shù)據(jù)如表3所示。

圖11 抽采流量為500 m3/min時的瓦斯分布Fig.11 Gas distribution when the pumping flow rate is 500 m3/min

圖12 抽采流量為600 m3/min時的瓦斯分布Fig.12 Gas distribution when the pumping flow rate is 600 m3/min

圖13 抽采流量為700 m3/min時的瓦斯分布Fig.13 Gas distribution when the pumping flow rate is 700 m3/min

由表3知，當抽采流量為500 m3/min時，上隅角瓦斯體積分數(shù)為2.1%，回風流面瓦斯體積分數(shù)為1.5%；當抽采流量為600 m3/min時，上隅角瓦斯體積分數(shù)為0.50%，回風流瓦斯體積分數(shù)為0.35%，達到國家安全生產(chǎn)標準；當抽采流量為700 m3/min時，上隅角瓦斯體積分數(shù)為0.34%，回風流瓦斯體積分數(shù)為0.2%。由此可見，上隅角等地的瓦斯體積分數(shù)隨著抽采流量的增大而減小，呈現(xiàn)負相關(guān)，抽采流量加大，低瓦斯體積分數(shù)區(qū)域不斷向采空區(qū)深部和上部擴散。抽采流量越大，上隅角瓦斯的治理效果越好，但是瓦斯抽采流量越大，越容易導致采空區(qū)漏風，采空區(qū)自燃“三帶”變大[14]，煤炭易自燃。綜合考慮保德煤礦現(xiàn)有抽采泵設備和經(jīng)濟效益等實際情況，抽采流量為600 m3/min時，治理效果最佳。

表3 不同抽采流量下瓦斯體積分數(shù)Table 3 Gas concentration under different pumping flow rates

4 現(xiàn)場試驗

4.1 施工方案

本次以保德煤礦81307工作面為試驗對象，采用采空區(qū)埋管抽采瓦斯的方法，治理上隅角瓦斯體積分數(shù)超限問題。試驗設置的采空區(qū)埋管抽采的參數(shù)為布置間距100 m、抽采負壓32 kPa、抽采流量600 m3/min。

試驗過程中由于要收集采空區(qū)和回采工作面上隅角瓦斯體積分數(shù)情況，因此在采空區(qū)采用提前鋪設束管。共采集2018年5～7月這3個月份的進風流瓦斯體積分數(shù)、工作面瓦斯體積分數(shù)、回風隅角流瓦斯體積分數(shù)和回風流瓦斯體積分數(shù)以此考察采空區(qū)埋管抽采治理上隅角瓦斯的效果。

4.2 效果考察

經(jīng)過監(jiān)測，整理出在設計的埋管抽采參數(shù)下81307工作面的5～7月份的瓦斯體積分數(shù)數(shù)據(jù)。為了便于研究埋管抽采對瓦斯體積分數(shù)的影響，對這些數(shù)據(jù)進行分析，采用曲線圖來分析瓦斯體積分數(shù)的變化趨勢。工作面、回風流及上隅角的瓦斯體積分數(shù)變化趨勢如圖14～16所示。

圖14 81037工作面瓦斯體積分數(shù)Fig.14 Gas concentration of 81037 working face

圖15 81037回風流瓦斯體積分數(shù)Fig.15 81037 return airflow gas concentration

圖16 81037上隅角瓦斯體積分數(shù)Fig.16 Gas concentration of the upper corner of 81037

通過圖14～16可以看出，采空區(qū)埋管抽采瓦斯后，81307工作面的瓦斯體積分數(shù)下降明顯。進風流中瓦斯體積分數(shù)試驗期間穩(wěn)定在0.1%左右，這是由于進風流不經(jīng)過采空區(qū)，埋管抽采影響較小。工作面瓦斯體積分數(shù)0.3%左右，回風流中瓦斯體積分數(shù)0.35%左右，上隅角瓦斯體積分數(shù)0.5%左右，小于1%，達到國家安全標準。采空區(qū)埋管抽采瓦斯，由于抽采管內(nèi)外壓差較大，導致采空區(qū)內(nèi)瓦斯體積分數(shù)降低，工作面漏風量下降，從而風流攜帶的瓦斯在上隅角處積聚的量也會下降，達到治理上隅角瓦斯體積分數(shù)超限問題的效果[15]。

5 結(jié)論

1)合理的采空區(qū)埋管抽采參數(shù)可以改變采空區(qū)內(nèi)漏風規(guī)律，減小采空區(qū)內(nèi)瓦斯含量，降低上隅角處瓦斯體積分數(shù)。模擬結(jié)果顯示81307工作面采空區(qū)埋管抽采布置間距為100 m、抽采負壓為32 kPa、抽采流量為600 m3/min時，安全和經(jīng)濟效益最佳。

2)現(xiàn)場抽采試驗，顯示采用模擬得到的抽采參數(shù)進行采空區(qū)瓦斯抽采，上隅角瓦斯體積分數(shù)降至0.5%左右，可以有效解決上隅角瓦斯超限問題。

3)對比分析數(shù)值模擬和現(xiàn)場試驗，利用COMSOL進行多物理場耦合模擬可以及時、準確的獲得采空區(qū)埋管抽采最佳參數(shù)，治理上隅角瓦斯體積分數(shù)超限問題，節(jié)省施工時間，保障生產(chǎn)安全，同時節(jié)約了資源，提高效益。