首页 / 专利库 / 信号处理 / 盲源分离 / 一种基于稀疏特征相似性的连轧机振动信号分离方法

一种基于稀疏特征相似性的连轧机振动信号分离方法

阅读:433发布:2020-05-14

专利汇可以提供一种基于稀疏特征相似性的连轧机振动信号分离方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于稀疏特征相似性的连 轧机 振动 信号 分离方法。该方法为:在连轧机的旋转部件敏感处表面一个或多个不同 位置 布置 加速 度 传感器 ,采集各振源工作时的振动信号即观测信号;通过基于时 频谱 分割的稀疏特征提取方法,获取各路采集的振动信号在最佳 原子 集合下的稀疏表示;再通过结构相似性对所有稀疏表示提取的原子进行聚类,估算出振源的个数,并采用各聚类中心对应的原子更新稀疏表示结果;最后,通过新的稀疏表示结果对应的聚类系数,估算出混合矩阵,进而计算出各振源的独立信号,实现连轧机振动信号的分离。该方法可有效滤除连轧机在运行过程中的干扰噪声,提高稀疏特征提取的 精度 ,能更精确地估算振源个数及混合矩阵。,下面是一种基于稀疏特征相似性的连轧机振动信号分离方法专利的具体信息内容。

1.一种基于稀疏特征相似性的连轧机振动信号分离方法,其特征在于,首先,在连轧机的旋转部件敏感处表面一个或多个不同位置布置加速传感器,采集各振源工作时的振动信号,即观测信号;然后,通过基于时频谱分割的稀疏特征提取方法,获取各路采集的振动信号在最佳原子集合下的稀疏表示,亦即,通过基于时频谱分割的稀疏分解方法对每一组观测信号进行稀疏分解,提取出各观测信号的稀疏特征,即观测信号在最佳原子集合下的稀疏表示;之后,再通过结构相似性对所有稀疏表示提取的原子进行聚类,估算出振源的个数,并采用各聚类中心对应的原子更新稀疏表示结果;最后,通过新的稀疏表示结果对应的聚类系数,估算出混合矩阵,进而计算出各振源的独立信号,实现连轧机振动信号的分离。
2.如权利要求1所述的基于稀疏特征相似性的连轧机振动信号分离方法,其特征在于,该方法具体包括以下步骤:
步骤一、振动信号的采集
在连轧机组中每一台轧机或其中至少两台轧机的旋转部件敏感处表面一个或多个不同位置,布置一个或多个加速度传感器,采集轧机轧时各振动源的振动信号,即观测信号S(t):
s(t)=[s1(t),s2(t),...,sk(t),...,sK(t)]T
其中,K为观测信号数目,k表示观测点,k=1,2,…,K;[]T表示矩阵的转置;
sk(t)为t时刻在第k个观测点采集到的振动信号,即观测信号;
步骤二、基于时频谱分割的稀疏分解方法的特征提取
利用基于时频谱分割的稀疏分解方法对每一组观测信号进行稀疏分解,提取出各观测信号的稀疏特征,即观测信号在最佳原子集合下的稀疏表示;具体包括多分辨率广义S变换、计算最优二值化时频谱、时频谱分割、稀疏表示四个处理过程;
步骤三、基于原子相似性的振源个数估计
将所有稀疏表示提取的原子以尺度因子为横坐标、频率因子为纵坐标做出原子分布图,然后根据K-均值聚类方法对其进行分类,聚类的个数即为振源个数;
步骤四、估算混合矩阵,分离振动信号中的各振源独立信号
根据上述步骤三的原子聚类信息,估算出混合矩阵,并分离出振动信号中的各振源独立信号。
3.如权利要求2所述的基于稀疏特征相似性的连轧机振动信号分离方法,其特征在于,步骤一中所述的旋转部件敏感处包括轧机的上支撑轴承座、下支撑辊轴承座、上工作辊轴承座、下工作辊轴承座,任选其中的一个或几个。
4.如权利要求2所述的基于稀疏特征相似性的连轧机振动信号分离方法,其特征在于,步骤二中,
(1)计算各观测信号不同尺度下的多分辨率广义S变换结果:
其中,
MGST——多分辨率广义S变换;
sk(t)——t时刻在第k个观测点采集到的振动信号,即观测信号;
u——位移因子;
f——频率因子;
λ——尺度因子。
5.如权利要求2所述的基于稀疏特征相似性的连轧机振动信号分离方法,其特征在于,步骤二中,
(2)计算最优二值化时频谱:
选择合适的阈值,将多分辨率广义S变换得到的时频谱转换成二值化时频谱:
其中,
T——二值化阈值,且T=2||aλ||2/M
aλ——尺度因子为λ的时频谱对应的数据矩阵;
M——矩阵aλ的数据总个数;
m——矩阵aλ的行索引值;
n——矩阵aλ的列索引值;
将所有分辨率下的二值化时频谱通过并集融合规则进行融合,计算出一个最优的二值化时频谱:
其中,
a——最优二值化时频谱对应的数据矩阵;
Π——所有元素的交运算;
λ1——第1个尺度因子;
λp——第p个尺度因子。
6.如权利要求2所述的基于稀疏特征相似性的连轧机振动信号分离方法,其特征在于,步骤二中,
(3)时频谱分割:
采用像素点8连通的原则,将最优二值化时频谱分割成若干个独立的连通域,并对每个连通域进行编号标注。
7.如权利要求2所述的基于稀疏特征相似性的连轧机振动信号分离方法,其特征在于,步骤二中,
(4)稀疏表示:
计算不同尺度下的广义S变换时频谱在每个标注后的连通域中的能量最大值,并根据能量最大值对应的时频因子提取各连通域的最佳表示原子,构成最佳原子集合,计算出观测信号sk(t)在最佳原子集合下的稀疏表示:
其中:
——振动信号sk(t)在最佳原子集合下的稀疏表示;
L——提取的最佳原子的个数;
gl——提取的第l个最佳原子,且
<sk(t),gl>——sk(t)与gl的內积运算;
——gl对应的时频因子。
8.如权利要求2所述的基于稀疏特征相似性的连轧机振动信号分离方法,其特征在于,步骤四具体包括更新稀疏表示、估算混合矩阵、信号盲源分离三个处理过程:
(1)更新稀疏表示:
将每个聚类下的稀疏表示原子用该聚类的聚类中心所对应的原子替代,并重新计算各通道信号在新的原子集合下的稀疏表示;
(2)估算混合矩阵:
计算每个通道各聚类信号对应的稀疏表示的幅值的平均值即聚类系数,并按照同一顺序排列,得到一个由聚类系数组成的矩阵,该矩阵描述了各振源信号在各通道的能量衰减情况,将矩阵每一列均进行归一化处理,计算出混合矩阵A;
(3)信号盲源分离:
根据计算出的混合矩阵,通过X(t)=A-1S(t)将源信号分离出来;
其中,X(t)为源信号;S(t)为观测信号。

