E0103工作面CO₂預(yù)裂技術(shù)數(shù)值模擬研究

唐強

摘 要：本文研究礦井為煤與瓦斯突出煤礦，E0103工作面存在煤與瓦斯突出危險性。為提高該工作面煤與瓦斯抽采效率，該煤礦在巷道掘進過程中運用CO2預(yù)裂技術(shù)進行增透，然后進行抽采工作。本文對CO2釋放過程中煤巖滲透率分布情況進行數(shù)值模擬研究，建立含瓦斯煤流-固耦合模型及CO2與瓦斯解析置換方程，運用軟件COMSOL并結(jié)合MATLAB對預(yù)裂技術(shù)實施過程進行數(shù)值模擬，得出工作面使用預(yù)裂技術(shù)過程中不同時間段煤巖滲透率分布云圖。結(jié)果表明：模擬結(jié)果同現(xiàn)場深孔快速取樣的樣本分析結(jié)果基本一致。

關(guān)鍵詞：瓦斯突出;CO2預(yù)裂;流-固耦合;滲透率

中圖分類號：TD712.6文獻標識號：A文章編號：1003-5168（2021）11-0073-03

The Numerical Simulation Study of CO2 Pre-Splitting Technique

on E0103 Working Face

TANG Qiang

（China Coal Technology Engineering Group Chongqing Research Institute， Chongqing 400037）

Abstract： This paper research mine was a coal and gas outburst mine， There is danger of coal and gas outburst in E0103 working face. The technology of CO2 pre-splitting antireflection was used in the process of the mine roadway drivage to improve the working face of coal and gas extraction efficiency， then begin extraction. This paper performed numerical simulations mainly for the CO2 release process of coal rock permeability distribution. We established a model of Flow - solid coupling model and CO2 and gas analytical displacement equation， using the software COMSOL combined with MATLAB to pre-splitting technology implementation process. We could get coal rock permeability distribution in the process of different period for working face with pre-splitting technique. Simulation results was basically identical with the field of deep hole fast sampling sample analysis results.

Keywords： gas outburst;CO2 pre-splitting;flow-solid coupling;permeability

近年來，煤礦形勢逐漸嚴峻，瓦斯突出事故頻發(fā)[1]。瓦斯抽采是防治瓦斯事故的關(guān)鍵措施，但是，部分礦井開采的煤層透氣性較差，導(dǎo)致瓦斯抽采效果不佳，瓦斯治理時間長[2]。為此，眾多學(xué)者對低透性煤層增透技術(shù)展開研究，主要技術(shù)措施包括深孔松動爆破、水力壓裂、水力割縫、CO2增透等，不同的增透技術(shù)有其自身的特點及應(yīng)用范圍[3-4]。本次研究礦井為煤與瓦斯突出礦井，為節(jié)約成本，同時消除該礦E0103工作面煤與瓦斯突出危險，礦井在掘進過程中運用CO2增透技術(shù)增加煤層透氣性，以達到快速抽采達標的目的。本文以此為背景[5-7]，建立該煤礦含瓦斯煤流-固耦合模型，模擬出運用CO2增透技術(shù)及未使用CO2增透技術(shù)兩種情況下，不同時間段里煤層中瓦斯壓力分布云圖，并將模擬結(jié)果同深孔快速取樣分析結(jié)果進行對比，發(fā)現(xiàn)二者的結(jié)果基本一致，故可推斷本文所建模型合理。

1 煤礦概況

研究礦井地處云南，為煤與瓦斯突出礦井，E0103工作面可采長度約700 m，工作面長約130 m，煤厚2.52 m，傾角15°，夾矸層數(shù)0～3層，采用綜合機械化開采工藝采煤;5煤直接頂為泥質(zhì)灰?guī)r、裂隙不發(fā)育，上部為粉砂質(zhì)泥巖和粉砂巖，直接底板為泥巖，底板下大概有三層煤線，煤線間夾泥巖。根據(jù)《地勘報告》，5煤瓦斯含量取7.9～15.45 m3/t，煤層硬度系數(shù)[f]為0.37，煤層透氣性系數(shù)為0.607～2.67 m2/（MPa2·d）。

