CN104808090A - 一种基于改进强跟踪滤波器的电力信号突变参数测量方法 - Google Patents

Abstract

本发明公开了一种基于改进实数型强跟踪滤波器(RSTF)的电力系统失真信号的参数估计方法，属于电力系统领域。该方法包括：步骤一：获取目标信号；步骤二：确定目标信号特征；步骤三：针对目标信号，建立滤波模型，此模型为连续三次采样正弦波电压的等式关系模型；步骤四：根据RSTF递推规则，得到RSTF的离散算法公式；步骤五：采用遗传算法对RSTF的参数设置进行在线优化求解，之后对离散目标信号进行滤波和参数估计。因为采用连续三次采样正弦波电压的等式关系模型，使得算法成为隐式双层滤波，提高了滤波能力。而RSTF的自适应和鲁棒性等，保证了此方法在线应用的实时性和有效性。另外，遗传算法优化了算法参数的求解，从而可以进一步提高算法的跟踪能力。

Description

一种基于改进强跟踪滤波器的电力信号突变参数测量方法

本发明属于电力系统领域，涉及一种基于改进实数型强跟踪滤波器的电力系统失真信号突变参数估计测量方法。

背景技术

参数估计对于现代电力系统的运行有着重要意义。本发明所述的参数估计主要包括幅值、频率和相位的估计，三者的估计结果分别反映了电力系统信号的变化的范围、变化的快慢以及信号开始时的情况，因此幅值、频率和相位的估计对于确定电力系统信号内部信息具有决定性作用。通常而言，在电力系统的参数估计中，频率估计是目前要解决的核心和难点。电力系统数字化保护中，继电器采用基于频率的信号处理算法，从采样得到的电压信号或电流信号中，可估计系统状态信息。在电能质量监控领域，电力信号的电压幅值和相位检测皆是基于频率测量而进行的。在分布式发电系统并网过程中，根据频率测量值可分析供电设备与电网之间的动态能量平衡关系。但是电力系统中电力电子装置的广泛使用，可控硅转换装置和高频逆变器注入非正弦电流信号，导致电力信号受到谐波和噪声信号的破坏而失真。因此快速并准确测量失真电压(电流)信号频率成为当前电力系统研究热点。

目前，国内外比较常用的方法有快速傅里叶变换，最小二乘法、牛顿递归法、自适应陷波滤波法、最小均方误差法、正交分量滤波法等方法。这些方法各具特点，只能针对某一特定的工程问题，在强噪声和谐波干扰条件下，大多表现不尽如人意。而卡尔曼滤波法能够克服噪声和谐波的影响，准确地估计信号的频率、相位和幅值。因此针对非线性系统正弦波频率估计，提出扩展复数卡尔曼滤波器应用于畸变电压信号的频率和幅值跟踪。然而由于扩展卡尔曼滤波在系统达到工况时，失去了对幅值和频率等电气参数突变工况的跟踪能力。为此提出了采用滞环比较的方法重置误差协方差阵，然而滞环比较的上下限阈值仅根据噪声的强度取经验值，自适应能力差。随着强跟踪滤波器思想的提出，采用强跟踪滤波器对电力系统失真信号参数估计，综合解决了滤波发散和跟踪突变能力不强等问题，但是在数值计算过程中，次优渐消因子的计算和分析比较需在复数域中完成，计算过程复杂，特别是针对高维过程。然而还需要注意的是，强跟踪滤波器的参数选择对跟踪结果有着极大地影响，而许多时候这些参数的设计均采用的是经验值，具有一定的盲目性。

因此，目前需要一种能够既降低运算复杂程度，又满足在线应用的实时性和有效性，并能提供快速准确跟踪的方法。

发明内容

有鉴于此，本发明的目的在于提供一种基于改进实数型强跟踪滤波器的电力系统失真信号突变参数估计测量方法，此算法利用三次采样正弦波关系模型对信号的平滑能力，以及强跟踪滤波器的自适应能力，改进了对信号突变跟踪的鲁棒性。并通过引入遗传算法对强跟踪滤波器进行改进，提高了整个算法的跟踪能力。另外，由于此算法考虑了实数型强跟踪滤波器，避免了复数型强跟踪滤波器的复杂运算过程。从而避免了复杂的方差重置规则设计，保证了此方法在线应用的实时性和有效性。

为达到上述目的，本发明提供如下技术方案：

一种基于改进实数型强跟踪滤波器的电力系统失真信号的频率、幅值、相位等参数估计方法，属于电力系统领域。该方法包括以下步骤：步骤一：获取滤波目标信号；步骤二：确定滤波目标信号特征；步骤三：针对滤波目标信号，建立滤波模型，此模型为连续三次采样正弦波电压的等式关系模型，从而得到相应的状态方程、观测方程，以及离散滤波目标信号；步骤四：根据实数型强跟踪滤波器递推规则，得到实数型强跟踪滤波器的离散算法公式；步骤五：采用遗传算法对实数型强跟着滤波器的参数设置进行在线优化求解，之后对离散滤波目标信号进行滤波和参数估计，得到相应的电力信号参数。

进一步，步骤一具体包括，通过电压互感器获取电力系统中的失真信号，并将该失真信号调制为弱电信号制，从而得到滤波目标信号y。

进一步，步骤二中，所述的滤波目标信号特征为，该信号为电力系统的所有正弦失真信号，包括三相或单相、电流或电压。

进一步，步骤三具体包括如下步骤：首先，将连续非线性的电力系统失真信号转换为离散时间信号y_k。之后，根据滤波目标信号的特征，采用连续三个采样时刻的电压等式关系构建滤波模型，包括状态方程和观测方程需要注意的是，当谐波和衰减直流分量导致信号失真时，建模所采用的等式关系以及构建得到的模型都需改变。

进一步，步骤四具体包括，在经典扩展卡尔曼滤波器算法公式的基础上，在递推状态预报误差协方差时引入次优渐消因子λ_k+1，保证每一步递推过程中，保持了不同时刻的残差序列处处正交。其工程意义在于,当存在模型不确定性或者电力信号参数突变时,在线调整增益矩阵K_k+1可使残差ε_k+1始终保持高斯白噪声的性质。该滤波器充分利用了残差序列中的有效信息,对电力信号参数变动具有较强的鲁棒性,而且在整个递推过程中所有的数据的数据类型均为实数型，从而得到实数型强跟踪滤波器的具体算法公式。

进一步，步骤五中利用改进的实数型强跟踪滤波器进行滤波和参数估计采用以下步骤：

(1)输入离散滤波目标信号y_k，设定遗忘因子ρ等各项相关参数值，其中通过遗传算法对弱化因子β进行优化求解。并根据该信号各项参数得到状态预报以及给定初始的预报误差协方差(此时设定离散时间k为1)；

(2)计算得到输出残差ε_k+1；

(3)通过此时的预报误差协方差进一步计算得到该时刻卡尔曼增益K_k，以及通过将此时的输出残差ε_k+1和遗忘因子ρ相结合运算得到该时刻的实际输出残差序列的协方差V_k+1。并进一步根据获得的卡尔曼增益K_k计算得出该时刻观测噪声协方差以及利用实际输出残差序列的协方差V_k+1和弱化因子β等数据，运算得到该时刻的次优渐消因子λ_k+1；

(4)根据当前的卡尔曼增益K_k和状态预报以及输入信号y_k等数据，计算得出该时刻的状态估计值接着利用前面几步得到的观测噪声协方差和次优渐消因子λ_k+1等数据运算得出当前的估计误差协方差

(5)时间参数k增加一个单位，之后重复(2)～(4)的步骤，直到滤波时间结束。

进一步，所述方法还包括以下步骤：为了甄别该方法的精准性，采用标准偏差STD和差异系数COV对其精准性进行评判；另外，为了确定该方法满足实时应用，通过获得该方法的计算时间来间接知道该方法的复杂程度，从而确定其是否满足实时应用要求。

本发明的有益效果在于：本发明所述的方法，利用三次采样正弦波等式关系对信号的平滑能力，以及强跟踪滤波器的自适应能力，改进了对信号突变跟踪的鲁棒性。并通过引入遗传算法对强跟踪滤波器进行改进，提高了整个算法的跟踪能力。另外，由于此算法考虑了实数型强跟踪滤波器，避免了复数型强跟踪滤波器的复杂运算过程。从而避免了复杂的方差重置规则设计，保证了此方法在线应用的实时性和有效性。

附图说明

为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行说明：

图1为本发明的流程示意图；

具体实施方式

下面将结合附图，对本发明的优选实施例进行详细的描述。

