CN115481657A - 一种基于电信号的风力发电机组通信滑环故障诊断方法 - Google Patents

一种基于电信号的风力发电机组通信滑环故障诊断方法 Download PDF

Info

Publication number
CN115481657A
CN115481657A CN202211026575.1A CN202211026575A CN115481657A CN 115481657 A CN115481657 A CN 115481657A CN 202211026575 A CN202211026575 A CN 202211026575A CN 115481657 A CN115481657 A CN 115481657A
Authority
CN
China
Prior art keywords
slip ring
domain
signal
fault diagnosis
value
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
Application number
CN202211026575.1A
Other languages
English (en)
Inventor
张轶东
刘博�
傅望安
张伟平
段周朝
李霄
薛文超
石壮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Thermal Power Research Institute Co Ltd
Clean Energy Branch of Huaneng Zhejiang Energy Development Co Ltd
Original Assignee
Xian Thermal Power Research Institute Co Ltd
Clean Energy Branch of Huaneng Zhejiang Energy Development Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Thermal Power Research Institute Co Ltd, Clean Energy Branch of Huaneng Zhejiang Energy Development Co Ltd filed Critical Xian Thermal Power Research Institute Co Ltd
Priority to CN202211026575.1A priority Critical patent/CN115481657A/zh
Publication of CN115481657A publication Critical patent/CN115481657A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial 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]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开一种基于电信号的风力发电机组通信滑环故障诊断方法,包括以下步骤:采集不同工况下通信滑环的电信号;对信号进行分解;利用多域特征向量对粒子群优化最小二乘支持向量机的故障诊断模型进行训练,实现故障诊断。本发明中,在信号降噪过程中,通过F‑GWO优化的VMD降噪方法将故障信号进行分解,减少模态混叠现象的出现;在信号特征提取过程中,根据自适应VMD的分解特点,通过计算奇异值及排列熵值分别提取了信号时域、频域和时频域特征,获取能够全面反映信号特点的多域特征向量;在故障诊断过程中,通过多域特征向量对PSO‑LSSVM进行训练,提高了故障诊断的准确度和速度,为风力发电机组通信滑环故障诊断提供了方向。

Description

一种基于电信号的风力发电机组通信滑环故障诊断方法
技术领域
本发明属于风力发电技术领域,具体涉及一种基于电信号的风力发电机组通信滑环故障诊断方法。
背景技术
风电滑环是用于风力发电机组的机舱和轮毂之间传输电能动力和电气信号的装置,滑环的故障现象主要表现为通信故障、环道烧毁、编码器故障等。其中,通信故障是一类较常见的故障,而其问题的根本原因是刷针与环道的接触不良,从而影响信号的传输。从问题本质上分析,主要分为二类:一类为滑环自身缺陷,导致接触不良,二类为未及时维护导致接触不良。因通信滑环故障所导致的风力发电机组运行不稳定现象屡有发生,造成巨大的经济损失。因此,通信滑环故障诊断具有重要的现实意义。
近年来,有学者进行了关于风力发电机组滑环故障诊断的研究,Chen等根据滑环、碳刷的热成像来判断碳刷-滑环的运行状态;Tang等根据滑环、碳刷的振动信号对碳刷-滑环的故障特性进行了研究;WURFEL等从电刷滑环装置运行温度的角度来分析判断其运行状况。但以上种种诊断及监测方法,主要还是从振动信号、图像监测等方面入手,依赖传感器、红外检测等外部设备,这些设备所需安装的数量多、所处环境恶劣时容易被损坏,因此,存在成本高、可适应性差的缺陷,且实际中正常运行的风力发电机组本身就具有复杂的振动特性,尤其在风速随机不稳定变化的冲击下,机组的振动情况更加复杂,使振动监测的准确度不高,在信号分析过程中虽然可以实现对噪声信号的处理,但均只考虑了时域的特征量,诊断结果的可靠性较差。
发明内容
为解决现有技术中存在的技术问题,本发明的目的在于提供一种基于电信号的风力发电机组通信滑环故障诊断方法。
为实现上述目的,达到上述技术效果,本发明采用的技术方案为:
一种基于电信号的风力发电机组通信滑环故障诊断方法,包括以下步骤:首先,采集不同工况下通信滑环的电信号,选取通信滑环故障时的电信号作为实验数据;随后,对信号进行分解,提取信号的时域、频域和时频域特征,得到能够全面反映信号特点的多域特征向量T;利用多域特征向量T对粒子群优化最小二乘支持向量机的故障诊断模型进行训练,完成故障模式识别,实现故障诊断。
进一步的,包括以下步骤:
1)采集不同工况下通信滑环的电信号,选取通信滑环故障时的电信号作为实验数据;
2)对信号进行自适应变分模态分解;
3)由步骤2)可得信号的K个固有模态分量IMF,即得到了原信号的时域信息,求取各个IMF时域矩阵的奇异值,得到时域矩阵奇异值构成的K维特征向量T1
4)对K个IMF进行快速傅里叶变换,得到原信号的频域信息,并求取各个IMF频域矩阵的奇异值,得到频域矩阵奇异值构成的K维特征向量T2
5)计算K个IMF时频域的排列熵Hpe,得到Hpe构成的K维特征向量T3
6)分别对时域、频域和时频域特征向量采用主成分分析法进行分析,选取贡献度最大的主成分分量,作为各域的特征量;
7)综合上述降维后的各域特征量,再将时域、频域、时频域的特征量综合,即得到能够全面反映信号特点的多域特征向量T;
8)采用PSO对LSSVM参数进行寻优,搭建粒子群优化最小二乘支持向量机的故障诊断模型PSO-LSSVM,根据设计计算的多域特征向量T对PSO-LSSVM进行训练;
9)在训练好的PSO-LSSVM中进行故障模式识别,即可实现故障诊断。
进一步的,步骤2)中,对信号进行自适应变分模态分解的步骤包括:
2.1)初始化F-GWO参数,设置迭代次数为n,灰狼ω的数量为m,将(K,α)设置为灰狼的位置向量,其中,K为VMD的模态数;
2.2)根据每匹灰狼的位置向量即(K,α),对信号进行自适应VMD,计算固有模态分量IMF的平均包络熵AEE;
2.3)更新最小平均包络熵MAEE,并把前三小的AEE的位置向量分别给α,β和δ,α代表最适解,β代表第二个最佳解决方案,δ代表第三个最佳解决方案;
2.4)根据α,β和δ灰狼的位置,更新m个灰狼ω的位置;
2.5)重复步骤2.2)至步骤2.4),直到迭代次数达到n;
2.6)输出α灰狼的位置向量,即得最佳分解参数
Figure BDA0003815996660000021
进一步的,步骤2.1)中,K的寻优范围为2-12,α的寻优范围为800-5000。
进一步的,步骤3)中,各个IMF时域矩阵的奇异值的求取步骤包括:
3.1)计算矩阵AAT的特征值;
3.2)对其特征值开方,即得到矩阵A的奇异值;
奇异值越大,其所代表的能量越大,小奇异值的奇异向量对原矩阵影响很小,大奇异值所对应的奇异向量更能代表原先的矩阵。
进一步的,步骤5)中,K个IMF时频域的排列熵Hpe的计算步骤包括;
5.1)将得到的各IMF分量看作一个时间序列X(i),i=1,2,…,N,进行相空间重构,得到矩阵Y;
5.2)将每一个重构分量按照升序重新排列,得到向量中各元素位置的列索引构成一组符号序列:S(l)={j1,j2,…,jm},l=1,2,…,n,且n≤m!;
5.3)计算每一种符号序列出现的次数,其除以m!种不同的符号序列出现的总次数作为该符号序列出现的概率,即{P1,P2,…,Pn};
5.4)时间序列X的排列熵Hpe的计算公式为:
Figure BDA0003815996660000031
5.5)排列熵的最大值为ln(d!),将排列熵值进行归一化处理,即:
Figure BDA0003815996660000032
进一步的,步骤5.1)中,矩阵Y为:
Figure BDA0003815996660000033
其中,m为嵌入维数,t为延迟时间,n=N-(m-1)t,矩阵Y中的每一行都是一个重构分量,共有n个重构分量。
进一步的,步骤8)中,根据设计计算的多域特征向量T对PSO-LSSVM进行训练,训练步骤包括:
8.1)将数据样本点划分为训练样本集和测试样本集两部分;
8.2)对最大迭代次数、最大和最小的加权因子wmax和wmin、学习因子c1,c2与径向基核函数中的宽度σ2进行初始化设置;
8.3)根据c1、c2和σ2的值,计算出目标函数的适应度值,目标函数为训练正确率;
8.4)将各个粒子当下获得的适应度值与该粒子之前的最优适应度值比较大小,如果目前获得的数值更优,则将其值作为个体的最佳位置pBest;
8.5)将各个粒子当前计算出的适应度值同该种群最优位置的适应度值进行比较,如果前者更优,则将其值作为群体的最优位置gBest;
8.6)根据如下两式更新定义的所有粒子的位置和速度大小,满足搜寻更匹配的参数c1、c2和σ2的值;
V=wV+c1rand()(pBest-p)+c2rand()(gBest-p)
p=p+V
其中,rand为(0,1)之间的随机函数值;c1,c2分别被称为学习因子;w被称为加权因子;V被称为粒子运行的速度;p被称为粒子当前所处的位置;
8.7)重复步骤8.3)至步骤8.6),直到满足最佳的匹配值。
进一步的,步骤8.3)中,正确率的计算公式如下:
accuracy=(TP+TN)/(P+N)
其中,TP是被正确地划分为正例的个数,即实际为正例且被分类器划分为正例的实例数(样本数);TN是被正确地划分为负例的个数,即实际为负例且被分类器划分为负例的实例数;P是实际为正的样本数;N是实际为负的样本数。
与现有技术相比,本发明的有益效果为:
本发明提出了一种基于电信号的风力发电机组通信滑环故障诊断方法,包括以下步骤:采集不同工况下通信滑环的电信号,采用快速灰狼算法(F-GWO)优化的变分模态分解(VMD)降噪方法对信号进行处理,在一定程度减少了模态混叠问题;接着,根据自适应VMD信号处理方法的特点,计算各模态时频域矩阵的奇异值及各模态的排列熵,并根据主成分分析法(PCA)简化特征信息,凸显出故障特征,构建多域特征向量;接着利用多域特征向量对PSO-LSSVM进行训练,提高了故障诊断的准确性。本发明中,从电信号入手采集故障信号,提高了信号采集精度,降低了采集成本,且不需要额外装设设备,还可以避免因检测设备故障造成的故障信号不稳定、可靠性差的问题;在信号降噪过程中,通过快速灰狼算法(F-GWO)优化的变分模态分解(VMD)降噪方法将故障信号进行分解来实现降噪处理,降低了复杂度高和非线性强的时间序列非平稳性,减少了模态混叠现象的出现;在信号特征提取过程中,根据自适应VMD的分解特点,通过计算奇异值及排列熵值分别提取了信号时域、频域和时频域特征,利用PCA简化特征信息,选取代表各域的主要特征量构成多域特征向量,提高故障诊断的准确度和效率,而且,所提取的时域、频域和时频域特征能够反应多域的特点,蕴含更加丰富的信息,能够使故障特征更具有独创性;在故障诊断过程中,采用PSO对LSSVM参数进行寻优,通过多域特征向量对粒子群算法(PSO)优化的最小二乘支持向量机(LSSVM)进行训练,本发明通过多域特征提取和智能算法优化参数,提高了故障诊断的准确度和速度,为风力发电机组通信滑环故障诊断提供了方向。
附图说明
图1为本发明的流程图;
图2为本发明的自适应VMD方法的流程图。
具体实施方式
下面对本发明进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
以下给出一个或多个方面的简要概述以提供对这些方面的基本理解。此概述不是所有构想到的方面的详尽综览,并且既非旨在指认出所有方面的关键性或决定性要素亦非试图界定任何或所有方面的范围。其唯一的目的是要以简化形式给出一个或多个方面的一些概念以为稍后给出的更加详细的描述之序。
如图1-2所示,一种基于电信号的风力发电机组通信滑环故障诊断方法,包括以下步骤:
1)采集不同工况下通信滑环的电信号,选取通信滑环故障时的电信号作为实验数据;
2)对原信号即原信号是步骤1)中选来作为实验数据的电信号进行自适应变分模态分解(VMD);
步骤2)中,对VMD分解中需要人为设定的参数,以最小平均包络熵作为适应度函数,引入快速灰狼算法(F-GWO)对VMD进行寻优,得到自适应VMD方法;
3)由步骤2)可得原信号的K个固有模态分量(IMF),即得到了原信号的时域信息,求取各个IMF时域矩阵的奇异值,得到时域矩阵奇异值构成的K维特征向量T1
4)对K个IMF进行快速傅里叶变换,得到原信号的频域信息,并求取各个IMF频域矩阵的奇异值,得到频域矩阵奇异值构成的K维特征向量T2
5)计算K个IMF时频域的排列熵(PE),得到PE构成的K维特征向量T3
6)分别对时域、频域和时频域特征向量采用主成分分析法(PCA)进行分析,选取贡献度最大的主成分分量,作为各域的特征量;
7)综合上述降维后的各域特征量,再将时域、频域、时频域的特征量综合,即得到能够全面反映信号特点的多域特征向量T;
8)采用PSO对LSSVM参数进行寻优,搭建了粒子群优化最小二乘支持向量机的故障诊断模型PSO-LSSVM,根据设计计算的多域特征向量T对PSO-LSSVM进行训练;
9)在训练好的PSO-LSSVM中进行故障模式识别,即可实现故障诊断。
实施例1
如图1-2所示,一种基于电信号的风力发电机组通信滑环故障诊断方法,包括以下步骤:
1)以某型号风力发电机组为分析对象,采集不同工况下通信滑环的电信号,选取通信滑环故障时的电信号作为实验数据,故障信号与正常信号的区别在于频率不同。
2)对原信号进行自适应变分模态分解;
需要说明的是,故障的特征频率幅值较小,且其往往与环境噪声一起混杂在复杂信号中,直接对电信号进行快速傅里叶变化很难观察到故障特征频率,为实现对原始信号的高分辨率时频分析并提取故障特征分量,需要详细了解组成复杂信号的各模态分量,将非线性、非平稳的时序信号分解为相对线性、平稳的各模态分量,清晰地显示各模态分量的规律,从而明确原始复杂信号的组成成分,因此,本步骤对原信号进行自适应变分模态分解(VMD),具体步骤包括:
2.1)初始化F-GWO参数,设置迭代次数为n,灰狼ω的数量为m,将(K,α)设置为灰狼的位置向量,其中,K为VMD的模态数,其寻优范围为2-12,α的寻优范围为800-5000;
2.2)根据每匹灰狼的位置向量即(K,α),对信号进行自适应VMD,计算固有模态分量IMF的平均包络熵(AEE);
2.3)更新最小平均包络熵(MAEE),并把前三小的AEE的位置向量分别给α,β和δ,α代表最适解,β代表第二个最佳解决方案,δ代表第三个最佳解决方案;
2.4)根据α,β和δ灰狼的位置,更新m个灰狼ω的位置;
2.5)重复步骤2.2)至步骤2.4),直到迭代次数达到n;
2.6)输出α灰狼的位置向量,即得最佳分解参数
Figure BDA0003815996660000071
3)由步骤2)可得原信号的K个固有模态分量(IMF),即得到了原信号的时域信息,求取各个IMF时域矩阵的奇异值,将各奇异值拼在一起得到时域矩阵奇异值构成的K维特征向量T1;其中,各个IMF时域矩阵的奇异值的求取步骤包括:
3.1)计算矩阵AAT的特征值;
3.2)对其特征值开方,即得到矩阵A的奇异值;
奇异值越大,其所代表的能量越大,小奇异值的奇异向量对原矩阵影响很小,大奇异值所对应的奇异向量更能代表原先的矩阵。因此,可选取前若干个奇异值和奇异向量,能够用少量的矩阵相乘来近乎无损表示出原矩阵的信息,如此通过奇异值分解(SVD)能得到代表矩阵最本质变化的特征元素。基于此,选取奇异值作为电信号时域和频域信息的特征量。
4)对K个IMF进行快速傅里叶变换,得到原信号的频域信息,并求取各个IMF频域矩阵的奇异值,得到频域矩阵奇异值构成的K维特征向量T2
5)故障信号和正常信号在分解后所得到的各分量复杂程度不同,不同类型的故障信号复杂程度也不同,通过计算K个IMF时频域的排列熵(PE),可以衡量各自的复杂度,得到PE构成的K维特征向量T3,直观放大了信号中的微小差别,各IMF的PE计算步骤如下;
5.1)将得到的各IMF分量看作一个时间序列X(i),i=1,2,…,N,进行相空间重构,得到矩阵Y为:
Figure BDA0003815996660000081
其中,m为嵌入维数,t为延迟时间,n=N-(m-1)t,矩阵Y中的每一行都是一个重构分量,共有n个重构分量;
5.2)将每一个重构分量按照升序重新排列,得到向量中各元素位置的列索引构成一组符号序列:S(l)={j1,j2,…,jm},l=1,2,…,n,且n≤m!;
5.3)计算每一种符号序列出现的次数,其除以m!种不同的符号序列出现的总次数作为该符号序列出现的概率,即{P1,P2,…,Pn};
5.4)时间序列X的排列熵的计算公式为:
Figure BDA0003815996660000082
5.5)排列熵的最大值为ln(d!),将排列熵值进行归一化处理,即:
Figure BDA0003815996660000083
6)分别对时域、频域和时频域特征向量采用主成分分析法(PCA)进行分析,选取贡献度最大的主成分分量,作为各域的特征量;
7)综合上述降维后的各域特征向量,再将时域、频域、时频域的特征向量综合,即得到能够全面反映信号特点的多域特征向量T;
8)采用PSO对LSSVM参数进行寻优,搭建了粒子群优化最小二乘支持向量机的故障诊断模型PSO-LSSVM,根据设计计算的多域特征向量T对PSO-LSSVM进行训练,具体步骤包括:
8.1)将数据样本点划分为训练样本集和测试样本集两部分;
8.2)对最大迭代次数、最大和最小的加权因子wmax和wmin、学习因子c1,c2与径向基核函数中的宽度σ2进行初始化设置;
8.3)根据C和σ2的值,计算出目标函数的适应度值,目标函数为训练正确率,计算公式如下所示:
accuracy=(TP+TN)/(P+N)
其中,TP是被正确地划分为正例的个数,即实际为正例且被分类器划分为正例的实例数(样本数);TN是被正确地划分为负例的个数,即实际为负例且被分类器划分为负例的实例数;P是实际为正的样本数;N是实际为负的样本数;
8.4)将各个粒子当下获得的适应度值与该粒子之前的最优适应度值比较大小,如果目前获得的数值更优,则将其值作为个体的最佳位置pBest;
8.5)将各个粒子当前计算出的适应度值同该种群最优位置的适应度值进行比较,如果前者更优,则将其值作为群体的最优位置gBest;
8.6)根据如下两式更新定义的所有粒子的位置和速度大小,满足搜寻更匹配的参数c1,c2和σ2的值;
V=wV+c1rand()(pBest-p)+c2rand()(gBest-p)
p=p+V
其中,rand为(0,1)之间的随机函数值;c1,c2分别被称为学习因子;w被称为加权因子;V被称为粒子运行的速度;p被称为粒子当前所处的位置;
8.7)重复步骤8.3)至步骤8.6),直到满足最佳的匹配值。
9)将测试样本集的特征向量输入到训练好的PSO-LSSVM中,在训练好的PSO-LSSVM中进行故障模式识别,输出诊断结果,实现故障诊断。
本发明未具体描述的部分或结构采用现有技术或现有产品即可,在此不做赘述。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (9)

