李 勃， 董小琴

（1.湖北省電力勘測(cè)設(shè)計(jì)院有限公司，武漢 430000； 2.湖北工業(yè)大學(xué)土木建筑與環(huán)境學(xué)院，武漢 430000）

隨著國(guó)內(nèi)建筑業(yè)和基礎(chǔ)設(shè)施建設(shè)的快速發(fā)展，社會(huì)對(duì)混凝土的使用量快速增長(zhǎng)，需要開(kāi)發(fā)大量采石場(chǎng)以開(kāi)采混凝土所需骨料［1-2］. 在采石場(chǎng)骨料資源開(kāi)采過(guò)程中，爆破作業(yè)是必不可少的工藝環(huán)節(jié)［3］，爆破過(guò)程中產(chǎn)生的大量煙塵，不僅會(huì)污染環(huán)境，還會(huì)給安全生產(chǎn)和現(xiàn)場(chǎng)作業(yè)人員的身體健康帶來(lái)極大的威脅［4-7］. 研究采石場(chǎng)爆破粉塵顆粒的運(yùn)移和擴(kuò)散規(guī)律對(duì)粉塵防治工作具有重要意義，可為制定防塵、降塵措施提供科學(xué)依據(jù).

國(guó)內(nèi)外學(xué)者對(duì)粉塵運(yùn)移及擴(kuò)散規(guī)律的研究頗多且研究方法多樣，既有數(shù)值模擬和實(shí)驗(yàn)分析的方法，也有采用現(xiàn)場(chǎng)試驗(yàn)的方法. 蔣仲安等［8］采用ANSYS軟件對(duì)卸礦站粉塵濃度進(jìn)行數(shù)值模擬，模擬結(jié)果與實(shí)測(cè)結(jié)果基本吻合. 馮博等［9］基于離散相模型探究了綜采工作面粉塵逸散規(guī)律，發(fā)現(xiàn)在不同區(qū)域的采煤機(jī)具有相似的粉塵擴(kuò)散規(guī)律. 譚聰?shù)龋?0］通過(guò)建立離散相模型對(duì)不同粉塵的含水率和風(fēng)速條件進(jìn)行了模擬，結(jié)果表明最佳控塵風(fēng)速為1.5 m/s. 李雨成等［11］對(duì)金川龍首礦長(zhǎng)距離掘進(jìn)巷道進(jìn)行定點(diǎn)監(jiān)測(cè)發(fā)現(xiàn)，粉塵顆粒比例隨粒徑增大而減小. 目前對(duì)于爆破粉塵的研究主要集中在爆破炮煙分布規(guī)律、粉塵排放量等方面. 張?jiān)迄i和孔令華［12］研究了爆破炮煙的運(yùn)移規(guī)律. 陶子夜等［13］建立了爆破拆除工程粉塵擴(kuò)散模型. 王嘉瑩等［14］根據(jù)爆破粉塵的親水性和疏水性進(jìn)行粉塵分離，確定了影響粉塵濕潤(rùn)性的主要因素為粉塵真密度和粒徑. 劉堅(jiān)等［15］運(yùn)用粉塵擴(kuò)散理論進(jìn)行現(xiàn)場(chǎng)試驗(yàn)發(fā)現(xiàn)，粉塵排放量隨著炸藥用量的增加而增加，同時(shí)采用現(xiàn)場(chǎng)數(shù)據(jù)驗(yàn)證了水封爆破技術(shù)對(duì)降塵的重要作用. 針對(duì)采場(chǎng)的研究，廖賢鑫等［16］通過(guò)對(duì)某鐵礦爆破煙塵源進(jìn)行特征分析得出，爆破粉塵運(yùn)移規(guī)律受到風(fēng)流流場(chǎng)的顯著影響；陳舉師等［17］針對(duì)采場(chǎng)爆破作業(yè)粉塵過(guò)多且不易排出的現(xiàn)象，建議安裝局部通風(fēng)設(shè)施，且指出聯(lián)絡(luò)巷通風(fēng)排塵的最佳風(fēng)速為2 m/s. 針對(duì)粉塵運(yùn)移的研究，研究成果主要集中在巷道型作業(yè)面、煤礦采場(chǎng)綜采、綜放工作面粉塵濃度分布規(guī)律等方面，針對(duì)采石場(chǎng)梯段深孔爆破粉塵的研究很少. 鑒于此，本研究以楊房溝采石場(chǎng)為研究對(duì)象，建立了三維幾何模型，并通過(guò)對(duì)數(shù)學(xué)模型的求解研究了料場(chǎng)粉塵的運(yùn)移規(guī)律，最后對(duì)模擬結(jié)果進(jìn)行了可視化處理.

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

基于氣固兩相流理論，模擬粉塵顆粒在風(fēng)流作用下運(yùn)移規(guī)律的數(shù)值計(jì)算方法主要有歐拉-歐拉法和歐拉-拉格朗日兩種方法. 本研究采用歐拉-拉格朗日方法設(shè)置的離散相模型，將空氣視為背景流體，將爆破粉塵顆粒視為離散流體，對(duì)采石場(chǎng)中爆破粉塵進(jìn)行數(shù)值模擬研究，分析其在風(fēng)流作用下的運(yùn)移規(guī)律. 假定空氣為不可壓縮流體，分別建立連續(xù)方程和動(dòng)量守恒方程：

式中：u、v、w分別表示速度矢量在x、y、z方向上的分量，m/s；ρ表示密度，kg/m3；t表示時(shí)間，s；ui、uj分別表示x、y方向上的速度，m/s；xi、xj分別表示x、y方向上的坐標(biāo)，m；τij表示應(yīng)力張量；p表示壓力，Pa；g表示重力，N；Fi表示控制體平均顆粒流體阻力，N.

根據(jù)標(biāo)準(zhǔn)k-ε湍流模型，顆粒的輸運(yùn)方程如式（3）和式（4）所示.

式中：k為湍流動(dòng)能，J；Gk為湍流動(dòng)能k的產(chǎn)生項(xiàng)；μ和μt分別為層流和湍流黏性系數(shù)，Pa·s；ε為湍流動(dòng)能耗散率，m2/s3；Sk為湍流動(dòng)能源項(xiàng)；Sε為湍流動(dòng)能耗散率源項(xiàng)；C1ε、C2ε為經(jīng)驗(yàn)常數(shù)，分別取值1.44、1.92.

在笛卡爾坐標(biāo)系下，粉塵顆粒作用力平衡方程為如式（5）所示.

式中：FD(uf-up)表示顆粒單位質(zhì)量曳力，N；ρp、ρ分別表示顆粒密度和流體密度，kg/m3；up、uf分別表示顆粒速度和流體速度，m/s.FD計(jì)算方法見(jiàn)式（6）.

式中：CD為曳力系數(shù)；dp為顆粒直徑，m.

本研究采用隨機(jī)軌道模型，故顆粒軌道控制方程為：

式中：ua為氣流瞬時(shí)速度，m/s；τp為顆粒松弛時(shí)間，s.

2 三維模型的建立

2.1 現(xiàn)場(chǎng)概況

楊房溝采石場(chǎng)位于山體斜坡左側(cè)，山體為35°～45°的斜坡地形；地表多基巖裸露，以堅(jiān)硬的弱風(fēng)化巖為主，內(nèi)部巖性為花崗閃長(zhǎng)巖，有用層儲(chǔ)量為1 782.18 萬(wàn)m3；骨料物理力學(xué)性質(zhì)良好，滿足人工骨料原巖質(zhì)量要求. 圖1 為楊房溝采石場(chǎng)附近氣象站點(diǎn)每月多年平均風(fēng)速.

圖1 采石場(chǎng)附近氣象站點(diǎn)每月多年平均風(fēng)速Fig.1 Monthly multi-year average wind speed at meteorological stations near the quarry

2.2 三維模型的建立

楊房溝采石場(chǎng)現(xiàn)場(chǎng)情況較為復(fù)雜，幾何模型建立過(guò)程中難以遵照所有現(xiàn)場(chǎng)情況，需要根據(jù)該采石場(chǎng)現(xiàn)場(chǎng)實(shí)際情況和實(shí)測(cè)數(shù)據(jù)對(duì)模型及其計(jì)算區(qū)間做一定簡(jiǎn)化和假設(shè)，具體如下：①采石場(chǎng)邊坡傾斜角度嚴(yán)格按照設(shè)計(jì)圖紙開(kāi)挖；②為保證現(xiàn)場(chǎng)人員安全，維護(hù)邊坡穩(wěn)定，邊坡采用梯級(jí)臺(tái)階方式開(kāi)挖，所有臺(tái)階寬度和高度均相等；③爆破面為平整水平爆破面；④料場(chǎng)外側(cè)山體表面為平整坡面. 基于上述簡(jiǎn)化和假設(shè)，對(duì)該采石場(chǎng)進(jìn)行幾何模型的建立，如圖2所示.

圖2 采石場(chǎng)幾何模型Fig.2 Geometry model of the quarry

根據(jù)模擬需要構(gòu)建模型計(jì)算域，并利用Mesh軟件對(duì)采石場(chǎng)和計(jì)算域進(jìn)行網(wǎng)格劃分，篩選網(wǎng)格質(zhì)量和數(shù)量. 完成網(wǎng)格劃分后，將三維模型導(dǎo)入FLUENT軟件中進(jìn)行模擬參數(shù)設(shè)定和數(shù)值模擬（圖3）.

圖3 計(jì)算域及采石場(chǎng)網(wǎng)格劃分Fig.3 Computational domain and grid division of quarry

2.3 參數(shù)設(shè)定

本研究采用離散相模型對(duì)采石場(chǎng)爆破粉塵顆粒進(jìn)行仿真分析，采用壓力基求解器，運(yùn)用SIMPLE壓力速度耦合算法和標(biāo)準(zhǔn)k-ε湍流模型求解連續(xù)相流場(chǎng). 結(jié)合采石場(chǎng)爆破實(shí)際情況創(chuàng)建離散相噴射源，設(shè)置材質(zhì)、密度. 根據(jù)炸藥用量與顆粒初始速度之間的函數(shù)關(guān)系，確定爆破顆粒初始速度為7 m/s. 通過(guò)爆破現(xiàn)場(chǎng)粉塵顆粒粒徑分布測(cè)量結(jié)果得知，粒徑分布符合Rosin-Rammler分布，且最小顆粒粒徑為1×10-6m，最大顆粒粒徑為1.16×10-4m，由此設(shè)置粒徑分布及質(zhì)量流率等參數(shù). 爆破粉塵顆粒在風(fēng)流作用下的速度等參數(shù)隨時(shí)間而變化，故時(shí)間設(shè)定選擇非穩(wěn)態(tài)；入口邊界類型和出口邊界類型分別為速度入口、自由出流，其他數(shù)值模擬參數(shù)設(shè)定見(jiàn)表1.

表1 計(jì)算模型參數(shù)Tab.1 The calculation model parameters

3 爆破粉塵運(yùn)移規(guī)律的數(shù)值模擬分析

3.1 不同視角下爆破粉塵顆粒的運(yùn)移軌跡

為研究采石場(chǎng)中爆破粉塵顆粒的運(yùn)移規(guī)律，在離散相模型中進(jìn)行塵源設(shè)置，爆破面為粉塵顆粒的噴射面. 分別從采石場(chǎng)正視圖和俯視圖視角下，截取爆破發(fā)生后0、10、100、200、300 s 時(shí)粉塵顆粒速度分布圖，研究粉塵顆粒速度大小在垂直方向和水平方向的分布情況，進(jìn)而探究采石場(chǎng)爆破后粉塵顆粒的運(yùn)移規(guī)律.

從圖4可以看出，0 s時(shí)，采石場(chǎng)爆破工作還未進(jìn)行，未產(chǎn)生爆破粉塵，其正視圖和俯視圖中均未出現(xiàn)粉塵顆粒.

圖4 采石場(chǎng)爆破粉塵顆粒速度大小及速度矢量分布情況（0 s）Fig.4 Velocity magnitude and velocity vector distribution of blasting dust particles in quarry（0 s）

由圖5a可知，爆破發(fā)生10 s后，由于爆破發(fā)生時(shí)間較短，粉塵顆粒在爆炸沖擊波的作用下從爆破面噴射而出，向上運(yùn)移距離較短，粉塵顆粒未在采石場(chǎng)內(nèi)大面積擴(kuò)散. 在重力和空氣阻力作用下，距離爆破面較遠(yuǎn)的粉塵顆粒速度較小，少部分較大粒徑的粉塵顆粒速度趨近于0 m/s；距離爆破面較近的粉塵顆粒速度偏大. 圖5b 為爆破發(fā)生10 s 后采石場(chǎng)正視視角下的粉塵顆粒速度矢量分布圖，從中可以看到，受風(fēng)流影響，采石場(chǎng)北側(cè)（圖中左側(cè)）粉塵顆粒運(yùn)移方向較為復(fù)雜，而南側(cè)（圖中右側(cè)）粉塵顆粒大都趨向于向南側(cè)運(yùn)移.圖5c為爆破發(fā)生10 s后采石場(chǎng)俯視視角下的粉塵顆粒速度矢量分布圖，由圖可知，采石場(chǎng)外側(cè)粉塵顆粒在風(fēng)流作用下被直接帶走，顆粒速度趨近于風(fēng)速；采石場(chǎng)內(nèi)側(cè)有三處速度較小，分別位于采石場(chǎng)北側(cè)、南側(cè)內(nèi)部和中前側(cè).

圖5 采石場(chǎng)爆破粉塵顆粒速度大小及速度矢量分布情況（10 s）Fig.5 Velocity magnitude and velocity vector distribution of blasting dust particles in quarry（10 s）

圖6展示了采石場(chǎng)爆破發(fā)生100 s后粉塵顆粒的運(yùn)移情況. 由圖6a可知，相較于爆破發(fā)生10 s后的運(yùn)移情況，粉塵顆粒在垂直方向上的速度大小分布較為均勻，分層不明顯. 從圖6b可以看出，速度較大的粉塵顆粒均受風(fēng)流流場(chǎng)的影響顯著，運(yùn)移方向較為統(tǒng)一，均朝向采石場(chǎng)南側(cè). 從圖6c可以看出，采石場(chǎng)北側(cè)粉塵顆粒速度較小，出現(xiàn)粉塵顆粒聚集現(xiàn)象. 整體來(lái)看，離開(kāi)采石場(chǎng)的粉塵顆粒速度比采石場(chǎng)內(nèi)部的粉塵顆粒速度大.

圖6 采石場(chǎng)爆破粉塵顆粒速度大小及速度矢量分布情況（100 s）Fig.6 Velocity magnitude and velocity vector distribution of blasting dust particles in quarry（100 s）

