潮流计算Word格式文档下载.docx
- 文档编号:21430138
- 上传时间:2023-01-30
- 格式:DOCX
- 页数:19
- 大小:166.42KB
潮流计算Word格式文档下载.docx
《潮流计算Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《潮流计算Word格式文档下载.docx(19页珍藏版)》请在冰豆网上搜索。
1.题目原始数据及其化简
原始数据:
节点数:
支路1:
支路2:
支路3:
支路4:
节点1:
节点2:
节点3:
节点4:
根据原始数据所画电路简化图如图1:
13
42
图1电路简化图
2.PQ分解法
2.1PQ分解法基本思想
PQ分解法是从改进和简化牛顿法潮流程序的基础上提出来的,它的基本思想是:
把节点功率表示为电压向量的极坐标形式,以有功功率误差作为修正电压向量角度的依据,以无功功率误差作为修正电压幅值的依据,这样,n-1+m阶的方程式便分解为一个n-1阶和一个m阶的方程,这两组方程分别进行轮流迭代,这就是所谓的有功-无功功率分解法。
牛顿法潮流程序的核心是求解修正方程式,当节点功率方程式采取极坐标系统时,修正方程式为:
或展开为:
(1)
以上方程式是从数学上推倒出来的,并没有考虑电力系统这个具体对象的特点。
在交流高压电网中,输电线路的电抗要比电阻大得多,系统中有功功率变化主要受电压相位的影响,无功功率则主要受母线电压幅值变化的影响。
在修正方程式的系数矩阵中,偏导数
和
的数值相对于偏导数
是相当小的,所以,矩阵N及J中各元素的数值相对是很小的,因此对牛顿法的第一步简化就是把有功功率和无功功率分开来进行迭代,即将式
(1)化简为:
(2)
这样,由于我们把2n阶的线性方程组变成了二个n阶的线性方程组,因而大大节省了机器内存和解题时间。
但是矩阵H和L都是节点电压幅值和相角差的函数,在迭代过程中仍然不断变化,而且又都是不对称矩阵。
对牛顿法的第二个化简,也是比较关键的一个化简,即把式
(2)中的系数矩阵简化为在迭代过程中不变的对称矩阵,即常数矩阵。
在一般情况下,线路两端电压的相角差是不大的(不超过10~20度),因此可以认为:
(3)
此外,与系统各节点无功功率相应的导纳
必定远远小于该节点自导纳的虚部,即:
或
(4)
考虑到以上关系后,式
(2)中系数矩阵中的元素表达式可以化简为:
(5)
这样,式(5)中系数矩阵可以表示为:
(6)
进一步可以把它们表示为以下矩阵的乘积:
(7)
将它代入
(2)中,并利用乘法结合率,可以把修正方程式变为:
(8)
及
(9)
将以上两式的左右两侧用以下矩阵左乘
就可以得到
(10)
(11)
以上两式就是P-Q分解法达到修正方程式,其中系数矩阵只不过是系统导纳矩阵的虚部,只是阶次不同,是对称矩阵,而且在迭代过程中维持不变。
它们与功率误差方程式
(12)
(13)
构成了P-Q分解法迭代过程中基本计算公式
2.2PQ分解法潮流计算基本步骤
1)形成系数矩阵
,并求其逆矩阵。
2)设PQ节点电压的初值和各节点相角初值为
(i=1,2,…,n,i≠s)和
(i=1,2,…,m,i≠s)。
3)通过(12)式计算各节点有功功率误差
从而求出
(i=1,2,…,n,i≠s)。
4)解修正方程式,求各节点电压相位角的变量
(i=1,2,…,n,i≠s)
5)求各节点电压相位角的新值
(i=1,2,…,n,i≠s)。
6)通过(13)式计算无功功率的不平衡量
,从而求出
7)解修正方程式,求各节点电压大小的变量
(i=1,2,…,m,i≠s)。
8)求各节点电压大小的新值
9)运用各节点电压的新值自第三步开始进入下一次迭代。
10)计算平衡节点功率和线路功率
3编程及运行
3.1PQ分解法潮流计算程序框图
是
否
图2PQ分解法潮流计算流程框图
3.2源程序代码
n=input('
请输入节点数:
n='
);
nl=input('
请输入支路数:
nl='
isb=input('
请输入平衡母线节点号:
isb='
pr=input('
请输入误差精度:
pr='
B1=input('
请输入由支路参数形成的矩阵:
B1='
%输入B1
B2=input('
B2='
%输入B2
X=input('
请输入由节点号及其对地阻抗形成的矩阵:
X='
%输入X
na=input('
请输入PQ节点数na='
Y=zeros(n);
YI=zeros(n);
e=zeros(1,n);
f=zeros(1,n);
V=zeros(1,n);
O=zeros(1,n);
fori=1:
n
ifX(i,2)~=0;
p=X(i,1);
Y(p,p)=1./X(i,2);
end
end
nl
ifB1(i,6)==0
p=B1(i,1);
q=B1(i,2);
elsep=B1(i,2);
q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
YI(p,q)=YI(p,q)-1./B1(i,3);
Y(q,p)=Y(p,q);
YI(q,p)=YI(p,q);
Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;
YI(q,q)=YI(q,q)+1./B1(i,3);
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
YI(p,p)=YI(p,p)+1./B1(i,3);
end%求导纳矩阵
G=real(Y);
B=imag(YI);
BI=imag(Y);
S(i)=B2(i,1)-B2(i,2);
BI(i,i)=BI(i,i)+B2(i,5);
end
P=real(S);
Q=imag(S);
e(i)=real(B2(i,3));
f(i)=imag(B2(i,3));
V(i)=B2(i,4);
n
ifB2(i,6)==2
V(i)=sqrt(e(i)^2+f(i)^2);
O(i)=atan(f(i)./e(i));
fori=2:
ifi==n
B(i,i)=1./B(i,i);
elseIC1=i+1;
forj1=IC1:
B(i,j1)=B(i,j1)./B(i,i);
fork=i+1:
forj1=i+1:
B(k,j1)=B(k,j1)-B(k,i)*B(i,j1);
p=0;
q=0;
ifB2(i,6)==2
p=p+1;
k=0;
forj1=1:
ifB2(j1,6)==2
k=k+1;
A(p,k)=BI(i,j1);
na
ifi==na
A(i,i)=1./A(i,i);
elsek=i+1;
forj1=k:
A(i,j1)=A(i,j1)./A(i,i);
na
A(k,j1)=A(k,j1)-A(k,i)*A(i,j1);
ICT2=1;
ICT1=0;
kp=1;
kq=1;
K=1;
DET=0;
ICT3=1;
whileICT2~=0|ICT3~=0
ICT2=0;
ICT3=0;
fori=1:
ifi~=isb
C(i)=0;
fork=1:
C(i)=C(i)+V(k)*(G(i,k)*cos(O(i)-O(k))+BI(i,k)*sin(O(i)-O(k)));
DP1(i)=P(i)-V(i)*C(i);
DP(i)=DP1(i)./V(i);
DET=abs(DP1(i));
ifDET>
=pr
ICT2=ICT2+1;
Np(K)=ICT2;
ifICT2~=0
fori=2:
DP(i)=B(i,i)*DP(i);
ifi~=n
IC1=i+1;
fork=IC1:
DP(k)=DP(k)-B(k,i)*DP(i);
else
forLZ=3:
i
L=i+3-LZ;
IC4=L-1;
forMZ=2:
IC4
I=IC4+2-MZ;
DP(I)=DP(I)-B(I,L)*DP(L);
O(i)=O(i)-DP(i);
kq=1;
L=0;
L=L+1;
C(i)=C(i)+V(k)*(G(i,k)*sin(O(i)-O(k))-BI(i,k)*cos(O(i)-O(k)));
DQ1(i)=Q(i)-V(i)*C(i);
DQ(L)=DQ1(i)./V(i);
DET=abs(DQ1(i));
ifDET>
ICT3=ICT3+1;
elsekp=0;
ifkq~=0;
L=0;
Nq(K)=ICT3;
ifICT3~=0
DQ(i)=A(i,i)*DQ(i);
forLZ=2:
i
L=i+2-LZ;
forMZ=1:
IC4
I=IC4+1-MZ;
DQ(I)=DQ(I)-A(I,L)*DQ(L);
DQ(k)=DQ(k)-A(k,i)*DQ(i);
L=L+1;
V(i)=V(i)-DQ(L);
kp=1;
K=K+1;
kq=0;
ifkp~=0
Dy(K-1,i)=V(i);
disp('
迭代次数'
disp(K);
每次没有达到精度要求的有功功率个数为'
disp(Np);
每次没有达到精度要求的无功功率个数为'
disp(Nq);
fork=1:
E(k)=V(k)*cos(O(k))+V(k)*sin(O(k))*j;
O(k)=O(k)*180./pi;
各节点的电压标幺值E为(节点号从小到大排):
'
disp(E);
各节点的电压V大小(节点号从小到大排)为:
disp(V);
各节点的电压相角O(节点号从小到大排)为:
disp(O);
forp=1:
C(p)=0;
forq=1:
C(p)=C(p)+conj(Y(p,q))*conj(E(q));
S(p)=E(p)*C(p);
各节点的功率S(节点号从小到大排)为:
disp(S);
各条支路的首端功率Sj(顺序同您输入B1时一样)为:
nl
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))));
disp(Si(p,q));
各条支路的末端功率Sj(顺序同您输入B1时一样)为:
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))));
disp(Sj(q,p));
各条支路的功率损耗DS(顺序同您输入B1时一样)为:
DS(i)=Si(p,q)+Sj(q,p);
disp(DS(i));
K
Cs(i)=i;
forj=1:
Dy(K,j)=Dy(K-1,j);
以下是每次迭代后各节点的电压值(如图所示)'
plot(Cs,Dy),xlabel('
),ylabel('
电压'
),title('
电压迭代次数曲线'
3.3运行程序及结果分析:
n=4
nl=4
isb=1
pr=0.00001
B1=[130.02+0.08i010;
140.04+0.12i010;
240.05+0.14i010;
340.04+0.12i010]
B2=[0-0.6-0.25i1.0002;
0-0.8-0.35i1.0002;
0.400.950.9503;
001.01.001]
X=[10;
20;
30;
40]
请输入PQ节点数na=2
迭代次数9
每次没有达到精度要求的有功功率个数为
333332220
每次没有达到精度要求的无功功率个数为
222211100
0.96411.0658+0.1659i0.9484+0.0546i0.9978+0.0666i
0.96411.07870.95001.0000
08.84563.29363.8196
-1.1374+0.2500i0.8000+0.3500i0.4000-0.6290i0.0000+0.2198i
-0.5746+0.3328i
-0.5628-0.0828i
0.8000+0.3500i
-0.1841-0.3341i
0.5841-0.2948i
0.5767+0.1246i
-0.7672-0.2583i
0.1905+0.3535i
0.0095+0.0379i
0.0139+0.0418i
0.0328+0.0917i
0.0065+0.0194i
每次迭代后各节点的电压值如图所示
图3电压迭代次数曲线
经过九轮迭代,节点功率不平衡量下降到0.000010以下,迭代到此结束,电压幅值和相角都能够满足计算精度的要求。
4.小结
通过本次电力系统分析实验,我更加熟悉地掌握了matlab的用法,在之前也用过matlab进行过实验的分析,都没能十分熟练地运用该软件,这次实践过后,加深了对此软件的理解以及更多的应用。
PQ潮流计算法在运用matlab时会显得比较简单,相比别的软件。
就自动化专业来讲,matlab会是一个很方便去解决问题的软件,因此在之后的学习过程中应该更加重视它,更近一步的掌握对以后的专业学习会有很大帮助。
电力系统分析是今年学习的一门比较重要的专业课,电力系统潮流计算又是这门课非常重要的内容,因此此次课设意义重大。
PQ法是进行潮流计算的一种方法,通过给定条件确定各母线的电压幅值和相角,各元件流过的功率和整个系统的功率损耗。
在做课设的过程中也遇到过很多问题,比如对课本知识掌握不透没能完全理解实验内容耽误了一些时间和精力。
但在同学的帮助以及查阅资料的情况下慢慢的克服了个个问题,完成了实验报告。
感悟到无论做什么事都应该平心静气,认认真真。
5.参考文献
[1]何仰赞等.电力系统分析[M].武汉:
华中理工大学出版社,2002.3
[2]西安交通大学等.电力系统计算[M].北京:
水利电力出版社,1993.12
[3]陈衍.电力系统稳态分析[M].北京:
水利电力出版社,2004.1
[4]李光琦.电力系统暂态分析[M].北京:
水利电力出版社,2002.5
[5]于永源,杨绮雯.电力系统分析(第二版)[M].北京:
中国电力出版社,2004.3
[6]华智明.张瑞林.电力系统分析.北京.重庆大学出版社.2002.5
[7]弋东方.电力工程电气设计手册.北京.中国电力出版社.1989.10
本科生课程设计成绩评定表
指导教师签字:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 潮流 计算