第一章电力系统潮流计算讲义.ppt

North China Electric Power University,Department of Electrical Engineering,Baoding 2008.11-2009.01,电力系统分析,第一章 电力系统潮流计算,一．概述,二．潮流计算问题的数学模型,三．潮流计算的几种基本方法,四．保留非线性潮流算法,五．最小化潮流算法,六．潮流计算中的自动调整,七．最优潮流问题,八．交直流电力系统的潮流计算,九．几种特殊性质的潮流计算问题简介,前面介绍的潮流计算可归结为求解非线性代数方程组问题，通过结合电力系统的物理特性，提出了多种求解该方程组的算法 但在实际计算中，对于一些病态系统（如重负荷系统等），往往出现计算过程振荡甚至不收敛的现象这时人们很难判定这是由于潮流算法不够完善而导致计算失败，还是从一定的初值出发，在给定的运行条件下，从数学上来讲，非线性的潮流方程组本来就是无解的（或者无实数解）五.最小化潮流算法,后来人们提出了潮流计算问题在数学上可以表示为求一个由潮流方程构成的函数（称为目标函数）最小值问题这就形成了采用数学规划的方法，称之为非线性规划潮流计算法这种方法的一个显著特点是从原理上保证了计算过程不会发散。

在给定的运行条件下，只要潮流问题有解，则目标函数最小值就迅速趋近于零；如果潮流问题无解，则目标函数先是逐渐减小，但最后却停留在一个不为零的正值上为给定条件下潮流问题的有解与无解提供了一个明确的判断途径五.最小化潮流算法,早期的应用数学规划方法的非线性规划潮流算法在内存需求量和计算速度方面都无法和前面介绍的各种潮流算法竞争，因而未得到实际推广应用以后，对非线性规划方法进行了改进，将数学规划原理和常规牛顿潮流算法有机结合起来，形成了一种新的潮流计算方法-带最优乘子的牛顿算法，通常简称为最优乘子法这种算法能有效地解决病态电力系统的潮流计算问题，并已得到广泛使用五.最小化潮流算法,一 潮流计算和非线性规划 将潮流计算问题概括为求解如下非线性代数方程组 (l-145) 构造标量函数 (l-147),五.最小化潮流算法,,,若 的解存在，则 最小值应该为零若此最小值不能为零，则说明不存在能满足原方程组 的解这样，就把原来的解代数方程组的问题转化为求 ，从而使 的问题从而将潮流计算问题转化为非线性规划问题由于没有附加的约束条件，因此在数学规划中属于无约束非线性规划的范畴五.最小化潮流算法,,,,,,按照数学规划的方法，通常由下述步骤求出 的极小点： (1)确定初始估计值 ； (2)置迭代次数 ； (3)从 出发，按照能使目标函数下降的原则，确定寻优方向 ；,五.最小化潮流算法,,,,,,(4)沿着 的方向确定使目标函数下降最多的一个点，即决定移动的步长。

由此得到新的迭代点 (l-149) 式中， 为步长因子，其数值的选择应使目标函数下降最多，即 (1-150) 可见,当 决定以后， 是 的一元函数通过求 对 的极值得到最优步长因子 五.最小化潮流算法,,,,,,,,(5) 校验 是否成立如成立，则 就是所求的解；否则，令 ，转向步骤(3)，重复循环计算五.最小化潮流算法,,,,,图1-11 求目标函数最小点示意图,五.最小化潮流算法,由上可见，为了求得问题的解，关键要解决两个问题： (1) 确定下一次迭代的搜索方向 ； (2) 确定下一次迭代的最优步长因子 确定迭代的搜索方向 常用的方法有： 梯度法(最速下降法) 、Powell方法 牛顿法、DFP算法、BFGS算法等.,五.最小化潮流算法,,确定最优步长因子 ： 确定最优步长因子用一维搜索法，通过一维搜索法确定最优步长因子常用的方法有： 解析法、插值法、黄金分割法等 早期的研究工作中，为确定 和 采用的非线性规划算法由于所需的内存量和计算速度都不能和牛顿法等常规潮流计算方法相比，因此作为一种潮流算法，没有被普遍采用五.最小化潮流算法,,非线性规划的计算过程能对收敛过程加以控制，迭代过程总是使目标函数下降，永远不发散，这些特点是牛顿法等常规潮流算法所没有的。

五.最小化潮流算法,二 带最优乘子的牛顿潮流算法 1)为了改进非线性规划潮流算法，首先在决定搜索方向 上，人们提出了利用常规牛顿潮流算法每次迭代所求出的修正量向量 (1-151) 作为搜索方向，称之为目标函数在 处的牛顿方向五.最小化潮流算法,,,,,由于牛顿法的雅可比矩阵高度稀疏并且已有了一套行之有效的求解修正方程式的方法，因此在决定 时可以充分利用原来牛顿潮流算法在内存和计算速度方面的优势五.最小化潮流算法,,,,,2)接着决定最优步长因子 已知对确定的 ，目标函数 是的一元函数 (l-152) 现在的问题是写出这个一元函数的解析表示式 如果有了这样的式子， 则 可以通过下式求得 (l-153),,五.最小化潮流算法,,,,,,,,,应用本章前面的式(1-77)，可以得到计算 的有效方法 由式(1-77)，采用直角坐标的潮流方程的泰勒展开式可以精确地表示为 (l-154) 引入标量乘子 调节变量 的修正步长，上式可写为 (l-155),五.最小化潮流算法,,,,,,为表达简明起见，分别定义三个向量 (1-156),五.最小化潮流算法,,,于是式(1-155)可简写为 (1-157) 代入式(1-146)，则目标函数可写为 (1-158),五.最小化潮流算法,,将 对 求导，并令其等于零，由此可以求得最优乘子 (1-159) 展开可得 (1-160) 其中 (1-161),五.最小化潮流算法,,,,,,可用牛顿法求解,以上介绍了从搜索方向 和最优步长因子 两个方面对非线性规划潮流算法作的改进。

不难看到经过改进的算法实质上是常规牛顿潮流算法和计算最优乘子这部分算法的结合因此对于现有的采用直角坐标的牛顿法潮流程序，只需增加计算最优乘子的部分，就可以改造成为应用非线性规划原理的算法，使得潮流计算的收敛过程能有效地得到控制五.最小化潮流算法,,,,,计算 的原理框图,,牛顿潮流法第k次迭代修正量的计算公式为： 上式左边为 ，等号右边为 ，也就是在k次迭代后， 已求得，为了求 ，只要计算 即可五.最小化潮流算法,3） 最优乘子法计算量的讨论,如果进一步推导，还可看出第（K+1）次潮流方程的偏差量（ ）可不必按 直接计算，可方便地用第k次迭代中已经求得的计算 而得即： （证明见书p36）,,,五.最小化潮流算法,这样，每次迭代，从原来要计算 简化为仅计算 进一步减少了计算量 分析可见，为了计算最优乘子而增加的计算量是很少的五.最小化潮流算法,三、 带有最优乘子的牛顿潮流算法具体应用： 可分为以下三种不同的情况讨论： (1) 从一定的初值出发，原来的潮流问题有解用带有最优乘子的牛顿潮流算法求解时，目标函数 下降为零， 经过几次迭代以后，稳定在1.0附近。

五.最小化潮流算法,,,(2) 从一定的初值出发，原来的潮流问题无解这种情况下使用这种算法求解时，目标函数开始时逐渐减小，但迭代到一定的次数以后即停滞在某一个不为零的正值上，不继续下降 的值则逐渐减小，最后趋近于零趋近于零是所给的潮流问题无解的标志，这说明 有异常变化，只是由于存在着一个趋于零的 ，才使得计算过程不致发散五.最小化潮流算法,,,(3) 有别于上两种情况，当采用这个方法计算时，不论迭代多少次， 的值始终在1.0附近摆动，但目标函数却不能降为零或不断波动 的值趋近于1.0说明了解的存在，而目标函数不能继续下降或产生波动可能是由于计算的精度不够所致，这时若改用双精度计算往往能解决问题五.最小化潮流算法,,,可见，采用带有最优乘子 的牛顿潮流算法以后，潮流计算不会发散，即从算法上保证了计算过程的收敛性，从而有效地解决了病态潮流的计算问题而通过 的具体数值，提供了在给定的运算条件下，潮流问题是否存在解的一个判断标志五.最小化潮流算法,,,前面介绍的各种潮流算法，构成了潮流程序的核心部分除此之外，一些实用的潮流程序往往还附有模拟实际系统运行控制特点的自动调整计算功能这些调整控制大都属于所谓的单一准则控制，即调整系统中单独的一个参数或变量以使系统的某一个准则得到满足。

