《华北电力大学 电力系统分析基础 课件 第四章复杂电力系统潮流的计算机算法讲义》由会员分享，可在线阅读，更多相关《华北电力大学 电力系统分析基础 课件 第四章复杂电力系统潮流的计算机算法讲义（92页珍藏版）》请在金锄头文库上搜索。
1、电力系统分析基础 Power System Analysis Basis （四）,主讲人：栗然,第四章 复杂电力系统潮流的计算机算法,基本要求：本章着重介绍运用电子计算机计算电力系统潮流分布的方法。它是复杂电力系统稳态和暂态运行的基础。 运用计算机计算的步骤，一般包括建立数学模型，确定解算方法，制定框图和编制程序，本章着重前两步。,第四章 复杂电力系统潮流的计算机算法,2. 功率方程、节点分类及约束条件,4.1 电力网络方程,电力网络方程指将网络的有关参数和变量及其相互关系归纳起来组成的，反映网络特性的数学方程式组。如节点电压方程、回路电流方程，割集电压方程。相应有： （1）节点导纳矩阵 （2
2、）节点阻抗矩阵 （3）回路阻抗矩阵,网络元件：恒定参数 发电机：电压源或电流源 负荷：恒定阻抗,代数方程,一、节点电压方程,一、节点电压方程,注意： 零电位是不编号的,负荷用阻抗表示,以母线电压作为待求量,电压源变为电流源,以零电位作为参考，根据基尔霍夫电流定律,一、节点电压方程,I2,一、节点电压方程,其中,一、节点电压方程,n 个独立节点的网络，n 个节点方程,一、节点电压方程,n 个独立节点的网络，n 个节点方程,一、节点电压方程,n 个独立节点的网络，n 个节点方程,Y 节点导纳矩阵 Yii 节点i的自导纳 Yij 节点i、j间的互导纳,一、节点电压方程,Y 矩阵元素的物理意义：,二、
3、节点导纳矩阵,自导纳,Y 矩阵元素的物理意义 互导纳,二、节点导纳矩阵,节点导纳矩阵中自导纳的确定,二、节点导纳矩阵,节点导纳矩阵中互导纳的确定,二、节点导纳矩阵,节点导纳矩阵Y 的特点,直观易得 稀疏矩阵 对称矩阵,二、节点导纳矩阵,三、节点导纳矩阵的修改,不同的运行状态，（如不同结线方式下的运行状况、变压器的投切或变比的调整等）,改变一个支路的参数或它的投切只影响该支路两端节点的自导纳和它们之间的互导纳，因此仅需对原有的矩阵作某些修改。,Y 矩阵的修改,不同的运行状态，（如不同结线方式下的运行状况、变压器的投切或变比的调整等）,三、节点导纳矩阵的修改,Y 矩阵的修改,三、节点导纳矩阵的修改
4、,电力网,Y 增加一行一列（n1）（n1）,（1）从原网络引出一条支路增加一个节点,Y 矩阵的修改,三、节点导纳矩阵的修改,Y 阶次不变,Y 矩阵的修改,（2）在原有网络节点i、j之间增加一条支路,三、节点导纳矩阵的修改,Y 阶次不变,（3）在原有网络的节点i、j之间切除一条支路,Y 矩阵的修改,三、节点导纳矩阵的修改,Y 矩阵的修改,（4）在原有网络的节点i、j之间的导纳由yij改变为yij,三、节点导纳矩阵的修改,Y 矩阵的修改,（5）在原有网络的节点i、j之间变压器的变比由k*改变为k*,三、节点导纳矩阵的修改,Y 矩阵的修改,（5）在原有网络的节点i、j之间变压器的变比由k*改变为k*
5、,三、节点导纳矩阵的修改,42 功率方程及其迭代解法,一、功率方程和变量、节点的分类,1、功率方程,等值电源功率,等值负荷功率,（a）简单系统,42 功率方程及其迭代解法,一、功率方程和变量、节点的分类,1、功率方程,（b）简单系统的等值网络,一、功率方程和变量、节点的分类,1、功率方程,（c）注入功率和注入电流,42 功率方程及其迭代解法,一、功率方程和变量、节点的分类,1、功率方程,42 功率方程及其迭代解法,一、功率方程和变量、节点的分类,1、功率方程,42 功率方程及其迭代解法,一、功率方程和变量、节点的分类,2、变量的分类,42 功率方程及其迭代解法,一个电力系统有n个节点，每个节点
6、可能有4个变量Pi,Qi ,ei, fi或Pi,Qi ,Ui, i,而上述功率方程只有2n个，所以需要事先给定2n个变量的值。根据各个节点的已知量的不同，将节点分成三类：PQ节点、PV 节点、平衡节点。,一、功率方程和变量、节点的分类,2、变量的分类,42 功率方程及其迭代解法,(1)、PQ节点(Load Buses) 已知Pi,Qi ，求,ei, fi（ Ui, i, ），负荷节点（或发固定功率的发电机节点），数量最多。 (2)、PV节点(Voltage Control Buses) 已知Pi, Ui ，求, Qi, i, ，对电压有严格要求的节点，如电压中枢点. (3)、平衡节点（Slac
7、k Bus or Voltage Reference bus） 已知Ui ， i,，求, Pi, Qi, ，只设一个。,一、功率方程和变量、节点的分类,2、变量的分类,设置平衡节点的目的,42 功率方程及其迭代解法,在结果未出来之前，网损是未知的，至少需要一个节点的功率不能给定，用来平衡全网功率。,电压计算需要参考节点。,一、功率方程和变量、节点的分类,3、约束条件,42 功率方程及其迭代解法,实际电力系统运行要求： 电能质量约束条件：Uimin Ui Uimax 电压相角约束条件 |ij|=| i - j | ijmax, 稳定运行的一个重要条件。 有功、无功约束条件 Pimin Pi Pi
8、max Qimin Qi Qimax,二、高斯赛德尔迭代法（既可解线性，也可解非线性方程）,42 功率方程及其迭代解法,二、高斯赛德尔迭代法（既可解线性，也可解非线性方程）,42 功率方程及其迭代解法,可改写为：,二、高斯赛德尔迭代法（既可解线性，也可解非线性方程）,42 功率方程及其迭代解法,假设变量（x1, x2, .,xn）的一组初值（ ） 将初值代入迭代格式,完成第一次迭代 将第一次迭代的结果作为初值，代入迭代公式，进行第二次迭代 检查是否满足收敛条件：,二、高斯赛德尔迭代法（既可解线性，也可解非线性方程）,42 功率方程及其迭代解法,求解过程：,迭代收敛条件：,同一道题可能存在多种迭
9、代格式，有的迭代格式收敛，有的迭代式不收敛。下面讨论收敛条件： 当迭代格式为 定理 如果 则迭代格式 对任意给定的初值都收敛。,42 功率方程及其迭代解法,例 已知方程组 用高斯-塞德尔求解（0.01）。 解：（1）将方程组 改写成迭代公式： （2）设初值 ；代入上述迭代公式,直到|x(k+1)-x(k)| ,42 功率方程及其迭代解法,二、高斯赛德尔迭代法（既可解线性，也可解非线性方程）,若式中的aij对于Yij、xi对应Ui，yi对应,42 功率方程及其迭代解法,二、高斯赛德尔迭代法（既可解线性，也可解非线性方程）,此时可用迭代法求解。如设节点1为平衡节点，其余为PQ节点，则有：,42 功
10、率方程及其迭代解法,(1),二、高斯赛德尔迭代法（既可解线性，也可解非线性方程）,此时可用迭代法求解。如设节点1为平衡节点，其余为PQ节点，则有：,42 功率方程及其迭代解法,二、高斯赛德尔迭代法（既可解线性，也可解非线性方程）,此时可用迭代法求解。如设节点1为平衡节点，其余为PQ节点，则有：,计算步骤为：,42 功率方程及其迭代解法,二、高斯赛德尔迭代法（既可解线性，也可解非线性方程）,对各类节点的计算和处理 由于节点的类型不同，已知条件和求解对象不同，约束条件不同，在计算过程中的处理不同。,（1）PQ节点：按标准迭代式直接迭代；,（2）PV节点：已知的式Pp和Up，求解的是Qp，p；按标准
11、迭代式算出Up (k)， p (k)后，首先修正:,然后修正,42 功率方程及其迭代解法,(2),二、高斯赛德尔迭代法（既可解线性，也可解非线性方程）,对各类节点的计算和处理,检查无功是否越限，如越限，取限值,此时：PVPQ,42 功率方程及其迭代解法,(3),例题：用G-S计算潮流分布,解：网络的节点导纳距阵为：,设 ，代入式（1）求,修正U3为 ，再用式（2）计算：,然后开始第二次迭代：,再修正U3为：,因此，第二次迭代结束时节点2的电压为 节点3的电压相位角为3=2.940,与之对应的节点3的无功功率为Q3=0.0596.,再计算,三、牛顿拉夫逊迭代法（常用于解非线性方程）,原理：,按泰
12、勒级数展开，并略去高次项,42 功率方程及其迭代解法,三、牛顿拉夫逊迭代法（常用于解非线性方程）,原理：,42 功率方程及其迭代解法,三、牛顿拉夫逊迭代法（常用于解非线性方程）,42 功率方程及其迭代解法,初值不当不收敛,三、牛顿拉夫逊迭代法（常用于解非线性方程）,42 功率方程及其迭代解法,三、牛顿拉夫逊迭代法（常用于解非线性方程）,42 功率方程及其迭代解法,三、牛顿拉夫逊迭代法（常用于解非线性方程）,42 功率方程及其迭代解法,三、牛顿拉夫逊迭代法（常用于解非线性方程）,42 功率方程及其迭代解法,三、牛顿拉夫逊迭代法（常用于解非线性方程）,42 功率方程及其迭代解法,三、牛顿拉夫逊迭代
13、法（常用于解非线性方程）,42 功率方程及其迭代解法,非线性代数方程的牛顿法迭代格式为：,三、牛顿拉夫逊迭代法（常用于解非线性方程）,（1）将xi(0)代入，算出f，J中各元素，代入上式方程组，解出xi(0)；,（2）修正xi(1) xi(0) xi(0) ，算出f，J中各元素，代入上式方程组，解出 xi(1) ；,计算步骤：,注意：xi的初值要选得接近其精确值，否则将不迭代。,42 功率方程及其迭代解法,4-3牛顿拉夫逊迭代法潮流计算,一、潮流计算时的修正方程式,节点电压用直角坐标表示：,4-3牛顿拉夫逊迭代法潮流计算,一、潮流计算时的修正方程式,首先对网络中各节点作如下约定： （1）网络中
14、共有n个节点，编号为1，2，3，n； （2）网络中（m1）个PQ节点，一个平衡节点，编号为1，2，m，其中1sm为平衡节点； （3）nm个PV节点，编号为m+1,m+2,，n.,一、潮流计算时的修正方程式,(m-1)个PQ节点(n-m)个PV节点，共n-1个,(m-1)个PQ节点,(n-m)个PV节点,4-3牛顿拉夫逊迭代法潮流计算,(4-36a),(4-36b),(4-36c),一、潮流计算时的修正方程式,用直角坐标表示的修正方程,PQ节点,PV节点,(4-37),4-3牛顿拉夫逊迭代法潮流计算,一、潮流计算时的修正方程式,相应的：,(4-38),4-3牛顿拉夫逊迭代法潮流计算,一、潮流计算
15、时的修正方程式,用直角坐标表示的修正方程,4-3牛顿拉夫逊迭代法潮流计算,一、潮流计算时的修正方程式,非对角元素(ij),雅可比矩阵元素值,4-3牛顿拉夫逊迭代法潮流计算,一、潮流计算时的修正方程式,对角元素(i=j),雅可比矩阵元素值,4-3牛顿拉夫逊迭代法潮流计算,70,雅可比矩阵的特点：,（1）雅可比矩阵各元素均是节点电压相量的函数，在迭代过程中，各元素的值将随着节点电压相量的变化而变化。因此，在迭代过程中要不断重新计算雅可比矩阵各元素的值； （2）雅可比矩阵各非对角元素均与YijGijjBij有关，当Yij0，这些非对角元素也为0，将雅可比矩阵进行分块，每块矩阵元素均为22阶子阵，分块矩阵与节点导纳矩阵有相同的稀疏性结构； （3）非对称矩阵。,71,分块雅可比矩阵：,一、潮流计算时的修正方程式,以极坐标表示:,4-3牛顿拉夫逊迭代法潮流计算,一、潮流计算时的修正方程式,以极坐标表示的另一种修正方程式为,PQ节点,PV节点,4-3牛顿拉夫逊迭代法潮流计算,(4-44),一、潮流计算时的修正方程式,以极坐标表示:,4-3牛顿拉夫逊迭代法潮流计算,用极坐标表示的修正方程式为,一、潮流计算时的修正方程式,4-3牛顿拉夫逊迭代法潮流计算,极坐标法系数推导,展开式,计及,(4-47a),(4-47b),(4-48),一、潮流计算时的修正方程
