不同封堵條件下巷道火災煙氣溫度場預測模型

馬礪 張鵬宇 李超華 劉順

摘 要：巷道火災蔓延速度快且溫度高，封堵控制過程中會對工作人員造成人身威脅，因此研究巷道火災煙氣封堵過程中溫度變化規(guī)律及預測對其安全控制有重要意義。文中通過分析煙氣質(zhì)量流率的卷吸特征，基于火災煙氣溫度控制方程建立了封堵條件下的巷道火災煙氣溫度預測模型。搭建1

∶10比例的巷道火災相似模擬實驗臺，開展火源功率為4，6，8 MW的巷道火災實驗。測試進風側單側封堵比例為0%，25%，50%，75%，100%的15種工況下頂棚煙氣溫度分布，得到了不同封堵條件下的煙氣溫度衰減規(guī)律。采用平均絕對誤差（MAE）、平均絕對百分誤差（MAPE）、均方根誤差（RMSE）和確定系數(shù)（R2）驗證實驗值與理論預測值的準確性。結果表明：當封堵比例大于50%時，隨著封堵比例增加，抑制火災效果明顯，在不同封堵下煙氣縱向溫度分布呈高斯函數(shù)規(guī)律;通過對比不同封堵條件下的溫度場的預測結果（MAE<0.075 1，MAPE<0132 1，MSE<0.077，R2=0.906 9），發(fā)現(xiàn)該模型對封閉過程火區(qū)溫度的預測結果較好。關鍵詞：安全科學與工程;巷道火災;溫度場;封堵;預測模型中圖分類號：TD 75+2

文獻標志碼：A

文章編號：1672-9315（2021）02-0213-08

DOI：10.13800/j.cnki.xakjdxxb.2021.0204開放科學（資源服務）標識碼（OSID）：

Prediction model of roadway fire temperature field under

different plugging conditions

MA Li1，2，ZHANG Pengyu1，2，LI Chaohua1，2，LIU Shun1，2

（1.College of Safety Science and Engineering，Xian University of Science and Technology，Xian 710054，China;

2.Shaanxi Key Laboratory of Prevention and Control of Coal Fire，Xian University of Science and Technology，Xian 710054，China）

Abstract：The roadway fire spreads sharply with the high temperature，which can threat the safety of staff during the plugging control.Therefore，it is of great significance to explore the temperature change and prediction in the process of roadway fire smoke plugging for its safety control.In this study，the entrainment characteristics of the mass flow rate of the flue gas were analyzed，and a prediction model of the flue gas temperature of the roadway fire under the plugging condition was established based on the control equation of the fire smoke temperature.A similar simulation test bench for the 1∶10 ratio of roadway fires was set up，and the roadway fire experiments with the fire source power of 4，6，8 MW were tested.The roof smoke temperatures under 15 conditions were tested on the inlet side at 0%，25%，50%，75%，and 100%，and the flue gas temperature variations under different plugging conditions were analyzed.In addition，the average absolute error（MAE），average absolute percentage error（MAPE），root mean square error（RMSE）and determination coefficient（R2）were applied to verify the similarities between the values of experimental and theories.The results indicated that：when the blocking ratio was greater than 50%，the fire suppression effect was remarkable with the increase of blocking ratio，and the vertical temperature distribution of the flue gas under different blocking ratios described a Gaussian function variation.By comparing the prediction results of the temperature field

under different plugging conditions（MAE<0.075 1，MAPE<0.132 1，MSE<0.077，and R2=0.906 9），the model turned to have a better prediction result for the fire zone temperature in the blocking process.

Key words：safety science and engineering;roadway fire;temperature field;plugging;predictive model

0 引 言

