电力系统动态潮流计算及网络拓扑分析

1、分类号：TM734单位代码：10422密级：学号：200413208论文题目：电力系统动态潮流计算及网络拓扑分析作者姓名张国衡专业电路与系统指导教师姓名专业技术职务王良副教授2007年5月15日目录摘要1Abstract2第1章绪论31.1 课题背景31.2 潮流计算的基本要求和要点31.3 潮流计算程序的发展41.4 动态潮流算法的提出5第2章潮流计算的数学模型61 节点网络方程式61 电力网络方程的求解方法81 潮流计算的定解条件11第3章P-Q分解法的基本潮流算法13(1) 牛顿一拉夫逊法的基本原理13(1) 极坐标下的牛顿-拉夫逊法潮流计算15(1) P-Q分解法的原理18(1) P-

2、Q分解法的特点20(1) P-Q分解法的潮流计算步骤21第4章基于电网频率计算的动态潮流221 电力系统的频率特性和一次调频231 频率计算271 微分方程的求解281 频率计算和潮流计算的联合30第5章基于面向对象的动态潮流程序32面向对象的编程思想32对象模型的建立32类的处理和实现34生成应用程序40算例分析42一次调频的手工算例46结论48第6章电力系统的网络拓扑分析49离线数据准备49网络拓扑分析50电网拓扑分析的例题53拓扑分析和潮流计算的接口56第7章动态潮流综合算例分析57程序流程图57I型考题综合算例59华北电网综合算例63结束语65参考文献66附录67致谢78攻读硕士学位期

3、间发表的学术论文79电力系统动态潮流计算及网络拓扑分析摘要电力系统潮流计算是电力系统规划设计与运行分析的基本工具。通过几十年的发展，潮流算法日趋成熟。但由于电网的复杂性，传统的潮流算法依然存在着一些方面的局限性。在电网各种运行方式中，节点注入功率的改变特别是节点注入停运将使系统节点有功、无功注入发生较大的变化，使系统功率出现严重不平衡，而使以往的潮流计算方法在计算这种情况下往往会出现收敛性差、计算结果与实际不符的情况。本文提出的动态潮流算法，主要是在常规潮流计算的基础上考虑了负荷和发电机的静态频率特性，其核心是潮流计算和频率计算。在动态潮流计算中，系统中由于功率扰动（切负荷、发电机增减出力）而

4、产生的不平衡功率按照各发电机和负荷的功频静特性系数在多台发电机及负荷之间进行分配，得到调整后的发电机出力和负荷的大小以及系统的频率连续变化的情况，这克服了常规潮流算法中由平衡节点独自承担不平衡功率而导致潮流收敛性差、结果和实际不符的情况。显然，相对于常规潮流算法来说，动态潮流算法能够对系统运行的实际情况进行更有效的模拟，是一个较大的进步。电力系统在正常运行情况下，可用电力系统状态参数之间的代数方程组来描述其某个特定运行状态。包括潮流计算的许多程序都是以导纳矩阵为基础的。而结点导纳矩阵又是随网络拓扑而变，若不能及时而准确地随着开关所处状态实时变化而修改网络拓扑结构，就会造成分析计算结果的错误。本

5、文将电力系统网络拓扑分析和电力系统的动态潮流两个问题衔结起来，构成一个整体。可实时地随着开关等电气元件信息的变化，进行动态潮流计算，这在工程实际中也有非常广泛的应用。论文还通过典型算例系统的仿真，分析了动态潮流算法的合理性，说明了新算法为准确分析电力系统随负荷变化后的潮流状态提供了一种比较理想的工具。并应用动态潮流和网络拓扑分析成功的模拟了电网的解列和并网操作。关键词：动态潮流；P-Q分解法；稀疏矩阵；节点导纳矩阵；因子表；频率静特性；网络拓扑分析InvestigationonDynamicLoadFlowandNetworkTopologicalAnalysisAbstractPowerfl

6、owcalculationisabasictoolforsystemprogrammingandoperationalanalyzinginpowersystem.Thealgorithmhasbeenputforwardbythedeepresearchofresearchersinandabroadforseveraldecadeyears.Butbecauseofthecomplexityofpowersystem,moreorlessthetraditionalalgorithmshavesomecertainlimitations.Thetraditionalloadflowisn&

7、#39;teasytoconvergewhenalargeunbalancepowerexistswhichcoursedbythechangeofnodeinjectionpower.Andtheresultdoesn'tconformtothereality.Thedynamicloadflowmodelbeputforwardinthispaperistohavebeenconsideredthefrequencystaticcharacteristicofloadandgeneratoronthefoundationoftraditionalloadflow,thecoreis

