CN103630792B - A kind of arrester resistive leakage current algorithm - Google Patents
A kind of arrester resistive leakage current algorithm Download PDFInfo
- Publication number
- CN103630792B CN103630792B CN201310300687.6A CN201310300687A CN103630792B CN 103630792 B CN103630792 B CN 103630792B CN 201310300687 A CN201310300687 A CN 201310300687A CN 103630792 B CN103630792 B CN 103630792B
- Authority
- CN
- China
- Prior art keywords
- current
- total current
- value
- resistance property
- compounent
- 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
Abstract
The invention discloses a kind of arrester resistive leakage current algorithm, including step:S1 samples the leakage total current of a cycle, S2 calculates total current DC component, S3 calculates total current AC compounent, S4 is to total current AC compounent bandpass filtering, S5 rearranges total current AC compounent waveform since 0 degree, S6 finds knee value and its position in preceding 1/4 cycle, S7 calculates current in resistance property, S8 finds peak value and its position of current in resistance property, S9 calculates the virtual value of current in resistance property, S10 calculates circuital virtual value, S11 proof of algorithm and S12 send out total current virtual value and current in resistance property peak value.The present invention can realize the high-acruracy survey of resistive leakage current, have the advantages that error is small, cost is low, fireballing, can cross preferably reflection arrester working status.
Description
Technical field
The present invention relates to a kind of lightning arrester resistance of a kind of arrester resistive leakage current algorithm, more particularly to electric system
Leakage current algorithm.
Background technology
Zinc-Oxide Arrester be the protection that is serially connected in beside important power equipment between high-voltage power line and the earth this set
For from overvoltage（Lightning stroke）A kind of important protective device of infringement.Due to the superior nonlinear wind vibration of Zinc-Oxide Arrester
（Piezo-resistive properties）With good through-current capability, quickly lightning stroke energy can be discharged by discharging over the ground, it is vast
Power department receives and widely uses.
Shown by operating experience for many years, since the operation of arrester is had been under high voltage, can made inside MOA
Resistor disc aging, simultaneously because the influence of environmental condition, its resistor disc makes moist deterioration, causes the increasing of resistive leakage current
Add.If cannot be out of service in time by abnormal arrester rapidly, within a period of time（Several month, days or a few hours）, when being
Occurs overvoltage in system, the possible thermal runaways of MOA, explode, trigger large area electric power accident.
The leakage total current of arrester contains two parts of capacitive leakage current and resistive leakage current.It is wherein resistive to let out
Leakage current is the technical parameter of real reflection arrester operating status.Since the detection of resistive leakage current must conventionally sample
Total current and voltage are leaked, by calculating voltage and leaking circuital angle, then extrapolates resistance using trigonometric function relation
Property leakage current.Resistive leakage current only accounts for leakage circuital 10%~15% under normal conditions, is calculated by multiplication twice
The resistive leakage current value error gone out is bigger.Cause two critical defects of common detection methods：Accuracy of detection is poor,
It can not truly reflect arrester working status.Testing cost is higher, due to needing the high pressure PT using sampled voltage, can not push away
Extensively.
Monitoring arrester leaks the size of current in resistance property in real time under circuit operating status, truly reflects under working status
Lightning arrester insulation characteristic, be the effective means for preventing surge arrester failure from causing large-area power-cuts.
Yet with the following characteristics of leakage current signal：It is difficult correct measurement to cause leakage current：
1st, have periodically and non-linear（It is sinusoidal）Characteristic, it is difficult to described in the time domain using mathematical formulae.
2nd, value range is wide and absolute figure is small（0.05~1.00mA）, it is difficult to amplified with standard amplifying circuit.
3rd, scene is subject to harmonic wave, electromagnetism and white noise acoustic jamming, it is difficult to is removed using common interference protection measure.
4th, the finite energy of generation, it is difficult to measured with common method, but still comply with law of conservation of energy.
How safe and reliable and accurate survey calculation arrester leaks current in resistance property, is real under circuit operating status
When monitor leakage current of an arrester basic means.
The content of the invention
It is an object of the invention to provide a kind of arrester resistive leakage current algorithm.Realize the height of resistive leakage current
Precision measure, has the advantages that error is small, cost is low, fireballing, can cross preferably reflection arrester working status
Realizing the technical solution of above-mentioned purpose is：A kind of arrester resistive leakage current algorithm, including step:S1 samplings one
The leakage total current in a cycle, S2 calculate total current DC component, S3 calculates total current AC compounent, S4 is to total current exchange point
Amount bandpass filtering, S5 rearranged since 0 degree total current AC compounent waveform, S6 found in preceding 1/4 cycle knee value and
Its position, S7 calculate current in resistance property, S8 finds the peak value of current in resistance property and its position, S9 calculate virtual value, the S10 of current in resistance property
Calculate circuital virtual value, S11 proof of algorithm and S12 and send out total current virtual value and current in resistance property peak value.
Above-mentioned steps S1 samples the leakage total current of a cycle, samples the waveform of a cycle N points altogether, N >=128, obtain
To discrete total current function Ix (n), 0~N of n values.
The formula that above-mentioned steps S2 calculates total current DC component is that Vdc=(∑ Ix (n))/N, Vdc is total current direct current point
Amount.
The formula that above-mentioned steps S3 calculates total current AC compounent is Ix (n)AC=Ix (n)-Vdc, Ix (n)ACFor total current
AC compounent.
Formula of the above-mentioned steps S4 to total current AC compounent bandpass filtering：Ix(n)BP=(∑Ix(n-i)AC×ai)/25,i
=(0~24), a0~a24Value be（4、5、7、11、16、22、28、34、39、44、48、50、51、50、48、44、39、34、
28、22、16、11、7、5、4）, Ix (n)BPFor the total current AC compounent after bandpass filtering.
It is Ix (n) that above-mentioned steps S5 rearranges total current AC compounent waveform since 0 degreeBP, rearranging waveform
Splice point carry out 1/2 smoothing processing.
Above-mentioned steps S6 finds knee value in preceding 1/4 cycle can be used as capacity current maximum for Ic and its position Nc, Ic
Value.
Above-mentioned steps S7 calculate current in resistance property formula be：Ir(n)=Ix(n)BP- Ic × Sin (N-Nc) × 2 π/N, Ir (n)
For resistive leakage current waveform.
Above-mentioned steps S8, which finds the peak I r (max) of current in resistance property and its position Nr, Ir (max), to be used as current in resistance property most
Big value.
The formula that above-mentioned steps S9 calculates the virtual value of current in resistance property is Ir2=(∑Ir(n)2)/N, Ir are current in resistance property
Virtual value.
The formula that above-mentioned steps S10 calculates circuital virtual value is Ix2=(∑Ix(n)BP 2)/N, Ix have to be circuital
Valid value.
The method of above-mentioned steps S11 proof of algorithm is to meet following two formula at the same time：Ir2+Ic2/2=Ix2, Nr-Nc=N/
4。
The beneficial effects of the invention are as follows：Realize the high-acruracy survey of resistive leakage current, have error is small, cost is low,
Fireballing advantage, can cross preferably reflection arrester working status.
Brief description of the drawings
Fig. 1 is the flow diagram of the present invention.
Embodiment
The present invention will be further described with reference to the accompanying drawings.
Referring to Fig. 1, a kind of arrester resistive leakage current algorithm is given in figure, including step:S1 samples a week
The leakage total current of phase, S2 calculate total current DC component, S3 calculates total current AC compounent, S4 to total current AC compounent band
Pass filter, S5 rearrange total current AC compounent waveform, S6 since 0 degree and find knee value and its position in preceding 1/4 cycle
Put, S7 calculates current in resistance property, S8 finds the peak value of current in resistance property and its position, S9 calculate the virtual value of current in resistance property, S10 is calculated
Circuital virtual value, S11 proof of algorithm and S12 send out total current virtual value and current in resistance property peak value.
Step S1 samples the leakage total current of a cycle, samples the waveform of a cycle N points altogether, and N takes 512, obtain from
Dissipate total current function Ix (n), 0~N of n values.
The formula that step S2 calculates total current DC component is that Vdc=(∑ Ix (n))/N, Vdc is total current DC component.
The formula that step S3 calculates total current AC compounent is Ix (n)AC=Ix (n)-Vdc, Ix (n)ACExchanged for total current
Component.
Formula of the step S4 to total current AC compounent bandpass filtering：Ix(n)BP=(∑Ix(n-i)AC×ai)/25,i=(0
~24), a0~a24Value be（4、5、7、11、16、22、28、34、39、44、48、50、51、50、48、44、39、34、28、
22、16、11、7、5、4）, Ix (n)BPFor the total current AC compounent after bandpass filtering.
It is Ix (n) that step S5 rearranges total current AC compounent waveform since 0 degreeBP, rearranging the spelling of waveform
Contact carries out 1/2 smoothing processing.
It is that Ic and its position Nc, Ic can be used as capacity current maximum that step S6 finds knee value in preceding 1/4 cycle.
Step S7 calculate current in resistance property formula be：Ir(n)=Ix(n)BP- Ic × Sin (N-Nc) × 2 π/N, Ir (n) is resistance
Property leakage current waveform.
Step S8, which finds the peak I r (max) of current in resistance property and its position Nr, Ir (max), can be used as current in resistance property maximum
Value.
The formula that step S9 calculates the virtual value of current in resistance property is Ir2=(∑Ir(n)2)/N, Ir are the effective of current in resistance property
Value.
The formula that step S10 calculates circuital virtual value is Ix2=(∑Ix(n)BP 2)/N, Ix are circuital effective
Value.
The method of step S11 proof of algorithm is to meet following two formula at the same time：Ir2+Ic2/2=Ix2, Nr-Nc=N/4.
The principle of the present invention is：Using basic trigonometric function replacement relational expression, corresponding calculation formula and algorithm are changed,
In the case where not using Voltage Reference amount, measurement detection resistive leakage current.Realize the high accuracy survey of resistive leakage current
Amount, has the advantages that error is small, cost is low, fireballing, can cross preferably reflection arrester working status.
The present invention is described in detail above in association with attached drawing embodiment, those skilled in the art can be according to upper
State and bright many variations example is made to the present invention.Thus, some details in embodiment should not form limitation of the invention, this
Invention will be used as protection scope of the present invention using the scope that the appended claims define.
Claims (1)
1. a kind of arrester resistive leakage current algorithm, it is characterised in that including step:The leakage of S1 sampling a cycles is entirely electric
Stream, S2 calculate total current DC component, S3 calculates total current AC compounent, S4 to total current AC compounent bandpass filtering, S5 from 0
Degree starts to rearrange total current AC compounent waveform, S6 finds knee value in preceding 1/4 cycle and its position, S7 calculate it is resistive
Electric current, S8 find the peak value of current in resistance property and its position, S9 calculate the virtual value of current in resistance property, S10 calculate it is circuital effectively
Value, S11 proof of algorithm and S12 send out total current virtual value and current in resistance property peak value；
The step:S1 sample a cycle leakage total current, altogether sample a cycle N points waveform, N >=128, obtain from
Dissipate total current function Ix (n), 0~N of n values；
The step：The formula that S2 calculates total current DC component is that Vdc=(∑ Ix (n))/N, Vdc is total current direct current point
Amount；
The step：The formula that S3 calculates total current AC compounent is Ix (n)AC=Ix (n)-Vdc, Ix (n)ACExchanged for total current
Component；
The step：Formula Ixs (n) of the S4 to total current AC compounent bandpass filteringBP=(∑ Ix (n-i)AC×ai)/25, i=
(0~24), a0~a24Value for (4,5,7,11,16,22,28,34,39,44,48,50,51,50,48,44,39,34,28,
22nd, 16,11,7,5,4), Ix (n)BPFor the total current AC compounent after bandpass filtering；
The step：It is Ix (n) that S5 rearranges total current AC compounent waveform since 0 degreeBP, rearranging the spelling of waveform
Contact carries out 1/2 smoothing processing；
The step：It is that Ic and its position Nc, Ic can be used as capacity current maximum that S6 finds knee value in preceding 1/4 cycle；
The step：The formula that S7 calculates current in resistance property is Ir (n)=Ix (n)BP- Ic × Sin (N-Nc) × 2 π/N, Ir (n) be
Resistive leakage current waveform；
The step：S8, which finds the peak I r (max) of current in resistance property and its position Nr, Ir (max), can be used as current in resistance property maximum
Value；
The step：The formula that S9 calculates the virtual value of current in resistance property is Ir2=(∑ Ir (n)2)/N, Ir are having for current in resistance property
Valid value；
The step：The formula that S10 calculates circuital virtual value is Ix2=(∑ Ix (n)BP 2)/N, Ix are circuital effective
Value；
The step：The method of S11 proof of algorithm is to meet following two formula at the same time：Ir2+Ic2/ 2=Ix2, Nr-Nc=N/
4。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310300687.6A CN103630792B (en) | 2013-07-17 | 2013-07-17 | A kind of arrester resistive leakage current algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310300687.6A CN103630792B (en) | 2013-07-17 | 2013-07-17 | A kind of arrester resistive leakage current algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103630792A CN103630792A (en) | 2014-03-12 |
| CN103630792B true CN103630792B (en) | 2018-04-17 |
Family
ID=50212046
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310300687.6A Expired - Fee Related CN103630792B (en) | 2013-07-17 | 2013-07-17 | A kind of arrester resistive leakage current algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103630792B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109085448A (en) * | 2018-07-31 | 2018-12-25 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | A kind of converter station equipment leakage current waveform measurement and composition analysis device |
| CN109900984A (en) * | 2019-03-01 | 2019-06-18 | 安徽天玄智能科技有限公司 | Monitoring method, system and the device of surge protector based on current in resistance property |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09196988A (en) * | 1996-01-19 | 1997-07-31 | Mitsubishi Electric Corp | Deterioration detector for zinc oxide type lightning arrester |
| CN1892240A (en) * | 2005-07-08 | 2007-01-10 | 上海市电力公司 | On-line monitoring method for dissipative current value in leakage current of lightning protector |
| CN101359016A (en) * | 2007-07-31 | 2009-02-04 | 上海电气自动化设计研究所有限公司 | Method for detecting resistance leakage current of lightning arrester without voltage reference |
| CN101533041A (en) * | 2009-04-10 | 2009-09-16 | 李新育 | Method for online separately detecting resistive leakage current of non-voltage reference source |
| CN102854431A (en) * | 2011-07-01 | 2013-01-02 | 上海思盛信息科技事务所 | Method for online detecting leakage resistive current of lightning arrestor |
| CN102901856A (en) * | 2012-09-17 | 2013-01-30 | 吉林省电力有限公司长春供电公司 | Cable line arrester resistive current detection method based on phase search |
-
2013
- 2013-07-17 CN CN201310300687.6A patent/CN103630792B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09196988A (en) * | 1996-01-19 | 1997-07-31 | Mitsubishi Electric Corp | Deterioration detector for zinc oxide type lightning arrester |
| CN1892240A (en) * | 2005-07-08 | 2007-01-10 | 上海市电力公司 | On-line monitoring method for dissipative current value in leakage current of lightning protector |
| CN101359016A (en) * | 2007-07-31 | 2009-02-04 | 上海电气自动化设计研究所有限公司 | Method for detecting resistance leakage current of lightning arrester without voltage reference |
| CN101533041A (en) * | 2009-04-10 | 2009-09-16 | 李新育 | Method for online separately detecting resistive leakage current of non-voltage reference source |
| CN102854431A (en) * | 2011-07-01 | 2013-01-02 | 上海思盛信息科技事务所 | Method for online detecting leakage resistive current of lightning arrestor |
| CN102901856A (en) * | 2012-09-17 | 2013-01-30 | 吉林省电力有限公司长春供电公司 | Cable line arrester resistive current detection method based on phase search |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103630792A (en) | 2014-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105137163B (en) | The resistive leakage current on-line monitoring method of metal oxide piezoresistive arrester or Surge Protector | |
| KR101086878B1 (en) | Method for measuring resistive leakage current of nonlinear device lightning arrester with differential and apparatus thereof | |
| CN102901856B (en) | Cable run arrester resistance current detection method based on phase search | |
| CN102830319B (en) | A kind of zinc oxide lightning arrester insulation state measuring device with electricity and method | |
| CN103364696A (en) | Intelligent distinguishing method of insulation characteristic of zinc oxide lightning arrester | |
| CN105954632A (en) | Zinc oxide lightning arrester on-line monitoring and diagnostic method | |
| CN104360197B (en) | Method for monitoring resistance of valve discs of arresters based on resistance-capacitance network | |
| CN108897951A (en) | Arrester nonlinear wind vibration " turn point " and its " bending coefficient " approximating method | |
| CN104820127A (en) | Metal oxide arrester resistance current extracting method and device | |
| CN102854431B (en) | Method for online detecting leakage resistive current of lightning arrestor | |
| CN205374588U (en) | Leakage current collection system | |
| CN103630792B (en) | A kind of arrester resistive leakage current algorithm | |
| KR100498927B1 (en) | Method and Device for measuring resistive leakage current by time-delay synthesis method in deciding the deterioration diagnosis of arrester | |
| KR100930552B1 (en) | Surge protection device and remote monitoring system for it | |
| Yan et al. | Study on the resistive leakage current characteristic of MOV surge arresters | |
| JP2011202956A (en) | Failure determination method of lightning arrester | |
| CN100460883C (en) | Detection method for testing current in resistance property of lightning arrester | |
| Bezerra et al. | Evaluation of surge arrester models for overvoltage studies | |
| KR20120008823A (en) | System for measuring ground impedance in high frequency band | |
| CN108872747B (en) | Surge protector resistive current extraction device and method based on correlation coefficient | |
| CN205792213U (en) | Turbine flow converter capable of resisting transient high voltage | |
| CN104849530A (en) | MOA resistive current cardinal wave measuring method | |
| CN104833841A (en) | Substation intrusion overvoltage detection method | |
| CN208207058U (en) | A kind of surge protector current in resistance property extraction element based on dephased current | |
| WO2009012550A2 (en) | Zinc oxide (zno) varistors degradation analysis system using parameters from the modified longevin model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20151126 Address after: 202172, No. 4, building 3366, 3116 north side highway, Chongming County, Shanghai Applicant after: Shanghai Sisheng Information Technology Office Address before: 201114, room 1, building 88, 1127 West Heng Road, Shanghai, Minhang District Applicant before: Shanghai Sisheng Electric Power Control Technology Co., Ltd. |
|
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180531 Address after: 200025 11 floor, 1 building, 430 Xujiahui Road, Huangpu District, Shanghai. Patentee after: Shanghai Electric Power Communication Co., Ltd. Address before: 202172 3116, room 4, 3366 North Rim highway, Chongming County, Shanghai. Patentee before: Shanghai Sisheng Information Technology Office |
|
| 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: 20180417 Termination date: 20200717 |