陳舉師　蔣仲安　王明

摘要：為了改善膠帶輸送巷道粉塵濃度超標的現狀，探索影響粉塵濃度分布的主要因素，根據相似原理，結合氣固兩相流的運動方程，導出了模擬膠帶輸送巷道粉塵運動的相似準則數，建立了膠帶輸送巷道相似模型，運用計算流體力學的Fluent軟件對膠帶輸送巷道相似模型粉塵濃度分布進行數值模擬，并與相似實驗數據對比分析，模擬結果與實驗數據基本吻合.研究結果表明，巷道平均風速及膠帶運行速度是影響粉塵濃度分布的兩大主要因素.巷道平均風速為0.15～0.60 m·s-1時，風速越大，粉塵濃度越低.膠帶運行速度為1～2.5 m·s-1時，運行速度越大，粉塵濃度越高.

關鍵詞：濃度分布；膠帶輸送巷道；粉塵；相似理論；數值模擬

中圖分類號：TD714.2 文獻標識碼：A

Abstract： In order to improve the current situation of dust concentration exceeding the limits in belt conveyer roadway，to explore the main influencing factors of dust concentration distribution， according to the similarity theory and the gassolid twophase flow equation motion， a similarity model device was designed for the specific situation in belt conveyer roadway. Fluent， which is a software of computational fluid dynamics， was used to analyze the dust concentration distribution in belt conveyor roadway of the similar model. Compared with similar experiment data analysis， the simulation result was consistent with the experimental data. The results show that the average wind speed in the roadway and the speed of the belt are the two main factors affecting the dust concentration distribution. When the average wind speed of the roadway is 0.15 m·s-1 ～ 0.60 m·s-1， the greater the wind speed， the lower the dust concentration. When the Tape speed is 1～ 2.5 m·s-1， the greater the speed， the higher the dust concentration.

Key words：concentration distribution； belt conveyer roadway； dust； similarity theory； numerical simulation

金屬礦山在井下開采過程中，礦巖的裝載與輸送是最為重要的環(huán)節(jié)之一，直接決定著礦山能否正常生產，影響著礦山的生產成本.礦巖在轉載及輸送過程中，均會產生大量的粉塵.金屬礦山粉塵中游離二氧化硅含量較高，長期吸入極易引發(fā)塵肺病，嚴重危害著井下工人的身體健康[1-3].

對此國內外的專家學者進行了大量的研究，一般通過控制礦石的含水率、進行塵源密閉、加強通風排塵效果、實施噴霧灑水、降低膠帶機高度勢能差及減小膠帶運行速度等措施來實現除塵目標.雖然這些措施能夠在一定程度上降低粉塵濃度，但其應用效果遠遠達不到國家衛(wèi)生標準的要求，需要在理論上和實踐上進行更加深入的研究[4-8].

因此，研究膠帶輸送巷道的粉塵濃度分布規(guī)律，確定出影響粉塵濃度分布的主要因素，對于探索適用于膠帶輸送巷道的粉塵控制新工藝、新技術、新設備、新方法具有十分重大的現實意義.

1粉塵產生機理

任何粉塵都要經過一定的傳播途徑，才能以空氣為媒介向周邊擴散.使粉塵顆粒從靜止狀態(tài)變成懸浮狀態(tài)的過程稱作“塵化”過程，弄清塵化機理，是治理粉塵的首要問題.在膠帶輸送巷道內，膠帶輸送機機頭及機尾轉載點處是主要塵源.此外，膠帶輸送機在運行過程中，由于膠帶自身的振動以及礦石與空氣的摩擦作用，也會產生少量的粉塵.根據現場調查及理論分析可知膠帶輸送巷道內粉塵塵化機理主要有：

