CN107330547B - 一种城市公交动态调度优化方法及系统 - Google Patents

Abstract

本发明属于数据处理技术领域，公开了一种城市公交动态调度优化方法及系统，城市公交动态调度优化方法包括：在获取公交车内拥挤度、道路拥挤状况和公交车实时定位等信息的基础上；建立包括时段车次数模型、发车时刻决策模型以及路径变更模型的三层调度决策模型；借助互联网技术实现人、车、路信息交互；时段车次数模型从平衡公交运营成本、乘客等车时间以及乘车舒适性出发，确定最优的时段车次数；发车时刻确定模型以区域总相遇次数最大和总相遇站点数最小为目标调整发车时刻，有效降低换乘等车概率；路径变更模型采用远离拥堵路段的思路，可有效提高行车速度和减少后面站点乘客等车时间，并设置响应站点，供越过路段乘客下车。

Description

一种城市公交动态调度优化方法及系统

技术领域

本发明属于数据处理技术领域，尤其涉及一种城市公交动态调度优化方法及系统。

背景技术

改革公交运营调度模式，来提高公交吸引力成为当前的热点，目前现有的调度技术主要有三种方式，一种就是调度员完全依靠个人经验进行调度，这是最普遍的调度模式，该模式具有调度效率低，反应速度慢以及针对性差等缺点；另一种是调度员在调度中心，在实时获取道路、客流信息的基础上，人工进行调度，目前国内少量应用，该模式同样具有主观性强、调度效率较低以及人工成本高的缺点；最后一种是利用计算机软件依靠历史数据进行制作公交调度计划，该种模式在国内几乎没有，该模式不具有以乘客换乘优先为目的确定发车时刻、实时调整路径，躲避拥堵的优点，针对性不强的缺点也较为突出。

综上所述，现有技术存在的问题是：

现有技术存在的技术问题：1、未考虑公交运营时段的差异，调度结果适应性差；(原因：未采用计算机技术进行调度，增加的工作量，人工难以完成。)2、没有实现定量评估调度方案；(其解决的难度在于定量评估模型的确定以及实时数据的获取)；3、没有最大化利用发车时刻来对公交换乘进行优化(其最大的解决难度在于相关模型的建立与无法保证公交车在途中的行驶时间，本方案的发车时刻确定模型和路径变更决策模型很好的解决了该问题，路径决策模型使公交车的行驶时间较为稳定，受外界的干扰少)；4、无法保证在公交线路出现交通拥堵等意外情况时，及时作出正确决策，保证公交服务水平和公交营运时间。(目前还未有考虑城市道路拥堵情况和因此作出变更路径的调度模式出现)。5、没有考虑设置响应站点，供乘客下车。

本方案相对于以上调度技术，能够解决如下问题：1、对公交营运时间进行了时段划分，实现了精准化调度，提高了调度效果；2、该方案综合考虑了公司、乘客利益，建立量化模型，确定时段车次数，实现了综合利益最大化；3、在时段车次数确定的情况下，调整发车时刻，提升乘客换乘效率。4、建立了变更行驶路径决策模型，极大地降低了交通拥堵对公交运行和服务质量的影响。5、设置响应站点，解决越过站点带来的车内到达越过站点乘客下车的问题。

发明内容

针对现有技术存在的问题，本发明提供了一种城市公交动态调度优化方法及系统。

本发明是这样实现的，一种城市公交动态调度优化方法，所述城市公交动态调度优化方法包括：

通过实时数据获取系统结合红外线和图像识别技术，实时采集乘客到站时刻和数量、公交车内拥挤度、道路拥挤状况和公交车实时定位信息；

建立时段车次数模型、发车时刻决策模型以及变更行驶路径模型的三层调度决策模型；利用信息共享系统实现人、车、路信息交互。

进一步，所述时段车次数模型从平衡公交运营成本、乘客等车时间以及乘车舒适性出发，确定最优的时段车次数；

发车时刻决策模型以区域总的相遇次数最大和总的相遇站点数最小为目标，降低乘客换乘等车概率，确定发车时刻；

变更行驶路径模型以城市道路经常发生交通拥堵和小区开放政策出台为背景，远离拥堵，实现车内乘客、变更路段站点候车乘客及后面站点候车乘客整体最优。

进一步，所述时段车次模型得建立方法包括：

由投资收益率指标以及相应的权重，构建目标函数如式(1)

maxZ＝ω₁×μ₁+ω₂×μ₂+ω₃×μ₃…(1)；

其中，μ₁,μ₂,μ₃分别表示投资收益率、候车时间满意率和拥挤度满意率，具体表达式如下式(2)、(3)、(4)；ω₁,ω₂,ω₃分别为三项指标的权重，由层次分析原理求得；

其中：

自变量为时段车次数正整数n，找到自变量的上下限后，求出最优解；确定忍受的最低满载率r_min和该线路乘客接受的最高满载率r_max；由式(5)得时段车次数n的上下限。