说明书全文

一种基于稀疏特征相似性的连轧机振动信号分离方法

技术领域

[0001] 本发明属于旋转机械状态监测与故障诊断技术领域,涉及一种多机架连轧机振动信号分离方法,尤其涉及一种基于稀疏特征相似性的连轧机振动信号分离方法。

背景技术

[0002] 冶金企业在带生产过程中,粗轧机组和精轧机组是整个轧制过程的核心设备,这些连轧机组中的轴承轧辊的工作状态将直接决定轧制产品的质量精度和性能指标。连轧机机组在长期连续、高强度工作中,经常会出现轧辊表层裂纹、轧辊表层局部剥落、轧辊支撑轴承失效等故障,进而影响产品质量精度,甚至导致设备停产及重大人员伤亡等事故,因此,连轧机的状态监测和故障诊断对保证生产稳定性、提高产品质量、降低生产成本具有重要意义。
[0003] 由于连轧机组是由多台连轧机组成,产生的机械振动会在轧机间进行传递,因此,通过采集振动信号对连轧机进行状态监测和故障诊断时,所采集的振动数据均为不同部件振动信号强耦合的混合数据,这样就会导致故障特征提取和故障定位存在极大的难度,特别是当多台机架同时出现故障时,不采用有效的信号分离方法,无法确定具体的故障点位置
[0004] 信号的盲源分离技术就是将实际检测的信号看成是各通道源信号与混合矩阵混合后的信号,盲源分离的关键在于准确地估算出混合矩阵,进而根据检测的信号逆推出各通道源信号。
[0005] 混合矩阵的估算方法比较常用的就是独立成分分析方法。文献检索查到相关专利:2018年3月20日公开的申请号为201711055592.7的发明专利《一种多振源机械设备故障信号分离装置及方法》,提供一种多振源机械设备故障信号分离装置及方法,将信号通过去均值、白化处理和独立成分分析三步去除各路振动信号之间的相互影响。但是,上述专利存在一个缺陷:在工程应用中,各机械振源之间往往会存在交叉频率成分,导致无法满足独立成分分析算法的独立性假设,从而造成算法失效。
[0006] 为了克服这一问题,文献检索查到相关论文:于刚,周以齐,张卫芬.一种新的相关机械振源盲分离方法[J].振动与冲击,2016,35(15):216-221.提出了一种新的相关机械振源盲分离方法,其中指出:当信号具有较好的稀疏性时,观测信号在某一时频点只有一个源的可能性增大。对于隶属同一源信号的时频点将会正比于源信号在测量信号中的混合衰减系数,因此信号的稀疏特性可以被用于估计源信号的混合矩阵。但是,该方法存在不足之处:通过时频散点图拟合各振源聚类直线并进而计算混合矩阵时,往往会因为噪声的干扰,对于直线方向的估计准确性和鲁棒性不好,无法达到对混合矩阵的良好估计。

