毕业论文-基于matlab的电力系统潮流计算

摘要潮流计算是电力系统的一项重要分析功能，是进行故障计算，继电保护整定，安全分析的必要工具。是对复杂电力系统正常和故障条件下稳态运行状态的计算。电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础。潮流计算是电力系统的各种计算的基础，同时它又是研究电力系统的一项重要分析功能，在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用电力系统潮流计算来定量的比较供电方案和现有电力系统运行方式的合理性，可靠性和经济性。实际电力系统的潮流技术那主要采用牛顿拉夫逊法。传统的潮流计算程序缺乏图形用户界面，结果显示不直观，难于与其他分析功能集成。网络原始数据输入工作量大且易于出错。随着计算机技术的飞速发展，MICROSOFTWINDOWS操作系统早已被大家所熟悉，其友好的图形用户界面已成为PC机的标准，而DOS操作系统下的应用程序因其界面不够友好，开发具有WINDOWS风格界面的电力系统分析软件已成为当前的主流趋势。另外，传统的程序设计方法是结构化程序设计方法，该方法基于功能分解，把整个软件工程看作是一个个对象的组合，由于对某个特定问题域来说，该对象组成基本不变，因此，这种基于对象分解方法设计的软件结构上比较稳定，易于维护和扩充。本文介绍了图形化潮流计算软件的开发设计思想和总体结构，阐述了该软件所具备的功能和特点。结合电力系统的特点，软件采用MATLAB语言运行于WINDOWS操作系统的图形化潮流计算软件。本系统的主要特点是操作简单，图形界面直观，运行稳定。计算准确。计算中，算法做了一些改进，提高了计算速度，各个类的有效封装又使程序具有很好的模块性。可维护性和可重用性。关键词电力系统潮流计算；牛顿拉夫逊法潮流计算；MATLABABSTRACTPOWERFLOWCALCULATIONISANIMPORTANTFUNCTIONOFTHEELECTRICPOWERSYSTEMISANALYZED,FAULTCALCULATION,RELAYPROTECTION,SAFETYANALYSISTOOLSNECESSARYISTHECALCULATIONOFTHESTEADYSTATEOPERATIONOFNORMALANDFAULTCONDITIONSOFCOMPLEXPOWERSYSTEMUNDERTHETHERESULTSOFPOWERFLOWCALCULATIONISTHEBASISOFCOMPUTATIONANDFAULTANALYSISOFPOWERSYSTEMSTABILITYPOWERFLOWCALCULATIONISTHEBASISOFALLKINDSOFPOWERSYSTEMCALCULATION,ANDITISALSOANIMPORTANTFUNCTIONOFPOWERSYSTEMANALYSIS,THERESEARCHONTHEOPERATINGMODEOFPOWERSYSTEMPLANNINGANDDESIGNOFTHEEXISTINGPOWERSYSTEM,NEEDTOUSEPOWERFLOWCALCULATIONTOQUANTITATIVECOMPARISONOFTHERATIONALITYOFPOWERSUPPLYSCHEMEANDOPERATIONMODEOFTHEEXISTINGPOWERSYSTEM,RELIABILITYANDECONOMYTHEACTUALPOWERFLOWTECHNOLOGYTHATMAINLYUSESTHENEWTONRAPHSONMETHODTHETRADITIONALFLOWCALCULATIONPROGRAMLACKSAGRAPHICALUSERINTERFACE,THEDISPLAYISNOTINTUITIVE,ITISDIFFICULTTOINTEGRATEWITHOTHERANALYSISFUNCTIONTHENETWORKINPUTDATAANDHEAVYWORKLOADANDERRORPRONEWITHTHERAPIDDEVELOPMENTOFCOMPUTERTECHNOLOGY,MICROSOFTWINDOWSOPERATINGSYSTEMHASLONGBEENFAMILIAR,ITSFRIENDLYGRAPHICALUSERINTERFACEHASBECOMEPCSTANDARD,ANDTHEAPPLICATIONOFTHEDOSOPERATINGSYSTEMBECAUSEOFITSINTERFACEISNOTFRIENDLYENOUGH,THEPOWERSYSTEMANALYSISSOFTWAREDEVELOPMENTWITHWINDOWSSTYLEINTERFACEHASBECOMETHEMAINTRENDINADDITION,THETRADITIONALDESIGNMETHODISASTRUCTUREDPROGRAMDESIGNMETHOD,THISMETHODISBASEDONTHEFUNCTIONDECOMPOSITION,THESOFTWAREENGINEERINGASACOMBINATIONOFOBJECTS,DUETOAPARTICULARPROBLEMDOMAIN,THECOMPOSITIONOFTHEOBJECTISESSENTIALLYTHESAME,THEREFORE,BASEDONTHESOFTWARESTRUCTUREOFTHEOBJECTDECOMPOSITIONMETHODONTHEDESIGNOFSTABLE,EASYTOMAINTAINANDEXPANSIONTHEDEVELOPMENTOFDESIGNIDEAANDOVERALLSTRUCTUREOFTHISPAPERINTRODUCESTHEGRAPHICALPOWERFLOWCALCULATIONSOFTWARE,EXPOUNDSTHEFUNCTIONANDFEATURESOFTHESOFTWAREACCORDINGTOTHECHARACTERISTICSOFTHEPOWERSYSTEM,THEGRAPHICALTRENDMATLABLANGUAGEOPERATINGSYSTEMRUNNINGONAWINDOWSCALCULATIONSOFTWARETHEMAINFEATUREOFTHISSYSTEMISSIMPLE,INTUITIVEGRAPHICALINTERFACE,ACCURATECALCULATIONOFSTABLEOPERATIONINTHECALCULATION,THEALGORITHMMADESOMEIMPROVEMENTS,IMPROVETHECALCULATIONSPEED,ANDTHEEFFECTIVEENCAPSULATIONOFCLASSPROGRAMMODULEHASGOODMAINTAINABILITYANDREUSABILITYKEYWORDSPOWERFLOWCALCULATIONNEWTONRAPHSONPOWERFLOWCALCULATIONMATLAB目录引言1第1章电力系统潮流计算概述211电力系统叙述212潮流计算简介213潮流计算的意义及其发展3第2章潮流计算的数学模型421导纳矩阵的原理及计算方法422潮流计算的基本方程723电力系统节点分类924潮流计算的约束条件10第3章牛顿拉夫逊法概述1231牛顿拉夫逊法基本原理1232牛顿拉夫逊法潮流求解过程1333牛顿拉夫逊法的程序框图17第4章MATLAB概述1841MATLAB简介1842MATLAB应用在潮流计算中的优势1843矩阵的运算19第5章潮流计算主界面设计实现2151主界面介绍2152数据初始化2153潮流计算2254数据处理2355数据传递的问题2356例某电网接线图及给定的参数2457运算结果24第6章系统潮流计算的前沿算法及发展前景3161保留非线性算法3162最优潮流分析法3163OPF分析法32结论与展望33致谢34参考文献35附录36附录A基于MATLAB的牛顿拉夫逊法潮流计算程序清单36附录B外文文献及译文39附录C参考文献的题录摘要45插图清单图21双绕组变压器原理图6图22变压器等值电路7图23潮流计算用的电网结构图8图24潮流计算等值网络8图31牛顿拉夫逊法的程序框图17图51某电网接线图24引言潮流计算是在给定电力系统网络结构、参数和决定系统运行状态的边界条件的情况下确定系统稳态运行状态的一种基本方法，是电力系统规划和运营中不可缺少的一个重要组成部分。可以说，它是电力系统分析中最基本、最重要的计算，是系统安全、经济分析和实时控制与调度的基础。是电力系统研究人员长期研究的一个课题。MATLAB自1980年问世以来，它的强大的矩阵处理功能给电力系统的分析、计算带来许多方便。在处理潮流计算时，其计算机软件的速度已无法满足大电网模拟和实时控制的仿真要求，而高效的潮流问题相关软件的研究已成为大规模电力系统仿真计算的关键。随着计算机技术的不断发展和成熟，对MATLAB潮流计算的研究为快速、详细地解决大电网的计算问题开辟了新思路。电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。在用数字计算机解电力系统潮流问题的开始阶段，普遍采取以节点导纳矩阵为基础的逐次代入法。这个方法的原理比较简单，要求的数字计算机内存量比较下，适应50年代电子计算机制造水平和当时电力系统理论水平。但它的收敛性较差，当系统规模变大时，迭代次数急剧上升，在计算中往往出现迭代不收敛的情况。这就迫使电力系统计算人员转向以阻抗矩阵为基础的逐次代入法。阻抗法改善了系统潮流计算问题的收敛性，解决了导纳法无法求解的一些系统的潮流计算，在60年代获得了广泛的应用。阻抗法的主要缺点是占用计算机内存大，每次迭代的计算量大。当系统不断扩大时，这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺点，60年代中期发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需要存储各个地区系统的阻抗矩阵及它们之间联络线的阻抗，这样不仅大幅度地节省了内存容量，同时也提高了计算速度。克服阻抗法缺点的另一途径是采用牛顿拉夫逊法。这是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从60年代中期，在牛顿法中利用了最佳顺序消去法以后，牛顿法在收敛性。内存要求。速度方面都超过了阻抗法，成为60年代末期以后广泛采用的优秀方法。第1章电力系统潮流计算概述11电力系统叙述电力工业发展初期，电能是直接在用户附近的发电站或称发电厂中生产的，各发电站孤立运行。随着工农业生产和城市的发展，电能的需要量迅速增加，而热能资源如煤田和水能资源丰富的地区又往往远离用电比较集中的城市和工矿区，为了解决这个矛盾，就需要在动力资源丰富的地区建立大型发电站，然后将电能远距离输送给电力用户。同时，为了提高供电可靠性以及资源利用的综合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来的整体，即称为电力系统。电力系统加上发电机的原动机如汽轮机、水轮机，原动机的力能部分如热力锅炉、水库、原子能电站的反应堆、供热和用热设备，则称为动力系统。现代电力系统提出了“灵活交流输电与新型直流输电”的概念。灵活交流输电技术是指运用固态电子器件与现代自动控制技术对交流电网的电压、相位角、阻抗、功率以及电路的通断进行实时闭环控制，从而提高高压输电线路的输送能力和电力系统的稳定水平。新型直流输电技术是指应用现电力电子技术的最新成果，改善和简化变流站的造价等。运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了电网在预想操作情况下电网的潮流分布以校验运行可靠性。在电力系统调度运行的多个领域都涉及到电网潮流计算。潮流是确定电力网络运行状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前提。12潮流计算简介电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和基础的计算。电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点1计算方法的可靠性或收敛性。2对计算机内存量的要求。3计算速度。4计算的方便性和灵活性。电力系统潮流计算问题在数学上是一组多元非线性方程式求解问题，其解法都离不开迭代。因此，对潮流计算方法，首先要求它能可靠地收敛，并给出正确答案。由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大，潮流计算的方程式阶数也越来越高，对这样的方程式并不是任何数学方法都能保证给出正确答案的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。13潮流计算的意义及其发展电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿拉夫逊法。在运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了多个在预想操作情况下电网的潮流分布以校验运行可靠性。在电力系统调度运行的多个领域都涉及到电网潮流计算。潮流是确定电力网络运行状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前提。在用数字计算机解电力系统潮流问题的开始阶段，普遍采取以节点导纳矩阵为基础的逐次代入法。这个方法的原理比较简单，要求的数字计算机内存量比较下，适应50年代电子计算机制造水平和当时电力系统理论水平。但它的收敛性较差，当系统规模变大时，迭代次数急剧上升，在计算中往往出现迭代不收敛的情况。这就迫使电力系统计算人员转向以阻抗矩阵为基础的逐次代入法。阻抗法改善了系统潮流计算问题的收敛性，解决了导纳法无法求解的一些系统的潮流计算，在60年代获得了广泛的应用。阻抗法的主要缺点是占用计算机内存大，每次迭代的计算量大。当系统不断扩大时，这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺点，60年代中期发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需要存储各个地区系统的阻抗矩阵及它们之间联络线的阻抗，这样不仅大幅度地节省了内存容量，同时也提高了计算速度。克服阻抗法缺点的另一途径是采用牛顿拉夫逊法。这是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从60年代中期，在牛顿法中利用了最佳顺序消去法以后，牛顿法在收敛性。内存要求。速度方面都超过了阻抗法，成为60年代末期以后广泛采用的优秀方法。第2章潮流计算的数学模型21导纳矩阵的原理及计算方法211自导纳和互导纳的确定方法电力网络的节点电压方程21BUYI为节点注入电流列向量，注入电流有正有负，注入网络的电流为正，流出网络BI的电流为负。根据这一规定，电源节点的注入电流为正，负荷节点为负。既无电源又无负荷的联络节点为零，带有地方负荷的电源节点为二者代数之和。为节点电压列向量，由于节点电压是对称于参考节点而言的，因而需先选定参BU考节点。在电力系统中一般以地为参考节点。如整个网络无接地支路，则需要选定某一节点为参考。设网络中节点数为（不含参考节点）,则，均为NN列向量。BIU为NN阶节点导纳矩阵。BY节电导纳矩阵的节点电压方程BIY展开为22121311122223123333123NNNNNNNYIUYYI是一个NN阶节点导纳矩阵，其阶数就等于网络中除参考节点外的节点数。节BY点导纳矩阵的对角元素I1，2，N成为自导纳。自导纳数值上就等于在I节IIY点施加单位电压，其他节点全部接地时，经节点I注入网络的电流，因此，它可以定义为23/0,IIJYIUI节点I的自导纳数值上就等于与节点直接连接的所有支路导纳的总和。I节点导纳矩阵的非对角元素J1，2，NI1，2，。，NJI称互导纳，JI由此可得互导纳数值上就等于在节点I施加单位电压，其他节点全部接地时，经节IJY点J注入网络的电流，因此可定义为24/0,JIJIJIUI节点J，I之间的互导纳数值上就等于连接节点J，I支路到导纳的负值。显然，IJY恒等于。互导纳的这些性质决定了节点导纳矩阵是一个对称稀疏矩阵。而且，由IJYJI于每个节点所连接的支路数总有一个限度，随着网络中节点数的增加非零元素相对愈来愈少，节点导纳矩阵的稀疏度，即零元素数与总元素的比值就愈来愈高。212节点导纳矩阵的性质及意义节点导纳矩阵的性质1为对称矩阵，。如网络中含有源元件，如移相变压器，则对称性不再BYIJYI成立。2对无接地支路的节点，其所在行列的元素之和均为零，即B,110NNIJJIJIY对于有接地支路的节点，其所在行列的元素之和等于该点接地支路的导纳。利用这一性质，可以检验所形成节点导纳矩阵的正确性。3具有强对角性对角元素的值不小于同一行或同一列中任一元素。BY4为稀疏矩阵，因节点I，J之间无支路直接相连时0，这种情况在实际电IJY力系统中非常普遍。矩阵的稀疏性用稀疏度表示，其定义为矩阵中的零元素与全部元素之比，即，式中为中的零元素。随节点数的增加而增加2ZNSZBYSN50，可达92；100，可达90；500，可达99，充分利用节点导纳矩NSN阵的稀疏性可节省计算机内存，加快计算速度，这种技巧称为稀疏技术。节点导纳矩阵的意义是NN阶方阵，其对角元素I1，2，N称为自导纳，非对角元素BYIYI，J1，2，N，称为互导纳。将节点电压方程展开为IJIJBUYI11212212NNNNIYYUI25/0,IIIJYUIJIJ表明，自导纳在数值上等于仅在节点I施加单位电压而其余节点电压均为零（即其I余节点全部接地）时，经节点I注入网络的电流。其显然等于与节点I直接相连的所有支路的导纳之和。同时可见。表明，互导纳/0,1,2,IJIJIYIJNJ在数值上等于仅在节点J施加单位电压而其余节点电压均为零时，经节点I注入网络的电流，其显然等于即。为支路的导纳，负号表示该电流流出网络。IJYIJIJYIJ如节点IJ之间无支路直接相连，则该电流为0，从而0。IJY注意字母几种不写法的不同意义粗体黑字表示导纳矩阵，大写字母代矩阵中的IJYB第I行第J列元素，即节点I和节点J之间的互导纳。小写字母I，J支路的导纳等于支路阻抗的倒数数。根据定义直接求取节点导纳矩阵时，注意以下几点1节点导纳矩阵是方阵，其阶数就等于网络中除去参考节点外的节点数。参考节点一般取大地，编号为零。2节点导纳矩阵是稀疏矩阵，其各行非零非对角元素就等于与该行相对应节点所连接的不接地支路数。3节点导纳矩阵的对角元素就等于各该节点所连接导纳的总和。因此，与没有接地支路的节点对应的行或列中，对角元素为非对角元素之和的负值。4节点导纳矩阵的非对角元素等于连接节点I，J支路导纳的负值。因此，一般情况下，节点导纳矩阵的对角元素往往大于非对角元素的负值。5节点导纳矩阵一般是对称矩阵，这是网络的互易特性所决定的。从而，一般只要求求取这个矩阵的上三角或下三角部分。213非标准变比变压器等值电路变压器型等值电路更便于计算机反复计算,更适宜于复杂网络的潮流计算双绕组变压器可用阻抗与一个理想变压器串联的电路表示理想变压器只是一个参数,那就是变比。现在变压器阻抗按实际变比归算到低压侧为例,推导出变压器型等值21UK电路图21双绕组变压器原理图流入和流出理想变压器的功率相等KIU/2126I1式中,是理想变压器的变比,和分别为变压器高,低绕组的实际电压21UK2从图B直接可得（212TUKIZ7）从而可得12T122TYUIZ（212T122TUI8）式中,又因节点电流方程应具有如下形式TY1/Z112IYU（229）将式（18）与（19）比较，得21TY/1TY/22因此可得各支路导纳为（212T10122TY/10）由此可得用导纳表示的变压器型等值电路图22变压器型等值电路22潮流计算的基本方程在潮流问题中，任何复杂的电力系统都可以归纳为以下元件（参数）组成（1）发电机（注入电流或功率）（2）负荷（注入负的电流或功率）（3）输电线支路（电阻，电抗）（4）变压器支路（电阻，电抗，变比）（5）母线上的对地支路（阻抗和导纳）（6）线路上的对地支路（一般为线路充电点容导纳）集中了以上各类型的元件的简单网络如图图23潮流计算用的电网结构图图24潮流计算等值网络采用导纳矩阵时，节点注入电流和节点电压构成以下线性方程组,其中YUI12NII12UN可展开如下形式212NIIJJ1IYU,2由于实际电网中测量的节点注入量一般不是电流而是功率，因此必须将式中的注入电流用节点注入功率来表示。节点功率与节点电流之间的关系为213ISIIPJQI式中，因此用导纳矩阵时，节点可以表示为IGILDIPIGILDIQP把这个关系代入式中得IS/IIIIJIU（21,2NIIIJPJYUN14）式（34）就是电力系统潮流计算的数学模型潮流方程。它具有如下特点1它是一组代数方程，因而表征的是电力系统的稳定运行特性。2它是一组非线性方程，因而只能用迭代方法求其数值解。3由于方程中的电压和导纳既可以表为直角坐标，又可表为极坐标，因而潮流方程有多种表达形式极坐标形式，直角坐标形式和混合坐标形式。取，得到潮流方程的极坐标形式IIU|IJIJIJYY2151NIIIIIJIJPJQUY取，得到潮流方程的直角坐标形式IIUEJFIJIJIJYGB2161111NNIIIJIJIIJIJJJNNIIIJIJIIJIJJJPEFGFBEQFBFEF取，得到潮流方程的混合坐标形式IIUIJIJIJYG2171COSINISNIIJIIJJJJIIJJIJIJJPUBQ不同坐标形式的潮流方程适用于不同的迭代解法。例如利用牛顿拉夫逊迭代法求解，以直角坐标和混合坐标形式的潮流方程为方便；而PQ解耦法是在混合坐标形式的基础上发展而成，故当然采用混合坐标形式。4它是一组N个复数方程，因而实数方程数为2N个但方程中共含4N个变量P，Q，U和，I1，2，N，故必须先指定2N个变量才能求解。23电力系统节点分类用一般的电路理论求解网络方程，目的是给出电压源或电流源研究网络内的电流或电压分布，作为基础的方程式，一般用线性代数方程式表示。然而在电力系统中，给出发电机或负荷连接母线上电压或电流都是向量的情况是很少的，一般是给出发电机母线上发电机的有功功率P和母线电压的幅值U，给出负荷母线上负荷消耗的有功功率P和无功功率Q。主要目的是由这些已知量去求电力系统内的各种电气量。所以，根据电力系统中各节点性质的不同，很自然地把节点分成三类1PQ节点对这一类点，事先给定的是节点功率P，Q，待求的未知量是节点电压向量U，所以叫PQ节点。通常变电所母线都是PQ节点，当某些发电机的输出功率P。Q给定时，也作为PQ节点。PQ节点上的发电机称之为PQ机或PQ给定型发电机。在潮流计算中，系统大部分节点属于PQ节点。2PU节点这类节点给出的参数是该节点的有功功率P及电压幅值U，待求量为该节点的无功功率Q及电压向量的相角。这类节点在运行中往往要有一定可调节的无功电源。用以维持给定的电压值。通常选择有一定无功功率储备的发电机母线或者变电所有无功补偿设备的母线做PU节点处理。PU节点上的发电机称为PU机或PU给定型发电机。3平衡节点在潮流计算中，这类节点一般只设一个。对该节点，给定其电压值，并在计算中取该节点电压向量的方向作为参考轴，相当于给定该点电压向量的角度为零。也就是说，对平衡节点给定的运行参数是U和，因此有城为U节点，而待求量是该节点的P。Q，整个系统的功率平衡由这一节点承担。关于平衡节点的选择，一般选择系统中担任调频调压的某一发电厂或发电机，有时也可能按其他原则选择，例如，为提高计算的收敛性。可以选择出线数多或者靠近电网中心的发电厂母线作平衡节点。以上三类节点4个运行参数P。Q。U。中，已知量都是两个，待求量也是两个，只是类型不同而已。24潮流计算的约束条件电力系统运行必须满足一定技术和经济上的要求。这些要求够成了潮流问题中某些变量的约束条件，常用的约束条件如下1节点电压应满足218MINMAX1,2IIUN2从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在额定电压附近。PU节点电压幅值必须按上述条件给定。因此，这一约束条件对PQ节点而言3节点的有功功率和无功功率应满足219MINMAXGIGIPQPQ节点的有功功率和无功功率，以及PU节点的有功功率，在给定是就必须满足上述条件，因此，对平衡节点的P和Q以及PU节点的Q应按上述条件进行检验。4节点之间电压的相位差应满足230MAX|IJIJIJ为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。这一约束的主要意义就在于此。因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束条件。常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用约束条件进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至修改系统的运行方式，重新进行计算。第3章牛顿拉夫逊法概述31牛顿拉夫逊法基本原理电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿拉夫逊法。牛顿拉夫逊法简称牛顿法在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程。即通常所称的逐次线性化过程。对于非线性代数方程组即310FX12,0INFX1,2IN在待求量X的某一个初始估计值附近，将上式展开成泰勒级数并略去二阶及以上的高阶项，得到如下的经线性化的方程组3200FXFX上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量330010FF将和相加，得到变量的第一次改进值。接着就从出发，重复上述计0X1X1X算过程。因此从一定的初值出发，应用牛顿法求解的迭代格式为0X34KKFFX351KX上两式中是函数对于变量X的一阶偏导数矩阵，即雅可比矩阵JK为FF迭代次数有上式可见，牛顿法的核心便是反复形式并求解修正方程式。牛顿法当初始估计值和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。0X牛顿潮流算法突出的优点是收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代45次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。牛顿法也具有良好的收敛可靠性，对于对以节点导纳矩阵为基础的高斯法呈病态的系统，牛顿法也能可靠收敛。牛顿法所需的内存量及每次迭代所需时间均较高斯法多。牛顿法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有可能根本不收敛或收敛到一个无法运行的节点上。对于正常运行的系统，各节点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值也称为平直电压，如假定或3601IU0I01IE0IF1,2IQNIS这样一般能得到满意的结果。但若系统因无功紧张或其它原因导致电压质量很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问题。解决这个问题的办法可以用高斯法迭代12次，以此迭代结果作为牛顿法的初值。也可以先用直流法潮流求解一次以求得一个较好的角度初值，然后转入牛顿法迭代。32牛顿拉夫逊法潮流求解过程以下讨论的是用直角坐标形式的牛顿拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量由于平12,NFFE衡节点的电压向量是给定的，因此待求两共需要2N1个方程式。事实上，除21N了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。对PQ节点来说，和是给定的，因而可以写出ISPISQ（370IIJIJIIJIJJJIJIJIJIJISIJJPFFEGEBF）对PV节点来说，给定量是和，因此可以列出ISPIV（38）2200ISIJIJIIJIJIJJIIFFEEB求解过程大致可以分为以下步骤1形成节点导纳矩阵。2将各节点电压设初值U。3将节点初值代入相关求式，求出修正方程式的常数项向量。4将节点电压初值代入求式，求出雅可比矩阵元素。5求解修正方程，求修正向量。6求取节点电压的新值。7检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步。8计算支路功率分布，PV节点无功功率和平衡节点柱入功率。以直角坐标系形式表示迭代推算式采用直角坐标时,节点电压相量及复数导纳可表示为39IIJJJVEFYGB将以上二关系式代入上式中,展开并分开实部和虚部假定系统中的第1,2,M号为PQ节点,第M1,M2,N1为PV节点,根据节点性质的不同,得到如下迭代推算式对于PQ节点31011NNIIIJIJIIJIJJJIIIJIJIIJIJJJEBFFBEQFGFEGF,2M对于PV节点3111122NNIIIJIJIIJIJJJIIIPEBFFBEV,IMN对于平衡节点平衡节点只设一个,电压为已知,不参见迭代,其电压为312NNVEJF修正方程式235和236两组迭代式工包括2N1个方程选定电压初值及变量修正量号之后代入式235和236,并将其按泰勒级数展开,略去二次方程及以各项,得到修IEF正方程如下313WJU11212MNPQWPU111MNEFFUEF314111111111111111111MMNNMMMMNPPPPEFEFEFEFQQQQEFEFEFEJ11111111112222211NMMMMNNFPPPPEFEFEFEFUU222111111111111222222NNNNNNNNMMUUEFEFEFEF22NN雅可比矩阵各元素的算式式328中,雅可比矩阵中的各元素可通过对式324和325进行偏导而求得。当时，雅可比矩阵中非对角元素为JI315220IIIJIJJJIIIJIJJJJJPQGEBFEFUEF当时,雅可比矩阵中对角元素为JI316112NIIJIJIIJIIJIJIIJJNIIJIJIIJIIJIJIIJJIIJIIPGEBFEFEFQFEFEEGBFUEFF由式329和3210看出,雅可比矩阵的特点1阵中各元素是节点电压的函数,在迭代过程中,这些元素随着节点电压的变化而变化。2导纳矩阵中的某些非对角元素为零时,雅可比矩阵中对应的元素也是为零。若,则必有。0IJY0IJJ3可比矩阵不是对称矩阵雅可比矩阵各元素的表示如下1,2IQNISIJIJIIJIJIJIIJJGEBFJIPHFEIJIJIIJIJIJIIJJEFJINFBFFIJIJIIJIJIJIIJJEGFJIQMFFEIJIJIIJIJIJIIJJEBFJILFF20IIJIJUREJ2IIJIJSFJF33牛顿拉夫逊法的程序框图YN图31牛顿拉夫逊法程序框图启动输入原始数据形成导纳矩阵给定电压初值0,EF置对于PU节点，计算。QP,对于PQ节点，计算2是否P,计算雅可比矩阵各元素求修正方程式，求FE,用，E1求修正FFF节点电压以FFE1,1以1E按系统的潮流分布计算平衡节点的功率及线路功率输出第4章MATLAB概述41MATLAB简介目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求1计算速度快。2内存需要少。3计算结果有良好的可靠性和可信性。4适应性好，亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强。5简单。MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。MATLAB程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB与C语言和FORTRAN语言相比更容易被掌握。通过M语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。另外，MATLAB提供了一种特殊的工具工具箱（TOOLBOXES）这些工具箱主要包括信号处理（SIGNALPROCESSING）、控制系统（CONTROLSYSTEMS）、神经网络（NEURALNETWORKS）、模糊逻辑FUZZYLOGIC、小波WAVELETS和模拟（SIMULATION）等等。不同领域、不同层次的用户通过相应工具的学习和应用，可以方便地进行计算、分析及设计工作。MATLAB设计中，原始数据的填写格式是很关键的一个环节，它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计，主要应从使用的角度出发，原则是简单明了，便于修改。42MATLAB应用在潮流计算中的优势MATLAB由一系列工具组成。这些工具方便用户使用MATLAB的函数和文件，其中许多工具采用的是图形用户界面。包括MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于用户浏览帮助、工作空间、文件的浏览器。随着MATLAB的商业化以及软件本身的不断升级，MATLAB的用户界面也越来越精致，更加接近WINDOWS的标准界面，人机交互性更强，操作更简单。而且新版本的MATLAB提供了完整的联机查询、帮助系统，极大的方便了用户的使用。简单的编程环境提供了比较完备的调试系统，程序不必经过编译就可以直接运行，而且能够及时地报告出现的错误及进行出错原因分析。MATLAB是一个高级的矩阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。用户可以在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用程序（M文件）后再一起运行。新版本的MATLAB语言是基于最为流行的C语言基础上的，因此语法特征与C语言极为相似，而且更加简单，更加符合科技人员对数学表达式的书写格式。使之更利于非计算机专业的科技人员使用。MATLAB是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。函数中所使用的算法都是科研和工程计算中的最新研究成果，而前经过了各种优化和容错处理。在通常情况下，可以用它来代替底层编程语言，如C和C。在计算要求相同的情况下，使用MATLAB的编程工作量会大大减少。基于MATLAB的电力系统潮流计算使计算机在计算、分析、研究复杂的电力系统潮流分布问题上又前进了一步。不管采用什么算法,所有的潮流计算都是基于矩阵的迭代运算。而MATLAB语言正是以处理矩阵见长,实践证明，MATLAB语言在电力系统潮流计算仿真研究中的应用是可行的,而且由于其强大的矩阵处理功能,完全可以应用于电力系统的其它分析计算中；用MATLAB语言编程效率高,程序调试十分方便,可大大缩减软件开发周期，如果像控制界一样开发出电力系统自己的专用工具箱,将系统分析用的一些基本计算以函数的形式直接调用,那么更高层次的系统软件也可以很容易地实现。43矩阵的运算矩阵是MATLAB数据存储的基本单元，而矩阵的运算是MATLAB语言的核心，在MATLAB语言系统中几乎一切运算均是以对矩阵的操作为基础的。矩阵的基本数学运算包括矩阵的四则运算、与常数的运算、逆运算、行列式运算、秩运算、特征值运算等基本函数运算，这里进行简单介绍。1四则运算矩阵的加、减、乘运算符分别为“，”，用法与数字运算几乎相同，但计算时要满足其数学要求在MATLAB中矩阵的除法有两种形式左除“”和右除“/”。在传统的MATLAB算法中，右除是先计算矩阵的逆再相乘，而左除则不需要计算逆矩阵直接进行除运算。通常右除要快一点，但左除可避免被除矩阵的奇异性所带来的麻烦。在MATLAB6中两者的区别不太大。2与常数的运算常数与矩阵的运算即是同该矩阵的每一元素进行运算。但需注意进行数除时，常数通常只能做除数。3基本函数运算矩阵的函数运算是矩阵运算中最实用的部分，常用的主要有以下几个DETA求矩阵A的行列式EIGA求矩阵A的特征值INVA或A1求矩阵A的逆矩阵RANKA求矩阵A的秩TRACEA求矩阵A的迹（对角线元素之和）我们在进行工程计算时常常遇到矩阵对应元素之间的运算。这种运算不同于前面讲的数学运算，为有所区别，我们称之为数组运算。4基本数学运算数组的加、减与矩阵的加、减运算完全相同。而乘除法运算有相当大的区别，数组的乘除法是指两同维数组对应元素之间的乘除法，它们的运算符为“”和“/”或“”。前面讲过常数与矩阵的除法运算中常数只能做除数。在数组运算中有了“对应关系”的规定，数组与常数之间的除法运算没有任何限制。另外，矩阵的数组运算中还有幂运算（运算符为）、指数运算（EXP）、对数运算（LOG）、和开方运算（SQRT）等。有了“对应元素”的规定，数组的运算实质上就是针对数组内部的每个元素进行的。矩阵的幂运算与数组的幂运算有很大的区别。5逻辑关系运算逻辑运算是MATLAB中数组运算所特有的一种运算形式，也是几乎所有的高级语言普遍适用的一种运算。第5章潮流计算主界面设计实现51主界面介绍主界面是进行潮流计算的的操作界面，能完成数据的初始化、进行潮流计算、并显示潮流计算的相关结果以供用户分析。具体来讲就是点击第一个弹出框选择潮流计算数据组，并在主界面上显示所选择数据组的电路接线图，单击计算按钮就能完成潮流计算并将潮流计算结果保存下来，单击第二个弹出框可以查看具体的结果，同时还将画出潮流计算电压迭代曲线。52数据初始化本设计将数据保存在EXCEL表格里，在主界面上单击弹出框选择数据组从而从EXCEL表格里调入数据完成数据初始化。本设计定义了7个数据供主程序调用，他们分别是节点数N，支路数N1，平衡节点标号ISB，精度值PR，节点参数矩阵B1，支路参数矩阵B2和节点号及其对地阻抗形成的举证X。从EXCEL表格中调入数据采用XLSREAD函数，这样EXCEL表格中的数据就以矩阵的形式存放在MATLAB中。但XLSREAD函数不能读入虚数，其将虚数标记符I当做字符处理，所以本设计把数据的实部和虚部分开分别读入，然后在MATLAB中将其转化回去。数据初始化程序如下所示VALGETHOBJECT,VALUESTRGETHOBJECT,STRINGSWITCHSTRVALCASE电力系统1M1CASE电力系统2M2ENDIFM1AXLSREADEFIRXLSHANDLESNA1HANDLESA2HANDLESISBA3HANDLESPRA4BXLSREADEBRXLSCXLSREADEBIXLSHANDLESB1B,16JC,16HANDLESB2B,712JC,712DXLSREADEXXLSHANDLESXD,12JD,34ELSEIFM2AXLSREADEFIR2XLSHANDLESNA1HANDLESN1A2HANDLESISBA3HANDLESPRA4BXLSREADEBR2XLSCXLSREADEBI2XLSHANDLESB1B,16JC,16HANDLESB2B,712JC,712DXLSREADEX2XLSHANDLESXD,12JD,34END53潮流计算本设计通过计算按钮调用潮流计算主程序完成潮流计算。潮流计算程序根据牛顿拉夫逊法编制，对主程序介绍如下1导纳矩阵的形成求导纳矩阵Y中的非对角元元素，若无变压器,则直接为线路阻抗分之一取负值，若有变压器，为线路阻抗乘以KT后分之一再取负值。求导纳矩阵Y中对角元元素，无变压器时为加上线路对地电导的一半乘J，有变压器时，对角元元素就与所输入的折算到哪一侧有关，如果支路起始端处于高压侧，支路起始节点的自导纳中要加上