基于空间特征参数的变电站多源局部放电检测方法及装置

摘要

本发明公开了基于空间特征参数的变电站多源局部放电检测方法及装置。该方法包括：在预设时间内采集UHF天线阵列接收的超高频信号；从所述UHF天线阵列获取所述超高频信号的脉冲波形的片段组，组成波形片段组序列；根据所述波形片段组序列计算各根天线的时延，计算空间特征向量并形成序列，对空间特征向量序列进行模糊聚类分析；将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型。采用本发明，可以通过基于不同局部放电源所产生超高频放电信号在时域与频率的特征不同，运用模糊聚类的方法对其进行聚类，以分离不同局部放电源所产生的局部放电UHF信号。

权利要求

1.一种基于空间特征参数的变电站多源局部放电信号检测分类方法，其特征在于，包括：在预设时间内采集UHF天线阵列接收的超高频信号；

从所述UHF天线阵列获取所述超高频信号的脉冲波形的片段组，组成波形片段组序列；

根据所述波形片段组序列计算各根天线的时延，确定空间特征向量并形成序列，以每个空间特征向量作为假定中心，根据其它空间特征向量与该假定中心的距离，获取该假定中心的可能性权值；

取可能性权值最大的假定中心作为第一聚类中心，并将该可能性权值作为最大权值Pmax；

设定可能性权值最小值Pmin＝ε×Pmax，其中，ε为预设的分类系数，由所有采集的UHF信号提取的空间特征向量组成一集合{x1,x2,…,xn},并假设集合中每一向量元素都为潜在的聚类中心，以如下公式(a1)评估各点为聚类中心的可能性Pi，取可能性最大的点为所述第一聚类中心：其中，Pi是集合中某点与其他所有点的距离的函数，ra用于定义“邻居”的有效半径距离；

* *

假设x1，P1 为计算出来的第一个聚类中心及其为聚类中心的可能性，根据以下公式(a2)重新计算各点为聚类中心的可能性：*

其中，rb用于定义需要减少为聚类中心可能性Pi的x1 的“邻居”的有效半径，这样*x1 附近的点为潜在聚类中心的可能性大大减少；

*

再依据以上公式(a2)的计算结果选取新的Pi最大值P2 对应的点为第二个聚类中心* * *x2，以此类推，求出第k个聚类中心xk、Pk 之后，根据以下公式(a3)求取各点新的Pi：取最大值对应的点定位为第k+1个聚类中心，根据所述可能性权值Pi>Pmin,确定分类数；

根据所述空间特征向量获取模糊聚类的目标函数的极小值；

在所述目标函数取极小值时，将各个聚类原型矢量组成最佳模糊分类矩阵；

根据所述最佳模糊分类矩阵，分类多源局部放电信号，并获取各类分类信号的聚类中心,对空间特征向量序列进行模糊聚类分析；

将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型。

2.根据权利要求1所述的基于空间特征参数的变电站多源局部放电信号检测分类方法，其特征在于：所述分类系数ε的范围包括，0.15≤ε≤0.2。

3.根据权利要求2所述的基于空间特征参数的变电站多源局部放电信号检测分类方法，其特征在于，将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型的步骤，包括：当所述聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征不相匹配时，确定该分类为干扰信号的波形。

4.一种基于空间特征参数的变电站多源局部放电信号检测分类装置，其特征在于，包括：信号采集单元，用于在预设时间内采集UHF天线阵列接收的超高频信号；

与所述信号采集单元相连的波形片段获取单元，用于从所述UHF天线阵列获取所述超高频信号的脉冲波形的片段组，组成波形片段组序列；

与所述波形片段获取单元相连的聚类分析单元，用于根据所述波形片段组序列计算各根天线的时延，确定空间特征向量并形成序列，对空间特征向量序列进行模糊聚类分析；

所述聚类分析单元包括：

参数设定单元，用于计算设定聚类分析的分类数，根据所述空间特征向量获取模糊聚类的目标函数的极小值；

与所述参数设定单元相连的最佳分类单元，用于在所述目标函数取极小值时，将各个聚类原型矢量组成最佳模糊分类矩阵；

