电力系统潮流计算实验报告

1、1 手算过程已知：节点1：PQ节点， s(1)= -0.5000-j0.3500 节点2：PV节点， p(2)=0.4000 v(2)=1.0500节点3：平衡节点，U(3)=1.00000.0000网络的连接图： 0.0500+j0.2000 1 0.0500+j0.2000231） 计算节点导纳矩阵由 ; ;导纳矩阵中的各元素：; ; ; ;形成导纳矩阵：2） 计算各PQ、PV节点功率的不平衡量，及PV节点电压的不平衡量：取:节点3是平衡节点，保持为定值。； ； ； ； ；于是： 3）计算雅可比矩阵中各元素 雅可比矩阵的各个元素分别为：又：； ； ； ； ； ； ； ； ；得到K=0时的雅

2、可比矩阵：4）建立修正方程组：解得：因为 ； ；所以 ；5）运用各节点电压的新值进行下一次迭代：即取: 节点3时平衡节点，保持为定值。6）计算各PQ、PV节点功率的不平衡量，及PV节点电压的不平衡量：；;于是： 7）计算雅可比矩阵中各元素：雅可比矩阵的各个元素分别为：又：； ； ； ； ； ； ； ；得到K=1时的雅可比矩阵：8）建立修正方程组如下：解得：因为 ； ；所以 ；9）运用各节点电压的新值进行下一次迭代：2. 程序流程图结束打印结果K=K+1计算平衡节点功率，计算各线路中功率计算新的电压初始值解修正方程式，求ei,fi计算雅可比矩阵各元素检验收敛Pimax，Qimax？计算Pi Qi

3、 UiK=0给定电压初始值形成导纳矩阵启动3. 源程序见电子版程序：余越.c4 输入数据3,2,0,1,1,0.000101,1,-0.5000,-0.3500,1.0000,0.00002,2,0.4000,1.0500,1.0000,0.00003,3,1.0000,0.00001,1,2,0.0500,0.20002,1,3,0.0500,0.20005 输出数据 潮流上机实习 华北电力大学 电自031 余越 7 * 原始数据 * = 节点数: 3 支路数: 2 对地支路数: 0 PQ节点数: 1 PV节点数: 1 精度:0. - PQ节点 节点 1 P1=-0. Q1=-0. PV节点

4、 节点 2 P2=0. V2=1. 平衡节点 节点 3 e3=1. f3=0. - 支路 1 相关节点: 1, 2 R=0. X=0. 支路 2 相关节点: 1, 3 R=0. X=0. = * 计算结果 * 节点导纳矩阵为: 2.35294+j-9.41177 -1.17647+j 4.70588 -1.17647+j 4.70588 -1.17647+j 4.70588 1.17647+j-4.70588 0.00000+j 0.00000 -1.17647+j 4.70588 0.00000+j 0.00000 1.17647+j-4.70588 =-输出雅可比矩阵和迭代次数-迭代的次数

5、为 1 9.41177 2.35294 -4.70588 -1.17647 -0.50000 -2.35294 9.41177 1.17647 -4.70588 -0.35000 -4.70588 -1.17647 4.70588 1.17647 0.40000 0.00000 0.00000 0.00000 2.00000 0.10250 =-输出 df,de-节点为 1 df=-0.01578 de=-0.02188节点为 2 df= 0.05094 de= 0.05125-输出迭代过程中的电压值- 节点为 1 f=-0.01578 e= 0.97812 节点为 2 f= 0.05094

6、e= 1.05125=-输出雅可比矩阵和迭代次数-迭代的次数为 2 9.51875 1.95000 -4.58437 -1.22500 -0.00541 -2.95000 8.81875 1.22500 -4.58437 -0.01555 -5.00699 -0.99706 4.74136 1.39706 -0.00697 0.00000 0.00000 0.10188 2.10250 -0.00522 =-输出 df,de-节点为 1 df=-0.00188 de=-0.00312节点为 2 df=-0.00343 de=-0.00232-输出迭代过程中的电压值- 节点为 1 f=-0.01

7、766 e= 0.97501 节点为 2 f= 0.04751 e= 1.04893=-输出雅可比矩阵和迭代次数-迭代的次数为 3 9.50304 1.95417 -4.56748 -1.23017 0.00001 -2.96650 8.76690 1.23017 -4.56748 -0.00007 -4.99205 -1.01046 4.72083 1.40412 -0.00003 0.00000 0.00000 0.09502 2.09787 -0.00002 =-输出 df,de-节点为 1 df=-0.00000 de=-0.00001节点为 2 df=-0.00001 de=-0.0

8、0001-输出迭代过程中的电压值- 节点为 1 f=-0.01766 e= 0.97499 节点为 2 f= 0.04750 e= 1.04892= 各节点电压为: U1=0. + j -0. U2=1. + j 0. U3=1. + j 0. = 平衡节点功率为: S3=0.+j 0. = 线路功率如下: 线路1-2的功率: -0. + j -0. 线路2-1的功率: 0. + j 0. 线路1-3的功率: -0. + j -0. 线路3-1的功率: 0. + j 0. = 线路上损耗的功率为: 线路1-2上的功率损耗: 0. + j 0. 线路1-3上的功率损耗: 0. + j 0. =

9、网络总损耗: 0. + j 0. *结束*6 程序说明1输入文件为 input.txt2输出文件为 output.txt3本程序采用直角坐标系,牛顿-拉夫逊法；4在输入文件中：1）节点数，总的支路数（包括对地支路）,对地支路数，PQ节点数，PV节点数，精度；用逗号隔开 2）节点编号，节点类型：1表示PQ节点：有功，无功，电压初值； 2表示PV节点：有功，电压，电压初值； 3表示平衡节点：电压值；用逗号隔开；3）支路编号，支路相关节点p1,p2;电阻和电抗;用逗号隔开；(有平行支路要连续输入)（对地支路要放在最后，若要输出对地支路的潮流则子程序还需适当的修改）；5在输出文件中输出了本程序计算的原

10、始数据、迭代过程及结果。7计算结果分析1由计算机所的结果：各节点电压为: U1=0. + j -0. U2=1. + j 0. U3=1. + j 0.平衡节点功率为: S3=0.+j 0.线路功率如下: 线路1-2的功率: -0. + j -0. 线路2-1的功率: 0. + j 0. 线路1-3的功率: -0. + j -0. 线路3-1的功率: 0. + j 0.网络总损耗: 0. + j 0.2结果分析1）各节点电压均在额定电压附近，平衡功率中无功功率为0.，即发出无功。2）所用线路有功都不超过1，没有过载，即转换功率比较大，运行较为安全。3）线路损耗中线路1-2较大，须提高电压等级，

11、以减少网损。3程序的优缺点1）优点：该程序采用牛顿拉夫逊算法，收敛较快。解方程用列主元素消去法，一般不会产生溢出。采用输入输出文件，可多次进行调用，且不因源程序多次修改数据而出错，减少数据输入的麻烦，方便看运行结果。本程序考虑到PQ、PV节点，计算的节点数，可根据电脑运行容量确定的计算方程的阶数来限定，一般情况下都能满足运算要求。运算精度可由用户自定义输入。通过多道题目的试验，本程序均可使用，即通用性良好。2）缺点：数据输入不美观，且输入格式要求严格。不能解决非标准变比变压器的潮流计算问题。没有建立良好的人机交互界面。8. 实验总结为期两周的电力系统潮流上机计算结束了，通过这两周的程序设计，手算题目等科目使我对运用计算机辅助计算电力系统潮流分布的原理和方法有了更深一步的了解，对电力系统潮流分布的概念有了更深刻的理解，并且通过程序来完成自己的理解。这套程序用VC+语言编程解决潮流计算问题，采用牛顿拉夫逊算法，计算速度快，收敛性比较好。生成的数据记录包括过程和结果，查看的时候比较方便，但是输入的格式比较死板，程序解决的问题还是比较简单，由于个人水平有限，所以没有