电力系统潮流和稳定计算设计.docVIP

目 录 . 一、N-R法潮流分析与计算 1 1.常规的潮流算法 1 2.牛顿-拉夫逊法的基本原理 1 3.极坐标形式的N-R法潮流计算 2 4.N-R法潮流计算过程 3 二、算法分析 3 1 牛顿算法数学原理： 3 2 开式网潮流计算 4 3 开式地方电力网潮流计算 6 三、程序设计 7 1. 牛顿—拉夫逊法潮流计算程序 7 四、算法实例验证 10 1.潮流方程 11 2.修正方程 11 3.收敛条件 13 4.结果分析 14 一、N-R法潮流分析与计算 1.常规的潮流算法 1）高斯-赛德尔法(GS) 内存需求量很少，但计算时间长。 2）牛顿-拉夫逊法(N-R) 具有较高的收敛可靠性和收敛速度，但是需较好的初始值，且内存占有量大。 3）PQ快速解耦法 计算时间少，内存占用少，但是对病态潮流敏感。 2.牛顿-拉夫逊法的基本原理 1）N-R法解单变量非线性方程 非线性方程： 设方程初始解为，初解与真解的偏差为 则真解为： 若修正量 很小，在 处展开为泰勒级数： 忽略高次项及更高项： 当 ，或者 时，无限接近真解。 初值要选择比较接近真解，否则迭代过程可能不收敛。 2）N-R法解多变量非线性方程组 非线性方程组： 初始值 按泰勒级数展开，并忽略高次项： 收敛条件： 3.极坐标形式的N-R法潮流计算 对节点编号：1~m为PQ节点，m+1~n-1 为PV节点， n为平衡节点 。 对PQ节点，功率偏差方程： . . . . 对PV节点，功率偏差方程： . PV节点注入无功功率与平衡节点注入功率不参加迭代计算。 其中： 雅可比矩阵： 4.N-R法潮流计算过程 1）、形成网络节点导纳矩阵 2）、设置节点电压初值 3)、求雅可比矩阵各元素 4)、求修正方程中功率的不平衡量 5)、解修正方程，求 6)、求各节点电压修正后值 7)、检查是否收敛 8)、若不收敛，重复步骤 3 ~7，直至收敛。 9)、迭代结束，计算PV节点注入无功、平衡节点功率、线路功率及网络总损耗。 二、算法分析 牛顿—拉夫逊算法 牛顿算法数学原理： 牛顿法 （Newton Method）：解非线性方程f(x)=0的牛顿(Newton) 法，就是将非线性方程线性化的一种方法。它是解代数方程和超越方程的有效方法之一。 设有单变量非线性方程，给出解的近似值，它与真解的误差为，则将满足，即 将上式左边的函数在附近展成泰勒级数，如果差值很小，二次及以上阶次的各项均可略去得： 这是对于变量的修正量的线性方程式，成为修正方程，解此方程可得修正量 用所求得的去修正近似解，便得 修正后的近似解同真解仍然有误差。为了进一步逼近真解，可以反复进行迭代计算，迭代计算通式是 迭代过程的收敛判据为 式中，和为预先给定的小正数。 牛顿-拉夫逊法实质上就是切线法，是一种逐步线性化的方法，此法不仅用于求单变量方程，也适用于多变量非线性代数方程的有效方法。 牛顿法至少是二阶收敛的，即牛顿法在单根附近至少是二阶收敛的，在重根附近是线性收敛的。 牛顿法收敛很快，而且可求复根，缺点是对重根收敛较慢，要求函数的一阶导数存在。 开式网潮流计算 2.1简单开式电力网潮流计算 第一种已知条件：已知同一点的功率和电压。 如已知和，求和；或由和，求和。 第二种已知条件：已知不同点的功率和电压。 如已知和，求和；或由和，求和。 （1） 已知和，求和 计算步骤： 变压器阻抗上的功率损耗： 变压器阻抗上的电压损失（忽略）： 变压器导纳支路上的功率损耗： 线路末端导纳支路上的功率损耗： 线路阻抗上的功率损耗： 线路阻抗上的电压损失： 线路首端导纳上的功率损耗： 与上述过程类似，可由和推出和。 (2) 已知首端电压和末端功率，求首端和 1 迭代法 1）假定末端电压为（一般取该网络的额定电压） 2）根据、，逐级推算功率损耗和电压损耗，求得、。 3）根据和给定的，从首端向末端逐级推算，求得和 4）若和假设的相近，则计算结束；否则用求得的和给定的重复第2）、3）步，直至推算的和前一次的迭代结果很接近为止。 该方法在实践中通常认为经过一次反复就可获得足够精确的结果。 2.2简化计算法 1)? 假设全网各节点电压为额定电压。 2)? 根据和，从末端向首端逐级推算功率分布，最终求得 。 3)? 根据推得的和已知的，从首端向末端逐级推算电压分布，求得。 步骤： 2） 3）