黃瓊，丁兆云

（1.中國消防救援學(xué)院 基礎(chǔ)部，北京 102202；2.國防科技大學(xué) 系統(tǒng)工程學(xué)院，長沙 410073）

0 引言

隨著公路、鐵路和城市地鐵隧道數(shù)量的劇增，隧道火災(zāi)事故發(fā)生頻率大幅增加，成為危害最嚴(yán)重的隧道交通事故類型之一。由于隧道屬于狹長密閉空間，火災(zāi)產(chǎn)生的熱量和濃煙大量聚集，會(huì)對(duì)隧道內(nèi)人員的生命安全造成嚴(yán)重威脅。目前關(guān)于隧道火災(zāi)的煙氣運(yùn)動(dòng)和控制已得到研究人員的廣泛關(guān)注，并對(duì)此進(jìn)行了大量研究［1-3］。作為火災(zāi)煙氣運(yùn)動(dòng)的重要參數(shù)，火災(zāi)煙氣速度的準(zhǔn)確估計(jì)具有重要意義。但由于火災(zāi)的復(fù)雜性，煙霧的傳播存在諸多不確定性，高溫?zé)煔馑俣裙烙?jì)是一項(xiàng)困難且具有挑戰(zhàn)性的任務(wù)。目前主要有兩種方法：一是基于物理模型設(shè)計(jì)測速設(shè)備進(jìn)行測量估計(jì)。比如利用皮托管根據(jù)壓差完成測量［4］，但該方法不適合低流速測量，誤差較大；采用熱線式設(shè)備測量［5］的準(zhǔn)確度較高，但探頭易被煙氣中的炭黑粒子和細(xì)小液滴堵塞損壞，測量成本較高。也有許多研究人員提出了低成本且簡單可靠的測量方法［6-8］，但由于需要預(yù)先布設(shè)相關(guān)測量設(shè)備，在實(shí)際火災(zāi)應(yīng)用中受到一定限制。二是基于火災(zāi)模型進(jìn)行有效的火災(zāi)蔓延和煙氣擴(kuò)散模擬估計(jì)。目前的火災(zāi)模型有經(jīng)驗(yàn)?zāi)Ｐ?、區(qū)域模型及計(jì)算流體動(dòng)力學(xué)（Computational Fluid Dynamics，CFD）模型。CFD 模型以火焰動(dòng)力學(xué)模擬器（Fire Dynamics Simulator，F(xiàn)DS）［9］為代表，在模擬火災(zāi)時(shí)，精確度高，但對(duì)計(jì)算資源要求也高，往往需要耗費(fèi)比真實(shí)火災(zāi)更多的時(shí)間實(shí)現(xiàn)模擬。區(qū)域模型，如綜合火災(zāi)增長和煙霧傳輸?shù)腃FAST（Consolidated model of Fire And Smoke Transport）［10］，將火災(zāi)區(qū)域分為兩區(qū)，且假設(shè)各區(qū)的物理量是均勻的，通過開口計(jì)算煙氣分布及溫度。在大多數(shù)情況下，現(xiàn)有的火災(zāi)模擬工具的準(zhǔn)確性很大程度上依賴于用戶預(yù)先設(shè)定的模型參數(shù)，如火源的熱釋放率、通風(fēng)條件、火源位置等［11-14］和系統(tǒng)的初始狀態(tài)。在實(shí)際火災(zāi)中，由于火災(zāi)的復(fù)雜性，影響火災(zāi)蔓延和煙氣擴(kuò)散的參數(shù)可能隨時(shí)間變化，導(dǎo)致傳統(tǒng)的火災(zāi)模型的模擬結(jié)果與實(shí)際情況有很大不同。許多學(xué)者試圖結(jié)合模型與傳感器數(shù)據(jù)［15-17］以改進(jìn)模擬模型，但仍然難以保證實(shí)時(shí)性。近年來，數(shù)據(jù)同化技術(shù)（Data Assimilation，DA）如集合卡爾曼濾波（Ensemble Kalman Filter，EnKF）、粒子濾波器等［18-20］在許多領(lǐng)域廣泛應(yīng)用。大部分研究［21-23］基于EnKF 結(jié)合觀測數(shù)據(jù)與模型，實(shí)時(shí)模擬火災(zāi)煙氣的溫度分布和煙層高度，能夠獲得較為精確的模擬估計(jì)性能，且對(duì)計(jì)算機(jī)資源的需求顯著降低。本文選擇粒子濾波算法模擬估計(jì)隧道火災(zāi)煙氣速度，以解決目前測量成本較高、模擬精度較低與較難保證實(shí)時(shí)性等問題。實(shí)驗(yàn)結(jié)果表明，本文方法能夠基本滿足實(shí)時(shí)性要求，且具有較高的模擬精度，能為消防救援人員提供有用的關(guān)鍵信息，也可為排煙系統(tǒng)和消防規(guī)劃策略提供理論依據(jù)。

