蔣仲安，鄧權(quán)龍，時訓(xùn)先，陳舉師

(1.北京科技大學(xué) 金屬礦山高效開采與安全教育部重點(diǎn)實驗室，北京 100083; 2.中國安全生產(chǎn)科學(xué)研究院 工業(yè)安全研究所，北京 100012)

石棉篩分車間粉塵質(zhì)量濃度分布規(guī)律的數(shù)值模擬*

蔣仲安1?，鄧權(quán)龍1，時訓(xùn)先2，陳舉師1

(1.北京科技大學(xué) 金屬礦山高效開采與安全教育部重點(diǎn)實驗室，北京 100083; 2.中國安全生產(chǎn)科學(xué)研究院 工業(yè)安全研究所，北京 100012)

為了掌握石棉選礦廠篩分車間內(nèi)部粉塵隨時間和空間運(yùn)移的規(guī)律，獲取通風(fēng)除塵優(yōu)化的參數(shù)，以茫崖石棉礦第一選礦廠篩分車間為研究背景，依據(jù)氣固兩相流理論，采用計算機(jī)流體力學(xué)的離散相(DPM)模型對篩分車間粉塵濃度分布規(guī)律數(shù)值模擬，并與現(xiàn)場粉塵濃度實測數(shù)據(jù)比較分析，模擬結(jié)果和實測數(shù)據(jù)相吻合.研究結(jié)果表明：篩分設(shè)備附近粉塵濃度較高，以篩分設(shè)備為中心隨著距離增加粉塵濃度降低；捕集邊界下粉塵濃度比反彈邊界較低，粉塵捕集效果更好；篩分車間進(jìn)風(fēng)風(fēng)速控制在 0.6 m/s左右時，粉塵沉降效果較好；安裝抽風(fēng)集塵罩后，主抽風(fēng)管風(fēng)速為 14 m/s時，平均控塵效率達(dá)到90%，石棉纖維濃度控制在2.52 f/mL以下.

石棉篩分車間；粉塵濃度；運(yùn)移規(guī)律；離散相；數(shù)值模擬

石棉因其具有獨(dú)特的耐腐蝕性、隔熱性和強(qiáng)抗拉性等優(yōu)良特性，廣泛應(yīng)用于航天、船舶、消防、特殊制造等領(lǐng)域[1-2].國內(nèi)石棉礦產(chǎn)資源豐富，工業(yè)使用量較大，2014年使用量達(dá)到50萬噸.近幾年歐美、德國、日本等國家禁止生產(chǎn)含石棉成分的產(chǎn)品；國內(nèi)采取全面禁止使用青石棉，安全使用溫石棉的方針[3].石棉在開采、選礦、加工等過程容易產(chǎn)生粉塵，污染新鮮空氣，威脅作業(yè)人員的健康，特別是粉塵含有微細(xì)纖維，人體器官組織吸入后易附著和沉降，導(dǎo)致石棉肺、肺癌、間皮瘤等疾病[4-6].篩分作業(yè)作為選礦工藝中重要的工序，主要采用振動篩分設(shè)備，通過網(wǎng)篩不停的振動，分離粒徑大于網(wǎng)篩孔徑的塊石，由于網(wǎng)篩循環(huán)的震蕩，致使礦石表面的粉塵揚(yáng)起，擴(kuò)散到車間污染空氣.目前，國內(nèi)外對粉塵控制研究主要采用理論分析、相似實驗、數(shù)值模擬等方法研究粉塵濃度分布規(guī)律.國內(nèi)學(xué)者的研究主要集中在礦山作業(yè)面[7-8]，在石棉篩分車間粉塵方面的研究較少.國外學(xué)者在工業(yè)粉塵控制做了一定的研究，例如Silvester等[9]針對壓碎破碎設(shè)備車間粉塵控制，通過計算流體動力學(xué)軟件預(yù)測車間粉塵分布；Constance[10]分別對稀釋通風(fēng)和局部通風(fēng)兩種控塵措施的優(yōu)缺點(diǎn)進(jìn)行概述.本文針對石棉篩分車間建立幾何模型，通過數(shù)值模擬車間不同條件下粉塵的運(yùn)移規(guī)律，結(jié)合現(xiàn)場實測數(shù)據(jù)，獲得篩分車間最優(yōu)通風(fēng)控塵參數(shù)，指導(dǎo)現(xiàn)場控塵改造，對于保障作業(yè)人員身體健康有重要意義.

