CN102842890A - Realizing method of current change rate protection - Google Patents

Realizing method of current change rate protection Download PDF

Info

Publication number
CN102842890A
CN102842890A CN201210357913XA CN201210357913A CN102842890A CN 102842890 A CN102842890 A CN 102842890A CN 201210357913X A CN201210357913X A CN 201210357913XA CN 201210357913 A CN201210357913 A CN 201210357913A CN 102842890 A CN102842890 A CN 102842890A
Authority
CN
China
Prior art keywords
current
rate
value
rising
protection
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.)
Granted
Application number
CN201210357913XA
Other languages
Chinese (zh)
Other versions
CN102842890B (en
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.)
Nari Technology Co Ltd
Original Assignee
Nari Technology 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 Nari Technology Co Ltd filed Critical Nari Technology Co Ltd
Priority to CN201210357913.XA priority Critical patent/CN102842890B/en
Publication of CN102842890A publication Critical patent/CN102842890A/en
Application granted granted Critical
Publication of CN102842890B publication Critical patent/CN102842890B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Emergency Protection Circuit Devices (AREA)
  • Electric Propulsion And Braking For Vehicles (AREA)

Abstract

The invention discloses a realizing method of current change rate protection, wherein a long data window mean algorithm is combined with a short data window mean algorithm, and the comprehensive logic judgment of mutually matching the two data window mean algorithms with different lengths well solves the contradiction of quickness and reliability of the current change rate protection to improve the integral performance of the rail direct current protection.

Description

