




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、课程设计(论文)题 目 名 称 电力系统潮流计算 课 程 名 称 电力系统稳态分析 学 生 姓 名 学 号 系 、专 业 电气工程系电气工程及其自动化电力方向 指 导 教 师 2009年1月6日前 言 在如今的社会,电力已经成为人们必不可少的需求,而建立结构合理的大型电力系统不仅便于电能生产与消费的集中管理、统一调度和分配,减少总装机容量,节省动力设施投资,且有利于地区能源资源的合理开发利用,更大限度地满足地区国民经济日益增长的用电需要。电力系统建设往往是国家及地区国民经济发展规划的重要组成部分。电力系统的出现,使高效、无污染、使用方便、易于调控的电能得到广泛应用,推动了社会生产各个领域的变化
2、,开创了电力时代,发生了第二次技术革命。电力系统的规模和技术水准已成为一个国家经济发展水平的标志之一。电力系统稳态分析包括潮流计算(或潮流分析)和静态安全分析。潮流计算针对电力革统各正常运行方式,而静态安全分析则要研究各种运行方式下个别系统元件退出运行后系统的状况。其目的是校验系统是否能安全运行,即是否有过负荷的元件或电压过低的母线等。原则上讲,静态安全分析也可U用潮流计算来代替。但是一般静态安全分析需要校验的状态数非常多,用严格的潮流计算来分析这些状态往往计算量过大,因此不得不寻求一些特殊的算法以满足要求。牛顿法是数学中解决非线性方程式的典型方法,有较好的收敛性。解决电力系统潮流计算问题是以
3、导纳距阵为基础的,因此,只要在迭代过程中尽可能保持方程式系数距阵的稀疏性,就可以大大提高牛顿法潮流程序的放率。自从20 世纪60 年代中期利用了最佳顺序消去法以后,牛顿法在收敛性、内存要求、速度方面都超过了阻抗法,成为直到目前仍在广泛采用的优秀方法。目 录第一章 系统概述1.1 设计目的与要求1.1.1 设计目的41.1.2 设计要求41.2 设计题目1.3 设计内容第二章 潮流计算设计题目2.1 潮流计算题目2.2 对课题的分析及求解思路第三章 潮流计算算法及手工计算3.1 潮流计算算法3.2 关于电力系统潮流计算手工计算3.2.1 节点导纳矩阵83.2.2 简化雅可比矩阵9修正、迭代10第
4、四章 Matlab概述4.1 Matlab简介4.2 矩阵的运算4.2.1 四则运算124.2.2 与常数的运算124.2.3 基本数学运算124.2.4 逻辑关系运算12第五章 潮流计算流程图及源程序5.1 潮流计算流程图5.2 潮流计算源程序5.3 运行计算结果总 结参考文献1.1 设计目的与要求 设计目的1.1.2 设计要求1.2 设计题目电力系统潮流计算(牛顿-拉夫逊法、P-Q 分解法)1.3 设计内容8.准备计算机演示答辩,书写该课程设计说明书(必须计算机打印)。2.1 潮流计算题目图2-1 电力系统接线图2.2 对课题的分析及求解思路此电力系统是一个5节点,4支路的电力网络。其中包
5、含3个PQ节点,一个PV节点,和一个平衡节点。综合比较牛顿拉夫逊法(直角坐标、极坐标)、PQ分解法等多种求解方法的特点,最后确定采用牛顿拉夫逊法(极坐标)。因为此方法所需解的方程组最少。3.1 潮流计算算法本题采用了题目要求的牛顿拉夫逊潮流计算的方法。牛顿-拉夫逊法潮流计算的公式。把牛顿法用于潮流计算,采用直角坐标形式表示的如式(1-3)所示的形式。其中电压和支路导纳可表示为: (1-2) 将上述表示式(1-2)代入(1-1)式的右端,展开并分出实部和虚部,便得: (1-3) 按照以上的分类,PQ节点的输出有功功率和无功功率是给定的,则第i节点的给定功率设为和(称为注入功率)。 假定系统中的第
6、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) 由上式不难看出,雅可比矩阵有
7、以下特点。 雅可比矩阵中的诸元素都是节点电压的函数,因此在迭代过程中,它们将随着节点电压的变化而不断的变化。 雅可比矩阵具有结构对称性,数据不对称。如非对角,。 由式(1-7)可以看出,当导纳矩阵中非对角元素为零时,。雅可比矩阵中相应的元素也为零,即矩阵是非常稀疏的。因此,修正方程的求解同样可以应用稀疏矩阵的求解技巧。正是由于这一点才使N-R法获得广泛的应用。3.2 关于电力系统潮流计算手工计算 节点导纳矩阵求得节点导纳矩阵Y = 各节点的导纳值如下:; ; ; ;; ;; ; ; ; ; ; ;. 简化雅可比矩阵形成有功迭代和无功迭代的简化雅可比矩阵B/和B/B/= B/= 将B/ 和B/进
8、行三角分解: -0.234654 -0.564073 -0.443749 -0.423068 2.357080 -1.02902 -0.234654 -0.564073 -0.423068 修正、迭代给定PQ节点初值和各节点电压相角初值V1=1.050。 ,V2(0)=V3(0)=1.0,V4=1.12(0)=3(0)=0, 4(0)=01 作第一次有功迭代,按公式计算节点有功功率不平衡量 P2(0)=-0.55-(-0.024037)=-0.525963 P3(0)=-0.30-(-0.022695)=-0.277305 P4(0)=0.500000 P1(0)/V1(0)=0.454545
9、 P2(0)/ V2(0)=-0.525963 P3(0)/V3(0)=-0.2773092做第一次无功迭代,按公式计算无功功率不平衡量,计算时电压相角最新的修正值。 Q2(0)=-0.13-(-0.001550)=-0.039594Q3(0)=-0.18-(-0.14406)=-0.039588Q2(0)/ V2(0)=-0.131553Q3(0)/V3(0)=-0.039588 解修正方程式,可得各节点电压幅值的修正量为 V3(0))=-0.014855 于是有: V2(1) = V2(0)+V2(1)=0.964776 V3(1) = V3(0)+V3(1)=0.985145 到这里为止
10、,第一轮有功迭代和无功迭代便做完了。3 按公式计算平衡节点功率,得: P1+jQ1=0.367885+j0.264696经过四轮迭代,节点不平衡功率也下降到10-5以下,迭代到此结束。4.1 Matlab简介目前电子计算机已广泛应用于电力系统的分析计算,潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求:(1)计算速度快(2)内存需要少(3)计算结果有良好的可靠性和可信性(4)适应性好,亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强(5)简单。 MATLAB是一种交互式、面向对象的程序设计语言,广泛应用于工业界与学术界,主要用于矩阵运算,同时
11、在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。MATLAB程序设计语言结构完整,且具有优良的移植性,它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题,特别是关于矩阵和矢量的计算。MATLAB与C语言和FORTRAN语言相比更容易被掌握。通过M语言,可以用类似数学公式的方式来编写算法,大大降低了程序所需的难度并节省了时间,从而可把主要的精力集中在算法的构思而不是编程上。另外,MATLAB提供了一种特殊的工具:工具箱(TOOLBOXES).这些工具箱主要包括:信号处理(SIGNAL PROCESSING)、控制系统(CONTROL SYSTEMS)
12、、神经网络(NEURAL NETWORKS)、模糊逻辑(FUZZY LOGIC)、小波(WAVELETS)和模拟(SIMULATION)等等。不同领域、不同层次的用户通过相应工具的学习和应用,可以方便地进行计算、分析及设计工作。MATLAB设计中,原始数据的填写格式是很关键的一个环节,它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计,主要应从使用的角度出发,原则是简单明了,便于修改。 与常数的运算 常数与矩阵的运算即是同该矩阵的每一元素进行运算。但需注意进行数除时,常数通常只能做除数。基本函数运算中,矩阵的函数运算是矩阵运算中最实用的部分,常用的主要有以下几个:det(a)
13、求矩阵a的行列式eig(a) 求矩阵a的特征值inv(a)或a (-1) 求矩阵a的逆矩阵rank(a) 求矩阵a的秩trace(a) 求矩阵a的迹(对角线元素之和)我们在进行工程计算时常常遇到矩阵对应元素之间的运算。这种运算不同于前面讲的数学运算,为有所区别,我们称之为数组运算。 基本数学运算数组的加、减与矩阵的加、减运算完全相同。而乘除法运算有相当大的区别,数组的乘除法是指两同维数组对应元素之间的乘除法,它们的运算符为“.*”和“./”或“.”。前面讲过常数与矩阵的除法运算中常数只能做除数。在数组运算中有了“对应关系”的规定,数组与常数之间的除法运算没有任何限制。另外,矩阵的数组运算中还有
14、幂运算(运算符为 . )、指数运算(exp)、对数运算(log)、和开方运算(sqrt)等。有了“对应元素”的规定,数组的运算实质上就是针对数组内部的每个元素进行的。矩阵的幂运算与数组的幂运算有很大的区别。 逻辑关系运算 逻辑运算是MATLAB中数组运算所特有的一种运算形式,也是几乎所有的高级语言普遍适用的一种运算。5.1 潮流计算流程图图5-1 潮流计算流程图5.2 潮流计算源程序据课题题目,本程序把节点1设为平衡节点,节点2、3、4为PQ节点,节点5为PV节点。G(1,1)=10.834;B(1,1)=-32.500;G(1,2)=-1.667;B(1,2)=5.000;G(1,3)=-1
15、.667;B(1,3)=5.000;G(1,4)=-2.500;B(1,4)=7.500;G(1,5)=-5.000;B(1,5)=15.000;G(2,1)=-1.667;B(2,1)=5.000;G(2,2)=12.917;B(2,2)=-38.750;G(2,3)=-10.000;B(2,3)=30.000;G(2,4)=0;B(2,4)=0;G(2,5)=-1.250;B(2,5)=3.750;G(3,1)=-1.667;B(3,1)=5.000;G(3,2)=-10.000;B(3,2)=30.000;G(3,3)=12.917;B(3,3)=-38.750;G(3,4)=-1.25
16、0;B(3,4)=3.750;G(3,5)=0;B(3,5)=0;G(4,1)=-2.500;B(4,1)=7.500;G(4,2)=0;B(4,2)=0;G(4,3)=-1.250;B(4,3)=3.750;G(4,4)=3.750;B(4,4)=-11.250;G(4,5)=0;B(4,5)=0;G(5,1)=-5.000;B(5,1)=15.000;G(5,2)=-1.250;B(5,2)=3.750;G(5,3)=0;B(5,3)=0;G(5,4)=0;B(5,4)=0;G(5,5)=6.250;B(5,5)=-18.750;Y=G+j*B;delt(1)=0;delt(2)=0;de
17、lt(3)=0;delt(4)=0;u(1)=1.0;u(2)=1.0;u(3)=1.0;u(4)=1.0;p(1)=0.20;q(1)=0.20;p(2)=-0.45;q(2)=-0.15;p(3)=-0.40;q(3)=-0.05;p(4)=-0.60;q(4)=-0.10;k=0;precision=1;N1=4; %the N1 is the amount of the PQ buswhile precision>0.00001 delt(5)=0; u(5)=1.06; for m=1:N1 for n=1:N1+1 pt(n)=u(m)*u(n)*(G(m,n)*cos(de
18、lt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n); qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); end pp(m)=p(m)-sum(pt); qq(m)=q(m)-sum(qt); end for m=1:N1 for n=1:N1+1 h0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(
19、n)+B(m,n)*sin(delt(m)-delt(n); j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n); L0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); end H(m,m)=sum(h0)-u(m)2*(G(m,m)*sin(delt(m)-delt(m)-B(m,m)*cos(delt(m)-delt(m); N(m,m)=sum(n0)-2*u(m)2*G(m,m)+u(m)2*(G(m,m
20、)*cos(delt(m)-delt(m)+B(m,m)*sin(delt(m)-delt(m); J(m,m)=sum(j0)+u(m)2*(G(m,m)*cos(delt(m)-delt(m)+B(m,m)*sin(delt(m)-delt(m); L(m,m)=sum(L0)+2*u(m)2*B(m,m)+u(m)2*(G(m,m)*sin(delt(m)-delt(m)-B(m,m)*cos(delt(m)-delt(m); end for m=1:N1 JJ(2*m-1,2*m-1)=H(m,m); JJ(2*m-1,2*m)=N(m,m); JJ(2*m,2*m-1)=J(m,m)
21、; JJ(2*m,2*m)=L(m,m); end for m=1:N1 for n=1:N1 if m=n else H(m,n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); J(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n); N(m,n)=-J(m,n); L(m,n)=H(m,n); JJ(2*m-1,2*n-1)=H(m,n); JJ(2*m-1,2*n)=N(m,n); JJ(2*m,2*n-1)=J(
22、m,n); JJ(2*m,2*n)=L(m,n); end end end for m=1:N1 PP(2*m-1)=pp(m); PP(2*m)=qq(m); end uu=-inv(JJ)*PP' precision=max(abs(uu); for n=1:N1 delt(n)=delt(n)+uu(2*n-1); u(n)=u(n)+uu(2*n); end k=k+1;endK=k-1,delt,u'%the following program is used to calculate the S5 and S(m,n)for n=1:N1+1 U(n)=u(n)*(
23、cos(delt(n)+j*sin(delt(n);end for m=1:N1+1 I(m)=Y(5,m)*U(m);end S5=U(5)*sum(conj(I) for m=1:N1+1 for n=1:N1+1 S(m,n)=U(m)*(conj(U(m)-conj(U(n)*conj(-Y(m,n); endendYS5.3 运行计算结果K =4;delt = -0.0461 -0.0839 -0.0896 -0.1044 0;U = 1.0365 1.0087 1.0073 1.0016 1.0600S5 =1.2982 + 0.2445iY = 10.8340 -32.5000i -1.6670 + 5.0000i -1.6670 + 5.0000i -2.5000 + 7.5000i -5.0000 +15.0000i -1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i 0 -1.2500 + 3.7500i -1.6
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 市总工会领导班子述职报告
- 斯玛特维商场互动导视系统解决方案v10(16-9)
- 2025年会计、审计及税务服务合作协议书
- 广东省广州市华侨、协和、增城中学等三校2024~2025学年高一下学期期中考试数学试卷(解析版)
- 安徽省鼎尖联考2024-2025学年高二下学期4月月考数学试题(解析)
- 2025年驾校学车项目建议书
- 2025年视听周边设备:耳机合作协议书
- 晚期肝癌护理措施
- 护理措施诊断
- 风疹患者护理规范
- JHA及SCL风险评价方法讲解(参考)
- DB11T 1933-2021 人乳库建立与运行规范
- 1.3.1动量守恒定律课件(共13张PPT)
- 国网北京市电力公司授权委托书(用电)
- 白黑白装饰画欣赏黑白装饰画的特点黑白装饰画的表现形式黑白装饰 bb
- 中小学教育惩戒规则(试行)全文解读ppt课件
- 调度指挥与统计分析课程教学设计
- 常暗之厢(7规则-简体修正)
- 终端塔基础预偏值(抬高值)计算表格
- 海外医疗服务委托合同协议书范本模板
- (完整版)研究者手册模板
评论
0/150
提交评论