2 CO2預(yù)裂過程數(shù)值模擬研究

鉆孔長度約為100 m，在鉆孔底部下放CO2爆破器主管，進行啟爆。CO2從泄能頭兩側(cè)的出氣孔急速沖出，時間約為40 ms，本文主要模擬這一時間段內(nèi)CO2對煤巖的影響。實驗研究表明，CO2發(fā)爆器啟動之后，會在鉆孔內(nèi)產(chǎn)生一個約20 MPa的瞬間壓力，該壓力會對含瓦斯煤層產(chǎn)生一種瞬態(tài)的能量，同時CO2會置換出含瓦斯煤中的瓦斯氣體。

2.1 模擬方程

通過現(xiàn)場實際考察可知，進行CO2預(yù)裂過程中，溫度對瓦斯解析的影響較小，故本文所建立的數(shù)學(xué)模型不考慮溫度對預(yù)裂過程的影響。在CO2釋放過程中，壓力逐漸衰減，而以往的模擬忽略了這一點，根據(jù)實驗研究得出壓力衰減方程[5]。泄能頭兩側(cè)的出氣孔噴出CO2壓力方程為：

[P=P0-tB] ? ? ? ? ? ? ? ? ? ? ? ? ? ? （1）

式中：[P]表示泄能頭兩側(cè)的出氣孔噴出CO2壓力;[P0]表示泄能頭兩側(cè)出氣孔噴出CO2初始壓力;[B]表示壓力衰減系數(shù);[t]表示時間。

煤層內(nèi)的瓦斯含量動態(tài)方程為[6-7]：

[MS=abcpρS1+bp] ? ? ? ? ? ? ? ? ? ? ? ? （2）

式中：[MS]表示瓦斯含量;[a]表示吸附常數(shù);[b]表示吸附常數(shù);[c]表示煤質(zhì)校正參數(shù)，計算公式如式（3）所示;[p]表示瓦斯壓力;[ρS]表示煤的密度。

[c=ρs11+0.31M×100-A-M100] ? ? ? ? ? ? ? ? ? （3）

式中：[M]表示煤的水分;[A]表示煤的灰分。

瓦斯運動方程[8-11]：

[V=-kμ1+mp?p] ? ? ? ? ? ? ? ? ? ? ? ?（4）

式中：[V]表示瓦斯?jié)B流速度;[k]表示滲透率;[μ]表示瓦斯黏度系數(shù);[m]表示Klinkerberg系數(shù)。

[?Ms?t+??ρV=0] ? ? ? ? ? ? ? ? ? ? （5）

式中：[V]表示瓦斯?jié)B流速度;[ρ]表示瓦斯密度。

2.2 模擬參數(shù)

將現(xiàn)場采集的煤樣帶到實驗室進行實驗參數(shù)測定，得到煤巖基本參數(shù)，具體見表1。

2.3 數(shù)值模擬

現(xiàn)有的基本設(shè)計圖如圖1所示。數(shù)值模擬所建立的幾何模型如圖2所示。

建立的幾何模型中，鉆孔長度為90 m，將建立的基本數(shù)學(xué)模型編譯成MATLAB（矩陣實驗室軟件）文件，鏈接到耦合模擬軟件COMSOL Multiphysics（康模數(shù)爾軟件）中，模擬時間為40 ms，步長為5 ms，現(xiàn)場進行測試得到頂板壓力為2.3 MPa，底板壓力為2 MPa，軸向壓力為0.8 MPa。根據(jù)實驗室測量的煤巖樣本基本參數(shù)級進行數(shù)值模擬，分別得出預(yù)裂之后10 ms、20 ms、30 ms和40 ms的滲透率如圖3至圖6所示。

由圖3至圖6可知，進行CO2預(yù)裂技術(shù)之后，鉆孔孔底滲透率逐漸增大，隨著時間的推移，滲透率增大，影響的范圍也逐漸增加。滲透率變化范圍較大的區(qū)域為鉆孔之間耦合區(qū)域，而鉆孔非耦合區(qū)域僅僅依靠孔底的CO2預(yù)裂仍舊無法到達致裂要求。

2.4 模擬結(jié)果分析

