CN109298362A - A kind of distribution line failure indicator three-phase synchronous method for testing precision, device and system - Google Patents
A kind of distribution line failure indicator three-phase synchronous method for testing precision, device and system Download PDFInfo
- Publication number
- CN109298362A CN109298362A CN201811051066.8A CN201811051066A CN109298362A CN 109298362 A CN109298362 A CN 109298362A CN 201811051066 A CN201811051066 A CN 201811051066A CN 109298362 A CN109298362 A CN 109298362A
- Authority
- CN
- China
- Prior art keywords
- acquisition unit
- recorder data
- benchmark
- data sequence
- recorder
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
Abstract
The present invention discloses a kind of distribution line failure indicator three-phase synchronous method for testing precision, device and system, when testing, by fault detector three-phase acquisition unit carry on the same phase line of testing stand, and test current waveform is passed through in testing stand primary side, then the recorder data of fault detector three-phase acquisition unit is obtained, related coefficient algorithm is recycled to achieve the purpose that compare time delay between three groups of recorder datas, and then the three-phase synchronous time error of fault detector is calculated, and then it can determine whether three-phase synchronous time error meets the accuracy rating of permission.Using the present invention without changing testing stand sample rate, whether three-phase acquisition unit error synchronization time can be judged in permitted range, and judging result is effectively reliable, testing cost is lower.
Description
Technical field
The present invention relates to power distribution automation equipment technical field, especially a kind of distribution line failure indicator three-phase synchronous
Method for testing precision, device and system.
Background technique
Distribution line failure indicator is mounted on distribution line, for detecting distribution line short and single-phase connecing
The device of earth fault and alert is a kind of simple, practical " two is distant " type power distribution automation solution.Distribution wire
Line fault indicator collects unit by three-phase acquisition unit and one and forms, and acquisition unit and collects between unit using the small nothing of short distance
Line communication mode is collected and is generally carried out data transmission using GPRS public network communication mode between unit and distant place data center.
The three-phase acquisition unit difference carry of distribution line failure indicator is responsible for acquisition, record is matched on three-phase line
Waveform is wirelessly transferred to and collects unit and then synthesize zero sequence transient state letter by electric line three-phase voltage and current waveform
Number, the correctness that three-phase synchronous performance directly affects the accuracy of voltage current waveform record, zero sequence transient signal calculates, into
And it is related to distribution line and occurs correctly trigger failure wave-recording when short circuit or ground fault, and failure accurately determine
Position.
Since three-phase fault indicator wireless channel conditions are not exactly the same, markers clock synchronization has error, therefore three-phase waveform
There are synchronous errors.To meet the zero sequence transient signal required precision for collecting unit synthesis, guarantee distribution line failure indicator
It can correctly indicate that fault location occurs and carries out for failure, it is desirable that error synchronization time of three-phase line acquisition unit should be less than 100 μ
s。
Therefore it needs by distribution line failure indicator test platform (hereinafter referred to as " testing stand "), to distribution line failure
The three-phase synchronous precision of indicator carries out detection test.Current testing stand, by each phase 10kV pressure-wire superimposed current source,
Injection high current waveform is simulated to distribution line failure indicator primary side while being applied in the case where not promoting actual power
High voltage, high current.But the sample rate of testing stand is generally 4000Hz to 5000Hz, and the sampling interval is about 250 μ s or 200 μ
S much larger than the three-phase synchronous time error 100 μ s of maximum value to be allowed, therefore can not correctly judge three-phase synchronous time error
Whether in permitted range.
Summary of the invention
The technical problem to be solved in the present invention are as follows: related coefficient algorithm is utilized, in the premise for not changing testing stand sample rate
Under, realize the three-phase synchronous accuracy test to distribution line failure indicator, and judging result is effectively reliable.
The technical scheme adopted by the invention is as follows: a kind of distribution line failure indicator three-phase synchronous method for testing precision, packet
It includes:
It obtains fault detector three-phase acquisition unit and responds the test current waveform on same route and collected, sampling
Moment corresponding three groups of recorder datas；
Using any acquisition unit as benchmark acquisition unit, other acquisition units are as control acquisition unit, based on setting
Data window width the recorder data of continuous sampling is chosen from the corresponding recorder data of benchmark acquisition unit, obtain benchmark recording
Data sequence；
Data window width based on aforementioned setting, using sliding window, in the sliding scale of setting, every time sliding respectively from
The recorder data for choosing continuous sampling in the corresponding recorder data of acquisition unit is compareed, multiple related recorder data sequences are obtained；
It is right that sampling time section corresponding to the sliding scale of the sliding window, covering and time span are greater than benchmark recorder data sequence institute
The sampling time section answered；
The corresponding multiple related recorder data sequences of each control acquisition unit are calculated separately, with benchmark recorder data sequence
Related coefficient；
Based on the corresponding related coefficient calculated result of each control acquisition unit, calculates corresponding control acquisition unit and adopted with benchmark
Collect the time synchronization error between unit, and then obtains the time synchronization error between any two acquisition unit.
Specifically, assuming there is acquisition unit A, B, C, using acquisition unit A as benchmark acquisition unit, control acquisition is calculated separately
Time synchronization error between unit B and C and A is tAB、tAC, then the time synchronization error compareed between acquisition unit B and C can
It is calculated, is tBC=tAC-tABTo get the time synchronization error arrived between any two acquisition unit.And then it can determine whether event
Whether in allowed limits to hinder indicator three-phase synchronous error, such as the synchronization time of the three-phase line acquisition unit of conventional need
Error should be less than 100 μ s.
The present invention " obtains fault detector three-phase acquisition unit to respond the test current waveform on same route and collect
, the corresponding three groups of recorder datas of sampling instant " in, each phase acquisition unit that " acquisition moment corresponding " finally obtains
Sampling periods corresponding to recorder data are identical, and for each sampling instant in sampling periods, in each group recorder data point
It Ju You not a corresponding sampled point.
Preferably, defining data window width is N, and benchmark recorder data sequence is X=(x [1], x [2] ..., x [N]), the
The corresponding related recorder data sequence of k sliding window is Yk=(yk[1],yk[2],...,yk[N]), then related recorder data sequence
Arrange YkWith the correlation coefficient r between benchmark recorder data sequence XkAre as follows:
Wherein n is the variable of N.Related coefficient algorithm is the prior art.
Preferably, described to be based on the corresponding related coefficient calculated result of each control acquisition unit, calculate corresponding control acquisition
Time synchronization error between unit and benchmark acquisition unit includes:
Obtain related coefficient maximum value, the maximum related recorder data sequence of related coefficient and benchmark recorder data sequence it
Between time shift relationship and sliding scale in, the related recording number adjacent to the maximum related recorder data sequence of related coefficient
According to the corresponding related coefficient of sequence；
Based on the related coefficient and the time shift relationship got, each control acquisition unit and benchmark acquisition unit are calculated
Time synchronization error.The time shift relationship is the corresponding sampling time section of correlation recorder data sequence and benchmark recorder data sequence
Sampling instant deviation between corresponding sampling time section.That is the present invention uses interpolation on the basis of related coefficient algorithm
Approximating method achievees the purpose that approach real time synchronous error, so that the calculated result of time synchronization error essence with higher
Degree.
Preferably, for each control acquisition unit, the sampling instant deviation of definition and benchmark recorder data sequence for m,
Correlation coefficient r between the corresponding related recorder data sequence of m-th of sliding window and benchmark recorder data sequencemFor related coefficient
Maximum value then compares the time synchronization error t of acquisition unit Yu benchmark acquisition unit are as follows:
Wherein, Δ T is the sampling period of testing stand, rm+1And rm-1Respectively m+1 and the corresponding phase of the m-1 sliding window
Close the related coefficient between recorder data sequence and benchmark recorded broadcast data sequence.In the present invention, the corresponding correlation of m-th of sliding window
The sampling time section of recorder data sequence, the time shift relationship between sampling time section corresponding with benchmark recorded broadcast data sequence are adopted
Sample Time of day offsets is m.
Preferably, the data window width is 80 sampled points.Then for sample frequency 4000Hz, (sampling period is 250 μ
S) testing stand, the sampling period of single recorder data sequence are 20ms, for current most testing stand, be can reach higher
Testing efficiency.
Preferably, sampling instant corresponding to the sliding scale both ends recorder data of the sliding window, it is respectively in advance and stagnant
The sampling instant corresponding to benchmark recorder data sequence both ends afterwards, and the time span of lead and lag is equal.Such as surpass respectively
Preceding and J sampling instant of lag, J is positive integer, then for single related recorder data sequence, the related coefficient being calculated has
(r-J,r-J+1,r-J+2,...,r0,r1,...,rJ-1,rJ), when for compareing existing between acquisition unit and benchmark acquisition unit
Between synchronize may be advanced, it is also possible to the case where lagging, it is clear that the setting of sliding scale of the present invention all can adapt to, and guarantee phase relation
Several maximum values can be obtained from the related recorder data sequence chosen, and then guarantee the standard that time synchronization error calculates
True reliability.
Preferably, benchmark recorder data is compared in sampling instant corresponding to the sliding scale both ends recorder data of sliding window
Sampling instant corresponding to the recorder data of sequence both ends, difference 20 sampled points of lead and lag.
Invention additionally discloses a kind of distribution line failure indicator three-phase synchronous accuracy test devices, comprising: recorder data
Module is obtained, the test current waveform on same route is responded for obtaining fault detector three-phase acquisition unit and collects
, the corresponding three groups of recorder datas of sampling instant；
Benchmark recorder data sequence chooses module, for the data window width based on setting from preset benchmark acquisition unit
The recorder data that continuous sampling is chosen in corresponding recorder data, obtains benchmark recorder data sequence；
Related recorder data sequence chooses module, is being set for the data window width based on aforementioned setting using sliding window
In fixed sliding scale, sliding is chosen and is continuously adopted respectively from the corresponding recorder data of preset each control acquisition unit every time
The recorder data of sample obtains multiple related recorder data sequences；Sampling time section corresponding to the sliding scale of the sliding window,
It covers and time span is greater than sampling time section corresponding to benchmark recorder data sequence；
Related coefficient computing module, for calculating separately the corresponding multiple related recorder data sequences of each control acquisition unit
Column, the related coefficient with benchmark recorder data sequence；
Time synchronization error computing module, for being based on the corresponding related coefficient calculated result of each control acquisition unit, meter
Calculate the time synchronization error between corresponding control acquisition unit and benchmark acquisition unit, so obtain any two acquisition unit it
Between time synchronization error.
Preferably, defining data window width is N, and benchmark recorder data sequence is X=(x [1], x [2] ..., x [N]), the
The corresponding related recorder data sequence of k sliding window is Yk=(yk[1],yk[2],...,yk[N]), then related coefficient calculates mould
When block carries out related coefficient calculating, related recorder data sequence YkWith the correlation coefficient r between benchmark recorder data sequence XkAre as follows:
Wherein n is the variable of N.
Preferably, time synchronization error computing module calculates the time between each control acquisition unit and benchmark acquisition unit
When synchronous error:
According to related coefficient calculated result, related coefficient maximum value, the maximum related recorder data sequence of related coefficient are obtained
Arrange between benchmark recorder data sequence time shift relationship and sliding scale in, with the maximum related recording number of related coefficient
According to the adjacent corresponding related coefficient of related recorder data sequence of sequence；
For each control acquisition unit, it is assumed that the sampling instant deviation of benchmark recorder data sequence is m, m-th of sliding
Correlation coefficient r between the corresponding related recorder data sequence of window and benchmark recorder data sequencemFor related coefficient maximum value, then
Compare the time synchronization error t of acquisition unit and benchmark acquisition unit are as follows:
Wherein, Δ T is the sampling period of testing stand, rm+1And rm-1Respectively m+1 and the corresponding phase of the m-1 sliding window
Close the related coefficient between recorder data sequence and benchmark recorded broadcast data sequence.In the present invention, the corresponding correlation of m-th of sliding window
Sampling instant deviation between the sampling time section of recorder data sequence, and the corresponding sampling time section of benchmark recorded broadcast data sequence
As m.
Invention additionally discloses a kind of distribution line failure indicator three-phase synchronous accuracy test systems, including testing stand and event
Hinder indicator three-phase synchronous accuracy test device；
The input terminal that fault detector to be measured corresponds to three acquisition units of three-phase is separately connected the same phase line of testing stand
Road, testing stand primary side are passed through test current waveform；
Distribution line failure indicator three-phase synchronous accuracy test device executes aforementioned distribution line failure indicator three-phase
Synchronization accuracy test method, to obtain the time synchronization error between fault detector any two acquisition unit.
Test current waveform of the present invention is the existing current waveforms such as square wave or sine wave.
Beneficial effect
(1) the three-phase synchronous accuracy test that distribution line failure indicator is carried out using the present invention, without changing testing stand
Sample rate, therefore using the testing stand of traditional routine sampling rate, it can be to three-phase acquisition unit error synchronization time whether in institute
Judged in the range of permission, and judging result is effectively reliable, testing cost is lower.
(2) test method of the invention is benchmark acquisition unit by the way that any acquisition unit is arranged, and other acquisitions are single
Member is control acquisition unit, and then chooses benchmark recorder data sequence using fixed window, multiple related to being chosen using sliding window
Recorder data sequence recycles related coefficient algorithm, calculates separately multiple related recorder data sequences and benchmark recorder data sequence
Related coefficient between column, and according to related coefficient calculated result, based on multiple related coefficient the larger value and corresponding related record
Each control acquisition unit and benchmark is calculated in sampling time section deviation between wave data sequence and benchmark recorder data sequence
Time synchronization error between acquisition unit achievees the purpose that compare three-phase acquisition unit waveform time delay, and then finally obtains and appoint
The time synchronization error anticipated between two acquisition units.
(3) present invention calculates time synchronization error in conjunction with interpolation fitting method, improves on the basis of correlation coefficient process
The precision of Conventional temporal measurement, has ensured the reliability of time synchronization error result, while algorithm redundancy is low, easy to accomplish.
Detailed description of the invention
Fig. 1 show carry scheme of installation of the fault detector acquisition unit on testing stand；
Fig. 2 show test current waveform and fault detector three-phase acquisition unit recording waveform diagram；
Fig. 3 show fixed window and sliding window chooses recorder data sequence diagram；
Fig. 4 show the flow diagram that fault detector three-phase synchronous accuracy test is carried out using the present invention.
Specific embodiment
It is further described below in conjunction with the drawings and specific embodiments.
Embodiment 1
Refering to what is shown in Fig. 3, a kind of distribution line failure indicator three-phase synchronous method for testing precision of the present embodiment, packet
It includes:
It obtains fault detector three-phase acquisition unit and responds the test current waveform on same route and collected, sampling
Moment corresponding three groups of recorder datas；
Using any acquisition unit as benchmark acquisition unit, other acquisition units are as control acquisition unit, based on setting
Data window width the recorder data of continuous sampling is chosen from the corresponding recorder data of benchmark acquisition unit, obtain benchmark recording
Data sequence；
Data window width based on aforementioned setting, using sliding window, in the sliding scale of setting, every time sliding respectively from
The recorder data for choosing continuous sampling in the corresponding recorder data of acquisition unit is compareed, multiple related recorder data sequences are obtained；
It is right that sampling time section corresponding to the sliding scale of the sliding window, covering and time span are greater than benchmark recorder data sequence institute
The sampling time section answered；
The corresponding multiple related recorder data sequences of each control acquisition unit are calculated separately, with benchmark recorder data sequence
Related coefficient；
Based on the corresponding related coefficient calculated result of each control acquisition unit, calculates corresponding control acquisition unit and adopted with benchmark
Collect the time synchronization error between unit, and then obtains the time synchronization error between any two acquisition unit.
Definition data window width is N, and benchmark recorder data sequence is X=(x [1], x [2] ..., x [N]), k-th of sliding
The corresponding related recorder data sequence of window is Yk=(yk[1],yk[2],...,yk[N]), then related recorder data sequence YkWith base
Correlation coefficient r between quasi- recorder data sequence XkAre as follows:
Wherein n is the variable of N.Related coefficient algorithm is the prior art.
It is described to be based on the corresponding related coefficient calculated result of each control acquisition unit, calculate corresponding control acquisition unit and base
Time synchronization error between quasi- acquisition unit includes:
Obtain related coefficient maximum value, the maximum related recorder data sequence of related coefficient and benchmark recorder data sequence it
Between time shift relationship and sliding scale in, the related recording number adjacent to the maximum related recorder data sequence of related coefficient
According to the corresponding related coefficient of sequence；
Based on the related coefficient and the time shift relationship got, each control acquisition unit and benchmark acquisition unit are calculated
Time synchronization error.The time shift relationship is the corresponding sampling time section of correlation recorder data sequence and benchmark recorder data sequence
Sampling instant deviation between corresponding sampling time section.That is the present invention uses interpolation on the basis of related coefficient algorithm
Approximating method achievees the purpose that approach real time synchronous error, so that the calculated result of time synchronization error essence with higher
Degree.
For each control acquisition unit, the sampling instant deviation of definition and benchmark recorder data sequence is m, m-th of cunning
Move the correlation coefficient r between the corresponding related recorder data sequence of window and benchmark recorder data sequencemFor related coefficient maximum value,
Then compare the time synchronization error t of acquisition unit Yu benchmark acquisition unit are as follows:
Wherein, Δ T is the sampling period of testing stand, rm+1And rm-1Respectively m+1 and the corresponding phase of the m-1 sliding window
Close the related coefficient between recorder data sequence and benchmark recorded broadcast data sequence.In the present invention, the corresponding correlation of m-th of sliding window
The sampling time section of recorder data sequence, the time shift relationship between sampling time section corresponding with benchmark recorded broadcast data sequence are adopted
Sample Time of day offsets is m.
Utilizing the present embodiment, it is assumed that fault detector has acquisition unit A, B, C, using acquisition unit A as benchmark acquisition unit,
Calculating separately the time synchronization error between control acquisition unit B and C and A is tAB、tAC, then compare between acquisition unit B and C
Time synchronization error can be calculated, be tBC=tAC-tABTo get the time synchronization arrived between any two acquisition unit
Error.And then whether in allowed limits to can determine whether fault detector three-phase synchronous error, such as the three-phase line of conventional need
Error synchronization time of acquisition unit should be less than 100 μ s.
Embodiment 2
Based on the test method of embodiment 1, as shown in Figure 1, the three-phase acquisition unit of fault detector to be measured is mounted on
On the same phase line of testing stand, test current waveform then is passed through to testing stand primary side.
The three-phase acquisition unit of fault detector distinguishes response test current waveform, and collected current waveform can be from failure
The unit that collects of indicator obtains, and with reference to Fig. 2, the acquisition waveform of three acquisition units is due to various reality factors, unavoidably
There are time synchronization errors.
The three-phase acquisition unit of failure definition indicator is acquisition unit 1/2/3, and the present embodiment is obtaining three-phase acquisition list
It after the recorder data of member, chooses 1 benchmark acquisition unit of acquisition unit and refers to phase, using fixed window from benchmark acquisition unit pair
80 sampled points that continuous sampling is chosen in the benchmark recorder data answered, form benchmark recorder data sequence.Control is acquired
Sliding window is respectively adopted from corresponding recorder data in unit 2 and acquisition unit 3, repeatedly slides, and sliding chooses 80 respectively every time
A continuous sampling point forms multiple related recorder data series.Refering to what is shown in Fig. 3, t2Moment is the starting sample point of fixed window
Sampling instant, the starting sample point sampling moment of sliding window is in t1And t3It is slided between moment.
It is 80 sampled points that data window width, which is arranged, in the present embodiment, then for the sample frequency 4000Hz (sampling period 250
μ s) testing stand, sampling period of single recorder data sequence is 20ms, for current most testing stand, can reach compared with
High testing efficiency.
If time synchronization error is not present in three acquisition units, the recording waveform for corresponding to three acquisition units should not had
It is differentiated.But regardless of whether there are time synchronization error, the sampling instant of the recorder data of corresponding three acquisition units is all
It is mutual corresponding.
In the present embodiment, sampling instant corresponding to the sliding scale both ends recorder data of the sliding window is advanced respectively
With lag behind sampling instant corresponding to benchmark recorder data sequence both ends, and the time span of lead and lag is equal.As divided
Other J sampling instant of lead and lag, J is positive integer, then for single related recorder data sequence, the phase relation being calculated
Number has (r-J,r-J+1,r-J+2,...,r0,r1,...,rJ-1,rJ), for existing between control acquisition unit and benchmark acquisition unit
Time synchronization may be advanced, it is also possible to the case where lagging, it is clear that the setting of the present embodiment sliding scale all can adapt to, and guarantee
The maximum value of related coefficient can be obtained from the related recorder data sequence chosen, and then guarantee time synchronization error meter
The accurate reliability of calculation.
When carrying out related recorder data sequence selection using sliding window, the sliding scale both ends recorder data institute of sliding window is right
The sampling instant answered, compared to sampling instant corresponding to the recorder data of benchmark recorder data sequence both ends, difference lead and lag
20 sampled points.Refering to what is shown in Fig. 3, the sliding scale of i.e. sliding window is to shift to an earlier date 20 sampled points after 20 with reference to fixed window
Sampled point, thus amounts to 41 sliding windows, and corresponding each control acquisition unit forms 41 related recorder data sequences, according to each
The sampling time section of related recorder data sequence and the time shift of benchmark recorder data sequential sampling period are numbered, i.e.,
{Y-20,Y-19,...,Y-1,Y0,Y1,...,Y20, wherein related recorder data sequence Y0With benchmark recorder data sequence X=(x
[1], [2] x ..., x [N]) in each sampled point sampling instant correspond.Then for acquisition unit 2 and acquisition unit 3, divide
41 related coefficient { r can not be calculated-20,r-19,...,r-1,r0,r1,...,r20}。
Assuming that in 41 related coefficients, maximum value and the sampling instant deviation of benchmark recorder data sequence be m, m
I.e. Y at a sliding windowmIt obtains, is rm, then further combined with YmAdjacent related recorder data sequence Ym-1And Ym+1Corresponding correlation
Coefficient rm-1And rm+1, the time synchronization error between control acquisition unit and benchmark acquisition unit is carried out using difference approximating method
Calculating, formula are as follows:
Assuming that the time synchronization error being calculated between acquisition unit 2/3 and acquisition unit 1 is tAB/tAC, then list is acquired
Time synchronization error between member 2 and acquisition unit 3 is tBC=tAC-tABTo get arrived between any two acquisition unit when
Between synchronous error.And then whether in allowed limits to can determine whether fault detector three-phase synchronous error,
Interpolation fitting algorithm is the prior art, in order to further increase the accuracy after fitting, be can be selected more related
Coefficient the larger value participates in interpolation fitting.
Embodiment 3
It is based on identical inventive concept with Examples 1 and 2, the present embodiment is that a kind of distribution line failure indicator three is identical
Walk accuracy test device, comprising: recorder data obtains module, responds same line for obtaining fault detector three-phase acquisition unit
The test current waveform of road and collected, the corresponding three groups of recorder datas of sampling instant；
Benchmark recorder data sequence chooses module, for the data window width based on setting from preset benchmark acquisition unit
The recorder data that continuous sampling is chosen in corresponding recorder data, obtains benchmark recorder data sequence；
Related recorder data sequence chooses module, is being set for the data window width based on aforementioned setting using sliding window
In fixed sliding scale, sliding is chosen and is continuously adopted respectively from the corresponding recorder data of preset each control acquisition unit every time
The recorder data of sample obtains multiple related recorder data sequences；Sampling time section corresponding to the sliding scale of the sliding window,
It covers and time span is greater than sampling time section corresponding to benchmark recorder data sequence；
Related coefficient computing module, for calculating separately the corresponding multiple related recorder data sequences of each control acquisition unit
Column, the related coefficient with benchmark recorder data sequence；
Time synchronization error computing module, for being based on the corresponding related coefficient calculated result of each control acquisition unit, meter
Calculate the time synchronization error between corresponding control acquisition unit and benchmark acquisition unit, so obtain any two acquisition unit it
Between time synchronization error.
Definition data window width is N, and benchmark recorder data sequence is X=(x [1], x [2] ..., x [N]), k-th of sliding
The corresponding related recorder data sequence of window is Yk=(yk[1],yk[2],...,yk[N]), then related coefficient computing module carries out phase
When relationship number calculates, related recorder data sequence YkWith the correlation coefficient r between benchmark recorder data sequence XkAre as follows:
Wherein n is the variable of N.
The time synchronization that time synchronization error computing module calculates between each control acquisition unit and benchmark acquisition unit is missed
When poor:
According to related coefficient calculated result, related coefficient maximum value, the maximum related recorder data sequence of related coefficient are obtained
Arrange between benchmark recorder data sequence time shift relationship and sliding scale in, with the maximum related recording number of related coefficient
According to the adjacent corresponding related coefficient of related recorder data sequence of sequence；
For each control acquisition unit, it is assumed that the sampling instant deviation of benchmark recorder data sequence is m, m-th of sliding
Correlation coefficient r between the corresponding related recorder data sequence of window and benchmark recorder data sequencemFor related coefficient maximum value, then
Compare the time synchronization error t of acquisition unit and benchmark acquisition unit are as follows:
Wherein, Δ T is the sampling period of testing stand, rm+1And rm-1Respectively m+1 and the corresponding phase of the m-1 sliding window
Close the related coefficient between recorder data sequence and benchmark recorded broadcast data sequence.In the present invention, the corresponding correlation of m-th of sliding window
Sampling instant deviation between the sampling time section of recorder data sequence, and the corresponding sampling time section of benchmark recorded broadcast data sequence
As m.
Embodiment 4
It is based on identical inventive concept with embodiment 1 to 3, the present embodiment is that a kind of distribution line failure indicator three is identical
Walk accuracy test system, including testing stand and fault detector three-phase synchronous accuracy test device；
The input terminal that fault detector to be measured corresponds to three acquisition units of three-phase is separately connected the same phase line of testing stand
Road, testing stand primary side are passed through test current waveform；
Distribution line failure indicator three-phase synchronous accuracy test device executes distribution line failure described in embodiment 1 and refers to
Show device three-phase synchronous method for testing precision, to obtain the time synchronization error between fault detector any two acquisition unit.
Test current waveform of the present invention is the existing current waveforms such as square wave or sine wave.
To sum up, the three-phase synchronous accuracy test that distribution line failure indicator is carried out using the present invention, without changing test
Platform sample rate, therefore using the testing stand of traditional routine sampling rate, can to three-phase acquisition unit error synchronization time whether
Judged in permitted range, and judging result is effectively reliable, testing cost is lower.
The present invention utilizes correlation coefficient process simultaneously, achievees the purpose that compare three-phase acquisition unit waveform time delay, and in correlation
On the basis of Y-factor method Y, time synchronization error is calculated in conjunction with interpolation fitting method, improves the precision of Conventional temporal measurement, is ensured
The reliability of time synchronization error result, while algorithm redundancy is low, it is easy to accomplish.
It should be understood by those skilled in the art that, embodiments herein can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the application, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
The embodiment of the present invention is described in conjunction with attached drawing above, but the invention is not limited to above-mentioned specific
Embodiment, the above mentioned embodiment is only schematical, rather than restrictive, those skilled in the art
Under the inspiration of the present invention, without breaking away from the scope protected by the purposes and claims of the present invention, it can also make very much
Form, all of these belong to the protection of the present invention.
Claims (11)
1. a kind of distribution line failure indicator three-phase synchronous method for testing precision, characterized in that include:
It obtains fault detector three-phase acquisition unit and responds the test current waveform on same route and collected, sampling instant
Corresponding three groups of recorder datas；
Using any acquisition unit as benchmark acquisition unit, other acquisition units are as control acquisition unit, the number based on setting
The recorder data for choosing continuous sampling from the corresponding recorder data of benchmark acquisition unit according to window width, obtains benchmark recorder data
Sequence；
Data window width based on aforementioned setting, using sliding window, in the sliding scale of setting, sliding is respectively from control every time
The recorder data that continuous sampling is chosen in the corresponding recorder data of acquisition unit obtains multiple related recorder data sequences；It is described
Sampling time section corresponding to the sliding scale of sliding window, covering and time span are greater than corresponding to benchmark recorder data sequence
Sampling time section；
The corresponding multiple related recorder data sequences of each control acquisition unit are calculated separately, it is related to benchmark recorder data sequence
Coefficient；
Based on the corresponding related coefficient calculated result of each control acquisition unit, calculates corresponding control acquisition unit and benchmark acquisition is single
Time synchronization error between member, and then obtain the time synchronization error between any two acquisition unit.
2. according to the method described in claim 1, it is characterized in that, definition data window width be N, benchmark recorder data sequence be X
=(x [1], x [2] ..., x [N]), the corresponding related recorder data sequence of k-th of sliding window is Yk=(yk[1],yk
[2],...,yk[N]), then related recorder data sequence YkWith the correlation coefficient r between benchmark recorder data sequence XkAre as follows:
Wherein n is the variable of N.
3. according to the method described in claim 1, it is characterized in that, it is described based on the corresponding related coefficient of each control acquisition unit
It calculates as a result, the time synchronization error calculated between corresponding control acquisition unit and benchmark acquisition unit includes:
Related coefficient maximum value is obtained, between the maximum related recorder data sequence of related coefficient and benchmark recorder data sequence
In time shift relationship and sliding scale, the related recorder data sequence adjacent to the maximum related recorder data sequence of related coefficient
Arrange corresponding related coefficient；
Based on the related coefficient and the time shift relationship got, the time of each control acquisition unit and benchmark acquisition unit is calculated
Synchronous error.
4. according to the method described in claim 3, it is characterized in that, for each control acquisition unit, definition and benchmark recorder data
The sampling instant deviation of sequence be m, the corresponding related recorder data sequence of m-th of sliding window and benchmark recorder data sequence it
Between correlation coefficient rmFor related coefficient maximum value, then the time synchronization error t of acquisition unit Yu benchmark acquisition unit is compareed are as follows:
Wherein, Δ T is the sampling period of testing stand, rm+1And rm-1Respectively m+1 and the corresponding related record of the m-1 sliding window
Related coefficient between wave data sequence and benchmark recorded broadcast data sequence.
5. according to the method described in claim 1, it is characterized in that, the data window width be 80 sampled points.
6. according to the method described in claim 1, it is characterized in that, corresponding to the sliding scale both ends recorder data of the sliding window
Sampling instant, lead and lag sampling instant corresponding to benchmark recorder data sequence both ends respectively, and lead and lag
Time span it is equal.
7. according to the method described in claim 6, it is characterized in that, adopted corresponding to the sliding scale both ends recorder data of sliding window
The sample moment adopts for lead and lag 20 respectively compared to sampling instant corresponding to the recorder data of benchmark recorder data sequence both ends
Sampling point.
8. a kind of distribution line failure indicator three-phase synchronous accuracy test device, characterized in that include:
Recorder data obtains module, responds test current wave on same route for obtaining fault detector three-phase acquisition unit
Shape and collected, the corresponding three groups of recorder datas of sampling instant；
Benchmark recorder data sequence chooses module, corresponding from preset benchmark acquisition unit for the data window width based on setting
Recorder data in choose continuous sampling recorder data, obtain benchmark recorder data sequence；
Related recorder data sequence chooses module, for the data window width based on aforementioned setting, using sliding window, in setting
In sliding scale, sliding chooses continuous sampling respectively from the corresponding recorder data of preset each control acquisition unit every time
Recorder data obtains multiple related recorder data sequences；Sampling time section corresponding to the sliding scale of the sliding window, covering
And time span is greater than sampling time section corresponding to benchmark recorder data sequence；
Related coefficient computing module, for calculating separately the corresponding multiple related recorder data sequences of each control acquisition unit, with
The related coefficient of benchmark recorder data sequence；
Time synchronization error computing module calculates phase for being based on the corresponding related coefficient calculated result of each control acquisition unit
Reply obtains between any two acquisition unit according to the time synchronization error between acquisition unit and benchmark acquisition unit
Time synchronization error.
9. distribution line failure indicator three-phase synchronous accuracy test device according to claim 8, characterized in that definition
Data window width is N, and benchmark recorder data sequence is X=(x [1], x [2] ..., x [N]), the corresponding correlation of k-th of sliding window
Recorder data sequence is Yk=(yk[1],yk[2],...,yk[N]), then related coefficient computing module carries out related coefficient calculating
When, related recorder data sequence YkWith the correlation coefficient r between benchmark recorder data sequence XkAre as follows:
Wherein n is the variable of N.
10. distribution line failure indicator three-phase synchronous accuracy test device according to claim 8, characterized in that when
Between synchronous error computing module calculate it is each control acquisition unit and benchmark acquisition unit between time synchronization error when:
According to related coefficient calculated result, obtain related coefficient maximum value, the maximum related recorder data sequence of related coefficient with
In time shift relationship and sliding scale between benchmark recorder data sequence, to the maximum related recorder data sequence of related coefficient
Arrange the adjacent corresponding related coefficient of related recorder data sequence；
For each control acquisition unit, it is assumed that the sampling instant deviation of benchmark recorder data sequence is m, m-th of sliding window pair
Correlation coefficient r between the related recorder data sequence answered and benchmark recorder data sequencemFor related coefficient maximum value, then compare
The time synchronization error t of acquisition unit and benchmark acquisition unit are as follows:
Wherein, Δ T is the sampling period of testing stand, rm+1And rm-1Respectively m+1 and the corresponding related record of the m-1 sliding window
Related coefficient between wave data sequence and benchmark recorded broadcast data sequence.In the present invention, the corresponding related recording of m-th of sliding window
Sampling instant deviation between the sampling time section of data sequence, and the corresponding sampling time section of benchmark recorded broadcast data sequence is
m。
11. a kind of based on any one of the claim 1-7 distribution line failure indicator three-phase synchronous method for testing precision
Distribution line failure indicator three-phase synchronous accuracy test system, characterized in that identical including testing stand and fault detector three
Walk accuracy test device；
The input terminal that fault detector to be measured corresponds to three acquisition units of three-phase is separately connected the same phase line of testing stand, examination
It tests platform primary side and is passed through test current waveform；
Distribution line failure indicator three-phase synchronous accuracy test device perform claim requires the described in any item distribution lines of 1-7
Fault detector three-phase synchronous method for testing precision, to obtain the time synchronization between fault detector any two acquisition unit
Error.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811051066.8A CN109298362B (en) | 2018-09-10 | 2018-09-10 | Three-phase synchronization precision testing method, device and system for distribution line fault indicator |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811051066.8A CN109298362B (en) | 2018-09-10 | 2018-09-10 | Three-phase synchronization precision testing method, device and system for distribution line fault indicator |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109298362A true CN109298362A (en) | 2019-02-01 |
| CN109298362B CN109298362B (en) | 2020-10-16 |
Family
ID=65166692
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811051066.8A Active CN109298362B (en) | 2018-09-10 | 2018-09-10 | Three-phase synchronization precision testing method, device and system for distribution line fault indicator |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109298362B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111896788A (en) * | 2020-08-07 | 2020-11-06 | 西安远测电力科技有限公司 | Trigger method for current recording of opening and closing coil of circuit breaker |
| CN113466615A (en) * | 2021-06-17 | 2021-10-01 | 三峡大学 | Drosophila optimization algorithm-based post-fault wave recording data synchronization method |
| CN114325555A (en) * | 2022-01-04 | 2022-04-12 | 国网上海市电力公司 | Metering equipment abnormity on-line monitoring model and error calculation method |
| CN115327371A (en) * | 2022-10-17 | 2022-11-11 | 广东电网有限责任公司佛山供电局 | Method for diagnosing state time point of medium-voltage switch and related device |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202256653U (en) * | 2011-10-12 | 2012-05-30 | 华北电力大学 | Overhead line type fault indicator performance test device |
| CN103698584A (en) * | 2014-01-08 | 2014-04-02 | 广州思泰信息技术有限公司 | Power distribution network overhead line three-phase current sampling synchronization method |
| CN104467186A (en) * | 2014-12-15 | 2015-03-25 | 国网安徽省电力公司 | Full network fault recording data parameter-free accurate synchronizing method |
| CN105721219A (en) * | 2016-03-07 | 2016-06-29 | 宁波迦南电子有限公司 | Power line fault indicator with time synchronization function |
| CN107015110A (en) * | 2017-04-01 | 2017-08-04 | 云南电网有限责任公司瑞丽供电局 | Synchronous recording fault location system, method |
| CN206515424U (en) * | 2017-01-22 | 2017-09-22 | 营口电力勘察设计院 | Novel electric power fault location system |
| CN107272525A (en) * | 2017-08-16 | 2017-10-20 | 杭州大有科技发展有限公司 | A kind of fault detector |
| CN206804790U (en) * | 2017-04-20 | 2017-12-26 | 广西电网有限责任公司电力科学研究院 | A kind of intelligent distribution network line fault monitoring system |
| US20180062553A1 (en) * | 2016-08-31 | 2018-03-01 | Intel Corporation | Monitoring health of electrical equipment |
| CN107959969A (en) * | 2017-12-11 | 2018-04-24 | 武汉中原电子信息有限公司 | A kind of method for synchronizing time applied to transient state recording type fault detector |
| CN108051708A (en) * | 2017-12-31 | 2018-05-18 | 杭州柯林电气股份有限公司 | A kind of transient state recording type fault detector |
| CN108107321A (en) * | 2017-12-14 | 2018-06-01 | 科大智能电气技术有限公司 | A kind of electric power system fault waveform comparison method |
| CN108226836A (en) * | 2016-12-31 | 2018-06-29 | 科大智能电气技术有限公司 | A kind of calibration method of fault detector synchronous waveform |
| CN108445279A (en) * | 2018-02-12 | 2018-08-24 | 北京和信瑞通电力技术股份有限公司 | Based on the wireless synchronous sampling method of short frequency in transient state recording type fault detector |
| CN108490310A (en) * | 2018-03-21 | 2018-09-04 | 国网河南省电力公司漯河供电公司 | Method for small electric current grounding system of distribution network earth fault line selection |
-
2018
- 2018-09-10 CN CN201811051066.8A patent/CN109298362B/en active Active
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202256653U (en) * | 2011-10-12 | 2012-05-30 | 华北电力大学 | Overhead line type fault indicator performance test device |
| CN103698584A (en) * | 2014-01-08 | 2014-04-02 | 广州思泰信息技术有限公司 | Power distribution network overhead line three-phase current sampling synchronization method |
| CN104467186A (en) * | 2014-12-15 | 2015-03-25 | 国网安徽省电力公司 | Full network fault recording data parameter-free accurate synchronizing method |
| CN105721219A (en) * | 2016-03-07 | 2016-06-29 | 宁波迦南电子有限公司 | Power line fault indicator with time synchronization function |
| US20180062553A1 (en) * | 2016-08-31 | 2018-03-01 | Intel Corporation | Monitoring health of electrical equipment |
| CN108226836A (en) * | 2016-12-31 | 2018-06-29 | 科大智能电气技术有限公司 | A kind of calibration method of fault detector synchronous waveform |
| CN206515424U (en) * | 2017-01-22 | 2017-09-22 | 营口电力勘察设计院 | Novel electric power fault location system |
| CN107015110A (en) * | 2017-04-01 | 2017-08-04 | 云南电网有限责任公司瑞丽供电局 | Synchronous recording fault location system, method |
| CN206804790U (en) * | 2017-04-20 | 2017-12-26 | 广西电网有限责任公司电力科学研究院 | A kind of intelligent distribution network line fault monitoring system |
| CN107272525A (en) * | 2017-08-16 | 2017-10-20 | 杭州大有科技发展有限公司 | A kind of fault detector |
| CN107959969A (en) * | 2017-12-11 | 2018-04-24 | 武汉中原电子信息有限公司 | A kind of method for synchronizing time applied to transient state recording type fault detector |
| CN108107321A (en) * | 2017-12-14 | 2018-06-01 | 科大智能电气技术有限公司 | A kind of electric power system fault waveform comparison method |
| CN108051708A (en) * | 2017-12-31 | 2018-05-18 | 杭州柯林电气股份有限公司 | A kind of transient state recording type fault detector |
| CN108445279A (en) * | 2018-02-12 | 2018-08-24 | 北京和信瑞通电力技术股份有限公司 | Based on the wireless synchronous sampling method of short frequency in transient state recording type fault detector |
| CN108490310A (en) * | 2018-03-21 | 2018-09-04 | 国网河南省电力公司漯河供电公司 | Method for small electric current grounding system of distribution network earth fault line selection |
Non-Patent Citations (5)
| Title |
|---|
| CHENGXI ZHANG: "Fault-tolerant adaptive finite-time attitude synchronization and tracking control for multi-spacecraft formation", 《AEROSPACE SCIENCE AND TECHNOLOGY》 * |
| 李振龙: "故障指示器信号录波同步和检测算法研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
| 杨增健: "电网故障脉冲到达时间精确测量方法研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
| 燕跃豪: "基于IEC61850的数字化故障录波器的建模与实现", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
| 谢成: "基于相关系数分析的配电网单相接地故障研判方法与试验研究", 《浙江电力》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111896788A (en) * | 2020-08-07 | 2020-11-06 | 西安远测电力科技有限公司 | Trigger method for current recording of opening and closing coil of circuit breaker |
| CN111896788B (en) * | 2020-08-07 | 2023-02-07 | 西安远测电力科技有限公司 | Trigger method for current recording of opening and closing coil of circuit breaker |
| CN113466615A (en) * | 2021-06-17 | 2021-10-01 | 三峡大学 | Drosophila optimization algorithm-based post-fault wave recording data synchronization method |
| CN113466615B (en) * | 2021-06-17 | 2022-06-03 | 三峡大学 | Drosophila optimization algorithm-based post-fault wave recording data synchronization method |
| CN114325555A (en) * | 2022-01-04 | 2022-04-12 | 国网上海市电力公司 | Metering equipment abnormity on-line monitoring model and error calculation method |
| CN114325555B (en) * | 2022-01-04 | 2023-11-14 | 国网上海市电力公司 | Metering equipment abnormality online monitoring model and error calculation method |
| CN115327371A (en) * | 2022-10-17 | 2022-11-11 | 广东电网有限责任公司佛山供电局 | Method for diagnosing state time point of medium-voltage switch and related device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109298362B (en) | 2020-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109298362A (en) | A kind of distribution line failure indicator three-phase synchronous method for testing precision, device and system | |
| CN106646139B (en) | A kind of electrical power distribution network fault location method based on three-phase current amplitude analysis | |
| Bíscaro et al. | Integrated fault location and power-quality analysis in electric power distribution systems | |
| WO2019232595A1 (en) | A method of estimating the location of a fault on an electrical distribution network and an associated system | |
| RU2632989C2 (en) | Method and device for determining location of single-phase-to-ground fault in distributing network based on wavelet transformation of transitional signals | |
| JP3352411B2 (en) | Control system, power system protection control system, and storage medium storing program | |
| DK1922550T3 (en) | Method and system for monitoring the flow of energy through an electric energitransmissonsledning | |
| CN107925240B (en) | Test system for traveling wave fault detector | |
| CN105116246B (en) | A kind of operation of power networks protection on-line analysis method | |
| CN102879639A (en) | Real-time frequency measuring method in power system | |
| CN104597320B (en) | A kind of method calculated suitable for multiple frequency alternating current signals | |
| CN107767722B (en) | Multifunctional digital electric energy metering training simulation system | |
| CN106771861A (en) | Complex electric network Fault Locating Method based on wide area traveling wave energy and time difference | |
| CN106324340B (en) | A kind of method of synchronized phasor and frequency measurement dynamic property | |
| US20110043045A1 (en) | Systems and methods for asynchronous sampling data conversion | |
| CN104215838A (en) | Remote nuclear phase method for intelligent substation | |
| CN104820133A (en) | Online measurement method for impedance parameter of three-phase asymmetric power transmission line | |
| CN102183709B (en) | Method of determining fault point of power grid and severity of fault | |
| US6924647B2 (en) | Fault location method and device | |
| Abdolkhalig et al. | Phasor measurement based on IEC 61850-9-2 and Kalman–Filtering | |
| KR101664010B1 (en) | An Estimation Method of Line Parameter based on Synchrophasor Measurements in Power System | |
| CN107179476B (en) | Distribution network fault distance measurement method | |
| CN105866576A (en) | Simulation detection system of intelligent transformer station secondary-side electric energy metering error influence, and detection analysis method | |
| CN106646138B (en) | Distribution net work earthing fault localization method based on the conversion of more sample frequency wavelet character energy | |
| Ren | SYNCHROPHASOR MEASUREMENT USING SUBSTATION INTELLIGENT |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |