CN109948194B - 一种高压断路器机械缺陷集成学习诊断方法 - Google Patents
一种高压断路器机械缺陷集成学习诊断方法 Download PDFInfo
- Publication number
- CN109948194B CN109948194B CN201910144353.1A CN201910144353A CN109948194B CN 109948194 B CN109948194 B CN 109948194B CN 201910144353 A CN201910144353 A CN 201910144353A CN 109948194 B CN109948194 B CN 109948194B
- Authority
- CN
- China
- Prior art keywords
- self
- encoder
- sub
- circuit breaker
- voltage circuit
- 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
- 238000003745 diagnosis Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000007547 defect Effects 0.000 title claims abstract description 30
- 239000002245 particle Substances 0.000 claims abstract description 52
- 230000006835 compression Effects 0.000 claims abstract description 51
- 238000007906 compression Methods 0.000 claims abstract description 51
- 230000009466 transformation Effects 0.000 claims abstract description 29
- 238000012549 training Methods 0.000 claims abstract description 19
- 238000005457 optimization Methods 0.000 claims abstract description 15
- 230000007246 mechanism Effects 0.000 claims abstract description 10
- 238000005070 sampling Methods 0.000 claims abstract description 10
- 238000003066 decision tree Methods 0.000 claims description 29
- 230000008569 process Effects 0.000 claims description 22
- 238000005259 measurement Methods 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000013507 mapping Methods 0.000 claims description 10
- 230000011218 segmentation Effects 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 3
- 230000000903 blocking effect Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000005056 compaction Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims description 2
- 238000009826 distribution Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 6
- 238000012512 characterization method Methods 0.000 abstract description 3
- 230000006872 improvement Effects 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 210000002569 neuron Anatomy 0.000 description 3
- 230000004913 activation Effects 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 102000002274 Matrix Metalloproteinases Human genes 0.000 description 1
- 108010000684 Matrix Metalloproteinases Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000007636 ensemble learning method Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种高压断路器机械缺陷集成学习诊断方法，属于高压断路器检测技术领域。所述方法首先获取机械振动信号数据样本及原始特征空间构建；然后基于有放回随机抽样方式组成多个特征子集；以每个特征子集和差集训练成最优旋转压缩决策子树模型；以投票机制构成最终的高压断路器机械缺陷诊断模型。本发明采用了粒子群算法对自编码器超参进行了迭代优化，避免人为调整参数的局限性，有利于高压断路器机械缺陷诊断精度的提高；本发明利用自编码器对原始特征空间进行了优化变换，使得集成学习中每个子分类器均表现出增强型决策，强化了各决策子树的表征能力，同时提高了集群诊断的辨识精度。
Description
技术领域
本发明涉及一种旋转压缩变换的集成学习方法，具体涉及一种优化特征变换空间的高压断路器机械缺陷集成学习诊断方法。
背景技术
高压断路器是区域性电网中最重要的开关设备。电力系统实际运行过程中，随着高压断路器运行时间的增加，电器与绝缘部件会有所老化，机械部件也会出现不同程度磨损，这可能会导致断路器功能失效。某种功能失效或者某种功能出现不允许的性能偏差都被称为故障，故障是设备存在的一种或者多种现象。最常见的显性故障，是通过人体的一般性感知或者正常的操作便能发现的故障，如储能不到位、拒分、拒合等一些隐性故障，是必须通过仪表或者其他设备的检测才能发现的故障，例如，分闸线圈匝间短路、分合闸速度超标以及主触头电阻变大等。对于高压断路器来讲，状态监测技术还不完善，是目前研究的热点之一。据国际大电网会议13.06工作组(参见参考文献《高压断路器机械特性监测及故障模式识别方法研究》邓帮飞，重庆大学硕士学位论文；参考文献[2]:《高压断路器事故调查》刘亚芳,国际电力1997年第3期)对包括22个国家的102个电力部门所作的第一次国际调查结果表明，在1974-1977年间，63kV及以上电压等级断路器故障中，机械性故障占70.3％。可见，对于高压断路器实施动机械特性的监测，及时了解其运行状况，掌握其运行特性变化及变化趋势，对提高其运行可靠性极为重要。
目前，能够有效反应高压断路器机械特性检测和状态诊断方法主要基于高压断路器的动触头行程信号、分合闸操作机构线圈电流信号以及机械振动信号三类。其中触头行程和线圈电流检测技术发展较早，较为成熟，但在现在应用测试过程中存在诸多不便，且部分机械缺陷难以辨识。高压断路器机械振动检测，相比于触头行程和线圈电流检测具有安装方便、检测机械缺陷种类多等优点，在近20年的发展过程中，振动信号采集系统的相关工作已比较成熟，但对振动信号特征提取与故障分类方法仍处于试验研究阶段，存在诊断准确率低、通用性差等问题，这也导致目前相关研究尚无成熟产品和应用方案，并未见在工程中成功应用的报道。因此，开展高压断路器机械状态评估与故障诊断技术的研究，研发高效、高精度的断路器机械状态检测方法具有重大而深远的意义。
发明内容
本发明的目的在于：提供一种优化特征变换空间的高压断路器机械缺陷集成学习诊断方法，该方法首先基于多类高压断路器机械缺陷的振动信号数据样本，并对其构建传统的小波包能量熵，组成振动信号原始特征空间；然后，对于振动信号原始特征空间，以有放回随机抽样方式形成多个与振动信号原始特征空间等容量的特征子集(袋内数据)以及相应的与振动信号原始特征空间的差集(袋外数据)；接着，对每个特征子集内的特征数据设计自编码器旋转压缩原始特征空间，并在变换后特征空间下训练决策树初步诊断模型，基于袋外数据评估诊断模型诊断结果，利用粒子群算法迭代优化旋转压缩变换及诊断模型参数，形成最优旋转压缩决策子树模型；最后，以投票机制融合全部最优旋转压缩决策子树模型的诊断结果，完成高压断路器机械缺陷最终诊断。
本发明采用的技术方案为：一种优化特征变换空间的高压断路器机械缺陷集成学习诊断方法，所述方法包括步骤如下：
步骤1：获取机械振动信号数据样本及原始特征空间构建；
采集高压断路器在不同机械缺陷下合闸过程的机械振动信号，基于高压断路器合闸过程表现为冲击性自由衰减信号，即信号为突变、频率随时间变化和非线性特点，因此利用小波包时频域分析，计算小波包能量熵，度量机械振动信号特征，组成机械振动信号原始特征空间。
步骤2：基于有放回随机抽样方式组成多个特征子集；
在原始特征空间下，有放回随机抽样数据样本，组成和原始特征空间中数据样本数一致的多个特征子集，特征子集中可能存在重复的特征数据，并定义每个特征子集为各自袋内数据，各特征子集与原始特征空间的差集为各自袋外数据，并且特征子集个数为最终形成最优旋转压缩决策子树的数量，即集成学习模型的容量。
步骤3：以每个特征子集和差集训练成最优旋转压缩决策子树模型；
步骤3.1基于每个特征子集建立自编码器，将原始特征空间旋转、压缩、非线性映射为新特征空间描述；
步骤3.1.1初始化自编码器的结构和相关参数；
步骤3.1.2将特征子集输入到自编码器中训练，使得自编码器的输出等于(近似等于)自编码器的输入，提取自编码器中间隐含层(自编码器一般只有1个输入层、1个隐含层和一个输出层)作为新特征集下的袋内数据，即完成原始特征空间的旋转、压缩以及非线性映射；所有新特征集构成新特征空间。
步骤3.1.3记录自编码器从输入到中间隐含层的映射关系，即自编码器旋转、压缩参数；
步骤3.2在新特征空间下以袋内数据构建决策树分类诊断模型；
步骤3.2.1将步骤3.1.2中得到的新特征空间下的袋内数据，放入决策树的树根节点；
步骤3.2.2以袋内数据中某个特征的某个数值点对袋内数据进行分割，基于基尼指数评估该特征在这一数值点的分割效果，寻找用于分割新特征空间下的袋内数据类别的最佳特征变量及数值点，将袋内数据分成两个子部分；
步骤3.2.3判断分割后的两个子部分的纯度(基尼指数)，若某个子部分的基尼指数小于要求值，则此子部分不再分割，并定义为叶节点，这个子部分样本最多的类别就是这个叶节点的输出类别；若某个子部分的基尼指数大于要求值则以步骤3.2.2的方式再次进行分割，直至满足要求，如此形成了基于新特征集下袋内数据的决策树分类诊断模型；
步骤3.3以袋外数据对生成的自编码器旋转压缩决策子树分类诊断模型进行评估，采用粒子群算法优化自编码器的超参和决策树分类诊断模型的参数，获得最优旋转压缩变换决策子树；所述超参包括自编码器的旋转和压缩参数。
步骤4：基于多个最优旋转压缩变换决策子树诊断模型，以投票机制构成最终的高压断路器机械缺陷诊断模型；
基于不同的抽样特征子集不断重复步骤3中训练-评估-优化过程，获得多个最优旋转压缩变换决策子树模型，并统计多个最优旋转压缩变换决策子树模型的分类结果，以少数服从多数的投票机制集成融合，完成高压断路器机械缺陷的最终诊断。
与现有技术相比，本发明的优点和有益效果是：
1、本发明技术方案中，采用了自编码器对原始特征空间进行了旋转压缩变换，增强相同故障的聚合性和不同故障之间的差异性，相比于其他的基于自编码器变换方法，本发明采用了粒子群算法对自编码器超参进行了迭代优化，避免人为调整参数的局限性，有利于高压断路器机械缺陷诊断精度的提高；
2、本发明技术方案中，采用样本袋装、决策树集成的诊断方法，有效降低了振动分散性对诊断结果的影响，相比于其他的基于集成学习方法，本发明利用自编码器对原始特征空间进行了优化变换，使得集成学习中每个子分类器(决策子树)均表现出增强型决策，强化了各决策子树的表征能力，同时提高了集群诊断的辨识精度；
3、本发明技术方案中，设计决策树分类诊断模型中旋转压缩映射变换的评估方法，为粒子群优化过程提供了评价标准，相比于其他的粒子群优化和单一模型诊断方法，本发明灵活运用了袋外数据不参与决策树分类诊断模型训练的优势，提出袋内训练-袋外评估的迭代优化方案，有效地度量了自编码器对原始特征空间的旋转、压缩映射变换程度的优劣，形成了最优旋转压缩变换决策子树，为进一步提高机器学习诊断能力提供发展方向。
附图说明
下面结合附图对本发明进一步说明。
图1是本发明提供的高压断路器机械缺陷集成学习诊断方法流程图；
图2是旋转压缩变换最优特征空间的集成学习方法结构图；
图3是最优旋转压缩优化特征空间的决策子树结构图；
图4是最优旋转压缩优化特征空间的决策子树生成流程图；
图5是旋转压缩映射的自编码器结构图；
图6是决策树结构图；
图7是粒子群算法流程图。
具体实施方式
下面将结合实施例和附图，对本发明实施例中的技术方案进行清楚、完整地描述。
本发明提供了一种优化特征变换空间的高压断路器机械缺陷集成学习诊断方法，图1展示出了本实施例中的所述诊断方法的流程图，图2为本实施例中所述方法的结构框图；所述诊断方法具体包括下述步骤：
步骤1：获取振动信号数据样本及原始特征空间构建；
步骤1.1使用振动信息测量设备采集高压断路器在不同机械缺陷下合闸过程的机械振动信号，得到数据样本数量为m个；
步骤1.2基于高压断路器合闸过程表现为冲击性自由衰减信号，即振动信号为突变、频率随时间变化和非线性特点，因此利用小波包时频域分析，计算小波包能量熵，度量机械振动信号特征，具体步骤如下：
步骤1.2.1对高压断路器机械振动信号Xq进行小波包时-频变换，其中q＝1,2,…,m表示测量数据样本(简称测量样本)编号，获得时-频能量矩阵Aq,t×f，并在时间轴t、频率轴f方向分别等分n1和n2份，每个子时间段元素个数为t/n1，每个子频率段元素个数为f/n2，计算每个时-频块的能量，如公式(1)所示。
其中，Eq,i,j表示第q个测量数据样本小波包时-频能量矩阵Aq,t×f进行n1×n2分块后第(i，j)个时-频块的能量，
表示第q个测量样本中第(i，j)个能量块中第(S1，S2)个元素，n1和n2表示为沿时间方向(行方向)和频率方向(列方向)对矩阵Aq,t×f的等分个数，t和f表示矩阵Aq,t×f的行数和列数。
步骤1.2.2根据公式(1)得到第q个测量样本的小波包时-频块的能量Eq,i,j，可以计算小波包能量熵，如公式(2)和(3)所示，由此，构成了高压断路器第q个测量样本的振动信号原始特征空间特征描述WPEq＝{WPEq,i,f,WPEq,t,j},i＝1,2,...,n1；j＝1,2,...,n2；q＝1,2,...,m。
步骤2：基于有放回随机抽样方式组成多组训练子集；
将步骤1中得到的原始特征空间中特征描述WPEq＝{WPEq,i,f,WPEq,t,j}定义为Cq＝{xq,1,xq,2,...,xq,w}，其中，i＝1,2,...,n1；j＝1,2,...,n2；q＝1,2,...,m，Cq表示第q个测量样本的小波包能量熵的特征集，w表示第q个测量样本的特征个数，即特征空间维度，定义每个特征集Cq的故障类型编号为Lg，g＝1,2,…,G，表示有G个机械故障的类别个数。
步骤2.1在原始特征空间下有放回随机抽样测量样本，组成和原始特征空间中样本数一致的多个特征子集
特征子集中可能存在重复的特征数据，定义每个特征子集为各自袋内数据，其中N表示特征子集个数，且为最终形成最优旋转压缩决策子树的数量，即集成学习模型的容量。
步骤2.2将各特征子集InCnum与原始特征空间中测量样本作差集，定义其为各自的袋外数据ExCnum，num＝1,2,…,N。
步骤3：以每个特征子集InCnum(袋内数据)和差集ExCnum(袋外数据)，训练评估出最优旋转压缩变换决策子树模型，模型结构如图3所示，生成流程如图4所示，具体步骤如下：
步骤3.1基于每个特征子集InCnum建立自编码器，将原始特征空间旋转、压缩、非线性映射为新特征空间描述；
步骤3.1.1初始化自编码器的旋转压缩参数，即超参，如图5，如隐含层神经元(y1,y2,…,yp)个数p、激活函数fθ(·)类型、正则化系数λ等；θ的取值是1～p，p为自编码器中编码器的输出个数。
步骤3.1.2将袋内数据特征子集
输入到自编码器中进行训练，使得自编码器中解码器输出向量
等于(近似等于)编码器输入向量
定义自编码器中隐含层(编码器)输出为
输入层-隐含层-输出层满足关系如公式(4)和(5)所示：
y＝fθ(Wx+b) (4)
公式(4)表示自编码器的编码过程，即输入层到隐含层输出过程，公式(5)表示自编码器的解码过程，即隐含层输出到输出层输出过程，其中编码解码激活函数分别为fθ(α)＝1/(1+e-α)，
当自编码器的超参数集
确定后即可得到自编码器输出向量
定义指标J，计算从输入到输出的误差，如公式(6)所示：
其中，
λ为正则项系数，自编码器的超参，初始化权值参数和偏执参数
根据计算公式(6)计算误差，基于神经网络误差反向原理，利用随机梯度下降方法更新权值参数和偏执参数，使得输入-输出的误差最小。显然，可以看出从公式(4)可以看出，当参数{W,b}确定后，从输入层到隐含层输出可以理解为
则可以理解为对原始特征空间进行了一定的旋转(坐标变换)、压缩(矩阵非满秩)以及非线性变换(激励函数为非线性)，形成新特征空间。
步骤3.1.3记录自编码器从输入层到中间隐含层的映射关系，即自编码器超参{W,b}；
步骤3.2在新特征空间下以袋内数据构建决策树分类诊断模型，决策树分类诊断模型结构如图6所示，其中集合S表示新特征空间下的袋内数据，Si，j，k表示第i-j-k次分割后的子集。
步骤3.2.1将步骤3.1.2中得到的新特征空间下的袋内数据集合S，放入决策树的树根节点；
步骤3.2.2以树根节点为例说明，袋内数据的样本个数是m，每个样本是由一个w维向量组成，可以理解为一个m×w的矩阵，矩阵行数表示样本个数，列数表示w个特征数，并且每一行(样本)对应一个机械缺陷的类别。选择某个一列(特征yi)，至少有N-1种数值Δi可以将样本分为2部分，即小于等于Δi的一部分S1，大于Δi的一部分S2，然后计算每个部分S1和S2中机械缺陷类别的出现概率，利用公式(7)计算出特征yi在Δi这个数值点上的Gini指数，最后根据公式(8)计算出两个部分S1和S2的加权Gini指数Gini(S,yi)，表示出特征yi在Δi这个数值点对袋内数据划分纯度的度量，其中|S|表示集合S中样本的数量。不断更换特征yi上的Δi值(最多N-1次)计算，以Gini(S,yi)最小判断出特征yi上的最优分割值。
再更换特征yi(w次)，根据上述过程计算出每个特征下的最优分割值，并比较每个特征的最优分割值下的Gini(S,yi)，选择最小值对应的特征和特征分割点为最佳特征变量及数值，将袋内数据集S分成两个子部分S1和S2。公式(7)、(8)如下：。
其中x表示一个变量，其存在K类状态值，Pk表示变量x属于第k类的概率值，公式(7)可以表示随机变量的纯度，当变量x仅存在一类状态值，则基尼指数为0，表示纯度最大；当变量x存在类别数K趋近于无穷大，且各类别出现概率一致，则基尼指数为1，表示纯度最小；利用公式(7)可以看出公式(8)表征了由特征变量yi的某个数值Δi对新袋内数据特征集进行分割得到两个子部分S1和S2的纯度度量，递归寻找用于分割新特征集下的袋内数据类别的最佳特征变量及数值，将新特征集下的袋内数据集S分成两个子部分S1和S2；
步骤3.2.3判断分割后的两个子部分的纯度(基尼指数)，若某个子部分的基尼指数小于要求值γ，则此子部分不再分割，并定义为叶节点，这个子部分样本最多的类别就是这个叶节点的输出类别；若某个子部分的基尼指数大于要求值γ，则以步骤3.2.2的方式继续进行分割，直至满足要求，形成了基于新特征集下袋内数据的决策树分类诊断模型；
步骤3.3以袋外数据对生成的自编码器和决策树分类诊断模型进行评估，采用粒子群算法优化自编码器的超参和决策树诊断模型的参数，流程如图7所示，获得最优旋转压缩变换决策子树，具体如下：
步骤3.3.1设置粒子群算法参数，如粒子数Np、最大迭代次数intermax、惯性因子ω及学习因子c1和c2等，定义每个粒子由自编码器超参构成，如隐含层神经元个数p、正则化系数λ以及决策树停止分割的纯度要求值γ等参数，由向量Pd＝[p,λ,γ,...]表示，初始化所有粒子向量以及每个粒子对应的自编码器参数
步骤3.3.2基于步骤2抽样的特征子集InCnum，根据步骤3.1训练所有粒子代表的自编码器，生成各粒子的旋转压缩特征，并利用步骤3.2生成各粒子旋转压缩特征下的决策树分类诊断模型；
步骤3.3.3将步骤2中特征子集InCnum对应的袋外数据ExCnum，放入所有粒子的自编码器和决策树分类诊断模型，获得并记录各粒子的诊断结果；
步骤3.3.4对袋外数据ExCnum的诊断结果进行排序，记录并更新所有粒子的历史全局最优位置Pgbest，以及每个粒子的历史自身最优位置Pi，best，其中i表示粒子编号。判断当前粒子是否优于历史自身最优值或是历史全局最优值：若是，则更新每个粒子的自身最优值及历史全局最优值，并判断是否到达迭代截止条件；若不是，则不改变每个粒子的历史自身最优及历史全局最优值，并判断是否到达迭代截止条件。若没有到达迭代截止条件，则以公式(9)对每个粒子进行更新，并以更新后的所有粒子参数，重新返回步骤3.3.2重新进行自编码器训练和决策树分类诊断模型构建，直至达到迭代截止条件停止优化过程。
其中，vi(inter)表示第i个粒子在第inter次迭代时的更新速度，ω表示粒子的惯性常数，c1和c2分别为历史自身最优学习因子和历史全局最优学习因子，r1和r2分别是两个服从正态分布的随机数；
步骤3.3.5根据历史全局最优的粒子参数，重新进行步骤3.3.2训练过程，获得最优旋转压缩变换决策子树；
步骤4：基于多个最优旋转压缩变换决策子树，以投票机制构成最终的高压断路器机械缺陷诊断模型；
基于不同的抽样子集不断重复步骤3中训练-评估-优化过程，获得多个最优旋转压缩变换子树ORCTTn(Optimal Rotation&Compression Transform Tree)，n＝1,2,…,N，并统计多个最优旋转压缩变换决策子树的分类结果，以少数服从多数的投票机制集成融合，如公式(10)所示，完成高压断路器机械缺陷的最终诊断。
其中，C表示测试样本小波包能量熵特征向量，Lg表示第n个最优旋转压缩变换决策子树的诊断结果(即缺陷类别)，I(ORCTTn,C,Lg)表示第g行为1的列向量，
表示获得列向量最大值索引的函数。
本发明应用实例如下，以某型号高压断路器为实验平台设置正常情况、分闸弹簧疲劳、合闸弹簧疲劳、油缓冲器漏油、传动轴阻尼增加及地脚螺栓松动6种运行工况，一共做了合闸试验各300次，每种运行工况50次，随机选用每种运行工况的70％数据，即35个数据进行训练，剩余数据进行测试验证。训练样本数m＝210，在小波包时频矩阵中，沿时间轴等分n1＝12，沿频率轴等分n2＝13，则共分为12×13个能量块，小波包能量熵构成的原始特征空间维度w为12+13等于25，设置最优旋转压缩子决策树数量N＝100，粒子群算法粒子数Np＝30，最大迭代次数intermax＝100，惯性因子ω＝0.5，学习因子c1和c2均等于2，每个粒子由隐含层神经元个数p、正则化系数λ以及决策树停止分割的纯度要求值γ三个参数构成。利用本专利所述方法，设计优化旋转压缩特征变换的决策树集成诊断方法，通过测试数据的诊断结果表明，该方法可以降低了振动分散性对诊断结果的影响，强化了高压断路器同类典型工况样本在特征空间的聚合程度以及不同故障间的分离程度，同时解决了诊断模型的超参选择问题，增强了各决策子树的表征能力，有效提升了故障分类的准确率。
最后应当说明的是：所描述的实施例仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
Claims (4)
1.一种高压断路器机械缺陷集成学习诊断方法，其特征在于：所述方法包括步骤如下，
步骤1：获取机械振动信号数据样本及原始特征空间构建；
采集高压断路器在不同机械缺陷下合闸过程的机械振动信号，利用小波包时频域分析，计算小波包能量熵，度量机械振动信号特征，组成机械振动信号原始特征空间；
步骤2：基于有放回随机抽样方式组成多个特征子集；
在原始特征空间下，有放回随机抽样数据样本，组成和原始特征空间中数据样本数一致的多个特征子集，特征子集中可能存在重复的特征数据，并定义每个特征子集为各自袋内数据，各特征子集与原始特征空间的差集为各自袋外数据，并且特征子集个数为最终形成最优旋转压缩决策子树的数量；
步骤3：以每个特征子集和差集训练成最优旋转压缩决策子树模型；
步骤3.1基于每个特征子集建立自编码器，将原始特征空间旋转、压缩、非线性映射为新特征空间描述；
步骤3.1.1初始化自编码器的结构和相关参数；
步骤3.1.2将特征子集输入到自编码器中训练，使得自编码器的编码器输出等于解码器的输入，提取自编码器中间隐含层作为新特征集下的袋内数据，即完成原始特征空间的旋转、压缩以及非线性映射；所有新特征集构成新特征空间；
步骤3.1.3记录自编码器从输入到中间隐含层的映射关系，即自编码器旋转、压缩参数；
步骤3.2在新特征空间下以袋内数据构建决策树分类诊断模型；
步骤3.2.1将步骤3.1.2中得到的新特征空间下的袋内数据，放入决策树的树根节点；
步骤3.2.2以袋内数据中某个特征的某个数值点对袋内数据进行分割，基于基尼指数评估该特征在这一数值点的分割效果，寻找用于分割新特征空间下的袋内数据类别的最佳特征及数值点，将袋内数据分成两个子部分；
步骤3.2.3判断分割后的两个子部分的基尼指数，若某个子部分的基尼指数小于要求值，则此子部分不再分割，并定义为叶节点，这个子部分样本最多的类别就是这个叶节点的输出类别；若某个子部分的基尼指数大于要求值则以步骤3.2.2的方式再次进行分割，直至满足要求，如此形成了基于新特征空间下袋内数据的决策树分类诊断模型；
步骤3.3以袋外数据对生成的自编码器旋转压缩决策子树分类诊断模型进行评估，采用粒子群算法优化自编码器的超参和决策树分类诊断模型，获得最优旋转压缩变换决策子树；所述超参包括自编码器的旋转和压缩参数；
步骤4：基于多个最优旋转压缩变换决策子树诊断模型，以投票机制构成最终的高压断路器机械缺陷诊断模型；
基于不同的抽样特征子集不断重复步骤3中训练-评估-优化过程，获得多个最优旋转压缩变换决策子树模型，并统计多个最优旋转压缩变换决策子树模型的分类结果，以少数服从多数的投票机制集成融合，完成高压断路器机械缺陷的最终诊断。
2.根据权利要求1所述的一种高压断路器机械缺陷集成学习诊断方法，其特征在于：步骤1具体如下，
步骤1.1采集高压断路器在不同机械缺陷下合闸过程的机械振动信号，得到数据样本数量为m个；
步骤1.2利用小波包时频域分析，计算小波包能量熵，度量机械振动信号特征，具体步骤如下：
步骤1.2.1对高压断路器机械振动信号Xq进行小波包时-频变换，其中q＝1,2,…,m表示测量样本编号，获得时-频能量矩阵Aq,t×f，并在时间轴t、频率轴f方向分别等分n1和n2份，每个子时间段元素个数为t/n1，每个子频率段元素个数为f/n2，计算每个时-频块的能量，如公式(1)所示：
其中，Eq,i,j表示第q个测量数据样本小波包时-频能量矩阵Aq,t×f进行n1×n2分块后第(i，j)个时-频块的能量，
表示第q个测量样本中第(i,j)个能量块中第(S1,S2)个元素，n1和n2表示为沿时间方向和频率方向对矩阵Aq,t×f的等分个数，t和f表示矩阵Aq,t×f的行数和列数；
步骤1.2.2根据公式(1)得到第q个测量样本的小波包时-频块的能量Eq,i,j，计算小波包能量熵，如公式(2)和(3)所示，由此，构成了高压断路器第q个测量样本的振动信号原始特征空间特征描述WPEq＝{WPEq,i,f,WPEq,t,j},i＝1,2,...,n1；j＝1,2,...,n2；q＝1,2,...,m：
3.根据权利要求1所述的一种高压断路器机械缺陷集成学习诊断方法，其特征在于：步骤3.3所述的优化和评估，具体如下，
步骤3.3.1设置粒子群算法参数，定义每个粒子由自编码器超参构成，由向量Pd＝[p,λ,γ,...]表示，初始化所有粒子向量以及每个粒子对应的自编码器参数
步骤3.3.2基于特征子集，训练所有粒子代表的自编码器，生成各粒子的旋转压缩参数，并构建各粒子旋转压缩参数下的决策树分类诊断模型；
步骤3.3.3将特征子集对应的袋外数据，放入所有粒子的自编码器和决策树分类诊断模型，获得并记录各粒子的诊断结果；
步骤3.3.4对袋外数据的诊断结果进行排序，记录并更新所有粒子的历史全局最优位置Pgbest，以及每个粒子的历史自身最优位置Pi，best，其中i表示粒子编号；判断当前粒子是否优于历史自身最优值或是历史全局最优值：若是，则更新每个粒子的历史自身最优值及历史全局最优值，并判断是否到达迭代截止条件；若不是，则不改变每个粒子的历史自身最优及历史全局最优值，并判断是否到达迭代截止条件；若没有到达迭代截止条件，则以公式(9)对每个粒子进行更新，并以更新后的所有粒子参数，重新返回步骤3.3.2重新进行自编码器训练和决策树分类诊断模型构建，直至达到迭代截止条件停止优化过程；
其中，vi(inter)表示第i个粒子在第inter次迭代时的更新速度，ω表示粒子的惯性常数，c1和c2分别为历史自身最优学习因子和历史全局最优学习因子，r1和r2分别是两个服从正态分布的随机数；
步骤3.3.5根据历史全局最优的粒子参数，重新进行步骤3.3.2训练过程，获得最优旋转压缩变换决策子树。
4.根据权利要求1所述的一种高压断路器机械缺陷集成学习诊断方法，其特征在于：步骤4所述的投票机制，如公式(10)所示，
其中，C表示测试样本小波包能量熵特征向量，Lg表示第n个最优旋转压缩变换决策子树的诊断结果，I(ORCTTn,C,Lg)表示第g行为1的列向量，
表示获得列向量最大值索引的函数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910144353.1A CN109948194B (zh) | 2019-02-27 | 2019-02-27 | 一种高压断路器机械缺陷集成学习诊断方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910144353.1A CN109948194B (zh) | 2019-02-27 | 2019-02-27 | 一种高压断路器机械缺陷集成学习诊断方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109948194A CN109948194A (zh) | 2019-06-28 |
| CN109948194B true CN109948194B (zh) | 2020-07-03 |
Family
ID=67007066
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910144353.1A Active CN109948194B (zh) | 2019-02-27 | 2019-02-27 | 一种高压断路器机械缺陷集成学习诊断方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109948194B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110595751B (zh) * | 2019-09-19 | 2020-12-18 | 华东理工大学 | 基尼指数引导的早期故障特征小波重构方法及其应用 |
| CN111026090A (zh) * | 2019-12-26 | 2020-04-17 | 浙江力石科技股份有限公司 | 物联网设备故障识别方法、系统、装置及可存储介质 |
| CN112485394A (zh) * | 2020-11-10 | 2021-03-12 | 浙江大学 | 一种基于稀疏自编码和极限学习机的水质软测量方法 |
| CN113687222B (zh) * | 2021-08-24 | 2024-07-02 | 青岛理工大学 | 一种sf6断路器弧触头状态评估方法及系统 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107633114A (zh) * | 2017-08-21 | 2018-01-26 | 江苏省电力试验研究院有限公司 | 一种高压断路器潜伏性机械故障振动信号仿真方法及装置 |
| CN108181107A (zh) * | 2018-01-12 | 2018-06-19 | 东北电力大学 | 计及多分类目标的风电机组轴承机械故障诊断方法 |
| CN108680348A (zh) * | 2018-05-14 | 2018-10-19 | 国网山东省电力公司莱芜供电公司 | 一种基于随机森林的断路器机械故障诊断方法及系统 |
| CN108734202A (zh) * | 2018-04-27 | 2018-11-02 | 西安工程大学 | 一种基于改进bp神经网络的高压断路器故障诊断方法 |
| CN108760031A (zh) * | 2018-04-20 | 2018-11-06 | 中国电力科学研究院有限公司 | 一种用于气体绝缘金属封闭开关设备振动模拟实验的系统及方法 |
| CN109102032A (zh) * | 2018-09-03 | 2018-12-28 | 中国水利水电科学研究院 | 一种基于深度森林和自编码的泵站机组诊断方法 |
| CN109214460A (zh) * | 2018-09-21 | 2019-01-15 | 西华大学 | 基于相对变换与核熵成分分析的电力变压器故障诊断方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103487749B (zh) * | 2013-09-18 | 2016-04-13 | 国家电网公司 | 高压断路器机械状态在线监测诊断系统及方法 |
| US11507064B2 (en) * | 2016-05-09 | 2022-11-22 | Strong Force Iot Portfolio 2016, Llc | Methods and systems for industrial internet of things data collection in downstream oil and gas environment |
| CN106372724A (zh) * | 2016-08-31 | 2017-02-01 | 西安西拓电气股份有限公司 | 一种人工神经网络算法 |
-
2019
- 2019-02-27 CN CN201910144353.1A patent/CN109948194B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107633114A (zh) * | 2017-08-21 | 2018-01-26 | 江苏省电力试验研究院有限公司 | 一种高压断路器潜伏性机械故障振动信号仿真方法及装置 |
| CN108181107A (zh) * | 2018-01-12 | 2018-06-19 | 东北电力大学 | 计及多分类目标的风电机组轴承机械故障诊断方法 |
| CN108760031A (zh) * | 2018-04-20 | 2018-11-06 | 中国电力科学研究院有限公司 | 一种用于气体绝缘金属封闭开关设备振动模拟实验的系统及方法 |
| CN108734202A (zh) * | 2018-04-27 | 2018-11-02 | 西安工程大学 | 一种基于改进bp神经网络的高压断路器故障诊断方法 |
| CN108680348A (zh) * | 2018-05-14 | 2018-10-19 | 国网山东省电力公司莱芜供电公司 | 一种基于随机森林的断路器机械故障诊断方法及系统 |
| CN109102032A (zh) * | 2018-09-03 | 2018-12-28 | 中国水利水电科学研究院 | 一种基于深度森林和自编码的泵站机组诊断方法 |
| CN109214460A (zh) * | 2018-09-21 | 2019-01-15 | 西华大学 | 基于相对变换与核熵成分分析的电力变压器故障诊断方法 |
Non-Patent Citations (2)
| Title |
|---|
| "Lung Nodule Classification Using Deep Features in CT Images";Devinder Kumar.etc.;《2015 12th Conference on Computer and Robot Vision》;20151231;全文 * |
| "基于神经网络集成模型在高压断路器机械故障诊断中的应用";赵科等;《高压电器》;20180716;第54卷(第7期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109948194A (zh) | 2019-06-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109948194B (zh) | 一种高压断路器机械缺陷集成学习诊断方法 | |
| CN110006645B (zh) | 一种多源融合的高压断路器机械故障诊断方法 | |
| CN110376522B (zh) | 一种数据融合的深度学习网络的电机故障诊断方法 | |
| CN111914883B (zh) | 一种基于深度融合网络的主轴轴承状态评估方法及装置 | |
| CN110132554B (zh) | 一种深度拉普拉斯自编码的旋转机械故障诊断方法 | |
| CN109344517A (zh) | 一种新能源汽车的高压绝缘故障诊断方法 | |
| CN116226646B (zh) | 轴承健康状态及剩余寿命的预测方法、系统、设备及介质 | |
| CN114623050A (zh) | 风电机组的运行数据异常检测方法、系统及介质 | |
| CN113538353B (zh) | 基于单通道图数据增强和迁移训练残差网络的五相异步电动机滚动轴承故障诊断方法 | |
| CN117056814B (zh) | 一种变压器声纹振动故障诊断方法 | |
| CN113987294A (zh) | 一种基于遗传优化gru神经网络的cvt在线故障诊断方法 | |
| CN113985156A (zh) | 一种基于变压器声纹大数据的智能故障识别方法 | |
| CN113533995A (zh) | 一种动力电池一致性检测方法 | |
| CN114236332A (zh) | 一种电力电缆绝缘状态判断方法及系统 | |
| CN115586406A (zh) | 一种基于特高频信号的gis局部放电故障诊断方法及系统 | |
| CN114705432A (zh) | 防爆电机轴承健康状态评估方法及系统 | |
| CN114487643A (zh) | 一种特高压gil设备现场交接验收综合测试平台 | |
| CN113611568A (zh) | 一种基于遗传卷积深度网络的真空断路器 | |
| CN112329825B (zh) | 基于信息分维及提升决策树的变压器机械故障诊断方法 | |
| Song et al. | Research on rolling bearing fault diagnosis method based on improved LMD and CMWPE | |
| CN115564075B (zh) | 一种城市电网主配一体故障协同诊断方法、系统 | |
| CN112069621A (zh) | 基于线性可靠度指标的滚动轴承剩余使用寿命的预测方法 | |
| CN113537010B (zh) | 基于单通道图数据增强和迁移训练残差网络的十五相异步电动机滚动轴承故障诊断方法 | |
| CN116008733A (zh) | 一种基于集成深度神经网络的单相接地故障诊断方法 | |
| CN114137915A (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 |