计算机潮流计算资料

合集下载

[理学]4 潮流计算

[理学]4 潮流计算

方程数: n 1 m n m 1 2(n 1)
i
未知量: ei , f i , i ( PQ PV ) , 2(n 1)
第四讲 电力系统潮流计算
(2)极坐标下的数学方程
将 Vi Vi e
j i
和 Yij Gij jBij 代入 Pi jQi Vi Yij V j
n个节点,m个PQ节点,n m 1个PV节点
编号
1 m
PQ
m 1 n 1
PV
n
平衡点
I YV
I i Yij V j
j 1 n
Vi I i Vi Y ij V j
j 1

n


Pi jQi Vi Yij V j
x(k )
x
讨论收敛区域和收敛条件。又称切线法。
第四讲 电力系统潮流计算
2、设初始点 xo , f ( xo ) 0
f ( xo x ) 0 df f ( xo ) dx f ( xo ) x df dx 1 d2 f x 2 dx 2 x 0
x0
x 2
( 0) x1 0 ( 0) x 2
(1) ( 0) ( 0) x1 x1 x1 (1) ( 0) ( 0) x2 x2 x2
第四讲 电力系统潮流计算
3、多维非线性方程组的迭代公式
基于同样的思想,我们可以得到 n维非线性方程 —牛顿 拉夫逊迭代公式
j 1
j i j i ( G jB ) V e ij ij i j 1 n

n
*

(完整word版)潮流计算方法

(完整word版)潮流计算方法

由于本人参加我们电气学院的电气小课堂,主讲的是计算机算法计算潮流这章,所以潜心玩了一个星期,下面整理给大家分享下。

