CN112699913B - 一种台区户变关系异常诊断方法及装置 - Google Patents

一种台区户变关系异常诊断方法及装置 Download PDF

Info

Publication number
CN112699913B
CN112699913B CN202011343787.3A CN202011343787A CN112699913B CN 112699913 B CN112699913 B CN 112699913B CN 202011343787 A CN202011343787 A CN 202011343787A CN 112699913 B CN112699913 B CN 112699913B
Authority
CN
China
Prior art keywords
abnormal
area
clustering
user
voltage data
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
CN202011343787.3A
Other languages
English (en)
Other versions
CN112699913A (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.)
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Metering Center of State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Metering Center of State Grid Hunan Electric Power Co Ltd
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 State Grid Corp of China SGCC, State Grid Hunan Electric Power Co Ltd, Metering Center of State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202011343787.3A priority Critical patent/CN112699913B/zh
Priority to PCT/CN2021/077141 priority patent/WO2022110557A1/zh
Publication of CN112699913A publication Critical patent/CN112699913A/zh
Application granted granted Critical
Publication of CN112699913B publication Critical patent/CN112699913B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/23—Clustering techniques
    • G06F18/232—Non-hierarchical techniques
    • G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
    • G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/04—Testing or calibrating of apparatus covered by the other groups of this subclass of instruments for measuring time integral of power or current
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/08—Probabilistic or stochastic CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Probability & Statistics with Applications (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Remote Monitoring And Control Of Power-Distribution Networks (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开一种台区户变关系异常诊断方法及装置，该方法步骤包括：S01.获取指定台区的总表和用户电表电压数据并进行降维处理，得到降维电压数据；S02.对降维电压数据采用K‑means聚类方法进行聚类，并在聚类过程中根据数据维度的最大和最小值选取初始聚类中心，查找出指定台区内户变关系存在异常的异常用户；S03.分别计算异常用户与各个台区总表之间的皮尔逊相关系数，将计算出的各所述皮尔逊相关系数进行排序，根据排序结果诊断出异常用户所属的正确台区。本发明能够实现同一台区一个及多个异常用户、不同台区多个异常用户的诊断，且具有实现方法简单、诊断效率以及精度高等优点。

Description

一种台区户变关系异常诊断方法及装置
技术领域
本发明涉及用电信息采集系统技术领域，尤其涉及一种台区户变关系异常诊断方法及装置。
背景技术
随着电网建设的快速发展，包含海量数据的用电信息采集系统不仅可直接反映配电网的运行状况，且能间接反映配电网的拓扑关系。现有配电网拓扑结构校验主要包括：线变关系校验、馈线拓扑校验、户变关系及相位校验、线户关系校验。正确的低压配电网拓扑，尤其正确的户变关系是当前配电网管理精细化和降耗减损的基础。然而用电信息采集系统易出现台区户变关系错误问题，存在部分台区线路临时改变使户变关系档案更新不及时或记录错误等原因，导致用户进线端和集中器归属关系记录不准确，户变关系与实际不符，造成台区线损计算时出现异常情况。目前低压台区户变关系错误主要形式为主站系统中记录的用户档案与实际信息不相符，且通常把某个台区用户错误挂接至相邻台区，或把某个馈线变压器错误挂接至相邻馈线，该情况会影响正确的户变关系，阻碍正常线损计算。若将实际属于台区A的用户1档案信息，错归为台区B所有，则线损计算时将造成台区A线损偏小，而台区B线损偏大结果发生。因此，台区户变关系的准确诊断与分析对当前配电网的精益化管理具有重要意义。
现有技术中校验台区户变关系主要有人工专用设备现场识别和在线自动识别两种方式，其中人工方式需安排工作人员至现场识别，识别效率低、成本高，在线自动识别可以解决人工方式存在的上述问题，目前在线自动识别方法主要是基于数据挖掘技术，通过提取台区异常用户特征来识别特征相符的台区，如通过提取用户电压曲线相似性特征实现台区户变关系异常诊，在低压配电网中由于各处负荷的不确定性，电压通常也随之波动，电气距离较近负荷，其电压波动曲线较相似，而电气距离较远负荷其电压波动曲线相似度较低，因此可选取用户电压曲线相似性作为台区户变关系异常诊断的依据。但是上述在线自动识别方法需要依赖于异常用户的特征提取，仅能针对台区中出现少量异常用户的情况，对于多达数百用户台区，若存在多个相邻台区存在数量较多的异常用户时，大量的异常用户特征提取困难，且提取的特征所能够表征的信息有限，就难以快速、准确的实现多用户台区和多个异常用户情况的诊断与分析。
发明内容
本发明要解决的技术问题就在于：针对现有技术存在的技术问题，本发明提供一种实现方法简单、诊断效率以及精度高的台区户变关系异常诊断方法及装置，能够实现同一台区一个及多个异常用户、不同台区多个异常用户的诊断。
为解决上述技术问题，本发明提出的技术方案为：
一种台区户变关系异常诊断方法，步骤包括：
S01.降维处理：获取指定台区的总表和用户电表电压数据并进行降维处理，得到降维电压数据；
S02.聚类处理：对所述降维电压数据采用K-means聚类方法进行聚类，并在聚类过程中根据数据维度的最大和最小值选取初始聚类中心，查找出指定台区内户变关系存在异常的异常用户；
S03.台区诊断：分别计算所述异常用户与各个台区总表之间的皮尔逊相关系数，将计算出的各所述皮尔逊相关系数进行排序，根据排序结果诊断出所述异常用户所属的正确台区。
进一步的，所述步骤S02中选取初始聚类中心的步骤包括：
S201.使用所述降维电压数据中每一维最大和最小值选取k个样本点作为k个初始聚类中心；
S202.分别计算所述降维电压数据中每个样本数据距离最近聚类中心的距离D(zij)，即：
其中，zij为所述降维电压数据；
S203.取所述降维电压数据中每个样本数据距离聚类中心的最大距离作为下一个聚类中心，以更新各簇的聚类中心，即：
μr＝argmax|D(zij)|
其中，μr为各簇更新的聚类中心，max|D(zij)|为每个样本与当前最近一个聚类中心的最大距离；
S204.重复执行步骤S202、步骤S203，直至聚类误差函数E收敛或达到最大步数，得到最终的k个聚类中心[μ1,μ2…,μk]。
进一步的，所述步骤S201中具体按照下式选取k个初始聚类中心；
μr＝ma(i)+(mi(i)-ma(i))+rand(),r＝1,2,...,k
其中，μr为初始聚类中心，ma为每一维最大数，mi为每一维最小数。
进一步的，所述步骤S03中具体取排序结果中最大皮尔逊相关系数所对应的台区作为所述异常用户所属的正确台区。
进一步的，所述步骤S03中还包括先分别计算查找出的所有异常电表的电压数据之间的相关性值，将相关性值小于预设阈值的异常电表判定为属于同一台区，从判定为同一台区的各异常电表中选取出一个作为待检测电表，然后分别计算待检测电表与各个台区总表之间的皮尔逊相关系数以诊断出所属的正确台区。
进一步的，具体对查找出的所有异常电表先进行聚类处理，形成多个聚类，获取作为每个聚类的聚类中心的异常电表作为待检测电表，然后分别计算待检测电表与各个台区总表之间的皮尔逊相关系数以诊断出目标异常电表所属的正确台区，由诊断出的目标异常电表所属的正确台区得到对应的整个聚类中各异常电表所属的正确台区。
进一步的，所述步骤S03诊断出所述异常用户所属的正确台区后，还包括以诊断出的所述正确台区作为目的台区，获取目的台区的用户电压数据并返回步骤S02以重新进行聚类处理，根据聚类结果最终确认异常用户是否属于目的台区。
进一步的：所述步骤S01中具体采用主成分分析方法进行降维处理。
一种台区户变关系异常诊断装置，包括：
降维处理模块，用于获取指定台区的各用户电压数据并进行降维处理，得到降维电压数据；
聚类处理模块，用于对所述降维电压数据采用K-means聚类方法进行聚类，并在聚类过程中根据数据维度的最大和最小值选取初始聚类中心，查找出指定台区内户变关系存在异常的异常用户作为待检测用户；
台区诊断模块，用于分别计算所述待检测用户与各个台区总表之间的皮尔逊相关系数，将计算出的各所述皮尔逊相关系数进行排序，根据排序结果诊断出所述异常用户所属的正确台区。
一种台区户变关系异常诊断装置，包括处理器以及存储器，所述存储器用于存储计算机程序，所述处理器用于执行所述计算机程序，所述处理器用于执行所述计算机程序，以执行如上述方法。
与现有技术相比，本发明的优点在于：
1、本发明通过获取台区总表和用户电表电压数据降维后，使用改进的K-means聚类提取电压数据特征，在聚类时根据数据维度的最大和最小值选取初始聚类中心，使得可以快速准确的查找出台区内所有的异常用户，同时对查找出的异常用户，使用改进皮尔逊相关系数方法，由异常用户与各个台区总表之间的皮尔逊相关系数之间的大小关系诊断出正确的台区，能够准确提取多个相邻台区异常用户间特征，实现多异常用户所属正确台区诊断，从而在同一台区一个及多个异常用户、不同台区多个异常用户情况下均能有效实现异常用户的准确检测与分析。
2、本发明能够高效实现电压数据间的特征信息挖掘，通过在使用K-means聚类方法提取电压数据特征时，根据数据维度的最大最小值选取更优聚类中心，使得可以得到不同特征异常用户，同时通过使用皮尔逊相关系数表征异常用户与台区之间相关程度的基础上，通过取待检测用户与各个台区总表间的最大相关系数作为所求的正确台区，可以有效避免皮尔逊相关系数参考值的选取，从而提高检测效率以及精度，尤其适用于同一个台区存在多台异常用户、相邻台区存在多台异常用户时的户变关系自动诊断。
3、本发明进一步先将异常台区按照相关性进行划分，将相关性值小于预设阈值的异常电表判定为属于同一台区，从判定为同一台区的各异常电表中再选取出一个作为待检测电表，然后分别计算待检测电表与各个台区总表之间的皮尔逊相关系数以诊断出所属的正确台区，使得仅需要少量的计算即可实现所有异常电表的台区诊断，无需一一对所有异常电表进行诊断，可以进一步提高诊断的效率，避免大量不必要的计算，尤其是在存在大规模的电表划分异常时，能够快速的诊断出所有电表所属的正确台区。
附图说明
图1是本实施例台区户变关系异常诊断方法的实现流程示意图。
图2是本实施例中电压数据聚类中心选取的原理示意图。
图3是在具体应用实施例中得到的仿真结果示意图。
图4是在具体应用实施例中得到的数据长度影响分析结果示意图。
具体实施方式
以下结合说明书附图和具体优选的实施例对本发明作进一步描述，但并不因此而限制本发明的保护范围。
如图1所示，本实施例台区户变关系异常诊断方法的步骤包括：
S01.降维处理：获取指定台区的总表和用户电表电压数据并进行降维处理，得到降维电压数据；
S02.聚类处理：对降维电压数据采用K-means聚类方法进行聚类，并在聚类过程中根据数据维度的最大和最小值选取初始聚类中心，查找出指定台区内户变关系存在异常的异常用户；
S03.台区诊断：分别计算异常用户与各个台区总表之间的皮尔逊相关系数，将计算出的各皮尔逊相关系数进行排序，根据排序结果诊断出异常用户所属的正确台区。
本实施例通过获取台区总表和用户电表电压数据降维后，使用改进的K-means聚类提取电压数据特征，在聚类时根据数据维度的最大和最小值选取初始聚类中心，而不是如传统的随机选取初始聚类中心，可以快速准确的查找出台区内所有的异常用户，同时对查找出的异常用户，使用改进皮尔逊相关系数方法，由异常用户与各个台区总表之间的皮尔逊相关系数之间的大小关系诊断出正确的台区，能够准确提取多个相邻台区异常用户间特征，实现多异常用户所属正确台区诊断，从而在同一台区一个及多个异常用户、不同台区多个异常用户情况下均能有效实现异常用户的准确检测与分析。
本实施例具体使用主成分分析PCA方法对原始电压数据进行降维处理。通常台区用户电表电压数据为一天24小时，每隔1小时采集一次，对于数据维度为24维的高维数据，传统聚类算法将面临高维数据包含海量冗余、不相干信息问题，直接对高维数据聚类极大降低性能，聚类算法难以实现高稳定性，本实施例通过采用PCA方法对原始数据进行降维处理，可以便于后续对电压数据进行聚类处理而查找出异常用户。
本实施例使用主成分分析PCA对原始电压数据进行降维处理的详细步骤为：
步骤S101.基于PCA将p维电压数据特征在m维上体现，该m维数据信息也被称为主成分，为在原有p维电压数据特征上构造出m维新特征，将台区n个用户电表p维电压数据表示为：
通过PCA得出的主成分之间互不相关，计算电压数据相关系数具体为：
式中，xi、xj为台区电压数据，为电压数据均值。
由式(2)可得电压数据相关系数矩阵为：
式中，rij(i,j＝1,2,…,p)为台区电压数据xi、xj的相关系数，rij＝rji。
根据式(3)解特征方程得到：
|λI-R|＝0 (4)
根据数值大小得特征值λ1≥λ2≥…≥λp≥0，同理可得对应于特征值λi的特征向量ei(i＝1,2,…,p)。
步骤S102.根据主成分贡献率法(CPV)确定主成分个数m，计算各主成分的贡献率CPVi与累计贡献率CPVa分别为：
式中，CPVi为第i个主成分的贡献率，λm为第m(m≤p)个主成分所对应的特征值，CPVa为控制限。
通过式(5)和(6)计算包含原始电压数据绝大部分信息的主成分个数m，使用式(1)中原始电压数据矩阵X的m个特征向量作线性组合得到主成分，则有：
式中，zij(i＝1,2,…,n,j＝1,2,…,m)为台区用户电表电压数据第i个样本第j个主成分。
相同台区用户电压数据波动相似性使其具有相同特征，本实施例基于该特性使用K-means聚类检测出异常用户。K-means聚类是根据迭代理论，求出最大迭代步数或使得聚类误差函数收敛得到聚类中心，该聚类方法为无监督学习算法，可以适用于检测台区中异常用户，但是传统的K-means聚类方法中是采用随机选取初始聚类中心的方式，若初始聚类中心选取不合适，会极大影响最终聚类结果。本实施例对传统的K-means聚类进行改进，通过根据数据维度的最大和最小值，基于尽可能得到相互距离远的初始聚类中心原则选取初始聚类中心，而不是如传统方法随机选取初始聚类中心，可以避免出现初始聚类中心选择错误或不合适的问题，从而大大提高异常用户的检测精度以及效率。
本实施例步骤S02中选取初始聚类中心的具体步骤包括：
S201.使用降维电压数据中每一维最大和最小值选取k个样本点作为k个初始聚类中心，具体按照下式(8)选取k个初始聚类中心；
μr＝ma(i)+(mi(i)-ma(i))+rand(),r＝1,2,...,k (8)
其中，μr为初始聚类中心，ma为每一维最大数，mi为每一维最小数；
S202.分别计算降维电压数据中每个样本数据距离最近聚类中心的距离D(zij)，即：
其中，zij为降维电压数据；
S203.取降维电压数据中每个样本数据距离聚类中心的最大距离作为下一个聚类中心，以更新各簇的聚类中心，即：
μr＝argmax|D(zij)| (10)
其中，μr为各簇更新的聚类中心，max|D(zij)|为每个样本与当前最近一个聚类中心的最大距离；
S204.重复执行步骤S202、步骤S203，直至聚类误差函数E收敛或达到最大步数，得到最终的k个聚类中心[μ1,μ2…,μk]。
本实施例先对于降维电压数据选取k个样本点作为k个初始聚类中心，针对电压中每个样本数据，分别计算其到已有最近聚类中心欧式距离，由式(9)所求距离值越大表示该样本被选取作为下一次聚类中心的概率越大，样本数据被分别归为最近聚类中心类别后得到k个簇，计算样本被选取作为聚类中心的概率可得：
式中，D(zij)为样本到聚类中心的距离。
在具体应用实施例中对台区用户电表电压数据通过PCA降维后，得到2维数据为A(0.1,0.1)，B(0.2,0.2)，C(0.2,0)，D(0.4,0.6)，E(0.5,0.6)，F(0.5,0.5)，G(0.6,0.5)，再对降维后电压数据进行聚类分析，并在首次聚类中选择C1(0.2,0.3)，C2(0.5,0.3)两点作为第1个和第2个聚类中心如图2所示。由图2可知，计算每一簇中其他电压数据样本与该簇当前已有最近聚类中心距离D(zij)，以及各簇中每个电压数据样本被选取作为下一个聚类中心的概率P(zij)分别如表1和2所示。
表1中心点1聚类距离与概率
表2中心点2聚类距离与概率
由表1可见，对于第一簇数据A点被选为下一个聚类中心的概率最大，对于第二簇数据D点被选为下一个聚类中心的概率最大。而由图2可知，A，D分别为距离初始聚类中心C1和C2点最远的两个点。
由式(9)求得每个样本与当前最近一个聚类中心的距离并选取最大值，基于相互距离尽可能远原则选取该样本为该簇新聚类中心点，则按照式(10)即可计算出最佳的新聚类中心点作为各簇下一个聚类中心；重复上述式(9)和(11)不断移动聚类中心直至聚类误差函数收敛或达到最大迭代步数，平方误差SSE函数为：
根据式(12)直至平方误差SSE收敛或达到最大步数，选出k个聚类中心[μ1,μ2…,μk]，并分别实现以μr为聚类中心的台区用户电表电压数据聚类求得台区内异常用户。
本实施例通过按照上述改进的K-means聚类分析步骤对相邻各台区内所有的用户电表电压数据进行聚类分析，即可查找出相邻各台区内所有的用户电表电压数据异常的用户。由于皮尔逊相关系数可以表征两个变量间的线性相关性，本实施例在查找出异常用户后，通过使用皮尔逊相关系数来判断异常用户与各台区总表之间用户电表电压数据的相关程度，以诊断出待检测异常用户所属的正确台区。
皮尔逊相关系数用来度量两个电压数据X与Y之间的相互关系，可度量两个变量线性相关的强弱，其总体相关系数为：
式中，μX，μY为期望值，σX、σY为总体标准差，cov(X,Y)为协方差。
式(13)对于样本的皮尔逊相关系数同样可用，计算电压样本相关系数为：
式中，为标准化变量；/>为样本均值；Sx、Sy为样本标准差。
皮尔逊相关系数因两个变量的位置和尺度的变化并不会引起该系数的改变，即把X移动到a+bX和把Y移动到c+dY，其中a、b、c和d是常数，对两个变量间相关系数毫无影响，且对于总体以及样本皮尔逊相关系数公式同时成立，由此可得出更一般的线性变换则会改变相关系数。
由于台区电表电压数据X的期望μX与方差σ2 X分别为：
同理，计算台区电表电压数据Y的期望μY与方差σ2 Y分别为：
根据两个电压数据X与Y间期望变换公式为：
E[(X-E(X))(Y-E(Y))]＝E(XY)-E(X)E(Y) (16)
可计算总体相关系数为
式中，ρX,Y为总体皮尔逊相关系数，E(X)、E(Y)为电压数据X、Y的期望值。
同理，由式(19)计算电压数据样本皮尔逊相关系数可得：
式中，rx,y为样本皮尔逊相关系数，Zx、Zy为标准化变量，Sx、Sy为X、Y样本标准差， 表示样本均值，n为样本数。
皮尔逊相关系数为|rx,y|≤1，rx,y大于零表示两者为正相关方向，小于零表示为负相关方向，评价皮尔逊相关系数的标准如表3所示。
表3皮尔逊相关系数关联度标准
由表3可见，若直接使用皮尔逊相关系数表征两个电压数据X与Y之间相关程度，皮尔逊相关系数的阈值取值会决定相关程度的判定结果，但实际应用中就难以准确的设定参考值以判断某用户是否属于某台区，且当通过聚类查找到台区内的异常用户后，由于该台区及邻近几个台区用户较多，若计算待检测用户与各个台区间之间的皮尔逊相关系数，即计算异常用户与台区中所有电表数据之间的皮尔逊相关系数，还会存在工作量巨大的问题。
本实施例在通过聚类查找到台区内的异常用户后，通过计算待检测用户与各台区总表之间的皮尔逊相关系数，再对计算出的各皮尔逊相关系数进行排序，依据排序结果来诊断异常用户所属的正确台区，无需计算异常用户与台区所有电压数据之间的皮尔逊相关系数，可以大大减少计算量，提高诊断效率，同时无需设置诊断的皮尔逊相关系数参考值，只需对待检测用户与各个台区总表间的皮尔逊相关系数进行大小排序，有效避免皮尔逊相关系数参考值的选取，可以充分利用相邻台区之间的皮尔逊相关系数，大大提高最终的诊断精度、降低诊断的复杂度。
本实施例步骤S03中具体取排序结果中最大皮尔逊相关系数所对应的台区作为异常用户所属的正确台区，即相关系数最大者对应台区为待检测用户所属正确台区。
本实施例上述方法，能够高效实现电压数据间的特征信息挖掘，通过在使用K-means聚类方法提取电压数据特征时，根据数据维度的最大最小值选取更优聚类中心，使得可以得到不同特征异常用户，同时通过使用皮尔逊相关系数表征异常用户与台区之间相关程度的基础上，通过取待检测用户与各个台区总表间的最大相关系数作为所求的正确台区，可以有效避免皮尔逊相关系数参考值的选取，从而提高检测效率以及精度，尤其适用于同一个台区存在多台异常用户、相邻台区存在多台异常用户时的户变关系自动诊断。
本实施例步骤S03中还包括先分别计算查找出的所有异常电表的电压数据之间的相关性值，将相关性值小于预设阈值的异常电表判定为属于同一台区，从判定为同一台区的各异常电表中选取出一个作为待检测电表，然后分别计算待检测电表与各个台区总表之间的皮尔逊相关系数以诊断出所属的正确台区。
当相邻台区中存在大量的异常用户时，若一一分别计算异常用户与各台区总表的相关系数进行台区诊断，依然会存在需要进行大量计算的问题。本实施例进一步先将异常台区按照相关性进行划分，将相关性值小于预设阈值的异常电表判定为属于同一台区，从判定为同一台区的各异常电表中再选取出一个作为待检测电表，然后分别计算待检测电表与各个台区总表之间的皮尔逊相关系数以诊断出所属的正确台区，使得仅需要少量的计算即可实现所有异常电表的台区诊断，无需一一对所有异常电表进行诊断，可以进一步提高诊断的效率，避免大量不必要的计算，尤其是在存在大规模的电表划分异常时，能够快速的诊断出所有电表所属的正确台区。
本实施例具体对查找出的所有异常电表先进行聚类处理，形成多个聚类，获取作为每个聚类的聚类中心的目标异常电表，然后分别计算目标异常电表与各个台区总表之间的皮尔逊相关系数以诊断出目标异常电表所属的正确台区，由诊断出的目标异常电表所属的正确台区得到对应的整个聚类中各异常电表所属的正确台区。即采用聚类的方式先将异常电表进行分类，后续仅需要计算各聚类中心与台区总表之间的皮尔逊相关系数，利用聚类的方式快速、精准的对所有异常电表进行诊断。
本实施例步骤S03诊断出异常用户所属的正确台区后，还包括以诊断出的正确台区作为目的台区，获取目的台区的用户电压数据并返回步骤S02以重新进行聚类处理，根据聚类结果最终确认异常用户是否属于目的台区。按照上述步骤诊断出的台区，也可能会存在诊断误差，如若两个台区的电压特性非常接近，所得到的皮尔逊相关系数也较为接近，直接依据皮尔逊相关系数来确定正确台区即可能会存在误差。本实施例在诊断出异常用户的台区后，进一步获取该台区的用户电压数据重新进行聚类处理，若诊断的台区正确则再次聚类时可以正确分类，因而可以依据再次聚类处理的结果验证所诊断的台区是否正确，进一步提高诊断的精度。
为验证本发明上述按照改进的K-means聚类分析查找出异常用户、以及使用最大皮尔逊相关系数诊断出正确台区的有效性，在具体实施例中使用数据为某个台区其中137个用户的每24个小时每隔1小时采集的电压数据进行仿真。先判断出异常用户，对于需校验用户，获取相邻台区总表电压数据，诊断待校验用户所属正确台区，分别采用以下各种情况对本发明上述方法进行验证，结果如图3所示。
(1)1台区1用户分析
首先针对该台区相邻1个台区中的1个用户错接入该台区的情况进行仿真和分析，将2台区1用户加入1台区中，并对该用户进行异常诊断，找出该用户所属正确台区，计算待检测用户与5个台区总表间的皮尔逊相关系数如图3(a)所示。由图3(a)可见，该用户与2台区总表间的皮尔逊相关系数为0.9910，在5个台区中为最大值，本实施例取相关系数最大者所对应台区为待检测用户所属的正确台区，即该用户属于2台区，与真实结果相同，由此验证本发明上述方法在相邻1个台区中1个用户错接入该台区情况下的准确性和可行性。
为验证本发明方法在1个台区多个用户错接入其他台区情况下的准确性和可行性，针对该台区相邻1个台区中的7个用户错接入该台区的情况进行仿真和分析，将3台区7个用户加入1台区中，并进行诊断分析，计算待检测7个用户与5个台区总表间的皮尔逊相关系数如图3(b)所示。计算待检测7个用户与5个台区总表间的皮尔逊相关系数，并将皮尔逊相关系数所对应台区与实际台区作比较结果如表4所示。
表4皮尔逊相关系数测试结果
由表4可见，可直观得出7个待检测用户与台区总表间最大皮尔逊相关系数对应的台区均为3台区，实际对应台区也均属于3台区，由此可见本发明方法在相邻1个台区中n个用户错接入该台区情况下的具有较高的准确度性。
为验证本发明上述方法在多个台区用户错接入其他台区情况下的准确性和可行性，针对该台区相邻2个台区中的各1个用户错接入该台区的情况进行仿真和分析，分别将2台区1个用户和3台区1个用户加入1台区，并进行诊断分析，计算待检测2个用户与5个台区总表间的皮尔逊相关系数如图3(c)所示。由图3(c)可知，待检测用户1仅与2台区总表间的皮尔逊相关系数超过0.9900，待检测用户2仅与3台区总表间的皮尔逊相关系数超过0.9800，根据相关系数最大者所对应的台区即为待检测用户所属的正确台区，表明2个待检测用户分别属于2台区、3台区，与真实结果相同，由此可验证本发明上述方法在相邻n个台区中1个用户错接入该台区情况下的准确性和有效性。
为验证本发明方法在多个台区多个用户错接入其他台区情况下的准确性和可行性，针对该台区相邻2个台区中各7个用户错接入该台区的情况进行仿真和分析，分别将2台区7个用户和3台区7个用户加入1台区，并进行诊断分析，计算待检测14个用户与5个台区总表之间的皮尔逊相关系数如图3(d)所示。由图3(d)可见，用户1到用户7与台区3皮尔逊相关系数最大，用户8、用户10到用户14与台区2皮尔逊相关系数最大，用户9与台区1皮尔逊相关系数最大。因此除用户9外，1到7个待检测用户均属于3台区，8到14个待检测用户均属于2台区。表明除用户9外，其余待检测用户诊断结果与真实结果相同，又因聚类结果中用户9已被识别为1台区异常用户，故从1台区排除，且用户9与台区2总表皮尔逊相关系数大于其他台区，故综合判断可得出用户9所属正确台区。由此可见本发明在相邻n个台区中n个用户错接入该台区情况下仍有较高的准确性。
为分析用户电能表电压数据长度对识别结果的影响，本实施例提出异常用户识别正确率指标，以反映台区户变关系识别准确性，该指标定义为诊断出异常用户数与待诊断总用户数比值，数值越大表明识别结果越准确，电压数据长度定义为用户电压时间维度，选取一天中2个时刻的电压值到24个时刻的电压值情况下，本发明上述方法的识别正确率如图4所示。由图4可见，当电压数据长度低于10维时，本发明方法识别正确率随着电压数据长度的增加，呈现先急剧上升后平缓上升的趋势，并在11维逐渐收敛，稳定在90％以上的正确率，表明电压数据长度越大，本发明方法对户变关系异常诊断结果越准确。
为比较不同户变关系异常诊断的方法，选取常用的改进灰色关联度分析法、BP神经网络法以及本发明上述方法进行仿真比较。以多个台区多个用户错接入其他台区为例，将该台区相邻2个台区中的各7个用户放入该台区进行诊断识别，其仿真对比如表5所示。
表5与常用算法诊断结果对比
由表5可见，3种检测方法异常诊断结果不同，其中改进灰色关联度分析法异常用户诊断数为7，正确识别率为50％，BP神经网络法识别异常用户数为9，正确识别率为64.29％，本文算法识别异常用户数为13，正确识别率为92.86％，综合判断可得到另外一个异常用户所属正确台区。即本发明上述通过改进K-means聚类和皮尔逊相关系数进行户变关系异常诊断的方法，相比常用检测方法具有更高的诊断准确率。
本实施例还提供台区户变关系异常诊断装置，该装置包括：
降维处理模块，用于获取指定台区的各用户电压数据并进行降维处理，得到降维电压数据；
聚类处理模块，用于对降维电压数据采用K-means聚类方法进行聚类，并在聚类过程中根据数据维度的最大和最小值选取初始聚类中心，查找出指定台区内户变关系存在异常的异常用户作为待检测用户；
台区诊断模块，用于分别计算待检测用户与各个台区总表之间的皮尔逊相关系数，将计算出的各皮尔逊相关系数进行排序，根据排序结果诊断出异常用户所属的正确台区。
本实施例中，台区诊断模块还包括先分别计算查找出的所有异常电表的电压数据之间的相关性值，将相关性值小于预设阈值的异常电表判定为属于同一台区，从判定为同一台区的各异常电表中选取出一个作为待检测电表，然后分别计算待检测电表与各个台区总表之间的皮尔逊相关系数以诊断出所属的正确台区。
本实施例中，台区诊断模块具体对查找出的所有异常电表先进行聚类处理，形成多个聚类，获取作为每个聚类的聚类中心的异常电表作为待检测电表，然后分别计算待检测电表与各个台区总表之间的皮尔逊相关系数以诊断出目标异常电表所属的正确台区，由诊断出的目标异常电表所属的正确台区得到对应的整个聚类中各异常电表所属的正确台区。
本实施例中，台区诊断模块诊断出异常用户所属的正确台区后，还包括以诊断出的正确台区作为目的台区，获取目的台区的用户电压数据并返回聚类处理模块以重新进行聚类处理，根据聚类结果最终确认异常用户是否属于目的台区。
本实施例台区户变关系异常诊断装置与上述台区户变关系异常诊断方法为一一对应，在此不再一一赘述。
在另一实施例中，本发明台区户变关系异常诊断装置还可以为包括处理器以及存储器，存储器用于存储计算机程序，处理器用于执行计算机程序，其特征在于，处理器用于执行计算机程序，以执行如上述台区户变关系异常诊断方法。
上述只是本发明的较佳实施例，并非对本发明作任何形式上的限制。虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明。因此，凡是未脱离本发明技术方案的内容，依据本发明技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均应落在本发明技术方案保护的范围内。

Claims (9)

1.一种台区户变关系异常诊断方法，其特征在于，步骤包括：
S01.降维处理：获取指定台区的总表和用户电表电压数据并进行降维处理，得到降维电压数据；
S02.聚类处理：对所述降维电压数据采用K-means聚类方法进行聚类，并在聚类过程中根据数据维度的最大和最小值选取初始聚类中心，查找出指定台区内户变关系存在异常的异常用户；
S03.台区诊断：分别计算所述异常用户与各个台区总表之间的皮尔逊相关系数，将计算出的各所述皮尔逊相关系数进行排序，根据排序结果诊断出所述异常用户所属的正确台区；
所述步骤S03中还包括先分别计算查找出的所有异常电表的电压数据之间的相关性值，将相关性值小于预设阈值的异常电表判定为属于同一台区，从判定为同一台区的各异常电表中选取出一个作为待检测电表，然后分别计算待检测电表与各个台区总表之间的皮尔逊相关系数以诊断出所属的正确台区。
2.根据权利要求1所述的台区户变关系异常诊断方法，其特征在于，所述步骤S02中选取初始聚类中心的步骤包括：
S201.使用所述降维电压数据中每一维最大和最小值选取k个样本点作为k个初始聚类中心；
S202.分别计算所述降维电压数据中每个样本数据距离最近聚类中心的距离D(zij)，即：
其中，zij为所述降维电压数据；
S203.取所述降维电压数据中每个样本数据距离聚类中心的最大距离作为下一个聚类中心，以更新各簇的聚类中心，即：
μr＝argmaxD(zij)
其中，μr为各簇更新的聚类中心，max|D(zij)|为每个样本与当前最近一个聚类中心的最大距离；
S204.重复执行步骤S202、步骤S203，直至聚类误差函数E收敛或达到最大步数，得到最终的k个聚类中心[μ1,μ2…,μk]。
3.根据权利要求2所述的台区户变关系异常诊断方法，其特征在于，所述步骤S201中具体按照下式选取k个初始聚类中心；
μr＝ma(i)+(mi(i)-ma(i))+rand(),r＝1,2,...,k
其中，μr为初始聚类中心，ma为每一维最大数，mi为每一维最小数。
4.根据权利要求1或2或3所述的台区户变关系异常诊断方法，其特征在于，所述步骤S03中具体取排序结果中最大皮尔逊相关系数所对应的台区作为所述异常用户所属的正确台区。
5.根据权利要求1所述的台区户变关系异常诊断方法，其特征在于，具体对查找出的所有异常电表先进行聚类处理，形成多个聚类，获取作为每个聚类的聚类中心的异常电表作为待检测电表，然后分别计算待检测电表与各个台区总表之间的皮尔逊相关系数以诊断出目标异常电表所属的正确台区，由诊断出的目标异常电表所属的正确台区得到对应的整个聚类中各异常电表所属的正确台区。
6.根据权利要求1或2或3所述的台区户变关系异常诊断方法，其特征在于，所述步骤S03诊断出所述异常用户所属的正确台区后，还包括以诊断出的所述正确台区作为目的台区，获取目的台区的用户电压数据并返回步骤S02以重新进行聚类处理，根据聚类结果最终确认异常用户是否属于目的台区。
7.根据权利要求1或2或3所述的台区户变关系异常诊断方法，其特征在于，所述步骤S01中具体采用主成分分析方法进行降维处理。
8.一种台区户变关系异常诊断装置，其特征在于，包括：
降维处理模块，用于获取指定台区的各用户电压数据并进行降维处理，得到降维电压数据；
聚类处理模块，用于对所述降维电压数据采用K-means聚类方法进行聚类，并在聚类过程中根据数据维度的最大和最小值选取初始聚类中心，查找出指定台区内户变关系存在异常的异常用户作为待检测用户；
台区诊断模块，用于分别计算所述待检测用户与各个台区总表之间的皮尔逊相关系数，将计算出的各所述皮尔逊相关系数进行排序，根据排序结果诊断出所述异常用户所属的正确台区；
所述台区诊断模块中还包括先分别计算查找出的所有异常电表的电压数据之间的相关性值，将相关性值小于预设阈值的异常电表判定为属于同一台区，从判定为同一台区的各异常电表中选取出一个作为待检测电表，然后分别计算待检测电表与各个台区总表之间的皮尔逊相关系数以诊断出所属的正确台区。
9.一种台区户变关系异常诊断装置，其特征在于，包括处理器以及存储器，所述存储器用于存储计算机程序，其特征在于，所述处理器用于执行所述计算机程序，以执行如权利要求1～7中任意一项所述方法。
CN202011343787.3A 2020-11-25 2020-11-25 一种台区户变关系异常诊断方法及装置 Active CN112699913B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202011343787.3A CN112699913B (zh) 2020-11-25 2020-11-25 一种台区户变关系异常诊断方法及装置
PCT/CN2021/077141 WO2022110557A1 (zh) 2020-11-25 2021-02-22 一种台区户变关系异常诊断方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011343787.3A CN112699913B (zh) 2020-11-25 2020-11-25 一种台区户变关系异常诊断方法及装置

Publications (2)

Publication Number Publication Date
CN112699913A CN112699913A (zh) 2021-04-23
CN112699913B true CN112699913B (zh) 2023-08-29

Family

ID=75506037

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011343787.3A Active CN112699913B (zh) 2020-11-25 2020-11-25 一种台区户变关系异常诊断方法及装置

Country Status (2)

Country Link
CN (1) CN112699913B (zh)
WO (1) WO2022110557A1 (zh)

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113297539B (zh) * 2021-05-07 2024-10-18 天云融创数据科技（北京）有限公司 户变关系异常识别方法、装置、电子设备和介质
CN113657171A (zh) * 2021-07-20 2021-11-16 国网上海市电力公司 基于图小波神经网络的低压配电网台区拓扑识别方法
CN113804965B (zh) * 2021-09-28 2024-05-31 北京来也网络科技有限公司 基于rpa和ai的异常计量点识别方法和装置
CN113744089B (zh) * 2021-11-08 2022-02-15 广东电网有限责任公司佛山供电局 一种台区户变关系辨识方法及设备
CN114152909A (zh) * 2021-11-29 2022-03-08 国网江苏省电力有限公司营销服务中心 基于大数据的中高压失准分析系统
CN114240690A (zh) * 2021-12-17 2022-03-25 国网辽宁省电力有限公司沈阳供电公司 基于多聚合用户电表数据的用户用电异常监测方法
CN115203286A (zh) * 2022-05-25 2022-10-18 广州番禺职业技术学院 一种基于知识驱动的台区户变关系识别方法
CN114825647B (zh) * 2022-06-30 2022-09-23 国网浙江省电力有限公司宁波供电公司 一种基于拓扑关系的台区线损监测方法及相关设备
CN115423250B (zh) * 2022-07-28 2023-07-28 国网浙江省电力有限公司营销服务中心 一种台区户变关系分析方法
CN115470858A (zh) * 2022-09-13 2022-12-13 国网天津市电力公司 基于ddtw距离和最小最大距离算法的配电台区拓扑辨识方法
CN115235676B (zh) * 2022-09-26 2022-12-09 山东济矿鲁能煤电股份有限公司阳城煤矿 一种矿井液压支架压力异常自动报警方法
CN115344567A (zh) * 2022-10-18 2022-11-15 国网天津市电力公司营销服务中心 一种适用于边缘计算的低压台区数据清洗治理方法及装置
CN115438035B (zh) * 2022-10-27 2023-04-07 江西师范大学 一种基于kpca和混合相似度的数据异常处理方法
CN115542062B (zh) * 2022-11-07 2024-01-09 北京志翔科技股份有限公司 户变关系异常的识别方法、装置、设备和存储介质
CN116303480B (zh) * 2023-05-22 2023-08-04 新乡市万新电气有限公司 基于云计算的电能表误差校验方法
CN116599055B (zh) * 2023-05-26 2023-12-01 联桥科技有限公司 一种低压配网台区的拓扑网络识别方法及系统
CN116482471B (zh) * 2023-06-21 2023-09-08 四川中电启明星信息技术有限公司 一种基于电压时空聚合曲线的户变关系识别方法
CN116499419B (zh) * 2023-06-27 2023-10-20 西安高商智能科技有限责任公司 一种舵机旋转角度异常检测方法及系统
CN116522181B (zh) * 2023-06-28 2024-03-29 广东电网有限责任公司 一种配电台区日线损特征提取方法及系统
CN116595426B (zh) * 2023-07-17 2023-09-26 济南大陆机电股份有限公司 一种工业物联网数据智能采集管理系统
CN116910596B (zh) * 2023-07-26 2024-06-21 江苏方天电力技术有限公司 一种基于改进dbscan聚类的用户窃电分析方法、设备及存储介质
CN116881745B (zh) * 2023-09-07 2023-12-19 深圳市特安工业科技有限公司 基于大数据的压力变送器异常监测方法
CN116954342B (zh) * 2023-09-18 2023-12-15 惠州市鑫晖源科技有限公司 基于物联网的全方位高效散热主机电源运行监测方法
CN117129790B (zh) * 2023-10-26 2024-01-23 山西思极科技有限公司 一种电力系统故障诊断系统
CN117148023B (zh) * 2023-10-31 2024-01-26 威海海泰电子有限公司 一种智能电源适配器生产检测方法及系统
CN117648657B (zh) * 2023-12-13 2024-05-14 青岛市建筑设计研究院集团股份有限公司 一种城市规划多源数据优化处理方法
CN118033249A (zh) * 2024-04-11 2024-05-14 国网江苏省电力有限公司常州供电分公司 台区相位识别方法和装置
CN118051862B (zh) * 2024-04-16 2024-06-21 洛阳禾安工程技术服务有限公司 一种建筑涂料的检测方法及系统
CN118070080B (zh) * 2024-04-17 2024-07-30 山东中电仪表有限公司 一种多功能电能表的用户电耗数据智能分析方法及系统
CN118536050A (zh) * 2024-07-29 2024-08-23 国网江西省电力有限公司南昌供电分公司 基于gis和电压数据融合的配电网线变关系识别方法
CN118551332A (zh) * 2024-07-30 2024-08-27 杭州致成电子科技有限公司 一种基于r方的台户关系识别方法、系统、设备和介质
CN118610903A (zh) * 2024-08-02 2024-09-06 国网山东省电力公司营销服务中心(计量中心) 一种上网关口电能计量装置及其计量误差判断方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104376057A (zh) * 2014-11-06 2015-02-25 南京邮电大学 一种基于最大最小距离和K-means的自适应聚类方法
CN109685128A (zh) * 2018-12-18 2019-04-26 电子科技大学 一种MB-kmeans++聚类方法及基于其的用户会话聚类方法
CN110084718A (zh) * 2019-04-25 2019-08-02 国网湖南省电力有限公司 一种台区户变关系精准校验方法及装置
CN110727662A (zh) * 2019-09-10 2020-01-24 国网浙江省电力有限公司电力科学研究院 基于相关性分析的低压台区用户相位识别方法和系统
CN111061821A (zh) * 2019-11-12 2020-04-24 国网山东省电力公司聊城供电公司 基于改进k值聚类算法的低压配电网拓扑校验方法及系统
CN111126429A (zh) * 2019-11-10 2020-05-08 国网浙江省电力有限公司 一种基于PCA降维和K-Means聚类的低压台区用户接入点识别方法
CN111444241A (zh) * 2020-03-26 2020-07-24 南京工程学院 一种基于数据挖掘的台区线损异常关联用户精准定位方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8095521B2 (en) * 2007-03-30 2012-01-10 Amazon Technologies, Inc. Recommendation system with cluster-based filtering of recommendations
CN106014858B (zh) * 2016-07-21 2019-11-22 浙江运达风电股份有限公司 一种风电机组对风误差自动校准方法及装置
CN109325545A (zh) * 2018-10-16 2019-02-12 南方电网科学研究院有限责任公司 低压配电网拓扑结构校验方法、装置、设备及存储介质
CN109461096A (zh) * 2018-11-05 2019-03-12 杭州致成电子科技有限公司 一种基于电量关联性的户变关系异常判别方法
CN110120668B (zh) * 2019-05-08 2021-09-07 许昌许继软件技术有限公司 一种台区拓扑自动识别方法及系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104376057A (zh) * 2014-11-06 2015-02-25 南京邮电大学 一种基于最大最小距离和K-means的自适应聚类方法
CN109685128A (zh) * 2018-12-18 2019-04-26 电子科技大学 一种MB-kmeans++聚类方法及基于其的用户会话聚类方法
CN110084718A (zh) * 2019-04-25 2019-08-02 国网湖南省电力有限公司 一种台区户变关系精准校验方法及装置
CN110727662A (zh) * 2019-09-10 2020-01-24 国网浙江省电力有限公司电力科学研究院 基于相关性分析的低压台区用户相位识别方法和系统
CN111126429A (zh) * 2019-11-10 2020-05-08 国网浙江省电力有限公司 一种基于PCA降维和K-Means聚类的低压台区用户接入点识别方法
CN111061821A (zh) * 2019-11-12 2020-04-24 国网山东省电力公司聊城供电公司 基于改进k值聚类算法的低压配电网拓扑校验方法及系统
CN111444241A (zh) * 2020-03-26 2020-07-24 南京工程学院 一种基于数据挖掘的台区线损异常关联用户精准定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
应用皮尔逊相关系数算法查找异常电能表用户;王涓等;《电力需求侧管理》;第16卷(第02期);52-54 *

Also Published As

Publication number Publication date
WO2022110557A1 (zh) 2022-06-02
CN112699913A (zh) 2021-04-23

Similar Documents

Publication Publication Date Title
CN112699913B (zh) 一种台区户变关系异常诊断方法及装置
CN111505433B (zh) 一种低压台区户变关系纠错及相位识别方法
CN109818349B (zh) 一种基于多维状态矩阵滑动匹配的电网鲁棒状态预测方法
CN106055918A (zh) 一种电力系统负荷数据辨识及修复方法
CN109670676A (zh) 基于支持向量数据描述的配网台区风险预警方法及系统
CN117214637B (zh) 充电枪耐压测试方法、装置、设备及存储介质
CN108334894B (zh) 基于无监督机器学习的变压器油温异常识别方法
CN116011158A (zh) 一种低压台区的拓扑识别方法、系统及装置
CN112418687B (zh) 基于用电特征的用户用电异常识别方法、装置和存储介质
CN111709668A (zh) 基于数据挖掘技术的电网设备参数风险识别方法及装置
CN111626559A (zh) 一种基于主因子分析法的中压配电网线损关键特征指标提取方法及系统
CN110991741B (zh) 一种基于深度学习的断面约束概率预警方法及系统
CN110544047A (zh) 一种不良数据辨识方法
CN112926686B (zh) 基于brb和lstm模型的电力大数据用电异常检测方法及装置
CN113204894B (zh) 一种电能计量异常诊断模型的构建方法及应用
CN117236022A (zh) 互感器剩余寿命预测模型训练方法、应用方法及电子设备
CN117744856A (zh) 锂电池剩余寿命预测方法、系统、计算机设备和存储介质
CN109784777B (zh) 基于时序信息片段云相似度度量的电网设备状态评估方法
CN112036718A (zh) 一种考虑新能源不确定性的电力系统安全风险评估方法
CN116662840A (zh) 一种基于机器学习的低压台区用户相位识别方法
CN115545422A (zh) 一种基于改进决策机制的台区户变关系识别方法
CN115598459A (zh) 一种配电网10kV馈线故障停电预测方法
CN113487080A (zh) 一种基于风速分类的风速动态场景生成方法、系统及终端
CN118171042B (zh) 时间加权的电压曲线相似度与不相似度综合拓扑识别方法
CN117371623B (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