电力系统潮流计算课程设计..

《电力系统潮流计算课程设计..》由会员分享，可在线阅读，更多相关《电力系统潮流计算课程设计..（14页珍藏版）》请在装配图网上搜索。

1、目录摘要-1-1. 设计意义与要求-2-1.1设计意义-2-1.2设计要求（具体题目）-2-2. 题目解析-3-2.1设计思路-3-2.2详细设计-4-2.2.1节点类型-4-2.2.2待求量-4-2.2.3导纳矩阵-4-2.2.4潮流方程-5-2.2.5牛顿一拉夫逊算法-6-2.2.5.1牛顿算法数学原理：-6-2.2.5.2修正方程-7-2.2.5.3收敛条件-9-3. 结果分析-10-4. 小结-11-参考文献-12-摘要电力系统的出现，使高效，无污染，使用方便，易于调控的电能得到广泛应 用，推动了社会生产各个领域的变化，开创了电力时代，发生于第二次技术革命。 电力系统的规模和技术水准已

2、经成为一个国家经济发展水平的标志之一。电力系统稳态分析包括潮流计算和静态安全分析。电力系统潮流计算是电力 系统最基本的计算，也是最重要的计算。所谓潮流计算，就是已知电网的接线方 式与参数及运行条件，计算电力系统稳态运行各母线电压、个支路电流与功率及 网损。对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电 流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划 的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮 流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等 提供原始数据。在数学上，潮流计算是多元非线性方程组的求解问题，

3、求解的方法有很多种。 牛顿一拉夫逊法是数学上解非线性方程式的有效方法，有较好的收敛性。将牛顿 法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性、稀疏性 及节点编号顺序优化等技巧，使牛顿法在收敛性、占用内存、计算速度等方面都 达到了一定的要求。本文以一个具体例子分析潮流计算的具体方法，并运用牛顿一拉夫逊算法求 解线性方程。关键词：电力系统 潮流计算 牛顿一拉夫逊算法1. 设计意义与要求1.1设计意义潮流计算是电力系统分析中的一种最基本的计算，它的任务是对给定的运行 条件确定系统的运行状态，如各母线上的电压(幅值及相角)、网络中的功率分 布及功率损耗等。潮流计算的结果是电力系统稳定计

4、算和故障分析的基础。具体表现在以下方面：(1) 在电网规划阶段,通过潮流计算，合理规划电源容量及接入点，合理规划 网架，选择无功补偿方案，满足规划水平的大、小方式下潮流交换控制、调峰、调 相、调压的要求。(2) 在编制年运行方式时，在预计负荷增长及新设备投运基础上，选择典型方 式进行潮流计算，发现电网中薄弱环节,供调度员日常调度控制参考，并对规划、 基建部门提出改进网架结构，加快基建进度的建议。(3) 正常检修及特殊运行方式下的潮流计算，用于日运行方式的编制,指导发 电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要 求及电压质量要求。(4) 预想事故、设备退出运行对静态

5、安全的影响分析及作出预想的运行方式 调整方案。总结为在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比 较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电 力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力 系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系 统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采 用在线潮流计算。1.2设计要求(具体题目)如图所示的电网:1）根据给定的运行条件，确定图中电力系统潮流计算时各节点的类型、待求量;2）求节点导纳矩阵Y；3）给出潮流方程或功率方程的表达式；4）当用牛

6、顿一拉夫逊法计算潮流时，给出修正方程和迭代收敛条件;系统如图1.1所示2. 题目解析2.1设计思路此题目为负载电力系统潮流计算模型。首先写出节点导纳矩阵，并分析各节点的类型，找出待求量。然后，确定潮 流方程。最后进行潮流计算，而最后一步，可利用牛顿一拉夫逊法进行潮流分析。2.2详细设计2.2.1节点类型电力系统潮流计算中，节点一般分为如下几种类型：PQ节点：节点注入的有功功率无功功率是已知的；PV节点：节点注入的有功功率已知，节点电压幅值恒定，一般由无功储备 比较充足的电厂和电站充当；平衡节点：节点的电压为1*exp(0)，其注入的有功无功功率可以任意调节， 一般由具有调频发电厂充当。更复杂的

7、潮流计算，还有其他节点，或者是这三种节点的组合，在一定条件 下可以相互转换。对于本题目，节点分析如下：节点1给出有功功率为2.，无功功率为1，PQ节点。节点2给出有功功率为0.5，电压幅值为1.0，PV节点。节点3电压相位是0，电压幅值为1，平衡节点。2.2.2待求量节点1待求量是P，Q；节点2待求量是Q，5 ；节点3待求量是U，5。2.2.3导纳矩阵导纳矩阵分为节点导纳矩阵、结点导纳矩阵、支路导纳矩阵、二端口导纳矩 阵。结点导纳矩阵：对于一个给定的电路(网络)，由其关联矩阵A与支路导纳矩 阵Y所确定的矩阵。支路导纳矩阵：表示一个电路中各支路导纳参数的矩阵。其行数和列数均为电路的支路总数。二端

8、口导纳矩阵：对应于二端口网络方程，由二端口参数组成。节点导纳矩阵：以导纳的形式描述电力网络节点注入电流和节点电压关系的 矩阵。它给出了电力网络连接关系和元件特性的全部信息，是潮流计算的基础方 程式。本例应用结点导纳矩阵具体计算时，根据如下公式：Yii= yi0 + 】yijYik= -yik由题给出的导纳可求的节点导纳矩阵如下：y = y + y =1.25-j5.5y = y =0.5 + j12 21y = y = 0.75 + 2.513 31y = y + y = 1.3 j 722 2123y = y =0.8 + j 423 32y = y + y = 1.55 j 6.53331

9、32进而节点导纳矩阵为：1.25-j5.5Y = -0.5 + j3-0.75 + j2.5-0.5 + j3 -0.75 + j2.5 1.3-j7-0.8 + j4-0.8 + j4 1.55-j6.52.2.4潮流方程网络方程是潮流计算的基础，如果给出电压源或电流源，便可解得电流、电 压分布。然而，潮流计算中，这些值都是无法准确给定的。这样，就需要列出潮 流方程。对n个节点的网络，电力系统的潮流方程一般形式是P jQn二七匕 （i=1,2,，n）v .j=1 ji其中P= P - P , Q = Q - Q ，即PQ分别为节点的有功功率无功功率。i Gi Ldi i Gi Ldi代入得潮

10、流方程：2 - j1U /511二 (1.25-j5.5) - U1/5 1 + (0.5-j3) - 1/52 + (0.75-j2.5) 1/0 0.5 jq1/52二 (0.5-j3) - U*8 + (1.3-j7) - 1Z52 + (0.8-j4) - 1/0 乙-q =(0.75-j2.5) - U Z6 +(0.8-j4) - 1/5 +(1.55-j6.5) - 1/0 1/01 122.2.5牛顿一拉夫逊算法2.2.5.1牛顿算法数学原理:牛顿法(Newton Method)：解非线性方程 f(x)=0 的牛顿(Newton)法， 就是将非线性方程线性化的一种方法。它是解代

