一种海上风电机组高速轴承剩余寿命预测方法

本发明属于基于振动监测分析的海上风电机组高速轴承剩余寿命预测领域，具体来讲，提出了一种全新的具有严格单调性的高速轴承健康指标构建方法，可实现仅利用在线历史数据进行高速轴承剩余寿命预测。

背景技术：

1、滚动轴承作为海上风电机组传动链中主要的运动部件之一，其一旦发生故障将会影响设备的生产效率，甚至导致整个机组停机维修，带来不必要的经济损失。海上风电机组传动链高速轴承长期处于高速、高温、变载的工作环境，故障率高，运维保养难度较大，因此对海上风电机组传动链高速轴承进行故障诊断以及寿命预测，以保证设备安全稳定运转是非常有必要的，这也是目前工业化进程的要求和发展趋势。相比于故障诊断，海上风电机组高速轴承剩余寿命预测的研究意义更大，可利用剩余寿命预测结果，分析设备未来的运行状态以及失效时刻，考虑是否提前准备设备部件，可以为海上风电机组预防性维修决策提供依据，从而延长相关机械设备的服役时间，提高生产效率，并且在一定程度上预防重大安全事故的发生等。

2、在滚动轴承剩余寿命预测领域，目前的预测手段主要可分为三类：基于物理失效模型、基于数据驱动模型及两者的融合模型。基于物理失效模型的方法通过建立相应的物理或数学模型来表征机械设备的衰退趋势，需要结合特定的先验知识或缺陷劣化方程，并根据观测数据对模型参数进行优化。基于数据驱动的方法通过采集的运行数据对机械设备的退化趋势进行预测，该方法不需要相应的专家经验知识，因此对于研究复杂机械系统较为实用。融合模型方法就是将物理失效模型与数据驱动模型结合，主要思想是通过实时的监测数据对给定的经验退化模型进行修正，但是由于融合了两种方法，使得构建的模型结构较为复杂且难以求解，因此该方法研究较少。基于模型的方法能够利用相关经验知识，同时可以结合观测数据潜在的故障信息，因此得到的剩余寿命预测结果更具说服力。不过，由于机械设备的复杂程度越来越高，再加上环境的干扰，使得其性能退化过程具有较强的随机性，因此构建合理的健康指标显得尤为重要。同时，现有的很多剩余寿命预测方法仅在试验台数据或训练集上表现良好，且往往都需要先利用全寿命数据进行训练，这在实际的工程应用中存在一定局限性。

3、高速轴承剩余寿命预测的重点之一在于健康指标构建，现有的健康指标构建方法大都存在以下不足：①利用神经网络、支持向量机等方法构建健康指标时，需要用到全寿命的振动数据进行训练，且训练结果往往仅适用于同故障类型的轴承。②现有方法主要依据特征信号不同时刻状态与初始时刻状态之间的差别进行健康指标构建，忽略了轴承的退化是不可逆过程，良好的健康指标应具有严格的单调性趋势。③有的方法能够避免使用全寿命数据进行训练，但构建的指标可预测性较差，难以进行准确的剩余寿命预测。

4、现有技术cn110909509a轴承寿命预测方法具体采用了一些技术手段，包括如下步骤，步骤一：将sdae、infogan和lsgan相结合，构建infolsgan，自动地从噪声数据中提取可解释的鲁棒特征，解决梯度消失问题；步骤二：采用基于ac的训练算法训练infolsgan，减少训练时间，加快收敛速度；步骤三：根据训练后的infolsgan，利用softmax分类器预测测试样本中风电齿轮箱轴承的剩余寿命。这种方法在小样本情况下能预测风电齿轮箱轴承的剩余寿命。但是这种算法过于复杂，且大样本情况下无法处理，仅可应用于实验用途。现有技术《风电机组传动链高速轴承剩余寿命预测》(作者黄乙珂)中公开了轴承的寿命预测方法，其中利用了信号的敏感特征，随后得到失效阈值，将敏感特征与失效阈值进行比较，从而在分析粒子滤波解决预测问题所存在优缺点的基础上，提出利用高斯过程回归方法预测结果作为观测值，同时通过粒子滤波对预测结果进行优化，从而构建出pf-gpr融合预测模型；确定了预测模型状态方程的具体形式及构建过程、核函数类型及剩余寿命计算方式。

5、现有技术的预测方法存在预测不够准确的问题。

技术实现思路

1、本发明通过高速轴承历史振动信号，构建了可预测性良好的健康指标，实现了较为精准的高速轴承剩余寿命预测。

