电力系统分析第二章 (4).ppt

1、第四章 电力系统潮流的计算机算法,第一节 电力网络的数学模型 第二节 等值变压器模型及其应用 第三节 节点导纳矩阵的形成和修改 第四节 功率方程和变量及节点分类 第五节 高斯-赛德尔法潮流计算 第六节 牛顿-拉夫逊法潮流计算 第七节 P-Q分解法潮流计算,第三章讨论简单电力网络的潮流分布计算，理解了与之相关的各种物理现象。对于复杂电力网络的潮流计算，一般必须借助电子计算机进行。 运用电子计算机，一般要完成以下步骤： 1、建立电力网络的数学模型 2、确定解算方法 3、制定计算流程和编制计算程序 本章将着重讨论前两项，主要阐述在电力系统潮流的实际计算中常用的、基本的方法。,第一节 电力网络的数学模

2、型,电力网络的数学模型指的是将网络有关参数及其相互关系归纳起来，组成可以反映网络性能的数学方程式组。也就是对电力系统的运行状态、变量和网络参数之间相互关系的一种数学描述。有： 节点电压方程 回路电流方程 割集电压方程等 节点电压方程又分为以节点导纳矩阵表示的节点电压方程和以节点阻抗矩阵表示的节点电压方程。,一、节点导纳矩阵的节点电压方程,在电路理论中，已经讲过了节点导纳矩阵的节点电压方程 对于n个节点的网络其展开为 上式中， 是节点注入电流的列向量。 是节点电压的列向量。 是一个nn阶节点导纳矩阵。,以网络节点导纳矩阵表示的节点电压方程在进行潮流计算时,可以减少计算机的内存,提高运算速度,因此

3、是最为常用的. 二、节点阻抗矩阵的节点电压方程 由 的两边都左乘 ，可得 ，而 ，则节点电压方程为,第二节 等值变压器模型及其应用,一、变压器为非标准变比时的修正 无论采用有名制或标么制，凡涉及多电压级网络的计算，在精确计算时都必须将网络中所有参数和变量按市价变比归算到同一电压等级。实际上，在电力系统计算中总是有些变压器的实际变比不等于变压器两侧所选电压基准值之比，也就是不等于标准变比，而且变压器的变比在运行中是可以改变的。这将使每改变一次变比都要从新计算元件参数，很不方便。下面将介绍另一种可等值地体现变压器电压变换功能的模型。 二、等值变压器模型,图一 等值双绕组变压器,（a）非标准变比时的

4、修正电路,（c）以变压器导纳表示,（b）以变压器导纳表示,由图一（a）电路可得到以下关系式：,解上联立方程得：,和,对于三绕组变压器，由于在高、中压两侧有分接头，其接入理想变压器的电路如图二所示：,（b）等值电路,（a）电路,图二 等值三绕组变压器模型,三、等值变压器模型的应用 （1）采用有名制，线路参数都未经归算，变压器参数则归在低压侧。变压器阻抗为：,相应的理想变压器的变比为：,其中， 、 分别为变压器高、低压绕组实际匝数相对应的电压。,（2）采用有名制，线路和变压器参数都已按选定的变比归算至高压侧这种情况下的线路阻抗分别为,相应的理想变压器的变比为：,其中， 、 分别为变压器高、低压侧的

5、电压。,变压器阻抗为,（3）采用标么制，线路和变压器参数都已按选定的基准电压折算为标么值。这种情况下的线路阻抗的标么值分别为,相应的理想变压器变比的标么值应取：,其中， 、 分别为折算参数时任选的变压器高、低压侧基准电压。,变压器阻抗标么值为,第三节 节点导纳矩阵的形成和修改,一、节点导纳矩阵的形成 节点导纳矩阵的计算归纳总结如下： 1、 节点导纳矩阵的阶数等于电力网络中除参考电（一般为大地）以外的节点数。 2、 节点导纳矩阵是稀疏矩阵，其各行非对角非零元素的个数等于对应节点所连的不接地支路数。 3、 节点导纳矩阵的对角元素，即各节点的自导纳等于相应节点所连支路的导纳之和，即,4、节点导纳矩阵

