交通网络中疏散路线设计

2015高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：重庆师范大学参赛队员(打印并签名)：1.王晨晨2.赵越3.彭穗军指导教师或指导教师组负责人(打印并签名)：日期：2015年7月29日赛区评阅编号（由赛区组委会评阅前进行编号）：2015高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：交通网络中疏散路线设计与调度方案摘要在发生对人们的健康和安全造成严重危害的自然或人为的紧急事件情况下，大规模的疏散和避难所避难是保护人口免受潜在危害的主要选择。疏散部署就是指把有限的救援力量投入到最需要救援的地方，使效率最高。针对问题1:首先利用GM(1，1)灰色模型法合理评估疏散区域的人口规模并转化为用于疏散的车辆数，引入符号ijN表示第i处疏散处到第j处避难处的最优车辆，ijC表示第i处疏散处到第j处避难处的运输次数，ijD表示第i处疏散处到第j处避难处的距离数。根据目标函数：VDzijmin与2ij1ijijminVDVDCz，利用Lingo软件，求得最优疏散时间为。针对问题2:首先简化模型，用Matlab计算出所有避难处和疏散处的坐标和最短距离，分析建筑物应急疏散空阔网络中任意节点的待疏散人员完成安全疏散的最优路线，并进行模拟，其结果表现为全局最优目标的实现。然后根据Warshall算法完成计算各疏散处到避难处的最短路径，把距离疏散处距离最短的避难处作为最佳匹配的避难处，构建了以最短路径为目标函数的整数规划模型，最后再考虑道路阻塞的情况下，分多次输送，得到最短疏散时间和整体疏散方案。由于重庆师范大学和重庆大学A、B、C区的建筑群较多，我们把大学看作一个特殊的具有一定的辐射范围的特殊避难处处理，然后进行数据处理与求解。随后我们结合实际路况并用调查得到的住宅的实际相关数据对模型进行验证，模型的疏散时间和疏散路程误差为：道路阻塞度为：结果证明了模型是科学、合理可操作的。针对问题3:我们采取0-1整体规划模型对疏散人员的行为偏好进行假设，得出更加符合实际情况的应急疏散复杂系统，利用0-1整数规划模型对此进行数据处理和求解，根据网络流原理和最优化理论，对疏散人群的行为进行有效假设，进行疏散性能的动态分析和疏散行动决策的全局优化。得到改进后更加贴近实际的模型。其性能指标主要包括任意时刻t预期能够实现安全疏散的人数，完成安全疏散所需的时间及最佳疏散路线的全局寻优等。本文中定义了两个评价原则：原则一：将某疏散处所有人员运送到避难处所需时间≦10min；原则二：保证道路阻塞密度不超过负荷峰值，且尽量接近于最优值；然后依据问题分析中两个评价原则，对所得方案性能进行评价。关键词：最优分配灰色人口预测道路阻塞模型Ford算法多目标规划一问题重述如飓风、洪水、火灾或化学泄漏等自然或人为的紧急事件，会对人们的健康和安全造成严重危害。在这种情况下，大规模的疏散和避难所避难是保护人口免受潜在危害的主要选择。请你以沙坪坝地区为例对交通网络中的疏散路线进行设计。首先，在疏散计划过程中，疏散处和避难处应该是确定的。从实际观点来看，住宅小区和公司商厦被确定为疏散处（见图1中红点）；运动场馆、学校、医院、广场和大学建筑等有容纳大数目人群的特点，应该被确定为避难处（见图中1绿点），其中重庆师范大学和重庆大学A、B、C区的建筑群较多，不能简单确定为一个避难处。其次，评估疏散区域的人口规模是必要的，这个人口规模应该转化为用于疏散的车辆数，考虑道路上的车辆速度会随着车辆数增多而下降，并制定出疏散方案。接着，当你确定了疏散路线方案后应该验证其性能，可以从疏散时间、疏散路程和阻塞规律等方面进行考虑。最后，对疏散人群的行为进行有效假设用来改进你的模型，这会使你的疏散计划更贴近实际，例如疏散人员离开疏散处的时间是一个泊松分布或疏散人员对疏散路径选择有一定的偏好等等。针对以上要求，我们要研究的问题如下：（1）合理评估疏散区域的人口规模及用于疏散的车辆数，制定出疏散方案；（2）从疏散时间，疏散路程，道路阻塞等方面，验证疏散路线方案的性能；（3）对疏散人群的行为进行有效假设来使模型更贴近实际；二问题分析针对问题1：首先利用GM(1，1)灰色模型法合理评估疏散区域的人口规模及转化为用于疏散的车辆数。灰色模型法不直接利用原始数据，而是通过累加生成灰色模型，滤去原始数据中可能混入随机量，从上下波动的时间寻找某种隐含规律，然后利用Warshall—Ford算法求最短路径，此时认为应急疏散系统中关于人的疏散行为的数学模型为假设所有的待疏散人员具有相同的疏散能力，井井有条地按预先制定的疏散计划，完成疏散行动。主要的任务是研究任意时刻，群集的疏散进展。最后制定出合理的疏散方案，计算出交通密度与用于疏散的车辆之间的关系。针对问题2：简化模型，用Matlab计算出所有避难处和疏散处的坐标和最短距离，分析建筑物应急疏散空阔网络中，任意节点的待疏散人员完成安全疏散的最优路线，并进行静态和动态模拟模拟结果表现为任意时刻，不同事故状态下，各节点完成安全疏散的全局最优目标的实现。最后完成最短矩阵距离与最佳速度之间的求解，考虑道路阻塞，分多次输送，得到最短疏散时间。其中疏散空间网络由节点和通道组成，其中各节点和通道均具有多个处于动态变化中的属性特征，如完成疏散的时间、疏散距离等，称为疏散成本属性。各节点和通道上的各种疏散成本属性值存在很大的差异，且随事故状态的发展而变化，表现为时间的函数可以将任一的应急疏散空间模化为G(U，E)网络。其中节点集u=|u1，u2......uN|节点可分为三类：源节点u1,传输节点u2和出口目标节点u3。源节点即只有流出群集，无流入群集的节点。出口目标节点u3属吸收式，即只存在由节点u2或节点u1指向u3的单向的群集流动。针对问题3：对疏散人群的行为进行有效假设，得到改进后更加贴近实际的模型。利用网络流原理和最优化理论，根据应急疏散空间事故状态的即时变化，考虑疏散人员的随机疏散行为特点，进行疏散性能的动态分析和疏散行动决策的全局优化。本文中定义了两个评价原则：原则一：将某疏散处所有人员运送到避难处所需时间≦10min；原则二：保证道路阻塞密度不超过负荷峰值，且尽量接近于最优值；然后依据问题分析中两个评价原则，对所得方案性能进行评价。三模型假设（1）假设避难处和疏散处所在位置固定不变。（2）用于疏散的车辆型号相同，即每辆所载人数和速度都相同。（3）假设用于疏散的车辆足够多（4）任何一个疏散处的人优先到达最近的避难所。（5）不考虑处用于疏散车之外的车辆造成的道路状况。（6）避难处与避难处，疏散处与疏散处之间相互独立，没有影响。（7）只考虑用于疏散的车辆走主干道路和次干道路。（8）假设每个疏散车的速度都是相同的，即疏散车到避难所的时间只疏散距离和道路阻塞程度有关。四符号说明Q：交通量（辆/h）；V：速度（km／h）；K：交通密度（辆／km）；ijN：第i处疏散处到第j处避难处的最优车辆数；fV：畅行速度；jK：阻塞密度；ijD：第i处疏散处到第j处避难处的距离；in：第i处疏散处所需车辆；iS：第i处疏散处总人数；x：每辆车的最大载人数；ijC：第i处疏散处到第j处避难处的运输次数；),(yxjj：第j个避难处的坐标；),(qpii：第i个疏散处坐标。五：模型的建立与求解5.1建模准备住宅小区和公司商厦被确定为疏散处（见图中红点）；运动场馆、学校、医院、广场和大学建筑等有容纳大数目人群的特点，应该被确定为避难处,为了找到疏散处到避难处的最短路径，先利用以沙坪坝地区交通网络中的疏散路线图作出疏散处到避难处的分布图，利用Matlab绘制出其分布图的各点坐标，然后计算出各个疏散处到避难处的距离，最终确定最短路径下的车辆调度方案。考虑重庆师范大学和重庆大学A、B、C区的建筑群较多，不能简单确定为一个避难处后，将所给图像预处理后得到的如图所示：图1：沙坪坝地区避难处与疏散处分布图5.2模型的建立5.2.1疏散区域的人口评估模型利用GM(1，1)灰色模型法，从灰色系统的建模，关联度以及残差辨识的思想出发，第一步设有原始数列：}...,{)0()()0()2()0()1()0(nxxxx，X对)0(x做一次累加，生成数列：}...,{)1()()1()2()1()1()1(nxxxx，式中ijjiXx1)0()()1()(，i=1,2,3...n。对数据使用这种累加生成技术后，使其变为较有规律的生成数列，然后再建立微分方程模型，因此灰色模型实际上是生成数列模型：uaxdxdx)1()1(。第二步求参运算，应用最小二乘法解a和u：NTTYBBBua1)(，其中，12/))()1((12/))3()2((12/))2()1(()1()1()1()1()1()1(nxnxxxxxB，那么微分方程的解，即时间响应函数为：aueauxkxxak)()1()1()1()1(，其中a,u为待估价的参数。需要说明的是此时预测出来的值是一个累加值，当需预测该值所在地点的数据时，要用这个值减去前一个预测值，即作累加还原，可得原始数据的估计值：)()1()1()1()1()0(kxkxkx5.2.2道路长度评估模型假设各节点之间是有连线的，在有连线的道路上根据题目中的已知条件和两点之间的长度计算公式，可以得到该城市内任意相邻两个避难处和疏散处的距离：22(,)()()iiiiliixxyy，综合计算，可求得任意两交叉路口间距离为：...22...22j()()(,)()()ijiiiiijijjjjxxyyijdijxxyyij当时；当时。对于i~：iiK~K，即(,)(,)iiiiKxyKxy；对于j~：jjK~K，即(,)(,)jjjjKxyKxy。其中，,iiKxy表示第i个交叉路口坐标；(,)iiKxy表示第i个交叉路口的邻接路口集合，坐标表示。可以得到各个路径之间的距离，再用Matlab编程使之在地图的道路上标出各自的距离。5.2.3交通密度评估模型对疏散区域的人口规模进行评估之后，将人口规模转化为用于疏散所需的车辆数，即in=[xiS]+1691i考虑道路上的车辆速度会随着车辆数增多而下降，根据交通量、速度和交通密度的关系式Q=V·KK=ijijDN691i，261j由格林希尔兹速度-密度线性模型V=fV(1-jKK)导出Q=fV（K-j2KK）=JK(V-f2VV)是二次函数关系，可用一条抛物线表示，故当流量最大时，共有两个约束条件，即：V=21fV，K=21jK当道路流量最大时可得最优车辆数，即ijN=ijj2DK691i，261j考虑疏散处所需车辆数与最优车辆数不相等，则有两种情况情况一：当inijN时，只需运输一次，并能以最大速度运输，可建立下面的约束条件：K=ijinD691i，261j情况二：当inijN时，需分批次运输，次数为：ijC+1=[ijinN]+1691i，261j前ijC次以最优车辆数运输，道路流量达到最大，速度为21fV，最后一次运输与情况一类似。由ij