2.793

2018影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于闭塞区间的高速列车运行时间与节能协同优化方法

赵辉 代学武

赵辉, 代学武. 基于闭塞区间的高速列车运行时间与节能协同优化方法. 自动化学报, 2020, 46(3): 471−481. doi: 10.16383/j.aas.c190200
引用本文: 赵辉, 代学武. 基于闭塞区间的高速列车运行时间与节能协同优化方法. 自动化学报, 2020, 46(3): 471−481. doi: 10.16383/j.aas.c190200
Zhao Hui, Dai Xue-Wu. Cooperative optimization method for high-speed trains running time and energy saving based on block sections. Acta Automatica Sinica, 2020, 46(3): 471−481. doi: 10.16383/j.aas.c190200
Citation: Zhao Hui, Dai Xue-Wu. Cooperative optimization method for high-speed trains running time and energy saving based on block sections. Acta Automatica Sinica, 2020, 46(3): 471−481. doi: 10.16383/j.aas.c190200

基于闭塞区间的高速列车运行时间与节能协同优化方法


DOI: 10.16383/j.aas.c190200
详细信息
    作者简介:

    东北大学博士研究生. 分别于2014年和2016年获得大连海事大学学士和硕士学位. 主要研究方向为智能交通系统, 高速列车调度与控制. E-mail: zhaohui_209@163.com

    东北大学流程工业综合自动化国家重点实验室教授. 主要研究方向为动态系统鲁棒状态估计, 无线传感测量与控制, 状态监测, 及其在工业物联网、高铁调度等领域的应用. 本文通信作者. E-mail: daixuewu@mail.neu.edu.cn

  • 基金项目:  国家自然科学基金(61790574, 61773111, U1834211), 中国国家铁路集团有限公司科技研究开发计划课题(N2019G020)资助

Cooperative Optimization Method for High-speed Trains Running Time and Energy Saving Based on Block Sections

More Information
  • Fund Project:  Supported by National Natural Science Foundation of China (61790574, 61773111, U1834211) and Science and Technology Project of China National Railway Group Co., Ltd. (N2019G020)
  • 摘要: 提出了一种高速列车运行时间与节能协同优化方法. 针对由动态调度层、优化控制层、跟踪控制层组成的列车运行控制与动态调度一体化结构, 设计了面向动态调度层和优化控制层的列车运行时间调整策略和节能速度位置曲线. 基于高速铁路闭塞区间, 建立了列车 − 区间模型和列车速度曲线节能优化模型. 利用模型预测控制方法对列车区间运行时间进行调整, 优化列车总延误时间; 根据调整后的区间运行时间设计列车运行优化速度位置曲线, 减少列车运行能耗. 仿真算例验证了设计的运行时间与节能协同优化策略的有效性.
  • 图  1  动态调度和运行控制一体化结构图

    Fig.  1  Integrated structure diagram of dynamic rescheduling and operation control

    图  2  列车11在区间分界点延误时间

    Fig.  2  Delay time of train 11 at section demarcation point

    图  4  列车16在区间分界点延误时间

    Fig.  4  Delay time of train 16 at section demarcation point

    图  3  列车13在区间分界点延误时间

    Fig.  3  Delay time of train 13 at section demarcation point

    图  5  列车13优化速度位置曲线

    Fig.  5  The optimal speed-position trajectory of train 13

    图  6  速度位置曲线跟踪效果图

    Fig.  6  Tracking performance of speed-position trajectory

    图  7  跟踪效果局部放大图

    Fig.  7  Partial enlarged detail of tracking performance

  • [1] 孙亚华. 高速列车运行控制与动态调度一体化的研究[硕士学位论文], 北京交通大学, 2015.

    Sun Ya-Hua. The Integration Study of High-Speed Train Operation Control and Scheduling [Master thesis], Beijing Jiaotong University, 2015.
    [2] 2 Gao S G, Hou Y H, Dong H R, Stichel S, Ning B. High-speed trains automatic operation with protection constraints: a resilient nonlinear gain-based feedback control approach. IEEE/CAA Journal of Automatica Sinica, 2019, 6(4): 992−999 doi:  10.1109/JAS.2019.1911582
    [3] 3 Gao S G, Dong H R, Ning B, Zhang Q. Cooperative prescribed performance tracking control for multiple high-speed trains in moving block signaling system. IEEE Transactions on Intelligent Transportation Systems, 2019, 20(7): 2740−2749 doi:  10.1109/TITS.2018.2877171
    [4] 4 Albrecht A R, Howlett P G, Pudney P J, Vu X. Energy-efficient train control: from local convexity to global optimization and uniqueness. Automatica, 2013, 49(10): 3072−3078 doi:  10.1016/j.automatica.2013.07.008
    [5] 5 Yan X H, Cai B G, Ning B, Wei S G. Moving horizon optimization of dynamic trajectory planning for high-speed train operation. IEEE Transactions on Intelligent Transportation Systems, 2016, 17(5): 1258−1270 doi:  10.1109/TITS.2015.2499254
    [6] 6 Yang H, Zhang K P, Liu H E. Online regulation of high speed train trajectory control based on T-S fuzzy bilinear model. IEEE Transactions on Intelligent Transportation Systems, 2016, 17(6): 1496−1508 doi:  10.1109/TITS.2015.2497320
    [7] 7 Li L, Dong W, Ji Y D, Zhang Z K, Lang T. Minimal-energy driving strategy for high-speed electric train with hybrid system model. IEEE Transactions on Intelligent Transportation Systems, 2013, 14(4): 1642−1653 doi:  10.1109/TITS.2013.2265395
    [8] 8 Fang W, Yang S X, Yao X. A survey on problem models and solution approaches to rescheduling in railway networks. IEEE Transactions on Intelligent Transportation Systems, 2016, 16(6): 2997−3016
    [9] 9 Binder S, Maknoon Y, Bierlaire M. The multi-objective railway timetable rescheduling problem. Transportation Research Part C Emerging Technologies, 2017, 78: 78−94 doi:  10.1016/j.trc.2017.02.001
    [10] 10 Dong H R, Zhu H N, Li Y D, Lv Y S, Gao S G, Zhang Q, et al. Parallel intelligent systems for integrated high-speed railway operation control and dynamic scheduling. IEEE Transactions on Cybernetics, 2018, 48(12): 3381−3389 doi:  10.1109/TCYB.2018.2852772
    [11] 11 Zheng Y J. Emergency train scheduling on Chinese high-speed railways. Transportation Science, 2018, 52(5): 1077−1091 doi:  10.1287/trsc.2017.0794
    [12] 12 Krasemann J T. Design of an effective algorithm for fast response to the re-scheduling of railway traffic during disturbances. Transportation Research Part C Emerging Technologies, 2012, 20(1): 62−78 doi:  10.1016/j.trc.2010.12.004
    [13] 13 Shakibayifar M, Sheikholeslami A, Corman F. A simulation-based optimization approach to rescheduling train traffic in uncertain conditions during disruptions. Scientia Iranica, 2018, 25(2): 646−662
    [14] 14 Yang L X, Qi J G, Li S K, Gao Y. Collaborative optimization for train scheduling and train stop planning on high-speed railways. Omega, 2016, 64: 57−76 doi:  10.1016/j.omega.2015.11.003
    [15] Zhou W L, Tian J L, Qin J, Deng L B, Wei T J. Optimization of multiperiod mixed train schedule on high-speed railway. Discrete Dynamics in Nature and Society, 2015, Article ID: 107048
    [16] 16 Hou Z P, Dong H R, Gao S G, Nicholson G, Chen L, Roberts C. Energy-saving metro train timetable rescheduling model considering ATO profiles and dynamic passenger flow. IEEE Transactions on Intelligent Transportation Systems, 2019, 20(7): 2774−2785 doi:  10.1109/TITS.2019.2906483
    [17] 17 Yan X H, Cai B, Ning B, Wei S G. Online distributed cooperative model predictive control of energy-saving trajectory planning for multiple high-speed train movements. Transportation Research Part C Emerging Technologies, 2016, 69: 60−78 doi:  10.1016/j.trc.2016.05.019
    [18] 范家璐, 张也维, 柴天佑. 一类工业过程运行反馈优化控制方法. 自动化学报, 2015, 41(10): 1754−1761

    18 Fan Jia-Lu, Zhang Ye-Wei, Chai Tian-You. Optimal operational feedback control for a class of industrial processes. Acta Automatica Sinica, 2015, 41(10): 1754−1761
    [19] 19 Li S K, Zhou X S, Yang L X, Gao Z Y. Automatic train regulation of complex metro networks with transfer coordination constraints: a distributed optimal control framework. Transportation Research Part B Methodological, 2018, 117: 228−253 doi:  10.1016/j.trb.2018.09.001
    [20] 20 Li S K, Yang L X, Gao Z Y. Efficient real-time control design for automatic train regulation of metro loop lines. IEEE Transactions on Intelligent Transportation Systems, 2019, 20(2): 485−496 doi:  10.1109/TITS.2018.2815528
    [21] 21 Wang X, Li S K, Su S, Tang T. Robust fuzzy predictive control for automatic train regulation in high-frequency metro lines. IEEE Transactions on Fuzzy Systems, 2019, 27(6): 1295−1308 doi:  10.1109/TFUZZ.2018.2877593
    [22] 22 Van Breusegem S, Campion G, Bastin G. Traffic modeling and state feed back control for metro lines. IEEE Transactions on Automatic Control, 1991, 36(7): 770−784 doi:  10.1109/9.85057
    [23] 23 Su S, Tang T, Li X, Gao Z Y. Optimization of multitrain operations in a subway system. IEEE Transactions on Intelligent Transportation Systems, 2014, 15(2): 673−684 doi:  10.1109/TITS.2013.2285737
    [24] Zhao H, Dai X W, Zhang M, Zhang Q. An event-triggered robust cruise control strategy for multiple high-speed trains. In: Proceedings of the 30th Chinese Control and Decision Conference. Shenyang, China: IEEE, 2018. 2041−2046
    [25] Tang H Y, Ge X C, Liu Q Q, Wang Q Y. Robust H control of high-speed trains with parameter uncertainties and unpredictable time-varying delays. In: Proceedings of the 35th Chinese Control Conference. Chengdu, China: IEEE, 2016. 10173−10178
    [26] 26 Dong H R, Gao S G, Ning B, Li L. Extended fuzzy logic controller for high speed train. Neural Computing and Applications, 2013, 22(2): 321−328 doi:  10.1007/s00521-011-0681-8
    [27] 27 Zhou Y H, Zhang Z L. High-speed train control based on multiple-model adaptive control with second-level adaptation. Vehicle System Dynamics, 2014, 52(5): 637−652 doi:  10.1080/00423114.2014.887209
    [28] 28 Ji H H, Hou Z S, Zhang R K. Adaptive iterative learning control for high-speed trains with unknown speed delays and input saturations. IEEE Transactions on Automation Science and Engineering, 2016, 13(1): 260−273 doi:  10.1109/TASE.2014.2371816
    [29] 29 Tang H Y, Wang Q Y, Feng X Y. Robust stochastic control for high-speed trains with nonlinearity, parametric uncertainty, and multiple time-varying delays. IEEE Transactions on Intelligent Transportation Systems, 2018, 19(4): 1027−1037 doi:  10.1109/TITS.2017.2710138
    [30] 徐君, 张国良, 曾静, 杜柏阳, 贾枭. 高阶离散多智能体系统在参数不确定和带外部干扰下的鲁棒H一致性控制. 自动化学报, 2017, 43(10): 1850−1857

    30 Xu Jun, Zhang Guo-Liang, Zeng Jing, Du Bo-Yang, Jia Xiao. Robust H consensus control for high-order discrete time multi-agent systems with parameter uncertainties and external disturbances. Acta Automatica Sinica, 2017, 43(10): 1850−1857
    [31] 俞立. 鲁棒控制——线性矩阵不等式处理方法. 北京: 清华大学出版社, 2002. 41−59

    Yu Li. Robust Control — Linear Matrix Inequality Method. Beijing: Tsinghua University Press, 2002. 41−59
  • [1] 贾兆红, 王燕, 张以文. 求解差异机器平行批调度的双目标协同蚁群算法[J]. 自动化学报, 2020, 46(6): 1121-1135. doi: 10.16383/j.aas.c180834
    [2] 王秀丽, 姜斌, 陆宁云. 基于相关向量机的高速列车牵引系统剩余寿命预测[J]. 自动化学报, 2019, 45(12): 2303-2311. doi: 10.16383/j.aas.c190204
    [3] 王伟, 任国恒, 陈立勇, 张效尉. 基于结构先验与协同优化的城市场景分段平面重建[J]. 自动化学报, 2019, 45(11): 2187-2198. doi: 10.16383/j.aas.2017.c170458
    [4] 王呈, 陈晶, 荀径, 李开成. 基于混合滤波最大期望算法的高速列车建模[J]. 自动化学报, 2019, 45(12): 2260-2267. doi: 10.16383/j.aas.c190193
    [5] 原豪男, 郭戈. 交通信息物理系统中的车辆协同运行优化调度[J]. 自动化学报, 2019, 45(1): 143-152. doi: 10.16383/j.aas.c180354
    [6] 关新平, 陈彩莲, 杨博, 华长春, 吕玲, 朱善迎. 工业网络系统的感知-传输-控制一体化:挑战和进展[J]. 自动化学报, 2019, 45(1): 25-36. doi: 10.16383/j.aas.c180484
    [7] 杨超, 彭涛, 阳春华, 陈志文, 桂卫华. 高速列车牵引传动系统故障测试与验证仿真平台研究[J]. 自动化学报, 2019, 45(12): 2218-2232. doi: 10.16383/j.aas.c190395
    [8] 宁滨, 董海荣, 郑伟, 荀径, 高士根, 王洪伟, 孟令云, 李浥东. 高速铁路运行控制与动态调度一体化的现状与展望[J]. 自动化学报, 2019, 45(12): 2208-2217. doi: 10.16383/j.aas.2019.y000004
    [9] 司方远, 汪晋宽, 韩英华, 赵强. 信息物理融合的智慧能源系统多级对等协同优化[J]. 自动化学报, 2019, 45(1): 84-97. doi: 10.16383/j.aas.2018.c180368
    [10] 陶宏伟, 彭涛, 杨超, 陈志文, 桂卫华. 高速列车牵引整流器多类故障联合诊断方法[J]. 自动化学报, 2019, 45(12): 2294-2302. doi: 10.16383/j.aas.c190258
    [11] 高开来, 丁进良. 蒸馏与换热协同的约束多目标在线操作优化方法[J]. 自动化学报, 2019, 45(9): 1679-1690. doi: 10.16383/j.aas.c180717
    [12] 周东华, 纪洪泉, 何潇. 高速列车信息控制系统的故障诊断技术[J]. 自动化学报, 2018, 44(7): 1153-1164. doi: 10.16383/j.aas.2018.c170392
    [13] 丁进良, 杨翠娥, 陈远东, 柴天佑. 复杂工业过程智能优化决策系统的现状与展望[J]. 自动化学报, 2018, 44(11): 1931-1943. doi: 10.16383/j.aas.2018.c180550
    [14] 伍铁斌, 阳春华, 李勇刚, 朱红求, 桂卫华. 基于模糊操作模式的砷盐除钴过程操作参数协同优化[J]. 自动化学报, 2014, 40(8): 1690-1698. doi: 10.3724/SP.J.1004.2014.01690
    [15] 衷路生, 李兵, 龚锦红, 张永贤, 祝振敏. 高速列车非线性模型的极大似然辨识[J]. 自动化学报, 2014, 40(12): 2950-2958. doi: 10.3724/SP.J.1004.2014.02950
    [16] 左斌, 李静, 胡云安. 基于反演设计的块控型极值搜索系统一体化控制方法研究[J]. 自动化学报, 2011, 37(9): 1114-1129. doi: 10.3724/SP.J.1004.2011.01114
    [17] 郑志刚, 汪增福. 基于区域间协同优化的立体匹配算法[J]. 自动化学报, 2009, 35(5): 469-477. doi: 10.3724/SP.J.1004.2009.00469
    [18] 王飞跃. 关于控制与机械一体化的数理计算体系结构研究[J]. 自动化学报, 2006, 32(2): 318-320.
    [19] 韩志刚. 关于建模与自适应控制的一体化途径[J]. 自动化学报, 2004, 30(3): 380-389.
    [20] 孙福权, 郑秉霖, 崔建江, 汪定伟, 刘新胜. 炼钢-热轧一体化管理的生产计划编制问题研究[J]. 自动化学报, 2000, 26(3): 409-413.
  • 加载中
图(7)
计量
  • 文章访问数:  1604
  • HTML全文浏览量:  615
  • PDF下载量:  95
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-20
  • 录用日期:  2019-08-15
  • 网络出版日期:  2020-03-30

基于闭塞区间的高速列车运行时间与节能协同优化方法

doi: 10.16383/j.aas.c190200
    基金项目:  国家自然科学基金(61790574, 61773111, U1834211), 中国国家铁路集团有限公司科技研究开发计划课题(N2019G020)资助
    作者简介:

    东北大学博士研究生. 分别于2014年和2016年获得大连海事大学学士和硕士学位. 主要研究方向为智能交通系统, 高速列车调度与控制. E-mail: zhaohui_209@163.com

    东北大学流程工业综合自动化国家重点实验室教授. 主要研究方向为动态系统鲁棒状态估计, 无线传感测量与控制, 状态监测, 及其在工业物联网、高铁调度等领域的应用. 本文通信作者. E-mail: daixuewu@mail.neu.edu.cn

摘要: 提出了一种高速列车运行时间与节能协同优化方法. 针对由动态调度层、优化控制层、跟踪控制层组成的列车运行控制与动态调度一体化结构, 设计了面向动态调度层和优化控制层的列车运行时间调整策略和节能速度位置曲线. 基于高速铁路闭塞区间, 建立了列车 − 区间模型和列车速度曲线节能优化模型. 利用模型预测控制方法对列车区间运行时间进行调整, 优化列车总延误时间; 根据调整后的区间运行时间设计列车运行优化速度位置曲线, 减少列车运行能耗. 仿真算例验证了设计的运行时间与节能协同优化策略的有效性.

English Abstract

