1、1 绪 论11.1 课题背景11.2 电力系统潮流计算的意义21.3 电力系统潮流计算的发展21.4 潮流计算的发展趋势42 潮流计算的数学模型52.1 电力线路的数学模型及其应用52.2 等值双绕组变压器模型及其应用62.3 电力网络的数学模型82.4 节点导纳矩阵92.4.1 节点导纳矩阵的形成92.4.2 节点导纳矩阵的修改102.5 潮流计算节点的类型112.6 节点功率方程1227 潮流计算的约束条件133 牛顿拉夫逊法潮流计算基本原理143.1 牛顿-拉夫逊法的基本原理143.2 牛顿-拉夫逊法潮流计算的修正方程173.3 潮流计算的基本特点203.4 节点功率方程214牛顿拉夫逊
2、法分解潮流程序2241 牛顿拉夫逊法分解潮流程序原理总框图224.2 形成节点导纳矩阵程序框图及代码234.2。1 形成节点导纳矩阵程序框图234.2.2 形成节点导纳矩阵的程序代码243 雅克比矩阵求取的程序框图及代码2531 形成节点导纳矩阵程序框图252 形成雅克比矩阵程序的代码254 求取Df、De的程序框图及代码281 求取Df、De的程序框图282 求取Df、De的程序代码285 实例与分析305.1 一个6节点算例3052 根据算例输入相应节点的线路参数315.3 算例运行335.3.1 原始数据的输入335.3.2 原始数据输入程序段345.3.3 导纳矩阵的形成355.3.4
3、 导纳矩阵Y为365.3.5 雅克比矩阵J（k=0）365.3.6 算例运行结果输出程序段365.3.7 算例运行结果输出385.4 潮流计算GUI界面40全文总结45参考文献46致谢47附录A 潮流计算程序48附录B“清除数据”按钮回调函数57附录C“关闭”按钮回调函数58河南城建学院本科毕业设计（论文） 1 绪论1 绪 论1.1 课题背景潮流计算是研究电力系统的一种最基本和最重要的计算，最初，电力系统潮流计算是通过人工手算的，后来为了适应电力系统日益发展的需要，采用了交流计算台。随着电子数字计算机的出现，1956年Ward等人编制了实际可行的计算机潮流计算程序。这样，就为日趋复杂的大规模电
4、力系统提供了极其有利的计算手段。经过几十年的时间，电力系统潮流计算已经发展的十分成熟。潮流计算是研究电力系统稳态运行情况的一种计算，是根据给定的运行的条件及系统接线情况确定整个电力系统各个部分运行的状态，如各母线的电压、各元件中流过的电流、系统的功率损耗等等。电力系统潮流计算是计算系统动态稳定和静态稳定的基础。在电力系统规划设计和现有电力系统运行的方式研究中，都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性。 电力系统潮流计算分为离线计算和在线计算，离线计算主要用于系统规划设计、安排系统的运行方式，在线计算则用于运算中系统的实时监测和实施控制。两种计算的原理在本
5、质上是相同的。实际电力系统的潮流技术主要采用牛顿拉夫逊法。牛顿拉夫逊法早在50年代末就已应用于求解电力系统潮流问题，但作为一种适用的、有竞争力的电力系统潮流计算方法，则是在应用了稀疏矩阵技巧和高斯消元法求修正方程式以后。牛顿拉夫逊法是求解非线性代数方程有效的迭代计算，本设计就是采用牛顿拉夫逊法计算电力系统潮流的。1.2 电力系统潮流计算的意义 （1）在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划网架，选择无功补偿方案，满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。 （2）在编制年运行方式时，在预计负荷增长及新设备投运基础上，选择典型方式进行潮流计算，发现电网
6、中薄弱环节，供调度员日常调度控制参考，并对规划、基建部门提出改进网架结构,加快基建进度的建议。 （3）正常检修及特殊运行方式下的潮流计算，用于日运行方式的编制，指导发电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求及电压质量要求。 （4）预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 总之在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划
7、设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。1.3 电力系统潮流计算的发展 利用电子计算机进行潮流计算从20世纪50年代中期就已经开始。此后，潮流计算曾采用了各种不同的方法，这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点： （1）算法的可靠性或收敛性 （2）计算速度和内存占用量 （3）计算的方便性和灵活性 电力系统潮流计算属于稳态分析范畴，不涉及系统元件的动态特性和过渡过程。因此其数学模型不包含微分方程，是一组高阶非线性方程。非线性代数方程组的解法离不开迭代，因此，潮流计算方法首先要求它是能可靠
8、的收敛，并给出正确答案。随着电力系统规模的不断扩大，潮流问题的方程式阶数越来越高，目前已达到几千阶甚至上万阶，对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的，这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法。 在用数字计算机求解电力系统潮流问题的开始阶段，人们普遍采用以节点导纳矩阵为基础的高斯-赛德尔迭代法（以下简称导纳法）。这个方法的原理比较简单，要求的数字计算机的内存量也比较小，适应当时的电子数字计算机制作水平和电力系统理论水平，于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法（简称阻抗法）。 20世纪60年代初，数字计算机已经发展到第二代，计算机的内存和计算速
9、度发生了很大的飞跃，从而为阻抗法的采用创造了条件。阻抗矩阵是满矩阵，阻抗法要求计算机储存表征系统接线和参数的阻抗矩阵。这就需要较大的内存量，而且阻抗法每迭代一次都要求顺次取阻抗矩阵中的每一个元素进行计算，因此，每次迭代的计算量很大。 阻抗法改善了电力系统潮流计算问题的收敛性，解决了导纳法无法解决的一些系统的潮流计算，在当时获得了广泛的应用，曾为我国电力系统设计、运行和研究作出了很大的贡献。但是，阻抗法的主要缺点就是占用计算机的内存很大，每次迭代的计算量很大。当系统不断扩大时，这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺点，后来发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分
10、割为几个小的地区系统，在计算机内只需存储各个地区系统的阻抗矩阵及它们之间的联络线的阻抗，这样不仅大幅度的节省了内存容量，同时也提高了节省速度。 克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法（以下简称牛顿法）。牛顿法是数学中求解非线性方程式的典型方法，有较好的收敛性。解决电力系统潮流计算问题是以导纳矩阵为基础的，因此，只要在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿潮流程序的计算效率。自从20世纪60年代中期采用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、计算速度方面都超过了阻抗法，成为直到目前仍被广泛采用的方法。 在牛顿法的基础上，根据电力系统的特点，抓住主要矛盾，对纯
11、数学的牛顿法进行了改造，得到了P-Q分解法。P-Q分解法在计算速度方面有显著的提高，迅速得到了推广。 牛顿法的特点是将非线性方程线性化。20世纪70年代后期，有人提出采用更精确的模型，即将泰勒级数的高阶项也包括进来，希望以此提高算法的性能，这便产生了保留非线性的潮流算法。另外，为了解决病态潮流计算，出现了将潮流计算表示为一个无约束非线性规划问题的模型，即非线性规划潮流算法。 近20多年来，潮流算法的研究仍然非常活跃，但是大多数研究都是围绕改进牛顿法和P-Q分解法进行的。此外，随着人工智能理论的发展，遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。但是，到目前为止这些新的模型和算法还不能取
12、代牛顿法和P-Q分解法的地位。由于电力系统规模的不断扩大，对计算速度的要求不断提高，计算机的并行计算技术也将在潮流计算中得到广泛的应用，成为重要的研究领域。1.4 潮流计算的发展趋势 通过几十年的发展，潮流算法日趋成熟。近几年，对潮流算法的研究仍然是如何改善传统的潮流算法，即高斯-塞德尔法、牛顿法和快速解耦法。牛顿法，由于其在求解非线性潮流方程时采用的是逐次线性化的方法，为了进一步提高算法的收敛性和计算速度，人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来，于是产生了二阶潮流算法。后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点，提出了采用直角坐标的保留非线性快速潮流算法。4河
13、南城建学院本科毕业设计（论文） 2 潮流计算的数学模型2 潮流计算的数学模型2.1 电力线路的数学模型及其应用 在电力系统稳态分析中的电力线路数学模型就是以电阻、电抗、电纳、电导表示的它们的等值电路。 式（2.1）式中为导线材料的电阻率（mm2/km）；s为导线的额定截面积（mm2）。 式（2.2）式中r为导线计算半径（mm或cm）；Dm为几何均距（mm或cm），其单位应与r的相同。 式（2.3） 式（2.4）式中b1导线单位长度的电纳（S/km）； g1导线单位长度的电导（S/km）；三相线路泄漏和电晕损耗功率（kW/km）；U线路线电压（kV）。按上式求得单位长度导线的电阻、电抗、电纳、电
14、导后，就可作最原始的电力线路等值电路图，如图2-1所示。这是单相等值电路。之所以可用单相等值电路代表三相，一方面由于本文中讨论的是三相对称运行方式，另一方面也因设架空线路都已经整循环换位。图2-1 中等线路等值模型以单相等值电路代表三相虽已简化了不少计算，但由于电力线路的长度往往有数十乃至数百公里，如将每公里的电阻、电抗、电纳、电导都一一绘于图上，所得的等值电路仍十分复杂。何况，严格说来，电力线路的参数并不是均匀分布的，即使是极短的一段线段，都有相应大小的电阻、电抗、电纳、电导。换言之，即使是如此复杂的等值电路，也不能认为精确。但好在电力线路一般都不长，需分析的又往往只是它们的端点状况两端电压
15、、电流、功率，通常可不考虑线路的这种分布参数特性，只是在个别情况下才要用双曲函数研究具有均匀分布参数的线路。以下，先讨论一般线路的等值电路。中等长度的线路通常指100km-300km之间的架空线路，这种线路的导纳一般不能略去，常用的是型等值电路。 当线路长度为（km）时：2.2 等值双绕组变压器模型及其应用无论采用有名制或标幺制，凡涉及多电压级网络的计算都必须将网络中所有参数和变量归算至同一电压等级。这是因为以型或T型等值电路做变压器模型时，这些等值电路模型并不能体现变压器实际具有的电压变换功能。以下将介绍另一种可等值的体现变压器电压变换功能的模型，它也是运用计算机进行电力系统分析时采用的变压
16、器模型，虽然运用这种模型时并不排斥手算。既然这种模型可体现电压变换，在多电压级网络计算中采用这种变压器模型后，就可不必进行参数和变量的归算，这正是这种变压器模型的主要特点之一。以下就介绍这种变压器模型。 如图2-2所示。 图2-2 等值双绕组变压器 首先，从一个未作电压级归算的简单网络入手。设图中变压器的导纳或励磁支路和线路的导纳支路都可略去；设变压器两侧线路的阻抗都未经归算，即分别为高、低电压侧、侧线路的实际阻抗，变压器本身的阻抗则归在低压侧；设变压器的变比k，其值为高、低压绕组电压之比。显然，在这些假设条件下，如在变压器阻抗ZT的左侧串联以变比为K的理想变压器如图（2-2c）所示，其效果就
17、如同将变压器及其低压侧线路的阻抗都归算至高压侧，或将高压侧线路的阻抗归算至低压侧，从而实际上获得将所有参数和变量都归算在同一侧的等值网络，只要变压器的变比取的是实际变比，这一等值网络无疑是严格的。由图（2-2c）可见流入理想变压器的功率为，流出理想变压器的功率为，流入流出变压器的功率应该相等，可得： 式（2.5）从而有： 式（2.6）另外由图2-2c可以直接得到： 式（2.7）联立解方程组： 式（2.8）可得： 式（2.9）即： 式（2.10）的成立体现了无源电路的互易特性，然后令，就可以作导纳支路表示的变压器模型如图（2-2e）所示以及以阻抗支路表示的变压器模型如图（2-2f）所示。其中，。
18、以下利用图2-2说明各种不同情况下等值变压器模型的应用，即多电压级网络中变压器和线路参数的计算，以及相应的理想变压器变比的取值。（1）有名值、线路参数都归算到低压侧。据以图2-2的情况，由图可见，此时线路阻抗分别为上图中，变压器阻抗则由；相应的理想变压器变比则为，这里所取得理想变压器的变比就是变压器实际变比。（2）有名值、线路参数都归算到高压侧。这种情况下的线路阻抗分别为从而理想变压器变比为： 2.3 电力网络的数学模型 有名制：所有参数和变量都以有名单位，如、S、kV（V）、kA（A）、MVA（VA）等表示。标幺制：所有参数和变量都以与他们同名基准值相对的标幺值表示，因此都没有单位。对多电压
19、级网络，变压器模型：采用等值变压器模型时，所有参数和变量可不进行归算；采用有名制或标幺制取决于习惯。在我国，电力工程界使用标幺值已有多年；但在国外，有名制的使用也很普遍。至于变压器模型的使用范围，则泾渭分明。手算时，都是用形或T型等值电路模型；计算机计算时，都是用等值变压器或型等值电路模型。 此外，在制定电力网络等值电路模型时，有时还同时作某些简化，常见的有：线路的电导通常都被略去；变压器的电导有时以具有定值的有功功率损耗的形式出现在电路中；100km以下架空线路的电纳被略去；100300km架空线路或变压器的电纳有时以具有定值的容性或感性无功功率损耗的形式出现在电路中。有时，整个元件，甚至部
20、分系统都可能不包括在等值电路中。例如，将某些发电厂的高压母线看作为可维持给定电压、输出给定功率的等值电源时，这些发电厂内部的元件就不再包括在等值电路中。2.4 节点导纳矩阵在电路原理课程中，已导出了运用节点导纳矩阵的节点电压方程 式（2.11）上式中，是节点注入电流的列向量，可理解为某个节点的电源电流与负荷电流之和，并规定电源流向网络的注入电流为正。因此，仅有负荷的负荷节点注入电流就具有负值。是节点电压的列向量。因通常以大地作参考节点，网络中有接地支路时，节点电压通常就指该节点的对地电压；网络中没有接地支路时，各节点电压可指各该节点与某一个被选定参考节点之间的电压差。是一个节点导纳矩阵，它的阶
21、数n等于网络中除参考节点外的节点数。它可展开为 式（2.12）2.4.1 节点导纳矩阵的形成 根据定义直接求取节点导纳矩阵时，注意以下几点：（1） 节点导纳矩阵是方阵，其阶数就等于网络中除去参考节点外的节点数。参考节点一般取大地，编号为零。 （2） 节点导纳矩阵是稀疏矩阵，其各行非零非对角元素就等于与该行相对应节点所连接的不接地支路数。 （3） 节点导纳矩阵的对角元素就等于该节点所连接导纳的总和。因此，与没有接地支路的节点对应的行或列中，对角元素为非对角元素之和的负值。 （4） 节点导纳矩阵的非对角元素等于连接节点i，j支路导纳的负值。因此，在一般情况下，节点导纳矩阵的对角元素往往大于非对角元
22、素的负值。 （5） 节点导纳矩阵一般是对称矩阵，这是网络的互易特性所决定的。从而，一般只要求取这个矩阵的上三角或下三角部分。（6） 网络中的变压器。2.4.2 节点导纳矩阵的修改 （1） 从原有网络引出一支路，同时增加一节点。设i为原有网络中的节点，j为新增加的节点，新增加支路导纳为。则因新增一节点，节点导纳矩阵将增加一阶。 新增的对角元，由于在节点j上只有一个支路，将为=；新增的非对角元；原有矩阵中的对角元将增加。（2） 在原有网络的节点i、j之间增加一支路。这时由于仅增加支路不增加节点，节点导纳矩阵阶数不变，但与节点i、j有关的元素应作一下修改，其增量为:（3） 在原有网络的节点i，j之间切除一支路。切除一导纳为的支路，相当于增加一导纳为的支路，从而与节点i、j有关的元素应作如下修改： （4） 原有网络
