基于“U”型下行通風的綜放工作面瓦斯治理數值模擬研究

張增輝，薛彥平

（1.國家能源神東煤炭集團保德煤礦，山西忻州 036600；2.中煤科工集團沈陽研究院有限公司，遼寧撫順 113122；3.煤礦安全技術國家重點實驗室，遼寧撫順 113122）

瓦斯是威脅煤礦井下安全生產最重要因素之一[1]，綜放開采工藝在高產高效高瓦斯礦井應用過程中常因回風隅角瓦斯積聚而影響工作面安全生產，而采空區是主要瓦斯來源，占30%～60%[2]。采空區屬于塌陷冒頂區域，無法通過人員直接檢測瓦斯分布參數，目前只能通過束管采樣的方式測定采空區瓦斯氣體情況，但束管采樣存在測量數據誤差大，測量區域有限等諸多不利條件。隨著流體流場模擬技術的發展和日趨成熟，對采空區瓦斯分布進行數字模擬成為一種行之有效的技術手段[3]。保德煤礦綜放工作面由于產量集中，加之采空區瓦斯抽采量和抽采間距依據性不強，回采期間經常出現回風隅角和回風流瓦斯體積分數瀕臨超限的問題，嚴重影響安全高效生產。為此，通過研究綜放工作面采空區瓦斯運移規律，多角度數值模擬不同條件下埋管抽采對回風隅角瓦斯運移和濃度分布情況的影響，確定最佳采空區埋管抽采參數，為有效解決采空區瓦斯治理提供行之有效的方案和科學合理的依據。

1 工作面概況

81308 綜放工作面位于8 號煤二號主、輔運大巷以南，井田邊界以北，工作面以東為81307 綜放工作面采空區，以西為81309 備用工作面。工作面寬240 m，可采長度2 545 m，煤層平均厚度6.4 m，傾角平均4°，平均含夾矸3 層，單層夾矸最大厚度0.8 m，屬于較穩定煤層。開采煤層為自燃煤層，煤塵具有爆炸性。采用綜合機械化放頂煤采煤工藝，全部垮落法管理頂板，工作面設計采高3.7 m，放煤高度2.6 m，采放比為1∶0.7，設計日產量1.3 萬t/d，回采率93%。回采煤量為567 萬t，設計回采周期14 個月。工作面區域煤層原始瓦斯含量5.95 m3/t，采取了預抽煤層瓦斯的治理措施。工作面采前煤層殘余可解吸瓦斯含量最大3.78 m3/t，采前瓦斯抽采達標。通過分源預測法測算并參考鄰近工作面回采期間瓦斯涌出數據，工作面鄰近層和開采層通過采空區涌出的瓦斯量將占工作面瓦斯涌出總量的60%以上。

2 工作面采空區瓦斯運移模型

2.1 數學模型建立的基本條件

1）采空區氣體流動是個復雜的滲流過程，受制于支護條件、推進速度、采放比例、頂板垮落、采空區漏風等因素。總體上看，采空區內氣體可視為恒溫不可壓縮氣體，適用質量和動量守恒定律。

2）煤巖充填空間視為各向同性均質的含孔隙、裂隙的多孔介質，氣體流動服從達西定律，主要以擴散形式運移[4]。采空區內氣體為氧氣、氮氣和甲烷組成的混合氣體，不考慮裂隙水；遺煤瓦斯壓力較低，考慮為理想氣體，采空區氣體流動為等溫過程。

3）采空區矸石垮落和遺煤充填空間視為各向同性均質的含孔隙裂隙的多孔介質，遺煤和鄰近煤層連續釋放瓦斯。

基于上述假設和基本條件，構建符合連續性方程、動量方程和成分守恒方程的三維空間不可壓縮氣體流動的數值模型。

式中：ρ 為氣體密度，kg/m3；φ 為各方向上速度矢量或瓦斯體積分數求解變量的通用變量；t 為時間，s；u 為平均流速向量的分量，m/s；Γ 為廣義擴散系數；S 為廣義源項[5]。

2.2 采空區滲透率

采空區氣體流動，孔隙率和滲透率隨著垮落和充填物質參數發生變化，基于Kozeny-Carman 方程，構建采空區滲透率和孔隙率的動態變化方程。

式中：κ 為滲透率，μm2；Dp為介質的粒徑，m；n為介質的孔隙率，其大小取決于介質碎脹參數。

