CN109901031A - 用于局部放电信号的信噪分离方法、信息数据处理终端 - Google Patents
用于局部放电信号的信噪分离方法、信息数据处理终端 Download PDFInfo
- Publication number
- CN109901031A CN109901031A CN201910146996.XA CN201910146996A CN109901031A CN 109901031 A CN109901031 A CN 109901031A CN 201910146996 A CN201910146996 A CN 201910146996A CN 109901031 A CN109901031 A CN 109901031A
- Authority
- CN
- China
- Prior art keywords
- signal
- pulse
- class
- partial discharge
- phase
- 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.)
- Granted
Links
Landscapes
- Testing Relating To Insulation (AREA)
Abstract
本发明属于高压电气设备局部放电检测技术领域,公开了一种用于局部放电信号的信噪分离方法、信息数据处理终端;采集局部放电波形信号;对局放信号进行提取,获取单个脉冲信号,记录脉冲的峰值及相位,并绘制相位分布谱图;将所提取脉冲通过多个带通滤波器,获得对应滤波器下脉冲的峰值信息,并通过主成分分析降维,获得特征参数;对特征参数进行聚类分析,将脉冲分为多类,并根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号。本发明能够高效准确提取各个脉冲波形;采用多个带通滤波器下的峰值信息作为特征向量,能够充分反映脉冲波形信息,提高了干扰信号与局放信号分离的精度。
Description
技术领域
本发明属于高压电气设备局部放电检测技术领域,尤其涉及一种用于局部放电信号的信噪分离方法、信息数据处理终端。
背景技术
目前,最接近的现有技术:随着社会现代化步伐的加快,我国电网的规模不断扩大,电力电缆也广泛应用到各级电力传输***中。局部放电的变化和特征在某种程度上显示了电缆存在的缺陷,能够反映电缆绝缘的早期劣化。与此同时,高压电气设备中的大型发电机在定子绝缘劣化过程中,局部放电活动是评判其劣化程度的重要标志,局部放电可作为定子绝缘劣化失效的先兆。电力变压器是电力***中最重要及最昂贵的设备之一,其安全运行的意义重大。在现场运行中,局部放电是导致电力变压器绝缘劣化的重要原因之一,特别是随着电力设备容量、电压等级不断增大的情况下,这个问题更为严重。因此,对高压电气设备施行局部放电在线检测有着重大意义。
然而,在对电缆局部放电进行在线检测时,总是不可避免地受到脉冲干扰的影响。这些脉冲干扰分为周期性脉冲干扰以及随机性脉冲干扰。周期性脉冲干扰主要包括可控硅设备在可控硅开、闭时产生的脉冲干扰信号,且周期性地出现在工频的某相位上;随机性脉冲干扰包括由电气设备的局部放电、开关操作及电焊操作等引起的无规律的随机性干扰。目前国内外主要通过将采集的信号与设置的局部放电阈值进行比较,从而判断高压电气设备是否产生局部放电信号。然而,脉冲型干扰与局部放电信号的时域波形十分相似,仅通过幅值阈值比较,或者其放电参数是否位于规定范围,很容易将脉冲干扰识别为局部放电信号,造成漏判或者误判。因此,有必要对采集的局部放电信号进行提取,获得单个脉冲信号,并根据脉冲波形特征进行逻辑判断,实现有用信号与噪声的分离。对局部放电信号进行在线检测时,由于背景噪声的影响,很难确定单个脉冲信号的起始以及结束边沿。现有技术存在着计算量大、脉冲边沿不清楚等问题,给后续的信号分离带来了不少困难。
目前对高压电气设备局部放电信号的研究方法主要有两种,第一种是事先通过大量的典型局部放电试验,建立局放以及噪声指纹库,然后通过人工神经网络比较放电实例与各种放电模型特征参量的相异度,确定识别结果。比如专利《一种变压器局部放电类型识别方法》(公布号:CN201810522817)分别搭建了油中气隙、油楔、针板、悬浮四种内部绝缘故障时的局部放电模拟试验平台,以上升沿t1、持续时间t2、能量频率分布熵Ed、幅值频率分布熵Ad作为特征参数,建立局放以及噪声指纹库,并构建优化概率神经网络实现脉冲信号的识别。这种方法属于监督识别算法,需要一定的先验知识,因此在对采集的信号进行噪声以及局放识别之前,需要进行大量的局放试验,建立放电特征数据库,这会消耗大量的人力物力资源;另一方面,由于现场环境的复杂性,往往同一典型局放在某一时刻获取的脉冲信号与另一时刻获得的信号有所差异,而局放数据库不可能涵盖所有时刻的脉冲信号,因此在识别过程中很容易出现误判。第二种方法是提取单个脉冲信号的时频域波形特征,采用聚类算法对这些特征进行聚类分析,从而实现信号与噪声的分离。采用脉冲波形特征表征脉冲信号的程度不同,最终干扰信号与局放信号分离的效果也有所不同,因此脉冲信号特征参数的选取十分重要。目前应用最为广泛的是等效时频法,即采用脉冲信号的时域等效宽度与频域等效宽度作为脉冲的特征参数,利用噪声与局放信号之间的特征差异,从而实现信号的分离。但是当两种类型的脉冲信号的时域等效宽度与频域等效宽度相同时,该方法很难完成分离。
综上所述,现有技术存在的问题是:
(1)目前对高压电气设备局部放电信号的研究方法对局部放电信号进行在线检测时,由于背景噪声的影响,很难确定单个脉冲信号的起始以及结束边沿。现有技术存在着计算量大、脉冲边沿不清楚等问题,给后续的信号分离带来了不少困难。
(2)监督识别算法需要一定的先验知识,因此在对采集的信号进行噪声以及局放识别之前,需要进行大量的局放试验,建立放电特征数据库,这会消耗大量的人力物力资源;另一方面,由于现场环境的复杂性,往往同一典型局放在某一时刻获取的脉冲信号与另一时刻获得的信号有所差异,而局放数据库不可能涵盖所有时刻的脉冲信号,因此在识别过程中很容易出现误判。
(3)目前对高压电气设备局部放电信号的研究方法脉冲波形特征表征脉冲信号的程度不同,最终干扰信号与局放信号分离的效果也有所不同,因此脉冲信号特征参数的选取十分重要。
解决上述技术问题的难度:在局部放电信号中提取单个脉冲信号时,需要确定信号的起点与终点,由于背景噪声的影响,单个脉冲信号的边沿很不明显。因此准确提取脉冲波形信号,确定脉冲的起始以及结束位置具有一定的难度。脉冲信号可以用脉冲的特征参数表示,特征参数表征脉冲信号的程度越高,信号分离的效果越好。特征参数的错误构造将会直接导致聚类效果的不理想。
解决上述技术问题的意义:本发明首先采用幅值-时间阈值法,通过确定单个脉冲第一个大于幅值阈值、最后一个大于幅值阈值的点初步确定脉冲的位置,然后根据时域能量法进一步确定脉冲的起始以及结束边沿,能够高效准确地提取各个脉冲波形。本发明提出采用多个带通滤波器下的峰值信息作为特征向量,能够充分反映脉冲波形信息,提高了干扰信号与局放信号分离的精度。本发明在对脉冲信号与局放信号进行分离时,无需建立局放以及噪声指纹库,结合聚类分离结果以及单类信号相位分布谱图特征,即可分辨脉冲干扰以及局部放电信号类型。
发明内容
针对现有技术存在的问题,本发明提供了一种用于局部放电信号的信噪分离方法、信息数据处理终端。
本发明是这样实现的,一种用于局部放电信号的信噪分离方法,所述用于局部放电信号的信噪分离方法采集局部放电波形信号;对局放信号进行提取,获取单个脉冲信号,记录脉冲的峰值及相位,并绘制相位分布谱图;将所提取脉冲通过多个带通滤波器,获得对应滤波器下脉冲的峰值信息,并通过主成分分析降维,获得特征参数;对特征参数进行聚类分析,将脉冲分为多类,并根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号。
进一步,所述对局放信号进行提取,获取单个脉冲信号,记录脉冲的峰值及相位,并绘制相位分布谱图具体包括:
步骤一,采集的局部放电波形信号为离散时间序列,将采集信号的幅值与设置的脉冲幅值阈值进行比较;若采集信号幅值大于所述脉冲幅值阈值,则存储该幅值以及所对应时间序列中的位置信息;否则,不进行存储;
步骤二,读取存储的采样点,将相邻采样点时间之差与设置的脉冲时间阈值进行比较,确定脉冲的第一个大于脉冲阈值的点P以及最后一个大于脉冲阈值的点Q,并分别进行存储;
步骤三,根据步骤二中存储的采样点,确定脉冲的起点以及终点,并存储整个脉冲的波形,取该脉冲幅值的最大值作为该脉冲的峰值;并将其所对应的时间序列进行转换确定局放相位,绘制相位分布谱图。
进一步,所述步骤三中确定脉冲起点以及终点方法具体实现如下:
第一步,计算背景噪声的均值以及方差,背景噪声可认为是服从于均值为μ,方差为σ2的正态分布的数字随机序列,均值以及方差计算公式:
式中,Xi为背景噪声序列,N为截取的背景噪声的点数;
第二步,计算背景噪声的时域能量,定义长度为M的信号序列时域能量如式;
第三步,确定脉冲的起始以及结束位置,分别在点P处进行前向搜索、点Q处进行后向搜索,计算长度为M的信号的时域能量;将时域能量与设置的能量阈值进行比较,当时域能量小于等于能量阈值时,完成搜索,确定脉冲的起始位置和结束位置。
进一步,所述将所提取脉冲通过多个带通滤波器,获得对应滤波器下脉冲的峰值信息,并通过主成分分析降维,获得特征参数具体包括以下步骤:
步骤一,对采集信号进行傅里叶变换得到频谱信号,获取采集信号的频率范围;
步骤二,在采集信号频率范围内,以0.1MHz为步进,计算提取的所有脉冲在相应频率下的方差;
步骤三,取方差较大的前三个值所在频率作为第一带通滤波器、第二带通滤波器、第三带通滤波器的中心频率fc1、fc2、fc3;沿方差曲线分别在所述中心频率两侧取3dB带宽作为第一带通滤波器、第二带通滤波器、第三带通滤波器的带宽BW1、BW2、BW3;
步骤四,将提取的单个脉冲信号分别通过三个带通滤波器,获得多个频带下脉冲的多个峰值,分别为第一峰值、第二峰值、第三峰值,并进行归一化处理;
步骤五,通过主成分分析降维算法对步骤四中所述的多个峰值进行降维,提取其前两阶X1、X2作为特征参数。
进一步,所述对特征参数进行聚类分析,将脉冲分为多类,并根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号具体包括:
步骤一,根据特征参数(X1,X2)绘制二维坐标散点图;
步骤二,采用模糊C均值聚类算法对特征参数进行聚类,得到分离的噪声信号与局放信号;
步骤三,根据分类类别号,使用不同图例在同一坐标系下绘制不同类别脉冲的二维坐标散点图;
步骤四,根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号。
进一步,采用聚类算法进行聚类时,聚类类别数有两种确定方式,一是人为确定,二是根据聚类有效性指标进行评分,计算聚类有效性函数,自适应确定最佳聚类类别数。
进一步,所述聚类有效性指标包括类内紧密性指标J以及类间分离性指标K,其中类内紧密性指标J表征类内样本的紧聚程度,类内紧密性指标J越小,类内样本的紧聚程度越高,各子类内样本点越相似,其计算公式如式:
类间分离性指标K表征类与类之间的分散程度,类间分离性指标K越大,类与类之间的分散程度越大,类间样本点的差异性越大,其计算公式如式:
其中,N为样本个数,c为聚类数,μij为第j个点与第i个分类簇之间的隶属关系;
聚类有效性函数f根据类内紧密性指标J以及类间分离性指标K,寻求紧致性最小且分离性最大的最佳聚类分类数c,其计算公式如式:
聚类有效性函数f的最小值对应最佳的聚类数c,即可自适应确定最佳聚类数。
本发明的另一目的在于提供一种执行所述用于局部放电信号的信噪分离方法的用于局部放电信号的信噪分离***,所述用于局部放电信号的信噪分离***包括:
信号采集模块,用于采集局部放电波形信号;
脉冲提取模块,用于对所述局放信号进行提取,获取单个脉冲信号,记录脉冲的峰值及相位并绘制相位分布谱图;
特征参数提取模块,用于将所提取脉冲通过多个带通滤波器,获得对应滤波器下脉冲的峰值信息,并进行主成分分析降维,获得特征参数;
聚类模块,用于对特征参数进行聚类分析,将脉冲分为多类,并根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号。
本发明的另一目的在于提供一种包含所述用于局部放电信号的信噪分离方法的信息数据处理终端。
本发明的另一目的在于提供一种包含所述用于局部放电信号的信噪分离方法的高压电气设备局部放电检测***。
综上所述,本发明的优点及积极效果为:本发明采用幅值-时间阈值法与时域能量法联合提取脉冲信号,首先通过幅值-时间阈值法初步确定脉冲起点P和终点Q,然后通过时域能量法精确定位单个脉冲的起始和结束边沿,能够高效准确提取各个脉冲波形。传统的时域能量法仅设定一个幅值阈值,很容易对一个脉冲的所有大于幅值阈值的点开展时域能量法确定脉冲的起点与终点,当P、Q之间有n个大于幅值阈值的点时,传统的时域能量法需要对(n+2)个点开展时域能量法求解,而本发明仅需在P、Q这2点处通过时域能量法搜索脉冲的起点与终点,不仅计算量明显减少,而且不会造成脉冲的重复提取,效率大大提高。采用多个带通滤波器下的峰值信息作为特征向量,能够充分反映脉冲波形信息,提高了干扰信号与局放信号分离的精度。常用的时频等效法以时域等效宽度、频域等效宽度作为脉冲信号的特征参数,通过计算聚类有效性函数比较两种算法信号分离的精度。根据相关公式,计算聚类个数c=3、特征量为时域等效宽度和频域等效宽度时,时频等效法的聚类有效性函数f=0.3144。本发明的聚类有效性函数f=0.1829。聚类有效性函数f越小,表明类内样本的紧聚程度越高,各子类内样本点越相似;类与类之间的分散程度越大,类间样本点的差异性越大。可以看出,本发明分离的精度更高。
附图说明
图1是本发明实施例提供的用于局部放电信号的信噪分离方法流程图。
图2是本发明实施例提供的用于局部放电信号的信噪分离方法实现流程图。
图3是本发明实施例提供的提取的单个脉冲波形示意图。
图4是本发明实施例提供的局部放电信号相位分布谱图。
图5是本发明实施例提供的聚类分析的二维坐标散点图。
图6是本发明实施例提供的聚类类别号为1的相位分布谱图。
图7是本发明实施例提供的聚类类别号为2的相位分布谱图。
图8是本发明实施例提供的聚类类别号为3的相位分布谱图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
针对目前对高压电气设备局部放电信号的研究方法需要进行大量的局放试验,耗费人力物力资源;脉冲信号提取计算量大、脉冲边沿不清楚;特征参数的错误构造将会直接导致聚类效果不理想的问题,本发明采用幅值-时间阈值法与时域能量法联合提取脉冲信号,能够高效准确提取各个脉冲波形;采用多个带通滤波器下的峰值信息作为特征向量,能够充分反映脉冲波形信息,提高了干扰信号与局放信号分离的精度。
下面结合附图对本发明的应用原理作详细的描述。
如图1所示,本发明实施例提供的用于局部放电信号的信噪分离方法包括以下步骤:
S101:采集局部放电波形信号;
S102:对局放信号进行提取,获取单个脉冲信号,记录脉冲的峰值及相位,并绘制相位分布谱图;
S103:将所提取脉冲通过多个带通滤波器,获得对应滤波器下脉冲的峰值信息,并通过主成分分析降维,获得特征参数;
S104:对特征参数进行聚类分析,将脉冲分为多类,并根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号。
在本发明的优选实施例中,对局放信号进行提取,获取单个脉冲信号,记录脉冲的峰值及相位,并绘制相位分布谱图具体包括:
步骤一,采集的局部放电波形信号为离散时间序列,将采集信号的幅值与设置的脉冲幅值阈值进行比较;若采集信号幅值大于所述脉冲幅值阈值,则存储该幅值以及所对应时间序列中的位置信息;否则,不进行存储;
步骤二,读取存储的采样点,将相邻采样点时间之差与设置的脉冲时间阈值进行比较,确定脉冲的第一个大于脉冲阈值的点P以及最后一个大于脉冲阈值的点Q,并分别进行存储;
步骤三,根据步骤二中存储的采样点,确定脉冲的起点以及终点,并存储整个脉冲的波形,取该脉冲幅值的最大值作为该脉冲的峰值;并将其所对应的时间序列进行转换确定局放相位,绘制相位分布谱图。
在本发明的优选实施例中,步骤三中确定脉冲起点以及终点方法具体实现如下:
第一步,计算背景噪声的均值以及方差,背景噪声可认为是服从于均值为μ,方差为σ2的正态分布的数字随机序列,均值以及方差计算公式:
式中,Xi为背景噪声序列,N为截取的背景噪声的点数;
第二步,计算背景噪声的时域能量,定义长度为M的信号序列时域能量如式;
第三步,确定脉冲的起始以及结束位置,分别在点P处进行前向搜索、点Q处进行后向搜索,计算长度为M的信号的时域能量;将时域能量与设置的能量阈值进行比较,当时域能量小于等于能量阈值时,完成搜索,确定脉冲的起始位置和结束位置。
在本发明的优选实施例中,将所提取脉冲通过多个带通滤波器,获得对应滤波器下脉冲的峰值信息,并通过主成分分析降维,获得特征参数具体包括以下步骤:
步骤一,对采集信号进行傅里叶变换得到频谱信号,获取采集信号的频率范围;
步骤二,在采集信号频率范围内,以0.1MHz为步进,计算提取的所有脉冲在相应频率下的方差;
步骤三,取方差较大的前三个值所在频率作为第一带通滤波器、第二带通滤波器、第三带通滤波器的中心频率fc1、fc2、fc3;沿方差曲线分别在所述中心频率两侧取3dB带宽作为第一带通滤波器、第二带通滤波器、第三带通滤波器的带宽BW1、BW2、BW3;
步骤四,将提取的单个脉冲信号分别通过三个带通滤波器,获得多个频带下脉冲的多个峰值,分别为第一峰值、第二峰值、第三峰值,并进行归一化处理;
步骤五,通过主成分分析降维算法对步骤四中所述的多个峰值进行降维,提取其前两阶X1、X2作为特征参数。
在本发明的优选实施例中,对特征参数进行聚类分析,将脉冲分为多类,并根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号具体包括:
步骤一,根据特征参数(X1,X2)绘制二维坐标散点图;
步骤二,采用模糊C均值聚类算法对特征参数进行聚类,得到分离的噪声信号与局放信号;
步骤三,根据分类类别号,使用不同图例在同一坐标系下绘制不同类别脉冲的二维坐标散点图;
步骤四,根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号。
在本发明的优选实施例中,采用聚类算法进行聚类时,聚类类别数有两种确定方式,一是人为确定,二是根据聚类有效性指标进行评分,计算聚类有效性函数,自适应确定最佳聚类类别数。
在本发明的优选实施例中,聚类有效性指标包括类内紧密性指标J以及类间分离性指标K,其中类内紧密性指标J表征类内样本的紧聚程度,类内紧密性指标J越小,类内样本的紧聚程度越高,各子类内样本点越相似,其计算公式如式:
类间分离性指标K表征类与类之间的分散程度,类间分离性指标K越大,类与类之间的分散程度越大,类间样本点的差异性越大,其计算公式如式:
其中,N为样本个数,c为聚类数,μij为第j个点与第i个分类簇之间的隶属关系;
聚类有效性函数f根据类内紧密性指标J以及类间分离性指标K,寻求紧致性最小且分离性最大的最佳聚类分类数c,其计算公式如式:
聚类有效性函数f的最小值对应最佳的聚类数c,即可自适应确定最佳聚类数。
图2为本发明实施例提供的用于局部放电信号的信噪分离方法的具体实现流程图。
下面结合具体实施例对本发明的应用原理作进一步的描述。
根据采集的局部放电信号,将幅值-时间阈值法与时域能量法相结合,提取脉冲信号,图3给出了提取的单个脉冲信号的一个波形实例。然后根据提取的所有脉冲的峰值及相位分布,绘制局部放电信号相位分布谱图,如图4所示。图4中脉冲分布在0~360°范围内,可以看出,仅通过阈值比较,无法将脉冲干扰和局部放电信号分离。
对提取的脉冲信号进行傅里叶变换,以0.1MHz为步进,计算提取的所有脉冲在相应频率下的方差,取方差较大的前三个值所在频率作为第一带通滤波器、第二带通滤波器和第三带通滤波器的中心频率。最终确定三个带通滤波器中心频率分别为1.1MHz、2.8MHz、19.5MHz,带宽分别为0.4MHz、0.3MHz、1.8MHz。然后将脉冲信号分别通过这三个带通滤波器,获取对应滤波器下脉冲的峰值信息,归一化后作为特征参数。对获取的三维特征参数进行主成分分析降维,映射到二维坐标系(X1,X2)中,并进行聚类分析。已知聚类个数c=3。
聚类分析后的二维坐标散点图如图5所示,脉冲信号分为三类,如图例所示,分别用圆点标记表示分类类别号为1的脉冲信号,星号标记表示分类类别号为2的脉冲信号,三角标记表示分类类别号为3的脉冲信号。根据相关公式,计算聚类有效性函数f=0.1829。
根据分类类别号、脉冲峰值、脉冲相位信息分别绘制各类信号的相位分布谱图,如图6、图7以及图8所示。由图6可以看出,分类类别号为1的脉冲分布在0~360°相位范围内,与工频信号的相位没有相关性;由图7可以看出,分类类别号为2的脉冲分布在工频信号正负半周的峰值附近,分布较为对称,幅值主要分布在70~300mV;由图8可以看出,分类类别号为3的脉冲分布在工频信号正半周的峰值附近以及负半周内,正负半周分布不对称,且正半周脉冲信号的峰值高于负半周脉冲信号的峰值。根据相位分布特征,可以初步确定分类类别号为1的脉冲信号为干扰脉冲信号,分类类别号为2的脉冲信号为局部放电中的气隙放电,分类类别号为3的脉冲信号为局部放电中的电晕放电。由此可以看出,本发明提出的一种用于局部放电信号的信噪分离方法能够有效地分离局部放电脉冲信号与干扰脉冲信号。
为了验证本发明算法的有效性,下面采用时频等效法进行聚类,以时域等效宽度、频域等效宽度作为脉冲信号的特征参数,通过聚类有效性函数进行对比分析。根据相关公式,计算聚类个数c=3、特征量为时域等效宽度和频域等效宽度时,聚类有效性函数f=0.3144。本发明的聚类有效性函数f=0.1829。聚类有效性函数f越小,表明类内样本的紧聚程度越高,各子类内样本点越相似;类与类之间的分散程度越大,类间样本点的差异性越大。因此,本发明提出的采用多个带通滤波器下的峰值信息作为特征向量,能够充分反映脉冲波形信息,提高了干扰信号与局放信号分离的精度。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种用于局部放电信号的信噪分离方法,其特征在于,所述用于局部放电信号的信噪分离方法采用幅值-时间阈值法与时域能量法联合提取脉冲信号,能够高效准确提取各个脉冲波形;采用多个带通滤波器下的峰值信息作为特征向量,能够充分反映脉冲波形信息,提高了干扰信号与局放信号分离的精度;在对脉冲信号与局放信号进行分离时,无需建立局放以及噪声指纹库,结合聚类分离结果以及单类信号相位分布谱图特征,即可分辨脉冲干扰以及局部放电信号类型。
2.如权利要求1所述的用于局部放电信号的信噪分离方法,其特征在于,所述用于局部放电信号的信噪分离方法包括:
(1)采集局部放电波形信号;
(2)对局放信号进行提取,获取单个脉冲信号,记录脉冲的峰值及相位,并绘制相位分布谱图;
(3)将所提取脉冲通过多个带通滤波器,获得对应滤波器下脉冲的峰值信息,并通过主成分分析降维,获得特征参数;
(4)对特征参数进行聚类分析,将脉冲分为多类,并根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号。
3.如权利要求2所述的用于局部放电信号的信噪分离方法,其特征在于,所述对局放信号进行提取,获取单个脉冲信号,记录脉冲的峰值及相位,并绘制相位分布谱图具体包括:
步骤一,采集的局部放电波形信号为离散时间序列,将采集信号的幅值与设置的脉冲幅值阈值进行比较;若采集信号幅值大于所述脉冲幅值阈值,则存储该幅值以及所对应时间序列中的位置信息;否则,不进行存储;
步骤二,读取存储的采样点,将相邻采样点时间之差与设置的脉冲时间阈值进行比较,确定脉冲的第一个大于脉冲阈值的点P以及最后一个大于脉冲阈值的点Q,并分别进行存储;
步骤三,根据步骤二中存储的采样点,确定脉冲的起点以及终点,并存储整个脉冲的波形,取该脉冲幅值的最大值作为该脉冲的峰值;并将其所对应的时间序列进行转换确定局放相位,绘制相位分布谱图。
4.如权利要求3所述的用于局部放电信号的信噪分离方法,其特征在于,所述步骤三中确定脉冲起点以及终点方法具体实现如下:
第一步,计算背景噪声的均值以及方差,背景噪声可认为是服从于均值为μ,方差为σ2的正态分布的数字随机序列,均值以及方差计算公式:
式中,Xi为背景噪声序列,N为截取的背景噪声的点数;
第二步,计算背景噪声的时域能量,定义长度为M的信号序列时域能量如式;
第三步,确定脉冲的起始以及结束位置,分别在点P处进行前向搜索、点Q处进行后向搜索,计算长度为M的信号的时域能量;将时域能量与设置的能量阈值进行比较,当时域能量小于等于能量阈值时,完成搜索,确定脉冲的起始位置和结束位置。
5.如权利要求2所述的用于局部放电信号的信噪分离方法,其特征在于,所述将所提取脉冲通过多个带通滤波器,获得对应滤波器下脉冲的峰值信息,并通过主成分分析降维,获得特征参数具体包括以下步骤:
步骤一,对采集信号进行傅里叶变换得到频谱信号,获取采集信号的频率范围;
步骤二,在采集信号频率范围内,以0.1MHz为步进,计算提取的所有脉冲在相应频率下的方差;
步骤三,取方差较大的前三个值所在频率作为第一带通滤波器、第二带通滤波器、第三带通滤波器的中心频率fc1、fc2、fc3;沿方差曲线分别在所述中心频率两侧取3dB带宽作为第一带通滤波器、第二带通滤波器、第三带通滤波器的带宽BW1、BW2、BW3;
步骤四,将提取的单个脉冲信号分别通过三个带通滤波器,获得多个频带下脉冲的多个峰值,分别为第一峰值、第二峰值、第三峰值,并进行归一化处理;
步骤五,通过主成分分析降维算法对步骤四中所述的多个峰值进行降维,提取其前两阶X1、X2作为特征参数。
6.如权利要求2所述的用于局部放电信号的信噪分离方法,其特征在于,所述对特征参数进行聚类分析,将脉冲分为多类,并根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号具体包括:
步骤一,根据特征参数(X1,X2)绘制二维坐标散点图;
步骤二,采用模糊C均值聚类算法对特征参数进行聚类,得到分离的噪声信号与局放信号;
步骤三,根据分类类别号,使用不同图例在同一坐标系下绘制不同类别脉冲的二维坐标散点图;
步骤四,根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号。
7.如权利要求6所述的用于局部放电信号的信噪分离方法,其特征在于,采用聚类算法进行聚类时,聚类类别数有两种确定方式,一是人为确定,二是根据聚类有效性指标进行评分,计算聚类有效性函数,自适应确定最佳聚类类别数;
所述聚类有效性指标包括类内紧密性指标J以及类间分离性指标K,其中类内紧密性指标J表征类内样本的紧聚程度,类内紧密性指标J越小,类内样本的紧聚程度越高,各子类内样本点越相似,其计算公式如式:
类间分离性指标K表征类与类之间的分散程度,类间分离性指标K越大,类与类之间的分散程度越大,类间样本点的差异性越大,其计算公式如式:
其中,N为样本个数,c为聚类数,μij为第j个点与第i个分类簇之间的隶属关系;
聚类有效性函数f根据类内紧密性指标J以及类间分离性指标K,寻求紧致性最小且分离性最大的最佳聚类分类数c,其计算公式如式:
聚类有效性函数f的最小值对应最佳的聚类数c,即可自适应确定最佳聚类数。
8.一种执行权利要求1所述用于局部放电信号的信噪分离方法的用于局部放电信号的信噪分离***,其特征在于,所述用于局部放电信号的信噪分离***包括:
信号采集模块,用于采集局部放电波形信号;
脉冲提取模块,用于对所述局放信号进行提取,获取单个脉冲信号,记录脉冲的峰值及相位,并绘制相位分布谱图;
特征参数提取模块,用于将所提取脉冲通过多个带通滤波器,获得对应滤波器下脉冲的峰值信息,并进行主成分分析降维,获得特征参数;
聚类模块,用于对特征参数进行聚类分析,将脉冲分为多类,并根据分类类别号、脉冲峰值、脉冲相位绘制单类脉冲相位分布谱图,确定局放信号与干扰信号。
9.一种包含权利要求1-7任意一项所述用于局部放电信号的信噪分离方法的信息数据处理终端。
10.一种包含权利要求1-7任意一项所述用于局部放电信号的信噪分离方法的高压电气设备局部放电检测***。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910146996.XA CN109901031B (zh) | 2019-02-27 | 2019-02-27 | 用于局部放电信号的信噪分离方法、信息数据处理终端 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910146996.XA CN109901031B (zh) | 2019-02-27 | 2019-02-27 | 用于局部放电信号的信噪分离方法、信息数据处理终端 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109901031A true CN109901031A (zh) | 2019-06-18 |
CN109901031B CN109901031B (zh) | 2021-06-11 |
Family
ID=66945734
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910146996.XA Active CN109901031B (zh) | 2019-02-27 | 2019-02-27 | 用于局部放电信号的信噪分离方法、信息数据处理终端 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109901031B (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111103517A (zh) * | 2020-01-20 | 2020-05-05 | 云南电网有限责任公司玉溪供电局 | 一种真空度局部放电脉冲群分类识别方法 |
CN111537850A (zh) * | 2020-05-21 | 2020-08-14 | 北京传动联普科技有限公司 | 一种局放信号智能分离分类诊断的大数据识别管理*** |
CN111579938A (zh) * | 2020-04-14 | 2020-08-25 | 杭州柯林电气股份有限公司 | 基于多传感器的多源局部放电信号分离方法、监测*** |
CN112307969A (zh) * | 2020-10-30 | 2021-02-02 | 全球能源互联网欧洲研究院 | 一种脉冲信号的分类辨识方法、装置及计算机设备 |
CN112698170A (zh) * | 2021-01-18 | 2021-04-23 | 上海格鲁布科技有限公司 | 基于时域波形的特高频局部放电检测抗干扰方法 |
CN112763871A (zh) * | 2020-12-30 | 2021-05-07 | 珠海华网科技有限责任公司 | 一种局部放电分类识别方法 |
CN112834875A (zh) * | 2020-12-31 | 2021-05-25 | 康威通信技术股份有限公司 | 一种局部放电脉冲分割方法及*** |
CN113092972A (zh) * | 2021-06-07 | 2021-07-09 | 杭州兆华电子有限公司 | 基于非接触式超声波传感器的局放类型识别方法和装置 |
CN113221615A (zh) * | 2020-12-31 | 2021-08-06 | 中国石油化工股份有限公司 | 一种基于降噪聚类的局部放电脉冲提取方法 |
CN113341280A (zh) * | 2021-04-29 | 2021-09-03 | 西安交通大学 | 一种变压器油中局部放电在线定位监测***及方法 |
CN113406445A (zh) * | 2021-03-24 | 2021-09-17 | 上海朋禾智能科技有限公司 | 一种变压器局部放电多维数据在线检测方法 |
CN113433437A (zh) * | 2021-06-28 | 2021-09-24 | 华北电力大学 | 一种提取放电电流脉冲的方法及*** |
CN113484700A (zh) * | 2021-06-03 | 2021-10-08 | 武汉慧测电力科技有限公司 | 基于室内智能巡检机器人的开关柜局放检测方法 |
CN113837003A (zh) * | 2021-08-19 | 2021-12-24 | 深圳供电局有限公司 | 识别局部放电类型方法、装置、存储介质及电子装置 |
CN114065824A (zh) * | 2022-01-17 | 2022-02-18 | 陕西公众电气股份有限公司 | 一种对局部放电信号类型快速识别的方法 |
CN114545167A (zh) * | 2022-02-23 | 2022-05-27 | 四川大学 | 一种基于t-SNE算法的电缆终端局放脉冲分类方法 |
CN114626423A (zh) * | 2022-05-12 | 2022-06-14 | 杭州兆华电子股份有限公司 | 一种改进Kmeans的局部放电信号成像方法 |
CN114859194A (zh) * | 2022-07-07 | 2022-08-05 | 杭州兆华电子股份有限公司 | 一种基于非接触式的局部放电检测方法及装置 |
CN115166453A (zh) * | 2022-09-08 | 2022-10-11 | 国网智能电网研究院有限公司 | 基于边缘实时射频脉冲分类的局放连续监测方法及装置 |
CN115980528A (zh) * | 2023-02-14 | 2023-04-18 | 湖南湘能智能配电设备有限公司 | 一种环网柜噪音信号快速诊断与消除方法与装置 |
CN117892067A (zh) * | 2024-03-15 | 2024-04-16 | 国网上海市电力公司 | 一种低频电流局放监测抗干扰方法、装置、设备及介质 |
CN118094264A (zh) * | 2024-04-29 | 2024-05-28 | 山东泰开电力电子有限公司 | 一种智能化电力电容器局部放电检测方法及*** |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102645620A (zh) * | 2012-05-17 | 2012-08-22 | 广东电网公司电力科学研究院 | 基于时频特征参数的变电站多源局部放电检测方法及装置 |
CN103149513A (zh) * | 2013-02-26 | 2013-06-12 | 中国科学院电工研究所 | 一种变压器局部放电超声源重建的定位方法及装置 |
CN103197209A (zh) * | 2013-03-27 | 2013-07-10 | 西安电子科技大学 | 局部放电工频同步信号传感装置 |
CN103338023A (zh) * | 2013-03-19 | 2013-10-02 | 甘肃省电力公司检修公司 | 高频局放测量时消除与局放无关的外部噪声的滤波器 |
KR101320433B1 (ko) * | 2012-05-30 | 2013-10-23 | 주식회사 와이즈오토모티브 | 절연 및 부분 방전 신호 검출 장치 |
CN103675610A (zh) * | 2013-09-29 | 2014-03-26 | 国家电网公司 | 局部放电在线检测中的特征因子提取方法 |
CN104635126A (zh) * | 2015-01-27 | 2015-05-20 | 国家电网公司 | 一种基于滑动窗口的局部放电单脉冲提取方法 |
CN106896306A (zh) * | 2017-04-26 | 2017-06-27 | 国网上海市电力公司 | 一种gis振荡冲击耐压试验用信号去噪方法 |
US10024903B2 (en) * | 2013-10-29 | 2018-07-17 | Prysmian S.P.A. | Partial discharge detection system and method employing a synthetized synchronization signal |
CN108446632A (zh) * | 2018-03-20 | 2018-08-24 | 珠海华网科技有限责任公司 | 一种局部放电脉冲边沿寻找与局部放电确认方法 |
-
2019
- 2019-02-27 CN CN201910146996.XA patent/CN109901031B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102645620A (zh) * | 2012-05-17 | 2012-08-22 | 广东电网公司电力科学研究院 | 基于时频特征参数的变电站多源局部放电检测方法及装置 |
KR101320433B1 (ko) * | 2012-05-30 | 2013-10-23 | 주식회사 와이즈오토모티브 | 절연 및 부분 방전 신호 검출 장치 |
CN103149513A (zh) * | 2013-02-26 | 2013-06-12 | 中国科学院电工研究所 | 一种变压器局部放电超声源重建的定位方法及装置 |
CN103338023A (zh) * | 2013-03-19 | 2013-10-02 | 甘肃省电力公司检修公司 | 高频局放测量时消除与局放无关的外部噪声的滤波器 |
CN103197209A (zh) * | 2013-03-27 | 2013-07-10 | 西安电子科技大学 | 局部放电工频同步信号传感装置 |
CN103675610A (zh) * | 2013-09-29 | 2014-03-26 | 国家电网公司 | 局部放电在线检测中的特征因子提取方法 |
US10024903B2 (en) * | 2013-10-29 | 2018-07-17 | Prysmian S.P.A. | Partial discharge detection system and method employing a synthetized synchronization signal |
CN104635126A (zh) * | 2015-01-27 | 2015-05-20 | 国家电网公司 | 一种基于滑动窗口的局部放电单脉冲提取方法 |
CN106896306A (zh) * | 2017-04-26 | 2017-06-27 | 国网上海市电力公司 | 一种gis振荡冲击耐压试验用信号去噪方法 |
CN108446632A (zh) * | 2018-03-20 | 2018-08-24 | 珠海华网科技有限责任公司 | 一种局部放电脉冲边沿寻找与局部放电确认方法 |
Non-Patent Citations (3)
Title |
---|
RAMY HUSSEIN,KHALED BASHIR SHABAN,AYMAN H. EL-HAG: "Robust feature extraction and classification of acoustic partial discharge signals corrupted with noise", 《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》 * |
XIAOXIN CHEN,YONG QIAN,YONGPENG XU,GEHAO SHENG,XIUCHEN JIANG: "Partial discharge and noise separation in", 《IET GENERATION, TRANSMISSION & DISTRIBUTION》 * |
王刘旺,朱永利,李莉,贾亚飞: "基于自适应双阈值的局部放电基本参数提取", 《高电压技术》 * |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111103517A (zh) * | 2020-01-20 | 2020-05-05 | 云南电网有限责任公司玉溪供电局 | 一种真空度局部放电脉冲群分类识别方法 |
CN111579938A (zh) * | 2020-04-14 | 2020-08-25 | 杭州柯林电气股份有限公司 | 基于多传感器的多源局部放电信号分离方法、监测*** |
CN111579938B (zh) * | 2020-04-14 | 2024-03-22 | 杭州柯林电气股份有限公司 | 基于多传感器的多源局部放电信号分离方法、监测*** |
CN111537850A (zh) * | 2020-05-21 | 2020-08-14 | 北京传动联普科技有限公司 | 一种局放信号智能分离分类诊断的大数据识别管理*** |
CN111537850B (zh) * | 2020-05-21 | 2022-05-27 | 北京传动联普科技有限公司 | 一种局放信号智能分离分类诊断的大数据识别管理*** |
CN112307969A (zh) * | 2020-10-30 | 2021-02-02 | 全球能源互联网欧洲研究院 | 一种脉冲信号的分类辨识方法、装置及计算机设备 |
CN112307969B (zh) * | 2020-10-30 | 2024-03-12 | 全球能源互联网欧洲研究院 | 一种脉冲信号的分类辨识方法、装置及计算机设备 |
CN112763871A (zh) * | 2020-12-30 | 2021-05-07 | 珠海华网科技有限责任公司 | 一种局部放电分类识别方法 |
CN112834875A (zh) * | 2020-12-31 | 2021-05-25 | 康威通信技术股份有限公司 | 一种局部放电脉冲分割方法及*** |
CN113221615A (zh) * | 2020-12-31 | 2021-08-06 | 中国石油化工股份有限公司 | 一种基于降噪聚类的局部放电脉冲提取方法 |
CN112698170A (zh) * | 2021-01-18 | 2021-04-23 | 上海格鲁布科技有限公司 | 基于时域波形的特高频局部放电检测抗干扰方法 |
CN113406445A (zh) * | 2021-03-24 | 2021-09-17 | 上海朋禾智能科技有限公司 | 一种变压器局部放电多维数据在线检测方法 |
CN113341280A (zh) * | 2021-04-29 | 2021-09-03 | 西安交通大学 | 一种变压器油中局部放电在线定位监测***及方法 |
CN113484700A (zh) * | 2021-06-03 | 2021-10-08 | 武汉慧测电力科技有限公司 | 基于室内智能巡检机器人的开关柜局放检测方法 |
CN113092972B (zh) * | 2021-06-07 | 2021-09-03 | 杭州兆华电子有限公司 | 基于非接触式超声波传感器的局放类型识别方法和装置 |
CN113092972A (zh) * | 2021-06-07 | 2021-07-09 | 杭州兆华电子有限公司 | 基于非接触式超声波传感器的局放类型识别方法和装置 |
CN113433437A (zh) * | 2021-06-28 | 2021-09-24 | 华北电力大学 | 一种提取放电电流脉冲的方法及*** |
CN113837003A (zh) * | 2021-08-19 | 2021-12-24 | 深圳供电局有限公司 | 识别局部放电类型方法、装置、存储介质及电子装置 |
CN113837003B (zh) * | 2021-08-19 | 2024-05-14 | 深圳供电局有限公司 | 识别局部放电类型方法、装置、存储介质及电子装置 |
CN114065824A (zh) * | 2022-01-17 | 2022-02-18 | 陕西公众电气股份有限公司 | 一种对局部放电信号类型快速识别的方法 |
CN114545167A (zh) * | 2022-02-23 | 2022-05-27 | 四川大学 | 一种基于t-SNE算法的电缆终端局放脉冲分类方法 |
CN114545167B (zh) * | 2022-02-23 | 2023-06-16 | 四川大学 | 一种基于t-SNE算法的电缆终端局放脉冲分类方法 |
CN114626423A (zh) * | 2022-05-12 | 2022-06-14 | 杭州兆华电子股份有限公司 | 一种改进Kmeans的局部放电信号成像方法 |
CN114859194B (zh) * | 2022-07-07 | 2022-09-23 | 杭州兆华电子股份有限公司 | 一种基于非接触式的局部放电检测方法及装置 |
CN114859194A (zh) * | 2022-07-07 | 2022-08-05 | 杭州兆华电子股份有限公司 | 一种基于非接触式的局部放电检测方法及装置 |
CN115166453A (zh) * | 2022-09-08 | 2022-10-11 | 国网智能电网研究院有限公司 | 基于边缘实时射频脉冲分类的局放连续监测方法及装置 |
CN115166453B (zh) * | 2022-09-08 | 2023-01-24 | 国网智能电网研究院有限公司 | 基于边缘实时射频脉冲分类的局放连续监测方法及装置 |
CN115980528A (zh) * | 2023-02-14 | 2023-04-18 | 湖南湘能智能配电设备有限公司 | 一种环网柜噪音信号快速诊断与消除方法与装置 |
CN115980528B (zh) * | 2023-02-14 | 2024-03-26 | 湖南湘能智能配电设备有限公司 | 一种环网柜噪音信号快速诊断与消除方法与装置 |
CN117892067A (zh) * | 2024-03-15 | 2024-04-16 | 国网上海市电力公司 | 一种低频电流局放监测抗干扰方法、装置、设备及介质 |
CN117892067B (zh) * | 2024-03-15 | 2024-05-28 | 国网上海市电力公司 | 一种低频电流局放监测抗干扰方法、装置、设备及介质 |
CN118094264A (zh) * | 2024-04-29 | 2024-05-28 | 山东泰开电力电子有限公司 | 一种智能化电力电容器局部放电检测方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN109901031B (zh) | 2021-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109901031A (zh) | 用于局部放电信号的信噪分离方法、信息数据处理终端 | |
CN106405339B (zh) | 基于高低频小波特征关联的输电线路故障原因辨识方法 | |
CN110175508B (zh) | 一种应用于超声波局放检测的特征值提取方法 | |
CN108469560B (zh) | 一种基于快速s变换时频空间模型的电磁干扰客观复杂度评估方法 | |
Seo et al. | Probabilistic wavelet transform for partial discharge measurement of transformer | |
Contin et al. | Classification and separation of partial discharge signals by means of their auto-correlation function evaluation | |
CN110247420B (zh) | 一种hvdc输电线路故障智能识别方法 | |
Firuzi et al. | A novel method for differentiating and clustering multiple partial discharge sources using S transform and bag of words feature | |
CN104198898A (zh) | 基于脉冲序列分析的局部放电发展过程诊断方法 | |
CN109856517A (zh) | 一种特高压设备局部放电检测数据的判别方法 | |
CN105701470A (zh) | 一种基于最优小波包分解的模拟电路故障特征提取方法 | |
CN103197218A (zh) | 一种高压电缆绝缘缺陷局放带电检测诊断方法 | |
CN106503612A (zh) | 基于压缩感知理论的直流电缆局部放电信号模式识别方法及*** | |
CN114325256A (zh) | 一种电力设备局部放电识别方法、***、设备及存储介质 | |
Yang et al. | Ship recognition via its radiated sound: The fractal based approaches | |
CN110176250A (zh) | 一种基于局部学习的鲁棒声学场景识别方法 | |
CN105893976B (zh) | 一种行波信号的参数化识别方法 | |
Malik et al. | Automatic threshold optimization in nonlinear energy operator based spike detection | |
CN115728612A (zh) | 变压器放电故障诊断方法及装置 | |
Du et al. | Feature extraction and pattern recognition algorithm of power cable partial discharge signal | |
Dang et al. | Fault diagnosis of power transformer by acoustic signals with deep learning | |
Seo et al. | A novel signal extraction technique for online partial discharge (PD) measurement of transformers | |
Ghosh et al. | A method for unambiguous identification of on-field recorded insulator leakage current waveforms portraying electrical activity on the surface | |
CN112581940A (zh) | 基于边缘计算与神经网络的放电声音检测方法 | |
Lu et al. | Convolutional neural network based transient earth voltage detection |
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 |