孫三祥，王 文，郭慧杰，雷鵬帥,2，張 健,3

(1.蘭州交通大學(xué) 環(huán)境與市政工程學(xué)院，甘肅 蘭州 730070；2.蘭州交通大學(xué) 寒旱地區(qū)水資源綜合利用教育部工程研究中心，甘肅 蘭州 730070；3.蘭州交通大學(xué) 甘肅省軌道交通力學(xué)應(yīng)用工程實(shí)驗(yàn)室，甘肅 蘭州 730070)

現(xiàn)場(chǎng)實(shí)測(cè)結(jié)果表明，高海拔隧道鉆爆法施工過(guò)程中，出碴階段CO排放量及濃度最高[1-2]。王耀等[1]采用CO檢測(cè)儀，測(cè)試分析了川藏公路雀兒山隧道(海拔4 000 m)爆破作業(yè)、內(nèi)燃機(jī)械的CO排放量以及隧道內(nèi)CO濃度變化趨勢(shì)，得出在出碴階段內(nèi)燃機(jī)械累計(jì)排放的CO總量多于爆破階段，施工機(jī)械，尤其是出蹅車(chē)排放的CO在出碴階段對(duì)隧道內(nèi)環(huán)境造成了嚴(yán)重污染。張玉偉等[2]采用便攜式CO檢測(cè)儀，跟蹤監(jiān)測(cè)了康定瓦九(瓦澤至九龍)公路雞丑山隧道(海拔4200m)不同施工階段CO濃度，得出在出碴階段工作區(qū)CO濃度最高。

國(guó)內(nèi)外專(zhuān)家對(duì)隧道施工期的通風(fēng)進(jìn)行了研究[3-5]，并對(duì)平原地區(qū)隧道爆破后CO的擴(kuò)散進(jìn)行了模擬研究。鄧祥輝等[6]采用RNGκ-ε湍流模型，基于ADINA軟件中的CFD模塊，對(duì)隧道內(nèi)爆破后的CO擴(kuò)散規(guī)律以及工作人員進(jìn)洞安全時(shí)間進(jìn)行了數(shù)值模擬；李孜軍等[7]采用RNGκ-ε湍流模型，基于Fluent軟件模擬分析了雙洞隧道獨(dú)頭掘進(jìn)爆破作業(yè)后CO擴(kuò)散規(guī)律；王曉玲等[8]采用高雷諾數(shù)κ-ε數(shù)學(xué)模型，基于STAR-CD軟件，模擬分析了爆破后不同通風(fēng)時(shí)刻隧洞內(nèi)CO遷移和分布規(guī)律。但高海拔地區(qū)，由于氣壓、空氣密度和含氧量均較低，內(nèi)燃機(jī)燃燒不充分，因此CO擴(kuò)散規(guī)律有別于平原地區(qū)。

針對(duì)高海拔隧道爆破后CO濃度擴(kuò)散的模擬，曹正卯等[9]依托新關(guān)角隧道(海拔3 300 m以上)，采用κ-ε湍流非穩(wěn)態(tài)模型，基于Fluent軟件，模擬分析不同海拔高度施工爆破后有害氣體濃度分布，得出了同一測(cè)點(diǎn)有害氣體質(zhì)量濃度與海拔高度的關(guān)系式，但未考慮施工機(jī)械及出蹅車(chē)的影響。

針對(duì)隧道出碴階段CO濃度擴(kuò)散的模擬，劉敦文等[10]采用RNGκ-ε湍流模型，基于Fluent軟件，應(yīng)用組分傳輸模型，得出了平原地區(qū)隧道出碴過(guò)程CO濃度場(chǎng)，提出了CO濃度限值及安全連續(xù)作業(yè)時(shí)間。專(zhuān)家還就隧道施工過(guò)程中車(chē)輛對(duì)隧道速度場(chǎng)及通風(fēng)效果的影響，且大多采用二維動(dòng)網(wǎng)格模型模擬[11-12]，而對(duì)高海拔隧道出碴車(chē)運(yùn)行過(guò)程中排放CO濃度擴(kuò)散的模擬目前未見(jiàn)報(bào)道。

