CN103116111B - 一种电力变压器绕组工作状态的诊断方法 - Google Patents

一种电力变压器绕组工作状态的诊断方法 Download PDF

Info

Publication number
CN103116111B
CN103116111B CN201310021637.4A CN201310021637A CN103116111B CN 103116111 B CN103116111 B CN 103116111B CN 201310021637 A CN201310021637 A CN 201310021637A CN 103116111 B CN103116111 B CN 103116111B
Authority
CN
China
Prior art keywords
signal
transformer winding
kth
vibration signal
winding
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.)
Active
Application number
CN201310021637.4A
Other languages
English (en)
Other versions
CN103116111A (zh
Inventor
李红雷
黄华
林一
王丰华
耿超
金之俭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
State Grid Corp of China SGCC
Shanghai Municipal Electric Power Co
East China Power Test and Research Institute Co Ltd
Original Assignee
Shanghai Jiaotong University
State Grid Corp of China SGCC
Shanghai Municipal Electric Power Co
East China Power Test and Research Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Jiaotong University, State Grid Corp of China SGCC, Shanghai Municipal Electric Power Co, East China Power Test and Research Institute Co Ltd filed Critical Shanghai Jiaotong University
Priority to CN201310021637.4A priority Critical patent/CN103116111B/zh
Publication of CN103116111A publication Critical patent/CN103116111A/zh
Application granted granted Critical
Publication of CN103116111B publication Critical patent/CN103116111B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)

Abstract

本发明涉及一种电力变压器绕组工作状态的诊断方法,包括以下步骤:(1)向激振器输入信号Vi对变压器绕组进行激振,通过变压器绕组表面的N个振动加速度传感器采集和记录各个测点的振动信号Voi;(2)分别对输入信号Vi和各个测点的振动信号Voi进行傅里叶变换,得到变压器绕组的振动频响曲线H(ω);(3)对振动频响曲线H(ω)进行傅里叶反变换,得到变压器绕组的自由振动信号H(t);(4)将自由振动信号H(t)分解为若干个固有模式分量;(5)将分解得到的全部固有模式分量进行希尔伯特变换;(6)构造解析信号;(7)得到变压器绕组的第k阶固有频率;(8)根据固有频率的变化对变压器绕组状态进行判别。与现有技术相比,本发明具有诊断准确、操作便捷、安全等优点。

Description

