实验一潮流计算

学院：电气工程学院专业班级：电气F1104

姓名：马路林学号：201123910725

潮流计算实验运行程序

>> %本程序的功能是用牛顿——拉夫逊法进行潮流计算

% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳

% 5、支路的变比；6、支路首端处于K侧为1，1侧为0

% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值

% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量

% 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；

% 3为PV节点；

clear;

n=10;

nl=10;

isb=1;

pr=0.00001;

B1=[1 2 0.03512+0.08306i 0.13455i 1 0;

2 3 0.0068+0.18375i 0 1.02381 1;

1 4 0.05620+0.13289i 0.05382i 1 0;

4 5 0.00811+0.24549i 0 1.02381 1;

1 6 0.05620+0.13289i 0.05382i 1 0;

4 6 0.04215+0.09967i 0.04037i 1 0;

6 7 0.0068+0.18375i 0 1.02381 1;

6 8 0.02810+0.06645i 0.10764i 1 0;

8 10 0.00811+0.24549i 0 1 1;

8 9 0.03512+0.08306i 0.13455i 1 0]

B2=[0 0 1.1 1.1 0 1;

0 0 1 0 0 2;

0 0.343+0.21256i 1 0 0 2;

0 0 1 0 0 2;

0 0.204+0.12638i 1 0 0 2;

0 0 1 0 0 2;

0 0.306+0.18962i 1 0 0 2;

0 0 1 0 0 2;

0.5 0 1.1 1.1 0 3;

0 0.343+0.21256i 1 0 0 2]

Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);

for i=1:nl

if B1(i,6)==0

p=B1(i,1);q=B1(i,2);

else

p=B1(i,2);q=B1(i,1);

end

Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));

Y(q,p)=Y(p,q);

Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;

Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;

end

disp

disp(Y)

%----------------------------------------------------------

G=real(Y);B=imag(Y);

for i=1:n

e(i)=real(B2(i,3));

f(i)=imag(B2(i,3));

V(i)=B2(i,4);

end

for i=1:n

S(i)=B2(i,1)-B2(i,2);

B(i,i)=B(i,i)+B2(i,5);

end

P=real(S);Q=imag(S);

ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;

while IT2~=0

IT2=0;a=a+1;

for i=1:n

if i~=isb

C(i)=0;D(i)=0;

for j1=1:n

C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)

D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bij*ej) end

P1=C(i)*e(i)+f(i)*D(i);

Q1=C(i)*f(i)-e(i)*D(i);

V2=e(i)^2+f(i)^2;

if B2(i,6)~=3

DP=P(i)-P1;

DQ=Q(i)-Q1;

for j1=1:n

if j1~=isb&j1~=i

X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de=-dQ/df

X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df=dQ/de

X3=X2; % X2=dp/df X3=dQ/de

X4=-X1; % X1=dP/de X4=dQ/df

p=2*i-1;q=2*j1-1;

J(p,q)=X3;J(p,N)=DQ;m=p+1;

J(m,q)=X1;J(m,N)=DP;q=q+1;

J(p,q)=X4;J(m,q)=X2;

elseif j1==i&j1~=isb

X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);

X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);

X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);

X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);

p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;

m=p+1;

J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;

J(m,q)=X2;

end

end

else

DP=P(i)-P1;

DV=V(i)^2-V2;

for j1=1:n

if j1~=isb&j1~=i

X1=-G(i,j1)*e(i)-B(i,j1)*f(i);

X2=B(i,j1)*e(i)-G(i,j1)*f(i);

X5=0;X6=0;

p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;

m=p+1;

J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;

J(m,q)=X2;

elseif j1==i&j1~=isb

X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de

X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df

X5=-2*e(i);

X6=-2*f(i);

p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;

m=p+1;

J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;

J(m,q)=X2;

end

end

end

end

end

for k=3:N0

k1=k+1;N1=N;

for k2=k1:N1

J(k,k2)=J(k,k2)./J(k,k);

end

J(k,k)=1;

if k~=3

k4=k-1;

for k3=3:k4

for k2=k1:N1

J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

if k==N0

break;

end

for k3=k1:N0

for k2=k1:N1

J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

else

for k3=k1:N0

for k2=k1:N1

J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

end

end

for k=3:2:N0-1

L=(k+1)./2;

e(L)=e(L)-J(k,N);

k1=k+1;

f(L)=f(L)-J(k1,N);

end

for k=3:N0

DET=abs(J(k,N));

if DET>=pr

IT2=IT2+1;

end

end

ICT2(a)=IT2;

ICT1=ICT1+1;

end

disp('迭代次数：');

disp(ICT1);

disp('没有达到精度要求的个数：');

disp(ICT2);

for k=1:n

V(k)=sqrt(e(k)^2+f(k)^2);

sida(k)=atan(f(k)./e(k))*180./pi;

E(k)=e(k)+f(k)*j;

end

disp('各节点的实际电压标幺值E为(节点号从小到大排列)：'); disp(E);

disp('-----------------------------------------------------');

disp('各节点的电压大小V为(节点号从小到大排列)：');

disp(V);

disp('-----------------------------------------------------');

disp('各节点的电压相角sida为(节点号从小到大排列)：'); disp(sida);

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为(节点号从小到大排列)：');

disp(S);

disp('-----------------------------------------------------');

disp('各条支路的首端功率Si为(顺序同您输入B1时一致)：');

for i=1:nl

p=B1(i,1);q=B1(i,2);

if B1(i,6)==0

Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))...

-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));

Siz(i)=Si(p,q);

else

Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))...

-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));

Siz(i)=Si(p,q);

end

disp(Si(p,q));

SSi(p,q)=Si(p,q);

ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];

disp(ZF);

disp('-----------------------------------------------------');

end

disp('各条支路的末端功率Sj为(顺序同您输入B1时一致)：');

for i=1:nl

