皮帶輸煤暗道通風(fēng)排塵改造方案優(yōu)化研究*

陳景序，荊德吉，葛少成，馬 恒

(遼寧工程技術(shù)大學(xué) 安全科學(xué)與工程學(xué)院，遼寧 阜新 123000)

0 引言

皮帶輸煤暗道內(nèi)部運(yùn)輸作業(yè)基本方式主要依賴皮帶輸煤系統(tǒng)，而隨著煤礦選煤或洗煤等生產(chǎn)環(huán)節(jié)產(chǎn)量的不斷增大，造成輸煤作業(yè)時暗道內(nèi)部粉塵污染嚴(yán)重[1-3]。暗道內(nèi)部粉塵濃度超過國家標(biāo)準(zhǔn)的主要原因是輸送過程中密閉導(dǎo)料槽內(nèi)煤料受運(yùn)輸過程中的碰撞和振動沖擊，及紊亂氣流所形成復(fù)雜卷吸混合作用下，內(nèi)部高速含塵氣流在出口或覆膜濾袋除塵器破口處發(fā)生逃逸現(xiàn)象。逃逸的微細(xì)粉塵長時間懸浮于密閉工作空間內(nèi)，極易誘發(fā)安全生產(chǎn)隱患和嚴(yán)重肺部疾病(CWP)[4-5]。

懸浮的粉塵由于其自身物理特性可長時間懸浮在暗道內(nèi)部空間，研究顯示懸浮粉塵的主要運(yùn)動特征，受風(fēng)流變化影響較為明顯，國內(nèi)外專家學(xué)者對此種類粉塵運(yùn)動規(guī)律展開了大量研究[6-10]。陳舉師，蔣仲安等[11]通過采用相似實(shí)驗(yàn)和數(shù)值模擬的方法，分析了多影響因子條件下煤塵在巷道內(nèi)部的分布規(guī)律。Camelli等[12]使用CFD模擬方法開發(fā)出更實(shí)際的計算流體動力學(xué)方法來驗(yàn)證地下粉塵污染物擴(kuò)散及其氣流動力學(xué)。Widiatmojo等[13]通過采用示蹤氣體結(jié)合數(shù)值仿真的方法描述了通風(fēng)系統(tǒng)擴(kuò)散特性，有利于檢測礦井通風(fēng)網(wǎng)的泄漏問題。因此，依托計算流體力學(xué)理論進(jìn)行數(shù)值模擬分析，探索皮帶輸煤暗道內(nèi)粉塵運(yùn)移規(guī)律及其濃度分布特點(diǎn)，給出合理通風(fēng)參數(shù)，為優(yōu)化輸煤暗道通風(fēng)除塵方案設(shè)計提供理論基礎(chǔ)。

1 現(xiàn)場工程概況

平朔公司木瓜界選煤廠給料機(jī)下料處264暗道輸煤，現(xiàn)有通風(fēng)方案采用一壓兩抽的通風(fēng)除塵方式。在機(jī)尾安全出口設(shè)置一臺zk25水暖型空調(diào)機(jī)組供風(fēng)制熱，機(jī)頭設(shè)置2臺T35/1.1 kW局扇進(jìn)行抽風(fēng)，來承擔(dān)輸煤暗道內(nèi)懸浮粉塵的清除工作。全負(fù)荷輸煤工作允許4臺給料機(jī)同時工作，現(xiàn)場產(chǎn)塵量較大，而人員出入安全門造成短路風(fēng)流，暗道內(nèi)部不能形成有效的排塵風(fēng)速，嚴(yán)重限制了懸浮粉塵的排除和稀釋，暗道內(nèi)部粉塵濃度嚴(yán)重超過國家標(biāo)準(zhǔn)。現(xiàn)有通風(fēng)除塵方案如圖1所示，圖1中箭頭方向代表該方案內(nèi)部風(fēng)速風(fēng)流場的運(yùn)移方向。表1為輸煤暗道工作面主要參數(shù)。

表1 輸煤暗道工作面參數(shù)

圖1 264輸煤暗道通風(fēng)除塵方案現(xiàn)狀Fig.1 264coal conveyer handling conduit ventilation dedusting program status quo

根據(jù)GBZ/T 192.1-2007《工作場所空氣中粉塵測定》分別在內(nèi)部呼吸帶高度1.6 m位置處共布置5個測點(diǎn)，在輸煤工作進(jìn)行時分別對暗道粉塵濃度、分散度進(jìn)行測量，暗道內(nèi)部測點(diǎn)布置如圖2布置。未采用改造方案密閉現(xiàn)場時懸浮粉塵分散度測定結(jié)果如表2所示。

圖2 暗道測點(diǎn)布置方案Fig.2 Underground track layout plan

結(jié)合表2實(shí)測數(shù)據(jù)顯示暗道未密閉時，現(xiàn)場粉塵污染嚴(yán)重，各測點(diǎn)數(shù)據(jù)顯示粒徑小于5%的粉塵分散度占比較大。測點(diǎn)1處小于5%的粉塵分散度比其他位置略低，是由于風(fēng)流卷吸作用下小粒徑粉塵隨風(fēng)流方向遷移，導(dǎo)致大粒徑粉塵占比增大；同理測點(diǎn)4處受尾部局扇影響，粒徑小于5%粉塵分散度有所降低，但不明顯；測點(diǎn)2,3處風(fēng)流場趨于靜止，粉塵大量集聚；測點(diǎn)5粒徑小于5%以下粉塵分散度反而上升，是由于機(jī)尾的粉塵集聚，風(fēng)流場短路而形成的。

表2 懸浮粉塵分散度測定Table.2 Determination of suspended dust dispersion

2 數(shù)值分析

2.1 數(shù)學(xué)模型的選定

為了準(zhǔn)確描述皮帶輸煤暗道中較小粒徑的粉塵受流場中拖拽力、壓力梯度、重力等影響而發(fā)生運(yùn)動軌跡的偏差。忽略粉塵運(yùn)動或碰撞而產(chǎn)生的邊緣流場變化對暗道內(nèi)氣相湍流場運(yùn)動狀態(tài)改變，故而只考慮流場的運(yùn)動狀態(tài)。用于描述運(yùn)動流體中不可壓縮流體的動量，熱和質(zhì)量傳遞過程的穩(wěn)態(tài)方程組被稱為Navier-Stokes方程。使用N-S偏微分方程通過離散化和數(shù)字化可用于分析解決氣固兩相流中顆粒相運(yùn)動。

連續(xù)性方程:

(1)

動量守恒方程:

(2)

相間曳力模型：

(3)

(4)

式中:ε=i時代表顆粒相，ε=j時代表氣相；α為質(zhì)量濃度，mg/m3；t為時間，s；τ為應(yīng)力張量；SM為動量源；ρ為密度，kg/m3；u為速度，m/s；dp為顆粒粒徑，um；μ為動力粘度，Pas；β為相間動力傳遞系數(shù)；Cd為曳力系數(shù)。

2.2 幾何模型的建立

本文主要對皮帶輸煤暗道內(nèi)通風(fēng)系統(tǒng)進(jìn)行模擬分析。為簡化風(fēng)流流場計算，對暗道內(nèi)設(shè)備模型簡化構(gòu)建，設(shè)橫截斷面底為4.2 m，壁高5.8 m，長度設(shè)為140 m的長方體，工作面上有振動給料機(jī)、皮帶、覆膜濾袋脈沖除塵器、密封帶等各種除塵設(shè)備均以標(biāo)準(zhǔn)幾何體進(jìn)行表示。采用solidworks繪圖軟件繪制皮帶輸煤暗道的三維幾何模型，并采用CFD軟件網(wǎng)格劃分程序進(jìn)行劃分，獲得的幾何圖形如圖3所示。

圖3 物理模型網(wǎng)格劃分Fig.3 Physical model meshing

2.3 邊界條件設(shè)定與模型假設(shè)

根據(jù)皮帶輸煤暗道內(nèi)現(xiàn)場情況及實(shí)測數(shù)據(jù)，給出數(shù)值模擬主要設(shè)置參數(shù)，設(shè)置參數(shù)如表3所示。

表3 顆粒元邊界條件參數(shù)設(shè)定Table.3 Particle element and boundary condition set

根據(jù)粉塵在氣相流場中的逸散規(guī)律，將氣體相流動視為不可壓縮的連續(xù)介質(zhì)流動或連續(xù)流體相穩(wěn)定流動，忽略重力影響；將逃逸的粉塵視為分布在暗道流域內(nèi)部球型粉塵顆粒無凝并、分裂或破碎，并考慮重力影響；計算過程中忽略煤塵顆粒與顆粒之間的作用力，以及粉塵對氣相流場的影響，只考慮氣相湍流夜里作用對粉塵運(yùn)動尺度的影響[14-15]；將機(jī)尾設(shè)為風(fēng)速入口，皮帶、支架等設(shè)備按照平面邊界處理，假定暗道內(nèi)部實(shí)現(xiàn)了全封閉，無漏風(fēng)問題。

2.4 數(shù)值模擬結(jié)果及分析

在實(shí)際模擬過程中根據(jù)《煤礦安全規(guī)程》及國家對運(yùn)輸機(jī)巷內(nèi)風(fēng)速相關(guān)標(biāo)準(zhǔn)規(guī)范規(guī)定，沿x軸正方向施加壓入式通風(fēng)速度矢量v，逐步增加風(fēng)速v由0.25 m/s至0.40 m/s，得出4種類通風(fēng)風(fēng)速條件下暗道空間內(nèi)粉塵濃度的沿程分布效果，如圖4所示。

圖4 輸煤暗道粉塵濃度分布Fig.4 Distribution of dust concentration in coal conveyer handling conduit

由圖4中可以看出，在暗道內(nèi)施加x正方向壓入式通風(fēng)風(fēng)速為0.25，0.30，0.35 ，0.40 m/s時，皮帶輸煤暗道內(nèi)部粉塵濃度分布規(guī)律。粉塵自各塵源位置發(fā)出后，受正壓風(fēng)速流場推動作用下沿x軸正方向運(yùn)移，粉塵運(yùn)移規(guī)律較為單一。當(dāng)排塵風(fēng)速較小時，粉塵受到的曳引作用及紊流脈動作用小于阻力等作用力，在內(nèi)部翻滾不易排出巷道，機(jī)尾漂浮的粉塵集聚在通道中部，中部粉塵濃度為最大值。由于出口條件設(shè)置為凍結(jié)，不發(fā)生反彈，在連續(xù)通風(fēng)的條件下，中后段粉塵濃度呈現(xiàn)出緩慢下降的趨勢。不同風(fēng)速條件下的粉塵濃度沿程變化均呈現(xiàn)一致的規(guī)律。從除塵效果來看當(dāng)風(fēng)速達(dá)到0.35 m/s時，最大粉塵濃度為3.2 mg/m3，降至國家標(biāo)準(zhǔn)的最大粉塵濃度4 mg/m3以下，然而隨著風(fēng)速的增加，粉塵濃度并未發(fā)生明顯下降，故選擇0.35 m/s為最優(yōu)化風(fēng)速。

3 通風(fēng)除塵方案優(yōu)化設(shè)計

為獲得較準(zhǔn)確的設(shè)備選型指導(dǎo)參數(shù)，依次對需風(fēng)量、通風(fēng)阻力分布、需熱量進(jìn)行計算?？紤]系統(tǒng)內(nèi)部設(shè)施的復(fù)雜性，本著便于實(shí)施、維護(hù)與使用經(jīng)濟(jì)性原則，設(shè)計過程中選擇可變頻風(fēng)機(jī)組對通風(fēng)性能進(jìn)行動態(tài)調(diào)節(jié)。

