最优潮流问题的应用与计算.doc

内蒙古工业大学电力学院本科毕业设计论文 1 （ 二二 一一年年 六六 月月 本本科科毕毕业业设设计计说说明明书书 题题 目目： 最优潮流问题的应用与计算最优潮流问题的应用与计算 学学生生姓姓名名 ： 史史 界界 铭铭 学学 院院 ： 电电 力力 学学 院院 系系 别别： 电电 力力 工工 程程 系系 专专 业业： 电电 气气 工工 程程 及及 其其 自自 动动 化化 班班 级级： 电电 气气0 0 6 6 3 3 班班 指指导导教教师师 ： 韩韩 如如 月月讲师讲师/硕士硕士 学校代码：学校代码： 1012810128 学学 号：号： 200621202089200621202089 内蒙古工业大学电力学院本科毕业设计论文 2 摘要摘要 电力系统的运行必须在任何情况下，都应该尽可能的保证电力系统运行的可靠 性，保证电能的良好质量，保证运行的最大经济性。因此，做好规划设计就十分重 要。 在规划领域，需要进行潮流分析验证规划方案的合理性。潮流是确定电力网络 运行状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前提。无论进行 电力系统规划设计，还是对各种运行状态的研究分析，都须进行潮流计算。同时， 为了求使电力系统的某一指标达到最优值的潮流分布，法国学者在 carpentier 在 20 世纪 60 年代提出了最优潮流的概念。 电力系统最优潮流，简称 opf（optimal power flow） 。最优潮流是一个复杂的 非线性规划问题，要求在满足特定的电力系统运行和安全约束条件下，通过调整系 统中可利用控制手段实现预定目标最优的系统稳定运行状态。发展到今天，最优潮 流应用领域以十分广泛，针对不同的应用，最优潮流模型可以选择不同的控制变量、 状态变量集合，不同的目标函数以及不同的约束条件。本次毕业设计将以媒耗和网 损为目标函数，进行研究。 关键词：电力系统潮流计算；牛顿拉夫逊法潮流计算； matlab；最优潮流； 煤耗；网损 内蒙古工业大学电力学院本科毕业设计论文 3 abstract power system operation must be under any circumstances, should as far as possible the operation reliability of power system, ensure the quality, ensure good electrical maximum efficiency. therefore, the plan design is very important in the planning area, the need for current analysis verify the rationality of the plan. tidal power network operation is to determine the basic factors, current state of research power system is the basis and premise of steady-state problems. whether for power system planning and design, or for the operational status of research and analysis, are calculated to tide. at the same time, in order to make the power system for an index to achieve optimal value of the trend, french scholars in carpentier distribution in the 1960s put forward the concept of optimal fashion. power system optimal fashion, referred to as optimal power flow opf (.). the optimal fashion is a complex nonlinear programming problem, in particular requirements of power system operation and safety conditions, and by adjusting control system can realize the optimal objectives of stable operation system. today, the development of optimal power flow in a very wide application field, according to the different application, the optimal power flow model can choose different control variables, the state variables set, different target function and constraints. the graduation design in coal consumption and net loss for target function, were studied. keywords: the power flow calculation system, thenewton-raphson method computation, matlab, opf, coal, network loss 内蒙古工业大学电力学院本科毕业设计论文 4 目 录 第一章第一章 电力系统常规潮流计算电力系统常规潮流计算6 6 1.1 潮流计算简介 6 1.2 潮流计算的意义及其发展 6 1.3 潮流计算的数学模型潮流方程 7 1.4 电力系统节点分类 8 1.5 潮流计算的约束条件 9 1.6 牛顿-拉夫逊法 10 1.7 电力系统潮流计算算例： .17 第二章第二章 电力系统最优潮流计算电力系统最优潮流计算2424 2.1 最优潮流概述 .24 2.2 算例的最优潮流计算： .25 2.3 最优潮流与常规潮流的比较： .37 第三章第三章 最优潮流的应用最优潮流的应用4040 3.1 最优潮流的应用 .40 总结总结4343 参考文献参考文献4444 总总 结结4545 参考文献参考文献4646 谢谢 辞辞4747 内蒙古工业大学电力学院本科毕业设计论文 5 第一章 电力系统常规潮流计算 1.1 潮流计算简介 电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运 行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元 件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运 行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合 理性。可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳 定的基础。所以潮流计算是研究电力系统的一种很重要和基础的计算。 电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设 计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。 利用电子数字计算机进行电力系统潮流计算从 50 年代中期就已经开始。在这 20 年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计 算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点： （1）计算方法的可靠性或收敛性； （2）对计算机内存量的要求； （3）计算速度； （4）计算的方便性和灵活性。 电力系统潮流计算问题在数学上是一组多元非线性方程式求解问题，其解法都 离不开迭代。因此，对潮流计算方法，首先要求它能可靠地收敛，并给出正确答案。 由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大，潮流计算的方 程式阶数也越来越高，对这样的方程式并不是任何数学方法都能保证给出正确答案 的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。 1.2 潮流计算的意义及其发展 电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统 正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运 行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电 压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的 内蒙古工业大学电力学院本科毕业设计论文 6 运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分 析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最 优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采 用牛顿-拉夫逊法。 在运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需 要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了多个 在预想操作情况下电网的潮流分布以校验运行可靠性。在电力系统调度运行的多个 领域都涉及到电网潮流计算。潮流是确定电力网络运行状态的基本因素，潮流问题 是研究电力系统稳态问题的基础和前提。 在用数字计算机解电力系统潮流问题的开始阶段，普遍采取以节点导纳矩阵为 基础的逐次代入法。这个方法的原理比较简单，要求的数字计算机内存量比较下， 适应 50 年代电子计算机制造水平和当时电力系统理论水平。但它的收敛性较差，当 系统规模变大时，迭代次数急剧上升，在计算中往往出现迭代不收敛的情况。这就 迫使电力系统计算人员转向以阻抗矩阵为基础的逐次代入法。阻抗法改善了系统潮 流计算问题的收敛性，解决了导纳法无法求解的一些系统的潮流计算，在 60 年代获 得了广泛的应用。阻抗法的主要缺点是占用计算机内存大，每次迭代的计算量大。 当系统不断扩大时，这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺 点，60 年代中期发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分 割为几个小的地区系统，在计算机内只需要存储各个地区系统的阻抗矩阵及它们之 间联络线的阻抗，这样不仅大幅度地节省了内存容量，同时也提高了计算速度。 克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法。这是数学中解决非线性方程 式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵 为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性， 就可以大大提高牛顿法潮流程序的效率。自从 60 年代中期，在牛顿法中利用了最佳 顺序消去法以后，牛顿法在收敛性。内存要求。速度方面都超过了阻抗法，成为 60 年代末期以后广泛采用的优秀方法。 1.3 潮流计算的数学模型潮流方程 它具有如下特点： （1）它是一组代数方程，因而表征的是电力系统的稳定运行特性。 内蒙古工业大学电力学院本科毕业设计论文 7 （2）它是一组非线性方程，因而只能用迭代方法求其数值解。 （3）由于方程中的电压和导纳既可以表为直角坐标，又可表为极坐标，因而潮流方 程有多种表达形式-极坐标形式，直角坐标形式。 取 ， ，得到潮流方程的直角坐标形式： iii uejf ijijij ygjb （1-1） 11 11 ()() ()() nn iiijjijjiijjijj jj nn iiijjijjiijjijj jj peg eb ffg fb e qfg eb feg fb e 取， ，得到潮流方程的极坐标形式： iii uu ijijij ygjb （1-2） 1 1 (cossin) (sincos) n iijijijijij j n iijijijijij j puugb quugb 不同坐标形式的潮流方程适用于不同的迭代解法。 1.4 电力系统节点分类 用一般的电路理论求解网络方程，目的是给出电压源(或电流源)研究网络内的 电流(或电压)分布，作为基础的方程式，一般用线性代数方程式表示。然而在电力 系统中，给出发电机或负荷连接母线上电压或电流(都是向量)的情况是很少的，一 般是给出发电机母线上发电机的有功功率(p)和母线电压的幅值(u)，给出负荷母线上 负荷消耗的有功功率(p)和无功功率(q)。主要目的是由这些已知量去求电力系统内的 各种电气量。所以，根据电力系统中各节点性质的不同，很自然地把节点分成三类： pq 节点 对这一类点，事先给定的是节点功率(p，q)，待求的未知量是节点电压向量 (u，)，所以叫 pq 节点。通常变电所母线都是 pq 节点，当某些发电机的输出功 率 p。q 给定时，也作为 pq 节点。pq 节点上的发电机称之为 pq 机(或 pq 给定型 发电机)。在潮流计算中，系统大部分节点属于 pq 节点。 内蒙古工业大学电力学院本科毕业设计论文 8 pu 节点 这类节点给出的参数是该节点的有功功率 p 及电压幅值 u，待求量为该节点的 无功功率 q 及电压向量的相角。这类节点在运行中往往要有一定可调节的无功电 源。用以维持给定的电压值。通常选择有一定无功功率储备的发电机母线或者变电 所有无功补偿设备的母线做 pu 节点处理。pu 节点上的发电机称为 pu 机(或 pu 给 定型发电机) 平衡节点 在潮流计算中，这类节点一般只设一个。对该节点，给定其电压值，并在计算 中取该节点电压向量的方向作为参考轴，相当于给定该点电压向量的角度为零。也 就是说，对平衡节点给定的运行参数是 u 和，因此有城为 u节点，而待求量是 该节点的 p。q，整个系统的功率平衡由这一节点承担。 关于平衡节点的选择，一般选择系统中担任调频调压的某一发电厂(或发电机)， 有时也可能按其他原则选择，例如，为提高计算的收敛性。可以选择出线数多或者 靠近电网中心的发电厂母线作平衡节点。 以上三类节点 4 个运行参数 p。q。u。中，已知量都是两个，待求量也是两个， 只是类型不同而已。 1.5 潮流计算的约束条件 电力系统运行必须满足一定技术和经济上的要求。这些要求够成了潮流问题中 某些变量的约束条件，常用的约束条件如下： 1. 节点电压应满足 (1-3) minmax( 1,2,) iii uuuin 从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行 在额定电压附近。pu 节点电压幅值必须按上述条件给定。因此，这一约束条件对 pq 节点而言。 2. 节点的有功功率和无功功率应满足 内蒙古工业大学电力学院本科毕业设计论文 9 (1-4) minmax minmax gigigi gigigi ppp qqq pq 节点的有功功率和无功功率，以及 pu 节点的有功功率，在给定是就必须满 足上述条件，因此，对平衡节点的 p 和 q 以及 pu 节点的 q 应按上述条件进行检验。 3. 节点之间电压的相位差应满足 (1-5)max | | | ijijij 为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数 值。这一约束的主要意义就在于此。 因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束 条件。常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用约束条件 进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至修改系统的运行 方式，重新进行计算。 1.6 牛顿-拉夫逊法 电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统 正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运 行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电 压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的 运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分 析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最 优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采 用牛顿-拉夫逊法。 以下讨论的是用直角坐标形式的牛顿拉夫逊法潮流的求解过程。当采用直角 坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量 由于平衡节点的电压向量是给定的，因此待求两共需 12 12 ,., n n fff eee 2(1)n 要 2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作 用以外，其余每个节点都可以列出两个方程式。对 pq 节点来说，是给定 is is q p 和 内蒙古工业大学电力学院本科毕业设计论文 10 的，因而可以写出 （1-6） ()()0 ()()0 iijijiijjijj isjjj j ij i ijijijjjijj iisijj j ij i pfff eg egepbb qqfff g eegebb 对 pv 节点来说，给定量是，因此可以列出 isisvp 和 （1-7） 2 222 ()()0 ()0 iisijijiijjijj jij j ij i iisi i fff eg egeppbb f vve 求解过程大致可以分为以下步骤： （1）形成节点导纳矩阵 （2）将各节点电压设初值 u， （3）将节点初值代入相关求式，求出修正方程式的常数项向量 （4）将节点电压初值代入求式，求出雅可比矩阵元素 （5）求解修正方程，求修正向量 （6）求取节点电压的新值 （7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第 3 步重新开始 进行狭义次迭代，否则转入下一步 （8）计算支路功率分布，pv 节点无功功率和平衡节点柱入功率。 以直角坐标系形式表示 .迭代推算式 采用直角坐标时,节点电压相量及复数导纳可表示为: (1-8) iii ijijij vejf ygjb 将以上二关系式代入上式中,展开并分开实部和虚部;假定系统中的第 1,2,m 号 为 pq 节点,第 m+1,m+2,n-1 为 pv 节点,根据节点性质的不同,得到如下迭 代推算式: 内蒙古工业大学电力学院本科毕业设计论文 11 对于 pq 节点 (1-9) 11 11 ()() ()() nn iiiijjijjiijjijj jj nn iiiijjijjiijjijj jj ppeg eb ffg fb e qqfg eb feg fb e 1,2,im 对于 pv 节点 (1-10) 11 2222 ()() () nn iiiijjijjiijjijj jj iiii ppeg eb ffg fb e vvef 1,2,1immn 对于平衡节点 平衡节点只设一个,电压为已知,不参见迭代,其电压为: (1-11)nnn vejf .修正方程 式(2-3-5)和(2-3-6)两组迭代式工包括 2(n-1)个方程.选定电压初值及变量修正量 符号之后代入式(2-3-5)和(2-3-6),并将其按泰勒级数展开,略去二次方程, ii ef 及以后各项,得到修正方程如下: (1-12) wj u 内蒙古工业大学电力学院本科毕业设计论文 12 1 1 1 2 1 1 2 1 m m m m n n p q p q w p u p u 1 1 1 1 1 1 m m m m n n e f e f u e f e f 11111111 111111 11111111 111111 1111 mmmmnn mmmmnn mmmmmmm mmmmn pppppppp efefefef qqqqqqqq efefefef ppppppp efefefe j 11 111111 11111111 111111 22222 11111 11 m n mmmmmmmm mmmmnn mmmmmmmm mmmmnn mmmmm mm p f qqqqqqqq efefefef pppppppp efefefef uuuuu efefe 222 111 1111 11111111 111111 222222 111111 1111 mmm mmnn nnnnnnnn mmmmnn nnnnnn mmmm uuu fef pppppppp efefefef uuuuuu efefef 22 11 11 nn nn uu ef (1-13) 内蒙古工业大学电力学院本科毕业设计论文 13 .雅可比矩阵各元素的算式 式(3-2-8)中, 雅可比矩阵中的各元素可通过对式(3-2-4)和(3-2-5)进行偏导而求 得.当时, 雅可比矩阵中非对角元素为ji (1-14) 22 () 0 ii ijiiji jj ii ijiiji jj jj pq g eb f ef pq b eg f fe uu ef 当时,雅可比矩阵中对角元素为:ji (1-15) 1 1 1 1 2 2 () () () () 2 2 n i ijjijjiiiiii j i n i ijjijjiiiiii j j n i ijjijjiiiiii j i n i ijjijjiiiiii j j i i j i i i p g eb fg eb f e p g fb eg fb e f q g fb eg fb e e q geb fg eb f f u e e u f f 由式(3-2-9 和(3-2-10)看出,雅可比矩阵的特点: 矩阵中各元素是节点电压的函数,在迭代过程中,这些元素随着节点电压的变化而 内蒙古工业大学电力学院本科毕业设计论文 14 变化; 导纳矩阵中的某些非对角元素为零时,雅可比矩阵中对应的元素也是为零.若 ,则必有;0 ij y 0 ij j 雅可比矩阵不是对称矩阵;(1,2, ;)iqn is 雅可比矩阵各元素的表示如下: () () ()() ijiiji i ij ijjijjii iiii j j i g eb f jip h g eb fg eb fjie ) () ()() ijiiji i ij ijjijjii iiii j j i b eg f jip n g fb eb eg fjif ) () ()() ijiiji i ij ijjijjii iiii j j i b eg f jiq m g fb eb eg fjie （1-16） ) () ()() ijiiji i ij ijjijjii iiii j j i g eb f jiq l g eb fg eb fjif 2 0 () 2() i ij ij jiu r ejie 2 0 () 2() i ij ij jiu s fjif 内蒙古工业大学电力学院本科毕业设计论文 15 牛顿牛顿拉夫逊法的程序框图拉夫逊法的程序框图 内蒙古工业大学电力学院本科毕业设计论文 16 1.7 电力系统潮流计算算例： 系统如下图所示 p4=5 v4=1.05 4 1:1.05 j0.015 235 v5=1.05 ? 0.08+j0.30 0.04+j0.25 0.1+j0.35 3.7+j1.32+j1 j0.25 1.05:1 j0.25 j0.25 j0.25 1.6+j0.8 1 计算节点导纳矩阵： 29166 . 6 37874 . 1 35 . 0 1 . 0 1 25 . 0 04 . 0 1 25 . 0 13121011 j jj jyyyy 90015 . 3 62402 . 0 25 . 0 04 . 0 1 121221 j j yyy 64150 . 2 75471 . 0 35 . 0 1 . 0 1 131331 j j yyy 98082.6615390 . 1 052 . 1 1 015. 0 1 30. 008 . 0 1 25 . 0 04 . 0 1 )25. 025 . 0 ( 2 42 42 23212022 j jjj jj k y yyyy 11203 . 3 82987 . 0 30 . 0 08 . 0 1 3223 j j yy 49206.63 05 . 1 1 015 . 0 1 42 42 4224 j jk y yy 内蒙古工业大学电力学院本科毕业设计论文 17 66667.66 015 . 0 1 44 j j y 33333.33 03 . 0 1 55 j j y 节点导纳矩阵： 1.37874 -0.62402 -0.75471 0 0 -j6.29166 j3.90015 j2.64150 0 0 -0.62402 1.45390 -0.82987 0 0 j3.90015 -j66.98082 j3.11203 j63.19206 0 -0.75471 -0.92987 1.58459 0 0 j2.64150 j3.11203 -j35.73786 -j66.66667 0 0 0 0 0 0 0 j63.49206 0 -j66.66667 0 0 0 0 0 0 0 0 j31.74603 0 -j33.33333 matlab 程序： clear g(1,1)=1.37874;b(1,1)=-6.29166;g(1,2)=-0.62402;b(1,2)=3.90015; g(1,3)=-0.754 71;b(1,3)=2.64150;g(1,4)=0;b(1,4)=0; g(1,5)=0;b(1,5)=0; g(2,1)=-0.62402;b(2,1)=3.90015;g(2,2)=1.45390;b(2,2)=-66.98082;g(2,3)=- 0.82987;b(2,3)=3.11203; g(2,4)=0;b(2,4)=63.49206;g(2,5)=0;b(2,5)=0; g(3,1)=-0.75471;b(3,1)=2.64150;g(3,2)=- 0.82987;b(3,2)=3.11203;g(3,3)=1.58459;b(3,3)=-35.73786; g(3,4)=0;b(3,4)=0;g(3,5)=0;b(3,5)=31.74603; g(4,1)=0;b(4,1)=0;g(4,2)=0;b(4,2)=63.49206;g(4,3)=0;b(4,3)=0; g(4,4)=0;b(4,4)=-66.66667;g(4,5)=0;b(4,5)=0; g(5,1)=0;b(5,1)=0;g(5,2)=0;b(5,2)=0;g(5,3)=0;b(5,3)=31.74603; g(5,4)=0;b(5,4)=0;g(5,5)=0;b(5,5)=-33.33333; 内蒙古工业大学电力学院本科毕业设计论文 18 y=g+j*b delt(1)=0;u(1)=1.0;delt(2)=0;u(2)=1.0;delt(3)=0;u(3)=1.0;delt(4)=0; p(1)=-1.69;p(2)=-2.2;p(3)=-4.5;p(4)=5.5; q(1)=-0.8;q(2)=-1;q(3)=-1.3; k=0;precision=1; n1=3;%the n1 is the amount of the pq bus n2=1; while precision0.00001 u(4)=1.05; delt(5)=0; u(5)=1.05; for m=1:n1+n2 for n=1:n1+n2+1 pt(n)=u(m)*u(n)*(g(m,n)*cos(delt(m)-delt(n)+b(m,n)*sin(delt(m)-delt(n); end pp(m)=p(m)-sum(pt); end for m=1:n1 for n=1:n1+n2+1 qt(n)=u(m)*u(n)*(g(m,n)*sin(delt(m)-delt(n)-b(m,n)*cos(delt(m)-delt(n); end qq(m)=q(m)-sum(qt); end for m=1:n1+n2 for n=1:n1+n2+1 h0(n)=u(m)*u(n)*(g(m,n)*sin(delt(m)-delt(n)-b(m,n)*cos(delt(m)-delt(n); end h(m,m)=sum(h0)-u(m)2*(g(m,n)*sin(delt(m)-delt(n)-b(m,n)*cos(delt(m)-delt(n); h(m,m)=sum(h0)-u(m)2*(g(m,m)*sin(delt(m)-delt(m)-b(m,m)*cos(delt(m)- delt(m); jj(2*m-1,2*m-1)=h(m,m); end % 内蒙古工业大学电力学院本科毕业设计论文 19 for m=1:n1 for n=1:n1+n2+1 n0(n)=-u(m)*u(n)*(g(m,n)*cos(delt(m)-delt(n)+b(m,n)*sin(delt(m)-delt(n); j0(n)=-u(m)*u(n)*(g(m,n)*cos(delt(m)-delt(n)+b(m,n)*sin(delt(m)-delt(n); l0(n)=-u(m)*u(n)*(g(m,n)*sin(delt(m)-delt(n)-b(m,n)*cos(delt(m)-delt(n); end n(m,m)=sum(n0)-2*u(m)2*g(m,m)+u(m)2*(g(m,m)*cos(delt(m)- delt(m)+b(m,m)*sin(delt(m)-delt(m); j(m,m)=sum(j0)+u(m)2*(g(m,m)*cos(delt(m)-delt(m)+b(m,m)*sin(delt(m)- delt(m); l(m,m)=sum(l0)+2*u(m)2*b(m,m)+u(m)2*(g(m,m)*sin(delt(m)-delt(m)- b(m,m)*cos(delt(m)-delt(m); jj(2*m-1,2*m)=n(m,m); jj(2*m,2*m-1)=j(m,m); jj(2*m,2*m)=l(m,m); end for m=1:n1+n2 for n=1:n1+n2 if m=n else h(m,n)=-u(m)*u(n)*(g(m,n)*sin(delt(m)-delt(n)-b(m,n)*cos(delt(m)- delt(n); jj(2*m-1,2*n-1)=h(m,n); end end end for m=1:n1+n2 for n=1:n1 if m=n else n(m,n)=-u(m)*u(n)*(g(m,n)*cos(delt(m)-delt(n)+b(m,n)*sin(delt(m)-delt(n); 内蒙古工业大学电力学院本科毕业设计论文 20 jj(2*m-1,2*n)=n(m,n); end end end for m=1:n1 for n=1:n1+n2 if m=n else j(m,n)=u(m)*u(n)*(g(m,n)*cos(delt(m)-delt(n)+b(m,n)*sin(delt(m)-delt(n); jj(2*m,2*n-1)=j(m,n); end end end for m=1:n1 for n=1:n1 if m=n else l(m,n)=-h(m,n); jj(2*m,2*n)=l(m,n); end end end % for m=1:n1+n2 pp(2*m-1)=pp(m); end for m=1:n1 pp(2*m)=qq(m); end uu=-inv(jj)*pp; precision=max(abs(uu); for n=1:n1 内蒙古工业大学电力学院本科毕业设计论文 21 delt(n)=delt(n)+uu(2*n-1); u(n)=u(n)+uu(2*n); end n=n1+n2; delt(n)=delt(n)+uu(n); delt(n)=delt(n)+uu(7); k=k+1; end k-1,delt,u % the following program is used to calculate the s5 for n=1:5 u(n)=u(n)*(cos(delt(n)+j*sin(delt(n); end for m=1:5 i(m)=conj(y(5,m)*conj(u(m); end s5=u(5)*sum(i) for m=1:5 i(m)=conj(y(4,m)*conj(u(m); end s4=u(4)*sum(i) % the following program is used to calculate the smn y0(1,2)=j*0.25;y0(2,1)=j*0.25;y0(1,3)=0;y0(3,1)=0;y0(1,4)=0;y0(1,5)=0; y0(2,3)=j*0.25;y0(3,2)=j*0.25;y0(2,4)=-j*0.3308;y0(4,2)=j*0.315;y0(2,2)=0; y0(3,5)=-j*0.6615;y0(5,3)=j*0.63;y0(2,4)=0;y0(2,5)=0;y0(3,3)=0;y0(3,4)=0; y0(4,1)=0;y0(4,3)=0;y0(4,4)=0;y0(4,5)=0;y0(5,1)=0;y0(5,2)=0;y0(5,4)=0;y0(5,5)=0; for m=1:5 for n=1:5 s(m,n)=u(m)*(conj(u(m)-conj(u(n)*conj(-y(m,n)+u(m)2*conj(y0(m,n); end end 内蒙古工业大学电力学院本科毕业设计论文 22 s pg1=real(s4); pg2=real(s5); f=(50.4395*pg12+200.4335*pg1+1200.6485)+(200.55*pg22+500.746*pg2+1857.20 1) 计算结果： 发电机有功及无功出力 各节点电压向量 发电机序号有功出力无功出力 45.49992.0711 53.23612.5965 母线序号电压幅值电压相角 10.8373-0.0956 21.07460.3448 31.0292-0.0945 41.05000.4217 51.05000 支路有功功率 支路号首末端母线号 pijpji 112-1.46621.5845 213-0.13380.1569 3231.4156-1.2774 424-55 535-2.57942.5794 内蒙古工业大学电力学院本科毕业设计论文 23 支路有功功率 第二章 电力系统最优潮流计算 2.1 最优潮流概述 电力系统最优潮流，简称 opf（optimal power flow） 。最优潮流是一个复杂的 非线性规划问题，要求在满足特定的电力系统运行和安全约束条件下，通过调整系 统中可利用控制手段实现预定目标最优的系统稳定运行状态。发展到今天，最优潮 流应用领域以十分广泛，针对不同的应用，最优潮流模型可以选择不同的控制变量、 状态变量集合，不同的目标函数以及不同的约束条件。 2.2 最优潮流模型 最优潮流模型是在以下前提条件下提出的： (1)各火电核电投入运行的机组己知(不解决机组开停问题)。 (2)各水电机组的出力己定(由水库经济调度确定)。 (3)电力网络结构确定(不受接线方式影响，不考虑网络重构问题)。 最优潮流在数学上是一个带约束条件的优化问题，其中主要构成包括变量集合， 约束条件和目标函数。 最优潮流模型中，变量主要分为两大类。一类是控制变量，是可以控制的自变 量，通常包括各机组的有功出力，各发电机或同步补偿机无功出力；移相器抽头位 置、可调变压器抽头位置、并联电抗器或电容器容量；在某些紧急状况下，水电机 内蒙古工业大学电力学院本科毕业设计论文 24 组快速启动、某些负荷的卸载也可以作为控制手段。另一类是状态变量，是控制变 量的因变量，通常包括各节点电压和各支路功率等。 最优潮流考虑的系统约束条件有： (1)各节点有功功率和无功功率平衡约束 (2) 各发电机有功出力上下界约束。 （3）各发电机/同步补偿机无功出力上下界约束。 (4)并联电抗器/电容器容量约束。 (5) 移相器抽头位置约束。 （6）可调变压器抽头位置约束。 (7)各节点电压幅值上下界约束。 (8) 各支路传输功率约束。 从数学观点来看，以上约束条件中， （1）为等式约束，其余为不等式约束；从 约束的物理特性来看， （2）-（6）为控制变量约束， （7） （8）称为状态变量约束。 最优潮流有各式各样的目标函数，最常用的形式有以下两种： （1） 系统运行成本最小。该目标函数一般表示为火电机组燃料费用最小，不考虑 机组的启动、停机等费用。其中机组成本耗费曲线是模型的关键问题，它不 仅影响解的最优性，还制约求解方法的选取。 （2） 有功传输损耗最小。无功优化潮流通常以有功传输损耗最小为目标函数，它 在减少系统有功损耗的同时，还能改善电压质量。 2.2 算例的最优潮流计算： 系统如图所示示：系统如图所示示： 内蒙古工业大学电力学院本科毕业设计论文 25 p4=5 v4=1.05 4 1:1.05 j0.015 235 v5=1.05 ? 0.08+j0.30 0.04+j0.25 0.1+j0.35 3.7+j1.32+j1 j0.25 1.05:1 j0.25 j0.25 j0.25 1.6+j0.8 1 线路传输功率边界 支路号首末端母线号线路传输功率边界 1122 2130.65 3232 4246 5355 发电机数据 状态变量 10 个： x=(1 v1 2 v2 3 v3 4 v4 5 v5 ) (2-1) 系统中两台发电机为控制变量： 出力上界出力下界燃料耗费参数 发电机 序号 母线 号有功无功有功无功 二次系 数 一次系数常数 14831-350.4395200.43351200.6485 25851-2.1200.55500.7461857.201 内蒙古工业大学电力学院本科毕业设计论文 26 u=（pg1 pg2 qg1 qg2） (2-2) 系统总变量为 14 个 以煤耗为目标函数的最优潮流： 目标函数： (2-3) )()( 02212 2 22201111 2 121 apapaapapaf gggg 约束条件: 非发电机节点 0)sincos( 5 1 ijijij j ijjdii bgvvpp (i=1,2,3) (2-4)0)cossin( 5 1 ijijij j ijjdii bgvvqq 发电机节点 0)sincos( 5 1 ijijij j ijjdi ik gki bgvvppp （i=4,5） (2-5)0)cossin( 5 1 ijijij j ijjdi ik gki bgvvqqq 式中：表示第 k 台发电机接在节点 i 上，ik. 52, 41 kk 不等式约束条件共有 14 个，分别为 （i=1,2） (2-6) gigigi ppp (i=1,2) (2-7) gigi gi qqq (i=1,2,3,4,5) (2-8) iii vvv (对所有