在火災的發(fā)展蔓延階段，高溫煙氣是阻礙人員疏散，造成人員傷亡的重要因素，研究煙氣蔓延的規(guī)律并準確預測煙氣溫度對火災的預防控制及人員的疏散救援有著重要意義[1-2]。巷道火災按照時間順序一般可劃分為3個階段，即火災初期增長階段、火災充分發(fā)展階段、火災衰減階段。在火災充分發(fā)展階段，根據(jù)火災煙氣蔓延規(guī)律分為5個階段[3]：火羽流上升階段、頂棚射流階段，二維徑向擴散階段、一維縱向蔓延轉化階段以及一維縱向蔓延階段[4-5]。對于火災煙氣溫度場的研究，MODIC利用IDS通風模擬軟件對公路隧道火災中的氣流溫度和速度進行研究，指導了消防安全疏散[6];HWANG等通過模擬得到不同通風條件下煙氣溫度和速度的分布規(guī)律以及對火災的影響[7];LEE等模擬了不同寬高比下的隧道火災煙氣溫度場、速度場等影響規(guī)律[8];BAILEY等通過LES大渦數(shù)值模擬得到了走廊煙氣溫度的縱向分布規(guī)律[9];

KHATTRI采用小尺寸數(shù)值模擬研究火災煙氣溫度分布規(guī)律，得到煙氣溫度的衰減并不一定是指數(shù)函數(shù)分布規(guī)律[10]。郭慶華等通過對海拔4 100 m的全尺寸隧道進行火災實驗，得到3種火源功率下最高溫度及縱向溫度分布規(guī)律[11];陶海軍等建立了1∶36小尺寸綜合管廊模型的火災實驗，得到煙氣溫度在管廊電力艙室的衰減規(guī)律，并分析煙氣溫度在不同高度下的跳躍現(xiàn)象[12];張鵬等基于1∶20地鐵換乘站相似實驗模型，利用高溫輻射流取代火源，得到煙氣溫度在頂棚沿站臺方向和樓梯處衰減規(guī)律[13];梁震寰等搭建圓形小尺寸管廊實驗模型，研究得到4種火源形狀的距離、煙氣溫度以及頂棚溫度的函數(shù)關系[14];劉斌等對復雜結構隧道火災進行數(shù)值模擬，分析得到不同縱向通風對煙氣溫度的分布規(guī)律及煙流的擴散影響[15];姜學鵬等通過FDS模擬設置不同坡度的人字坡隧道模型，得到最高煙氣溫度與火源功率呈現(xiàn)3/4次冪的關系且沿隧道長度方向呈指數(shù)函數(shù)衰減[16];姜濤利用CFD方法設計不同工況下的地鐵過江隧道模型，得到火災煙氣溫度場及濃度場，為救援工作提供技術指導[17]。對于火災封堵控制研究，

YAO等通過建立小尺寸實驗臺研究火災封堵窒息，結果表明完全密封的時間越早，燃燒時間越短，溫度越低，抑制火災效果越好[18-19];紀道溪通過數(shù)值模擬研究封堵控制因素對滅火效率影響規(guī)律[20];董炳燕等通過數(shù)值模擬得出火源功率為50 MW的臨界封堵比例為75%，火源功率為100，150 MW的臨界封堵比例為50%[21].

綜上所述，當前對于火災煙氣溫度的研究中，多數(shù)學者通過建立小尺寸實驗模型來研究煙氣溫度的分布影響規(guī)律，對于封堵過程研究的模型較少。文中主要研究在巷道火災封閉過程中煙氣一維蔓延階段溫度預測模型，為井下封閉火區(qū)的蔓延控制提供理論基礎。

1 實驗裝置及條件

1.1 實驗系統(tǒng)實驗系統(tǒng)利用西安科技大學煤火災害重點實驗室巷道裝置模型[22]，建立了1

∶10比例的巷道火災發(fā)展封堵過程實驗模型。巷道模型頂部為拱形，底部為矩形，如圖1所示。該系統(tǒng)采用分段式設計，單節(jié)長度0.7 m，共11節(jié)，巷道模型截面寬度0.48 m，高度0.42 m，為便于觀察巷道火災煙氣溫度的變化趨勢，巷道一面采用耐高溫玻璃制成，另一面采用角鋼制成。

實驗臺采用自然通風，風速以實際測量為準，本實驗為0.3 m/s，火源距進風側250 cm。實驗布

置19個K型鎧裝熱電偶，且相對火源中心對稱布置，分別距火源中心點距離為12.5，25，56.5，75，93.5，140.5，159，199.5和240 cm.火源及熱電偶分布如圖2所示，采用C7108溫度數(shù)據(jù)模塊采集溫度。

