羅俊峰 蔣仲安 王洪勝 張中意

(北京科技大學土木與環(huán)境工程學院)

綜放采場的瓦斯來源主要有工作面、采空區(qū)、鄰近層瓦斯涌出等［1-2］，采空區(qū)瓦斯涌出占到采掘空間瓦斯涌出量約40%，甚至可以達到70%。在高瓦斯礦井中，上隅角瓦斯積聚和超限問題比較普遍，對煤礦安全生產(chǎn)造成了嚴重威脅，有可能引起瓦斯爆炸，造成巨大的財產(chǎn)損失和重大的人員傷亡。瓦斯治理主要采用的措施有:加強通風、加強瓦斯抽放以及綜合通風［3］。U+I是治理綜放采場采空區(qū)瓦斯的一種常用通風方式，通過對U+I型和U型兩種通風方式在不用風量情況下的采空區(qū)瓦斯運移規(guī)律的比較，為U+I型通風方式及其調(diào)控技術提供更充分的理論基礎。

1 采空區(qū)瓦斯流動方程

根據(jù)上覆巖層在采場中的移動情況，分析采空區(qū)頂板巖石性質(zhì)和垮落巖石的破壞情況。在橫向上可以將采空區(qū)劃分為自然堆積區(qū)、載荷影響區(qū)及壓實穩(wěn)定區(qū);在豎直方向上可以劃分為冒落帶、裂隙帶和彎曲下沉帶。其中，自然堆積區(qū)距工作面0～20 m，載荷影響區(qū)距工作面20～80 m，重新壓實區(qū)距工作面大于80 m。瓦斯在采空區(qū)中的流速極低，盡管壓力梯度很大，仍可以將瓦斯在采空區(qū)中的流動視為不可壓流動［4-5］。因此，可將采空區(qū)瓦斯流動視為近似層流運動，這也是大量研究基于Darcy定律的原因。

1.1 流體力學控制方程

CFD的基礎是建立Navier-Stocks方程。流體運動遵循最基本的3個守恒定律，即質(zhì)量守恒定律、動量守恒定律和能量守恒定律。具體在流體力學中，就是流體力學中的連續(xù)性方程、動量守恒方程和能量方程［6-8］。

為了研究的方便，將采空區(qū)視作各向同性的多孔介質(zhì)。采空區(qū)多孔介質(zhì)的動量損失特性可以用下式描述:

式中，Si為第i個(x，y或z)動量方程的源;D、C為預定義矩陣;μ為分子黏度為速度絕對值;vj為在x，y和z方向上的速度分量。

式(1)中第一項為黏性損失項，即Darcy定律，是壓力梯度驅(qū)動下的黏性損失;第二項為慣性損失項，是高速流動中的動量損失，對于采空區(qū)多孔介質(zhì)的低速流動，第二項可以忽略不計。因此可用下式描述采空區(qū)多孔介質(zhì)的動量損失源項:

式中，α為多孔介質(zhì)的滲透率，m2。

1.2 k-ε 方程

瓦斯在采空區(qū)中的流動視為層流，在巷道中的流動雷諾數(shù)Re較大，并且在巷道中存在拐角，應視作湍流。k-ε模型為FLUENT中自帶的可處理氣體流動湍流狀態(tài)的流動模型，源于嚴格的數(shù)學統(tǒng)計技術，其湍動能k和耗散率ε方程可以用下式表示:

式中，ρ為流體密度;k為湍動能;ε為湍動能耗散率;μ為動力黏度系數(shù);μt為湍流黏度;σk為流體流動應力;Gk為平均速度梯度引起的湍流動能的產(chǎn)生項;Gb為由于浮力引起的湍流動能的產(chǎn)生項;YM為可壓湍流中脈動擴張的貢獻;σε為導致應變的應力;C1ε、C2ε、C3ε為模型常數(shù)。

2 工程背景和數(shù)值模型

2.1 工作面概況