8、flowcalculationandfrequencycalculation.Indynamicloadflowcalculationtheunbalancepowerwhichcoursedbythechangeofnodeinjectionpowerisdistributedaccordingtothefrequencystaticcharacteristicofeachgeneratorbetweenallgenerators,andgetsthechangeofsystematicfrequencyfollowloadaswellasgeneratorsactivepower.So,i

9、thassurmountedtheconditionthatflowconvergencetobenotgoodandresultisinconsistentwithrealitywhichsincebalancednodeundertakesunbalancedpoweraloneintraditionalloadflow.Obviously,relativetotraditionalloadflow,dynamicloadflowalgorithmcangoonfortheactualconditionthatsystemruns,ismoreeffectiveandsimulated,i

10、tisagreatadvance.Alotofpowerprogramsthatincludeflowcalculationarebasicwithbusadmittancematrix.Andbusadmittancematrixisalsotomakerubbingwiththechangeofnetworktopologicalstructure.Ifwecannotmodifynetworktopologicalstructureintimeaccuratelyalongwiththeswitchlocatedstatechangeofrealtimetomakerubbingwill

11、causeresultmistake.Thethesislinksupdynamicloadflowanalysisandnetworktopologicalanalysisinpowersystemasawholetotakedynamicloadflowcalculationofrealtime.Thisprojecthasbroadapplicationinpractice.Thethesisalsovalidatestherationalityofthedynamicpowerflowalgorithmthroughemulatingonatypicalsystem.Theresult

12、indicatesthatthenewalgorithmhasprovidedamoreidealtoolforanalyzingthepowerflowstateofpowersystemafterthechangeofsystemloadswellandtrulyFinallywehavesimulateduntierowandparalleloperationofelectricalnetworksuccessfullyaccordingtothetheoryofdynamicloadflowandnetworktopologicalanalysis.Keywords:Dynamiclo

13、adflow,P-Qdecompositionmethod,Sparsematrix,Busadmittancematrix,Factortable,Frequencystaticcharacteristic;Networktopologicalanalysis第1章绪论课题背景当今，电能以其清洁、高效、便于输送等突出优点，己经成为全球广泛使用的最主要能源，在社会经济的发展中起到举足轻重的作用。电力系统是当今世界最庞大的人工系统，它包括发电、输电、配电、用电四个环节，对其运行的最基本要求有三点：保证安全可靠的供电、要有合乎要求的电能质量、要有良好的经济性。随着国民经济的进一步发展，社会各部门对

14、电力系统的发展提出了更新、更高的要求，当今电力系统正向着超高压、大容量、远距离的输电方式发展，全国电网区域互联也成为一种必然的趋势，这样的发展趋势在很大程度上提高了系统运行效率，增加了经济效益，促进了对能源的合理开发和利用，减轻了对环境保护的压力，但同时也给电力系统的安全运行带来了新的问题。经过近几年的技术改造和升级，大多数电网数据采集与监控（SCADA）系统都已通过实用化验收。在此基础上，对电力系统调度自动化的要求也不断提高，能量管理系统（EMS）在调度中心的应用逐步走向实用化。一个实用化的能量管理系统包括：数据库管理系统、人机管理系统、经济调度系统、网络管理系统和高级应用软件系统，五大部分

15、有机地结合在一起，为用户提供服务1。其中，高级应用软件系统是整个能量管理系统的灵魂，只有它的正常工作才可以为用户提供分析与控制服务。高级应用软件系统又由实时网络状态估计、在线潮流、安全分析、最优潮流等可以不断延伸的各种分析与控制功能组成。在这些众多功能中，潮流计算是基本的工具，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态，是电力系统运行、规划以及安全性、可靠性和优化的基础，也是各种暂态分析的基础和出发点。现代电力系统规模庞大，结构复杂，其规划、设计和运行均需借助电子计算机进行潮流计算分析。近几十年来，随着电力系统规模的不断扩大与结构的日益复杂、电力系统自动化水平的提高和计

16、算机技术的日新月异，电力系统潮流的计算机算法也在不断地进步与更新。但是，从计算效率、收敛性以及对实际系统的有效模拟程度等各个方面综合起来看，现有的诸多潮流算法仍然存在不少尚待解决的问题。所以，还需对电力系统潮流进行更深入的研究，发展更加完善的潮流算法，以满足处于不断扩大与更新中的当代电力系统的需求。潮流计算的基本要求和要点潮流计算随计算性质不同而有不同的要求，如长距离输电、区域性网络、城市配电网络等都有不尽相同的要求，但仍有其共同的基本要求。首先是不同类型的网络在各种运行方式下，网络各节点的电压水平应符合有关规定。其次如网络中各线路的潮流分布不应有线路过载等。对潮流计算的分析主要根据计算的目的