1 數(shù)學(xué)模型與控制方程

1.1 數(shù)學(xué)模型

粉塵粒子在空氣運(yùn)移的理論基礎(chǔ)是氣固兩相流理論，目前，在數(shù)值計算方面主要分兩種方法：歐拉-歐拉方法和歐拉-拉格朗日方法.其中，歐拉-拉格朗日方法是運(yùn)用歐拉法、拉格朗日法分別處理背景、離散流體，將氣體作為背景流體，粉塵粒子視為離散流體，離散相模型(DPM)的理論基礎(chǔ)就是基于該方法.本文運(yùn)用離散相模型研究石棉選廠篩分車間粉塵運(yùn)動規(guī)律，將車間內(nèi)部空氣視為背景流體，采用歐拉法進(jìn)行研究；將塵顆粒作為離散流體，采用拉格朗日法研究粉塵粒子運(yùn)移線路，最終得出粉塵在車間空氣中的濃度分布.

1.2 控制方程

假定流體為不可壓縮的牛頓流體，忽略體積力，則連續(xù)相的運(yùn)動方程選用三維不可壓縮Navier-Stokes方程：

連續(xù)方程：

(1)

動量守恒方程：

(2)

式中:ρ為氣體密度，kg/m3；τij為應(yīng)力張量；p為流體相壓力，Pa；g為重力加速度，m/s2；xi,xj分別為x,y方向上的坐標(biāo)，m；ui,uj分別為流體在x,y方向上的速度(m/s)；Fi為控制體平均的顆粒流體阻力，N.

k-ε湍流動能方程：

Gk-ρε+Sk

(3)

(4)

通過積分拉格朗日坐標(biāo)下顆粒作用力的微分方程來求解顆粒的軌跡，在求解過程中，只考慮阻力和重力及顆粒相的作用力[11-13].

(5)

式中:FD(u-up)為顆粒的單位質(zhì)量阻力；FC為顆粒與顆粒(邊界)碰撞力；FD可根據(jù)式(4)計算：

(6)

2 幾何建模與參數(shù)設(shè)置

2.1 幾何模型與網(wǎng)格劃分

石棉選廠篩分車間主要有供料皮帶間、篩分設(shè)備間、物料分類收集間3部分，幾何模型圖如圖1所示.設(shè)備間長12 m、寬9 m、高4 m;車間安裝2臺篩分機(jī)，篩分機(jī)尺寸為長4 m、寬1 m、高1.2 m;開有兩扇門，高3 m、寬2 m;給料管段斷面為0.8 m×0.6 m.集塵罩整體形狀為規(guī)則的六面體，大罩口斷面1.0 m×0.6 m，集塵罩小口直接與抽風(fēng)通連接.采用Gambit軟件劃分模型網(wǎng)格，詳見圖2.

圖1 篩分車間幾何模型圖Fig.1 The geometrical model of screening workshop

圖2 篩分車間網(wǎng)格劃分Fig. 2 The meshing of screening workshop

2.2 計算模型參數(shù)設(shè)置

根據(jù)茫崖石棉礦第一選廠篩分車間的具體情況和相關(guān)實測數(shù)據(jù)，結(jié)合CFD相關(guān)理論，對FLUENT數(shù)值模擬進(jìn)行參數(shù)設(shè)定[14-16]，詳見表1.

表1 計算模型參數(shù)設(shè)定表Tab.1 The calculation model parameters

3 數(shù)值模擬結(jié)果與對比分析

3.1 風(fēng)流流場分析

設(shè)置車間進(jìn)風(fēng)口風(fēng)速為 0.25 m/s，通過對車間風(fēng)流流場進(jìn)行數(shù)值模擬，得到風(fēng)流場速度矢量圖和風(fēng)流速度云圖，如圖3，圖4所示.從圖中看出：①風(fēng)流進(jìn)入車間先沿著風(fēng)流初始方向向前運(yùn)動，隨后前方受到車間壁面邊界的阻擋后，風(fēng)向幾乎發(fā)生90°變化流向車間出口.②風(fēng)流流場影響區(qū)域呈現(xiàn)不斷擴(kuò)大的趨勢，流場風(fēng)速則是呈現(xiàn)減小的趨勢，風(fēng)流受附壁效應(yīng)的影響，流向車間中央風(fēng)流較少，形成漩渦區(qū)，風(fēng)流交換周期較長.③根據(jù)風(fēng)流狀態(tài)不同分為：射流區(qū)、渦流區(qū)、回流區(qū)、出口區(qū)，射流區(qū)和出口區(qū)風(fēng)速在0.20～0.25 m/s，回流區(qū)風(fēng)速保持在 0.10～0.20 m/s，渦流區(qū)風(fēng)速范圍 0～0.10 m/s.

