课程设计任务书1PQ法潮流计算程序设计电气1201乌通.docx

1、课程设计任务书1PQ法潮流计算程序设计电气1201乌通稳态分析课设目录1.任务书 22.模型简介及等值电路 33.修正方程的建立 44程序流程图. 55. MATLAB程序编写. 76.结果分析 97.设计总结 108.参考文献 11 电力系统分析课程设计任务书题 目利用PQ法实现电力系统潮流计算程序设计学生姓名 乌通学号201205080810专业班级电气1201设计内容与要求1. 设计要求：掌握MATLAB语言编程方法；理解和掌握运用计算机进行潮流计算的基本算法原理；针对某一具体电网，进行潮流计算程序设计。其目的在于加深学生对电力系统稳态分析中课程中基本概念和计算方法的理解，培养学生运用所

2、学知识分析和解决问题的能力。2. 内容1）学习并掌握MATLAB语言。2）掌握变压器非标准变比概念及非标准变比变压器的等值电路。掌握节点导纳矩阵的概念及导纳矩阵的形成和修改方法。3）掌握电力系统功率方程、变量和节点分类。4）掌握利用PQ法进行潮流计算的方法和步骤。5）选择一个某一具体电网，编制程序流程框图，用MATLAB语言进行PQ法潮流计算程序设计。6）利用MATLAB语言编制潮流计算的程序，并上机调试程序，对计算结果进行验证分析。7）整理课程设计论文。起止时间2015 年 5 月 25日 至 2015 年 5 月 29日指导教师签名2015年 5 月 20 日 系（教研室）主任签名 201

3、5年 5 月 21日学生签名2015年 5 月 22 日 二 模型简介图节点导纳之路等值电路闭式电力网结线图如图所示，发电厂G2为一基载厂。承担固定负荷，包括高压母线上的负荷其运算功率为117-j22MVA，变电所H的运算负荷（包括高压母线负荷）为440+j136MVA。发电厂G1承担除G2厂供给功率之外的系统所需负荷及网络损耗。线路L1，L2，L3的阻抗分别为3.17+j20.7，5.13+j27.2及7.33+j48.0。已知发电厂G1的高压母线电压U1为238kV，试计算潮流分布。三 修正方程式的建立PQ分解法潮流计算时的修正方程是计电力系统的特点后对牛拉法修正方程的简化。 由牛拉法修正

4、方程式，可得出P-Q分解法潮流计算的修正方程式：或者简写为： 对修正方程式的第一个化简是：计电力网络的各元件的阻抗一般远大于电阻，以致各节点电压相位角的改变主要影响各元件中的有功功率潮流从而各节点的注入功率流；各节点电压大小的改变主要影响各元件中的无功功率潮流，从而各节点的注入无功功率；可将修正方程式简化为：电力网络上的特点，可以由电力系统分析第三章中对纵，横向附加电势的讨论得到证实。对修正方程式的第二个简化基于对状态变量i的约束条件，即不宜过大。计及这一条件，在计及可以认为： 于是则有再按自导纳的定义，上式中得项应为各元件电抗远大于电阻的前提下，除节点i外，其他节点都接地时，由节点i注入无功

5、功率。这功率必远大于正常运行时节点的注入无功功率。雅可比矩阵中的两个子矩阵H，L的元素将有相同的表达式，但阶数不同。展开如下：将上式代入（2.1）式将上面两式左右都乘以下式：上两式可以写成：为了加速收敛，使它们更有效的进行修正，可以考虑在B中尽量去掉那些有功功率及电压向量角度无关或影响较小的因素。为此，我们以电力系统导纳矩阵的虚部作为B，但是去掉了充电电容和变压器非标准变化的影响。具体的说，B的非对角元素和对角元素分别按下式计算：式中：rij和xij分别为支路ij的电阻和感抗。 与迭代用的B按式计算： (4-14)式中bi0为节点i的接地支路导纳。四 程序流程图五 程序clc;clear;g(

6、2,1)=0.0067; b(2,1)=-0.0355; g(2,3)=0.0072; b(2,3)=-0.0472;g(1,2)=0.0067; b(1,2)=-0.0355;g(1,3)=0.0031; b(1,3)=-0.0204;g(3,2)=0.0072; b(3,2)=-0.0472;g(3,1)=0.0031; b(3,1)=-0.0204;delt(1)=0.0;delt(2)=0.0; u(2)=220;u(1)=220; p(2)=117;q(2)=-22; p(1)=-440;q(1)=-136;k=0;precision=1;delt,k,u;N1=2; for m=1

7、:N1+1 for n=1:N1+1 if m=n G(m,m)=g(m,1)+g(m,2)+g(m,3); B(m,m)=b(m,1)+b(m,2)+b(m,3); else G(m,n)=-g(m,n); B(m,n)=-b(m,n); end endendY=G+1i*B;B1=B(1:2,1:2); B11=B(1:2,1:2); while precision0.00001 delt(3)=0;u(3)=238; for m=1:N1 for n=1:N1+1pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-

8、delt(n); end pp(m)=p(m)-sum(pt); end for m=1:N1 PU(m)=pp(m)/u(m); PU1(m)=pp(m)/(u(m)2); end ddelt=-inv(B1)*PU1; precision=max(abs(PU); for n=1:N1 delt(n)=delt(n)+ddelt(n); end for m=1:N1 for n=1:N1+1 qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); end qq(m)=q(m)-sum(qt); end

9、for m=1:N1 QU(m)=qq(m)/u(m); end uu=-inv(B11)*QU; precision=max(abs(QU); for n=1:N1 u(n)=u(n)+uu(n); end k=k+1; k,delt,u end for n=1:N1+1 U(n)=u(n)*(cos(delt(n)+j*sin(delt(n); endfor m=1:N1+1I(m)=Y(3,m)*U(m);endS3=U(3)*sum(conj(I); for m=1:N1+1 for n=1:N1+1S(m,n)=U(m)*(conj(U(m)-conj(U(n)*conj(-Y(m,

10、n); endend YB1B11ppqq udeltUkSS3 六 结果分析delt = 0 0k = 0k = 1delt = -1.8560e-001 -4.3317e-002 0u = 2.0974e+002 2.2521e+002 2.3800e+002k = 6delt = -1.7888e-001 -4.4556e-002 0u = 2.0949e+002 2.2498e+002 2.3800e+002Y = Columns 1 through 2 Column 3 9.8000e-003 -5.5900e-002i -6.7000e-003 +3.5500e-002i -3.1

11、000e-003 +2.0400e-002i -6.7000e-003 +3.5500e-002i 1.3900e-002 -8.2700e-002i -7.2000e-003 +4.7200e-002i -3.1000e-003 +2.0400e-002i -7.2000e-003 +4.7200e-002i 1.0300e-002 -6.7600e-002i B1 = -5.5900e-002 3.5500e-002 3.5500e-002 -8.2700e-002B11 = -5.5900e-002 3.5500e-002 3.5500e-002 -8.2700e-002pp = -5.

12、3147e-004 1.1721e-003qq = -4.2952e-004 6.1704e-005u = 2.0949e+002 2.2498e+002 2.3800e+002delt = -1.7888e-001 -4.4556e-002 0U = Columns 1 through 2 Column 3 2.0614e+002 -3.7274e+001i 2.2476e+002 -1.0021e+001i 2.3800e+002 k = 6S = Columns 1 through 2 Column 3 0 -2.4298e+002 -5.7879e+001i -1.9702e+002

13、-7.8121e+001i 2.5028e+002 +9.6548e+001i 0 -1.3328e+002 -1.1855e+002i 2.0448e+002 +1.2716e+002i 1.3526e+002 +1.3156e+002i 0 S3 = 3.3974e+002 +2.5873e+002i七 设计总结电力系统潮流计算是研究电力系统稳态运行情况的计算方法，它根据给定的运行条件及系统接线情况等确定整个电力系统的运行状态、母线电压、 元件中流过的功率、系统的功率损耗等。其主要目的在于计算电力网络中设备和用户运行状态的参数。 计算的结果可用来考查所研究的运行状态下，网络特性对各种电力设

14、备和用电设备的适用性。在这次课程设计中，通过积极查阅各种相关资料，再结合老师课堂上讲解的相关知识，对于用PQ法计算潮流分布有了更深的认识。P-Q分解法是牛顿-拉夫逊潮流计算的一种简化方法，改进和提高的运算速度。在拿到老师给的课题之后，跟小组同学经过认真的讨论、交流，画出等值电路模型，建立了修正方程，运用MATLAB软件编制了潮流计算的程序。其中也遇到了很多困难，程序也出现了很多错误，但在经过我们努力的修改，相互配合之下，终于成功的得出最终的结论，收获了很多。非常感谢老师的细心指导和我们这个小组成员的不懈努力。 八 参考文献1.陈珩 电力系统温态分析（第三版） 中国电力出版社2.中国自动化网3.电力系统温态分析 水利水电出版社