杜晨曦 魯長(zhǎng)江 毛珊珊 程 奕

?

中小學(xué)校疏散演練人員疏散模型研究1

杜晨曦1）魯長(zhǎng)江1）毛珊珊2）程 奕1）

1）四川省地震局，成都 610200 2）成都理工大學(xué)，成都 610059

本文以元胞自動(dòng)機(jī)模型理論為指導(dǎo)，以人員疏散時(shí)間最短為目標(biāo)，構(gòu)建中小學(xué)校疏散演練人員疏散模型。本文從開(kāi)展地震應(yīng)急疏散模型研究出發(fā)，為地震應(yīng)急演練疏散訓(xùn)練軟件提供核心算法。通過(guò)對(duì)空間基礎(chǔ)構(gòu)建元素及搭建方法的分析研究，設(shè)計(jì)完成疏散虛擬場(chǎng)景的建模方法。最后開(kāi)展人員疏散算法構(gòu)建方式的研究，并結(jié)合應(yīng)急疏散演練人員行為規(guī)范的研究成果，完成疏散模型的總體設(shè)計(jì)，達(dá)到為應(yīng)急演練疏散軟件提供人員疏散演化模型以及為人員疏散過(guò)程合理性評(píng)判提供依據(jù)的目的。

地震應(yīng)急疏散模型 元胞自動(dòng)機(jī) Dijkstra算法 遞推 疏散策略

引言

地震是群災(zāi)之首，全球地震災(zāi)害多發(fā)，給人類(lèi)造成了巨大損失。我國(guó)屬于地震災(zāi)害多發(fā)國(guó)家，頻發(fā)的地震災(zāi)害對(duì)中國(guó)經(jīng)濟(jì)建設(shè)、社會(huì)穩(wěn)定和城市化進(jìn)程等構(gòu)成了嚴(yán)重威脅。特別是2008年汶川地震和2010年玉樹(shù)地震，造成了嚴(yán)重的人員傷亡和財(cái)產(chǎn)損失。而在汶川地震中四川省安縣桑棗中學(xué)2000多名師生在1分36秒全體成功逃生又為我們提供了寶貴的經(jīng)驗(yàn)，即科學(xué)、合理的日常疏散演練能大大提高地震逃生概率，為生命安全提供有力保障。

因此，開(kāi)展地震應(yīng)急疏散研究，特別是針對(duì)中小學(xué)的地震應(yīng)急疏散研究顯得尤為重要。全面分析和研究地震應(yīng)急疏散模型可以為減少地震災(zāi)害造成的人員和財(cái)產(chǎn)損失打下良好基礎(chǔ)。

目前，國(guó)內(nèi)外在人員疏散模型的研究上主要采用連續(xù)和離散兩類(lèi)模擬方法。其中連續(xù)模型多是把行人運(yùn)動(dòng)看作氣體、水流、顆粒流等物理介質(zhì)，借鑒一些物理方程來(lái)建立模型。如Pauls（1978）通過(guò)對(duì)火災(zāi)情況下人員疏散演習(xí)實(shí)測(cè)得出的多層建筑人員疏散時(shí)間經(jīng)驗(yàn)公式。而離散模型則把時(shí)間和空間劃分為不連續(xù)的“時(shí)間步”和“網(wǎng)格”，通過(guò)一定局部或整體的行為運(yùn)動(dòng)規(guī)律來(lái)模擬人員疏散過(guò)程，代表模型為元胞自動(dòng)機(jī)。國(guó)內(nèi)對(duì)人員疏散的研究起步較晚，但經(jīng)過(guò)近幾年的努力，已經(jīng)取得了不少成績(jī)。張培紅（2002）在東北大學(xué)校園等處對(duì)人員流動(dòng)狀態(tài)進(jìn)行了實(shí)際觀測(cè)，從中歸納出了彎道90°以及下樓梯時(shí)人員流動(dòng)速度衰減系數(shù)，建立了離散狀態(tài)人員疏散行動(dòng)的計(jì)算機(jī)模擬模型SHEBR-SGEM；張樹(shù)平（2004）調(diào)查了一百多起火災(zāi)，對(duì)火災(zāi)發(fā)生時(shí)人員所處的環(huán)境狀態(tài)、人員疏散行為等進(jìn)行了理論研究；中國(guó)科學(xué)技術(shù)大學(xué)的楊立中教授等利用元胞自動(dòng)機(jī)模型對(duì)疏散過(guò)程中的從眾行為、親情行為、建筑物出口動(dòng)力學(xué)特征、信息傳播等進(jìn)行了深入的探討（Yang等，2002）。

