電力系統(tǒng)分析P-Q分解法潮流計算

華南理工大學(xué)廣州學(xué)院

專業(yè)課程設(shè)計報告

P-Q分解法潮流計算

系別電氣工程系專業(yè)班級09級電氣4班學(xué)生姓名鐘劍帆學(xué)號200930084778指導(dǎo)教師房大中提交日期2011年11月12日

華南理工大學(xué)廣州學(xué)院

目錄

P-Q分解法潮流計算....................................1

一、原理分析.........................................1

二、程序流程.........................................1

三、設(shè)計內(nèi)容.........................................2

3.1程序設(shè)計......................................2

3.2.1輸入數(shù)據(jù)....................................7

3.2.2輸出數(shù)據(jù)....................................8

3.2.3結(jié)果分析...................................11

四、心得總結(jié)........................................12

華南理工大學(xué)廣州學(xué)院

一、原理分析

從潮流計算的基本方程出發(fā)，采用PQ分解法并通過建立矩陣的修正方程來依次迭代，逐步逼近真值來計算電網(wǎng)的電壓和功率分布。

二、程序流程

1

華南理工大學(xué)廣州學(xué)院

三、設(shè)計Acline(i).fbno=Bus(k).no;

end

ifstrcmp(Acline(i).tbname,Bus(k).name)Acline(i).tbno=Bus(k).no;

end

end

2

華南理工大學(xué)廣州學(xué)院

end

fid=fopen(‘Transdata.txt’);

T=textscan(fid,‘%s%f%f%s%f%f%f%f’)

fclose(fid);

tansnumber=size(T{1},1)

fori=1:tansnumber

Trans(i).fbname=T{1}(i);

Trans(i).fbBase_KV=T{2}(i);

Trans(i).fbrated_KV=T{3}(i);

Trans(i).tbname=T{4}(i);

Trans(i).tbBase_KV=T{5}(i);

Trans(i).tbrated_KV=T{6}(i);

Trans(i).R=T{7}(i);

Trans(i).X=T{8}(i);

fork=1:busnumber

ifstrcmp(Trans(i).fbname,Bus(k).name)

Trans(i).fbno=Bus(k).no;

end

ifstrcmp(Trans(i).tbname,Bus(k).name)

Trans(i).tbno=Bus(k).no;

end

end

Trans(i).k=Trans(i).tbrated_KV*Trans(i).fbBase_KV/Trans(i).fbrated_KV/Trans(i).tbBase_KV;

tempx=Trans(i).fbrated_KV/Trans(i).fbBase_KV;

Trans(i).X=tempx*Trans(i).X;

Trans(i).R=tempx*Trans(i).R;

end

%N=0

%Trans(1)

%Trans(2)

%forY=G+jBmatrix

[G,B,B2]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber);%B:=B’;B2:=B"

dlmwrite(‘Gmatrix.txt’,G,‘delimiter’,‘\t’,’precision’,6);

dlmwrite(‘Bmatrix.txt’,B,‘delimiter’,‘\t’,’precision’,6);

G

B

B2

pause

[JP,JQ]=FormJPQmatrix(Bus,B,B2,busnumber);

JP

iJP=-inv(JP)

3

華南理工大學(xué)廣州學(xué)院

JQ

iJQ=-inv(JQ)

pause

%maxiteration=0

fori=1:busnumber

NodeV(i)=Bus(i).V;

Nodea(i)=Bus(i).angle;

VX(i)=Bus(i).V*cos(Bus(i).angle);

VY(i)=Bus(i).V*sin(Bus(i).angle);

dQGQL(i)=Bus(i).QG-Bus(i).QL;

dPGPL(i)=Bus(i).PG-Bus(i).PL;

end

NodeV=NodeV’

Nodea=Nodea’

%VX=VX’

%VY=VY’

dQGQL=dQGQL’

dPGPL=dPGPL’

pause

%fornointer=1:10

maxdP=1.;

maxdQ=1.;

epsilon=0.000001;

noiteration=0;

while(maxdP>epsilon)&(maxdP>epsilon)

[deltaP,deltaQ,maxdP,maxdQ]=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B,G,busnumber);

deltaP;

deltaQ;

maxdP;

maxdQ;

da=iJP*deltaP;

dV=iJQ*deltaQ;

Nodea=Nodea+da;

NodeV=NodeV+dV;

noiteration=noiteration+1;

ifnoiteration>20

break

end

end

fori=1:busnumber

Bus(i).V=NodeV(i);

NodeV(i)=NodeV(i)*Bus(i).Base_KV;

Bus(i).angle=Nodea(i);

4

華南理工大學(xué)廣州學(xué)院

Nodea(i)=Nodea(i)*180/pi;

end

noiteration

Nodea=Nodea’

NodeV=NodeV’

Clear

子函數(shù)

%生成G、B矩陣

function

[G,B,X]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber)

Y=zeros(busnumber);

X=zeros(busnumber);

fori=1:busnumber

Y(i,i)=Y(i,i)+Bus(i).pb*j;

end

fori=1:aclinenumber

f=Acline(i).fbno;

t=Acline(i).tbno;

Y(f,f)=Y(f,f)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);

Y(t,t)=Y(t,t)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);

Y(f,t)=Y(f,t)-1/(Acline(i).R+Acline(i).X*j);

Y(t,f)=Y(t,f)-1/(Acline(i).R+Acline(i).X*j);

X(f,f)=X(f,f)-1/Acline(i).X;

X(t,t)=X(t,t)-1/Acline(i).X;

X(f,t)=1/Acline(i).X;

X(t,f)=1/Acline(i).X;

end

fori=1:tansnumber

f=Trans(i).fbno;

t=Trans(i).tbno;