进一步，所述发车时刻优化模型的建立方法包括：

构造换乘站点集

不同的公交线路相互交汇，其交汇点往往就是换乘站点，这些换乘站点组成粗换乘站点集。但在现实生活中，相同两条公交线路可能存在2个及以上的换乘站点，但两线路同一班次公交车不可能在所有站点都相遇，现实情况是基本只能在一个换乘站点相遇。考虑乘客选择换乘站点的特性，当公交车相遇随机发生在各个换乘站点时，乘客选择对换乘站点的选择也呈随机分布；当将公交车相遇控制在同一换乘站点时，乘客也会自发的选择到该换乘站点换乘。如说明书内仅存2个换乘站点的例子，单个换乘站点方案与随机多个换乘站点的方案分别的等车概率为

和

因此，固定换乘站点可有效提高乘客换乘的便利性。本发明在粗换乘站点集的基础上，筛选出在其所在的两条公交线路的换乘集中，所能容纳同时停靠多辆公交车的站点，将该站点设置为该线路的换乘站点，这些换乘站点组成细换乘站点集。

在确定了细换乘站点集之后，先以使在该时段、该区域内两辆公交车在细换乘站点集内的站点总的相遇次数最大为目标，初步确定发车时刻；接着以使相遇点数最小(即使尽可能多的公交车在同一站点、相近时间相遇)调整公交发车时刻，确定最优发车时刻。具体模型如下：

以公交车区域内站点总相遇次数和所有线路总相遇点数作为目标函数，通过公式：

计算每两辆车总的相遇次数最大值；

其中：Z—总的相遇次数；

X_ki—线路k的第i辆车的发车时间；

X_qj—线路q上第j辆车的发车时间；

T_kn—线路k从起点到站点n的时间；

X_ki+T_kn+T_kin—线路k上第i辆车离开站点n的时间；

D_ntki—线路k第i辆车在时间t到达终点n取1，否则取0；

线路上总的相遇点的最小值确定每条线路上的最佳发车时间，可根据公式

计算总相遇点数，Y_nt表示时间t到达站点n的车辆数

其中，M为该区域公交线路的条数；

F_k为该时段线路k的时段车次数；

该模型存在如下的约束条件：

1).时段第一辆公交车发车时间限制：

T_k1≤S_k；

其中，T_k1为线路k在该时段第一辆车的发车时间；

S_k为线路k所规定的该时段第一辆车的最迟发车时间；

2).时段最后一辆公交车发车时间限制：

其中，T为该时段结束时的时间点；

H_k该线路k所允许的最大发车间隔；

为该时段最后一辆车的发车时间；

F_k为该时段的发车车次数；

3).时段内相邻公交发车时间间隔限制：

h_k≤X_k(i+1)-X_ki≤H_k(k＝1,2,…,M,i＝1,2,…,F_k-1)；

其中，h_k表示线路k所允许的最小发车间隔；

X_ki表示线路k在该时段第i车次的发车时刻；

4).站点所能接受的最大相遇车辆数限制：

Y_nt≤R_n；

其中，R_n表示站点n所允许同时停靠的公交车辆数；

5).判断公交车是否在某时刻到达站点：

D_ntki＝max{0,(0.5-|t-(X_ki+T_kn)|)}；

其中：(n＝1,2,……,N；t＝1,2,……,T；k＝1,2,……M；i＝1,2,……,F_k)。

公交车实际到达的时刻与t相比，误差小于0.5分钟，则认为公交车到达。

进一步，路径决策模型的建立方法包括：

该路径变更决策模式是在公交线路上发生的交通拥堵时，进行一定程度的行驶路径变更的一种新型城市公交调度模式。他在保留常规公交低成本性的基础上，提高了公交巴士运行的灵活性。远离拥堵路段可以提高公交巴士的行驶速度和降低后面站点(公交车还未到达的公交站点)乘客等待时间，满足乘客公交出行快捷性和及时性的需求，大大提高公交吸引力。