与所述最佳分类单元相连的聚类中心获取单元，用于根据所述最佳模糊分类矩阵，分类多源局部放电信号获取各类分类信号的聚类中心；

所述参数设定单元还包括：

可能性分析单元，用于以每个空间特征向量作为假定中心，根据其它空间特征向量与该假定中心的距离，获取该假定中心的可能性权值；

与所述可能性分析单元相连的第一聚类中心单元，用于取可能性权值最大的假定中心作为第一聚类中心，并将该可能性权值作为最大权值Pmax；

与所述第一聚类中心单元相连的聚类中心确定单元，用于设定可能性权值最小值Pmin＝ε×Pmax，其中，ε为预设的分类系数，由所有采集的UHF信号提取的空间特征向量组成一集合{x1,x2,…,xn},并假设集合中每一向量元素都为潜在的聚类中心，以如下公式(a4)评估各点为聚类中心的可能性Pi，取可能性最大的点为所述第一聚类中心：其中，Pi是集合中某点与其他所有点的距离的函数，ra用于定义“邻居”的有效半径距离；

* *

假设x1，P1 为计算出来的第一个聚类中心及其为聚类中心的可能性，根据以下公式(a5)重新计算各点为聚类中心的可能性：*

其中，rb用于定义需要减少为聚类中心可能性Pi的x1 的“邻居”的有效半径，这样*x1 附近的点为潜在聚类中心的可能性大大减少；

*

再依据以上公式(a5)的计算结果选取新的Pi最大值P2 对应的点为第二个聚类中心* * *x2，以此类推，求出第k个聚类中心xk,Pk 之后，根据以下公式(a6)求取各点新的Pi：取最大值对应的点定位为第k+1个聚类中心；

与所述第一聚类中心单元相连的聚类中心确定单元，还用于根据所述可能性权值Pi>Pmin,确定分类数；

与所述聚类分析单元相连的对比识别单元，用于将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型。

5.根据权利要求4所述的基于空间特征参数的变电站多源局部放电信号检测分类装置，其特征在于，所述聚类中心确定单元包括：分类系数设定单元，用于预设所述分类系数ε。

6.根据权利要求5所述的基于空间特征参数的变电站多源局部放电信号检测分类装置，其特征在于，还包括：与所述对比识别单元相连的干扰判断单元，用于在所述聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征不相匹配时，确定该分类为干扰信号的波形。

说明书

基于空间特征参数的变电站多源局部放电检测方法及装置

技术领域

[0001] 本发明涉及电力系统高电压与绝缘技术领域，尤其涉及基于空间特征参数的变电站多源局部放电检测方法及装置。

背景技术

[0002] 国内现有对变电站的局部放电进行监测及定位的方法基本都是先针对变电站内单个电力设备进行局部放电监测与定位，然后综合各个电力设备的监测情况评估整个变电站电力设备的局部放电状况。针对单个变电站设备(GIS地理信息系统、变压器、容性设备等)的局部放电检测，一般通过监测局部放电所产生的电信号、射频信号、声信号或者其他化学量或物理量来实现。在长期的检测与维修过程中，因为设备一直处于维护过程中，单个设备多处同时发生局部放电的可能性比较小，目前应用到单个设备的局部放电监测系统基本都不考虑多源局部放电的情况。

[0003] 而针对整个变电站进行局部放电监测，变电站中的任何高压电力设备均可能会发生局部放电故障，要想对全站的一次电气设备实施监测，采用上述监测方式，就需要在每一个设备上都安装局部放电监测装置，这样可以自动区分不同电力设备的局部放电信号，但是这样需要安装多套局部放电监测装置，需要耗费大量的时间、财力进行设备采购及安装；测试时需携带多套不同类型仪器，操作不便；众多监测装置的维护和管理也需耗费大量的时间和人力。目前的监测形式难以适应今后智能变电站远程管理的发展需要。

发明内容

[0004] 为了解决上述技术问题，本发明提供了基于空间特征参数的变电站多源局部放电检测方法及装置，能够检测由不同局部放电源所产生的局部放电UHF(Ultra High Frequency，超高频)信号。

[0005] 本发明提供了基于空间特征参数的变电站多源局部放电信号检测及分类方法，包括：

[0006] 在预设时间内采集UHF天线阵列接收的超高频信号；

