最小化潮流算法知识分享

一、潮流计算方法之间的区别联系高斯-赛德尔法：原理简单，导纳矩阵对称且高度稀疏，占用内存小。

收敛速度很慢，迭代次数随节点数直接上升，计算量急剧增加，不适用大规模系统。

牛顿-拉夫逊法：收敛速度快，迭代次数和网络规模基本无关。

相对高斯-赛德尔法，内存量和每次迭代所需时间较多，其可靠的收敛还取决于一个良好的启动初值。

PQ 分解法（快速解耦法）：PQ 分解法实际上是在极坐标形式的牛顿法的基础上，在交流高压电网中，输电线路等元件的R<<X ，即有功功率主要取决于电压相角，而无功功率主要取决于电压幅值，根据这种特性对方程组进行简化，从而实现了有功和无功的解耦。

两大条件：(1)线路两端的相角相差不大(小于10°～20°)，而且||||ij ij G B ≤，于是可以认为：cos 1;sin ij ij ij ij G B θθ≈≤； (2)与节点无功功率相对应的导纳2/i i Q U 通常远小于节点的自导纳ii B ，也即2i i ii Q U B <<。

1. PQ 分解法用一个1n -阶和一个1n m --阶的方程组代替牛顿法中22n m --阶方程组，显著减少了内存需量和计算量。

2. 计算过程中B '、B ''保持不变，不同于牛顿法每次迭代都要重新形成雅可比矩阵，因此显著提高了计算速度。

3.雅可比矩阵J 不对称，而B '、B ''都是对称的，使求逆等运算量和所需的存储容量都大为减少。

4. PQ 分解法的迭代次数要比牛顿法多，但是每次迭代所需时间比牛顿法少，所以总的计算速度仍是PQ 分解法快。

在低压配电网中PQ 分解法不适用。

交流高压电网的输电线路的元件满足R<<X ，PQ 分解法正是基于此条件简化而来；而低电压配电网络一般R/X 比值很大，大R/X 比值病态问题也正是PQ 分解法应用中的一个最大障碍。

第四章潮流计算的计算机算法第一节概述潮流计算是电力系统最基本、最常用的计算。

根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值及相角），各元件中流过的功率、整个系统的功率损耗等。

潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。

因此潮流计算在电力系统的规划设计、生产运行、调度管理及科学研究中都有着广泛的应用。

电力系统潮流计算分为离线潮流计算和在线潮流计算。

前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。

本章主要讨论离线潮流计算问题，它的基本算法同样适用于在线潮流计算。

潮流计算在数学上是多元非线性方程组的求解问题，求解的方法有很多种。

自从五十年代计算机应用于电力系统以来，当时求解潮流的方法是以节点导纳矩阵为基础的逐次代入法（导纳法），后来为解决导纳法的收敛性较差的问题，出现了以阻抗矩阵为基础的逐次代入法（阻抗法）。

到六十年代，针对阻抗法占用计算机内存大的问题又出现了分块阻抗法及牛顿－拉夫逊（Newton－Raphson）法。

Newton —Raphson法是数学上解非线形方程式的有效方法，有较好的收敛性。

将N－R法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧，使N－R法在收敛性、占用内存、计算速度方面的优点都超过了阻抗法，成为六十年代末期以后普遍采用的方法。

同时国内外广泛研究了诸如非线形规划法、直流法、交流法等各种不同的潮流计算方法。

七十年代以来，又涌现出了更新的潮流计算方法。

其中有1974年由B、Stott、O、Alsac 提出的快速分解法以及1978年由岩本伸一等提出的保留非线性的高129速潮流计算法。

其中快速分解法（Fast decoupled load flow）从1975年开始已在国内使用，并习惯称之为PQ分解法。

由于PQ分解法在计算速度上大大超过N－R法，不但能应用于离线潮流计算，而且也能应用于在线潮流计算。

潮流计算的基本算法及使用方法Company number：【0089WT-8898YT-W8CCB-BUUT-202108】潮流计算的基本算法及使用方法一、 潮流计算的基本算法1.牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。

这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。

牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。

因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。

而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。

1．2 一般概念对于非线性代数方程组即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。

由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。

接着再从()1x 出发，重复上述计算过程。

因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。

潮流计算的概念和基本原理一、 潮流计算的意义电力系统潮流的计算和分析是电力系统运行和规划工作的基础。

运行中的电力系统，通过潮流计算可以预知，随着各种电源和负荷的变化以及网络结构的改变,网络所有母线的电压是否能保持在允许范围内,各种元件是否会出现过负荷而危及系统的安全,从而进一步研究和制订相应的安全措施.规划中的电力系统,通过潮流计算,可以检验所提出的网络规划方案能否满足各种运行方式的要求，以便制定出既满足未来供电负荷增长的需求,又保证安全稳定运行的网络规划方案。

二、 潮流计算的基本概念潮流计算的一般提法是:已知电力网络的结构和参数，已知各负荷点、电源点吸取或发出的有功功率和无功功率(ＰＱ节点),给定电压控制点的电压幅值和有功功率(P Ｖ节点),对指定的一个平衡节点给定其电压幅值和相位角（V θ点）,求解全网各节点电压幅值和相位角，并进一步算出各支路的功率分布和网络损耗．求解潮流问题的基本方程式是节点功率平衡方程。

三、 潮流计算的基本原理1. 潮流计算的基本模型1．1潮流方程电力系统是由发电机、变压器、输电线路及负荷等组成,其中发电机及负荷是非线性元件,但在进行潮流计算时,一般可以用接在相应节点上的一个电流注入量来代表。