1.一种基于电信号的风力发电机组通信滑环故障诊断方法,其特征在于,包括以下步骤:首先,采集不同工况下通信滑环的电信号,选取通信滑环故障时的电信号作为实验数据;随后,对信号进行分解,提取信号的时域、频域和时频域特征,得到能够全面反映信号特点的多域特征向量T;利用多域特征向量T对粒子群优化最小二乘支持向量机的故障诊断模型进行训练,完成故障模式识别,实现故障诊断。
2.根据权利要求1所述的一种基于电信号的风力发电机组通信滑环故障诊断方法,其特征在于,包括以下步骤:
1)采集不同工况下通信滑环的电信号,选取通信滑环故障时的电信号作为实验数据;
2)对信号进行自适应变分模态分解;
3)由步骤2)可得信号的K个固有模态分量IMF,即得到了原信号的时域信息,求取各个IMF时域矩阵的奇异值,得到时域矩阵奇异值构成的K维特征向量T1
4)对K个IMF进行快速傅里叶变换,得到原信号的频域信息,并求取各个IMF频域矩阵的奇异值,得到频域矩阵奇异值构成的K维特征向量T2
5)计算K个IMF时频域的排列熵Hpe,得到Hpe构成的K维特征向量T3
6)分别对时域、频域和时频域特征向量采用主成分分析法进行分析,选取贡献度最大的主成分分量,作为各域的特征量;
7)综合上述降维后的各域特征量,再将时域、频域、时频域的特征量综合,即得到能够全面反映信号特点的多域特征向量T;
8)采用PSO对LSSVM参数进行寻优,搭建粒子群优化最小二乘支持向量机的故障诊断模型PSO-LSSVM,根据设计计算的多域特征向量T对PSO-LSSVM进行训练;
9)在训练好的PSO-LSSVM中进行故障模式识别,即可实现故障诊断。
3.根据权利要求1所述的一种基于电信号的风力发电机组通信滑环故障诊断方法,其特征在于,步骤2)中,对信号进行自适应变分模态分解的步骤包括:
2.1)初始化F-GWO参数,设置迭代次数为n,灰狼ω的数量为m,将(K,α)设置为灰狼的位置向量,其中,K为VMD的模态数;
2.2)根据每匹灰狼的位置向量即(K,α),对信号进行自适应VMD,计算固有模态分量IMF的平均包络熵AEE;
2.3)更新最小平均包络熵MAEE,并把前三小的AEE的位置向量分别给α,β和δ,α代表最适解,β代表第二个最佳解决方案,δ代表第三个最佳解决方案;
2.4)根据α,β和δ灰狼的位置,更新m个灰狼ω的位置;
2.5)重复步骤2.2)至步骤2.4),直到迭代次数达到n;
2.6)输出α灰狼的位置向量,即得最佳分解参数
Figure FDA0003815996650000021
4.根据权利要求3所述的一种基于电信号的风力发电机组通信滑环故障诊断方法,其特征在于,步骤2.1)中,K的寻优范围为2-12,α的寻优范围为800-5000。
5.根据权利要求1所述的一种基于电信号的风力发电机组通信滑环故障诊断方法,其特征在于,步骤3)中,各个IMF时域矩阵的奇异值的求取步骤包括:
3.1)计算矩阵AAT的特征值;
3.2)对其特征值开方,即得到矩阵A的奇异值;
奇异值越大,其所代表的能量越大,小奇异值的奇异向量对原矩阵影响很小,大奇异值所对应的奇异向量更能代表原先的矩阵。
6.根据权利要求1所述的一种基于电信号的风力发电机组通信滑环故障诊断方法,其特征在于,步骤5)中,K个IMF时频域的排列熵Hpe的计算步骤包括;
5.1)将得到的各IMF分量看作一个时间序列X(i),i=1,2,…,N,进行相空间重构,得到矩阵Y;
5.2)将每一个重构分量按照升序重新排列,得到向量中各元素位置的列索引构成一组符号序列:S(l)={j1,j2,…,jm},l=1,2,…,n,且n≤m!;
5.3)计算每一种符号序列出现的次数,其除以m!种不同的符号序列出现的总次数作为该符号序列出现的概率,即{P1,P2,…,Pn};
5.4)时间序列X的排列熵Hpe的计算公式为:
Figure FDA0003815996650000022
5.5)排列熵的最大值为ln(d!),将排列熵值进行归一化处理,即:
Figure FDA0003815996650000031
7.根据权利要求6所述的一种基于电信号的风力发电机组通信滑环故障诊断方法,其特征在于,步骤5.1)中,矩阵Y为:
Figure FDA0003815996650000032
其中,m为嵌入维数,t为延迟时间,n=N-(m-1)t,矩阵Y中的每一行都是一个重构分量,共有n个重构分量。
8.根据权利要求1所述的一种基于电信号的风力发电机组通信滑环故障诊断方法,其特征在于,步骤8)中,根据设计计算的多域特征向量T对PSO-LSSVM进行训练,训练步骤包括:
8.1)将数据样本点划分为训练样本集和测试样本集两部分;
8.2)对最大迭代次数、最大和最小的加权因子wmax和wmin、学习因子c1,c2与径向基核函数中的宽度σ2进行初始化设置;
8.3)根据c1、c2和σ2的值,计算出目标函数的适应度值,目标函数为训练正确率;
8.4)将各个粒子当下获得的适应度值与该粒子之前的最优适应度值比较大小,如果目前获得的数值更优,则将其值作为个体的最佳位置pBest;
8.5)将各个粒子当前计算出的适应度值同该种群最优位置的适应度值进行比较,如果前者更优,则将其值作为群体的最优位置gBest;
8.6)根据如下两式更新定义的所有粒子的位置和速度大小,满足搜寻更匹配的参数c1、c2和σ2的值;
V=wV+c1rand()(pBest-p)+c2rand()(gBest-p)
p=p+V
其中,rand为(0,1)之间的随机函数值;w被称为加权因子;V被称为粒子运行的速度;p被称为粒子当前所处的位置;
8.7)重复步骤8.3)至步骤8.6),直到满足最佳的匹配值。
9.根据权利要求1所述的一种基于电信号的风力发电机组通信滑环故障诊断方法,其特征在于,步骤8.3)中,正确率的计算公式如下:
accuracy=(TP+TN)/(P+N)
其中,TP是被正确地划分为正例的个数,即实际为正例且被分类器划分为正例的实例数(样本数);TN是被正确地划分为负例的个数,即实际为负例且被分类器划分为负例的实例数;P是实际为正的样本数;N是实际为负的样本数。
CN202211026575.1A 2022-08-25 2022-08-25 一种基于电信号的风力发电机组通信滑环故障诊断方法 Pending CN115481657A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211026575.1A CN115481657A (zh) 2022-08-25 2022-08-25 一种基于电信号的风力发电机组通信滑环故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211026575.1A CN115481657A (zh) 2022-08-25 2022-08-25 一种基于电信号的风力发电机组通信滑环故障诊断方法