1）剪切作用造成的塵化.礦石在機尾轉載過程中，上游工序中的礦石由具有一定高度的漏斗口自由落下至膠帶輸送機表面，在該過程中，礦石中所攜帶的礦粉受空氣的迎面阻力作用發(fā)生了剪切效應，導致空氣被卷入礦粉流中，礦粉流逐漸擴散，相互的卷吸作用使粉塵不斷地向外飛揚，并長時間懸浮在空氣中，造成膠帶機尾處粉塵飛揚.

2）誘導空氣造成的塵化.礦石在輸送過程中，由于塊度大小不均，膠帶上礦堆表面凹凸不平，當其以一定速度在空氣中運動時，會帶動礦石表面附近區(qū)域內的空氣一起流動，產生了所謂的誘導空氣.誘導空氣與礦堆表面的礦粉相互混合，將部分礦粉卷吸入空氣中，隨風流的運動擴散開來，造成膠帶輸送機表面粉塵飛揚.

3）設備運動造成的塵化.膠帶輸送機運行過程中，由于膠帶托輥的作用，膠帶自身會劇烈振動，導致細小粉塵顆粒自礦石堆中脫離出來，并與礦石堆表面附近區(qū)域內空氣發(fā)生混合，最終隨風流飄散.

4）裝入礦石造成的塵化.礦石在機頭轉載時，受自身重力作用自由下落，會排擠出與下落礦石相同體積的空氣，這些空氣會由漏斗口向上逸出.在該過程中，下落礦石中所含的礦粉顆粒將與空氣發(fā)生混合，并隨上逸空氣排出漏斗口外，造成機頭轉載點處粉塵彌漫.

2相似模型的建立

膠帶輸送巷道粉塵濃度分布的相似實驗主要研究粉塵在氣體中的運動規(guī)律以及分布情況，因此相似準則數的導出要同時考慮氣體和粉塵.由于氣流是推動粉塵運動的基本動力之一，因此，應使用氣體的運動方程和粉塵的運動方程來描述氣固兩相流動過程的全部物理現象[9-12].

2.1相似準則數的導出

根據氣固兩相流的運動方程，其中有量綱的物理量為ρg，ρp，Ug，Up，μg，dp，g，l，t和P等10個.

通過量綱分析，可導出包括斯托克斯準則數（Stk=d2pρp|Ug-Up|lμg）、均時性準則數（Ho=Ugtl）、弗勞德準則數（Fr=Uggl）、雷諾準則數（Re=lρgUgμg）、密度準則數（ρpρg）、運動準則數（UpUg）及歐拉準則數（Eu=PρgU2g）等7個相似準則數.其中，歐拉準則數Eu是非定型相似準則數，予以排除.

根據流動的單值條件，包括幾何條件、物理條件（ρg，μg，g）、氣流進出口速度分布等，可導出幾何相似準則數（δl=DL）、顆粒雷諾準則數（Rep=dpρg|Ug-Up|μg）和界面的相對粗糙度（ΔD）.

因此，本次實驗中可導出的獨立相似準則數共有9個：Stk，Ho，Fr，Re，ρpρg，UpUg，δl，Rep和ΔD.

2.2相似準則數的簡化

在進行相似模型實驗設計時，由于受到比例尺制約關系的限制，要想使實驗模型與原型之間各個相似準則數均相等是比較困難的，但在保證足夠準確度的情況下，運用近似模型法來進行相似模型實驗設計，以保持局部相似或近似是完全可以實現的.相似實驗模型與原型的流場均處于穩(wěn)定狀態(tài)，可以不考慮Ho；模型與原型中的粉塵顆粒都很細小，其重力可忽略不計，可以不考慮Fr；粉塵顆粒能在極短的時間內充分加速到氣流速度的99%，氣固兩相之間相對速度可忽略不計，則Up/Ug可近似地取1；粘性流體在流動時具有“自模性”，只要模型和原型的雷諾數Re均位于同一自模區(qū)，即可認為Re相似；由于本次實驗過程中所使用的礦石均取自于西石門鐵礦，其產生的粉塵與現場完全相同，則ρp/ρg可以滿足.最終，相似模型實驗滿足與原型相似的相似準則數只有Stk，δl，Rep和Δ/D等4個.

2.3相似模型的建立

模型以空氣為介質，其密度為ρ，粘性系數μ在對應點上與實型保持相等，顆粒的物質密度相等，即：ρp=ρ′p，ρg=ρ′g，μg=μ′g.

當ρg=ρ′g，μg=μ′g，l=2l′時，由雷諾數Re=lρgUgμg相等得出U′g=2Ug；

當ρg=ρ′g，l=2l′，U′g=2Ug時，由Stokes數Stk=d2pρp|Ug-Up|lμg相等得出d=2d′；由U′g=2Ug可以得出模擬實驗時的風速應該是實際風速的2倍，但根據公式Re=lvυ算出Re>2 300，此時流體流動處于完全紊流狀態(tài)，位于阻力平方區(qū)或第二自模區(qū).根據尼古拉茲實驗，在阻力平方區(qū)，只要滿足動力相似，也就保證了流體運動相似.所以風速不必滿足上述相似條件，只需按照幾何相似準則進行模型設計.

膠帶輸送巷道相似模型以西石門鐵礦11/96膠帶斜井為原型，按照幾何尺寸比為Cl=l∶l=2∶1建立而成.巷道原型高3.7 m，寬4 m，取巷道長度l=24 m，膠帶輸送機寬1.4 m，高1.4 m，膠帶寬1 m.對應建立的巷道模型高1.85 m，寬2 m，長度為l=12 m，膠帶輸送機模型寬0.7 m，高0.7 m，膠帶寬0.5 m.

該模型為獨頭結構，獨頭端以及頂板底板均為混凝土結構，兩側壁面為木質板結構.膠帶輸送機長5 m，布置在巷道中部，距前端4.5 m，后端2.5 m.將壓入式風機及送風管道安裝在相似實驗模型內右側，管道中心位置高1.3 m，距右壁0.2 m，實現壓入式通風，風流撞擊獨頭端后返回，沿著巷道流動，將膠帶輸送機安置在風流流場較為均勻的區(qū)域.膠帶輸送巷道相似實驗模型裝置示意圖，如圖1所示.

3數值模擬結果及分析

3.1幾何模型的建立及求解

3.1.1幾何模型的建立及網格劃分

根據相似實驗模型裝置的實際情況，結合數值模擬的需要，對膠帶輸送巷道相似模型做出以下簡化及假設：

1）巷道模型視為長方體，送風管道視為規(guī)則的圓柱體，放料漏斗視為規(guī)則的倒梯形臺；

2）平面、支架等簡化認為為平面邊界；

3）實驗用電機與變頻器由于體積較小，不考慮其對風流的影響；

4）不考慮相似實驗模型內部的漏風；

5）不考慮膠帶機運行過程中的震動.

根據上述簡化及假設條件，按照相似模型的實際尺寸進行建模，并進行網格劃分，得膠帶輸送巷道相似模型的幾何模型建立及網格劃分如圖2所示.幾何模型中巷道斷面寬2.0 m，高1.85 m，長12 m；送風管道直徑為0.4 m，出風端中心高度1.3 m，進風端中心高度0.3 m.網格的基本尺寸取0.1～0.2 m.

3.1.2邊界條件的設定及求解

通過查閱相關文獻資料，根據相似模型實驗裝置的實際情況，結合膠帶輸送巷道粉塵運動三維數學模型和Fluent數值模擬方法，并對區(qū)域網格進行自適應等調試，最終求解出膠帶輸送巷道相似模型粉塵濃度分布[13-15].

3.2風流流場分布

根據膠帶輸送巷道相似模型的實驗結果，在Fluent中進行參數設置，得出巷道平均風速為0.3 m·s-1時模型空間風流流場分布如圖3及圖4所示.

從圖中可以看出：

