CN114609483A - 一种基于Hilbert变换的GIS局部放电信号特征提取方法 - Google Patents
一种基于Hilbert变换的GIS局部放电信号特征提取方法 Download PDFInfo
- Publication number
- CN114609483A CN114609483A CN202210073612.8A CN202210073612A CN114609483A CN 114609483 A CN114609483 A CN 114609483A CN 202210073612 A CN202210073612 A CN 202210073612A CN 114609483 A CN114609483 A CN 114609483A
- Authority
- CN
- China
- Prior art keywords
- partial discharge
- characteristic
- signal
- feature
- 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.)
- Pending
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 35
- 238000012847 principal component analysis method Methods 0.000 claims abstract description 6
- 239000013598 vector Substances 0.000 claims description 113
- 230000009467 reduction Effects 0.000 claims description 34
- 238000001228 spectrum Methods 0.000 claims description 28
- 238000004422 calculation algorithm Methods 0.000 claims description 21
- 230000000694 effects Effects 0.000 claims description 19
- 238000000513 principal component analysis Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 239000006185 dispersion Substances 0.000 claims description 7
- 230000005484 gravity Effects 0.000 claims description 7
- 230000002238 attenuated effect Effects 0.000 claims description 6
- 238000012549 training Methods 0.000 claims description 6
- 230000002776 aggregation Effects 0.000 claims description 5
- 238000004220 aggregation Methods 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000009413 insulation Methods 0.000 abstract description 8
- 239000000284 extract Substances 0.000 abstract description 5
- 230000006866 deterioration Effects 0.000 abstract 1
- 230000002349 favourable effect Effects 0.000 abstract 1
- 230000000875 corresponding effect Effects 0.000 description 28
- 230000007547 defect Effects 0.000 description 7
- 230000000630 rising effect Effects 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008092 positive effect Effects 0.000 description 2
- 238000011946 reduction process Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000002923 metal particle Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
- G01R31/1227—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials
- G01R31/1254—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials of gas-insulated power appliances or vacuum gaps
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- 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
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
- G06F2218/16—Classification; Matching by matching signal segments
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Discrete Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Testing Relating To Insulation (AREA)
Abstract
本发明请求保护一种基于Hilbert变换的GIS局部放电信号特征提取方法。该方法利用Hilbert变换求取GIS局部放电信号的包络,并基于包络提取局部放电脉冲波形的统计特征参数、时域特征参数和频域特征参数,采用改进的核主成分分析法对提取的特征参数进行降维,最终得到维数较低的且能正确执行局部放电分类的特征量。相比通常的GIS局部放电信号特征提取方法,本发明提取的局部放电特征量不仅维数较低,而且物理意义更为明确,更能准确地表征局部放电信号。本发明对GIS局部放电信号特征进行有效提取,有利于对GIS内产生的局部放电进行检测识别,及时了解设备内部的绝缘劣化情况,对于GIS设备的安全运行具有重要价值。
Description
技术领域
本发明GIS局部放电在线监测技术领域,具体涉及一种基于Hilbert变换的GIS局部放电信号特征提取方法。
背景技术
随着我国国民经济的快速发展,电力工业建设的突飞猛进,大电网、大容量、超高压、大机组己经成为现代电力***的特点和发展趋势,为保障电力***的安全性、稳定性和可靠性,对电力设备运行提出了更高的要求。由于气体绝缘组合电器(Gas InsulatedSwitchgear,简称GIS)具有占地面积小、运行安全可靠、维护工作量少、检修周期长等一系列优点,近些年来在现代电力***中得到了广泛的应用。
GIS在制造、运输、安装和运行的过程中容易出现金属毛刺、电极电位浮动、自由金属颗粒和绝缘老化等各种缺陷,这些缺陷将导致不同程度的局部放电,长期的局部放电将使设备绝缘进一步劣化,严重时将造成整个绝缘击穿,从而对设备的安全运行造成严重的威胁。因此,有必要对GIS内产生的局部放电进行检测识别,判断引发局部放电的绝缘缺陷类型,及时针对不同的缺陷类型进行相应的检修维护工作,防止出现绝缘整体击穿的重大事故。其中局部放电信号特征的提取是局部放电信号检测识别的关键。
当前国内外学者对局部放电信号特征的提取开展了大量研究,主要形成了以下四类特征提取方法:(1)统计特征提取法;(2)分形特征提取法;(3)小波特征提取法;(4)图像矩特征提取法,这四类特征提取方法主要存在以下问题:
统计特征提取法所需的样本数量较大,得到的特征量维数较高,会产生一定的信息冗余,不利于后续的模式识别;分形特征提取法的分形维数计算容易受到局部放电信号长度和有效信号点数的影响,并且提取的特征量维数较高;小波特征提取法在计算特征量时很难选择合适的小波基和分解层数;图像矩特征法计算复杂度较高,且没有明确的物理意义,在实际应用中受到一定限制。
因此目前急需一种能简单、有效提取GIS局部放电信号特征的方法,以便对GIS内产生的局部放电信号进行快速、有效地识别。
经过检索,申请公开号为CN110231548A,一种基于超声信号包络谱的GIS局部放电特征提取方法,该方法可以帮助识别不同类型的局部放电绝缘缺陷。包括:设计三种GIS典型绝缘缺陷,采集局部放电超声信号;利用FIR带通滤波器对局放超声信号进行滤波处理;利用Hilbert变换对滤波后的局放超声信号进行包络提取;对信号包络进行工频相关性及相位模式分析,得到不同缺陷类型的测量特征,
上述发明基于超声信号包络谱只提取了局部放电信号的峰峰值、峰值、有效值、偏斜度和峭度等波形统计特征参数,却忽略了方差、波形因子、峰值因子、脉冲因子和裕度因子等包含局部放电信号重要波形特征信息的特征参数;上述发明基于超声信号包络谱只提取了局部放电信号的50Hz频率相关性和100Hz频率相关性等频域特征参数,却忽略了频率幅值均值、重心频率、均方根频率、频率方差和第一主频频率等包含局部放电信号重要频域特征信息的特征参数。因此上述发明提取的局部放电特征没有包含局部放电信号的大部分信息,并不能很好地表征局部放电信号,会影响后续局部放电模式识别的精度。
本发明分别基于信号的时域波形包络和频谱包络来提取时域特征参数和频域特征参数,并联合由局部放电有效数据提取得到的统计特征参数共同构建放电特征数据,由上述三类特征参数构建的局部放电特征数据总共包含了24组信号特征指标,基本包含了局部放电信号的全部特征信息,解决了前述发明存在的局部放电特征包含的信息量不足等问题。同时本发明还创新性地提出了改进的KPCA算法,并利用改进的KPCA算法对提取的局部放电特征数据进行降维处理,得到了维数较低且能正确执行局部放电分类的特征量,其较低的特征维数也使得后续识别模型的训练和识别时间大大缩减,提高了局部放电模式识别效率
因此相比前述发明,本发明提取的局部放电特征参数包含更多的有效放电信息,且特征维数较低,分类效果更好,更有利于后续的局部放电模式识别。
发明内容
本发明旨在解决以上现有技术的问题。提出了一种能简单、有效地基于Hilbert变换的GIS局部放电信号特征提取方法。本发明的技术方案如下:
一种基于Hilbert变换的GIS局部放电信号特征提取方法,其包括以下步骤:
采集现场检测的气体绝缘组合电器GIS局部放电信号;
对采集的GIS局部放电信号进行有效数据的提取,去除不包含放电信息的无效信号数据;
利用Hilbert变换求取放电信号的包络;
基于放电信号的包络提取局部放电脉冲波形的统计特征参数、时域特征参数和频域特征参数,利用这三类特征参数构建局部放电特征数据;
利用改进的核主成分分析法对提取的局部放电特征数据进行降维处理,核主成分分析KPCA是通过一个非线性映射将原始输入空间的特征数据投影到高维特征空间,然后在高维特征空间中利用主成分分析PCA对数据进行降维,最终得到维数较低的且能正确执行局部放电分类的特征数据。
进一步的,所述采集现场检测的气体绝缘组合电器GIS局部放电信号是利用Keysight品牌的DSOX6002A型号示波器进行采集,其带宽为和采样率最高分别可达4GHz和20.0GSa/s,可以实现对GIS局部放电信号的有效采集。
进一步的,所述对局部放电信号进行有效数据的提取,具体包括:只截取局部放电信号中间的一段放电脉冲数据进行特征参数的提取,将放电脉冲出现前和放电脉冲衰减完后的数据进行去除。
进一步的,所述利用Hilbert变换求取放电信号的包络,具体包括:
首先对局部放电信号时域波形的电压值取绝对值,得到局部放电信号的单极性波形,然后利用Hilbert变换求取该单极性波形的包络,得到局部放电信号的时域波形包络;在完成时域波形包络的提取后,对放电信号波形进行快速傅里叶变换,将信号的波形由时域变换到频域中,然后利用Hilbert变换求取放电信号的频谱包络。
进一步的,所述Hilbert变换是数学和信号处理的一种同域变换,提供了实现给定函数或傅里叶序列的调和共轭的具体方法,对于连续时域信号x(t),定义其Hilbert变换为:
其中t是时间自变量,是用于作卷积运算的时间变换量。x(t)的解析信号a(t)定义为
其中j是复数的虚部单位。解析信号的模值E(t)即为x(t)的包络
而对于长度为m的离散信号序列x(n)及其FFT序列X(k),有
其中A(k)为x(n)的离散解析信号a(n)对应的FFT序列,则x(n)的包络E(n)为
E(n)=|a(n)|=|IFFT[A(k)]| (5)
IFFT表示快速傅里叶逆变换。
进一步的,所述基于信号的包络提取局部放电脉冲波形的统计特征参数、时域特征参数和频域特征参数,具体包括:
分别基于信号的时域波形包络和频谱包络来提取时域特征参数和频域特征参数,并联合由局部放电有效数据提取得到的统计特征参数共同构建放电特征数据;其中时域特征参数主要包括上升时间、下降时间、峰值时间和脉冲宽度,频域特征参数主要包括频率幅值均值、重心频率、均方根频率、频率方差、第一主频频率和第一主频峰值,统计特征参数主要包括均值、最大值、最小值、峰峰值、绝对均值、方根幅值、均方根值、方差、偏斜度、峭度、波形因子、峰值因子、脉冲因子和裕度因子;因此由上述三类特征参数构建的局部放电特征数据包含了24组信号特征指标,其特征维数为24维,几乎包含了局部放电信号的全部特征信息。
进一步的,构建局部放电特征数据的部分特征参数定义如下:
偏斜度:偏斜度α主要反映了信号的不对称情况,当偏斜度值越大,不对称性越严重,xi是信号第i个数据点的幅值,N为数据点的数目;
峭度:峭度β主要反映了信号的陡峭程度,对信号中的脉冲比较敏感;
波形因子:波形因子Sf是信号的均方根值与绝对均值的比值,其敏感性不好,但是稳定性较好;
峰值因子:峰值因子Cf是信号的最大值与均方根值的比值,反映了峰值在波形中的极端程度,xmax为信号的幅值最大值;
裕度因子:裕度因子CLf是信号的最大值与方根幅值的比值;
脉冲因子:脉冲因子If是信号的最大值与绝对均值的比值;
上升时间:上升时间tr为脉冲从峰值的10%处上升到峰值的90%处所需的时间,tr_90%和tr_10%分别为脉冲上升沿中峰值的90%和10%处所对应的时刻;
tr=tr_90%-tr_10% (12)
下降时间:下降时间td为脉冲从峰值的90%处下降到峰值的10%处所需的时间,td_90%和td_10%分别为脉冲下降沿中峰值的90%处和10%处所对应的时刻;
td=td_10%-td_90% (13)
峰值时间:峰值时间tp为脉冲从起始处上升到峰值处所需的时间,tpeak和t0分别为脉冲峰值处和起始处所对应的时刻;
tp=tpeak-t0 (14)
脉冲宽度:脉冲宽度tw为脉冲从上升沿中峰值的50%处到下降沿中峰值的50%处所需的时间,tr_50%和td_50%分别为脉冲上升沿中峰值的50%处和下降沿中峰值的50%处所对应的时刻;
tw=td_50%-tr_50% (15)
重心频率:重心频率FC能够描述信号在频谱中分量较大的信号成分的频率,反映信号功率谱的分布情况,fi是时域信号经快速傅里叶变换后产生的频谱中第i个信号成分的频率,ai为第i个信号成分对应的幅值,K为频谱中信号成分的数目;
频率方差:频率方差VF描述信号功率谱能量分布的分散程度;
第一主频频率:第一主频频率fi为频谱中幅值最大的信号成分对应的频率。
第一主频峰值:第一主频峰值为频谱中所有信号成分的最大的幅值。
进一步的,所述利用改进的KPCA算法对提取的局部放电特征数据进行降维处理,最终得到维数较低的且能正确执行局部放电分类的特征数据,具体包括:
一般的KPCA算法步骤如下:
对于一般PCA方法,即通过求解特征方程
λv=Cv (19)
获得方差贡献率大的特征值(对应较大的特征值)及与之对应的特征向量。现引入非线性映射函数Φ,使输入空间中的特征数据x1,x2,...,xM变换为高维特征空间中的特征数据Φ(x1),Φ(x2),...,Φ(xM),并假设
因此,高维特征空间中的PCA是求解方程
中的特征值λ和特征向量ν,进而有
注意到上式中特征向量ν可以由Φ(xk)(k=1,2,...,M)线性表示出,即
其中αi为Φ(xi)的系数,由式(22)~式(24)得
定义M×M矩阵K:
Kij≡Φ(xi)·Φ(xj) (26)
式(25)简化为
MλKα=K2α (27)
其中α是M×1的向量,第j个元素为αj,显然满足
Mλα=Kα (28)
则必然满足式(27),通过对式(28)的求解,以及联合式(24)即可获得方差贡献率大的特征值及与之对应的特征向量;
式中1ij=1(对于所有的i,j);
通过分别计算高维特征空间中的特征数据Φ(x1),Φ(x2),...,Φ(xM)在前n个最大的特征值λ1,λ2,...,λn所对应的特征向量ν1,ν2,…,νn上的投影值,即可得到降维后维数为n的特征数据。
上述一般的KPCA算法通过使用方差贡献率来选择特征向量,从而完成对高维特征空间中的特征数据的降维,具有较好的数据降维效果,但是在整个降维过程中只考虑最大化地保留样本的总特征信息量(特征值代表了对应特征向量方向上包含的信息量)而没有考虑样本的类别信息,即只从数学角度考虑矩阵主要成分而没有考虑实际样本的分类问题,这会导致降维后的特征量对于异常数据的分类并没有很好的效果。
而对于分类问题来说,特征提取的目的就是要提取出有利于后续分类的特征组合,为此本发明提出了一种改进的KPCA算法,通过各特征向量的类内聚集程度和类间离散程度来度量类别信息,不仅保留了较好的降维效果,而且使得降维后的特征数据更有利于后续的模式分类。每个特征向量的类内聚集程度描述为:
根据公式(32)可知,Jk越小就说明类别信息量越大,那么就越容易进行分类,即这个特征向量方向上包含的类别信息量越大。但并不是类别信息量越大就越能提高分类准确率,有些特征向量方向上包含的类别信息量大,但不一定会对分类效果起到积极的作用。因此,本发明使用相似度来衡量提取的各特征向量与样本的类别向量之间的相关性,通过特征向量与样本类别向量间的相关系数来选取特征向量的组合,去除对分类效果不明显的特征向量,其公式为:
式中:ξ为特征向量与样本类别属性向量的相关系数;vi为特征向量(其中i=1,2,3,…,n);vz是原样本中的类别属性向量。若ξ的绝对值大于零或等于1,则说明两特征向量是相关的;若ξ相似趋于或等于0,则说明两向量是不相关的。通过计算特征向量与样本的类别向量间的相关系数ξ,排除不相关或相关系数非常小的特征向量,去掉后会减小后续分类器的计算量,从而减小了分类的训练和识别时间。
改进的KPCA的具体降维步骤总结如下:
假设输入空间中有M组特征数据,每组特征数据有N个特征属性(不包含类别特征),根据改进的KPCA原理,特征提取的具体实现过程为:
(1)列出一个(M×N)维数据矩阵,其中M为样本个数,N为每个样本的属性个数;
(3)求取各特征向量的信息度量Jk,按照信息度量Jk从小到大的顺序,对相应的特征向量进行排序;
(4)提取前m个特征向量ν1,ν2,…,νm,从第一个特征向量ν1开始,依次根据式(33)计算ν1,ν2,...,νm与原样本中的类别属性向量vz的相关系数ξ;
(5)从步骤(4)中找到与类别向量vz的相关系数ξ为0或趋于0的特征向量,将其从特征向量组中去除,得到剩余的n个特征向量;
(6)对剩余特征向量按照相关系数从大到小的顺序重新排列;
(7)通过分别计算高维特征空间中的特征数据在剩余的特征向量上的投影值,即可得到降维后维数为n的且分类效果较好的特征数据。
本发明的优点及有益效果如下:
本发明的创新主要是提出了一种能简单、有效地提取GIS局部放电信号特征的方法。通常的局部放电特征提取方法(例如小波特征提取法、图像矩特征提取法等)存在计算复杂度高、计算参数设置困难和物理意义不明确等问题,而本发明采用基于Hilbert变换的特征提取方法,通过局部放电信号包络来提取放电特征量,不仅计算复杂度低,而且提取的特征量物理意义明确,基本包含了局部放电的全部特征信息,能准确地表征局部放电信号;通常的局部放电特征提取方法存在提取的特征量维数较高且分类效果不理想等问题,而本发明利用改进的KPCA算法对提取的局部放电特征数据进行降维处理,得到了维数较低且能正确执行局部放电分类的特征量,其较低的特征维数也使得后续识别模型的训练和识别时间大大缩减,提高了局部放电模式识别效率;相比一般的KPCA算法通过使用方差贡献率来选择特征向量,改进的KPCA算法通过求取各特征向量的信息度量以及各特征向量与原样本中类别属性向量的相关系数来选择特征向量,不仅保留了一般KPCA算法的良好降维效果,而且其降维后的特征数据分类效果更好。
相比前述发明(申请公开号为CN110231548A,一种基于超声信号包络谱的GIS局部放电特征提取方法),本发明提取的特征参数更多,分别提取了统计特征参数、时域特征参数和频域特征参数等共计24组特征参数来构建局部放电特征数据,基本包含了局部放电信号的全部特征信息,解决了前述发明存在的局部放电特征包含的信息量不足等问题。同时本发明还创新性地提出了改进的KPCA算法,并利用改进的KPCA算法对提取的局部放电特征数据进行降维处理,得到了维数较低且能正确执行局部放电分类的特征量,其较低的特征维数也使得后续识别模型的训练和识别时间大大缩减,提高了局部放电模式识别效率。
因此相比前述发明,本发明提取的局部放电特征参数包含更多的有效放电信息,且特征维数较低,分类效果更好,更有利于后续的局部放电模式识别。
附图说明
图1是本发明提供优选实施例的步骤流程示意图;
图2是现场测得的原始GIS局部放电信号波形图;
图3是去除无效信号数据后的GIS局部放电信号波形图;
图4是提取的放电信号时域波形包络;
图5是提取的放电信号频谱包络。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、详细地描述。所描述的实施例仅仅是本发明的一部分实施例。
本发明解决上述技术问题的技术方案是:
如图1所示,一种基于Hilbert变换的GIS局部放电信号特征提取方法,其包括以下步骤:
1)利用局部放电信号采集单元采集现场检测的GIS局部放电信号。
2)对采集的GIS局部放电信号进行有效数据的提取,去除不包含放电信息的无效信号数据。
3)利用Hilbert变换求取放电信号的包络。
4)基于信号的包络提取局部放电脉冲波形的统计特征参数、时域特征参数和频域特征参数,利用这三类特征参数构建局部放电特征数据。
5)利用改进的核主成分分析法对提取的局部放电特征数据进行降维处理,最终得到维数较低的且能正确执行局部放电分类的特征数据。
进一步的,所述步骤2)对局部放电信号进行有效数据的提取,由于GIS局部放电脉冲出现前和衰减完后的数千个采样数据点基本不包含有效的放电信息,因此只截取信号中间的一段放电脉冲数据进行特征参数的提取,将放电脉冲出现前和衰减完后的数据进行去除。
进一步的,所述步骤3)中首先对局部放电信号时域波形的电压值取绝对值,得到局部放电信号的单极性波形,然后利用Hilbert变换求取该单极性波形的包络,得到局部放电信号的时域波形包络;在完成时域波形包络的提取后,对放电信号波形进行快速傅里叶变换,将信号的波形由时域变换到频域中,然后利用Hilbert变换求取放电信号的频谱包络;Hilbert变换是数学和信号处理的一种同域变换,提供了实现给定函数或傅里叶序列的调和共轭的具体方法,对于连续时域信号x(t),定义其Hilbert变换为:
其中t是时间自变量,τ是用于作卷积运算的时间变换量。
x(t)的解析信号a(t)定义为
其中j是复数的虚部单位。
解析信号的模值E(t)即为x(t)的包络
而对于长度为m的离散信号序列x(n)及其FFT序列X(k),有
其中A(k)为x(n)的离散解析信号a(n)对应的FFT序列。则x(n)的包络E(n)为
E(n)=|a(n)|=|IFFT[A(k)]| (5)
IFFT表示快速傅里叶逆变换。
进一步的,所述步骤4)中分别基于信号的时域波形包络和频谱包络来提取时域特征参数和频域特征参数,并联合统计特征参数(由局部放电有效数据提取得到)共同构建放电特征数据;其中时域特征参数主要包括上升时间、下降时间、峰值时间和脉冲宽度,频域特征参数主要包括频率幅值均值、重心频率、均方根频率、频率方差、第一主频频率和第一主频峰值,统计特征参数主要包括均值、最大值、最小值、峰峰值、绝对均值、方根幅值、均方根值、方差、偏斜度、峭度、波形因子、峰值因子、脉冲因子和裕度因子;因此由上述三类特征参数构建的局部放电特征数据包含了24组信号特征指标,其特征维数为24维,几乎包含了局部放电信号的全部特征信息;部分特征参数的定义如下:
偏斜度:偏斜度α主要反映了信号的不对称情况,当偏斜度值越大,不对称性越严重(xi是信号第i个数据点的幅值,N为数据点的数目)。
峭度:峭度β主要反映了信号的陡峭程度,对信号中的脉冲比较敏感。
波形因子:波形因子Sf是信号的均方根值与绝对均值的比值,其敏感性不好,但是稳定性较好。
峰值因子:峰值因子Cf是信号的最大值与均方根值的比值,反映了峰值在波形中的极端程度(xmax为信号的幅值最大值)。
裕度因子:裕度因子CLf是信号的最大值与方根幅值的比值,其敏感性很好,但稳定性一般。
脉冲因子:脉冲因子If是信号的最大值与绝对均值的比值,其拥有较好的敏感性,但是稳定性一般
上升时间:上升时间tr为脉冲从峰值的10%处上升到峰值的90%处所需的时间,tr_90%和tr_10%分别为脉冲上升沿中峰值的90%和10%处所对应的时刻。
tr=tr_90%-tr_10% (12)
下降时间:下降时间td为脉冲从峰值的90%处下降到峰值的10%处所需的时间,td_90%和td_10%分别为脉冲下降沿中峰值的90%处和10%处所对应的时刻。
td=td_10%-td_90% (13)
峰值时间:峰值时间tp为脉冲从起始处上升到峰值处所需的时间,tpeak和t0分别为脉冲峰值处和起始处所对应的时刻。
tp=tpeak-t0 (14)
脉冲宽度:脉冲宽度tw为脉冲从上升沿中峰值的50%处到下降沿中峰值的50%处所需的时间,tr_50%和td_50%分别为脉冲上升沿中峰值的50%处和下降沿中峰值的50%处所对应的时刻。
tw=td_50%-tr_50% (15)
重心频率:重心频率FC能够描述信号在频谱中分量较大的信号成分的频率,反映信号功率谱的分布情况(fi是时域信号经快速傅里叶变换后产生的频谱中第i个信号成分的频率,ai为第i个信号成分对应的幅值,K为频谱中信号成分的数目)。
频率方差:频率方差VF描述信号功率谱能量分布的分散程度。
第一主频频率:第一主频频率fi为频谱中幅值最大的信号成分对应的频率。
第一主频峰值:第一主频峰值为频谱中所有信号成分的最大的幅值
进一步的,所述步骤5)中利用改进的KPCA算法对提取的局部放电特征数据进行降维处理,最终得到维数较低的且能正确执行局部放电分类的特征数据;改进的KPCA算法降维的主要原理如下:核主成分分析(KPCA)是一种非线性数据处理方法,其核心思想是通过一个非线性映射将原始输入空间的特征数据投影到高维特征空间,然后在高维特征空间中利用主成分分析(PCA)对数据进行降维。
一般的KPCA算法步骤如下:
对于一般PCA方法,即通过求解特征方程
λv=Cv (19)
获得方差贡献率大的特征值(对应较大的特征值)及与之对应的特征向量。现引入非线性映射函数Φ,使输入空间中的特征数据x1,x2,...,xM变换为高维特征空间中的特征数据Φ(x1),Φ(x2),...,Φ(xM),并假设
因此,高维特征空间中的PCA是求解方程
中的特征值λ和特征向量ν,进而有
注意到上式中特征向量ν可以由Φ(xk)(k=1,2,...,M)线性表示出,即
其中αi为Φ(xi)的系数,由式(22)~式(24)得
定义M×M矩阵K:
Kij≡Φ(xi)·Φ(xj) (26)
式(25)简化为
MλKα=K2α (27)
其中α是M×1的向量,第j个元素为αj,显然满足
Mλα=Kα (28)
则必然满足式(27)。通过对式(28)的求解,以及联合式(24)即可获得方差贡献率大的特征值(对应较大的特征值)及与之对应的特征向量。
式中1ij=1(对于所有的i,j)。
通过分别计算高维特征空间中的特征数据Φ(x1),Φ(x2),...,Φ(xM)在前n个最大的特征值λ1,λ2,...,λn(前n个最大特征值的累计方差贡献率应达到90%以上)所对应的特征向量ν1,ν2,…,νn上的投影值,即可得到降维后维数为n的特征数据。
上述一般的KPCA算法通过使用方差贡献率来选择特征向量,从而完成对高维特征空间中的特征数据的降维,具有较好的数据降维效果,但是在整个降维过程中只考虑最大化地保留样本的总特征信息量(特征值代表了对应特征向量方向上包含的信息量)而没有考虑样本的类别信息,即只从数学角度考虑矩阵主要成分而没有考虑实际样本的分类问题,这会导致降维后的特征量对于异常数据的分类并没有很好的效果。
而对于分类问题来说,特征提取的目的就是要提取出有利于后续分类的特征组合,为此本发明提出了一种改进的KPCA算法,通过各特征向量的类内聚集程度和类间离散程度来度量类别信息,不仅保留了较好的降维效果,而且使得降维后的特征数据更有利于后续的模式分类。每个特征向量的类内聚集程度描述为:
根据公式(32)可知,Jk越小就说明类别信息量越大,那么就越容易进行分类,即这个特征向量方向上包含的类别信息量越大。但并不是类别信息量越大就越能提高分类准确率,有些特征向量方向上包含的类别信息量大,但不一定会对分类效果起到积极的作用。因此,本发明使用相似度来衡量提取的各特征向量与样本的类别向量之间的相关性,通过特征向量与样本类别向量间的相关系数来选取特征向量的组合,去除对分类效果不明显的特征向量,其公式为:
式中:ξ为特征向量与样本类别属性向量的相关系数;vi为特征向量(其中i=1,2,3,…,n);vz是原样本中的类别属性向量。若ξ的绝对值大于零或等于1,则说明两特征向量是相关的;若ξ相似趋于或等于0,则说明两向量是不相关的。通过计算特征向量与样本的类别向量间的相关系数ξ,排除不相关或相关系数非常小的特征向量,去掉后会减小后续分类器的计算量,从而减小了分类的训练和识别时间。
改进的KPCA的具体降维步骤总结如下:
假设输入空间中有M组特征数据,每组特征数据有N个特征属性(不包含类别特征),根据改进的KPCA原理,特征提取的具体实现过程为:
(1)列出一个(M×N)维数据矩阵,其中M为样本个数,N为每个样本的属性个数;
(3)求取各特征向量的信息度量Jk,按照信息度量Jk从小到大的顺序,对相应的特征向量进行排序;
(4)提取前m个特征向量ν1,ν2,...,νm,从第一个特征向量ν1开始,依次根据式(33)计算ν1,ν2,...,νm与原样本中的类别属性向量vz的相关系数ξ;
(5)从步骤(4)中找到与类别向量vz的相关系数ξ为0或趋于0的特征向量,将其从特征向量组中去除,得到剩余的n个特征向量;
(6)对剩余特征向量按照相关系数从大到小的顺序重新排列;
(7)通过分别计算高维特征空间中的特征数据在剩余的特征向量上的投影值,即可得到降维后维数为n的且分类效果较好的特征数据。
本发明的具体实现过程如图1所示。首先利用局部放电信号采集单元采集现场检测的GIS局部放电信号;然后对采集的GIS局部放电信号进行有效数据的提取,由于GIS局部放电脉冲出现前和衰减完后的数千个采样数据点基本不包含有效的放电信息,因此只截取信号中间的一段放电脉冲数据进行特征参数的提取,将放电脉冲出现前和衰减完后的数据进行去除;接着对局部放电信号时域波形的电压值取绝对值,得到局部放电信号的单极性波形,然后利用Hilbert变换求取该单极性波形的包络,得到局部放电信号的时域波形包络;在完成时域波形包络的提取后,对放电信号波形进行快速傅里叶变换,将信号的波形由时域变换到频域中,然后利用Hilbert变换求取放电信号的频谱包络;基于信号的时域波形包络和频谱包络来提取时域特征参数和频域特征参数,并联合统计特征参数共同构建放电特征数据;最后利用改进的核主成分分析法对提取的局部放电特征数据进行降维处理,最终得到维数较低的且能正确执行局部放电分类的特征数据。
以某次现场试验测得的GIS局部放电信号的特征提取为例,其原始数据如图2所示。对该放电信号进行有效数据的提取,去除不包含放电信息的无效信号数据,结果如图3所示。利用Hilbert变换求取该放电信号的时域波形包络和频谱包络,结果如图4、图5所示。基于信号的时域波形包络和频谱包络来提取时域特征参数和频域特征参数,并联合统计特征参数共同构建放电特征数据,利用改进的核主成分分析法对提取的放电特征数据进行降维处理,最终得到维数较低的且能正确执行局部放电分类的特征数据。
上述实施例阐明的***、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的,计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上这些实施例应理解为仅用于说明本发明而不用于限制本发明的保护范围。在阅读了本发明的记载的内容之后,技术人员可以对本发明作各种改动或修改,这些等效变化和修饰同样落入本发明权利要求所限定的范围。
Claims (7)
1.一种基于Hilbert变换的GIS局部放电信号特征提取方法,其特征在于,包括以下步骤:
采集现场检测的气体绝缘组合电器GIS局部放电信号;
对采集的GIS局部放电信号进行有效数据的提取,去除不包含放电信息的无效信号数据;
利用Hilbert变换求取放电信号的包络;
基于放电信号的包络提取局部放电脉冲波形的统计特征参数、时域特征参数和频域特征参数,利用这三类特征参数构建局部放电特征数据;
利用改进的核主成分分析法对提取的局部放电特征数据进行降维处理,核主成分分析KPCA是通过一个非线性映射将原始输入空间的特征数据投影到高维特征空间,然后在高维特征空间中利用主成分分析PCA对数据进行降维,最终得到维数较低的且能正确执行局部放电分类的特征数据。
2.根据权利要求1所述的一种基于Hilbert变换的GIS局部放电信号特征提取方法,其特征在于,所述采集现场检测的气体绝缘组合电器GIS局部放电信号是利用Keysight品牌的DSOX6002A型号示波器进行采集,其带宽为和采样率最高分别可达4GHz和20.0GSa/s,可以实现对GIS局部放电信号的有效采集。
3.根据权利要求1所述的一种基于Hilbert变换的GIS局部放电信号特征提取方法,其特征在于,所述对局部放电信号进行有效数据的提取,具体包括:只截取局部放电信号中间的一段放电脉冲数据进行特征参数的提取,将放电脉冲出现前和放电脉冲衰减完后的数据进行去除。
4.根据权利要求3所述的一种基于Hilbert变换的GIS局部放电信号特征提取方法,其特征在于,所述利用Hilbert变换求取放电信号的包络,具体包括:
首先对局部放电信号时域波形的电压值取绝对值,得到局部放电信号的单极性波形,然后利用Hilbert变换求取该单极性波形的包络,得到局部放电信号的时域波形包络;在完成时域波形包络的提取后,对放电信号波形进行快速傅里叶变换,将信号的波形由时域变换到频域中,然后利用Hilbert变换求取放电信号的频谱包络。
5.根据权利要求4所述的一种基于Hilbert变换的GIS局部放电信号特征提取方法,其特征在于,所述Hilbert变换是数学和信号处理的一种同域变换,提供了实现给定函数或傅里叶序列的调和共轭的具体方法,对于连续时域信号x(t),定义其Hilbert变换为:
其中t是时间自变量,τ是用于作卷积运算的时间变换量;
x(t)的解析信号a(t)定义为
其中j是复数的虚部单位;
解析信号的模值E(t)即为x(t)的包络
而对于长度为m的离散信号序列x(n)及其FFT序列X(k),有
其中A(k)为x(n)的离散解析信号a(n)对应的FFT序列,则x(n)的包络E(n)为
E(n)=|a(n)|=|IFFT[A(k)]| (5)
IFFT表示快速傅里叶逆变换。
6.根据权利要求5所述的一种基于Hilbert变换的GIS局部放电信号特征提取方法,其特征在于,所述基于信号的包络提取局部放电脉冲波形的统计特征参数、时域特征参数和频域特征参数,具体包括:
分别基于信号的时域波形包络和频谱包络来提取时域特征参数和频域特征参数,并联合由局部放电有效数据提取得到的统计特征参数共同构建放电特征数据;其中时域特征参数主要包括上升时间、下降时间、峰值时间和脉冲宽度,频域特征参数主要包括频率幅值均值、重心频率、均方根频率、频率方差、第一主频频率和第一主频峰值,统计特征参数主要包括均值、最大值、最小值、峰峰值、绝对均值、方根幅值、均方根值、方差、偏斜度、峭度、波形因子、峰值因子、脉冲因子和裕度因子;因此由上述三类特征参数构建的局部放电特征数据包含了24组信号特征指标,其特征维数为24维,几乎包含了局部放电信号的全部特征信息。
7.根据权利要求6所述的一种基于Hilbert变换的GIS局部放电信号特征提取方法,其特征在于,所述利用改进的KPCA算法对提取的局部放电特征数据进行降维处理,最终得到维数较低的且能正确执行局部放电分类的特征数据,具体包括:
根据公式(20)可知,Jk越小就说明类别信息量越大,那么就越容易进行分类,即这个特征向量方向上包含的类别信息量越大。但并不是类别信息量越大就越能提高分类准确率,有些特征向量方向上包含的类别信息量大,但不一定会对分类效果起到积极的作用,因此,本发明使用相似度来衡量提取的各特征向量与样本的类别向量之间的相关性,通过特征向量与样本类别向量间的相关系数来选取特征向量的组合,去除对分类效果不明显的特征向量,其公式为:
式中:ξ为特征向量与样本类别属性向量的相关系数;vi为特征向量(其中i=1,2,3,…,n);vz是原样本中的类别属性向量。若ξ的绝对值大于零或等于1,则说明两特征向量是相关的;若ξ相似趋于或等于0,则说明两向量是不相关的。通过计算特征向量与样本的类别向量间的相关系数ξ,排除不相关或相关系数非常小的特征向量,去掉后会减小后续分类器的计算量,从而减小了分类的训练和识别时间。
改进的KPCA的具体降维步骤总结如下:
假设输入空间中有M组特征数据,每组特征数据有N个特征属性(不包含类别特征),根据改进的KPCA原理,特征提取的具体实现过程为:
(1)列出一个(M×N)维数据矩阵,其中M为样本个数,N为每个样本的属性个数;
(2)通过一般的KPCA算法求解式(16)并联合式(12),求得高维特征空间中的协方差矩阵C的特征值及对应的特征向量;
(3)求取各特征向量的信息度量Jk,按照信息度量Jk从小到大的顺序,对相应的特征向量进行排序;
(4)提取前m个特征向量ν1,ν2,...,νm,从第一个特征向量ν1开始,依次根据式(21)计算ν1,ν2,...,νm与原样本中的类别属性向量vz的相关系数ξ;
(5)从步骤(4)中找到与类别向量vz的相关系数ξ为0或趋于0的特征向量,将其从特征向量组中去除,得到剩余的n个特征向量;
(6)对剩余特征向量按照相关系数从大到小的顺序重新排列;(7)通过分别计算高维特征空间中的特征数据在剩余的特征向量上的投影值,即可得到降维后维数为n的且分类效果较好的特征数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210073612.8A CN114609483A (zh) | 2022-01-21 | 2022-01-21 | 一种基于Hilbert变换的GIS局部放电信号特征提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210073612.8A CN114609483A (zh) | 2022-01-21 | 2022-01-21 | 一种基于Hilbert变换的GIS局部放电信号特征提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114609483A true CN114609483A (zh) | 2022-06-10 |
Family
ID=81857643
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210073612.8A Pending CN114609483A (zh) | 2022-01-21 | 2022-01-21 | 一种基于Hilbert变换的GIS局部放电信号特征提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114609483A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114719731A (zh) * | 2022-06-08 | 2022-07-08 | 中国航发四川燃气涡轮研究院 | 一种叶尖间隙峰峰值提取方法及叶片转速计算方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106324405A (zh) * | 2016-09-07 | 2017-01-11 | 南京工程学院 | 一种基于改进主成分分析的变压器故障诊断方法 |
CN107101813A (zh) * | 2017-04-26 | 2017-08-29 | 河北工业大学 | 一种基于振动信号的框架式断路器机械故障程度评估方法 |
CN110824310A (zh) * | 2019-11-21 | 2020-02-21 | 深圳供电局有限公司 | 局部放电类型判断方法及装置 |
-
2022
- 2022-01-21 CN CN202210073612.8A patent/CN114609483A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106324405A (zh) * | 2016-09-07 | 2017-01-11 | 南京工程学院 | 一种基于改进主成分分析的变压器故障诊断方法 |
CN107101813A (zh) * | 2017-04-26 | 2017-08-29 | 河北工业大学 | 一种基于振动信号的框架式断路器机械故障程度评估方法 |
CN110824310A (zh) * | 2019-11-21 | 2020-02-21 | 深圳供电局有限公司 | 局部放电类型判断方法及装置 |
Non-Patent Citations (2)
Title |
---|
LINGLIANG ZHOU ET: "Feature extraction and classification of partial discharge signal in GIS based on Hilbert transform", 2021 INTERNATIONAL CONFERENCE ON INFORMATION CONTROL, ELECTRICAL ENGINEERING AND RAIL TRANSIT (ICEERT), 10 January 2022 (2022-01-10), pages 208 - 213 * |
何新 等: "特征提取算法KPCA的改进与设计", 高师理科学刊, vol. 35, no. 6, 30 June 2015 (2015-06-30), pages 34 - 37 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114719731A (zh) * | 2022-06-08 | 2022-07-08 | 中国航发四川燃气涡轮研究院 | 一种叶尖间隙峰峰值提取方法及叶片转速计算方法和装置 |
CN114719731B (zh) * | 2022-06-08 | 2022-09-23 | 中国航发四川燃气涡轮研究院 | 一种叶尖间隙峰峰值提取方法及叶片转速计算方法和装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111626153B (zh) | 一种基于集成学习的局放故障状态识别方法 | |
Wang et al. | A novel weighted sparse representation classification strategy based on dictionary learning for rotating machinery | |
CN106443316B (zh) | 一种电力变压器绕组形变状态多信息检测方法及装置 | |
Lin et al. | Detection and classification of multiple power-quality disturbances with wavelet multiclass SVM | |
WO2014101579A1 (zh) | 基于奇异值分解算法的变压器局部放电模式识别方法 | |
CN112881942B (zh) | 基于小波分解与经验模态分解的异常电流诊断方法和*** | |
CN114237046B (zh) | 一种基于sift数据特征提取算法与bp神经网络模型的局部放电模式识别方法 | |
CN113155462A (zh) | 基于辛几何模式分解与图结构增强动态时间规整的轴承故障诊断方法 | |
CN114609483A (zh) | 一种基于Hilbert变换的GIS局部放电信号特征提取方法 | |
Li et al. | Pumping unit fault analysis method based on wavelet transform time-frequency diagram and cnn | |
Yang et al. | Fault Classification in Distribution Systems Using Deep Learning With Data Preprocessing Methods Based on Fast Dynamic Time Warping and Short-Time Fourier Transform | |
CN113076354A (zh) | 一种基于非侵入式负荷监测的用户用电数据分析方法和装置 | |
CN115186772B (zh) | 一种电力设备的局部放电的检测方法、装置及设备 | |
CN116432054A (zh) | 一种基于聚类分析的局部放电严重程度识别方法 | |
He et al. | Robust locally linear embedding and its application in analogue circuit fault diagnosis | |
Chen et al. | Application of EMD-AR and MTS for hydraulic pump fault diagnosis | |
Wu et al. | Classification of complex power quality disturbances based on modified empirical wavelet transform and light gradient boosting machine | |
CN114091593A (zh) | 一种基于多尺度特征融合的网络级电弧故障诊断方法 | |
Fu et al. | A Partial Discharge Signal Separation Method Applicable for Various Sensors Based on Time-Frequency Feature Extraction of t-SNE | |
Xiong et al. | Application of multi-kernel relevance vector machine and data pre-processing by complementary ensemble empirical mode decomposition and mutual dimensionless in fault diagnosis | |
Xiang et al. | Research on Transformer Voiceprint Feature Extraction Oriented to Complex Noise Environment | |
Anguswamy et al. | Distribution Network Power Quality Insights With Optimally Placed Micro-PMUs Incorporating Synthetic and Real Field Data | |
He et al. | Gas-insulated Switchgear Discharge Pattern Recognition Based on KPCA and Multilayer Neural Network | |
Sun et al. | Identification method of gas-liquid two-phase flow regime based on distance evaluation | |
Shang et al. | Transformer partial discharge fault diagnosis based on improved adaptive local iterative filtering‐bidirectional long short‐term memory |
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 |