CN113553545B - 小电流电器的事件检测方法及系统、设备、存储介质 - Google Patents
小电流电器的事件检测方法及系统、设备、存储介质 Download PDFInfo
- Publication number
- CN113553545B CN113553545B CN202110794271.9A CN202110794271A CN113553545B CN 113553545 B CN113553545 B CN 113553545B CN 202110794271 A CN202110794271 A CN 202110794271A CN 113553545 B CN113553545 B CN 113553545B
- Authority
- CN
- China
- Prior art keywords
- window
- event
- variance
- value
- calculation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 86
- 238000004364 calculation method Methods 0.000 claims abstract description 119
- 230000001052 transient effect Effects 0.000 claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000009825 accumulation Methods 0.000 claims description 76
- 230000001186 cumulative effect Effects 0.000 claims description 25
- 238000004590 computer program Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 11
- NRNCYVBFPDDJNE-UHFFFAOYSA-N pemoline Chemical compound O1C(N)=NC(=O)C1C1=CC=CC=C1 NRNCYVBFPDDJNE-UHFFFAOYSA-N 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 abstract description 23
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 230000008859 change Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000000605 extraction Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/02—Measuring effective values, i.e. root-mean-square values
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Complex Calculations (AREA)
- Measurement Of Current Or Voltage (AREA)
Abstract
本发明公开了一种小电流电器的事件检测方法及系统、设备、存储介质,该方法通过在电流有效值时间序列上依次定义均值计算窗口、暂态检测窗口和方差计算窗口,相比于现有的CUSUM算法不仅增加了方差计算窗口来判断序列的波动情水平,而且将方差计算窗口内的方差值与预设的方差阈值进行比对,当计算的方差值小于方差阈值时,对累积和的计算加上了系数(1+δd),针对小电流电器,虽然其电流有效值和波动值较小,但是通过系数(1+δd)可以显著提高小电流电器的投入\或切出事件累积和,有利于准确识别出小电流电器的投入/切出事件,有效降低了小电流电器事件被漏检的概率。
Description
技术领域
本发明涉及非侵入式负荷识别技术领域,特别地,涉及一种小电流电器的事件检测方法及系统、设备、计算机可读取的存储介质。
背景技术
目前,用户负荷监测按数据的获取方式分为侵入式和非侵入式。20世纪80年代,Hart提出了非侵入式负荷监测(Non-Intrusive Load Monitoring,NILM),相比于侵入式监测,NILM具有低安装成本、用户接受程度高等优势。NILM主要包含3个步骤,分别是事件检测、特征提取和负荷识别。其中特征提取和负荷识别两部分均以事件检测为基础,提高事件检测的准确性可有效增强负荷识别率。目前变点检测是更具有发展前景的检测方法,主要有4种概率方法用于检测电气量信号的变化:广义似然比、卡方拟合优度、贝叶斯信息准则和累积和(Cumulative Sum,CUSUM)方法。其中,CUSUM算法简单实用,稳定性强。
但是,在对小电流电器的投入/切出事件进行检测时,由于小电流电器的电流有效值较小,且电流波动值也小,现有的CUSUM算法经常会出现小电流电器事件被漏检的问题。并且,现有的CUSUM算法只能检测出投入/切出事件发生的暂态时刻,而无法检测出投入/切出事件发生后进入稳态的时刻,容易将暂态过程认定为稳态,对于后续的特征提取部分,无法保证稳态特征提取的准确性。
发明内容
本发明提供了一种小电流电器的事件检测方法及系统、设备、计算机可读取的存储介质,以解决现有的CUSUM算法容易漏检小电流电器的投入/切出事件的技术问题。
根据本发明的一个方面,提供一种小电流电器的事件检测方法,包括以下内容:
获取电流有效值时间序列;
在电流有效值时间序列上定义三个连续的滑动窗口,沿滑动方向排列依次为均值计算窗口、暂态检测窗口和方差计算窗口;
定义投入事件累积和和切出事件累积和的计算公式如下:
其中,为第k次计算时的投入事件累积和,/>为第k-1次计算时的投入事件累积和,/>为第k次计算时的切出事件累积和,/>为第k-1次计算时的切出事件累积和,Mm为第k次计算时的均值计算窗口的电流有效值均值,Md为第k次计算时的暂态检测窗口的电流有效值均值,β为噪声参数,d为时间延迟因子,δ为权重参数,V为第k次计算时的方差计算窗口内的电流有效值的方差值;
设置投入事件累积和的初始值和切出事件累积和的初始值均为0,并预设方差阈值Δmin、确定投入事件发生的累积和阈值H1、确定切出事件发生的累积和阈值H2;
利用三个滑动窗口在电流有效值时间序列上不断同步滑动,且每滑动一次分别计算均值计算窗口、暂态检测窗口和方差计算窗口内的电流有效值均值,并计算方差计算窗口内的电流有效值的方差值V;
基于每次窗口滑动的计算结果对投入事件累积和和切出事件累积和进行分析,从而判断出是否存在负荷投入事件和/或负荷切出事件。
进一步地,在所述基于每次窗口滑动的计算结果对投入事件累积和和切出事件累积和进行分析,从而判断出是否存在负荷投入事件和/或负荷切出事件的过程中:
当Md>Mm+β时,则进行投入事件累积和计算,直至时确定投入事件发生;当Md<Mm+β时,则进行切出事件累积和计算,直至/>时确定切出事件发生。
进一步地,当确定投入事件或切出事件发生后还包括以下内容:
继续滑动窗口并计算累积和,直至累积和不再增加,且此时方差计算窗口内的方差值小于方差阈值,则判定电流有效值时间序列进入稳态。
进一步地,所述方差阈值Δmin通过以下过程确定:
利用稳态时间窗口对小电流电器的稳态运行数据进行分析,计算得到小电流电器稳态运行时电流有效值的方差,则所述方差阈值应大于小电流电器稳态运行时的方差值。
进一步地,对于电流波动越小的电器,其权重参数δ值越大。
另外,本发明还提供一种小电流电器的事件检测系统,采用如上所述的事件检测方法,包括
数据采集单元,用于获取电流有效值时间序列;
滑动窗口设置单元,用于在电流有效值时间序列上定义三个连续的滑动窗口,沿滑动方向排列依次为均值计算窗口、暂态检测窗口和方差计算窗口;
累积和定义单元,用于定义投入事件累积和和切出事件累积和的计算公式;
初始化单元,用于设置投入事件累积和的初始值和切出事件累积和的初始值均为0,并预设方差阈值Δmin、确定投入事件发生的累积和阈值H1、确定切出事件发生的累积和阈值H2;
计算单元,用于在窗口滑动过程中计算均值计算窗口、暂态检测窗口和方差计算窗口内的电流有效值均值,并计算方差计算窗口内的电流有效值的方差值V;
分析单元,用于基于每次窗口滑动的计算结果对投入事件累积和和切出事件累积和进行分析,从而判断出是否存在负荷投入事件和/或负荷切出事件。
进一步地,还包括
稳态时刻识别单元,用于在所述分析单元确定投入事件或切出事件发生后,继续滑动窗口并计算累积和,直至累积和不再增加,且此时方差计算窗口内的方差值小于方差阈值,则判定电流有效值时间序列进入稳态。
进一步地,还包括
方差阈值设置单元,用于利用稳态时间窗口对小电流电器的稳态运行数据进行分析,计算得到小电流电器稳态运行时电流有效值的方差,则所述方差阈值应大于小电流电器稳态运行时的方差值。
另外,本发明还提供一种设备,包括处理器和存储器,所述存储器中存储有计算机程序,所述处理器通过调用所述存储器中存储的所述计算机程序,用于执行如上所述的方法的步骤。
另外,本发明还提供一种计算机可读取的存储介质,用于存储对小电流电器进行事件检测的计算机程序,所述计算机程序在计算机上运行时执行如上所述的方法的步骤。
本发明具有以下效果:
本发明的小电流电器的事件检测方法,通过在电流有效值时间序列上依次定义均值计算窗口、暂态检测窗口和方差计算窗口,相比于现有的CUSUM算法不仅增加了方差计算窗口来判断序列的波动情水平,而且将方差计算窗口内的方差值与预设的方差阈值进行比对,当计算的方差值小于方差阈值时,对累积和的计算加上了系数(1+δd),针对小电流电器,虽然其电流有效值和波动值较小,但是通过系数(1+δd)可以显著提高小电流电器的投入\或切出事件累积和,有利于准确识别出小电流电器的投入/切出事件,有效降低了小电流电器事件被漏检的概率。
另外,本发明的小电流电器的事件检测系统、设备、计算机可读取的存储介质同样具有上述优点。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明优选实施例的小电流电器的事件检测方法的流程示意图。
图2是本发明优选实施例的滑动窗滑过投入事件的发生时刻的各阶段示意图。
图3是本发明优选实施例的小电流电器的事件检测方法与现有的CUSUM算法在事件检测过程中累积和的对比示意图。
图4是是本发明另一实施例的小电流电器的事件检测系统的单元结构示意图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明,但是本发明可以由下述所限定和覆盖的多种不同方式实施。
如图1所示,本发明的优选实施例提供一种小电流电器的事件检测方法,包括以下内容:
步骤S1:获取电流有效值时间序列;
步骤S2:在电流有效值时间序列上定义三个连续的滑动窗口,沿滑动方向排列依次为均值计算窗口、暂态检测窗口和方差计算窗口;
步骤S3:定义投入事件累积和和切出事件累积和的计算公式如下:
其中,为第k次计算时的投入事件累积和,/>为第k-1次计算时的投入事件累积和,/>为第k次计算时的切出事件累积和,/>为第k-1次计算时的切出事件累积和,Mm为第k次计算时的均值计算窗口的电流有效值均值,Md为第k次计算时的暂态检测窗口的电流有效值均值,β为噪声参数,d为时间延迟因子,δ为权重参数,V为第k次计算时的方差计算窗口内的电流有效值的方差值;
步骤S4:设置投入事件累积和的初始值和切出事件累积和的初始值均为0,并预设方差阈值Δmin、确定投入事件发生的累积和阈值H1、确定切出事件发生的累积和阈值H2;
步骤S5:利用三个滑动窗口在电流有效值时间序列上不断同步滑动,且每滑动一次分别计算均值计算窗口、暂态检测窗口和方差计算窗口内的电流有效值均值,并计算方差计算窗口内的电流有效值的方差值V;
步骤S6:基于每次窗口滑动的计算结果对投入事件累积和和切出事件累积和进行分析,从而判断出是否存在负荷投入事件和/或负荷切出事件。
可以理解,本实施例的小电流电器的事件检测方法,通过在电流有效值时间序列上依次定义均值计算窗口、暂态检测窗口和方差计算窗口,相比于现有的CUSUM算法不仅增加了方差计算窗口来判断序列的波动情水平,而且将方差计算窗口内的方差值与预设的方差阈值进行比对,当计算的方差值小于方差阈值时,对累积和的计算加上了系数(1+δd),针对小电流电器,虽然其电流有效值和波动值较小,但是通过系数(1+δd)可以显著提高小电流电器的投入或切出事件累积和,有利于准确识别出小电流电器的投入/切出事件,有效降低了小电流电器事件被漏检的概率。
可以理解,所述步骤S6具体为:
当Md>Mm+β时,则进行投入事件累积和计算,直至时确定投入事件发生;当Md<Mm+β时,则进行切出事件累积和计算,直至/>时确定切出事件发生。
在每次窗口滑动的计算结果中,当暂态检测窗口的电流有效值均值Md大于均值计算窗口的电流有效值均值Mm与噪声参数β之和时,证明电流负荷数据是增加的,则可能存在电器投入,故而进行投入事件累积和计算,直至时确定投入事件发生。相反地,若暂态检测窗口的电流有效值均值Md小于均值计算窗口的电流有效值均值Mm与噪声参数β之和时,则可能存在电器切出,故而进行切出事件累积和计算,直至/>时确定切出事件发生。
可以理解,投入事件的检测原理与切出事件的检测原理相同,故本发明中后续的说明以投入事件的检测来做示范性说明,切出事件的检测不再赘述。
具体地,首先获取电流有效值时间序列其中,i(j)表示的是电流有效值时间序列I中的第j个电流有效值。
然后,在电流有效值时间序列I上定义三个连续的滑动窗口,沿滑动方向排列依次为均值计算窗口Wm、暂态检测窗口Wd和方差计算窗口Wv,其中,均值计算窗口Wm的窗口长度为m,暂态检测窗口Wd的窗口长度为n,方差计算窗口Wv的窗口长度为v。
为了检测电器投入和切出,定义投入事件累积和为:
通过上式可以有效提高电流波动小的小电流电器的累积和,并且,对于电流波动越小的电器,其权重参数δ值越大,在中当前统计值所占的比重越大,则/>累计的能力也就越强。
并且,权重参数δ可以根据实际检测效果进行优化,检测效果=正确检测个数/总个数*100%,其中,总个数=正确检测个数+错误检测个数+漏测个数,根据大量训练结果可以得到检测效果最好时的权重参数δ。
另外,所述方差阈值Δmin通过以下过程确定:
利用稳态时间窗口对小电流电器的稳态运行数据进行分析,计算得到小电流电器稳态运行时电流有效值的方差,则所述方差阈值应大于小电流电器稳态运行时的方差值。其中,窗口的长度优选为50个点,当然在本发明的其它实施例中,窗口的长度也可以选择其它数值,例如60个点、100个点等,在此不做具体限定。
然后,设置投入事件累积和的初始值为0,并预设方差阈值Δmin和确定投入事件发生的累积和阈值H1。利用三个滑动窗口在电流有效值时间序列I上进行同步滑动,每次滑动均需计算均值计算窗口、暂态检测窗口和方差计算窗口内的电流有效值均值,并计算方差计算窗口内的电流有效值的方差值V,具体计算公式如下:
其中,Md为暂态检测窗口内的电流有效值均值,Mm为均值计算窗口内的电流有效值均值,Mv为方差计算窗口内的电流有效值均值,V为方差计算窗口内电流有效值的方差值,k为均值计算窗口的第一个采样点。
当暂态检测窗口内的电流均值Md>均值计算窗口内的电流均值Mm与噪声参数β之和时,值开始从0开始累加。当/>时,说明此时存在投入事件发生的可能性,但是累积和尚未达到确定投入事件发生的程度。此时,引入时间延迟因子d,设置其初始值为0,当/>时,令d=d+1,同时判断方差计算窗口内的方差值V与方差阈值Δmin的大小。当方差值V≤方差阈值Δmin时,Md-(Mm+β)将乘上(1+δd)的系数,从而计算/>直至时确定投入事件发生,而发生时刻可以根据t=k-d逆变换得出。为避免序列震荡导致的投入事件或者切出事件的误检测,引入时间延迟因子d后,时间延迟因子每增加1,都要对/>和/>进行比较,如果/>则认为此时引起电流变化的是一个波动并非事件发生,并令/>从而避免由数据波动导致的多识别事件。
可以理解,如图2所示,在电流波动小的小电流电器的投入事件发生后,本发明的滑动窗滑过一个投入事件的发生时刻时可以分为五个阶段。
(1)阶段一:方差计算窗口Wv尚未滑动到投入事件发生时刻,此时三个窗口的电流均值均保持不变,此时,Md-Mm=0,则
(2)阶段二:事件发生时刻位于方差计算窗口Wv内,方差计算窗口内的Mv和V不断变化,而均值计算窗口和暂态检测窗口的Mm和Md保持不变,此时Md-Mm=0,
(3)阶段三:事件发生时刻位于暂态检测窗口Wd内,t1时刻为投入事件发生时刻,令I1=I0+ΔI,其中,I1为电器投入后的电流值,I0为电器投入前的电流值,ΔI为电器投入前后的电流差值,并设td=t-t1,在该阶段内,1≤td≤n,n为暂态检测窗口Wd的窗口长度。则每个时刻对应的Md和Mm为:
Mm=I0
此时存在三种情况:
1)、情况一:方差计算窗口Wv已经进入稳态,暂态检测窗口的Md不断变化,均值计算窗口和方差计算窗口的Mm、Mv、V保持不变。若则当td=n时,
2)、情况二:方差计算窗口仍有部分位于暂态过程,方差计算窗的Mv和V不断变化,直至方差计算窗口全部进入稳态,若则
此时,大于情况一。当td=n时,
3)、情况三:当td>n0(1<n0<n)后才满足此时/>介于情况一和情况二之间。
(4)阶段四:即当1≤td-n≤m-1时,每个时刻对应的Md和Mm为:
在该阶段内td对应的时刻均有:
所以在该阶段内的累计和如下式所示:
同样地,也存在三种情况:
1)、情况1:若则
当td=m+n-1时,计算出此时的累计和为
2)、情况2:若则
当td=m+n-1时,
3)、情况3:若td>n0后才满足则此时/>介于情况1和情况2之间。
(5)阶段五:三个窗口均处于稳态阶段,此时累积和不再变化。
综上所述,分别采用本发明的事件检测算法和现有的CUSUM算法,的变化曲线如图3所示。
可以理解,基于暂态特征的负荷识别算法,不仅要确定负荷启动时刻,还要确定负荷何时进入稳态,因为大部分负荷如微波炉、电磁炉等,并非启动后立刻进入稳定工作状态。而现有的CUSUM算法无法准确识别出进入稳态的时刻。作为优选的,所述事件检测方法在确定投入事件或切出事件发生后还包括以下内容:
继续滑动窗口并计算累积和,直至累积和不再增加,且此时方差计算窗口内的方差值小于方差阈值,则判定电流有效值时间序列进入稳态。
窗口在检测出有事件发生后,累积和还会继续增加,直至t=ts1时刻,累计和不再增加,且方差计算窗的方差V小于方差阈值Δmin。即当事件发生后,若在ts1至ts1+m+n+v时刻保持累积和不发生变化且电流波动较小时,即电流达到稳定,则认为ts1时刻为进入稳态的时刻。
方差计算窗口的方差值V表示电流时间序列的波动水平,因此本发明的算法可精准确定进入稳态的时刻。以检测投入事件为例,当累积和时,说明此时有进入稳态的可能,但是方差计算窗口内的电流尚未稳定,令k=k+1,计算/>和Vk,直至/>且Vk≤Δmin时,认为进入稳态的时刻ts为k。
为了验证本发明的事件检测算法与现有的CUSUM算法的区别,如下表一所示,本申请发明人进行了大量的实例分析。
表一、本发明的事件检测算法与现有的CUSUM算法的试验结果比对
由上表可知,本发明的事件检测算法与现有的滑动窗双边CUSUM事件检测算法对于投入事件/切出事件的发生时刻的检测结果相同,但现有的滑动窗双边CUSUM事件检测算法无法正确检测出投入事件/切出事件发生后进入稳态的时刻,易误将暂态过程认定为稳态,而本发明的事件检测算法则可以正确检测到投入事件/切出事件进入稳态的时刻。而正确检测稳态时刻是正确提取稳态和暂态特征的基础,也是所提取特征有效性的保证,间接提高了识别算法的训练效果。
另外,如图4所示,本发明的另一实施例还提供一种小电流电器的事件检测系统,采用如上所述的事件检测方法,所述系统包括
数据采集单元,用于获取电流有效值时间序列;
滑动窗口设置单元,用于在电流有效值时间序列上定义三个连续的滑动窗口,沿滑动方向排列依次为均值计算窗口、暂态检测窗口和方差计算窗口;
累积和定义单元,用于定义投入事件累积和和切出事件累积和的计算公式;
初始化单元,用于设置投入事件累积和的初始值和切出事件累积和的初始值均为0,并预设方差阈值Δmin、确定投入事件发生的累积和阈值H1、确定切出事件发生的累积和阈值H2;
计算单元,用于在窗口滑动过程中计算均值计算窗口、暂态检测窗口和方差计算窗口内的电流有效值均值,并计算方差计算窗口内的电流有效值的方差值V;
分析单元,用于基于每次窗口滑动的计算结果对投入事件累积和和切出事件累积和进行分析,从而判断出是否存在负荷投入事件和/或负荷切出事件。
可以理解,本实施例的小电流电器的事件检测系统,通过在电流有效值时间序列上依次定义均值计算窗口、暂态检测窗口和方差计算窗口,相比于现有的CUSUM算法不仅增加了方差计算窗口来判断序列的波动情水平,而且将方差计算窗口内的方差值与预设的方差阈值进行比对,当计算的方差值小于方差阈值时,对累积和的计算加上了系数(1+δd),针对小电流电器,虽然其电流有效值和波动值较小,但是通过系数(1+δd)可以显著提高小电流电器的投入事件累积和或者显著降低小电流电器的切出事件累积和,有利于准确识别出小电流电器的投入/切出事件,有效降低了小电流电器事件被漏检的概率。
另外,所述事件检测系统还包括稳态时刻识别单元,用于在所述分析单元确定投入事件或切出事件发生后,继续滑动窗口并计算累积和,直至累积和不再增加,且此时方差计算窗口内的方差值小于方差阈值,则判定电流有效值时间序列进入稳态。
另外,所述事件检测系统还包括方差阈值设置单元,用于利用稳态时间窗口对小电流电器的稳态运行数据进行分析,计算得到小电流电器稳态运行时电流有效值的方差,则所述方差阈值应大于小电流电器稳态运行时的方差值。
可以理解,本实施例的系统中的各个模块与上述方法实施例的各个步骤的内容相对应,故各个模块的工作原理在此不再赘述,参考上述方法实施例即可。
另外,本发明的另一实施例还提供一种设备,包括处理器和存储器,所述存储器中存储有计算机程序,所述处理器通过调用所述存储器中存储的所述计算机程序,用于执行如上所述的方法的步骤。
另外,本发明的另一实施例还提供一种计算机可读取的存储介质,用于存储对小电流电器进行事件检测的计算机程序,所述计算机程序在计算机上运行时执行如上所述的方法的步骤。
一般计算机可读取介质的形式包括:软盘(floppy disk)、可挠性盘片(flexibledisk)、硬盘、磁带、任何其与的磁性介质、CD-ROM、任何其余的光学介质、打孔卡片(punchcards)、纸带(paper tape)、任何其余的带有洞的图案的物理介质、随机存取存储器(RAM)、可编程只读存储器(PROM)、可抹除可编程只读存储器(EPROM)、快闪可抹除可编程只读存储器(FLASH-EPROM)、其余任何存储器芯片或卡匣、或任何其余可让计算机读取的介质。指令可进一步被一传输介质所传送或接收。传输介质这一术语可包含任何有形或无形的介质,其可用来存储、编码或承载用来给机器执行的指令,并且包含数字或模拟通信信号或其与促进上述指令的通信的无形介质。传输介质包含同轴电缆、铜线以及光纤,其包含了用来传输一计算机数据信号的总线的导线。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种小电流电器的事件检测方法,其特征在于,包括以下内容:
获取电流有效值时间序列;
在电流有效值时间序列上定义三个连续的滑动窗口,沿滑动方向排列依次为均值计算窗口、暂态检测窗口和方差计算窗口;
定义投入事件累积和和切出事件累积和的计算公式如下:
,
其中,为第k次计算时的投入事件累积和,/>为第k-1次计算时的投入事件累积和,/>为第k次计算时的切出事件累积和,/>为第k-1次计算时的切出事件累积和,Mm为第k次计算时的均值计算窗口的电流有效值均值,Md为第k次计算时的暂态检测窗口的电流有效值均值,β为噪声参数,d为时间延迟因子,δ为权重参数,V为第k次计算时的方差计算窗口内的电流有效值的方差值;
设置投入事件累积和的初始值和切出事件累积和的初始值均为0,并预设方差阈值Δmin、确定投入事件发生的累积和阈值H1、确定切出事件发生的累积和阈值H2;
利用三个滑动窗口在电流有效值时间序列上不断同步滑动,且每滑动一次分别计算均值计算窗口、暂态检测窗口和方差计算窗口内的电流有效值均值,并计算方差计算窗口内的电流有效值的方差值V;
基于每次窗口滑动的计算结果对投入事件累积和和切出事件累积和进行分析,从而判断出是否存在负荷投入事件和/或负荷切出事件。
2.如权利要求1所述的小电流电器的事件检测方法,其特征在于,在所述基于每次窗口滑动的计算结果对投入事件累积和和切出事件累积和进行分析,从而判断出是否存在负荷投入事件和/或负荷切出事件的过程中:
当Md>Mm+β时,则进行投入事件累积和计算,直至>H1时确定投入事件发生;当Md<Mm +β时,则进行切出事件累积和计算,直至/>>H2时确定切出事件发生。
3.如权利要求2所述的小电流电器的事件检测方法,其特征在于,当确定投入事件或切出事件发生后还包括以下内容:
继续滑动窗口并计算累积和,直至累积和不再增加,且此时方差计算窗口内的方差值小于方差阈值,则判定电流有效值时间序列进入稳态。
4.如权利要求1所述的小电流电器的事件检测方法,其特征在于,所述方差阈值Δmin通过以下过程确定:
利用稳态时间窗口对小电流电器的稳态运行数据进行分析,计算得到小电流电器稳态运行时电流有效值的方差,则所述方差阈值应大于小电流电器稳态运行时的方差值。
5.如权利要求1所述的小电流电器的事件检测方法,其特征在于,对于电流波动越小的电器,其权重参数δ值越大。
6.一种小电流电器的事件检测系统,采用如权利要求1~5任一项所述的事件检测方法,其特征在于,包括
数据采集单元,用于获取电流有效值时间序列;
滑动窗口设置单元,用于在电流有效值时间序列上定义三个连续的滑动窗口,沿滑动方向排列依次为均值计算窗口、暂态检测窗口和方差计算窗口;
累积和定义单元,用于定义投入事件累积和和切出事件累积和的计算公式;
初始化单元,用于设置投入事件累积和的初始值和切出事件累积和的初始值均为0,并预设方差阈值Δmin、确定投入事件发生的累积和阈值H1、确定切出事件发生的累积和阈值H2;
计算单元,用于在窗口滑动过程中计算均值计算窗口、暂态检测窗口和方差计算窗口内的电流有效值均值,并计算方差计算窗口内的电流有效值的方差值V;
分析单元,用于基于每次窗口滑动的计算结果对投入事件累积和和切出事件累积和进行分析,从而判断出是否存在负荷投入事件和/或负荷切出事件。
7.如权利要求6所述的小电流电器的事件检测系统,其特征在于,还包括
稳态时刻识别单元,用于在所述分析单元确定投入事件或切出事件发生后,继续滑动窗口并计算累积和,直至累积和不再增加,且此时方差计算窗口内的方差值小于方差阈值,则判定电流有效值时间序列进入稳态。
8.如权利要求6所述的小电流电器的事件检测系统,其特征在于,还包括
方差阈值设置单元,用于利用稳态时间窗口对小电流电器的稳态运行数据进行分析,计算得到小电流电器稳态运行时电流有效值的方差,则所述方差阈值应大于小电流电器稳态运行时的方差值。
9.一种电子设备,其特征在于,包括处理器和存储器,所述存储器中存储有计算机程序,所述处理器通过调用所述存储器中存储的所述计算机程序,用于执行如权利要求1~5任一项所述的方法的步骤。
10.一种计算机可读取的存储介质,用于存储对小电流电器进行事件检测的计算机程序,其特征在于,所述计算机程序在计算机上运行时执行如权利要求1~5任一项所述的方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110794271.9A CN113553545B (zh) | 2021-07-14 | 2021-07-14 | 小电流电器的事件检测方法及系统、设备、存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110794271.9A CN113553545B (zh) | 2021-07-14 | 2021-07-14 | 小电流电器的事件检测方法及系统、设备、存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113553545A CN113553545A (zh) | 2021-10-26 |
CN113553545B true CN113553545B (zh) | 2023-09-01 |
Family
ID=78131747
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110794271.9A Active CN113553545B (zh) | 2021-07-14 | 2021-07-14 | 小电流电器的事件检测方法及系统、设备、存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113553545B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109813978A (zh) * | 2018-12-25 | 2019-05-28 | 武汉中原电子信息有限公司 | 一种综合暂态特征和稳态间的变化特征的非侵入负荷式识别方法 |
CN110488128A (zh) * | 2019-09-12 | 2019-11-22 | 广东电网有限责任公司佛山供电局 | 双边累积和事件检测方法 |
CN111208351A (zh) * | 2020-01-17 | 2020-05-29 | 北京市腾河电子技术有限公司 | 基于负荷跳变计算供电线路阻抗的方法、存储介质 |
CN111308185A (zh) * | 2020-02-28 | 2020-06-19 | 宁波三星医疗电气股份有限公司 | 一种非侵入式负荷识别方法 |
CN112396007A (zh) * | 2020-11-24 | 2021-02-23 | 国网江苏省电力有限公司营销服务中心 | 一种非侵入式居民负荷突变事件三阈值检测方法和系统 |
CN112782465A (zh) * | 2020-12-01 | 2021-05-11 | 国网浙江省电力有限公司营销服务中心 | 一种基于滑动窗口技术的暂态负荷启停事件快速捕捉方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10139437B2 (en) * | 2014-02-18 | 2018-11-27 | Encored Technologies, Inc. | Apparatus, server, system and method for energy measuring |
-
2021
- 2021-07-14 CN CN202110794271.9A patent/CN113553545B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109813978A (zh) * | 2018-12-25 | 2019-05-28 | 武汉中原电子信息有限公司 | 一种综合暂态特征和稳态间的变化特征的非侵入负荷式识别方法 |
CN110488128A (zh) * | 2019-09-12 | 2019-11-22 | 广东电网有限责任公司佛山供电局 | 双边累积和事件检测方法 |
CN111208351A (zh) * | 2020-01-17 | 2020-05-29 | 北京市腾河电子技术有限公司 | 基于负荷跳变计算供电线路阻抗的方法、存储介质 |
CN111308185A (zh) * | 2020-02-28 | 2020-06-19 | 宁波三星医疗电气股份有限公司 | 一种非侵入式负荷识别方法 |
CN112396007A (zh) * | 2020-11-24 | 2021-02-23 | 国网江苏省电力有限公司营销服务中心 | 一种非侵入式居民负荷突变事件三阈值检测方法和系统 |
CN112782465A (zh) * | 2020-12-01 | 2021-05-11 | 国网浙江省电力有限公司营销服务中心 | 一种基于滑动窗口技术的暂态负荷启停事件快速捕捉方法 |
Non-Patent Citations (1)
Title |
---|
基于小电流电器的滑动窗双边CUSUM事件检测改进算法;黄伟; 郑明杨;水电能源科学;第40卷(第2期);206-210 * |
Also Published As
Publication number | Publication date |
---|---|
CN113553545A (zh) | 2021-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7216063B2 (en) | Method and apparatus for comparing a data set to a baseline value | |
CN110750429A (zh) | 运维管理系统的异常检测方法、装置、设备及存储介质 | |
CN111444168A (zh) | 一种配电室变压器日最大负荷异常数据检测处理方法 | |
CN113225209B (zh) | 一种基于时间序列相似检索的网络监控实时预警方法 | |
CN117235557B (zh) | 基于大数据分析的电气设备故障快速诊断方法 | |
CN117407700A (zh) | 一种带电作业过程中作业环境的监测方法 | |
CN107276779B (zh) | 一种监控方法、系统及设备 | |
CN110175353A (zh) | 电梯运行指标趋势分析方法 | |
CN111161097A (zh) | 基于假设检验的事件检测算法检测开关事件的方法和装置 | |
CN113553545B (zh) | 小电流电器的事件检测方法及系统、设备、存储介质 | |
CN109214318B (zh) | 一种寻找非稳态时间序列微弱尖峰的方法 | |
CN114076841A (zh) | 基于用电信息数据的窃电行为识别方法及系统 | |
CN113642421A (zh) | 周期性多状态电器的事件检测方法及系统、设备、存储介质 | |
CN114896302B (zh) | 小负荷用户的户变关系识别方法及系统、设备、存储介质 | |
CN110298767A (zh) | 一种火力发电厂时间序列变量异常监测方法及系统 | |
CN115144172A (zh) | 一种基于声音的gis断路器故障在线监测系统及方法 | |
CN111103490B (zh) | 一种负荷辨识事件检测方法、装置及设备 | |
CN114005256A (zh) | 基于报警持续时间概率函数的报警延迟器设计方法及系统 | |
CN116595338B (zh) | 基于物联网工程信息采集处理系统 | |
CN117913999B (zh) | 智能配电柜及智能配电方法 | |
CN112232580A (zh) | 一种供电中断损失分析方法和装置 | |
CN110907544A (zh) | 一种变压器油中溶解气体含量异常阶跃数据识别方法 | |
CN118033409B (zh) | 一种gcb灭弧室开关电阻测试方法 | |
CN117825899B (zh) | 一种电容器组的绝缘状态监测方法 | |
CN117150325B (zh) | 一种异常残差数据的处理方法、装置、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |