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

辽宁工业大学2012年数学建模（论文）题目：火灾紧急撤离问题院（系）：机械工程及自动化专业班级：机械1106班学生姓名：王哲、郭爽、吴建彬起止时间：2012.5.21—2012.5.27摘要本文借用流体动力学中的微分关系，通过将离散的人员转化为连续的人流，以人流密度为研究主体，建立了人员撤离的动态微分方程优化模型，分析了地震发生时人员紧急撤离的问题。并根据我们所在教学楼的楼层建筑的数据分别估算了混乱状况下与有组织时人员撤离的时间，为人员的紧急撤离提供了参考方案。第一，本文分析了在无组织的状态下，人员撤离的一般情形。一方面，无组织下人员的运动具有随机性，故此引入人流密度作为基本研究对象。另一方面，流量的变化率是人流密度对距离积分后对时间的导数，人流量对时间的积分即为撤离人员的数量。由此几方面关系，可以列出整个动态过程的微分方程。经分析发现，单位时间的人流量与密度和速度成正比关系，而整体的人流速度与密度之间又是成一次线性关系，恰好符合流体力学中的流量、流速与密度之间的关系。根据实际情况对整求解过程做了简化，以楼道中的平均人流量为研究主体，最终以数值解求得全部人员逃离所需时间大约为420s.第二，利用得出的人流量随时间变化的图像可知，由于人员无组织的涌出教室，导致人流密度很大，人群得不到有效的移动，从而使流量达到最大值后又迅速减小。故最好的撤离方式是在达到流量最大的时候，保持住一定的人流密度从而来维持最大的流量。结合数据后可知，在撤离开始一分钟的时候应该有人组织撤离，这样可以避免由于人员的过多涌入楼道而导致的拥堵现象。这样子调控后最佳的撤离时间可以降到240秒左右。第三，除去人为堵塞的因素对撤离时间影响较大外，改变楼层的设计同样可以缩短撤离所用时间。于是，文章讨论了实际楼层中的参数，如楼层中疏散通道的宽度、教室门的宽度以及疏散口的数量等，对紧急撤离时间的影响。并得出结论疏散口的增加与疏散通道的加宽对撤离时间的缩短有明显的提高。最后，由于不同的楼层人员速度不一样会导致在楼道中的互相推挤现象，此举对人员在楼道中人员的有效流动有较大影响。故我们引入混乱时间的概念，用来具体量化由此导致的时间的浪费情况。分析后可知混乱时间主要决定于相临两层人员的速度差，由于混乱时间与速度差成正比关系，而且在速度差为正值的时候时间较大，而为负值时时间较小，故利用指数函数来表示两者的关系。由此建立了以总的混乱时间最小为目标的优化模型。利用atlab对各种指派情形进行比较，得出最了优解。关键词：人流量动态微分方程最佳撤离混乱时间目录一、问题的提出.................................................................................1二、基本假设及符号说明......................................................................12.1基本假设....................................................................................12.2符号说明....................................................................................1三、问题分析.........................................................................................23.1问题一.........................................................................................23.2问题二.........................................................................................23.3问题三.........................................................................................33.4问题四........................................................................................3四、模型准备.........................................................................................3五、模型建立.........................................................................................4六、参考文献.......................................................................................11七、附录...............................................................................................121一、问题的提出近年来，诸如地震、火灾等突发事件时常发生。虽然人们在很多情况下还不能准确预报这些突发事件，但当灾难发生时，尽可能在灾难中减少伤亡人数是人们应对突发事件的首选。在突发事件中，身处灾难环境的人员快速撤离灾难地点可以有效减少伤亡人数。本着居安思危的态度，假设某一天上午，学生正在我校C教学楼上课，突然该楼发生火灾，请你完成任务:1.用数学建模方法，给出一种使学生快速撤离C教学楼的方案；2.针对我校C教学楼，用你的方案给出其第一、二楼学生快速撤离的具体方案和所用的时间。二、基本假设及符号说明2.1基本假设疏散过程中，人群的流量与疏散通道的宽度、行走速度有关；所有人员在突发事件发生后同时疏散，中途不退后；所有人员在疏散过程中不发生踩踏事件；每个年级在同一个楼层；2.2符号说明符号说明ρi(x,t)第x位置，第t时刻，第i层楼的人员密度；Qi(x,t)第x位置，第t时刻，第i层楼的人流量；Lw走廊宽度；2Ls楼梯宽度；б(ρ)人员密度为时的拥挤调控系数；Vw撤离人员在走廊的平均速度；VwMax撤离人员在走廊的最大速度；Vs撤离人员在楼梯的平均速度；VsMax撤离人员在楼梯的最大速度；R(t)第i层楼进入楼梯间的人数；M单位时间内从教室进入走廊的人员数；Ni第i个教室的总人数；三、问题分析3.1问题一根据人流运动的特点，建立基于流体动力学的微分方程模型。将每一个楼层分为教室出口处和非教室出口处，由于不同位置流入流出的人流量不同，故可以动态的分析出不同时刻不同位置的人流量密度。以任意小区间段的人流量为考虑对象，该区间两端人流量之差即为该区间人数变化率。根据此关系建立微分方程，接触各个位置的人流量密度函数。此时，用出口处的人流量对时间积分即为已撤离出的人员数目。求解积分方程，得到当撤离人数为教学楼总人数时所经历的时间，即为人全部撤离出所消耗的时间。3.2问题二通过分析问题一结果发现影响撤离时间的主要因素为人流量密度，当人流量3密度过大时会导致人员移动区间变小，使得撤离速度的下降。故当撤离人流量最大时，应使撤离人数与从教室流入走廊的人数相同，以此保持人流量密度一直保持在最好的水平，使得单位时间内撤离的人数最多，此为最佳撤离方案。3.3问题三分析问题一所列出的微分方程式，结合实际，认为走廊、楼道宽度以及楼梯数量为制约撤离时间的主要因素，应适当改进。并结合实际情况中的其他因素，提出若干可行性建议。3.4问题四认为每一个楼层只安排一个年级，将撤离时间定为无干扰时的撤离时间以及相邻两层互相制约而产生的混乱时间之和。无论楼层如何安排，无干扰撤离时间不变。而对于混乱时间，其数值正比于下层撤离速度与该层速度之差。如果下一层的速度比本层大，则混乱时间很小，如果下一层的撤离速度大于该层，则会产生很大的混乱时间。基于此情况，定义混乱时间是关于相邻两层人员撤离速度差的指数函数。通过编程，对各种情况进行遍历，可以求出混乱时间最小时的楼层安排方案，即为最合理的教室分配方案。四、模型准备4.1基本公式准备撤离人员从走廊、楼梯撤离，其情形就像在湍急的江河中奔腾的流水一样。故我们运用流体动力学中的概念去解决该问题。根据流体动力学中速度与密度的定义，得到速度与流体密度关系公式：其中：ρm表示单位平面内能够容纳的最多人数；VwMax表示人在撤离稳定时期能够行走的最大速度；此公式表示在人流密度增大初期，人行走速度逐渐降低，当人流密度到达最大值的时候由于没有行走空间，所以人行走速度降为零。由于在稳定的时候，撤离人员处于稳定连续状态。所以根据流体动力学知识有如下方程：4表示：时刻t，区间[a，b]内的撤离人员的数量为单位时间内通过a，b点的流量Qi(a,t),Qi(b,t)之差等于撤离人员数量的变化率。4.2楼梯间中行走长度的确定根据实际情况，假设在楼梯间行走的距离为，如图1：图楼梯长度示意图五、模型建立5.1针对问题一5.1.1模型一将整栋楼看作一个整体，设其密度均匀，人员从20个教室流入，从楼出口处流出，根据3.1节中基本公式列微分方程求解。将整个撤离过程分为稳定前阶段和稳定阶段两个时间段。1）稳定前阶段设定稳定前状态即为距离楼梯间最近的教室的人员到达下一楼层之前的时间段。此时，由于整个教学楼均处于畅通的状态，所以设此阶段的撤离速度为Vbefor。由此可以确定出到达稳定阶段时，楼梯间以及走廊的初始密度。5可以得到稳定前阶段所消耗的时间：2）稳定阶段根据3.1节中的公式，可以列微分方程：其中：表示整个教学楼的平均密度；表示在楼梯间中的平均撤离速度；表示走廊长度总和；表示楼梯间的长度总和；该方程未考虑由于走廊拥堵而造成的教室中人员无法到达走廊的情况，为考虑此情况，引入拥挤调控系数表示在教室出口处的密度与能够从教室出去人数的调控函数，如果出口处密度为零时,即教室里面的撤离人员均能到达走廊。如果出口处的密度已为最大值，即教室里面的撤离人员无法到达走廊。此时的微分方程为：根据3.1节，将速度与密度的关系代入上式有：解微分方程可以得到人员撤离密度关于时间的函数，设经过T时间，所有人员撤离完毕，有：6反解出时间T即为稳定时期撤离完毕所花费的时间。故全部人员撤离消耗的时间为：由于此模型过于简化，并没有考虑到各层楼走廊和楼梯间的拥堵情况，分析问题过于粗略。故我们将每层的走廊和楼梯间分开考虑，引入模型二。5.1.2模型二将每一个楼层看作一个整体，将其密度视为均匀，记为ρi(t)。此时依旧从稳定前与稳定两个阶段去分析整个撤离过程。稳定前阶段与模型一一致，在此不详细论述。对于稳定时期，对每个楼层进行分析，根据模型一的思想，列方程如下：代入人流量与密度方程，化简得到：如模型一考虑，引入拥挤调控系数б(ρ)，得到方程如下：解微分方程可以得到，得到第i个楼层进入楼梯间的人数：每一个楼层均为此，考虑楼梯间的情况。将每一个楼层视为教室，整个楼梯间为走廊，列类似微分方程：解该微分方程，可以得到楼梯间的人流密度ρ(t)。此时，得到最终从出口撤离的人员数量与时间的关系：5.1.3模型三对整个走廊和楼梯间的整个长度分成无限小段，基于流体动力学知识，建立更具有一般性的模型，计算出全部人员撤离完毕需要时间。依旧将整个撤离过程分为稳定前、稳定两个阶段进行分析。1）针对稳定前阶段设定稳定前状态即为距离楼梯间最近的教师的人员到达下一楼层之前的时间段。此时，由于整个教学楼均处于畅通的状态，所以设此阶段的撤离速度为Vbefor。7由此可以确定出到达稳定阶段时，楼梯间以及各层走廊的初始密度。可以得到稳定前阶段所消耗的时间：