本发明涉及有载分接开关故障诊断
技术领域：
，具体涉及一种基于itd与svm的有载分接开关oltc机械故障诊断方法。
背景技术：
：随着对电能质量要求的提高，电网大量应用自动电压控制等系统，现有有载分接开关调节相当频繁，故障发生率很高。据国外资料统计，分接开关故障占有载调压变压器故障的40％以上，国内平均统计数据表明，分接开关的故障占变压器故障的20％以上，且主要为机械故障，若不及时发现和处理，其故障会严重破坏oltc和变压器的固有结构，影响电力设备和系统的正常安全运行并造成严重后果。因此，为了确保分接开关安全可靠地运行，有必要开展分接开关机械故障诊断方法的相关研究。作为变压器中唯一可进行机械动作的部件，分接开关操作过程中包含一系列动作事件，这些事件中包含丰富的振动信号，其中触头动作过程中的碰撞、摩擦等都伴有机械振动信号的产生。利用加速度传感器监测设备操作过程中的振动信号，对振动信号进行分析，提取故障特征，是目前比较有效的监测和诊断方法。已有的振动信号分析方法有小波奇异性检测、自组织映射法、经验模态分解(emd)和小波包等。这些方法大多是基于线性理论，然而，研究表明，oltc切换过程中的振动信号表现出明显的非线性行为。因此，将振动信号假设为平稳或分段平稳信号进行时频分析的方法，其分析效果不是很明显。而emd分解是近年来提出的一种信号时频分析方法，因此更适合处理非平稳信号，但是emd分解存在比较严重的端点效应和虚假分量的问题。时间固有尺度分解(intrinsictime-scaledecomposition,itd)是在研究emd和局部均值分解(localmeandecomposition,lmd)的基础上提出的一种自适应时频分析方法，因此，本发明采用时间固有尺度分解itd方法提取oltc机械故障特征，可以克服emd的缺点，再结合支持向量机svm，可以取得较好的oltc机械故障诊断效果。技术实现要素：为解决现有技术中的不足，本发明提供一种基于itd与svm的有载分接开关oltc机械故障诊断方法，其诊断结果确实精度高，结构简单，可操作性强。为了实现上述目标，本发明采用如下技术方案：一种基于itd与svm的oltc机械故障诊断方法，包括以下步骤：1)通过加速度传感器对有载分接开关正常状态下的振动信号、故障状态下的振动信号进行采集，并对振动信号做预处理；2)对预处理后的振动信号进行固有时间尺度分解itd分析，构造特征向量，作为svm的输入；3)对svm进行训练，将步骤2)构造的特征向量输入到svm中，对svm进行训练，将测试数据输入到训练好的svm从而判断oltc的故障模式。前述的一种基于itd与svm的oltc机械故障诊断方法，其特征是：所述预处理具体为：对振动信号进行去噪处理。前述的一种基于itd与svm的oltc机械故障诊断方法，其特征是：所述故障状态包括有载分接开关oltc触头松动以及弹簧性能下降。前述的一种基于itd与svm的oltc机械故障诊断方法，其特征是：所述加速度传感器安装在分接开关的顶端。前述的一种基于itd与svm的oltc机械故障诊断方法，其特征是：所述固有时间尺度分解itd分析，具体步骤包括：①确定原始信号xt(t≥0)的极值x(t)及相应的时刻tk'，k'＝1,2,…,m，m为极值点的个数，并且计算：式中，t表示原始信号的序列时刻，k＝1,2,3,…,m-2，0＜α＜1，x(k)、x(k+1)、x(k+2)分别表示原始信号xt在第k、k+1、k+2个极值点的极值，tk、tk+1和tk+2表示极值x(k)、x(k+1)、x(k+2)分别对应的时刻，lk+1表示各基线控制点；②式(1)中可以得到l2,l3,…,lm-1，两端点l1、lm的值需要估计计算，采用镜像延拓向左右两端各延拓一个极值点，得到左右两端极值点分别为(t0,x(0))、(tm+1,x(m+1))，其中x(0)为向左端延拓的极值，t0为极值x(0)对应的时刻，x(m+1)也为向右端延拓的极值，tm+1为极值x(m+1)对应的时刻，令k分别等于0和m-1，按照式(1)求出l1、lm的值，采用有理样条插值函数拟合所有的l1，l2，l3，...lm，得到基线信号lt；③将基线信号lt从原始信号中分离出来，得到剩余信号，即信号经itd分解，第一次分解的结果记为h1(t)：h1(t)＝xt-lt(2)若h1(t)满足内禀尺度分量isc的条件，输出h1(t)，令isc1＝h1(t)，isc1是第一个分解出来的内禀尺度分量isc；否则，h1(t)作为原始数据，继续重复上述步骤①到步骤③，直到h1s(t)为内禀尺度分量，记isc1＝h1s(t)，h1s(t)为h1(t)重复步骤①～③s次后满足isc条件后h1(t)的值；④将isc1分量从原始信号中分离出来，则得到新的剩余信号u1(t)＝xt-isc1；⑤新的剩余信号u1(t)作为原始数据，重复上述步骤①～④循环n-1次，直到剩余信号un(t)为单调信号或常信号，即为itd分解终止条件，将信号能量添加到分解中作为终止条件判据，即：式中，l的范围为1～原始信号的序列长度；uzn(t)为剩余信号un(t)的均值，xz(t)为x(t)的均值，x(t)为原始信号，ε为阈值；itd分解结束后固有时间尺度分解itd将非平稳信号分解成若干个表征信号特征的内禀尺度分量isc和一个单调信号之和，即：式中，isci为第i个分解出来的内禀尺度分量；⑥采用itd分解得到的前5个分量isc1、isc2、…、isc5，分别计算出它们的能量e1、e2、e3、e4和e5，pm为第m个isc的能量占总能量e中的比重，m＝1，2，…，5，能量为信号各点幅度值平方后的求和，将p1、p2、p3、p4和p5作为有载分接开关机械故障诊断的特征向量。前述的一种基于itd与svm的oltc机械故障诊断方法，其特征是：所述振动数据分成两组，每组包括正常振动信号与故障信号，一组用于svm的训练，另一组用于svm的测试；步骤3)中，先将特征向量归一化处理到[0，1]之间，将归一化后的训练特征向量输入到支持向量机svm中，对支持向量机进行训练，然后将归一化后的测试特征向量输入到训练好的svm，从而判断oltc的故障模式。本发明所达到的有益效果：1、本发明采用的itd算法可以克服经验模态分解(emd)、局部均值分解(lmd)中存在比较严重的端点效应和虚假分量的问题，适合提取oltc机械故障信号的特征量；2、本发明不需要大量数据进行svm的训练，诊断精度更高；3、本发明对oltc的诊断效果明显优于bp网络。附图说明图1是正常信号的itd分解结果图；图2是正常信号的emd分解结果图。具体实施方式下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。一种基于itd与svm的oltc机械故障诊断方法，包括以下步骤：1)通过加速度传感器对有载分接开关正常状态下的振动信号、故障状态下的振动信号进行采集，并对振动信号做预处理；2)对预处理后的振动信号进行固有时间尺度分解itd分析，构造特征向量，作为svm的输入；3)对svm进行训练，将步骤2)构造的特征向量输入到svm中，对svm进行训练，将测试数据输入到训练好的svm从而判断oltc的故障模式。步骤1)中振动传感器采用永磁体吸附在测试点的表面的安装方式，考虑到振动信号的传播介质以及传播过程的阻尼，本发明将振动传感器安装在分接开关的顶端，此位置所拾取的振动信号高频衰减比较少，信号较完整；将采集到的振动数据分成两组，每组都包括正常振动信号与故障信号，一组用于svm的训练，另一组用于svm的测试。故障状态是指有载分接开关oltc触头松动以及弹簧性能下降两种情况。振动信号做预处理包括：对振动信号进行去噪处理。步骤2)中，对预处理后的振动信号进行固有时间尺度分解itd分析，提取特征量，itd具体分解步骤如下：①确定原始信号xt(t≥0)的极值x(t)及相应的时刻tk'，k'＝1,2,…,m，m为极值点的个数，并且计算：式中，t表示原始信号的序列时刻，k＝1,2,3,…,m-2，0＜α＜1，一般α＝0.5，x(k)、x(k+1)、x(k+2)分别表示原始信号xt在第k、k+1、k+2个极值点的极值，tk、tk+1和tk+2表示极值x(k)、x(k+1)、x(k+2)分别对应的时刻，lk+1表示各基线控制点；②式(1)中可以得到l2,l3,…,lm-1，两端点l1、lm的值需要估计计算，采用镜像延拓向左右两端各延拓一个极值点，得到左右两端极值点分别为(t0,x(0))、(tm+1,x(m+1))，其中x(0)为向左端延拓的极值，t0为极值x(0)对应的时刻，x(m+1)为向右端延拓的极值，tm+1为极值x(m+1)对应的时刻，令k分别等于0和m-1，按照式(1)求出l1、lm的值，采用有理样条插值函数拟合所有的l1，l2，l3，...lm，得到基线信号lt；③将基线信号lt从原始信号中分离出来，得到剩余信号，即信号经itd分解第一次分解的结果记为h1(t)：h1(t)＝xt-lt(2)若h1(t)满足内禀尺度分量isc的条件，输出h1(t)，令isc1＝h1(t)，isc1是第一个分解出来的内禀尺度分量isc；否则，h1(t)作为原始数据，继续重复上述步骤①到步骤③，直到h1s(t)为内禀尺度分量，记isc1＝h1s(t)，h1s(t)为h1(t)重复步骤①～③s次后满足内禀尺度分量isc条件后h1(t)的值；在实际应用中，设置一个变量h，当|lt|小于等于h时，则h1(t)满足内禀尺度分量isc条件；④将isc1分量从原始信号中分离出来，则得到新的剩余信号u1(t)＝xt-isc1；⑤新的剩余信号u1(t)作为原始数据，重复上述步骤①～④循环n-1次，直到剩余信号un(t)为单调信号或常信号，剩余信号un(t)为单调信号或常信号即为itd分解终止条件；在实际工程应用中，在对分析结果没有影响的情况下，提高分解速度可减少分解时间，因此itd分解终止条件可适当的放宽，将信号能量添加到分解中作为终止条件判据，即：式中，l的范围为1～原始信号的序列长度；uzn(t)为剩余信号un(t)的均值，xz(t)为x(t)的均值，x(t)为原始信号，ε为阈值，本发明ε取0.001，即剩余信号能量小于原始信号能量的0.1％时，itd分解结束，固有时间尺度分解itd将非平稳信号分解成若干个表征信号特征的内禀尺度分量isc和一个单调信号之和，即：式中，isci为第i个分解出来的内禀尺度分量；⑥有载分接开关处于不同的机械状态下，产生的振动信号能量随频率的分布也不同。当信号经固有时间尺度分解itd后，n-1个isc分量包含的频率成分也不一样，经分析，itd分解得到的内禀尺度分量其主要信息集中在前5个分量，因此，本发明采用itd分解得到isc1、isc2、…、isc5，分别计算出它们的能量e1、e2、e3、e4和e5，pm为第m个isc的能量占总能量e中的比重，这个比重记为能量熵，m＝1，2，…，5，能量为信号各点幅度值平方后的求和。本发明将p1、p2、p3、p4和p5作为有载分接开关机械故障诊断的特征向量。由于经验模态分解过程中常出现虚假分量，因此，本文针对虚假分量问题，比较itd和emd的分解效果。由于虚假分量与原信号的相关性很小，通常采用相关系数来识别虚假分量，但相关系数法存在误识别情况，因此，本文采用k-l散度来识别虚假分量。k-l散度是2个概率分布差异的度量，设信号x＝{x1,x2,…,xn}和y＝{y1,y2,…,yn}的概率分布分别为p(x)、q(x)，则信号x和y的k_l散度d(p,q)为：d(p,q)＝δ(p,q)+δ(q,p)(5)式中，δ(p,q)以及δ(q,p)为信号x和y的k_l距离，n为信号点总数。步骤3)中，先使用matlab自带的数据预处理-归一化函数mapminmax将特征量归一化处理[0，1]之间，将归一化后的训练特征向量输入到svm中，对svm进行训练，然后将归一化后的测试特征向量输入到训练好的svm，从而判断oltc的故障模式。训练特征向量和测试特征向量分别利用采集到的训练信号和测试信号进行itd分解得到。实施例：对cmiii-500-63b-10193w型分接开关模拟实验。此分接开关为三相y连接，最大的分接位置数为19。振动传感器采用分辨率高且抗干扰能力强的lc0151型压电式加速度传感器。实验时振动传感器采用永磁体吸附在测试点的表面的安装方式，把振动传感器安装在分接开关的顶端，通过加速度传感器对oltc正常状态下的振动信号、故障状态(oltc触头松动以及弹簧性能下降)下的振动信号进行采集，将采集到的振动数据分成2组，每组都包括正常振动信号与故障信号，一组用于svm的训练，另一组用于svm的测试；为了体现itd方法的优势，本发明对oltc正常状态下的振动信号分别进行emd、itd分解，比较两者性能。其中itd算法分解得到7个内禀尺度分量，而emd算法分解得到10个分量，其主要分量波形如图1、图2。由图1、图2可以看出(图1中仅展示了itd分解得到的前五个分量，图2仅展示了emd分解得到的前六个分量)，itd分解得到的分量与emd相比，前者毛刺部分较少。由于emd分解过程中常出现虚假分量，因此，本发明针对虚假分量问题，比较itd和emd的分解效果。针对虚假分量问题，itd和emd的分解效果很难从图1、图2中直观看出。由于虚假分量与原信号的相关性很小，通常采用相关系数来识别虚假分量，但相关系数法存在误识别情况，因此，本文采用k-l散度来识别虚假分量。k-l散度值越小，表明该分量与原信号的关系越紧，是真实信号；反之，表明该分量与原始信号的关系越远，是信号的虚假分量，应当去除。真实分量与虚假分量具有数量级差异的散度值，因此，经多次仿真分析，本发明选取0.1为阈值，k-l散度大于0.1时，表明该分量是虚假分量，反之，则说明该分量为真实信号。如表1、表2为原始信号与各分量间的k-l散度值。由表1、表2可知，表1中的后两个为虚假分量，itd的虚假分量占2/7*100＝28.6％，表2中的后4个为虚假分量，emd虚假分量占比40％，由此可知，itd算法改善了emd虚假分量问题。表1原信号与itd各分量间的k-l散度值sci1sci2sci3sci4sci5sci6sci70.00100.00530.00640.01070.07110.37621表2原信号与emd各分量间的k-l散度值imf1imf2imf3imf4imf5imf6imf7imf8imf9imf100.00100.00400.00570.00320.01050.07060.10240.27310.55371注：imf为分接开关振动信号经emd分解后得到的10个分量，分别记为imf1，imf2，...imf10；由于itd分解由于emd分解，因此，本发明使用itd分解得到的各个分量能量熵构造特征向量。由表1、表2分析可知，分接开关经itd分解后，其机械运行状态信息主要集中在前5个分量，因此本发明使用itd分解得到的前5个分量的能量熵构造特征向量，作为svm的输入。对svm进行训练，将构造的特征向量输入到svm中，对svm进行训练，将测试数据输入到训练好的svm从而判断oltc的故障模式。采用同样的15组数据对bp网络进行训练，并对这15组数据进行预测，表3为svm的预测结果与bp网络预测结果的比较：表3svm的预测结果与bp网络预测结果的比较注：1表示分接开关正常状态，2表示分接开关处于触头松动故障，3表示分接开关弹簧性能下降故障，每种状态5组数据，共15组，表1给出的是统计后的结果比较。由表3中可以看出，基于itd和svm的oltc机械故障诊断方法可以有效地识别出故障类型，从而说明了本发明的可行性，由表3也知，本发明诊断效果优于bp网络。以上所述仅是本发明的优选实施方式，应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。当前第1页12