1、电 力 系 统 分 析 课程 设 计学 院 电气工程学院 班 级 学 号 姓 名 指导教师 时 间 电力系统分析课程设计1前言电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和很基础的计算。随着科学技术的发展，电力系统变得越来越复杂，电气工程师
2、掌握一种好的能对电力系统进行仿真的软件是学习和研究的需要。与众多专门的电力系统仿真软件相比，MATLAB 软件具有易学、功能强大和开放性好，是电力系统仿真研究的有力工具。电力系统分析课程设计2目录1 设计题目 .31.1 系统图的确定 .31.2 各节点的初值及阻抗参数 .42 潮流计算 .52.1 潮流计算概述与发展 .52.2 复杂电力系统潮流计算 .52.3 MATLAB 概述 .62.4 牛顿-拉夫逊法原理 .72.5 牛顿-拉夫逊法解决潮流计算问题 .82.6 计算机潮流计算的步骤 .92.7 计算机潮流计算流程图 .113 手算潮流计算 .123.1 确定节点类型 .123.2 求
3、节点导纳矩阵 Yb .123.3 计算各节点功率的修正方程的初始值（不平衡量） .134 计算机算法潮流计算 .154.1 计算机法潮流计算过程 .154.2 计算机法潮流计算结果 .234.3 系统功率分布图 .25个人心得 .29附录：源程序 .30参考文献 .40电力系统分析课程设计31 设计题目1.1 系统图的确定选择六节点、环网、两电源和多引出的电力系统，简化电力系统图如图 1-1 所示，等值导纳图如图 1-2 所示。运用以直角坐标表示的牛顿-拉夫逊计算如图 1-1 所示系统中的潮流分布。计算精度要求各节点电压的误差或修正量不大于 。 501 . 0 5 : 1 1 . 0 5 :
4、1U = 1 . 0 5 U = 1 . 0 52 + j 13 . 7 + j 1 . 31 . 8 + j 0 . 4 1 . 6 + j 0 . 80 . 0 8 + j 0 . 30 . 0 6 + j 0 . 0 2 50 . 0 4 + j 0 . 2 50 . 1 + j 0 . 3 5L 2L 3L 4L 5图 1-1 电力系统图电力系统分析课程设计41 . 8 + j 4 . 01 . 6 + j 0 . 8j 0 . 2 5j 0 . 2 5j 0 . 2 5j 0 . 2 5j 0 . 2 5j 0 . 2 53 . 7 + j 1 . 32 + j 11 . 0 5 :
5、 11 . 0 5 : 1oUj05.130.1+j0.350.06+j0.0250 . 0 4 + j 0 . 2 5图 1-2 电力系统等值导纳图1.2 各节点的初值及阻抗参数该系统中，节点为平衡节点，保持 =1.05+j0 为定值，节点为 PV 节1U点，其他四个节点都是 PQ 节点。给定的注入电压标幺值如表 1-1、线路阻抗标幺值如表 1-2、输出功率标幺值如表 1-3。表 1-1 各节点电压标幺值参数U1U 2U3U 4U5U 61.05 1.00 1.00 1.00 1.00 1.05表 1-2 线路、变压器阻抗标幺值 线路 L2 L3 L4 L5 T1 T2 Y/2阻抗 0.06
6、+j0.25 0.04+j0.25 0.08+j0.30 0.1+j0.35 j0.03 j0.015 j0.25电力系统分析课程设计5表 1-3 节点输出功率节点 功率 2+j1 1.8+j0.4 1.6+j0.8 3.7+j1.3 5注：各 PQ 节点的电压取 1 是为了方便计算和最后验证程序的正确性。2 潮流计算2.1 潮流计算概述与发展电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。利用电子数字计算机进行电力系统潮流计算从 50 年代中期就已经开始。在这 20 年内，潮流计算曾采用了各种不同的方法，
7、这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。牛顿-拉夫逊法作为一种实用的，有竞争力的电力系统潮流计算方法，是在应用了稀疏矩阵技巧和高斯消去法求修正方程后。牛顿-拉夫逊法是求解非线性代数方程有效的迭代计算。2.2 复杂电力系统潮流计算电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行方式下的节点电压和功率分布，用以检查系统各元件是否过负荷、各点电压是否满足要求、功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计都是以潮流计算为基础。潮流计算结果的用途，例如用于电力系统稳定研究、安全估计
8、或最优潮流等也对潮流计算的模型和方法有直接影响。节点类型：电力系统分析课程设计6（1）PV 节点：柱入有功功率 P 为给定值，电压也保持在给定数值。（2）PQ 节点：诸如有功功率和无功功率是给定的。（3）平衡节点：用来平衡全电网的功率。选一容量足够大的发电机担任平衡全电网功率的职责。平衡节点的电压大小与相位是给定的，通常以它的相角为参考量，即取其电压相角为 0。一个独立的电力网中只设一个平衡点。基本步骤：（1）形成节点导纳矩阵；（2）将各节点电压设初值 U；（3）将节点初值代入相关求式，求出修正方程式的常数项向量；（4）将节点电压初值代入求式，求出雅可比矩阵元素；（5）求解修正方程，求修正向量
9、；（6）求取节点电压的新值；（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第 3 步重新开始进行狭义次迭代，否则转入下一步；（8）计算支路功率分布，PV 节点无功功率和平衡节点柱入功率。2.3 MATLAB 概述目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1）计算速度快；（2）内存需要少；（3）计算结果有良好的可靠性和可信性；（4）适应性好，亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强；（5）简单。MATLAB 是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，
10、主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。MATLAB 程序设计语言结构完整，且具有优良的移植性，它的基本数据元素电力系统分析课程设计7是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB 与 C 语言和 FORTRAN 语言相比更容易被掌握。通过 M 语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。2.4 牛顿-拉夫逊法原理假设有 n 个联立的非线性代数方程： 0).,(.,)(2121nnnxff假设以给出各变量的
11、初值 ， ， ，令其分别为个变量的修)0()0(nx正量，使满足以上方程，所以： 0),.,(,., )( ()0()0(2)(01) ()()()()2 002011 nn nxxxff将上式中的 n 个多元函数在初始值附近分别展开成泰勒级数，并略去含有， ， 的二次及以上阶次的各项，便得：)0(1x)0(2)0(n 0.),.( .),.( 0)(30)0(2)0(10()0(21 )(302)0(2)0(120()0(212 )(1)(1)(1)()( xfxfxfxxf ffff xxxx nnnnnn 方程可写成：电力系统分析课程设计8 )0(2)(10021220101)0()0(
12、212)()( *.,.,. nnnnnn xfxffxfxxff 以上方程是对于修正量 ， ， 的线性方程组，称为牛)0()()(顿法的修正方程，可解出 ， ， 。对初始近似解进行修正：)(1x)(2)0(nx（i=1,2,，n）)0()()1(iiixx反复迭代，在进行 k+1 次迭代时，从求解修正方程式： )(2)(1212211)()(212)()( *.,.,. knkknknknkknknkk xfxffxfxxff 得到修正量 , , ,对各量进行修正 )(1k)(2k)(kn)()()1( kikii（i=1,2,，n）迭代过程一直进行到满足收敛判据 1)()(21|,max|
13、 knki xf，2.5 牛顿-拉夫逊法解决潮流计算问题节点总数为 n；PQ 节点有 m， ；PV 节点有 n-m-1，平衡节点有 1 个，节点编号按照先 PQ 节点，再 PV 节点，最后平衡节点的顺序进行编号，即:1，2，m 为 PQ 节点；m+1，m+2，n-1 为 PV 节点；n 为平衡节点。可形成结点导纳矩阵。导纳矩阵元素可表示为 ，本文中节点电压以ijijijBgY直角坐标形式表示，即 。由此下列公式可求出 Pi，QiiiijfeV.电力系统分析课程设计9nj jijiinj jijiii j jijiij jijiii eBfGefBeGfQP11 )()(假设系统中的第 1，2，
14、m 号节点为 PQ 节点，第 i 个节点的给定功率为和 ，对该节点可列方程：isPis ),2,1()()(11 mieBfGefBeGfQPnj nj jijiijijiiisiisi j j jijiijijiiisiisi 假设系统中的第 m+1，m+2，n-1 号节点为 PV 节点，则对其中每一个节点可列方程:)1,2,1()()(22211 nmifeVeBGBGePiiisiisinj nj jijiijijiiiiisi第 n 号节点为平衡节点，其电压为是给定的 ，故不参加迭代。nnjfeV.修正方程 可写成分块矩阵的形式：JW13211,2,1, ,331221,1321 *nnnnn VJJJ通过反复求解修正方程，解出各节点的未知量，再通过收敛判据判定是否已为真值。从而求得 PQ 节点的电压 V 及相角 的真值，PV 节点的 Q、 真值，平衡节点的 P、Q 真值，以上即为牛顿-拉夫逊迭代法的潮流计算过程，其优点为计算精确，运行速度快。其中的各个环节都可通过 MATLAB 程序来实现。2.6 计算机潮流计算的步骤（1）对电力网络的所有参数设初值,包括电压、相角、有功、无功等。（2）处理非标准变比支路,使其变成标准变比为 1 的变压器支路。（3）形成节点导纳矩阵 Y。
