电力系统潮流计算的计算机仿真

1、毕业设计（论文）材料之二（1）安徽工程大学本科安徽工程大学本科毕业设计（论文）毕业设计（论文）专专 业：业： 电气工程及其自动化电气工程及其自动化 题题 目：目： 电力系统潮流计算的电力系统潮流计算的 计算机仿真计算机仿真 作作 者者 姓姓 名名： 张三张三 导导师师及及职职称称 ： 黄梅初（讲师）黄梅初（讲师） 导师所在单位：导师所在单位： 电气工程学院电气工程学院 2015 年年 6 月月 2 日日-安徽工程大学安徽工程大学本科毕业设计（论文）任务书本科毕业设计（论文）任务书 20152015 届届 电气工程电气工程 学院学院 电气工程及其自动化电气工程及其自动化 专业专业学生姓名：学生姓

2、名： 张三张三 毕业设计（论文）题目毕业设计（论文）题目中文：电力系统潮流计算的计算机仿真中文：电力系统潮流计算的计算机仿真英文：英文：computer analysis of load flow calculation for power system 原始资料原始资料 （也可以附参考文献）（也可以附参考文献）五节点等值电力网络如上图所示。图中，节点 1 为平衡节点，其它节gg523141s2s3s4s点都是 pq 节点，保持为定值，给定各支路阻抗：006. 1.5ju，18. 006. 02112jzz18. 006. 03113jzz12. 004. 04114jzz，06. 002.

3、05115jzz03. 001. 03223jzz24. 008. 05225jzz，也给定各节点输出功率：，24. 008. 04334jzz2 . 02 . 01js，。15. 045. 02js05. 04 . 03js1 . 06 . 04js试运用以极坐标表示的牛顿拉夫逊法计算此五节点等值网络中的潮流分布。计算精确度要求各节点不平衡量不大于。510 毕业设计（论文）任务内容毕业设计（论文）任务内容1、课题研究的意义电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算，又是研究电力系统的一项重要分析功能，同时也是进行故障计算、继电保护鉴定、安全分析的工具。电力系统潮流计算

4、是计算系统动态稳定和静态稳定的基础。在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性。2、本课题研究的主要内容：掌握简单电力网络潮流的手工计算，包括简单支路、辐射形网络和闭合环网的潮流计算。掌握复杂电力系统潮流计算的数学模型，包括电力网络方程的基本概念，功率方程及变量、节点的分类。熟练掌握 matlab软件的基本功能应用，包括基本运算和程序设计。熟悉目前流行的复杂电力系统潮流的计算机算法，掌握高斯赛德尔迭代法、牛顿拉夫逊法和pq 分解法的基本原理。结合具体算例，画出潮流计算的流程框图，运用matlab 软件编制应

5、用程序进行计算机分析，并对程序进行调试校验，运算结果包含各节点注入功率、节点注入电流、节点电压、线路流动功率、网络功率损耗和迭代次数等。-3、提交的成果：（1）毕业设计（论文）正文；（2）主要源程序：（3）至少一篇引用的外文文献及其译文；（4）附不少于 10 篇主要参考文献的题录及摘要。指导教师（签字） 教研室主任 （ 签 字 ） 批 准 日 期 2015年1月4日 接 受 任 务 书 日 期 2015年1月10日 完 成 日 期 2015年6月2日 接受任务书学生（签字） 安徽工程大学毕业设计（论文）-i-电力系统潮流计算的计算机仿真电力系统潮流计算的计算机仿真摘摘 要要（摘要写成三段式，具

6、体见论文格式的要求）（摘要写成三段式，具体见论文格式的要求）电力系统潮流计算是电力系统各种计算的基础，又是研究电力系统的一项重要分析功能，同时也是进行故障计算、继电保护鉴定、安全分析的工具。电力系统潮流计算是计算系统动态稳定和静态稳定的基础。在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用电力系统潮流计算来定量地比较供电方案或运行方式的合理性、可靠性和经济性。本文首先介绍了简单电力网络潮流的手工计算和复杂电力系统潮流计算的数学模型，包括简单支路、辐射形网络和闭合环网的潮流计算，电力网络方程的基本概念，功率方程以及变量、节点的分类。然后介绍到目前流行的复杂电力系统潮流的计算机算法，重

7、点详述高斯赛德尔法、牛顿拉夫逊法和 pq 分解法的基本原理。最后，本文采用一个五节点等值电力网络进行实例分析，画出潮流计算的流程框图，运用matlab 软件编制应用程序进行计算机分析，运算内容包含各节点注入功率、节点注入电流、节点电压、线路流动功率、网络功率损耗和迭代次数等，计算结果表明程序的算法具有良好的收敛性和实用性。关键词：关键词：潮流计算；复杂电力系统；matlab；牛顿拉夫逊法张三：电力系统潮流计算的计算机仿真-ii-computer analysis of load flow calculation for power systemabstractpower flow calcul

8、ation is the foundation of power system calculation.it is not only an important analysis function for power system, but also a tool for fault calculation, identification of relay protection and security analysis. power flow calculation is the basic calculation of dynamic and static stabilization. in

