CN111830405A - 一种利用频率差异量的负荷开关事件检测方法和*** - Google Patents

一种利用频率差异量的负荷开关事件检测方法和*** Download PDF

Info

Publication number
CN111830405A
CN111830405A CN202010671056.5A CN202010671056A CN111830405A CN 111830405 A CN111830405 A CN 111830405A CN 202010671056 A CN202010671056 A CN 202010671056A CN 111830405 A CN111830405 A CN 111830405A
Authority
CN
China
Prior art keywords
signal
sequence
nth
transform coefficient
coefficient vector
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.)
Withdrawn
Application number
CN202010671056.5A
Other languages
English (en)
Inventor
翟明岳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Petrochemical Technology
Original Assignee
Guangdong University of Petrochemical Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong University of Petrochemical Technology filed Critical Guangdong University of Petrochemical Technology
Priority to CN202010671056.5A priority Critical patent/CN111830405A/zh
Publication of CN111830405A publication Critical patent/CN111830405A/zh
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/327—Testing of circuit interrupters, switches or circuit-breakers
    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J13/00—Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
    • H02J13/00032—Systems characterised by the controlled or operated power network elements or equipment, the power network elements or equipment not otherwise provided for
    • H02J13/00036—Systems characterised by the controlled or operated power network elements or equipment, the power network elements or equipment not otherwise provided for the elements or equipment being or involving switches, relays or circuit breakers

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Remote Monitoring And Control Of Power-Distribution Networks (AREA)

Abstract

本发明的实施例公开一种利用频率差异量的负荷开关事件检测方法和***，所述方法包括：步骤101获取按时间顺序采集的信号序列S；步骤102生成N个信号一阶差分序列；步骤103生成N个信号二阶差分序列；步骤104求取DFT变换系数矢量；步骤105求取归一化DFT变换系数矢量；步骤106求取N个频率差异量；步骤107求取事件判断阈值；步骤108判断负荷开关事件。

Description

一种利用频率差异量的负荷开关事件检测方法和***
技术领域
本发明涉及电力领域，尤其涉及一种负荷开关事件检测方法和***。
背景技术
随着智能电网的发展，家庭用电负荷的分析变得越来越重要。通过用电负荷的分析，家庭用户可以及时获得每个电器的用电信息，以及电费的精细化清单；电力部门可以获得更详尽的用户用电信息，并可以提高用电负荷预测的准确度，为电力部门提供统筹规划的依据。同时，利用每个电器的用电信息，可获知用户的用电行为，这对于家庭能耗评估和节能策略的研究具有指导意义。
当前用电负荷分解主要分为侵入式负荷分解和非侵入式负荷分解两种方法。非侵入式负荷分解方法不需要在负荷的内部用电设备上安装监测设备，只需要根据用电负荷总信息即可获得每个用电设备的负荷信息。非侵入式负荷分解方法具有投入少、方便使用等特点，因此，该方法适用于家庭负荷用电的分解。
非侵入式负荷分解算法中，电气设备的开关事件检测是其中最重要的环节。最初的事件检测以有功功率P的变化值作为事件检测的判断依据，方便且直观。这是因为任何一个用电设备的运行状态发生变化，其所消耗的功率值也必然发生改变，并且该改变也将会在所有电器所消耗的总功率中体现出来。这种方法除了需要设置功率变化值的合理阈值，还需要解决事件检测方法在实际应用中存在的问题：某些电器启动时刻的瞬时功率值会出现较大的尖峰(例如，马达启动电流远大于额定电流)，会造成电器稳态功率变化值不准确，从而影响对开关事件的判断，这种尖峰其实就是脉冲噪声；而且不同家用电器的暂态过程或长或短(脉冲噪声的持续时间和发生频率相差较大)，因此功率变化值的确定变得较为困难；由于电能质量的变化(如电压突降)有功功率会出现突变的情况，这样很可能会出现误判。(脉冲)噪声的强度很大，并且背景噪声对开关事件的正确检测有较大的影响。
现在常用的负载开关事件往往利用功率数据的变化来确定：当功率变化值超过预置的阈值时，则认为有负载开关事件发生。这种方法虽然简单易行，但由于脉冲噪声以及非线性负载的普遍使用，造成开关事件检测精度的极大下降。
因此，开关事件检测过程中，如何提高开关事件检测精度是非常重要的。负载开关事件检测是能量分解中最为重要的一步，既要检测到事件发生，同时还能确定事件发生的时刻。但是开关事件检测的精度受功率信号(功率序列)中噪声的影响较大，尤其是功率信号中普遍存在着脉冲噪声，进一步影响了检测精度。因此有效提高负载开关事件的检测精度是当前非常重要的一项工作。
发明内容
现在常用的负载开关事件往往利用功率数据的变化来确定：当功率变化值超过预置的阈值时，则认为有负载开关事件发生。这种方法虽然简单易行，但由于脉冲噪声以及非线性负载的普遍使用，造成开关事件检测精度的极大下降。
本发明的目的是提供一种利用频率差异量的负荷开关事件检测方法和***，所提出的方法利用了不同负荷开关动作时所产生的暂态功率信号差异，通过频率差异量性质区分不同负荷的不同开关事件以及异常功率信号造成的异常事件。所提出的方法具有较好的开关事件检测性能，计算也非常简单。
为实现上述目的，本发明提供了如下方案：
一种利用频率差异量的负荷开关事件检测方法，包括：
步骤101获取按时间顺序采集的信号序列S；
步骤102生成N个信号一阶差分序列，具体为：第n个信号一阶差分序列记为
Figure BDA0002582314070000021
其第i个元素的值为s|i+1|N-si；其中，i为元素信号，其取值范围为i＝1,···,N；n为差分序列序号，其取值范围为n＝1,2,···,N；N为所述信号序列S的长度；si为所述信号序列S的第i个元素；s|i+1|N为所述信号序列S的第|i+1|N个元素；|i+1|N表示以N为模取余数运算；
步骤103生成N个信号二阶差分序列，具体为：第n个信号二阶差分序列记为
Figure BDA0002582314070000022
其第i个元素的值为s|i+2|N-si；其中，s|i+2|N为所述信号序列S的第|i+2|N个元素；
步骤104求取DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure BDA0002582314070000023
的DFT变换系数矢量为
Figure BDA0002582314070000024
求取公式为
Figure BDA0002582314070000025
第n个所述信号二阶差分序列
Figure BDA0002582314070000026
的DFT变换系数矢量为
Figure BDA0002582314070000027
求取公式为
Figure BDA0002582314070000028
步骤105求取归一化DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure BDA0002582314070000029
的DFT变换系数矢量
Figure BDA00025823140700000210
的归一化DFT变换系数矢量为
Figure BDA00025823140700000211
其计算公式为
Figure BDA00025823140700000212
第n个所述信号二阶差分序列
Figure BDA00025823140700000213
的DFT变换系数矢量
Figure BDA00025823140700000214
的归一化DFT变换系数矢量为
Figure BDA00025823140700000215
其计算公式为
Figure BDA00025823140700000216
步骤106求取N个频率差异量，具体为：第n个频率差异量记为Hn，其计算公式为其中，
Figure BDA00025823140700000218
为所述信号一阶差分序列
Figure BDA00025823140700000219
的中心频率；
Figure BDA00025823140700000220
为所述信号二阶差分序列
Figure BDA00025823140700000221
的中心频率；D0为J差异因子，其计算公式为
Figure BDA00025823140700000222
其中，In是维数为n的全1列矢量；Δf为所述信号序列的采样频率；
步骤107求取事件判断阈值，具体为：事件判断阈值记为ε0，其计算公式为
Figure BDA00025823140700000223
Figure BDA00025823140700000224
步骤108判断负荷开关事件，具体为：如果第n个频率差异量Hn大于或者等于所述事件判断阈值ε0，则在所述信号序列的第n点处检测到负荷开关事件；否则，未检测到负荷开关事件。
一种利用频率差异量的负荷开关事件检测***，包括：
模块201获取按时间顺序采集的信号序列S；
模块202生成N个信号一阶差分序列，具体为：第n个信号一阶差分序列记为
Figure BDA0002582314070000031
其第i个元素的值为s|i+1|N-si；其中，i为元素信号，其取值范围为i＝1,···,N；n为差分序列序号，其取值范围为n＝1,2,···,N；N为所述信号序列S的长度；si为所述信号序列S的第i个元素；s|i+1|N为所述信号序列S的第|i+1|N个元素；|i+1|N表示以N为模取余数运算；
模块203生成N个信号二阶差分序列，具体为：第n个信号二阶差分序列记为
Figure BDA0002582314070000032
其第i个元素的值为s|i+2|N-si；其中，s|i+2|N为所述信号序列S的第|i+2|N个元素；
模块204求取DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure BDA0002582314070000033
的DFT变换系数矢量为
Figure BDA0002582314070000034
求取公式为
Figure BDA0002582314070000035
第n个所述信号二阶差分序列
Figure BDA0002582314070000036
的DFT变换系数矢量为
Figure BDA0002582314070000037
求取公式为
Figure BDA0002582314070000038
模块205求取归一化DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure BDA0002582314070000039
的DFT变换系数矢量
Figure BDA00025823140700000310
的归一化DFT变换系数矢量为
Figure BDA00025823140700000311
其计算公式为
Figure BDA00025823140700000312
第n个所述信号二阶差分序列
Figure BDA00025823140700000313
的DFT变换系数矢量
Figure BDA00025823140700000314
的归一化DFT变换系数矢量为
Figure BDA00025823140700000315
其计算公式为
Figure BDA00025823140700000316
模块206求取N个频率差异量，具体为：第n个频率差异量记为Hn，其计算公式为
Figure BDA00025823140700000317
其中，
Figure BDA00025823140700000318
为所述信号一阶差分序列
Figure BDA00025823140700000319
的中心频率；
Figure BDA00025823140700000320
为所述信号二阶差分序列
Figure BDA00025823140700000321
的中心频率；D0为J差异因子，其计算公式为
Figure BDA00025823140700000322
其中，In是维数为n的全1列矢量；Δf为所述信号序列的采样频率；
模块207求取事件判断阈值，具体为：事件判断阈值记为ε0，其计算公式为
Figure BDA00025823140700000323
Figure BDA00025823140700000324
模块208判断负荷开关事件，具体为：如果第n个频率差异量Hn大于或者等于所述事件判断阈值ε0，则在所述信号序列的第n点处检测到负荷开关事件；否则，未检测到负荷开关事件。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
现在常用的负载开关事件往往利用功率数据的变化来确定：当功率变化值超过预置的阈值时，则认为有负载开关事件发生。这种方法虽然简单易行，但由于脉冲噪声以及非线性负载的普遍使用，造成开关事件检测精度的极大下降。
本发明的目的是提供一种利用频率差异量的负荷开关事件检测方法和***，所提出的方法利用了不同负荷开关动作时所产生的暂态功率信号差异，通过频率差异量性质区分不同负荷的不同开关事件以及异常功率信号造成的异常事件。所提出的方法具有较好的开关事件检测性能，计算也非常简单。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍。显而易见，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的方法流程示意图；
图2为本发明的***流程示意图；
图3为本发明的具体实施案例流程示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1一种利用频率差异量的负荷开关事件检测方法的流程示意图
图1为本发明一种利用频率差异量的负荷开关事件检测方法的流程示意图。如图1所示，所述的一种利用频率差异量的负荷开关事件检测方法具体包括以下步骤：
步骤101获取按时间顺序采集的信号序列S；
步骤102生成N个信号一阶差分序列，具体为：第n个信号一阶差分序列记为
Figure BDA0002582314070000041
其第i个元素的值为s|i+1|N-si；其中，i为元素信号，其取值范围为i＝1,···,N；n为差分序列序号，其取值范围为n＝1,2,···,N；N为所述信号序列S的长度；si为所述信号序列S的第i个元素；s|i+1|N为所述信号序列S的第|i+1|N个元素；|i+1|N表示以N为模取余数运算；
步骤103生成N个信号二阶差分序列，具体为：第n个信号二阶差分序列记为
Figure BDA0002582314070000042
其第i个元素的值为s|i+2|N-si；其中，s|i+2|N为所述信号序列S的第|i+2|N个元素；
步骤104求取DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure BDA0002582314070000043
的DFT变换系数矢量为
Figure BDA0002582314070000044
求取公式为
Figure BDA0002582314070000045
第n个所述信号二阶差分序列
Figure BDA0002582314070000046
的DFT变换系数矢量为
Figure BDA0002582314070000047
求取公式为
Figure BDA0002582314070000048
步骤105求取归一化DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure BDA0002582314070000049
的DFT变换系数矢量
Figure BDA00025823140700000410
的归一化DFT变换系数矢量为
Figure BDA00025823140700000411
其计算公式为
Figure BDA00025823140700000412
第n个所述信号二阶差分序列
Figure BDA00025823140700000413
的DFT变换系数矢量
Figure BDA00025823140700000414
的归一化DFT变换系数矢量为
Figure BDA00025823140700000415
其计算公式为
Figure BDA00025823140700000416
步骤106求取N个频率差异量，具体为：第n个频率差异量记为Hn，其计算公式为
Figure BDA0002582314070000051
其中，
Figure BDA0002582314070000052
为所述信号一阶差分序列
Figure BDA0002582314070000053
的中心频率；
Figure BDA0002582314070000054
为所述信号二阶差分序列
Figure BDA0002582314070000055
的中心频率；D0为J差异因子，其计算公式为
Figure BDA0002582314070000056
其中，In是维数为n的全1列矢量；Δf为所述信号序列的采样频率；
步骤107求取事件判断阈值，具体为：事件判断阈值记为ε0，其计算公式为
Figure BDA0002582314070000057
Figure BDA0002582314070000058
步骤108判断负荷开关事件，具体为：如果第n个频率差异量Hn大于或者等于所述事件判断阈值ε0，则在所述信号序列的第n点处检测到负荷开关事件；否则，未检测到负荷开关事件。
图2一种利用频率差异量的负荷开关事件检测***的结构意图
图2为本发明一种利用频率差异量的负荷开关事件检测***的结构示意图。如图2所示，所述一种利用频率差异量的负荷开关事件检测***包括以下结构：
模块201获取按时间顺序采集的信号序列S；
模块202生成N个信号一阶差分序列，具体为：第n个信号一阶差分序列记为
Figure BDA0002582314070000059
其第i个元素的值为s|i+1|N-si；其中，i为元素信号，其取值范围为i＝1,···,N；n为差分序列序号，其取值范围为n＝1,2,···,N；N为所述信号序列S的长度；si为所述信号序列S的第i个元素；s|i+1|N为所述信号序列S的第|i+1|N个元素；|i+1|N表示以N为模取余数运算；
模块203生成N个信号二阶差分序列，具体为：第n个信号二阶差分序列记为
Figure BDA00025823140700000510
其第i个元素的值为s|i+2|N-si；其中，s|i+2|N为所述信号序列S的第|i+2|N个元素；
模块204求取DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure BDA00025823140700000511
的DFT变换系数矢量为
Figure BDA00025823140700000512
求取公式为
Figure BDA00025823140700000513
第n个所述信号二阶差分序列
Figure BDA00025823140700000514
的DFT变换系数矢量为
Figure BDA00025823140700000515
求取公式为
Figure BDA00025823140700000516
模块205求取归一化DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure BDA00025823140700000517
的DFT变换系数矢量
Figure BDA00025823140700000518
的归一化DFT变换系数矢量为
Figure BDA00025823140700000519
其计算公式为
Figure BDA00025823140700000520
第n个所述信号二阶差分序列
Figure BDA00025823140700000521
的DFT变换系数矢量
Figure BDA00025823140700000522
的归一化DFT变换系数矢量为
Figure BDA00025823140700000523
其计算公式为
Figure BDA00025823140700000524
模块206求取N个频率差异量，具体为：第n个频率差异量记为Hn，其计算公式为
Figure BDA00025823140700000525
其中，
Figure BDA00025823140700000526
为所述信号一阶差分序列
Figure BDA00025823140700000527
的中心频率；
Figure BDA00025823140700000528
为所述信号二阶差分序列
Figure BDA00025823140700000529
的中心频率；D0为J差异因子，其计算公式为
Figure BDA00025823140700000530
其中，In是维数为n的全1列矢量；Δf为所述信号序列的采样频率；
模块207求取事件判断阈值，具体为：事件判断阈值记为ε0，其计算公式为
Figure BDA0002582314070000061
Figure BDA0002582314070000062
模块208判断负荷开关事件，具体为：如果第n个频率差异量Hn大于或者等于所述事件判断阈值ε0，则在所述信号序列的第n点处检测到负荷开关事件；否则，未检测到负荷开关事件。
下面提供一个具体实施案例，进一步说明本发明的方案
图3为本发明具体实施案例的流程示意图。如图3所示，具体包括以下步骤：
步骤301获取按时间顺序采集的信号序列S；
步骤302生成N个信号一阶差分序列，具体为：第n个信号一阶差分序列记为
Figure BDA0002582314070000063
其第i个元素的值为s|i+1|N-si；其中，i为元素信号，其取值范围为i＝1,···,N；n为差分序列序号，其取值范围为n＝1,2,···,N；N为所述信号序列S的长度；si为所述信号序列S的第i个元素；s|i+1|N为所述信号序列S的第|i+1|N个元素；|i+1|N表示以N为模取余数运算；
步骤303生成N个信号二阶差分序列，具体为：第n个信号二阶差分序列记为
Figure BDA0002582314070000064
其第i个元素的值为s|i+2|N-si；其中，s|i+2|N为所述信号序列S的第|i+2|N个元素；
步骤304求取DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure BDA0002582314070000065
的DFT变换系数矢量为
Figure BDA0002582314070000066
求取公式为
Figure BDA0002582314070000067
第n个所述信号二阶差分序列
Figure BDA0002582314070000068
的DFT变换系数矢量为
Figure BDA0002582314070000069
求取公式为
Figure BDA00025823140700000610
步骤305求取归一化DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure BDA00025823140700000611
的DFT变换系数矢量
Figure BDA00025823140700000612
的归一化DFT变换系数矢量为
Figure BDA00025823140700000613
其计算公式为
Figure BDA00025823140700000614
第n个所述信号二阶差分序列
Figure BDA00025823140700000615
的DFT变换系数矢量
Figure BDA00025823140700000616
的归一化DFT变换系数矢量为
Figure BDA00025823140700000617
其计算公式为
Figure BDA00025823140700000618
步骤306求取N个频率差异量，具体为：第n个频率差异量记为Hn，其计算公式为
Figure BDA00025823140700000619
其中，
Figure BDA00025823140700000620
为所述信号一阶差分序列
Figure BDA00025823140700000621
的中心频率；
Figure BDA00025823140700000622
为所述信号二阶差分序列
Figure BDA00025823140700000623
的中心频率；D0为J差异因子，其计算公式为
Figure BDA00025823140700000624
其中，In是维数为n的全1列矢量；Δf为所述信号序列的采样频率；
步骤307求取事件判断阈值，具体为：事件判断阈值记为ε0，其计算公式为
Figure BDA00025823140700000625
Figure BDA00025823140700000626
步骤308判断负荷开关事件，具体为：如果第n个频率差异量Hn大于或者等于所述事件判断阈值ε0，则在所述信号序列的第n点处检测到负荷开关事件；否则，未检测到负荷开关事件。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言，由于其与实施例公开的方法相对应，所以描述较为简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (2)

1.一种利用频率差异量的负荷开关事件检测方法，其特征在于，包括：
步骤101获取按时间顺序采集的信号序列S；
步骤102生成N个信号一阶差分序列，具体为：第n个信号一阶差分序列记为
Figure FDA0002582314060000011
其第i个元素的值为
Figure FDA0002582314060000012
其中，i为元素信号，其取值范围为i＝1,···,N；n为差分序列序号，其取值范围为n＝1,2,···,N；N为所述信号序列S的长度；si为所述信号序列S的第i个元素；
Figure FDA0002582314060000013
为所述信号序列S的第|i+1|N个元素；|i+1|N表示以N为模取余数运算；
步骤103生成N个信号二阶差分序列，具体为：第n个信号二阶差分序列记为
Figure FDA0002582314060000014
其第i个元素的值为
Figure FDA0002582314060000015
其中，
Figure FDA0002582314060000016
为所述信号序列S的第|i+2|N个元素；
步骤104求取DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure FDA0002582314060000017
的DFT变换系数矢量为
Figure FDA0002582314060000018
求取公式为
Figure FDA0002582314060000019
第n个所述信号二阶差分序列
Figure FDA00025823140600000110
的DFT变换系数矢量为
Figure FDA00025823140600000111
求取公式为
Figure FDA00025823140600000112
步骤105求取归一化DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure FDA00025823140600000113
的DFT变换系数矢量
Figure FDA00025823140600000114
的归一化DFT变换系数矢量为
Figure FDA00025823140600000115
其计算公式为
Figure FDA00025823140600000116
第n个所述信号二阶差分序列
Figure FDA00025823140600000117
的DFT变换系数矢量
Figure FDA00025823140600000118
的归一化DFT变换系数矢量为
Figure FDA00025823140600000119
其计算公式为
Figure FDA00025823140600000120
步骤106求取N个频率差异量，具体为：第n个频率差异量记为Hn，其计算公式为
Figure FDA00025823140600000121
其中，
Figure FDA00025823140600000122
为所述信号一阶差分序列
Figure FDA00025823140600000123
的中心频率；
Figure FDA00025823140600000124
为所述信号二阶差分序列
Figure FDA00025823140600000125
的中心频率；D0为J差异因子，其计算公式为
Figure FDA00025823140600000126
其中，In是维数为n的全1列矢量；Δf为所述信号序列的采样频率；
步骤107求取事件判断阈值，具体为：事件判断阈值记为ε0，其计算公式为
Figure FDA00025823140600000127
步骤108判断负荷开关事件，具体为：如果第n个频率差异量Hn大于或者等于所述事件判断阈值ε0，则在所述信号序列的第n点处检测到负荷开关事件；否则，未检测到负荷开关事件。
2.一种利用频率差异量的负荷开关事件检测***，其特征在于，包括：
模块201获取按时间顺序采集的信号序列S；
模块202生成N个信号一阶差分序列，具体为：第n个信号一阶差分序列记为
Figure FDA00025823140600000128
其第i个元素的值为
Figure FDA00025823140600000129
其中，i为元素信号，其取值范围为i＝1,···,N；n为差分序列序号，其取值范围为n＝1,2,···,N；N为所述信号序列S的长度；si为所述信号序列S的第i个元素；
Figure FDA0002582314060000021
为所述信号序列S的第|i+1|N个元素；|i+1|N表示以N为模取余数运算；
模块203生成N个信号二阶差分序列，具体为：第n个信号二阶差分序列记为
Figure FDA0002582314060000022
其第i个元素的值为
Figure FDA0002582314060000023
其中，
Figure FDA0002582314060000024
为所述信号序列S的第|i+2|N个元素；
模块204求取DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure FDA0002582314060000025
的DFT变换系数矢量为
Figure FDA0002582314060000026
求取公式为
Figure FDA0002582314060000027
第n个所述信号二阶差分序列
Figure FDA0002582314060000028
的DFT变换系数矢量为
Figure FDA0002582314060000029
求取公式为
Figure FDA00025823140600000210
模块205求取归一化DFT变换系数矢量，具体为：第n个所述信号一阶差分序列
Figure FDA00025823140600000211
的DFT变换系数矢量
Figure FDA00025823140600000212
的归一化DFT变换系数矢量为
Figure FDA00025823140600000213
其计算公式为
Figure FDA00025823140600000214
第n个所述信号二阶差分序列
Figure FDA00025823140600000215
的DFT变换系数矢量
Figure FDA00025823140600000216
的归一化DFT变换系数矢量为
Figure FDA00025823140600000217
其计算公式为
Figure FDA00025823140600000218
模块206求取N个频率差异量，具体为：第n个频率差异量记为Hn，其计算公式为
Figure FDA00025823140600000219
其中，
Figure FDA00025823140600000220
为所述信号一阶差分序列
Figure FDA00025823140600000221
的中心频率；
Figure FDA00025823140600000222
为所述信号二阶差分序列
Figure FDA00025823140600000223
的中心频率；D0为J差异因子，其计算公式为
Figure FDA00025823140600000224
其中，In是维数为n的全1列矢量；Δf为所述信号序列的采样频率；
模块207求取事件判断阈值，具体为：事件判断阈值记为ε0，其计算公式为
Figure FDA00025823140600000225
模块208判断负荷开关事件，具体为：如果第n个频率差异量Hn大于或者等于所述事件判断阈值ε0，则在所述信号序列的第n点处检测到负荷开关事件；否则，未检测到负荷开关事件。
CN202010671056.5A 2020-07-14 2020-07-14 一种利用频率差异量的负荷开关事件检测方法和*** Withdrawn CN111830405A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010671056.5A CN111830405A (zh) 2020-07-14 2020-07-14 一种利用频率差异量的负荷开关事件检测方法和***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010671056.5A CN111830405A (zh) 2020-07-14 2020-07-14 一种利用频率差异量的负荷开关事件检测方法和***