本文所提及的人員疏散模型是地震應(yīng)急疏散演練訓(xùn)練器軟件的核心組成部件之一。地震應(yīng)急疏散演練訓(xùn)練器軟件的應(yīng)用范圍是中小學(xué)校，軟件的研制目的是幫助中小學(xué)校通過(guò)計(jì)算機(jī)模擬的方式快速得出適合于本校實(shí)際情況的疏散演練方案，學(xué)校基于疏散演練方案開(kāi)展震前疏散演練訓(xùn)練，達(dá)到有序、高效、安全疏散的目的，做好震前疏散準(zhǔn)備工作，一旦發(fā)生地震，學(xué)校師生可以按照優(yōu)化后的疏散演練方案有序、快速地撤離。

1 模型的建立

本文基于城鎮(zhèn)中小學(xué)疏散場(chǎng)景模型，通過(guò)對(duì)迪科斯徹算法、元胞自動(dòng)機(jī)疏散模型以及人員調(diào)度策略的研究分析，結(jié)合人員分布情況及建筑物結(jié)構(gòu)特點(diǎn)，得出適合中小學(xué)校人員的疏散模型。

1.1 疏散場(chǎng)景模型的建立

把已有的疏散軟件及常用建模軟件建筑物內(nèi)部空間的搭建方式，與四川省中小學(xué)校典型建筑結(jié)構(gòu)類(lèi)型及內(nèi)部空間特點(diǎn)的走訪調(diào)研數(shù)據(jù)相結(jié)合，得出合理可行的疏散場(chǎng)景模型搭建方法。

1.1.1 疏散場(chǎng)景模型的分類(lèi)

一個(gè)疏散場(chǎng)景可用2種方法表示，即精細(xì)網(wǎng)絡(luò)模型和粗糙網(wǎng)絡(luò)模型。精細(xì)網(wǎng)絡(luò)模型使用大量的網(wǎng)絡(luò)節(jié)點(diǎn)表示一個(gè)建筑物空間，每個(gè)節(jié)點(diǎn)對(duì)應(yīng)一小塊面積的疏散空間，并且允許在疏散模型中對(duì)其進(jìn)行詳細(xì)的描述。這種模型通常是對(duì)建筑物幾何的精確描述，并且在建模時(shí)通常需要建筑物的圖紙。粗糙網(wǎng)絡(luò)模型使用節(jié)點(diǎn)和與之相連接的邊來(lái)表示一個(gè)疏散場(chǎng)景，一條邊可以表示疏散場(chǎng)景內(nèi)一個(gè)真實(shí)的連接關(guān)系，節(jié)點(diǎn)代表建筑物空間中一個(gè)獨(dú)立的疏散空間，比如一個(gè)教室、一個(gè)辦公室等。

1.1.2 疏散場(chǎng)景模型的搭建

通過(guò)對(duì)四川省中小學(xué)校的走訪調(diào)研選定建立模型的基礎(chǔ)環(huán)境，根據(jù)人員所在建筑物的位置和影響人員快速撤離的空間結(jié)構(gòu)等因素，梳理出疏散場(chǎng)景搭建所需的點(diǎn)、線、面以及障礙物、樓梯等構(gòu)建元素類(lèi)型，分析得出精細(xì)網(wǎng)絡(luò)模型是更為合理的疏散場(chǎng)景模型搭建方法。模型搭建步驟如下：

（1）通過(guò)自主研發(fā)編輯器編輯空間模型（未使用AutoCAD、Visio等軟件）。

（2）空間模型搭建由整體到局部、由框架到細(xì)節(jié)。

圖1 疏散場(chǎng)景模型的搭建流程

通過(guò)精細(xì)網(wǎng)絡(luò)模型搭建的教學(xué)樓空間結(jié)構(gòu)圖如圖2所示。

圖2 教學(xué)樓空間結(jié)構(gòu)示意圖

1.2 人員疏散模型的建立

建筑物內(nèi)的疏散人員可以從宏觀和微觀兩個(gè)角度表示。從宏觀角度，疏散人員可看作是一個(gè)群體，在這個(gè)群體中每個(gè)疏散人員是完全相同的，按照同樣的行為進(jìn)行疏散，沒(méi)有個(gè)體的特質(zhì)。從微觀角度，每個(gè)疏散人員都會(huì)有自己的個(gè)體特征，或者被隨機(jī)分配不同的屬性，在疏散過(guò)程中，每個(gè)疏散人員在不同的條件下可能表現(xiàn)出不同的疏散行為。近年來(lái)，在建筑物人員疏散的微觀模型中，元胞自動(dòng)機(jī)模型以其簡(jiǎn)單性、靈活性和高效性在疏散流仿真方面受到越來(lái)越多的關(guān)注。

