潮流上机实验报告.docx
潮流上机实验报告电力系统潮流上机课程设计报告 院系 班级： 学号： 学生姓名： 指导教师： 设计周数 成绩： 日期：年月日课程课程设计报告 一、课程设计的目的与要求 培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识 二、设计正文 1.手算： 要求应用牛顿-拉夫逊法或P-Q分解法手算求解，要求精度为0.001MW。 &=1.0Ð0°，节点2为PQ节点，负荷功率S=0.8+j0.6，节节点1为平衡节点，电压U12点3是PV节点，P3=0.4,U3=1.1，两条支路分别为Z13=0.01+j0.04，Z12=0.05+j0.2，对地支路y30=j0.33。 1 课程课程设计报告 2.计算机计算：编写潮流计算程序，要求如下： 2.1据给定的潮流计算任务书整理潮流计算的基础数据：节点的分类，线路模型，等值变压器模型，电压等级的归算，标幺值的计算； 2.2基础数据的计算机存储：节点数据，支路数据； 2.3用牛顿-拉夫逊法计算； 2.4根据所选潮流计算方法画流程图，划分出功能模块，有数据输入模块，导纳阵形成模块，解线性方程组模块，计算不平衡功率模块，形成雅可比矩阵模块，解修正方程模块，计算线路潮流，网损，PV节点无功功率和平衡节点功率，数据输出模块； 2.5据上述模块编制程序并上机调试程序，得出潮流计算结果； 2.6源程序及其程序中的符号说明集、程序流图 简单系统如下图所示，支路数据如下： z12=0.1+j0.41，z13=j0.3，z14=0.12+j0.5，z24=0.08+j0.40 ，y10,4=y40,1=j0.0192，y20,4=y40,2=j0.01413 y10,2=y20,1=j0.01528k=1.1 节点数据如下： &=-0.30-j0.18，S&=-0.55-j0.13， S21&=0.5，U3=1.10，U4=1.05Ð0o S31)节点导纳阵 #include <stdio.h> #include <math.h> #include <fstream.h> #include "LF.h" /form node conductance matrix intMakeY( intnB, intnL, Line* sL, double* YG, double* YB ) inti,j,l; double r,x,d1,g,b,t; for(i=0;i<nB;i+) for(j=0;j<nB;j+) YGij=0.0; YBij=0.0; for(i=0;i<nL;i+) r=sLi.R; x=sLi.X; g=r/(r*r+x*x); b=-x/(r*r+x*x); switch(sLi.Type) case 1:/Line 2 课程课程设计报告 break; case 2:/Transformer b*=1/sLi.K; break; g*=1/sLi.K; YGsLi.NumIsLi.NumI+=g; YGsLi.NumJsLi.NumJ+=g; YGsLi.NumIsLi.NumJ-=g; YGsLi.NumJsLi.NumI-=g; YBsLi.NumIsLi.NumI+=b+sLi.B; YBsLi.NumIsLi.NumJ-=b; YBsLi.NumJsLi.NumI-=b; YBsLi.NumJsLi.NumJ+=b+sLi.B; printf("实部：n"); for(i=0;i<nB;i+) for(j=0;j<nB;j+) printf("%lft",YGij); printf("n"); for(i=0;i<nB;i+) for(j=0;j<nB;j+) printf("%lft",YBij); printf("n"); printf("虚部：n"); /* Check the Y matrix */ ofstreamfout("out.txt"); fout<< "-Y Matrix-" <<endl; for(i=0;i<nB;i+) for(j=0;j<nB;j+) fout<< YGij << "+j" << YBij << "t" fout<<endl; fout.close; return 0; 2）计算功率不平衡量 #include <stdio.h> #include <math.h> #include <fstream.h> #include "LF.h" /form delta p and delta q intCalDeltaPQ( intnpv, intnpq, Bus* bus, double* YG, double* YB, int* p_Jtobus, double* deltaf ) intk,i,j; for(k=0;k<npv+npq*2;k+) i=p_Jtobusk; if(k<npv) 3 课程课程设计报告 deltafk=busi.GenP-busi.LoadP; for(j=0;j<npv+npq+1;j+) deltafk-=busi.Volt*busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Phase-busj.Phase); printf("PV节点%d的有功功率是%lfn",i,deltafk); if(k<npq+npv&&k>=npv) deltafk=busi.GenP-busi.LoadP; for(j=0;j<npv+npq+1;j+) deltafk-=busi.Volt*busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Phase-busj.Phase); printf("PQ节点%d的有功功率是%lfn",i,deltafk); if(k<npq*2+npv&&k>=npv+npq) deltafk=busi.GenQ-busi.LoadQ; for(j=0;j<npv+npq+1;j+) deltafk-=busi.Volt*busj.Volt*(YGij*sin(busi.Phase-busj.Phase)-YBij*cos(busi.Phase-busj.Phase); 3）雅各比矩阵的计算 /* Purpose: for undergraduate course Task: Load Flow Copyright NCEPU, Liu Chongru */ #include <stdio.h> #include <math.h> #include <fstream.h> #include "LF.h" /form Jacobian matrix intFormJacobian( intnpv, intnpq, Bus* bus, double* YG, double* YB, int* p_Jtobus, double* Jac ) intnp = npv+npq,j,k,i,m; /TODO for(k=0;k<npv+npq*2;k+) i=p_Jtobusk; qi=0; 4 printf("PQ节点%d的无功功率是 %lfn",i,deltafk); return 0; double a14,q14; ai=0; 课程课程设计报告 if(k<np)/H N if(j!=i) for(j=0;j<np+1;j+) ai+=busj.Volt*(YGij*sin(busi.Phase-busj.Phase)-YBij*cos(busi.Phase-busj.Phase); qi+=busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Phase-busj.Phase); j=p_Jtobusm; if(j!=i) for(m=0;m<npv+npq*2;m+) if(m<np) Jackm=busi.Volt*busj.Volt*(YGij*sin(busi.Phase-busj.Phase)-YBij*cos(busi.Phase-busj.Phase);/Form H else Jackm=busi.Volt*busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Phase-busj.Phase);/Form N else for(j=0;j<np+1;j+) if(m<np) Jackm=-busi.Volt*ai;/Form H else if(j=i) else Jackm=busi.Volt*qi+2*busi.Volt*busi.Volt*YGij;/Form N if(j!=i) ai+=busj.Volt*(YGij*sin(busi.Phase-busj.Phase)-YBij*cos(busi.Phase-busj.Phase); qi+=busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Phase-busj.Phase); j=p_Jtobusm; if(j!=i) for(m=0;m<npv+npq*2;m+) if(m<np) Jackm=-busi.Volt*busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Phase-busj.Phase); /Form J else Jackm=busi.Volt*busj.Volt*(YGij*sin(busi.Phase-busj.Phase)-YBij*cos(busi.Phase-busj.Phase); /Form L else if(j=i) if(m<np) 5 课程课程设计报告 Jackm=busi.Volt*qi; Jackm=busi.Volt*ai-2*busi.Volt*busi.Volt*YBij; else for(i=0;i<np+npq;i+) for(int j=0;j<np+npq;j+) printf("%d %d %f ",i,j,Jacij); printf("n"); /Output the matrix to check the Jacobian matrix ofstreamfout("out.txt",ios:app); fout<< "-Jacobian Matrix-" <<endl; for(i=0; i<np+npq;i+ ) for(j=0; j<np+npq; j+ ) fout<<Jacij << "t" fout<<endl; fout.close; return 0; 4）线路损耗 /8.calculate the power flow double* p_Pij, *p_Qij, *p_Pji, *p_Qji; p_Pij = new doublenL; p_Qij = new doublenL; p_Pji = new doublenL; p_Qji = new doublenL; int x1,x2; for( i=0; i<nL; i+ ) x1=linei.NumI; x2=linei.NumJ; if(linei.Type=1) p_Piji=busx1.Volt*busx1.Volt*(-YGx1x2)-busx1.Volt*busx2.Volt*(-YGx1x2)*cos(busx1.Phase-busx2.Phase)+(-YBx1x2)*sin(busx1.Phase-busx2.Phase); p_Qiji=-busx1.Volt*busx1.Volt*(linei.B+(-YBx1x2)-busx1.Volt*busx2.Volt*(-YGx1x2)*sin(busx1.Phase-busx2.Phase)-(-YBx1x2)*cos(busx1.Phase-busx2.Phase); p_Pjii=busx2.Volt*busx2.Volt*(-YGx2x1)-busx2.Volt*busx1.Volt*(-YGx2x1)*cos(busx2.Phase-busx1.Phase)+(-YBx2x1)*sin(busx2.Phase-busx1.Phase); p_Qjii=-busx2.Volt*busx2.Volt*(linei.B+(-YBx2x1)-busx2.Volt*busx1.Volt*(-YGx2x1)*sin(busx2.Phase-busx1.Phase)-(-YBx2x1)*cos(busx2.Phase-busx1.Phase); else p_Piji=busx1.Volt*busx1.Volt*(-YGx1x2)/linei.K-busx1.Volt*busx2.Volt*(-YGx1x2)*cos(busx1.Phase-busx2.Phase)+(-YBx1x2)*sin(busx1.Phase-busx2.Phase); p_Qiji=-busx1.Volt*busx1.Volt*(-YBx1x2)/linei.K+linei.B)-busx1.Volt*bu6 课程课程设计报告 sx2.Volt*(-YGx1x2)*sin(busx1.Phase-busx2.Phase)-(-YBx1x2)*cos(busx1.Phase-busx2.Phase); p_Pjii=busx2.Volt*busx2.Volt*(-YGx2x1*linei.K)-busx2.Volt*busx1.Volt*(-YGx2x1)*cos(busx2.Phase-busx1.Phase)+(-YBx2x1)*sin(busx2.Phase-busx1.Phase); p_Qjii=-busx2.Volt*busx2.Volt*(-YBx2x1)*linei.K+linei.B)-busx2.Volt*busx1.Volt*(-YGx2x1)*sin(busx2.Phase-busx1.Phase)-(-YBx2x1)*cos(busx2.Phase-busx1.Phase); /p and q of PH bus and PV bus int s=0; double p9,q9,Ps9,Qs9,PS=0,QS=0; pi=0; qi=0; for( i=0; i<nB; i+ ) for(int j=0; j<nB; j+ ) pi+=(busj.Volt*(YGij)*cos(busj.Phase)-busj.Volt*(YBij)*sin(busj.Phase); qi-=(busj.Volt*(YGij)*sin(busj.Phase)+busj.Volt*(YBij)*cos(busj.Phase); Psi=busi.Volt*cos(busi.Phase)*pi-busi.Volt*sin(busi.Phase)*qi; Qsi=busi.Volt*cos(busi.Phase)*qi+busi.Volt*sin(busi.Phase)*pi; for(i=0;i<nB;i+) PS+=Psi; QS+=Qsi; printf("PS=%7.7f,QS=%7.7fn",PS,QS); /loss doublePloss=0, Qloss=0; for( i=0; i<nB; i+ ) Ploss+=p_Piji+p_Pjii; Qloss+=p_Qiji+p_Qjii; 5）程序流图如下 启动 输入原始数据 形成节点导纳矩阵 设非平衡节点电压初值ei(0)fi(0); 令迭代次数k=0 对PQ节点计算DPi(k).DQi(k)对PV节点计算DPi(k),DUi2(k) 令节点号i=1 7 课程课程设计报告 迭代是否收敛 解修正方程，由DPi(k).DQi(k)DUi2(k)及雅客比矩阵用高斯法求各节点的电压增量De(k)i雅客比矩阵是否形成，i>n? 计算雅客比矩阵各元素(k)(k)(k)(k)(k) Hij,Nij,Jij,L(ijk),Rij,Sij增加迭代次数k=k+1 增加节点号i=i+1 计算节点的新电压 ,Dfi(k)ei(k+1)=ei(k)+Dei(k) fi(k+1)=fi(k)+Dfi(k)求出De(k)max,Df(k)maxDe(k),Df(k)maxmax<e计算平衡节点的功率及线路功率 停止 6）得到的数据 -Y Matrix- 0+j-17.3611 0+j0 0+j0 0+j17.3611 0+j0 0+j-16 0+j17.3611 0+j0 0+j0 0+j0 0+j0 0+j0 0+j0 0+j0 0+j0 0+j0 0+j0 0+j0 0+j16 0+j0 0+j0 0+j-17.0648 0+j0 0+j0 0+j0 0+j0 0+j0 0+j17.0648 0+j0 0+j0 3.30738+j-39.3089 -1.36519+j11.6041 -1.94219+j10.5107 0+j0 0+j0 0+j0 0+j0 0+j0 0+j0 -1.36519+j11.6041 2.55279+j-17.3382 0+j0 -1.1876+j5.97513 0+j0 0+j0 0+j0 0+j0 0+j0 -1.94219+j10.5107 0+j0 3.2242+j-15.8409 0+j0 0+j0 -1.28201+j5.58824 0+j0 0+j16 0+j0 0+j0 -1.1876+j5.97513 0+j0 2.80473+j-35.4456 -1.61712+j13.698 0+j0 0+j0 0+j0 -1.28201+j5.58824 0+j0 -1.15509+j9.78427 2.4371+j-32.1539 -16.4 0 0 0 0 0 0 0 0 0 0 8 0+j0 0+j0 0+j0 0+j0 0+j0 0+j0 -1.61712+j13.698 2.77221+j-23.3032 -1.15509+j9.78427 0+j0 0+j0 0+j17.0648 16.4 0 0 0 0 0 0 0 0 0 0 -Jacobian Matrix- 0 0 0 0 0 0 0 0 0 0 0 0 -1.1876 0 17.4915 0 -17.4915 0 40.1703 -11.6041 -10.5107 0 -11.6041 17.5792 0 3.30738 -1.36519 -1.94219 0 -1.36519 2.55279 0 -5.97513 0 课程课程设计报告 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -10.5107 0 0 0 0 0 0 0 16.0989 0 0 -5.58824 -1.94219 0 0 0 0 0 3.2242 0 0 -1.28201 -16.4 -5.97513 0 36.0731 -13.698 0 -1.1876 0 2.80473 -1.61712 0 -13.698 23.4822 -9.78427 0 -5.58824 0 0 0 0 -1.61712 2.77221 -1.15509 -1.28201 0 -1.15509 2.4371 0 0 0 -17.4915 0 -9.78427 32.864 0 -3.30738 1.36519 1.94219 0 1.36519 -2.55279 0 1.94219 0 0 0 0 0 0 0 0 1.1876 0 -3.2242 0 0 38.4474 -11.6041 -10.5107 0 -11.6041 17.0972 0 0 -5.97513 0 0 0 0 0 0 15.5829 0 0 1.1876 0 -5.97513 0 1.28201 -10.5107 0 -5.58824 -2.80473 1.61712 0 34.8181 -13.698 0 1.61712 -2.77221 1.15509 0 1.15509 -2.4371 0 0 0 0 0 0 0 0 -13.698 23.1242 -9.78427 -9.78427 31.4437 0 0 0 0 1.28201 0 0 0 0 0 -5.58824 0 1.68793 0 0 0 -Jacobian Matrix- 16.9269 0 0 0 0 0 0 0 0 -16.9269 0 0 0 18.1691 0 -18.1691 0 0.883627 -2.10262 0 0 0 41.9297 -12.1301 -11.1536 0 -12.0455 18.0609 0 -11.0484 0 0 0 -6.36224 0 0 0 16.8144 0 0 3.54272 -1.0628 -1.76646 0 -1.78138 1.30819 0 0 0 -0.357534 0 0 -2.509 2.42598 0 -6.01539 0 -5.76607 -2.33608 0 -1.97778 -16.9269 0 0 0 0 0 0 0 0 0 0 0 37.9476 -14.6585 0 3.05959 -0.930027-14.4721 24.8873 -10.4152 0 -6.05157 0 1.86088 -1.47389 0 -0.9919730 0 0 -18.1691 0 2.6627 0 0 0 0 -10.4721 34.6928 0 0 0 0 -0.733327 -3.52149 1.0628 1.76646 0 1.78138 -3.884 2.33608 0 0 0 0 0 0 0 0 0 -4.31386 0 2.509 0 42.0299 -12.1301 -11.1536 0 -12.0455 17.2037 0 0 0 -6.36224 0 0 0 16.2993 0 0 2.10262 0 -6.01539 0 1.97778 -11.0484 0 -5.76607 1.68793 0 0.357534 0 -2.97549 0.930027 0 -3.98289 1.47389 0 38.3226 -14.6585 0 -14.4721 24.2355 -10.4152 -6.05157 0 0 -10.4721 34.8585 0 0 0 0.883627 0 0 0.733327 0 0 0 0.991973 -2.60893 0 0 0 0 0 0 0 0 0 0 -Jacobian Matrix- 16.7457 0 0 0 0 0 0 0 0 -16.7457 0 0 0 1.63043 0 0 0 18.0388 0 -18.0388 0 0.850196 -2.0217 0 0 0 41.3695 -11.8919 -10.9686 0 -11.8057 17.6918 0 -10.8651 0 0 0 -6.21183 0 0 0 16.5476 0 0 3.48069 -1.02775 -1.73712 0 -1.7602 1.28091 0 0 0 -0.382862 0 0 2.40655 0 -5.8861 0 -5.68251 -2.29737 0 -1.91027 -16.7457 0 0 0 0 0 0 0 0 0 0 0 37.3041 -14.3465 0 2.95313 -0.937485-14.1704 24.4052 -10.2348 0 -5.94693 0 -2.42909 1.86079 -1.43353 -0.757656 0 -0.9892080 0 0 -18.0388 0 2.59847 0 0 0 0 -10.2872 34.273 0 0 0 0 -3.48089 1.02775 1.73712 0 1.7602 -3.78189 0 2.29737 0 0 0 0 0 0 0 0 0.382862 0 -4.20764 0 0 41.3703 -11.8919 -10.9686 0 -11.8057 16.6941 0 0 0 -6.21183 0 0 0 15.9488 0 0 2.0217 0 -5.8861 0 1.91027 -10.8651 0 -5.68251 1.63043 0 -2.95077 0.937485 0 37.3083 -14.3465 0 2.42909 -3.86262 1.43353 0 0.757656 0 0 0 -14.1704 23.7059 -10.2348 -5.94693 0 0 -10.2872 34.2743 0.850196 0 0 0.989208 -2.59706 0 0 0 9 -Jacobian Matrix- 16.7435 0 0 -16.7435 0 0 0 0 0 0 0 1.63 0 0 0 18.0374 0 -18.0374 0 0.85 课程课程设计报告 0 0 0 0 0 0 41.3625 -11.8888 -10.9664 0 -11.8026 17.6871 0 -10.8628 0 0 0 -6.20987 0 0 0 16.5444 0 0 0 0 0 3.48016 -1.02713 -1.73662 0 -1.76008 1.28053 0 0 0 -0.383399 0 0 2.40632 0 0 0 0 0 0 -5.8845 0 -2.02045 0 -5.68158 -2.29703 0 -1.90929 -16.7435 0 0 0 0 0 0 0 0 0 0 0 37.296 -14.3426 0 2.95114 -0.937742-14.1667 24.3994 -10.2326 0 -5.94567 0 -2.42794 1.86097 -1.43302 -0.758139 0 -0.9892030 0 0 -18.0374 0 2.59734 0 0 0 0 -10.285 34.2681 0 0 0 0 -3.48016 1.02713 1.73662 0 1.76008 -3.78053 0 2.29703 0 0 0 0 0 0 0.383399 0 -4.20632 0 0 41.3625 -11.8888 -10.9664 0 -11.8026 16.6871 0 0 -6.20987 0 0 0 0 15.9444 0 0 2.02045 0 -5.8845 0 1.90929 -10.8628 0 -5.68158 1.63 0 -2.95114 0.937742 0 37.296 -14.3426 0 2.42794 -3.86097 1.43302 0 0.989203 -2.59734 0 -14.1667 23.6994 -10.2326 -10.285 34.2681 0.85 0 0.758139 0 -5.94567 0 -iteration- iteration = 4 -voltage magnitude and angle- 1.04 0 1.025 1.025 0 0.161967 9.28001 0.0814153 4.66476 -2.21679 -3.98881 -3.6874 0.727537 1.96672 1.02579 -0.0386902 0.995631 -0.0696178 1.01265 -0.0643572 1.01588 0.0126979 1.03235 0.0343257 1 2 3 4 5 6 7 8 9 1.02577 0.064921 3.7197 -bus P and Q- 0.71641 0.270459 1.63 0.0665366 0.85 -0.108597 0 0 -0.5 -1.25 -0.9 -0.3 0 0 0 0 -1 -0.35 -line flow- -NUM-i-j-begin-end- 1 2 3 4 5 6 7 8 4 7 9 7 9 7 9 5 1 2 3 8 8 5 6 4 -0.71641+j-0.239232 -0.85+j0.149553 0.71641+j0.27046 -1.63+j0.0917816 1.63+j0.0665365 0.85+j-0.108597 0.763799+j-0.00797398 -0.759046+j-0.107041 0.241834+j0.0311946 -0.240954+j-0.242958 0.866201+j-0.0838079 -0.843202+j-0.113128 0.608166+j-0.180748 -0.594627+j-0.134566 10 -0.406798+j-0.386872 0.409373+j0.228931 课程课程设计报告 9 6 4 -0.305372+j-0.165433 0.307037+j0.0102993 -Ploss and Qloss
