劉夢(mèng)婷，朱 偉

(北京市科學(xué)技術(shù)研究院 北京城市系統(tǒng)工程研究中心 燃?xì)狻⒐峒暗叵鹿芫W(wǎng)運(yùn)行安全北京市重點(diǎn)實(shí)驗(yàn)室，北京100035)

0 引言

城市快速發(fā)展帶來(lái)了人口的迅速增加，地鐵站、超市、場(chǎng)館、寫(xiě)字樓等公共建筑內(nèi)聚集大量人員，工廠(chǎng)與居住區(qū)、人口密集區(qū)相鄰，且隨著我國(guó)國(guó)際化程度不斷提高，傳統(tǒng)和非傳統(tǒng)威脅并存，建筑物內(nèi)的人員疏散除了面臨火災(zāi)等安全問(wèn)題，還可能面臨生產(chǎn)安全事故、恐怖主義或犯罪等導(dǎo)致的毒氣事件，因此，緊急情況下人員快速、安全地撤離具有重要的現(xiàn)實(shí)意義。例如，1995年日本邪教組織人員在東京地下鐵三線(xiàn)共5列列車(chē)上投放沙林毒氣，造成13人遇難及5 510人以上受傷；2013年吉林寶源豐禽業(yè)公司發(fā)生特大火災(zāi)爆炸事故，造成121人遇難、76人受傷，受傷致死的原因主要為氨氣中毒和燒傷；2017年德國(guó)漢堡機(jī)場(chǎng)逾50人疑吸入有毒氣體導(dǎo)致身體不適。

現(xiàn)有的研究多側(cè)重于的大尺度毒氣擴(kuò)散導(dǎo)致的區(qū)域路網(wǎng)疏散。熊立春等[1]對(duì)液氨泄漏擴(kuò)散濃度進(jìn)行了高斯擴(kuò)散模擬，在此基礎(chǔ)上建立了風(fēng)險(xiǎn)矩陣和人口密度矩陣，對(duì)人員疏散進(jìn)行了最優(yōu)路徑分析；Shen等[2]提出了毒氣泄露事故的疏散范圍的計(jì)算方法，以毒負(fù)荷最小為目標(biāo)函數(shù)提出了路徑選擇優(yōu)化方法；賀政綱等[3]提出了與氣體濃度區(qū)間、離事故點(diǎn)的距離以及距離事故發(fā)生時(shí)間3個(gè)因素有關(guān)的恐慌因子，對(duì)有毒氣體泄漏事故的人員疏散范圍進(jìn)行修正，以疏散路徑總長(zhǎng)度最短為目標(biāo)函數(shù)，建立了疏散路徑優(yōu)化模型；Gai等[4]提出了毒氣泄漏事故中的風(fēng)險(xiǎn)評(píng)估模型，并以疏散時(shí)間和路徑風(fēng)險(xiǎn)最小為目標(biāo)函數(shù)，建立了多目標(biāo)疏散路徑優(yōu)化模型。

建筑內(nèi)的緊急疏散研究側(cè)重于火災(zāi)場(chǎng)景，對(duì)毒氣事故場(chǎng)景下的人員疏散研究較少。Karafyllidis等[5]建立了基于元胞自動(dòng)機(jī)的火災(zāi)擴(kuò)散模型，元胞的狀態(tài)為元胞過(guò)火面積與總面積的比值，根據(jù)相鄰元胞狀態(tài)轉(zhuǎn)移規(guī)則確定下一時(shí)刻元胞狀態(tài)；Georgoudas等[6]利用Karafyllidis & Thanailakis的火災(zāi)擴(kuò)散模型，模擬了行人接近火源的疏散過(guò)程；Zheng等[7]在元胞自動(dòng)機(jī)地面場(chǎng)的基礎(chǔ)上擴(kuò)展了火場(chǎng)以模擬火災(zāi)中行人的疏散過(guò)程；Jirasingha等[8]利用FDS+Evac進(jìn)行了模擬；Liu等[9]利用Kundu等[10]的方法構(gòu)建一個(gè)簡(jiǎn)化的火災(zāi)擴(kuò)散模型，疏散模型中行人根據(jù)火勢(shì)是否進(jìn)入視野區(qū)域以及障礙物的位置決定是否改變路徑；Fang等[11]考慮到了火災(zāi)產(chǎn)物如煙氣和一氧化碳對(duì)人員生理和行為的影響，建立了基于多格子的疏散模型，其中出口選擇取決于人員避開(kāi)能見(jiàn)度低和溫度高區(qū)域的期望值；Wan等[12]將社會(huì)力模型和高斯煙團(tuán)模式結(jié)合，用于模擬突發(fā)毒氣事件的人員疏散，其中毒氣對(duì)行人個(gè)體產(chǎn)生的排斥力由行人距離毒氣源的距離決定，且不同距離的毒負(fù)荷也決定了行人的運(yùn)動(dòng)速度?，F(xiàn)有研究中僅考慮了行人運(yùn)動(dòng)過(guò)程中對(duì)危險(xiǎn)的避讓行為，而出口選擇是行人運(yùn)動(dòng)前和過(guò)程中根據(jù)風(fēng)險(xiǎn)信息的決策和判斷，相關(guān)研究較少，且沒(méi)有考慮行人對(duì)危險(xiǎn)信息的獲取和傳遞。

本文分析了人群分布和毒氣分布對(duì)個(gè)體出口選擇的影響，建立了疏散出口選擇模型，并基于地面場(chǎng)疏散模型擴(kuò)展了毒氣場(chǎng)表達(dá)和毒氣信息傳遞模式。將出口選擇模型引入擴(kuò)展的疏散模型進(jìn)行人群疏散仿真，對(duì)比了不同模式下的疏散模擬結(jié)果，并分析了出口間隔、毒氣場(chǎng)參數(shù)、毒氣源位置和擴(kuò)散范圍等對(duì)疏散結(jié)果的影響。

1 疏散動(dòng)力學(xué)模型

1.1 出口選擇影響因素

空間壓力主要體現(xiàn)在有毒氣體擴(kuò)散對(duì)行人的危害。決策者主要通過(guò)2個(gè)方面感知危險(xiǎn)：一是通過(guò)視覺(jué)和嗅覺(jué)感知毒氣濃度；二是從其他行人處獲取毒氣信息。因此，一方面為了模擬行人對(duì)毒氣濃度的感知模式，在地面場(chǎng)的基礎(chǔ)上引入毒氣場(chǎng)，作用包括3方面：一是行人的避險(xiǎn)行為使得個(gè)體移動(dòng)的方向是由毒氣濃度大的位置向濃度小的位置；二是通過(guò)某個(gè)時(shí)間位置上毒氣濃度和人體毒負(fù)荷的對(duì)比，確定行人的生存狀態(tài)；三是個(gè)體進(jìn)入致傷區(qū)域感知到毒氣濃度，從而獲得毒氣信息。另一方面，為了模擬毒氣信息的傳遞模式，假設(shè)感知到毒氣濃度的行人會(huì)通過(guò)聲音和手勢(shì)等警示他人，因此，未進(jìn)入毒氣擴(kuò)散范圍內(nèi)的行人也可能提前獲取他人的危險(xiǎn)信息改變路徑選擇。

1.2 考慮擁堵和毒氣避免的出口選擇模型

出口選擇模型基于多元Logit離散選擇模型，效用函數(shù)根據(jù)疏散時(shí)間和毒氣源位置構(gòu)建。假設(shè)建筑內(nèi)有K個(gè)出口(k=1,2，…,K)，行人i(位置(ix,iy))在時(shí)刻t選擇出口k的概率由下式計(jì)算：

(1)

(2)

效用函數(shù)中的行走時(shí)間因素為：

(3)

效用函數(shù)中的等待時(shí)間因素為：

(4)

效用函數(shù)中的毒氣源因素為：

(5)

1.3 考慮毒氣危害擴(kuò)散的行人疏散模型

本文使用二維元胞自動(dòng)機(jī)(Cellular Automata，CA)模型來(lái)模擬逃生者的運(yùn)動(dòng)。將室內(nèi)空間劃分為L(zhǎng)×W的元胞網(wǎng)格，每個(gè)元胞對(duì)應(yīng)的空間是0.5 m×0.5 m，取值1(有人)或0(無(wú)人)2種狀態(tài)。每時(shí)步逃生者移動(dòng)到相鄰的未被占用的元胞中，或者保持原地不動(dòng)，采用Von Neumann型鄰居。

行人的方向選擇由靜態(tài)地面場(chǎng)(Static Floor Field，SFF)[14-15]決定?？紤]到行人對(duì)毒氣危害避免行為，加入了毒氣地面場(chǎng)(Concentration Floor Field，CFF)表示毒氣釋放后的濃度分布。因此，逃生者選擇方向位置(x,y)的轉(zhuǎn)移概率表示為：

P(x,y)∝exp(-kSS(x,y))exp(-kCC(x,y))

(6)

式中：參數(shù)kS和kC分別表示行人對(duì)出口路徑和毒氣分布的掌握程度。