17、而定。在电力系统运行方式中一般含高峰负荷和低谷负荷时运行方式下，在具有水力发电厂的电力系统中根据水电厂水文特点又有丰水期、平水期、枯水期的运行方式，止匕外，也需要研究事故运行方式和各种特殊运行方式。在潮流计算中首先应效验网络枢纽点的电压水平及网络各节点的电压是否满足要求，其次效验各发电厂发电机的有功及无功出力是否符合技术要求，另外根据计算的要求对各线路、变压器的潮流进行分析。潮流计算程序的发展电力系统潮流计算是研究电力系统稳态运行情况的一项基本运算，其数学本质是一组多元非线性方程，主要采用迭代的方法求解。电力系统潮流计算从提出至今，经历了一个由手工，利用交、直流计算台到应用数字电子计算机的发展

18、过程，现有的潮流算法都以计算机的应用为前提。由于潮流计算在电力系统分析研究中具有重要的地位，吸引了大量的专家学者对其进行了研究，针对各种实际情况以及特殊需求，发展了多种用于电力系统潮流计算的计算机算法。在现有的潮流算法之中，最早出现的是常规潮流算法，其它潮流算法都是根据不同的实际需求在常规潮流的基础上发展起来的。利用电子计算机进行电力系统潮流计算始于上个世纪50年代中期，最初以节点导纳矩阵为基础进行迭代，其原理简单，易于编程实现，同时由于导纳矩阵是稀疏矩阵，对计算机内存需求不大。但是此法收敛性较差，其迭代次数会随着系统规模的扩大而急剧增加，易出现不收敛的情况。在这种情况下，出现了基于阻抗矩阵的

19、迭代方法，大大改善了潮流计算的收敛性，可以求解一些用导纳法无法收敛的潮流问题。但是，阻抗矩阵是满秩矩阵，不但占用内存大，而且每次迭代所需的计算量也比较大，这就引入了新的问题。随后出现的分块阻抗法，它将一个大系统分为若干小系统，只需存储各个小系统的阻抗矩阵以及它们之间联络线的阻抗，此法能够在一定程度上克服阻抗法对内存需求大以及计算效率低的缺点。为了使潮流算法得到进一步的完善，数学中求解非线性问题的经典方法一牛顿-拉夫逊方法也被引入到了电力系统潮流计算当中，它以导纳矩阵为基础，其方程有直角坐标和极坐标两种形式，在不同的应用情况下各有所长。相对于阻抗法来说，它在保证良好收敛性以及计算精度的前提下，降

20、低了对计算机内存的需求，提高了运算速度。正因如此，牛顿一拉夫逊法至今仍然是使用最为广泛、效果最好的一种潮流计算方法，也是目前所有潮流计算机算法中最为成熟的一种方法。此后，由牛顿-拉夫逊法的极坐标形式经过一定的简化和改进而得到的PQ分解法（又称改进牛顿法），也是一种性能比较优越的潮流计算方法，它根据电力系统的特点，抓住主要矛盾，以有功功率误差作为修正电压相角的依据，以无功功率误差作为修正电压幅值的依据，使有功功率和无功功率迭代分开进行，不但降低了修正方程组的阶数，而且使雅可比矩阵的元素在整个迭代过程中维持常数，不必在每次迭代时重新求解，因而在运算速度方面较以前的潮流算法有了很大的突破。由于速度上

21、的明显优势，PQ分解法还可以用于在线计算2。目前，常规潮流算法仍然大量地应用于电力系统各个领域，但由于其模型过于简单,不能全面考虑系统运行时多方面的实际情况，同时选择不同的发电机节点作为平衡节点亦会使所得的潮流结果存在差异，因而在一些特殊的场合以及特定需求下，产生了以常规潮流为基础，而又在某些方面具有特殊功能的其它潮流算法。动态潮流算法的提出常规潮流截取某一时间断面进行计算，其前提是假设系统中功率绝对平衡，全部发电机的输出功率正好等于所有负荷功率与网损之和。然而，实际的电力系统是一个动态的系统，各处的负荷时刻都在发生变化，为了达到供需平衡，系统中发电机的有功输出总体上跟随负荷的变化而变化。在电

22、力系统中，供需恰巧平衡，不存在不平衡功率，频率不发生变化的情况是极为罕见的。通常情况下，都是供需大体平衡，系统存在着不平衡（有功）功率，这将导致系统频率发生变化。本文试图在常规潮流计算的基础上引进频率计算模块，通过发电机和负荷的一次调频来动态分配电网的功率扰动所产生的不平衡功率，改造后的潮流算法称之为动态潮流算法。动态潮流主要是在常规潮流计算的基础上考虑了负荷和发电机的频率动态特性，具核心是潮流计算和频率计算。在动态潮流计算中，系统中由于功率扰动（切负荷、发电机增减出力）而产生的不平衡功率按照各发电机和负荷的功频静特性系数在多台发电机及负荷之间进行分配，得到调整后的发电机出力和负荷的大小以及系