该模式主要由站台单元、车载单元、路侧单元以及调度控制中心等构成(这4者之间的关系。站台单元、车载单元以及路侧单元分别采集站点乘客到达数、公交车位置以及道路交通状况等信息，并汇总形成数据库，调度控制中心根据信息，做出决策，实现公共巴士合理调度，降低交通拥堵产生的影响，提升公共巴士吸引力。

构造路径集与响应站点：

以是否存在可替代备选路径为标准进行路段划分，交叉口是公交车驶入、驶出规定线路的关键节点。因此，本发明首先将该公交线路上的所有交叉口汇总形成交叉口集；交叉口之间的路径组成粗路径集，接着剔除粗路径集中高于3倍规定路径长度的备选路径，剩余路径组成路径集(路径越长，被采用的概率就越低，限制长度可提高运算效率)。

当公共巴士做出变更行驶路径的决策时，对应的原路径可能没有、可能存在1个或者多个公交车站的情况，没有公交车站的情况不必考虑。当存在公交车站时，公交车内可能有在被越过站点(原路径内的公交车站)下车的乘客，我们采用了设置响应站点的方案，公交车辆在该响应站点停靠，乘客在响应站点下车，乘客可以在响应站点等待下一班车(仅允许变更行驶路径公交车的后一辆公交车在响应站点停靠)。乘客也可以选择步行或自行车出行到达目的站点，因为公交车变更行驶路径的决策是在原路径出现因交通拥堵的情况下做出的，步行等慢性交通与公交车行驶的时间相差不大，甚至会低于公交车行驶时间，响应站点也尽可能选在离被越过站点近的地方。

分析路径变化的影响，建立反应调度综合影响的目标函数如下式：

minz＝w_a×t_a+w_b×t_b+w_c×t_c…(8)；

其中t_a,t_b,t_c分别为车内乘客乘车时间、变更路段站点乘客的等车时间及后面站点乘客的等车时间，具体表达式如式(9),(10),(11)；w_a,w_b,w_c为对应变量权重，由调度员依照层次分析法原理确定；

t_a＝m₁×(t_z+t_a2)…(9)；

其中，m₁车内乘客数，t_z为变更路段的交通阻抗，t_a2为后面路段的交通阻抗；T_g为下一公交车到达变更路段站点的时刻，T_ci为第i个乘客到达变更路段站点时刻，m₂为在公交车到达站点前到达的乘客数；T_gj为公交车到达后面第j个站点的时刻，T_ji为第i乘客到达第j站点的时刻，n为后面站点数，m_j为在公交车到达后面第j个站点到站的乘客数；

交通阻抗采用美国的BPR函数和交叉口影响修正系数确定,取α＝0.15,β＝4；

由上可得，表示调度综合影响的目标函数的自变量与交通阻抗有关，设该路段可供选择的行驶路径为自变量；并将其设为x＝(x₁,x₂,,…,x_H),H表示有H条路径可供选择；

其中x₁+x₂+…+x_H＝1,

自变量对应的交通阻抗向量为：t＝(t₁,t₂,…,t_H),

则变更路径的阻抗为：t_z＝x×t^T,

在表示出变更路径的交通阻抗后，结合预测的其他路段阻抗，乘客到站等信息，将目标函数表示为关于x的函数，利用模拟退火算法求出最优路径，并将其发送给驾驶员。

进一步，路径决策配套的方法包括：

路段划分与路径集的构建：

路段划分以是否存在可替代备选路径为标准进行，交叉口是公交车驶入、驶出规定线路的关键节点；该模型首先将该公交线路上的所有交叉口汇总形成交叉口集；交叉口之间的路径组成粗路径集，接着剔除粗路径集中高于3倍规定路径长度的备选路径，剩余路径组成路径集；

响应站点的设置：

采取变更路径决策时，若有站点被越过，则公交车辆在该响应站点停靠，被越过站点乘客在响应站点下车，乘客可以在响应站点等待下一班车，仅允许变更行驶路径公交车的后一辆公交车在响应站点停靠；响应站点以离被越过站点距离最近为原则进行设置；

当前一辆公交车越过该站点时，强制要求下一辆公交车必须途径该站点，避免该站点乘客等待过久。

本发明的另一目的在于提供城市公交动态调度优化系统包括：

同于结合红外线和图像识别技术，实时采集乘客到站时刻和数量、公交车内拥挤度、道路拥挤状况和公交车实时定位信息的实时数据获取系统；

用于建立时段车次数模型、发车时刻决策模型以及变更行驶路径模型的三层调度决策模型的调度决策系统；

用于实现人、车、路信息交互的信息共享系统。

本发明的优点及积极效果为：

本发明的公交调度模式以数据为导向，具有成本低廉、调度科学以及地区适应性强等优点，能有效提高公交的运营调度效率，具有广阔的应用前景。

本发明的时段车次数模型投资收益率、拥挤度满意率分别高0.8％、9.4％，候车时间满意率低3.5％。

本发明的发车时刻确定模型可有效提高换乘效率，在下述实例中，本模型至少降低了

的等车概率。

本发明的变更路径决策模型，在较少增加被越过站点乘客候车时间的基础上，极大地提高了公交车的运行速度和减少了后面站点乘客的候车时间。而且，发车间隔越密集、原始路径越拥堵、被越过站点乘客到达数越少、后面乘客到达数越多以及后面站点数越多，变更行驶路径的负面影响就越小、正面效应就越大。

本发明采集数据涵盖站台、车内及道路等客流数据，更加准确地了解客流分布状况。人、车、路实现信息交互，根据反馈的信息乘客可自主调节出行方式，公交企业实现精准调度。调度模式的选择的考虑了地区之间的差异，拓宽了调度的适用范围。建立了路径变更模型，可有效降低道路拥堵造成的危害。

本发明提供的模型特征应该包括相较于其他可变线路和越站模型，是否越过站点取决于道路拥堵状况和乘客到达情况。为解决越过站点乘客下车问题，提出响应站点的方案。路径变更也有限制，不管道路拥堵状况如何，不能进行连续变更路径。

附图说明

图1是本发明实施例提供的城市公交动态调度优化方法流程图。

图2是本发明实施例提供的城市公交动态调度优化系统示意图。

图中：1、实时数据获取系统；2、调度决策系统；3、信息共享系统。

图3是本发明实施例提供的车内及上下车客流数据获取原理流程图。

图4是本发明实施例提供的换乘站点分布示意图。

图5是本发明提供的城市公交动态调度优化系统中路径变更部分示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

下面结合附图对本发明的应用原理作详细描述。

如图1所示，本发明实施例提供的城市公交动态调度优化方法，包括：

S101：通过实时数据获取系统结合红外线和图像识别技术，实时采集乘客到站时刻和数量、公交车内拥挤度、道路拥挤状况和公交车实时定位信息；

S102：建立时段车次数模型、发车时刻决策模型以及变更行驶路径模型的三层调度决策模型；实现人、车、路信息交互。

如图2所示，本发明实施例提供的城市公交动态调度优化系统，包括：

实时数据获取系统1、调度决策系统2及信息共享系统3；实时数据获取系统同时结合红外线和图像识别技术，实时采集乘客到站时刻和数量、公交车内拥挤度、道路拥挤状况和公交车实时定位信息；

调度决策系统用于建立时段车次数模型、发车时刻决策模型以及变更行驶路径模型的三层调度决策模型；

信息共享系统用于实现人、车、路信息交互。

1)实时数据获取系统1：

