基于小波多尺度变换快速辨识CVT故障的方法与流程

本发明涉及故障识别与电力系统测量装置领域，具体涉及一种基于小波多尺度变换快速辨识cvt故障的方法。

背景技术：

近年来，随着我国现代电网的不断发展，已呈现出智能化、特高压化的特点和趋势，对于电力互感器提出了更高的要求来适应这种系统电压等级和输送容量不断提高、网络智能化程度不断加深的发展形势。

电容式电压互感器(capacitorvoltagetransformer，cvt)因具有结构简单、经济安全、维护工作容易、绝缘可靠性高等特性广泛使用于高电压监测工作当中。但受限于cvt严格密封的制作工艺等因素，在实际运行过程中，除了发生漏油、产生异常声响等易直观发现的故障外，内部故障比如电容器击穿、二次端子松动等运行缺陷时有发生且不易发现，导致在线路投切、补偿装置等投入或退出过程中cvt故障率居高不下，并造成监督排查工作的复杂性，威胁电网的安全运行。

从大量电气工程实践中，人们认识到要使cvt设备可靠有效的运行，充分发挥其效益，必须发展设备监测和故障诊断技术，因而cvt故障识别具有重要的研究价值。

技术实现要素：

针对现有技术中的上述不足，本发明提供的基于小波多尺度变换快速辨识cvt故障的方法通过对cvt故障检测可以避免cvt采集不准后加速失效、完全击穿等缺陷。

为了达到上述发明目的，本发明采用的技术方案为：

提供一种基于小波多尺度变换快速辨识cvt故障的方法，其包括：

a1、采集在线电力系统中各cvt装置的电压信号；

a2、获取已选取的小波基和阈值选择方法，并根据选取的小波基，采用小波多尺度变换将电压信号分解得到尺度系数和小波系数；

a3、采用阈值选择方法对小波系数进行阈值调整，并将阈值调整后的各级小波系数按小波变换的反演算法进行信号重构，得到消噪信号；

a4、根据选取的小波基，采用小波多尺度变换对消噪信号进行分解；

a5、选取2、3、4三个尺度，采用模极大值法对单尺度重构信号进行处理，得到三个电压信号；

a6、根据重构后的三个电压信号的波形，判断cvt装置是否发生故障，若发生，则输出故障的位置点，否则返回步骤a1。

本发明的有益效果为：本方案相比传统算法对一、二次测量电压对比分析及引入其他判定量容易限于数据处理量繁琐、时效性不够好等问题，本方案利用小波多尺度分析可实现实时可视化监测、快速辨别故障波形、对比判定故障点等，为调度监测人员进行下一步cvt故障诊断提供了有效、可行的分析途径。

利用模极大值法捕捉到异常信号的模极大值点，能很快辨识到cvt异常信号发生的位置；同时结合多尺度重构算法捕捉到cvt异常波形。消噪信号具体实现过程中削弱了信号的高频噪声，有利于提高故障识别和定位能力。

进一步地，所述采用小波多尺度变换将电压信号分解得到尺度系数和小波系数为采用一维离散小波将电压信号f(x)分解为尺度系数ck和小波系数dk，f(x)的表达式为：

其中，k＝0，1，2，…，n-1；n为总样本数；j为分解的层数；构成l²空间的一个子空间vj+1的正交基；和的计算式分别为：

其中，l(n)为低通滤波器；h(n)为高通滤波器；z为正整数。

采用上述进一步方案的有益效果是：本方案通过已选取的小波基对电压信号进行分解得到的上述表达式，可以提高后续故障识别和定位能力。

进一步地，所述小波基和阈值选择方法的选取方法进一步包括：

s1、获取cvt装置的历史电压信号，并设置不同时序位置处的电压波动分别为1％和2％的击穿故障；

s2、随机加入固定长度的噪声，采用多种候选小波基将加入击穿故障的历史电压信号分解得到尺度系数和小波系数；

s3、采用多种阈值选择方法分别对同一个小波系数进行阈值调整，并将阈值调整后的各级小波系数按小波变换的反演算法进行信号重构，得到消噪信号；

s4、计算每个小波系数得到的多个消噪信号的信噪比，并选取同一个小波系数对应的多个消噪信号中信噪比最大的消噪信号，同时将信噪比最大的消噪信号对应的阈值选择方法的选择次数累加一次；

s5、采用选取的消噪信号对应的候选小波基将消噪信号分解得到尺度系数和小波系数；

s6、采用模极大值法对步骤s5得到的尺度系数和小波系数进行处理，得到每个尺度信号模极大值点的坐标，并与击穿故障处的位置信息进行比较；

s7、当比较结果小于预设阈值时，则将对应的候选小波基的有效次数累加一次，并将重复执行步骤s1至步骤s7的循环次数累加一次；

s8、当循环次数等于设定次数时，进入步骤s9，否则返回步骤s1；

s9、对比所有候选小波基的有效次数，选取有效次数最大对应的候选小波基作为在线故障识别时的小波基；

