电力系统综合负荷模型的辨识算法研究

电力系统综合负荷模型的辨识算法研究Study on A ggregate Pow er L oad Param eter E sti m ati on A lgo rithm
长沙电力学院(湖南长沙410077) 王　进　苏　盛
湖南大学(湖南长沙410082) 李欣然
摘　要:研究了一种具有全电压范围适应性的动态负荷模型,提出了一种用于模型参数辨识的改进模式搜索算法,仿真结果表明了所提出的模型辨识算法的可行性和有效性。

关键词:电力系统综合负荷模型;动态负荷模型;模型参数辨识;改进模式搜索算法
中图分类号:TM715
文献标识码:A
文章编号:100329171(2002)1120008205
电力系统运行状况的仿真计算,是电力系统规划设计和运行决策的重要依据。

电力综合负荷建模在电力系统仿真计算中占有极为重要的地位。

不正确或不恰当的负荷模型,难以描述负荷的真实运行行为,使电力系统仿真计算结果的可信程度大为下降,甚至可能导致决策的错误。

如参考文献[1]所述,在对加拿大渥太华地区一个局部系统从互联大系统解列后的系统电压稳定研究中发现,采用动态负荷模型与静态负荷模型得出的结果差别巨大,静态模型甚至导致完全错误的结论。

这说明电力负荷的运行特性或运行行为的描述对电力系统电压稳定性有着极其重要的影响,动态负荷模型对电压稳定分析的重要性越来越受到广泛的重视。

本文在研究一种具有全电压范围性的综合负荷模型的基础上,提出了一种改进模式搜索优化算法,详细阐述了模型的参数辨识原理及算法实现过程。

通过将综合负荷模型用于动模实验的成功建模,得出了一些有意义的结论。

1　系统模型[2]
111　综合负荷等值模型
设综合负荷由静态负荷与动态负荷两部分组成。

静态部分用一恒定导纳(Y s=G s-j B s)等值;动态部分(Y d=G d-j B d)的等值模型由式(1)～式(5)描述。

负荷电纳B>0时为感性。

G1(V)=k g0V k g1,ΠV≥V cr
G2(V)=g0,ΠV<V cr
(1)
B(V)=b0+b12-(G(V) b1)2
1-1-(G(V) b1)2+
b2G b3(V),ΠV≥V cr
B(V)=b0+b12-(G(V) b1)2
1+1-(G(V) b1)2+
b2G b3(V),ΠV<V cr
(2)
T g(d G d t)=f G1(V,G)・
f G2(G)=f G(V,G)
(3)
f G1(V,G)=P ls-P ld=P ls-V2G(t)(4)
f G2(G)=+K G21-(G G cr)2
1-1-(G G cr)2
f G2(G)=-K G21-(G G cr)2
1+1+(G G cr)2
(5)
(1)～(5)式中各有关部分所代表的意义如下[1]:
(1)式(1)中,V为综合负荷的运行电压,G (V)为静态等值电导;G1、G2分别对应于稳定区和失稳定区,V cr是其失稳临界电压,系数k g0,k g1为稳定区或稳定模式下的模型参数;g0为失稳区或失稳模式下的模型参数。

(2)式(2)中,B(V)是综合负荷的静态等值电纳,且其中之上式对应于稳定区或稳定模式,下式对应于失稳区或失稳模式;b0、b1、b2、b3为模型参数;G(V)由(1)式确定。

分析式(2)可以看出,G 满足关系:0<G≤b1,因此,由综合负荷等值导纳的变化规律即知,模型参数b1正是综合负荷等值临界电导G cr,即:b1=G cr。

(3)式(3)～式(5)中,T g为综合负荷的惯性时间常数,G cr和K是模型参数;f G1描述了负荷与
系统间的功率平衡关系,其中P ls是负荷为完成其能量转换功能而要求从电网吸收的功率,P ld是综合负荷在受扰动后实际能够从电网得到的功率;
f G2描述了综合负荷电导对时间的变化率与电导自身间的关系,其中式(5)的上、下两式分别对应于负荷的稳定区和失稳区,它所描述的就是负荷
的稳定行为模式(稳定模式和失稳模式)。

112　系统状态方程
待辨识系统为一阶系统,以动态等值电导为状态变量,其状态方程即式(3),且设其中P ls= P ls0恒定。

动态等值电纳则由(2)式确定。

113　系统输出方程
以综合负荷从电网吸收的总功率为输出变量,负荷电压为输入变量,系统输出方程如式(6)所示。

P=V2(G s+G d)
Q=V2(B s+B d)
(6) 114　系统稳态方程
由式(1)、(3)和(4),系统受扰动前的稳态条件下,将有
P ls0-V20G d0=0,G d0=k g0V k g10(7)由(2)式
B d0=b0+b12-(G d0 b1)21-1-(G
d0 b1)2
+
b2G b3d0;b1=G cr(8)由(6)式
P0=V20(G s+G d0)
Q0=V20(B s+B d0)
(9)
由此解出
G s=(P0 V20)-G d0
B s=(Q0 V20)-B d0
(10)
式(7)、(8)、(10)即构成辨识模型的稳态方程。

2　模型辨识原理
211　辨识原理框图[3,4,5]
本文采用如图1所示的辨识原理。

以系统输入向量u(t)的时间序列u(0)、u(1)、…、u(N)作为被辨识的实际系统和等效模型系统的输入;同时将实际系统在t=0时刻的输出y(0)反馈给模型系统作为其附加输入,以保证模型输出满足实际系统的初始稳态条件。

图1所对应的辨识过程分两步进行
:
(1)参数辨识
开关K2之021接通,022断开,同时开关K3接通,辨识模型参数,确定独立待辨识参数向量Α,并根据系统的初始稳态条件确定非独立的导出辨识参数向量Β。

因为模型响应不仅取决于模型参数及系统激励,而且也与综合负荷的初始负荷功率及电压水平有关。

本文把初始稳态条件称为初始等效条件。

与此相对应,图1中之开关K1当且仅当t=0时闭合。

图1　辨识原理示意图
(2)模型检验2动态回响测试
开关K2之021断开,022接通,同时开关K3断开,由系统输入量u(t)的时间序列求解模型响应。

这一过程也相当于综合负荷建模时的系统动态回响测试[5]。

212　辨识准则
以系统实测响应与模型响应之差的平方和最小作为辨识准则。

目标函数如式(11)所示。

m in J[x(t),u(t),Α,Β]=m in6k=N k=0[y(k)-
y m(k)]T[y(k)-y m(k)](11)上式中,k是实测响应的采样时刻,N为系统实测数据的数据长度;x(t)=G d为状态向量;u(t) =V(t)为输入向量;y为系统实测响应或系统输出响应向量,可以取为综合负荷从电网吸收的总功率或总负荷电流,本文中一律取综合负荷功率P、Q,即有:y=[P　Q]T,它们由现场或动模实验的实测给出;y m是与系统实测响应(即y)相对应的模型输出响应向量,由输出方程确定。

待辨识的模型参数为:
(1)Α——独立的待辨识参数向量,且Α= [k g0k g1b0b1b2b3T g K]T,且b1=G cr;
(2)Β——非独立的导出待辨识参数向量,由稳态方程解出,且Β=[P ls0G s B s]T。

3　辨识算法
311　模式搜索法(Hook-Jeeves法)[6]
动态辨识本质上是在满足一定约束条件下的
参数优化问题。

用于负荷模型辨识的数学优化方法很多,本文采用模式搜索直接优化算法(Hook2 Jeeves法)。

这种方法的每一步迭代中,由两种移动组成,一种称为探测移动,目的在于探测一个参考点的周围是否有更好的点,从而找出下降方向;另一种称为模式移动,目的在于沿有利的下降方向进行加速。

两种移动交替进行,逐步逼近最优点。

该算法无需计算目标函数导数,虽然计算速度较慢,但数值稳定性好,使用方便可靠,加之此处是应用于一阶模型,且待辨识参数少,实践证明其计算速度对于离线辨识来说是完全可以接受的。

