一种考虑疏散引导的室内多出口行人流仿真方法

摘要

本发明涉及一种考虑疏散引导的室内多出口行人流仿真方法，首先将室内空间做格网化处理，接着初始化疏散人员，以均匀随机分布将疏散行人安排在随机的空白元胞中，且以与引导人员不同颜色的圆点表示；然后搜索各元胞内的行人，解析各行人所受到的引导作用力、行人之间作用力、可见出口吸引力、障碍物斥力，确定各行人可能前往的四个方向，以此得到方向矩阵A；最后解析方向矩阵，各行人参与竞争，得到实际移动方向；更新行人位置，直至所有行人到达出口。本发明能够完善现有基于元胞自动机的行人流疏散模型，揭示多出口引导疏散场景中行人流疏散效果与引导作用之间的内在规律，为行人流疏散引导策略的评估与优化提供理论支持。

权利要求

1.一种考虑疏散引导的室内多出口行人流仿真方法，其特征在于：包括以下步骤：步骤S1：将室内空间做格网化处理，构成(W+2)·(W+2)的二维离散元胞网格系统；根据室内疏散策略对格网化后的室内空间进行场景设置，以不同颜色方格或圆点分别代表墙壁、出口、障碍物、引导人员；

步骤S2：初始化疏散人员，以均匀随机分布将疏散行人安排在随机的空白元胞中，且以与引导人员不同颜色的圆点表示；

步骤S3:搜索各元胞内的行人，解析各行人所受到的引导作用力、行人之间作用力、可见出口吸引力、障碍物斥力，确定各行人可能前往的四个方向，以此得到方向矩阵A；

步骤S4：解析方向矩阵，各行人参与竞争，得到实际移动方向；

步骤S5：更新行人位置，直至所有行人到达出口；

其中，所述步骤S3具体包括以下步骤：

步骤S31:计算引导作用力

其中，D为引导人员的最大引导系数；φ为个体引导时目标视轴与干扰视轴的夹角，群体引导时φ＝0；n为单独引导集中程度系数，该程度系数越大则单独引导效果越明显；r为正确引导人员到疏散行人的距离；r′i为其他作为干扰的引导人员到疏散行人的距离；

为引导作用力的方向，其方向从疏散行人指向正确引导人员的所在的出口；

步骤S32：计算行人之间作用力

其中， 表示行人排斥力； 表示行人吸引力；kab为行人斥力系数；gab为行人引力系数；rab为疏散行人a与疏散行人b之间的距离；dab为疏散行人a与疏散行人b之间的最小步长；w为视野范围步长；qa为行人a的带电系数；qb为行人b的带电系数； 为行人斥力的方向， 为行人斥力的方向；

步骤S33：计算障碍物斥力

其中，kwb为障碍物斥力系数；rwb为障碍物w与疏散行人b之间的距离；dwb为障碍物w与疏散行人b之间的最小步长；w为视野范围步长；qw为障碍物w的带电系数；qb为行人b的带电系数； 为障碍物斥力的方向；

步骤S34：计算可见出口吸引力

其中，E为可见出口引力系数；dexit为可见出口与疏散行人之间的最小步长；w为视野范围步长； 为可见出口吸引力的方向，其方向从疏散行人指向可见出口；

步骤S35：得到方向矩阵A：

其中，Address是该行人原先的位置，m1,m2,m3,m4∈{±1,±(W+1),±(W+2),±(W+3)}。

2.根据权利要求1所述的一种考虑疏散引导的室内多出口行人流仿真方法，其特征在于：所述步骤S4具体为：各行人按照方向矩阵A中决定参与方向竞争的顺序，吸引力大的方向优先考虑；当多个疏散人员同时竞争同一移动目标，此位置将被系统以等概率随机分配给其中一个行人作为下一时间步移动位置；当疏散行人在竞争中未被成功分配某一目标网格，则选择次优网格作为竞争方向；当疏散行人在所有可能移动网格的竞争中落选，该行人保持不动以等概率的根据竞争结果分配疏散行人下一步移动元胞。

说明书

一种考虑疏散引导的室内多出口行人流仿真方法

技术领域

[0001] 本发明涉及行人流特征研究领域，特别是一种考虑疏散引导的室内多出口行人流仿真方法。

背景技术

[0002] 复杂情况下大型集会的突发事件可能造成拥挤踩踏事故，疏散场所内的行人个体往往会因不熟悉环境而找不到最佳的逃生方向。通过行人流疏散模型来评价突发事件发生时集会疏散能力并提供有效的疏散策略是行人流领域的研究重点。

