电力系统分析课程设计

电力系统分析课程设计学院电气工程学院班级学号姓名指导教师时间电力系统分析课程设计1前言电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和很基础的计算。随着科学技术的发展，电力系统变得越来越复杂，电气工程师掌握一种好的能对电力系统进行仿真的软件是学习和研究的需要。与众多专门的电力系统仿真软件相比，MATLAB软件具有易学、功能强大和开放性好，是电力系统仿真研究的有力工具。电力系统分析课程设计2目录1设计题目..........................................................31.1系统图的确定.................................................31.2各节点的初值及阻抗参数.......................................42潮流计算..........................................................52.1潮流计算概述与发展...........................................52.2复杂电力系统潮流计算.........................................52.3MATLAB概述..................................................62.4牛顿-拉夫逊法原理............................................72.5牛顿-拉夫逊法解决潮流计算问题................................82.6计算机潮流计算的步骤.........................................92.7计算机潮流计算流程图........................................113手算潮流计算.....................................................123.1确定节点类型................................................123.2求节点导纳矩阵Yb...........................................123.3计算各节点功率的修正方程的初始值（不平衡量）................134计算机算法潮流计算...............................................154.1计算机法潮流计算过程........................................154.2计算机法潮流计算结果........................................234.3系统功率分布图..............................................25个人心得...........................................................29附录：源程序.......................................................29参考文献...........................................................40电力系统分析课程设计31设计题目1.1系统图的确定选择六节点、环网、两电源和多引出的电力系统，简化电力系统图如图1-1所示，等值导纳图如图1-2所示。运用以直角坐标表示的牛顿-拉夫逊计算如图1-1所示系统中的潮流分布。计算精度要求各节点电压的误差或修正量不大于510。1.05:11.05:1U=1.05U=1.052+j13.7+j1.3①②③④⑤⑥1.8+j0.41.6+j0.80.08+j0.30.06+j0.0250.04+j0.250.1+j0.35L2L3L4L5图1-1电力系统图电力系统分析课程设计41.8+j4.01.6+j0.8j0.25j0.25j0.25j0.25j0.25j0.253.7+j1.32+j11.05:11.05:1oUj005.103.01①⑥⑤②③④0.1+j0.350.06+j0.0250.04+j0.25图1-2电力系统等值导纳图1.2各节点的初值及阻抗参数该系统中，节点①为平衡节点，保持1U=1.05+j0为定值，节点⑥为PV节点，其他四个节点都是PQ节点。给定的注入电压标幺值如表1-1、线路阻抗标幺值如表1-2、输出功率标幺值如表1-3。表1-1各节点电压标幺值参数U1U2U3U4U5U61.051.001.001.001.001.05表1-2线路、变压器阻抗标幺值线路L2L3L4L5T1T2Y/2阻抗0.06+j0.250.04+j0.250.08+j0.300.1+j0.35j0.03j0.015j0.25电力系统分析课程设计5表1-3节点输出功率节点②③④⑤⑥功率2+j11.8+j0.41.6+j0.83.7+j1.35注：各PQ节点的电压取1是为了方便计算和最后验证程序的正确性。2潮流计算2.1潮流计算概述与发展电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。牛顿-拉夫逊法作为一种实用的，有竞争力的电力系统潮流计算方法，是在应用了稀疏矩阵技巧和高斯消去法求修正方程后。牛顿-拉夫逊法是求解非线性代数方程有效的迭代计算。2.2复杂电力系统潮流计算电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行方式下的节点电压和功率分布，用以检查系统各元件是否过负荷、各点电压是否满足要求、功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计都是以潮流计算为基础。潮流计算结果的用途，例如用于电力系统稳定研究、安全估计或最优潮流等也对潮流计算的模型和方法有直接影响。节点类型：电力系统分析课程设计6（1）PV节点：柱入有功功率P为给定值，电压也保持在给定数值。（2）PQ节点：诸如有功功率和无功功率是给定的。（3）平衡节点：用来平衡全电网的功率。选一容量足够大的发电机担任平衡全电网功率的职责。平衡节点的电压大小与相位是给定的，通常以它的相角为参考量，即取其电压相角为0。一个独立的电力网中只设一个平衡点。基本步骤：（1）形成节点导纳矩阵；（2）将各节点电压设初值U；（3）将节点初值代入相关求式，求出修正方程式的常数项向量；（4）将节点电压初值代入求式，求出雅可比矩阵元素；（5）求解修正方程，求修正向量；（6）求取节点电压的新值；（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步；（8）计算支路功率分布，PV节点无功功率和平衡节点柱入功率。2.3MATLAB概述目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1）计算速度快；（2）内存需要少；（3）计算结果有良好的可靠性和可信性；（4）适应性好，亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强；（5）简单。MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。电力系统分析课程设计7MATLAB程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB与C语言和FORTRAN语言相比更容易被掌握。通过M语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。2.4牛顿-拉夫逊法原理假设有n个联立的非线性代数方程：0)......,(0)......,(0)......,(21212211nnnnxxxfxxxfxxxf假设以给出各变量的初值)0(1x，)0(2x，……，)0(nx，令其分别为个变量的修正量，使满足以上方程，所以：0),......,(0),......,(0),......,()0()0()0(2)0(2)0(1)0(1)0()0()0(2)0(2)0(1)0(12)0()0()0(2)0(2)0(1)0(11nnnnnnnxxxxxxfxxxxxxfxxxxxxf将上式中的n个多元函数在初始值附近分别展开成泰勒级数，并略去含有)0(1x，)0(2x，……，)0(nx的二次及以上阶次的各项，便得：0......),......,(0......),......,(0......),......,()0(303)0(202)0(101)0()0(2)0(1)0(3032)0(2022)0(1012)0()0(2)0(12)0(3031)0(2021)0(1011)0()0(2)0(11xxfxxfxxfxxxfxxfxxfxxfxxxfxxfxxfxxfxxxfnnnnnnn电力系统分析课程设计8方程可写成：)0()0(2)0(1002010202201201021011)0()0(2)0(1)0()0(2)0(12)0()0(2)0(11*..................),......,(),......,(),......,(nnnnnnnnnnnxxxxfxfxfxfxfxfxfxfxfxxxfxxxfxxxf以上方程是对于修正量)0(1x，)0(2x，……，)0(nx的线性方程组，称为牛顿法的修正方程，可解出)0(1x，)0(2x，……，)0(nx。对初始近似解进行修正：)0()0()1(iiixxx（i=1,2,……，n）反复迭代，在进行k+1次迭代时，从求解修正方程式：)()(2)(1212221212111)()(2)(1)()(2)(12)()(2)(11*..................),......,(),......,(),......,(knkkknnknknknkkknkkknkknknkkknkkxxxxfxfxfxfxfxfxfxfxfxxxfxxxfxxxf得到修正量)(1kx,)(2kx,……,)(knx,对各量进行修正)()()1(kikikixxx（i=1,2,……，n）迭代过程一直进行到满足收敛判据1)()(2)(1|}),,(max{|knkkixxxf，2.5牛顿-拉夫逊法解决潮流计算问题节点总