因此，本文選取RNGκ-ε湍流模型，利用動(dòng)網(wǎng)格及用戶(hù)自定義函數(shù)，選用彈簧近似光滑模型、局部重劃模型計(jì)算網(wǎng)格的動(dòng)態(tài)變化；考慮海拔對(duì)氣壓和空氣密度的修正；由通用有限速率模型反映CO與空氣的耦合；僅考慮出碴車(chē)排放的CO，其他機(jī)械排放的CO作為固定源；基于Fluent軟件，進(jìn)行通風(fēng)工況下出碴車(chē)運(yùn)行過(guò)程中高海拔隧道內(nèi)氣流速度場(chǎng)和壓力場(chǎng)、CO濃度場(chǎng)的三維非穩(wěn)態(tài)數(shù)值模擬?；诔霾赀^(guò)程現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)驗(yàn)證模擬結(jié)果。根據(jù)模擬結(jié)果分析出碴車(chē)排放CO擴(kuò)散規(guī)律及行駛間隔時(shí)間，為高海拔隧道施工通風(fēng)方案的制訂提供一定的理論支持。

1 模型建立

1.1 隧道三維模型

以在建的采用鉆爆法施工的敦格鐵路當(dāng)金山單洞單線(xiàn)特長(zhǎng)隧道(平均海拔3 000 m)為例建立隧道三維有限元模型和網(wǎng)格簡(jiǎn)圖，取正洞長(zhǎng)度為1 600 m，平導(dǎo)長(zhǎng)度為2 000 m，動(dòng)網(wǎng)格計(jì)算長(zhǎng)度為100 m。隧道洞身坡度為12.3‰，單面上坡。出碴車(chē)簡(jiǎn)化為長(zhǎng)×寬×高為8.55 m×2.49 m×3.45 m的六面體，底面距地面高0.6 m，出碴車(chē)尾氣管口位于車(chē)體左壁，距車(chē)頭3 m，距地面0.8 m，出口直徑為0.08 m。隧道風(fēng)筒直徑為1.2 m，位于隧道側(cè)壁中上部，正洞風(fēng)筒距掌子面30 m，壓入式通風(fēng)。采用非結(jié)構(gòu)化四面體網(wǎng)格對(duì)流場(chǎng)區(qū)域進(jìn)行網(wǎng)格劃分，對(duì)出碴車(chē)車(chē)體周?chē)败?chē)頭、車(chē)位部網(wǎng)格進(jìn)行加密，如圖1所示。通過(guò)網(wǎng)格敏感性分析，確定最終節(jié)點(diǎn)數(shù)62 905個(gè)，網(wǎng)格數(shù)337 448個(gè)。出碴車(chē)行駛及CO排放對(duì)應(yīng)流場(chǎng)采用動(dòng)網(wǎng)格劃分。

1.2 動(dòng)網(wǎng)格模型的建立

將隧道內(nèi)氣流場(chǎng)分為出碴車(chē)運(yùn)動(dòng)區(qū)域和不運(yùn)動(dòng)區(qū)域，基于Fluent軟件提供的DEFINE_CG_MOTION和DEFINE_PROFILE宏，應(yīng)用C語(yǔ)言編譯用戶(hù)自定義函數(shù)，建立出蹅車(chē)運(yùn)動(dòng)和CO排放速度變化的動(dòng)網(wǎng)格模型。

1.3 氣流流動(dòng)與CO擴(kuò)散的方程

1.3.1 基本假設(shè)

(1)通風(fēng)氣流為不可壓縮流體。

(2)氣體在稀釋擴(kuò)散過(guò)程中不發(fā)生化學(xué)反應(yīng)及相變反應(yīng)。

(3)隧道壁面為等溫壁面。

(4)僅考慮出碴車(chē)排放CO的擴(kuò)散，將其他機(jī)械排放的CO統(tǒng)一作為固定源，即將出蹅車(chē)通過(guò)監(jiān)測(cè)點(diǎn)前的CO濃度值作為模擬計(jì)算中的背景值。

1.3.2 基本方程

