建筑物出口的人員疏散研究1

趙士達(dá) 張 楠 趙 穎

?

建筑物出口的人員疏散研究1

趙士達(dá) 張 楠 趙 穎

（天津市地震局，天津 300201）

根據(jù)建筑物出口疏散的特點(diǎn)，提出了一種改進(jìn)型的元胞自動(dòng)機(jī)模型。研究了影響人員疏散速度的兩個(gè)因素：人員之間的相互擁擠產(chǎn)生的沖突和人員移動(dòng)到出口時(shí)由于轉(zhuǎn)彎導(dǎo)致的移動(dòng)速度減慢。使用數(shù)學(xué)方法推導(dǎo)得出了摩擦沖突函數(shù)和轉(zhuǎn)彎函數(shù)。通過(guò)計(jì)算機(jī)仿真與實(shí)際疏散實(shí)驗(yàn)對(duì)比發(fā)現(xiàn)，在不考慮轉(zhuǎn)彎因素的情況下，仿真結(jié)果與實(shí)驗(yàn)結(jié)果偏差較大，而在同時(shí)考慮相互沖突和轉(zhuǎn)彎因素時(shí)的仿真結(jié)果與實(shí)驗(yàn)結(jié)果偏差很小，驗(yàn)證了這兩個(gè)因素是疏散仿真不可忽略的因素。通過(guò)分析仿真結(jié)果，驗(yàn)證了本文所采用的模型和函數(shù)能夠準(zhǔn)確地反映實(shí)際疏散過(guò)程，具有較高的實(shí)用價(jià)值。

元胞自動(dòng)機(jī) 人群疏散 摩擦函數(shù) 轉(zhuǎn)彎函數(shù) 疏散仿真

引言

隨著社會(huì)的發(fā)展，人們對(duì)地震、火災(zāi)等突發(fā)事件越來(lái)越重視。為了減少人員傷亡和財(cái)產(chǎn)損失，各級(jí)地方政府先后制定了各種應(yīng)對(duì)突發(fā)事件的應(yīng)急預(yù)案。其中，制定人員疏散方案是應(yīng)急預(yù)案的一項(xiàng)重點(diǎn)內(nèi)容（趙士達(dá)等，2014），需要對(duì)人的運(yùn)動(dòng)行為進(jìn)行系統(tǒng)研究。在過(guò)去的十幾年里，各國(guó)學(xué)者對(duì)人員疏散行為進(jìn)行了廣泛深入的研究。其中人員疏散中的決策過(guò)程主要通過(guò)元胞自動(dòng)機(jī)（Cellular Automata，CA）的規(guī)則進(jìn)行演化并表現(xiàn)出來(lái)（朱剛等，2006；楊兆升等，2011）。自1966年由Von Neumann首次提出元胞自動(dòng)機(jī)的概念以來(lái)，因?yàn)樵詣?dòng)機(jī)能將局部變化的反饋?zhàn)饔帽憩F(xiàn)出來(lái)，并在模擬系統(tǒng)中具有較好的突變和自組織等特性，已被廣泛應(yīng)用于人群疏散、城市交通控制等領(lǐng)域（岳昊，2008）。為了減輕災(zāi)害發(fā)生時(shí)的損失，國(guó)內(nèi)外諸多學(xué)者以CA模型為基礎(chǔ)，從不同的角度對(duì)人員疏散過(guò)程進(jìn)行了深入的研究，如：Heather等（2009）和Gong等（2000）將模糊規(guī)則引入到CA模型，建立了模糊元胞自動(dòng)機(jī)（Fuzzy Cellular Automata，F(xiàn)CA）模型，并廣泛應(yīng)用于人員流動(dòng)、城市發(fā)展等研究領(lǐng)域；朱藝等（2007）針對(duì)人口密度和出口條件對(duì)疏散過(guò)程進(jìn)行了研究；張俊娜等（2012）將人工勢(shì)能場(chǎng)與CA模型相結(jié)合，對(duì)人員疏散過(guò)程進(jìn)行了研究。經(jīng)過(guò)統(tǒng)計(jì)調(diào)查發(fā)現(xiàn)，在緊急情況下人員疏散時(shí)間主要花費(fèi)在建筑物出口附近（Steffen等，2009）。所以，本文將對(duì)建筑物出口人員疏散進(jìn)行研究?？紤]到出口附近等待被疏散的人員通常是后排人與前排人交錯(cuò)排列，為此還對(duì)經(jīng)典的二維元胞自動(dòng)機(jī)空間模型進(jìn)行了改進(jìn)，并在此模型基礎(chǔ)上分析了人之間的摩擦力和人運(yùn)動(dòng)到出口附近轉(zhuǎn)彎對(duì)他人移動(dòng)速度的影響。

