基于MATLAB软件的P-Q分解法潮流计算

1、南京工程学院基于MATLA欧件的P-Q分解法潮流计算摘要电力系统潮流计算是研究电力系统稳态运行情况的一种重要的分析计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性，可靠性和经济性。所以，电力系统潮流计算是进行电力系统故障计算，继电保护整定，安全分析的必要工具。随着电力系统网络的急剧扩大和不断复杂，运用手算进行潮流计算已经不现实。但是，伴随着计算机技术的飞速发展，基于计算机的潮流计算也就应运而生了。这样，通过潮

2、流计算，实现对系统的分析成为可能。本文结合潮流计算的三个基本要求，紧跟该领域的发展，介绍了基于MATLAB软件P-Q分解法潮流计算的程序，该程序用于粗略的计算中小型电力网络的潮流，实现对其的分析。本文所设计的程序，在计算中，所用的算法通俗易懂并对以往的主流算法做了一些改进，提高了计算速度。同时，该程序采用了GUI人机对话，将Excel表格、TXT文档与MATLAB程序紧密联系起来，使输入输出界面更加人性化。关键词：电力系统潮流计算；P-Q分解法；MATLAB软件南京工程学院PowerflowcalculationofP-QmodebasedonMATLABsoftwareAbstractPow

3、erflowcalculationisoneoftheimportantcalculationswhicharetostudytheoperationofpowersystemsteadystateanalysis.Itisbasedonthegivenoperatingconditionsandsystemwiringtoidentifythevariouspartsofthepowersystemoperatingstate:thebusesvoltage,thestreamcomponentspower,systempowerloss.bothpowersystemplanningdes

4、ignandoperationofexistingpowersystemmodeofstudyareneedtousethepowerflowcalculationtoquantitativelycomparetheprogramorrunmodepowersupplyreasonable,reliabilityandeconomy.Therefore,thepowerflowcalculationisanessentialtoolforacalculationofpowersystemfaults,protectionsetting,securityanalysis.withtherapid

5、expansionofpowersystemnetworkandcontinuingtobemorecomplex,usinghandcalculationforflowcalculationhasbeenunrealistic.But,withtheceleritydevelopmentincomputertechnology,computer-basedpowerflowcalculationhasalsoemerged.Inthisway,Itispossibletoanalysispowersystemthroughthepowerflowcalculation.Basedonthet

6、hreebasicrequirementsofpowerflowcalculationandfollowedbythedevelopmentofthefield,ThispaperintroducesthePQmodepowerflowcalculationprocedurebasedonMATLABsoftware.Itisusedforaroughcalculationofthesmallandmediumpowernetworktoachieveitsanalysis.Thealgorithmusedintheprocedurementionedinthispaperismoreeasy

7、tounderstandandmadesomeimprovementstoenhancethecomputing-II-南京工程学院speedratherthanthepast.Atthesametime,theprogramusestheGUImanmachinedialogue.SoExceltable,TXTdocumentsiscloselylinkedwiththeMATLABprogramtomaketheinputandoutputinterfacesmorehumanity.Keywords:powerflowcalculation；P-Qdecompositionmode；M

8、ATLABsoftware-III-南京工程学院目录摘要IAbstractII第1章绪论1.1.1 课题背景1.1.2 电力系统潮流计算 电力系统潮流计算简介 电力系统潮流计算的基本要求21.3 潮流计算的意义及其发展31.4 本次毕业设计主要工作4第2章潮流计算的原理及具体算法过程62.1 电力网络的数学模型.62.1.1 电力网络的基本方程62.1.2 导纳矩阵的形成 电力网络中几种特殊的数学模型82.2 电力系统潮流计算 电力系统潮流计算数学模型112.2.2 电力系统节点分类 潮流计算的约束条件1.32.3 牛顿

9、-拉夫逊法求解潮流计算 牛顿-拉夫逊法原理 P-Q分解法潮流计算1.5第3章基于MATLA歆件P-Q法潮流计算203.1 P-Q分解法程序框图203.2 计算步骤及实现各部分功能的程序213.2.1 原始数据的输入213.2.2 导纳矩阵及B,B”形成233.2.3 计算不平衡功率出及修正相角A9253.2.4 计算不平衡功率i及修正相电压AVi263.2.5 程序运行结果的输出27第4章算例验证与分析284.1 算例说明及分析284.1.1 算例说明284.1.2 算例分析28-IV-南京工程学院4.2 根据算例输入相应节点线路参数284.3 算例运行结果29

10、结论32致谢33参考文献.34附录A36附录B46附录C63南京工程学院第1章绪论1.1 课题背景电力是衡量一个国家经济发展的主要指标，也是反映人民生活水平的重要标志，它已成为现代工农业生产、交通运输以及城乡生活等许多方面不可或缺的能源和动力。电力系统是由发电、输电、变电、配电和用电等环节组成的电能生产与消费系统。它的功能是将自然界的一次能源通过发电动力装置转化成电能，再经输电、变电和配电将电能供应到各用户。为实现这一功能，电力系统在各个环节和不同层次还具有相应的信息与控制系统，对电能的生产过程进行测量、调节、控制、保护、通信和调度，以保证用户获得安全、经济、优质的电能如图1-1所示。图1-1

11、电力系统示意图电力系统的出现，使电能得到广泛应用，推动了社会生产各个领域的变化，开创了电力时代，出现了近代史上的第二次技术革命。20世纪以来,电力系统的发展使动力资源得到更充分的开发，工业布局也更为合理，使电能的应用不仅深刻地影响着社会物质生产的各个侧面，也越来越广地渗透到人类日常生活的各个层面。电力系统的发展程度和技术水准已成为各国经济发展水平的标志之一。1.2 电力系统潮流计算1.2.1 电力系统潮流计算简介电力系统潮流计算1是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母南京工程学院线的电压，各元件中流过的功率，系统的功率损耗

12、等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性，可靠性和经济性。止匕外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以，潮流计算是研究电力系统的一种很重要和基础的计算。电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。1.2.2 电力系统潮流计算的基本要求对潮流计算的要求可以归纳为下面三点：计算方法的可靠性或收敛性；占用内存少，计算速度快；计算的方便性和灵活性。 收敛性一直以来，为了解决潮流计算中的收敛问题，人们

13、做了各种努力，首先要保证输入正确的数据，这样才可以使潮流收敛2。有不少研究是对基本的牛顿拉夫逊法和P-Q分解法的改进。在牛顿拉夫逊算法中，应用矩阵分块求逆方法对阶数较高的雅可比阵求逆计算进行改进，使阶数较高的雅可比矩阵的求逆变为阶数较低的四个子阵的求逆3,可以提高收敛速度。也有文章讨论通过对雅克比矩阵的简化4-5,提高计算速度。在P-Q分解法中，B和B“的不同组合形式会直接影响潮流的收敛性，文献6详细分析了对地并联支路导纳在形成B*时的作用，并且通过多个系统的潮流计算结果验证了：在形成B时考虑支路电阻，在形成B”时忽略支路电阻，并且将节点的并联电纳和线路的并联电纳做二倍处理时，P-Q分解法收敛

14、最快。针对病态系统，出现了最优乘子7和线性规划法网，有些文章还讨论了潮流计算中对负荷电压静特性的考虑，迭代中PV节点无功越限问题，以及不平衡功率的分配等问题网。 减少内存占用量自1967年美国学者WF.Tinney将稀疏矩阵技术引入电力系统潮流计算，用于牛顿潮流计算中，不仅减少了内存占用量，还大大提高了潮流计算的速度。虽然随着计算机技术的不断发展，计算机的内存不断增加，但是稀疏存储在节省运算量和计算机内存还是有其必要性的。文献10提到的三元素牛顿拉夫逊法是基于这样的假设，各节点只与2个相邻节点连接，且在取节点号时可假设节点1仅和节点2,3连接，节点2仅和节点3,4连接，以此类推，

15、其它节点之间的连接被忽略。这样进行高斯消去时，只需对偶数行南京工程学院消元，且只需消去一个元素，消元计算量将大大减少，但是这种算法节点编号麻烦，需要重新进行节点编号。动态形成十字链表，同时存储网络的拓扑结构信息和参数信息，适合于电力系统中运行方式变化及故障等情况下的对网络结构的修改11,与之相应的要增加存储空间，对于网络结构不经常变化的系统，没有必要采用这种存储方式。 计算的方便性和灵活性完善的潮流程序应该可以重复应用，可移植性高。为了使程序的可复用性最大化，以使软件的维护和升级所需费用和时间降低，文献12采用了CBD（基于构件的开发）和OOP（面向对象程序设计）技术。OOP技术

16、的应用主要集中在通过其继承特性重复应用已经存在的类。另外，还出现了一些对传统算法的综合和改进的算法。电力系统松弛算法13就是综合了时间增量松弛法与波形松弛法的优点，得到的一种改进的新方法。对初值要求严格，迭代速度快的特点，利用电力网的结构特点，使用高斯塞德尔迭代法的第一次迭代结果作为牛顿拉夫逊法的计算初值140这样既解决了牛顿-拉夫逊法对初值要求高的问题，又提高了收敛速度。随着GPS技术的出现，出现了一些应用PMU的潮流计算的方法。根据PMU测量精度和配置的不同，可以采用不同的方式将PMU的测量结果应用于潮流计算中15-16。但是由于PMU的配置还没有普遍，故该方法并不实用。1.3 潮流计算的

17、意义及其发展电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算的结果对电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术主要采用牛顿-拉夫逊法及其衍生的P-Q分解法。在运行方式管理中，潮流是确定电网运行方式的基本出发点

18、；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了多个在预想操作情况下电网的潮流分布以校验运行可靠性。在电力系统调度运行的多个领域都涉及到电网潮流计算。潮流是确定电力网络运行南京工程学院状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前提。在用数字计算机解电力系统潮流问题的开始阶段，普遍采取以节点导纳矩阵为基础的逐次代入法。这个方法的原理比较简单，要求的数字计算机内存量比较低，适应20世纪50年代电子计算机制造水平和当时电力系统理论水平。但它的收敛性较差，当系统规模变大时，迭代次数急剧上升，在计算中往往出现迭代不收敛的情况。这就迫使电力系统计算人员转向以

19、阻抗矩阵为基础的逐次代入法。阻抗法改善了系统潮流计算问题的收敛性，解决了导纳法无法求解的一些系统的潮流计算，在60年代获得了广泛的应用。阻抗法的主要缺点是占用计算机内存大，每次迭代的计算量大。当系统不断扩大时，这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺点，60年代中期发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需要存储各个地区系统的阻抗矩阵及它们之间联络线的阻抗，这样不仅大幅度地节省了内存容量，同时也提高了计算速度。克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法。这是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计

20、算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从60年代中期，在牛顿法中利用了最佳顺序消去法以后，牛顿法在收敛性。内存要求。速度方面都超过了阻抗法，成为60年代末期广泛采用的优秀方法。但随着网络规模的扩大以及计算机从离线计算向在线计算的发展，牛顿-拉夫逊法的内存需要量及计算速度方面越来越不适应要求。70年代中期出现的快速分解法比较成功的解决了上述问题，使潮流计算在牛顿-拉夫逊法的基础上向前迈进了一大步，成为取代牛顿-拉夫逊法的算法之一。近30年来，潮流问题算法的研究仍非常活跃，但是大多数研究是围绕着改进牛顿拉

21、夫逊法和P-Q分解法进行的。止匕外，随着人工智能理论的发展，遗传算法、人工神经网络、模糊算法也逐渐引入潮流计算。但是，到目前为止，这些新模型和算法还不能取代牛顿法和P-Q分解法的地位。由于电力系统的不断扩大和对计算速度要求的不断提高，计算机的并行技术也引起一些研究人员的兴趣。1.4本次毕业设计主要工作本文致力于研究分析电力网络的运行情况。结合电力系统潮流计算及暂态稳定计算的特点，设计MATLAB程序实现较复杂电力网络的潮流计算。探讨学习潮流计算的基本原理。P-Q法潮流计算的基础是牛顿-拉夫逊法潮流计算。所以，要实现对前者的了解，就要对后者有深刻的理解。本文较详细的说明了对各种电力系统网络导纳矩

22、阵形成过程，阐述了牛顿-拉夫逊法潮流计算公式的形成过程、原理。-4-南京工程学院基于MATLAB的编程。本文对MATLAB软件编程中所用的M语言有一定的介绍，说明了常用函数的作用，并对M语言中所涉及到的逻辑关系符，运算符，矩阵的正确输入以及简单的人机对话功能有了初步的解释。举例验证。真实的电力网络是复杂的又是简单的，复杂在于其网络的结构各式各样十分复杂，但大多结构是不同简单结构的不断重复；简单在于其所包含的器件基本上相同。对于本文来说，寻找到一个包含所有类型元件，并包含少许节点和线路的例子就可证明改程序对所有类型的电力系统网络适用。人机对话的形成。本文介绍了将程序的运算移到了幕后”，展现在用户

23、面前是个人性化界面的人机对话功能。本设计通过简单的程序段将MATLAB软件与Excel表格和TXT文档联系在一起，使用户更好的完成计算工作。南京工程学院第2章潮流计算的原理及具体算法过程2.1 电力网络的数学模型2.1.1 电力网络的基本方程在潮流问题中，任何复杂的电力系统都可以归纳为以下元件（参数）组成。发电机（注入电流或功率）负荷（注入负的电流或功率）输电线支路（电阻，电抗）变压器支路（电阻，电抗，变比）母线上的对地支路（阻抗和导纳）线路上的对地支路（一般为线路充电点容导纳）集中了以上各类型的元件的简单网络如图2-1所示。图2-1潮流计算用的电网结构图将图2-1中的发电机和负荷节点用无阻抗

24、线从网络中抽出（为不失一般性，将非发电机又非负荷的浮动节点当作零注入功率的母线抽出网络之外），剩下的部分即由接地和不接地支路组成一个无源线性网络如图2-2所示。南京工程学院4-国一磔无源线性网络|4/盘1图2-2潮流计算等值网络由图2-2,根据电路理论里边的戴维南定理，电力网络的节点电压方程它展开为：可表示为17:-IJJY1丫2丫3IIIYnlw|2Y21Y22Y23II1Y2nV2II3Y31Y32Y33II1Y3nHl1+FilV3(2-2)JnI1i，YYn2Yn3IIIYnn下Mn一I=yV(2-1)公式(2-1)中，Y是一个nxn阶节点导纳矩阵，其阶数就等于网络中节点总数。2.1.

25、2 导纳矩阵的形成节点导纳矩阵的对角元素Yi(i=1,2,n)成为自导纳。自导纳Yii数值上就等于在i节点施加单位电压，其他节点全部接地时，经节点i注入网络的电流，因此，它可以定义为：Yi=I/M(IVj0,j=)(2-3)节点i的自导纳Yi数值上就等于与节点直接连接的所有支路导纳的总和。节点导纳矩阵的非对角元素Yij(j=1,2,，n;i=1,2,，n;j错误！未找到引用源。i)称互导纳，由此可得互导纳Yj数值上就等于在节点i施加单位电压，其他节点全部接地时，经节点j注入网络的电流，因此可定义Yj=l|j/V(Vj0,j#i)(2-4)节点j,i之间的互导纳Yj数值上就等于连接节点j,i支路

26、到导纳的负值。显然，Yj恒等于Yji。互导纳的这些性质决定了节点导纳矩阵是一个对南京工程学院称稀疏矩阵。而且，由于每个节点所连接的支路数总有一个限度，随着网络中节点数的增加非零元素相对愈来愈少，节点导纳矩阵的稀疏度，即零元素数与总元素的比值就愈来愈高。注意字母几种不写法的不同意义：粗体黑字表示导纳矩阵，大写字母Yj为矩阵Y中的第i行第j列元素，即节点i和节点j之间的互导纳。小写字母i,j支路的导纳等于支路阻抗的倒数，yj=1/Zj。根据定义直接求取节点导纳矩阵时，注意以下几点：节点导纳矩阵是方阵，其阶数就等于网络中除去参考节点外的节点数。参考节点一般取大地，编号为零。节点导纳矩阵是稀疏矩阵，其

27、各行非零非对角元素就等于与该行相对应节点所连接的不接地支路数。节点导纳矩阵的对角元素就等于该节点所连接导纳的总和。因此，在没有接地支路的节点对应的行或列中，对角元素为非对角元素之和为负值。节点导纳矩阵的非对角元素等于连接节点i,j支路导纳的负值。因此，在一般情况下，节点导纳矩阵的对角元素往往大于非对角元素的负值。节点导纳矩阵一般是对称矩阵，这是网络的互易特性所决定的。从而，一般只要求取这个矩阵的上三角或下三角部分。2.1.3电力网络中几种特殊的数学模型输电线路的参数有四个：反映线路通过电流时产生有功功率损失效应的电阻；反映载流导线产生磁场效应的电感；反映线路带电时绝缘介质中产生泄漏电流及导线附

28、近空气游离而产生有功功率损失的电导；反映带电导线周围电场效应的电容。输电线路的这些参数通常可以认为是沿全长均匀分布的，每单位长度的参数为电阻P0O电感L0,电导g0及电容C0,其一项等值电路如图2-3所示。r0+j3Lg。jcoC0图2-3架空线等值电路在工程计算中，既要保证必要的精度，又要尽可能的简化计算，采用近似参数时，遵守以下规则:南京工程学院当线路长度l100km,线路可用”字型等值电路代替。r0+jlL0O图2-4架空线字型等值电路当线路长度100kml300km,线路可用错误！未找到引用源。”型等值电路代替。r0+jlL01久Ntj3c0/2j3c0/2图2-5架空线错误！未找到引

29、用源。”型等值电路当线路长度300km1,线路可用用型连接的多个错误！未找到引用源。”型电路来模拟，每个错误！未找到引用源。”型电路代替长度为200km-300km的一段线路。还需指出这里所讲的处理方法仅适应于工频下的稳态计算。本论文中，涉及到架空线导线问题的，均为前两种情况，所以第三种情况的示意图略。非标准变比变压器等值电路双绕组变压器如图2-6可用阻抗与一个理想变压器串联的电路表示理想变压器只是一个参数，那就是变比k=Vi/V2o本论文涉及到的变压器阻抗按实际变比分别归算到低压侧和高压侧，现将变压器型等值电路分别介绍如下：.111k:112.ViV2Zt图2-6双绕组变压器原

30、理图1、变压器阻抗归算到低压侧等值模型如图2-7所示。-9-南京工程学院图2-7阻抗归算到低压侧变压器原理图流入和流出理想变压器的功率相等，即：V1II1=MI/2k式中，k=Vi/V2是理想变压器的变比Vi和V2分别为变压器高、实际电压，从图2-7直接可得：V1k=VI12Tz从而可得:式中Yt=1/Zt,，V1V2ytV1YTV2Il二一2二2k2ZTkZTk2kVV2II2-YtV?kZTZtk又因节点电流方程应具有如下形式：；I1=YM+丫12V2、-|2=丫21V1+丫22V2将式(2-7)与(2-8)比较，得：-2Y1=Yt/k2丫12=-Yr/k，产-Yr/k工丫22=丫r因此可

31、得各支路导纳为：*2=-丫2=YT/kY21=-丫=YT/k1-k、0=丫11-丫2=tYtkk1丫20=丫22-，1=Yrlk由此可得用导纳表示的变压器型等值电路如图2-8所示(2-5)低绕组的(2-6)(2-7)(2-8)(2-9)-10-南京工程学院Ytk1-kk2YT图2-8阻抗归算到低压侧变压器等值电路2、变压器阻抗归算到高压侧等值模型如图2-9所示ViI1Ztk:1V2图2-9阻抗归算到高压侧变压器原理图同理,可得用导纳表示的变压器型等值电路如图2-10所示。kYT31-k)YTLl(k2-k)YT图2-10阻抗归算到高压侧变压器等值电路2.2 电力系统潮流计算2.2.1 电力系统