我們結(jié)合Dijkstra算法、人員調(diào)度策略以及中小學(xué)教學(xué)樓這一特定的疏散環(huán)境對(duì)疏散流進(jìn)行建模仿真。

1.2.1 元胞自動(dòng)機(jī)模型原理

元胞自動(dòng)機(jī)模型（Kirchner，2002；Burstedde，2001；岳昊等，2010）是在一個(gè)由離散、有限狀態(tài)的元胞組成的元胞空間中，按照一定的局部規(guī)則，在離散的時(shí)間維上演化的動(dòng)力學(xué)系統(tǒng)。元胞自動(dòng)機(jī)系統(tǒng)狀態(tài)用一種常規(guī)的格子組成單元。每個(gè)單元可以處于幾種狀態(tài)中的一種（一般有兩種，0和1），元胞的狀態(tài)為0時(shí)表示空出，為1時(shí)表示有人或障礙物，采用二維Von-Neumann鄰居或Moore鄰居。

元胞自動(dòng)機(jī)模型通常采用并行更新規(guī)則，在同一時(shí)步，系統(tǒng)中的所有行人同步并行更新運(yùn)動(dòng)到相對(duì)應(yīng)的格點(diǎn)上。若出現(xiàn)多位行人同時(shí)選擇一個(gè)格點(diǎn)，則發(fā)生沖突，解決的辦法通常是采用等概率隨機(jī)地選擇其中一位到達(dá)該格點(diǎn)，其余保持在原來(lái)位置，等到下一個(gè)時(shí)間步再做出選擇。

1.2.2 Dijkstra算法原理

Dijkstra算法（Dijkstra's algorithm）由荷蘭計(jì)算機(jī)科學(xué)家Dijkstra（1959）提出。該算法使用了廣度優(yōu)先搜索解決非負(fù)權(quán)有向圖的單源最短路徑問(wèn)題，算法最終得到一個(gè)最短路徑樹(shù)。該算法常用于路由算法或者作為其他圖算法的一個(gè)子模塊。

圖3 Dijkstra算法示意圖

第三步 D尚未著色，繼續(xù)執(zhí)行第二步。

第三步 D尚未著色，繼續(xù)執(zhí)行第二步。

第三步 D尚未著色，繼續(xù)執(zhí)行第二步。

第三步 D尚未著色，繼續(xù)執(zhí)行第二步。

第三步 D尚未著色，繼續(xù)執(zhí)行第二步。

第三步 D尚未著色，繼續(xù)執(zhí)行第二步。

得到從A到D的最短距離10。

1.2.3 人員疏散的前提條件

人員在疏散時(shí)遵循以下前提條件：

（1）網(wǎng)格尺寸為：0.4m×0.4m。（相關(guān)文獻(xiàn)資料顯示我國(guó)人體平均體厚與肩寬投影面積為0.146m2，且元胞自動(dòng)機(jī)網(wǎng)格化通用網(wǎng)格大小劃分為0.4m×0.4m）。

（2）每個(gè)元胞有8個(gè)運(yùn)動(dòng)方向，每一個(gè)時(shí)間步內(nèi)人員只能移動(dòng)一個(gè)網(wǎng)格。

（4）有疏散人員的節(jié)點(diǎn)稱(chēng)為疏散源。每次遞推時(shí)每個(gè)疏散人員可能位于不同的疏散源，所以到每個(gè)出口的最短距離不同。

1.2.4 教室內(nèi)疏散人員模型設(shè)計(jì)

教室內(nèi)疏散人員模型規(guī)則：

（1）利用迪科斯切算法求出出口到每個(gè)格子的距離。

（2）利用元胞自動(dòng)機(jī)思想，在格點(diǎn)上的元胞選取離出口距離最小的格點(diǎn)前進(jìn)。

（3）如果離出口距離最小的格點(diǎn)位置被占據(jù)，元胞則自動(dòng)選擇次好位置，但絕不能走差位置（即距離出口位置比當(dāng)前位置還遠(yuǎn)的格點(diǎn)）。

（4）移動(dòng)時(shí)，先讓在出口位置的元胞移動(dòng)，離開(kāi)教室，再以出口為起點(diǎn)，由近及遠(yuǎn)調(diào)度其他元胞移動(dòng)。

