[工学]电力系统稳态分析-第四章 复杂电力系统潮流的计算机算法宋芸-ok.pptVIP

随着计算机技术的迅速发展和普及，电子计算机已成为分析计算复杂电力系统各种运行情况的主要工具。在本章中将介绍应用计算机计算复杂电力系统潮流的数学模型和计算方法。 本章主要内容： 建立电力网络的网络方程，形成网络的节点导纳矩阵及节点阻抗矩阵。 导出复杂网络的功率方程（潮流方程）。 求解潮流方程，得出全网的潮流发布即电压、功率分布。 一、电力系统的等值模型 电力系统的等值模型实际上是系统中各元件等值模型按它们的相关关系而组成的，主要有： 1) 发电机模型：由它的端电压和输出功率来表示； 2) 负荷模型：由一个恒功率或负荷电压静态特性表示； 3) 输电线模型：是一个分布参数的电路，可用一个集中参数的∏型等值电路表示； 4) 变压器模型：通常用集中参数的г型等值电路表示。 各节点的净注入功率 节点导纳矩阵Y 和节点阻抗矩阵Z区别 实际上，由于所有的PQ节点和PV节点的注入有功功率都以给定，而网络中的总有功功率损耗是未知的，因此平衡节点的注入有功功率 必须平衡全系统的有功功率和有功损耗而不加以给定，这就是为什么称它为平衡节点的原因。 二、 高斯-塞德尔法潮流计算 高斯迭代法 如此反复迭代： 确定数列{xi}的极限 则称迭代过程收敛，极限值x*为方程的根。 上述迭代法是一种逐次逼近迭代法，称为高斯迭代法。 高斯-塞德尔迭代法计算潮流的步骤： 设定各节点电压的初值，并给定迭代误差判据； 对每一个PQ节点，以最新迭代得到的节点电压值代入功率迭代方程式求出新值； 对于PV节点，求出其无功功率，并判断是否越限，如越限则将PV节点转化为PQ节点； 判别各节点电压前后二次迭代值相量差的模是否小于给定误差，如不小于，则回到第2步，继续进行计算，否则转到第5步； 根据功率方程求出平衡节点注入功率Ps和Qs； 求支路功率分布和支路功率损耗。 由于非线性代数方程式可能有多个解，要得出不同的解需要给定不同的初值。但有时在求解过程中永远得不出所需要的精度，甚至迭代过程中得到的解离准确解越来越远，这种情况称为迭代不收敛。是因为所给定的初值离解太远。 因此要想得到某个希望的解，应该给定一个距离这一解较近的初值。因此初值的给定对于牛顿法来说十分重要。 牛顿-拉夫逊法计算潮流 节点功率方程式： 节点电压以直角坐标形式表示的牛顿-拉夫逊法潮流计算 即节点电压表示为： 功率方程可分成实部和虚部两个方程： 对功率方程求导，得到修正方程为： 不平衡量为： 修正方程中对各类节点的处理： PQ节点：每个PQ节点有两个变量 待求，都要参加联立求解； PV节点：节点电压有效值给定，它们之间的关系为： ，用这个关系式来代替该节点无功功率表达式，并改变雅可比矩阵中对应该节点相应的部分JL换成RS；见p130公式（4-37） 对功率方程求导，得到修正方程为： 3.节点电压以完全极坐标形式表示的牛顿-拉夫逊法潮流计 算方法，即节点电压和节点导纳矩阵都以极坐标形式表示。 功率方程为： 修正方程中对各类节点的处理： PQ节点：都要参加联立求解； PV节点：该类节点只有有功部分参加联立求解，而雅可比矩阵中该类节点无功部分则除去相应的行和列，但每次迭代完成需计算该节点的无功功率，以校验是否越限； 平衡节点：因其电压大小、相位均为已知，所以不需要参加联立求解，一般处理为，在雅可比矩阵中对应该节点的对角元素为一大数，其他部分为0，当迭代结束后再求该节点的有功功率和无功功率。 雅可比矩阵的特点 雅可比矩阵为一非奇异方阵。传统的，当节点电压以极坐标表示时，该矩阵为2(n-1)-m阶方阵（m为PV节点数）；当节点电压以直角坐标表示时，该矩阵为2(n-1)阶方阵。现在为了便于编程，一般为经过处理的2n阶。 矩阵元素与节点电压有关，故每次迭代时都要重新计算。 与导纳矩阵具有相似的结构，当Yij=0，Hij、Nij、Jij、Lij均为0，因此也是高度稀疏的矩阵。 具有结构对称性，但数值不对称 注意：当在计算过程中发生PV节点的无功功率越限时，PV节点要转化为PQ节点 牛顿-拉夫逊法计算电力系统潮流的基本步骤： 形成节点导纳矩阵； 给各节点电压设初值； 将节点电压初值代入，求出修正方程式的常数项向量； 将节点电压初值代入，求出雅可比矩阵元素； 求解修正方程式，求出变量的修正向量； 求出节点电压的新值； 如有PV节点，则检查该类节点的无功功率是否越限；