直角坐标系下牛顿法潮流计算.docxVIP

直角坐标系下牛顿法潮流计算 电力系统潮流计算 潮流计算是电力系统分析中的一种最基本的计算， 它的任务是对给定的运行 条件确定系统的运行状态，如母线上的电压 （幅值及相角 ）、网络中的功率分布及 功率损耗等。 在电力系统规划设计和现有电力系统运行方式的研究中， 都需要利 用潮流计算来定量地分析比较供电方案或运行方式的合理性 .可靠性和经济性。 此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。 节点导纳矩阵的形成 在图 1（a）的简单电力系统中，若略去变压器的励磁功率和线路电容，负 荷用阻抗表示， 便可以得到一个有 5个节点（包括零电位点） 和 7 条支路的等值 网络，如图 1（b）所示。将接于节点 1 和 4 的电势源和阻抗的串联组合变换成 等值的电流源和导纳的并联组合， 变得到图 1（c）的等值网络， 其中 I1 y10 E1和 I4 y40E4分别称为节点 1和 4 的注入电流源。 3(a)?1y24?11 y12y′102 y23 3 (a) ?1 y24 ?1 1 y12 y′10 2 y23 y34 2 y23 3 y20 y40 ?4 y12 y20 y101 (b) 4 y40 y34 (c) 图 1 电力系统及其网络 以零电位点作为计算节点电压的参考点，根据基尔霍夫定律，可以写出 4 个独立节点的电流平衡方程如下： y10U1 y12 (U1 U2) I1 y12 (U2 U1) y20U2 y23 (U2 U3) y24 (U2 U4 ) 0 y23(U3 U2 ) y34(U3 U4) 0 2-1)y24(U4 U2 ) y34(U4 U3) y40U4 I 2-1) y34 y34 上述方程组经过整理可以写成 Y11U1 Y12U2 I1 Y21U1 Y22U2 Y23U3 Y24U4 0 Y32U2 Y33U3 Y34U4 0 Y42U2 Y43U3 Y44U4 I I4 （2-2） 式 中 ， Y11 y10 y12 ； Y22 y20 y23 y24 y12 ； Y33 y23 y34 ； Y44 y 40 y24 y34 ； Y12 Y21 y12 ； Y23 Y32 y 23 ； Y24 Y42 y24 般的，对于有 n 个独立节点的网络，可以列写 n个节点方程 Y11U1 Y12U2 L Y1nUn I1 Y21U1 Y22U2 L Y2nUn I2 M Yn1U1 Yn2U2 L YnnUn In （2-3） 也可以用矩阵写成 Y11 Y12L Y1n U1 I1 Y21 Y22L Y2n U2 I2 M M M Yn1 Yn2 L Ynn Un In n n （2-4） 或缩写为 YU I （2-5） 矩阵 Y 称为节点导纳矩阵。它的对角线元素 Yii称为节点 i的自导纳，其值等 于接于节点 i 的所有支路导纳之和。 非对角线元素 Yij 称为节点 i 、 j 间的互导纳， 它等于直接接于节点 i 、j间的支路导纳的负值。 若节点i 、j间不存在直接支路，则有 Yij 0 。由此可知节点导纳矩阵是一个稀疏的对称矩阵。 牛顿 -拉夫逊法潮流计算 牛顿 -拉夫逊法的基本原理 牛顿—拉夫逊法( Newton —Raphson 法)是求解非线性方程代数方程组的 有效迭代计算方法。 在牛顿—拉夫逊法的每一次迭代过程中， 对非线性方程通过 线性化处理逐步近似。下面以单变量加以说明。 设有单变量非线性方程 (0) (0) (0) (0) f( (0) (0) (0) (0) f(x x ) f (x ) f (x ) 2 x(0) f (x(0))( x2(!0)) (n) (0) f (x ) n ( x(0))nn! f(x) 0 (3-1) (0) (0) 求解此方程时。先给出解的近似值 x 它与真解的误差为 x ，则 (0) (0) x x x 将满足方程，即 (0) (0) f(x x(0) ) 0 (3-2) (0) 将(3-8)式左边的函数在 x 附近展成泰勒级数，于是便得 3-3) x(0) 处的一阶导数，? .，n式中 f (x(0) ) ,?? f ( x(0) 处的一阶导数，? .，n 阶导数。 如果差值 x(0) 很小， 3-9式右端 x(0) 的二次及以上阶次的各项均可略去。 于是， 3-9 便简化为 f(x(0)(0)(0) (0) f(x(0) (0) (0) (0) ＝0 (0) (0) (0) x(0)) f (x(0)) f (x(0)) x (3-4) (0) 这是对于变量的修正量 x 的现行方程式，亦称修正方程式。解此方程可得 修正量 (0)xff( (0) x ff((xx(0(0)))) (3-5) 用所求的 x(0) 去修正近似解，变