发明内容
本发明的目的在于提供一种电网频率动态过程仿真轨迹的校正方法。本发明在电网模型和参数准确条件下,以事故下系统真实的频率轨迹作为参照依据,对调速器参数在其可行的范围内进行调整,使得仿真轨迹逼近真实的频率轨迹;建立误差评价模型,通过轨迹误差评价函数,来评价频率轨迹的逼近程度;通过禁忌搜索,来调整相关参数,从而实现对频率仿真轨迹的有效校正。
为实现上述目的本发明采用的的技术方案为:一种电网频率动态过程仿真轨迹的校正方法,包括以下步骤:
步骤1)通过故障录波器获取多个电网功率脱落事件的频率实测轨迹;
步骤2)在电力系统仿真软件中建立仿真模型和输入参数,获得电网动态过程频率仿真轨迹;
步骤3)根据步骤1)和步骤2)建立频率仿真轨迹和实测轨迹相近度的数学评价体系;所述数学评价体系概括为如下内容:
①定义频率动态过程仿真轨迹和实测轨迹接近度的特征指标,包括初始频降差值,频降斜率差值,回升频率差值,频率轨迹最大差值,频率轨迹平均差值;
②依据上述特征指标,建立频率仿真轨迹和实测轨迹误差评价函数;
③对多个电网功率脱落事件的实测轨迹和频率仿真轨迹进行上述误差计算,取最大值作为该仿真条件下误差值;
步骤4)确定要调整的调速器参数及其可调范围,形成参数变量集合;
步骤5)基于禁忌搜索算法,调整参数,通过如下步骤使频率仿真轨迹逼近频率实测轨迹;
第一步对调速器参数变量集合进行编码;
第二步产生初始解;
第三步将单个移动和交换移动分别作用于初始解,得到一组可行的试验解;
第四步对上述可行的试验解分别进行频率仿真计算和误差评价计算,从中保留误差最小的试验解;
第五步更新禁忌搜索表;
第六步更新初始解;
第七步判断终止条件:若连续几次迭代目标函数值没有改进或达到最大允许迭代次数,则停止计算,输出结果;
第八步若没有达到终止条件,重复第三步;
步骤6)输出结果,得到调整后的调速器仿真参数,用于频率控制或频率恢复策略。
所述获得电网动态过程频率仿真轨迹的过程如下:
a.在电力系统仿真软件中建立目标网架;
b.根据电网实际情况,选取和建立电气元件模型和录入相关参数;
c.在仿真软件中设置运行方式;
d.在仿真软件中设置故障;
e.运行电力系统仿真软件,得到频率仿真轨迹。
所述的数学评价体系为:
a.将各特征指标定义如下:
初始频降差值:为发生故障后频率初始下降到最低点,频率实测轨迹和仿真轨迹的频率差值;
为保证频率实测轨迹频降最低点高于仿真轨迹频率,定义差值函数为f:
频降斜率:频降斜率指频率开始下降到最低点的曲线斜率,这个值由初始频降与初始频降时间计算得到;
频降斜率差值:为仿真轨迹频降斜率和实测轨迹频率斜率差值的绝对值;
回升频率:回升频率指频率下降到最低点后频率回升所达到的最大值;
回升频率差值:为频率仿真轨迹和实测轨迹回升频率差值的绝对值;
频率轨迹最大差值:为频率仿真轨迹和实测轨迹在时间轴上的最大差值的绝对值;
频率轨迹平均差值:为频率曲线所有时间点上频率仿真轨迹和实测轨迹频率差值的均方差;
上述特征量通过存储在计算机中的频率曲线数据,进行计算得到;
b.某一故障条件下的频率实测轨迹和仿真轨迹的误差评价函数定义如下:
fe=∑(λ1T1+λ2T2+…λnTn) 式一
其中:fe为误差评价函数;T代表某个特征指标;λ表征相应特征指标在误差评价函数中的权重系数;n为所考虑的特征指标的个数;
c.在具有多个实际电力系统功率脱落事件下的实测频率轨迹的条件下,需要分别对多个功率脱落事件的实测轨迹和仿真轨迹进行上述误差计算,取数值最大者作为仿真轨迹和实测轨迹的误差值;
该误差值具体计算公式如下:
Fe=max(fe1,fe2,...fem) 式二
其中:Fe为在一组确定的调速器参数条件下,仿真轨迹和实测轨迹的误差;m为考虑的系统功率脱落事件下的实测频率轨迹个数。
所述的调速器仿真参数的调整,是指在电网模型和参数准确的前提下,对调速器相关参数进行小范围调整,改变频率仿真曲线的特性和形状,得到与实际电力系统情况相符的仿真轨迹。
所述的对参数变量集合进行编码,包括:
a.对某一可调参数,根据参数的可调节范围,把该参数间隔的进行离散化,形成该参数的待选集和;
b.对该参数待选集和中的每个元素,进行二进制编码;若某参数T的二进制编码对应的十进制数为d,则有:
式三
其中:Tmax和Tmin为该参数的最大和最小取值;n为二进制编码位数;
c.对所有参数进行上述编码工作;
d.把所有参数集合的二进制编码串联,形成一个(0-1)码串。
所述的产生初始解,包括:
a.对可调参数变量集合,赋予初值;
b.把该参数集合在电力系统动态仿真软件中进行修正;
c.针对不同的运行方式和故障条件,进行计算,得到不同功率脱落事件下的频率仿真轨迹;
d.计算在该参数值下,频率仿真轨迹和实测轨迹的误差函数值;
e.把参数的初值向量写成二进制编码。
所述的将单个移动和交换移动分别作用于初始解,包括:
a.单个移动:随机选取码串的某一位进行取反操作;
b.交换移动:随机选取码串的某两位进行取反操作;
c.将两种移动分别作用于当前解,从而得到一组调速器参数向量的试验解。
所述的可行试验解,包括:
a.该试验解不在禁忌搜索表中;
b.该试验解虽然在禁忌表中,但已满足释放准则。
对上述可行的试验解分别进行频率仿真计算和误差评价计算,从中保留误差最小的试验解;包括:
a.取某一个可行的试验解,通过其二进制编码,计算参数数值;
b.在电力系统仿真程序的调速器模型参数中,按上述数值进行修改;
c.在仿真程序中设置和功率脱落事件对应的运行方式和故障,进行仿真计算,得到相应的频率仿真轨迹;
d.计算误差评价函数;
e.与原有的评价函数值作比较,保留函数值最小者和相对应的参数值;
f.重复至a,直到所有试验解都计算完毕。
所述的更新禁忌搜索表,包括:
a.将实现了的移动的反方向移动记录到禁忌表中;
b.禁忌表采用先进先出的管理方式。
本发明的优点为:
1.本发明在电网模型和参数准确的条件下,以事故下系统真实的频率轨迹作为参照依据,对调速器参数在其可行的范围内进行调整,使得仿真轨迹逼近真实的频率轨迹。这种方法符合电力系统运行的实际情况。
2.本发明建立了频率仿真轨迹和实测轨迹误差的数学评价体系。与简单的通过轨迹形状进行二者接近度进行判断相比,这种评价方法更为科学和有效。
3.本发明采用禁忌搜索的方法来对参数进行调整,可以有效的实现频率仿真轨迹逼近实测轨迹。
4.本发明能够使电网频率仿真轨迹真实的体现系统的频率动态特性,为电力系统的频率控制和低频减载方案提供依据,提高电力系统运行的稳定性和可靠性。
具体实施方式
如图1、2、3、4、5-1、5-2所示,一种电网频率动态过程仿真轨迹的校正方法,包括以下步骤:
步骤1)通过故障录波器(市购产品)获取多个电网功率脱落事件的频率实测轨迹;
步骤2)在电力系统仿真软件(即电力系统分析综合程序(Power SystemAnalysis Software Package)简称PSASP)中建立仿真模型和输入参数,获得电网动态过程频率仿真轨迹;
步骤3)根据步骤1)和步骤2)建立频率仿真轨迹和实测轨迹相近度的数学评价体系;所述数学评价体系概括为如下内容:
①定义频率动态过程仿真轨迹和实测轨迹接近度的特征指标,包括初始频降差值,频降斜率差值,回升频率差值,频率轨迹最大差值,频率轨迹平均差值;
②依据上述特征指标,建立频率仿真轨迹和实测轨迹误差评价函数;
③对多个电网功率脱落事件的实测轨迹和频率仿真轨迹进行上述误差计算,取最大值作为该仿真条件下误差值;
步骤4)确定要调整的调速器参数及其可调范围,形成参数变量集合;
步骤5)基于禁忌搜索算法,调整参数,通过如下步骤使频率仿真轨迹逼近频率实测轨迹;
第一步对调速器参数变量集合进行编码;
第二步产生初始解;
第三步将单个移动和交换移动分别作用于初始解,得到一组可行的试验解;
第四步对上述可行的试验解分别进行频率仿真计算和误差评价计算,从中保留误差最小的试验解;
第五步更新禁忌搜索表;
第六步更新初始解;
第七步判断终止条件:若连续几次迭代目标函数值没有改进或达到最大允许迭代次数,则停止计算,输出结果;
第八步若没有达到终止条件,重复第三步;
步骤6)输出结果,得到调整后的调速器仿真参数,用于频率控制或频率恢复策略。
所述获得电网动态过程频率仿真轨迹的过程如下:
a.在电力系统仿真软件中建立目标网架;
b.根据电网实际情况,选取和建立电气元件模型和录入相关参数;
c.在仿真软件中设置运行方式;
d.在仿真软件中设置故障;
e.运行电力系统仿真软件,得到频率仿真轨迹。
所述的数学评价体系为:
a.将各特征指标定义如下:
初始频降差值:为发生故障后频率初始下降到最低点,频率实测轨迹和仿真轨迹的频率差值;
为保证频率实测轨迹频降最低点高于仿真轨迹频率,定义差值函数为f:
频降斜率:频降斜率指频率开始下降到最低点的曲线斜率,这个值由初始频降与初始频降时间计算得到;
频降斜率差值:为仿真轨迹频降斜率和实测轨迹频率斜率差值的绝对值;
回升频率:回升频率指频率下降到最低点后频率回升所达到的最大值;
回升频率差值:为频率仿真轨迹和实测轨迹回升频率差值的绝对值;
频率轨迹最大差值:为频率仿真轨迹和实测轨迹在时间轴上的最大差值的绝对值;
频率轨迹平均差值:为频率曲线所有时间点上频率仿真轨迹和实测轨迹频率差值的均方差;
上述特征量通过存储在计算机中的频率曲线数据,进行计算得到;
b.某一故障条件下的频率实测轨迹和仿真轨迹的误差评价函数定义如下:
fe=∑(λ1T1+λ2T2+…λnTn) 式一
式一中:fe为误差评价函数;T代表某个特征指标;λ表征相应特征指标在误差评价函数中的权重系数;n为所考虑的特征指标的个数;
c.在具有多个实际电力系统功率脱落事件下的实测频率轨迹的条件下,需要分别对多个功率脱落事件的实测轨迹和仿真轨迹进行上述误差计算,取数值最大者作为仿真轨迹和实测轨迹的误差值;
该误差值具体计算公式如下:
Fe=max(fe1,fe2,...fem) 式二
式二中:Fe为在一组确定的调速器参数条件下,仿真轨迹和实测轨迹的误差;m为考虑的系统功率脱落事件下的实测频率轨迹个数。
所述的调速器仿真参数的调整,是指在电网模型和参数准确的前提下,对调速器相关参数进行小范围调整,改变频率仿真曲线的特性和形状,得到与实际电力系统情况相符的仿真轨迹。
所述的对参数变量集合进行编码,包括:
a.对某一可调参数,根据参数的可调节范围,把该参数间隔的进行离散化,形成该参数的待选集和;
b.对该参数待选集和中的每个元素,进行二进制编码;若某参数T的二进制编码对应的十进制数为d,则有:
式三
式三中:Tmax和Tmin为该参数的最大和最小取值;n为二进制编码位数;
c.对所有参数进行上述编码工作;
d.把所有参数集合的二进制编码串联,形成一个(0-1)码串。
所述的产生初始解,包括:
a.对可调参数变量集合,赋予初值;
b.把该参数集合在电力系统动态仿真软件中进行修正;
c.针对不同的运行方式和故障条件,进行计算,得到不同功率脱落事件下的频率仿真轨迹;
d.计算在该参数值下,频率仿真轨迹和实测轨迹的误差函数值;
e.把参数的初值向量写成二进制编码。
所述的将单个移动和交换移动分别作用于初始解,包括:
a.单个移动:随机选取码串的某一位进行取反操作;
b.交换移动:随机选取码串的某两位进行取反操作;
c.将两种移动分别作用于当前解,从而得到一组调速器参数向量的试验解。
8、根据权利要求1所述的一种电网频率动态过程仿真轨迹的校正方法,其特征在于:所述的可行试验解,包括:
a.该试验解不在禁忌搜索表中;
b.该试验解虽然在禁忌表中,但已满足释放准则。
对上述可行的试验解分别进行频率仿真计算和误差评价计算,从中保留误差最小的试验解;包括:
a.取某一个可行的试验解,通过其二进制编码,计算参数数值;
b.在电力系统仿真程序的调速器模型参数中,按上述数值进行修改;
c.在仿真程序中设置和功率脱落事件对应的运行方式和故障,进行仿真计算,得到相应的频率仿真轨迹;
d.计算误差评价函数;
e.与原有的评价函数值作比较,保留函数值最小者和相对应的参数值;
f.重复至a,直到所有试验解都计算完毕。
所述的更新禁忌搜索表,包括:
a.将实现了的移动的反方向移动记录到禁忌表中;
b.禁忌表采用先进先出的管理方式。
本发明采用禁忌搜索技术来对参数进行适当调整,达到对频率仿真轨迹校正的目的。禁忌搜索又称为TS法,是一种扩展邻域的启发式搜索方法,最早由Fred Glover于1977年提出,作为一个最优工具来求解非线性覆盖问题。它利用灵活记忆的特殊形式能在搜索过程中获得知识,来避免搜索陷入局部最优。已成功地应用于求解大规模的组合最优问题。TS方法的基本原理为:从初始解X(为n维向量)开始搜索。通过采用一组操作符“移动”,对当前解进行操作,随机产生一组当前解的邻域N(X)试验解X1,X2,...,XK。在产生的这些邻域试验解中,选择其中最能改善评价函数的那个解X*作为当前局部最优解,即令X=X*,重复迭代,直到所有的‘移动’均不能改善评价函数,表明当前的解为局部最优解。为了避免陷入局部最优解,TS法中采用了一种灵活的“记忆技术”。另外,为了尽可能不错过产生最优解的“移动”,TS方法还采用“释放准则”策略。
采用禁忌搜索进行计算的一般步骤如下:
(1)读入原始变量集合进行编码。
(2)产生初始解:在控制变量约束范围内随机产生一个初始解X,求得目标函数值f(X),置最好解向量Xopt=f(X)
(3)产生一组试验解:将单个移动和交换移动分别作用于X,得到一组可行的试验解X1,X2,...,XK,并求得相应的f(X1),f(X2),...,f(XK)。
(4)搜索邻域:从上述试验解中寻优,得到一个最好的试验解X*,如果X*不在禁忌表中,或者虽在禁忌表中但已满足释放准则,则用X*更新X,即令X=X*;如果该解在禁忌表中,同时又没有满足释放准则,则寻找次好的解,并重复此过程。
(5)更新禁忌表:将实现了的移动的反方向移动记录到禁忌表中,如果禁忌表已满,则首先排除最先记录的移动。
(6)更新Xopt:若f(X*)<f(Xopt),则用X*更新Xopt,否则Xopt保持不变。
(7)判断终止条件:若连续几次迭代目标函数值没有改进或达到最大允许迭代次数,则停止优化,输出结果;否则,转到步骤(3)继续迭代。
与其它搜索算法相比,禁忌搜索法具有较强的跳出局部最优解的能力,因而收敛特性好、解的质量高,已经被成功应用于大规模的电力系统参数优化或调整的问题中。本发明采用禁忌搜索技术对影响频率动态特性的相关参数进行适当调整,来达到对频率仿真轨迹校正的目的。由于本方法的目的是实现对频率仿真轨迹的有效调整,使之逼近实测轨迹,更真实地体现电力系统的频率动态特性,因此搜索过程中,不必一定到达最优点,一定的迭代次数后,也能有效的改进频率仿真轨迹。
本发明在电网模型和参数基本准确条件下,以事故下系统真实的频率轨迹作为参照依据,对调速器参数在其可行的范围内进行微调,使得仿真轨迹逼近真实的频率轨迹;建立误差评价模型,通过轨迹误差评价函数,来评价频率轨迹的逼近程度;通过禁忌搜索,来调整相关参数,从而实现对频率仿真轨迹的有效校正。本发明采用的方法,可以避免电力系统中频率动态过程研究和控制无据可依,大量凭借经验的现状,为电网频率轨迹的校正提供了有效而切实的途径。本发明采用的方法通过对频率仿真轨迹进行有效校正,使电网频率仿真轨迹真实的体现系统的频率动态特性,为电力系统的频率控制和低频减载方案提供依据,提高电力系统运行的稳定性和可靠性。
图1是通过参数调整使频率仿真轨迹接近实测轨迹的总体流程图。其中,实际系统指的是要研究的某个电网,对应着实际的物理网络;仿真系统,是指在电力系统仿真程序中建立的与实际系统对应的数字模拟系统。无论是频率实测曲线,还是频率仿真曲线,都是在系统中产生功率脱落事件的条件下得到的。
①、频率量测。通过现代电网的能量管理系统和装设在变电站内的故障录波器,在系统发生功率脱落事件时,可以得到频率实测轨迹的数据。每次发生功率脱落事件的情况都不尽相同,可以选择对当前电网最有参考价值的若干个事件,作为后期频率仿真轨迹调整的依据。在记录频率实测轨迹的同时,也要通过能量管理系统,记录发生功率脱落事件时,所对应的系统运行方式、控制方式等,以便在仿真计算时,与其相对应。
②、功率电压稳态量量测。通过能量管理系统的数据采集与监控系统,获取发生功率脱落时刻,实际电网的运行状态
③、潮流调整、扰动设置。在仿真系统中对系统稳态运行方式进行调整,使其与发生功率脱落事件时刻系统的运行状态相对应;根据功率脱落情况,在仿真软件中进行扰动设置,例如切机、切负荷和故障等;对于选定的多个功率脱落事件,应分别进行相应的潮流调整和扰动设置,逐个进行计算,进而得到多个功率脱落事件下的频率仿真轨迹。
④、轨迹接近度评价计算。用于确定仿真轨迹和实测轨迹的接近程度。首先定义了一些反映轨迹接近程度的特征指标,这些指标可以通过对轨迹曲线的简单计算而得到;通过综合考虑这些特征指标,来衡量频率仿真轨迹和实测轨迹的接近程度(具体定义和计算方式见说明部分);在计算时,各特征量的权重系数的取值可以根据调整的要求而定。一般而言,前三个特征量:初始频降差值、频降斜率差值和回升频率差值决定了轨迹的主要形状,应赋予较大的比重;其他的特征量所占比重略小。
⑤、禁忌搜索、调整参数。如果频率仿真轨迹和实测轨迹吻合度不好,误差较大,则需要对相关参数进行调整,重新仿真计算,再次校核。这里在电网模型和参数基本准确条件下,对调速器参数在其可行的范围内进行调整,来改进仿真轨迹。根据确定要调整的调速器参数及其可调范围,形成参数变量集合,作为禁忌搜索的变量集合;然后依据禁忌搜索的计算规律,不断进行仿真计算,进而得到调整后的调速参数。更具体的步骤将在下面结合图3做进一部的说明。
⑥、确定参数。通过搜索计算,最后可以得到一组比较满意的调速参数;在这组参数下进行仿真计算,得到的频率仿真轨迹与实测频率轨迹具有较好的吻合度。基于仿真系统调整后的系统模型和参数,就可以在此基础上来研究系统的一些预想的频率事件和频率动态特性,为电力系统的频率控制奠定基础。正是由于对仿真系统中的相关参数进行了有效校正,提高了仿真轨迹的准确度,才使基于频率仿真的研究工作更加可信和可靠。
如图2所示,仿真轨迹和实测轨迹特征指标示意图。
其中的实线曲线是某大区电网实测的频率轨迹,图中的虚线曲线是实测轨迹相对应的运行状态和扰动下,频率仿真曲线。图中的黑色圆圈点是轨迹的初始频降点;图中的直线表征轨迹的频降斜率;图中的白色圆圈点是轨迹的回升频率。通过这些特征指标,可以有效的表述轨迹的接近程度。
如图4所示,程序模块和数据调用关系图。
为实现本发明所述校正方法,采用电网动态仿真程序软件,误差评价计算和禁忌搜索;禁忌搜索程序对要调整的参数变量集合参数进行优化调整,在调整过程中,对某一组既定的参数进行评价的时候,就需要调用轨迹误差评价计算程序,来决定该组参数的取舍;而轨迹误差评价计算程序在计算的时候,需要知道测试频率轨迹和仿真频率轨迹的数据。其中频率实测轨迹数据可以事先得到,保存在存储介质上;仿真频率轨迹则需要根据相应的参数值,进行仿真计算得到。因此轨迹误差评价计算程序需要调用电力系统动态仿真程序,得到频率仿真轨迹数据。当前,电力系统仿真计算软件种类较多,在我国电力系统广泛应用的电力系统仿真计算软件主要有中国电力科学研究院开发的电力系统分系综合程序PSASP,邦纳维尔电力局开发的电力系统离线分析程序BPA(Bonneville Pow er A dm in ist rat ion)等。一般这类程序都提供了数据和调用的接口,因而上述整个过程可以通过程序来实现。
如图5-1、5-2所示,国内某区域电网实施参数调整前后频率轨迹对比图。
本发明所述方法已经成功应用于国内某个区域电网的频率仿真轨迹校正工作中。该电网总装机容量为6300万千瓦,500KV线路76条,220KV线路768条,500KV变电站30座,220KV变电站358座。对这个系统近一年半时间,收集了若干个功率脱落事件和相应的频率实测轨迹。通过频率实测轨迹和仿真轨迹的比对发现,二者之间存在较大的误差。这种频率仿真轨迹无法用于系统频率动态过程的研究和频率控制。采用本发明提出的方法,进行了计算分析。禁忌搜索部分和轨迹误差评价部分自行编制程序,电力系统仿真软件采用PSASP,整个计算在一台PC机上完成。禁忌搜索经过30次迭代后的计算结果(如图5-1、5-2所示)。
其中黑色虚线为实测频率轨迹曲线,黑色实线为经过调整计算后的频率仿真曲线。可见,本方法可以对频率仿真曲线进行有效的调整,更加逼近真实的频率轨迹。