[0007] 从所述UHF天线阵列获取所述超高频信号的脉冲波形的片段组，组成波形片段组序列；

[0008] 根据所述波形片段组序列计算各根天线的时延，确定空间特征向量并形成序列，对空间特征向量序列进行模糊聚类分析；

[0009] 将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型。

[0010] 相应地，本发明还提供了基于空间特征参数的变电站多源局部放电信号检测及分类装置，包括：

[0011] 信号采集单元，用于在预设时间内采集UHF天线阵列接收的超高频信号；

[0012] 与所述信号采集单元相连的波形片段获取单元，用于从所述UHF天线阵列获取所述超高频信号的脉冲波形的片段组，组成波形片段组序列；

[0013] 与所述波形片段获取单元相连的聚类分析单元，用于根据所述波形片段组序列计算各根天线的时延，确定空间特征向量并形成序列，对空间特征向量序列进行模糊聚类分析；

[0014] 与所述聚类分析单元相连的对比识别单元，用于将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型。

[0015] 实施本发明，具有如下有益效果：

[0016] 本发明在变电站场地内安装若干个宽带全向UHF天线传感器阵列接收局部放电发出的电磁波信号，可用一套装置对全站设备进行局部放电的监测，获取和定位整个变电站站域的放电情况。通过所采集波形片段能量积累函数的拐点计算天线阵列各天线所接收的局部放电脉冲信号的起始点，从而计算局部放电信号到达天线阵列各天线的时延。基于步时延及天线坐标计算所接收局部放电信号源的空间特征参数。运用模糊聚类的方法对其进行聚类，以分离不同局部放电源所产生的局部放电UHF信号。进一步地，可将其分类结果与系统定位结果的分类结果进行交叉对比，以确认多局部放电源的存在与定位的准确性。

附图说明

[0017] 图1是典型放电UHF信号波形及其频谱；

[0018] 图2是本发明基于空间特征参数的变电站多源局部放电信号检测及分类方法的流程图；

[0019] 图3是本发明基于空间特征参数的变电站多源局部放电信号检测及分类方法的示意图；

[0020] 图4是本发明基于空间特征参数的变电站多源局部放电信号检测及分类方法的第一实施例流程图；

[0021] 图5是本发明基于空间特征参数的变电站多源局部放电信号检测分类方法的第二实施例流程图；

[0022] 图6是本发明基于空间特征参数的变电站多源局部放电信号检测及分类装置的示意图；

[0023] 图7是本发明基于空间特征参数的变电站多源局部放电信号检测及分类装置的第一实施例的聚类分析单元示意图；

[0024] 图8是本发明基于空间特征参数的变电站多源局部放电信号检测及分类装置的第二实施例的参数设定单元示意图。

具体实施方式

[0025] 为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步地详细描述。

[0026] 绝缘故障是电力设备在运行中最主要的可能故障之一，电力设备发生绝缘故障前，一般都会有一个逐渐发展的局部放电过程，最终导致绝缘击穿。如果在这个过程能够对运行设备进行局部放电监测和诊断，及时发现局部放电信号，提前对缺陷进行处理，就能有效避免绝缘击穿故障的发生。对局部放电位置的定位，也有助于制定更有针对性的检修处理方案，减少停电时间，提高检修效率。

[0027] 超高频(UHF)电磁波法是局部放电检测的一种新方法，该方法通过UHF天线传感器接收当电力设备内发生局部放电过程辐射的电磁波中的超高频段(300～3000MHz)信号来检测局部放电。超高频电磁波检测的优点为：检测频段较高，可以有效地避开常规局部放电测量中的电晕、开关操作等多种电气干扰；检测频带宽，所以其检测灵敏度很高，且电磁波在空气中的传播速度近似光速，通过构建一个UHF天线矩阵，利用传播到天线矩阵内各个UHF天线的UHF信号的时差可以用来计算局部放电源的位置。相对于单个电力设备的局部放电监测而言，由于整个变电站的局部放电监测的电力设备对象多，所以多个电力设备同时发生局部放电的可能性大，在监测过程中区分不同局部放电源所激发的UHF信号，以便更加准确的定位。

