Optimal Feeder Flex-route Transit Schedule Design of Urban Rail Transit
-
摘要:
为尽可能满足乘客出行需求,制定科学合理的公交运行时刻表能够保持公交稳定运行的同时降低运行成本. 首先,针对可变线路公交区别于常规定点定线公交的运行特征,提出一种将常规公交转变为可变线路公交接驳地铁的策略,在无需开设新接驳线路的前提下服务原有公交车出行需求乘客和短途出行需求的地铁乘客;然后,通过分析乘客出行行为的规律与特点,根据可变线路公交和地铁的时空信息以及乘客换乘行为的时空参数,利用混合整数非线性规划构建可变线路公交接驳地铁的路径、时间、调度协同优化模型和算法;最后,结合案例,分析车辆运行速度与出行需求水平以及不同时间成本侧重情况对发车间隔、车辆运行时间、运行成本等系统各项指标的影响. 研究表明:在运行周期内适当增加发车频率可以有效降低乘客出行时间成本,且不会增加系统总成本;为可变线路公交接驳地铁的时刻表构建提供了基本理论和自动化编制方法,提高了公交运行效率和服务质量.
Abstract:A scientific schedule should be able to meet the needs of most passengers, maintain stable operation, and reduce operating costs. Aiming at the operation characteristics of the flex-route transit different from the traditional fixed-route transit, a method of making the schedule of the feeder flex-route transit of urban rail transit is proposed. A strategy of traditional fixed-route transit into flex-route transit to connect urban rail transit is proposed, which can serve the original traditional fixed-route passengers and urban rail transit passengers with short-distance travel demand at the same time without the need to open new feeder routes. The rules and characteristics of passenger travel behavior are analyzed. According to the space-time parameters of transfer behavior of passengers and the space-time information of flex-route transit and urban rail transit, the mixed integer nonlinear programming is used to construct the path, time, and scheduling collaborative optimization model and algorithm of the feeder flex-route transit of urban rail transit. Combined with the case analysis, the influence of vehicle driving speed and travel demand level and different time cost emphasis on the system indicators such as departure interval, vehicle running time, and operation cost are analyzed. The results indicate that increasing the departure frequency appropriately during the operating cycle can effectively reduce the travel time cost of passengers without increasing the total system cost. It provides basic theory and automatic compilation method for the feeder flex-route transit schedule design of urban rail transit, reduces the travel cost and improves the operation efficiency and service quality.
-
Key words:
- Flex-route transit /
- Urban rail transit /
- Schedule /
- Transfer analysis /
- Operations coordination
-
随着出行需求的不断增长和出行选择日渐向快捷、舒适方向发展,私家车或网约车使用率提高,造成交通拥堵[1],同时也导致公共交通特别是公交车的使用率逐渐下降. 因此,如何最大程度地发挥公共交通的优势,并满足各类灵活出行需求,是迫切需要深入研究的问题[2-3].
公交车接驳大型轨道交通线路或站点不仅可以解决乘客的“最后一公里”问题,还能够维持公交车一定的运力. 在公交车接驳轨道交通时刻表编制方面,Wang等[4]提出了一种公交接驳地铁发车时刻表和车辆调度一体化优化模型;窦雪萍等[5]构建了公交接驳轨道交通时刻表并与车队规模联合的优化模型,以减少乘客排队等待出行成本与公交企业运营成本;Xiong等[6]提出了关于社区班车与地铁服务协同优化的研究问题,使用遗传算法和Frank-Wolfe算法相结合的混合算法进行求解,得出最优同步调度方案;Liu等[7]针对因轨道交通网络中断而导致大量乘客滞留的问题,构建了公交接驳轨道交通时刻表优化模型,并提出相关接运方案;Chen等[8]为应对轨道交通中断问题,将快速和短途路线合并,提出时变需求下的公交桥接路线和时刻表综合优化框架.
由于地铁辐射范围越来越广,一些公交车线路因需求量变化导致改道或停运,还有少数公交车线路出于各方面原因与地铁线路保持较高程度的重合. 如果能够充分利用这部分公交车特性,在满足原有公交车出行需求的同时服务有短途出行需求的地铁乘客,无需开设新接驳线路,具有高效率低成本的特性. 可变线路公交(flex-route transit, FRT)[9]的特点恰好可以满足上述需求,在与地铁重合的线路和站点中既能服务原有公交出行需求,又能在固定站点间改变线路,接送有短途出行需求的地铁乘客. 在可变线路公交时刻表编制方面,赵展轮[10]结合可变线路式公交双重特性,构建车辆发车间隔模型和发车时间窗宽度确定方法;姚尹杰[11]从发车时刻等方面对需求响应公交进行系统设计,综合考虑乘客上下车站点时空约束、车辆不可逆行约束等条件,构建了通道型需求响应公交时刻表与路径规划协同优化模型.
综上所述,由于可变线路公交基于乘客预约出行需求进行运行组织,因此运行方案设计需要与站点和线路同步调整. 目前针对可变线路公交运行周期、发车间隔和车辆配置与接驳站点、线路同步优化的研究相对较少,特别是考虑不同时空属性下的灵活出行特征(乘客类型及预约时空维度信息)和接驳站点处时间衔接的影响,因此有待进一步研究和探索. 乘客搭乘地铁和短途出行时会考虑地铁的到站、发车时间,以安排自己的出行计划. 为有效提供接驳服务,需要将可变线路公交与地铁的运营组织方案相配合. 目前,国内大部分城市的公交行车计划仍然依靠人工编制[12],因此需要在考虑运营需求的前提下,提出一种智能化公交行车计划编制方法.
1. 时刻表编制模型与方法
针对接驳地铁站点“最后一公里”出行需求,本研究建立了可变线路公交协同地铁运行时刻表编制模型,编制最优化的可变线路公交时刻表,以高效完成接驳服务,提高公交系统服务水平.
1.1 发车时间间隔
由于可变线路公交具有服务合乘的特点和车辆运行不确定性,既要考虑固定站点的发车间隔和发车时刻,也要考虑站点间灵活接送乘客时的松弛服务时间,因此,传统固定线路公交时刻表编制方式[13-14]和发车频率计算方法[15]并不完全适用于服务灵活出行需求的接驳可变线路公交. 本节提出了一种基于地铁时刻表、乘客预约需求和发车时间窗约束的可变线路公交车接驳地铁发车时间间隔计算方法.
在长(L)*宽(W)服务区域内存在I类乘客:地铁站出站(O)、站外/线外到达需求点(D)和II类乘客:站外/线外出发需求点(O)、地铁站进站(D),如图1所示,现存的常规公交转变成可变线路公交,完成乘客在站点和基准线外之间的需求的接驳服务.
地铁换乘接驳可变线路公交和接驳可变线路公交换乘地铁2种情况,最佳换乘的情况分别为:地铁换乘接驳可变线路公交的乘客抵达公交站台时,可变线路公交车辆恰好在站;接驳可变线路公交换乘地铁的乘客到达地铁站台时,地铁恰好在站,具体换乘过程用时空图表示,如图2所示. 图中:RTa、RTd分别为地铁到站和发车时间;ε为乘客进出地铁站换乘所需时间,min;tp、td分别为上、下可变线路公交时刻;ta,k、tb,k分别为可变线路公交班次k在固定站点的到达、发车时间;Pp,k、Pd,k分别为可变线路公交班次k理想时在固定站点接上、送达乘客时刻;K为接驳可变线路公交车辆发车次数集合,k∈K.
左侧乘客从地铁换乘可变线路公交时,需要在RTa时刻完成下车,经过换乘时间ε在tp时刻(tb,k时刻前)到达接驳站公交点. 右侧乘客自可变线路公交换乘地铁需要在td时刻或之前到接驳站,经过换乘时间ε在RTd时刻或之前到达地铁站点进行换乘.
乘客从地铁换乘可变线路公交时如果过早到达公交站,需要等待一定时间才能被服务;若乘客晚到公交站点,若上一班次被安排服务该乘客,车辆需等待乘客到达站点后接上乘客再发车,这一延误同样也记为等待时间成本. 乘客从可变线路公交换乘地铁时如果早到换乘地铁站则不记为延误;若实际到达时刻晚于预计到达时刻则乘客只能乘坐下一班次地铁,这一晚到延误记为晚到时间成本.
1.1.1 目标函数
接驳地铁可变线路公交时刻表编制模型的目标函数为系统成本U最小化,如式(1)所示.
minU=λ1UW+λ2UL+λ3UF+λ4UO, (1) 式中:系统成本包括乘客/可变线路公交车辆等待时间成本UW、乘客乘坐可变线路公交晚到时间成本UL、换乘失败成本UF和企业运营成本UO,λ1、λ2、λ3和λ4为权重系数.
等待时间涉及车辆发车间隔和排班,故定义为式(2);晚到时间成本根据上述规则定义为式(3);换乘过程中换乘失败成本与换乘失败乘客人数有关,见式(4);企业运营成本取决于车辆班次的运行时间成本,见式(5).
UW=θW∑q∈QI∑k∈Kyq,k|Pp,q,k−tp,q|, (2) UL={0,Pd,q,k⩽td,q,θL∑q∈QII∑k∈Kyq,k(Pd,q,k−td,q),Pd,q,k>td,q, (3) UF=θF(qm−∑k∈Kqk), (4) UO=θO∑α∈S∑β∈S∑k∈Kxαβ,kdαβ,k/v, (5) 式中:θW为地铁换乘可变线路公交的乘客或可变线路公交车辆在公交站的等待时间成本单位值,元/min;θL为可变线路公交换乘地铁的乘客乘坐车辆晚到公交站的时间成本单位值,元/min;θF为换乘失败(乘客未被服务)成本的单位转换值(人数转换为时间成本),元/min;θO为车辆/企业运营成本的单位值,元/min;QI、QII分别为I、II类乘客的集合;Q为所有乘客(需求)集合,Q=QI∪QII,q∈Q,qm为服务区域内总乘客数量;qk为实际乘坐可变线路公交班次k的乘客数量;RC为研究区域内接驳站的集合,RC={r1, r2,…, rn},rn为研究区域内接驳站的数量;RN为接驳可变线路公交车辆需要经过的临时站点(灵活需求点)集合;S为所有站点的集合,S = RC∪RN;Pp,q,k、Pd,q,k分别为可变线路公交班次k接上、送达乘客q的时间;tp,q为地铁换乘可变线路公交的乘客q到达公交车站的时间;td,q为可变线路公交换乘地铁的乘客q预定到达公交车站的时间; dαβ,k为可变线路公交班次k行驶自站点α至站点β的Manhattan距离,km;v为可变线路公交车辆行驶速度,千米/小时;xαβ,k为0-1变量,可变线路公交班次k经过弧(α, β)时取1,否则取0;yq,k为0-1变量,乘客q乘坐可变线路公交班次k时取1,否则取0.
1.1.2 约束条件
式(6)和式(7)表示班次k在运行过程中都必须从可变线路公交首站(第1个接驳站)出发,最后到达末站(最后1个接驳站),同时,站点最多只能被访问一次,即被访问的站点只能有一次到达和一次出发.
∑α∈Sxαβ,k=1,∀β∈S∖{r1},k∈K, (6) ∑β∈Sxαβ,k=1,∀α∈S∖{rn},k∈K. (7) 式(8)表示当xαβ,k = 1时,班次k车辆到达站点β的时间将不早于从站点α出发的时间加上两站之间的行驶时间,当xαβ,k = 0时,此约束无效.
ta,β,k⩾tb,α,k+xαβ,kdαβ,k/v−M(1−xαβ,k),∀(α,β)∈S,k∈K, (8) 式中:ta,β,k 、tb,α,k分别为可变线路公交班次k在站点β和α的到达、发车时间,M为一个足够大的数.
式(9)确保班次k在非接驳站的发车时间等于车辆到达非接驳站的时间加上停车服务时间. 类似地,式(10)确保班次k在接驳站的发车时间等于车辆在接驳站的到达时间加上服务时间.
tb,α,k=ta,α,k+trn,∀α∈RN,k∈K, (9) tb,α,k=ta,α,k+tr,∀α∈RC∖{r1,rn},k∈K, (10) 式中:trn为可变线路公交在临时站点(灵活需求点)的服务时间,tr为在接驳站的停靠时间.
式(11)表示消除使其他点不相交的所有循环,确保在最优解中不包含无意义的子循环
∑α∈ψ∑α∈ψxαβ,k⩽|ψ|−1,∀ψ⊆S,ψ≠∅,k∈K, (11) 式中:ψ为集合S除空集外的真子集.
式(12)表示每个乘客最多被服务一次.
∑k∈Kyq,k⩽1,∀q∈Q. (12) 式(13)和(14)用于判断换乘乘客能否顺利完成换乘,换乘成功(使用接驳可变线路公交)yq,k则取1,反之为0.
M(yq,k−1)⩽|Pp,q,k−tp,q|<Myq,k,∀q∈QI,k∈K, (13) M(yq,k−1)⩽|Pd,q,k−td,q|<Myq,k,∀q∈QII,k∈K. (14) 当换乘乘客q乘坐可变线路公交班次k的等待时间或晚到时间大于τ时为1,否则为0,如式(15)、(16)所示.
M(zq,k−1)⩽|Pp,q,k−tp,q|−τ<Mzq,k,∀q∈QI,k∈K, (15) M(zq,k−1)⩽|Pd,q,k−td,q|−τ<Mzq,k,∀q∈QII,k∈K, (16) 式中:τ为乘客最大可容忍等待或晚到时间, min;zq,k为0-1变量.
式(17)、(18)表示乘客下车时刻要大于上车时刻加上接驳站点或临时站点的服务时间.
Pp,q,k+tr<Pd,q,k,∀q∈QI,k∈K, (17) Pp,q,k+trn<Pd,q,k,∀q∈QII,k∈K. (18) 式(19)表示可变线路公交运行中的发车间隔受时间区间约束.
Hmin⩽H⩽Hmax,∀H∈N+, (19) 式中:Hmin、Hmax分别为车辆发车间隔下、上限,min;H为可变线路公交运营周期内发车间隔,min.
运营周期内发车班次次数由下式(20)得到.
km=hH,∀h∈N+,H|h (20) 式中:h为可变线路公交运营周期,分钟;km为运营周期内发车班次次数,次.
通过相关标准及经验确定发车间隔的上下限以及潜在的最优发车间隔(如:5 min、10 min),再根据上述提出的模型依次计算和对比不同发车间隔的各项系统成本,从而确定出最优发车间隔和发车班次.
1.2 站点发车时间
由于需求的不确定性和不均匀性,分配和制定站点之间的松弛时间(制定车辆运行时刻表时额外分配给服务灵活需求点的时间,其值等于系统所分配的运行周期时间减去车辆沿基准路线的直接行驶时间)和发车时间并不是简单的求平均值,要根据实际需求来进行合理的规划. 对于每个班次车辆来说,车辆到站和发车时间受各个接驳站之间的松弛时间影响,根据车辆行驶时间、车辆在每个临时站点和每个接驳站的服务时间可计算得到车辆每班次单程的运行周期时间,见式(21).
ςk⩾∑α∈S∑β∈Sxαβ,k×dαβ,kv+qi,ktrn+(rm−2)tr,∀i∈RC∖{rn},k∈K, (21) 式中:ςk为可变线路公交班次k运行周期,min;qi,k为可变线路公交班次k在接驳站ri与ri + 1之间的临时站点上、下车乘客的数量.
式(22)表示运行周期减去车辆沿基准线行驶的时间即为总松弛时间;式(23)表示各个接驳站区间的松弛时间之和,即各班次车辆完整运行的总松弛时间. 各个区间内的松弛时间主要由区间内车辆所需服务的临时站点数来分配和确定,不再受不同区间的服务区域大小的影响,如式(24).
Lk=ςk−∑i∈RC∖{rn}dr,i,k/v,∀k∈K, (22) Lk=∑i∈RC∖{rn}Li,k,∀k∈K, (23) Li,k=qi,k∑i∈RC∖{rn}qi,kLk,∀i∈RC∖{rn},k∈K, (24) 式中:Lk为可变线路公交班次k运行周期内总的松弛时间,min;Li,k为可变线路公交班次k在接驳站ri与ri + 1之间分配的松弛时间,min;dr,i,k为班次k沿基准线路行驶时在接驳站ri与ri + 1之间的行驶距离,km.
各个接驳站的发车时间通常由前一个接驳站的发车时间和2个接驳站间的行程时间和松弛时间共同决定,见式(25).
ξi,k=ξi−1,k+Li−1,k+dr,i−1,kv,∀i∈RC∖{r1,rn},k∈K, (25) 式中:ξi,k为可变线路公交班次k在接驳站ri的发车时间.
2. 算法设计
在编制接驳可变线路公交时刻表时,可以将问题细分为“发车间隔设置算法设计”“班次设置算法设计”和“发车时间设置算法设计”3个子问题.
2.1 发车间隔设置阶段
设置发车间隔时需要考虑公交时刻表的特性,制定合理高效的方案. 具体求解算法伪代码如表1所示.
表 1 发车间隔设置算法设计Table 1. Algorithm design of departure interval setting步骤 指令 输入 发车间隔下限Hmin,发车间隔上限Hmax,发车间隔Hν,其中:Hν = Hmin,⋯, Hmax,ν = 1, 2, ⋯, n 步骤 取初始H1 = Hmin,令H* = H1,计算当前U(H1),令U*(H*) = U(H1) 循环 当ν < n时,执行步骤1、2 步骤1 在上述基础上,获取新的Hν + 1,判断是否满足约束,若Hν + 1不满足约束,则重复步骤1 步骤2 计算Hν + 1下的U(Hν + 1),若U(Hν + 1) < U*(H*),则接受Hν + 1作为新的最优发车间隔,H* = Hν + 1,否则保留原来的H* 步骤3 若满足终止条件Hn = Hmax,则输出H*和U*(H*),结束程序,否则返回步骤1. 输出 最优发车间隔H*,最优系统成本U*(H*) 2.2 班次设置阶段
每个不同的发车间隔对应了不同的乘客分配和可变线路公交车辆班次设置,通过对分配和设置后的具体接驳情况和班次安排计算得到系统成本,因此,这一部分也需要设计相应算法进行求解. 乘客分配和班次设置决定每个班次具体服务哪些乘客,可以运用遗传算法求解,在乘客分配和车辆班次设置的时候,尝试将被拒绝的乘客重新安排回运行班次,从而提高解的质量,具体求解算法伪代码如表2所示.
表 2 班次设置算法设计Table 2. Algorithm design of shift setting步骤 指令 输入 种群规模,交叉概率,变异概率,最大迭代数Imax 步骤0 输入运营参数,设置种群数量、交叉、变异概率,最大迭代次数Imax,初始化种群迭代次数In = 0 循环 当In < Imax时,执行步骤1~8 步骤1 编码:染色体采用由0-1变量构成的矩阵体现乘客分配和班次设置的二维关系,并假设存在一个新的虚拟班次存放被拒绝的乘客 步骤2 适应度函数:当完成运营周期内一次完整班次设置之后,计算染色体对应的适应度 步骤3 染色体更新:1) 在每个乘客和班次完成分配和设置后,将被拒绝乘客的信息反馈给染色体,对相应的基因进行更改,实现染色体的更新;2) 尝试将被拒绝的乘客重新安排回接驳可变线路公交中,提高染色体质量,跳出局部最优 步骤4 选择:精英保留策略和锦标赛法 步骤5 交叉:选用单点交叉方式 步骤6 变异:选用单点变异方式 步骤7 种群适应度更新并记录最优解 步骤8 迭代次数更新,In = In + 1 步骤9 算法终止条件:In ≥ Imax 输出 最优解u,最优系统成本U 2.3 发车时间设置阶段
乘客分配和班次设置后需要对各班次的路径进行规划可求出系统成本U,同时路径规划后才能求解得到接驳可变线路公交系统的运行周期时间、松弛时间和站点发车时间. 站点发车时间主要指固定站点的发车时间,为避免过高的拒绝率或过长的延误时间,需要考虑整体运行周期,对每个运行区间分配合适的松弛时间. 具体求解算法如表3所示.
表 3 发车时间设置算法设计Table 3. Algorithm design of departure time setting步骤 指令 输入 初始运行站点、线路,出行需求信息,乘客分配情况,车辆班次情况,初始车辆状态 步骤0 输入出行需求行程时间分布和接驳可变线路公交线路运营各项参数 步骤1 将乘客起终点在各分配的车辆班次中依次插入 步骤2 判断约束,若乘客起点或终点的插入会导致可变线路公交违背接驳站的时间约束,则拒绝请求 步骤3 被拒绝乘客重新安排,明确班次的选择顺序,分别将乘客实际上、下车时间与预计/预约上、下车时间比较,时间成本越小,班次排序越靠前,乘客优先尝试插入到该班次中,不再尝试安排到其他班次; 步骤4 计算得到各班次的运行周期时间,松弛时间和固定站点发车时间; 步骤5 编制接驳地铁的可变线路公交时刻表. 输出 最优各班次车辆运行路径,运行周期时间,松弛时间,接驳站发车时间 3. 案例研究
3.1 案例描述
本文选取南京市地铁2号线和公交5路作为研究对象,根据流量、线路、时刻表等运营数据验证模型及算法,布局走向如图3所示. 图中蓝色实线为南京市地铁2号线,红色实线为南京市公交5路,位于鼓楼区、秦淮区和栖霞区的交界处,承担着通勤、科研教育、旅游休闲等出行需求.
线路周围POI点非常多,如图4所示,图中紫色小点为各类型POI点,吸引着大量乘客出行往返,地铁和公交承担着高水平的出行需求服务. 乘客出行乘坐地铁后大概率要到达某个POI点,同样地,乘客乘坐地铁前会从某个POI点出发. 本文主要解决地铁乘客“最后一公里”出行问题,因此筛选得到线路和站点附近的POI点(潜在出行需求点),在图中显示,下述出行需求集合会在这部分POI点中随机生成,同时考虑线路和车辆可达性,形成可供可变线路公交接驳的潜在出行需求集合.
公交车辆行驶速度和乘客出行需求水平在一天中是分时段变化的,具体随时间段变化情况如图5所示. 在高峰期时,出行需求水平较高,同时由于道路拥堵严重,公交车辆行驶速度较低;在平峰时段,出行需求较少,道路较为通畅,公交车辆行驶速度随之提高. 各时段出行需求从乘客出行需求集合中随机选取,用作系统仿真分析数据基础.
南京公交5路站点和线路与地铁2号线新街口-马群段高度重合,可将5路常规公交转换为可变线路公交,在保证原有公交定点定线乘客出行需求的同时接驳服务地铁乘客短途出行需求. 将重合度较高的地铁-公交站点作为接驳可变线路公交的固定站点,包括新街口-新街口东站、大行宫-大行宫东站、西安门•浦发银行-西安门站、明故宫-明故宫东站、下马坊-小卫街站、孝陵卫-孝陵卫站、钟灵街-钟灵街站、马群-马群站. 5路公交转变为可变线路公交后,无站外需求时按照原始线路(基准线)运行,如有站外需求点,包括站外到地铁乘车或地铁站下车到站外需求点,可变线路公交可以偏离基准线到站外需求点接、送乘客,方便各类型乘客的出行.
对系统进行仿真分析时,由于系统道路沿线路网较为复杂,因此根据前面章节的相关描述,将系统抽象为一个LW的服务区域,可变线路公交在服务区域内接驳乘客出行. L为接驳站两端间的路线长度,考虑到“最后一公里”短途出行的目的,将服务区域宽度W设为1公里. 接驳站相应位置信息如表4所示.
表 4 固定站点地铁运行时刻表Table 4. Urban rail transit running timetable data of fixed stations编号 接驳站
(地铁-公交)固定站点
相对位置始发
时间末班
时间r1 新街口(东) 0 06:11 23:33 r2 大行宫(东) 1.00 06:13 23:35 r3 西安门•浦发银行 2.00 06:15 23:37 r4 明故宫(东) 3.40 06:17 23:40 r5 下马坊-小卫街 6.20 06:21 23:44 r6 孝陵卫 7.20 06:23 23:46 r7 钟灵街 8.30 06:25 23:48 r8 马群 11.00 06:00 23:52 接驳地铁的可变线路公交需要配合地铁运行时刻表才能更好地完成乘客接送服务,地铁到站/发车时刻表目前有始发时间和末班时间数据,如表4所示,具体运行周期(如1天)内的详细时刻表未被公布. 因此本文根据现有信息和相关运行规则,合理推算服务周期内的地铁到/发站时刻表,并依此仿真生成短途出行需求的时空信息,接驳可变线路公交根据出行需求时空信息合理编制运行时刻表,制定切实可行的发车和运行计划.
研究系统中相关参数参考[6,13-15]. θW、θL和θO都设为1 元/min,考虑到公共交通的公益性,尽可能地满足出行需求,减少拒绝率,因此将θF设为15 元/min[15],Hmin为3 min, Hmax为30 min,ε为3 min,τ为5 min,tr、trn分别为1 min和0.3 min. 与相关研究[14,15]类似,应用直线运动来再现真实的道路网络.
本节基于Matlab R2021a平台,通过编程对模型进行求解. 对于每一次的仿真模拟实验,遗传算法的参数设置如下:染色体种群规模为80,发生交叉行为的概率为0.9,发生变异的概率为0.1,最大迭代次数为500次.
3.2 案例分析
根据系统相关参数值、车辆行程速度和乘客各时段需求时空数据,可得到不同需求水平下不同车速所对应的最优发车间隔,如表5所示. 这里的车辆发车间隔主要是针对转变后的可变线路公交车,并非原始5路公交车全部按照这一发车间隔和方案运营.
表 5 不同需求水平下不同车速所对应的最优发车间隔Table 5. Optimal departure interval corresponding to different driving speeds under different demand levels需求/(人·h−1) 发车间隔H/(分钟) v=20千米/
小时v=30千米/
小时v=40千米/
小时10 30 30 30 20 30 30 30 30 20 20 30 40 15 20 20 50 12 15 15 60 10 12 15 根据系统中车辆行程时间和乘客需求随时间的变化情况,可得到运营周期内不同时段下接驳地铁可变线路公交发车时间间隔的变化情况,如表6所示.
表 6 不同时段下接驳地铁可变线路公交发车时间间隔Table 6. Departure time interval of feeder FRT at different periods时段 V/(km·h−1) 需求/(人·h−1) H /min 7:00-8:59 20 60 10 9:00-10:59 40 40 20 11:00-12:59 30 50 15 13:00-14:59 30 50 15 15:00-16:59 40 30 30 17:00-18:59 20 60 10 19:00-20:59 30 20 30 在计算系统成本时,默认各项时间成本的权重系数都为1,若重视乘客等待和晚到时间成本,增加其相应的权重,如设置等待时间成本、晚到时间成本、换乘失败成本和企业运营成本的权重系数比为2∶2∶1∶1,最优发车间隔随之发生变化,如图6所示.
从图6可以看出,若重视乘客等待和晚到时间成本,相应权重系数增大后,最优发车间隔则发生改变,需要减小发车间隔,增加发车频率,降低乘客相关出行时间成本的同时得到最优系统运行成本.
根据表6,在早晚高峰时段(7:00-8:59、17:00-18:59),出行需求较大的同时车辆行驶速度缓慢,因此,发车时间间隔相应缩短;当出行需求较低车速较快时,发车间隔时间可以相应增大,也能满足乘客出行需求. 同时,计算得到各运行周期内的系统运行成本U、各运行周期内的平均车辆运行时间ςk、各运行周期内的发车次数和所需车辆数,如表7所示. 其中,所需车辆数是根据每趟班车单程运行时间确定的,因本文主要对接驳系统和方法进行研究和验证,计算单程时刻表即可,假设回程所需时间与去程相等,从而计算各运行周期内所需总车辆数.
表 7 系统各运行时段的仿真结果Table 7. Simulation results of each operating period of system时段 H/min U/元 平均ςk/min 发车次数
/次车辆数
/辆7:00-8:59 10 801.22 53.71 12 11 9:00-10:59 20 279.18 33.73 6 4 11:00-12:59 15 430.76 42.17 8 6 13:00-14:59 15 493.44 43.48 8 6 15:00-16:59 30 192.75 35.12 4 3 17:00-18:59 10 895.30 57.99 12 11 19:00-20:59 30 210.46 41.27 4 3 从表7中可以发现,车辆在早晚高峰时段运行因受制于行程速度限制以及高出行需求水平,其单程运行周期较长,导致在这2个时段内的运营费用相对较高. 为满足出行需求和发车时间间隔要求,需要在这2个时段内配置更多的车辆. 相比之下,在平峰时段,随着需求量减少和车辆行驶速度增加,较少车辆即可满足出行需求,系统总成本相应降低. 因此,在实际运营中,应根据不同时段行程速度和乘客需求的变化,设定不同的发车时间间隔,不仅有助于满足不同时段乘客出行需求,还能有效控制接驳地铁可变线路公交系统运营成本,从而实现更好的运营效果.
若重视乘客等待和晚到时间成本,运行周期内的平均车辆运行时间和系统运行成本也会随之发生变化,如图7、8所示.
从图7可以看出,重视乘客出行时间则系统发车间隔减小,运行周期内车辆数随之增加,使得各个班次车辆接送乘客平均数减少,对应的平均运行距离减少,因此,运行周期内平均车辆运行时间相应减少.
图8中常规系统运行成本指默认权重值情况下系统的总运行成本,与表7一致;常规发车间隔重视乘客系统运行成本指重视乘客时间成本情况下原发车间隔时的系统运行成本;重视乘客系统运行成本指发车间隔和平均车辆运行时间相应减少时的系统运行成本. 可以看出,当乘客权重系数增大时常规情况下的系统运行成本随之增大,此时降低发车间隔,增大发车频率,增加运行车辆数则可以更大程度上降低乘客等待和晚到时间,使系统总成本相应减少. 由于车辆数增加,相应的运营成本也随着增高,所以系统运行成本不会有大幅度的降低.
4. 结 论
1) 分析了部分使用地铁出行乘客的“最后一公里”出行需求特征和换乘行为特征,分析了公交车和地铁的运行组织特征和衔接特征,提出一种将常规公交转变为可变线路公交接驳地铁的策略,在不需要开设新接驳线路的前提下可以同时服务原有公交车出行需求乘客和有短途出行需求的地铁乘客.
2) 针对可变线路公交区别于常规定点定线公交的运行特征,以公共交通运行数据和服务范围内POI数据为基础,综合考虑了乘客出行成本和公交系统车辆运营成本,构建可变线路公交运行时刻表编制方法和接驳地铁的路径规划与运行时刻表编制协同优化模型与算法.
3) 以实际路网中地铁和公交系统为例,分析研究了车辆运行速度与出行需求水平对系统各项指标的影响,并分析不同时间成本侧重情况下对系统发车间隔、车辆运行时间、运行成本等指标的影响,结果表明在运行周期内适当增加发车频率可以有效降低乘客出行时间成本同时不会过多增加系统运营成本,并编制车辆行车时刻表. 可为企业实际运营提供一定的参考,对具有短途出行需求的地铁乘客经济、便捷出行和公交系统低成本、高效性、公益性具有一定的促进意义.
4) 在本文给出的时刻表构建方法和案例分析的基础上,进一步扩大协同服务网络规模、引入不同车辆车型组合是开展后续相关研究的重要方向.
致谢:合肥大学城市轨道交通安全与应急管理重点实验室开放课题(2024GD0007).
-
表 1 发车间隔设置算法设计
Table 1. Algorithm design of departure interval setting
步骤 指令 输入 发车间隔下限Hmin,发车间隔上限Hmax,发车间隔Hν,其中:Hν = Hmin,⋯, Hmax,ν = 1, 2, ⋯, n 步骤 取初始H1 = Hmin,令H* = H1,计算当前U(H1),令U*(H*) = U(H1) 循环 当ν < n时,执行步骤1、2 步骤1 在上述基础上,获取新的Hν + 1,判断是否满足约束,若Hν + 1不满足约束,则重复步骤1 步骤2 计算Hν + 1下的U(Hν + 1),若U(Hν + 1) < U*(H*),则接受Hν + 1作为新的最优发车间隔,H* = Hν + 1,否则保留原来的H* 步骤3 若满足终止条件Hn = Hmax,则输出H*和U*(H*),结束程序,否则返回步骤1. 输出 最优发车间隔H*,最优系统成本U*(H*) 表 2 班次设置算法设计
Table 2. Algorithm design of shift setting
步骤 指令 输入 种群规模,交叉概率,变异概率,最大迭代数Imax 步骤0 输入运营参数,设置种群数量、交叉、变异概率,最大迭代次数Imax,初始化种群迭代次数In = 0 循环 当In < Imax时,执行步骤1~8 步骤1 编码:染色体采用由0-1变量构成的矩阵体现乘客分配和班次设置的二维关系,并假设存在一个新的虚拟班次存放被拒绝的乘客 步骤2 适应度函数:当完成运营周期内一次完整班次设置之后,计算染色体对应的适应度 步骤3 染色体更新:1) 在每个乘客和班次完成分配和设置后,将被拒绝乘客的信息反馈给染色体,对相应的基因进行更改,实现染色体的更新;2) 尝试将被拒绝的乘客重新安排回接驳可变线路公交中,提高染色体质量,跳出局部最优 步骤4 选择:精英保留策略和锦标赛法 步骤5 交叉:选用单点交叉方式 步骤6 变异:选用单点变异方式 步骤7 种群适应度更新并记录最优解 步骤8 迭代次数更新,In = In + 1 步骤9 算法终止条件:In ≥ Imax 输出 最优解u,最优系统成本U 表 3 发车时间设置算法设计
Table 3. Algorithm design of departure time setting
步骤 指令 输入 初始运行站点、线路,出行需求信息,乘客分配情况,车辆班次情况,初始车辆状态 步骤0 输入出行需求行程时间分布和接驳可变线路公交线路运营各项参数 步骤1 将乘客起终点在各分配的车辆班次中依次插入 步骤2 判断约束,若乘客起点或终点的插入会导致可变线路公交违背接驳站的时间约束,则拒绝请求 步骤3 被拒绝乘客重新安排,明确班次的选择顺序,分别将乘客实际上、下车时间与预计/预约上、下车时间比较,时间成本越小,班次排序越靠前,乘客优先尝试插入到该班次中,不再尝试安排到其他班次; 步骤4 计算得到各班次的运行周期时间,松弛时间和固定站点发车时间; 步骤5 编制接驳地铁的可变线路公交时刻表. 输出 最优各班次车辆运行路径,运行周期时间,松弛时间,接驳站发车时间 表 4 固定站点地铁运行时刻表
Table 4. Urban rail transit running timetable data of fixed stations
编号 接驳站
(地铁-公交)固定站点
相对位置始发
时间末班
时间r1 新街口(东) 0 06:11 23:33 r2 大行宫(东) 1.00 06:13 23:35 r3 西安门•浦发银行 2.00 06:15 23:37 r4 明故宫(东) 3.40 06:17 23:40 r5 下马坊-小卫街 6.20 06:21 23:44 r6 孝陵卫 7.20 06:23 23:46 r7 钟灵街 8.30 06:25 23:48 r8 马群 11.00 06:00 23:52 表 5 不同需求水平下不同车速所对应的最优发车间隔
Table 5. Optimal departure interval corresponding to different driving speeds under different demand levels
需求/(人·h−1) 发车间隔H/(分钟) v=20千米/
小时v=30千米/
小时v=40千米/
小时10 30 30 30 20 30 30 30 30 20 20 30 40 15 20 20 50 12 15 15 60 10 12 15 表 6 不同时段下接驳地铁可变线路公交发车时间间隔
Table 6. Departure time interval of feeder FRT at different periods
时段 V/(km·h−1) 需求/(人·h−1) H /min 7:00-8:59 20 60 10 9:00-10:59 40 40 20 11:00-12:59 30 50 15 13:00-14:59 30 50 15 15:00-16:59 40 30 30 17:00-18:59 20 60 10 19:00-20:59 30 20 30 表 7 系统各运行时段的仿真结果
Table 7. Simulation results of each operating period of system
时段 H/min U/元 平均ςk/min 发车次数
/次车辆数
/辆7:00-8:59 10 801.22 53.71 12 11 9:00-10:59 20 279.18 33.73 6 4 11:00-12:59 15 430.76 42.17 8 6 13:00-14:59 15 493.44 43.48 8 6 15:00-16:59 30 192.75 35.12 4 3 17:00-18:59 10 895.30 57.99 12 11 19:00-20:59 30 210.46 41.27 4 3 -
[1] 夏亮,江欣国,范英飞. 基于均值-方差理论的多阶段公交走廊设计[J]. 西南交通大学学报,2023,58(6): 1294-1302. doi: 10.3969/j.issn.0258-2724.20210874XIA Liang, JIANG Xinguo, FAN Yingfei. Phased transit service design based on mean-variance theory[J]. Journal of Southwest Jiaotong University, 2023, 58(6): 1294-1302. doi: 10.3969/j.issn.0258-2724.20210874 [2] 蒋阳升,李衍,李皓,等. 基于模块化仿真的共享汽车联合调度优化[J]. 西南交通大学学报,2023,58(1): 74-82. doi: 10.3969/j.issn.0258-2724.20210083JIANG Yangsheng, LI Yan, LI Hao, et al. Optimization for joint relocation of carsharing based on modular simulation[J]. Journal of Southwest Jiaotong University, 2023, 58(1): 74-82. doi: 10.3969/j.issn.0258-2724.20210083 [3] 胡路,梁志梅,蒋阳升. 拥挤传播对汽车共享系统运行影响的仿真分析[J]. 西南交通大学学报,2023,58(3): 499-510. doi: 10.3969/j.issn.0258-2724.20220231HU Lu, LIANG Zhimei, JIANG Yangsheng. Simulation analysis on influence of congestion propagation on operation of carsharing systems[J]. Journal of Southwest Jiaotong University, 2023, 58(3): 499-510. doi: 10.3969/j.issn.0258-2724.20220231 [4] WANG J D, YUAN Z Z, CAO Z C, et al. Optimal bus bridging schedule with transfer passenger demand during disruptions of urban rail transit[J]. Journal of Transportation Engineering, Part A: Systems, 2021, 147(10): 04021071.1-04021071.12. [5] 窦雪萍,董冉,李同飞,等. 考虑乘客排队行为的接驳公交行车计划优化[J]. 系统科学与数学,2022,42(12): 3355-3367. doi: 10.12341/jssms22368DOU Xueping, DONG Ran, LI Tongfei, et al. Feeder bus schedule optimization considering passenger queuing behavior[J]. Journal of Systems Science and Mathematical Sciences, 2022, 42(12): 3355-3367. doi: 10.12341/jssms22368 [6] XIONG J, HE Z B, GUAN W, et al. Optimal timetable development for community shuttle network with metro stations[J]. Transportation Research Part C: Emerging Technologies, 2015, 60: 540-565. doi: 10.1016/j.trc.2015.10.007 [7] LIU J R, WU X Y. Research of emergency feeder bus scheme for urban rail transit under unexpected incident[J]. IOP Conference Series: Materials Science and Engineering, 2019, 688(4): 044040.1-044040.7. [8] CHEN Y, AN K. Integrated optimization of bus bridging routes and timetables for rail disruptions[J]. European Journal of Operational Research, 2021, 295(2): 484-498. doi: 10.1016/j.ejor.2021.03.014 [9] 邱丰,李文权,沈金星. 可变线路式公交的两阶段车辆调度模型[J]. 东南大学学报(自然科学版),2014(5): 1078-1084.QIU Feng, LI Wenquan, SHEN Jinxing. Two-stage model for flex-route transit scheduling[J]. Journal of Southeast University (Natural Science Edition), 2014(5): 1078-1084. [10] 赵展轮. 可变线路公交运行时刻表编制方法研究[D]. 南京:东南大学,2019. [11] 姚尹杰. 通道型需求响应公交的调度优化研究[D]. 广州:华南理工大学,2021. [12] SHANG H Y, HUANG H J, WU W X. Bus timetabling considering passenger satisfaction: an empirical study in Beijing[J]. Computers & Industrial Engineering, 2019, 135: 1155-1166. [13] Associates I. TRB’s Transit Cooperative Research Program (TCRP) report 100: transit capacity and quality of service manual, second edition[M]. Washington, D. C. : TRB, 2003. [14] 美国交通运输研究委员会. 公共交通通行能力和服务质量手册[M]. 北京:中国建筑工业出版社,2010. [15] Avishai C. 公共交通规划与运营──理论、建模及应用[M]. 北京:清华大学出版社,2010. [16] ZHENG Y, LI W Q, QIU F. A slack arrival strategy to promote flex-route transit services[J]. Transportation Research Part C: Emerging Technologies, 2018, 92: 442-455. doi: 10.1016/j.trc.2018.05.015 [17] 茧敏,王卓,陈哲轩,等. 基于时空路径推算的城轨网络客流分配方法[J]. 西南交通大学学报,2023,58(5): 1117-1125. doi: 10.3969/j.issn.0258-2724.20220545JIAN Min, WANG Zhuo, CHEN Zhexuan, et al. Passenger flow assignment method for urban rail transit networks based on inference of spatiotemporal path[J]. Journal of Southwest Jiaotong University, 2023, 58(5): 1117-1125. doi: 10.3969/j.issn.0258-2724.20220545 [18] QUADRIFOGLIO L, DESSOUKY M M, PALMER K. An insertion heuristic for scheduling Mobility Allowance Shuttle Transit (MAST) services[J]. Journal of Scheduling, 2007, 10(1): 25-40. doi: 10.1007/s10951-006-0324-6 [19] ZHENG Y, LI W Q, QIU F, et al. The benefits of introducing meeting points into flex-route transit services[J]. Transportation Research Part C: Emerging Technologies, 2019, 106: 98-112. doi: 10.1016/j.trc.2019.07.012 -