一种在线辨识电力系统低频振荡特征参数的方法与系统技术方案

技术编号:27066120 阅读:41 留言:0更新日期:2021-01-15 14:47
本发明专利技术公开了一种在线辨识电力系统低频振荡特征参数的方法与系统。本发明专利技术的方法包括步骤:输入卡尔曼滤波器的初始值;根据k时刻状态空间变量值计算下一时刻状态空间变量预测值;更新k+1时刻的预测误差协方差矩阵;计算k+1时刻的卡尔曼增益;计算k+1时刻的误差协方差矩阵;计算k+1时刻的状态空间变量值;输出功率振荡特征参数辨识结果。本发明专利技术针对不同机组功率信号采样频率不同造成的收敛速度慢的问题,提出变速率扩展卡尔曼滤波方法,能够有效提高辨识收敛速度,降低辨识误差。

【技术实现步骤摘要】
一种在线辨识电力系统低频振荡特征参数的方法与系统
本专利技术属于电力系统运行与控制领域,涉及一种在线辨识电力系统低频振荡特征参数的方法与系统。
技术介绍
电力系统低频振荡是指电力系统在受到干扰的情况下发生的一种功角稳定性问题,通常表现为有功功率在0.1-2.5Hz的频率范围内振荡,如振荡幅值得不到收敛,它将威胁电力系统安全,严重时会导致电力系统的崩溃。电力系统低频振荡对汽轮发电机组的安全稳定运行有较大的影响,当电力系统功率晃动幅值超出一定范围或不能自行衰减,将会给机组的安全稳定运行带来严重威胁。随着发电机组跟随电网功率振荡,机组转速超过死区,一次调频回路发挥作用,汽轮机调节阀根据一次调频回路指令随着转速的振荡而反复开关晃动,此时应尽快采取措施平息转速和功率振荡,减小低频振荡对发电机组造成的危害。但是目前发电机组控制系统尚无法检测识别低频振荡,机组运行人员不能在发生低频振荡时及时采取措施,长时间振荡会应发跳机危险。电力系统低频振荡主要特征参数包括:振荡频率以及振荡阻尼比。由于发电机组低频振荡特征参数无法直接测量,因此加大了对发电机组功率振荡的识别难度。
技术实现思路
为解决上述现有技术存在的问题,及时检测识别汽轮发电机组强迫振荡,并针对性的采取应对措施,使发电机组快速恢复平稳运行,本专利技术提供一种基于变速率扩展卡尔曼滤波器的在线辨识电力系统低频振荡特征参数的方法与系统,其用于识别汽轮发电机组低频振荡并及时采取措施抑制汽轮发电机组长时间振荡,防止低频振荡事件恶化而危及电力系统安全。本专利技术采用的一种技术方案为:一种在线辨识电力系统低频振荡特征参数的方法,其包括步骤:输入卡尔曼滤波器的初始值;根据k时刻状态空间变量值xk计算下一时刻状态空间变量预测值x-k+1;更新k+1时刻的预测误差协方差矩阵P-k+1;计算k+1时刻的卡尔曼增益Kk+1;计算k+1时刻的误差协方差矩阵Pk+1;计算k+1时刻的状态空间变量值xk+1;输出功率振荡特征参数辨识结果。本专利技术采用扩展卡尔曼滤波方法对汽轮发电机组功率信号进行处理,当机组发生功率振荡时,通过卡尔曼滤波器处理后输出功率实时振荡频率以及阻尼比,并将该信号输出至控制器中;考虑到不同机组功率信号采样频率不同,本专利技术提出采用变速率扩展卡尔曼滤波的方法,解决功率信号采样频率低带来的误差。为解决不同发电机组功率信号采样频率不同导致滤波器输出存在误差的问题,本专利技术采用可变速率扩展卡尔曼滤波的方法,在采样点之间增加计算点数来加快拟合速度,降低误差。进一步地,所述卡尔曼滤波器的初始值包括状态空间变量初值x0、误差协方差矩阵初值P0、功率振荡频率初值ω0、功率振荡阻尼比初值σ0、过程噪声协方差初值Q0以及测量噪声协方差初值R0。进一步地,将所述的发电机组功率振荡视为多个带指数增益余弦信号之和,其数学模型用下式表示:式中:N表示组成功率振荡余弦信号个数;A表示余弦信号幅值;ω表示振荡频率;σ表示振荡阻尼比;φ表示信号初相位;n(t)表示零均值白噪声信号;对单个带指数增益余弦信号进行分析,对功率信号以频率fs进行采样,k采样时刻发电机组功率如下式所示:式中:A表示振荡信号幅值;k表示采样次数;fs表示采样频率;σ表示功率振荡阻尼比;φ表示功率信号初相位;ω表示振荡频率;为便于采用卡尔曼滤波器处理功率信号,定义状态空间变量值xk如下式所示:式中,x1、x2、x3、x4分别表示状态空间变量xk中的第1、2、3、4行参数;根据上述状态空间变量建立用于卡尔曼滤波器的状态空间方程以及测量方程;根据公式(3)搭建如下状态空间方程:式中:wk为一个四行一列矩阵,其表示过程噪声,该噪声为均值为零的白噪声;将公式(3)代入到公式(2)中可得到k采样时刻的测量方程如下:yk=h(xk)+nk,(5)其中:h(xk)=x1,k·cos(φ)-x2,k·sin(φ)(6)其中:nk表示测量噪声,该噪声为均值为零的白噪声;x1,k、x2,k分别表示在k时刻状态空间变量xk的第1、2行参数。进一步地,根据k时刻状态空间变量值xk计算k+1时刻的状态空间变量预测值x-k+1;如下式所示:其中,函数f(xk,wk)对应于公式(4)所示状态空间方程,表示k+1时刻状态空间变量的预测值。进一步地,更新k+1时刻预测误差协方差矩阵P-k+1的公式如下:其中:Fk、Lk表示状态空间方程的雅克比矩阵,计算公式分别为:根据式(4)知Lk为一个4阶单位矩阵,即Lk=I4×4;Qk为k时刻过程噪声协方差;Pk为k时刻误差协方差矩阵。进一步地,k+1时刻卡尔曼增益的计算公式如下所示:式中:Hk表示测量方程的雅克比矩阵,计算公式为函数h(xk)对应于公式(6)所示测量方程,Rk+1表示测量噪声协方差。进一步地,根据k+1时刻预测误差协方差矩阵P-k+1以及卡尔曼增益Kk+1更新k+1时刻的误差协方差矩阵,计算公式如下所示:式中:I为对应维度单位矩阵。进一步地,k+1时刻的状态空间变量值xk+1的计算公式如下所示:上式中:Yk+1表示发电机组功率信号在k+1时刻的实际测量值;对应于公式(6)所示方程,表示根据k+1时刻状态空间变量的预测值计算得到的功率信号预测值;在发电机组运行时,通过对实时测量的功率数据的滤波处理,得到该时刻发电机组功率振荡特征参数,即振荡频率以及阻尼比。本专利技术采用的另一种技术方案为:一种在线辨识电力系统低频振荡特征参数的系统,其包括:初始值输入单元,用于输入卡尔曼滤波器的初始值;状态空间变量预测值计算单元,用于根据k时刻状态空间变量值xk计算下一时刻状态空间变量预测值x-k+1;预测误差协方差矩阵更新单元,用于更新k+1时刻的预测误差协方差矩阵P-k+1;卡尔曼增益计算单元,用于计算k+1时刻的卡尔曼增益Kk+1;误差协方差矩阵计算单元,用于计算k+1时刻的误差协方差矩阵Pk+1;状态空间变量值计算单元,计算k+1时刻的状态空间变量值xk+1;辨识结果输出单元,用于输出功率振荡特征参数辨识结果。基于上述在线辨识电力系统低频振荡特征参数的方法,本专利技术提供一种在线辨识电力系统低频振荡特征参数的系统。将上述变速率扩展卡尔曼滤波算法布置在系统中,以功率信号作为输入测量信号,输出功率振荡特征参数振荡频率以及阻尼比。本专利技术具有的有益效果如下:本专利技术能实时在线辨识并输出发电机组功率振荡特征参数;当发电机组发生功率振荡时,可通过特征参数的变化向机组发出报警信号,解决了无法直接测量功率振荡特征参数的难题,为快速识别发电机组功率振荡提供了理论依据,可防止汽轮发电机组长时间功率振荡而危害发电机组以及电网安全。附图说明本文档来自技高网...