1.2 實驗樣品為了模擬井下巷道火災，火源大小設計為4，6，8 MW.火源類型為柴油，熱值為42 kJ/g，火源相似計算值見表1。

2 實驗數(shù)據(jù)分析

2.1 頂棚最高煙氣溫度分析3種火源功率下頂棚煙氣溫度封堵曲線如圖3所示。結果表明：隨著火源功率增加，封堵比例升高頂棚煙氣溫度增加，這是由于在巷道內(nèi)燃燒消耗的氧氣增多，造成熱量聚積。0～50%封堵火災充分燃燒時間和溫度情況相近，50%封堵的燃燒時間更長溫度更高。表明在封堵比例小于50%時，隨著新鮮風流進入巷道，氧氣的供給足夠維持柴油的充分燃燒，使得熱量難以聚積，但由于火災充分燃燒時間長，不利于人員的疏散。75%封堵比例時，火災的充分燃燒時間相比0～50%封堵明顯縮短，且隨著火源功率的增加，達到充分燃燒階段的時間提前，燃燒所需的時間隨之減少，且衰減階段溫度下降較快，表明75%封堵比例對氧氣的供給造成了極大阻礙，火源功率越高消耗氧氣越多，供氧不足，火災得到了明顯的抑制。完全封堵（比例100%）時火災充分燃燒階段最短，溫度上升和衰減的幅度也最大，這是由于阻斷了火災燃燒消耗所需的氧氣來源，導致頂棚熱量聚集，燃燒不充分，火災熄滅的時間快。圖4表示的是在熄滅時，各封堵比例的縱向最高煙氣溫度變化曲線，從圖中可以看出隨著封堵比例升高，溫度衰減的幅度越大;綜合圖3圖4可知，當封堵比例大于50%時，在距火源240 cm處的下風側各封堵比例的溫度相近，當時間到達500 s時，100%封堵的溫度最低，75%封堵的溫度與0～25%的溫度接近。綜上所述，由于100%封堵導致熱量聚集過快，溫度快速上升至600 ℃，對巷道頂棚破壞嚴重，考慮到保護救援人員，當封堵50%以上時，隔氧滅火效果明顯。

2.2 不同封堵條件對縱向煙氣溫度的衰減規(guī)律KUNSCH通過理論分析和實驗研究得到了頂棚煙氣溫度縱向分布模型[23]，經(jīng)過國內(nèi)外學者驗證該模型服從指數(shù)衰減規(guī)律，如公式（1）所示。文中取距火源中心點距離x-x0為自變量，煙氣溫ΔT/ΔT0為因變量，利用該模型對實驗數(shù)據(jù)進行擬合，擬合后各參數(shù)數(shù)據(jù)見表2。

ΔTΔT0

=Ta-TTa-T0

=a·exp（-b（x-x0））

（1）

分析表中擬合數(shù)據(jù)得到封堵條件對參數(shù)a和b的影響。

1）當封堵比例η=0時，煙氣溫度與距離的關系

ΔTΔT0

=

1.1·exp（-0.006 7（x-x0））

（2）2）當封堵比例η>0時，煙氣溫度會隨著封堵比例呈線性變化

ΔTΔT0

=

1·exp（-（0.003 932+0.004 246）（x-x0））

（3）

3 巷道火災封堵煙氣溫度預測模型3.1 巷道火災一維溫度場分布通過分析火災煙氣溫度場控制方程和質(zhì)量流率的卷吸方程[3]，建立符合礦井巷道條件下的一維溫度場分布模型。該模型增加了封堵條件，巷道截面如圖5所示。

邊界條件：x=0時，

=0;x=x0時，T=T0。假設條件：在一維蔓延階段，巷道火災的煙氣層在橫向上溫度一致;煙氣為不可壓縮流體。巷道火災的煙氣層在x方向質(zhì)量流率變化，等

于煙氣層卷吸空氣質(zhì)量流率變化，得到連續(xù)性方程

ddx=w

（4）巷道火災煙氣層在x方向上的焓值變化，等于煙氣層卷吸空氣的焓值增量與對流輻射損失的焓值量總和，得到能量方程

d（cPT）

