电力系统的潮流分布计算计划书

1 电力系统的潮流分布计算计划书 第一章 电力系统潮流计算简述 电力系统的潮流分布计算（牛顿 拉夫逊直角坐标法） 第一节 潮流计算的意义 (1)在电网规划阶段 ,通过潮流计算 ,合理规划电源容量及接入点 ,合理规划网架 ,选择无功补偿方案 ,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。 (2)在编制年运行方式时 ,在预计负荷增长及新设备投运基础上 ,选择典型方式进行潮流计算 ,发现电网中薄弱环节 ,供调度员日常调度控制参考 ,并对规划、基建部门提出改进网架结构 ,加快基建进度的建议。 (3)正常检修及 特殊运行方式 下的潮流计算 ,用于日运行方式的编制 ,指导发电厂开机方式 ,有功、无功调整方案及负荷调整方案 ,满足线路、变压器热稳定要求及电压质量要求。 (4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 总结为在 电力系统运行方式 和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了 实时监控 电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在 电力系统运行状态 的实时监控中，则采用在线潮流计算。 第二节 潮流计算的发展史 利用电子计算机进行潮流计算从 20 世纪 50 年代中期就已经开始。此后，潮流计算曾采用了各种不同的方法，这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点： （ 1）算法的可靠性或收敛性 （ 2）计算速度和内存占用量 （ 3）计算的方便性和灵活性 电力系统潮流计算属于稳态分析范畴，不涉及系 统元件的动态特性和过渡过程。因此其数学模型不包含微分方程，是一组高阶非线性方程。非线性代数方程组的解法离不开迭代，因此，潮流计算方法首先要求它是能可靠的收敛，并给出正确答案。随着电力系统规模的不断扩大，潮流问题的方程式阶数越来越高，目前已达到几千阶甚至上万阶，对这样规模的方 2 程式并不是采用任何数学方法都能保证给出正确答案的。这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法。 在用数字计算机求解电力系统潮流问题的开始阶段，人们普遍采用以节点导纳 矩阵 为基础的高斯 个方法的原理比较简单，要求的数字计算机的内存量也比较小，适应当时的 电子数字计算机 制作水平和电力系统理论水平，于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法。 20 世纪 60 年代初，数字计算机已经发展到第二代，计算机的内存和计算速度发生了很大的飞跃，从而为阻 抗法的采用创造了条件。阻抗矩阵是满矩阵，阻抗法要求计算机储存表征系统接线和参数的阻抗矩阵。这就需要较大的内存量。而且阻抗法每迭代一次都要求顺次取阻抗矩阵中的每一个元素进行计算，因此，每次迭代的计算量很大。 阻抗法改善了电力系统潮流计算问题的收敛性，解决了导纳法无法解决的一些系统的潮流计算，在当时获得了广泛的应用，曾为我国电力系统设计、运行和研究作出了很大的贡献。但是，阻抗法的主要缺点就是占用计算机的内存很大，每次迭代的计算量很大。当系统不断扩大时，这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺 点，后来发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需存储各个地区系统的阻抗矩阵及它们之间的联络线的阻抗，这样不仅大幅度的节省了内存容量，同时也提高了计算速度。 克服阻抗法缺点的另一途径是采用牛顿 下简称 牛顿法 ）。牛顿法是数学中求解非线性方程式的典型方法，有较好的收敛性。解决电力系统潮流计算问题是以导纳矩阵为 基础的，因此，只要在迭代过程中尽可能保持方程式 系数矩阵 的稀疏性，就可以大大提高牛顿潮流程序的计算效率。自从 20 世纪 60 年代中期采用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、计算速度方面都超过了阻抗法，成为直到目前仍被广泛采用的方法。 在牛顿法的基础上，根据电力系统的特点，抓住主要矛盾，对纯数学的牛顿法进行了改造，得到了 解法。 解法在计算速度方面有显著的提高，迅速得到了推广。 牛顿法的特点是将非线性方程线性化。 20 世纪 70 年代后期，有人提出采用更精确的模型，即将泰勒级数的高阶项也包括进来，希望以此提高算法的性能，这便产生了保留非线性的潮流算法。另外，为了解决病态潮流计算，出现了将潮流计算表示为一个无约束非线性规划问题的模型，即非线性规划潮流算法。 近 20 多年来，潮流算法的研究仍然非常活跃，但是大多数研究都是围绕改进牛顿法和 解法进行的。此外，随着人工智能理论的发展，遗传算法、 人工神经网络 、模糊算法也逐渐被引入潮流计算。但是，到目前 3 为止这些新的模型和算法还不能取代牛顿法和 解法的地位。由于电力系统规模的不断扩大，对计算速度的要求不断提高，计算机的并行计算技术也将在潮流计算中得到广泛的应用，成为重要的研究领域。 第三节 潮流计算的发展趋势 通过几十年的发展，潮流算法日趋成熟。近几年，对潮流算法的研究仍然是如何改善传统的潮流算法，即高斯 顿法和快速解耦法。牛顿法，由于其在求解非线性潮流方程 时采用的是逐次线性化的方法，为了进一步提高算法的收敛性和计算速度，人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来，于是产生了二阶潮流算法。后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点，提出了采用直角坐标的保留非线性快速潮流算法。 对于一些病态系统，应用非线性潮流计算方法往往会造成计算过程的振荡或者不收敛，从数学上讲，非线性的潮流计算方程组本来就是无解的。这样，人们提出来了将潮流方程构造成一个函数，求此函数的最小值问题，称之为非线性规划潮流的计算方法。优点是原理上保证了计算过程永远 不会发散。如果将数学规划原理和牛顿潮流算法有机结合一起就是最优乘子法。另外，为了优化系统的运行，从所有以上的可行潮流解中挑选出满足一定指标要求的一个最佳方案就是最优潮流问题。最优潮流是一种同时考虑经济性和安全性的 电力网络 分析优化问题。 电力系统的安全运行、经济调度、可靠性分析、能量管理以及电力定价等方面得到了广泛的应用。 第二章 程序设计 第一节 潮流计算的基本步骤 点导纳矩阵 上图可得相应的节点导纳矩阵 Y= 5i 0 0; 5i i i i 0i 0; 4 0 i 0i 0 0 由于节点 1 为平衡节点，其他的节点均为 点，系数矩阵 B ， 应对该等值网络，不存在除去与有功功率和电压相位或无功功率和电压大小关系较小的因素的可能性， 这两个矩阵 B ， B 完全相同。他们就由导纳矩阵的虚部部分中除第一行第一列外的各个元素所组成 . 2 计算各节点有功功率不平衡量5U=1 =0； 02U = 03U= 04U = 01U=02 = 03= 04 = 05=0, 计算各节点有功功率不平衡量 0 0 00 s o s G 各节点电压的相位角i(弧度 ) 由矩阵方程式 - 1 B ( 00 = 00 U 3 计算各节点无功功率不平衡量下式计算各节点无功功率不平衡量 0 0 11 c o ss in G 由矩阵方程式计算各节点电压大小 - 1 B )/( 00 = 0U 5 4 计算平衡节点功率 51*= *iU*第二节 潮流计算 设计资料及参数 题目：（一）设计原始资料网络接线如图所示，各支路导纳均以标幺值标于图中。节点注入功率分别为： 1 0 . 2 0 0 . 2 0 ，2 0 . 4 5 0 . 1 5 ， 3 0 . 4 0 0 . 0 5 ， 4 0 . 6 0 0 . 1 0 ，其中节点 1 连接的实际上相当于给定功率的发电厂。设节点 5 电压保持定值， 5 。试运用以极坐标表示的 高斯 塞德尔法、牛顿 拉夫逊法、 算该系统的潮流分布。计算精度要求各节点电压修正量不大于510 。 0 8 + j 0 . 2 4 0 . 0 1 + j 0 . 0 0 4 + j 0 . 1 4 5 + j 0 . 1 5 0 . 4 + j 0 . 0 50 . 6 + j 0 . 1- ( 0 . 2 + j 0 . 2 )123 45 6 第三节 解法的程序流程图 ： 第四节 算设计题目潮流分布的完整的 序 节点总数为： ); N=5 平衡节点为： ); 1 点为： ); 2,3,4,5 7 Y= 5i 0 0; 5i i i i 0i 0; 0 i 0i 0 0 G=);B=); (2:5, 2:5); 1; b1=1); b2=2); )=0; )=0; )=0; )=0; u(2)=u(3)= 8 u(4)=u(5)=p(2)=p(3)=p(4)=p(5)=q(2)=q(3)=q(4)=q(5)=k=0; ; )=0;u(1)=m=2:5 n=1:5 pt(n)=u(m)*u(n)*(G(m,n)*dt(m)n)+B(m,n)*dt(m)n); 各节点有功率不平衡量为 ： ) dp(m)=p(m) 9 m=2:5 u1=u); u2=:5,2:5); u3=u(2:5); u4=:4,:); :5); :4,:); ,1:1); 各节点电压的相位角 ： ) dt(m)=dt(m)+m=2:5 n=1:5 qt(n)=u(m)*u(n)*(G(m,n)*dt(m)n)-B(m,n)*dt(m)n); 各节点无功率不平衡量为 ： ) qq(m)=q(m) 10 m=2:5 :5); :4,:); b2* 各节点电压的修正量 ： ) ut 各节点电压的大小为 ： ) u(m)=u(m)+ut( n=1:5 U(n)=u(n)*(dt(n)+j*dt(n); n=1:5 I(n)=Y(1,n)*U(n); 各节点功率 s 为 ： ) (1)*) m=1:5 n=1:5 S(m,n)=U(m)*(m)(n)*Y(m,n); 11 各节点功率 ： ) S k=k+1 五节 程序运行结果显示 节点总数为 ： N = 5 平衡节点为： 1 2 3 4 5 Y = 0 0 0 12 0 0 节点有功率不平衡量为： 0 节点有功率不平衡量为： 0 节点有功率不平衡量为： 0 节点有功率不平衡量为： 0 节点电压的相位角 0 0 0 0 节点电压的相位角 0 0 0 节点电压的相位角 0 0 13 各节点电压的相位角 0 节点无功率不平衡量为： 0 节点无功率不平衡量为： 0 节点无功率不平衡量为： 0 节点无功率不平衡量为： 0 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量 ： 0 0 0 0 0 0 14 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点功率 节点功率 S = 0 15 0 0 0 0 0 0 0 0 0 k = 1 各节点有功率不平衡量为： 0 节点有功率不平衡量为： 0 节点有功率不平衡量为： 0 节点有功率不平衡量为： 0 节点电压的相位角 16 0 节点电压的相位角 0 节点电压的相位角 0 节点电压的相位角 0 节点无功率不平衡量为： 0 节点无功率不平衡量为： 0 节点无功率不平衡量为： 0 节点无功率不平衡量为： 0 节点电压的修正量： 0 0 0 0 0 0 0 0 0 17 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 18 各节点电压的大小为： u = 节点功率 节点功率 S = 0 0 0 0 0 0 0 0 0 0 k =2 各节点有功率不平衡量为： 0 节点有功率不平衡量为： 0 19 各节点有功率不平衡量为： 0 节点有功率不平衡量为： 0 节点电压的相位角 0 节点电压的相位角 0 节点电压的相位角 0 节点电压的相位角 0 节点无功率不平衡 量为： 0 节点无功率不平衡量为： 0 节点无功率不平衡量为： 0 20 各节点无功率不平衡量为： 0 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 0 0 0 21 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点功率 节点功率 S = 0 0 0 0 0 0 22 0 0 0 0 k = 3 各节点有功率不平 衡量为： 0 节点有功率不平衡量为： 0 节点有功率不平衡量为： 0 节点有功率不平衡量为： 0 节点电压的相位角 0 节点电压的相位角 23 0 节点电压的相位角 0 节点电压的相位角 0 节点无功率不平衡量为： 0 节点无功率不平衡量为： 0 节点无功率不平衡量为： 0 节点无功率不平衡量为： 0 节点电压的修正量： 0 0 0 0 0 0 24 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点电压的修正量： 25 0 0 0 0 0 0 0 0 0 0 0 0 各节点电压的大小为： u = 节点功率 节点功率 S = 0 0 0 0 0 0