基于动态情绪驱动的人员应急疏散模型

2021-11-25李佳洋唐加福郭恩铭

东北大学学报（自然科学版） 2021年11期

李佳洋, 唐加福,2, 郭恩铭

(1. 东北大学工商管理学院, 辽宁沈阳 110819; 2. 东北财经大学管理科学与工程学院, 辽宁大连 116025; 3. 东北大学计算机科学与工程学院, 辽宁沈阳 110169)

人员动态疏散模型是应急疏散研究的重要工具,模拟突发事件的真实情境是评估模型的一个重要指标[1-2].在微观模型中,元胞自动机(cellular automata,CA)具有模拟复杂系统时空动态演变的能力,并且原理简单、计算效率高,因此在运动主体动态行为分析等方面得到了广泛应用. Burstedde等[3]提出“静态场”和“动态场”的概念,由此构建场域CA模型用以研究个体的疏散特征. Zheng等[4]提出了一种扩展的CA模型以分析三维烟雾扩散对个体疏散行为的影响. Zhang等[5]考虑行人对距离的认知差异并对背景场规则进行修正,构建了基于背景场的粒子群优化算法的出口选择策略.

情绪等心理状态在个体行为决策中起着非常重要的作用,从而影响个体运动.近年来,研究人员逐渐关注疏散模型如何表达个体的情绪状态. Helbing等[6]建立了一个模拟个体在紧急情况下紧张情绪和逃离动态行为的模型. Zheng等[7]提出了一种改进的地板场CA模型,分析了危机的扩散和个体情绪对疏散行为的影响. Cocking等[8]通过实证分析指出,紧急疏散时情绪能在人群中快速传播,即便是陌生人之间依旧可以传播.情绪的模型化是一个难题,如何表述情绪的更迭和动态演化,并将情绪状态的动态性嵌入行人运动模型,是目前疏散模型构建的一个难题.Mao等[9]从研究同伴行为互动的角度,提出了一个基于亲密关系的情绪传染模型来模拟紧急疏散中的同伴决策机制. Fu等[10]将流行病学SIR(susceptible-infectious-removal)模型与CA结合,认为在某种程度上,情绪传染类似于疾病的传播,每个个体可能处于三种状态之一:易感(S)、感染(I)或恢复(R),据此提出了一个新的情绪传染模型.梳理现有文献,情绪传播及更新是一个复杂的过程,当前疏散建模中考虑情绪因素的动态疏散行为方面研究是有限的,情绪的动态更新和行为决策结合的模型仍需要研讨.

在突发事件下,情绪会影响人们的判断和行为决策[11].而传统传染病SIS(susceptible-infectious-susceptible)模型考虑到个体很难达到完全免疫,仅存在易感者和感染者这两种人群,更加符合行人的情绪传播特性,因此本文在SIS模型基础上提出一种改进的情绪传染模型.本文基于SIS传染模型定义疏散个体状态,考虑个体速度变化和方向感知域,构建一种符合行人疏散特征的情绪更新机制,并将情绪传播动力学与行人运动结合,改进转移概率,建立动态情绪驱动的行人运动模型,最后对模型进行疏散过程仿真,分析微观情绪传染演化模型、情绪传染机理对总疏散时间和安全疏散人数的影响.该模型有助于应急管理者了解紧张情绪在人与人之间传播特征,便于制定有效的管控策略.

1 经典传染病模型

不同的情绪强度不仅对个体产生影响,其本身可以通过疏散个体之间的情绪感染而蔓延,而情绪传染类似于传染病的传播[10].传统的传染病学模型主要有两类:SIR模型和SIS模型.

SIR模型是指个体状态为“易感者”、“感染者”、“免疫者”三种类型.S,I,R分别表示三类群体占总群体的比例.SIR模型在传染过程中不计算新出生率和死亡率,故S+I+R=1.该模型反映了人群传染的宏观过程:当传染疾病发生时,三类人群总数量恒定,“易感者”一旦与“感染者”接触,以β概率被传染;“感染者”以一定的概率γ转变到“免疫者”;“免疫者”等同于疫苗注射者,对疾病免疫，不再被传染.

另一种经典SIS模型考虑某些疾病传染时个体很难达到完全免疫,因此在传染过程中仅存在“易感者”和“感染者”.该模型的微分系统动力学方程如式(1)所示,状态转移如图1所示.

图1 SIS模型状态转移示意图

(1)

式中:β为日接触传染率;γ为日治愈率.

由于SIS模型采用常微分方程表示传染动力学宏观模型,是一个宏观、静态数学模型,并不适合用于描述疏散个体的微观情绪感染和自由移动个体状态;并且SIS模型假定人群总数在传染过程中不发生变化,而疏散时行人一旦离开出口即为安全,在危险区的人数每时每刻都是变动的.综上,经典SIS并不完全适用.本文结合CA构建突发事件下考虑情绪感染的行人疏散动态模型.

2 情绪传染机制

2.1 状态定义

当室内突发异常事件,如火灾、有毒气体扩散、非法侵害等,行人均带有紧张情绪.因室内空间有限,个体对危机认知及心理承受能力的差异,具有不同情绪状态.模型中疏散个体根据其紧张情绪值分为潜在情绪状态和显著情绪状态.