本人一个星期以来的汗水,弄清楚了计算机算法计算潮流的基础,如果有什么不懂的可以发信息到邮箱:zenghao616@接下来开始弄潮流的优化问题,吼吼!电力系统的潮流计算的计算机算法:以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);%由于电路网络的互易性,导纳矩阵为对称的矩阵for i=1:1:5for j=1:1:5y(j,i)=y(i,j);endend%节点导纳矩阵的形成Y=zeros(5,5);%求互导纳for i=1:1:5for j=1:1:5if i~=jY(i,j)=-y(i,j);endendend%求自导纳for i=1:1:5%这句话是说将y矩阵的第i行的所有元素相加,得到自导纳的值 Y(i,i)=sum(y(i,:));end%上面求得的自导纳不包含该节点的对地导纳数值,需要加上Y(2,2)=Y(2,2)+0.067*1i;Y(3,3)=Y(3,3)+0.022*1i;Y(4,4)=Y(4,4)+0.0187*1i;Y(5,5)=Y(5,5)+0.0246*1i;%导纳矩阵的实部和虚部G = real(Y);B = imag(Y);Qc2=0;Qc3=0;%原始节点功率%这里电源功率为正,负荷功率为负S(1)=0;S(2)=-0.217-0.121*1i+Qc2*1i;S(3)=-0.749-0.19*1i+Qc3*1i;S(4)=-0.658+0.039*1i;S(5)=-0.076-0.016*1i;%节点功率的P QP = real(S);Q = imag(S);%下面是两个PV节点的无功初始值Q(2) = 0;Q(3) = 0;U=ones(5,1); %1列5行的‘1’矩阵%节点电压初始值U(1)=1.06;U(2)=1.045;U(3)=1.01;U_reg=U;Sum_YU0=0;%中间变量Sum_YU1=0;%中间变量for cont=1:1:6 %这里的cont是迭代次数for i=2:1:5for j=1:1:iif i~=jSum_YU0 = Sum_YU0 + Y(i,j)*U_reg(j);endendfor j=i+1:1:5Sum_YU1 = Sum_YU1 + Y(i,j)*U(j);endU(i)=( (P(i)-Q(i)*1i ) / conj(U(i)) - Sum_YU0 - Sum_YU1 ) / Y(i,i); U_reg(i)=U(i);%PV节点计算%下面是把求出的U2、U3只保留其相位,幅值不变if i==2angle_U2 = angle(U(2));U(2)=1.045*cos(angle_U2)+1.045*sin(angle_U2)*1i;Q(2)=imag( U(2)*( conj(Sum_YU0) + conj(Sum_YU1) + conj(Y(2,2)*U(2)) ) );endif i==3angle_U3 = angle(U(3));U(3)=1.01*cos(angle_U3)+1.01*sin(angle_U3)*1i;Q(3)=imag( U(3)*( conj(Sum_YU0) + conj(Sum_YU1) + conj(Y(3,3)*U(3)) ) );end% 下面做越界检查%if Q(4)>Q_Max% Q(4) = Q_Max;%end%if Q(4)<Q_Min% Q(4) = Q_Min;%end%下面可以做PV节点收敛判断Sum_YU0 = 0;Sum_YU1 = 0;endend%节点注入无功,流入为正,流出为负Qc2=Q(2)+0.121-1.045^2 * 0.067;Qc3=Q(3)+0.19-1.01^2 * 0.022;%电压幅值和相角angle_U=angle(U)*180/pi;U=abs(U);S_Line=zeros(5,5);%计算平衡节点功率S_BalanceNode=0;for j=1:1:5S_BalanceNode = S_BalanceNode + U(1) * conj(Y(1,j)*U(j));end%下面由上面算出的电压值求线路的功率%这里计算出来的线路功率的有功、无功%for i=1:1:5% for j=i:1:5% if i~=j% S_Line(i,j)=U(i)*( conj(U(i))-conj(U(j)) ) * conj(y(i,j));% end% if i==2% %S_Line(2,j)=S_Line(2,j)+U(2)*conj(0.067*1i);% end% if i==3% %S_Line(3,j)=S_Line(3,j)+U(3)*conj(0.022*1i);% end% end%end计算网络的潮流分布 ---- Newton算法(直角坐标)clear;clc;%电稳书Page 102 例题3-5%计算网络的潮流分布 ---- Newton算法(直角坐标)%其中节点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);%由于电路网络的互易性,导纳矩阵为对称的矩阵for i=1:1:5for j=1:1:5y(j,i)=y(i,j);endend%节点导纳矩阵的形成Y=zeros(5,5);%求互导纳for i=1:1:5for j=1:1:5if i~=jY(i,j)=-y(i,j);endendend%求自导纳for i=1:1:5%这句话是说将y矩阵的第i行的所有元素相加,得到自导纳的值Y(i,i)=sum(y(i,:));end%上面求得的自导纳不包含该节点的对地导纳数值,需要加上Y(2,2)=Y(2,2)+0.067*1i;Y(3,3)=Y(3,3)+0.022*1i;Y(4,4)=Y(4,4)+0.0187*1i;Y(5,5)=Y(5,5)+0.0246*1i;%导纳矩阵的实部和虚部G = real(Y);B = imag(Y);%节点2、3需补偿的无功Qc2=0;Qc3=0;%原始节点功率%这里电源功率为正,负荷功率为负S(1)=0;S(2)=-0.217-0.121*1i+Qc2*1i;S(3)=-0.749-0.19*1i+Qc3*1i;S(4)=-0.658+0.039*1i;S(5)=-0.076-0.016*1i;%节点功率的P QP = real(S);Q = imag(S);%下面是两个PV节点的无功初始值Q(2) = 0;Q(3) = 0;%给点电压初始值e=[1.06,1.045,1.01,1,1];f=[0,0,0,0,0];U=e+f*1i;delta_U=zeros(1,5);delta_P=zeros(1,5);delta_Q=zeros(1,5);delta_PQV=ones(8,1);Sum_GB1=0;Sum_GB2=0;cont=0;while max(delta_PQV > 1e-6),cont=cont+1;%for cont=1:1:3%下面开始计算delta_P/delta_Q/delta_Ufor i=2:1:5for j=1:1:5Sum_GB1=Sum_GB1 + ( G(i,j)*e(j) - B(i,j)*f(j) );Sum_GB2=Sum_GB2 + ( G(i,j)*f(j) + B(i,j)*e(j) );enddelta_P(i)=P(i)-e(i)*Sum_GB1-f(i)*Sum_GB2;if i~=2 && i~=3 %不为节点2,3则计算无功delta_Q(i)=Q(i)-f(i)*Sum_GB1+e(i)*Sum_GB2;endif i==2 || i==3 %这里计算delta_U的值,始终为零delta_U(i)=U(i)^2-( e(i)^2 + f(i)^2 );endSum_GB1=0;Sum_GB2=0;end%___________________________________%%下面计算雅克比矩阵J=zeros(8,8);for ii=2:1:5i=ii-1;for j=1:1:5Sum_GB1=Sum_GB1 + ( G(ii,j)*e(j) - B(ii,j)*f(j) );Sum_GB2=Sum_GB2 + ( G(ii,j)*f(j) + B(ii,j)*e(j) );endfor jj=2:1:5j=jj-1;if ii~=2 && ii~=3 %PQ节点if ii==jjJ(2*i-1,2*i-1)=-Sum_GB1-G(ii,ii)*e(ii)-B(ii,ii)*f(ii);J(2*i-1,2*i)=-Sum_GB2+B(ii,ii)*e(ii)-G(ii,ii)*f(ii);J(2*i,2*i-1)=Sum_GB2+B(ii,ii)*e(ii)-G(ii,ii)*f(ii);J(2*i,2*i)=-Sum_GB1+G(ii,ii)*e(ii)+B(ii,ii)*f(ii);elseJ(2*i-1,2*j-1)=-(G(ii,jj)*e(ii)+B(ii,jj)*f(ii)); J(2*i-1,2*j)=B(ii,jj)*e(ii)-G(ii,jj)*f(ii);J(2*i,2*j-1)=B(ii,jj)*e(ii)-G(ii,jj)*f(ii);J(2*i,2*j)=(G(ii,jj)*e(ii)+B(ii,jj)*f(ii));endelse%PV节点if ii==jjJ(2*i-1,2*i-1)=-Sum_GB1-G(ii,ii)*e(ii)-B(ii,ii)*f(ii);J(2*i-1,2*i)=-Sum_GB2+B(ii,ii)*e(ii)-G(ii,ii)*f(ii);J(2*i,2*i-1)=-2*e(ii);J(2*i,2*i)=-2*f(ii);elseJ(2*i-1,2*j-1)=-(G(ii,jj)*e(ii)+B(ii,jj)*f(ii)); J(2*i-1,2*j)=B(ii,jj)*e(ii)-G(ii,jj)*f(ii);J(2*i,2*j-1)=0;J(2*i,2*j)=0;endendendSum_GB1=0;Sum_GB2=0;end%在求解修正方程之前建议把delta_P和delta_Q,delta_U全部放在一个矩阵delta_PQV=[delta_P(2);delta_U(2);delta_P(3);delta_U(3);delta_P(4) ;delta_Q(4);delta_P(5);delta_Q(5)];%下面求解修正方程;注意矩阵运算时候的左除和右除的区别delta_ef=-J\delta_PQV;%下面修正各个节点的电压for i=2:1:5e(i)=e(i)+delta_ef(2*(i-1)-1);f(i)=f(i)+delta_ef(2*(i-1));end%到这里第一轮迭代完成end%电压幅值和相角U=e+f*1i;angle_U=angle(U)*180/pi;%节点注入无功,流入为正,流出为负Sum_YU=0;for i=2:1:3for j=1:1:5Sum_YU = Sum_YU + Y(i,j)*U(j);endQ(i)=imag( U(i)*conj( Sum_YU ) );Sum_YU=0;endQc2=Q(2)+0.121-1.045^2 * 0.067;Qc3=Q(3)+0.19-1.01^2 * 0.022;U=abs(U);disp(['Iteration times : ' num2str(cont)]);%显示最终的迭代次数牛顿算法求解潮流 (极坐标):clear;clc;%牛顿算法求解潮流 (极坐标)%计算网络的潮流分布%其中节点5是平衡节点%节点1、2、3是PQ节点,节点4是PV节点% 如果节点有对地导纳支路%需将对地导纳支路算到自导纳里面%------------------------------------------------%%输入原始数据,每条支路的导纳数值,包括自导和互导纳;Y=[0.8381-3.7899*1i,-0.4044+1.6203*1i,0,0,-0.4337+2.2586*1i;...-0.4044+1.6203*1i,0.7769-3.3970*1i,-0.3726+1.8557*1i,0,0;...0,-0.3726+1.8557*1i,1.1428-7.0210*1i,-0.5224+4.1792*1i,-0.2739+1. 2670*1i;...0,0,-0.5224+4.1792*1i,0.5499-4.3591*1i,0;...-0.4337+2.2586*1i,0,-0.2739+1.2670*1i,0,0.7077-3.4437*1i];%导纳矩阵的实部和虚部G = real(Y);B = imag(Y);%给点电压初始值U = [1,1,1,1,1.05];angle_U=[0,0,0,0,0];%for i=1:1:5% U(i)=U_abs(i)*cos(angle_U(i))+U_abs(i)*sin(angle_U(i))*1i;%end%原始节点功率%这里电源功率为正,负荷功率为负%下面给点PQ PV节点功率值S=[-0.22-0.14*1i,-0.18-0.09*1i,-0.27-0.13*1i,0.35,0];%节点功率的P QP = real(S);Q = imag(S);%下面是PV节点的无功初始值Q(4) = 0;delta_P=zeros(1,5);delta_Q=zeros(1,5);%delta_angleU=zeros(1,4);%delta_absU=zeros(1,4);delta_PQ=ones(8,1);Sum_GB1=0;Sum_GB2=0;cont=0;%最外层循环,cont代表迭代的次数,这里可以用约束条件来代替%for cont=1:1:4while max(delta_PQ)>1e-6,%下面计算delta_P/delta_Q/delta_Ucont=cont+1;for i=1:1:4for j=1:1:5Sum_GB1=Sum_GB1 + U(j)*( G(i,j)*cos(angle_U(i)-angle_U(j)) + B(i,j)*sin(angle_U(i)-angle_U(j)) );Sum_GB2=Sum_GB2 + U(j)*( G(i,j)*sin(angle_U(i)-angle_U(j)) - B(i,j)*cos(angle_U(i)-angle_U(j)) );enddelta_P(i)=P(i)-U(i)*Sum_GB1;if i~=4 %不为节点四则计算无功delta_Q(i)=Q(i)-U(i)*Sum_GB2;endSum_GB1=0;Sum_GB2=0;end%_______________________________________________________%%下面计算雅克比矩阵J=zeros(7,7);for ii=1:1:4for jj=1:1:4if ii ~= 4 %PQ节点if ii==jjJ(2*ii-1,2*ii-1)=U(ii)^2*B(ii,ii)+Q(ii);J(2*ii-1,2*ii)=-U(ii)^2*G(ii,ii)-P(ii);J(2*ii,2*ii-1)=U(ii)^2*G(ii,ii)-P(ii);J(2*ii,2*ii)=U(ii)^2*B(ii,ii)-Q(ii);elseJ(2*ii-1,2*jj-1)=-U(ii)*U(jj)*( G(ii,jj)*sin(angle_U(ii)-angle_U( jj)) - B(ii,jj)*cos(angle_U(ii)-angle_U(jj)) );J(2*ii-1,2*jj)=-U(ii)*U(jj)*( G(ii,jj)*cos(angle_U(ii)-angle_U(jj)) + B(ii,jj)*sin(angle_U(ii)-angle_U(jj)) );J(2*ii,2*jj-1)=U(ii)*U(jj)*( G(ii,jj)*cos(angle_U(ii)-angle_U(jj)) + B(ii,jj)*sin(angle_U(ii)-angle_U(jj)) );J(2*ii,2*jj)=-U(ii)*U(jj)*( G(ii,jj)*sin(angle_U(ii)-angle_U(jj)) - B(ii,jj)*cos(angle_U(ii)-angle_U(jj)) );endelse%PV节点if ii==jjJ(2*ii-1,2*ii-1)=U(ii)^2*B(ii,ii)+Q(ii);J(2*ii-1,2*ii)=-U(ii)^2*G(ii,ii)-P(ii);elseJ(2*ii-1,2*jj-1)=-U(ii)*U(jj)*( G(ii,jj)*sin(angle_U(ii)-angle_U( jj)) - B(ii,jj)*cos(angle_U(ii)-angle_U(jj)) );J(2*ii-1,2*jj)=-U(ii)*U(jj)*( G(ii,jj)*cos(angle_U(ii)-angle_U(jj)) + B(ii,jj)*sin(angle_U(ii)-angle_U(jj)) );endendendend%在求解修正方程之前建议把delta_ef和delta_ef全部放在一个矩阵delta_PQ=[delta_P(1);delta_Q(1);delta_P(2);delta_Q(2);delta_P(3); delta_Q(3);delta_P(4)];%下面求解修正方程;注意矩阵运算时候的左除和右除的区别J=J(1:7,1:7);delta_ef=-J\delta_PQ;%下面修正各个节点的电压for i=1:1:4if i~=4U(i)=U(i)+delta_ef(2*i)*U(i);endangle_U(i)=angle_U(i)+delta_ef(2*i-1);end%到这里第一轮迭代完成end%下面显示出满足条件后的迭代的次数disp(['Iteration times : ' num2str(cont)]);%下面计算平衡节点5的功率PQfor j=1:1:5Sum_GB1=Sum_GB1 + U(j)*( G(5,j)*cos(angle_U(5)-angle_U(j)) + B(5,j)*sin(angle_U(5)-angle_U(j)) );Sum_GB2=Sum_GB2 + U(j)*( G(5,j)*sin(angle_U(5)-angle_U(j)) - B(5,j)*cos(angle_U(5)-angle_U(j)) );endP(5)=U(5)*Sum_GB1;Q(5)=U(5)*Sum_GB2;%下面将相角用角度表示for i=1:1:5angle_U(i)=angle_U(i)*180/pi;End计算计算法P-Q算法计算潮流:这个算法是由牛顿算法的极坐标形式简化而来。