9、 the research of power system planning and the existing operation mode, it is based on power flow calculation to quantificationally compare rationality, reliability and economy of power supply scheme or operation mode.first of all,this paper introduces the manual calculation of simple power flow and

10、 the mathematical model of complex power flow calculation, including load flow calculation of simple branches, radial network and closed loop network, the basic concept of power network equation, power equation, as well as the classification of variable and node. then this paper introduces popular c

11、omputer algorithm for complex power flow at present, which focus on the details of the basic theory of gauss - seidel method, newton - raphson method and pq decomposition method. finally, this paper uses a power network equivalence of five nodes to carry on the example analysis, which has drawn the

12、flow diagram of load flow calculation, using matlab software applications for computer analysis. its operation contains the nodal injection power, nodal injection current, node voltage, line flow power, network power loss and the number of iterations. and the calculation results show that the algori

13、thm has good convergence and practicability.keywords: load flow calculation; complex power system; matlab; newton - raphson method安徽工程大学毕业设计（论文）-iii-目目 录录(每章和每节标题按要求设置好后，目录最后能自动生成每章和每节标题按要求设置好后，目录最后能自动生成)引引 言言.1第第 1 章章 绪论绪论.21.1 课题的学术背景与研究意义.21.2 潮流计算的国内外发展概况.21.3 课题研究的主要内容.3第第 2 章章 简单电力网络潮流的手工计算简单电

14、力网络潮流的手工计算.42.1 简单辐射网络的潮流计算.42.1.1 简单支路的潮流分布.42.1.2 辐射型网络的潮流计算方法.62.2 简单环形网络的潮流计算.92.2.1 两端电压相等.92.2.2 两端电压不相等.92.3 手工潮流算例分析.10第第 3 章章 复杂电力系统潮流计算的数学模型复杂电力系统潮流计算的数学模型.133.1 电力网络的数学模型.133.1.1 节点电压方程.133.1.2 节点导纳矩阵的形成和修改.143.2 功率方程及变量、节点的分类.153.2.1 功率方程.153.2.2 变量的分类.163.2.3 节点的分类.18第第 4 章章 复杂电力系统潮流的计算

15、机算法复杂电力系统潮流的计算机算法.194.1 高斯赛德尔法潮流计算.194.2 牛顿拉夫逊法潮流计算.224.2.1 牛顿拉夫逊法的基本原理.224.2.2 直角坐标形式的牛顿拉夫逊法.244.2.3 极坐标形式的牛顿拉夫逊法.274.3 pq 分解法潮流计算.29第第 5 章章 潮流计算实例的计算机辅助分析潮流计算实例的计算机辅助分析.335.1 matlab 软件简介.335.2 给定的等值电力网络.345.3 牛顿拉夫逊算法的程序设计.345.4 潮流计算的运算结果.35结论与展望结论与展望.38致致 谢谢.39参考文献参考文献.40附附 录录.42附录 a 主要源程序.42张三：电力

16、系统潮流计算的计算机仿真-iv-附录 b 外文文献及其译文.46附录 c 主要参考文献的题录及摘要.58安徽工程大学毕业设计（论文）-v-插图清单插图清单图 2-1 简单支路示意图.4图 2-2 两端电压相量示意图.5图 2-3 单电源辐射型电力系统.6图 2-4 等值电路 i.7图 2-5 等值电路 ii.8图 2-6 等值电路 iii.8图 2-7 简单环形网等效两端供电网.9图 2-8 两端电压不等的网与环网等值.9图 2-9 10kv 配电网络.10图 3-1 电力系统的等值网络图.13图 3-2 简单系统及其等值网络.15图 4-1 线路上流通的电流和功率.21图 4-2 直角坐标形

17、式的牛顿拉夫逊潮流计算流程图.27图 4-3 pq 分解法潮流计算的流程框图.32图 5-1 五节点等值电力网络.34图 5-2 极坐标形式的牛顿拉夫逊法潮流计算流程图.35张三：电力系统潮流计算的计算机仿真-vi-表格清单表格清单表 5-1 节点注入功率.36表 5-2 节点注入电流.36表 5-3 节点电压.36表 5-4 各线路上的流动功率.36表 5-5 各线路上的功率损耗.37安徽工程大学毕业设计（论文）-1-引引 言言电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算，又是研究电力系统的一项重要分析功能，同时也是进行故障计算、继电保护鉴定、安全分析的工具。电力系统潮

18、流计算是计算系统动态稳定和静态稳定的基础。在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性。潮流计算的目的在于：确定电力系统的运行方式；检查系统中的各元件是否过压或过载；为电力系统继电保护的整定提供依据；为电力系统的稳定计算提供初值，为电力系统规划和经济运行提供分析的基础。因此，电力系统潮流计算是电力系统中一项最基本的计算，既具有一定的独立性，又是研究其他问题的基础。随着计算机技术的不断发展和成熟，基于 matlab 的计算机潮流计算研究近年来得到了长足的发展，本文最后的实例分析将在 matlab 软件的开发环境

19、下设计出一个能够在电力系统中完成潮流分析的计算程序软件，在保证电力系统供电可靠性和电能质量的前提下，根据给定的电力系统网络图及网络基本参数，运用一种具有较好的收敛性和实用性的潮流算法自行完成电力系统潮流计算这种电力系统基本计算，尽可能提高潮流计算的效率，降低人力资源消耗，从而提高电力系统运行的经济性。从引言开始，就可以在全文标注参考文献，如：随着电力系统容量日益增大，范围越来越广，仅设置系统各元件的继电保护装置，远不能防止发生全电力系统长期大面积停电的严重事故。为此必须从电力系统全局出发，研究故障元件被相应继电保护装置的动作切除后，系统将呈现何种工况，系统失去稳定时将出现何种特征，如何尽快恢复

20、其正常运行等。系统保护的任务就是当大电力系统正常运行被破坏时，尽可能将其影响范围限制到最小，负荷停电时间减到最短15。并且从序号为 1 的参考文献开始标注起，后面依次是 2,3,4 等等，在你引用内容的右上角标注，如上所示。张三：电力系统潮流计算的计算机仿真-2-第第 1 章章 绪论绪论1.1 课题的学术背景与研究意义电力系统潮流计算的分析与比较是建立在计算机潮流计算的基础上进行分析与比较的。潮流计算的计算机求解始于 20 世纪 50 年代。50 年代中期，潮流计算主要采用以节点导纳拒阵为基础的高斯赛德尔迭代法。该方法原理简单，占用计算机内存少，适应当时的电子数字计算机制作水平和电力系统理论水

21、平，但高斯赛德尔迭代法收敛性差，当系统规模变大时，迭代次数急剧上升，且常有不收敏的情况。于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法（以下简称阻抗法） 。20 世纪 60 年代初，数字计算机已经发展到第二代，计算机的内存和计算速度发生了很大的飞跃，从而为阻抗法的采用创造了条件。阻抗法改善了收敛性，但因占用内存多，使解题规模受到一定限制。60 年代初期即开始研究牛顿拉夫逊法。研究表明，牛顿拉夫逊法具有很好的收敛性。直到 60 年代末期，优化节点编号和稀疏矩阵程序技巧的高斯消去法的实际应用，才使牛顿拉夫逊法潮流计算在收敏性、内存需求、计算速度等方面都超过其他方法，成为广泛采用的优秀方法。70

22、 年代初，在牛顿拉夫逊法的基础上，根据电力系统的特点发展了潮流计算 pq 分解法。该方法所占内存约为牛顿拉夫逊法的1/21/4，计算速度也明显加快。由于牛顿拉夫逊法和 pq 分解法的显著优点，使得到 90 年代末期为止，它们仍然是实际应用的电力系统潮流计算的主要方法。此外，作为方法的研究和探讨，还提出了非线性快速潮流计算法、最优乘子法、非线性规划法、网流法等。为适应电力网调度自动化的需要，在线潮流计算方法及其应用也得到重视和发展。70 年代后期至 80 年代以来，大型商用电力系统分析软件包得到广泛应用，其中不少潮流计算程序包，同时备有几种算法供用户选择，以便有助于解决各类潮流计算的收敛性问题。

23、电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算。其任务是要在给定某些运行参数的情况下，计算出系统中的全部参数，包括各母线电压的大小和相位、各发电机和负荷的功率和电流以及各个变压器和线路等元件所通过的功率、电流和损耗。潮流计算又是研究电力系统的一项重要分析功能，同时也是进行故障计算、继电保护鉴定、安全分析的工具。电力系统潮流计算是计算系统动态稳定和静态稳定的基础。在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性。潮流计算的目的在于：确定电力系统的运行方式；检查系统中的各元件是否过压或过载；为电力

24、系统继电保护的整定提供依据；为电力系统的稳定计算提供初值，为电力系统规划和经济运行提供分析的基础。因此，电力系统潮流计算是电力系统中一项最基本的计算，既具有一定的独立性，又是研究其他问题的基础。1.2 潮流计算的国内外发展概况近 20 多年来，电力系统潮流算法的研究仍然非常活跃，但是大多数研究都是围绕改进牛顿拉夫逊法和 pq 分解法进行的。此外，随着人工智能理论的发展，遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。但是，到目前为止这些新的模型和算法还不能取代牛顿拉夫逊法和 pq 分解法的地位。由于电力系统规模的不断扩大，对计算速度的要求不断提高，计算机的并行计算技术也将在潮流计算中得到

25、广泛的应用，成为重要的研究领域。此外，通过几十年的发展，潮流算法日趋成熟，但对潮流安徽工程大学毕业设计（论文）-3-算法的研究仍然是如何改善传统的潮流算法，即高斯塞德尔法、牛顿拉夫逊法和快速解耦法。另外针对潮流计算在不同领域的应用，研究人员提出了更为专业，适应性强的潮流计算算法。文献1提出一种以牛顿拉夫逊法为基础的 zbus 潮流计算算法，经过推导得出电流不匹配量，并把网络方程虚部和实部的分开写，最后再通过牛顿拉夫逊法求解非线性方程，即通过修改雅可比矩阵元素来处理 pv 节点的算法。文献2提出电流型牛顿拉夫逊法的一般形式，雅可比矩阵(y+d)部分在每次迭代过程中只有对角块发生不大的变化，因此易