本发明将实时数据获取分为三部分：(1)乘客到站时刻和数量；(2)公交车内拥挤度；(3)道路拥挤状况和公交车实时定位。数据获取的关键技术以图像识别为主，红外线、车载GPS等为辅，多种方法综合运用。

本发明结合红外线和图像识别技术：在公交站台除广告牌的其它三面分别安装红外线设备，顶棚一角安装摄像头。当有乘客到达站台时，红外线设备被触动，激发站台顶棚的摄像头立即拍摄此刻画面并记录时间，拍摄画面经过处理便能识别出站台内的候车人数。当有乘客持续挡在红外线发射器前，摄像设备不会连续拍摄画面，造成不必要的内存浪费。

车内数据的采集主要依靠双目摄像设备，分别在公交车前门和后门相应位置安装摄像头。当乘客上、下车时，摄像设备可清晰记录所有乘客的上下车数量和时间，运用计算机视觉技术，采集视频中人体头部形状为分析目标，得到车内拥挤度。

确定车内拥挤度的流程为：(1)双目摄像机捕捉到视频图像，发送到服务器；(2)服务器将运动物体从视频背景中检测、分离出来，进行筛选、过滤，获得真实的运动物体；(3)服务器将识别出来的人体精确地累加进出的人数，并叠加日期和时间，生成一条进出记录，可得到车内拥挤度。

当公交车到达起点站或终点站后，通过车内摄像头拍摄车内状况，进行图像识别，如若识别出车内有乘客，则进行乘客数的分析确认，否则将这一数据清零，继续下一班次的数据采集。其流程图3所示。

公交车载GPS每一时间间隔就发送当前位置数据至服务器。服务器先根据相邻两次经纬度，算出该间隔内车辆的行驶路程，再由公式：路程＝速度×时间，得出平均速度，时间间隔很小时，可将其视为瞬时速度。

路网数据信息的获取不仅对公交调度有重要作用，对于紧急状况发生的预测和处理也尤为重要，通过测得车辆速度可预测所在路线路网交通量，更有利于合理调度。路况数据的监测流程：(1)公交车将车载终端采集到的速度信息发送给处理器；(2)服务器计算出车辆的实时速度和位置，并将线路状况进行处理，通过红(非常拥堵)、黄(拥堵)、浅绿(常规)、绿色(畅通)来表示；(3)服务器将拥堵信息通过App等展示给乘客，并将车辆实时位置等反馈给调度系统。如表1.

表1路况信息采集反馈

2)时段车次确定模型：

公交客流量具有时间特性，故可由历史客流数据，将公交营运时间分为多个时段。为使企业和乘客整体利益最大化，本发明构建了投资收益率等三个指标。各地区存在差异，参照层次分析法原理，调度员可由实际情况确定指标权重。最后由指标和权重确定目标函数，自变量为时段发车车次数。如表2。

表2主要变量说明

由投资收益率等指标以及相应的权重，构建了目标函数如式(1)

max Z＝ω₁×μ₁+ω₂×μ₂+ω₃×μ₃…(1)

其中，μ₁,μ₂,μ₃分别表示投资收益率、候车时间满意率和拥挤度满意率，其具体表达式如下三式；ω₁,ω₂,ω₃分别为三项指标的权重，由层次分析原理求得。