采空區多孔介質滲透率受到工作面推進長度、傾向長度、頂底板巖性及垮落、鄰近煤層及保安煤柱等相關參數影響，不同研究人員和不同計算方式確定的參數有很大差異[6]。

綜放工作面瓦斯涌出主要來自于煤壁、落煤和采空區[7]。依據鄰近已采81307 綜放工作面實際統計數據推測，81308 綜放工作面絕對瓦斯涌出量約為16.2 m3/min，采空區瓦斯涌出約10.9 m3/min。

2.3 邊界條件及物理模型

81308 綜放工作面通風方式為“U”型下行通風，即81308 1#進風巷和81308 運輸巷進風、81308 運輸巷靠近切眼處聯巷回風，81308 綜放工作面布置示意圖如圖1。

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

忽略其他因素的影響，僅考慮瓦斯抽采對采空區自燃危險區域分布的動態影響。參考81307 綜放工作面現場觀測數據，81308 綜放工作面采空區“兩道”浮煤厚度約4.3 m，中部采空區浮煤厚度約1.2 m，采空區深部距離工作面500 m 處按照壁面處理，采空區滲流范圍為工作面2 道之間的松散煤巖體以及跨落帶和斷裂帶。在高出煤層一定范圍后的漏風量很小，故假定影響滲流的裂隙帶在采空區范圍內各處高度一樣。因此，81308 工作面采空區物理模型范圍為：工作面長度240 m，采空區深部距工作面500 m；中部浮煤厚度為1.2 m，“兩道”厚度為4.3 m，將“兩道”煤柱和采空區深部500 m 設為漏風邊界。以采空區進風側外邦為坐標原點，向回風側方向為x 方向，垂直工作面到采空區深部方向為y 方向。

2.4 幾何模型及網格劃分

采煤工作面配備型號、尺寸各異的設備，難以完全按照實際條件建立幾何模型，在設計模型時理想化參數條件，忽略采煤工作面各種設備的影響。將采煤工作面及其進回風巷都按矩形斷面處理[8-9]。最終設計簡化模型尺寸為：巷道寬5.2 m，高3.5 m，工作面長240 m，采空區長500 m，高30 m，同時對模型進行網格劃分。

3 模擬結果

3.1 埋管布置間距對回風隅角瓦斯分布規律影響

采用Fluent 流體數值模擬軟件對工作面及回風隅角瓦斯體積分數分布規律進行數值模擬，不同埋管間距下采空區瓦斯體積分數云圖如圖2，不同埋管間距下工作面、回風隅角及回風流瓦斯體積分數如圖3。

圖2 不同埋管間距下采空區瓦斯體積分數Fig.2 Goaf gas concentration diagrams under different pipe spacings

圖3 不同埋管間距下工作面、回風隅角及回風流瓦斯體積分數Fig.3 Gas concentration of working face, return air corner and return air flow under different pipe spacings

通過圖2 和圖3 可知，當布置間距為130 m 時，埋管布置間距較大時，有效抽采半徑無法覆蓋回風隅角區域，治理效果差，回風隅角瓦斯體積分數為1.12%，回風流瓦斯體積分數為0.94%，回風隅角和回風流瓦斯體積分數瀕臨報警值。當布置間距小于等于120 m 時，有效抽采半徑覆蓋回風隅角區域，治理效果明顯，回風隅角瓦斯體積分數為0.80%，回風流瓦斯體積分數為0.70%。綜合考慮巷道條件和經濟效益等因素，建議埋管間距設置在90～120 m范圍內均比較合適。

3.2 抽采流量對回風隅角瓦斯分布規律影響

基于采空區埋管抽采流量分別為500、600、700 m3/min 條件下，對工作面、回風隅角及回風流瓦斯體積分數分布規律進行數值模擬，不同埋管抽采量下的采空區瓦斯體積分數云圖如圖4，不同瓦斯抽采流量下的工作面、回風隅角及回風流瓦斯體積分數如圖5。

圖4 不同埋管抽采量下采空區瓦斯體積分數Fig.4 Goaf gas concentration diagrams under different gas drainage amount of buried pipe

圖5 不同抽采流量下的工作面、回風隅角及回風流瓦斯體積分數Fig.5 Gas concentration of working face, return air corner and return air flow under different gas drainage flow

