高速铁路区间完全中断场景双向到发线运用优化

高博文，欧冬秀，董德存; GAO Bowen; OU Dongxiu; DONG Decun

网刊加载中。。。

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

高速铁路区间完全中断场景双向到发线运用优化 PDF

- ORCID：
高博文 ^1,2
✉
- ORCID：
欧冬秀 ^1,2
✉
- ORCID：
董德存 ^1,2

1. 同济大学道路与交通工程教育部重点实验室，上海 201804； 2. 同济大学上海市轨道交通结构耐久与系统安全重点实验室，上海201804

中图分类号： U292.4

最近更新：2024-01-03

DOI：10.11908/j.issn.0253-374x.23386

目录contents

摘要

针对高速铁路区间完全中断场景下的列车运行调整问题，提出了双向到发线运用优化策略。从宏观角度将列车运行过程抽象为由事件和活动构成的网络，对提出的双向到发线运用策略进行建模，建立高速铁路区间完全中断场景下列车运行调整混合整数线性规划（MILP）模型。根据现场实时性要求，采用自编算法与商业优化软件Gurobi相结合的方式对模型进行求解。运用模型和算法对构建的9个现场中断案例进行算例测试。结果表明，双向到发线运用优化策略可以有效减少中断事件对列车运行的影响，相比于目前到发线运用策略，最多可减少10%的列车延误时间。

关键词

高速铁路; 区间中断; 到发线运用; 时刻表调整; 混合整数线性规划（MILP）

随着铁路智能化进程的不断推进，越来越多的国家如中国、瑞士、荷兰、日本等相继提出了一系列调度辅助决策系统^［

1］。各类突发事件下的列车智能调整方法是近年来铁路运输组织领域研究的重点问题之一。铁路系统作为一个开环系统，列车运行过程不可避免地受到各类干扰因素的影响，从而产生列车延误。在这种情况下，调度员需要根据历史经验对延误列车的作业计划进行及时调整，但很难应对大范围的延误干扰。国内外学者运用各类数学模型模拟列车运行作业过程，并尝试利用计算机算法求解列车运行调整方案。Cacchiani等^{［参考文献 2

百度学术}2］总结了列车运行实时调整的模型和算法。目前，对于小干扰事件的列车运行调整研究已经十分成熟，但缺乏对大干扰事件如中断事件的研究。

区间中断场景下，一旦中断事件发生，在中断事件恢复前列车将不被允许通过中断路段。Ghaemi等^［

3］引入浴缸模型来描述交通水平在中断事件下的变化。现有研究给出了解决列车运行调整问题的基本模型，按照模型类型大致可分为宏观和微观2类。宏观模型侧重于整个运输网络，应用“事件‒活动”网络重构列车时刻表，未考虑信号、列车速度和道岔切换等微观特征，而是将列车运行过程抽象为各个列车在各个车站到达和出发的事件集合以及连接任意2个事件的活动集合^{［参考文献 4-6}4-6］。宏观模型具有较低的计算复杂度，能够更有效地解决大范围列车运行调整问题。与此相反，微观模型关注运输网络的局部路段，需要获取列车和线路的详细信息，将运输网络抽象为一系列离散的行车资源，并通过构建等效的路网替代图来描述列车的运行过程^{［参考文献 7-9}7-9］。微观模型着重描述列车的微观行为，更适用于解决局部路网中的列车运行调整问题。

无论采用宏观模型还是微观模型，选择合适的调度策略至关重要。科学合理的调度策略能显著减少中断事件对时空的影响，缩短列车恢复到正常运行的时间。城市轨道交通领域已有多种适用于中断场景的调度策略，包括灵活停站^［

10］、临时折返变更交路^{［参考文献 10-11}10-11］、单线双向运行^{［参考文献 12

百度学术}12］、渡线反向行车^{［参考文献 13

百度学术}13］等。然而，由于旅客席位的预先预定，因此目前铁路领域的研究主要关注以下2种策略：第1种，针对单线铁路和双线高速铁路，在完全中断场景下将列车依次扣停在中断事件上游车站的策略^{［参考文献 14-15}14-15］；第2种，针对区间部分中断场景，利用区间正常轨道列车交替通过中断区间的策略^{［参考文献 16

百度学术}16］。相比之下，铁路领域中可选择的调度策略十分有限，因此铁路领域中的调度策略设计和优化需要得到更多的关注。

针对高速铁路区间完全中断场景，提出了双向到发线运用优化策略。对于高速铁路小干扰^［

17-18］和中断场景^{［参考文献 14-16}14-16］下的到发线运用研究均视双方向列车在车站内的作业相对独立，即每个列车仅使用车站内的正线和同侧到发线进行作业。考虑到大部分车站内的到发线具备双向作业条件，在文献［6］的基础上，运用“事件‒活动”网络对区间完全中断场景下的列车运行调整问题和到发线运用策略进行建模，以列车作业间隔、车站到发线数量、调度调整策略等为约束，构建混合整数线性规划（MILP）模型。最后，通过实际案例展示双向到发线运用优化策略在列车运行调整方案中的优化效果。

1 问题描述

1.1　问题定义

我国高速铁路均为双线铁路，列车运行方向分为上行和下行方向。通常情况下，上、下行列车在区间和车站的运行相对独立，互不干扰。然而，若发生区间完全中断事件，则双向列车均无法通过中断区间，需在上游车站内等候修复。此时，调度员需根据列车的运行特征和位置，在中断事件期间为各列车分配停靠站点；待事件修复后，需合理安排列车发车次序，缩短作业时间，促使各列车尽快恢复正点运行。

图1a展示了一个区间完全中断场景下的列车运行调整案例，该线路共有6个车站（ $s_{1}$ 、 $s_{2}$ ，…， $s_{6}$ ），每个车站到发线布置方式相同且共有4条到发线供列车作业使用。图1b展示了6个车站内的到发线布置情况。在正常情况下，2个方向上的列车数量相近。为避免交叉作业的影响，任意运行方向的列车仅使用车站内的正线及外侧到发线，即如图1b实线所示的2条列车作业径路。图1a中实线表示实际运行时刻表，虚线表示计划时刻表。中断发生后，已经进入中断区间的列车可以继续行驶，如列车 $g_{1}$ ；其余列车需要根据作业特征，由调度员按照经验扣停在中断区间上游车站，如列车 $g_{7} 、 g_{9}$ 扣停在车站 $s_{2}$ ，列车 $g_{3} 、 g_{5}$ 扣停在车站 $s_{3}$ 。

