一种基于极值点的自适应事件检测方法与流程

1.本发明涉及非侵入式负荷监测领域，具体涉及一种基于极值点的自适应事件检测方法。


背景技术：

2.能源是当今社会运作和发展不可或缺的资源，而随着全球经济的快速发展，社会各个行业和家庭用户对电能数量和质量的需求日渐增长。为了节约电能，合理安排用电计划，缓解社会能源压力、实现能源的可持续发展，提出了非侵入式负荷监测。现有的大部分基于事件的nilm算法的第一步就是对总用电曲线进行事件检测。
3.近些年，国内外学者对提出了许多针对用电负荷事件检测的方法。负荷工作状态的变化可称之为事件。事件的发生可以由设备用电曲线的工作特征变化反应出来，比如有功功率、无功功率、电压和电流等。相应的，这些变化特征也可以从总用电曲线中反映出来。根据发生时的负荷特征可以实现负荷识别和监测工作。常用的事件检测算法是基于窗口的算法，比如广义似然比检验glr、拟合优度检验gof和滑动窗的双边cusum算法。此类算法均需要设置事件检测窗口的大小和较多的参数变量，对于不同采样频率的数据需要相应更改窗口大小。


技术实现要素：

4.本发明针对现有技术存在的问题提供一种可以自适应检测事件的启停时刻的基于极值点的自适应事件检测方法。
5.本发明采用的技术方案是：一种基于极值点的自适应事件检测方法，包括以下步骤：
6.步骤1：获取总功率信号，对总功率信号进行去噪；
7.步骤2：在步骤1获取的信号对应曲线中获取最大极值点和最小极值点；
8.步骤3：根据步骤2得到的极值点得到上升区间和下降区间；
9.步骤4：在步骤3得到上升区间中对上升沿进行正反向搜索，在下降区间中下降沿进行正反向搜索；
10.步骤5：若步骤4同一个区间的正反向搜索均具有满足设定条件的点，则表示检测到事件；否则退出。
11.进一步的，所述步骤1中采用mallat塔式分解方法对总功率信号进行去噪；
12.首先对噪声进行强度估计：
[0013][0014]
式中：σj为第j层的噪声强度估计，δ
*
为经验系数，dj(i)为第j层的高频系数，n为高频系数的维数；
[0015]
计算第j层的噪声阈值tj：
[0016][0017]
若第j层的高频系数的绝对值小于阈值tj则将其置零，反之收缩至tj；
[0018]
对信号进行小波分解得到高频系数和低频系数，采用上述计算得到的阈值对高频系数进行修正；对处理后的吸波系数进行逆变换得到重构的信号，即可完成信号去噪处理。
[0019]
进一步的，所述步骤2中极大值点满足以下要求：
[0020][0021]
极小值点满足以下要求：
[0022][0023]
其中：xn、x
n-1
，x
n+1
均为曲线中的点。
[0024]
进一步的，步骤3中上升区间和下降区间的确定方法如下：
[0025]
设xi、xj、xk为三个相邻的极值点，i《j《k且xi和xk为极小值点，xj为极大值点；
[0026]
若x
j-xi＞g
p
，那么上升区间为
[0027]
若x
j-xk＞g
p
，那么下降区间为
[0028]
其中，g
p
为给定的功率门限。
[0029]
进一步的，所述步骤4中在上升区间对上升沿进行正反向搜索过程如下：
[0030]
正向搜索：从n＝i搜索到n＝j，若xn满足x
n+1-xn＞g
σ
，则停止搜索记录该点对应的时间为上升沿的时间起点，否则此次搜索无果；
[0031]
反向搜索：从n＝j搜索到n＝i，若xn满足x
n-x
n-1
＞g
σ
，则停止搜索记录该点对应的时间为上升沿的时间终点，否则此次搜索无果；
[0032]
其中g
σ
为噪声门限。
[0033]
进一步的，所述步骤4中在下降区间对下降沿进行正反向搜索过程如下：
[0034]
正向搜索：从n＝j搜索到n＝k，若xn满足x
n-x
n+1
＞g
σ
，则停止搜索记录该点对应的时间为上升沿的时间起点，否则此次搜索无果；
[0035]
反向搜索：从n＝k搜索到n＝j，若xn满足x
n-1-xn＞g
σ
，则停止搜索记录该点对应的时间为上升沿的时间终点，否则此次搜索无果；
[0036]
其中g
σ
为噪声门限。
[0037]
本发明的有益效果是：
[0038]
(1)本发明首先对信号进行小波软阈值去噪处理降低噪声干扰，有效避免了噪声波动可能引起的事件误判。
[0039]
(2)本发明不需要设置事件检测窗口的大小，只需要设置功率参数门限和噪声门限即可从总用电曲线中自适应地检测事件的启停时刻，剔除了缓变功率的信号的干扰，有较强的有效性和实用性。
附图说明
[0040]
图1为本发明的流程示意图。
[0041]
图2为实施例中某家庭总用电曲线。
[0042]
图3为实施例中事件检测输出结果。
具体实施方式
[0043]
下面结合附图和具体实施例对本发明做进一步说明。
[0044]
如图1所示，一种基于极值点的自适应事件检测方法，包括以下步骤：
[0045]
步骤1：获取总功率信号，对总功率信号进行去噪；利用mallat塔式分解方法对总功率信号进行小波软阈值去噪。
[0046]
首先对信号进行小波分解，得到一组高频系数d＝{d1,d2,
…
,dm}和一组低频系数。由于高频部分不仅包含了噪声信号，还隐藏了许多重要的时频特征，需要合理的选择阈值将高频系数有选择性地减少或置零，从而降低对信号细节特征的影响，减少原始负荷特征的丢失。针对每一层高频系数的修正。
[0047]
首先对噪声进行强度估计：
[0048][0049]
式中：σj为第j层的噪声强度估计，δ
*
为经验系数，dj(i)为第j层的高频系数，n为高频系数的维数；
[0050]
计算第j层的噪声阈值tj：
[0051][0052]
若第j层的高频系数的绝对值小于阈值tj则将其置零，反之收缩至tj；
[0053]
对信号进行小波分解得到高频系数和低频系数，采用上述计算得到的阈值对高频系数进行修正；对处理后的吸波系数进行逆变换得到重构的信号，即可完成信号去噪处理。
[0054]
步骤2：在步骤1获取的信号对应曲线中获取最大极值点和最小极值点；
[0055]
极大值点满足以下要求：
[0056]
[0057]
极小值点满足以下要求：
[0058][0059]
其中：xn、x
n-1
，x
n+1
均为曲线中的点。
[0060]
步骤3：根据步骤2得到的极值点得到上升区间和下降区间；
[0061]
上升区间和下降区间的确定方法如下：
[0062]
设xi、xj、xk为三个相邻的极值点，i《j《k且xi和xk为极小值点，xj为极大值点；
[0063]
若x
j-xi＞g
p
，那么上升区间为
[0064]
若x
j-xk＞g
p
，那么下降区间为
[0065]
其中，g
p
为给定的功率门限。
[0066]
步骤4：在步骤3得到上升区间中对上升沿进行正反向搜索，在下降区间中下降沿进行正反向搜索；
[0067]
在上升区间对上升沿进行正反向搜索过程如下：
[0068]
正向搜索：从n＝i搜索到n＝j，若xn满足x
n+1-xn＞g
σ
，则停止搜索记录该点对应的时间为上升沿的时间起点，否则此次搜索无果；
[0069]
反向搜索：从n＝j搜索到n＝i，若xn满足x
n-x
n-1
＞g
σ
，则停止搜索记录该点对应的时间为上升沿的时间终点，否则此次搜索无果；
[0070]
其中g
σ
为噪声门限。
[0071]
在下降区间对下降沿进行正反向搜索过程如下：
[0072]
正向搜索：从n＝j搜索到n＝k，若xn满足x
n-x
n+1
＞g
σ
，则停止搜索记录该点对应的时间为上升沿的时间起点，否则此次搜索无果；
[0073]
反向搜索：从n＝k搜索到n＝j，若xn满足x
n-1-xn＞g
σ
，则停止搜索记录该点对应的时间为上升沿的时间终点，否则此次搜索无果；
[0074]
其中g
σ
为噪声门限。
[0075]
步骤5：若步骤4同一个区间的正反向搜索均具有满足设定条件的点，则表示检测到事件；否则退出。
[0076]
若同一个区间的正反向搜索均有满足条件的点表示成果检测到事件，根据上述搜索结果得到事件的起点和终点时刻。否则只是用电器功率缓变过程，对所有事件区间完成搜索，即成功总用电曲线事件检测。事件检测就是负荷投切的过程，也是功率跳变的过程，要记录起点和终点也就是检测整个投切事件的开始和结束过程。
[0077]
选取开源redd数据库中house 2某一天内的用电数据进行总用电器曲线的事件检测，检测前数据如图2所示，经过上述方法检测后曲线如图3所示。从图中可以看出，可以看出，算法从上升区间中搜索到了事件的开启时刻，下降区间搜索到了事件的结束时刻，完整地从总用电曲线中检测到了负荷投切的过程，无需设置检测窗口的大小，不需要较高的计
算复杂度，参数设置也较为简单。
[0078]
本发明首先对信号进行了去噪处理降低噪声干扰，且不需要设置事件检测窗口的大小。只需提前设置两个参数便可以从总用电曲线中自适应地检测事件的启停时刻，有较强的有效性和实用性。