1 粒子濾波原理

火災(zāi)煙氣運(yùn)動(dòng)作為復(fù)雜的非線性動(dòng)態(tài)系統(tǒng)，受多種不確定因素的影響。粒子濾波算法采用貫序的蒙特卡洛方法，能夠?qū)Ψ蔷€性的狀態(tài)空間進(jìn)行有效估計(jì)，適用于非線性非高斯系統(tǒng)。動(dòng)態(tài)時(shí)變系統(tǒng)的狀態(tài)方程和觀測方程可以描述為：

其中：Xk是k時(shí)刻系統(tǒng)的狀態(tài)量；Zk是k時(shí)刻的系統(tǒng)觀測量；Uk和Rk分別是系統(tǒng)的狀態(tài)噪聲和觀測噪聲，表征系統(tǒng)狀態(tài)估計(jì)及觀測過程的不確定性程度；Fk和Hk分別表示狀態(tài)轉(zhuǎn)移函數(shù)和系統(tǒng)觀測方程，用于描述系統(tǒng)在狀態(tài)空間中的狀態(tài)變化過程、系統(tǒng)狀態(tài)與觀測量間的映射關(guān)系。

1.1 火災(zāi)煙氣運(yùn)動(dòng)的狀態(tài)方程

當(dāng)隧道發(fā)生火災(zāi)時(shí)，熱煙氣上升和卷吸的空氣混合形成火羽流，受到頂棚撞擊后，發(fā)生水平流動(dòng)。紀(jì)杰等［24］將隧道火災(zāi)煙氣的運(yùn)動(dòng)區(qū)域劃分為：自由上升羽流階段、頂棚徑向蔓延階段、過渡階段和一維水平蔓延階段。本文主要針對(duì)煙氣一維水平蔓延過程中的速度估計(jì)展開。已知煙氣速度在火羽流內(nèi)幾乎一致，但超出火羽流，急劇下降為0［25］。與火羽流中軸線距離r的氣流速度可表示為：

其中：Q為熱釋放速率；H為隧道頂部到火源的高度。

基于離散時(shí)間序列，將k時(shí)刻的熱釋放速率Qk泰勒展開，如式（2）所示，同時(shí)記為常數(shù)。

結(jié)合式（2）～（4），建立煙氣運(yùn)動(dòng)狀態(tài)方程：

其中：dt為k時(shí)刻和k-1 時(shí)刻的時(shí)間間隔；系統(tǒng)狀態(tài)Xk為通過添加噪聲參數(shù)U=(U0，U1，U2)T使?fàn)顟B(tài)方程能夠更好地反映煙氣運(yùn)動(dòng)過程中的不確定性，提高煙氣速度估計(jì)的可靠性和準(zhǔn)確性。

1.2 基于溫度的觀測方程

在離開火羽流中軸線的任意距離r處，溫度分布體現(xiàn)為：隨著垂直頂棚的距離增加，溫度升高，當(dāng)達(dá)到最高值后，又逐漸降至環(huán)境溫度T0。根據(jù)一系列全尺寸火災(zāi)實(shí)驗(yàn)結(jié)果，導(dǎo)出如下溫度分布關(guān)系式［25］：

將式（1）代入式（6），得到系統(tǒng)觀測方程，如式（7）所示，并記為常數(shù)。

