CN103259251B - 一种基于权重数学形态学的变压器励磁涌流识别方法 - Google Patents
一种基于权重数学形态学的变压器励磁涌流识别方法 Download PDFInfo
- Publication number
- CN103259251B CN103259251B CN201310159424.8A CN201310159424A CN103259251B CN 103259251 B CN103259251 B CN 103259251B CN 201310159424 A CN201310159424 A CN 201310159424A CN 103259251 B CN103259251 B CN 103259251B
- Authority
- CN
- China
- Prior art keywords
- transformer
- mathematical morphology
- excitation flow
- sigma
- recognition method
- 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
Links
Landscapes
- Protection Of Transformers (AREA)
Abstract
本发明公开了一种基于权重数学形态学的变压器励磁涌流识别方法,包括以下步骤:(1)按照一定采样频率采集变压器差动保护两侧电流互感器的差动电流信号;(2)判断连续若干次采样值之间的突变量是否大于预设的突变量限定值,如果是则经延时采满半个周期数据后进入步骤(3),否则继续采样;(3)选取数据窗口为半个基频周期长的差动电流信号值,进行权重数学形态学的处理,得到处理后的信号波形;(4)求取波形相关系数,判断波形相关系数是否小于预设的波形相关系数整定值,如果是,则认为变压器发生内部故障;否则判别为变压器励磁涌流。本发明方法具有步骤简单、计算量小、延时小、实现简单、通用性强等优点。
Description
技术领域
本发明涉及变压器继电保护技术领域,特别涉及一种基于权重数学形态学的变压器励磁涌流识别方法。
背景技术
变压器作为电能传输和电压变换的重要电气设备,其运行状态直接影响供电的安全与稳定。同时,变压器造价昂贵,一旦因故障而遭到破坏,检修难度大、时间长,将造成严重的经济损失。因此研究新的性能可靠、技术先进的变压器保护方案具有较大的工程和理论价值。
差动保护因灵敏度高和选择性好而作为变压器的主保护。但是变压器差动保护一直受励磁涌流问题的困扰,在空载合闸或外部故障切除后恢复供电时,差动回路会流入可与内部故障电流相比拟的励磁电流(又称为励磁涌流),从而引起变压器差动保护的误动,因此正确区分励磁涌流和内部故障电流是保证变压器保护可靠动作的关键所在。
由于励磁涌流通常包含有大量以二次谐波为主的高次谐波并且波形之间出现间断,通常采用二次谐波制动原理和间断角原理区分励磁涌流与内部故障电流。二次谐波制动原理是根据检测到变压器励磁涌流的二次谐波含量较高而闭锁保护。但是近年来随着变压器铁芯材料的磁特性变化,励磁涌流中二次谐波含量降低,而在系统带有长线路以及分布电容存在的情况下,内部故障电流的二次谐波含量有时会接近甚至大于励磁涌流的二次谐波含量,从而导致了变压器差动保护的拒动或误动。间断角原理是利用励磁涌流波形之间有较大间断角的特性,通过检测差动电流间断角的大小实现涌流识别,但是由于对硬件的要求较高,并且在CT(current transformer,电流互感器)饱和时因反向电流的出现导致间断角减少,所以一般较少使用。总结以上分析可知传统的识别方法已不能满足现代保护技术高可靠性的要求,因此,进一步探索快速、准确的区分变压器励磁涌流和内部故障电流的新方法以提高变压器差动保护的性能,是十分必要的。
在目前已公开的专利、期刊和会议文献中,许多国内外学者针对如何识别变压器励磁涌流和内部故障电流进行了广泛的研究,并提出了许多新原理、新方法用于励磁涌流识别,也取得了一定成效,主要包括波形对称原理、磁通特性制动原理以及小波理论、模糊数学、数学形态学、神经网络等新手段,但这些方法中仍存在着诸如数据窗口长、保护整定复杂、特殊情况不可靠或计算负担较重等缺陷。
中国华北电力大学的Sun Yang等人在《Novel Theory of Identifying CurrentBased on Half-cycle Sine Waveform》(International Conference on Computer,Mechatronics,Control and Electronic Engineering,2010)(《基于正弦半波的励磁涌流识别新原理》(2010年IEEE国际计算机、机电、控制和电子工程学术会议))中利用输入波形与标准半波正弦比较。该方法需要的时间窗口为半个基频周期,延迟较小,但是该方法精度不高,因为构造的标准半波正弦受最大采样值的影响。当最大采样值与实际正弦波形峰值误差较大时,将不能可靠识别励磁涌流与内部故障电流。
伊朗伊斯兰阿扎德大学的G.Mokryani等人在《Detection of Inrush CurrentBased on Wavelet Transform and LVQ Neural Network》(Transmission andDistribution Conference and Exposition,2010IEEE PES)(《基于小波变换和LVQ神经网络的励磁涌流检测》(2010IEEE电力与能源国际会议))中将小波变换与LVQ神经网络相结合用于鉴别励磁涌流,该方法能可靠识别励磁涌流,但利用神经网络需要大量的训练模型,并且小波变换计算较为复杂和需要的一个周波的数据窗口,这些缺点使得其实现起来比较困难。
英国利物浦大学的Z.Lu,W.H.Tang,T.Y.Ji和Q.H.Wu等人在《AMorphological Scheme for Inrush Identification in Transformer Protection》(IEEETransactions on Power Delivery,Vol.24,No.2,2009)(《一种基于形态学方法的变压器励磁涌流识别》(IEEE输电汇刊,2009年第2期))提出了一种基于数学形态学的分解方法提取励磁涌流的波形特征来进行识别,该算法具有计算负担不大、数据窗口短等特点。但这种方法在一些情况下并不可靠,如当故障出现不是在0度合闸角时出现,在电流的前一个半波将不能可靠识别。
公告号为CN 101567552 A的中国发明专利公开了一种利用形态结构的电力变压器励磁涌流与内部故障识别方法,该发明专利提出了利用正弦半波结构元素分别对变压器三相差流进行形态学分析,求取原始信号与形态学分析结果波形的一致性,其本质为反应原波形信号与正弦波的相似程度,通过比较波形一致性系数的大小实现电力变压器励磁涌流与内部故障电流的识别;该方法虽然克服了一些现有技术的不足,但其采样数据窗口长达25ms,造成了较长时间的延迟,不能满足高压变压器需要快速动作的要求,较长的动作时间也加剧了内部故障电流对变压器的损害。
发明内容
本发明的主要目的在于克服现有技术的缺点与不足,提供一种基于权重数学形态学的变压器励磁涌流识别方法,该方法具有延时小、可靠性高及计算量小的优点。
本发明的目的通过以下的技术方案实现:一种基于权重数学形态学的变压器励磁涌流识别方法,包括以下步骤:
(1)按照一定采样频率采集变压器差动保护两侧电流互感器的差动电流信号;
(2)判断连续若干次采样值之间的突变量是否大于预设的突变量限定值,如果是则经延时采满半个周期数据后进入步骤(3),否则继续采样;
(3)选取数据窗口为半个基频周期长的差动电流信号值,进行权重数学形态学的处理,经过处理后的信号波形表示如下:
其中:
I为差动电流信号,b=[cosmφ,…,cos2φ,cosφ,1,cosφ,cos2φ,…,cosmφ]为结构元素,φ=ωΔt,ω为系统基波角频率,Δt为采样时间间隔,m为结构元素b的长度,s为结构元素b的定义域,k为差动电流信号I的采样点;
(4)求取波形相关系数J,计算方法如下:
Cov(I,D)=E(ID)-E(I)E(D)
(5)判断J是否小于预设的波形相关系数整定值,如果是,则认为变压器发生内部故障;否则判别为变压器励磁涌流。
优选的,所述步骤(2)中,连续3次采样得到的差动电流信号I(n-2N)、I(n-N)、I(n)符合如下启动判据时,即认为发生了故障或空载合闸:
|I(n)-I(n-N)|-|I(n-N)-I(n-2N)|>Iqd;
其中,Iqd为预设的突变量限定值,N为一个周期的采样点数。
更进一步的,所述预设的突变量限定值Iqd设定为0.2In,In为额定电流。
优选的,所述步骤(5)中,预设的波形相关系数整定值为0.9~1.5。
更进一步的,所述预设的波形相关系数整定值优选为1.2。
优选的,结构元素的长度m的值为0.075N,N为一个周期的采样点数。
本发明与现有技术相比,具有如下优点和有益效果:
1、本发明方法仅选取半个基波周期长的数据窗口即可实现励磁涌流的识别,因此本发明方法的延时小。
2、本发明通过计算经权重数学形态学处理后的波形与原始波形的波形相关系数的值,实现了对励磁涌流和内部故障电流的判别,所利用的权重数学形态学算法只涉及到少量的运算,因此本发明方法的步骤简单、计算量小,通过简单的硬件即可实现。
3、本发明对于不同类型的励磁涌流与内部故障电流的识别具有较好的通用性。
附图说明
图1是本发明方法的流程图。
图2(a)和(b)分别是本发明方法中励磁涌流和内部故障电流的波形图。
图3是本发明方法的励磁涌流在数据窗口内的波形及其数学形态学算法处理后的结果。
图4是本发明方法的内部故障电流在数据窗口内的波形及其数学形态学算法处理后的结果。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例1
如图1所示,一种基于权重数学形态学的变压器励磁涌流识别方法,包括以下步骤:
(1)采集变压器差动保护两侧电流互感器的差动电流信号I;
(2)采用fs为4kHZ的采样频率对采集到的差动电流信号I进行采样,得到差动电流信号I在各采样点的值I(k);
(3)当步骤(2)中采样得到的差动电流信号I采样值连续3点的突变量大于预设的突变量限定值Iqd,即符合如下启动判据时即认为发生了故障或空载合闸:
|I(n)-I(n-N)|-|I(n-N)-I(n-2N)|>Iqd;
式中,N为一个周期的采样点数。经10ms延时采满半个周期数据后进行保护程序的计算。本实施例突变量限定值Iqd=0.2A,In为额定电流,为1A,N=80。
(4)选取数据窗口为半个基频周期长10ms的差动电流信号I,对该信号进行权重数学形态学的处理,权重数学形态学的膨胀和腐蚀定义如下:
其中,b=[cosmφ,…,cos2φ,cosφ,1,cosφ,cos2φ,…,cosmφ]为结构元素,φ=ωΔt,ω为系统基波角频率,Δt为采样时间间隔,m为结构元素b的长度,本实施例优选取m=6。定义:
表示经过权重数学形态学处理后的信号波形。在数据窗口内,变压器励磁涌流和内部故障电流经数学形态学处理后得到的波形分别如图3、4所示。
(5)定义波形相关系数J,波形相关系数的实质是反映权重数学形态学处理后的信号波形与原始信号波形的相似程度。对于内部故障电流J≈1,而对于励磁涌流J>>1。波形相关系数J的计算方法如下:
Cov(I,D)=E(ID)-E(I)E(D)
(6)变压器励磁涌流和内部故障的识别判据:根据步骤(1)至步骤(3)计算处理变压器的差流,得到波形相关系数J,对于内部故障电流J≈1,而对于励磁涌流J>>1;由此可得变压器励磁涌流和内部故障电流识别判据如下式所示:
J≥Jset;
式中,Jset为预设的变压器励磁涌流和内部故障电流识别的波形相关系数整定值,本发明Jset整定值优选1.2。当变压器差流满足J<Jset时,判别为变压器发生内部故障;否则判别为变压器励磁涌流。
本发明依据原理:发生内部短路故障时,变压器差流近似于正弦波,如图2(b)所示,而出现励磁涌流时,涌流是非正弦的,如图2(a)所示。
若电流为正弦波,相应的第k个采样值为:
I(k)=Acos(ω·kΔt+θ);
以I(k)为中心点,其左右的采样值可以表示为:
I(k+n)=Acos[ω·(k+n)Δt+θ]
;
I(k-n)=Acos[ω·(k-n)Δt+θ]
则I(k+n)+I(k-n)=2Acos(ω·kΔt+θ)cos(ω·nΔt)=2I(k)cos(ω·nΔt)。
而 所以当电流信号为正弦波时,Dn(k)≈I(k),则J≈1。然而由于变压器励磁涌流波形是非正弦的,故不符合。
上述方法所采用的数学形态学是图象分析中形状表示的重要方法,是形状的定量描述。数学形态学在信号处理中的应用主要是灰度形态学,其两种基本的形态函数是灰度膨胀和灰度腐蚀。假设待处理信号f(n)是采样得到的一维信号,其定义域为D[f]={1,2,3,...,N};N为f(n)的序列长度。g(k)为一维结构元素序列,其定义域为D[g]={1,2,3,...,P};其中,P为g(k)的序列长度,P和N都是整数,N≥P。f(n)在以下公式中记为f,g(k)在以下公式中均记为g。则f(n)的灰度膨胀(f⊕g)(n)和灰度腐蚀(fΘg)(n)分别定义为:
(fΘg)(n)=min{f(n+x)-g(x)|(n+x)∈Df且x∈Dg};
表示膨胀运算,Θ表示腐蚀运算,腐蚀运算是求取极小值,膨胀运算是求取极大值。max表示集合{f(n-x)+g(x)|(n-x)∈Df且x∈Dg}中的最大元素,min表示集合{f(n+x)-g(x)|(n+x)∈Df且x∈Dg}中的最小元素,x为平移变量。
而权重数学形态学是根据标准数学形态学重新定义了膨胀运算与腐蚀运算,分别为:
(fΘg)(n)=min{f(n+x)/g(x)|(n+x)∈Df且x∈Dg};
表示膨胀运算,Θ表示腐蚀运算,腐蚀是求取极小值,膨胀运算是求取极大值。max表示集合中{f(n-x)/g(x)|(n-x)∈Df且x∈Dg}中的最大元素,min表示集合{f(n+x)/g(x)|(n+x)∈Df且x∈Dg}中的最小元素,x为平移变量。
本发明所用的权重数学形态学与标准数学形态学相比,由于其结构元素为正弦结构并且将标准数学形态学的加减运算改为除法运算,因此其能更好地提取出波形信号的正弦特征。
而励磁涌流和内部故障电流的波形有明显的区别,励磁涌流是因变压器铁芯饱和而产生的,其波形是非正弦的,并且每一个波形之间都有一个间断,而内部故障电流基本是正弦波形。因此,利用本发明的权重数学形态学可以提取信号波形的正弦特征,用以区分励磁涌流和内部故障电流。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (6)
1.一种基于权重数学形态学的变压器励磁涌流识别方法,其特征在于,包括以下步骤:
(1)按照一定采样频率采集变压器差动保护两侧电流互感器的差动电流信号;
(2)判断连续若干次采样值之间的突变量是否大于预设的突变量限定值,如果是则经延时采满半个周期数据后进入步骤(3),否则继续采样;
(3)选取数据窗口为半个基频周期长的差动电流信号值,进行权重数学形态学的处理,经过处理后的信号波形表示如下:
其中:
I为差动电流信号,b=[cosmφ,…,cos2φ,cosφ,1,cosφ,cos2φ,…,cosmφ]为结构元素,φ=ωΔt,ω为系统基波角频率,Δt为采样时间间隔,m为结构元素b的长度,s是结构元素b的定义域,k是差动电流信号I的采样点;
(4)求取波形相关系数J,计算方法如下:
Cov(I,D)=E(ID)-E(I)E(D)
(5)判断J是否小于预设的波形相关系数整定值,如果是,则认为变压器发生内部故障;否则判别为变压器励磁涌流。
2.根据权利要求1所述的基于权重数学形态学的变压器励磁涌流识别方法,其特征在于,所述步骤(2)中,连续3次采样得到的差动电流信号I(n-2N)、I(n-N)、I(n)符合如下启动判据时,即认为发生了故障或空载合闸:
|I(n)-I(n-N)|-|I(n-N)-I(n-2N)|>Iqd;
其中,Iqd为预设的突变量限定值,N为一个周期的采样点数。
3.根据权利要求2所述的基于权重数学形态学的变压器励磁涌流识别方法,其特征在于,所述预设的突变量限定值Iqd设定为0.2In,In为额定电流。
4.根据权利要求1所述的基于权重数学形态学的变压器励磁涌流识别方法,其特征在于,所述步骤(5)中,预设的波形相关系数整定值为0.9~1.5。
5.根据权利要求4所述的基于权重数学形态学的变压器励磁涌流识别方法,其特征在于,所述预设的波形相关系数整定值为1.2。
6.根据权利要求2所述的基于权重数学形态学的变压器励磁涌流识别方法,其特征在于,所述步骤(3)中,结构元素的长度m的值为0.075N,N为一个周期的采样点数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310159424.8A CN103259251B (zh) | 2013-05-02 | 2013-05-02 | 一种基于权重数学形态学的变压器励磁涌流识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310159424.8A CN103259251B (zh) | 2013-05-02 | 2013-05-02 | 一种基于权重数学形态学的变压器励磁涌流识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103259251A CN103259251A (zh) | 2013-08-21 |
CN103259251B true CN103259251B (zh) | 2015-10-28 |
Family
ID=48963025
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310159424.8A Active CN103259251B (zh) | 2013-05-02 | 2013-05-02 | 一种基于权重数学形态学的变压器励磁涌流识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103259251B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103545789B (zh) * | 2013-08-26 | 2016-12-28 | 江苏科技大学 | 变压器差动保护的励磁涌流模糊识别方法 |
CN103490394B (zh) * | 2013-09-30 | 2016-07-06 | 山东大学 | 有源配电网的自同步正序故障分量电流差动保护方法 |
CN103683198B (zh) * | 2013-12-03 | 2017-01-04 | 昆明理工大学 | 一种基于差动电流相邻阶次差分构成的平面上相邻点距离的励磁涌流快速识别方法 |
CN104200055B (zh) * | 2014-09-26 | 2017-06-20 | 国家电网公司 | 特高压变压器空载合闸情况下的励磁涌流仿真方法及装置 |
CN104993455A (zh) * | 2015-07-28 | 2015-10-21 | 株洲南车时代电气股份有限公司 | 一种牵引变压器过流保护方法 |
CN106058810B (zh) * | 2016-07-08 | 2018-12-25 | 山东鲁能智能技术有限公司 | 一种基于功率故障分量判据的励磁涌流识别方法 |
CN106324328B (zh) * | 2016-08-09 | 2019-01-29 | 华南理工大学 | 一种基于形态学级联腐蚀运算的变压器励磁涌流识别方法 |
CN109638780B (zh) * | 2019-01-08 | 2020-04-07 | 广东电网有限责任公司 | 一种保护装置控制方法 |
CN110445095B (zh) * | 2019-07-29 | 2021-06-04 | 天津大学 | 基于复合环流与复合零序电流波形相关度的励磁涌流识别方法 |
CN111273108B (zh) * | 2020-03-17 | 2022-06-21 | 深圳供电局有限公司 | 变压器空充跳闸原因的判断方法 |
CN112039021B (zh) * | 2020-09-08 | 2022-04-12 | 河南理工大学 | 基于差分波形参量的变压器励磁涌流鉴别方法 |
CN114167117A (zh) * | 2021-12-02 | 2022-03-11 | 合肥工业大学 | 一种双绕组变压器差动保护励磁涌流的识别方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101257208A (zh) * | 2007-12-21 | 2008-09-03 | 清华大学 | 一种变压器励磁涌流的鉴别方法 |
CN101567552A (zh) * | 2009-06-03 | 2009-10-28 | 昆明理工大学 | 一种利用形态结构的电力变压器励磁涌流与内部短路识别方法 |
CN103050941A (zh) * | 2012-12-19 | 2013-04-17 | 华南理工大学 | 一种基于形态学梯度的变压器励磁涌流识别方法 |
-
2013
- 2013-05-02 CN CN201310159424.8A patent/CN103259251B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101257208A (zh) * | 2007-12-21 | 2008-09-03 | 清华大学 | 一种变压器励磁涌流的鉴别方法 |
CN101567552A (zh) * | 2009-06-03 | 2009-10-28 | 昆明理工大学 | 一种利用形态结构的电力变压器励磁涌流与内部短路识别方法 |
CN103050941A (zh) * | 2012-12-19 | 2013-04-17 | 华南理工大学 | 一种基于形态学梯度的变压器励磁涌流识别方法 |
Non-Patent Citations (2)
Title |
---|
利用数学形态学提取暂态量的变压器保护新原理;马静 等;《中国电机工程学报》;20060330;第26卷(第6期);19-23 * |
采用改进数学形态学识别变压器励磁涌流的新方法;黄家栋 罗伟强;《中国电机工程学报》;20090305;第29卷(第7期);98-105 * |
Also Published As
Publication number | Publication date |
---|---|
CN103259251A (zh) | 2013-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103259251B (zh) | 一种基于权重数学形态学的变压器励磁涌流识别方法 | |
CN103050941B (zh) | 一种基于形态学梯度的变压器励磁涌流识别方法 | |
Dash et al. | Hybrid S-transform and Kalman filtering approach for detection and measurement of short duration disturbances in power networks | |
Liu et al. | Application of the variational mode decomposition-based time and time–frequency domain analysis on series DC arc fault detection of photovoltaic arrays | |
CN103399257B (zh) | 一种中性点不接地系统的铁磁谐振故障检测方法 | |
CN105262051B (zh) | 基于采样序列绝对值偏态分布的变压器励磁涌流鉴别方法 | |
CN103323698A (zh) | 一种新型电流互感器饱和识别方法 | |
Zuhaib et al. | A novel method for locating the source of sustained oscillation in power system using synchrophasors data | |
Li et al. | A fault pattern and convolutional neural network based single-phase earth fault identification method for distribution network | |
Chuanli et al. | A new approach to detect transformer inrush current by applying wavelet transform | |
CN105067906A (zh) | 一种基于三阶中心矩的电流互感器饱和检测方法 | |
CN102435896A (zh) | 船舶中压电力系统中间歇性接地故障的快速识别方法 | |
CN101820166B (zh) | 基于变压器等效励磁阻抗频域特征的励磁涌流判别方法 | |
CN105606892A (zh) | 一种基于sst变换的电网谐波与间谐波分析方法 | |
CN105116251A (zh) | 变压器故障鉴别保护方法 | |
CN104410044A (zh) | 基于峭度和偏度的变压器励磁涌流识别方法 | |
CN105044458B (zh) | 一种基于tt变换的电力系统谐波检测方法 | |
CN109088402A (zh) | 输电线路时域全波形保护装置及其相突变量方向判断方法 | |
CN103560509B (zh) | 一种基于小波分析的电压下陷检测装置及该装置的控制方法 | |
Yongli et al. | Fault location method based on the periodicity of the transient voltage traveling wave | |
CN108092234B (zh) | 适用于变频电动机的差动保护方法和装置 | |
Samantaray et al. | Spectral energy function for fault detection during power swing | |
CN113866628B (zh) | 一种基于电流、电压双相序比对的开关状态分析方法 | |
Zhang et al. | Fast Fault Detection and Location System for Distribution Network Lines Based on Power Electronic Disturbance Signals | |
CN203368041U (zh) | 一种在线检测自励式mcr中励磁控制晶闸管状态的装置 |
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 |