第七章 潮流计算的数学模型及基本解法
第7章潮流计算的数学模型及基本解法
代入(7-10)式，经整理后有
Vn D I n Ys Vs L Vn U Vn 1
2013-1-10
（7-11）
12
考虑到电流和功率的关系式，上式写成迭代格式为
（ k 1）
Vi
i -1 n （ k） （ k） 1 Si Ys Vs Yij V j Yij V j i 1,2, ,n Yii j1 ji 1 k Vi
2013-1-10
7
综上所述，若选第N个节点为平衡节点，剩下n个节点（n=N-1） 中有r个节点是PV节点，则有n-r个节点是PQ节点。因此除了平衡 节点外，有n个节点注入有功功率，n-r个节点注入无功功率以及r 个节点的电压幅值是已知量。 在直角坐标系，待求的状态变量共2n个，用
x e
T
f
T T
2013-1-10 13
0
（7-12）
考所以，用 V j 代替 V j 可出得到更好的收敛果。 这就是高斯，赛德尔(Gauss-Seidel)选代的基本思 想，即一旦求出电压新值，在随后的迭代中立即使 用。这种方法的选代格式是
（ k 1）
（ k 1）
（ k）
Vi
（ k） （ k） i -1 n 1 Si Ys Vs Yij V j Yij V j i 1,2, ,n Yii j1 ji 1 k Vi
（7-13）
考高斯一赛德尔法比高斯迭代法收敛性要好。 考在导纳矩阵法的迭代公式中，导纳矩阵高度稀 疏，每行只有少数几个是非零元素，非对角非零元 素个数与和节点j相联的支路数相等。所以，上一次
潮流计算的基本算法及使用方法之欧阳音创编
潮流计算的基本算法及使用方法一、二、潮流计算的基本算法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 为迭代次数。
由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。
现代电力系统分析理论与方法 第7章 电力系统最优潮流
最优潮流计算
在系统的结构参数及负荷情况给定情况下，通过控制变量的优选， 找到能够满足所有给定的约束条件，并使系统的某一技术指标达到 最优（如网损、煤耗）时的潮流分布。
注：u为待选变量 约束条件分为等式约束条件和不等式约束条件。 采用的方法为：非线性规划
4
第一节
概述
随着电力系统规模扩大，对计算速度和系统安全性提出了更高要求，这 些经典调度理论已不能满足要求。将电力系统的潮流计算和优化理论结合， 并且计及系统的各种约束条件和电能质量，即形成了经典的优化理论—— 最优潮流(OPF)。OPF已在电力市场很多经济理论中广泛应用。
11
第二节
最优潮流的数学模型
考虑电力系统的经济因素，20世纪60年代末出现了一些经济调度理论， 例如最优分配有功负荷分布的等耗量微增率和无功电源最优分布的等网损 微增率。等耗量微增率准则是指系统所有发电机组具有同样的耗量微增率 时，系统运行所需要的费用最小，等网损微增率是指系统所有无功电源配 置具有相同的网损微增率时，系统网损最小。
最优潮 流的目 标函数
全系统火电机组燃料总费用，即 f Ki (PGi ) inG
式中：nG 为全系统所有发电机的集合，Ki (PGi ) 为第i台发 电机的耗量特性，一般用二次多项式表示，PGi 为第i台发电
机的有功出力。
有功网损，即 f (Pij Pji ) (i, j )nl 式中，nl 表示所有支路的集合。 9
可以证明最优潮流包含了等耗量微增率和等网损微增率，是这2个准则 在电力系统中的进一步发展运用（通过对目标函数的比较、约束条件的比 较、物理含义的分析等等）。
12
第三节
最优潮流的简化梯度算法
13
第三节
潮流计算的基本算法及使用方法
潮流计算的基本算法及使用方法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 为迭代次数。
潮流计算的基本算法及使用方法之欧阳家百创编
潮流计算的基本算法及使用方法一、欧阳家百（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.牛顿－拉夫逊法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 为迭代次数。
由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。
第七章 潮流计算的数学模型及基本解法
3 关于高斯法的讨论 非线性代数方程组 高斯迭代公式：
高斯法迭代的收敛性主要由
的谱半径[或矩阵 φ ( x* ) 的最大特征值]决定。当φ ( x* )的谱 半径小于1时高斯法可以收敛，φ ( x* ) 的谱半径越小收敛 性越好。
7.3 牛顿—拉夫逊法潮流计算
1. 牛顿拉夫逊法的一般描述
节点功率方程, x是节点电压。
2. 基于阻抗矩阵的方法 以平衡节点为电压给定节点的阻抗矩阵法
解法： (1) Yn 的稀疏因子表法+前代回代
(2) Z n法
(3) 高斯—赛德尔迭代
2．将接地支路用等效注入电流代替的阻抗矩阵法 节点导纳矩阵拆成不包含接地支路和只包含接地支 路的两部分。把平衡节点列写在最后
高斯迭代格式：
高斯—赛德尔迭代格式：
主要内容：从数学上说，潮流计算是要求解一组由 主要内容 潮流方程描述的非线性代数方程组。潮流计算问题 的数学模型；高斯迭代法(Gauss法)为基础的潮流计 算方法；牛顿—拉夫逊法潮流计算。
研究意义： 研究意义 确定电力系统稳态运行状态的方法之一。是电 力系统运行、规划以及安全性、可靠性分析和优化 的基础，也是各种电磁暂态和机电暂态分析的基础 和出发点。 对潮流计算方法的基本要求 基本要求： 基本要求 (1) 要有可靠的收敛性，对不同的系统及不同的运 行条件都能收敛； (2) 占用内存少、计算速度快； (3) 调整和修改容易，能满足工程上提出的各种要 求。
潮流计算结束时若平衡节点的有功功率无功功率和实际情况不符就要调整其它节点给定的边界条件以使平衡节点的功率在实际允许的范围之内
第二篇 电力系统潮流计算
主要内容：潮流计算的数学模型及基本解法；潮流方 主要内容 程的特殊解法；潮流计算中的特殊问题；潮流计算问 题的扩展。
电力系统潮流计算
电力系统潮流计算
Power System Load Flow Calculation 主要内容：
一． 二． 三． 四．
潮流计算的数学模型及解算方法 潮流方程的特殊解法 潮流计算中的特殊问题 潮流计算问题的扩展
一、潮流计算的数学模型及解算方法
什么是潮流计算？ 对电力系统正常运行状况的分析和 计算。通常需要已知系统参数和条件， 给定一些初始条件，从而计算出系统运 行的电压和功率等。 潮流计算有什么作用？ 确定系统的运行方式；系统规划设 计；稳定计算；故障计算等等。
四、潮流计算问题的扩展
3、约束方程 对控制变量的约束
– – – – – – –
发电机有功出力上、下限 发电机机端电压上、下限 无功电源控制电压上、下限 可投切电容、电抗容量上、下限 变压器变比上、下限 移相器可调相角上、下限 允许切除负荷容量上、下限
u
min
≤u≤u
max
四、潮流计算问题的扩展
3、约束方程 对依从变量的约束
– – –
负荷母线电压 发电机无功出力 线路有、无功潮流或电流
min
h
≤ h( x, u , D, p, A) ≤ h
max
四、潮流计算问题的扩展
4、常规潮流
f ( x, u , D, p, A) = 0 u h
minx
≤ h( x, u , D, p, A) ≤ h
四、潮流计算问题的扩展
四、潮流计算问题的扩展
1、变量的划分
网络结构关联阵A 网络结构关联阵A 网络元件参数p 干扰变量D 控制变量u 依从变量x
四、潮流计算问题的扩展
2、潮流方程
PGi − PDi − ∑ Pij (θ ,U , Y ) = 0 j∈i QGi − QDi − ∑ Qij (θ , U , Y ) = 0 j∈i f ( x, u, D, p, A) = 0
牛顿拉夫逊法潮流计算
摘要本文，首先简单介绍了基于在MALAB中行潮流计算的原理、意义，然后用具体的实例，简单介绍了如何利用MALAB去进行电力系统中的潮流计算。
众所周知，电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各线的电压、各元件中流过的功率、系统的功率损耗等等。
在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济性。
此外，在进行电力系统静态及暂态稳定计算时，要利用潮流计算的结果作为其计算的基础；一些故障分析以及优化计算也需要有相应的潮流计算作配合；潮流计算往往成为上述计算程序的一个重要组成部分。
以上这些，主要是在系统规划设计及运行方式安排中的应用，属于离线计算范畴。
牛顿－拉夫逊法在电力系统潮流计算的常用算法之一，它收敛性好，迭代次数少。
本文介绍了电力系统潮流计算机辅助分析的基本知识及潮流计算牛顿－拉夫逊法，最后介绍了利用MTALAB程序运行的结果。
关键词：电力系统潮流计算，牛顿－拉夫逊法，MATLABABSTRACTThis article first introduces the flow calculation based on the principle of MALAB Bank of China, meaning, and then use specific examples, a brief introduction, how to use MALAB to the flow calculation in power systems.As we all know, is the study of power flow calculation of power system steady-state operation of a calculation, which according to the given operating conditions and system wiring the entire power system to determine the operational status of each part: the bus voltage flowing through the components power, system power loss and so on. In power system planning power system design and operation mode of the current study, are required to quantitatively calculated using the trend analysis and comparison of the program or run mode power supply reasonable, reliability and economy.In addition, during the power system static and transient stability calculation, the results of calculation to take advantage of the trend as its basis of calculation; number of fault analysis and optimization also requires a corresponding flow calculation for cooperation; power flow calculation program often become the an important part. These, mainly in the way of system design and operation arrangements in the application areas are off-line calculation.Newton - Raphson power flow calculation in power system is one commonly used method, it is good convergence of the iteration number of small, introduce the trend of computer-aided power system analysis of the basic knowledge and power flow Newton - Raphson method, introduced by the last matlab run results.Keywords：power system flow calculation, Newton – Raphson method, matlab目录1 绪论 (1)1.1 课题背景 (1)1.2 电力系统潮流计算的意义 (2)1.3 电力系统潮流计算的发展 (2)1.4 潮流计算的发展趋势 (4)2 潮流计算的数学模型 (5)2.1 电力线路的数学模型及其应用 (5)2.2 等值双绕组变压器模型及其应用 (6)2.3 电力网络的数学模型 (8)2.4 节点导纳矩阵 (9)2.4.1 节点导纳矩阵的形成 (9)2.4.2 节点导纳矩阵的修改 (10)2.5 潮流计算节点的类型 (11)2.6 节点功率方程 (12)2·7 潮流计算的约束条件 (13)3 牛顿－拉夫逊法潮流计算基本原理 (14)3.1 牛顿-拉夫逊法的基本原理 (14)3.2 牛顿-拉夫逊法潮流计算的修正方程 (17)3.3 潮流计算的基本特点 (20)3.4 节点功率方程 (21)4牛顿－拉夫逊法分解潮流程序 (22)4·1 牛顿－拉夫逊法分解潮流程序原理总框图 (22)4.2 形成节点导纳矩阵程序框图及代码 (23)4.2。
潮流计算的基本算法及使用方法之欧阳数创编
潮流计算的基本算法及使用方法一、二、潮流计算的基本算法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 为迭代次数。
由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程式。
潮流计算
Sb SG STc S0c jQB 2 jQB3
1 b Tb 2 c Tc 3
A
d Td
SLDb
G
SG
SL D d
14
二、两级电压的开式电力网计算 计算方法一：包含理想变压器，计算时，经过理 想变压器功率保持不变，两侧电压之比等于实际 变比k。 T b d c L-1 L-2 SLD A
V1 arctg V1 V1
4
网络元件的功率损耗
功率损耗包括：电阻和等值电抗上的损耗 对地等值导纳上产生的损耗
V1S1 , I1 S ' I
jQB1
B j 2
R jX
S '', I S 2 , I 2 V2
jQB 2
B j 2
线路
VS1 , I1
线路
S0 (GT jBT )V 2
I0% S0 P0 jQ0 P0 j SN 100
开式网络的电压和功率分布计算
一、已知供电点电压和负荷点功率时的计算方法 已知末端的功率和电压：从末端开始依次计算出 电压降落和功率损耗。 已知电源点的电压和负荷的功率：采取近似的方 法通过叠代计算求得满足一定精度的结果
X2 k2 X2
T
A
A
B2 B2 / k 2 d c L-2 SLD
R'2+ j X'2 j B'2/2
16
R1+ jX1
j B1/2 j B1/2
b ΔS0
Z'T
c' j B'2/2
d'
SLD
二、两级电压的开式电力网计算 计算方法三：用π型等值电路代表变压器
(完整)潮流计算的概念和基本原理
潮流计算的概念和基本原理一、 潮流计算的意义电力系统潮流的计算和分析是电力系统运行和规划工作的基础。
运行中的电力系统，通过潮流计算可以预知,随着各种电源和负荷的变化以及网络结构的改变,网络所有母线的电压是否能保持在允许范围内,各种元件是否会出现过负荷而危及系统的安全，从而进一步研究和制订相应的安全措施。
规划中的电力系统，通过潮流计算,可以检验所提出的网络规划方案能否满足各种运行方式的要求，以便制定出既满足未来供电负荷增长的需求，又保证安全稳定运行的网络规划方案。
二、 潮流计算的基本概念潮流计算的一般提法是：已知电力网络的结构和参数，已知各负荷点、电源点吸取或发出的有功功率和无功功率（PQ 节点）,给定电压控制点的电压幅值和有功功率(PV 节点），对指定的一个平衡节点给定其电压幅值和相位角(V θ点）,求解全网各节点电压幅值和相位角，并进一步算出各支路的功率分布和网络损耗。
求解潮流问题的基本方程式是节点功率平衡方程。
三、 潮流计算的基本原理1. 潮流计算的基本模型1。
1潮流方程电力系统是由发电机、变压器、输电线路及负荷等组成，其中发电机及负荷是非线性元件，但在进行潮流计算时，一般可以用接在相应节点上的一个电流注入量来代表。
因此潮流计算所用的电力网络系由变压器、输电线路、电容器、电抗器等静止线性元件所构成，并用集中参数表示的串联或并联等值支路来模拟。
结合电力系统的特点，对这样的线性网络进行分析，普通采用的是节点法，节点电压与节点电流之间的关系V Y I= (1－1)其展开式为j nj ij i V Y I ∑==1),,3,2,1(n i = （1－2）在工程实际中,已经的节点注入量往往不是节点电流而是节点功率，为此必须应用联系节点电流和节点功率的关系式ii i i V jQ P I *-= ),,3,2,1(n i = （1－3） 将式(1－3）代入式（1－2）得到jnj ij iii V Y V jQ P ∑=*=-1),,3,2,1(n i = （1－4）交流电力系统中的复数电压变量可以用两种极坐标来表示i j ii e V V θ= （1－5）或ii i jf e V += （1－6）而复数导纳为ij ij ij jB G Y += （1－7)将式（1－6）、式(1－7）代入以导纳矩阵为基础的式(1－4），并将实部与虚部分开，可以得到以下两种形式的潮流方程。
第七章 潮流计算的数学模型及基本解法20110409概论
Eˆ I Sˆ Eˆ diag{Uˆ }是节点电压共轭组成的N N 阶对角线矩阵. i
Sˆ Eˆ YU
这是N阶非线性代数方程组.
上式展开后可写成:
P i
jQ i
Uˆ
i
Y U ij
j
ji
i 1,2,, N
❖ 在直角坐标系中的潮流方程
Pi jQi Uˆi YijU j ji
i 1,2,, N
❖ 潮流方程分析－1
对于N个节点的电力系统，选第N个节点为平衡节点，其余（n n N 1)
个节点中有r个PV节点，则有n r个PQ节点。
在直角坐标系中，待求的变量为2n个，
用x [eT f T ] [e1 e2
en f1 f2
直角坐标系中的潮流方程是：
fn ]表示。
Pi PiSP (eiai fibi ) 0 Qi QiSP ( fiai eibi ) 0
fibi eibi
i 1, 2, , N
❖ 极坐标系下的潮流方程
Pi jQi Uˆi YijU j
i 1,2,, N
ji
令Ui Uii代入基本方程得
Pi jQi Ui i (Gij jBij )U j j
ji
Ui (Gij jBij )U j (cosij jsinij )
i 1,2,, n r
i n r 1,, n
ΔP(e, f) Psp P(e, f)
f(x)
ΔQ(e,
f)
Qsp Q(e, f)
ΔU2(e, f) (U sp )2 U2 (e, f)
n维 n r维
r维
ΔP(e, f) Psp P(e, f)
f(x)
i 1, 2, , n i 1, 2, , n r
第7章潮流计算的数学模型及基本解法
2013-1-10
7
综上所述，若选第N个节点为平衡节点，剩下n个节点（n=N-1） 中有r个节点是PV节点，则有n-r个节点是PQ节点。因此除了平衡 节点外，有n个节点注入有功功率，n-r个节点注入无功功率以及r 个节点的电压幅值是已知量。 在直角坐标系，待求的状态变量共2n个，用
x e
T
f
T T
代入(7-10)式，经整理后有
Vn D I n Ys Vs L Vn U Vn 1
2013-1-10
（7-11）
12
考虑到电流和功率的关系式，上式写成迭代格式为
（ k 1）
Vi
i -1 n （ k） （ k） 1 Si Ys Vs Yij V j Yij V j i 1,2, ,n Yii j1 ji 1 k Vi
Vn Y（I n Ys Vs ）
-1 n
（7-14）
上式也可以写成
Vn Z（I n Ys Vs ） n
~
~
（7-15）
其中 Z n 是 Yn 的逆矩阵，即以平衡节点为电压 给定节点建立的节点阻抗矩阵。
二、关于高斯法的讨论
对于形如
f x 0
（7-16）
的非线性代数方程组，总可以写成 x x
故有
（7-4）
Pi ei ai f i bi Q f a e b i i i i i
Vi Vi
i 1,2,……, N
（7-5）
式（7-4）和式（7-5）是直角坐标系表示的潮流方程。 如果节点电压用极坐标表示，即令
高等电力系统分析-第七章-潮流计算的数学模型和基本解法
第七章潮流计算的数学模型和基本解法作业：7－4，7－8，7－91潮流计算物理上给定电力系统的网络结构、参数和决定系统运行状态的边界条件，求解电力系统的稳态运行状态。
数学上求解一组非线性代数方程组（潮流方程），与计算工具的发展息息相关。
目标收敛性好、占用内存少、计算速度快、易于调整和修改、使用灵活2潮流计算的发展历史50年代Y矩阵法；60年代初Z矩阵法；60年代Newton－Raphson法；60年代中Tinney稀疏矩阵技术；1974年 B Stott 提出快速分解法（Fast Decoupled Load Flow）；347.1 潮流问题的数学模型7.1.1 潮流方程YV =IˆˆiS =diag{V }YV ˆ()1,2,.....,i i i ij ij jj iP jQ V G jB V i N节点导纳方程实部虚部展开，共2N 个方程节点类型划分N个节点，每个节点有4个运行变量（P、Q、V、θ），全系统共有4N个变量，潮流方程2N个，所以要给定2N个变量，求解另外2N个。
PQ节点——给定PQ，求解Vθ负荷节点、联络节点、无AVR的发电机节点PV节点——给定PV，求解Qθ配备AVR的发电机节点Vθ节点——给定Vθ ，求解PQ计算上需要一般选择调节余量较大的发电机56lossP GP DP电力网络11,cos N Nloss iijijiji i j iP V P V V G11,cos NNloss ii j ij iji i j iQ V Q V V B,,,V Q V P loss loss 事先不知，所以，至少有一个节点的P 、Q 不能给定，来平衡总的损耗，该节点Vθ需要给定，叫Slack bus 。
松弛节点？——松弛什么？平衡节点？——平衡什么？系统中的网损都由平衡节点吸收了么？参考节点？——Vθ节点=参考节点么？Vθ节点如果计算得到的功率明显与物理不符合怎么办？扰动时的功率缺额应该由谁承担？按照什么原则选取Vθ节点？78直角坐标潮流方程ij j j ij ij i i i i jf e jB G jf e jQ Pi i i i ii ij j ij j j iP e a f b a G e B fi i i i ii ij j ij j j iQ f a e b b G f B e i = 1,2,…,Ni = 1,2,…,N指定一个节点为平衡节点, 例如N ，其给定。
潮流计算的基本算法及使用方法
潮流计算的基本算法及使用方法潮流计算的基本算法及使用方法一、潮流计算的基本算法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）可见，牛顿法的核心便是反复形成求解修正方程式。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
研究意义： 研究意义 确定电力系统稳态运行状态的方法之一。是电 力系统运行、规划以及安全性、可靠性分析和优化 的基础，也是各种电磁暂态和机电暂态分析的基础 和出发点。 对潮流计算方法的基本要求 基本要求： 基本要求 (1) 要有可靠的收敛性，对不同的系统及不同的运 行条件都能收敛； (2) 占用内存少、计算速度快； (3) 调整和修改容易，能满足工程上提出的各种要 求。
2. 基于阻抗矩阵的方法 以平衡节点为电压给定节点的阻抗矩阵法
解法： (1) Yn 的稀疏因子表法+前代回代
(2) Z n法
(3) 高斯—赛德尔迭代
2．将接地支路用等效注入电流代替的阻抗矩阵法 节点导纳矩阵拆成不包含接地支路和只包含接地支 路的两部分。把平衡节点列写在最后
高斯迭代格式：
高斯—赛德尔迭代格式：
在 x0 处将上式进行一阶台劳展开
f 潮流雅可比矩阵： J = T x
一般形式：
2. 直角坐标的牛顿—拉夫逊法
修正x( k ) 得x( k +1) 的新值。 重复上述过程直至max fi ( x( k ) ) < ε 。
3. 极坐标的牛顿—拉夫逊法
x = [θ U ]
T T T
雅可比矩阵的各部分子矩阵具有公式一致的形式
3 关于高斯法的讨论 非线性代数方程组 高斯迭代公式：
高斯法迭代的收敛性主要由
的谱半径[或矩阵 φ ( x* ) 的最大特征值]决定。当φ ( x* )的谱 半径小于1时高斯法可以收敛，φ ( x* ) 的谱半径越小收敛 性越好。
7.3 牛顿—拉夫逊法潮流计算
1. 牛顿拉夫逊法的一般描述
节点功率方程, x是节点电压。
第二篇 电力系统潮流计算
主要内容：潮流计算的数学模型及基本解法；潮流方 主要内容 程的特殊解法；潮流计算中的特殊问题；潮流计算问 题的扩展。
第七章 潮流计算的数学模型及基本解法
主要内容：从数学上说，潮流计算是要求解一组由 主要内容 潮流方程描述的非线性代数方程组。潮流计算问题 的数学模型；高斯迭代法(Gauss法)为基础的潮流计 算方法；牛顿—拉夫逊法潮流计算。
U i = ei + jf i
极坐标形式： 极坐标形式
2. 潮流方程的讨论和节点类型的划分 对于N个节点的电力系统，每个节点有四个运行 变量。共有4N个变量。实数方程为2N个。给定2N个， 求其余2N个。 全系统还应当满足功率平衡条件，即全网注入 功率之和应等于网络损耗
ห้องสมุดไป่ตู้
节点类型： 1．PQ节点：负荷 2．PV节点：发电机 3．节点：平衡节点（松弛节点） 平衡节点的选取是一种计算上的需要，有一定的 任意性。为使潮流计算结果符合实际，常把平衡节点 选在有较大调节余量的发电机节点。 潮流计算结束时若平衡节点的有功功率、无功功 率和实际情况不符，就要调整其它节点给定的边界条 件以使平衡节点的功率在实际允许的范围之内。
V N个节点的系统中： θ 节点1个，剩下共 n = N 1 个
节点， PV 节点
r 个，PQ节点n r 个。
直角坐标形式:待求量共2n个
极坐标形式：待求量共2n r 个
7.2 高斯迭代法为基础的潮流计算方法
1. 基于导纳矩阵的方法 将平衡节点s排在最后，取出前n个方程
高斯—赛德尔法 高斯 赛德尔法： 赛德尔法
发展状况： 发展状况 （1）高斯迭代法（导纳矩阵），收敛性极差； （2）牛顿—拉夫逊（阻抗矩阵），收敛性好， 内存占用量大； （3）牛顿—拉夫逊+Tinney提出的稀疏矩阵技 术和节点优化编号技术 （4）快速分解法，计算速度快，在线
7.1 潮流计算问题的数学模型
1. 潮流方程
直角坐标形式： 直角坐标形式 令