北京交通大学电力系统课程设计.doc

北京交通大学 电气工程学院 潮流计算机算法 班级： 电气0909班 学号： 姓名： 熊飞 指导教师： 吴俊勇 目录 一、在极坐标下的牛顿-拉夫逊法 3 二、PQ分解法 5 三、牛顿-拉夫逊法原理 8 四、测试系统数据IEEE14节点 9 五、牛顿拉夫逊法程序详解以及运算结果比较 12 六、PQ分解法程序以及运算结果比较 24 一、在极坐标下的牛顿-拉夫逊法： 采用极坐标表示时，节点功率方程可以表示为： 其中若1---m号为PQ节点，m+1---n-1号为PV节点，第n号为平衡节点。Vn与 是给定的，PV节点电压也是给定的。因此只剩下n-1个节点电压相角和m个节点电压幅值是未知变量。 实际上，对于每一个PQ节点或每一个PV节点都可以列写一个有功功率不平衡方程式： 总共包含了n-1+m个方程式。 将以上方程转化为矩阵形式： 即： 雅克比矩阵的各个元素： 当时： 当时 一般地，把节点不平衡功率对节点电压幅值的偏导数乘以该节点电压，相应地把节点电压修正量都除以该节点的电压复制，这样雅克比矩阵元素比较整齐。 二、极坐标下的PQ分解法 PQ分解法是由极坐标形式的牛顿法演化而来，以有功功率作为修正电压向量角度的依据，以无功功率作为修正电压幅值的依据，把有功功率和无功功率迭代分开进行。 中的N与H均可消去为0. 同时，引入以下近似： (1)一般情况下，线路两端电压的相角差不大(不超过10°~20°)，因此可以认为 (2)与系统各节点无功功率相对应的导纳 通常远小于该节点自导纳的虚部，即 ， 考虑到上述关系，略去相关项可将系数矩阵化简为： 得：P-Q分解法的修正方程式为 框图： PQ分解法： 极坐标牛顿拉夫逊法： 三、牛顿-拉夫逊法原理 潮流计算的结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响的分析。实际电力系统的潮流技术主要采用牛顿-拉夫逊法及其衍生的P-Q分解法。以下为牛顿-拉夫逊法为例分析潮流计算的基本原理。 牛顿-拉夫逊方法的核心是将一个非线性代数方程组的求解问题转化为线性代数方程组的迭代求解问题，先给出接的一个近似值，它与真值的误差为修正量，加原非线性方程组在近似值处线性化，得到线性化的修正方程，修正方程的系数矩阵为雅克比矩阵，反复求解修正方程，用所得到的修正量去不断修正近似解，直到迭代收敛为止。即通常所称的逐次线性化过程。 非线性代数方程组： QUOTE 在待求量x的某一个初始估计值x(0)附近，将上式展开成泰勒级数并略去二阶及以上的高阶项，得到如下的经线性化的方程组： 上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量： 将Δx(0)和x(0)相加，得到变量的第一次改进值x(1)。接着就从x(1)出发，重复上述计算过程。因此从一定的初值x(0)出发，应用牛顿法求解的迭代格式为： 上两式中：是函数对于变量x的一阶偏导数矩阵，即雅可比矩阵J，k为迭代次数。 迭代计算的收敛判据为： 由上可见，牛顿法的核心便是反复形式并求解修正方程式。应用牛顿拉夫逊法时，当初始估计值x(0)和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代4~5次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。 牛顿-拉夫逊法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有可能根本不收敛或收敛到一个无法运行的节点上。对于正常运行的系统，各节点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值(也称为平直电压)，如假定： QUOTE 在牛顿拉夫逊算法中，应用矩阵分块求逆方法对阶数较高的雅可比阵求逆计算进行改进，使阶数较高的雅可比矩阵的求逆变为阶数较低的四个子阵的求逆[3]，可以提高收敛速度。 四、计算步骤，及各部分功能的程序： 采用IEEE14节点测试系统以及数据： 网络拓扑结构： 该IEEE14系统特点： 1、5个电源，一个作为平衡节点，其余4个作为PV节点。 2、含有三个变压器支路，要考虑调压变压器对潮流计算的影响。 3、符合个数较多，分布广泛。、 4、9号支路有对地静止并联补偿电容。要考虑对地并联支路的影响。 5、该系统总共有14个节点，17条线路，与3条变压器支路。 将系统参数整理如下： %线路和变压器参数矩阵： 入端 出端 电阻 电抗 电导 电纳