選取RNGκ-ε湍流模型，建立氣流流動(dòng)與CO擴(kuò)散的方程，包括基本控制方程(質(zhì)量守恒方程、動(dòng)量方程、能量方程)、湍動(dòng)能方程、湍動(dòng)能耗散率方程和組分輸運(yùn)方程，分別如下。

(1)基本控制方程通式為

(1)

式中：t為時(shí)間，s；ρ為流體密度，kg·m-3；φ為待求通用物理量；U為速度矢量；Γφ和Sφ分別為對(duì)應(yīng)變量的廣義擴(kuò)散系數(shù)和廣義源項(xiàng)。

(2)湍動(dòng)能κ方程為

(2)

其中

式中：ui和uj為速度分量，m·s-1；κ為紊流動(dòng)能，m2·s-2；xi和xj為坐標(biāo)；μt為紊流動(dòng)力黏性系數(shù)，Pa·s；σκ為湍流普朗特?cái)?shù)，σκ=1.00；ε為紊流的動(dòng)能耗散率，m2·s-3；μeff為有效黏性系數(shù)；μ為層流動(dòng)力黏性系數(shù)，Pa·s；cμ為經(jīng)驗(yàn)常數(shù)，取值0.09；G為紊流脈動(dòng)動(dòng)能產(chǎn)生項(xiàng)。

(3)湍動(dòng)能耗散率ε方程為

(3)

式中：σε，c1，c2均為經(jīng)驗(yàn)常數(shù)，σε=1.30，c1=1.44，c2=1.92。

(4)組分輸運(yùn)方程為

(4)

1.3.3 邊界條件

(1)進(jìn)口邊界：出碴車(chē)排氣管為CO進(jìn)口邊界?；诂F(xiàn)場(chǎng)測(cè)試結(jié)果，確定尾氣管出口CO濃度最大值為427 mg·m-3，對(duì)應(yīng)尾氣溫度為600 K，尾氣設(shè)置為空氣與CO的混合體，其中CO在尾氣中的質(zhì)量百分比取0.042 7%。出碴車(chē)加速、勻速運(yùn)行階段CO的排放速度分別為30和23 m·s-1?，F(xiàn)場(chǎng)實(shí)測(cè)隧道內(nèi)風(fēng)速為0.3 m·s-1，該風(fēng)速設(shè)為隧道模型入口的風(fēng)速。

(2)出口邊界: 設(shè)為自由出流邊界條件。

(3)壁面邊界: 出碴車(chē)及隧道壁面設(shè)為壁面邊界，出蹅車(chē)為滑移壁面，隧道壁面為無(wú)滑移絕熱邊界。

(4)正洞工作面需風(fēng)量為1 000 m3·min-1，對(duì)應(yīng)正洞風(fēng)管出口風(fēng)速為14.8 m·s-1，入口輸送新鮮空氣(即CO質(zhì)量分?jǐn)?shù)為0),空氣溫度為293 K。

2 數(shù)值模擬算法及算法驗(yàn)證

2.1 數(shù)值模擬算法

采用非定常隱式解法求解上述方程。其中速度迭代采用PISO算法；壓力采用標(biāo)準(zhǔn)離散方式；其他參數(shù)的離散采用二階迎風(fēng)格式。網(wǎng)格的動(dòng)態(tài)變化過(guò)程選用彈簧近似光滑模型和局部重劃模型模擬。時(shí)間步長(zhǎng)設(shè)置為0.01 s。

2.2 算法驗(yàn)證

采用現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)對(duì)模擬結(jié)果進(jìn)行驗(yàn)證?，F(xiàn)場(chǎng)實(shí)測(cè)期間，隧道內(nèi)氣溫為13.8～18.5 ℃，大氣壓為70.7～71.5 kPa。出碴車(chē)為陜汽德隆F3000型自卸車(chē)，自4號(hào)橫通道運(yùn)行到3號(hào)橫通道。分別在距棧橋50，100，150 m的斷面布置測(cè)點(diǎn)，并分別標(biāo)記為測(cè)點(diǎn)1、測(cè)點(diǎn)2、測(cè)點(diǎn)3。測(cè)點(diǎn)位于出碴車(chē)右側(cè)，且距右側(cè)隧道壁面1.2 m，距隧道底面1.5 m。使用CD4(B)型多參數(shù)氣體測(cè)定器檢測(cè)CO濃度，數(shù)據(jù)采集時(shí)間間隔為15 s，持續(xù)監(jiān)測(cè)300 s。