6、的非对角元素 等于节点 和 间支路导纳的负值，即 5、节点导纳矩阵是对称方阵，因此一般只需要求取这个矩阵的上三角或下三角部分。 6、对网络中的变压器，采用计及非标准变比时以导纳表示的等值电路，并将之接入网络中。然后按此等职电路用前述方法很方便地形成节点导纳矩阵。在实际程序中，往往直接计算变压器支路对节点导纳矩阵的影响。即当新接入非标准变比的变压器支路 、 时，对原来的节点导纳矩阵修正如下：,1）增加非零非对角元素为 2）节点 的自导纳，增加一个改变量为 3）节点 的自导纳，也增加一个改变量为,二、节点导纳矩阵的修改 在电力系统计算中，对于已知网络，其节点导纳矩阵已经形成。如果网络接线发生局部变

7、化，此时不必重新计算节点导纳矩阵。仅仅需要在原节点导纳矩阵的基础上进行必要的局部修改就可以得到所求节点导纳矩阵。下面介绍几种情况。,图三 电力网络接线变更示意图,(a),(b),(c),(d),（1）从原有网络中引出一条新的支路，图三（a）。同时增加一个新的节点。 新增加节点的对角元素为： 新增加非对角元素为： 原有节点的自导纳增量为：,（2）在原有节点 和 间增加一条支路，图三（b）。此情况下节点导纳矩阵的阶数不变。有关元素修改如下：,（3）在原有节点间切除一条阻抗为 的支路，见图三（c）这种情况下，相当于在节点 和 间增加阻抗为 的支路，此时，节点导纳矩阵的阶数不变，其元素修正如下：,（4

8、）原有网络节点 和 之间支路阻抗由 改变为 ，这种情况下，可以看作是在节点 和 间切除阻抗为 的支路，并在节点 和 间增加阻抗为 的支路，如图三（d）。此时，节点导纳矩阵的阶数不变，其元素修正如下：,（5）原有网络节点 和 之间变压器的变比由 变为 时，相当于在原网络节点 和 之间切除一变比为 的变压器支路，而又增加一个变比为 的变压器支路。其元素修正如下：,第四节 功率方程和变量及节点分类,一、功率方程 每节点的注入功率方程式为： 其中： 对于N个节点的电力网络，可以列出2N个功率方程。每个 节点具有四个变量，N个节点有4N个变量，但只有2N个关 系方程式。,二、变量的分类,1、负荷消耗的有

9、功、无功功率（ 、 ）取决于用户，因而是无法控制的，故称为不可控变量或扰动变量。一般以列向量 表示，即,2、电源发出的有功、无功功率（ 、 ）是可以控制的变量，故称为控制变量，以列向量 表示，即,3、母线或节点电压和相位角（ 、 ），是受控制变量控制的因变量。其中 主要受 的控制， 主要受 的控制。故 、 称为系统的状态变量，以列向量 表示，即,三、节点的分类 1、PQ节点：已知P、Q 负荷、过渡节点，PQ给定的 发电机节点，为大部分节点 2、 PV节点：已知P、V 给定PV的发电机节点， 具有可调电源的变电所， 为少量节点 3、 平衡节点基准节点： 也称为松弛节点，摇摆节点,PQ节点,平衡节

10、点,第五节 高斯-塞德尔法潮流计算,迭代法 考察下列形式的方程： 这种方程是隐式的，因而不能直接得出它的根，但如果给出根的某个猜测值，代入上式的右端，即可求得： 再进一步得到：,如此反复迭代： 确定数列xk有极限 则称迭代过程收敛，极限值x*为方程的根。 上述迭代法是一种逐次逼近迭代法，称为高斯迭代法。,高斯-塞德尔迭代法 在高斯法的每一次迭代过程中是用上一次迭代的全部分量来计算本次的所有分量，显然在计算第i个分量时，已经计算出来的最新分量并没有被利用，从直观上看，最新计算出来的分量可能比旧的分量要好些。因此，对这些最新计算出来的第k+1次近似分量加以利用，就是高斯-塞德尔迭代法。 高斯-塞德