发明内容

[0007] 本发明的目的在于,克服现有技术的不足,解决现有技术中存在的问题,提供一种基于稀疏特征相似性的连轧机振动信号分离方法,以提高盲源分离时混合矩阵的计算精度,进而达到连轧机多振源振动信号的准确分离目的。
[0008] 本发明的目的是通过以下技术方案实现的:
[0009] 本发明针对现有技术中存在的以上问题,提出一种新的基于稀疏特征相似性的连轧机振动信号分离方法:首先,在连轧机的旋转部件敏感处表面一个或多个不同位置布置加速传感器,采集各振源(振动源)工作时的振动信号,即观测信号;然后,通过基于时频谱分割的稀疏特征提取方法,获取各路采集的振动信号在最佳原子集合下的稀疏表示,亦即,通过基于时频谱分割的稀疏分解方法对每一组观测信号进行稀疏分解,提取出各观测信号的稀疏特征,即观测信号在最佳原子集合下的稀疏表示;之后,再通过结构相似性对所有稀疏表示提取的原子进行聚类,估算出振源的个数,并采用各聚类中心对应的原子更新稀疏表示结果;最后,通过新的稀疏表示结果对应的聚类系数,估算出混合矩阵,进而计算出各振源的独立信号,实现连轧机振动信号的分离。
[0010] 本发明的基于稀疏特征相似性的连轧机振动信号分离方法,具体包括以下步骤:
[0011] 步骤一、振动信号的采集
[0012] 在连轧机组中每一台轧机或其中至少两台轧机的旋转部件敏感处(该旋转部件敏感处包括轧机的上支撑辊轴承座、下支撑辊轴承座、上工作辊轴承座、下工作辊轴承座,任选其中的一个或几个)表面一个或多个不同位置,布置一个或多个加速度传感器,采集轧机轧钢时各振动源的振动信号,即观测信号S(t):
[0013] s(t)=[s1(t),s2(t),...,sk(t),...,sK(t)]T
[0014] 其中,K为观测信号数目,k表示观测点,k=1,2,…,K;[]T表示矩阵的转置;
[0015] sk(t)为t时刻在第k个观测点采集到的振动信号,即观测信号。
[0016] 步骤二、基于时频谱分割的稀疏分解方法的特征提取
[0017] 利用基于时频谱分割的稀疏分解方法对每一组观测信号进行稀疏分解,提取出各观测信号的稀疏特征,即观测信号在最佳原子集合下的稀疏表示;
[0018] 具体包括多分辨率广义S变换、计算最优二值化时频谱、时频谱分割、稀疏表示四个处理过程:
[0019] (1)计算各观测信号不同尺度下的多分辨率广义S变换结果:
[0020]
[0021] 其中,
[0022] MGST——多分辨率广义S变换;
[0023] sk(t)——t时刻在第k个观测点采集到的振动信号,即观测信号;
[0024] u——位移因子;
[0025] f——频率因子;
[0026] λ——尺度因子。
[0027] (2)计算最优二值化时频谱:
[0028] 选择合适的阈值,将多分辨率广义S变换得到的时频谱转换成二值化时频谱:
[0029]
[0030] 其中,
[0031] T——二值化阈值,且T=2||aλ||2/M
[0032] aλ——尺度因子为λ的时频谱对应的数据矩阵;
[0033] M——矩阵aλ的数据总个数;
[0034] m——矩阵aλ的行索引值;
[0035] n——矩阵aλ的列索引值。
[0036] 将所有分辨率下的二值化时频谱通过并集融合规则进行融合,计算出一个最优的二值化时频谱:
[0037]
[0038] 其中,
[0039] a——最优二值化时频谱对应的数据矩阵;
[0040] Π——所有元素的交运算;
[0041] λ1——第1个尺度因子;
[0042] λp——第p个尺度因子。
[0043] (3)时频谱分割:
[0044] 采用像素点8连通的原则,将最优二值化时频谱分割成若干个独立的连通域,并对每个连通域进行编号标注。
[0045] (4)稀疏表示:
[0046] 计算不同尺度下的广义S变换时频谱在每个标注后的连通域中的能量最大值,并根据能量最大值对应的时频因子提取各连通域的最佳表示原子,构成最佳原子集合,计算出观测信号sk(t)在最佳原子集合下的稀疏表示:
[0047]
[0048] 其中:
[0049] ——振动信号sk(t)在最佳原子集合下的稀疏表示;
[0050] L——提取的最佳原子的个数;
[0051] gl——提取的第l个最佳原子,且
[0052] <sk(t),gl>——sk(t)与gl的內积运算;
[0053] ——gl对应的时频因子。
[0054] 步骤三、基于原子相似性的振源个数估计
[0055] 不同结构的设备其振动信号振荡衰减特性不同,所产生的冲击信号的结构不同,特别是尺度因子和频率因子不同,即只有同一振源产生的冲击信号具有相似性。因此,将所有稀疏表示提取的原子以尺度因子为横坐标、频率因子为纵坐标做出原子分布图,然后根据K-均值聚类方法对其进行分类,聚类的个数即为振源个数。
[0056] 步骤四、估算混合矩阵,分离振动信号中的各振源独立信号
[0057] 根据上述步骤三的原子聚类信息,估算出混合矩阵,并分离出振动信号中的各振源独立信号;具体包括更新稀疏表示、估算混合矩阵、信号盲源分离三个处理过程:
[0058] (1)更新稀疏表示:
[0059] 将每个聚类下的稀疏表示原子用该聚类的聚类中心所对应的原子替代,并重新计算各通道信号在新的原子集合下的稀疏表示。
[0060] (2)估算混合矩阵:
[0061] 混合信号中的各通道信号均包含各振源信号,即各聚类信号。计算每个通道各聚类信号对应的稀疏表示的幅值的平均值(即聚类系数),并按照同一顺序排列,得到一个由聚类系数组成的矩阵,该矩阵描述了各振源信号在各通道的能量衰减情况,将矩阵每一列均进行归一化处理,计算出混合矩阵A。
[0062] (3)信号盲源分离:
[0063] 根据计算出的混合矩阵,通过X(t)=A-1S(t)即可将源信号分离出来;
[0064] 其中,X(t)为源信号;S(t)为观测信号。
[0065] 本发明的有益效果:
[0066] 本发明提供了一种基于稀疏特征相似性的连轧机振动信号分离方法,其优点是:①基于时频谱分割的稀疏分解方法可以有效滤除连轧机在运行过程中的干扰噪声,提高稀疏特征提取的精度;②通过稀疏特征相似性的策略对提取原子进行聚类,可以更加精确地估算振源个数及混合矩阵。由于原子结构仅与尺度因子和频率因子相关,与盲源个数无关,因此,适用于多机架连轧机振动信号的盲源分离。
[0067] 随着用户对连轧机设备稳定性和安全性的日益提高,连轧机的状态监测和故障诊断显得越来越重要,提前发现连轧机关键部位的早期故障并提供可靠的解决方案是国内外各大钢厂评价连轧机设备预知维修可靠度的主要指标。本发明的技术主要围绕如何提高连轧机多源振动信号的分离精度,从而保证后续信号分析和故障诊断的精度,具有较广泛的推广价值。附图说明
[0068] 图1是本发明一种基于时频谱分割的变速旋转机械故障诊断方法的流程图
[0069] 图2是2个通道(2个观测点)的观测信号图;
[0070] 图3是第一组观测信号的稀疏表示结果图;
[0071] 图4是第二组观测信号的稀疏表示结果图;
[0072] 图5是根据提取原子的结构信息产生的散点图及聚类中心;
[0073] 图6是振源1的振动信号;图7是振源2的振动信号;
[0074] 图8是振源3的振动信号;图9是振源4的振动信号。