A kind of implementation method of current changing rate protection
Technical field
The implementation method of current changing rate protection belongs to electrical technology field.
Background technology
Along with developing rapidly of urban rail transit in China in recent years, DC traction power-supply system has obtained application more and more widely.In DC traction power-supply system; Because the direct-current short circuit electric current is different from alternating current, dc switch disjunction prospective short circuit current stationary value is difficulty very, therefore; Require dc switch before fault current reaches the expection stationary value, just not answer the disjunction fault current, and more little good more.Because the current changing rate protection is distinguish short-circuit fault current and locomotive starting current effectively, can protect and be used as in fault early period of origination fast detecting to fault in the first-selection of subway direct current protecting.The sensitivity and the reliability that how to guarantee the current changing rate protection are one of key points of protection.
Summary of the invention
The object of the present invention is to provide a kind of current changing rate protection implementation method, solve the rapidity of current changing rate protection and the contradiction between the reliability preferably, improved the overall performance of track direct current protecting.
The current changing rate protection is divided into current increment protection and current-rising-rate protection.
The current increment protection if monitoring current increment Delta I is higher than corresponding setting value, and continues the corresponding time, then protection action outlet.In-plant nonmetal character short trouble (metallicity direct short-circuit fault is tripped by the electromagnetic tripping apparatus of circuit breaker self) during current increment Δ I protection is primarily aimed at.
The current-rising-rate protection starts Δ T (for setting value) afterwards if monitor the current changing rate protection, and current increment value Δ I then protects the action outlet greater than minimum threshold definite value.Be primarily aimed at remote or resistive type fault.
The general average value filtering algorithm that adopts in direct current, i.e. the R point data window algorithm of averaging.If it is bigger that the R value is chosen, filter effect is good, and computational accuracy is higher, can cause the current changing rate protection to delay start-up time, thereby has influence on choosing of current reference value in the current changing rate protection, is difficult to satisfy the sensitivity demand of protection.Less but if the R value is chosen, protection has higher rapidity and sensitivity, but has the malfunction possibility.
The present invention provides a kind of implementation method of current changing rate protection, it is characterized in that, comprises following steps:
1) data acquisition: the acquired original by FPGA completion data, after the first Filtering Processing of the intermediate value method of average, data are stored among the data buffer area i [N+1], N is that the data sampling of 20ms is counted;
2) calculate current-rising-rate K value:
Current-rising-rate K computing formula K=di/dt;
K 1[k]=(i[k]-i[k-z 1])/t 1
K 2[k]=(i[k]-i[k-z 2])/t 2
Wherein:
K representes the current sampling point counting, between 0 to N+1, circulates;
The current sampling data that i [k] expression k is ordered;
I [k-z 1] expression k-z 1The current sampling data of point;
I [k-z 2] expression k-z 2The current sampling data of point;
t 1Expression k point and k-z 1The time interval between the point;
t 2Expression k point and k-z 2The time interval between the point;
K 1[k], K 2[k] is illustrated respectively in t 1The time interval, t 2First current-rising-rate under the time interval, the instantaneous value of second current-rising-rate.
3) data filtering:
31) data window R: be used for the number of data points of calculating mean value, R can not be greater than N+1;
32) data window R chooses smaller value R 1, calculate current average Calculate the mean value of first current-rising-rate K 1 R = 1 R 1 Σ m = 0 R 1 K ( k - m ) ;
33) data window R chooses higher value R 2, calculate the mean value of second current-rising-rate K 2 R = 1 R 2 Σ m = 0 R 2 K ( k - m ) .
4) current changing rate protection starting component: as the mean value K of first current-rising-rate 1RMean value K with second current-rising-rate 2RAll, then be judged to current changing rate protection starting component and start greater than starting threshold value; If current changing rate protection starting component starts, as the mean value K of first current-rising-rate 1RMean value K with second current-rising-rate 2RAll less than returning threshold value, current changing rate protection starting component just returns.
5) choosing of current reference value: do not start at current changing rate protection starting component, and the mean value K of first current-rising-rate 1RMean value K with second current-rising-rate 2RDuring all less than the startup threshold value, the current reference value of real-time update current changing rate protection, current reference value power taking levelling average I FRInstantaneous value; Mean value K when first current-rising-rate 1RMean value K with second current-rising-rate 2RAmong both any one stops to upgrade when starting threshold value.
6) current increment protection: when current changing rate protection starting component starts, get into increment protection logic determines, increment protection logic determines adopts current average I FRInstantaneous value subtracts the current reference value as the protection input value; If the protection input value is protected definite value greater than increment, and satisfied time-delay demand, then increment protection action; If the increment protection is moved, then no longer judge starting component, need only judge that whether input value protects definite value less than increment, if satisfy condition, then returns.
7) current-rising-rate protection: when current changing rate protection starting component starts; Get into current-rising-rate protection logic determines; Current-rising-rate picks up counting, if current changing rate protection starting component starts and last till the time-delay of adjusting of current-rising-rate protection always, then adopts current average I FRInstantaneous value subtracts the current reference value as input value, if the smallest incremental definite value that input value is protected greater than current-rising-rate, then current-rising-rate protection action.
The software that operates in the CPU that DSP is a core CPU has been realized the implementation method that this current changing rate is protected.Experimental test shows that the current changing rate protection is after this method is handled, and definite value precision, operate time and stability satisfy the requirement of track direct current protecting system fully.
The beneficial effect that the present invention reached: this method combines the advantage of long data window mean value algorithm and short data window mean value algorithm; The integrated logic that adopts two kinds of different length data window mean value algorithms to cooperatively interact is judged; Solve the rapidity of current changing rate protection and the contradiction between the reliability preferably, improved the overall performance of track direct current protecting.
Embodiment
Further describe in the face of the present invention down.Following examples only are used for technical scheme of the present invention more clearly is described, and can not limit protection scope of the present invention with this.
Method of measurement proposed by the invention mainly embodies on the software inner with operating in DSP, and concrete steps are following:
1) data acquisition: the acquired original by FPGA accomplishes data after the first Filtering Processing of the intermediate value method of average, is stored in data among the data buffer area i [N+1]; Accomplished the acquired original of data earlier by FPGA, frequency acquisition F can reach 40K, and median filtering method is got 4 point data, removes maximum and minimum value, and remaining two values are asked on average, and this mean value is stored among the data buffer area i [N+1], and wherein N gets 200.
2) calculate current-rising-rate K value:
Current-rising-rate K computing formula K=di/dt;
K 1[k]=(i[k]-i[k-z 1])/t 1
K 2[k]=(i[k]-i[k-z 2])/t 2
Wherein: z 1=10, z 2=30;
K representes the current sampling point counting, circulation between 0 to 201;
The current sampling data that i [k] expression k is ordered;
I [k-z 1] expression k-z 1The current sampling data of point;
I [k-z 2] expression k-z 2The current sampling data of point;
t 1Expression k point and k-z 1The time interval between the point;
t 2Expression k point and k-z 2The time interval between the point;
K 1[k], K 2[k] is illustrated respectively in t 1The time interval, t 2First current-rising-rate under the time interval, the instantaneous value of second current-rising-rate.
3) data filtering:
31) data window R chooses smaller value R 1=10, calculate current average I FR, calculate the mean value K of first current-rising-rate 1R
32) data window R chooses higher value R 2=30, calculate the mean value K of second current-rising-rate 2R
4) current changing rate protection starting component: as the mean value K of first current-rising-rate 1RMean value K with second current-rising-rate 2RAll, then be judged to current changing rate protection starting component and start greater than starting threshold value; If current changing rate protection starting component starts, as the mean value K of first current-rising-rate 1RMean value K with second current-rising-rate 2RAll less than returning threshold value, current changing rate protection starting component just returns;
5) choosing of current reference value: do not start at current changing rate protection starting component, and the mean value K of first current-rising-rate 1RMean value K with second current-rising-rate 2RDuring all less than the startup threshold value, the current reference value of real-time update current changing rate protection, current reference value power taking levelling average I FRInstantaneous value; Mean value K when first current-rising-rate 1RMean value K with second current-rising-rate 2RAmong both any one stops to upgrade when starting threshold value;
6) current increment protection: when current changing rate protection starting component starts, get into increment protection logic determines, increment protection logic determines adopts current average I FRInstantaneous value subtracts the current reference value as the protection input value; If the protection input value is protected definite value greater than increment, and satisfied time-delay demand, then increment protection action; If the increment protection is moved, judge no longer then whether current changing rate protection starting component starts, need only judge that whether input value protects definite value less than increment, if satisfy condition, then returns;
7) current-rising-rate protection: when current changing rate protection starting component starts; Get into current-rising-rate protection logic determines; Current-rising-rate picks up counting, if current changing rate protection starting component starts and last till the time-delay of adjusting of current-rising-rate protection always, then adopts current average I FRInstantaneous value subtracts the current reference value as input value, if the smallest incremental definite value that input value is protected greater than current-rising-rate, then current-rising-rate protection action.
Above embodiment is merely explanation technological thought of the present invention, can not limit protection scope of the present invention with this, every technological thought that proposes according to the present invention, and any change of on the technical scheme basis, being done all falls within the protection range of the present invention.

