潮流计算的计算机算法
高等电力系统分析
（潮流计算的计算机算法）PQ分解法潮流计算（IEEE14)
目录
一、MATLAB源程序
二、对支路参数（B1）、节点参数（B2）的说明
三、带入数据，运行结果
一、MATLAB源程序
clear
close all
n=input('请输入节点数:n=');
n1=input('请输入支路数:n1=');
isb=input('请输入平衡节点号:isb=');
pr=input('请输入误差精度:pr=');
B1=input('请输入支路参数:B1=');
B2=input('请输入节点参数:B2=');
n2=input('请输入PQ节点个数:n2=');
Y=zeros(n);
for i=1:n1
p=B1(i,1);
q=B1(i,2);
Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4)*1j); %非对角元
Y(q,p)=Y(p,q);
Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j; %对角元 Y(q,q)=Y(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;
end
disp('导纳矩阵Y=');
disp(Y)
%---------------------------------------------
%---------------下面是求P,Q,V,O矩阵---------------
V=zeros(1,n);O=zeros(1,n);P=zeros(1,n);Q=zeros(1,n);
G=real(Y);B=imag(Y);
for i=1:n
P(i)=B2(i,3);
Q(i)=B2(i,4);
V(i)=B2(i,5);
O(i)=B2(i,6);
end
B3=B(1:n-1,1:n-1); %不含平衡节点，由节点导纳虚部构成
B4=B(1:n2,1:n2); %所有PQ节点
%----------------------------------------------
%---------------下面是求ΔP,ΔQ矩阵---------------
DX=0;ICT=1;Mp=1;Mq=1;
while ICT~=0
m1=1;m2=1;
for i=1:n
if i~=isb
C(i)=0;
D(i)=0;
for j1=1:n
C(i)=C(i)+V(i)*V(j1)*(G(i,j1)*cos(O(i)-O(j1))+B(i,j1)*sin(O(i)-O( j1)));
D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*sin(O(i)-O(j1))-B(i,j1)*cos(O(i)-O( j1)));
end
DP(m1)=P(i)-C(i);
m1=m1+1;
if B2(i,2)==1
DQ(m2)=Q(i)-D(i);
m2=m2+1;
end
end
end
m1=m1-1; %所有节点数
m2=m2-1; %PQ节点数
DPQ=[DP';DQ']; %求DP,DQ
V1=V(:,1:m1);
V2=diag(V1);
V3=inv(V2); %对V矩阵求逆
H=V3*DP'; %ΔP/V
K=-inv(B3)*H; %-ΔP/V/B3
deltO=V3*K; %Δ角=-ΔP/V/V/B3
max1=max(abs(DP));
for i=1:m1
if max1 Mp=0; else O(i)=O(i)+deltO(i)'; Mq=1; end end V4=V(:,1:m2); V5=diag(V4); V6=inv(V5); L=V6*DQ'; N=-inv(B4)*L; deltV=N; %ΔV=-ΔQ/V/B max2=max(abs(DQ)); for i=1:m2 if max2 Mq=0; else if B2(i,2)==1; V(i)=V(i)+deltV(i)'; Mp=1; end end end if Mp==0&&Mq==0 ICT=0; else ICT=1; end DX=DX+1; end %------------------------------------------- %----------------迭代结束，开始输出结果---------------- disp('----------------------------------------------'); disp('迭代次数为：'); disp(DX); for i=1:n E(i)=V(i)*cos(O(i))+1j*V(i)*sin(O(i)); o(i)= 180*angle(E(i))/pi; end disp('---------------------------------------------'); disp('修正后各节点电压标么值为（节点号从小到大排列）:'); disp(V); disp('---------------------------------------------'); disp('修正后各节点电压相角为（节点号从小到大排列）:'); disp(o); %------------计算各个节点的功率---------------------- disp('----------------------------------------------'); disp('各节点的功率为:'); for p=1:n C(p)=0; for q=1:n C(p)=C(p)+conj(Y(p,q)*conj(E(q))); end S(p)=E(p)*C(p); end disp(S); %------------计算各支路的功率---------------------- for i=1:n1 p=B1(i,1);q=B1(i,2); Si(p,q)=E(p)*(conj(E(p))*conj(Y(p,p)-Y(p,q))+(conj(E(p))-conj(E(q )))*conj(Y(p,q))); disp('----------------------------------------------'); disp('各条支路的首端功率为:'); disp(Si(p,q)); Si(q,p)=E(q)*(conj(E(q))*conj(Y(q,q)-Y(p,q))+(conj(E(q))-conj(E(p )))*conj(Y(p,q))); disp('----------------------------------------------'); disp('各条支路的末端功率为:'); disp(Si(q,p)); DS(i)=Si(p,q)+Si(q,p); disp('----------------------------------------------'); disp('各条支路的功率损耗为:'); disp(DS(i)); end %----------------计算平衡节点功率------------- Sp=0; for i=1:n Sp=Sp+V(n)*conj(Y(n,i))*conj(V(i)); end disp('----------------------------------------------'); disp('平衡节点功率为:'); disp(Sp); 二、对支路参数（B1）、节点参数（B2）的说明 1.节点数：14 2.支路数：20 3.支路矩阵B1的各支路参数： 起点编号，终点编号，电阻，电抗，电导，电纳 [1 2 0.01335 0.04211 0 0 ; 1 3 0 0.2091 2 0 0 ; 1 4 0 0.55618 0 0 ; 1 10 0.05811 0.1763 2 0 0.034 ; 1 11 0.06701 0.17103 0 0.0128 ; 2 10 0.05695 0.17388 0 0.0346 ; 2 12 0 0.25202 0 0 ; 2 14 0.0540 3 0.2230 4 0 0.0492 ; 3 4 0 0.11001 0 0 ; 3 13 0 0.17615 0 0 ; 4 5 0.03181 0.0845 0 0 ; 4 9 0.12711 0.27038 0 0 ; 5 6 0.08205 0.19207 0 0 ; 6 12 0.09498 0.1989 0 0 ; 7 8 0.22092 0.19988 0 0 ; 7 12 0.12291 0.25581 0 0 ; 8 9 0.17093 0.34802 0 0 ; 8 12 0.06615 0.13027 0 0 ; 10 11 0.04699 0.19797 0 0.0438 ; 10 14 0.01938 0.05917 0 0.0528 ;] 4. 节点参数矩阵B2的各节点参数：（对应的每一列为） 节点编号，类型，注入有功，注入无功，电压幅值，电压相位其中节点类型：1=PQ节点，2=PV节点，0=平衡节点 [1 1 -0.478 0.039 1 0; 2 1 -0.076 -0.016 1 0; 3 1 0 0 1 0; 4 1 -0.29 5 -0.16 6 1 0; 5 1 -0.09 -0.058 1 0; 6 1 -0.035 -0.018 1 0; 7 1 -0.061 -0.016 1 0; 8 1 -0.135 -0.058 1 0; 9 1 -0.149 -0.05 1 0; 10 2 0.183 0 1.045 0; 11 2 -0.942 0 1.01 0; 12 2 -0.112 0.047 1.7 0; 13 2 0 0.174 1.9 0; 14 0 0 0 1.06 0;] 三、带入数据，运行结果 >> clear close all n=input('请输入节点数:n='); n1=input('请输入支路数:n1='); isb=input('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入支路参数:B1='); B2=input('请输入节点参数:B2='); n2=input('请输入PQ节点个数:n2='); Y=zeros(n); for i=1:n1 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4)*1j); %非对角元 Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j; %对角元 Y(q,q)=Y(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j; end disp('导纳矩阵Y='); disp(Y) %--------------------------------------------- %---------------下面是求P,Q,V,O矩阵--------------- V=zeros(1,n);O=zeros(1,n);P=zeros(1,n);Q=zeros(1,n); G=real(Y);B=imag(Y); for i=1:n P(i)=B2(i,3); Q(i)=B2(i,4); V(i)=B2(i,5); O(i)=B2(i,6); end B3=B(1:n-1,1:n-1); %不含平衡节点，由节点导纳虚部构成 B4=B(1:n2,1:n2); %所有PQ节点 %---------------------------------------------- %---------------下面是求ΔP,ΔQ矩阵--------------- DX=0;ICT=1;Mp=1;Mq=1; while ICT~=0 m1=1;m2=1; for i=1:n if i~=isb C(i)=0; D(i)=0; for j1=1:n C(i)=C(i)+V(i)*V(j1)*(G(i,j1)*cos(O(i)-O(j1))+B(i,j1)*sin(O(i)-O(j1))); D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*sin(O(i)-O(j1))-B(i,j1)*cos(O(i)-O(j1))); end DP(m1)=P(i)-C(i); m1=m1+1; if B2(i,2)==1 DQ(m2)=Q(i)-D(i); m2=m2+1; end end end m1=m1-1; %所有节点数 m2=m2-1; %PQ节点数 DPQ=[DP';DQ']; %求DP,DQ V1=V(:,1:m1); V2=diag(V1); V3=inv(V2); %对V矩阵求逆 H=V3*DP'; %ΔP/V K=-inv(B3)*H; %-ΔP/V/B3 deltO=V3*K; %Δ角=-ΔP/V/V/B3 max1=max(abs(DP)); for i=1:m1 if max1 Mp=0; else O(i)=O(i)+deltO(i)'; Mq=1; end end V4=V(:,1:m2); V5=diag(V4); V6=inv(V5); L=V6*DQ'; N=-inv(B4)*L; deltV=N; %ΔV=-ΔQ/V/B max2=max(abs(DQ)); for i=1:m2 if max2 Mq=0; else if B2(i,2)==1; V(i)=V(i)+deltV(i)'; Mp=1; end end end if Mp==0&&Mq==0 ICT=0; else ICT=1; end DX=DX+1; end %------------------------------------------- %----------------迭代结束，开始输出结果---------------- disp('----------------------------------------------'); disp('迭代次数为：'); disp(DX); for i=1:n E(i)=V(i)*cos(O(i))+1j*V(i)*sin(O(i)); o(i)= 180*angle(E(i))/pi; end disp('---------------------------------------------'); disp('修正后各节点电压标么值为（节点号从小到大排列）:'); disp(V); disp('---------------------------------------------'); disp('修正后各节点电压相角为（节点号从小到大排列）:'); disp(o); %------------计算各个节点的功率---------------------- disp('----------------------------------------------'); disp('各节点的功率为:'); for p=1:n C(p)=0; for q=1:n C(p)=C(p)+conj(Y(p,q)*conj(E(q))); end S(p)=E(p)*C(p); end disp(S); %------------计算各支路的功率---------------------- for i=1:n1 p=B1(i,1);q=B1(i,2); Si(p,q)=E(p)*(conj(E(p))*conj(Y(p,p)-Y(p,q))+(conj(E(p))-conj(E(q)))*conj(Y(p,q))); disp('----------------------------------------------'); disp('各条支路的首端功率为:'); disp(Si(p,q)); Si(q,p)=E(q)*(conj(E(q))*conj(Y(q,q)-Y(p,q))+(conj(E(q))-conj(E(p)))*conj(Y(p,q))); disp('----------------------------------------------'); disp('各条支路的末端功率为:'); disp(Si(q,p)); DS(i)=Si(p,q)+Si(q,p); disp('----------------------------------------------'); disp('各条支路的功率损耗为:'); disp(DS(i)); end %----------------计算平衡节点功率------------- Sp=0; for i=1:n Sp=Sp+V(n)*conj(Y(n,i))*conj(V(i)); end disp('----------------------------------------------'); disp('平衡节点功率为:'); disp(Sp); 请输入节点数:n=14 请输入支路数:n1=20 请输入平衡节点号:isb=14 请输入误差精度:pr=0.00001 请输入支路参数:B1=[1 2 0.01335 0.04211 0 0 ; 1 3 0 0.2091 2 0 0 ; 1 4 0 0.55618 0 0 ; 1 10 0.05811 0.1763 2 0 0.034 ; 1 11 0.06701 0.17103 0 0.0128 ; 2 10 0.05695 0.17388 0 0.0346 ; 2 12 0 0.25202 0 0 ; 2 14 0.0540 3 0.2230 4 0 0.0492 ; 3 4 0 0.11001 0 0 ; 3 13 0 0.17615 0 0 ; 4 5 0.03181 0.0845 0 0 ; 4 9 0.12711 0.27038 0 0 ; 5 6 0.08205 0.19207 0 0 ; 6 12 0.09498 0.1989 0 0 ; 7 8 0.22092 0.19988 0 0 ; 7 12 0.12291 0.25581 0 0 ; 8 9 0.17093 0.34802 0 0 ; 8 12 0.06615 0.13027 0 0 ; 10 11 0.04699 0.19797 0 0.0438 ; 10 14 0.01938 0.05917 0 0.0528 ;] 请输入节点参数:B2=[1 1 -0.478 0.039 1 0; 2 1 -0.076 -0.016 1 0; 3 1 0 0 1 0; 4 1 -0.29 5 -0.16 6 1 0; 5 1 -0.09 -0.058 1 0; 6 1 -0.035 -0.018 1 0; 7 1 -0.061 -0.016 1 0; 8 1 -0.135 -0.058 1 0; 9 1 -0.149 -0.05 1 0; 10 2 0.183 0 1.045 0; 11 2 -0.942 0 1.01 0; 12 2 -0.112 0.047 1.7 0; 13 2 0 0.174 1.9 0; 14 0 0 0 1.06 0;] 请输入PQ节点个数:n2=9 导纳矩阵Y= Columns 1 through 5 10.5130 -38.2963i -6.8410 +21.5786i 0.0000 + 4.7819i 0.0000 + 1.7980i 0.0000 + 0.0000i -6.8410 +21.5786i 9.5680 -34.8916i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 4.7819i 0.0000 + 0.0000i 0.0000 -19.5490i 0.0000 + 9.0901i 0.0000 + 0.0000i 0.0000 + 1.7980i 0.0000 + 0.0000i 0.0000 + 9.0901i 5.3261 -24.2825i -3.9020 +10.3654i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -3.9020 +10.3654i 5.7829 -14.7683i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.8809 + 4.4029i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.4240 + 3.0291i 0.0000 + 0.0000i -1.6860 + 5.1158i -1.7011 + 5.1939i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.9860 + 5.0688i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 3.9679i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 5.6770i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.0259 + 4.2350i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i Columns 6 through 10 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.6860 + 5.1158i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.7011 + 5.1939i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.4240 + 3.0291i 0.0000 + 0.0000i -1.8809 + 4.4029i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 3.8359 - 8.4970i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 4.0150 - 5.4279i -2.4890 + 2.2520i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -2.4890 + 2.2520i 6.7249 -10.6697i -1.1370 + 2.3150i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.1370 + 2.3150i 2.5610 - 5.3440i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 9.5213 -30.1895i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.1350 + 4.7819i -1.9550 + 4.0941i -1.5260 + 3.1760i -3.0989 + 6.1028i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -4.9991 +15.2631i Columns 11 through 14 -1.9860 + 5.0688i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 3.9679i 0.0000 + 0.0000i -1.0259 + 4.2350i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 5.6770i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.9550 + 4.0941i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.5260 + 3.1760i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -3.0989 + 6.1028i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.1350 + 4.7819i 0.0000 + 0.0000i 0.0000 + 0.0000i -4.9991 +15.2631i 3.1210 - 9.7941i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 6.5799 -17.3407i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 - 5.6770i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 6.0250 -19.3961i ---------------------------------------------- 迭代次数为： 42 --------------------------------------------- 修正后各节点电压标么值为（节点号从小到大排列）: Columns 1 through 10 1.2128 1.2148 1.5627 1.5364 1.5602 1.6264 1.6792 1.6654 1.5809 1.0450 Columns 11 through 14 1.0100 1.7000 1.9000 1.0600 --------------------------------------------- 修正后各节点电压相角为（节点号从小到大排列）: Columns 1 through 10 -13.1377 -11.9998 -14.4540 -15.0006 -15.3011 -15.7590 -16.4954 -16.3362 -16.0220 -5.2081 Columns 11 through 14 -12.1568 -16.1918 -14.4540 0 ---------------------------------------------- 各节点的功率为: Columns 1 through 5 1.4718 - 1.8800i 0.9041 - 1.1423i 0.0000 - 0.0000i -0.1314 + 0.1399i -0.0475 - 0.0646i Columns 6 through 10 -0.0300 + 0.0035i 0.0039 - 0.1228i -0.0656 - 0.0724i -0.0352 - 0.1437i -1.5950 - 1.3905i Columns 11 through 14 -0.4870 - 1.2396i 3.2005 + 3.0414i 1.7586 + 3.1847i -2.7057 + 0.3992i ---------------------------------------------- 各条支路的首端功率为: 26.1711 +87.9128i ---------------------------------------------- 各条支路的末端功率为: 23.5667 +83.4812i ---------------------------------------------- 各条支路的功率损耗为: 4.9738e+01 + 1.7139e+02i ---------------------------------------------- 各条支路的首端功率为: 15.2545 +65.3872i ---------------------------------------------- 各条支路的末端功率为: 0.2082 +56.8008i ---------------------------------------------- 各条支路的功率损耗为: 1.5463e+01 + 1.2219e+02i ---------------------------------------------- 各条支路的首端功率为: 15.3538 +59.6752i ---------------------------------------------- 各条支路的末端功率为: 12.6813 +60.6685i ---------------------------------------------- 各条支路的功率损耗为: 2.8035e+01 + 1.2034e+02i ---------------------------------------------- 各条支路的首端功率为: 18.4735 +62.4536i ---------------------------------------------- 各条支路的末端功率为: 11.6194 +39.6840i ---------------------------------------------- 各条支路的功率损耗为: 3.0093e+01 + 1.0214e+02i ---------------------------------------------- 各条支路的首端功率为: 18.0012 +62.4930i ---------------------------------------------- 各条支路的末端功率为: 5.5097 +1 6.2405i ---------------------------------------------- 各条支路的功率损耗为: 23.5109 +78.7335i ---------------------------------------------- 各条支路的首端功率为: 17.0452 +57.7870i ---------------------------------------------- 各条支路的末端功率为: 11.7622 +39.7706i ---------------------------------------------- 各条支路的功率损耗为: 28.8075 +97.5576i ---------------------------------------------- 各条支路的首端功率为: 13.5220 +59.6677i ---------------------------------------------- 各条支路的末端功率为: 19.6150 +58.2875i ---------------------------------------------- 各条支路的功率损耗为: 3.3137e+01 + 1.1796e+02i ---------------------------------------------- 各条支路的首端功率为: 16.5471 +56.5546i ---------------------------------------------- 各条支路的末端功率为: 6.9281 +2 7.4025i ---------------------------------------------- 各条支路的功率损耗为: 23.4752 +83.9571i ---------------------------------------------- 各条支路的首端功率为: -0.2082 +69.5644i ---------------------------------------------- 各条支路的末端功率为: 12.7806 +79.1441i ---------------------------------------------- 各条支路的功率损耗为: 1.2572e+01 + 1.4871e+02i ---------------------------------------------- 各条支路的首端功率为: 0.0000 +64.5962i ---------------------------------------------- 各条支路的末端功率为: -0.0000 +37.3498i ---------------------------------------------- 各条支路的功率损耗为: -1.7764e-15 + 1.0195e+02i ---------------------------------------------- 各条支路的首端功率为: 21.7957 +82.2160i ---------------------------------------------- 各条支路的末端功率为: 23.5612 +60.7482i ---------------------------------------------- 各条支路的功率损耗为: 4.5357e+01 + 1.4296e+02i ---------------------------------------------- 各条支路的首端功率为: 15.8994 +64.7377i ---------------------------------------------- 各条支路的末端功率为: 9.9896 +20.6498i ---------------------------------------------- 各条支路的功率损耗为: 25.8891 +85.3875i ---------------------------------------------- 各条支路的首端功率为: 18.7607 +47.1609i ---------------------------------------------- 各条支路的末端功率为: 15.0088 +33.6107i ---------------------------------------------- 各条支路的功率损耗为: 33.7695 +80.7716i ---------------------------------------------- 各条支路的首端功率为: 15.4666 +33.8366i ---------------------------------------------- 各条支路的末端功率为: 24.5068 +61.3933i ---------------------------------------------- 各条支路的功率损耗为: 39.9734 +95.2299i ---------------------------------------------- 各条支路的首端功率为: 18.2997 +21.5842i ---------------------------------------------- 各条支路的末端功率为: 25.5956 +35.9107i ---------------------------------------------- 各条支路的功率损耗为: 43.8954 +57.4949i ---------------------------------------------- 各条支路的首端功率为: 15.7256 +24.3488i ---------------------------------------------- 各条支路的末端功率为: 23.3240 +59.2040i ---------------------------------------------- 各条支路的功率损耗为: 39.0496 +83.5528i ---------------------------------------------- 各条支路的首端功率为: 21.6792 +35.6719i ---------------------------------------------- 各条支路的末端功率为: 9.3603 +19.4667i ---------------------------------------------- 各条支路的功率损耗为: 31.0395 +55.1385i 第四章潮流计算的计算机算法 第一节概述 潮流计算是电力系统最基本、最常用的计算。根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值及相角），各元件中流过的功率、整个系统的功率损耗等。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。因此潮流计算在电力系统的规划设计、生产运行、调度管理及科学研究中都有着广泛的应用。 电力系统潮流计算分为离线潮流计算和在线潮流计算。前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。本章主要讨论离线潮流计算问题，它的基本算法同样适用于在线潮流计算。 潮流计算在数学上是多元非线性方程组的求解问题，求解的方法有很多种。自从五十年代计算机应用于电力系统以来，当时求解潮流的方法是以节点导纳矩阵为基础的逐次代入法（导纳法），后来为解决导纳法的收敛性较差的问题，出现了以阻抗矩阵为基础的逐次代入法（阻抗法）。到六十年代，针对阻抗法占用计算机内存大的问题又出现了分块阻抗法及牛顿－拉夫逊（Newton－Raphson）法。Newton —Raphson法是数学上解非线形方程式的有效方法，有较好的收敛性。将N－R法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧，使N－R法在收敛性、占用内存、计算速度方面的优点都超过了阻抗法，成为六十年代末期以后普遍采用的方法。同时国内外广泛研究了诸如非线形规划法、直流法、交流法等各种不同的潮流计算方法。七十年代以来，又涌现出了更新的潮流计算方法。其中有1974年由B、Stott、O、Alsac 提出的快速分解法以及1978年由岩本伸一等提出的保留非线性的高 129 关于电力系统经济调度的潮流计算分析 发表时间：2016-05-24T15:57:29.347Z 来源：《电力设备》2016年第2期作者：秦先威 [导读] （国网山东省电力公司烟台市牟平区供电公司山东烟台 264100）随着经济的快速发展和科技的不断进步，社会各行业对电力资源的需求量越来越大，我国的电力系统建设规模也越来越大。 （国网山东省电力公司烟台市牟平区供电公司山东烟台 264100） 摘要:潮流计算是电力调度中最重要也是最基本的计算之一,它应用于电力系统中实时电价计算、输电权分配、网络阻塞管理等多方面。 关键词:电力系统；经济调度；潮流计算 前言 随着经济的快速发展和科技的不断进步，社会各行业对电力资源的需求量越来越大，我国的电力系统建设规模也越来越大。电力调度对电力系统的正常运行有很大的影响，而潮流计算则是电力调度中最重要的基本计算方法，潮流计算对电价计算、输电分配、电网线路管理有十分重要的影响。随着经济的快速发展，我国的电力企业得到了飞速的发展，与此同时，人们对供电质量的要求也越来越高，为满足人们的用电需求，电力系统在运行过程中，必须保证电力调度的合理性、科学性，潮流计算是电力系统经济调度最重要的计算方法之一，潮流计算的结果准确性很高，科学性很强，潮流计算对电力系统经济调度有十分重要的作用。 一、潮流计算的概述 1.1 潮流计算的概述 潮流计算是指利用已知的电网接线方式、参数、运行条件，将电力系统的各个母线电压、支路电流、功率、网损计算出来。通过潮流计算能判断出正在运行的电力系统的母线电压、支路电流、功率是否在允许范围内运行，如果超出允许范围，就需要采用合理的措施，对电力系统的进行方式进行调整。在电力系统规划过程中，采用潮流计算，能为电网供电方案、电气设备的选择提供科学的依据，同时潮流计算还能为自动装置定整计算、继电保护、电力系统稳定计算、故障计算提供原始数据。 1.2 潮流计算的电气量 潮流计算是根据电力系统接线方式、运行条件、参数等已知条件，将稳定状态下电力系统的电气量计算出来。一般情况下，给出的条件有电源、负荷节点的功率、平衡节点的电压、相位角、枢纽点的电压，需要计算的电气量有各节点的电压、相位角、各支路通过的电流、功率、网络的功率损耗等。 1.3 传统的潮流计算方法 传统的潮流计算方法，包括很多不同的内容，具有一定的优点和缺点。例如，传统的潮流计算方法，包括非线性规划法、二次规划法和线性规划法等。在电力系统经济调度的过程中，应用传统的潮流计算方法，优点是：可以根据目标函数的导数信息，确定需要进行搜索的方向，因此在计算的时候，具有较快的速度和清晰的计算过程。而且，可信度比较高。 1.5 智能的潮流计算方法 潮流计算中人工智能方法的优点是：随机性：属于全局优化算法，跳出局部极值点比较容易；与导数无关性：在工程中，一些优化问题的目标函数处于不可导状态。如果进行近似和假设，会对求解的真实性造成影响；内在并行性：操作对象为一组可行解，在一定程度上可以克服内在并发性开放中性能的不足。而其缺点，主要是：需要按照概率进行操作，不能保证可以完全获取最优解；算法中的一些控制参数需要根据经验人文地给出，对专家经验和一定量的试验要求比较高；表现不稳定，在同一问题的不同实例中应用算法会出现不同的效果。 二、潮流计算的分类 根据电力系统的运行状态，潮流计算可以分为离线计算和在线计算两种方法，离线计算主要用于电力系统规划设计和电力系统运行方式安排中；在线计算主要用于电力系统运行监控和控制中；根据潮流计算的发展，潮流计算可以分为传统方法和人工智能方法两种情况，下面分别对这两种方法进行分析。 2.1 潮流计算的传统方法 潮流计算的传统方法有非线性规划法、线性规划法、二次规划法等几种情况，潮流计算的传统方法具有计算速度快、解析过程清晰、结果真实可靠等优点，但传统方法对目标函数有一定的限制，需要简化处理，这样求出来的值有可能不是最优值。 2.2 潮流计算的人工智能方法 潮流计算的人工智能方法是一种新兴的方法，人工智能方法不会过于依赖精确的数学模型，它有粒子群优化算法、遗传法、模拟退火法等几种情况，人工智能方法的计算结果和导数没有关系，其操作对象是一组可行解，能克服内在并行性存在的问题，但人工智能方法表现不太稳定，在计算过程中，有的控制参数需要根据经验得出，因此，采用人工智能方法进行计算时，需要计算人员有丰富的经验。 三、潮流计算在电力系统经济调度中的应用 3.1 在输电线路线损计算的应用 在进行输电线路线损计算过程中，通过潮流计算能得出经济潮流数据。潮流程度能根据线路的功率因数、有功负荷、无功负荷等参数，计算出潮流线损，例如一条长为38.1km，型号为LGJ—150的导线，当潮流为20MW、功率因数为0.9时，该线路线损为0.24MW，线损率为1.18%；当潮流为30MW、功率因数为0.9时，该线路线损为0.57MW，线损率为1.91%；潮流为50MW、功率因数为0.9时，该线路线损为1.95MW，线损率为3.90%；由此可以看出，潮流小于30MW时，线损率小于2%，潮流超过50MW时，线损率将超过4%，因此，该输电线路的经济输送潮流为30MW以下。调度人员可以根据计算结果，编制线路经济运行方案，从而实现节能调度。 3.2 在变压器变损中的应用 调度人员可以利用潮流计算程序，将变压器在不同负荷下的损耗、变损率计算出来，从而为变压器控制提供依据。例如一台40MVA双 潮流计算的基本算法及使用方法 一、 潮流计算的基本算法 1. 牛顿－拉夫逊法 1．1 概述 牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把对非线 性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。 牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏 导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。 1．2 一般概念 对于非线性代数方程组 ()0=x f 即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1) 在待求量x 的某一个初始计算值() 0x 附件，将上式展开泰勒级数并略去二阶及以上的高 阶项，得到如下的线性化的方程组 ()()()() ()0000=?'+x x f x f (1－2) 上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 ()() ()[]()()0 1 00x f x f x -'-=? (1－3) 将() 0x ?和() 0x 相加，得到变量的第一次改进值()1x 。接着再从() 1x 出发，重复上述计算 过程。因此从一定的初值() 0x 出发，应用牛顿法求解的迭代格式为 ()()()()() k k k x f x x f -=?' (1－4) ()()()k k k x x x ?+=+1 (1－5) 上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代 1.节能效果好，环境负面影响小 2.可以提高供电的安全性和可靠性 3.调峰性能好 4.服务多样化 5.可以满足特殊场合的需才 6.节省投资，经济效益好 在传统集中式供电模式中，功率从高电压等级的电源侧传送到低电压等级的用户侧， 配电系统是单端供电。分布式发电的引入，将改变传统配电网络拓扑，使得潮流特性发生变化，此外接入电网的电气元件和用户侧的电压也将产生变化，给系统的稳定性可能造成很大影响。分布式发电并网会产生两个方面的问题:一是并网系统本身的结构和性能;另 一个是分布式发电并网后对电力系统运行、控制、保护等各方面产生的影响。 <1)对配电系统负荷预测和规划的影响。在传统配电系统引入分布式发电加大了电 力系统负荷预测的不确定性，使配电系统规划者难于准确预测负荷增长情况，从而影响配电系统的规划设计。另外，在配电系统引入分布式发电，对降低网损起到很大作用。但是配电系统本身节点数非常多，分布式发电节点的出现，使得在所有可能网络结构中寻找分布式发电的最优网络布置方案更加困难。 (2)对继电保护的影响。分布式发电接入配电网后短路电流将会增大，加之整个配 电网变成多电源的网络，网络潮流的流向具有不确定性，而传统系统的潮流是从电源到拜户单向流动，故保护系统的设计基础应该发生相应的变化。 C3)对系统可靠性的影响。若分布式发电仅作为备用电源，则可提高系统供电的可 靠性;若分布式发电与电网并联运行，就有可能降低系统的可靠性。分布式发电接入配网，改变了电网拓扑结构，使得短路电流发生变化，导致按原有设计的保护装置误动作，破坏保护设备间的协调运行，妨碍了自动重合闸动作。 t4)对系统电能质量的影响。比如引起电压闪烁，还可能引入谐波，造成谐波污染; (5)对系统稳定性的影响。配电系统中引入少量的分布式发电对整个电网不会构成 太大的影响，但是当电网中存在较多的分布式发电单元或者存在大容量的分布式发电单元时，分布式发电的引入将会带来系统稳定问题。 t6)对电力市场走向和最后格局将产生深远影响。电力公司和用户间将形成新型关系，用户不仅可以从电力公司买电，也可用自己的分布式发电向其卖电或为其提供有偿削峰、紧急功率支持等服务，分布式发电也为其它行业(如天然气公司)进入市场打开了方 便之门，因此未来电力市场的竞争将更加激烈。 近年来，随着分布式发电技术的兴起，有许多小型的水电站、热电联产电厂、风电场 直接连接到配电系统，辐射式的网络将变为一遍布电源和用户互联的网络，潮流也不一定单向地从变电站母线流向各负荷，有可能会出现回流和复杂的电压变化，同时可能增大或减小系统损耗，这取决于分布式发电的位置、与负荷量的相对大小以及网络的拓扑结构等因素y o}，给潮流分析带来了影响。另外，如果配电网中含有风力发电或太阳能光伏发电系统，由于它们的输出受天气的影响很大，具有随机变化的特性。上述负荷的变化、风电和光伏电池输出的变化可能会在配电网中交替出现，使系统的潮流具有随机性分布式发电的定义 目前，国际上对于分布式发电并没有一个统一的规范化定义，但基本概念是一致的。 相对于传统集中式发电而言，分布式发电是指分布安置在需求侧的能源梯级综合利用，其 容量很小(大多在几十KW一几十MW之间)。通过在需求现场根据用户对能源的不同 高等电力系统分析 （潮流计算的计算机算法）PQ分解法潮流计算（IEEE14) 目录 一、MATLAB源程序 二、对支路参数（B1）、节点参数（B2）的说明 三、带入数据，运行结果 一、MATLAB源程序 clear close all n=input('请输入节点数:n='); n1=input('请输入支路数:n1='); isb=input('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入支路参数:B1='); B2=input('请输入节点参数:B2='); n2=input('请输入PQ节点个数:n2='); Y=zeros(n); for i=1:n1 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4)*1j); %非对角元 Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j; %对角元 Y(q,q)=Y(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j; end disp('导纳矩阵Y='); disp(Y) %--------------------------------------------- %---------------下面是求P,Q,V,O矩阵--------------- V=zeros(1,n);O=zeros(1,n);P=zeros(1,n);Q=zeros(1,n); G=real(Y);B=imag(Y); for i=1:n P(i)=B2(i,3); Q(i)=B2(i,4); V(i)=B2(i,5); O(i)=B2(i,6); end B3=B(1:n-1,1:n-1); %不含平衡节点，由节点导纳虚部构成 B4=B(1:n2,1:n2); %所有PQ节点 %---------------------------------------------- %---------------下面是求ΔP,ΔQ矩阵--------------- DX=0;ICT=1;Mp=1;Mq=1; while ICT~=0 m1=1;m2=1; for i=1:n 随机潮流研究现状 摘要：随机潮流可用于分析线路潮流、节点电压的概率分布、期望值、方差和极限值，以期对整个电网在各种运行条件下的性能有一个全面、综合的评价，并对电网存在的薄弱环节做出量化分析，这些信息对规划和调度部门的决策极具参考价值，因而广泛应用于中长期电网规划和短期运行规划、状态估计及量测点布置、输电系统输送容量和无功规划等。在电力市场环境下，由于发电竞价上网、输电转运等因素，潮流分布的不确定性增大，随机潮流计算将成为日常和必备的分析工具。分布式发电的出现也使电力系统的负荷预测、规划和运行与过去相比有更大的不确定性。由于大量的用户会安装分布式发电为其提供电能，使得配电网规划人员更加难于准确预测负荷的增长情况，从而影响后续的规划。电力系统的随机潮流计算的研究对于整个电网运行与规划都具有重要意义。 关键词：电力系统；随机潮流；潮流计算方法；随机潮流现状 1. 电力系统随机潮流计算的产生 1.1 随机潮流研究背景 潮流计算主要目的是对电力系统中节点电压分布和支路功率分布进行计算，计算得到的结果能够为指挥电网运行、选择导线截面积和输变电设备、检验和确定网络拓扑结构等提供依据，因此它是电力系统中最重要最基本的计算。目前大部分的潮流计算是确定性的计算，也就是说它只能针对网络结构、节点注入量以及部分节点电压和变压器变比均为已知的确定的值,这种确定的运行方式来计算各节点电压和支路潮流以及其它待求量。但是在实际的电力系统中存在着许多的随机因素，这就使得网络结构和节点注入量等都具有一定的随机性。在这种情况下，如果仍然使用确定性的潮流计算就需要对系统众可能发生的情况分别进行统计与分析研究。容易看出，想要把所有可能发生的情况都做相应的计算不现实也没有必要，即便如此，不仅计算量相当大，计算的结果也不一定能够令人满意。目前电网中大量不确定因素主要有： （1）在实际运行环境中，当前系统运行状态是通过仪表测到的值来描述的，在使用仪表量测过程中难免有量测误差的存在； （2）在电网设计和规划过程中，要规划设计若干年以后的电源和电网的发展，因此系统负荷预测值不可能很准确；系统负荷也不再是一个通常意义下的已知的确定的量，而是一个随机变量。 （3）从严格意义上说，有时也需要将发电机输出功率作为一个随机变量来处理，因为发电机也不是百分之百可靠的，也会出现因故障而退出运行的可能性。 （4）网络中的线路和变压器等输电设备并不是完全可靠地，它们的随机故障导致系统运行方式发生变化，因此网络拓扑结构并不确定。（5）近年来随着国家能源结构的调整，可再生性的清洁能源得到了重视并在电力系统中所占的比例不断地增加。由于很多新能源发电的有功出力受自然天气条件的影响很大而具有随机性、间歇性和不可调度性的缺点。它们的接入给电网带来随机性的扰动，对电网的影响也日益突出。（6）电力用户中，尤其是工业用户，对电价的变化十分敏感。他们的工作进度可能安排在合适的电价下，例如，在远离高峰负荷的时间里安排工作生产时间。由于这种调整用电时间的行为，使电力系统规划运行 由于本人参加我们电气学院的电气小课堂，主讲的是计算机算法计算潮流这章，所以潜心玩了一个星期，下面整理给大家分享下。 本人一个星期以来的汗水，弄清楚了计算机算法计算潮流的基础，如果有什么不懂的可以发信息到邮箱：zenghao616@https://www.wendangku.net/doc/263818452.html, 接下来开始弄潮流的优化问题，吼吼！ 电力系统的潮流计算的计算机算法：以MATLAB为环境 这里理论不做过多介绍，推荐一本专门讲解电力系统分析的计算机算法的书籍---------《电力系统分析的计算机算法》—邱晓燕、刘天琪编著。 这里以这本书上的例题【2-1】说明计算机算法计算的过程，分别是牛顿拉弗逊算法的直角坐标和极坐标算法、P-Q分解算法。主要是简单的网络的潮流计算，其实简单网络计算和大型网络计算并无本质区别，代码里面只需要修改循环迭代的N即可，这里旨在弄清计算机算法计算潮流的本质。代码均有详细的注释. 其中简单的高斯赛德尔迭代法是以我们的电稳教材为例子讲，其实都差不多，只要把导纳矩阵Y给你，节点的编号和分类给你，就可以进行计算了，不必要找到原始的电气接线图。 理论不多说，直接上代码： 简单的高斯赛德尔迭代法： 这里我们只是迭代算出各个节点的电压值，支路功率并没有计算。 S_ij=P_ij+Q_ij=V_i(V_i* - V_j*) * y_ij* 可以计算出各个线路的功率 在显示最终电压幅角的时候注意在MATLAB里面默认的是弧度的形式，需要转化成角度显示。 clear;clc; %电稳书Page 102 例题3-5 %计算网络的潮流分布 --- 高斯-赛德尔算法 %其中节点1是平衡节点 %节点2、3是PV节点，其余是PQ节点 % 如果节点有对地导纳支路 %需将对地导纳支路算到自导纳里面 %------------------------------------------------% %输入原始数据，每条支路的导纳数值，包括自导和互导纳； y=zeros(5,5); y(1,2)=1/(0.0194+0.0592*1i); y(1,5)=1/(0.054+0.223*1i); y(2,3)=1/(0.04699+0.198*1i); y(2,4)=1/(0.0581+0.1763*1i); %由于电路网络的互易性，导纳矩阵为对称的矩阵 电 力 系 统 三 种 潮 流 计 算 方 法 的 比 较 一、高斯 -赛德尔迭代法： 以导纳矩阵为基础， 并应用高斯 -- 塞德尔迭代的算法是在电力系统中最早得到应用的潮流计算方法，目前高斯一塞德尔法已很少使用。 将所求方程 f ( x ) 0 改写为 x ( x ) 不能直接得出方程的根，给一个猜测值 x 0 得 x 1( x 0 ) 又可取 x1 为猜测值，进一步得： x 2 ( x 1 ) 反复猜测 x k 1 迭代 则方程的根 ( x k ) 优点： 1. 原理简单，程序设计十分容易。 2. 导纳矩阵是一个对称且高度稀疏的矩阵，因此占用内存非常节省。 3. 就每次迭代所需的计算量而言，是各种潮流算法中最小的，并且和网络所包 含的节点数成正比关系。 缺点： 1. 收敛速度很慢。 2. 对病态条件系统，计算往往会发生收敛困难：如节点间相位角差很大的重负 荷系统、包含有负电抗支路 (如某些三绕组变压器或线路串联电容等 )的系统、具有较长的辐射形线路的系统、长线路与短线路接在同一节点上，而且长短 线路的长度比值又很大的系统。 3. 平衡节点所在位置的不同选择，也会影响到收敛性能。 二、牛顿 -拉夫逊法： 求解 f ( x ) 0 设 x x 0 x ，则 按牛顿二项式展开： 当 △x 不大，则取线性化（仅取一次项） 则可得修正量 对 得： 作变量修正： x k 1x k x k ，求解修正方程 20 世纪 牛顿法是数学中求解非线性方程式的典型方法，有较好的收敛性。自从 60 年代中期采用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、计算速度方面都超过了其他方法，成为直到目前仍被广泛采用的方法。 优点： 1. 收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭 代 4—5 次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。 2. 具有良好的收敛可靠性， 对于前面提到的对以节点导纳矩阵为基础的高斯一 塞德尔法呈病态的系统，牛顿法均能可靠地收敛。 3. 牛顿法所需的内存量及每次迭代所需时间均较前述的高斯一塞德尔法为多， 并与程序设计技巧有密切关系。 缺点： 电力系统潮流分析 —基于牛拉法和保留非线性的随机潮流 ， 姓名：*** 学号：*** 1 潮流算法简介 常规潮流计算 常规的潮流计算是在确定的状态下。即：通过已知运行条件（比如节点功率或网络结构等）得到系统的运行状态（比如所有节点的电压值与相角、所有支路上的功率分布和损耗等）。 常规潮流算法中的一种普遍采用的方法是牛顿-拉夫逊法。当初始值和方程的精确解足够接近时，该方法可以在很短时间内收敛。下面简要介绍该方法。 牛顿拉夫逊方法原理 对于非线性代数方程组式（1-1），在待求量x 初次的估计值(0)x 附近，用泰勒级数（忽略二阶和以上的高阶项）表示它，可获得如式（1-2）的线性化变换后的方程组，该方程组被称为修正方程组。'()f x 是()f x 对于x 的一阶偏导数矩阵，这个矩阵便是重要的雅可比矩阵J 。 12(,,,)01,2, ,i n f x x x i n == （1-1） (0)'(0)(0)()()0f x f x x +?= （1-2） ' 由修正方程式可求出经过第一次迭代之后的修正量(0)x ?，并用修正量(0)x ?与估计值(0) x 之和，表示修正后的估计值(1)x ，表示如下（1-4）。 (0)'(0)1(0)[()]()x f x f x -?=- （1-3） (1)(0)(0)x x x =+? （1-4） 重复上述步骤。第k 次的迭代公式为： '()()()()()k k k f x x f x ?=- （1-5） (1)()()k k k x x x +=+? （1-6） 当采用直角坐标系解决潮流方程，此时待解电压和导纳如下式： i i i ij ij ij V e jf Y G jB =+=+ （1-7） 假设系统的网络中一共设有n 个节点，平衡节点的电压是已知的，平衡节点表示如下。 n n n V e jf =+ （1-8） } 毕业设计（论文）题目配电网潮流计算与程序设计 完成日期年月日 目录 1绪论 (4) 1.1潮流计算的目的和意义 (4) 1.2潮流计算的发展 (4) 1.3本文的主要工作 (6) 2 配电网络模型 (6) 2.1线路模型 (6) 2.2变压器的模型 (9) 2.3 负荷模型 (13) 3基于前推回代法的潮流计算 (14) 4基于matlab的程序设计 (15) 4.1 MATLAB简介 (15) 4.2 程序设计 (16) 4.3 算例 (20) 5 总结 (24) 配电网潮流计算与程序设计 摘要：在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。为了保证电力系统的稳定运行，要进行潮流调节。 Abstract:In the normal operation of the power system, with the change of power load and the change of system operation mode, the network losses will also change. To strictly ensure that all users at any time have a rated voltage is impossible, so the node appears in the operation of the system voltage deviation is inevitable. In order to guarantee the stable operation of power system, must carry on the trend of the regulation. 关键词: 电网潮流计算前推后代法MATLAB Keywords:Power flow calculation Push the offspring before MATLAB 前言 电力是以电能作为动力的能源。发明于19世纪70 年代，电力的发明和应用掀起了第二次工业化高潮。成为人类历史18世纪以来，世界发生的三次科技革命之一，从此科技改变了人们的生活。既是是当今的互联网时代我们仍然对电力有着持续增长的需求，因为我们发明了电脑、家电等更多使用电力的产品。不可否认新技术的不断出现使得电力成为人们的必需品。 20世纪出现的大规模电力系统是人类工程科学史上最重要的成就之一，是由发电、输电、变电、配电和用电等环节组成的电力生产与消费系统。它的功能是将自然界的一次能源通过发电动力装置转化成电能，再经输电、变电和配电将电能供应到各用户。为实现这一功能，电力系统在各个环节和不同层次还具有相应的信息与控制系统，对电能的生产过程进行测量、调节、控制、保护、通信和调度，以保证用户获得安全、经济、优质的电能。 电力工业发展初期，电能是直接在用户附近的发电站（或称发电厂）中生产的，各发电站孤立运行。随着工农业生产和城市的发展，电能的需要量迅速增加，而热能资源和水能资源丰富的地区又往往远离用电比较集中的城市和工矿区，为了解决这个矛盾，就需要在动力资源丰富的地区建立大型发电站，然后将电能远距离输送给电力用户。同时，为了提高供电的可靠性以及资源利用的综合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来的整体，即称为电力系统。 随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。电力系统潮流计算是电力系统分析计算中最基本的内容，也是电力系统运行及设计中必不可少的工具。根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线电压的幅值及相角、各元件中流过 潮流计算数学模型与数值方法 1. 什么是潮流计算潮流计算的主要作用有哪些 潮流计算，电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。 潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。 2. 潮流计算有哪些待求量、已知量 （已知量：1、电力系统网络结构、参数 2、决定系统运行状态的边界条件 待求量：系统稳态运行状态 例如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等） 3. 潮流计算节点分成哪几类分类根据是什么 （分成三类：PQ 节点、PV 节点和平衡节点，分类依据是给定变量的不同） 4. 教材牛顿-拉夫逊法及有功-无功分解法是基于何种电路方程可否采用其它类型方程 答:基于节点电压方程，还可以采用回路电流方程和割集电压方程等。但是后两者不常用。 5. 教材牛顿-拉夫逊法是基于节点阻抗方程、还是基于节点导纳方程进行迭代计算的试阐述这两种方程的优点与缺点。 1.不能由等值电路直接求出 2.满秩矩阵内存量大 3.对角占优矩阵。。 节点导纳矩阵的特点：1.直观容易形成2.对称阵3.稀疏矩阵（零元素多）：每一行的零元素个数=该节点直接连出的支路数。 6. 说出至少两种建立节点导纳矩阵的方法，阐述其中一种方法的原理与过程。 方法：1.根据自导纳和互导纳的定义直接求取2.运用一节点关联矩阵计算3.阻抗矩阵的逆矩阵 节点导纳矩阵的形成：1.对角线元素ii Y 的求解)1,,0(=≠==i j I i ii U i j U U I Y 【除i 外的其他节点接地，0=j U ，只在i 节点加单位电压值】解析ii Y 等于与i 节点直接相连的的所有支路导纳和2.互导纳),0,1(j k U U U I Y k j j i ij ≠===，ji ij Y Y =（无源网络导纳之间是对称的）解析：ij Y 等于j i ,节点之间直接相连的支路导纳的负值。 7. 潮流计算需要考虑哪些约束条件 答: 为了保证系统的正常运行必须满足以下的约束条件： 对控制变量 潮流计算的基本算法及使用方法 一、 欧阳光明（2021.03.07） 二、 潮流计算的基本算法 1. 牛顿－拉夫逊法 1．1 概述 牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。 牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。 1．2一般概念 对于非线性代数方程组 即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1) 在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并 略去二阶及以上的高阶项，得到如下的线性化的方程组 ()()()() ()0000=?'+x x f x f (1－2) 上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 ()() ()[]()()0 1 00x f x f x -'-=? (1－3) 将()0x ?和()0x 相加，得到变量的第一次改进值()1x 。接着再从()1x 出 发，重复上述计算过程。因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为 ()()()() ()k k k x f x x f -=?' (1－4) ()()()k k k x x x ?+=+1 (1－5) 上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。 由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成 求解修正方程式。牛顿法当初始估计值()0x 和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。 1．3潮流计算的修正方程 运用牛顿－拉夫逊法计算潮流分布时，首先要找出描述电力系统的非线性方程。这里仍从节点电压方程入手，设电力系统导纳矩阵已知，则系统中某节点（i 节点）电压方程为 从而得 ∑=* *? ?=n j j ij i i U Y U S 1 进而有 ()01 =-+* =* ? ∑j n j ij i i i U Y U jQ P （1 －6） 潮流计算的定义（课后题） 各种潮流计算模型和算法的特点、适用范围以及相互之间的区别和联系（课后题） 影响潮流收敛性的因素，以及如何改善潮流计算的收敛性（课后题） 通过功率方程说明为什么潮流计算的数学模型是非线性的？应该采用什么样的数学方法求解？（03A 、05A ） 电力系统的潮流计算有哪些常规算法？有哪些扩展算法？（05B ） 潮流计算的目的是什么？其数学模型是什么？有何特点？（06B ） 简要说明潮流计算的概念、模型及计算方法。（07B ） 高斯赛德尔迭代法和牛顿拉夫逊迭代法是常规的潮流计算方法，请介绍一下最优潮流（OPF ）算法的原理及其应用。（04电科院） 潮流计算的目的： 常规潮流计算的目的是在已知电力网络参数和各节点的注入量的条件下，求解各节点电压。 目的1： 1. 在电网规划阶段，通过潮流计算，合理规划电源容量和接入点，合理规划网架，选择无 功补偿方案，满足规划水平年的大小方式下潮流交换控制、调峰、调相、调压的要求。 2. 在编制年运行方式，在预计复合增长及新设备投运基础上，选择典型方式进行潮流计算， 发现电网中的薄弱环节，供调度人员异常调度控制参考，并对规划、基建部门提出改进网架结构，加快基建进度的建议。 3. 正常检修及特殊运行方式下的潮流计算，用于日常运行方式的编制，指导发电厂开机方 式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求及电压质量要求。 4. 预想事故、设备退出运行对静态安全分析的影响及做出预想的运行方式调整方案。 目的2： A. 检查电力系统各元件是否过负荷； B. 检查电力系统各节点的电压是否满足电压质量的要求； C. 根据对各种运行方式的潮流分布计算，可以正确的选择系统接线方式，合理调整负荷， 以保证电力系统安全、可靠地运行，向用户供给高质量的电能； D. 根据功率分布，可以选择电力系统的电气设备和导线截面积，可以为电力系统继电保护 整定计算提供必要的数据等； E. 为电力系统扩建和规划提供依据； F. 为调压计算、经济运行计算、短路计算等提供必要的数据。 数学模型：数学模型为：潮流计算所用的电力网络系由变压器、输电线路、电容器、电抗器等静止线性元件所构成，并用集中参数表示的串联或并联等值支路来模拟。普遍采用节点法，I YU =来建立潮流计算的数学模型。在实际工程中，节点注入量不是电流，而是节点功率，因此节点电压方程要进行修改：*,(1,2,...,)i i i P jQ I i n U -==，进一步得到 * *1 ,(1,2,...,)n i i ij j j i P jQ Y U i n U =-==∑，上式为电压的非线性隐函数，无法直接求解，必须通过一定的算法求近似解。这是潮流计算问题最基本的方程式，是一个以节点电压U ?为变量 3.3复杂电力网潮流计算的计算机解法 3.3.1 导纳矩阵的形成 1．自导纳 节点i的自导纳，亦称输入导纳，在数值上等于在节点i施加单位电压，其他节点全部接地时，经节点i注入网络的电流。 主对角线元素，更具体地说，就等于与节点连接的所有支路导纳的和。 2．互导纳 节点i、j间的互导纳，在数值上等于在节点i施加单位电压，其他节点全部接地时，经节点j注入网络的电流。非对角线元素。 更具体地说，是连接节点j和节点i支路的导纳之和再加上负号而得。 3．导钠矩阵的特点： （1）因为，导纳矩阵Y是对称矩阵； （2）导纳矩阵是稀疏矩阵，每一非对角元素是节点i和j间支路导纳的负值，当i和j间没有直接 相连的支路时，即为零，根据一般电力系统的特点，每一节点平均与3-5个相邻节点有直接联系，所以导纳矩阵是一高度稀疏的矩阵； （3）导纳矩阵能从系统网络接线图直观地求出。 4．节点导纳矩阵的修改 （1）从原有网络引出一支路，同时增加一节点，设i为原有网络结点，j为新增节点，新增支路ij的导纳为y ij。如图3-17（a）所示。 因新增一节点，新的节点导纳阵需增加一阶。且新增对角元Y jj=y ij，新增非对角元Y ij=Y ji=－y ij，同时对原阵中的对角元Y ii进行修改，增加ΔY ii＝y ij。 （2）在原有网络节点i、j间增加一支路。如图3-17（b）所示。 设在节点i增加一条支路，由于没有增加节点数，节点导纳矩阵Y阶次不变，节点的自导纳Y ii、Y jj和互导纳Y ij分别变化量为 （3-57） 图 3-17 网络接线的变化图 （a）网络引出一支路，（b）节点间增加一支路，（c）节点间切除一支路，（d）节点间导纳改变 （3）在原有网络节点i、j间切除一支路。如图3-17（c）所示。 设在节点i切除一条支路，由于没有增加节点数，节点导纳矩阵Y阶次不变，节点的自导纳Y ii、Y jj和互导纳Y ij分别发生变化，其变化量为 毕业设计（论文）任务书 信息与电气工程系电工电子基础教研室 系（教研室）主任:（签名）年月日 学生姓名:学号:专业:电气工程及其自动化 1 设计（论文）题目及专题：电力系统概率潮流计算的计算方法与比较 2 学生设计（论文）时间：自 2020年1月9日开始至 2020年5月 25日止 3 设计（论文）所用资源和参考资料： [1]陈倪.电力系统概率潮流计算[D].东南大学,1990. [2]戴小青. 电力系统概率潮流新算法及其应用[D].华北电力大学(北京),2006. [3] 张建芬,王克文,宗秀红,谢志棠.几种概率潮流模型的准确性比较分析[J].郑州大学学报(工学版),2003(04):32-36. [4]DING Ming, LI Shenghu, HUANG Kai. Probabilistic load flow analysis based on Monte-Carlo simulation [J]. Power System Technology, 2007, 25(11):10-14. [5]MORALES J M,BARINGO L,CONEJO A J, et al. Probabilistic power flowwith correlated wind sources[J]. IET Generation, Transmission & Distribution, 2010, 4(5):641-651. [6]代景龙,韦化,鲍海波,等. 基于无迹变换含分布式电源系统的随机潮流[J]. 电力自动化设备, 2016, 36(3):86-93. [7]张衡,程浩忠,曾平良,等. 分位数拟合的点估计法随机潮流在输电网规划中的应用[J]. 电力自动化设备, 2018, 38(11):43-49. [8]方斯顿,程浩忠,徐国栋,等. 基于Nataf变换和准蒙特卡洛模拟的随机潮流计算[J]. 电力自动化设备, 2015, 35(8):38-44. 4 设计（论文）应完成的主要内容： （1）电力系统概率潮流概述；（2）基于蒙特卡洛法的概率潮流计算；（3）基于累积量法的概率潮流计算；（4）基于点估计法的概率潮流计算；（5）基于无迹变换法的概率潮流计算；（6）各种计算方法的比较分析。 近年来，大多数研究都是围绕改进牛顿法和P-Q分解法进行。此外，随着人工智能理论的发展，遗传算法、人工神经网络、模糊算法也逐渐引入潮流计算。但是，到目前为止这些新的模型和算法还不能取代牛顿法和P-Q分解法的地位。由于电力系统规模不断扩大，对计算速度要求不断提高，计算机的并行计算技术也将在潮流计算中得到广泛应用，成为重要的研究领域。 经过三十多年的发展，潮流算法已经比较成熟，但是仍存在不少尚待解决的问题。例如各种牛顿法潮流算法，对于某些条件可能导致不收敛。潮流计算的多解现象及其机理在重负荷情况下，临近多根与电压不稳定问题的关联。当前无论在实践上还是在理论上，均有许多问题需待解决，特别是如何快速求解成千上万个变量的大规模非线性规划问题。 近几年，对潮流计算的研究仍然是如何改善传统的潮流算法。牛顿拉夫逊法，由于其在求解非线性潮流方程时采用逐次线性化方法，为了进一步提高算法的收敛性和计算速度，人们考虑采用将泰勒级数高阶项或非线性项也考虑进来，于是产生了二阶潮流算法。后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点，提出了采用直角坐标的保留非线性快速潮流算法。 在这种情况下，进行电力系统规划和运行条件分析时，若不考虑随机变化因素，就要对众多可能发生的情况作大量的方案计算，计算时间是难以承受的，并且很难反映系统整体的状况。随机潮流计算是解决上述问题的有效方法和手段。应用随机理论来描述这种不确定性，探讨相应的数学建模，计算机算法和实际应用，称为随机潮流（Probabilistic Load Flow,简写为PLF）研究，也称为随机潮流。采用随机潮流计算方法，输入数据为已知的随机变量，给定的是它们的随机统计特性（例如，给定节点注入功率的期望和方差或随机密度函数等），输出数据则是节点电压和支路潮流的统计特性，有期望值和方差或随机密度函数等。由这些结果，可以知道节点电压、支路功率、PV节点无功功率及平衡节点功率的平均值、取值范围以及其随机等。这样，只要通过一次计算就能为电力系统的运行条件提供更完备的信息，减少了大量的计算工作量。根据这些信息，可以更深刻地揭示系统运行状况、存在问题和薄弱环节，为规划与运行决策提供更全面的信息，可以更恰当地确定输电线和无功补偿装置的容量以及系统的备用容量等，从而提高了电力系统的安全运行水平。 1.3.1国外关于随机潮流计算的研究现状 把随机分析方法应用在电力系统的潮流研究上来最初是 B.Borkowska在1974年提出来的。自从那以后，就有两种方法采用了随机分析方法来研究潮流问题：随机潮流方法和随机潮流方法。在随机潮流研究中，负荷和发电量在ti 瞬间被看成随机变量。这种方法研究了这种不确定性在每个瞬间给传统的潮流计潮流计算的计算机算法
关于电力系统经济调度的潮流计算分析
潮流计算的基本算法及使用方法Word版
潮流计算
潮流计算的计算机算法
随机潮流研究现状
(完整word版)潮流计算方法
有关电力系统三种潮流计算方法的比较.docx
电力系统潮流计算方法分析
潮流计算设计
潮流计算简答题
2021年潮流计算的基本算法及使用方法
潮流计算问题
33复杂电力网潮流计算的计算机解法
电力系统概率潮流计算的计算方法与比较毕业设计任务书
潮流计算国内外现状
- 潮流计算的计算机方法
- 潮流计算的计算机算法
- 潮流计算的计算机算法课程设计范文 (2)
- 电力系统分析潮流计算的计算机算法
- 潮流计算的计算机算法
- 电力系统潮流的计算机算法
- 复杂电力系统潮流的计算机算法分享
- 潮流计算的计算机算法资料
- 电力系统分析潮流计算的计算机算法
- 潮流计算的计算机算法
- 第四章复杂电力系统潮流的计算机算法ppt-PowerPo
- 电力系统潮流的计算机算法夏道止版
- 潮流计算的计算机算法课程设计范文
- 潮流计算的计算机算法
- 第四章复杂电力系统潮流的计算机算法ppt-PowerPo
- 电力系统分析潮流计算的计算机算法
- 潮流的计算机算法ppt课件
- 电力系统潮流计算的计算机算法概述
- 电力系统潮流的计算机算法
- 潮流的计算机算法课件PPT