第四章电力系统潮流的计算机算法.pptVIP

第四章 复杂电力系统潮流的计算机算法 第一节 网络方程式 基本概念 网络方程：反映系统中电流与电压之间相互关系的数学方程； 如节点电压方程、回路电流方程、割集电压方程等，一般来讲，由于系统的等值电路中的接地支路较多，采用节点电压方程时的方程数比回路电流少，故在电力系统潮流计算中大都采用节点电压方程。 用节点导纳矩阵表示的网络方程式 1、网络方程的形成 用节点阻抗矩阵表示的网络方程式 第二节 潮流计算的节点功率方程和节点分类 直角标表示的节点功率方程 节点的分类 ——————一般节点———————— 负荷节点：该节点上没有发电机而只有负荷； 发电机节点：该节点上只有发电机而没有负荷； 联络节点：该节点上既没有发电机，也没有负荷； ———潮流计算中的节点分类—— 电力系统有n个节点，每个节点可能有4个变量，则共有4n个变量，而上述功率方程只有2n个，所以需要事先给定2n个变量的值。根据各个节点的已知量的不同，分成三类：PQ节点、PV 节点、平衡节点。 PQ节点：给定的是注入有功功率P和注入无功功率Q；待求量是节点电压有效值U和电压的相位?。 PV节点：给定的是注入有功功率P和节点电压有效值U；待求量是节点的注入无功功率和电压的相位?。 平衡节点：给定的是节点电压有效值U和电压的相位? ；待求量是节点的注入无功功率Q和注入有功功率P。也称为V- ?节点、松弛节点或电压参考节点 第三节 潮流计算的牛顿-拉夫逊法 极座标法系数推导 当i≠j ，对特定的j，只有特定节点的δj，从而δij= δi- δj 是变量 对特定的j，只有该特定节点的Uj是变量 极座标法系数推导 当i=j ，由于δi是变量，从而所有δij= δi- δj 都是变量，可得 相似地，由于Ui是变量，可得 极座标法潮流计算的基本步骤： 1. 输入原始数据和信息：y、C、Pis、Qis、Uis、约束条件 2. 形成节点导纳矩阵YB=CTyC 设置各节点电压初值Ui(0), δi(0) 将初始值代入(4-48) 求不平衡量?Pi(0), ?Qi(0) 5. 计算雅可比矩阵各元素（Hij、Lij、Nij、Jij） 6. 解修正方程(4-53) ，求? Ui(k), ? δi(k) 7. 求节点电压新值Ui(k+1) = Ui(k) -? Ui(k), δi(k+1) = δi(k) -? δi(k) 极座标法潮流计算的基本步骤： 8. 判断是否收敛：Max| ? Ui(k) |≤ε, Max| ? δi(k) |≤ε 9. 重复迭代第4、5、6、7步，直到满足第8步的条件 求平衡节点的功率和PV节点的Qi及各支路的功率 三、直角坐标形式的N-R潮流算法 1. 直角座标法： PQ节点 PU节点 直角坐标的缩写形式： 直角座标法矩阵表示 雅可比矩阵元素值 非对角元素(i≠j) 对角元素(i=j) 雅可比矩阵的特点： 各元素是各节点电压的函数 不是对称矩阵 ∵ Yij=0，∴ Hij= Nij= Jij= Lij= 0，另Rij= Sij= 0，故稀疏 直角坐标潮流计算的基本步骤： 1. 输入原始数据和信息：y、C、Pis、Qis、Uis、约束条件 2. 形成节点导纳矩阵YB=CTyC 设置各节点电压初值ei(0), fi(0) 4. 将初始值代入(4-38)求不平衡量?Pi(0), ?Qi(0), ?Ui2(0) 5. 计算雅可比矩阵各元素（Hij、Lij、Nij、Jij、Rij、Sij） 6. 解修正方程(4-37) ，求? ei(k), ? fi(k) 7. 求节点电压新值ei(k+1) =ei(k) -? ei(k), fi(k+1) = fi(k) -? fi(k) 8. 判断是否收敛：Max| ? fi(k) |≤ε, Max| ? ei(k) |≤ε 9. 重复迭代第4、5、6、7步，直到满足第8步的条件 求平衡节点的功率和PV节点的Qi及各支路的功率 返回 第四节 N-R潮流计算中的收敛性和稀疏技术 一、稀疏矩阵的存储 15 10 8 11 5 4 20 7 14 9 12 3 2 22 17 13 16 18 1、按坐标存储的方案 18 22 12 20 11 15 DIAG 6 5 4 3 2 1 顺序号L 对角元素 15 10 8 11 5 4 20 7 14 9 12 3 2 22 17 13 16 18 非对角元素 ICOL IROW OFFD 顺序号m 5 3 6 4 5 3 6 4 1 1 3 2 6 6 5 5 4 4 3 3 3 2 1 1 16 13 17 2 3 9 14 7 5 8 4 10 12 11 10 9 8 7 6 5 4 3 2 1 特点：按坐标位置存储，简单、