dx=cPwTa+

（5）其中Kunsch提出一維水平蔓延的煙氣卷吸的質(zhì)量流率可以表示為[23]

w=βρawΔux

（6）巷道火災煙氣的能量損失可以看成是煙氣層與環(huán)境溫度之間的能量交換（忽略壁面厚度）

=αD（Ta-T）

（7）

在一維蔓延階段，卷吸系數(shù)β的值很小，在數(shù)量級上可以趨近于0

βρaΔuρs

x→0

（8）為了建立封堵比例η與空氣流速u的關系，文中簡化在封堵情況下的封堵兩側的體積流量變化，即自然通風風量減去進入巷道風量等于被封堵處風量的損失值

Q進-Q=Q損

（9）損失值可以用進入巷道風流體積流量減去開口處進入風流體積流量

Q損=uA-uA1

，帶入化簡可以得到

u=um1+η

（10）綜合以上因素，聯(lián)立化簡方程，根據(jù)參考點位置x0，T0可以積分得到

當封堵比例η=01）當h≤R時

ΔTΔT0

=exp

-

απR180

arccos

R-h0R

cPwρsum

（x-x0）

（11）

2）當h≥R時

ΔTΔT0

=exp

-

απR180

arccos

h0-RR+2α（h0-R）

cPwρsum

（x-x0）

（12）

當封堵比例η>01）當h≤R時

ΔTΔT0

=exp

-

α（1+η）πR180

arccos

R-h0R

cPwρsum

（x-x0）

（13）

2）當h≥R時

ΔTΔT0

=

exp

-

α（1+η）πR180

arccos

R-h0R+2α（1+η）（h0-R）

cPwρsum

（x-x0）

（14）

其中 為巷道火災煙氣層質(zhì)量流率，kg/s;w為煙氣層卷吸空氣的質(zhì)量流率，kg/s;cp為比熱容，kJ/（kg·K）;

T為煙氣層溫度，℃;

Ta為環(huán)境溫度，℃;為熱交換損失的能量流率，kW;

β為水平卷吸系數(shù);ρa為空氣密度，kg/m3;w為煙氣層與空氣之間的分層寬度，m;

Δu為煙氣層和空氣層的相對速度，m/s;h為煙氣層厚度，m;ρs為煙氣層密度，kg/m3;h0為穩(wěn)定階段初始煙氣層高度，m;R為

圓弧拱頂半徑，m;h1為圓弧拱頂截面的最高高度，m。

3.2 理論模型預測結果根據(jù)巷道火災實驗臺的實際尺寸和環(huán)境條件，確定以下基礎數(shù)據(jù)：表面?zhèn)鳠嵯禂?shù)

α取45 W/（m·K）;煙氣比熱容cp為1.02 kJ/（kg·K）;穩(wěn)定階段煙氣初始質(zhì)量流量根據(jù)煙氣密度ps，乘以煙氣層與空氣之間的分層寬度w來表示=wρsu。得到封堵條件下巷道火災一維煙氣溫度預測值。采用平均絕對誤差（MAE）、平均絕對百分誤差（MAPE）、根均方誤差（RMSE）和確定系數(shù)（R2）等指標來評價預測模型的準確性[24-25]。

MAE=1nni=1

|fi-yi|

（15）

MAPE=1nni=1

fi-yiyi×100%

（16）

RMSE=1nni=1

（fi-yi）2

（17）

R2=1-ni=1

（fi-yi）2

ni=1

（-yi）2

（18）

式中 fi為預測值;yi為實際值;為實際值的平均值。

3.3理論模型的預測值與實驗數(shù)據(jù)的對比圖6表示的是5種封堵條件下實驗數(shù)據(jù)擬合結果和理論模型預測結果，從圖中可以看出兩者曲線相近，計算得到各曲線之間相關系數(shù)R2均在0.9以上，表明實驗值與預測值存在較強的相關性，理論模型預測值較好。文中將不同封堵條件下理論值與實驗值進行對比，為了量化模型的預測性能，將其誤差指標（MAE<0.075 1，MAPE<0.132 1，RMSE<0.077，R2>0.906 9）列于表3中。結果表明，封堵條件下的巷道火災煙氣理論模型能夠準確預測巷道火災煙氣溫度。

