楊文宇，張成法，逄崇晨，馬 輝

（1.山東科技大學(xué) 礦業(yè)與安全工程學(xué)院，山東 青島 266590；2.招金礦業(yè)股份有限公司 夏甸金礦，山東 招遠(yuǎn) 265418）

高地溫問題已成為礦業(yè)工程常見的地質(zhì)災(zāi)害，在深部采礦過程、采礦工程活動(dòng)必然要受到礦井熱環(huán)境的制約。隨著我國(guó)淺部礦產(chǎn)資源逐漸減少，地下開采的深度也越來(lái)越大，隨之而來(lái)的熱害問題也變得愈加嚴(yán)峻[1]。為滿足工業(yè)衛(wèi)生標(biāo)準(zhǔn)對(duì)深部礦井要求的熱環(huán)境條件，礦山企業(yè)一般采用2個(gè)手段來(lái)控制獨(dú)頭巷道的熱環(huán)境，即增大風(fēng)量與降低風(fēng)溫[2]。

隨著上世紀(jì)60年代計(jì)算流體力學(xué)（CFD）的出現(xiàn)，熱害礦井掘進(jìn)工作面的數(shù)值模擬研究取得了快速發(fā)展。例如，日本學(xué)者中山伸介等運(yùn)用CFD數(shù)值模擬研究了掘進(jìn)工作面風(fēng)流流場(chǎng)，高建良同樣運(yùn)用數(shù)值模擬的方法對(duì)掘進(jìn)工作面換熱系數(shù)的分布與風(fēng)溫預(yù)測(cè)開展了研究[3]。之后，F(xiàn)LUENT與Airpak等通用商業(yè)數(shù)值模擬軟件與CLIMSIM，MIVENA，Ventsim，VUMA，Multiflux與Kduct等 面向礦山熱害的專用數(shù)值模擬軟件在礦井掘進(jìn)工作面的研究工作中得到了廣泛應(yīng)用，取得了較多的研究成果[3]。

對(duì)于金屬礦山，隨著金屬礦山機(jī)械化程度的提高，以鏟運(yùn)機(jī)為代表的地下無(wú)軌設(shè)備得到了前所未有的發(fā)展，內(nèi)燃機(jī)產(chǎn)生的熱量隨之加大目前對(duì)獨(dú)頭巷道的熱環(huán)境數(shù)值模擬研究較少考慮通風(fēng)障礙物的影響，為研究通風(fēng)降溫的最優(yōu)能耗運(yùn)作方式，本文以金屬礦山獨(dú)頭巷道為研究對(duì)象，在考慮通風(fēng)。

1 幾何模型的建立

障礙物如鏟運(yùn)機(jī)存在情況下對(duì)其進(jìn)行熱環(huán)境數(shù)值模擬研究。

本研究建立的幾何模型以實(shí)測(cè)巷道為參照，建立等比例的幾何模型，將其計(jì)算區(qū)域定為如圖1所示的長(zhǎng)方體區(qū)域，其幾何尺寸為：X×Y×Z=40×4.9×4.6m3，為真實(shí)再現(xiàn)作業(yè)環(huán)境，在模型中還包括作業(yè)人員2人，礦工身高定為1.75m，距離迎頭3m，鏟運(yùn)機(jī)1臺(tái)尺寸為8.0×2.5×2.0m3。獨(dú)頭巷道中配備膠皮風(fēng)筒，其邊長(zhǎng)為800mm，風(fēng)筒長(zhǎng)30m，因此風(fēng)筒口距離迎頭10m，風(fēng)筒位于巷道頂板中部。

2 假設(shè)及邊界條件

（1）假設(shè)風(fēng)流為低速不可壓縮氣體，密度符合近似假設(shè)，圍巖表面熱量均勻分布，且熱物性參數(shù)為常數(shù)，壁面粗糙度均勻；

（2）出入口邊界：其中入口邊界的風(fēng)流初始速度為v=12m/s，送風(fēng)溫度t=25℃，相對(duì)濕度為70%；出口邊界為自由回流；

（3）壁面邊界：給所有壁面施加無(wú)滑動(dòng)邊界條件，假設(shè)表面溫度均勻分布，巷道兩幫及頂板熱量按照實(shí)測(cè)數(shù)據(jù)取平均值：134.8W/m2；