p=B1(i,1);q=B1(i,2);

if B1(i,6)==0

Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))...

-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));

Sjy(i)=Sj(q,p);

else

Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))...

-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));

Sjy(i)=Sj(q,p);

end

disp(Sj(q,p));

SSj(q,p)=Sj(q,p);

ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];

disp(ZF);

disp('-----------------------------------------------------');

end

disp('各条支路的功率损耗DS为(顺序同您输入B1时一致)：');

for i=1:nl

p=B1(i,1);q=B1(i,2);

DS(i)=Si(p,q)+Sj(q,p);

disp(DS(i));

DDS(i)=DS(i);

ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];

disp(ZF);

disp('-----------------------------------------------------'); end

figure(1);

subplot(1,2,1);

plot(V);

xlabel('节点号');ylabel('电压标幺值');

grid on;

subplot(1,2,2);

plot(sida);

xlabel('节点号');ylabel('电压角度');

grid on;

figure(2);

subplot(2,2,1);

P=real(S);Q=imag(S);

bar(P);

xlabel('节点号');ylabel('节点注入有功');

grid on;

subplot(2,2,2);

bar(Q);

xlabel('节点号');ylabel('节点注入无功');

grid on;

subplot(2,2,3);

P1=real(Siz);Q1=imag(Siz);

bar(P1);

xlabel('支路号');ylabel('支路首端注入有功');

grid on;

subplot(2,2,4);

bar(Q1);

xlabel('支路号');ylabel('支路首端注入无功');

grid on;

程序运行结果

B1 =

1.0000

2.0000 0.0351 + 0.0831i 0 + 0.1346i

1.0000 0

2.0000

3.0000 0.0068 + 0.1838i 0 1.0238 1.0000

1.0000 4.0000 0.0562 + 0.1329i 0 + 0.0538i 1.0000 0

4.0000

5.0000 0.0081 + 0.2455i 0 1.0238 1.0000

1.0000 6.0000 0.0562 + 0.1329i 0 + 0.0538i 1.0000 0

4.0000 6.0000 0.0422 + 0.0997i 0 + 0.0404i 1.0000 0

6.0000

7.0000 0.0068 + 0.1838i 0 1.0238 1.0000

6.0000 8.0000 0.0281 + 0.0664i 0 + 0.1076i 1.0000 0

8.0000 10.0000 0.0081 + 0.2455i 0 1.0000 1.0000

8.0000 9.0000 0.0351 + 0.0831i 0 + 0.1346i 1.0000 0

B2 =

0 0 1.1000 1.1000 0 1.0000

0 0 1.0000 0 0 2.0000

0 0.3430 + 0.2126i 1.0000 0 0 2.0000

0 0 1.0000 0 0 2.0000

0 0.2040 + 0.1264i 1.0000 0 0 2.0000

0 0 1.0000 0 0 2.0000

0 0.3060 + 0.1896i 1.0000 0 0 2.0000

0 0 1.0000 0

0 2.0000

0.5000 0 1.1000 1.1000 0 3.0000

0 0.3430 + 0.2126i 1.0000 0 0 2.0000

导纳矩阵Y=

9.7177 -22.8591i -4.3185 +10.2135i 0 -2.6996 + 6.3834i 0 -2.6996 + 6.3834i 0 0 0 0

-4.3185 +10.2135i 4.5104 -15.3311i -0.1964 + 5.3083i 0 0 0 0 0 0 0

0 -0.1964 + 5.3083i 0.2011 - 5.4347i 0 0 0 0 0 0 0

-2.6996 + 6.3834i 0 0 6.4271 -18.7292i -0.1313 + 3.9744i -3.5993 + 8.5110i 0 0 0 0

0 0 0 -0.1313 + 3.9744i 0.1344 - 4.0690i 0 0 0 0 0

-2.6996 + 6.3834i 0 0 -3.5993 + 8.5110i 0 11.8891 -32.7444i -0.1964 + 5.3083i -5.3984 +12.7660i 0 0

0 0 0 0 0 -0.1964 + 5.3083i 0.2011 - 5.4347i 0 0 0

0 0 0 0 0 -5.3984 +12.7660i 0 9.8514 -26.9275i -4.3185 +10.2135i -0.1344 + 4.0690i

0 0 0 0 0 0 0 -4.3185 +10.2135i 4.3185 -10.1462i 0

0 0 0 0 0 0 0 -0.1344 + 4.0690i 0 0.1344 - 4.0690i

迭代次数：

4

没有达到精度要求的个数：

17 18 17 0

各节点的实际电压标幺值E为(节点号从小到大排列)：

1.1000 1.0757 - 0.0207i 1.0050 - 0.0780i 1.0772 - 0.0175i 1.0171 -

0.0631i 1.0762 - 0.0152i 1.0112 - 0.0666i 1.0778 - 0.0051i 1.0996 + 0.0304i

1.0177 - 0.0814i

-----------------------------------------------------

各节点的电压大小V为(节点号从小到大排列)：

1.1000 1.0759 1.0080 1.0773 1.0191 1.0763 1.0134 1.0778 1.1000 1.0209

-----------------------------------------------------

各节点的电压相角sida为(节点号从小到大排列)：

0 -1.1046 -4.4373 -0.9283 -3.5503 -0.8106 -3.7665 -0.2718

1.5822 -4.5707

各节点的功率S为(节点号从小到大排列)：

0.7165 + 0.2587i 0.0000 + 0.0000i -0.3430 - 0.2126i-0.0000 + 0.0000i -0.2040 - 0.1264i -0.0000 + 0.0000i -0.3060 - 0.1896i -0.0000 + 0.0000i 0.5000 + 0.0089i -0.3430 - 0.2126i

-----------------------------------------------------

各条支路的首端功率Si为(顺序同您输入B1时一致)：

0.3485 + 0.0932i