23、统的频率连续变化的情况，这完全克服了常规潮流算法中由于平衡节点选取的差异而导致潮流结果不同的情况。显然，相对于常规潮流算法来说，动态潮流算法能够在一段时间范围内对系统运行的实际情况进行更有效的模拟，是一个较大的进步。程序为了进一步仿真电力系统的扰动，还添加了拓扑分析模块，模仿大电网由于支路开关的断合而导致的系统解列、并网的过程，并能对各个子系统进行潮流计算和频率计算，通过对小系统分析处理完成子网的并网。第2章潮流计算的数学模型应用电子计算机对电力系统进行分析计算时，需要掌握电力系统的数学模型，计算方法和程序设计三个方面的知识，在这一章里我们将介绍潮流计算的数学模型和计算方法。电力系统的数学模型

24、是对电力系统运行状态的一种数学描述。通过数学模型可以把电力系统中物理现象的分析归结为某种形式的数学问题。2.1节点网络方程式电力网络的运行状态可用节点方程来描述，节点方程以母线电压为待求量，母线电压能唯一地确定网络的运行状态。知道了母线电压，就可以算出母线功率、支路功率和电流。电力系统计算通常采用节点方程。在图2-1(a)的简单电力系统中，若略去变压器励磁功率和线路电容，负荷用阻抗表示，便可得到一个有5个节点(包括零电位点)和7条支路的等值网络，如图2-1(b)所示。将接于节点1和4的电势源和阻抗的串联组合变换成等值的电流源和导纳的并联组合，便得到图2-1(c)所示的等值网络，其中&二

25、y10合和除y40层4分别称为节点1和4的注入电流源。2410-R3D1广一0F1(a)以零电位作为计算节点电压的参考点，根据基尔霍夫电流定理，可以写出4个独立节点的电流平衡方程如下:y10V1+y12(V1-V2)=I1yi2（V2-V1）+y20V2+y23（V2-V3）+y24（V2-V4）=0y23（V3-V2）+y34（V3-V4）=0（2-2）y24（V4-V2）+y34（V4-V3）+y40V4=I4上述方程组经过整理可以写成：YiiVi+Y12V2=I1Yl2Vl+Y22V2+Y23V3+Y24V4=0Y23V2+Y33V3+Y34V4=0（2-2）Y42V2+Y43V3+Y4

26、4V4=I4式中Yii=yi0+yi2；Y22=y20+y23+y24+yi2；丫33=y23+y34；Y44=y40+y24+y34；Yi2=Y2i=-yi2；Y23=Y32=-y23；Y24=Y42=-y24；Y34=Y43=-y34；一般地，对于有n个独立节点的网络，可以列写n个节点方程YiiVi+Y12V2+YinVn=Ii+YnnVn=In（2-4）Y2iVi+Y22V2+Y2nVn=I2（2-3）也可以用矩阵写成?YiY12LYinVi?Ii丫21丫22L，nV2?I2MMMMMYniYn2LYnn?Vn?InYniVi+Yn2V2+或缩写为YV&&矩阵Y称为导钠矩

27、阵。导纳矩阵的形成可以归纳以下几点：1）导纳矩阵的阶数等于电力电力网络的节点数。2）导纳矩阵各行非零非对角元个数等于对应节点所连接的不接地支路数。3）它的对角线元素Yii称为节点的自导纳，其值等于接于节i的所有支路导纳之和。4）非对角线元素Yij称为节点i、j间的互导纳，它等于直接接于节点i、j间支路导纳的负值，若节点i、j间不存在直接支路，则有Yij=0（由此可知节点导纳矩阵是一个稀疏的对称矩阵)。按照以上原则，则无论电力网络如何复杂，都可以根据给定的输电线路参数和接线拓扑直接求出导纳矩阵。对含变压器的支路，根据兀型等值电路网，可以求出节点p、q的自导纳和互导纳2-2变压器支路等值电路分别为

28、：Ypp=1/kz+(k-1)/kz=1/zYqq=1/kz+(1-k)/k2z=1/k2z(2-5)Ypq=Yqp=-1/kz电力网络通常是由相应的节点导纳矩阵来描述的。在现代电力系统分析中，我们需要面对成千上万个节点及电力网络所连接的电力系统。对电力网络的描述和处理往往成为解决有关问题的关键。电力网络的导纳矩阵具有良好的稀疏特性，可以用来高效处理电力网络方程，是现代电力系统分析中广泛应用的数学模型。2.2电力网络方程的求解方法用高斯消去法解网络方程目前电力网络方程主要用高斯消去法求解50高斯消去法求解线性方程组由消去运算和回代运算两部分组成。消去运算又叫前代运算，可按行也可按列进行，同样回