（5）在每個(gè)元胞等待移動(dòng)時(shí)，均需計(jì)算等待時(shí)間，調(diào)度時(shí)優(yōu)先調(diào)度等待時(shí)間長(zhǎng)的元胞，讓該元胞選擇移動(dòng)格點(diǎn)。

（6）在元胞移動(dòng)時(shí)，為每個(gè)元胞增加一個(gè)行動(dòng)力，每次加1，總的行動(dòng)力不超過(guò)2，每走一步數(shù)據(jù)更新一次，以此解決速度問(wèn)題。

（7）把人員所在位置到出口的距離與教室內(nèi)總?cè)藬?shù)相結(jié)合，對(duì)教室內(nèi)人員選擇不同出口進(jìn)行總體調(diào)度。

需要注意的是，教室內(nèi)疏散模型設(shè)計(jì)的第（7）步是為了對(duì)前面6步進(jìn)行優(yōu)化而提出的調(diào)度策略。在模型仿真實(shí)驗(yàn)中，我們把沒(méi)有第（7）步（即未優(yōu)化版本）和有第（7）步（即優(yōu)化版本）的結(jié)果進(jìn)行對(duì)比，進(jìn)一步證明調(diào)度策略的作用。

1.2.5 整棟教學(xué)樓疏散人員模型設(shè)計(jì)

（3）將疏散時(shí)間數(shù)組，疏散人員分配矩陣，中間結(jié)果數(shù)組初始化成0：

（4）計(jì)算每次遞推時(shí)每個(gè)出口的疏散時(shí)間的中間值：

（6）更新每個(gè)出口對(duì)應(yīng)的疏散時(shí)間，完成第次遞推：

圖4 疏散模型程序流程圖

2 仿真

本文僅針對(duì)一個(gè)教室內(nèi)的人員進(jìn)行疏散模型仿真實(shí)驗(yàn)。仿真分為未優(yōu)化版本和優(yōu)化版本，每個(gè)版本又分別考慮了人員均勻和不均勻分布情況（表1）。每次仿真實(shí)驗(yàn)，都進(jìn)行了相關(guān)數(shù)據(jù)采集，下面將對(duì)各仿真結(jié)果進(jìn)行說(shuō)明。

表1 不同版本疏散完成數(shù)據(jù)的對(duì)比

（1）未優(yōu)化版本的人員均勻分布的疏散

如圖5所示，疏散開(kāi)始前，教室內(nèi)人員均勻分布。疏散開(kāi)始后，教室內(nèi)的人員開(kāi)始分別往左右出口移動(dòng)，到疏散即將完成時(shí)，左右出口的人數(shù)相差很小，疏散基本同時(shí)完成。

圖5 未優(yōu)化版本的人員均勻分布情況

（2）未優(yōu)化版本的人員不均勻分布的疏散

如圖6所示，疏散開(kāi)始前，教室內(nèi)人員是不均勻分布的。疏散開(kāi)始后，教室內(nèi)的人員開(kāi)始分別往左右出口移動(dòng)，到疏散即將完成時(shí)，左右出口的人數(shù)相差很大，很明顯，在這種情況下，右邊出口擁堵情況變得更加嚴(yán)重。疏散步數(shù)相對(duì)人員均勻分布的情況有所增加。

圖6 未優(yōu)化版本的人員不均勻分布情況

（3）優(yōu)化版本的人員均勻分布的疏散

在優(yōu)化版本里加入教室內(nèi)人員調(diào)度模型，把疏散人員數(shù)量與出口的距離結(jié)合后進(jìn)行算法優(yōu)化，優(yōu)化后的程序有一個(gè)控制臺(tái)窗口，其中記錄了給每個(gè)出口分配的人數(shù)以及總步數(shù)。如圖7所示，疏散開(kāi)始前，教室內(nèi)人員是均勻分布的。疏散開(kāi)始后，教室內(nèi)的人員開(kāi)始分別往左右出口移動(dòng)，到疏散即將完成時(shí)，左右出口的人數(shù)相差很小，疏散基本同時(shí)完成。

圖7 優(yōu)化版本的人員均勻分布情況

（4）優(yōu)化版本的人員不均勻分布的疏散

優(yōu)化后的不均勻分布狀態(tài)，如圖8所示。疏散開(kāi)始前，教室內(nèi)人員是不均勻分布的。疏散開(kāi)始后，教室內(nèi)的人員開(kāi)始分別往左右出口移動(dòng)，到疏散即將完成時(shí)，很明顯，左右兩個(gè)出口的分布更加平均，相對(duì)未優(yōu)化的版本疏散完成步數(shù)減少6步，說(shuō)明經(jīng)優(yōu)化后疏散時(shí)間減少。