这方面的具体例子有： (1) 自动调整有载调压变压器的分抽头以保持变压器某侧节点或某个远方节点的电压为规定的数值六．潮流计算中的自动调整,(2) 自动调整移相变压器的移相抽头以保持通过该移相变压器的有功功率为规定值 (3) 自动调整互联系统中某一个区域的一个（或数个）节点的有功出力以保持本区域和其它区域间的净交换有功功率为规定的数值六．潮流计算中的自动调整,（4）此外，节点的无功功率越界、节点的电压越界的自动处理，负荷静态特性的考虑等也属于潮流计算中自动调整的范畴六．潮流计算中的自动调整,为了在潮流计算中引入自动调整，对于单一准则控制问题，通常有两类方法： 第一类方法：按照所要保持的系统状态量 和当前的计算值 的差值大小，不断地在迭代中改变控制参数 的大小 大小的改变按照偏差反馈的原理进行，即 (1-164) 式中， 对减少迭代次数，保证收敛有很大影响六．潮流计算中的自动调整,这一类方法不改变原来的潮流计算方程，算法的迭代矩阵以及变量的组成均无变化 由于加入了调整，往往使得达到收敛所需的迭代次数和无调整的潮流计算相比有较多的增加，有的达到2-3倍六．潮流计算中的自动调整,第二类方法：则要改变原来潮流方程的构成，如增加或改写其中的一些方程式，为此待求变量的组成以及迭代矩阵（如雅可比矩阵等）的结构也有变化。

属于这一类的一些比较成功的自动算法能使达到收敛所需的迭代次数非常接近无调整的算法六．潮流计算中的自动调整,各种潮流计算方法，往往要根据算法本身的特点，以不同的方式引入自动调整本节介绍在牛顿法潮流算法中实现自动调整的有关方法六．潮流计算中的自动调整,一 节点无功功率越界和 节点电压越界的处理 发电机节点及具有可调无功电源的节点，常被指定为 节点在潮流计算过程中，它们的无功出力 可能会超出其出力限制值 （包括上界及下界）为此，潮流程序必须对 节点的无功出力加以监视并在出现越界时加以处理六．潮流计算中的自动调整,,,对于用牛顿算法的程序，当在迭代过程中发现无功功率越界时，即将这一节点转化为给定无功功率 的 节点显然，这种节点类型的改换将导致修正方程结构的变化对采用极坐标形式的修正方程将增加一个与 对应的方程式而在采用直角坐标形式时，则用与 对应的方程式代替原来与 对应的方程式六．潮流计算中的自动调整,,,,,,,,由干牛顿法每次迭代都要重新形成雅可比矩阵，因此就每一次迭代来说，采用这种节点形式转换的处理方法并不增加多少计算量在随后的迭代过程中，若出现该节点的电压又高于（对应于原来 越上界）或低于（对应于原来 越下界） 节点的规定电压值 时，则该节点在下一次迭代中应重新转换成 节点。

六．潮流计算中的自动调整,,,,,,,,,,,,节点的电压越界（包括越上界及下界）可以通过将该节点转换成 节点的办法来处理，即将该节点的电压固定在电压的上界或下界上但这种处理方式的前提是该节点必须具有足够的无功调节能力（即有可调的无功电源，包括无功补偿设备），因而不是所有的节点都可以这样处理六．潮流计算中的自动调整,,,,在迭代过程中，这种节点由 节点再复原为 节点的判据是节点的实际无功功率计算值 和原来给定的无功功率 的差出现正或负值（分别对应于原来节点电压越上界和越下界）。