一种基于希尔伯特黄变换算法对断路器故障诊断的方法技术

本发明专利技术公开了一种基于希尔伯特黄变换算法对断路器故障诊断的方法，对数据信号进行经验模态分解，信号经过经验模态分解得到各固有模态函数分量，对每一个固有模态函数分量进行Hilbert变换，得到瞬时幅值和瞬时频率参数，将振幅和瞬时频率绘制成能量、频率、时间分布的Hilbert谱，应用于非线性、非平稳系列的分析中，并采用镜像延拓方法抑制Hilbert变换端点效应。测试结果证实了本发明专利技术在电力系统断路器振动监测分析中具有灵活性和有效性，使用EMD可将信号分解为n个IMF分量，振动信号所含的不同频率成分应包含在所分解的IMF分量中，其中IMF2中的频率分布恰好包含了分离开关的基础频率部分，提高了分析的精度，简化了分析的步骤。

A method for fault diagnosis of circuit breaker based on Hilbert Huang transform method

The invention discloses a method based on Hilbert Huang transform method for fault diagnosis of circuit breaker, the data signal, empirical mode decomposition, EMD signal after each IMFs, Hilbert transformation to every IMFs, get the instantaneous amplitude and instantaneous frequency, amplitude and the instantaneous frequency of draw energy, frequency and time distribution of the Hilbert spectrum is applied to analyze the nonlinear and non-stationary series, and the mirror extension method to restrain end effects of Hilbert transform. The test results confirmed that the invention has the advantages of flexibility and effectiveness in the analysis of power system circuit breaker vibration monitoring, the use of EMD signals can be decomposed into n IMF components, different frequency components contained in vibration signal should be included in the decomposition of the IMF component, the IMF2 frequency distribution just includes the basic frequency separation the switch part, improve the analysis accuracy, simplifies the analysis steps.

【技术实现步骤摘要】