圖3 篩分車間風(fēng)流場速度矢量分布圖Fig. 3 The distribution of velocity vector in screening workshop

圖4 篩分車間風(fēng)流速度云圖Fig. 4 The speed cloud in screening workshop

3.2 粉塵運(yùn)動規(guī)律

圖5 粉塵顆粒運(yùn)動軌跡圖Fig. 5 The particle trajectory of dust particles

圖6 篩分車間粉塵濃度空間分布圖Fig. 6 The dust concentration distribution in screening workshop

為直觀觀察粉塵粒子在篩分車間內(nèi)的運(yùn)動軌跡，在離散相模型設(shè)置篩分機(jī)物料出口處隨機(jī)產(chǎn)生120個粉塵粒子，運(yùn)用隨機(jī)游走模型(DRW)追蹤粒子運(yùn)動軌跡，得到圖5.圖6，圖7分別為車間粉塵質(zhì)量濃度空間分布圖和分布云圖.從圖5～圖7得出：①粉塵從塵源產(chǎn)生后隨風(fēng)流運(yùn)動方向進(jìn)行逸散，粉塵受到壁面邊界阻擋和碰撞后停止運(yùn)動.②在篩分設(shè)備物料出口產(chǎn)生的粉塵濃度達(dá)到最大值，隨著距離該區(qū)域距離增大粉塵濃度呈降低趨勢.③粉塵受車間壁面邊界阻擋和碰撞等效應(yīng)，同時部分粉塵隨風(fēng)流排出，在不同程度上起到降低粉塵的作用，但車間內(nèi)粉塵濃度依然較大，平均濃度約 65 mg/m3，篩分設(shè)備物料出口濃度高達(dá) 186.5 mg/m3，中央漩渦區(qū)濃度低于 28 mg/m3且粉塵停留時間較長，排塵與沉降效果差.

3.3 不同邊界條件粉塵濃度

為比較不同邊界條件篩分車間粉塵運(yùn)動和濃度分布規(guī)律的差異性，分別模擬了邊界條件設(shè)置為捕集和反彈兩種情況，圖7，圖8分別為篩分車間邊界設(shè)置為捕集條件、反彈條件粉塵濃度分布云圖.經(jīng)對比模擬結(jié)果發(fā)現(xiàn)：在粉塵濃度分布特征上兩者基本一致，在數(shù)值上反彈邊界情況較大，由此可見，在篩分車間日常作業(yè)過程中，現(xiàn)有工藝允許條件下，采取適當(dāng)?shù)貪櫇駢γ?、地板的措施，加?qiáng)對壁面的捕集效率，有助于減輕篩分車間粉塵污染.

圖7 篩分車間粉塵濃度分布云圖Fig. 7 The dust concentration cloud in screening workshop

圖8 邊界設(shè)置反彈時粉塵濃度云圖Fig. 8 The dust concentration cloud under reflect wall condition

3.4 不同車間進(jìn)風(fēng)風(fēng)速粉塵濃度

研究篩分車間不同進(jìn)風(fēng)風(fēng)速條件粉塵濃度空間分布情況，選取了車間進(jìn)風(fēng)風(fēng)速為0.20 m/s,0.40 m/s,0.60 m/s,0.80 m/s 4種情況下，換算成進(jìn)風(fēng)風(fēng)量分別對應(yīng)為1.2 m3/s,2.4 m3/s,3.6 m3/s,4.8 m3/s，分別進(jìn)行模擬，模擬結(jié)果詳見圖9.為了量化不同進(jìn)風(fēng)風(fēng)速排塵效果，對車間內(nèi)人體呼吸帶高度(H=1.5 m)平面粉塵濃度取平均值，進(jìn)風(fēng)風(fēng)速為0.20 m/s,0.40 m/s,0.60 m/s,0.80 m/s 4種情況下對應(yīng)的呼吸帶高度粉塵濃度平均值分別為97.58 mg/m3,82.42 mg/m3,70.87 mg/m3,67.21 mg/m3.從模擬結(jié)果看出：當(dāng)進(jìn)風(fēng)風(fēng)速在 0.20～0.80 m/s區(qū)間，篩分車間排塵效果隨著風(fēng)速變大而增強(qiáng)；當(dāng)進(jìn)風(fēng)風(fēng)速從0.60 m/s增加到0.80 m/s時，排塵效果增幅較小.考慮到當(dāng)風(fēng)速較大時，風(fēng)流流動容易揚(yáng)起車間地面積塵，產(chǎn)生二次揚(yáng)塵；當(dāng)風(fēng)速較小時，車間內(nèi)部排塵效果不好.因此，綜合考慮，建議篩分車間進(jìn)風(fēng)風(fēng)速控制在0.6 m/s左右.

圖9 不同進(jìn)風(fēng)風(fēng)速下粉塵濃度分布云圖Fig. 9 The dust concentration cloud at different inlet wind speeds

3.5 安裝抽風(fēng)集塵罩粉塵濃度分布

為改善車間作業(yè)環(huán)境，在兩臺篩分設(shè)備的物料出口處安裝集塵罩，用直徑為400 mm的風(fēng)管與集塵罩連通，模擬主抽風(fēng)管風(fēng)速為10 m/s,12 m/s,14 m/s,16 m/s 4種情況車間粉塵濃度分布，車間進(jìn)風(fēng)風(fēng)速為0.6 m/s，控塵效果如圖10所示.對比分析得出：安裝抽風(fēng)集塵罩后，大部分粉塵經(jīng)集塵罩和抽風(fēng)管排出，車間內(nèi)部粉塵濃度下降明顯，粉塵污染得到良好的控制.隨著主抽風(fēng)管風(fēng)速增加，控塵效果先顯著增強(qiáng)，后趨于穩(wěn)定，當(dāng)抽風(fēng)風(fēng)速達(dá)到14 m/s時，抽風(fēng)集塵罩基本能控制篩分設(shè)備物料出口處粉塵的外逸，控塵效果較好，當(dāng)增加抽風(fēng)風(fēng)速至16 m/s時，控塵效果趨于穩(wěn)定.綜合篩分車間控塵效果和風(fēng)機(jī)耗能等因素，最終確定抽風(fēng)集成罩主抽風(fēng)管適宜風(fēng)速為14 m/s.

圖10 不同抽風(fēng)風(fēng)速篩分車間控塵效果圖Fig. 10 The dust removal effect at different extraction wind speeds

4 實測結(jié)果與模擬結(jié)果對比分析

依據(jù)GBZ/T192.1—2007，GBZ/T192.5—2007以及相關(guān)文獻(xiàn)[17-18]，結(jié)合篩分車間空間結(jié)構(gòu)和產(chǎn)塵源位置情況，對車間內(nèi)呼吸帶高度布置10個測點(diǎn)，測點(diǎn)布置如圖11所示.現(xiàn)場采用IFC-2防爆型粉塵采樣儀，選用過氯乙烯纖維濾膜(孔徑0.8 μm)分別對集塵罩安裝前后粉塵進(jìn)行采樣.實驗室使用BK5000相差顯微鏡對樣本石棉纖維濃度進(jìn)行測定，數(shù)據(jù)結(jié)果如表2所示，采樣濾膜處理后在顯微鏡下觀察石棉纖維形態(tài)，如圖12所示.運(yùn)用濾膜稱重法對樣本進(jìn)行粉塵濃度測定，每個測點(diǎn)采用3次以上重復(fù)操作取平均值，并將實測數(shù)據(jù)與模擬結(jié)果進(jìn)行對比，對集塵罩安裝前后的控塵效率進(jìn)行計算，如圖13所示.

圖11 粉塵測點(diǎn)布置圖Fig. 11 The layout of dust measuring points

測定數(shù)據(jù)測點(diǎn)編號12345678910集塵罩安裝前石棉纖維濃度/(f·mL-1)23.2120.8314.1612.4417.6019.826.229.864.113.26集塵罩安裝后石棉纖維濃度/(f·mL-1)2.522.341.551.421.862.100.961.110.640.58

圖12 采樣濾膜石棉纖維在顯微鏡下放大視野圖Fig.12 The pictures of asbestos fibers on the sampling membrane under microscope

圖13 抽風(fēng)集塵罩安裝前后粉塵濃度實測與模擬數(shù)據(jù)Fig.13 The measured and simulated data of dust concentration before and after installation of collector

如圖12為石棉纖維在顯微鏡下放大的視野圖，石棉纖維的形態(tài)呈現(xiàn)細(xì)條針絲狀，而巖塵顆粒為細(xì)小顆粒狀，石棉粉塵中還有微細(xì)纖維，人體吸入體內(nèi)后，容易附著于器官組織引起病變.由表2，圖13可以看出：①篩分車間安裝集塵罩前后模擬結(jié)果和實測數(shù)據(jù)基本吻合，說明幾何建模和離散相模型(DPM)等相關(guān)的設(shè)置均合理，數(shù)值模擬結(jié)果具有一定的可靠性.②經(jīng)過對比集塵罩安裝前后的各測點(diǎn)粉塵濃度數(shù)據(jù)，發(fā)現(xiàn)安裝集塵罩后，粉塵通過風(fēng)管抽走，篩分車間內(nèi)部粉塵急劇下降，各測點(diǎn)的平均控塵效率達(dá)到90%.③集塵罩安裝后，篩分車間粉塵濃度降低的同時，粉塵中石棉纖維濃度也得到了良好的控制，1號測點(diǎn)石棉纖維濃度由原來最高的23.21 f/mL下降至2.52 f/mL.圖14為抽風(fēng)集塵罩安裝前后篩分車間的現(xiàn)場照片.

(a)安裝前 (b)安裝后圖14 抽風(fēng)集塵罩安裝前后車間作業(yè)時現(xiàn)場照片F(xiàn)ig. 14 The workshop photos before and after the installation of the exhaust dust collector

5 結(jié)論與建議

通過對茫崖石棉礦第一選礦廠篩分車間粉塵濃度分布進(jìn)行數(shù)值模擬，并與現(xiàn)場實測粉塵濃度比較分析，得到如下結(jié)論與建議：

1)兩臺篩分設(shè)備附近粉塵濃度較大，并以篩分設(shè)備為中心隨著距離增大粉塵濃度呈現(xiàn)降低的趨勢.建議對篩分設(shè)備采取防塵降塵措施，從塵源位置遏止粉塵逸散.

2)與反彈邊界相比，捕集邊界下車間內(nèi)粉塵濃度更低，捕集粉塵顆粒效果更好.在工藝允許前提下，適當(dāng)采取地濕潤車間墻壁、地面、頂板等措施.

3)根據(jù)車間進(jìn)風(fēng)口和出風(fēng)口的現(xiàn)有布局情況，通過數(shù)值模擬得出：當(dāng)車間進(jìn)風(fēng)風(fēng)速控制在0.6 m/s左右時，能夠較好地排塵車間粉塵.建議將篩分車間進(jìn)風(fēng)風(fēng)速設(shè)為 0.6 m/s左右.

4)在兩臺篩分設(shè)備物料出口處安裝負(fù)壓抽風(fēng)集塵罩后，各測點(diǎn)的平均控塵效率達(dá)到90%，車間石棉纖維濃度控制在 2.52 f/mL以下，主抽風(fēng)管風(fēng)速建議控制在 14 m/s時為宜.

抽風(fēng)控塵研究的內(nèi)容還有很多，其中集塵罩的尺寸、形狀、布置方式等均可做進(jìn)一步深入的研究.

[1] 王瑜.我國石棉開采和加工企業(yè)粉塵危害研究[J].中國安全生產(chǎn)科學(xué)技術(shù),2011,7(4):144-147.

WANG Yu. Study on dust hazard of asbestos mining and processing enterprises in China[J].Journal of Safety Science and Technology,2011,7(4):144-147.(In Chinese)

[2] 王勇毅,姜亢,郭建中,等.石棉相關(guān)產(chǎn)品生產(chǎn)過程粉塵危害與控制對策[J].中國安全生產(chǎn)科學(xué)技術(shù),2012,8(12):156-160

WANG Yongyi,JIANG Kang,GUO Jiangzhong,etal.Dust hazards and control countermeasures in the process of asbestos-related production[J].Journal of Safety Science and Technology,2012,8(12):156-160.(In Chinese)

[3] 樊晶光.石棉粉塵控制現(xiàn)狀分析[J].勞動保護(hù),2005(1):26-27.

FAN Jingguang. The analysis of asbestos dust control[J].Labour Protection,2005(1):26-27.(In Chinese)

[4] 陳照亮,楊俊杰,金春姬.我國石棉的安全性分析及對策研究[J].中國安全生產(chǎn)科學(xué)技術(shù),2007,3(2):36-42.

CHEN Zhaoliang,YANG Junjie,JIN Chunji.Safety analysis and countermeasures of asbestos[J].Journal of Safety Science and Technology,2007,3(2):36-42.(In Chinese)

[5] 陳志霞,陳水平,張林忠,等.某石棉礦接塵工人惡性腫瘤10年回顧性調(diào)查[J].中國職業(yè)醫(yī)學(xué),2008,35(5):391-393.

CHEN Zhixia,CHEN Shuiping,ZHANG Linzhong,etal.Ten-year retrospective study on malignant tumor in workers exposed to asbestos dust in an asbestos mine[J].China Occupational Medicine,2008,35(5):391-393.(In Chinese)

[6] 甘四洋,董發(fā)勤,曾婭莉,等.溫石棉、納米SiO2、硅灰石及人造纖維粉塵的細(xì)胞毒性研究[J].安全與環(huán)境學(xué)報,2009,9(4):13-16.

GAN Siyang,DONG Faqin,ZENG Yali,etal.Research on toxicity of chrysotile,nana-silica,wollastonite and 3 kinds of artificial mineral fibers[J].Journal of Safety Science and Environment,2009,9(4):13-16.(In Chinese)

[7] 譚聰，蔣仲安，王明，等.綜放工作面多塵源粉塵擴(kuò)散規(guī)律的相似實驗[J].煤炭學(xué)報，2015，40(1)：122-127.

TAN Cong，JIANG Zhongan，WANG Ming，etal．Similarity experiment on multi-source dust diffusion law in fully mechanized caving face [J].Journal of China Coal Society,2015，40(1)：122-127.(In Chinese)

[8] 周剛，程衛(wèi)民，陳連軍，等.綜放工作面粉塵濃度空間分布規(guī)律的數(shù)值模擬及其應(yīng)用[J].煤炭學(xué)報，2010，35(12)：2094-2099.

ZHOU Gang,CHENG Weimin,CHEN Lianjun,etal.Numerical simulation and its application of dust concentration spatial distribution regularities in fully-mechanized caving face [J].Journal of China Coal Society,2010,35(12):2094-2099.(In Chinese)

[9] SILVESTER S A,LOWNDES I S,KINGMAN S W,etal.Improved ventilation and dust capture in underground crushing plants[J].Applied Mathematical Modelling,2007,31(2):311-331.

[10] CONSTANCE J A. Ways to control airborne contaminants in your plant [J].Powder and Bulk Engineering,2004,18(5):118-119.

[11] 黃優(yōu),桂祥友,夢絮屹,等.篩分車間粉塵顆粒在氣流中的受力分析及其運(yùn)動軌跡的研究[J].煤炭工程,2008(8):88-90.

HUANG You,GUI Xiangyou,MENG Xuyi,etal.Research on stressed analysis and moving trace on dust particles in air flow of screening shop[J].Coal Engineering,2008(8):88-90.(In Chinese)

[12] 荊俊山,張仁健,傅剛.風(fēng)速對生產(chǎn)車間納米粉塵濃度影響的模擬[J].中國分體技術(shù),2009,15(2):88-90.

JING Junshan,ZHANG Renjian,FU Gang.Numerical simulation of nano-sized particles dispersion inside production workshop under different wind speed[J].China Powder Science and Technology,2009,15(2):88-90.(In Chinese)

[13] 陳舉師,姜蘭,蔣仲安.邊坡鉆孔作業(yè)中粉塵分布及其影響因素的數(shù)值模擬[J].工程科學(xué)學(xué)報,2015,37(6):685-692.

CHEN Jushi,JIANG Lan,JIANG Zhongan. Numerical simulation of dust distribution and influencing factors in slope drilling[J].Chinese Journal of Engineering,2015,37(6):685-692.(In Chinese)

[14] 馬云東,羅根華,郭昭華.轉(zhuǎn)載點(diǎn)粉塵顆粒擴(kuò)散運(yùn)動的數(shù)值模擬[J].安全與環(huán)境學(xué)報,2006,6(2):16-18.

MA Yundong,LUO Genhua,GUO Zhaohua. Numerical simulation on application of diffuse regulation with powre dust in transshipping site[J].Journal of Safety and Environment,2006,6(2):16-18.(In Chinese)

[15] 譚聰,蔣仲安,陳舉師,等.綜采割煤粉塵運(yùn)移影響因素的數(shù)值模擬[J].北京科技大學(xué)學(xué)報,2014,36(6):716-721.

TAN Cong,JIANG Zhongan,CHEN Jushi,etal.Numerical simulation of influencing factors on dust movement during coal cutting at fully mechanized working faces[J].Journal of University of Science and Technology Bengjing,2014,36(6):716-721.(In Chinese)

[16] 陳舉師,蔣仲安,王 明.膠帶輸送巷道粉塵濃度分布的數(shù)值模擬及實驗研究[J].湖南大學(xué)學(xué)報：自然科學(xué)版,2015,42(6):127-134.

CHEN Jushi,JIANG Zhongan,WANG Ming.Numerical simulation and experimental research on dust concentration distribution in belt conveyer roadway[J].Journal of Hunan University :Natural Sciences,2015,42(6):127-134.(In Chinese)

[17] GBZ/192.1—2007 工業(yè)場所空氣中粉塵測定(第1部分)：總粉塵濃度[S].北京：人民衛(wèi)生出版社,2010:1-8.

GBZ/192.1—2007 Method for determination of dustin the air of workplace Part 1:Total dust concentration[S].Beijing: People’s Medical Publishing House,2010:1-8.(In Chinese)

[18] GBZ/192.5—2007 工業(yè)場所空氣中粉塵測定(第5部分)：石棉纖維濃度[S].北京：人民衛(wèi)生出版社,2010:26-35.

GBZ/192.5—2007 Method for determination of dustin the air of workplace Part 5:Asbestos fiber concentration[S].Beijing:People’s Medical Publishing House,2010:26-35.(In Chinese)

Numerical Simulation of Dust Mass Concentration Distribution inScreening Workshop of Asbestos Concentrator

JIANG Zhongan1?,DENG Quanlong1,SHI Xunxian2,CHEN Jushi1

(1.Key Laboratory of Minitry of Education for Efficient Mining and Safety of Metal Mine,University ofScience and Technology Beijing,Beijing 100083,China;2.Institute of Industrial Safety,China Academy of Safety Science & Technology,Beijing 100012,China)

In order to master the regularities of dust movement with time and space in the screening workshop of asbestos concentrator and obtain the optimized parameters of dust removal by ventilation,screening workshop of No. 1 concentrator in Mangya Asbestos Mine was considered as the research background. Based on the theory of gas-solid two-phase flow,the discrete phase model (DPM) of computational fluid mechanics was used to simulate the dust concentration in screening workshop. By comparing with field measurement of dust concentration distribution,simulation results are essentially coincident with the measured data. The results indicate that the dust concentration near the screening equipment is higher,and the value decreases as the distance from the screening equipment increases. The dust concentration is lower under trap wall condition than the reflect,and the dust trapping effect is better. The dust settling effect is better when the inlet wind speed is about 0.6 m/s. When the exhaust dust cover is installed and the wind speed of main exhaust pipes is maintained at 14 m/s,the average dust-removal efficiency can reach 90%,and the asbestos fiber concentration can be controlled below 2.52 f/mL.

asbestos screening workshop; dust concentration; migration regularity; discrete phase; numerical simulation

1674-2974(2017)12-0135-07

10.16339/j.cnki.hdxbzkb.2017.12.021

2017-02-07

國家科技支撐計劃資助項目(2015BAK40B02),National Science and Technology Support Program(2015BAK40B02)；國家自然科學(xué)基金資助項目(51574016),National Natural Science Foundation of China(51574016)；北京市自然科學(xué)基金青年科學(xué)基金資助項目(8164010),the Young Scientists Fund of the Beijing Municipal Natural Science Foundation(8164010);中央高校基本科研業(yè)務(wù)費(fèi)專項資金資助項目(FRF-BD-17-007A),Fundamental Research Funds for the Central Universities(FRF-BD-17-007A)

蔣仲安(1963-)，男，浙江諸暨人，北京科技大學(xué)教授，博士?通訊聯(lián)系人，E-mail: jza1963@263.net

X964

A