地震紧急撤离问题数学建模

辽宁工业大学2012年数学建模（论文）题目： 火灾紧急撤离问题院（系）：机械工程及自动化专业班级：机械口06班学生姓名：王哲、郭爽、吴建彬起止时间：二摘要本文借用流体动力学中的微分关系，通过将离散的人员转化为连续的人流，以人流密度为研究主体，建立了人员撤离的动态微分方程优化模型，分析了地震发生时人员紧急撤离的问题。并根据我们所在教学楼的楼层建筑的数据分别估算了混乱状况下与有组织时人员撤离的时间，为人员的紧急撤离提供了参考方案。第一，本文分析了在无组织的状态下，人员撤离的一般悄形。一方面，无组织下人员的运动具有随机性，故此引入人流密度作为基本研究对象。另一方面，流量的变化率是人流密度对距离积分后对时间的导数，人流量对时间的积分即为撤离人员的数量。山此儿方面关系，可以列出整个动态过程的微分方程。经分析发现，单位时间的人流量与密度和速度成正比关系，而整体的人流速度与密度之间乂是成一次线性关系，恰好符合流体力学中的流量、流速与密度之间的关系。根据实际悄况对整求解过程做了简化，以楼道中的平均人流量为研究主体，最终以数值解求得全部人员逃离所需时间大约为420s.笫二，利用得出的人流量随时间变化的图像可知，山于人员无组织的涌出教室，导致人流密度很大，人群得不到有效的移动，从而使流量达到最大值后乂迅速减小。故最好的撤离方式是在达到流量最大的时候，保持住一定的人流密度从而来维持最大的流量。结合数据后可知，在撤离开始一分钟的时候应该有人组织撤离，这样可以避免山于人员的过多涌入楼道而导致的拥堵现象。这样子调控后最佳的撤离时间可以降到240秒左右。第三，除去人为堵塞的因素对撤离时间影响较大外，改变楼层的设讣同样可以缩短撤离所用时间。于是，文章讨论了实际楼层中的参数，如楼层中疏散通道的宽度、教室门的宽度以及疏散口的数量等，对紧急撤离时间的影响。并得出结论疏散口的增加与疏散通道的加宽对撤离时间的缩短有明显的提高。最后，III于不同的楼层人员速度不一样会导致在楼道中的互相推挤现象,此举对人员在楼道中人员的有效流动有较大影响。故我们引入混乱时间的概念，用来具体量化山此导致的时间的浪费情况。分析后可知混乱时间主要决定于相临两层人员的速度差，山于混乱时间与速度差成正比关系，而且在速度差为正值的时候时间较大，而为负值时时间较小，故利用指数函数来表示两者的关系。山此建立了以总的混乱时间最小为LI标的优化模型。利用atlab对各种指派情形进行比较，得出最了优解。关键词：人流量动态微分方程最佳撤离混乱时间目录TOC\o"1-5"\h\z\o"CurrentDocument"一、 问题的提出 1\o"CurrentDocument"二、 基木假设及符号说明 2基木假设 2符号说明 2\o"CurrentDocument"三、 问题分析 3问题一 3问题二 3问题三 3问题四 3\o"CurrentDocument"四、 模型准备 4\o"CurrentDocument"五、 模型建立 5\o"CurrentDocument"六、 参考文献 121七、 附录 13一、问题的提出近年来，诸如地震、火灾等突发事件时常发生。虽然人们在很多情况下还不能准确预报这些突发事件，但当灾难发生时，尽可能在灾难中减少伤亡人数是人们应对突发事件的首选。在突发事件中，身处灾难环境的人员快速撤离灾难地点可以有效减少伤亡人数。木着居安思危的态度，假设某一天上午，学生正在我校C教学楼上课,突然该楼发生火灾，请你完成任务：1•用数学建模方法，给出一种使学生快速撤离C教学楼的方案；2.针对我校C教学楼，用你的方案给出其第一、二楼学生快速撤离的具体方案和所用的时间。二、基本假设及符号说明基本假设疏散过程中，人群的流量与疏散通道的宽度、行走速度有关；所有人员在突发事件发生后同时疏散，中途不退后；所有人员在疏散过程中不发生踩踏事件；每个年级在同一个楼层；符号说明符号 说明Pi(x,t) 第x位置，第t时刻，第i层楼的人员密度;Qi(x,t) 第x位置，第t时刻，第i层楼的人流量;Lw 走廊宽度；Ls6(P)VwVwMaxVsVsMaxR(t)MNi楼梯宽度；人员密度为时的拥挤调控系数；撤离人员在走廊的平均速度；撤离人员在走廊的最大速度；撤离人员在楼梯的平均速度；撤离人员在楼梯的最大速度；第i层楼进入楼梯间的人数；单位时间内从教室进入走廊的人员数;第i个教室的总人数；三、问题分析问题一根据人流运动的特点，建立基于流体动力学的微分方程模型。将每一个楼层分为教室出口处和非教室出口处，由于不同位置流入流出的人流量不同，故可以动态的分析出不同时刻不同位置的人流量密度。以任意小区间段的人流量为考虑对象，该区间两端人流量之差即为该区间人数变化率。根据此关系建立微分方程，接触各个位置的人流量密度函数。此时，用出口处的人流量对时间积分即为已撤离出的人员数LI。求解积分方程，得到当撤离人数为教学楼总人数时所经历的时间，即为人全部撤离出所消耗的时间。问题二通过分析问题一结果发现影响撤离时间的主要因素为人流量密度，当人流量密度过大时会导致人员移动区间变小，使得撤离速度的下降。故当撤离人流量最大时，应使撤离人数与从教室流入走廊的人数相同，以此保持人流量密度一直保持在最好的水平，使得单位时间内撤离的人数最多，此为最佳撤离方案。问题三分析问题一所列出的微分方程式，结合实际，认为走廊、楼道宽度以及楼梯数量为制约撤离时间的主要因素，应适当改进。并结合实际情况中的其他因素，提出若干可行性建议。问题四认为每一个楼层只安排一个年级，将撤离时间定为无干扰时的撤离时间以及相邻两层互相制约而产生的混乱时间之和。无论楼层如何安排，无干扰撤离时间不变。而对于混乱时间，其数值正比于下层撤离速度与该层速度之差。如果下一层的速度比本层大，则混乱时间很小，如果下一层的撤离速度大于该层，则会产生很大的混乱时间。基于此情况，定义混乱时间是关于相邻两层人员撤离速度差的指数函数。通过编程，对各种情况进行遍历，可以求出混乱时间最小时的楼层安排方案，即为最合理的教室分配方案。四、模型准备基本公式准备撤离人员从走廊、楼梯撤离，其情形就像在湍急的江河中奔腾的流水一样。故我们运用流体动力学中的概念去解决该问题。根据流体动力学中速度与密度的定义，得到速度与流体密度关系公式：pm其中：Pm表示单位平面内能够容纳的最多人数；VwMax表示人在撤离稳定时期能够行走的最大速度；此公式表示在人流密度增大初期，人行走速度逐渐降低，当人流密度到达最大值的时候由于没有行走空间，所以人行走速度降为零。由于在稳定的时候，撤离人员处于稳定连续状态。所以根据流体动力学知识有如下方程：Qi(g,/)=—[|/?j(x,z)xQi(g,/)=—[|/?j(x,z)xLw\cbi表示：时刻t,区间［a,b］内的撤离人员的数量为［“(卸)&单位时间内通过a,b点的流量Qi(a,t),Qi(b,t)之差等于撤离人员数量的变化率。楼梯间中行走长度的确定根据实际悄况，假设在楼梯间行走的距离为2 ,如图1：2 ,如图1：图楼梯长度示意图五、模型建立针对问题一模型一将整栋楼看作一个整体，设其密度均匀，人员从20个教室流入，从楼出口处流出，根据节中基本公式列微分方程求解。将整个撤离过程分为稳定前阶段和稳定阶段两个时间段。1)稳定前阶段设定稳定前状态即为距离楼梯间最近的教室的人员到达下一楼层之前的时间段。此时，山于整个教学楼均处于畅通的状态，所以设此阶段的撤离速度为Vbeforo山此可以确定出到达稳定阶段时，楼梯间以及走廊的初始密度。— 壬(厶wxL\)IbeforeQ找就何二Q找就何二A/X可以得到稳定前阶段所消耗的时间:品H品H+L\Ibefore2)稳定阶段根据节中的公式，可以列微分方程：— —d—20jW+p(/)xLsxK=—|p(/)x(L»-x5^+LsxX9)|dt其中:3(/)表示整个教学楼的平均密度；v■-表示在楼梯间中的平均撤离速度；弘表示走廊长度总和；表示楼梯间的长度总和；该方程未考虑山于走廊拥堵而造成的教室中人员无法到达走廊的情况，为考虑此情况，引入拥挤调控系数6(。)表示在教室出口处的密度与能够从教室出去人数的调控函数，如果出口处密度为零时，即教室里面的撤离人员均能到达走廊。如果出口处的密度已为最大值，即教室里面的撤离人员无法到达走廊。此时的微分方程为：20jV/x5(p)+p(f)xZsxI/=—|p(/)x(Lj^x5w+AsxX?)|dt根据节，将速度与密度的关系代入上式有：20A7xc(p)-t-x?(r)xL.$x (I-—)=—[p(z)x(L^xSV+LsxSs)]Pmdt解微分方程可以得到人员撤离密度关于时间的函数，设经过T时间，所有人员撤离完毕，有：⑴乂心衣厶*卩一也|=ya；Jo rtr反解出时间T即为稳定时期撤离完毕所花费的时间。故全部人员撤离消耗的时间为：Time-兀如+T由于此模型过于简化，并没有考虑到各层楼走廊和楼梯间的拥堵情况，分析问题过于粗略。故我们将每层的走廊和楼梯间分开考虑，引入模型二。模型二将每一个楼层看作一个整体，将其密度视为均匀，记为Pi(t)o此时依旧从稳定前与稳定两个阶段去分析整个撤离过程。稳定前阶段与模型一一致，在此不详细论述。对于稳定时期，对每个楼层进行分析，根据模型一的思想，列方程如下：4M—0(4/?+厶1+L2+厶3+厶■»)=G(f)x(4£?+L\+L2+£s4-Z.4)xLw代入人流量与密度方程，化简得到：4AY—pi(t)X刃亿如x[l-P⑴]XLW=/?;(/)X(4D+厶1+厶2+厶3+La)XLwQ/w如模型一考虑，引入拥挤调控系数6(p),得到方程如下：4XIx6(pi)—/?/(/)xMaxx11———xLwx橈)=—Q(f)x(4Z?+Zj+Li+厶3+厶4)xLwr)m (it解微分方程可以得到，得到第i个楼层进入楼梯间的人数:Ri(t)=pxJ'irwxX(1-—)xZjfxQ(q如)每一个楼层均为此，考虑楼梯间的情况。将每一个楼层视为教室，整个楼梯间为走廊，列类似微分方程：&⑴X6(p)-p(t)X x[l——XLs-=—/Xt)X x厶pw dt解该微分方程，可以得到楼梯间的人流密度P(t)o此时，得到最终从出口撤离的人员数量与时间的关系：/?喽⑴=0⑴x厶xKsa^x[1-Pm模型三对整个走廊和楼梯间的整个长度分成无限小段，基于流体动力学知识，建立更具有一般性的模型，计算出全部人员撤离完毕需要时间。依旧将整个撤离过程分为稳定前、稳定两个阶段进行分析。针对稳定前阶段设定稳定前状态即为距离楼梯间最近的教师的人员到达下一楼层之前的时间段。此时，山于整个教学楼均处于畅通的状态，所以设此阶段的撤离速度为Vbeforo山此可以确定出到达稳定阶段时，楼梯间以及各层走廊的初始密度。L\Vbejor^L\Vbejor^m(LwxAi)Ibefore可以得到稳定前阶段所消耗的时间：丁冋+L1Sian二 ['before2）针对稳定阶段对每一层进行考虑，根据实际位置分布建立坐标系，如图图任意楼层走廊坐标系的建立确定每一楼层每一个时刻不同位置的撤离人员密度。由于在教室的出口处会有撤离人员流向走廊，所以在每个教室的出口处，密度变化率大于其他位置，如图： ►■ ► D D图教室出口处的密度变化示意图i教室出口处的密度变化率：取平均位置，设撤离人员流出位置为教室门处，列微分方程如下:=| 创X」W・X+ -&/)x丄沪xjk(x- +"X△/-a(.x」)x厶沪xl+(x)xA/= a(為r曲丰［/Xx-Ax,r)xLwx14(x-Ar)4M一p(x,/)>(丄护x%(x)1xAr」r简化得:

—Ax.t)xLwxVw(x-Ax)—Ax.t)xLwxVw(x-Ax)+M—p(jrj)xLwxPV(x)令AtfO有：原式二lima* A/二q(jv-心J)x£眇xViv二lima* A/二q(jv-心J)x£眇xViv(j-Aj)+A/-pi(xj)x厶妒x)iv楼梯间里面各个位置密度的变化率：针对于楼梯间内部，将其类比为某一楼层的走廊，各个楼层进入楼梯间的撤离人员数量相当于该楼层各个教室进入走廊的数量。基于此假设，列出类似微分方程：dry-与—\ /）dy y-勒j）*厶％心（》一纫）+忆⑴朋（Q”）-以”/）xisx取刃atJy丿为交接口处•sdr尸为— Q（川旳=刊y-为；“X厶乂临y）-p{j,r）xLsxVs{y}atJyb为剩余位置最终只需计算出从出口撤离的人员数量与时间的关系即可：R楼⑴二Q樓//,r)xLsxPs.VtoX[l-pR楼⑴二Q樓//,r)xLsxPs.VtoX[l-——二 1Qiw此时，经过时间T,所有人员撤离完毕，有如下方程:TJ尺楼（？）力=工入1/=!反解出T即可。此时所有人员撤离需要的时间为:11加£=1start十J但是山于模型三方程过于复杂，无法给出解析解，只能利用微分转差分进行数值分析。鉴于前两个模型已经给出相应的撤离时间，故此模型并没有作相应的计算，仅作为一更详细的基本方程。问题二模型建立题LI要求给出撤离时间最短的时候的人员撤离方案，山图可以清楚的看出，撤离人员的总人数即为图像所围成的面积。在总人数不变的情况下，要缩短撤离总时间，只能使得各个时刻出口处的人流量保持最大值，增大图像面积，如图8竖线所示。根据上题分析，出口处人流量降低的主要原因是楼道人流量密度过大，导

致撤离速度降低、人流量变小。为解决上述问题，应当在撤离过程进行了50秒左右时，保持楼梯间的人员数量，使得人流量密度为最恰当的数值。由图8可知，人流量最大为每秒人，要保持人流量密度则要从教室进入走廊的人数等于撤离教学楼的人数。模型检验根据参考文献中，给出的相应的疏散时间经验公式:Trine=Trine=丁曲“+其中:TmartNLsE表示稳定前时间；表示总人数；表示逃生口宽度；表示有效流动系数；根据经验取E=人/秒，根据我们的基本数据，计算出来的撤离总时间为秒。与计算得到的撤离时间相似，结果较为符合实际。问题三结合问题一中模型三的方程，考虑楼房设计的最合理方案，给出该微分方程:

5石 Q楼（叩心）R楼（r）二°建（ Vs.\iaxx[l ]2 pm根据上述方程，可以明显看出，要是某时刻撤离的人员最多能采取如下方法:1） 增大走廊的宽度以及楼梯间的宽度；2） 一般应在幕近建筑标准层的两端设置楼梯间。楼梯间的布置应均匀。且有多方向的出口。3） 高楼层应配套增加紧急通道，可以直接从5楼到达1楼；4） 增设广播设备，当险悄发生时能够及时通知，缩短人员的反应时间；问题四模型准备根据节中的用作模型检验的撤离时间经验公式，定义撤离时间如下：LsxQLsxQVSMax其中：T 为需要的总时间；Tmart为达到稳定状态之前需要的准备时间；N 为整栋楼的总人数；Q 平均每秒钟的撤离的人数；考虑拥堵时候的撤离时间如果不考虑拥堵情况，第i层人员的撤离时间满足上述定义，总时间即为其求和。现在考虑拥堵情况，故在此基础上加入混乱时间表示笫i层的人员和第i+1层的人员交汇之后因为拥堵产生的附加时间。时整个撤离时间即为:7-=fL+fAL1=1 1=1i认为与相邻两层人员的速度差成正比关系，||如果Vi>Vi+l,说明下层人的速度大于