动态模型参数估计的初始化

申请号 CN200680040903.6 申请日 2006-10-31 公开(公告)号 CN101300532B 公开(公告)日 2010-12-15
申请人 ABB研究有限公司; 发明人 P·科尔巴; M·齐马; M·拉森;
摘要 本 发明 涉及电 力 系统中机电振荡的监控,以及如EP-A 1 489 714中详细描述的根据重复测量并评估的 信号 而由适应性 算法 进行的它们的识别。为了估计合理收敛的电力系统的模型参数,要求递归计算的适当初始化,包括限 制模 型和计算的调谐参数的定义。然后通过复制之前为不同信号调谐的这组调谐参数可将所用的第二信号的初始化简化。调节增益乘以第二信号建立了不同信号之间的兼容性,信号的预先滤波又去掉了包括典型电磁振荡的频带之外的影响。
权利要求

1.一种初始化由电系统(1)的参数模型的估计的模型参数(a1,a2,...)推导电力系统(1)的机电振荡模式的频率和阻尼(f,ξ)的方法,其中模型参数(a1,a2,...)的估计基于电力系统(1)的第二系统量(y2)的一系列测量的值(y21,y22)并且其中每次测量第二系统量(y2)的新值(y2k)时适应性地估计所述模型参数(a1,a2,...),其中该初始化方法包括:
-调谐一组调谐参数(tp2),用于模型参数(a1,a2,...)的随后的估计,
其特征在于该方法还包括:
-通过根据电力系统(1)的第一系统量(y1)复制先前为估计模型参数(a1,a2,...)而调谐的调谐参数(tp1)来调谐这组调谐参数(tp2),以及
-在模型参数(a1,a2,...)的每次适应性估计之前确定对第二系统量(y2)的测量的值(y21,y22,...)进行缩放的调节增益(G2)。
2.根据权利要求1的方法,其特征在于确定调节增益(G2)包括比较有关第一和第二系统量(y1,y2)的统计度量(Sr1,Sr2)。
3.根据权利要求1的方法,其特征在于其包括提供带通滤波器(f1,f2),用于在借助调节增益(G2)进行缩放之前对电力系统(1)的第二系统量(y2)的该系列测量的值(y21,y22,...)进行滤波。
4.根据权利要求3的方法,其特征在于其包括通过带通滤波器(f1,f2)截止不代表机电电力系统振荡的频率。
5.根据权利要求3的方法,其特征在于其包括提供阈值,用于如果最后滤波和缩放的值(y2k)低于所述阈值那么就忽略它。
6.一种用于由电力系统(1)的参数模型的估计的模型参数(a1,a2,...)推导电力系统(1)的机电振荡模式的频率和阻尼(f,ξ)的系统,包括用于测量第一和第二系统量(y1,y2)的两个测量单元(20)以及用于根据电力系统(1)的第二系统量(y2)的一系列测量的值(y21,y22,...)估计模型参数(a1,a2,...)的监控中心(21),其中每次测量第二系统量(y2)的新值(y2k)时适应性地估计所述模型参数(a1,a2,...),并且其中调谐一组调谐参数(tp2)用以初始化模型参数(a1,a2,...)的随后的估计,其特征在于,该系统包括
-用于调谐该组调谐参数(tp2)的装置,其中该调谐是通过根据电力系统(1)的第一系统量(y1)复制之前为估计模型参数(a1,a2,...)调谐的调谐参数(tp1)而进行的,以及
-用于确定调节增益(G2)的装置,该调节增益(G2)用于在模型参数(a1,a2,...)的每次适应性估计之前对第二系统量(y2)的测量的值(y21,y22,...)进行缩放。
7.根据权利要求1到4中的一个权利要求的方法的使用,用于根据由卡尔曼滤波技术估计的模型参数(a1,a2,...)推导电力系统(1)中的机电振荡的频率和阻尼(f,ξ)。
8.根据权利要求7的使用,其特征在于对缩放参数G进行在线调节。
9.一种用于使用计算机程序控制电力系统(2)中的功率流以及阻尼电磁振荡的方法,其中可将计算机程序加载到数字计算机的内部存储器中,并且该计算机程序包括计算机程序代码装置,当所述程序加载到所述内部存储器中时,计算机程序代码装置使得计算机根据权利要求8执行控制器(1)的功能。