29、代运算也可按行或列进行。通常采用“消去按列，回代按行”的方式进行。设有n阶线性方程组AX=B。其中矩阵A和向量B的元素可以是实数也可以是复数由于消去运算只对A和B进行，因此可以把B作为第n+1列附在A之后，形成nx(n+1)阶增广矩阵Lana21a22La2nAABMMMMan1an2Lannba12Lana1,n1b2a21a22La2na2,n1MMMMMMbnan1an2Lannan,n1(2-6)为了讨论方便就用aj,n+1替代了bj(j=1,2,n)。首先讨论按列消去过程，它的运算步骤如下:第一步消去第一列首先把增广矩阵的第一列规格化为:式中a1ja12a13a1，n+1(2-7)=

30、a1j/an(j=2,3，,n+1)然后用式(2-7)所表示的行消去A的第一列对角线下各元素a21a3-an1元素,结果使A的第2n行其它元素化为aij(1)=aijwa1j(j=2,3，,n+1;i=2,3，,n)式中：上标表示该元素第一次运算的结果。这时矩阵A变为Ai：ABiaii互L都a22(1)La2nMMMMan2Lann(1)ai,na2,n1an,n与之对应的方程组是A1X=B1,它与AX=B同解。矩阵未标出的元素为零第二步，消去第二列，步骤同上一般地，在消去第k列时要做以下的运算akj的=akj(k-1)/akk(k-1)(j=k+1,n+1)(2-8)aij(k)=aij(k

31、-1)-aik(k-1)akj(k)(j=k+1,，n+1;i=k+1,n)(2-9)经过对矩阵A的n次消去运算，即k从1依次取到n按式(2-8)、(2-9)运算使矩阵A对角线以下的元素全部转化为零，从而得到增广矩阵AnAnBna(1)a1na2nMa1a2与之相对应的方程组是AnX=Bn,即anl,n1,n1M(n),n1(2-10)a1a(1),n12,n1(1)a12x1X2(1)a13X3(2)a23X3(1)anXn(2)a2nXnXn1an1(n1)XnXn(n1)an1,n1a(n)an,n1它与AX=B同解。现在来讨论按行回代过程。对于方程组(2-11)回代运算自下而上进行。首

32、先由第n个方程可知Xn=an,n+1。然后将Xn代入第n-1个方程，解出Xn-1=an-1,n+1(n-1)-an-1,n(n-1)Xn再将Xn-1和Xn代入第n-2个方程可解出Xn-2。一般地把已求出的Xi+1,Xi+2,Xn代入第i个方程，即可求出Xi=ai,n+1Wai,j(2-11)Xj(i=n,-,2,1)(2-这就是回代的一般公式。利用因子表法在实际计算中，常常遇到这种情况：对于方程组以=8需要多次求解，每次仅改变其常数项B,而系数矩阵A通常是不变的。这时，为了提高计算速度，可以利用因子表对线形方程组求解。因子表可以理解为高斯消去法解线性方程过程中对常数项B全部运算的一种记录表格。

33、高斯消去法分为消去过程和回代过程。回代过程的运算由对系数矩阵进行消去后得到的上三角矩阵元素确定。见式(2-10)。为了对常数项进行消去运算(又称回代运算)，还必须记录消去过程运算所需要的运算因子。消去过程又分为规格化运算和消去运算。由式(2-8)、(2-9)可知，消去过程对常数项B中的第i个元素Bi的运算包括：bia.(j1)ijij不难看出，因子表式(4-9)中下三角部分的元素就是系数矩阵在消去过程中曾出现的元素，因此只要把它们保留在原来的位置，并把对角线取倒数就可以得到因子表的下三角部分。而因子表的上三角部分的元素就是系数矩阵在消去完成后的结果。对于方程组，需要多次求解，每次仅改变其常数项

34、B而系数矩阵A是不变的情况,应首先对其系数矩阵A进行消去运算，形成因子表。有了因子表，就可以对不同的常数=bi(i-1)/aii(i-1)(i=1,2-.n)(2-13)bi(k)=bi(k-1)-aik(k-1)bk(k)(k=1,2.i-1)(2-14)将上式中的运算因子ai1,ai2,aik(k-1)-ai,i-1(i-2)以及1/aii(i-1)逐行放在下三角部分和式(2-10)的上三角矩阵元素合在一起，就得到因子表：(1)a14LL(1)an1a12a13a11a211c(T)a23a22a24La2na31a32a33a34a3n(2-15)a41a42a43a4na44(1)an

35、1an2an3an41T(n1)ann其中下三角及对角线元素可用来对常数项运算。因子表也可写成如下的形式：B进行(消去)前代运算，上三角用来进行回代(2-16)ij)ji)DuU12U13U14LLU1nL21D22U23U24LLU2nL31L32D33U34LLU3nL41L42L43D44LLU4nLLLLLn1Ln2Ln3L04LLDnn其中d.A.dii(i1)aiiua(i)ijij项B求解。这时可以直接应用因子表中的元素，用下面的公式代替式(2-13)、(2-14)进行消去运算：bi=bi(i-1)/dii(i-1)(i=1,2-.n)(2-17)bi(k)=bi(k-1)-li

