Summary of the invention
For when dependence impedance method simple in solution electrical power system transmission circuit is found range, distance accuracy can not meet range finding requirement, be subject to the factor impacts such as transition resistance, the easy error starting of simple dependence Fault Location With Traveling Wave timer, and during single end distance measurement, there is perception reflex ripple difficulty, exist during both-end distance measuring unstable networks, on time the factor such as dropout impact, propose a kind of hybrid ranging method that impedance method and traveling wave method combine, both avoid the error starting phenomenon existed in simple Fault Location With Traveling Wave, turn improve the reliability of measuring distance of transmission line fault.
For achieving the above object, concrete scheme of the present invention is as follows:
The row ripple hybrid ranging method that impedance method and traveling wave method combine, comprises the following steps:
The power-frequency voltage at power frequency data acquisition unit Real-time Collection circuit two ends and power current data, the high-frequency current channel data at high-frequency data collecting unit Real-time Collection circuit two ends, judges whether to break down according to the power-frequency voltage of Real-Time Monitoring and power current data;
Power frequency startup amount carries out power-frequency voltage, the operation of power current data record ripple at once after judging fault, carry out the operation of high-frequency traveling wave current record ripple simultaneously, obtain power frequency recorded wave file and the high-frequency current recorded wave file of the correspondence at circuit two ends after record ripple completes;
Judge mutually by starting the power frequency recorder data of record, the calculating of voltage, current effective value, utilizes impedance method analysis to obtain faulty line corresponding to circuit two ends number, fault moment, trouble spot respectively apart from the both-end resistance algorithm length of circuit and trouble spot apart from the single-ended impedance range finding length of circuit;
The faulty line number obtained with impedance method, fault moment are for reference, traveling wave method analysis is carried out to the high-frequency traveling wave current data corresponding respectively with power frequency recorder data of record, obtain both-end distance measuring length and the single end distance measurement length of circuit one end, and the both-end distance measuring length of the circuit other end and single end distance measurement length;
Relation between the range finding length at the range finding length at the circuit two ends obtained utilizing impedance method and the circuit two ends that utilize traveling wave method to obtain and line length provides the row ripple hybrid ranging length of trouble spot apart from circuit one end.
Further, the high frequency sample frequency that power frequency sample frequency is 5KHz, high-frequency data collecting unit is corresponding that power frequency data acquisition unit is corresponding is 2MHz; Power frequency judges that mode that record ripple starts comprises that voltage and current passage is out-of-limit, sudden change starts, and zero sequence and negative-sequence current sudden change, out-of-limit startup, zero sequence and negative sequence voltage sudden change, out-of-limit startup and switching value start.
Further, power frequency recorder data and high-frequency current recorder data are all with unified standard, and be the absolute timestamp information of GPS or the Big Dipper, absolute time when power frequency recorded wave file title and high frequency recorded wave file title are all to start is as file name.
Further, its record wavelength degree of power frequency recorded wave file and high frequency recorded wave file is configurable.
Further, faulty line number judges that the fault phase starting the operation of record ripple obtains according to power frequency.
Further, fault moment obtains after carrying out the detection of sudden change value to the voltage and current data in the power frequency recorder data obtained.
Further, both-end resistance algorithm length obtains by writing the differential equation to the voltage and current amount row in the power frequency recorder data file of both sides; Single-ended impedance range finding length is then by calculating the voltage and current amount in the power frequency recorded wave file of both sides respectively, obtains the resistance value of circuit when fault occurs and obtains compared with the line impedance of unit length.
The concrete grammar of Fault Location With Traveling Wave is:
A. establish and intercept time window t
m: with the fault moment t provided
0mfor benchmark, frontly push away Ams, pusher Bms, establish t
m∈ [t
0m-A, t
0m+ B] absolute time section as intercepting time window, wherein A and B is adjustable parameters, and is arithmetic number;
B. Mintrop wave, reflection wave absolute time is established: to high frequency recorded wave file y
min be in and intercept time window t
minterior high-frequency current data carry out phase-model transformation, wavelet transformation, obtain wavelet modulus maxima, and carry out amplitude screening to the modulus maximum that wavelet transformation obtains, and determine that the Mintrop wave head of initial row ripple arrives the absolute time t of local terminal
1the absolute time t of local terminal is arrived with trouble spot reflection wave wave head
2; In like manner to high frequency recorded wave file y
nin be in intercepting time t
nin data process equally, obtain this fault produce initial row ripple Mintrop wave head arrive N end absolute time t
3and trouble spot reflection wave wave head arrives the absolute time t of N end
4.
C. the length of the both-end distance measuring length of trouble spot distance M, N end and the Single Terminal Traveling Wave Fault Location of trouble spot distance M, N end is calculated according to absolute time difference.
Further, in described step C, the calculating of length is specially:
According to described absolute time t
1and t
3mistiming calculate trouble spot distance M hold both-end distance measuring length be:
Wherein L
mnfor circuit MN total length, v is the speed of row wave traveling;
According to described absolute time t
1and t
3mistiming calculate trouble spot distance N hold both-end distance measuring length be:
According to described absolute time t
1and t
2mistiming calculate trouble spot and apart from the length of Single Terminal Traveling Wave Fault Location of M end be:
According to described absolute time t
3and t
4mistiming calculate trouble spot and apart from the length of Single Terminal Traveling Wave Fault Location of N end be:
lmn is the distance between m and the n of circuit two ends, and v is traveling wave speed.
Further, power frequency recorded wave file x
mwith corresponding high frequency recorded wave file y
mbetween mapping relations establish in the following way:
First with power frequency recorded wave file time and x
mfor benchmark, frontly push away setting-up time T, pusher setting-up time T, form an absolute time sweep interval [x
m-T, x
m+ T];
Then scanning search is positioned at all high frequency recorded wave files in this interval, if there is multiple high speed recorded wave file, then gets and x
mthe minimum high frequency recorded wave file of the absolute value of mistiming as y
m.
Further, range finding length and the range finding length at the circuit two ends utilizing traveling wave method to obtain at the circuit two ends obtained utilizing impedance method carry out comprehensively analyzing the row ripple hybrid ranging length providing trouble spot and hold apart from M, and concrete steps are as follows:
1) if both-end resistance algorithm length, both-end travelling wave ranging length all exist, and the length all L within line length that find range
mdz∈ [0, L
mn], L
mdx∈ [0, L
mn]:
Calculate the absolute value L of the difference of both-end resistance algorithm length and both-end travelling wave ranging length
mdfor:
L
md=|L
mdz-L
mdx|
If 0≤L
md≤ (L
mn× i%): then row ripple hybrid ranging length A
mget row ripple both-end distance measuring length L
mdx, otherwise row ripple hybrid ranging length A
mget impedance both-end distance measuring length L
mdz; I is adjustable parameters;
2) if both-end resistance algorithm length L
mdz∈ [0, L
mn], both-end travelling wave ranging length L
mdxfor outside district, then row ripple hybrid ranging length A
mget both-end resistance algorithm length L
mdz;
3) if both-end travelling wave ranging length L
mdx∈ [0, L
mn], both-end resistance algorithm length L
mdzfor outside district, then row ripple hybrid ranging length A
mget both-end travelling wave ranging length L
mdx;
4) if both-end resistance algorithm length L
mdzfor outside district, both-end travelling wave ranging length L
mdxfor outside district, and single-ended impedance range finding length L
mszmeet L
msz∈ [0, L
mn], Single Terminal Traveling Wave Fault Location length L
msxmeet L
msx∈ [0, L
mn], then carry out following judgement:
If L
msz≤ (L
mn× j%), then row ripple hybrid ranging length A
mget L
msz, otherwise row ripple hybrid ranging length A
mget L
msx; Wherein, j is adjustable parameters;
5) if both-end resistance algorithm length L
mdzfor outside district, both-end travelling wave ranging length L
mdxfor outside district, and single-ended impedance range finding length L
mszmeet L
msz∈ [0, L
mn], Single Terminal Traveling Wave Fault Location length L
msxfor outside district, then row ripple hybrid ranging length A
mget single-ended impedance range finding length L
msz;
6) if both-end resistance algorithm length L
mdzfor outside district, both-end travelling wave ranging length L
mdxfor outside district, and Single Terminal Traveling Wave Fault Location length L
msxmeet L
msx∈ [0, L
mn], single-ended impedance range finding length L
mszfor outside district, then row ripple hybrid ranging length A
mget Single Terminal Traveling Wave Fault Location length L
msx;
7) if both-end resistance algorithm length L
mdzfor outside district, both-end travelling wave ranging length L
mdxfor outside district, Single Terminal Traveling Wave Fault Location length L
msxfor outside district, single-ended impedance range finding length L
mszfor outside district, then row ripple hybrid ranging length A
mfor outside district.
Beneficial effect of the present invention:
1, power frequency sampled data judges the mode that record ripple starts, and starts definite value and quantizes, effectively can reduce the error starting of circuit.
2, power frequency recorder data and the comparative analysis of high frequency recorder data, solve in simple Fault Location With Traveling Wave unanalyzable problem of can only finding range.
3, various location algorithm combines and ensure that the reliability of range finding, decreases the situation that range finding is failed.
Embodiment
Below in conjunction with accompanying drawing, embodiments of the present invention are illustrated:
Fig. 1 is specific embodiment of the invention process flow diagram, and its specific implementation process is as follows:
Carry out record ripple by the power-frequency voltage and power current of monitoring transmission line of electricity and start judgement.First, after transmission line of electricity breaks down, power frequency startup amount can trigger the record ripple record operation that record ripple logging software module carries out power-frequency voltage, power current data and high-frequency current data, obtains power frequency recorded wave file and high frequency recorded wave file;
Then impedance method is by calculating resistance algorithm length to the voltage and current amount in power frequency recorded wave file, the faulty line that traveling wave method provides with impedance method, fault moment are reference, the interval intercept operation of record ripple is carried out to high frequency recorder data, carries out the process such as wavelet transformation, phase-model transformation afterwards and obtain Fault Location With Traveling Wave length;
Row ripple hybrid ranging length is obtained finally by the comprehensive analysis of resistance algorithm length and travelling wave ranging length.
The row ripple hybrid ranging method that a kind of impedance method of the present invention and traveling wave method are combined being described as specifically:, comprise the following steps:
Step 1: the power-frequency voltage at power frequency data acquisition unit Real-time Collection circuit MN two ends and power current data also judge whether startup separator record ripple, the high-frequency current channel data at high-frequency data collecting unit Real-time Collection circuit MN two ends;
Step 2: power frequency startup amount at once starts power frequency logging software module and carries out power-frequency voltage, the operation of power current data record ripple after judging fault, trigger higher frequency logging software module carries out the operation of high-frequency traveling wave current record ripple simultaneously, and after record ripple completes, MN both sides obtain power frequency recorded wave file x
m, x
nand high-frequency current recorded wave file y
m, y
n;
Step 3: by the power frequency recorder data x to record in step 2
m, x
ncarry out starting judging mutually, voltage, current effective value calculate, and obtain M and hold faulty line m, M to hold fault moment t
0m, trouble spot apart from M end both-end resistance algorithm result L
mdzand trouble spot is apart from the single-ended impedance range measurement L of M end
msz; In like manner, faulty line n, the N that can obtain N end hold fault moment t
0n, trouble spot apart from N end both-end resistance algorithm result L
ndzand trouble spot is apart from the single-ended impedance range measurement L of N end
nsz;
Step 4: with the faulty line provided in step 3 number, fault moment for reference, to record in step 2 and x
mand x
nhigh-frequency traveling wave current data y corresponding respectively
mand y
ncarry out traveling wave method analysis, obtain M and hold traveling wave method both-end distance measuring result L
mdxwith single end distance measurement result L
msx, and N holds the both-end distance measuring result L of traveling wave method
ndxwith single end distance measurement result L
nsx;
Step 5: M in step 3 is held to the impedance method range measurement L obtained
mdzand L
mszfault Location With Traveling Wave result L is held with the M obtained in step 4
mdxand L
msxcarry out comprehensive analysis and provide the row ripple hybrid ranging result A of trouble spot apart from M end
m, in like manner obtain the row ripple hybrid ranging result A of trouble spot apart from N end
n.
More specifically, the power frequency sample frequency in described step 1 is 5KHz, high frequency sample frequency is 2MHz; Power frequency judges that mode that record ripple starts comprises that voltage and current passage is out-of-limit, sudden change starts, zero sequence and negative-sequence current sudden change, out-of-limit startup, and zero sequence and negative sequence voltage sudden change, out-of-limit startup and switching value start;
Power frequency recorder data in described step 2 and high-frequency current recorder data are all with the absolute timestamp information of unified standard (GPS or the Big Dipper), absolute time (during * * * * month * * day * *, * * divides * * * second * millisecond * * microsecond * * nanosecond) when power frequency recorded wave file title and high frequency recorded wave file title are all to start is as file name, as file is called the power frequency recorded wave file x of 2014-11-05-16-33-48-511-000,33 points of 48 seconds 511 milliseconds of 000 microseconds when the absolute time representing when power frequency amount judges fault initiating is 5 days 16 November in 2014; In like manner, file is called the high frequency recorded wave file y of 2014-11-05-16-33-48-518-170-000,33 points of 48 seconds 518 milliseconds of 170 microsecond 000 nanoseconds when absolute time when representing high frequency startup is 5 days 16 November in 2014.
Power frequency recorded wave file x in step 2
m, x
nand high frequency recorded wave file y
m, y
nits record wavelength degree is configurable;
According to power frequency, faulty line m and n in step 3 judges that the fault phase starting the operation of record ripple obtains;
Fault moment (t in step 3
0m, t
0n) be obtain after the detection of sudden change value is carried out to the voltage and current data in the power frequency recorder data obtained, if namely there are continuous print 3 sampled points, the effective value parameter of its corresponding point is compared with effective value parameter of corresponding moment in last cycle, if when its sudden change value all exceedes the threshold value of setting, then using the absolute time of first sampled point in these 3 sampled points interim for this week as fault moment, if do not find 3 sampled points meeting above-mentioned condition, then will start the absolute time of power frequency record ripple operation as fault moment;
Both-end resistance algorithm result L in step 3
mdz, L
ndzobtain by writing the differential equation to the voltage and current amount row in the power frequency recorder data file of both sides; And single-ended impedance range measurement L
msz, L
nszbeing then by calculating the voltage and current amount in the power frequency recorded wave file of both sides respectively, obtaining the resistance value of circuit when fault occurs and obtaining compared with the line impedance of unit length;
The concrete grammar (for M end) of the Fault Location With Traveling Wave in step 4 is:
1, intercepting time window t is established
m: with the fault moment t provided in step 3
0mfor benchmark, frontly push away Ams, pusher Bms, establish t
m∈ [t
0m-A, t
0m+ B] absolute time section as intercepting time window, wherein A and B is adjustable parameters, and is arithmetic number, and default setting is A=15, B=10;
2, Mintrop wave, reflection wave absolute time is established: to high frequency recorded wave file y
min be in and intercept time window t
minterior high-frequency current data carry out phase-model transformation, wavelet transformation, obtain wavelet modulus maxima, and carry out amplitude screening to the modulus maximum that wavelet transformation obtains, and determine that the Mintrop wave head of initial row ripple arrives the absolute time t of local terminal
1the absolute time t of local terminal is arrived with trouble spot reflection wave wave head
2; In like manner to high frequency recorded wave file y
nin be in intercepting time t
nin data process equally, obtain this fault produce initial row ripple Mintrop wave head arrive N end absolute time t
3and trouble spot reflection wave wave head arrives the absolute time t of N end
4;
3, according to described absolute time t
1and t
3mistiming calculate trouble spot distance M hold both-end distance measuring result be:
Wherein L
mnfor circuit MN total length, v is the speed of row wave traveling;
4, according to described absolute time t
1and t
3mistiming calculate trouble spot distance N hold both-end distance measuring result be:
5, according to described absolute time t
1and t
2mistiming calculate trouble spot and apart from the result of Single Terminal Traveling Wave Fault Location of M end be:
6, according to described absolute time t
3and t
4mistiming calculate trouble spot and apart from the result of Single Terminal Traveling Wave Fault Location of N end be:
Power frequency recorded wave file x in described step 4
mwith corresponding high frequency recorded wave file y
mbetween mapping relations establish in the following way:
First with power frequency recorded wave file time and x
mfor benchmark, frontly push away 20ms, pusher 20ms, form an absolute time sweep interval [x
m-20, x
m+ 20];
Then scanning search is positioned at all high frequency recorded wave files in this interval, if there is multiple high speed recorded wave file, then gets and x
mthe minimum high frequency recorded wave file of the absolute value of mistiming as y
m.
In described step 5, M end is according to impedance method range measurement L
mdz, L
mszwith Fault Location With Traveling Wave result L
mdx, L
msxthe row ripple hybrid ranging result A of M end is obtained after comprehensively analyzing
mconcrete steps as follows:
1. if both-end resistance algorithm result L
mdz, both-end travelling wave ranging result L
mdxall exist, and range measurement all L within line length
mdz∈ [0, L
mn], L
mdx∈ [0, L
mn]:
Calculate the absolute value L of the difference of both-end resistance algorithm result and both-end travelling wave ranging result
mdfor:
L
md=|L
mdz-L
mdx|
If 0≤L
md≤ (L
mn× i%): then row ripple hybrid ranging result A
mget row ripple both-end distance measuring result L
mdx, otherwise row ripple hybrid ranging result A
mget impedance both-end distance measuring result L
mdz;
Note: i is adjustable parameters, acquiescence i=2.
2. if both-end resistance algorithm result L
mdz∈ [0, L
mn], both-end travelling wave ranging result L
mdxfor outside district, then row ripple hybrid ranging result A
mget both-end resistance algorithm result L
mdz;
3. if both-end travelling wave ranging result L
mdx∈ [0, L
mn], both-end resistance algorithm result L
mdzfor outside district, then row ripple hybrid ranging result A
mget both-end travelling wave ranging result L
mdx;
4. if both-end resistance algorithm result L
mdzfor outside district, both-end travelling wave ranging result L
mdxfor outside district, and single-ended impedance range measurement L
mszmeet L
msz∈ [0, L
mn], Single Terminal Traveling Wave Fault Location result L
msxmeet L
msx∈ [0, L
mn], then carry out following judgement:
If L
msz≤ (L
mn× j%), then row ripple hybrid ranging result A
mget L
msz, otherwise row ripple hybrid ranging result A
mget L
msx;
Note: j is adjustable parameters, acquiescence j=10.
5. if both-end resistance algorithm result L
mdzfor outside district, both-end travelling wave ranging result L
mdxfor outside district, and single-ended impedance range measurement L
mszmeet L
msz∈ [0, L
mn], Single Terminal Traveling Wave Fault Location result L
msxfor outside district, then row ripple hybrid ranging result A
mget single-ended impedance range measurement L
msz;
6. if both-end resistance algorithm result L
mdzfor outside district, both-end travelling wave ranging result L
mdxfor outside district, and Single Terminal Traveling Wave Fault Location result L
msxmeet L
msx∈ [0, L
mn], single-ended impedance range measurement L
mszfor outside district, then row ripple hybrid ranging result A
mget Single Terminal Traveling Wave Fault Location result L
msx;
7. if both-end resistance algorithm result L
mdzfor outside district, both-end travelling wave ranging result L
mdxfor outside district, Single Terminal Traveling Wave Fault Location result L
msxfor outside district, single-ended impedance range measurement L
mszfor outside district, then row ripple hybrid ranging result A
mfor outside district.
In order to method of the present invention is better described, carry out simulation operations, specific as follows:
Fig. 2 is emulation test system model schematic.
Transmission line of electricity MN as shown in the figure, line length is 41.634km, suppose (to hold 29.867km apart from M at circuit F, 11.767km is held apart from N) place's generation C phase short circuit grounding fault, power frequency startup amount triggers the operation of record ripple record, obtain power frequency recorded wave file and high frequency recorded wave file, the recorder data of both sides exchange afterwards completes the preliminary work to this fault data source, and its corresponding waveform is respectively as shown in Fig. 3, Fig. 4, Fig. 5, Fig. 6.
For M station, be illustrated in figure 3 M to stand this side power frequency recorder data waveform, judge that obtaining M holds faulty line to be M line according to power frequency startup amount, according to the analysis of impedance method to this side power frequency recorder data, to obtain this side faulty line be the time is 2014/09/0814:46:50:618:199, then this time of providing with impedance method of traveling wave method is for benchmark, and the front 15ms of pushing away, pusher 10ms obtain this side and intercept time window t
minitial time be 2014/09/0814:46:50:603:199, be 2014/09/0814:46:50:628:199 by the time, between two verniers, be corresponding intercepting time window t as shown in Figure 4
m.In like manner M end to the startup amount of the offside power frequency recorder data that offside sends carry out judgement obtain N hold faulty line be N line, carrying out impedance method analysis to offside power frequency recorder data, to obtain offside fault-time be 2014/09/0814:46:50:617:997, this fault-time that traveling wave method is given with impedance method is for benchmark, before push away 15ms, pusher 10ms obtains the intercepting time window t of offside high frequency fault data
ninitial time be 2014/09/0814:46:50:602:997, be 2014/09/0814:46:50:627:997 by the time, between two verniers, be corresponding intercepting time window t as shown in Figure 6
n.
If Fig. 7 is that M holds both-end resistance algorithm schematic diagram, M holds both-end resistance algorithm length to be 29.402km as seen from the figure, and error is 0.465km; Single-ended impedance range finding length is 30.212km, and error is 0.345km.
If Fig. 8 is that N holds both-end resistance algorithm schematic diagram, N holds both-end resistance algorithm length to be 12.232km as seen from the figure, and error is 0.465km; Single-ended impedance range finding length is 9.656km, and error is 2.111km.
M end is establishing the intercepting time window t of this side and offside high frequency fault data
mand t
nafter, traveling wave method with respective intercepting time window for scope, the both sides high frequency fault data got local terminal respectively obtain M after carrying out the process such as wavelet transformation, phase-model transformation and hold row ripple Mintrop wave head and reflection wave wave head time in this side and offside fault data, calculate M side both-end travelling wave ranging length afterwards, in like manner N end also carries out similar process and obtains N and hold both-end distance measuring length.
Note: when both-end travelling wave ranging length exists (outside Bu Wei district), automatic algorithms does not carry out Single Terminal Traveling Wave Fault Location flow process.
Being illustrated in figure 9 M and holding manual row ripple single end distance measurement schematic diagram, is 30.1638km/11.4702km by the single-ended manual range finding length obtaining shown in figure holding apart from M trouble spot, and actual fault point position is 29.867km apart from M end, and error is 0.2968km.
As shown in Figure 10 for N holds manual Single Terminal Traveling Wave Fault Location schematic diagram, be 11.1363km/30.4977km by the single-ended manual range finding length obtaining shown in figure holding apart from N trouble spot, actual fault point position is 11.767km apart from N end, and error is 0.6307km.
As shown in figure 11 for M holds both-end travelling wave ranging schematic diagram, both-end travelling wave ranging length is hold 30.1095km apart from M as seen from the figure, and actual fault point is hold 29.867km apart from M, and error is 0.2425km.
As shown in figure 12 for N holds both-end travelling wave ranging schematic diagram, both-end travelling wave ranging length is hold 11.5245km apart from N as seen from the figure, and actual fault point is hold 11.767km apart from N, and error is 0.2425km.
Comprehensive the above, for M end, both-end resistance algorithm length L
mdzfor 29.402km, both-end travelling wave ranging length L
mdxfor 30.1095km, then the difference absolute value L both
mdfor 0.7075km, be less than 2% of total track length, then row ripple hybrid ranging length gets both-end travelling wave ranging length L
mdxvalue, i.e. A
mvalue is 30.1095km, and error is 0.2425km; In like manner N holds hybrid ranging length An value to be 11.5245km, and error is 0.2425km.
By reference to the accompanying drawings the specific embodiment of the present invention is described although above-mentioned; but not limiting the scope of the invention; one of ordinary skill in the art should be understood that; on the basis of technical scheme of the present invention, those skilled in the art do not need to pay various amendment or distortion that creative work can make still within protection scope of the present invention.