CN112327219A - 特征自动挖掘和参数自动寻优的配电变压器故障诊断方法 - Google Patents
特征自动挖掘和参数自动寻优的配电变压器故障诊断方法 Download PDFInfo
- Publication number
- CN112327219A CN112327219A CN202011178015.9A CN202011178015A CN112327219A CN 112327219 A CN112327219 A CN 112327219A CN 202011178015 A CN202011178015 A CN 202011178015A CN 112327219 A CN112327219 A CN 112327219A
- Authority
- CN
- China
- Prior art keywords
- distribution transformer
- encoder
- function
- fault diagnosis
- automatic
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005457 optimization Methods 0.000 title claims abstract description 46
- 238000003745 diagnosis Methods 0.000 title claims abstract description 42
- 238000005065 mining Methods 0.000 title claims abstract description 12
- 238000012549 training Methods 0.000 claims abstract description 34
- 238000007637 random forest analysis Methods 0.000 claims abstract description 31
- 239000013598 vector Substances 0.000 claims abstract description 24
- 238000000605 extraction Methods 0.000 claims abstract description 17
- 238000012360 testing method Methods 0.000 claims abstract description 14
- 238000002372 labelling Methods 0.000 claims abstract description 5
- 230000006870 function Effects 0.000 claims description 54
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 18
- 238000013528 artificial neural network Methods 0.000 claims description 8
- 230000004913 activation Effects 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000003066 decision tree Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 4
- 239000003795 chemical substances by application Substances 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 238000011478 gradient descent method Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 3
- 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 claims description 2
- 238000001514 detection method Methods 0.000 description 7
- 238000004587 chromatography analysis Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000004804 winding Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013145 classification model Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000013021 overheating Methods 0.000 description 1
- 230000003595 spectral effect Effects 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/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
- G01R31/62—Testing of transformers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- 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/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- 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
-
- 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
Abstract
本发明涉及一种特征自动挖掘和参数自动寻优的配电变压器故障诊断方法。包括：将振动信号采集装置安装在配电变压器上，并从振动信号采集装置中采集配电变压器运行时的振动波形；构建基于二次调优后的堆栈自编码器的配电变压器故障特征提取模型；利用二次调优后的堆栈自编码器提取振动信号特征向量yn，并对相应的特征打上标签，建立一个包含正常及各类故障的数据库；分割数据集，即将数据集按X1：X2的比例拆分为训练集和测试集；利用训练集训练随机森林分类器；基于训练的随机森林分类器的网络参数，建立配电变压器故障诊断模型，实现配电变压器故障诊断。本发明能够显著提高配电变压器故障诊断的准确性，并具有良好的鲁棒性和突出的诊断性能。
Description
技术领域
本发明涉及一种特征自动挖掘和参数自动寻优的配电变压器故障诊断方法。
背景技术
配电系统中的变压器是我国电力系统变电传输环节中一个极为重要的输变电设备，而用户使用的配电变压器就是面向用户的一种终端输变电设备，数量多应用范围广，其可靠的运行对整个电网的稳定和安全运行至关重要，一旦变压器发生了故障就有可能给用户造成巨大的经济损失。而根据相关的统计显示,变压器绕组的故障在我国变压器总事故中约占一半以上。因此极为有必要对配电系统变压器绕组的正常工作运行状态和可靠性进行实时有效地的监测，及时发现可能造成故障的原因和隐患。
目前，电力变压器的运行状态检测的方法可以大致分为两个主要的大类：分别是离线运行状态检测和在线运行状态监测。离线监测的方法主要有频率响应法、短路阻抗法和低压脉冲法等；在线监测的方法很多，但目前研究最多的是油色谱分析法和振动分析法，且对油色谱分析法的研究更加成熟，应用也更加广泛。但油色谱分析法能检测的故障类型主要是过热和放电等引起的绝缘油性质发生变化的故障，对变压器的机械故障反应比较迟钝甚至无法检测，且用于油色谱分析的设备大多价格昂贵、安装不便。因此，相对于油色谱分析法，基于振动分析法的变压器在线检测方法不仅对变压器主要部件的机械结构变化或故障反应灵敏(主要是变压器铁芯、绕组以及分接开关等部件发生位移、松动、磨损、变形等)，而且基于振动分析法的变压器检测系统与整个电力系统之间没有电气连接，可以再电力系统正常运行时进行监测，从而可以快速高效、安全可靠地达到变压器在线状态检测、故障诊断以及潜伏性故障诊断的目的。振动信号主要有时域特征量和频域特征量这两类特征量。时域上的特征量包括波形平均值、方差、幅值、波形相似度、动态时间距离等：频域上的特征量包括频谱包络、时频能量熵、能-频分布等。这些特征量在一定程度上能够刻画振动信号的形态特征，但实际应用的效果不理想，易受扰动和噪声等因素的影响。
多年来经过研究人员的努力，已经在振动信号获取、特征信号提取和故障检测模型建立等方面取得了一些进展，但是基于振动分析法的变压器在线检测方法尚不成熟，仍然需要做更深入的研究。
为有效解决配电变压器故障诊断中面临的数据特征提取、局部最优、实时故障分类等问题，本发明提出了一种基于自动调参的堆栈自编码器(auto-encoders,SAE)和随机深林(Random Forest,RF)分类器的配电变压器故障诊断方法。
发明内容
本发明的目的在于提供一种特征自动挖掘和参数自动寻优的配电变压器故障诊断方法，能够显著提高配电变压器故障诊断的准确性，并具有良好的鲁棒性和突出的诊断性能。
为实现上述目的，本发明的技术方案是：一种特征自动挖掘和参数自动寻优的配电变压器故障诊断方法，包括如下步骤：
步骤S1、将振动信号采集装置安装在配电变压器上，并从振动信号采集装置中采集配电变压器运行时的振动波形；
步骤S2、构建基于二次调优后的堆栈自编码器的配电变压器故障特征提取模型；
步骤S3、利用二次调优后的堆栈自编码器提取振动信号特征向量yn，并对相应的特征打上标签，建立一个包含正常及各类故障的数据库；
步骤S4、分割数据集，即将数据集按X1：X2的比例拆分为训练集和测试集；
步骤S5、利用训练集训练随机森林分类器；
步骤S6、基于步骤S5训练的随机森林分类器的网络参数，建立配电变压器故障诊断模型，实现配电变压器故障诊断。
在本发明一实施例中，所述步骤S2具体实现如下：
S21、创建一个堆栈自编码器模型，利用测量得到的振动数据逐层学习数据，提取数据特征，得到每层自编码器的偏差和权重；再通过前馈反向神经网络来缩小网络损失值，从而不断更新权重和偏差，得到堆栈自编码器中每个编码器的权重和偏差，具体实现方式如下：
假设一训练样本x＝[x1,x2,...,xn]T∈Rn，由公式(1)可将其向隐含层进行映射，映射后的隐含层表示为y1＝[y1,y2,...,ym]T∈Rm；
y1＝fθ(x)＝se(Wx+b) (1)
其中，se(·)为编码阶段的激活函数；θ＝{W,b}为编码阶段的参数，其中包含权值矩阵W∈Rm×n、偏置矩阵b∈Rm；
接着，采用解码器公式(2)将隐含表示变换回x的重构向量
其中，sd(·)为解码阶段的激活函数；
为解码阶段的参数，包含权值矩阵
偏置矩阵
近似于x的程度越高，自编码器网络的性能越优良；因此，要不断最小化
和x之间的误差，即最小化损失函数：
对于网络参数
的调整，同BP神经网络的调节机制，采用梯度下降法来最小化平均重构误差；具体的损失函数为：
其中，
为交叉熵，用于度量输入向量与重构向量之间的差异，两者差异越小，交叉熵的值越小，即损失函数的值越小；
随后将此编码器提取出的隐含层表示特征y1作为下一个编码器的输入，同样利用公式(1)可得y1映射后的隐含层表示为y2；再采用公式(2)可得y2的重构向量
此时同样利用公式(4)来不断最小化y1和
之间的误差；利用此步骤依次对n个编码器进行训练，可得每个自编码器的权重和偏差；
S22、设置堆栈自编码器特征提取网络中的迭代次数和学习率为可调超参数，并设定可调超参数调整范围，将这些可调超参数传递给步骤S23；
S23、利用贝叶斯优化进行超参数自动寻优，贝叶斯优化公式表示为：
其中，D＝{(x1,y1),(x2,y2),...,(xn,yn)}表示已采集样本点，p(f)为f的先验概率，p(f|D)为f的后验概率，f表示目标函数或模型中的参数，后验概率分布是通过已采集样本点对先验概率进行修正后目标函数或模型中参数的置信度；采用高斯过程作为概率代理模型，该过程是多元高斯概率分布的泛化，由协方差函数和均值函数构成，表达式为：
y＝gp(m(x),k(x,x')) (6)
其中，m(x)为均值函数，k(x,x')为协方差函数，该过程对一系列离散数据点(xi,yi)进行函数拟合，得到一个逼近真实值的目标函数；
采用PI函数作为采集函数，其形式为：
其中f(X)为目标函数，f(X+)是目前为止最优的目标函数，μ(x)是f(X)后验分布的均值，δ(x)是协差矩阵，ν是权衡系数，用于避免陷入局部最优；采集函数求解的过程就是找到一个新的x，使PI(x)最大；
S24、判断超参数寻优的目标函数是否达到预设的期望值，若是则转至S25，否则转至S23；
S25、得到基于堆栈自编码器的配电变压器故障初始特征提取模型；
S26、对堆栈自编码器进行二次调优，得到基于二次调优后的堆栈自编码器的配电变压器故障特征提取模型。
在本发明一实施例中，所述步骤S26具体实现如下：
S261、根据式(8)和式(9)计算各层的激励响应；
Z(l)＝W(l-1)a(l-1)+b(l-1) (8)
a(l)＝f(Z(l)) (9)
S262、根据式(10)计算残差；
S263、根据式(11)、式(12)计算权重和偏差的更新公式；
S264、根据式(13)、(14)计算权重和偏差的改变值；
其中ΔW(l)和W(l)是大小相同的矩阵，Δb(l)和b(l)是行数相同的向量，初始值都设置为0；
S265、根据式(15)、(16)得到更新后的权重和偏差；
其中，m表示样本数量，μ为学习速率。
在本发明一实施例中，所述步骤S5具体实现如下：
S51、利用Bootstrap抽样从原始数据中有放回的随机抽取X个样本构成N个子训练集，每个样本的样本容量与原始数据集一样；
S52、根据N个子训练集训练生成N个决策树模型，从而得到N个分类结果；
S53、通过投票表决得出最终的分类结果，此处投票表决服从“少数服从多数”的原则，即最终的分类结果为N个结果中获得票数最多的那个类别；
S54、判断随机森林分类器网络是否达到预设的期望值，若是则转至步骤S57，否则转至步骤S55；
S55、设置随机森林分类器中迭代次数、学习率以及决策树棵树为可调超参数，给定初值及调整范围，然后对预设模型参数采用贝叶斯优化算法进行优化；
S56、贝叶斯优化后，更新随机森林模型的权重和偏差，并转入S52；
S57、保存训练好的随机森林分类器网络参数，并转至步骤S61。
在本发明一实施例中，所述步骤S6具体实现如下：
S61、基于步骤S5训练的随机森林分类器的网络参数，建立配电变压器故障诊断模型，并将测试集输入到该配电变压器故障诊断模型；
S62、实现配电变压器故障诊断。
相较于现有技术，本发明具有以下有益效果：
1、本发明利用贝叶斯优化算法自动调整网络的超参数，获得网络的最优参数，这样极大节省了人力和时间；
2、将堆栈自编码器与随机森林结合起来，利用堆栈自编码器逐层学习数据、提取数据特征，再利用随机森林对已提取的特征进行深入学习，从而可将配电变压器振动信号中蕴含的更为丰富的故障特征提取出来，进而显著提高配电变压器故障诊断的准确性；
3、本发明建立的配电变压器故障诊断系统可应用于工程实践中，其同时实现振动信号的自动特征提取和故障诊断并具有良好的鲁棒性和突出的诊断性能。
附图说明
图1为本发明配电变压器故障诊断过程示意图。
图2为堆栈自编码构建流程图。
图3为贝叶斯优化流程。
图4为随机森林建模过程。
图5为本发明配电变压器故障诊断方法的诊断过程流程图。
图6为本发明堆栈自编码器模型提取特征展示图。
具体实施方式
下面结合附图，对本发明的技术方案进行具体说明。
如图1所示，本发明提供了一种特征自动挖掘和参数自动寻优的配电变压器故障诊断方法，包括如下步骤：
步骤S1、将振动信号采集装置安装在配电变压器上，并从振动信号采集装置中采集配电变压器运行时的振动波形；
步骤S2、构建基于二次调优后的堆栈自编码器的配电变压器故障特征提取模型；
步骤S3、利用二次调优后的堆栈自编码器提取振动信号特征向量yn，并对相应的特征打上标签，建立一个包含正常及各类故障的数据库；
步骤S4、分割数据集，即将数据集按X1：X2的比例拆分为训练集和测试集；
步骤S5、利用训练集训练随机森林分类器；
步骤S6、基于步骤S5训练的随机森林分类器的网络参数，建立配电变压器故障诊断模型，实现配电变压器故障诊断。
如图5所示，本发明配电变压器故障诊断方法的诊断过程，具体如下：
S1：将振动信号采集装置安装在配电变压器上，并从振动信号采集装置中采集配电变压器运行时的振动波形；
S2：构建基于堆栈自编码器的配电变压器故障特征提取模型；
S21：创建一个堆栈自编码器模型，利用采集得到的振动数据逐层学习数据，提取数据特征，得到每层自编码器的偏差和权重；再通过前馈反向神经网络来缩小网络损失值，从而不断更新权重和偏差，得到堆栈自编码器中每个编码器的权重和偏差，具体实现方法如下：
假设一训练样本x＝[x1,x2,...,xn]T∈Rn，由公式(1)可将其向隐含层进行映射，映射后的隐含层表示为y1＝[y1,y2,...,ym]T∈Rm；
y1＝fθ(x)＝se(Wx+b) (1)
其中，se(·)为编码阶段的激活函数；θ＝{W,b}为编码阶段的参数，其中包含权值矩阵W∈Rm×n、偏置矩阵b∈Rm；
接着，采用解码器公式(2)将隐含表示变换回x的重构向量
其中，sd(·)为解码阶段的激活函数；
为解码阶段的参数，包含权值矩阵
偏置矩阵
近似于x的程度越高，自编码器网络的性能越优良；因此，要不断最小化
和x之间的误差，即最小化损失函数：
对于网络参数
的调整，同BP神经网络的调节机制，采用梯度下降法来最小化平均重构误差；本申请采用的损失函数为：
其中，
为交叉熵，用于度量输入向量与重构向量之间的差异，两者差异越小，交叉熵的值越小；
随后将此编码器提取出的隐含层表示特征y1作为下一个编码器的输入，同样利用公式(1)可得y1映射后的隐含层表示为y2；再采用公式(2)可得y2的重构向量
此时同样利用公式(4)来不断最小化y1和
之间的误差；利用此步骤依次对n个编码器进行训练，可得每个自编码器的权重和偏差；该过程的流程图如图2所示。
S22：设置堆栈自编码器特征提取网络中的迭代次数和学习率为可调超参数，并设定参数可调范围，将这些可调超参数传递给S23的贝叶斯优化网络中。
S23：利用贝叶斯优化进行超参数自动寻优。
贝叶斯优化算法是基于贝叶斯定理而提出的，表示为：
其中，D＝{(x1,y1),(x2,y2),...,(xn,yn)}表示已采集样本点，p(f)为f的先验概率，p(f|D)为f的后验概率，f表示目标函数或模型中的参数，后验概率分布是通过已采集样本点对先验概率进行修正后目标函数或模型中参数的置信度；本申请采用的概率代理模型为高斯过程，该过程是多元高斯概率分布的泛化，由协方差函数和均值函数构成，表达式为：
y＝gp(m(x),k(x,x')) (6)
其中，m(x)为均值函数，k(x,x')为协方差函数，该过程对一系列离散数据点(xi,yi)进行函数拟合，得到一个逼近真实值的目标函数；
常见的采集函数有上置信边界(UCB)、改进期望(EI)和改进概率(PI)三种，本申请采用PI函数，其形式为：
其中f(X)为目标函数，f(X+)是目前为止最优的目标函数，μ(x)是f(X)后验分布的均值，δ(x)是协差矩阵，ν是权衡系数，用于避免陷入局部最优；采集函数求解的过程就是找到一个新的x，使PI(x)最大；
贝叶斯优化与其他调参方法不同之处在于，它考虑了上一次参数信息，通过不断地更新先验，使迭代次数更少、速度更快。实现流程如图3所示。
S24：判断超参数寻优的目标函数是否达到预设的期望值，若是则转至S25，否则转至S23；
S25：得到基于堆栈自编码器的配电变压器故障初始特征提取模型；
S26：堆栈自编码器二次调优；
二次调优是指将堆栈自编码器看作一个多层神经网络，利用有标签的训练样本集，对该神经网络的权重和偏差进行调整。二次调优可以很好地提高堆栈自编码器的性能。本研究采用对误差进行反向传播来进行二次调优，具体步骤如下：
步骤一：根据式(8)和式(9)计算各层的激励响应；
Z(l)＝W(l-1)a(l-1)+b(l-1) (8)
a(l)＝f(Z(l)) (9)
步骤二：根据式(10)计算残差；
步骤三：根据式(11)、式(12)计算权重和偏差的更新公式；
步骤四：根据式(13)、(14)计算权重和偏差的改变值；
其中ΔW(l)和W(l)是大小相同的矩阵，Δb(l)和b(l)是行数相同的向量，初始值都设置为0；
步骤五：根据式(15)、(16)得到更新后的权重和偏差；
其中，m表示样本数量，μ为学习速率。
S3：利用二次调优后的堆栈自编码器提取振动信号特征向量yn，并对相应的特征打上标签，建立一个包含正常及各类故障的数据库；
S4：分割数据集；
S41：将数据集按X1：X2的比例拆分为训练集和测试集；
S42：将训练集随机打乱后，作为S5的数据输入；
S43：将测试集作为S6的数据输入；
S5：利用训练集训练随机森林分类器，随机森林算法思想如图4所示；
S51：利用Bootstrap抽样从原始数据中有放回的随机抽取X个样本构成N个子训练集，每个样本的样本容量与原始数据集一样；
S52：根据N个子训练集训练生成N个决策树模型，从而得到N个分类结果；
S53：通过投票表决得出最终的分类结果，此处投票表决服从“少数服从多数”的原则，即最终的分类结果为N个结果中获得票数最多的那个类别。
S54：判断随机森林分类器网络是否达到预设的期望值，若是则转至S57，否则转至S55；
S55：设置随机森林分类器中迭代次数、学习率以及决策树棵树为可调超参数，给定初值及调整范围，然后对预设模型参数进行S23中的贝叶斯优化算法。
S56：贝叶斯优化后，更新随机森林分类器的权重和偏差，并转入执行S52；
S57：保存训练好的随机森林分类器网络参数，并转至S61。
S6：故障诊断阶段；
S61：结合S57保存的网络参数，建立配电变压器故障诊断模型，并将S42制作好的测试样本输入到该模型；
S62：实现配电变压器故障诊断。
下面结合具体实例对本发明作进一步说明：
在本实施例中，选用型号为S11-M-315/10的配电变压器；选用型号为LC0156A的压电加速度传感器；传感器采样频率设置为10kHz，采样时间为1s，因此每个振动信号中包含10000个采样值。每类故障的样本各采集200个，其中150个用作训练集，50个用作测试集。
在本实例中，研究的配电变压器的数据集为实测的配电变压器振动数据集，如表1所示。
表1 实测的配电变压器振动数据集
由表1可知，实测数据中配电变压器状态包括：类别1：绕组变形、类别2：绕组螺丝松动、类别3：正常。鉴于配电变压器故障数据难以获取，本实施例中，通过人为制造故障的方式获得类别1和类别2的故障数据。
为了验证本发明方法的有效性，构建一个3步(1000→500→250→100)降维的堆栈自编码器模型。根据上文S2步骤所述方法使用测得的数据对该模型进行预训练，然后使用有标签数据对该模型进行二次调整；再根据S3、S4步骤的方法将提取到的特征值建立数据库，并按3:1的比例分为训练集和测试集；然后根据S5步骤来训练RF分类器模型。最终得到网络的最优参数，并以此作为故障诊断模型的参数，其参数如表2所示。振动信号经堆栈自编码器提取特征信息后，可得到100维的特征向量，由于100维向量无法展示，此处利用T-SNE将100维降为2维进行展示。图6中展示了堆栈自编码器模型有无二次调优时提取到的特征量。
表2 网络参数
图6(a)为无调整情况下特征分布图。每个类别已经形成多个簇聚集在一起，但同一个类别的簇分散在不同的区域，图中还存在许多散乱分布的点游离于分布空间中，这种分布结果给分类器带来一定的挑战。图6(b)是经过调整后的效果。可以看出每个类别聚集程度相对于图6(a)有明显的改善，每个类别已经占据一定的区域，但图上依然有个别点存在发散现象，原因是采样环境存在噪声干扰，导致这些数据质量不好。总之，使用本发明所提的堆栈自编码器网络提取的特征分类更为清晰明确，效果较为突出。因而，将此特征值重新构造数据库，划分训练集和测试集，使用训练集训练随机森林分类器，最终利用测试集对该分类模型进行验证，得到的测试集准确度可达96.67％。结果证明，本发明所提方法对配电变压器故障信号特征提取及分类具有良好的效果，证实了该方法的有效性。
以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超出本发明技术方案的范围时，均属于本发明的保护范围。
Claims (5)
1.一种特征自动挖掘和参数自动寻优的配电变压器故障诊断方法，其特征在于，包括如下步骤：
步骤S1、将振动信号采集装置安装在配电变压器上，并从振动信号采集装置中采集配电变压器运行时的振动波形；
步骤S2、构建基于二次调优后的堆栈自编码器的配电变压器故障特征提取模型；
步骤S3、利用二次调优后的堆栈自编码器提取振动信号特征向量yn，并对相应的特征打上标签，建立一个包含正常及各类故障的数据库；
步骤S4、分割数据集，即将数据集按X1：X2的比例拆分为训练集和测试集；
步骤S5、利用训练集训练随机森林分类器；
步骤S6、基于步骤S5训练的随机森林分类器的网络参数，建立配电变压器故障诊断模型，实现配电变压器故障诊断。
2.根据权利要求1所述的特征自动挖掘和参数自动寻优的配电变压器故障诊断方法，其特征在于，所述步骤S2具体实现如下：
S21、创建一个堆栈自编码器模型，利用测量得到的振动数据逐层学习数据，提取数据特征，得到每层自编码器的偏差和权重；再通过前馈反向神经网络来缩小网络损失值，从而不断更新权重和偏差，得到堆栈自编码器中每个编码器的权重和偏差，具体实现方式如下：
假设一训练样本x＝[x1,x2,...,xn]T∈Rn，由公式(1)可将其向隐含层进行映射，映射后的隐含层表示为y1＝[y1,y2,...,ym]T∈Rm；
y1＝fθ(x)＝se(Wx+b) (1)
其中，se(·)为编码阶段的激活函数；θ＝{W,b}为编码阶段的参数，其中包含权值矩阵W∈Rm×n、偏置矩阵b∈Rm；
接着，采用解码器公式(2)将隐含表示变换回x的重构向量
其中，sd(·)为解码阶段的激活函数；
为解码阶段的参数，包含权值矩阵
偏置矩阵
近似于x的程度越高，自编码器网络的性能越优良；因此，要不断最小化
和x之间的误差，即最小化损失函数：
对于网络参数
的调整，同BP神经网络的调节机制，采用梯度下降法来最小化平均重构误差；具体的损失函数可表示为：
事实上，
是交叉熵，用于度量输入向量与重构向量之间的差异，两者差异越小，交叉熵的值越小；
随后将此编码器提取出的隐含层表示特征y1作为下一个编码器的输入，同样利用公式(1)可得y1映射后的隐含层表示为y2；再采用公式(2)可得y2的重构向量
此时同样利用公式(4)来不断最小化y1和
之间的误差；利用此步骤依次对n个编码器进行训练，可得每个自编码器的权重和偏差；
S22、设置堆栈自编码器特征提取网络中的迭代次数和学习率为可调超参数，并设定可调超参数调整范围，将这些可调超参数传递给步骤S23；
S23、利用贝叶斯优化进行超参数自动寻优，贝叶斯优化公式表示为：
其中，D＝{(x1,y1),(x2,y2),...,(xn,yn)}表示已采集样本点，p(f)为f的先验概率，p(f|D)为f的后验概率，f表示目标函数或模型中的参数，后验概率分布是通过已采集样本点对先验概率进行修正后目标函数或模型中参数的置信度；采用高斯过程作为概率代理模型，该过程是多元高斯概率分布的泛化，由协方差函数和均值函数构成，表达式为：
y＝gp(m(x),k(x,x')) (6)
其中，m(x)为均值函数，k(x,x')为协方差函数，该过程对一系列离散数据点(xi,yi)进行函数拟合，得到一个逼近真实值的目标函数；
采用PI函数作为采集函数，其形式为：
其中f(X)为目标函数，f(X+)是目前为止最优的目标函数，μ(x)是f(X)后验分布的均值，δ(x)是协差矩阵，ν是权衡系数，用于避免陷入局部最优；采集函数求解的过程就是找到一个新的x，使PI(x)最大；
S24、判断超参数寻优的目标函数是否达到预设的期望值，若是则转至S25，否则转至S23；
S25、得到基于堆栈自编码器的配电变压器故障初始特征提取模型；
S26、对堆栈自编码器进行二次调优，得到基于二次调优后的堆栈自编码器的配电变压器故障特征提取模型。
3.根据权利要求2所述的特征自动挖掘和参数自动寻优的配电变压器故障诊断方法，其特征在于，所述步骤S26具体实现如下：
S261、根据式(8)和式(9)计算各层的激励响应；
Z(l)＝W(l-1)a(l-1)+b(l-1) (8)
a(l)＝f(Z(l)) (9)
S262、根据式(10)计算残差；
S263、根据式(11)、式(12)计算权重和偏差的更新公式；
S264、根据式(13)、(14)计算权重和偏差的改变值；
其中ΔW(l)和W(l)是大小相同的矩阵，Δb(l)和b(l)是行数相同的向量，初始值都设置为0；
S265、根据式(15)、(16)得到更新后的权重和偏差；
其中，m表示样本数量，μ为学习速率。
4.根据权利要求1所述的特征自动挖掘和参数自动寻优的配电变压器故障诊断方法，其特征在于，所述步骤S5具体实现如下：
S51、利用Bootstrap抽样从原始数据中有放回的随机抽取X个样本构成N个子训练集，每个样本的样本容量与原始数据集一样；
S52、根据N个子训练集训练生成N个决策树模型，从而得到N个分类结果；
S53、通过投票表决得出最终的分类结果，此处投票表决服从“少数服从多数”的原则，即最终的分类结果为N个结果中获得票数最多的那个类别；
S54、判断随机森林分类器网络是否达到预设的期望值，若是则转至步骤S57，否则转至步骤S55；
S55、设置随机森林分类器中迭代次数、学习率以及决策树棵树为可调超参数，给定初值及调整范围，然后对预设模型参数采用贝叶斯优化算法进行优化；
S56、贝叶斯优化后，更新随机森林分类器网络中的权重和偏差；并转入步骤S52；
S57、保存训练好的随机森林分类器网络参数，并转至步骤S6。
5.根据权利要求1所述的特征自动挖掘和参数自动寻优的配电变压器故障诊断方法，其特征在于，所述步骤S6具体实现如下：
S61、基于步骤S5训练的随机森林分类器的网络参数，建立配电变压器故障诊断模型，并将测试集输入到该配电变压器故障诊断模型；
S62、实现配电变压器故障诊断。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011178015.9A CN112327219B (zh) | 2020-10-29 | 2020-10-29 | 特征自动挖掘和参数自动寻优的配电变压器故障诊断方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011178015.9A CN112327219B (zh) | 2020-10-29 | 2020-10-29 | 特征自动挖掘和参数自动寻优的配电变压器故障诊断方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112327219A true CN112327219A (zh) | 2021-02-05 |
| CN112327219B CN112327219B (zh) | 2024-03-12 |
Family
ID=74297777
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011178015.9A Active CN112327219B (zh) | 2020-10-29 | 2020-10-29 | 特征自动挖掘和参数自动寻优的配电变压器故障诊断方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112327219B (zh) |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113033674A (zh) * | 2021-03-25 | 2021-06-25 | 安徽理工大学 | 一种贝叶斯优化随机森林算法的苹果多光谱图像无损检测方法 |
| CN113049035A (zh) * | 2021-03-12 | 2021-06-29 | 辽宁工程技术大学 | 一种基于物联网的变压器状态监测系统 |
| CN113177594A (zh) * | 2021-04-29 | 2021-07-27 | 浙江大学 | 基于贝叶斯优化的pca-极限随机树的空调故障诊断方法 |
| CN113219374A (zh) * | 2021-03-05 | 2021-08-06 | 国网江苏省电力有限公司徐州供电分公司 | 一种变压器电压故障检测方法 |
| CN113283546A (zh) * | 2021-07-20 | 2021-08-20 | 深圳市佳运通电子有限公司 | 加热炉完整性管理集控装置的炉况异常报警方法及系统 |
| CN113313049A (zh) * | 2021-06-11 | 2021-08-27 | 北京百度网讯科技有限公司 | 超参数的确定方法、装置、设备、存储介质以及计算机程序产品 |
| CN113505445A (zh) * | 2021-09-10 | 2021-10-15 | 北京航空航天大学 | 基于序贯随机森林的实时故障诊断方法和系统 |
| CN113566953A (zh) * | 2021-09-23 | 2021-10-29 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种柔壁喷管的在线监测方法 |
| CN113570070A (zh) * | 2021-09-23 | 2021-10-29 | 深圳市信润富联数字科技有限公司 | 流式数据采样与模型更新方法、装置、系统与存储介质 |
| CN113935406A (zh) * | 2021-09-27 | 2022-01-14 | 苏州大学 | 基于对抗流模型的机械设备无监督故障诊断方法 |
| CN114295981A (zh) * | 2022-01-05 | 2022-04-08 | 湖南师范大学 | 一种基于rn网络电流多特征形态识别的电机故障诊断方法 |
| CN114325236A (zh) * | 2021-12-28 | 2022-04-12 | 广东电网有限责任公司 | 一种基于频谱熵和随机森林的配电网故障识别方法及系统 |
| CN114357663A (zh) * | 2022-03-15 | 2022-04-15 | 中国长江三峡集团有限公司 | 一种训练齿轮箱故障诊断模型方法、齿轮箱故障诊断方法 |
| CN117056734A (zh) * | 2023-10-12 | 2023-11-14 | 山东能源数智云科技有限公司 | 基于数据驱动的设备故障诊断模型的构建方法及装置 |
| CN117092445A (zh) * | 2023-10-19 | 2023-11-21 | 盛隆电气集团有限公司 | 基于大数据的配电系统的故障检测方法及系统 |
| CN117150243A (zh) * | 2023-10-27 | 2023-12-01 | 湘江实验室 | 一种基于故障影响解耦网络的故障隔离与估计方法 |
| CN117571901A (zh) * | 2023-11-17 | 2024-02-20 | 承德神源太阳能发电有限公司 | 一种光伏电站变压器故障预警及检修的方法、系统及设备 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102221655A (zh) * | 2011-06-16 | 2011-10-19 | 河南省电力公司济源供电公司 | 基于随机森林模型的电力变压器故障诊断方法 |
| US20160245686A1 (en) * | 2015-02-23 | 2016-08-25 | Biplab Pal | Fault detection in rotor driven equipment using rotational invariant transform of sub-sampled 3-axis vibrational data |
| CN107167702A (zh) * | 2017-05-04 | 2017-09-15 | 国网福建省电力有限公司 | 一种配电网馈线故障类型识别方法及装置 |
| CN108153982A (zh) * | 2017-12-26 | 2018-06-12 | 哈尔滨工业大学 | 基于堆叠自编码深度学习网络的航空发动机修后性能预测方法 |
| CN109613351A (zh) * | 2018-11-21 | 2019-04-12 | 北京国网富达科技发展有限责任公司 | 一种变压器的故障诊断方法、设备及系统 |
| CN109669087A (zh) * | 2019-01-31 | 2019-04-23 | 国网河南省电力公司 | 一种基于多源信息融合的电力变压器故障诊断方法 |
| CN110110768A (zh) * | 2019-04-24 | 2019-08-09 | 西安电子科技大学 | 基于并行特征学习和多分类器的滚动轴承故障诊断方法 |
| CN110348489A (zh) * | 2019-06-19 | 2019-10-18 | 西安理工大学 | 一种基于自编码网络的变压器局部放电模式识别方法 |
| CN110596492A (zh) * | 2019-09-17 | 2019-12-20 | 昆明理工大学 | 一种基于粒子群算法优化随机森林模型的变压器故障诊断方法 |
| CN111751763A (zh) * | 2020-06-08 | 2020-10-09 | 武汉大学 | 一种基于GSMallat-NIN-CNN网络的电力变压器绕组故障诊断方法 |
-
2020
- 2020-10-29 CN CN202011178015.9A patent/CN112327219B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102221655A (zh) * | 2011-06-16 | 2011-10-19 | 河南省电力公司济源供电公司 | 基于随机森林模型的电力变压器故障诊断方法 |
| US20160245686A1 (en) * | 2015-02-23 | 2016-08-25 | Biplab Pal | Fault detection in rotor driven equipment using rotational invariant transform of sub-sampled 3-axis vibrational data |
| CN107167702A (zh) * | 2017-05-04 | 2017-09-15 | 国网福建省电力有限公司 | 一种配电网馈线故障类型识别方法及装置 |
| CN108153982A (zh) * | 2017-12-26 | 2018-06-12 | 哈尔滨工业大学 | 基于堆叠自编码深度学习网络的航空发动机修后性能预测方法 |
| CN109613351A (zh) * | 2018-11-21 | 2019-04-12 | 北京国网富达科技发展有限责任公司 | 一种变压器的故障诊断方法、设备及系统 |
| CN109669087A (zh) * | 2019-01-31 | 2019-04-23 | 国网河南省电力公司 | 一种基于多源信息融合的电力变压器故障诊断方法 |
| CN110110768A (zh) * | 2019-04-24 | 2019-08-09 | 西安电子科技大学 | 基于并行特征学习和多分类器的滚动轴承故障诊断方法 |
| CN110348489A (zh) * | 2019-06-19 | 2019-10-18 | 西安理工大学 | 一种基于自编码网络的变压器局部放电模式识别方法 |
| CN110596492A (zh) * | 2019-09-17 | 2019-12-20 | 昆明理工大学 | 一种基于粒子群算法优化随机森林模型的变压器故障诊断方法 |
| CN111751763A (zh) * | 2020-06-08 | 2020-10-09 | 武汉大学 | 一种基于GSMallat-NIN-CNN网络的电力变压器绕组故障诊断方法 |
Non-Patent Citations (1)
| Title |
|---|
| 武立平: "基于振动信号指标能量的变压器机械故障检测", 《电测与仪表》, vol. 57, no. 6, pages 126 - 131 * |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113219374A (zh) * | 2021-03-05 | 2021-08-06 | 国网江苏省电力有限公司徐州供电分公司 | 一种变压器电压故障检测方法 |
| CN113219374B (zh) * | 2021-03-05 | 2024-04-30 | 国网江苏省电力有限公司徐州供电分公司 | 一种变压器电压故障检测方法 |
| CN113049035B (zh) * | 2021-03-12 | 2022-05-27 | 辽宁工程技术大学 | 一种基于物联网的变压器状态监测系统 |
| CN113049035A (zh) * | 2021-03-12 | 2021-06-29 | 辽宁工程技术大学 | 一种基于物联网的变压器状态监测系统 |
| CN113033674A (zh) * | 2021-03-25 | 2021-06-25 | 安徽理工大学 | 一种贝叶斯优化随机森林算法的苹果多光谱图像无损检测方法 |
| CN113177594A (zh) * | 2021-04-29 | 2021-07-27 | 浙江大学 | 基于贝叶斯优化的pca-极限随机树的空调故障诊断方法 |
| CN113313049A (zh) * | 2021-06-11 | 2021-08-27 | 北京百度网讯科技有限公司 | 超参数的确定方法、装置、设备、存储介质以及计算机程序产品 |
| CN113283546A (zh) * | 2021-07-20 | 2021-08-20 | 深圳市佳运通电子有限公司 | 加热炉完整性管理集控装置的炉况异常报警方法及系统 |
| CN113283546B (zh) * | 2021-07-20 | 2021-11-02 | 深圳市佳运通电子有限公司 | 加热炉完整性管理集控装置的炉况异常报警方法及系统 |
| CN113505445A (zh) * | 2021-09-10 | 2021-10-15 | 北京航空航天大学 | 基于序贯随机森林的实时故障诊断方法和系统 |
| CN113566953A (zh) * | 2021-09-23 | 2021-10-29 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种柔壁喷管的在线监测方法 |
| CN113570070A (zh) * | 2021-09-23 | 2021-10-29 | 深圳市信润富联数字科技有限公司 | 流式数据采样与模型更新方法、装置、系统与存储介质 |
| CN113935406A (zh) * | 2021-09-27 | 2022-01-14 | 苏州大学 | 基于对抗流模型的机械设备无监督故障诊断方法 |
| CN114325236A (zh) * | 2021-12-28 | 2022-04-12 | 广东电网有限责任公司 | 一种基于频谱熵和随机森林的配电网故障识别方法及系统 |
| CN114325236B (zh) * | 2021-12-28 | 2023-08-18 | 广东电网有限责任公司 | 一种基于频谱熵和随机森林的配电网故障识别方法及系统 |
| CN114295981A (zh) * | 2022-01-05 | 2022-04-08 | 湖南师范大学 | 一种基于rn网络电流多特征形态识别的电机故障诊断方法 |
| CN114357663A (zh) * | 2022-03-15 | 2022-04-15 | 中国长江三峡集团有限公司 | 一种训练齿轮箱故障诊断模型方法、齿轮箱故障诊断方法 |
| CN117056734A (zh) * | 2023-10-12 | 2023-11-14 | 山东能源数智云科技有限公司 | 基于数据驱动的设备故障诊断模型的构建方法及装置 |
| CN117056734B (zh) * | 2023-10-12 | 2024-02-06 | 山东能源数智云科技有限公司 | 基于数据驱动的设备故障诊断模型的构建方法及装置 |
| CN117092445A (zh) * | 2023-10-19 | 2023-11-21 | 盛隆电气集团有限公司 | 基于大数据的配电系统的故障检测方法及系统 |
| CN117150243A (zh) * | 2023-10-27 | 2023-12-01 | 湘江实验室 | 一种基于故障影响解耦网络的故障隔离与估计方法 |
| CN117150243B (zh) * | 2023-10-27 | 2024-01-30 | 湘江实验室 | 一种基于故障影响解耦网络的故障隔离与估计方法 |
| CN117571901A (zh) * | 2023-11-17 | 2024-02-20 | 承德神源太阳能发电有限公司 | 一种光伏电站变压器故障预警及检修的方法、系统及设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112327219B (zh) | 2024-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112327219B (zh) | 特征自动挖掘和参数自动寻优的配电变压器故障诊断方法 | |
| CN109685138B (zh) | 一种xlpe电力电缆局部放电类型识别方法 | |
| CN105550700B (zh) | 一种基于关联分析和主成分分析的时间序列数据清洗方法 | |
| CN103728551B (zh) | 一种基于级联集成分类器的模拟电路故障诊断方法 | |
| CN109948597B (zh) | 一种高压断路器机械故障诊断方法 | |
| CN110929847A (zh) | 一种基于深度卷积神经网络的换流变压器故障诊断方法 | |
| Tong et al. | Detection and classification of transmission line transient faults based on graph convolutional neural network | |
| CN109214416A (zh) | 一种基于深度学习的多维信息融合变压器故障诊断方法 | |
| CN101231672A (zh) | 基于改进型bp神经网络的模拟电路软故障诊断方法 | |
| CN110082082B (zh) | 一种基于振动信号主成分分析法的gis状态识别方法 | |
| CN110298085A (zh) | 基于XGBoost和随机森林算法的模拟电路故障诊断方法 | |
| Wang et al. | Few-shot transfer learning with attention mechanism for high-voltage circuit breaker fault diagnosis | |
| CN112130026B (zh) | 基于小波包能量熵和dbn的mmc-hvdc单极接地故障定位方法 | |
| CN111680725A (zh) | 基于重构贡献的气体传感器阵列多故障隔离算法 | |
| CN116154972B (zh) | 一种分布式电网电能质量监测方法及系统 | |
| CN116380445A (zh) | 基于振动波形的设备状态诊断方法及相关装置 | |
| CN116089857A (zh) | 一种基于ceemdan-dbn的变压器故障识别方法 | |
| CN112486137A (zh) | 有源配电网故障特征库构建方法、系统及故障诊断方法 | |
| CN115600088A (zh) | 一种基于振动信号的配电变压器故障诊断方法 | |
| CN113514743B (zh) | 一种基于多维特征的gis局部放电模式识别系统构建方法 | |
| CN113379116A (zh) | 基于聚类和卷积神经网络的台区线损预测方法 | |
| CN110110784B (zh) | 一种基于变压器相关运行数据的变压器故障辨识方法 | |
| CN117520809A (zh) | 一种基于EEMD-KPCA-CNN-BiLSTM的变压器故障诊断方法 | |
| CN116244600A (zh) | 一种gis间歇性放电模式识别模型的构建方法、系统及设备 | |
| CN115840116A (zh) | 一种基于gan-cnn的开关柜局放故障模式识别方法 |
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 |