电力系统分析课程设计
学院电气工程学院
班级
学号
姓名
指导教师
时间
电力系统潮流计算是研究电力系统稳态运行情况的一种计算,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态:各母线的电压,各元件中流过的功率,系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中,都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。可靠性和经济性。此外,电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和很基础的计算。
随着科学技术的发展,电力系统变得越来越复杂,电气工程师掌握一种好的能对电力系统进行仿真的软件是学习和研究的需要。与众多专门的电力系统仿真软件相比,MATLAB软件具有易学、功能强大和开放性好,是电力系统仿真研究的有力工具。
1 设计题目 (3)
1.1系统图的确定 (3)
1.2各节点的初值及阻抗参数 (4)
2 潮流计算 (5)
2.1潮流计算概述与发展 (5)
2.2复杂电力系统潮流计算 (5)
2.3 MATLAB概述 (6)
2.4牛顿-拉夫逊法原理 (7)
2.5牛顿-拉夫逊法解决潮流计算问题 (8)
2.6计算机潮流计算的步骤 (9)
2.7计算机潮流计算流程图 (11)
3 手算潮流计算 (12)
3.1确定节点类型 (12)
3.2求节点导纳矩阵Yb (12)
3.3计算各节点功率的修正方程的初始值(不平衡量) (13)
4 计算机算法潮流计算 (15)
4.1计算机法潮流计算过程 (15)
4.2计算机法潮流计算结果 (23)
4.3系统功率分布图 (25)
个人心得 (29)
附录:源程序 (30)
参考文献 (40)
1 设计题目
1.1系统图的确定
选择六节点、环网、两电源和多引出的电力系统,简化电力系统图如图1-1所示,等值导纳图如图1-2所示。运用以直角坐标表示的牛顿-拉夫逊计算如图1-1所示系统中的潮流分布。计算精度要求各节点电压的误差或修正量不大于5
ε。
10-
=
1.8+j0.4 1.6+j0.8
图1-1 电力系统图
图1-2 电力系统等值导纳图
1.2各节点的初值及阻抗参数
U=1.05+j0为定值,节点⑥为PV节点,该系统中,节点①为平衡节点,保持
1
其他四个节点都是PQ节点。给定的注入电压标幺值如表1-1、线路阻抗标幺值如表1-2、输出功率标幺值如表1-3。
表1-1 各节点电压标幺值参数
表1-2 线路、变压器阻抗标幺值
表1-3 节点输出功率
注:各PQ节点的电压取1是为了方便计算和最后验证程序的正确性。
2 潮流计算
2.1潮流计算概述与发展
电力系统潮流计算也分为离线计算和在线计算两种,前者主要用于系统规划设计和安排系统的运行方式,后者则用于正在运行系统的经常监视及实时控制。
利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。在这20年内,潮流计算曾采用了各种不同的方法,这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。
牛顿-拉夫逊法作为一种实用的,有竞争力的电力系统潮流计算方法,是在应用了稀疏矩阵技巧和高斯消去法求修正方程后。牛顿-拉夫逊法是求解非线性代数方程有效的迭代计算。
2.2复杂电力系统潮流计算
电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行方式下的节点电压和功率分布,用以检查系统各元件是否过负荷、各点电压是否满足要求、功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建,对新的电力系统进行规划设计都是以潮流计算为基础。
潮流计算结果的用途,例如用于电力系统稳定研究、安全估计或最优潮流等也对潮流计算的模型和方法有直接影响。
节点类型:
(1)PV节点:柱入有功功率P为给定值,电压也保持在给定数值。
(2)PQ节点:诸如有功功率和无功功率是给定的。
(3)平衡节点:用来平衡全电网的功率。选一容量足够大的发电机担任平衡全电网功率的职责。平衡节点的电压大小与相位是给定的,通常以它的相角为参考量,即取其电压相角为0。一个独立的电力网中只设一个平衡点。
基本步骤:
(1)形成节点导纳矩阵;
(2)将各节点电压设初值U;
(3)将节点初值代入相关求式,求出修正方程式的常数项向量;
(4)将节点电压初值代入求式,求出雅可比矩阵元素;
(5)求解修正方程,求修正向量;
(6)求取节点电压的新值;
(7)检查是否收敛,如不收敛,则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代,否则转入下一步;
(8)计算支路功率分布,PV节点无功功率和平衡节点柱入功率。
2.3 MATLAB概述
目前电子计算机已广泛应用于电力系统的分析计算,潮流计算是其基本应用之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求:(1)计算速度快;
(2)内存需要少;
(3)计算结果有良好的可靠性和可信性;
(4)适应性好,亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强;
(5)简单。
MATLAB是一种交互式、面向对象的程序设计语言,广泛应用于工业界与学术界,主要用于矩阵运算,同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。
MATLAB程序设计语言结构完整,且具有优良的移植性,它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题,特别是关于矩阵和矢量的计算。MATLAB与C语言和FORTRAN语言相比更容易被掌握。通过M语言,
可以用类似数学公式的方式来编写算法,大大降低了程序所需的难度并节省了时间,从而可把主要的精力集中在算法的构思而不是编程上。
2.4牛顿-拉夫逊法原理
假设有n 个联立的非线性代数方程:
??
??
???
===0)......,(0)......,(0)......,(21212211n n n n x x x f x x x f x x x f
假设以给出各变量的初值)0(1x ,)
0(2x ,……,)
0(n x ,令其分别为个变量的修正
量,使满足以上方程,所以:
??
?
?
???
=?+?+?+=?+?+?+=?+?+?+0),......,(0),......,(0),......,()0()0()0(2)0(2)0(1)0(1)
0()0()0(2)0(2)0(1)0(12)
0()0()0(2)0(2)0(1)0(11n n n n n n n x x x x x x f x x x x x x f x x x x x x f 将上式中的n 个多元函数在初始值附近分别展开成泰勒级数,并略去含有
)0(1x ?,)
0(2
x ?,……,)
0(n x ?的二次及以上阶次的各项,便得:
??????
??
??
?=?+?+?+
=?+?+?+=?+?+?+0......
),......,(0......),......,(0......),......,()
0(3
03
)
0(202
)
0(1
01)0()0(2)0(1)
0(3032)0(2022)0(1012)0()0(2)0(12)
0(303
1)0(2021)0(1011)
0()0(2)0(11x x f x x f x x f x x x f x x f x x f x x f x x x f x x f x x f x x f x x x f n n n n n n n αααααααααααααααααα
方程可写成:
??
???????????????????
?
?????????=????
??????)0()0(2)0(100
201
02
2201201
21011)0()0(2)0(1)0()0(2)0(12)
0()0(2)0(11*..................),......,(),......,(),......,(n n
n
n n n n n n n n x x x x f x f x f x f x f x
f x f x f x f x x x f x x x f x x x f αααααααααααααααααα
以上方程是对于修正量)0(1x ?,)
0(2x ?,……,)
0(n x ?的线性方程组,称为牛顿法的修正方程,可解出)0(1x ?,)
0(2x ?,……,)
0(n x ?。对初始近似解进行修正:
)
0()0()1(i i i x x x ?+= (i=1,2,……,n )
反复迭代,在进行k+1次迭代时,从求解修正方程式:
??
???????????????????
?
?????????-=????
??????)()(2)(121
2
22121
2111
)()(2)(1)()(2)(12)
()(2)(11*..................),......,(),......,(),......,(k n k k k n
n
k
n k
n k n k
k
k n k
k
k n k k n k n k k k n k k x x x x f x f x f x f x f x
f x f x f x f x x x f x x x f x x x f αααααααααααααααααα
得到修正量)(1k x ?,)
(2k x ?,……,)
(k n x ?,对各量进行修正)
()()
1(k i k i k i
x x x ?+=+
(i=1,2,……,n )迭代过程一直进行到满足收敛判据
1
)
()(2)(1|}),,(max {|λ?k n k k i x x x f ,
2.5牛顿-拉夫逊法解决潮流计算问题
节点总数为n ;PQ 节点有m ,;PV 节点有n-m-1,平衡节点有1个,节点编号按照先PQ 节点,再PV 节点,最后平衡节点的顺序进行编号,即:1,2,…,m 为PQ 节点;m+1,m+2,…,n-1为PV 节点;n 为平衡节点。可形成结点导纳矩阵。导纳矩阵元素可表示为ij
ij ij jB g Y +=,本文中节点电压以直角坐标形式表
示,即
i
i i jf e V +=.
。由此下列公式可求出Pi ,Qi
??
???
??
+--=++-=∑∑∑∑====n
j j ij j ij i n j j ij j ij i i n
j j ij j ij i n
j j ij j ij i i e B f G e f B e G f Q e B f G f f B e G e P 111
1)()()()(
假设系统中的第1,2,…,m 号节点为PQ 节点,第i 个节点的给定功率为is
P 和
is
Q ,对该节点可列方程:
),,2,1()()()()(1111
m i e B f G e f B e G f Q Q Q Q e B f G f f B e G e P P P P n j n
j j ij j i i j i j ij i is i is i n
j n
j j ij j i i j i j ij i is i is i ??=??
???
??++--=-=?+---=-=?∑∑∑∑====
假设系统中的第m+1,m+2,…,n-1号节点为PV 节点,则对其中每一个节点可列方程:
)1,,2,1()()()(22222211
-??++=??
?
??
+-=-=?+---=-=?∑∑==n m m i f e V V V V e B f G f f B e G e P P P P i i is i is i n
j n
j j ij j i i j i j ij i is i is i
第n 号节点为平衡节点,其电压为是给定的
n
n n jf e V +=.
,故不参加迭代。
修正方程V J W ?-=?可写成分块矩阵的形式:
?????????
????
??????????????????????
?-=??
?
?????
?????
???????---------13211,12
,11,11
,33231
1,222211,112
111321*n n n n n n n n n V V V V J J J J J J J J J J J J W W W W 通过反复求解修正方程,解出各节点的未知量,再通过收敛判据判定是否已为真值。从而求得PQ 节点的电压V 及相角δ的真值,PV 节点的Q 、δ真值,平衡节点的P 、Q 真值,以上即为牛顿-拉夫逊迭代法的潮流计算过程,其优点为计算精确,运行速度快。其中的各个环节都可通过MATLAB 程序来实现。
2.6计算机潮流计算的步骤
(1)对电力网络的所有参数设初值,包括电压、相角、有功、无功等。 (2)处理非标准变比支路,使其变成标准变比为1的变压器支路。 (3)形成节点导纳矩阵Y 。
(4)计算有功功率的不平衡量ΔPi,从而求出s)
i n,1,2,3...,(i /≠=?i i U P 。
(5)根据节点的类型形成J 。
(6)解修正方程式,求各节点的电压的变化量Δei(i=1,2,3...n,i ≠s)