測(cè)點(diǎn)1不同時(shí)間間隔CO濃度的測(cè)定值與數(shù)值模擬值最大差值為0.12 mg·m-3，相對(duì)誤差小于5%，可見(jiàn)兩者吻合很好，由此表明，本文提出的模擬算法可以有效地揭示出碴車(chē)行駛過(guò)程中排放CO擴(kuò)散規(guī)律。

3 數(shù)值模擬結(jié)果

3.1 出碴車(chē)運(yùn)行速度曲線(xiàn)

根據(jù)現(xiàn)場(chǎng)監(jiān)測(cè)，出碴車(chē)運(yùn)行速度曲線(xiàn)取3種，如圖2所示，即出碴車(chē)的加速時(shí)間均取4 s, 加速度a取1.00，1.25，1.50 m·s-2，對(duì)應(yīng)的勻速運(yùn)行速度v取4，5，6 m·s-1。

圖2 出碴車(chē)運(yùn)行速度曲線(xiàn)

3.2 出碴車(chē)環(huán)隙流特性

出渣車(chē)分別以圖2中的3條速度曲線(xiàn)行駛至8 s時(shí)，車(chē)頭環(huán)隙流速度分布如圖3所示。其中，以加速度為1.50 m·s-2，對(duì)應(yīng)速度為6 m·s-1的速度曲線(xiàn)行駛至8 s時(shí)，出碴車(chē)不同部位環(huán)隙流速度分布如圖4所示。圖中縱坐標(biāo)環(huán)隙速度方向與車(chē)體同向?yàn)檎聪驗(yàn)樨?fù)；橫坐標(biāo)以車(chē)體軸線(xiàn)為基準(zhǔn)，運(yùn)行方向右側(cè)至隧道壁面的距離為正，左側(cè)為負(fù)。由圖3和圖4可知：出碴車(chē)行駛速度越大，環(huán)隙流速度越大；出碴車(chē)行駛速度一定時(shí)，越靠近出碴車(chē)尾部，環(huán)隙流速度越大。

圖3 出碴車(chē)以3種速度曲線(xiàn)行駛時(shí)環(huán)隙流速度分布

圖4 出碴車(chē)不同部位環(huán)隙流速度分布

3.3 出碴過(guò)程中隧道內(nèi)氣流的壓力場(chǎng)和速度場(chǎng)

考慮出渣車(chē)基本沿隧道中線(xiàn)行駛，以距隧道底面1.5 m高度的隧道中線(xiàn)為例，區(qū)分相同車(chē)速不同時(shí)刻、相同時(shí)刻不同車(chē)速2種情況，根據(jù)模擬結(jié)果繪制不同時(shí)刻氣流的壓力和速度的沿程(沿隧道縱向距離)分布圖。

1)相同車(chē)速，不同時(shí)刻

以圖2中出碴車(chē)的加速度1.50 m·s-2，勻速運(yùn)行速度6 m·s-1的速度曲線(xiàn)為例，不同時(shí)刻氣流的壓力和速度的沿程分布曲線(xiàn)如圖5和圖6所示。

圖5 不同時(shí)刻氣流壓力的沿程分布曲線(xiàn)

由圖5可知：出碴車(chē)自啟動(dòng)(t=0)到加速結(jié)束時(shí)(t=4 s)，車(chē)頭前端氣流壓力由0增大為90 Pa，車(chē)尾增大為36 Pa，車(chē)頭為正壓，車(chē)尾為負(fù)壓，且氣流壓力均達(dá)到最大值；勻速行駛時(shí)，車(chē)頭前端區(qū)間氣流壓力由0增大為18 Pa，車(chē)尾壓力由0增大為10 Pa；加速結(jié)束時(shí)車(chē)頭氣流壓力最大值是勻速行駛時(shí)最大值的5倍，車(chē)尾氣流壓力最大值是勻速行駛時(shí)最大值的3.6倍。

圖6 不同時(shí)刻氣流速度的沿程分布曲線(xiàn)