1）風流自送風管道口射出后，以逐漸擴大的風流射向獨頭端，風速沿管道軸線方向逐漸減小，導致到達獨頭端的風速較低；風流在獨頭的阻礙下發(fā)生反向回流，并受模型空間及射流卷吸作用的影響，形成渦流，最后沿模型走向經過膠帶輸送機排出.

2）從送風管道口到獨頭端，管道軸線上的風速從6.5 m·s-1衰減到2 m·s-1；在靠近獨頭端的區(qū)域，風速較大，且變化較快；距離管道口越近，風速越大，變化也越劇烈.氣流渦旋中心的風速較低，邊緣風速較高；巷道模型內膠帶輸送機所處區(qū)域平均風速約為0.3 m·s-1.

3）風流在流動至膠帶輸送機機尾處（距獨頭端4.5 m）基本達到穩(wěn)定，且在膠帶機所在巷道區(qū)域內始終保持穩(wěn)定.說明相似模型實驗中管道布置參數選擇合理，用來模擬膠帶輸送巷道內穩(wěn)定的風流流場可行.

3.3粉塵濃度空間分布

在離散相模型中進行粉塵源相關參數設置，令膠帶輸送機機頭及機尾各隨機產生420個粉塵顆粒，并跟蹤其運動軌跡得到膠帶運行速度為1 m·s-1，巷道平均風速為0.3 m·s-1時粉塵運動軌跡如圖5所示.通過大批量地進行粉塵顆粒跟蹤，并對其運動軌跡進行統(tǒng)計分析，得相似巷道模型內粉塵濃度分布如圖6、圖7所示.從圖中可知：

1）粉塵顆粒在放料口和膠帶表面產生后，在風流的作用下往巷道出口方向擴散，橫向隨機脈動.隨著擴散距離的增大，由于脈動量的累加，粉塵顆粒覆蓋范圍也逐漸增加，最終充滿整個巷道斷面.

2）在膠帶輸送機所在巷道區(qū)域內，隨著距獨頭端距離的增加，粉塵濃度沿程先上升至最大值，然后逐步緩慢降低；

3）在巷道斷面內，水平方向上粉塵濃度在膠帶中部達到最大值，并以膠帶中部為中心向兩側降低，膠帶左側粉塵濃度低于膠帶右側；豎直方向上粉塵濃度以上下部膠帶表面為中心向上下兩側逐步降低.

3.4不同邊界條件下粉塵濃度分布

通過現場調查及理論分析可知，巷道平均風速及膠帶運行速度是膠帶輸送巷道內影響粉塵濃度分布的兩大主要因素.因此，取人行道中央斷面（W=0.3 m）與H=1 m平面的交線作為基準線研究這兩個參數對粉塵濃度分布的影響.

3.4.1巷道平均風速

在膠帶運行速度為1 m·s-1，巷道平均風速分別為0.15 m·s-1，0.3 m·s-1，0.45 m·s-1及0.6 m·s-1條件下粉塵濃度沿程分布如圖8所示.

從圖8中可以看出：

1）在巷道平均風速為0.15～0.6 m·s-1的范圍內，風速越大，巷道空間內粉塵濃度整體越低.隨著風速的增大，粉塵濃度達到最大值的位置距獨頭端越遠.

2）較大的風速有利于粉塵顆粒的稀釋及排出，但同時也會造成已沉降粉塵的二次飛揚，在本次模擬過程中，所選擇的風速均未達到揚塵風速，風流對粉塵顆粒主要起著稀釋及排出的效果.

3）在現場應用中，在滿足規(guī)程要求的情況下，為達到一個較好的排塵效果，應將巷道風速保持在一個比較合理的數值范圍，這樣能實現風流的排塵效果，也能保證已沉降粉塵不被風流揚起.

3.4.2膠帶運行速度

在巷道平均風速為0.3 m·s-1，膠帶運行速度分別為1 m·s-1，1.5 m·s-1，2 m·s-1及2.5 m·s-1條件下粉塵濃度沿程分布如圖9所示.從圖9中可以看出：

1）膠帶運行速度對巷道內粉塵濃度分布規(guī)律的影響較弱，其變化趨勢基本保持一致.在膠帶運行速度為1～2.5 m·s-1的范圍內，隨著膠帶運行速度的增大，巷道內粉塵濃度整體升高.

2）膠帶運行速度越大，所引起的膠帶自身的振動就越強烈，且礦石與空氣的摩擦效應也更強，導致膠帶表面所產生的粉塵顆粒逃脫束縛的能力越強，模型空間粉塵濃度相應升高.

3）實際生產過程中，在保證完成生產任務的前提下，應盡量減小膠帶輸送機的運行速度，以降低膠帶輸送巷道內粉塵的濃度.

2）粉塵濃度測點布置

在巷道模型內部選取6個待測斷面，每個斷面布置3×4個粉塵濃度測點，測點布置參照風速測點布置圖.沿著膠帶輸送機運料方向在輸送機左側人行道中心線、右側送風管道中心線及輸送機中軸線上布置測點，從機尾放料口處開始，每隔0.8 m布置1個測點，直至機頭處.相似實驗模型粉塵濃度測點布置如圖11所示.

4.2風流流場分布

取巷道平均風速為0.3 m·s-1研究巷道模型內風流流場分布情況.由于篇幅有限，此處只選取H= 1.0 m平面作為代表進行風流流場分析.圖12所示為H=1.0 m平面分別與W=0.3 m，W=0.9 m，W=1.6 m斷面的交線上風速的沿程變化規(guī)律，從圖12中可看出：

在H=1.0 m平面內，W=0.9 m，W=1.6 m斷面上風速沿著回風方向先急劇下降，在距獨頭端約3.5 m處趨于穩(wěn)定；W=0.3 m斷面風速沿程先減小，后增大，隨后又迅速減小，在距獨頭端約3.5 m處也趨于穩(wěn)定.在距獨頭端3.5 m以外區(qū)域內，風速基本保持在0.32 m·s-1左右.

4.3粉塵濃度分布

取巷道平均風速為0.3 m·s-1對巷道模型空間粉塵濃度分布進行詳細測定，實驗過程中膠帶機運行速度為1.5 m·s-1.圖13、圖14分別為不同巷道斷面、不同高度平面內粉塵濃度沿程分布，從圖13、圖14中可以看出：

隨著距膠帶機放料口距離的增加，粉塵濃度沿程先逐步上升至一個最大值，后逐步緩慢降低.在水平方向上，粉塵濃度在膠帶中部分布最高，膠帶右側居中，膠帶左側最低.在高度方向上，粉塵濃度在H=1.0 m高度平面內達到最大值，并以該平面為中心向上下兩側逐步降低.

①膠帶運行速度在1.0～2.5 m·s-1范圍內，隨著膠帶運行速度的增加，粉塵濃度在人行道中央斷面與H=1.0 m平面大交線上整體分布也相應升高，且最大值出現的位置距放料口越近.

②膠帶運行速度越大，膠帶機自身的振動頻率也就越快，振動幅度隨之增大，導致產生的粉塵顆粒越多，粉塵濃度越高.

③較大的膠帶運行速度使得粉塵顆粒的產生速率大大超出了巷道風流的排塵能力，導致粉塵顆粒不能及時排出，故其最大值出現的位置距放料口要近一些.

4.5模擬結果的驗證

為了驗證模擬結果的準確性，使數值模擬與相似實驗過程中所有參數及邊界條件設置均保持一致，取膠帶運行速度為1 m·s-1，巷道平均風速為0.3 m·s-1時人行道中央斷面與H=1 m平面的交線作為基準線，將沿程風速及粉塵濃度分布作為比較對象，對模擬結果及實驗數據進行對比分析結果如圖17，圖18所示.