1 改進(jìn)的出口模型

經(jīng)典的元胞自動(dòng)機(jī)二維網(wǎng)格劃分方法是將整個(gè)空間劃分為若干個(gè)大小相等、整齊排列的正方形，它的馮·諾依曼型鄰域如圖1（a）所示。但是在人員疏散時(shí)，后排人為了獲得更好的視野會(huì)選擇交叉排列的方式，所以本文采用交叉式的網(wǎng)格劃分方式，其鄰域如圖1（b）所示。從圖1（b）可以看出，人會(huì)根據(jù)運(yùn)動(dòng)規(guī)律移動(dòng)到其鄰域元胞中，但是在實(shí)際疏散中出口附近的人不會(huì)選擇后退，即不會(huì)向遠(yuǎn)離出口的方向運(yùn)動(dòng)。所以，本文對(duì)出口附近元胞自動(dòng)機(jī)模型進(jìn)行了改進(jìn)，改進(jìn)后的模型空間網(wǎng)格劃分方式如圖2所示。元胞E是緊鄰出口的元胞，人需要經(jīng)過(guò)元胞E到達(dá)出口。筆者將元胞1、2、3、4的人移動(dòng)到元胞E的概率定義為∈［0，1］，并將在一個(gè)時(shí)間步長(zhǎng)內(nèi)通過(guò)元胞E離開(kāi)出口的概率定義為∈［0，1］。當(dāng)元胞1、2、3、4只有1個(gè)人要進(jìn)入元胞E時(shí)，該人進(jìn)入元胞E的概率為1。當(dāng)有2個(gè)或者2個(gè)以上人要進(jìn)入元胞E時(shí)，就會(huì)產(chǎn)生沖突，即最多只有1個(gè)人可以進(jìn)入元胞E，他人將保持在原來(lái)的位置不動(dòng)。

圖1 整齊排列和交叉排列兩種模型的鄰域定義方式

圖2 改進(jìn)的出口模型

同時(shí)將全部人員保持在原來(lái)位置的概率定義為，那么人員成功移動(dòng)到元胞E的概率就是1-。人員保持在原來(lái)位置的概率φ(k)可以用公式（1）來(lái)表示，其中k為同時(shí)向元胞E移動(dòng)的人員數(shù)量。在下文中筆者將φ稱(chēng)為摩擦系數(shù)。

（1）

2 疏散人員之間的摩擦函數(shù)

通過(guò)公式（1）可以看出，在圖2模型中當(dāng)ke≥2時(shí)，摩擦系數(shù)是一個(gè)常量，它并不能反映出擁擠時(shí)人員之間的相互作用。所以本文引入摩擦函數(shù)ζ∈［0，1］。其中，ζ∈［0，1］表示在多個(gè)人員同時(shí)向同一個(gè)元胞移動(dòng)時(shí)，不給其他人員讓行的概率；摩擦函數(shù)ζ表示的物理意義是疏散時(shí)發(fā)生沖突并且沖突不能夠解決的概率。這里，首先分析沖突發(fā)生時(shí)沖突被解決的兩種情況：一是當(dāng)沖突發(fā)生時(shí)，所有人員都選擇保持在原來(lái)位置上給其他人讓路，不存在多人競(jìng)爭(zhēng)同一個(gè)位置的情況；二是只有一個(gè)人選擇不給其他人讓路，其他的人員選擇保持在原來(lái)的位置不動(dòng)，也不存在人員相互競(jìng)爭(zhēng)的情況。摩擦函數(shù)的數(shù)學(xué)表達(dá)式如公式（2）：

在公式（2）中，（1-ζ）k表示沖突發(fā)生時(shí)沖突被解決的第一種情況；kζ（1-ζ）k-1表示沖突被解決的第二種情況。由公式（2）可以得出：當(dāng)k=1時(shí)，ζ=0；當(dāng)k=∞時(shí)，ζ=1。其中，k=1時(shí)的物理意義是只有一個(gè)人試圖移動(dòng)到出口元胞時(shí)，不會(huì)發(fā)生沖突，移動(dòng)到出口元胞的概率為1；k=∞時(shí)的物理意義是非常多的人試圖同時(shí)移動(dòng)到出口元胞時(shí)，沖突十分嚴(yán)重，幾乎沒(méi)有人能移動(dòng)到出口元胞。從Daichi等（2007）的實(shí)驗(yàn)中可以看出，當(dāng)k≥2時(shí)摩擦系數(shù)始終為常量，而摩擦函數(shù)隨著k的增大而增大。這說(shuō)明摩擦函數(shù)反映了不同k的阻塞強(qiáng)度。這里，筆者引用Daichi等（2007）的摩擦系數(shù)和摩擦函數(shù)隨k變化的函數(shù)圖，如圖3所示。

圖3 摩擦系數(shù)和摩擦函數(shù)隨k改變的變化

Fig. 3 Values of the frictional parameter and the frictional function vs.e

從圖3可以看出，摩擦函數(shù)隨著k的增大而增大。同時(shí)從公式（2）也可以看出，摩擦函數(shù)的大小與產(chǎn)生沖突的人員數(shù)量有關(guān)，與元胞自動(dòng)機(jī)網(wǎng)格劃分方式無(wú)關(guān)。也就是說(shuō)，本文所采用的摩擦函數(shù)不僅適用于本文所采用的元胞空間劃分方式，而且也可以用于其它的元胞空間劃分方式。

3 人員運(yùn)動(dòng)時(shí)的轉(zhuǎn)彎函數(shù)

在疏散過(guò)程中人員的運(yùn)動(dòng)方向會(huì)隨時(shí)發(fā)生變化。在每次運(yùn)動(dòng)方向發(fā)生變化時(shí)，人員會(huì)做出轉(zhuǎn)彎的動(dòng)作，這導(dǎo)致了人員的移動(dòng)速度減慢。本文引入轉(zhuǎn)彎函數(shù)用以表示行人在運(yùn)動(dòng)中由于轉(zhuǎn)彎導(dǎo)致的減速效果，其數(shù)學(xué)表達(dá)式如公式（3）：

（3）

其中θ∈［0，π］為人員移動(dòng)轉(zhuǎn)彎的角度。

例如在圖2中，元胞1和4的人員先要移動(dòng)到元胞E，再?gòu)脑鸈移動(dòng)到出口，這需要轉(zhuǎn)彎90°，即θ=90°；元胞2和3的人員要移動(dòng)到出口需要轉(zhuǎn)彎30°。Η≥0為慣性系數(shù)，它反應(yīng)人員慣性的大小。在本文第一部分提到了人員移動(dòng)到元胞E的概率和在一個(gè)時(shí)間步長(zhǎng)移動(dòng)到出口的概率，這兩個(gè)概率的關(guān)系可以通過(guò)一個(gè)簡(jiǎn)化的指數(shù)函數(shù)來(lái)表示（Ansgar等，2003），其數(shù)學(xué)表達(dá)式如公式（4）：

（4）

由公式（4）可以看出，在一定時(shí)轉(zhuǎn)彎的角度越大，人員移動(dòng)速度降低的越快，離開(kāi)出口的概率也就越小。

4 平均疏散速度

人員要移動(dòng)到出口需要先移動(dòng)到出口的鄰域元胞，假設(shè)出口的鄰域元胞個(gè)數(shù)為n，鄰域元胞中同時(shí)想移動(dòng)到出口元胞的人員個(gè)數(shù)為k。那么n個(gè)鄰域中，k個(gè)人員試圖移動(dòng)到出口元胞的概率（k）如公式（5）：

一個(gè)人成功移動(dòng)到出口的概率就是k個(gè)人同時(shí)試圖移動(dòng)到出口并且沖突被解決的所有可能，因此，一個(gè)人成功移動(dòng)到出口的概率(n)可以用公式（6）表示：

定義在時(shí)刻出口元胞沒(méi)有被疏散人員占據(jù)的概率為π（0）和出口元胞被其鄰域元胞中的人員占據(jù)的概率為π（），其中為出口元胞的第個(gè)鄰域，那么在+1時(shí)刻，出口元胞的狀態(tài)可以通過(guò)以下方程式表示：

（7）

在→∞時(shí)，通過(guò)方程式（7）可以推導(dǎo)出π∞()的數(shù)學(xué)表達(dá)式（8）：

其中∈［1，n］。因?yàn)椋杂晒剑?）和（8）可以得出π∞（0）的數(shù)學(xué)表達(dá)式（9）：

假設(shè)出口寬度為一個(gè)元胞的寬度，那么一個(gè)步長(zhǎng)時(shí)間內(nèi)的平均疏散人數(shù)可以表示為：

（10）

5 計(jì)算機(jī)仿真與疏散實(shí)驗(yàn)

本文對(duì)人員疏散進(jìn)行計(jì)算機(jī)仿真并與實(shí)際疏散仿真實(shí)驗(yàn)進(jìn)行對(duì)比。在仿真中需要對(duì)18名人員按9種方式從建筑物中疏散。其中建筑物為長(zhǎng)4m，寬3.5m的長(zhǎng)方形，出口位于長(zhǎng)邊中央處，寬度為0.5m，9種疏散方式如圖4所示。

在（H）和（I）兩種疏散方式中人員可以自由移動(dòng)，在其它疏散方式中人員按照?qǐng)D4中所示的方式移動(dòng)。平均疏散人數(shù)可以通過(guò)以下公式計(jì)算得出：

（11）

式中，和分別表示第個(gè)被疏散的人員和第個(gè)被疏散的人員；t和t分別表示第個(gè)人員和第個(gè)人員被疏散的時(shí)刻。

在（A）到（D）的疏散方式中，仿真初期人員運(yùn)動(dòng)相對(duì)較快，出口附近沒(méi)有沖突，計(jì)算平均疏散人數(shù)時(shí)要等人員運(yùn)動(dòng)狀態(tài)穩(wěn)定后進(jìn)行計(jì)算，所以本文取=5，=18。在（E）到（I）的疏散方式中，疏散末期出口附近人員數(shù)量減少，人員移動(dòng)速度相對(duì)較快，在計(jì)算平均疏散人數(shù)時(shí)不應(yīng)包括疏散末期，所以本文取=1，=14。

對(duì)這9種疏散方式分別進(jìn)行2次仿真。在第一次仿真時(shí)，只考慮人員之間的摩擦力作用，不考慮轉(zhuǎn)彎對(duì)疏散的影響，可計(jì)算得到平均疏散人數(shù)；在第二次仿真時(shí)，同時(shí)考慮摩擦力和轉(zhuǎn)彎對(duì)疏散的影響，可計(jì)算得到平均疏散人數(shù)。在仿真中將元胞的邊長(zhǎng)設(shè)定為0.5m，仿真步長(zhǎng)為0.3s。Daichi等（2009）使用最小二乘法，得到的=0.25、=0.26、=0.09，筆者直接使用這些參數(shù)進(jìn)行仿真，2次仿真的結(jié)果與實(shí)際試驗(yàn)結(jié)果的對(duì)比如圖5所示。

通過(guò)計(jì)算機(jī)仿真結(jié)果和實(shí)際試驗(yàn)結(jié)果綜合分型，筆者得出了以下幾點(diǎn)結(jié)論：

（1）對(duì)比（B）、（E）和（G）三種疏散方式可以發(fā)現(xiàn)，隨著n（n≥2）的增大，出口附近發(fā)生沖突也越嚴(yán)重，疏散速度也越慢。

（2）當(dāng)n一定時(shí)，在不考慮轉(zhuǎn)彎函數(shù)的仿真時(shí)，（B）、（C）和（D）三種疏散方式的仿真疏散速度基本相同，但其仿真速度與實(shí)際試驗(yàn)的速度相差較大；在考慮轉(zhuǎn)彎函數(shù)的仿真和實(shí)際試驗(yàn)時(shí)，（B）、（C）和（D）三種疏散方式的仿真疏散速度隨著θ增大而降低，與實(shí)際試驗(yàn)的疏散速度相差不大，驗(yàn)證了轉(zhuǎn)彎過(guò)程對(duì)人員疏散的減速效果。

（3）觀察（G）、（H）和（I）三種疏散方式發(fā)現(xiàn)，其疏散速度非常相似。這是因?yàn)椋℉）和（I）兩種疏散方式在人員之間會(huì)自由運(yùn)動(dòng)，當(dāng)達(dá)到穩(wěn)定狀態(tài)時(shí)人員之間會(huì)分別占據(jù)出口的4個(gè)鄰域，其疏散方式與（G）類(lèi)似。

（4）通過(guò)圖5可以發(fā)現(xiàn)，相比更接近實(shí)驗(yàn)結(jié)果，這說(shuō)明轉(zhuǎn)彎函數(shù)對(duì)疏散速度的影響十分重要，是疏散仿真不可忽略的因素之一。

（5）仿真中的結(jié)果與實(shí)驗(yàn)結(jié)果十分接近，但是在（B）方式中卻存在一定差異，這主要是由于在實(shí)際疏散中人員會(huì)有意識(shí)地按照交替順序進(jìn)行移動(dòng)，使得沖突得到一定程度的抵消從而提高了疏散效率，這種現(xiàn)象被稱(chēng)為“拉鏈”效應(yīng)（Yushi等，2012）。

6 結(jié)論與展望

根據(jù)建筑物出口人員疏散的特點(diǎn)，本文提出了一種改進(jìn)型的元胞自動(dòng)機(jī)空間劃分模型。介紹了影響人員疏散速度的兩個(gè)函數(shù)：摩擦函數(shù)和轉(zhuǎn)彎函數(shù)。根據(jù)數(shù)學(xué)理論推導(dǎo)得出了相關(guān)函數(shù)的數(shù)學(xué)表達(dá)式。同時(shí)，使用文中的數(shù)學(xué)公式和相關(guān)參數(shù)對(duì)人員疏散進(jìn)行了計(jì)算機(jī)仿真，分別得出了考慮轉(zhuǎn)彎因素和不考慮轉(zhuǎn)彎因素的人員疏散時(shí)間，并將這一仿真結(jié)果與Daichi等（2009）的疏散實(shí)驗(yàn)結(jié)果進(jìn)行了對(duì)比。通過(guò)對(duì)比發(fā)現(xiàn)，考慮轉(zhuǎn)彎函數(shù)的仿真結(jié)果與實(shí)際疏散結(jié)果更接近，更能充分地反映人員的疏散過(guò)程。

楊兆升，高學(xué)英，孫迪，2011. 城市交通疏散救援的元胞自動(dòng)機(jī)模型. 交通運(yùn)輸工程學(xué)報(bào), 2（4）：114—120.

岳昊，2008. 基于元胞自動(dòng)機(jī)的行人流仿真模型研究. 北京：北京交通大學(xué), 14—25.

張俊娜，范海菊，2012. 基于人工勢(shì)場(chǎng)和Agent的人員疏散系統(tǒng)仿真. 計(jì)算機(jī)應(yīng)用，32（6）：1753—1756.

趙士達(dá)，張楠，楊爽，2014. 基于Android系統(tǒng)的避難場(chǎng)所查詢(xún)軟件開(kāi)發(fā). 電子技術(shù)應(yīng)用，40（3）：133—136.

朱剛，馬良，2006. 基于元胞自動(dòng)機(jī)的物流系統(tǒng)選址模型. 上海理工大學(xué)學(xué)報(bào)，（1）：19—22.

朱藝，楊立中，李健，2007. 不同房間結(jié)構(gòu)下人員疏散的CA模擬研究. 火災(zāi)科學(xué)，16（3）：175—179.

Steffen B., Seyfried A., 2009. Methods for measuring pedestrian density, flow, speed and direction with minimal scatter. Physica A: Statistical Mechanics and its Applications, (12): 015.

Kirchen A., Nishinani K., Schadschneider A., 2003. Friction effects and clogging in a cellular automaton model for pedestrian dynamics. Physical Review E-PHYS REVE, 67: 056122.

Yanagisawa D., Nishinar K., 2007a. Mean-field theory for pedestrian outflow through an exit. Physical Review E-PHYS REVE, 76 (6): 75-80.

Yanagisawa D., Kimuna A., 2009b. Analysis on Pedestrian Outflow through an Exit with an Obstacle. ICROS-SICE International Joint Conference, 8: 5040—5045.

Gong Q.L., Mendel J.M., 2000. Interval type-2 fuzzy logic systems: theory and design. IEEE Transactions on Fuzzy Systems, 8 (5): 535—550.

Heather Betel, Paola Flocchini, 2009. On the asymptoti behavior of fuzzy cellular Automata. Electronic Notes in Theoretical Computer Science, 252: 23—40.

Yushi Suma, Daichi Yanagisawa., Katsuhiro Nishinari, 2012. Anticipation effect in pedestrian dynamics: Modeling and experiments. Physica A: Statistical Mechanics and its Applications, 391 (1): 248—263.

Research on the Safe Evacuation near the Building Exit

Zhao Shida, Zhang Nan and Zhao Ying

(Tianjin Seismological Bureau, Tianjin 300201, China)

The paper presents an improved cellular automaton model according to the feature ofevacuation near the outlet. We studied friction and turning factors that affect pedestrian evacuation speed, and derived expression offriction function and turning function by using mathematical methods. The average pedestrian outflow of the simulation that includes the effects of both the frictional function and the turning function agrees well with experiment result. In contrast, the simulation results that only include the effect of the frictional function are not corresponding to the experimental ones well. Simulation results show that friction and turning should not be ignored. By analyzing the simulation results, we verified that the model can accurately reflect the actual evacuation process and has practical value.

Cellular automaton; Pedestrian evacuation; Friction function; Turning function; Evacuation simulation

天津市地震局青年基金項(xiàng)目（121010）

2014-07-30

趙士達(dá)，男，生于1983年。工程師。主要研究方向：地震應(yīng)急指揮、地震現(xiàn)場(chǎng)通訊和災(zāi)害評(píng)估。E-mail：zhaoshida@163.com