4 結 論

1）根據(jù)封堵過程中巷道內(nèi)外體積流量變化的平衡條件，得到流速與封堵比例之間的函數(shù)關系，基于巷道火災實驗臺建立了不同封堵條件下的煙氣溫度預測模型。

2）通過分析不同封堵比例下巷道火災煙氣溫度的預測模型，得到溫度的衰減主要受巷道尺寸、空氣流速和煙流質(zhì)量流率影響。

3）實驗分析得到頂棚煙氣溫度隨著火源功率的增加而升高;當封堵比例為50%時，火災充分發(fā)展時間最長;當封堵比例大于50%時，隨著封堵比例增加，煙氣溫度下降的時間越早。

參考文獻（References）：

[1]

郭軍，劉蔭，金永飛，等.礦井膠帶火災巷道環(huán)境多參數(shù)時空演化規(guī)律[J].西安科技大學學報，2019，39（1）：21-27.GUO Jun，LIU Yin，JIN Yongfei，et al.Spatiotemporal evolution law of multiple parameters of roadway environment for rubberbelt fire in mine[J].Journal of Xian University of Science and Technology，2019，39（1）：21-27.

[2]馬礪，范晶，任立峰，等.地鐵站臺不同排煙模式下的煙氣流動數(shù)值模擬[J].西安科技大學學報，2017，37（5）：643-648.MA Li，F(xiàn)AN Jing，REN Lifeng，et al.Numerical simulation of smoke flow on subway platform under different extraction modes[J].Journal of Xian University of Science and Technology，2017，37（5）：643-648.

[3]胡隆華.隧道火災煙氣蔓延的熱物理特性研究[D].合肥：中國科學技術大學，2006.HU Longhua.Studies on thermal physics of smoke movement in tunnel fires[D].Hefei：University of Science and Technology of China，2006.

[4]陽東，胡隆華，霍然，等.縱向風對通道火災煙氣豎向分層特性的影響[J].燃燒科學與技術，2010，16（3）：252-256.YANG Dong，HU Longhua，HUO Ran，et al.Effects of longitudinal air flow on smoke stratification in channel fire[J].Journal of Combustion Science and Technology，2010，16（3）：252-256.

[5]HU L H，HUO R，LI Y Z，et al.Fullscale burning tests on studying smoke temperature and velocity along a corridor[J].Tunnelling and Underground Space Technonlogy，2005，20（3）：233-229.

[6]MODIC J.Fire simulation in road tunnels[J].Tunnelling and Underground Space Technology，2003，18（5）：525-530.

[7]HWANG C C，EDWARDS J C.The critical ventilation velocity in tunnel firesacomputer simulation[J].Fire Safety Journal，2005，40：213-244.

[8]LEE S R，RYOU H S.A numerical study on smoke movement in longitudinal ventilation tunnel fires for different aspect ratio[J].Building and Environment，2006，40（6）：719-725.

[9]BAILEY J L，F(xiàn)ORNEY G P，TATEM P A，et al.Development and validation of corridor flow submodel for CFAST[J].Journal of Fire Protection Engineering，2002，12（3）：139-161.

[10]KHATTRI S K.From smallscale tunnel fire simulations to predicting fire dynamics in realistic tunnels[J].Tunnelling and Underground Space Technology，2017，61：198-204.

[11]郭慶華，閆治國，朱合華.高海拔隧道全尺寸火災煙氣及溫度場特征試驗研究[J].土木工程學報，2017，50（8）：114-120.GUO Qinghua，YAN Zhiguo，ZHU Hehua.Fullscale experimental study on smoke and temperature distribution characteristics of the tunnel fire at high altitude[J].China Civil Engineering Journal，2017，50（8）：114-120.

[12]陶海軍，朱國慶，趙永昌，等.綜合管廊電纜燃燒煙氣溫度實驗研究[J].消防科學與技術，2017，36（6）：745-748.TAO Haijun，ZHU Guoqing，ZHAO Yongchang，et al.Experimental study on fire smoke temperature of urban utility tunnel cable[J].Fire Science and Technology，2017，36（6）：745-748.

