电力系统分析课程设计-潮流计算和短路计算的程序实现

2012-2013学年度下学期电力系统分析课程设计电力系统的潮流计算和短路故障的计算机算法程序设计专业电气工程及其自动化姓名学号班级0310405指导教师钟建伟2013年4月14日信息工程学院课程设计任务书目录1电力系统图及初步分析 ........................ 错误!未指定书签。

1.1电力系统图及设计任务................... 错误!未指定书签。

1.2初步分析 ................................ 错误!未指定书签。

2牛顿-拉夫逊法简介............................ 错误!未指定书签。

2.1概述..................................... 错误!未指定书签。

2．2一般概念............................... 错误!未指定书签。

2.3潮流计算的修正方程..................... 错误!未指定书签。

2.4直角坐标表示的修正方程 ................ 错误!未指定书签。

3程序设计...................................... 错误!未指定书签。

3.1程序流程图.............................. 错误!未指定书签。

错误!未指定书签。

3.2潮流计算程序运行结果如下：............ 错误!未指定书签。

二．三相短路计算 .................................. 错误!未指定书签。

2.1计算原理：利用节点阻抗矩阵计算短路电流 . 错误!未指定书签。

注意：上述计算方法以及公式来源于电力系统分析上册P136-P137 ..................................... 错误!未指定书签。

2.2三相短路计算流程图：...................... 错误!未指定书签。

2.3习题实例.................................... 错误!未指定书签。

【例6-3】在如图2-3所示的电力系统中分别在节点1和节点5接入发电机支路，其标幺值参数为：............. 错误!未指定书签。

2.4三相短路计算程序及结果如下：............. 错误!未指定书签。

三．不对称短路计算................................ 错误!未指定书签。

3.1不对称短路课程设计的题目 ................. 错误!未指定书签。

3.2课程设计的设计任务及设计大纲............. 错误!未指定书签。

3.3电力系统不对称故障时元件的序参数和等值电路错误!未指定书签。

3.3.1电力系统不对称故障时用标幺值表示的各序等值电路错误!未指定书签。

3.4电力系统不对称故障时各序等值电路的化简与计算错误!未指定书签。

3.4.1正序等值电路的化简计算.............. 错误!未指定书签。

3.4.2负序等值电路的化简计算 ................. 错误!未指定书签。

3.4.3零序等值电路的化简计算 ................. 错误!未指定书签。

3.5电力系统不对称故障时元件参数的计算...... 错误!未指定书签。

3.5.1理论分析.............................. 错误!未指定书签。

3.5.2各元件各序等值电路电抗标幺值的计算错误!未指定书签。

3.6电力系统不对称故障分析与计算 ................ 错误!未指定书签。

3.6.1单相接地短路............................. 错误!未指定书签。

3.6.2两相直接接地短路 ........................ 错误!未指定书签。

3.6.3两相短路.............................. 错误!未指定书签。

注释：以上程序中的计算公式都是根据正序等效定则得到的。

错误!未指定书签。

3.7正序等效定则的内容 ........................ 错误!未指定书签。

3.8短路计算的matlab/simulink模型如下：.... 错误!未指定书签。

3.9.1变压器和线路参数设置：.............. 错误!未指定书签。

3.9.2短路模块和负载模块的参数设置....... 错误!未指定书签。

3.9.3故障相短路相电流和相电压波形....... 错误!未指定书签。

设计总结............................................ 错误!未指定书签。

参考文献............................................ 错误!未指定书签。

附录................................................ 错误!未指定书签。

一．潮流计算1电力系统图及初步分析1.1电力系统图及设计任务此电力系统图有AutoCAD2012软件画出网络各元件参数的标幺值如下：Z12=0.1+j0.4;y120=y210=j0.01538;z13=j0.13;k=1.1;z14= 0.12+j0.5;y140=y410=j0.01920;z24=0.08+j0.4;y240=y420=j0.01413系统中节点1，2为PQ节点，节点3为P节点，节点4为平衡节点，已给定P1s+jQ1s=-0.3-j0.18,P2s+jQ2s=-0.55-j0.13P3s=0.5,V3s=1.10,V4s=o0.1∠05容许误差为510-。

