智曉峰，鞠興軍，張周愛，繆衛(wèi)峰，杜勇志，張 洋，楊小彬，程虹銘

（1.大雁公司（神寶能源）寶日希勒露天煤礦，內(nèi)蒙古 呼倫貝爾 021008；2.中國礦業(yè)大學(xué)（北京）應(yīng)急管理與安全工程學(xué)院，北京 100083）

霧是空氣中水蒸氣接近飽和時(shí)，氣象參數(shù)（溫度、濕度、氣壓等）發(fā)生變化，多余水蒸氣凝結(jié)成細(xì)小水滴析出、匯聚而成［1］。工程中，霧害主要是降低作業(yè)空間能見度，給作業(yè)空間內(nèi)的人員帶來安全隱患。張清琳［2］分析了城市管廊工程中空氣霧化原因，提出一種新型的綜合管廊通風(fēng)除霧系統(tǒng)；薛永慶等［3］采用實(shí)測和數(shù)值模擬的方法分析冬季引水隧道霧氣產(chǎn)生原因，提出升溫除濕的方式減小引水隧道霧氣的產(chǎn)生；孫寶芝等［4］基于熱力學(xué)原理分析了地下海水坑道霧氣形成機(jī)制，提出升溫、降溫、通風(fēng)等幾種消霧方法并進(jìn)行了現(xiàn)場應(yīng)用；李增華等［5］、吳吉男等［6］、秦憲禮等［7］、祁學(xué)謙等［8］對(duì)井工煤礦巷道產(chǎn)霧原因進(jìn)行了分析并提出了相應(yīng)的治理措施。神寶能源寶日希勒露天煤礦冬季鏟裝作業(yè)中不同程度地發(fā)生了霧氣，加之，采場地勢低，空氣對(duì)流慢，霧氣縈繞在電鏟鏟斗與煤壁之間，消散慢，能見度低，電鏟作業(yè)司機(jī)視線受阻，既影響鏟裝效率又留下安全隱患。

為此，現(xiàn)場測定采場煤壁及鏟裝空間的溫度、濕度和風(fēng)速，基于熱力學(xué)原理分析露天采場霧氣產(chǎn)生原因，采用數(shù)值模擬的方法分析采場霧氣分布情況，提出增設(shè)風(fēng)機(jī)干預(yù)采場流場的方法進(jìn)行除霧，并模擬分析了不同風(fēng)機(jī)組合下采場流場分布及霧氣消散情況。

1 現(xiàn)場測試及結(jié)果分析

結(jié)合露天煤礦采場冬季霧氣形成機(jī)理，對(duì)寶日希勒露天煤礦采場在冬季12 月份鏟裝過程中形成霧氣區(qū)域的空氣溫度、濕度、風(fēng)速進(jìn)行現(xiàn)場測試。測點(diǎn)布置在采場煤壁及煤壁前方15 m 范圍內(nèi)，測點(diǎn)布置間隔為5 m，其中煤壁處測點(diǎn)（A1、A2、A3、A4、B1、C1、D1）主要測定煤壁剛暴露時(shí)涌出空氣的溫度、濕度，但考慮到煤壁的不穩(wěn)定性，可適當(dāng)遠(yuǎn)離煤壁測定，其余9 個(gè)測點(diǎn)測定采場空氣的溫度、濕度和風(fēng)速。測定編號(hào)和布置如圖1 所示。各測點(diǎn)參數(shù)測5 次，取5 次測定結(jié)果的平均值作為測定的最終結(jié)果。測定的最終結(jié)果見表1～表3。

圖1 采場空間測點(diǎn)布置

表1 采場空間溫度測定結(jié)果 單位：℃

表3 采場空間風(fēng)速測定結(jié)果 單位：m/s

霧是空氣相對(duì)濕度達(dá)到100%時(shí)空氣中的水蒸氣凝結(jié)形成細(xì)小水滴漂浮在空氣中形成的。空氣相對(duì)濕度的大小由濕空氣中水蒸氣的分壓與相同溫度下飽和水蒸氣的分壓比值決定，空氣相對(duì)濕度的增大可以是濕空氣中水蒸氣的分壓增大或空氣溫度降低所造成。從表1 中可以看出，新暴露的煤壁溫度與采場空間溫度相差較大，采場空間飽和水蒸氣分壓小于新暴露煤壁處的飽和水蒸氣分壓，煤壁瞬時(shí)涌出的濕空氣在采場空間達(dá)到飽和狀態(tài)析出細(xì)小水滴成霧；并且采場空間空氣相對(duì)濕度低于煤壁附近空氣（表2），相對(duì)濕度較大的濕空氣涌入相對(duì)濕度較小的空間，為霧氣源源不斷生產(chǎn)提供了條件；采場空間地勢低，兩面為煤壁，采場流場微弱（表3），亦不利于霧氣消散；同時(shí)，鏟裝過程中形成的大量煤塵作為凝結(jié)核，也為采場空間的汽-水轉(zhuǎn)化提供了依托。

