基于粒子濾波的隧道火災(zāi)煙氣速度估計(jì)方法

黃瓊，丁兆云

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

0 引言

隨著公路、鐵路和城市地鐵隧道數(shù)量的劇增，隧道火災(zāi)事故發(fā)生頻率大幅增加，成為危害最嚴(yán)重的隧道交通事故類(lèi)型之一。由于隧道屬于狹長(zhǎng)密閉空間，火災(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ì)測(cè)速設(shè)備進(jìn)行測(cè)量估計(jì)。比如利用皮托管根據(jù)壓差完成測(cè)量［4］，但該方法不適合低流速測(cè)量，誤差較大；采用熱線式設(shè)備測(cè)量［5］的準(zhǔn)確度較高，但探頭易被煙氣中的炭黑粒子和細(xì)小液滴堵塞損壞，測(cè)量成本較高。也有許多研究人員提出了低成本且簡(jiǎn)單可靠的測(cè)量方法［6-8］，但由于需要預(yù)先布設(shè)相關(guān)測(cè)量設(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)增長(zhǎng)和煙霧傳輸?shù)腃FAST（Consolidated model of Fire And Smoke Transport）［10］，將火災(zāi)區(qū)域分為兩區(qū)，且假設(shè)各區(qū)的物理量是均勻的，通過(guò)開(kāi)口計(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í)性。近年來(lái)，數(shù)據(jù)同化技術(shù)（Data Assimilation，DA）如集合卡爾曼濾波（Ensemble Kalman Filter，EnKF）、粒子濾波器等［18-20］在許多領(lǐng)域廣泛應(yīng)用。大部分研究［21-23］基于EnKF 結(jié)合觀測(cè)數(shù)據(jù)與模型，實(shí)時(shí)模擬火災(zāi)煙氣的溫度分布和煙層高度，能夠獲得較為精確的模擬估計(jì)性能，且對(duì)計(jì)算機(jī)資源的需求顯著降低。本文選擇粒子濾波算法模擬估計(jì)隧道火災(zāi)煙氣速度，以解決目前測(cè)量成本較高、模擬精度較低與較難保證實(shí)時(shí)性等問(wèn)題。實(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)方程和觀測(cè)方程可以描述為：

其中：Xk是k時(shí)刻系統(tǒng)的狀態(tài)量；Zk是k時(shí)刻的系統(tǒng)觀測(cè)量；Uk和Rk分別是系統(tǒng)的狀態(tài)噪聲和觀測(cè)噪聲，表征系統(tǒng)狀態(tài)估計(jì)及觀測(cè)過(guò)程的不確定性程度；Fk和Hk分別表示狀態(tài)轉(zhuǎn)移函數(shù)和系統(tǒng)觀測(cè)方程，用于描述系統(tǒng)在狀態(tài)空間中的狀態(tài)變化過(guò)程、系統(tǒng)狀態(tài)與觀測(cè)量間的映射關(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ū)域劃分為：自由上升羽流階段、頂棚徑向蔓延階段、過(guò)渡階段和一維水平蔓延階段。本文主要針對(duì)煙氣一維水平蔓延過(guò)程中的速度估計(jì)展開(kāi)。已知煙氣速度在火羽流內(nèi)幾乎一致，但超出火羽流，急劇下降為0［25］。與火羽流中軸線距離r的氣流速度可表示為：

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

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

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

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

1.2 基于溫度的觀測(cè)方程

在離開(kāi)火羽流中軸線的任意距離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)觀測(cè)方程，如式（7）所示，并記為常數(shù)。

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

基于上述建立的狀態(tài)方程和觀測(cè)方程，運(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 隧道模型和測(cè)點(diǎn)布置

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

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

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

2.2 結(jié)果分析

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

2.2.1 敏感性分析

在粒子濾波中，粒子數(shù)目越高，粒子多樣性越好，能有效延緩?fù)嘶?，但必然增加算法的時(shí)間復(fù)雜度，在實(shí)時(shí)性要求較高的應(yīng)用領(lǐng)域無(wú)法容忍。因此有必要對(duì)粒子數(shù)n進(jìn)行討論，以確定n對(duì)煙氣速度估計(jì)的影響。在場(chǎng)景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 中清楚地反映出來(lái)，實(shí)驗(yàn)中利用均方根誤差（Root Mean Square Error，RMSE）衡量估計(jì)值與真實(shí)值間的偏差。RMSE 越低，估計(jì)精度越高。可以看出n從100 增加到150，RMSE 急劇下降，而150 增加到200，RMSE 基本相當(dāng)，均處于較低值。若總模擬時(shí)間為3 000 s，每10 s 獲取一次觀測(cè)數(shù)據(jù)，為滿足實(shí)時(shí)性要求，每次獲取觀測(cè)數(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í)間超過(guò)3 s；當(dāng)n=100，120，150時(shí)，計(jì)算時(shí)間少于3 s，可在每次獲取觀測(cè)數(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 觀測(cè)數(shù)據(jù)分析

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

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

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

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

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

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

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

表3 A1與A5場(chǎng)景下的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 不同場(chǎng)景下的煙氣速度估計(jì)結(jié)果Fig.5 Smoke velocity estimation results in different scenarios

3 結(jié)語(yǔ)

本文提出一種基于粒子濾波算法的隧道火災(zāi)煙氣速度估計(jì)方法，實(shí)驗(yàn)結(jié)果表明，該方法可通過(guò)傳感器的觀測(cè)值修正實(shí)現(xiàn)煙氣速度的精確估計(jì)。同時(shí)討論了影響模擬估計(jì)精度的主要因素，主要結(jié)論如下：1）粒子數(shù)目反映了粒子多樣性程度，粒子數(shù)越高，煙氣速度估計(jì)的精確度越高，但帶來(lái)的是計(jì)算時(shí)間的大幅增加，無(wú)法滿足實(shí)時(shí)性要求；2）觀測(cè)位置與觀測(cè)的不確定性同樣也會(huì)影響估計(jì)精度。觀測(cè)位置越接近估計(jì)點(diǎn)，觀測(cè)的不確定程度越低，估計(jì)性能越好，因此可對(duì)觀測(cè)數(shù)據(jù)進(jìn)行預(yù)處理，以提高模擬估計(jì)效果。3）在火災(zāi)穩(wěn)定燃燒階段的估計(jì)效果優(yōu)于燃燒初期，但可通過(guò)調(diào)整方程中的噪聲參數(shù)，使模擬估計(jì)效果達(dá)到較好水平。本文方法可用于一般隧道火災(zāi)場(chǎng)景，但由于實(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ù)。