《电力系统动态潮流计算及网络拓扑分析.pdf》由会员分享，可在线阅读，更多相关《电力系统动态潮流计算及网络拓扑分析.pdf（81页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、山东大学硕士学位论文电力系统动态潮流计算及网络拓扑分析姓名：张国衡申请学位级别：硕士专业：电路与系统指导教师：王良20070515原创性声明本人郑重声明：所呈交的学位论文，是本人在导师的指导下，独立进行研究所取得的成果。除文中已经注明引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写过的科研成果。对本文的研究作出重要贡献的个人和集体，均已在文中以明确方式标明。本声明的法律责任由本人承担。论文作者签名：；础日关于学位论文使用授权的声明本人完全了解山东大学有关保留、使用学位论文的规定，同意学校保留或向国家有关部门或机构送交论文的复印件和电子版，l 允许论文被查阅和借阅；本人授权山东大学可以
2、将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或其他复制手段保存论文和汇编本学位论文。(保密论文在解密后应遵守此规定)论文作者签名：知红四母师签名：二五t 日期：1 2 丘f 山东大学硕士学位论文电力系统动态潮流计算及网络拓扑分析摘要电力系统潮流计算是电力系统规划设计与运行分析的基本工具。通过几十年的发展，潮流算法日趋成熟。但由于电网的复杂性，传统的潮流算法依然存在着一些方面的局限性。在电网各种运行方式中，节点注入功率的改变特别是节点注入停运将使系统节点有功、无功注入发生较大的变化，使系统功率出现严重不平衡，而使以往的潮流计算方法在计算这种情况下往往会出现收敛性差、计算
3、结果与实际不符的情况。本文提出的动态潮流算法，主要是在常规潮流计算的基础上考虑了负荷和发电机的静态频率特性，其核心是潮流计算和频率计算。在动态潮流计算中，系统中由于功率扰动(切负荷、发电机增减出力)而产生的不平衡功率按照各发电机和负荷的功频静特性系数在多台发电机及负荷之间进行分配，得到调整后的发电机出力和负荷的大小以及系统的频率连续变化的情况，这克服了常规潮流算法中由平衡节点独自承担不平衡功率而导致潮流收敛性差、结果和实际不符的情况。显然，相对于常规潮流算法来说，动态潮流算法能够对系统运行的实际情况进行更有效的模拟，是一个较大的进步。电力系统在正常运行情况下，可用电力系统状态参数之间的代数方程
4、组来描述其某个特定运行状态。包括潮流计算的许多程序都是以导纳矩阵为基础的。而结点导纳矩阵又是随网络拓扑而变，若不能及时而准确地随着开关所处状态实时变化而修改网络拓扑结构，就会造成分析计算结果的错误。本文将电力系统网络拓扑分析和电力系统的动态潮流两个问题衔结起来，构成一个整体。可实时地随着开关等电气元件信息的变化，进行动态潮流计算，这在工程实际中也有非常广泛的应用。论文还通过典型算例系统的仿真，分析了动态潮流算法的合理性，说明了新算法为准确分析电力系统随负荷变化后的潮流状态提供了一种比较理想的工具。并应用动态潮流和网络拓扑分析成功的模拟了电网的解列和并网操作。关键词：动态潮流；P-Q 分解法；稀
5、疏矩阵；节点导纳矩阵；因子表；频率静特性；网络拓扑分析山东大学硕士学位论文I n v e s t i g a t i o no nD y n a m i cL o a dF l o wa n dN e t w o r kT o p o l o g i c a lA n a l y s i sA b s t r a c tP o w e rf l o wc a l c u l a t i o ni sab a s i ct o o lf o rs v s t e mp r o g r a m m i n ga n do p e r a t i o n a la n a l y z i n gi
6、np o w e rs y s t e m T h ea l g o r i t h mh a sb e e np u tf o r w a r db yt h ed e e pr e s e a r c ho fr e s e a r c h e r si na n da b r o a df o rs e v e r a ld e c a d ey e a r s B u tb e c a a s eo ft h ec o m p l e x i t yo fp o w f f fs y s t e r n，m o r eo rl e s st h et r a d i t i o n a
7、 la l g o r i t h m sh a v es o m ec e r t a i nl i m i t a t i o n s T h et r a d i t i o n a l1 0 a df l o wi s n te a s yt oc o n v e r g ew h e nal a r g eu n b a l a n c ep o w e re x i s t sw h i c hc o u r s e db yt h ec h a n g eo fn o d ei n j e c t i o np o w e r A n dt h er e s u l td o e
8、 s n tc o n f o r mt ot h er e a l i t y T h ed y n a m i cl o a df l o wm o d e lb ep u tf o r w a r di nt h i sp a p e ri st oh a v eb e e nc o n s i d e r e dt h ef r e q u e n c ys t a t i ec h a r a e t e r i s t i co fl o a da n dg e n e r a t o ro nt h ef o u n d a t i o no ft r a d i t i o n
9、 a ll o a df l o w,t l l ec o r ei sf l o we a l c u l a t i o na n df r e q u e n c yc a l c u l a t i o n I nd y n a m i cl o a df l o wc a l c u l a t i o nt h eu n b a l R a c ep o w e rw h i c hc o u r s e db yt h ec h a n g eo f n o d ei n j e c t i o np o w e ri sd i s t r i b u t e da c c o
10、r d i n gt ot h ef r e q u e n c ys t a t i cc h a r a c t e r i s t i eo fe a c hg e n e r a t o rb e t w e e na l lg e n e r a t o r s，a n dg e t st h ec h a n g eo f s y s t e m a t i cf r e q u e n c yf o l l o w1 0 a da sw e l la sg e n e r a t o r s。a c t i v ep o w e r S o，i th a ss u r m o u
11、 n t e dt h ec o n d i t i o nt h a tf l o wc o n v e r g e n c et ob en o tg o o da n dr e s u l ti si n e o n s i s t e n tw i t hr e a l i t yw h i c hs i n c eb a l a n t e dn o d eu n d e r t a k e su n b a l a n c e dp o w e ra l o n ei nt r a d i t i o n a lI o a df l o w O b v i o u s l y,r
12、e l a t i v et ot r a d i t i o n a ll o a df l o w。d y n a m i c1 0 a df l o wa l g o r i t h mc a ng oo nf o rt h ea c t u a lc o n d i t i o nt h a ts y s t e mI n n s，i sm o r ee f f e c t i v ea n ds i m u l a t e d,i ti sag r e a ta d v a n c e Al o to fp o w e rp r o g r a m st h a ti n c l u
13、 d ef l o we a l c u l a t i o na r eb a s i ew i t hb u sa d m i t t a n c em a t r i x A n db u sa d m i t t a n c em a t r i xi sa l s ot on a a k em b b m gw i t ht h ec h a n g eo fn e t w o r kt o p o l o 西c a ls t r u c t u r e I fw eC a l ln o tm o d i f yn e t w o r kt o p o l o g i c a ls
14、t r u c t u r ei nt i m ea c c u r a t e l ya l o n gw i t l lt h es w i t c hl o c a t e ds t a t ec h a n g eo f r e a lt i m et om a k er u b b i n g，w i l lc a u s er e s u l tm i s t a k e T h et h e s i sl i n k su Dd y n a m i cl o a df l o wa n a l y s i sa n dn e t w o r kt o p o l o g i c
15、a la n a l y s i si np o w e rs y s t e ma saw h o l et ot a k ed y n a m i cl o a df l O Wc a l c u l a t i o no fr e a lt i m e T h i sp r o j e c th a sb r o a da p p l i c a t i o ni np r a c t i c e T h et h e s i sa l s ov a l i d a t e st h er a t i o n a l i t yo ft h ed y n a m i cp o w e r
16、f l o wa l g o r i t h mt h r o u g he m u l a t i n go nat y p i c a ls y s t e m T h er e s u l ti n d i e a t 酷t h a tt h en e wa l g o r i t h mh a sp r o v i d e dam o r ei d e a lt o o lf o ra n a l y z i n gt h ep o w e rf l o ws t a t eo fD 0 w e rs y s t e ma f t e rt h ec h a n g eo fs y s
17、 t e ml o a d sw e l la n dt r u l y。F i n a l l yw eh a v es i m u l a t e du n t i er o wa n dp a r a l l e lo p e r a t i o no fd e c t r i c a ln e t w o r ks u c c e s s f u l l ya c c o r d i n gt ot h et h e o r yo f d y n a m i cl o a df l o wa n dn e t w o r kt o p o l o g i c a la n a l y
18、s i s K e yw o r d s：D y n a m i cl o a df l o w,p-Qd e c o m p o s i t i o nm e t h o d，S p a r s em a r x，B u sa d m i t t a n c em a t r i x，F r e q u e n c ys t a t i cc h a r a c t e r i s t i c；N e t w o r kt o p o l o g i c a la n a l y s i s2山东大学硕士学位论文第1 章绪论1 1课题背景当今，电能以其清洁、高效、便于输送等突出优点，己经成为
19、全球广泛使用的最主要能源，在社会经济的发展中起到举足轻重的作用。电力系统是当今世界最庞大的人工系统，它包括发电、输电、配电、用电四个环节，对其运行的最基本要求有三点：保证安全可靠的供电、要有合乎要求的电能质量、要有良好的经济性。随着国民经济的进一步发展，社会各部门对电力系统的发展提出了更新、更高的要求，当今电力系统正向着超高压、大容量、远距离的输电方式发展，全国电网区域互联也成为一种必然的趋势，这样的发展趋势在很大程度上提高了系统运行效率，增加了经济效益，促进了对能源的合理开发和利用，减轻了对环境保护的压力，但同时也给电力系统的安全运行带来了新的问题。经过近几年的技术改造和升级，大多数电网数据
20、采集与监控(S C A D A)系统都已通过实用化验收。在此基础上，对电力系统调度自动化的要求也不断提高，能量管理系统(E M S)在调度中心的应用逐步走向实用化。一个实用化的能量管理系统包括：数据库管理系统、人机管理系统、经济调度系统、网络管理系统和高级应用软件系统，五大部分有机地结合在一起，为用户提供服务【l】。其中，高级应用软件系统是整个能量管理系统的灵魂，只有它的正常工作才可以为用户提供分析与控制服务。高级应用软件系统又由实时网络状态估计、在线潮流、安全分析、最优潮流等可以不断延伸的各种分析与控制功能组成。在这些众多功能中，潮流计算是基本的工具，它根据给定的运行条件及系统接线情况确定整
21、个电力系统各部分的运行状态，是电力系统运行、规划以及安全性、可靠性和优化的基础，也是各种暂态分析的基础和出发点。现代电力系统规模庞大，结构复杂，其规划、设计和运行均需借助电子计算机进行潮流计算分析。近几十年来，随着电力系统规模的不断扩大与结构的日益复杂、电力系统自动化水平的提高和计算机技术的日新月异，电力系统潮流的计算机算法也在不断地进步与更新。但是，从计算效率、收敛性以及对实际系统的有效模拟程度等各个方面综合起来看，现有的诸多潮流算法仍然存在不少尚待解决的问题。所以，还需对电力系统潮流进行更深入的研究，发展更加完善的潮流算法，以满足处于不断扩大与更新中的当代电力系统的需求。1。2潮流计算的基
22、本要求和要点潮流计算随计算性质不同而有不同的要求，如长距离输电、区域性网络、城市配电网络等都有不尽相同的要求，但仍有其共同的基本要求。首先是不同类型的网络在各种运行方式下，网络各节点的电压水平应符合有关规定。其次如网络中各线路的潮流分布山东大学硕士学位论文不应有线路过载等。对潮流计算的分析主要根据计算的目的而定。在电力系统运行方式中一般含高峰负荷和低谷负荷时运行方式下，在具有水力发电厂的电力系统中根据水电厂水文特点又有丰水期、平水期、枯水期的运行方式，此外，也需要研究事故运行方式和各种特殊运行方式。在潮流计算中首先应效验网络枢纽点的电压水平及网络各节点的电压是否满足要求，其次效验各发电厂发电机
23、的有功及无功出力是否符合技术要求，另外根据计算的要求对各线路、变压器的潮流进行分析。1 3潮流计算程序的发展电力系统潮流计算是研究电力系统稳态运行情况的一项基本运算，其数学本质是一组多元非线性方程，主要采用迭代的方法求解。电力系统潮流计算从提出至今，经历了一个由手工，利用交、直流计算台到应用数字电子计算机的发展过程，现有的潮流算法都以计算机的应用为前提。由于潮流计算在电力系统分析研究中具有重要的地位，吸引了大量的专家学者对其进行了研究，针对各种实际情况以及特殊需求，发展了多种用于电力系统潮流计算的计算机算法。在现有的潮流算法之中，最早出现的是常规潮流算法，其它潮流算法都是根据不同的实际需求在常
24、规潮流的基础上发展起来的。利用电子计算机进行电力系统潮流计算始于上个世纪5 0 年代中期，最初以节点导纳矩阵为基础进行迭代，其原理简单，易于编程实现，同时由于导纳矩阵是稀疏矩阵，对计算机内存需求不大。但是此法收敛性较差，其迭代次数会随着系统规模的扩大而急剧增加，易出现不收敛的情况。在这种情况下，出现了基于阻抗矩阵的迭代方法，大大改善了潮流计算的收敛性，可以求解一些用导纳法无法收敛的潮流问题。但是，阻抗矩阵是满秩矩阵，不但占用内存大，而且每次迭代所需的计算量也比较大，这就引入了新的问题。随后出现的分块阻抗法，它将一个大系统分为若干小系统，只需存储各个小系统的阻抗矩阵以及它们之间联络线的阻抗，此法
25、能够在一定程度上克服阻抗法对内存需求大以及计算效率低的缺点。为了使潮流算法得到进一步的完善，数学中求解非线性问题的经典方法一牛顿拉夫逊方法也被引入到了电力系统潮流计算当中，它以导纳矩阵为基础，其方程有直角坐标和极坐标两种形式，在不同的应用情况下各有所长。相对于阻抗法来说，它在保证良好收敛性以及计算精度的前提下，降低了对计算机内存的需求，提高了运算速度。正因如此，牛顿一拉夫逊法至今仍然是使用最为广泛、效果最好的一种潮流计算方法，也是目前所有潮流计算机算法中最为成熟的一种方法。此后，由牛顿拉夫逊法的极坐标形式经过一定的简化和改进而得到的P Q 分解法(又称改进牛顿法)，也是一种性能比较优越的潮流计
26、算方法，它根据电力系统的特点，抓住主要矛盾，以有功功率误差作为修正电压相角的依据，以无功功率误差作为修正电压幅值的依据，使有功功率和无功功率迭代分开进行，不但降低4山东大学硕十学位论文了修正方程组的阶数，而且使雅可比矩阵的元素在整个迭代过程中维持常数，不必在每次迭代时重新求解，因而在运算速度方面较以前的潮流算法有了很大的突破。由于速度上的明显优势，P Q 分解法还可以用于在线计算【2】。目前，常规潮流算法仍然大量地应用于电力系统各个领域，但由于其模型过于简单，不能全面考虑系统运行时多方面的实际情况，同时选择不同的发电机节点作为平衡节点亦会使所得的潮流结果存在差异，因而在一些特殊的场合以及特定需
27、求下，产生了以常规潮流为基础，而又在某些方面具有特殊功能的其它潮流算法。1 4动态潮流算法的提出常规潮流截取某一时间断面进行计算，其前提是假设系统中功率绝对平衡，全部发电机的输出功率正好等于所有负荷功率与网损之和。然而，实际的电力系统是一个动态的系统，各处的负荷时刻都在发生变化，为了达到供需平衡，系统中发电机的有功输出总体上跟随负荷的变化而变化。在电力系统中，供需恰巧平衡，不存在不平衡功率，频率不发生变化的情况是极为罕见的。通常情况下，都是供需大体平衡，系统存在着不平衡(有功)功率，这将导致系统频率发生变化。本文试图在常规潮流计算的基础上引进频率计算模块，通过发电机和负荷的一次调频来动态分配电
28、网的功率扰动所产生的不平衡功率，改造后的潮流算法称之为动态潮流算法。动态潮流主要是在常规潮流计算的基础上考虑了负荷和发电机的频率动态特性，其核心是潮流计算和频率计算。在动态潮流计算中，系统中由于功率扰动(切负荷、发电机增减出力)而产生的不平衡功率按照各发电机和负荷的功频静特性系数在多台发电机及负荷之间进行分配，得到调整后的发电机出力和负荷的大小以及系统的频率连续变化的情况，这完全克服了常规潮流算法中由于平衡节点选取的差异而导致潮流结果不同的情况。显然，相对于常规潮流算法来说，动态潮流算法能够在一段时间范围内对系统运行的实际情况进行更有效的模拟，是一个较大的进步。程序为了进一步仿真电力系统的扰动
29、，还添加了拓扑分析模块，模仿大电网由于支路开关的断合而导致的系统解列、并网的过程，并能对各个子系统进行潮流计算和频率计算，通过对小系统分析处理完成子网的并网。山东大学硕士学位论文第2 章潮流计算的数学模型应用电子计算机对电力系统进行分析计算时，需要掌握电力系统的数学模型，计算方法和程序设计三个方面的知识，在这一章里我们将介绍潮流计算的数学模型和计算方法。电力系统的数学模型是对电力系统运行状态的一种数学描述。通过数学模型可以把电力系统中物理现象的分析归结为某种形式的数学问题。2 1节点网络方程式电力网络的运行状态可用节点方程来描述，节点方程以母线电压为待求量，母线电压能唯一地确定网络的运行状态。
30、知道了母线电压，就可以算出母线功率、支路功率和电流。电力系统计算通常采用节点方程。在图2-1(a)的简单电力系统中，若略去变压器励磁功率和线路电容，负荷用阻抗表示，便可得到一个有5 个节点(包括零电位点)和7 条支路的等值网络，如图2-l(b)所示。将接于节点1 和4 的电势源和阻抗的串联组合变换成等值的电流源和导纳的并联组合，便得到图2-l(c)所示的等值网络，其中i。=y。毛和|4=y 4 0 I!。分别称为节点1 和4 的注入电流源。24(a)O246山东大学硕士学位论文Y 1 2(V 2-V I)+Y e o V 2+y 2 3(V2 一V3)+y 2 4(V2-V4)=0y 2 3(
31、V3 一V2)+y 3 4(V3 一V 4 户0y 2 4(V 4 一V2)+y 3 4(V 4-V3)+y 4 0 V 4-1 4上述方程组经过整理可以写成：y l l V I+h 2 V j=I，H 2 审j+玛2 也+b 3 审j+，幺审j=0兄3+b 3 v，+蚝4 一0K 2 V Y 4 3 V 3+“V 4=1 4式中Y l l=Y l o+Y 1 2；Y 2 fY 2 0+Y 2 3+y 2 4+Y 1 2：Y 3 3=y 2 3+y 3 4；Y 4 4=y 4 0+Y 2 4+Y 3 4；Y I fY 2 1 2 Y 1 2；Y 2 3 一Y 3 2Y 2 3；Y 矿Y 4 产
32、-m；Y 矿Y 4 3=一y 3 4：一般地，对于有n 个独立节点的网络，可以列写n 个节点方程Y l l V I+Y 1 2 V 2+Y l n V n=IlY j l V l+Y 2 2 V 2+Y 2 n V n=I2Y n l V I+Y n 2 V 寸+Y n n v n=I n也可以用矩阵写成巧，巧：X。E。E。瓦：KK：K厶j l：(2 2)(2 2)(2 3)(2 4)或缩写为讨=t矩阵Y 称为导钠矩阵。导纳矩阵的形成可以归纳以下几点：1)导纳矩阵的阶数等于电力电力网络的节点数。2)导纳矩阵各行非零非对角元个数等于对应节点所连接的不接地支路数。3)它的对角线元素Y i i 称为
33、节点的自导纳，其值等于接于节f 的所有支路导纳之和。4)非对角线元素j 称为节点i、j 间的互导纳，它等于直接接于节点i、j 间支路导纳的负值，若节点i、J 间不存在直接支路，则有M j _ O(由此可知节点导纳矩7山东大学硕士学位论文阵是一个稀疏的对称矩阵)。按照以上原则，则无论电力网络如何复杂，都可以根据给定的输电线路参数和接线拓扑直接求出导纳矩阵。对含变压器的支路，根据丌型等值电路。1，可以求出节点p、q 的自导纳和互导纳分别为：Y p p=l k z+(k-1)I 肛-l zY q q=l k z+(1 助肥F l 七乏(2 5)r p q r q pl k z电力网络通常是由相应的节
34、点导纳矩阵来描述的。在现代电力系统分析中，我们需要面对成千上万个节点及电力网络所连接的电力系统。对电力网络的描述和处理往往成为解决有关问题的关键。电力网络PZl：kq仁 叫驴2-2 变压器支路等值电路的导纳矩阵具有良好的稀疏特性，可以用来高效处理电力网络方程，是现代电力系统分析中广泛应用的数学模型。2 2电力网络方程的求解方法2 2 1 用高斯消去法解网络方程目前电力网络方程主要用高斯消去法求解口l。高斯消去法求解线性方程组由消去运算和回代运算两部分组成。消去运算又叫前代运算，可按行也可按列进行，同样回代运算也可按行或列进行。通常采用“消去按列，回代按行”的方式进行。设有刀阶线性方程组4 x
35、釜B。其中矩阵彳和向量矗的元素可以是实数也可以是复数。由于消去运算只对彳和占进行，因此可以把丑作为第肿l 列附在彳之后，形成n X(时1)阶增广矩阵互=4 曰】=q Iq 2 q。6 l吃l 口2 2 口2。6 21 2 吒q Ia 1 2q Hq m+la 2 1 口巴 a 2 m“1 2为了讨论方便就用吩一1 替代了巧(，=l，2，)。步骤如下：(2 6)，“J首先讨论按列消去过程，它的运算第一步消去第一列。首先把增广矩阵的第-N 规格化为：1a 1 2 1 口1 3 a 1 开+l 1(2-7)式中a l j t=口I J a I l0=2，3,-，I l+1)然后用式(2 7)所表示的
36、行消去j 的第一列对角线下各元素a 2 la 3 1 Q n l 元素，结果使j 的第2 n 行其它元素化为口1 5 呀口f，a o(1(，=2，3，九+l；i=-2，3，功式中：上标1 表示该元素第一次运算的结果。这时矩阵j 变为互：山东大学硕士学位论文五=【4 马】=a t Iq 2(na l (1)q+l 1 色2(n口2 n(1)啦埘+l(”2(1 l a n n(1)a n 月+1 0 1与之对应的方程组是一l 苍啦l，它与A X=B 同解。矩阵未标出的元素为零。第二步，消去第二列，步骤同上。一般地，在消去第k 列时要做以下的运算硝毋=a 勾k-1“柚0=k+l，妒1)(2 8)口妒
37、铆，M 吲1 口(j=k+l,-,n+l；i=k+l，n)、(2 9)经过对矩阵j 的n 次消去运算，即k 从l 依次取到n 按式(2 8)、(2 9)运算使矩阵A 对角线以下的元素全部转化为零，从而得到增广矩阵互=【4 最】_1a 1 2(n口l (1 q H+1(1)1啦H(2 口2 m 十l(2)1a n。，与之相对应的方程组是A B X=B。，即五 l a 1 2(1)x l+a 1 3(1)x 3+q 月(1)矗=q 坩“(1)恐+(2)x 3+a 2 n(2)X n=a 2 J+l(2)毛一l+一l(”q)=一I m+I(”一1)=a n+I(”)(2 1 0)(2 1 1)它与A
38、 X=B 同解。现在来讨论按行回代过程。对于方程组(2 1 1)回代运算自下而上进行。首先由第n 个方程可知X n=a n,州俐。然后将粕代入第疗1 个方程，解出X n J=a n 1,n+l(n a n-l,n 似力h再将X n-I 和而代入第n-2 个方程可解出而o。一般地把已求出的而小船2 而代入第i个方程，即可求出x i=a i,。+1 似一Z a 哥(f=拧，2，1)(2-1 2)这就是回代的一般公式。2 2 2 利用因子表法在实际计算中，常常遇到这种情况：对于方程组A X=B 需要多次求解，每次仅改变其常数项曰，而系数矩阵彳通常是不变的。这时，为了提高计算速度，可以利用因子表对线形
39、方程组求解。因子表可以理解为高斯消去法解线性方程过程中对常数项B 全部运算的一种记录表格。高斯消去法分为消去过程和回代过程。回代过程的运算由对系数矩阵进行消去后得到的上三角矩阵元素确定。见式(2 1 0)。为了对常数项进行消去运算(又称回代运9山东大学硕士学位论文算)，还必须记录消去过程运算所需要的运算因子。消去过程又分为规格化运算和消去运算。由式(2-8)、(2-9)可知，消去过程对常数项B 中的第i 个元素B i 的运算包括：6 p=6“)a i j O I(i _ 1，2 n)(2 1 3)6。)-b k。1)口j 产。酣”(k=1，2 i 1)(2 1 4)将上式中的运算因子嘶，n，“
40、)a L ，(以以及1 口户1 逐行放在下三角部分和式(2-1 0)的上三角矩阵元素合在一起，就得到因子表：1 q 2 na t 3(”q 4(q 月(1)a l t呸-专舻砰“a 2 n(2)码ta 3 2 0)矿1(3)“a 3 n O)甜1 7 4 3(2)专轳口J 舻a n 3(2)群)士(2-1 5)其中下三角及对角线元素可用来对常数项口进行(消去)前代运算，上三角用来进行回代运算。因子表也可写成如下的形式：U。U 2。以。U。厶t 厶：厶，厶巩(2-1 6)其中，1巩2 a f t O-1)=(f ，)岛=沪1(，f)不难看出，因子表式(4 9)中下三角部分的元素就是系数矩阵在消去
41、过程中曾出现的元素，因此只要把它们保留在原来的位置，并把对角线取倒数就可以得到因子表的下三角部分。而因子表的上三角部分的元素就是系数矩阵在消去完成后的结果。对于方程组，需要多次求解，每次仅改变其常数项占而系数矩阵A 是不变的情况，应首先对其系数矩阵A 进行消去运算，形成因子表。有了因子表，就可以对不同的常数1 0如巩kkn 易厶厶山东大学硕士学位论文项B 求解。这时可以直接应用因子表中的元素，用下面的公式代替式(2 1 3)、(2 1 4)进行消去运算：6)-b j(“1)甜。1(i=1，2 n)(2 1 7)6 p)-6 产n 厶也(k=l，2 i-I)(2 1 8)用以下公式代替式进行回代
42、运算：珩；x i=b f 一u t j x j(i=l l ，2，1)(2 1 9)2 3潮流计算的定解条件图2 3 简单电力系统及其等值电路电力系统由发动机、变压器、输电线路及负荷等构成。如上图2-3 表示2 1 节的简单电力系统的等值电路，其网络方程为t=(砖+吃+圪+匕吃)(扛1，2，3，4)节点电流可以用节点功率和电压表示：，：善：一S G i-S l l-：!生=戈2 二!丝二鱼垡2KKK把(2 2 1)代入(2 2 0)可得：(2 2 0)(2 2 1)量墨L 二!鱼立三；!1 2 趔：巧吆-I-r：吃+鬈。吃o：1，2，3)(2 2 2)，这是一组复数方程式，而且是对于V 的非线
43、形方程，如果把实部和虚部分开便得到6个实数方程，但是每一个节点都有6 个变量：发电机发出的有功功率和无功功率，负荷需要的有功功率和无功功率，以及节点电压的幅值和相位(或对应于某一选定参考直角坐标的实部和虚部)。对于1 1 个节点的网络，可列写2 n 个方程，但是有6 n 个变量。通常把发电机与负荷功率作为已知量，并把节点注入功率P，=P c r P L D f 和Q i=Q o，Q 切引入网络方程，就成为4 n 个变量。这样，n 个节点电力系统的潮流方程一般形式是：山东大学硕士学位论文e,-尥：I：窆巧以(f：1，2，n)J=l或只+_，Q=矿毛吃(i=1，2，)J-I(2 2 3)(2 2
44、4)将上述方程式的实部和虚部分开，对每一个节点可得两个实数方程，但是变量仍有4个，既P、Q、K 占。我们必须给定其中的2 个，而留下两个作为待求变量，方程组式可以求解。根据电力系统运行条件，按给定变量的不同一般将节点分为以下三种类型。1 P Q 节点：这类节点有功功率P 和无功功率Q 给定，节点电压(以J)是待定量，通常变电所都是这一种类型的节点，由于没有发电设备，其发电功率为零。在一些情况下，系统中某些发电厂送出的功率在一定时间内为固定时，该发电厂母线也作为P Q 节点，因此，电力系统的大多数节点属于P Q 节点。2 P V 节点：节点的P、V 给定，Q，占待求，这类节点必须得有足够的可调无
45、功功率，用以维持给定的电压幅值，因此又称为电压控制节点。一般是选择有一定无功储备的发电厂和具有可调无功电源设备的变电所作为P V 节点。在电力系统中这一类节点很少。3 平衡节点：节点的几万给定，P、Q 待求。在潮流分布算出以前，网络中的功率损耗是末知的。因此，网络中至少有一个节点的功率不能给定，这个节点承担了系统的功率平衡，故称之为平衡节点旧。另外必须选定一节点，其电压相位均为零，作为各节点电压的参考，这个节点称之为基准节点(其电压幅值给定)。为了计算方便，常将平衡节点和基准节点选为同一节点，可称之为平衡节点，平衡节点只有一个。它的电压幅值和相位已经给定，而其有功功率和无功功率待求。一般选择调
46、频发电厂为平衡节点比较合理，但在进行潮流计算时也可按照别的原则来选择。从以上的讨论可以看到，尽管网络方程是线形方程，但是由于在定解条件中不能给定节点电流，只能给出节点功率，这就使潮流方程变为非线形方程了。由于平衡节点的电压已经给定，所以平衡节点不参加求解。1 2山东大学硕士学位论文第3 章P-0 分解法的基本潮流算法3 1牛顿一拉夫逊法的基本原理首先介绍牛顿法，这是解非线性方程式的有效方法。这个方法把非线性方程式的求解过程变成反复对相应的线性方程式的求解过程，通常称为逐次线性化过程。设有单变量非线性方程m)=O(3 1)求解此方程时，先给出解的近似值x(o)，它与真解的误差为A x(o)，则将
47、满足方程(3 1)，即_，酽+“o 佃(3-2)将式(3 2)左边的函数在x(o 附近展成泰勒级数，于是便得刑+耐0)=删o)+厂酽+厂铲)尘孚+町一)尘笋+(3-3)式中，I x(o)，贝功p o)分别为函数在x(o 处的一阶导数，n 阶导数。如果差值x(o 很小，式(3 3)右端的二次及以上阶次的各项均可略去。于是，式(3 3)可简化为f(x o+A x o)=，(x o)+，(x o)A x o=0(3-4)这是关于修正量o 的线性方程，亦称为修正方程式。解此方程式可得修正量舻一怨，s，用所求得的0 去修正近似解，便得b o)舻犁)一怨由于式(3 4)是略去了高次项的简化式，因此7所解出
48、的修正量也只是近似值。修正后的近似解0 1 与真解仍然有误差。但是，这样的迭代计算可以反复进行下去，迭代计算的通式是+1)“虬爹器 6，迭代过程的收敛判据为l f(x)I 毛(3-7)D或I A x I 岛(3-8)图3-1 牛顿法的几何解释山东大学硕士学位论文式中的8 1，9 2 是预先给定的小正数。这种解法的几何意义可以从m(3 1)中得到证明。函数y=f(x)为图中的曲线，坟x)=O的解相当于曲线与X 轴的交点。如果第k 次迭代中得到x(k)，则过【x，y(k)=f(x 啦瑚点做一切线，此切线同x 轴的交点便确定了下一个近似解x(k+1)。由此可见，牛顿拉夫逊法实质上就是切线法，是一种逐
49、步线性化的方法。牛顿法不仅用于求解单变量方程，也是求解多变量非线性方程的有效方法。设有n 个联立的非线性代数方程：石(，而，毛)=0五(而，X 2，j)2o(3-9)：五(而，恐，x o)=0应用牛顿法求解多变量非线性方程组(3 9)时，假定已给出各变量的初值x 1(o)，X 2(m，而(0)，令A x 尸，彩0)，粕(0 分别为各变量的修正量，使其满足方程，即f,(x f o+日”，拶+“，掣+坪)=0正(+酬”，吐+“，掣+坪)=01Z(研+：0”，霉+谬，王(0+拶)=0(3 1 0)将上式中的n 个多元函数在初始值附近分别展成泰勒级数，并略去含有缸p，血夕)，(o 的二次及以上阶次的各
50、项，使得石(”，掣，x(O)+正(墨”，拶，掣)+；z(妒拶，)+詈I o 酬+薏I。衅)+剖。衅)_。l oI ol o方程式可以写成矩阵形式1 4氯盟8 x 2|o 钆I oI o烈烈；鼍I。篆l o，(0)，缸m，(3 1 1)(3 1 2)线性方程组，称为牛顿法的修正=甜衅盟钒盟钆斗一+峭峭钎一篮+科科盟甄鲠一钆衅峭；衅川H川川J刈均盟魄欲一钆；盟瓯一h卅弘)、，J蚋帕蚋删拶谔；坤尉mOO)钟“皤北lllll【式石五Z0稚山东大学硕士学位论文方程式。利用高斯消去法或三角分解法可以解出修正量缸(0)，x 2(，缸夕)。然后对初始近似解进行修正0 1)-0 0)+彰们0=1，2，曲(3 1