CN106504516A - 一种基于公交站点信息化的多形式公交动态调度方法 - Google Patents

Abstract

本发明公开了一种基于公交站点信息化的多形式公交动态调度方法，包括如下步骤：(1)通过公交站点信息与乘客手机终端的信息连接和交互获得各种数据；(2)建立初始公交调度模型，对模型进行求解，得出包含多形式的公交车最佳初始发车时间间隔；(3)计算各种系数，确定每相邻两辆公交车的调度形式和运行区间；(4)建立实时多形式公交调度模型，对模型进行求解，得出每相邻两车间最佳的站点个数和站点间的最佳运行时间；(5)每相邻两车间的最佳运行时间间隔通过时间之差来求和。本发明实现了公交站点信息化技术与乘客手机终端的信息融合，实现了动态的不同公交形式和公交车最佳车间运行间隔的调度，提高了公交的运行效率。

Description

一种基于公交站点信息化的多形式公交动态调度方法

技术领域

本发明涉及智能公交控制领域，尤其是一种基于公交站点信息化的多形式公交动态调度方法。

背景技术

随着经济和社会的发展，机动车保有量迅速增加，城市交通拥堵现象日益严重，发展公共交通作为缓解交通压力的有效手段受到越来越多的重视，但是目前的公共交通系统开发多集中于单一调度形式，或者一成不变的调度模式，不能针对客流的实时变化，考虑不同的调度形式的组合以及高级的调度形式，还有车间的时间间隔也是公交的核心之一。其能够达到实时调度公交形式和车间的时间间隔的方案主要依据主要研究交线路客流信息，即乘客OD分布，为制定准确、完整的客流信息进行预测使用许多模型和方法，但其效果等还有待提高。

发明内容

本发明所要解决的技术问题在于，提供一种基于公交站点信息化的多形式公交动态调度方法，可以实现动态的公交形式和公交车最佳车间间隔的调度，提高了公交的运行效率的同时实现城市的公共交通系统良性循环发展。

为解决上述技术问题，本发明提供一种基于公交站点信息化的多形式公交动态调度方法，包括如下步骤：

(1)通过公交站点信息与乘客手机终端的信息连接和交互，得到实时乘客OD分布数据，通过公交车的GPS定位模块获取公交车的运行位置和时间数据；

(2)根据乘客OD分布得出每个站点上下车人数，计算公交车的站点停靠时间，结合相邻站点间运行时间，建立初始公交调度模型，其中公交车共有常用的三种不同运营形式，以单位时间内所有乘客总出行成本和企业运营成本之和最小为目标对模型进行求解，得出不同形式下公交车的最佳初始发车时间间隔；

(3)在每相邻两车间位于靠后位置的公交车离开某站点时刻，统计线路客流OD分布，计算线路客流方向不均衡系数、断面不均衡系数和站点不均衡系数，根据车辆调度形式的判断准则和方法确定出公交车的调度形式和运行区间；

(4)在每相邻两车间位于靠后位置的公交车离开某站点时刻，统计线路每个站点上下车人数，与位于靠前位置公交车之间达到最佳时间间隔的运行时段内每站点新增上下车人数，建立实时公交调度模型，以单位时间内所有乘客总出行成本和企业运营成本之和最小为目标对模型进行求解，得出每相邻两车间最佳的站点个数和站点间的最佳运行时间；

(5)每相邻两车间的最佳运行时间间隔是在位于靠后位置的公交车离开某站点时刻的原有车间运行时间间隔，与靠前位置公交车需调整的运行时间与靠后位置公交车需调整的运行时间之差的求和。

优选的，步骤(2)中初始公交调度模型的目标函数为：

其中：n—线路上站点总数目；k—线路上运营公交车总数目；—第j辆公交车离开i站的时刻；p—调度形式，p＝1为全程车(普通公交车)、p＝2为大站快车和p＝3区间车；—0-1变量，第j辆公交车在i站停靠，则为1，否则为0；—0-1变量，第j辆公交车经过从i-1站到i站区间，则为1，否则为0；—第j辆公交车在t时刻，乘客从i站至i+s站的人数(人/分钟)；一第j辆公交车离i-1站时车上的乘客数；Rⁱ—公交车从i-1站到i站的运行时间，开始输入历史数据的均值；—第j辆公交车在i站的停靠时间，取决于i站点的上车和下车的人数；—第j辆公交车到达i站点时下车的人数；—第j辆公交车到达i站点时上车的人数；β—单位乘客平均上下车时间(分钟/人)；C₁—乘客等车单位时间成本；C₂—乘客在车单位时间成本(元/分钟)；C₃—单位公交车运营时间成本(元/分钟)。

优选的，步骤(4)中公交调度模型的目标函数为：

其中：

优选的，步骤(2)中采用遗传算法对初始公交调度模型进行求解，包括如下步骤：

(2.1)初始化参数，参数包括：乘客出行OD分布、站间运行时间、各时间成本、初始种群大小、交叉概率、变异概率以及循环次数；

(2.2)对模型进行初始编码，令K′为最大的发车次数，将发车时间间隔和发车类型这两个决策变量对应的二进制串进行组合，得到一条8K’位二进制的染色体，前6K’位二进制染色体为时间间隔，后2K’位二进制代表公交车调度形式；根据最大发车间隔和最小发车间隔的约束条件，确定编码的长度区间；

(2.3)根据初始化确定的发车次数、发车间隔和发车形式，结合输入的参数计算目标函数值及每个个体的适应值；

