《电力系统分析第11章课件》由会员分享，可在线阅读，更多相关《电力系统分析第11章课件（76页珍藏版）》请在金锄头文库上搜索。
1、第三节复杂电力系统潮流的计算机算法,TCSC、STATCOM、 UPFC、 FACTS,本章主要内容：,2. 功率方程、节点分类及约束条件,第一节 电力网的数学模型,一、节点电压方程,二、导纳矩阵的形成,三、导纳矩阵的修改,YB节点导纳矩阵,第一节 电力网的数学模型,1、节点电压方程（示例）, 参考节点的选取接地点,第一节 电力网的数学模型,注：Y距阵的维数（n-1）,第一节 电力网的数学模型,2、导纳矩阵的形成,自导纳,互导纳,节点导纳距阵的特点： 1、阶数 2、对称性 3、稀疏性,第一节 电力网的数学模型,3、导纳矩阵的修改,增加一节点,增加一条支路,第一节 电力网的数学模型,切除一条支路
2、,修改一条支路的导纳值（ yij 改变为yij ）,第一节 电力网的数学模型,修改一条支路的变压器变比值（ k*改变为k* ）,第二节 功率方程、节点分类及约束条件,一、功率方程,n个节点，2n个方程，6n个变量,n个节点，2n个方程，4n个变量:,二、节点分类,一将节点分成三类：PQ节点、PV 节点、平衡节点。 1、PQ节点(Load Buses) 已知Pi,Qi ，求Vi, i,负荷节点（或发固定功率的发电机节点），数量最多。 2、PV节点(Voltage Control Buses) 已知Pi, Vi ，求,Qi, i, ，对电压有严格要求的节点，如电压中枢点。具有无功储备的发电厂，可调
3、无功设备的变电所，数目较少,3、平衡节点 已知Vi ， i,，求, Pi, Qi, ，只设一个。 设置平衡节点的目的 在结果未出来之前，网损是未知的，至少需要一个节点的功率不能给定，用来平衡全网功率。 电压计算需要参考节点。,三、约束条件,实际电力系统运行要求： 电能质量约束条件：Uimin Ui Uimax 电压相角约束条件 |ij|=| i - j | ijmax, 稳定运行的一个重要条件。 有功、无功约束条件 Pimin Pi Pimax Qimin Qi Qimax,第四节 牛顿拉夫逊法潮流计算,一、N-R原理,1. 非线性方程的求解：,f(x)=0 设：x(0)为方程的初始近似解，x
4、(0)为与真实解的偏差 则：x= x(0) +x(0) f(x(0) +x(0)=0 按Taylors展开 f(x(0) +x(0)= f(x(0)+f(x(0)x(0)+.+ fn(x(0) (x(0)n/n!+.=0,第四节 牛顿拉夫逊法潮流计算,由于x(0）较小，故忽略高次项后： f(x(0) +x(0)= f(x(0) +f(x(0)x(0) =0 x(0) = -f(x(0)/ f(x(0) x(1) = x(0) +x(0) = x(0) - f(x(0)/ f(x(0) k次迭代时修正方程为： f(x(k)+ f(x(k)x(k) =0 x(k) = -f(x(k)/ f(x(k
5、) x(k+1) = x(k) - f(x(k)/ f(x(k),第四节 牛顿拉夫逊法潮流计算,结束迭代的条件（收敛）：|f(x(k)|1 或 |x(k)|2,物理意义,初值不当不收敛,x(k+1) = x(k) - f(x(k)/ f(x(k),第四节 牛顿拉夫逊法潮流计算,2. 非线性方程组的求解：,推广于多变量非线性方程组,第四节 牛顿拉夫逊法潮流计算,(4-32),式中 为函数 fi(x1,x2,.xn) 对自变量 xj 的偏导在初始值处的值,第四节 牛顿拉夫逊法潮流计算,(4-33),用矩阵表示,得到新的近似解：,第四节 牛顿拉夫逊法潮流计算,(4-33a),更一般的表示,第k+1次
6、迭代后的解为：,(4-33b),第四节 牛顿拉夫逊法潮流计算,可简写为： F(x(k)=J(k) x(k) (4-34) J(k) 为nn阶雅可比矩阵,其元素 为函数 fi(x1,x2,.xn) 对自变量 xj 的偏导在点（ x(k) ）的值 式(4-33b)可简写为： x(k+1)= x(k) x(k) (4-34b) 第k次迭代后用下面的公式检查是否收敛,(4-35a),(4-35b),第三节 高斯塞德尔迭代法潮流计算,两种常见的求解非线性方程的方法：,高斯-塞德尔迭代法 牛顿-拉夫逊迭代法,牛顿拉夫逊法潮流计算,潮流计算时的修正方程式（直角座标法）,(4-38a),(4-38b),1.
7、直角座标法：,PQ节点,PU节点,第四节 牛顿拉夫逊法潮流计算,二、潮流计算时的修正方程式（直角座标法）,(4-38a),(4-38b),1. 直角座标法：,PQ节点,PU节点,第四节 牛顿拉夫逊法潮流计算,直角坐标的缩写形式：,第四节 牛顿拉夫逊法潮流计算,直角座标法矩阵表示,(4-37),第四节 牛顿拉夫逊法潮流计算,雅可比矩阵元素值,非对角元素(ij),对角元素(i=j),(4-41a),(4-41b),第四节 牛顿拉夫逊法潮流计算,矩阵的特点及计算步骤,雅可比矩阵的特点：,各元素是各节点电压的函数 不是对称矩阵 Yij=0， Hij= Nij= Jij= Lij= 0，另Rij= Si
8、j= 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|
9、fi(k) |, Max| ei(k) | 9. 重复迭代第4、5、6、7步，直到满足第8步的条件 求平衡节点的功率和PV节点的Qi及各支路的功率,潮流计算时的修正方程式（框图）,见书上P165,第四节 牛顿拉夫逊法潮流计算,2. 极座标法：,第四节 牛顿拉夫逊法潮流计算,极座标法矩阵表示,(4-44),第四节 牛顿拉夫逊法潮流计算,极座标法系数推导,展开式,计及,(4-45a),(4-45b),(4-48),第四节 牛顿拉夫逊法潮流计算,极座标法系数推导,(4-49a),(4-49b),当ij ，对特定的j，只有特定节点的j，从而ij= i- j 是变量,对特定的j，只有该特定节点的Uj是变
10、量,第四节 牛顿拉夫逊法潮流计算,极座标法系数推导,(4-49c),(4-49d),当i=j ，由于i是变量，从而所有ij= i- j 都是变量，可得,相似地，由于Ui是变量，可得,第四节 牛顿拉夫逊法潮流计算,极座标法潮流计算的基本步骤：,1. 输入原始数据和信息：y、C、Pis、Qis、Uis、约束条件 2. 形成节点导纳矩阵YB=CTyC 设置各节点电压初值Ui(0), i(0) 将初始值代入(4-45) 求不平衡量Pi(0), Qi(0) 5. 计算雅可比矩阵各元素（Hij、Lij、Nij、Jij） 6. 解修正方程(4-44) ，求 Ui(k), i(k) 7. 求节点电压新值Ui(
11、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及各支路的功率,第三节 高斯塞德尔迭代法潮流计算,一、功率方程的非线性,非线性方程组，不能用常规代数求解方程方法求解 。,直角坐标形式：,极坐标形式：,第三节 高斯塞德尔迭代法潮流计算,二、高斯-塞德尔迭代法原理及求解步骤,例6-1 已知方程组 用高斯-塞德尔求解（0.01）。 解：（1）将方
12、程组 改写成迭代公式： （2）设初值 ；代入上述迭代公式,直到|x(k+1)-x(k)| ,第三节 高斯塞德尔迭代法潮流计算,二、高斯-塞德尔迭代法原理及求解步骤,设有非线性方程组的一般形式：,将其改写成下述便于迭代的形式：,第三节 高斯塞德尔迭代法潮流计算,二、高斯-塞德尔迭代法原理及求解步骤,假设变量（x1, x2, .,xn）的一组初值（ ） 将初值代入迭代格式（6-18），完成第一次迭代 将第一次迭代的结果作为初值，代入迭代公式，进行第二次迭代 检查是否满足收敛条件：,第三节 高斯塞德尔迭代法潮流计算,二、高斯-塞德尔迭代法原理及求解步骤,迭代公式：,更一般的形式：,第三节 高斯塞德尔
13、迭代法潮流计算,二、高斯-塞德尔迭代法原理及求解步骤,简化形式：,迭代收敛条件：,第三节 高斯塞德尔迭代法潮流计算,二、高斯-塞德尔迭代法原理及求解步骤,同一道题可能存在多种迭代格式，有的迭代格式收敛，有的迭代式不收敛。下面讨论收敛条件： 当迭代格式为 定理 如果 则迭代格式 对任意给定的初值都收敛。,第三节 高斯塞德尔迭代法潮流计算,三、高斯-塞德尔迭代法潮流计算,用高斯-塞德尔法计算电力系统潮流首先要将功率方程改写成能收敛的迭代形式 Q : 设系统有n个节点，其中 m个是PQ节点，n-(m+1)个是PV节点，一个平衡节点，且假设节点1为平衡节点（电压参考节点） 功率方程改写成：,1. 方程
14、表示：,第三节 高斯塞德尔迭代法潮流计算,三、高斯-塞德尔迭代法潮流计算,或更具体的形式为：,第三节 高斯塞德尔迭代法潮流计算,三、高斯-塞德尔迭代法潮流计算,上述迭代公式假设n 个节点全部为PQ节点。 式中等号右边采用第k次迭代结果，当ji时，采用第k次迭代结果。 用G-S迭代法求解的步骤： 第一步：形成节点导纳距阵； 第二步：设除平衡节点外的其它节点的初值，一般都设 ； 第三步：迭代求解，判断收敛与否？若满足收敛条件， 则迭代停止,2. 求解的步骤：,第三节 高斯塞德尔迭代法潮流计算,三、高斯-塞德尔迭代法潮流计算,3. PV节点的处理：,由于该类节点的V已知，Q未知，故在给定初值时，对该
15、类节点增加初值 ； 增加计算无功的迭代公式： 对于PV节点的计算步骤： 除了完成（6-24）的迭代计算外，还要执行（6-25）的迭代计算 对（6-25）得到的结果要进行下列三种情况的校核：,第三节 高斯塞德尔迭代法潮流计算,三、高斯-塞德尔迭代法潮流计算,(a) 这种情况由于计算得到的结果比允许的最小值还小，所以不允许以计算得到的结果再代入进行迭代，以 作为PV节点的无功功率，此时，PV节点就转化为PQ节点 (b) 这种情况由于计算得到的结果比允许的最大值还大，所以不允许以计算得到的结果再代入进行迭代，而是以 作为PV节点的无功功率，此时，PV节点就转化为PQ节点。 (c) 因求出的无功功率满足要求，所以迭代得到的结果继续代入公式（6-25）进行计算,第三节 高斯塞德尔迭代法潮流计算,三、高斯-塞德尔迭代法潮流计算,对于PV节点，由于它的U值是给定的，每次用公式（6-24）得到的结果 中的 一般不等于给定的值，这种情况要用给定的U代替计算得到的幅值，用 组成新的电压初值。如果通过迭代得到的与限值比较已经越限，则转化为PQ节点后，就不必做电压幅值的更换了。,第三节 高斯塞德尔迭代法潮流计算,三、高斯-塞德尔迭代法潮流计算,4. 潮流计算：,平衡节点的功率：,支路功率：,支路功率损耗：,第三节 高斯塞德尔
