《电力系统潮流计算(共18页).docx》由会员分享，可在线阅读，更多相关《电力系统潮流计算(共18页).docx（18页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上系统的潮流及三相短路电流计算班级：电气1班学号：姓名：杨鹏摘要潮流计算，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。牛顿-拉夫逊法是电力系统潮流计算的常用算法之一，收敛性好，迭代次数少。本文基于MATLAB计算系统的潮流及三相短路电流。关
2、键词：潮流计算 matlab 牛顿-拉夫逊 1电力系统的潮流计算电力系统常规潮流计算的任务是根据给定电网机构、发电计划及负荷分布情况，求出整个电网的运行状态，其中包括各节点母线电压、相角、线路传输的有功功率和无功功率等。在电网的潮流计算中，一般给定的运行参数有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。待求的参数包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。1.1节点的功率方程对n节点电力系统，节点i注入的有功功率Si：极坐标形式的节点功率方程：直角坐标形式的节点功率方程：1.2节点分类:根据实际运行条件，节点可分成三类：PQ节点、PV节点和平衡
3、节点PQ节点：节点注入的有功P和无功Q皆为给定量的节点。一般负荷节点，联络节点和给定有功和无功的发电机节点在潮流计算中都视作PQ节点，PQ节点的节点电压（其幅值U和相角，或其实部e和虚部f）为待求变量。PV节点：节点注入的有功P和无功Q皆为给定量的节点。一般负荷节点，联络节点和给定有功和无功的发电机节点在潮流计算中都视作PQ节点，PQ节点的节点电压（其幅值U和相角，或其实部e和虚部f）为待求变量。平衡节点：平衡节点的节点电压是给定值，对极坐标形式的节点功率方程，平衡节点的电压幅值一般情况下取作U=1.0，相角取作 ，对直角坐标形式的节点功率方程，平衡节点的实部和虚部一般分别取作e=1.0和 f
4、=0.0。计算中所得其它节点电压的相角以平衡节点的相角为参考。平衡节点提供的有功和无功注入除了需要平衡整个电网发电和负荷的不平衡功率，还要平衡整个电网的有功和无功损耗，其值只有在潮流计算后才能确定。潮流计算中原则上可以选择任意发电机节点作平衡节点，但通常以选择容量较大，离负荷中心电气距离较近的发电机节点作平衡节点。1.2电力网络的潮流方程当电力网络的节点类型确定后，便可以列写全网络的潮流方程。1.2.1极坐标形式的潮流方程对PQ，PV节点：对PQ节点：上式可联立表示成非线性向量的形式：1.2. 2直角坐标形式的潮流方程对PQ，PV节点：对PQ节点：还要对所有PV节点增补如下方程，其中，是PV节
5、点给定的电压幅值，上式可联立表示成非线性向量的形式：2牛顿-拉夫逊潮流算法2.1牛顿-拉夫逊算法牛顿法的核心是把非线性方程f(x)=0的求解过程，转变为反复求解对应的线性增量方程，并用线性增量方程的解修正非线性方程解的过程。在x(k)附近将f(x)=0左端向量函数展开成泰勒级数，并略去二阶及以上的高阶项，得到如下线性化方程，称作牛顿法的修正方程。其中，称作非线性向量函数在 x(k)点处的雅可比矩阵。解修正方程组，求得迭代修正量，即：则，第k+1步迭代的解向量判断收敛:，若收敛，则为非线性代数方程组的解，退出迭代。否则输出迭代不收敛信息，退出迭代。2.2 牛顿-拉夫逊算法算法流程3系统的潮流及三
6、相短路电流算例潮流计算：方法不限，计算系统的节点电压和相角。短路电流：4号母线发生金属性三相短路时（zf=0），分别按照精确算法和近似算法计算短路电流、系统中各节点电压以及网络中各支路的电流分布，并对两种情况下的计算结果进行比较。近似算法是指：形成节点导纳矩阵时，所有节点的负荷都略去不计，短路前网络处于空载状态，短路前各节点电压均为1.0。计算步骤：(1) 进行系统正常运行状态的潮流计算，求得(2) 形成不含发电机和负荷的节点导纳矩阵YN；(3) 将发电机表示为电流源()和导纳()的并联组合；节点负荷用恒阻抗的接地支路表示；形成包括所有发电机支路和负荷支路的节点导纳矩阵Y，即在YN中的发电机节
7、点和负荷节点的自导纳上分别增加发电机导纳和负荷导纳（）；(4) 利用，计算节点阻抗矩阵，从而得到阻抗矩阵中的第f列；(5) 利用公式(6-7)或(6-10)计算短路电流；(6) 利用公式(6-8)或(6-11)计算系统中各节点电压；(7) 利用公式(6-9)计算变压器支路的电流；对输电线路利用型等值电路计算支路电流。表1 9节点系统支路参数支路R（p.u.）X（p.u.）B/2(TK) （p.u.）1400.05761.02700.06251.03900.05861.0450.010.0850.088460.0170.0920.079570.0320.1610.153690.0390.170.
8、179780.00850.0720.0745890.01190.10080.1045表2 9节点系统发电机参数发电机编号节点类型PG（p.u.）VG（p.u.）（p.u.）（p.u.）1平衡1.040.31.1372PV1.631.0250.31.2113PV0.851.0250.31.043表3 9节点系统负荷参数节点编号节点类型Pi（p.u.）Qi（p.u.）4PQ005PQ1.250.56PQ0.90.37PQ008PQ10.359PQ00表4 短路前9节点潮流计算结果（供参考）节点编号电压幅值(p.u.)电压相角(rad)节点有功(p.u.)节点无功(p.u.)11.040.00.71
9、640.270521.0259.281.630.066531.0254.660.85-0.108641.0258-2.2170050.9956-3.989-1.25-0.561.0127-3.687-0.90-0.371.02583.7200081.01590.728-1.0-0.3591.03241.967003.1程序运行结果迭代次数：12节点导纳矩阵:0-17.36i000+17.36i0000000-16i00000+16i00000-17.06i000000+17.064i0+17.36i003.30-39i-1.36+11.6i-1.94+10.51i000000-1.36+11.
10、6i2.55-17.33i0-1.18+5.97i00000-1.94+10.51i03.22-15.84i00-1.28+5.58i00+16i00-1.18+5.97i02.8-35.44i-1.6+13.69i0000000-1.6+13.69i2.77-23.30i-1.15+9.78i000+17.36i00-1.28+5.58i0-1.15+9.78i2.43-32.15i电压幅值U 1.0000 1.0000 1.0000 1.0704 0.8482 1.6001 1.6665 1.3389 1.0702电压相角a = 0 9.8974 4.4486 -2.9480 -3.538
11、6 -3.1649 3.2647 0.3083 1.7942节点有功P =0.5904 1.4109 0.2977 0.1643 -1.1511 -0.0598 0.9936 -1.2437 -0.1567节点无功Q =0.4458 0.5927 -0.7510 0.9473 -0.8220 -0.5932 0.7520 -0.2676 0.6408输出修正后的发电机节点和负荷节点的自导纳 1点的自导纳= 0.0000 -20.4443i 2点的自导纳= 0.0000 -19.3332i 3点的自导纳= 0.0000 -20.5585i 5点的自导纳= 3.2561 -17.4375i 6点的
12、自导纳= 4.9498 -16.8225i 8点的自导纳= 3.4869 -23.1242i精确计算短路电流If的模值和角度 模值ans =6.5819 相角ans =-73.8255精确计算短路各节点电压模值ans = 0.6688 0.7274 0.2331 0.0000 0.8082 0.1374 0.5896 0.7166 0.7706近似计算短路电流If的模值和角度 模值ans =6.7770 相角ans =-71.8791近似计算短路各节点电压模值ans = 1.8855 1.8146 1.3719 1.2158 1.6528 1.7934 1.3461 1.2421 1.5864
13、电流幅值误差wuchai0 =0.6910电流相角误差wuchai1 =0.4827电压幅值误差ans = 1.8855 1.8146 1.3719 1.2158 1.6528 1.7934 1.3461 1.2421 1.5864电压相角误差ans = 1.8855 1.8146 1.3719 1.2158 1.6528 1.7934 1.3461 1.2421 1.5864短路时各支路电流 a =i.j. Iij实部 Iij虚部 Iji实部 Iji虚部4.5. -0. 1. 0. -1.94.6. -0. 1. 0. -1.45.7. -0. 1. 0. -1.86.9. -0. 1. 0
14、. -1.07.8. 0. 0. -0. -0.08.9. 0. 0. -0. -0.61.4. 0. -3. -0. 3.62.7. 1. -1. -1. 1.33.9. 0. -1. -0. 1.7i j Iij模值 Iij相角 Iji模值 Iji相角4.5 0. 1. 0. -0.94.6 0. 1. 0. -0.35.7 0. 1. 0. -0.06.9 0. 1. 0. -0.87.8 0. 0. 0. -1.68.9 0. 0. 0. -1.31.4 0. -0. 0. 1.72.7 0. -0. 0. 1.23.9 0. -0. 0. 1.6附件:matlab代码format l
15、ongA=1 4 0 0.0576 0; 2 7 0 0.0625 0; 3 9 0 0.0586 0; 4 5 0.01 0.085 0.088; 4 6 0.017 0.092 0.079; 5 7 0.032 0.161 0.153; 6 9 0.039 0.17 0.179; 7 8 0.0085 0.072 0.0745; 8 9 0.0119 0.1008 0.1045;y=zeros(9,9);for k=1:9 m=A(k,1); n=A(k,2); y(m,m)=y(m,m)+1/(A(k,3)+1i*A(k,4)+ 1i*A(k,5); y(n,n)=y(n,n)+1/(A
16、(k,3)+1i*A(k,4)+ 1i*A(k,5); y(m,n)=y(m,n)-1/(A(k,3)+1i*A(k,4); y(n,m)=y(n,m)-1/(A(k,3)+1i*A(k,4);endU=1.04,1.025,1.025,1,1,1,1,1,1;a=0,0,0,0,0,0,0,0,0;P=1,1.63,0.85,0,-1.25,-0.9,0,-1,0;Q=0,0,0,0,-0.5,-0.3,0,-0.35,0;dP=0,0,0,0,0,0,0,0,0;dQ=0,0,0,0,0,0,0,0,0;G=real(y);B=imag(y);Pt=zeros(1,9);Qt=zeros(
17、1,9);H=zeros(8,8);N=zeros(8,6);M=zeros(6,8);L=zeros(6,6);Ai=zeros(1,9);Bi=zeros(1,9);JJ=zeros(14,14);I=zeros(1,9);k=0;precison=1;while precison0. for m=2:9 for n=1:9 pt(n)=U(m)*U(n)*(G(m,n)*cos(a(m)-a(n)+B(m,n)*sin(a(m)-a(n); end dP(m)=P(m)-sum(pt); end for m=4:9 for n=1:9 qt(n)=U(m)*U(n)*(G(m,n)*si
18、n(a(m)-a(n)-B(m,n)*cos(a(m)-a(n); end dQ(m)=Q(m)-sum(qt); end for m=1:8 for n=1:8 if m=n else H(m,n)=-U(m+1)*U(n+1)*(G(m+1,n+1)*sin(a(m+1)-a(n+1)-B(m+1,n+1)*cos(a(m+1)-a(n+1); end end end for m=1:8 for n=1:6 if m=n+2 elseN(m,n)=-U(m+1)*U(n+3)*(G(m+1,n+3)*cos(a(m+1)-a(n+3)+B(m+1,n+3)*sin(a(m+1)-a(n+3
19、); end end end for m=1:6 for n=1:8 if m+2=n else M(m,n)=U(m+3)*U(n+1)*(G(m+3,n+1)*cos(a(m+3)-a(n+1)+B(m+3,n+1)*sin(a(m+3)-a(n+1); end end end for m=1:6 for n=1:6 if m=n else L(m,n)=-U(m+3)*U(n+3)*(G(m+3,n+3)*sin(a(m+3)-a(n+3)-B(m+3,n+3)*cos(a(m+3)-a(n+3); end end end %雅可比矩阵i不等于j for m=1:8 for n=1:9
20、Ai(n)=U(n)*(G(m+1,n)*sin(a(m+1)-a(n)-B(m+1,n)*cos(a(m+1)-a(n); end H(m,m)=U(m+1)*U(m+1)*B(m+1,m+1)+U(m+1)*sum(Ai); end for m=1:6 for n=1:9 Bi(n)=U(n)*(G(m+1,n)*cos(a(m+1)-a(n)+B(m+1,n)*sin(a(m+1)-a(n); end N(m+2,m)=-U(m+1)*U(m+1)*G(m+1,m+1)-U(m+1)*sum(Bi); end for m=1:6 for n=1:9 Bi(n)=U(n)*(G(m+3,n
21、)*cos(a(m+3)-a(n)+B(m+3,n)*sin(a(m+3)-a(n); end M(m,m+2)=-U(m+3)*U(m+3)*G(m+3,m+3)-U(m+3)*sum(Bi); end for m=1:6 for n=1:9 Ai(n)=U(n)*(G(m+3,n)*sin(a(m+3)-a(n)-B(m+3,n)*cos(a(m+3)-a(n); end L(m,m)=U(m+3)*U(m+3)*B(m+3,m+3)-U(m+3)*sum(Ai); end;%雅可比矩阵i等于j%Ai(n)=U(n)*(G(m+1,n)*sin(a(m+1)-a(n)-B(m+1,n)*c
22、os(a(m+1)-a(n);%Bi(n)=U(n)*(G(m+1,n)*cos(a(m+1)-a(n)+B(m+1,n)*sin(a(m+1)-a(n); for m=1:8 for n=1:8 JJ(m,n)=H(m,n); endendfor m=1:8 for n=1:6 JJ(m,n+8)=N(m,n); endendfor m=1:6 for n=1:8 JJ(m+8,n)=M(m,n); endendfor m=1:6 for n=1:6 JJ(m+8,n+8)=L(m,n); endend%雅可比矩阵JJ形成for m=1:8 PQ(m)=dP(m+1);endfor m=9:
23、14 PQ(m)=dQ(m-5);enddUa=-inv(JJ)*PQ;precison=max(abs(dUa);for m=1:8 a(m+1)=a(m+1)+dUa(m);endfor m=9:14 U(m-5)=U(m-5)+dUa(m);endk=k+1;endu=U.*cos(a)+i*U.*sin(a);i=y*u.;I=conj(i);pq=u.*I.;P=real(pq);Q=imag(pq);format longXE=1 0.3 1.137; 2 0.3 1.211; 3 0.3 1.043; 4 0 0; 5 0 0; 6 0 0; 7 0 0; 8 0 0; 9 0
24、0;disp(大作业)disp(迭代次数)kdisp(节点导纳矩阵)ydisp(电压幅值)Udisp(电压相角)a=a*360/(2*pi)disp(节点有功)Pdisp(节点无功)Qdisp(输出修正后的发电机节点和负荷节点的自导纳)for m=1:9 if abs(P(m)0. & abs(Q(m)0. if P(m)0 y(m,m)=y(m,m)-1i*(1/XE(m,2); else y(m,m)=y(m,m)-conj(pq(m)/U(m)/U(m); end disp(m); disp( 点的自导纳=); disp(y(m,m); endendz=inv(y)line=z(:,4)
25、;If=u(4)/z(4,4); for m=1:9 uu(m)=u(m)-line(m)*If; end If1=1/z(4,4); for m=1:9 uu1(m)=1-line(m)*z(4,4); end wuchai0=(abs(If)-abs(If1)/abs(If); wuchai1=(angle(If)-angle(If1)/angle(If); for m=1:9 wuchau0(m)=(abs(uu(m)-abs(uu1(m)/abs(uu(m); end for m=1:9 wuchau1(m)=(angle(uu(m)-angle(uu1(m)/angle(uu(m);
26、 end for m=1:9 hang=A(m,:); y=zeros(2,2); if hang(1)3 & hang(2)3 y(1,1)=y(1,1)+1/(A(m,3)+1i*A(m,4)+1i*A(m,5); y(1,2)=y(1,2)-1/(A(m,3)+1i*A(m,4); y(2,1)=y(2,1)-1/(A(m,3)+1i*A(m,4); y(2,2)=y(2,2)+1/(A(m,3)+1i*A(m,4)+1i*A(m,5); ii(m,1)=y(1,1)*uu(hang(1)+y(1,2)*uu(hang(2); ii(m,2)=y(2,1)*uu(hang(1)+y(2,
27、2)*uu(hang(2); else y(1,1)=y(1,1)+1/(A(m,3)+1i*A(m,4); y(1,2)=y(1,2)-1/(A(m,3)+1i*A(m,4); y(2,1)=y(2,1)-1/(A(m,3)+1i*A(m,4); y(2,21)=y(2,2)+1/(A(m,3)+1i*A(m,4); ii(m,1)=(uu(hang(2)-uu(hang(1)*y(1,2); ii(m,2)=(uu(hang(1)-uu(hang(2)*y(2,1); end end a=zeros(9,6); for m=1:9 a(m,1)=A(m,1); a(m,2)=A(m,2);
28、 a(m,3)=real(ii(m,1); a(m,4)=imag(ii(m,1); a(m,5)=real(ii(m,2); a(m,6)=imag(ii(m,2); end b=zeros(9,6); for m=1:9 b(m,1)=A(m,1); b(m,2)=A(m,2); b(m,3)=abs(ii(m,1); b(m,4)=angle(ii(m,1)*360/(2*pi); b(m,5)=abs(ii(m,2); b(m,6)=angle(ii(m,2)*360/(2*pi); end for n=1:3 for m=1:5 s=a(m,:); a(m,:)=a(m+1,:);
29、a(m+1,:)=s; m=m+1; end n=n+1; end for n=1:3 for m=1:5 s=a(m+3,:); a(m+3,:)=a(m+4,:); a(m+4,:)=s; m=m+1; end n=n+1; end for n=1:3 for m=1:5 s=b(m,:); b(m,:)=b(m+1,:); b(m+1,:)=s; m=m+1; end n=n+1; end for n=1:3 for m=1:5 s=b(m+3,:); b(m+3,:)=b(m+4,:); b(m+4,:)=s; m=m+1; end n=n+1; end disp(输出节点阻抗矩阵的第
30、四列)linedisp(精确计算短路电流If的模值和角度)disp( 模值)abs(If)disp( 相角)angle(If)*360/(2*pi)disp(精确计算短路各节点电压模值)sb=abs(uu);sb.disp(近似计算短路电流If的模值和角度)disp( 模值)abs(If1)disp( 相角)angle(If1)*360/(2*pi)disp(近似计算短路各节点电压模值)sb=abs(uu1);sb.disp(电流幅值误差)wuchai0disp(电流相角误差)wuchai1disp(电压幅值误差)sb2=wuchau0(m);sb.disp(电压相角误差)sb2=wuchau1(m);sb.disp(短路时各支路电流)disp( i j Iij实部 Iij虚部 Iji实部 Iji虚部)adisp( i j Iij模值 Iij相角 Iji模值 Iji相角)b专心-专注-专业