城市轨道交通能量运行控制系统摘 要：城市轨道交通和干线铁路的高速增长带来的能耗、碳排放等问题愈发引起社会的关注为提升城市轨道交通系统能量利用效率，文章提出一种用于城市轨道交通牵引供电网的能量运行控制系统该系统采用 4G 通信获取每辆列车的瞬时功率和位置、牵引供电系统的参数以及变电站和双向变流器的实时状态信息，依此迭代计算网络空间矩阵，完成牵引供电网络的数字重构基于重构的牵引供电网络，能量运行控制系统根据目标函数以及系统物理约束，采用智能算法寻优牵引变电站直流输出电压指令，并下发给全线路双向变流器，实现整条地铁线路的潮流分配，从而使地铁牵引电力系统具有更好的供电质量和更高的效率通过对采用该能量运行控制系统的某城市地铁 4 号线进行现场测试，结果显示系统总体效率提升了 0.767%，验证了所提城市轨道交通能量运行控制系统的有效性 本文源自李红波; 张超; 黄子昊; 王雄; 罗文广; 张志学， 控制与信息技术 发表时间：2022-06-07 关键词：轨道交通牵引供电;能量运行控制;系统节能;负荷预测;运行图优化;数字重构 0 引言 2022 年 12 月获批的?巴黎协定?其中心目标是将温度上升限制在比工业化前水平高 2 ℃以下，并努力将温度上升限制在比工业化前水平高 1.5 ℃ [1] 。
公路运输在温室气体排放中占有主要地位 [2]要减少 80%的总体碳足迹，就必须减少 95%的公路运输 [3] ;然而，经济增长需要加强运输能力，这将导致温室气体排放量上升针对减少温室气体排放和增加运输这两种要求相互矛盾的难题，城市电气化轨道交通由于具备良好的零碳排放和较强的运输能力，其重要性日益凸显 [4] 现有城市轨道交通电网通常采用多脉波(12 或 24 脉波)整流器，其将 35 kV 或 10 kV 的中压交流电转换为 750 V 或 1 500 V 的直流电给列车供电 [5] 由于二极管整流器具有单向工作特性，再生制动能量无法被回收，导致直流母线电压升高，直接危害列车以及车载电气设备平安运行 在过去的几十年中，轨道交通系统通常通过配置额外的制动电阻器来耗散多余的再生制动能量，以防止直流母线电压过冲城市电气轨道交通是兆瓦级系统，再生制动能量巨大，所以制动电阻产生的热量会导致地铁隧道和车站的温度升高 [6]制动电阻器解决方案不仅浪费再生制动能量，同时增加了环境控制系统的能耗 为了解决这个问题，业界开展了再生制动能量回收技术的研究再生制动能量回收技术可以分为两种：一种是通过能量回馈逆变器将再生制动能量回馈至 35 kV 电网，以用于其他负载，例如电站中的辅助设备[7];另一种是通过储能系统回收再生制动能量[8-9] 。
再生制动能量回馈逆变器具有两种运行状态：当线路制动能量小于逆变器的容量时，能量回馈逆变器将所有制动能量有效地引入直流母线电压;当制动能量大于逆变器的容量时，能量回馈逆变器将以额定功率作为恒定电源运行，尽可能多地回收制动能量并抑制直流母线中的过电压但如果能量回馈逆变器的电压指令设置不正确，由于相邻的列车无法充分吸收该能量，那么会导致再生能量的利用效率降低此外，在某些工况下，回收的能量被回馈至 110 kV 公共电网，不能被轨道交通系统负荷利用 [7]这种复杂的潮流增加了电网失稳的风险，且较长的传输路径会增加功率损耗并降低效率而能量存储方案因为所有制动能量都被就近存储，防止了由复杂的潮流引起的功率损耗，因此具有比能量回馈逆变器更高的能量利用率;但再生制动功率和列车起动、加速等运行状态的功率曲线会呈现顶峰值的脉冲形状为了使储能系统具有较快的响应速度，应采用超级电容器或飞轮等功率型储能元件;但这些低能量密度的存储元件将增加系统体积和本钱 [10]另外，城轨交通供电系统目前采用电力数据采集和监控系统(power supervisory control and data acquisition，PSCADA)，只对电网局部状态进行数据监测及界面数据、图形显示、报表等简单数据处理 [11]，不具备利用数据对整个供电网络进行潮流优化调度的功能。
针对上述问题，本文提出一种新型的城市轨道交通智能牵引供电能量运行控制系统(简称“运控系统〞)与传统轨道交通牵引供电系统不同，本文牵引供电系统中的双向脉冲整流器和能量回馈逆变器被双向变换器代替 [12]另外，通过配置车载和地面 4G 无线通信模块(简称“4G 模块〞)，在列车、双向变流器和能量运控系统(ECMS)之间建立通信环路能量运控系统通过 4G 模块收集列车、双向变流器、变压器、牵引变电站以及牵引线和钢轨的参数，对双向变流器的输出电压指令进行优化，最终到达降低电能传输损耗、抑制列车电压波动、降低钢轨电位的目的 1 牵引供电能量运控系统架构 牵引供电系统主要包括基于双向变流器的变电站、列车、地面和车载 4G 模块以及能量运控系统(图 1)本文所提能量运控系统(ECMS)是可逆型智慧地铁牵引电网的上层控制单元 4G 模块被分别安装在双向变流器、列车和能量运控系统内部，在变流器和能量运控系统之间进行相互数据通信出于系统平安性考虑，列车仅向能量运控系统发送状态量，但不接收能量运控系统指令，因此列车必须按照规划的运行图行驶，在地铁运营期间无法通过能量运控系统对运行图进行实时调整。
能源运控单元是能量运控系统的核心，包括效劳器、数据库、4G 通信模块、监视器、键盘和其他外围输入及输出设备，提供数据发送和接收、数据分析以及控制策略制定的功能根据目标函数和系统物理约束条件，能量运控系统通过智能算法确定直流链的功率分配，以实现节能、稳定电压和优化线路质量的目标 作为物理执行单元，双向变流器跟踪来自能量运控系统的输出电压命令，以完成列车和牵引电源系统之间的实际潮流分配 牵引变电站的位置、直流输出电压、电流、列车功率和运行位置以及牵引网导线阻抗、钢轨阻抗等参数信息通过 4G 模块传输至能量运控系统根据接收到的数据，能量运控系统迭代计算牵引供电系统的导纳矩阵、节点电压和节点电流，以完成实际牵引供电系统的数字化重构然后在确保牵引供电系统和列车平安稳定运行的前提下，考虑输电损耗、交流 110 kV 电网功率因数以及变电站容量、电网电压波动和钢轨电位等约束条件，通过启发式算法搜索直流输出电压、无功和有功电流的最优指令另外，能量运控系统可以动态调整牵引网全线双向变流器的输出特性，灵活地分配各个双向变流器的输入、输出功率，从而实现对临近故障变电站的能量救援，以确保牵引供电系统的供电可靠性。
2 能量运控系统工作原理 牵引供电能量运控系统优化调度策略有两个关键：(1)每个采样周期内重建真实的供电网络;(2)针对多个目标和约束条件，能量运控系统通过遗传算法求解双向变流器的直流输出电压 2.1 牵引供电系统的数字重构 2.1.1 系统建模 由于每个牵引变电站都采用双向变流器与列车交换能量，为方便表述，本文用双向变流器替代牵引变电站作为牵引供电系统模型中的电压源或功率源节点 [14] 牵引供电网的等效电路如图 2 所示，其中 Vx， rcx 和 (x=1, 2, …, n) 分别表示第 x 个双向变流器的直流输出电压、电流和等效输出阻抗;VTx 表示第 x 个列车的直流输入电压;Rx 表示第 x 个和第 (x+1) 个双向变流器之间传输线的阻抗整个运行过程中，双向变流器的状态函数 s 为(逆变 / 整流模式) (1) 式中：Pmax——双向变流器峰值功率，其中 Pmax_T 表示整流模式下峰值功率，Pmax_F 表示逆变模式下峰值功率;PLoad——列车功率 当双向变流器的最大功率大于负载功率需求时，将变流器等效为电压源;而当双向变流器的最大功率小于负载功率需求时，将变流器等效为恒功率源。
两种模式的功率输出方式不同，Pmax_F是Pmax_T的两倍，这是因为在工作周期内，在逆变模式下，双向变流器先运行 30 s 回收制动能量，然后再停止运行 90 s;而在整流模式下，双向变流器连续输出功率 近年来，随着车载制动电阻器的取消，线路上运行的列车可被等效为一个时变功率源 2.1.2 牵引供电系统节点电压和潮流计算 牵引供电系统节点电压计算流程如下： (1)根据牵引线和钢轨的阻抗、每列列车的实时位置以及每个变电站的位置来制定节点导纳矩阵，然后设置每个双向变流器和每个列车的初始电压值 (2)利用 Picard 算法计算节点电压，并将其用于更新双向变流器和列车的状态如果迭代值收敛，那么迭代过程完成;否那么，再更新来自列车的数据并修改节点导纳矩阵之后，计算程序将跳至 Picard 迭代步骤，并重复之前的步骤图 3 示出节点电压计算具体流程 为方便定义节点导纳矩阵，图 4 示出牵引供电系统的简化网络拓扑 导纳矩阵的节点根据距地铁线路始发站的距离进行升序排列由于列车的移动性，通常需要在下一个采样周期内对导纳矩阵的节点进行重新排序以形成新矩阵而双向变流器的位置不变，因此不必每次都重新配置导纳矩阵中双向变流器的节点位置。
于是，可依次按双向变流器、上行列车、下行列车顺序排布，构建节点导纳矩阵 Y： (2) 式中：Yss——双向变流器之间的节点导纳矩阵; Ytt——列车之间的节点导纳矩阵;Yst——双向变流器和列车之间的节点导纳矩阵;Yts——列车与双向变流器之间的节点导纳矩阵，由于矩阵的对称性，因此Yts 和 Yst 相同 以 Yss 为例说明计算过程假设无列车在地铁线上行驶，那么每个供电区间的导纳矩阵 Yss 中的牵引线导纳(g i T)、钢轨导纳(g i R)和泄露导纳(g i X)依次按式 (3) 计算获得 (3) 式中：li——第 i 个电源局部的长度;r1, r2——牵引线缆和钢轨的单位电阻;g0——钢轨的单位泄漏导纳 考虑第 i 个电源局部中的上行列车的情况，应将第 i 个双向变流器和第(i + 1)个双向变流器对应的节点导纳修改为 (4) 同样，如果第 i 个电源 s 中有一个下行列车，那么 Yss 也可以根据式 (4) 进行修正Yss 和 Yst 需要根据每列列车的实时位置进行更新Yss 中的对角元素可以由 Ytt 和 Yst 获得，而无须额外进行计算 在获得节点导纳矩阵之后，可以通过 Picard 法来计算整个系统节点电压和系统潮流，从而完成牵引供电系统的数字重构。
2.2 基于双向变流器的能量管理策略 基于双向变流器的能量管理策略旨在通过动态调整双向变流器的直流输出电压，在保证整个地铁线路正常运行的前提下尽可能地降低牵引电网能耗 牵引电网能耗的目标函数如下： (5) 式中：Vi——第 i 个双向变流器的直流输出电压; Ii——第 i 个双向变流器的直流输出电流;PBC， PBC_rated——双向变流器的实时功率和额定功率;Vtrain——列车电压;Vsteel_rail——节点轨电位;n——地铁线路中变电站的总数 根据标准 EN50122-1:1998?铁路应用——地面装置?，轨电位不得高于 110 V另外，制开工况下，当 Vtrain>1 900 V 时，列车将启动闸瓦制动，以减少再生制动能量，阻止电压继续升高;但此时，列车制动能量将被转化为热能耗散到空气中，造成能量浪费因此在优化目标函数中，增加列车电压不超过 1 900。