考慮公共建筑空間內(nèi)化學(xué)恐怖襲擊的場(chǎng)景，使用簡(jiǎn)化的高斯煙團(tuán)模型模擬毒氣擴(kuò)散模式。假定毒氣為氯氣，泄漏總量為5 kg，則靜風(fēng)條件下，在30 min內(nèi)，點(diǎn)位(x,y)的毒氣最大濃度表示為[16]：

C(x,y)=26.15×exp(-2.263 5×(x2+y2))

(7)

根據(jù)英國(guó)衛(wèi)生安全執(zhí)行局(Health and Safety Executive，HSE)提出的毒負(fù)荷概念進(jìn)行處理，當(dāng)氯氣瞬間泄露時(shí)，氯氣濃度在2 500 mg以上為致死區(qū)；氯氣濃度在100 mg到2 500 mg之間為重傷區(qū)；氯氣濃度在40 mg到100 mg之間為輕傷區(qū)。

2 仿真實(shí)驗(yàn)與討論

2.1 仿真步驟與實(shí)驗(yàn)場(chǎng)景

個(gè)體決策的仿真流程如圖1所示。VCAk是出口k的可見(jiàn)范圍(Visibility Catchment Area，VCA)，這里假設(shè)為圓形，半徑等于行人的可見(jiàn)距離，m。CWT(Cumulative Waiting Time)為累積等待時(shí)間，s，每當(dāng)行人的忍耐時(shí)間耗盡后CWT重置為0。

圖1 個(gè)體的疏散出口選擇流程Fig.1 Flow chart of exit choice for one occupant

1)初始化行人位置，所有行人根據(jù)出口選擇模型確定初始的目標(biāo)出口。

2)根據(jù)行人疏散模型動(dòng)態(tài)更新行人位置。在此過(guò)程中，一旦行人第一次進(jìn)入某一出口的可見(jiàn)區(qū)域VCA、累積等待時(shí)間大于其忍耐時(shí)間CWT>Pa、行人獲取了毒氣源的信息即OGK由0變1，則根據(jù)出口選擇模型對(duì)目標(biāo)出口進(jìn)行一次重新選擇。

3)更新所有逃生者的位置，即更新元胞狀態(tài)。當(dāng)行人進(jìn)入毒氣濃度致死區(qū)后，行人位置不再變化，元胞狀態(tài)始終為1，等同于障礙物。到達(dá)出口后的行人將從系統(tǒng)中移除，當(dāng)所有行人都離開(kāi)疏散空間仿真結(jié)束，否則返回到步驟2。

2.2 疏散模型驗(yàn)證和參數(shù)設(shè)定

通過(guò)真人實(shí)驗(yàn)數(shù)據(jù)對(duì)疏散模型進(jìn)行驗(yàn)證和參數(shù)設(shè)定。實(shí)驗(yàn)構(gòu)造了1個(gè)具有3個(gè)出口的房間，面積為8.5 m×5.5 m，志愿者為50人，進(jìn)行了2種情景的多次實(shí)驗(yàn)[17]，實(shí)驗(yàn)與模擬情景如圖2與圖3所示。情景1：行人初始隨機(jī)分布，開(kāi)燈環(huán)境完全可見(jiàn)，房間開(kāi)雙出口如圖2(a)所示；情景2：行人初始集中分布在遠(yuǎn)離出口的位置，關(guān)燈環(huán)境不可見(jiàn)，房間開(kāi)單個(gè)出口如圖3(a)所示。利用元胞自動(dòng)機(jī)疏散模型分別對(duì)這2種情景的行人疏散進(jìn)行了計(jì)算機(jī)仿真，結(jié)果為100次仿真的平均值。

圖2 可見(jiàn)環(huán)境下的疏散實(shí)驗(yàn)和模擬過(guò)程快照Fig.2 Snapshot of the real and simulated evacuation process in a visible room

圖3 不可見(jiàn)環(huán)境下的疏散實(shí)驗(yàn)和模擬過(guò)程快照Fig.3 Snapshot of the real and simulated evacuation process in an invisible room

真人實(shí)驗(yàn)和模擬實(shí)驗(yàn)的結(jié)果擬合較好，說(shuō)明疏散模型的有效性。在完全可見(jiàn)環(huán)境下，模型的靜態(tài)場(chǎng)參數(shù)設(shè)定為kS=10，計(jì)算得到出口流率f=1.92人/(m·s-1)。不可見(jiàn)環(huán)境下，模型的靜態(tài)場(chǎng)參數(shù)設(shè)定為kS=5，計(jì)算得到出口流率f=1.75 人/(m·s-1)。因此，以下模擬實(shí)驗(yàn)采用kS=5；出口流率f=1.75 人/(m·s-1)；行人速度為1 m/s，每時(shí)步行人移動(dòng)1個(gè)元胞，因此，每時(shí)步耗時(shí)為0.5 m÷1 m/s = 0.5 s。2種不同的情景下，疏散人數(shù)與時(shí)間的關(guān)系對(duì)比如圖4和圖5所示。

