电力系统专业课程设计牛顿拉夫逊法潮流计算.doc

课程设计阐明书题目电力系统分析系(部)专业(班级)姓名学号指引教师起止日期电力系统分析课程设计任务书系（部）： 专业： 指引教师：课题名称电力系统专家潮流初步设计设计内容及规定1. 理解电力系统专家潮流计算基本原则2. 潮流计算不收敛因素分析3. 潮流计算收敛性分析4. 电力系统专家潮流计算流程图设计5. 分析成果设计工作量1. 掌握有关基本概念2. 理解潮流计算不收敛数学解释3. 对潮流计算收敛分析4. 设计计算流程图进度安排起止日期（或时间量）设计内容（或预期目的）备注第1天课题简介，收集有关材料，分析原始数据第2天学习有关基本理论第3天初步理解潮流计算收敛问题第4天流程图设计第5天编写设计阐明书教研室意见年 月 日系（部）主管领导意见年 月 日目录一、潮流计算基本原理1.1潮流方程基本模型1.2潮流方程讨论和节点类型划分1.3、潮流计算意义二、牛顿拉夫逊法2.1牛顿-拉夫逊法基本原理2.2节点功率方程2.3修正方程2.4牛顿法潮流计算重要流程三、收敛性分析四、算例分析总结参照文献电力系统分析潮流计算一、潮流计算基本原理1.1潮流方程基本模型电力系统是由发电机、变压器、输电线路及负荷等构成，其中发电机及负荷是非线性元件，但在进行潮流计算时，普通可以用接在相应节点上一种电流注入量来代表。因而潮流计算所用电力网络系由变压器、输电线路、电容器、电抗器等静止线性元件所构成，并用集中参数表达串联或并联等值支路来模仿。结合电力系统特点，对这样线性网络进行分析，普通采用是节点法，节点电压与节点电流之间关系 （11）其展开式为 （12）在工程实际中，已经节点注入量往往不是节点电流而是节点功率，为此必要应用联系节点电流和节点功率关系式 （13）将式（13）代入式（12）得到 （14）交流电力系统中复数电压变量可以用两种极坐标来表达 （15）或 （16）而复数导纳为 （17）将式（16）、式（17）代入以导纳矩阵为基本式（14），并将实部与虚某些开，可以得到如下两种形式潮流方程。潮流方程直角坐标形式为（18）（19）潮流方程极坐标形式为（110） （111）以上各式中，表达号后标号节点必要直接和节点相联，并涉及状况。这两种形式潮流方程普通称为节点功率方程，实牛顿拉夫逊等潮流算法所采用重要数学模型。1.2潮流方程讨论和节点类型划分对于电力系统中每个节点，要拟定其运营状态，需要由四个变量：有功注入注入有功、无功注入、电压幅值及电压相角。对于有个独立节点网络，其潮流方程有个，变量数为个。依照电力系统实际运营状况，普通每个节点4个变量中总有两个是已知，两个是未知。按各个节点所已经变量不同，可把节点提成三种类型。(1) 节点。此类节点已知节点注入有功功率、无功功率，待求未知量是节点电压值及相位角，因此称此类节点为节点。普通电力系统中没有发电设备变电所母线、发固定功率发电厂母线可作为节点，此类节点在电力系统中占大某些。(2) 节点。此类节点已经节点注入有功功率和电压值，待求未知量是节点注入无功功率及相位角，因此称此类节点为节点。此类节点普通为有一定无功功率储备发电厂母线和有一定无功功率电源变电所母线，此类节点在电力系统中位数不多，甚至可有可无。(3)平衡节点。潮流计算时，普通只设一种平衡节点，全网功率由平衡节点作为平衡机来平衡。平衡节点电压幅值及相位角是已知，如果给定、，待求则是注入功率、。1.3潮流计算意义早在20世纪50年代中期，就已开始使用数字计算机进行电力系统潮流计算。时至今日，潮流计算曾采用过各种不同办法，这些办法形成和发展都环绕着潮流计算某些基本规定进行。这些规定基本上可以归纳为如下几种方面：算法可靠性和收敛性、成果可信性；满足计算速度和内存占用量规定；计算以便灵活、适应性好。电力系统潮流计算和分析是电力系统运营和规划工作基本。运营中电力系统，通过潮流计算可以预知,随着各种电源和负荷变化以及网络构造变化,网络所有母线电压与否能保持在容许范畴内，各种元件与否会浮现过负荷而危及系统安全，从而进一步研究和制定相应安全办法。规划中电力系统，通过潮流计算，可以检查所提出网络规划方案能否满足各种运营方式规定，以便制定出既满足将来供电负荷增长需求，又保证安全稳定运营网络规划方案。二、牛顿拉夫逊法2.1牛顿-拉夫逊法基本原理设有单变量非线性方程 (4-1) 求解此方程时。先给出解近似值它与真解误差为，则将满足方程，即 (4-2)将(3-8)式左边函数在附近展成泰勒级数，于是便得 (4-3)式中，,分别为函数在处一阶导数，.，n阶导数。 如果差值很小，(3-9)式右端二次及以上阶次各项均可略去。于是，(3-9)便简化为 0 (4-4)这是对于变量修正量现行方程式，亦称修正方程式。解此方程可得修正量 (4-5)用所求去修正近似解，变得 (4-6)由于(3-10)是略去高次项简化式，因而所解出修正量也只是近似值。修正后近似解同真解依然有误差。但是，这样迭代计算可以重复进行下去，迭代计算通式是 (4-7)迭代过程收敛判据为 (4-8)或 (4-9)式中，为预先给定小正数。 这种解法几何意义可以从图31得到阐明。函数yf(x)为图中曲线。f(x)0解相称于曲线与x轴交点。如果第k次迭代中得到，则过点作一切线，此切线同x轴交点便拟定了下一种近似值。由此可见，牛顿拉夫逊法实质上就是切线法，是一种逐渐线性化办法。 应用牛顿法求解多变量非线性方程组(3-1)时，假定已给出各变量初值，. ，令，. 分别为各变量修正量，使其满足方程(3-1)即(4-10)将上式中n个多元函数在初始值附近分别展成泰勒级数,并略去具有,，二次及以上阶次各项，便得 (4-11)方程式(3-17)也可以写成矩阵形式 (4-12)方程式(3-18)是对于修正量,， 线性方程组,称为牛顿法修正方程式.运用高斯消去法或三角分解法可以解出修正量,，。然后对初始近似值进行修正 (i=1,2,.,n) (4-13)如此重复迭代，在进行k1次迭代时，从求解修正方程式 (4-14)得到修正量，并对各变量进行修正 (i=1,2,n) (4-15)式(3-20)和(3-21)也可以缩写为 (4-16)和 (4-17) 式中X和分别是由n个变量和修正量构成n维列向量；F(X)是由n个多元函数构成n维列项量；J是n阶方阵，称为雅可比矩阵，它第i、j个元素是第n个函数对第j个变量偏导数；上角标(k)表达阵每一种元素都在点处取值。 迭代过程始终到满足收敛判据 (4-18)或 (4-19)为止。和为预先给定小正数。2.2节点功率方程 电力系统负荷习惯用功率表达，对于有n个节点电力系统，系统中各节点注入电流与注入功率以标幺值表达关系为 i=1，2，n （3-20）式中表达其共轭复数。将此关系式代入节点电压方程通式，可得到以节点注入功率表达节点电压方程: (3-21) 上述方程式，普通称为功率方程。依照方程中节点电压向量表达不同，可以得到不同形式功率方程。 若节点电压向量以直角坐标表达，即以复数平面上实轴与虚轴上投影表达可写成 （3-22）其共轭值为 （3-23）导纳表达为 （3-24）把这两关系式代回式（3-21）功率方程中，展开后再将功率方程实部和虚某些别写成有功、无功功率分离节点方功率方程： （3-25）式中：i=1，2，n为各节点编号。若节点电压以极坐标表达，则或写成 （3-26） 将其同导纳复数表达式一起代入式（3-21）功率方程，进整顿可以得到 （3-27）式中：i与j节点电压相角差。由式（3-25）和（3-27）给出功率方程表达办法避免了复数运算，因而，在潮流计算中普遍采用。2.3修正方程采用牛顿法计算潮流时，需要对功率方程进行修改。下面将依照在不同坐标内修改进行讨论：（1）在直角坐标系内时，由PQ节点功率方程（3-25）可知：节点i注入功率是各点电压函数，设节点电压已知，代入式（3-25），可以求出节点i有功及无功功率，它们与给定PQ 节点注入功率差值应满足如下方程 （3-28）对于PV 节点，已知节点注入有功功率及节点电压大小，记作，其节点有功功率应满方程： （3-29）对于平衡节点，由于其电压给定，故不需要迭代求解。通过以上分析可见，式（3-28）和式（3-29）共2（n-1）个方程，待求量共2（n-1）个。将上述2（n-1）个方程按泰勒级数展开，并略去修正量高次方项后得到修正方程如下： （3-30）其中雅克比矩阵各元素可以对式（3-28）和式（3-29）求偏导数获得。对于非对角元素（）有 （3-31）对于对角元素（有 （3-32）由上述表达式可以看到，雅克比矩阵具备如下特点：（1）各元素是各节点电压函数，迭代过程中每迭代一次各节点电压都要变化，因而各元素每次也变化；（2）雅克比矩阵不具备对称性；（3）互导纳，与之相应非对角元素亦为零，此外因非对角元素，故雅克比矩阵是稀疏矩。当在极坐标系内时，由功率方程（3-27）可知节点i注入功率是各节点电压幅值和相角函数。代入式（3-27）可以求出节点i有功功率和无功功率，它们与给定PQ节点注入功率差值满足下面方程： （3-33）式中：i与j节点电压相角差。 在有n个节点系统中，假定第号节点为PQ节点，第m+1n-1号节点为PV节点，第n号节点为平衡节点。和是给定，PV节点电压幅值也是给定，因而，只剩余n-1个节点电压相角和m个节点电压幅值是未知量。由（3-33）可知一共包括了n-1+m方程式，正好同未知量数目相等，而直角坐标形式方程少了n-1-m个。由方程（3-33）可以写出修正方程 (3-34)式中 （3-35）其中：H是阶方阵，其元素为；N是阶矩阵，其元素为；K是阶矩阵，其元素为；L是阶矩阵，其元素为。 对式（3-33）求偏导数，可得雅克比矩阵元素表达式如下：非对角元素（） （3-36）对角元素（） （3-37）2.4牛顿法潮流计算重要流程（1）形成节点导纳矩阵；（2）给各节点电压设初值；（3）将节点电压初值代入（3-28）（3-29），求出修正方程式常数项向量；（4）将节点电压初值代入（3-31），（3-32），求出雅可比矩阵元素；（5）求解修正方程式（3-30），求出变量修正向量；（6）求出节点电压新值；（7）如有PV节点，则检查该类节点无功功率与否越限；（8）检查与否收敛，由式（3-19）可知，若电压趋近于真解时，功率偏移量将趋于零。如不收敛，则以各节点电压新值作为初值自第3步重新开始下一次迭代，否则转入下一步。（9）计算支路功率分布，PV节点无功功率和平衡节点注入功率，最后输出成果，并结束。牛顿-拉夫逊潮流计算程序框图如图3-2所示 三、收敛性分析小阻抗支路两端都是PQ节点时，其功率不平衡方程为：小阻抗支路电抗x非常小，它电纳（1/x）非常大。且1/x数值远不不大于上述方程组中出当代数量 ，为了以便阐明，不妨称这些代数量为小代数量。（1）第一次迭代忽视式（4.9）中小代数量。设立初始电压 ，可得：整顿式（4.13）可得这样可得第一次迭代后节点电压虚部关系：式（4.9）乘以k后，加上式(4.10)，设立初始电压为 ,可得:整顿式（4.16）可得第一次迭代后节点电压实部关系：式（4.9）乘以k后，加上式(4.10)，设立初始电压为 ,可得:式（4.11）乘以k后，加上式(4.12)，设立初始电压为 ,可得:综合式（4.15）、（4.17）可得，节点电压通过第一次迭代后满足收敛电压关系（3.14）。且可得新功率不平衡方程组（4.14）、（4.17）、（4.18）和（4.19）。这些方程中并不存在小阻抗问题，潮流计算可正常收敛。（2）第二次迭代忽视式（4.9）中小代数量，并将第一次迭代成果式（4.15）、（4.17）代入，可得：将式（4.15）代入式（4.20），可得忽视式（4.11）中小代数量，并将第一次迭代成果式（4.15）、（4.17）代入，可得：或者由式（4.21）、（4.23），可得四、算例分析采用3.4节5节点系统算例来论证以上结论，收敛精度为 。3种情形下迭代成果分别如各表所示，表格中 表达各个节点初给定值与计算值之间最大差值。计算成果表白：采用改进算法后，3种情形下，节点5和节点2电压值在各次迭代中保持关系： 情形1相应是小阻抗支路两端都为PQ节点状况，采用常规直角坐标牛顿法计算该情形潮流时，发散。采用改进办法，通过4次迭代后，收敛。情形2相应是小阻抗支路处在PQ节点和PV节点间状况，常规算法和改进算法都可收敛，切所需迭代次数相似。情形3相应是小阻抗支路处在PQ节点和平衡节点间状况，常规算法和改进算法都可收敛，且所需迭代次数相似。总结电力系统潮流计算分布计算，是指电力系统在某一稳定状态正常运营方式下，电力网络各节点电压和功率分布计算。它重要目：(1) 依照功率分布，可以选折电力系统电气设备和导线截面积，可觉得电力系统继电保护整定计算提供必要数据等。(2)检查电力系统各节点电压与否满足电压质量规定。(3)依照对各种运营方式潮流分布计算，可以协助咱们对的地选取系统接线方式，合理调节负荷，以保证电力系统安全、可靠地运营，向顾客供应高质量电能。(4) 检查电力系统各元件与否过负荷。(5)为电力系统规划和扩建提供根据。(6)为调节计算、经济运营计算、短路计算和稳定计算提供必要数据。潮流计算是电力系统分析中一种最基本计算，它任务是对给定运营条件拟定系统运营状态，如各母线上电压、网络中功率分布以及功率损耗等。潮流计算数学模型是以节点方程为基本，推导出相应功率方程。当电力系统中必须已知条件给定后潮流分布，取决于网络构造，而网络构造在功率方程中反映是节点导纳矩阵或节点阻抗矩阵。在复杂电力网中，在各个节点中没有直接相连节点诸多，从而使矩阵中有诸多零元素、节点导纳矩阵成为稀疏矩阵。普通来说，对角元素绝对值不不大于非对角元素绝对值，使节点导纳矩阵成为具备对角线优势矩阵。因而，节点导纳矩阵是一种对称、稀疏且具备对角线优势方阵。这给后来分析计算带来了很大以便，它有助于节约内存、提高计算速度以及改进收敛等。功率方程是非线性代数方程组，必要采用数值求解办法进行计算。在一方面讨论了电力系统节点分类以及潮流计算成果约束条件后，详细简介了惯用于潮流计算重要办法牛顿拉夫逊法。在实例计算这一章节中验证，从计算量来看，计算量不是特别庞大，计算成果可以迅速收敛，可以迅速精确计算出各个节点待求量。从效果来看，牛顿拉夫逊法迭代次数较少。从计算速度来看，速度比较快。在计算初始时，要选定比较适合初始值才干满足计算成果迅速收敛，如果选值不适当，计算成果不会收敛，也许成为发散型算式。参照文献1电力系统分析第二版，孟祥萍，高等教诲出版社，2当代电力系统分析，王锡凡，科学出版社，3电力系统分析，孙淑琴，机械出版社，4电力系统分析课程设计与综合实验，祝淑萍，中华人民共和国电力出版社，5 电力系统分析题解，何仰赞、温增银，华中科技大学出版社，