华侨大学工学院
电力系统综合设计课程设计报告
题 目：复杂电力系统潮流计算
专业、班级：10级电气⑵班 学生： _______________________ 学 号： ________________________ 指导教师： _________ 黄永杰 _________ 分 数： _____________________________
2013年_6—月26 日
摘要 ..... 一、任务书
、基础资料 二、计算
3.1 节点导纳矩阵 ............................ 5 3.2 设定所求变量的初值 ......................... 6 3.3 计算修正方程 ............................ 7 3.4 形成雅可比矩阵 .......................... 9 3.5 求解修正方程 .......................... 10 3.6 进行修正和迭代 .......................... 10 3.7 迭代精度的确认 ........................... 11 3.8
各节点电压计算功率分布 ...................... 11
四、 结论 ................................................................. 五、 ...................................... 致谢
13
六、 参考文献 ............................................................. 13 14
本次的课程设计主要针对复杂电力系统进行潮流计算。对电力网络的各种设 计方案及各种运行方式进行潮流计算， 可以得到各种电网各节点的电压，并求得 网络的潮流。采用牛顿-拉夫逊算法，牛顿一拉夫逊法是数学上解非线性方程式 的有效方法，有较好的收敛性。将牛顿法用于潮流计是以导纳矩阵为基础， 由于 利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧，使牛顿法在收敛 性、占用存、计算速度等方面都达到了一定的要求。
关键词：潮流分布 迭代 牛顿-拉夫逊算法
、任务书
题目二：如图二所示电力系统接线图，系统额定电压为 110KV各元件参数为
LGJ-120，r仁0.21 Q /km,x1=0.4 Q /km,b1=2.85 x 10-6s/km,线路长度分别为 I仁 150km,l2=100km,l3=75km.变压器容量为 63000KVA额定电压为 110/38.5KV , 短路电
压百分数为10.5，变压器的实际变比为1.1282，电容器导纳为j0.05。
取 SB=100MVA,UB=UN.
取节点4为平衡节点，节点3为PV节点，节点1,2均为PQ节点。
1.试用直角坐标表示的牛顿一拉夫逊计算系统中的潮流分布。 0.001)
(迭代精度为
、基础资料
牛顿迭代法(Newton's method )又称为牛顿-拉夫逊方法，它是牛顿在
17
世纪提出的一种在实数域和复数域上近似求解方程的方法。方法使用函数 f(x) 的泰勒级数的前面几项来寻找方程 f(x) = 0 的根。牛顿迭代法是求方程根的重 要方法之一，其最大优点是在方程 f(x) = 0 的单根附近具有平方收敛，而且该 法还可以用来求方程的重根、复根。
线性网络的常用解法有节点电压法和回路法，前者须列写节点电流平衡方
程，后者则须列写回路方程。
一般的，对于有n个节点的网络，可以列写n个节点方程
丫1厲 丫必 L YnUn &
WnU& & 冷出Y22U& L
Yn1® Yn2l&2 L 绻说 &
丫
也可以用矩阵写成
Y1n 丫.1 Yn2 L Ynn 观 或缩写为
1&
2
Y2i
YU I
对潮流计算的要求可以归纳为下面几点: (1) 算法的可靠性或收敛性 (2) 计算速度和存占用量
(3) 计算的方便性和灵活性 牛顿法，由于其在求解非线性潮流方程时采用的是逐次线性化的方法， 为了进 步提高算法的收敛性和计算速度，人们考虑采用将泰勒级数的高阶项或非线性项 也考虑进来，于是产生了二阶潮流算法。后来又提出了根据直角坐标形式的潮流 方程是一个二次代数方程的特点，提出了采用直角坐标的保留非线性快速潮流算 法。
三、计算
解：1、( 1)线路参数的标幺值：
Ze1 Ze2 Ze3
0.2603 j 0.4959 0.1736 j 0.3306 0.1302 j 0.2479
yel / 2 j 0 .0259 ye2 / 2
j 0.0172
ye3 /2 j 0.0129
(2) 变压器参数的标幺值：
Zn Yn2
kZT 1 k kZ
j0.188 j 0.6044
Yni
k 1 kZT
j0.6818
2、各串联支路导纳:
y12
y31
j5.3182 ；
1.66 j3.1619 ；
y41 1.245 j2.3714 ；
1
y43
0.2603 j 0.4959
0.00 j 1.5809
自导纳:
yn
y
2.905
j11.503
;
33
2.49 j4.7 ；
Y22
；
j 4.6638
互导纳:
y 2.075 j 3.9092
0.83 j1.5809 ； 1.245
j2.3714
；
Y
Y34 Y43 Y41 Y14 Y13 Y31
Y42 Y24 0 ；
23 Y 32
1.66 j3.1619 ；
Y12 Y21 j5.3182
2.905 j 11.5031 i
j5.3182
Y
B
1.660 j3.1619
1.2450 j 2.3714
j n
P
(0)
e(0)
Gije
(0)
j 1 j n
Qi
(0)
fi(0)
Gje(
0)
j 1
经计算得：
R
(0)
0.14525 ； P2(0)
Qi
(0)
0.37494 ； Q2
0)
又 P(0)
P R
(0)
；
P(0) 0.14525
P(0) 0.11290 Q10) 0.37494 Q30)
0.12330
4、计算各节点功率的不平衡量
(0)
3、初值：e（0）
i
i
（0）
Gij fj
（0）
Bij e
(0)
j5.3182 1.66
j 3.1619 1.245 j2.3714
j 4.6638
0
0
0
2.49 j 4.7040 0.83 j1.5809 0
0.83 j 1.5809
2.075 j3.9092
Bj f j
(0)
Bj f j
(0)
)
)
Go fj
(0
Bij e
(0
0 ； P3
(0)
0.08710
0.65434 ； Q3
0)
0.12330
Qi(0)
Qi Qi(0)
P2(0) 0.50000
Q20)
0.35434
取 U,0) i jo ； U
0)
20)
i jo ；
u3
0)
1.05 j0 ； u41.05 j0
5、计算雅克比矩阵中各元素:
先计算各节点注入电流
(0)
I
R(0) jQi(0) * (0)
i
ai( jbi(
0)
0)
Ui
1
(0) 1
P
(0)
jQi(0)
(0)
0.1453 j0.3749 1.0 j0
(0)
(
0.1453 j0.3749 a1 jbi0)
相似地可得:
•••计算雅克比矩阵各元素
BueT ；11 1(0)
10)
GneT ) ) B f(°a(°
11 11 11
10) (0)
!0) 11 1(0) 110) 11
(0) G f(0) (0) 11 11 h
0.6544
； bT
11.5031 1 2.905 1 2.905
0.1174
2.905 0 11.5031) 0 (11.5031)
(0.3749) (0.1453) 0 ( 0.1453)
11.1282 2.7597
3.0503
：
HNJ
Gf
；
GneBfa
Bne
22
(0) J22
30)
H2
Be2
G22e22)
0)
G f
(°)
22 2
11.5032.905 0.3749) 11.878
1
4.6638 1.0 0 0.6544 5.3182 0 1.0
(4.6638) 0 0
0
4.0094 0.1174 5.0566 0.083 2.6975 0 0
? N2)2
H3
(0) (0) 22 f 22
0
G22e2
0)
B
a22
0 1.0
G22
(0) b(0) 2 S
4.6638 1.0 0 0.6544 4.704 1.05 2.49 0 2.49 1.05 (
S
3N)
B33事
G33e3?
f
(0) (0) 3 5
B33 f3
(0)
4.704) 0
2 f3(0)
B12eT B13e
(0)
33
2.1
5.3182
(1.66) 0
3.1619
G12 f1G13 f1
(0)
5.3182 3.1619
(0)
Nj G12
e1
0)
NJ
G^e⑼
(0) J12 B,⑼
(0) (0) J13 B f
13 M
(0) <0) 12 G2 f1
(0) 13
G”
(0) H 21
B21&6
(0) H 23 B23e20) N(0) G21
e2
0)
21
N(0) G23e2
0)
23 (0) 21 B f
(0)
2112
(0) 23 B23f
2(0)
(0) L21 G2J2⑼
|
(0) G f
(°)
L23 23
1 2
(0) H 31 B31e30)
H(o) H
32 B32e30) N 3
Gaie3
0)
0N3)
2)(0)
J
31 B3/3⑼ (0) J32
G/y B12 f1
(0)
0 1 5.3182
0 0
B13 f1 ⑼
1.66 1 3.1619 0 1.66 G2
辭
5.3182
0 0 1 0
G1
3e1(0)
3.1619 0 (1 .66) 1 1.66 Bg® 0 0 5.3182 1
5.3182
B13e
(0)
1.66 0
3.1619 1
3.1619 G21 f2
(0)
5.3182 1.0 0 0
5.3182
0)
G23 f2
0 1.0 0 0 0
B21 f V 0
1.0 5.3182
0 0
B23 f
2(0>
0 0 0
0 0
G21
e2
0)
5.3182 0 0
1 0
G23e20：
)0
B21e20)
0
5.3182 1
5.3182
B23e2
0)
0
G31 f. 0)
3.1619 1.05
( 1.660) 0
3.3199
G32 f3(
：0)
0
B31 f J
1.66 1.05 3.1619 0
1.743
G31&
0)
3.1619 0 ( 1.66) 1.05 1.743
31e3
0)
1.66 0 3.1619 1.05 3.320
BR
31
0
； 為
0
；
R32 0 ； S32 0
•••列出k=0时的雅克比矩阵：
11.1282 2.7597 5.3182 3.0503 11.878 0 (0)5.3182 j⑼
0 5.3182 0
5.3182 0 3.320 1.743 0
0
0
0
6. 逆矩阵为:
0.2304
0.0713 0.2304 0.0837 0.1594 0.0837
(J(0)) 1
0.2304
0.0713 0.4184
0.1110 0.21115 0.1110 0.1801 0.0081 0.1801
0 0
0
0 5.3182 0 4.0094 0 0
0.0946 0.2115 0.0946 0.5299 0.0108 0
3.1619 1.66 0 0 5.0566 0
0.1675 0 0.1675 0 0.3077 0
1.66 3.1619 0 0 2.6975 2.1
0.1403 0.3062 0.1403 0.4062 0.2406 0.4762
7.
R(0)
0.3749 0.1453 0.50 0.3544 0.1129 0.1233
i
(0) (0)
i
(0) f1
e
P2
QQ
11
（0）
(0)
f (0) f
2
230)
e2e(
0)
P3
(1)
（0）
(0) f 3
0)
0)
0.05230 0.05224
0.14630
0.01905 0.00131 0.05872
f1⑴ e⑴ f2
(1)
0.05230 0
0.94776 0.14631.01905 0.00131 0.99128
e2f3e(
1)
(1)
1)
迭代过程中各节点功率的不平衡量:
k=0 时：
R
(k)j Qi
(k)
R j Q1
(0)
(0)
0.1453
j0.3749
P
2
(0)
0
j Q2
0)
0.5 j0.3544
j0.1233
P3
j Q3 0.1129
0)
8.求得各节点电压的新值后，就开始第二次迭代。每次迭代所得示于表 1~4
表1迭代过程中各节点功率的不平衡量 K
P1(k) j Q1k) P2(k) j Q2(k) P3(k) j Q3k) 0 0.14525-j0.37494 -0.50000+j0.35434 0.11290-j0.12330 1
0.02520-j0.01318 -0.00016+j0.01465 0.04603-j-0.06592 -0.11269-j0.01217 0.00668-j0.00545 0.00092-j0.007 2
表2迭代过程雅克比矩阵各对角元 k 0 1
H仆) H11 I (k) L11 H() H 22 kI (k) L22 H33 )33 11.1282 10.5985 10.8023 11.8780 10.9019 10.7712 5.3182 5.0415 5.0235 4.0094 4.4678 4.36 5.0566 4.63 4.7006 2.6975 2.6571 2.4654 2
表3迭代过程中各节点电压的修正量 k 0 1
e(k) 1 J f1 . f (h) e(k) 2 j f. r ( h) 2 e3k) j f3(h) -0.05224-j0.0523 -0.06030+j0.0272 -0.01120-j0.0115 0.01905-j0.1463 -0.05-j0.0740 -0.0218-j0.0354 -0.05872+j0.00131 -0.01170+j0.06660 0.00440 +j0.00520 2
表4迭代过程中各节点电压 k 0 1
e1 (k) jf* (k) 1 e2k) jf2(k) e3k) f(k) 1.00000+j0.00000 0.94776-j0.05230 0.88746-j0.02510 1.00000+j0.00000 1.01905-j0.14630 0.92955-j0.22030 1.05000+j0.00000 0.99128+j0.00131 0.97958+j0.06791 2
9、计算平衡节点功率S4和线路功率Sj
j n *
*
平衡节点功率：
S4 U4 丫4jUj
j 1
(1.05 j0)[( 1.245 j2.3714)(0.94336 (0.83 j1.5809)(0.98848 j0.0223) 0
j3.909)
j0.02301) (2.075
线路功率：
S
12
S21
同理：
S41
S31
S43
网络总损耗:
i n
S
Si i 1
(1.05 0j)]
0.21042 0.3108j
* * * * *
U1[U1 y10 (U1 U 2) y12]
0.5 j0.307 *
* *
* *
U2[U2
y
20
(U 2
0.5 j 0.246
0.271
j0.101 ；
S14
0.257 j0.110 ；
0.256 j0.195 S31
0.243 j 0.197 ；；
0.057 j
瓦4
0.05j0.002
0.057
；
6
0.21042 j 0.3108 0.5 j0.3 0.2 0.15 j0.1 1.06042 j 0.7108
四、结论
这次的电力系统分析课程设计让我对平时所学的专业知识有了更深刻更具
体的了解，明白了理论知识必须与实践相结合才能更好的发挥作用。 在不停的翻 书上网查资料的过程中，我积累了大量的潮流计算以及电力系统的知识， 全面透 彻的了解了相关知识。
本设计采用直角坐标形式的牛顿一拉夫逊法作常规潮流计算。
P-Q分解法利
用了电力系统的一些特有的运行特点， 对牛顿一拉夫逊法作了简化，可提高计算 速度，但较难理解，牛顿
一拉夫逊法的雅克比矩阵在每次迭代过程中都有变化， 需要重新形成和求解，这占用了计算的大部分时间，成为牛顿一拉夫逊计算速度 不能提高的主要因素，但收敛性好，物理概念也较为清晰
通过本次设计，我最大的感慨是学习这些知识的时候没有投入全部的精力， 以致课设刚开始时花费了大量时间补充基础知识。计算机的操作还需提升，英语 也要加强。从这些事我下决心以后上课时必须认真听讲，不能虚度光阴。
五、致谢
本次课程设计是在我的老师和的亲切关怀和悉心指导下完成的。她们严肃 的科学态度，严谨的治学精神，精益求精的工作作风，深深地感染和激励着我。 我坚信这次的课设为我的毕业论文和未来打下了很好的基础，
使我能更从容的面
对以后生活和学习上的种种挑战。 同学与老师的帮助也是必不可少的， 在我遇到 困难的时候，得到了来自同学的协助与老师的指导， 使我克服了困难，冲破了阻 碍，如期完成了任务，在此特别感谢他们。我会从此走的更远。
参考文献
[1] 胡学林编著•可编程控制器原理及应用(第二版)•北京：电子工业， 2012.7.
[2] 珩等.电力系统稳态分析第三版 北京：中国电力，2007.6 [3] 周全仁等.电网计算与程序设计[M].:科学技术，1983. [4] 伯明,寿.高等电力网络分析[M].北京：清华大学,1996. [5] 珩.电力系统稳态分析.水利电力,1994
⑹ 少兵,骆平.电力系统分析的教学软件开发，电力系统潮流分析.华北电
力技术,2000(10)