Method of Determining Design Value of Rockfall Block Size Based on the Return Period
-
摘要:
山区交通落石事件频发且具有突发性和随机性,基于概率的防护结构可靠性设计是降低灾害风险的必要和有效手段,然而,作为设计基本参数的落石尺寸设计值的确定方法还不明确,使得可靠性设计不能有效进行. 介绍基于泊松-广义帕累托分布(GPD)复合模型的落石尺寸-重现期预测方法,结合宝成铁路白—上段落石记录数据,进行GPD阈值选取,并分析记录中缺失的最大、最小落石尺寸对阈值选取和对预测结果的影响;与地震作用标准值超越概率对应,提出以“小落石不坏、中落石可修、大落石不塌”三水准设防为标准的落石尺寸设计值确定方法. 研究表明:记录中缺失的最大落石尺寸对阈值选取影响小,但对落石尺寸预测结果影响大,与重现期对应的落石尺寸预测值随设定的最大落石尺寸呈指数型增大趋势,而记录的最小落石尺寸对阈值选取和预测结果的影响都较小;以宝成铁路白—上段落石记录数据为例,根据选取出的合理阈值,设定不同最大落石尺寸值,利用本文方法,预测得到宝成铁路白—上段50、100、
1000 、10000 a重现期落石尺寸区间范围,计算得到50、100 a设计基准期三水准设防标准对应的落石尺寸标准值即设计值,可为类似工程设计提供参考.Abstract:The rockfall events on transport lines in mountainous areas occur frequently and are sudden and random. Probability-based reliability design of protective structures is a necessary and effective means to reduce the hazard risk. However, the method of determining the design value of rockfall size, as a basic parameter of the design, is still unclear, making the reliability design fail to be carried out effectively. The prediction method of rockfall size and the return period was introduced based on the Poisson-generalized Pareto distribution (GPD) composite model. According to the rockfall record data of the Bai−Upper section of Baoji−Chengdu Railway, the GPD threshold was selected, and the impact of the missing maximum and minimum rockfall sizes (vmax and vmin) in the record on the threshold selection and prediction results was analyzed. Corresponding to the exceedance probability of the standard value of earthquake action, a method of determining the design value of rockfall size was proposed, which was based on the three-level protection standard featuring “the structure doesn’t break under small rockfalls, can be repaired under medium rockfalls, and doesn’t collapse under large rockfalls”. The results show that the missing maximum rockfall size in the record has a small impact on the threshold selection but a large impact on the prediction results of the rockfall size. The prediction value of rockfall size corresponding to the return period increases exponentially with the set maximum rockfall size, while the recorded minimum rockfall size has less impact on both the threshold selection and the prediction results. By taking the record data of the Bai−Upper section of Baoji−Chengdu Railway as an example, according to the selected reasonable threshold, different maximum rockfall size values are set. By using the method presented in this paper, the interval range of rockfall size in return periods of 50, 100, 1 000, and 10 000 years on the Bai−Upper section of Baoji−Chengdu Railway is predicted, and the standard value (design value) of rockfall size corresponding to the three-level protection standard in design reference periods of 50 and 100 years is calculated, which can provide a reference for similar engineering designs.
-
Key words:
- protective structure /
- rockfall size /
- return period /
- generalized Pareto distribution /
- design value
-
随着出行需求的不断增长和出行选择日渐向快捷、舒适方向发展,私家车或网约车使用率提高,造成交通拥堵[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,$\cdots $, Hmax,ν = 1, 2, $\cdots $, 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).
-
年份 v>3.0 m3 2.0 m3<v≤3.0 m3 1.5 m3<v≤2.0 m3 1.0 m3<v≤1.5 m3 0.5 m3<v≤1.0 m3 v≤0.5 m3 合计 1961 年 1 1 0 10 12 1962 年 4 1 2 2 17 26 1963 年 4 2 1 4 30 41 1964 年 3 4 2 3 5 54 71 1965 年 1 2 2 45 50 1966 年 5 2 3 12 99 121 1967 年 1 1 1 5 8 1969 年 3 1 5 26 35 1970 年 1 1 5 16 51 74 1971 年 1 1 2 3 3 30 40 1972 年 3 2 3 3 2 39 52 合计 23 15 10 24 52 406 530 分占比例% 4.34 2.83 1.89 4.53 9.81 76.6 100.00 表 2 GPD参数σ、ξ估计
Table 2. Estimated results of GPD parameters ξ and σ
条件 ξ σ 90% 置信区间 [0.40,0.60] [0.86,1.02] 95% 置信区间 [0.38,0.63] [0.84,1.04] 均值 0.50 0.94 表 3 各重现期对应落石尺寸均值及其置信区间
Table 3. Mean rockfall sizes and their confidence intervals during different return periods
m3 重现期/a 均值 90% 置信区间 95% 置信区间 vmax=4.0 m3 vmax=6.0 m3 vmax=4.0 m3 vmax=6.0 m3 vmax=4.0 m3 vmax=6.0 m3 50 4.93 13.42 [4.28,5.89] [11.30,14.92] [4.13, 6.12] [11.03, 15.26] 100 5.09 15.91 [4.39,6.12] [13.34,17.86] [4.23, 6.38] [12.89,18.72] 1000 5.41 26.75 [4.60,6.69] [21.23,31.83] [4.42, 7.01] [20.28,33.38] 10000 5.58 42.99 [4.70,7.02] [31.95,53.43] [4.51, 7.39] [30.21,56.86] 表 4 各设计基准期落石尺寸标准值
Table 4. Standard values of rockfall size in different design reference periods
m3 设计基准期/a 设防
水准T/a 1/T/a− 1 1−1/T N 年不超越
概率
((1− 1/T)N)N 年超越
概率
(1−(1−1/T)N)与 T 对应的
落石尺寸标准值vmax=4.0 m3 vmax=6.0 m3 50 Ⅰ 50 0.0200 0.9800 0.36 0.64 4.93 13.4 Ⅱ 475 0.0021 0.9979 0.90 0.10 5.33 22.8 Ⅲ 2475 0.0004 0.9996 0.98 0.02 5.49 32.4 100 Ⅰ 100 0.0100 0.9900 0.37 0.63 5.09 15.9 Ⅱ 900 0.0011 0.9989 0.89 0.11 5.40 26.2 Ⅲ 5000 0.0002 0.9998 0.98 0.02 5.54 37.4 -
[1] MARCHELLI M, COLTRINARI G, ALFARO DEGAN G, et al. Towards a procedure to manage safety on construction sites of rockfall protective measures[J]. Safety Science, 2023, 168: 106307.1-106307.8. [2] 中华人民共和国住房和城乡建设部. 铁路工程结构可靠性设计统一标准:GB 50216—2019[S]. 北京:中国计划出版社,2019. [3] 中华人民共和国交通运输部. 公路工程结构可靠性统一标准:JTG 2120—2020[S]. 北京:人民交通出版社,2020. [4] 国家铁路局. 铁路隧道设计规范:TB 10003—2016[S]. 北京:中国铁道出版社, 2017. [5] 中华人民共和国交通运输部. 公路隧道设计规范 第一册 土建工程:JTG 3370.1—2018[S]. 北京:人民交通出版社, 2019. [6] 中华人民共和国住房和城乡建设部. 建筑结构可靠性设计统一标准:GB 50068—2018[S]. 北京:中国建筑工业出版社. [7] 铁道第二勘测设计院. 铁路工程设计技术手册. 隧道[M]. 北京:中国铁道出版社,1999. [8] 曾廉. 明洞顶设计荷载的研究[J]. 铁路标准设计通讯,1974(7): 3-17. [9] ILLEDITSCH M, PREH A. Determination of meaningful block sizes for rockfall modelling[J]. Natural Hazards, 2024, 120(6): 5685-5710. doi: 10.1007/s11069-024-06432-4 [10] STELZER G, BICHLER A. ONR 24810–A comprehensive guideline for building better rockfall protection structures[C]//64th Annual Highway Geology Symposium. New Hampshire: Highway Geology Symposium, 2013. [11] MARCHELLI M, DE BIAGI V, PEILA D. Reliability-based design of rockfall passive systems height[J]. International Journal of Rock Mechanics and Mining Sciences, 2021, 139: 104664.1-104664.9. [12] MARCHELLI M, DE BIAGI V, PEILA D. Reliability-based design of protection net fences: influence of rockfall uncertainties through a statistical analysis[J]. Geosciences, 2020, 10(8): 280.1-280.24. [13] DE BIAGI V, MARCHELLI M, PEILA D. Reliability analysis and partial safety factors approach for rockfall protection structures[J]. Engineering Structures, 2020, 213: 110553.1-110553.1. [14] Austrian Standards Institute. Technical protection against rockfall-terms and definitions, effects of actions, design, monitoring and maintenance: ONR 24810[S]. Vienna: [s.n.], 2017. [15] MELZNER S, ROSSI M, GUZZETTI F. Impact of mapping strategies on rockfall frequency-size distributions[J]. Engineering Geology, 2020, 272: 105639.1-105639.11. [16] LAIMER H J. Determination of rockfall design blocks in Upper Triassic limestones and Dolomites (Dachstein Formation, Northern Calcareous Alps)[J]. Bulletin of Engineering Geology and the Environment, 2020, 79(3): 1581-1590. doi: 10.1007/s10064-019-01640-w [17] DE BIAGI V, LIA NAPOLI M, BARBERO M, et al. Estimation of the return period of rockfall blocks according to their size[J]. Natural Hazards and Earth System Sciences, 2017, 17(1): 103-113. doi: 10.5194/nhess-17-103-2017 [18] CURCEAC S, ATKINSON P M, MILNE A, et al. An evaluation of automated GPD threshold selection methods for hydrological extremes across different scales[J]. Journal of Hydrology, 2020, 585: 124845.1-124845.15. [19] STORMY A. MATLAB[M]. Fifth edition. Boston: [s.n.], 2019. [20] COLES S. An Introduction to Statistical Modeling of Extreme Values[M]. London: Springer, 2001. [21] 史道济. 实用极值统计方法[M]. 天津:天津科学技术出版社,2006. [22] 仇学艳,王超,秦崇仁. 阈值法在河海工程设计中的应用[J]. 水利学报,2001(8): 32-37. doi: 10.3321/j.issn:0559-9350.2001.08.006QIU Xueyan, WANG Chao, QIN Chongren. The application of threshold procedure in design of river and coastal engineering[J]. Shuili Xuebao, 2001(8): 32-37. doi: 10.3321/j.issn:0559-9350.2001.08.006 [23] PHILIP J, DAVID R, JENNY W, et al. Uncertainties in return from extreme value analysis of peaks over threshold using the generalised Pareto distribution[J]. Ocean Engineering, 2021, 220: 107725.1-107725.17. [24] THOMAS R K, ILARIA P. Use of peak over threshold data for flood frequency estimation: An application at the UK national scale.[J]. Jounal of Hydrology, 2023, 626(5): 130235.1-130235.31. [25] 中华人民共和国住房和城乡建设部. 建筑抗震设计规范GB 50011—2010(2016版)[S]. 北京:中国建筑工业出版社,2016. -