[0003] 行人流特征的研究主要由行人微观行为特征和宏观群体特征两方面组成。行人微观行为特征研究主要是针对行人个体行为及其影响因素进行分析研究，例如分析行人的行人恐惧度、行人慢启动效应、超越行为、羊群效应、团体行为效应等。行人流宏观行为特征研究主要是针对行人群体行为及其评价进行分析研究，分析研究行人群体拥挤机理以及群集行为特性等。

[0004] 目前，典型的行人流疏散模型包括元胞自动机模型、格子气模型、社会力模型、流体动力模型等。其中，元胞自动机作为一种能模拟行人个体行为并能呈现行人微观特征及宏观行为的离散模型，已经被广泛地应用于行人流仿真研究中，对于减少室内行人疏散拥堵现象，避免发生行人踩踏事故，完善集会突发情况下快速疏散的基础理论及其应用水平均具有重要作用。

[0005] Schultz结合路径查找算法及摩擦效应，利用元胞自动机模拟了机场旅客疏散情况；永贵等通过改进的层次域元胞自动机模型进行了行人疏散问题的研究，考虑了多种逃生受阻因素，并利用正菱形网格避免行人与墙壁的碰撞；陈亮等详细阐述了行人之间排斥力、路径通行能力等因素的影响，深入研究了行人疏散过程中的路径选择行为。

[0006] 然而上述基于元胞自动机的行人仿真研究中，都没有关注引导作用对行人疏散过程的影响。而在考虑紧急情况下行人的心理行为及相互作用的条件下,引导人员可有效提高公共场所的疏散效率，是大型集会疏散方案中必不可少的一部分。因此，有必要将引导人员的引导作用引入仿真模型，构建考虑多出口疏散引导情况的行人流仿真模型，量化多出口疏散场景中的引导作用力，以更为综合地考察室内环境行人流运动演化，更好地揭示多出口引导疏散场景中行人流的内在规律。

发明内容

[0007] 有鉴于此，本发明的目的是提出一种考虑疏散引导的室内多出口行人流仿真方法，完善现有基于元胞自动机的行人流疏散模型，揭示多出口引导疏散场景中行人流疏散效果与引导作用之间的内在规律，为行人流疏散引导策略的评估与优化提供理论支持。

[0008] 本发明采用以下方案实现：一种考虑疏散引导的室内多出口行人流仿真方法，具体包括以下步骤：

[0009] 步骤S1：将室内空间做格网化处理，构成(W+2)·(W+2)的二维离散元胞网格系统；根据室内疏散策略对格网化后的室内空间进行场景设置，以不同颜色方格或圆点分别代表墙壁、出口、障碍物、引导人员；

[0010] 步骤S2：初始化疏散人员，以均匀随机分布将疏散行人安排在随机的空白元胞中，且以与引导人员不同颜色的圆点表示；

[0011] 步骤S3:搜索各元胞内的行人，解析各行人所受到的引导作用力、行人之间作用力、可见出口吸引力、障碍物斥力，确定各行人可能前往的四个方向，以此得到方向矩阵A；

[0012] 步骤S4：解析方向矩阵，各行人参与竞争，得到实际移动方向；

[0013] 步骤S5：更新行人位置，直至所有行人到达出口。

[0014] 进一步地，所述步骤S3具体包括以下步骤：

[0015] 步骤S31:计算引导作用力

[0016]

[0017] 其中，D为引导人员的最大引导系数；φ为个体引导时目标视轴与干扰视轴的夹角，群体引导时φ＝0；n为单独引导集中程度系数，该程度系数越大则单独引导效果越明显；r为正确引导人员到疏散行人的距离；ri′为其他作为干扰的引导人员到疏散行人的距离； 为引导作用力的方向，其方向从疏散行人指向正确引导人员的所在的出口；

[0018] 步骤S32：计算行人之间作用力

[0019]

[0020] 其中， 表示行人排斥力； 表示行人吸引力；kab为行人斥力系数；gab为行人引力系数；rab为疏散行人a与疏散行人b之间的距离；dab为疏散行人a与疏散行人b之间的最小步长；w为视野范围步长；qa为行人a的带电系数；qw为行人b的带电系数； 为行人斥力的方向， 为行人斥力的方向；

[0021] 步骤S33：计算障碍物斥力

[0022]

[0023] 其中，kwb为障碍物斥力系数；rwb为障碍物w与疏散行人b之间的距离；dwb为障碍物w与疏散行人b之间的最小步长；w为视野范围步长；qw为障碍物w的带电系数；qb为行人b的带电系数； 为障碍物斥力的方向；