圖4 可見(jiàn)環(huán)境下的疏散實(shí)驗(yàn)和模擬結(jié)果對(duì)比Fig.4 Comparison of the results of experiment and simulation in a visible room

圖5 不可見(jiàn)環(huán)境下的疏散實(shí)驗(yàn)和模擬結(jié)果對(duì)比Fig.5 Comparison of the results of experiment and simulation in an invisible room

2.3 毒氣泄漏場(chǎng)景下的仿真效果

以1個(gè)面積為25 m×25 m(50個(gè)×50個(gè)元胞網(wǎng)格)的雙出口大廳為疏散場(chǎng)景，2個(gè)出口位于同一面墻，出口寬度為1個(gè)元胞，即wk=0.5 m，VCA半徑為5 m。行人總數(shù)為500人，為了反映出口選擇的現(xiàn)象，行人的初始位置隨機(jī)分布于偏向1個(gè)出口的半個(gè)房間內(nèi)。圖6和圖7分別為無(wú)氯氣和有氯氣泄漏的疏散仿真過(guò)程，反映了所構(gòu)建出口模型的有效性。

2.4 出口間距對(duì)疏散結(jié)果的影響

觀察2個(gè)出口之間的間距對(duì)疏散結(jié)果的影響，分別設(shè)置出口間距為5，10，15，20和24 m，VCA半徑為5 m，得到2個(gè)出口的疏散人數(shù)和時(shí)間分別如圖8和圖9所示。當(dāng)出口間距為15 m時(shí)，總疏散時(shí)間最短，2個(gè)出口的疏散人數(shù)和疏散時(shí)間較平衡且不確定性較小；當(dāng)出口間距為20 m時(shí)，2個(gè)出口對(duì)應(yīng)的疏散人數(shù)的方差比較大，說(shuō)明行人在2個(gè)出口之間較頻繁地往返，不確定性大。因此，出口間距的設(shè)置存在臨界值，應(yīng)大于或小于該臨界值，否則容易造成行人頻繁改變目標(biāo)出口，影響到疏散的安全性。

圖6 無(wú)毒氣泄漏場(chǎng)景的疏散過(guò)程Fig.6 Evacuation process in the scenario without toxic gas leakage

圖7 毒氣泄漏場(chǎng)景的疏散過(guò)程Fig.7 Evacuation process in the scenario with toxic gas leakage

圖8 疏散人數(shù)與出口間隔的關(guān)系Fig.8 Relationship of number of the occupants evacuated and exit separation

圖9 疏散時(shí)間與出口間隔的關(guān)系Fig.9 Relationship of evacuation time and exit separation

該臨界值與人的視野半徑有關(guān)。圖10顯示了VCA半徑R分別為2.5，5和7.5 m時(shí)2個(gè)出口的通過(guò)人數(shù)。視野半徑為2.5 m且出口間距為10 m、視野半徑為5 m且出口間距為20 m、視野半徑為7.5 m且出口間距為24 m時(shí)，2個(gè)出口對(duì)應(yīng)的疏散人數(shù)的方差相比較最大。隨著視野范圍的增大，所獲得的不同出口路徑上的擁堵和毒氣等危險(xiǎn)信息越多，使得行人不會(huì)盲目改變選擇。

圖10 不同視野半徑下疏散人數(shù)與出口間隔的關(guān)系Fig.10 Relationship of number of the occupants evacuated and exit separation for different visual radius

2.5 毒氣場(chǎng)參數(shù)對(duì)疏散結(jié)果的影響

考慮不同的毒氣擴(kuò)散場(chǎng)景，如圖11和圖12所示。行人的初始位置隨機(jī)分布于整個(gè)房間，2個(gè)出口分別位于相對(duì)的2面墻上，2個(gè)場(chǎng)景的毒氣源位置相同，但場(chǎng)景二的毒氣擴(kuò)散范圍比場(chǎng)景一大。由于毒氣泄漏為穩(wěn)態(tài)分布，傷亡人數(shù)隨時(shí)間變化較小，為了便于計(jì)算和對(duì)比，此處不設(shè)定死亡模式。

圖11 場(chǎng)景一：毒氣擴(kuò)散范圍較小Fig.11 Scenario 1: the smaller diffusion range of toxic gas