從圖17及圖18中可以看出，在進行比較的基準線上，風速及粉塵濃度的模擬結果與實驗數據基本吻合，分布及變化規(guī)律基本保持一致.但相比較之下兩者數值上還是略微有所偏差，其中，風速模擬結果與實驗數據相比其相對誤差在所有測點的平均值為17%，粉塵濃度模擬結果相對誤差的平均值為25%.這是由于在實驗數據測定、模型建立及參數設置過程中均會存在一定誤差所造成的.通過對比分析，說明采用離散相模型對膠帶輸送巷道相似模型粉塵運動進行模擬是可行的，模擬結果可信.

5結論

1）膠帶輸送巷道風速及粉塵濃度的模擬結果與實驗數據基本吻合，分布及變化規(guī)律基本保持一致.說明采用離散相模型對膠帶輸送巷道相似模型粉塵濃度分布及其影響因素進行模擬是可行的，模擬結果可信.

2）風流自送風管道口射出后，以逐漸擴大的風流射向獨頭端，風速沿管道軸線方向逐漸減小，在獨頭的阻礙下發(fā)生反向回流，并在膠帶輸送機所處區(qū)域保持穩(wěn)定.說明相似模型實驗中送風管道布置參數選擇合理，用來模擬膠帶輸送巷道內穩(wěn)定的風流流場可行.

3）隨著距獨頭端距離的增加，巷道區(qū)域內粉塵濃度沿程先上升至最大值，然后逐步緩慢降低；水平方向上，粉塵濃度在膠帶中部達到最大值，并以膠帶中部為中心向兩側逐步降低；豎直方向上，粉塵濃度以上下部膠帶表面為中心向上下兩側逐步降低.

4）巷道平均風速及膠帶運行速度是膠帶輸送巷道相似模型內影響粉塵濃度分布的兩大主要因素.在巷道平均風速為0.15～0.6 m·s-1的范圍內，風速越大，巷道空間內粉塵濃度整體越低.在膠帶運行速度為1～2.5 m·s-1的范圍內，隨著膠帶運行速度的增大，巷道內粉塵濃度整體升高.

參考文獻

[1]周江.西石門鐵礦皮帶運輸粉塵運動規(guī)律數值模擬及控制的研究[D].北京：北京科技大學土木環(huán)境工程學院，2010：28-31.

ZHOU Jiang.Study on numerical simulation of dust movement law and dust control of belt conveyers in Xishimen iron mine[D].Beijing：School of Civil & Environmental Cngineering. University of Science and Technology Beijing，2010：28-31.（In Chinese）

[2]朱虹，常立偉，姜歡.關于選礦廠皮帶運輸機粉塵控制的探討[J].有色礦冶，2009，25（4）：57-58.

ZHU Hong，CHANG Liwei，JIANG Huan.Dust controlling of belt conveyers in concentrating mills[J].NonFerrous Mining and Metallurgy， 2009， 25（4）： 57-58.（In Chinese）

[3]鐘新民.物流公司皮帶運輸系統(tǒng)除塵方案的可行性分析[J].甘肅冶金，2010，32（6）：154-155.

ZHONG Xinmin.Feasibility analysis of the dusting project transportation system in logistics[J].Gansu Metallurgy，2010，32（6）：154-155.（In Chinese）

[4]呂太，田志葉，李想.電廠輸煤皮帶落煤點吸塵裝置的數值模擬與實踐[J].東北電力學院學報，2006，26（2）：5-9.

LV Tai，TIAN Zhiye，LI Xiang.Numerical simulation and improvement of a dust collecting device at the coal dropping point of a power plant coaltransport belt[J].Journal of Northeast Dianli Universuty，2006，26（2）：5-9.（In Chinese）