赵辉, 代学武. 基于闭塞区间的高速列车运行时间与节能协同优化方法. 自动化学报, 2020, 46(3): 471−481. doi: 10.16383/j.aas.c190200
引用本文: 赵辉, 代学武. 基于闭塞区间的高速列车运行时间与节能协同优化方法. 自动化学报, 2020, 46(3): 471−481. doi: 10.16383/j.aas.c190200
Zhao Hui, Dai Xue-Wu. Cooperative optimization method for high-speed trains running time and energy saving based on block sections. Acta Automatica Sinica, 2020, 46(3): 471−481. doi: 10.16383/j.aas.c190200
Citation: Zhao Hui, Dai Xue-Wu. Cooperative optimization method for high-speed trains running time and energy saving based on block sections. Acta Automatica Sinica, 2020, 46(3): 471−481. doi: 10.16383/j.aas.c190200
  • 近年来, 高速铁路飞速发展, 以其快速、安全、准时等特点极大地改变了人们的出行方式. 高速铁路逐渐在中国轨道交通运输体系中居于核心地位. 随着高速铁路网络的快速发展, 其已呈现出运营里程长、路网复杂、区域分布广、客流量大等特点, 对高速列车运行控制和动态调度策略提出了更高的要求. 这就要求我们不断地改进高速列车的运行控制和动态调度方法, 实现高速列车运行快速、运营安全、旅客舒适、到站准时等服务目标[1-3].

    高速列车运行控制的目的在于优化列车运行能耗, 控制列车按期望的速度位置曲线安全运行. 列车节能优化运行一直是列车运行控制领域的研究热点. Albrecht等[4]指出列车节能优化运行主要包括最大牵引、巡航、惰行、最大制动四个过程. Yan等[5]研究了具有不确定参数影响下的列车运行轨迹的动态优化问题, 将列车运行中能量消耗作为性能指标, 采用滚动时域优化策略优化了列车速度位置曲线. 文献[6]则基于列车运行的非线性动态特征和列车历史运行数据, 构建了一个模糊双线性模型来描述高速列车的动态运行过程, 设计了自适应预测控制方法来实现高速列车的安全、节能运行. 文献[7]提出了一种混合系统模型来描述高速列车的特性, 将列车运行能耗转化为非线性优化问题, 设计了一种固定运行时间下的最小能耗运行策略.

    列车动态调度主要是指: 由于灾害天气、线路和设备故障等突发事件导致列车实绩运行图偏离计划运行图时, 需要通过一次或多次制订列车运行阶段调整计划, 以尽快恢复计划运行图的动态调整过程. 动态性主要体现在需要根据不断变化的高铁路网运行状态不断调整列车运行, 因此要求动态调度策略具有实时性、稳定性和可靠性的特点. 已有研究表明, 列车动态调度问题具有多目标、多约束、非线性等特征[8]. 列车动态调度目标通常设置为: 总晚点列车数最少、列车总晚点时间最小、加权总晚点时间最小、列车正点率最高、列车到站时间偏差和列车离站时间偏差加权和最小、旅客满意度高等一个或多个目标[9]. 针对此类多目标优化问题, 列车动态调度手段主要包括改变列车到发时间、压缩列车区间运行时分、缩减列车车站停车时间、改变列车运行线路、改变列车到发顺序等[10]. 综上所述, 调整列车运行时间是列车动态调度的重要组成部分. 当列车延误不是特别严重时, 利用基本图中预留的缓冲时间余量, 对列车运行时间进行调整就能够恢复列车正常运行. 列车运行时间调整问题的优化求解方法主要有三类: 数学规划方法[11]、人工智能算法[12]和仿真方法[13]. 例如, 文献[14]将车站停车时间作为约束考虑到列车运行时间调整中, 将停车时间和延误时间作为优化目标, 考虑了高速列车运行时间调整和车站停车时间的协同优化问题, 给出了列车最优运行方式. Zhou等[15]采用启发式算法研究了多周期混合条件下高速列车运行时间优化问题, 在运行时间调整问题中考虑了旅客满意度, 得到了旅客满意度较高的优化解.

    当列车运行偏离基本运行图时, 恢复按图运行的过程实质上与列车运行控制和动态调度同时相关, 动态调度提供策略, 运行控制去执行. 运行控制与动态调度是存在内在联系的, 当前的研究成果多是只考虑运行控制或是只考虑动态调度, 综合考虑列车运行控制和动态调度的学者还不多[16]. 随着技术的进步, 相关理论的完善以及高铁基础设施的建设, 高速铁路的智能化已经成为未来高铁路网的发展趋势. 而在高速列车运行控制和动态调度研究领域, 高速列车的运行控制与动态调度一体化理论和技术正在成为研究重点. 如文献[17]采用分布式协同模型预测控制和蚁群算法, 基于高速列车运行控制与动态调度一体化的思想, 通过微调列车运行时间, 在相邻列车运行距离约束、相邻列车运行间隔约束、列车运行速度约束等约束条件下, 重新设计了运行在线路上列车的速度位置曲线, 在约束允许的范围内优化了列车运行能耗.

    目前, 高速列车的运行控制与动态调度还相对独立, 针对列车运行控制与动态调度一体化的研究, 还处于初步研究阶段, 研究成果也相对较少. 高速列车运行控制和动态调度一体化主要问题在于如何将业务流程中相对独立、具有不同优化目标、但又密切相关的运行控制系统和调度指挥系统更好地结合起来, 通过信息共享和信息融合, 提高运行控制系统和动态调度系统之间的协调能力. 由于动态调度与运行控制的优化目标不尽相同且各有侧重, 有时还相互冲突, 如何利用自动化和信息技术对动态调度和运行控制进行多目标协同优化也是一体化研究中的难点.

    针对上述问题, 本文以调度系统中运行时间调整问题和控制系统中节能优化运行问题为研究对象, 借鉴工业自动化系统中复杂系统的分层控制体系和全流程一体化运行反馈控制思想[18], 提出了高速列车运行时间和节能协同优化方法. 以列车繁忙时段小范围延误为研究背景, 借鉴城市轨道交通离散事件建模理论[19-22], 建立了基于列车运行时间的列车 − 区间模型. 在此基础上, 采用模型预测控制设计列车运行时间调整策略, 减少列车总延误时间和区间运行时间压缩量, 恢复列车按图运行. 根据调整后的列车运行时间, 采用数值分析和二分法迭代相结合的方法, 在给定时间和速度约束条件下, 优化列车运行能耗, 生成优化速度位置曲线. 利用跟踪控制跟踪优化速度位置曲线. 同时, 将列车的实时速度、位置反馈回上层动态调度层, 构成列车运行时间和节能协同优化反馈闭环结构, 减少列车总延误时间、优化列车运行能耗.

    • 高速列车运行控制与动态调度一体化是指运用先进的感知、传输、控制方法和技术, 实现列车运行优化控制和智能动态调度, 进而深度融合列车运行控制和调度指挥, 实现高铁路网所有列车整体运行效率全局最优. 高速列车运行控制与动态调度一体化通过将高速铁路运营中调度指挥和运行控制的分层模式有机整合起来, 实现列车总延误时间最小、列车节能优化运行并能够准确跟踪优化速度位置曲线等优化目标.

      结合高速列车运行控制和动态调度的业务流程, 高速列车运行控制与动态调度一体化的体系结构如图1所示. 该一体化结构分为三层, 分别是动态调度层(上层)、优化控制层(中层)和跟踪控制层(下层). 动态调度层对应调度集中系统和调度员, 主要对列车延误进行动态调度. 通过调整列车的区间运行时间、列车运行路径、运行次序等, 实现减少列车延误时间、车次数等优化目标, 并将调整后的列车运行图传送给优化控制层. 优化控制层根据调度层给出的运行图, 在给定运行图的约束下实现对列车运行能耗的优化, 并将优化的速度位置曲线传送给跟踪控制层. 跟踪控制层主要目标为跟踪优化控制层所给出的速度位置曲线, 实现对速度位置曲线的有效跟踪. 同时, 高速列车通过铁路综合数字移动通信系统将列车速度、位置信息传输给列车调度指挥系统, 也即将信息反馈给调度层. 将列车的实时速度、位置反馈回上层调度层, 可计算出列车延误时间, 构成反馈闭环, 为列车运行时间调整提供可靠的参考信息. 通过列车运行控制与动态调度一体化反馈优化结构, 实现列车优化、抗扰、准点运行.

      图  1  动态调度和运行控制一体化结构图

      Figure 1.  Integrated structure diagram of dynamic rescheduling and operation control

      本文主要针对该一体化结构的一部分, 即运行时间调整和节能速度位置曲线进行协同优化设计. 高速列车运行时间调整和节能协同优化要求列车运行节能优化、区间运行时间智能最优、节能运行与区间运行时间调整全方位协同. 接下来将建立基于闭塞区间的列车 − 区间模型和列车速度曲线节能优化模型, 并给出目标函数.

    • 考虑由$ N $列高速列车组成的列车群单向运行在有$ M-1 $个闭塞区间的高铁线路上, 将始发站和终点站作为线路的起止点. 列车按发车时间先后顺序分别标号1到$ N $, 闭塞区间按顺序标号1到$ M-1 $, 闭塞区间分界点按顺序标号1到$ M $(始发站视为1号区间分界点, 终点站视为$ M $号区间分界点). 目前, 中国高速铁路运营广泛采用了CTCS-3级列车控制系统, 该系统可获取运行在线路上列车的实时位置和速度信息. 基于该系统, 可以获得列车出清区间的时刻. 因此, 列车$ i,\; i\in{1,2,\cdots,N}$出清区间$ m,\;m\in{1,2,\cdots,M-1} $的实际时刻可以表示为

      $$ t_{m+1}^i = t_{m}^i+r_{m}^i $$ (1)

      其中, $ t_{m}^i $表示第$ i $列列车出清区间$ m-1 $的实际时刻, $ r_{m}^i $则表示列车在闭塞区间$ m $内的实际占用时间. 若列车在区间分界点$ m+1 $不停车, 则列车在闭塞区间$ m $内的实际占用时间$ r_{m}^i $即是列车实际区间运行时间; 若停车, $ r_{m}^i $表示列车实际区间运行时间与车站停车时间之和. 本文假设列车在车站实际停车时间与计划停车时间相同, 不作调整, 因而列车在车站停车时, $ r_{m}^i $也可以表示列车区间运行时间.

      列车$ i $在闭塞区间$ m $内的实际占用时间与计划占用时间关系为

      $$ r_{m}^i = R_{m}^i+u_{m}^i+\omega_{m}^i $$ (2)

      其中, $ R_m^i $表示列车$ i $在闭塞区间$ m $内的计划占用时间; $ u_{m}^i $表示列车$ i $在闭塞区间$ m $内的运行时间调整量; $ \omega_m^i $表示列车$ i $在闭塞区间$ m $运行过程中由于各种原因导致的区间运行时间发生改变, 在时间域上体现为一系列扰动. 该扰动对列车的正点运行造成影响.

      由于本文考虑的是列车晚点运行状态, 优化列车运行时间实质上就是压缩列车区间运行时间, $ u_{m}^i $通常取负值, 表示调度员要求缩短列车$ i $在区间$ m $的运行时间. 受限于线路条件、列车动力系统执行能力和运行安全, 运行时间的调整量$ u_{m}^i $是有限的, 以保证列车能够生成安全可行的速度位置曲线. 因此$ u_{m}^i $满足

      $$ u_{m \min}^i\leq u_{m}^i\leq 0 $$ (3)

      其中, $ u_{m \min}^i $为负值, 表示列车$ i $在第$ m $个闭塞区间运行时间调整量的最小值, 即运行时间调整量的最小值为压缩区间运行时间最大值.

      由式(1)和式(2)可得

      $$ t_{m+1}^i = t_{m}^i+R_{m}^i+u_{m}^i+\omega_{m}^i $$ (4)

      与式(1)相似, 列车$ i $出清区间$ m $的计划时刻可以表示为

      $$ T_{m+1}^i = T_{m}^i+R_{m}^i $$ (5)

      其中, $ T_{m}^i $表示列车$ i $出清区间$ m-1 $的计划时刻.

      将式(5)代入式(4), 可得

      $$ t_{m+1}^i-T_{m+1}^i = t_{m}^i-T_{m}^i+u_{m}^i+\omega_{m}^i $$ (6)

      定义$ e_{m}^i $表示第$ i $列列车出清区间$ m-1 $的实际时刻与计划时刻之差, 则

      $$ e_{m}^i = t_{m}^i-T_{m}^i $$ (7)

      也即列车运行偏离计划时刻表的时间且有$ e_{m}^i\geq 0 $. 因此, 根据式(7), 式(6)可转化为

      $$ e_{m+1}^i = e_{m}^i+u_{m}^i+\omega_{m}^i $$ (8)

      除了时间调整量$ u_{m}^i $的限制条件(3), 其他限制条件有: 为防止列车在运行过程中发生碰撞等事故, 相邻列车出清同一闭塞区间的时间间隔应满足

      $$ t_{m}^i-t_{m}^{i-1}\geq t_{m}^{\min} $$ (9)

      其中, $ t_{m}^{\min} $表示列车运行的最小时间间隔.

      相邻列车出清同一闭塞区间的计划时间间隔为

      $$ H_{m}^i = T_{m}^i-T_{m}^{i-1} $$ (10)

      综合式(7), (9)和(10)可得, 式(9)中限制条件等价于

      $$ e_{m}^{i-1}-e_m^i\leq H_{m}^i-t_{m}^{\min} $$ (11)

      假设高速列车运行在繁忙时段, 每个闭塞区间均被占用. 但由于每一个闭塞区间内仅允许一列高速列车运行, 定义$ t(k) = \left[t_{1}^{k-1},t_{2}^{k-2},\cdots,t_{M}^{k-M}\right]^{\rm T} $, 其中, $ M $表示区间分界点编号, $ k $为列车运行状态编号. 因此, $ t(k) $可以视为每列列车出清闭塞区间时刻的集合. 可以根据$ t(k) $构建高速列车在线路上运行的列车 − 区间模型. 则由式(4)可得, 列车 − 区间模型如下式所示

      $$ t(k+1) = At(k)+R(k)+u(k)+\omega(k)+t_1(k) $$ (12)

      其中,

      $$ \begin{split} & u(k) = \left[u_{0}^k,u_{1}^{k-1},\cdots,u_{M-1}^{k-M+1}\right]^{\rm T}, \; u_{0}^k = 0\\ &R(k) = \left[R_{0}^{k},R_{1}^{k-1},\cdots,R_{M-1}^{k-M+1}\right]^{\rm T}, \; R_{0}^{k} = 0\\ &\omega(k) = \left[\omega_{0}^k,\omega_{1}^{k-1},\cdots,\omega_{M-1}^{k-M+1} \right]^{\rm T},\; \omega_{0}^k = 0\\ &t_1(k) = \left[t_1^k,0,\cdots,0\right]^{\rm T}\\ &A = \left [\begin{matrix} 0&0&\cdots&0\\ 1&0&\cdots&0\\ \vdots&\vdots&\ddots&\vdots\\ 0&\cdots&1&0\end{matrix}\right]_{M\times M} \end{split} $$

      此外, 定义$T(k) =[T_{1}^{k-1},T_{2}^{k-2}, \cdots ,T_{M}^{k-M}]^{\rm T} $, $e(k) = $$ [e_{1}^{k-1},e_{2}^{k-2},\cdots ,e_{M}^{k-M}]^{\rm T} $, $ e_1(k) = [e_1^k,0,\cdots, $$ 0]^{\rm T} $, 可以得到列车 − 区间偏差模型

      $$ e(k+1) = Ae(k)+u(k)+\omega(k)+e_1(k) $$ (13)

      约束条件(3)和(11)可以转化为

      $$ u_{\min}(k)\leq u(k)\leq 0 \qquad\qquad\qquad$$ (14)
      $$ e(k-1)-e(k)\leq H(k)-t_{\min}(k) $$ (15)

      其中,

      $$ \begin{split} & H(k) = \left[H_{1}^{k-1},H_{2}^{k-2},\cdots,H_{M}^{k-M}\right]^{\rm T} \\ &t_{\min} (k) = \left[t_{1}^{\min},t_{2}^{\min},\cdots,t_{M}^{\min}\right]^{\rm T} \\ &u_{\min}(k) = \left[u_{0 \min}^k,u_{1 \min}^{k-1},\cdots,u_{M-1 \min}^{k-M+1}\right]^{\rm T} \end{split} $$

      注1. 本文中所建立的高速铁路列车 − 区间模型是基于城市轨道交通离散事件模型[19-22]扩展而来, 主要有以下不同:

      1)本文中所建立的列车 − 区间模型是基于广义闭塞区间的, 放松了对闭塞区间的限制, 站间可以有任意数目闭塞区间; 而城市轨道交通离散事件模型是站间闭塞的, 要求站站停车, 即每个区间分界点都要停车. 本文的模型对停车区间没作限制, 可以任意指定. 此外, 该模型可以充分利用线路, 提高线路容量.

      2)本文中列车车站停车时间没有具体考虑, 被包含在列车区间占用时间中. 这是因为, 与地铁车站停车时间与客流量相关不同, 高铁列车按座位售票乘车, 乘客数量有限且相对固定. 高速列车车站停车时间在列车延误不是极为严重的情况下是不变的. 因此, 针对本文研究的小范围列车延误, 不对列车车站停车时间进行调整, 将列车车站停车时间整合到区间占用时间中.

      列车运行时间调整的优化目标为尽量减少高速列车的总延误时间, 同时也期望列车运行时间调整量的绝对值尽量小, 即尽可能用较小的运行时间压缩量来最大程度的恢复列车按图运行. 因此, 根据模型(13), 上述优化目标可用目标函数(16)来表征

      $$ J = \sum\limits_{k = k_{0}}^{K_{D}}\left\{e^{\rm T}(k+1)Qe(k+1)+u^{\rm T}(k)Ru(k)\right\} $$ (16)

      其中, $ Q $, $ R $为正定权重矩阵, $ k_{0} $表示系统(13)初始状态编号, $ K_{D} $表示系统(13)终端状态编号. 该目标函数第1部分物理意义为使列车尽量恢复正点运行, 所有列车总晚点之和最小. 第2部分物理意义为使列车运行时间的压缩总量最小, 也即尽量小范围的调整运行时间, 保证旅客的满意度以及满足列车本身运行控制的可执行条件. 权重矩阵$ Q $, $ R $的选取主要与两个优化目标的优先性有关. 若优先使列车恢复正点运行, 则选取一个较大的$ Q $, 目标函数第1部分越小, 列车总晚点之和越少; 若优先最小化列车运行时间压缩量, 则选取一个较大的$ R $, 目标函数第2部分越小, 列车运行时间压缩量越少. 实质上, 优化目标函数(16)是列车总晚点时间和列车运行时间压缩量的综合最优. 因此, 通过最小化目标函数(16), 可使列车尽量恢复正点运行并满足运行时间压缩量尽可能小. 综上所述, 针对线路上运行的所有高速列车, 该优化问题可以描述为

      $$ \begin{split}&\min J=\sum\limits_{k=k_{0}}^{K_{D}}\left\{e^{\rm T}(k+1)Qe(k+1)+u^{\rm T}(k)Ru(k)\right\} \\ &{\rm s.t.}\quad (13)\sim (15) \\[-12pt] \end{split} $$ (17)

      值得注意的是, 目标函数第2部分优化目标为最小化列车运行时间压缩量. 同时, 由于区间运行时间越短, 能耗越大[23], 尽可能少地压缩运行时间可以减少列车运行能耗. 因此, 目标函数第2部分也表示列车运行能耗. 优化目标函数第2部分, 可以节约列车运行能耗. 但优化问题(17)不能得到节能优化速度曲线, 接下来将给出列车速度曲线优化模型, 实现列车运行时间与节能的协同优化.

    • 求解优化问题(17)可以求得每列列车在每个闭塞区间内的运行时间压缩量, 再由列车在每个区间的计划运行时间, 可以得到列车在每个闭塞区间的实际运行时间.

      这里只考虑列车牵引能耗, 则高速列车$ i $速度曲线节能优化问题可描述为[7]

      $$ \min \Pi = \int_{t_{m}^{i}}^{t_{m+1}^{i}}v(t)F(v){\rm d}t $$ (18)

      其中, $ \Pi $表示列车在闭塞区间内的牵引能耗; $ v(t) $表示列车的运行速度; $ F(v) $表示列车的牵引力. 对列车运行的每个区间分别求解优化问题(18), 得到列车在一个区间的优化牵引能耗值. 再将列车在各个区间内的优化牵引能耗相加, 最后将每列车牵引能耗求和, 得到总的优化牵引能耗, 即

      $$ \min \sum\limits_{i = 1}^N\sum\limits_{m = 1}^{M-1}\Pi = \sum\limits_{i = 1}^N\sum\limits_{m = 1}^{M-1}\int_{t_{m}^{i}}^{t_{m+1}^{i}}v(t)F(v){\rm d}t $$ (19)

      注2. 给出列车在每个区间的运行时间后, 将列车全程运行能耗优化问题转化为优化列车在每个区间的能耗, 再进行求和. 同时, 将所有列车运行能耗相加, 即可求得所有列车运行总能耗. 通过以上描述可知, 固定时刻表下的多列车节能优化运行问题可以转化为单列车区间节能优化问题, 然后进行求和. 因此, 在本文中, 多列车的节能优化问题实质上与单列车区间节能优化问题是等价的.

      列车节能速度曲线优化问题实质上为多约束优化问题. 首先, 列车运行符合列车的动力学特征, 即

      $$ M_g\frac{{\rm d}v(t)}{{\rm d}t} = F(v)-B(v)-f_{1}(v)-f_{2}(p) $$ (20)

      其中, $ M_g $表示列车质量; $ B(v) $表示列车的制动力; $ f_{1}(v) $表示列车的基本运行阻力; $ f_{2} (p) $表示列车的附加运行阻力.

      列车的基本运行阻力主要由机械阻力和空气阻力组成, 可以由戴维斯方程表示[24]

      $$ f_{1}(v) = c_{0}+c_{1}v(t)+c_{2}v^2(t) $$ (21)

      其中, $ c_{0} $, $ c_{1} $, $ c_{2} $为阻力系数.

      列车的附加运行阻力主要与线路的坡度、线路中的弯道和隧道有关. 列车的附加运行阻力可由式(22)表示[25]

      $$ f_{2}(p) = G_{i}{\rm {sin}}\theta_{p}(p)+r_{ci}(p)+r_{ti}(p) $$ (22)

      其中, $ G_{i}{\rm {sin}}\theta_{p}(p) $表示由轨道坡度造成的运行阻力, $ G_i $表示列车重力, $ \theta_{p}(p) $表示轨道的坡度; $ r_{ci}(p) $表示由线路中弯道造成的附加阻力; $ r_{ti}(p) $表示由隧道引起的附加阻力.

      其次, 列车要在规定的时间进出闭塞区间, 并完成区间运行的全程, 满足列车运行速度约束, 即

      $$\begin{split}&{S_m} = \displaystyle\int_{t_m^i}^{t_{m + 1}^i} v (t){\rm{d}}t\\ &v(t)\leq V_{m} \\ &v \left(t_{m}^i \right) = v_{d} \left(t_{m}^i\right) \\ &v\left(t_{m+1}^i \right) = v_{d} \left(t_{m+1}^i\right) \end{split} $$ (23)

      其中, $ S_{m} $表示闭塞区间的长度; $ v_{d} (t_{m}^i) $表示列车占用闭塞区间$ m $起始时刻的期望速度; $ v_{d} (t_{m+1}^i) $表示列车出清闭塞区间$ m $时刻的期望速度; $ V_{m} $为区间限速.

      综上所述, 列车节能速度位置曲线优化问题可以描述为

      $$ \begin{split}& \min\Pi = \displaystyle\int_{t_{m}^{i}}^{t_{m+1}^{i}}v(t)F(v){\rm d}t \\ &{\rm s.t.}\; \; (20) \sim (23) \end{split} $$ (24)

      离散化问题(24), 将区间长度$ S_{m} $等分成$ L $个部分, $ L+1 $个离散的数据点, 可得

      $$ \begin{split} &{\min\Pi = \displaystyle\sum\limits_{l = 1}^L {({S_{m,l}} - {S_{m,l - 1}})\times F({v_{l - 1}})} }\\ &{{\rm{s}}.{\rm{t}}.}\quad{v_{l + 1}^2 = v_l^2 + 2({S_{m,l + 1}} - {S_{m,l}})\times}\\ &\quad\quad\quad\quad\quad\ \dfrac{F({v_l}) - B({v_l}) - {f_1}({v_l}) - {f_2}({S_{m,l}})}{{M_g}}\\ &\quad\quad{{v_l} \le {V_m}}\\ &\quad\quad{{v_0} = {v_d}\left(t_m^i\right)}\\ &\quad\quad{{v_L} = {v_d}\left(t_{m + 1}^i\right)}\\ &\quad\quad{{f_1}({v_l}) = {c_0} + {c_1}{v_l} + {c_2}v_l^2}\\ &\quad\quad{f_2}({S_{m,l}}) = {G_i}{\rm {sin}}{\theta _{{S_{m,l}}}}({S_{m,l}}) + \\ &\quad\quad\quad\quad\quad\quad\quad{r_{ci}}({S_{m,l}}) + {r_{ti}}({S_{m,l}})\\ &\quad\quad{\displaystyle\sum\limits_{l = 1}^{L-1} \dfrac{ {{S_{m,l + 1}} - {S_{m,l}}} } {v_l} = t_{m + 1}^i - t_m^i}\\[-15pt] \end{split}$$ (25)

      其中, $ F(v_{l}) $, $ B(v_{l}) $, $ f_{1}(v_{l}) $, $ f_{2}(S_{m,l} ) $分别表示离散化的牵引力、制动力、基本运行阻力、附加运行阻力.

      求解优化问题(25), 就是在给定区间运行时分、起止速度、线路限速、列车基本参数等约束条件下寻求节能的列车速度位置曲线. 接下来, 将根据建立的运行时间优化模型(17)和列车速度曲线节能优化模型(25), 设计列车运行时间调整策略、求解节能优化速度位置曲线, 实现对列车运行时间和节能的协同优化.

    • 本节针对优化问题(17), 采用模型预测控制方法, 设计列车运行调整策略. 由于模型(13)存在扰动, 首先对扰动进行估计. 在$ k-1 $步, 扰动的估计值可以表示为

      $$ \hat{\omega}(k-1) = e(k)-Ae(k-1)-u(k-1)-e_1(k-1) $$ (26)

      因此, 在第$ k $步, 需要对扰动进行补偿. 令

      $$ \omega(k) = \hat{\omega}(k-1) $$ (27)

      相对于以往大部分相关研究[19-20]中直接将扰动认为是已知的情况, 本文通过上述补偿策略, 可以估计列车在运行过程中所受到的扰动影响, 进而在下一步的运行调整中进行补偿.

      设定预测时域和控制时域均为$ K_{f} $, $ k+t|k $表示第$ k $步对第$ k+t $步的预测, 则优化问题(17)转化为基于补偿策略的模型预测控制问题

      $$ \begin{split} & \min J = \displaystyle\sum\limits_{t = 1}^{K_{f}}\{e^{\rm T}(k+t|k)Qe(k+t|k)+\\ &\qquad u^{\rm T}(k+t-1|k)Ru(k+t-1|k)\} \\ &{\rm s.t.}\; \; (26), \; (27)\\ &\qquad e(k+1|k) = Ae(k|k)+u(k|k)+\omega(k)+e_1(k|k)\\ &\qquad e(k+t+1|k) = Ae(k+t|k)+u(k+t|k)+\\ &\qquad \qquad\quad e_1(k+t|k),\; t = 1,2,\cdots,K_{f}-1\\ &\qquad e(k+t|k)\geq 0, \;t = 0,1,\cdots,K_{f}\\ &\qquad u_{\min}(k+t)\leq u(k+t|k)\leq 0\\ &\qquad e(k+t-1|k)\!-\!e(k+t|k)\!\le\! H(k+t)\!-\!t_{\min}(k+t) \end{split}$$ (28)

      采用式(28)基于补偿策略的预测控制方法, 可以解决由于未知随机扰动对列车运行时间调整造成的影响. 接下来, 将基于补偿策略的预测控制问题转化为标准的二次规划问题进行求解.

      定义

      $$\begin{split} E(k) =\;& [e^{\rm T}(k+1|k),e^{\rm T}(k+2|k),\cdots, \\&e^{\rm T}(k+K_{f}|k)]^{\rm T} \end{split} \hspace{60pt}$$
      $$\begin{split} U(k) =\;& [u^{\rm T}(k|k),u^{\rm T}(k+1|k),\cdots,\\&u^{\rm T}(k+K_{f}- 1|k)]^{\rm T} \end{split} \hspace{73pt}$$
      $$ \bar{e}(k) = [e_1^{\rm T}(k|k),e_1^{\rm T}(k+1|k), \cdots,e_1^{\rm T}(k+{K_f}-1|k)]^{\rm T} $$

      对未来$ K_{f} $步的预测可由下式表示:

      $$ E(k) = A_{x}e(k)+B_{u}U(k)+D_{d}\omega(k)+B_{u}\bar{e}(k) $$ (29)

      其中,

      $$ \begin{array}{l} A_{x} = [A,A^2,\cdots,A^{K_{f}}]^{\rm T} \\ B_{u} = \left [\begin{matrix}I&0&\cdots&0\\A&I&0&0\\ \vdots&\vdots&\ddots&0\\A^{K_{f}-1}&\cdots&A&I\end{matrix}\right] \\ D_{d} = [I,A,\cdots,A^{K_f-1}]^{\rm T} \end{array} $$

      式(28)中目标函数转化为

      $$ \min J = E^{\rm T}(k)\bar{Q}E(k)+U^{\rm T}(k)\bar{R}U(k) $$ (30)

      其中, $ \bar{Q} = {\rm {diag}}\{Q,Q,\cdots,Q\}_{MK_{f}\times MK_{f}} $, $ \bar{R} = {\rm {diag}}\{R, $$R,\cdots,R\}_{MK_{f}\times MK_{f}} $.

      列车运行时间调整量约束(14)可以转化为

      $$ \left[ {\begin{array}{*{20}{c}} { - T}\\ T \end{array}} \right]\times U(k) \ge \left[ {\begin{array}{*{20}{c}} 0\\ \vdots \\ 0\\ {{u_{\min }}(k)}\\ \vdots \\ {{u_{\min }}(k + {K_f} - 1)} \end{array}} \right]$$ (31)

      其中, $ T = {\rm {diag}}\{I,I,\cdots,I\}_{MK_{f}\times MK_{f}} $.

      相邻列车间的运行间隔约束(15)可以转化为

      $$ \begin{array}{l} {T_1}e(k) - {T_2}E(k) \le \\ \qquad\quad\left[ {\begin{array}{*{20}{c}} {H(k + 1) - {t_{\min }}(k + 1)}\\ \vdots \\ {H(k + {K_f}) - {t_{\min }}(k + {K_f})} \end{array}} \right] \end{array} $$ (32)

      其中,

      $T_1 = [I,0,\cdots,0]_{M\times MK_{f}}^{\rm T}$

      $$ T_{2} = \left [\begin{matrix}I&0&\cdots&0\\ -I&I&\cdots&0\\ \vdots&\vdots&\ddots&\vdots\\ 0&\cdots&-I&I\end{matrix}\right]_{MK_{f}\times MK_{f}}\qquad\quad $$

      根据式(29), 上式等价为

      $$ \begin{array}{l} {T_2}{B_u}U(k) \ge ({T_1} - {T_2}{A_x})e(k) - {T_2}{D_d}\omega (k) - \\ \qquad{T_2}{B_u}\bar e(k) - \left[ \!\!\!\!{\begin{array}{*{20}{c}} {H(k + 1) - {t_{\min }}(k + 1)}\\ \vdots \\ {H(k + {K_f}) - {t_{\min }}(k + {K_f})} \end{array}}\!\!\! \right] \end{array} $$ (33)

      根据式(29)和式(30), 可得

      $$ \begin{aligned} J = \;&{{{[{A_x}e(k) + {B_u}U(k) + {D_d}\omega (k) + {B_u}\bar e(k)]}^{\rm{T}}}\bar Q}\times\\ &{[{A_x}e(k) + {B_u}U(k) + {D_d}\omega (k) + {B_u}\bar e(k)] + }\\ &{{U^{\rm{T}}}(k)\bar RU(k)}=\\ & {{U^{\rm{T}}}(k)\bar RU(k) + {{[{B_u}U(k)]}^{\rm{T}}}\bar Q[{B_u}U(k)] + 2}\times \\ &{{{[{A_x}e(k) + {D_d}\omega (k) + {B_u}\bar e(k)]}^{\rm{T}}}\bar Q[{B_u}U(k)] + \Xi } \end{aligned}$$ (34)

      其中,

      $$\begin{split} \Xi =\;& [A_{x}e(k)]^{\rm T}\bar{Q}[A_{x}e(k)]+2[D_d\omega(k)]^{\rm T} \times \\&\bar{Q}[A_{x}e(k)]+[D_{d} \omega(k)]^{\rm T}\bar{Q}[D_{d}\omega(k)]+[B_u\bar{e}(k)]^{\rm T}\times \\& \bar{Q}[B_u\bar{e}(k)] +2[D_d\omega(k)]^{\rm T} \bar{Q}[B_u\bar{e}(k)]+\\& 2[A_xe(k)]^{\rm T}\bar{Q}[B_u\bar{e}(k)] \end{split}$$

      为已知. 因此, 优化问题(28)转化为

      $$ \begin{split}&\min J=\,U^{\rm T}(k)\bar{R}U(k)+[B_{u}U(k)]^{\rm T}\bar{Q}[B_{u}U(k)]+\\ &\qquad\quad\;\;\; 2[A_{x}e(k)+D_{d}\omega(k)+B_u\bar{e}(k)]^{\rm T}\bar{Q}[B_{u}U(k)] \\ &{\rm s.t.}\; \; (26), \; (27), \; (29), \; (31), \; (33) \\[-10pt] \end{split}$$ (35)

      由式(35)可知, 多列车运行时间调整问题转化为二次规划问题. 根据预测控制的基本原理, 得到的控制序列的第1步将作用于系统. 因此, 求解式(35), 多列车在区间内的运行时间调整量为

      $$ u(k) = [I,0,\cdots,0]U(k) $$ (36)

      重复进行上述过程, 每一次预测都进行计算求解并得到运行时间调整量, 可得运行时间调整量序列. 进而计算出列车调整后的区间运行时间.

    • 根据第2.1节优化后的列车区间运行时间, 线路限速、线路长度、运行时间、坡度以及列车进入区间的初速度限制、末速度限制等限制条件, 求解列车速度曲线节能优化模型, 得到每列车在每个区间内的节能优化速度位置曲线. 一般情况下, 列车节能运行策略将列车运行过程分为最大牵引、巡航、惰行、最大制动四个阶段. 根据区间运行情况, 也可能只包含部分工况. 基于上述研究, 本文采用二分法和数值分析方法计算给定运行时间、给定区间初速度末速度限制下的列车优化速度位置曲线. 具体方案如算法1所示.

      算法1.

      步骤1. 输入区间基本条件: 运行时间约束、速度约束、线路路况、列车基本参数等条件; 给定能耗迭代初值, 包括列车区间能耗设定最小值$ \Pi_{\min} $, 列车区间能耗设定最大值$ \Pi_{\max} $.

      步骤2. 计算$ \Pi_{\rm{mid}} = (\Pi_{\min}+\Pi_{\max})/2 $.

      步骤3. 求出在能耗为$ \Pi_{\rm{mid}} $条件下列车运行时间$ T(\Pi_{\rm{mid}}) $和最优速度位置曲线. $ T(\Pi_{\rm{mid}}) $和最优速度位置曲线可利用算法2中给出的方法求得.

      步骤4. 将运行时间$ T(\Pi_{\rm{mid}}) $与给定列车运行时间相比较. 若$ T(\Pi_{\rm{mid}}) $大于给定列车运行时间, 则用$ \Pi_{\rm{mid}} $代替$ \Pi_{\min} $, $ \Pi_{\max} $不变, 返回步骤2; 若$ T(\Pi_{\rm{mid}}) $小于给定列车运行时间, 则用$ \Pi_{\rm{mid}} $代替$ \Pi_{\max} $, $ \Pi_{\min} $不变, 返回步骤2; 若$ T(\Pi_{\rm{mid}}) $等于给定列车运行时间, 则$ \Pi_{\min} $为所求解的优化能耗, 同时得到最优速度位置曲线.

      注3. 这里$ T(\Pi_{\rm{mid}}) $与给定列车运行时间偏差小于一个极小的时间, 即认为$ T(\Pi_{\rm{mid}}) $与给定列车运行时间相同. 算法1步骤3中$ T(\Pi_{\rm{mid}}) $和最优速度位置曲线可由文献[23]给出的方法获得. 本文根据高速列车节能优化运行实际情况对其进行应用, 具体内容如算法2所示.

      算法2.

      步骤1. 将区间$ S_{m} $划分为$ L $个部分, 即将速度位置曲线分为$ L+1 $个离散的数据点. 构成速度序列: $ \{v_{l}:0\!\leq \!l\!\leq \!L,l\!\in \!{\bf N}\} $, 位置序列: $ \{S_{m,l}:0\!\leq\! l\!\leq\! L,l\!\in\! {\bf N}\} $. 令$ l = 0 $, $ v_{l} = v_{d}(t_{m}^i) $.

      步骤2. 令$ l = l+1 $, 计算列车加速度、速度和能耗剩余值: $ \Pi_{\rm{mid}} = \Pi_{\rm{mid}}-F(v_{l-1} )\times (S_{m,l}-S_{m,l-1}) $, 列车处于最大牵引状态.

      步骤3. 若$ v_{l}<V_{k} $$ \Pi_{\rm{mid}}>0 $, 返回步骤2; 否则, 若$ \Pi_{\rm{mid}}>0 $, $ v_{l+1} = v_{l} $, $ \Pi_{\rm{mid}} = \Pi_{\rm{mid}}-[f_{2}(S_{m,l})+ $$ f_{1}(v_{l})]\times (S_{m,l+1}-S_{m,l}) $, 列车处于巡航状态. 若$ \Pi_{\rm{mid}}>0 $, 令$ l = l+1 $, 重复步骤3.

      步骤4. 若$ \Pi_{\rm{mid}} = 0 $, $ l<L $, 列车进入惰行状态, $ v_{l+1}^2 = v_{l}^2-2\times[f_{1}(v_{l})+f_{2}(S_{m,l})]\times(S_{m,l+1}-$$ S_{m,l})/M_g. $$ l = l+1, $ 重复步骤4.

      步骤5. 令$ j = L $, $ v_{j} = v_{d}(t_{m+1}^i ) $, 如果$ j>l $, 反向计算刹车曲线: $ v_{j-1}^2 = v_{j}^2+2\times[B(v_{j})+f_{1}(v_{j})+ $$f_{2} (S_{m,j})]\times (S_{m,j}-S_{m,j-1})/M_g $. 令$ j = j-1 $, 重复步骤5, 直到与惰行曲线交汇, 得到刹车曲线与惰行曲线交汇点.

      步骤6. 输出速度位置曲线, 列车区间运行时间: $ T(\Pi_{\rm{mid}}) = \displaystyle\sum\nolimits_{l = 1}^{L-1}[S_{m,l+1}-S_{m,l}]/v_{l} $.

      本文中采取对每列列车在每个运行区间逐个计算节能优化速度位置曲线策略, 以达到高铁线路上所有列车运行总能耗优化的目标. 本文采用的分层优化策略, 将列车运行时间与节能协同优化这一多目标优化问题进行分解. 首先优化列车运行时间, 使列车尽可能恢复正点运行; 其次, 在恢复正点运行的基础上优化列车速度位置曲线.

    • 作为高速列车运行控制和动态调度一体化控制结构的三个核心层的底层, 跟踪控制是根据协同优化后生成的节能速度位置曲线, 实现速度位置曲线的准确跟踪. 在一体化整体框架中, 列车跟踪控制器所起的作用相当于具备自主抗扰能力的执行器. 这部分不是本文的主要贡献, 为保持本文所涉及的一体化框架的完整性和有助于理解仿真实验的结果, 本节对列车跟踪速度位置曲线问题作简要介绍.

      针对高速列车跟踪给定速度位置曲线问题, 大量学者进行了深入研究, 并将一系列的控制理论和方法应用到跟踪控制中, 如模糊控制[26]、自适应控制[27]、迭代学习控制[28]、鲁棒控制[29]等.

      根据第2.2节生成的列车优化速度位置曲线以及区间运行时间, 采用文献[30-31]中给出的方法设计鲁棒控制器跟踪优化速度位置曲线, 这里不具体介绍设计过程, 只简要给出结论. 跟踪控制器形式为

      $$ \bar{u}(t) = \bar{K}x(t) $$ (37)

      其中, $ x(t) $表示列车状态, $ \bar{K} $为控制器参数. 通过跟踪控制器, 可实现对速度位置曲线的跟踪.

    • 为了验证本文所提出的高速列车运行时间与节能协同优化策略的有效性, 仿真分析选择京津城际高速铁路. 京津城际高速铁路全线共设北京南站、亦庄站、永乐站、武清站、天津站5座车站, 全程120 km. 仿真实验中, 将北京南到亦庄间分为2个闭塞区间, 亦庄站到永乐站之间分为2个闭塞区间, 永乐站到武清站之间分为2个闭塞区间, 武清站到天津站之间分为3个闭塞区间, 即全线共分为9个闭塞区间, 共计10个闭塞区间分界点. 实验中假设列车从天津站出发上行至北京南站. 为了便于描述, 线路途中8个区间分界点编号2到9, 天津站编号为1, 北京南站编号为10.

    • 根据京津城际高速铁路实际运营中在繁忙时段列车运行间隔为5分钟这一实际情况. 设定相邻列车在每个区间分界点的计划运行间隔和最小运行间隔分别为$ H_m^i = 300\;{\rm{s}}$, $ t_m^{\min} = 180\;{\rm{s}}$. 预测时域和控制时域$ K_f = 5 $. 列车区间运行时间调整量最小值为$-30\,{\rm s} $, 即$ u_{m \min}^i = -30\;{\rm{s}} $. 分别选取$ Q = 0.7I $, $ R = $$ 0.4I $; $Q = $$ 0.12I $, $ R = 3I $以及$ Q = 0.07I $, $ R = 0.4I $进行3组实验. 第1组实验选取列车11为例, 第2组实验选取列车13, 第3组实验选取列车16. 设定列车11, 13, 16在天津站发车时初始运行偏差分别为100 s, 100 s, 260 s. 列车在运行过程中扰动分别在区间2和区间6对列车运行造成5 ~ 20 s随机延误. 求解式(35)中优化问题, 得到线路上运行所有列车运行时间调整结果. 根据模型预测控制基本原理, 这里只取控制序列中的第1步作为调整结果并应用. 之后重复进行上述过程, 得到列车运行时间调整序列. 列车到达每个区间分界点实际时间与计划时间的偏差如图2$\sim $4所示.

      图  2  列车11在区间分界点延误时间

      Figure 2.  Delay time of train 11 at section demarcation point

      图  4  列车16在区间分界点延误时间

      Figure 4.  Delay time of train 16 at section demarcation point

      图2$\sim $4可知, 列车在运行过程中, 虽然受外部扰动的影响, 但列车在逐渐恢复正点运行, 当列车从天津出发到北京南站, 列车基本能够恢复正点运行. 图2图3中, 列车11、13初始延误相同, 但由于两组实验选取权重矩阵不同, 优化结果不同. 图2中, 选取了较大的$ Q $, 列车11优先恢复正点运行. 图3中, 选取了较大的$ R $, 列车13优先优化运行时间压缩量. 而图4中, 列车16初始晚点时间为260 s, 相对于京津城际全程运行仅需35分钟左右且列车运行频繁的实际情况, 已经是较大的晚点, 经过运行时间调整后, 列车延误时间逐渐在减小, 到达终点站时已仅有20 s延误, 基本恢复正点到站.

      图  3  列车13在区间分界点延误时间

      Figure 3.  Delay time of train 13 at section demarcation point

      同时, 本文设计的运行时间调整策略可以在每个闭塞区间都压缩部分运行时间, 而不是在一个区间内大幅度压缩运行时间. 采用本文中的运行时间调整策略可以保证列车的执行能力, 也可以节约能耗, 更加有利于下一步设计列车优化速度位置曲线来进一步优化能耗. 对于经停武清站的列车, 由于京津城际高速列车只有少数列车经停武清站, 且列车至武清站时延误时间相对于初始列车延误时间已经减少. 因此, 本文设计的列车运行时间调整策略能够减少列车延误时间. 但由于本文优化目标是综合考虑列车晚点时间和运行时间压缩量, 因此, 对于经停武清站的列车, 列车到达武清站时, 延误时间减小, 但不能完全恢复正点到达武清站.

    • 本小节基于第3.1节给出的列车运行调整时间, 以列车13为例, 生成速度位置曲线. 由于实验过程相同, 这里仅选取一个区间进行仿真实验. 假设列车从天津站至区间分界点2正常运行时间为280 s, 经运行调整后运行时间为252 s. 运行距离13 km, 列车区间限速250 km/h, 初速度为0, 末速度为230 km/h. 列车选取CRH-3型高速列车为例, 列车质量为$ 4.08\times 10^5\;{\rm{kg}}$. 列车在线路上运行的平均附加运行阻力取2.4 kN. 列车最大牵引力(kN)和最大制动力(kN)可按下式计算[7]:

      $$ F(v) =\! \left\{ {\begin{aligned} &{300 - 0.284v},&\!\!{0 \le v \le 119.7\;{\rm{km/h}}\quad\qquad\;\;}\\ &{266\times \frac{{119.7}}{v}},&\!\!{119.7\;{\rm{km/h}} < v \le 300\;{\rm{km/h}}} \end{aligned}} \right. $$
      $$ B(v) =\!\left\{ {\begin{aligned} &{300 - 0.281v},&\!\!{0 \le v \le 106.7\;{\rm{km/h}}}\quad\qquad\;\;\\ &{270\times \frac{{106.7}}{v}},&\!\!{106.7\;{\rm{km/h}} < v \le 300\;{\rm{km/h}}} \end{aligned}} \right.$$

      列车基本运行阻力为

      $$ f_1(v) = 6.7744+0.05719v+0.0008235v^2\; \rm{kN} $$

      设初始$ \Pi_{\min} = 1.5\times 10^6\;{\rm kJ}, $$ \Pi_{\max} = 2.5\ \times$$10^6\,{\rm kJ} $. 采用算法1和算法2求解式(25)中优化问题, 求得优化运行能耗为$ 1.672\times 10^6 $ kJ. 得到列车13从天津站到区间分界点2优化速度位置曲线如图5所示.

      图  5  列车13优化速度位置曲线

      Figure 5.  The optimal speed-position trajectory of train 13

      为了保持一体化结构的完整性, 图6给出了跟踪控制器的跟踪效果. 由实验结果可知, 跟踪控制器能够很好地跟踪优化位置速度曲线. 由于列车运行距离较长, 图6只能表示出跟踪的整体效果. 接下来, 对图6的部分细节进行放大处理, 如图7所示. 图7给出了跟踪效果的局部放大图, 以运行过程中11 km到13 km为例, 包括部分巡航、惰行以及制动工况. 由图7可以看出, 相对于高速行驶的列车, 总体上还是能够很好地跟踪期望的速度位置曲线.

      图  6  速度位置曲线跟踪效果图

      Figure 6.  Tracking performance of speed-position trajectory

      图  7  跟踪效果局部放大图

      Figure 7.  Partial enlarged detail of tracking performance

    • 本文针对小范围延误情况下的高速列车运行时间调整和节能运行问题, 在运行控制和动态调度一体化的框架下, 提出了一种基于闭塞区间的高速列车运行时间和节能协同优化方法. 基于列车 − 区间模型和列车速度曲线节能优化模型, 并结合分层优化和设定值优化的思想, 设计了列车运行时间调整和节能运行协同优化策略. 该协同优化方法能够实现对列车小范围延误的运行时间调整, 减小列车延误时间, 恢复列车正常运行, 实现高速列车“按秒控制, 按秒调整”的运营需求, 并通过优化速度位置曲线节约列车运行能耗. 在现有方案中, 动态调度与速度位置曲线节能优化还没有实现真正意义上的联合优化, 所设计的列车 − 区间模型应用场景有限. 在未来的研究中, 在计算复杂度允许的情况下, 将考虑更加符合高铁实际运营的突发事件下列车运行控制与动态调度的多目标联合优化.

参考文献 (31)

目录

    /

    返回文章
    返回