不同通風(fēng)參數(shù)對掘進工作面降溫的影響

辛 嵩，劉尚校，張 逍，朱曉鎮(zhèn)

（山東科技大學(xué) 安全與環(huán)境工程學(xué)院，山東 青島266590）

隨著煤炭開采強度增加以及開采工藝的改進，淺部資源日益減少，開采深部礦井資源已成為必然選擇。當(dāng)?shù)V井開采到一定深度，圍巖溫度達到幾十?dāng)z氏度，高溫、高濕的作業(yè)環(huán)境嚴重影響著開采進度[1]。我國2016 年版《煤礦安全規(guī)程》規(guī)定：當(dāng)掘進工作面空氣溫度超過26 ℃，必須縮短超溫地點作業(yè)人員的工作時間，當(dāng)掘進工作面的空氣溫度超過30℃必須停止作業(yè)[2]。獨頭巷道由于缺少回風(fēng)巷，往往是礦井熱害最為嚴重的地方，如何有效改善采掘工作面的熱環(huán)境狀況，成為眾多學(xué)者重點研究的方向。目前，國內(nèi)外學(xué)者已經(jīng)在數(shù)值模擬方面對掘進工作面的風(fēng)流場、溫度場進行了研究。姬建虎[3]得出掘進工作面射流沖擊換熱系數(shù)的關(guān)聯(lián)式，并對不同組合方式下?lián)Q熱情況進行了數(shù)值模擬，李勇[4]對巷道風(fēng)流流場及降溫前后溫度場的分布及變化進行了模擬分析。但很少有學(xué)者對掘進工作面不同通風(fēng)參數(shù)下進行較為全面的數(shù)值模擬，為此選取了趙樓煤礦7302 軌道巷掘進工作面60 m 為研究對象，使用ANSYS_FLUENT 數(shù)值模擬軟件對掘進工作面不同通風(fēng)條件下的熱環(huán)境進行數(shù)值模擬，并與現(xiàn)場實測數(shù)據(jù)進行對比分析。

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

湍流模型是表征地下熱環(huán)境流動特性的關(guān)鍵組成部分。使用商用CFD 程序（Fluent）求解質(zhì)量守恒方程（連續(xù)方程）、動量方程（Navier-Stokes 方程）、能量守恒方程以及紊流模型方程。計算了巷道內(nèi)的空氣流動和傳熱，將SIMPLE 算法用于對流項的二階格式[5-6]。對比分析了幾種常用的湍流模型選擇了工程中最常用的標(biāo)準(zhǔn)K-Epsilon 模型，該模型考慮了湍流動能k 及其與湍流黏性耦合的耗散率ε 2 個方程，在預(yù)測巷道內(nèi)的空氣速度和溫度場方面表現(xiàn)出最佳效果。

連續(xù)性方程：

式中：E 為總能量，J；hj為j 組分的生成焓；Jj為j 組分的擴散通量，kg/(m2·s)；Sh為化學(xué)反應(yīng)熱和自定義的任何其他體積熱源，J；k 為湍流動能，J。ε 為耗散率，%；μ 為流體的動態(tài)黏度，Pa·s；Gk為由于平均速度梯度產(chǎn)生的湍流動能；C1ε、C2ε分別為模型常數(shù)；σk、σε分別為k 方程和ε 方程對應(yīng)的湍流普朗特數(shù)；μt為湍流黏度，Pa·s。

2 物理模擬的構(gòu)建

1）巷道幾何模型。根據(jù)趙樓煤礦7302 軌道巷掘進工作面現(xiàn)場測量的尺寸，運用Rhino6.0 建立了長抽短壓式通風(fēng)系統(tǒng)下的掘進巷道簡化后的幾何模型，幾何模型由掘進機、巷道、壓入式風(fēng)筒、抽出式風(fēng)筒4 部分組成。巷道為梯形，凈高4 m，上凈寬4.8 m，下凈寬5.2 m，選用φ0.8 m 的風(fēng)筒，掘進機尺寸為11 m×3.5 m×1.8 m。

2）邊界條件命名與網(wǎng)格劃分。運用CFD－ICEM對趙樓煤礦7302 掘進工作面幾何模型進行邊界命名和網(wǎng)格劃分，將壓入式風(fēng)筒入風(fēng)口設(shè)定為入口邊界，巷道的末端斷面以及抽出式風(fēng)筒的出風(fēng)口設(shè)定為出口邊界。3Dinflation 提供了靠近墻邊界的高質(zhì)量網(wǎng)格生成功能，以解決物理特性的變化。為了得到更精確的模擬結(jié)果，對受風(fēng)流擾動和熱量交換影響的較為強烈的風(fēng)筒、掘進機附近區(qū)域進行了網(wǎng)格加密，共得到2 033 160 個網(wǎng)格，網(wǎng)格平均質(zhì)量0.88。