26、于通过编程来实现；同时在处理 pv 节点的时候把无功功率 q 作为状态变量被引入节点功率方程,该方法方便了 pv、pq 两种类型节点之间的相互转换。该算法具有较好的收敛性，不过在修正方程式的时候，有可能维数增加带来内存开销增加和运算量也会有所增加，但系数矩阵增加的部分仍然保持高度稀疏性，因此该算法在使用稀疏技术和节点优化编号后，内存的开销增加量和运算过程中的计算量增加就不会太明显。文献3提出改进 pq 分解法，将小阻抗、零阻抗和 r/x 比值较大的问题用节点等效附加注入法进行统一处理。改进后的 pq 分解法能够满足多种状态下的潮流计算，却能够保持着 pq 分解法的高效特点。文献4提出一种以配电

27、网络特有的层次结构特性为基础的新颖的分层前推回代算法。该算法采用支路分层计算，得出各支路的电压降以及其功率的损耗，从而提高了潮流计算的运算速度。文献5提出树状网的交替迭代的潮流计算算法，即以节点电压和线路功率损耗量交替迭代，在迭代过程中首先计算出各节点的电压，随后算出每条线路支路的功率损耗量，从而最终求出每条线路的支路功率。该算法对小阻抗输电元件、r/x 值较大的元件，以及重负荷输电线路都具有相当好的收敛性。文献6分析了全局潮流计算的特殊困难以及要满足的基本要求，并建立了潮流计算的数学模型，使用了较为严格的全局潮流计算主从分裂法。该算法将庞大的电力系统潮流分成为发电、输电潮流以及一系列规模较小

28、的配电馈线潮流子系统，因此通过该算法使得潮流计算规模得到很大的降低。文献7在静态电压稳定的连续潮流算法的基础上然后结合 pq 分解法潮流计算的特点，提出了以 pq 分解法为基础的连续潮流法，该算法使得连续潮流法的计算速度得到进一步提高。文献8提出以雅可比矩阵和快速分解法为基础，同时与 pq 分解连续潮流算法相结合的算法，即在预测环节采用完整的雅可比矩阵来求取切向量，避免使用 b与 b时产生的误差，不过在初始潮流计算、校正环节仍然采用 pq 分解法进行求解。1.3 课题研究的主要内容掌握简单电力网络潮流的手工计算，包括简单支路、辐射形网络和闭合环网的潮流计算。掌握复杂电力系统潮流计算的数学模型，

29、包括电力网络方程的基本概念，功率方程及变量、节点的分类；熟悉节点电压方程、节点导纳矩阵的特点及其应用。熟练掌握 matlab 软件的基本功能应用，包括基本运算和程序设计。熟悉目前流行的复杂电力系统潮流的计算机算法，掌握高斯赛德尔迭代法、牛顿拉夫逊法和 pq 分解法的基本原理。最后，结合具体算例，画出了潮流计算的流程框图，运用 matlab 软件编制应用程序进行计算机分析，并对程序进行调试校验，运算结果包含各节点注入功率、节点注入电流、节点电压、线路流动功率、网络功率损耗和迭代次数等。张三：电力系统潮流计算的计算机仿真-4-安徽工程大学毕业设计（论文）-5-第第 2 章章 简单电力网络潮流的手工

30、简单电力网络潮流的手工计算计算2.1 简单辐射网络的潮流计算大约半个多世纪以前，数字计算机还没有出现的时候，潮流计算都是采用手工的计算方法。虽然潮流计算的本质是解电力系统的节点功率方程，然而手工的计算方法是不可能用解上述节点功率方程的方法来进行潮流计算的。手工潮流计算是根据一个简单支路的电压和功率传输关系，将较为复杂的电力系统分解为若干个简单支路来进行潮流计算的。因此，任何复杂的潮流计算都是从一个简单支路的潮流分布和电压降落的计算开始的。2.1.1 简单支路的潮流分布如图 2-1 所示的简单支路，节点 1 和 2 之间的阻抗为已知；两端的电压jxrz分别为和，从节点 1 注入该支路的复功率为，