11、尔迭代法计算潮流 功率方程的特点：描述电力系统功率与电压关系的方程式是一组关于电压的非线性代数方程式，不能用解析法直接求解 。,假设有n个节点的电力系统，没有PV节点，平衡节点编号为s，功率方程可写成下列复数方程式： 对每一个PQ节点都可列出一个方程式，因而有n-1个方程式。在这些方程式中，注入功率Pi和Qi都是给定的，平衡节点电压也是已知的，因而只有n-1个节点的电压为未知量，从而有可能求得唯一解。,高斯-塞德尔迭代法解潮流如下： 如系统内存在PV节点，假设节点p为PV节点，设定的节点电压为Up0。假定高斯-塞德尔迭代法已完成第k次迭代，接着要做第k+1次迭代前，先按下式求出节点p的注入无功

12、功率：,然后将其代入下式，求出节点p的电压： 在迭代过程中，按上式求得的节点p的电压大小不一定等于设定的节点电压Up0，所有在下一次的迭代中，应以设定的Up0对电压进行修正，但其相角仍保持上式所求得的值，使得 如果所求得PV节点的无功功率越限，则无功功率在限，该 PV节点转化为PQ节点。,高斯-塞德尔迭代法计算潮流的步骤： 设定各节点电压的初值，并给定迭代误差判据； 对每一个PQ节点，以前一次迭代的节点电压值代入功率迭代方程式求出新值； 对于PV节点，求出其无功功率，并判断是否越限，如越限则将PV节点转化为PQ节点； 判别各节点电压前后二次迭代值相量差的模是否小于给定误差，如不小于，则回到第2

13、步，继续进行计算，否则转到第5步； 根据功率方程求出平衡节点注入功率； 求支路功率分布和支路功率损耗。,第六节 牛顿-拉夫逊法潮流计算,一、牛顿-拉夫逊法的基本原理 设有单变量非线性方程 (4-1) 给出解的近似值 ，它与真解的误差为 ，则 可得 将上式左边的函数在 附近展成泰勒级数，便得 （4-2）,如果差值 很小， 的二次及以上阶次的各项可略去，式（4-2）便简化成 上式是修正量 的线性方程式，也称为修正方程式，解此方程可得修正量 用所求的 去修正近似解，便得 修正后的近似解 同真解仍有误差，为进一步逼近真解，这样的迭代计算反复进行下去，迭代计算通式是,（4-3） 迭代过程的收敛判剧为 （

14、4-4） 或 （4-5） 式中， 和 为预先给定的小正数。 下图为这种解法的几何意义，函数 为图中曲线。 的解相当于曲线与 轴的交点。如果第 次迭代中得到 ，则过点 点作一切线，此切线同 轴的交点便确定了下一个近似解 。由此可见，牛顿-拉夫逊法实质上就是切线法，是一种逐步线性的方法。,牛顿法的几何解释 牛顿法也适用于多变量非线性代数方程的求解。 设有 个联立的非线性代数方程 （4-6）,假定已给出各变量的初值 ，令 分别为各变量的修正量，使其满足方程（4-6），即 （4-7） 将上式的 个多元函数在初始值附近分别展成泰勒级数，并略去含有 的二次及以上阶次各项，使得 （4-8）,方程式（4-8）

15、也可以写成矩阵形式 （4-9） 方程式（4-9）是对于修正量 的线性方程组，称为牛顿法的修正方程式，利用高斯消去法或三角分解法可解出 ，然后对初始近似解进行修正 （4-10） 如此反复迭代，在进行第 次迭代时，从而求得修正方程式,（4-11） 得到修正量 ，并对各修正量进行修正 （4-12） 式（4-11）和式（4-12）也可缩写为 （4-13） 和 （4-14）,式中， 和 分别是由 个变量和修正量组成的 维列向量； 是由 个多元函数组成的 维列向量； 是 阶方阵，称为雅可比矩阵，它的第 ， 个元素 是第 个函数 对第 个变量 的偏导数；上角 标 表示 阵每一个元素都在点 处取值。 迭代过程