3）邊界條件具體參數(shù)設(shè)定。根據(jù)對趙樓煤礦7302 掘進工作面現(xiàn)場的實測數(shù)據(jù)分析，然后對邊界條件相關(guān)參數(shù)進行了設(shè)定。假設(shè)風(fēng)流為低速不可壓縮流體，密度符合近似假設(shè)，可忽略由流體的黏性力做功所引起的耗散熱。圍巖壁面溫度均勻分布，且熱物性參數(shù)為常數(shù)。假設(shè)巷道壁面粗糙度均勻，所有壁面施加無滑移邊界條件，假設(shè)壁面溫度均勻分布，巷道圍巖壁面及迎頭壁溫按照實測數(shù)據(jù)取平均值。計算湍流模型為Realizablek-epsilon 模型，開啟能量方程，采用SIMPLE 算法，待迭代平衡后，運用CFD-Post 對結(jié)果進行分析[7-8]。具體邊界條件設(shè)定為：①入口溫度：23 ℃；②入口風(fēng)量：300 m3/min；③抽出風(fēng)管出口風(fēng)速：－15 m/s；④巷道出口：自由回流；⑤圍巖壁面溫度：35 ℃；⑥掘進迎頭溫度：39 ℃；⑦掘進機壁面溫度：50 ℃。

3 模擬結(jié)果

研究為了量化壓入風(fēng)管位置對趙樓煤礦7302掘進工作面長抽短壓式通風(fēng)降溫效果的影響，結(jié)合礦井實際情況，固定變量壓風(fēng)筒距離迎頭距離為20 m，通過改變壓風(fēng)管出風(fēng)口與掘進巷道迎頭端面之間的距離Lf（Lf分別設(shè)為2.5、5、7.5、10 m），對巷道內(nèi)的風(fēng)流速度、溫度進行數(shù)值模擬，得到壓風(fēng)管不同布置方式下，巷道內(nèi)的風(fēng)流速度、溫度分布規(guī)律，通過比較可以確定壓風(fēng)管的合理位置。

為明確抽壓風(fēng)量配比對趙樓煤礦7302 軌道巷掘進工作面通風(fēng)降溫效果的影響，結(jié)合礦井實際情況，將壓入風(fēng)量Q 設(shè)為300 m3/min，通過改變抽壓風(fēng)量配比B（B 分別設(shè)為1、1.5、2）分析7302 掘進工作面通風(fēng)系統(tǒng)在不同抽壓風(fēng)量配比情況下的通風(fēng)降溫效果，從而得到通風(fēng)效果最佳的抽壓風(fēng)量配比。

3.1 壓風(fēng)管出風(fēng)口位置對風(fēng)流速度的影響

不同Lf條件下掘進巷道內(nèi)風(fēng)流流線分布如圖1。由圖1 可以看出，長抽短壓式通風(fēng)條件下巷道內(nèi)風(fēng)流的運移狀況要比單一壓入式通風(fēng)條件下更加復(fù)雜。設(shè)定距離迎頭前30 m 區(qū)域為工作區(qū)域，當(dāng)Lf=2.5 m 時，工作區(qū)域的局部區(qū)域的風(fēng)流速度都保持在0.3 m/s 以下，風(fēng)流經(jīng)自由擴張后快速進入到作業(yè)區(qū)域，由于受到掘進機的阻擋作用，反向回風(fēng)不暢導(dǎo)致作業(yè)區(qū)域的回流速度過低，不符合煤礦安全規(guī)程對空氣的供給量規(guī)定的最小風(fēng)速。當(dāng)Lf=10 m時，由于掘進機的阻擋作用，限制了射流邊界層的發(fā)展擴散，這部分風(fēng)流由于掘進機的阻擋作用形成了回流，從而保證了大部分區(qū)域的風(fēng)速在0.3 m/s以上。由于壓風(fēng)口、吸風(fēng)口處的卷吸作用，在距離掘進迎頭20 m 以內(nèi)的區(qū)域內(nèi)，形成了多個不規(guī)則渦流風(fēng)流場，當(dāng)Lf=7.5、10 m 時，明顯風(fēng)流的速度要高于其他條件下的風(fēng)流速度，考慮到獨頭巷道的除塵效果，工作區(qū)域內(nèi)的風(fēng)流速度不宜過大[9-10]。

