电力系统潮流计算课程设计Word下载.doc
- 文档编号:13194664
- 上传时间:2022-10-08
- 格式:DOC
- 页数:20
- 大小:675.50KB
电力系统潮流计算课程设计Word下载.doc
《电力系统潮流计算课程设计Word下载.doc》由会员分享,可在线阅读,更多相关《电力系统潮流计算课程设计Word下载.doc(20页珍藏版)》请在冰豆网上搜索。
Us﹪=10.5,I0﹪=1.1;
变压器T4:
SFL1-16000/110,(110±
2.5﹪)/10.5,ΔPs=62kW,ΔP0=11.6kW,Us﹪=10.5,I0﹪=1.10。
导线型号均为LGJ-150,参数r0=0.21Ω/km,x0=0.4Ω/km,b0=2.8×
10-6S/km。
电网潮流计算
1.计算各元件参数,画出等值电路;
2.进行网络潮流计算;
3.不满足供电要求,进行调压计算。
2思路分析
这是一道潮流计算题,按照一般潮流计算的步骤将元件转换为等值参数,这里我们进行真实值的直接计算,并用近似计算计算。
由于负载给出,线路长度已知,我们可以将如图闭环的潮流计算分解成4个开环单电源的潮流问题进行计算,并计算是否有调压的必要。
3潮流计算过程
3.1各元件参数计算
①120Km线路
②100Km线路
③70Km线路
④变压器T1,T2
⑤变压器T3
⑥变压器T4
3.2绘制等效电路
图3-1等效电路
3.3功率分布计算
3.3.1各元件功率损耗
①两台T4变压器并联损耗
②T3变压器损耗
③100Km与70Km线路交点4末端功率损耗
④120Km与100Km线路交点3末端功率损耗
⑤1.4间100Km线路损耗
⑥1.3间120Km线路损耗
⑦2.4间70Km线路损耗
⑧2.3间100Km线路损耗
位置1点总损耗
位置2点总损耗
3.4调压计算
3.4.1计算1.4线路上的电压值
位置4由G1提供的电压为
由于119.84的输入电压大于110额定值,所以调压关系不满足。
(经验证其他3路均不满足关系)
所以需要降低121KV端的输出电压和提高110的输入额定值
T1的变压器取输出
T4的变压器取输入
调压后:
所以满足调压关系。
经验证其他三路均满足,调压成功。
结论:
将变压器T1,T2置于和
将变压器T3,T4置于调压档可满足条件。
4程序
clear;
n=8;
%input('
请输入节点数:
n='
);
nl=8;
请输入支路数:
nl='
isb=1;
请输入平衡母线节点号:
isb='
pr=0.00001;
请输入误差精度:
pr='
B1=[128.5+20.1i0.000556i10;
1413.6+32.16i0.0002224i10;
1613.6+32.16i0.0002224i10;
231.495+40.335i01.051;
451.78+53.885i01.0251;
4610.2+24.12i0.0001668i10;
671.495+40.335i01.0251;
686.8+16.08i0.0004448i10;
891.78+53.885i01.0251;
8108.5+20.1i0.000556i10];
请输入由支路参数形成的矩阵:
B1='
B2=[0022922901;
00220002;
050+30.987i220002;
040+27.79i220002;
060+37.18i220002;
200022922903];
请输入各节点参数形成的矩阵:
B2='
Y=zeros(n);
e=zeros(1,n);
f=zeros(1,n);
V=zeros(1,n);
sida=zeros(1,n);
S1=zeros(nl);
%-------修改部分------------
ym=1;
SB=100;
UB=220;
%ym=input('
您输入的参数是标么值?
(若不是则输入一个不为零的数值)'
ifym~=0
%SB=input('
请输入功率基准值:
SB='
%UB=input('
请输入电压基准值:
UB='
YB=SB./UB./UB;
BB1=B1;
BB2=B2;
fori=1:
nl
B1(i,3)=B1(i,3)*YB;
B1(i,4)=B1(i,4)./YB;
end
disp('
B1矩阵B1='
disp(B1)
n
B2(i,1)=B2(i,1)./SB;
B2(i,2)=B2(i,2)./SB;
B2(i,3)=B2(i,3)./UB;
B2(i,4)=B2(i,4)./UB;
B2(i,5)=B2(i,5)./SB;
B2矩阵B2='
disp(B2)
end
%%%---------------------------------------------------
fori=1:
nl %支路数
ifB1(i,6)==0 %左节点处于低压侧
p=B1(i,1);
q=B1(i,2);
else
p=B1(i,2);
q=B1(i,1);
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;
%对角元K侧
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
%对角元1侧
%求导纳矩阵
disp('
导纳矩阵Y='
disp(Y)
%----------------------------------------------------------
G=real(Y);
B=imag(Y);
%分解出导纳阵的实部和虚部
n %给定各节点初始电压的实部和虚部
e(i)=real(B2(i,3));
f(i)=imag(B2(i,3));
V(i)=B2(i,4);
%PV节点电压给定模值
n %给定各节点注入功率
S(i)=B2(i,1)-B2(i,2);
%i节点注入功率SG-SL
B(i,i)=B(i,i)+B2(i,5);
%i节点无功补偿量
%===================================================================
P=real(S);
Q=imag(S);
ICT1=0;
IT2=1;
N0=2*n;
N=N0+1;
a=0;
whileIT2~=0
IT2=0;
a=a+1;
fori=1:
ifi~=isb %非平衡节点
C(i)=0;
D(i)=0;
forj1=1:
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);
%节点功率P计算eiΣ(Gij*ej-Bij*fj)+fiΣ(Gij*fj+Bij*ej)
Q1=C(i)*f(i)-e(i)*D(i);
%节点功率Q计算fiΣ(Gij*ej-Bij*fj)-eiΣ(Gij*fj+Bij*ej)
%求P'
Q'
V2=e(i)^2+f(i)^2;
%电压模平方
%=========以下针对非PV节点来求取功率差及Jacobi矩阵元素=========
ifB2(i,6)~=3 %非PV节点
DP=P(i)-P1;
%节点有功功率差
DQ=Q(i)-Q1;
%节点无功功率差
%===============以上为除平衡节点外其它节点的功率计算=================
%=================求取Jacobi矩阵===================
forj1=1:
ifj1~=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/dfX3=dQ/de
X4=-X1;
%X1=dP/deX4=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;
elseifj1==i
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 潮流 计算 课程设计
![提示](https://static.bdocx.com/images/bang_tan.gif)