关键词 信号优先；快速公交；公交调度；行车时刻表；优化算法
1 引言
运营调度是公交企业的核心工作，一般分为行车时刻表编制、现场调度和统计分析三个阶段。其中，行车时刻表也叫行车计划，是组织和指导公交企业运营生产全过程的生产作业性计划，是公交运营调度的基础。对于行车时刻表优化，国内外许多学者做过相关研究。Scheele[1]提出了基于最小乘客出行时间的公交线路发车频率优化模型，Koutsopoulos等[2]提出了一个以乘客等待费用、运营成本和拥挤费用最小为目标的非线性规划模型。孙芙灵[3]探讨了不同客流状态下确定时段配车数和发车频率的方法，杨兆升[4]提出了一个最大化社会效益的发车频率优化模型。
研究表明，交叉口延误占公交整个行程时间的10%～20%，占整个延误时间的50%[5]。因此，只有降低了交叉口延误，才能显著降低公交车辆的行程延误，真正实现快速公交的快速。因此，本文提出了一种旨在降低公交车辆在交叉路口等待时间的快速公交行车时刻表优化模型，并在模型的基础上建立了行车时刻表优化算法。
2 快速公交运行描述
快速公交车辆有序地运行在一系列的站台与路口之间，在每个站台停车上下乘客，在需要的路口停车等待绿灯。考虑一条典型的快速公交线路（如图1所示），与传统公交线路相比，快速公交线路具有如下特征：
1) 公交专用道
公交专用道将其他交通与快速公交相隔离，给与了快速公交道路优先权，从而避免了其他交通流对快速公交运营速度的影响。
2) 交叉口公交信号优先
快速公交的交叉口一般设置公交专用信号灯，并对快速公交给与信号优先，从而提高快速公交的运营速度。
图1 快速公交线路示意图
对线路上的所有站台和交叉路口分别进行编号，在AB方向，站台顺序为<1>，<2>，……，<M-1>，<M>，在BA方向站台编号依次为<M+1>，<M+2>，……，<2M-2>，<2M-1>，<2M>。路口顺序在AB方向为<1>，<2>，……，<N-1>，<N>，在BA方向依次为<N+1>，<N+2>，……，<2N-2>，<2N-1>，<2N>。定义站台k和k+1之间的路段为路段k，则在AB方向上路段编号为<1>，<2>，……，<M-1>，在BA方向上路段编号为<M+1>，<M+2>，……，<2M-1>。
同样地，对行车时刻表中的所有车次按照发车时刻的顺序进行编号，假设共有K次行车计划，则按照发车时刻将行车时刻表中的所有计划依次编号为<1>，<2>，……，<K>。以上标来表示车次编号，下标表示站台编号，则第i次公交车到达第k个站台的时刻可以表示为（如图2所示）：
其中，Tk-1i为第i次车到达k-1站台的时刻；tsk-1i为第i次车在站台k-1的停靠时间；trk-1i为第i次车在路段k-1的行驶时间；twk-1i为第i次车通过路段k-1上的交叉口所用的时间。
图2 快速公交车辆运行示意图
即公交车辆从到达某一站点到到达下一站点的行程时间分为三部分，分别为站点服务乘客的停靠时间，所经路段上的行驶时间和途中通过交叉口的时间。即：
行程时间 = 站台停靠作业时间 + 路段行驶时间 + 通过交叉口时间。
3 快速公交行车时刻表优化模型
模型的目标是使公交车辆在交叉口的等待时间最小，但应同时兼顾乘客的利益和公交公司的收入。因此，选定考虑的因素为车辆在交叉口的等待时间、乘客的候车时间、车辆的满载率和企业的运营收入。下面首先分别对这四个因素建立数学模型，最终建立快速公交行车时刻表的优化模型。
3.1 车辆交叉口等待时间建模
对车辆交叉口等待时间，使用隶属度函数来表示。假设当公交车辆通过路口不需要等待绿灯，即等待时间为0时，隶属度为1；当车辆需要在路口停车等待绿灯时，隶属度随着车辆在路口等待时间的增加而降低。因此，使用函数μj(x)来表示第j辆车交叉口等待时间的隶属度：
式中：x为车辆的交叉口等待时间；qwj为第j次车的交叉口等待时间的期望值；m0为等待时间小于期望值qwj时的车辆交叉口等待时间隶属度值；Tsj为第j次车通过的所有路口的信号机红灯时长之和。其中，车辆的交叉口等待时间的计算公式如下：
其中，kj为第j次车通过的路口数量，wij为第j次车通过第i个路口的等待时间。
3.2 乘客候车时间建模
乘客在等待公交车辆时，都希望候车时间越短越好，因此，在模型中应考虑乘客的候车时间。为了计算的方便，采用候车时间小于耐受时间的乘客数与总乘客人数的比例，作为乘客候车时间的度量。
当公交车辆到达站台时，如果候车人数超过在耐受时间内到达的乘客数，则候车乘客中必定存在候车时间超过耐受时间的乘客，这一部分人数等于该站候车总人数减去在耐受时间内到达的乘客人数。计算公式如下：
其中，uij为第j次车在i站停靠时的候车超过耐受时间的乘客数；pij为第j次车到达i站时的候车人数；tij为第j次车到达i站的时间；L为乘客能够接受的最大候车时间；λi(t)为i站t时刻的乘客到达率(人/分钟)。
在给定初始状态的情况下，可以逐站计算出候车时间超过耐受时间的乘客数，从而得到乘客候车时间的隶属度计算公式如下，其中mj为第j次车经过的站台数量：
满载率是公交企业制定运营计划时的一个重要指标，如果满载率过高，将造成车上过于拥挤，服务质量下降，从而引起乘客不满；如果满载率过低，则浪费运力，增加了企业的运营成本。使用函数μj(x)来表示第j辆车的满载率隶属度：
式中：x为车辆的平均满载率；lej为第j次车的满载率指标；其中，车辆平均满载率x的计算公式如下：
其中，mj为第j次车经过的站台数量，lij为第j次车在i站台上下客后的满载率。
3.4 企业运营收入建模[6]
公交企业在运营满足客流需求的前提下，总是希望运营成本较小，企业收入较高。假设每发一辆车企业对收入都有一个最低期望值和一个满意期望值，当实际收入低于最低期望值时，企业收入隶属度为0；当实际收入大于满意期望值时，企业收入隶属度为1。因此，企业运营收入隶属度使用如下公式表示：
基于蒙特卡洛方法的高斯混合采样粒子滤波算法研究
朴素贝叶斯分类在入侵检测中的应用