其中：

自变量为时段车次数n(正整数)，找到了他的上下限过后，便可以编写程序求出最优解。通过财务分析和交通调查，确定公交公司可以忍受的最低满载率r_min该和线路乘客可以接受的最高满载率r_max。由式(5)可得时段车次数n的上下限，最后运用李发智的发车间隔模型和本模型对一实例进行计算，本模型投资收益率、拥挤度满意率分别高0.8％、9.4％，候车时间满意率低3.5％。

3)发车时刻优化模型：

考虑乘客换乘，本模型以总相遇次数最大和总相遇站点数最小为目标。相遇次数一定，减少相遇站点数可降低换乘时间。

构造换乘站点集：

不同的公交线路相互交汇，其交汇点往往就是换乘站点，这些换乘站点组成粗换乘站点集。但在现实生活中，相同两条公交线路可能存在2个及以上的换乘站点，但两线路同一班次公交车不可能在所有站点都相遇，现实情况是基本只能在一个换乘站点相遇。考虑乘客选择换乘站点的特性，当公交车相遇随机发生在各个换乘站点时，乘客选择对换乘站点的选择也呈随机分布；当将公交车相遇控制在同一换乘站点时，乘客也会自发的选择到该换乘站点换乘。如图4，1路乘客在换乘站点仅为A站点和A,B站点均分相遇次数情况下的等车概率分别为

因此，固定换乘站点可有效提高乘客换乘的便利性。本发明在粗换乘站点集的基础上，筛选出在其所在的两条公交线路的换乘集中，所能容纳同时停靠多辆公交车的站点，将该站点设置为该线路的换乘站点，这些换乘站点组成细换乘站点集。

在确定了细换乘站点集之后，先以使在该时段、该区域内两辆公交车在细换乘站点集内的站点总的相遇次数最大为目标，初步确定发车时刻；接着以使相遇点数最小(即使尽可能多的公交车在同一站点、相近时间相遇)调整公交发车时刻，确定最优发车时刻。具体模型如下：

以公交车区域内站点总相遇次数和所有线路总相遇点数作为目标函数，通过公式：

计算每两辆车总的相遇次数最大值；

其中：Z—总的相遇次数；

X_ki—线路k的第i辆车的发车时间；

X_qj—线路q上第j辆车的发车时间；

T_kn—线路k从起点到站点n的时间；

X_ki+T_kn+T_kin—线路k上第i辆车离开站点n的时间；

D_ntki—线路k第i辆车在时间t到达终点n取1，否则取0；

线路上总的相遇点的最小值确定每条线路上的最佳发车时间，根据公式

计算总相遇点数，Y_nt表示时间t到达站点n的车辆数

其中，M为该区域公交线路的条数；

F_k为该时段线路k的时段车次数；

该模型的约束条件包括：

1)时段第一辆公交车发车时间限制：

T_k1≤S_k；

其中，T_k1为线路k在该时段第一辆车的发车时间；

S_k为线路k所规定的该时段第一辆车的最迟发车时间；

2)时段最后一辆公交车发车时间限制：

其中，T为该时段结束时的时间点；

H_k该线路k所允许的最大发车间隔；

为该时段最后一辆车的发车时间；

F_k为该时段的发车车次数；

3)时段内相邻公交发车时间间隔限制：

h_k≤X_k(i+1)-X_ki≤H_k(k＝1,2,…,M,i＝1,2,…,F_k-1)；

其中，h_k表示线路k所允许的最小发车间隔；

X_ki表示线路k在该时段第i车次的发车时刻；

4)站点所能接受的最大相遇车辆数限制：

Y_nt≤R_n；

其中，R_n表示站点n所允许同时停靠的公交车辆数；

5)判断公交车是否在某时刻到达站点：

D_ntki＝max{0,(0.5-|t-(X_ki+T_kn)|)}；

其中：(n＝1,2,……,N；t＝1,2,……,T；k＝1,2,……M；i＝1,2,……,F_k)；

公交车实际到达的时刻与t相比，误差小于0.5分钟，则认为公交车到达。

对四条不同线路Ⅰ、Ⅱ、Ⅲ、Ⅳ)进行模拟分析，结果为选择线路Ⅰ在节点2换乘将会花费更久的时间，等车概率为

结合时段车次数模型，确定科学的发车时刻，缩小了乘客的换乘等待时间，提高了公交车调度效率，使得公交调度更加合理化。

4)路径决策模型：

该路径变更决策模式是在公交线路上发生的交通拥堵时，进行一定程度的行驶路径变更的一种新型城市公交调度模式。该调度模式在保留常规公交低成本性的基础上，提高了公交巴士运行的灵活性。远离拥堵路段可以提高公交巴士的行驶速度和降低后面站点(公交车还未到达的公交站点)乘客等待时间，满足乘客公交出行快捷性和及时性的需求，大大提高公交吸引力。