2、本发明提出一种递推型的高速轴承健康指标构建方法，既满足轴承真实退化规律，且具有良好可预测性；同时在寿命预测阶段，提出自约束状态空间估计器，利用自约束曲线为状态空间估计器在未来时刻的迭代更新提供伪观测，提高预测精度，获得更准确的剩余寿命预测结果。

3、本发明可以通过以下技术方案来实现：

4、一种海上风电机组高速轴承剩余寿命预测方法，包括：

5、步骤一：在海上风电机组的齿轮箱高速轴两侧轴承对应齿轮箱体的轴承座上以及发电机轴承座上的若干个采样点分别安装n个振动传感器，采集各个采样点振动信号，其中n为大于2的自然数；随后将振动信号由模拟值转换成数字信号，并进行预处理；

6、步骤二：将步骤一预处理后的采样点振动信号进行转换，获得每个采样点振动信号的时域特征和频域特征，并将这些特征保存为高速轴承的振动基础特征集；

7、步骤三：对步骤二获得的基础特征集分别进行单调性度量计算，将具有较高单调性的特征确定为轴承退化的敏感特征；

8、步骤四：基于步骤三所获得的敏感特征计算特征退化率，并构建递增型健康指标；

9、步骤五：通过指数核空间转换，将健康指标转换成核空间指标，确定核空间指标阈值，进而确定健康指标的失效阈值；

10、步骤六：基于自约束状态空间估计器预测未来时刻的健康指标变化趋势，将健康指标预测结果与失效阈值进行比较，确定轴承剩余寿命。

11、进一步的，步骤三中，单调性度量的表达式：

12、

13、

14、式中，mon是单调性度量，f(t)是t时刻的特征值。

15、进一步的，步骤四中，特征退化率的表达式：

16、

17、

18、式中，fi(t)是t时刻第i个敏感特征的特征值，di(t)是t时刻第i个敏感特征的特征退化率，i是t时刻敏感特征的数量，d(t)是t时刻所有敏感特征的特征退化率，tert是磨合结束时刻；

19、健康指标构建的表达式：

20、

21、式中，fb是其中一个基础特征，表示系统的基础健康状态；h(t)是t时刻的健康指标。

22、进一步的，步骤五中，指数核空间转换的表达式：

23、

24、式中，p是超参数，ksi(t)是t时刻核空间指标，tfdt是退化起始点。

25、进一步的，步骤五中，核空间指标阈值和核空间指标残差的表达式：

26、

27、

28、j(x,y)＝|time span of xi＜y| (5)

29、式中，t(l)end是ksi(l)的最后时刻，ksit是核空间指标阈值，ksie是核空间指标残差。

30、进一步的，步骤六中，自约束算法的表达式：

31、

32、式中，g(.)是约束曲线，w是约束曲线的参数，t是计算时刻，tnow是当前时刻，yt是t时刻历史观测，是约束曲线的拟合参数，τ是未来区间，是(tnow+τ)时刻的伪观测。

33、进一步的，将自约束曲线融入状态空间估计器，利用历史数据回归得到与历史数据发展趋势相符合的自约束曲线，为状态空间估计器在未来时刻的迭代更新提供伪观测，提高预测精度，获得更准确的剩余寿命预测结果。

34、本发明的有益效果在于：

35、1)现有健康指标构建方法没有关注到轴承退化过程为一不可逆过程，因此构建的指标不具有严格单调性。本发明针对这一问题，提出一种新的特征单调性趋势评价指标，即满足轴承真实退化情况，同时具有良好的趋势性及单调性，易于预测。

36、2)健康指标构建的难点之一是在避免对数据进行归一化处理的情况下消除不同特征间的量纲差异对指标构建结果产生的影响，本发明基于特征退化率有效提取不同特征包含的退化信息；在此基础上，考虑到轴承退化为一连续过程，建立了一种递推型的健康指标构建方法，并以正常运行阶段振动信号均方根的均值作为初始健康指标值表征不同轴承的正常状态；为解决失效阈值难确定问题，提出指数核空间转换方法，并基于核空间指标失效阈值确定健康指标的失效阈值，从而实现了仅利用历史振动数据进行轴承剩余寿命预测，具有较大工程实际意义。

37、3)在构建趋势性良好的轴承退化健康指标的基础上，基于对现有状态空间估计器的优缺点分析，提出了自约束状态空间估计器预测模型，其主要思想是利用自约束曲线为状态空间估计器在未来时刻的迭代更新提供伪观测，提高预测精度，获得更准确的剩余寿命预测结果。