基于自搜索特征峰值与局部极差的离心泵故障诊断方法与流程

1.本发明属于离心泵故障诊断技术领域，具体涉及一种基于自搜索特征峰值与局部极差的离心泵故障诊断方法。


背景技术：

2.离心泵以其结构简单、性能稳定、维修方便、价格便宜等优点，在机械、石油化工、国防工业等领域得到了广泛应用。然而离心泵往往在高温、高速的恶劣工况下工作，随着服役时间的增长，离心泵性能会出现明显退化，从而导致故障频发，而离心泵在工程实际中运送的介质又一般具有毒性、腐蚀性、易燃易爆等特点，若无法对离心泵故障进行及时有效地诊断，将会造成环境污染、经济损失甚至人员伤亡等重大事故。因此，对离心泵故障诊断方法进行研究，进而保障其安全平稳运行显得尤为重要。
3.离心泵的常见故障主要包括轴承故障、叶轮故障和气穴。近年来，相关人员对离心泵的故障诊断主要是在提取特征之后，通过诊断专家直接观察、分析所提取的特征，或依靠专业知识在时域、频域和时频域内对信号谱图进行人工分析，人为地对离心泵的健康状态进行判断，从而实现故障诊断。上述方法需要人工判断且过于依赖专家及先验知识，费时费力且易发生状态误判，极大地增大了维护难度，难以实现对离心泵故障的高效准确识别。


技术实现要素：