[0024] 步骤S34：计算可见出口吸引力

[0025]

[0026] 其中，E为可见出口引力系数；dexit为可见出口与疏散行人之间的最小步长；w为视野范围步长； 为可见出口吸引力的方向，其方向从疏散行人指向可见出口；

[0027] 步骤S35：得到方向矩阵A：

[0028]

[0029] 其中，Address是该行人原先的位置，m1,m2,m3,m4∈{±1,±(W+1),±(W+2),±(W+3)}。

[0030] 进一步地，所述步骤S4具体为：各行人按照方向矩阵A中决定参与方向竞争的顺序，吸引力大的方向优先考虑；当多个疏散人员同时竞争同一移动目标，此位置将被系统以等概率随机分配给其中一个行人作为下一时间步移动位置；当疏散行人在竞争中未被成功分配某一目标网格，则选择次优网格作为竞争方向；当疏散行人在所有可能移动网格的竞争中落选，该行人保持不动以等概率的根据竞争结果分配疏散行人下一步移动元胞。

[0031] 与现有技术相比，本发明有以下有益效果：本发明针对多出口疏散环境中引导人员在固定位置通过声音、动作、警示灯等途径为疏散行人提供相应疏散方向的引导方式，并考虑引导信号之间的相互影响，提出一种考虑疏散引导的室内多出口行人流仿真方法。本发明通过引入无线通信理论计算引导受信度，量化引导作用力，结合可见出口吸引力、行人吸引力、行人排斥力和障碍物排斥力进一步反映行人流的微观特征，能弥补现有离散技术对描述引导作用影响的不足，能更真实地反映人员疏散的现实状况，揭示多出口引导疏散场景中行人流疏散效果与引导作用强度的内在规律，可为行人流疏散引导策略的评估与优化提供理论支持。本发明规则简单、运算速度快。

附图说明

[0032] 图1为本发明实施例中的方法流程示意图。

[0033] 图2为本发明实施例中的Moore型元胞自动机示意图。

[0034] 图3为本发明实施例中的W＝19的室内疏散仿真环境示意图。

[0035] 图4为本发明实施例中的群体引导及个体引导的信号及干扰分析示意图。

[0036] 图5为本发明实施例中的行人之间作用力示意图。

[0037] 图6为本发明实施例中的w＝3时的行人之间的作用力示意图。

[0038] 图7为本发明实施例中的w＝3时的障碍物斥力示意图。

[0039] 图8为本发明实施例中的某疏散行人的合力分解示意图。

[0040] 图9为本发明实施例中的疏散过程模拟结果图。

[0041] 图10为本发明实施例中的不同行人密度下引导作用系数对系统疏散时间的影响仿真结果示意图。

具体实施方式

[0042] 下面结合附图及实施例对本发明做进一步说明。

[0043] 如图1所示，本实施例提供了一种考虑疏散引导的室内多出口行人流仿真方法，具体包括以下步骤：

[0044] 步骤S1：将室内空间做格网化处理，构成(W+2)·(W+2)的二维离散元胞网格系统；根据室内疏散策略对格网化后的室内空间进行场景设置，以不同颜色方格或圆点分别代表墙壁、出口、障碍物、引导人员；

[0045] 步骤S2：初始化疏散人员，以均匀随机分布将疏散行人安排在随机的空白元胞中，且以与引导人员不同颜色的圆点表示；

[0046] 步骤S3:搜索各元胞内的行人，解析各行人所受到的引导作用力、行人之间作用力、可见出口吸引力、障碍物斥力，确定各行人可能前往的四个方向，以此得到方向矩阵A；

[0047] 步骤S4：解析方向矩阵，各行人参与竞争，得到实际移动方向；

[0048] 步骤S5：更新行人位置，直至所有行人到达出口。

[0049] 下面结合图2至图10具体对本实施例进行阐述。

[0050] 第一步，将室内空间做格网化处理。采用Moore型的元胞自动机邻居结构，将疏散空间均匀划分为等宽的元胞网格，每个元胞网格视为一个基本单位，其对应的8个元胞移动方向如图2所示。任意一个元胞网格都只有三种可能状态，分别为被疏散行人占据、被障碍物占据及未被占据。选取场景W＝19并建立如图3的疏散空间，其中浅灰色方格代表墙壁，黑色方格代表出口，深灰色方格代表障碍物，浅灰色圆点代表引导人员。设置元胞尺寸为0.5m×0.5m，每个元胞只能容纳一个行人，每个出口宽度为1个元胞大小，整个疏散空间为9.5m×9.5m。出口是行人离开疏散空间的唯一途径。

