CN104865400B - 一种风电机组转速的检测识别方法及系统 - Google Patents
一种风电机组转速的检测识别方法及系统 Download PDFInfo
- Publication number
- CN104865400B CN104865400B CN201510175876.4A CN201510175876A CN104865400B CN 104865400 B CN104865400 B CN 104865400B CN 201510175876 A CN201510175876 A CN 201510175876A CN 104865400 B CN104865400 B CN 104865400B
- Authority
- CN
- China
- Prior art keywords
- signal
- speed
- frequency
- rotating speed
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Abstract
本发明公开了一种风电机组转速的检测识别方法及系统,通过测取风电机组发电机轴承座振动加速度信号,首先去除信号中的趋势项和进行带通滤波,滤波通带的上下截止频率分别取双馈式风电机组发电机轴的最低和最高旋转频率;对滤波后的信号做奇异值降噪处理,提出带通信号中的突出旋转频率成分;最后对滤波降噪后的信号进行离散傅里叶变换,利用内插傅里叶变换法得到转子旋转频率的估计值。该方法可以给出发电机转子旋转频率的精确估计,为不同部件的故障特征提取提供必要的转速信息,不需要通过安装转速测量系统来获取转子转速信号,简化了振动监测系统,降低了系统成本,易于工程实施。
Description
技术领域
本发明机械设备运行参数的检测技术领域,尤其涉及一种风电机组转速的检测识别方法及系统。
背景技术
风力发电机组是实现将风能转换成电能的设备,其核心部分是由风轮、传动链、发电机构成的能量转换和传递系统。风轮将风能转换成旋转机械能,经过传动链传给发电机转子,发电机实现旋转机械能到电能的转换。风力发电机组通常安装在较差的地区,而且风速、载荷和风轮转速连续变化,导致机组运行工况复杂,极易造成风轮、传动链、发电机等承受载荷部件出现故障甚至损坏,而且风电场机组分布范围广,出现故障不易及时发现和维修,影响机组的安全经济性。随着风电机组投运数量、单机容量的不断增长,以及海上风电规模的逐渐扩大,提升风电机组检修技术水平成为风电运行企业的迫切需求,因此风电机组运行状态监测与故障诊断技术逐渐引起高度关注。风电机组的机械故障主要发生在风轮至发电机的整个传动链上的各个零部件,包括主轴轴承、增速齿轮箱内的各级齿轮及轴承、发电机转子轴承等。与一般机械设备中的传动链相比,风电机组传动链的载荷状态、运行工况、环境条件和结构布局等方面均比较复杂,导致传动链中齿轮、轴承等主要零部件的故障机理和故障发展模式等存在特殊性,因此故障率高于其它行业的同类设备,实际使用寿命远远低于设计寿命。
转速是旋转机械设备的主要参数之一,也是旋转机械设备故障诊断的一个重要参数,特别是对于风电机组这类变速变载旋转机械设备,其振动特性受转速的影响非常明显。有无转速信息,对于机组设备运行状态的判断及故障特征提取的效果影响非常大。大型风电机组设备通常在风轮低速主轴和发电机高速转子末端安装有转速测点,将转速信号接入机组监控系统。振动监测系统一般为独立运行的系统,主要监测各个部件的振动加速度信号和必要的运行参数,如机组的输出功率、转速等。为了避免影响机组监控系统的正常运行,振动监测系统的信号通常取自独立安装的各种传感器测点。振动加速度传感器的现场安装比较方便,而测量转速的传感器安装相对困难,受到机组现场条件的限制。对于便携式振动监测系统,转速传感器的安装和信号获取问题更加突出。
发明内容
为了解决上述问题,本发明提出一种风电机组转速的检测识别方法,包括:
步骤1、通过安装在发电机转子前后两端支撑轴承座上的振动加速度传感器来检测发电机结构振动加速度信号;对振动加速度信号进行预处理,去除信号中的趋势项和直流分量;
步骤2、对信号进行窄带滤波处理,窄带滤波的上下截止频率分别取双馈式风电机组发电机的最低和最高旋转频率;
步骤3、对滤波后的信号做奇异值降噪处理,提高窄带信号的信噪比,突出振动信号中的旋转频率成分;
步骤4、对降噪后的信号进行离散傅里叶变换,得到信号的离散频谱,利用内插法得到转子旋转频率的估计值;
步骤5、根据旋转频率估计值,计算发电机转子的转速,并根据增速齿轮箱内的传动齿轮结构和齿轮齿数,计算齿轮箱中各级齿轮轴的转速以及风轮低速主轴的转速。
所述步骤1中某一端支撑轴承座上的振动加速度信号x(n)经过窄带滤波处理后得到带通信号将其当做一个噪声干扰的正弦信号,即:
式中N为有限长信号的数据点数,u(n)为干扰噪声,f0为信号频率,A0为正弦信号的幅值、为正弦信号的初始相位
所述步骤3中包括对滤波信号进行相空间重构,得到重构相空间矩阵A,然后对A进行奇异值分解,奇异值分解公式为:UTAV=Σ,U、V分别为分解矩阵;分解得到信号的奇异值矩阵Σ,其中包含N个按大小排序的奇异值;保留Σ中最大的K个奇异值,K<N,然后做反变换,即A′=UΣVT,由反变换得到的相空间矩阵A′进行重构信号,从而获得了降噪信号。
所述步骤4对经过窄带滤波和奇异值降噪后的信号进行离散傅里叶变换,得到信号的离散频谱|X(m)|;如果信号不满足整周期截取的条件,离散幅值谱中的峰值点m偏离信号频率f0,存在偏差θ;采用内插法求出偏差的估计值偏差估计值的计算公式为:
式中,如果|X(m+1)|>|X(m-1)|,取变量α=1,否则取变量α=-1,m为正整数。
所述步骤4中根据偏差估计值求出旋转频率的估计值计算公式为:
其中,频率估计值的最大误差Δf为离散频谱的频率间隔,fs为离散信号的采样频率,m为正整数,为偏差的估计值;N为有限长信号的数据点数。
所述整周期截取的条件为:式中,B为截短窗宽度,T0为信号的周期,N为有限长信号的数据点数,f0为信号的频率,fs为离散信号的采样频率,m为正整数,θ为偏差。
一种风电机组转速的检测识别系统,包括:振动传感器、数据采集器、振动监测工作站,振动传感器安装在发电机转子前后两端支撑轴承座上,数据采集器分别与振动传感器和振动监测工作站相连;
所述振动监测工作站包括依次相连的:窄带滤波模块、奇异值降噪模块、频率估算模块、转速计算模块;
其中,窄带滤波模块对采集到的振动加速度信号进行窄带滤波处理,窄带滤波的上下截止频率分别取双馈式风电机组发电机的最低和最高旋转频率;
奇异值降噪模块对滤波后的信号做奇异值降噪处理,提高窄带信号的信噪比,突出振动信号中的旋转频率成分;
频率估算模块对降噪后的信号进行离散傅里叶变换,得到信号的离散频谱,利用内插法得到转子旋转频率的估计值;
转速计算模块根据旋转频率估计值,计算发电机转子的转速,并根据增速齿轮箱内的传动齿轮结构和齿轮齿数,计算齿轮箱中各级齿轮轴的转速以及风轮低速主轴的转速。
本发明的有益效果在于:
1)本方法直接使用振动监测系统中的振动加速度测量信号提取转速信息,不需要通过安装转速测量系统来获取转子转速信号,简化了振动监测系统,降低了系统成本。
2)本方法只对风电机组结构振动加速度信号进行分析处理,就能得到旋转频率的比较精确估计值,并进一步换算成转速。旋转频率估计值可以达到较高的精度,满足振动监测与故障诊断的需求。
3)风电机组的结构振动是由许多激励源共同作用的结果,振动加速度信号包含复杂的频率成分,对旋转频率识别效果产生影响。采用窄带频率技术可以滤除与旋转频率无关的频率成分,而将旋转频率成分突出出来。提高转速估计的精度。
4)由于振动监测系统安装在实际运行的机组设备上,现场条件复杂,振动监测信号往往受到较强的干扰噪声的影响。其中随机干扰噪声分布在较宽的频域范围,对旋转频率估计精度的影响较大。采用奇异值分解降噪可以有效地抑制宽频带干扰噪声的影响,进一步突出旋转频率成分,提供旋转频率估计精度。
5)本发明的方法不需要在被检测设备上另外安装硬件设备,因此不会对运行造成影响,简单方便,易于工程实施。
附图说明
图1为转速检测方法的硬件系统示意图
图2为原始信号、带通滤波信号、奇异值降噪信号对比
图3为信号的离散幅值谱示意图
具体实施方式
下面结合附图,对优选实施例作详细说明。
一种风电机组转速的检测识别系统,如图1所示,包括:振动传感器、数据采集器、振动监测工作站,振动传感器安装在发电机转子前后两端支撑轴承座上,数据采集器分别与振动传感器和振动监测工作站相连;
所述振动监测工作站包括依次相连的:窄带滤波模块、奇异值降噪模块、频率估算模块、转速计算模块;
其中,窄带滤波模块对采集到的振动加速度信号进行窄带滤波处理,窄带滤波的上下截止频率分别取双馈式风电机组发电机的最低和最高旋转频率;
奇异值降噪模块对滤波后的信号做奇异值降噪处理,提高窄带信号的信噪比,突出振动信号中的旋转频率成分;
频率估算模块对降噪后的信号进行离散傅里叶变换,得到信号的离散频谱,利用内插法得到转子旋转频率的估计值;
转速计算模块根据旋转频率估计值,计算发电机转子的转速,并根据增速齿轮箱内的传动齿轮结构和齿轮齿数,计算齿轮箱中各级齿轮轴的转速以及风轮低速主轴的转速。
本发明的转速检测与估计的原理是,转子不平衡是发电机结构振动的主要激励源之一,由于转子不平衡或多或少总是存在,因此发电机结构振动中总是包含转子不平衡引起的周期振动成分,根据振动理论,该周期振动成分的基频即为发电机转子的旋转频率。通过安装在发电机转子两端支撑轴承座上的振动加速度传感器,可以拾取发电机结构振动信号,其中就包含转子旋转频率的信息。通过对测量振动信号的分析处理,可以给出发电机转子旋转频率的精确估计,从而进一步换算出发电机转子转速,以及传动链上各个旋转部件的转速,为不同部件的故障特征提取提供转速必要的转速信息。
本发明提供了一种风电机组转速检测方法,具体包括:
(一)在风电机组的发电机转子两端靠近支撑轴承部位分别安装压电振动加速度传感器。传感器输出信号经过信号调理电路进行放大和抗混叠滤波后,进行模数转换,得到发电机转子前后两端的结构振动加速度信号x(n),y(n)。采集得到的数字信号经过预处理后,存入计算机进行分析处理。
(二)对经过预处理的振动加速度信号x(n)进行数字带通滤波,得到带通滤波信号理想带通滤波器的频响函数表达式为:
式中,f1为下限截止频率,f2为上限截止频率。滤波截止频率根据风电机组发电机的工作转速范围进行设定。经过带通滤波后,振动加速度信号中的旋转频率成分凸显出来,可以视为一个噪声干扰的正弦信号,即:
式中N为有限长信号的数据点数,u(n)为干扰噪声。
(三)对带通滤波信号进行奇异值分解降噪,提高带通信号的信噪比。具体过程为,对信号在一定的延迟时间和嵌入维数下进行相空间重构,得到重构相空间矩阵A,然后对A进行奇异值分解。奇异值分解公式为:
UTAV=Σ
分解得到信号的奇异值矩阵Σ,其中包含若干个按大小排序的奇异值。保留Σ中最大的几个奇异值,然后做反变换,即
A′=UΣVT
有反变换得到的相空间矩阵A′进行重构信号,从而获得了降噪信号x′(n)。
图2示出实测风电机组发电机驱动端轴承座振动加速度信号x(n)、带通滤波信号奇异值降噪后信号x′(n)的对比,可见经过带通滤波和奇异值降噪后,信号中的旋转频率成分被提取出来。
(四)对经过带通滤波和奇异值降噪后的信号x′(n)进行离散傅里叶变换,得到信号的离散频谱。如果带通滤波降噪信号不满足整周期截取的条件,即:
式中,B为截短窗宽度,T0为信号的周期,f0为信号的频率,fs为离散信号的采样频率,m为正整数。
N点非整周期截取信号x′(n)的幅值谱表达式如下:
图3为信号x′(n)的离散幅值谱示意图。由于非整周期采样,使得频域离散峰值点m偏离f0,偏差为θ。如果直接采用峰值点m对应的频率作为旋转频率的估计值,由于存在频域泄漏问题,使得频率估计值的误差较大,最大误差可以达到:
对于非整周期采样的正弦信号,其离散幅值谱|X(k)|的能量主要集中在峰值点|X(m)|及其两个相邻点|X(m+1)|和|X(m-1)|上。利用这三点的幅值信息,通过内插法可以实现旋转频率的精确估计。内插法频率估计分以下三步实现:
(1)找出信号离散幅值谱中的最大值|X(m)|及其对应的点m。
(2)根据下式m点与旋转频率f0对应点偏差θ的估计值:
式中,如果|X(m+1)|>|X(m-1)|,取α=1,否则取α=-1。
(3)根据下式求出旋转频率的估计值
根据旋转频率估计值,可以按照下式计算得到发电机转子的转速估计值,实现转速识别
(五)对连续采集的振动信号,取每N点按照上述(一)至(四)的步骤进行转速估计,可以得到风电机组转速的连续变化曲线。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (1)
1.一种风电机组转速的检测识别方法,其特征在于,包括:
步骤1、通过安装在发电机转子前后两端支撑轴承座上的振动加速度传感器来检测发电机结构振动加速度信号;对振动加速度信号进行预处理,去除信号中的趋势项和直流分量;
步骤2、对信号进行窄带滤波处理,窄带滤波的上下截止频率分别取双馈式风电机组发电机的最低和最高旋转频率;
步骤3、对滤波后的信号做奇异值降噪处理,提高窄带信号的信噪比,突出振动信号中的旋转频率成分;
步骤4、对降噪后的信号进行离散傅里叶变换,得到信号的离散频谱,利用内插法得到转子旋转频率的估计值;
步骤5、根据旋转频率估计值,计算发电机转子的转速,并根据增速齿轮箱内的传动齿轮结构和齿轮齿数,计算齿轮箱中各级齿轮轴的转速以及风轮低速主轴的转速;
所述步骤1中某一端支撑轴承座上的振动加速度信号x(n)经过窄带滤波处理后得到带通信号将其当做一个噪声干扰的正弦信号,即:
式中N为有限长信号的数据点数,u(n)为干扰噪声,f0为信号频率,A0为正弦信号的幅值,为正弦信号的初始相位;
所述步骤3中包括对滤波信号进行相空间重构,得到重构相空间矩阵A,然后对A进行奇异值分解,奇异值分解公式为:UTAV=Σ,U、V分别为分解矩阵;分解得到信号的奇异值矩阵Σ,其中包含N个按大小排序的奇异值;保留Σ中最大的K个奇异值,K<N,然后做反变换,即A′=UΣVT,由反变换得到的相空间矩阵A′进行重构信号,从而获得了降噪信号;
所述步骤4对经过窄带滤波和奇异值降噪后的信号进行离散傅里叶变换,得到信号的离散频谱|X(m)|;如果信号不满足整周期截取的条件,离散幅值谱中的峰值点m偏离信号频率f0,存在偏差θ;采用内插法求出偏差的估计值偏差估计值的计算公式为:
式中,如果|X(m+1)|>|X(m-1)|,取变量α=1,否则取变量α=-1,m为正整数;
所述步骤4中根据偏差估计值求出旋转频率的估计值计算公式为:
其中,频率估计值的最大误差Δf为离散频谱的频率间隔,fs为离散信号的采样频率,m为正整数,为偏差的估计值;N为有限长信号的数据点数;
所述整周期截取的条件为:式中,B为截短窗宽度,T0为信号的周期,N为有限长信号的数据点数,f0为信号的频率,fs为离散信号的采样频率,m为正整数,θ为偏差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510175876.4A CN104865400B (zh) | 2015-04-14 | 2015-04-14 | 一种风电机组转速的检测识别方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510175876.4A CN104865400B (zh) | 2015-04-14 | 2015-04-14 | 一种风电机组转速的检测识别方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104865400A CN104865400A (zh) | 2015-08-26 |
CN104865400B true CN104865400B (zh) | 2018-08-28 |
Family
ID=53911377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510175876.4A Expired - Fee Related CN104865400B (zh) | 2015-04-14 | 2015-04-14 | 一种风电机组转速的检测识别方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104865400B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105487010A (zh) * | 2015-11-21 | 2016-04-13 | 成都科瑞信科技有限责任公司 | 基于带通滤波的高效测试转速自动调整型电机测试系统 |
CN105548595B (zh) * | 2015-12-18 | 2019-04-30 | 国网河北能源技术服务有限公司 | 一种提取风电齿轮箱各级轴转速检测方法 |
CN106092524B (zh) * | 2016-05-13 | 2018-11-13 | 长兴昇阳科技有限公司 | 一种使用振动信号精确提取转速信号的方法 |
DE102016215635B4 (de) * | 2016-08-19 | 2022-02-03 | Robert Bosch Gmbh | Vorrichtung und Verfahren zur Bestimmung einer Drehzahl eines rotierenden Walzenkörpers |
CN107037327A (zh) * | 2016-10-09 | 2017-08-11 | 中国电力科学研究院 | 局部放电故障判定特征提取方法及判定方法 |
US10266278B2 (en) * | 2016-10-11 | 2019-04-23 | Unison Industries, Llc | Starter issue detection |
CN107192843B (zh) * | 2017-06-28 | 2019-05-24 | 北京航天控制仪器研究所 | 一种高精度转台的低速速率精度检测方法 |
CN108680247B (zh) * | 2018-05-16 | 2020-07-31 | 江苏大学镇江流体工程装备技术研究院 | 基于振动烈度低频滤波修正的振动信号转换方法 |
CN108932940B (zh) * | 2018-07-05 | 2021-04-27 | 浙江众邦机电科技有限公司 | 一种基于运转工况的工业缝纫机降噪方法、装置及设备 |
CN109813547B (zh) * | 2019-01-08 | 2020-09-22 | 华南理工大学 | 基于稀疏分解优化算法的旋转机械局部型故障诊断方法 |
EP3736970A1 (en) * | 2019-05-06 | 2020-11-11 | Siemens Gamesa Renewable Energy A/S | Position and speed calculation for an electric generator |
CN110987404B (zh) * | 2019-12-23 | 2022-06-10 | 安徽容知日新科技股份有限公司 | 一种转速识别方法及计算设备 |
CN112684200A (zh) * | 2020-12-09 | 2021-04-20 | 江苏大学 | 一种应用于海水淡化系统的转速监测仪 |
CN113252929B (zh) * | 2021-07-05 | 2022-02-01 | 格创东智(深圳)科技有限公司 | 转速确定方法、装置、电子设备及计算机可读存储介质 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003259695A (ja) * | 2002-03-06 | 2003-09-12 | Katsuhiro Sakai | 小型風力発電機の制御方法 |
JP4007268B2 (ja) * | 2003-07-22 | 2007-11-14 | 株式会社日立製作所 | 風力発電装置 |
CN201288636Y (zh) * | 2008-10-23 | 2009-08-12 | 皇明太阳能集团有限公司 | 小型风力发电机转速检测与控制装置 |
CN102590543A (zh) * | 2012-02-24 | 2012-07-18 | 中国汽车技术研究中心 | 一种转速测量装置及其测量方法 |
CN102661786B (zh) * | 2012-05-22 | 2013-09-04 | 国电联合动力技术有限公司 | 变转速风电机组转子振动信号的采集分析方法 |
CN103116032A (zh) * | 2013-01-18 | 2013-05-22 | 北京金风科创风电设备有限公司 | 用于获取风力发电机组的转速的方法和设备 |
CN203479392U (zh) * | 2013-09-26 | 2014-03-12 | 云南电网公司临沧供电局 | 一种输电塔振动监测系统 |
CN103743922B (zh) * | 2014-01-26 | 2016-03-02 | 深圳科瑞技术股份有限公司 | 一种基于dsp的发动机转速测量方法和装置 |
CN104374939B (zh) * | 2014-11-06 | 2017-02-15 | 西安交通大学 | 基于振动信号同步压缩变换的旋转机械瞬时转速估测方法 |
-
2015
- 2015-04-14 CN CN201510175876.4A patent/CN104865400B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104865400A (zh) | 2015-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104865400B (zh) | 一种风电机组转速的检测识别方法及系统 | |
CN102156043B (zh) | 风力发电机组在线状态监测与故障诊断系统 | |
Zhao et al. | Health assessment of rotating machinery using a rotary encoder | |
Peng et al. | Sensor fault detection and isolation for a wireless sensor network-based remote wind turbine condition monitoring system | |
Yoon et al. | On the use of a single piezoelectric strain sensor for wind turbine planetary gearbox fault diagnosis | |
Urbanek et al. | Comparison of amplitude-based and phase-based method for speed tracking in application to wind turbines | |
Daneshi-Far et al. | Review of failures and condition monitoring in wind turbine generators | |
CN105784366A (zh) | 一种变转速下的风电机组轴承故障诊断方法 | |
CN104819841B (zh) | 基于内置编码信息的单传感柔性角度域平均方法 | |
Guo et al. | Gear fault diagnosis of wind turbine based on discrete wavelet transform | |
Lin et al. | A review and strategy for the diagnosis of speed-varying machinery | |
Peng et al. | Speed estimation in planetary gearboxes: A method for reducing impulsive noise | |
Teng et al. | Detection and quantization of bearing fault in direct drive wind turbine via comparative analysis | |
US9146100B2 (en) | Wind turbine condition monitoring method | |
Hong et al. | Application of spectral kurtosis and Vold-Kalman filter based order tracking in wind turbine gearbox fault diagnosis | |
Wang et al. | Tacholess order-tracking approach for wind turbine gearbox fault detection | |
Guo et al. | Fault diagnosis of wind turbine bearing using synchrosqueezing wavelet transform and order analysis | |
CN112465068A (zh) | 一种基于多传感器数据融合的旋转设备故障特征提取方法 | |
CN110219816A (zh) | 用于风机故障诊断的方法和系统 | |
Moghadam et al. | Experimental validation of angular velocity measurements for wind turbines drivetrain condition monitoring | |
Liu et al. | Bispectrum analysis for feature extraction of pitting fault in wind turbine gearbox | |
CN104459186A (zh) | 基于稀疏分段拟合与积分逼近的无转速计阶比分析方法 | |
CN113218653B (zh) | 一种基于消除低频扰动的海上风电齿轮箱监测方法及装置 | |
Crabtree et al. | Condition monitoring algorithm suitable for wind turbine use | |
CN112664379A (zh) | 水轮机组故障预判方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180828 Termination date: 20210414 |
|
CF01 | Termination of patent right due to non-payment of annual fee |