对比所有阈值选择方法对应的选择次数，选取选择次数最大对应的阈值选择方法作为在线故障识别时的阈值选择方法。

采用上述进一步方案的有益效果是：采用统计实验确定最优小波基和阈值选择方法，利用模极大值法捕捉到异常信号的模极大值点，能很快辨识到cvt异常信号发生的位置；由于在在线检测之前对最优小波基和阈值选择方法进行确认，可以简化在线检测步骤，提高检测效率。

进一步地，所述信噪比的计算公式为：

其中，qsnr为信噪比；ps为历史电压信号的功率；pn为噪声的功率；si为历史电压信号；fi为消噪信号。

采用上述进一步方案的有益效果是：基于噪声评价确定各阈值选择方式，将不确定的去噪效果反应到具体物理指标上，利于确定各阈值选择方式的特点。

进一步地，每次循环时采用的历史电压信号均不相同。由于选取的历史电压信号均不相同，这样可以提高统计最优小波基和阈值选择方法的准确性。

进一步地，所述候选小波基包括选择haar小波、daubechies系列小波、coiflet系列小波、symlets系列小波、reversebior系列小波和fejerkorovkin系列小波；多种阈值选择方法分别为无偏风险估计阀值、启发式阈值、固定门限阈值和最大最小准则阈值。

附图说明

图1为基于小波多尺度变换快速辨识cvt故障的方法一个实施例的原理框图。

图2为实验例中cvt电气结构原理图。

图3为实验例中小波系数第2、3、4三个尺度来重构采集的电压信号的波形图。

具体实施方式

下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

电容式电压互感器是由电容分压器和电磁单位(中间变压器)组合而成的具有独特结构的电气设备。图2为cvt电气结构简图，电容分压器分高压电容c1和中压电容c2两部分，其中高压电容c1由c11、c12和c133部分串联而成。c11、c12、c13分布在cvt的上节瓷套管中，中压电容c2分布在下节瓷套管中。c11、c12和c13都可以看成一定数量的电容原件串联而成；图1中l为补偿电抗器，f为保护装置，a1为主二次绕组，da和ds构成剩余电压绕组。电容原件经真空干燥处理后浸在电容器油中，串连在一起，由绝缘支撑杆固定。

cvt内部故障通常表现在电容元件被击穿，这会引起总电容的变化，出现分压比变化和一、二次电压异常等异常现象。设cvt内部各电容原件电容量相等，则有电容变化率：

其中，n0个电容元件串联，n1个电容被击穿。电容δc％变化一一对应cvt测量电压变化δu％，根据各厂站的测试和理论计算给δu％设定一个阈值数值(常用1％、2％等)，若超过该数值，则表示该处cvt异常。

但由于实际操作过程cvt击穿元件较少，即n1远小于n0，因此电容变化率不大，电压变化率不大，无法直接通过监测二次电压发现监测电压的变化情况。通过实用有效的方法敏锐地捕捉到电压值的变化是表征cvt采集波形异常的关键。

参考图1，图1示出了基于小波多尺度变换快速辨识cvt故障的方法一个实施例的原理框图；如图1所示，该方法包括步骤a1至步骤a6。

在步骤a1中，采集在线电力系统中各cvt装置的电压信号。

在步骤a2中，获取已选取的小波基和阈值选择方法，并根据选取的小波基，采用小波多尺度变换将电压信号分解得到尺度系数和小波系数。

实施时，本方案优选所述采用小波多尺度变换将电压信号分解得到尺度系数和小波系数为采用一维离散小波将电压信号f(x)分解为尺度系数ck和小波系数dk，f(x)的表达式为：

其中，k＝0，1，2，…，n-1；n为总样本数；j为分解的层数；构成l²空间的一个子空间vj+1的正交基；和的计算式分别为：

其中，l(n)为低通滤波器；h(n)为高通滤波器；z为正整数。

在步骤a3中，采用阈值选择方法对小波系数进行阈值调整，并将阈值调整后的各级小波系数按小波变换的反演算法进行信号重构，得到消噪信号。

其中，阈值选择方法的阈值函数主要分为硬阈值函数和软阈值函数；硬阈值处理对小波分解得到的较小函数进行置零，其函数形式为：

软阈值的计算方法是将较小的系数置零但对较大的小波系数趋向零收缩，其函数形式为：

目前，工程中常用阈值选择方法有无偏风险估计阀值(rigrsur)、启发式阈值(heursure)、固定门限阈值(sqtwolog)、最大最小准则阈值(minimaxi)共四种。

在步骤a4中，根据选取的小波基，采用小波多尺度变换对消噪信号进行分解；此处分解所采用的算例及达到的效果与步骤a2完全相同，此处就不再赘述。

在步骤a5中，选取2、3、4三个尺度，采用模极大值法对单尺度重构信号进行处理，得到三个电压信号。

模极大值法是基于函数奇异性理论来检测信号在某一序列中是否具有突变性的，并以来lipschitzα来描述，小波变换的模极大值被用来检测信号的奇异性，设wψf(s,x)(s＝2^j)是信号f(x)在尺度j下的小波变换，对于任意x，x∈(x0-δ,x0+δ)，有：

|wψf(s,x)|≤|wψf(s,x0)|

则x0称为小波变换在尺度j下的模极大值点，|wψf(s,x0)|为小波变换的模极大值。模极大值点与信号的奇异点一一对应，模极大值点对应信号的奇异点跳变的边缘，模极大值的幅值对应信号奇变的强度。

在步骤a6中，根据重构后的三个电压信号的波形，判断cvt装置是否发生故障，若发生，则输出故障的位置点，否则返回步骤a1。由于发生故障处的波形会发生突变，呈现出突然上升并快速下降，通过波形能够直观地发现这一现象，而快速确定有故障的点。

在本发明的一个实施例中，小波基和阈值选择方法的选取方法进一步包括：

s1、获取cvt装置的历史电压信号，并设置不同时序位置处的电压波动分别为1％和2％的击穿故障；

s2、随机加入固定长度的噪声，采用多种候选小波基将加入击穿故障的历史电压信号分解得到尺度系数和小波系数；

s3、采用多种阈值选择方法分别对同一个小波系数进行阈值调整，并将阈值调整后的各级小波系数按小波变换的反演算法进行信号重构，得到消噪信号；

s4、计算每个小波系数得到的多个消噪信号的信噪比：

其中，qsnr为信噪比；ps为历史电压信号的功率；pn为噪声的功率；si为历史电压信号；fi为消噪信号。

并选取同一个小波系数对应的多个消噪信号中信噪比最大的消噪信号，同时将信噪比最大的消噪信号对应的阈值选择方法的选择次数累加一次。

s5、采用选取的消噪信号对应的候选小波基将消噪信号分解得到尺度系数和小波系数；

s6、采用模极大值法对步骤s5得到的尺度系数和小波系数进行处理，得到每个尺度信号模极大值点的坐标，并与击穿故障处的位置信息进行比较；

s7、当比较结果小于预设阈值时，则将对应的候选小波基的有效次数累加一次，并将重复执行步骤s1至步骤s7的循环次数累加一次；

s8、当循环次数等于设定次数时，进入步骤s9，否则返回步骤s1，在进入下一次循环时，每次循环时采用的历史电压信号均不相同。

s9、对比所有候选小波基的有效次数，选取有效次数最大对应的候选小波基作为在线故障识别时的小波基；

对比所有阈值选择方法对应的选择次数，选取选择次数最大对应的阈值选择方法作为在线故障识别时的阈值选择方法。

本方案运用小波多分辨分析方法对阈值去噪进行了一定的探讨，以信噪比的大小作为评价，并结合cvt运行和采样的特点找到了适用于cvt电压去噪的阈值选择规则。

本方案首先通过步骤s1至s9所建立的评价指标，结合统计实验的思想确定适合的阈值选择方法或小波基，确定最优小波参数后，再对实时数据取样后进行去噪和波形检测，这样躲避了由于小波簇理论并不完善而常在应用问题上缺乏系统的最优参数选择问题，提高了cvt监测的可靠性、运行的安全性、计量电压的公平性等优点，同时本发明具备的在线识别和可视化特点可提高变电站等高电压场所无人看守、自主运行等能力。

实施时，本方案所采用的候选小波基包括选择haar小波、daubechies系列小波、coiflet系列小波、symlets系列小波、reversebior系列小波和fejerkorovkin系列小波；多种阈值选择方法分别为无偏风险估计阀值、启发式阈值、固定门限阈值和最大最小准则阈值。

下面结合具体的实例对本方案快速辨识cvt故障的方法的效果进行说明：

选择某500kv变电站数据库采样100个点，随机选择5个点20、40、60、80、99设置δu＝1％的突变信号，去噪后使用rbio3.1小波基进行模极大值辨别，计算出三种情况下的模极大值，其结果如表1所示，从表1中可以看出三种分解尺度下求取的模极大值点与预设故障点位置十分接近，说明本方法的准确性非常高，同时注意到不同分解尺度对应求出的模极大值点存在一定的差异，但至少有两个点非常准确，保证了本方法的可靠性。

表1不同分解尺度下的模极大值点

选择上述一个预设故障点80进行第2、3、4三个层次进行重构，通过重构后的波形参见图3，图3提取出了故障信号三种分解尺度下重构的波形，容易看到三种尺度下采样位置均在80左右时检测到了信号的模极大值，而且突变能力很强，表现在于采样点很小的区域内就出现显著的奇变。

输入信号由于采样间隔较长，波形变化幅度较大，基于数据比较的传统监测方法难保证分辨准确度。采用模极大值法，首先将不同尺度下的故障波形提取出来，并将非明显变化区间“归零”处理，可以快速“聚焦”到异常信号段，这使得监测工作更加直观和容易。

而三种不同的单支重构信号代表着三种不同的分辨率，三个波形在模极大值点附近都产生异常波形，且第二层和第四层波形在模极大值点附近都能明显观测到信号突变的特点，且两者波形相似，可以确保判定该模极大值点出现时cvt发生异常。