计算机潮流计算.docx
- 文档编号:525565
- 上传时间:2023-04-29
- 格式:DOCX
- 页数:31
- 大小:271.59KB
计算机潮流计算.docx
《计算机潮流计算.docx》由会员分享,可在线阅读,更多相关《计算机潮流计算.docx(31页珍藏版)》请在冰点文库上搜索。
学号1350803233
《工厂供电》
课程设计
(2013级本科)
题目:
某工厂变电所一次系统设计
系(部)院:
物理与机电工程学院
专业:
电气工程及其自动化
作者姓名:
杨兴海
指导教师:
田娜职称:
助教
完成日期:
2016年1月5日
河西学院本科生课程设计任务书
设计题目
计算机潮流计算
作者姓名
杨兴海
学院、专业、年级
物电学院电气工程专业13级
指导教师姓名、职称
田娜助教
任务下达日期
2015年6月1日
一、设计内容:
结合《电力系统分析》课程,进一步学习并理解电力系统潮流计算的物理及数学模型,利用PQ分解方法对某电力系统进行潮流计算,根据潮流计算结果进行电压调整分析。
1.网络接线图及计算数据如下:
2.线路数据:
L1、L2和L3型号均为LGJQ4×400,长度分别为340km、280km和300km;L4和L5型号分别为LGJ240和LGJ185,长度分别为80km和8km。
变压器数据:
自行选择变压器型号。
母线2、10、6、11、8的电压等级分别为:
500kV、35kV、220kV、110kV和10kV。
变压器标幺值变比取1.0。
节点类型:
节点1、4分别为平衡节点和PV节点,其余为PQ节点。
发电机组数据:
等值发电机G2有功出力700MW,等值发电机G1有功最大出力700MW。
负荷数据:
节点2:
450+j100MVA;节点3:
300+j120MVA;
节点5:
350+j130MVA;节点7:
120+j80MVA;
节点9:
2.5+j0.9MVA;节点10:
30+j10MVA。
。
二、任务要求:
(1)查设计手册,确定变压器、线路等的电气参数;
(2)选择线路、变压器等模型,进行有名值和标幺值参数计算,并归算形成等值网络;
(3)手算Jaccobi矩阵的部分元素;
(4)应用程序进行潮流计算(机端母线电压取为1.0p.u.);
(5)分析不同节点电压水平的差别及其原因;
(6)选择适当调压方式进行电压调整;
(7)撰写设计报告。
三、论文(设计)进度安排
阶段
论文(设计)各阶段名称
起止日期
1
熟悉设计任务书、设计题目及设计背景资料
6.1~6.5
2
查阅有关资料
6.6~6.7
3
阅读设计要求必读的参考资料
6.8~6.10
4
书写设计说明书
6.11~6.19
5
小组答辩质疑
6.21~6.22
6
上交设计成果
6.30
四、需收集和阅读的资料及参考文献(指导教师指定)
[1]:
陈珩.电力系统稳态分析(第三版)[M],北京,中国电力出版社,2007
[2]:
何仰赞.温增银.《电力系统分析》第三版[M],武汉,华中科技大学出版社,2002
[3]:
陈悦.《电气工程毕业设计指南电力系统分册》[M],北京,中国水利水电出版社,2008
教研室意见
负责人签名:
年月日
目录
摘要:
1
第一章电力系统潮流计算概述 2
1.1潮流计算简介 2
1.2潮流计算的意义及其发展 2
第二章电气参数的确定 3
2.1变压器、线路参数确定 3
2.2电力线路、变压器模型的建立 4
第三章PQ分解方法潮流计算分析 5
3.1PQ分解法的极坐标表示及简化算法 5
3.1.1潮流计算的定义 5
3.1.2潮流计算的约束条件 6
3.1.3节点电压用极坐标表示时的牛顿-拉夫逊潮流计算 6
3.1.4对牛顿—拉夫逊法潮流计算的数学模型进行简化修正 8
3.2PQ分解法潮流计算的简化算法 10
3.3PQ分解法潮流计算的基本步骤 13
第四章程序编写及结果分析 13
4.1程序编程 13
4.2执行结果 20
第五章课程设计心得 22
参考文献 22
摘要
潮流计算是电力系统最基本最常用的计算。
根据系统给定的运行条件,网络接线及元件参数,通过潮流计算可以确定各母线的电压(幅值和相角),各支路流过的功率,整个系统的功率损耗。
潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。
因此,潮流计算在电力系统的规划计算,生产运行,调度管理及科学计算中都有着广泛的应用。
潮流计算在数学上是多元非线性方程组的求解问题,PQ分解法是数学上解非线性方程组的有效方法,有较好的收敛性。
运用电子计算机计算一般要完成以下几个步骤:
建立数学模型,确定解算方法,制订计算流程,编制计算程序。
关键字:
PQ分解法计算机潮流计算MATLAB
第一章电力系统潮流计算
1.1潮流计算简介
电力系统潮流计算是研究电力系统稳态运行情况的一种计算,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态:
各母线的电压,各元件中流过的功率,系统的功率损耗等等。
在电力系统规划的设计和现有电力系统运行方式的研究中,都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。
可靠性和经济性。
此外,电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。
所以潮流计算是研究电力系统的一种很重要和基础的计算。
电力系统潮流计算也分为离线计算和在线计算两种,前者主要用于系统规划设计和安排系统的运行方式,后者则用于正在运行系统的经常监视及实时控制。
利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。
在这20年内,潮流计算曾采用了各种不同的方法,这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。
对潮流计算的要求可以归纳为下面几点:
(1)计算方法的可靠性或收敛性;
(2)对计算机内存量的要求;
(3)计算速度;
(4)计算的方便性和灵活性。
电力系统潮流计算问题在数学上是一组多元非线性方程式求解问题,其解法都离不开迭代。
因此,对潮流计算方法,首先要求它能可靠地收敛,并给出正确答案。
由于电力系统结构及参数的一些特点,并且随着电力系统不断扩大,潮流计算的方程式阶数也越来越高,对这样的方程式并不是任何数学方法都能保证给出正确答案的。
这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。
1.2潮流计算的意义及其发展
(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。
(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。
(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。
(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。
总结为在电力系统运行方式和规划方案的研究中,都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。
同时,为了实时监控电力系统的运行状态,也需要进行大量而快速的潮流计算。
因此,潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。
在系统规划设计和安排系统的运行方式时,采用离线潮流计算;在电力系统运行状态的实时监控中则采用在线潮流计算。
近20多年来,潮流算法的研究仍然非常活跃,但是大多数研究都是围绕改进牛顿法和P-Q分解法进行的。
此外,随着人工智能理论的发展遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。
但是,到目前为止这些新的模型和算法还不能取代牛顿法和P-Q分解法的地位。
由于电力系统规模的不断扩大,对计算速度的要求不断提高,计算机的并行计算技术也将在潮流计算中得到广泛的应用,成为重要的研究领域。
第二章电气参数的确定
2.1变压器、线路参数确定
查设计手册:
LGJQ4×400型导线:
=0.020Ω/km,=0.276Ω/km,=4.02S/km;LGJ240型导线:
=0.131Ω/km,=0.432Ω/km;LGJ185型导线:
=0.170Ω/km,=0.440Ω/km。
变压器型号:
T1,T2为SF-360000/220;T3部分参数为:
额定容量240/120/120MVA,额定电压500/242/38.5kv。
T4部分参数为:
额定容量120/120/60MVA,,额定电压242/121/10.5kv,等值电路中所有参数都归算到高压侧。
线路L1,L2的电导,L4,L5的导纳都可略去。
表2-1短路电压(未经归算)
绕组
高压-中压
高压-低压
中压-低压
短路电压百分值(%)
10
30
20
2.1.1电力变压器及线路参数的计算
1.采用有名值:
2.采用标幺值:
取基准功率为1000MVA,基准电压为500kv
2.2.2电力变压器及线路参数的计算结果
电力线路参数见表2-1,电力变压器参数见表2-2
表2-1电力线路参数
线路
参数
LGJQ4×400
LGJ240
L4
LGJ185
L5
L1
L2
L3
6.8
5.6
6.0
10.48
1.36
0.0272
0.0224
0.0240
0.0419
0.0054
93.84
77.28
82.8
34.56
3.52
0.3754
0.3091
0.3312
0.1382
0.0141
6.8+j93.84
5.6+j77.28
6.0+j82.8
10.48+j34.56
1.36+j3.52
0.0272+j0.3754
0.0224+j0.3091
0.024+j0.3312
0.0419+j0.1382
0.0054+j0.0141
S
1.3668
1.1256
1.206
/
/
3.417
2.814
3.015
/
/
表2-2电力变压器参数
变压器
参数
T1
T2
T3
T4
25.54
25.54
0.243
0.142
0.486
0.142
0.486
0.284
0.1022
0.1022
1.92
1.92
104.17
48.8
0
0
208.33
97.6
0.00768
0.00768
0.4167
0.1952
0
0
0.8333
0.3904
25.54+j1.92
25.54+j1.92
0.243+j104.17
0.142+j48.8
0.486
0.142
0.486+j208.33
0.284+j97.6
0.1022+j0.00768
0.1022+j0.00768
+j0.4167
+j0.1952
+j0.8333
+j0.3904
×S
17.77
17.77
9.56
5.46
0.00444
0.00444
0.00239
0.001365
×S
20.83
20.83
16.39
12.29
0.005208
0.005208
0.004098
0.003073
Ω
17.77+j20.83
17.77+j20.83
9.56+j16.39
5.46+12.29
4.44+j5.208
4.44+j5.208
2.39+j4.098
1.365+j3.073
2.2电力线路、变压器模型的建立
线路采用Π型等值电路,变压器采用Γ型等值电路,有名值和标幺值等值电路分别见:
图2-1、图2-2。
图2-1有名值等值网络
图2-2标幺值等值网络
第三章PQ分解方法潮流计算分析
3.1PQ分解法的极坐标表示及简化算法
3.1.1潮流计算的定义
潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件,确定电力系统各部分稳态运行状态参数的计算。
通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。
待求的运行状态参量包括电网各母线节点的电压幅值和相角,以及各支路的功率分布、网络的功率损耗等。
PQ分解法的极坐标表示是派生于以极坐标表示的牛顿—拉夫逊法。
3.1.2潮流计算的约束条件
为了保证电力系统的正常运行潮流问题中某些变量应满足一定的约束条件,常用的约束条件有:
(1)所有节点电压必须满足
(i=1,2,…,n)
从保证电能质量和供电安全的要求看,电力系统的所有电气设备都必须运行在额定电压附近。
PV节点的电压幅值必须按上述条件给定。
因此,这一约束主要是对PQ节点而言。
(2)所有电源节点的有功功率和无功功率必须满足
PQ节点的有功功率和无功功率以及PV节点的有功功率,在给定时就必须满足上述条件。
因此,对平衡节点的P和Q以及PV节点的Q应按上述条件进行检验。
(1)某些节点之间电压的相位差应满足
为了保证系统运行的稳定性,要求某些输电线路两端的电压相位差不超过一定的数值。
因此,潮流计算可以归纳为求解一组非线性方程组,并使其解答满足一定的约束条件。
如果不满足,则应修改某些变量的给定值,甚至修改系统的运行方式,重新进行计算。
3.1.3节点电压用极坐标表示时的牛顿-拉夫逊潮流计算
采用极坐标时,节点电压表示为
节点功率方程(11-25)将写成
公式1
式中,,是i,j两节点电压的相差角。
方程式(公式1)把节点功率表示为节点电压的幅值和相角的函数。
在有n节点的系统中,假定第1~m号节点PQ节点,第m+1~n-1号节点的PV节点,第n号节点为平衡节点。
和是给定的,PV节点的电压幅值Vm+1~Vn-1也是给定的。
因此,只剩下n-1个节点的电压相角和m个节点的电压幅值,,…,是未知量。
实际上,对于每一个PQ节点或每一个PV节点都可以列写一个有功功率不平衡量方程式
=-=-(i=1,2,…,n-1)公式2
而对于每一个PQ节点还可以再列写一个无功功率不平衡量方程式
=-=-(i=1,2,…,m)
公式3
式2和式3一共包含了n-1+m个方程式,正好同未知量的数目相等,而比直角坐标形式的方程式少了n-1-m个。
对于方程式2和式3可以写出修正方程式如下:
公式4
式中
;;;公式5
H是(n-1)×(n-1)阶方阵,其元素为;N是(n-1)×m阶矩阵,其元素为;K是m×(n-1)阶矩阵,其元素为;L是m×m阶方阵,其元素为。
3.1.4对牛顿—拉夫逊法潮流计算的数学模型进行简化修正
在交流高压电网中,输电线路的电抗要比电阻大得多,系统中母线有功功率的变化主要受电压相位的影响,无功功率的变化则主要受母线电压幅值变化的影响。
在修正方程式的系数矩阵中,偏导数和的数值相对于偏导数和是相当小的。
作为简化的第一步,可以将上式子块N和K略去不计,即认为它们的元素都等于零。
这样,n-1+m阶的方程式便分解为一个n-1阶和一个m阶的方程。
这一简化大大地节省了机器内存和解题时间。
以上方程式表明,节点的有功功率不平衡量只用于修正电压的相位,节点的无功功率不平衡量只用于修正电压的幅值。
这两组方程分别轮流进行迭代,这就是所谓有功-无功功率分解法。
但矩阵H和L的元素都是节点电压幅值和相角差的函数,其数值在迭代过程中是不断变化的。
因此,最关键的一步简化就在于,把系数矩阵H和L简化成常数矩阵。
它的根据是什么呢?
在一般情况下,线路两端电压的相角差是不大的(不超过100~200),因此可以认为
此外,与系统各节点无功功率相适应的导纳必远小于该节点自导纳的
或
考虑到以上的关系,矩阵H和L的元素的表达式便被简化成
(i,j=1,2,…,n-1)公式6
(i,j=1,2,…,m)公式7
而系数矩阵H和L则可以分别写成
H=
=×
=公式8
L=
=×
=公式9
代入得:
公式10
公式11
这就是简化了的修正方程,也可展开为
=-公式12
=-公式13
在这两个修正方程式中,系数矩阵都由节点导纳矩阵的虚部构成,只是阶次不同,矩阵为n-1阶,不含平衡节点对应的行和列,矩阵为m阶,不含平衡节点和PV节点所对应的行和列。
由于修正方程的系数矩阵为常数矩阵,只要作一次三角分解,即可反复使用,结合采用稀疏技巧,还可进一步的节省机器内存和计算时间。
3.2PQ分解法潮流计算的简化算法
PQ分解法潮流计算修正方程的简化根据极坐标表示的牛顿-拉夫逊法可得:
每一个PQ或PV节点的有功功率不平衡方程:
=-=-(i=1,2,…,n-1)
公式14
每一个PQ节点的无功功率不平衡方程:
=-=-(i=1,2,…,m)
公式15
由此可以写出修正方程如下:
公式16
对式24和式25求偏导数,可以得到雅克比矩阵元素的表达如下:
a)当i≠j时
b)当i=j时
对修正方程的第一个简化是:
将式26中的N、K子阵略去而将其简化为
公式17
上式可分别写成以下两式
公式18
公式19
对修正方程式的第二个简化式:
在一般情况下,线路两端电压的相角差是不大的(不超过100~200),因此可以认为
因此可得:
(i,j=1,2,…,n-1)
(i,j=1,2,…,m)
对修正方程的第三个简化是:
当采用标幺值计算是,可近似地认为各节点电压的大小,且当忽略接地支路时,导纳元素的对角元为非对角元之和,这样可得:
且此时式29也可以简化为
公式20
这是雅克比矩阵的两个子阵H、L具有相同的表达式,只是其阶数不同,前者为(n-1)阶,后者为m阶,两个子阵的展开式如下:
将式28代入式30可得:
=-公式21
=-公式22
式31和式32可简写为
公式23
式33为简化P—Q分解法的修正方程式,原P—Q分解法的修正方程的简化形式
通过两者的比较可知,简化P—Q分解法的修正方程式比原P—Q分解法的修正方程更为简单。
3.3PQ分解法潮流计算的基本步骤
(1)形成系数矩阵,,并求其逆矩阵。
(2)设各节点电压的初值为(i=1,2,…,n,i≠s)和(i=1,2,…,m,i≠s)。
(3)通过有功功率的不平衡方程计算有功功率的不平衡量,从而求出(i=1,2,…,n,i≠s)。
(4)解修正方程式,求各节点电压相位角的变量(i=1,2,…,n,i≠s)。
(5)求各节点电压相位角的新值(i=1,2,…,n,i≠s)。
(6)通过无功功率的不平衡方程计算无功功率的不平衡量,从而求出(i=1,2,…,m,i≠s)。
(7)解修正方程式,求各节点电压大小的变量(i=1,2,…,m,i≠s)。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算机 潮流 计算