3.第二章 力系统潮流计算.doc

第六节 保留非线性潮流算法 牛顿法求解非线性潮流方程时采用了逐次线性化的方法。70年代后期，人们开始研究这样的问题，即如果采用更加精确的数学模型，将泰勒级数的高阶项或非线性项也考虑进来，也许会进一步提高算法的收敛性能及计算速度，于是便产生了一类称之为保留非线性的潮流算法。又因为其中大部算法主要包括了泰勒级数的前三项即取到泰勒级数的二阶项，所以也称为二阶潮流算法，实现这种想法的第一个尝试是在极坐标形式的牛顿法修正方程式中增加了泰勒级数的二阶项，所得到的算法对收敛性能略有改善，但计算速度无显著提高。 后来，日本早稻田大学岩本伸一及田村康男根据直角坐标形式的潮流方程是一个二次代数方程组的这一特点，提出了采用直角坐标的保留非线性快速潮流算法，在速度上比牛顿法有较多的提高，引起了广泛的重视。在此以后，又出现了一些计入非线性的其它潮流算法。这些算法除了作为常规的潮流计算工具之外，也已经在状态估计、最优潮流等其它计算中得到应用。 一、二次型函数的泰勒展开 考虑只包含二次项的代数函数为如下形式 （2-78） 该函数三阶以上的导数为零，对各自变量的一、二次偏导数为 （2-79） 该函数在，处泰勒展开为 （2-79） 即二次型函数的泰勒展开二次项部分与原函数有同样的形式，只是变量换成了变量的增量。由于三阶以上的导数为零，所以（2-79）是函数的精确表达。 因此从理论上讲，如（2-78）是方程的左端，方程的右端为一常数，则如能够从展开式（2-78）设法求得变量的修正量，并将它对估计初值加以修正，则只要一步就可求得方程组的解。而牛顿法由于线性近似，略去了高阶项，因此用每次迭代所求得的修正量对上一次的估计值加以改进后;仅是向真值接近了一步而已。 二、潮流方程中的二次型及求解 直角坐标形式的潮流方程为仅包含二次项的非线性代数方程组 （2-80） 表示为 （2-81） 其中为列矢量，为节点设定量，是已知的；为潮流方程的函数矢量；为待求变量，是维矢量。 按照二次型函数的特点，将（2-81）右侧在解的初值处展开式为 （2-82） 是在的雅可比矩阵。（2-82）是（2-81）的精确表达，没任何近似，它是一个以作为变量的二次代数方程组。以平启动电压为初值，由（2-82）解出并用修正可得原方程的解。但（2-82）也是一个非线性的，其形式与原方程一致，如使用原方程（2-80）同样的解法就失去了本研究的意义。 式（2-82）可以写成 （2-83） 由此可得到一个迭代式 （2-84） 上式中k表示迭代次数；J为方程（2-80）解的初值处的雅可比矩阵，在迭代过程中保持不变，所以是一种定雅可比的算法。k=0时,如令，于是和牛顿法的第一次迭代计算完全相同。 收敛条件为： （2-85） 或 （2-86） 三、算法特点及性能估计 保留非线性快速潮流算法的特点可以通过和牛顿法进行比较而得以揭示。 设求解的方程是 ，则牛顿法的迭代公式是 （2-87） 此算法的迭代公式是 （2-88） 由两组公式可见，与牛顿法的在迭代过程中变化的雅可比矩阵不同，保留非线性快速潮流算法采用的是用初值计算而得的恒定雅可比矩阵，整个计算过程只需一次形成，可使用因子表法，所以每次迭代所需时间可以大大节省。 就每一次迭代所需的计算量而言，牛顿法要重新计算，而保留非线性快速潮流算法则要计算，由于计算函数式完全相同，仅变量不同、所以这部分的计算量是完全相同的。 两种算法的△x的含义是不同的。牛顿法的是相对于上一次迭代所得到的迭代点的修正量；而保留非线性快速潮流算法的△x(k)则是相对于始终不变的初始估计值的修正量。 保留非线性快速潮流算法达到收敛所需的迭代次数比牛顿法要多，在半对数坐标纸上收敛特性近似为一条直线。但由于每次迭代所需的计算量比牛顿法节省很多，所以总的计算速度比牛顿法可提高很多。 由于不具对称性质的雅可比矩阵经三角分解后，其上