為了更好地模擬隧道火災(zāi)環(huán)境下的觀測過程，式（7）中引入測量噪聲以表征傳感器的不確定性觀測誤差分別為正態(tài)分布的均值和方差。監(jiān)測離開火羽流中軸線任意距離處的溫度值，根據(jù)式（7）建立煙氣運(yùn)動(dòng)過程中某點(diǎn)溫度與該點(diǎn)煙氣速度之間的特征關(guān)系，以準(zhǔn)確反映煙氣特性間的變化規(guī)律。

基于上述建立的狀態(tài)方程和觀測方程，運(yùn)用粒子濾波的隧道火災(zāi)煙氣速度估計(jì)方法實(shí)現(xiàn)流程如下。

步驟6 若滿足終止條件，結(jié)束流程；否則k=k+1，轉(zhuǎn)步驟2，繼續(xù)估計(jì)k+1 時(shí)刻的煙氣速度。

2 實(shí)驗(yàn)與結(jié)果分析

2.1 隧道模型和測點(diǎn)布置

本文通過FDS 火災(zāi)模擬器獲取實(shí)驗(yàn)數(shù)據(jù)。模擬隧道長50 m，寬7 m，高5 m，隧道兩端口的屬性設(shè)置為“OPEN”，即隧道端口直接連通開放空間。通過網(wǎng)格獨(dú)立性分析，確定合適的網(wǎng)格尺寸，以保證模擬的可靠性和準(zhǔn)確性。McGrattan等［26］建議取值為［4，16］，其中：δx是網(wǎng)格尺寸；D*為火源特征長度，因此網(wǎng)格尺寸設(shè)為0.5 m。隧道頂棚和側(cè)壁的材料屬性設(shè)置為“CONCRETE”。隧道內(nèi)初始環(huán)境溫度設(shè)為20 ℃。著火點(diǎn)布置在隧道縱向中心線處。選擇方形作為模擬火源形狀，與隧道地面齊平，燃燒面尺寸為1 m×1 m，燃料屬性為正庚烷?？紤]到公路隧道中小型汽車和載貨卡車的火災(zāi)熱釋放速率（Heat Release Rate，HRR）的范圍約為3～5 MW 和10～20 MW［27］，為使模擬實(shí)驗(yàn)更貼近真實(shí)火災(zāi)，令最大熱釋放速率分別為3 MW 和15 MW。隧道內(nèi)設(shè)置了溫度和氣流速度的測點(diǎn)。由于豎直分布的溫度最大值在距離頂棚0.02H（H為隧道頂部到火源的高度）處［24］，且系統(tǒng)狀態(tài)方程和觀測方程的建立條件為至少距離火羽流中軸線0.18H外的任意處，因此為保證算法的精確性和有效性，在隧道頂棚下方10 cm，氣流速度測點(diǎn)布置在距離火源3 m 和6 m 處，溫度測點(diǎn)則在火源一側(cè)6 m 范圍內(nèi)每0.5 m 為間隔縱向排列。隧道模型和測點(diǎn)布置如圖1 所示。在該模型的基礎(chǔ)上進(jìn)行4 個(gè)場景的實(shí)驗(yàn)，如表1 所示。

圖1 隧道模型和測點(diǎn)布置Fig.1 Tunnel model and measuring point layout

表1 實(shí)驗(yàn)場景設(shè)置Tab.1 Experimental scene setting

2.2 結(jié)果分析

將隧道火災(zāi)模擬實(shí)驗(yàn)產(chǎn)生的數(shù)據(jù)，通過插值法進(jìn)行預(yù)處理，并依據(jù)隧道模型參數(shù)建立系統(tǒng)狀態(tài)方程和觀測方程，再基于粒子濾波算法實(shí)現(xiàn)煙氣速度估計(jì)驗(yàn)證與分析。

2.2.1 敏感性分析

