CN112302976B - 一种基于熵权法的风机叶片故障预警方法 - Google Patents
一种基于熵权法的风机叶片故障预警方法 Download PDFInfo
- Publication number
- CN112302976B CN112302976B CN202010987002.XA CN202010987002A CN112302976B CN 112302976 B CN112302976 B CN 112302976B CN 202010987002 A CN202010987002 A CN 202010987002A CN 112302976 B CN112302976 B CN 112302976B
- Authority
- CN
- China
- Prior art keywords
- fan
- signal
- early warning
- characteristic
- fan blade
- 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
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D27/00—Control, e.g. regulation, of pumps, pumping installations or pumping systems specially adapted for elastic fluids
- F04D27/001—Testing thereof; Determination or simulation of flow characteristics; Stall or surge detection, e.g. condition monitoring
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D25/00—Pumping installations or systems
- F04D25/02—Units comprising pumps and their driving means
- F04D25/08—Units comprising pumps and their driving means the working fluid being air, e.g. for ventilation
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/26—Rotors specially for elastic fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
本发明涉及风机故障诊断技术领域，公开了一种基于熵权法的风机叶片故障预警方法，包括：建立一数据存储单元；提取待测信号X和参考信号Y，并分别对所述待测信号X和所述参考信号Y进行待测信号特征矩阵A和参考信号特征矩阵B的建立；基于熵权法分别计算A和B中同一位置特征向量的数值；构建预警指标a及门槛值矩阵alfa，并根据预警门槛值aij和预警指标a来判断风机叶片是否正常并建立状态表Tab。本发明技术方案的风机叶片故障预警方法规避了采用以神经网络等建立整个风场故障识别模型过程中建立完整的数据集的困难度，提高了风机叶片故障预警的精确度和稳定性。
Description
技术领域
本发明涉及风机故障诊断技术领域，具体而言，涉及一种基于熵权法的风机叶片故障预警方法。
背景技术
风机长期运行在露天的环境中，叶片由于较恶劣的天气环境常常发生腐蚀、开裂等故障，如果不能够及时的将故障排除，不但影响风机的发电效率，长期的运行会造成叶片的故障程度加深，甚至影响叶片的使用寿命，而随着人工智能在工业各方面的渗透，针对叶片的故障诊断涌现出多种智能算法，主要以神经网络等建立整个风场的故障识别模型为主，但是这种算法需要叶片的信号声音库足够完善，故障类型足够多，考虑到实际情况，完整数据集的建立过程困难较大，且这类方法多处于实验室阶段，实际应用案例极少。
发明内容
本发明目的在于提出一种基于熵权法的风机叶片故障预警方法，考虑到一系列实际问题，采用风机叶片正常时的信号，建立不同风速、不同风机的叶片正常时的最大可预警指标的门槛值矩阵。在对风机进行诊断预警时，计算该风机的预警指标值，并通过对比分析该风机的相同风速下的门槛值，来对风机作出是否预警的指示。进一步而言，利用小波变换从风机叶片正常的声音信号中提取能量特征矩阵，并基于熵权法建立故障预警指标的门槛值矩阵，将计算得到的待测风机的预警指标与相对应的门槛值进行比较分析，并改变小波分解的层数重复上述步骤，最后综合分析实现待测风机是否进行故障预警的任务。
本发明提供的技术方案，核心构思是通过待测信号与参考信号计算预警指标，并构建预警指标的门槛值矩阵，最后通过对比分析预警指标及其门槛值的大小综合判断风机叶片是否故障。在本发明的多个具体实施例中所述的方法包括，
步骤1.建立一数据存储单元，所述存储单元定期存储和更新一风机组在无故障情况下提取的所有风机叶片声音信号；所述声音信号包括所述风机组所有的风机处于不同风速情况下发出的声音；提取所述声音信号作为所述风机组中所有风机叶片故障预警时的备用参考信号。
步骤2.从所述风机组某风机中提取待测信号X和参考信号Y，并分别对所述待测信号X和所述参考信号Y进行小波分解并提取各自能量系数特征属性值，各个所述能量系数特征属性值组成能量系数特征属性特征维，进而分别构建待测信号特征矩阵A和参考信号特征矩阵B。
步骤3.分别提取所述待测信号特征矩阵A和所述参考信号特征矩阵B的第n列同一位置的所述能量系数特征属性值，组成m×2的特征矩阵Q且共计K个，其中n＝{1，2，3，……K}，K为所述小波分解的层数，m为所述信号划分的段数；利用熵权法分别计算K个所述特征矩阵Q中第一维和第二维所述能量系数特征维的特征权重值，分别为bn1和bn2。
步骤4.构建所述风机组风机叶片故障预警的指标a，用以定量描述叶片的运行状态，公式如下：
并由所述风机叶片故障预警指标a构建门槛值矩阵alfa。
步骤5.根据所述待测风机的编号和风速，遍历所述门槛值矩阵alfa，确定所述待测风机的预警门槛值aii和所述风机组风机叶片故障预警的指标a，并建立所述风机组中风机叶片状态信息表Tab。
步骤6.调整所述小波分解层数K值大小及其变化次数，并将每次计算结果记录在所述风机组中风机叶片状态信息表Tab中，判断所述风机组中风机的运行状况。
进一步的，所述步骤2中所述待测信号X和所述参考信号Y设置一定数目的采样点作为固定帧长，定义为wlen，则待测信号X和所述参考信号Y分别被分成m段信号，即
其中m为式中结果的整数部分，wlen的采样点个数要求至少包含所述风机旋转一周的声音采样点，X、Y代表所述待测信号和所述参考信号，wlen满足
进一步的，所述步骤2中所述小波分解，利用小波能量公式计算所述参考信号Y和所述待测信号X分解后的能量系数，即En＝∑|xn|2其中xn为被分解后的每段信号对应参考值，n＝{1，2，3......K}，K为所述信号小波分解的层数；所述待测信号X和参考信号Y分别被分成m段信号，每段信号分别被进行K层小波分解，每段每一层小波分解后的所述信号能量系数为能量系数特征属性值，每段K层小波分解后的所述信号能量系数特征属性组成能量系数特征维，所述m段能量系数特征维组成待测信号特征矩阵A和参考信号特征矩阵B。
进一步的，所述步骤3分别计算K个所述特征矩阵Q中第一维和第二维所述特征属性的特征权重值采用的熵权法，即首先根据所述特征矩阵Q中各个所述特征属性值进行归一化即
其中i＝{1，2，3，......m}，j＝{1，2}，min(Xi)为第i段所述声音信号对应的能量系数特征属性值的最小值，max(Xi)为第i段所述声音信号对应的能量系数特征属性值的最大值，Xij为第i段第j维度所述声音信号对应的能量系数特征属性值，以获得所述特征矩阵Q的归一化矩阵Q′。
进一步的，所述步骤4中所述构建叶片故障预警指标a即
其中bn1和bn2分别表示所述特征矩阵Q中第一维和第二维所述能量系数特征属性的特征权重值。
进一步的，所述步骤4所述的门槛值矩阵alfa，分别从三个时间段的声音信号存储单元中提取同风速、同一个风机的声音信号，且相互之间所述风机叶片故障预警指标a按照步骤2到4进行计算，分别用a1、a2、a3表示，则所述风机在此风速下的故障预警指标a为上述三个值中的最大值；针对所述风机组中其他风机做同样的处理即得到预警指标a的门槛值矩阵alfa：
其中u是所述风机组中风机编号总数，v是所述风机组中风机不同风速的数目，i的取值范围是[1，u]间的整数，j的取值范围是[1，v]的整数。
进一步的，所述步骤5中所述待测风机的预警门槛值aij和所述风机叶片的故障预警的指标a二者对比，当a＜aij时，则所述待测风机的叶片正常，反之说明所述待测风机叶片预警指标超过所述风机叶片预警门槛值并建立与之对应的所述风机叶片状态信息表Tab，记录诊断结果。
进一步的，根据所述归一化矩阵Q′中的能量系数通过计算公式
其中
如果pij＝0则定义
得出所述归一化矩阵Q′中的特征权重值。
进一步的，所述步骤2中所述待测信号X为采集所述风机组中某风机在一固定时刻风机叶片发出的一段声音信号；所述参考信号Y为遍历所述数据存储单元中的所述备用参考信号，提取与所述待测信号X对应的风机在同一风速下的声音信号，且所述待测信号X和所述参考信号Y两种信号的保持时间和帧长一致。
进一步的，所述步骤2针对所述待测信号X和所述参考信号Y两种信号中同一段声音信号分别进行K层小波分解，K是人为设定的一个常数，并利用小波能量公式提取分解后各自的能量系数，各自生成对应的待测信号特征维X1和参考信号特征维Y1；剩余段所述声音信号做上述同样的处理，各自生成待测信号特征矩阵A和参考信号特征矩阵B且大小均为m×K，其中m为所述特征矩阵行数，即声音的分段数，K为所述特征矩阵的列数，即小波分解的层数。
本发明技术方案带来的技术效果的一个方面在于，规避了采用以神经网络等建立整个风场故障识别模型过程中建立完整的数据集的困难度，提高了风机叶片故障预警的精确度和稳定性，相比于建立故障识别模型方法而言，操作简便，计算量较小，不要求完备的故障叶片数据库，更贴合于实际工程中的应用。
附图说明
图1为本发明一种基于熵权法的风机叶片故障预警方法的流程图；
图2为本发明一个实施例中风机K层小波分解波形分析示意图。
具体实施方式
实施例1
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
具体实现方法如下所述：
如图1所示，为本发明一种基于熵权法的风机叶片故障预警方法的流程图，其步骤包括：
步骤1.针对一个风场，建立风机叶片声音信号的数据存储单元，采集的声音包括风场所有风机不同风速尺度下的声音信号，将这些声音信号作为风机叶片故障预警时的备用参考信号，这个存储单元需要每隔一定的时间进行全部更新，时间段以月为单位来计量；
步骤2.在一固定时刻，对一风机的叶片进行故障诊断时，采集叶片在该时刻发出的一段声音信号作为待测信号X，设置一定数目的采样点作为固定帧长，定义为wlen，则待测信号可以被分成m段声音信号，具体的计算公式如下：
m是式(1)中结果的整数部分，wlen的采样点个数要求至少包含风机旋转一周声音采样点个数，低速轴的转速为w，单位是r/min，声音的采样频率是f，单位是Hz，则wlen需满足以下关系式：
在参考信号存储单元中寻找与待测信号同台风机，同样风速下的声音信号，且参考信号与待测信号保持时间、帧长一致；
针对两种信号的每段声音信号进行K层小波分解，K是人为设定的一个常数，并利用小波能量公式提取分解后的能量系数，计算公式如下：
En＝∑|xn|2 (3)
其中xn为被分解后的每段信号对应参考值，n＝{1，2，3......K}，本篇幅中K为所述小波分解的层数，优选的，K∈[3，4，5，6，7]。
于是每段声音信号都可以用能量系数生成的向量来表示，则对于待测信号和参考信号，可以通过上述的变换提取出能够表征两种信号本质的能量特征矩阵A和B，其中两个矩阵大大小都为m×K，即矩阵的行数为声音的分段数，矩阵的列为小波分解的层数；
步骤3.分别提取所述待测信号特征矩阵A和所述参考信号特征矩阵B的第n列同一位置的所述能量系数特征属性值，组成m×2的特征矩阵Q且共计K个，其中n＝{1，2，3，……K}，K为所述小波分解的层数，m为所述信号划分的段数；利用熵权法分别计算K个所述特征矩阵Q中第一维和第二维所述能量系数特征维的特征权重值，分别为bn1和bn2。
所述步骤3分别计算K个矩阵Q中第一维和第二维的特征权重值采用的所述熵权法，即为根据所述特征矩阵Q中各个所述能量系数值进行归一化，即
其中i＝{1，2，3，......m}，j＝{1，2}，min(Xi)为第i段所述声音信号对应的能量系数的最小值，max(Xi)为第i段所述声音信号对应的能量系数的最大值，Xij为第i段第j维度所述声音信号对应的能量系数，进行标准化后可以得到数据标准化表，以获得所述特征矩阵Q的归一化矩阵Q′。根据所述归一化矩阵Q′中的能量系数通过计算公式
其中
如果pij＝0则定义
得出所述归一化矩阵Q′中的特征权重值。
步骤4.构建所述风机组风机叶片故障预警的指标a，用以定量描述叶片的运行状态，公式如下：
并由所述风机叶片故障预警指标a构建门槛值矩阵alfa。
步骤5.利用声音信号的存储单元，构建预警指标a的门槛值矩阵alfa，
建立三个时间段上的声音信号的存储单元，分别从三个存储单元中取出同风速下，同一个风机的声音信号，相互之间的a指标按照步骤2到4进行计算，针对这个风机固定工况下可以计算出三个指标的值，分别用a1、a2、a3表示，则该个风机在这一风速下的故障预警指标为取三个指标中的最大值。针对不同的风机做同样的处理即可以得到预警指标a的门槛值矩阵alfa：
式(5)中，u是某一风场的风机的编号总数，v是风场不同的风速尺度数目，i的取值范围是[1，u]间的整数，j的取值范围是[1，v]的整数，
步骤6.根据待测风机的编号和风速情况，从alfa中确定该个风机的预警门槛值，当计算出来的a＜aij时，则叶片正常，否则就预警，提醒风机设备管理人风机叶片可能存在故障，并将本次得到的结果记录在叶片状态信息表Tab中；
步骤7.改变小波分解层数K值，则步骤2到步骤6中涉及到的小波分解层数必须保持一致，通过设置K值的变化次数，即可以实现整个流程的循环的次数，将每次结果都记录在叶片状态信息表Tab中，如果表中的预警次数大于正常次数，则该个风机预警叶片故障，否则显示叶片正常，从而实现了风机叶片的故障预警。
| 小波分解层数 | 1号风机 | 2号风机 | 3号风机 | …… | u号风机 |
| 1 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 |
| 2 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 |
| 3 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 |
| 4 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 |
| 5 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 |
| 6 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 |
| 7 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 |
| …… | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 |
| K | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 | 正常/故障 |
实施例2
本实施例结合现场一实例进一步阐述一种基于熵权法的风机叶片故障预警方法的实施过程及注意事项。
某一风机现场共计5台风机组，现需要对上述5台风机组中所有的风机利用熵权法建立风机叶片故障预警系统，用于现场风机设备的管理。其具体方法如下：
1.设在某一风速下提取待测声音信号Y和参考声音信号X且将待测声音信号Y和参考声音信号X分成4段进行3层小波分解，故测得参考声音信号Y和待测声音信号X如下所示：
所述待测参考声音信号
所述参考声音信号
2.分别提取所述待测信号特征矩阵A和所述参考信号特征矩阵B的第n列所述能量系数并组合成m×2的矩阵Q，其中n＝{1,2,3，……K}，共计K个矩阵Q，此处m为4，K为3，故所述特征矩阵Q分别为：
对Q1、Q2、Q3利用熵权法分别计算第一维和第二维的特征权重，分别为bn1、bn2，其中n＝{1,2,3,4}。首先根据归一化公式
得到各自对应的归一化Q11、Q21、Q31、Q41，即：
根据计算公式
其中
如果pij＝0则定义
得出所述归一化矩阵Q′中的特征权重值，即bn1、bn2。
构建叶片故障预警的指标a，用以定量描述叶片的运行状态，计算公式如下所示：
a＝|b11-b12|+|b21-b22|+|b31-b32|＝0.012+0.278+0.16＝0.45，风机叶片故障预警指标为0.45。
3.构建预警指标a的门槛值矩阵alfa。建立三个时间段上的声音信号的存储单元并采取同风速下、同一个风机的声音信号，预警指标a按照上述步骤进行计算，针对该风机固定工况下可以计算出三个指标的值，分别用a1、a2、a3表示，则该个风机在这一风速下的故障预警指标为取三个指标中的最大值。具体过程如下，假设该风机为1号风机，分别在时间段1、时间段2和时间段3等三个时间段，在风机无故障情况下和实际运行状态下各提取了一段声音信号，且此段声音信号均被分成4段进行3层小波分解，进行一系列运算得出相对应的预警指标a1、a2、a3。测得参考声音信号Y和待测声音信号X如下所示：
上述各组声音信号对应的Q矩阵即分别提取参考声音信号X和待测声音信号Y同一列能量系数，组合成Q矩阵，如下：
对上述各个Q矩阵进行熵权法归一化处理，得出各自的归一化矩阵，具体如下：
根据计算公式
其中
如果pij＝0则定义
得出所述归一化矩阵Q′中的特征权重值，即bi1、bi2。具体如下表所述：
故障预警的指标a，用以定量描述叶片的运行状态，计算公式如下所示：
其中bn1和bn2分别表示所述特征矩阵Q中第一维和第二维所述能量系数特征属性的特征权重值，故a1、a2、a3分别为：0.12、0.66、0.5。该1号风机在这一风速下故障预警指标为取三个指标中的最大值，即为0.66。
其他风机也做上述运算，由于实例上述已设定对现场5台风机利用熵权法建立风机叶片故障预警系统，用于现场风机设备的管理。设在某一风速下提取参考声音信号Y和待测声音信号X且将参考声音信号Y和待测声音信号X分成4段进行3层小波分解，且对风机在三个速度段进行取样，故v为3，u为5，具体如下：
| 风机组 | 风速1 | 风速2 | 风速3 |
| 1号风机 | 0.66 | 0.45 | 0.42 |
| 2号风机 | 0.64 | 0.42 | 0.27 |
| 3号风机 | 0.78 | 0.38 | 0.56 |
| 4号风机 | 0.45 | 0.31 | 0.69 |
| 5号风机 | 0.38 | 0.39 | 0.36 |
根据本实施例第2步计算所述风机组风机叶片的故障预警指标0.45和上表中风机组各个风机在不同风速下的故障预警实际值做比较，即根据待测风机的编号和风速情况，从alfa中确定该个风机的预警门槛值，当计算出来的a<aij时，则叶片正常，否则就预警，提醒风机设备管理人员风机叶片可能存在故障，并将本次得到的结果记录在叶片状态信息表Tab中，具体如下：
最后尝试调整K值大小和次数，即小波分解的层数，将每次结果都记录在叶片状态信息表Tab中，如下表所示，用以风机组叶片故障的预警。
| 小波分解层数 | 1号风机 | 2号风机 | 3号风机 | 4号风机 | 5号风机 |
| 1 | NA/故障 | NA/故障 | 正常/NA | NA/故障 | NA/故障 |
| 2 | NA/故障 | NA/故障 | 正常/NA | NA/故障 | NA/故障 |
| 3 | NA/故障 | NA/故障 | 正常/NA | NA/故障 | NA/故障 |
实施例3
本实施例和实施例2的区别在于利用熵权法对归一化矩阵Q′中的特征属性权重进行二次熵权，比较下最终结果和实施例1之间的差异和联系。下表为实施例2求得的归一化矩阵Q′中的特征权重：
根据所述风机样本特征属性各个特征维的权重计算公式
可以得到所述特征矩阵Y中各个特征维的权重如下表所示：
此时对于故障预警的指标a利用公式
计算所得a1、a2、a3分别为：0.05、0.08、0.03，此时1号风机在该风速下故障预警指标为取三个指标中的最大值，即为0.08。对于风机同样做上述运算处理，，由于实例上述已设定对现场5台风机利用熵权法建立风机叶片故障预警系统，用于现场风机设备的管理。设在某一风速下提取参考声音信号Y和待测声音信号X且将参考声音信号Y和待测声音信号X分成4段进行3层小波分解，且对风机在三个速度段进行取样，故v为3，u为5，具体如下：
具体数据以下列表格体现，如下：
| 风机组 | 风速1 | 风速2 | 风速3 |
| 1号风机 | 0.35 | 0.02 | 0.06 |
| 2号风机 | 0.78 | 0.05 | 0.03 |
| 3号风机 | 0.09 | 0.02 | 0.08 |
| 4号风机 | 0.02 | 0.03 | 0.35 |
| 5号风机 | 0.01 | 0.06 | 0.08 |
根据本实施例所述风机组风机叶片的故障预警指标0.08，参照表中风机组各个风机在不同风速下的故障预警预警实际值做比较，即根据待测风机的编号和风速情况，从alfa中确定该个风机的预警门槛值，当计算出来的a<aij时，则叶片正常，否则就预警，提醒风机设备管理人员风机叶片可能存在故障，并将本次得到的结果记录在叶片状态信息表Tab中，具体如下：
最后尝试调整K值大小和次数，即小波分解的层数，将每次结果都记录在叶片状态信息表Tab中，如下表所示，用以风机组叶片故障的预警。
| 小波分解层数 | 1号风机 | 2号风机 | 3号风机 | 4号风机 | 5号风机 |
| 1 | NA/故障 | NA/故障 | 正常/NA | NA/故障 | NA/故障 |
| 2 | NA/故障 | NA/故障 | 正常/NA | NA/故障 | NA/故障 |
| 3 | NA/故障 | NA/故障 | 正常/NA | NA/故障 | NA/故障 |
通过本实施例的计算结果可知，对基于熵权法提取所述风机组所有风机叶片能量系数矩阵特征维的二次权重相对于实施例2结果一致，即采用熵权法求取所述归一化矩阵特征属性权重或者对归一化矩阵特征维的二次权重均可以对风机组叶片故障进行预警判断。
实施例4
本实施例对某一风机组中所有风机进行的K层小波分解得到一系列离散平稳的小波进行波形分析,此实施例中会涉及Matlab软件的使用。
在对离散型波形分析所用到的函数有swt小波分解和iswt小波重构。假设某一风机场风机组共计5台风机，对所有风机提取组件所述风机特征矩阵D，并利用熵权法求取所述风机的特征属性特征权重值和一系列风机叶片故障预警值进行3层小波分解，对得出的波形进一步分析，和上述实施例2和实施例3中所述风机叶片预警值Tab表做横向比较。具体结果如下：
首先使用命令行实现一维分析，即
1.加载信号(所述风机组某一风机某段声音信号)load noisdopp；
2.设置变量s＝noisdopp；对于SWT变换，如果需要在第k层分解信号,k∈[1,2,……K]，那么原始信号需要能够平分成2^k份。所以如果原始信号的长度不满足要求，需要使用Signal Extension GUI工具或使用wextend函数来扩展它。
3.执行一层平稳小波分解[swa,swd]＝swt(s,1,'db1')；函数执行将产生1层近似和细节的系数，二者信号的长度相等，这也是平稳小波和普通小波两者不同的地方。
4.显示近似、细节小波系数
Subplot(1,2,1),plot(swa)；title('Approximation cfs')
subplot(1,2,2),plot(swd)；title('Detail cfs')
5.使用平稳小波逆变换恢复原始信号
A0＝iswt(swa,swd,'db1')；
重构的误差为err＝norm(s-A0),err＝2.1450e-014
6.从系数中构建近似和细节，输入如下代码
nulcfs＝zeros(size(swa))；
A1＝iswt(swa,nulcfs,'db1')；
D1＝iswt(nulcfs,swd,'db1')；
显示结果如下
Subplot(1,2,1),plot(A1)；title('Approximation A1')；
Subplot(1,2,2),plot(D1)；title('Detail D1')；
7.执行多层平稳小波分解
使用db1小波做3层小波分解[swa,swd]＝swt(s,3,'db1')；产生的近似系数在swa中，细节系数在swd中，而且它们有相同的长度。
8.显示近似和细节的系数
kp＝0；
for i＝1:3
subplot(3,2,kp+1),plot(swa(i,:))；
title(['Approx.cfs level',num2str(i)])
subplot(3,2,kp+2),plot(swd(i,:))；
title(['Detail cfs level',num2str(i)])
kp＝kp+2；
end
9.从系数中重建第3层的近似
mzero＝zeros(size(swd))；
A＝mzero；
A(3,:)＝iswt(swa,mzero,'db1')；
10.从系数中重建细节
D＝mzero；
for i＝1:3
swcfs＝mzero；
swcfs(i,:)＝swd(i,:)；
D(i,:)＝iswt(mzero,swcfs,'db1')；
End
11.从第3层近似和第2、3层细节中重建第1、2层近似
重建2、3层的近似
A(2,:)＝A(3,:)+D(3,:)
A(1,:)＝A(2,:)+D(2,:)；
显示第1、2、3层的近似和细节
kp＝0；
for i＝1:3
subplot(3,2,kp+1),plot(A(i,:))；
title(['Approx.level',num2str(i)]；
subplot(3,2,kp+2),plot(D(i,:))；
title(['Detail level',num2str(i)])
kp＝kp+2；
end
执行后，得到波形如图2所示，根据1、2、3层波形的分析，提取某个风机在某一风速下的趋势，其离散频次和离散幅度代表着再次瞬间风机叶片的稳定状态，映射于上述实施例2和实施例3中所述风机叶片预警值Tab表记录的相关信息，最终结果毫无二致。
实施例5
本实施例将对所述风机组中所有风机提取的所述待测信号特征矩阵A和所述参考信号特征矩阵B的第n列所述能量系数并组合成m×2的矩阵Q，第一特征维和第二特征维对应的参考声音信号Y或者待测声音信号X的先后顺序是否影响最终判断结果做进一步说明。同样假设某一风机场某风机组所有风机共计5台，对5台风机所有的风机叶片提取的特征矩阵做3层小波分解处理，其结果如下：
所述待测声音信号
所述参考声音信号
分别提取所述待测信号特征矩阵和所述参考信号特征矩阵的第n列所述能量系数并组合成m×2的矩阵Q，其中n＝{1,2,3，……K}，共计K个矩阵Q，此处m为4，K为3，故所述特征矩阵Q分别为：
此时待测信号特征矩阵特征维为特征矩阵Q的第一特征维，参考信号特征矩阵特征维为特征矩阵Q的第二特征维。而实际在求取故障预警的指标a，用以定量描述叶片的运行状态，计算公式如下所示：
其中bn1和bn2分别表示所述特征矩阵Q中第一维和第二维所述能量系数特征属性的特征权重值,可知特征矩阵Q的第一特征维和第二特征维分别对应于参考声音信号Y或者待测声音信号X的先后顺序对风机叶片故障预警指标a无影响。
Claims (9)
1.一种基于熵权法的风机叶片故障预警方法，其特征在于，包括:
步骤1.建立一数据存储单元，所述存储单元定期存储和更新一风机组在无故障情况下提取的所有风机叶片声音信号；所述声音信号包括所述风机组所有的风机处于不同风速情况下发出的声音；提取所述声音信号作为所述风机组中所有风机叶片故障预警时的备用参考信号；
步骤2.从所述风机组某风机中提取待测信号X和参考信号Y，并分别对所述待测信号X和所述参考信号Y进行小波分解并提取各自能量系数特征属性值，各个所述能量系数特征属性值组成能量系数特征属性特征维，进而分别构建待测信号特征矩阵A和参考信号特征矩阵B；
步骤3.分别提取所述待测信号特征矩阵A和所述参考信号特征矩阵B的第n列同一位置的所述能量系数特征属性值，组成m×2的特征矩阵Q且共计K个，其中n＝{1,2,3，……K}，K为所述小波分解的层数，m为所述信号划分的段数；利用熵权法分别计算K个所述特征矩阵Q中第一维和第二维的能量系数特征维的特征权重值，分别为bn1和bn2；
步骤4.构建所述风机组风机叶片故障预警的指标a，用以定量描述叶片的运行状态，公式如下：
并由所述风机叶片故障预警指标a构建门槛值矩阵alfa；
步骤5.根据待测风机的编号和风速，遍历所述门槛值矩阵alfa，确定待测风机的预警门槛值aij和所述风机组风机叶片故障预警的指标a，并建立所述风机组中风机叶片状态信息表Tab；
步骤6.调整所述小波分解层数K值大小及其变化次数，并将每次计算结果记录在所述风机组中风机叶片状态信息表Tab中，判断所述风机组中风机的运行状况。
2.根据权利要求1所述的一种基于熵权法的风机叶片故障预警方法，其特征在于：所述步骤2中所述待测信号X和所述参考信号Y设置一定数目的采样点作为固定帧长，定义为wlen，则待测信号X和所述参考信号Y分别被分成m段信号，即
其中m为式中结果的整数部分，wlen的采样点个数要求至少包含所述风机旋转一周的声音采样点，X、Y代表所述待测信号和所述参考信号，wlen满足
低速轴的转速为w，单位是r/min；声音的采样频率是f，单位是Hz。
3.根据权利要求1所述的一种基于熵权法的风机叶片故障预警方法，其特征在于：所述步骤2中所述小波分解，利用小波能量公式计算所述参考信号Y和所述待测信号X分解后的能量系数，即En＝∑|xn|2，其中xn为被分解后的每段信号对应参考值，n＝{1,2,3……K}，K为所述信号小波分解的层数；所述待测信号X和参考信号Y分别被分成m段信号，每段信号分别被进行K层小波分解，每段每一层小波分解后的信号能量系数为能量系数特征属性值，每段K层小波分解后的信号能量系数特征属性组成能量系数特征维，m段能量系数特征维组成待测信号特征矩阵A和参考信号特征矩阵B。
4.根据权利要求1所述的一种基于熵权法的风机叶片故障预警方法，其特征在于：所述步骤3分别计算K个所述特征矩阵Q中第一维和第二维特征属性的特征权重值采用的熵权法，即首先根据所述特征矩阵Q中各个所述特征属性值进行归一化即
其中i＝{1，2，3，......m}，j＝{1，2}，min(Xi)为第i段所述声音信号对应的能量系数特征属性值的最小值，max(Xi)为第i段所述声音信号对应的能量系数特征属性值的最大值，Xij为第i段第j维度所述声音信号对应的能量系数特征属性值,以获得所述特征矩阵Q的归一化矩阵Q′。
5.根据权利要求1所述的一种基于熵权法的风机叶片故障预警方法，其特征在于：所述步骤4中所述构建叶片故障预警指标a即
其中bn1和bn2分别表示所述特征矩阵Q中第一维和第二维所述能量系数特征属性的特征权重值。
6.根据权利要求1所述的一种基于熵权法的风机叶片故障预警方法，其特征在于：所述步骤4门槛值矩阵alfa，分别从三个时间段的信号存储单元中提取同一风速、同一风机的声音信号，单个时间段的声音信号进行所述风机叶片故障预警指标a的计算，分别用a1、a2、a3表示，则所述该风机在该风速下的故障预警指标a为上述三个值中的最大值；针对所述风机组中其他风机做同样的处理即得到预警指标a的门槛值矩阵alfa：
其中u是所述风机组中风机编号总数，v是所述风机组中风机不同风速的数目，i的取值范围是[1,u]间的整数，j的取值范围是[1,v]的整数。
7.根据权利要求1所述的一种基于熵权法的风机叶片故障预警方法，其特征在于：所述步骤5中所述待测风机的预警门槛值aij和所述风机叶片的故障预警的指标a二者对比，当a<aij时，则所述待测风机的叶片正常，反之说明所述待测风机叶片预警指标超过所述风机叶片预警门槛值并建立与之对应的所述风机叶片状态信息表Tab，记录诊断结果。
8.根据权利要求1所述的一种基于熵权法的风机叶片故障预警方法，其特征在于：所述步骤2中所述待测信号X为采集所述风机组中某风机在一固定时刻风机叶片发出的一段声音信号；所述参考信号Y为遍历所述数据存储单元中的所述备用参考信号，提取与所述待测信号X对应的风机在同一风速下的声音信号，且所述待测信号X和所述参考信号Y两种信号保持时间和帧长一致。
9.根据权利要求1所述的一种基于熵权法的风机叶片故障预警方法，其特征在于：所述步骤2针对所述待测信号X和所述参考信号Y两种信号中同一段声音信号分别进行K层小波分解，K是人为设定的一个常数，并利用小波能量公式提取分解后各自的能量系数，生成对应的待测信号特征维和参考信号特征维；剩余段所述声音信号做上述处理，各自生成所述待测信号特征矩阵A和所述参考信号特征矩阵B且大小均为m×K，其中m为所述特征矩阵行数，即声音的分段数，K为所述特征矩阵的列数，即小波分解的层数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010987002.XA CN112302976B (zh) | 2020-09-18 | 2020-09-18 | 一种基于熵权法的风机叶片故障预警方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010987002.XA CN112302976B (zh) | 2020-09-18 | 2020-09-18 | 一种基于熵权法的风机叶片故障预警方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112302976A CN112302976A (zh) | 2021-02-02 |
| CN112302976B true CN112302976B (zh) | 2022-02-01 |
Family
ID=74483297
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010987002.XA Active CN112302976B (zh) | 2020-09-18 | 2020-09-18 | 一种基于熵权法的风机叶片故障预警方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112302976B (zh) |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104155650B (zh) * | 2014-08-14 | 2016-11-23 | 电子科技大学 | 一种基于熵权值法点迹质量评估的目标跟踪方法 |
| CN107133394A (zh) * | 2017-04-14 | 2017-09-05 | 浙江经贸职业技术学院 | 基于响应面法和熵权法相结合的风扇多目标性能优化方法 |
| CN107346541B (zh) * | 2017-06-02 | 2020-02-18 | 华南理工大学 | 一种基于超声射频时间序列小波分析的组织定征方法 |
| CN108131321B (zh) * | 2018-02-12 | 2019-06-28 | 山东理工大学 | 一种轴流式风机失速故障监测方法 |
| CN110057586B (zh) * | 2019-04-25 | 2020-10-30 | 长江大学 | 轴承故障振动信号Schatten改进小波包与重构降噪方法 |
| CN110905735A (zh) * | 2019-12-17 | 2020-03-24 | 天津市津能风电有限责任公司 | 一种基于声音周期性的叶片故障诊断方法 |
| CN111461564A (zh) * | 2020-04-08 | 2020-07-28 | 湖南大学 | 基于云模型和最优组合赋权的风电机组功率特性评估方法 |
-
2020
- 2020-09-18 CN CN202010987002.XA patent/CN112302976B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112302976A (zh) | 2021-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109597401B (zh) | 一种基于数据驱动的设备故障诊断方法 | |
| CN111562358B (zh) | 基于联合模型的变压器油中气体含量预测方法及系统 | |
| CN105224872B (zh) | 一种基于神经网络聚类的用户异常行为检测方法 | |
| Bagnall et al. | A run length transformation for discriminating between auto regressive time series | |
| CN112417028B (zh) | 一种风速时序特征挖掘方法及短期风电功率预测方法 | |
| DE112015005427B4 (de) | Qualitätssteuermaschine für komplexe physikalische Systeme | |
| CN108399434B (zh) | 基于特征提取的高维时间序列数据的分析预测方法 | |
| CN111027629A (zh) | 基于改进随机森林的配电网故障停电率预测方法及系统 | |
| CN116861354A (zh) | 电力变压器监测数据的异常识别及模式区分方法及系统 | |
| CN115564114A (zh) | 一种基于图神经网络的空域碳排放短期预测方法及系统 | |
| CN111553186A (zh) | 一种基于深度长短时记忆网络的电磁信号识别方法 | |
| CN112434743A (zh) | 基于gil金属微粒局部放电时域波形图像的故障识别方法 | |
| Adeyemo et al. | Effects of normalization techniques on logistic regression in data science | |
| Gowtham Sethupathi et al. | Efficient rainfall prediction and analysis using machine learning techniques | |
| CN112302976B (zh) | 一种基于熵权法的风机叶片故障预警方法 | |
| CN112445690B (zh) | 信息获取方法、装置及电子设备 | |
| CN112651499A (zh) | 一种基于蚁群优化算法和层间信息的结构化模型剪枝方法 | |
| CN115907079B (zh) | 一种基于注意力时空图卷积网络的空域交通流量预测方法 | |
| Sahraoui et al. | New criteria for wrapper feature selection to enhance bearing fault classification | |
| CN111028086A (zh) | 基于聚类与lstm网络的增强指数跟踪方法 | |
| CN113657726B (zh) | 基于随机森林的人员的危险性分析方法 | |
| CN115452362A (zh) | 一种齿轮箱故障诊断方法 | |
| CN115879369A (zh) | 一种基于优化LightGBM算法的磨煤机故障预警方法 | |
| CN115079052A (zh) | 一种变压器故障诊断方法和系统 | |
| CN113296947A (zh) | 基于改进XGBoost模型的资源需求预测方法 |
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 |