电力系统潮流计算

武汉理工大学《电力系统分析》课程设计说明书1目录摘要...........................................................................................................2电力系统潮流计算........................................................................................31任务及题目要求.........................................................................................32原理介绍...................................................................................................32.1节点导纳矩阵...................................................................................42.2牛顿-拉夫逊法.................................................................................52.2.1牛顿-拉夫逊法基本原理...........................................................52.2.2牛顿--拉夫逊法潮流求解过程介绍.............................................73分析计算.................................................................................................114结果分析.................................................................................................165总结.......................................................................................................17参考资料...................................................................................................18武汉理工大学《电力系统分析》课程设计说明书2摘要电力网的运行状态可用节点方程或回路方程来描述。节点导纳矩阵是以系统元件的等值导纳为基础所建立的、描述电力网络各节点电压和注入电流之间关系的线性方程。潮流计算是电力系统分析中的一种最基本的计算，它的任务是对给定的运行条件确定系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布及功率损耗等。本文就节点导纳矩阵和潮流进行分析和计算。关键词：潮流计算，导纳矩阵，等值导纳，电力系统武汉理工大学《电力系统分析》课程设计说明书3电力系统潮流计算1任务及题目要求题目初始条件：如图所示电网。其元件导纳参数为：y12=0.5-j3,y23=0.8-j4,y13=0.75-j2.5任务及要求:1）根据给定的运行条件，确定图2所示电力系统潮流计算时各节点的类型和待求量；2）求节点导纳矩阵Y；3）给出潮流方程或功率方程的表达式；4）当用牛顿-拉夫逊法计算潮流时，给出修正方程和迭代收敛条件。2原理介绍110.51.02+j1j1.123y13y23y12武汉理工大学《电力系统分析》课程设计说明书42.1节点导纳矩阵节点导纳矩阵既可根据自导纳和互导纳的定义直接求取，也可根据电路知识中找出改网络的关联矩阵，在节点电压方程的矩阵形式进行求解。本章节我们主要讨论的是直接求解导纳矩阵。根据节点电压方程章节我们知道，在利用电子数字计算机计算电力系统运行情况时，多采用IYV形式的节点方程式。其中阶数等于电力网络的节点数。从而可以得到n个节点时的节点导纳矩阵方程组：11112211211222221122nnnnnnnYYYnYYYnYYYnVVVIVVVIVVVI（2-1）由此可以得到n个节点导纳矩阵：111212212212nnnnnnYYYYYYYYYY（2-2）它反映了网络的参数及接线情况，因此导纳矩阵可以看成是对电力网络电气特性的一种数学抽象。由导纳短阵所联系的节点方程式是电力网络广泛应用的一种数学模型。通过上面的讨论，可以看出节点导纳矩阵的有以下特点：（1）导纳矩阵的元素很容易根据网络接线图和支路参数直观地求得，形成节点导纳矩阵的程序比较简单。（2）导纳矩阵为对称矩阵。由网络的互易特性易知jjiYiY。武汉理工大学《电力系统分析》课程设计说明书5（3）导纳矩阵是稀疏矩阵。它的对角线元素一般不为零，但在非对角线元素中则存在不少零元素。在电力系统的接线图中，一般每个节点与平均不超过3~4个其他节点有直接的支路连接。因此，在导纳矩阵的非对角线元素中每行仅有3~4个非零元素，其余的都是零元素，而且网络的规模越大，这种现象越显著。节点导纳矩阵的形式可归纳如下：（1）导纳矩阵的阶数等于电力网络（2）导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连得不接地支路数。（3）导纳矩阵各对角元素，即节点的自导纳等于相应节点之间的支路导纳之和。（4）导纳矩阵非对角元素，即节点之间的互导纳等于相应节点之间的支路导纳的负值。2.2牛顿-拉夫逊法2.2.1牛顿-拉夫逊法基本原理牛顿--拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程。即通常所称的逐次线性化过程。对于非线性代数方程组：()0fx即12(,,,)0infxxx(1,2,,)in(2-3)在待求量x的某一个初始估计值(0)x附近，将上式展开成泰勒级数并略去二阶及以上的高阶项，得到如下的经线性化的方程组：(0)'(0)(0)()()0fxfxx(2-4)武汉理工大学《电力系统分析》课程设计说明书6上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量(0)'(0)1(0)[()]()xfxfx(2-5)将(0)x和(0)x相加，得到变量的第一次改进值(1)x。接着就从(1)x出发，重复上述计算过程。因此从一定的初值(0)x出发，应用牛顿法求解的迭代格式为：'()()()()()kkkfxxfx(2-6)(1)()()kkkxxx(2-7)上两式中：'()fx是函数()fx对于变量x的一阶偏导数矩阵，即雅可比矩阵J；k为迭代次数。由上式可见，牛顿法的核心便是反复形式并求解修正方程式。牛顿法当初始估计值(0)x和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。牛顿潮流算法突出的优点是收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代4~5次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。牛顿法也具有良好的收敛可靠性，对于对以节点导纳矩阵为基础的高斯法呈病态的系统，牛顿法也能可靠收敛。牛顿法所需的内存量及每次迭代所需时间均较高斯法多。牛顿法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有可能根本不收敛或收敛到一个无法运行的节点上。对于正常运行的系统，各节点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值(也称为平直电压)，如假定：(0)1iU(0)0i或(0)1ie(0)0if(1,2,,;)iqnis(2-8)这样一般能得到满意的结果。但若系统因无功紧张或其它原因导致电压质量很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问题。解决这个问题的办法可以用高斯法迭代1~2次，以此迭代结果作为牛顿法武汉理工大学《电力系统分析》课程设计说明书7的初值。也可以先用直流法潮流求解一次以求得一个较好的角度初值，然后转入牛顿法迭代。2.2.2牛顿--拉夫逊法潮流求解过程介绍以下讨论的是用直角坐标形式的牛顿—拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量1212,,,...,nnfffeee由于平衡节点的电压向量是给定的，因此待求两共2(n-1)需要2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。对PQ节点来说，isP和isQ是给定的，因而可以写出()()0()()0iijijiijjijjisjjjjijiijijijjjijjiisijjjijipfffeGeGePBBQQfffGeeGeBB（2-9）对PV节点来说，给定量是isisVP和，因此可以列出2222()()0()0iisijijiijjijjjijjijiiisiifffeGeGePPBBfVVe（2-10）求解过程大致可以分为以下步骤：（1）形成节点导纳矩阵；（2）将各节点电压设初值U（3）将节点初值代入相关求式，求出修正方程式的常数项向量；（4）将节点电压初值代入求式，求出雅可比矩阵元素；（5）求解修正方程，求修正向量；（6）求取节点电压的新值；（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步；武汉理工大学《电力系统分析》课程设计说明书8（8）计算支路功率分布，PV节点无功功率和平衡节点注入功率。以直角坐标系形式表示：○1迭代推算式采用直角坐标时,节点电压相量及复数导纳可表示为:iiiijijijVejfYGjB(2-11)将以上二关系式代入上式中,展开并分开实部和虚部;假定系统中的第1,2,,m号为P—Q节点,第m+1,m+2,,n-1为P—V节点,根据节点性质的不同,得到如下迭代推算式:⑴于PQ节点1111()()()()nniiiijjijjiijjijjjjnniiiijjijjiijjijjjjPPeGeBffGfBeQQfGeBfeGfBe1,2,,im(2-12)⑵对于PV节点112222()()()nniiiijjijjiijjijjjjIiiiPPeGeBffGfBeVVef1,2,,1immn(2-13)⑶对于平衡节点平衡节点只设一个,电压为已知,不参见迭代,其电压为:nnnVejf(2-14)○2修正方程两组迭代式中包括2(n-1)个方程.选定电压初值及变量修正量符号之后代入,并将其按泰勒级数展开,略去,iief二次方程及以后各项,得到修正方程如下:WJU(2-15)武汉理工大学《电力系统分析》课程设计说明书9其中，11121121mmmmnnPQPQWPUPU；111111mmmmnnefefUefef11111111111111111111111111111111mmmmnnmmmmnnmmmmmmmmmmmnPPPPPPPPefefefefQQQQQQQQefefefefPPPPPPPefefefeJ