电力系统稳定分析与计算课程设计.docx

《电力系统稳定分析与计算课程设计.docx》由会员分享，可在线阅读，更多相关《电力系统稳定分析与计算课程设计.docx（46页珍藏版）》请在三一办公上搜索。

1、课程设计报告电力系统稳定分析与计算毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作 者 签 名： 日 期： 指导教师签名： 日期： 使用授权说明本人完全了解 大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权

2、保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。作者签名： 日 期： 学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。作者签名： 日期： 年 月 日学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意

3、学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权 大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。涉密论文按学校规定处理。作者签名：日期： 年 月 日导师签名： 日期： 年 月 日注 意 事 项1.设计（论文）的内容包括：1）封面（按教务处制定的标准封面格式制作）2）原创性声明3）中文摘要（300字左右）、关键词4）外文摘要、关键词 5）目次页（附件不统一编入）6）论文主体部分：引言（或绪论）、正文、结论7）参考文献8）致谢9）附录（对论文支持必要时）2.论文字数要求：理工类设计（论文）正

4、文字数不少于1万字（不包括图纸、程序清单等），文科类论文正文字数不少于1.2万字。3.附件包括：任务书、开题报告、外文译文、译文原文（复印件）。4.文字、图表要求：1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准请他人代写2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符合国家技术标准规范。图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒手画3）毕业论文须用A4单面打印，论文50页以上的双面打印4）图表应绘制于无格子的页面上5）软件工程类课题应有程序清单，并提供电子文档5.装订顺序1）设计（论文）2）附件：按照任务书、开题报告、外文

5、译文、译文原文（复印件）次序装订3）其它目 录课程设计题目C：电力系统稳定分析和计算31.电网参数计算31.1 线路参数的计算和标幺化31.2 节点信息统计32.电网潮流计算32.1 采用Matlab计算电网潮流32.1.1 网络信息处理32.1.2 基于Matlab的PQ分解法程序32.1.3 PQ分解法潮流计算的结果分析32.2采用Power World计算潮流32.2.1 Power world仿真软件简介32.2.2 基于Power World软件构建潮流模型32.2.3 采用power world仿真软件的潮流计算结果32.3 Matlab仿真计算结果与Power World仿真软件

6、计算结果对比34. 转移阻抗的求解34.1 负荷、发电机、变压器、线路的处理34.1.1 负荷的处理34.1.2 发电机的处理34.1.3 变压器的处理34.1.4 输电线路的处理34.2 网络变换（稳态时）34.2.1 消去节点5、634.2.2 消去节点2、3、4、734.3 各种运行方式下转移阻抗的求解34.4 功角特性计算35. 故障仿真35.1 采用Matlab软件分析故障35.1.1 三相短路故障35.1.2 单相接地故障35.2采用Power World仿真软件分析故障（二阶）35.2.1 线路2-4中点三相短路时的故障分析35.2.2 线路2-4中点三相短路时的故障分析35.2

7、.3 线路2-4中点单相接地故障时的仿真分析35.2.4 Matlab与Power World仿真结果的比较35.3 采用Power World仿真软件分析故障（三阶）35.3.1 三相短路时的故障分析35.3.2 线路2-4中点发生单相接地时的故障分析36. 发电机模型选择对稳定计算结果的影响3参考文献3附录：3附1. PQ分解法的Matlab程序3附2. Matlab潮流计算结果截图3附3. Power World潮流计算截图3附4. 零序网（左）负序网（右）3附5. 单相接地网架等效图3课程设计题目C：电力系统稳定分析和计算姓名： 指导教师： 李晓华 一、 一个220kV分网结构和参数如

8、下：18km12km12km12km25km16km15km#1#2#4#5#6#3G10km500kV站（#1）的220kV母线视为无穷大母线，电压恒定在230kV。图中，各变电站负荷均等效接于220kV母线。变电站参数和负荷如下表：编号类型220kV最大负荷，MVA#1500kV站平衡节点#2220kV站180+j110#3220kV站350+j180#4220kV站260+j120#5220kV站250+j120#6220kV站150+j70各线路长度如图所示。所有线路型号均为LGJ-2*300，基本电气参数为：正序参数： r = 0.054/km, x = 0.308/km, C =