由模擬結(jié)果可知：進行CO2預(yù)裂之后，10 ms時滲透率的影響范圍約為5 m;20 ms時滲透率的影響范圍約為30 m;30 ms時滲透率的影響范圍約為40 m;40 ms時滲透率的影響范圍約為60 m;隨著時間的推移，滲透率的影響范圍會不斷增大，滲透率提高會提升瓦斯抽采效率。可見，影響煤層滲透率較大的區(qū)域為鉆孔耦合區(qū)域，而雖然非耦合區(qū)域的煤層滲透率發(fā)生了變化，但變化并不明顯，不能滿足預(yù)裂要求。通過深孔快速取樣對耦合區(qū)域和非耦合區(qū)域的煤樣進行滲透率參數(shù)檢測，檢測結(jié)果與數(shù)值模擬結(jié)果基本一致。因此，在實際施工過程中，需要對原方案進行優(yōu)化，不僅要在孔底放置CO2預(yù)裂爆破裝置，還要在鉆孔中間段一定位置放置CO2預(yù)裂爆破裝置;同時，對于滲透率變化不大的非耦合區(qū)域，仍需要進行CO2預(yù)裂，消除低滲透率盲區(qū)。

3 結(jié)語

本文針對礦井高瓦斯低透氣性煤層的特點，建立含瓦斯煤層流-固耦合模型。針對現(xiàn)場工程實踐進行數(shù)值模擬，得出進行CO2預(yù)裂技術(shù)后不同時間段的煤層滲透率分布云圖。同時，進行深孔快速定點取樣，將樣本帶到實驗室中進行測定，測得的滲透率與數(shù)值模擬對比，二者結(jié)果基本一致。研究結(jié)果對推進山西省高瓦斯低透氣性煤層CO2預(yù)裂技術(shù)具有一定的指導(dǎo)意義和現(xiàn)實意義。在實際施工中，要不斷優(yōu)化CO2預(yù)裂技術(shù)，以滿足低成本消除礦井煤與瓦斯突出危險的要求。

參考文獻：

[1]趙丹，吳禹默，王明玉，等.低滲透煤層中復(fù)合射孔技術(shù)的應(yīng)用研究[J].能源技術(shù)與管理，2021（1）：4-7.

[2]王正帥.水力化增透技術(shù)在碎軟煤層瓦斯抽采中的應(yīng)用研究[J].煤炭工程，2021（2）：85-89.

[3]李喜員，榮海，范超軍.低透氣性煤層水力沖孔和抽采鉆孔瓦斯治理模擬研究[J].礦業(yè)安全與環(huán)保，2021（1）：11-16.

[4]張曉剛，姜文忠，都鋒.高瓦斯低透氣性煤層增透技術(shù)發(fā)展現(xiàn)狀及前景展望[J].煤礦安全，2021（2）：169-176.

[5]張悅，張民波，朱天玲，等.低透氣性煤層CO2增透預(yù)裂技術(shù)應(yīng)用[J].科技導(dǎo)報，2013（23）：36-39.

[6]王惠明.CO2預(yù)裂增透瓦斯抽采技術(shù)的應(yīng)用：以宏遠煤業(yè)150103運輸順槽掘進工作面為例[J].技術(shù)與市場，2014（12）：53-55.

[7]楊曉國.CO2預(yù)裂增透技術(shù)在高瓦斯低透氣性煤層中的應(yīng)用[J].中國煤炭，2015（4）：122-126.

[8]吳錦旗.液態(tài)CO2預(yù)裂強化預(yù)抽消突技術(shù)在突出煤層揭煤過程中的應(yīng)用[J].煤炭與化工，2015（7）：105-106.

[9]李寶華.CO2增透預(yù)裂技術(shù)在常村煤礦瓦斯抽采中的應(yīng)用[J].煤，2015（9）：53-54.

[10]郭有慧，孫銳.寺家莊礦CO2預(yù)裂爆破強化抽采技術(shù)試驗研究[J].科技與創(chuàng)新，2014（1）：6-8.

[11]陳運.基于CO2預(yù)裂爆破的斜井快速揭煤技術(shù)[J].能源技術(shù)與管理，2016（1）：34-35，43.