通過圖4 和圖5 可知，抽采流量越大，工作面及回風隅角瓦斯的治理效果越好。在抽采流量為500 m3/min 條件下，回風隅角瓦斯體積分數為1.1%，回風流瓦斯體積分數為1.04%；在抽采流量為600 m3/min條件下，回風隅角瓦斯體積分數為0.68%，回風流瓦斯體積分數為0.59%；在抽采流量為700 m3/min 條件下，回風隅角瓦斯體積分數為0.32%，回風流瓦斯體積分數為0.28%。考慮到瓦斯抽采流量越大，采空區漏風越嚴重，易造成采空區遺煤自燃[10-11]，建議抽采流量為600 m3/min，并采取一定的安全措施。

4 現場應用

保德煤礦81308 綜放工作面2019 年5 月開始回采，采用“U”型下行通風，配風量2 100 m3/min，沿81308 輔運巷設置1 趟DN800 mm 抽采主管，運輸巷和輔運巷之間通過聯巷或煤柱螺旋鉆孔設置DN800 mm 抽采支管。在間距大于100 m 的聯巷中間施工大孔徑螺旋鉆孔布置抽采支管，保證抽采支管間距在70～120 m 之間。為確保采空區瓦斯抽采安全，采取了自動監控和人工監控自然發火標志性氣體的安全措施，并在距最近在用抽采支管口100 m 范圍內安設自動噴粉抑爆裝置。

4.1 抽采參數與回風隅角瓦斯體積分數關系

1）抽采量是回風隅角瓦斯體積分數主要控制因素。抽采量在500 m3/min 左右時，回風隅角周平均瓦斯體積分數為0.8%，實測最大值超過1%；抽采量在610 m3/min 左右時，回風隅角周平均瓦斯體積分數在0.4%左右，實測最大值0.6%；抽采量在680 m3/min 左右時，回風隅角周平均瓦斯體積分數在0.24%左右，實測最大值0.42%。

2）埋管間距為回風隅角瓦斯體積分數的次要控制因素。在相同抽采量的前提下，回風隅角瓦斯體積分數隨埋管間距增大而增大，在間距70～120 m 范圍內，回風隅角瓦斯體積分數在±0.18%的范圍內波動。2019 年4 月至11 月期間，抽采參數及回風隅角瓦斯體積分數統計圖如圖6。

圖6 抽采參數及回風隅角瓦斯體積分數統計圖Fig.6 Statistics of drainage parameters and return air corner gas concentration

4.2 抽采參數與回風流瓦斯體積分數關系

1）抽采量是回風流瓦斯體積分數主要控制因素。抽采量在500 m3/min 左右時，回風流周平均瓦斯體積分數為0.6%，實測最大值超過0.92%；抽采量在610 m3/min 左右時，回風流周平均瓦斯體積分數為0.3%，實測最大值0.5%；抽采量在680 m3/min 左右時，回風流周平均瓦斯體積分數在0.18%左右，實測最大值0.36%。

2）埋管間距為回風隅角瓦斯體積分數的次要控制因素。在相同抽采量的前提下，回風隅角瓦斯體積分數隨埋管間距增大而增大，在間距70～120 m范圍內，回風流瓦斯體積分數在±0.11%的范圍內波動。2019 年4 月至11 月期間，抽采參數及回風流瓦斯體積分數統計圖如圖7。

圖7 抽采參數及回風流瓦斯體積分數統計圖Fig.7 Statistics of drainage parameters and return air flow gas concentration

81308 綜放工作面現場長達7 個月的試驗結果與數字模擬分析結果基本吻合，證明采取的數字建模方式是科學的、可靠的，具有超前指導礦井未來綜放工作面采空區埋管抽采參數的科學價值。

5 結 語

1）高瓦斯綜放工作面“U”型通風條件下，回風隅角瓦斯治理尤為關鍵。通過建立采空區瓦斯滲流數學模型，模擬不同抽采參數條件下工作面及采空區瓦斯運移及瓦斯體積分數分布狀況，可為合理選擇回風隅角瓦斯治理措施提供依據。

2）保德煤礦綜放工作面在“U”型下行通風條件下，依據模擬結果并通過現場試驗，采取合理的采空區埋管抽采參數有效解決了回風隅角和回風流瓦斯超限問題。在類似煤層瓦斯賦存及開采條件下，綜合考慮經濟和采空區防滅火等因素，埋管抽采間距為100 m，抽采量為580～610 m3/min。3）抽采量是影響工作面回風隅角及回風流瓦斯積聚的主控因素。隨著工作面瓦斯涌出量的增大，若繼續采取采空區埋管抽采的治理措施，勢必需相應增大抽采量，這樣將對采空區防滅火帶來不利影響。屆時，建議采取綜合瓦斯治理措施。