CN104036362A - 一种台变用电负荷异常数据的快速检测方法 - Google Patents
一种台变用电负荷异常数据的快速检测方法 Download PDFInfo
- Publication number
- CN104036362A CN104036362A CN201410283769.9A CN201410283769A CN104036362A CN 104036362 A CN104036362 A CN 104036362A CN 201410283769 A CN201410283769 A CN 201410283769A CN 104036362 A CN104036362 A CN 104036362A
- Authority
- CN
- China
- Prior art keywords
- power load
- data
- abnormal data
- platform
- sigma
- 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.)
- Pending
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种台变用电负荷异常数据的快速检测方法,该方法首先对所采集到的数据资料进行预处理,得到初始样本集;对样本集进行一阶向后差分;使用差分序列建立Gaussian分布模型,采用最小二乘估计法求解模型参数,在给定置信水平下求解置信区间;将待检测用电负荷数据的差分序列代入模型,当超出置信区间时判定该时间点为异常值。本发明提供一种准确性高、可操作性高的台变用电负荷的异常数据检测方法,实现台变用电负荷异常数据的快速检测及挖掘。
Description
技术领域
本发明涉及一种台变用电负荷异常数据的快速检测方法,属台变运行检测技术领域。
背景技术
在正常情况下,电网的负荷呈现出一定的有规律的趋势。但是,在某些特殊的情况下,负荷曲线会出现巨大的波动,而这种负荷数据上的巨大波动有可能是采集与传输上误差所致,也可能标示着负荷的剧变。因此对负荷数据的剧变进行分析、挖掘与预测是一个非常重要的内容,它可以给调度管理部门提前给一个估计,以便制定紧急事故下的运行管理方式及拉闸限电序位表,以防止电网崩溃与瓦解。
开展异常数据检测与挖掘是比较棘手的工作,所得的结果少而散,并且由于异常数据的多样性与复杂性导致并不存在着异常数据检测的通用算法,需要针对不同类型的数据提出有效的检测方法。
发明内容
本发明的目的是,为了实现台变用电负荷异常数据的快速检测及挖掘,本发明提出一种台变用电负荷异常数据的快速检测方法。
本发明方法包括下述步骤:
(1)对实测数据进行数据预处理,剔除实测数据中的缺损数据,建立训练样本集;
(2)对所建立的样本数据集进行差分处理,得到台变用电负荷数据的差分序列;
(3)使用步骤(2)得到的台变用电负荷数据的差分序列,构造Gaussian分布模型,确定Gaussian分布模型参数;
(4)计算给定置信度水平下步骤(3)中Gaussian分布的置信区间,作为台变用电负荷异常数据的判定阈值;
(5)利用步骤(4)所建立的分布模型进行台变用电负荷的异常数据检测。
本发明步骤(2)中,差分处理是一阶向后差分,其数学表达式为:
Δx(t)=x(t)-x(t-1),t=2,3,...
其中,x(t)是t时刻台变用电负荷值,Δx(t)是t时刻台变用电负荷的一阶向后差分值。
本发明步骤(4)中,Gaussian分布模型的数学表达式为:
其中f为Gaussian分布的概率密度函数,μ是位置参数,σ是尺度参数。
本发明步骤(4)中,确定Gaussian分布参数的方法为最小二乘估计法。
本发明步骤(4)的计算给定置信水平下Gaussian分布的置信区间的方法包括:
(4.1)估计样本标准差σ,其数学表达式为:
其中,s为总体标准差,n为样本个数。
(4.2)当抽取的样本数量足够大时,根据中心极限定理,认为样本均值近似地服从Gaussian分布。其Z统计量的计算公式为:
其中为样本均值,μ为总体均值。
(4.3)在给定1-α的置信水平下,总体均值μ的置信区间为:
其中Zα/2的值通过查表得到。
本发明所述异常数据的判断方法为:对待判断的台变用电负荷数据z(t),t=1,2,...,M,计算其一阶向前差分序列Δz(t),t=2,3...,M;当t时刻的差分值Δz(t)落入步骤4所求置信区间内时,则该时刻台变用电负荷数据为正常值;反之,当t时刻的差分值Δz(t)超出步骤4所求置信区间时,则该时刻台变用电负荷数据为异常值。
本发明的效果在于:本发明所述的用于台变用电负荷数据的异常数据检测方法,能够提供一种准确性高、可操作性高的台变用电负荷的异常数据检测方法,实现台变用电负荷异常数据的快速检测及挖掘。
本发明适用于台变用电负荷异常数据的快速检测。
附图说明
图1是台变用电负荷异常数据检测流程图;
图2是采用本方法的实施例进行的一个实际应用试验的台变用电负荷异常数据检测的结果图。
具体实施方式
下面结合附图,对优选实施例作详细说明。
图1是本发明台变用电负荷异常数据检测的流程图。图1中,本发明提供的方法包括以下步骤:
步骤1:数据预处理,建立训练样本集和测试样本集。
对于所选用的实施例数据,首先剔除采集资料中的缺损数据,将整个样本集按照时间序列顺序排列,建立初始样本集。将初始样本集划分为训练样本集和测试样本集。其中,训练样本集用于台变用电负荷的异常数据检测模型的建立和模型参数估计,测试样本集为待检测的样本数据段。
本实施例中,数据采集间隔为半小时,则一天的数据量为48点。初始样本集为***中所有采集数据经过剔除缺损数据得到,测试样本集为48×7天=336点。
步骤2:对训练样本集进行一阶向后差分处理。
假设训练样本集的时间序列为x(t),t=1,2,...,N,则其一阶向后差分序列的数学表达式为:
y(t)=Δx(t)=x(t)-x(t-1),t=2,3,...N
步骤3:对差分序列y(t),t=2,3,...,N建立Gaussian分布模型。
步骤3.1:假设差分的总体均值为μ,则差分序列y(t),t=2,3,...,N可分解为均值μ与误差e(t)之和:
y(t)=μ+e(t),t=2,3,...,N
步骤3.2:使用最小二乘法估计总体平均数μ,就是使y(t)与μ之间的误差平方和为最小,即使得下式最小:
求Q对μ的偏导数,并令其等于0,得:
则得到总体平均数的估计量为
步骤3.3:估计总体样本方差σ2,数学表达式如下:
因而,σ2的估计为:
步骤3.4:差分序列Δx(t),t=2,3,...,N建立Gaussian分布模型为
本实施例中,所求得的Gaussian分布模型为N(0.1278,0.16602)。
步骤4:给定置信水平1-α,计算步骤3中所求Gaussian分布对应的置信区间。
步骤4.1:由步骤3,假设则有:
步骤4.2:根据给定的置信水平1-α,查表得zα/2的值。
步骤4.3:由Gaussian概率分布的定义,可得如下数学表达式
则有
随机区间即为置信水平为1-α的置信区间。本实施例中,选取1-α=0.95,求得步骤3中所得Gaussian分布的置信区间为(0.1222,0.1334)。
步骤5:测试样本集的异常数据检测。假设待检测的测试样本集为z(t),t=1,2,...,M,首先计算其一阶向后差分序列Δz(t)=z(t)-z(t-1),t=2,3,...,M。当t时刻的差分值Δz(t)落入步骤4所求置信区间内时,则该时刻台变用电负荷数据为正常值;反之,当t时刻的差分值Δz(t)超出步骤4所求置信区间时,则该时刻台变用电负荷数据为异常值。
本实施例中,待检测的测试样本集选取为24×7天=168点,所得到的异常数据检测结果如图2所示。
Claims (6)
1.一种台变用电负荷异常数据的快速检测方法,其特征在于,所述方法包括下述步骤:
(1)对实测数据进行数据预处理,剔除实测数据中的缺损数据,建立训练样本集;
(2)对所建立的样本数据集进行差分处理,得到台变用电负荷数据的差分序列;
(3)使用步骤(2)得到的台变用电负荷数据的差分序列,构造Gaussian分布模型,确定Gaussian分布模型参数;
(4)计算给定置信度水平下步骤(3)中Gaussian分布的置信区间,作为台变用电负荷异常数据的判定阈值;
(5)利用步骤(4)所建立的分布模型进行台变用电负荷的异常数据检测。
2.根据权利要求1所述的一种台变用电负荷异常数据的快速检测方法,其特征在于,所述差分处理是一阶向后差分,其数学表达式为:
Δx(t)=x(t)-x(t-1),t=2,3,...
其中,x(t)是t时刻台变用电负荷值,Δx(t)是t时刻台变用电负荷的一阶向后差分值。
3.根据权利要求1所述的一种台变用电负荷异常数据的快速检测方法,其特征在于,所述Gaussian分布模型的数学表达式为:
其中f为Gaussian分布的概率密度函数,μ是位置参数,σ是尺度参数。
4.根据权利要求1所述的一种台变用电负荷异常数据的快速检测方法,其特征在于,所述确定Gaussian分布参数的方法为最小二乘估计法。
5.根据权利要求1所述的一种台变用电负荷异常数据的快速检测方法,其特征在于,所述计算给定置信水平下Gaussian分布的置信区间的方法包括:
(4.1)估计样本标准差σ,其数学表达式为:
其中,s为总体标准差,n为样本个数;
(4.2)当抽取的样本数量足够大时,根据中心极限定理,认为样本均值近似地服从Gaussian分布;其Z统计量的计算公式为:
其中为样本均值,μ为总体均值;
(4.3)在给定1-α的置信水平下,总体均值μ的置信区间为:
其中Zα/2的值通过查表得到。
6.根据权利要求1所述的一种台变用电负荷异常数据的快速检测方法,其特征在于,所述异常数据的判断方法为:对待判断的台变用电负荷数据z(t),t=1,2,...,M,计算其一阶向前差分序列Δz(t),t=2,3...,M当t时刻的差分值Δz(t)落入步骤4所求置信区间内时,则该时刻台变用电负荷数据为正常值;反之,当t时刻的差分值Δz(t)超出步骤4所求置信区间时,则该时刻台变用电负荷数据为异常值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410283769.9A CN104036362A (zh) | 2014-06-23 | 2014-06-23 | 一种台变用电负荷异常数据的快速检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410283769.9A CN104036362A (zh) | 2014-06-23 | 2014-06-23 | 一种台变用电负荷异常数据的快速检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104036362A true CN104036362A (zh) | 2014-09-10 |
Family
ID=51467125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410283769.9A Pending CN104036362A (zh) | 2014-06-23 | 2014-06-23 | 一种台变用电负荷异常数据的快速检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104036362A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105205736A (zh) * | 2015-10-14 | 2015-12-30 | 国家电网公司 | 基于经验模态分解的电力负荷异常数据的快速检测方法 |
WO2016101690A1 (zh) * | 2014-12-22 | 2016-06-30 | 国家电网公司 | 基于时间序列分析的输变电设备的状态监测数据清洗方法 |
CN107271764A (zh) * | 2017-06-19 | 2017-10-20 | 宁波三星医疗电气股份有限公司 | 一种用电器耗电异常检测方法及装置 |
CN107478988A (zh) * | 2017-09-21 | 2017-12-15 | 山东大学 | 基于非精确贝叶斯模型的断路器异常判别方法和*** |
CN109242321A (zh) * | 2018-09-17 | 2019-01-18 | 国网河北省电力有限公司电力科学研究院 | 用户电力负荷在线分析方法及终端设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0991045A (ja) * | 1995-09-25 | 1997-04-04 | Nec Corp | 電源電圧監視装置 |
CN102073785A (zh) * | 2010-11-26 | 2011-05-25 | 哈尔滨工程大学 | 基于广义动态模糊神经网络的燃气日负荷组合预报方法 |
CN102930178A (zh) * | 2012-11-26 | 2013-02-13 | 慈溪市供电局 | 一种负荷数据异常检测方法 |
-
2014
- 2014-06-23 CN CN201410283769.9A patent/CN104036362A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0991045A (ja) * | 1995-09-25 | 1997-04-04 | Nec Corp | 電源電圧監視装置 |
CN102073785A (zh) * | 2010-11-26 | 2011-05-25 | 哈尔滨工程大学 | 基于广义动态模糊神经网络的燃气日负荷组合预报方法 |
CN102930178A (zh) * | 2012-11-26 | 2013-02-13 | 慈溪市供电局 | 一种负荷数据异常检测方法 |
Non-Patent Citations (2)
Title |
---|
康丽峰: "基于数据挖掘的短期负荷预测方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
陈斌等: "异常检测综述", 《山东大学学报(工学版)》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016101690A1 (zh) * | 2014-12-22 | 2016-06-30 | 国家电网公司 | 基于时间序列分析的输变电设备的状态监测数据清洗方法 |
CN105205736A (zh) * | 2015-10-14 | 2015-12-30 | 国家电网公司 | 基于经验模态分解的电力负荷异常数据的快速检测方法 |
CN107271764A (zh) * | 2017-06-19 | 2017-10-20 | 宁波三星医疗电气股份有限公司 | 一种用电器耗电异常检测方法及装置 |
CN107271764B (zh) * | 2017-06-19 | 2020-05-12 | 宁波三星医疗电气股份有限公司 | 一种用电器耗电异常检测方法及装置 |
CN107478988A (zh) * | 2017-09-21 | 2017-12-15 | 山东大学 | 基于非精确贝叶斯模型的断路器异常判别方法和*** |
CN109242321A (zh) * | 2018-09-17 | 2019-01-18 | 国网河北省电力有限公司电力科学研究院 | 用户电力负荷在线分析方法及终端设备 |
CN109242321B (zh) * | 2018-09-17 | 2020-10-09 | 国网河北省电力有限公司电力科学研究院 | 用户电力负荷在线分析方法及终端设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104036362A (zh) | 一种台变用电负荷异常数据的快速检测方法 | |
CN105205736A (zh) | 基于经验模态分解的电力负荷异常数据的快速检测方法 | |
EP3575960B1 (en) | System upgrade assessment method based on system parameter correlation coefficients | |
CN104766175A (zh) | 一种基于时间序列分析的电力***异常数据辨识与修正方法 | |
CN103020166B (zh) | 一种电力实时数据异常检测方法 | |
CN110648480B (zh) | 一种基于变化速率的单变量报警***及方法 | |
CN105279365A (zh) | 用于学习异常检测的样本的方法 | |
CN111563524A (zh) | 一种多站融合***运行态势异常监测和告警合并方法 | |
CN108304348B (zh) | 一种基于二元维纳过程的轴承剩余寿命预测方法 | |
CN107424441B (zh) | 基于Hotelling′s T2统计量的航空器航迹变点检测与估计方法 | |
CN109815627B (zh) | 一种径流变异影响因子确定方法及*** | |
CN101841155B (zh) | 用于电力***暂态稳定分析的典型故障集识别方法 | |
CN102721987A (zh) | 多普勒雷达遥感强风暴的预警方法 | |
CN104899448B (zh) | 一种集合卡尔曼滤波静态局地化方案的自适应补偿方法 | |
CN104462794A (zh) | 一种基于对比统计分析发现环境监测异常数据的算法 | |
CN110334406B (zh) | 一种考虑风速特大值的极值风速重现期确定方法和装置 | |
CN113889252B (zh) | 基于生命体征大数据聚类核心算法和区块链的远程互联网大数据智慧医疗*** | |
CN110990788A (zh) | 一种基于三元维纳过程的轴承剩余寿命预测方法 | |
CN110889554A (zh) | 基于重现时间间隔分析法的电力负荷波动分析和风险预警方法 | |
CN103678938B (zh) | 一种面向空间形状和误差范围的退化模型一致性检验方法 | |
CN117394337A (zh) | 一种电网负荷预警方法及其*** | |
CN103134433B (zh) | 一种利用位移监测鉴别边坡失稳致滑因子的方法 | |
CN109142979B (zh) | 一种配电网异常状态的检测方法及装置 | |
CN106442830B (zh) | 变压器油中气体含量告警值的检测方法和*** | |
CN102710349B (zh) | 一种脉冲干扰环境下基于数据挑选的频谱感知方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140910 |
|
RJ01 | Rejection of invention patent application after publication |