36、kbk的(k=1,2.i-1)(2-18)用以下公式代替式进行回代运算：xn=bn(n);xi=biWui,jxj(i=n2,1)(2-19)2.3潮流计算的定解条件图2-3简单电力系统及其等值电路电力系统由发动机、变压器、输电线路及负荷等构成。如上图2-3表示2.1节的简单电力系统的等值电路，其网络方程为由(Yi冷丫2通Y泊丫/)(i1,2,3,4)(2-20)节点电流可以用节点功率和电压表示:(PGiPLDi)j(QGiQLDi)(2-21)把(2-21)代入(2-20)可得:(PGiPLDi)j(QGiQLDi)(i1,2,3)(2-22)这是一组复数方程式，而且是对于V的非线形方程，如

37、果把实部和虚部分开便得到6个实数方程，但是每一个节点都有6个变量：发电机发出的有功功率和无功功率，负荷需要的有功功率和无功功率，以及节点电压的幅值和相位(或对应于某一选定参考直角坐标的实部和虚部)。对于n个节点的网络，可列写2n个方程，但是有6n个变量。通常把发电机与负荷功率作为已知量，并把节点注入功率Pi=PGi-PLDi和Qi=QGi-QLdi引入网络方程，就成为4n个变量。这样，n个节点电力系统的潮流方程一般形式是：n.RjQiV?YjVj(i1,2,L,n)(2-23)j1n或RjQiViY?V?j(i=1,2,，n)(2-24)ji将上述方程式的实部和虚部分开，对每一个节点可得两个实

38、数方程，但是变量仍有4个，既P、Q、V、6。我们必须给定其中的2个，而留下两个作为待求变量，方程组式可以求解。根据电力系统运行条件，按给定变量的不同一般将节点分为以下三种类型。PQ节点：这类节点有功功率P和无功功率Q给定，节点电压(V、6)是待定量，通常变电所都是这一种类型的节点，由于没有发电设备，其发电功率为零。在一些情况下，系统中某些发电厂送出的功率在一定时间内为固定时，该发电厂母线也作为PQ节点，因此,电力系统的大多数节点属于PQ节点。PV节点：节点的P、V给定，Q、6待求，这类节点必须得有足够的可调无功功率，用以维持给定的电压幅值，因此又称为电压控制节点。一般是选择有一定无功储备的发电

39、厂和具有可调无功电源设备的变电所作为PV节点。在电力系统中这一类节点很少。平衡节点：节点的V、6给定，P、Q待求。在潮流分布算出以前，网络中的功率损耗是末知的。因此，网络中至少有一个节点的功率不能给定，这个节点承担了系统的功率平衡，故称之为平衡节点60另外必须选定一节点，其电压相位均为零，作为各节点电压的参考，这个节点称之为基准节点(其电压幅值给定)。为了计算方便，常将平衡节点和基准节点选为同一节点，可称之为平衡节点，平衡节点只有一个。它的电压幅值和相位已经给定，而具有功功率和无功功率待求。一般选择调频发电厂为平衡节点比较合理，但在进行潮流计算时也可按照别的原则来选择。从以上的讨论可以看到，尽

40、管网络方程是线形方程，但是由于在定解条件中不能给定节点电流，只能给出节点功率，这就使潮流方程变为非线形方程了。由于平衡节点的电压已经给定，所以平衡节点不参加求解。第3章P-Q分解法的基本潮流算法牛顿一拉夫逊法的基本原理首先介绍牛顿法，这是解非线性方程式的有效方法。这个方法把非线性方程式的求解过程变成反复对相应的线性方程式的求解过程，通常称为逐次线性化过程。设有单变量非线性方程f(x)=0(3-1)求解此方程时，先给出解的近似值x(0),它与真解的误差为A。)，则将满足方程(3-1),即f(x(0)+及)=0(3-2)将式(3-2)左边的函数在x(0)附近展成泰勒级数，于是使得/(0)x2(0)

