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 PDF

Info

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
Application number
CN201811051066.8A
Other languages
Chinese (zh)
Other versions
CN109298362B (en
Inventor
曾飞
袁栋
黄哲忱
袁宇波
周建华
孙健
高磊
卜强生
杨雄
彭志强
李虎成
吴楠
方鑫
张刘冬
陈实
李娟�
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd, Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201811051066.8A priority Critical patent/CN109298362B/en
Publication of CN109298362A publication Critical patent/CN109298362A/en
Application granted granted Critical
Publication of CN109298362B publication Critical patent/CN109298362B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00Testing 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

A kind of distribution line failure indicator three-phase synchronous method for testing precision, device and System
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.
CN201811051066.8A 2018-09-10 2018-09-10 Three-phase synchronization precision testing method, device and system for distribution line fault indicator Active CN109298362B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (15)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
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
CN105116246B (en) A kind of operation of power networks protection on-line analysis method
CN102879639A (en) Real-time frequency measuring method in power system
CN102377180A (en) Power system load modeling method based on electric energy quality monitoring system
CN104597320B (en) A kind of method calculated suitable for multiple frequency alternating current signals
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
CN109324228A (en) It is a kind of based on multi-frequency-Taylor's dynamic model synchronous phasor measuring method
CN103257273A (en) Method for measuring common-frequency periodic signal phase differences
CN104391177B (en) EMUs side harmonics test system and method
CN104215838A (en) Remote nuclear phase method for intelligent substation
CN104820133A (en) Online measurement method for impedance parameter of three-phase asymmetric power transmission line
CN102778634B (en) Distribution network voltage fluctuation positioning method and device
US6924647B2 (en) Fault location method and device
CN110297208A (en) Digitize live various dimensions test method and test macro across interval metering system
Abdolkhalig et al. Phasor measurement based on IEC 61850-9-2 and Kalman–Filtering
CN105866576A (en) Simulation detection system of intelligent transformer station secondary-side electric energy metering error influence, and detection analysis method
CN107179476B (en) Distribution network fault distance measurement method
CN109359822A (en) Electronic type voltage transformer measuring state appraisal procedure and system
CN103605904A (en) Self-compensating amplitude calculating method applied to electrical power system and based on error estimation

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