32、潮流计算数学模型电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。采用导纳矩阵时，节点注入电流和节点电压构成如式(2-1)所示线性方程组可展开如下形式：nIi八YjV(i=12Wn)(2-10)jm由于实际电网中测量的节点注入量一般不是电流而是功率，因此，必须-11-南

33、京工程学院将式中的注入电流用节点注入功率来表示。节点功率与节点电流之间的关系为：Si=Pi-jQi=Villi(2-11)式中Pi=巳PLDi,Qi-Qci_QLDi因此用导纳矩阵时，PQ节点可以表示为lL=Si/ViP-jQi*Vi把这个关系代入式（2-10）中，得到如下公式:P-jQi*Vin八YM(i=1,2M|n)j1(2-12)2.2.2 电力系统节点分类用一般的电路理论求解网络方程，目的是给出电压源（或电流源）研究网络内的电流（或电压）分布，作为基础的方程式，一般用线性代数方程式表示。然而在电力系统中，给出发电机或负荷连接母线上电压或电流（都是向量）的情况是很少的，一般是给出发电机

34、母线上发电机的有功功率（P）和母线电压的幅值（V）,给出负荷母线上负荷消耗的有功功率（P）和无功功率（Q）。主要目的是由这些已知量去求电力系统内的各种电气量。所以，根据电力系统中各节点性质的不同，很自然地把节点分成三类18:PQ节点对这一类点，事先给定的是节点功率（P,Q）,待求的未知量是节点电压向量（V,9）,所以叫PQ节点。通常变电所母线都是PQ节点，当某些发电机的输出功率P,Q给定时，也作为PQ节点。PQ节点上的发电机称之为PQ机（或PQ给定型发电机）。在潮流计算中，系统大部分节点属于PQ节点。PV节点这类节点给出的参数是该节点的有功功率P及电压幅值V,待求量为该节点的无功功率Q及电压向