3.1 皮帶輸煤暗道通風(fēng)計算

3.1.1 現(xiàn)場風(fēng)機(jī)需風(fēng)量計算

依據(jù)模擬結(jié)果計算現(xiàn)場通風(fēng)需風(fēng)量，依據(jù)GB21455-2008《通風(fēng)除塵風(fēng)機(jī)參數(shù)計算標(biāo)準(zhǔn)》，264暗道優(yōu)化排塵風(fēng)速為0.35 m/s，獲得暗道最優(yōu)需風(fēng)量。

依據(jù)最優(yōu)排塵風(fēng)速對暗道需風(fēng)量計算：

Q=V×S

(5)

式中：Q為風(fēng)量，m3/s；S為掘進(jìn)巷道的斷面積，m2；V為巷道斷面風(fēng)速，m/s。

通過以上計算，求得輸煤暗道所需風(fēng)量為30 693.6 m3/h。

3.1.2 暗道通風(fēng)阻力分布

本文中輸煤暗道通風(fēng)阻力實(shí)地測量較為復(fù)雜，為了解礦井通風(fēng)系統(tǒng)的阻力分布情況通過假定暗道內(nèi)部完全密閉條件下獲得暗道通風(fēng)阻力圖，獲得的優(yōu)化排塵風(fēng)速下的264暗道阻力分布如圖5所示，可以看出264暗道最大通風(fēng)阻力為245.94 Pa。

圖5 優(yōu)化排塵風(fēng)速下的264輸煤暗道阻力分布Fig.5 The optimizes resistance distribution of ventilation dedusting on 264 coal conveyer handling conduit

3.1.3 換熱量計算

當(dāng)?shù)貧庀缶謿庀蠼y(tǒng)計數(shù)據(jù)顯示冬季最低氣溫為-26℃，為防止冬季輸煤過程中煤塊發(fā)生結(jié)凍，送風(fēng)機(jī)需要采用空調(diào)風(fēng)機(jī)實(shí)現(xiàn)冬季供暖。根據(jù)供暖設(shè)計可知風(fēng)機(jī)出風(fēng)口溫度需達(dá)40℃，計算該巷道的換熱量：

P=QρcΔt/3 600=747.68 kW

(6)

式中：Q為巷道通風(fēng)量，Q=30 693.6m3/h；ρ為空氣的密度，1.29 kg/m3；c為空氣的比熱，1.03 kJ/kg×℃；Δt為空氣的溫差，℃。

通過以上計算，根據(jù)輸煤暗道所換熱量747.68 kW，為現(xiàn)場選得最佳熱風(fēng)機(jī)型號。

3.2 通風(fēng)除塵風(fēng)機(jī)組現(xiàn)場布置方案設(shè)計

依托現(xiàn)場情況通過設(shè)立密閉風(fēng)門、降低通風(fēng)阻力來使風(fēng)機(jī)高效運(yùn)行。依托現(xiàn)場漏塵特點(diǎn)初步制定3種通風(fēng)除塵方案。

方案①：采用抽出式通風(fēng)。將存在內(nèi)部漏風(fēng)的相鄰廢棄通道封閉，將機(jī)頭房人員通道和設(shè)備通道用雙扇風(fēng)門封閉以杜絕漏風(fēng)，依據(jù)暗道需風(fēng)量和阻力選擇FBN05/7.5 kW抽出式局扇，風(fēng)機(jī)工況點(diǎn)采用變頻調(diào)節(jié)，適時調(diào)整工況。布置方案如圖6所示。

圖6 方案①Fig.6 Scheme①

方案②：采用壓入式通風(fēng)。考慮暗道需要保暖，機(jī)尾采用zk60變頻空調(diào)機(jī)組送風(fēng)，機(jī)尾暗道設(shè)置風(fēng)門，風(fēng)筒穿過風(fēng)門供風(fēng)，機(jī)頭房風(fēng)門處于常開狀態(tài)，如圖7所示。

