电力系统课程设计潮流计算.docx
《电力系统课程设计潮流计算.docx》由会员分享,可在线阅读,更多相关《电力系统课程设计潮流计算.docx(36页珍藏版)》请在冰豆网上搜索。
电力系统课程设计潮流计算
电力系统课程设计潮流计算
一、概述
1.1设计目的与要求..................................................3
1.1.1设计目的.....................................................3
1.1.2设计要求.....................................................3
1.2设计题目.....................................................3
1.3设计内容.....................................................3
二电力系统潮流运算概述......................................4
2.1电力系统简介...................................................4
2.2潮流运算简介...................................................4
2.3潮流运算的意义及其进展........................................5
三潮流运算设计题目...........................................6
3.1潮流运算题目.................................................6
3.2对课题的分析及求解思路.......................................7
四潮流运算算法及手工运算....................................7
4.1极坐标下P-Q法的算法..........................................7
4.2节点电压方程..................................................8
4.3节点导纳矩阵..................................................9
4.4导纳矩阵在潮流运算............................................10
4.5潮流运算的手工运算............................................12
五Matlab概述..................................................13
5.1Matlab简介...................................................14
5.2Matlab的应用..................................................14
5.3矩阵的运算....................................................14
5.3.1与常数的运算..............................................14
5.3.2差不多数学运算...................................................14
5.3.3逻辑关系运算...................................................14
5.4Matlab中的一些命令.............................................15
六潮流运算流程图及源程序......................................18
6.1潮流运算流程图..................................................18
6.2潮流运算源程序图................................................19
6.3运行运算结果....................................................27
七总结...................................................29
八参考文献...............................................29
第一章系统概述
1.1设计目的与要求
1.1.1设计目的
1.把握电力系统潮流运算的差不多原理;
2.把握并能熟练运用一门运算机语言(MATLAB语言或C语言或C++语言);
3.采纳运算机语言对潮流运算进行运算机编程。
1.1.2设计要求
1.程序源代码;
2.给定题目的输入,输出文件;
3.程序说明;
4.给定系统的程序运算过程;
5.给定系统的手算过程(至少迭代2次)。
1.1设计题目
电力系统潮流运算(牛顿-拉夫逊法、P-Q分解法)
1.2设计内容
1.依照电力系统网络推导电力网络数学模型,写出节点导纳矩阵;
2.给予各节点电压变量(直角坐标系形式)初值后,求解不平稳量;
3.形成雅可比矩阵;
4.求解修正量后,重新修改初值,从2开始重新循环运算;
5.求解的电压变量达到所要求的精度时,再运算各支路功率分布、功率损耗和平稳节点功率;
6.上机编程调试;
7.运算分析给定系统潮流分析并与手工运算结果做比较分析;
8.书写课程设计说明书。
第二章电力系统潮流运算概述
2.1电力系统叙述
电力工业进展初期,电能是直截了当在用户邻近的发电站(或称发电厂)中生产的,各发电站孤立运行。
随着工农业生产和都市的进展,电能的需要量迅速增加,而热能资源和水能资源丰富的地区又往往远离用电比较集中的都市和工矿区,为了解决那个矛盾,就需要在动力资源丰富的地区建立大型发电站,然后将电能远距离输送给电力用户。
同时,为了提高供电的可靠性以及资源利用的综合经济性,又把许多分散的各种形式的发电站,通过送电线路和变电所联系起来。
这种由发电机、升压和降压变电所,送电线路以及用电设备有机连接起来的整体,即称为电力系统。
现代电力系统提出了“灵活交流输电和新型直流输电”的概念。
灵活交流输电技术是指运用固态电子器件与现代自动操纵技术对交流电网的电压、相位角、阻抗、功率以及电路的通断进行实时闭环操纵,从而提高高压输电线路的诉讼能力和电力系统的稳态水平。
新型直流输电技术是指应用现电力电子技术的最新成果,改善和简化变流站的造价等。
运营方式治理中,潮流是确定电网运行方式的差不多动身点:
在规划领域,需要进行潮流分析验证规划方案的合理性;在实时运行环境,调度员潮流提供了电网在预想操作预想下的电网的潮流分布以及校验运行的可靠性。
在电力系统调度运行的多个领域都涉及到电网潮流运算。
潮流是确定电力网咯运行状态的差不多因素,潮流问题是研究电力系统稳态问题的基础和前提。
2.2潮流运算简介
电力系统潮流运确实是研究电力系统稳态运行情形的一种运算,它依照给定的运行条件及系统接线情形确定整个电力系统各部分的运行状态:
各母线的电压。
各元件中流过的功率,系统的功率损耗等等。
在电力系统规划的设计和现有电力系统运行方式的研究中,都需要利用潮流运算来定量的分析比较供电方案或运行方式的合理性。
可靠性和经济性。
此外,电力系统的潮流运算也是运算机系统动态稳固和静态稳固的基础,因此潮流运确实是研究电力系统的一种和重要和基础的运算。
电力系统潮流运算也分为离线运算和在线运算两种,前者要紧用于系统规划设计和安排系统的运行方式,后者则用于正在运行系统的经常监视及实时操纵。
利用电子数字运算机进行潮流运算从50年代中期就差不多开始了。
在这20年内,潮流运算曾采纳了各种不同的方法,这些方法的进展要紧围绕着对潮流运算的一些差不多要求进行的,对潮流运算的要求能够归纳为以下几点:
(1)运算方法的可靠性或收敛性;
(2)对运算机内存量的要求;
(3)运算速度;
(4)运算的方便性和灵活性。
2.3潮流运算的意义及其进展
电力系统潮流运确实是电力系统分析中的一种最差不多的运算,是对复杂电力系统正常和故障条件下稳态运行状态的运算。
潮流运算的目标是求取电力系统在给定运行状态的运算,即节点电压和功率分布,用以检查系统各元件是否过负荷。
各点电压是否满足要求,功率的分布和分配是否合理以及功率损耗等。
对现有的电力系统的运行和扩建,对新的电力系统进行规划设计以及对电力系统进行静态和稳态分析差不多上以潮流运算为基础。
潮流运算结果可用如电力系统稳态研究,安全估量或最优潮流等对潮流运算的模型和方法有直截了当阻碍。
实际电力系统的潮流技术那要紧采纳牛顿—拉夫逊法。
运行方式治理中,潮流是确定电网运行方式的差不多动身点;在规划领域,需要进行潮流分析验证规划方案的合理性;在实时运行环境,调度员潮流提供了多个在预想操作情形下电网的潮流分布以及校验运行可靠性。
在电力系统调度运行的多个领域问题是研究电力系统稳态问题的基础和前提。
在用数字解算运算机解电力系统潮流问题的开始时期,普遍采取以节点导纳矩阵为基础的逐次代入法。
那个方法的原理比较简单,要求的数字运算机内存量比较差下,适应50年代电子运算机制造水平和当时电力系统理论水平,但它的收敛性较差,当系统规模变大时,迭代次数急剧上升,在运算中往往显现迭代不收敛的情形。
这就迫使电力系统的运算人员转向以阻抗矩阵为基础的逐次代入法。
阻抗法改善了系统潮流运算问题的收敛性,解决了导纳无法求解的一些系统的潮流运算,在60年代获得了广泛的应用,阻抗法德要紧缺点是占用运算机内存大,每次迭代的运算量大。
当系统不断扩大时,这些缺点就更加突出,为了克服这些缺点,60年代中期进展了以阻抗矩阵为基础的分块阻抗法。
那个方法把一个大系统分割为几个小的地区系统,在运算机内只需要储备各个地区系统的阻抗矩阵及它们之间联络的阻抗,如此不仅大幅度的节约了内存容量,同时也提高了运算速度。
克服阻抗法缺点是另一个途径是采纳牛顿-拉夫逊法。
这是数学中解决非线性方程式的典型方法,有较好的收敛性。
在解决电力系统潮流运算问题时,是以导纳矩阵为基础的,因此,只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性,就能够大大提高牛顿法潮流程序的效率。
自从60年代中期,牛顿法中利用了最佳顺序消去法以后,牛顿法在收敛性。
内存要求。
速度方面都超过了阻抗法,成为了60年代末期以后广泛采纳的优秀方法。
第三章潮流运算设计题目
3.1潮流运算课题
题目:
在图1所示的简单电力系统中,系统中节点1、2为
节点,节点3为
节点,节点4为平稳节点,已给定
,
,
,
,
,
,网络各元件参数的标幺值如表2所示,给定电压的初始值如表2所示,收敛系数
。
试求:
图3-1简单电力系统
表3-1网络各元件参数的标幺值
支路
电阻
电抗
输电线路
变压器变比k
1—2
0.03
0.09
0.02
—
1—3
0.02
0.05
0.02
—
2—3
0.04
0.08
—
—
2—4
0.0
0.05
—
0.9625
3—4
0.03
0.07
—
—
表2各节点电压(初值)标幺值参数
节点i
1
2
3
4
1.00+j0.0
1.0+j0.0
1.0+j0.0
1.05+j0.0
3.2对课题的分析及求解思路
此电力系统是一个4节点,5支路的电力网络。
综合比较牛顿拉夫逊法(直角坐标、极坐标)、PQ分解法等多种求解方法的特点,最后确定采纳牛顿拉夫逊法(极坐标)。
因为此方法所需解的方程组最少。
第四章潮流运算算法及手工运算
4.1极坐标下P-Q法的算法
4.1.1节点导纳矩阵Y
依照题目提供的各节点的参数,求得节点导纳矩阵
=
4.1.2简化雅可比矩阵B/和B//
通过上一步的导纳矩阵,形成有功迭代和无功迭代的简化雅可比矩阵B/和B//
对雅可比矩阵进行三角分解,形成因子表,为后面进行修正方程运算作好预备。
4.1.3修正和迭代
第一步,给定PQ节点初值和各节点电压相角初值。
第二步,作第一次有功迭代,按公式运算节点有功功率不平稳量。
第三步,做第一次无功迭代,按公式运算无功功率不平稳量,运算时电压相角最新的修正值。
解修正方程式,可得各节点电压幅值的修正量。
第四步,第一轮有功迭代和无功迭代便做完了。
第五步,按公式运算平稳节点功率。
直到节点不平稳功率下降到10-5以下,迭代便能够终止。
4.2潮流运算算法
本题采纳了题目要求的牛顿-拉夫逊潮流运算的方法。
牛顿-拉夫逊法潮流运算的公式。
把牛顿法用于潮流运算,采纳极坐标形式表示的如式(1-3)所示的形式。
其中电压和支路导纳可表示为:
将上述表示式(1-2)代入(1-1)式的右端,展开并分出实部和虚部,便得:
(1-3)
按照以上的分类,PQ节点的输出有功功率和无功功率是给定的,则第i节点的给定功率设为
和
(称为注入功率)。
假定系统中的第1、2、…、m节点为PQ节点,对其中每一个节点的N-R法表达式F(x)=0[如
、
、
]形式有些下列方程:
(1-4)
=(1、2、…、m)
PV节点的有功功率和节点电压幅值是给定的。
假定系统中的第m+1、m+2、…、n-1节点为PV节点,则对其中每一PV节点能够列写方程:
(1-5)
=(m+1、m+2、…、n-1)
(6)形成雅可比矩阵。
N-R法的思想是
;本例
;对F(x)求偏导的式(1-6)、式(1-7),即式(1-4)、式(1-5)中的
、
、
是多维变量的函数,对多维变量求偏导(
、
、
、
、
、
、
、…),并以矩阵的形式表达称为雅可比矩阵。
当j=i时,对角元素为
(1-6)
当
时,矩阵非对角元素为:
(1-7)
由上式不难看出,雅可比矩阵有以下特点。
1雅可比矩阵中的诸元素差不多上节点电压的函数,因此在迭代过程中,它们将随着节点电压的变化而不断的变化。
2雅可比矩阵具有结构对称性,数据不对称。
如非对角
,
,
。
3由式(1-7)能够看出,当导纳矩阵中非对角元素
为零时,。
雅可比矩阵中相应的元素也为零,即矩阵是专门稀疏的。
因此,修正方程的求解同样能够应用稀疏矩阵的求解技巧。
正是由于这一点才使N-R法获得广泛的应用。
4.3手工运算
4.3.1节点导纳矩阵
求得节点导纳矩阵Y
=
各节点的导纳值如下:
Y11=10.2299-27.2214i
Y12=-3.3333+10.0000i
Y13=-6.8966+17.2414i
Y14=0
Y21=-3.3333+10.0000i
Y22=8.3333-38.5181i
Y23=-5.0000+10.0000i
Y24=0+19.2500i
Y31=-6.8966+17.2414i
Y32=-5.0000+10.0000i
Y33=17.0690-39.3003i
Y34=-5.1724+12.0690i
Y41=0
Y42=0+19.2500i
Y43=-5.1724+12.0690i
Y44=5.1724-32.0690i
4.3.2简化雅可比矩阵
迭代中的雅克比矩阵:
-27.586210.000017.5862-10.09203.3333
10.0000-40.412510.20003.3333-8.2333
17.586210.2000-40.88817.03455.1000
10.3678-3.3333-7.0345-26.856610.0000
-3.33338.4333-5.100010.0000-36.6237
-28.388610.653017.7355-10.15633.4957
10.6198-42.338910.55743.5953-8.8273
17.989210.7970-40.88816.82785.0990
10.9577-3.4957-7.4620-27.795210.6530
-3.59539.4444-5.578110.6198-42.1163
-28.342910.619017.7239-10.13863.4922
10.5905-42.249110.54163.5776-8.7973
17.970010.7682-40.88816.83125.1009
10.9386-3.4922-7.4464-27.742910.6190
-3.57769.3973-5.554010.5905-41.8495
-28.342610.618917.7238-10.13853.4921
10.5903-42.248610.54143.5777-8.7971
17.970010.7682-40.88816.83105.1007
10.9385-3.4921-7.4465-27.742610.6189
-3.57779.3971-5.554110.5903-41.8486
-28.342710.618917.7238-10.13853.4921
10.5903-42.248610.54143.5777-8.7971
17.970010.7682-40.88816.83105.1007
10.9385-3.4921-7.4465-27.742710.6189
-3.57779.3971-5.554110.5903-41.8486
4.3.3修正、迭代
给定PQ节点初值和各节点电压相角初值
V1=1.0,V2(0)=V3(0)=1.0,V4=1.05
δ2(0)=δ3(0)=0,δ4(0)=0
1作第一次有功迭代,按公式运算节点有功功率不平稳量
迭代中的△P:
-0.2621-0.20000.3156
0.00070.0086-0.0167
-0.00000.00000.0008
0.0000-0.0000-0.0000
0.0000-0.00000.0000
2做第一次无功迭代,按公式运算无功功率不平稳量,运算时电压相角最新的修正值。
迭代中的△Q:
0.06481.6944
-0.0033-0.0887
-0.0000-0.0002
-0.0000-0.0000
-0.0000-0.0000
解修正方程式,可得各节点电压幅值的修正量为
迭代中电压的模:
1.01591.0470
1.01501.0448
1.01501.0448
1.01501.0448
1.01501.0448
到那个地点为止,第一轮有功迭代和无功迭代便做完了。
3按公式运算平稳节点功率,得:
P1+jQ1=0.3159+1.3621i
通过四轮迭代,节点不平稳功率也下降到10-5以下,迭代到此终止。
4.4输出功率的手工运算
全线路各个点的功率分配如下:
0-0.0582-0.3378i-0.3418+0.0172i0
0.0612-0.4727i0-0.0958-0.4238i-0.2654-0.9014i
0.3441-0.0424i0.1013-0.3746i0-0.0455-0.4275i
00.2654+0.9332i0.0505+1.2558i0
第五章Matlab概述
5.1Matlab简介
MATLAB是由美国mathworks公司公布的要紧面对科学运算、可视化以及交互式程序设计的高科技运算环境。
它将数值分析、矩阵运算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值运算的众多科学领域提供了一种全面的解决方案,并在专门大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学运算软件的先进水平
MATLAB和Mathematica、Maple并称为三大数学软件。
它在数学类科技应用软件中在数值运算方面首屈一指。
MATLAB能够进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等,要紧应用于工程运算、操纵设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。
5.2Matlab的应用
MATLAB的差不多数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的情况简捷得多,同时mathwork也吸取了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。
在新的版本中也加入了对C,FORTRAN,C++,JAVA的支持。
能够直截了当调用,用户也能够将自己编写的有用程序导入到MATLAB函数库中方便自己以后调用,此外许多的MATLAB爱好者都编写了一些经典的程序,用户能够直截了当进行下载就能够用。
MALAB产品族能够用来进行以下各种工作:
●数值分析
●数值和符号运算
●工程与科学绘图
●操纵系统的设计与仿真
●数字图像处理技术
●数字信号处理技术
●通讯系统设计与仿真
●财务与金融工程
MATLAB的应用范畴专门广,包括信号和图像处理、通讯、操纵系统设计、测试和测量、财务建模和分析以及运算生物学等众多应用领域。
附加的工具箱(单独提供的专用MATLAB函数集)扩展了MATLAB环境,以解决这些应用领域内特定类型的问题。
5.3与常数的运算
常数与矩阵的运算即是同该矩阵的每一元素进行运算。
但需注意进行数除时,常数通常只能做除量。
差不多函数运算中,矩阵的函数运确实是矩阵预算中最有用的部分,常用的要紧有以下几个:
det(a)求矩阵a的行列式
eig(a)求矩阵a的特点值
inv(a)或a^(-1)求矩阵a的逆矩阵
rank(a)求矩阵a的秩
trace(a)求矩阵a的迹(对角线元素之和)
我们进行工程运算经常常遇到矩阵对应元素之间的运算。
这种运算不同于前面讲的数学运算,为有所区别,我们称之为数组运算。
5.4差不多数学运算
数组的加、减与矩阵的加、减运算完全相同。
而乘除法运算有相当大的区别,数组的乘除法是指两同维数组对应元素之间的乘除法,它们的运算符为“.*”和“./”或“.\。
”前面讲过常数与矩阵的除法运算中常数只能做除数。
在数组运算中有了“对应关系”的规定,数组与常数之间的除法运算没有任何限制。
另外,矩阵的数组