一种基于听觉感知的空间人群疏散模拟方法
【专利摘要】一种基于听觉感知的空间人群疏散模拟方法，本发明涉及基于听觉感知的空间人群疏散模拟方法。本发明是为了解决现有疏散模拟方法中没有考虑听觉要素如警报器的数量和位置的影响的问题。本发明步骤为：步骤一、将建筑物平面进行网格划分，划分后的网格分为两种情况，一种为被人或障碍物占据，一种为空；步骤二、设置网格中人员移动速度和方向；步骤三、依据听觉感知元胞自动机模型得到人员i在下一时间步移动到网格j的概率；步骤四、对网格中所有人员的移动概率进行计算，对人员在当前时间位置至到达出口的整个疏散过程进行模拟。本发明应用于建筑安全性评价领域。
【专利说明】
一种基于听觉感知的空间人群疏散模拟方法
技术领域
[0001 ]本发明涉及基于听觉感知的空间人群疏散模拟方法。
【背景技术】
[0002] 目前，人员疏散模型的建模方法大体上可分为两种:一种是宏观的方法，即把行人 视为连续流动介质，因为现代人员疏散研究是从交通流的研究中分化出来的，因而也就很 自然地继承了流体研究中己经完善和成熟的方法。最早的宏观模型是由Henderson提出的， 他认为行人的运动行为类似于气体或液体的流动，行人行为的气态动力学方程与 Boltzmann方程相似，不过它考虑了行人之间的相互影响和行人的目的。Hughes采用连续介 质理论(continuum theory)研究大型人群的运动特征，并根据Navier-Stokes方程进一步 推导出大型人群流动的控制方程和行人避免向高密度人群运动的方程，该模型成功地解释 了麦加朝圣中的人群的运动状况。然而人群是不遵守动量和能量守恒的，宏观模型的缺点 在于没有将人员的相互作用考虑在内，因此不适于紧急情况下对行人流进行研究。宏观模 型忽略了个体之间的差异，于是研究者又提出一种新的模型，即微观模型。
[0003] 另一种微观的方法，它把行人视为相互作用的粒子，其中最为著名的就是Helbing 的社会力模型。微观模型可以对行人流的具体行为进行描述，近年来引起了极大的关注。该 模型的研究方法主要有连续型和离散型两种。其中具有代表性的模型有社会力模型、磁场 力模型、元胞自动机模型等和格子气(Lattice Gas)模型。
【发明内容】
[0004] 本发明是为了解决现有疏散方法中没有考虑听觉要素如警报器的数量和位置的 影响的问题，而提出的一种基于听觉感知的空间人群疏散模拟方法。
[0005] -种基于听觉感知的空间人群疏散模拟方法按以下步骤实现：
[0006] 步骤一:将建筑物平面进行网格划分，划分后的网格分为两种情况，一种为被人或 障碍物占据，一种为空；
[0007] 步骤二:设置网格中人员移动速度和方向；
[0008 ]步骤三:依据听觉感知元胞自动机模型得到人员i在下一时间步移动到网格j的概 率为：
[0010 ] nmajPn j为网格j最大容纳人数和当前已容纳人数，nmax = n j时表示网格j被障碍物 占据，ks为静态信息吸引力系数，kD为动态信息吸引力系数；
[0011] Si, j为网格j对人员i的静态信息吸引力，Di,j为网格j对人员i的动态信息吸引力，Μ 为建筑物平面分割的网格数；
[0012] 步骤四、对网格中所有人员的移动概率进行计算，对人员在当前时间位置至到达 出口的整个疏散过程进行模拟。
[0013] 发明效果：
[0014] 本发明提出了基于听觉感知的疏散模型，模型原型为元胞自动机模型。本发明考 虑了听觉要素如警报器的数量和位置的影响，并通过声场测试和疏散实验来验证确定模型 的变量。模型的精确度通过体育场内观测实验来比对得到确定。模型不仅能用来预计疏散 时间，而且还能帮助分析和改进建筑的扬声器布局对疏散的影响，评估建筑在紧急情况发 生时人员疏散的安全性，利用模型给出相应的疏散设计和预警方案。本发明算法与测量值 的误差小于10%。
【附图说明】
[0015] 图1为听觉感知元胞自动机模型移动概率分布示意图；
[0016] 图2为疏散观测实验出口位置示意图；
[0017] 图3为各出口疏散人数方差比对图；
[0018]图4为疏散初始人员位置不意图；
[0019]图5为疏散进彳丁至20s时的人员位置不意图；
[0020]图6为疏散进彳丁至60s时的人员位置不意图；
[0021 ]图7为疏散进行至100s时的人员位置示意图；
[0022]图8为疏散进彳丁至140s时的人员位置不意图；
[0023]图9为疏散进彳丁至180s时的人员位置不意图；
[0024] 图10为计算机模拟实验中通过各个出口的人员数量的方差对比图；
[0025] 图11为单位时间内通过各个出口的疏散人数平均值图。
【具体实施方式】
【具体实施方式】 [0026] 一:一种基于听觉感知的空间人群疏散模拟方法包括以下步骤：
[0027] 步骤一:将建筑物平面进行网格划分，划分后的网格分为两种情况，一种为被人或 障碍物占据，一种为空；
[0028] 模型基于同元胞自动机基本模型，对建筑物的平面进行均匀的网格划分，每个网 格被障碍物占据、或被人员占据、或为空。元胞自动机模型中每个网格对应空间为边长为 〇.5m的正方形，但应用在大空间建筑内时，网格大小可根据建筑物尺寸做出相应调整，每个 网格也可以根据人员密度不再限定为一人，但需要注意的是为了简化计算处在同一个网格 (元胞)中的人员彼此不受影响即可。人员每次移动最多只能行进一格。
[0029] 步骤二:设置网格中人员移动速度和方向；
[0030] 步骤三:依据听觉感知元胞自动机模型得到人员i在下一时间步移动到网格j的概 率为：
[0032]式中Pi,』表示人员i在下一时间步移动到网格j的概率，在人员所处邻域内有
[0033] nmax和η」为网格j最大容纳人数和当前已容纳人数，nmax = nj时表示网格j被墙壁等 障碍物占据，ks为静态信息吸引力系数，kD为动态信息吸引力系数；
[0034] Si, j为网格j对人员i的静态信息吸引力，Di,j为网格j对人员i的动态信息吸引力，Μ 为建筑物平面分割的网格数，根据建筑物平面尺寸确定。
[0035]步骤四、对网格中所有人员的移动概率进行计算，对人员在当前时间位置至到达 出口的整个疏散过程进行模拟。
[0036]【具体实施方式】二:本实施方式与【具体实施方式】一不同的是:所述网格为0.5m的正 方形，每个网格中大于等于1人小于等于5人。其它步骤及参数与【具体实施方式】一相同。
[0037] 其它步骤及参数与【具体实施方式】一相同。
【具体实施方式】 [0038] 三:本实施方式与一或二不同的是:所述步骤二中网 格中人员移动速度和方向具体为：
[0039] 在模型中采用Von Neumann邻域，即人员能向自身周围前后左右四个方向运动，如 图1所示。
[0040] 设置人员移动速度为lm/s~2m/s，人员在所在位置向6土_、Pij3、四个 位置移动。
[0041] 人正常行走的速度约为lm/s，但在紧急情况下快速行走速度可达1.5m/s，而小跑 的速度可以达到2m/s。另外个人的身体素质差异较大，年轻人的行动速度高于老人，男性高 于女性;行李和随行的老人孩子都会使人员行动能力减慢。为了体现这些差别，模型里人员 的最大速度可以设置两档：2m/s和lm/s。分别代表一个时间步移动一次和两个时间步移动 一次。
[0042] 其它步骤及参数与【具体实施方式】一或二相同。
[0043]【具体实施方式】四：本实施方式与【具体实施方式】一至三之一不同的是：所述步骤三 中kD具体为：
[0044] 高能见度无声音信息指令情况下kD = 0;高能见度有声音信息指令情况下kD = 〇. 84;低能见度无声音信息指令情况下kD= 1.43;低能见度有声音信息指令情况下kD= 1.2。 高能见度为大于等于5米，低能见度为小于5米。
[0045] 其它步骤及参数与【具体实施方式】一至三之一相同。
【具体实施方式】 [0046] 五:本实施方式与一至四之一不同的是:所述步骤三 中匕具体为：
[0047] 当网格j是出口时，ks=l表示该因素对疏散个体有吸引力；当网格j是墙壁等阻碍 物时，ks = -l表示该因素对疏散个体有排斥力;ks = 0表示不考虑该因素。
[0048] 其它步骤及参数与【具体实施方式】一至四之一相同。
【具体实施方式】 [0049] 六:本实施方式与一至五之一不同的是:所述步骤三 中Si, j的表达式为：
[0050] hi = ai,ek^J
[0051 ] 为阻碍系数，，表示人员i到出口 ％这段距离中墙和楼梯对人员的阻碍程度；
[0053] α和β为阻碍系数，η是网格j到出口 ek这段距离中的平地长度，rs是网格j到出口 ek 这段距离中的楼梯长度，元胞所在位置距离出口越近，对疏散人员的吸引力就越大，是 网格j到出口 ％的距离，采用坐标差值法通过网格j到出口改的距离来确定；
[0054] ri,j 表示为：
[0055] ru = max(ij^mme.[\xj - xe.\ + |y； - tmneii\\x} - xc,;\ + |y； -yeJ]
[0056] 设建筑物中共有k个出口，式中-Λ·,)表示出口改的位置坐标;（Xj，yj)表示网格j 的位置坐标；
[0057]式中等式右侧第一项表示所有格点距各出口最远距离，右侧第二项表示元胞j到 各出口的最小值。这样能够保证距离最近出口最远处的元胞对人员的吸引力为0。
[0058]其它步骤及参数与【具体实施方式】一至五之一相同。
[0059]【具体实施方式】七:本实施方式与【具体实施方式】一至六之一不同的是：所述步骤三 中Di,j的表达式为：
[0060] ^i,j _ b.i.j · fi.'w.ifi)
[0061] bi,j体现影响人员对报警信号做出的反应，其取值采用Weibull分布；表示 人员i在网格j内所受到的噪音干扰，表示为：
[0062] fi.Wk(t) = gLj(t) [θυ ·
[0063] gi, j (t)为能见度系数，0 <gi, j (t) < 1，使用 Gauss ian 分布：
[0065] 为人员i在当前时刻能够接收到的警报声信息的权重，警报声信息权重取决 于当前位置的声压级LP，警报器为点声源，声功率级为LW，LP为当前位置的声压级，当前位置 声压级随测点距声源距离的变化为：
[0067]其中，Q是指向性因子，r是声源距测点距离，R是房间常数， ，其中s是室内总 表面积，S室内平均吸声系数;因此。fc(〇表示为：
[0069] 其它步骤及参数与【具体实施方式】一至六之一相同。
[0070] 实施例一：
[0071] 实验在哈尔滨工业大学体育馆中进行。体育馆一层东西两侧共有A、B、C、D四个出 口用来疏散，如图2所示。本次实验共有1127名实验对象参加，均匀分布在一层体育场大厅 内。根据经验疏散开始后的预判时间设置为l〇s，表示在接收到疏散指令到开始疏散的时间 是 10s 〇
[0072] 通过各出口疏散的人数以及各出口疏散时间被记录下来，疏散结果表明，总疏散 用时是3分52秒，每个出口成功疏散人数非常接近。A出口疏散270人，疏散时间3分40秒;B出 口疏散284人，疏散时间3分10秒;C出口疏散330人，疏散时间4分;D出口疏散274人，疏散时 间3分30秒。出口B时间略短于其他出口，出口D的波动性显示高于其他出口。图3显示了各个 出口通行人员数量的方差。虽然每个出口的疏散过程是相似的，但其疏散能力是不均衡的。 疏散过程分为四个阶段。第一阶段是疏散开始至20s，通过出口的人数较少，但随着时间逐 渐增加。随着越来越多的人往出口方向聚集以及出口通行能力的限制，疏散人数增加的趋 势变得缓慢，达到3人/秒后不再增加。第三阶段是100s以后，疏散人员数量的方差开始保持 不变，接近1.5人/秒。180s后进入第四阶段，疏散人员数量开始下降直至疏散过程结束。
[0073] 计算机模拟的实验数据录入与疏散实验相同，首先用CAD画出体育馆平面图，接下 来录入人员位置，人数1127人平均分布在体育馆一层，预判时间设置为10s，人员移动速度 lm/s。图4-图9显示了利用听觉感知元胞自动机模型模拟不同时间段的疏散过程，分别是疏 散过程刚开始和进行到20s、60s、100s、140s和180s时的人员位置。
[0074] 从过程图中可以看出，各个出口疏散的人数基本平均。疏散总用时为3分19秒。模 拟结果同时还得出了每10秒内各个出口通行的人数已经各出口的总通行人数，详见表1。 [0075]表1计算机模拟实验每10s内各出口通行人数
[0076]
[0077]图10显示了模拟实验中通过各个出口的人员数量的方差，这说明了各个出口疏散 的人员数量越接近平均值，疏散过程所用的时间越短。疏散过程是在疏散的起始时间疏散 的人数很低，经过一个快速的上升期后持续一段时间，再平稳的下降。
[0078] 观测实验所得的数据，每十秒疏散人数的方差被计算用以与听觉感知元胞自动机 模型模拟所得数据进行对比。数值小于1的占22%，数值大于1小于3的占21%，数值大于3小 于5的占27%，数值大于5小于10的占26%，而数值大于10的为3%。最大值是15,最小值是0 (意味着等于平均值）。方差值大于5的说明在该单位时间内疏散的人数明显的不同于平均 值。观测实验与模拟实验结果的最大差值用以与平均疏散人数进行比较，其所占百分比越 小表示精确度越高。通过出口 A的人数最大差值出现在190s，差值人数为9人，占平均疏散人 数的9.4%。通过出口 B的人数中最大差值出现在80s，差值人数为7人，占平均疏散人数的 15.2 %。通过出口 C的人数中最大差值出现在120s，差值人数为10人，占平均疏散人数的 11.1 %。通过出口 D的人数中，最大差值出现在40s，差值人数为15人，占平均疏散人数的
[0079] 图11所示为观测实验和模拟实验中单位时间内通过各个出口的疏散人数平均值 比较，最大差值出现在100s，为25人，所占平均疏散人数的10.8%。其他值较为接近，差值低 于所占平均疏散人数的10%的比例为80%。说明利用听觉感知元胞自动机模型模拟的各个 出口的模拟结果接近平均值。听觉感知元胞自动机模拟的结果中，模拟显示的疏散过程在 20s以前和160s以后与观测实验的疏散过程是极为相近的。其区别在于元胞自动机没能很 好的模拟在出口处出现拥堵的情况，出现较大差距的地方主要原因是对于瓶颈现象缺乏考 虑造成的。然而，这些差距是在一个可接受的范围内。因此利用听觉感知元胞自动机模型的 模拟结果是可信的。
【主权项】
1. 一种基于听觉感知的空间人群疏散模拟方法，其特征在于，所述基于听觉感知的空 间人群疏散方法包括W下步骤： 步骤一:将建筑物平面进行网格划分，划分后的网格分为两种情况，一种为被人或障碍 物占据，一种为空； 步骤二:设置网格中人员移动速度和方向； 步骤Ξ:依据听觉感知元胞自动机模型得到人员i在下一时间步移动到网格j的概率 为：nmax和η j为网格j最大容纳人数和当前已容纳人数，nmax = η j时表示网格j被障碍物占据， ks为静态信息吸引力系数，kD为动态信息吸引力系数； Si, j为网格j对人员i的静态信息吸引力，Di,j为网格j对人员i的动态信息吸引力，Μ为建 筑物平面分割的网格数； 步骤四、对网格中所有人员的移动概率进行计算，对人员在当前时间位置至到达出口 的整个疏散过程进行模拟。2. 根据权利要求1所述的一种基于听觉感知的空间人群疏散模拟方法，其特征在于，所 述网格为0.5m的正方形，每个网格中大于等于1人小于等于5人。3. 根据权利要求2所述的一种基于听觉感知的空间人群疏散模拟方法，其特征在于，所 述步骤二中网格中人员移动速度和方向具体为： 设置人员移动速度为Im/s~2m/s，人员在巧,Λ所在位置向巧知、％、巧成和％四个位置 移动。4. 根据权利要求3所述的一种基于听觉感知的空间人群疏散模拟方法，其特征在于，所 述步骤Ξ中kD具体为： 高能见度无声音信息指令情况下kD = 0;高能见度有声音信息指令情况下kD = 0.84;低 能见度无声音信息指令情况下kD= 1.43;低能见度有声音信息指令情况下kD= 1.2。5. 根据权利要求4所述的一种基于听觉感知的空间人群疏散模拟方法，其特征在于，所 述步骤Ξ中ks具体为： 当网格j是出口时，ks = 1表示对疏散个体有吸引力；当网格j是被阻碍物占据时，ks = -1 表示对疏散个体有排斥力;ks = 0表示不考虑静态信息吸引力。6. 根据权利要求5所述的一种基于听觉感知的空间人群疏散模拟方法，其特征在于， 所述步骤Ξ中Sij的表达式为：awfc为阻碍系数，，表示人员i到出口 ek的距离中墙和楼梯对人员的阻碍程度；α和峽J阻碍系数，rf是网格j到出口 ek的平地长度，。是网格j到出口 ek的楼梯长度，ri,j 是网格j到出口 ek的距离，采用坐标差值法通过网格巧Ij出口 ek的距离来确定； ri,j表示为：设建筑物中共有k个出口，式中(Λ-.τ,)表示出口 ek的位置坐标；（Xj，yj)表示网格j的位 置坐标。7.根据权利要求6所述的一种基于听觉感知的空间人群疏散模拟方法，其特征在于，所 述步骤Ξ中Di,j的表达式为：bi,j表示人员对报警信号做出的反应，其取值采用We化ull分布;孩wfc货)表示人员i在网 格j内所受到的噪音干扰，表示为：gi,j(t)为能见度系数，〇<gi,j(t)<l，使用Gaussian分布：kwh脚为人员i在当前时刻能够接收到的警报声信息的权重，表示为：U为声功率级，Q是指向性因子，r是声源距测点距离，R是房间常数，R = || :，其中S是室 内总表面积，S室内平均吸声系数。
【文档编号】G06F17/50GK106096168SQ201610452004
【公开日】2016年11月9日
【申请日】2016年6月21日
【发明人】武悦, 康健, 张姗姗, 朱丽玮, 王超, 薛明辉, 付本臣
【申请人】哈尔滨工业大学