图1为本发明的流程示意图，该电力系统失真信号参数估计方法包括五个步骤：步骤一：获取滤波目标信号；步骤二：确定滤波目标信号特征；步骤三：针对滤波目标信号，建立滤波模型，此模型为连续三次采样正弦波电压的等式关系模型，从而得到相应的状态方程、观测方程，以及离散滤波目标信号；步骤四：根据实数型强跟踪滤波器递推规则，得到实数型强跟踪滤波器的离散算法公式；步骤五：采用遗传算法对实数型强跟着滤波器的参数设置进行在线优化求解，之后对离散滤波目标信号进行滤波和参数估计，得到相应的电力信号参数。

1.由于本发明针对的是电力系统的失真信号，则步骤一中所述的获取滤波目标信号具体步骤为，通过电压互感器获取电力系统中的失真信号，并将该失真信号调制为弱电信号制，从而得到滤波目标信号y。

2.在本实施例中，步骤二中所述的滤波目标信号特征，首先需明确的一点便是该信号中混有附加噪声或谐波信号，本发明主要考虑电力系统信号混有零均值的高斯白噪声情况。其次，该滤波目标信号的特征还包括，该信号可以为单相或者三相非线性正弦波电压或电流信号。

3.在本实施例中，步骤三采用的具体内容为：首先将连续非线性的电力系统失真信号转换为离散时间信号y_k。由于三次采样正弦波关系模型对信号具有较好的平滑能力，可弥补强跟踪滤波平滑度差的性质，从而本发明采用如

${\hat{y}}_k-2\cos \mathrm{\ω}{T}_s{\hat{y}}_{k-1}+{\hat{y}}_{k-2}=0---\left(1\right)$

所示的连续三个采样时刻的电压等式模型对离散滤波目标信号y_k进行构建模型，得到的信号模型为：

y_k＝Acos(kωT_s+φ)+ε_k    (2)

其中：y_k为瞬时信号值；A为信号幅值；k为采样时刻；T_s为采样周期；ω为弧频率；φ为相位；ε_k为附加噪声(假设为零均值的高斯白噪声，方差为)。

该模型包括的状态方程和观测方程分别如下所示:

${\hat{x}}_k={\left[\begin{array}{ccc}2\cos \mathrm{\ω}{T}_s& {\hat{y}}_{k-1}& {\hat{y}}_{k-2}\end{array}\right]}^T---\left(3\right)$

${\hat{y}}_k={\left[\begin{array}{ccc}0& 2\cos \mathrm{\ω}{T}_s& -1\end{array}\right]}^T{\hat{x}}_k+{\mathrm{\ϵ}}_k---\left(4\right)$

需要注意的是，当谐波和衰减直流分量导致信号失真时，建模所采用的等式关系以及构建得到的模型都将发生变化，分别变形为：

${\hat{y}}_k=2\cos \mathrm{\ω}{T}_s{\hat{y}}_{k-1}-{\hat{y}}_{k-2}+\underset{m=1}{\overset{M}{\mathrm{\Σ}}}[A_m\·\left(2\cos \mathrm{m\ω}{T}_s\right)\·\sin \{\mathrm{m\ω}(k-1)T_s+{\mathrm{\φ}}_m\left\}\right]---\left(5\right)$

$y_k=\underset{m=1}{\overset{M}{\mathrm{\Σ}}}{A}_m\cos \left(\mathrm{m\ω}{T}_s{\mathrm{\φ}}_m\right)+{\mathrm{\ϵ}}_k---\left(6\right)$

式中：y_k为k时刻信号采样值；m为谐波阶次；M为信号中谐波最高阶次；A_m为m阶谐波幅值；φ_m为m阶谐波相位；ε_k为附加噪声。

此外，需要特别指出的是，此模型中的所有参数的数据类型均为实数。

4.本发明所述参数估计方法的步骤四主要包含的内容即为根据实数型强跟踪滤波器递推规则，得到实数型强跟踪滤波器的离散算法公式。根据步骤三中给出的电力系统失真信号采样后的离散信号关系式(2)，推导出满足如式(3)所示的状态方程，以及与之对应的测量方程(4)。之后借助经典扩展卡尔曼滤波器算法公式，在递推状态预报误差协方差时引入次优渐消因子λ_k+1，并且在整个递推过程中要求所有的数据的数据类型均为实数型，从而得到实数型强跟踪滤波器的具体算法公式。相比与已有的复数型强跟踪滤波器算法，整个实数型强跟踪滤波器算法具有更低的计算复杂程度，更加有利于本发明在实际电力系统中的应用。

