一种高噪声背景下电能质量扰动信号降噪方法与流程

本发明涉及电能质量扰动源识别与治理技术领域，特别是一种高噪声背景下电能质量扰动信号降噪方法。

背景技术：

在现代科技发展的今天，各种电力电子器件的应用，给电网带来巨大的挑战。同时，各种精密仪器以及高精度控制装置的使用，也对电网的供电质量提出了更高的要求。各种家用电器的使用，给电网带来了嘈杂的背景，作为电能质量扰动源识别以及治理的基础，电能质量复合扰动降噪就显得格外重要。

在降噪方面有：中值滤波器、平均值滤波器、灰度级形态学、双曲s变换降噪等。中值滤波器、平均值滤波器均不适合处理条纹状噪声。灰度级形态学对处理条纹状噪声效果较好，但算法复杂，并且并不能消除信号的时频通域内噪声对信号影响。双曲s变换降噪，需要将含噪声的信号与原始信号进行对比，以此确定信号的时频通域，在实际应用中并不可行，并且不能减少噪声对基频和谐波的信号干扰。

技术实现要素：

本发明要解决的技术问题是提供一种高噪声背景下电能质量扰动信号降噪方法，消除了噪声对信号的影响，同时最大程度的压制了噪声，并且减少了降噪时给信号带来的损失。

为解决上述技术问题，本发明所采取的技术方案如下。

一种高噪声背景下电能质量扰动信号降噪方法，其特征在于包括以下步骤：

a、利用变电站中的电压、电流互感器及二次设备对相关扰动信号进行在线采集和记录；

b、将扰动信号中的中高频噪声进行抑制；

c、通过区分噪声和暂态扰动，消除噪声对信号的影响；

d、通过区分噪声和信号，对噪声进行压制。

作为优选，步骤b中，将扰动信号利用多分辨率双曲s变换，令双曲s变换的离散表达式如下

当n≠0时

当n＝0(相当于零频率)时

式中j,m,n＝0,1,...,n-1；gh(m,n)的离散化表达式：

n为进行双曲s变换的离散信号的点数；t采集信号相邻两点的时间，即采样周期；f为频率；t是时间；m用来实现在不同频率变换；参数τ用来控制时间轴上双曲窗口的位置；信号通过快速傅里叶变换得到加入参数m得到gh(m,n)为双曲窗口函数；决定了双曲窗口前半曲线的衰减度；决定了双曲窗口后半曲线的衰减度；λhy确定了双曲窗口的曲率；的增加会使双曲窗口的前后两半曲线衰减，从而改变了双曲窗口的形状；通过调整值的和，使窗宽随频率成反比变化的速度，达到调节时频精度的目的；

多分辨率双曲s变换的结果为复数矩阵，对矩阵每一个元素取模得到多分辨率双曲s变换的模矩阵；模矩阵列表示采样时间，行表示频率；基频幅值曲线，即模矩阵频率为50hz处的行向量；频率幅值曲线，即分模矩阵频中每行的最大值。

作为优选，低频部分频率范围为1hz≤fl≤100hz，多分辨率双曲s变换离散表达式的nl的取值范围为nt≤nl≤100nt时的参数为λhy＝3。

作为优选，中频部分频率范围为101hz≤fm≤900hz，多分辨率双曲s变换离散表达式的nm的取值范围为101nt≤nm≤900nt时的参数为λhy＝16。

作为优选，高频部分频率范围为901hz≤fh≤1300hz，离散表达式的nh的取值范围为901nt≤nh≤1300nt时的参数为λhy＝32。

作为优选，步骤c中，根据噪声和扰动在时频域表现不同，以模矩阵行的平均值作为基准线，噪声以及扰动使得模矩阵的行向量沿着采样点方向上下波动，在模矩阵的行向量上由于扰动的能量大于噪声的能量，由此可以通过波动能量密度来区分噪声和暂态扰动。

作为优选，