试用牛顿法计算潮流分布1.2初步分析潮流计算在数学上可归结为求解非线性方程组，其数学模型简写如下：2牛顿-拉夫逊法简介2.1概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。

这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。

牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵J ，朝减小方程的误差的方向前进一步，在新的点上再计算误差和雅可比矩阵，重复这一过程直到误差达到收敛标准，即得到了非线性方程组的解。

因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。

2．2一般概念对于非线性代数方程组即()0,,,21=n i x x x f Λ()n i Λ,2,1=(2－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (2－2) 上式称之为牛顿法的修正方程式。

由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆(2－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。

接着再从()1x 出发，重复上述计算过程。

因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆'(2－4)()()()k k k x x x ∆+=+1(2－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。

由式（2－4）和式子（2－5）可见，牛顿法的核心便是反复形成求解修正方程式。

牛顿法当初始估计值()0x 和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。

2.3潮流计算的修正方程运用牛顿－拉夫逊法计算潮流分布时，首先要找出描述电力系统的非线性方程。

这里仍从节点电压方程入手，设电力系统导纳矩阵已知，则系统中某节点（i 节点）电压方程为 从而得∑=**••=nj j ij i i U Y U S 1进而有()01=-+*=*•∑j nj ij i i i U Y U jQ P （2－6）式（2－6）中，左边第一项为给定的节点注入功率，第二项为由节点电压求得的节点注入功率。

他们二者之差就是节点功率的不平衡量。

现在有待解决的问题就是各节点功率的不平衡量都趋近于零时，各节点电压应具有的价值。

由此可见，如将式（2－6）作为牛顿－拉夫逊中的非线性函数()0=X F ，其中节点电压就相当于变量X 。

建立了这种对应关系，就可列出修正方程式，并迭代求解。

但由于节点电压可有两种表示方式——以直角做表或者极坐标表示，因而列出的迭代方程相应地也有两种，下面分别讨论。

2.4直角坐标表示的修正方程节点电压以直角坐标表示时，令i i i jf e U +=•、j j j jf e U +=•，且将导纳矩阵中元素表示为ij ij ij jB G Y +=，则式（2－7）改变为 ()()()()01=--+-+∑=n j j j ij ij i i i i jf e jB G jf e jQ P （2－7）再将实部和虚部分开，可得()()[]()()[]⎪⎪⎭⎪⎪⎬⎫=+---=++--∑∑==0011nj j ij j ij i j ij j ij i i nj j ij j ij i j ij j ij i i e B f G e f B e G f Q e B f G f f B e G e P （2－8）这就是直角坐标下的功率方程。

可见，一个节点列出了有功和无功两个方程。

对于PQ 节点（1,,21-=m i Λ，），给定量为节点注入功率，记为i P '、i Q '，则由式（2－8）可得功率的不平衡量，作为非线性方程 ()()[]()()[]⎪⎪⎭⎪⎪⎬⎫+---'=∆++--'=∆∑∑==nj j ij j ij i j ij j ij i i i nj j ij j ij i j ij j ij i i i e B f G e f B e G f Q Q e B f G f f B e G e P P 11 （2－9）式中i P ∆、i Q ∆——分别表示第i 节点的有功功率的不平衡量和无功功率的不平衡量。

对于PV 节点（n m m i ,,2,1Λ++=），给定量为节点注入有功功率及电压数值，记为i P '、i U '，因此，可以利用有功功率的不平衡量和电压的不平衡量表示出非线性方程，即有()()[]()⎪⎭⎪⎬⎫+-'=∆++--'=∆∑=22221i i i i nj j ij j ij i j ij j ij i i i f e U U e B f G f f B e G e P P（2－10）式中i U ∆为电压的不平衡量。