S(1,2)=0.3485+0.093157i

-----------------------------------------------------

0.3441 + 0.2420i

S(2,3)=0.34409+0.24201i

-----------------------------------------------------

0.1904 + 0.0760i

S(1,4)=0.19038+0.07599i

-----------------------------------------------------

0.2044 + 0.1400i

S(4,5)=0.20445+0.13999i

-----------------------------------------------------

0.1777 + 0.0895i

S(1,6)=0.17767+0.089525i

-----------------------------------------------------

-0.0163 - 0.0055i

S(4,6)=-0.016305-0.0054856i

-----------------------------------------------------

0.3069 + 0.2128i

S(6,7)=0.30686+0.21281i

-----------------------------------------------------

-0.1477 - 0.0234i

S(6,8)=-0.14767-0.02338i

-----------------------------------------------------

0.3443 + 0.2509i

S(8,10)=0.34427+0.25091i

-----------------------------------------------------

-0.4925 - 0.1508i

S(8,9)=-0.49251-0.15077i

-----------------------------------------------------

各条支路的末端功率Sj为(顺序同您输入B1时一致)：-0.3441 - 0.2420i

S(2,1)=-0.34409-0.24201i

-----------------------------------------------------

-0.3430 - 0.2126i

S(3,2)=-0.343-0.21256i

-----------------------------------------------------

-0.1881 - 0.1345i

S(4,1)=-0.18815-0.13451i

-----------------------------------------------------

-0.2040 - 0.1264i

S(5,4)=-0.204-0.12638i

-----------------------------------------------------

-0.1755 - 0.1482i

S(6,1)=-0.17551-0.14815i

-----------------------------------------------------

0.0163 - 0.0413i

S(6,4)=0.016326-0.041272i

-----------------------------------------------------

-0.3060 - 0.1896i

S(7,6)=-0.306-0.18962i

-----------------------------------------------------

0.1482 - 0.1001i

S(8,6)=0.14824-0.10014i

-----------------------------------------------------

-0.3430 - 0.2126i

S(10,8)=-0.343-0.21256i

-----------------------------------------------------

0.5000 + 0.0089i

S(9,8)=0.5+0.0089402i

-----------------------------------------------------

各条支路的功率损耗DS为(顺序同您输入B1时一致)：

0.0044 - 0.1488i

DS(1,2)=0.0044095-0.14885i

-----------------------------------------------------

0.0011 + 0.0294i

DS(2,3)=0.0010897+0.029445i

-----------------------------------------------------

0.0022 - 0.0585i

DS(1,4)=0.0022306-0.058518i

-----------------------------------------------------

0.0004 + 0.0136i

DS(4,5)=0.00044972+0.013613i

-----------------------------------------------------

0.0022 - 0.0586i

DS(1,6)=0.0021584-0.058629i

-----------------------------------------------------

0.0000 - 0.0468i

DS(4,6)=2.1344e-005-0.046758i

-----------------------------------------------------

0.0009 + 0.0232i

DS(6,7)=0.00085804+0.023186i

-----------------------------------------------------

0.0006 - 0.1235i

DS(6,8)=0.00056584-0.12352i

-----------------------------------------------------

0.0013 + 0.0384i

DS(8,10)=0.001267+0.038353i

-----------------------------------------------------

0.0075 - 0.1418i

DS(8,9)=0.0074931-0.14183i

-----------------------------------------------------

潮流计算(matlab)实例计算

潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新

的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图

电力系统运行方式及潮流分析实验报告

电力系统运行方式及潮 流分析实验报告 Company Document number：WUUT-WUUY-WBBGB-BWYTT-1982GT

电力系统第一次实验报告——电力系统运行方式及潮流分析实验

实验1 电力系统运行方式及潮流分析实验 一、实验目的 1、掌握电力系统主接线电路的建立方法 2、掌握辐射形网络的潮流计算方法； 3、比较计算机潮流计算与手算潮流的差异； 4、掌握不同运行方式下潮流分布的特点。 二、实验内容 1、辐射形网络的潮流计算； 2、不同运行方式下潮流分布的比较分析 三、实验方法和步骤 1．辐射形网络主接线系统的建立 输入参数（系统图如下）： G1：300+j180MVA（平衡节点） 变压器B1：Sn=360MVA，变比=18/121，Uk％=％，Pk=230KW，P0=150KW，I0/In=1％； 变压器B2、B3：Sn=15MVA，变比=110/11 KV，Uk％=％，Pk=128KW， P0=，I0/In=％； 负荷F1：20+j15MVA；负荷F2：28+j10MVA； 线路L1、L2：长度：80km，电阻：Ω/km，电抗：Ω/km，电纳：×10-6S/km。 辐射形网络主接线图 （1）在DDRTS中绘出辐射形网络主接线图如下所示： （2）设置各项设备参数： G1：300+j180MVA（平衡节点） 变压器B1：Sn=360MVA，变比=18/121，Uk％=％，Pk=230KW，P0=150KW，I0/In=1％；

变压器B2、B3：Sn=15MVA，变比=110/11 KV，Uk％=％，Pk=128KW， P0=，I0/In=％； 负荷F1：20+j15MVA；负荷F2：28+j10MVA； 线路L1、L2：长度：80km，电阻：Ω/km，电抗：Ω/km，电纳：×10-6S/km。2．辐射形网络的潮流计算 （1）调节发电机输出电压，使母线A的电压为115KV，运行DDRTS进行系统潮流计算，在监控图页上观察计算结果 项目DDRTS潮流计算结果 变压器B2输入功率+ 变压器B2输出功率+ 变压器B3输入功率+ 变压器B3输出功率+ 线路L1输入功率+ 线路L1输出功率+ 线路L2输入功率+ 线路L2输出功率+ （2）手算潮流： （3）计算比较误差分析 通过比较可以看出，手算结果与计算机仿真结果相差不大。产生误差原因：手算时是已知首端电压、末端功率的潮流计算，计算过程中要将输电线路对地电容吸收的功率以及变压器励磁回路吸收的功率归算到运算负荷中，并且在每一轮的潮流计算中都用上一轮的电压或功率的值（第一轮电压用额定电压）。 3．不同运行方式下潮流比较分析 （1）实验网络结构图如上。由线路上的断路器切换以下实验运行方式： ①双回线运行（L1、L2均投入运行） ②单回线运行（L1投入运行，L2退出）将断路器断开 对上述两种运行方式分别运行潮流计算功能，将潮流计算结果填入下表：

电力系统分析潮流计算例题

电力系统的潮流计算 西安交通大学自动化学院 2012.10 3．1 电网结构如图3—11所示，其额定电压为10KV 。已知各节点的负荷功率及参数： MVA j S )2.03.0(2 +=， MVA j S )3.05.0(3+=， MVA j S )15.02.0(4+= Ω+=)4.22.1(12j Z ，Ω+=)0.20.1(23j Z ，Ω+=)0.35.1(24j Z 试求电压和功率分布。 解：（1）先假设各节点电压均为额定电压，求线路始端功率。 0068.00034.0)21(103.05.0)(2 2223232232323j j jX R V Q P S N +=++=++=?0019.00009.0)35.1(10 15.02.0)(2 2 224242242424j j jX R V Q P S N +=++=++=?

则： 3068.05034.023323j S S S +=?+= 1519.02009.024424j S S S +=?+= 6587.00043.122423' 12 j S S S S +=++= 又 0346 .00173.0)4.22.1(106587.00043.1)(2 2 212122'12'1212j j jX R V Q P S N +=++=++=? 故： 6933.00216.112'1212 j S S S +=?+= （2） 再用已知的线路始端电压kV V 5.101 =及上述求得的线路始端功率 12 S ，求出线 路 各 点 电 压 。

kV V X Q R P V 2752.05 .104.26933.02.10216.1)(11212121212=?+?=+=? kV V V V 2248.101212=?-≈ kV V V V kV V X Q R P V 1508.100740.0) (24242 2424242424=?-≈?=+=? kV V V V kV V X Q R P V 1156.101092.0) (23232 2323232323=?-≈?=+=? （3）根据上述求得的线路各点电压，重新计算各线路的功率损耗和线路始端功率。 0066.00033.0)21(12.103.05.02 2 223j j S +=++=? 0018.00009.0)35.1(15 .1015.02.02 2 224j j S +=++=? 故 3066.05033.023323j S S S +=?+= 1518.02009.024424j S S S +=?+= 则 6584.00042.122423' 12 j S S S S +=++= 又 0331.00166.0)4.22.1(22 .106584.00042.12 2 212j j S +=++=? 从而可得线路始端功率 6915.00208.112 j S +=

武汉大学电力系统分析实验报告

电气工程学院 《电力系统分析综合实验》2017年度PSASP实验报告 学号： 姓名： 班级：

实验目的： 通过电力系统分析的课程学习，我们都对简单电力系统的正常和故障运行状态有了大致的了解。但电力系统结构较为复杂，对电力系统极性分析计算量大，如果手工计算，将花费 大量的时间和精力，且容易发生错误。而通过使用电力系统分析程序PSASP，我们能对电 力系统潮流以及故障状态进行快速、准确的分析和计算。在实验过程中，我们能够加深对电力系统分析的了解，并学会了如何使用计算机软件等工具进行电力系统分析计算，这对我们以后的学习和工作都是有帮助的。 潮流计算部分： 本次实验潮流计算部分包括使用牛顿法对常规运行方式下的潮流进行计算，以及应用PQ分解法规划运行方式下的潮流计算。在规划潮流运行方式下，增加STNC-230母线负荷的有功至1.5.p.u，无功保持不变，计算潮流。潮流计算中，需要添加母线并输入所有母线 的数据，然后再添加发电机、负荷、交流线、变压器、支路，输入这些元件的数据。对运行方案和潮流计算作业进行定义，就可以定义的潮流计算作业进行潮流计算。 因为软件存在安装存在问题，无法使用图形支持模式，故只能使用文本支持模式，所以 无法使用PSASP绘制网络拓扑结构图，实验报告中的网络拓扑结构图均使用Visio绘制， 请见谅。 常规潮流计算： 下图是常规模式下的网络拓扑结构图，并在各节点标注电压大小以及相位。 下图为利用复数功率形式表示的各支路功率（参考方向选择数据表格中各支路的i侧母

线至j侧），因为无法使用图形支持模式，故只能通过文本支持环境计算出个交流线功率，下图为计算结果。

潮流计算问答题

1.什么是潮流计算？潮流计算的主要作用有哪些？ 潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。 对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置整定计算、电力系统故障计算和稳定计算等提供原始数据。 2.潮流计算有哪些待求量、已知量？ （已知量： 电力系统网络结构、参数； 决定系统运行状态的边界条件 待求量：系统稳态运行状态 例如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等）通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。 待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。 3.潮流计算节点分成哪几类？分类根据是什么？ （分成三类：PQ节点、PV节点和平衡节点，分类依据是给定变量的不同） PV节点（电压控制母线）：有功功率Pi和电压幅值Ui为给定。这种类型节点相当于发电机母线节点，或者相当于一个装有调相机或静止补偿器的变电所母线。 PQ节点：注入有功功率Pi和无功功率Qi是给定的。相当于实际电力系统中的一个负荷节点，或有功和无功功率给定的发电机母线。 平衡节点：用来平衡全电网的功率。平衡节点的电压幅值Ui和相角δi是给定的，通常以它的相角为参考点，即取其电压相角为零。 一个独立的电力网中只设一个平衡节点。 4.教材牛顿-拉夫逊法及有功-无功分解法是基于何种电路方程？可否采用其它类型方程？ 基于节点电压方程，还可以采用回路电流方程和割集电压方程等。但是后两者不常用。

电力系统分析潮流实验报告(优选.)

最新文件---------------- 仅供参考--------------------已改成-----------word文本 --------------------- 方便更改 赠人玫瑰，手留余香。 南昌大学实验报告 学生姓名：学号：专业班级： 实验类型：□验证□综合■设计□创新实验日期：实验成绩： 电力系统潮流计算实验 一、实验目的： 本实验通过对电力系统潮流计算的计算机程序的编制与调试，获得对复杂电力系统进行潮流计算的计算机程序，使系统潮流计算能够由计算机自行完成，即根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，由计算程序运行完成该电力系统的潮流计算。通过实验教学加深学生对复杂电力系统潮流计算计算方法的理解，学会运用电力系统的数学模型，掌握潮流计算的过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。 二、实验内容： 编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，完成该电力系统的潮流计算，要求计算出节点电压、功率等参数。 1、在各种潮流计算的算法中选择一种，按照计算方法编制程序。 2、将事先编制好的电力系统潮流计算的计算程序原代码由自备移动存储设备导入计算机。

3、在相应的编程环境下对程序进行组织调试。 4、应用计算例题验证程序的计算效果。 三、实验程序： function [e,f,p,q]=flow_out(g,b,kind,e,f) %计算潮流后efpq的终值 s=flow(g,b,kind,e,f); k=0; while max(abs(s))>10^-5 J=J_out(g,b,kind,e,f); J_ni=inv(J); dv=J_ni*s; l=length(dv)/2; for i=1:l e(i)=e(i)-dv(2*i-1); f(i)=f(i)-dv(2*i); end s=flow(g,b,kind,e,f); end l=length(e); for i=1:l s1=0; s2=0; for j=1:l s1=s1+g(i,j)*e(j)-b(i,j)*f(j); s2=s2+g(i,j)*f(j)+b(i,j)*e(j); end p(i)=e(i)*s1+f(i)*s2; q(i)=f(i)*s1-e(i)*s2;

潮流计算实验

电力系统分析实验报告 实验一：潮流计算的计算机算法 >> clear; n=10; nl=10; isb=1; pr=0.00001; B1=[120.03512+0.08306i0.13455i10; 230.0068+0.18375i0 1.023811; 140.05620+0.13289i0.05382i10; 450.00811+0.24549i0 1.023811; 160.05620+0.13289i0.05382i10; 460.04215+0.09967i0.04037i10; 670.0068+0.18375i0 1.023811; 680.02810+0.06645i0.10764i10; 8100.00811+0.24549i011; 890.03512+0.08306i0.13455i10] B2=[00 1.1 1.101; 001002; 00.343+0.21256i1002; 001002; 00.204+0.12638i1002; 001002; 00.306+0.18962i1002; 001002; 0.50 1.1 1.103; 00.343+0.21256i1002] Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); for i=1:nl if B1(i,6)==0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;

电力系统潮流计算实验报告

1． 手算过程 已知： 节点1：PQ 节点， s(1)= -0.5000-j0.3500 节点2：PV 节点， p(2)=0.4000 v(2)=1.0500 节点3：平衡节点，U(3)=1.0000∠0.0000 网络的连接图： 0.0500+j0.2000 1 0.0500+j0.2000 2 3 1）计算节点导纳矩阵 由2000.00500.012j Z += ? 71.418.112j y -=; 2000.00500.013j Z += ? 71.418.113j y -=; ∴导纳矩阵中的各元素： 42.936.271.418.171.418.1131211j j j y y Y -=-+-=+=; 71.418.11212j y Y +-=-=; 71.418.11313j y Y +-=-=; =21Y 71.418.11212j y Y +-=-=; 71.418.12122j y Y -==; 002323j y Y +=-=; =31Y 71.418.11313j y Y +-=-=; =32Y 002323j y Y +=-=; 71.418.13133j y Y -==; ∴形成导纳矩阵B Y ： ?? ?? ? ?????-++-+-+-+-+--=71.418.10071.418.10071.418.171.418.171.418.171.418.142.936.2j j j j j j j j j Y B 2）计算各PQ 、PV 节点功率的不平衡量，及PV 节点电压的不平衡量： 取:000.0000.1)0(1)0(1)0(1j jf e U +=+= 000.0000.1)0(2) 0(2)0(2j jf e U +=+= 节点3是平衡节点，保持000.0000.1333j jf e U +=+=为定值。 ()()[] ∑==++-=n j j j ij j ij i j ij j ij i i e B f G f f B e G e P 1 )0()0()0()0()0()0() 0(；

第三章简单电力系统的潮流计算汇总

第一章 简单电力系统的分析和计算 一、 基本要求 掌握电力线路中的电压降落和功率损耗的计算、变压器中的电压降落和功率损耗的计 算；掌握辐射形网络的潮流分布计算；掌握简单环形网络的潮流分布计算；了解电力网络的简化。 二、 重点内容 1、电力线路中的电压降落和功率损耗 图3-1中，设线路末端电压为2U 、末端功率为222~jQ P S +=，则 （1）计算电力线路中的功率损耗 ① 线路末端导纳支路的功率损耗： 222 2* 222~U B j U Y S Y -=?? ? ??=? ……………（3-1） 则阻抗支路末端的功率为： 222~~~Y S S S ?+=' ② 线路阻抗支路中的功率损耗： ()jX R U Q P Z I S Z +'+'==?2 2 22222 ~ ……（3-2） 则阻抗支路始端的功率为： Z S S S ~ ~~21?+'=' ③ 线路始端导纳支路的功率损耗： 2121* 122~U B j U Y S Y -=?? ? ??=? …………（3-3） 则线路始端的功率为： 111~ ~~Y S S S ?+'= ~~~图3-3 变压器的电压和功率 ~2 ? U （2）计算电力线路中的电压降落 选取2U 为参考向量，如图3-2。线路始端电压 U j U U U δ+?+=2 1 其中 2 2 2U X Q R P U '+'= ? ； 222U R Q X P U '-'=δ ……………（3-4） 则线路始端电压的大小： ()()2 221U U U U δ+?+= ………………（3-5） 一般可采用近似计算： 2 2 2221U X Q R P U U U U '+'+ =?+≈ ………………（3-6）

PQ分解法潮流计算实验

xxxx实验报告 学生姓名：学号：专业班级： 实验类型：□验证□综合■设计□创新实验日期：实验成绩： 一、实验项目名称 P-Q分解法潮流计算实验 二、实验目的与要求： 目的：电力系统分析的潮流计算是电力系统分析的一个重要的部分。通过对电力系统潮流分布的分析和计算，可进一步对系统运行的安全性，经济性进行分析、评估，提出改进措施。电力系统潮流的计算和分析是电力系统运行和规划工作的基础。 潮流计算是指对电力系统正常运行状况的分析和计算。通常需要已知系统参数和条件，给定一些初始条件，从而计算出系统运行的电压和功率等；潮流计算方法很多：高斯-塞德尔法、牛顿-拉夫逊法、P-Q分解法、直流潮流法，以及由高斯-塞德尔法、牛顿-拉夫逊法演变的各种潮流计算方法。 本实验采用P-Q分解法进行电力系统分析的潮流计算程序的编制与调试，获得电力系统中各节点电压，为进一步进行电力系统分析作准备。通过实验教学加深学生对电力系统潮流计算原理的理解和计算，初步学会运用计算机知识解决电力系统的问题，掌握潮流计算的过程及其特点。熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。 要求：编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，完成该电力系统的潮流计算，要求计算出节点电压、功率等参数。 三、主要仪器设备及耗材 每组计算机1台、相关计算软件1套 四、实验内容： 1.理论分析： P-Q分解法是从改进和简化牛顿法潮流程序的基础上提出来的，它的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓住主要矛盾，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功功率和无功功率迭代分开来进行。 牛顿法潮流程序的核心是求解修正方程式，当节点功率方程式采取极坐标系统时，

电力系统分析潮流计算大作业

电力系统分析潮流计算大作业（源程序及实验报告）

源程序如下： 采用直角坐标系的牛顿-拉夫逊迭代 function chaoliujisuan() m=3; %m=PQ节点个数 v=1;%v=PV节点个数 P=[-0.8055 -0.18 0]; %P=PQ节点的P值 Q=[-0.5320 -0.12 0]; %Q=PQ节点的Q值 PP=[0.5];%PP=PV节点的P值 V=[1.0];%V=PV节点的U值 E=[1 1 1 1.0 1.0]'; %E=PQ,PV,Vθ节点e的初值 F=[0 0 0 0 0]'; %F=PQ,PV,Vθ节点f的初值 G=[ 6.3110 -3.5587 -2.7523 0 0; -3.5587 8.5587 -5 0 0; -2.7523 -5 7.7523 0 0; 0 0 0 0 0; 0 0 0 0 0 ]; B=[ -20.4022 11.3879 9.1743 0 0; 11.3879 -31.00937 15 4.9889 0; 9.1743 15 -28.7757 0 4.9889; 0 4.9889 0 5.2493 0; 0 0 4.9889 0 -5.2493 ]; Y=G+j*B; X=[]; %X=△X n=m+v+1;%总的节点数 FX=ones(2*n-2,1);%F(x)矩阵 F1=zeros(n-1,n-1);%F(x)导数矩阵 a=0;%记录迭代次数 EF=zeros(n-1,n-1);%最后的节点电压矩阵 while max(FX)>=10^(-5) for i=1:m %PQ节点 FX(i)=P(i);%△P FX(n+i-1)=Q(i);%△Q for w=1:n FX(i)= FX(i)-E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)-F(i)*G(i,w)*F(w)-F(i)*B(i,w)*E(w); %△P FX(n+i-1)=FX(n+i-1)-F(i)*G(i,w)*E(w)+F(i)*B(i,w)*F(w)+E(i)*G(i,w)*F(w)+E(i)*B(i ,w)*E(w); %△Q end