[0051] 第二步，初始化疏散人员，以均匀随机分布将疏散行人安排在随机的空白元胞中，并用白色圆点表示。

[0052] 第三步，解析行人受力，计算方向矩阵A：

[0053] 1)引导作用力

[0054] 在实际的疏散过程中，引导人员的引导作用会影响行人流的疏散效果。由于疏散人员可接收到多个引导人员的疏散指令，疏散行人存在方向犹豫性。疏散信号的受信度将直接影响引导作用力的大小。本文引入无线通信中的天线发送及接收模型，把引导人员的引导信号视为发射天线，疏散人员视为接收天线，通过计算载波及干扰信号的接收功率得到正确疏散信号的受信度。

[0055] 同时，本实施例将引导作用分为群体引导和个体引导两个阶段，群体引导及个体引导中的引导人员分别视为全向天线及定向天线，如图4所示。当疏散行为从群体引导进入个体引导阶段时，疏散行人所受到的干扰进一步缩小，更易被正确的引导人员疏导。引导作用力为

[0056]

[0057] 式中：D为引导人员的最大引导系数； 为个体引导时目标视轴与干扰视轴的夹角，群体引导时 n为单独引导集中程度系数，系数越大则单独引导效果越明显；r为正确引导人员到疏散行人的距离；ri′为其他作为干扰的引导人员到疏散行人的距离； 为引导作用力的方向，其方向从疏散行人指向正确引导人员的所在的出口。

[0058] 2)行人之间作用力

[0059] 疏散行人将产生行人之间作用力，并引起行人流的羊群效应及冲突效应。其大小与方向取决于疏散行人之间的距离，并包括行人排斥力、行人吸引力、无作用力三种情况，如图5所示。行人之间作用力为

[0060]

[0061] 式中：kab为行人斥力系数；gab为行人引力系数；rab为疏散行人a与疏散行人b之间的距离；dab为疏散行人a与疏散行人b之间的最小步长；w为视野范围步长；qa为行人a的带电系数；qw为行人b的带电系数； 为行人斥力的方向， 为行人斥力的方向，其方向如图6所示。

[0062] 3)障碍物斥力

[0063] 障碍物斥力为

[0064]

[0065] 式中：kwb为障碍物斥力系数；rwb为障碍物w与疏散行人b之间的距离；dwb为障碍物w与疏散行人b之间的最小步长；w为视野范围步长；qw为障碍物w的带电系数；qb为行人b的带电系数； 为障碍物斥力的方向，其方向如图7所示。

[0066] 4)可见出口吸引力

[0067] 当疏散行人位于某处元胞网格，且其可视范围内包含出口时，则该出口即可见出口，并产生可见出口吸引力。当人离开疏散空间或出口不可见，将不存在可见出口吸引力。可见出口吸引力为

[0068]

[0069] 式中：E为可见出口引力系数；dexit为可见出口与疏散行人之间的最小步长；w为视野范围步长； 为可见出口吸引力的方向，其方向从疏散行人指向可见出口。

[0070] 5)方向矩阵

[0071] 方向矩阵的元素是各力的合力沿邻居方向的分力方向所对应的元胞位置，并以分力大小依次排列。例如图8，某行人所受合力 的四个分力为垂直向上、水平向右、右上、右下四个方向，且 因此方向矩阵为

[0072]

[0073] 第四步，解析方向矩阵，各行人参与竞争，得到实际移动方向。各行人按照方向矩阵A中决定参与方向竞争的顺序，吸引力大的方向优先考虑；当多个疏散人员同时竞争同一移动目标，此位置将被系统以等概率随机分配给其中一个行人作为下一时间步移动位置；当疏散行人在竞争中未被成功分配某一目标网格，则选择次优网格作为竞争方向；当疏散行人在所有可能移动网格的竞争中落选，该行人保持不动以等概率的根据竞争结果分配疏散行人下一步移动元胞。

[0074] 第五步更新行人位置，回到第三步，行人继续在疏散空间内运动，直至所有行人到达出口完成疏散。

[0075] 本方法的疏散模拟过程结果如图9所示，可揭示多出口引导疏散场景中行人流疏散效果与引导作用强度的内在规律，仿真结果如图10所示。

[0076] 以上所述仅为本发明的较佳实施例，凡依本发明申请专利范围所做的均等变化与修饰，皆应属本发明的涵盖范围。