表2 采場空間相對(duì)濕度測定結(jié)果 %

2 采場霧氣產(chǎn)生及分布數(shù)值模擬

2.1 采場霧氣分布數(shù)學(xué)模型

采用Ansys Fluent 軟件對(duì)冬季露天采場霧氣產(chǎn)生及分布規(guī)律進(jìn)行數(shù)值計(jì)算與分析。霧氣為細(xì)小水滴的混合體，露天采場霧氣形成中又包裹著煤塵，可以用離散相模型（DPM）模擬分析細(xì)小水滴在采場空間中的分布規(guī)律，假定露天采場霧氣顆粒生成過程中滿足牛頓第二定律，進(jìn)一步可運(yùn)用拉格朗日法對(duì)霧氣顆粒運(yùn)動(dòng)方程進(jìn)行積分求解，可求得外界流場作用下霧氣顆粒的分布規(guī)律，霧氣顆粒的受力方程［9］可表示為

式中：u為流體速度；up為顆粒速度；ρ為流體密度；ρp為顆粒密度；Fx為附加加速度；FD（u-up）為顆粒的單位質(zhì)量曳力；gx為重力加速度。

相對(duì)雷諾數(shù)（顆粒雷諾數(shù)）

式中：dp為顆粒直徑；CD為曳力系數(shù)；α1、α2、α3根據(jù)相對(duì)雷諾數(shù)的范圍取不同的值。

2.2 采場幾何模型

依據(jù)寶日希勒露天煤礦采場布置情況及相關(guān)參數(shù)，建立三維幾何計(jì)算模型。采場空間長40 m、寬20 m、高30 m，鏟車位置在長度方向20 m 處，煤壁1、2、3 為電鏟司機(jī)可視區(qū)域，鏟車與煤壁1、2、3 間為霧氣生成區(qū)域和鏟車作業(yè)區(qū)，幾何模擬具體尺寸及布置情況如圖2 所示。

圖2 采場物理模型

2.3 采場霧氣分布規(guī)律

圖3 為Fluent 模擬得到的采場霧氣分布規(guī)律?？梢钥闯?，在近煤壁處生成的霧氣濃度相對(duì)較低，離煤壁距離增大，采場空氣溫度、濕度與煤壁涌出的濕空氣差異增大，霧氣生成濃度增加，最大濃度達(dá)17.7 g/m3。由于鏟車司機(jī)的視線高度大約在10 m，截取了采場10 m 高度處的霧氣濃度分布情況，如圖3（b）所示，可以看出，鏟車司機(jī)的視線平面內(nèi)，霧氣的分布密集范圍在10 m 左右，也就是整個(gè)鏟斗的工作區(qū)域均分布著濃霧。

圖3 采場霧氣分布規(guī)律

3 風(fēng)機(jī)干預(yù)下采場流場分布及除霧效果分析

霧氣的異常生成常出現(xiàn)在空氣流通不暢、空氣溫度、濕度急劇變化的管廊、地下空間工程等，一般可采取增溫、加大通風(fēng)等方法降低空氣的相對(duì)濕度，達(dá)到除霧的目的。露天采場為半封閉空間，空間范圍大，增溫除霧達(dá)不到良好的效果，為此，提出增設(shè)風(fēng)機(jī)以改變采場流場的除霧方法。

3.1 采場流場模型

露天采場空間大，要形成有效的流場風(fēng)機(jī)需要較大的出口流速，并且受煤壁的影響，風(fēng)流在近煤壁區(qū)域會(huì)形成渦流場，F(xiàn)luent 軟件中的RNGk-ε模型是使用重整化群理論推導(dǎo)出來的，相比標(biāo)準(zhǔn)k-ε模型，RNGk-ε模型在模擬高速流和漩渦流方面提高了精度［10-11］，更適合模擬本文風(fēng)機(jī)作用下采場流場分布情況。采場流場數(shù)學(xué)模型［12］及幾何模型如下：

1）流場連續(xù)方程：

式中：ui為速度張量；xi為位移張量。

2）流場動(dòng)量方程：

式中：i=x、y、z，作為下標(biāo)時(shí)代表X、Y、Z軸，否則代表位移張量；μ為動(dòng)力黏度。

3）RNGk-ε湍流模型的動(dòng)能方程：

4）RNGk-ε湍流模型的擴(kuò)散方程：

式中：k為湍動(dòng)能，m2/s2；ε為耗 散率，m2/s3；Gk為 由層流速度梯度產(chǎn)生的湍流動(dòng)能；Gb為由浮力產(chǎn)生的湍流動(dòng)能；αε、αk分別為ε方程和k方程的湍流普朗特?cái)?shù)；Sk、Sε為用戶自定義；μi為湍流動(dòng)力黏滯系數(shù)，Pa·s；YM為可壓縮湍流中過渡的擴(kuò)散產(chǎn)生的波動(dòng)。其中，μt為湍流黏性系數(shù)，模型常數(shù)分別為C1ε=1.42，C2ε=1.68。

5）幾何模型：根據(jù)電鏟情況，在電鏟不同位置增設(shè)3 臺(tái)風(fēng)機(jī)以改變采場空間流場。如圖2 所示，在離地面3 m 處，電鏟兩側(cè)分別布置兩臺(tái)風(fēng)機(jī)（風(fēng)機(jī)1、2），再在風(fēng)機(jī)2 上方2 m 處內(nèi)錯(cuò)1 m 布置風(fēng)機(jī)3。

3.2 風(fēng)機(jī)干預(yù)下采場流場分布情況

為探究風(fēng)機(jī)作用下采場流場的分布情況，分別模擬了單臺(tái)風(fēng)機(jī)、兩臺(tái)風(fēng)機(jī)和3 臺(tái)風(fēng)機(jī)的作用效果，風(fēng)機(jī)出口風(fēng)速設(shè)置為30 m/s。截取三維流場中10 m高度處平面，如圖4 所示，圖4（a）為單臺(tái)風(fēng)機(jī)1 作用效果，可以看出流場流速由風(fēng)機(jī)出口逐漸向煤壁衰減，煤壁處最小風(fēng)速為8.5 m/s，圖4（b）、（c）、（d）分別為水平兩臺(tái)風(fēng)機(jī)1 和2、垂直兩臺(tái)風(fēng)機(jī)2 和3 及3 臺(tái)風(fēng)機(jī)的作用效果，煤壁處的風(fēng)速分別為10.4、11.0、11.6 m/s。從流場的空間分布上看（以流場流速11 m/s 等值線為例），單臺(tái)風(fēng)機(jī)的作用范圍最小，3 臺(tái)風(fēng)機(jī)的作用范圍最大，垂直兩臺(tái)風(fēng)機(jī)減小風(fēng)速的衰減，但其作用范圍小于水平布置的兩臺(tái)風(fēng)機(jī)。

圖4 風(fēng)機(jī)干預(yù)下采場流場分布

3.3 風(fēng)機(jī)干預(yù)下除霧效果

DPM 模型中顆粒的運(yùn)移分布受曳力的影響，見式（1），外在流場作用下可以增大顆粒的單位質(zhì)量曳力，促進(jìn)顆粒的運(yùn)移。通過增設(shè)風(fēng)機(jī)改變采場流場分布，進(jìn)而促進(jìn)霧氣顆粒運(yùn)移，達(dá)到除霧效果。不同風(fēng)機(jī)組合下流場分布及影響范圍不同，本節(jié)模擬分析不同風(fēng)機(jī)組合下除霧效果，如圖5 所示。