漳村煤礦位于山西省長治市候堡鎮(zhèn)境內(nèi)，該礦2306工作面開采煤層屬二疊系下統(tǒng)山西組的3#煤層，位于23采區(qū)的南部，東面位于23采區(qū)大巷，北為2305工作面采空區(qū)，西距常村煤礦井田邊界188 m，南距井田邊界20～24 m。2306工作面煤層厚度5.7 ～6.9 m，煤層賦存穩(wěn)定，煤厚變異較小，結(jié)構(gòu)較簡單，兩層夾矸穩(wěn)定發(fā)育。煤層直接頂板為0～7 m泥巖，堅固性系數(shù)f約2.5，直接底板為0～6.4 m 泥巖;煤層傾角 2°～6°;工業(yè)儲量 207 萬 t，可采儲量為192萬t;煤層瓦斯涌出量較大，回采瓦斯涌出量為10～15 m3/min。

2.2 數(shù)值模型

2.2.1 幾何模型

2306工作面傾斜長度167.7 m，工作面可采長度1 443 m，面積 2 082 249 m2，煤層厚 6.34 m。進風巷為矩形截面，寬5 m，高3.3 m。回風巷為矩形截面，寬 5.0 m，高 3.5 m。瓦斯巷矩形截面寬4.4 m，高3.5 m，沿煤層頂板1 m 布置。采空區(qū)寬177.7 m，長取200 m，高度由經(jīng)驗公式計算。

根據(jù)現(xiàn)場實際情況和數(shù)值模擬的需要，將2306工作面采空區(qū)幾何模型進行了一定程度的簡化:①工作面、進風巷、回風巷、瓦斯巷和采空區(qū)均視作長方體，不考慮工作面設備情況;②將采空區(qū)分為冒落帶1、冒落帶2、冒落帶3，裂隙帶1和裂隙帶2，以區(qū)別不同的孔隙率;③采空區(qū)走向長度取200 m，傾斜長度取178 m(包括進、回風巷長度)，垂直高度取40 m;④進風巷、回風巷長度均取15 m;⑤進風巷寬度5 m，高度3 m，回風巷寬度5 m，高度3 m，瓦斯巷寬度4.4 m，高度3.5 m，工作面長度取168 m，寬度取5 m，高度取3 m。

2.2.2 基本假設

為簡化分析，對采空區(qū)以及采空區(qū)氣體進行如下假設:①本煤層和煤層上覆巖層是均勻連續(xù)的;②因采空區(qū)內(nèi)漏風風速較小，而采空區(qū)兩側(cè)壓差也較小，近似認為采空區(qū)氣體不可壓縮;③采空區(qū)溫度恒定假設，空氣進入采空區(qū)要和采空區(qū)內(nèi)的煤體發(fā)生氧化反應，各處的溫度并非恒定一致，為方便研究，不考慮化學反應和傳熱效應;④采空區(qū)內(nèi)滲透系數(shù)各向相同性，因為采空區(qū)內(nèi)冒落巖石堆積的無序性，無法區(qū)分哪個方向的孔隙率較大，故將采空區(qū)各處不同方向上的滲透系數(shù)考慮為一致;⑤工作面進風不考慮水蒸氣等小組分氣體，假設空氣中只含有氧氣、氮氣，無瓦斯氣體。

2.2.3 邊界條件設置

采空區(qū)滲流選取 Realizable k-模型［9］;在采空區(qū)的冒落帶和裂隙帶分別設置瓦斯源，模擬采空區(qū)瓦斯涌出;進風口設置為速度進口，即velocity-inlet;回風口設置為壓力出口，即pressure-outlet;選擇組分運輸模型，只考慮3種主要氣體，即 N2，O2，CH4;其余設置為墻［10］。

根據(jù)實際情況和前人的研究結(jié)果，對采空區(qū)巖體的孔隙率、黏性阻力系數(shù)和瓦斯源項按照表1設置［4，6，11］。

3 數(shù)值模擬研究結(jié)果

按照幾何模型建模劃分網(wǎng)格，根據(jù)表1的邊界條件設置模型參數(shù)。

表1 邊界條件的設置

3.1 不同風量下的數(shù)值模擬

在進風巷道風速為 1.0，1.5，2.0，2.5 和 3.0 m/s時，對U型通風方式下的采空區(qū)瓦斯運移規(guī)律進行數(shù)值模擬。在不同風量下，取z=0～3 m(x=199 m，y=178 m)作瓦斯?jié)舛茸兓€，見圖1。

圖1 U型上隅角瓦斯?jié)舛茸兓?/p>

分析圖1可知:

