第四复杂电力系统潮流的计算机算法剖析.ppt

基本要求：着重介绍运用电子计算机计算电力系统潮流分布的方法。它是复杂电力系统稳态和暂态运行的基础。 运用计算机计算的步骤，一般包括建立数学模型，确定解算方法，制定框图和编制程序，本章着重前两步。 §4.1 电力网络方程 电力网络方程指将网络的有关参数和变量及其相互关系归纳起来组成的，反映网络特性的数学方程式组。如节点电压方程、回路电流方程，割集电压方程。相应有： （1）节点导纳矩阵 （2）节点阻抗矩阵 （3）回路阻抗矩阵 u10=un1 u12=un1-un2 u23=un2-un3 u20=un2 u30=un3 4.1.1 节点电压方程和节点导纳矩阵的形成 几点说明 1、IB：节点注入电流的列向量，和节点注入功率类似，节点注入电流=节点电源电流和节点负荷电流之和，电源电流流入网络，其方向为正，负荷电流流出网络，其方向为负。 2、UB是节点电压的列向量，各节点电压是指各该节点与某一个被选定的参考节点之间的电压差。通常以大地为参考节点，编号为“0”。 3、YB是一个nxn阶几点导纳矩阵，其阶数n等于网络几点中除参考节点外的节点数。 二、节点导纳矩阵 按定义形成的节点导纳矩阵 因此，我们得到以下几点更为实用的结论： 节点i的自导纳Yii在数值上就等于与该节点直接相连的所有支路导纳的总和。 节点j，i之间的互导纳Yji在数值上就等于连接节点j和i之间支路导纳的负值，显然，经过简单的运算发现，Yji=Yij 另外，若节点J和i之间没有电路的连接，也不考虑两电路之间的互感，则Yji=Yij=0. 4.1.2 节点导纳矩阵的修改 修改节点导纳矩阵的原因： YB的形成与电网的结构和参数有关，它随电网的结构和参数变化而变化。（如两个节点间的支路断开，电网结构发生变化，该支路对应的互导纳变为零），因此我们要考虑结构、参数发生变化以后的节点导纳矩阵的形成。 修改节点导纳矩阵的典型方法 例题 解： 4.2 功率方程及其迭代法 1、直角坐标表示 2、极坐标表示 三、节点的编号 四、约束条件 现假设有一非线性方程组 在x（0）的附近按泰勒级数展开并忽略二阶及以上的高次项 二、计算步骤 4.3牛顿-拉夫逊潮流计算 三、修正方程 雅克比矩阵各元素 雅克比各元素具体值（证明过程略） 用直角坐标系表示时雅克比矩阵的特点 当i≠j时，Jij=-Nij ，Hij=Lij 但当i≠j时， Hij ≠Hji，Nij ≠Nji，Jij ≠Jji，Lij ≠Lji，因此雅克比矩阵并不是一个对称矩阵，这一点务必与节点导纳矩阵要区分开来。 当i≠j时，Rij=0，Sij=0，又因为电网中有很多节点之间并没有直接的连接，互导纳Yij=Gij+jBij=0，与之对应的很多非对角元素都为0，所以雅克比矩阵是一个稀疏矩阵。 雅克比矩阵的各个元素 （2）当i=j 四、牛拉潮流计算的基本步骤 五、PV节点向PQ节点的转化 4.4 P-Q分解法潮流计算 先将修正方程式4-44重新排列 简写成为： 可将上式中的N、J略去，得到。 另外，对无功注入功率（4-43b）也可以进行简化 在考虑到 说明： 以一个n-1阶和一个m-1阶系数矩阵B’、B”代替原有的n+m-2阶系数矩阵J，提高了计算速度，降低了对存储容量的要求； 修正方程的系数矩阵B’和B”为对称常数矩阵，且在迭代过程中保持不变； P-Q分解法具有线性收敛特性，与牛顿-拉夫逊法相比，当收敛到同样的精度时需要的迭代次数较多，但所用的时间少； P-Q分解法一般只适用于110KV及以上电网的计算。因为35KV及以下电压等级的线路r/x比值很大，不满足上述简化条件，可能出现迭代计算不收敛的情况。所以若不满足RX，则不能用P-Q解耦法计算潮流。 本章完毕 谢谢同学们的支持 代入： 并且等号两边同时前乘以逆矩阵，整理得出P-Q分解法修正方程： 简化为： 这两个矩阵不仅阶数不同，而且矩阵元素值也不相同。 P-Q分解法的特点： 由此可知，初值选择适当迭代收敛，选择不当将不收敛。 选择不当 选择恰当 由上节内容可知，牛顿型潮流计算的核心问题是修正方程式的建立和求解，在修正方程的建立过程中，需要对各类节点的编号做一个约定（在上一节中已讲到关于节点编号，这里在重复一下）： 一、节点的编号 二、潮流节点实际功率方程的个数 从理论上讲，假设一个n个节点的电力网络，除去一个平衡节点外，还剩下n-1各节点，所以功率方程的总数为2(n-1)个【包括有功和无功方程】。 网络中状态变量的个数也是2(n-1)个。 但现在根据上述编号约定原则，得出： ?PQ节点个数为m-1个，那么PQ节点对应的有功和无功方程共有 2(m-1)个。 ?PV节点个数为n-m个，那么PV节点对应的有功功率方程共有 n-m个。（因P