圖5（a）為單臺(tái)風(fēng)機(jī)1 作用下霧氣顆粒濃度分布情況。可以看出，單臺(tái)風(fēng)機(jī)作用下，采場霧氣濃度開始下降，但是采場拐角處霧氣濃度依然很高，最大霧氣濃度約為15.2 g/m3。圖5（b）為兩臺(tái)水平風(fēng)機(jī)1、2 作用下霧氣顆粒濃度分布情況，兩臺(tái)水平風(fēng)機(jī)流場作用范圍大，采場霧氣濃度顯著下降，采場拐角處積聚的最大霧氣濃度為4.53 g/m3。圖5（c）為兩臺(tái)垂直風(fēng)機(jī)2、3 作用下霧氣顆粒濃度分布情況，兩臺(tái)垂直風(fēng)機(jī)減小了流場流速的衰減，但是影響范圍小于兩臺(tái)水平風(fēng)機(jī)，采場霧氣消散情況劣于兩臺(tái)水平風(fēng)機(jī)作用效果，采場大部分區(qū)域仍然遍布霧氣，霧氣最大濃度為9.94 g/m3。圖5（d）為模擬試驗(yàn)3 臺(tái)風(fēng)機(jī)作用下霧氣顆粒濃度分布情況，3 臺(tái)風(fēng)機(jī)作用下流場影響范圍大，流速衰減小，除霧效果顯著，只在采場拐角處有殘霧存在，最大濃度在3.5 g/m3左右。

圖5 風(fēng)機(jī)干預(yù)下采場霧氣濃度分布

以上數(shù)值模擬試驗(yàn)表明：露天采場增設(shè)風(fēng)機(jī)改變流場的方法可有效消除采場霧氣，且不同的風(fēng)機(jī)組合下除霧效果不同。實(shí)際運(yùn)用中，可根據(jù)現(xiàn)場霧氣分布情況選擇不同的風(fēng)機(jī)組合形式，例如：霧氣濃度較小時(shí)，可開設(shè)一臺(tái)風(fēng)機(jī)，既可達(dá)到除霧效果又減小電能消耗；霧氣濃度較大時(shí)，初期開設(shè)3 臺(tái)風(fēng)機(jī)，可快速消除鏟車司機(jī)前方霧氣，隨著鏟車司機(jī)視野的開闊可減少風(fēng)機(jī)數(shù)量，通過水平兩臺(tái)或垂直兩臺(tái)風(fēng)機(jī)進(jìn)一步消除采場霧氣。

4 結(jié)論

通過對(duì)神寶能源寶日希勒露天煤礦冬季采場霧氣進(jìn)行現(xiàn)場實(shí)測，結(jié)合理論分析和數(shù)值模擬，對(duì)該露天礦冬季采場霧氣形成機(jī)理及風(fēng)機(jī)干預(yù)除霧效果進(jìn)行分析，得到以下結(jié)論：

1）該露天礦冬季采場空間內(nèi)以電鏟鏟裝過程中形成的大量煤塵作為凝結(jié)核依托，煤中濕空氣與采場中低溫、低濕空氣混合，降溫增濕及混合增濕兩種方式下，采場空間局部空氣相對(duì)濕度達(dá)到100%后形成霧氣，霧氣形成過程中關(guān)鍵參數(shù)是空氣溫度和濕度。

2）煤中濕空氣與采場中低溫、低濕空氣混合后，電鏟與煤壁中間是霧氣形成的主要區(qū)域，數(shù)值模擬顯示霧氣生成的范圍在10 m 左右，生成的霧氣最大濃度為17.7 g/m3，已嚴(yán)重阻擋鏟車司機(jī)的作業(yè)視線。

3）增設(shè)風(fēng)機(jī)的方法可有效改善采場空間的流場，數(shù)值模擬得到：單臺(tái)風(fēng)機(jī)作用下，煤壁拐角處風(fēng)速為8.5 m/s；兩臺(tái)風(fēng)機(jī)作用下，煤壁拐角處風(fēng)速有所增加，相比兩臺(tái)風(fēng)機(jī)水平布置，垂直布置的兩臺(tái)風(fēng)機(jī)作用范圍相對(duì)較?。? 臺(tái)風(fēng)機(jī)同時(shí)布置，作用范圍最大。

4）流場的作用改變了霧氣顆粒的運(yùn)移，可有效地消散霧氣。數(shù)值模擬得出：當(dāng)增設(shè)1 臺(tái)風(fēng)機(jī)時(shí)，煤壁拐角處的霧氣的最大濃度約為15.2 g/m3；增設(shè)兩臺(tái)水平風(fēng)機(jī)時(shí)，煤壁拐角處的霧氣的最大濃度約為4.53 g/m3；增設(shè)垂直兩臺(tái)水平風(fēng)機(jī)時(shí)，煤壁拐角處的霧氣的最大濃度約為9.94 g/m3；增設(shè)3 臺(tái)風(fēng)機(jī)時(shí)，煤壁拐角處的霧氣幾乎全部消散。