文 | 威猛贝塔编辑 | 威猛贝塔●○前言○●由于风电机组受到变速变载的影响，对风电机组振动信号直接进行频谱或包络分析会受到频率和幅值调制的影响而无法进行有效的故障诊断，因此需要提取振动信号的瞬时转速对振动信号进行阶次跟踪以进行阶次分析。
为降低安置转速传感器的安装成本以及数据传输成本，所以提岀了一种基于迭代频域能量算子的风电机组无键相瞬时转速信息提取方法，弥补了常规时频分析提取时频脊线方法无法对感兴趣的振动成分，进行有效的阶次跟踪的不足。
通过风电场实际采集的振动数据验证了该方法的有效性，并且与窄带解调以及多阶次概率加权法进行了对比，验证了所提岀方法的有效性。
●○风电机组振动信号○●随着风电机组设备趋向于大型化、海洋化，风电机组运行状态监测成为了保障风电机组安全可靠运行的重要支撑技术，风电机组传动系统一般都配置了振动监测系统，对主轴承、增速齿轮箱、发电机轴承等机械部件进行状态监测。
由于风电机组的变风速、变载荷工况特点，使得振动信号包含幅值及相位调制现象，影响了状态监测的效果，所以强调阶次跟踪是解决转速波动问题最为直接有效的方法，由于阶次跟踪的实施需要与振动信号进行同步采集转速信息。
风电机组现有的转速信号由包含发电机编码器的SCADA系统，进行采集铁由于SCADA与CMS系统存在时间差和采样频率差异，通过发电机编码器所提供的转速信号对大转速波动情况下，振动信号进行阶次跟踪将受到极大影响。
为提高振动信号分析效果，通过振动信号本身，提取风电机组的瞬时转速成为了目前的研究重点，叫根据振动信号提取瞬时转速进行阶次跟踪的优点还在于，瞬时转速信息的获取无需额外安置转速编码器。
在降低了硬件成本及数据传输、存储成本的同时，保证了通过CMS系统采集的振动信号进行状态监测的可靠性，通过振动信号本身进行瞬时转速提取主要有两种方法，分别为提取时频脊线和通过对相位解调获取的相位进行微分冋。
其中对于提出用同步提取转换的方法来提高时频分析结果的时频分辨率，从而提升瞬时频率估计的准确性将Wigner-Ville分布与小波尺度谱进行融合，之后用峰值搜索提取时频脊线，有效地提取出振动信号的瞬时频率。
以上方法均选取一个感兴趣的振动成分并提取其时频脊线来进行瞬时频率估计，而部分场合由于噪声以及传递函数的影响，选取单一的时频脊线估计瞬时频率容易出现较大的误差。
所以专家提出了一种概率融合，提取短时变换中多个振动成分的时频脊线，并进行多阶次加权平均来估计瞬时频率的多阶次概率加权法。
通过Hilbert变换窄带解调估计出风电机组10s瞬时转速，并进行后续角度域同步平均分析，同时还提出了一种通过Hilbert变换进行迭代相位解调的瞬时频率估计方法。
相较于未迭代的窄带解调瞬时频率估计方法，该方法能够更有效地提取感兴趣的振动成分的瞬时频率并通过其进行阶次跟踪，提高了故障诊断的准确度。
风力涡轮机齿轮箱在不同运行条件下的状态监测，通常齿轮箱系统是非线性的，目前已开发了简化的非线性齿轮模型，在此基础上首次应用提出的时频分析方法。
在模拟中及实际风力，涡轮机齿轮箱实验数据中检验了不同负载的影响，为对风电机组传动系统振动信号的瞬时频率进行精确的估计，本实验提出了一种基于迭代频域能量算子方法来处理大波动变转速风电机组振动信号，并通过实测数据对风电机组进行验证。
验证结果表明，尽管风电机组实测数据的振动成分复杂，但通过IFDEO进行迭代更新瞬时频率并进行阶次跟踪能够将高阶振动成分与其边频带进行有效地分离，在此基础上进行进一步的信号处理能够有效地识别风电机组传动系统的故障。
●○迭代频域能量算子○●由于频率能量算子进行相位解调获取信号的瞬时频率需要单分量信号，因此为降低频率混叠的影响，选取的振动成分一般应为低阶频率成分，并对其进行窄带解调，假设一阶振动成分的频率为0且转速波动为±W%，则其在频谱上的带宽为1±W%，因此其高阶谐波的带宽为1±W%。
如选取一阶振动成分进行瞬时频率估计，转速波动应限制在±33%以内，以保证一、二阶振动成分在频域内能够有效地分离，选取更高阶的振动成分进行瞬时频率估计将导致转速波动更为有限，如选取二阶振动成分进行瞬时频率估计的转速波动将限制在±20%以内。
而选取三阶振动成分进行瞬时频率估计的转速波动将限制在±14.3%以内，尽管通过提取低阶振动成分进行瞬时频率估计，并进行阶次跟踪能够降低绝大部分转速波动造成的频率混叠的影响，使得绝大部分振动成分与其倍频在阶次谱上不再混叠。
由于低阶振动成分转速波动缓慢，通过低阶振动成分进行瞬时频率估计仍会产生较大的误差，使得高阶振动成分极易与其边频带混叠在一起，降低了故障诊断的准确度。
因此，可以通过选取阶次跟踪后的角度域振动信号中更高阶的振动成分进行进一步的瞬时频率估计，并进行迭代阶次跟踪，通过提高瞬时频率精确度以实现细化阶次谱的目的。
本实验提出的通过IFDEO进行振动信号无键相阶次跟踪的流程图所示，在进行阶次跟踪后的角度域振动信号可进行后续的包络分析或更进一步的信号处理以实现故障诊断的目的。
通过便携式多通道风电数据采集系统以12.8kHz在某风场采集了2台1.5MW机组齿轮箱的加速度信号，齿轮箱传动比为94.995，加速度传感器安置于低速轴后轴承，灵敏度为100mV/g，同时在高速轴安置电涡流传感器采集高速轴瞬时速度信息。
●○作者观点○●每台机组采集信号时长为200s，信号采集期间，机组均处于正常并网发电状态，并未因数据采集做出任何特殊控制处理，后续维修发现1号机组输入轴出现故障，图2、3分别为2台机组的时域、2000Hz以内频谱图和时频分析，由于两台机组在200s内转速波动大，所以我们可以看出图2，3频谱图出现了频率混叠的情况，导致无法进行后续的信号处理与分析。
图4为1号机组振动信号通过IFDEO估计瞬时转速并进行阶次跟踪后的阶次谱图，由于高速轴及其谐波振动振幅明显，因此以高速轴转速为参考阶次，依次提取高速轴0.5阶，19阶高速轴齿轮啮合成分以及38阶振动成分2阶高速轴齿轮啮合成分进行阶次跟踪。
由于整段信号的振动成分复杂，转速波动剧烈约为士13%，且多个振动成分阶次接近，因此在第一次进行提取0.5阶振动成分进行转速信息估计的过程中将整段信号等时间分成10段，每段有5s的重叠以进行后续的拼接，每段信号的转速波动平均约为±3%。
在此基础上，估计每段信号的0.5阶振动成分，最终通过汉宁窗将每段信号的瞬时频率拼接得到整段信号的瞬时频率估计。
随着迭代次数以及选取振动成分阶次的增加，阶次谱振动成分的模糊程度明显降低，随着瞬时频率精确度的提升，齿轮啮合频率的边频带与啮合频率分离，图为2号机组振动信号通过IFDEO进行瞬时转速估计并通过阶次跟踪得到的阶次谱图，相较于2号机组，1号机组的阶次谱中出现了更多的频率成分。
通过后续对2台机组的角度域振动信号提取包络得到各自的平方包络谱，其阶次成分为0.0053以及其倍频，因此可推算出为风电传动系统输入轴松动，出现1/2X及其倍频，第0.5阶振动成分进行窄带解调获取的瞬时转速进行阶次跟踪后的2台机组的阶次谱图。
由于振动信号的信噪比低和振动成分复杂的影响，相较于IFDEO的瞬时频率，估计岀的瞬时频率误差更大，导致阶次谱中高阶振动成分频率与其边频带混淆，影响后续故障诊断结果。
通过MOPA提取时频脊线后的2台机组的阶次谱图，图中选取了第1阶、2.5阶、18阶、36阶以及78阶时频脊线进行瞬时频率估计，虽然其准确度相较于提取单一的时频脊线有所提升，但因其仍具有信噪比低、时频分辨率低的特点，导致通过MOPA估计的瞬时频率误差较大，阶次谱图模糊，影响后续故障诊断的有效实行。
图分别为2台机组通过电涡流传感器采集的瞬时速度信号进行阶次跟踪的阶次谱，通过对比振动信号可以看岀，本文提岀的方法所得到的阶次谱，能够将高阶振动频率分开阶次谱频带60耀100阶，并与图中的阶次谱最为相近，因此IFDEO能够保证后续解调获得的平方包络谱有效识别1号机组故障。
●○作者观点○●本文提岀了基于IFDEO的无键相瞬时频率估计方法，通过实际采集2台风电机组信号进行验证，确定该方法能够在变转速的情况下，对风电机组传动系统输入轴松动进行有效识别，利用能量算子解调的方法提取瞬时频率的时频分辨率比通过提取时频脊线进行瞬时频率估计的时频分辨率高。
与相位解调法相比，通过迭代提取更高阶的振动成分瞬时频率能够进一步提咼瞬时频率的精确度，使得高阶振动成分在阶次谱中能够和临近振动成分有效分离，该方法能够更有效地提取感兴趣的振动成分的瞬时频率，并进行阶次跟踪。
该方法不受硬件的限制，无需安置额外的传感器获取键相信号，因此能够降低风电机组的监测成本，同时能够对长时间尧大转速波动情况下采集的振动信号进行瞬时频率估计，作为需要长时间采集的信号的前期处理。
