山东交通学院电力系统分析课程设计报告书院(部)别信息科学与电气工程学院班级学号姓名指导教师时间 2014.06.9-2013.06.13课程设计任务书题目复杂网络牛顿—拉夫逊法潮流分析院(部) 信息科学与电气工程学院专业电气工程及其自动化班级学生姓名学号6 月9 日至 6 月13 日共 1 周指导教师(签字)负责人(签字)年月日成绩评定表目录摘要 (2)1.1 潮流计算 (3)1.1.1潮流计算要求 (4)1.2 MATLAB简介 (5)2.1 实验设计要求 (5)2.2 电力系统等值网络 (5)2.2.1牛顿拉夫逊法程序流程图 (5)2.2.2 实验设计步骤 (7)2.2.3程序输出结果 (10)心得体会 (13)参考文献 (14)附录 (15)..摘要本课程设计通过对地区电网的设计，巩固和运用前面所学到的基础理论知识，掌握电力系统设计的一般原则和方法，培养分析问题和解决问题的能力。

在此过程中，要考虑到各方面的相互关系和相互影响，综合地运用课程中所学到的知识，进行独立思考。

潮流计算，电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。

潮流计算是电力系统分析中的一种最基本的计算，它的任务是对给定的运行条件确定系统的运行状态，如各母线上的电压（幅值及其相角）、网络中的功率分布及功率损耗等。

通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。

待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。

传运用MATLAB软件进行仿真潮流计算，图形界面直观，运行稳定，计算准确，提高了计算速度，各个类的有效封装又使程序具有很好的模块性.可维护性和可重用性。

因此，潮流计算在电力系统的规划计算、生产运行、调度管理及科学计算中都有着广泛的应用。

也就是说，对于电气工程及其自动化专业的学生来说，掌握潮流计算是非常重要和必要的。

关键字：潮流计算牛拉法Matlab1潮流计算1.1潮流计算概述潮流计算是研究电力系统稳态运行情况的一种基本电气计算，常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。

潮流计算的结果是电力系统稳定计算和故障分析的基础。

通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。

对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。

潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。