一种电力变压器绕组工作状态的诊断方法
技术领域
本发明涉及一种基于信号监测的诊断方法,尤其是涉及一种电力变压器绕组工作状态的诊断方法。
背景技术
变压器是电力系统中最重要的设备之一,其运行的稳定性对电力系统安全影响重大。随着我国电网容量的日益增大,短路容量亦随之不断增大,变压器出口短路形成的冲击电流产生的巨大电磁作用力对变压器绕组的机械强度和动稳定性构成了严重的威胁。目前变电站设备及线路的运行环境始终不容乐观,因外部短路造成变压器绕组受冲击而引发的变形是变压器运行过程中较为常见的故障,对系统的安全运行造成了很大的威胁。
变压器遭受突发短路后,其绕组可能首先发生松动或轻微变形,通过大量的实验研究分析变压器绕组变形具有累积效应,如果对于松动或变形不能及时发现和修复,那么在变压器的松动或变形累积到一定程度后会使变压器的抗短路能力大幅下降而在遭受较小的冲击电流下也会引发大的事故发生。绕组的变形一方面会导致机械抗短路电流冲击能力的下降,另一方面也会导致线圈内部局部绝缘距离发生变化,使局部出现绝缘薄弱点,当遇到过电压作用时,绕组有可能发生饼间或匝间短路导致变压器绝缘击穿事故,或者由于局部场强增大而引起局部放电,绝缘损伤部位会逐渐扩大,最终导致变压器发生绝缘击穿事故而引发进一步的事态扩大。
因此,在运行过程中当变压器经历了外部短路事故后或运行一段时间后的常规检修中,如何有效地检测出变压器绕组是否存在松动和变形,从而判断变压器是否需要检修处理显得十分重要,是保障变压器安全运行的一个重要手段,因此变压器绕组变形的检测是目前变压器常规试验项目之一。
目前实际应用的对变压器绕组状态的检测方法主要有以下三种:
1、短路阻抗法
变压器短路阻抗是当负载阻抗为零时变压器内部的等效阻抗,短路阻抗是变压器绕组的漏抗和电阻的矢量和,由于变压器直流电阻相对于漏抗数值很小,因此变压器的短路阻抗反映的主要是变压器绕组的漏抗。由变压器的理论分析可知,变压器漏抗值是由绕组的几何尺寸所决定的,或者说是由绕组的结构决定的,一旦变压器绕组发生变形,从理论上来说变压器的漏抗相应也会发生变化,因此通过对变压器短路阻抗的检测可以间接地反映变压器绕组内部是否发生了变形。
一般情况下,运行中的变压器受到了短路电流的冲击后,或在定期常规检查时要将测得的短路阻抗值与原有的记录进行比较来判断绕组是否发生了变形,如果短路阻抗值变化较大,例如国标中设定为变化超过3%,则可确认绕组有显著变形。
按照有关标准规定,变压器在短路阻抗测试试验中,要求测量每一相的短路阻抗,并把试验后所测量的短路阻抗值与以往试验的数据加以比较,根据其变化的程度,作为判断被试变压器绕组是否合格的重要依据之一。
从实际应用情况来看,短路阻抗法在长期的生产实践中已建立了标准,判据较为明确,在国际电工标准IEC60076-5和GB1095-85中均明确给出了线圈变形程度的判据。但很多情况下这种方法的灵敏度很低,故障的检出率较低,只有在线圈整体变形情况较为严重时才能够得到较明确的反映。
2、频响分析法
频响分析法的基本原理是将变压器绕组视为一个分布参数网络,它由对地电容C、纵向电容K、电感L等分布参数构成一无源线性双端口网络,该网络的特性在频域上可以用传递函数H(jω)来描述。
绕组发生局部机械变形后,其内部的分布电感L、纵向电容K和对地电容C等分布参数会发生相应的变化,从而在网络的传递函数上得到反映。因此分析变压器绕组的网络传递函数曲线的变化情况就可以分析内部的网络电参数是否发生变化,从而推断相应的机械结构是否发生了变形,这是频响分析法测试变压器绕组变形的依据和基础。
频响法测试首先将一稳定的正弦扫频电压信号Vi施加到被试变压器绕组的一端,然后同时记录该端口和其它输出端口上的电压Vo,从而得到被试绕组的一组频响特性曲线,其表达式为
H(jω)=Vo/Vi
频响法的测试灵敏度较短路阻抗法高,但由于其频响波形的复杂性,对绕组状况的判别需要较多的经验,较难形成明确的定量判据,因此至今没有形成判别标准。
上述两种方法是目前判别变压器绕组状况最常用的,两种方法都是采用电测方法,出发点都是基于变压器绕组发生明显变形的状况下模型中对应的元件电参数发生变化来进行测量判别,这对变压器绕组发生较明显的变形情况较为适宜,但对绕组发生轻微变形,尤其是对变压器绕组存在的相对松动和扭曲变形的状态不能给出较明确的判断,因为这些情况下反映在等效电路模型中的电参数几乎没有变化,其传递函数的变化也就非常小。然而变压器绕组松动或扭曲变形对其抗短路能力有很大的影响,因此研究绕组的状况需有灵敏度更高的方法来进行判别。
3、振动分析法
振动分析法的基本原理是把变压器绕组看作一个机械结构体,则当绕组结构或受力发生任何变化时,都可以从它的机械振动特性变化上得到反映。因此,可通过分析箱壁上的振动信号来对绕组的工作状态进行检测。与前述方法相比较,振动分析法的最大优点是可通过吸附在变压器箱壁上的振动传感器来获得变压器的振动信号,通过分析振动特性的变化来判断绕组状态的变化情况,只要绕组的机械特性(如结构变形、预紧力松动等)发生变化,都可以从它的机械振动特性变化上得到反映,从而大大提高了检测的灵敏度。此外,将振动传感器置于箱壁上的振动检测与整个强电系统没有直接的连接,对于整个电气系统的正常运行没有任何影响,因此,可发展成为一种较准确、便捷、安全的在线监测方法。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种诊断准确、操作便捷、安全的电力变压器绕组工作状态的诊断方法。
本发明的目的可以通过以下技术方案来实现:
一种电力变压器绕组工作状态的诊断方法,所述的变压器绕组上连接有激振器,且变压器绕组表面上放置有N个振动加速度传感器,该方法包括以下步骤:
(1)向激振器输入信号Vi对变压器绕组进行激振,通过变压器绕组表面的N个振动加速度传感器采集和记录各个测点的振动信号Voi(i=1,2,…,N);
(2)分别对输入信号Vi和各个测点的振动信号Voi进行傅里叶变换,得到变压器绕组的振动频响曲线H(ω):
H ( ω ) = ( Σ i = 1 N V oi ( ω ) ) / V i ( ω )
式中,Vi(ω)为输入信号Vi的傅里叶变换,为Voi(ω)各个测点振动信号Voi的傅里叶变换;
(3)对振动频响曲线H(ω)进行傅里叶反变换,得到变压器绕组的自由振动信号H(t);
(4)将自由振动信号H(t)分解为若干个固有模式分量;
(5)将分解得到的全部固有模式分量通过以下公式进行希尔伯特变换:
d k ( t ) = 1 π ∫ - ∞ + ∞ c k ( t ) t - τ dτ , k = 1 , . . . , p
式中,p为分解得到的固有模式分量的数目,dk(t)为第k个固有模式分量的希尔伯特变换值,ck(t)为分解得到的第k个固有模式分量,t为时间;
(6)根据以下公式构造解析信号:
Z k ( t ) = c k ( t ) + j d k ( t ) = a k ( t ) e j θ k ( t )
式中,Zk(t)为根据第k个固有模式分量及其希尔伯特变换值得到的第k个解析信号,ak(t)为幅值函数,θk(t)为相位函数,且有下述关系式
θ k ( t ) = arctan [ d k ( t ) c k ( t ) ] = ω dk t + φ k
式中,Ak为第k个解析信号的幅值;为第k个解析信号的阻尼比;ωk为第k个解析信号对应的固有频率;ωdk为第k个解析信号的固有阻尼频率,且有
(7)计算第k个相位函数θk(t)的导数,得到变压器绕组的第k阶固有频率;
(8)根据固有频率的变化对变压器绕组状态进行判别:当变压器绕组的各阶固有频率向低频方向偏移且数值减小为原来的5%及以上时,判定变压器绕组发生松动或变形。
所述的步骤(4)具体包括以下步骤:
4a)对自由振动信号H(t)求导,得到时间序列y(t);
4b)计算时间序列y(t)相邻两点的乘积
pyi(t)=yi(t)×yi-1(t)
其中,i=2,3,…,n-1,n为自由振动信号的点数;
4c)根据乘积pyi(t)和时间序列y(t)的正负,依次找寻自由振动信号H(t)的所有局部极大值点eb(t)和所有局部极小值点es(t):
当pyi(t)<0时,若yi-1(t)<0,则Hi-1(t)为局部极小值点;若yi-1(t)>0,则Hi-1(t)为局部极大值点;
当pyi(t)>0时,Hi-1(t)为非极值点;
当pyi(t)=0时,若yi-1(t)=0,计算两点yi(t)和yi-2(t)的乘积,令pyi(t)'=yi(t)×yi-2(t),若pyi(t)'<0且yi-2(t)<0,则Hi-1(t)为局部极小值点;若pyi(t)'<0且yi-2(t)>0,则Hi-1(t)为局部极大值点;若yi-2(t)=0,则Hi-1(t)为非极值点;
4d)将所述所有局部极大值点eb(t)和所有局部极小值点es(t)用三次样条插值函数s(t)连接起来分别求出上包络线emax(t)和下包络线emin(t);
4e)根据求得的上包络线emax(t)和下包络线emin(t)计算上、下包络线的均值m(t)=(emax(t)+emin(t))/2,将自由振动信号H(t)减去m(t),得到一个新的时间序列Y(t);
4f)判断上述时间序列Y(t)是否同时满足下述两个条件:
A.在整个信号长度上,极值点和过零点的数目必须相等或者至多只相差一个;
B.在任意时刻,由极大值点定义的上包络线和由极小值点定义的下包络线的平均值为零;
若同时满足上述两个条件,则Y(t)为固有模式分量;若不能同时满足上述两个条件,则将Y(t)作为一个原始分量,重复前述步骤4a)~步骤4e),直到Y(t)同时满足上述两个条件,将Y(t)记为ci(t),则ci(t)为自由振动信号H(t)的一个固有模式分量,i=1,2,…,p;
4g)将ci(t)从自由振动信号H(t)中分离出来,得到差值信号ri(t),将差值信号ri(t)作为待处理的振动信号H1(t),ri(t)=H(t)-ci(t);
4h)重复上述步骤4a)~步骤4g),直至满足迭代终止准则,得到全部p个固有模式分量,所述的迭代终止准则为:得到的新的时间序列Y(t)为窄带信号。
所述的步骤4d)中的三次样条插值函数s(t)是在自由振动信号H(t)的每一个小区间[ti,ti+1]上不超过三次的多项式,i=1,2,…,n-1,其表达式为
s ( t ) = m i ( t i + 1 - t ) 3 6 ( t i + 1 - t i ) + m i + 1 ( t - t i ) 3 6 ( t i + 1 - t i ) + x i + 1 ( t ) - x i ( t ) t i + 1 - t i - t i + 1 - t i 6 ( m i + 1 - m i ) + x i ( t ) - m i ( t i + 1 - t i ) 2 6
式中,mi和mi+1为三次样条插值函数s(t)在区间[ti,ti+1]两端点处的二阶导数值。
所述的输入信号Vi为白噪声信号。
与现有技术相比,本发明通过分析变压器绕组的振动特性获得固有频率的变化,进而地变压器绕组的工作状态进行诊断,具有诊断准确、操作便捷、安全等优点。
附图说明
图1为本发明的流程示意图;
图2为实施例中变压器绕组状态良好时的振动频响曲线;
图3为实施例中变压器绕组状态恶化时的振动频响曲线;
图4为实施例中将自由振动信号分解为4个IMF分量的结果示意图;
图5为实施例中变压器绕组状态良好时的4个IMF的相位函数曲线示意图;
图6为实施例中变压器绕组状态恶化时的4个IMF的相位函数曲线示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例
如图1所示,一种电力变压器绕组工作状态的诊断方法,将其应用于一10kV的变压器绕组进行监测诊断,具体包括以下步骤:
(1)在该变压器绕组上放置激振器,将经过功率放大器放大的20kHz的白噪声信号Vi输入激振器对变压器绕组进行激振,在该变压器绕组表面放置20个振动加速度传感器采集和记录各个测点的振动信号Voi(i=1,2,…,20),采集时间为0.04s。
(2)分别对输入的白噪声信号Vi和各个测点的振动信号Voi进行傅里叶变换,得到变压器绕组的振动频响曲线H(ω):
H ( &omega; ) = ( &Sigma; i = 1 N V oi ( &omega; ) ) / V i ( &omega; )
式中,Vi(ω)为输入的白噪声信号的傅里叶变换;Voi(ω)为各个测点振动信号的傅里叶变换;本实施例中,N=20,图2显示了该变压器绕组状态良好时的振动频响曲线,图3显示了该变压器绕组状态恶化时的振动频响曲线。
(3)对振动频响曲线H(ω)进行傅里叶反变换,得到变压器绕组的自由振动信号H(t);
(4)根据以下步骤将自由振动信号H(t)分解为若干个固有模式分量(IntrinsicMode Function,简称为IMF分量):
4a)对自由振动信号H(t)求导,得到时间序列y(t);
4b)计算时间序列y(t)相邻两点的乘积
pyi(t)=yi(t)×yi-1(t)
其中,i=2,3,…,n-1,n为自由振动信号的点数;
4c)根据乘积pyi(t)和时间序列y(t)的正负,依次找寻自由振动信号H(t)的所有局部极大值点eb(t)和所有局部极小值点es(t):
当pyi(t)<0时,若yi-1(t)<0,则Hi-1(t)为局部极小值点;若yi-1(t)>0,则Hi-1(t)为局部极大值点;
当pyi(t)>0时,Hi-1(t)为非极值点;
当pyi(t)=0时,若yi-1(t)=0,计算两点yi(t)和yi-2(t)的乘积,令pyi(t)'=yi(t)×yi-2(t),若pyi(t)'<0且yi-2(t)<0,则Hi-1(t)为局部极小值点;若pyi(t)'<0且yi-2(t)>0,则Hi-1(t)为局部极大值点;若yi-2(t)=0,则Hi-1(t)为非极值点;
4d)将所述所有局部极大值点eb(t)和所有局部极小值点es(t)用三次样条插值函数s(t)连接起来分别求出上包络线emax(t)和下包络线emin(t);
所述的三次样条插值函数s(t)是在自由振动信号H(t)的每一个小区间[ti,ti+1]上不超过三次的多项式,i=1,2,…,n-1,其表达式为
s ( t ) = m i ( t i + 1 - t ) 3 6 ( t i + 1 - t i ) + m i + 1 ( t - t i ) 3 6 ( t i + 1 - t i ) + x i + 1 ( t ) - x i ( t ) t i + 1 - t i - t i + 1 - t i 6 ( m i + 1 - m i ) + x i ( t ) - m i ( t i + 1 - t i ) 2 6
式中,mi和mi+1为三次样条插值函数s(t)在区间[ti,ti+1]两端点处的二阶导数值;
4e)根据求得的上包络线emax(t)和下包络线emin(t)计算上、下包络线的均值m(t)=(emax(t)+emin(t))/2,将自由振动信号H(t)减去m(t),得到一个新的时间序列Y(t);
4f)判断上述时间序列Y(t)是否同时满足下述两个条件:
A.在整个信号长度上,极值点和过零点的数目必须相等或者至多只相差一个;
B.在任意时刻,由极大值点定义的上包络线和由极小值点定义的下包络线的平均值为零;
若同时满足上述两个条件,则Y(t)为固有模式分量;若不能同时满足上述两个条件,则将Y(t)作为一个原始分量,重复前述步骤4a)~步骤4e),直到Y(t)同时满足上述两个条件,将Y(t)记为ci(t),则ci(t)为自由振动信号H(t)的一个固有模式分量,i=1,2,…,p;
4g)将ci(t)从自由振动信号H(t)中分离出来,得到差值信号ri(t),将差值信号ri(t)作为待处理的振动信号H1(t),ri(t)=H(t)-ci(t);
4h)重复上述步骤4a)~步骤4g),直至满足迭代终止准则,得到全部p个固有模式分量,所述的迭代终止准则为:得到的新的时间序列Y(t)为窄带信号,窄带信号是指信号的带宽Δf远小于中心频率fC的信号。
经过上述步骤,最初的自由振动信号H(t)已经被分解为4个IMF分量,图4显示了将自由振动信号分解为4个IMF分量的结果,其中,图(4a)为分量IMF1的示意图,图(4b)为分量IMF2的示意图,图(4c)为分量IMF3的示意图,图(4d)为分量IMF4的示意图。
(5)将分解得到的全部固有模式分量通过以下公式进行希尔伯特变换:
d k ( t ) = 1 &pi; &Integral; - &infin; + &infin; c k ( t ) t - &tau; d&tau; , k = 1 , . . . , p
式中,p为分解得到的固有模式分量的数目,dk(t)为第k个固有模式分量的希尔伯特变换值,ck(t)为分解得到的第k个固有模式分量,t为时间。
(6)根据以下公式构造解析信号:
Z k ( t ) = c k ( t ) + j d k ( t ) = a k ( t ) e j &theta; k ( t )
式中,Zk(t)为根据第k个固有模式分量及其希尔伯特变换值得到的第k个解析信号,ak(t)为幅值函数,θk(t)为相位函数,且有下述关系式
&theta; k ( t ) = arctan [ d k ( t ) c k ( t ) ] = &omega; dk t + &phi; k
式中,Ak为第k个解析信号的幅值;为第k个解析信号的阻尼比;ωk为第k个解析信号对应的固有频率;ωdk为第k个解析信号的固有阻尼频率,且有
(7)计算第k个相位函数θk(t)的导数,得到变压器绕组的第k阶固有频率。图5显示了该变压器绕组状态良好时的4个IMF的相位函数曲线,图(5a)~(5d)分别为分量IMF1~分量IMF4的相位函数曲线;图6显示了该变压器绕组状态恶化时的4个IMF的相位函数曲线,图(6a)~(6d)分别为分量IMF1~分量IMF4的相位函数曲线。
(8)根据固有频率的变化对变压器绕组状态进行判别:当变压器绕组的各阶固有频率向低频方向偏移且数值减小为原来的5%及以上时,判定变压器绕组发生松动或变形,此时需要及时进行处理,避免形成重大故障。

Claims (4)

1.一种电力变压器绕组工作状态的诊断方法,所述的变压器绕组上连接有激振器,且变压器绕组表面上放置有N个振动加速度传感器,其特征在于,该方法包括以下步骤:
(1)向激振器输入信号Vi对变压器绕组进行激振,通过变压器绕组表面的N个振动加速度传感器采集和记录各个测点的振动信号Voi,其中,i=1,2,…,N;
(2)分别对输入信号Vi和各个测点的振动信号Voi进行傅里叶变换,得到变压器绕组的振动频响曲线H(ω):
H ( &omega; ) = ( &Sigma; i = 1 N V oi ( &omega; ) ) / V i ( &omega; )
式中,Vi(ω)为输入信号Vi的傅里叶变换,为Voi(ω)各个测点振动信号Voi的傅里叶变换;
(3)对振动频响曲线H(ω)进行傅里叶反变换,得到变压器绕组的自由振动信号H(t);
(4)将自由振动信号H(t)分解为若干个固有模式分量;
(5)将分解得到的全部固有模式分量通过以下公式进行希尔伯特变换:
d k ( t ) = 1 &pi; &Integral; - &infin; + &infin; c k ( t ) t - &tau; d&tau; , k = 1 , &CenterDot; &CenterDot; &CenterDot; , p
式中,p为分解得到的固有模式分量的数目,dk(t)为第k个固有模式分量的希尔伯特变换值,ck(t)为分解得到的第k个固有模式分量,t为时间;
(6)根据以下公式构造解析信号:
Z k ( t ) = c k ( t ) + jd k ( t ) = a k ( t ) e j&theta; k ( t )
式中,Zk(t)为根据第k个固有模式分量及其希尔伯特变换值得到的第k个解析信号,ak(t)为幅值函数,θk(t)为相位函数,且有下述关系式
&theta; k ( t ) = arctan [ d k ( t ) c k ( t ) ] = &omega; dk t + &phi; k
式中,Ak为第k个解析信号的幅值;为第k个解析信号的阻尼比;ωk为第k个解析信号对应的固有频率;ωdk为第k个解析信号的固有阻尼频率,且有
(7)计算第k个相位函数θk(t)的导数,得到变压器绕组的第k阶固有频率;
(8)根据固有频率的变化对变压器绕组状态进行判别:当变压器绕组的各阶固有频率向低频方向偏移且数值减小为原来的5%及以上时,判定变压器绕组发生松动或变形。
2.根据权利要求1所述的一种电力变压器绕组工作状态的诊断方法,其特征在于,所述的步骤(4)具体包括以下步骤:
4a)对自由振动信号H(t)求导,得到时间序列y(t);
4b)计算时间序列y(t)相邻两点的乘积
pyi(t)=yi(t)×yi-1(t)
i=2,3,…,n-1,其中,n为自由振动信号的点数;
4c)根据乘积pyi(t)和时间序列y(t)的正负,依次找寻自由振动信号H(t)的所有局部极大值点eb(t)和所有局部极小值点es(t):
当pyi(t)<0时,若yi-1(t)<0,则Hi-1(t)为局部极小值点;若yi-1(t)>0,则Hi-1(t)为局部极大值点;
当pyi(t)>0时,Hi-1(t)为非极值点;
当pyi(t)=0时,若yi-1(t)=0,计算两点yi(t)和yi-2(t)的乘积,令pyi(t)'=yi(t)×yi-2(t),若pyi(t)'<0且yi-2(t)<0,则Hi-1(t)为局部极小值点;若pyi(t)'<0且yi-2(t)>0,则Hi-1(t)为局部极大值点;若yi-2(t)=0,则Hi-1(t)为非极值点;
4d)将所述所有局部极大值点eb(t)和所有局部极小值点es(t)用三次样条插值函数s(t)连接起来分别求出上包络线emax(t)和下包络线emin(t);
4e)根据求得的上包络线emax(t)和下包络线emin(t)计算上、下包络线的均值m(t)=(emax(t)+emin(t))/2,将自由振动信号H(t)减去m(t),得到一个新的时间序列Y(t);
4f)判断上述时间序列Y(t)是否同时满足下述两个条件:
A.在整个信号长度上,极值点和过零点的数目必须相等或者至多只相差一个;
B.在任意时刻,由极大值点定义的上包络线和由极小值点定义的下包络线的平均值为零;
若同时满足上述两个条件,则Y(t)为固有模式分量;若不能同时满足上述两个条件,则将Y(t)作为一个原始分量,重复前述步骤4a)~步骤4e),直到Y(t)同时满足上述两个条件,将Y(t)记为ci(t),则ci(t)为自由振动信号H(t)的一个固有模式分量,i=1,2,…,p;
4g)将ci(t)从自由振动信号H(t)中分离出来,得到差值信号ri(t),将差值信号ri(t)作为待处理的振动信号H1(t),ri(t)=H(t)-ci(t);
4h)重复上述步骤4a)~步骤4g),直至满足迭代终止准则,得到全部p个固有模式分量,所述的迭代终止准则为:得到的新的时间序列Y(t)为窄带信号。
3.根据权利要求2所述的一种电力变压器绕组工作状态的诊断方法,其特征在于,所述的步骤4d)中的三次样条插值函数s(t)是在自由振动信号H(t)的每一个小区间[ti,ti+1]上不超过三次的多项式,i=1,2,…,n-1,其表达式为
s ( t ) = m i ( t i + 1 - t ) 3 6 ( t i + 1 - t i ) + m i + 1 ( t - t i ) 3 6 ( t i + 1 - t i ) + x i + 1 ( t ) - x i ( t ) t i + 1 - t i - t i + 1 t i 6 ( m i + 1 - m i ) + x i ( t ) - m i ( t i + 1 - t i ) 2 6
式中,mi和mi+1为三次样条插值函数s(t)在区间[ti,ti+1]两端点处的二阶导数值。
4.根据权利要求1所述的一种电力变压器绕组工作状态的诊断方法,其特征在于,所述的输入信号Vi为白噪声信号。
CN201310021637.4A 2013-01-21 2013-01-21 一种电力变压器绕组工作状态的诊断方法 Active CN103116111B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310021637.4A CN103116111B (zh) 2013-01-21 2013-01-21 一种电力变压器绕组工作状态的诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310021637.4A CN103116111B (zh) 2013-01-21 2013-01-21 一种电力变压器绕组工作状态的诊断方法

Publications (2)

Publication Number Publication Date
CN103116111A CN103116111A (zh) 2013-05-22
CN103116111B true CN103116111B (zh) 2014-12-24

Family

ID=48414537

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310021637.4A Active CN103116111B (zh) 2013-01-21 2013-01-21 一种电力变压器绕组工作状态的诊断方法

Country Status (1)

Country Link
CN (1) CN103116111B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103176104A (zh) * 2013-02-04 2013-06-26 上海交通大学 一种变压器绕组运行状态的诊断方法
CN103454520A (zh) * 2013-08-02 2013-12-18 国家电网公司 基于在线频率响应法的变压器绕组变形在线监测方法
CN103454526B (zh) * 2013-08-23 2016-09-07 上海交通大学 一种基于klc模型的电力变压器绕组故障类型判定方法
CN104345233A (zh) * 2014-10-22 2015-02-11 国家电网公司 一种用于配电变压器故障检测的在线监测终端及其检测方法
CN105203914B (zh) * 2015-07-22 2018-09-04 广东电网有限责任公司电力科学研究院 一种突发短路下的变压器绕组状态诊断方法
CN105093059B (zh) * 2015-08-26 2018-12-11 云南电网有限责任公司电力科学研究院 一种基于归一化奇异谱熵的变压器绕组工作状态检测方法
CN105352471A (zh) * 2015-11-30 2016-02-24 云南电网有限责任公司电力科学研究院 一种变压器绕组状态诊断方法
CN105627904B (zh) * 2016-02-01 2018-08-07 国网浙江省电力公司电力科学研究院 一种变压器绕组变形的判定方法
CN107239739B (zh) * 2017-05-05 2020-10-27 北京化工大学 一种尺度参数控制可调的信号包络线提取方法
CN108917906B (zh) * 2018-04-02 2019-08-02 西南交通大学 一种牵引变压器绕组变形故障检测方法
CN109029959B (zh) * 2018-08-27 2020-08-21 深圳供电局有限公司 一种变压器绕组的机械状态检测方法
CN110133146A (zh) * 2019-05-28 2019-08-16 国网上海市电力公司 一种考虑不平衡数据样本的变压器故障诊断方法及系统
CN110632477A (zh) * 2019-11-02 2019-12-31 广东石油化工学院 一种利用Hilbert空间因子的变压器运行状态振声检测方法和系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4897794A (en) * 1987-07-21 1990-01-30 Egg Co., Ltd. Impulse coil tester
JP2000097982A (ja) * 1998-09-21 2000-04-07 Ikd:Kk コイル試験評価装置
CN1776441A (zh) * 2005-11-30 2006-05-24 上海市电力公司 利用扫频电源激振检测变压器绕组状态的装置
CN101556299A (zh) * 2009-05-18 2009-10-14 重庆电力科学试验研究院 检测变压器绕组微小形变的方法
CN101738567A (zh) * 2008-11-25 2010-06-16 上海市电力公司 利用恒流扫频电源激振检测变压器绕组状态的系统和方法
CN101937047A (zh) * 2009-06-30 2011-01-05 上海市电力公司 利用振动波形检测变压器绕组状态的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4897794A (en) * 1987-07-21 1990-01-30 Egg Co., Ltd. Impulse coil tester
JP2000097982A (ja) * 1998-09-21 2000-04-07 Ikd:Kk コイル試験評価装置
CN1776441A (zh) * 2005-11-30 2006-05-24 上海市电力公司 利用扫频电源激振检测变压器绕组状态的装置
CN101738567A (zh) * 2008-11-25 2010-06-16 上海市电力公司 利用恒流扫频电源激振检测变压器绕组状态的系统和方法
CN101556299A (zh) * 2009-05-18 2009-10-14 重庆电力科学试验研究院 检测变压器绕组微小形变的方法
CN101937047A (zh) * 2009-06-30 2011-01-05 上海市电力公司 利用振动波形检测变压器绕组状态的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
振动法在线监测突发短路时变压器绕组状态;王丰华 等;《控制工程》;20110731;第18卷(第4期);第596-599页 *
振动频响分析法在大型变压器绕组变形检测中的应用;邵宇鹰 等;《上海交通大学学报》;20100930;第44卷(第9期);第1223-1228页 *

