基于MATLAB的电力系统复杂潮流分析课程设计报告书.doc
- 文档编号:30578683
- 上传时间:2023-08-17
- 格式:DOC
- 页数:60
- 大小:515.50KB
基于MATLAB的电力系统复杂潮流分析课程设计报告书.doc
《基于MATLAB的电力系统复杂潮流分析课程设计报告书.doc》由会员分享,可在线阅读,更多相关《基于MATLAB的电力系统复杂潮流分析课程设计报告书.doc(60页珍藏版)》请在冰豆网上搜索。
科技大学
电力系统稳态分析课程设计
题目:
基于MATLAB的电力系统 复杂潮流分析
学生:
建峰
学号:
1167130207
专业:
电气工程及其自动化
班级:
电气2011—2班
指导教师:
景霞
摘要
电力系统潮流计算是电网分析的基础应用,是对复杂电力系统正常和故障条件下稳态运行状态的计算。
给定电力系统的网络结构、参数和决定电力系统运行状况的边界条件,确定电力系统运行的方法之一是朝流计算。
MATLAB是一种交互式、面向对象的程序设计语言,广泛应用于工业界与学术界,主要用于矩阵运算.采用迭代法,通过建立矩阵的修正方程来依次迭代,逐步逼近真值来计算出电力网的电压,功率分布。
PQ分解法是极坐标形式牛顿-拉潮流计算的一种简化计算方法,。
P—Q分解法通过对电力系统具体特点的分析,对牛顿法修正方程式的雅可比矩阵进行了有效的简化和改进。
由于这些简化只涉及修正方程式的系数矩阵,并未改变节点功率平衡方程和收敛判据,因不会降低计算结果的精度。
用手算和计算机算法对其进行设计。
使用MATLAB软件进行编程,在很大程度上节省了存,减少了计算量。
通过对本题计算我们了解了一些工程计算和解决工程问题的方法。
关键词:
潮流计算,PQ分解法,MATLAB
ElectricalpowersystemcomplextidalcurrentanalysisbasedonMATLAB
PowerFlowAnalysisGridcomputingisthebasisofapplications,thecomplexpowersystemundernormalandfaultconditionsforthecalculationofsteadystateoperation.Giventhepowersystemnetworkstructure,parametersanddecisionsoperationofthepowersystemboundaryconditions,todeterminethemethodofoperationofthepowersystemisoneofNorthKoreaflowcalculation.
MATLABisaninteractive,object-orientedprogramminglanguage,widelyusedinindustryandacademia,mainlyformatrixcalculation.Usingiteration,theamendmentthroughtheestablishmentofmatrixiterativeequationtoturn,graduallymovingtowardsatruevaluetocalculatethevoltageelectricitygrid,powerdistribution.
PQdecompositionmethodistheformofpolarcoordinatesNewton-thewideningtrendofasimplifiedcalculationmethod.P-Qdecompositionmethodadoptedbythespecificcharacteristicsofthepowersystemanalysis,Newton’sLawoftheJacobianmatrixformulahaseffectivelysimplifiedandimproved.Asaresultofthesesimplifiedformulathatinvolvesonlythecoefficientmatrix,thebalanceofpowerhasnotchangednodeequationsandtheconvergencecriterion,becausetheresultswillnotreducetheaccuracy.
UseMATLABsoftwareprogramming,savingmemorytoalargeextent,reducetheamountofcomputation.Bythiscalculationweunderstandthatanumberofengineeringcalculationandsolveengineeringproblems
Keywords:
Thetrend,thePQdecompositionmethod,MATLAB
目录
科技大学课程设计任务书 5
第一章引言 7
1.1研究背景及意义 7
1.2潮流计算的意义 8
1.3电力系统稳态分析潮流计算总结 9
1.4MATLAB的概述 9
第二章理论计算 11
2.1P-Q法潮流计算的基本步骤 11
2.2PQ分解法潮流计算流程图 12
2.3两机五节点网络潮流计算 13
第三章程序设计 18
3.1设计流程 18
3.2程序设计 18
3.3程序运行结果 22
第四章设计感想 44
参考文献 45
科技大学课程设计任务书
课程名称
电力系统稳态分析课程设计
设计题目
两机五节点网络潮流计算—P-Q法
指导教师
景霞
时间
1周
一、教学要求
电力系统稳态分析课程设计以设计和优化电力系统的潮流分析为重点,提高学生综合能力为目标,尽可能结合实际工程进行。
设计容的安排要充分考虑学校现有的设备,设计时间及工程实际需要,并使学生初步学会运用所学知识解决工程中的实际问题。
二、设计资料及参数
(一)设计原始资料
1、待设计电气设备系统图
2、电力系统网络各元件参数
3、电力系统电气元件的使用规
4、电力工程电气设计手册
(二)设计参考资料
1、《电力系统稳态分析》,珩,中国电力,2007,第三版
2、《电力系统分析》,祯祥,大学,2005,第三版
3、《电力系统分析课程实际设计与综合实验》,祝书萍,中国电力,2007,第一版
三、设计要求及成果
1.根据给定的参数或工程具体要求,收集和查阅资料;学习相关软件(软件自选)。
2.在给定的电力网络上画出等值电路图。
3.运用计算机进行潮流计算。
4.编写设计说明书。
基本要求:
1.编写潮流计算程序;
2.在计算机上调试通过(?
);
3.运行程序并计算出正确结果(?
);
4.写出课程设计报告(包括以下容)(1份)
(1)程序框图;
(2)源程序;(3)符号说明表;(4)算例及计算结果
5.编写计算说明书(1份)。
四、进度安排
根据给定的参数或工程具体要求,收集和查阅资料(半天)
学习软件(MATLAB或C语言等)(一天半)
编程计算复杂系统潮流计算(三天)
编写计算设计书(一天)
五、评分标准
课程设计成绩评定依据包括以下几点:
1)工作态度(占10%);
2)基本技能的掌握程度(占20%);
3)程序编写是否合理是否有运行结果(40%);
4)课程设计说明书编写水平(占30%)。
5)分为优、良、中、合格、不合格五个等级。
考核方式:
设计期间教师现场检查;评阅设计报告
第一章引言
1.1研究背景及意义
电力系统在运行时,在电源电势激励作用下,电流或功率从电源通过系统各元件流入负荷,分布于电力网各处,称为潮流分布。
潮流计算是电力系统分析中的一种最基本的计算,即根据给定电力系统的网络结构和参数,在满足电力系统运行状态参数的边界条件情况下,确定电力系统稳态运行状态的最基本的方法。
它的任务是在给定的接线方式和运行条件下,确定系统的运行状态,如各母线上的电压(幅值和相角)、网络中的功率分布及功率损耗等,是电力系统的稳态计算。
潮流计算是对电力系统正常运行状况的分析和计算,即电力系统中的电压、电流、功率的计算,即潮流计算;潮流计算方法很多:
高斯—塞德尔法、牛顿—拉夫逊法、P-Q分解法、直流潮流法,以及由高斯—塞德尔法、牛顿—拉夫逊法演变的各种潮流计算方法。
潮流计算可以用传统的手工方式进行,也可以计算机为工具通过软件完成。
两种方法各有优缺点。
前者物理概念清晰,可用来计算一些接线较简单的电力网,但若将其用于接线复杂的电力网则计算量过大,难于保证计算准确性。
后者从数学上看可归结为用数值方法解非线性代数方程,数学逻辑简单完整,借助计算机可快速精确地完成计算,但其缺点是物理概念不明显,物理规律被埋没在循环往复的数值求解过程中。
潮流计算是电力系统非常重要的分析计算,用以研究系统规划和运行中提出的各种问题。
对规划中的电力系统,通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求:
对运行中的电力系统,通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全,系统中所有母线的电压是否在允许的围以,系统中各元件(线路、变压器等)是否会出现过负荷,以及可能出现过负荷时应事先采取哪些预防措施等。
因此潮流计算的目的是:
①为电力系统规划设计提供接线、电气设备选择和导线截面选择的依据。
②提供电力线运行方式和制定检修计划的依据。
③提供继电保护、自动装置设计和整定计算的依据。
④为调压计算、经济运行计算、短路和稳定计算提供必要的数据。
潮流计算的研究是从20世纪50年代随着电网的产生而开始的,涌现出各种算法都是围绕着算法的可靠性或收敛性;对计算机存量的要求和计算速度;计算的方便性和灵活性等。
1.2潮流计算的意义
潮流计算是研究电力系统稳态运行情况的一种基本电气计算,常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态,如各母线上的电压(幅值及相角)、网络中的功率分布以及功率损耗等。
潮流计算的结果是电力系统稳定计算和故障分析的基础。
通过潮流计算可以判断电网母线电压、支路电流和功率是否越限,如果有越限,就应采取措施,调整运行方式。
对于正在规划的电力系统,通过潮流计算,可以为选择电网供电方案和电气设备提供依据。
潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。
具体表现在以下方面:
(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。
(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。
(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。
(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。
总结为在电力系统运行方式和规划方案的研究中,都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。
同时,为了实时监控电力系统的运行状态,也需要进行大量而快速的潮流计算。
因此,潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。
在系统规划设计和安排系统的运行方式时,采用离线潮流计算;在电力系统运行状态的实时监控中,则采用在线潮流计算。
1.3电力系统稳态分析潮流计算总结
潮流计算针对电力革统各正常运行方式,而静态安全分析则要研究各种运行方式下个别系统元件退出运行后系统的状况。
其目的是校验系统是否能安全运行,即是否有过负荷的元件或电压过低的母线等。
原则上讲,静态安全分析也可U用潮流计算来代替。
但是一般静态安全分析需要校验的状态数非常多,用严格的潮流计算来分析这些状态往往计算量过大,因此不得不寻求一些特殊的算法以满足要求。
利用电子数字计算机进行电力系统潮流计算从20世纪50年代中期就己开始,此后,潮流计算曾采用了各种不同的方法,这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的,对潮流计算的要求可以归纳为下面几点:
(1)计算方法的可靠性或收敛性
(2)对计算速度和存量的要求
(3)计算的方便性和灵活性
一般情况下,采用P-Q分解法计算时要求的迭代次数较采用牛顿-拉夫逊法时多,但每次迭代所需的时间则较采用牛顿-拉夫逊法时少,以致总的计算速度仍是P-Q分解法快。
1.4MATLAB的概述
目前电子计算机已广泛应用于电力系统的分析计算,潮流计算是其基本应用软件之一。
现有很多潮流计算方法。
对潮流计算方法有五方面的要求:
(1)计算速度快
(2)存需要少(3)计算结果有良好的可靠性和可信性(4)适应性好,亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强(5)简单。
MATLAB是一种交互式、面向对象的程序设计语言,广泛的应用于工业界和学术界,主要用于矩阵运算,同时在数值分析、自动控制模拟、MATLAB程序设计语言结构完整,且具有优良的移植性,它的基本数据元素是不需要定义的数组。
它可以高效地解决工业计算的问题,特别是关于矩阵和矢量的计算。
MATLAB与C语言和FORTRAN语言相比更容易掌握。
通过MATLAB语言,可以用类似数学公式的方式来编写算法,大大降低了程序所需的难度并节省了时间,从而可以把主要的精力集中在算法的构思而不是编程上。
另外,MATLAB提供了一种特殊的工具:
工具箱(TOOLBOXES).这些工具箱主要包括:
信号处理(SIGNALPROCESSING)、控制系统(CONTROLSYSTEMS)、神经网络(NEURALNETWORKS)、模糊逻辑(FUZZYLOGIC)、小波(WAVELETS)和模拟(SIMULATION)等等。
不同领域、不同层次的用户通过相应工具的学习和应用,可以方便地进行计算、分析及设计工作。
MATLAB设计中,原始数据的填写格式是很关键的一个环节,它与程序使用的方便性和灵活性有着直接的关系。
原始数据输入格式的设计,主要应从使用的角度出发,原则是简单明了,便于修改。
第二章理论计算
2.1P-Q法潮流计算的基本步骤
运用P-Q分解法计算潮流分布时的基本步骤是:
(1)形成系数矩阵Bˊ、B〞,并求其逆阵。
(2)设各节点电压的初值(i=1,2,…,n,i≠s)和(i=1,2,…,m,i≠s)。
(3)按式=-计算有功功率的不平衡量,从而求出(i=1,2,…,n≠s)。
(4)解修正方程式()=﹣BˊU,求各节点电压相位角的变量(i=1,2,…,n,i≠s)。
(5)求各节点电压相位角的新值=+(i=1,2,…,n,i≠s)。
(6)按式=-计算无功功率的不平衡量,从而求出(i=1,2,…,m,i≠s)。
(7)解修正方程式()=﹣B〞U(i=1,2,…,m,i≠s)。
(8)求各节点电压大小的新值=+(i=1,2,…,m,i≠s)。
(9)运用各节点电压的新值自第三步开始进入下一次迭代。
(10)计算平衡节点功率和线路功率。
2.2PQ分解法潮流计算流程图
启动
输入原始数据
形成节点导纳矩阵
形成BˊB"并求其逆阵
计算计算电压新值
新值是否满足要求
满
足
停止
置迭代次数k=0
不满足
k+1
2.3两机五节点网络潮流计算
系统接线图
解:
1,.形成系数矩阵Bˊ、B〞并求他们的逆阵
由于节点1为平衡节点,其他的节点均为PQ节点,系数矩阵,阶数相同。
又应对该等值网络,不存在除去与有功功率和电压相位或无功功率和电压大小关系较小的因素的可能性,这两个矩阵Bˊ、B〞完全相同。
他们就由导纳矩阵的虚部部分中除第一行第一列外的各个元素所组成,即
Y=[6.25-18.75j-5+15j-1.25+3.75j00
-5+15j10.834-32.5j-1.667+5j-1.667+5j-2.5+7.5j
-1.25+3.75j-1.667+5j12.917-38.75j-10+30j0
0-1.667+5j-10+30j12.917-38.75j-1.25+3.75j
0-2.5+7.5j0-1.25+3.75j3.75-11.25j]
Bˊ=B〞=
[-32.5005.0005.0007.500
5.000-38.75030.0000
5.00030.000-38.7503.750
7.50003.750-11.250]
由此可见,网络的节点导纳矩阵虽为奇异矩阵,但它的虚数部分的子阵Bˊ或B〞则是非奇异矩阵,可以求逆,其逆阵为
==
[-0.056190-0.041905-0.44762-0.052381
-0.041905-0.099048-0.087619-0.057143
-0.044762-0.087619-0.105714-0.065079
-0.052381-0。
057143-0.065079-0.145503]
2.计算各节点有功功率的不平衡量
取=1.06,=0;====1.0;====0,按下式计算各节点有功功率不平衡量
=-
=0.20—1.0×1.06(-5.000cos0+15.000sin0)—1.0×1.0(10.384cos0—32.500sin0)—1.0×1.0(-1.667cos0+5.000sin0)—1.0×1.(-1.667cos0+5.000sin0)—1.0×1.0(-2.500cos0+7.500sin0)=0.50000
相似地可得
=-0.375000;=-0.400000;=-0.600000
3.计算各节点电压的相位角(弧度)
由方程式—()=U
可得
=—0.036952;=—0.085524;=—0.091810;=—0.108571;
所以
=—0.036952;=—0.085524;=—0.091810;=—0.108571;
4.计算各节点无功功率不平衡量
按下式计算各节点无功功率不平衡量
=-
=0.20—1.0×1.06[—5.00sin(-0.036952—0)—15.000cos(-0.036952—0)]—1.0×1.0[10.834sin0+32.5cos0]—1.0×1.0[-1.667sin(-0.036952+0.085524)—5.000cos(-0.036952+0.085524)]—1.0×1.0[-1.667sin(—0.036952+0,091810)—5.000cos(-0,036952+0.091810)]—1.0×1.0[-2.500sin(—0.036952+0.108571)—7.500cos(-0.036952+0.108571)]=1.211930
相似地可得
=—0.077279;=—0.191947;=—0.319599;
5.计算各节点电压的大小
由方程式
—()=U
可得
=1.039528;=0.008050;=0.006386;=0.000072;
所以
=1.039528;=0.008050;=0.006386;=0.000072;
求得各节点电压的新值后,就开始第二次迭代。
每次迭代所得适于表1~表3。
由表1可见经6次迭代就可满足ε≤的要求。
表1迭代过程中各节点功率的不平衡量
K
0
-0.500000
1.211930
-0.375000
-0.077279
0.400000
-0.191947
-0.600000
-0.319599
1
-0.411720
-0.136497
0.049802
0.018300
0.076248
0.024849
0.120899
0.040218
2
0.043223
0.014875
-0.005899
-0.001936
-0.007856
-0.002664
-0.012588
-0.004265
3
-0.004833
-0.001658
-0.000650
0.000215
0.000872
0.000295
0.001376
0.000466
4
0.000537
0.000184
-0.000072
-0,000024
-0.000096
-0.000033
-0.000150
-0.000051
5
-0.000060
-0.000021
0.000008
0.000003
0.000011
0.000004
0.000016
0.000006
6
0.000007
0.000000
-0.000001
0.000000
-0.000001
0.0000000
-0.000002
0.000000
表2迭代过程中各节点电压的修正量
k
0
-0.036952
0.039528
-0.085524
0.008050
-0.091810
0.006386
-0.108571
0.000072
1
-0.010063
-0.003406
0.001828
0.000757
0.002461
0.000941
0.004597
0.001618
2
0.001054
0.000385
-0.000230
-0.000064
-0.000285
-0.000082
-0.000484
-0.000149
3
-0.000120
-0.000043
0.000023
0.0
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 MATLAB 电力系统 复杂 潮流 分析 课程设计 报告书