由圖6可知：出碴車(chē)自啟動(dòng)(t=0)到加速結(jié)束時(shí)(t=4 s)，車(chē)體周?chē)鷼饬魉俣茸兓沧畲螅俣茸畲笾党霈F(xiàn)在車(chē)尾；車(chē)頭前端速度變化影響區(qū)域長(zhǎng)度約5 m，車(chē)尾渦流區(qū)影響長(zhǎng)度較大，且隨時(shí)間的延續(xù)影響區(qū)域長(zhǎng)度逐漸增大；出碴車(chē)加速階段車(chē)體不同部位氣流速度均遠(yuǎn)大于勻速行駛階段，其中車(chē)尾渦流區(qū)速度最大，為勻速運(yùn)行階段車(chē)尾渦流區(qū)最大速度的2～3倍。

2)相同時(shí)刻，不同車(chē)速

以出碴車(chē)穩(wěn)定運(yùn)行12 s時(shí)刻為例，不同車(chē)速時(shí)氣流壓力和速度的沿程分布曲線(xiàn)如圖7和圖8所示。

圖7 不同車(chē)速時(shí)氣流壓力沿程分布曲線(xiàn)

圖8 不同車(chē)速時(shí)氣流速度沿程分布曲線(xiàn)

由圖7可知：車(chē)頭前端氣流正壓影響區(qū)域長(zhǎng)度約5 m；遠(yuǎn)小于車(chē)尾負(fù)壓區(qū)域長(zhǎng)度；隨著車(chē)速增大，車(chē)頭前端壓力及車(chē)尾負(fù)壓均增大，即車(chē)頭對(duì)前面空氣的擠壓、車(chē)尾對(duì)后部空氣的抽吸越強(qiáng)，同時(shí)壓力變化影響區(qū)域也增大；車(chē)頭前端5 m以遠(yuǎn)區(qū)域壓力沿程變化趨于穩(wěn)定，不同車(chē)速形成的壓差基本相同，約為1.2 Pa。

由圖8可知：車(chē)頭前端氣流速度變化影響區(qū)域長(zhǎng)度約5 m；車(chē)尾渦流區(qū)影響區(qū)域長(zhǎng)度較大，且隨時(shí)間的延續(xù)影響長(zhǎng)度逐漸增大；車(chē)速越大，車(chē)尾渦流速度越大，影響區(qū)域長(zhǎng)度越大；運(yùn)行歷時(shí)12 s后，3種車(chē)速下的氣流速度場(chǎng)均趨于穩(wěn)定。

3.4 出碴車(chē)行駛時(shí)CO濃度場(chǎng)

鑒于加速階段歷時(shí)短，重點(diǎn)分析出渣車(chē)勻速行駛階段的CO濃度場(chǎng)。出渣車(chē)以圖2中加速度1.50 m·s-2，勻速運(yùn)行速度6 m·s-1的速度曲線(xiàn)行駛至8 s時(shí)，尾氣管出口速度矢量圖如圖9所示。由圖9可以看出：尾氣管出口下方形成渦流區(qū)；受車(chē)速和環(huán)隙流的影響，出碴車(chē)尾部氣流速度變大，紊動(dòng)作用加強(qiáng)，從而增強(qiáng)了CO的擴(kuò)散。

圖9 尾氣管出口速度矢量圖

由于尾氣管安裝在車(chē)體左側(cè)，車(chē)體行駛過(guò)程中尾氣排放應(yīng)屬于浮射流范疇，其沖擊浮射流基本特性與有限空間受橫流影響的沖擊射流有相似之處，二者間的區(qū)別另文討論。

以圖2中加速度1.50 m·s-2，勻速運(yùn)行速度6 m·s-1的速度曲線(xiàn)為例，出渣車(chē)行駛至8 s時(shí)，距尾氣管不同距離斷面CO濃度分布云圖如圖10所示。由圖10可知： CO主要集中在車(chē)尾渦流區(qū)；出碴車(chē)勻速行駛時(shí)，距尾氣管出口越近，CO濃度越高，但超標(biāo)(高于標(biāo)準(zhǔn)TB10204—2002《鐵路隧道施工規(guī)范》中的限值30 mg·m-3)區(qū)域極小；單車(chē)勻速行駛時(shí)CO濃度超標(biāo)，但超標(biāo)區(qū)域的高度未達(dá)到人體呼吸高度。