④應(yīng)將統(tǒng)籌城鄉(xiāng)水利基礎(chǔ)設(shè)施和公共服務(wù)作為城鄉(xiāng)經(jīng)濟社會發(fā)展一體化的重要環(huán)節(jié)。統(tǒng)籌考慮城鄉(xiāng)水利基礎(chǔ)設(shè)施建設(shè)和水利公共服務(wù)，同步納入政府公共財政預(yù)算支出項目；建立城鄉(xiāng)經(jīng)濟社會一體化發(fā)展制度，并將加快農(nóng)村飲水安全建設(shè)和實施農(nóng)村水環(huán)境治理作為加快發(fā)展農(nóng)村公共事業(yè)的重要內(nèi)容。

圖1 不同Lf 條件下掘進巷道內(nèi)風(fēng)流流線分布Fig.1 Distribution of wind flow lines in tunneling under different Lf conditions

3.2 壓風(fēng)管出風(fēng)口位置對風(fēng)流溫度的影響

不同的Lf值下，巷道后端的溫度分布沒有明顯的差異?？梢哉J為，巷道后端的風(fēng)流主要是由回流組成的，風(fēng)流與機械設(shè)備和圍巖之間發(fā)生了相同的熱傳遞，且大部分回流通過抽出風(fēng)管排出巷道，只有少部分污風(fēng)回流到巷道后端。相反，巷道工作區(qū)域的溫度分布存在明顯差異。不同Lf條件下巷道內(nèi)風(fēng)流速度在y=0 平面上的變化情況如圖2。從圖2可以看出，較小的Lf在作業(yè)區(qū)域具有更好的冷卻性能，然而，Lf過小時則會起到相反的效果，因為Lf過小時，由于掘進機的阻礙作用，新鮮冷風(fēng)流反向回風(fēng)不暢導(dǎo)致作業(yè)區(qū)域的部分區(qū)域溫度較高，而Lf設(shè)置為5 m 時，工作區(qū)域的平均溫度最低，冷風(fēng)的擴散面積最大。

綜上所述，當(dāng)Lf設(shè)置為5 m 時，工作區(qū)域的大部分風(fēng)速符合指標(biāo)，溫度約為27 °C，對礦工來說較為舒適。因此，在后續(xù)分析中應(yīng)采用Lf=5 m 進行研究。在優(yōu)化通風(fēng)降溫系統(tǒng)時，應(yīng)注意風(fēng)管位置對通風(fēng)降溫性能的影響。

圖2 不同Lf 條件下巷道內(nèi)風(fēng)流速度y=0 平面上的變化情況Fig.2 Changes of the air velocity in the tunnel on the y=0 plane under different Lf conditions

3.3 抽壓風(fēng)量配比對風(fēng)流速度的影響

不同Lf值條件下巷道內(nèi)風(fēng)流溫度在y=0 平面上的變化情況如圖3。不同B 值條件下巷道內(nèi)風(fēng)流速度在y=0 平面上的變化情況如圖4。

圖3 不同Lf 值條件下巷道內(nèi)風(fēng)流溫度在y=0 平面上的變化情況Fig.3 The change of the airflow temperature in the tunnel on the y=0 plane under different Lf conditions

圖4 不同B 值條件下巷道內(nèi)風(fēng)流速度在y=0 平面上的變化情況Fig.4 B variation of the air velocity in the tunnel on the y=0 plane under different conditions

3 種不同抽壓風(fēng)量配比下，作業(yè)區(qū)域的風(fēng)流速度變化不大，相反，在巷道后端風(fēng)流速度產(chǎn)生了明顯變化，當(dāng)抽出風(fēng)量等于壓入風(fēng)量（即B=1）時，巷道后端壓抽重疊段風(fēng)流平均速度低于0.2 m/s，當(dāng)抽出風(fēng)量大于壓入風(fēng)量（即B=1.5 及B=2）時，整個巷道中的風(fēng)速基本滿足相關(guān)規(guī)定，絕大部分區(qū)域大于0.3 m/s，然而當(dāng)B=2 時抽風(fēng)口附近得風(fēng)流較紊亂，整個巷道的風(fēng)速較大，使得部分污風(fēng)擴散至整個巷道，對通風(fēng)效果起到反作用，且會加大通風(fēng)能耗。