简单环网潮流计算算例

简单环网潮流计算算例 图中所示为110kV 闭式电力网，A 为某发电厂的高压母线，其运行电压为117kV 。网络各元件的参数如下： 线路每公里的参数为 线路I 、II r 0 = 0.27 Ω， x 0 = 0.423 Ω， b 0 = 2.69×10-6S Ω III r 0 = 0.45 Ω， x 0 = 0.44 Ω， b 0 = 2.58×10-6S Ω 线路I 的长度为60km ，线路II 为50km ，线路III 为40km 。 各变电所每台变压器的额定容量、励磁功率和归算到110kV 电压级的阻抗分别为 变电所b S N =20MV A ，ΔS 0=0.05+j0.6MV A ，R T =4.84Ω，X T =63.5Ω 变电所c S N =10MV A ，ΔS 0=0.03+j0.35MV A ，R T =11.4Ω，X T =127Ω 负荷功率 S LDb =24+j18MV A ，S LDc =12+j9MV A 试求电力网的功率分布及最大电压损耗。 解： 1． 计算网络参数及制订等值电路 线路I ： var M 975.0var M 1101061.12 1Q S 1061.1S 60102.6938.25j 2.1660)423.0j 27.0(24BI 46I I -=???-=??=??=Ω +=Ω?+=---B Z 线路II ： var M 815.0var M 1101035.12 1Q S 1035.1S 50102.6915.21j 5.1350)423.0j 27.0(24BII 46II II -=???-=??=??=Ω +=Ω?+=---B Z 图（a ）

psasp潮流计算实验说明-2014

PSASP 潮流计算实验 一、实验目的 理解电力系统分析中潮流计算的相关概念，掌握用PSASP 软件对系统潮流进行计算的过程。学会在文本方式下和图形方式下的对潮流计算结果进行分析。 二、预习要求 复习《电力系统分析》中有关潮流计算的内容，了解有关潮流计算的功能，理解常用潮流计算方法，了解PQ 、PV 和V θ（平衡节点，在PSASP 中称为Slack 节点）的设置。 三、实验内容 （一） PSASP 潮流计算概述 潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。 PSASP 潮流计算的流程和结构如下图所示： 潮流计算 各种计算公共部分 图形方式 文本方式

以一个图所示9节点系统为例，计算其在常规、规划两种运行方式下的潮流。规划运行方式即在常规运行方式下，其中接于一条母线（STNB-230）处的负荷增加，对原有电网进行改造后的运行方式，具体方法为：在母线GEN3-230和STNB-230之间增加一回输电线，增加发电机3的出力及其出口变压器的容量，新增或改造的元件如下图虚线所示。 （二）数据准备 1. 指定数据目录及基准容量 双击PSASP图标，弹出PSASP封面后，按任意键，即进入PSASP主画面：

在该画面中，要完成的工作如下： (1)指定数据目录 第一次可通过“创建数据目录”按钮，建立新目录，如：F:\CLJS。以后可通过“选择数据目录”按钮，选择该目录。 (2)给定系统基准容量 系统基准容量项中，键入该系统基准容量，如100MVA。建立了数据之后，该数不要轻易改动。 2. 录入系统潮流计算数据 基础数据(系统参数)如下：

电力系统分析潮流计算的计算机算法

潮流计算的计算机算法实验报告 姓名：学号：班级： 一、实验目的 掌握潮流计算的计算机算法。 熟悉MATLAB，并掌握MATLAB程序的基本调试方法。 二、实验准备 根据课程内容，熟悉MATLAB软件的使用方法，自行学习MATLAB程序的基础语法，并根据所学知识编写潮流计算牛顿拉夫逊法(或PQ分解法) 的计算程序，用相应的算例在MATLAB上进行计算、调试和验证。 三、实验要求 每人一组，在实验课时内，用MATLAB调试和修改运行程序，用算例计算输出潮流结果。 四、实验程序 clear; %清空内存 n=input('请输入节点数:n='); n1=input('请输入支路数:n1='); isb=input('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入支路参数:B1='); B2=input('请输入节点参数:B2='); X=input('节点号和对地参数:X='); Y=zeros(n);

Times=1; %一：创建节点导纳矩阵 for i=1:n1 if B1(i,6)==0 %不含变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/B1(i,3); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4); Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4); else %含有变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3); Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3)); end end Y; %将OrgS、DetaS初始化 OrgS=zeros(2*n-2,1);

