中小学校疏散演练人员疏散模型研究1

杜晨曦 鲁长江 毛珊珊 程 奕



中小学校疏散演练人员疏散模型研究1

杜晨曦1）鲁长江1）毛珊珊2）程 奕1）

1）四川省地震局，成都 610200 2）成都理工大学，成都 610059

本文以元胞自动机模型理论为指导，以人员疏散时间最短为目标，构建中小学校疏散演练人员疏散模型。本文从开展地震应急疏散模型研究出发，为地震应急演练疏散训练软件提供核心算法。通过对空间基础构建元素及搭建方法的分析研究，设计完成疏散虚拟场景的建模方法。最后开展人员疏散算法构建方式的研究，并结合应急疏散演练人员行为规范的研究成果，完成疏散模型的总体设计，达到为应急演练疏散软件提供人员疏散演化模型以及为人员疏散过程合理性评判提供依据的目的。

地震应急疏散模型 元胞自动机 Dijkstra算法 递推 疏散策略

引言

地震是群灾之首，全球地震灾害多发，给人类造成了巨大损失。我国属于地震灾害多发国家，频发的地震灾害对中国经济建设、社会稳定和城市化进程等构成了严重威胁。特别是2008年汶川地震和2010年玉树地震，造成了严重的人员伤亡和财产损失。而在汶川地震中四川省安县桑枣中学2000多名师生在1分36秒全体成功逃生又为我们提供了宝贵的经验，即科学、合理的日常疏散演练能大大提高地震逃生概率，为生命安全提供有力保障。

因此，开展地震应急疏散研究，特别是针对中小学的地震应急疏散研究显得尤为重要。全面分析和研究地震应急疏散模型可以为减少地震灾害造成的人员和财产损失打下良好基础。

目前，国内外在人员疏散模型的研究上主要采用连续和离散两类模拟方法。其中连续模型多是把行人运动看作气体、水流、颗粒流等物理介质，借鉴一些物理方程来建立模型。如Pauls（1978）通过对火灾情况下人员疏散演习实测得出的多层建筑人员疏散时间经验公式。而离散模型则把时间和空间划分为不连续的“时间步”和“网格”，通过一定局部或整体的行为运动规律来模拟人员疏散过程，代表模型为元胞自动机。国内对人员疏散的研究起步较晚，但经过近几年的努力，已经取得了不少成绩。张培红（2002）在东北大学校园等处对人员流动状态进行了实际观测，从中归纳出了弯道90°以及下楼梯时人员流动速度衰减系数，建立了离散状态人员疏散行动的计算机模拟模型SHEBR-SGEM；张树平（2004）调查了一百多起火灾，对火灾发生时人员所处的环境状态、人员疏散行为等进行了理论研究；中国科学技术大学的杨立中教授等利用元胞自动机模型对疏散过程中的从众行为、亲情行为、建筑物出口动力学特征、信息传播等进行了深入的探讨（Yang等，2002）。

本文所提及的人员疏散模型是地震应急疏散演练训练器软件的核心组成部件之一。地震应急疏散演练训练器软件的应用范围是中小学校，软件的研制目的是帮助中小学校通过计算机模拟的方式快速得出适合于本校实际情况的疏散演练方案，学校基于疏散演练方案开展震前疏散演练训练，达到有序、高效、安全疏散的目的，做好震前疏散准备工作，一旦发生地震，学校师生可以按照优化后的疏散演练方案有序、快速地撤离。

1 模型的建立

本文基于城镇中小学疏散场景模型，通过对迪科斯彻算法、元胞自动机疏散模型以及人员调度策略的研究分析，结合人员分布情况及建筑物结构特点，得出适合中小学校人员的疏散模型。

1.1 疏散场景模型的建立

把已有的疏散软件及常用建模软件建筑物内部空间的搭建方式，与四川省中小学校典型建筑结构类型及内部空间特点的走访调研数据相结合，得出合理可行的疏散场景模型搭建方法。

1.1.1 疏散场景模型的分类

一个疏散场景可用2种方法表示，即精细网络模型和粗糙网络模型。精细网络模型使用大量的网络节点表示一个建筑物空间，每个节点对应一小块面积的疏散空间，并且允许在疏散模型中对其进行详细的描述。这种模型通常是对建筑物几何的精确描述，并且在建模时通常需要建筑物的图纸。粗糙网络模型使用节点和与之相连接的边来表示一个疏散场景，一条边可以表示疏散场景内一个真实的连接关系，节点代表建筑物空间中一个独立的疏散空间，比如一个教室、一个办公室等。