(2.4)对种群进行繁殖，每繁殖一代包括选择、交叉、变异和基于约束条件进行修正；

(2.5)重复繁殖过程直到达到设置的迭代次数，选择适应值最大的个体代表的组合调度方案为初始发车间隔方案。

优选的，步骤(3)中对公交车的调度形式和运行区间进行求解，包括如下步骤：

(3.1)当断面不均衡系数大于一般规定值，即(一般取1.2-1.5)，路段客流量差ΔQ_i＞＞nq(n一般取2-4，q为公交车额定容量)，根据相应的准则，初步定出区间车和大站快车运行的路段及站点。

(3.2)当线路上站点不均衡系数大于一般规定值，即(一般取1.4-2.0)；方向不均衡系数大于一般规定值，即(一般取1.2-1.4)，根据相应的准则，初步定出大站快车运行的路段及站点。

优选的，步骤(4)中采用遗传算法对实时公交调度模型进行求解，包括如下步骤：

(4.1)初始化参数，参数包括：乘客出行OD分布，包括公交车离开某站点时刻、运营形式、以及该时刻线路乘客OD分布和每相邻公交车间运行时间间隔内新增乘客的OD分布、总站点数目、区间车和大站快车运行路段和站点、站间历史运行时间、各时间成本、初始种群大小、交叉概率、变异概率以及循环次数；

(4.2)对模型进行初始编码，令K′为最大的发车次数，将公交车的站间运行时间这个决策变量对应的二进制串进行组合，得到一条10K’位二进制的染色体，前六位表示公交车的站间运行时间和后四位表示最佳站点个数的取值；

(4.3)根据初始化确定的发车次数和发车间隔和发车形式，结合输入的参数计算目标函数值及每个个体的适应值；

(4.4)对种群进行繁殖，每繁殖一代包括选择、交叉、变异和修正；

(4.5)重复繁殖过程直到达到设置的迭代次数，选择适应值最大的个体代表的调度方案为实时动态调度策略，即选取每相邻两车间最佳的站点个数和站点间的最佳运行时间。

优选的，约束条件包括公交车满载率约束和发车时间间隔约束；所述满载率约束为高峰时段满载率不超过120％，平峰时段满载率不低于50％；所述发车时间间隔约束为不超过设定的最小发车间隔和最大发车间隔范围，且在研究时段T内满足公交车发车间隔X_j的总和不大于T。

优选的，步骤(5)中通过实时调度方法求解每相邻公交车间保持最佳的车间间隔运行的具体方法为：在j与j+1公交车间最佳运行时间间隔进行调整的模型中，需求解出j与j+1公交车停靠站点时间和站点间的运行时间Rⁱ，并求解出j+1与j+2公交车停靠站点时间和站点间的运行时间Rⁱ′，所以分别调整到最佳运行时间间隔的时刻为和所以j与j+1公交车需调整到的最佳运行时间间隔为时刻两车间的时间间隔与L^j+1和L^j作差之和，即

本发明的有益效果为：本发明实现了公交站点信息化技术与乘客手机终端的信息融合，避免其对乘客OD分布的预测，主动获取到实时乘客OD分布，实现了动态的公交形式和公交车最佳车间间隔的调度；通过车路协同保障其保持最佳的车间间隔运行，提高了公交的运行效率，同时也通过公交站点网络技术的广泛应用，即IOT(Internet-oriented-transit)，以公交站点网络技术来主动引导公共交通智能化发展，实现城市的公共交通系统有机良性循环发展，使居民或乘客能方便的选用公交的优先方式出行，最终增强公交吸引力。

附图说明

图1为本发明的求解初始发车间隔的原理图。

图2为本发明的求解路段上车间运行时间间隔的原理图。

图3为本发明的显示遗传算法的流程示意图。

图4为本发明的公交站点信息传递原理图。

具体实施方式

如图1、2、3和4所示，一种基于公交站点信息化的多形式公交动态调度方法，包括如下步骤：

(1)通过公交站点信息与乘客手机终端的信息连接和交互，得到实时乘客OD分布数据，通过公交车的GPS定位模块获取公交车的运行位置和时间数据；

(2)根据乘客OD分布得出每个站点上下车人数，计算公交车的站点停靠时间，结合相邻站点间运行时间，建立初始公交调度模型，其中公交车共有常用的三种不同运营形式，以单位时间内所有乘客总出行成本和企业运营成本之和最小为目标对模型进行求解，得出不同形式下公交车的最佳初始发车时间间隔；

(3)在每相邻两车间位于靠后位置的公交车离开某站点时刻，统计线路客流OD分布，计算线路客流方向不均衡系数、断面不均衡系数和站点不均衡系数，根据车辆调度形式的判断准则和方法确定出公交车的调度形式和运行区间；

(4)在每相邻两车间位于靠后位置的公交车离开某站点时刻，统计线路每个站点上下车人数，与位于靠前位置公交车之间达到最佳时间间隔的运行时段内每站点新增上下车人数，建立实时公交调度模型，以单位时间内所有乘客总出行成本和企业运营成本之和最小为目标对模型进行求解，得出每相邻两车间最佳的站点个数和站点间的最佳运行时间；

(5)每相邻两车间的最佳运行时间间隔是在位于靠后位置的公交车离开某站点时刻的原有车间运行时间间隔，与靠前位置公交车需调整的运行时间与靠后位置公交车需调整的运行时间之差的求和。