從圖7可以看出，爆破發(fā)生200 s后，部分粉塵顆粒在重力和空氣阻力作用下發(fā)生自然沉淀，被地面及采石場(chǎng)三面臺(tái)階邊坡吸附，部分粒徑較小的粉塵顆粒由于重力較小，所以繼續(xù)向上方及四周運(yùn)移，此時(shí)粉塵顆粒已運(yùn)移至采石場(chǎng)垂直方向上一半以上區(qū)域. 圖7a顯示，采石場(chǎng)內(nèi)部大多為速度較小的粉塵顆粒. 由圖7b可以得出，爆破發(fā)生200 s后，速度較大的粉塵顆粒依然是受風(fēng)流影響較大的粉塵顆粒，其運(yùn)移方向與風(fēng)流方向緊密相關(guān). 在風(fēng)流流場(chǎng)作用下仍停留在采石場(chǎng)內(nèi)部的粉塵顆粒速度均偏小. 圖7c反映了粉塵顆粒的水平運(yùn)移情況，此時(shí)，較小粒徑的粉塵顆粒布滿整個(gè)采石場(chǎng)區(qū)域，仍存在北側(cè)粉塵顆粒聚集現(xiàn)象.

圖7 采石場(chǎng)爆破粉塵顆粒速度大小及速度矢量分布情況（200 s）Fig.7 Velocity magnitude and velocity vector distribution of blasting dust particles in quarry（200 s）

圖8為采石場(chǎng)爆破發(fā)生300 s后粉塵顆粒的運(yùn)移情況. 從圖中可以看出，相較于之前的粉塵顆粒分布情況，在風(fēng)流的夾帶作用、顆粒自身重力及空氣阻力的共同作用下，采石場(chǎng)內(nèi)部粉塵顆粒數(shù)量明顯減少. 停留在采石場(chǎng)內(nèi)部的粉塵顆粒粒徑較小，受自身重力影響較小，由于風(fēng)流流場(chǎng)作用，這些粉塵顆粒受風(fēng)流影響不顯著.

圖8 采石場(chǎng)爆破粉塵顆粒速度大小及速度矢量分布情況（300 s）Fig.8 Velocity magnitude and velocity vector distribution of blasting dust particles in quarry（300 s）

3.2 不同高度下爆破粉塵的質(zhì)量濃度分布情況

在采石場(chǎng)現(xiàn)場(chǎng)作業(yè)中，有工人（呼吸帶高度約為1.5 m）、挖掘機(jī)及運(yùn)輸卡車司機(jī)（呼吸帶高度約為2 m）等進(jìn)行石料的開(kāi)采和運(yùn)輸，因此研究不同高度下粉塵的質(zhì)量濃度分布規(guī)律，可為制定更加有效的降塵措施提供科學(xué)依據(jù). 本研究利用FLUENT后處理軟件TECPLOT，分別對(duì)z=0 m（噴射面）、z=1.5 m（工人呼吸帶高度）和z=2.0 m（挖掘機(jī)及運(yùn)輸卡車司機(jī)呼吸帶高度）的粉塵質(zhì)量濃度云圖進(jìn)行對(duì)比分析，以研究其分布規(guī)律.

圖9分別展示了采石場(chǎng)爆破發(fā)生后1、50、100、300 s時(shí)在z=0 m、z=1.5 m和z=2.0 m高度下的粉塵質(zhì)量濃度.

