一种断路器操动机构振动特征提取方法

1.本发明涉及电气设备故障诊断技术领域，尤其是一种断路器操动机构振动特征提取方法。


背景技术：

2.高压断路器作为动作型控制和保护设备，其可靠性是保证电力系统安全稳定运行的基础。一旦发生故障，首先影响其保护的电缆线路、电器设备运行安全，其次会造成电网事故不断扩大，给生产、生活带来巨大影响。国际大电网会议调查表明，断路器运行中操动机构故障占比61％。振动信号蕴含大量的机械部件状态信息，利用非侵入式振动信号辨识断路器机械故障取得良好的实践效果。经典振动特征提取方法包括短时能量法、包络分析、功率谱和希尔伯特
‑
黄变换等。断路器动作产生的振动波通过透射、反射和衍射等传输过程呈现出非平稳、非线性特点，在时域和频域均有繁杂的结构特征。传统方法在精准描述信号细节特征和复杂故障辨识方面受到限制。
3.对振动信号进行分解处理，是将复杂信号化繁为简的有效方法。ceemdan(自适应噪声集合经验模态分解complete ensemble empirical mode decomposition with adaptive noise)是在emd(经验模态分解empirical mode decomposition)基础上改进的一种自适应分解算法，能有效解决模态混叠问题。拆解后信号的局部变化可以较好反映细节特征。目前的特征提取方法是在量化细节特征方面缺乏针对性，难以刻画非平稳信号的波形突变、功率分布和频谱形态等关键特征指标，无法充分利用细节变化蕴含的机械部件信息。


技术实现要素：

