CN103513212B - 基于重构的bh曲线特征ct状态识别以及不饱和度计算方法 - Google Patents

基于重构的bh曲线特征ct状态识别以及不饱和度计算方法 Download PDF

Info

Publication number
CN103513212B
CN103513212B CN201310414726.5A CN201310414726A CN103513212B CN 103513212 B CN103513212 B CN 103513212B CN 201310414726 A CN201310414726 A CN 201310414726A CN 103513212 B CN103513212 B CN 103513212B
Authority
CN
China
Prior art keywords
degree
unsaturation
iron core
reconstruct
current
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
CN201310414726.5A
Other languages
English (en)
Other versions
CN103513212A (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201310414726.5A priority Critical patent/CN103513212B/zh
Publication of CN103513212A publication Critical patent/CN103513212A/zh
Application granted granted Critical
Publication of CN103513212B publication Critical patent/CN103513212B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明提供了一种基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，用于电力系统差动保护两侧CT工作状态的识别，具体是通过对差动两侧电流附带的有效信息关键特征的提取与识别，在不增加任何硬件情况下，利用差动两侧电流的采样数据，计算得出磁场强度H以及磁感应强度B的变化趋势，构造出反映铁芯内部磁通变化过程的BH曲线，并提供量化的CT铁芯状态信息，根据不同状态下的BH曲线特征，自动调整识别门槛，有效识别外部穿越性电流导致的CT饱和，并且能够有效区分CT开路、区内故障及区外故障导致的饱和，降低由于CT饱和程度差异导致的差动误动作概率，提高差动保护动作的可靠性。

Description

基于重构的BH曲线特征CT状态识别以及不饱和度计算方法
技术领域
本发明涉及电力系统差动保护两侧CT工作状态的识别方法，具体是通过两侧电流附带的有效信息关键特征的提取与识别，提供量化后的CT铁芯状态信息，为差动保护提供判据，属于电力系统自动化控制技术领域。
背景技术
电力系统中差动保护通过对被保护区域两侧电流的连续交流采样，实时计算差动电流以及制动电流，有效识别区内及区外故障，保护设备的可靠运行。对于差动保护CT而言，由于剩磁、二次回路阻抗、一次电流中的非周期分量等因素导致CT铁芯进入饱和，饱和后的CT二次电流出现畸变失真，实时计算得出的差动电流变大并容易导致差动保护误动作，严重影响运行的可靠性。
在电力系统中，由于电动机、发电机、电抗器、线路保护两侧CT特性差异导致的差动误动作频繁发生（《电力系统继电保护动作实例分析》景敏慧2012），影响保护的可靠性，在工程设计阶段，虽已根据可能出现的最大故障电流进行了保护用CT的计算选型，但电力系统中CT稳态饱和较少出现，而对于突发故障或瞬间合闸导致的非周期分量，特别是两侧差动CT二次回路的阻抗大小存在明显差异时更易导致两侧CT的饱和程度出现差异，不能准确传变一次电流，从而产生了传变误差，导致差动保护的错误判断，从而影响保护动作的可靠性。
由于电力系统呈感性，合闸瞬间的励磁涌流以及故障电流均存在大量非周期分量，这些电流对于目前大量应用的P类保护级CT而言，是个严峻的考验，往往会导致CT暂态饱和，因此需要采取有效识别此类饱和，避免保护的误动作，确保电力系统的稳定可靠运行。
国内外针对电流互感器饱和的问题提出了多种识别方法，有时差法、异步法、差分法、导数法、谐波比法、小波法、形态学法等，其中以南瑞继保提出的异步法和国电南自提出的时差法应用最为广泛，并且已应用到实际工程中。上述各种识别方法都是在CT饱和出现后根据饱和电流的某个或多个特征提出，识别原理相对复杂，并且无法反映CT铁芯的变化过程，只能定性判断，无法提供量化的铁芯饱和程度信息。
发明内容
使用两侧差动CT电流的瞬时采样数据，计算得出铁芯中磁场强度H以及磁感应强度B的变化趋势，重构出CT铁芯的BH曲线，反映铁芯磁通变化过程。
根据CT铁芯不同特性，将CT伏安特性曲线中拐点电流放大k_sat倍后作为不饱和区域阈值,设置合适的不饱和区域。
通过实时计算T_samp内滞留在不饱和区的采样点数与T_samp内采样总点数的比值，得出铁芯不饱和度S_nsat，将原本无法量化的铁芯磁滞过程数字化。
根据各种运行状态下重构的BH曲线特征，设计了一个识别不同状态的不饱和度阈值曲线。
根据本时刻之前T_dis内电流大小自动调整不饱和度阈值，适应各种工作状态。
根据权利2所述的基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，其特征在于：所述k_sat取1.7。
根据权利3所述的基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，其特征在于：所述T_samp取20ms。
根据权利3所述的基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，其特征在于：识别为铁芯饱和后，通过不饱和度可以计算得出饱和度。
根据权利4所述的基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，其特征在于：包括但不限于本曲线，基于此原理衍生的其他曲线方式也在保护权利之中，并且参数的调整不影响权利的主张。
根据权利5所述的基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，其特征在于：T_dis可取20ms或10ms，20ms可以保证在励磁涌流等电流偏向一侧时的灵敏性，10ms则可以缩短识别时间。
有益效果
在不增加任何硬件情况下，本发明利用差动两侧电流的采样数据，可以有效识别外部穿越性电流导致的CT饱和，并且能够有效区分CT开路、区内故障及区内故障导致的饱和，降低了由于CT饱和因素导致的差动误动作概率，提高了差动保护动作的可靠性。
附图说明
图1为在正常运行、一侧CT开路、区内故障、区外故障不同工作状态下，使用两侧电流重构出来的BH曲线；
图2为区外故障两侧CT饱和程度不一致时，使用两侧电流重构出来的BH曲线；
图3为不饱和度阈值曲线；
图4为“正常运行→一侧CT开路→正常运行→区内故障（未饱和）→区内故障（饱和）→区内故障消除恢复正常运行→区外故障导致轻微饱和→区外故障导致严重饱和”不同状态下的实时识别效果；
图5为区外故障两侧CT饱和程度不一致的识别效果。
图中曲线说明：
A.正常运行无故障CT未饱和；
B.区外故障导致CT饱和；
C.区内故障CT未饱和；
D.区内故障CT饱和；
E.一侧CT开路；
F.CT铁芯拐点对应的饱和识别阈值。
具体实施方式
在差动保护CT回路中，当差动保护区域外部故障时，巨大的电流穿过差动保护区域形成穿越性电流，而由于两侧CT的特性存在差异，往往饱和程度存在差异，形成差流导致差动保护装置错误判断，在非本保护区域内故障时错误的跳闸降低了差动保护的可靠性。
根据CT的工作原理，一次电流产生交变磁通，经过铁芯耦合到二次线圈形成感应电势并产生二次电流去平衡一次电流的磁通变化影响，但由于漏磁、铁芯铁损、外部回路阻抗电压降等因素，一次电流与二次电流形成的磁通并不能完全抵消，对CT二次回路而言，微机保护装置消耗负载很小，在电流互感器的二次负载中电缆电阻占有绝大部分，故二次阻抗呈电阻性（《电力系统继电保护动作实例分析》景敏慧2012中国电力出版社），而电阻在CT饱和暂态过程中不会像电感那样因频率变化而阻抗发生变化，是个固定的数值。
根据大量的工程实际动作录波以及软件仿真结果，发现当两侧CT二次回路阻抗存在差异时，往往二次回路阻抗较大一侧先出现明显的饱和过程，而较小的另一侧则仍能够正确进行电流的变送，鉴于此，将准确变送的一侧视为一次电流，两侧CT二次回路电流瞬时采样值的差值为等效励磁电流，对畸变大的一侧二次电流进行的实时积分得出铁芯中对应的磁感应强度变化趋势。由于CT漏磁以及铁损所占比重较小，铁芯内部的剩余磁通Ф主要与CT两侧的二次回路阻抗电压U密切相关，CT两端的二次回路的电压降，从而得出铁芯内部的磁通，铁芯内部的磁感应强度，而外部磁场强度可以通过两侧电流计算得出：，其中为计算中使用的系数，s为铁芯的截面积，为了简化分析，在不影响分析结果的前提下，本次上述三个参数取1，为差动保护两侧CT的电流采样数据。
使用上述计算得出的磁场强度H以及磁感应强度B进行饱和一侧CT铁芯的BH曲线的重构，图1为各种典型运行状态下的重构BH曲线，图2为区外故障时电流导致的两侧CT饱和程度不一致时重构的BH曲线。通过图1、图2可以看出，未饱和时励磁电流较小，而饱和后的励磁电流变化明显，并且不同状态下重构的BH曲线具有明显的特征，CT铁芯饱和后按照磁滞回线特性往返，需要经过非饱和区并滞留一段时间方能到达反向饱和区，而CT开路以及区内短路则直接穿过不饱和区域抵达对侧，滞留时间短。
根据不同状态下在不饱和阈值区内滞留时间差异，通过对滞留在图1、图2中不饱和识别阈值F区域内的采样点数的统计，利用电流采样点数存在差异的特征，设置合理的比例阈值来识别CT是否进入饱和状态。参考磁滞回线曲线簇特性，不同大小的励磁电流对应的BH曲线在不饱和阈值区域内滞留的采样点比例不同，故此阈值不能为恒定值，需要根据励磁电流的峰值大小而动态变化。
由于不同CT铁芯特性不同，为了选择合适的不饱和区域F所对应的识别阈值，可以根据CT伏安特性测试报告数据中伏安特性曲线的拐点电流来设置合适的不饱和阈值，由于工频下伏安特性曲线测试已考虑CT铁芯的铁损因素，并且该数值为工频下基波的模值，而不饱和阈值需要根据电流的瞬时值进行识别，故需要将模值转换为对应的峰值，并且实际应用中为了可靠识别起见，还需将其放大1.2倍，k_sat=1.414*1.2=1.6968，本次取1.7。
为了能够判断铁芯的工作状态，通过实时计算当前时刻之前T_samp内滞留在不饱和区的采样点数与T_samp内采样总点数的比值，可以计算得出铁芯当前时刻的实时不饱和度S_nsat，将原本无法量化的非线性的铁芯磁滞过程数字化，如识别为饱和状态，则在此基础上亦可得出实时饱和度（1-S_nsat），本次T_samp取50HZ时一个周期时间20ms。
为了能够通过不饱和度有效识别不同状态下重构出来的BH曲线特征，设计了如图3所示的不饱和度阈值曲线，该曲线基于大量的实际故障录波数据分析以及仿真验证后归纳而出，实时计算得出的不饱和度如进入曲线的上部饱和区并且延时满足要求，则可以判断为区外故障导致的CT饱和，从而进行差动保护的闭锁，避免差动保护的误动作。
由于电流变化的不可预期性，需要不断的自动调整不饱和度识别阈值，根据电流的连续不可突变性，采用本时刻之前T_dis时间窗口内的励磁电流峰值的绝对值来自动调整识别阈值，故本算法具有一定的滞后性，但识别时间不大于20ms，考虑到差动保护的判断往往也需要1～2个周波时间，可以将此判断结果作为差动输出的一个闭锁条件，避免误动。T_dis可取20ms或10ms，20ms可以保证在励磁涌流等电流偏向一侧时的灵敏性，10ms则可以缩短识别时间。
该算法较好的解决了目前差动CT饱和识别问题，并且通过采用励磁电流与不饱和区阈值的比较，判断是否进入饱和区，将曲线转换为简单便于程序实现的时间数字序列，同时给出了能够量化反映铁芯状态的不饱和度的计算方法，便于微机保护装置的实现。
图4为对应于与图1中的“正常运行→一侧CT开路→正常运行→区内故障（未饱和）→区内故障（饱和）→区内故障消除恢复正常运行→区外故障导致轻微饱和→区外故障导致严重饱和”不同状态下，T_dis取20ms时本方法的实际识别结果，从图中可以看出，区外故障导致CT饱和的状态B被正确识别出来。
图5为对应于与图2中的区外故障时两侧CT饱和程度不一致时，T_dis取20ms时本方法的实际识别结果，从图中可以看出，识别正确。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何不经过创造性劳动想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书所限定的保护范围为准。

Claims (5)

1.一种基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，其特征在于：使用两侧差动CT电流和CT两侧的二次回路阻抗电压的瞬时采样数据，计算得出铁芯中磁场强度H以及磁感应强度B的变化趋势，重构出CT铁芯的BH曲线，反映铁芯磁通变化过程；根据CT铁芯不同特性，将CT伏安特性曲线中拐点电流放大1.7倍后作为不饱和区域阈值，设置合适的不饱和区域；通过实时计算T_samp内滞留在不饱和区的采样点数与T_samp内采样总点数的比值，得出铁芯不饱和度S_nsat，将原本无法量化的铁芯磁滞过程数字化，其中，T_samp表示一个周期的采样时间；根据各种运行状态下重构的BH曲线特征，设计了一个识别不同状态的不饱和度阈值曲线；根据当前时刻之前T_dis内电流大小自动调整不饱和度阈值，适应各种工作状态，其中，T_dis表示一个周期的识别窗口时间。
2.根据权利要求1所述的基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，其特征在于：所述T_samp取20ms。
3.根据权利要求1所述的基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，其特征在于：识别为铁芯饱和后，通过不饱和度可以计算得出饱和度。
4.根据权利要求1所述的基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，其特征在于：参数是可以调整的。
5.根据权利要求1所述的基于重构的BH曲线特征CT状态识别以及不饱和度计算方法，其特征在于：T_dis可取20ms或10ms；20ms可以保证在励磁涌流等电流偏向一侧时的灵敏性，10ms则可以缩短识别时间。
CN201310414726.5A 2013-09-12 2013-09-12 基于重构的bh曲线特征ct状态识别以及不饱和度计算方法 Active CN103513212B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310414726.5A CN103513212B (zh) 2013-09-12 2013-09-12 基于重构的bh曲线特征ct状态识别以及不饱和度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310414726.5A CN103513212B (zh) 2013-09-12 2013-09-12 基于重构的bh曲线特征ct状态识别以及不饱和度计算方法

Publications (2)

Publication Number Publication Date
CN103513212A CN103513212A (zh) 2014-01-15
CN103513212B true CN103513212B (zh) 2016-01-27

Family

ID=49896235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310414726.5A Active CN103513212B (zh) 2013-09-12 2013-09-12 基于重构的bh曲线特征ct状态识别以及不饱和度计算方法

Country Status (1)

Country Link
CN (1) CN103513212B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105186931B (zh) * 2015-09-09 2018-01-23 乐晓蓉 异步电机启动过程故障电流的识别方法
CN105137166B (zh) * 2015-09-09 2017-11-28 乐晓蓉 电力系统励磁涌流识别方法
CN105223450B (zh) * 2015-11-02 2018-06-26 乐晓蓉 变压器铁芯状态在线监测方法
CN105182261B (zh) * 2015-11-02 2018-06-22 乐晓蓉 线圈内部铁芯磁场强度非接触式测量方法
CN106199480B (zh) * 2016-07-16 2018-11-13 太原理工大学 基于b样条小波变换的电流互感器饱和检测方法
CN106291057B (zh) * 2016-08-14 2018-11-16 卢庆港 基于多ct下重构bh曲线特征的暂态饱和电流识别方法
CN106932747B (zh) * 2017-05-08 2019-08-13 中国南方电网有限责任公司 一种全局和局部特征相结合的电流互感器饱和检测方法
CN114035138A (zh) * 2021-11-08 2022-02-11 许昌许继软件技术有限公司 一种适用于输电线路的电流互感器饱和判别方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5993568A (en) * 1998-03-25 1999-11-30 Nkk Corporation Soft magnetic alloy sheet having low residual magnetic flux density
CN103105594A (zh) * 2013-01-24 2013-05-15 江苏省电力公司电力科学研究院 基于小信号斜率映射的电流互感器剩磁检测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0713625B2 (ja) * 1987-12-04 1995-02-15 株式会社日立製作所 劣化検出方法及びその装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5993568A (en) * 1998-03-25 1999-11-30 Nkk Corporation Soft magnetic alloy sheet having low residual magnetic flux density
CN103105594A (zh) * 2013-01-24 2013-05-15 江苏省电力公司电力科学研究院 基于小信号斜率映射的电流互感器剩磁检测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
P类电流互感器饱和原因分析及对策;景敏慧等;《电力系统自动化》;20071110;第31卷(第21期);第95-96页第2.1节，图5 *
切除外部故障时电流互感器局部暂态饱和对变压器差动保护的影响及对策;袁宇波等;《中国电机工程学报》;20050531;第25卷(第10期);第12-17页 *
基于堆栈操作的电流互感器铁芯磁化模型;崔迎宾等;《电力系统自动化》;20110510;第35卷(第9期);第86-91页 *
电流互感器稳态及暂态饱和特性的分析研究;黄鸣宇;《宁夏电力》;20091031(第5期);第1-4页 *

Also Published As

Publication number Publication date
CN103513212A (zh) 2014-01-15

Similar Documents

Publication Publication Date Title
CN103513212B (zh) 基于重构的bh曲线特征ct状态识别以及不饱和度计算方法
Kang et al. A CT saturation detection algorithm
Hajipour et al. Current-transformer saturation compensation for transformer differential relays
Bo et al. A new technique for transformer protection based on transient detection
Lu et al. A morphological scheme for inrush identification in transformer protection
CN103323698B (zh) 一种新型电流互感器饱和识别方法
Sahebi et al. Efficient method for discrimination between inrush current and internal faults in power transformers based on the non‐saturation zone
CN105262051B (zh) 基于采样序列绝对值偏态分布的变压器励磁涌流鉴别方法
CN104237731A (zh) 基于eemd与能量法的谐振接地配电网单相接地故障选线方法
CN105653882A (zh) 一种利用电流波形特征识别励磁涌流的方法
CN105140893A (zh) 一种差动保护ct饱和识别方法
CN105353271A (zh) 一种电流差动保护饱和判别方法
Altuve et al. Sizing current transformers for line protection applications
Sahebi et al. Identifying internal fault from magnetizing conditions in power transformer using the cascaded implementation of wavelet transform and empirical mode decomposition
US11435411B2 (en) Time domain distance protection based on polarity comparison for power transmission lines
Saravanan et al. Inrush blocking scheme in transformer differential protection
CN109740105B (zh) 一种用于母线保护装置的抗工频磁场干扰方法
Rahmati et al. A fast WT‐based algorithm to distinguish between transformer internal faults and inrush currents
CN102082422A (zh) 防止非故障相二次电流突变导致线路纵差保护误动的方法
CN103983890A (zh) 交流串扰对继电器影响产生的故障点查找方法和测试电路
Kang et al. Compensating algorithm suitable for use with measurement-type current transformers for protection
De Alencar et al. Transformer inrush currents and internal faults identification in power transformers using wavelet energy gradient
Moosavi et al. A new fast method for improvement of power transformer differential protection based on discrete energy separation algorithm
CN107748836B (zh) 一种故障时电流互感器铁芯饱和时间计算方法
Abniki et al. A novel technique for internal fault detection of power transformers based on moving windows

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