在粒子濾波中，粒子數(shù)目越高，粒子多樣性越好，能有效延緩?fù)嘶?，但必然增加算法的時(shí)間復(fù)雜度，在實(shí)時(shí)性要求較高的應(yīng)用領(lǐng)域無法容忍。因此有必要對(duì)粒子數(shù)n進(jìn)行討論，以確定n對(duì)煙氣速度估計(jì)的影響。在場景A1 下，粒子數(shù)分別為100、120，150、180、200 時(shí)，煙氣速度估計(jì)的精度及計(jì)算時(shí)間，結(jié)果如圖2、3 所示。

從圖2 可看出，煙氣速度估計(jì)曲線在n=150，180，200 時(shí)，基本趨于一致，接近于真實(shí)狀態(tài)，而n=100，120 時(shí)，估計(jì)狀態(tài)與真實(shí)狀態(tài)偏離較遠(yuǎn)。這一點(diǎn)同樣可以在圖3 中清楚地反映出來，實(shí)驗(yàn)中利用均方根誤差（Root Mean Square Error，RMSE）衡量估計(jì)值與真實(shí)值間的偏差。RMSE 越低，估計(jì)精度越高?？梢钥闯鰊從100 增加到150，RMSE 急劇下降，而150 增加到200，RMSE 基本相當(dāng)，均處于較低值。若總模擬時(shí)間為3 000 s，每10 s 獲取一次觀測數(shù)據(jù)，為滿足實(shí)時(shí)性要求，每次獲取觀測數(shù)據(jù)后，將響應(yīng)時(shí)間控制在10 ms 以內(nèi)，即在本實(shí)驗(yàn)中的計(jì)算時(shí)間應(yīng)小于3 s，是可接受的。由表2 可知，當(dāng)n=180，200 時(shí)，計(jì)算時(shí)間超過3 s；當(dāng)n=100，120，150時(shí)，計(jì)算時(shí)間少于3 s，可在每次獲取觀測數(shù)據(jù)后進(jìn)行實(shí)時(shí)估計(jì)。綜合考慮估計(jì)精度與計(jì)算時(shí)間，當(dāng)n=150 時(shí)，實(shí)時(shí)估計(jì)隧道火災(zāi)煙氣速度的效果最佳。

圖2 粒子數(shù)對(duì)估計(jì)精度的影響Fig.2 Influence of particle number on estimation accuracy

圖3 不同煙氣速度的RMSE對(duì)比Fig.3 Comparison of RMSE for different smoke velocities

表2 粒子數(shù)對(duì)計(jì)算時(shí)間的影響Tab.2 Influence of particle number on calculation time

2.2.2 觀測數(shù)據(jù)分析

本節(jié)從兩個(gè)方面研究觀測數(shù)據(jù)對(duì)煙氣速度估計(jì)效果的影響如圖4 所示。

圖4 觀測數(shù)據(jù)對(duì)估計(jì)精度的影響Fig.4 Influence of observation data on estimation accuracy

1）由于觀測位置的偏差，必須對(duì)觀測數(shù)據(jù)進(jìn)行近似處理。如在場景A2 和A3 下，溫度Ts1分別根據(jù)Ts4和Ts5、Ts3和Ts6通過線性插值方法獲得，討論近似誤差對(duì)Vs1估計(jì)精度的影響，如圖4（a）所示。

2）觀測數(shù)據(jù)的不確定性對(duì)估計(jì)結(jié)果的影響至關(guān)重要，主要來源于傳感器的測量誤差。因此為解釋傳感器的內(nèi)在特性對(duì)估計(jì)效果的影響，在場景A4 下，往觀測數(shù)據(jù)中分別加入?yún)f(xié)方差為Rt1=10和Rt2=30 的噪聲，這是影響粒子濾波效果的一個(gè)重要參數(shù)，如圖4（b）所示。