4.本发明要解决的技术问题是提供一种断路器操动机构振动特征提取方法，能够解决现有技术的不足，提高了对断路器操动故障识别准确率。
5.为解决上述技术问题，本发明所采取的技术方案如下。
6.一种断路器操动机构振动特征提取方法，包括以下步骤：
7.a、采集断路器操动过程的振动信号；
8.b、对振动信号进行分解得到不同频段的模态分量信号imf，并求取各imf的功率谱；
9.c、计算故障敏感因子λ
k
，利用故障敏感因子λ
k
对模态分量进行筛选，得到有效模态分量；
10.d、在极坐标下，根据极径尺度ω和极角尺度b，将极坐标以极点为中心、呈辐射状划分为若干等面积子区域；
11.e、将笛卡尔坐标下功率谱波形的频率和幅值(x，y)经发散因子δ在极坐标(r，θ)中发散；
12.f、以波形散布在子区域频数构建信息熵的概率函数，定义为极坐标下感知波形变
异和功率主峰分布的谱形状熵特征。
13.作为优选，步骤a中，用ni usb
‑
4431采集卡和ad50s压电加速度传感器采集断路器振动信号，ad50s压电加速度传感器采用磁吸方式固定在断路器箱体，以线圈电流信号作为ni usb
‑
4431采集卡的采集触发源。
14.作为优选，步骤b中，计算模态分量信号imf和imf的功率谱包括以下步骤，
15.b1、设e
j
(
·
)为emd分解得到的第j个imf分量，ceemdan分解得到的第k个模态分量记为imf
k
(t)，v
i
(t)表示第i次试验中添加具有标准正态分布的白噪声序列(i＝1，2，
…
i)，ε
j
‑1代表第j个模态分量信噪比系数；
16.b2、原始信号s(t)添加噪声分量构成s(t)+ε0v
i
(t)进行i次试验，ceeemdan分解得到第一阶模态分量imf1(t)及第一阶段的残余信号r
i
(t)为，
[0017][0018]
r1(t)＝s(t)
‑
imf1(t)；
[0019]
b3、由emd对信号r1(t)+ε1e1(v
i
(t))进行i次试验，得到第二阶模态分量imf2(t)为，
[0020][0021]
b4、重复步骤b3，得第i个残余信号r
k
(t)和第i+1阶模态分量imf
k+1
(t)为，
[0022]
r
k
(t)＝r
k
‑1(t)
‑
imf
k
(t)
[0023][0024]
b5、重复步骤b4，直至残余信号极点个数不超过2个时停止分解，最终残余信号r(t)为，
[0025][0026]
对各阶模态分量求功率谱得到pimf
k
(f)，
[0027][0028]
其中，表示取复共轭，m
k
表示第k阶分量的采样点数。
[0029]
作为优选，步骤c中，α
k
为故障信号与其各imf分量的相关系数，β
k
为该故障imf分量与同工况下正常信号的相关系数(k＝1，2，
…
，n)。
[0030]
作为优选，步骤d中，
[0031][0032]
n＝ω
×
b b∈z
[0033]
a为在极径方向划分子区域的基值，ω表示在极径方向划分的段数，a
ω
为极径方向每段划分的长度，b为等分极角的段数，n为极坐标下划分子区域的个数。
[0034]
作为优选，步骤e中，
[0035][0036]
θ＝δ
×
arctan(y/x)
[0037]
δ为4时原本只存在于0～90
°
极角范围的功率谱波形扩展到0～360
°
范围，极径维持不变，起到发散波形的功能。
[0038]
作为优选，步骤f中，概率函数的计算公式为，
[0039][0040]
其中，p
g
为功率谱波形散布在第g个子区域的频数，n为极坐标下划分子区域的个数。
[0041]
采用上述技术方案所带来的有益效果在于：本发明根据振动信号经ceemdan分解后频域谱线聚集性特点，以感知功率谱波形变化、主峰分布、描述非平稳信号细节特征，提出极坐标下定义的谱形状熵特征提取方法。利用ceemdan将频率结构复杂的非平稳振动逐步拆解并求取各分量功率谱，经故障敏感因子筛选后，由发散因子δ对波形进行发散处理，提高对主峰区域的敏感度，根据信息熵理论中发生概率较小事件所含信息量反而较大的理念，在极坐标下定义感知波形变化、功率分布和频谱形态的谱形状熵特征。
[0042]
断路器操动故障源于机械部件本身或部件间能量传递特性变化，伴随振动信号体现为频率成分、强度及其关键时序发生变化。由于断路器动作过程的瞬时性，振动信号具有非线性、非平稳的特点，其频域谱线呈现出繁杂的结构特征，直接分析难以挖掘细节变化特征。ceemdan作为自适应分解算法，是将复杂信号化繁为简的有效方法。ceemdan具有良好的抗混叠效果，频域谱线聚集性强。随着分解阶数增加，逐步剥离出原信号高频信息，实现各阶pimf包含主峰区域的有效划分。
[0043]
故障多是由正常状态受外环境影响逐步演化形成，除包含故障信息外，仍含有部分正常状态特征。故障敏感因子可用于去除与正常状态相关而保留与故障特征紧密联系的imf分量，突出故障状态本质特征。现有的故障敏感因子方法计算值范围为整个实数，数值分布较为分散。当计算值为负数时，负号不代表任何物理含义，不便于进行数值间的比较。本发明对现有的故障敏感因子进行改进，得到二次表达式，使计算结果均为正数，消除符号影响，便于展现不同故障的显著差异。敏感因子数值越大，表征其对应的imf分量蕴含不同与正常状态的特征越重要，由此可筛选出有效的模态分量。
[0044]
本发明开创性的提出了发散因子，能够将笛卡尔坐标下功率谱波形的频率和幅值在极坐标中进行发散。由于ceemdan分解后各频段的功率谱具有较强的聚集性，多个峰值相互叠加，难以提取进一步的细节变化特征。经过发散处理，功率谱中幅值较小的点被限定在
零极轴附近，幅值较大的点沿零轴逆时针旋转发散，使主峰区域数据的散布范围扩大，使功率谱主峰波形落在同一子区域的频数降低，有利于感知功率谱波形变化、主峰分布、描述的细节特征。谱形状熵数值主要取决于波形分布较为分散子区域中的频数，即功率谱主峰区域。因此，谱形状熵特征对功率谱波形和主峰分布的变化具有较强敏感程度，能够反映细节变化特征。
附图说明
[0045]
图1是谱形状熵特征提取流程图。
[0046]
图2是ceemdan分解时域波形图。
[0047]
图3是ceemdan分解功率谱波形图。
[0048]
图4是正常及三种故障下各阶分量故障敏感因子对比图。
[0049]
图5是正常合闸状态下有效模态分量极坐标功率谱图。
[0050]
图6是弹簧疲劳状态下有效模态分量极坐标功率谱图。
[0051]
图7是铁芯卡涩状态下有效模态分量极坐标功率谱图。
[0052]
图8是轴销脱落状态下有效模态分量极坐标功率谱图。
[0053]
图9是基于谱形状熵特征的四种算法辨识结果。
具体实施方式
[0054]
参照图1，本发明的一个具体实施方式包括以下步骤：
[0055]
a、采集断路器操动过程的振动信号；
[0056]
b、对振动信号进行分解得到不同频段的模态分量信号imf，并求取各imf的功率谱；
[0057]
c、计算故障敏感因子λ
k
，利用故障敏感因子λ
k
对模态分量进行筛选，得到有效模态分量；
[0058]
d、在极坐标下，根据极径尺度ω和极角尺度b，将极坐标以极点为中心、呈辐射状划分为若干等面积子区域；
[0059]
e、将笛卡尔坐标下功率谱波形的频率和幅值(x，y)经发散因子δ在极坐标(r，θ)中发散；
[0060]
f、以波形散布在子区域频数构建信息熵的概率函数，定义为极坐标下感知波形变异和功率主峰分布的谱形状熵特征。
[0061]
步骤a中，用ni usb
‑
4431采集卡和ad50s压电加速度传感器采集断路器振动信号，ad50s压电加速度传感器采用磁吸方式固定在断路器箱体，以线圈电流信号作为ni usb
‑
4431采集卡的采集触发源。
[0062]
步骤b中，计算模态分量信号imf和imf的功率谱包括以下步骤，
[0063]
b1、设e
j
(
·
)为emd分解得到的第j个imf分量，ceemdan分解得到的第k个模态分量记为imf
k
(t)，v
i
(t)表示第i次试验中添加具有标准正态分布的白噪声序列(i＝1，2，
…
i)，ε
j
‑1代表第j个模态分量信噪比系数；
[0064]
b2、原始信号s(t)添加噪声分量构成s(t)+ε0v
i
(t)进行i次试验，ceeemdan分解得到第一阶模态分量imf1(t)及第一阶段的残余信号r1(t)为，
[0065][0066]
r1(t)＝s(t)
‑
imf1(t)；
[0067]
b3、由emd对信号r1(t)+ε1e1(v
i
t))进行i次试验，得到第二阶模态分量imf2(t)为，
[0068][0069]
b4、重复步骤b3，得第i个残余信号r
k
(t)和第i+1阶模态分量imf
k+1
(t)为，
[0070]
r
k
(t)＝r
k
‑1(t)
‑
imf
k
(t)
[0071][0072]
b5、重复步骤b4，直至残余信号极点个数不超过2个时停止分解，最终残余信号r(t)为，
[0073][0074]
对各阶模态分量求功率谱得到pimf
k
(f)，
[0075][0076]
其中，表示取复共轭，m
k
表示第k阶分量的采样点数。
[0077]
以k最大值等于9为例进行分解，得到不同特征尺度的9个模态分量imf和1个残余分量，如图2所示。计算各分量功率谱，得到pimf波形如图3所示。
[0078]
步骤c中，α
k
为故障信号与其各imf分量的相关系数，β
k
为该故障imf分量与同工况下正常信号的相关系数(k＝1，2，
…
，n)。
[0079]
各分量故障敏感因子如图4所示，根据图4，选取第1和第3阶imf分量为特征分析对象。
[0080]
步骤d中，
[0081][0082]
n＝ω
×
b b∈z
[0083]
a为在极径方向划分子区域的基值，ω表示在极径方向划分的段数，a
ω
为极径方向每段划分的长度，b为等分极角的段数，n为极坐标下划分子区域的个数。
[0084]
步骤e中，
[0085][0086]
θ＝δ
×
arctan(y/x)
[0087]
δ为4时原本只存在于0～90
°
极角范围的功率谱波形扩展到0～360
°
范围，极径维持不变，起到发散波形的功能。
[0088]
图5
‑
8展示发散处理后各状态类型第1和3阶imf功率谱极坐标波形。功率谱主峰代表该频段信号的主要频率成分，与断路器运行过程中的参与部件以及部件的机械特性相关。
[0089]
步骤f中，概率函数的计算公式为，
[0090][0091]
其中，p
g
为功率谱波形散布在第g个子区域的频数。
[0092]
取正常及三种模拟故障状态的振动信号各50组，提取各组数据的谱形状熵特征后，选用支持向量机(svm)、粒子群优化支持向量机(pso
‑
svm)、反向神经网络(bpnn)、k近邻(knn)四分类方法分别进行状态辨识，结果如图6所示。
[0093]
由图9可知，基于谱形状熵的各种辨识方法均取得较好的结果。谱形状熵能够刻画振动信号波形变异和功率主峰分布特征，对断路器操动故障识别准确率高。
[0094]
在本发明的描述中，需要理解的是，术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
[0095]
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。