高 海,程傳興

(1.中煤科工集團(tuán)沈陽研究院有限公司 煤礦安全技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,遼寧 撫順 113122;2.山東科技大學(xué) 安全與環(huán)境工程學(xué)院,山東 青島 266590)

煤炭作為世界三大能源之首,為國民經(jīng)濟(jì)發(fā)展做出了突出貢獻(xiàn)[1]。近年來,我國鼓勵(lì)新能源發(fā)展,但煤炭仍是我國能源消費(fèi)的重要組成部分[2-3]。然而,隨著機(jī)械化水平不斷提高,煤巷掘進(jìn)速度得到提升,同時(shí)產(chǎn)生的粉塵量也大大增加[4]。其中,粉塵中含有的大量SiO2是導(dǎo)致塵肺病的重要原因之一。塵肺病是我國目前最嚴(yán)重的職業(yè)病,2020年全國共報(bào)告各類職業(yè)病新病例中職業(yè)性塵肺病14 367例,占84.2%。因此改善綜掘工作面作業(yè)環(huán)境、保障作業(yè)工人職業(yè)衛(wèi)生健康尤為重要。

20世紀(jì)初,美國、日本等國外學(xué)者對(duì)粉塵在巷道中的運(yùn)動(dòng)進(jìn)行了相關(guān)研究,并得出了與粉塵運(yùn)移相關(guān)的數(shù)學(xué)模型[5],從此揭開了對(duì)粉塵運(yùn)移研究的新篇章。MORENO等[6]研究發(fā)現(xiàn)巷道通風(fēng)可以有效降低可吸入粉塵的質(zhì)量和濃度。HAMDANI等[7]驗(yàn)證了粉塵沉積和表面粗糙度的關(guān)系,進(jìn)一步驗(yàn)證了粉塵的沉積規(guī)律。WITT等[8]利用計(jì)算機(jī)流體動(dòng)力學(xué)(CFD)模型預(yù)測了傳輸裝置附近揚(yáng)塵的運(yùn)動(dòng)軌跡。TORANO等[9]根據(jù)現(xiàn)場測量結(jié)果驗(yàn)證了基于CFD的作業(yè)過程中粉塵擴(kuò)散行為模擬的可靠性,并發(fā)現(xiàn)改變風(fēng)筒與底板和工作面之間的距離會(huì)極大地影響粉塵控制性能。COLINET等[10]詳細(xì)介紹了粉塵控制領(lǐng)域的研究成果,匯總了各種粉塵控制的方法,為礦井粉塵防治工作提供了有力幫助。

除此之外,國內(nèi)學(xué)者也對(duì)掘進(jìn)工作面風(fēng)流及粉塵的運(yùn)動(dòng)進(jìn)行了詳細(xì)研究。陳榮策等[11]在現(xiàn)場實(shí)驗(yàn)的基礎(chǔ)上,研究了沒有局部通風(fēng)的條件下,掘進(jìn)鑿巖時(shí)粉塵的產(chǎn)生及分布規(guī)律。李雨成等[12]分析了掘進(jìn)工作面采用壓入式、抽出式及長壓短抽式3種通風(fēng)方式下,風(fēng)筒距工作面不同距離時(shí)的粉塵分布特征。蔣仲安等[13]研究了掘進(jìn)工作面采取長壓短抽通風(fēng)方式下,風(fēng)筒高度與直徑對(duì)渦流的作用,并對(duì)長壓短抽通風(fēng)系統(tǒng)參數(shù)進(jìn)行了優(yōu)化。秦躍平等[14]利用FLUENT模擬軟件研究了“長壓短抽”通風(fēng)方式下,掘進(jìn)巷道內(nèi)粉塵的運(yùn)移和分布規(guī)律,提出了“壓風(fēng)分流”通風(fēng)方式,并對(duì)其降塵效果進(jìn)行了模擬。程衛(wèi)民等[15]利用FLUENT軟件對(duì)綜掘工作面旋流氣幕抽吸控塵進(jìn)行了數(shù)值模擬,并設(shè)計(jì)了由附壁風(fēng)筒和抽塵凈化裝置構(gòu)成的綜掘工作面旋流氣幕抽吸控塵系統(tǒng)。

