一种变桨轴承失效的预警方法及系统与流程

本发明涉及风力发电机技术技术领域，具体地说是一种变桨轴承失效的预警方法及系统。

背景技术：

风电行业作为清洁能源在最近时期得到了快速发展，风力发电已经成为重要的电力能源解决方案。伴随着风电行业快速发展的同时，风机部件的可靠性也更加重要，因部件失效带来的损失也在逐年上升。在这方面，变桨轴承开裂问题尤其突出，已出现多起批量性质量问题，甚至导致风机倒塌，带来重大损失以及人员伤亡。由于变桨轴承结构以及运行工况复杂，导致其受力状况复杂，难以直接分析，外置振动检测传感器施工难度大，成本高，并且实施效果不理想。

如何建立有效便捷的预警机制，通过风机参数的计算来评估变桨轴承运行状态来尽早发现变桨轴承的异常状态，为变桨轴承预防性维护提供数据依据，减少机组带病运行的时间，从而避免机组的变桨轴承开裂导致的风电机组损坏和恶性事故发生，成为风电行业急需解决的问题。

技术实现要素：

本发明为了解决上述问题，提供了一种变桨轴承失效的预警方法及系统。

本发明解决其技术问题所采取的技术方案是：

本发明第一方面提供了一种变桨轴承失效的预警方法，该方法包括以下步骤：

对机组数据进行采集并存储；

判断机组数据采集周期是否结束；

如果结束，则对变点进行评估；

根据评估结果与预警阈值的对比判断是否报警。

结合第一方面，在第一方面第一种可能实现的方式中，所述的机组数据包括时间、风机编号、变桨过程五度接近开关值。

结合第一方面，在第一方面第二种可能实现的方式中，对变点进行评估具体包括对均值变点评估、对最小二乘法线性回归系数变点评估、对中位数分组回归的分组均值差与系数差变点评估。

结合第一方面，在第一方面第三种可能实现的方式中，对均值变点评估具体包括以下方法：

计算五度接近开关值均值式中，vk为五度接近开关值，k＝1，2，3...n。

结合第一方面，在第一方面第四种可能实现的方式中，对最小二乘法线性回归系数变点评估具体包括以下方法：

在风机变桨过程中获得五度接近开关值为vk；

设定一元线性回归函数：y＝a+bx；

第i个数据与拟合直线上相应的值之间残差为：

vi＝δyi＝yi-(a+bxi)；

vi＝δyi＝yi-(a+bxi)

最小二乘法拟合原理为为最小值，即

设置对a和b的一阶偏导数等于0，即：

其中

其中

结合第一方面，在第一方面第五种可能实现的方式中，对中位数分组回归的分组均值差与系数差变点评估具体包括以下方法：

将风机的在变桨过程中获得五度接近开关值按照由小到大排序获得样本数组为v′k，样本数量为n，k＝1，2，3......；v′n为中位数，则有中位数分组数组x1＝(v′1，v′2，...，v′n)；x2＝(vn+1，vn+2，...，vn)；

针对数组x1进行最小二乘法线性回归，设定一元线性回归函数为：

y＝a1+b1x

其中第i个数据vi与拟合直线上相应的值之间残差为：

vi＝δyi＝yi-(a1+b1xi)

最小二乘法拟合原理为为最小值，即

该结果需要对a1和b1的一阶偏导数等于0，即：

其中

其中

针对数组x2进行最小二乘法线性回归，设定一元线性回归函数为：

y＝a2+b2x

其中第i个数据vi与拟合直线上相应的值之间残差为：

vi＝δyi＝yi-(a2+b2xi)

最小二乘法拟合原理为为最小值，即：

该结果需要对a2和b2的一阶偏导数等于0，即

其中

其中

结合第一方面，在第一方面第六种可能实现的方式中，判断是否报警具体包括：在均值变点评估中，当e(vk)＞t1时(t1为均值变点预警阈值)，表明风机的变桨过程中五度接近开关值出现均值变点，进行变桨轴承失效预警报警；

在对最小二乘法线性回归系数变点评估中，当b＞t2(t2为最小二乘法线性回归系数变点预警阈值)，表明风机的变桨过程中五度接近开关值出现回归变点，进行变桨轴承失效预警报警；

在对中位数分组回归的分组均值差与系数差变点评估中，(t3为中位数分组均值差变点预警阈值)，表明风机的变桨过程中五度接近开关值出现变点，进行变桨轴承失效预警报警；b2-b1＞t4(t4为中位数分组系数差变点预警阈值)，表明风机的变桨过程中五度接近开关值出现变点，进行变桨轴承失效预警报警。