3.4 抽壓風(fēng)流配比對風(fēng)流溫度的影響

當(dāng)抽出風(fēng)量大于壓入風(fēng)量時，可對巷道內(nèi)的高溫環(huán)境產(chǎn)生有利影響，因此混合通風(fēng)可以通過適當(dāng)提高抽出風(fēng)量來改善巷道內(nèi)的高溫環(huán)境。然而，過高的抽出風(fēng)量并不能有效改善巷道內(nèi)的熱環(huán)境，不同B 值條件下巷道內(nèi)風(fēng)流溫度在y=0 平面上的變化情況如圖5，不同B 值條件下巷道內(nèi)各斷面上風(fēng)流溫度的變化情況如圖6。

圖5 不同B 值條件下巷道內(nèi)風(fēng)流溫度在y=0 平面上的變化情況Fig.5 B change of airflow temperature in the tunnel on the y=0 plane under different conditions

圖6 不同B 值條件下巷道內(nèi)各斷面上風(fēng)流溫度的變化情況Fig.6 B changes of airflow temperature on each section in the tunnel under different conditions

當(dāng)B=2 時，壓入風(fēng)管與抽出風(fēng)管正負壓作用下的渦流區(qū)內(nèi)的風(fēng)流紊亂，使得大部分高溫風(fēng)流擴散到巷道后端，在距離掘進迎頭10 m 處的斷面平均溫度27.8 ℃，巷道內(nèi)平均溫度較高，對降溫過程產(chǎn)生負面影響。

當(dāng)B=1.5 時，抽壓風(fēng)流配合合理，距離掘進迎頭1、5、10 m 處斷面的平均溫度分別為26.4、27.0、27.3℃時達到了較好的通風(fēng)降溫效果。

4 數(shù)值模擬與實測驗證

為了全面掌握巷道中的風(fēng)流參數(shù)以及驗證數(shù)值模擬的有效性，根據(jù)斷面與掘進迎頭的選取了8 個斷面，分別距離迎頭1、5、10、15、20、25、30、60 m，每個巷道斷面設(shè)置了4 個測量點，巷道斷面測點布置圖如圖7。

待風(fēng)機開啟5 min 左右，風(fēng)流場穩(wěn)定后，使用風(fēng)速計測量巷道內(nèi)的風(fēng)流速度，采用用干式溫度計、紅外線測溫儀測量巷道內(nèi)空氣溫度。

各測點風(fēng)流速度大小、方向及風(fēng)流溫度大小實測數(shù)據(jù)與數(shù)值模擬對比見表1，表2。

圖7 巷道斷面測點布置圖Fig.7 Arrangement of measuring points of tunnel cross section

表1 現(xiàn)場實測風(fēng)流速度、方向與數(shù)值模擬對比Table 2 Field measured wind speed, direction and numerical simulation comparison

表2 現(xiàn)場實測風(fēng)流速度與數(shù)值模擬對比Table 2 Comparison of field measured wind speed and numerical simulation

5 結(jié) 論

1）通過改變壓風(fēng)管出風(fēng)口與掘進巷道迎頭端面之間的距離Lf（分別設(shè)為2.5、5、7.5、10 m），對巷道內(nèi)的風(fēng)流速度、溫度進行數(shù)值模擬，結(jié)果顯示當(dāng)Lf設(shè)置不同值時，工作區(qū)域形成的渦流場大小不一，對通風(fēng)降溫系統(tǒng)影響較大。

2）將壓入風(fēng)量設(shè)為Q=300 m3/min，通過改變抽壓風(fēng)量配比B（B 分別設(shè)為1、1.5、2）分析掘進工作面通風(fēng)系統(tǒng)在不同抽壓風(fēng)量配比工況下的通風(fēng)降溫效果，結(jié)果顯示當(dāng)壓入風(fēng)管距離掘進迎頭5 m，抽壓風(fēng)量配比為1.5 時，工作區(qū)域平均溫度27 ℃，長抽短壓式通風(fēng)系統(tǒng)可達到最佳的通風(fēng)降溫效果。

3）將數(shù)值模擬結(jié)果應(yīng)用到趙樓煤礦7302 軌道巷掘進工作面通風(fēng)優(yōu)化中，經(jīng)過現(xiàn)場實測，實測風(fēng)速與數(shù)值模擬結(jié)果平均相對誤差10.07%，實測風(fēng)溫與數(shù)值模擬結(jié)果平均相對誤差1.60%，數(shù)值模擬結(jié)果較為準(zhǔn)確。