1-3及4高斯法和牛顿法.

任课教师：王守相天津大学电气与自动化工程学院2008年现代电力系统分析现代电力系统分析1.3高斯一塞德尔法潮流以导纳矩阵为基础，并应用高斯--塞德尔迭代的算法是在电力系统中最早得到应用的潮流计算方法。优点：原理简单，程序设计十分容易。导纳矩阵是一个对称且高度稀疏的矩阵，因此占用内存非常节省。就每次迭代所需的计算量而言，是各种潮流算法中最小的，并且和网络所包含的节点数成正比关系。现代电力系统分析缺点：本算法的主要缺点是收敛速度很慢。病态条件系统，计算往往会发生收敛困难节点间相位角差很大的重负荷系统；包含有负电抗支路(如某些三绕组变压器或线路串联电容等)的系统；具有较长的辐射形线路的系统；长线路与短线路接在同一节点上，而且长短线路的长度比值又很大的系统。此外，平衡节点所在位置的不同选择，也会影响到收敛性能。目前高斯一塞德尔法已很少使用现代电力系统分析1.4牛顿一拉夫逊法1.4.1牛顿一拉夫逊法的一般概念牛顿一拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程，即通常所称的逐次线性化过程。现代电力系统分析对于非线性代数方程组即在待求量x的某一个初始估计值x(0)附近，将上式展开成泰勒级数并略去二阶及以上的高阶项，得到如下的经线性化的方程组0)(xf),...2,1(0),...,,(21nixxxfni0)()()0()0()0(xxfxf现代电力系统分析上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量将相加，得到变量的第一次改进值x(1)。接着就从x(1)出发，重复上述计算过程。因此从一定的初值x(0)出发，应用牛顿法求解的迭代格式为:)()]([)0(1)0()0(xfxfx)0()0(xx和)()(1)()()()()('kkkkkkxxxxfxxf）（现代电力系统分析上两式中：f’(x)是函数f(x)对于变量x的一阶偏导数矩阵，即雅可比矩阵J；k为迭代次数。由上两式可见，牛顿法的核心便是反复形成并求解修正方程式。牛顿法当初始估计值x(0)和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。现代电力系统分析1.4.2牛顿潮流算法的修正方程式在将牛顿法用于求解电力系统潮流计算问题时，由于所采用f(x)的数学表达式以及复数电压变量采用的坐标形式的不同，可以形成牛顿潮流算法的不同形式。以下讨论用得最为广泛的f(x)采用功率方程式模型，而电压变量则分别采用极坐标和直角坐标的两种形式。现代电力系统分析(一)极坐标形式令则采用极坐标形式的潮流方程是：对每个PQ节点及PV节点对每个PQ节点iiiUU现代电力系统分析将上述方程式在某个近似解附近用泰勒级数展开，并略去二阶及以上的高阶项后，得到以矩阵形式表示的修正方程式为:式中：n为节点总数；m为PV节点数，雅可比矩阵是(2n-m-2)阶非奇异方阵。现代电力系统分析(二)直角坐标形式令在这里，潮流方程的组成与上不同，对每个节点，都有二个方程式，所以在不计入平衡节点方程式的情况下，总共有2(n-1)个方程式。iiijfeU现代电力系统分析对每个PQ节点，根据式(1—11)和式(1—12)有：现代电力系统分析对每个PV节点，除了有与式(1-39)相同的有功功率方程式之外，还有采用直角坐标形式的修正方程式为0)()(2222iiisiUfeU现代电力系统分析仔细分析以上两种类型的修正方程式，可以看出两者具有以下的共同特点。(1)修正方程式的数目分别为2(n-1)-m及2(n-1)个，在PV节点所占比例不大时，两者的方程式数目基本接近2(n-1)个。(2)雅可比矩阵的元素都是节点电压的函数，每次迭代，雅可比矩阵都需要重新形成。现代电力系统分析(3)分析雅可比矩阵的非对角元素的表示式可见，某个非对角元素是否为零决定于相应的节点导纳矩阵元素Yij是否为零。因此如将修正方程式按节点号的次序排列，并将雅可比矩阵分块，把每个2X2阶子阵；作为分块矩阵的元素，则按节点号顺序而构成的分块雅可比矩阵将和节点导纳矩阵具有同样的稀疏结构，是一个高度稀疏的矩阵。(4)和节点导纳矩阵具有相同稀疏结构的分块雅可比矩阵在位置上对称，但雅可比矩阵不是对称阵。复习并分析这些特点非常重要，因为正是修正方程式的这些特点决定了牛顿法潮流程序的主要轮廓及程序特色。现代电力系统分析1.4.3修正方程式的处理和求解在本节的开头就已提到，牛顿算法的核心就是反复形成并求解修正方程式。因此如何有效地处理修正方程式就成为提高牛顿法潮流程序计算速度并降低内存需量的关键所在。?现代电力系统分析从算法的发展过程来看，在50年代末就已经提出了牛顿法潮流的雏形。先是用迭代法求解修正方程式，但遇到迭代法本身不收敛的问题。用高斯消去法等直接法求解，但如前所分析，修正方程式的数目在2(n-1)左右，如果不利用雅可比矩阵的稀疏特性，当网络节点数增加为N倍，存储雅可比矩阵的内存量将正比于N2倍,利用直接法求解修正方程的计算量将正比于N3倍地增长。这就限制了牛顿法潮流程序的解题规模，从而使得这种方法的推广应用一度止步不前。现代电力系统分析其后正是人们注意到了雅可比矩阵高度稀疏的特点，求解修正方程式时采用了稀疏程序设计技巧，并且发展了一套在消元过程中旨在尽量保持其稀疏性、以减少内存需量并提高计算速度的有效方法(即著名的最优顺序消去法)，才使牛顿法真正得到了突破，因而在60年代中期以后被普遍采用。现代电力系统分析结合修正方程式的求解，目前在实用的牛顿法潮流程序中所包含的程序特点主要有以下三个方面，这些程序特点对牛顿法潮流程序性能的提高起着决定性的作用。(1)对于稀疏矩阵，在计算机中以“压缩”方式只储存其非零元素，且只有非零元素才参加运算。现代电力系统分析(2)修正方程式的求解过程，采用对包括了修正方程常数项的增广矩阵以按行消去而不是传统的按列消去的方式进行消元运算。由于消元运算系按行进行，因此可以不需先形成整个增广矩阵，然后进行消元运算，而是采取边形成、边消元、边存储的方式，即每形成增广矩阵的一行便马上进行消元，并且消元结束后便随即将结果送内存存储。现代电力系统分析图1-1是增广矩阵按行消元的示意图，图中表示了五阶增广矩阵的前四行，其中1-3行已完成了消元运算且已经存放在内存中，接着要进行的是第四行的消元运算，即消去对角元以左的三个元素。在具体的程序中，待消行是放在一个专用的工作数组中进行消元运算的。现代电力系统分析这种按行消元做法的好处:是对于消元过程中新注入的非零元素，当采用“压缩”存储方式时，可以方便地按序送入内存，不需要预留它们的存放位置。特别值得注意的是由于不必一次形成整个雅可比矩阵，且常数项的消元运算已和矩阵的消元过程同时进行，因此这种牛顿潮流算法求解修正方程式时，所需的矩阵存储量只是消元运算结束时所得到的用以进行回代的上三角矩阵而已。现代电力系统分析(3)消元的最优顺序或节点编号优化经过消元运算得到的上三角矩阵一般仍属稀疏阵，但由于消元过程中在原来是零元素的位置上有新元素注入，使得它的稀疏度比原来雅可比矩阵的上三角有所降低。但分析表明，注入元素的多少和消元的顺序或节点编号有关。节点编号优化的作用即在于找到一种网络节点的重新编号方案，使得按此构成的节点导纳矩阵以及和它相应的雅可比矩阵在高斯消元或三角分解过程中出现的注入元素数目能大大减少。节点编号优化通常有三种方法：现代电力系统分析静态法——按各节点静态连接支路数的多少顺序编号；半动态法——按各节点动态连接支路数的多少顺序编号；动态法——按各节点动态增加支路数的多少顺序编号。三种节点编号优化方法：动态法效果最好，但优化本身所需计算量也最多，而静态法则反之。对于牛顿法潮流计算来说，一般认为，采用半动态法似乎是较好的选择。现代电力系统分析1.4.4牛顿潮流算法的性能和特点牛顿潮流算法突出的优点是收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代4—5次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。牛顿法也具有良好的收敛可靠性，对于上节中提到的对以节点导纳矩阵为基础的高斯一塞德尔法呈病态的系统，牛顿法均能可靠地收敛。现代电力系统分析牛顿法所需的内存量及每次迭代所需时间均较前述的高斯一塞德尔法为多，并与程序设计技巧有密切关系。牛顿法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有可能根本不收敛或收敛到一个无法运行的解点上。对于正常运行的系统，各节点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值(也称为“平直电压”)，现代电力系统分析“平直电压”法假定：或这样一般能得到满意的结果。但若系统因无功紧张或其它原因导致电压质量很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问题。0100iiU);,...,2,1(0100sinifeii现代电力系统分析解决这个问题的办法可以先用上一节的高斯一塞德尔法迭代1-2次；以此迭代结果作为牛顿法的初值。也可以先用直流法潮流求解一次以求得一个较好的角度初值，然后转入牛顿法迭代。