《电力系统暂态分析课程设计.doc》由会员分享，可在线阅读，更多相关《电力系统暂态分析课程设计.doc（48页珍藏版）》请在三一办公上搜索。
1、内蒙古科技大学电力系统稳态分析课程设计题 目：基于MALB的电力系统复杂潮流分析学生姓名:刘建峰学 号：1173020专 业:电气工程及其自动化班 级：电气201班指导教师：刘景霞摘要电力系统潮流计算是电网分析的基础应用，是对复杂电力系统正常和故障条件下稳态运行状态的计算。给定电力系统的网络结构、参数和决定电力系统运行状况的边界条件,确定电力系统运行的方法之一是朝流计算。MATLAB是一种交互式、面向对象的程序设计语言,广泛应用于工业界与学术界，主要用于矩阵运算.采用迭代法,通过建立矩阵的修正方程来依次迭代,逐步逼近真值来计算出电力网的电压，功率分布。PQ分解法是极坐标形式牛顿-拉潮流计算的一
2、种简化计算方法，。Q分解法通过对电力系统具体特点的分析，对牛顿法修正方程式的雅可比矩阵进行了有效的简化和改进。由于这些简化只涉及修正方程式的系数矩阵，并未改变节点功率平衡方程和收敛判据，因不会降低计算结果的精度。用手算和计算机算法对其进行设计。使用MATLB软件进行编程,在很大程度上节省了内存,减少了计算量。通过对本题计算我们了解了一些工程计算和解决工程问题的方法。关键词：潮流计算，PQ分解法,MATLAB Electriclpowr ste ompletidal curren alyss asd ATLAowe Flow Anays Gidmpting is thebasiof applic
3、atns， th omlexpowersystem udr norma and faultcdtions fr h caulato f teady st orato. Gin e poer yt netoksructu， pmeters ad dcisions oprtion of hpowerystm bdary cnditions,to dtermine t mthod fratinoft pwr system isone of Nort Korealow cculo.ML is an eracive,object-oriente progming angug, wely usdinind
4、utry nd cadeia, maiyforatrixaculation. sing ieration, eamenm trouh he ealisment of matrix iteratie equion totur,gradallymoving towrd a truelto calate the voltage lecrc grid, oer distrbuin.PQ decmpositio mthod is thform opolr coordinates Newton thewideg rend of silifi cculatinmehod. P-Qdcomposin metd
5、 apted by thsecific cractistisof the power systm aays, Nwos L of the Jacbian mtrxformulahs effectvel simlified nd iprovd. Aa result of the imlifed formutt involeonly the coefficien marx，the baceooerhas notcaged node eqatonsan he ovrc rieron， becaue theelts wll no rdce e acuac.Use MATLABsoftwar proga
6、mmg, saving memoryto are xtent, rduc teaount of computtion. y thscalculati we undrstand tht numbe engierng calcultiand sleengineringroblemKeord ： Th red， P decomposii metod,MATLAB目录内蒙古科技大学课程设计任务书5第一章引言1.1 研究背景及意义71潮流计算的意义81.3电力系统稳态分析潮流计算总结91.4 MATAB的概述9第二章 理论计算1.1-Q法潮流计算的基本步骤112. PQ分解法潮流计算流程图12.3 两机
7、五节点网络潮流计算1第三章 程序设计1831设计流程13.2 程序设计183.3 程序运行结果22第四章 设计感想参考文献45内蒙古科技大学课程设计任务书课程名称电力系统稳态分析课程设计设计题目两机五节点网络潮流计算 PQ法指导教师刘景霞时间1周一、教学要求电力系统稳态分析课程设计以设计和优化电力系统的潮流分析为重点，提高学生综合能力为目标，尽可能结合实际工程进行。设计内容的安排要充分考虑学校现有的设备，设计时间及工程实际需要，并使学生初步学会运用所学知识解决工程中的实际问题。二、设计资料及参数(一）设计原始资料1、待设计电气设备系统图2、电力系统网络各元件参数3、电力系统电气元件的使用规范4
8、、电力工程电气设计手册(二)设计参考资料、电力系统稳态分析，陈珩,中国电力出版社，2007,第三版、电力系统分析，韩祯祥，浙江大学出版社，005,第三版3、电力系统分析课程实际设计与综合实验,祝书萍，中国电力出版社，07,第一版三、设计要求及成果.根据给定的参数或工程具体要求,收集和查阅资料；学习相关软件（软件自选)。2在给定的电力网络上画出等值电路图。.运用计算机进行潮流计算。.编写设计说明书。基本要求：1.编写潮流计算程序;在计算机上调试通过(？）；3.运行程序并计算出正确结果（？）；4.写出课程设计报告(包括以下内容）(1份）(1)程序框图；(2)源程序;（）符号说明表;（）算例及计算结
9、果.编写计算说明书(1份)。四、进度安排根据给定的参数或工程具体要求,收集和查阅资料(半天）学习软件（MATLA或C语言等)（一天半)编程计算复杂系统潮流计算(三天)编写计算设计书（一天）五、评分标准课程设计成绩评定依据包括以下几点:1） 工作态度（占1%);2） 基本技能的掌握程度（占0%）；3） 程序编写是否合理是否有运行结果(0%);4） 课程设计说明书编写水平(占30%）。5） 分为优、良、中、合格、不合格五个等级。考核方式:设计期间教师现场检查；评阅设计报告第一章 引言11 研究背景及意义电力系统在运行时，在电源电势激励作用下,电流或功率从电源通过系统各元件流入负荷，分布于电力网各处
10、,称为潮流分布。潮流计算是电力系统分析中的一种最基本的计算，即根据给定电力系统的网络结构和参数,在满足电力系统运行状态参数的边界条件情况下,确定电力系统稳态运行状态的最基本的方法。它的任务是在给定的接线方式和运行条件下,确定系统的运行状态,如各母线上的电压（幅值和相角）、网络中的功率分布及功率损耗等，是电力系统的稳态计算。潮流计算是对电力系统正常运行状况的分析和计算，即电力系统中的电压、电流、功率的计算，即潮流计算;潮流计算方法很多：高斯塞德尔法、牛顿拉夫逊法、-Q分解法、直流潮流法,以及由高斯塞德尔法、牛顿拉夫逊法演变的各种潮流计算方法。潮流计算可以用传统的手工方式进行，也可以计算机为工具通
11、过软件完成。两种方法各有优缺点。前者物理概念清晰，可用来计算一些接线较简单的电力网,但若将其用于接线复杂的电力网则计算量过大，难于保证计算准确性。后者从数学上看可归结为用数值方法解非线性代数方程，数学逻辑简单完整,借助计算机可快速精确地完成计算,但其缺点是物理概念不明显,物理规律被埋没在循环往复的数值求解过程中。潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。对规划中的电力系统,通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求:对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在
12、允许的范围以内，系统中各元件（线路、变压器等）是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。因此潮流计算的目的是：为电力系统规划设计提供接线、电气设备选择和导线截面选择的依据。 提供电力线运行方式和制定检修计划的依据。 提供继电保护、自动装置设计和整定计算的依据。为调压计算、经济运行计算、短路和稳定计算提供必要的数据。潮流计算的研究是从20世纪50年代随着电网的产生而开始的,涌现出各种算法都是围绕着算法的可靠性或收敛性;对计算机内存量的要求和计算速度;计算的方便性和灵活性等。1. 潮流计算的意义潮流计算是研究电力系统稳态运行情况的一种基本电气计算,常规潮流计算的任务是根据给定
13、的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统,通过潮流计算,可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。具体表现在以下方面：(1)在电网规划阶段，通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案，满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。
14、（)在编制年运行方式时,在预计负荷增长及新设备投运基础上，选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。(）正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制，指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。总结为在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时,为了实时监控电力系统的运行状态,也需要进行大量而快速的潮
15、流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时,采用离线潮流计算;在电力系统运行状态的实时监控中，则采用在线潮流计算。1.3电力系统稳态分析潮流计算总结潮流计算针对电力革统各正常运行方式,而静态安全分析则要研究各种运行方式下个别系统元件退出运行后系统的状况。其目的是校验系统是否能安全运行,即是否有过负荷的元件或电压过低的母线等。原则上讲,静态安全分析也可U用潮流计算来代替。但是一般静态安全分析需要校验的状态数非常多，用严格的潮流计算来分析这些状态往往计算量过大,因此不得不寻求一些特殊的算法以满足要求。利用电子数字计算机进行电力系
16、统潮流计算从20 世纪0 年代中期就己开始，此后，潮流计算曾采用了各种不同的方法,这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的，对潮流计算的要求可以归纳为下面几点:(1）计算方法的可靠性或收敛性 (2)对计算速度和内存量的要求 （3)计算的方便性和灵活性一般情况下，采用P-Q分解法计算时要求的迭代次数较采用牛顿-拉夫逊法时多，但每次迭代所需的时间则较采用牛顿-拉夫逊法时少，以致总的计算速度仍是PQ分解法快。.4 MATAB的概述目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1)计算速度快（2)内存需
17、要少（3)计算结果有良好的可靠性和可信性（）适应性好,亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强(5）简单。 MATLAB是一种交互式、面向对象的程序设计语言，广泛的应用于工业界和学术界,主要用于矩阵运算,同时在数值分析、自动控制模拟、MAT程序设计语言结构完整,且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效地解决工业计算的问题，特别是关于矩阵和矢量的计算。MATB与C语言和FORRN语言相比更容易掌握。通过ATLA语言，可以用类似数学公式的方式来编写算法,大大降低了程序所需的难度并节省了时间，从而可以把主要的精力集中在算法的构思而不是编程上。另外
18、,MATLAB提供了一种特殊的工具:工具箱（TOOBXES）.这些工具箱主要包括：信号处理（SINAL PROESSIN)、控制系统(CNTOL SYSTEMS）、神经网络（EUR NEWOS)、模糊逻辑(FUZZY LOGIC)、小波(AVELETS)和模拟（MULTION）等等。不同领域、不同层次的用户通过相应工具的学习和应用,可以方便地进行计算、分析及设计工作。 AAB设计中,原始数据的填写格式是很关键的一个环节,它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计,主要应从使用的角度出发,原则是简单明了，便于修改。第二章 理论计算21 Q法潮流计算的基本步骤 运用P-Q分
19、解法计算潮流分布时的基本步骤是：（1） 形成系数矩阵B、，并求其逆阵。（2） 设各节点电压的初值(i1,2，n，is)和(i=，2，m，is)。（3） 按式=-计算有功功率的不平衡量,从而求出(i=，2，ns)。（4） 解修正方程式 （）=BU，求各节点电压相位角的变量（=1，2，,n，s)。（5） 求各节点电压相位角的新值=(i=,2，n，i)。（6） 按式=- 计算无功功率的不平衡量,从而求出（i1,2，,m,is）。（7） 解修正方程式()=BU（i=,2,，m，is）。（8） 求各节点电压大小的新值+（i=1,2，i）。（9） 运用各节点电压的新值自第三步开始进入下一次迭代。（10）
20、计算平衡节点功率和线路功率。2.2 PQ分解法潮流计算流程图 启动输入原始数据形成节点导纳矩阵形成BB并求其逆阵计算 计算电压新值新值是否满足要求满足 停止置迭代次数k=0不满足k+12.3两机五节点网络潮流计算系统接线图解: 1，.形成系数矩阵B、B并求他们的逆阵由于节点1为平衡节点,其他的节点均为PQ节点，系数矩阵，阶数相同。又应对该等值网络，不存在除去与有功功率和电压相位或无功功率和电压大小关系较小的因素的可能性,这两个矩阵、B完全相同。他们就由导纳矩阵的虚部部分中除第一行第一列外的各个元素所组成,即6.5-18.5j -+15j -1.2+375j 0 -55j 10.834325j
21、-1+5 1.6+5j -.5+7.5 125+3.75j -1.67+5j .17-38.75j -10+3j 0 0 -16+5j 10+30j 12.917-8.5j 125+.5 2.55j 0 -.25+3.75j 3.75-15jB=B= 2.00 5.000 000 .500 500 -38.70 30.00 0 .000 3.000 -38.750 3.70 .50 0 375 11.250 由此可见,网络的节点导纳矩阵虽为奇异矩阵,但它的虚数部分的子阵B或则是非奇异矩阵,可以求逆,其逆阵为= -005619 -0.19 .44762 -0.0381 -0.04195 -099
22、48 .0769 -0.057143 -0.04762 -0819 -0.0514 -0.6507 -0.05238 -0。05713 -0.065079 -.4503 2计算各节点有功功率的不平衡量 取=106,=0;=.0;=0,按下式计算各节点有功功率不平衡量 =-=0201.01.06（-5.000cos+15.0in0)1.01.0(1.3cos032.500sn0)1.00(-1.667o05.000in0）1.01（1.co+5.0in）0.0（-2.0cos+7.500si0)0.50000相似地可得 -0.350；-0.0000；-0.60000 3计算各节点电压的相位角（弧
23、度） 由方程式 （)=U可得0.03652;=0.554；=0910；.108;所以=.392； =0.08554； =0.09180； 0.185； 4.计算各节点无功功率不平衡量按下式计算各节点无功功率不平衡量 =-=.01.01.065.0sin(369520）.000cos（-0.09520).1.010.8sin02.5cos1.010-1.667in（-.36952+0.8524）5.00os（-0.3692.085524)101.7sin(.0952+0,01810）5.000cs（0,036920.0981)0.0-2.00(.036920.108571）7.50co(-.03
24、6952+.108571)=.211930相似地可得 =0.077279; =0.11947; .1959； .计算各节点电压的大小 由方程式 ()U可得=.039528； =.008050; 0.0638； =0.0072;所以=.03928； =00080； =0.0038； =0.0002;求得各节点电压的新值后,就开始第二次迭代。每次迭代所得适于表1表。由表1可见经6次迭代就可满足的要求。表1 迭代过程中各节点功率的不平衡量-0.000001.1930-0.370-.077790.40000.19470.600000-.3195991-0.10-01367.42.01800.0728.
25、24849.12089904021820.432230.0148750005890.0196-0.00785600026640.1258-0042653.004833-006580.0500.00050.0007200029.03760.04640.00037.000184-0.000720,0002-.00096-.0003.001500.000515-000006-.0000.00000.0003.001.000040.0006000000660.00000000-0000010.0000-0000.0000-0000002000000表2 迭代过程中各节点电压的修正量k000369520
26、.03928-.8524.0000.09180063860.105710.00071-.1063-0.0034060.00128.07700024.0009410.0045970.0016120.10540385-00300.004-00085-.0002-.484-0.00014930.00120-.0000430.000.0000-00028.00000.00490.0001540.00030005-0.000002-.00000-003-0.00001-0.00005-.00-0.00002-.0000000000.00000.000000.000000.000010.00000表3 迭
27、代过程中各节点电压k00.000000.0000000.00000.000000.0000001.00000.001.000010.0369521.039528-.05241080-0.08100060，0857.000072-047061.03612-0.083961.008-0,0893481.0073270,10397410016903-0.59206507-.8392100-.896331.00725-0104581.001544-0.0402103643-0.83903.00850-.0860510723-0.10410100155-.04606136468-0.03961.000-
28、.096081.00252-0.104151.005540.04600103648-0,083906008750-0.08901.00722-0.1014.00554 =.036468; =1.000 =1.007252; =1.0154 计算平衡节点功率和线路功率 迭代收敛后,就可以计算平衡节点功率和线路功率。结果如下 =1.29812j0.2472各线路功率如表4表4 各线路功率j i134508905+j0138620.40867+j0.105802-087509-j.09385.24688.81450.27919+0.00600.54870+0.1323-0.39596-j0.0674
29、-4308-j.070130.8979-j0.012204-024598-j0.06645-018726+j0.0137806332j0.0036750599-0.097684-0308j00021计算完毕。第三章 程序设计3. 设计流程2程序设计针对两机五节点网络潮流计算基于MATLAB的程序如下：C 清除当前屏幕上显示的内容,但并不清楚工作空间的数据Clar %清除工作空间内的所有变量disp（节点总数为:)； %显示内容N=5 dip(平衡节点为:）;1=inut（请输入导纳矩阵：Y=);pr=t(请输入误差精度:pr=);Greal(Y)；mag（Y）; real实部，iae虚部B1=
30、B（:5, 2：5)； %去掉平衡节点B2=B1; b=inv(B1); %取逆=inv（B）;ure=npu（请输入平衡节点1电压初值的实部，e=);dm=input(请输入平衡节点电压初值的虚部,I=);fo m=2：5dt(m)0; %电压虚部初值u（m)1； 电压实部初值edor m=2:5p(m)=inpt(请输入PQ节点的有功功率：);（)npu(请输入P节点的无功功率：=);edk=0;wa1; %迭代次数,hle ha（1）p;u（1)ure；dt(1)dim； %平衡节点电压实部fr m=2：5 %平衡节点电压虚部fr =:5p()u(m)u(n）*(m,n)os(d(m)-
31、t()+B(m,)*si(t(m)-t(); %enddisp(各节点有功率不平衡量为：) dp(m)=p(m)-um(t) =-endr m=：u1=ag（u); 对角矩阵nv（u(2:5,2：）); %取逆u3=ipd(ro90（u（2：5）); %逆时针90上下对称u4=3(:5-,：）; p1flid(rot90（dp(2：5)； %（)=Up2=dp(：,:)；t=(-b1*dp2u4)*u2；tt=(:,:）；tt=fipd(ot90(tt)dp（各节点电压的相位角dt为:) d(）=dt(m)dtt(m-)edfor m=：5for =1:5 qt()=u()*u(n)(G(m,
32、n)*sin(d（)-t(）-B(m，n)*cos(d()dt();end %p(各节点无功率不平衡量为:)q(m)=()-sum（qt) % =-endfr m=:q=lip(rot0(q（:5）); q2=qq1(:5-1,:)；ut=-b2*qq2/u4;dis( 各节点电压的修正量：)utdp(各节点电压的大小为:） u（m）=u(m)+ut(m1) edr n=1:5 U(n)=u(n)*(cos（d（)+*si(dt(n); %电压endfo =：5 (n)Y(1,n)*U（n)；endisp（各节点功率s为:)1U(1)sum（conj(I) %电流fo=1：5for =:5 S
33、(m，）U(m)*（conj(U（m)-coj(（）)）*cnj（-Y(，n）); %功率endendsp(各节点功率sij为：）Sucha=x(abs(ut)）=k+1 %迭代+nd.3 程序运行结果节点总数为:N = 5平衡节点为：as = 请输入导纳矩阵:=6.25-18.75i -5+1 -.25+3.5i 0 0 -+15i .34-3. -.667+5i -1675i 2.i-1.5.5i -1.667+5i 2.917-3.75i -0+0i 0 0 -1667+5 -10+30 12913875i -1.25+3.7i0 -+7.5i 0 -.25+3.75i .75-1125
34、请输入误差精度：pr=0.00001请输入平衡节点1电压初值的实部，Re.06请输入平衡节点电压初值的虚部，Im=0请输入P节点的有功功率：p=0.2请输入PQ节点的无功功率：=0.请输入节点的有功功率：p=45请输入PQ节点的无功功率：q-0.15请输入PQ节点的有功功率:p=-04请输入PQ节点的无功功率:=0.05请输入Q节点的有功功率：=-.6请输入PQ节点的无功功率:q=-.1各节点有功率不平衡量为:d= 0.000各节点有功率不平衡量为: 0.00 -0.750各节点有功率不平衡量为:dp= 0.500 -0.3750 -0.000各节点有功率不平衡量为:dp 0 0.5000 0
35、3750 -0.400 -0.000dt = -0030 0.085 -0.0918 -0.1086各节点电压的相位角d为：dt = -.0370 0 0 0t = 0.0370 .855 -0.918 -0.1086各节点电压的相位角dt为：d = 0 070 0.05 0dtt = -0.0370 0.0855 0.0918 -0.108各节点电压的相位角为：t = 0 -0.37 0.055 0.018 dt = -0.070 0.05 -0.918 0.186各节点电压的相位角dt为：d 0 -.70 -0.05 -0.098 -0.1086各节点无功率不平衡量为:qq = 0 .2119各节点无功率不平衡量为：qq= 1.2119 -0.0773各节点无功率不平衡量为:qq= 0 1.219 -0.073 -0.19各节点无功率不平衡量为:qq = 0 1219 -0.7 -0.199 -0.316各节点