從圖9可以看出，爆破發(fā)生1 s后，粉塵在爆炸沖擊波作用下從噴射面上的爆破孔噴出，噴射面上粉塵顆粒密布，但此時(shí)爆破剛發(fā)生，垂直方向上的粉塵顆粒運(yùn)移距離很短，故相較于噴射面z=0 m，在z=1.5 m平面上的粉塵顆粒較少，在z=2.0 m 平面上幾乎沒(méi)有粉塵顆粒出現(xiàn). 爆破發(fā)生50 s 后，由于風(fēng)流流場(chǎng)作用，粉塵顆粒在采石場(chǎng)北側(cè)發(fā)生聚集現(xiàn)象，并在中部位置形成一條高質(zhì)量濃度分布帶，最高處的粉塵質(zhì)量濃度達(dá)到4.8×10-3kg/m3，但從z=0 m到z=2.0 m的平面上，由于粉塵顆粒粒徑不同引起的自身重力差異，不同粉塵顆粒的受力大小不同，故隨著與噴射面的距離越來(lái)越遠(yuǎn)，粉塵質(zhì)量濃度整體上出現(xiàn)逐漸下降的現(xiàn)象. 爆破發(fā)生后50 s到100 s時(shí)段內(nèi)，在自身重力影響下，部分粉塵顆粒發(fā)生自然沉降，部分粉塵顆粒在風(fēng)流的夾帶作用下逐漸向其他方向運(yùn)移，所以高質(zhì)量濃度分布帶上的粉塵質(zhì)量濃度逐漸下降，但此時(shí)，采石場(chǎng)北側(cè)出現(xiàn)粉塵顆粒聚集的現(xiàn)象，故其粉塵質(zhì)量濃度相對(duì)較高. 爆破發(fā)生100 s后，采石場(chǎng)北側(cè)越靠近噴射面，粉塵質(zhì)量濃度越高.爆破發(fā)生300 s后，高質(zhì)量濃度分布帶已完全消失，粒徑較小的粉塵顆粒較為均勻地分布在噴射面、現(xiàn)場(chǎng)工人呼吸帶高度、挖掘機(jī)及運(yùn)輸卡車司機(jī)呼吸帶高度上，在z=1.5 m高度上的粉塵質(zhì)量濃度約為6×10-5kg/m3.

4 結(jié)論

基于氣固兩相流理論，采用流體力學(xué)軟件FLUENT，建立了離散相模型，并對(duì)楊房溝采石場(chǎng)爆破粉塵的運(yùn)移規(guī)律進(jìn)行數(shù)值模擬，同時(shí)采用TECPLOT軟件對(duì)模擬結(jié)果進(jìn)行可視化處理分析. 得出主要結(jié)論如下：

1）爆破發(fā)生后，在風(fēng)流的夾帶作用下，采石場(chǎng)外側(cè)的爆破粉塵顆粒沿著風(fēng)流方向運(yùn)移，被風(fēng)流帶出采石場(chǎng)，離開(kāi)采石場(chǎng)后，粉塵顆粒運(yùn)移速度趨近于風(fēng)速.

2）由于受風(fēng)流流場(chǎng)的顯著影響，采石場(chǎng)內(nèi)部的爆破粉塵顆粒速度較小，在采石場(chǎng)北側(cè)出現(xiàn)粉塵顆粒聚集現(xiàn)象，粉塵質(zhì)量濃度較高，并在采石場(chǎng)中部位置形成一條高質(zhì)量濃度分布帶. 部分粉塵顆粒因自身重力和空氣阻力的影響發(fā)生自然沉降，部分粉塵顆粒在風(fēng)流持續(xù)作用下沿著風(fēng)流方向運(yùn)移，在爆破發(fā)生50 s后高質(zhì)量濃度分布帶上的粉塵質(zhì)量濃度逐漸下降，爆破發(fā)生300 s后，高質(zhì)量濃度分布帶已完全消失.

3）爆破發(fā)生后300 s內(nèi)，采石場(chǎng)內(nèi)不同高度的粉塵質(zhì)量濃度分布存在差異. 爆破發(fā)生1 s后，z=0 m的噴射面上粉塵的質(zhì)量濃度較大，而z=2 m的平面上幾乎沒(méi)有粉塵顆粒出現(xiàn)；爆破發(fā)生50 s后，z=0 m的噴射面上形成一條高質(zhì)量濃度分布帶，而z=2 m的平面上的粉塵顆粒分布規(guī)律與z=0 m和z=1.5 m平面上的不相同，且粉塵質(zhì)量濃度偏低. 爆破發(fā)生300 s后，停留在采石場(chǎng)內(nèi)的粉塵顆粒粒徑較小，不同高度上的粉塵質(zhì)量濃度差異也較小，粉塵顆粒分布較為均勻.