圖10 距尾氣管不同距離斷面CO濃度分布云圖

出渣車(chē)以圖2中加速度1.50 m·s-2，勻速運(yùn)行速度6 m·s-1的速度曲線(xiàn)行駛時(shí)，不同時(shí)刻、不同車(chē)速時(shí)隧道內(nèi)CO濃度分布云圖如圖11所示。出渣車(chē)以圖2中的3條速度曲線(xiàn)行駛至8 s時(shí)，CO濃度分布云圖如圖12所示。由圖11和圖12可得如下結(jié)論。

圖11 不同時(shí)刻隧道內(nèi)CO濃度分布云圖

圖12 不同車(chē)速時(shí)隧道內(nèi)CO濃度分布云圖

(1)CO自車(chē)體左側(cè)的尾氣管排出后，車(chē)體左側(cè)尾部渦流區(qū)的CO濃度較高。

(2)出碴車(chē)行駛速度越大，尾部渦流影響區(qū)也越長(zhǎng)，CO擴(kuò)散作用越強(qiáng)。因此建議出碴車(chē)車(chē)速根據(jù)現(xiàn)場(chǎng)情況盡可能以較大車(chē)速勻速行駛。

(3)出碴車(chē)加速運(yùn)行區(qū)段，CO累積嚴(yán)重，車(chē)體尾部濃度高，濃度超過(guò)標(biāo)準(zhǔn)限值30 mg·m-3。因此建議出碴車(chē)選擇出碴路線(xiàn)時(shí)，盡量減少作業(yè)點(diǎn)，從而減少出碴車(chē)加減速的次數(shù)，以減少CO排放量；作業(yè)人員選擇避讓點(diǎn)時(shí)，應(yīng)盡量選擇在出碴車(chē)勻速行駛區(qū)段、遠(yuǎn)離車(chē)輛尾氣管的一側(cè)，避免選擇在出碴車(chē)加速區(qū)段(例如在通過(guò)棧橋的區(qū)段)。

(4)出蹅車(chē)經(jīng)過(guò)測(cè)點(diǎn)較長(zhǎng)一段時(shí)間后，由于環(huán)隙流的作用，隧道頂部的高濃度CO向隧道底部運(yùn)移，致使底部CO濃度升高。根據(jù)模擬結(jié)果，車(chē)速為4，5，6 m·s-1時(shí)，隧道內(nèi)前、后2輛出碴車(chē)輛間隔時(shí)間大于5 min時(shí)，CO累積濃度不超過(guò)標(biāo)準(zhǔn)限值30 mg·m-3。因此，建議前、后2輛出碴車(chē)間隔時(shí)間盡量大于5 min。

3.5 出碴車(chē)尾部CO擴(kuò)散規(guī)律

基于出碴車(chē)尾部CO濃度的分布云圖，選取距車(chē)體中心線(xiàn)左側(cè)1.2 m，距隧道底面1.5 m的點(diǎn)為觀(guān)測(cè)點(diǎn)(該點(diǎn)斷面濃度最大)，出渣車(chē)分別以圖2中的3條速度曲線(xiàn)行駛至8 s時(shí)， CO濃度沿程變化曲線(xiàn)如圖13所示。由圖13可知，該曲線(xiàn)反映的擴(kuò)散規(guī)律符合瞬時(shí)點(diǎn)源一維隨流擴(kuò)散方程基本解(不同車(chē)速擬合優(yōu)度R2均大于0.8)，即按照e指數(shù)衰減。這也與文獻(xiàn)[7]基于實(shí)測(cè)結(jié)果所得結(jié)論一致。

圖13 CO濃度沿程變化曲線(xiàn)

隧道內(nèi)污染物瞬時(shí)點(diǎn)源一維隨流擴(kuò)散方程基本解為[13]

(5)

