潮流计算的计算机方法.docx

潮流计算的计算机方法
对于复杂网络的潮流计算，一般必须借助电子计算机进展。其计算步骤是：建立电力网 络的数学模型，确定计算方法、制定框图和编制程序。本章重点介绍前两局部，并着重阐述 在电力系统潮流实际计算中常用的、根本的方法。
1, 电力网络的数学模型
电力网络的数学模型指的是将网络有关参数相变量及其相互关系财纳起来所组成的.可 以反映网络性能的数学方程式组。也就是对电力系统的运行状态、变量和网络参数之间相互 关系的一种数学描述。电力网络的数学模型有节点电压方程和回路电流方程等，前者在电力 系统潮流计算中广泛采用。节点电压方程又分为以节点导纳矩阵表示的节点电压方程和以节 点阻抗矩阵表示的节点电压方程。
(1) 节点导纳矩阵在电路理论课中。已讲过了用节点导纳矩阵表示的节点电压方程： 对于n个节点的网络其展开为：
上式中，I是节点注入电流的列向量。在电力系统计算中，节点注入电流可理解为节点电源 电流与负荷电流之和，并规定电源向网络节点的注人电流为正。那么，只有负荷节点的注入 电流为负，而仅起联络作用的联络节点的注入电流为零。U是节点电压的列向量。网络中有 接地支路时，通常以大地作参考点，节点电压就是各节点的对地电压。井规定地节点的编号 为0。y是一个nXn阶节点导纳矩阵，其阶数n就等于网络中除参考节点外的节点数。
物理意义：节点i单位电压，其余节点接地，此时各节点向网络注入的电流就是节点i 的自导纳和其余节点的与节点i之间的互导纳。
特点：对称矩阵，稀疏矩阵，对角占优
⑵节点阻抗矩阵
对导纳阵求逆，得:
Z “8 = Ub其中
称为节点阻抗矩阵，是节点导纳矩阵的逆阵，
物理意义：节点i注入单位电流，其余节点不注入电流，此时各节点的电压就是节点i 的自阻抗和其余节点的与节点i之间的互阻抗。
特点：满阵，对称，对角占优
2, 功率方程、变量和节点分类
(1) 功率方程
的是节点的注入功率，因此，需要重新列写方程：
其展开式为：
所以：
展开写成极坐标方程的形式：
所以节点的功率方程为：
(2) 变量分类
负荷消耗的有功、无功功率取决于用户，因而是无法控制的，故称为不可控变童或扰动 变量。一般以列向量d表示，即
电源发出的有功、无功功率是可以控制的变量，故称为控制变量，以列向量u表示，即 母线或节点电压和相位角是受控制变量控制的因变量，故称为状态向童。
—般对于有n个节点的电力系统(除接地点外)，扰动变量d,控制变量u,状态变量x 皆是2n阶列向量，共有变量6n个.对于实际的电力系统仍然不好求解。于是对于实际的电 力系统作了某些符合实际的规定：出于节点负荷己知.于是给定2n个扰动变量。其次，又 给定2(n — 2)个控制变量，余下2个控制变量待定，以便平衡系统中的有功和无功功率，最 后给定2个状态变量，要求确定2(n—1)个状态变量。
由上述的规定.就确定了 4n个变量、只剩下2n个变量是待求的。这样就可以从2n个 方程式中解出2n个未知变量。但实际上，这个解还应满足一些约束条件。这些约束条件足 保证系统正常运行不可少的。
系统中的节点因给定的变量不同分为三类
(3) 节点分类
笫一类称PQ节点。对于这类节点，等值负荷功率和等值电源功率是给定的，从而注入 功率也是给定的，待求的那么是节点电压的大小。属于这一类节点的有按给定有功、无功功 率发电的发电母线和没有电源的变电所母线。
第二类称PV节点。对这类节点，等值负荷和等值电源的有功功率是给定的，从而注入 有功功率是给定的。等值负荷的无功功率和节点电压的大小是给定的。特求的那么是等值电 源的无功功率和节点电压的相位角。有一定无功储藏的发电厂和有一定无功功率电源的变电 所母线都可选作PV节点。
第三类称平衡节点。潮流计算时、一般都只设一个平衡节点。对这个节点，等值负荷功 率是给定的，节点电压的大小和相位角也是给定的，待求的那么是等值电源功率。担负调整 系统频率任务的发电厂母线往往被选作平衡节点。
进展计算时，平衡节点是不可少的，一般只有一个；PQ节点是大量的，PV节点少， 甚至可以不设。
3, 高斯——塞德尔方法
(1) 雅可比迭代法雅可比迭代法的根本思想：
以导纳矩阵为根底的潮流计算的根本方程式是：
展开为：
再改写为以节点电压为求解对象的形式：
那么雅可比迭代法求解潮流方程的迭代格式为：
收敛条件为：
4, 牛顿一拉夫逊法潮流计算
是目前求解非线性方程最好的方法，根本思想是把非线性方程的求解过程变成反复对线 性方程组的求解，通常称为逐次线性化过程。这里先从一维方程式的解来说明它的意义和推 导过程，然后再推广到n维的情况。
设有非线性方程式：
求解此方程，设&为近似值，△刈为近似值与真解的误差，那么有： 台劳展开有：
略去高次项有：
这是对于变量的修正量的线/生方程式，称修正方程式，用它可以求出修正量：
由于是修正量的近似值，故用它修正后的xl并不是方程的真解，只是向真解更逼近了 一些。
得到更逼近的解：
这种迭代继续进展下去，直至：
方程的解为：
牛顿——拉夫逊法可以推广到多变量非线性方程组的情况，设有非线性方程组：
用近似解和修正量表示如下：
求偏导数，略去高次项，写为矩阵的形式有：
缩写为：
迭代格式为：
收敛条件为：
从以上分析看出：牛顿•拉夫逊法求解非线性方程组的过程，实际上是反复求解修正方 程式的过程。因此，牛顿一拉夫逊法的收敛性比较好，但要求其初值选择得较为接近它们的 准确解、当初值选择得不当，可能出现不收敛或收敛到无实际工程意义的解的情况，这种现 象。为此，应用牛顿一拉夫逊法计算潮流分布的某些程序中，采用对初值不太敏感的高斯- 塞得尔法迭代一、二次后再转入牛顿一拉夫逊法继续迭代这样就能收到比较好的效果。
下面来看一下，如何通过牛顿一拉夫逊法求解潮流方程。潮流方程的根本形式：
△Pi =七-Pdi -UjX vj(Gjj cos丹 + B . sin 介.)尸：i=1、2、〔公式 4-85]
△0 = QGi~Qdi-U力Jj(G.. sin弓 一 B. cos^..)j=i
这样的方程一共有2n个。然而由于节点类型的不同，参加迭代求解的方程也不同。
(1) 对于PQ节点，R和Q,所以两个方程全部参加迭代，待求状态量为&和L对于PV节点，P,而Q未知，所以只有有功方程参加迭代；由于电压幅值已确 定，故待求状态量为&
(2) 对于平衡节点，R和Q.都未知，所以都不参加迭代。
假设系统中节点数为n, PV节点数为m,那么PQ节点数为n-m-1 ,参加迭代的方程为 m+2(n-m-1)个。待求的状态变量也为m+2(n-m-1)。具体方程如下：
整理得：
其中：
〔公式 4-90 和 4-91〕
求得到各待求的状态变量后，再通过节点功率方程计算得到平衡节点功率和PV节点得 无功。
解算步骤：
(1) 稔入原始数据和信息〔网络参数，负荷功率，PV节点有功和电压幅值，PQ 节点有功和无功，平衡节点电压〕
(2) 形成节点导纳矩阵
(3) 给定待求状态变量初值
(4) 迭代次数k=l
(5) 求方程的不平衡量APiK和△Q/<
(6) 判断是否收敛？ max(Z\P和左Q)