CN113704698B - 一种非介入式负荷辨识的事件检测方法及系统 - Google Patents

一种非介入式负荷辨识的事件检测方法及系统 Download PDF

Info

Publication number
CN113704698B
CN113704698B CN202111265840.7A CN202111265840A CN113704698B CN 113704698 B CN113704698 B CN 113704698B CN 202111265840 A CN202111265840 A CN 202111265840A CN 113704698 B CN113704698 B CN 113704698B
Authority
CN
China
Prior art keywords
window
active power
point
fluctuation
time point
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
Application number
CN202111265840.7A
Other languages
English (en)
Other versions
CN113704698A (zh
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.)
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
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 China Electric Power Research Institute Co Ltd CEPRI filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN202111265840.7A priority Critical patent/CN113704698B/zh
Publication of CN113704698A publication Critical patent/CN113704698A/zh
Application granted granted Critical
Publication of CN113704698B publication Critical patent/CN113704698B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • 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
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Mathematical Optimization (AREA)
  • Operations Research (AREA)
  • Mathematical Analysis (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Development Economics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Quality & Reliability (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Probability & Statistics with Applications (AREA)
  • Game Theory and Decision Science (AREA)
  • Algebra (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Remote Monitoring And Control Of Power-Distribution Networks (AREA)

Abstract

本申请公开了一种非介入式负荷辨识的事件检测方法及系统。其中，该方法包括：基于实时采集的计量点电压电流波形计算实时的有功功率波形，并获取所述有功功率波形的波动点以及所述波点的波动时间点集合；基于所述波动点的波动时间点集合，计算基态窗内的有功功率均值，并根据所述基态窗内的有功功率均值，确定检测窗内的有功功率累加值，当累加和大于判定阈值时，确定当前时间点为启动/停机暂态开启点；当暂态为负荷状态时，根据所述窗宽，计算检测窗与稳态窗内的差值累加和以及稳态窗方差，当所述差值累加和大于所述判定阈值，并且所述稳态窗方差小于预定参数与稳态窗内的有功功率稳态值的乘积时，确定当前时间点为启动/停机暂态结束点。

Description

一种非介入式负荷辨识的事件检测方法及系统
技术领域
本申请涉及电力负荷用电监测特别是涉及一种非介入式负荷辨识的事件检测方法及系统。
背景技术
工商业用户和居民用户的用电负荷相比，其用电设备种类繁多，且启动方式相差较大，稳态运行时波动较大，启动的暂态时间差异也比较明显，适用于居民用户的非介入式负荷辨识事件检测方案并不完全适用于工商业用户，容易出现漏判和误判现象。此外，居民负荷的启停阶段的暂态过程通常在5s以内，工商业负荷的暂态过程多在1min以上，需要准确判断设备开启和关闭状态中暂态开始和暂态结束四类状态。
现有工商业用户的负荷的事件检测方法，鲁棒性不强、不可以支撑长暂态过程、不能准确判断暂态开始与结束。
发明内容
本公开的实施例提供了一种非介入式负荷辨识的事件检测方法及系统，以至少解决现有技术中存在的现有工商业用户的负荷的事件检测方法，鲁棒性不强、不可以支撑长暂态过程、不能准确判断暂态开始与结束的技术问题。
根据本公开实施例的一个方面，提供了一种非介入式负荷辨识的事件检测方法，包括：基于实时采集的计量点电压电流波形计算实时的有功功率波形，并获取所述有功功率波形的波动点以及所述波动点的波动时间点集合；
基于所述波动点的波动时间点集合，确定波动时间点之间的时间差，去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值，并将所述波动点均值作为窗宽；
当暂态为开始状态时，根据所述窗宽，计算基态窗内的有功功率均值，并根据所述基态窗内的有功功率均值，确定检测窗内的有功功率累加值，当累加和大于判定阈值时，确定当前时间点为启动/停机暂态开启点；
当暂态为负荷状态时，根据所述窗宽，计算检测窗与稳态窗内的差值累加和以及稳态窗方差，当所述差值累加和大于所述判定阈值，并且所述稳态窗方差小于预定参数与稳态窗内的有功功率稳态值的乘积时，确定当前时间点为启动/停机暂态结束点。
根据本公开实施例的另一方面，还提供了一种非介入式负荷辨识的事件检测系统，包括：
获取波动点模块，用于基于实时采集的计量点电压电流波形计算实时的有功功率波形，并获取所述有功功率波形的波动点以及所述波动点的波动时间点集合；
确定窗宽模块，用于基于所述波动点的波动时间点集合，确定波动时间点之间的时间差，去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值，并将所述波动点均值作为窗宽；
判定开始状态模块，用于当暂态为开始状态时，根据所述窗宽，计算基态窗内的有功功率均值，并根据所述基态窗内的有功功率均值，确定检测窗内的有功功率累加值，当累加和大于判定阈值时，确定当前时间点为启动/停机暂态开启点；
判定负荷暂态模块，用于当暂态为负荷状态时，根据所述窗宽，计算检测窗与稳态窗内的差值累加和以及稳态窗方差，当所述差值累加和大于所述判定阈值，并且所述稳态窗方差小于预定参数与稳态窗内的有功功率稳态值的乘积时，确定当前时间点为启动/停机暂态结束点。
在本发明中，通过新增稳态窗，首次提出开启暂态结束点和关闭暂态结束点的判断方法。在波动点的基础上进行状态判断，无需实时计算，资源占用率较低。改进偏移累加量判定方法，降低对于负荷高频波动情况下的误判率。
附图说明
此处所说明的附图用来提供对本公开的进一步理解，构成本申请的一部分，本公开的示意性实施例及其说明用于解释本公开，并不构成对本公开的不当限定。在附图中：
图1是根据本公开实施例所述的一种非介入式负荷辨识的事件检测方法的流程示意图；
图2是根据本公开实施例所述的一种非介入式负荷辨识的事件检测系统的示意图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式，然而，本发明可以用许多不同的形式来实施，并且不局限于此处描述的实施例，提供这些实施例是为了详尽地且完全地公开本发明，并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中，相同的单元/元件使用相同的附图标记。
除非另有说明，此处使用的术语（包括科技术语）对所属技术领域的技术人员具有通常的理解含义。另外，可以理解的是，以通常使用的词典限定的术语，应当被理解为与其相关领域的语境具有一致的含义，而不应该被理解为理想化的或过于正式的意义。
根据本实施例的第一个方面，提供了一种非介入式负荷辨识的事件检测方法100。参考图1所示，该方法100包括：
S102: 基于实时采集的计量点电压电流波形计算实时的有功功率波形，并获取所述有功功率波形的波动点以及所述波动点的波动时间点集合；
S104: 基于所述波动点的波动时间点集合，确定波动时间点之间的时间差，去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值，并将所述波动点均值作为窗宽；
S106: 当暂态为开始状态时，根据所述窗宽，计算基态窗内的有功功率均值，并根据所述基态窗内的有功功率均值，确定检测窗内的有功功率累加值，当累加和大于判定阈值时，确定当前时间点为启动/停机暂态开启点；
S108: 当暂态为负荷状态时，根据所述窗宽，计算检测窗与稳态窗内的差值累加和以及稳态窗方差，当所述差值累加和大于所述判定阈值，并且所述稳态窗方差小于预定参数与稳态窗内的有功功率稳态值的乘积时，确定当前时间点为启动/停机暂态结束点。
具体地，本实施例首先基于实时采集的计量点电压电流波形计算实时的有功功率波形，应用自适应抽样方法获取波形的波动点，设定基态窗、监测窗和稳态窗，比较基态窗和监测窗之间的差值，确认启动/停机暂态开始点，比较监测窗和稳态窗之间的差值，确认启动/停机暂态结束点。其中暂态开始点与暂态结束点成一一对应的关系，完成暂态开始点的判断后再进入暂态结束点的判断。
具体步骤如下：
步骤1：基于实时采集的计量点电压电流波形计算实时的有功功率波形，并应用自适应抽样方法获取有功功率波形的波动点；
步骤2：基于波动点序列，计算波动时间点之间的时间差，去除最大最最小的干扰项，计算剩余差值的均值
Figure 727015DEST_PATH_IMAGE001
，并将其作为基态窗、监测窗和稳态窗的窗宽；
步骤3：进入暂态开始状态的时间点判断流程，设定基态窗和检测窗，计算基态窗内的有功功率均值，改进CUSUM方法，计算检测窗与基态窗内的差值累加和，当累加和大于判定阈值时，当前时间点为启动/停机暂态开启点；
步骤4：当判定出现负荷暂态开始时间点，进入暂态结束时间点判断流程，当完成暂态结束时间点判断后，才能进入下一次的负荷暂态开始时间点判断，设定检测窗和稳态窗，计算稳态窗内的有功功率均值和方差，计算检测窗与稳态窗内的差值累加和，当累加和稳态窗内有功功率方差均满足判定条件时，当前时间点为启动/停机暂态结束点。
进一步，所述步骤1包括以下的步骤：
步骤1.1：计算实时的有功功率时序序列
Figure 3276DEST_PATH_IMAGE002
：
Figure 971232DEST_PATH_IMAGE003
（1）
其中
Figure 332943DEST_PATH_IMAGE004
和
Figure 310127DEST_PATH_IMAGE005
为瞬时电压电流采样点值，
Figure 440894DEST_PATH_IMAGE006
为工频周波的采样点数，根据具体需求，可将参与计算有功功率的实时电压电流采样点数扩大至
Figure 845330DEST_PATH_IMAGE007
的整数倍；
步骤1.2：对有功功率波形进行中值滤波，设定滤波窗宽为q，保留波动上升沿与下降沿，消除噪声干扰；
步骤1.3：初始化异动点判定条件，判定当前点是否属于临界点和异动点，并根据当前时间点的有功功率值实时更新双重判定阈值；
步骤1.4：完成总线侧多维特征量序列的异动点判定，保留序列中的异动点，完成自适应抽样，在保留原始波形轮廓信息的基础上，完成波动信号的有效抽样；
步骤1.5：获取有功功率时序序列中的波动时间点集合
Figure 428758DEST_PATH_IMAGE008
；
进一步，所述步骤2包括以下的步骤：
步骤2.1：计算波动时间点之间的时间差，并进行排序，得到差值集合
Figure 944053DEST_PATH_IMAGE009
（2）
Figure 444174DEST_PATH_IMAGE010
（3）
步骤2.2：去除最大最小的干扰项，计算剩余差值的均值
Figure 19511DEST_PATH_IMAGE011
，并将其作为基态窗、监测窗和稳态窗的窗宽：
Figure 355815DEST_PATH_IMAGE012
（4）
其中a为设定的判定阈值。
进一步，所述步骤3包括以下的步骤：
步骤3.1：对于波动时间点
Figure 409221DEST_PATH_IMAGE013
，基态窗和检测窗内的有功功率序列为：
Figure 514581DEST_PATH_IMAGE014
（5）
Figure 260820DEST_PATH_IMAGE015
（6）
步骤3.2：计算基态窗内的有功功率基准值：
Figure 818840DEST_PATH_IMAGE016
（7）
步骤3.3：设定参数
Figure 941517DEST_PATH_IMAGE017
，向后遍历检测窗内的有功功率值，基于传统的CUSUM方法进行改进，计算累加和：
Figure 635803DEST_PATH_IMAGE018
（8）
Figure 552944DEST_PATH_IMAGE019
（9）
步骤3.4：设定判定阈值h，当累加和
Figure 614572DEST_PATH_IMAGE020
满足以下条件时，
Figure 275360DEST_PATH_IMAGE020
所对应的时间点即为启动/停机暂态开启点：
Figure 824153DEST_PATH_IMAGE021
（10）
步骤3.5：当检测窗内完成遍历，累加和
Figure 177774DEST_PATH_IMAGE020
未满足式（10）时，基态窗和检测窗移至离
Figure 710387DEST_PATH_IMAGE022
最近的波动点进行下一轮判定；
进一步，所述步骤4包括以下的步骤：
步骤4.1：当检测到启动/停机暂态开始点后，检测窗和稳态窗移至最近的波动点进行启动/停机暂态结束点的判定。对于波动时间点
Figure 174866DEST_PATH_IMAGE023
，检测窗和稳态窗内的有功功率序列为：
Figure 578166DEST_PATH_IMAGE024
（11）
Figure DEST_PATH_IMAGE025
（12）
步骤4.2：计算稳态窗内的有功功率稳态值和方差：
Figure 368267DEST_PATH_IMAGE026
（13）
Figure 388176DEST_PATH_IMAGE027
（14）
步骤4.3：设定参数
Figure 125188DEST_PATH_IMAGE028
，向前遍历检测窗内的有功功率值，计算累加和：
Figure 490568DEST_PATH_IMAGE029
（15）
Figure 185992DEST_PATH_IMAGE030
（16）
步骤4.4：设定参数
Figure 427617DEST_PATH_IMAGE031
，当累加和
Figure 233899DEST_PATH_IMAGE032
和稳态窗方差S满足以下条件时，
Figure 877370DEST_PATH_IMAGE032
所对应的时间点即为启动/停机暂态结束点：
Figure 478116DEST_PATH_IMAGE033
（17）
Figure 472616DEST_PATH_IMAGE034
（18）
步骤4.5：当检测窗内完成遍历，累加和和稳态窗方差S满足未满足式（17）和式（18）时，检测窗和稳态窗移至离
Figure 817010DEST_PATH_IMAGE035
最近的波动点进行下一轮判定。
可选地，基于实时采集的计量点电压电流波形计算实时的有功功率波形，并获取所述有功功率波形的波动点以及所述波动点的波动时间点集合，包括：
计算实时的有功功率时序序列
Figure 314988DEST_PATH_IMAGE036
：
Figure 86635DEST_PATH_IMAGE037
（1）
其中
Figure 319164DEST_PATH_IMAGE038
和
Figure 467248DEST_PATH_IMAGE039
为瞬时电压电流采样点值，
Figure 819732DEST_PATH_IMAGE040
为工频周波的采样点数；
对所述有功功率时序序列的有功功率波形进行中值滤波，设定滤波窗宽为q，保留波动上升沿与下降沿，根据所述上升沿与所述下降沿，消除噪声干扰；
获取有功功率时序序列中的波动时间点集合
Figure 27860DEST_PATH_IMAGE041
，其中，n为波动时间点的序列号。
可选地，基于所述波动点的波动时间点集合，确定波动时间点之间的时间差，去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值，并将所述波动点均值作为窗宽，包括：
基于所述波动点的波动点序列，确定波动时间点之间的时间差，并对所述时间差进行排序，得到差值集合
Figure 731374DEST_PATH_IMAGE042
Figure 417570DEST_PATH_IMAGE043
（2）
Figure 890140DEST_PATH_IMAGE044
（3）
其中，n为波动时间点的序列号；
去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值
Figure 3589DEST_PATH_IMAGE045
，并将所述波动点均值作为基态窗、监测窗和稳态窗的窗宽：
Figure 194399DEST_PATH_IMAGE046
（4）
其中a为设定的判定阈值。
可选地，当暂态为开始状态时，根据所述窗宽，计算基态窗内的有功功率均值，并根据所述基态窗内的有功功率均值，确定检测窗内的有功功率累加值，当累加和大于判定阈值时，确定当前时间点为启动/停机暂态开启点，包括：
对于第一波动时间点
Figure 684286DEST_PATH_IMAGE047
，确定所述第一波动时间点的基态窗内的有功功率和所述第一波动时间点的检测窗内的有功功率为：
Figure 526209DEST_PATH_IMAGE048
（5）
Figure 810560DEST_PATH_IMAGE049
（6）
其中，
Figure 488666DEST_PATH_IMAGE050
为所述第一波动时间点的基态窗内的有功功率
Figure 516665DEST_PATH_IMAGE051
为基态窗以及检测窗的窗宽，
Figure 963827DEST_PATH_IMAGE052
为在所述第一波动时间点往前数窗宽
Figure 419079DEST_PATH_IMAGE053
的时间，
Figure 584481DEST_PATH_IMAGE054
为在所述第一波动时间点往前数窗宽
Figure 416171DEST_PATH_IMAGE055
再加1秒的时间，
Figure 717839DEST_PATH_IMAGE056
为在所述第一波动时间点往前数窗宽
Figure 609572DEST_PATH_IMAGE057
对应的有功功率，
Figure 996691DEST_PATH_IMAGE058
为与在所述第一波动时间点往前数窗宽
Figure 117224DEST_PATH_IMAGE053
再加1秒的时间对应的功率，
Figure 538979DEST_PATH_IMAGE059
为所述第一波动时间点的有功功率；
Figure 601612DEST_PATH_IMAGE060
为所述第一波动时间点的检测窗内的有功功率，
Figure 210448DEST_PATH_IMAGE061
为第一波动时间点加1秒的时间，
Figure 649520DEST_PATH_IMAGE062
为第一波动时间点加2秒的时间，
Figure 925780DEST_PATH_IMAGE063
为在所述第一波动时间点往后数窗宽
Figure 893736DEST_PATH_IMAGE055
的时间，
Figure 255448DEST_PATH_IMAGE064
为第一波动时间点加1秒的时间对应的有功功率，
Figure 967052DEST_PATH_IMAGE065
为为第一波动时间点加2秒的时间对应的有功功率，
Figure 363398DEST_PATH_IMAGE066
所述第一波动时间点往后数窗宽
Figure 20032DEST_PATH_IMAGE055
的时间对应的有功功率；
根据所述第一波动时间点的基态窗内的有功功率，计算所述第一波动时间点的基态窗内的有功功率基准值：
Figure 869039DEST_PATH_IMAGE067
（7）
其中，
Figure 118755DEST_PATH_IMAGE068
为基态窗内的有功功率基准值， Pk为基态窗内k时间点的瞬时有功功率值；
遍历所述检测窗内的有功功率值，基于CUSUM方法进行改进，计算累加所述检测窗内的有功功率值，确定检测窗内的有功功率累加值：
Figure 369608DEST_PATH_IMAGE069
（8）
Figure 944945DEST_PATH_IMAGE070
（9）
其中，
Figure 281249DEST_PATH_IMAGE071
为检测窗内的有功功率累加值的初始值，
Figure 334656DEST_PATH_IMAGE072
为检测窗内的有功功率累加值；
设定判定有功功率阈值h，当有功功率累加值
Figure 440015DEST_PATH_IMAGE072
满足以下条件时，确定
Figure 186254DEST_PATH_IMAGE072
所对应的时间点为启动/停机暂态开启点：
Figure 744274DEST_PATH_IMAGE073
（10）
当检测窗内完成遍历，并且所述有功功率累加值
Figure 617683DEST_PATH_IMAGE074
未满足式（10）时，基态窗和检测窗移至离
Figure 311970DEST_PATH_IMAGE075
最近的波动点进行下一轮判定。
可选地，当暂态为负荷状态时，根据所述窗宽，计算检测窗与稳态窗内的差值累加和以及稳态窗方差，当所述差值累加和大于所述判定阈值，并且所述稳态窗方差小于预定参数与稳态窗内的有功功率稳态值的乘积时，确定当前时间点为启动/停机暂态结束点，包括：
对于第二波动时间点
Figure 229110DEST_PATH_IMAGE076
，确定所述第二波动时间点的检测窗内的有功功率和所述第二波动时间点的稳态窗内的有功功率为：
Figure 540006DEST_PATH_IMAGE077
（11）
Figure 200794DEST_PATH_IMAGE025
（12）
其中，
Figure 15167DEST_PATH_IMAGE078
为所述第二波动时间点的检测窗内的有功功率，
Figure 103208DEST_PATH_IMAGE051
为检测窗以及稳态窗的窗宽，
Figure 635821DEST_PATH_IMAGE079
为在所述第二波动时间点往前数窗宽
Figure 834721DEST_PATH_IMAGE055
的时间，
Figure 503600DEST_PATH_IMAGE080
为在所述第二波动时间点往前数窗宽
Figure 28122DEST_PATH_IMAGE053
再加1秒的时间，
Figure 31719DEST_PATH_IMAGE081
为在所述第二波动时间点往前数窗宽
Figure 299889DEST_PATH_IMAGE057
对应的有功功率，
Figure 88854DEST_PATH_IMAGE082
为与在所述第二波动时间点往前数窗宽
Figure 253119DEST_PATH_IMAGE055
再加1秒的时间对应的功率，
Figure 760324DEST_PATH_IMAGE083
为所述第二波动时间点的有功功率；
Figure 566606DEST_PATH_IMAGE084
为所述第二波动时间点的稳态窗内的有功功率，
Figure 944497DEST_PATH_IMAGE085
为第二波动时间点加1秒的时间，
Figure 76401DEST_PATH_IMAGE086
为第二波动时间点加2秒的时间，
Figure 539744DEST_PATH_IMAGE087
为在所述第二波动时间点往后数窗宽
Figure 149717DEST_PATH_IMAGE055
的时间，
Figure 664006DEST_PATH_IMAGE088
为第二波动时间点加1秒的时间对应的有功功率，
Figure 435653DEST_PATH_IMAGE089
为第二波动时间点加2秒的时间对应的有功功率，
Figure 917450DEST_PATH_IMAGE090
所述第二波动时间点往后数窗宽
Figure 799955DEST_PATH_IMAGE055
的时间对应的有功功率；
根据所述第二波动时间点的稳态窗内的有功功率，计算稳态窗内的有功功率稳态值和稳态窗内的方差：
Figure 418018DEST_PATH_IMAGE091
（13）
Figure 360566DEST_PATH_IMAGE092
（14）
其中，
Figure 798501DEST_PATH_IMAGE093
为稳态窗内的有功功率稳态值，
Figure 750276DEST_PATH_IMAGE094
为稳态窗内的方差；
设定参数
Figure 222846DEST_PATH_IMAGE095
，向前遍历检测窗内的有功功率值，计算检测窗与稳态窗内的差值累加和：
Figure 601875DEST_PATH_IMAGE096
（15）
Figure 44882DEST_PATH_IMAGE097
（16）
设定参数
Figure 534769DEST_PATH_IMAGE098
，当所述差值累加和
Figure 861845DEST_PATH_IMAGE099
和稳态窗内的方差S满足以下条件时，确定所述差值累加和
Figure 411775DEST_PATH_IMAGE100
所对应的时间点为启动/停机暂态结束点：
Figure 824302DEST_PATH_IMAGE101
（17）
Figure 852301DEST_PATH_IMAGE102
（18）
当检测窗内完成遍历，并且所述差值累加
Figure 299463DEST_PATH_IMAGE103
和和稳态窗方差S满足未满足式（17）和式（18）时，检测窗和稳态窗移至离
Figure 20294DEST_PATH_IMAGE104
最近的波动点进行下一轮判定。
从而，通过新增稳态窗，首次提出开启暂态结束点和关闭暂态结束点的判断方法。在波动点的基础上进行状态判断，无需实时计算，资源占用率较低。改进偏移累加量判定方法，降低对于负荷高频波动情况下的误判率。
根据本申请的另一个方面，还提供了一种非介入式负荷辨识的事件检测系统200。参考图2所示，该系统200包括：
获取波动点模块210，用于基于实时采集的计量点电压电流波形计算实时的有功功率波形，并获取所述有功功率波形的波动点以及所述波动点 的波动时间点集合；
确定窗宽模块220，用于基于所述波动点的波动时间点集合，确定波动时间点之间的时间差，去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值，并将所述波动点均值作为窗宽；
判定开始状态模块230，用于当暂态为开始状态时，根据所述窗宽，计算基态窗内的有功功率均值，并根据所述基态窗内的有功功率均值，确定检测窗内的有功功率累加值，当累加和大于判定阈值时，确定当前时间点为启动/停机暂态开启点；
判定负荷暂态模块240，用于当暂态为负荷状态时，根据所述窗宽，计算检测窗与稳态窗内的差值累加和以及稳态窗方差，当所述差值累加和大于所述判定阈值，并且所述稳态窗方差小于预定参数与稳态窗内的有功功率稳态值的乘积时，确定当前时间点为启动/停机暂态结束点。
可选地，获取波动点模块210，包括：
计算有功功率时序子模块，用于计算实时的有功功率时序序列
Figure 185696DEST_PATH_IMAGE105
：
Figure 751807DEST_PATH_IMAGE106
（1）
其中
Figure 53475DEST_PATH_IMAGE107
和
Figure 695940DEST_PATH_IMAGE108
为瞬时电压电流采样点值，
Figure 83059DEST_PATH_IMAGE109
为工频周波的采样点数；
消除噪声干扰子模块，用于对所述有功功率时序序列的有功功率波形进行中值滤波，设定滤波窗宽为q，保留波动上升沿与下降沿，根据所述上升沿与所述下降沿，消除噪声干扰；
获取波动时间点集合子模块，用于获取有功功率时序序列中的波动时间点集合
Figure 452861DEST_PATH_IMAGE110
，其中，n为波动时间点的序列号。
可选地，确定窗宽模块220，包括：
得到差值集合子模块，用于基于所述波动点的波动点序列，确定波动时间点之间的时间差，并对所述时间差进行排序，得到差值集合
Figure 874615DEST_PATH_IMAGE111
Figure 937249DEST_PATH_IMAGE112
（2）
Figure 811664DEST_PATH_IMAGE113
（3）
其中，n为波动时间点的序列号
确定窗宽子模块，用于去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值
Figure 985156DEST_PATH_IMAGE114
，并将所述波动点均值作为基态窗、监测窗和稳态窗的窗宽：
Figure 995837DEST_PATH_IMAGE115
（4）
其中a为设定的判定阈值。
可选地，判定开始状态模块230，包括：
确定第一波动时间点有功功率子模块，用于对于第一波动时间点
Figure 229373DEST_PATH_IMAGE116
，确定所述第一波动时间点的基态窗内的有功功率和所述第一波动时间点的检测窗内的有功功率为：
Figure 591084DEST_PATH_IMAGE117
（5）
Figure 817535DEST_PATH_IMAGE118
（6）
其中，
Figure 948302DEST_PATH_IMAGE119
为所述第一波动时间点的基态窗内的有功功率，
Figure 352738DEST_PATH_IMAGE120
为基态窗以及检测窗的窗宽，
Figure 936166DEST_PATH_IMAGE052
为在所述第一波动时间点往前数窗宽
Figure 717041DEST_PATH_IMAGE053
的时间，
Figure 967893DEST_PATH_IMAGE121
为在所述第一波动时间点往前数窗宽
Figure 277652DEST_PATH_IMAGE055
再加1秒的时间，
Figure 348376DEST_PATH_IMAGE122
为在所述第一波动时间点往前数窗宽
Figure 667362DEST_PATH_IMAGE057
对应的有功功率，
Figure 772721DEST_PATH_IMAGE058
为与在所述第一波动时间点往前数窗宽
Figure 269693DEST_PATH_IMAGE053
再加1秒的时间对应的功率，
Figure 827713DEST_PATH_IMAGE059
为所述第一波动时间点的有功功率；
Figure 684811DEST_PATH_IMAGE123
为所述第一波动时间点的检测窗内的有功功率，
Figure 910256DEST_PATH_IMAGE124
为第一波动时间点加1秒的时间，
Figure 827396DEST_PATH_IMAGE062
为第一波动时间点加2秒的时间，
Figure 607133DEST_PATH_IMAGE125
为在所述第一波动时间点往后数窗宽
Figure 267922DEST_PATH_IMAGE055
的时间，
Figure 347873DEST_PATH_IMAGE064
为第一波动时间点加1秒的时间对应的有功功率，
Figure 435915DEST_PATH_IMAGE126
为第一波动时间点加2秒的时间对应的有功功率，
Figure 968527DEST_PATH_IMAGE066
所述第一波动时间点往后数窗宽
Figure 685204DEST_PATH_IMAGE053
的时间对应的有功功率；
计算有功功率基准值子模块，用于根据所述第一波动时间点的基态窗内的有功功率，计算所述第一波动时间点的基态窗内的有功功率基准值：
Figure 619662DEST_PATH_IMAGE127
（7）
其中，
Figure 613026DEST_PATH_IMAGE128
为基态窗内的有功功率基准值， Pk为基态窗内k时间点的瞬时有功功率值；
确定有功功率累加值子模块，用于遍历所述检测窗内的有功功率值，基于CUSUM方法进行改进，计算累加所述检测窗内的有功功率值，确定检测窗内的有功功率累加值：
Figure 632934DEST_PATH_IMAGE129
（8）
Figure 901105DEST_PATH_IMAGE130
（9）
其中，
Figure 690069DEST_PATH_IMAGE131
为检测窗内的有功功率累加值的初始值，
Figure 854334DEST_PATH_IMAGE132
为检测窗内的有功功率累加值；
确定暂态开启点子模块，用于设定判定有功功率阈值h，当有功功率累加值
Figure 361539DEST_PATH_IMAGE132
满足以下条件时，确定
Figure 167821DEST_PATH_IMAGE132
所对应的时间点为启动/停机暂态开启点：
Figure 811292DEST_PATH_IMAGE133
（10）
开始状态下一轮判定子模块，用于当检测窗内完成遍历，并且所述有功功率累加值
Figure 428349DEST_PATH_IMAGE132
未满足式（10）时，基态窗和检测窗移至离
Figure 891692DEST_PATH_IMAGE134
最近的波动点进行下一轮判定。
可选地，判定负荷暂态模块240，包括：
确定第二波动时间点有功功率子模块，用于对于第二波动时间点
Figure 501664DEST_PATH_IMAGE135
，确定所述第二波动时间点的检测窗内的有功功率和所述第二波动时间点的稳态窗内的有功功率为：
Figure 999642DEST_PATH_IMAGE136
（11）
Figure 771289DEST_PATH_IMAGE137
（12）
其中，
Figure 987506DEST_PATH_IMAGE138
为所述第二波动时间点的检测窗内的有功功率，
Figure 135591DEST_PATH_IMAGE120
为检测窗以及稳态窗的窗宽，
Figure 753654DEST_PATH_IMAGE079
为在所述第二波动时间点往前数窗宽
Figure 961782DEST_PATH_IMAGE053
的时间，
Figure 665295DEST_PATH_IMAGE080
为在所述第二波动时间点往前数窗宽
Figure 351492DEST_PATH_IMAGE055
再加1秒的时间，
Figure 73329DEST_PATH_IMAGE081
为在所述第二波动时间点往前数窗宽
Figure 452358DEST_PATH_IMAGE057
对应的有功功率，
Figure 643168DEST_PATH_IMAGE139
为与在所述第二波动时间点往前数窗宽
Figure 867476DEST_PATH_IMAGE053
再加1秒的时间对应的功率，
Figure 460131DEST_PATH_IMAGE083
为所述第二波动时间点的有功功率；
Figure 744482DEST_PATH_IMAGE140
为所述第二波动时间点的稳态窗内的有功功率，
Figure 422588DEST_PATH_IMAGE085
为第二波动时间点加1秒的时间，
Figure 716166DEST_PATH_IMAGE086
为第二波动时间点加2秒的时间，
Figure 897748DEST_PATH_IMAGE141
为在所述第二波动时间点往后数窗宽
Figure 618580DEST_PATH_IMAGE055
的时间，
Figure 269135DEST_PATH_IMAGE088
为第二波动时间点加1秒的时间对应的有功功率，
Figure 100825DEST_PATH_IMAGE142
为第二波动时间点加2秒的时间对应的有功功率，
Figure 668072DEST_PATH_IMAGE090
所述第二波动时间点往后数窗宽
Figure 294226DEST_PATH_IMAGE055
的时间对应的有功功率；
计算稳态窗方差子模块，用于根据所述第二波动时间点的稳态窗内的有功功率，计算稳态窗内的有功功率稳态值和稳态窗内的方差：
Figure 681345DEST_PATH_IMAGE143
（13）
Figure 316725DEST_PATH_IMAGE144
（14）
其中，
Figure 738480DEST_PATH_IMAGE145
为稳态窗内的有功功率稳态值，
Figure 535534DEST_PATH_IMAGE146
为稳态窗内的方差；
计算差值累加和子模块，用于设定参数
Figure 409949DEST_PATH_IMAGE147
，向前遍历检测窗内的有功功率值，计算检测窗与稳态窗内的差值累加和：
Figure 583442DEST_PATH_IMAGE148
（15）
Figure 100181DEST_PATH_IMAGE149
（16）
确定暂态结束点子模块，用于设定参数
Figure 333716DEST_PATH_IMAGE150
，当所述差值累加和
Figure 695427DEST_PATH_IMAGE151
和稳态窗内的方差S满足以下条件时，确定所述差值累加和
Figure 407031DEST_PATH_IMAGE151
所对应的时间点为启动/停机暂态结束点：
Figure 537798DEST_PATH_IMAGE152
（17）
Figure 942235DEST_PATH_IMAGE153
（18）
负荷暂态下一轮判定子模块，用于当检测窗内完成遍历，并且所述差值累加和
Figure 791242DEST_PATH_IMAGE154
和稳态窗方差S满足未满足式（17）和式（18）时，检测窗和稳态窗移至离
Figure 572116DEST_PATH_IMAGE155
最近的波动点进行下一轮判定。
本发明的实施例的一种非介入式负荷辨识的事件检测系统200与本发明的另一个实施例的一种非介入式负荷辨识的事件检测方法100相对应，在此不再赘述。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产品的形式。本申请实施例中的方案可以采用各种计算机语言实现，例如，面向对象的程序设计语言Java和直译式脚本语言JavaScript等。
本发明是参照根据本申请实施例的方法、设备（系统）、和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本申请的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。
显然，本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样，倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内，则本申请也意图包含这些改动和变型在内。

Claims (10)

1.一种非介入式负荷辨识的事件检测方法，其特征在于，包括：
基于实时采集的计量点电压电流波形计算实时的有功功率波形，并获取所述有功功率波形的波动点以及所述波动点的波动时间点集合；
基于所述波动点的波动时间点集合，确定波动时间点之间的时间差，去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值，并将所述波动点均值作为窗宽；
当暂态为开始状态时，根据所述窗宽，计算基态窗内的有功功率均值，并根据所述基态窗内的有功功率均值，确定检测窗内的有功功率累加值，当累加和大于判定阈值时，确定当前时间点为启动/停机暂态开启点；
当暂态为负荷状态时，根据所述窗宽，计算检测窗与稳态窗内的差值累加和以及稳态窗方差，当所述差值累加和大于所述判定阈值，并且所述稳态窗方差小于预定参数与稳态窗内的有功功率稳态值的乘积时，确定当前时间点为启动/停机暂态结束点。
2.根据权利要求1所述的方法，其特征在于，基于实时采集的计量点电压电流波形计算实时的有功功率波形，并获取所述有功功率波形的波动点以及所述波动点 的波动时间点集合，包括：
计算实时的有功功率时序序列
Figure 563455DEST_PATH_IMAGE001
：
Figure 266969DEST_PATH_IMAGE002
（1）
其中
Figure 703897DEST_PATH_IMAGE003
和
Figure 442046DEST_PATH_IMAGE004
为瞬时电压电流采样点值，
Figure 821075DEST_PATH_IMAGE005
为工频周波的采样点数；
对所述有功功率时序序列的有功功率波形进行中值滤波，设定滤波窗宽为q，保留波动上升沿与下降沿，根据所述上升沿与所述下降沿，消除噪声干扰；
获取有功功率时序序列中的波动时间点集合
Figure 277464DEST_PATH_IMAGE006
，其中，n为波动时间点的序列号。
3.根据权利要求1所述的方法，其特征在于，基于所述波动点的波动时间点集合，确定波动时间点之间的时间差，去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值，并将所述波动点均值作为窗宽，包括：
基于所述波动点的波动点序列，确定波动时间点之间的时间差，并对所述时间差进行排序，得到差值集合
Figure 501772DEST_PATH_IMAGE007
Figure 828848DEST_PATH_IMAGE008
（2）
Figure 378778DEST_PATH_IMAGE009
（3）
其中，n为波动时间点的序列号；
去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值
Figure 56884DEST_PATH_IMAGE010
，并将所述波动点均值作为基态窗、监测窗和稳态窗的窗宽：
Figure 84883DEST_PATH_IMAGE011
（4）
其中a为设定的判定阈值。
4.根据权利要求1所述的方法，其特征在于，当暂态为开始状态时，根据所述窗宽，计算基态窗内的有功功率均值，并根据所述基态窗内的有功功率均值，确定检测窗内的有功功率累加值，当累加和大于判定阈值时，确定当前时间点为启动/停机暂态开启点，包括：
对于第一波动时间点
Figure 781312DEST_PATH_IMAGE012
，确定所述第一波动时间点的基态窗内的有功功率和所述第一波动时间点的检测窗内的有功功率为：
Figure 236565DEST_PATH_IMAGE013
（5）
Figure 401967DEST_PATH_IMAGE014
（6）
其中，
Figure 233656DEST_PATH_IMAGE015
为所述第一波动时间点的基态窗内的有功功率，
Figure 800904DEST_PATH_IMAGE016
为基态窗以及检测窗的窗宽，
Figure 427057DEST_PATH_IMAGE017
为在所述第一波动时间点往前数窗宽
Figure 79756DEST_PATH_IMAGE018
的时间，
Figure 449557DEST_PATH_IMAGE019
为在所述第一波动时间点往前数窗宽
Figure 605732DEST_PATH_IMAGE018
再加1秒的时间，
Figure 419098DEST_PATH_IMAGE020
为在所述第一波动时间点往前数窗宽
Figure 293513DEST_PATH_IMAGE021
对应的有功功率，
Figure 467006DEST_PATH_IMAGE022
为与在所述第一波动时间点往前数窗宽
Figure 743266DEST_PATH_IMAGE018
再加1秒的时间对应的功率，
Figure 976802DEST_PATH_IMAGE023
为所述第一波动时间点的有功功率；
Figure 338513DEST_PATH_IMAGE024
为所述第一波动时间点的检测窗内的有功功率，
Figure 315696DEST_PATH_IMAGE025
为第一波动时间点加1秒的时间，
Figure 180884DEST_PATH_IMAGE026
为第一波动时间点加2秒的时间，
Figure 585320DEST_PATH_IMAGE027
为在所述第一波动时间点往后数窗宽
Figure 434328DEST_PATH_IMAGE018
的时间，
Figure 201820DEST_PATH_IMAGE028
为第一波动时间点加1秒的时间对应的有功功率，
Figure 452673DEST_PATH_IMAGE029
为为第一波动时间点加2秒的时间对应的有功功率，
Figure 28010DEST_PATH_IMAGE030
所述第一波动时间点往后数窗宽
Figure 364314DEST_PATH_IMAGE018
的时间对应的有功功率；
根据所述第一波动时间点的基态窗内的有功功率，计算所述第一波动时间点的基态窗内的有功功率基准值：
Figure 683300DEST_PATH_IMAGE031
（7）
其中，
Figure 788659DEST_PATH_IMAGE032
为基态窗内的有功功率基准值，Pk为基态窗内k时间点的瞬时有功功率值；
遍历所述检测窗内的有功功率值，基于CUSUM方法进行改进，计算累加所述检测窗内的有功功率值，确定检测窗内的有功功率累加值：
Figure 269319DEST_PATH_IMAGE033
（8）
Figure 827339DEST_PATH_IMAGE034
（9）
其中，
Figure 950016DEST_PATH_IMAGE035
为检测窗内的有功功率累加值的初始值，
Figure 660614DEST_PATH_IMAGE036
为检测窗内的有功功率累加值；
设定判定有功功率阈值h，当有功功率累加值
Figure 577755DEST_PATH_IMAGE036
满足以下条件时，确定
Figure 623071DEST_PATH_IMAGE036
所对应的时间点为启动/停机暂态开启点：
Figure 283859DEST_PATH_IMAGE037
（10）
当检测窗内完成遍历，并且所述有功功率累加值
Figure 363811DEST_PATH_IMAGE036
未满足式（10）时，基态窗和检测窗移至离
Figure 451853DEST_PATH_IMAGE038
最近的波动点进行下一轮判定。
5.根据权利要求1所述的方法，其特征在于，当暂态为负荷状态时，根据所述窗宽，计算检测窗与稳态窗内的差值累加和以及稳态窗方差，当所述差值累加和大于所述判定阈值，并且所述稳态窗方差小于预定参数与稳态窗内的有功功率稳态值的乘积时，确定当前时间点为启动/停机暂态结束点，包括：
对于第二波动时间点
Figure 984465DEST_PATH_IMAGE039
，确定所述第二波动时间点的检测窗内的有功功率和所述第二波动时间点的稳态窗内的有功功率为：
Figure 448945DEST_PATH_IMAGE040
（11）
Figure 117823DEST_PATH_IMAGE041
（12）
其中，
Figure 891613DEST_PATH_IMAGE042
为所述第二波动时间点的检测窗内的有功功率，
Figure 645943DEST_PATH_IMAGE043
为所述第二波动时间点的稳态窗内的有功功率；
根据所述第二波动时间点的稳态窗内的有功功率，计算稳态窗内的有功功率稳态值和稳态窗内的方差：
Figure 648534DEST_PATH_IMAGE044
（13）
Figure 437498DEST_PATH_IMAGE045
（14）
其中，
Figure 132922DEST_PATH_IMAGE046
为稳态窗内的有功功率稳态值，
Figure 374547DEST_PATH_IMAGE047
为稳态窗内的方差；
设定参数，向前遍历检测窗内的有功功率值，计算检测窗与稳态窗内的差值累加和：
Figure 180829DEST_PATH_IMAGE048
（15）
Figure 89879DEST_PATH_IMAGE049
（16）
设定参数
Figure 425046DEST_PATH_IMAGE050
，当所述差值累加和
Figure 419546DEST_PATH_IMAGE051
和稳态窗内的方差S满足以下条件时，确定所述差值累加和
Figure 514673DEST_PATH_IMAGE052
所对应的时间点为启动/停机暂态结束点：
Figure 278229DEST_PATH_IMAGE053
（17）
Figure 49876DEST_PATH_IMAGE054
（18）
当检测窗内完成遍历，并且所述差值累加
Figure 531673DEST_PATH_IMAGE055
和和稳态窗方差S满足未满足式（17）和式（18）时，检测窗和稳态窗移至离
Figure 679758DEST_PATH_IMAGE056
最近的波动点进行下一轮判定。
6.一种非介入式负荷辨识的事件检测系统，其特征在于，包括：
获取波动点模块，用于基于实时采集的计量点电压电流波形计算实时的有功功率波形，并获取所述有功功率波形的波动点以及所述波动点的波动时间点集合；
确定窗宽模块，用于基于所述波动点的波动时间点集合，确定波动时间点之间的时间差，去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值，并将所述波动点均值作为窗宽；
判定开始状态模块，用于当暂态为开始状态时，根据所述窗宽，计算基态窗内的有功功率均值，并根据所述基态窗内的有功功率均值，确定检测窗内的有功功率累加值，当累加和大于判定阈值时，确定当前时间点为启动/停机暂态开启点；
判定负荷暂态模块，用于当暂态为负荷状态时，根据所述窗宽，计算检测窗与稳态窗内的差值累加和以及稳态窗方差，当所述差值累加和大于所述判定阈值，并且所述稳态窗方差小于预定参数与稳态窗内的有功功率稳态值的乘积时，确定当前时间点为启动/停机暂态结束点。
7.根据权利要求6所述的系统，其特征在于，获取波动点模块，包括：
计算有功功率时序子模块，用于计算实时的有功功率时序序列
Figure 32242DEST_PATH_IMAGE057
：
Figure 240369DEST_PATH_IMAGE058
（1）
其中
Figure 943883DEST_PATH_IMAGE059
和
Figure 630079DEST_PATH_IMAGE060
为瞬时电压电流采样点值，
Figure 632144DEST_PATH_IMAGE061
为工频周波的采样点数；
消除噪声干扰子模块，用于对所述有功功率时序序列的有功功率波形进行中值滤波，设定滤波窗宽为q，保留波动上升沿与下降沿，根据所述上升沿与所述下降沿，消除噪声干扰；
获取波动时间点集合子模块，用于获取有功功率时序序列中的波动时间点集合
Figure 745593DEST_PATH_IMAGE062
，其中，n为波动时间点的序列号。
8.根据权利要求6所述的系统，其特征在于，确定窗宽模块，包括：
得到差值集合子模块，用于基于所述波动点的波动点序列，确定波动时间点之间的时间差，并对所述时间差进行排序，得到差值集合
Figure 936403DEST_PATH_IMAGE063
；
Figure 160711DEST_PATH_IMAGE064
（2）
Figure 753367DEST_PATH_IMAGE065
（3）
其中，n为波动时间点的序列号；
确定窗宽子模块，用于去除最大波动点和最小波动点的干扰项，计算剩余波动点差值的波动点均值
Figure 303297DEST_PATH_IMAGE066
，并将所述波动点均值作为基态窗、监测窗和稳态窗的窗宽：
Figure 981403DEST_PATH_IMAGE067
（4）
其中a为设定的判定阈值。
9.根据权利要求6所述的系统，其特征在于，判定开始状态模块，包括：
确定第一波动时间点有功功率子模块，用于对于第一波动时间点
Figure 9402DEST_PATH_IMAGE068
，确定所述第一波动时间点的基态窗内的有功功率和所述第一波动时间点的检测窗内的有功功率为：
Figure 456563DEST_PATH_IMAGE069
（5）
Figure 928127DEST_PATH_IMAGE070
（6）
其中，
Figure 93529DEST_PATH_IMAGE071
为所述第一波动时间点的基态窗内的有功功率，
Figure 925219DEST_PATH_IMAGE016
为基态窗以及检测窗的窗宽，
Figure 226887DEST_PATH_IMAGE072
为在所述第一波动时间点往前数窗宽
Figure 118620DEST_PATH_IMAGE073
的时间，
Figure 505739DEST_PATH_IMAGE019
为在所述第一波动时间点往前数窗宽
Figure 875540DEST_PATH_IMAGE018
再加1秒的时间，
Figure 297295DEST_PATH_IMAGE020
为在所述第一波动时间点往前数窗宽
Figure 359928DEST_PATH_IMAGE021
对应的有功功率，
Figure 968764DEST_PATH_IMAGE022
为与在所述第一波动时间点往前数窗宽
Figure 657104DEST_PATH_IMAGE073
再加1秒的时间对应的功率，
Figure 933364DEST_PATH_IMAGE023
为所述第一波动时间点的有功功率；
Figure 166899DEST_PATH_IMAGE074
为所述第一波动时间点的检测窗内的有功功率，
Figure 528611DEST_PATH_IMAGE025
为第一波动时间点加1秒的时间，
Figure 974635DEST_PATH_IMAGE026
为第一波动时间点加2秒的时间，
Figure 370982DEST_PATH_IMAGE075
为在所述第一波动时间点往后数窗宽
Figure 775418DEST_PATH_IMAGE018
的时间，
Figure 624425DEST_PATH_IMAGE076
为第一波动时间点加1秒的时间对应的有功功率，
Figure 139720DEST_PATH_IMAGE029
为为第一波动时间点加2秒的时间对应的有功功率，
Figure 141306DEST_PATH_IMAGE030
所述第一波动时间点往后数窗宽
Figure 716643DEST_PATH_IMAGE018
的时间对应的有功功率；
计算有功功率基准值子模块，用于根据所述第一波动时间点的基态窗内的有功功率，计算所述第一波动时间点的基态窗内的有功功率基准值：
Figure 787368DEST_PATH_IMAGE077
（7）
其中，
Figure 106354DEST_PATH_IMAGE078
为基态窗内的有功功率基准值， Pk为基态窗内k时间点的瞬时有功功率值；
确定有功功率累加值子模块，用于遍历所述检测窗内的有功功率值，基于CUSUM方法进行改进，计算累加所述检测窗内的有功功率值，确定检测窗内的有功功率累加值：
Figure 211713DEST_PATH_IMAGE079
（8）
Figure 957952DEST_PATH_IMAGE080
（9）
其中，
Figure 515972DEST_PATH_IMAGE081
为检测窗内的有功功率累加值的初始值，
Figure 638649DEST_PATH_IMAGE082
为检测窗内的有功功率累加值；
确定暂态开启点子模块，用于设定判定有功功率阈值h，当有功功率累加值
Figure 598515DEST_PATH_IMAGE082
满足以下条件时，确定
Figure 767852DEST_PATH_IMAGE082
所对应的时间点为启动/停机暂态开启点：
Figure 78748DEST_PATH_IMAGE083
（10）
开始状态下一轮判定子模块，用于当检测窗内完成遍历，并且所述有功功率累加值
Figure 473957DEST_PATH_IMAGE083
未满足式（10）时，基态窗和检测窗移至离
Figure 288330DEST_PATH_IMAGE084
最近的波动点进行下一轮判定。
10.根据权利要求6所述的系统，其特征在于，判定负荷暂态模块，包括：
确定第二波动时间点有功功率子模块，用于对于第二波动时间点
Figure 641950DEST_PATH_IMAGE085
，确定所述第二波动时间点的检测窗内的有功功率和所述第二波动时间点的稳态窗内的有功功率为：
Figure 174563DEST_PATH_IMAGE086
（11）
Figure 373463DEST_PATH_IMAGE087
（12）
其中，
Figure 307921DEST_PATH_IMAGE088
为所述第二波动时间点的检测窗内的有功功率，
Figure 832443DEST_PATH_IMAGE089
为所述第二波动时间点的稳态窗内的有功功率；
计算稳态窗方差子模块，用于根据所述第二波动时间点的稳态窗内的有功功率，计算稳态窗内的有功功率稳态值和稳态窗内的方差：
Figure 321194DEST_PATH_IMAGE090
（13）
Figure 340096DEST_PATH_IMAGE091
（14）
其中，
Figure 129061DEST_PATH_IMAGE092
为稳态窗内的有功功率稳态值，
Figure 558905DEST_PATH_IMAGE093
为稳态窗内的方差；
计算差值累加和子模块，用于设定参数
Figure 66110DEST_PATH_IMAGE094
，向前遍历检测窗内的有功功率值，计算检测窗与稳态窗内的差值累加和：
Figure 872392DEST_PATH_IMAGE095
（15）
Figure 515863DEST_PATH_IMAGE096
（16）
确定暂态结束点子模块，用于设定参数
Figure 382187DEST_PATH_IMAGE097
，当所述差值累加和
Figure 845530DEST_PATH_IMAGE098
和稳态窗内的方差S满足以下条件时，确定所述差值累加和
Figure 455503DEST_PATH_IMAGE098
所对应的时间点为启动/停机暂态结束点：
Figure 468327DEST_PATH_IMAGE099
（17）
Figure 239974DEST_PATH_IMAGE100
（18）
负荷暂态下一轮判定子模块，用于当检测窗内完成遍历，并且所述差值累加和
Figure 721771DEST_PATH_IMAGE101
和稳态窗方差S满足未满足式（17）和式（18）时，检测窗和稳态窗移至离
Figure 869855DEST_PATH_IMAGE102
最近的波动点进行下一轮判定。
CN202111265840.7A 2021-10-29 2021-10-29 一种非介入式负荷辨识的事件检测方法及系统 Active CN113704698B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111265840.7A CN113704698B (zh) 2021-10-29 2021-10-29 一种非介入式负荷辨识的事件检测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111265840.7A CN113704698B (zh) 2021-10-29 2021-10-29 一种非介入式负荷辨识的事件检测方法及系统

Publications (2)

Publication Number Publication Date
CN113704698A CN113704698A (zh) 2021-11-26
CN113704698B true CN113704698B (zh) 2022-02-11

Family

ID=78647335

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111265840.7A Active CN113704698B (zh) 2021-10-29 2021-10-29 一种非介入式负荷辨识的事件检测方法及系统

Country Status (1)

Country Link
CN (1) CN113704698B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114996182B (zh) * 2022-05-23 2024-04-26 中国计量大学 一种用于负荷识别模型训练的稳态数据流输出装置
CN115907567B (zh) * 2023-02-21 2023-05-09 浙江大学 基于鲁棒性随机分割森林算法的负荷事件检测方法及系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11300616B2 (en) * 2017-07-12 2022-04-12 Mayo Foundation For Medical Education And Research Systems and methods for non-invasive current estimation
CN110504679B (zh) * 2019-07-25 2021-05-18 深圳供电局有限公司 一种基于km匹配算法的非侵入式负荷辨识方法
CN110488128A (zh) * 2019-09-12 2019-11-22 广东电网有限责任公司佛山供电局 双边累积和事件检测方法
CN110954744A (zh) * 2019-11-18 2020-04-03 浙江工业大学 一种基于事件检测的非侵入式负荷监测方法
CN112782465A (zh) * 2020-12-01 2021-05-11 国网浙江省电力有限公司营销服务中心 一种基于滑动窗口技术的暂态负荷启停事件快速捕捉方法
CN112881793B (zh) * 2021-01-13 2022-02-01 吉林大学 一种结合时间阈值的非侵入式负荷事件检测方法

Also Published As

Publication number Publication date
CN113704698A (zh) 2021-11-26

Similar Documents

Publication Publication Date Title
CN113704698B (zh) 一种非介入式负荷辨识的事件检测方法及系统
CN109213654B (zh) 一种异常检测方法及装置
CN107871190B (zh) 一种业务指标监控方法及装置
CN106708694B (zh) 一种定位毛刺异常点的方法、装置及计算设备
CN109115257B (zh) 传感器特性曲线的修正方法、装置、设备及存储介质
CN112149860A (zh) 一种自动异常检测方法和系统
CN114978956B (zh) 智慧城市网络设备性能异常突变点检测方法及装置
CN112396007B (zh) 一种非侵入式居民负荷突变事件三阈值检测方法和系统
CN111161097B (zh) 基于假设检验的事件检测算法检测开关事件的方法和装置
CN110601157B (zh) 一种间歇性单相接地故障保护方法及装置
CN108762908B (zh) 系统调用异常检测方法及装置
CN110858177A (zh) 内存使用量监测方法及基板处理装置
CN115167281A (zh) 直流氩弧焊工艺执行告警方法、装置、电子设备及介质
CN110907843A (zh) 电池阻抗计算方法及装置
CN112445679B (zh) 一种信息检测方法、装置、服务器及存储介质
CN110196354B (zh) 一种负载的开关事件的检测方法和装置
CN105425683B (zh) 一种基于阶梯式动态死区的模拟量采集方法
CN114184326A (zh) 天然气泄漏传感器状态识别方法及系统
CN113810743A (zh) 一种直播间流量监控方法和装置
KR101946826B1 (ko) 신호 처리 장치 및 방법
CN114083987A (zh) 电池监测参数的修正方法、装置和计算机设备
CN106709251A (zh) 一种评估方法及装置
CN115112962A (zh) 一种非侵入式负荷辨识方法、系统、存储介质及设备
CN114367547B (zh) 一种轧制数据的统计方法和装置
CN118091325B (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