一种风力发电机组塔筒振动异常检测方法与流程

本发明涉及一种风力发电机组塔筒振动异常检测方法。

背景技术：

塔筒是风力发电机组的关键部件。机舱底座固定在塔筒顶端，并与轮毂连接，后者用于固定叶片。风力发电机组几乎所有部件重力载荷以及动力载荷都会作用到塔筒，因此塔筒振动保护涉及到整机的安全，也是发电企业重点关注的问题。

目前，传统的风力发电机组塔筒振动异常检测方法主要是通过判断振动测量值是否超过固定的限值来实现。另外还有通过加装应力传感器，利用力载荷测量监测塔筒安全的方法。

中国专利201010565456.4提供了一种风电机组振动检测装置及方法。该装置包括数个振动传感器和至少一个振动开关，以及振动评估模块，用于根据预设的加权值和加权公式对振动信号进行处理，以获取风电机组振动状态数据，并与预设的故障值进行比较，判断是否触发振动故障保护逻辑。

中国专利201210302806.7提供一种风力发电机组低频振动的监测方法，首先利用加速度传感器采集塔筒振动信号，然后通过振动分析仪中的滤波器和检波器进行相应的处理，得到各种低频振动监测信号，并与PLC控制器中对应各个监测信号设置的报警阀值曲线的阀值进行对比，进而保护机组。

中国专利201610505126.3公开了一种风电机组的预警机构及预警方法，其中预警机构包括三向加速度传感器、振动信号实时调理系统、模态实时采集与谱分析系统、多个力传感器组件、多个应变信号实时调理系统、多个应变实时采集分析系统、实时频谱比较系统、实时应变比较系统、判断系统以及数据输出系统。

然而上述专利主要是基于设置的固定的振动保护阈值，忽略了现场实际机组由于维护时间、风况特点、地形地质、以及传感器安装或测量等因素引起的塔筒振动差异性。

技术实现要素：

为克服上述现有技术的缺点，本发明提出一种风力发电机组塔筒振动异常检测方法。

本发明利用风电场运行数据生成各台同型号机组的风速-振动表，绘制各相同风速段下全风电场机组振动值箱线图，以发现塔筒振动异常的风电机组；再将塔筒振动正常机组的风速-振动表合并成全风电场该型号机组的风速振动特征表，然后将各台机组当日生成的风速-振动表与全风电场该型号机组的风速振动特征表进行差异对比，实现当日对机组塔筒振动异常检测。

本发明包括以下步骤：

1、收集风电场各台风力发电机组近3个月的运行数据，提取其中风速5分钟均值、塔筒振动有效值的5分钟均值。其中塔筒振动有效值是塔筒振动瞬时值在一个塔筒振动周期内的均方根，塔筒振动周期可以从机组设计手册中查到。

2、对步骤1得到的每台风力发电机组风速5分钟均值、塔筒振动有效值的5分钟均值进行风速-振动散点拟合，首先对风速5分钟均值进行分仓，间隔为1m/s，然后对相同仓内振动数据进行均值计算，公式为：

和

其中，X_i为风速在bin区间i内的平均值，Y_i为振动在bin区间i内的平均值，x_i,j为风速在bin区间i内的第j个数值，y_i,j为振动在bin区间i内的第j个数值，N_i为数据在bin区间i内的有效点个数。

最后形成每台机组的风速-振动表，该表有8行1列，每行代表1个风速仓，分别对应3m/s至10m/s，第1列为X_i，第2列为Y_i。

3、采用箱线图异常值计算方法求取步骤2形成的全场各台机组风速-振动表在各风速段的振动均值正常范围，超出振动均值正常范围的风力发电机组则视为塔筒振动异常，至此完成对风场机组塔筒历史振动异常的检测。

箱线图异常值的计算方法为：首先将相同风速段下的所有机组振动均值从小到大排序，然后选出上四分位数Q3，中位数，下四分位数Q1，再计算上四分位数和下四分位数之间的差值Q3-Q1，又称IQR，绘制箱线图的矩形，上边缘为上四分位数，下边缘为下四分位数。在矩形内部中位数的位置画一条横线，即中位线，在矩形上边缘上方的1.5IQR范围内出现的最大数据点处画一条线，作为正常值上限，在矩形下边缘下方的1.5IQR范围内出现的最小数据点处画一条线，作为正常值下限。超出正常值限值的数据为异常点，用“+”表示。

4、剔除步骤3检测出的塔筒振动异常的风力发电机组，将风电场其余塔筒振动正常机组的风速-振动表合并，具体方法为将各风力发电机组风速-振动表相同位置表格内的数值相加，并除以机组数量，形成全风电场机组的风速振动特征表。

5、在风力发电机组运行时，每日采集0点至24点的风速5分钟均值、塔筒振动有效值的5分钟均值，并按照步骤2的方法形成机组当天的风速-振动表。

6、将步骤5形成的风力发电机组当天的风速-振动表与步骤4形成的全风电场机组的风速振动特征表进行对比，计算得到机组当天的振动水平指标，公式如下：

其中，g为机组当天的振动水平，vib_i,2为机组当天的风速-振动表第i行第2列的振动值，为全场机组的风速振动特征表第i行第2列的振动值。

如果机组当天的振动水平g大于阈值，则触发塔筒当日振动异常告警，至此完成对各台机组当日运行过程中塔筒振动异常检测。所述阈值通常取机组振动保护值的50％，例如某机组振动保护值为0.8m/s²，g_i则取0.4m/s²。

附图说明

图1本发明方法示意图；

图2全场机组在各风速段的振动均值的箱线图；

图3某机组当日风速-振动表与全场机组风速振动特征表绘制的对比曲线。

具体实施方式

下面结合附图和具体实施方式对本发明做进一步说明。

本发明方法步骤如下：

1、收集风电场各台风力发电机组近3个月的运行数据，提取其中风速5分钟均值、塔筒振动有效值的5分钟均值。其中塔筒振动有效值是塔筒振动瞬时值在一个塔筒振动周期内的均方根，塔筒振动周期可以从机组设计手册中查到。

2、对步骤1形成的每台风力发电机组风速5分钟均值、塔筒振动有效值的5分钟均值进行风速-振动散点拟合，最后形成每台机组的风速-振动表。

3、采用箱线图异常值计算方法求取步骤2形成的全场各台机组风速-振动表在各风速段的振动均值正常范围，超出振动均值正常范围的风力发电机组则视为塔筒振动异常，至此完成方法中对风场机组塔筒历史振动异常的检测。

4、剔除掉步骤3检测出的塔筒振动异常的风力发电机组，将风电场其余塔筒振动正常机组的风速-振动表合并，具体方法为将各风力发电机组风速-振动表相同位置表格内的数值相加，并除以机组数量，形成全风电场机组的风速振动特征表。

5、在风力发电机组运行时，每日采集0点至24点的风速5分钟均值、塔筒振动有效值的5分钟均值，并按照步骤2的方法形成机组当天的风速-振动表。

6、将风力发电机组当日形成的风速-振动表与全风电场机组的风速振动特征表进行比较，差异较大则触发塔筒当日振动异常告警。

以下结合实施例进一步说明本发明，如图1所示，

1、在步骤1中，对某容量5万千瓦风场33台风力发电机组的近3个月的运行数据进行收集，提取其中风速5分钟均值、塔筒振动有效值的5分钟均值。其中塔筒振动有效值是塔筒振动瞬时值在一个塔筒振动周期内的均方根，塔筒振动周期从机组设计手册中查到为3.18秒。

2、对步骤1形成的每台风力发电机组数据进行风速-振动散点拟合，首先对风速5分钟均值进行分仓，间隔为1m/s，然后对相同仓内振动数据进行均值计算，公式为：

和

其中，X_i为风速在bin区间i内的平均值，Y_i为振动在bin区间i内的平均值，x_i,j为风速在bin区间i内的第j个数值，y_i,j为振动在bin区间i内的第j个数值，N_i为数据在bin区间i内的有效点个数。

最后形成33台机组的风速-振动表，该表有8行1列，每行代表1个风速仓，分别对应3m/s至10m/s，第1列为X_i，第2列为Y_i。

3、采用箱线图异常值计算方法求取步骤2形成的全场各台机组风速-振动表在各风速段的振动均值正常范围，超出振动均值正常范围的2号和8号机组被检测出塔筒振动异常，至此完成对风场机组塔筒历史振动异常的检测。

全风电场机组在各风速段的振动均值的箱线图如图2所示。以风速段8.5m/s-9.5m/s的振动均值正常范围计算为例，首先将相同风速段下的所有机组振动均值从小到大排序，然后选出上四分位数Q3为0.0771m/s²，中位数为0.0503m/s²，下四分位数Q1为0.0427m/s²，接下来计算上四分位数和下四分位数之间的差值Q3-Q1为0.0344m/s²，又称IQR，绘制箱线图的矩形，上边缘为上四分位数，下边缘为下四分位数。在矩形内部中位数的位置画一条横线，即中位线，在矩形上边缘上方的1.5IQR范围内出现的最大数据点0.0872m/s²处画一条线，作为正常值上限，在矩形下边缘下方的1.5IQR范围内出现的最小数据点0.0313m/s²处画一条线，作为正常值下限。超出正常值范围的数据为异常点，用“+”表示。

4、剔除步骤3检测出的塔筒振动异常的2号和8号机组，将风电场其余塔筒振动正常机组的风速-振动表合并，具体方法为将各风力发电机组风速-振动表相同位置表格内的数值相加，并除以机组数量，形成全风电场机组的风速振动特征表。

5、在风力发电机组运行时，每日采集0点至24点的风速5分钟均值、塔筒振动有效值的5分钟均值，并按照步骤2的方法形成机组当天的风速-振动表。

6、将步骤5形成的风力发电机组当天的风速-振动表与步骤4形成的全场机组的风速振动特征表进行对比计算，得到机组当天的振动水平指标，公式如下：

其中，g为机组当天的振动水平，vib_i,2为机组当天的风速-振动表第i行第2列的振动值，为全场机组的风速振动特征表第i行第2列的振动值。

机组当日振动水平大于阈值，则触发塔筒当日振动异常告警。如图3所示，17号机组g等于0.0313，大于该风场设定阈值0.02，因此该机组于凌晨0点触发塔筒当日振动异常告警，并标识出对应类别号和对比波形，供专业人员查看分析。