安装在站点的设备(设备与WiFi兼容)很容易和站点候车乘客的智能终端设备进行连接和交互，使其以最低成本获得大量信息。

能够在手机终端弹出在此站点经过的公交线路，能够提供目前公交车是否停靠或是途径的信息让乘客选择出行班次和方式。

此相关软件(app)或是通过手机短信的方式让乘客输入信息统计各站点的候车乘客OD分布，能够通过手机输入自己的OD数据，通过在停靠站附近的路侧设备传输到数据中心。乘客也可以通过智能终端获取公交车实时运行情况，例如下一班次及后续班次的位置及到站时间、车内剩余座位数目等。

在已得到乘客的OD分布，从而得出公交线路上由站点i上车至站点j下车的客流量。在公交车离开站点时会再次统计线路上乘客OD分布，通过在数据监测中心中进行储存和处理，即可得到实时公交线路动态OD数据。

因为不同车站和不同断面上客流不均衡性的客流的空间分布特性。并且包括公交网络中不同线路和一条线路不同断面的也有不同客流分布特性，本专利着重分析一条线路不同断面的客流空间分布特征。

城市快速公交线路的断面客流量，是指单位时间内，通过公交线路某一地点的客流量。在单位时间内，客流量最大的断面称为最大客流断面。为定量描述快速公交客流在空间上分布的不均衡性，选择断面不均衡系数、方向不均衡系数、站点不均衡系数三个指标来进行分析。

断面不均衡系数K_i和路段客流量差ΔQ_i，是单位时间内快速公交线路某断面客流量与平均断面客流量之比。可以显示出该线路在这个时刻各路段的客流动态，表示这种断面的不均衡性。

断面不均衡系数K_i，是单位时间内快速公交线路某断面客流量与平均断面客流量之比。K_i＝Q_i/Q_n，式中：K_s—为线路分时段不均衡系数，Q_si—为线路单位时间内第i个断面客流量，Q_n—为线路单位时间内平均断面客流量，路段客流量差ΔQ_i，指统计时间内运营线路某路段客流量与沿线各路段平均客流量之差。ΔQ_i＝Q_i-Q_n，式中：ΔQ_i—为线路分时段不均衡系数，Q_si—为线路单位时间内第i个断面客流量，Q_n—为线路单位时间内平均断面客流量。

快速公交线路上下行的客流量在同一时间内一般是不相等的。有的线路两个方向的运量几乎相等，而有的线路则差异很大。客流在方向上的不均衡性一般用方向不均衡系数表示，其计算方法是单位时间内某条线路最大单向断面客流量与双向平均客流量之比。

即：客流在站点上的不均衡系数K_j，指线路在j站乘客集散量与沿线各停车站平均乘客集散量之比。用来表示线路中各个站点客流量的不均衡性。

即：K_j＝Q_j/Q_n，式中：K_s—为线路各个站点不均衡系数，Q_si—为统计时间内站点i进、出站客流量，Q_n—为线路卑位时间内各站平均客流集散量。

为解决出现高峰时段乘客过分拥挤、平峰时段公交车满载率不高的现象，公交车的发车在符合客流分布需求，采用多种公交调度形式去满足需求，故此专利主要分析全程车、区间车、大站快车调度形式。通过上述可知公交线路的客流在空间上都具有不均衡性，为提高快速公交的运营效率针对不同的客流分布规律制定相应的公交车组合调度形式。

全程车、区间车、大站快车的具体概念，全程车：指快速公交车从线路起点发车直到终点站为止，必须在沿线按所有的固定停车站点依次停靠，并行驶满全程的一种基本调度形式，也称为慢车。区间车：公交车仅行驶线路上某一客流量较高的路段或区段的一种辅助调度形式。线路各断面的客流量往往是不均衡的，当客流断面的不均衡程度达到某一临界条件时,就应采用区间车的调度方法，从而使得断面客流量趋于均衡。大站快车：为适应公交线路乘距较长乘客的乘车需求，仅在沿线乘客集散量较大的停车站点停靠，是一种越站快速运行的调度形式。大站快车的调度方法即组织部分公交车只在大站进行停站，以减少停站次数，提高平均技术速度，从而节省时间。

首先初步分析客流分布，计算线路客流方向不均衡系数、断面不均衡系数和站点不均衡系数，根据公交车调度形式的判断准则和方法确定出每种公交车的运行区间。组合调度是由全程车、区间车和大站快车组成，根据客流情况在高峰时段高流量的线路上采用不同的组合方式(如全程车、全程车和区间车、全程车和大站快车、以及三种线路的组合)，提高系统运营效率。

具体步骤如下：根据公式分别计算线路上行和下行各断面不均衡系数或路段客流量差，线路上站点不均衡系数或方向不均衡系数。

当线路上站点不均衡系数大于一般规定值，即(一般取1.4-2.0)；断面不均衡系数大于一般规定值，即(一般取1.2-1.5)，路段客流量差ΔQ_i＞＞nq(n一般取2-4，q为公交车额定容量)，根据相应的准则，初步定出区间车和大站快车运行的路段及站点。综合考虑线路站距、折返点位置、道路条件是否满足有超车道以及调度工作方便等因素，拟定可行的运行路段及站点。

在运营过程中，系统的总成本主要体现在乘客和公交公司两方面。使得乘客出行成本和运营商成本的总和最小，来确定以多大的发车间隔、在路段以多大的时间间隔行驶和发何种类型的公交车的组合实现最优。