16、一直进行到满足收敛判剧 （4-15） 或 （4-16） 为止， 和 为预先给定的小正数。,将牛顿-拉夫逊法用于潮流计算，要求将潮流方程写成形如式（4-16）的形式。由于节点电压可以采用不同的坐标系表示，牛顿-拉夫逊法潮流计算也就相应地采取不同的计算公式。 二、节点电压用直角坐标表示时的牛顿-拉夫逊法潮流计算 采用直角坐标时，节点电压可表示为 导纳矩阵元素则可表示为,将上述表达式代入 ，展开并分出实部和虚部，便得 （4-17） 假定系统中的第 号节点为 节点，第 个节点的给定功率设为 和 ，对该节点可列写方程 （4-18）,假定系统中的第 号节点为 节点，则对其中每一个节点可列写方程,（4-19

17、）,第 号节点为平衡节点，其电压 是给定的，不参加迭代。,式（4-18）和式（4-19）总共包含了 个方程，待求的变量有 也是 个 。我们还可看到，方程（4-18）和式（4-19）已经具备了方程组（4-16）的形式，因此，不难写出如下的修正方程式,（4-20）,式中,上述方程中雅可比矩阵的各元素，可以对式（4-18）和式（4-19）求偏导获得，当 时 （4-21） 当 时,见式（4-22）,（4-22） 修正方程式（4-20）还可以写成矩阵的形式，如下,式中， 和 都是二维列向量； 是 阶方阵。 对于 节点,（4-23）,（4-24）,对于 节点,（4-25）,从表达式（4-21）（4-25）

18、得出雅可比矩阵的以下特点： （1）雅可比矩阵各元素都是节点电压的函数，它们的数值将在迭代过程中不断地改变。 （2）雅可比矩阵的子块 中的元素的表达式只用到导纳矩阵中的对应元素 ，则必有 。因此，（4-23）式中分块形式的雅可比矩阵同节点导纳矩阵一样稀疏，修正方程的求解同样可以应用稀疏矩阵的技巧。 （3）无论在式（4-20）或式（4-23）中雅可比矩阵的元素或子块都不具有对称性。 用牛顿-拉夫逊计算潮流的程序框图示于下图,输入原始数据,形成节点导纳矩阵,按公式（11-49）和（11-50）计算雅可比矩阵各元素,计算平衡节点功率及全部线路功率,输出,输电线路功率的计算公式如下 （4-28） 三、节

19、点电压用极坐标表示时的牛顿-拉夫逊潮流计算 采用极坐标时，节点电压表示为 节点功率方程将写成 （4-29） 式中 ，是 两节点电压的相角差。,方程式（4-29）把节点功率表示为节点电压的幅值和相角的函数。在有 个节点的系统中，假定第 号节点为 节点，第 号节点为 节点，第 号节点为平衡节点。 和 是给定的， 节点的电压幅值 也是给定的。因此，只剩下 个节点的电压相角 和第 个节点的电压幅值 是未知量。 对于每一个 或每一个 节点都可以列写一个有功功率不平衡量方程式,（4-30）,而对于每一个 节点还可以再列写一个无功功率不平衡量方程式 （4-31） 式（4-30）和式（4-31）一共包含了 个

20、方程式，正好同未知量数目相等，而比直角坐标形式的方程式少了 个。 对于方程式（4-30）和式（4-31）可以写出修正方程式如下 （4-32） 式中,； ； （4-33） ； 是 阶方阵 ，其元为 ； 是 阶矩阵，其元素为 ； 是 阶矩阵，其元素为 ； 是 阶矩阵，其元 素为 。,在上式中把节点不平衡功率对节点电压幅值的偏导数都乘以该节点电压，相应地把节点电压的修正量都除以该节点的电压幅值，这样，雅可比矩阵元素的表达式就具有比较整齐的形式。 对式（4-30）和（4-31）求偏导数，可以得到雅可比矩阵元素的表达式如下 当 时 （4-34）,当 时,计算的步骤和程序框图与直角坐标形式的相似。,（4-35）,第七节 P-Q分解法潮流计算,P-Q分解法是牛顿-拉夫逊法潮流计算的一种简化方法。 牛顿-拉夫逊法的缺点：牛顿-拉夫逊法的雅可比矩阵在每一次迭代过程中都有变化，需要重新形成和求解，这占据了计算的大部分时间，成为牛顿-拉夫逊法计算速度不能提高的主要原因。 P-Q分解法利用了电力系统的一些特有的运行特性，对牛顿-拉夫逊法做了简化，以改进和提