节点导纳矩阵及潮流计算Word格式.docx
《节点导纳矩阵及潮流计算Word格式.docx》由会员分享,可在线阅读,更多相关《节点导纳矩阵及潮流计算Word格式.docx(15页珍藏版)》请在冰豆网上搜索。
从而可以得到
n个节点时的节点导纳矩阵方程
组:
Y11V1Y12V2LY1nVnI1
Y21V1Y22V2LY2nVnI2
L
(2-1)
Yn1V1Yn2V2LYnnVnIn
由此可以得到n个节点导纳矩阵:
它反映了网络的参数及接线情况,因此导纳矩阵可以看成是对电力网络电气特性的一种数学抽象。
由导纳短阵所联系的节点方程式是电力网络广泛应用的一种数学模型。
通过上面的讨论,可以看出节点导纳矩阵的有以下特点:
(1)导纳矩阵的元素很容易根据网络接线图和支路参数直观地求得,形成节点导纳矩阵的程序比较简单。
(2)导纳矩阵为对称矩阵。
由网络的互易特性易知YijYji。
(3)导纳矩阵是稀疏矩阵。
它的对角线元素一般不为零,但在非对角线元素中则存在不少零元素。
在电力系统的接线图中,一般每个节点与平均不超过3~4个其他节点有直接的支路连接。
因此,在导纳矩阵的非对角线元素中每行仅有3~4个非零元素,其余的都是零元素,而且网络的规模越大,这种现象越显著。
节点导纳矩阵的形式可归纳如下:
(1)导纳矩阵的阶数等于电力网络
(2)导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连得不接地支路数。
(3)导纳矩阵各对角元素,即节点的自导纳等于相应节点之间的支路导纳之和。
(4)导纳矩阵非对角元素,即节点之间的互导纳等于相应节点之间的支路导纳的负值。
2.2牛顿-拉夫逊法
2.2.1牛顿-拉夫逊法基本原理
牛顿--拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法。
其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程。
即通常所称的逐次线性化过程。
对于非线性代数方程组:
在待求量x的某一个初始估计值x(0)附近,将上式展开成泰勒级数并略去二阶及以上的高阶项,得到如下的经线性化的方程组:
f(x(0))f'
(x(0))x(0)0(2-4)
上式称之为牛顿法的修正方程式。
由此可以求得第一次迭代的修正量
x(0)[f'
(x(0))]1f(x(0))(2-5)
将x(0)和x(0)相加,得到变量的第一次改进值x
(1)。
接着就从x
(1)出发,重复上述计算过程。
因此从一定的初值x(0)出发,应用牛顿法求解的迭代格式为:
'
(k)(k)(k)
f(x)xf(x)(2-6)
上两式中:
f'
(x)是函数f(x)对于变量x的一阶偏导数矩阵,即雅可比矩阵J;
k为迭代次数。
由上式可见,牛顿法的核心便是反复形式并求解修正方程式。
牛顿法当初始估计值x(0)和方程的精确解足够接近时,收敛速度非常快,具有平方收敛特性。
牛顿潮流算法突出的优点是收敛速度快,若选择到一个较好的初值,算法将具有平方收敛特性,一般迭代4~5次便可以收敛到一个非常精确的解。
而且其迭代次数与所计算网络的规模基本无关。
牛顿法也具有良好的收敛可靠性,对于对以节点导纳矩阵为基础的高斯法呈病态的系统,牛顿法也能可靠收敛。
牛顿法所需的内存量及每次迭代所需时间均较高斯法多。
牛顿法的可靠收敛取决于有一个良好的启动初值。
如果初值选择不当,算
法有可能根本不收敛或收敛到一个无法运行的节点上。
对于正常运行的系统,
各节点电压一般均在额定值附近,偏移不会太大,并且各节点间的相位角差也不大,所以对各节点可以采用统一的电压初值(也称为平直电压),如假定:
Ui(0)1i(0)0或ei(0)1fi(0)0(iq1,2,L,n;
is)(2-8)
这样一般能得到满意的结果。
但若系统因无功紧张或其它原因导致电压质量很差或有重载线路而节点间角差很大时,仍用上述初始电压就有可能出现问题。
解决这个问题的办法可以用高斯法迭代1~2次,以此迭代结果作为牛顿法的初值。
也可以先用直流法潮流求解一次以求得一个较好的角度初值,然后转入牛顿法迭代。
2.2.2牛顿--拉夫逊法潮流求解过程介绍
以下讨论的是用直角坐标形式的牛顿—拉夫逊法潮流的求解过程。
当采用直角坐标时,潮流问题的待求量为各节点电压的实部和虚部两个分量e1,f1,e2,f2...en,fn由于平衡节点的电压向量是给定的,因此待求两共2(n-1)需
要2(n-1)个方程式。
事实上,除了平衡节点的功率方程式在迭代过程中没有约束作用以外,其余每个节点都可以列出两个方程式。
对PQ节点来说,Pis和Q是给定的,因而可以写出
Pipisei
(GijejBijfj)fjjij
i(Gijf
i
Bijej)0
(2-
QiQisf
i(GijejBijfj)ejji
(Gijf
ji
jBijej)0
9)
对PV
节点来说,给
定量
是Pis和Vis,
因此可以列出
PiPisei
j
(GijejBijfj)fi
iji
(Gijfj
(2-10)
2222
Vi2Vi2s(ei2fi)0求解过程大致可以分为以下步骤:
(1)形成节点导纳矩阵;
(2)将各节点电压设初值U
3)将节点初值代入相关求式,求出修正方程式的常数项向量;
(4)将节点电压初值代入求式,求出雅可比矩阵元素;
(5)求解修正方程,求修正向量;
(6)求取节点电压的新值;
(7)检查是否收敛,如不收敛,则以各节点电压的新值作为初值自第3步
重新开始进行狭义次迭代,否则转入下一步;
8)计算支路功率分布,PV节点无功功率和平衡节点注入功率以直角坐标系形式表示:
○1迭代推算式
采用直角坐标时,节点电压相量及复数导纳可表示为:
V&
ieijfi
iii(2-11)
YijGijjBij
将以上二关系式代入上式中,展开并分开实部和虚部;
假定系统中的第
1,2,L,m号为P—Q节点,第m+1,m+2,L,n-1为P—V节点,根据节点性质的不同,得到如下迭代推算式:
⑴于PQ节点
⑵对于PV节点
⑶对于平衡节点
平衡节点只设一个,电压为已知,不参见迭代,其电压为:
Vnenjfn(2-14)
○2修正方程
两组迭代式中包括2(n-1)个方程.选定电压初值及变量修正量符号之
后代入,并将其按泰勒级数展开,略去ei,fi二次方程及以后各项,得到修
正方程如下:
W
JU
(2-15)
P1
e1
Q1
f1
M
其中,
Pm;
em
Qm
Ufm
Pm1
em1
U2
m1
fm1
Pn1
en1
U2n1
fn1
fm
em1
fm1
fn1
Pm
U2m1
Um1
fn1
U2n1
(2-16)
○3雅可比矩阵各元素的算式
式(2-12)中,雅可比矩阵中的各元素可通过对式(2-8)和(2-9)进行偏导而求得.当ji时,雅可比矩阵中非对角元素为
ej
当ji时,雅可比矩阵中对角元素为
由式(2-13)和(2-18)看出,雅可比矩阵的特点:
○1矩阵中各元素是节点电压的函数,在迭代过程中,这些元素随着节点电压的变化而变化;
○2导纳矩阵中的某些非对角元素为零时,雅可比矩阵中对应的元素也是为零若Yij0,则必有Jij0;
雅可比矩阵各元素的表示如下:
3分析计算
1.根据给定的运行条件,确定图中所示电力系统潮流计算时各节点的类型和待求量
根据图中可以看出各节点的类型和待求量分别为:
节点1:
PQ节点待求量:
U,
节点2:
PV节点待求量:
Q,
节点3:
平衡节点待求量:
P,Q
2.求节点导纳矩阵Y
Y22y12y231.3j7
Y33
y23y131.55j6.5
Y21
Y12
0.5
j3
Y32
Y23
0.8
j4
Y31
Y13
0.75
j2.5
所以节点导纳矩阵为:
1.25
j5.5
Y0.5
1.3
j7
1.55
j6.5
3.潮流方程或功率方程的表达式
因为对n个节点的网络,电力系统的潮流方程一般形式是:
PijQin.(i=1,2,⋯,n)
*YijVj
其中Pi=PGi-PLdi,Qi=QGi-QLdi,即PQ分别为节点的有功功率无功功率。
所以代入得潮流方程:
1)修正方程
计算1、2节点的不平衡量Pi、Qi和Vi
20202
V220V2SV200节点3是平衡节点,其电压Vieijfi是给定的,故不参加迭代。
根据给定的容许误差105,按收敛判据maxPik,Qik,Vi2k进
行校验,以上节点1、2的不平衡量都未满足收敛条件,于是继续以下计算。
修正方程式为WJV(n=3)
W[P1
2T
P2V22]T
Ve1
e2f2T
e2f2
Je1
JP2
P2
V22
以上雅可比矩阵J中的各元素值是通过求偏导数获得的,对PQ节点来说,
Pis和Qis是给定的,因而可以写出
Pipisei(GijejBijfj)fj(GijfjBijej)0
jiji
QiQisfi(GijejBijfj)ej(GijfjBijej)0
对PV节点来说,给定量是Pis和Vis,因此可以列出
PiPisei(GijejBijfj)fi(GijfjBijej)0
Vi2Vei2s(i2fi)0
当ji时,雅可比矩阵中非对角元素为
当ji时,雅可比矩阵中对角元素为:
iei
(Gijej
j1
Bijfj)
Giiei
Biifi
Pifj
n
(Gijfj
Bijej)
Giifi
Biiei
Qiei
(Gijfj
Qifj
(Gijej
Biif
Ui2ej
2ei
Ui2
2fi
fi
代入数值后的修正方程为
5.5
3
2
1
7
e2
f2
求解修正方程得
e10.2547
0.3611
e20
0.1015
2)收敛条件
一轮迭代结束,根据收敛条件收敛判据maxPik,Qik,Vi2k,若等
式成立,结果收敛,迭代结束,计算平衡节点的功率和线路潮流计算,否则继续计算雅可比矩阵,解修正方程,直到满足收敛判据。
4结果分析
给定节点电压初值e10e20e301.0,f10f20f300,经过四次笔算迭代过程后,得到节点电压和不平衡功率的变化情况分别于表4.1和表4.2所示(取105):
迭代计数k
节点电压
?
V1e1jf1
V2e2jf2
V3e3jf3
0.7453-j0.3611
1-j0.1015
0.4131-j0.3510
0.9901-j0.1479
1.2973-j0.3797
1.0083-j0.0185
4
0.8217-j0.3644
0.9986-j0.0880
表4.1迭代过程中节点电压变化情况
节点不平衡量
-2
-1
-0.1482
-0.9769
-0.0726
-0.0103
-0.0902
-0.6071
-0.0480
-0.0022
-0.6272
-4.3251
-0.3610
-0.0171
-0.1816
-1.2510
-0.1042
-0.0049
表4.2迭代过程中节点不平衡量变化情况
结果值与我的小组同学基本一样,也在预期之内。
得到了基本一致的结果。
并且确定牛顿法具有很好的二次收敛性,是求解多元非线性方程的正确算法。
5总结
这次的电力系统分析课程设计让我对平时所学的专业知识有了更深刻更具体的了解,明白了理论知识必须与实践相结合才能更好的发挥作用。
在不停的翻书上网查资料的过程中,我积累了大量的导纳矩阵和潮流计算以及电力系统的知识,全面透彻的了解了相关知识的应用。
使自己的知识更加牢固,并且有了更深的理解。
通过这次毕业设计,我才明白学习是一个长期积累的过程,在以后的工作、生活中都应该不断的学习,努力提高自己的知识和综合素质。
总之,万事开头难,知识必须通过应用才能实现其价值!
有些东西以为学会了,但真正到
用的时候才发现是两回事,所以我认为只有到真正会用的时候才是真的学会
参考资料
[1]何仰赞温増银.电力系统分析(第三版).华中科技大学出版社
[2]吴国炎.电力系统分析.浙江大学出版社
[3]华智明岳湖山.电力系统稳态计算.重庆大学出版社