（4）人體熱舒適邊界：井下工作人員為繁重體力勞動(dòng)作業(yè)，新陳代謝率為400W/m2，鏟運(yùn)機(jī)平均散熱功率為66.8kW。

3 網(wǎng)格劃分

由于網(wǎng)格的質(zhì)量直接影響計(jì)算結(jié)果和精度，網(wǎng)格劃分是數(shù)值模擬中非常重要的一環(huán)，本模擬首先采用非結(jié)構(gòu)化網(wǎng)劃分器的默認(rèn)設(shè)置進(jìn)行粗糙網(wǎng)格的劃分。判斷劃分完的粗糙網(wǎng)格是否可以充分表示巷道的物理模型，并確定其是否能夠滿足網(wǎng)格劃分規(guī)則的最小網(wǎng)格量。運(yùn)用得到的粗糙網(wǎng)格進(jìn)行一次求解，并查看求解出來(lái)的結(jié)果是否合理，預(yù)測(cè)所需要的計(jì)算時(shí)間。將其作為后續(xù)進(jìn)行更為細(xì)致的網(wǎng)格劃分，求解更精確結(jié)果的評(píng)估依據(jù)。對(duì)劃分網(wǎng)格進(jìn)行長(zhǎng)寬比、扭曲度的網(wǎng)格質(zhì)量檢查。最終生成粗糙網(wǎng)格數(shù)33304，節(jié)點(diǎn)數(shù)38283，細(xì)化網(wǎng)格數(shù)403725，節(jié)點(diǎn)數(shù)426579。各方向網(wǎng)格單元對(duì)應(yīng)最大尺寸為總尺寸的1/20，對(duì)風(fēng)筒送風(fēng)口、巷道出風(fēng)口、與鏟運(yùn)機(jī)周邊溫度梯度和速度梯度比較大的地方網(wǎng)格局部加密細(xì)化。

圖1 獨(dú)頭巷道幾何模型

4 改變送風(fēng)速度對(duì)巷道熱環(huán)境的影響

為了選擇更準(zhǔn)確的湍流模型，在對(duì)獨(dú)頭巷道進(jìn)行數(shù)值模擬研究之前，用對(duì)比分析的方法，在各類參數(shù)均相同的情況下，將幾個(gè)常用的湍流模型，如零方程模型、標(biāo)準(zhǔn)k-ε模型和RNG k-ε模型所計(jì)算出的巷道風(fēng)速與風(fēng)溫與利用風(fēng)速儀與溫度儀進(jìn)行實(shí)測(cè)的結(jié)果進(jìn)行比較與分析，本文選用RNG k-ε模型作為計(jì)算模型，進(jìn)行下一步的數(shù)值模擬實(shí)驗(yàn)。

將風(fēng)管出風(fēng)口的送風(fēng)速度分別設(shè)定為v1=12m/s，v2=13m/s，v3=14m/s，v4=15m/s，其溫度場(chǎng)數(shù)值模擬結(jié)果分別如圖2所示。

圖2 XY面Z=2.3m處速度分布（（a）v1=15m/s（b）v2=14m/s（c）v3=13m/s（d）v4=12m/s）