4.为了克服上述现有技术存在的缺点，本发明的目的在于提出一种基于自搜索特征峰值与局部极差的离心泵故障诊断方法，提高离心泵故障诊断的效率和准确性。
5.为达到上述目的，本发明采取的技术方案如下：
6.基于自搜索特征峰值与局部极差的离心泵故障诊断方法，包括以下步骤：
7.步骤1：获取离心泵的振动信号样本及三相电流信号样本；
8.步骤2：自动搜索离心泵转频，对获取的离心泵振动信号样本做功率谱，在离心泵额定转速内自动搜索功率谱峰值，谱峰值对应的频率即为自动搜索所得的离心泵转频f
search
；
9.步骤3：构建自搜索转频峰值指标rp
v
，首先对获取的离心泵振动信号样本进行低通滤波，得到频率在1280hz以下的信号成分，并计算滤波后振动信号时域峰峰值x
pp
：
10.x
pp
＝max(x)
‑
min(x)
11.其中，x表示低通滤波后的振动信号；
12.然后对低通滤波后的信号进行快速傅里叶变换，在获得的频谱中，对步骤2所得的离心泵转频f
search
左右各2hz频率区间内进行最大幅值自动搜索，得到转频附近最大幅值f
max
，进而得到自搜索转频峰值指标rp
v
：
13.14.步骤4：构建自搜索特征频率峰值指标fp
v
，首先对获取的离心泵振动信号样本进行带通滤波，得到频率在[2500hz,10000hz]区间的信号成分，并计算滤波后振动信号时域峰峰值：
[0015]
x'
pp
＝max(x')
‑
min(x')
[0016]
其中，x'表示带通滤波后的振动信号；
[0017]
利用db10小波对带通滤波后的振动信号进行四层小波分解，并对各层小波分解的高频细节信号进行平方包络，得到各层平方包络谱其中k表示层数，k∈[1,4]；其次将中转频f
search
左右各0.5hz频率区间内的信号幅值取反，得到转频负谱即：
[0018][0019]
通过计算：
[0020][0021]
得到各层消除转频后的平方包络谱根据转频f
search
分别计算离心泵中轴承的外圈、内圈、滚动体、保持架故障特征频率，在中的各故障特征频率附近频段内搜索最高峰值及其对应频率，并输出：
[0022][0023][0024][0025][0026]
其中，分别为外圈、内圈、滚动体、保持架故障特征频率的频段输出，分别为1、2、3层小波分解中的各故障特征频段最高峰值对应频率；最终计算得到自搜索特征频率峰值指标fp
v
：
[0027][0028]
步骤5：构建电流局部极差指标pd
c
，将三相电流信号样本分别分割为长度为1s的小段，每相电流信号包含j段样本，计算分段后三相电流信号每段的峰峰值和最小值其中，j＝1,2,...,j，并依此得到u、v、w相电流局部极差指标：
[0029][0030][0031][0032]
将获得的u、v、w相电流局部极差指标取平均，得到电流局部极差指标pd
c
：
[0033][0034]
步骤6：对步骤3、步骤4、步骤5中构造的指标根据幅值分别设定相应的阈值thh
r
、thh
f
、thh
p
，并根据所构建的指标和相应阈值对离心泵健康状态进行划分：若rp
v
＞thh
r
，则诊断为轴承故障，否则，如果满足fp
v
＞thh
f
，此时判断故障类型为叶轮故障；每组数据均须判断是否有气穴现象，若pd
c
＞thh
p
，则离心泵入口流量不足，存在气穴，若pd
c
≤thh
p
，则离心泵无气穴。
[0035]
本发明的有益效果为：本发明提出一种基于自搜索特征峰值与电流局部极差的离心泵故障诊断方法，首先利用振动信号和电流信号分别构建自搜索转频峰值指标rp
v
、自搜索特征频率峰值指标fp
v
以及电流局部极差指标pd
c
，并根据所构建的不同故障指标分别设定不同阈值，通过所构建指标与阈值的比较实现离心泵故障的自动诊断，克服了传统方法对诊断专家和先验知识的过分依赖，提高了诊断效率和诊断准确性。
附图说明
[0036]
图1为本发明的流程图。
[0037]
图2为实施例离心泵故障的诊断结果。
具体实施方式
[0038]
下面结合实施例和附图对本发明做进一步阐述。
[0039]
如图1所示，基于自搜索特征峰值与局部极差的离心泵故障诊断方法，包括以下步
骤：
[0040]
步骤1：获取离心泵的振动信号样本及三相电流信号样本；
[0041]
步骤2：自动搜索离心泵转频；对获取的离心泵振动信号样本做功率谱，在离心泵额定转速内自动搜索功率谱峰值，谱峰值对应的频率即为自动搜索所得的离心泵转频f
search
；
[0042]
步骤3：构建自搜索转频峰值指标rp
v
；首先对获取的离心泵振动信号样本进行低通滤波，得到频率在1280hz以下的信号成分，并计算滤波后振动信号时域峰峰值x
pp
：
[0043]
x
pp
＝max(x)
‑
min(x)
[0044]
其中，x表示低通滤波后的振动信号；
[0045]
然后对低通滤波后的信号进行快速傅里叶变换，在获得的频谱中，对步骤2所得的离心泵转频f
search
左右各2hz频率区间内进行最大幅值自动搜索，得到转频附近最大幅值f
max
，进而得到自搜索转频峰值指标rp
v
：
[0046][0047]
步骤4：构建自搜索特征频率峰值指标fp
v
；首先对获取的离心泵振动信号样本进行带通滤波，得到频率在[2500hz,10000hz]区间的信号成分，并计算滤波后振动信号时域峰峰值：
[0048]
x'
pp
＝max(x')
‑
min(x')
[0049]
其中，x'表示带通滤波后的振动信号；
[0050]
利用db10小波对带通滤波后的振动信号进行四层小波分解，并对各层小波分解的高频细节信号进行平方包络，得到各层平方包络谱其中k表示层数，k∈[1,4]；其次将中转频f
search
左右各0.5hz频率区间内的信号幅值取反，得到转频负谱即：
[0051][0052]
通过计算：
[0053][0054]
得到各层消除转频后的平方包络谱根据转频f
search
分别计算离心泵中轴承的外圈、内圈、滚动体、保持架故障特征频率，在中的各故障特征频率附近频段内搜索最高峰值及其对应频率，并输出：
[0055]
[0056][0057][0058][0059]
其中，分别为外圈、内圈、滚动体、保持架故障特征频率的频段输出，分别为1、2、3层小波分解中的各故障特征频段最高峰值对应频率；最终计算得到自搜索特征频率峰值指标fp
v
：
[0060][0061]
步骤5：构建电流局部极差指标pd
c
；将三相电流信号样本分别分割为长度约为1s的小段，每相电流信号包含j段样本，计算分段后三相电流信号每段的峰峰值和最小值其中，j＝1,2,...,j，并依此得到u、v、w相电流局部极差指标：
[0062][0063][0064][0065]
将获得的u、v、w相电流局部极差指标取平均，得到电流局部极差指标pd
c
：
[0066][0067]
步骤6：对步骤3、步骤4、步骤5中构造的指标根据幅值分别设定相应的阈值thh
r
、thh
f
、thh
p
，并根据所构建的指标和相应阈值对离心泵健康状态进行划分：若rp
v
＞thh
r
，则诊断为轴承故障，否则，如果满足fp
v
＞thh
f
，此时判断故障类型为叶轮故障；由于气穴可能伴随轴承故障、叶轮故障同时存在，因此每组数据均须判断是否有气穴现象，若pd
c
＞thh
p
，则离心泵入口流量不足，存在气穴，若pd
c
≤thh
p
，则离心泵无气穴。
[0068]
实施例：以应用离心泵故障实验数据为例，验证本发明的可行性。
[0069]
获得的离心泵故障实验数据是分别在30hz、40hz、50hz三种不同转速下测得的离心泵不同故障类型实验数据，其中离心泵故障包括轴承故障、叶轮故障以及气穴，由于气穴可能伴随轴承故障、叶轮故障同时存在，因此故障类型又可进一步细分为轴承故障 气穴、单一轴承故障、叶轮故障 气穴、单一叶轮故障和单一气穴，每组数据均需判断是否有气穴。离心泵轴承型号为6306和6307，包含nsk和skf两种，两类轴承在转频f下的理论故障特征频率分别如表1和表2所示：
[0070]
表1 nsk 6306/skf 6306轴承故障特征频率(转频f)
[0071][0072]
表2 nsk 6307/skf 6307轴承故障特征频率(转频f)
[0073][0074]
实验过程中，分别使用加速度传感器和电流钳获取离心泵的振动信号和三相电流信号。
[0075]
使用本发明方法对离心泵进行故障诊断，首先针对振动信号和三相电流信号，提取本发明提出的自搜索转频峰值指标rp
v
、自搜索特征频率峰值指标fp
v
和电流局部极差指标pd
c
；其次分别为每个指标设定相应的阈值，其中thh
r
设定为0.5，thh
f
设定为1000，thh
p
设定为1000；然后对不同转速、不同故障类型的离心泵进行故障诊断，诊断结果如图2所示，从图2中可以看出，对于不同的离心泵故障类型，针对不同转速下的监测信号，本发明方法能实现对离心泵故障的自动诊断，且诊断准确率高，表明本发明有效地克服了传统方法对诊断专家和先验知识的过分依赖，提高了离心泵故障诊断的效率和性能。
[0076]
在具体应用中，实施者可根据实际情况，通过对本发明中的阈值进行适当调整，便可用本发明方法实现对离心泵的自动故障诊断，有助于提高离心泵故障诊断的效率和性能。应当指出，在不脱离本发明构想的前提下，对本发明方法所做的调整和变形，也应视为本发明的保护范围。