9、0.0116 F/km；零序参数： r0 = 0.204/km, x0 = 0.968/km, C0 = 0.0078 F/km；40C长期运行允许的最大电流：1190A。燃煤发电厂G有三台机组，发电机与升压变之间采用单元接线。电厂220kV侧采用单母分段接线，正常运行时分段开关闭合。发电机组主要参数如下表：机组台数单台容 量（MW）额定电压（EV）功率因数升压变容量MVA升压变Vs%XdXdXqTd0TJ=2H330010.50.8535010.52.00.251.7411当发电机采用三阶模型时，励磁环节（含励磁机和励磁调节器）模型如下（不考虑PSS）:上图中参数如下：TR=0，KA=25，

10、 TA=0.1，Te=0.15，KE=1， KF=0.05，TF=0.7 发电厂升压变均采用Y/Y0接线，变比10.5kV/242kV。不计内阻和空载损耗。发电厂按PV方式运行，高压母线电压定值为1.05VN。发电厂厂用电均按出力的7%考虑。稳定仿真中不考虑发电厂的调速器和原动机模型。负荷采用恒阻抗模型，负序阻抗标幺制取0.2。二、 设计的主要内容：1、手工进行参数计算和标幺化，形成潮流计算的网架参数；2、用Matlab编制潮流计算程序，要求采用P-Q分解潮流计算方法。3. 考虑该电厂开机三台，均为额定容量输出的方式，用编制的程序计算潮流。3、用PowerWorld软件进行潮流计算并与自己编制

11、的软件计算结果进行校核和分析；4、设#2和#4母线之间双回线路中一回的中点分别发生以下2种故障：1） 1s时发生三相短路，1.1s同时切除故障线路三相；2） 1s时单相接地短路，1.1s时同时切除故障线路三相。1.9s时三相重合闸。因重合于永久性故障，2.1s时再次切除故障线路。试手工计算序网以及用网络变换法求解转移阻抗；5、针对问题4，用Matlab编制稳定计算程序（三台机可并联等值成一台机），发电机采用二阶经典模型（注：用ode45函数既可求解），用软件计算出摇摆曲线，要求输出发电机功角，角速度。 6、用PowerWorld软件的分析问题4，并与编程计算结果进行比较校核。7、发电机采用三阶

12、模型，用PowerWorld作为分析工具，对问题4的两种故障方案进行稳定计算，给出摇摆曲线，并计算故障的极限切除时间。8、比较两种模型的仿真结果，分析发电机模型选择对于稳定计算结果的影响。9、编制课程设计报告。三、 设计要求和设计成果：1、2位同学为一组，自行分工，但任务不能重复；2、每位同学对自己的设计任务编写课程设计说明书一份；3、一组同学共同完成一份完整的设计报告；2、设计说明和报告应包含： 以上设计任务每一部分的计算过程和结果分析； 所编制的潮流和稳定计算源程序（主要语句应加注释）； 潮流计算结果（潮流图） 稳定计算的功角曲线等； 网络变换法求解转移阻抗的变换过程图。附注： ODE函数

13、说明Matlab提供了一阶常微分方程组求解的系列函数：ode*。包括：ode45, ode23, ode113等，还有针对刚性系统的ode15s，ode23s等。这里可采用ode45编程（大家也可选择和对比其它函数，不同编号采用的数值积分算法不同）。函数形式：t, y=ode45(odefun,t1,tf,x0,options);说明：odefun：列向量1*n，通过函数计算柄输出的微分方程的右端项；t1,tf: 分别制定积分的时间起点和终点；x0：列向量1*n，状态变量初值options：微分优化参数，是一个结构体，使用odeset可以设置其具体参数，详细内容查看帮助。t：为时间列向量1*m

14、y为状态变量计算结果矩阵，m行代表时间点，n列代表n个状态变量的时间序列值。例如，求解如下微分方程：初值为x1=1,x2=0从0s积分到3s，步长：0.1s则，首先定义函数myfunc，计算微分方程右端项的值：function dx=myfunc(x) dx=x(2) 2*sin(x(1) ;Ode45函数引用如下：x0= 1, 0options=odeset;options.reltol=1e-8; t,y=ode45(myfunc,0,3,x0,options)华南理工大学电力系 电气工程和及其自动化 专业课程设计(论文)任务书兹发给2010级电气工程及其自动化1班学生 课程设计任务书,内

15、容如下：课程设计题目： 电力系统运行方式分析和计算 应完成的项目：用Matlab编制PQ分解法潮流计算程序，完成典型运行方式的潮流计算并进行分析；用PowerWord软件对自己编制的软件计算结果进行校核和分析；用Matlab编制稳定计算程序，发电机采用二阶经典模型，要求给出网络变换法求解转移阻抗的变换过程图；选择2-3种故障方案，计算故障的极限切除时间和极限切除角；用Powerworld作为分析工具，发电机采用三阶模型，对上面的2-3种故障方案进行稳定计算，计算故障的极限切除时间，分析发电机模型选择对于稳定计算结果的影响，并且分析励磁调节系统参数变化对于稳定计算结果的影响。参考资料以及说明电力

16、系统分析（上、下册）华中科技大学出版发电厂电气部分高等学校教材电网调度运用技术东北大学出版社PowerWorld 15使用手册基于MATLAB/Simulink的系统仿真技术与应用本毕业设计（论文）任务书于2013年12月23日发出，应于2014年 1 月 3 日前完成，然后提交课程考试委员会进行答辩。系主任 批准 年 月 日教员组主任 审核 年 月 日指导老师 签 发 2013年 12月 23日1.电网参数计算1.1 线路参数的计算和标幺化220kV分网结构的型等值电路图1如下图所示：图1.电网等效线路图待求解电网的基准功率取为100MVA，基准电压220kV，则有 （1）线路的正序参数：

17、（2）线路阻抗标幺值计算公式： （3）线路等值电纳标幺值计算公式： （4）当线路为双回线路时，线路阻抗为单回线线路的一半，导纳为单回线路的两倍。各线路参数的标幺化计算如下：线路L12长度15km，则：Z12=0.515（0.054+j0.308）=（0.405+j2.32） B12=22500.01161510-6=1.093310-4S 标幺值：Z12*=（0.405+j2.32）100/2202=8.367810-4+j0.004793 B12*=1.093310-42202/100=0.052916线路L13长度18km，则：Z13=0.518（0.054+j0.308）=（0.486+

18、j2.772） B13=22500.01161810-6=1.311310-4S 标幺值：Z13*=（0.486+j2.77）100/2202=1.00410-3+j0.005723 B12*=1.311310-42202/100=0.063467线路L24长度12km，则：Z24=0.512（0.054+j0.308）=（0.324+j1.848） B24=22500.01161210-6=8.741810-5S 标幺值：Z24*=（0.324+j1.848）100/2202=6.6910-4+j0.003818 B24*=8.741810-52202/100=0.04231线路L45长度1

19、0km，则：Z45=0.510（0.054+j0.308）=（0.27+j1.54） B45=22500.01161010-6=7.284810-5S 标幺值：Z45*=（0.27+j1.54）100/2202=5.5810-4+j0.003182 B45*=7.284810-52202/100=0.03526线路L46长度12km，则：Z46=12（0.054+j0.308）=（0.648+j3.696） B46=2500.01161210-6=4.373110-5S 标幺值：Z46*=（0.648+j3.696）100/2202=1.33910-3+j0.007636 B46*=4.373

20、110-52202/100=0.021166线路L3G长度16km，则：Z3G=0.516（0.054+j0.308）=（0.432+j2.464） B3G=22500.01161610-6=1.1655710-4S 标幺值：Z3G*=（0.432+j2.464）100/2202=8.9310-4+j0.005091 B3G*=1.1655710-42202/100=0.05641线路L5G长度25km，则：Z5G=0.525（0.054+j0.308）=（0.675+j3.85） B5G=22500.01162510-6=1.821210-4S 标幺值：Z5G*=（0.675+j3.85）1

21、00/2202=1.39510-3+j0.007955 B5G*=1.821210-42202/100=0.08815线路L6G长度12km，则：Z6G=0.512（0.054+j0.308）=（0.324+j1.848） B6G=22500.01161210-6=8.741810-5S 标幺值：Z6G*=（0.324+j1.848）100/2202=6.6910-4+j0.003818 B6G*=8.741810-52202/100=0.04231综上，电网的线路参数的有名值和标幺值如下表所示：表1.线路参数的标幺值和有名值统计线路名长度/km有名值 标幺值阻抗Z/导纳B/S阻抗Z*/导纳B

22、*/SL122150.405+j2.321.093310-48.36810-4+j0.004790.05292L132180.486+j2.7721.311310-41.00410-3+j0.0057230.06347L242120.324+j1.8488.741810-56.6910-4+j0.0038180.04231L452100.27+j1.547.284810-55.5810-4+j0.0031820.03526L46120.648+j3.6964.373110-51.33910-3+j0.0076360.02117L3G2160.432+j2.4641.165610-48.9310

23、-4+j0.0050910.05641L5G2250.675+j3.851.821210-41.39510-3+j0.0079550.08815L6G2120.324+j1.8488.741810-56.6910-4+j0.0038180.04231线路零序参数的标幺值计算与正序相似，如表2.所示。表2.线路零序参数的计算线路名长度/km 标幺值阻抗Z*/导纳B*/SL122153.16110-3+j0.0150.03558L132183.79310-3+j0.0180.04270L242122.52910-3+j0.0120.02846L452102.10710-3+j0.0100.0237

24、2L46125.05810-3+j0.0240.01423L3G2163.37210-3+j0.0160.03795L5G2255.26910-3+j0.0250.05930L6G2122.52910-3+j0.0120.02846故障L24125.05810-3+j0.0240.014231.2 节点信息统计在给定的220kV网架中共有7个节点，其中节点1是500kV变电站，为平衡节点，视为无穷大系统，电压稳定在230kV，即1.045455VN。节点2、3、4、5、6均为PQ节点，并且带有一定的负荷。此外，节点6为PV节点，由于发电机并不总是满载运行，在正常运行时，考虑到7%的厂用电，因而

25、发电机机组发出的总有功功率为：满载运行时， （6）于是，当发电机机组满载时，每台机组的出力为837/3=279MW，发电机机端额定电压为10.5kV，出线侧的高压母线电压稳定在1.05VN，各个节点的参数信息如表3所示。表3.各节点参数一览表编号123456G节点类型平衡节点PQ节点PQ节点PQ节点PQ节点PQ节点PV节点有功功率P/MW180350260250150837无功功率Q/MVar10018012012070电压V/kV2302422.电网潮流计算2.1 采用Matlab计算电网潮流通过对待求解电网的等效线路参数的计算与标幺化，形成节点导纳矩阵。然后利用各个节点的类型信息，利用PQ

26、分解法进行电网潮流计算。本文编写了基于Matlab平台的PQ分解法计算程序，对电网进行了潮流计算，得出了系统稳定时各节点的电压，以及线路损耗等，分析了系统稳定运行时各线路的输送功率。2.1.1 网络信息处理形成电网线路的节点导纳矩阵是进行潮流计算前必须要做的准备。首先利用线路的标幺值参数整合成矩阵，使其包含线路的所有信息，再形成节点导纳矩阵。具体如下：构造支路信息矩阵B1，每一行代表一条线路，每行中的各元素表示线路不同的信息（注：节点G记为节点7）：1、支路首端节点号；2、支路末端节点号；3、支路阻抗标幺值；4、支路对地电纳标幺值；5、支路变比；6、支路存在变压器时，若首段：末端=k:1，则为

27、1，反之则为0，如果支路不存在变压器则默认为0；7、是否为变压器支路。本题形成的B1矩阵见附录。构造节点信息矩阵B2，题目给定的网络共有7个节点，本文编程中的B2矩阵为76的矩阵，每一行代表一个节点的信息，每一行中的各元素的含义：1、该节点的发电机功率（注入）；2、该节点的负荷功率（输出）；3、空行（保持矩阵格式）；4、PV节点电压给定值和PQ节点的初始值；5、各节点电压相角；6、节点类型（其中，1表示平衡节点，2表示PQ节点，3表示PV节点）。待求解网架中，1号节点为平衡节点，26号节点为PQ节点，7号节点（即母线G）为PV节点。本题形成的B2矩阵见附录。2.1.2 基于Matlab的PQ分

28、解法程序编程思想如下：1）根据支路信息矩阵B1和节点信息矩阵B2，先形成待求解电网的节点导纳矩阵：2）先对空矩阵加入对角线元素，即先加入每个节点的对地导纳y0；3）根据支路信息矩阵B1，对导纳矩阵再加入互导纳Ypq；4）将互导纳Ypq加入到导纳矩阵的对角线元素，形成自导纳。至此，形成节点导纳矩阵Y。5）初始化各节点电压（幅值和相角）的存储矩阵和各节点注入功率（实部和虚部）的存储矩阵，以存放迭代过程中的中间值和输出最终值；6）构造各节点的潮流方程，进入潮流计算程序，计算不平衡功率P，判断P是否小于可接受误差。若P大于可接受误差，则进行相角修正。计算不平衡功率Q判断Q是否小于可接受误差。若Q大于可

29、接受误差，则进行电压幅值修正。7）当和都满足要求时，计算结束。计算平衡节点功率和全部线路功率。8）基于MATLAB平台，得出计算结果，采用误差精度为=1.010-5。由结果可知，共用6次迭代便能满足精度要求，共用时间仅为0.073687s；2.1.3 PQ分解法潮流计算的结果分析利用2.1.2节的PQ分解法的matlab程序计算的电网潮流结果如附录中附2图所示。将matlab程序计算的结果整理如下：1）各节点的电压和注入功率如表4.所示。表4. 各节点电压和注入功率节点编号1234567基准电压220220220220220220220标幺电压1.04551.03751.04171.03581

30、.03631.04291.05实际电压（kV）230.01228.25229.17227.88227.99229.44231.00节点相角（度）0-0.6033-0.2606-0.7560-0.7119-0.19320.3608节点注入功率（标幺值）3.5671+j1.7501-1.8000+j1.0000-3.5000-j1.8000-2.6000-j1.2000-2.5000-j1.2000-1.5000-j0.70008.3700+j3.9248节点注入功率（实际值）356.71+j175.01-180.00-j100.00-350.00-j180.00-260.00-j120.00-2

31、50.00-j120.00-150.00-j70.00837.00+j392.48其中，matlab中各节点的相角是以弧度的形式表示的，其换算关系为：角度=弧度180，表3.中的功角是角度的形式表示的。2）各支路的状态如表5.所示。表中的功率值均为标幺值，基准功率SB=100MW。表.5 各支路的状态支路名称线路首端有功功率（MW）线路首端无功功率（Mvar）线路末端有功功率（MW）线路末端无功功率（Mvar）有功损耗（MW）无功损耗（Mvar）L122.61171.2732-2.6052-1.29330.0065-0.0201L130.95530.4769-0.9543-0.53990.00

32、11-0.0630L240.80520.2933-0.8047-0.33610.0005-0.0428L45-0.2829-0.15240.28300.11480.0001-0.0376L46-1.5123-0.71151.51580.70840.0035-0.0031L37-2.5457-1.26012.55231.23590.0066-0.0242L57-2.7830-1.31482.79511.28820.0121-0.0267L67-3.0158-1.40843.02261.40080.0068-0.0077全网的功率损耗可按下式（5）计算： （5） 即：网损率2.2采用Power W

33、orld计算潮流2.2.1 Power world仿真软件简介Power World Simulator（仿真器）是一个电力系统仿真软件包，其设计界面友好，并有高度的交互性。该仿真软件能够进行专业的工程分析。而且由于其可交互性和可绘图性，它也可以用于向非专业用户解释电力系统的运行操作。该仿真器是一个集成的产品，其核心是一个全面、强大的潮流计算程序。它能够有效地计算高达10,0000个节点的电力网络，因此当它作为一个独立的潮流分析软件包时，性非常实用。与其它商业潮流计算软件包不同，该软件可以让用户通过生动详细的全景图来观察电力系统。此外，系统模型可以通过使用仿真软件的图形编辑工具很容易地进行修改

34、，用户只需轻轻点击几下鼠标就可以在检修期间切换线路、增加新的线路或发电机、确定新的交易容量。仿真器广泛地使用了图形和动画功能，大大地增强了用户对系统特性、问题和约束的理解，以便于用户对系统进行维护。它基本的工具包括经济调度、区域功率经济分配分析、功率传输分配因子计算（PTDF）、短路分析以及事故分析等功能的工具。2.2.2 基于Power World软件构建潮流模型在编辑模式下，画出网架的模型，并且将节点1设置为平衡节点，电压恒定为VN=1.0455VB，节点1连接一台无穷大的发电机。具体参数设置如图2.所示。图2. 平衡机的参数设置节点6为PV节点，P=837MW，V=1.05VN，每台发电

35、机的输出功率P=297MW。具体参数设置如图3.所示。其余各节点只设置基准电压VB=220kV即可。图3. 发电机的参数设置采用power world画出的单项接线图如图4.所示。图4. 单项接线图2.2.3 采用power world仿真软件的潮流计算结果运行模式下，点击tools中的slove选项即可对电网进行潮流计算，在case information中点击Power Flow List可以查看各个节点的电压和线路的潮流。Power World仿真软件中的潮流计算结果如附录三所示。整理后的潮流计算结果如下。1）各节点的电压如表6所示。表6、各节点的电压节点编号基准电压kV标幺电压实际电压

36、(kV)相角 (度)12201.0455230.01022201.0375228.25-0.6032201.0418229.20-0.2642201.0358227.88-0.7652201.0364228.00-0.7162201.0430229.46-0.1972201.0500231.000.362）发电机的状态如表7.所示。节点编号发电机编号有功出力（MW）无功出力（MW）设定电压1#1356.70175.741.04557#1279.00130.571.05007#2279.00130.571.05007#3279.00130.571.05003）各支路的状态如表8.所示。表8. 各

37、支路的状态线路名称首端有功功率（MW）首端无功功率（Mvar）末端有功功率（MVA）末端无功功率（Mvar）有功损耗（MW）无功损耗（Mvar）L12261.16127.62-260.52-129.640.64-2.02L1395.5448.10-95.42-54.400.12-6.30L2480.5229.64-80.46-33.920.06-4.28L45-28.32-15.0828.3211.320-3.76L46-151.22-70.96151.5670.650.34-0.31L37-254.58-125.60255.24123.160.66-2.44L57-278.32-131.32

38、279.52128.661.20-2.66L67-301.56-140.66302.24139.880.68-0.78注：Lij表示i为线路首端，j为线路末端4）电网的网损率计算电网的网损为发电机发出的功率与平衡节点1与外界的交换功率的之和减去负荷功率之和。即公式（7） （7）其中，PD=发电机发出的有功功率+平衡节点与外界交换的有功功率，PL=负荷的有功率代数和，PD=3279+356.7=1193.7MW，总负荷PL=180+350+260+250+150=1190MW。于是网损率=（1193.7-1190）/1193.7=0.31%，与Matlab计算结果形同。2.3 Matlab仿真计

39、算结果与Power World仿真软件计算结果对比通过以上两种计算结果的比较，我们将其整理如下1）各节点电压的对比如表9所示：表9.各节点电压的对比节点编号标幺电压实际电压(kV)相角 (度)MatlabPower WorldMatlabPower WorldMatlabPower World11.04551.0455230.01230.010021.03751.0375228.25228.25-0.6033-0.6031.04171.0418229.17229.20-0.2606-0.2641.03581.0358227.88227.88-0.7560-0.7651.03631.036422

40、7.99228.00-0.7119-0.7161.04291.0430229.44229.46-0.1932-0.1971.05001.0500231.00231.000.36080.36比较两种计算结果的各节点电压可以发现，节点电压的最大偏差为0.03kV，误差极小。两种计算潮流得出的节点电压具有很好的一致性。说明两种计算方法都是可行的。2）各支路状态的对比如表10所示。表10、各支路状态的对比支路区间首端有功功率（MW）首端无功功率（Mvar）有功损耗（MW）无功损耗（Mvar）PowerWorldMatlabPowerWorldMatlabPowerWorldMatlabPowerWorldMatlab1，2261.16261.17127.620.47690.640.0065-2.02-0.02011，395.5495.5348.1047.690.120.0011-6.30-0.06302，480.5280.5229.6429.330.060.0005-4.28-0.0428