1、电力系统分析课程报告题目：基于Matlab电力系统短路潮流计算系 别 电气工程系 专业班级 学生姓名 学生学号 指导教师 提交日期 2012年6月15日 目录一短路的原因及其后果1.1 短路的原因31.2 短路的后果3二短路类型2.1 对称短路-三相短路 32.2 不对称短路3三电力系统短路潮流计算的意义和目的3.1 短路潮流计算的意义43.2 短路潮流计算的目的5四短路潮流计算的基本原理五短路计算机计算的流程5.1 短路计算流程图65.2 短路计算源程序图75.3 运行计算结果及分析12六总结七参考文献八附录一、短路的原因及其后果 1.1短路的原因1、电气设备及载流导体因绝缘老化，或遭受机械
2、损伤，或因雷击、过电压引起绝缘损坏；2、架空线路因大风或导线履冰引起电杆倒塌等，或因鸟兽跨接裸露导体等；3、电气设备因设计、安装及维护不良所致的设备缺陷引发的短路；4、运行人员违反安全操作规程而误操作，如带负荷拉隔离开关，线路或设备检修后未拆除接地线就加上电压等。 1.2短路的后果1、强大的短路电流通过电气设备使发热急剧增加，短路持续时间较长时，足以使设备因过热而损坏甚至烧毁；2、巨大的短路电流将在电气设备的导体间产生很大的电动力，可能使导体变形、扭曲或损坏；3、短路将引起系统电压的突然大幅度下降，系统中主要负荷异步电动机将因转矩下降而减速或停转，造成产品报废甚至设备损坏；4、短路将引起系统中
3、功率分布的突然变化，可能导致并列运行的发电厂失去同步，破坏系统的稳定性，造成大面积停电。这是短路所导致的最严重的后果。5、巨大的短路电流将在周围空间产生很强的电磁场，尤其是不对称短路时，不平衡电流所产生的不平衡交变磁场，对周围的通信网络、信号系统、晶闸管触发系统及自动控制系统产生干扰。二、短路类型2.1、对称短路 三相短路 三相电流和电压仍是对称的 k(3) 2.2、不对称短路 ：两相短路： k(2)单相接地短路： k(1)两相短路接地： k(1, 1)相间短路与接地短路 ：相间短路：三相短路、两相短路接地短路：单相接地短路、两相短路接地三、电力系统短路潮流计算的意义和目的3.1短路潮流计算的
4、意义短路潮流计算是电力系统的一项重要分析功能，是进行故障计算，继电保护整定，安全分析的必要工具。MATLAB自1980年问世以来，它的强大的矩阵处理功能给电力系统的分析、计算带来许多方便。在处理短路潮流计算时，其计算机软件的速度已无法满足大电网模拟和实时控制的仿真要求，而高效的短路问题相关软件的研究已成为大规模电力系统仿真计算的关键。随着计算机技术的不断发展和成熟，对MATLAB短路潮流计算的研究为快速、详细地解决大电网的计算问题开辟了新思路。短路潮流计算是电力系统分析中的一种基本的计算，其任务时要根据给定的网络接线和其他已知条件，计算网络中的功率分布、功率损耗和未知的节点电压。采用计算机计算
5、，可以减轻计算人员的工作量，提高计算准确度和计算速度。短路潮流计算是电力系统分析中的一种基本的计算，其任务时要根据给定的网络接线和其他已知条件，计算网络中的功率分布、功率损耗和未知的节点电压。采用计算机计算，可以减轻计算人员的工作量，提高计算准确度和计算速度。3.2短路潮流计算的目的1、为选择和校验各种电气设备的机械稳定性和热稳定性提供依据，为此，计算短路冲击电流以校验设备的机械稳定性，计算短路电流的周期分量以校验设备的热稳定性；2、为设计和选择发电厂和变电所的电气主接线提供必要的数据；3、为合理配置电力系统中各种继电保护和自动装置并正确整定其参数提供可靠的依据。四、短路潮流计算的基本原理短路
6、电流计算的关键是化简网络以求得各电源点相对短路点的转移阻抗。那么根据上面介绍的数学模型，在程序设计时，关键也就在于如何化简节点导纳矩阵，求得各电源节点相对于短路点的互导纳，进而求得转移电抗，查曲线数字表求得短路电流。对总节点数为，电源节点数为的网络，设其节点导纳矩阵为，我们可以利用高斯消元法消去个中间节点和联络节点，将其化简为仅含有个电源节点和一个短路节点的等效网络。化简后的矩阵为转移导纳矩阵，转移导纳矩阵仍然是对称矩阵，其阶数为阶，矩阵中各元素的意义与化简前相同，即：对角线元素称为节点的自导纳，其值等于接于节点的所有支路导纳之和。非对角线元素称为节点、间的互导纳，它等于直接联接于节点、间的支
7、路导纳的负值。所以各电源节点与电源节点之间的转移阻抗就是转移导纳矩阵中对应的互导纳的倒数，于是便可求得各电源供给的短路电流。各电源供给的短路电流之和就是短路点的总短路电流。这一系列过程都可以由微机自动完成。其步骤如下。）准备工作：将网络中的各节点用从起的正整数编号，先电源节点后其它节点。）支路阻抗计算：各支路按其元件性质不同分别计算阻抗。）节点导纳矩阵形成：根据数学模型中介绍的节点导纳矩阵中各元素的意义，可由上一步中的各支路阻抗值形成节点导纳矩阵。）节点导纳矩阵化简：利用高斯消元法进行化简。）查曲线：运算曲线以曲线数字表的形式存放在磁盘中。应用时按电源的性质（水轮发电机、汽轮发电机或系统）区分
8、，调用不同的运算曲线。如果数值介于两数之间，则用插值法求结果。实践证明，采用线性插值得到的结果已经能够满足一般的计算要求。 五、短路计算机计算的流程5.1短路计算流程图输入数据用公式(6-10)计算短路电流İf输出结果用公式（6-9）计算指定支路电流用公式（6-11）计算各节点电压计算节点阻抗矩阵第f列元素选择故障点f形成节点导纳矩阵并进行三角分解 5.2短路计算源程序图Sbase_MVA=100.fid=fopen(Nodedata.txt);N=textscan(fid, %s %u %d %f %f %f %f %f %f);fclose(fid);busnumber=size(N1,1
9、)for i=1:busnumber Bus(i).name=N1(i); Bus(i).type=N2(i); Bus(i).no=i; Bus(i).Base_KV=N3(i); Bus(i).PG=N4(i); Bus(i).QG=N5(i); Bus(i).PL=N6(i); Bus(i).QL=N7(i); Bus(i).pb=N8(i); Bus(i).V=1.0; Bus(i).angle=0;endfid=fopen(Aclinedata.txt);A=textscan(fid, %s %s %f %f %f %f);fclose(fid);aclinenumber=size(
10、A1,1)for i=1:aclinenumber Acline(i).fbname=A1(i); Acline(i).tbname=A2(i); Acline(i).Base_KV=A3(i); Acline(i).R=A4(i); Acline(i).X=A5(i); Acline(i).hB=A6(i); for k=1:busnumber if strcmp(Acline(i).fbname, Bus(k).name) Acline(i).fbno=Bus(k).no; end if strcmp(Acline(i).tbname, Bus(k).name) Acline(i).tbn
11、o=Bus(k).no; end end endfid=fopen(Transdata.txt);T=textscan(fid, %s %f %f %s %f %f %f %f);fclose(fid);tansnumber=size(T1,1)for i=1:tansnumber Trans(i).fbname=T1(i); Trans(i).fbBase_KV=T2(i); Trans(i).fbrated_KV=T3(i); Trans(i).tbname=T4(i); Trans(i).tbBase_KV=T5(i); Trans(i).tbrated_KV=T6(i); Trans(
12、i).R=T7(i); Trans(i).X=T8(i); for k=1:busnumber if strcmp(Trans(i).fbname, Bus(k).name) Trans(i).fbno=Bus(k).no; end if strcmp(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).fbrat
13、ed_KV2/Trans(i).fbBase_KV2; Trans(i).X=tempx*Trans(i).X; Trans(i).R=tempx*Trans(i).R;end%N=0%Trans(1)%Trans(2)% for Y=G+jB matrixG,B,B2=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber); %B:=B;B2:=Bdlmwrite(Gmatrix.txt, G, delimiter, t,precision, 6);dlmwrite(Bmatrix.txt, B, delimiter,
14、t,precision, 6);%G%B%B2YN=G+j*B %不包括发电机导纳和负荷导纳的节点导纳矩阵pauseJP,JQ=FormJPQmatrix(Bus,B,B2,busnumber);JP;iJP=-inv(JP);JQ;iJQ=-inv(JQ);%pause%maxiteration=0for i=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-Bu
15、s(i).QL; dPGPL(i)=Bus(i).PG-Bus(i).PL; endNodeV=NodeV;Nodea=Nodea;%VX=VX%VY=VYdQGQL=dQGQL;dPGPL=dPGPL;%pause;%for nointer=1:10maxdP=1.;maxdQ=1.;epsilon=0.000001;noiteration=0;while (maxdPepsilon)&(maxdQepsilon)deltaP,deltaQ,maxdP,maxdQ=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B,G,busnumber);deltaP;
16、deltaQ;maxdPmaxdQda=iJP*deltaP;dV=iJQ*deltaQ;Nodea=Nodea+da;NodeV=NodeV+dV;noiteration=noiteration+1pause if noiteration20 break endendfor i=1:busnumber Bus(i).V=NodeV(i); NodeV(i)=NodeV(i)*Bus(i).Base_KV; Bus(i).angle=Nodea(i); Nodea(i)=Nodea(i)*180/pi;endnoiteration%潮流计算结束Nodea=NodeaNodeV=NodeV%潮流
17、计算节点电压模及相角y=1.0/(j*0.3)pausefor i=1:3 YN(i,i)=YN(i,i)+y; %将发电机导纳加入节点导纳矩阵endfor i=4:9 y=(Bus(i).PL-j* Bus(i).QL)/Bus(i).V2 YN(i,i)=YN(i,i)+y; %将负荷导纳加入节点导纳矩阵endYN %包括发电机导纳和负荷导纳的节点导纳矩阵uM=0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0; %以第四个元素为1的单位向量Z4=inv(YN)*uM %Z4:节点阻抗矩阵的第四列If=(Bus(4).V*cos(Bus(4).angle)+j*Bus(
18、4).V*sin(Bus(4).angle)/Z4(4) %母线对地短路电流标幺值IB230=100.0/sqrt(3)/230.0Ifymzh=abs(If*IB230) %母线对地短路电流有名值 1.6208kAfor i=1:busnumber V0(i)=Bus(i).V*cos(Bus(i).angle)+j*Bus(i).V*sin(Bus(i).angle); Vf(i)=V0(i)-Z4(i)*If;endV0Vf%由V0和Vf计算各支路电流程序 补上for i=1:busnumberfor k=1:busnumberI(i,k)=(Vf(k)-Vf(i)*YN(i,j); %
19、求出各支路电流 endend 5.3运行计算结果及分析Sbase_MVA =100busnumber =9aclinenumber =6tansnumber =3YN = 0 -17.3611i 0 0 0 +17.3611i 0 0 0 0 0 0 0 -16.0000i 0 0 0 0 0 +16.0000i 0 0 0 0 0 -17.0648i 0 0 0 0 0 0 +17.0648i 0 +17.3611i 0 0 3.3074 -39.3089i -1.3652 +11.6041i -1.9422 +10.5107i 0 0 0 0 0 0 -1.3652 +11.6041i 2
20、.5528 -17.3382i 0 -1.1876 + 5.9751i 0 0 0 0 0 -1.9422 +10.5107i 0 3.2242 -15.8409i 0 0 -1.2820 + 5.5882i 0 0 +16.0000i 0 0 -1.1876 + 5.9751i 0 2.8047 -35.4456i -1.6171 +13.6980i 0 0 0 0 0 0 0 -1.6171 +13.6980i 2.7722 -23.3032i -1.1551 + 9.7843i 0 0 0 +17.0648i 0 0 -1.2820 + 5.5882i 0 -1.1551 + 9.784
21、3i 2.4371 -32.1539imaxdP =1.9161e-007maxdQ =7.6235e-007noiteration =9Nodea = 0 9.6687 4.7711 -2.4066 -4.3499 -4.0173 3.7991 0.6215 1.9256NodeV = 18 18 18 227 220 224 229 227 231YN = 0 -20.6944i 0 0 0 +17.3611i 0 0 0 0 0 0 0 -19.3333i 0 0 0 0 0 +16.0000i 0 0 0 0 0 -20.3982i 0 0 0 0 0 0 +17.0648i 0 +17.3611i 0 0 3.3074 -39.3089i -1.3652 +11.6041i -1.9422 +10.5107i 0 0 0 0 0 0 -1.3652 +11.6041i 3.9159 -17.8835i 0 -1.1876 + 5.9751i 0 0 0 0 0 -1.9422 +10.5107i 0 4.1700 -16.1562i 0 0 -1.2820 + 5.5882i 0 0 +16.0000i 0 0 -1.1876 + 5.9751i 0 2.8047 -35.4456i -1.6171 +13.6980i