Publications (1)

Publication Number Publication Date
CN115481657A true CN115481657A (zh) 2022-12-16

Family

ID=84421956

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211026575.1A Pending CN115481657A (zh) 2022-08-25 2022-08-25 一种基于电信号的风力发电机组通信滑环故障诊断方法

Country Status (1)

Country Link
CN (1) CN115481657A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116340728A (zh) * 2023-04-07 2023-06-27 平顶山天安煤业股份有限公司 输电杆塔的倾斜测量信号的处理方法及计算机存储介质
CN116838947A (zh) * 2023-06-30 2023-10-03 中国人民解放军总医院第二医学中心 一种供节氧系统的故障诊断方法及供节氧系统

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116340728A (zh) * 2023-04-07 2023-06-27 平顶山天安煤业股份有限公司 输电杆塔的倾斜测量信号的处理方法及计算机存储介质
CN116838947A (zh) * 2023-06-30 2023-10-03 中国人民解放军总医院第二医学中心 一种供节氧系统的故障诊断方法及供节氧系统
CN116838947B (zh) * 2023-06-30 2024-02-13 中国人民解放军总医院第二医学中心 一种供节氧系统的故障诊断方法及供节氧系统

Similar Documents

Publication Publication Date Title
CN109187025B (zh) 一种集成kelm的滚动轴承剩余使用寿命预测方法
Wang et al. A novel weighted sparse representation classification strategy based on dictionary learning for rotating machinery
CN115481657A (zh) 一种基于电信号的风力发电机组通信滑环故障诊断方法
CN111914883B (zh) 一种基于深度融合网络的主轴轴承状态评估方法及装置
CN111060337B (zh) 一种基于专家系统的运转设备实时故障诊断方法
CN105678343B (zh) 基于自适应加权组稀疏表达的水电机组噪声异常诊断方法
CN111881594B (zh) 一种核动力设备的非平稳信号状态监测方法及系统
CN113240016A (zh) 一种基于深度学习的滚动轴承故障诊断方法
CN112069930A (zh) 提升gis设备故障诊断准确率的振动信号处理方法及装置
CN110595778B (zh) 一种基于mmf与igra风电机组轴承故障诊断方法
CN112507479B (zh) 一种基于流形学习和softmax的石油钻机健康状态评估方法
CN114091504A (zh) 一种基于生成对抗网络的旋转机械小样本故障诊断方法
CN111898644A (zh) 一种无故障样本下航天液体发动机健康状态智能识别方法
CN117272143A (zh) 基于格拉姆角场与残差网络的输电线故障辨识方法和装置
CN109726770A (zh) 一种模拟电路故障测试诊断方法
CN116861320A (zh) 基于短时傅里叶同步压缩变换的转子故障诊断方法
CN116089857A (zh) 一种基于ceemdan-dbn的变压器故障识别方法
Zhao et al. Research on constructing a degradation index and predicting the remaining useful life for rolling element bearings of complex equipment
CN117009794B (zh) 一种基于无监督子域自适应的机器故障诊断方法及系统
CN117992777A (zh) 一种基于图像融合的电机电气故障诊断方法
CN117592543A (zh) 一种基于自监督学习的航空发动机滚动轴承故障诊断方法
CN115618202A (zh) 一种基于流形嵌入和关键特征选取的机械故障诊断方法
CN115561575A (zh) 基于多维矩阵轮廓的海上风电场电气异常状态辨别方法
CN114593917A (zh) 一种基于三元组模型的小样本轴承故障诊断方法
Jiang et al. Research on the effect of wind turbine bearing fault diagnosis method based on multi-feature calculation and Bayesian optimized machine learning method

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