电力系统分析基础课件-04复杂电力系统潮流的计算机算法.ppt

* * * * * * * * * 一、潮流计算时的修正方程式 (m-1)个PQ节点＋(n-m)个PV节点，共n-1个 (m-1)个PQ节点 (n-m)个PV节点 4-3牛顿－拉夫逊迭代法潮流计算 (4-36a) (4-36b) (4-36c) 一、潮流计算时的修正方程式 用直角坐标表示的修正方程 PQ节点 PV节点 2(n－m) 2(m－1) 2(n－m) 2(m－1) (4-37) 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 相应的： (4-38) 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 用直角坐标表示的修正方程 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 非对角元素(i≠j) 雅可比矩阵元素值 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 对角元素(i=j) 雅可比矩阵元素值 4-3牛顿－拉夫逊迭代法潮流计算 * 雅可比矩阵的特点： （1）雅可比矩阵各元素均是节点电压相量的函数，在迭代过程中，各元素的值将随着节点电压相量的变化而变化。因此，在迭代过程中要不断重新计算雅可比矩阵各元素的值； （2）雅可比矩阵各非对角元素均与Yij＝Gij＋jBij有关，当Yij＝0，这些非对角元素也为0，将雅可比矩阵进行分块，每块矩阵元素均为2×2阶子阵，分块矩阵与节点导纳矩阵有相同的稀疏性结构； （3）非对称矩阵。 * 分块雅可比矩阵： 一、潮流计算时的修正方程式 以极坐标表示: 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 以极坐标表示的另一种修正方程式为 PQ节点 PV节点 2(n－m) 2(m－1) 2(n－m) 2(m－1) 4-3牛顿－拉夫逊迭代法潮流计算 (4-44) 一、潮流计算时的修正方程式 以极坐标表示: 4-3牛顿－拉夫逊迭代法潮流计算 用极坐标表示的修正方程式为 一、潮流计算时的修正方程式 4-3牛顿－拉夫逊迭代法潮流计算 极坐标法系数推导 展开式 计及 (4-47a) (4-47b) (4-48) 一、潮流计算时的修正方程式 极坐标法系数推导 (4-49a) (4-49b) 当i≠j ，对特定的j，只有特定节点的δj，从而δij= δi- δj 是变量 对特定的j，只有该特定节点的Uj是变量 一、潮流计算时的修正方程式 极坐标法系数推导 (4-49c) (4-49d) 当i=j ，由于δi是变量，从而所有δij= δi- δj 都是变量，可得 相似地，由于Ui是变量，可得 二、潮流计算基本步骤 1. 输入原始数据和信息：y、Pis、Qis、Uis、约束条件 2. 形成节点导纳矩阵YB 3. 设置各节点电压初值ei(0), fi(0) 或Ui(0), δi(0) 4. 将初始值代入(4-38)或（4-45）求不平衡量?Pi(0), ?Qi(0), ?Ui2(0) 5. 计算雅可比矩阵各元素（Hij、Lij、Nij、Jij、Rij、Sij） 6. 解修正方程(4-37) ，求? ei（k）， ? fi（k)或（4-44）求? Ui(k), ? δi(k) 4-3牛顿－拉夫逊迭代法潮流计算 7. 求节点电压新值ei(k+1) =ei(k)＋? ei(k), fi(k+1) = fi(k) ＋? fi(k)或Ui(k+1)= Ui(k)＋? Ui(k), δi(k+1) = δi(k)＋? δi(k+1) 8.判断是否收敛：Max| ? fi(k) |≤ε, Max| ? ei(k) |≤ε或Max| ? Ui(k |≤ε, Max| ? δi(k+1) |≤ε 9.重复迭代第4、5、6、7步，直到满足第8步的条件 10. 求平衡节点的功率和PV节点的Qi及各支路的功率 二、潮流计算基本步骤 4-3牛顿－拉夫逊迭代法潮流计算 二、潮流计算基本步骤 4-3牛顿－拉夫逊迭代法潮流计算 牛顿-拉夫逊法的缺点：牛顿-拉夫逊法的雅可比矩阵在每一次迭代过程中都有变化，需要重新形成和求解，这占据了计算的大部分时间，成为牛顿-拉夫逊法计算速度不能提高的主要原因。 P-Q分解法利用了电力系统的一些特有的运行特性，对牛顿-拉夫逊法做了简化，以改进和提高计算速度。 一、潮流计算时的修正方程式 (m-1)×(n-1) (m-1)×(m-1) (n-1)×(n-1) (n-1)×(m-1) 4-4 P－Q分解法潮流计算 一、潮流计算时的修正方程式 1、对修正方程式的第一步简化 高压网络中，各元件的XR，δ→P，相应的J≈0；U →Q，N ≈0。 4-4 P－Q分解法潮流计算 一、潮流计算时的修正方程式 2、对修正方程式的第二步简化 高压网络中，各元件的XR，使G