微分方程求解采用四阶龙格库塔法,这是一种普遍使用的求解微分方程的数值积分算法。

312　算法的改进
在负荷建模过程中,笔者发现有时目标函数的变化区域内存在着一个狭长的区域——深沟。

利用模式搜索法往往会跨过这个区域,而找不到这个狭长区域的走向,从而无法搜索到最优解。

针对算法用于负荷测辨时存在的不足之处,笔者提出了一种改进模型搜索法,该算法是在模式搜索方法的基础上对步长加以改进而形成的,具体过程如下:
(1)进行探测搜索时,步长S除初始时给定以外,后面的步长依赖于前面2个基点的距离。

设前2个基点为b k-1,b k,则令步长为:
S=Ξ～‖b k-b k-1‖
其中,Ξ为步长收缩因子,0<Ξ<1,一般取011～015,参考点为b。

(2)在探测搜索结束后,设基点为b(k)0,参考点为b(k),然后取基点r=b(k)+Α(b(k)-b(k)0)(Α为步长加速因子,1<Α<2)开始进行模式移动,模式移动失败时,常规的方法是取前次搜索的终点为参考点开始下一次探测搜索过程。

笔者在使用该算法时,对此进行了改进,模式移动失败时,延伸方向并不立即放弃,而是取一个数Β(0<Β<1)令
r1=b(k)+Β(b(k)-b(k)0),r2=b(k)-Β(b(k)-b(k)0),先后在r1和r2处以缩小的步长做新的延伸,若成功,则以此次延伸的终点作为参考点开始新的探测过程;若均失败,则按常规方法进行下一次探测,求新的延伸方向。

对常规方法的改进避免了因模式移动的速度过快而越过了目标函数变化的狭长区域,从而保证搜索方向的可行性。

(3)常规算法在给定初始值后即开始探测搜索,因而对初始值的选定要求较高,改进算法中,在给定初始值后首先进行直线搜索获取最优初始值,将其作为参考点进行探测搜索,从而降低了算法对初始值的敏感度。

4　辨识实例
411　实测数据来源
实测数据取自在清华大学国家重点实验室动模实验所获得的实测动态综合负荷记录数据[2],辨识所用记录对应的负荷构成相同(带机械负载的感应电动机、电阻电炉、白炽灯照明负荷),电压激励扰动形式相同(阶跃下降而后恢复),但激励强度(电压下降幅度和持续时间)不同,共使用该类动模实验动态负荷数据记录9个。

辨识算法采用上述第3节中提出的改进模式搜索法。

412　辨识结果
为了考察模型对综合负荷的描述能力,笔者对9个记录均进行了动态辨识。

这里给出其中3个记录的部分辨识结果。

分别如表1和图2～图5所示。

表1中,第一列是负荷记录编号,第二列为电压扰动幅度,其它各列为辨识所得的模型参数。

我们把在某一记录对应的激励下辨识所得之模型参数对于本记录的建模激励下的模型响应,称为自辨识响应;否则为内插或外推模型响应。

在图2(a)、(b)及图3(a)中,分别给出了2条记录使用常规算法辨识的有功及无功的实测及模型响应曲线比较。

在图2(a)、(b)及图3(a)、(b)中,分别给出了2条记录使用改进算法辨识的有功及无功的实测及模型响应曲线比较。

其中,在响应曲线比较图中,虚线代表模型的实测曲线,实线代表模型响应曲线。