圖4（a）的結(jié)果表明，總體上估計(jì)精度能夠滿足要求，尤其在火災(zāi)穩(wěn)定燃燒階段；觀測點(diǎn)位置距離估計(jì)點(diǎn)越近，估計(jì)值越接近真實(shí)狀態(tài)值。主要原因在于觀測數(shù)據(jù)的誤差削弱了它在模擬估計(jì)過程中對(duì)系統(tǒng)狀態(tài)的修正作用。下一步將研究更為有效的近似方法，以提高觀測數(shù)據(jù)質(zhì)量，保證較好的估計(jì)效果。從圖4（b）可看出，觀測不確定性程度越小，估計(jì)效果越好。由于觀測器件的內(nèi)在特性，遠(yuǎn)離火源的傳感器的觀測誤差相較于火源附近的探測器的觀測誤差更小。因此以后的研究將對(duì)火源附近的估計(jì)點(diǎn)獲取的觀測數(shù)據(jù)進(jìn)行有效合理的降噪處理，以降低它對(duì)估計(jì)精度的影響。

2.2.3 煙氣速度估計(jì)效果測試

在A1 和A5 場景下實(shí)現(xiàn)煙氣速度估計(jì)測試，如圖5 所示?？梢钥闯?，在t=250 s和t=560 s 之前，估計(jì)誤差相對(duì)較大，之后煙氣速度的估計(jì)曲線與真實(shí)狀態(tài)隨時(shí)間的變化趨勢基本一致，估計(jì)值在真實(shí)煙氣速度的±20%誤差范圍內(nèi)。使用平均絕對(duì)誤差（Mean Absolute Error，MAE）進(jìn)行評(píng)估，如表3 所示，各場景下煙氣速度的MAE 均低于真實(shí)狀態(tài)值的20%（ts0.2）?？傮w模擬估計(jì)效果較好，而且誤差基本能夠控制在真實(shí)值的20%以內(nèi)。因?yàn)橛^測方程基于穩(wěn)態(tài)方程建立，火災(zāi)燃燒的初始階段為非穩(wěn)態(tài)過程，在一定程度上影響了估計(jì)精度。為解決這一問題，可以通過調(diào)整觀測方程中的噪聲參數(shù)以減少對(duì)系統(tǒng)狀態(tài)修正的影響。

表3 A1與A5場景下的MAE與真實(shí)狀態(tài)值20%的平均值的對(duì)比 單位：m/sTab.3 Comparison of MAE and 20% of true state value in scenario A1 and A5 unit：m/s

圖5 不同場景下的煙氣速度估計(jì)結(jié)果Fig.5 Smoke velocity estimation results in different scenarios

3 結(jié)語

本文提出一種基于粒子濾波算法的隧道火災(zāi)煙氣速度估計(jì)方法，實(shí)驗(yàn)結(jié)果表明，該方法可通過傳感器的觀測值修正實(shí)現(xiàn)煙氣速度的精確估計(jì)。同時(shí)討論了影響模擬估計(jì)精度的主要因素，主要結(jié)論如下：1）粒子數(shù)目反映了粒子多樣性程度，粒子數(shù)越高，煙氣速度估計(jì)的精確度越高，但帶來的是計(jì)算時(shí)間的大幅增加，無法滿足實(shí)時(shí)性要求；2）觀測位置與觀測的不確定性同樣也會(huì)影響估計(jì)精度。觀測位置越接近估計(jì)點(diǎn)，觀測的不確定程度越低，估計(jì)性能越好，因此可對(duì)觀測數(shù)據(jù)進(jìn)行預(yù)處理，以提高模擬估計(jì)效果。3）在火災(zāi)穩(wěn)定燃燒階段的估計(jì)效果優(yōu)于燃燒初期，但可通過調(diào)整方程中的噪聲參數(shù)，使模擬估計(jì)效果達(dá)到較好水平。本文方法可用于一般隧道火災(zāi)場景，但由于實(shí)驗(yàn)環(huán)境建立在FDS 模擬基礎(chǔ)上，為進(jìn)一步驗(yàn)證估計(jì)方法在實(shí)際火災(zāi)中的模擬效果，需要再進(jìn)行小尺寸或全尺寸的火災(zāi)實(shí)驗(yàn)，將模擬估計(jì)值與真實(shí)火災(zāi)狀態(tài)值進(jìn)行比較，以應(yīng)用于實(shí)際工作，為消防救援和人員疏散提供有用的關(guān)鍵信息，同時(shí)也為排煙系統(tǒng)和消防規(guī)劃策略提供理論依據(jù)。