35、量的相角日。这类节点在运行中往往要有一定可调节的无功电源，用以维持给定的电压值。通常选择有一定无功功率储备的发电机母线或者变电所有无功补偿设备的母线做PV节点处理，PV节点上的发电机称为PV机（或PV给定型发电机）。平衡节点在潮流计算中，这类节点一般只设一个。对该节点，给定其电压值，并在计算中取该节点电压向量的方向作为参考轴，相当于给定该点电压向量的角度为零。也就是说，对平衡节点给定的运行参数是V和8,因此有城为VB节点，而待求量是该节点的巳Q,整个系统的功率平衡由这一节点承-12-南京工程学院担。关于平衡节点的选择，一般选择系统中担任调频调压的某一发电厂(:发电机)，有时也可能按其他原则选择

36、，例如，为提高计算的收敛性。可以选择出线数多或者靠近电网中心的发电厂母线作平衡节点。以上三类节点4个运行参数P,Q,V,8中，已知量都是两个，待求量也是两个，只是类型不同而已。2.2.3 潮流计算的约束条件电力系统运行必须满足一定的技术和经济上的要求。这些要求构成了潮流问题中某些变量的约束条件，常用的约束条件如下：节点电压应满足小于节点最大额定电压并大于最小额定电压，即：Vimin-V00=1,2,，n)(2-16)在待求量x的某一个初始估计值x(0)附近，将上式展开成泰勒级数并略去二阶及以上的高阶项，得到如下的经线性化的方程组：错误！未找到引用源。fx(0)f(x(0)；：x(0)=0xxx

37、(2-17)上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量：错误！未找到引用源。x(0)=f(x(0)fx(2-18)将错误！未找到引用源。x(0)和错误！未找到引用源。相加，得到变量的第一次改进值x错误！未找到引用源。接着就从错误！未找到引用源。出发，重复上述计算过程。因此从一定的初值x(0)错误！未找到引用源。出发，应用牛顿法求解的迭代格式为:(,(k)(k)(k)(2-19)I”3)Ax=-f(x)(k粕(k)(k)x=x+Ax上两式中：错误！未找到引用源。f(x)是函数错误！未找到引用源。对于变量x的一阶偏导数矩阵，即雅可比矩阵J,k为迭代次数。由上式可见，牛顿法的核心便

38、是反复形式并求解修正方程式。应用牛顿拉夫逊法时，当初始估计值x(0)错误！未找到引用源。和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。牛顿拉夫逊潮流算法突出的优点是收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代45次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。牛顿拉夫逊法也具有良好的收敛可靠性，对于以节点导纳矩阵为基础的高斯法呈病态的系统，牛顿法拉夫逊也能可靠收敛。牛顿法所需的内存量及每次迭代所需时间均较高-14-南京工程学院斯法多。牛顿-拉夫逊法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有可能根本不收敛或收敛到

39、一个无法运行的节点上。对于正常运行的系统，各节点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值(也称为平直电压)，如假定：V：0)=i,0(0)=0错误！未找到引用源。这样一般能得到满意的结果。但若系统因无功紧张或其它原因导致电压质量很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问题。解决这个问题的办法可以用高斯法迭代12次，以此迭代结果作为牛顿-拉夫逊法的初值。也可以先用直流法潮流求解一次以求得一个较好的角度初值，然后转入牛顿-拉夫逊法迭代。2.3.2P-Q分解法潮流计算 极坐标下的潮流计算模型以下讨论的

40、是用极坐标形式叙述牛顿-拉夫逊法潮流的求解过程。当采用极坐标时，潮流问题的待求量为各节点电压的幅值和相角Vi,9错误!未找到引用源。由于平衡节点的电压向量是给定的，不参加迭代；系统有m个PQ节点，则PV节点有n-(m+1)个。对于PV节点，因为电压幅值给定，这就减少了n-(m+1)个未知数，同时，PV节点的注入无功功率为可调节量，不能预先给定，AQ也就失去了约束作用。因此，对于PV节点，将式(2-20)按泰勒级数展开并略去AQ,AV二次及以后各项，修正方程为式(2-22)。取错误！未找到引用源。VVz0,y产G+jB得到潮流方程的极坐标形式：Pi=ViZVjGijcos0ij+Bijsin0i

41、j)(2-20)1j土、nQi”EVjGijsinbjBijcosbj)j士式中9ij=0i-9j错误！未找到引用源。表示i,j两节点电压的相角差。则每个节点的功率差表示为以下方程组：Pi=Pis-Pi(2-21)”=Qisg方程式把节点功率差AP(8,V),AQ(8,V)错误！未找到引用源。表示为节点电压V和相角错误！未找到引用源。的函数，其中APis错误！未找到引用源。，AQis错误！未找到引用源。为外加功率，表示发电机等注入功率时符号为正，表示负荷等流出功率时符号为负，节点功率差为两者的代数和。则由式(2-20)、(2-21)可知，对各节点的错误！未找到引用源。P,错误！未找到引用源。Q

