复杂电力系统潮流计算课程设计.docx
《复杂电力系统潮流计算课程设计.docx》由会员分享，可在线阅读，更多相关《复杂电力系统潮流计算课程设计.docx（13页珍藏版）》请在冰豆网上搜索。
复杂电力系统潮流计算课程设计
华侨大学厦门工学院
电力系统综合设计课程设计报告
题目：
复杂电力系统潮流计算专业、班级：
10级电气
（2）班学生姓名：
学号：
指导教师：
黄永杰
分数:
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，r1=0.21Ω/km,x1=0.4Ω/km,b1=2.85×10-6s/km,线路长度分别为l1=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'smethod）又称为牛顿-拉夫逊方法，它是牛顿在17
世纪提出的一种在实数域和复数域上近似求解方程的方法。
方法使用函数f（x）的泰勒级数的前面几项来寻找方程f（x）=0的根。
牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f（x）=0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。
线性网络的常用解法有节点电压法和回路法，前者须列写节点电流平衡方程，后者则须列写回路方程。
一般的，对于有n个独立节点的网络，可以列写n个节点方程
Y11U1Y12U2Y1nUnI1
Y21U1Y22U2Y2nUnI2
Yn1U1Yn2U2YnnUnIn
也可以用矩阵写成
或缩写为
YUI
对潮流计算的要求可以归纳为下面几点：
（1）算法的可靠性或收敛性
（2）计算速度和内存占用量
（3）计算的方便性和灵活性牛顿法，由于其在求解非线性潮流方程时采用的是逐次线性化的方法，为了进步提高算法的收敛性和计算速度，人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来，于是产生了二阶潮流算法。
后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点，提出了采用直角坐标的保留非线性快速潮流算法。
三、计算
解：
1、
（1）线路参数的标幺值：
自导纳：
互导纳：
Y34Y430.83j1.5809；
Y41Y141.245j2.3714；
Y13Y311.66j3.1619；
Y12Y21j5.3182
2.905j11.5031
j5.3182
1.66j3.1619
1.245j2.3714
YB
j5.3182
j4.6638
0
0
1.660j3.1619
0
2.49j4.7040
0.83j1.5809
1.2450j2.3714
0
0.83j1.5809
2.075j3.9092
（0）
3、初值：
ei1；
（0）i
4、
计算各节点功率的不平衡量
取
U
1（0）1j0；
U2（0）
1j0；
U
（0）
3
1.05j0；
U4（0）
1.05j0
jn
P
（0）
ei（0）Gi
（0）ej
Bijfj（0）ijj
fi（0）Gijfj（0）Bije（j0）
j1
Qi（0）fi（0）Gije（j0）Bijfj（0）ei（0）Gijfj（0）Bije（j0）
j1
经计算得：
P1（0）0.14525；P2（0）0；P3（0）0.08710
Q1（0）0.37494；Q2（0）0.65434；Q3（0）0.12330
又Pi（0）PiPi（0）;
Qi（0）QiQi（0）
P1（0）0.14525
；P2（0）0.50000
P3（0）0.11290
Q1（0）0.37494
;Q2（0）0.35434
Q3（0）0.12330
5、计算雅克比矩阵中各元素：
先计算各节点注入电流
0.1453j0.3749a1（10）jb1（10）相似地可得：
a2（20）0；a3（30）0；
b2（20）0.654；4b3（30）0.1174
∴计算雅克比矩阵各元素
H1（10）B11e1（0）G11f1（0）b1（10）11.503112.9050（0.3749）11.1282
N1（10）G11e1（0）B11f1（0）a1（10）2.9051（11.5031）0（0.1453）2.7597
J1（10）G11e1（0）B11f1（0）a1（10）2.9051（11.5031）0（0.1453）3.0503
L1（10）B11e1（0）G11f1（0）b1（10）11.503112.9050（0.3749）11.878
H2（20）B22e2（0）G22f2（0）b2（20）4.66381.000.65445.3182
N2（20）G22e2（20）B22f2（0）a2（20）01.0（4.6638）000
J2（20）G22e2（0）B22f2（20）a2（20）01.0000
L（202）B22e2（0）G22f2（0）b2（20）4.66381.000.65444.0094
H3（30）B33e3（0）G33f3（0）b3（30）4.7041.052.4900.11745.0566
N3（30）G33e3（30）B33f3（0）a3（30）2.491.05（4.704）00.0832.6975
R332f3（0）0；S332e3（0）2.1
H1（20）B12e1（0）G12f1（0）5.31821005.3182
H1（30）B13e1（0）G13f1（0）3.16191（1.66）03.1619
N1（20）G12e1（0）B12f1（0）015.318200
N1（30）G13e1（0）B13f1（0）1.6613.161901.66
J1（20）B12f1（0）G12e1（0）5.31820010
J1（30）B13f1（0）G13e1（0）3.16190（1.66）11.66
L1（20）G12f1（0）B12e1（0）005.318215.3182
L1（30）G13f1（0）B13e1（0）1.6603.161913.1619
H2（10）B21e2（0）G21f2（0）5.31821.0005.3182
H2（30）B23e2（0）G23f2（0）01.0000
N2（10）G21e2（0）B21f2（0）01.05.318200
N2（30）G23e2（0）B23f2（0）00000
J2（10）B21f2（0）G21e2（0）5.31820010
（0）
J23
B23f2（0）G23e2（0）0
G21f2（0）B21e2（0）05.318215.3182
H3（10）B31e3（0）G31f3（0）3.16191.05（1.660）03.3199
H3（20）B32e3（0）G32f3（0）0
N3（10）G31e3（0）B31f3（0）1.661.053.161901.743
N3（20）0
J3（10）B31f3（0）G31e3（0）3.16190（1.66）1.051.743
（0）
32
R310；
S31
0；
R320；
S32
0
列出k=0时的雅克比矩阵
11.1282
3.0503
2.7597
11.878
5.31820
J（0）
5.3182
0
5.3182
0
5.3182
0
3.320
1.743
0
0
0
0
0
3.1619
1.66
5.3182
1.66
3.1619
0
0
0
4.0094
0
0
0
5.0566
2.6975
0
0
2.1
6.逆矩阵为：
0.2304
0.0713
0.2304
0.0946
0.1675
0.1403
0.0837
0.1594
0.0837
0.2115
0
0.3062
0.2304
0.0713
0.4184
0.0946
0.1675
0.1403
0.1110
0.21115
0.1110
0.5299
0
0.4062
0.1801
0.0081
0.1801
0.0108
0.3077
0.2406
0
0
0
0
0
0.4762
（J（0））1
7.
迭代过程中各节点功率的不平衡量：
k=0时：
Pi（k）jQi（k）P1（0）jQ1（0）0.1453j0.3749
P2（0）jQ2（0）0.5j0.3544
P3（0）jQ3（0）0.1129j0.1233
8.求得各节点电压的新值后，就开始第二次迭代。
每次迭代所得示于表1~4
表1迭代过程中各节点功率的不平衡量
K
P1（k）jQ（1k）
P2（k）jQ2（k）
P3（k）jQ3（k）
0
0.14525-j0.37494
-0.50000+j0.35434
0.11290-j0.12330
1
0.02520-j0.01318
0.04603-j-0.06592
0.00668-j0.00545
2
-0.00016+j0.01465
-0.11269-j0.01217
0.00092-j0.00897
表2迭代过程雅克比矩阵各对角元
k
H（k）
H11
（k）L11
H（k）
H22
（k）
L22
H3（k3）
33
S（3k3）
0
11.1282
11.8780
5.3182
4.0094
5.0566
2.6975
1
10.5985
10.9019
5.0415
4.4678
4.6893
2.6571
2
10.8023
10.7712
5.0235
4.3689
4.7006
2.4654
表3迭代过程中各节点电压的修正量
k
e1（k）jf1（h）
（k）（h）e2jf2
e3（k）jf3（h）
0
-0.05224-j0.0523
0.01905-j0.1463
-0.05872+j0.00131
1
-0.06030+j0.0272
-0.0895-j0.0740
-0.01170+j0.06660
2
-0.01120-j0.0115
-0.0218-j0.0354
0.00440+j0.00520
表4迭代过程中各节点电压
k
（k）（k）e1jf1
e（2k）jf2（k）
e（3k）jf3（k）
0
1.00000+j0.00000
1.00000+j0.00000
1.05000+j0.00000
1
0.94776-j0.05230
1.01905-j0.14630
0.99128+j0.00131
2
0.88746-j0.02510
0.92955-j0.22030
0.97958+j0.06791
9、计算平衡节点功率S~4和线路功率S~ij
jn**
~
平衡节点功率：
S4U4Y4jUjj1
（1.05j0）[（1.245j2.3714）（0.94336j0.0223）0
（0.83j1.5809）（0.98848j0.02301）（2.075j3.909）（1.050j）]
0.210420.3108j
线路功率：
S~12U1[U1y10（U1U2）y12]0.5j0.307
S~21U2[U2y20（U2U1）y21]0.5j0.246
S310.256j0.195；S310.243j0.197；
S~430.057j0.057；S~340.056j0.002
网络总损耗：
in
S~S~i0.21042j0.31080.5j0.30.20.15j0.11.06042j0.7108
i1
四、结论
这次的电力系统分析课程设计让我对平时所学的专业知识有了更深刻更具体的了解，明白了理论知识必须与实践相结合才能更好的发挥作用。
在不停的翻书上网查资料的过程中，我积累了大量的潮流计算以及电力系统的知识，全面透彻的了解了相关知识。
本设计采用直角坐标形式的牛顿—拉夫逊法作常规潮流计算。
P-Q分解法利用了电力系统的一些特有的运行特点，对牛顿—拉夫逊法作了简化，可提高计算速度，但较难理解，牛顿—拉夫逊法的雅克比矩阵在每次迭代过程中都有变化，需要重新形成和求解，这占用了计算的大部分时间，成为牛顿—拉夫逊计算速度不能提高的主要因素，但收敛性好，物理概念也较为清晰
通过本次设计，我最大的感慨是学习这些知识的时候没有投入全部的精力，以致课设刚开始时花费了大量时间补充基础知识。
计算机的操作还需提升，英语也要加强。
从这些事我下决心以后上课时必须认真听讲，不能虚度光阴。
五、致谢本次课程设计是在我的老师和的亲切关怀和悉心指导下完成的。
她们严肃的科学态度，严谨的治学精神，精益求精的工作作风，深深地感染和激励着我。
我坚信这次的课设为我的毕业论文和未来打下了很好的基础，使我能更从容的面对以后生活和学习上的种种挑战。
同学与老师的帮助也是必不可少的，在我遇到困难的时候，得到了来自同学的协助与老师的指导，使我克服了困难，冲破了阻碍，如期完成了任务，在此特别感谢他们。
我会从此走的更远。
参考文献
[1]胡学林编著.可编程控制器原理及应用（第二版）．北京：
电子工业出版社，2012.7.
[2]陈珩等.电力系统稳态分析第三版北京：
中国电力出版社，2007.6
[3]周全仁等.电网计算与程序设计[M].长沙：
湖南科学技术出版社,1983.
[4]张伯明,陈寿孙.高等电力网络分析[M].北京：
清华大学出版社,1996.
[5]陈珩.电力系统稳态分析.水利电力出版社,1994
[6]杨少兵,骆平.电力系统分析的教学软件开发，电力系统潮流分析.华北电力技术,2000（10）
升级会员 