乘客的出行成本主要体现在公交完成出行过程所花费的时间上，即乘客从起始站点到目的站点的时间，分为车站等车时间、在车时间(包括在此过程中公交车在公交站点的停靠时间和交叉口延误时间、换乘时间)。

选取研究时段内所有公交车的运行时间来表示公交公司的成本，包括公交车在站间的运行时间和站点停靠时间。

乘客等车时间成本，本文快速公交车的发车形式分全程车、区间车、大站快车三种,当第j辆车是大站快车时，在非停靠站点没有乘客上下车、公交车没有停靠时间；当第j辆车是区间车时，在非运行区段没有站间运行时间，因此在的计算式中引入0-1变量，计算公式如

式中：—第j辆车离开i站的时刻；

—第j辆车在首站的发车时刻；

rⁱ公交车从i-1站到i站的运行时间，假定各种形式的公交车在相同站间的运行速度是一定的；

第j辆车在i站停留的时间；

p调度形式，p＝1为全程车、p＝2为大站快车、p＝3为区间车；

—0-1变量，第j辆车如果经过i站，则为1，否则为0。

公交车在停靠站点的停留时间与在站点上、下车的人数有关，乘客上下车在前后门打开，因此公交车的停站时间取决于人数多的一方。计算方法为：

式中：—第j辆车到达第i站时下车的人数；

—第j辆车到达第i站时上车的人数；

β—单位乘客平均上下车时间为(分钟/人)；

N—线路上站点数。

第j辆车在i站上车的人数即为在一定时间间隔内，到达i站点,前往线路下游运行区段各站的乘客数总和，即i站上车的人数同样，在i站下车的人数

式中：—0-1变量，第j辆车在i站停靠，则为1，否则为0。

乘客的在车时间通常分为两部分:乘客在乘车区间内公交车运行时间和在乘车区间内不下车的乘客在车站的停靠时间，其表达式ω₂＝C₂(t_y+t_z)

ω₂—乘客在车时间成本(元)；

C₂—乘客单位在车时间成本(元/min)；

t_y—公交车在运行过程中乘客的在车时间(min)；

t_z—不下车的乘客在车站的停靠时间(min)；

公交车j在运行过程中乘客的在车时间，由站间运行时间和车上的乘客数共同决定，表达式为：

式中:一第j辆车离i-1站时车上的乘客数

Rⁱ一公交车从i-1站到i站的运行时间

第j辆车离开i站时车上的乘客数等于第j辆车离开i-1站时车上的乘客数减去第j辆在第i站时下车的人数加上在第i站时上车的人数，计算公式如

在乘车区间内不下车的乘客在车站的停靠时间，由不下车的乘客数量和在该站点的停靠时间共同决定，表达式如

式中：—第j辆车在i站时没有下车的乘客数，其计算式

综上分析，考虑模型中是全程车、区间车和大站快车三种公交车的组合，当第j辆车是大站快车时，在非停靠站点公交车没有停靠时间；当第j辆车是区间车时，在非运行区段没有站间运行时间，因此，乘客的在车时间成本表达式如

公交车的运营时间h包括两部分，站点间运行时间和在站点的停靠时间。当发出公交车的类型不同时，沿途经过站点和停靠的站点不同。运营时段内所有公交车的总运行时间h可表示如当第j辆车是全程车即p＝1时，公交车运行经过所有站点并且在每个站点都停靠，所有的站间运行时间和站点的停靠时间都计算在单程运营时间内；第j辆车是大站快车即p＝2时，公交车经过所有站点并且只在特定的大站停靠，因此单程运营时间为站间运行时间与停靠站点的停车时间之和，对于i不属于M，M为除了大站停靠的其他站点集合，当第j辆车是区间车即p＝3时，公交车在公交线路的部分路段运行，只经过部分站点，因此单程运营时间只是站间运行时间和站点的停车时间之和，区间车不经过的站点其区间运行时间为零，因没有乘客上下车故停靠时间也为零。设单位公交车运营时间成本C₃(元/min)，则公交车运营时间成本表达式

乘客出行时间成本和运营商旳运营成本两个目标都是要达到最小值，因此，对两个目标进行线性加和，以总成本最小作为目标函数，在固定需求条件下，单位时间内所有乘客总出行成本和企业运营成本之和最小。考虑两部分成本是同等重要的，取相同的权重，目标函数min Z＝C_ck+C_y。

乘客出行时间成本由两部分组成，包括乘客等车时间成本(C₁)和在车时间成本(C₂)，对于公交公司成本以公交车在单程调度过程中所花费的时间成本(C_y)衡量。所以初始的目标函数为：

在明确组合调度模型的目标函数后，还需要设定模型的约束条件，才能得到合理的调度组合方案。在该模型中约束条件主要考虑公交车乘客的满载率方面、出行时间约束方面和发车间隔方面等因素。

公交车满载率是指公交车在线路某断面载运乘客的平均满载程度，标志着城市公交管理的状况和运行质量的高低。满载率计算公式如下：满载率＝公交车j离开车站i时车上乘客数/公交车额定载客人数*100％。一般情况下从乘客舒适性和安全性角度考虑，高峰时段满载率不应超过120％，同时考虑公交企业的经济效益平峰时满载率一般也不能低于50％。

发车间隔的限制是指相邻两车从首站发出的间隔时间需要在一定范围内，如果发车间隔过小，发车次数增加，公交车的运营效益无法满足。发车间隔过大，容易导致服务水平降低，引起乘客不满意。本模型考虑乘客与高峰时段发车间隔可接受程度，并结合以往历史发车时间间隔，确定最大发车间隔和最大发车间隔。

在研究时段T内需要满足公交车发车间隔的总和不大于T，则在研究时段内公交公司的发车总次数K′的公式为：从建立的公交调度的模型来看，需要确定的是在研究时间内，选择的调度形式，保证公交出行乘客的出行成本和公交公司的运营成本最低。需要确定的初始发车时间间隔与发车总数是相关的，即K是不确定的，与决策变量相关，无法将两者脱离开来。模型中涉及到了很多的参数，这些参数之间的关系错综复杂，并且参数的确定过程中用到了决策变量，所以选取遗传算法进行求解，求解过程需要将决策变量用编码的形式表示出来，由此可以获得决策变量的可能解，从而推出相关的参数。

初始化参数，其一为模型固有参数，包括：乘客出行OD分布(包括公交车离开站点时刻统计线路上乘客的OD分布)、站间历史运行时间、乘客等车单位时间成本、运行时间成本等，其二为模型运行参数，如：初始种群大小、交叉概率、变异概率以及循环次数等。

对模型进行初始编码，令K′为最大的发车次数，将发车时间间隔和发车类型这两个决策变量对应的二进制串进行组合，得到一条8K′位二进制的染色体，前6K′位二进制染色体为时间间隔，后2K′位二进制代表公交车调度形式；根据最大发车间隔和最小发车间隔的约束条件，确定编码的长度区间。

根据初始化确定的发车次数和发车间隔，结合输入的参数计算目标函数值及每个个体的适应值。

对种群进行繁殖，每繁殖一代包括选择、交叉、变异和基于约束条件进行修正。

重复繁殖过程直到达到设置的迭代次数，趋于稳定，最终输出满意解，选择适应值最大的个体代表的组合调度方案为初始发车间隔和发车组合调度方案。由于公交线路上的公交客流是动态变化的，所以对每辆公交车离开某一个公交站点后需记录该时刻线路全部乘客OD分布，实时更新客流数据，从而研究客流分布与变化规律，所以每相邻公交车的发车间隔需要做出调整。

具体实时调度过程如下：

(1)在每相邻j和j+1两辆公交车间位于靠后位置的公交车离开某站点时刻，统计线路客流OD分布，计算线路客流方向不均衡系数、断面不均衡系数和站点不均衡系数三个指标，根据车辆调度形式的判断准则和方法确定出公交车的调度形式和运行区间。当断面不均衡系数大于一般规定值，即(一般取1.2-1.5)，路段客流量差ΔQ_i＞＞nq(n一般取2-4，q为公交车额定容量)，根据相应的准则，初步定出区间车运行的路段及站点。当线路上站点不均衡系数大于一般规定值，即(一般取1.4-2.0)；方向不均衡系数大于一般规定值，即(一般取1.2-1.4)，根据相应的准则，初步定出大站快车运行的路段及站点。80综合考虑线路站距、折返点位置、道路条件是否满足有超车道以及调度工作方便等因素，拟定可行的运行路段及站点。

(2)当线路判断存在大车站点时，需记录此时的大车站点顺序号，数目为1,2,3…W，并且记录W个站点的乘客OD分布，计算每相邻j和j+1公交车间最佳的运行时间间隔只需输入在两车之间的大车站点顺序号和其站点OD分布，其他站点为乘客数为0，也为0。

(3)当线路判断存在区间站点时，划分出相应的区间，需记录此时区间内站点顺序号，数目为1,2,3…Q，并且记录Q个站点的乘客OD分布，当有公交车两车之间的存在有区间站点时，则公交车实施区间车的调度形式，只在区间站点停靠。需计算每相邻j和j+1公交车间最佳的运行时间间隔只需输入在两车间的区间站点顺序号和其站点OD分布。

(4)当线路判断不存在区间站点和大车站点时，每相邻j和j+1公交车不实施区间快车和大车快车的调度方式，则公交车只能实施普通公交的调度形式。

(5)此时刻需对j和j+1两辆公交车的调度形式做出调整，如果下一时刻断面不均衡系数、方向不均衡系数、站点不均衡系数发生改变，按此时刻的条件调整，两辆公交车的调度形式将被重新调整，以满足该时刻的公交车调度计划。

(6)把j与j+1公交车作为相邻研究的公交车，当j+1公交车离开第i+s站点时刻j+2公交车行到i+s+k至i+s+k+1站点区间(包括刚好离开i+s+k站点时)，记录j+1与j+2公交车在此时刻的位置，并重新记录线路上每站点的乘客OD分布，若j+1与j+2公交车之间相距的站点为k个，线路上每站点新的乘客OD分布为j+1公交车从此时刻位置至j+2公交车在此时刻位置的运行时间段内，第i+s+1至i+s+k站点(包括i+s+k站点)间k个站点的乘客OD分布数目与每站点新增乘客数目之和，当高峰时公交车内会出现乘客站立的情况，所以公交车承载的乘客数目(乘客承载能力)扩大1.2倍。

(7)若该运行时间段内各站点的新增乘客数目为0，最新计算的公交车承载乘客数目如果大于公交车的剩余座位数目(剩余承载力)，说明j+1公交车没有剩余承载能力，则取公交车的剩余座位数目作为上限约束条件，计算满足j+1公交车最大的承载率时需运行经过站点k2个公交站点，应到达的站点i+s+k2，所以k2个站点可以看作目标函数求解数值的上限。

(8)若取公交车的最新计算的公交车承载乘客数目作为上限约束条件，计算满足j+1公交车最大的承载率时需运行经过站点k22个公交站点，应到达的站点i+s+k22，所以k22个站点可以看作目标函数求解数值的上限。

(9)若该运行时间段内各站点新增乘客数目取历史数据的最大值，最新计算的公交车承载乘客数目如果大于公交车的剩余座位数目(剩余承载力)，说明j+1公交车没有剩余承载能力，则取公交车的剩余座位数目作为下限约束条件，计算满足j+1公交车最大的承载率时需运行经过站点k1个公交站点，应到达的站点i+s+k1，所以k1个站点可以看作目标函数求解数值的下限。

(10)若取公交车的最新计算的公交车承载乘客数目作为下限约束条件，计算满足j+1公交车最大的承载率时需运行经过站点k12个公交站点，应到达的站点i+s+k12，所以k12个站点可以看作目标函数求解数值的下限。

(11)实际情况下，只要新增乘客OD分布的始点或终点落在i+s+1至i+s+k区间(包括i+s+k)的任意站点，则根据目标函数求解公交车相距的站点数目k的最优值，根据该运行时间段内各站点新增乘客数目，判断最新计算的公交车承载乘客数目是否大于公交车的剩余座位数目(剩余承载力)，如果最新计算的公交车承载乘客数目大于公交车的剩余座位数目，可以得到j+1与j+2公交车相距的站点数目K的取值范围k1≤h≤k2，如果最新计算的公交车承载乘客数目小于公交车的剩余座位数目，可以得到j+1与j+2公交车相距的站点数目K的取值范围k12≤h≤k22。

(12)在j+1公交车离开第i+s站点时刻j车运行到i-1至i站点区间(包括刚好停靠在i站点时)，记录j与j+1公交车在此时刻的位置，并重新记录线路上每站点的乘客OD分布，若j与j+1公交车之间相距的站点为s个，线路上每站点新的乘客OD分布为j公交车从此时刻位置至j+1公交车在此时刻位置的运行时间段内，第i+1至i+s站点(包括i+s站点)间s个站点的乘客OD分布数目与每站点新增乘客数目之和，当高峰时公交车内会出现乘客站立的情况，所以公交车承载的乘客数目(乘客承载能力)扩大1.2倍。

(13)若该运行时间段内各站点的新增乘客数目为0，最新计算的公交车承载乘客数目如果大于公交车的剩余座位数目(剩余承载力)，说明j公交车没有剩余承载能力，则取公交车的剩余座位数目作为上限约束条件，计算满足j公交车最大的承载率时需运行经过站点h2个公交站点，应到达的站点i+h2，所以h2个站点可以看作目标函数求解数值的上限。

(14)若取公交车的最新计算的公交车承载乘客数目作为上限约束条件，计算满足j公交车最大的承载率时需运行经过站点h22个公交站点，应到达的站点i+h22，所以h22个站点可以看作目标函数求解数值的上限。

(15)若该运行时间段内各站点新增乘客数目取历史数据的最大值，最新计算的公交车承载乘客数目如果大于公交车的剩余座位数目(剩余承载力)，说明j公交车没有剩余承载能力，则取公交车的剩余座位数目作为下限约束条件，计算满足j公交车最大的承载率时需运行经过站点h1个公交站点，应到达的站点i+h1，所以h1个站点可以看作目标函数求解数值的下限。

(16)若取公交车的最新计算的公交车承载乘客数目作为下限约束条件，计算满足j公交车最大的承载率时需运行经过站点h12个公交站点，应到达的站点i+h12，所以h12个站点可以看作目标函数求解数值的下限。

(17)实际情况下，只要新增乘客OD分布的始点或终点落在i+1至i+k区间(包括i+k)的任意站点，则根据目标函数求解公交车相距的站点数目h的最优值，根据该运行时间段内各站点新增乘客数目，判断最新计算的公交车承载乘客数目是否大于公交车的剩余座位数目(剩余承载力)，如果最新计算的公交车承载乘客数目大于公交车的剩余座位数目，可以得到j与j+1公交车相距的站点数目s的取值范围h1≤h≤h2，如果最新计算的公交车承载乘客数目小于公交车的剩余座位数目，可以得到j与j+1公交车相距的站点数目h的取值范围h12≤h≤h22。

第j+1辆公交车从此时的位置至i+s+k站点运行时间段内，考虑乘客出行时间成本和公交公司运营成本两个目标都是要达到最小值。因此，对两个目标进行合并，采用线性加和以总成本最小作为目标函数。即在约束条件下，单位时间内所有乘客总出行成本和企业运营成本之和最小。考虑两部分成本是同等重要的，取相同的权重，目标函数min Z＝C_ck+C_y。

乘客出行时间成本由两部分组成，包括乘客等车时间成本(ω₁)和在车时间成本(ω₂)，对于运营商成本本文以公交车在单程调度过程中所花费的时间成本(C_y)衡量。所以最终的目标函数为：

式中：

同理，第j辆公交车从此时的位置至i+h站点运行时间段内，考虑乘客出行时间成本和公交公司运营成本两个目标都是要达到最小值。所以最终的目标函数为：

在明确第j辆与第j+1辆公交车调度模型的目标函数后，还需要设定模型的约束条件，才能得到合理的调度组合方案。在该模型中约束条件主要考虑公交车乘客的满载率方面、出行时间约束方面和发车间隔方面等因素。

参数设置包含两方面，其一为模型固有参数，如:乘客出行OD分布(包括公交车离开站点时刻统计的线路乘客的OD分布和调整车间时间间隔的新增乘客的OD分布)、总站点数目、、区间车和大站快车运行路段和站点、站间运行时间、乘客等车单位时间成本，公交车在单程调度过程中所花费的时间成本等，其二为模型运行参数，如:初始种群大小、交叉概率、变异概率以及循环次数等。

通过模型制定出j和j+1公交车与j+1和j+2的在调整时间间隔中，通过上述遗传算法求解出j和j+1公交车需要停靠站点区间运行时间Rⁱ，是公交站点间的运行时间的最佳取值，调整所需的时间所以j与j+1的最佳运行时间间隔是在时刻的车间时间再加上L^j+1与L^j之差的和，即

尽管本发明就优选实施方式进行了示意和描述，但本领域的技术人员应当理解，只要不超出本发明的权利要求所限定的范围，可以对本发明进行各种变化和修改。

Claims (8)

1.一种基于公交站点信息化的多形式公交动态调度方法，其特征在于，包括如下步骤：

(1)通过公交站点信息与乘客手机终端的信息连接和交互，得到实时乘客OD分布数据，通过公交车的GPS定位模块获取公交车的运行位置和时间数据；

(2)根据乘客OD分布得出每个站点上下车人数，计算公交车的站点停靠时间，结合相邻站点间运行时间，建立初始公交调度模型，其中公交车共有常用的三种不同运营形式，以单位时间内所有乘客总出行成本和企业运营成本之和最小为目标对模型进行求解，得出不同形式下公交车的最佳初始发车时间间隔；

(3)在每相邻两车间位于靠后位置的公交车离开某站点时刻，统计线路客流OD分布，计算线路客流方向不均衡系数、断面不均衡系数和站点不均衡系数，根据车辆调度形式的判断准则和方法确定出公交车的调度形式和运行区间；

(4)在每相邻两车间位于靠后位置的公交车离开某站点时刻，统计线路每个站点上下车人数，与位于靠前位置公交车之间达到最佳时间间隔的运行时段内每站点新增上下车人数，建立实时公交调度模型，以单位时间内所有乘客总出行成本和企业运营成本之和最小为目标对模型进行求解，得出每相邻两车间最佳的站点个数和站点间的最佳运行时间；

(5)每相邻两车间的最佳运行时间间隔是在位于靠后位置的公交车离开某站点时刻的原有车间运行时间间隔，与靠前位置公交车需调整的运行时间与靠后位置公交车需调整的运行时间之差的求和。

2.如权利要求1所述的基于公交站点信息化的多形式公交动态调度方法，其特征在于，步骤(2)中初始公交调度模型的目标函数为：

$\begin{array}{c}\min Z=C_1\underset{j=1}{\overset{k}{\mathrm{\Σ}}}\underset{i=1}{\overset{n-1}{\mathrm{\Σ}}}\underset{s=1}{\overset{n-i}{\mathrm{\Σ}}}{\∫}_{T_i^{j1}}^{T_i^j}{\mathrm{\γ}}_{i,j}^p{\mathrm{\γ}}_{i+s,j}^p{D}_{i,i+s}^j\left(t\right)(T_i^j+t)dt+C_2\underset{j=1}{\overset{k}{\mathrm{\Σ}}}\underset{i=2}{\overset{n}{\mathrm{\Σ}}}({\mathrm{\δ}}_{i,j}^p{S}_{i-1}^j{R}^i\\ +{\mathrm{\γ}}_{i,j}^p{\mathrm{\τ}}_i^j(S_{i-1}^j-\underset{u=1}{\overset{i-1}{\mathrm{\Σ}}}{B}_{d,i}^j))+C_3\underset{j=1}{\overset{k}{\mathrm{\Σ}}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}({\mathrm{\δ}}_{i,j}^p{R}^i+{\mathrm{\γ}}_{i,j}^p{\mathrm{\τ}}_i^j)\end{array}$

其中：n—线路上站点总数目；k—线路上运营公交车总数目；—第j辆公交车离开i站的时刻；p—调度形式，p＝1为全程车(普通公交车)、p＝2为大站快车和p＝3区间车；—0-1变量，第j辆公交车在i站停靠，则为1，否则为0；—0-1变量，第j辆公交车经过从i-1站到i站区间，则为1，否则为0；—第j辆公交车在t时刻，乘客从i站至i+s站的人数(人/分钟)；一第j辆公交车离i-1站时车上的乘客数；Rⁱ—公交车从i-1站到i站的运行时间，开始输入历史数据的均值；—第j辆公交车在i站的停靠时间，取决于i站点的上车和下车的人数；—第j辆公交车到达i站点时下车的人数；—第j辆公交车到达i站点时上车的人数；β—单位乘客平均上下车时间(分钟/人)；C₁—乘客等车单位时间成本；C₂—乘客在车单位时间成本(元/分钟)；C₃—单位公交车运营时间成本(元/分钟)。

3.如权利要求1所述的基于公交站点信息化的多形式公交动态调度方法，其特征在于，步骤(4)中公交调度模型的目标函数为：

$\begin{array}{c}\min Z=C_1\underset{a=i+s+1}{\overset{i+s+x-1}{\mathrm{\Σ}}}\underset{s1=1}{\overset{n-i-x}{\mathrm{\Σ}}}{\∫}_{T_a^{(j+1)1}}^{T_a^{j+1}}{\mathrm{\γ}}_{a,j}^p{\mathrm{\γ}}_{a+s1,j}^p{D}_{a,a+s1}^{j+1}\left(t\right)(T_a^{j+1}-t)dt\\ +C_2\underset{a=i+s+1}{\overset{i+s+x}{\mathrm{\Σ}}}({\mathrm{\δ}}_{a,j+1}^p{S}_{a-1}^{j+1}{R}^a+{\mathrm{\γ}}_{a,j+1}^p{\mathrm{\τ}}_a^{j+1}(S_{a-1}^{j+1}-\underset{d=1}{\overset{i+s}{\mathrm{\Σ}}}{B}_{d,a}^{j+1}))\\ +C_3\underset{a=i+s+1}{\overset{i+s+x}{\mathrm{\Σ}}}({\mathrm{\δ}}_{a,j+1}^p{R}^a+{\mathrm{\γ}}_{a,j+1}^p{\mathrm{\τ}}_a^{j+1})\end{array}$

其中：

$T_a^{j+1}={\mathrm{\γ}}_{a,j+1}^p\{{\mathrm{\Σ}}_{i=2}^a{R}^i+T_1^{j+1}+{\mathrm{\Σ}}_{i=2}^a{\mathrm{\τ}}_i^{j+1}\}.$

4.如权利要求1所述的基于公交站点信息化的多形式公交动态调度方法，其特征在于，步骤(2)中采用遗传算法对初始公交调度模型进行求解，包括如下步骤：

(2.1)初始化参数，参数包括：乘客出行OD分布、站间运行时间、各时间成本、初始种群大小、交叉概率、变异概率以及循环次数；

(2.2)对模型进行初始编码，令K′为最大的发车次数，将发车时间间隔和发车类型这两个决策变量对应的二进制串进行组合，得到一条8K’位二进制的染色体，前6K’位二进制染色体为时间间隔，后2K’位二进制代表公交车调度形式；根据最大发车间隔和最小发车间隔的约束条件，确定编码的长度区间；

(2.3)根据初始化确定的发车次数、发车间隔和发车形式，结合输入的参数计算目标函数值及每个个体的适应值；

(2.4)对种群进行繁殖，每繁殖一代包括选择、交叉、变异和基于约束条件进行修正；

(2.5)重复繁殖过程直到达到设置的迭代次数，选择适应值最大的个体代表的组合调度方案为初始发车间隔方案。

5.如权利要求1所述的基于公交站点信息化的多形式公交动态调度方法，其特征在于，步骤(3)中对公交车的调度形式和运行区间进行求解，包括如下步骤：

(3.1)当断面不均衡系数大于一般规定值，即(一般取1.2-1.5)，路段客流量差ΔQ_i＞＞nq(n一般取2-4，q为公交车额定容量)，根据相应的准则，初步定出区间车和大站快车运行的路段及站点；

(3.2)当线路上站点不均衡系数大于一般规定值，即(一般取1.4-2.0)；方向不均衡系数大于一般规定值，即(一般取1.2-1.4)，根据相应的准则，初步定出大站快车运行的路段及站点。

6.如权利要求1所述的基于公交站点信息化的多形式公交动态调度方法，其特征在于，步骤(4)中采用遗传算法对实时公交调度模型进行求解，包括如下步骤：

(4.1)初始化参数，参数包括：乘客出行OD分布，包括公交车离开某站点时刻、运营形式、以及该时刻线路乘客OD分布和每相邻公交车间运行时间间隔内新增乘客的OD分布、总站点数目、区间车和大站快车运行路段和站点、站间历史运行时间、各时间成本、初始种群大小、交叉概率、变异概率以及循环次数；

(4.2)对模型进行初始编码，令K′为最大的发车次数，将公交车的站间运行时间这个决策变量对应的二进制串进行组合，得到一条10K’位二进制的染色体，前六位表示公交车的站间运行时间和后四位表示最佳站点个数的取值；

(4.3)根据初始化确定的发车次数和发车间隔和发车形式，结合输入的参数计算目标函数值及每个个体的适应值；

(4.4)对种群进行繁殖，每繁殖一代包括选择、交叉、变异和修正；

(4.5)重复繁殖过程直到达到设置的迭代次数，选择适应值最大的个体代表的调度方案为实时动态调度策略，即选取每相邻两车间最佳的站点个数和站点间的最佳运行时间。

7.如权利要求3或5所述的基于公交站点信息化的多形式公交动态调度方法，其特征在于，约束条件包括公交车满载率约束和发车时间间隔约束；所述满载率约束为高峰时段满载率不超过120％，平峰时段满载率不低于50％；所述发车时间间隔约束为不超过设定的最小发车间隔和最大发车间隔范围，且在研究时段T内满足公交车发车间隔X_j的总和不大于T。

8.如权利要求1所述的基于公交站点信息化的多形式公交动态调度方法，其特征在于，步骤(5)中通过实时调度方法求解每相邻公交车间保持最佳的车间间隔运行的具体方法为：在j与j+1公交车间最佳运行时间间隔进行调整的模型中，需求解出j与j+1公交车停靠站点时间和站点间的运行时间Rⁱ，并求解出j+1与j+2公交车停靠站点时间和站点间的运行时间Rⁱ′，所以分别调整到最佳运行时间间隔的时刻为 和所以j与j+1公交车需调整到的最佳运行时间间隔为时刻两车间的时间间隔与L^j+1和L^j作差之和，即