圖7 方案②Fig.7 Scheme②

方案③：聯(lián)合式通風(fēng)?；诜桨涪谕L(fēng)動力不足的缺點(diǎn)，根據(jù)“風(fēng)機(jī)串聯(lián)，風(fēng)量不變風(fēng)壓增加”的原則，對機(jī)頭房進(jìn)行封閉，增設(shè)1臺抽出式局扇，型號為FBNO5/7.5 kW，如圖8所示。

圖8 方案③Fig.8 Scheme③

由于木瓜界輸煤暗道沒有供暖系統(tǒng)，暗道距地表較近，暗道溫度受地表溫度影響較大，冬季溫度較低情況下給料機(jī)易發(fā)生結(jié)凍。暗道主要需要依靠空調(diào)機(jī)組為作業(yè)設(shè)備提供熱量來避免設(shè)備結(jié)凍，所以該方案①排除；但由于空調(diào)機(jī)組內(nèi)部有過濾裝置和換熱裝置，隨著機(jī)器設(shè)備增大，過濾裝置和換熱裝置阻力也隨之增加，現(xiàn)有設(shè)備不足以克服暗道245.94 Pa的阻力，排除方案②?？紤]單獨(dú)使用空調(diào)機(jī)組動力不足，采用一壓一抽方式布置FBN05/7.5 kW局扇進(jìn)行補(bǔ)償除塵，布置的空調(diào)送風(fēng)機(jī)組和抽出式局扇均采用變頻調(diào)節(jié)，既能保證聯(lián)合風(fēng)機(jī)工況點(diǎn)處于合理區(qū)域，又起到節(jié)約能耗的作用?；谏鲜?，方案③為優(yōu)選通風(fēng)除塵方案。

3.3 實(shí)測數(shù)據(jù)與模擬結(jié)果對比分析

采用上述通風(fēng)除塵方案③一壓一抽聯(lián)合式通風(fēng)系統(tǒng)進(jìn)行改造，對漏風(fēng)處實(shí)施封閉處理，可使暗道內(nèi)部形成有效除塵風(fēng)流場。采用粉塵濃度測試儀在標(biāo)定測點(diǎn)進(jìn)行測量，每個測點(diǎn)的測量次數(shù)均不少于3次，通過調(diào)節(jié)變頻風(fēng)機(jī)風(fēng)速，獲得不同風(fēng)速條件下測量結(jié)果，并結(jié)合模擬結(jié)果進(jìn)行對比分析。實(shí)測粉塵濃度數(shù)據(jù)見表4。

表4 不同風(fēng)速條件下各測點(diǎn)實(shí)測粉塵濃度Table.4 Dust concentration under different wind speedconditions at each measuring point

通過表4現(xiàn)場實(shí)測各測點(diǎn)位置數(shù)據(jù)顯示，輸煤暗道內(nèi)部機(jī)尾處風(fēng)門測點(diǎn)1處，由壓入式風(fēng)機(jī)提供穩(wěn)定的正壓風(fēng)速，卷吸機(jī)尾處逃逸或懸浮空氣中粉塵向機(jī)頭處吸塵風(fēng)機(jī)運(yùn)移；然而粉塵在測點(diǎn)2，3處形成集聚，造成粉塵濃度隨風(fēng)速移動方向形成逐步上升趨勢，導(dǎo)致測點(diǎn)3處粉塵濃度上升為最大；上述方案③中機(jī)頭處設(shè)置2臺變頻抽出式局扇，在吸塵口處形成穩(wěn)定負(fù)壓場，懸浮在暗道內(nèi)部的小粒徑粉塵受風(fēng)流場曳力作用下再次發(fā)生明顯增速移動，跟隨風(fēng)流由局扇排出暗道，使暗道內(nèi)部測點(diǎn)4，5處粉塵濃度降低。