Also Published As

Publication number Publication date
CN103116111A (zh) 2013-05-22

Similar Documents

Publication Publication Date Title
CN103116111B (zh) 一种电力变压器绕组工作状态的诊断方法
CN103176104A (zh) 一种变压器绕组运行状态的诊断方法
CN102998545B (zh) 一种变压器绕组工作状态的在线监测方法
CN103822696B (zh) 一种变压器绕组状态的诊断方法
CN105182172B (zh) 基于振动信号形态谱的变压器突发短路下绕组状态诊断方法
CN101937047A (zh) 利用振动波形检测变压器绕组状态的方法
US9404957B2 (en) Fault diagnosis and preliminary location system and method for transformer core looseness
CN104236702B (zh) 电力变压器内部松动判断系统及其方法
CN202735425U (zh) 基于振动的电力变压器故障检测系统
CN100373166C (zh) 利用扫频电源激振检测变压器绕组状态的装置
CN102998544A (zh) 变压器短路时绕组工作状态诊断方法
CN206114822U (zh) 一种电力变压器绕组形变状态多信息检测装置
CN105092024B (zh) 一种电力变压器绕组轴向压紧状态在线监测方法及系统
CN105203914B (zh) 一种突发短路下的变压器绕组状态诊断方法
CN105093059B (zh) 一种基于归一化奇异谱熵的变压器绕组工作状态检测方法
CN101738567A (zh) 利用恒流扫频电源激振检测变压器绕组状态的系统和方法
CN106019107B (zh) 一种空心电抗器局部放电的检测方法及系统
CN103968939B (zh) 基于平均位移法的变压器绕组松动故障检测方法
CN104061851A (zh) 基于过电压响应的变压器绕组变形在线监测方法
WO2016065959A1 (zh) 中性点不接地的10kV系统中铁磁谐振的诊断方法
CN102081136A (zh) 一种适用于现场的冲击电压下gis局部放电检测方法
CN104655967B (zh) 配电变压器绕组振动信号特征量提取方法
CN205091401U (zh) 一种变压器绕组振动频响测试系统
CN205300880U (zh) 一种变压器绕组振动法测试系统
CN104237727A (zh) 一种变压器近区短路信号监测装置及短路录波分析方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant