CN112782421A - 一种基于音频的转速识别方法 - Google Patents
一种基于音频的转速识别方法 Download PDFInfo
- Publication number
- CN112782421A CN112782421A CN202011478552.5A CN202011478552A CN112782421A CN 112782421 A CN112782421 A CN 112782421A CN 202011478552 A CN202011478552 A CN 202011478552A CN 112782421 A CN112782421 A CN 112782421A
- Authority
- CN
- China
- Prior art keywords
- signal
- rotating speed
- value
- autocorrelation function
- zero
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
Abstract
本申请及一种基于音频的转速识别方法,仅需通过音频采集装置采集声音信号,记录风机叶片扫风声音,然后分别通过过零次数、自相关函数、功率谱密度算法迅速计算当前叶片的实时转速,互相之间进行验证并对计算错误的样本进行过滤,整体资源消耗极低,可以在边缘设备运行。
Description
技术领域
本申请涉及一种基于音频的转速识别方法,适用于转速识别的技术领域。
背景技术
在旋转机械中,对于转速的识别对于生产和维护来说具有重要的意义。例如,对于风机叶片来说,无论是从风机安全运行方面还是从发电综合性能上来看,叶片转速都是风电机组最基本也是较为重要的控制量,从其他角度来看,出于对叶片运行状况的监控,叶片转速可以为后续的数据分析算法提供可靠输入,是重要的参考量。但是现实中,出于硬件条件或安全性、隐私性考虑,并不一定可以实时通过SCADA(数据采集与监视控制)系统采集到叶片实时转速,此时需要通过其他方法计算得到近似的实时转速,目前自动计算叶片实时转速的方法还是空白。
公开号为CN110658519A的中国专利中,通过对外辐射源雷达接收到的直达波信号和雷达回波信号进行预处理,获取距离谱数据,再根据自相关谱的谱峰时间间隔估计叶片转速。该专利中,计算叶片转速为中间步骤,在实施时需要对风电机组加装外辐射源雷达,根据两类波形信号进而计算谱数据,在工程实施上较为复杂,在多风机条件下如果外辐射源雷达系统参考通道中出现多径干扰,会影响谱数据的计算精度进而影响叶片转速的计算准确度。
公开号为CN110470265A的中国专利中,基于微波传感器对动叶片叶尖间隙与定时值的实时在线测量,可计算获得动叶片叶尖间隙、振动及应力等关键参数,进而通过相位峰值与转速计算模块计算相位峰值和叶片转速。该专利中需要加装精密的传感设备,具有一定的复杂度。
上述方法中,都需要加装额外的传感设备,具有一定的复杂度和高成本。计算转速为中间结果,难以为后续分析过程提供精准输入。
发明内容
本申请的目的是提供一种基于音频的转速识别方法,其能够根据声音计算转速,能够脱离对数据采集和监控系统的依赖,可独立运行、效率高、运算准确,可以和后续算法配合,为后续算法输入转速信息。
本申请涉及一种基于音频的转速识别方法,包括以下步骤:
(1)通过音频采集设备采集音频数据;
(2)将采集到的音频数据转换为时频域矩阵,再转换为对数能量谱;
(3)通过滤波截取中频段的信号,将此频段范围内的对数能量谱按时域方向求均值得到中频能量向量;
(4)通过中值滤波对得到的中频能量向量进行降噪处理,选定零值位置并计算过零次数,根据过零次数估算转速;
(5)通过时间序列分析得到信号的自相关函数,通过定位自相关函数的最高点来确定信号的周期,根据周期来估算转速;
(6)计算功率谱密度,通过定位功率谱密度的峰值来确定信号的周期,根据周期来估算转速;
(7)将上述步骤(4)、步骤(5)和步骤(6)中估算得到的转速求平均值,作为最终确定的转速数值。
其中,在步骤(4)中,以中频能量向量均值为起点,上下一个标准差为波动范围定位零值位置;在步骤(5)中,对自相关函数值的峰值定位方法为:通过识别信号中的任一样本点,其直接的临近点取值均低于该点,定位到自相关函数中所有的极大值点;在步骤(7)中,对根据过零次数、自相关函数和功率谱密度确定的转速分别计算相对误差,若任两个计算得到的相对误差超过设定误差阈值,则不输出计算结果。
本申请仅需加装音频采集装置,就可以通过采集声音信号,记录风机叶片扫风声音,然后通过至少两种算法迅速计算当前叶片的实时转速,互相之间进行验证并对计算错误的样本进行过滤,整体资源消耗极低,可以在边缘设备运行。本申请还可以在无法接入SCADA系统时,作为补充手段迅速准确计算出风机当前的实时转速,并可以作为信息输入后续需要转速信息的故障检测算法中。
附图说明
图1为本申请实施例中的叶片扫风声音信号的示意图。
图2为本申请实施例中得到的对数能量谱的示意图。
图3为本申请实施例中截取的频率范围。
图4为本申请实施例中得到的中频能量特征。
图5为本申请实施例中对中频能量特征向量进行平滑去噪后确定的零值位置。
图6为本申请实施例中对滤波后信号计算滞后阶数的自相关函数的示意图。
图7为本申请实施例中根据极值点坐标得到滤波后信号周期的示意图。
图8为本申请实施例中对滤波后信号计算功率谱函数确定功率谱密度最大值及其坐标。
图9为本申请的转速识别方法的步骤示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚明白,下文中将结合附图对本申请的实施例进行详细说明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互任意组合。
参见图9,根据本申请的一种基于音频的转速识别方法,包括以下步骤:
(1)通过音频采集设备采集音频数据;
(2)然后将采集到的音频数据采用例如短时傅里叶变换的方法转换为时频域矩阵,再转换为对数能量谱,此时矩阵内元素为分贝形式表示;
(3)通过滤波截取中频段的信号,将此频段范围内对数能量谱按时域方向求均值得到中频能量向量;这是因为传感器接收风机叶片扫风声音信号为周期性信号,为了避免低频噪声和高频故障特征对信号质量的影响,可以通过带通滤波截取中频段的信号,通常为2048-4096Hz之间,将此频段范围内对数能量谱按时域方向求均值得到中频能量向量;
(4)计算信号的过零次数
首先通过中值滤波对上述中频能量向量进行降噪处理,通过选定合适的邻域大小可以去除部分离群值对算法稳定性的影响。
选定零值位置用于计算过零次数,通常在稳定的扫风情况下,零值位置的定位比较灵活,为了避免叶片间不对称扫风形态的影响,通常以中频能量向量均值为起点,上下一个标准差为波动范围定位零值位置。
本实施例中,设滤波后能量信号为x,x为周期性信号,则零值位置t可表示为:
t=mean(x)+std(x)
对信号x进行去零值处理,减去零值后得到信号x_可表示为:
x_=x-t
过零次数的计算方法为,将序列x_分解为由x_第二个元素至末尾的子序列x_1和由x_起始至倒数第二个元素的子序列x_0,将两个子序列相乘后计算小于零元素的个数,即为信号x_的过零次数。这是因为,如果一个元素相邻两点乘积小于0,即相邻两点分别位于正负位置,说明此处经过了0点。由于叶片扫风信号的特殊性,每一支叶片转动一次会产生两个过零点,对于三支叶片风机来说,此时根据声音计算得到的估计转速rs可表示为:
rs=sum((x1×x0)<0)//6
(5)计算自相关函数
通过时间序列分析得到信号x的自相关函数,其中自相关函数ACF用来度量时间序列中每隔k个时间单位(xt和xt-k)的观测值之间的相关,可使用不同滞后长度观测值的子序列之间的相关系数表示。因此,在选定一定的滞后阶数k后,会计算k个不同滞后长度的观测值子序列与原始序列的相关系数,得到长度为k的ACF序列,一个序列和其不同滞后阶数k的历史观测值之间的相似程度可表示为:
在自相关函数值第一次达到最高点时,代表此时原序列经过了一个扫风周期,因此通过定位自相关函数最高点即可确定原始信号x的周期。对自相关函数值的峰值定位方法为:通过识别信号中的任一样本点,其直接的临近点取值均低于该点定位到自相关函数中所有的极大值点;可以设置0.8为阈值,仅保留超过阈值的最大极大值点,为最终定位的第一次最高点坐标h,如图6所示。h代表着信号x中一支叶片和临近的一支叶片扫风之间的间隔周期,对于三支叶片风机来说此时根据声音计算得到的估计转速rs可表示为:
其中,length(x)表示信号的长度;
(6)计算功率谱密度
同样基于信号x,对其进行离散傅里叶变换(DFT),序列x(n),n=1,2,…,N-1的DFT是一个复数序列X(f):
其中N表示傅里叶变换的点数,k表示傅里叶变换的第k个频谱,k/N表示每个系数捕捉到的频率。基于DFT计算功率谱密度P,可表示为:
通过定位功率谱密度的峰值即可定位到原始信号的周期,功率谱密度中每个元素表示在频率k/N的能量,也即周期N/k的能量,因此功率谱密度峰值所在频率的倒数即为原始信号的周期。
对于三支叶片风机来说,此时根据声音计算得到的估计转速rs可表示为:
其中,length(x)表示信号的长度,argmax(P)表示功率谱密度最大值对应的坐标。
(7)最终确认转速
设根据过零次数、自相关函数和功率谱密度定位的转速分别为rs1,rs2,rs3,则首先对rs1,rs2,rs3计算相对误差,以5%为预先设定的误差阈值,对每两组转速计算误差百分比,若有一次计算误差超过5%则代表此组数据转速计算有误,不输出计算结果。在保证组间相对误差都小于5%时,以rs1,rs2,rs3的均值输出作为最终的估计转速,可表示为:
rs=mean(rs1,rs2,rs3)
本申请针对已有的依据声音进行故障诊断,且依赖实时转速输入的算法,提出了无需接入SCADA系统即可计算转速的方法,能够降低整体系统的复杂度,使故障检测系统能够脱离风机自身数据采集和监控系统独立运行。
实施例
图1所示为通过麦克风采集风机运行时长为1分钟的叶片扫风声音信号。对信号进行短时傅里叶变换,FFT点数选择8192,交叠长度选择7168,得到STFTs谱图矩阵,根据对数变换公式20log10(STFTs)对谱图矩阵计算得到对数能量谱,如图2所示。如图3中虚线所示,截取频率范围在2048Hz-4096Hz的信号,按照时域方向计算对数能量均值,得到如图4所示的中频能量特征。使用中值滤波,邻域长度选择为11,对中频能量特征向量进行平滑去噪。如图5所示,根据滤波后信号x_均值确定零值位置,计算此时信号x_过零次数为66,则根据过零次数得到的叶片转速为66/6=11RPM。对滤波后信号x_计算滞后阶数1000的自相关函数,如图6所示。以0.8为阈值,通过峰值定位算法确定超过最大值0.8的极值点,坐标为86,如图7所示。根据极值点坐标得到滤波后信号周期为86,信号长度为3001,此时得到转速为3001/86//3=11。对滤波后信号x_计算功率谱函数,信号时长为60秒,长度为3001,采样率为len(x_)/60=3001/60≈50Hz,分别得到功率谱密度Pxx_den和对应的频率值f,如图8所示。确定功率谱密度最大值及其坐标,最大值所在坐标为argmax(Pxx_den)=35,滤波后信号x_在第35个元素时达到最大能量,在f中对应的频率为0.583,则信号x_的周期为1/f[argmax(Pxx_den)]*(len(x_)/60)=0.583*50=85.74285714285715,此时叶片转速为3001/85.74285714285715//3=11。由于根据过零次数、自相关函数和功率谱计算得到的转速都为11,则最终转速为11RPM。
虽然本申请所揭露的实施方式如上,但所述的内容只是为了便于理解本申请而采用的实施方式,并非用以限定本申请。任何本申请所属技术领域内的技术人员,在不脱离本申请所揭露的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本申请的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (4)
1.一种基于音频的转速识别方法,其特征在于,包括以下步骤:
(1)通过音频采集设备采集音频数据;
(2)将采集到的音频数据转换为时频域矩阵,再转换为对数能量谱;
(3)通过滤波截取中频段的信号,将此频段范围内的对数能量谱按时域方向求均值得到中频能量向量;
(4)通过中值滤波对得到的中频能量向量进行降噪处理,选定零值位置并计算过零次数,根据过零次数估算转速;
(5)通过时间序列分析得到信号的自相关函数,通过定位自相关函数的最高点来确定信号的周期,根据周期来估算转速;
(6)计算功率谱密度,通过定位功率谱密度的峰值来确定信号的周期,根据周期来估算转速;
(7)将上述步骤(4)、步骤(5)和步骤(6)中估算得到的转速求平均值,作为最终确定的转速数值。
2.根据权利要求1所述的转速识别方法,其特征在于,在步骤(4)中,以中频能量向量均值为起点,上下一个标准差为波动范围定位零值位置。
3.根据权利要求1或2所述的转速识别方法,其特征在于,在步骤(5)中,对自相关函数值的峰值定位方法为:通过识别信号中的任一样本点,其直接的临近点取值均低于该点,定位到自相关函数中所有的极大值点。
4.根据权利要求1-3中任一项所述的转速识别方法,其特征在于,在步骤(7)中,对根据过零次数、自相关函数和功率谱密度确定的转速分别计算相对误差,若任两个计算得到的相对误差超过设定误差阈值,则不输出计算结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011478552.5A CN112782421B (zh) | 2020-12-15 | 2020-12-15 | 一种基于音频的转速识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011478552.5A CN112782421B (zh) | 2020-12-15 | 2020-12-15 | 一种基于音频的转速识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112782421A true CN112782421A (zh) | 2021-05-11 |
CN112782421B CN112782421B (zh) | 2022-06-17 |
Family
ID=75750956
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011478552.5A Active CN112782421B (zh) | 2020-12-15 | 2020-12-15 | 一种基于音频的转速识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112782421B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113738632A (zh) * | 2021-09-07 | 2021-12-03 | 中国商用飞机有限责任公司 | 用于测量液压泵转速的方法和系统 |
CN113915078A (zh) * | 2021-11-26 | 2022-01-11 | 龙源(北京)风电工程技术有限公司 | 光纤光栅式风电机组叶片工况智能识别方法及系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2172819A2 (en) * | 2008-08-27 | 2010-04-07 | China Steel Corporation | Method for evaluating grinding state |
JP2011252715A (ja) * | 2010-05-31 | 2011-12-15 | Rion Co Ltd | エンジン回転数計測方法及びその装置 |
CN102682765A (zh) * | 2012-04-27 | 2012-09-19 | 中咨泰克交通工程集团有限公司 | 高速公路音频车辆检测装置及其方法 |
CN103743922A (zh) * | 2014-01-26 | 2014-04-23 | 深圳科瑞技术股份有限公司 | 一种基于dsp的发动机转速测量方法和装置 |
EP3015866A1 (en) * | 2014-10-30 | 2016-05-04 | ABB Technology AG | A method of automatic determination of rotational speed of a shaft in a rotating machine |
CN107843740A (zh) * | 2017-10-31 | 2018-03-27 | 湖南大学 | 一种融合振动与声音信号频谱特征的转速测量方法 |
CN107918032A (zh) * | 2017-11-14 | 2018-04-17 | 湖南大学 | 一种空间多源声音信号融合的转速测量方法 |
CN109357749A (zh) * | 2018-09-04 | 2019-02-19 | 南京理工大学 | 一种基于dnn算法的电力设备音频信号分析方法 |
CN110792563A (zh) * | 2019-11-04 | 2020-02-14 | 北京天泽智云科技有限公司 | 基于卷积生成对抗网络的风电机组叶片故障音频监测方法 |
CN111161756A (zh) * | 2020-02-13 | 2020-05-15 | 北京天泽智云科技有限公司 | 风机叶片扫风声音信号中异常哨声轮廓的提取及识别方法 |
-
2020
- 2020-12-15 CN CN202011478552.5A patent/CN112782421B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2172819A2 (en) * | 2008-08-27 | 2010-04-07 | China Steel Corporation | Method for evaluating grinding state |
JP2011252715A (ja) * | 2010-05-31 | 2011-12-15 | Rion Co Ltd | エンジン回転数計測方法及びその装置 |
CN102682765A (zh) * | 2012-04-27 | 2012-09-19 | 中咨泰克交通工程集团有限公司 | 高速公路音频车辆检测装置及其方法 |
CN103743922A (zh) * | 2014-01-26 | 2014-04-23 | 深圳科瑞技术股份有限公司 | 一种基于dsp的发动机转速测量方法和装置 |
EP3015866A1 (en) * | 2014-10-30 | 2016-05-04 | ABB Technology AG | A method of automatic determination of rotational speed of a shaft in a rotating machine |
CN107843740A (zh) * | 2017-10-31 | 2018-03-27 | 湖南大学 | 一种融合振动与声音信号频谱特征的转速测量方法 |
CN107918032A (zh) * | 2017-11-14 | 2018-04-17 | 湖南大学 | 一种空间多源声音信号融合的转速测量方法 |
CN109357749A (zh) * | 2018-09-04 | 2019-02-19 | 南京理工大学 | 一种基于dnn算法的电力设备音频信号分析方法 |
CN110792563A (zh) * | 2019-11-04 | 2020-02-14 | 北京天泽智云科技有限公司 | 基于卷积生成对抗网络的风电机组叶片故障音频监测方法 |
CN111161756A (zh) * | 2020-02-13 | 2020-05-15 | 北京天泽智云科技有限公司 | 风机叶片扫风声音信号中异常哨声轮廓的提取及识别方法 |
Non-Patent Citations (2)
Title |
---|
童亦斌等: "MW级风电机组最大功率追踪复合控制策略", 《农业机械学报》, vol. 42, no. 9, 30 September 2011 (2011-09-30), pages 129 - 132 * |
魏毅立等: "音频法测速研究与探索", 《电声技术》, 31 December 2007 (2007-12-31), pages 4 - 7 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113738632A (zh) * | 2021-09-07 | 2021-12-03 | 中国商用飞机有限责任公司 | 用于测量液压泵转速的方法和系统 |
CN113915078A (zh) * | 2021-11-26 | 2022-01-11 | 龙源(北京)风电工程技术有限公司 | 光纤光栅式风电机组叶片工况智能识别方法及系统 |
CN113915078B (zh) * | 2021-11-26 | 2023-06-06 | 龙源(北京)风电工程技术有限公司 | 光纤光栅式风电机组叶片工况智能识别方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112782421B (zh) | 2022-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105547698B (zh) | 滚动轴承的故障诊断方法及装置 | |
JP7199608B2 (ja) | 風力タービン翼を検査するための方法および装置、ならびにそれらの機器および記憶媒体 | |
CN112782421B (zh) | 一种基于音频的转速识别方法 | |
CN110838302B (zh) | 基于信号能量尖峰识别的音频分割方法 | |
CN104502099A (zh) | 齿轮箱瞬变工况特征分量循频提取方法 | |
CN109029999B (zh) | 基于增强调制双谱分析的滚动轴承故障诊断方法 | |
Liu et al. | An online bearing fault diagnosis technique via improved demodulation spectrum analysis under variable speed conditions | |
Zhang et al. | Improved local cepstrum and its applications for gearbox and rolling bearing fault detection | |
CN108398260B (zh) | 基于混合概率方法的齿轮箱瞬时角速度的快速评估方法 | |
RU2705070C1 (ru) | Способ распознавания типа самолёта с турбореактивным двигателем в импульсно-доплеровской радиолокационной станции | |
CN115954017A (zh) | 一种基于hht的发动机小样本声音异常故障识别方法及系统 | |
CN117606782B (zh) | 风力发电机叶片故障检测方法、装置、设备及存储介质 | |
CN113139430A (zh) | 用于故障检测的声信号分析方法、故障检测的方法、装置 | |
CN110376437B (zh) | 一种克服非阶次频率成分干扰的阶次分析方法 | |
EP4254261A1 (en) | Blade fault diagnosis method, apparatus and system, and storage medium | |
CN114994650A (zh) | 一种改进的螺旋桨转速估计方法 | |
CN115683644A (zh) | 航空发动机双源拍振特征识别方法 | |
CN116256158A (zh) | 一种基于深度信号分离的旋转机械瞬时相位自适应提取方法 | |
CN108106717A (zh) | 一种基于声音信号识别机组状态的方法 | |
CN113090471B (zh) | 风力发电机组的塔架净空音频监测系统、方法及装置 | |
CN113887450A (zh) | 基于调制谱特征的水中目标在线筛选方法 | |
CN117725539B (zh) | 一种复杂工况下风机旋转部件的故障特征提取及分析方法 | |
CN112784702A (zh) | 基于自适应旋转同步提取调频变换的信号处理方法 | |
CN115855461B (zh) | 一种风电机组风机叶片噪音异常系统和方法 | |
CN112863541B (zh) | 一种基于聚类和中值收敛的音频切割方法和系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |