基于定子电流数据驱动的双馈风电机组桨叶不平衡检测方法技术

一种基于定子电流数据驱动的双馈风电机组桨叶不平衡检测方法，在桨叶不平衡条件下，采集多组恒转速运行阶段的定子相电流iA和转速r信息，并进行希尔伯特变换，得到幅值包络信号，对包络信号进行FFT变换，提取叶轮1倍转频处幅值；利用最小二乘法拟合桨叶平衡时1倍转频处幅值与定子相电流的关系；机组运行中，采集恒转速运行阶段定子相电流iA和转速r，并计算出1倍频的幅值Q(1)；根据采集的iA并结合拟合的曲线，计算出Q_base(1)，若满足Q(1)>Q_base(1)+offset，则判断不平衡存在，且根据差值大小给出不平衡程度。本发明专利技术简单易行、准确性良好、适用性较好。

【技术实现步骤摘要】

本专利技术属于风电机组的在线故障诊断
，涉及一种基于定子电流数据驱动的双馈风电机组桨叶不平衡检测方法。
技术介绍
风电机组的运行环境通常都比较恶劣，叶片更会因风沙、结冰等造成质量不平衡，特别是近年来机组容量的增大及低风速机组的推广应用，使得叶片越来越长。叶轮的不平衡，包括气动不平衡，易传递到传动链，造成整机的振动加剧，严重时甚至导致事故的发生。考虑到不同叶片的差异性、故障机理的复杂性，难以建立较准确地诊断模型。采用传统类似固定阈值的判断方法，因阈值的确定依赖经验及少量试验，难以满足对不同风场、不同叶片的适应性。目前风力发电机组桨叶不平衡的检测可通过视频分析、振动检测、电信息等手段，但是这些方法大都需要增加设备，不仅成本大，而且可靠性不高、适应性较差。近来各研究机构及风机制造商也对桨叶不平衡进行了相关研究，以及公布了一些专利。如“叶轮不平衡故障的双馈风力发电机定子电流诊断方法”，该方法对特征分量的提取需经二次求导重构定子电流和二次FFT运算，处理过程复杂。又如“一种双馈风力发电机叶片不平衡在线故障诊断方法”，是对电功率进行处理，需要采集电压和电流量，并且这些方法理论上可行，但工程上面对不同地域的风场、不同厂家的叶片、不同的风机运行工况，适应性较差。
技术实现思路
为了克服已有双馈风力发电机组叶片不平衡检测方法在准确性、适应性方面的不足，本专利技术提供一种简单易行、准确性良好、适用性较好的基于定子电流数据驱动的双馈风电机组桨叶不平衡检测方法。本专利技术解决其技术问题所采用的技术方案是：一种基于定子电流数据驱动的双馈风电机组桨叶不平衡检测方法，所述检测方法包括如下步骤：1)桨叶平衡条件下，采集多组机组恒转速运行阶段的定子相电流iA和转速r信息，并进行希尔伯特变换，得到幅值包络信号，要求数据采集期间未进行变桨动作；2)对步骤1中的包络信号进行FFT变换，并根据发电机转速r计算叶轮1倍转频f(1)；3)根据步骤1)、2)的计算，得出叶轮1倍转频处的幅值大小，利用最小二乘法拟合出桨叶平衡时，1倍频处幅值与定子相电流的关系曲线；4)机组运行中，采集恒转速运行期间定子相电流iA和转速r，并按照步骤1)和2)中方法计算出1倍频的幅值Q(1)；5)根据步骤4)中采集的电流，并结合步骤3)中拟合的曲线，计算出Q_base(1)；6)比较Q(1)和Q_base(1)及offset的关系，若其满足Q(1)>Q_base(1)+offset，则判断不平衡存在，且根据差值大小给出不平衡程度，其中offset为残差补偿。进一步，所述步骤1)中，所述的恒转速运行阶段，转速方均根值在设定的合理范围内，所述合理范围为0-15之间。再进一步，在步骤2)中，根据发电机转速计算叶轮转速1倍频f(1)的方法如下：式(1)中，r为发电机转速，b为齿轮箱速比。更进一步，在步骤3)中，采用最小二乘法进行曲线拟合，方法如下：式中为待求函数，为某一函数类，a0,a1,…an为待求系数。式中，P表示求得的数据与实际数据满足误差的平方和的最小值。本专利技术的有益效果主要表现在：1)本专利技术涉及的信息的提取无需增加新的设备，硬件成本较低；2)因故障特征量本身很小，且不同工况下亦有不同。本专利技术通过拟合的方式，得出该条件下基准值，通过与基准值的比较得出诊断结论。进一步排除了噪声的干扰，效果显著，且通过误差补偿的方式，提高了检测率。附图说明图1为本专利技术流程图。图2为信号分析流程图。图3为机组正常运行模式下的转矩－转速曲线。图4为桨叶平衡时发电机定子A相电流和转速曲线。图5为桨叶平衡与不平衡时A相电流FFT分析对比曲线。图6为桨叶平衡与不平衡时Hilbert处理后频谱对比曲线。图7为桨叶平衡与不平衡时特征幅值与相电流有效值关系。图8为三种桨叶状态下，诊断结果的对比图。具体实施方式下面结合附图对本专利技术作进一步描述。参照图1～图8，一种基于定子电流数据驱动的双馈风电机组桨叶不平衡检测方法，该方法步骤如下：1)风电机组运行于恒转速附近，且叶片平衡时采集定子相电流iA和发电机转速r。在叶片不平衡时，会在定子电流中有呈现，如以质量不平衡为例进行介绍。假设某一叶片在距轮毂R处增加了一质量块m，此时的气动转矩可表示成：式中Tm为风力机输出转矩，ωm为叶轮转速，为初始相位。1.1)对(4)式，结合转子运动方程，可以得到电磁转矩中也含有振荡分量。式中Te为发电机电磁转矩，ΔTe不平衡转矩的幅值，为初始相位。1.2)在电网电压定向下，定子侧发电机惯例、转子侧电动机惯例，进行坐标变换，可以得到dq坐标下电磁转矩的表达式，进而可求得定子d轴电流isd。式中，isd为d轴电流，Δisd为震荡分量幅值，为初始相位。1.3)在1.2)条件下，dq轴系下定子电压方程和磁链可分别表示为式(7)和(8)，再有(6)式，可求出定子q轴电流。1.4)由定子dq轴电流，根据坐标变换即可求出abc坐标系下的值iA。从式(10)中，可以看到不平衡时会在定子电流中产生f1±fm次的谐波分量。中采集的数据进行Hilbert变换，对幅值包络信号FFT处理。对定子相电流iA做Hilbert变换，得到变换之后的信号iB，进一步得到解析信号Z(t)，对解析信号求幅值包络A(t)式中iB(t)是iA(t)的希尔伯特变换输出值。Z(t)＝iA(t)+jiB(t)(12)式中Z(t)为解析函数。式中A(t)为iA的包络幅值。2)对步骤1中的包络信号进行FFT变换，并根据发电机转速r计算叶轮1倍转频f(1)；方法如下：3)对包络信号A(t)进行FFT变换得到频谱信号，从中提取叶轮转速的1倍频幅值；利用最小二乘法拟合出桨叶平衡时，1倍频处特征幅值与定子相电流的关系：式中为待求函数，为某一函数类，a0,a1,…an为待求系数。式中，P表示求得的数据与实际数据满足误差的平方和的最小值；4)风机运行时，采集定子侧相电流和转速信息，按照步骤1)和2)中方法求出叶轮转速1倍频处特征幅值Q(1)；5)根据步骤4)中采集的定子电流数据，结合步骤3)中拟合的结果，求出该工况下对应的基准幅值Q_base(1)。6)比较Q(1)和Q_base(1)关系，若满足Q(1)>Q_base(1)+offset，则判断桨叶不平衡存在。本实施例涉及的参数均可从风力发电机的现有传感器中获取，无需增加新设备，监测成本低；因故障特征量本身很小，且不同工况下亦有不同。本专利技术通过拟合的方式，得出该条件下基准值，通过与基准值的比较得出诊断结论。进一步排除了噪声的干扰，效果显著，且通过误差补偿的方式，提高了检测率。以某三桨叶的1.5MW双馈风力发电机组为研究对象，机组运行于恒转速阶段时，即转速为1750rpm附近，运行区域如图3中所示，采集双馈机组恒转速运行附近的转速和定子相电流并进行信号分析，电流和转速时域波形如图4所示。图5是叶片平衡与不平衡时电流的频域图，可以明显看到在频率f＝50±0.275Hz处两者的区别，叶片不平衡时该处幅值偏大。为了进一步消除噪声、使特征量的提取更加准确，进行Hilbert变换，对幅值包络FFT分析，频谱如图6所示。图7是桨叶平衡与不平衡时特征幅值与定子A相电流的关系曲线，其中实线表示不平衡时，虚线表示平衡时，很明显看到随着定子侧输出功率的本文档来自技高网...

【技术保护点】
一种基于定子电流数据驱动的双馈风电机组桨叶不平衡检测方法，其特征在于：所述检测方法包括如下步骤：1)桨叶平衡条件下，采集多组机组恒转速运行阶段的定子相电流iA和转速r信息，并进行希尔伯特变换，得到幅值包络信号，要求数据采集期间未进行变桨动作；2)对步骤1中的包络信号进行FFT变换，并根据发电机转速r计算叶轮1倍转频f(1)；3)根据步骤1)、2)的计算，得出叶轮1倍转频处的幅值大小，利用最小二乘法拟合出桨叶平衡时，1倍频处幅值与定子相电流的关系曲线；4)机组运行中，采集恒转速运行期间定子相电流iA和转速r，并按照步骤1)和2)中方法计算出1倍频的幅值Q(1)；5)根据步骤4)中采集的电流，并结合步骤3)中拟合的曲线，计算出Q_base(1)；6)比较Q(1)和Q_base(1)及offset的关系，若其满足Q(1)>Q_base(1)+offset，则判断不平衡存在，且根据差值大小给出不平衡程度，其中offset为残差补偿。

【技术特征摘要】
1.一种基于定子电流数据驱动的双馈风电机组桨叶不平衡检测方法，其特征在于：所述检测方法包括如下步骤：1)桨叶平衡条件下，采集多组机组恒转速运行阶段的定子相电流iA和转速r信息，并进行希尔伯特变换，得到幅值包络信号，要求数据采集期间未进行变桨动作；2)对步骤1中的包络信号进行FFT变换，并根据发电机转速r计算叶轮1倍转频f(1)；3)根据步骤1)、2)的计算，得出叶轮1倍转频处的幅值大小，利用最小二乘法拟合出桨叶平衡时，1倍频处幅值与定子相电流的关系曲线；4)机组运行中，采集恒转速运行期间定子相电流iA和转速r，并按照步骤1)和2)中方法计算出1倍频的幅值Q(1)；5)根据步骤4)中采集的电流，并结合步骤3)中拟合的曲线，计算出Q_base(1)；6)比较Q(1)和Q_base(1)及offset的关系，若其满足Q(1)>Q_base(1)+o...

【专利技术属性】
技术研发人员：邱纪星，许国东，陶钧炳，韩小良，
申请(专利权)人：浙江运达风电股份有限公司，
类型：发明
国别省市：浙江;33