1、基于排队论的校园自行车排队问题建模与仿真研究摘要:随着校园内自行车数目不断递增,校园内人流较多的地方的交通拥堵问题变得越来越严重。本文利用排队论的相关理论,建立自行车排队系统模型对拥堵区域进行了相关理论分析。同时利用模拟仿真软件 ProModel 对校园局部自行车流动情况进行模拟并作出相关优化方案。这样理论结合软件模拟的系统优化研究方法对于实际的交通系统优化具有一定的参考价值。 关键词:排队论;交通流;ProModel 仿真;系统优化; Research on Modeling and Simulation of campus bicycle problems based on queuing
2、 theory Yang Fan (College of civil engineeringSoutheast UniversityNanjing211189) Abstract: As the number of the school bicycles increasing, traffic congestion problems are becoming more and more serious. This paper uses the queuing theory to build the model to analyze the congestion area. Also, the
3、paper uses the ProModel software to build the model of the traffic of the campus and work out the optimum proposal. This method used in this paper by calculating and modeling to get the optimum proposal has the reference value for actual traffic system optimization. Keywords: queuing theory, traffic
4、 flow, ProModel simulation, system optimization 中图分类号:C29 文献标识码:A 0.引言 由于高校扩招,导致各个高校的新生数目不断增加,而学生在校园内的主要交通工具自行车的数目也在逐年递增,在食堂、教学楼等区域的高峰时期往往会存在较为严重的局部交通拥堵现象。部分道路上的排队现象既影响了通行的效率,也存在一定的安全隐患,因此利用排队论的理论结合模拟软件仿真的相关内容对其进行优化研究具有非常重要的实际意义。 1.排队论基本概念 排队论(queuing theory), 或称随机服务系统理论, 是通过对服务对象到来及服务时间的统计研究,得出这些数量
5、指标(等待时间、排队长度、忙期长短等)的统计规律,然后根据这些规律来改进服务系统的结构或重新组织被服务对象,使得服务系统既能满足服务对象的需要,又能使机构的费用最经济或某些指标最优。 排队论是为了解决是为了解决排队问题而发展起来的学科。它建立在系统排队模型的基础上,分析系统的内在规律,使系统的性能得到优化。一般地,排队论模型如下图所示: 由上图可知,排队系统的运行过程是:顾客按照一定的规则进入系统;如果服务台忙,则顾客按一定的排队规则进入排队队列,等待服务;当服务台空闲时,将按照一定的服务规则从队列中选择顾客,并为顾客提供服务,接受服务的顾客在服务完成后离开系统。1 排队论通过研究人员在各个阶
6、段的时间和规则的概率分布情况来实现对整个排队系统的研究,通过对相关参数的调整和控制从而达到优化系统的目标。 2.校园自行车交通系统建模 根据排队论的基本理论,针对校园自行车的系统,我们也可以利用相关的理论来研究,如下图: 自行车通过一定的到达规则进入系统,如果自行车太多,则自行车需要排队来到达车位进行停车,当交通没有发生堵塞,或者轮到排队的自行车时,那么自行车将会按照一定的停车规则来停车,自行车停放时间符合相应的停车规律,到达时间后自行车便会离开到达下一个停车点。由以上的一些初步的认识我们可以得到,建立自行车排队论系统,需要的元素至少有:自行车(bicycle,可能有不同的种类) 、到达规则(
7、rules for arrival) 、排队规则(rules for queuing) 、停车规则(rules for parking) 、排队队列(queue) 、车位(places for parking) 、停车时间(parking time)等,而这些元素需要进行统计观测以及数据分析得出。 3.模型分析 校园自行车排队系统主要是针对校园自行车车流较大的局部区域,可以通过前期的考察对实际的地形进行简化并划分排队区域和停车区域。排队现象因为在排队区域中车流过大而使得自行车无法及时到达停车区域。其中自行车的到达一般符合泊松分布,停车时间一般符合负指数分布。因此校园自行车排队系统可以看成是 M
8、/M/1(N)排队系统。 以多通道排队 M/M/N(单通道排队系统可以取 N=1)来分析,系统假设自行车平均到达率为 ,平均行车时间为 ,记 =/N 为交通强度或利用系数,当 1 时,系统处于不稳定状态,排队会越来越长;1 时,系统处于稳定状态,排队队伍可以消散。因此,要想使得交通顺畅的条件是 1,即 N。2 设随机变量 其概率密度函数为: (x)=e-x(x)n-1n-1!x0 其中 n 为正整数 E=0xe-x(x)n-1n-1!dx 积分后得: E=n,由于 E 就是队列的平均停车时间 1/N,所以, Nn=。 系统中有车辆的概率 P0=1k=0N-1kk!+NN!(1-/N)2 系统中
9、有 k 辆车 Pk=kk!P0(kN) kN!Nk-N(kN) 系统中平均排队长度 Lq=N+1N!NP0(1-/N)2 系统中平均车辆数n=Lq+ 自行车在系统中平均消耗时间d=Lq+1 平均等待时间=Lq3 4.实例分析 以东南大学九龙湖校区教学楼楼群之间的区域进行分析,教学楼平面简图如下: 图 1 教学楼平面图 从图中可以看出,八个教学楼四周都有相应的自行车停车区域,将整个教学楼围合内部作为分析系统。实际情况显示在高峰阶段,95%以上的自行车都是经过道路 1 进入系统,而通过道路 3 进入系统的自行车非常少,可以忽略。因此,在道路 2 的阴影部分十分容易形成排队,而这个排队区域也是研究的
10、主要对象。 通过实地观测统计道路 2 入口处车流情况,在上课期间、午休期间以及晚上 10 点以后自行车车流忽略不计(由于自行车车流太小,不会造成排队,故不作为研究对象) ,在上下课的高峰时段自行车单向流量约为 3300 辆/h,考虑各项因素,正常情况下每辆车步行时速度为 1.5m/s,骑行时速度为 8m/s。 可以计算出 =33003600=0.92 辆/s, 由于形成排队时自行车绝大多数都是步行,且由于行人影响,故取自行车步行速度 1.0m/s,自行车长取 2m/辆, =1.0/2=0.5 辆/s, 对于单通道排队系统,=1.841,所以系统不稳定,排队长度会越来越长,严重后果阻碍交通。 对
11、于 N2 的排队系统,N1,故对于多通道排队系统系统都是稳定的,可以列表来分析各项指标,见表 1: 表 1 自行车排队系统排队指标 指标 系统 M/M/2 M/M/3 M/M/4 10.58 0.68 0.022 12.42 2.52 1.862 13.5 2.73 2.024 11.5 0.73 0.024 0.043 0.161 0.183 0.08 0.296 0.337 0.074 0.272 0.31 0.068 0.167 0.17 0.062 0.102 0.057 0.002 0.053 0.049 0.044 P(8) 0.47 从表 1 中可知,车道的数目对于车流量有积极的
12、影响,而车道的数目反映到校园内道路上就是道路的宽度,从上述计算可知,当道路宽度在单方向上容许 3 辆以上自行车同时走过时排队现象会大幅减轻,因此学校在平衡经济和地理环境因素后可以适当拓宽道路宽度使得排队系统更加优化。 5.软件模拟 相比于上述排队系统的简化计算,实际情况显然要复杂得多。首先自行车的到达和停车并不一定服从泊松分布和负指数分布,真实的概率分布情况需要根据实际的统计得出,不仅如此,仍选用教学楼为例,教学楼处排队情况并不是仅在入口处存在排队,而是在每栋教学楼的入口处都可能有排队发生。如果综合考虑这些实际情况现在通过计算来建模太过复杂,因此可以利用仿真模拟软件来模拟实际的排队系统的状态。
13、 针对校园自行车排队系统模拟仿真采用的模拟软件是ProModel,ProModel 是一款离散事件仿真软件,它可以构造多种生产、物流和服务系统模型,而对于本系统可以采用简化、等效等方法模拟出教学楼围合处自行车的排队情况,最终对于系统进行优化调整。 5.1 区域划分及数据处理 将教学楼围合处停车车位按照实际情况进行划分为七大停车区域: 据统计:区 221 个车位,区 156 个车位,区 210 个车位,区 176 个车位,区 173 个车位,区 562 个车位,区 200 个车位 统计主要测量上午 7:408:00,9:359:45,下午 13:4014:00,15:3515:45,时间段自行车
14、的流动情况,因为这些时间段为自行车流动的可能高峰期,容易体现存在的关键问题,最后将所有时间整合在一起作为一套连续的自行车高峰时段排队情况模拟。 图 2 教学楼停车区域划分图 通过统计及数据拟合,得出: 自行车到达入口的概率分布为 L(0.0159,4.15,1.66)min; 自行车在教一教二之间行动的分布概率为 N(25.7,11.1)sec; 自行车到达停车点所用的锁车时间分布概率为 T(6,12,54.1)sec。 5.2 基本假设 由于教学楼的特殊性质,自行车在教学楼的停留时间绝大多数是45 分钟的倍数,而下课可能有前后 5 分钟的容错时间,因此以 45n 表示停留时间,通过调查可知:
15、 n=1,1%; n=2,35.4%; n=3,13.5%; n=4,46.7%; n=5,3.4% (此比例根据几个较大院系学生课表及各班人数确定,教学楼内部移动不计入自行车移动); 假定教一和教二之间,教三和教四之间的马路最多同时容许 30 辆自行车通过,教一、教二门口的道路允许同时通过 15 辆自行车; 假定当教一教二之间马路上车辆数目小于 35 时,后来的自行车有15%的比例停放在教一教二南部的停车点,否则比例为 25%; 假定自行车从入口走到教一教二的马路上用时 3 秒; 假定自行车从教一教二之间的马路走到其他三个马路上用时 7 秒,在其他三个马路上平均用时 15 秒; 假定从入口到教一教二南部停车点用时 10 秒,从道路到停车点用时 10 秒; 非停在教一教二南部的自行车在各教学楼分区内的停车概率如下:教一:22.7%; 教二:20% 教三:22.1% 教四:14.6%