《潮流计算的计算机方法》由会员分享，可在线阅读，更多相关《潮流计算的计算机方法(11页珍藏版)》请在人人文库网上搜索。

1、精品文库 一、潮流计算的计算机方法 对于复杂网络的潮流计算，一般必须借助电子计算机进行。其计算步骤是：建立电力网 络的数学模型，确定计算方法、制定框图和编制程序。本章重点介绍前两部分，并着重阐述 在电力系统潮流实际计算中常用的、基本的方法。 1,电力网络的数学模型 电力网络的数学模型指的是将网络有关参数相变量及其相互关系归纳起来所组成的.可 以反映网络性能的数学方程式组。也就是对电力系统的运行状态、变量和网络参数之间相互 关系的一种数学描述。 电力网络的数学模型有节点电压方程和回路电流方程等，前者在电力 系统潮流计算中广泛采用。 节点电压方程又分为以节点导纳矩阵表示的节点电压方程和以节 点阻抗。

2、矩阵表示的节点电压方程。 (1) 节点导纳矩阵 在电路理论课中。已讲过了用节点导纳矩阵表示的节点电压方程： -n * a 对于n个节点的网络其展开为: 扎 f I是节点注入电流的列向量。在电力系统计算中，节点注入电流可理解为节点电源 上式中， 电流与负荷电流之和，并规定电源向网络节点的注人电流为正。那么，只有负荷节点的注入 电流为负，而仅起联络作用的联络节点的注入电流为零。U是节点电压的列向量。 网络中有 接地支路时，通常以大地作参考点， 节点电压就是各节点的对地电压。并规定地节点的编号 为0。y是一个n X n阶节点导纳矩阵，其阶数 n就等于网络中除参考节点外的节点数。 物理意义：节点i单位。

3、电压，其余节点接地，此时各节点向网络注入的电流就是节点 的自导纳和其余节点的与节点i之间的互导纳。 特点：对称矩阵，稀疏矩阵，对角占优 (2)节点阻抗矩阵 得: 对导纳阵求逆， 其中 称为节点阻抗矩阵，是节点导纳矩阵的逆阵。 物理意义：节点i注入单位电流，其余节点不注入电流，此时各节点的电压就是节点 的自阻抗和其余节点的与节点 i之间的互阻抗。 特点：满阵，对称，对角占优 2,功率方程、变量和节点分类 (1) 功率方程 已知的是节点的注入功率，因此，需要重新列写方程: 其展开式为: n ZYijU j丄 P - jQi Ui 欢迎下载 所以: P +jQi 展开写成极坐标方程的形式: n P 。

4、=UiZ U j (Gij cosgj + Bij sin 6ij) j4 n Qi =U i 2 U j (Gij sin 6 Bj cosj) j4 所以节点的功率方程为： n AR =PGi Pdi -UiS U j(Gij cos + Bij sin) j 二 n AQi =QGi Qdi U i 送 U j (Gij sin 6 - Bij cos ) (2) 变量分类 负荷消耗的有功、 。一般以列向量 j经 无功功率取决于用户， 因而是无法控制的， 故称为不可控变量或扰动 d表示，即 J =九P计屮q 33 电源发出的有功、 无功功率是可以控制的变量，故称为控制变量, 以列向量U表。

5、示，即 =尸门屜&縮血ar 母线或节点电压和相位角是受控制变量控制的因变量，故称为状态向量。 般对于有n个节点的电力系统(除接地点外)，扰动变量d,控制变量U，状态变量x 皆是2n阶列向量，共有变量6n个.对于实际的电力系统仍然不好求解。于是对于实际的电 力系统作了某些符合实际的规定：出于节点负荷己知.于是给定2n个扰动变量。其次，又 给定2(n 一 2)个控制变量，余下2个控制变量待定，以便平衡系统中的有功和无功功率，最 后给定2个状态变量，要求确定 2(n 1)个状态变量。 由上述的规定就确定了4n个变量、只剩下2n个变量是待求的。这样就可以从2n个 方程式中解出2n个未知变量。但实际上，。

6、这个解还应满足一些约束条件。这些约束条件足 保证系统正常运行不可少的。 系统中的节点因给定的变量不同分为三类 (3) 节点分类 第一类称PQ节点。对于这类节点，等值负荷功率和等值电源功率是给定的，从而注入 功率也是给定的，待求的则是节点电压的大小。属于这一类节点的有按给定有功、无功功率 发电的发电母线和没有电源的变电所母线。 第二类称PV节点。对这类节点，等值负荷和等值电源的有功功率是给定的，从而注入 有功功率是给定的。等值负荷的无功功率和节点电压的大小是给定的。待求的则是等值电源 的无功功率和节点电压的相位角。 有一定无功储备的发电厂和有一定无功功率电源的变电所 母线都可选作PV节点。 第三。

