CN103245846B - Dynamic zero drift filtering algorithm for relay protection - Google Patents
Dynamic zero drift filtering algorithm for relay protection Download PDFInfo
- Publication number
- CN103245846B CN103245846B CN201210029910.3A CN201210029910A CN103245846B CN 103245846 B CN103245846 B CN 103245846B CN 201210029910 A CN201210029910 A CN 201210029910A CN 103245846 B CN103245846 B CN 103245846B
- Authority
- CN
- China
- Prior art keywords
- drift
- value
- drift value
- zero drift
- alternation
- 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
- 238000001914 filtration Methods 0.000 title abstract description 4
- 238000005070 sampling Methods 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 14
- 230000001681 protective effect Effects 0.000 claims description 15
- 238000005259 measurement Methods 0.000 description 4
- 230000002411 adverse Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Landscapes
- Emergency Protection Circuit Devices (AREA)
Abstract
The invention provides a dynamic zero drift filtering algorithm for relay protection, which comprises the following steps: (1), inputting a sampling passage into a short circuit, so as to measure an initial zero drift value which is then solidified into a memorizer of the protection device; (2), calculating the sum of sampling point data; (3), calculating an average zero drift values according to the sum of the sampling point data; (4), calculating a zero drift variation between the initial zero drift value and the average zero drift value; (5), calculating a zero drift gradation step length according to a gradation period; and (6), progressively adjusting the zero drift value of each sampling point according to the gradation period. Due to the adoption of the algorithm , the zero drift value can be automatically calculated during the operation of a device, so as to eliminate the error caused by that the zero drift varies with time and temperature variation in the practical operation process, and further to improve the protection and metering precision.
Description
Technical field
The present invention relates to technical field of electric power, particularly relate to the algorithm of a kind of kinetic measurement in relay protection filtering drift.
Background technology
Drift is the abbreviation of zero point drift (Zero drift), and in relay protection, drift refers to when additional interchange is zero, the size of protective device sampling quantity.It is because the temperature drift of the devices such as sensor, metering circuit, attenuator, amplifier, wave filter and time drift cause.
Operational amplifier is all adopt direct-coupled mode, and direct coupling type amplifying circuit Q point at different levels is interactional, and due to amplification at different levels, the faint change of the first order, can make the change that output stage generation is very large.When input short, export and will slowly change in time.
The parameter of transistor, by the impact of temperature, also can produce zero point drift.
The precision of metering can be affected zero point when there is drift, also can bring error to protection.Drift is the one of the main reasons causing systematic error, and some fast algorithms and instantaneous value protection algorism are more easily subject to the impact of drift.
Existing protective relaying device, generally measures drift at the device before measurement that dispatches from the factory, is solidificated in protective device as a definite value.During actual motion, when protective device is sampled at every turn, all deduct the drift value that this is fixing.But the drift of protective device is not a changeless value, its can fluctuate along with the change of environment temperature, working time up and down, and comparatively greatly, degree of accuracy is lower for the sampled data therefore causing deducting fixing drift value and True Data error.
Summary of the invention
The invention provides a kind of dynamic filter drift algorithm for relay protection, in plant running process, automatically calculate drift value, to eliminate in actual moving process drift in time with temperature variation and change the error brought, improve the precision of protection and metering.
In order to achieve the above object, the invention provides following technical scheme:
For a dynamic filter drift algorithm for relay protection, comprise the following steps:
(1) sampling channel is inputted short circuit, measure the initial drift value of passage and be cured in the storer of protective device;
(2) calculating sampling point data sum;
(3) according to sampling number according to sum, obtain average drift value;
(4) the drift variable quantity between average drift value and initial drift value is calculated;
(5) drift alternation step-length is gone out according to alternation computation of Period;
(6) according to the described alternation cycle, alternation adjusts the drift value of each sampled point.
Improve as one, further comprising the steps of:
(7) the drift value after adjustment is cured in the storer of protective device by timing.
The process of described calculating sampling point data sum is: calculate
individual sampling number according to sum,
, wherein
for single sampled data,
,
for each cycle sampling number,
for the multiple of total Cycle Length divided by single Cycle Length of sampling.
Preferably,
value be integer.
Preferably, the process obtaining average drift value described in is:
, wherein
for average drift value.
As preferably, the process of described calculating drift variable quantity is:
, wherein
for initial drift value,
for drift value variable quantity.
As preferably, the process of described calculating drift alternation step-length is:
, wherein
for drift alternation step-length,
,
for needing the amount of cycles of adjustment.
As preferably,
value be integer.
As preferably, the process that described alternation formula adjusts the drift value of each sampled point is: the drift value adjusting each sampled point, makes
.
Improve as one, when breaking down, close described dynamic filter drift algorithm.
Dynamic filter drift algorithm is simply effective, fast response time, and reliability is high, not by the impact of frequency change, has good applicability.In addition without the need to hardware costs, overcome the deficiency of hard ware measure method, eliminate the error that hardware frequency measurement is introduced due to the factor of device, improve the precision of frequency measurement.
Accompanying drawing explanation
Fig. 1 is dynamic filter drift algorithm flow schematic diagram provided by the invention.
Fig. 2 is that algorithm application provided by the invention is at the protection of protective relaying device and the schematic flow sheet in measuring.
specific implementation method
Below with reference to specific embodiment, technical scheme provided by the invention is described in detail:
Embodiment:
As shown in Figure 1, dynamic conditioning drift process is realized by following step:
(1) sampling channel is inputted short circuit, measure the initial drift value of passage
and be cured in the storer ROM of protective device, define initial drift value
.
(2) be AC signal due to what measure in relay protection, the integrated value of this signal within each cycle is zero, and therefore within each cycle, signal data sum is drift value sum.Get the integral multiple that computation period length is individual signals Cycle Length during calculating, multiple is
, this ensure that in computation period, interchange value accumulated value is zero,
can need to arrange according to the precision calculated.If
be the sampling number in the signal period, definition
, suppose that drift value is
, then in whole computation period, calculating sampling point AC signal data sum
.
(3) according to sampling number according to sum
with sampled point quantity
, average drift value in whole computation period can be calculated
.
(4) the drift variable quantity between average drift value and initial drift value is calculated: establish
for drift variable quantity, calculate
.
(5) go out drift alternation step-length according to alternation computation of Period: by adjustment drift period definition be the alternation cycle, owing to will progressively adjust drift value within the alternation cycle, preferably get the integral multiple that alternation Cycle Length is single Cycle Length, multiple will be
.?
when being worth constant,
value determine each amplitude of adjustment and the time of adjustment,
larger, the time of adjustment is longer, and the number of times of adjustment is also more, and the amplitude of corresponding each adjustment is less,
less, the time of adjustment is shorter, and the amplitude of each adjustment is also corresponding larger,
can arrange the different demands of adjusting range and regulation time according in practical operation.If
if,
for drift alternation step-length, i.e. the amplitude of each adjustment, draws
.
(6) according to the described alternation cycle, alternation formula adjusts the drift value of each sampled point: the object of alternation formula adjustment drift is the change of level and smooth drift, reduces the impact on protection and metering.Although drift slowly changes with temperature in time, the variable quantity in computation period is very little, average through the alternation cycle, and adjusting range is each time very little, does not almost affect calculating.In this example, in former drift value
on to add up drift alternation step-length according to following formula:
, from above formula, pass through
secondary adjustment, each adjusting range is
, make when alternation end cycle, the drift value set be adjusted to the actual drift value calculated in former internal memory
.
(7) preferably the drift value after adjustment is cured in the storer of protective device by timing, has new drift initial value when making again Pneumatic protecting device in storer ROM.The frequency of timing write drift value can be arranged as required.
Because when breaking down, waveform can produce sudden change, the sampled data sum now in computation period no longer only comprises drift value, and now dynamically filtering drift algorithm will be no longer accurate, preferably close dynamically filter drift algorithm when therefore fault occurs.
As shown in Figure 2, after the above-mentioned dynamic filter drift algorithm of connected applications, the protection in protective relaying device and measuring process as follows:
First, the initialization of protective relay, sampling module, drift module is carried out;
Secondly application dynamic filter drift algorithm adjustment drift provided by the invention;
New drift value after adjustment be applied to protection and measure.If as broken down in protection and measuring process, then close dynamically filter drift algorithm, use the drift value be solidificated in protective device to calculate.
Apply algorithm provided by the invention, can on the basis not affecting protection and metering, the drift of dynamic tracking means, eliminates the adverse effect that drift causes system with temperature variation in time.
Claims (3)
1., for a dynamic filter drift algorithm for relay protection, it is characterized in that comprising the following steps:
(1) sampling channel is inputted short circuit, measure the initial drift value of passage and be cured in the storer of protective device;
(2) calculating sampling point data sum, detailed process is:
Calculate
individual sampling number according to sum,
, wherein
for single sampled data,
,
for each cycle sampling number,
for sample total Cycle Length divided by single Cycle Length multiple and for integer;
(3) according to sampling number according to sum, obtain average drift value according to following formula:
，
Wherein
for average drift value,
for sampling number is according to sum,
for sampled point quantity;
(4) by the drift variable quantity between the average drift value of following formulae discovery and initial drift value:
Wherein
for initial drift value,
for drift value variable quantity;
(5) go out drift alternation step-length according to alternation computation of Period, computing formula is:
Wherein
for drift alternation step-length,
,
for need adjustment amount of cycles and be integer;
(6) according to the described alternation cycle, alternation formula adjusts the drift value of each sampled point, makes
.
2. the dynamic filter drift algorithm for relay protection according to claim 1, characterized by further comprising following steps:
(7) the drift value after adjustment is cured in the storer of protective device by timing.
3. the dynamic filter drift algorithm for relay protection according to claim 1 and 2, is characterized in that: when breaking down, and closes described dynamic filter drift algorithm.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210029910.3A CN103245846B (en) | 2012-02-10 | 2012-02-10 | Dynamic zero drift filtering algorithm for relay protection |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210029910.3A CN103245846B (en) | 2012-02-10 | 2012-02-10 | Dynamic zero drift filtering algorithm for relay protection |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103245846A CN103245846A (en) | 2013-08-14 |
| CN103245846B true CN103245846B (en) | 2015-03-11 |
Family
ID=48925491
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210029910.3A Active CN103245846B (en) | 2012-02-10 | 2012-02-10 | Dynamic zero drift filtering algorithm for relay protection |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103245846B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104678154B (en) * | 2015-03-13 | 2017-10-20 | 南京大全自动化科技有限公司 | A kind of method that regulation power distribution equipment zero based on adaptive nulling technology floats |
| CN106546934A (en) * | 2015-09-18 | 2017-03-29 | 比亚迪股份有限公司 | Zero shift rectifying method and apparatus when current of electric is sampled |
| CN107134760B (en) * | 2017-07-13 | 2019-03-12 | 积成软件有限公司 | It is a kind of be biased to based on bus protection spill current zero simple drift filtering method |
| CN113238178A (en) * | 2021-05-14 | 2021-08-10 | 武汉理工大学 | Zero drift signal elimination method of alternating current signal sensor |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU809493A1 (en) * | 1979-02-06 | 1981-02-28 | Ивано-Франковский Институт Нефтии Газа | M-dm dc amplifier |
| CN1080726A (en) * | 1993-05-08 | 1994-01-12 | 浙江大学 | The zero blanking float method of micro current instrument and device thereof |
| CN2434684Y (en) * | 2000-08-10 | 2001-06-13 | 中国科学院等离子体物理研究所 | Long time non-zero floation photoelectric isolator |
-
2012
- 2012-02-10 CN CN201210029910.3A patent/CN103245846B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU809493A1 (en) * | 1979-02-06 | 1981-02-28 | Ивано-Франковский Институт Нефтии Газа | M-dm dc amplifier |
| CN1080726A (en) * | 1993-05-08 | 1994-01-12 | 浙江大学 | The zero blanking float method of micro current instrument and device thereof |
| CN2434684Y (en) * | 2000-08-10 | 2001-06-13 | 中国科学院等离子体物理研究所 | Long time non-zero floation photoelectric isolator |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103245846A (en) | 2013-08-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102998525B (en) | Single-phase multifunctional electric energy metering circuit and ammeter automatic calibrating method thereof | |
| CN103245846B (en) | Dynamic zero drift filtering algorithm for relay protection | |
| CN103472259B (en) | Method for silicon micro-resonant type accelerometer temperature compensation | |
| CN101488031A (en) | High-precision magnetic bearing axial control method based on interference observer | |
| CN105043348A (en) | Accelerometer gyroscope horizontal angle measurement method based on Kalman filtering | |
| CN103814276A (en) | Method and apparatus for determining differential flow characteristics of multiple-meter fluid flow system | |
| CN207502605U (en) | A kind of frequency measuring system based on VCXO references | |
| CN105387859A (en) | Temperature drift error compensation method of MEMS (Micro-electromechanical System) sensor group | |
| CN107063311A (en) | A kind of dynamic measurement system method of single coil vibrating sensor | |
| CN110186522B (en) | Moisture overreading compensation and flow measurement method combining vortex street amplitude characteristic | |
| CN101706563A (en) | Digital calibration system of partial discharge test equipment | |
| CN104199307B (en) | Hardware-in-loop simulation method and system | |
| CN204694288U (en) | A kind of gas-solid diphasic stream parameter measurement mechanism | |
| CN105298889A (en) | Gas compressor surge detection method | |
| CN202886471U (en) | Multi-channel leakage current phase difference detection equipment | |
| CN106772193B (en) | Measuring method using current transformer frequency characteristic measuring device | |
| CN105319434B (en) | A kind of method and device of electric energy meter measurement electric current | |
| CN103206997A (en) | Gas Roots flow measuring device and flow correction method for same | |
| CN103206995B (en) | Be applicable to gas flow surveying instrument and the method for testing of factory building condition | |
| CN102737145A (en) | Measurement-based modeling method by prediction on electromagnetic emission broadband behavioral level of electronic component | |
| CN103163924A (en) | System on chip (SoC) chip and compensation method and compensation device of operational amplifier offset voltage thereof | |
| CN102128860A (en) | Dynamic parameter measuring system for humidity sensor | |
| CN110098610A (en) | The real-time identification method and system of power system oscillation dominant pattern under fault disturbance | |
| CN102306249B (en) | Equipment efficiency simulation method and system | |
| CN104678154B (en) | A kind of method that regulation power distribution equipment zero based on adaptive nulling technology floats |
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 | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Dynamic zero drift filtering algorithm for relay protection Effective date of registration: 20190717 Granted publication date: 20150311 Pledgee: Nanjing Bank Co.,Ltd. Nanjing North Branch Pledgor: Nanjing Hongyi Electrical Automation Co.,Ltd. Registration number: 2019320000333 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Date of cancellation: 20200813 Granted publication date: 20150311 Pledgee: Nanjing Bank Co.,Ltd. Nanjing North Branch Pledgor: NANJING HONGYI ELECTRICAL AUTOMATION Co.,Ltd. Registration number: 2019320000333 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A Dynamic Zero Drift Filter Algorithm for Relay Protection Effective date of registration: 20231229 Granted publication date: 20150311 Pledgee: Zheshang Bank Co.,Ltd. Nanjing Branch Pledgor: Nanjing Hongyi Electrical Automation Co.,Ltd. Registration number: Y2023980075267 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right |