Claims (1)

1. the implementation method of a current changing rate protection is characterized in that, may further comprise the steps:
1) data acquisition: the acquired original by FPGA completion data, after the first Filtering Processing of the intermediate value method of average, data are stored among the data buffer area i [N+1], N is that the data sampling of 20ms is counted;
2) calculate current-rising-rate K value:
Current-rising-rate K computing formula K=di/dt;
K 1[k]=(i[k]-i[k-z 1])/t 1
K 2[k]=(i[k]-i[k-z 2])/t 2
Wherein:
K representes the current sampling point counting, between 0 to N+1, circulates;
The current sampling data that i [k] expression k is ordered;
I [k-z 1] expression k-z 1The current sampling data of point;
I [k-z 2] expression k-z 2The current sampling data of point;
t 1Expression k point and k-z 1The time interval between the point;
t 2Expression k point and k-z 2The time interval between the point;
K 1[k], K 2[k] is illustrated respectively in t 1The time interval, t 2First current-rising-rate under the time interval, the instantaneous value of second current-rising-rate.
3) data filtering:
31) data window R: be used for the number of data points of calculating mean value, R can not be greater than N+1;
32) data window R chooses smaller value R 1, calculate current average
Figure FDA00002185075600011
, calculate the mean value of first current-rising-rate K 1 R = 1 R 1 Σ m = 0 R 1 K ( k - m )
33) data window R chooses higher value R 2, calculate the mean value of second current-rising-rate K 2 R = 1 R 2 Σ m = 0 R 2 K ( k - m ) .
4) current changing rate protection starting component: as the mean value K of first current-rising-rate 1RMean value K with second current-rising-rate 2RAll, then be judged to current changing rate protection starting component and start greater than starting threshold value; If current changing rate protection starting component starts, as the mean value K of first current-rising-rate 1RMean value K with second current-rising-rate 2RAll less than returning threshold value, current changing rate protection starting component just returns.
5) choosing of current reference value: do not start at current changing rate protection starting component, and the mean value K of first current-rising-rate 1RMean value K with second current-rising-rate 2RDuring all less than the startup threshold value, the current reference value of real-time update current changing rate protection, current reference value power taking levelling average I FRInstantaneous value; Mean value K when first current-rising-rate 1RMean value K with second current-rising-rate 2RAmong both any one stops to upgrade when starting threshold value.
6) current increment protection: when current changing rate protection starting component starts, get into increment protection logic determines, increment protection logic determines adopts current average I FRInstantaneous value subtracts the current reference value as the protection input value; If the protection input value is protected definite value greater than increment, and satisfied time-delay demand, then increment protection action; If the increment protection is moved, then no longer judge starting component, need only judge that whether input value protects definite value less than increment, if satisfy condition, then returns.
7) current-rising-rate protection: when current changing rate protection starting component starts; Get into current-rising-rate protection logic determines; Current-rising-rate picks up counting, if current changing rate protection starting component starts and last till the time-delay of adjusting of current-rising-rate protection always, then adopts current average I FRInstantaneous value subtracts the current reference value as input value, if the smallest incremental definite value that input value is protected greater than current-rising-rate, then current-rising-rate protection action.
CN201210357913.XA 2012-09-24 2012-09-24 Realizing method of current change rate protection Active CN102842890B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210357913.XA CN102842890B (en) 2012-09-24 2012-09-24 Realizing method of current change rate protection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210357913.XA CN102842890B (en) 2012-09-24 2012-09-24 Realizing method of current change rate protection