[5]賈惠艷，馬云東.選煤廠輸煤系統(tǒng)轉載點粉塵產出控制技術[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.（In Chinese）

[6]吳應豪.巷道粉塵沉降規(guī)律與轉載點噴霧降塵系統(tǒng)研究[D].太原：太原理工大學機械工程學院，2007：65-72.

WU Yinghao.Research on tunnel dust settlement rule and reprinting spraying dust removing system[D].Taiyuan： School of Mechanical Engineering， Taiyuan University of Technology，2007：65-72.（In Chinese）

[7]張大明.輸煤巷道煤塵運移規(guī)律及治理技術研究[D].阜新：遼寧工程技術大學資源與環(huán)境工程學院，2010：88-91.

ZHANG Daming.Study on regularity of coal dust transport and dust treatment technologies in coal conveyer roadway[D].Fuxin： College of Resource and Environment Engineering， Liaoning Engineering Technology University，2010：88-91.（In Chinese）

[8]施春紅，歐盛南，金龍哲.礦井粉塵運移規(guī)律性的試驗研究[J]. 北京科技大學學報，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 and Technology Beijing，2007，29（2）：1-5.（In Chinese）

[9]賈惠艷.皮帶輸煤系統(tǒng)轉載點粉塵析出逸散規(guī)律及數值模擬研究[D].阜新：遼寧工程技術大學資源與環(huán)境工程學院，2007：35-40.

JIA Huiyan.Study on the regularity of dust precipitation and dissipation and numerical simulation of transpersite in the coal conveyer belt system[M].Fuxin：College of Resource and Environment Engineering，Liaoning Engineering Technology University，2007： 35-40.（In Chinese）

[10]馬云東，羅根華，郭昭華.轉載點粉塵顆粒擴散運動的數值模擬[J].安全與環(huán)境學報，2006，6（2）：16-18.

MA Yundong，LUO Genhua，GUO Shaohua.Numerical simulation on application of diffuse regulation with power dust in transshipping site[J].Journal of Safety and Environment，2006，6（2）：16-18.（In Chinese）

[11]劉少虹，許振良.輸煤巷道煤塵擴散的數值模擬[J].科學技術與工程，2007，7（21）：5672-5674.

LIU Shaohong，XU Zhenliang.Numerical simulation of diffusion of coal dust in haulage roadway[J].Science Technology and Engineering，2007，7（21）： 5672-5674.（In Chinese）

[12]劉毅，蔣仲安，蔡衛(wèi)，等.綜采工作面粉塵運動規(guī)律的數值模擬[J].北京科技大學學報，2007，29（4）：351-353.

LIU Yi，JIANG Zhongan，CAI Wei，et al.Numerical simulation of the dust movement rule in fullymechanized coal faces[J].Journal of University of Science and Technology Beijing，2007，29（4）：351-353.（In Chinese）

[13]王曉珍，蔣仲安，王善文，等.煤巷掘進過程中粉塵濃度分布的數值模擬[J].煤炭學報，2007，32（4）：386-390.

WANG Xiaozhen，JIANG Zhongan，WANG Shanwen， et al.Numerical simulation of distribution regularities of dust concentration during the ventilation process of coal roadway driving[J].Journal of China Coal Society，2007，32（4）：386-390.（In Chinese）

[14]牛偉，蔣仲安，劉毅.綜采工作面粉塵運動規(guī)律數值模擬及應用[J].遼寧工程技術大學學報：自然科學版，2010，29（3）：357-360.

NIU Wei，JIANG Zhongan，LIU Yi.Numerical simulation on dust movement regularities at fullymechanized coal faces and its utilization[J].Journal of Liaoning Technical University：Natural Science，2010，29（3）：357-360.（In Chinese）

[15]杜翠鳳，王輝，蔣仲安，等.長壓短抽式通風綜掘工作面粉塵分布規(guī)律的數值模擬[J].北京科技大學學報，2010，32（8）：957-962.

DU Cuifeng，WAMG Hui，JIANG Zhongan，et al. Numerical simulations of dust distribution in a fully mechanized excavation face with farpressingnearabsorption ventilation[J]. Journal of University of Science and Technology Beijing，2010，32（8）：957-962.（In Chinese）