[0028] 图1是典型放电UHF信号波形及其频谱。

[0029] 如图1所示，对变电站电力设备的局部放电信号进行分类整理，可以将其归纳为针尖放电，沿面放电、悬浮放电、金属微粒、绝缘内部放电等。另外电力设备中局部放电发生空间的绝缘介质可能为空气，SF6和绝缘油等。理论上而言，不同的放电类型，局部放电形成的机理不同，其放电的过程不尽相同。其所激发的局部放电UHF信号在时域内与频域内的特性则不尽相同。另外局部放电所激发UHF信号的特性受到局部放电空间的绝缘介质有关，同时也和局部放电周围空间电力设备的几何结构与所用材料，接受UHF信号的传感器与局部放电源之间的空间距离与路径有关，因为上述因素会影响UHF信号的传播。

[0030] 试验研究也证实了上述理论分析。以下数据给出了相同的放电形式在不同的绝缘介质中所激发的UHF信号频谱分布不同，以及不同放电形式在相同绝缘介质中所激发的UHF信号频谱分布不同的示例。基于GIS实物模型的四种缺陷模拟试验结果也反映出不同的放电类型其放电波形及频谱不同，试验结果如图1所示。而下表则总结了不同放电类型所产生的UHF信号的频带分布特征。表1：

[0031]

[0032] 以上试验结果充分证明，对于局部放电：相同的绝缘介质不同的放电类型，相同的放电类型不同的绝缘介质，其激发的UHF信号波形及其频谱特征存在很大的差异。对变电站局部放电进行整体监测，可以通过分析所监测到的局部放电UHF信号的波形及频谱特性，用以区分多源局部放电信号。

[0033] 图2是本发明基于空间特征参数的变电站多源局部放电信号检测及分类方法的流程图，包括：

[0034] S101：在预设时间内采集UHF天线阵列接收的超高频信号；

[0035] S102：从所述UHF天线阵列获取所述超高频信号的脉冲波形的片段，组成波形片段组序列；

[0036] S103：根据所述波形片段组序列计算各根天线的时延，确定空间特征向量并形成序列，对空间特征向量序列进行模糊聚类分析；

[0037] S104：将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型。

[0038] 变电站局部放电监测与定位系统所监测的对象为所监测变电站整体。变电站中存在众多的电力设备，如变压器，断路器，互感器，避雷器等，不排除整个变电站存在多个局部放电源，同时产生局部放电的可能性，如果存在这样的情况，判断是否存在多源局部放电信号，并对其进行定位与识别成为变电站整体监测与定位系统的关键性能。本发明在变电站场地内安装宽带全向UHF天线传感器阵列接收局部放电发出的电磁波信号，可用一套装置对全站设备进行局部放电的监测，获取和定位整个变电站站域的放电情况。

[0039] 图3是本发明基于空间特征参数的变电站多源局部放电信号检测及分类方法的示意图。

[0040] 首先，分段同步采集UHF天线阵列接收的UHF脉冲信号组，并持续采集一段时间；将获取的多个局部放电UHF脉冲波形的片段组，组成波形片段组序列；如图3所示，所述UHF天线阵列由四根天线组成时，所述波形片段组包括由这些天线接收的四个UHF脉冲信号。最后，从所述波形片段组序列中计算空间特征参数确定空间特征向量，以该空间特征向量进行模糊聚类分析，根据聚类分析的结果与典型的UHF脉冲波形的特征进行对比，确定所述局部放电UHF脉冲波形的类型。通过将聚类分析的结果与预设的超高频脉冲波形的特征进行对比，可以分离不同局部放电源所产生的局部放电UHF信号，从而对其进行定位与识别。通过该装置首先粗选有缺陷的设备或设备部件，确定设备或设备部件有了缺陷和故障风险后再进行深入分析和定位，成本低，效率高。

[0041] 需要补充说明的是，从所述波形片段组序列中提取空间特征参数确定空间特征向量的过程包括：首先，基于波形片段组，采用能量积累法，通过所采集波形片段能量积累函数的拐点计算天线阵列各天线所接收的局部放电脉冲信号的起始点(t1,t2,t3,t4…)。然后，从而计算局部放电信号到达天线阵列各天线的时延。最后，基于计算的各天线接收的信号的时延，及天线坐标P1(x1,y1,z1)、P2(x2,y2,z2)、P3(x3,y3,z3)、P4(x4,y4,z4)，采用迭代法计算所接收局部放电信号源的的空间特征参数。

[0042] 具体地，空间特征向量的计算原理：

[0043] 设局部放电源为P(x,y,z),以球面波的形式向四周辐射超高频信号。以图3的四根接收天线为例，4个射频传感器接收天线坐标为S1(x1,y1,z1)、S2(x2,y2,z2)、S3(x3,y3,z3)、S4(x4,y4,z4),可固定安装在变电站某一位置，也可安装在可移动的设备间，各传感器间相隔一定距离，避免时延太小，采样设备无法区分。

[0044] 假定超高频局部放电信号以光速v＝3x108m/s传播，设tj为信号到达Sj的时间，易得局部放电源到传感器距离Li(i、j取1、2、3、4)和时间差为：

[0045]

[0046] ΔTj1＝tj-t1 (3-9)

[0047] Li＝v(t1+ΔTj1) (3-10)

[0048] 由(3-8)(3-10)可建立以下球面方程组：

[0049]

[0050] 使残差向量平方和满足一定精度的最优解。即令各残差为f1、f2、f3、f4，取一定精度tol其表达式为：

[0051] f(x,y,z,t)＝(x-xi)2+(y-yi)2+(z-zi)2-v2(t1+ΔTj1)2 (3-12)[0052]

[0053] 满足式(3-13)的(x,y,z,t)即为方程组(3-11)基于最小二乘法牛顿迭代计算的最优解。

[0054] 图4是本发明基于空间特征参数的变电站多源局部放电信号检测及分类方法的第一实施例流程图，与图2相比，本发明实施例的采用传统的模糊聚类方法对空间特征向量进行分类。

[0055] S201：在预设时间内采集UHF天线阵列接收的超高频信号；

[0056] S202：从所述UHF天线阵列获取所述超高频信号的脉冲波形的片段组，组成波形片段组序列；

[0057] S203：根据所述波形片段组序列计算各根天线的时延，确定空间特征向量并形成序列；

[0058] S204：计算设定聚类分析的分类数，根据所述空间特征向量序列获取模糊聚类的目标函数的极小值；

[0059] S205：在所述目标函数取极小值时，将各个聚类原型矢量组成最佳模糊分类矩阵；

[0060] S206：根据所述最佳模糊分类矩阵，分类多源局部放电信号获取各类分类信号的聚类中心；

[0061] S207：将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型。

[0062] 在变电站安装数个宽带全向UHF天线传感器，构成UHF天线阵列，接收全站范围内激发的UHF信号。局部放电监测系统采用脉冲触发方式检测UHF天线接收的UHF信号，并对检测的信号分段进行采集与存储。整个检测过程是连续的，天线阵内n个天线接收的UHF信号是同步采集并存储的，该存储的信号为一含局部放电UHF脉冲的波形片段组序列a。每个局部放电UHF脉冲的波形片段组所含的波形数与UHF天线阵所含UHF天线数相等。系统采样率在1.5GS/s以上。以记录的局部放电UHF脉冲的波形片段组序列中提取出的同一时刻同步采集的UHF脉冲波形片段组(如图3所示)，计算确定时延，并以此时延计算满足式(3-13)的空间特征向量组(x,y,z,t)。

[0063] 采用传统的模糊聚类方法，以提取的空间特征向量为参数，对记录的变电站局部放电UHF脉冲的波形片段进行分类。确定分类数，以及每类信号的典型波形及其空间参数。优选地，本发明使用模糊聚类分析，其原理如下：

[0064] 聚类分析是多元统计分析的一种，它把一个没有类别标记的样本集按某种准则划分成若干个子集，相似的样本尽可能归为一类，不相似的样本划分到不同的类中。模糊数学是用来处理人的经验与认知先验信息的不确定问题，处理含有大量未知信息的局部放电现场检测数据即是这样的问题。利用模糊聚类对局部放电的数据进行处理是有意义的。

