电力系统分析(-).pptVIP

现代电力系统分析 任课教师：葛少云 第六节 保留非线性潮流算法 牛顿法求解非线性潮流方程时采用了逐次线性化的方法。70年代后期，人们开始研究这样的问题，即如果采用更加精确的数学模型，将泰勒级数的高阶项或非线性项也考虑进来，也许会进一步提高算法的收敛性能及计算速度，于是便产生了一类称之为保留非线性的潮流算法。又因为其中大部算法主要包括了泰勒级数的前三项即取到泰勒级数的二阶项，所以也称为二阶潮流算法，实现这种想法的第一个尝试是在极坐标形式的牛顿法修正方程式中增加了泰勒级数的二阶项，所得到的算法对收敛性能略有改善，但计算速度无显著提高。 后来，参考文献[14]（岩本伸一及田村康男）根据直角坐标形式的潮流方程是一个二次代数方程组的这一特点，提出了采用直角坐标的保留非线性快速潮流算法，在速度上比牛顿法有较多的提高，引起了广泛的重视。在此以后，又出现了一些计入非线性的其它潮流算法。这些算法除了作为常规的潮流计算工具之外，也已经在状态估计、最优潮流等其它计算中得到应用。 一、保留非线性快速潮流算法 采用直角坐标形式的潮流方程为 该潮流问题实际上就是求解一个不含变量一次项的二次代数方程组。 对这样的方程组用泰勒级数展开，则二阶项系数已是常数，没有二阶以上的高阶项，所以泰勒级数只要取三项就能够得到一个没有截断误差的精确展开式。因此从理论上，假若能够从这个展开式设法求得变量的修正量，并将它对估计初值加以修正，则只要一步就可求得方程组的解。而牛顿法由于线性近似，略去了高阶项，因此用每次迭代所求得的修正量对上一次的估计值加以改进后;仅是向真值接近了一步而已。 （一）数学模型表达 我们先定义 为n 维函数给定值相量 x为n 维未知变量相量 一个具有n个变量的齐次二次代数方程式的普遍形式为: 于是潮流方程组就可以写成如下的矩阵形式 或 式中系数矩阵为 对上式在初值x（0）附近展开，可得到如下没有截断误差的精确展开式 于是与上式对应的精确的泰勒展开式（1-74）为 式中 为修正量相量 H是一个常数矩阵，其阶数很高，但高度稀疏。注意若式(1—74)中略去第三项，就成为通常的牛顿法的展开式。 研究表明可以进一步将式(1—74)改写为 上式的推出促成了本算法的突破，因为可以非常方便地计算二阶项。值得顺便指出的是，该式是一个很重要的关系式，在研究其它算法时将多次引用。 （二）数值计算迭代公式 上式是一个以 作为变量的二次代数方程组，从一定的初值 x（0）出发，求解满足该式的 仍然要采用迭代的方法。上式可改写成： 于是算法的具体迭代公式为： k表示迭代次数；J为按x=x(0)估计得到；在进行第一次迭代时，k=0,令 ，于是和牛顿法的第一次迭代计算完全相同。 算法的收敛判据为算法的收敛判据为 也可以采用相继二次迭代的二阶项之差，作为收敛判据，即 (四)算法特点及性能估计 保留非线性快速潮流算法的特点可以通过和牛顿法进行比较而得以揭示。 设求解的方程是 ，则牛顿法的迭代公式是 此算法的迭代公式是 图1—9表示了两种算法的迭代过程。 由两组公式可见，与牛顿法的在迭代过程中变化的雅可比矩阵不同，保留非线性快速潮流算法采用的是用初值 x（0）计算而得的恒定雅可比矩阵，整个计算过程只需一次形成，并三角分解构成因子表。 就每一次迭代所需的计算量而言，牛顿法要重新计算y(x(k))，而保留非线性快速潮流算法则要计算y(△x(k))，由于计算函数式完全相同，仅变量不同、所以这部分的计算量是完全相同的，但由于保留非线性快速潮流算法不需重新形成雅可比矩阵并三角分解，所以每次迭代所需时间可以大大节省。 两种算法的△x的含义是不同的。牛顿法的△x(k)是相对于上一次迭代所得到的迭代点x(k)的修正量；而保留非线性快速潮流算法的△x(k)则是相对于始终不变的初始估计值x(0)的修正量。图l-9中AA1：对应于y(x(0))-ys)，A1A2、A1A3 、…、 A1An对应于逐次迭代中变化着的二阶项y(△x)；逐次迭代就对应于求解一系列相似三角形，平行的斜边说明用的是和第一次迭代同样的恒定不变的雅可比矩阵。 保留非线性快速潮流算法达到收敛所需的迭代次数比牛顿法要多，在半对数坐标纸上收敛特性近似为一条直线。但由于每次迭代所需的计算量比牛顿法节省很多，所以总的计算速度比牛顿法可提高很多。 由于不具对称性质的雅可比矩阵经三角分解后，其上下三角元素都需要保存,和牛顿法的仅需保存上三角元素相比较，此算法所需的矩阵存储量将较后者要增加35％一40％ 另外，由于利用以初始值计算得到的恒定雅可比矩阵进行迭代