本专利技术属于电力系统监测领域，尤其涉及一种基于希尔伯特黄变换算法对断路器故障诊断的方法。
技术介绍
电力系统中，高压断路器数量众多、检修量大、费用高。有关统计表明，变电站维护费用的一半以上是用在高压断路器上，而其中60％又是用于断路器的小修和例行检修上。另外据统计，10％的断路器故障是由于不正确的检修所致，断路器的大修完全解体，既费时间，费用也很高，而且解体和重新装配有可能会引起很多新的缺陷。在目前相对保守的计划检修中，检修缺乏一定的针对性。对断路器的重要参数进行长期连续的在线监测，不仅可以准确、实时地反映设备现在的运行状态，而且还能分析各种重要参数的变化趋势，预测使用寿命，判断是否存在故障先兆，为检修决策提供依据，实现设备预知性维修，从而延长设备的维修保养周期，提高设备的利用率，减少维修保养费用，形成一套融故障检修、定期检修、状态检修和主动检修为一体的、优化的综合检修方式，以保证电力系统安全可靠运行，并可以大大减少盲目定期检修带来的资金浪费，降低电力成本。随着电网建设规模的扩大，电网对于一次设备的要求也越来越高。断路器是电力系统中最重要的控制与保护设备，它在电网中的作用至关重要，对其实施在线监测，提前发现故障征兆，必将是今后发展的趋势。基于FFT技术的经典频谱估计由于是对信号真实功率谱的渐近无偏估计，而非一致估计。当分析样本长度增加时，频谱估计的方差、偏差和分辨率并不能同时得到改善；另外，由于实时工程检测信号的非平稳性，必然会存在非同步采样误差，加窗截断处理往往会导致频谱能量泄漏和频谱估计偏差，最终影响机电设备故障诊断的准确性。本文采用基于希尔伯特黄变换(HHT)的现代谱估计方法对断路器的非平稳信号进行频谱估计，以期进一步提高其分析精度和准确性。
技术实现思路
本专利技术的目的在于提供一种基于希尔伯特黄变换算法对断路器故障诊断的方法，旨在解决现有的断路器的计划检修缺乏针对性，对断路器的非平稳信号的分析精度和准确性不够高的问题。本专利技术是这样实现的，一种基于希尔伯特黄变换算法对断路器故障诊断的方法对数据信号进行经验模态分解，信号经过经验模态分解得到各固有模态函数分量，对每一个固有模态函数分量进行Hilbert变换，得到瞬时幅值和瞬时频率参数，将振幅和瞬时频率绘制成能量、频率、时间分布的Hilbert谱，应用于非线性、非平稳系列的分析中，并采用镜像延拓方法抑制Hilbert变换端点效应。进一步，对于一个任意的信号X(t)，采用下式进行Hilbert变换：Y(t)=1πP∫-∞+∞X(τ)t-τdτ]]>其中，P是Cauchy主值，当时间序列X(t)的频率f＞0时，P＝-1；根据式上式得到解析信号：Z(t)＝X(t)+iY(t)＝a(t)eiθ(t)其中，a(t)=[X2(t)+Y2(t)]1/2θ(t)=arctan(Y(t)X(t))]]>a(t)为瞬时振幅，θ(t)为瞬时角速度，ω(t)瞬时角频率为：ω(t)=dθ(t)dt]]>则瞬时频率：f(t)=ω(t)2π=12π·ω(t)]]>瞬时频率只有唯一的与其对应的值。进一步，经验模态分解的过程为：求取原信号在时间长度内的所有的极值点，将极大值点经过三次样条插值函数拟合形成原信号的上包络线；将所有的极小值点通过三次样条插值函数拟合形成原信号的下包络线，上、下包络线的均值记作m1(t)；原信号的IMF分量h1(t)由下式得：h1(t)＝X(t)-m1(t)将第一次筛选出来的信号h1(t)作为被筛选的原始信号，按下式进行重复筛选：h11(t)＝h1(t)-m11(t)重复筛选过程k次，直到得到的第一个h1(t)符合固有模态函数的要求，得到原始信号经过经验模态分解的第一个固有模态函数分量c1(t)，即：h1k(t)＝h1(k-1)(t)-m1k(t)c1(t)＝h1k(t)当两次筛选的结果满足下式时，并设定SD值下式作为停止筛选的判据，SD表示成下式：SD=Σt=0τ[(h1k(t)-h1(k-1)(t))h1(k-1)(t)]2,k=1,2...]]>SD的限值可以根据经验来设定，前后两次筛选结果的SD值达到限值，即可以停止此次筛选过程；当一次筛选过程结束时，将c1(t)从原信号X(t)分离出来，得到一个去掉高频分量c1(t)的剩余信号r1(t)：r1(t)＝X(t)-c1(t)当一次筛选过程结束时，要将r1(t)作为新信号施行重复的筛选过程，得到n个固有模态函数分量：r2(t)＝r1(t)-c2(t)，…，rn(t)＝rn-1(t)-cn(t)得到一系列的固有模态函数分量：X(t)=Σi=1n+1ci(t)=Σi=1nci(t)+rn(t)]]>cn+1(t)＝rn(t)经过分解原始信号X(t)经逐级的分解后，分解成得到n个固有模态函数分量c1(t)，c2(t)，…cn(t)和一个趋势分量rn(t)。进一步，经过希尔伯特黄变换，将原始信号x(t)用各个固有模态函数的瞬时频率ωj(t)和瞬时幅值aj(t)近似的表示成下式：x(t)=Σj=1naj(t)ej∫ωj(t)dt]]>其中aj(t)和ωj(t)分别表示第j个固有模态函数分量的瞬时幅值和频率。进一步，对原始信号采用下式进行Fourier展开：x(t)=Σj=1∞aj(t)ejω]]>其中，aj和ωj分别表示第j个谐波分量的幅值和频率，它们不随着时间而变，是一个固定值。进一步，对Hilbert谱在时间上进行积分得到边际谱：h(ω)=∫0τH(ω,t)dt]]>其中，H(ω，t)表示Hilbert谱。测试结果证实了本专利技术在电力系统断路器振动监测分析中具有灵活性和有效性。使用EMD可将信号分解为n个IMF分量，振动信号所含的不同频率成分应包含在所分解的IMF分量中，其中IMF2中的频率分布恰好包含了分离开关的基础频率部分，提高了分析的精度，简化了分析的步骤。附图说明图1是本专利技术实施例提供的标准端点延拓措施IMF分量时域图；图2是本专利技术实施例提供的端点镜像延拓IMF分量时域图；图3是本专利技术实施例提供的标准端点延拓措施IMF分量频谱图；图4是本专利技术实施例提供的端点镜像延拓IMF分量频谱图；图5是本专利技术实施例提供的断路器处于未合闸状态，基频为2.7hz时使用EMD分解的IMF2分量进行分析的结果图；图6是本专利技术实施例提供的断路器处于合闸中间状态，基频为3.3hz时使用EMD分解的IMF2分量进行分析的结果图；图7是本专利技术实施例提供的断路器处于合紧状态，基频为3.4hz时使用EMD分解的IMF2分量进行分析的结果图。具体实施方式为能进一步了解本专利技术的
技术实现思路
、特点及功效，兹例举以下实施例，并配合附图详细说明如下。请参阅图1-图7：本专利技术是这样实现的，一种基于希尔伯特黄变换算法对断路器故障诊断的方法对数据信号进行经验模态分解，信号经过经验模态分解得到各固有模态函本文档来自技高网...

【技术保护点】
一种基于希尔伯特黄变换算法对断路器故障诊断的方法，其特征在于，所述的基于希尔伯特黄变换算法对断路器故障诊断的方法对数据信号进行经验模态分解，信号经过经验模态分解得到各固有模态函数分量，对每一个固有模态函数分量进行Hilbert变换，得到瞬时幅值和瞬时频率参数，将振幅和瞬时频率绘制成能量、频率、时间分布的Hilbert谱，应用于非线性、非平稳系列的分析中，并采用镜像延拓方法抑制Hilbert变换端点效应。

【技术特征摘要】
1.一种基于希尔伯特黄变换算法对断路器故障诊断的方法，其特征在于，所述的基于希尔伯特黄变换算法对断路器故障诊断的方法对数据信号进行经验模态分解，信号经过经验模态分解得到各固有模态函数分量，对每一个固有模态函数分量进行Hilbert变换，得到瞬时幅值和瞬时频率参数，将振幅和瞬时频率绘制成能量、频率、时间分布的Hilbert谱，应用于非线性、非平稳系列的分析中，并采用镜像延拓方法抑制Hilbert变换端点效应。2.如权利要求1所述的基于希尔伯特黄变换算法对断路器故障诊断的方法，其特征在于，对于一个任意的信号X(t)，采用下式进行Hilbert变换：Y(t)=1πP∫-∞+∞X(τ)t-τdτ]]>其中，P是Cauchy主值，当时间序列X(t)的频率f＞0时，P＝-1；根据式上式得到解析信号：Z(t)＝X(t)+iY(t)＝a(t)eiθ(t)其中，a(t)=[X2(t)+Y2(t)]1/2θ(t)=arctan(Y(t)X(t))]]>a(t)为瞬时振幅，θ(t)为瞬时角速度，ω(t)瞬时角频率为：ω(t)=dθ(t)dt]]>则瞬时频率：f(t)=ω(t)2π=12π·ω(t)]]>瞬时频率只有唯一的与其对应的值。3.如权利要求1所述的基于希尔伯特黄变换算法对断路器故障诊断的方法，其特征在于，经验模态分解的过程为：求取原信号在时间长度内的所有的极值点，将极大值点经过三次样条插值函数拟合形成原信号的上包络线；将所有的极小值点通过三次样条插值函数拟合形成原信号的下包络线，上、下包络线的均值记作m1(t)；原信号的IMF分量h1(t)由下式得：h1(t)＝X(t)-m1(t)将第一次筛选出来的信号h1(t)作为被筛选的原始信号，按下式进行重复筛选：h11(t)＝h1(t)-m11(t)重复筛选过程k次，直到得到的第一个h1(t)符合固有模态函数的要求，得到原始信号经过经验模态分解的第一个固有模态函数分量c1(t)，即：h1k(t)＝h1(k-1...

【专利技术属性】
技术研发人员：唐远，李健，徐珂航，朱敏，徐智，
申请(专利权)人：国网四川省电力公司眉山供电公司，南京优能特电力科技发展有限公司，
类型：发明
国别省市：四川;51