5.在本实施例中，步骤五采用遗传算法进行改进的实数型强跟踪滤波器，算法具体步骤如下所述：

(1)输入滤波目标信号，设定遗忘因子ρ和弱化因子β，并根据该信号各项参数得到状态预报以及给定初始的预报误差协方差(此时设定离散时间k为1)，其中遗忘因子ρ取值需满足0≤ρ≤1的要求，而弱化因子β＜1则通过将

$\mathrm{\β}:\underset{\mathrm{\β}}{\min }\left(\underset{k=0}{\overset{L}{\mathrm{\Σ}}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\right|x_{i_k}-{\hat{x}}_{i_k}\left|\right)---\left(7\right)$

设为遗传算法的适应度函数来进行在线优化求解。引入此弱化因子β的目的是为了使状态估计值更为平滑，从而将弱化因子β引入到次优渐消因子λ_k+1的计算中去，削弱λ_k+1的调节作用，避免过调节情况的出现；

(2)计算得到输出残差ε_k+1，即为瞬间信号y_k+1与估计信号之间的差值；

(3)通过预报误差协方差进一步计算得到该时刻卡尔曼增益K_k，以及通过将输出残差ε_k+1和遗忘因子ρ相结合得到该时刻的实际输出残差序列的协方差V_k+1。并进一步根据获得的卡尔曼增益K_k计算得出该时刻观测噪声协方差以及利用实际输出残差序列的协方差V_k+1和弱化因子β等数据，运算得到该时刻的次优渐消因子λ_k+1。

其中，V_k+1和λ_k+1的具体运算公式如下所示：

$V_{k+1}=\left\{\begin{array}{cc}{\mathrm{\ϵ}}_1{\mathrm{\ϵ}}_1^T,& k=0\\ \frac{\mathrm{\ρ}{V}_k+{\mathrm{\ϵ}}_{k+1}{\mathrm{\ϵ}}_{k+1}^T}{1+\mathrm{\ρ}},& k\≥1\end{array}\right.---\left(7\right)$

N_k+1＝V_k+1-βR_k+1    (8)

$M_{k+1}=H_{k+1}{F}_{k+1}{\hat{P}}_{k+1|k+1}{F}_{k+1}^T{H}_{k+1}^T---\left(9\right)$

${\mathrm{\&lambda;}}_{k+1}=\left\{\begin{array}{cc}{\mathrm{\&lambda;}}_0,& {\mathrm{\&lambda;}}_0\&GreaterEqual;1\\ 1,& {\mathrm{\&lambda;}}_0<1\end{array}\right.,{\mathrm{\&lambda;}}_0=\frac{\mathrm{tr}\left[N_{k+1}\right]}{\mathrm{tr}\left[M_{k+1}\right]}---\left(10\right)$

式(10)中，tr[·]为求矩阵迹的算子；

(4)根据当前的卡尔曼增益K_k和状态预报以及输入信号y_k等数据，计算得出该时刻的状态估计值接着利用前面几步得到的观测噪声协方差和次优渐消因子λ_k+1等数据运算得出当前的估计误差协方差

(5)时间参数k增加一个单位，之后重复(2)～(4)的步骤，直到滤波时间结束。

6.在本实施例中，所述方法还包括以下步骤：为了甄别该方法的精准性，采用标准偏差STD和差异系数COV对其精准性进行评判，二者的计算公式分别如下所示：

$\mathrm{STD}=\sqrt{\frac{\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{(x_k-\stackrel{\&OverBar;}{X})}^2}{N-1}}---\left(8\right)$

$\mathrm{COV}=\frac{\mathrm{STD}}{\stackrel{\&OverBar;}{X}}\&times;100\%---\left(9\right)$

式中：N为最大采样步数；x_k为第k个采样步数时的测量值，而则为测量值的平均值。

另外，为了确定该方法满足实时应用，通过获得该方法的计算时间来间接知道该方法的复杂程度，从而确定其是否满足实时应用要求。

通过以上六个步骤，能够为电力系统的失真信号提供有效的参数估计，特别是在信号发生突变时，能够提供快速准确的参数估计测量，从而为电力系统运行时的安全监控和有效保护提供数据基础。

最后说明的是，以上优选实施例仅用以说明本发明的技术方案而非限制，尽管通过上述优选实施例已经对本发明进行了详细的描述，但本领域技术人员应当理解，可以在形式上和细节上对其作出各种各样的改变，而不偏离本发明权利要求书所限定的范围。

Claims (7)

1.一种基于改进实数型强跟踪滤波器的电力系统失真信号突变参数估计测量方法，其特征在于：包括以下步骤：

步骤一：获取滤波目标信号；

步骤二：确定滤波目标信号特征；

步骤三：针对滤波目标信号，建立滤波模型，此模型为连续三次采样正弦波电压的等式关系模型，从而得到相应的状态方程、观测方程，以及离散滤波目标信号；

步骤四：根据实数型强跟踪滤波器递推规则，得到实数型强跟踪滤波器的离散算法公式；

步骤五：采用遗传算法对实数型强跟着滤波器的参数设置进行在线优化求解，之后对离散滤波目标信号进行滤波和参数估计，得到相应的电力信号参数。

2.根据权利要求1所述的基于遗传算法的电力系统失真信号突变参数估计测量方法，其特征在于：步骤一具体包括，通过电压互感器获取电力系统中的失真信号，并将该失真信号调制为弱电信号制，从而得到估计目标信号y。

3.根据权利要求2所述的基于改进实数型强跟踪滤波器的电力系统失真信号突变参数估计测量方法，其特征在于：步骤二中，所述的滤波目标信号特征为，该信号为电力系统的所有正弦失真信号，包括三相或单相、电流或电压。

4.根据权利要求3所述的基于改进实数型强跟踪滤波器的电力系统失真信号突变参数估计测量方法，其特征在于：步骤三具体包括如下步骤：首先，将连续非线性的电力系统失真信号转换为离散时间信号y_k。之后，根据滤波目标信号的特征，采用连续三个采样时刻的电压等式关系构建滤波模型，包括状态方程和观测方程需要注意的是，当谐波和衰减直流分量导致信号失真时，建模所采用的等式关系以及构建得到的模型都需改变。

5.根据权利要求4所述的基于改进实数型强跟踪滤波器的电力系统失真信号突变参数估计测量方法，其特征在于：步骤四具体包括，在经典扩展卡尔曼滤波器算法公式的基础上，在递推状态预报误差协方差时引入次优渐消因子λ_k+1，保证每一步递推过程中，保持了不同时刻的残差序列处处正交。其工程意义在于,当存在模型不确定性或者电力信号参数突变时,在线调整增益矩阵K_k+1可使残差ε_k+1始终保持高斯白噪声的性质。该滤波器充分利用了残差序列中的有效信息,对电力信号参数变动具有较强的鲁棒性,而且在整个递推过程中所有的数据的数据类型均为实数型，从而得到实数型强跟踪滤波器的具体算法公式。

6.根据权利要求5所述的基于改进实数型强跟踪滤波器的电力系统失真信号突变参数估计测量方法，其特征在于：步骤五中利用改进实数型强跟踪滤波器进行参数估计采用以下步骤：

(1)输入离散滤波目标信号y_k，设定遗忘因子ρ等各项相关参数值，其中通过遗传算法对弱化因子β进行优化求解。并根据该信号各项参数得到状态预报以及给定初始的预报误差协方差(此时设定离散时间k为1)；

(2)计算得到输出残差ε_k+1；

(3)通过此时的预报误差协方差进一步计算得到该时刻卡尔曼增益K_k，以及通过将此时的输出残差ε_k+1和遗忘因子ρ相结合运算得到该时刻的实际输出残差序列的协方差V_k+1。并进一步根据获得的卡尔曼增益K_k计算得出该时刻观测噪声协方差以及利用实际输出残差序列的协方差V_k+1和弱化因子β等数据，运算得到该时刻的次优渐消因子λ_k+1；

(4)根据当前的卡尔曼增益K_k和状态预报以及输入信号y_k等数据，计算得出该时刻的状态估计值接着利用前面几步得到的观测噪声协方差和次优渐消因子λ_k+1等数据运算得出当前的估计误差协方差

(5)时间参数k增加一个单位，之后重复(2)～(4)的步骤，直到滤波时间结束。

7.根据权利要求6所述的基于改进实数型强跟踪滤波器的电力系统失真信号突变参数估计测量方法，其特征在于：所述方法还包括以下步骤：为了甄别该方法的精准性，采用标准偏差STD和差异系数COV对其精准性进行评判；另外，为了确定该方法满足实时应用，通过获得该方法的计算时间来间接知道该方法的复杂程度，从而确定其是否满足实时应用要求。