式中：C為隧道內(nèi)采樣點(diǎn)空氣污染物濃度，是時(shí)間t和位置x的函數(shù)，mg·m-3；C0為隧道內(nèi)測(cè)點(diǎn)初始時(shí)刻空氣污染物濃度，mg·m-3；u為隧道內(nèi)氣流速度，m·s-1；D為CO擴(kuò)散系數(shù)，m2·s。

將式(5)進(jìn)行變換可得

(6)

將圖13對(duì)應(yīng)的實(shí)測(cè)數(shù)據(jù)代入式(6)，可得出碴車(chē)不同車(chē)速時(shí)CO擴(kuò)散系數(shù)沿程變化曲線(xiàn)，如圖14所示。由圖14可知：車(chē)速越大，車(chē)體尾部CO擴(kuò)散系數(shù)越大；受尾部渦流的影響，距車(chē)尾近區(qū)，CO擴(kuò)散系數(shù)較大，擴(kuò)散較快，遠(yuǎn)離尾部渦流區(qū)，CO濃度擴(kuò)散系數(shù)趨于穩(wěn)定。

圖14 CO擴(kuò)散系數(shù)隨距離的變化

為便于計(jì)算，考慮CO擴(kuò)散系數(shù)與車(chē)速之間的關(guān)系，車(chē)速分別為4，5，6 m·s-1時(shí)對(duì)應(yīng)的CO擴(kuò)散系數(shù)均值分別為6.20，5.97，5.91 m2·s。

4 結(jié) 論

(1)出碴車(chē)分別以圖2中3條速度曲線(xiàn)運(yùn)行時(shí)，車(chē)頭前端氣流速度和壓力變化影響的區(qū)域長(zhǎng)度約為5 m，車(chē)尾影響的區(qū)域較長(zhǎng)；出碴車(chē)運(yùn)行速度越大，越靠近車(chē)輛尾部，環(huán)隙流速度越大；運(yùn)行歷時(shí)12 s后速度場(chǎng)均趨于穩(wěn)定；加速階段車(chē)體不同部位氣流速度均遠(yuǎn)大于勻速行駛階段，其中車(chē)尾渦流區(qū)氣流速度最大，為勻速運(yùn)行階段車(chē)尾渦流區(qū)氣流最大速度的2～3倍；隨著車(chē)速增大，車(chē)頭前端壓力及車(chē)尾負(fù)壓均增大，加速結(jié)束時(shí)，車(chē)頭氣流壓力最大值是勻速行駛時(shí)最大值的5倍，車(chē)尾氣流壓力最大值是勻速行駛時(shí)最大值的3.6倍。

(2)CO主要集中在車(chē)尾渦流區(qū)；出碴車(chē)勻速行駛時(shí)，距尾氣管出口越近，CO濃度越高，但超過(guò)標(biāo)準(zhǔn)限值30 mg·m-3的區(qū)域極小，且未達(dá)到人體呼吸高度；出碴車(chē)加速運(yùn)行區(qū)段為CO高濃度積聚區(qū)，其濃度超過(guò)標(biāo)準(zhǔn)限值30 mg·m-3。因此建議：出碴車(chē)盡量減少作業(yè)點(diǎn)，并以較高的速度勻速行駛，以減少CO排放；作業(yè)人員選擇避讓點(diǎn)時(shí)，盡量選擇在出碴車(chē)勻速行駛區(qū)段，且遠(yuǎn)離車(chē)輛尾氣管的一側(cè)。

(3)為了盡可能減小前后2輛出碴車(chē)行駛排放CO濃度的疊加，使CO累積濃度小于標(biāo)準(zhǔn)限值30 mg·m-3，建議前后2輛出碴車(chē)行駛間隔時(shí)間大于5 min。

(4)出碴車(chē)運(yùn)行過(guò)程中，CO擴(kuò)散規(guī)律符合瞬時(shí)點(diǎn)源一維擴(kuò)散方程基本解，即按照e指數(shù)衰減。車(chē)速越大，車(chē)體尾部CO擴(kuò)散系數(shù)越大；受尾部渦流的影響，距車(chē)尾近區(qū)，CO擴(kuò)散系數(shù)較大，遠(yuǎn)離尾部渦流區(qū)，CO擴(kuò)散系數(shù)趨于穩(wěn)定。