说明书全文

技术领域

发明涉及监控包括多个发电机和用户的电系统中的电磁振荡的领域。其从权利要求1的前序部分中所述的初始化电力系统的参数模型的模型参数估计的方法中出发。

背景技术

随着正在进行的电力市场的反调节(deregulation),从遥远的发电机到本地用户的负载输送以及电力运转已经变成惯例。作为优化资产的新兴需要和公用事业(utility)之间竞争的结果,增长的电功率量通过现有的网络进行输送,这样总是造成部分输电系统的过剩,输送瓶颈和/或振荡。在这个方面上,输电网络非常动态。总的来说,包括几个交流发电机的电力系统中的电磁振荡频率小于几赫兹,并且只要它们可以衰减,就认为可以接受。它们由系统负载中正常的小变化发起,并且它们是任何电力系统的特性。然而,在电力系统的工作点改变时,例如,由于随着发电机,负载和/或输电线的连接或断开所带来的功率流的新分配而引起的改变,将发生不充分的阻尼振荡。同样地,即使它们互连之前后者不单独出现任何严重阻尼振荡,几种现有电网的互连将发生不足够的阻尼振荡。在这些情况下,几兆瓦的输电功率的增加在稳定振荡和不稳定振荡之间形成差别,其中不稳定振荡的电位使得系统崩溃或者导致失同步,失互连,并且最终不能提供电功率给用户。对电力系统进行适当的监控可有助于网络工作人员精确地评定电力系统状态,通过像连接特定设计的振荡阻尼设备那样采取适当的动作来避免总停电。
专利申请EP-A 1 489 714中,电力系统中机电振荡的适应性检测基于线性时变模型。采样像例如网络选择节点上的电压电流的幅值或度之类的系统量或信号,并且通过卡尔曼滤波技术估计表示电力系统性能(behavior)的线性模型参数。这个过程以递归方式进行,也就是说,每次测量系统量的新值时,更新模型参数。最后,从模型的估计参数中推导出例如频率或阻尼之类的振荡模型的参数,并将这些参数呈现给工作人员。这种适应性识别过程可以进行电力系统当前状态的实时分析,尤其包括主要的功率振荡模式的阻尼ξ以及频率f,该主要的功率振荡模式就是具有最低相对阻尼率的模式。
为了使得这种动态模型参数的估计正确运行,必须由一组适当选择的调谐参数对该估计进行初始化,这些调谐参数例如动态模型的阶,过程和测量噪声,信号预先滤波的截止频率等等。总的来说,调谐参数的值取决于监控的特定电力系统以及选择成监视算法输入的特定信号。然后通过调试工程师进行调节或调谐,其中调试工程师分析每个输入信号并确定下一步的估计过程的输出,即,被估计的主频率以及阻尼,足够快速地响应并且相对于测量噪声不是太敏感。特别地,调试工程师不得不这样调节调谐参数的值:由被测信号和例如由上述线性时变模型预期的信号之间的差所给定的估计误差最小,并且采用尽可能小阶的动态模型来足以精确估计感兴趣的捕获振荡模式。结果是这种调谐程序时间密集,并且要求调试工程师具有一定程度的知识和经验。
为了识别电力系统中的振荡,诸如电压,电流和功率流的幅值或相角之类的不同系统量可用作所提出的识别程序的输入。然而,这种信号相对于它们的诸如幅值和信号方差的统计性质不同。为了简化调谐程序,即,为了找到要求启动估计算法的调谐参数的最佳初始值,上述欧洲专利申请提出引入信号调整用于从监控的电力系统得到的所有可容许测量值。

发明内容

本发明的目的是增加检测和监控电力系统中的机电功率振荡的灵活性,而不增加调试的工程复杂性或者工作量。这些目的可通过根据权利要求1和5的初始化电力系统参数模型的模型参数估计的方法以及监控电力系统的系统来实现。进一步的优选实施例从专利的从属权利要求中可以明显得到。
根据本发明,从以下事实得到优势:对于同一个电力系统而言,多个不同的系统量可用,即,不同位置上测量的不同信号。然而,在这些不同的输入信号中,可观察到相同的动态现象,例如机电振荡现象。因此,例如,为了地理或者时间上跟踪某个振荡模式,可从一个信号交换到另一个信号上,其中某个振荡模式例如中间区域模式可从第一位置的测量电压并且可从电力系统的另一位置的电流信号中可观察到,或者例如可随着输电线或发电机的连接或断开而切换。
为了避免在采用例如最初所述的检测机电振荡的方法时对每个系统量单独进行调谐,仅单独确定第一或者参考系统量的调谐参数。然后可通过复制或者重新使用全部或部分上述的调谐参数并通过确定足够的重新缩放因子作为信号调节增益而简化第二或者其他的系统量的初始化程序,其中第二或其他的系统量可例如为某个振荡模式提供更好的可观察性。后者可通过比较第一和第二系统量来确定,其能承受兼容的不同输入信号并且是源自第一系统量的调谐参数值能成功再使用的前提。然后应用包括所述调节增益的调谐参数值复制组,以根据第二系统量的一系列测量值来识别代表电力系统性能的模型参数。
优选的是,通过比较包含在被测信号中的统计信息来确定调节增益,这些被测信号例如第一和第二系统量的多个测量值的最大信号功率,平均值或者均方根值。可以实时设置缩放因子的适应性。
在有利的实施例中,在上述信号调节之前为第二系统量的测量值提供带通滤波器。该滤波器可基于被跟踪的振荡的一般知识,或者基于(即围绕)第一或参考系统量从以前的分析中得到的主要机电振荡频率进行限定。
上述简化的参数估计调谐过程的使用涉及诸如电力系统中主要振荡模式的频率或阻尼之类的信息从被估计的动态模型参数中的推导。为此,动态模型参数优选地由卡尔曼滤波技术确定。
附图说明
下文中参考附图示出的优选示意性实施例详细解释了本发明的主题,其中:
图1示意性地示出电力系统,
图2示出估计模型参数的过程的流程图
图3是根据第一系统量y1分析电力系统的结果,
图4是根据第二系统量y2并且采用相同的调谐参数初始化分析相同的电力系统的结果。
附图中所用的附图标记以及它们的含义在附图标记列表中以概括的形式列出。原则上,相同的部件在附图中采用相同的参考标记。

具体实施方式

