CN115859083A - 一种水电机组中多振源耦合振动信号的解耦方法 - Google Patents
一种水电机组中多振源耦合振动信号的解耦方法 Download PDFInfo
- Publication number
- CN115859083A CN115859083A CN202211558002.3A CN202211558002A CN115859083A CN 115859083 A CN115859083 A CN 115859083A CN 202211558002 A CN202211558002 A CN 202211558002A CN 115859083 A CN115859083 A CN 115859083A
- Authority
- CN
- China
- Prior art keywords
- vibration
- vibration signal
- signal
- generating set
- signals
- 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.)
- Pending
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种水电机组中多振源耦合振动信号的解耦方法,包括如下步骤:根据预设频率采集水电机组水平方向和垂直方向的振动信号并进行预处理,得到实测振动信号;利用构建的时域卷积窗对实测振动信号进行加权后再进行傅里叶变换,得到恒定频率;对幅值和相位表示形式的实测振动信号的瞬时相量进行二阶泰勒展开,得到多振源耦合振动信号模型;求解实测振动信号与该模型的误差方程,得到实测振动信号的瞬时向量即为原振动信号的调制信号,实现振动信号的解耦。本发明通过构建新型窗函数实现振动信号原信号频率准确测量的基础上,基于泰勒展开式建立多振源耦合下的振动信号模型,然后通过求解模型参数以实现对快时变、强调频下振动信号的解耦。
Description
技术领域
本发明属于水电机组振动信号处理领域,具体涉及一种水电机组中多振源耦合振动信号的解耦方法。
背景技术
水电机组运行可靠性对水电站的安全与稳定有着至关重要的影响。然而水电机组的运动状态会受到水力、机械和电气三个方面因素的耦合作用和影响,如机组的旋转机械部分的机械运动会导致机组中的水流和电机气隙间的磁场发生变化;机组中水流流场的变化也会导致机组中电机气隙间磁场的变化,同时也会引起机械的振动。因而水电机组是水-机-电耦合的多源非线性系统,振动信号也会表现出强烈的多振源耦合特性。
现有的振动信号处理方法主要分为时域、频域和时频域。时域方法通常计算信号自身的数学特征或概率分布作为检测依据,如互相关函数、概率密度函数等,但易受噪声干扰;以傅里叶变换为基础的频域方法可以获得振动信号的特征频谱,但缺乏对非平稳信号的时域描述信息;时频域方法可以从时间、频率两个维度刻画信号特征,在处理时变信号方面具有优势。
尽管目前对水电机组非平稳振动信号的故障特征提取方面的研究较多,但解决水电机组多振源耦合振动信号的辨识与分离方面并不多见。
发明内容
本发明要解决的技术问题是针对现有技术的不足而提供水电机组中多振源耦合振动信号的解耦方法,对水电机组多振源耦合振动信号进行辨识与分离。
为解决上述技术问题,本发明的内容包括:
一种水电机组中多振源耦合振动信号的解耦方法,包括如下步骤:
S1、根据预设频率采集水电机组水平方向和垂直方向的振动信号;
S2、对采集的振动信号进行初步去噪处理,得到实测振动信号;
S3、构建时域卷积窗;
S4、利用构建的时域卷积窗对实测振动信号进行加权,并对加权后的实测振动信号进行傅里叶变换,得到恒定频率;
S5、对幅值和相位表示形式的实测振动信号的瞬时相量进行二阶泰勒展开,得到的近似信号作为多振源耦合振动信号模型;
S6、求解实测振动信号与多振源耦合振动信号模型的误差方程,得到实测振动信号的瞬时向量,即为原振动信号的调制信号,从而实现振动信号的解耦。
进一步的,所述步骤S2中,采用滤波器及低秩稀疏分解法对采集的振动信号进行初步去噪处理。
进一步的,所述步骤S2中,实测振动信号用x(n)表示,其中n=0,1,…,N,N为奇数。
进一步的,所述步骤S3中,将长度为M1的三阶最大旁瓣衰减窗和长度为M2的三阶最小旁瓣Nuttall窗进行时域卷积并补零后,得到长度为M的时域卷积窗w(m),其中,M=M1+M2,m=0,1,…,M,M<N。
进一步的,所述步骤S4中,加权后的实测振动信号为xw(m)=w(m)x(m),对加权后的实测振动信号xw(m)进行傅里叶变换后得到恒定频率f0。
进一步的,所述步骤S4中,采用二谱线插值离散傅里叶变换或三谱线插值离散傅里叶变换求解得到恒定频率f0。
进一步的,所述步骤S5中,实测振动信号x(n)的幅值和相位表示形式为:
进一步的,根据最小二乘法求得使多组误差方程en最小时的p(n0)、p′(n0)、p″(n0),得到瞬时向量p(n)。
本发明的有益效果是:
1)构建新型窗函数实现振动信号原信号频率的测量,有效降低了离散频率校正方法中的频谱泄露;2)建立多振源耦合下振动信号模型,用泰勒展开式(即高阶多项式)实现对信号瞬时频率、瞬时幅值的精确估计;3)通过准确求解泰勒展开式系数,建立可读性强的动态信号特征表征方法,完成多振源耦合下振动信号动态特征的映射。
具体实施方式
为便于理解本发明,下面结合具体实施方式对本发明作进一步详细的说明。本领域技术人员应该明了,所述实施例仅仅是帮助理解本发明,不应视为对本发明的具体限制。
本发明提供了一种水电机组中多振源耦合振动信号的解耦方法,包括如下步骤:
S1、根据预设频率fs采集水电机组水平方向和垂直方向的振动信号。
S2、对采集的振动信号进行初步去噪处理,得到实测振动信号x(n),其中n=0,1,…,N, N为奇数。
本步骤中,对于采集到的含噪程度高的振动信号,预先采用滤波器、低秩稀疏分解法等方法进行初步去噪处理。
S3、构建时域卷积窗:将长度为M1的三阶最大旁瓣衰减窗和长度为M2的三阶最小旁瓣 Nuttall窗进行时域卷积并补零后,得到长度为M的时域卷积窗w(m),其中,M=M1+M2, m=0,1,…,M,M<N。
S4、利用构建的时域卷积窗w(m)对实测振动信号x(n)进行加权,并对加权后的实测振动信号xw(m)=w(m)x(m)进行傅里叶变换,得到恒定频率f0。
本步骤中,恒定频率f0的求解可采用二谱线插值离散傅里叶变换、三谱线插值离散傅里叶变换等离散频谱校正方法获得。
本步骤中,实测振动信号x(n)的幅值和相位表示形式为:
S6、根据最小二乘法求得使实测振动信号x(n)与多振源耦合振动信号模型的多组误差方程/>最小时的p(n0)、p′(n0)、p″(n0),得到实测振动信号x(n)的瞬时向量p(n),即为原振动信号的调制信号,从而实现振动信号的解耦。
本发明通过构建新型时域卷积窗函数实现振动信号原信号频率准确测量的基础上,基于泰勒展开式建立多振源耦合下的振动信号模型,然后通过求解模型参数以实现对快时变、强调频下振动信号的解耦。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (9)
1.一种水电机组中多振源耦合振动信号的解耦方法,其特征在于,包括如下步骤:
S1、根据预设频率采集水电机组水平方向和垂直方向的振动信号;
S2、对采集的振动信号进行初步去噪处理,得到实测振动信号;
S3、构建时域卷积窗;
S4、利用构建的时域卷积窗对实测振动信号进行加权,并对加权后的实测振动信号进行傅里叶变换,得到恒定频率;
S5、对幅值和相位表示形式的实测振动信号的瞬时相量进行二阶泰勒展开,得到的近似信号作为多振源耦合振动信号模型;
S6、求解实测振动信号与多振源耦合振动信号模型的误差方程,得到实测振动信号的瞬时向量,即为原振动信号的调制信号,从而实现振动信号的解耦。
2.根据权利要求1所述的水电机组中多振源耦合振动信号的解耦方法,其特征在于:所述步骤S2中,采用滤波器及低秩稀疏分解法对采集的振动信号进行初步去噪处理。
3.根据权利要求1所述的水电机组中多振源耦合振动信号的解耦方法,其特征在于:所述步骤S2中,实测振动信号用x(n)表示,其中n=0,1,…,N,N为奇数。
4.根据权利要求3所述的水电机组中多振源耦合振动信号的解耦方法,其特征在于:所述步骤S3中,将长度为M1的三阶最大旁瓣衰减窗和长度为M2的三阶最小旁瓣Nuttall窗进行时域卷积并补零后,得到长度为M的时域卷积窗w(m),其中,M=M1+M2,m=0,1,…,M,M<N。
5.根据权利要求4所述的水电机组中多振源耦合振动信号的解耦方法,其特征在于:所述步骤S4中,加权后的实测振动信号为xw(m)=w(m)x(m),对加权后的实测振动信号xw(m)进行傅里叶变换后得到恒定频率f0。
6.根据权利要求5所述的水电机组中多振源耦合振动信号的解耦方法,其特征在于:所述步骤S4中,采用二谱线插值离散傅里叶变换或三谱线插值离散傅里叶变换求解得到恒定频率f0。
9.根据权利要求8所述的水电机组中多振源耦合振动信号的解耦方法,其特征在于:根据最小二乘法求得使多组误差方程en最小时的p(n0)、p′(n0)、p″(n0),得到瞬时向量p(n)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211558002.3A CN115859083A (zh) | 2022-12-06 | 2022-12-06 | 一种水电机组中多振源耦合振动信号的解耦方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211558002.3A CN115859083A (zh) | 2022-12-06 | 2022-12-06 | 一种水电机组中多振源耦合振动信号的解耦方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115859083A true CN115859083A (zh) | 2023-03-28 |
Family
ID=85670415
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211558002.3A Pending CN115859083A (zh) | 2022-12-06 | 2022-12-06 | 一种水电机组中多振源耦合振动信号的解耦方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115859083A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117040339A (zh) * | 2023-10-09 | 2023-11-10 | 潍柴动力股份有限公司 | 一种永磁同步电机矢量解耦方法、装置、设备和存储介质 |
-
2022
- 2022-12-06 CN CN202211558002.3A patent/CN115859083A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117040339A (zh) * | 2023-10-09 | 2023-11-10 | 潍柴动力股份有限公司 | 一种永磁同步电机矢量解耦方法、装置、设备和存储介质 |
CN117040339B (zh) * | 2023-10-09 | 2024-01-16 | 潍柴动力股份有限公司 | 一种永磁同步电机矢量解耦方法、装置、设备和存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107451557B (zh) | 基于经验小波变换与局部能量的输电线路短路故障诊断方法 | |
CN103245832B (zh) | 基于快速s变换的谐波时频特性参数估计方法及分析仪 | |
CN108152674B (zh) | 一种基于特征点辨识和线性插值的故障行波滤波方法 | |
CN107247182B (zh) | 一种基于量测相量数据的间谐波分量还原方法 | |
CN110389312B (zh) | 一种适用于现场pmu测试的校准器相量测量方法 | |
CN103576002A (zh) | 一种容性绝缘设备介质损耗角的计算方法 | |
CN110221116B (zh) | 基于加窗插值和解析模式分解的电压闪变包络检测方法 | |
CN108535613B (zh) | 一种基于组合窗函数的电压闪变参数检测方法 | |
CN109541312A (zh) | 一种新能源汇集地区次同步谐波检测方法 | |
CN115859083A (zh) | 一种水电机组中多振源耦合振动信号的解耦方法 | |
Liu et al. | Improved processing of harmonics and interharmonics by time-domain averaging | |
CN105137180A (zh) | 基于六项余弦窗四谱线插值的高精度谐波分析方法 | |
CN109540560B (zh) | 旋转机械结构复谐动态过程的绝对抗混叠多尺度滤波方法 | |
CN109669072B (zh) | 一种配电网的自适应同步相量量测方法 | |
CN102253282A (zh) | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 | |
CN112018784B (zh) | 一种基于同步相量测量数据的次同步谐振溯源方法 | |
CN104502707A (zh) | 一种基于三次样条插值的电力系统同步相量测量方法 | |
CN104714075A (zh) | 一种电网电压闪变包络参数提取方法 | |
CN103941090A (zh) | 基于谱线能量插值的谐波测量方法 | |
Xu et al. | A novel frequency estimation method based on complex bandpass filters for P-class PMUs with short reporting latency | |
CN109992748B (zh) | 一种河口潮位调和分析方法 | |
CN103592514A (zh) | 一种新型谐波高精度检测方法 | |
CN104849552A (zh) | 一种基于自适应陷波滤波器的谐波提取方法 | |
CN104076203B (zh) | 一种考虑负频率影响的超低频间谐波检测方法 | |
CN111523231B (zh) | 一种基于EEMD和Prony方法的次同步振荡分析方法 |
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 |