暂态稳定分析仿真
考虑如图 2 所示的 9 节点电力系统。该系统有三台发电机、3 个负荷以及 9
条支路。支路数据和发电机参数分别列于表 1 和表 2,系统频率为 60Hz。干扰是
在零秒线路 5—7 靠近母线 7 处发生三相短路,故障在 5 个周波(约 0.08333)
由断开线路 5—7 而被消除。
2
3
2
7
8
9
3
5
6
4
1
1
图 2
表 1 支路数据
首端母线名 末端母线名
电阻(标幺
电抗(标幺
容纳之半(标
变压器非标
值)
0.010
0.017
0.032
0.039
0.0085
0.0119
0
0
0
值)
0.085
0.092
0.161
0.170
0.072
0.1008
0.0576
0.0625
0.0586
幺值)
0.088
0.079
0.153
0.179
0.0745
0.1045
准变比
1.0
1.0
1.0
5
6
7
9
8
9
4
7
9
4
4
5
6
7
8
1
2
3
发电 母线
JT
aR
dX
'dX
qX
'qX
0'dT
0'qT
D
表 2 发电机数据
机
名
1
2
3
1
2
3
47.28
12.80
6.02
0
0
0
0.1460
0.0608
0.0969
0.0969
0.8958
0.1198
0.8645
0.1969
1.3125
0.1813
1.2578
0.2500
8.96
6.00
8.59
0.535
0.600
0
0
0
2.1 初值计算结果
利用 Matlab 进行潮流计算得到系统正常运行情况下的参数,运行结果如下
所示:
整理可得到表 3:
表 3 正常运行情况下的系统潮流
母线名
电压
发电机
负荷
幅值
相角
有功功率 无功功率 有功功率 无功功率
0.7164
1.6300
0.8500
0.2705
0.0665
-0.1086
1
2
3
4
5
6
7
1.0400
1.0250
1.0250
1.0258
0.9956
1.0127
1.0258
0
9.2800
4.6647
-2.2168
-3.9888
-3.6874
3.7197
1.2500
0.9000
0.5000
0.3000
8
9
1.0159
1.0324
0.7275
1.9667
1.0000
0.3500
各发电机的暂态电动势,功角和输入机械功率初值计算结果如下:
计凸极效应:
不计凸极效应(令
X
q
X ):
'
d
整理得到表 4:
表 4 发电机的 '
q
E 、 及P
(0)
(0)
m
发电机
不计凸极效应
计凸极效应
'
qE
1.0566
1.0502
1.0170
(0)
2.2716
19.7316
13.1664
'
qE
1.0564
0.7882
0.7679
(0)
3.5857
61.0984
54.1366
1
2
3
mP
(0)
0.7164
1.6300
0.8500
各节点的等值导纳计算结果如下:
即,负荷(节点 5):1.2610—j0.5044
负荷(节点 6):0.8776—j0.2925
负荷(节点 8):0.9690—j0.3391
2.2 故障系统和故障后系统描述
故障期间的电网络相当于在 7 号母线处并联一条阻抗为零的接地支路,这时
只要将正常情况下导纳矩阵Y 中的对角元素 77Y 改为无穷大(实际计算中取 2010 ),
即可得到故障期间系统的节点导纳矩阵。
故障后的电网络是切除线路 5—7 后的情况,由于线路 5—7 对导纳的贡献为:
Y
l
(5 7)
5
7
1
r
jx
1
r
jx
5
7
jb
r
1
jx
jb
1
r
jx
式中: 0.032,
r
x
0.161,
b
0.153
Y
,因此故障后导纳矩阵为:
P
Y Y
(5 7)
l
。
2.3 相对摇摆角计算结果
根据系统的初值,利用 Matlab 求解网络方程和微分方程得到相对摇摆角计
算结果如图 3 和图 4 所示。
160
140
120
1
2
/
度
角
摆
摇
对
相
100
80
60
40
0
90
80
70
60
50
40
30
20
10
1
2
/
度
角
摆
摇
对
相
0.5
1
时间/s
1.5
2
图 3 计凸极效应
0
0
0.5
1
时间/s
1.5
2
图 4 不计凸极效应
由图 3 和图 4 可以看出,无论计凸极效应还是不计凸极效应,系统都是暂态
稳定的。在计凸极效应的情况下,最大相对摇摆角为 149.8768 度。在不计凸极
效应的情况下,最大相对摇摆角为 85.3312 度,而第二摆的角度 85.1118 度比第
一摆的角度小。
最后,对不计凸极效应和计凸极效应情况下的临界切除时间进行仿真模拟。
得到前者所对应的临界切除时间在 0.163s—0.164s 之间,后者所对应的临界切
除时间在 0.086s—0.087s 之间。它们对应的摇摆曲线分别如图 5 和图 6。
1
2
/
度
角
摆
摇
对
相
160
140
120
100
80
60
40
20
0
-20
-40
0
160
140
120
1
2
/
度
角
摆
摇
对
相
100
80
60
40
0
1200
1000
800
1
2
/
度
角
摆
摇
对
相
600
400
200
0
0
0.5
1
时间/s
1.5
2
0.5
1
时间/s
1.5
2
图 5 不计凸极效应相对摇摆角与时间的关系曲线
(左切除时间 0.163s,右切除时间 0.164s)
900
800
700
600
1
2
/
度
角
摆
摇
对
相
500
400
300
200
100
0
0
0.5
1
时间/s
1.5
2
0.5
1
时间/s
1.5
2
图 6 计凸极效应相对摇摆角与时间的关系曲线
(左切除时间 0.086s,右切除时间 0.087s)