7、类称平衡节点。 潮流计算时、一般都只设一个平衡节点。对这个节点， 率是给定的，节点电压的大小和相位角也是给定的，待求的则是等值电源功率。 统频率任务的发电厂母线往往被选作平衡节点。 进行计算时，平衡节点是不可少的，一般只有一个； 甚至可以不设。 3,高斯一一塞德尔方法 (1) 雅可比迭代法 雅可比迭代法的基本思想： PQ节点是大量的, xK + = f (xK) 以导纳矩阵为基础的潮流计算的基本方程式是: I =YU 展开为: . n. IYjU j j壬 i =1.2.n 等值负荷功 担负调整系 PV节点少, iH Ui n . -=YiiUi + Z YjU j i =12.n j知式 再。

8、改写为以节点电压为求解对象的形式: 则雅可比迭代法求解潮流方程的迭代格式为: UjK Y(P 严- YjUjK) i =1.2.n YiiU ij土j# 收敛条件为: =maXu K 十-u K 4,牛顿一拉夫逊法潮流计算 是目前求解非线性方程最好的方法， 基本思想是把非线性方程的求解过程变成反复对线 性方程组的求解，通常称为逐次线性化过程。这里先从一维方程式的解来阐明它的意义和推 导过程，然后再推广到 n维的情况。 设有非线性方程式： 求解此方程，设X0为近似值，AX0为近似值与真解的误差，则有: 台劳展开有: 心、+()=+ fa)色認“j十尸3 (甞严 ! +*“ +严y今尹 略去高次项。

9、有: 这是对于变量的修正量的线性方程式，称修正方程式，用它可以求出修正量: 沈丈閃) 由于AX0是修正量的近似值，故用它修正后的 一些。 X1并不是方程的真解，只是向真解更逼近了 工-黑I舄 得到更逼近的解： 这种迭代继续进行下去，直至: I fa I 1 或！ Zu： I 方程的解为： 工=工+ 工 牛顿一一拉夫逊法可以推广到多变量非线 性方程组的情况，设有非线性方程组： 人(工1 以2/tr) = o /(旳于工2 *F工；)=0 X /m空， * * * 也0 用近似解和修正量表示如下: 3严+ 近严声严十工严、，龙严 几(卫严+ 叭山严+ 对叫山； 几 7 R L 品 1T 工2T 、。

10、吐 0去2 * 0如 打 =0 缩写为: 迭代格式为: 收敛条件为: max |0严山严, max 坯 严I , 从以上分析看出：牛顿-拉夫逊法求解非线性方程组的过程， 程式的过程。因此，牛顿一拉夫逊法的收敛性比较好，但要求其初值选择得较为接近它们的 精确解、当初值选择得不当，可能出现不收敛或收敛到无实际工程意义的解的情况， 象。为此，应用牛顿一拉夫逊法计算潮流分布的某些程序中，采用对初值不太敏感的高斯 塞得尔法迭代一、二次后，再转入牛顿一拉夫逊法继续迭代这样就能收到比较好的效果。 下面来看一下，如何通过牛顿一拉夫逊法求解潮流方程。潮流方程的基本形式： n 的=PGi -Pdi U 正 U j。

11、(Gij cosj +Bij sinj) j n 心 Qi =QGi Qdi Ui 送 U j (Gij si n 6 - Bij cos) j# 实际上是反复求解修正方 i=1、2、n (公式 这种现 4-85) 这样的方程一共有 2n个。然而由于节点类型的不同，参加迭代求解的方程也不同。 (1) 对于PQ节点，Pi和Qi已知，所以两个方程全部参加迭代，待求状态量为 Ui 对于PV节点，Pi已知而Qi未知，所以只有有功方程参加迭代；由于电压幅值 已确定，故待求状态量为 3 (3)对于平衡节点，Ps和Qs都未知，所以都不参加迭代。 假设系统中节点数为 n, PV节点数为m,则PQ节点数为n-m。

12、-1，参加迭代的方程为 m+2(n-m-1)个。 心R AQi iP2 AQ2 H11 N11 H12 N12 Il H 1p H1n J11 L11 J 12 L12 - J1p J 1n H21 N21 H22 N22 Il H2p N2p J 21 I L21 J 22 L22 J2p N2p Hp1 1 Np1 Hp2 Np2 - Il H pp Hpn LHn1 Nn1 Hn2 Nn2 - Il H np Hnn 待求的状态变量也为 m+2(n-m-1)。具体方程如下: 3 U1/U1 込 au2/u2 整理得: Hij签小厂第Ujjj 餌丄厂誥Uj (公式4-90和4-91) 其中: 求得到各待求的状态变量后, 无功。 解算步骤： (1) 再通过节点功率方程计算得到平衡节点功率和PV节点得 输入原始数据和信息(网络参数，负荷功率，PV节点有功和电压幅值，PQ节 点有功和无功，平衡节点电压) 形成节点导纳矩阵 给定待求状态变量初值 迭代次数k=1 求方程的不平衡量 APiK和AQiK 判断是否收敛？ max( AP和AQ)e?若是，转11 求雅可比矩阵各元素 解修正方程，得 aS和auk 计算节点电压新值：扌+1= +A K和UK+1=UK+ AUK k=k+1 ，转 5 计算节点功率，计算输电线路功率 结束。