11、数方程和超越方程的有效方法之设有单变量非线性方程f G)=0，给出解的近似值x(。)，它与真解的误差为Ax(0)，则 x = x(0 + Ax(0 将满足 f (x)= 0，即：+ Ax(0将上式左边的函数在x(。)附近展成泰勒级数，如果差值Ax(。)很小，Ax(。)二次 及以上阶次的各项均可略去得：f C (0) + Ax (0 )= f C (0)+ f、( (0)Xx (0) = 0修正量：Ax(0 ) = - f用所求得的Ax(0)去修正近似解，便得 叫x Q = x (0 + Ax (0 = x (0 一 f f(0这是对于变量的修正量Ax(0)的线性方程式，成为修正方程，解此方程可

12、得 (0 )修正后的近似解xG)同真解仍然有误差。为了进一步逼近真解，可以反复进行迭代计算，迭代计算通式是 x +1)= x G)- fV *) f)迭代过程的收敛判据为或AxQ 式中，E和2为预先给定的小正数。牛顿-拉夫逊法实质上就是切线法，是一种逐步线性化的方法，此法不仅用于求单变量方程，也适用于多变量非线性代数方程的有效方法。牛顿法至少是二阶收敛的，即牛顿法在单根附近至少是二阶收敛的，在重根附近是线性收敛的。牛顿法收敛很快，而且可求复根，缺点是对重根收敛较慢，要求函数的一阶导数存在。2.2.5.2修正方程计算1、2节点的不平衡量A、AQ和AVii iAP（0）= P - P（0）= P

13、-11s 11s=-2-（G+ G + G）=-2-0 = -2111213e（0）E G e（0）- B f（0）+ f（0）E G f（0）+ B e（0）11j j1j j11j j1j jj=1j=1AQ（0） = Q -Q（0） = Q - f（0）G e（0）-B f（0）-e（0）G f（0）+ B e（0） 11j j1j j11j j1j jj=111s 11s=-1 -（B+ B + B ）=-1 - 0 = -11112132 = 0AV2（0 = V21 2Sj=1节点3是平衡节点，其电压Vi = ei + jf.是给定的，故不参加迭代。PQ),AQ(k),AV2(k)

14、 Le 进根据给定的容许误差e = 10-5，按收敛判据maxii i行校验，以上节点1、2的不平衡量都未满足收敛条件，于是继续以下计算。修正方程式为:AW = -JAV（n=3）AW = APAQ1Af11AV = lAe1AP2Ae2AV2t2Af 12SAP1SesaQ1SesaPSAP1SesaQ1SesaPSAP1七saQiSedAV/22-Sei2 卢 1 SAV 2 2- 副1Se-2 SAV 2 2-Se22户2SAV 22-Sf2以上雅可比矩阵J中的各元素值是通过求偏导数获得的，对PQ节点来说，P和Q是给定的，因而可以写出代 isa P = P - e (G e - B f)

15、- f (G f + Be)= 0、i is iij j ij j j . ij j ij jaq = Q - f (G e - B f)+e (G f + B e)= 0i is iij j ij j jij j对pv节点来说，给定量是P和V.,因此可以列出a P = P - e (Ge - B f)- f (G f + Be)= 0i is iij j ij j iAV 2 = V2 - ( 2 + f 2) = 0 iis i iSAP当j卫i时，雅可比矩阵中非对角元素为SAQi = -(G e + B f )SAeSAfij iij iSAP _ saQ _Sf = M = Be 一

16、GjfjjSAU 2 SAU 2 八=0SeSf当j=，时，雅可比矩阵中对角元素为:dAP8 (G e B f ) G e B fdeij j ij j ii i ii iij=1dAPVi =-乙(G f + B e ) G f + B edfij j ij j ii i ii ijj =1dAQr = 8 (G f + B e ) G f + B edeij j ij j ii i ii iij =1= 8 (G Ae- B f ) + G e + B fij j ij j ii i ii iij j ij j ii j=1dAU 2 仁i = 2edeijdAU 2AeiAAe2Af2

17、I0.25470.3611 0.1015i代入数值后的修正方程为:-1.255.50.55.51.2530.53 1.33 一0.57Ae 就 Ae=- 2-10.520020f_ 0 _求解修正方程得:计算节点电压的第一次近似值2.2.5.3收敛条件el = e(0)+ Ae(0)= 1 0.2547 = 0.7453 f(i)= f(0)+ f0)= 0 0.3611 = 0.3611 eG) = e(0)+ Ae(0)= 1 0 = 1f1 = fj) + f) = 0 0.1015 = 0.1015一轮迭代结束，根据收敛条件收敛判据max AP(k),AQk),AV2。) ii i式成

18、立，结果收敛，迭代结束，计算平衡节点的功率和线路潮流计算，否则继续 计算雅可比矩阵，解修正方程，直到满足收敛判据。3. 结果分析给定节点电压初值e（0）= e（0）= e（0）= 1.0, f（0）= f（0）= f（0）= 0 ,经过四次笔算123123迭代过程后，得到节点电压和不平衡功率的变化情况分别于表1和表2所示（取 = 10 -5）:迭代计数k节点电压口 1= + f 2 =e2 + jf2宁 3 =e+jf 3221110.7453-j0.3611221-j0.101533120.4131-j0.35100.9901-j0.1479131.2973-j0.37971.0083-j0

19、.0185140.8217-j0.36440.9986-j0.08801表3.1迭代过程中节点电压变化情况迭代计数k节点不平衡量AP1AQ1AP 2A 2 20-2-10.501-0.1482-0.9769-0.0726-0.01032-0.0902-0.6071-0.0480-0.00223-0.6272-4.3251-0.3610-0.01714-0.1816-1.2510-0.1042-0.0049表3.2迭代过程中节点不平衡量变化情况又利用C语言计算矩阵，进行了四次迭代，结果仍然没有收敛。因此，我经过分析认为，次组初值，并不合适，无法收敛。我又阅读了有关文章，其中有如下结论：“多年的实

20、践证明，牛顿法具有很好的二次收敛性，是求解多元非线性方程 的经典算法，至今仍是电力系统潮流计算的主流。因此，一般认为算法不是导致 不收敛的原因，潮流不收敛产生的主要原因是计算的初始条件给得不合理，导 致潮流方程无解。”中国自动化网.改善调度员潮流计算收敛性的措施4. 小结通过这次电力系统分析课程设计，我对自己以往的学习有了很深的认识。电 力系统分析专业课的系统学习刚刚结束，但是由于电力系统分析课程本身较难， 在学习过程中没有做到足够的理解学习以及在以后的时间没有做到充分的温习， 因此在这次设计作业的解题过程中感遇到了很大的阻力。首先拿到图就感到不知所云，只有第二问能轻松做出，其他几文完全犹如丈

21、 二和尚摸不着头脑。已经超出我们平时作业的水平，也超出了我们的能力范围， 以至好长时间都没有动笔，只是每次拿出来看几眼。慢慢看出了些眉头后，又通 过请教同学和老师把部分不懂的细节弄清了。完成这个艰难的第一步后，对整个 题目解析宏观的把握，对解题方法也有了大致思路。第二步就是详细分析，计算 线路参数作出等值电路，这一步相对进行的比较顺利。接下来就是理清思路认真 做潮流计算，在做潮流计算的时候由于计算量比较大，线路参数比较多，所花的 时间很多。更困难的是，其中很多计算如果用笔算，需要大量时间精力，而选择 编程计算机计算同样需要大量时间，经过斟酌，还是决定笔算。然而，意外总是 在最糟糕的情况下发生，

22、经过多次计算，我们仍无法得出答案，只好用计算机编 程计算。于是采用了我比较熟练的C语言编程，利用MicrosoftVisualCV60 进行计算。最终，计算机的结果与笔算结果一致，因此，我们也有了上述的的结 论。完成整个潮流计算后，面临这一个更艰巨的任务，就是手写稿转化成电子稿。 由于公式比较复杂，上下标多，对公式编辑及画图软件应用甚少，操作起来很费 时。通过此次的课程设计，让我认识到一个很严重的问题：我们的学习不扎实。 考完之后或是经过一段时间我们就对前面学过的知识没有映像了，这主要是因为 没有理解性学习。所以，在以后的学习中，一定要学的扎实，学的牢靠，并且能 学以致用，活学活用。同时，作为电气工程及其自动化专业本科生，我们对计算 机编程等使用还不够熟练，毕竟，以后的计算，基本由计算机完成，掌握好编程， 就如同海上航行把握住了风向，必然一帆风顺。另外，这次课程设计也为明年的毕业论文打下了很好的基础。参考文献1 何仰赞温增银电力系统分析华中科技大学出版社20062 何仰赞温增银电力系统分析题解华中科技大学出版社20063 中国自动化网4 于永源 杨绮雯 电力系统分析（第二版）中国电力出版社20075 谭浩强C程序设计第三版