CN107016352B - 电力变压器有载分接开关振动脉冲数自动识别方法 - Google Patents
电力变压器有载分接开关振动脉冲数自动识别方法 Download PDFInfo
- Publication number
- CN107016352B CN107016352B CN201710146756.0A CN201710146756A CN107016352B CN 107016352 B CN107016352 B CN 107016352B CN 201710146756 A CN201710146756 A CN 201710146756A CN 107016352 B CN107016352 B CN 107016352B
- Authority
- CN
- China
- Prior art keywords
- signal
- vibratory impulse
- value
- model
- mixed signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
- G06F2218/16—Classification; Matching by matching signal segments
- G06F2218/20—Classification; Matching by matching signal segments by applying autoregressive analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- General Engineering & Computer Science (AREA)
- Housings And Mounting Of Transformers (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种电力变压器有载分接开关振动脉冲数自动识别方法,步骤为:在变压器油箱表面获得振动信号,并进行零均值处理;以一个采样间隔为分割点将所述混合信号分为前后两个统计时间段内的自回归模型,并逐次计算对应的回归系数和模型阶数;利用最大似然估计公式对自回归模型的误差进行最大似然估计;计算信号回归模型误差的信息熵A值;求取信息熵A值的最小值和最大值;根据有载分接开关的类型,选择有效的振动脉冲信号长度区间,求取信号长度区间对应的信息熵A值的梯度;应用梯度法求取局部最大值的经典算法,获得每个包络的最大值,最大值得个数即对应于振动脉冲的数目。本发明能够实现对有载分接开关多个振动脉冲信号的自动识别。
Description
技术领域
本发明属于电气设备监测技术领域,,特别涉及到电力变压器有载分接开关振动脉冲信号与变压器油箱表面振动信号的自适应提取方法。
背景技术
有载调压分接开关(On-Load Tap Changer,OLTC)是变压器完成调压的核心部件,不仅能稳定负载中心电压,而且也是联络电网、调整负载潮流、改善无功分配等不可缺少的重要设备,它的性能状况直接关系到有载调压变压器的安全运行。发达国家对容量在10MVA及以上的变压器大都安装了有载分接开关,我国也愈来愈多地采用了有载调压变压器。有载分接开关由选择器、切换开关和电动机构组成,机械故障是其主要故障类型,包括选择开关触头接触不良、操作机构失灵、限位开关失灵、切换开关拒切、中止或动作滞后、内部紧固件松动和脱落等。这些故障可以通过对有载分接开关振动信号的分析进行有效的检测和判断。由于电力变压器在运行时,铁心硅钢片的磁致伸缩引起铁心振动,绕组在负载电流电场力作用下产生绕组振动,这些振动都将传递到变压器的油箱表面,形成变压器的主要振动。在变压器油箱表面获得的振动信号除了有载分接开关振动以外,还包括铁心和绕组的振动。因此要对有载分接开关振动信号进行分析和处理,首先需要对油箱表面振动信号进行信号提取,从中提取出有载分接开关振动信脉冲号成分。
有载分接开关动作振动波形主要是由动静触头的接触摩擦引起的振动冲击脉冲信号,而铁心和绕组的振动信号主要为周期信号,同时有载分接开关在一次动作过程中,多个触点将发生接触与提取,产生的振动波形将存在多个振动冲击脉冲,振动冲击脉冲信号与分接开关触点的接触状态相关,能够反映出触点接触的状态。
发明内容
本发明的目的是提供一种电力变压器有载分接开关振动脉冲数自动识别方法,以实现对有载分接开关多个振动脉冲信号的自动识别。
为实现上述目的,本发明采用的技术方案为:
一种电力变压器有载分接开关振动脉冲数自动识别方法,包括以下步骤:
(1)在变压器油箱表面获得振动信号,该振动信号是包括有载分接开关振动脉冲信号、环境噪声、变压器铁心振动信号的混合信号;并对所述混合信号进行零均值处理;
(2)根据所述混合信号的统计特性,以一个采样间隔为分割点将所述混合信号分为前后两个统计时间段内的自回归模型,并逐次计算对应的回归系数和模型阶数;
(3)利用最大似然估计公式对自回归模型的误差进行最大似然估计;
(4)计算信号回归模型误差的信息熵A值;
(5)求取信息熵A值的最小值和最大值;该最小值对应噪声自回归模型和振动脉冲信号自回归模型的最佳分割点,是在最小二乘意义下两个模型的最佳匹配,是分接开关触头机械动作产生的振动脉冲信号的起始点;
(6)根据有载分接开关的类型,选择有效的振动脉冲信号长度区间,求取信号长度区间对应的信息熵A值的梯度;则可对有效信号长度内的多个振动脉冲信号形成统计意义上的包络;从第二个包络开始,包络的最大值表示振动脉冲信号的信息熵变化率最大,即前后两个振动脉冲信号在统计意义上的分割点,标志着后一个振动脉冲信号起始点;
(7)应用梯度法求取局部最大值的经典算法,获得每个包络的最大值,最大值得个数即对应于振动脉冲的数目。
有益效果:本发明的方法利用振动冲击脉冲信号和环境噪声信号的统计独立特性,在对数据进行处理及分段的基础上,通过计算信号的信息熵来判断脉冲信号的起始时间,进一步计算信息熵的梯度,并在运用求取最大值和最小值经典算法的基础上,结合设定的阈值来实现对有载分接开关多个振动脉冲信号的自动识别。
附图说明
图1为分接开关振动信号;
图2为分接开关振动信号信息熵曲线;
图3为分接开关振动信号脉冲数识别图形。
具体实施方式
下面结合附图对本发明做更进一步的解释。
一种电力变压器有载分接开关振动脉冲数自动识别方法,包括以下步骤:
(1)在变压器油箱表面获得振动信号,该振动信号是包括有载分接开关振动脉冲信号、环境噪声、变压器铁心振动信号的混合信号;并对所述混合信号进行零均值处理;
变压器油箱表面采集得到的混合信号为xn={x1,...,xN},则混合信号去除均值的方法为:
式中,表示去掉均值后的第n个信号的幅值,xn表示采集得到的混合信号时间序列中第n个信号的幅值,N表示采集得到的混合信号时间序列的总个数,xm表示在第m个混合信号的幅值,下标n、m分别表示采集到的混合信号中的第n、m个信号。
(2)根据所述混合信号的统计特性,以一个采样间隔为分割点将所述混合信号分为前后两个统计时间段内的自回归模型,并逐次计算对应的回归系数和模型阶数;
具体方法为:
变压器表面采集得到的混合信号的时间序列为xn={x1,...,xN}。
首先进行第一次数据分段,从第一个数据x1开始,把混合信号的时间序列分成{x1}和{x2,...,xN}两个时间段,如图2所示。
然后进行第二次分段,第一个时间段增加一个数据x2,第二个时间段减少一个数据x2,即两个时间段分别为{x1,x2}和{x3,...,xN}。依次递推,对于总共有N个数据的混合信号,则总共需要N-1次数据分段。
对于第j次分段后的两段数据{x1,x2,...,xj}和{xj+1,xj+2,...,xN},分别用下式建立第一段数据和第二段数据的自回归模型:
式中,为自回归模型的系数,Mi(i=1,2)分别为噪声模型的阶数和信号模型的阶数。该模型将时间序列在模型窗内分为确定的和预测误差若误差为高斯分布,则其均值方差且与时间序列确定的部分不相关
(3)利用最大似然估计公式对自回归模型的误差进行最大似然估计;具体方法为:
把j作为分割点,使用自回归模型的系数来提取时间序列[M1+1,j]和[j+1,N-M2]中的预测误差部分。假设误差为高斯分布,则这两个时段误差的似然函数L可用下式表示:
式中分别表示时段1和时段2的模型参数(依赖于j),而且p1=M1+1,p2=j+1,q1=j,q2=N-M2,n1=j-M1,n2=N-M2-j;
若把K作为分割点,使用自回归模型的系数来提取时间序列[M1+1,K]和[K+1,N-M2]中的预测误差部分。
对方程最大似然函数取对数并求偏导数,寻找模型参数的最大似然估计,如下式所示:
求解上式,可求得模型参数:
通过实验研究发现,对于变压器表面测得的混合信号自回归模型误差参数的最大似然估计,可以用下式代替:
上式中,下标k表示采集得到的混合信号中的第k个信号。
(4)对于第j次分段,通过下式计算信号回归模型误差的信息熵A值,计算结果如附图2所示:
上式中,C=(M1+M2-N)(log2π+1),为常数。
(5)采用MINMAX算法获取信息熵A值的最小值和最大值。
(6)根据有载分接开关的类型,选择有效的振动脉冲信号长度区间,采用DIFF算法求取信号长度区间对应的信息熵A值的梯度,则可对有效信号长度内的多个振动脉冲信号形成统计意义上的包络。
(7)应用梯度法求取局部最大值的经典算法,可以获得每个包络的最大值,最大值得个数对应于脉冲的数目,如图3所示。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (7)
1.一种电力变压器有载分接开关振动脉冲数自动识别方法,其特征在于:包括以下步骤:
(1)在变压器油箱表面获得振动信号,该振动信号是包括有载分接开关振动脉冲信号、环境噪声、变压器铁心振动信号的混合信号;并对所述混合信号进行零均值处理;
(2)根据所述混合信号的统计特性,以一个采样间隔为分割点将所述混合信号分为前后两个统计时间段内的自回归模型,并逐次计算对应的回归系数和模型阶数;
(3)利用最大似然估计公式对自回归模型的误差进行最大似然估计;
(4)计算信号回归模型误差的信息熵A值;
(5)求取信息熵A值的最小值和最大值;
(6)根据有载分接开关的类型,选择有效的振动脉冲信号长度区间,求取信号长度区间对应的信息熵A值的梯度;则可对有效信号长度内的多个振动脉冲信号形成统计意义上的包络;
(7)应用梯度法求取局部最大值的经典算法,获得每个包络的最大值,最大值得个数即对应于振动脉冲的数目。
2.根据权利要求1所述的电力变压器有载分接开关振动脉冲数自动识别方法,其特征在于:所述步骤(1)中,变压器油箱表面采集得到的混合信号为xn={x1,...,xN},则通过下式对混合信号进行零均值处理:
式中,表示去掉均值后的第n个信号的幅值,xn表示采集得到的混合信号时间序列中第n个信号的幅值,N表示采集得到的混合信号时间序列的总个数,xm表示在第m个混合信号的幅值,下标n、m分别表示采集到的混合信号中的第n、m个信号。
3.根据权利要求1所述的电力变压器有载分接开关振动脉冲数自动识别方法,其特征在于:所述步骤(2)的具体方法为:
变压器油箱表面采集得到的混合信号的时间序列为xn={x1,...,xN};
首先进行第一次数据分段,从第一个数据x1开始,把混合信号的时间序列分成{x1}和{x2,...,xN}两个时间段;
然后进行第二次分段,第一个时间段增加一个数据x2,第二个时间段减少一个数据x2,即两个时间段分别为{x1,x2}和{x3,...,xN};依次递推,对于总共有N个数据的混合信号,则总共需要N-1次数据分段;
对于第j次分段后的两段数据{x1,x2,...,xj}和{xj+1,xj+2,...,xN},分别用下式建立第一段数据和第二段数据的自回归模型:
其中t=1,2;
式中,为自回归模型的系数,其中m=1,...,Mi,Mi分别为噪声模型的阶数和信号模型的阶数,其中i=1,2;该模型将时间序列在模型窗内分为确定的和预测误差若误差为高斯分布,则其均值方差且与时间序列确定的部分不相关
4.根据权利要求1所述的电力变压器有载分接开关振动脉冲数自动识别方法,其特征在于:所述步骤(3)的具体方法为:
把j作为分割点,使用自回归模型的系数来提取时间序列[M1+1,j]和[j+1,N-M2]中的预测误差部分;假设误差为高斯分布,则这两个时段误差的似然函数L可用下式表示:
式中分别表示时段1和时段2的模型参数,其中依赖于j,而且p1=M1+1,p2=j+1,q1=j,q2=N-M2,n1=j-M1,n2=N-M2-j;
若把K作为分割点,使用自回归模型的系数来提取时间序列[M1+1,K]和[K+1,N-M2]中的预测误差部分;
对方程最大似然函数取对数并求偏导数,寻找模型参数的最大似然估计,如下式所示:
求解上式,求得模型参数:
对于变压器油箱表面测得的混合信号自回归模型误差参数的最大似然估计,用下式代替:
上式中,下标k表示采集得到的混合信号中的第k个信号。
5.根据权利要求1所述的电力变压器有载分接开关振动脉冲数自动识别方法,其特征在于:所述步骤(4)中,对于第j次分段,通过下式计算信号回归模型误差的信息熵A值:
上式中,C=(M1+M2-N)(log2π+1),为常数。
6.根据权利要求1所述的电力变压器有载分接开关振动脉冲数自动识别方法,其特征在于:所述步骤(5)中,采用MINMAX算法获取信息熵A值的最小值和最大值。
7.根据权利要求1所述的电力变压器有载分接开关振动脉冲数自动识别方法,其特征在于:所述的步骤(6)中,采用DIFF算法求取信号长度区间对应的信息熵A值的梯度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710146756.0A CN107016352B (zh) | 2017-03-13 | 2017-03-13 | 电力变压器有载分接开关振动脉冲数自动识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710146756.0A CN107016352B (zh) | 2017-03-13 | 2017-03-13 | 电力变压器有载分接开关振动脉冲数自动识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107016352A CN107016352A (zh) | 2017-08-04 |
CN107016352B true CN107016352B (zh) | 2019-06-18 |
Family
ID=59440679
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710146756.0A Expired - Fee Related CN107016352B (zh) | 2017-03-13 | 2017-03-13 | 电力变压器有载分接开关振动脉冲数自动识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107016352B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102018105087B3 (de) * | 2018-03-06 | 2019-06-27 | Maschinenfabrik Reinhausen Gmbh | System, Verfahren und Computerprogrammprodukt zur Erfassung physikalischer Größen von mindestens einer Komponente eines Stufentransformators und zum Monitoring der Komponenten eines Stufentransformators |
CN109270441B (zh) * | 2018-07-17 | 2021-03-19 | 国网浙江省电力有限公司检修分公司 | 基于振动信号的高压断路器分闸特性参数在线检测方法 |
CN110702215B (zh) * | 2019-10-19 | 2021-04-06 | 广东石油化工学院 | 一种利用回归树的变压器运行状态振声检测方法和系统 |
CN112257640A (zh) * | 2020-10-30 | 2021-01-22 | 国网上海市电力公司 | 基于Hibert变换和信息熵的脉冲特征参数提取方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103308156A (zh) * | 2013-05-16 | 2013-09-18 | 国家电网公司 | 电力变压器有载调压开关及本体振动信号的分离方法 |
CN104215324A (zh) * | 2014-08-18 | 2014-12-17 | 国家电网公司 | 电力变压器有载分接开关振动信号提取方法 |
CN204988500U (zh) * | 2015-08-20 | 2016-01-20 | 国网天津武清供电有限公司 | 一种有载分接开关机械振动信号调理电路 |
CN105510812A (zh) * | 2015-10-12 | 2016-04-20 | 国家电网公司 | 一种换流变分接开关状态监测系统及非平稳信号频谱估计方法 |
CN106840379A (zh) * | 2017-03-08 | 2017-06-13 | 潘小胜 | 一种机械振动信号的故障分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009120765A1 (en) * | 2008-03-25 | 2009-10-01 | Abb Research Ltd. | Method and apparatus for analyzing waveform signals of a power system |
-
2017
- 2017-03-13 CN CN201710146756.0A patent/CN107016352B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103308156A (zh) * | 2013-05-16 | 2013-09-18 | 国家电网公司 | 电力变压器有载调压开关及本体振动信号的分离方法 |
CN104215324A (zh) * | 2014-08-18 | 2014-12-17 | 国家电网公司 | 电力变压器有载分接开关振动信号提取方法 |
CN204988500U (zh) * | 2015-08-20 | 2016-01-20 | 国网天津武清供电有限公司 | 一种有载分接开关机械振动信号调理电路 |
CN105510812A (zh) * | 2015-10-12 | 2016-04-20 | 国家电网公司 | 一种换流变分接开关状态监测系统及非平稳信号频谱估计方法 |
CN106840379A (zh) * | 2017-03-08 | 2017-06-13 | 潘小胜 | 一种机械振动信号的故障分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107016352A (zh) | 2017-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107016352B (zh) | 电力变压器有载分接开关振动脉冲数自动识别方法 | |
CN107305774B (zh) | 语音检测方法和装置 | |
CN101430360B (zh) | 供电系统中二次设备的错误数据鉴别方法 | |
Wang et al. | A new method with Hilbert transform and slip-SVD-based noise-suppression algorithm for noisy power quality monitoring | |
US20160322064A1 (en) | Method and apparatus for signal extraction of audio signal | |
KR101505341B1 (ko) | 로우 비트 레이트 신호 코더 및 디코더 | |
CA3190353A1 (en) | Method and apparatus for compressing and decompressing a higher order ambisonics representation | |
CN106199267B (zh) | 一种电力设备故障特征分析方法 | |
CN109932053B (zh) | 一种用于高压并联电抗器的状态监测装置及方法 | |
CN112014047A (zh) | 一种有载分接开关机械故障诊断方法 | |
CN109541455A (zh) | 一种基于s变换时频谱svd降噪的oltc冲击特性提取方法 | |
CN112345213A (zh) | 一种低压直流断路器机械故障诊断方法 | |
CN108872666B (zh) | 一种负载开关事件检测方法及系统 | |
CN1232834C (zh) | 基于振动分析的真空断路器触头关合时刻的在线检测方法 | |
CN110658445B (zh) | 一种有载分接开关机械故障的分析诊断方法 | |
CN102279364A (zh) | 一种基于music与psa的笼型异步电动机转子断条故障检测方法 | |
Konstantinopoulos et al. | Synchrophasor missing data recovery via data-driven filtering | |
CN110207974A (zh) | 基于振动信号时频能量分布特征的断路器故障识别方法 | |
CN107703433A (zh) | 基于异或逻辑的hbsm‑mmc内部igbt开路故障识别方法 | |
CN106501602A (zh) | 一种基于滑窗频谱分离的基波参数测量方法 | |
CN111079647A (zh) | 一种断路器缺陷识别方法 | |
CN106053936B (zh) | 一种获取电学信号瞬时频率的方法及系统 | |
Zhang et al. | Vibration signal processing and state analysis technology for on-load tap-changer | |
Wen-Xi et al. | Automatic segmentation method for voltage sag detection and characterization | |
Khurram et al. | Detection of oscillatory modes in power systems using empirical wavelet transform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190618 |