具体表现在以下方面：(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。

(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。

(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。

(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。

总结为在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。

同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。

在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。

此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。

所以潮流计算是研究电力系统的一种很重要和基础的计算。

1.1.1潮流计算的要求电力系统运行必须满足一定技术和经济上的要求。

这些要求够成了潮流问题中某些变量的约束条件，常用的约束条件如下：1.节点电压应满足min max (1,2,)i i i U U U i n ≤≤=从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在额定电压附近。

PU 节点电压幅值必须按上述条件给定。

因此，这一约束条件对PQ 节点而言。

2.节点的有功功率和无功功率应满足min max min max Gi Gi Gi Gi Gi Gi P P P Q Q Q ≤≤⎫⎬≤≤⎭PQ 节点的有功功率和无功功率，以及PU 节点的有功功率，在给定是就必须满足上述条件，因此，对平衡节点的P 和Q 以及PU 节点的Q 应按上述条件进行检验。

3.节点之间电压的相位差应满足max ||||||ij i j i j θθθθθ=-<-为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。

这一约束的主要意义就在于此。

因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束条件。

常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用约束条件进行检验。

如果不能满足要求，则应修改某些变量的给定值，甚至修改系统的运行方式，重新进行计算。

1 MATLAB 简介1.2 MATLAB 概述 MATLAB 是矩阵实验室（Matrix Laboratory ）的简称，是美国MathWorks 公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB 和Simulink 两大部分。

基本功能MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。

它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究必须进行有效数值计算的众多科学领域提供了一种全面的解决方案。

2.1实验设计要求下图示网络中，变压器的变比、各支路阻抗和一半的对地电纳均以标么值标于图中，设5节点为平衡节点，电压为1，节点4为PV节点，电压为1,P=0.5。

试求该网络的潮流分布，方法不限，求解精度为10e-5。

2.2电力系统等值网络2.2.1牛顿拉夫逊法程序流程图2.2.2实验设计步骤第一步，节点编号。

由已知条件可知G1为平衡节点，编号为5；G2为PV节点，编号为4；其余为PQ节点，编号为1，2，3；列出支路表与节点数据表如表1-1和1-2所示。

表1-1 支路数据表1-2 节点数据第二步，求节点导纳矩阵Y 11=j0.09+1/(0.03+j0.1)+j0.07+1/(0.025+j0.08)=6.311-j20.4652 Y 12=Y 21=-1/(0.025+j0.08)=-3.5587+j11.3879 Y 13=Y 31=-1/(0.03+j0.1)=-2.7523+j9.1743Y 22=1/(1.05*1.05*j0.1905)+j0.05+j0.07+1/(0.025+j0.08)+1/(0.02+j0.06)=8.5587-j31.0292Y 23=Y 32=-1/(0.02+j0.06)=-5+j15 Y 24=Y 42=-1/(1.05*j0.1905)=j4.9993Y 33=1/(1.05*1.05*j0.1905)+j0.09+1/(0.03+j0.1)+j0.05+1/(0.02+j0.06)=7.7523-j28.7956Y 35=Y 53=-1/(1.05*j0.1905)=j4.9993 Y 44=1/(j0.1905)=-j5.2493 Y 55=1/(j0.1905)=-j5.2493Y 14=Y 41=Y 15=Y 51=Y 25=Y 52=Y 34=Y 43=Y 45=Y 54=0 则节点导纳矩阵为B B B I Y U =第三步，设定初值。

给定个节点电压初始值(0)(0)e f 、 第四步，根据节点功率方程()()()()112222()j ni i i ij j ij j i ij j ij j j j ni i i ij j ij j i ij j ij j j i i i P P e G e B f f G f B e Q Q f G f B f e G f B e U U e f ====⎡⎤=--++⎣⎦⎡⎤=---+⎣⎦=-+∑∑ (1)式，求出修正方程式常数项向量(0)(0)2(0)P Q U 、、 △ P1=-0.8055 △Q1=-0.372;△P2=-0.18 △Q2=0.2475;△ P3=0 △Q3=-0.3875 △ P4=0.5 △U4=-2第五步：将电压初始值代入下式，求雅克比矩阵方程。

()()()()()2211=-=-+=-=-==0==i iij i ij i j j i iij i ij ij j i i j jj ni i i ij j ij j i ij j ij j j j ni i i ij j ij j i ij j ij j j i i j i P Q G e B f e f P Q B e G f f e U U e f j i P P e G e B f f G f B e Q Q f G f B f e G f B e P e ====≠∂∂∂∂∂∂∂∂∂∂∂∂⎡⎤=--++⎣⎦⎡⎤=---+⎣⎦∂∂∑∑当时，矩阵中非对角元素是当时，对角元素是()()()()1==1==1122---=-+-=+-=-++=-2e =-2j nij j ij j ii i ii i j j ni ij j ij j ii i ii i j i j ni ij j ij j ii i ii ij i j ni ij j ij j ii i ii ij i i ii i ijG e B f G e B f P G f B e B e G f f Q G f B e B e G f e Q G e B f G e B f f U e U f f ====-∂+∂∂+∂∂-∂∂∂∂∂∑∑∑∑(2)第六步：解修正方程式，求出修正量(0)(0)e f 、第七步：修正各节点电压(1)(0)(0)(1)(0)(0)=+=+e e ef f f⎧⎪⎨⎪⎩ 第八步：将(1)(1)e f 、再代入（1）式，求出(1)(1)2(1)P Q U 、、第九步：校验是否收敛，其收敛条件为(k)()()(x)=,k k f P Q ε﹤式中，()(),k k P Q 为向量()()k k P Q 、中大分量的绝对值。

这个收敛条件比较直观，它可以直接显示出最终结果的功率误差第十步：如果收敛，就进一步计算各段电力线路潮流和平衡节点功率，并打印出计算结果；如果不收敛，转回3步进行下一次迭代计算，直到收敛为止。