[13]張鵬，李嘉欣，李雪進，等.基于相似分析的地鐵換乘站縮尺寸火災實驗研究[J].消防科學與技術，2018，37（11）：1457-1460.

ZHANG Peng，LI Jiaxin，LI Xuejin，et al.Experiment on scaled fire of metro transfer station based on similarity analysis[J].Fire Science and Technology，2018，37（11）：1457-1460.

[14]梁震寰，朱國慶，田承飛.綜合管廊內(nèi)線型火源煙氣溫度場研究[J].消防科學與技術，2018，37（12）：1635-1637，1641.LIANG Zhenhuan，ZHU Guoqing，TIAN Chengfei.Research on temperature field of linear fire source in integrated pipe gallery[J].Fire Science and Technology，2018，37（12）：1635-1637，1641.

[15]劉斌，安軍，郭汝杰，等.復雜結構隧道火災煙氣逆流與溫度分布的數(shù)值模擬[J].安全與環(huán)境學報，2019，19（3）：818-823.LIU Bin，AN Jun，GUO Rujie，et al.Numerical simulation study on the firesmoke reverse flow and tempera ture distribution in the complicated structure tunnel[J].Journal of Safety and Environment，2019，19（3）：818-823.

[16]姜學鵬，廖湘娟，何振華.人字坡山嶺隧道火災頂板下方煙氣最高溫度的研究[J].安全與環(huán)境學報，2018，18（3）：925-929.JIANG Xuepeng，LIAO Xiangjuan，HE Zhenhua.Maximum temperature of smoke beneath the ceiling in the herringbone slope tunnel fire[J].Journal of Safety and Environment，2018，18（3）：925-929.

[17]姜濤.地鐵過江隧道火災CFD數(shù)值模擬研究[J].湖北大學學報（自然科學版），2015，37（5）：477-483.JIANG Tao.CFD simulation study of fire in a subway tunnel under Yangtze river[J].Journal of Hubei University（Natural Science），2015，37（5）：477-483.

[18]YAO Y，CHENG X，SHI L，et al.Experimental study on the effects of initial sealing time on fire behaviors in channel fires[J].International Journal of Thermal Sciences，2018，125：273-282.

[19]YAO Y，HE K，PENG M，et al.Maximum gas temperature rise beneath the ceiling in a portalssealed tunnel fire[J].Tunnelling and Underground Space Technology，2018，80：10-15.

[20]紀道溪.重載鐵路隧道極富燃料火災特性及封堵控制研究[D].長沙：中南大學，2012.

JI Daoxi.Reserch on fire characteristics and control by sealing portal in heavy haul railway tunnel with extremely abundant fuel[D].Changsha：Central South University，2012.

[21]董炳燕，吳晉湘，黃有波.封堵比例對隧道火災溫度影響的數(shù)值模擬研究[J].消防科學與技術，2018，37（8）：1051-1054.DONG Bingyan，WU Jinxiang，HUANG Youbo.Numerical investigation on the effect of tunnel portal sealing ratio on gas temperature during tunnel fire[J].Fire Science and Technology，2018，37（8）：1051-1054.

[22]李士戎.移動火源對隧道溫度場分布及煙氣流動影響規(guī)律研究[D].西安：西安科技大學，2013.

LI Shirong.Research on influence of moving fire on temperature field distribution and smoke spread in tunnel[D].Xian：Xian University of Science and Technology，2013.

[23]KUNSCH J P.Simple model for control of fire gases in a ventilated tunnel[J].Fire Safety Journal，2002，37（1）：67-81.

[24]LEI C K，DENG J，CAO K，et al.A random forest approach for predicting coal spontaneous combustion[J].Fuel，2018，233：63-67.

[25]劉倩穎，阮應君，時翔，等.基于Kmeans聚類與BP神經(jīng)網(wǎng)絡算法的辦公建筑逐時電負荷預測[J].熱能動力工程，2018，33（3）：138-144.

LIU Qianying，RUAN Yingjun，SHI Xiang，et al.Hourly power load forecasting of office building based on Kmeans clustering and BP neural network algorithm[J].Journal of Engineering for Thermal Energy and Power，2018，33（3）：138-144.