[0065] 经过特征提取的信号，可借助模糊聚类对局部放电脉冲信号的特征参数进行处理。假设c为需要划分的类数，则有2≤c≤K，Pi(i＝1，2，…，c)表示第i类的聚类原型矢量，则模糊聚类的目标函数Jm为

[0066]

[0067] 其中μik为第i个脉冲的第k个特征参数，m称为加权指数，dik为某一脉冲样本xk与第i类聚类原型Pi之间的欧几里得距离，表示该样本与聚类原型之间的失真度，即2

[0068] (dik) ＝||xk-pi|| (6)

[0069] 聚类的准则即为取Jm(U,P)的极小值min{Jm(U,P)}。

[0070] 在给定聚类类别数c，初始化聚类原型P以及加权指数m后，可求得使得Jm(U,P)为最小的μik值为

[0071]

[0072] 以及Jm(U,P)为最小值时pi的值为

[0073]

[0074] 通过上述迭代计算即可得到最佳模糊分类矩阵和聚类中心，随后根据最佳模糊分类矩阵得到c类信号，即完成了对局部放电脉冲的模糊聚类处理。将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型。

[0075] 当所述聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征不相匹配时，确定该分类为干扰信号的波形。

[0076] 图5是本发明基于空间特征参数的变电站多源局部放电信号检测及分类方法的第二实施例流程图，与图3的第一实施例相比，图4的流程图根据空间特征向量的迭代分类，自动获取分类数以及各个分类的聚类中心。它的优点在于提高了模糊聚类的智能性和自适应性。

[0077] S301：在预设时间内采集UHF天线阵列接收的超高频信号；

[0078] S302：从所述UHF天线阵列获取所述超高频信号的脉冲波形的片段组，组成波形片组序列；

[0079] S303：根据所述波形片段组序列计算各根天线的时延，确定空间特征向量并形成序列；

[0080] S304：以每个空间特征向量作为假定中心，根据其它空间特征向量与该假定中心的距离，获取该假定中心的可能性权值；

[0081] S305：取可能性权值最大的假定中心作为第一聚类中心，并将该可能性权值作为最大权值Pmax；

[0082] S306：设定可能性权值Pmin＝ε×Pmax，根据所述可能性权值Pi>Pmin,确定分类数；其中，ε为预设的分类系数；

[0083] S307：将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型；

[0084] S308：当所述聚类分析的结果中中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征不相匹配时，确定该分类为干扰信号的波形。

[0085] 传统的聚类的分类数根据数据事先选定，在计算过程中，为了进一步提高运算的自动性，本实施例采用下面算法对聚类数及聚类中心进行自动计算。它是考虑局部放电UHF信号特征参数向量在多维空间分布密度的一种自适应算法。

[0086] 下面结合本发明的空间特征参数，具体介绍局部放电UHF脉冲的分类数及聚类中心的求取：

[0087] 将采集的局部放电UHF信号提取的空间特征参数x，y，z，t构建成一个空间特征向量Xi(i为所采集UHF信号的序号)。由所有采集的UHF信号提取的空间特征向量组成一集合{x1,x2,…,xn}，并假设集合中每一向量元素都为潜在的聚类中心。为按照式(9)，评估各点为聚类中心的可能性Pi，取可能性最大的点为第一聚类中心。

[0088]

[0089] Pi是集合中某点与其他所有点的距离的函数。一个含有更多“邻居”的点，其为聚类中心的可能性更大。ra是一个常数，用来定义“邻居”的有效半径距离。

[0090] 假设x1*，P1*为计算出来的第一个聚类中心及其为聚类中心的可能性，根据式(10)重新计算各点为聚类中心的可能性。

[0091]

[0092] 式中rb也为一常数，用来定义需要减少为聚类中心可能性Pi的x1*的“邻居”的*有效半径。这样x1 附近的点为潜在聚类中心的可能性大大减少。再依据式(10)的计算* *
结果选取新的Pi最大值P2 对应的点为第二个聚类中心x2。以此类推，求出第k个聚类* *
中心xk,Pk 之后可按式(11)求取各点新的Pi，

[0093]

[0094] 取最大值对应的点定位为第k+1个聚类中心，直至所求的新的聚类中心对应的Pi*小于εP1，ε为一设定极小数，实际计算中，可根据情况在0.15～0.2之间选取。

[0095] 本发明基于不同局部放电源所产生超高频放电信号在时域与频率的特征不同，运用模糊聚类的方法对其进行聚类，以分离不同局部放电源所产生的局部放电UHF信号。并将其分类结果与系统定位结果的分类结果进行交叉对比，以确认多局部放电源的存在与定位的准确性。

[0096] 当根据聚类分析的结果中分好类的各类脉冲波形的特征与典型的UHF脉冲波形的特征对比，无法确定所述局部放电UHF脉冲波形的类型时，确定该类型为干扰信号的波形。

[0097] 另外，各种类似放电的干扰信号，如电力电子器件、发动机等，也会更有可能被射频天线接收，本发明具备将干扰分离的能力。当一类放电UHF脉冲波形的片段没有找到相应的定位结果时，可以确定该类型为干扰信号的波形。

[0098] 图6是本发明基于空间特征参数的变电站多源局部放电信号检测及分类装置的示意图，包括：

[0099] 信号采集单元，用于在预设时间内采集UHF天线阵列接收的超高频信号；

[0100] 与所述信号采集单元相连的波形片段获取单元，用于从所述UHF天线阵列获取所述超高频信号的脉冲波形的片段组，组成波形片段组序列；

[0101] 与所述波形片段获取单元相连的聚类分析单元，用于根据所述波形片段组计算各根天线的时延，确定空间特征向量并形成序列，对空间特征向量序列进行模糊聚类分析；

[0102] 与所述聚类分析单元相连的对比识别单元，用于将聚类分析的结果中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征进行对比，确定所述脉冲波形的类型。

[0103] 图6与图2相对应，上述各个单元的运行方式与方法中的相同。

[0104] 图7是本发明基于空间特征参数的变电站多源局部放电信号检测及分类装置的第一实施例的聚类分析单元示意图。在本实施例中，采用传统的聚类分析方法对空间特征向量进行聚类分析。

[0105] 在第一实施例中，所述聚类分析单元包括：

[0106] 参数设定单元，用于计算设定聚类分析的分类数，根据所述空间特征向量获取模糊聚类的目标函数的极小值；

[0107] 与所述参数设定单元相连的最佳分类单元，用于在所述目标函数取极小值时，将各个聚类原型矢量组成最佳模糊分类矩阵；

[0108] 与所述最佳分类单元相连的聚类中心获取单元，用于根据所述最佳模糊分类矩阵，分类多源局部放电信号获取各类分类信号的聚类中心。

[0109] 图7与图4相对应，图中各个单元的运行方式与方法中的相同。

[0110] 图8是本发明基于空间特征参数的变电站多源局部放电信号检测及分类装置的第二实施例的参数设定单元示意图。与图7相比，图8实施例采用智能性和自适应性更强的聚类分析方法对空间特征向量进行聚类分析。

[0111] 在第二实施例中，所述参数设定单元包括：

[0112] 可能性分析单元，用于以每个空间特征向量作为假定中心，根据其它空间特征向量与该假定中心的距离，获取该假定中心的可能性权值；

[0113] 与所述可能性分析单元相连的第一聚类中心单元，用于取可能性权值最大的假定中心作为第一聚类中心，并将该可能性权值作为最大权值Pmax；

[0114] 与所述第一聚类中心单元相连的聚类中心确定单元，用于设定可能性权值Pmin＝ε×Pmax，根据所述可能性权值Pi>Pmin,确定分类数，其中，ε为预设的分类系数。

[0115] 图8与图5相对应，图中各个单元的运行方式与方法中的相同。

[0116] 如图8所示，所述聚类中心确定单元包括：分类系数设定单元，用于预设所述分类系数ε。

[0117] 在其中一个实施例中，本发明基于空间特征参数的变电站多源局部放电信号检测分类装置，还包括：

[0118] 与所述对比识别单元相连的干扰判断单元，用于在所述聚类分析的结果中中分好类的各类脉冲波形的特征与预设的超高频脉冲波形的特征不相匹配时，确定该分类为干扰信号的波形。

[0119] 以上所述的本发明实施方式，并不构成对本发明保护范围的限定。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明的权利要求保护范围之内。