粉塵在流場中的運(yùn)移狀態(tài)呈現(xiàn)出很強(qiáng)的“跟隨性”，易被氣流輸運(yùn)而不易實(shí)現(xiàn)重力分離或沉降，在機(jī)頭機(jī)尾處形成的穩(wěn)定流場，可有效降低粉塵濃度，而在暗道中段風(fēng)速有所衰減，粉塵集聚，濃度最高。

4 結(jié)論

1)基于氣固兩相流理論模型研究暗道內(nèi)部風(fēng)速流場干涉懸浮粉塵運(yùn)動，獲得不同排塵風(fēng)速條件下粉塵濃度分布。模擬結(jié)果顯示雖然較大風(fēng)速有利于粉塵的稀釋及排出，然而隨著風(fēng)速的增大，降塵效果改善不明顯，為避免引發(fā)粉塵二次飛揚(yáng)，依據(jù)模擬結(jié)果選擇0.35 m/s時為最優(yōu)排塵風(fēng)速，為優(yōu)化暗道通風(fēng)改造方案提供了理論依據(jù)。

2)通過計算分析了暗道冬季最低需熱量，暗道內(nèi)部通風(fēng)阻力，依據(jù)計算結(jié)果對除塵方案提供設(shè)備選型基礎(chǔ)，對比分析各除塵方案，優(yōu)選設(shè)計方案③一壓一抽聯(lián)合式通風(fēng)除塵方案進(jìn)行工程改造。

3)采用一壓一抽聯(lián)合式通風(fēng)進(jìn)行通風(fēng)改造后，獲得了良好的降塵效果，滿足除塵技術(shù)要求，該設(shè)計即保證除塵性能，又能起到節(jié)約能耗的作用，對治理連續(xù)輸煤作業(yè)過程中的粉塵污染具有重要實(shí)踐指導(dǎo)意義。然而單純使用通風(fēng)除塵并不能從根本解決粉塵污染問題，還應(yīng)該深入現(xiàn)場粉塵污染源頭，分析粉塵污染特征，為治理粉塵污染提出更合理的改進(jìn)意見。

[1] 郭雁飛, 周建新, 陳衛(wèi)紅. 粉塵危害事件調(diào)查分析[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2009, 5(4):80-84.

GUO Yanfei， ZHOU Jianxin， CHEN Weihong. Investigation and analysis of dust hazard incidents[J]. Journal of Safety Science and Technology, 2009, 5(4): 80-84.

[2] 孫百存, 甘小根. 帶式輸送機(jī)轉(zhuǎn)載裝置的研究與應(yīng)用[J]. 煤礦機(jī)械, 2009 (5):120-121.

SUN Baicun, GAN Xiaogen. Adhesive tape conveyer reprint installment research and application [J]. Coal Mine Machinery, 2009(5): 120-121.

[3] 賈惠艷, 馬云東. 選煤廠輸煤系統(tǒng)轉(zhuǎn)載點(diǎn)粉塵產(chǎn)出控制技術(shù)[J]. 環(huán)境污染與防治, 2007, 29(10):767-770.

JIA Huiyan, MA Yundong. Generation control techniques of coal dust of transferring points of the coal transportation system in the coal selection plant [J]. Environmental Pollution and Control, 2007, 29(10): 767-770.

[4] 葛少成, 邵良杉, 齊慶杰. 選煤廠轉(zhuǎn)運(yùn)點(diǎn)除塵方案模擬優(yōu)化設(shè)計[J]. 遼寧工程技術(shù)大學(xué)學(xué)報, 2007, 26(6):805-808.

GE Shaocheng, SHAO Liangshan, QI Qingjie. Simulation and optimal design of dust removal scheme of transshipment point in a coal washery [J]. Journal of Liaoning Technical University, 2007, 26 (6): 805-808.

[5] 施春紅,歐盛南,金龍哲．礦井粉塵運(yùn)移規(guī)律性的試驗(yàn)研究[J]. 北京科技大學(xué)學(xué)報, 2007, 29(2): 1-5.

SHI Chunhong, OU Shengnan, JIN Longzhe. A study and analysis on the law of motion of the coal dust [J]. Journal of University of Science Technology Beijing, 2007, 29(2): 1-5.

[6] 張茂宇, 賈福音. 煤礦井下粉塵監(jiān)測、防治技術(shù)的現(xiàn)狀與趨勢[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2011, 7 (12):188-191.

ZHANG Maoyu, JIA Fuyin. Situation and trend on dust monitoring & control technologies in underground mine [J]. Journal of Safety Science and Technology, 2011, 7(12): 188-191.

[7] 王曉珍. 煤巷掘進(jìn)過程中粉塵濃度影響因素分析[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2011, 7(4): 75-79.

WANG Xiaozhen. Analysis of dust concentration influence factor in coal roadway driving [J]. Journal of Safety Science and Technology, 2011, 7(4): 75-79.

[8] 李雨成, 李智, 高倫. 基于風(fēng)流及粉塵分布規(guī)律的機(jī)掘工作面風(fēng)筒布置[J]. 煤炭學(xué)報, 2014, 39(S1):130-135.

LI Yucheng, LI Zhi, GAO Lun. Arrangement of air duct in tunneling working face based on the distribution laws of airflow and dust [J]. Journal of China Coal Society, 2014, 39(S1): 130-135.

[9] 張建卓, 祝天姿, 高猛, 等. 綜掘面風(fēng)幕集塵除塵系統(tǒng)的參數(shù)優(yōu)化與數(shù)值模擬[J]. 中國安全科學(xué)學(xué)報, 2011, 21(4): 128-131.

ZHANG Jianzhuo, ZHU Tianzi, GAO Meng, et al. Parameter optimization and numerical simulation of dust-collecting and dedusting system with air curtain in fully mechanized excavation face [J]. China Safety Science Journal, 2011,21(4): 128-131.

[10] Yuqing FENG, TIM SWENSERSMITH, PETER J WITT, et al. CFD modeling of gas-solid flow in an internally circulating fluidized bed [J]. Powder Technology, 2012, 219: 78-85.

[11] 陳舉師, 蔣仲安, 姜蘭. 膠帶輸送巷道粉塵分布及其影響因素的實(shí)驗(yàn)研究[J]. 煤炭學(xué)報, 2014, 39(1): 135-140.

CHEN Jushi, JIANG Zhongan, JIANG Lan. Experimental research on dust distribution and its influencing factors in belt conveyer roadway[J]. Journal of China Coal Society, 2014, 39(1): 135-140.

[12] CAMELLI F E,BYRNE G, LOHNER R.Modeling subway air flow using CFD[J]. Tunneling and Underground Space Technology, 2014,(43): 20-31.

[13] WIDIATMOJO A, SASAKI K, SUGAI Y. Assessment of air dispersion characteristic in underground mine ventilation: Field measurement and numerical evaluation[J]. Process Safety& Environmental Protection, 2015,93(13): 173-181.

[14] 荊德吉, 葛少成. 閉環(huán)回旋控塵技術(shù)的三維數(shù)值模擬研究[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2015, 11(3): 79-84.

JING Deji, GE Shaocheng. Study on 3D numerical simulation of closed-loop dust control technology [J]. Journal of Safety Science and Technology, 2015, 11(3): 79-84.

[15] 張大明, 馬云東. 輸煤巷道通風(fēng)排塵特征參數(shù)數(shù)值模擬優(yōu)化研究[J]. 安全與環(huán)境學(xué)報, 2012, 2(1): 178-181.

ZHANG Daming, MA Yundong. New numerical simulation for optimizing the characteristic is parameters of ventilation dedusting in the coal conveyer roadway [J]. Journal of Safety and Environment, 2012, 12(1): 178-181.