电力系统潮流计算课程设计终极版.docx

《电力系统潮流计算课程设计终极版.docx》由会员分享，可在线阅读，更多相关《电力系统潮流计算课程设计终极版.docx（13页珍藏版）》请在第一文库网上搜索。

1、目录摘要-11 .设计意义与要求-2-1.1 设计意义-2-12设计要求（具体题目）-2-2.题目解析一3-1 .1设计思路-3-2 .2详细设计-4-2.1.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）当用牛顿一拉夫逊法计算潮流时，给出修正方程和迭代收敛条件；系统如图1

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

7、可以相互转换。对于本题目，节点分析如下：节点1给出有功功率为2.,无功功率为1,PQ节点。节点2给出有功功率为0.5,电压幅值为10,PV节点。节点3电压相位是0,电压幅值为1,平衡节点。2.2.2待求量节点1待求量是P,Q；节点2待求量是Q,i；节点3待求量是U,io2.2.3导纳矩阵导纳矩阵分为节点导纳矩阵、结点导纳矩阵、支路导纳矩阵、二端口导纳矩阵。结点导纳矩阵：对于一个给定的电路（网络），由其关联矩阵A与支路导纳矩阵Y所确定的矩阵。支路导纳矩阵：表示一个电路中各支路导纳参数的矩阵。其行数和列数均为电路的支路总数。二端口导纳矩阵：对应于二端口网络方程，由二端口参数组成。节点导纳矩阵：以导

8、纳的形式描述电力网络节点注入电流和节点电压关系的矩阵。它给出了电力网络连接关系和元件特性的全部信息，是潮流计算的基础方程式。本例应用结点导纳矩阵具体计算时，根据如下公式：Yii=yio+Eyij丫汰由题给出的导纳可求的节点导纳矩阵如下：),=y12+13=1.25-j5.5进而节点导纳矩阵为：1.2j.5O.3-0.72./-O.31.守7jO.752.S.超42.2.4潮流方程网络方程是潮流计算的基础，如果给出电压源或电流源，便可解得电流、电压分布。然而，潮流计算中，这些值都是无法准确给定的。这样，就需要列出潮流方程。耳TQ*Vi对n个节点的网络，电力系统的潮流方程一般形式是/7(i=1,2

9、,n)=%户其中Pi=Pgi-P1di,Qi=Qc,i-Q,ai,即PQ分别为节点的有功功率无功功率。代入得潮流方程：上-=(1.25-j5.5)UG+(O.5j3)+(0.75-j2.5)IZOU1N1用2=(o.5-j3)Ug+(1.3-j7)1Z2+(0.8-j4)IZOoz24=(0.75-j2.5)-t1Z+(0.8-j4)1Z2+(1.55-j6.5)IZOo2.2.5牛顿一拉夫逊算法2.2.5.1牛顿算法数学原理:牛顿法（NewtonMethod）:解非线性方程f（x）-0的牛顿（NeWtOn）法,就是将非线性方程线性化的一种方法。它是解代数方程和超越方程的有效方法之设有单变量非

10、线性方程/=0,给出解的近似值.*）,它与真解的误差为心，则4）七更）将满足CH=o,即：将上式左边的函数在/）附近展成泰勒级数，如果差值Ax（）很小，-（）二次及以上阶次的各项均可略去得:这是对于变量的修正量以的线性方程式，成为修正方程，解此方程可得修正量：AS=/）用所求得的心（）去修正近似解，便得修正后的近似解同真解仍然有误差。为了进一步逼近真解，可以反复进行迭代计算，迭代计算通式是乂I)=乂Q迭代过程的收敛判据为I*卜卒幽y式中，Z和兀为预先给定的小正数。牛顿-拉夫逊法实质上就是切线法，是一种逐步线性化的方法，此法不仅用于求单变量方程，也适用于多变量非线性代数方程的有效方法。牛顿法至少

11、是二阶收敛的，即牛顿法在单根附近至少是二阶收敛的，在重根附近是线性收敛的。牛顿法收敛很快，而且可求复根，缺点是对重根收敛较慢，要求函数的一阶导数存在。2.2.5.2修正方程计算1、2节点的不平衡量为-4f节点3是平衡节点，其电压Yf十是给定的，故不参加迭代。根据给定的容许误差E=IO5,按收敛判据m咫令哈匕纤捻进行校验，以上节点1、2的不平衡量都未满足收敛条件，于是继续以下计算。修正方程式为：AVX(n=3)孙aP,aaraP,de泌0朗de,*皿,deiP2以de、也WP)HAAHA%de况e2明5V225V225V225V22de用e2也以上雅可比矩阵J中的各元素值是通过求偏导数获得的，对

12、PQ节点来说,A.、和Q,.是给定的，因而可以写出对PV节点来说，给定量是不因此可以列出当i时，雅可比矩阵中非对角元素为=TGijei+BJJ2P-MQRpcf0,8kj,1,j,U25(2=。组Sfj当/=i时,雅可比矩阵中对角元素为：1.25-5505Xg七与3-Q5-7。3-13-252301-55Q5。(GiJeJ-8)GiR-RJj/=I(Gijfj+BvJ-Gji+B；鸟j=1(Gjj+-GJ+Be03i=A1=(GjjAej-ijj)+GiR+BjcJjj=公=-2f.M,代入数值后的修正方程为:求解修正方程得:Ae1-0.25471-0.3611e20A九-0.1015计算节点

13、电压的第一次近似值2 .2.5.3收敛条件M)=雪)+.。)=it254空Q7453#=#)+邱)=O-Q361-Q361S)=雾)+/)=1-0=1)=。)+加。)=O-Q1O1与-QIO1式成立，结果收敛，迭代结束，计算平衡节点的功率和线路潮流计算，否则继续计算雅可比矩阵，解修正方程，直到满足收敛判据。3 .结果分析给定节点电压初值辛叁号经过四次笔算迭代过程后，得到节点电压和不平衡功率的变化情况分别于表1和表2所示（取W=IO5）:迭代计数k节点电压立=q=j1=G+7%10.7453-j0.36111-j0.1015120.4131-j0.35100.9901-j0.1479131.29

14、73-j0.37971.0083-j0.0185140.8217-j0.36440.9986-j0.08801表3.1迭代过程中节点电压变化情况迭代计数k节点不平衡量Q1P2v220-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语言计算矩阵，进行了四次迭代，结果仍然没有收敛。因此，我经过分析认为，次组初值，并不合适，无法收敛。我又阅读了有关文章，其中有如下结论：“多年的实践证明，牛顿法具有很好的二次收敛性，是求解多元非线性方程的经典算法，至今仍是电力系统潮流计算的主流。因此，一般认为算法不是导致不收敛的原因，潮流不收敛产生的主要原因是计算的初始条件给得不