(1)工作面上隅角瓦斯?jié)舛仁芄ぷ髅骘L度的影響產(chǎn)生變化:在工作面風速較小時(小于2.0 m/s)，瓦斯?jié)舛茸兓^小;工作面風速較大時(大于2.0 m/s)，瓦斯?jié)舛茸兓^大，相比較低風速，瓦斯?jié)舛认陆盗?0%左右。

(2)U型通風方式下，工作面上隅角的瓦斯?jié)舛却笥?%，在高瓦斯?jié)舛鹊墓ぷ髅妫琔型通風不能滿足治理上隅角瓦斯超限的要求［12］。

3.2 U+I型與U型通風方式數(shù)值模擬對比

U+I通風方式的幾何模型如圖2所示，幾何模型尺寸如前所述，瓦斯巷和回風巷之間的水平距離為15 m，豎直距離6 m。根據(jù)表1的邊界條件設置模型參數(shù)，U型通風方式的幾何模型同U+I通風方式相比，沒有瓦斯巷。

圖2 采空區(qū)U+I幾何模型

比較兩種通風方式在工作面風速1.5 m/s時采空區(qū)冒落帶的瓦斯?jié)舛确植?，?shù)值模擬見圖3。

分析圖3，比較不同通風方式下采空區(qū)瓦斯?jié)舛仍茍D可知:

圖3 采空區(qū)瓦斯?jié)舛仍茍D

(1)自然堆積區(qū)內(nèi)，U+I型通風方式下的風速和漏風量相比于U型通風方式要大，該區(qū)域內(nèi)的瓦斯在U+I通風方式下被稀釋的程度較大，瓦斯?jié)舛容^小，一般低于1%。

(2)載荷影響區(qū)內(nèi)孔隙率較小，黏性阻力系數(shù)較大，漏風量相對較小。在U型通風方式下，瓦斯?jié)舛葹?0% ～60%;在U+I通風方式下，濃度為20%～50%。載荷影響區(qū)內(nèi)的瓦斯是造成采空區(qū)回風巷瓦斯?jié)舛容^高的原因之一。

(3)壓實穩(wěn)定區(qū)內(nèi)的瓦斯孔隙率很小，黏性阻力系數(shù)很大，通風對該區(qū)域的影響較小，瓦斯運移的動力主要是瓦斯?jié)舛葔毫?，此區(qū)域內(nèi)的瓦斯?jié)舛容^高，一般在70%以上。

(4)在進行采空區(qū)瓦斯抽放時，抽放孔應布置在瓦斯?jié)舛认鄬^高的區(qū)域內(nèi)，即回風側(cè)的壓實穩(wěn)定區(qū)內(nèi)或者載荷影響區(qū)靠近壓實穩(wěn)定區(qū)的范圍內(nèi)。

3.3 上隅角瓦斯?jié)舛鹊臄?shù)值模擬對比

采煤工作面上隅角處風流拐彎，易形成渦流區(qū)，瓦斯難以被風流帶走，容易形成瓦斯積聚超限。降低上隅角瓦斯?jié)舛鹊拇胧┲饕峭ㄟ^改變工作面通風方式、瓦斯抽放和稀釋上隅角瓦斯?jié)舛?加裝局扇、導風簾等)等措施。

U+I和U型通風方式下回風側(cè)上隅角瓦斯?jié)舛茸兓鐖D4所示。其中U+I(高位)通風方式是指相對于U+I型的通風方式，在豎直高度上提高瓦斯巷5 m，這是為了避免由于頂板變形或者底板塌透，導致瓦斯巷失效，發(fā)揮不了治理瓦斯的作用。

分析圖5，對比兩種不同通風方式下上隅角瓦斯?jié)舛瓤芍?

圖4 上隅角瓦斯?jié)舛葘Ρ?/p>

(1)U型通風方式下，上隅角瓦斯?jié)舛入S著高度的增加而逐漸增加，由于瓦斯密度小于空氣密度，導致上層瓦斯?jié)舛却笥谙聦油咚節(jié)舛取?/p>

(2)U+I和U+I(高位)通風方式下，上隅角瓦斯?jié)舛认啾扔赨型通風方式大大降低，冒落帶高度在12 m以下時，瓦斯?jié)舛鹊陀?.8%。

(3)對比U+I和U+I(高位)通風方式可以看到，U+I(高位)通風方式對于冒落帶12 m以上的區(qū)域瓦斯抽放效果要好于U+I型通風方式，瓦斯?jié)舛认鄬^低。因此采用瓦斯巷治理瓦斯時，建議布置高位瓦斯巷，有助于提高瓦斯抽放效果。

4 結(jié)論

(1)采空區(qū)3個分區(qū)受到工作面通風的影響具有差異，自然堆積區(qū)的孔隙率較大，黏性阻力系數(shù)較小，靠近工作面流速較大，瓦斯?jié)舛群苄?載荷影響區(qū)的流速次之，自載荷影響區(qū)至壓實穩(wěn)定區(qū)的瓦斯?jié)舛瘸试黾于厔?而壓實穩(wěn)定區(qū)的瓦斯在數(shù)值模擬收斂后基本無瓦斯運移發(fā)生，受工作面通風影響有限，瓦斯?jié)舛群芨摺?/p>

(2)U型通風方式下，改變工作面風速對上隅角瓦斯?jié)舛扔绊懹邢蓿诟咄咚共擅汗ぷ髅骓毑捎闷渌L方式，防止瓦斯?jié)舛瘸蕖?/p>

(3)U+I和U+I(高位)型通風方式對治理上隅角瓦斯?jié)舛瘸蘧哂忻黠@的效果，有效降低上隅角瓦斯?jié)舛戎?.8%以下。在這兩種通風方式下，采空區(qū)瓦斯隨漏風經(jīng)瓦斯巷排出，大大降低了采空區(qū)自然堆積區(qū)的瓦斯?jié)舛?。但是由于大部分瓦斯?jīng)由瓦斯巷排出，要注意瓦斯巷的瓦斯?jié)舛瘸迒栴}。

(4)比較U+I和U+I(高位)兩種通風方式，在條件允許的情況下，采空區(qū)U+I(高位)通風方式，對于降低采空區(qū)上部瓦斯?jié)舛染哂懈玫男Ч?/p>

(5)U+I和U+I(高位)適用于不易自燃高瓦斯采煤工作面。針對不同的工作面條件，具體分析采空區(qū)瓦斯運移變化，運用以上兩種通風方式，提高瓦斯抽放效果。

［1］周世寧，林柏泉.煤層瓦斯流動賦存理論［M］.北京:煤炭工業(yè)出版社，1998.

［2］胡千庭.煤礦瓦斯抽采與瓦斯災害防治［M］.徐州:中國礦業(yè)大學出版社，2007.

［3］俞啟香，王 凱，楊勝強.中國采煤工作面瓦斯涌出規(guī)律及其控制研究［J］.中國礦業(yè)大學學報，2000(1):9-14.

［4］李宗翔.綜放工作面采空區(qū)瓦斯涌出規(guī)律的數(shù)值模擬研究［J］.煤炭學報，2002(2):173-178.

［5］胡千庭，梁運培，劉見中.采空區(qū)瓦斯流動規(guī)律的CFD模擬［J］.煤炭學報，2007(7):719-723.

［6］王福軍.計算流體動力學分析［M］.北京:清華大學出版社，2004.

［7］高 魁，劉澤功，劉 健，等.工作面漏風對采空區(qū)瓦斯流動規(guī)律影響的數(shù)值模擬［J］.煤礦安全，2012(7):8-11.

［8］胡方坤，時國慶，王德明，等.高瓦斯采空區(qū)多組分氣體場分布數(shù)值模擬［J］.煤礦安全，2013(12):14-16.

［9］王 凱，蔣曙光，張衛(wèi)清，等.尾巷改變采空區(qū)瓦斯流場的數(shù)值模擬研究［J］.采礦與安全工程學報，2012(1):124-130.

［10］金龍哲，姚 偉，張 君.采空區(qū)瓦斯?jié)B流規(guī)律的CFD模擬［J］.煤炭學報，2010，35(9):1476-1480.

［11］高建良，王海生.采空區(qū)滲透率分布對流場的影響［J］.中國安全科學學報，2010(3):81-85.

［12］李宗翔，孫廣義，秦書玉.回采采空區(qū)上隅角瓦斯治理的數(shù)值模擬與參數(shù)確定［J］.中國地質(zhì)災害與防治學報，2001(4):9-12.