具体实施方式

[0075] 以下结合附图和具体实施例对本发明作进一步的说明。
[0076] 实施例1
[0077] 通过加速度传感器采集4机架连轧机组(包括4台连轧机)的振动信号,并应用本发明一种基于稀疏特征相似性的连轧机振动信号分离方法,对振动信号进行有效分离。
[0078] 如图1所示,具体方法步骤如下:
[0079] (1)通过加速度传感器采集连轧机组的振动信号:
[0080] 分别在连轧机组其中2台连轧机F1和F2的下工作辊轴承座(2个观测点)布置加速度传感器,采集轧机轧钢时各振动源的振动信号;图2所示为分别在上述2个观测点采集到的2个通道的观测信号(两组观测信号),其中,S1为第1个观测点采集到的振动信号,S2为第2个观测点采集到的振动信号。
[0081] (2)利用基于时频谱分割的稀疏分解方法对每一组观测信号进行稀疏分解,提取出各观测信号的稀疏特征(即观测信号在最佳原子集合下的稀疏表示):
[0082] 具体包括多分辨率广义S变换、计算最优二值化时频谱、时频谱分割、稀疏表示四个操作:
[0083] (1)计算各观测信号不同尺度下的多分辨率广义S变换结果:
[0084]
[0085] 其中,
[0086] MGST——多分辨率广义S变换;
[0087] sk(t)——t时刻在第k个观测点采集到的振动信号,即观测信号;
[0088] u——位移因子;
[0089] f——频率因子;
[0090] λ——尺度因子。
[0091] (2)计算最优二值化时频谱:
[0092] 选择合适的阈值,将多分辨率广义S变换得到的时频谱转换成二值化时频谱:
[0093]
[0094] 其中,
[0095] T——二值化阈值,且T=2||aλ||2/M
[0096] aλ——尺度因子为λ的时频谱对应的数据矩阵;
[0097] M——矩阵aλ的数据总个数;
[0098] m——矩阵aλ的行索引值;
[0099] n——矩阵aλ的列索引值。
[0100] 将所有分辨率下的二值化时频谱通过并集融合规则进行融合,计算出一个最优的二值化时频谱:
[0101]
[0102] 其中,
[0103] a——最优二值化时频谱对应的数据矩阵;
[0104] Π——所有元素的交运算;
[0105] λ1——第1个尺度因子;
[0106] λp——第p个尺度因子。
[0107] (3)时频谱分割:
[0108] 采用像素点8连通的原则,将最优二值化时频谱分割成若干个独立的连通域,并对每个连通域进行编号标注。
[0109] (4)稀疏表示:
[0110] 计算不同尺度下的广义S变换时频谱在每个标注后的连通域中的能量最大值,并根据能量最大值对应的时频因子提取各连通域的最佳表示原子,构成最佳原子集合,计算出观测信号sk(t)在最佳原子集合下的稀疏表示:
[0111]
[0112] 其中:
[0113] ——振动信号sk(t)在最佳原子集合下的稀疏表示;
[0114] L——提取的最佳原子的个数;
[0115] gl——提取的第l个最佳原子,且
[0116] <sk(t),gl>——sk(t)与gl的內积运算;
[0117] ——gl对应的时频因子。
[0118] 按上述方法,分别对S1、S2两组观测信号进行稀疏分解,分别提取出每一组观测信号的稀疏特征。图3所示为第一组观测信号的稀疏表示结果,图4所示为第二组观测信号的稀疏表示结果,图中的S1为第1个观测点采集到的振动信号(观测信号),S2为第2个观测点采集到的振动信号(观测信号),图中的稀疏表示,即稀疏表示信号;图中的残差,即残差信号;残差信号=观测信号-稀疏表示信号。
[0119] 步骤三、基于原子相似性的振源个数估计
[0120] 不同结构的设备其振动信号振荡衰减特性不同,所产生的冲击信号的结构不同,特别是尺度因子和频率因子不同,即只有同一振源产生的冲击信号具有相似性。因此,将所有稀疏表示提取的原子以尺度因子为横坐标、频率因子为纵坐标做出原子分布图,然后根据K-均值聚类方法对其进行分类,图5所示为根据提取原子的结构信息产生的散点图及聚类中心;聚类的个数即为振源个数,从图5可看出,该4机架连轧机组有4个振源。
[0121] 步骤四、估算混合矩阵,分离振动信号中的各振源独立信号
[0122] 根据上述步骤三的原子聚类信息,估算出混合矩阵,并分离出振动信号中的各振源独立信号;具体包括更新稀疏表示、估算混合矩阵、信号盲源分离三个处理过程:
[0123] (1)更新稀疏表示:
[0124] 将每个聚类下的稀疏表示原子用该聚类的聚类中心所对应的原子替代,并重新计算各通道信号在新的原子集合下的稀疏表示。
[0125] (2)估算混合矩阵:
[0126] 混合信号中的各通道信号均包含各振源信号,即各聚类信号。计算每个通道各聚类信号对应的稀疏表示的幅值的平均值(即聚类系数),并按照同一顺序排列,得到一个由聚类系数组成的矩阵,该矩阵描述了各振源信号在各通道的能量衰减情况,将矩阵每一列均进行归一化处理,计算出混合矩阵A。
[0127]
[0128] (3)信号盲源分离:
[0129] 根据计算出的混合矩阵,通过X(t)=A-1S(t)将四个源信号即四个振源的振动信号分离出来;其中,X(t)为源信号;S(t)为观测信号。
[0130] 图6所示为振源1的振动信号,图7所示为振源2的振动信号,图8所示为振源3的振动信号,图9所示为振源4的振动信号。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