城市轨道交通直流牵引供电系统动态仿真

1、城市轨道交通直流牵引供电系统动态仿真1 引言城轨交通在现代城市中对缓解交通运输压力方面有着不可 替代的作用。 城市轨道交通牵引供电系统的正常运转是保障城轨 列车以及整个城轨系统正常且安全运行的基础。 因而在设计供电 系统之初就需要对其进行潮流分析， 以选择或校验牵引变电所的 位置和相关设备的额定容量， 同时对列车运行时供电系统能否提 供所需功率等进行检验。对直流潮流进行解析的方法不断发展， 从算法上： 国内从最 初的平均运量法、 列车运行图法 1-3 进行估算到引进节点法 4 、 5 、梯形法 6 等应用精确的电路理论进行计算。从模型上：随 着整流技术的不断发展， 整流机组的模型从简单的单折线

2、模型发 展到了如今的 24 脉波整流的多折线模型 7 ，列车也从电流源模 型逐渐调整为更为精确的功率源模型。 另外， 考虑到钢轨对大地 的漏电阻，有些模型采用的是上下行割裂计算 8 ，有些模型将 上下行钢轨合二为一 9 ，这都降低了模型的精度。本文中将应用并改进相关计算模型并应用较为精确的数学 模型对直流牵引供电系统进行动态的仿真研究， 以期得到牵引供 电系统的电气量动态变化结果为进一步的分析与研究提供更为 精确的计算基础数据。2 直流牵引供电系统各部分模型在牵引供电系统中的直流侧是直接为电动机车提供电能的 部分。主要包括轨道，接触网 （或第三轨 ） ，馈电线和回流线。另 外，地铁列车的运行分

3、为上行和下行，上行和下行的接触网（ 或第三轨 ）在变电所母线处是连在一起的，为了建模方便通常认为 上行和下行轨道的回流线仅在变电所处连通。 如图 1 所示，牵引 变电站 （Tss） 将高压的交流电降压整流成 750V 或其它的电压等 级的直流电，通过第三轨或接触网供给列车，然后经轨道回流。图 1 直流侧牵引供电系统示意图2.1 整流机组的多折线外特性模型整流机组的外特性是指， 直流输出端电压随其负荷电流或短 路电流变化的关系曲线。 整流机组在直流牵引供电系统仿真中常 处理成带内阻的恒定理想电压源。 实际上， 整流机组的外特性是 呈非线性的。 本文采用文献 7 中的 2 个并联 12脉波整流电路

4、的 输出特性曲线等效 24 脉波整流电路的外特性曲线。2.2 列车建模由于列车所受取的网功率与网压波动无关， 所受取的牵引网 电流却受网压波动的影响较大， 所以列车模型应该采用功率源模 型更为精确。 列车在某一时刻的功率可由牵引计算和列车运行图 得到。然而，由于列车采用功率源模型，电路求解方程变为非线 性方程组，需要通过迭代求解。2.3 牵引网建模 由于牵引电流并非全部沿走行轨流回牵引变电所， 而是有一 部分泄漏入大地， 再由大地流回钢轨并回到牵引变电所。 计算等 值电阻时， 必须考虑走行轨对大地的泄漏电阻， 按均匀传输线理 论，对每一段走行轨都采用精确等值型电路 10 。这里考虑建模 的精确

5、行，对上下行钢轨应该分开建模如图 2 所示。图 2 直流牵引供电网模型3 算法流程3.1 动态仿真预处理 通过导入的列车牵引计算结果和对运行参数的设置，可获 得系统进行动态仿真时所需的列车每时刻位置、 功率等信息并进 行处理，从而得到当前时刻牵引网的等效电路。3.2 算法整体流程 算法流程如图 3 所示。图 3 直流牵引供电动态求解流程图（1）根据当前时刻牵引网状态和列车位置生成牵引等效电 路；（2）设定整流机组初始工作区间或更新工作区间；（3）设定各列车初始电压；（4）根据牵引网当前的等效电路生成相应的节点导纳矩阵；（5）对得到的节点导纳矩阵进行 LU 分解；（6）利用节点电压方程 YV=I

6、,以及（5）中得到的LU分解 结果计算出新的节点电压， 进而得到各列车的电压更新值并进行 收敛判断：其中为第 i 辆车在 t 时刻的电压值， e 为设定的误差限， 若每辆车都满足收敛条件那么该时刻的LU计算迭代结束，否则更新各列车电压，继续该步循环；（7）各列车电压收敛后计算各变电所电流，并判断各变电 所的工作区间是否出现变化，如果出现变化进入第（2）步修正整流机组参数， 否则该时刻的迭代结束， 输出该时刻各个变电所 的电气量变化；（ 8）更新时间， 返回步骤 （1），进入下一时刻的潮流计算。4 仿真实例4.1 仿真条件 这里以国内某条地铁线路为例进行仿真。 对线路进行动态仿真需要获取每时刻的

7、列车信息以得到牵 引网的等效电路拓扑， 因而列车的运行数据要作为程序输入的一 部分。另外，列车运行时间设置比如发车间隔，上下车首车发车 时间等决定了线路上运行列车的密度以及列车位置的动态相对 关系，也是重要的运行参数。这里将发车间隔设置为300s，上下行首列车发车时间都设置为 0s，那么在首列车运行完全线后， 线路上的电气量将以 300s 为周期变化。由于列车运行全线的时 间接近3000s，为进行全线列车仿真，这里将仿真时间设置为 3000-3300s 。电气参数方面， 该线路为 750V 电压等级， 空载电压为 836V。 接触网单位长度电阻0.007 Q/km,钢轨单位长度电阻0.009

8、Q/km,对地漏电导取 0.005S/km。4.2 仿真结果与分析下面为 300s 发车间隔下 3000-3300s 的仿真结果。 以下是该线路上某变电所的电压、电流以及功率变化情况。图 4. 某站 3000-3300s 电压变化曲线图 5. 某站 3000-3300s 电流变化曲线图 6. 某站 3000-3300s 功率变化曲线 下表是线路上各变电所在 3000-3300s 内输出的电能 :表 1. 3000-3300s 各变电所电能输出变电站 1 2 3 4 5 6 7电能 /kwh 38.01 41.14 43.79 42.01 43.71 45.80 46.41变电站 8 9 10 11 12 13电能 /kwh 44.86 44.44