圖8 優(yōu)化版本的人員不均勻分布情況

3 結(jié)論

根據(jù)本文的研究目標(biāo)有針對(duì)性地對(duì)四川省部分中小學(xué)校進(jìn)行了調(diào)研，在收集基礎(chǔ)資料的同時(shí)也了解了學(xué)校疏散工作開(kāi)展情況。

國(guó)外現(xiàn)有人員疏散模型的理論思想多是集中于還原在某一特定場(chǎng)景下人員真實(shí)的疏散過(guò)程。本文針對(duì)于學(xué)校這一特定場(chǎng)景，結(jié)合疏散演練人員疏散過(guò)程的特殊性（有人員引導(dǎo)且有序疏散），設(shè)計(jì)了基于元胞自動(dòng)機(jī)理論基礎(chǔ)的人員調(diào)度疏散模型。在完成了教室內(nèi)人員疏散模型的搭建工作之后，對(duì)模型開(kāi)展了仿真實(shí)驗(yàn)，并且對(duì)比了無(wú)人員調(diào)度和有人員調(diào)度的仿真過(guò)程，結(jié)果顯示有人員調(diào)度的疏散模型在人員分布不均勻的情況下能更高效地將人員疏散。

劉歡，2013．建筑物內(nèi)人員疏散算法與仿真．上海：復(fù)旦大學(xué)．

岳昊，邵春福，關(guān)宏志等，2010．基于元胞自動(dòng)機(jī)的行人視線受影響的疏散流仿真研究．物理學(xué)報(bào)，59（7）：4499—4507．

張福浩，劉紀(jì)平，李青元，2004．基于Dijkstra算法的一種最短路徑優(yōu)化算法．遙感信息，（2）：38-41．

張培紅，2002．建筑物火災(zāi)時(shí)人員疏散行為的研究．沈陽(yáng)：東北大學(xué)．

張樹(shù)平，2004．建筑火災(zāi)中人的行為反應(yīng)研究．西安：西安建筑科技大學(xué)．

Burstedde C., Klauck K., Schadschneider A., et al., 2001. Simulation of pedestrian dynamics using a two-dimensional cellular automaton. Physica A: Statistical Mechanics and its Applications, 295(3-4): 507—525.

Dijkstra E. W., 1959. A note on two problems in connection with graphs. Numerische Mathematik, 1: 269—271.

Kirchner A., Schadschneider A., 2002. Simulation of evacuation processes using a bionics-inspired cellular automaton model for pedestrian dynamics. Physica A: Statistical Mechanics and its Applications, 312(1-2): 260—276.

Pauls J. L., 1978. "Evacuation of high rise office buildings." Buildings, 72(5): 84—88.

Yang L. Z., Fang W. F., Huang R., et al., 2002. Occupant evacuation model based on cellular automata in fire. Chinese Science Bulletin, 47(17): 1484—1488.

The Evacuation Model of Evacuating the Drill Personnelin Primary and Secondary School

Du Chenxi1), Lu Changjiang1),Mao Shanshan2)and Cheng Yi1)

1) Earthquake Administration of Sichuan Province, Chengdu 610200, China 2) Chengdu University of Technology, Chengdu 610059, China

On the basis of cellular automata models, we build the evacuation model of evacuating the drill personnel in primary and secondary school with the aim of shortest evacuation time. This model provides the core algorithm for earthquake emergency training evacuation training software. Through analysis of spatial foundation construction elements and construction methods, we design the virtual scene modeling method. We also construct a method of personnel evacuation algorithm and combine the research results of emergency evacuation practice personnel to complete the overall design of evacuation model. The model provides the basis for judging the rationality of personnel evacuation process in practise.

Earthquake emergency evacuation model; Cellular automata; Dijkstra's algorithm; Recursive；Evacuation strategy

10.11899/zzfy20170421

“十二五”國(guó)家科技支撐計(jì)劃項(xiàng)目“城鎮(zhèn)地震防災(zāi)與應(yīng)急處置一體化服務(wù)系統(tǒng)及其應(yīng)用示范”課題三“城鎮(zhèn)地震災(zāi)害應(yīng)急處置關(guān)鍵技術(shù)研究”（2015BAK18B03）

2017-04-12

杜晨曦，男，生于1988年。助理工程師。主要從事地震現(xiàn)場(chǎng)應(yīng)急救援工作。E-mail：184778316@qq.com

杜晨曦，魯長(zhǎng)江，毛珊珊，程奕，2017．中小學(xué)校疏散演練人員疏散模型研究．震災(zāi)防御技術(shù)，12（4）：936—946．