基于改进多胞型观测方法的风力机桨距执行器故障诊断方法技术

本发明专利技术提出一种基于改进多胞型观测方法的风力机桨距执行器故障诊断方法，主要包括：(1)建立风力机线性变参数模型，在MATLAB中进行模型搭建；（2）在无故障条件下运行MATLAB中的风力机模型，采集正常运行时的桨距角数据；（3）基于风力机线性变参数模型，利用多胞型未知输入观测方法实现干扰项的解耦，通过重构干扰矢量解耦干扰，从而消除未知干扰的影响；（4）通过配置状态反馈矩阵来减小噪声信号对系统的影响；（5）观测系统正常运行时的输出残差，设计阈值进行残差的评估，利用一段时间里残差信号的均方根值是否超过阈值来检测系统是否发生故障。本发明专利技术能够有效判断出桨距执行器故障，利于及时发现故障，确保系统稳定运行。

【技术实现步骤摘要】
基于改进多胞型观测方法的风力机桨距执行器故障诊断方法
本专利技术涉及到风力机系统故障诊断
，具体为一种利用改进多胞型未知输入观测方法实现桨距执行器故障诊断的方法。
技术介绍
风能作为资源丰富的可再生能源，清洁无污染，因为它具有取之不尽用之不竭的特性，所以风能利用已经成为世界各国用于解决能源危机、环境污染问题的有效方法。风力机是实现风能转换成电能的有效途径。桨距系统属于风力机系统中重要的组成部分，而桨距执行器又是桨距系统中的关键部件，该部件主要负责调整叶片桨距角。桨距执行器是一个特别容易出现故障的部件，一旦桨距执行器出现故障，就会造成桨距角得不到及时的调整或者无法调整，导致桨距系统的动态滞后甚至造成系统失控，对快速随机变化的风速反应不及时，造成工作在额定风速以上时输出功率不稳定。为确保风力机系统安全稳定的运行，有效阻止故障的传播和扩散，减少风力机系统运行维护的成本，及时的诊断和处理桨距执行器故障是必要的。利用观测方法实现故障诊断属于基于分析模型的方法，在对象模型的基础上根据输入信号和输出信号，结合相应的方法，诊断出系统当前存在的故障。
技术实现思路
本专利技术针对风力机桨距执行器故障诊断受干扰和噪声影响的问题，提出一种基于改进多胞型未知输入观测方法的故障诊断方法，实现对桨距执行器的故障诊断。所设计的改进多胞型观测方法能够在能在解耦干扰的同时抑制噪声信号的影响，并能准确诊断出桨距执行器的故障。按照本专利技术提供的技术方案，所述基于改进多胞型观测方法的风力机桨距执行器故障诊断方法包括5个步骤：步骤1、在风力机各个子模型的基础上建立线性变参数模型，在MATLAB中进行模型搭建；步骤2、在无故障条件下运行MATLAB中的风力机模型，采集风力机模型正常运行时的桨距角数据，并将采集到的桨距角运行数据离线保存；步骤3、基于风力机线性变参数模型，利用多胞型未知输入观测方法实现干扰项的解耦，通过重构干扰矢量解耦干扰，从而消除未知干扰的影响；步骤4、考虑到实际运行中噪声的影响，通过配置状态反馈矩阵来减小噪声信号对系统的影响；步骤5、观测系统正常运行时的输出残差，设计阈值进行残差的评估，利用一段时间里残差信号的均方根值是否超过阈值来检测系统是否发生故障。具体的，步骤1中，风力机系统包括空气动力子系统、传动子系统、桨距子系统和功率子系统，在风力机各个子模型的基础上建立风力机线性变参数模型；风机的空气动力学模型为传动系统模型为功率系统的模型为桨距系统模型为其中，ρ为空气密度常数，R为风轮转子半径，vr(t)为有效风速，在风速vr(t)的驱动下带动风轮转子以速度ωr(t)旋转，从而将风能转化为机械能；所转化的效率由功率系数Cp(λ(t),β(t))决定，功率系数Cp(λ(t),β(t))是叶尖速比λ(t)和桨距角β(t)的函数；传动机构由风轮转子侧的低速轴和发电机侧的高速轴组成，风轮捕获风中的动能，使低速轴转动产生机械转矩，然后通过变比为Ng的齿轮箱与高速轴连接，将机械转矩传递到与之相连的发电机上带动其旋转；Jr、Jg和Br、Bg分别为高、低速轴的转动惯量和摩擦系数，Bdt和Kdt为传动机构的扭力衰减系数和扭转刚度，θΔ(t)为传动系统的输出扭矩角，ωg(t)和Tg(t)分别为发电机速度和力矩；功率系统能够将机械能转换成电能并输送给电网，通过控制发电机力矩Tg(t)来跟踪其参考输出Tg,ref(t)，获得相应的功率；发电机产生的电能Pg(t)由发电机速度和负载大小决定，ηg是发电机的效率；β(t)和βref(t)分别为桨距角实际值和参考值，ωn和ζ分别为桨距执行器的自然频率和阻尼系数；风力机线性变参数模型为各矩阵在凸包范围内变化，即：其中ωn0为无故障时的自然频率，kr,β(θ)，Br(θ)，kr,v(θ)均为和力矩相关的气动参数；按照上述模型在MATLAB中进行风力机模型搭建。具体的，步骤3中，基于风力机线性变参数模型，利用多胞型未知输入观测方法实现干扰项的解耦，通过重构干扰矢量解耦干扰，从而消除未知干扰的影响；多胞型未知输入观测模型为z(t),x(t)∈Rn分别为观测模型的状态向量和系统状态向量，u(t)∈Rm为系统的控制输入向量，y(t)∈Rl为系统输出向量，d(t)∈Rr为系统未知输入，分别为x(t),y(t)的估计向量，ex(t)为状态偏差，N(θ),G(θ),L(θ)均为适维的参变矩阵存在矩阵M，MCD＝I，使得干扰项估计值的表达式为：令K＝I-HC，整理可得多胞型未知输入观测结构的表达式为式中H＝DM，N(θ)＝KA(θ)-L(θ)C，G(θ)＝KB(θ)，为待设计矩阵。具体的，步骤4中，考虑到实际运行中噪声的影响，通过配置状态反馈矩阵来减小噪声信号对系统的影响；状态估计偏差表示成下列形式状态估计偏差的协方差为p(t)＝ε{ex(t)ex(t)T}协方差的更新公式为p(t+1)＝N(θ)p(t)N(θ)T+KQKT+HRHT+L(θ)RL(θ)T,令E(θ)＝W(θ)p(t)CT(VT)-1当p(t+1)最小时滤波效果达到最优，此时L(θ)V-E(θ)＝0；求得反馈矩阵表示如下L(θ)＝W(θ)p(t)CT(Cp(t)CT+R)-1。步骤5中，采用前面4步搭建的观测器观测正常运行时桨距系统模型输出与它的故障模型输出对比产生的残差信号理想情况下，当系统正常运行时输出残差应该为零，故障情况下残差偏离零值。但实际运行中，由于建模偏差的影响，即使在正常情况下残差也不是严格收敛于零值的。因此，设计阈值进行残差的评估。阈值的设定采用RMS值的方法，即利用一段时间里残差信号的均方根值来检测系统的故障。具体的，将设定时间段残差信号均方根值的上确界定义为阈值；在实际运行中，当残差的均方根大于阈值时判断系统出现故障，当均方根值小于阈值时判断系统正常运行。本专利技术的优点如下：1、所设计的观测方法能够实现未知干扰解耦；2、所设计的观测方法对噪声具有鲁棒性；3、所设计的观测方法能够准确诊断出风力机的桨距执行器故障。附图说明图1为仿真风速。图2为仿真风速下的转矩。图3为仿真风速下的桨距角。图4为正常运行时的残差曲线。图5为白噪声下的桨距角输出。图6为泵堵塞时的观测残差。图7是本专利技术总体流程图。具体实施方式下面结合附图和实施例对本专利技术做进一步的说明。如图7所示，本专利技术总体包括以下步骤：步骤1、在风力机各个子模型的基础上建立线性变参数模型，在MATLAB中进行模型搭建；步骤2、在无故障条件下运行MATLAB中的风力机模型，采集风力机模型正常运行时的桨距角数据，并将采集到的桨距角运行数据离线保存；步骤3、基于风力机线性变参数模型，利用多胞型未知输入观测方法实现干扰项的解耦，通过重构干扰矢量解耦干扰，从而消除未知干扰的影响；步骤4、考虑到实际运行中噪声的影响，通过配置状态反馈矩阵来减小噪声信号对系统的影响；步骤5、观测系统正常运行时的输出残差，设计阈值进行残差的评估，利用一段时间里残差信号的均方根值是否超过阈值来检测系统的故障。风力机系统包括空气动力子系统、传动子系统、桨距子系统和功率子系统，在风力机各个子模型的基础上建立风力机线性变参数模型。实施例中以4.8MW三叶片水平轴风力机模型在MATLAB中进行仿真，仿真参数如表1所示：表1风力机仿真参数仿本文档来自技高网...

【技术保护点】
1.基于改进多胞型观测方法的风力机桨距执行器故障诊断方法，其特征是，包括以下步骤：步骤1、在风力机各个子模型的基础上建立线性变参数模型，在MATLAB中进行模型搭建；步骤2、在无故障条件下运行MATLAB中的风力机模型，采集风力机模型正常运行时的桨距角数据，并将采集到的桨距角运行数据离线保存；步骤3、基于风力机线性变参数模型，利用多胞型未知输入观测方法实现干扰项的解耦，通过重构干扰矢量解耦干扰，从而消除未知干扰的影响；步骤4、考虑到实际运行中噪声的影响，通过配置状态反馈矩阵来减小噪声信号对系统的影响；步骤5、观测系统正常运行时的输出残差，设计阈值进行残差的评估，利用一段时间里残差信号的均方根值是否超过阈值来检测系统是否发生故障。

【技术特征摘要】
1.基于改进多胞型观测方法的风力机桨距执行器故障诊断方法，其特征是，包括以下步骤：步骤1、在风力机各个子模型的基础上建立线性变参数模型，在MATLAB中进行模型搭建；步骤2、在无故障条件下运行MATLAB中的风力机模型，采集风力机模型正常运行时的桨距角数据，并将采集到的桨距角运行数据离线保存；步骤3、基于风力机线性变参数模型，利用多胞型未知输入观测方法实现干扰项的解耦，通过重构干扰矢量解耦干扰，从而消除未知干扰的影响；步骤4、考虑到实际运行中噪声的影响，通过配置状态反馈矩阵来减小噪声信号对系统的影响；步骤5、观测系统正常运行时的输出残差，设计阈值进行残差的评估，利用一段时间里残差信号的均方根值是否超过阈值来检测系统是否发生故障。2.如权利要求1所述基于改进多胞型观测方法的风力机桨距执行器故障诊断方法，其特征是，步骤1中，风力机系统包括空气动力子系统、传动子系统、桨距子系统和功率子系统，在风力机各个子模型的基础上建立风力机线性变参数模型；风机的空气动力学模型为传动系统模型为功率系统的模型为桨距系统模型为其中，ρ为空气密度常数，R为风轮转子半径，vr(t)为有效风速，在风速vr(t)的驱动下带动风轮转子以速度ωr(t)旋转，从而将风能转化为机械能；所转化的效率由功率系数Cp(λ(t),β(t))决定，功率系数Cp(λ(t),β(t))是叶尖速比λ(t)和桨距角β(t)的函数；传动机构由风轮转子侧的低速轴和发电机侧的高速轴组成，风轮捕获风中的动能，使低速轴转动产生机械转矩，然后通过变比为Ng的齿轮箱与高速轴连接，将机械转矩传递到与之相连的发电机上带动其旋转；Jr、Jg和Br、Bg分别为高、低速轴的转动惯量和摩擦系数，Bdt和Kdt为传动机构的扭力衰减系数和扭转刚度，θΔ(t)为传动系统的输出扭矩角，ωg(t)和Tg(t)分别为发电机速度和力矩；功率系统能够将机械能转换成电能并输送给电网，通过控制发电机力矩Tg(t)来跟踪其参考输出Tg,ref(t)，获得相应的功率；发电机产生的电能Pg(t)由发电机速度和负载大小决定，ηg是发电机的效率；β(t)和βref(t)分别为桨距角实际值和参考值...

【专利技术属性】
技术研发人员：吴定会，刘稳，郑洋，肖仁，黄海波，王佳宇，潘庭龙，沈艳霞，
申请(专利权)人：江南大学，
类型：发明
国别省市：江苏,32