黃 旭，趙 勉

(1.隴東學院能源工程學院，甘肅慶陽745099；2.中煤能源山東有限公司，山東日照276826)

基于平面紊動射流的風幕控塵裝置研究

黃 旭1，趙 勉2

(1.隴東學院能源工程學院，甘肅慶陽745099；2.中煤能源山東有限公司，山東日照276826)

礦井綜掘工作面因其作業(yè)空間狹小，作業(yè)強度大，而使粉塵難以擴散、沉降，濃度遠超國家規(guī)定，嚴重影響礦井安全和作業(yè)人員的身心健康。風幕控塵技術與傳統(tǒng)的掘進面除塵技術相比有著諸多優(yōu)勢。依據(jù)平面自由紊動射流理論，對研制出的一種新型風幕控塵裝置進行礦井適應性改進，應用CFD技術模擬巷道在加設該裝置前后的粉塵濃度變化，證明了這種改進型風幕控塵裝置的實用性。

紊動射流；CFD技術；風幕控塵；控塵

1 風幕控塵原理及裝置結構

綜掘面的風幕裝置是將風幕射流裝置安設工作面與掘進司機之間的掘進機上，利用動力裝置將空氣壓入，通過出射條縫形成一道風幕墻，相當于在綜掘司機前方加設一個附加阻力層[1]。風幕控塵裝置需要設置壓入和抽出兩個模塊，在壓入式風筒提供新鮮風流的同時，抽出式風筒將工作面附近高濃度的粉塵抽出，并使污染區(qū)內(nèi)形成負壓，在風幕裝置出口處形成一個曲面風幕[2]，更好地抑制粉塵的外散。對粒徑小于5μm的呼吸性粉塵[3]，其效果更佳。

遼寧工程技術大學李雨成、楊靖等人[2,3,4]根據(jù)相應原則，對風幕控塵裝置的結構和參數(shù)進行了相關設計。

在對風幕控塵射流裝置的研究過程中，射流條縫的長度遠遠大于寬度，噴射流能夠發(fā)展充分，流體所受壓力不足以壓縮流體，所以風幕控塵條縫射流應為平面自由紊動不可壓射流(以下簡稱平面自由紊動射流)。

2 平面自由紊動射流風速衰減規(guī)律

射流的整體結構是相當復雜的。射流主要可以分為初始段、過渡段及主體段三部分，但過渡段較短，可忽略不計。

表1 半擴展厚度系數(shù)與軸向最大時均速度衰變系數(shù)統(tǒng)計表

注：此處k1為be時半擴展厚度系數(shù)

表2 風幕射流裝置上側和右側k2的擬合值

對軸向最大時均速度衰變系數(shù)k2取均值，上側擬合均值為4.25，右側擬合均值為3.63。

右側擬合的軸向最大時均速度衰變系數(shù)k2在理論值的范圍內(nèi)，上側擬合的軸向最大時均速度衰變系數(shù)k2考慮誤差影響及測量位置與射流軸向最大速度處有偏差，故k2略大是允許的。應用擬合k2進行實際巷道中的射流裝置風量衰減情況的預測是完全可行的。

圖1 風幕射流上側擬合圖

圖2 風幕射流右側擬合圖

3 山西某礦新增輔助運輸巷風幕控塵裝置設計及布局

根據(jù)山西某礦掘進工作面中所使用的S200MJ掘進機的外形及布置位置參數(shù)：長14.3m、寬3.6m、高1.8m、搖臂直徑1030mm、搖臂長度5m、操作臺(掘進司機位置)距煤壁距離14.3m、掘進機左側至巷道左壁距離0.7m、機體前側至掘進面距離5m[3]。對原有設備尺寸進行改進，主要設計參數(shù)如表3所示。

表3 風幕控塵裝置主要參數(shù)指標

此巷道為矩形，寬5.2m，巷道高3.8m，風筒出口軸心至底板距離3m，巷道布置符合相關規(guī)定和標準。為了得到更好的除塵效果，就要對現(xiàn)有工作面加裝風幕控塵裝置，并對輔助運輸巷做如下改造：

(1)該裝置安設在掘進機司機前方。把射流箱安設在離工作面7m的掘進機上；

(2)風幕壓風風機應放置在特制的拖車上，根據(jù)工作面的推進而隨掘進機前移；

(3)安設抽出式風筒。在原巷道加設帶有除塵裝置的抽出式風筒，抽風量要大于壓風量。距工作面3～4m范圍內(nèi)是高濃度粉塵最集中的地方，所以把抽出式通風的風筒吸風口的位置定在離工作面4m的位置為宜。根據(jù)礦中現(xiàn)有材料，風筒直徑選擇600mm，軸心離底板位置2.1m或3.0m(下文確定最佳位置)，緊貼右側巷道壁。

射流末端的風流風速大小足以控制住粉塵即可。在進行最低控塵風速計算時，考慮到風速富余問題，半擴展厚度(比例)系數(shù)k1應選取較小值，根據(jù)相關學者的研究，取Abramovich給定的實驗結果k1=0.21進行末端最小控塵風速的計算。

對最小控塵風速的計算時，采用極端物理模型假設的方法：

設定出射的風流形成了一個垂直向上規(guī)整的平面矩形風幕，粉塵顆粒的運動速度v顆粒入射為巷道的實測風速1.15m/s，并在底部水平入射進入風幕區(qū)域，粉塵顆粒的出射速度v顆粒出射的縱向分量剛好達到風幕內(nèi)風流v射流的速度。

通過數(shù)學計算，風幕控塵裝置的最低控塵風速為5.48m/s，當風幕的射流到達巷道壁時的末速度為5.48m/s時，粉塵顆粒不能逃逸出風幕區(qū)域。依據(jù)射流軸向最大時均速度衰變公式和軸向最大時均速度衰變系數(shù)k2的擬合值，得到上側和兩側的最大時均速度衰變公式，并通過風量公式計算總需風量Q總=Q上+Q左+Q右=149.4m3/min。

實測壓入式風機的出口風量為1224.4m3/min。壓入式風機風量一部分用于工作面通風，另一部分需要供給風幕控塵裝置，所以壓入式風機實際用于工作面通風的風量為1075m3/min。為了使出射風流偏向工作面一側，理論抽出式風筒風量要大于1075m3/min。

4 風幕控塵裝置數(shù)值模擬研究

4.1 數(shù)學模型的建立

本文在FLUENT模擬中采用拉格朗日法所對相應的Discrete Phase Model模型。在連續(xù)相中，做溫度與壓強對流場影響不大；不考慮張力作用；不存在熱交換的假設(條件)。

對離散相做如下假設：把粉塵顆粒為密度均勻的球體，只考慮與巷道壁的作用及顆粒的重力。掘進面看成一個面塵源，根據(jù)對粉塵分散度的測定，其粒度范圍在0.5～21μm之間，中位徑為1.5μm。

由于掘進面的作業(yè)強度主要密集在工作面附近，且主要考慮風幕控塵裝置的效果，所以巷道取30m范圍內(nèi)進行研究。連續(xù)相為巷內(nèi)空氣，根據(jù)巷道內(nèi)實測數(shù)據(jù)(干溫度=濕溫度=14℃，壓力=102.954kPa)，計算出密度為1.242kg/m3，溫度為287K，粘性系數(shù)為1.759×10-5kg/(m·s)。粉塵顆粒密度為1200kg/m3，質量流量為0.0004kg/s。實測巷道內(nèi)風速為1.15m/s。壓入式風筒風量為1075m3/min，風速為35.67m/s，風幕出射口風速為12.04m/s。風筒及風幕射流出口雷諾數(shù)Re、水力直徑和湍流強度I見表4。

表4 風筒及風幕射流出口雷諾數(shù)Re、水力直徑和湍流強度I

通過Gambit軟件完成巷道幾何模型構建，并劃分網(wǎng)格。

4.2 軟件重要條件的設定及說明

將壓入式風筒出口和風幕條縫出射口為velocity-inlet；抽出式風筒前端風流入口和巷道后部截斷面為outflow，并為“escape”；巷道四周壁，風筒壁為wall，并為“reflect”，并使用RNGk-ε模型[6]。

4.3 CFD模擬與精度分析

取濃度測量范圍120m長的巷道，將氣固兩相流基本參數(shù)設置完畢后，利用FLUENT軟件開始進行迭代，經(jīng)過200步左右的計算，達到給定的收斂標準。

圖3 粉塵濃度數(shù)據(jù)模擬精度對比

將呼吸帶高度(距地面1.5m)實測總粉塵濃度與模擬值對比(如圖3所示)，由結果顯示實際測量濃度與模擬值之間存在偏差，這主要是由于測試與掘進機工作的相對時間不確定，粉塵在巷道中做不規(guī)則的紊動運動，不排除模擬時所取點渦流較大，造成偏差，但趨勢與實際相符。該模擬簡化了現(xiàn)場實際條件，真實條件下會有很多影響粉塵濃度的因素，但模擬與實測具有一致性，這說明應用給定的邊界條件進行粉塵濃度的模擬研究是可行的。

圖4 粉塵濃度相對數(shù)量級對比

圖5 無風幕裝置時粉塵的跡線圖與粉塵濃度云圖(a)無風幕裝置時粉塵的跡線圖；(b)z=5m、z=10m、z=15m、z=20m、z=25m、z=30m、y=1.5m；(c)z=5m、z=10m、z=15m、z=20m、z=25m、z=30m、x=2.6m

4.4 模擬結果及分析

在未加風幕控塵裝置兩種情況下，依據(jù)上文所建立的幾何模型及所設邊界條件，模擬得到粉塵顆粒跡線圖(如圖5(a))和z=5m、z=10m、z=15m、z=20m、z=25m、z=30m、y=1.5m、x=2.6m處濃度分布的切片云圖(如圖5(b)、圖5(c))。