本发明第二方面提供了一种变桨轴承失效的预警系统，该系统包括：数据采集模块，用于采集机组数据；数据存储模块，用于存储机组数据；评估模块，用于对变点进行评估；比对模块，用于根据评估结果与预警阈值进行比对；判断模块，用于判断数据采集周期是否结束以及根据比对结果判断是否报警。

本发明第二方面的所述系统能够实现第一方面及第一方面的各实现方式中的方法，并取得相同的效果。

由以上技术方案可见，本发明通过分析计算风力发电机变桨过程中五度接近开关值的变化趋势，判断风力发电机变桨轴承运行异常情况，及时发现变桨轴承的异常，能够提前判断隐患，避免故障停机的发生。通过对变桨轴承变桨过程中的运行参数进行分析计算，实现难度小，成本低，而且评估效果有很大改善。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为一种变桨轴承失效的预警方法流程示意图

图2为本发明实施例所应用的一种变桨轴承失效的预警系统结构示意图。

具体实施方式

为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。

如图1所示，一种变桨轴承失效的预警方法，该方法包括以下步骤：

对机组数据进行采集；机组数据包括时间、风机编号、变桨过程五度接近开关值；

对采集后的数据进行存储；

判断机组数据采集周期是否结束；如果没有结束，则继续采集机组数据；

如果结束，则对变点进行评估；

根据评估结果与预警阈值的对比判断是否报警。

对变点进行评估包括对均值变点评估、对最小二乘法线性回归系数变点评估、对中位数分组回归的分组均值差与系数差变点评估。

对均值变点评估具体包括以下方法：

计算五度接近开关值均值式中，vk为五度接近开关值，k＝1，2，3...n。

对最小二乘法线性回归系数变点评估具体包括以下方法：

在风机变桨过程中获得五度接近开关值为vk；

设定一元线性回归函数：y＝a+bx；

第i个数据与拟合直线上相应的值之间残差为：

vi＝δyi＝yi-(a+bxi)；

vi＝δyi＝yi-(a+bxi)

最小二乘法拟合原理为为最小值，即

设置对a和b的一阶偏导数等于0，即：

其中

其中

对中位数分组回归的分组均值差与系数差变点评估具体包括以下方法：

将风机的在变桨过程中获得五度接近开关值按照由小到大排序获得样本数组为v′k，样本数量为n，k＝1，2，3......；v′n为中位数，则有中位数分组数组x1＝(v′1，v′2，...，v′n)；x2＝(v′n+1，v′n+2，...，v′n)；

针对数组x1进行最小二乘法线性回归，设定一元线性回归函数为：

y＝a1+b1x

其中第i个数据vi与拟合直线上相应的值之间残差为：

vi＝δyi＝yi-(a1+b1xi)

最小二乘法拟合原理为为最小值，即

该结果需要对a1和b1的一阶偏导数等于0，即：

其中

其中

针对数组x2进行最小二乘法线性回归，设定一元线性回归函数为：

y＝a2+b2x

其中第i个数据vi与拟合直线上相应的值之间残差为：

vi＝δyi＝yi-(a2+b2xi)

最小二乘法拟合原理为为最小值，即：

该结果需要对a2和b2的一阶偏导数等于0，即

其中

其中

判断是否报警具体包括：在均值变点评估中，当e(vk)＞t1时(t1为均值变点预警阈值)，表明风机的变桨过程中五度接近开关值出现均值变点，进行变桨轴承失效预警报警。

在对最小二乘法线性回归系数变点评估中，当b＞t2(t2为最小二乘法线性回归系数变点预警阈值)，表明风机的变桨过程中五度接近开关值出现回归变点，进行变桨轴承失效预警报警。

在对中位数分组回归的分组均值差与系数差变点评估中，(t3为中位数分组均值差变点预警阈值)，表明风机的变桨过程中五度接近开关值出现变点，进行变桨轴承失效预警报警；b2-b1＞t4(t4为中位数分组系数差变点预警阈值)，表明风机的变桨过程中五度接近开关值出现变点，进行变桨轴承失效预警报警。

如图2所示，一种变桨轴承失效的预警系统，该系统包括：数据采集模块，用于采集机组数据；数据存储模块，用于存储机组数据；评估模块，用于对变点进行评估；比对模块，用于根据评估结果与预警阈值进行比对；判断模块，用于判断数据采集周期是否结束以及根据比对结果判断是否报警。

以上所述仅是本发明的具体实施方式，使本领域技术人员能够理解或实现本发明。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。