Y(f,f)=Y(f,f)+1/(Trans(i).R+Trans(i).X*j);

Y(t,t)=Y(t,t)+1/(Trans(i).R+Trans(i).X*j)/Trans(i).k;

Y(f,t)=Y(f,t)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k;

Y(t,f)=Y(t,f)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k;

X(f,f)=X(f,f)-1/Trans(i).X;

X(t,t)=X(t,t)-1/Trans(i).X;

X(f,t)=1/Trans(i).X;

X(t,f)=1/Trans(i).X;

end

G=real(Y);

B=imag(Y);

end

%生成JP、JQ矩陣

5

華南理工大學(xué)廣州學(xué)院

function[JP,JQ]=FormJPQmatrix(Bus,B,B2,busnumber)

JP=B;

JQ=B2;

fori=1:busnumber

ifBus(i).type==1

fork=1:busnumber

JQ(i,k)=0.;

JQ(k,i)=0.;

JP(i,k)=0.;

JP(k,i)=0.;

end

JQ(i,i)=1.;

JP(i,i)=1.;

end

ifBus(i).type==3

fork=1:busnumber

JQ(i,k)=0.;

JQ(k,i)=0.;

end

JQ(i,i)=1.;

end

end

end

%計算偏節(jié)點PQ差量

function

[deltaP,deltaQ,maxdP,maxdQ]=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B,G,busnumber)

deltaQ=dQGQL;

deltaP=dPGPL;

maxdP=0.;

maxdQ=0.;

fori=1:busnumber

ifBus(i).type==1

deltaQ(i)=0.;

deltaP(i)=0.;

end

ifBus(i).type==3

deltaQ(i)=0.;

%y1=0;

%y2=0;

y3=0;

fork=1:busnumber

if(B(i,k)~=0|G(i,k)~=0)

6

華南理工大學(xué)廣州學(xué)院

%y1=y1+(G(i,k)*VX(k)-B(i,k)*VY(k));

%y2=y2+(G(i,k)*VY(k)+B(i,k)*VX(k));

y3=y3+NodeV(k)*(G(i,k)*cos(Nodea(i)-Nodea(k))+B(i,k)*sin(Nodea(i)-Nodea(k)));

end

end

deltaP(i)=deltaP(i)-y3*NodeV(i);

%deltaP2(i)=(deltaP2(i)-(y1*VX(i)+y2*VY(i)))/Bus(i).V;

end

ifBus(i).type==2

%y1=0;

%y2=0;

y3=0;

y4=0;

fork=1:busnumber

if(B(i,k)~=0|G(i,k)~=0)

%y1=y1+(G(i,k)*VX(k)-B(i,k)*VY(k));

%y2=y2+(G(i,k)*VY(k)+B(i,k)*VX(k));

y3=y3+NodeV(k)*(G(i,k)*cos(Nodea(i)-Nodea(k))+B(i,k)*sin(Nodea(i)-Nodea(k)));

y4=y4+NodeV(k)*(G(i,k)*sin(Nodea(i)-Nodea(k))-B(i,k)*cos(Nodea(i)-Nodea(k)));

end

end

deltaP(i)=deltaP(i)-y3*NodeV(i);

%deltaP2(i)=(deltaP2(i)-(y1*VX(i)+y2*VY(i)))/Bus(i).V;

deltaQ(i)=deltaQ(i)-y4*NodeV(i);

%deltaQ2(i)=(deltaQ2(i)-(y1*VY(i)-y2*VX(i)))/Bus(i).V;

end

ifmaxdP<abs(deltaP(i));

maxdP=abs(deltaP(i));

end

ifmaxdQ<abs(deltaQ(i));

maxdQ=abs(deltaQ(i));

end

deltaP(i)=deltaP(i)/NodeV(i);

deltaQ(i)=deltaQ(i)/NodeV(i);

end

end

3.2.1輸入數(shù)據(jù)

節(jié)點數(shù)據(jù)（Nodedata.txt）

bus1.0.

70.

華南理工大學(xué)廣州學(xué)院

bus2318.1.6.

bus3318.0.8.

bus42.0.0.

bus52.250.50.

bus62.90.30.

bus72.0.0.

bus82.00.350.

bus92.0.0.

支路數(shù)據(jù)（Aclinedata.txt）

bus4bus5230.0.010.0850.088

bus4bus6230.0.0170.0920.079

bus5bus7230.0.0320.1610.153

bus6bus9230.0.0390.170.179

bus7bus8230.0.00850.0720.0745

bus8bus9230.0.01190.10080.1045

變壓器數(shù)據(jù)（Transdata.txt）

bus118.018.0bus40.0576bus218.018.0bus70.0625bus318.018.0bus90.0586

3.2.2輸出數(shù)據(jù)

Sbase_MVA=

100

N=

{9x1cell}[9x1uint32][9x1int32][9x1double]

[9x1double][9x1double][9x1double][8x1double]

busnumber=

9

A=

{6x1cell}{6x1cell}[6x1double][6x1double]

[6x1double]

aclinenumber=

6

T=

{3x1cell}[3x1double][3x1double]{3x1cell}

[3x1double][3x1double][3x1double]

tansnumber=

3

8double]double]double][9x1[6x1[3x1

華南理工大學(xué)廣州學(xué)院

9

華南理工大學(xué)廣州學(xué)院

NodeV=1

1

1

1

1

1

1

1

1

Nodea=

dQGQL=0000-0.5000-0.3000

10

華南理工大學(xué)廣州學(xué)院

-0.3500

dPGPL=

1.6300

0.8500

-1.2500

-0.9000

-1.0000

noiteration=

9

Nodea=

09.66874.7711-2.4066-4.3499-4.01733.79910.62151.9256

NodeV=

1818