Publications (2)

Publication Number Publication Date
CN102842890A true CN102842890A (en) 2012-12-26
CN102842890B CN102842890B (en) 2015-03-25

Family

ID=47370027

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210357913.XA Active CN102842890B (en) 2012-09-24 2012-09-24 Realizing method of current change rate protection

Country Status (1)

Country Link
CN (1) CN102842890B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103364622A (en) * 2013-08-08 2013-10-23 株洲硬质合金集团有限公司 Catching method for breakdown points in electrolytic capacitor breakdown test process
CN103746357A (en) * 2013-12-27 2014-04-23 北京石油化工学院 A relay protection method of a DC traction power grid
CN105896564A (en) * 2016-06-03 2016-08-24 国网冀北节能服务有限公司 Series compensation device and safety protection method of same
CN106786430A (en) * 2017-01-19 2017-05-31 上海电气集团股份有限公司 A kind of discriminating conduct of short circuit malfunction
CN112332381A (en) * 2020-09-21 2021-02-05 天津中铁电气化设计研究院有限公司 Overhead earth wire short-circuit protection method based on current change rate and increment
CN112467688A (en) * 2020-11-11 2021-03-09 广州市扬新技术研究有限责任公司 Current increment protection method of direct current traction protection measurement and control device
CN115166340A (en) * 2022-09-06 2022-10-11 中铁电气化勘测设计研究院有限公司 Processing method of sampling data of subway direct current protection device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050057875A1 (en) * 2003-09-15 2005-03-17 Joseph Kover Programmable system to dynamically control the DC current through a load circuit and/or the DC voltage across a load circuit
CN102185288A (en) * 2011-05-18 2011-09-14 北京四方继保自动化股份有限公司 Current increment protection method capable of reflecting mean current change rate
CN102195275A (en) * 2011-05-18 2011-09-21 北京四方继保自动化股份有限公司 Protecting method of current change rate with self-adaptive output time-delay characteristic

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050057875A1 (en) * 2003-09-15 2005-03-17 Joseph Kover Programmable system to dynamically control the DC current through a load circuit and/or the DC voltage across a load circuit
CN102185288A (en) * 2011-05-18 2011-09-14 北京四方继保自动化股份有限公司 Current increment protection method capable of reflecting mean current change rate
CN102195275A (en) * 2011-05-18 2011-09-21 北京四方继保自动化股份有限公司 Protecting method of current change rate with self-adaptive output time-delay characteristic

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张目然: "电流上升率和电流增量保护的定量测量", 《都市快轨交通》, vol. 20, no. 3, 30 June 2007 (2007-06-30), pages 46 - 48 *
王宏: "地铁牵引供电系统直流馈线保护", 《电气化铁道》, no. 4, 30 April 2002 (2002-04-30), pages 41 - 44 *
郎兵: "直流牵引供电系统馈线保护的研究", 《北京交通大学学报》, vol. 33, no. 5, 31 October 2009 (2009-10-31), pages 65 - 68 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103364622B (en) * 2013-08-08 2015-09-16 株洲硬质合金集团有限公司 A kind of method for catching of electrochemical capacitor flash test breakdown point
CN103364622A (en) * 2013-08-08 2013-10-23 株洲硬质合金集团有限公司 Catching method for breakdown points in electrolytic capacitor breakdown test process
CN103746357A (en) * 2013-12-27 2014-04-23 北京石油化工学院 A relay protection method of a DC traction power grid
CN105896564B (en) * 2016-06-03 2018-06-01 国网冀北节能服务有限公司 The method for security protection of series compensation device and series compensation device
CN105896564A (en) * 2016-06-03 2016-08-24 国网冀北节能服务有限公司 Series compensation device and safety protection method of same
CN106786430B (en) * 2017-01-19 2018-12-18 上海电气集团股份有限公司 A kind of discriminating conduct of short circuit malfunction
CN106786430A (en) * 2017-01-19 2017-05-31 上海电气集团股份有限公司 A kind of discriminating conduct of short circuit malfunction
CN112332381A (en) * 2020-09-21 2021-02-05 天津中铁电气化设计研究院有限公司 Overhead earth wire short-circuit protection method based on current change rate and increment
CN112332381B (en) * 2020-09-21 2022-08-30 中铁电气化勘测设计研究院有限公司 Overhead earth wire short-circuit protection method based on current change rate and increment
CN112467688A (en) * 2020-11-11 2021-03-09 广州市扬新技术研究有限责任公司 Current increment protection method of direct current traction protection measurement and control device
CN112467688B (en) * 2020-11-11 2023-03-24 广州市扬新技术研究有限责任公司 Current increment protection method of direct current traction protection measurement and control device
CN115166340A (en) * 2022-09-06 2022-10-11 中铁电气化勘测设计研究院有限公司 Processing method of sampling data of subway direct current protection device
CN115166340B (en) * 2022-09-06 2023-01-10 中铁电气化勘测设计研究院有限公司 Processing method of sampling data of subway direct current protection device

