A 题 高铁牵引供电系统运行数据分析及等值建模
我国是世界上电气化铁路运营里程最长、服役电力机车型号最多、运营最
繁忙的国家。截至 2020 年底,我国铁路年消耗电量约 800 亿千瓦时,约占三峡
年总发电量的 80%。
图 1 为高速铁路牵引供电系统示意图,主要包括牵引变电所、牵引网和动车
组。牵引变电所将 220kV 三相电压转换为满足电力机车运行的 27.5kV 单相电压,
其核心元件是牵引变压器。牵引网由馈线、接触网、钢轨、回流线等组成,完成
对电力机车的送电任务。
高速动车组在运行过程中将消耗大量电能,且其功率随动车组运行工况不同
而变化,动车组运行也将带来区域电网电压、电流的不平衡和谐波问题,恶化电
能质量,给电网运行调控带来不良影响。由此,如何对含高铁牵引负荷的区域电
网电能质量进行评估,并制定节能降耗策略及对牵引负荷进行精准预测尤为重要。
对此,请结合上述背景知识及附录提供的相关实测数据,解决如下问题。
图 1 高速铁路牵引供电系统示意图
问题 1:电能质量评估
作为牵引供电系统的主要负荷,高速动车组负荷具有非线性、单相性、冲击
性等特点,导致了区域电网的谐波、负序等电能质量问题。请结合附录 1 提供的
牵引供电系统高压侧(220kV)实测三相电压、三相电流数据(采样间隔 50 微
秒),分析计算如下问题(提供可运行的源程序):
1) 高速动车组的主要运行工况可分为空载、牵引、制动,在三种运行工况下,
各取 0.2 秒时长的数据:【1,4000】(空载),【1600000,1604000】(牵引),
【4000000,4004000】(制动)
a. 计算各运行工况下电压和电流的正序、负序、零序分量,并绘制序分量图,
220kV外部电网外部电网动车组27.5kV牵引网牵引变电所能量传递谐波负序钢轨牵引变压器电能质量污染牵引能量消耗再生能量反馈回流线馈线
分析本段内电压和电流不平衡程度的差异。
b.绘制各运行工况电流频谱,计算幅值最大的 5 个谐波分量幅值及频率。
c. 请结合附录 1 中的国家电能质量评估标准,计算各运行工况下的电压、电
流不平衡度和电流总谐波畸变率,据此评价动车组运行工况对电压、电流
不平衡程度和电流谐波的影响。
2) 取全部数据为分析对象:
a. 高速动车组运行工况多变,可能带来电流不平衡程度和谐波含有率的变化,
试绘制电流正序、负序、零序分量幅值和 3 次谐波含有率随时间变化的曲
线,分析不平衡最严重时段的电流不平衡度和 3 次谐波含有率。
b.计算列车瞬时功率,并绘制其随时间变化曲线,试找出列车牵引功率、制
动功率最大的时刻,分析牵引功率、制动功率与三相不平衡度和谐波的关
系。
问题 2:节能降耗方案
高速动车组广泛应用了再生制动技术,在制动工况中将动能转换为电能,并
无偿地回馈至电力系统。结合附录 1 提供的现场实测数据,分析计算如下问题
(提供可运行的源程序):
a. 根据第 1 题中瞬时功率计算结果,统计高速动车组消耗和回馈的电量。
b. 为了将再生制动能量用于高速动车组牵引,请设计置于车载储能的再生制动
能量利用方案(含储能系统的功率电量配置及电力电子接入与控制方案),并
结合附录 2 所给电价、储能装置参数,测算其在降低动车组能耗方面的收益。
c. 请设计置于牵引变电所的再生制动能量利用储能方案(含储能系统的功率电
量配置及电力电子接入与控制方案),并结合附录 2 所给电价、储能装置参
数,测算其在降低动车组能耗方面的收益。
从经济、技术、安全可靠、储能容量利用率等方面比较 b、c 两方案的优劣。
问题 3:动态牵引负荷精确预测
研究表明,高铁牵引变电所动态牵引负荷(有功功率)需求与动车组在区间
运行过程(是否停靠站)、编组(8 编组和 16 编组)、运行方向(上行、下行)、
车型(不同车型额定功率各异)等因素密切相关。具体表征为:1)具有相同的
运行过程、编组、运行方向、车型的动车组,其在区间运行过程中动态功率需求
趋势相同、电度消耗/返送相近;2)对于同一车型,16 编组动车组的额定功率为
8 编组动车组额定功率的 2 倍;3)停靠站前动车组需进行短时间减速制动,将
产生大量的再生制动能量。附录 3 给出了某高铁牵引变电所全天的有功功率实测
数据,采样间隔为 1 秒,共计 86400 个数据。假定牵引变电所供电区间内不存在
同时给两列及以上动车组供电的场景。请结合实测数据,回答如下问题(提供可
运行的源程序):
a. 提取各列动车组在区间运行过程中的动态有功功率需求序列(动态牵引负
荷),并形成相应的动态牵引负荷数据库。
b. 应用相关性等分析方法,结合附录 3 提供的动车组运行方向、车型等信息,
实现所提取动态牵引负荷的分类和辨识。
c. 针对归类后的动态牵引负荷,应用线性回归等方法建立各类动态牵引负荷模
型,进而构建相应的动态牵引负荷模型库。
d. 结合附录 4 提供的列车时刻表,实现高铁牵引变电所动态牵引负荷的模拟和
电量的精确预测。
问题 4:高铁牵引变电所-牵引负荷等值建模
为了深入研究高速动车组运行对电网的影响,有必要建立反映高速动车组运
行的数学模型。
精准的数学模型的仿真结果应该能够与附录 1 所给数据相吻合。
图 2 给出了高铁牵引变电所-牵引负荷的示意图,请构建反映高速动车组运
行电气特性的数学模型,分析模型的电压、电流与附录 1 的实测电压、电流之间
的误差(提供可运行的源程序);说明模型有效性及改进的前景。
图 2 高铁牵引变电所-牵引负荷示意图
atvatibtvctvbticti钢轨回流线馈线
附录 1 三相电压、三相电流实测数据及电能质量标准
1. 三相电压、三相电流实测数据文件 Data.mat
Data.mat 文件中:Va、Vb、Vc,Ia、Ib、Ic 分别表示 A、B、C 三相瞬时电压值
(单位:V),三相瞬时电流值(单位:A);fs=20kHz,表示采样频率。
2. 相关电能质量标准及电能质量计算方法
[1]中华人民共和国国家标准. GB/T 14549-1993 电能质量 公用电网谐波.
1993.
[2]中华人民共和国国家标准. GB/T 15543-2008 电能质量 三相电压允许不
平衡度. 2008.
第 h 次谐波电压含有率
计算公式:
式中: 为第 次谐波电压(方均根值); 为基波电压(方均根值)。
第 h 次谐波电流含有率
计算公式:
式中: 为第 次谐波电流(方均根值); 为基波电流(方均根值)。
谐波电压含量 计算公式:
谐波电流含量 计算公式:
电压总谐波畸变率
计算公式:
电流总谐波畸变率
计算公式:
hHRUhh1100%UHRU=UhUh1UhHRIhh1100%IHRI=IhIh1IHU2Hhh2=U=UHI2Hh2hI=IuTHDHu1100%UTHDUiTHD
负序电压不平衡度 、零序电压不平衡度 的计算公式:
式中:U1 为三相电压的正序分量方均根值,单位为伏(V);U2 为三相电压
的负序分量方均根值,单位为伏(V);U0 为三相电压的零序分量方均根值,单
位为伏(V)。
将式中 U1、U2、U0 换为 I1、I2、I0 则为相应的电流不平衡度 和 的表达
式。
附录 2 高速动车组用电电价及用于节能降耗储能技术经济参数
高速动车组用电电价取 0.63 元/ kWh。用于节能降耗储能技术经济参数如附
表 1。
附表 1 储能参数
名称
单价/万元
系统集成/MW
隔离变压器/MW
并网逆变器/MW
四象限变流器/MW
双向DC/DC/MW
超级电容/MW
5
3
20
20
15
60
附录 3 高铁牵引变电所典型日负荷及列车运行信息
1. 高铁牵引变电所典型日负荷详见 P_TPSN.mat 文件
P_TPSN.mat 文件中:P_TPSN 表示负荷有功功率,单位为 W。
2. 列车运行信息详见附表 2
附表 2 给出了每趟高速动车组的编号、运行方向及车型。由表可知,该线路
每天运行动车组 29 趟,运行动车组车型为 CRH380D(车型 1)、CRH2A(车型
2)、CRH3C(车型 3)或 CRH380B(车型 4)。需要注意的是,实际运行过程中,
由于功率传输损耗等因素,导致动车组在额定功率需求下,牵引变电所的实际功
Hi1100%ITHDIU2U02U210U01100%100%U=UU=UI2I0
率较额定功率高。
附表 2 列车运行信息表
编号
1
2
3
4
5
6
7
8
运行方向 下行 上行 上行 上行 下行 下行 下行 下行
车型
编号
3
9
1
10
1
11
3
12
4
13
1
14
1
15
1
16
运行方向 上行 上行 下行 上行 下行 下行 下行 上行
车型
编号
1
17
1
18
1
19
1
20
2
21
1
22
3
23
1
24
运行方向 下行 上行 上行 下行 上行 下行 下行 下行
车型
编号
3
25
3
26
1
27
1
28
1
29
运行方向 下行 下行 下行 上行 下行
车型
4
2
1
1
1
2
3
1
附录 4 列车时刻表
车次编号
时间
车型
编组
停站
上下行
7:17
8:04
8:27
8:43
8:48
9:02
9:27
9:43
10:01
10:21
10:38
10:58
11:13
11:32
11:53
12:01
12:16
3
1
3
1
1
3
1
3
1
1
1
3
1
1
1
1
3
8
8
8
8
8
16
8
8
8
8
8
8
16
16
16
16
8
否
否
否
否
否
否
否
否
否
否
否
否
否
否
否
否
否
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
12:30
12:48
13:12
13:36
13:58
14:13
14:23
14:32
14:46
14:55
15:35
15:53
15:58
16:10
16:24
16:44
16:59
17:22
17:55
18:11
18:19
18:31
19:05
19:23
19:30
19:46
20:14
20:30
20:41
21:28
21:40
22:32
7:38
7:46
8:55
9:19
1
3
3
3
1
1
1
1
1
1
3
1
3
1
1
1
3
3
3
3
1
3
1
3
1
3
1
3
1
1
3
3
1
1
1
2
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
16
16
8
8
8
8
8
16
16
16
8
8
8
8
16
8
16
8
8
8
16
16
8
16
16
8
8
8
16
8
8
8
16
8
16
8
否
否
否
否
否
否
否
否
否
否
否
否
否
否
否
否
否
否
否
否
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
否
上行
否
否
否
否
否
否
否
否
否
否
否
否
否
否
否
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
上行
下行
下行
下行
下行
9:26
10:02
10:37
10:48
10:59
11:36
11:51
12:18
12:28
12:41
13:02
13:13
13:33
13:54
14:12
14:25
14:40
14:55
15:05
15:22
15:35
15:49
16:03
16:26
16:50
17:04
17:11
17:21
17:33
17:59
18:09
18:21
18:54
19:05
19:17
19:24
2
2
4
4
4
1
4
4
2
1
3
4
1
4
4
2
1
1
1
4
2
2
4
4
1
1
2
4
4
1
2
2
1
2
1
4
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
8
8
16
8
8
8
8
16
8
8
16
8
16
8
16
16
16
8
16
8
8
8
16
16
16
16
8
8
8
8
8
8
8
8
16
8
否
否
是
是
是
否
是
是
否
否
是
是
否
是
是
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
否
下行
否
否
否
是
否
否
是
是
否
否
否
是
是
否
否
否
否
否
否
是
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行
下行