Publications (1)

Publication Number Publication Date
CN111830405A true CN111830405A (zh) 2020-10-27

Family

ID=72922772

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010671056.5A Withdrawn CN111830405A (zh) 2020-07-14 2020-07-14 一种利用频率差异量的负荷开关事件检测方法和***

Country Status (1)

Country Link
CN (1) CN111830405A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113721103A (zh) * 2021-11-01 2021-11-30 广东电网有限责任公司东莞供电局 一种高压开关柜在线故障诊断方法及***

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113721103A (zh) * 2021-11-01 2021-11-30 广东电网有限责任公司东莞供电局 一种高压开关柜在线故障诊断方法及***
CN113721103B (zh) * 2021-11-01 2022-02-11 广东电网有限责任公司东莞供电局 一种高压开关柜在线故障诊断方法及***

Similar Documents

Publication Publication Date Title
CN109145825B (zh) 一种相干噪声滤除方法及***
CN111680590A (zh) 一种利用收缩梯度的功率信号滤波方法和***
CN111666870A (zh) 一种利用二次约束的功率信号重构方法和***
CN111830405A (zh) 一种利用频率差异量的负荷开关事件检测方法和***
CN109241874B (zh) 能量分解中功率信号滤波方法
CN110221119B (zh) 基于功率与akie融合信息的负载开关事件检测方法和***
CN112434567B (zh) 一种利用噪声抖动性质的功率信号滤波方法和***
CN108918929B (zh) 一种负荷分解中功率信号自适应滤波方法
CN110542855B (zh) 基于离散余弦变换的负荷开关事件检测方法和***
CN110196354B (zh) 一种负载的开关事件的检测方法和装置
CN112329637B (zh) 一种利用模式特征的负荷开关事件检测方法和***
CN110244115B (zh) 一种基于信号连接性的负载开关事件检测方法及***
CN111832474A (zh) 一种利用能量尺度的功率信号滤波方法和***
CN111639606A (zh) 一种利用Dantzig总梯度最小化的功率信号滤波方法和***
CN112307986B (zh) 一种利用高斯梯度的负荷开关事件检测方法和***
CN112257576B (zh) 一种利用Maha距离测度的负荷开关事件检测方法和***
CN110702981A (zh) 一种利用分类树的负荷开关事件检测方法和***
CN110749841A (zh) 一种利用转换空间因子的负荷开关事件检测方法和***
CN111948477A (zh) 一种利用固定b抽样的负荷开关事件检测方法和***
CN111737645A (zh) 一种利用预测矩阵的功率信号重构方法和***
CN112180154A (zh) 一种利用置信度优化的负荷开关事件检测方法和***
CN112180155A (zh) 一种利用紧支集的负荷开关事件检测方法和***
CN112180153A (zh) 一种利用KUllback-Leibler距离的负荷开关事件检测方法和***
CN112347922B (zh) 一种利用Hankerl矩阵的功率信号滤波方法和***
CN112270282B (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
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20201027