式(1)中，s^*(i,j)表示mhst模矩阵，a，b表示模矩阵的行向量与基准线相交同一波峰或波谷的端点；通过公式(1)和阈值的比较，以此确定暂态扰动的时频通域；阈值的设置和采样频率以及mhst的窗口函数相关；

在确定暂态信号扰动的时频通域，通过设置滤波器，消除噪声对信号的干扰，主要消除对谐波以及基频信号的干扰，对噪声起到一定的压制作用；滤波器表达式为，

式(2)中，为mhst的时频复数矩阵，为消噪时频滤波器，为滤除噪声后的时频复数矩阵；滤波器ff表达式如下：

式(3)中q为暂态扰动的时频通域；用公式(3)对暂态扰动赋权值1进行保留，对非暂态扰动区域的点赋权值矫正该区域点的幅值至该点所在模矩阵的行向量均值；对复数矩阵低频部分中的基频发生暂降的区域、高频部分中的暂态震荡区域赋权值1进行保留；主要对复数矩阵低频部分中的基频未发生扰动的区域、中频部分中谐波区域点赋权值消除噪声的影响；对非信号的时频通域内的噪声点赋权值对噪声起到一定的压制作用。

作为优选，步骤d中，根据噪声和信号在时频域的表现不同，含有信号的通域的能量密度是明显大于不含信号通域，由确定信号的时频通域和非信号的时频通域。

作为优选，

式(4)中，a、b间隔为3个采样点，其中n采样点，中括号为向下取整；由式(4)与阈值的比较可以确定信号时频通域与非信号时频通域，需要压制非信号的时频通域内的噪声，除噪时频滤波器ft的表达式为，

式(5)中r为信号的时频通域，用公式(5)对噪声赋权值0进行压制，对信号区域的点赋权值1进行保留。

采用上述技术方案所带来的有益效果在于：本发明在降噪时，消除了噪声对信号的影响，同时最大程度的压制了噪声，并且减少了降噪时给信号带来的损失。

附图说明

图1为本发明的流程图。

图2多分辨率双曲s变换抑制中高频部分噪声效果图。

图3为消除噪声对信号的影响的基频幅值效果图。

图4为压制噪声的模矩阵灰度效果图。

具体实施方式

参照图1-4，本发明一个具体实施方式包括以下步骤，

(1)电能质量扰动信号采集

利用变电站中的电压、电流互感器及二次设备对相关扰动信号进行采集和记录

(2)将扰动信号利用多分辨率双曲s变换及其抑制中高频部分噪声

令双曲s变换的离散表达式如下

当n≠0时

当n＝0(相当于零频率)时

式中j,m,n＝0,1,...,n-1。gh(m,n)的离散化表达式：

n为进行双曲s变换的离散信号的点数；t采集信号相邻两点的时间，即采样周期；f为频率；t是时间；m用来实现在不同频率变换；参数τ用来控制时间轴上双曲窗口的位置。信号通过快速傅里叶变换得到加入参数m得到gh(m,n)为双曲窗口函数；决定了双曲窗口前半曲线的衰减度；决定了双曲窗口后半曲线的衰减度；λhy确定了双曲窗口的曲率。的增加会使双曲窗口的前后两半曲线衰减，从而改变了双曲窗口的形状。通过调整值的和，使窗宽随频率成反比变化的速度，达到调节时频精度的目的。

多分辨率双曲s变换的结果为复数矩阵，对矩阵每一个元素取模得到多分辨率双曲s变换的模矩阵。模矩阵列表示采样时间，行表示频率。基频幅值曲线，即模矩阵频率为50hz处的行向量。频率幅值曲线，即分模矩阵频中每行的最大值。

低频部分频率范围为1hz≤fl≤100hz，多分辨率双曲s变换离散表达式的nl的取值范围为nt≤nl≤100nt的参数为λhy＝3。

中频部分频率范围为101hz≤fm≤900hz，多分辨率双曲s变换离散表达式的nm的取值范围为101nt≤nm≤900nt的参数为λhy＝16。采用单一分辨率窗口函数的广义s变换(gst)以及双曲s变换(hst)，由于设置参数时要兼顾到时域的精度，会导致噪声幅值在模矩阵中沿着频率方向不断增加，给降噪带来困难。在当前窗口函数下，多分辨率双曲s变换(mhst)，为含高噪声(信噪比为20db)的无扰动信号变换的频率幅值曲线中频沿着频率方向幅值大小基本没有增长趋势，而单一分辨率窗口函数的gst以及双曲s变换hst的频率幅值曲线中频部分沿着频率方向幅值大小增长趋势明显。在中频部分，mhst的基频幅值曲线在各频率的幅值明显小于hst、gst基频幅值曲线，体现了mhst在中频部分优异的抑制噪声的能力。

高频部分频率范围为901hz≤fh≤1300hz，离散表达式的nh的取值范围为901nt≤nh≤1300nt的参数为λhy＝32。在当前窗口函数下，多分辨率双曲s变换(mhst)，为含高噪声(信噪比为20db)的无扰动信号变换的频率幅值曲线中频沿着频率方向幅值大小基本没有增长趋势，而单一分辨率窗口函数的gst以及双曲s变换hst的频率幅值曲线高频部分沿着频率方向幅值大小增长趋势明显。在高频部分，mhst的基频幅值曲线在各频率的幅值明显小于hst、gst基频幅值曲线，体现了mhst在高频部分优异的抑制噪声的能力。

(3)消除噪声对信号的影响

mhst优异的抑噪性能，但是噪声给信号带来的干扰并没有消除。下面根据噪声和扰动在时频域表现不同，提出了基于波动能量密度降噪的方法。由于噪声具有随机性，在某一频率下，沿着采样点方向的能量是小于由扰动产生的能量。以模矩阵行的平均值作为基准线，噪声以及扰动使得模矩阵的行向量会沿着采样点方向上下波动，在模矩阵的行向量上由于扰动的能量大于噪声的能量，由此可以通过波动能量密度来区分mhst复数矩阵内的噪声和扰动。

式(1)中，s^*(i,j)表示mhst模矩阵，a，b表示模矩阵的行向量与基准线相交同一波峰(波谷)的端点。通过公式(1)和阈值的比较，以此确定暂态扰动的时频通域。阈值的设置和采样频率以及mhst的窗口函数有关。

在确定暂态信号扰动的时频通域，需要设置滤波器，消除噪声对信号的干扰，主要消除对谐波以及基频信号的干扰，对噪声起到一定的压制作用。滤波器表达式如式(2)。

式(2)中，为mhst的时频复数矩阵，为消噪时频滤波器，为滤除噪声后的时频复数矩阵。滤波器ff表达式如下：

式(3)中q为暂态扰动的时频通域。用公式(3)对暂态扰动赋权值1进行保留，对非暂态扰动区域的点赋权值矫正该区域点的幅值至该点所在模矩阵的行向量均值。对复数矩阵低频部分中的基频发生暂降(中断、暂升等)的区域、高频部分中的暂态震荡区域赋权值1进行保留。主要对复数矩阵低频部分中的基频未发生扰动的区域、中频部分中谐波区域点赋权值消除噪声的影响。对非信号的时频通域(由式4确定)内的噪声点赋权值对噪声起到一定的压制作用。

(4)压制噪声

通过上述滤波之后，并不能最大限度的滤除噪声，根据噪声和信号在时频域的表现不同，提出了基于能量密度的降噪方法。由mhst时频模矩阵可知，在含有信号的通域的能量密度是明显大于不含信号通域，由此可以确定信号的时频通域和非信号的时频通域。

式(13)中，a、b间隔为3个采样点，其中n采样点，中括号为向下取整。由式(4)与阈值的比较可以确定信号时频通域与非信号时频通域，需要压制非信号的时频通域内的噪声，除噪时频滤波器ft的表达式如式(5)。

式(5)中r为信号的时频通域，用公式(5)对噪声赋权值0进行压制，对信号区域的点赋权值1进行保留。

以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。