國內(nèi)外學(xué)者都對(duì)巷道內(nèi)粉塵運(yùn)移規(guī)律進(jìn)行了詳細(xì)的研究,但現(xiàn)有研究對(duì)單壓入式通風(fēng)條件下掘進(jìn)巷道風(fēng)流-粉塵運(yùn)移規(guī)律及最佳通風(fēng)控塵風(fēng)量的研究并不充分。基于此,本文對(duì)壓入式掘進(jìn)巷道風(fēng)流及粉塵的運(yùn)移規(guī)律及最佳控塵風(fēng)量進(jìn)行了研究,進(jìn)而客觀地預(yù)測了通風(fēng)除塵的效果,為解決巷道粉塵嚴(yán)重污染問題提供理論支持。

1 工作面概況

3206回風(fēng)順槽工作面地質(zhì)條件簡單,褶曲寬緩,煤層賦存平穩(wěn),約在1 400 m處有一處地質(zhì)構(gòu)造,其他隱伏地質(zhì)構(gòu)造需在掘進(jìn)過程中進(jìn)一步探明。3#煤層位于二疊系下統(tǒng)山西組下部,所采煤層平均埋藏深度約為500 m;煤層平均厚度為5.76 m,黑色,亮煤為主,煤層穩(wěn)定,全區(qū)可采,煤層傾角為2°～10°,平均為6°;煤層單軸抗壓強(qiáng)度平均為18.09 MPa。絕對(duì)瓦斯涌出量為2.45 m3/min,平均二氧化碳絕對(duì)涌出量為0.20 m3/min,煤層無自然發(fā)火傾向,無地溫異常區(qū)。該工作面水文地質(zhì)條件較簡單,涌水來源主要為3#煤層上覆砂巖、粉砂巖等裂隙水。巷道設(shè)計(jì)長度2 134.68 m,斷面寬為5.00 m,高為3.40 m。3206回順掘進(jìn)期間采用前探梁臨時(shí)支護(hù)或錨桿鉆車臨時(shí)支護(hù)。煤層頂?shù)装迩闆r如表1所示。

表1 煤層頂?shù)装迩闆rTable 1 Roof and floor of coal seam

2 數(shù)學(xué)模型的選取

2.1 風(fēng)流相數(shù)學(xué)模型的選取

掘進(jìn)工作面壓入式通風(fēng)屬于有限空間受限附壁射流,選用Realizableκ-ε模型能較好地處理模擬中的流動(dòng),模型求解的κ和ε方程如下:

(1)

式中:ρ為流體密度,kg/m3;κ為湍動(dòng)能,m2/s3;t為時(shí)間,s;ui為瞬間速度,m/s;xi、xj為張量坐標(biāo)x方向的速度,m/s;μ為流體的分子黏性系數(shù);μt為渦黏性;σκ為κ方程的紊流普特朗數(shù),默認(rèn)值為1.0;Gκ為平均速度梯度引起的湍動(dòng)能產(chǎn)生;Gb為浮力影響引起的湍動(dòng)能產(chǎn)生;ε為紊流耗散率,m2/s3。

(2)

式中:σε為湍動(dòng)耗散率對(duì)應(yīng)的普朗特?cái)?shù),取1.2;C1ε、C2ε、C3ε為經(jīng)驗(yàn)常數(shù),FLUENT中默認(rèn)值為C1ε=1.44,C2ε=1.92,C3ε=0.09。

(3)

式中:C1為常數(shù);η為氣體動(dòng)力黏度,Pa·s;S為顆粒實(shí)際表面積,m2。

2.2 粉塵相數(shù)學(xué)模型的選取

根據(jù)“作用在塵埃顆粒上的力是平衡的”原理,拉格朗日坐標(biāo)系中的運(yùn)動(dòng)方程可導(dǎo)出下式:

(4)

式中:up為顆粒速度,m/s;FD為曳力,N;u為相續(xù)連速度,m/s;gx為X方向重力加速度,m/s2;ρp為顆粒密度,kg/m3;Fx為X方向的其他作用力,N。

流體對(duì)粉塵的拖曳阻力計(jì)算公式為:

(5)

式中:dp為顆粒直徑,m;CD為曳力系數(shù);Re為相對(duì)雷諾數(shù)。

相對(duì)雷諾數(shù)Re表達(dá)式為:

(6)

曳力系數(shù)CD的表達(dá)式為:

(7)

式中:α1,α2,α3為常數(shù)。

3 物理模型的建立與網(wǎng)格劃分

3.1 物理模型的建立

通過對(duì)王坡煤礦掘進(jìn)工作面進(jìn)行實(shí)地調(diào)查,并結(jié)合該礦實(shí)際布置情況,利用Solidworks軟件建立了1∶1仿真物理模型。物理模型主要由掘進(jìn)工作面、壓風(fēng)筒、掘進(jìn)機(jī)、轉(zhuǎn)載機(jī)、輸送機(jī)5部分組成。掘進(jìn)機(jī)為EBZ160型掘進(jìn)機(jī),轉(zhuǎn)載機(jī)為DZQ-80/30/11橋式轉(zhuǎn)載機(jī),輸送機(jī)為DTL80/20/22帶式輸送機(jī)。其中綜掘巷道為長80 m、寬5 m、高3.4 m的長方體,掘進(jìn)機(jī)長寬高分別為9.3 m、2.9 m、1.68 m,風(fēng)筒為直徑1 m的圓柱體。掘進(jìn)工作面物理模型如圖1所示。

圖1 掘進(jìn)工作面物理模型Fig.1 Physical model of excavation working face

3.2 網(wǎng)格劃分及獨(dú)立性測試

在數(shù)值模擬過程中,網(wǎng)格獨(dú)立性對(duì)模擬結(jié)果的準(zhǔn)確性非常重要,在進(jìn)行數(shù)值模擬之前應(yīng)先進(jìn)行網(wǎng)格獨(dú)立性測試[16-17]。利用ICEM CFD生成了3種不同質(zhì)量的非結(jié)構(gòu)網(wǎng)格,分別命名為網(wǎng)格I、II、III。利用FLUENT軟件分別對(duì)3種不同質(zhì)量的網(wǎng)格進(jìn)行了風(fēng)流模擬,結(jié)果如圖2(a)所示。網(wǎng)格II和網(wǎng)格III產(chǎn)生了幾乎相同的模擬結(jié)果,當(dāng)使用網(wǎng)格I時(shí)獲得的數(shù)據(jù)差異較大?？紤]模擬結(jié)果、模型形成程度和美觀程度后選擇了網(wǎng)格III進(jìn)行后續(xù)研究。進(jìn)行網(wǎng)格劃分后的物理模型如圖2(b)所示。劃分后的網(wǎng)格總數(shù)為1 638 017,最大網(wǎng)格質(zhì)量為0.999 75,最小網(wǎng)格質(zhì)量為0.292 237,平均網(wǎng)格質(zhì)量為0.745。其中,99.99%的網(wǎng)格質(zhì)量在0.4以上,網(wǎng)格沒有負(fù)值,網(wǎng)格質(zhì)量良好,符合模擬計(jì)算精度和收斂速度的條件。網(wǎng)格質(zhì)量分布如圖3所示。氣流和粉塵的計(jì)算可以分別看作是穩(wěn)態(tài)解和非穩(wěn)態(tài)解,具體參數(shù)設(shè)置如表2所示。

圖2 網(wǎng)格劃分及獨(dú)立性測試Fig.2 Grid division and independence test

圖3 網(wǎng)格質(zhì)量分布Fig.3 Grid mass distribution

表2 數(shù)值模擬的參數(shù)設(shè)置Table 2 Parameters of numerical simulation

4 綜掘工作面壓入式通風(fēng)模擬結(jié)果及分析

4.1 風(fēng)流流場數(shù)值模擬

運(yùn)用FLUENT數(shù)值模擬軟件對(duì)掘進(jìn)工作面的風(fēng)流流場運(yùn)移情況進(jìn)行模擬,并運(yùn)用CFD-POST軟件對(duì)模擬結(jié)果進(jìn)行后處理,風(fēng)流流場數(shù)值模擬結(jié)果如圖4所示。由圖4可知:

1)風(fēng)流以27.6 m/s的速度從壓風(fēng)口排出,在射流卷吸作用下沿著巷道壁面向前運(yùn)移,在與工作面相距0～5 m處形成附壁射流。射流在前進(jìn)過程中不斷卷吸周邊空氣,使向工作面運(yùn)移的射流區(qū)域范圍呈錐形不斷擴(kuò)大。其中,風(fēng)筒出口處的風(fēng)流不受空間限制呈自由射流狀態(tài),靠近壁面?zhèn)鹊娘L(fēng)流受到壁面阻擋沿壁面向前運(yùn)移,遠(yuǎn)離壁面?zhèn)鹊娘L(fēng)流卷吸周圍空氣不斷向外擴(kuò)散。

2)由于射流卷吸作用及受空氣阻力的影響,風(fēng)流從壓風(fēng)口向工作面運(yùn)移過程中速度不斷減小,在到達(dá)工作面時(shí)速度衰減為13.6 m/s,該速度滿足工作面除塵所需風(fēng)速要求。

3)當(dāng)射流到達(dá)迎頭時(shí)撞擊工作面,因受工作面的阻擋導(dǎo)致風(fēng)流方向發(fā)生改變,改變后的風(fēng)流方向大致分為3部分:一部分流向下部底板,一部分流向工作面頂板,絕大部分風(fēng)流在巷道風(fēng)筒側(cè)流向巷道出口方向,進(jìn)而在巷道遠(yuǎn)離風(fēng)筒側(cè)產(chǎn)生回流區(qū)域。

4)由于綜掘巷道空間有限,在射流卷吸作用和回流區(qū)的共同作用下,風(fēng)流流場在距工作面0～5 m、10～16 m及18～24 m處形成了3處閉合狀渦流區(qū):第一處特征明顯的渦流區(qū)出現(xiàn)在掘進(jìn)機(jī)截割頭附近,由于掘進(jìn)機(jī)體積龐大,回流區(qū)風(fēng)流在移動(dòng)過程中受到掘進(jìn)機(jī)機(jī)身的阻擋作用,使得巷道1.7 m以下的風(fēng)流流場遭到破環(huán),從而造成截割頭與工作面之間風(fēng)流場紊亂,且形成了較為明顯的渦流區(qū);第二處渦流區(qū)出現(xiàn)在距工作面10～16 m的掘進(jìn)機(jī)尾部處;因受轉(zhuǎn)載機(jī)及運(yùn)輸機(jī)的影響,在距工作面18～24 m處,即轉(zhuǎn)載機(jī)和運(yùn)輸機(jī)之間形成了第三處范圍較大且形狀不規(guī)則的渦流區(qū)。

5)在距工作面40 m后,風(fēng)流己基本脫離前方壓風(fēng)射流場及渦形風(fēng)流場的影響,風(fēng)流跡線逐漸平滑且保持穩(wěn)定,形成風(fēng)速較小(1.27 m/s左右)、風(fēng)流跡線垂直于巷道出口的平穩(wěn)流。

圖4 風(fēng)流流場分布Fig.4 Air flow field distribution

4.2 粉塵運(yùn)移數(shù)值模擬及時(shí)空演化規(guī)律

對(duì)壓入式通風(fēng)條件下巷道的粉塵運(yùn)移情況進(jìn)行了模擬,并對(duì)模擬結(jié)果進(jìn)行后處理。圖5為不同時(shí)刻下粉塵的擴(kuò)散情況。由圖5可知:

1)單壓入式通風(fēng)條件下掘進(jìn)巷道的風(fēng)流流場可分為射流區(qū)、回流區(qū)、渦流區(qū)3個(gè)區(qū)域,而巷道內(nèi)粉塵質(zhì)量濃度分布和風(fēng)流流場分布一致。高速射流從壓風(fēng)筒流出后迅速將工作面產(chǎn)生的粉塵沖散,粉塵在風(fēng)流的作用下開始向巷道出口方向擴(kuò)散。

2)工作面附近風(fēng)流流場復(fù)雜,存在一處風(fēng)速較小的渦流區(qū)域。顆粒較大的粉塵在渦流及自身重力的作用下發(fā)生沉降,造成掘進(jìn)機(jī)前方粉塵積聚。當(dāng)t=4 s時(shí),部分粉塵已逐漸擴(kuò)散到掘進(jìn)機(jī)尾部的渦流區(qū)域,并開始運(yùn)移至轉(zhuǎn)載機(jī)與運(yùn)輸機(jī)附近的渦流區(qū)。由于渦流區(qū)域攜帶能力差,導(dǎo)致顆粒較大的粉塵在渦流區(qū)沉降,造成該區(qū)域粉塵質(zhì)量濃度較高,達(dá)到613 mg/m3左右。

3)在t=8 s時(shí),粉塵開始脫離渦流區(qū)域,進(jìn)入40～80 m的平穩(wěn)風(fēng)場,并向風(fēng)場較穩(wěn)定的區(qū)域擴(kuò)散,此時(shí)粉塵的運(yùn)動(dòng)軌跡相對(duì)平緩,質(zhì)量濃度超過350 mg/m3的較高濃度粉塵主要集中在遠(yuǎn)離風(fēng)筒的一側(cè)及掘進(jìn)機(jī)運(yùn)輸機(jī)之間的渦流區(qū)域。在t=48 s時(shí),粉塵在風(fēng)流的作用下初次到達(dá)巷道出口,且在t=145 s時(shí)粉塵質(zhì)量濃度在巷道中達(dá)到穩(wěn)定狀態(tài)。

從整體上看,隨著距地板距離的增加,粉塵質(zhì)量濃度變小,這是由于大顆粒粉塵受重力和阻力等因素的影響在逐漸沉積。巷道壓風(fēng)筒一側(cè)粉塵質(zhì)量濃度較低,質(zhì)量濃度超過350 mg/m3的高濃度粉塵主要集中在巷道遠(yuǎn)離風(fēng)筒的一側(cè),說明粉塵在運(yùn)移過程中主要受風(fēng)流的影響。

4.3 不同壓風(fēng)量下粉塵運(yùn)移規(guī)律

為了研究不同壓風(fēng)量對(duì)通風(fēng)控塵的影響,從而尋找通風(fēng)控塵的最佳通風(fēng)量,選取了1 300 m3/min、1 350 m3/min、1 400 m3/min、1 450 m3/min、1 500 m3/min共計(jì)5組壓風(fēng)量進(jìn)行模擬,觀察巷道內(nèi)粉塵運(yùn)移情況。5種不同通風(fēng)量下巷道內(nèi)的粉塵情況,如圖6所示。

圖5 不同時(shí)刻下粉塵運(yùn)移情況Fig.5 Dust transportation at different times

圖6 不同壓風(fēng)量下巷道內(nèi)粉塵濃度分布Fig.6 Dust concentration distribution in roadway under different pressure air volume

采用原始方案時(shí),即壓風(fēng)量為1 300 m3/min時(shí),整個(gè)巷道內(nèi)充滿了高質(zhì)量濃度粉塵,給掘進(jìn)巷道清潔生產(chǎn)工作造成了阻礙。在壓風(fēng)量從1 300 m3/min提升到1 400 m3/min的過程中,巷道內(nèi)整體粉塵質(zhì)量濃度明顯降低,這是由于隨著風(fēng)量的增加,巷道內(nèi)風(fēng)速增大,從而粉塵攜帶能力增大,粉塵在風(fēng)流的卷吸作用下攜帶到巷道后方,此時(shí)巷道內(nèi)空氣流通較好。但壓風(fēng)量繼續(xù)增加時(shí),即壓風(fēng)量從1 400 m3/min增加到1 500 m3/min過程中,巷道內(nèi)粉塵質(zhì)量濃度明顯增高。這是由于壓風(fēng)量不斷增加,導(dǎo)致巷道內(nèi)風(fēng)速過大,造成巷道二次揚(yáng)塵,所以巷道內(nèi)粉塵質(zhì)量濃度明顯增高。當(dāng)壓風(fēng)量為1 400 m3/min時(shí)通風(fēng)控塵效果相對(duì)較好,巷道內(nèi)粉塵質(zhì)量濃度較低。所以選取壓風(fēng)量1 400 m3/min為最佳控塵風(fēng)量。

5 現(xiàn)場實(shí)測結(jié)果及分析

為了驗(yàn)證模擬結(jié)果的準(zhǔn)確性和科學(xué)性,根據(jù)實(shí)際情況,在王坡煤礦3206工作面現(xiàn)場選取測試點(diǎn),測試點(diǎn)設(shè)置(X、Y、Z)坐標(biāo),Y值為測量距巷道地板的垂直距離,Y值選取呼吸區(qū)高度,即Y=1.55 m;Z值為測量點(diǎn)距風(fēng)筒側(cè)巷道壁面水平距離,Z值選取兩個(gè)距離,分別為1.5 m和4.5 m,測量點(diǎn)坐標(biāo)可表示為#A(X,1.55,1.50)、#B(X,1.55,4.50);X為測量點(diǎn)距工作面的水平距離,X值分別取5 m、10 m、30 m、50 m、70 m。選取最佳壓風(fēng)量1 400 m3/min進(jìn)行現(xiàn)場測量,待巷道風(fēng)流場穩(wěn)定后,對(duì)各個(gè)測點(diǎn)的風(fēng)流大小及粉塵質(zhì)量濃度進(jìn)行了實(shí)測,并與數(shù)值分析的結(jié)果進(jìn)行比較,計(jì)算出兩者的誤差范圍。測試點(diǎn)分布如圖7所示,各測點(diǎn)風(fēng)速、粉塵質(zhì)量濃度實(shí)測值與模擬值對(duì)比如表3和表4所示。由表中數(shù)據(jù)可知,各測點(diǎn)風(fēng)速及粉塵質(zhì)量濃度實(shí)測值與模擬值的相對(duì)誤差控制在10%以內(nèi),證明了模擬的有效性。

表3 #A(X,1.55,1.5)各測點(diǎn)風(fēng)速及粉塵質(zhì)量濃度實(shí)測值與模擬值對(duì)比Table 3 Comparison of measured and simulated values of wind velocity and dust concentration at measuring points in section A (X, 1.55, 1.5)

表4 #B(X,1.55,4.5)各測點(diǎn)風(fēng)速及粉塵質(zhì)量濃度實(shí)測值與模擬值對(duì)比Table 4 Comparison of the measured and simulated values of wind velocity and dust concentration at measuring points in section B (X, 1.55, 1.5)

6 結(jié)論

1)壓入式通風(fēng)條件下巷道內(nèi)風(fēng)流流場主要分為射流區(qū)、渦流區(qū)、回流區(qū)3個(gè)區(qū)域。射流區(qū)風(fēng)速最高且風(fēng)速降低較快;巷道內(nèi)存在3處特征比較明顯的渦流,分別位于工作面附近、掘進(jìn)機(jī)后方、轉(zhuǎn)載機(jī)與運(yùn)輸機(jī)連接處,風(fēng)流流場在經(jīng)過渦流區(qū)域后速度呈大幅度衰減趨勢(shì);回流區(qū)風(fēng)速較為穩(wěn)定,速度衰減較慢,最后風(fēng)流以1.26 m/s的速度流出巷道。

2)壓入式通風(fēng)條件下巷道內(nèi)質(zhì)量濃度超過350 mg/m3的粉塵主要積聚在掘進(jìn)機(jī)與工作面之間以及渦流區(qū)域及巷道內(nèi)遠(yuǎn)離風(fēng)筒的一側(cè),與風(fēng)流流場射流區(qū)、渦流區(qū)、回流區(qū)3個(gè)區(qū)域發(fā)生位置保持一致,說明了粉塵運(yùn)動(dòng)主要受風(fēng)流的影響。

3)提高風(fēng)筒壓風(fēng)量在一定程度上可以提高通風(fēng)控塵的效果,但壓風(fēng)量過大會(huì)造成巷道內(nèi)二次揚(yáng)塵。研究發(fā)現(xiàn)當(dāng)壓風(fēng)量為1 400 m3/min時(shí)通風(fēng)控塵效果最佳。