该系统主要由站台单元、车载单元、路侧单元以及调度控制中心等构成(这4者之间的关系如图5。站台单元、车载单元以及路侧单元分别采集站点乘客到达数、公交车位置以及道路交通状况等信息，并汇总形成数据库，调度控制中心根据信息，做出决策，实现公共巴士合理调度，降低交通拥堵产生的影响，提升公共巴士吸引力。

构造路径集与响应站点：

以是否存在可替代备选路径为标准进行路段划分，交叉口是公交车驶入、驶出规定线路的关键节点。因此，本文首先将该公交线路上的所有交叉口汇总形成交叉口集；交叉口之间的路径组成粗路径集，接着剔除粗路径集中高于3倍规定路径长度的备选路径，剩余路径组成路径集(路径越长，被采用的概率就越低，限制长度可提高运算效率)。

当公共巴士做出变更行驶路径的决策时，对应的原路径可能没有、可能存在1个或者多个公交车站的情况，没有公交车站的情况不必考虑。当存在公交车站时，公交车内可能有在被越过站点(原路径内的公交车站)下车的乘客，我们采用了设置响应站点的方案，公交车辆在该响应站点停靠，乘客在响应站点下车，乘客可以在响应站点等待下一班车(仅允许变更行驶路径公交车的后一辆公交车在响应站点停靠，不允许同一公交站点被连续两次越过)。乘客也可以选择步行或自行车出行到达目的站点，因为公交车变更行驶路径的决策是在原路径出现因交通拥堵的情况下做出的，步行等慢性交通与公交车行驶的时间相差不大，甚至会低于公交车行驶时间，响应站点也尽可能选在离被越过站点近的地方。

通过分析路径变化的影响，建立目标函数如下式：

minz＝w_a×t_a+w_b×t_b+w_c×t_c…(8)

其中t_a,t_b,t_c分别为车内乘客乘车时间、变更路段站点乘客的等车时间及后面站点乘客的等车时间，具体表达式如式(9),(10),(11)；w_a,w_b,w_c为对应变量权重，由调度员依照层次分析法原理确定。

t_a＝m₁×(t_z+t_a2)…(9)；

其中，m₁车内乘客数，t_z为变更路段的交通阻抗，t_a2为后面路段的交通阻抗；T_g为下一公交车到达变更路段站点的时刻，T_ci为第i个乘客到达变更路段站点时刻，m₂为在公交车到达站点前到达的乘客数；T_gj为公交车到达后面第j个站点的时刻，T_ji为第i乘客到达第j站点的时刻，n为后面站点数，m_j为在公交车到达后面第j个站点到站的乘客数。交通阻抗采用美国的BPR函数(取α＝0.15β＝4)和交叉口影响修正系数确定。

由上可得，目标函数的自变量与交通阻抗有关，故设该路段可供选择的行驶路径为自变量。并将其设为x＝(x₁,x_2,,…,x_H),H表示有H条路径可供选择。其中x₁+x₂+…+x_H＝1,

自变量对应的交通阻抗向量为：t＝(t₁,t₂,…,t_H),则变更路径的阻抗为：t_z＝x×t^T,在表示出变更路径的交通阻抗后，结合预测的其他路段阻抗，乘客到站等信息，可以将目标函数表示为关于x的函数，利用智能算法求出最优路径，并将其发送给驾驶员。

变更行驶路径的实例分析：

重庆607路公交车是由鸳鸯开往盘龙转盘的一票制常规公交线路。本文选取由金渝路公园1号至人和实验校这一路段(该路段包括龙宁路、金福路两站点)为例进行分析，龙健路口为公交路径与金渝大道的交叉口，位于金渝路公园1号与龙宁路，公交车可龙健路口驶出线路经金渝大道、龙政路越过龙宁路驶入线路，该路径设为备选路径；由龙健路口沿公交线路至驶入交叉口的路径设为原路径；驶入交叉口至金福路的路段称为后面路径1，金福路至人和实验校的路段称为后面路段2。

对该路径进行高峰时期的交通调查，获取信息如下：至龙健路口时公交车的载客量为21人，假设本班公交车较上一辆公交车在原路径的交通阻抗增加5min，龙健路口至龙宁路的交通阻抗约为原路径的15，其他路径交通阻抗基本不变。若上一辆公交车在原路径、后面路径1和后面路径2的平均交通阻抗为8min、6min和5min。备选路径的交通阻抗为8min。高峰时期公交车的平均发车间隔为5min。各站点乘客到达服从泊松分布，龙宁路、金福路以及人和实验校三个站点的平均乘客到达率分别为1(人/min)、1.2(人/min)和1.6(人/min)。利用以上数据和公式(1)、(2)和(3)求得是否变更行驶路径的三个决策指标的取值如下表：

表1两种方案决策指标取值

若各车内乘客乘车时间、越过站点候车时间以及后面站点候车时间的权重分别为0.2,0.4,0.4，则采用原路径和备选路径行驶的z值分别为225.6和147，在该决策指标中，z值越小，则乘客的整体利益最大，故在该情况下，采取变更行驶路径的决策方案。

由上可知，采用变更行驶路径的方案对乘客在公交站点的候车时间影响较大，大幅度增加了变更路段站点乘客的候车时间，但也大幅度减少了后面站点乘客的候车时间。分析表达式发现，影响程度与发车间隔、城市道路拥挤状况、公交站点乘客的到达率相关以及涉及站点数相关。发车间隔越密集、原始路径越拥堵、被越过站点乘客到达数越少、后面乘客到达数越多以及后面站点数越多，变更行驶路径的负面影响就越小、正面效应就越大。

本发明(1)采集数据涵盖站台、车内及道路等客流数据，更加准确地了解客流分布状况。(2)人、车、路实现信息交互，根据反馈的信息乘客可自主调节出行方式，公交企业实现精准调度。(3)调度模式的选择的考虑了地区之间的差异，拓宽了调度的适用范围。(4)基于换乘优先的发车时刻确定模型可有效降低换乘等车概率。(5)建立了路径变更模型，可有效降低道路拥堵造成的危害。

本发明采集数据种类全面、采集设备成本低，可有效反映客流情况；指标权重的确定考虑了地区差异，拓宽了应用范围，利于推广；作品以乘客为导向，可提高公交吸引力，并减少了调度员的工作量，提高运营调度水平，降低了人工成本；本作品还考虑了城市道路拥挤对公交运营的影响，对大中城市有较大的吸引力，具有良好的应用前景。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (4)

1.一种城市公交动态调度优化方法，其特征在于，所述城市公交动态调度优化方法包括：

在实时采集乘客到站时刻和数量、公交车内拥挤度、道路拥挤状况和公交车实时定位信息的基础上；建立时段车次数模型、发车时刻决策模型以及变更行驶路径模型的三层调度决策模型；并利用互联网通信技术，实现人、车、路信息交互；

所述时段车次数模型从平衡公交运营成本、乘客等车时间以及乘车舒适性出发，确定最优的时段车次数；

所述时段车次数模型得建立方法包括：

由投资收益率指标以及相应的权重，构建目标函数如式(1)

max Z＝ω₁×μ₁+ω₂×μ₂+ω₃×μ₃…(1)；

其中，μ₁,μ₂,μ₃分别表示投资收益率、候车时间满意率和拥挤度满意率，具体表达式如下式(2)、(3)、(4)；ω₁,ω₂,ω₃分别为三项指标的权重，由层次分析原理求得；

其中：

其中：ρ表示公交车票价，θ表示平均车次成本，m表示站点数，r_j表示第j路段满载率，p_j表示第j站点到达乘客总数，t_j表示公交车在j站点停靠的时间间隔，t_u,j表示上辆车到达j站点所花时间，a_j+1表示第j+1站点下车乘客数，s该时段的乘客总数，n表示时段车次数，d表示公交车标准载客量，q_j表示第j路段车内乘客数，t_jh表示第j站点乘客等待时间上限，r_h表示该线路乘客所能忍受的满载率上限，t_d,j表示公交车到达j站点所花的时间，T表示该时段长度；

自变量为时段车次数正整数n，找到自变量的上下限后，求出最优解；确定忍受的最低满载率r_min和该线路乘客接受的最高满载率r_max；由式(5)得时段车次数n的上下限；

在时段车次数确定的情况下，发车时刻决策模型以区域总的相遇次数最大和总的相遇站点数最小为目标，降低换乘等车概率，确定发车时刻；

变更行驶路径模型以城市道路经常发生交通拥堵为背景，绕开拥堵路段，提高车辆运行速度和降低后面站点乘客等车时间，实现车内乘客、变更路段站点候车乘客及后面站点候车乘客整体最优；

所述发车时刻决策模型的建立方法包括：

构造换乘站点集：不同的公交线路相互交汇，其交汇点往往就是换乘站点，这些换乘站点组成粗换乘站点集；在粗换乘站点集中，筛选出某两条线路中所能容纳同时停靠最多公交车辆数的站点，并将该站点设置为两条线路的换乘站点，这些换乘站点组成细换乘站点集；

在确定细换乘站点集之后，先以使在该时段、该区域内两辆公交车在细换乘站点集内的站点总的相遇次数最大为目标，初步确定发车时刻；接着以使相遇点数最小调整公交发车时刻，确定最优发车时刻；具体模型如下：

以公交车区域内站点总相遇次数和所有线路总相遇点数作为目标函数，通过公式：

计算每两辆车总的相遇次数最大值；

其中：Z—总的相遇次数；

X_ki—线路k的第i辆车的发车时间；

X_qj—线路q上第j辆车的发车时间；

T_kn—线路k从起点到站点n的时间；

X_ki+T_kn+T_kin—线路k上第i辆车离开站点n的时间；

D_ntki—线路k第i辆车在时间t到达终点n取1，否则取0；

线路上总的相遇点的最小值确定每条线路上的最佳发车时间，根据公式

计算总相遇点数，Y_nt表示时间t到达站点n的车辆数

其中，M为该区域公交线路的条数；

F_k为该时段线路k的时段车次数；

该模型的约束条件包括：

1)时段第一辆公交车发车时间限制：

