1、电力系统分析课程设计基于 Matlab 的电力系统潮流计算专 业 电气工程及其自动化 班 级 1402 班 姓 名 郭娟 学 号 14039208 同 组 成 员 孙帆 郭娟 陈朝辉 苏新波 吕应发 指 导 教 师 楚冰清 完 成 时 间 2016 .12 .30 1目 录1 概述 .11.1 设计目的 .11.2 设计题目 .11.3 设计内容 .12 电力系统潮流计算概述 .22.1 电力系统简介 .22.2 潮流计算简介 .22.3 潮流计算意义及其发展 .23 Matlab 概述 .33.1 Matlab 简介 .33.2 Matlab 应用 .44 潮流计算过程 .44.1 系统图及
2、参数 .54.1.1 系统图 .54.1.2 各元件参数 .54.2 电网潮流计算思路 .54.3 潮流计算过程 .54.3.1 各元件参数计算 .54.3.2 绘制等效电路图 .74.3.3 功率分布计算 .75 Matlab 程序及运行结果 .105.1Matlab 程序 .105.2 运行结果 .206 设计心得与体会 .23参考文献 .2521 概述1.1 设计目的在如今的社会，电力已经成为人们必不可少的需求，而建立结构合理的大型电力系统不仅便于电能生产与消费的集中管理、统一调度和分配，减少总装机容量节省动力设施投资，且有利于地区能源资源的合理开发利用，更大限度地满足地区国民经济日益增
3、长的用电需要。电力系统建设往往是国家及地区国民经济发展规划的重要组成部分。电力系统的出现，使高效、无污染、使用方便、易于调控的电能得到广泛应用，推动了社会生产各个领域的变化，开创了电力时代，发生了第二次技术革命。电力系统的规模和技术水准已成为一个国家经济发展水平的标志之一。电力系统稳态分析包括潮流计算（或潮流分析)和静态安全分析。潮流计算针对电力革统各正常运行方式，而静态安全分析则要研究各种运行方式下个别系统元件退出运行后系统的状况。其目的是校验系统是否能安全运行，即是否有过负荷的元件或电压过低的母线等。原则上讲，静态安全分析也可 U 用潮流计算来代替。但是一般静态安全分析需要校验的状态数非常
4、多，用严格的潮流计算来分析这些状态往往计算量过大，因此不得不寻求一些特殊的算法以满足要求。牛顿法是数学中解决非线性方程式的典型方法，有较好的收敛性。解决电力系统潮流计算问题是以导纳距阵为基础的，因此，只要在迭代过程中尽可能保持方程式系数距阵的稀疏性，就可以大大提高牛顿法潮流程序的放率。自从 20 世纪 60 年代中期利用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、速度方面都超过了阻抗法，成为直到目前仍在广泛采用的优秀方法。1.2 设计题目电力系统潮流计算（牛顿-拉夫逊法）1.3 设计内容（1）根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；（2）赋予各节点电压变量（直角坐标系形式）
5、初值后，求解不平衡量；（3）形成雅可比矩阵；（4）求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗3和平衡节点功率；（5）上机编程调试；（6）计算分析给定系统潮流分析并与手工计算结果做比较分析；（7）书写课程设计报告。2 电力系统潮流计算概述2.1 电力系统简介 电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。电力系统潮流计算是电力系统稳定计算和故障分析的基础。潮流计算可以用传统的手工方式进行，也可以计算机软件完成。两种方法各有优缺点。手工
6、方式可用来计算一些接线较简单的电力网，但若将其用于接线复杂的电力网则计算量过大，难于保证计算准确性；计算机方式从数学上看可归结为用数值方法解非线性代数方程，数学逻辑简单完整，可快速精确地完成计算，但其缺点是物理概念不明显，物理规律被埋没在循环往复的数值求解过程中，基本原理不太明显。2.2 潮流计算简介潮流计算，电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡
7、点的电压和相位角。待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。2.3 潮流计算的意义及其发展电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算，即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。4对现有的电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和稳态分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法
8、有直接影响。实际电力系统的潮流技术那主要采用牛顿拉夫逊法。运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了多个在预想操作情况下电网的潮流分布以及校验运行可靠性。在电力系统调度运行的多个领域问题是研究电力系统稳态问题的基础和前提。在用数字见算计算机解电力系统潮流问题的开始阶段，普遍采取以节点导纳矩阵为基础的逐次代入法。这个方法的原理比较简单，要求的数字计算机内存量比较差下，适应 50 年代电子计算机制造水平和当时电力系统理论水平，但它的收敛性较差，当系统规模变大时，迭代次数急剧上升，在计算中往往出现迭代不收敛的
9、情况。这就迫使电力系统的计算人员转向以阻抗矩阵为基础的逐次代入法。阻抗法改善了系统潮流计算问题的收敛性，解决了导纳无法求解的一些系统的潮流计算，在 60 年代获得了广泛的应用，阻抗法德主要缺点是占用计算机内存大，每次迭代的计算量大。当系统不断扩大时，这些缺点就更加突出，为了克服这些缺点，60 年代中期发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需要存储各个地区系统的阻抗矩阵及它们之间联络的阻抗，这样不仅大幅度的节省了内存容量，同时也提高了计算速度。克服阻抗法缺点是另一个途径是采用牛顿-拉夫逊法。这是数学中解决非线性方程式的典型方法，有较好的收敛性
10、。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从 60 年代中期，牛顿法中利用了最佳顺序消去法以后，牛顿法在收敛性。内存要求。速度方面都超过了阻抗法，成为了 60 年代末期以后广泛采用的优秀方法。3 Matlab 概述3.1 Matlab 简介MATLAB 是由美国 mathworks 公司发布的主要面对科学计算、可视化以及5交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设
11、计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如 C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB 和 Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB 可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 3.2 Matlab 的应用MATLAB 的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形
12、式十分相似，故用 MATLAB 来解算问题要比用 C，FORTRAN 等语言完成相同的事情简捷得多，并且 mathwork 也吸收了像 Maple 等软件的优点,使MATLAB 成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN， C+ ，JAVA 的支持。可以直接调用 ,用户也可以将自己编写的实用程序导入到 MATLAB 函数库中方便自己以后调用，此外许多的 MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。MATLAB 产品族可以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理技术 数字信号处理技术 通
13、讯系统设计与仿真 财务与金融工程 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工6具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特定类型的问题 。 4 潮流计算过程4.1 系统图及参数4.1.1 系统图图 1 潮流计算设计图4.1.2 各元件参数变压器 T1、T2：SFL1-16000/110，（12122.5）/6.3，Ps=110kW，P0=10.5kW，U0=10.5，I0=0.9；变压器 T3：SFL1-8000/110 ，（1105 ）/6.6，Ps
14、=52kW，P0=12.76kW,Us=10.5，I0=1.1；变压器 T4：2SFL1-16000/110 ，（11022.5）/10.5，Ps=62kW，P0=11.6kW，Us=10.5，I0=1.10 ；导线型号均为 LGJ-150，参数 r0=0.21/km，x0=0.4/km，b0=2.810-6S/km。4.2 电网潮流计算思路（1）计算各元件参数，画出等值电路；（2）进行网络潮流计算；7（3）不满足供电要求，进行调压计算。4.3 潮流计算过程4.3.1 各元件参数计算120Km 线路SSlbBxXlrR461010 1036.128.4.5. （1）100Km 线路SSlbBx
15、XlrR462020 108.218.40. （2）70Km 线路SSlbBxXlrR463030 109.718.24. （3）变压器 T1，T2MVAjKVAjSIjPSUXSPRNkTNk )14.05.()1069.5.(10%82.962.1.01601 3322 32323 （4）变压器 T38MVAjKVAjSIjPSUXPRNKTN )08.1276.0()108.76.2(10%3.55.81.9081203 3232333 （5）变压器 T4 KVAjKVAjSIjPSUXRNKTN )176.0.()106.(10%4.795.3.21062104323324 （6）4.
16、3.2 绘制等效电路4.3.3 功率分布计算图 2 环形网络的等效图（1）各元件功率损耗1 234T1 T2T3T46.3/2 1/6.36./1010/.5S0=.185+j.4 S0=.185+j.4G1 G2S0=.16+j.7S0=.16+j.7LD4=1.2+j6LD3=5.6+j4S0=.976+j0.8RX R+jXT2RT3+jXT3RT4XT4RXT4R1+jXR3+jX2 4B1/2B/ B/jB/2jjB2/ j429两台 T4 变压器并联损耗MVAjjXRUSjTNT )468.0173.()05.346.1(02.)21(35.04240 （7）T 3 变压器损耗jj
17、STN )6.0412.()83.15.9(106.5)(232（8）100Km 与 70Km 线路交点 4 末端功率损耗 MVAjSSLDT 352.()0.73.(044 ）MVAj)2.705.1（9）120Km 与 100Km 线路交点 3 末端功率损耗 VAjSSLDT )4.08.6.()5097.412.(033 VAj)154.6.(（10）1.4 间 100Km 线路损耗 MVAjjUBjXRUSS MVAjBj jSNNLDLD )6305.1.0()197.2.4()()197.2.4)81.2.(0.5.1(7221414)(12414 MVAj)3.025.（11）1.3 间 120Km 线路损耗 )026.139.()302.56.2()()310.569.2)43.2.(4.1.(021121313)2(133 jUBjXRUSS MVAjBj jjS NNLDLD