潮流计算(matlab)实例计算

潮流计算（matlab）实例计算 潮流例题：根据给定的参数或工程具体要求（如图） 本设计选择Matlab 进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿- 拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到 更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f（x） = 0 的单根 附 近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不 ，收集和查阅资料；学习相关软件（软件自选：

潮流计算（matlab）实例计算 平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y 。 （2）设个节点电压的初始值U 和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也 ２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化（增加，减少或保持不变）。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“ Matrix Laboratory ”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图

潮流计算 例题

例题：如图1所示的简单电力网中，已知变压器的参数为S N =31.5MV A ，0S S 031kW,190kW,%=10.5,%=0.7P P U I ?=?=；线路单位长度的参数为61110.21/km,0.416/km, 2.7410S/km r x b -=Ω=Ω=?。 如图所示的简单电力网中，当线路首端电压U A =120kV 时，试求： （1）线路和变压器的电压损耗； （2）变压器运行在额定变比时的低压侧电压及电压偏移。 说明：以上计算忽略电压降落的横分量。 图1 解：如题画等值电路图如下： 线路参数为： 0.21408.40.4164016.64l l l l R rl X x l ==?=Ω==?=Ω 变压器参数为 Ω=???=??=317.210) 105.31(110190103232 322N N S T S U P R Ω=?????=33.4010105.311001105.1010100%332 32N N S T S U U X （1） 变压器的功率损耗和励磁功率为 22 2T S 02N 2515()19031193.760.1937631.5S P P P kW MW S +?=?+?=?+== 222S 0N T N %%10.5(2515)0.731.5 3.0538var 10010010031.5100U S I S Q M S ?+??=+=+=?

1点处线路的充电功率 var 66308.01104074.22121222M lU b Q N l B =???== 计算L S 2 为：MVA j j Q Q Q j P P S B T LD T LD L 39.1719.25)66308.00538.315(19376.025)(22+=-+++=-?++?+= 线路阻抗中的功率损耗为： MW R U Q P P l L L l 65044.0104021.0110 173902519010322 23222222=???+=?+=?-- v a r 2885.11040416.0110 173902519010322 23222222M X U Q P Q l L L l =???+=?+=?-- 计算功率1S 为 M V A j j Q Q j P P S l L l I 68.1884.25)2885.139.17(65044.019.25)(221+=+++=?++?+= 线路电压损耗（忽略电压降落的横分量） 1125.848.418.6816.64 4.40120 =l l l A PR Q X U kV U +?+??== 1点电压为：1120 4.40115.60-A l U U U kV =?=-= 计算功率T S 2 为 M V A j j Q Q j P P S T LD T LD T 833.17163.25)833.215(163.025)(2+=+++='?++'?+= 变压器电压损耗 22125.163 2.3217.83340.33 6.73115.60 =T T T T T P R Q X U kV U +?+??== （2） 变压器低压侧折算到高压侧的电压为 2 1115.60 6.73108.87=-T U U U kV '?=-= 变压器低压侧的实际电压 22108.8710.8910=U U kV k '== 电压偏移为 2210.8910%100%8.9%10N N U U m U --= ?==

北京交通大学电力系统运行方式与潮流计算实验报告

电力系统分析 电力系统运行方式及潮流分析实验实验报告 姓名： 班级： 学号： 北京交通大学电气工程学院电气工程综合实验中心

实验1 电力系统运行方式及潮流分析实验 一、实验目的 1、掌握电力系统主接线电路的建立方法 2、掌握辐射形网络的潮流计算方法； 3、比较计算机潮流计算与手算潮流的差异； 4、掌握不同运行方式下潮流分布的特点。 二、实验内容 1、辐射形网络的潮流计算； 2、不同运行方式下潮流分布的比较分析 三、实验方法和步骤 1．辐射形网络主接线系统的建立 输入参数（系统图如下）： G1：300+j180MVA（平衡节点） 变压器B1：Sn=360MVA，变比=18/121，Uk％=14.3％，Pk=230kW，P0=150kW，I0/In=1％； 变压器B2、B3：Sn=15MVA，变比=110/11 kV，Uk％=10.5％，Pk=128kW，P0=40.5kW，I0/In=3.5％； 负荷F1：20+j15MVA；负荷F2：28+j10MVA； 线路L1、L2：长度：80km，电阻：0.21Ω/km，电抗：0.416Ω/km，电纳：2.74×10-6S/km（注意转换）。

辐射形网络主接线图 2．辐射形网络的潮流计算 （1）调节发电机输出电压，使母线A的电压为115kV，运行DDRTS进行系统潮流计算，在监控图页上观察计算结果，并填入下表： （2）手算潮流： 线路参数计算： =1/2*0.21*80=8.4Ω R L =1/2*0.416*80=16.64Ω X L =2*2.74*10-6 * 80=4.384*10-4s B L

电力系统运行方式及潮流分析实验报告.doc

电力系统第一次实验报告 ——电力系统运行方式及潮流分析实验

实验 1电力系统运行方式及潮流分析实验 一、实验目的 1、掌握电力系统主接线电路的建立方法 2、掌握辐射形网络的潮流计算方法； 3、比较计算机潮流计算与手算潮流的差异； 4、掌握不同运行方式下潮流分布的特点。 二、实验内容 1、辐射形网络的潮流计算； 2、不同运行方式下潮流分布的比较分析 三、实验方法和步骤 1．辐射形网络主接线系统的建立 输入参数（系统图如下）： G1：300+j180MVA （平衡节点） 变压器 B1：Sn=360MVA ，变比 =18/121，Uk％ =14.3％，Pk=230KW，P0=150KW，I0/In=1 ％； 变压器B2、B3： Sn=15MVA ，变比 =110/11 KV ，Uk ％=10.5％， Pk=128KW ，P0=40.5KW，I0/In=3.5 ％； 负荷 F1：20+j15MV A ；负荷 F2：28+j10MVA ； 线路 L1、L2 ：长度： 80km，电阻： 0.21 Ω/km，电抗： 0.416 Ω/km，电纳： 2.74 ×10-6 S/km。 辐射形网络主接线图

（ 1）在 DDRTS中绘出辐射形网络主接线图如下所示： （ 2）设置各项设备参数： G1：300+j180MVA （平衡节点）

变压器 B1：Sn=360MVA ，变比 =18/121，Uk％ =14.3％，Pk=230KW，P0=150KW，I0/In=1 ％； 变压器B2、B3： Sn=15MVA ，变比 =110/11 KV ，Uk ％=10.5％， Pk=128KW ，P0=40.5KW，I0/In=3.5 ％；

简单环形网络的潮流计算

银川能源学院 课程设计 课程名称：电力系统分析 设计题目：简单环形网络的潮流计算 学院：电力学院 专业：电气工程及其自动化 班级：电气（本）1202班 姓名：罗通 学号：1210240073 成绩： 指导教师：李莉、张彦迪 日期：2014年12月8日—2014年12月19日

潮流计算是在给定电力系统网络结构、参数和决定系统运行状态的边界条件的情况下确定系统稳态运行状态的一种基本方法，是电力系统规划和运营中不可缺少的一个重要组成部分。可以说，它是电力系统分析中最基本、最重要的计算，是系统安全、经济分析和实时控制与调度的基础。常规潮流计算的任务是根据给定的运行条件和网络结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。是电力系统研究人员长期研究的一个课题。它既是对电力系统规划设计和运行方式的合理性、可靠性及经济性进行定量分析的依据，又是电力系统静态和暂态稳定计算的基础。

前言------------------------------------------------------------------------------------------2 第一章：简单环形网络的潮流计算原理--------------------------------------4 1.1 电力线路和变压器上的功率损耗、电压降落及电能损耗--------------- 4 1.2电压降落、电压损耗、电压偏移及电压调整的概念---------------------- 5 1.3闭环网的潮流计算步骤---------------------------------------------------------- 6第二章：简单环形网络的潮流计算过程-------------------------------------- 7 2.1参数整理---------------------------------------------------------------------------- 7 2.2计算网络参数及等效电路------------------------------------------------------- 8 2.3电力系统潮流计算的运用------------------------------------------------------- 10 2.4注意事项---------------------------------------------------------------------------- 10 第三章：P-Q分解法的基本潮流算法-------------------------------------------11 3.1 P-Q分解法的原理----------------------------------------------------------------11 3.2 P-Q分解法的特点 ------------------------------------------ 13 3.3 P-Q分解法的潮流计算步骤 --------------------------------- 14 总结-------------------------------------------------------------------------------------------16谢辞-------------------------------------------------------------------------------------------17参考文献------------------------------------------------------------------------------------18