T_k1≤S_k；

其中，T_k1为该时段线路k在该时段第一辆车的发车时间；

S_k为该时段线路k所规定的该时段第一辆车的最迟发车时间；

2)时段最后一辆公交车发车时间限制：

其中，T为该时段结束时的时间点；

H_k该时段线路k所允许的最大发车间隔；

为该时段线路k最后一辆车的发车时间；

F_k为该时段线路k的发车车次数；

3)时段内相邻公交发车时间间隔限制：

h_k≤X_k(i+1)-X_ki≤H_k(k＝1,2,…,M,i＝1,2,…,F_k-1)；

其中，h_k表示该时段线路k所允许的最小发车间隔；

X_ki表示该时段线路k在该时段第i车次的发车时刻；

4)站点所能接受的最大相遇车辆数限制：

Y_nt≤R_n；

其中，R_n表示站点n所允许同时停靠的公交车辆数；

5)判断公交车是否在某时刻到达站点：

D_ntki＝max{0,(0.5-|t-(X_ki+T_kn)|)}；

其中：(n＝1,2,……,N；t＝1,2,……,T；k＝1,2,……M；i＝1,2,……,F_k)；

公交车实际到达的时刻与t相比，误差小于0.5分钟，则认为公交车到达。

2.如权利要求1所述的城市公交动态调度优化方法，其特征在于，变更行驶路径模型的建立方法包括：

分析路径变化的影响，建立反应调度综合影响的目标函数如下式：

min z＝w_a×t_a+w_b×t_b+w_c×t_c…(8)

其中t_a,t_b,t_c分别为车内乘客乘车时间、变更路段站点乘客的等车时间及后面站点乘客的等车时间，具体表达式如式(9),(10),(11)；w_a,w_b,w_c为对应变量权重，由调度员依照层次分析法原理确定；

t_a＝m₁×(t_z+t_a2)…(9)；

其中，m₁车内乘客数，t_z为变更路段的交通阻抗，t_a2为后面路段的交通阻抗；T_g为下一公交车到达变更路段站点的时刻，T_ci为第i个乘客到达变更路段站点时刻，m₂为在公交车到达站点前到达的乘客数；T_gj为公交车到达后面第j个站点的时刻，T_ji为第i乘客到达第j站点的时刻，n为后面站点数，m_j为在公交车到达后面第j个站点到站的乘客数；

交通阻抗采用BPR函数和交叉口影响修正系数确定,取α＝0.15,β＝4；

表示调度综合影响的目标函数的自变量与交通阻抗有关，设该路段可供选择的行驶路径为自变量；并将其设为x＝(x₁,x_2,,…,x_H),H表示有H条路径可供选择；

其中x₁+x₂+…+x_H＝1,

自变量对应的交通阻抗向量为：t＝(t₁,t₂,…,t_H),

则变更路径的阻抗为：t_z＝x×t^T,

在表示出变更路径的交通阻抗后，结合预测的其他路段阻抗，乘客到站信息，将目标函数表示为关于x的函数，利用模拟退火算法求出最优路径，并将其发送给驾驶员。

3.如权利要求2所述的城市公交动态调度优化方法，其特征在于，变更行驶路径模型还包括：

路段划分与路径集的构建：

路段划分以是否存在可替代备选路径为标准进行，交叉口是公交车驶入、驶出规定线路的关键节点；该模型首先将该公交线路上的所有交叉口汇总形成交叉口集；交叉口之间的路径组成粗路径集，接着剔除粗路径集中高于3倍规定路径长度的备选路径，剩余路径组成路径集；

响应站点的设置：

采取变更行驶路径模型决策时，若有站点被越过，则公交车辆在该响应站点停靠，被越过站点乘客在响应站点下车，乘客可以在响应站点等待下一班车，仅允许变更行驶路径公交车的后一辆公交车在响应站点停靠；响应站点以离被越过站点距离最近为原则进行设置；

当前一辆公交车越过该站点时，强制要求下一辆公交车必须途径该站点，避免该站点乘客等待过久。

4.一种如权利要求1所述的城市公交动态调度优化方法的城市公交动态调度优化系统，其特征在于，所述城市公交动态调度优化系统包括：

实时数据获取系统、调度决策系统及信息共享系统；实时数据获取系统结合红外线和图像识别技术，实时采集乘客到站时刻和数量、公交车内拥挤度、道路拥挤状况和公交车实时定位信息；

调度决策系统用于建立时段车次数模型、发车时刻决策模型以及变更行驶路径模型的三层调度决策模型；

信息共享系统用于实现人、车、路信息交互。

Images