圖12 場(chǎng)景二：毒氣擴(kuò)散范圍較大Fig.12 Scenario 2: the larger diffusion range of toxic gas

圖13 疏散人數(shù)與毒氣場(chǎng)參數(shù)kC的關(guān)系Fig.13 Relationship of number of the occupants evacuated and kC

圖14 疏散時(shí)間與毒氣場(chǎng)參數(shù)kC的關(guān)系Fig.14 Relationship of evacuation time and kC

疏散人數(shù)和時(shí)間與毒氣場(chǎng)參數(shù)kC的關(guān)系如圖13和圖14所示。同一場(chǎng)景下，改變毒氣場(chǎng)參數(shù)kC的取值，發(fā)現(xiàn)kC不同的情況下，2個(gè)出口的疏散人數(shù)和時(shí)間差別較??；不同場(chǎng)景下，即毒氣擴(kuò)散范圍不同的情況下，2個(gè)出口的疏散人數(shù)和時(shí)間差別也較小。因此，毒氣的擴(kuò)散范圍和行人對(duì)毒氣分布的掌握程度對(duì)疏散效果的影響較小。

但是從圖13(b)和圖14(b)發(fā)現(xiàn)，相比較場(chǎng)景二，場(chǎng)景一的疏散人數(shù)和時(shí)間的方差較大，說(shuō)明場(chǎng)景一中的行人排隊(duì)到達(dá)出口1(距離毒氣源較遠(yuǎn)的出口)的情況下，很可能冒險(xiǎn)選擇出口2(距離毒氣源較近的出口)。主要是因?yàn)閳?chǎng)景二的毒氣擴(kuò)散范圍較大，初始狀態(tài)下獲取毒氣信息的人數(shù)較多，通過(guò)信息傳遞規(guī)則，更多的人獲取了毒氣信息而趨向于選擇距離毒氣源較遠(yuǎn)的出口。

2.6 毒氣源位置對(duì)疏散結(jié)果的影響

仿真場(chǎng)景出口間距為15 m(見(jiàn)圖7)。由于出口的對(duì)稱(chēng)性，僅在1/4的房間內(nèi)改變毒氣源的位置，假設(shè)房間的左下角的坐標(biāo)為(0 m,0 m)，則出口的坐標(biāo)為(5 m,25 m)和(20 m,25 m)，分別設(shè)置毒氣源的位置如圖15所示，進(jìn)行疏散仿真，結(jié)果如圖16和圖17所示，說(shuō)明毒氣源位置對(duì)疏散結(jié)果具有顯著影響。

圖15 毒氣源的位置坐標(biāo)(單位：m)Fig.15 Coordinates of the cell position of toxic gas source(unit：m)

圖16 不同毒氣源位置下的2個(gè)出口的疏散人數(shù)Fig.16 Numbers of the occupants evacuated through two exits for different positions of toxic gas source

圖17 不同毒氣源位置下的2個(gè)出口的疏散時(shí)間Fig.17 Evacuation time of two exits for different positions of toxic gas source

從圖17可以看出，毒氣源越接近出口1和出口2的中界線(xiàn)上，2個(gè)出口的疏散時(shí)間越短且越接近相同；毒氣源越接近出口1，出口1和出口2的疏散時(shí)間都逐漸增加。主要因?yàn)樾腥嗽诒苊鈸矶潞捅苊舛練庵g權(quán)衡，極度擁堵的情況下有人冒險(xiǎn)回到出口1，由于行人有一定的忍耐時(shí)間，所以出口重新選擇往往發(fā)生在疏散最后階段，從而2個(gè)出口的疏散時(shí)間相差不大。

3 結(jié)論

1)出口間距的設(shè)置存在臨界值，應(yīng)大于或小于該臨界值，否則容易造成行人頻繁變換目標(biāo)出口，反而容易延誤疏散時(shí)機(jī)，且人的視野半徑越大，該臨界值越大。

2)毒氣源位置對(duì)疏散結(jié)果具有顯著影響。毒氣的擴(kuò)散范圍和行人對(duì)毒氣分布的掌握程度對(duì)疏散結(jié)果的影響較小，但是毒氣擴(kuò)散范圍由于影響到毒氣信息攜帶者的數(shù)量進(jìn)而影響到疏散結(jié)果的不確定性。

3)設(shè)定毒氣泄漏模式為一定時(shí)間內(nèi)的穩(wěn)態(tài)分布，后續(xù)將進(jìn)一步探究有風(fēng)條件下毒氣動(dòng)態(tài)擴(kuò)散對(duì)行人出口選擇的影響。