电力系统潮流计算.docx
《电力系统潮流计算.docx》由会员分享,可在线阅读,更多相关《电力系统潮流计算.docx(33页珍藏版)》请在冰豆网上搜索。
电力系统潮流计算
第四章电力系统潮流分析与计算
电力系统潮流计算是电力系统稳态运行分析与控制的基础,同时也是安全性分析、稳定性分析电磁暂态分析的基础(稳定性分析和电磁暂态分析需要首先计算初始状态,而初始状态需要进行潮流计算)。
其根本任务是根据给定的运行参数,例如节点的注入功率,计算电网各个节点的电压、相角以及各个支路的有功功率和无功功率的分布及损耗。
潮流计算的本质是求解节点功率方程,系统的节点功率方程是节点电压方程乘以节点电压构成的。
要想计算各个支路的功率潮流,首先根据节点的注入功率计算节点电压,即求解节点功率方程。
节点功率方程是一组高维的非线性代数方程,需要借助数字迭代的计算方法来完成。
简单辐射型网络和环形网络的潮流估算是以单支路的潮流计算为基础的。
本章主要介绍电力系统的节点功率方程的形成,潮流计算的数值计算方法,包括高斯迭代法、牛顿拉夫逊法以及PQ解藕法等。
介绍单电源辐射型网络和双端电源环形网络的潮流估算方法。
4-1潮流计算方程--节点功率方程
1.支路潮流
所谓潮流计算就是计算电力系统的功率在各个支路的分布、各个支路的功率损耗以及
各个节点的电压和各个支路的电压损耗。
由于电力系统可以用等值电路来模拟,从本质上
说,电力系统的潮流计算首先是根据各个节点的注入功率求解电力系统各个节点的电压,当各个节点的电压相量已知时,就很容易计算出各个支路的功率损耗和功率分布。
假设支路的两个节点分别为k和丨,支路导纳为ykl,两个节点的电压已知,分别为
V&和V&,如图4-1所示。
图4-1支路功率及其分布
那么从节点k流向节点I的复功率为(变量上面的“一”表示复共扼)
从节点I流向节点k的复功率为:
功率损耗为:
(4-3)
Si二&iSk=(V&-V&WMF)=-yki£
因此,潮流计算的第一步是求解节点的电压和相位,根据电路理论,可以采用节点导纳方程求解各个节点的电压。
2.节点功率方程
根据电路理论,要想求系统各个节点的电压,需要利用系统的节点导纳方程。
I刃
*
t1
L:
f
1测
尸
V
图4-2电网络示意图
如图4-2所示的电网络,有N个节点,假如已知各个节点的注入电流源的电流,以及各个支路的支路导纳,那么可以根据节点导纳方程求出电网各个节点的电压:
YV=
Is
(4-4)
其中
Yii
Yl2
L
YiN-
Y
Y22
21
L
Y2N
Y=
L
L
L
L
Yni
Yn2
L
Ynn
为电网络的节点导纳矩阵,Ykk(k=1,2,LN)为自导纳,是与k节点所有连接支路
导纳之和,Ykl(k=1)为互导纳,等于负的连接k和|节点的所有支路导纳之和。
V=[Vi,V2,L,Vn]T为各个节点的电压相量,Is=[lsi,ls2,L,Isn]T为注入到各个节点的总电流。
2.1节点复功率方程
要想计算各个节点电压,除了需要知道系统参数及节点导纳矩阵以外,还需要知道节点的注入电流源的电流。
然而电力系统中,节点的注入电流是不知道的,已知的是各个节点的注入功率。
这就需要将节点电压方程转化为节点功率方程。
方程4-4中第k(k=1,2,L,N)个节点的方程可以写作:
N
、YklV&二YdV&Yk2V&2LYkkV&kLYrnV&n=l(4-5)
&Sk
在方程4-5两端乘以Vk,得到:
(4-6)
Vk"YkiV&iRkl&sk三Ssk=Psk一jQsk
假如在电力系统中,各个节点的注入复功率都已知,那么就可以用方程4-6组成的方
程组求解各个节点的电压。
然而实际情况并非如此,已知的条件是:
有的节点的注入复功率S是已知的,有的节点的电压幅值和注入有功功率是已知的,有的节点的电压和相角是已知的。
根据这三种不同的情况,电力系统中各个节点分为三种类型:
PQ节点、PV节点
和V节点。
所谓PQ节点,就是该节点的注入复功率S是已知的,这样的节点一般为中间节点或者是负荷节点。
PV节点,指该节点已知的条件是注入节点的有功功率P和该节点的电压幅值V,这
样的节点通常是发电机节点。
V■:
.节点指的是该节点的电压幅值和相角是已知的,这样的节点通常是平衡节点,在每个局部电网中只有一个这样的节点。
当然,PQ节点和PV节点在一定条件下还可以互相转化,例如,当发电机节点无法维持该节点电压时,发电机运行于功率极限时,发电机节点的有功和无功变成了已知量,而电压幅值则未知,此时,该节点由PV节点转化为PQ节点。
再比如某个负荷节点,运行要
求电压不能越限,当该节点的电压幅值处于极限位置,或者电力系统调压要求该节点的电压恒定,此时该负荷节点就由PQ节点转化为PV节点。
假如全系统有N个节点,其中有M个PQ节点,N-M-1个PV节点,1个平衡节点,每个节点有四个参数:
电压幅值V、相位角、:
(用极坐标表示电压,如果用直角坐标表示
电压相量则是e和f)注入有功功率Ps和无功功率Qs,任何一个节点的四个参数中总有两个是已知的,因此N个节点,有2N个未知变量,N个复数方程(即2N个实数方程,实部和虚部各一个),通过解这个复数方程就可得到另外2N个参数。
这就是潮流计算的本质。
但在实际求解过程中,由于我们求解的对象是电压,因此,实际上不需要2N个功率
方程,对于M个PQ节点,有2M个功率方程(M个实部有功功率方程,M个虚部无功功率方程);对于N-M-1个PV节点,由于电压有效值V已知,因此只有N-M-1个有功功率方程;对于平衡节点,由于电压和相角已知,不需要功率方程。
因此总计有2M+N-M-
仁N+M-1个功率方程。
如果电压相量用极坐标表示,即V&=VkE£:
k,则M个PQ节点有
2M个未知数(M个电压有效值,M个电压相角),N-M-1个PV节点有N-M-1个未知数
(电压有效值已知,未知数为电压相角),平衡节点没有未知数,因此未知数的个数也是N+M-1个,与方程数一致。
如果复电压用直角坐标表示,V&=ek亠jfk,则有2(N-1)个未
知数,还需要增加N-M-1个电压方程,即22。
kk/feV
+=
2.2用直角坐标表示的电力系统节点功率方程
对于PQ节点,已知的是注入节点的功率p和Q,将Ykm=Gkm-jBkm和V&=ekjfk带入节点功率方程的复数表示式中,可以得到有功功率和无功功率两个方程:
NJN1
Psk二PGk-PLk二ek=(Gkmem-Bkmfm)"fk*(Gkmfm'BkmOm)
(4-7)
mAm二
-N_1N_1
Qsk二QGk-QLk二fk'(Gkmem-Bkmfm)-e^(Gkmfm'Bkmem)
mAm二
上式中psk和QSk为注入到节点k的净功率,即注入和消耗的代数和。
PGk、QGk表示
注入的功率,pLk和QLk为消耗的功率。
对于PV节点,除了有功功率方程外,因为已知该节点的电压幅值,还有一个电压方
程:
222
Vkek■fk(4-8)
方程4-7可以抽象的表示为:
APk(ei,fi,L丄f”」)=0
丿(4-9)
人Qk(e1,f1,L,3丄fN」)=0
方程4-8可以抽象的表示为
Vk(e1,f1,L,eN」,fN」)=0(4-10)
因此,对于一个具有N个节点的电力系统,其中M个PQ节点,N-M-1个PV节点,1个平衡节点,有方程如下:
也P(e1,f1,L,eN」,fN」)=0
.1
-Q1(e1,f1,L,eN」,f”」)=0
LLLLLL>2M个PQ节点的方程
匚PmG,£丄◎丄f”v)=0
lQmG,f1,L,eN」,fn」)=0
八Pm1(e1,f1,L,eNv,fN4)=0
八Vm1(e1,f1,L,eN/,fN4)=0
lVnj(e1,f1,L,eN4,fnj)
N个节点,平衡节点的电压幅值和相角已知,即其横分量和纵分量已知,因此平衡节点不参与计算。
N-1个节点的电压的横分量和纵分量为未知数,共2N-2个未知数。
2M个
PQ节点方程,2(N-M-1)个PV节点方程,共计2N-2个方程。
解这个方程组,就可以得到电力系统N个节点的电压相量,根据各个节点的电压相量
和已知的注入功率,就可以计算出各个支路的潮流分布,及各个支路的功率损耗。
2.3极坐标表示的节点功率方程
对于PQ节点,已知的是注入节点的功率p和Q,将Ykm二GkmjBkm和
V&二Vk/l:
k带入节点功率方程的复数表示式中,可以得到实部和虚部两个方程:
N-.-
(4-12)
Psk二PGk-PLk二Vk=Vm(GkmCOSkmBkmSin门km)
Qsk=QGk-QLk=VkgVm(GkmSin§km—BkmCOS§km)
m-1
上式中,V代表电压幅值,亦二、*-5。
对于PV节点,由于节点的电压幅值已知,因此只有有功功率方程而没有无功功率方程。
同样,方程4-12可以抽象的表示为:
1
-Pk(V1L,Vm,1,L,;!
n」)=0(4-13a)
Qk(V1L,Vm,、1,L,、n」)=0(4-13b)
lQm(V1,L,Vm,--1,L,':
n4)0
因此,对于一个具有N个节点的电力系统,其中M个PQ节点,N-M-1个PV节点,1个平衡节点,有方程如下:
P(V1,L,Vm,1,L,:
nj^0N_M-1个PV节点方程
Q(V1,L,Vm,1,L=0
LLLLLL眩M个PQ节点方程
:
Pm(V1,L,Vm,、1,L
J
匚Pm1(V1,L,Vm,--1,L,':
n4)=0
LLLLLL(4-14)
•:
Pn-1(V1,L,Vm,、1,L,、nG=0
1
除了平衡节点外,N-1个节点中,有W彳!
个QT功点的率方幅值和相角都是未知数,N-
M-1个PV节点的相角为未知数,因此共有2M+N-M-仁N+M-1个未知数,2M+N-M-
仁N+M-1个方程。
在方程4-14中,可以把N-1个有功功率方程放在一起,M个无功功率方程放在一起:
P(V1,L,Vm,1,L,、N4)=0
LLLLLL
Pn4(V1,L,Vm,1,L,':
Nd)=0|l
Q(V1,L,Vm,1,L,、nv)=0Qm(V1,L,Vm,
LLLLLL:
M个无功功率方程61,L,6N_)=0U
(4-15)解上述方程组,就可以得到电力系统中各个节点的电压幅值和相角,进而可以计算出各个支路的潮流分布和损耗。
3.小结
潮流计算是计算电力网各个支路的功率潮流分布和功率损耗,同时也计算各个支路的电压损耗。
首先要求电力网各个节点的电压相量。
根据电网络理论,节点电压通常采用节点导纳方程来求解,即已知电网络的节点导纳矩阵和各个节点的注入电流源的电流,求解节点导纳方程。
然而通常电力系统各个节点的注入电流是未知的,已知的是各个节点的注入功率,因此需要将节点电压方程转化为节点功率方程。
实际电力系统的节点注入功率并非都已知,有的已知注入有功功率P和无功功率Q称
为PQ节点;有的已知注入有功功率P和节点电压有效值V,称为PV节点;有的已知节
点电压V和相角d,称为平衡节点或V、•节点。
无论哪种类型节点,每一个节点均含有4个
参量P、Q、V、、(或e、f)已知的是其中的两个,故而可以利用节点功率方程(4-6)求
解出另外两个参量。
假设系统有N个节点,必然有2N个未知数,同样有2N个节点功率
方程(4-17中的实部和虚部各一个)。
实际上,我们求解的目标是电压,对于PV节点和V、•节点来说,前者电压有效值已知,
后者电压相量已知,因此不存在2N个未知数,当然也不需要2N个方程。
假设系统有N
个节点,M个PQ节点,1个平衡节点,对于直角坐标表示的节点电压来说,有2(N-1)个
未知数,2M+N-M-1个功率方程,只需要再补充N-M-1个电压方程就可以了;对于极坐标表示的电压来说,只有N-1个「未知数,M个V的未知数,因此只需要N+M-1个功率方程就足够了。
无论怎样,潮流计算是解决这样的一组非线性代数方程组:
F(X,C,U)=0(4-16)
其中,X代表系统状态,包括电压V和相角J.;C代表参数,包括电导G和电纳B;U表示系统激励,即注入的功率。
求解这样的多维非线性代数方程组,需要利用计算机进行辅助迭代计算,即先给定一个初值,然后不断迭代,逼近真实解。
方法有:
高斯-赛德尔迭代法,牛顿-拉夫逊法和PQ
解耦法。
4-2高斯—赛德尔叠代法
1.基本原理
为了方便理解这个n维方程组的叠代求解方法,先从一元非线性方程的求解开始。
假设有一维方程f(x)=0,高斯法的基本原理是,先将方程转化为:
x=g(x)(4-17)
那么给定一个初值x[0],代入就可以得到一个新值x⑴二g(x[0]),第k次叠代的值为:
(4-18)
(4-19)
[k1]/叫
x^g(x)
一直叠代到误差满足要求为止,即
[N][N-4]
x-X
这个解方程的方法称为高斯叠代法。
这个叠代求解的过程可以这样来理解:
x=g(x)的解可以认为是两个曲线y=x和讨二g(x)的交点的横坐标x,首先给定一个
初值x[0],g(x[0])与斜线y=x的交点的横坐标即为叠代后的新解x[1],g(x⑴)与斜线
y=x的交点的横坐标即为叠代后的新解x[2],如此围绕交点往复循环,不断地逼近方程
的解,如图4-4所示。
n为方程组为:
(4-20)
(4-21)
高斯迭代法可以推广到n维非线性代数方程组,假设fl(xi,x2L,xn)=0
f2(xi,x2L,Xn)=0
LLLLLL
fn(Xi,x2L,Xn)=0
首先将方程组4-20转化为:
Xi=g(Xi,x2,L,xn)
刈=g(Xi,x2,L,Xn)
」LLLLLL
M=g(Xi,X2,L,Xn)
给定一组初始值X[0]=[』,』丄,Xn0]]T,带入上式,得到一组新值x[1]=g(X[0]),
不断叠代,循环往复,第k次叠代为:
(4-22)
X[k1]二g(X[k])
其中第j个方程为
[k1]/k]k]k]
xgj(xi,X2,L,Xn)(4-23)
[[[
直到叠代前后的解的最大误差不超过允许的误差为止,即
[N1][N]|,、
max{xj-Xj}(4-24)
为了提高高斯叠代法的收敛速度,赛德尔提出将已经叠代出的新值代替旧值参与叠代计算,如在第k次叠代中,第j个方程为
[k1]/k1][1][k]k]、
X=gj(Xi,L,x_,X,L,Xn)(4-25)
[k[
第1至j-1个元素已经叠代出k+1次的值,因此代替第k次的值参与第j个元素的叠代,就可以提高收敛速度。
2.电力系统潮流计算的高斯-赛德尔迭代法
电力系统潮流计算需要求解节点功率方程,其中第m(m=1,2,…,N个节点功率方
程为:
Np2Nq
(4-26)
(4-27)
Vm'YmlV&1二YmmVmfm、YmV&I=PSm—jQsm14l=1
如上式变换为x=g(x)的形式,可以得到如下的方程:
\/&_1(PSm]jQSm_、NYmlV&l)
mYmmVmy
lrm
根据高斯-赛德尔迭代法,首先给定电压相量的初值,对于PQ节点,不仅需要给定
电压幅值的初值,还要给出相角的初值(设为零)。
假如第m号节点为PQ节点,第k次叠代公式为(第m个节点以前的节点第k次叠代已经完毕,因此用k+1次的值取代k次的值,而在第m个节点以后的节点尚未进行第k
次叠代):
对于PV节点,给定的初值的电压幅值为给定的电压,相角初值设为零。
可是对于PV
节点来说,注入该节点的无功功率未知,因此第k次叠代时,首先按照下式计算注入PV
节点(假设第m个节点是PV节点)的无功功率:
k[]&[klk[1&国N【k_1]
Qsmlm[VmISm]=lm[Vm(YmlVlYmlVl(4-29)
[k])]l=1lorn
如果在叠代计算过程中,任意节点的电压和无功功率必须满足不等约束条件:
如果在叠代过程中,PQ节点的电压幅值超出允许的范围,则该节点的电压幅值就固
定为允许电压的上限(如果超出上限)或下限(如果越过下限),PQ节点就变为PV节点
继续进行叠代。
同样,对于PV节点来说,如果在叠代过程中,无功功率Q超出了允许的
范围,则PV节点就变为PQ节点继续参与叠代。
高斯-赛德尔叠代法的计算过程如下:
(1)第一步:
设置初始值,对于PQ节点,由于其电压相量的幅值和相角都未知,因
此初始的电压相量的幅值可以设定为各个点的额定电压,相角选择为零;对于PV节点,
由于其电压相量的幅值已知,因此幅值用已知的设定电压,初始相角设定为零。
(2)第二步:
对于PQ节点,直接将设定的初始值代入,用4-28求得下一次迭代的电压值,然后判断是否电压越限,如果越限,则用其限值(越过上限用上限值,越过下限则用下限值),该节点在下一次迭代过程中转化为PV节点;对于PV节点,则首先利用式4-29求出注入的无功功率,然后校验无功功率是否越限,如果越限则采用上限值或者下限
值,下一次迭代时该节点转化为PQ节点,将求得的注入无功功率和已知的有功功率代入
4-28求解下一次迭代的电压相量值。
(3)第三步:
判断误差是否满足要求,用第k次迭代的结果和k-1次迭代的结果进行比较,如果其最大的误差满足事先设定的误差要求,则输出计算结果,如果不满足要求,
则返回第二步继续迭代。
其计算流程图如图4-5所示。
没宜初值
訂PQ召点?
卜理一]计算无功Q(k)
*
Q越限?
f
Q駆限值,菸为PQ节点
计算电叭gv.k)
i
"取限值.转为PVP点
亠
T赵遍¥
计算电
图4-5高斯赛德尔迭代法求解电力系统潮流的计算流程图
4-3牛顿—拉夫逊法
因此,可以得到这个一元非线性方程的求解步骤为:
首先给定解的初值x[0],然后根
据公式4-32求出初始值的修正值.-;x[0],由此可以得到该方程的新的解X[1]=X[0]•AX®,如此反复叠代,直到误差满足要求Ax【N]|£g。
迭代计算流程如图4-6所示。
楠出结END
图4-6牛顿拉夫逊法计算流程
其迭代求解过程的几何意义如图4-7所示。
图4-7牛顿拉夫逊法的几何解释
可以把上述求解一元非线性代数方程的方法推广到n维非线性代数方程(如4-20)的
求解。
非线性代数方程组4-20可以表示为矩阵形式:
F(X)=0
同样假定X0是该方程组的近似解,与真实解之间的误差为
勒级数:
(4-33)
X,在X0处展开一阶泰
F(X0:
X):
F(X0)JX=0
(4-34)
其中:
L
cX
CX2
1
Cf2
L
依1
CX2
L
L
L
Cfn
L
01
CX2
:
fl
£Xn
cf2
-Xn
L
Cfn
€Xn_x_x0
(4-35)
XJ‘F(X0)
(4-36)
被称为雅克比矩阵。
4-34称为修正方程,修正方程可得到修正值X:
计算过程与一维方程的牛顿法求解类似,首先给定初值X[0]=[X1[0],L,Xn0]]T,
并计算出在初始值处的雅克比矩阵J0,利用4-36式计算初始值的修正值
X[0]二J:
F(X[0]),根据这个差值可以得到修正后的解X⑴二X[0]•厶X[0]。
如此循环
往复,在第k次叠代时,计算雅克比矩阵Jk,根据4-34计算修正值X[k]=-J:
F(X[k]),
得到第k+1次修正后的解:
x[k1]=x[k]•厶x[k],重复上述过程,直到误差满足要求为止。
可见,牛顿拉夫逊法的关键在于求解雅克比矩阵J,由于直角坐标表示和极坐标表示
电压相量的节点功率方程有所不同,因此其雅克比矩阵也有很大的差异。
(4-37)
2.直角坐标节点功率方程的牛顿-拉夫逊法
仍然假设系统有N个节点,其中M个PQ节点,N-M-1个PV节点,1个平衡节点。
则M个PQ节点方程为(假设1号节点至M号节点为PQ节点):
NN
-PSk-ekV(Gklel-Bkifi)-fkX
k-1,2,L,M。
(4-38)
\7
222
-Vk-Vk-^e-"fk0
1
1■_
1
1
Lm
L
Mm
Lm|
1
Mm1
Lmn1
L
Mn1
II
「—-
1佔
L
Nrrtm
1
Hkfm|
Hkltm1
L
Ntn!
II
%]
Rh11
Sn
L
Rmm
Srr|m|
IRr1m1
Stpi!
L
FTn!
^nlfr
1
LL
L
L|LLL
L
L
1
1
H11
L
Nm
HUm|
N1m1
H^1m1
L
Nt!
11
1
1
Snn
L
FR1m
S1m|
S1dr
L
FR4n1
11
(4-39)
其中:
Nkj=「:
Pk--Gkjek-Bkjfk(j-k)
:
:
Pk--Gkkek-Bkkfk-'(Gklel-Bklfl),:
eki-1
nA
=Bkkek—Gkkfk-、(Gkifi亠Bkie)l丄
n」
=Bkkek-Gkkfk+迟(Gklfi+Bkie)
i=1
Lkj
;:
:
Qk
f
=BkjfkGkjek
Lkk
■Ql=GkkekBkkfk
fk
n.1
(Gkiei
id
-Bkifi)
Rkj
-k)
_ej
基于直角坐标的牛顿-拉夫逊法求解潮流计算的步骤如下:
(1)第一步:
设定初值,对于PQ节点,其电压幅值的初值设定为