图1 区间完全中断情况下列车运行调整示例

Fig.1 Example of timetable rescheduling in case of a complete blockage

目前针对我国高速铁路区间中断场景的列车运行调整，均是将车站内的上、下行列车分开作业，每个方向上的列车只能使用图1b中的实线径路。实际上，每个车站的进、出站咽喉区具备多类渡线，使得一个方向上的列车能够利用对向到发线进行车站作业（如图1b中对应的虚线径路）。由图1a的调整结果可见，若列车能够更多地停靠在靠近中断区间的车站，则可以有效减缓事件影响向上游车站的传播速度。同时，由于在中断持续阶段受影响的列车（如列车 $g_{2} 、 g_{4} 、 g_{6} 、 g_{8}$ ）无法通过中断区间，因此中断区间的下游车站（如车站 $s_{1} 、 s_{2} 、 s_{3}$ ）有大量到发线空闲。若这部分的到发线被合理利用起来，供对向受中断影响的列车（列车 $g_{3} 、 g_{5} 、 g_{7} 、 g_{9}$ ）停靠或者通过作业，将有效减少中断事件引发的列车延误。

在区间中断事件发生后，考虑在中断区间2个运行方向上的上游车站内，将对向到发线使用权限开放给受中断影响的列车使用。由于使用对向到发线有可能造成2个方向上的列车在到发线上的交叉运行，从而引发新的延误事件，因此需要对区间中断场景下双向到发线运用方式进行优化。提出了3个前提条件（参见3节），在确保列车总延误不增加的情况下，令受影响的列车优先使用本运行方向的到发线；若本方向到发线均被占满，则根据列车总延误情况，酌情利用对向到发线进行作业。

1.2　问题假设

区间中断场景下的列车运行调整问题相对复杂，为简化模型，明确应用场景，作出如下假设：

（1）所有车站内的到发线均为配有乘降站台的双向到发线，并支持列车在任意到发线进行作业。

（2）所有区间均为双线区间，上、下行列车在区间内的作业相互独立、互不干扰。

（3）中断事件前，列车按计划时刻表准点运行。发生中断后，已进入中断区间的列车可正常通过，未进入中断区间的列车需在上游车站停车等待，称为“中断列车”。

（4）由于当日车次席位已售完，调整过程仅涉及作业次序和时间的变更，不考虑取消列车运行计划。

2 列车运行调整模型构建

2.1　问题抽象与参数定义

为了模拟列车在高速线路上的运行过程，从宏观角度将高速铁路网络抽象为由事件和活动构成的网络 $N = (E, A)$ 。每个列车的运行过程被抽象为列车在各个车站 $s_{i}$ 的到达事件 $d \in E_{a r r, s_{i}}$ 和出发事件 $f \in E_{d e p, s_{i}}$ 的集合；活动 $a \in A$ 连接网络中的任意2个事件 $e_{1}$ 、 $e_{2}$ ，记录着事件间状态变化的时间间隔。按照2个事件发生的列车主体 $t_{e_{1}}$ 、 $t_{e_{2}}$ 的不同，将活动集合 $A$ 划分为单车运行活动集合 $A_{t r a i n}$ 和多车间隔活动集合 $A_{h e a d w a y}$ ，集合 $A_{t r a i n}$ 表示同一列车主体在任意车站 $A_{t r a i n, s_{i}}$ 或任意区间 $A_{t r a i n, c_{i}}$ 的运行活动集合，集合 $A_{h e a d w a y}$ 表示相邻的不同列车主体发生的事件间隔活动集合。集合 $S 、 C 、 G$ 分别表示选定线路上的车站集合、区间集合和列车集合， $s_{i} \in S, c_{i} \in C$ 表示下行列车依次经过的第 $i$ 个车站与第 $i$ 个区间。对于任意事件 $e$ ， $s_{e}$ 、 $g_{e}$ 和 $r_{e}$ 分别表示对应的作业车站、列车主体和运行方向。模型的输入参数和决策变量如表1所示。

表1 参数定义

Tab.1 Parameter definitions