Also Published As

Publication number Publication date
CN102842890B (en) 2015-03-25

Similar Documents

Publication Publication Date Title
CN102842890A (en) Realizing method of current change rate protection
US10310005B2 (en) Time-domain distance line protection of electric power delivery systems
CN104466903B (en) Current transformer saturation recognition method for differential protection
CN103235213A (en) Differentiation method of lightening back flashover and lightening shielding failure on electric transmission line
EP2502319A1 (en) Method and apparatus for ground distance protection
CN104242267A (en) Distance protection method for wind power generation transmission line
CN105914714A (en) Method based on current waveform recognition for preventing failure protection maloperation of circuit breaker
CN106771877A (en) The determination method and apparatus of the position of failure point of system with non effectively earth ed neutral
CN105743073A (en) Improved inverse time limit over-current protection method
CN102707193A (en) Ultra-low-power implementing method for fault indicator
CN104466920A (en) Circuit breaker failure protection method
CN103986133B (en) The time delay over-current guard method started based on reactive power Sudden Changing Rate
CN103219696A (en) Overcurrent protection method and device
CN201266832Y (en) Protection observe and control device for high-voltage motor
CN104198889A (en) Successive action based single-terminal location method for instant ground fault of high-voltage line
CN101800414A (en) Method for judging and protecting low voltage in transformer short circuit fault
CN104391191B (en) A kind of high resistive fault detection method of utilization line mould current-differencing quadratic sum
CN103715667B (en) A kind of electrical power trans mission/distribution system guard method and transmission & distribution electrical protective device
US11162994B2 (en) Fault current calculation during transformer saturation using the waveform unsaturated region
CN101534009B (en) Analysis method for operating characteristics of power-frequency variable directional protection of AC-DC interconnected power grid
CN104578006B (en) The accident overload blocking method for domain distance protection III sections of standing
CN112564065A (en) Pseudo-synchronous differential protection method and system for direct-current transmission line
CN104953565B (en) Line interphase fault relay protection method is surveyed based on transition resistance
CN102856884A (en) Heat accumulation based frequency converter overload method
CN105071356A (en) Method for preventing sampling channel AD over-limit differential misoperation

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