由圖5可以看出，由于壓入式風筒的作用，新鮮風流將風筒前部的粉塵吹散、稀釋，所以在風筒側粉塵濃度較低，而另一側粉塵會在掘進機前方形成小渦旋，同時在壓入式風筒出風口與工作面之間會出現(xiàn)大渦旋，大部分粉塵顆粒會通過掘進機與巷道壁的間隙向整個巷道擴散，這是由于巷道與掘進機之間形成了窄縫，風流經(jīng)過此時，風速會變大，隨之引起壓強降低，形成低壓區(qū)，隨后會有更多的粉塵顆粒被壓入此處。當粉塵顆粒通過這個段后會由巷道中部分別向上和向下逐漸運移。

在加裝風幕控塵裝置兩種情況下，所設邊界條件與未加風幕時完全相同，將抽出式風機出口中心位置設置在距底板2.1m和3.0m處，得到粉塵顆粒跡線圖(如圖6)和z=5m、z=10m、z=15m、z=20m、z=25m、z=30m、y=1.5處濃度分布的切片云圖(如圖7)。

由圖6和圖7可以看出，由于抽出式風筒的作用，工作面前部粉塵濃度有所降低，在抽出式風筒一側形成的小渦旋變小，并向抽出式風筒移動，說明抽出式風機帶走了大量的粉塵顆粒。根據(jù)粉塵顆粒運動軌跡可以看出，一部分粉塵會經(jīng)風幕射流后折返。加裝風幕控塵裝置后，整個巷道粉塵濃度明顯降低，說明風幕控塵裝置起到了應有的作用，但仍難以避免一些高速粉塵顆粒穿過風幕裝置向整個巷道擴散。

加裝風幕控塵裝置后，對抽出式風機出口中心位置距底板2.1m和3.0m時巷道中粉塵濃度相同數(shù)量級的對比(如圖4所示)可以發(fā)現(xiàn)，當抽出式風機出口中心位置距底板2.1m時的效果要比3.0m時的控塵效果更好。

圖6 有風幕裝置時粉塵跡線圖

(a)抽出式風機風筒出口中心位置距底板2.1m，z=5m、z=10m、z=15m、z=20m、z=25m、z=30m、y=1.5m；(b)抽出式風機風筒出口中心位置距底板2.1m，z=5m、z=10m、z=15m、z=20m、z=25m、z=30m、x=2.6m；(c)抽出式風機風筒出口中心位置距底板3m，z=5m、z=10m、z=15m、z=20m、z=25m、z=30m、y=1.5m；(d)抽出式風機風筒出口中心位置距底板3m，z=5m、z=10、z=15m、z=20m、z=25m、z=30m、x=2.6m圖7 有風幕裝置時粉塵濃度云圖

[1]李帥君，熊守美，Mei Li，等.應用兩相流模型模擬壓鑄充型過程的卷氣現(xiàn)象[J].金屬學報，2009，45(10)：1153-1158.

[2]楊靖，李雨成.風幕集塵射流裝置研究[J].世界科技研究與發(fā)展，2013，35(5)：591-593.

[3]趙勉.綜掘面風幕射流形成條件及數(shù)值模擬研究[D].遼寧工程技術大學，2015：2-5.

[4]楊靖.綜掘面風幕控塵綜合實驗平臺設計研究[D].遼寧工程技術大學，2014：3-5.

[5]Kramer T. J. and Depew C.A.，Analysis of mena flow characteristics of gas-solid suspension[J].Tuns.ASME л basic Engng，1992(94)：731.

[6]趙國景，丁廣攘.采煤工作面含塵風流流動的數(shù)值模擬[C].第六屆全國計算流體力學討論會，1993：126-127.

【責任編輯 趙建萍】

A Research on Air Curtain Dust Control Devices Based on Planar Turbulent Jet

HUANG Xu1，ZHAO Mian2

(1.SchoolofEnergyEngineering,LongdongUniversity,Qingyang745000,Gansu; 2.ChinaCoalEnergyCo.,Ltd.Shandong,Rizhao276826,Shandong)

Because of its narrow working space,large amounts of work intensity,it is difficult to spread the dust concentration,sedimentation in mine fully mechanized workface,and the dust concentration is far exceeding the relevant state regulations,and seriously affects the safety of mine workers. Compared with the traditional driving surface dust removal technology,air curtain dust control technology has many advantages. In this paper,based on the free plane turbulent jet theory,a new type of air curtain dust control is developed which means mine adaptive improvement. Application of CFD simulation technology of roadway in Canada provided the dust concentration variation of the device before and after which shows the practicability of the improved air curtain dust control device.

turbulent jet;CFD technology;wind screen control dust;dust control

1674-1730(2017)03-0085-05

2016-05-24

2014年隴東學院青年基金項目《劉園子煤礦涌水量預測模型研究》(XYBY140201)

黃 旭(1986—)，男，遼寧阜新人，助教，碩士，主要從事采礦工程教學與研究。

TD714

A