2.2 情绪感染机制

疏散个体情绪值与自身持有情绪及周围他人情绪相关[9],定义PEj(t)表示t时刻个体j的紧张情绪强度值,取值范围(0,1],即在疏散全过程中个体不存在紧张情绪消失的时刻.个体j的紧张情绪强度值随着时间变化的动态方程如式(2)所示.

PEj(t)=PEj(t-1)+(1-ξ)·Δej(t)+τ.

(2)

其中:Δej(t)表示t时刻个体j情绪感知增量;ξ是感知衰减系数,取值范围[0,1],体现了情绪感知非理性;τ为随机扰动.

此处,情绪感知增量Δej(t)是t时刻个体j的速度引起的情绪变化量Δfj(v,t)及周围行人的情绪均值Δφj(t)的相关函数,如式(3)所示.

Δej(t)=Δfj(v,t)+Δφj(t) .

(3)

该函数是情绪强度值的动态演化的关键,是本文的研究重点,下面对其构成开展研究.

1) 个体速度变化引起的情绪增量Δfj(v,t).在疏散过程中,若个体行走速度低于最大期望速度,他的不良情绪(烦躁、紧张、焦虑)增长较快.那么,t时刻个体j由速度引起的情绪感知增量函数如式(4)所示[12]:

(4)

该函数可以清楚表示移动速度的变化对个体紧张情绪的波动.若实际速度远小于最大期望速度,Δfj(v,t)值将趋近于1,即个体原地不动时,紧张情绪接近峰值.

2) 情绪感知域内他人引起的情绪增量Δφj(t).行人疏散对周围他人情绪有认知且受其影响,文献[4,10]定义R为情绪感知域半径,行人j的情绪感知范围离散化后是以该行人元胞为中心,2R为边长的正方形.然而,文献[12]指出,疏散个体视线内的信息传播与非视线内信息传播影响力存在明显差异.据此,本文认为情绪感知域分为可视区域RV和不可视区域RI,且可视区域的情绪传染性高于不可视区域.

t时刻个体j的情绪感知域定义为:t时刻个体j所处元胞坐标cellj(x,y)为中心,半径为R的圆,进行元胞离散化后区域为2R边长的正方形.

t时刻个体j的不可视感知域RI定义为:以个体j所处元胞以坐标cellj(x,y)为中心,2R为边长的正方形中除去RV的部分,如图2深色区域所示.

图2 t时刻个体j的情绪感知域

个体j的情绪感知域内,情绪增量是与距离有关的函数,那么,t时刻个体元胞j感知到的他人情绪变化量函数为式(5)所示:

(5)

其中,距离权重ddk为式(6)所示:

(6)

3 疏散模型

基于情绪感染的疏散模型分为情绪传染、情绪更新和行为决策三个步骤[9],本文构建的动态情绪驱动的行人运动疏散模型以元胞自动机为基础的行人运动模型(crowd emotional cellular automata, CECA),行人运动时因情绪状态不同,采取不同行动策略,构成情绪驱动行为机制.该疏散模型的基本流程是:t时刻疏散个体j更新情绪强度值,情绪强度达到一定阈值后发生状态的更新,进而改变行为策略[3，10],其过程如图3所示.

图3 t时刻疏散个体j一次迭代过程

3.1 情绪更新

t时刻疏散个体j疏散策略迭代过程由情绪更新判别状态更新,执行状态对应的疏散策略.t时刻疏散个体j紧张情绪强度值依据t-1时刻的紧张情绪强度值、行动策略、状态、前进方向有关,受到情绪感知范围内其他人紧张情绪强度值影响,情绪变动将引起疏散个体j状态及行为策略的动态变化,具体状态更新流程如图4所示.

图4 t时刻疏散个体j紧张情绪更新机制

3.2 状态更新

Fu等[10]提出基于情绪阈值的状态改变实现疏散人群状态的动态变化.若潜在情绪状态(S)的紧张情绪强度高于阈值λemotion时,以β概率转成显著情绪状态 (I),β称感染系数.若显著情绪

状态(I)的紧张情绪强度低于阈值λemotion时,以γ概率转为潜在情绪状态(S),γ称平静系数.

3.3 行为策略执行

CA行为策略是个体元胞j下一步的移动策略,个体j情绪状态决定了j的期望速度和期望目标元胞位置,评估期望目标元胞的转移概率,确定最终移动的目标元胞,实施移动策略.

根据个体情绪状态定义不同的行为策略,即:

(7)

其中:ks为到出口距离影响系数；kρ为密度影响系数;cpji(t)表示t时刻疏散个体j的邻域元胞i的状态参数,若被占据,则该值为0,否则为1;coji(t)表示t时刻疏散个体j的邻域元胞i是否为障碍物,若是障碍物,则该值为0，否则为1;ρji(t)表示t时刻疏散个体j的邻域元胞i的周围的密度.Sji(t)表示t时刻个体元胞j的邻域元胞i到出口的最短距离.设个体j的位置坐标为C(xj,yj),邻域元胞i的位置坐标为C(xi,yi),疏散空间中有m个有效出口,每个出口的中点坐标为C(xek,yek),则Sji(t)如式(8)所示.

(8)

4 结果与讨论

将动态情绪驱动行人运动模型与经典CA模型进行比较,验证该模型是否符合疏散时个体运动特征.在此基础上,分析情绪强度阈值λemotion,情绪感染系数β,情绪平静系数γ对总疏散时间的影响.

4.1 模型有效性分析

模拟的场景取文献[12]的基础场景,单位元胞为cl×cl,cl=0.4 m,18cl×14cl的单出口房间,出口宽度为el1=2cl.初始状态T=0时,150位疏散个体随机分布在空间位置上,每人占据一个元胞.疏散时间以时间步Time step为单位，用T表示，每一时间步定义为0.5 s,疏散个体移动时采用Moore邻居规则.CECA模型的其他参数设置为:λemotion=0.6,ξ=0,β=0.5,γ=0.5.图5和图6分别记录两个模型在相同仿真场景下不同时间步的密度图.

由图5和图6可以看出,随着时间推移CECA模型与CA模型中的人流特性相同,两种模型的行人均逐渐聚集在出口位置,在出口处表现出行人流拱形特征,但由于CA模型整体疏散较慢,其拱形特征更加明显.在相同时间步,CECA模型中个体移动速度更快,可以表明情绪波动及感染对减少总体疏散时间是有利的,这结论与文献[8]一致.

图5 经典CA模型不同时刻人员密度图

图6 CECA模型不同时刻人员密度图

为了验证CECA模型的动态特性,在同一模拟场景下,改变情绪强度阈值λemotion,得到每一时间步的系统平均速度,如图7所示.由图可见,CECA模型中情绪动态变化,驱动了行为策略的动态变化,可以反映实际疏散过程中速度的动态性[1].而由于经典CA模型中的每一个行人元胞遵循同样的规则进行同步更新,行人状态单一,不能体现不同行人的差异化的运动特征.

图7 系统平均速度

4.2 情绪强度阈值 λemotion的影响

文献[10]将情绪强度阈值设定为λemotion=0.6开展了仿真实验,然而并未讨论阈值对疏散效果的影响.为了分析情绪强度阈值对总疏散时间的影响,文中模拟一个12 m×12 m的单出口房间,出口宽度为el1=3cl,设参数ξ=0,β=0.5,γ=0.5.初始人数分别为N=100,200,300,400,500时,个体随机分布在空间中,每次仿真程序运行100次得到平均总疏散时间,结果如图8所示.

图8 不同情绪阈值和不同初始人数时的系统总疏散时间

由图8所示,情绪强度阈值对系统总疏散时间是有影响的,其影响可以分为三个区间:λemotion∈[0,0.4]的疏散时间最快;λemotion∈[0.6,1]时的系统疏散时间最慢;而λemotion∈[0.4,0.6]时，是总疏散时间发生显著变化的区间.

进一步分析情绪阈值影响疏散时间的机理,对两类状态群体疏散过程动态演化进行了仿真实验.在上述同一场景中以200人为例,相同的初始位置分布和初始情绪值,情绪阈值分别取λemotion=0.3,0.5,0.7,两类状态群体在疏散场地内的位置如图9所示,蓝色表示S状态个体,红色表示I状态个体,绿色表示出口.

情绪阈值λemotion=0.3时,个体分布如图9a和图 9b 所示, 初始时系统内I的人数为 155,S的人数为45,I状态的人数多于S;T=5时系统中绝大

图9 不同情绪阈值初始时刻和T=5时刻群体分布

多数是I状态个体,I状态的人数为185,S状态的人数为2.由图9e和图9f可见,当λemotion=0.7时,初始时系统内S状态的比例较大,此时S状态的人数为156,I状态的人数为44;而T=5时系统内绝大多数是S状态个体,S状态和I状态的人数分别为183和6.λemotion=0.5时,个体分布如图9c和图9d所示,初始时,S状态的人数为101,I状态的人数为99,两种状态人数相当,T=5时系统中S状态的人数为100,I状态的人数为87,短时间内并未明显出现哪类群体多于其他群体.

4.3 感染系数β和平静系数γ的影响

感染系数β和平静系数γ是疏散过程中个体情绪驱动行为改变的重要因子,主要分析其对总疏散时间的影响.以上文场景作为模拟场地,初始人数为200,考虑情绪阈值分别为λemotion=0.5和λemotion=0.6两种情况,分别探讨感染系数β和平静系数γ单独变化时系统中显著情绪状态(I)的比例,仿真程序每种情况执行100次取平均值,结果如图10所示.

由图10可见，γ=0.1时,λemotion=0.5,感染系数β值越大,系统中I的比例由0.5增加到1的速度越快.而λemotion=0.6时,无论β为何值时,系统中I群体的数量逐渐减少.β=0.1时,λemotion=0.5,系统中I群体比例最终都趋于1.I群体比例经过一个小幅下降后,逐步上升,特别是在T=40之后振荡上升至1.而λemotion=0.6时,系统最终均为S群体,平静系数γ数值越大,系统变为全S群体的速度越快.