电力系统分析教学课件作者刘学军第4章节电力系统潮流的计算机算法课件幻灯片.ppt

电力系统分析（刘学军主编）机械工业出版社 解：第一步，确定节点类型：由已知条件发电机G1发出的有功、无功可调，可知其满足平衡节点的要求，因此节点5为平衡节点；发电机G2所在的节点4为PV节点；其余节点1、2、3均为PQ节点。 电力系统分析（刘学军主编）机械工业出版社 第二步，形成节点导纳矩阵： 由给出的系统图中的数据可以列出节点导纳矩阵如下 解：第一步，确定节点类型：由已知条件发电机G1发出的有功、无功可调可知其满足平衡节点的要求，因此节点5为平衡节点；G2所在的节点4为PV节点；其余节点1、2、3则为PQ节点。 因YB为对称阵，只示出了上三角部分，下同。 第三步，设定初值： = = = ， , 。 电力系统分析（刘学军主编）机械工业出版社 第四步，计算失配功率 显然， 电力系统分析（刘学军主编）机械工业出版社 第五步，形成雅可比矩阵（阶数为7×7） 第六步，解修正方程，得到 ， 电力系统分析（刘学军主编）机械工业出版社 从而 然后转入下一次迭代。经三次迭代后， 。 电力系统分析（刘学军主编）机械工业出版社 4.3.3牛顿-拉夫逊法潮流计算的相关问题 1．稀疏矩阵表示法 电力网络的节点导纳矩阵是一个高度稀疏的矩阵，也就是矩阵中很多非对角元素是零。如果不考虑导纳矩阵的稀疏性，导纳矩阵要用二维数组存放，那么矩阵中的大量零元素不仅要占用大量的计算机内存，而且还要由于这些零元素进行大量的不必要的运算，从而降低了计算效率。而在用牛顿-拉夫逊计算潮流计算时，雅可比矩阵具有与节点导纳矩阵相似的结构，如果不考虑雅可比矩阵的稀疏性，在每次迭代的雅可比矩阵的形成与分解计算中也必将发生对零元素的不必要的存取和计算。因此，不考虑稀疏性的牛顿-拉夫逊法不仅需要庞大的存贮空间，而且计算效率很低，速度很慢。 电力系统分析（刘学军主编）机械工业出版社 将稀疏矩阵技巧应用于导纳矩阵或雅可比矩阵时可将这些矩阵的对角元素存放在一个数组A中，而将矩阵的上三角矩阵的非对角元素的非零元素按行按列紧凑的压缩在一个一维数组B中由导纳或雅可比矩阵的性质可知，A数组元素的个数等于系统结点个数，B组数组原数的个数与系统不接地支路数相应（为支路数与消去过程中出现的非零注入元素之和）。因此，矩阵的稀疏存放极大地减少了所需存贮单元，所以，在以一个稀疏方式存贮一个n阶的稀疏方阵时大约需要6b+3n个存贮单元，其中b是系统非接地支路数，n为系统节点数。在典型的电力系统中，b1≈1.5n，所以总的存贮单元大约需要12n，与利用二维数组存放所需单元相比，对500个节点的系统，大约减少了40倍。同时。相应于零元素的额外计算工作也大大减少。 电力系统分析（刘学军主编）机械工业出版社 2．节点的优化编号 在电力系统的节点导纳矩阵或牛顿-拉夫逊法计算潮流的雅可比矩阵中的非对角元素，相应与矩阵的行、列号相应节点之间直接相连的支路。?如果电网节点编号改变，对应的节点导纳矩阵或雅可比矩阵的行，列号也将改变，所以出现非零注入元素的可能性和数目也将发生变化，或者说注入元与节点编号的排序有关。 显然，用不同的节点编号方案，所得到的注入元数目是不同的。为了在计算过程中减少内存容量和提高计算速度，要求找到一种或数种使节点注入元最少的节点编号方案，或最优的节点编号方案。事实上由于电力网节点数目很大，要找到注入元最少的最优的节点编号方案非常困难，实用上常常一准优化的节点编号方案来代替。准优化节点排序常用的方法有两种： （1）静态按最少出线支路数排序 （2）动态按最少出线支路排序 电力系统分析（刘学军主编）机械工业出版社 （1）静态按最少出线支路数排序 这种排序的原则是： 1）按节点所连支路数（不包括接地支路，两个节点间有多条并联支路时，只算作一条）的多少排序，连接之路数少的节点排在前面； 2）非零元素尽量安排在导纳矩阵的准对角线上或靠近对角线 应用这种静态优化法时，在导纳矩阵形成前，通过各节点连接支路多少的分析，可一次定出节点编号，因而计算工作量小。 （2）动态按最少出线支路排序 在高斯消去法按列消元的过程中，每消去一列所增加的注入元个数，相当于网络中消去与该列同号的节点，给剩余节点增加的相互连线数。 电力系统分析（刘学军主编）机械工业出版社 3．牛顿-拉夫逊法的收敛特性 牛顿-拉夫逊法具有平方收敛特性。它在开始时收敛得比较慢，而在几次迭代以后，收敛得特别快。由于他的这种收敛特性，在使用稀疏矩阵技巧和优化节点编号以后，在使用平直电压启动（设置全网电压初值幅值均为1.0pu，相角为0°）时，牛顿-拉夫逊法的迭代次数实际上与系统规模无关。一个设计良好的牛顿拉夫逊潮流计算程序的每一次迭代时间