CN105574669B - 时空联合数据聚类分析的输变电设备状态异常检测方法 - Google Patents
时空联合数据聚类分析的输变电设备状态异常检测方法 Download PDFInfo
- Publication number
- CN105574669B CN105574669B CN201510943246.7A CN201510943246A CN105574669B CN 105574669 B CN105574669 B CN 105574669B CN 201510943246 A CN201510943246 A CN 201510943246A CN 105574669 B CN105574669 B CN 105574669B
- Authority
- CN
- China
- Prior art keywords
- time
- time window
- data
- space
- class
- 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
- 230000005540 biological transmission Effects 0.000 title claims abstract description 30
- 230000009466 transformation Effects 0.000 title claims abstract description 27
- 238000004458 analytical method Methods 0.000 title claims abstract description 8
- 238000001514 detection method Methods 0.000 title abstract description 17
- 230000002159 abnormal effect Effects 0.000 claims abstract description 51
- 239000011159 matrix material Substances 0.000 claims abstract description 39
- 238000011161 development Methods 0.000 claims abstract description 7
- 238000000034 method Methods 0.000 claims description 16
- 230000005856 abnormality Effects 0.000 claims description 11
- 238000007621 cluster analysis Methods 0.000 claims description 5
- 230000002123 temporal effect Effects 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 13
- 238000012544 monitoring process Methods 0.000 description 11
- 238000011156 evaluation Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 230000018109 developmental process Effects 0.000 description 4
- 238000009413 insulation Methods 0.000 description 4
- 230000000903 blocking effect Effects 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 239000011248 coating agent Substances 0.000 description 2
- 238000000576 coating method Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000013210 evaluation model Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 241000238413 Octopus Species 0.000 description 1
- 240000008866 Ziziphus nummularia Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011208 chromatographic data Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004868 gas analysis Methods 0.000 description 1
- 238000009421 internal insulation Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/28—Databases characterised by their database models, e.g. relational or object models
- G06F16/284—Relational databases
- G06F16/285—Clustering or classification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/30—Information retrieval; Database structures therefor; File system structures therefor of unstructured textual data
- G06F16/35—Clustering; Classification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/30—Information retrieval; Database structures therefor; File system structures therefor of unstructured textual data
- G06F16/35—Clustering; Classification
- G06F16/355—Class or cluster creation or modification
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- General Engineering & Computer Science (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Data Mining & Analysis (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
本发明公开了时空联合数据聚类分析的输变电设备状态异常检测方法,包括以下步骤:对输变电设备状态数据的时间分量添加时窗获得时间子序列,通过结合空间分量以及生成的时间子序列,得到时空子序列;通过使用c均值模糊聚类FCM对时空子序列按不同时窗进行聚类,得到分块矩阵;每个分块矩阵描述分块矩阵对应的时窗中存在的类;根据输变电设备状态历史正常数据计算单个时空子序列的异常度值,并利用分块矩阵对每个时窗的类赋予异常度值,根据异常度值大小判断设备状态的异常情况;如果异常度值高,则说明设备状态存在异常情况;利用分块矩阵建立不同时窗中的各个类之间的关系,从而形象化地观察异常类的发展情况。
Description
技术领域
本发明涉及输变电设备检测技术领域,具体是一种时空联合数据聚类分析的输变电设备状态异常检测方法。
背景技术
输变电设备在实际运行过程中会受到过负荷、过电压、内部绝缘老化、自然环境等异常事件影响,这些异常运行状态会导致设备缺陷、故障的发生,所以对输变电设备状态进行异常检测具有很强的必要性。
国内外文献提出了许多输变电设备的异常检测与状态评估技术。《基于MCD稳健统计分析的变压器油色谱异常值检测及分布特性》采用MCD稳健统计方法进行异常值检测,并将异常值和正常数据进行单独分析,给出正常油色谱数据的分布规律,《基于贝叶斯网络的电力变压器状态评估》建立贝叶斯网络对变压器进行状态评估,《基于可拓分析法的电力变压器本体绝缘状态评估》对变压器本体绝缘状态进行等级划分的基础上,通过构造适当的评估指标关联函数,进而建立基于物元可拓评价法的变压器本体绝缘状态评估模型,《基于云理论的电力变压器绝缘状态评估方法》结合油中气体分析、红外色像仪技术与云模型方法对变压器绝缘状态进行评价,《人工神经网络和信息融合技术在变压器状态评估中的应用》综合应用人工神经网络和Dempster-Shafer证据理论构建多信息融合的变压器综合状态评估模型。
目前的输变电设备异常检测方法利用的数据量小,没有充分考虑到历史数据的作用,并且只考虑了时间维度的数据,没有将时间大数据与地理空间信息结合。设备的状态数据(包括在线监测、带电检测、预防性试验数据等)具有体量大,种类繁多等特点,常规的数据处理技术难以对设备状态数据进行全面、深入的分析。近年来大数据技术不断发展,将大数据技术应用于输变电设备异常检测,可大大提高设备状态数据的利用率,有利于更加全面,深刻地分析设备状态。输变电设备的某些状态数据与地理空间位置存在一定的关联性,将空间维度信息与时间维度信息结合考虑,可以更加全面地分析设备的异常情况。
发明内容
本发明的目的就是为了解决上述问题,提供一种时空联合数据聚类分析的输变电设备状态异常检测方法,它具有实现状态数据时间与空间维度信息的关联性挖掘的优点。
为了实现上述目的,本发明采用如下技术方案:
时空联合数据聚类分析的输变电设备状态异常检测方法,包括以下步骤:
步骤S1:对输变电设备状态数据的时间分量添加时窗获得时间子序列,通过结合空间分量以及生成的时间子序列,得到时空子序列;
步骤S2:通过使用c均值模糊聚类FCM对时空子序列按不同时窗进行聚类,得到分块矩阵;进入步骤S3;
步骤S3:根据输变电设备状态历史正常数据计算单个时空子序列的异常度值,并利用步骤S2的分块矩阵对每个时窗的类赋予异常度值,根据异常度值大小判断设备状态的异常情况。
时空联合数据聚类分析的输变电设备状态异常检测方法,还包括:
步骤S4:利用步骤S2的分块矩阵建立不同时窗中的各个类之间的模糊关系,并获得相邻两个时窗之间的模糊关系矩阵,通过模糊关系矩阵观察异常类的发展情况。
所述步骤S2:每个分块矩阵中的元素描述分块矩阵对应的时窗中存在的类;
所述步骤S3:如果异常度值越高,则说明设备状态存在异常情况的可能性越大。
所述步骤S4:每个时窗中都会进行聚类,将每个时窗中的时空子序列分为若干类,每个时窗对应着分块矩阵,每个时窗中的类对应着分块矩阵中的元素,利用分块矩阵建立不同时窗中的各个类之间的模糊关系,并获得相邻两个时窗之间的模糊关系矩阵R;
模糊关系矩阵R的第i行第j列的元素ri,j代表第一个时窗W1中的第i类与第二个时窗W2中的第j类的相关程度,
rs,t的值越大,则说明第一个时窗W1中的第i类与第二个时窗W2中的第j类的相关程度越大,
相关程度大说明第二个时窗W2中的第j类是由第一个时窗W1中的第i类发展而来的,由此得出一个时窗中的异常类与相邻时窗中对应的相关类,从而了解异常类的发展情况。
所述步骤S1具体包括:
假设有N个数据x1,x2,...,xN,每个数据由一个空间分量和一个时间分量组成。
对于第k个数据xk,有:
xk=[xk(s)|xk(t)], (1)
其中,xk(s)代表数据的空间分量,xk(t)代表数据的时间分量。
假设空间分量的维数为r(通常r=2),时间分量的维数为q,则得到维数为n=r+q的第k个数据的表达式如公式(2)所示:
xk=[xk1(s),...,xkr(s)|xk1(t),...,xkq(t)]T, (2)
其中,xkr(s)表示第k个数据空间分量的第r个元素,xkq(t)表示第k个数据时间分量的第q个元素。
假设时窗的长度为l,则第j个时窗中第k个数据的时间子序列xkj(t)的表达式如公式(3)所示:
xkj(t)=[xk((j-1)*l+1)(t),...,xk(j*l)(t)]T, (3)
将第j个时窗中第k个数据的时间子序列与空间分量结合,得到第j个时窗中第k个数据的时空子序列xkj的表达式如公式(4)所示:
xkj=[xk1(s),...,xkr(s)|xk((j-1)*l+1)(t),...,xk(j*l)(t)]T, (4)
所述步骤S2具体包括:
利用c均值模糊聚类方法,计算聚类中心以及分块矩阵:
其中,vi为第i个聚类中心,uik为代表分块矩阵的元素,表示第k个数据到第i个聚类中心的隶属度,
dλ(vi,xk)表示第i个聚类中心到第k个数据的时空复合距离。λ表示时间分量所占的权重,vi(s)表示第i个聚类中心的空间分量,vi(t)表示第i个聚类中心的时间分量。
通过求得的聚类中心以及分块矩阵对第k个数据进行重构:
在设定范围内尝试多个不同的λ值,利用重建误差E(λ)来评估聚类的质量,选择一个使重建误差最小的λ值;
所述步骤S3具体包括:
假设xk是待检测数据,yk是历史正常数据,xkj是时空序列xk在时窗Wj中的子序列;ykj是时空序列yk在时窗Wj中的子序列,用公式(9)计算子序列xkj的异常度值:
fkj=||xkj-ykj||2, (9)
其中,fkj表示第j个时窗中第k个数据时空子序列的异常度值。
假设U是对时窗Wj中时空子序列聚类得到的分块矩阵,对时窗Wj中的每个子序列的异常度值进行加权平均得到在时窗Wj中每个聚类中心对应的类的异常度值sj={si,i=1,2,...,cj},其中
cj为在时窗Wj中类的数目。
si的值越大,则在时窗Wj中从属于第i类的子序列异常度越高;异常度值越高说明设备出现异常的可能性越大。
si的值越小,则在时窗Wj中从属于第i类的子序列异常度越低;异常度值越低说明设备出现异常的可能性越小。
所述步骤S4具体包括:
假设U1为时窗W1的分块矩阵,U1的维数为c1×N,U2为时窗W2的分块矩阵,U2的维数为c2×N,其中N是数据的个数,c1是U1中类的个数,c2是U2中类的个数。
构建一个维数为c1×c2的模糊关系矩阵
R=[ri,j],i=1,2,...,c1,j=1,2,...,c2, (11)
ri,j表示时窗W1中的第i类与时窗W2中的第j类的相关程度。
用迭代的方式更新R中的项,R中的第s行第t列的元素rs,t表示为:
其中,rs,t(iter+1)表示元素rs,t示第iter+1次的迭代结果。
< >表示rs,t(iter+1)的值限制在[0,1]的区间内,α为正的学习速度,iter为迭代序号。
通过选用极大极小复合算子,得到:
其中:
uyj,k表示分块矩阵U2第j行第k列的元素,uxi,k表示分块矩阵U1第i行第k列的元素。
rs,t表示W1中的第s类与W2中的第t类之间的相关程度,rs,t的值越大,则说明两个类越相关。
本发明的有益效果:
结合了设备状态数据的时间和空间信息,综合考虑其历史数据,可以更加全面地分析状态数据的异常情况,具有较高的准确性。模糊关系可以从时间发展的维度上观察异常类的演化情况,能将异常类别的演化情况通过矩阵的形式表征,使得状态数据的结构形象化、可视化,有利于更好地理解数据,分析数据。
附图说明
图1(a)为济南地区在线监测导线温度数据图;
图1(b)为安丘地区在线监测导线温度数据图;
图1(c)为济南地区在线监测环境温度数据图;
图1(d)为安丘地区在线监测环境温度数据图;
图1(e)为济南地区在线监测平均风速数据图;
图1(f)为安丘地区在线监测平均风速数据图;
图1(g)为济南地区在线监测湿度数据图;
图1(h)为安丘地区在线监测湿度数据图;
图2(a)为时窗1时空联合数据聚类结果图;
图2(b)为时窗2时空联合数据聚类结果图;
图2(c)为时窗3时空联合数据聚类结果图;
图2(d)为时窗4时空联合数据聚类结果图;
图2(e)为时窗5时空联合数据聚类结果图;
图2(f)为时窗6时空联合数据聚类结果图;
图2(g)为时窗7时空联合数据聚类结果图;
图2(h)为时窗8时空联合数据聚类结果图;
图3为异常度值与模糊关系的图形表示图;
图4为本发明的方法流程图。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
如图4所示,时空联合数据聚类分析的输变电设备状态异常检测方法,包括以下步骤:
步骤S1:对输变电设备状态数据的时间分量添加时窗获得时间子序列,通过结合空间分量以及生成的时间子序列,得到时空子序列;
步骤S2:通过使用c均值模糊聚类FCM对时空子序列按不同时窗进行聚类,得到分块矩阵;每个分块矩阵描述分块矩阵对应的时窗中存在的类;分别进入步骤S3和步骤S4;
步骤S3:根据输变电设备状态历史正常数据计算单个时空子序列的异常度值,并利用步骤S2的分块矩阵对每个时窗的类赋予异常度值,根据异常度值大小判断设备状态的异常情况;如果异常度值高,则说明设备状态存在异常情况;
步骤S4:利用步骤S2的分块矩阵建立不同时窗中的各个类之间的关系,从而形象化地观察异常类的发展情况。
以山东11个地区(济南,安丘,滨州,昌邑,德州,枣庄,章丘,招远,诸城,淄博,邹城)某些输电线路为例,选取其2013年12月1日-2014年2月28日90天的环境温度,导线温度,湿度,平均风速的在线监测数据作为历史正常数据,用于计算异常度值,同时选取其2014年12月1日-2015年2月28日的在线监测数据作为异常检测的样本,综合地理位置信息,判断线路的覆冰异常情况。部分地区在线监测数据如图1(a)-图1(h)所示。从图中可直接看出,在第44-58个数据点,安丘地区的导线温度与环境温度都有明显降低,平均风速与湿度有比较大的升高。
采用时窗长度为20,移动步长为10,第一个时窗为1-20天,第二个时窗为11-30天,以此类推,共8个时窗。将各个地区的经纬度坐标转化到笛卡尔坐标系,在每一个时间序列之后添加空间坐标,构成时空联合数据。在每个时窗中,通过使用1.2中的重构标准方法,选取[0,40]范围内的不同λ值,得到不同的重构误差,选取各时窗的最优λ值如表1所示:
表1不同时窗中的最优λ值
使用c均值模糊聚类方法对各时窗中的时空联合数据进行聚类,得到的聚类结果如图2(a)-图2(h)所示:
由图2(a)-图2(h)可知,在各时窗中,在线监测时空联合数据被分为了3类,位于山东较北部的安丘,昌邑,招远,诸城分为了一类,位于山东纬度靠近中部的滨州和淄博分为了一类,位于山东较南部的济南,德州,枣庄,章丘,邹城分为了一类。将异常检测样本数据与其前一年相对应时窗中的时空联合数据的欧式距离作为每一个序列的异常度值,进而加权得到每一类的异常度值如表2所示。
表2不同时窗中各类的异常度值
由表2可知,第四个时窗中的第二类,第五个时窗中的第三类以及第六个时窗中的第二类均有较高的异常度值,可知数据出现异常的地点为安丘、昌邑、招远、诸城,发生异常的时间点为第41-60个数据点,通过观察异常地区第五个时窗的数据,发现第44-58个数据点的导线温度较低,环境温度较低,平均风速较高,湿度较高,很有可能出现覆冰情况。
由各个时窗中的分块矩阵U1~U8可以得到相邻时窗中各类的模糊关系如表3所示。由时窗W1到W2的模糊关系矩阵表明在时窗W1中的第一、二、三类分别演化为时窗W2中的第三、一、二类。
表3相邻时窗中各类的模糊关系
不同时窗中各类的演化进程及异常度值如图3所示,其中节点表示类,箭头表示不同类之间的关系,每一列节点表示一个时窗,节点上方的数字表示该类的异常度值。
图3可以用来追踪正常与异常类的发展。例如时窗W3中的第三类演化为了时窗W4的第三类中的第二类,该类发生了异常,进而演化出了时窗W5中的具有很高异常度的第三类,时窗W5中的第三类演化出时窗W6中的第二类,该类同样具有较高异常度,时窗W6中的第二类演化出了时窗W7中的第一类,该类已经不再异常,可见覆冰异常情况在时窗W4(对应第31到50个数据点)开始出现,在时窗W7(对应第61到80个数据点)恢复正常。
由此可以得出异常检测的结论:在2014年12月31日到2015年1月29日这段时间中,安丘、昌邑、招远、诸城地区的环境温度、导线温度、平均风速与湿度发生异常;由不同时窗中类之间的模糊关系可知,在2014年12月1日到12月30日,上述地区的微气象参数比较正常,但在12月31日到2015年1月29日这段时间,以上地区全部出现微气象参数异常的情况,2015年1月30日之后,微气象参数恢复正常。因此,在2014年12月31日到2015年1月29日这段时间中,安丘、昌邑、招远、诸城地区很有可能出现覆冰异常的情况。
根据上述地区的输电线路巡视记录可知,在2015年1月13日到1月27日这段时间,上述地区输电线路存在覆冰情况。由此算例可得,本发明方法的异常检测结论与实践情况大致相符,具有一定的实用性。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (4)
1.时空联合数据聚类分析的输变电设备状态异常检测方法,其特征是,包括以下步骤:
步骤S1:对输变电设备状态数据的时间分量添加时窗获得时间子序列,通过结合空间分量以及生成的时间子序列,得到时空子序列;
所述步骤S1具体包括:
假设有N个数据x1,x2,...,xN,每个数据由一个空间分量和一个时间分量组成;
对于第k个数据xk,有:
xk=[xk(s)|xk(t)], (1)
其中,xk(s)代表数据的空间分量,xk(t)代表数据的时间分量;
假设空间分量的维数为r,时间分量的维数为q,则得到维数为n=r+q的第k个数据的表达式如公式(2)所示:
xk=[xk1(s),...,xkr(s)|xk1(t),...,xkq(t)]T, (2)
其中,xkr(s)表示第k个数据空间分量的第r个元素,xkq(t)表示第k个数据时间分量的第q个元素;
假设时窗的长度为l,则第j个时窗中第k个数据的时间子序列xkj(t)的表达式如公式(3)所示:
xkj(t)=[xk((j-1)*l+1)(t),...,xk(j*l)(t)]T, (3)
将第j个时窗中第k个数据的时间子序列与空间分量结合,得到第j个时窗中第k个数据的时空子序列xkj的表达式如公式(4)所示:
xkj=[xk1(s),...,xkr(s)|xk((j-1)*l+1)(t),...,xk(j*l)(t)]T, (4)
步骤S2:通过使用c均值模糊聚类FCM对时空子序列按不同时窗进行聚类,得到分块矩阵;进入步骤S3;
步骤S3:根据输变电设备状态历史正常数据计算单个时空子序列的异常度值,并利用步骤S2的分块矩阵对每个时窗的类赋予异常度值,根据异常度值大小判断设备状态的异常情况;
步骤S4:利用步骤S2的分块矩阵建立不同时窗中的各个类之间的模糊关系,并获得相邻两个时窗之间的模糊关系矩阵,通过模糊关系矩阵观察异常类的发展情况。
2.如权利要求1所述的时空联合数据聚类分析的输变电设备状态异常检测方法,其特征是,所述步骤S2:每个分块矩阵中的元素描述分块矩阵对应的时窗中存在的类。
3.如权利要求1所述的时空联合数据聚类分析的输变电设备状态异常检测方法,其特征是,所述步骤S3:如果异常度值越高,则说明设备状态存在异常情况的可能性越大。
4.如权利要求1所述的时空联合数据聚类分析的输变电设备状态异常检测方法,其特征是,所述步骤S4:每个时窗中都会进行聚类,将每个时窗中的时空子序列分为若干类,每个时窗对应着分块矩阵,每个时窗中的类对应着分块矩阵中的元素,利用分块矩阵建立不同时窗中的各个类之间的模糊关系,并获得相邻两个时窗之间的模糊关系矩阵R;
模糊关系矩阵R的第i行第j列的元素ri,j代表第一个时窗W1中的第i类与第二个时窗W2中的第j类的相关程度,
ri,j的值越大,则说明第一个时窗W1中的第i类与第二个时窗W2中的第j类的相关程度越大,
相关程度大说明第二个时窗W2中的第j类是由第一个时窗W1中的第i类发展而来的,由此得出一个时窗中的异常类与相邻时窗中对应的相关类,从而了解异常类的发展情况。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510943246.7A CN105574669B (zh) | 2015-12-16 | 2015-12-16 | 时空联合数据聚类分析的输变电设备状态异常检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510943246.7A CN105574669B (zh) | 2015-12-16 | 2015-12-16 | 时空联合数据聚类分析的输变电设备状态异常检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105574669A CN105574669A (zh) | 2016-05-11 |
CN105574669B true CN105574669B (zh) | 2020-02-14 |
Family
ID=55884772
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510943246.7A Active CN105574669B (zh) | 2015-12-16 | 2015-12-16 | 时空联合数据聚类分析的输变电设备状态异常检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105574669B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107248004B (zh) * | 2016-07-20 | 2020-02-14 | 国网山东省电力公司电力科学研究院 | 一种用于线路故障预测的时序数据粒度统一转换方法 |
CN107202852A (zh) * | 2017-05-23 | 2017-09-26 | 国家电网公司 | 一种基于可变阈值的油色谱在线监测数据异常值检测方法 |
CN107368556B (zh) * | 2017-07-04 | 2020-10-20 | 广西电网有限责任公司电力科学研究院 | 一种输电线路多源地理信息一致性匹配系统 |
CN110187242A (zh) * | 2019-07-02 | 2019-08-30 | 上海交通大学 | 一种基于时空相关性分析的配电网异常检测及定位方法 |
CN111612082B (zh) * | 2020-05-26 | 2023-06-23 | 河北小企鹅医疗科技有限公司 | 检测时间序列中异常子序列的方法及设备 |
CN112527788B (zh) * | 2020-12-17 | 2024-04-30 | 北京中恒博瑞数字电力科技有限公司 | 变压器监测数据异常值检测与清洗的方法及装置 |
CN114861729A (zh) * | 2022-05-20 | 2022-08-05 | 西安邮电大学 | 无线传感器网络中时间序列异常的检测方法及装置 |
CN116258281B (zh) * | 2023-05-12 | 2023-07-25 | 欣灵电气股份有限公司 | 基于云平台管理的物联网消防监测及调控系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101149760A (zh) * | 2007-11-08 | 2008-03-26 | 太原科技大学 | 一种天体光谱离群数据挖掘系统和方法 |
CN104700321A (zh) * | 2015-03-16 | 2015-06-10 | 国家电网公司 | 一种输变电设备状态运行趋势分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9500710B2 (en) * | 2012-10-15 | 2016-11-22 | Siemens Energy, Inc. | Generator neutral ground monitoring system and method |
-
2015
- 2015-12-16 CN CN201510943246.7A patent/CN105574669B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101149760A (zh) * | 2007-11-08 | 2008-03-26 | 太原科技大学 | 一种天体光谱离群数据挖掘系统和方法 |
CN104700321A (zh) * | 2015-03-16 | 2015-06-10 | 国家电网公司 | 一种输变电设备状态运行趋势分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105574669A (zh) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105574669B (zh) | 时空联合数据聚类分析的输变电设备状态异常检测方法 | |
CN110361176B (zh) | 一种基于多任务特征共享神经网络的智能故障诊断方法 | |
CN109446187B (zh) | 基于注意力机制与神经网络的复杂装备健康状态监测方法 | |
CN110087207B (zh) | 无线传感器网络缺失数据重建方法 | |
CN110619479B (zh) | 光伏组件的故障判别模型建模方法及故障判别方法 | |
Rao et al. | Dropout and pruned neural networks for fault classification in photovoltaic arrays | |
Gibson et al. | Evaluating synoptic systems in the CMIP5 climate models over the Australian region | |
CN112149887A (zh) | 一种基于数据时空特征的pm2.5浓度预测方法 | |
CN113988210A (zh) | 结构监测传感网失真数据修复方法、装置及存储介质 | |
CN112668526A (zh) | 基于深度学习和压电主动传感的螺栓群松动定位监测方法 | |
Kelder et al. | Interpreting extreme climate impacts from large ensemble simulations—are they unseen or unrealistic? | |
Moezi et al. | Simultaneous fault localization and detection of analog circuits using deep learning approach | |
García et al. | A Bayesian hierarchical spatio-temporal model for extreme rainfall in Extremadura (Spain) | |
Mahajan et al. | Incremental outlier detection in air quality data using statistical methods | |
CN117791856B (zh) | 一种基于巡检机器人的电网故障预警方法及装置 | |
Eo et al. | Deep learning framework with essential pre-processing techniques for improving mixed-gas concentration prediction | |
CN111931738B (zh) | 用于遥感影像的神经网络模型预训练方法、装置 | |
Richards et al. | A unifying partially-interpretable framework for neural network-based extreme quantile regression | |
Li et al. | Meteorological radar fault diagnosis based on deep learning | |
CN116756503A (zh) | 一种结构健康监测地震事件识别与异常数据诊断方法 | |
CN114598627A (zh) | 一种基于知识图谱的异常网络信息检测方法 | |
Hassan et al. | Fault classification of power plants using artificial neural network | |
CN116739996A (zh) | 基于深度学习的输电线路绝缘子故障诊断方法 | |
CN114282608B (zh) | 一种电流互感器隐性故障诊断预警方法及系统 | |
CN115510763A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |