电气工程基础--第三节 复杂电力系统的潮流计算.pptxVIP

第三节复杂电力系统的潮流计算;

一、高斯塞德尔法潮流计算

描述电力系统功率与电压关系的方程式(3-15)是一组关于电压U的非线性代数方程式，不能用解析法直接求解。高斯迭代法是一种简单可行的求解方法。

先假设有n个节点的电力系统，没有PV节点，平衡节点编号为s，1≤s≤n，则式(3-15)可写成下列复数方程式：;

高斯法的基本思想是用迭代计算来求解式(3-45)，其等号右边是前一次迭代的计算值，等号左边为新值。;

式中，k为迭代的次数。在给定节点电压的初值后，对所有的PQ节点逐个进行式(3-46)的迭代计算，求得所有PQ节点的电压新值，然后以新值代入式(3-46)右边，进行下一次迭代。这样反复迭代，直至所有节点电压前一次的迭代值与后一次迭代值相量差的模小于给定的允许误差值ε后，结束迭代，即;?;;?;;

二、牛顿拉夫逊法潮流计算

1.牛顿拉夫逊算法原理

牛顿拉夫逊(Newton-Raphson)算法是求解非线性代数方程有效的迭代计算方法。在每一次的迭代过程中，非线性问题通过线性化逐步近似。下面以一个变量为x的非线性函

数求解过程为例加以说明。设一维非线性方程为

求解x，设真值为x*。;;;;?;;

展开成泰勒级数，并略去二阶以上项：;

整理成为如下的矩阵方程：;;;

2.直角坐标系下的牛顿拉夫逊算法

运用牛顿拉夫逊法计算潮流时，节点导纳矩阵的形成、平衡节点和线路功率的计算与运用高斯塞德尔法时相同，不同的只是迭代过程。迭代过程中，两种方法应用的基本方

程都是运用高斯赛德尔法时将其展开为电压方程，而运用牛顿拉夫逊法时将其展开为功率方程，即;

在采用直角坐标系下，节点电压和导纳可表示成

将式(3-68)代入式(3-67)，展开取出实部和虚部，得到;

根据节点分类，若第i个节点为PQ节点，给定功率设为Pis和Qis，则功率误差方程可列为;

若第i个节点为PV节点，Pis和Uis给定，则功率和电压的误差方程可列为;;;;;

当j≠i时，非对角线元素为;

由以上表达式可得出雅可??矩阵的特点：

(1)矩阵中的元素是节点电压的函数，在迭代过程中将随着节点电压的变化而改变。

(2)矩阵是不对称的。

(3)当导纳矩阵中的非对角线元素Yij为零时，雅克比矩阵中相应的元素也为零。矩阵是稀疏的，可以应用稀疏矩阵的求解技巧。;

3.极坐标系中的牛顿拉夫逊法潮流算法;

对一个具有n个独立节点，其中有n-m-1个PV节点的网络，式(3-76)组成的方程组共有n-1+m个方程式。采用极坐标时，方程组个数较采用直角坐标表示时少了n-m-1个。因为PV节点采用极坐标时，待求的只有电压的相位和注入的无功功率；而采用直角坐标时，待求量为电压的实数部分、虚数部分和注入的无功功率，因此采用极坐标可使未知变量少了n-m-1个，方程数也少了n-m-1个，这样建立修正方程式的矩阵形式为;

式中：;

H是(n-1)×(n-1)阶方阵，N是(n-1)×m阶矩阵，M是m×(n-1)阶矩阵，L是m×m阶矩阵。各矩阵种元素分别为;;

4.牛顿拉夫逊法的程序计算步骤及流程图

本节以直角坐标系为例对牛顿拉夫逊法的计算步骤说明如下：;

牛顿-拉夫逊法潮流计算流程如图3-30所示。利用极坐标系计算潮流的过程与此类似。;;

三、P-Q分解法潮流计算

1.原理分析

P-Q分解法是从简化牛顿拉夫逊法极坐标的形式上提出来的。它的基本思想是根据电力系统的实际运行特点(通常网络上各支路的电抗远大于电阻值)。因此，节点功率方程在用极坐标形式表示时，牛顿拉夫逊法的修正方程为式(3-77)，由该式可知系统母线电压幅值的微小变化ΔU对母线有功功率的改变ΔP影响很小。同样，母线电压相位的少许变化Δδ也对母线无功功率的变化ΔQ影响很小。;

因此，节点功率方程在采用极坐标形式表示时，修正方程可简化为

这就是把2(n-1)阶的线性方程组变成了两个n-1阶的线性方程组，将P和Q分开来进行迭代计算，因而大大减少了计算工作量。但是H、L在迭代过程中仍然在不断地变化，而且又都是不对称的矩阵。对牛顿拉夫逊法进一步简化(也是最关键的一步)，即把式(3-80)中的系数矩阵简化为在迭代过程中不变的对称矩阵。;;;;

将式(3-84)和式(3-85)代入式(3-80)中，得到;

经进一步整理得到简化后的修正方程：;

简记为

简记为;

P-Q分解法迭代公式的特点是，P-δ和Q-V迭代分别交替进行，功率偏差计算时使用最近修正过的电压值，且有功和无功偏差都用电