电力网潮流电压计算例题与matlab程序

1电力网潮流电压计算例题与MATLAB程序编写佘名寰本文介绍了两个电力网潮流计算例题。一例为5个节点5条支路；另一例有6个节点7条支路，有PQ负荷节点也有PV发电机节点，变压器支路标么变比不为1。本文给出了完整的计算程序和计算结果。程序包括计算网络节点导纳矩阵和潮流电压两部分，程序中引用的M函数见参考文献1。例题选自研究生教材，比较典型，计算结果与教材核对基本一致。本文可供电力专业师生和基层技术人员潮流计算时参考。【例1.3】如图1-4所示，五个节点三条母线简单系统，系统参数均用标么值（pu.）表示。节点4、5通过变压器与节点2、3相连，变比为1.05:1，节点1、2、3为PQ负荷节点，节点4为PV发电机节点，节点5为参考节点，节点4、5电压均为1.05，图中箭头表示各支路的正方向，各个节点功率如表2.1所示,计算各节点电压、发电机无功功率和系统潮流分布：图1-4[例1.3]网络接线图表1-1母线有功负载无功负载发电机有功发电机无功1PD1=1.6QD1=0.82PD2=2.0QD2=1.03PD3=3.7QD3=1.345PG4=5.0PG5=？QG4=未定QG5=？4Z23Z13Z21YC2YC3YC1Z42Z531:1.051.05:142315S3S2S1Z21=0.04+J0.25Z23=0.08+J0.30Z13=0.1+J0.35Z42=J0.015Z53=J0.03YC1=J0.25YC2=J0.50YC3=J0.25S1=1.6+J0.8S2=2+J1S3=3.7+J1.3P4=5.0V4=1.05V5=1.05P4=5.0V4=1.05vV5=1.052【例1.3】程序说明：本例题与参考文献2的【例2.3.1】不同之处在于有两个变比不为1的变压器支路，并且引入了PV发电机节点。由于网络结构与例图1-1完全相同，网络节点导纳矩阵已在文献1之1.3节求出不再赘述。潮流计算程序参考【例2.3.1】节点电压符号矩阵：一共五个节点，独立节点数N1=4，PV节点数N2=1节点电压幅值符号矩阵为u=sym('[u1,u2,u3,u4,u5]');节点电压相角符号矩阵为delt=sym('[d1,d2,d3,d4,d5]');节点功率给定值为：三个PQ节点，一个PV节点p=[-1.6,-2,-3.7,5];q=[-0.8,-1.0,-1.3];【例1.3】源程序NU114N.m%*******NU114N.m,example2-3-2,fig2-2********%%ThefollowingProgramforloadflowcalculationisbasedonMATLAB2007clear%bus1,2,3isPQbus,bus4isPVbus,bus5isslackbusglobalNpNbbtpmk%Npisnumberofnodepoint,Nbisnumberofbraches,%pthetransformerturnsratio,withoff-nominaltap-setting,fig2-4,PSCLFp6,%btistwo-dimensionalarray,lineoneispoint1,linetwoispoint2,%misnumberofnodepoint,kisnumberoftransformersNp=5;Nb=5;bt=[2,3;4,5];p=[1/1.05,1/1.05];m=5;k=2;%nstart--thestartpointofbranches,nend--theendpoint,%mm--networkincidencematrixnstart=[2,2,1,4,5];nend=[1,3,3,2,3];mm=ffm(nstart,nend);%zb1,theseriesimpedancesoftransmissionline%yb1,theseriesadmittancesoftransmissionlinezb1r=[0.04,0.08,0.1,0.0,0.0];zb1i=[0.25,0.30,0.35,0.015,0.03];zb1=zb1r+zb1i*j;yb1=zb1.^(-1);yb=diag(conj(yb1'));y=mm*yb*(mm)';%yb0,theshuntadmittancesoftransmissionlineyb0i=[0.25,0.25,0.0,0.0,0.0];yb0=0+yb0i*j;y0=diag(conj(yb0'));yg=mm*y0*(mm)';yn=diag(diag(yg));yy=yn-yg;y=y+yn;3%yn0,thenode-admittanceofshuntcapacitor%yn0i=[];%yn0=0+yn0i*j;%ync=diag(conj(yn0'));%y=y+ync;[Y,YY]=fdt1(y,yy);G=real(Y);B=imag(Y);%YYistheshundadmittancesofthelineandtransformersateachend,%Yisnode-admittancematrixfornetwork%*********************************************************************u=sym('[u1,u2,u3,u4,u5]');delt=sym('[d1,d2,d3,d4,d5]');%u--nodevoltagemagnitude,delt--anglep=[-1.6,-2,-3.7,5];q=[-0.8,-1.0,-1.3];k=0;precision=1;N1=4;%theN1istheamountofthePQandPVbusN2=1;%theN2istheamountofthePVbusform=1:N1forn=1:N1+1pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));endpp(m)=p(m)-sum(pt);endform=1:N1-N2forn=1:N1+1qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));endqq(m)=q(m)-sum(qt);endJ1=jacobian(pp,['d1,d2,d3,d4,u1,u2,u3']);J2=jacobian(qq,['d1,d2,d3,d4,u1,u2,u3']);J=vertcat(J1,J2);uu=[1.0,1.0,1.0,1.05,1.05];dd=[0,0,0,0,0];whileprecision0.00001u1=uu(1);u2=uu(2);u3=uu(3);u4=uu(4);u5=uu(5);d1=dd(1);d2=dd(2);d3=dd(3);d4=dd(4);d5=dd(5);form=1:N1PP(m)=eval(pp(m));end4form=1:N1-N2PP(N1+m)=eval(qq(m));endJJ=eval(J);du=-inv(JJ)*PP';precision=max(abs(du));forn=1:N1dd(n)=dd(n)+du(n);endforn=1:N1-N2uu(n)=uu(n)+du(N1+n);endk=k+1;endk-1,dd*180/pi,uu%***********************************************%thefollowingprogramisusedtocalculatetheSmandSmnforn=1:N1+1U(n)=uu(n)*(cos(dd(n))+j*sin(dd(n)));endUm=conj(U');I=Y*Um;Sm=diag(Um)*conj(I)form=1:N1+1forn=1:N1+1Smn(m,n)=U(m)*(conj(U(m))-conj(U(n)))*conj(-Y(m,n))+U(m)*conj(U(m)*YY(m,n));endendSmn【例1.3】潮流计算结果：迭代次数ans=4节点1至5电压相角（度）ans=-4.778517.8535-4.281921.84330节点1至5电压幅值（pu）uu=0.86221.07791.03641.05001.0500节点1至5功率（pu）Sm=-1.6000-0.8000i5-2.0000-1.0000i-3.7000-1.3000i5.0000+1.8131i2.5794+2.2994i支路功率（行列编号均为节点1至5）Smn=0-1.4662-0.4091i-0.1338-0.3909i001.5845+0.6726i01.4155-0.2443i-5.0000-1.4282i00.1568+0.4713i-1.2774+0.2032i00-2.5794-1.9745i05.0000+1.8131i000002.5794+2.2994i00【例1.4】如图1-5所示，六个节点网络。节点2、1间和节点4、3间通过变压器相连，变比分别为0.909:1和0.976:1，节点1、2、3、4为PQ负荷节点，节点5为PV发电机节点，节点电压为1.10;节点6为参考节点,节点电压为1.05。输电线路和变压器数据见表2.2,各个节点功率如表2.3所示,计算各节点电压、发电机无功功率和系统潮流分布：图1-5[例1.4]网络接线图表1-2输电线路和变压器数据汇编线路首端节点编号线路末端节点编号充电功率换算导纳BC/2线路阻抗换算电导G线路阻抗换算电纳B变压器变比620.00700.2791-1.2910640.00900.2170-0.9137510.00.2224-0.3230530.00.2383-0.6542210.00.0-3.75940.909240.00760.2771-1.1625430.00.0-1.66670.976GG①⑤②③④⑥87563214910V5=1.10V6=1.056表1-3各节点功率母线有功负载无功负载发电机有功发电机无功1PD1=0.275QD1=0.0652PD2=0.0QD2=0.03PD3=0.15QD3=0.09456PD4=0.25QD4=0.025PG5=0.25PG6=？QG4=未定QG6=？【例1.4】程序说明：本例题计算程序可分为形成节点导纳矩阵和计算节点电压潮流两部分。计算节点导纳矩阵参考文献1第1.3节。这里有6个节点7条支路，节点编号原则是先PQ节点，后PV节点，平衡节点编号最大，支路序号如图所示。各支路的首末端节点编号用数组nstart和nend表示，各支路导纳用数组yb1r,yb1i表示，各支路两端充电电容导纳为数组yb0i:Np=6;Nb=7;nstart=[6,6,2,3,5,5,1];nend=[4,2,4,4,3,1,2];yb1r=[0.2170,0.2791,0.2771,0.0,0.2383,0.2224,0.0];yb1i=[-0.9137,-1.2910,-1.1625,-1.6667,-0.6542,-0.3230,-3.7594];yb0i=[0.0099,0.0070,0.0076,0.0,0.0,0.0,0.0];两台变压器跨接在节点2-1和4-3之间，参见图1-2，i节点对应2、4，j节点对应1、3，变比为1：1/0.909和1:1/0.976,程序中用数组bt[]和p[]表示，m=6,k=2表明用6个节点两条变压器支路：bt=[2,4;1,3];p=[1/0.909,1/0.976];m=6;k=2;计算潮流电压参照例题[1.3]。这里独立节点数N1=5，PV节点数N2=1，节点电压初始给定值幅值为数组uu,相角为dd，节点有功功率为p,无功功率为q:uu=[1.000,