31、从节点 2 流出的功率为，阻1u2u1s2s抗消耗的功率为。根据电路理论，、和、这四个变量，任何两个变量已s1u1s2u2s知都可以求出另外两个变量9,10。jx12r1s2s1u2u图 2-1 简单支路示意图1、已知一侧的电压和功率求另一侧的电压和功率假设已知节点 2 的电压和流出的功率，可知道流过该支路的电流为：2u2s （2-1）22usi如果以作为参考相量，阻抗引起的电压降落和功率损耗分别为：2uz （2-2）222)()(ujqpjxrud （2-3）2222222)(uqpjxrzis因此另一端节点 1 的电压为： （2-4）uduu21张三：电力系统潮流计算的计算机仿真-6-流过

32、节点 1 的复功率为： （2-5）sss21两端电压的关系还可以从如图 2-2 所示的相量图中得到（以为参考相量） ，为2.u末端电压和电流的夹角，称为功率因数角。从相量图中，不难得到阻抗引起的电压z降落的横分量和纵分量分别为： （2-6）222222222222sincossincossincossincosurqxpuiruixurixiuuxqrpuixuiruxiriu可得到首端的电压幅值和相角分别为： （2-7）2221)(uuuu （2-8）uuuarctg2如果已知首端（节点 1）的电压和功率，求末端的电压和功率，其基本原理同上，读者可以自行推导分析。2uiirijx1u图 2-

33、2 两端电压相量示意图2、已知一端的电压和流过另一端的复功率假如已知首端电压和末端的功率，要求首端的功率和末端的电压，我1u2s1s2u们可以利用两端电压的关系以及两端功率的关系列出如下方程组（以为参考相量）：1.u （2-9）jxruqpss22222221 （2-10）11111112urqxpjuxqrpuu直接求解上面这个相量方程组是很麻烦的，可以通过迭代法来求解：先给定一个末端电压的初值，这个初值可以设定为该节点的平均额定电压，然后将之代入式（2-安徽工程大学毕业设计（论文）-7-9） ，得到，然后再利用根据式（2-10）得到，重复上面的过程，直到误差满足1s1s2u要求为止。由于潮

34、流计算通常是在电力系统的稳态运行条件下，此时节点电压与平均额定电压差别不大。因此，在手工近似计算中，将上述的迭代过程只进行一次，即先设定未知的电压为平均额定电压，利用式（2-3） ，根据末端的功率计算支路的功率损耗，然后利用式（2-5）计算出首端的功率，再利用首端的功率和首端的电压计算系统的电压损耗，最后计算出末端的电压。2.1.2 辐射型网络的潮流计算方法所谓辐射型网络就是单电源供电的非环形网络，系统中所有的负荷都由一个电源供电，辐射型网络是由若干个简单支路树枝状串级联接而成的。对于辐射型网络中的接地支路可以做如下处理：（1）将对电力系统中的接地支路等效为该支路消耗的功率，对地支路的电压用额

35、定电压来替代。例如，对地支路的导纳为，那么这个对地支路消耗的功率为jbg ；2)(nujbgs（2）将同一节点消耗的功率进行合并。通过这样处理，辐射型网络就化减为若干简单支路的级联，可以利用简单支路的潮流和电压计算方法逐级进行潮流计算。辐射型网络的手工潮流计算一般从系统末端开始，因为通常辐射型网络的末端的负荷为已知，首先计算潮流的近似分布，然后再从电源端开始根据潮流分布计算出各个节点的电压。因此，辐射型网络的手工潮流估算仅包含三步：第一步，根据电力系统各个元件的电气参数，建立电力系统的等值计算电路，接着将对地支路等效为支路消耗的功率，并将各个节点消耗的功率进行合并。第二步，首先将系统中各个节点

36、的未知电压设为系统平均额定电压，然后从辐射型网络的末端开始，依次计算各个支路的功率损耗，最后得到潮流在辐射型网络中的近似分布。第三步，根据估算出的潮流分布，从电源端开始，根据前面简单支路的电压计算公式依次计算各个节点的电压。下面通过一个实例来说明潮流计算的过程，如图 2-3 所示的辐射型单电源的简单电力系统，已知节点 1（发电机节点）的电压和各个节点的负荷、au1ls2ls3ls，求该系统的功率和电压的分布9,10。4lst-1g2t-2134l1u1ls2ls3ls4ls图 2-3 单电源辐射型电力系统已知电力系统的各个元件的参数如下所示：张三：电力系统潮流计算的计算机仿真-8-变压器 t-

37、1：额定容量，额定变比，空载损耗，空载电流百nsniinituuk/10p分数，短路损耗，短路电压百分数；%0ikp%ku输电线路 l：每公里长的正序阻抗，每公里长的对地电纳，线路长度 ；1z0bl变压器 t-2：额定容量，额定变比，空载损耗，空载电流百nsniinituuk/20p分数，短路损耗，短路电压百分数；%0ikp%ku第一步作出等效电路及其参数：首先做电力系统的等值电路，根据上述各个元件的参数，我们可以得到各个元件的等效电路及其电路参数，等效电路 i 如图 2-4 所示。在计算等值电路中各个元件参数之前，先选择功率和电压的基准值、。bs1bu2bu3bu变压器 t1（根据等值电路，

38、变压器参数都归算到高压侧）： （2-11）211*1*1*1*120*12220*1*1*1*1222*12222*1/kkk100%100%bbniinibtttttniintbbniittttbbnniiktbbnniiktuuuujbgyusibsuupgjxrzussuuxussupr；输电线路： （2-12）bblbblsulbbuslzz220*0221*；变压器 t2（根据等值电路，变压器参数都归算到高压侧）： （2-13）322*2*2*2*220*22220*2*2*2*2222*22222*2/kkk100%100%bbniinibtttttnintbbnittttbbnn

39、iktbbnniktuuuujbgyusibsuupgjxrzussuuxussupr；1u1ls2ls3ls4ls123as1tk2tk1tz2tzlz1ty2ty20lbj20lbj4图 2-4 等值电路 i第二步，将对地支路简化为对地功率损耗：安徽工程大学毕业设计（论文）-9-如果电压基准值的选取与变压器的实际变比相匹配，那么，如果不1*2*1ttkk匹配，则需要将变压器的变比的标么值等效到电路中，把变压器的阻抗支路，变为型等效电路。为了说明问题，我们假设电压基准值选取与变压器实际变比匹配，或者忽略非标准变比的影响。对地支路假设为对地损耗功率，其对地支路的损耗用该点的额定电压来计算，等效

40、电路变为如图 2-5 所示。1u1ls2ls3ls4ls1234as1tz2tzlz1ts1ls2ts2ls图 2-5 等值电路 ii其中：*22*4*2*02*3*2*02*2*1*22*2*1)2/()2/(tntlnllnltntyusjbusjbusyus；第三步，节点功率合并：然后，将 1、2、3、4 各个节点上的所有功率合并，如图 2-6 所示：1u1234as1tz2tzlz1s2s3s4s2s3s4s2s3s4s2s 3s 4s 图 2-6 等值电路 iii其中：*1*1*1*1*2*2*2*2*3*3*4*4ltlltlllssssssssssss；第四步，从末端开始，根据末

41、端功率计算功率分布：先用各个节点的额定电压以及流出支路的功率来计算各个支路损耗以及功率分布：*1*2*2*2*2*1*12*22*2*2*2*3*2*3*3*3*2*32*3*3*3*4*3*4*4*4*2*22*42*4*4)()()(ssssssjxrussssssssjxrussssssssjxrussattnllnttn ；这样，就求得了功率的分布和节点 1 的注入功率。as张三：电力系统潮流计算的计算机仿真-10-第五步，从首端开始，根据首端电压计算电压损耗和各个节点的电压：*4*3*4*3*2*4*2*4*3*2*4*2*4*4*3*2*3*2*3*3*2*3*3*3*2*1*2*

42、1*1*2*1*2*1*1*2*1*2*2uuuurqxpjuxqrpuuuuurqxpjuxqrpuuuuurqxpjuxqrputtttlllltttt ；2.2 简单环形网络的潮流计算环形网可以等效成两端供电网，两端供电网也可以等效成环形网10,11。2.2.1 两端电压相等如图 2-7 所示，可将（a）图等效成（b）图。12311231u11uuasbs12z31z23z2s3sas12z2s3sbs23z31z（a） （b） 图 2-7 简单环形网等效两端供电网（a）环形网；（b）两端供电网 （2-14）zzszzzzszzssmma)(312312313312322.2.2 两端电

43、压不相等两端电压不相等的网络，可以等效成回路电压不为零的单一环网，如图 2-8 所示。ascs2s3sud4u14uu12z23z34z231uas12z2s3s23z141u34z（a） （b） 图 2-8 两端电压不等的网与环网等值安徽工程大学毕业设计（论文）-11-（a）两端供电网；（b）环形网 （2-15）34231234231234334232)(zzzuudzzzzszzssna其中称为循环功率。cnnszuudzzzuud342312对环形网的潮流分布，首先求出、，然后求各支路上的流动功率，即初步的asbs潮流分布，没有计及网络各段的电压降落、功率损耗。初步潮流分布的目的，在于找

44、出功率分点，以便在功率分点把闭环网打开成两个辐射网。然后，以功率分点为末端，对这两个辐射网分别用逐段推算法进行潮流分布计算。从中要计及各段的电压降落和功率损耗，所运用的公式与计算辐射网时完全相同。在两端供电网中，当两端电压相量不等，不论是模值还是相位不等都将产生循环功率。在环网中，循环功率是由于环网中有多台变压器，而变压器的变比不匹配引起的。所谓变比不匹配则是指环网中有两台及以上变压器时，由于变压器变比的不同使得网络空载且开环时开口两侧有电压差，即开口两侧感应电势不同，因而闭环后，即使空载也有环路电流，产生循环功率。应该特别注意正确地确定环网中循环功率的方向。循环功率的正方向取决于电压降落的正

45、方向。环网和两端供电网中的循环功率可改变网络中功率的分布。2.3 手工潮流算例分析10kv 配电网络的电网结构如图 2-9 所示。已知各节点的负荷功率及线路参数如下：=1.2+j2.4，=1.0+j2.0，=1.5+j3.0，=0.3+j0.2mva，=0.5+j0.3mva12z23z24z2s3s，=0.2+j0.15mva。设母线 1 的电压为 10.5kv，线路始端功率容许误差为 0.3%。4s110.5kvu 12444jqp 24z22jqp 33jqp 12z23z2u3u12s23s24s4u312s图 2-9 10kv 配电网络首先，假设各节点电压均为额定电压，功率损耗计算的

46、支路顺序为 3-2、4-2、2-1，第一轮计算依上列支路顺序计算各支路的功率损耗和功率分布。张三：电力系统潮流计算的计算机仿真-12-mvajjjxruqpsmvajjjxruqpsnn0019. 00009. 0)35 . 1 (1015. 02 . 0)(0068. 00034. 0)21 (103 . 05 . 0)(2222424224242422223232232323则mvajssssmvajsssmvajsss6587. 00043. 11519. 02009. 03068. 05034. 022423122442423323又mvajjjxruqpsn0346. 00173.

47、0)4 . 22 . 1 (106587. 00043. 1)(2221212221221212mvajsss6933. 00216. 1121212然后，用已知的线路始端电压=10.5kv 及上述求得的线路功率、，1u12s24s23s按上列相反的顺序求出线路各点电压，计算中忽略电压降落横分量。kvuuuuxqrpukvuuuuxqrpukvuuuuxqrpu1912.100336. 01508.100740. 02248.102752. 0232322323232323242422424242424121211212121212根据上述求得的线路各点电压，重新计算各线路的功率损耗和功率分布

48、。mvajjsmvajjs0018. 00009. 0)35 . 1 (1508.1015. 02 . 00065. 00033. 0)21 (1912.103 . 05 . 02222422223故mvajsssmvajsss1518. 02009. 03065. 05033. 02442423323则mvajssss6583. 00042. 12242312又mvajjs0331. 00166. 0)4 . 22 . 1 (2248.106583. 00042. 122212安徽工程大学毕业设计（论文）-13-从而可得线路始端功率mvajs6914. 00208. 112经过上述两轮迭代计

49、算，结果与第一步所得的计算结果比较相差小于 0.3%，计算到此结束。最后一次迭代结果可作为最终计算结果。张三：电力系统潮流计算的计算机仿真-14-第第 3 章章 复杂电力系统潮流计算的数学模型复杂电力系统潮流计算的数学模型3.1 电力网络的数学模型电力网络的数学模型指的是将网络的有关参数和变量及其相互关系归纳起来所组成的、可以反映网络性能的数学方程式组。在利用计算机的复杂电力系统潮流计算中，用的最多的是节点电压方程。节点电压方程用节点电压和导纳来表示支路电流，根据基尔霍夫电流定律列出方程组。3.1.1 节点电压方程如图 3-1 所示是一个电力系统的等值网络图，它共有三个节点。1231u2u3u

50、1i2i10y20y30y13y23y图 3-1 电力系统的等值网络图根据基尔霍夫电流定律，对该电路图列写节点电压方程得： （3-1）333232131332322212123132121111uyuyuyiuyuyuyiuyuyuyi式（3-1）就是图 3-1 所示电力网络等值电路的数学模型。将其用矩阵形式表式为： （3-2）bbbuyi式（3-2）可以展开为 （3-3）321333231232221131211321uuuyyyyyyyyyiii在式（3.2）中是节点注入电流的列向量。在电力系统计算中，节点注入电流可bi以理解为与该节点相连的正电流源与负电流源之和，其中规定注入该节点的电流

51、为正，流出该节点的电流为负。有的节点不与电流源相连，则其注入电流为零，如图 3-1 中的节点 3。是节点电压列向量，节点电压是该节点对参考地的电压。是一个buby阶的节点导纳矩阵，其阶数就等于网络中除参考地之外的节点数。nnn安徽工程大学毕业设计（论文）-15-3.1.2 节点导纳矩阵的形成和修改1、节点导纳矩阵的形成将电力网络的拓扑化为如图 3-1 所示的等值电路后，就可以进行节点导纳矩阵的计算了。求取节点导纳矩阵时，需注意以下几点9,12：（1）节点导纳矩阵是方阵，其阶数就等于网络中除参考地之外的节点数。n（2）节点导纳矩阵是稀疏矩阵，其各行非零非对角元数就等于与该行相对应节点所连接的不接

52、地支路数。如图 3-1 所示，与节点 2 对应的第二行非零非对角元数为 2。（3）节点导纳矩阵的对角元，称为自导纳，其数值上就等于与该节点所连接iiy导纳的总和。如图 3-1 中，与节点 2 对应的对角元。23212022yyyy（4）节点导纳矩阵的非对角元，称为互导纳，其数值上就等于连接节点 、ijyi支路导纳的负值。如图 3-1 所示，、。j2121yy2323yy（5）节点导纳矩阵一般是对称矩阵，这是网络的互易特性所决定的。2、节点导纳矩阵的修改在电力系统计算中，往往要计算不同结线方式下的运行状况。以下介绍几种典型的修改方法9,13：（1）从原有网络引出一支路，同时增加一节点，节点导纳矩

53、阵将增加一阶。新增的对角元，；新增的非对角元，；原有矩阵中的jjyijjjyyijyijjiijyyy对角元将增加，。iiyiiyijiiyy （2）在原有网络的节点 、之间增加一支路，应作如下修改ij；ijiiyy ijjjyyijjiijyyy（3）在原有网络的节点 、之间切除一支路，应作如下修改ij；ijiiyyijjjyyijjiijyyy（4）原有网络的节点 、之间的导纳由改变为，应作如下修改ijijyijy；ijijiiyyyijijjjyyyijijjiijyyyy（5）原有网络节点 、之间变压器的变比由改变为，应作如下修改ijk k；tiiykky22110jjytjiijyk

54、kyy11张三：电力系统潮流计算的计算机仿真-16-3.2 功率方程及变量、节点的分类3.2.1 功率方程设有简单系统如图 3-2 所示9。图中，、分别为母线 1、2 的等值电源功率；1gs2gs、分别为母线 1、2 的等值负荷功率；它们的合成、1ls2ls111lgsss分别为母线 1、2 的注入功率，与之对应的电流、222lgsss111lgiii分别为母线 1、2 的注入电流。于是222lgiii （3-4）222222212112121111usuyuyiusuyuyi （3-5）121222222212111111uyuuyusuyuuyusg1u2u12yg210y20yg1u2u

55、g21111llljqps222llljqps111gggjqps222gggjqps111llljqps222llljqps111gggjqps222gggjqps（a） （b） 1u2u12y210y20y1111lgsss222lgsss111lgiii222lgiii（c）图 3-2 简单系统及其等值网络（a）简单系统；（b）简单系统的等值网络；（c）注入功率和注入电流如令安徽工程大学毕业设计（论文）-17- （3-6）212211)90(21122112)90(212012102211jjjmjseuueuueyyyyyeyyyyyyyms；并将式（3-6）代入式（3-5）展开，将有

56、功功率和无功功率分列，可得这个简单系统的功率方程为 （3-7）)cos(cos)cos(cos)sin(sin)sin(sin121222222212121111121222222212121111mmsslgmmsslgmmsslgmmsslguuyuyqqquuyuyqqquuyuypppuuyuyppp将式（3-7）中的第一、二式相加，第三、四式相加，又可得这个系统的有功、无功功率平衡关系 （3-mmssllggmmssllgguuyuuyqqqquuyuuyppppcos)cos(2cos)(sin)cos(2sin)(2121222121212121222121218）由式（3-7）

57、可见，在功率方程中，母线电压的相位角是以差值的形式出)(21现的，亦即决定功率大小的是相对相位角或相对功率角，而不是绝对相位1221角或绝对功率角或。12由式（3-8）可见，这个简单系统的有功、无功功率损耗分别为mmssmmssuuyuuyquuyuuypcos)cos(2cos)(sin)cos(2sin)(2121222121212221它们都是母线电压，和相位角，相对相位角或的非线性函数。1u2u12123.2.2 变量的分类由式（3-4）还可见，在这四个一组的功率方程组中，除网络参数、symys外，共有十二个变量，它们是14：m负荷消耗的有功、无功功率、；1lp1lq2lp2lq电源发

58、出的有功、无功功率、；1gp1gq2gp2gq母线或节点的大小和相位角、。1u2u12因此，除非已知其中的八个变量，否则将无法求解。在这十二个变量中，负荷消耗的有功、无功功率无法控制，因它们取决于用户。张三：电力系统潮流计算的计算机仿真-18-它们就称为不可控变量或扰动变量。之所以称扰动变量是由于这些变量出现事先没有预计的变动时，系统将偏离它们的原始运行状况。不可控变量或扰动变量用列向量表示。余下的八个变量中，电源发出的有功、无功功率是可以控制的自变量，因而d它们就称为控制变量。控制变量用列向量表示。最后余下的四个变量母线电压u或节点电压的大小和相位角是受控制变量控制的因变量。其中、主要受、1

59、u2u1gq的控制，、主要受、的控制。这四个变量就是这简单系统的状态变量。2gq121gp2gp状态变量一般用列向量表示。x变量的这种分类也适用于个节点的复杂系统。只是对于这种复杂的系统，变量n数将增加到个，其中扰动变量、控制变量、状态变量各个。换言之，扰动向量、n6n2d控制向量、状态向量，都是阶的列向量。uxn2看来似乎将变量做如上分类后，只要已知或给定扰动变量和控制变量，就可以运用功率方程式（3-7）解出状态变量。其实不然，因已如上述，功率方程中，母线或节点电压的相位角是以相对值出现的，以致式（3-7）中的、变化同样大小时，功率12的数值不变，从而不可能运用它们求取绝对相位角。也如上所述

60、，系统的功率损耗本身就是状态变量的函数，在解得状态变量前，不可能确定这些功率损耗，从而也不可能按功率平衡关系式（3-8）给定所有控制变量，因它们的总和，如式（3-8）中的（+） 、 （+）尚属未知。1gp2gp1gq2gq为克服上述困难，可对变量的给定稍作调整：在一个有个节点的系统中，只给定对控制变量、，余下的一对控n) 1( ngipgiq制变量、待定。这一对控制变量、将使系统功率，包括电源功率、负荷gspgsqgspgsq功率和损耗功率保持平衡。在这系统中，给定一对状态变量、，只要求确定sus对状态变量、。给定的通常为零。这实际上就是相当于取节点 的电压) 1( niuiss相量为参考轴。