图1示出电力系统1,包括两个发电机10,10’和分别由汇流排(busbar)11,11’,11”表示并由多个传输线互连的几个变电站。诸如电压或电流的相角和/或幅值,频率,功率流等之类的系统量y1,y2,y3由位于电力系统1上不同变电站或节点上的合适测量单元20,20’,20”进行测量。将测量单元20,20’,20”测量的信号发送到振荡监控中心21中并在其中进行使用。总的来说,几个测量单元20,20’,20”可在一个单独的设备中实现,此外其不需要成为专用设备,相应的测量功能同样可由为系统1中的保护和控制任务所设的智能电子设备来执行。此外,监控中心21可以与其中一个测量单元20相同。
如上所述,模型参数的适应性估计的适当初始化要求递归计算所用的调谐参数的调谐或者离线调节。举例来说,在上述欧洲专利申请EP-A 1 489 714中所提出的程序中,选择离散时间自回归模型的动态阶n是最重要的一个方面,其中这个阶等于将要估计的参数数量。如果这个阶太低,在频域中得到的频谱将高度平滑,并且在频谱中具有低值尖峰的感兴趣的振荡被掩盖。另一方面,如果阶n太高,伪低值尖峰将被引入到频谱之中。此外,测量噪声Qm和过程噪声Qp的相关矩阵表示其他较不敏感的调谐参数。其他调谐参数是系统量y的连续测量值之间的采样时间Ts,信号预先滤波的截止频率f1,f2以及信号调节因子或增益G,如下面将要详细描述的。
图2示出了上述欧洲专利申请EP-A 1 489 714中所述的电力系统振荡监控的适应性实时算法的有利改进,其中该欧洲专利申请公开的内容在此引入以供参考。在初始化步骤30中,根据本发明确定与系统量y2一起使用的调谐参数tp2,即,为不同的系统量y1事先确定的调谐参数tp1进行复制。步骤30包括例如基于考虑中的系统量y1,y2的有限数量的测量值{y1’},{y2’}的离线分析最初确定调节增益G2,并包括下面所述的滤波和统计步骤。在该算法的重复执行过程中,在测量步骤31中采用1/Ts的采样或更新频率测量第二系统量y2的新值y2(k)。然后在滤波步骤32中对这一系列测量值y2(k)进行带通滤波,其中作为调谐参数的截止频率f1,f2已经在上面介绍,从而产生y2(k)的一系列滤波值。这一系列滤波值的统计度量在统计步骤33中针对调节增益G2的最后更新而确定。最后,在缩放步骤34中,y2(k)的该系列滤波值采用调节增益G2的实际值重新缩放。如果最后测量的,滤波的并且缩放的值y2(k)超过了特定阈值,和/或如果某个计数器显示出这样,在模型参数更新步骤35中进一步采用y2(k)的该系列滤波的缩放的值,如本领域中已知的那样。
更详细地,重新缩放之前的带通滤波步骤32去掉了例如0.1赫兹的较低截止频率f1以下的DC分量以及例如2赫兹的较高截止频率f2之上的高频。电力系统振荡的一般频率是已知的事实,允许限定如所指示的带通范围,然而,如果例如递归算法的结果显示这样做,那么任何时刻可以调节截止频率f1,f2。
在被认为初始化或更新调节增益G2的周期T中测量的值的带通滤波序列s(k)的统计度量包括例如最大信号功率,平均值,或者均方根值,如下:
Sr=max0<k<Ts(k)(最大值)
Sr=1TΣk=1T|s(k)|(平均值)
Sr=1TΣk=1Ts(k)2(均方根值)
调节增益G2可根据考虑的第一和第二系统量y1,y2的相应统计度量Sr1,Sr2通过除法:G2=Sr1/Sr2而计算得到。
在低信噪比的情形下,例如在测量单元故障的情况下,输入信号(即系统量的一些随后的测量值)可临时构成噪声,其平均值接近零而不是实际数据。这样考虑所有测量值确切等于零是有益的,否则估计的是噪声的主频率而不是被测信号的被测信号的主频率。基于平均信号功率的观察,固定阈值,并且如果实际信号功率低于该阈值,将估计的模型参数冻结(不更新)。
下面,一个实例示出所提出的程序的效果,其中选择两个完全不同的信号并且采用所开发的用于检测振荡的工具进行分析。实际测量的数据包括在Ts=0.05秒的时间间隔采样的1600个值y1k,y2k的两个序列,对应于80秒的短数据采集时间间隔。
第一系统量y1:输入信号是RMS幅值为400kV±2kV或者以传统符号表示为1p.u.±0.005的交流电压,其中1p.u.=400kV。这在图3中的第一曲线图上示出。第二个曲线图上,示出经过滤波的信号。采用特定组的总共19个调谐参数tp1,随后的适应性程序导致主频率f以及其相对阻尼ξ的估计,如图3的第三曲线图和第四曲线图示出,在所示的时间间隔内收敛到值f≈0.45Hz,ξ≈17%。当事先没有包括其他的信息并且估计的全部模型参数(a1,a2,...)从任何初始值(这里为零)开始并且快速收敛到它们的正确值时,两个底部曲线图中的初始尖峰由模型参数估计算法的瞬态性能所导致。
第二系统量y2:输入信号是功率线路中的功率流,其值为图4中的第一曲线图所示的1350MW±60。这种信息在收集几个采样值并运行第一分析之后立即对调试工程师可用。根据本发明,基于功率流测量的第二系统量的调谐参数tp2根据电压测量从第一组进行复制。这种情况中经过滤波的第二信号所用的调节增益G2可计算成G2=0.005/60=8.3e-5。结果,采用第二系统量y2,所估计的振荡参数频率和相对阻尼,显然以与基于第一系统量y1(图3)根据递归计算的参数类似的速度收敛(图4,第三和第四曲线图)。因此,这种简化的初始化实际上使得工作时间以及调谐过程最小化,而不负面地影响结果的质量
附图标记列表:
1     电力系统
10    发电机
11    汇流排
20    测量单元
21    振荡监控中心
QQ群二维码
意见反馈