本发明涉及风电机组故障预警，具体涉及一种基于集群分析的风电机组齿轮箱轴承故障预警方法。
背景技术：
风电机组运行环境较为恶劣，需要对风电机组定期检修和故障后维修，由于风电机组的特殊性，定期检修选择在小风时停机，对风电机组进行简单维护或更换耗材，无法发现大部件故障隐患，而在设计寿命周期内，风电机组一旦发生故障，需要停机维修，这会造成较大的经济损失。
目前，风电机组的故障统计数据表明，发生故障次数较多的部件主要是变流器、齿轮箱、变桨系统、控制柜、偏航系统等，根据部分风电场的运行数据统计，故障诊断、检修以及风电机组受累停运时间占故障损失运行时间的79％，而定期维护时间仅占12％。对于大部件如果仅采用故障后被动的维修方式，容易造成缺陷扩大，维修时间延长，长期停机，加剧经济损失，对于海上风电机组，海上运输、吊装、维修成本更加高昂。
我国风电整机装备技术起步较晚，早期生产和运行的风电机组主要采用国外技术直接引进或购买生产许可证生产，在技术引进、消化以及国产化的研发进程中，恰遇我国风电产业大发展，注重产量和产值粗犷式的发展同时，对技术细节关注不够，导致目前运行的风电机组故障率远远超过国外同类型风电机组。
目前风电机组均配有监控与数据采集系统(SCADA，supervisory control and data acquisition)，能够监控风电机组功率、转速，以及齿轮箱、发电机等大部件的温度等参数，对主要部件的运行状态进行简单的判定，但是一般只针对单一参数设定固定阈值，时效性较差。
本发明利用风电机组运行SCADA数据，对同一风电场内同型号的全部运行风机进行集群分析，通过风电机组发电功率、风轮转速、齿轮箱油温等参数与齿轮箱轴承温度之间的相对关系，对齿轮箱轴承故障进行预警，对故障预警具有更高的准确性和及时性，并且直接利用SCADA数据，不额外增加成本，占用运算资源少。
技术实现要素：
本发明提供一种风电机组齿轮箱轴承故障预警方法，该方法包括如下步骤：
步骤1：SCADA运行数据的采集；
步骤2：参数选择；
步骤3：参数相关性分析；
步骤4：确定目标变量和分析变量，进行目标变量的区间分析；
步骤5:目标变量的集群统计分析；
步骤6：风电机组状态的判断。
所述步骤1SCADA运行数据的采集包括：
采集同一风电场内同型号的单次运算数据周期不少于1天的至少10台风电机组的数据，数据选用10min平均值；
验证采集的SCADA运行数据，保证最终信号值的正确性。
所述步骤2参数选择包括：
从采集的SCADA运行数据中选出所需参数及其对应的数据；
所选参数包括：齿轮箱轴承温度、齿轮箱油温、平均风速、发电机转速、有功功率；
所述参数对应的数据的有效性核查包括：检查数据一致性，处理无效值、重复值和缺失值。
所述步骤3参数相关性分析包括：
对选出的参数进行相关性分析，根据下述算法中的任意一种计算参数间的相关系数r：
算法一：
算法二：
算法三：
其中，x、y为需要计算相关系数的2个参数，r为这2个参数之间的相关系数，n为参与计算的数据点个数，和均表示平均数，σx、σy、σxy、Lxx、Lyy和Lxy均为函数；
所述步骤4，选择齿轮箱轴承温度为目标变量Y，同一风电机组内与齿轮箱轴承温度相关系数最高的参数为分析变量X；对同一风电场内的m个风电机组进行目标变量区间分析包括：
将一个运算周期内分析变量X的最小值和最大值分别作为区间的起点和终点，将该区间等值划分为f个小区间x1,x2…xf；
分别计算每一个风电机组在每个小区间内的目标变量Y的平均值y1，y2…yf，其中的1、2……f为小区间标号，形成X-Y关系曲线。
所述步骤5目标变量集群统计分析包括：
每个小区间内的m个风电机组的目标变量Y的平均值如下所示：
第1个小区间的目标变量Y分析结果为：y1,1、y1,2……y1,m，其中的1、2……m为风电机组标号；
第f个小区间的目标变量Y分析结果为：yf,1,、yf,2……yf,m，其中的1、2……m为风电机组标号。
所述步骤5包括：
按下式计算同一小区间内m个风电机组目标变量Y的分析结果的平均值Af：
按下式计算同一小区间内m个风电机组目标变量Y分析结果的标准差Df：
其中，m为风电机组的个数，yf,i为第f个小区间内第i台风电机组的目标变量Y的平均值。
所述步骤6状态判断包括：
按下式计算各个区间的故障报警阈值δf：
δf＝Af+ε·Df
其中，ε为阈值设定系数，Af为同一小区间内m个风电机组目标变量Y的分析结果的平均值，Df为同一小区间内m个风电机组目标变量Y分析结果的标准差。
所述步骤6包括：
如小区间k内的风电机组j的目标变量Y的平均值yk,j超过该区间的报警阈值δk，则风电机组j在小区间k发生状态预警；
当风电机组j在f个小区间中有超过30％的小区间发生状态预警，则确定风电机组j发生齿轮箱轴承故障报警，30％为一般典型数值，实际应用中可根据情况调整。
风电机组齿轮箱轴承故障的预警结果的准确性与风电机组数量、运算数据周期、数据可靠性相关，同一风电场中用于分析的风电机组数量越多，预警结果越准确，运算数据周期选择1天-1个月之内较为合适，时间太短或太长均会影响预警结果。SCADA数据可靠性越高，预警结果约准确。
与最接近的现有技术比，本发明提供的技术方案具有以下优异效果：
1、使用SCADA数据对风电机组齿轮箱轴承故障进行预警，运算方法简明高效，不额外增加成本，极少占用系统资源。
2、可对SCADA数据进行有效选择，选取能够准确表达齿轮箱轴承运行状态的参数，对比一般监控系统的单一阈值判定，更为准确，时效性强。
3、使用区间分析法对风电机组的不同运行状态分别计算，更贴切的考虑到风电机组运行中的复杂情况。
4、可针对不同风电场，不同项目的实际情况对预警参数进行调整，具有很强的灵活性。
附图说明
图1为本发明参数相关性分析的结果示例；
图2为本发明区间分析的结果示例；
图3为本发明集群统计分析的结果示例；
图4为本发明状态判定的结果示例。
具体实施方式
下面结合附图对本发明做进一步详细的的说明：
本发明提供一种风电机组齿轮箱轴承故障预警方法，该方法包括如下步骤：
步骤1：SCADA运行数据的采集；
步骤2：参数选择；
步骤3：参数相关性分析；
步骤4：确定目标变量和分析变量，进行目标变量的区间分析；
步骤5:目标变量的集群统计分析；
步骤6：风电机组状态的判断。
所述步骤1SCADA运行数据的采集包括：
采集同一风电场内同型号的单次运算数据周期不少于1天的至少10台风电机组的数据，数据选用10min平均值；
验证采集的SCADA运行数据，保证最终信号值的正确性。
所述步骤2参数选择包括：
从采集的SCADA运行数据中选出所需参数及其对应的数据；
所选参数包括：齿轮箱轴承温度、齿轮箱油温、平均风速、发电机转速、有功功率；
所述参数对应的数据的有效性核查包括：检查数据一致性，处理无效值、重复值和缺失值。
如图1所示为本发明步骤3参数相关性分析的结果示例图，所述步骤3包括：
对选出的参数进行相关性分析，根据下述算法中的一种计算相关系数：
算法一：
算法二：
算法三：
其中，x、y为需要计算相关系数的2个参数，r为这2个参数之间的相关系数，n为参与计算的数据点个数，和均表示平均数，σx、σy、σxy、Lxx、Lyy和Lxy均为函数。
如图2所示为本发明步骤4区间分析的结果示例图，所述步骤4包括：选择齿轮箱轴承温度为目标变量Y，同一风电机组内与齿轮箱轴承温度相关系数最高的参数为分析变量X；本实施例选择齿轮箱轴承温度为目标变量Y，齿轮箱油温为分析变量X。
以一个运算周期内分析变量X的最小值和最大值分别作为区间的起点和终点，将该区间等值划分为f个小区间x1,x2…xf；
分别计算每一个风电机组在每个小区间内的目标变量Y的平均值y1，y2…yf，形成X-Y关系曲线，其中的1、2……f为小区间标号。
如图3所示为本发明步骤5集群统计分析的结果示例图，以齿轮箱轴承温度为目标变量Y，齿轮箱油温为分析变量X，风电机组总数m为10，所述步骤5包括：
分别对每个小区间内的m个目标变量Y的平均值进行分析，
第1个小区间的目标变量Y分析结果为：y1,1、y1,2……y1,m，其中1、2……m为风电机组标号；
第f个小区间的目标变量Y分析结果为：yf,1,、yf,2……yf,m，其中1、2……m为风电机组标号。
所述步骤5包括：
根据公式(7)计算同一小区间内m个风电机组目标变量Y分析结果的平均值A：
根据公式(8)计算同一小区间内m个风电机组目标变量Y分析结果的标准差D：
其中，Af为第f个小区间内m个风电机组目标变量Y的平均值的平均值，Df为第f个小区间内m个风电机组目标变量Y的平均值的标准差。
如图4所示为本发明步骤6的结果示例图，以齿轮箱轴承温度为目标变量Y，齿轮箱油温为分析变量X，风电机组总数m为10，所述步骤6包括：
设定各个区间的故障报警阈值为δ，如下式(9)所示：
δf＝Af+ε·Df (9)
其中，ε为阈值设定系数，典型值为2，可根据实际情况调整；δf为第f个小区间的故障报警阈值。
所述步骤6包括：
如果某个小区间k内的风电机组j的目标变量Y分析结果yk,j超过该区间的报警阈值δk，则风电机组j在区间k发生状态预警；
当同一风电机组j有超过30％的小区间发生预警，则确定风电机组j发生齿轮箱轴承故障报警，30％为一般典型数值，实际应用中可根据情况调整。
风电机组齿轮箱轴承故障的预警结果的准确性与风电机组数量、运算数据周期、数据可靠性相关，同一风电场中用于分析的风电机组数量越多，预警结果越准确，运算数据周期选择1天-1个月之内较为合适，时间太短或太长均会影响预警结果。SCADA数据可靠性越高，预警结果约准确。
所述方法所用装置包括：
用于数据采集的SCADA系统，用于数据分析的数据处理模块和用于报警的报警模块。
最后应当说明的是：以上实施例仅用于说明本发明的技术方案而非对其保护范围的限制,尽管参照上述实施例对本申请进行了详细的说明,所属领域的普通技术人员应当理解：本领域技术人员阅读本申请后依然可对申请的具体实施方式进行种种变更、修改或者等同替换，但这些变更、修改或者等同替换，均在申请待批的权利要求保护范围之内。