【技术保护点】
1.一种在线辨识电力系统低频振荡特征参数的方法,其特征在于,包括步骤:/n输入卡尔曼滤波器的初始值;/n根据k时刻状态空间变量值x

【技术特征摘要】
1.一种在线辨识电力系统低频振荡特征参数的方法,其特征在于,包括步骤:
输入卡尔曼滤波器的初始值;
根据k时刻状态空间变量值xk计算下一时刻状态空间变量预测值x-k+1;
更新k+1时刻的预测误差协方差矩阵P-k+1;
计算k+1时刻的卡尔曼增益Kk+1;
计算k+1时刻的误差协方差矩阵Pk+1;
计算k+1时刻的状态空间变量值xk+1;
输出功率振荡特征参数辨识结果。


2.根据权利要求1所述的一种在线辨识电力系统低频振荡特征参数的方法,其特征在于,所述卡尔曼滤波器的初始值包括状态空间变量初值x0、误差协方差矩阵初值P0、功率振荡频率初值ω0、功率振荡阻尼比初值σ0、过程噪声协方差初值Q0以及测量噪声协方差初值R0。


3.根据权利要求1或2所述的一种在线辨识电力系统低频振荡特征参数的方法,其特征在于,将所述的发电机组功率振荡视为多个带指数增益余弦信号之和,其数学模型用下式表示:



式中:N表示组成功率振荡余弦信号个数;A表示余弦信号幅值;ω表示振荡频率;σ表示振荡阻尼比;φ表示信号初相位;n(t)表示零均值白噪声信号;
对单个带指数增益余弦信号进行分析,对功率信号以频率fs进行采样,k采样时刻发电机组功率如下式所示:



式中:A表示振荡信号幅值;k表示采样次数;fs表示采样频率;σ表示功率振荡阻尼比;φ表示功率信号初相位;ω表示振荡频率;
为便于采用卡尔曼滤波器处理功率信号,定义状态空间变量值xk如下式所示:



式中,x1、x2、x3、x4表示状态空间变量xk中的第1、2、3、4行参数;
根据上述状态空间变量建立用于卡尔曼滤波器的状态空间方程以及测量方程;
根据公式(3)搭建如下状态空间方程:



式中:wk为一个四行一列矩阵,其表示过程噪声,该噪声为均值为零的白噪声;
将公式(3)代入到公式(2)中可得到k采样时刻的测量方程如下:
yk=h(xk)+nk,(5)
其中:
h(xk)=x1,k·cos(φ)-x2,k·sin(φ)(6)
其中:nk表示测量噪声,该噪声为均值为零的白噪声;x1,k、x2,k分别表示在k时刻状态空间变量xk的第1、2行参数。


4.根据权利要求3所述的一种在线辨识电力系统低频振荡特征参数的方法,其特征在于,根据k时刻状态空间变量值xk计算k+1时刻的状态空间变量预测值x-k+1;如下式所示:



其中,函数f(xk,wk)对应于公式(4)所示状态空间方程,表示k+1时刻状...

【专利技术属性】
技术研发人员:鲍文龙顾正皓张宝丁阳俊高宽吴瑞康应光耀
申请(专利权)人:国网浙江省电力有限公司电力科学研究院杭州意能电力技术有限公司
类型:发明
国别省市:浙江;33

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1