因此潮流计算所用的电力网络系由变压器、输电线路、电容器、电抗器等静止线性元件所构成，并用集中参数表示的串联或并联等值支路来模拟。

结合电力系统的特点，对这样的线性网络进行分析,普通采用的是节点法，节点电压与节点电流之间的关系V Y I= ﻩ (1－1)ﻩ其展开式为 j n j ij i V Y I ∑==1 ),,3,2,1(n i =(1－2)ﻩ在工程实际中，已经的节点注入量往往不是节点电流而是节点功率，为此必须应用联系节点电流和节点功率的关系式ﻩﻩi ii i V jQ P I *-= ),,3,2,1(n i = （1－３)将式(１-3)代入式（１－２)得到j n j ij iii V Y V jQ P ∑=*=-1 ),,3,2,1(n i =ﻩ (1－4) 交流电力系统中的复数电压变量可以用两种极坐标来表示ﻩ i j i i e V V θ= （1－5）或ﻩﻩﻩ ii i jf e V += ﻩ （1－６） 而复数导纳为ij ij ij jB G Y += （1-7)ﻩ将式（1－6)、式(1－７)代入以导纳矩阵为基础的式（１－4),并将实部与虚部分开，可以得到以下两种形式的潮流方程。

潮流计算的基本算法及使用方法一、欧阳家百（2021.03.07）二、潮流计算的基本算法1.牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。

这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。

牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。

因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。

而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。

1．2一般概念对于非线性代数方程组即()0,,,21=nixxxf ()ni,2,1= (1－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。

由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。

接着再从()1x 出发，重复上述计算过程。

因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。

由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。

配电网的前推回推法潮流计算的计算机算法一、用途潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。

对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求：对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危机系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件（线路、变压器等）是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。

潮流计算是电力系统分析最基本的计算。

除他自身的重要作用之外，潮流计算还是网损计算、静态安全分析、暂态稳定计算、小干扰静态稳定计算、短路计算、静态和动态等值计算的基础。

二、技术特点在辐射状配电网的潮流计算中，采用基于支路电流或支路功率的前推回推法进行潮流计算，其中在寻找节点的计算顺序时采用关于节点邻接表的广度优先搜索方法，或者其他更加简洁的方法。

1、原理配电网潮流有很多算法，主要包括牛顿法、快速解耦法、Z法、bus 回路阻抗法和前推回推法，其中前推回推法具有更大的优势，更适于求解配电网的潮流计算。

它具有编程简单、数值稳定性好、计算效率高等优点。

辐射型配电网的接线方式可以分为辐射式、链式、干线式三种网络。

辐射型配电网潮流计算有如下特点:（1）辐射型配电网支路数一定小于节点数。

因此，网络节点导纳矩阵稀疏度很高。

（2）低压配电网由于线路电阻较大，一般不满足R<<X，因此在配电网中采用P-Q解耦法进行网络潮流计算难以收敛；由于配电网络直接面向用户，所以网络节点众多，如采用传统的潮流算法（牛顿拉夫逊法、快速解耦法）会导致导纳矩阵非常庞大，处理的工作量较大，占用的资源也较多。

（3）对于末端负荷节点前的支路电流就只是由于末端运算负荷功率产生的，所以可以直接用于末端支路的电流，依次往前推。

因此，可以采用前推回推法进行网络潮流计算。

前推回推法有基于支路电流和基于支路功率两种形式。

潮流计算的基本算法及使用方法一、欧阳歌谷（2021.02.01）二、潮流计算的基本算法1.牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。

这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。

牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。

因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。

而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。

1．2一般概念对于非线性代数方程组即()0,,,21=nixxxf ()ni,2,1= (1－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。

由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。

接着再从()1x 出发，重复上述计算过程。

因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()k k k x x x ∆+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。

由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。

潮流计算的基本算法及使用方法潮流计算的基本算法及使用方法一、潮流计算的基本算法1. 牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。

这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。

牛顿-拉夫逊法的基本原理是在解的某一邻域的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。

因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。

而所谓“某一邻域”是指雅可比方向均指向解的围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域。

1．2 一般概念对于非线性代数方程组()0=x f即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1)在待求量x 的某一个初始计算值()0x附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=?'+x x f x f (1－2)上式称之为牛顿法的修正方程式。

由此可以求得第一次迭代的修正量()[]()()0100x f x f x -'-=? (1－3)将()0x ?和()0x相加，得到变量的第一次改进值()1x 。

接着再从()1x 出发，重复上述计算过程。

因此从一定的初值()0x出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=?' (1－4)()()()k k k x x x ?+=+1 (1－5)上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。

由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。

潮流计算的基本算法及使用方法一、二、潮流计算的基本算法1.牛顿－拉夫逊法1．1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。

这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。

牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。

因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。

而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。

1．2一般概念对于非线性代数方程组即 ()0,,,21=n i x x x f ()n i ,2,1=(1－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (1－2)上式称之为牛顿法的修正方程式。

由此可以求得第一次迭代的修正量()()()[]()()0100x f x f x -'-=∆ (1－3)将()0x ∆和()0x 相加，得到变量的第一次改进值()1x 。

接着再从()1x 出发，重复上述计算过程。

因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为()()()()()k k k x f x x f -=∆' (1－4)()()()kkk xxx∆+=+1 (1－5)上两式中：()xf'是函数()x f对于变量x的一阶偏导数矩阵，即雅可比矩阵J；k为迭代次数。

由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。