第四章-潮流的计算机算法.pptVIP

* 【例4-1】 如图所示系统中 * 用导纳表示的等值电路图 试形成系统的节点导纳矩阵。 * 解：自导纳元素为 * 互导纳为 * 最终形成的节点导纳矩阵为 * 【例4-3】 图4-7所示系统中各元件的参数与例4-1相同(为了使节点按照先PQ节点，再PV节点，最后为平衡节点的次序进行编号，以便于与公式相对照，节点的编号与例4-1不同)。取节点4为平衡节点，节点3为PV节点，节点1为联络节点，节点2为PQ节点。试用极坐标形式的牛顿-拉夫逊法计算系统的潮流分布。 平衡节点 PV节点 PQ节点 PQ节点 * 【分析】4节点系统： 2个PQ节点，分别对应两个未知量：电压幅值(或有效值)和相角； 可列写有功功率平衡方程和无功功率平衡方程两个方程式 （2*2个方程） 1个PV节点，对应一个未知量：电压相角； 极坐标形式下只需要列写一个有功功率平衡方程式 平衡节点，不需要列写方程 (1个方程) 5个方程：3个有功功率平衡方程，2个无功功率平衡方程(ΔP1、ΔP2、ΔP3、ΔQ1、ΔQ2) 5个未知量：U1、U2、θ1、θ2、θ3 * 【解】1. 计算导纳矩阵。对应于图4-7中的节点编号，导纳矩阵为 2. 给定节点电压初值， (PQ) (PQ) (PV) (Vθ) * 3. 置k=0，计算功率不平衡量 * 4．形成修正方程式。雅可比矩阵的形式为： * 经计算得修正方程式为 5. 求解修正方程，得 * 6. 进行修正，得 * 最后，用这些节点电压相位和有效值代入式(4.3.17)中的无功功率方程式便可以得出各个PV节点的注入无功功率，而代入式(4.3.18)便可以得出平衡节点的注入有功功率和注入无功功率。 于是，潮流计算问题便转化为求解非线性方程式 从而求出其中的n+m-1个未知量，即PQ节点和PV节点的电压相位和PQ节点的电压有效值。 它们和给定的平衡节点电压相位、PV节点和平衡节点给定的电压有效值一起，便是系统中全部n个节点的电压相位和电压有效值。 * 牛顿法求解功率平衡方程(多变量非线性方程) 2．修正方程式及其求解 如何应用？ * (2) f (x): 功率平衡方程式 (1) x ：待求的电压相位和有效值，写成向量形式 * 于是，便可以直接应用前面所介绍的解多变量非线性方程组的牛顿法中的方法和迭代公式来对它进行求解。 f (x): * (3)可以得出修正方程式(略去表示迭代次数的上标k)为 注意，式中采用了 来代替变量 的修正量 。以后将会看到，这样做的目的是为了使雅可比矩阵中各元素的计算式在形式上一致，以便简化雅可比矩阵元素的计算，而这并不影响计算的收敛性和计算结果的精度。 * 其中雅可比矩阵的各元素分别为： * (a) 对于式的雅可比矩阵中由虚线分开的各分块矩阵之非对角元素(j≠i)有 雅可比矩阵中的各元素： * 常数项 导数为0 当k≠j时，导数为0； 当k=j时，导数不为0 * 常数项 导数为0 当k≠j时，导数为0； 当k=j时，导数不为0 * (b) 对于各分块矩阵的对角元素(j≠i)有： * 常数项导数为0 j≠i、 j=i分开 * 常数项导数为0 j≠i、 j=i分开 * 雅可比矩阵 中的分块矩阵 、 、 和 的元素是在未知量的取值为 和 时，由上两式计算而得出的。 (4)对修正方程式进行求解后，可以得出修正量 和 ，从而可以得出 每次迭代之后,因为U,θ都有更新,所以雅各比矩阵的各个元素都需要重新计算 * 牛顿法计算潮流的迭代格式 可见，形成雅可比矩阵和求解修正方程式是牛顿法潮流计算中的主体。 * (1) 雅可比矩阵的阶数为 (2) 如果节点i和节点j之间的互导纳 ，则雅可比矩阵各子阵中的相应元素也等于零，因此雅可比矩阵也是稀疏矩阵。 (3) 雅可比矩阵中的两个对角子矩阵是不对称的，因为 修正方程式有以下特点： 但 * 因此，每次迭代要重新计算雅可比矩阵中的各个元素，从而增加了计算工作量，这是影响牛顿法潮流计算速度最重要的因素。 (4) 雅可比矩阵中各元素都是节点电压有效值和相位的函数， 因此在整个迭代过程中，所有元素都将随着节点电压相量的逐次修正而不断变化。 * 牛顿法潮流计算的收敛判据一般取 其中的ε为节点功率不平衡量(或称失配量)的容许误差，其取值范围为10-3~10-7。 迭代解出各节点的电压和相位后，便可以按式(4.3.18) 和式(4.3.17)求出平衡节点的注入有功功率