1.1.2 疏散场景模型的搭建

通过对四川省中小学校的走访调研选定建立模型的基础环境，根据人员所在建筑物的位置和影响人员快速撤离的空间结构等因素，梳理出疏散场景搭建所需的点、线、面以及障碍物、楼梯等构建元素类型，分析得出精细网络模型是更为合理的疏散场景模型搭建方法。模型搭建步骤如下：

（1）通过自主研发编辑器编辑空间模型（未使用AutoCAD、Visio等软件）。

（2）空间模型搭建由整体到局部、由框架到细节。

图1 疏散场景模型的搭建流程

通过精细网络模型搭建的教学楼空间结构图如图2所示。

图2 教学楼空间结构示意图

1.2 人员疏散模型的建立

建筑物内的疏散人员可以从宏观和微观两个角度表示。从宏观角度，疏散人员可看作是一个群体，在这个群体中每个疏散人员是完全相同的，按照同样的行为进行疏散，没有个体的特质。从微观角度，每个疏散人员都会有自己的个体特征，或者被随机分配不同的属性，在疏散过程中，每个疏散人员在不同的条件下可能表现出不同的疏散行为。近年来，在建筑物人员疏散的微观模型中，元胞自动机模型以其简单性、灵活性和高效性在疏散流仿真方面受到越来越多的关注。

我们结合Dijkstra算法、人员调度策略以及中小学教学楼这一特定的疏散环境对疏散流进行建模仿真。

1.2.1 元胞自动机模型原理

元胞自动机模型（Kirchner，2002；Burstedde，2001；岳昊等，2010）是在一个由离散、有限状态的元胞组成的元胞空间中，按照一定的局部规则，在离散的时间维上演化的动力学系统。元胞自动机系统状态用一种常规的格子组成单元。每个单元可以处于几种状态中的一种（一般有两种，0和1），元胞的状态为0时表示空出，为1时表示有人或障碍物，采用二维Von-Neumann邻居或Moore邻居。

元胞自动机模型通常采用并行更新规则，在同一时步，系统中的所有行人同步并行更新运动到相对应的格点上。若出现多位行人同时选择一个格点，则发生冲突，解决的办法通常是采用等概率随机地选择其中一位到达该格点，其余保持在原来位置，等到下一个时间步再做出选择。

1.2.2 Dijkstra算法原理

Dijkstra算法（Dijkstra's algorithm）由荷兰计算机科学家Dijkstra（1959）提出。该算法使用了广度优先搜索解决非负权有向图的单源最短路径问题，算法最终得到一个最短路径树。该算法常用于路由算法或者作为其他图算法的一个子模块。

图3 Dijkstra算法示意图

第三步 D尚未着色，继续执行第二步。

第三步 D尚未着色，继续执行第二步。

第三步 D尚未着色，继续执行第二步。

第三步 D尚未着色，继续执行第二步。

第三步 D尚未着色，继续执行第二步。

第三步 D尚未着色，继续执行第二步。

得到从A到D的最短距离10。

1.2.3 人员疏散的前提条件

人员在疏散时遵循以下前提条件：

（1）网格尺寸为：0.4m×0.4m。（相关文献资料显示我国人体平均体厚与肩宽投影面积为0.146m2，且元胞自动机网格化通用网格大小划分为0.4m×0.4m）。

（2）每个元胞有8个运动方向，每一个时间步内人员只能移动一个网格。

（4）有疏散人员的节点称为疏散源。每次递推时每个疏散人员可能位于不同的疏散源，所以到每个出口的最短距离不同。

1.2.4 教室内疏散人员模型设计

教室内疏散人员模型规则：

（1）利用迪科斯切算法求出出口到每个格子的距离。

（2）利用元胞自动机思想，在格点上的元胞选取离出口距离最小的格点前进。

（3）如果离出口距离最小的格点位置被占据，元胞则自动选择次好位置，但绝不能走差位置（即距离出口位置比当前位置还远的格点）。

（4）移动时，先让在出口位置的元胞移动，离开教室，再以出口为起点，由近及远调度其他元胞移动。

（5）在每个元胞等待移动时，均需计算等待时间，调度时优先调度等待时间长的元胞，让该元胞选择移动格点。

（6）在元胞移动时，为每个元胞增加一个行动力，每次加1，总的行动力不超过2，每走一步数据更新一次，以此解决速度问题。