在图4(a)、(b)及图5(a)中,分别给出了使用常规算法辨识的有功和无功的实测、外推及内插模型响应曲线比较,在图4(a)、(b)及
表1　模型对清华动模实测负荷记录的辨识结果(标么参数:S B=P0,V B=380V)
序号∃V%k g0k g1b0b1(G cr)b2b3T g K P ls0G s B s 145013543-21644051961819195300134140144533011007018040013628016858-514622 267014346014942-1415878151382901279711649026108792171770138660167391512420 377014240-415760-1415827161024001411301363546117872172400144120157631419856
图2　表1之记录1的有功实测及模型响应(V0=0196,∃V=45%)图3　表1之记录2的无功实测及模型响应(V0=0198,∃V=77%)图4　表1中记录2对记录1的有功实测及外推模型响应比较图5　表1中记录2对记录3的内插无功实测及内插模型响应比较
图5(a)、(b)中,分别给出了使用改进算法辨识的有功和无功的实测、外推及内插模型响应曲线比较,其中,粗实线代表实测响应,细实线代表自辨识响应,虚线代表外推或内插响应。

5　结果分析及其评价
综合上述辨识结果,可对综合负荷模型的描述能力及相应的改进算法做出如下基本评价:
(1)模型揭示了综合负荷G、B,进而P、Q之间的内在联系,从而实现了对有功功率和无功功率的综合描述,这是除感应电动机机理模型外的其他非机理模型所不具备的优点。

尤其当无功功率很小甚至为负时,一般非机理模型很难描述[7],而本文所采用的模型则具有较好的描述能力(如图3、5)。

(2)由图2～图5可知,模型响应与实测响应间的稳态误差很小,虽然存在一定的暂态误差,但就其总体效果而言,模型对有功和无功功率均具有较强的综合描述能力。

但也应注意,模型对无功功率的逼近能力次于有功功率。

这主要是因为模型是一阶的,相当于忽略了动态负荷成份的电磁暂态过程(只计及机械暂态)。

但由于模型的稳态逼近效果好,故在一般情况下,使用本文模型不会改变电压稳定性分析的基本结论。

(3)模型对有功功率具有较强的内插和外推能力[5]。

其中图4中自辨识响应曲线(电压扰动激励67%),与外推响应曲线(建模激励扰动幅度45%)以及图5中自辨识响应曲线(电压扰动激励67%)与内插响应曲线(建模激励扰动幅度77%)相比,均有较好的逼近效果。

(4)从图2～图5可明显看出,改进算法无论在模型的自辨识响应,或是模型的外推和内插能力描述上,均有不同程度的改善,因而有力地证明了该算法的有效性和可行性。

参考文献
[1]Eb rah i m V aahedi,et al.D ynam ic L oad M odeling in L arge Scale Stab ility Studies.IEEE T ran s.on Pow er System s,1983,3(3):1039～1045
[2]李欣然.考虑负荷特性的电压稳定及其负荷建模的研究[博士论文].北京:华北电力大学,1998
[3]鞠平,马大强.电力系统负荷建模.北京:水利电力出版社,199516
[4]沈善德.电力系统辨识.北京:清华大学出版社, 1993112
[5]贺仁睦,魏孝铭等.电力负荷动特性实测建模的外推和内插.中国电机工程学报,1996,16(3)
[6]薛嘉庆.最优化原理与方法.北京:冶金工业出版社, 1983
[7]张伶俐.电力负荷动特性综合的研究[硕士论文].北京:华北电力大学,1997
(收稿日期:2002204220)
下　期
　要　目　
京津唐东部地区电网安全稳定分析
电铁供电站高次谐波滤波兼无功补偿装置
　的设计
一种改进的三相电路谐波和负序电流检测
　方法研究
超导储能装置变流器控制策略的分析
滦河电厂5号主变色谱异常分析及局部放
　电试验
新、老《火电厂污染物排放标准》的对比与探讨
提高电除尘器除尘效率的研讨
在线监测介质损耗因数的影响因素与结果
　分析
20C r1M o1V T i B螺栓断裂原因分析
基于一种新方法的煤灰结渣预测
四角切圆锅炉浓淡燃烧器的应用研究
大同二电厂一次调频的实现方式
关于俄制发电机定子接地保护的改造
下花园发电厂3号炉排粉机改造
秦电2号炉汽包焊缝缺陷的检验与处理
220kV开关非全相故障原因分析
滦河发电厂汽轮发电机组“风温逆转”
风机主要调速节能技术的比较。