從圖2可以看出從風(fēng)筒口射出的有限空間射流，在獨(dú)頭工作面會(huì)形成射流、回流和渦流三種氣流形態(tài)。國(guó)際隧道協(xié)會(huì)對(duì)空氣的供給量最小風(fēng)速應(yīng)大于0.3m/s，當(dāng)風(fēng)管出風(fēng)口風(fēng)速為12m/s時(shí)在X軸方向上存在大范圍的風(fēng)速不達(dá)標(biāo)區(qū)域（風(fēng)速小于0.3m/s），這是由于受到鏟運(yùn)機(jī)等通風(fēng)障礙物阻塞作用，反向回風(fēng)不暢導(dǎo)致回流速度過低。當(dāng)風(fēng)管出風(fēng)口風(fēng)速為15m/s時(shí)，鏟運(yùn)機(jī)的后部邊壁限制了射流邊界層的發(fā)展擴(kuò)散，這部分風(fēng)流由于鏟運(yùn)機(jī)的阻塞作用形成了回流，從而保證了大部分區(qū)域的風(fēng)速在0.3m/s。使得15m/s時(shí)風(fēng)速不達(dá)標(biāo)區(qū)域最小。由圖3還可以看到，當(dāng)風(fēng)速?gòu)?5m/s改變到12m/s時(shí)，風(fēng)流在鏟運(yùn)機(jī)與迎頭區(qū)域形成了較強(qiáng)的回流區(qū)域，使得此區(qū)域的空氣齡較無(wú)鏟運(yùn)機(jī)時(shí)明顯增大。由此可見，在實(shí)施通風(fēng)降溫工程時(shí)必須考慮鏟運(yùn)機(jī)或掘進(jìn)機(jī)對(duì)風(fēng)流流場(chǎng)的阻塞作用，在鏟運(yùn)機(jī)與巷道壁圍成的區(qū)域內(nèi)還會(huì)形成風(fēng)流風(fēng)道，起到促進(jìn)通風(fēng)的作用。

當(dāng)風(fēng)流從15m/s至12m/s的速度從風(fēng)筒射出，巷道內(nèi)空氣的平均溫度與送風(fēng)溫度相差3.2℃左右。隨著送風(fēng)速度的逐漸降低，巷道內(nèi)平均溫度呈明顯增加趨勢(shì)。巷道空間內(nèi)的風(fēng)流受熱上浮，靠近頂板區(qū)域溫度比靠近底板區(qū)域的溫度升高。頂板區(qū)域的熱氣流受熱向外回流，靠近底板區(qū)域的一小部分空氣受風(fēng)筒口氣流的卷吸作用沿X方向運(yùn)動(dòng)。

由于鏟運(yùn)機(jī)在工作過程中會(huì)產(chǎn)生大量的熱量，導(dǎo)致鏟運(yùn)機(jī)附近的溫度比巷道內(nèi)其他地點(diǎn)溫度較高，巷道熱環(huán)境明顯受到鏟運(yùn)機(jī)運(yùn)轉(zhuǎn)放熱的影響。因此，從通風(fēng)降溫效果來(lái)看，在一定程度內(nèi)通過提高送風(fēng)速度，能夠改善巷道熱環(huán)境，這將有利于在較低成本下改善井下工作人員的工作環(huán)境，從會(huì)提高工作人員工作效率。

由此可見，在考慮經(jīng)濟(jì)成本的前提下并考慮風(fēng)管的合理布置，適當(dāng)提高送風(fēng)速度，對(duì)巷道熱環(huán)境進(jìn)行優(yōu)化的方法是可取的。但需要注意的是，過高的送風(fēng)速度會(huì)導(dǎo)致獨(dú)頭巷道風(fēng)速超過相關(guān)安全規(guī)程規(guī)定的上限，并且由于過高的送風(fēng)速度配合鏟運(yùn)機(jī)的阻塞作用，更容易在獨(dú)頭巷道內(nèi)形成渦流，從而抑制風(fēng)流的回流，因此，持續(xù)增大送風(fēng)速度并不能達(dá)到預(yù)期的降溫效果，對(duì)獨(dú)頭巷道的熱環(huán)境優(yōu)化存在一個(gè)最優(yōu)送風(fēng)速度區(qū)間。

5 結(jié)語(yǔ)

在金屬礦山采場(chǎng)獨(dú)頭巷道，通過增加風(fēng)管出風(fēng)口風(fēng)速，可在一定程度上改善巷道熱環(huán)境，但單純依靠增大風(fēng)速進(jìn)行巷道降溫時(shí)，需要注意其存在的局限性，有必要根據(jù)實(shí)際情況探索出最優(yōu)風(fēng)速值。巷道內(nèi)存在的鏟運(yùn)機(jī)、掘進(jìn)機(jī)等通風(fēng)障礙物的阻塞作用會(huì)對(duì)獨(dú)頭巷道的風(fēng)流場(chǎng)與溫度場(chǎng)產(chǎn)生較大影響，在優(yōu)化通風(fēng)降溫工程時(shí)必須考慮鏟運(yùn)機(jī)或掘進(jìn)機(jī)對(duì)熱環(huán)境的影響。