（7）把人员所在位置到出口的距离与教室内总人数相结合，对教室内人员选择不同出口进行总体调度。

需要注意的是，教室内疏散模型设计的第（7）步是为了对前面6步进行优化而提出的调度策略。在模型仿真实验中，我们把没有第（7）步（即未优化版本）和有第（7）步（即优化版本）的结果进行对比，进一步证明调度策略的作用。

1.2.5 整栋教学楼疏散人员模型设计

（3）将疏散时间数组，疏散人员分配矩阵，中间结果数组初始化成0：

（4）计算每次递推时每个出口的疏散时间的中间值：

（6）更新每个出口对应的疏散时间，完成第次递推：

图4 疏散模型程序流程图

2 仿真

本文仅针对一个教室内的人员进行疏散模型仿真实验。仿真分为未优化版本和优化版本，每个版本又分别考虑了人员均匀和不均匀分布情况（表1）。每次仿真实验，都进行了相关数据采集，下面将对各仿真结果进行说明。

表1 不同版本疏散完成数据的对比

（1）未优化版本的人员均匀分布的疏散

如图5所示，疏散开始前，教室内人员均匀分布。疏散开始后，教室内的人员开始分别往左右出口移动，到疏散即将完成时，左右出口的人数相差很小，疏散基本同时完成。

图5 未优化版本的人员均匀分布情况

（2）未优化版本的人员不均匀分布的疏散

如图6所示，疏散开始前，教室内人员是不均匀分布的。疏散开始后，教室内的人员开始分别往左右出口移动，到疏散即将完成时，左右出口的人数相差很大，很明显，在这种情况下，右边出口拥堵情况变得更加严重。疏散步数相对人员均匀分布的情况有所增加。

图6 未优化版本的人员不均匀分布情况

（3）优化版本的人员均匀分布的疏散

在优化版本里加入教室内人员调度模型，把疏散人员数量与出口的距离结合后进行算法优化，优化后的程序有一个控制台窗口，其中记录了给每个出口分配的人数以及总步数。如图7所示，疏散开始前，教室内人员是均匀分布的。疏散开始后，教室内的人员开始分别往左右出口移动，到疏散即将完成时，左右出口的人数相差很小，疏散基本同时完成。

图7 优化版本的人员均匀分布情况

（4）优化版本的人员不均匀分布的疏散

优化后的不均匀分布状态，如图8所示。疏散开始前，教室内人员是不均匀分布的。疏散开始后，教室内的人员开始分别往左右出口移动，到疏散即将完成时，很明显，左右两个出口的分布更加平均，相对未优化的版本疏散完成步数减少6步，说明经优化后疏散时间减少。

图8 优化版本的人员不均匀分布情况

3 结论

根据本文的研究目标有针对性地对四川省部分中小学校进行了调研，在收集基础资料的同时也了解了学校疏散工作开展情况。

国外现有人员疏散模型的理论思想多是集中于还原在某一特定场景下人员真实的疏散过程。本文针对于学校这一特定场景，结合疏散演练人员疏散过程的特殊性（有人员引导且有序疏散），设计了基于元胞自动机理论基础的人员调度疏散模型。在完成了教室内人员疏散模型的搭建工作之后，对模型开展了仿真实验，并且对比了无人员调度和有人员调度的仿真过程，结果显示有人员调度的疏散模型在人员分布不均匀的情况下能更高效地将人员疏散。

刘欢，2013．建筑物内人员疏散算法与仿真．上海：复旦大学．

岳昊，邵春福，关宏志等，2010．基于元胞自动机的行人视线受影响的疏散流仿真研究．物理学报，59（7）：4499—4507．

张福浩，刘纪平，李青元，2004．基于Dijkstra算法的一种最短路径优化算法．遥感信息，（2）：38-41．

张培红，2002．建筑物火灾时人员疏散行为的研究．沈阳：东北大学．

张树平，2004．建筑火灾中人的行为反应研究．西安：西安建筑科技大学．

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

“十二五”国家科技支撑计划项目“城镇地震防灾与应急处置一体化服务系统及其应用示范”课题三“城镇地震灾害应急处置关键技术研究”（2015BAK18B03）

2017-04-12

杜晨曦，男，生于1988年。助理工程师。主要从事地震现场应急救援工作。E-mail：184778316@qq.com

杜晨曦，鲁长江，毛珊珊，程奕，2017．中小学校疏散演练人员疏散模型研究．震灾防御技术，12（4）：936—946．