41、n(0)(0)(0)(0)、(0)(0)、(xJ(n),(0)、(x)f(xx)f(x)f(x)xf(x)Lf(x)L(3-3)2!n!式中f'(x(0),，f(n)(x(0)分别为函数f(x)在x(0)处的一阶导数，n阶导数。如果差值Ax很小，式(3-3)右端的二次及以上阶次的各项均可略去。于是，式(3-3)可简化为f(x(0)x(0)f(x(0)f(x(0)x(0)0(3-4)这是关于修正量Ax(0)的线性方程，亦称为修正方程式。解此方程式可得修正量式中的日，2是预先给定的小正数。(0)f(x(0)f(x(0)(3-5)用所求得的Ax去修正近似解，使得x(1)x(0)x(0)x(0

42、)f(x(0)f(x(0)由于式(3-4)是略去了高次项的简化式，因此所解出的修正量Ax也只是近似值。修正后的近似解x与真解仍然有误差。但是，这样的迭代计算可以反复进行下去，迭代计算的通式是f(x(k)f(x(k)(3-6)迭代过程的收敛判据为|f(x(k)|1(3-7)或|x(k)|2(3-8)图3-1牛顿法的几何解释这种解法的几何意义可以从图(3-1)中得到证明。函数y=f(x)为图中的曲线，f(x)=0的解相当于曲线与X轴的交点。如果第k次迭代中得到X(k),则过x(k),y(k)=f(x(k)点做一切线，此切线同X轴的交点便确定了下一个近似解x(k+1)0由此可见，牛顿-拉夫逊法实质上

43、就是切线法，是一种逐步线性化的方法。牛顿法不仅用于求解单变量方程，也是求解多变量非线性方程的有效方法。设有n个联立的非线性代数方程：(3-9)f1(X,X2,LXn)f2(X1,X2,LXn)Mfn(X1,X2,LXn)应用牛顿法求解多变量非线性方程组(3-9)时，假定已给出各变量的初值X1(0),X2(0),-，Xnf1(x；0)f2(X1(0),令AX1(0),及2(0),Axn分别为各变量的修正量，使其满足方程，即(0)(0)I(0)X1,X2X2,L,Xn(0)(0)I(0)X1,X2X2,L,XnM(0)(0)(0)I(0)X1,X2X2,XnX0)xn0)xn0)(3-10)将上式

44、中的n个多元函数在初始值附近分别展成泰勒级数，并略去含有AX1(0),及2(0),，加的二次及以上阶次的各项，使得fi(x,x20),l父)fi(0)X1f1X2fiXn(0)Xn(0)(0).(0)、f2(X1,X2,L,Xn)f2K(0)X1f2X2(0)X2(0)XnXn(3-11)fn(Xi(0),x20),L,Xn0)fnfnX2fnXn(0)Xn方程式可以写成矩阵形式f1(x1(0),x20),Lf2(X1(0),x20),LMfn(研,X20),Lxn0)X?)x?)方程式(3-12)是对于修正量f1X1f2X1MfnX1f1XnX2MfnX2AX1(0),Ax2(0),f2Xn

45、MfnXn(3-12)&n的线性方程组，称为牛顿法的修正方程式。利用高斯消去法或三角分解法可以解出修正量Axi(0),及2(0),，Axn(0)o然后对初始近似解进行修正Xi=xi+Axi(i=1,2,n)如此反复迭代，在进行第k次迭代时，求解修正方程式fi(Xi(k),x2k),Lxnk)f2(Xi(k),x2k),Lxnk)Mfn(X(k),x2k),LXnk)flXif2XiMfiMfnXiX2MfnXnk得到修正量Axi的，AX2(k),，AXn(k),并对各变量进行修正Xi(k+i)=Xi(k)+AXi(k)(i=i,2,n)式(3-I4)和(3-I5)也可以缩写为F(X(k

46、)=-J(k)AX(k)和X(k+i)=X(k)+ZX(k)式中，X和似分别是有n个变量和修正量组成的n维列向量;(3-I3)(3-I4)(3-I5)(3-I6)(3-I7)F(x)是由n个多元函数组成的n维列向量；J是nxn阶方阵，称为雅可比矩阵，它的第i、j个元素JjfjXj是第i个函数fi(XI,X2，,Xn)对第j个变量的偏导数，上角标(k)代表示J阵每一个元素都在点(XI吗X2(k),，刈(k)处取值。迭代过程一直进行到满足收敛判据MaX|fi(Xi(k),X2k),LXnk)|<i(3-I8)或MaX|AXi(k)|<2(3-I9)为止，©和2为预先给定的小正

47、数。将牛顿-拉夫逊法进行潮流计算，要求将潮流方程写成形如(3-9)的形式。由于节点电压可以采用不同的坐标系来表示，牛顿-拉夫逊法潮流计算也将相应的采用不同的计算公式。极坐标下的牛顿-拉夫逊法潮流计算采用极坐标时，节点电压表示为V&ViiVi(cosijsinJ,导纳矩阵元素则表示为YjGjjBj将其和导纳矩阵表示式带入节点的功率方程(2-24)右端，展开并分出实部和虚部，变得R=ViVGjcos斗+马sin斗)j=in(3-20)Qi=ViVj(Gijsin即-Bjcos弓)j=i式中，&=a-6是两节点电压的相角差。方程式(3-20)把节点功率表示为节点电压的幅值和相角的函数

48、。在有n个节点的系统中，假定第1m号节点为PQ节点，第m+1n-1号节点为PV节点，第n号节点为平衡节点。在极坐标系中Vn和&是给定的，PV节点的电压幅值Vm+1Vn-1也是给定的。因此，只剩下n-1个节点的电压相角&,金，,吊-1和m个节点的电压幅值5,72,Vm是未知量，一共包含了n-1+m个方程式，正好同未知数的数目相同。实际上，对于每一个PQ节点或每一个PV节点都可以列写一个有功功率不平衡量方程式nPPsPPisViVj(GijCosijBijsinj)0(i=1,2,n-1)(3-21)j1而对于每一个PQ节点还可以在列写一个无功功率不平衡量方程式nQiQisQiQi

49、sViVj(GjsinjBjcosj)0(i=1,2,m)(3-22)j1对于方程式(3-21)和(3-22)可以写出修正方程式如下RH11H12LH1pH1,n1N11N12N1m1P2H21H22LH2P1H2,n1N21N22LN1m2ML；11LLMPpHp1Hp2LHppHp,n1Np1Np2N1mpP,1Hn1,1Hn1,2LHn1,pHn1,n1,Nn1,1Nn1,2Nn1,mn1Q1J11J12LJ1pJ1,n1L11L12LLm711MQ2J21J22J2pJ2,n1-L21L22L2mV2/V2MLLILLMQmJm,1Jm,2Jm,pJm,n1&Lm,1Lm,2L

50、m,mVm/Vm方程式可简写为式中PHNQKLnQ1P2；QQ2MMPn1QmPV1V2MVmV1(3-23)(3-24)H是(n-l)X(n-l)阶方阵,；VdV2其元素为Hp三；N是（n-1）xm阶矩阵，其兀素为旦；L是mXm阶方阵，其元素p一NjVj；K是mX（n-1）阶矩阵，其兀素为KjVjQ,为LjVP"。Vj在这里把节点不平衡功率对节点电压幅值的偏导数都乘以该节点电压，相应地把节点电压的修正量都除以该节点的电压幅值，这样，雅可比矩阵的表达式就具有比较整齐的形式。iwj时:对式（3-21）和式（3-22）求偏导数，可以得到雅可比矩阵元素的表达式如下：当当i=j时:HijP一

51、-VVj(Gjsin°-Bjcosj)HhP-2_iQi吊力NijKij3Qi-VVj(GcosVVj(GjcosjjBjsinj)(3-25)Bjsinij)KiiVTQi-P-V2Gi(3-26)-PVIG.LijQi八Vji-VVj(GjsinVjij-Bjcosj)-QiV2Bh阵。计算步骤如下，首先要输入网络的原始数据以及各节点的给定值并形成节点导纳矩输入节点电压初值Vi和8置迭代计数k=0。然后开始进入牛顿法的迭代过程。在进入第k+1次迭代时其计算步骤如下:（1）按上一次迭代算出的个节点电压，利用式（3-21）和（3-22）计算各类节点的不平衡量APi(k)、AQi(k)

52、和AVi2(k)(2)按条件(3-18)校验收敛，即max|APi的,AQi(k),泉(k)|<。计算雅克比矩阵的各元素(4)解修正方程式求各节点的修正量AVi的和A*)修正各节点的电压Vi(k+1)=Vi的+AVi的，泸+1)=Xk)+A泄(6)迭代计数加一，返回第一步继续迭代过程。3.3P-Q分解法的原理采用极坐标形式表示节点电压，能够根据电力系统实际运行状态的物理特点，对牛顿潮流计算的数学模型进行合理的简化。在交流高压电网中，输电线路的电抗要比电阻大得多，系统中母线有功功率的变化P一?P则主要受母线电压幅值变化的影响。在修正方程式的系数矩阵中,偏导数和一数V?6值是相当小的。作为简

53、化的第一步，可以将方程式(3-24)中的子块N和K略去不计，即认为它们的元素为零。这样，n-1+m阶的方程式(3-24)便分解为一个n-1阶和一个m阶的方程AP=-HA6(3-27)郎-LVd-1AV(3-28)这一简化大大地节省了机器内存和解题时间。方程式(3-27)和(3-28)表明，节点的有功功率不平衡量只用于修正电压的相位，节点的无功功率不平衡量只用于修正电压的幅值。这两组方程轮流迭代，这就是所谓的有功-无功功率分解法。但是矩阵H和L的元素都是节点电压幅值和相角差的函数，其数值在迭代过程中是不断变化的。因此，最关键的一步简化就在于，把系数矩阵H和L简化为常数矩阵。它的根据是什么呢？在一般情况下，线路两端电压的相角差是不大的(不超过10。20。)，因此可以认为cos6产1,Gijsin6ij<<Bij此外，与系统各节点无功功率相适应的导纳BLDi必远小于该节点自导纳的虚部，即BLDi=Qi/Vi2<<Bii和Qi<<Vi2Bii考虑到以上的关系，矩阵H和L的元素的表达