42、分别对错误！未找到引用源。求偏导数并将各节点的-15-南京工程学院错误！未找到引用源。值代入，可得到雅可比矩阵J。其具体形式如下式:PlcPPiaPiaPiaPi&Pic8C020n_LRi仅2机工给P2给P2&P2P2&P2辔P2c0c02石QnXN必Ml&PnA&PnLPnA&Pn-龄PnA出PnLc02CQnX叫J=咫咫gg注Qica君92c&LRi匆2CVn&Q2Q2&Q2&Q2松Q2Q2CaC02C9nXRi12Vln_L小Qn,Qn双工RnL凶n.caC02CRi旬nL则可得修正方程组:一ARIHiiHi2HinNliNi2PP2H.2iH22HuN?iN22*HnjHn.J2Hn

43、JniNnnNnJ2+1.F.QJiJi2JimLiL2J2iJ”J2mL2ibIJmiJm2JmmLmLm2(2-22)该修正方程可用分块矩阵的形式简化如下:(2-23)NA91LAV/V_在式(2-23)中，电压幅值的修正量采用AVi/Vi的形式，是为了使雅克比矩阵中各元素有比较相似的形式。雅克比矩阵中各9,V元素对式(2-24)取偏导数求得，其计算式为:P1=-ViVjGijsin-ij-Bucos-ji=jHuD,n.TP-ZVjGijSin3BijCOfi=j-jj1jTLZ-16-南京工程学院F-Vj-WjGijsin-BijCosi=jV=i二jL，j11QnFVi=VVjGsi

44、niij-BijCosiij)2V：BiiJVijijIrPVj=JvVjfeijcoseij+BijsineJi卢jz-1。Nij-二Dn2dVi=VjGijco汨uBijsinaj-2ViGiiJVij与上|西rVjGcosG+Bijsin6)ijI二Jiji=jQj-=Vi二VjGjcosftBijsinfti-0jm士需要指出的是，在上述的迭代过程中，由于某个PV节点因无功越线而转为PQ节点时，修正方程中也应更换或增加相应的行。在极坐标系下，应增加一行对应于节点的无功功率不平衡量错误！未找到引用源。Qi的关系式，在列向量中，若Q（v）2Q时，取Q（v）=Q；若Q（v）EQ时，取iima

45、sxiimaxiimin(v)Q二Q。iiminPg解法潮流计算P-Q分解法是从简化牛顿法极坐标形式计算潮流程序的基础上提出来的。由于交流高压电网中输电线路等元件的x错误！未找到引用源。r,因此电力系统呈现了这样的物理特性，即电力系统中有功功率主要与各节点电压相角有关，无功功率则主要受各节点电压幅值的影响。因此，P-Q分解法的基本思想是：把节点功率表示成电压向量的极坐标方程式，以有功功率误差作为修正电压相角的依据，以无功功率误差作为修正电压幅值的依据,(2-24)把有功功率和无功功率迭代分开来进行。即将式（2-23）简化为:一AP_H0一A9IJaqjLLJ（AV/V_这样，由于

46、把2n阶的线性方程组变成了两个n阶的线性方程组，因而，计算量和内存方面都有很大的改善。但是，H、L在迭代过程中仍然不断变化，而且又都是不对称矩阵，因此，牛顿法的第二个简化，就是把式（2-24）中的系统矩阵简化为在迭代过程中不变的对称矩阵。一般情况下，线路两端电压相角差是不大的（通常不超过10-20。），因此可以认为：-17-南京工程学院(2-25)Gijsin6冬B，j此外，与系统各节点无功功率相应的导纳BLI必定远远小于该节点自导纳的虚部，即BbR：bV，(2-26)因此,Q处沦考虑到以上的关系，式（2-24）的系数矩阵中的各元素可表示为:(2-27)Hj=VVjBj(-1,2,n)LijW

47、iVjBijQJ-Zim）而系数矩阵H和L则可以分别写成：-ViBiVV2B2V1*V1B12V2VB1nNnLV2B22V2V2B2nNn工irVnlBnJlVlVnrBnW2VnBnnlVniV1VVnjB21BnqB12B1n.LB22B2n工Bn12Bn_1V2Vn。(2-28)%1BVd1V1B11V1L=V2BUV1VmBm1VV1B12V2V2B22V2VmBm2V2VBmVmV2B2mVmV1=V2.L，D2BVD2B11B12B21B22Bm1Bm2BmmVmV2Vm_(2-29)将式(2-28)和(2-29)代入式(2-24)中，得到l-Apl=-VdJbVd1hd1-aq

48、LVd2bhv1-18-南京工程学院用v/口Vd2别左乘以上两式，便得D1VD1由揖I-BWJa1VD21*1=-b这就是简化了的修正方程式,一卬1:它们也可展开写成ViP2V2Pr-BnDB12BgV16B21口一B229B2nJ.V2aBjjBn也Bn山Vni0n4(2-30)Viy2BllB21B12B22BlmlVll(2-31)V2Qm.Vm_|miBm2V.2Bmm_V在这两个修正方程式中系数矩阵元素就是系数导纳矩阵的虚部，因而系数矩阵是对称矩阵，且在迭代过程中保持不变。这就大大减少了计算工作量。式(2-30)和(2-31)构成了P-Q分解法迭代过程的基本方程。-19-南京工程学院

49、第3章基于MATLA欧件P-Q法潮流计算目前，电子计算机已广泛应用于电力系统的分析计算尤其是在潮流计算之中，MATLAB是其基本应用软件之一。3.1 P-Q分解法程序框图P-Q分解法的实现过程应用到循环嵌套原理，其程序框图如图3-1所小。图3-1P-Q分解法程序框图-20-南京工程学院3.2 计算步骤及实现各部分功能的程序3.2.1 原始数据的输入 原始数据的输入为了让程序更有通用性，以及方便初学用户使用，近几年，MATLAB编程大多都用到GUI功能。本次设计也不例外，用到了简单的人机对话界面作为原始数据输入界面。在这次设计中，用到了一个Excel表格作为原始数据输入界面，通过该

50、界面，用户不用依照矩阵的形式，将一连串的数据输入，而是，按照图表的提示，在图表中填入要求的电力系统节点、支路的参数，节点、支路个数以及要求精度即可。数据输入界面如图3-2所示。图3-2Excel输入界面在此，对数据的输入有以下几点说明：在节点信息里，节点电压及相角为迭代计算时所设的初值。在节点信息里，节点类型一栏中，“1蓑示PQ节点，“2裳示PV节点，“3表示平衡节点。KT一栏要求输入的是变压器的变比，非标准变比变压器，KT=k（k错误！未找到引用源。1）,标准变压器KT=1,若该线路无变压器，KT=0o输入变压器电阻、电抗时，如无特殊说明，均采用归算到低压侧的数值，再进行计算。折算到哪一侧的

51、标志一栏，表示如果系统简图的支路起始节点处于高压侧则请输入“1；否则输入“0：本条程序默认节点数，支路数均在在100以下，可以解决绝大部分电力系统的潮流问题，若遇到超大系统，可对程序做稍加调整，仍然使用。-21-南京工程学院原始数据的输入程序段1、输入参数：point为节点信息矩阵，zhilu为支路信息矩阵x=xlsread(pqinput.xls,A2:A2);y=xlsread(pqinput.xls,B2:B2);e=xlsread(pqinput.xls,B4:B4);point=xlsread(pqinput.xls,D3:H100);zhilu=xlsread(pqinput.xls,J3:R100);TYPE=zeros(x,1);U=zeros(x,1);a=zeros(x,1);P=zeros(x,1);Q=zeros(x,1);I=zeros(y,1);J=zeros(y,1);Rij=zeros(y,1);Xij=zeros(y,1);Zij=Rij+j*Xij;Y=zeros(x);G=zeros(x);B=zeros(x);B0=zeros(y,1);RT=zeros(y,1);XT=zeros(y,1);ZT=RT+j*XT;KT=zeros(y,1);2、矩阵赋初值TYPE=p