潮流计算方法

潮流计算方法

计算电流及功率
for n=1:nbus Sc = conj(conj(V(n))*(Y(n,n)*V(n))); DP(n) = P(n) - real(Sc); DQ(n) = Q(n) - imag(Sc); end

返回
Qi Lii Vi Vi V j (Gij sin ij Bij cos ij ) 2Vi 2 Bii Vi 2 Bii Qi ji i j i
直流潮流与牛拉法的结合
1.通过直流潮流计算各节点的相角。 2.将结果作为牛拉法初值进行计算。
J ij
Qi Vi V j (Gij cos ij Bij sin ij ) j
Qi Lij V j Vi V j (Gij sin ij Bij cos ij ) Vj
雅克比矩阵的形成
当j=i 时
H ii
Pi Vi V j (Gij sin ij Bij cos ij ) Vi 2 Bii Qi ji i j i

数据输入:
将各节点的初始值及线路数据输入excel文 件中,并通过xlsread这一函数输入到MATLAB 中。
形成导纳矩阵
1.形成一个零矩阵Y 2.检测线路参数前两列,将Y对应元素加上相应 参数。
形成导纳矩阵


1.非对角线元素:
for n=1:nbr Y(n1(n),n2(n))= Y(n1(n),n2(n))-Y1(n)/(a(k)^2); Y(n2(n),n1(n))= Y(n1(n),n2(n)); end
总体思路:
1.直接使用“牛顿——拉夫逊”法,通过图像 显示每次迭代后的误差并记录程序运行时间。 2.直接使用“PQ分解法”,通过图像显示每次 迭代后的误差并记录程序运行时间。 3.首先使用直流潮流计算,计算各节点电压的 相角,并将计算结果作为 “牛顿——拉夫逊” 法的初值,通过图像显示每次迭代后的误差并 记录程序运行时间。

第四章 潮流的计算机算法new-教材

第四章 潮流的计算机算法new-教材

y34 4 y440
7
将接在同一节点上的接 地导纳进行并联,得
I1
y12 y10 y13 y23 y20
I2
y10 = y120 + y130
⎫ ⎪ y20 = y210 + y230 ⎪ ⎬ y30 = y310 + y320 + y340 ⎪ ⎪ y40 = y440 + y430 ⎭
3 y30 y34 4 y40 I4
& ⎤ L Y1i L Y1 j L Y1n ⎤ ⎡U 1 ⎢ ⎥ ⎥ O M O M O M ⎥⎢ M ⎥ & ⎥ L Yii L Yij L Yin ⎥ ⎢U i ⎥⎢ ⎥ O M O M O M ⎥⎢ M ⎥ & ⎥ L Y ji L Y jj L Y jn ⎥ ⎢U ⎥⎢ j⎥ O M O M O M ⎥⎢ M ⎥ ⎢U ⎥ & L Yni L Ynj L Ynn ⎥ ⎥ ⎦⎢ ⎣ n⎥ ⎦
&⎤ L Z1i L Z1 j L Z1n ⎤ ⎡ I 1 ⎢ ⎥ O M O M O M ⎥ ⎥⎢ M ⎥ &⎥ L Z ii L Z ij L Z in ⎥ ⎢ I i ⎥⎢ ⎥ O M O M O M ⎥⎢ M ⎥ &⎥ L Z ji L Z jj L Z jn ⎥ ⎢ I ⎥⎢ j⎥ O M O M O M ⎥⎢ M ⎥ ⎢I &⎥ L Z ni L Z nj L Z nn ⎥ ⎥ ⎢ ⎦⎣ n⎥ ⎦
22
节点阻抗矩阵
Z =Y
−1
节点阻抗矩阵
Z ii
Z ij
对角元素,称为节点i的自阻抗 非对角元素称为节点i和节点j之间的互阻抗
⎡ Z11 ⎢ M ⎢ ⎢ Z i1 ⎢ ⎢ M ⎢ Z j1 ⎢ ⎢ M ⎢Z ⎣ n1

潮流计算的计算机方法

潮流计算的计算机方法

一、潮流计算的计算机方法对于复杂网络的潮流计算,一般必须借助电子计算机进行。

其计算步骤是:建立电力网络的数学模型,确定计算方法、制定框图和编制程序。

本章重点介绍前两部分,并着重阐述在电力系统潮流实际计算中常用的、基本的方法。

1,电力网络的数学模型电力网络的数学模型指的是将网络有关参数相变量及其相互关系归纳起来所组成的.可以反映网络性能的数学方程式组。

也就是对电力系统的运行状态、变量和网络参数之间相互关系的—种数学描述。

电力网络的数学模型有节点电压方程和回路电流方程等,前者在电力系统潮流计算中广泛采用。

节点电压方程又分为以节点导纳矩阵表示的节点电压方程和以节点阻抗矩阵表示的节点电压方程。

(1)节点导纳矩阵在电路理论课中。

已讲过了用节点导纳矩阵表示的节点电压方程:对于n个节点的网络其展开为:上式中,I是节点注入电流的列向量。

在电力系统计算中,节点注入电流可理解为节点电源电流与负荷电流之和,并规定电源向网络节点的注人电流为正。

那么,只有负荷节点的注入电流为负,而仅起联络作用的联络节点的注入电流为零。

U是节点电压的列向量。

网络中有接地支路时,通常以大地作参考点,节点电压就是各节点的对地电压。

并规定地节点的编号为0。

y是一个n×n阶节点导纳矩阵,其阶数n就等于网络中除参考节点外的节点数。

物理意义:节点i单位电压,其余节点接地,此时各节点向网络注入的电流就是节点i 的自导纳和其余节点的与节点i之间的互导纳。

特点:对称矩阵,稀疏矩阵,对角占优(2) 节点阻抗矩阵对导纳阵求逆,得:其中称为节点阻抗矩阵,是节点导纳矩阵的逆阵。

物理意义:节点i注入单位电流,其余节点不注入电流,此时各节点的电压就是节点i 的自阻抗和其余节点的与节点i之间的互阻抗。

特点:满阵,对称,对角占优2,功率方程、变量和节点分类(1)功率方程已知的是节点的注入功率,因此,需要重新列写方程: **==B B B B B U S I U Y其展开式为: i i i nj j ij U jQ P U Y ~1-=∑= 所以:∑=**=+nj jij i i i U Y U jQ P 1 展开写成极坐标方程的形式:)cos sin ()sin cos (11ij ij ij ij n j j i i ij ij ij ij n j j i i B G U U Q B G U U P δδδδ-=+=∑∑==所以节点的功率方程为:)cos sin ()sin cos (11ij ij ij ij n j j i di Gi i ij ij ij ij nj j i di Gi i B G U U Q Q Q B G U U P P P δδδδ---=∆+--=∆∑∑==(2) 变量分类负荷消耗的有功、无功功率取决于用户,因而是无法控制的,故称为不可控变量或扰动变量。

潮流计算的计算机算法资料

潮流计算的计算机算法资料

第四章潮流计算的计算机算法第一节概述潮流计算是电力系统最基本、最常用的计算。

根据系统给定的运行条件、网络接线及元件参数,通过潮流计算可以确定各母线的电压(幅值及相角),各元件中流过的功率、整个系统的功率损耗等。

潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。

因此潮流计算在电力系统的规划设计、生产运行、调度管理及科学研究中都有着广泛的应用。

电力系统潮流计算分为离线潮流计算和在线潮流计算。

前者主要用于系统规划设计和安排系统的运行方式,后者则用于正在运行系统的经常监视及实时控制。

本章主要讨论离线潮流计算问题,它的基本算法同样适用于在线潮流计算。

潮流计算在数学上是多元非线性方程组的求解问题,求解的方法有很多种。

自从五十年代计算机应用于电力系统以来,当时求解潮流的方法是以节点导纳矩阵为基础的逐次代入法(导纳法),后来为解决导纳法的收敛性较差的问题,出现了以阻抗矩阵为基础的逐次代入法(阻抗法)。

到六十年代,针对阻抗法占用计算机内存大的问题又出现了分块阻抗法及牛顿-拉夫逊(Newton-Raphson)法。

Newton —Raphson法是数学上解非线形方程式的有效方法,有较好的收敛性。

将N-R法用于潮流计算是以导纳矩阵为基础的,由于利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧,使N-R法在收敛性、占用内存、计算速度方面的优点都超过了阻抗法,成为六十年代末期以后普遍采用的方法。

同时国内外广泛研究了诸如非线形规划法、直流法、交流法等各种不同的潮流计算方法。

七十年代以来,又涌现出了更新的潮流计算方法。

其中有1974年由B、Stott、O、Alsac 提出的快速分解法以及1978年由岩本伸一等提出的保留非线性的高129速潮流计算法。

其中快速分解法(Fast decoupled load flow)从1975年开始已在国内使用,并习惯称之为PQ分解法。

由于PQ分解法在计算速度上大大超过N-R法,不但能应用于离线潮流计算,而且也能应用于在线潮流计算。

潮流计算的主要方法

潮流计算的主要方法

潮流计算的主要方法
最近几年,随着计算机仿真技术和复杂系统全面发展,潮流计算也受到越来越多的重视。

潮流计算是研究不同电力网络的物理特性和操作规律的一项重要工作。

针对潮流计算的主要方法,总结如下:
一、基于动力学的方法
1. 碰撞模型:根据动力学方法,计算电力系统的运行稳定性。

基于动力学的碰撞模型能够快速而精确地预测两个潮流的变化情况。

2. 时变快速收敛:在碰撞模型的基础上,为快速求解电力系统潮流,提出了时变快速收敛算法。

可以更快地获得潮流解。

二、基于牛顿迭代法的方法
1.牛顿迭代潮流计算方法:根据牛顿迭代法,采用迭代算法,求解电力系统潮流运行状态。

2. 功率流计算方法:计算机基于牛顿迭代法,快速求解节点电能的功率流公式。

可以有效的缩短潮流计算的时间,提高计算效率。

三、基于模糊聚类算法的方法
1. 基于模糊聚类的潮流计算方法:采用模糊聚类算法,对潮流计算进行多维度分析,可以得出最优的潮流结果。

2. 基于模糊划分的多目标模糊控制:根据模糊聚类理论,对潮流算法进行最佳控制,以满足电力网不同优化目标。

四、基于期望最大化的方法
1、基于粒子群优化的潮流计算方法:采用粒子群优化算法,将电力网潮流计算定义为多目标最优化问题,以期望最大化来求解潮流值,提高计算效率。

2、基于遗传算法的潮流计算方法:遗传算法利用进化过程来搜索全局最优解,使用遗传变异原则来改变候选解,以期望最大化来求解潮流计算问题。

潮流计算的计算机算法

潮流计算的计算机算法

高等电力系统分析(潮流计算的计算机算法)PQ分解法潮流计算(IEEE14)目录一、MATLAB源程序二、对支路参数(B1)、节点参数(B2)的说明三、带入数据,运行结果非对角元 对角元MATLAB 源程序 clearclose alln=in put(' 请输入节点数:n=');n 仁in put(' 请输入支路数:n 1=') )isb=in put(' 请输入平衡节点号 :isb=')pr=i nput(' 请输入误差精度:pr=');B1=i nput(' 请输入支路参数:B1 =');B2=i nput(' 请输入节点参数:B2=');n2=in put(' 请输入PQ 节点个数: :n2=');Y=zeros (n); for i=1: n1P=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; enddisp(' 导纳矩阵Y=');disp(Y)% --------------------------------------------% -------------- 下面是求P,Q,V,0矩阵 -----------------------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: nP(i)=B2(i,3);Q(i)=B2(i,4);V(i)=B2(i,5);O(i)=B2(i,6);endB3=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~=0m1=1;m2=1;for i=1:n所有节点数 节点数 求 DP,DQ 对V 矩阵求逆 △ P/V △ P/V/B3 △角=-△ P/V/V/B3△ V=- △ Q/V/Bif i~=isbC(i)=O;D(i)=O;for j1=1:n C(i)=C(i)+V(i)*V(j1)*(G(i,j1)*cos(O(i)-O(j1))+B(i,j1)*si n( O(i)-O(j1)));D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*si n(O(i)-O(j1))-B(i,j1)*cos(O(i)-O( j1)));endDP(m1)=P(i)-C(i); m1=m1+1;if B2(i,2)==1 DQ(m2)=Q(i)-D(i); m2=m2+1;endendendm1=m1-1;% m2=m2-1; %PQ DPQ=[DP';DQ']; % V 仁 V(:,1:m1);V2=diag(V1);V3=i nv(V2);% H=V3*DP'; %K=-i nv(B3)*H; %- deltO=V3*K;% max1=max(abs(DP));for i=1:m1if max1<prMp=0;elseO(i)=O(i)+deltO(i)'; Mq=1;endendV4=V(:,1:m2);V5=diag(V4);V6=i nv(V5);L=V6*DQ';N=-i nv(B4)*L; deltV=N; %max2=max(abs(DQ));for i=1:m2if max2<prMq=0;elseif B2(i,2)==1;V(i)=V(i)+deltV(i):Mp=1;endendendif Mp==0&&Mq==0ICT=0;elseICT=1;endDX=DX+1;end% ------------------------------------------% --------------- 迭代结束,开始输出结果disp(' --------------------------------------------');disp(' 迭代次数为:');disp(DX);for i=1: nE(i)=V(i)*cos(O(i))+1j*V(i)*si n(O(i));o(i)= 180*a ngle(E(i))/pi;enddisp(' --------------------------------------------');disp(' 修正后各节点电压标么值为(节点号从小到大排列) :');disp(V);disp(' --------------------------------------------');disp(' 修正后各节点电压相角为(节点号从小到大排列) :');disp(o);% ----------- 计算各个节点的功率disp(' --------------------------------------------');disp(' 各节点的功率为:');for p=1: nC(p)=0;for q=1: nC(p)=C(p)+conj(Y(p,q)*conj(E(q)));endS(p)=E(p)*C(p);enddisp(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: nSp=Sp+V( n)*conj(Y( n, i))*conj(V(i));enddisp(' -------------------------------------------- ');disp(' 平衡节点功率为:');disp(Sp);二、对支路参数(B1)、节点参数(B2)的说明1•节点数:142•支路数:203•支路矩阵B1的各支路参数:起点编号,终点编号,电阻,电抗,电导,电纳[1 2 0.01335 0.04211 0 0 ;1 3 0 0.20912 0 0 J1 4 0 0.55618 0 0 J1 10 0.05811 0.17632 0 0.0341 11 0.06701 0.17103 0 0.01282 10 0.05695 0.17388 0 0.0346 2 12 0 0.25202 0 0 J2 14 0.05403 0.22304 0 0.04923 4 0 0.11001 0 0 J3 13 0 0.17615 0 0 J4 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 07 8 0.22092 0.19988 0 0 ;[1 1 -0.4782 1 -0.0763 10 04 1 -0.2955 1 -0.096 1 -0.0357 1 -0.0618 1 -0.1359 1 -0.14910 2 0.18311 2 -0.94212 2 -0.11213 214 00.039 1 0;-0.016 1 0;1 0;-0.166 1 0;-0.058 1 0;-0.018 1 0;-0.016 1 0;-0.058 1 0;-0.05 1 0;0 1.045 0;0 1.01 0;0.047 1.70;0 0.174 1.9 0;0 0 1.06 0;]7 12 0.12291 0.25581 0 08 9 0.17093 0.34802 0 0 ;8 12 0.06615 0.13027 0 0 ;10 11 0.04699 0.19797 0 0.043810 14 0.01938 0.05917 0 0.05284.节点参数矩阵B2的各节点参数:(对应的每一列为)节点编号,类型,注入有功,注入无功,电压幅值,电压相位其中节点类型:仁PQ节点,2=PV节点,0二平衡节点三、带入数据,运行结果>> clearclose alln=input('请输入节点数:n=');n1=input('请输入支路数:n仁');isb=i nput('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr=');B1= input('请输入支路参数:B仁');B2=input('请输入节点参数:B2=');n2=input('请输入PQ节点个数:n2=');Y=zeros (n);for i=1: n1p=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;enddisp('导纳矩阵Y=');disp(Y)% --------------------------------------------% -------------- 下面是求PQ,V,0矩阵-------------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: nP(i)=B2(i,3);Q(i)=B2(i,4);V(i)=B2(i,5);O(i)=B2(i,6);endB3=B(1: n-1,1: n-1); %不含平衡节点,由节点导纳虚部构成B4=B(1:n2,1:n2); %所有PQ 节点% --------------------------------------------endendend m1=m1-1;m2=m2-1; DPQ=[DP';DQ']; V1=V(:,1:m1); V2=diag(V1);V3=i nv(V2);H=V3*DP';K=-i nv(B3)*H; deltO=V3*K;% 所有节点数%PQ节点数%求DPQQ%对V矩阵求逆%A P/V%-A P/V/B3%A 角=-A% -------------- 下面是求△ P,A Q矩阵------------DX=0;ICT=1;Mp=1;Mq=1;while ICT~=0m1=1;m2=1;for i=1: nif i~=isbC(i)=0;D(i)=0;for j1=1: nC(i)=C(i)+V(i)*V(j1)*(G(i,j1)*cos(O(i)-O(j1))+B(i,j1)*si n(O(i)-O(j1)));D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*si n(O(i)-O(j1))-B(i,j1)*cos(O(i)-O(j1)));endDP(m1)=P(i)-C(i); m1=m1+1;if B2(i,2)==1DQ(m2)=Q(i)-D(i); m2=m2+1;max1=max(abs(DP));for i=1:m1if max1<prMp=0;elseO(i)=O(i)+deltO(i)';Mq=1;endendV4=V(:,1:m2);V5=diag(V4); V6=i nv(V5);L=V6*DQ';N=-i nv(B4)*L;deltV=N; %△V=-A Q/V/Bmax2=max(abs(DQ));for i=1:m2if max2<prMq=0;elseif B2(i,2)==1;V(i)=V(i)+deltV(i): Mp=1;endendendif Mp==0&&M q==0ICT=0;elseICT=1;endDX=DX+1;end% ------------------------------------------% --------------- 迭代结束,开始输出结果disp(' --------------------------------------------');disp('迭代次数为:');disp(DX);for i=1: nE(i)=V(i)*cos(O(i))+1j*V(i)*si n(0 (i));o(i)= 180*a ngle(E(i))/pi;enddisp(' --------------------------------------------');disp('修正后各节点电压标么值为(节点号从小到大排列) :');disp(V);disp(' --------------------------------------------');disp('修正后各节点电压相角为(节点号从小到大排列) :');disp(o);% ----------- 计算各个节点的功率-----------------disp(' --------------------------------------------');disp('各节点的功率为:');for p=1: nC(p)=0;for q=1: nC(p)=C(p)+conj(Y(p,q)*conj(E(q)));endS(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: nSp=Sp+V( n)*conj(Y( n,i))*conj(V(i));enddisp(' -------------------------------------------- ');disp('平衡节点功率为:');disp(Sp);请输入节点数:n=14请输入支路数:n 1=20请输入平衡节点号:isb=14请输入误差精度:pr=0.00001请输入支路参数:B1=[1 2 0.01335 0.04211 0 0 ;1 3 0 0.20912 0 0 J1 4 0 0.55618 0 0 J1 10 0.05811 0.17632 0 0.0341 11 0.06701 0.17103 0 0.01282 10 0.05695 0.17388 0 0.03462 12 0 0.25202 0 0 J2 14 0.05403 0.22304 0 0.049210.5130 -38.2963i -6.8410 +21.5786i 0.0000 + 4.7819i0.0000 + 1.7980i0.0000 +0.0000i-6.8410 +21.5786i 9.5680 -34.8916i0.0000 + 0.0000i 0.0000 + 0.0000i0.0000 +0.0000i0.0000 + 4.7819i 0.0000 + 0.0000i0.0000 -19.5490i 0.0000 + 9.0901i0.0000 +0.0000i0.0000 + 1.7980i0.0000 + 0.0000i0.0000 + 9.0901i5.3261 -24.2825i-3.9020 +10.3654i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-3.9020 +10.3654i 5.7829-14.7683i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i 0.0000 + 0.0000i-1.8809 +4.4029i0.0000 + 0.0000i 0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 +O.OOOOi3 4 0 0.11001 00 ; 3 13 0 0.17615 00 J 4 5 0.03181 0.08450 ; 4 9 0.12711 0.27038 0 0 ; 5 60.08205 0.19207 00 ; 6 12 0.09498 0.1989 0 0 ; 7 8 0.22092 0.19988 0 0 ; 7 12 0.12291 0.25581 00 ; 8 9 0.17093 0.34802 0 0 ; 8 12 0.06615 0.13027 0 0 ;10 11 0.04699 0.19797 00.0438 ; 10 14 0.01938 0.05917 0 0.0528;请输入节点参数:B2=[1 1 -0.478 0.039 2 1 -0.076 -0.016 10;3 1 0 0 1 0;4 1 -0.295 -0.166 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 10; 10 2 0.183 01.0450;11 2 -0.9420 1.01 0;12 2 -0.112 0.047 1.7 0;13 20 0.1741.9 0;14 0 0 0 1.06 0;]请输入PQ 节点个数:n2=9导纳矩阵 Y=1Columns 1 through 5;0.0000 +O.OOOOiO.OOOOi0.0000 + 0.0000i 0.0000i-1.6860 + 5.1158i 0.0000i-1.9860 + 5.0688i 0.0000i0.0000 + 0.0000i 0.0000i0.0000 + 0.0000i 0.0000i0.0000 + 0.0000i 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.7011 + 5.1939i0.0000 + 0.0000i0.0000 + 3.9679i0.0000 + 0.0000i-1.0259 + 4.2350i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 5.6770i0.0000 + 0.0000i0.0000 + 0.0000i-1.4240 + 3.0291i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000iColu mns 6 through 100.0000 + 0.0000i 5.115 8i0.0000 + 0.0000i 5.1939i0.0000 + 0.0000i 0.0000i0.0000 + 0.0000i 0.0000i-1.8809 + 4.4029i 0.0000i3.8359 - 8.4970i 0.0000i0.0000 + 0.0000i 0.0000i0.0000 + 0.0000i 0.0000i0.0000 + 0.0000i 0.0000i0.0000 + 0.0000i -30.1895i0.0000 + 0.0000i 4.7819i-1.9550 + 4.0941i 0.0000i0.0000 + 0.0000i 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i4.0150 -5.4279i-2.4890 + 2.2520i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.5260 + 3.1760i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-2.4890 + 2.2520i6.7249 -10.6697i-1.1370 + 2.3150i0.0000 + 0.0000i0.0000 + 0.0000i-3.0989 + 6.1028i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.4240 + 3.0291i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.1370 + 2.3150i2.5610 - 5.3440i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +-1.6860 +-1.7011 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +9.5213-1.1350 +0.0000 +0.0000 +0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -4.9991-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 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.1350 + 4.7819i 3.1210 - 9.7941i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 3.9679i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.9550 + 4.0941i-1.5260 + 3.1760i-3.0989 + 6.1028i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i6.5799 -17.3407i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 5.6770i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 - 5.6770i0.0000 + 0.0000i0.0000 + 0.0000i-1.0259 + 4.2350i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-4.9991 +15.2631i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i6.0250 -19.3961i+15.2631iColu mns 11 through 14迭代次数为:42修正后各节点电压标么值为Colu mns 1 through 10(节点号从小到大排列)1.2128 1.21481.5809 1.04501.5627 1.5364 1.5602 1.6264 1.6792 1.6654Colu mns 11 through 141.0100 1.7000 1.90001.0600修正后各节点电压相角为(节点号从小到大排列)Colu mns 1 through 10-13.1377 -11.9998 -14.4540 -15.0006 -15.3011 -15.7590 -16.4954 -16.3362 -16.0220 -5.2081Colu mns 11 through 14-12.1568 -16.1918 -14.4540-0.0475 --1.5950 -各节点的功率为:Columns 1 through 5 1.4718 - 1.8800i 0.9041 - 1.1423i0.0000 - O.OOOOi -0.1314 + 0.1399i0.0646iColu mns 6 through 10 -0.0300 + 0.0035i 0.0039 - 0.1228i -0.0656 - 0.0724i-0.0352 - 0.1437i1.3905iColu mns 11 through 14 -0.4870 - 1.2396i3.2005 + 3.0414i1.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 +16.2405i各条支路的功率损耗为23.5109 +78.7335i各条支路的首端功率为17.0452 +57.7870i各条支路的末端功率为11.7622 +39.7706i各条支路的功率损耗为28.8075 +97.5576i各条支路的首端功率为13.5220 +59.6677i19.6150 +58.2875i各条支路的功率损耗为:3.3137e+01 + 1.1796e+02i各条支路的首端功率为16.5471 +56.5546i各条支路的末端功率为6.9281 +27.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+02i21.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.3933i39.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各条支路的首端功率为27.4696 +46.8495i各条支路的末端功率为27.7461 +67.4149i各条支路的功率损耗为:5.5216e+01 + 1.1426e+02i各条支路的首端功率为10.9761 +38.1226i各条支路的末端功率为4.9835 +14.8560i各条支路的功率损耗为15.9596 +52.9785i各条支路的首端功率为17.4469 +49.3022i各条支路的末端功率为10.7497 +39.1332i各条支路的功率损耗为28.1966 +88.4354i平衡节点功率为:-0.0889 - 0.5670i。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

学号1350803233《工厂供电》课程设计(2013级本科)系(部)院:物理与机电工程学院专业:电气工程及其自动化作者姓名:杨兴海指导教师:田娜职称:助教完成日期:2016 年 1 月 5 日目录摘要: (1)第一章电力系统潮流计算概述 (2)1.1 潮流计算简介 (2)1.2 潮流计算的意义及其发展 (2)第二章电气参数的确定 (3)2.1变压器、线路参数确定 (3)2.2 电力线路、变压器模型的建立 (4)第三章 PQ分解方法潮流计算分析 (5)3.1 PQ分解法的极坐标表示及简化算法 (5)3.1.1潮流计算的定义 (5)3.1.2潮流计算的约束条件 (6)3.1.3节点电压用极坐标表示时的牛顿-拉夫逊潮流计算 (6)3.1.4对牛顿—拉夫逊法潮流计算的数学模型进行简化修正 (8)3.2 PQ分解法潮流计算的简化算法 (10)3.3 PQ 分解法潮流计算的基本步骤 (13)第四章程序编写及结果分析 (13)4.1程序编程 (13)4.2执行结果 (20)第五章课程设计心得 (22)参考文献 (22)摘要潮流计算是电力系统最基本最常用的计算。

根据系统给定的运行条件,网络接线及元件参数,通过潮流计算可以确定各母线的电压(幅值和相角),各支路流过的功率,整个系统的功率损耗。

潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。

因此,潮流计算在电力系统的规划计算,生产运行,调度管理及科学计算中都有着广泛的应用。

潮流计算在数学上是多元非线性方程组的求解问题,PQ分解法是数学上解非线性方程组的有效方法,有较好的收敛性。

运用电子计算机计算一般要完成以下几个步骤:建立数学模型,确定解算方法,制订计算流程,编制计算程序。

关键字:PQ分解法计算机潮流计算 MATLAB第一章电力系统潮流计算1.1 潮流计算简介电力系统潮流计算是研究电力系统稳态运行情况的一种计算,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态:各母线的电压,各元件中流过的功率,系统的功率损耗等等。

在电力系统规划的设计和现有电力系统运行方式的研究中,都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。

可靠性和经济性。

此外,电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。

所以潮流计算是研究电力系统的一种很重要和基础的计算。

电力系统潮流计算也分为离线计算和在线计算两种,前者主要用于系统规划设计和安排系统的运行方式,后者则用于正在运行系统的经常监视及实时控制。

利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。

在这20年内,潮流计算曾采用了各种不同的方法,这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。

对潮流计算的要求可以归纳为下面几点:(1)计算方法的可靠性或收敛性;(2)对计算机内存量的要求;(3)计算速度;(4)计算的方便性和灵活性。

电力系统潮流计算问题在数学上是一组多元非线性方程式求解问题,其解法都离不开迭代。

因此,对潮流计算方法,首先要求它能可靠地收敛,并给出正确答案。

由于电力系统结构及参数的一些特点,并且随着电力系统不断扩大,潮流计算的方程式阶数也越来越高,对这样的方程式并不是任何数学方法都能保证给出正确答案的。

这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。

1.2 潮流计算的意义及其发展(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。

(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。

(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。

(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。

总结为在电力系统运行方式和规划方案的研究中,都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。

同时,为了实时监控电力系统的运行状态,也需要进行大量而快速的潮流计算。

因此,潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。

在系统规划设计和安排系统的运行方式时,采用离线潮流计算;在电力系统运行状态的实时监控中则采用在线潮流计算。

近20多年来,潮流算法的研究仍然非常活跃,但是大多数研究都是围绕改进牛顿法和P-Q分解法进行的。

此外,随着人工智能理论的发展遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。

但是,到目前为止这些新的模型和算法还不能取代牛顿法和P-Q分解法的地位。

由于电力系统规模的不断扩大,对计算速度的要求不断提高,计算机的并行计算技术也将在潮流计算中得到广泛的应用,成为重要的研究领域。

第二章电气参数的确定2.1变压器、线路参数确定查设计手册:LGJQ4×400型导线:=0.020Ω/km,=0.276Ω/km,=4.02S/km; LGJ240型导线:=0.131Ω/km,=0.432Ω/km; LGJ185型导线:=0.170Ω/km,=0.440Ω/km。

变压器型号:T1,T2为SF-360000/220;T3部分参数为:额定容量240/120/120MVA,额定电压500/242/38.5kv。

T4部分参数为:额定容量120/120/60MVA, ,额定电压242/121/10.5kv,等值电路中所有参数都归算到高压侧。

线路L1,L2的电导,L4,L5的导纳都可略去。

表2-1短路电压(未经归算1.采用有名值:2.采用标幺值:取基准功率为1000MVA,基准电压为500kv2.2.2电力变压器及线路参数的计算结果电力线路参数见表2-1,电力变压器参数见表2-2表2-1电力线路参数表2-2电力变压器参数2.2 电力线路、变压器模型的建立线路采用Π型等值电路,变压器采用Γ型等值电路,有名值和标幺值等值电路分别见:图2-1、图2-2。

图2-1 有名值等值网络图2-2标幺值等值网络第三章PQ分解方法潮流计算分析3.1 PQ分解法的极坐标表示及简化算法3.1.1潮流计算的定义潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件,确定电力系统各部分稳态运行状态参数的计算。

通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。

待求的运行状态参量包括电网各母线节点的电压幅值和相角,以及各支路的功率分布、网络的功率损耗等。

PQ分解法的极坐标表示是派生于以极坐标表示的牛顿—拉夫逊法。

3.1.2潮流计算的约束条件为了保证电力系统的正常运行潮流问题中某些变量应满足一定的约束条件,常用的约束条件有:(1)所有节点电压必须满足(i=1,2,…,n)从保证电能质量和供电安全的要求看,电力系统的所有电气设备都必须运行在额定电压附近。

PV节点的电压幅值必须按上述条件给定。

因此,这一约束主要是对PQ节点而言。

(2)所有电源节点的有功功率和无功功率必须满足PQ节点的有功功率和无功功率以及PV节点的有功功率,在给定时就必须满足上述条件。

因此,对平衡节点的P和Q以及PV节点的Q应按上述条件进行检验。

(1)某些节点之间电压的相位差应满足为了保证系统运行的稳定性,要求某些输电线路两端的电压相位差不超过一定的数值。

因此,潮流计算可以归纳为求解一组非线性方程组,并使其解答满足一定的约束条件。

如果不满足,则应修改某些变量的给定值,甚至修改系统的运行方式,重新进行计算。

3.1.3节点电压用极坐标表示时的牛顿-拉夫逊潮流计算采用极坐标时,节点电压表示为)(i i i i i jsin cos V i V V δδδ+=∠= 节点功率方程(11-25)将写成⎪⎪⎭⎪⎪⎬⎫=+=∑∑==n1j ij ij ij ij j i i ij ij ij ij n1j j i i cos B -sin G V V Q sin B cos G V V P δδδδ 公式 1 式中,j i ij δδδ-=,是i ,j 两节点电压的相差角。

方程式(公式1)把节点功率表示为节点电压的幅值和相角的函数。

在有n 节点的系统中,假定第1~m 号节点PQ 节点,第m+1~n-1号节点的PV 节点,第n 号节点为平衡节点。

n V 和n δ是给定的,PV 节点的电压幅值Vm+1~Vn-1也是给定的。

因此,只剩下n-1个节点的电压相角1-n 21δδδ,,, 和m 个节点的电压幅值1V ,2V ,…,m V 是未知量。

实际上,对于每一个PQ 节点或每一个PV 节点都可以列写一个有功功率不平衡量方程式i P ∆=is P -i P =is P -iV 0)sin cos (GV ij ij ij ijn1j j =+∑=δδB (i=1,2,…,n -1)公式 2而对于每一个PQ 节点还可以再列写一个无功功率不平衡量方程式i Q ∆=is Q -i Q =is Q -iV 0)cos sin (ij ij ij ij1j =-∑=δδB GV nj (i=1,2,…,m )公式 3式2和式3一共包含了n-1+m 个方程式,正好同未知量的数目相等,而比直角坐标形式的方程式少了n-1-m 个。

对于方程式2和式3可以写出修正方程式如下:⎥⎥⎦⎤⎢⎢⎣⎡∆∆⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡∆∆-V L KN HQ P V D 12δ公式 4式中⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡∆∆∆=∆-P P P n P 121 ; ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡∆∆∆=∆1-n Q Q Q Q M ; ⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡∆∆∆=∆-δδδδ121n ⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡∆∆∆=∆-V V V n V 121 ; ⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=V V V V m D O 212 公式 5H 是(n-1)×(n-1)阶方阵,其元素为jiij P H δ∂∆∂=;N 是(n-1)×m 阶矩阵,其元素为jij ij V P V N ∂∆∂=;K 是m ×(n -1)阶矩阵,其元素为j i ijQ K δ∂∆∂=;L 是m × m 阶方阵,其元素为jijij V Q V L ∂∆∂=。

3.1.4对牛顿—拉夫逊法潮流计算的数学模型进行简化修正在交流高压电网中,输电线路的电抗要比电阻大得多,系统中母线有功功率的变化主要受电压相位的影响,无功功率的变化则主要受母线电压幅值变化的影响。

在修正方程式的系数矩阵中,偏导数V P ∂∆∂和δ∂∆∂Q 的数值相对于偏导数δ∂∆∂P 和V Q∂∆∂是相当小的。

相关文档
最新文档