客船人员疏散仿真模型的验证方法.pdf

第 卷第期 年 月武汉理工大学学报(交通科学与工程版)J o u r n a l o fWu h a nU n i v e r s i t yo fT e c h n o l o g y(T r a n s p o r t a t i o nS c i e n c e&E n g i n e e r i n g)V o l N o O c t 客船人员疏散仿真模型的验证方法李亚朋)蔡薇)胡敏)陈祺)(江苏科技大学海洋装备研究院)镇江 )(武汉理工大学绿色智能江海直达船舶与邮轮游艇研究中心)武汉 )(中船邮轮科技发展有限公司)上海 )摘要:文中探讨了两类可对仿真模型进行验证的疏散场景,分别为疏散分析指南中的 个测试场景(T e s tT e s t)和 实 船 疏 散 试 验 项 目S A F E GUA R D中 的 两 个 实 船 场 景(S GV D S 和S G V D S)利用软件P a t h f i n d e r对各场景进行仿真,并说明完整的验证过程关键词:客船;疏散仿真;疏散场景;模型验证;P a t h f i n d e r中图法分类号:U d o i:/j i s s n 收稿日期:第一作者:李亚朋(),男,博士,主要研究领域为船舶人员疏散安全通信作者:蔡薇(),女,博士,教授,主要研究领域为船舶智能设计方法、游艇邮轮设计基金项目:工信部高技术船舶项目(G C Z )、工信部高技术船舶项目(工信部装函 )引言国际 海 事 组 织(i n t e r n a t i o n a lm a r i t i m eo r g a n i z a t i o n,I MO)下属机构海事安全委员会(m a r i t i m es a f e t yc o mm i t t e e,M S C)在 年发布第一版客船疏散分析指南M S C/C i r c ,用于辅助海上人命安全公约(s a f e t yo f l i f ea ts e a,S O L A S)对脱险通道的评估,而后经过多次修改,M S C在 年通过 修订的新造和现有客船疏散分析指南M S C/C i r c (后文简写为 指南),于 年强制执行 指南 中给出了两种疏散分析的方法,简化方法和高级方法,前者基于宏观的水力网络模型,后者基于微观的疏散仿真模型简化方法较容易使用,但当船舶布局形式变得复杂时,其评估结果的可信度会大大降低;高级方法需要建立船舶疏散空间的完整几何模型,赋予每个疏散个体不同的行为特征(如行走速度、反应时间等),因此能反映更真实的疏散过程,而且该方法对疏散过程和拥堵点的可视化呈现效果也更好在理论研究和实际工程中,一般更青睐运用仿真方法评估客船的疏散性能 评估过程中,可能需要构建一种新的仿真模型然后加以运用,或者直接使用一种现有的仿真软件然而,许多研究在使用仿真工具时并未分析其在船舶领域的适用性,未经验证便将其应用于指南中的疏散仿真场景,这将对疏散性能的计算和评估结果产生直接的影响文中探讨了两类可以用于验证疏散仿真工具的场景,包括 指南 中对仿真工具提出的 个测试场景(T e s tT e s t)和实船疏散试验项目S A F E GUA R D 中的两个实船场景(S GV D S 和S GV D S),详细展现完整的验证过程,并对本文的研究进行讨论和总结 I MO疏散 指南 中的测试场景 指南 对疏散仿真工具提出了四个方面的验证需求)模块测试(c o m p o n e n tt e s t i n g)用于检验工具的各个组成部分是否按照预期运行,该部分包括七个测试T e s t T e s t)功能测试(f u n c t i o n a lv e r i f i c a t i o n)用于检验模型是否能实现其各项功能或任务,任务不同测试则不同,所以未规定具体测试,但要求测试过程应包含在该工具的技术文档中)定性测试(q u a l i t a t i v ev e r i f i c a t i o n)用于检验模型是否能够反映人在疏散中的典型行为或现象,该部分包括五个测试T e s t T e s t )定量测试(q u a n t i t a t i v ev e r i f i c a t i o n),用于检验模型的预测结果与真实疏散演习的结果是否吻合,由于缺乏足够的试验数据,未规定具体的测试图为P a t h f i n d e r对 个测试场景的仿真过程截图图P a t h f i n d e r对 个测试场景的仿真T e s t 在走廊上维持设定的行走速度一个人以m/s的行走速度通过长 m、宽m的走道,则其所花的时间应是 s T e s t 在上楼梯时维持设定的行走速度一个人以m/s的行走速度通过长 m、宽m的向上斜坡或楼梯,则其所花的时间应是 s T e s t 在下楼梯时维持设定的行走速度一个人以m/s的行走速度通过长 m、宽m的向下斜坡或楼梯,则其所花的时间应是 s T e s t 出口的流率不能超过 人/(ms)个人从mm的房间里逃生,出口m宽,位于m边的中间,则整个过程中出口的流率应不超过 人/(ms)(这里出口宽度是m,所以也可是 人/s)T e s t 在设定的反应时间之后才开始移动 个人从mm的房间里逃生,出口m宽,位于m边的中间,这 个人的反应时间服从在 s,s 上的均匀分布,则每个人应经过适当的反应时间才开始移动T e s t 在走过拐角时身体不穿过边界 个人初始均匀分布在mm的矩形中,他们需要经过一个拐角,则每个人的身体边界都不应与拐角的边界有重合(即不会出现穿墙现象)T e s t 各年龄段人群的速度参数与 指南中一致将 指南 中男性、岁年龄段的速度参数赋予给 个人,即在 m/s,m/s 上的均匀分布则每个人应按照适当的行走速度移动T e s t 对向流中的人越多,疏散时间越长在这个场景中,两个边长 m的正方形房间由一条长 m、宽m的通道连接,左边房间放置 个人,他们向右边房间移动,同时右边房间放置N人,他们向左边房间移动,这样两边的人群将逆向流动 N分别为、,分别记录次实验中左边房间 人全部到达右边房间的时间,则随着逆流的增加(右边房间人数的增加),所记录的时间应增加T e s t 总出口宽度越短,疏散时间越长 个人均匀分布在长 m、宽 m的大房间中,每个长边上都有两个宽m的出口,人群选择离自己最近的出口逃出房间做两次实验,一次个出口全开,一次只开个出口,记录清空房间所需的时间,则由于总出口宽度减半,第二次的时间应大概是第一次的倍T e s t 按照分配的路径(或出口)逃生构建图T e s t 的舱室区域,其中包括 个舱室,个出口,个人分布于左侧八个舱室,令其向较宽的主要出口移动;八个人分布于右侧四个舱室,令其向较窄的次要出口移动观察这 个人的行走路径,则他们应按照分配的出口逃生T e s t 稳定的人群流通过楼梯时在楼梯底部会出现拥堵在这个场景中,一间长m、宽m的房间与一个宽m的向上楼梯由一条长 m、宽m的走道相连,人均匀分布于房间中,他们需要依次走出房间、通过走道、走上楼梯观察人群在各个位置的密度,则可以发现在房间的出口处将产生拥堵,人群流在走道上逐渐趋于稳定然后在楼梯底部再次产生拥堵T e s t 高密度的人群流通常比中等密度的人群流更慢 指南 中对该测试未设置具体的测试场景本文构建了一个长 m、宽m的走道,N人均匀分布于走道m,m 的区域内,人群向右移动,区域 m,m 作为密度(p/m)和流率(p/s)的测量区域分别设置N为、进行仿真实验,记录测量密度和流率,则当人群密度非常高时(N ,密度大于人/m),流率比中等密度时(N ,密度约为武汉理工大学学报(交通科学与工程版)年第 卷人/m)的流率小P a t h f i n d e r的模拟结果均符合预期,因此其通过了该项验证(在P a t h f i n d e r官网上,可查阅到该软件的开发者对 指南 中的T e s tT e s t 做了更加详细的测试)S A F E GUA R D项目中的实船疏散数据集 指南 中提到了定量测试,即与实船疏散演习的数据进行比较,但由于数据缺乏,未强制要求进行该测试欧盟第七框架计划(E UF P)中的实船疏散试验项目S A F E GAUR D()公开了两条船的试验数据集(S A F E GAUR D V a l i d a t i o nD a t aS e t s,S GV D S 和S GV D S),数据是根据在船上进行的集合演习得到,可考虑作为船舶疏散仿真工具的补充验证各个数据集提供了构建疏散场景及运行和验证仿真工具所需的所有信息,包括船舶几何(g e o m e t r y)、人口参数(p o p u l a t i o n)、反应时间分布(r e s p o n s et i m ed i s t r i b u t i o n)、到达曲 线(a r r i v a lc u r v e s,反映疏散时间与到达集合站人数之间关系的曲线)、衡量指标(v a l i d a t i o nm e t r i c)、接受准则(a c c e p t a n c ec r i t e r i a)等验证的基本流程为:根据A u t o C A Dd x f图纸文件构建船舶的舱室与通道;根据人群的初始位置将其分布于船舶上的对应位置,设置合适的目标点、反应时间、年龄段等参数;将该场景仿真 次,每次需要改变一次人群的部分参数;得到 次仿真的整体到达曲线,根据衡量指标找出与试验数据吻合度最高的一次仿真;对该仿真的整体到达曲线和各集合站到达曲线进一步计算衡量指标,检查各指标数值是否通过接受准则S GV D S主要使用到达曲线上的数据点(即某个到达人数对应的时间点)计算衡量指标,指标共有四个,分别为:)欧式相对差(t h ee u c l i d e a nr e l a t i v ed i f f e r e n c e,E R D)用于评估模型数据和试验数据之间的距离,理论上可能取值的范围是,值为表示模型数据与试验数据完全一致,计算公式见式()欧氏投影系数(t h ee u c l i d e a np r o j e c t i o nc o e f f i c i e n t,E P C)用于评估模型数据和试验数据之间的吻合程度,理论上可能取值的范围是(,值为表示吻合程度最好(但不能说明两者完全一致),计算公式见式()正割余弦(t h es e c a n tc o s i n e,S C)用于评估模型和试验的到达曲线形状的相似程度,该指标利用了各曲线的一阶导数,理论上可能取值的范围是(,值为表示两者曲线形状完全一致,计算公式见式()总集合时间百分差异(p e r c e n t a g ed i f f e r e n c eo f t o t a la s s e m b l yT i m e,d i f fT AT)用于衡量总集合时间的差异,与前面三个指标不同,该指标的计算仅需用到一个数据点(最后一个人到达的时间点),计算公式见式()EmE ni(Eimi)niEi()E,mmniEiminimi()E,mEmnis(EiEis)(mimis)s(titi)nis(EiEis)s(titi)nis(mimis)s(titi)()T(TETm)TE ()式中:Ei为试验的第i个数据点;mi为模型的第i个数据点;n为数据点总数;s为光顺项,用于移除试验和模型数据的噪点,取值根据数据点总数确定;TE为试验的总集合时间;Tm为模型的总集合时间 数据集S G V D S 该数据集源自一艘客滚渡轮,船上总乘客数 人,总疏散人数 人,计入到达曲线中的人数 人(有些人起始和终止位置在同一个集合站,有些人没有被标记但仍然参与了演习,绘制到达曲线时需排除这些人),人群分布于酒吧、餐厅、零售店、座位区、集合站等区域,涉及三层甲板,总共四个集合站该场景建模三维图见图从 次仿真中选择整体到达曲线E R D值最小的一次仿真,对其整体和各集合站到达曲线进行进一步评估各指标的接受准则为:E R D ;E P C;当s/n 时,S C;d i f fT A T 评估过程分两步:()对整体到 达 曲 线 计 算E R D,E P C,S C和d i f f第期李亚朋,等:客船人员疏散仿真模型的验证方法图S G V D S 客滚船疏散建模三维图T A T,当所有指标均满足接受准则时,才能进行接下来的评估;()分别对四个集合站的到达曲线计算E R D、E P C和S C,当 个值中的至少个值满足接受准则时,评估通过P a t h f i n d e r对S GV D S 的验证结果见表可以看到其通过了第一步和第二步的评估,因此该模型通过了S GV D S 的验证 P a t h f i n d e r仿真的和S GV D S 试验的整体到达曲线见图 P a t h f i n d e r只是勉强通过了评估(第二步中,刚好 个值中的个满足条件),而且仿真与试验的曲线差别较大,其中的原因可能是:实际的演习中,大量未被标记的人也参与了集合过程,影响到计入到达曲线的人的集合,从而影响试验曲线;另外,各集合站的到达曲线的数据点较少,都只有不到 个表P a t h f i n d e r对S G V D S 的验证结果nE R D E P C S C(s/n )d i f fT A T整体 集合站A 集合站B 集合站C 集合站D 注:n数据点个数;各指标未标灰色背景通过接受准则;标灰未通过图P a t h f i n d e r仿真曲线与S G V D S 试验曲线(整体)的比较 数据集S G V D S S GV D S 中的船舶是一艘大型邮轮,船上总乘客数 人,总疏散人数 人,计入到达曲线中的人数 人(需要排除本来就在集合站没有移动的人),人群分布于居住舱室、餐厅、休闲区、集合站等区域,涉及 层甲板,总共个集合站,该场景建模三维图见图图S G V D S 邮轮疏散建模三维图相比S GV D S,该数据集中各指标的接受准则更为严苛,为:E R D ;E P C;当s/n 时,S C;d i f fT A T 由于该数据集中的集合站也是个,所以评估过程与S GV D S 完全一样 P a t h f i n d e r对S GV D S 的验证结果见表由表可知:其通过了第一步和第二步的评估,因此该模型通过了S GV D S 的验证 P a t h f i n d e r仿真的和S GV D S 试验的整体到达曲线见图由表可知:P a t h f i n d e r仿真得到的结果满足了所有指标的接受准则,而且仿真的和试验的到达曲线吻合度非常高表P a t h f i n d e r对S G V D S 的验证结果nE R D E P C S C(s/n )d i f fT A T整体 集合站A 集合站B 集合站C 集合站D 图P a t h f i n d e r仿真曲线与S G V D S 试验曲线(整体)的比较结论)I MO仅要求进行其 指南 给出的验证,因此,即使没有验证S A F E GUA R D数据集也能进行 指南 规定的疏散分析/评估;但后一项验证中的数据集是笔者目前已知的最为完整的且规模最大的实船演习数据集,这些数据集也得到了其他学者的认可 ,可以考虑作为船舶人员疏散仿真模型的补充验证场景)从 指南 中的验证场景来看,I MO并未要求船舶仿真工具具备不同于陆上仿真工具的功武汉理工大学学报(交通科学与工程版)年第 卷能或特性,因此不是一定需要用m a r i t i m e E X O DU S、E V I或A E N E A S等船舶专用软件才能进行 指南 规定的疏散分析,只要仿真工具通过适当的验证即可)S A F E GUA R D数据集给出了详细的疏散场景和验证过程,能够较好地支持新的船舶人员疏散仿真模型的验证然而,由于实船试验的复杂性,如许多未被标记的乘客也参与了集合过程,可能使得数据集难以反映真实的疏散结果参 考 文 献I MO I n t e r i m g u i d e l i n e sf o ras i m p l i f i e de v a c u a t i o na n a l y s i so f r o r op a s s e n g e r s h i p s:M S C/C i r c SN e wY o r k:M a r i t i m eS a f e t yC o mm i t t e e,I MOR e v i s e dg u i d e l i n e so ne v a c u a t i o na n a l y s e sf o rn e wa n de x i s t i n gp a s s e n g e rs h i p s:M S C/C i r c SN e wY o r k:M a r i t i m eS a f e t yC o mm i t t e e,K L P F E LH,ME Y E R KN I G T,WAHL EJ,e ta l M i c r o s c o p i c s i m u l a t i o no f e v a c u a t i o np r o c e s s e so np a s s e n g e rs h i p sJT h e o r y&P r a c t i c a lI s s u e so nC e l l u l a rA u t o m a t a,():HAS,KUN,R OH M,e t a l C e l l b a s e de v a c u a t i o ns i m u l a t i o nc o n s i d e r i n gh u m a nb e h a v i o r i nap a s s e n g e rs h i pJ O c e a nE n g i n e e r i n g,():HU M,C A IW,Z HAO H S i m u l a t i o no fP a s s e n g e rE v a c u a t i o nP r o c e s s i nC r u i s eS h i p sB a s e do nA M u l t i G r i dM o d e lJ S y mm e t r y,():张尧考虑水上环境的客船人员应急疏散过程分析与建模D成都:西南交通大学,孙军基于人员疏散仿真的客船人员疏散问题研究D大连:大连海事大学,GA L E AER,D E E R ES,B R OWN R,e ta l A ne v a c u a t i o nv a l i d a t i o nd a t a s e t f o r l a r g ep a s s e n g e r s h i p sJ P e d e s t r i a na n dE v a c u a t i o nD y n a m i c s,():B R OWNR,GA L E AER,D E E R ES,e t a l P a s s e n g e r r e s p o n s e t i m ed a t a s e t s f o r l a r g ep a s s e n g e r f e r r i e sa n dc r u i s es h i p sd e r i v e df r o ms e at r i a l sJ I n t e r n a t i o n a lJ o u r n a lo f M a r i t i m e E n g i n e e r i n g,():D E E R ESJ,G A L E A E R,F I L I P P I D I SL,e ta l D a t ac o l l e c t i o n m e t h o d o l o g i e su s e di nt h eS A F E GUA R Dp r o j e c t t oc o l l e c th u m a nf a c t o r sd a t aCS a f e g u a r dP a s s e n g e rE v a c u a t i o nS e m i n a r,L o n d o n,WANG W L,L I USB,L OSM,e ta l P a s s e n g e rS h i pE v a c u a t i o nS i m u l a t i o na n d V a l i d a t i o nb yE x p e r i m e n t a lD a t a S e t sJP r o c e d i a E n g i n e e r i n g,():P A R K KP,HAM SH,HAS V a l i d a t i o no fa d v a n c e de v a c u a t i o na n a l y s i so np a s s e n g e rs h i p su s i n ge x p e r i m e n t a l s c e n a r i oa n dd a t ao ff u l l s c a l ee v a c u a t i o nJC o m p u t e r si nI n d u s t r y,():V a l i d a t i o nM e t h o do fP a s s e n g e rS h i pE v a c u a t i o nS i m u l a t i o nM o d e lL IY a p e n g)C A IW e i)H U M i n)C H E NQ i)(G r e e n&S m a r tR i v e r s e a g o i n gS h i p,C r u i s e a n dY a c h tR e s e a r c hC e n t e r,W u h a nU n i v e r s i t yo fT e c h n o l o g y,W u h a n ,C h i n a)(C S S CC r u i s eT e c h n o l o g yD e v e l o p m e n tC oL t d,S h a n g h a i ,C h i n a)A b s t r a c t:T w ok i n d so fe v a c u a t i o ns c e n a r i o st h a tc a nv e r i f yt h es i m u l a t i o n m o d e lw e r ed i s c u s s e d,n a m e l y,t e s t s c e n a r i o s i nt h ee v a c u a t i o na n a l y s i sg u i d e(T e s t T e s t )a n dr e a l s h i ps c e n a r i o si nt h er e a l s h i pe v a c u a t i o nt e s tp r o j e c tS A F E GUA R D(S GV D S a n dS GV D S)T h es i m u l a t i o ns o f t w a r eP a t h f i n d e rw a su s e d t os i m u l a t e e a c hs c e n e,a n d t h e c o m p l e t ev e r i f i c a t i o np r o c e s sw a s e x p l a i n e d K e yw o r d s:p a s s e n g e rs h i p;e v a c u a t i o ns i m u l a t i o n;e v a c u a t i o ns c e n a r i o s;m o d e l v a l i d a t i o n;P a t h f i n d e r第期李亚朋,等:客船人员疏散仿真模型的验证方法