CN106768933A - 一种基于遗传算法的风电机组叶片故障诊断方法 - Google Patents
一种基于遗传算法的风电机组叶片故障诊断方法 Download PDFInfo
- Publication number
- CN106768933A CN106768933A CN201611097540.1A CN201611097540A CN106768933A CN 106768933 A CN106768933 A CN 106768933A CN 201611097540 A CN201611097540 A CN 201611097540A CN 106768933 A CN106768933 A CN 106768933A
- Authority
- CN
- China
- Prior art keywords
- genetic algorithm
- wind power
- generation unit
- power generation
- fault
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Genetics & Genomics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Physiology (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种基于遗传算法的可变K均值风电机组叶片故障诊断方法,其主要针对典型的故障情况,对采集的风电机组叶片故障数据进行预处理,得到故障特征信息;运用主成分分析法PCA对数据进行标准化,并提取出导致风电机组叶片故障的主要特征信息;将所提取的主要故障特征信息,利用基于遗传算法的可变K均值聚类算法,对训练样本进行聚类分析。本发明结合了遗传算法,在传统K均值算法的基础上进行K值的计算选定,缩短了故障信息分析时间,提高了数据聚类分析精确度,降低了成本,从而使故障诊断的过程变得简单可靠,有效可行。
Description
技术领域
本发明涉及故障诊断方法领域,具体地说,特别涉及到一种基于遗传算法的风电机组叶片故障诊断方法。
背景技术
近年来,环境污染和能源危机的态势已经愈演愈烈,急需一种绿色环保新能源的出现。风电作为一种可再生、无排放的新能源,已经越来越受到全球各国的重视,这也使得风力发电技术在近些年来得到了蓬勃发展。在我国,风电场的大量建立和投入使用,将是响应政府改善能源结构、应对气候变化的号召,实现传统能源向新能源逐步过渡的具体体现。
除了可再生、无污染以外,风力发电在我国还具有资源丰富,占地面积广,单机容量小等优点。但是,由于风电场通常都建设在偏远地区,还存在技术条件差,环境恶劣等问题。同时,大规模风电机组的投入应用也使其安全稳定性得到了人们的高度重视。由于风电机组长期工作在恶劣的自然环境中,使得风场风况复杂多变,容易引发各种故障的产生。
因此,故障诊断环节是必不可少的。而风机的叶片作为风能捕捉的主要组件之一,除了受到恶劣自然环境的影响,还容易受到因负荷波动和快速变化而产生的疲劳损伤,叶轮偏载等故障,从而导致捕风效率下降,严重者甚至无法继续投入使用。
目前,对于叶片故障诊断,大多采用神经网络的方法。虽然神经网络算法具有准确度高,并行分布处理能力强,分布存储及学习能力强,能充分逼近复杂的非线性关系等优点,但是同时也存在着明显的劣势,如需要提前设定大量阈值和权值,学习时间过长,样本数量庞大时网络过于复杂等。
发明内容
本发明的目的在于针对现有技术中的不足,提供一种简单易行,成本低廉,过程时间短的新型风电机组叶片故障诊断方法,即基于遗传算法的可变K均值风电机组叶片故障诊断方法,从而使故障诊断的过程变得简单可行提高风电机组叶片的安全、可靠性。
本发明所解决的技术问题可以采用以下技术方案来实现:
一种基于遗传算法的风电机组叶片故障诊断方法,包括如下步骤:
1)采集风电机组叶片的故障数据,然后对采集到的风电机组叶片的故障数据进行预处理,得到故障特征信息,并对故障特征信息进行标准化处理;
2)运用主成分分析法PCA对标准化处理后的故障特征信息进行分析,提取出导致风电机组叶片故障的主要特征信息;
3)将提取出的主要特征信息,利用基于遗传算法的可变K均值聚类算法,对要特征信息进行聚类分析,获得风电故障分析结果;
4)将上述风电故障分析结果与故障数据库或专家知识库进行比对,得到最终故障类型的分析结果,最后将分析结果显示在人机交互界面上。
进一步的,所述主成分分析法PCA的步骤如下:
a1)数据标准化,公式为其中则新的输入变量的数据矩阵为
a2)求相关系数矩阵,
其中,rij(i,j=1,2,......,p)为标准化后变量与的相关系数,其计算公式为
a3)求系数矩阵R的特征值和特征向量,记为λi=(i=1,2,...,p),对应特征向量为ei(i=1,2,...,p);
a4)计算主成份贡献率、累计贡献率并确定主成分个数。总方差中属于第i个主成分zi的比例称为贡献率:前i个主成分的贡献率之和称为w1,w2,......,wi的累计贡献率一般取累计贡献率达到85%至95%的特征值所对应的λ1,λ2,...,λm的第1主成分,第2主成分,……第m(m≤p)个主成分。
a5)计算基于主成分的输入变量数据矩阵z,
进一步的,所述聚类分析的步骤如下:
b1)利用遗传算法的过程,生成初始种群,设定终止条件,即遗传算法迭代次数;
b2)利用K均值算法对筛选出来的个体进行优化,并用优化后的个体取代原来的个体;
b3)对种群中的个体进行选择、交叉、变异操作,并在完成后重新计算K值;
b4)重复步骤b2)和b3)至满足终止条件。
与现有技术相比,本发明的有益效果如下:
1、对采集的风电机组叶片故障数据进行标准化预处理,减少了故障特征信息数值上的差异,使聚类分析的过程变得简便。
2、通过PCA来提取故障信息中的主要特征信息,降低了数据的复杂度,实现了数据不相关和模型简化。
3、使用基于遗传算法的可变K均值聚类算法,提高了故障特征信息聚类分析的有效性,可行性和效率。
附图说明
图1为本发明所述的风电机组叶片故障诊断方法流程示意图。
图2为本发明所述的可变K均值算法的流程示意图。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。
参见图1和图2,本发明所述的yiz基于遗传算法的可变K均值风电机组叶片故障诊断方法,具体实施过程如下:
步骤一为故障数据监测与采集,即先进行数据采集与标准化,即针对典型的故障情况,对风电机组叶片故障数据进行监测和采集,并对故障数据进行标准化和初步预处理,得到故障特征信息。
步骤二为运用PCA,即主成分分析,提取出导致风电机组叶片故障的主要特征信息。PCA具体操作步骤为:
数据标准化,公式为其中
则新的输入变量的数据矩阵为
(2)求相关系数矩阵,
其中rij(i,j=1,2,......,p)为标准化后变量与的相关系数,其计算公式为
(3)求系数矩阵R的特征值和特征向量,记为λi=(i=1,2,...,p),对应特征向量为ei(i=1,2,...,p)
(4)计算主成份贡献率、累计贡献率并确定主成分个数。总方差中属于第i个主成分zi的比例称为贡献率:前i个主成分的贡献率之和称为w1,w2,……,wi的累计贡献率一般取累计贡献率达到85%至95%的特征值所对应的λ1,λ2,...,λm的第1主成分,第2主成分,……第m(m≤p)个主成分。
(5)计算基于主成分的输入变量数据矩阵z,
步骤三为运用故障聚类模型,即利用基于遗传算法的可变K均值算法进行聚类分析。
步骤四为故障解释机制,即将风电故障分析结果与故障数据库或专家知识库进行比对,得到最终故障类型的分析结果。
步骤五为将结果显示在人机交互界面上。
上述过程为对故障数据聚类分析主过程。
图2是故障信息聚类分析模型,即基于遗传算法的可变K均值算法具体操作图。步骤如下:
步骤一为用PCA进行故障信息的数据预处理。
步骤二为进行染色体编码设计和适应度函数设计。由于风电机组叶片故障具有多维度,数量大的特点,因此,为结合实际情况,本发明采用基于聚类中心的浮点数编码。染色体结构为其中K为基因的长度,为随机产生的聚类中心的个数,为染色体第l+1位的基因,它有n个维度。xl∈(1,2,...,c),l=1,2,...,K。适应度函数是用来评价个体,区别个体优劣的标准。本发明的适应度函数采用聚类中的目标函数。
步骤三为初始化信息设定,包括种群生成和终止条件设定。其中,终止条件设定即为规定遗传算法迭代次数,超过即终止。
步骤四为进行选择,交叉和变异操作,并结合K均值算法进行聚类分析。
(1)选择。为了保证适应度函数越大的染色体保留到下一代的概率越大,本文采用了遗传算法中经典的轮盘赌算法选择算子。具体做法如下:
1)进行K均值操作,更新染色体;
2)根据适应度函数,计算当前种群染色体的适应度,并记录其中适应度最大的个体;
3)再根据个体的适应度值来计算各个体的选择概率。
(2)交叉。随机从种群中选择M/2对染色体作为父代,对每对染色体,产生随机数P(0<P<1),当<Pc时,在两个染色体i,j中随机产生交叉点,对交叉点后的基因进行交叉运算,并重新计算交叉后的染色体长度,直到所有个体进行过交叉为止。
(3)变异。对每一个个体,每一个基因位,产生随机数P,当P<Pm时,对该基因位进行随机变异运算,生成下一代种群。(其中Pm为交叉概率)。
步骤五为检查染色体长度,当染色体长度大于最大聚类个数时,将这条染色体从种群中删去,否则执行步骤七。
步骤六为以变异并通过检查后产生的染色体基因为中心,重新确K值,进行K均值聚类,然后用K均值聚类的结果来取代原来的染色体中的基因。同时,判断是否满足终止条件,即是否达到预先设定的遗传算法迭代次数。若未达到,则重复到步骤四进行选择,交叉和变异,用K均值优化个体。若到达,输出结果。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (3)
1.一种基于遗传算法的风电机组叶片故障诊断方法,其特征在于,包括如下步骤:
1)采集风电机组叶片的故障数据,然后对采集到的风电机组叶片的故障数据进行预处理,得到故障特征信息,并对故障特征信息进行标准化处理;
2)运用主成分分析法PCA对标准化处理后的故障特征信息进行分析,提取出导致风电机组叶片故障的主要特征信息;
3)将提取出的主要特征信息,利用基于遗传算法的可变K均值聚类算法,对要特征信息进行聚类分析,获得风电故障分析结果;
4)将上述风电故障分析结果与故障数据库或专家知识库进行比对,得到最终故障类型的分析结果,最后将分析结果显示在人机交互界面上。
2.根据权利要求1所述的基于遗传算法的风电机组叶片故障诊断方法,其特征在于,所述主成分分析法PCA的步骤如下:
a1)数据标准化,公式为其中则新的输入变量的数据矩阵为
a2)求相关系数矩阵,
其中,rij(i,j=1,2,......,p)为标准化后变量与的相关系数,其计算公式为
a3)求系数矩阵R的特征值和特征向量,记为λi=(i=1,2,...,p),对应特征向量为ei(i=1,2,...,p);
a4)计算主成份贡献率、累计贡献率并确定主成分个数。总方差中属于第i个主成分zi的比例称为贡献率:前i个主成分的贡献率之和称为w1,w2,......,wi的累计贡献率一般取累计贡献率达到85%至95%的特征值所对应的λ1,λ2,...,λm的第1主成分,第2主成分,……第m(m≤p)个主成分。
a5)计算基于主成分的输入变量数据矩阵z,
3.根据权利要求1所述的基于遗传算法的风电机组叶片故障诊断方法,其特征在于,所述聚类分析的步骤如下:
b1)利用遗传算法的过程,生成初始种群,设定终止条件,即遗传算法迭代次数;
b2)利用K均值算法对筛选出来的个体进行优化,并用优化后的个体取代原来的个体;
b3)对种群中的个体进行选择、交叉、变异操作,并在完成后重新计算K值;
b4)重复步骤b2)和b3)至满足终止条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611097540.1A CN106768933A (zh) | 2016-12-02 | 2016-12-02 | 一种基于遗传算法的风电机组叶片故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611097540.1A CN106768933A (zh) | 2016-12-02 | 2016-12-02 | 一种基于遗传算法的风电机组叶片故障诊断方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106768933A true CN106768933A (zh) | 2017-05-31 |
Family
ID=58884313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611097540.1A Pending CN106768933A (zh) | 2016-12-02 | 2016-12-02 | 一种基于遗传算法的风电机组叶片故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106768933A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107178477A (zh) * | 2017-07-10 | 2017-09-19 | 龙源(北京)风电工程技术有限公司 | 基于深度自编码模型的风电机组叶片故障监控方法及*** |
CN107256546A (zh) * | 2017-05-23 | 2017-10-17 | 上海海事大学 | 基于pca卷积池化softmax的海流机叶片附着物故障诊断方法 |
CN109657795A (zh) * | 2018-12-12 | 2019-04-19 | 华中科技大学 | 一种基于属性选择的硬盘故障预测方法 |
CN110147808A (zh) * | 2019-03-26 | 2019-08-20 | 张锐明 | 一种新型电池成组筛选方法 |
CN110259648A (zh) * | 2019-07-05 | 2019-09-20 | 河北工业大学 | 一种基于优化K-means聚类的风机叶片故障诊断方法 |
CN115791142A (zh) * | 2023-02-09 | 2023-03-14 | 中国航发四川燃气涡轮研究院 | 轴向限位叶片结构及配置方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101782475A (zh) * | 2010-02-08 | 2010-07-21 | 天津工业大学 | 基于风力发电机组振动的叶片故障诊断方法 |
CN102175449A (zh) * | 2011-03-18 | 2011-09-07 | 天津工业大学 | 基于风力机应变能响应的叶片故障诊断方法 |
CN102507205A (zh) * | 2011-10-27 | 2012-06-20 | 西北工业大学 | 一种检测航空发动机风扇叶片颤振故障的方法 |
CN102706562A (zh) * | 2012-06-01 | 2012-10-03 | 河海大学 | 基于无线传感器网络的海上风电机组状态监测***及方法 |
CN102944416A (zh) * | 2012-12-06 | 2013-02-27 | 南京匹瑞电气科技有限公司 | 基于多传感器信号融合技术的风电机组叶片故障诊断方法 |
CN102944418A (zh) * | 2012-12-11 | 2013-02-27 | 东南大学 | 一种风电机组叶片故障诊断方法 |
CN103323461A (zh) * | 2013-06-14 | 2013-09-25 | 上海大学 | 一种非接触式的风力发电机叶片运动在线检测方法 |
CN104515677A (zh) * | 2015-01-12 | 2015-04-15 | 华北电力大学 | 风电机组叶片故障诊断与状态监控*** |
CN106124982A (zh) * | 2016-06-14 | 2016-11-16 | 都城绿色能源有限公司 | 一种风电机组的自动专家综合故障诊断***及诊断方法 |
-
2016
- 2016-12-02 CN CN201611097540.1A patent/CN106768933A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101782475A (zh) * | 2010-02-08 | 2010-07-21 | 天津工业大学 | 基于风力发电机组振动的叶片故障诊断方法 |
CN102175449A (zh) * | 2011-03-18 | 2011-09-07 | 天津工业大学 | 基于风力机应变能响应的叶片故障诊断方法 |
CN102507205A (zh) * | 2011-10-27 | 2012-06-20 | 西北工业大学 | 一种检测航空发动机风扇叶片颤振故障的方法 |
CN102706562A (zh) * | 2012-06-01 | 2012-10-03 | 河海大学 | 基于无线传感器网络的海上风电机组状态监测***及方法 |
CN102944416A (zh) * | 2012-12-06 | 2013-02-27 | 南京匹瑞电气科技有限公司 | 基于多传感器信号融合技术的风电机组叶片故障诊断方法 |
CN102944418A (zh) * | 2012-12-11 | 2013-02-27 | 东南大学 | 一种风电机组叶片故障诊断方法 |
CN103323461A (zh) * | 2013-06-14 | 2013-09-25 | 上海大学 | 一种非接触式的风力发电机叶片运动在线检测方法 |
CN104515677A (zh) * | 2015-01-12 | 2015-04-15 | 华北电力大学 | 风电机组叶片故障诊断与状态监控*** |
CN106124982A (zh) * | 2016-06-14 | 2016-11-16 | 都城绿色能源有限公司 | 一种风电机组的自动专家综合故障诊断***及诊断方法 |
Non-Patent Citations (2)
Title |
---|
王辉 等: ""基于遗传算法的可变K均值聚类算法研究"", 《中国信息经济学会2012学术年会暨第四届博士生论坛》 * |
甄立敬: ""风力发电机组发电机和齿轮箱故障诊断方法研究"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107256546A (zh) * | 2017-05-23 | 2017-10-17 | 上海海事大学 | 基于pca卷积池化softmax的海流机叶片附着物故障诊断方法 |
CN107178477A (zh) * | 2017-07-10 | 2017-09-19 | 龙源(北京)风电工程技术有限公司 | 基于深度自编码模型的风电机组叶片故障监控方法及*** |
CN109657795A (zh) * | 2018-12-12 | 2019-04-19 | 华中科技大学 | 一种基于属性选择的硬盘故障预测方法 |
CN109657795B (zh) * | 2018-12-12 | 2022-06-21 | 华中科技大学 | 一种基于属性选择的硬盘故障预测方法 |
CN110147808A (zh) * | 2019-03-26 | 2019-08-20 | 张锐明 | 一种新型电池成组筛选方法 |
CN110259648A (zh) * | 2019-07-05 | 2019-09-20 | 河北工业大学 | 一种基于优化K-means聚类的风机叶片故障诊断方法 |
CN110259648B (zh) * | 2019-07-05 | 2020-10-09 | 河北工业大学 | 一种基于优化K-means聚类的风机叶片故障诊断方法 |
CN115791142A (zh) * | 2023-02-09 | 2023-03-14 | 中国航发四川燃气涡轮研究院 | 轴向限位叶片结构及配置方法 |
CN115791142B (zh) * | 2023-02-09 | 2023-06-13 | 中国航发四川燃气涡轮研究院 | 轴向限位叶片结构及配置方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106768933A (zh) | 一种基于遗传算法的风电机组叶片故障诊断方法 | |
CN106712713B (zh) | 一种光伏电站的监测***以及监测和异常定位方法 | |
CN110929918B (zh) | 一种基于CNN和LightGBM的10kV馈线故障预测方法 | |
CN102707256B (zh) | 基于BP-AdaBoost神经网络的电子式电能表故障诊断方法 | |
CN109670553B (zh) | 基于自适应神经模糊推理***的光伏阵列故障诊断方法 | |
CN111426950B (zh) | 多尺度时空卷积深度信念网络的风力发电机故障诊断方法 | |
CN105354587A (zh) | 一种风力发电机组齿轮箱的故障诊断方法 | |
CN104200032B (zh) | 广义负荷建模中基于时段性的横向时间轴聚类方法 | |
CN103488869A (zh) | 一种最小二乘支持向量机的风力发电短期负荷预测方法 | |
CN110689069A (zh) | 一种基于半监督bp网络的变压器故障类型诊断方法 | |
CN104408562A (zh) | 一种基于bp神经网络的光伏***发电效率综合评估方法 | |
CN110879377B (zh) | 基于深度信念网络的计量装置故障溯源方法 | |
CN110705879A (zh) | 一种高比例可再生能源接入下的电网脆弱性评估方法 | |
CN108876163A (zh) | 综合因果分析与机器学习的暂态功角稳定快速评估方法 | |
CN111091141B (zh) | 一种基于分层Softmax的光伏背板故障诊断方法 | |
CN104218571B (zh) | 一种风力发电设备的运行状态评估方法 | |
CN110705887A (zh) | 一种基于神经网络模型的低压台区运行状态综合评价方法 | |
CN108667069A (zh) | 一种基于偏最小二乘法回归的短期风电功率预测方法 | |
CN106611243A (zh) | 一种基于garch模型的风速预测残差修正方法 | |
CN110889565B (zh) | 一种基于多维矩阵决策的配网巡检周期计算方法 | |
CN116842337A (zh) | 基于LightGBM优选特征与COA-CNN模型的变压器故障诊断方法 | |
CN113379116A (zh) | 基于聚类和卷积神经网络的台区线损预测方法 | |
Yan et al. | Deep learning based total transfer capability calculation model | |
CN112085108A (zh) | 基于自动编码器及k均值聚类的光伏电站故障诊断算法 | |
CN116317937A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170531 |
|
RJ01 | Rejection of invention patent application after publication |