符号	定义
$o_{e}$	事件 $e$ 的计划作业时刻
$G_{s_{i}, r}$	车站 $s_{i}$ 在运行方向 $r$ 上的到发线数量
$T_{a r r, a r r}, T_{d e p, d e p}$	同一运行方向上任意到达事件间或任意出发事件间的最小间隔时间
$T_{a r r, a r r}^{'}, T_{d e p, d e p}^{'}$	不同运行方向上任意到达事件间或任意出发事件间的最小间隔时间
$T_{d e p, a r r}, T_{d e p, a r r}^{'}$	发生在同一到发线的同一或不同运行方向上的到达事件与出发事件间的最小间隔时间
$T_{c_{i}}$ ， $T_{s_{i}}$	区间 $c_{i}$ 的最小纯运行时间与车站 $s_{i}$ 的最小停站时间
$Δ_{a c c e}, Δ_{d e c e}$	列车在车站启、停作业的额外附加时间
$H_{s t a r t}, H_{e n d}$	中断事件的开始时间与结束时间
$H_{o p e n}$	对向到发线开放使用权限的结束时间
$c_{k}$	中断事件发生的区间位置
$M$	足够大的正整数
$H_{G u r o b i}$	Gurobi求解器的最大求解时间
$x_{e}$	事件 $e$ 经调整后的发生时刻
$p_{g_{e}, s_{i}}$	车站作业类型0‒1变量（若列车 $g_{e}$ 需要在车站 $s_{i}$ 停站作业，则值为1；否则，值为0）
$λ_{e_{1}, e_{2}}$	事件次序0‒1变量（若事件 $e_{1}$ 在事件 $e_{2}$ 之前发生，则值为1；否则，值为0）
$φ_{e, h}$	事件时序0‒1变量（若事件 $e$ 在时刻 $h$ 之前发生，则值为1；否则，值为0）
$μ_{g_{e}, s_{i}}$	事件到发线0‒1变量（若列车 $g_{e}$ 在车站 $s_{i}$ 对向到发，则值为1；否则，值为0）
$γ_{g_{e}, s_{i}}$	对向到发线使用状态0‒1变量（若列车 $g_{e}$ 在中断期间停留在车站 $s_{i}$ 对向到发线，则值为1；否则，值为0）

2.2　基本模型

中断发生后，双向列车均无法通过中断区间，调度员需逐一扣停在线运行的列车。事件结束后，调度员需迅速安排扣停的列车发车，并采用压缩区间运行时间、车站作业时间以及变更越行关系等手段，使列车尽快恢复计划时刻。为减缓突发事件对列车运行的影响，优化目标主要为最小化列车到发事件与原定计划的时间偏移量。基于2.1节定义的符号，得到不考虑事件影响的列车运行调整基本模型，如下所示：

m i n \sum_{e \in E} (x_{e} - o_{e})

（1）

s.t.

0 \leq x_{e} - o_{e}, \forall e \in E

（2）

x_{e} = o_{e}, \forall e \in E : o_{e} < H_{s t a r t}

（3）

p_{g_{d}, s_{i}} = 1, \forall o_{d} \neq o_{f}, s_{i} = s_{d} = s_{f}, g_{d} = g_{f}

（4）

x_{f} - x_{d} \geq T_{s_{i}} p_{g_{d}, s_{i}}, \forall s_{i} \in S, (f, d) \in A_{t r a i n, s_{i}}

（5）

Δ_{a c c e} p_{g_{d}, s_{i}} + T_{c_{i}, p l a n} + Δ_{d e c e} p_{g_{d}, s_{i + 1}} \geq x_{d} - x_{f} \geq Δ_{a c c e} p_{g_{d}, s_{i}} + T_{c_{i}} + Δ_{d e c e} p_{g_{d}, s_{i + 1}}, \forall (d, f) \in A_{t r a i n, c_{i}}

（6）

x_{e_{2}} - x_{e_{1}} + M (1 - λ_{e_{1}, e_{2}}) \geq T_{e_{1}, e_{2}}, \forall e_{1}, e_{2} \in E_{a r r, s_{i}} / e_{1}, e_{2} \in E_{d e p, s_{i}}, s_{i} \in S, (e_{1}, e_{2}) \in A_{h e a d w a y}

（7）

λ_{e_{1}, e_{2}} + λ_{e_{2}, e_{1}} = 1, \forall e_{1}, e_{2} \in E_{a r r, s_{i}} / e_{1}, e_{2} \in E_{d e p, s_{i}}, s_{i} \in S, (e_{1}, e_{2}) \in A_{h e a d w a y}

（8）

x_{f} - x_{d} + M λ_{f, d} \geq 1 - T_{d e p, a r r}, \forall d \in E_{a r r, s_{i}}, f \in E_{d e p, s_{i}}, s_{i} \in S, (f, d) \in A_{h e a d w a y}

（9）

x_{d} - x_{f} + M (1 - λ_{f, d}) \geq T_{d e p, a r r}, \forall d \in E_{a r r, s_{i}}, f \in E_{d e p, s_{i}}, s_{i} \in S, (f, d) \in A_{h e a d w a y}

（10）

\sum_{d \in E_{a r r, s_{i}} \ d^{'} : r_{d} = r_{d^{'}}} λ_{d, d^{'}} - \sum_{f \in E_{d e p, s_{i}} : g_{d} = g_{f}} λ_{f, d^{'}} \leq G_{s_{i}, r_{d^{'}}} - 1, \forall s_{i} \in S

（11）

λ_{f, f^{'}} = λ_{d, d^{'}}, \forall f, f^{'} \in E_{d e p, s_{i}}, d, d^{'} \in E_{a r r, s_{i, r_{d}}}, g_{d} = g_{f}, g_{d^{'}} = g_{f^{'}}, r_{d} = r_{d^{'}}

（12）

p_{g_{e}, s_{i}}, λ_{e_{1}, e_{2}} \in {0,1}, x_{e} \in N

（13）

目标函数式（1）表示最小化列车事件的调整时刻与计划时刻的偏差。式（2）表示各列车事件调整时间不得早于其计划时间；式（3）表示在突发事件发生前，各列车按照计划时刻表准点运行。式（4）和式（5）约束了计划时刻表中列车的作业类型。若计划中列车需要停站，则 $p_{g_{d}, s_{i}}$ 等于1；若计划中列车不需要停站，则 $p_{g_{d}, s_{i}}$ 根据调整情况取0或1。式（5）表示列车在车站作业需要满足的最小停站时间约束。式（6）表示列车在区间作业时需要满足的最小运行时间约束（为避免列车在区间长时间运行，将区间纯运行时间最大值设置为各列车计划时刻表中的纯运行时间 $T_{c_{i}, p l a n}$ ）。式（7）和式（8）共同约束了同一车站同一运行方向上的任意两列车的到达/出发事件间的最小间隔时间 $T_{a r r, a r r}$ / $T_{d e p, d e p}$ ；式（9）和式（10）共同约束了同一车站同一运行方向上同一到发线作业的任意两列车的出发事件和到达事件的最小间隔时间 $T_{d e p, a r r}$ 。式（11）表示当列车到达事件 $d^{'}$ 发生时，车站至少存在一条空闲到发线可供列车使用，即车站容量约束。式（12）表示区间列车越行约束，车站 $s_{i, r_{d}}$ 表示车站 $s_{i}$ 在运行方向 $r_{d}$ 的下一站。式（13）给出了各决策变量的取值范围。

2.3　调度策略建模

各方向上的列车在中断期间需依次扣停在车站，当且仅当中断事件恢复后才可发出。为实现这一控制策略，引入事件时序变量 $φ_{e, h}$ 来记录事件 $e$ 是否在某一时刻 $h$ 前发生，其与 $x_{e}$ 的关系如下所示：

h - x_{e} + M (1 - φ_{e, h}) \geq 1, e \in E

（14）

x_{e} - h + M φ_{e, h} \geq 0, e \in E

（15）

由事件时序变量 $φ_{e, h}$ 可得区间完全中断场景下的列车调度策略约束，如下所示：

x_{f} \geq H_{e n d}, \forall f \in E_{d e p, s_{c_{k}, r_{f}}} : H_{s t a r t} \leq o_{f} \leq H_{e n d}

（16）

\sum_{f \in E_{d e p, s_{i}} : o_{f} \leq H_{e n d}} φ_{f, H_{e n d}} - \sum_{d \in E_{a r r, s_{i, r}} : g_{d} = g_{f}} φ_{d, H_{e n d}} + \sum_{d^{'} \in E_{a r r, s_{i, r}} : o_{d^{'}} \leq H_{e n d}} φ_{d^{'}, H_{e n d}} - \sum_{f^{'} \in E_{d e p, s_{i, r}} : g_{d^{'}} = g_{f^{'}}} φ_{f^{'}, H_{e n d}} \leq G_{s_{i, r}, r}, \forall s_{i} \in S_{c_{k}, r} \ s_{c_{k}, r}

（17）

式（16）表示2个运行方向上的列车无法在中断期间通过中断区间（车站 $s_{c_{k}, r_{f}}$ 表示中断区间 $c_{k}$ 在运行方向 $r_{f}$ 的入口站）。式（17）中， $S_{c_{k}, r}$ 表示区间 $c_{k}$ 在运行方向 $r$ 上的上游车站集合， $s_{i, r}$ 表示车站 $s_{i}$ 在运行方向 $r$ 的下一个车站。 $(\sum φ_{f, H_{e n d}} - \sum φ_{d, H_{e n d}})$ 表示中断期间车站 $s_{i}$ 向下一站发出的且还未到站的列车数量， $(\sum φ_{d^{'}, H_{e n d}} - \sum φ_{f^{'}, H_{e n d}})$ 表示中断期间在下一站扣留的列车数量。由于中断期间需要根据车站内的到发线数量进行扣车作业，因此 $(\sum φ_{f, H_{e n d}} - \sum φ_{d, H_{e n d}})$ 和 $(\sum φ_{d^{'}, H_{e n d}} - \sum φ_{f^{'}, H_{e n d}})$ 要小于等于下一站的到发线数量。

3 双向到发线运用优化建模

所提出的双向到发线运用优化策略允许在中断事件发生后的一定时间内，各方向受到影响的中断列车可以借助对向到发线进行扣车或通过作业。为避免对向到发线的使用给对向列车带来过多的延误影响，列车只有在本方向到发线均被占用且对向车站到发线有空余时，才可使用对向到发线进行作业。因此，使用对向到发线的列车到达事件 $d$ 需要满足以下3个前提条件：

（1）事件 $d$ 的发生时刻在规定的可以使用对向到发线的时间范围内。

（2）在事件 $d$ 的发生时刻，车站内供本方向列车作业的到发线均被占用。

（3）在事件 $d$ 的发生时刻，车站内供对向列车作业的到发线至少有一个可供使用。

为方便表述，假设车站 $s_{i}$ 上运行方向为 $r$ 的列车可以使用对向到发线，将另一个运行方向表示为 $r^{'}$ 。

3.1　新的活动间隔约束

需要增加的活动间隔约束主要由2类构成。第1类，2个运行方向 $r$ 、 $r^{'}$ 上的列车有可能在对向到发线上交替作业，需要增加同一车站不同运行方向上列车到达、出发事件间的活动间隔约束；第2类，运行方向 $r$ 上的部分列车可能使用对向到发线，需要增加使用对向到发线的列车与其同运行方向使用本方向到发线的列车到达、出发事件间的活动间隔约束。新增的活动间隔约束如下所示：

M (1 - μ_{g_{\hat{e}}, s_{i}}) + x_{e_{2}} - x_{e_{1}} + M (1 - λ_{e_{1}, e_{2}}^{'}) \geq T_{e_{1}, e_{2}}^{'}

（18）

M (1 - μ_{g_{\hat{e}}, s_{i}}) + x_{e_{1}} - x_{e_{2}} + M λ_{e_{1}, e_{2}}^{'} \geq 1 - T_{e_{1}, e_{2}}^{'}

（19）

式中： $g_{\hat{e}}$ 为 $g_{e_{1}}$ 和 $g_{e_{2}}$ 中的可能在对向到发线作业的列车主体； $μ_{g_{\hat{e}}, s_{i}}$ 为事件到发线0‒1变量。

对于第1类活动间隔约束（ $r_{e_{1}} \neq r_{e_{2}}$ ），若2个事件 $e_{1}$ 、 $e_{2}$ 均为到达事件或者出发事件，则式（18）和式（19）共同约束了双向列车的到达间隔（ $T_{e_{1}, e_{2}}^{'} = T_{a r r, a r r}^{'}$ ）和出发间隔（ $T_{e_{1}, e_{2}}^{'} = T_{d e p, d e p}^{'}$ ）。当列车 $g_{\hat{e}}$ 不在车站对向到发线作业时，为使 $λ_{e_{1}, e_{2}}^{'}$ 或 $λ_{e_{2}, e_{1}}^{'}$ 都等于零，需要增加约束，如下所示：

λ_{e_{1}, e_{2}}^{'} + λ_{e_{2}, e_{1}}^{'} = μ_{g_{\hat{e}}, s_{i}}

（20）

若2个事件分别为到达事件和出发事件，则式（18）和式（19）共同约束了双向列车的到、发间隔或发、到间隔。当列车 $g_{\hat{e}}$ 不在对向到发线作业时，为使 $λ_{e_{1}, e_{2}}^{'}$ 或 $λ_{e_{2}, e_{1}}^{'}$ 都等于零，需要增加约束，如下所示：

λ_{e_{1}, e_{2}}^{'} \leq μ_{g_{\hat{e}}, s_{i}}

（21）

λ_{e_{2}, e_{1}}^{'} \leq μ_{g_{\hat{e}}, s_{i}}

（22）

对于第2类活动间隔约束（ $r_{e_{1}} = r_{e_{2}}$ ），增加约束的方式与第1类活动间隔约束相同，主要的区别是将运行方向为 $r^{'}$ 的列车主体替换成运行方向为 $r$ 的列车主体， $g_{\hat{e}}$ 列车主体保持不变。

3.2　前提条件约束

为满足设定的3个前提条件，需要对列车停站作业进行限制。首先，对于任意一个可以使用对向到发线的列车到达事件 $d$ ，需要满足前提条件（1）。式（23）约束了当事件 $d$ 发生在 $H_{o p e n}$ 后，即对向到发线不再开放给列车使用时， $μ_{g_{d}, s_{i}}$ 为零。

μ_{g_{d}, s_{i}} - φ_{d, H_{o p e n}} \leq 0, \forall d \in E_{a r r, s_{i}} : o_{d} \leq H_{o p e n}

（23）

其次，为满足前提条件（2）和前提条件（3），需要对事件 $d$ 发生时刻 $x_{d}$ 时车站内2个运行方向 $r$ 、 $r^{'}$ 上的车站到发线使用情况进行统计。2.2节基础模型中的车站容量约束式（11）被替换为：

M (1 - μ_{g_{d}, s_{i}}) \geq G_{s_{i}, r} - (m - p) \geq 1 - μ_{g_{d}, s_{i}}

（24）

M (1 - μ_{g_{d}, s_{i}}) + G_{s_{i}, r^{'}} μ_{g_{d}, s_{i}} - (p + n) \geq 1

（25）

G_{s_{i}, r^{'}} - m^{'} - p^{'} \geq 1

（26）

m = \sum_{d_{1} \in E_{a r r, s_{i}} \ d : r_{d_{1}} = r} λ_{d_{1}, d} - \sum_{f_{1} \in E_{d e p, s_{i}} : g_{f_{1}} = g_{d_{1}}} λ_{f_{1}, d}

n = \sum_{d_{2} \in E_{a r r, s_{i}} : r_{d_{2}} = r^{'}} λ_{d_{2}, d}^{'} - \sum_{f_{2} \in E_{d e p, s_{i}} : g_{f_{2}} = g_{d_{2}}} λ_{f_{2}, d}^{'}

p = \sum_{d_{3} \in E_{a r r, s_{i}} \ d : r_{d_{1}} = r, o_{d_{3}} \leq H_{o p e n}} λ_{d_{3}, d}^{'} - \sum_{f_{3} \in E_{d e p, s_{i}} : g_{f_{3}} = g_{d_{3}}} λ_{f_{3}, d}^{'}

式中： $m$ 为时刻 $x_{d}$ 停在车站 $s_{i}$ 内运行方向为 $r$ 的列车数量； $n$ 为时刻 $x_{d}$ 停在车站 $s_{i}$ 内运行方向为 $r^{'}$ 的列车数量； $p$ 为时刻 $x_{d}$ 停在车站 $s_{i}$ 对向到发线上的运行方向为 $r$ 的列车数量。式（24）、（25）约束了任意运行方向 $r$ 上的到达事件 $d$ 的车站容量。当 $μ_{g_{d}, s_{i}}$ 为0时，即列车停在车站本方向到发线， $(G_{s_{i}, r} - (m - p))$ 大于等于1，确保有到发线空闲；当 $μ_{g_{d}, s_{i}}$ 为1时，即列车停在车站的对向到发线， $(G_{s_{i}, r} - (m - p))$ 等于0， $(G_{s_{i}, r^{'}} - (p + n))$ 大于等于1，确保前提条件（2）、（3）成立。式（26）共同约束了任意运行方向 $r^{'}$ 上的到达事件 $d^{'}$ 的车站作业。 $m^{'}$ 、 $p^{'}$ 的计算方法如下所示：

m^{'} = \sum_{d_{1} \in E_{a r r, s_{i}} \ d^{'} : r_{d_{1}} = r^{'}} λ_{d_{1}, d^{'}} - \sum_{f_{1}^{'} \in E_{d e p, s_{i}} : g_{f_{1}} = g_{d_{1}}} λ_{f_{1}, d^{'}}

p^{'} = \sum_{d_{2} \in E_{a r r, s_{i}} : r_{d_{2}} = r, o_{d_{2}} \leq H_{o p e n}} λ_{d_{2}, d^{'}}^{'} - \sum_{f_{2} \in E_{d e p, s_{i}} : g_{f_{2}} = g_{d_{2}}} λ_{f_{2}, d^{'}}^{'}

式中： $m^{'}$ 表示时刻 $x_{d^{'}}$ 停在车站 $s_{i}$ 内运行方向为 $r^{'}$ 的列车数量； $p^{'}$ 表示时刻 $x_{d^{'}}$ 停在车站 $s_{i}$ 内对向到发线上的运行方向为 $r$ 的列车数量。为保证车站容量约束，当事件 $d^{'}$ 发生时，车站 $s_{i}$ 内供给运行方向 $r^{'}$ 的到发线应至少有一个空闲，即 $(G_{s_{i}, r^{'}} - m^{'} - p^{'})$ 大于等于1。

3.3　调度策略约束替换

由于开放了对向到发线作业的权限，部分列车 $g_{d}$ 可以停在下一站 $s_{i, r}$ 的对向到发线，因此本站 $s_{i}$ 可以发送更多的列车。2.3节中的式（17）需要进行调整。为了统计中断列车在中断期间是否停靠在车站 $s_{i, r}$ 的对向到发线，引入0‒1变量 $γ_{g_{d}, s_{i}}$ 。 $γ_{g_{d}, s_{i}}$ 与 $μ_{g_{d}, s_{i}}$ 、 $φ_{d, H_{e n d}}$ 、 $φ_{f, H_{e n d}}$ 的约束关系如下所示：

γ_{g_{d}, s_{i}} \leq μ_{g_{d}, s_{i}}

（27）

μ_{g_{d}, s_{i}} - 1 \leq γ_{g_{d}, s_{i}} - (φ_{d, H_{e n d}} - φ_{f, H_{e n d}}) \leq 1 - μ_{g_{d}, s_{i}}, g_{d} =

g_{f}, s_{d} = s_{f} = s_{i}

（28）

式（27）、（28）共同约束了当且仅当 $μ_{g_{d}, s_{i}}$ 为1、 $(φ_{d, H_{e n d}} - φ_{f, H_{e n d}})$ 为1（表示列车 $g_{d}$ 在中断期间未从车站 $s_{i}$ 发出）时， $γ_{g_{d}, s_{i}}$ 取1。根据变量 $γ_{g_{d}, s_{i}}$ ，式（17）需修改为

\sum_{f \in E_{d e p, s_{i}} : o_{f} \leq H_{e n d}} φ_{f, H_{e n d}} - \sum_{d \in E_{a r r, s_{i, r}} : g_{d} = g_{f}} φ_{d, H_{e n d}} + \sum_{d^{'} \in E_{a r r, s_{i, r}} : o_{d^{'}} \leq H_{e n d}} φ_{d^{'}, H_{e n d}} - \sum_{f^{'} \in E_{d e p, s_{i, r}} : g_{d^{'}} = g_{f^{'}}} φ_{f^{'}, H_{e n d}} - \sum_{d^{'} \in E_{a r r, s_{i, r}} : o_{d^{'}} \leq m i n (H_{o p e n}, H_{e n d})} γ_{g_{d^{'}}, s_{i, r}} \leq G_{s_{i, r}, r}

（29）

式（29）与式（17）的区别在于，若车站 $s_{i}$ 发出的列车即将停靠在车站 $s_{i, r}$ 内的对向到发线，或车站 $s_{i, r}$ 中的中断列车已经停靠在对向到发线上，则不计入车站 $s_{i}$ 的发车数量。

4 模型求解方法与有效不等式

构建的高速铁路列车运行调整模型属于混合整数线性规划模型的范畴并属于NP‒hard问题。目前对于此类问题的求解，商业求解软件如Gurobi、Cplex等可在短时间内提供有效的解决方案。考虑到优化软件Gurobi可以利用多核处理器进行并行计算，采用Gurobi求解器对构建的模型进行求解。

利用求解器求解需要消耗一些时间（假设为 $H_{G u r o b i}$ ），这使得在时段 $[H_{s t a r t}, H_{s t a r t} + H_{G u r o b i})$ 内无法提供可行列车调整方案。为贴合现场实时性的需求，利用自编算法对2个方向上的在线列车按照计划时刻表中的作业次序依次进行扣停作业（耗时在几毫秒内），输出的调整方案将供该时段使用，构建的列车运行调整模型仅对 $(H_{s t a r t} + H_{G u r o b i})$ 时刻后的列车作业进行调整。

为了进一步加快列车运行调整模型的求解速度，在求解过程中添加了一些有效不等式，这些不等式均是由同一列车的到达、出发事件的有序性而产生的。添加的有效不等式如下所示：

λ_{d, d^{'}} \geq λ_{f, d^{'}}, \forall d, d^{'} \in E_{a r r, s_{i}} : r_{d} = r_{d^{'}}, f \in E_{d e p, s_{i}} : g_{d} = g_{f}, s_{i} \in S

（30）

φ_{d, h} \geq φ_{f, h}, \forall d \in E_{a r r, s_{i}}, f \in E_{d e p, s_{i}} : g_{d} = g_{f}, s_{i} \in S

（31）

φ_{f, h} \geq φ_{d^{'}, h}, \forall d^{'} \in E_{a r r, s_{i, r}}, f \in E_{d e p, s_{i}} : g_{d} = g_{f}, s_{i} \in S

（32）

5 算例分析

5.1　场景描述

算例数据来源于2019年广珠城际高速铁路深湛线部分路段的运行图数据。整个线路从新会站至茂名站，沿线经过13个车站、12个区间。计划时刻表中包含29对上、下行列车数据，列车速度等级相同。各车站到发线数量配置如表2所示，其中 $(a, b)$ 表示该站正常情况下供上行列车作业的到发线数为 $a$ ，供下行列车作业的到发线数为 $b$ 。

表2 车站到发线数量

Tab.2 Number of station tracks

序号	站名	到发线数量	序号	站名	到发线数量	序号	站名	到发线数量
1	新会	（2，2）	6	恩平	（2，3）	11	马踏	（2，2）
2	江门深茂场	（1，1）	7	大槐	（1，1）	12	电白	（2，2）
3	双水镇	（2，2）	8	阳东	（2，2）	13	茂名	（3，3）
4	台山	（3，3）	9	阳江	（4，3）
5	开平南	（2，2）	10	阳西	（3，2）

各区间最小纯运行时间如表3所示，其中 $(a, b)$ 依次表示上、下行的最小纯运行时间，单位为min。

表3 区间最小纯运行时间

Tab.3 Pure running time of segments

序号	区间名	纯运行时间	序号	区间名	纯运行时间	序号	区间名	纯运行时间
1	新会―江门深茂场	（5，4）	5	开平南―恩平	（10，10）	9	阳江―阳西	（14，14）
2	江门深茂场―双水镇	（8，8）	6	恩平―大槐	（6，6）	10	阳西―马踏	（10，11）
3	双水镇―台山	（7，7）	7	大槐―阳东	（6，6）	11	马踏―电白	（7，7）
4	台山―开平南	（5，5）	8	阳东―阳江	（9，8）	12	电白―茂名	（11，11）

结合线路不同时段内列车开行密度的特征，通过设定中断事件的发生区间 $c_{k}$ 和开始时间 $H_{s t a r t}$ ，构建了3个中断场景；每个中断场景下，通过调整中断事件持续时间又衍生了3个中断案例。详细的场景信息如表4所示。其中，中断场景 $(a, b, c)$ 表示中断事件发生在 $a$ 区间、 $b$ 时刻，共持续 $c$ min。

表4 完全中断场景信息

Tab.4 Information of complete blockage scenarios

序号	中断场景	中断区间	序号	中断场景	中断区间	序号	中断场景	中断区间
1	（5，9，60）	开平南―恩平	4	（7，13，60）	大槐―阳东	7	（9，16，60）	阳江―阳西
2	（5，9，120）	开平南―恩平	5	（7，13， 20）	大槐―阳东	8	（9，16，120）	阳江―阳西
3	（5，9，180）	开平南―恩平	6	（7，13，180）	大槐―阳东	9	（9，16，180）	阳江―阳西

参照文献［

15-16］，同方向列车作业间隔

T_{a r r, a r r}

、

T_{d e p, d e p}

、

T_{d e p, a r r}

分别设置为3、2、3 min。停靠在对向到发线的列车与其同运行方向的其他列车的额外间隔时间与同向列车间隔时间设置相同。停靠在对向到发线的列车与其对向列车间的作业间隔

T_{a r r, a r r}^{'}

、

T_{d e p, d e p}^{'}

、

T_{d e p, a r r}^{'}

设置为0、0、3 min。列车最小停站时间

T_{s_{i}}

设置为2 min，列车启、停附加时间

Δ_{a c c e} 、 Δ_{d e c e}

均设置为1 min。对向到发线使用权限开放时间默认设为

[H_{s t a r t}, H_{e n d} + 60 m i n]

，参数

M

取1 440 min。

所有算例均在CPU为Inter（R） Core（TM） i7‒12700K、内存为32 GB、操作系统为Windows11的64位台式机上运用Python平台构建列车运行调整模型，并调用Gurobi 9.5.1求解（最大求解时间 $H_{G u r o b i}$ 取5 min）。

5.2　到发线运用优化前后结果对比

为了比较双向到发线运用优化前后的结果，针对每个中断案例构建了2类列车运行调整模型。第1类列车运行调整模型中不添加3节中的优化约束，每列车仅使用本运行方向到发线进行作业；第2类列车运行调整模型在中断事件发生后使用双向到发线运用优化策略。表5展示了9个中断案例的调整结果。

表5 不同中断案例下的列车运行调整结果

Tab.5 Results for different blockage cases

中断场景	仅使用本方向到发线作业			开放对向到发线作业				差值占比/%
中断场景	列车总延误时间/min	求解时间/s	Gap	列车总延误时间/min	求解时间/s	Gap	使用次数	差值占比/%
（5，9，60）	2 026	0.6	0	1 990	1.3	0	1	1.8
（5，9，120）	8 743	2.6	0	8 216	118.7	0	4	6.0
（5，9，180）	20 620	300.2	0.01	18 560	300.2	0.02	5	10.0
（7，13，60）	1 689	0.5	0	1 603	1.6	0	1	5.1
（7，13，120）	8 827	24.0	0	8 186	300.2	0.01	4	7.3
（7，13，180）	19 784	300.2	0.01	18 382	300.3	0.03	3	7.1
（9，16，60）	3 123	0.5	0	3 123	1.2	0	0	0.0
（9，16，120）	11 287	6.7	0	10 682	300.1	0	4	5.4
（9，16，180）	24 801	75.7	0	22 627	300.3	0.02	7	8.8

注：差值占比表示第一类模型结果与第二类模型结果的差值占第一类模型结果的百分比；Gap表示当前最优解与最优可能解之间的差距。

表5中每个中断案例下2类模型的求解Gap均小于等于0.03，说明所提出的模型均可在5 min内快速获得较优甚至最优的列车运行调整方案。对比2类模型的求解结果，除了中断场景（9，16，60）外，其余案例均通过使用对向到发线减少了列车总延误时间。随着中断事件影响的增大，列车对于对向到发线的使用次数也逐渐增加，并且列车每使用一次对向到发线对总延误的优化效果也随之变得显著。在这8个案例中，平均每使用一次对向到发线将有效减少219.4 min的列车延误。在案例（9，16，60）中，中断事件的影响相对较小，仅使用本方向到发线便可以很好地控制延误的传播，因此2类模型输出了同一种列车运行调整方案。为了更形象地显示2类模型输出的列车运行调整方案的差异，图2a和图2b（圆点表示对向到发线作业上的到达事件和出发事件）展示了案例（9，16，180）的2类模型输出结果。

图2 中断场景（9，16，180）的列车运行调整方案

Fig.2 Final timetable for blockage scenario (9,16,180)

由图2a可知，中断事件的影响在下行和上行方向已经传播到了序号为6和13的车站。从图2b来看，由于开放了对向到发线，中断事件的影响在下行和上行方向仅传播到了序号为9和11的车站。综合求解时间和求解方案的优劣，9个案例的测试结果反映出所提出的双向到发线运用优化策略可以根据中断事件的影响范围，在必要时充分利用对向到发线为中断列车提供更多的作业选择，在不影响对向列车作业的前提下，极大程度地减少中断事件造成的列车延误影响。

5.3　对向到发线开放时长灵敏度分析

在前文案例中，默认将对向到发线开放的时长设定为60 min。为了研究对向到发线开放时长对最终结果的影响，保持其他模型参数不变，将对向到发线的开放时长分别修改为0、30、90 min，获得的最终列车延误结果如图3所示。图3中，以开放时长0 min下的调整结果为基线，展示了各个案例其他开放时长下的调整结果与该基线的差值。

图3 不同对向到发线开放时长下的目标函数值

Fig.3 Objective function values for different opening times

随着开放时长由 $H_{e n d}$ 逐渐增加至 $H_{e n d}$ +60 min，案例2、3、9的列车延误数值有小幅下降，但延误改善效果并不明显。当开放时长继续增加至 $H_{e n d}$ +90 min时，案例3、6、9的列车总延误数值有所上升。3个案例均是3个中断场景中影响范围最大的，延长对向到发线的开放时间造成解空间的进一步扩大，加上将最大求解时间设置为5 min，有限的算力使得求解器求得了一个较差的解。综合测试结果可知，在小范围中断事件中可适度延长开放时间，而在大范围中断事件中宜缩短开放时间，以确保算力充足。

6 结语

提出了双向到发线运用优化策略，允许中断列车在部分车站使用对向到发线，以最小化列车总延误时间为目标构建混合整数线性规划模型。以广珠城际高速铁路为例，构建了9个中断案例。结果表明，基于双向到发线运用优化策略的模型输出的调整方案受中断影响范围更小，列车总延误时间更少。中断列车每使用一次对向到发线，平均减少219.4 min的列车延误。

作者贡献声明

高博文：模型设计与程序实现，最终版本修订。

欧冬秀：提出双向到发线运用优化方案。

董德存：区间中断场景下列车运行调整方案可行性验证。

参考文献

文超，冯永泰，胡瑞，等. 高速铁路智能调度辅助决策系统功能分析［J］. 中国铁路， 2020（7）：9. [百度学术]

WEN Chao， FENG Yongtai， HU Rui， et al. Functional analysis on intelligent traffic control auxiliary decision-making system of high speed railway ［J］. Chinese Railways， 2020（7）：9. [百度学术]

CACCHIANI V， HUISMAN D， KIDD M， et al. An overview of recovery models and algorithms for real-time railway rescheduling ［J］. Transportation Research， Part B： Methodological， 2014， 63： 15. [百度学术]

GHAEMI N， CATS O， GOVERDE R M P. Railway disruption management challenges and possible solution directions ［J］. Public Transport， 2017， 9： 343. [百度学术]

CADARSO L， MARIN A. Integration of timetable planning and rolling stock in rapid transit networks ［J］. Annals of Operations Research， 2012， 199： 113. [百度学术]

LOUWERSE I， HUISMAN D. Adjusting a railway timetable in case of partial or complete blockades ［J］. European Journal of Operational Research， 2014， 235（3）： 583. [百度学术]

GAO B， DONG D， WU Y， et al. Macroscopic real-time timetable rescheduling approach for high-speed railway under complete blockages using a three-stage algorithm ［J］. Transportation Research Record， 2022， 2676（3）： 535. [百度学术]

CORMAN F， D’ARIANO A， PACCIARELLI D， et al. Centralized versus distributed systems to reschedule trains in two dispatching areas ［J］. Public Transport， 2010， 2（3）： 219. [百度学术]

PELLEGRINI P， MARLIERE G， RODRIGUEZ J. Optimal train routing and scheduling for managing traffic perturbations in complex junctions ［J］. Transportation Research， Part B： Methodological， 2014， 59： 58. [百度学术]

GHAEMI N， CATS O， GOVERDE R M P. A microscopic model for optimal train short-turnings during complete blockages ［J］. Transportation Research， Part B： Methodological， 2017， 105： 423. [百度学术]

ZHU Y， GOVERDE R M P. Railway timetable rescheduling with flexible stopping and flexible short-turning during disruptions ［J］. Transportation Research， Part B： Methodological， 2019， 123： 149. [百度学术]

周玮腾，钱蕾，韩宝明，等. 运营区间双向中断下城市轨道交通列车运行调整［J］. 东南大学学报（自然科学版）， 2022， 52（4）： 770. [百度学术]

ZHOU Weiteng， QIAN Lei， HAN Baoming， et al. Train rescheduling of urban rail transit under bi-directional interruption in operation section ［J］. Journal of Southeast University（Natural Science Edition）， 2022， 52（4）： 770. [百度学术]

刘浩江. 单线双向运行在地铁行车组织中的应用［J］. 都市快轨交通， 2008， 21（3）： 17. [百度学术]

LIU Haojiang. Application of bidirectional operation on single track in subway train operation organization ［J］. Urban Rapid Rail Transit， 2008， 21（3）： 17. [百度学术]

XU X， LI K， YANG L. Rescheduling subway trains by a discrete event model considering service balance performance ［J］. Applied Mathematical Modelling， 2016， 40（2）： 1446. [百度学术]

孟令云，杨肇夏，李海鹰. 单线铁路区间能力失效条件下列车运行调整模型［J］. 系统工程理论与实践， 2012， 32（4）： 885. [百度学术]

MENG Lingyun， YANG Zhaoxia， LI Haiying. Train dispatching models under field capacity breakdowns on single-track railway lines ［J］. Systems Engineering： Theory & Practice， 2012， 32（4）： 885. [百度学术]

占曙光，赵军，彭其渊，等. 高速铁路区间能力全失效条件下列车运行实时调整研究［J］. 铁道学报， 2015， 37（11）： 1. [百度学术]

ZHAN Shuguang， ZHAO Jun， PENG Qiyuan， et al. Real-time train rescheduling on high-speed railway under complete segment blockages ［J］. Journal of the China Railway Society， 2015， 37（11）： 1. [百度学术]

占曙光，赵军，彭其渊. 高速铁路区间能力部分失效情况下列车运行实时调整研究［J］. 铁道学报， 2016， 38（10）： 1. [百度学术]

ZHAN Shuguang， ZHAO Jun， PENG Qiyuan. Real-time train rescheduling on high-speed railway under partial segment blockages ［J］. Journal of the China Railway Society， 2016， 38（10）： 1. [百度学术]

彭其渊，宁佳，鲁工圆. 大型高铁客运站到发线运用调整模型及算法［J］. 铁道学报， 2019， 41（1）： 10. [百度学术]

PENG Qiyuan， NING Jia， LU Gongyuan. Model and algorithm for train platform scheme rescheduling at large high-speed railway station ［J］. Journal of the China Railway Society， 2019， 41（1）： 10. [百度学术]

王宇晴，查伟雄，万平. 列车运行调整与到发线运用的协同优化研究［J］. 科学技术与工程， 2022， 22： 9797. [百度学术]

WANG Yuqing， ZHA Weixiong， WAN Ping. Coordinated optimization of train operation adjustment and arrival departure-track operation ［J］. Science Technology and Engineering， 2022， 22： 9797. [百度学术]

高速铁路区间完全中断场景双向到发线运用优化 PDF

摘要

关键词

1 问题描述

1.1 问题定义

1.2 问题假设

2 列车运行调整模型构建

2.1 问题抽象与参数定义

2.2 基本模型

2.3 调度策略建模