CN107390100B - A kind of power cable partial discharge positioning method based on time reversal - Google Patents
A kind of power cable partial discharge positioning method based on time reversal Download PDFInfo
- Publication number
- CN107390100B CN107390100B CN201710657493.XA CN201710657493A CN107390100B CN 107390100 B CN107390100 B CN 107390100B CN 201710657493 A CN201710657493 A CN 201710657493A CN 107390100 B CN107390100 B CN 107390100B
- Authority
- CN
- China
- Prior art keywords
- power cable
- partial discharge
- discharge signal
- tested
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
- G01R31/1227—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials
- G01R31/1263—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials of solid or fluid materials, e.g. insulation films, bulk material; of semiconductors or LV electronic components or parts; of cable, line or wire insulation
- G01R31/1272—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials of solid or fluid materials, e.g. insulation films, bulk material; of semiconductors or LV electronic components or parts; of cable, line or wire insulation of cable, line or wire insulation, e.g. using partial discharge measurements
Abstract
The invention discloses a kind of power cable partial discharge positioning method based on time reversal carries out transmission function test to tested power cable first, and the attenuation constant and phase constant of power cable are found out according to the transmission function that test obtains；Shelf depreciation test is carried out to tested power cable again, obtains incident Partial discharge signal and reflection Partial discharge signal；Incident Partial discharge signal and reflection Partial discharge signal are modified using attenuation constant and phase constant later, according to revised incident Partial discharge signal and reflection Partial discharge signal, entire power cable position-energy spectrogram is obtained, then the energy maximum point in spectrogram is the partial discharge position of power cable.The present invention is without estimating the time delay of local discharge signal, and the influence by the frequency dependent characteristic of power cable phase velocity to power cable breakdown location is taken into account, to improve power cable shelf depreciation positioning accuracy.
Description
Technical field
The invention belongs to telecommunication diagnostic techniques fields, are related to the running shelf depreciation diagnostic techniques of power cable, more
Specifically, being related to a kind of power cable partial discharge positioning method based on time reversal.
Background technique
Cross-inked polyethylene power cable is largely used in power grid because of its reliable electrical, mechanical performance；But manufacture,
The complexity of abnormal activity present in transport and installation process and power cable service condition, will lead to power cable
There is power cable defect in local location；In this way under forceful electric power field action, it is existing that power cable defective locations can generate shelf depreciation
As.Shelf depreciation is one of the form of expression of power cable insulation aging, while also making to insulate and further deteriorating and lead to electric power
One of the main reason for cable insulation fails.Therefore, in order to ensure power cable reliability of operation, need to power cable into
Row Partial Discharge Detection.
In the prior art, based on mature time domain reflectometry, by detection power cable partial discharge pulse signal come real
Existing power cable shelf depreciation positioning.Number of patent application is that the application documents of CN201510173365.9 disclose one kind based on vibration
The distribution power cable local discharge measuring method for swinging wave test initially sets up the power cable based on oscillation wave test and locally puts
Then electrical testing and positioning system utilize above-mentioned test macro, the method by being forced into preset value at tested power cable end
It charges to the power cable, and resonance is occurred in tested electric power electricity by inductance in system and tested power cables capacitance
Cable-end generates damped oscillation voltage；Two pulses of incidence wave and back wave according to power cable discharge position reach test lead
Time difference determines the position that shelf depreciation occurs.Number of patent application is that the application documents of CN201510092868.3 disclose
A kind of power cable partial discharge positioning method analyzed based on auto-correlation-Wavelet Modulus Maxima, first acquisition shelf depreciation letter
Number, incidence wave-reflected wave pulses data segment is then extracted in the local discharge signal of acquisition, then anti-to the incidence wave-of extraction
Ejected wave pulse data section carries out autocorrelation analysis and Wavelet Modulus Maxima is analyzed, and obtains power cable shelf depreciation point position.On
State two kinds of power cable breakdown location methods be based on incidence wave and back wave time delay estimadon combine fixed spread speed into
Row positioning, but due to accurately can not carry out time delay estimation and do not consider power cable phase velocity frequency dependent characteristic influence, and
Two kinds of power cable breakdown location methods are made to be unable to reach higher positioning accuracy, to increase the operation and maintenance of power cable
Cost.
Therefore, the power cable shelf depreciation location technology with higher positioning accuracy is developed, power cable is reduced
Operation expense is significant in booming power industry.
Summary of the invention
For the not high problem of current power cable partial discharge positioning method positioning accuracy, the purpose of the present invention is intended to mention
For a kind of power cable partial discharge positioning method based on time reversal, to realize to the accurate of power cable Partial Discharge Sources
Positioning.
Thinking of the invention is handled local discharge signal using time reversal, and from treated, shelf depreciation is believed
Incident Partial discharge signal and reflection Partial discharge signal number are isolated, and incident Partial discharge signal and reflection Partial discharge signal are utilized in frequency domain
The attenuation constant and phase constant of tested power cable are modified, so that the frequency dependent characteristic influence of power cable phase velocity be examined
Including worry；Revised incident Partial discharge signal and reflection Partial discharge signal are converted to time domain again,
Find out the signal energy that tested power cable assumes partial discharge source position；Traversal is entirely tested power cable, obtains it
Position-energy spectrogram, the energy maximum point in spectrogram are the partial discharge position of power cable.
Power cable partial discharge positioning method provided by the invention based on time reversal, the composition packet of technical solution
Include following steps:
The test of S1 power cable transmission function, finds out the attenuation constant α (ω) and phase constant β (ω) of power cable
Transmission function is carried out with the reference power cable of batch with model to tested power cable or with tested power cable
e-(α(ω)+jβ(ω))lTest, l be tested power cable or with reference to power cable length,And the biography obtained according to test
Delivery function finds out the attenuation constant α (ω) and phase constant β (ω) of power cable；
S2 is tested the test of power cable shelf depreciation
Shelf depreciation test is carried out to tested power cable, and obtained local discharge signal s (t) the progress time is anti-
It drills, obtains time reversal local discharge signal s (T-t), then the local discharge signal s (T-t) after time reversal is separated,
Obtain incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t)；Wherein t is the time, and T is appointing no more than testing time length
Meaning constant；
S3 is tested the positioning of power cable shelf depreciation
S31 utilizes Fast Fourier Transform (FFT) by incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t) frequency domain is transformed to,
Obtain the incident Partial discharge signal V of frequency domain1(ω) and reflection Partial discharge signal V2(ω)；
S32 assumes a partial discharge source position x, the attenuation constant α (ω) obtained using step S1 on tested power cable
With phase constant β (ω) to incident Partial discharge signal V1(ω) and reflection Partial discharge signal V2(ω) is modified, obtain it is revised enter
Penetrate Partial discharge signal V3(ω) and reflection Partial discharge signal V4(ω):
V3(ω)=V1(ω)e-(α(ω)+jβ(ω))x
V4(ω)=V2(ω)e-(α(ω)+jβ(ω))(2l-x)
Wherein, V3(ω)、V4(ω) is respectively incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t) x is propagated respectively
Frequency-domain expression with 2l-x apart from rear signal；L is tested power cable length；e-(α(ω)+jβ(ω))xFor incoming signal v1(t) it passes
Transmission function expression formula after broadcasting x distance, e-(α(ω)+jβ(ω))(2l-x)To reflect signal v2(t) the transmitting letter after propagating 2l-x distance
Number expression formula；
S33 utilizes Fast Fourier Transform Inverse by revised incident Partial discharge signal V3(ω) and reflection Partial discharge signal V4
(ω) transforms to time domain, obtains the incident Partial discharge signal v that power cable assumes time domain at the x of partial discharge source position3(t) and reflection partial discharge
Signal v4(t)；
S34 assumes the incident Partial discharge signal v of time domain at the x of partial discharge source position based on tested power cable3(t) and reflection partial discharge
Signal v4(t), the signal energy p (x) that power cable in testing time length assumes partial discharge source position x is obtained:
P (x)=∫ (| v3(t)|+|v4(t)|)2dt；
S35 repeats step S32 to step S34, and the entire power cable that traversal length is l obtains entire power cable position
- energy spectrogram is set, then the energy maximum point in spectrogram is the partial discharge position of power cable.
The above-mentioned power cable partial discharge positioning method based on time reversal can use time domain method or frequency in step S1
Domain method carries out transmission function test with the reference power cable of batch with model to tested power cable or with tested power cable,
Specific test can be with reference to the test method that this field has disclosed (referring to Papazyan R, Eriksson
R.Calibration for time domain propagation constant measurements on power
Cables [J] .IEEE Transactions on Instrumentation and Measurement, 2003,52 (2):
415-418)；The positioning of power cable shelf depreciation, frequency domain test equipment commonly used in the art are Network Analyzer (such as Agilent
E5061B), before being tested using the equipment, passed to eliminate the test coaxial cable connecting with power cable to power cable
The influence for passing characteristic needs to calibrate (referring to Joel P D.Handbook of Microwave test coaxial cable
Component Measurements with Advanced VNA Techniques [M] .United Kingdom:John
Wiley&Sons, 2012).The transmitting letter that will further can be tested after calibration using the time-domain analysis function of Network Analyzer
Number is transformed into time domain, and the time domain gating function of being carried by Network Analyzer selects reflection power cable transmitting in the time domain
The useful signal of function is influenced with reducing the impedance mismatch of test coaxial cable and power cable junction.It is surveyed in the present invention
Trying obtained transmission function is e-(α(ω)+jβ(ω))l, l is to be tested power cable or refer to power cable length, then foundation through time domain
Attenuation constant α (ω) and phase constant β (ω) that obtained transmission function finds out power cable are gated, wherein ω is angular frequency.
The above-mentioned power cable partial discharge positioning method based on time reversal can use this field in step s 2
The partial discharge test system (such as Germany match bar OWTS HV250) having disclosed carries out shelf depreciation survey to tested power cable
Examination.The method that obtained local discharge signal is carried out time reversal processing, can be with reference to the place disclosed in the prior art
Reason method, such as Fink M, Prada C, Wu F, et al, Self focusing with time reversal mirror in
Inhomogeneous media [J] .IEEE Ultrasonic Symp., 1989:681-686.
The above-mentioned power cable partial discharge positioning method based on time reversal in step S35, is selected using equidistant mode
Tested power cable is taken to assume partial discharge source position.It in a preferred embodiment, can be first using the entire electricity of the first spacing traversal
Power cable chooses tested power cable and assumes partial discharge source position, repeats step S32 to step S34, obtains entire power cable position
Set-energy spectrogram；Power cable length is tested before and after very big position in spectrogram using less than the second spacing of the first spacing again
Tested power cable is further chosen in 1%~5% range and assumes partial discharge source position, is repeated step S32 to step S34, is obtained
Very big position front and back is tested position-energy spectrogram in 1%~5% range of power cable length, energy in the map in spectrogram
Maximum point is the partial discharge position of power cable.
Power cable partial discharge positioning method at present is based primarily upon time delay estimadon technology, since there are time delays to estimate
Meter error and the frequency dependent characteristic influence for not considering power cable phase velocity, and lead to power cable shelf depreciation positioning accuracy not
It is high.
Compared with prior art, the power cable partial discharge positioning method provided by the invention based on time reversal, tool
There are following very prominent and advantageous effects:
1, the present invention determines power cable Partial Discharge Sources using the method for solving power cable signal energy extreme value
Position is estimated without the time delay to local discharge signal, to avoid because error existing for time delay estimation is to positioning accuracy
Influence；
2, the present invention combines the transmission function of power cable to carry out propagation emulation to power cable Partial discharge signal, by electric power electricity
Influence of the frequency dependent characteristic of cable phase velocity to power cable breakdown location is taken into account, so that it is fixed to improve power cable shelf depreciation
Position precision；
3, the present invention carries out time reversal processing to the power cable partial discharge signal of test, by local discharge signal
Dispersion phenomenon is inverted, and the dispersion phenomenon at partial discharge source position is offset, to reduce dispersion to power cable part
The influence of breakdown location precision；
4, the present invention using the incident Partial discharge signal isolated from local discharge signal and reflection Partial discharge signal be superimposed into
Row positioning, can effectively improve the power cable shelf depreciation positioning accuracy under Low SNR；
5, the present invention is sampled due to needing to carry out positioning analysis to entire power cable by Partial Discharge Measurement System
Rate influence is smaller, and the accurate positioning of power cable shelf depreciation is still able to achieve in the lower situation of sample rate.
Detailed description of the invention
Fig. 1 is that the present invention is based on the power cable partial discharge positioning method flow charts of time reversal.
Fig. 2 is the tested power cable shelf depreciation positioning step flow chart of the present invention.
Fig. 3 is test coaxial cable calibration front and back transmission characteristic in the embodiment of the present invention 1.
Fig. 4 is the attenuation constant and phase constant curve that obtained reference power cable is tested in the embodiment of the present invention 1,
Middle solid line indicates that attenuation constant curve, dotted line indicate phase constant curve.
Fig. 5 is to obtain partial discharge pulse using oscillation wave partial discharge test system and test in the embodiment of the present invention 1 to believe
Number.
Fig. 6 is to obtain partial discharge pulse using oscillation wave partial discharge test system and test in the embodiment of the present invention 1 to believe
The partial discharge pulse's signal obtained after number time reversal processing.
Fig. 7 is the power cable position-energy spectrogram drawn in the embodiment of the present invention 1.
Fig. 8 is the power cable position-energy spectrogram drawn in the embodiment of the present invention 2.
Specific embodiment
The embodiment of the present invention is provided below with reference to attached drawing, and technical solution of the present invention is carried out into one by embodiment
Clear, the complete explanation of step.Obviously, the embodiment is only a part of the embodiments of the present invention, rather than whole implementation
Example.Based on the content of present invention, those of ordinary skill in the art are obtained all without making creative work
Other embodiments belong to the range that the present invention is protected.
Embodiment 1
The tested power cable that the present embodiment is directed to is that the shelf depreciation of 498m ZR-YJV22-8.7/15 power cable is fixed
Position, power cable defect type are that semi-conductive layer overlaps bad defect.
The present embodiment carries out transmission function test to reference power cable using Agilent E5061B network tester.
The present embodiment using Germany match bar OWTS HV250 oscillation wave partial discharge test macro to tested power cable into
The test of row shelf depreciation.
Power cable partial discharge positioning method provided in this embodiment based on time reversal, as shown in Figures 1 and 2,
The following steps are included:
The test of S1 power cable transmission function
By with the same model of tested power cable, with batch length be l0=10 meters of reference power cable both ends lead to respectively
The signal output end and signal receiving end for crossing test coaxial cable and network tester connect, under the conditions of 100kHz-50MHz,
It is tested using transmission function of the Network Analyzer to reference power cable.School is carried out to test coaxial cable first before test
Standard eliminates influence of the test coaxial cable to reference power cable transmission characteristic；Test coaxial cable transmission characteristic calibration front and back
As shown in Figure 3；Coaxial cable is tested before calibration, in test result and shows apparent transmission characteristic, after calibration, tests coaxial electrical
The transmission characteristic of cable is eliminated.Then reference power cable is once tested using Network Analyzer, it is coaxial due to testing
Cable and refer to power cable impedance mismatch, the two junction exist reflection signal, it is therefore desirable to will test obtain just
Pacing test result is selected using the time-domain analysis power and energy of Network Analyzer to time domain, and using the time domain that Network Analyzer carries
Logical function chooses the useful signal that reflection refers to power cable transmission function, reduces coaxial cable and refers to power cable junction
Impedance mismatch influence.Finally utilize the power cable transmission function obtained through time domain gating(l0For ginseng
Examine power cable length) find out attenuation constant α (ω) and phase constant β (ω), the attenuation constant and phase with reference to power cable
Bit constant and the attenuation constant of tested power cable are identical as phase constant, with reference to the attenuation constant and phase constant of power cable
Curve is as shown in Figure 4.
S2 is tested the test of power cable shelf depreciation
Shelf depreciation test, the office tested are carried out to tested power cable using oscillation wave partial discharge test macro
Portion's Discharge pulse waveform is as shown in figure 5, the local discharge signal s (t) that then will test carries out time reversal processing t → T-
T, local discharge signal s (T-t) that inverting obtains as shown in fig. 6, again by the local discharge signal s (T-t) after time reversal into
Row separation, obtains incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t)；Power cable partial discharge signal is through time reversal
After processing, the dispersion phenomenon of local discharge signal is inverted, and is offset the dispersion phenomenon at partial discharge source position, in this way may be used
To reduce influence of the dispersion to power cable shelf depreciation positioning accuracy.
S3 is tested the positioning of power cable shelf depreciation
S31 utilizes Fast Fourier Transform (FFT) by incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t) frequency domain is transformed to,
Obtain the incident Partial discharge signal V of frequency domain1(ω) and reflection Partial discharge signal V2(ω)；
V1(ω)=FFT [v1(t)]
V2(ω)=FFT [v2(t)]
S32 assumes a partial discharge source position x, the attenuation constant α (ω) obtained using step S1 on tested power cable
With phase constant β (ω) to incident Partial discharge signal V1(ω) and reflection Partial discharge signal V2(ω) is modified, obtain it is revised enter
Penetrate Partial discharge signal V3(ω) and reflection Partial discharge signal V4(ω):
V3(ω)=V1(ω)e-(α(ω)+jβ(ω))x
V4(ω)=V2(ω)e-(α(ω)+jβ(ω))(2l-x)
Wherein, V3(ω)、V4(ω) is respectively incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t) x is propagated respectively
Frequency-domain expression with 2l-x apart from rear signal；L is tested power cable length；e-(α(ω)+jβ(ω))xFor incoming signal v1(t) it passes
Transmission function expression formula after broadcasting x distance, e-(α(ω)+jβ(ω))(2l-x)To reflect signal v2(t) the transmitting letter after propagating 2l-x distance
Number expression formula；
S33 utilizes Fast Fourier Transform Inverse by revised incident Partial discharge signal V3(ω) and reflection Partial discharge signal V4
(ω) transforms to time domain, obtains the incident Partial discharge signal v that power cable assumes time domain at the x of partial discharge source position3(t) and reflection partial discharge
Signal v4(t)；
v3(t)=IFFT [V3(ω)]
v4(t)=IFFT [V4(ω)]
S34 assumes the incident Partial discharge signal v of time domain at the x of partial discharge source position based on tested power cable3(t) and reflection partial discharge
Signal v4(t), the signal energy p (x) that power cable in testing time length assumes partial discharge source position x is obtained:
P (x)=∫ (| v3(t)|+|v4(t)|)2dt；
S35 repeats step S32 to step S34, using 0.5 meter of the first spacing, changes hypothesis office on tested power cable
Source position x is put, step S32 to step S34 is repeated, so that it is traversed length l=498 meters of entire tested power cable, seeks
Each energy p (x) for assuming partial discharge source position, and make mutually it will be assumed that the entire power cable of partial discharge source position and energy relationship
Position-energy spectrogram；For the office within the scope of 5 meters before and after position very big in entire power cable position-energy spectrogram of acquisition
Portion is tested power cable, using 0.01 meter of the second spacing, changes on locally tested power cable and assumes partial discharge source position x, weight
Multiple step S32 to step S34, the part for choosing its traversal are tested power cable, seek each energy for assuming partial discharge source position
It measures p (x), and makes mutually it will be assumed that partial discharge source position and the part of energy relationship are tested power cable position-energy spectrogram, then office
The energy maximum point that portion is tested in power cable position-energy spectrogram is the partial discharge position of power cable.
Power cable position-energy spectrogram that the present embodiment obtains is as shown in fig. 7, as can be seen from the figure in tested electric power
There are an energy maximum at cable 249m, exactly tested power cable partial discharge position at this.
Embodiment 2
The tested power cable that the present embodiment is directed to is that the shelf depreciation of 498m ZR-YJV22-8.7/15 power cable is fixed
Position, power cable defect type are longitudinal tool marks defect.
The present embodiment carries out transmission function test to reference power cable using Agilent E5061B network tester.
The present embodiment using Germany match bar OWTS HV250 oscillation wave partial discharge test macro to tested power cable into
The test of row shelf depreciation.
Power cable partial discharge positioning method provided in this embodiment based on time reversal, as shown in Figures 1 and 2,
The following steps are included:
The test of S1 power cable transmission function
By with the same model of tested power cable, with batch length be l0=10 meters of reference power cable both ends lead to respectively
The signal output end and signal receiving end for crossing test coaxial cable and network tester connect, under the conditions of 100kHz-50MHz,
It is tested using transmission function of the Network Analyzer to reference power cable.School is carried out to test coaxial cable first before test
Standard eliminates influence of the test coaxial cable to reference power cable transmission characteristic.Then using Network Analyzer to reference electric power
Cable is once tested, and due to test coaxial cable and refers to power cable impedance mismatch, exists in the two junction anti-
Penetrate signal, it is therefore desirable to will the obtained the tested results of test using Network Analyzer time-domain analysis power and energy to when
Domain, and the useful signal that reflection refers to power cable transmission function is chosen using the time domain gating function that Network Analyzer carries,
Reduce coaxial cable and is influenced with reference to the impedance mismatch of power cable junction.Finally utilize the electric power obtained through time domain gating
Cable transmission function(l0To refer to power cable length) find out attenuation constant α (ω) with reference to power cable
With phase constant β (ω), the attenuation constant and phase constant and the attenuation constant of tested power cable are identical as phase constant.
S2 is tested the test of power cable shelf depreciation
Shelf depreciation test is carried out to tested power cable using oscillation wave partial discharge test macro, then will test
The local discharge signal s (t) that arrives carries out time reversal processing t → T-t, then by the local discharge signal s (T-t) after time reversal
It is separated, obtains incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t)；Power cable partial discharge signal is anti-through the time
After drilling processing, the dispersion phenomenon of local discharge signal is inverted, and offsets the dispersion phenomenon at partial discharge source position, in this way
It can reduce influence of the dispersion to power cable shelf depreciation positioning accuracy.
S3 is tested the positioning of power cable shelf depreciation
S31 utilizes Fast Fourier Transform (FFT) by incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t) frequency domain is transformed to,
Obtain the incident Partial discharge signal V of frequency domain1(ω) and reflection Partial discharge signal V2(ω)；
V1(ω)=FFT [v1(t)]
V2(ω)=FFT [v2(t)]
S32 assumes a partial discharge source position x, the attenuation constant α (ω) obtained using step S1 on tested power cable
With phase constant β (ω) to incident Partial discharge signal V1(ω) and reflection Partial discharge signal V2(ω) is modified, obtain it is revised enter
Penetrate Partial discharge signal V3(ω) and reflection Partial discharge signal V4(ω):
V3(ω)=V1(ω)e-(α(ω)+jβ(ω))x
V4(ω)=V2(ω)e-(α(ω)+jβ(ω))(2l-x)
Wherein, V3(ω)、V4(ω) is respectively incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t) x is propagated respectively
Frequency-domain expression with 2l-x apart from rear signal；L is tested power cable length；e-(α(ω)+jβ(ω))xFor incoming signal v1(t) it passes
Transmission function expression formula after broadcasting x distance, e-(α(ω)+jβ(ω))(2l-x)To reflect signal v2(t) the transmitting letter after propagating 2l-x distance
Number expression formula；
S33 utilizes Fast Fourier Transform Inverse by revised incident Partial discharge signal V3(ω) and reflection Partial discharge signal V4
(ω) transforms to time domain, obtains the incident Partial discharge signal v that power cable assumes time domain at the x of partial discharge source position3(t) and reflection partial discharge
Signal v4(t)；
v3(t)=IFFT [V3(ω)]
v4(t)=IFFT [V4(ω)]
S34 assumes the incident Partial discharge signal v of time domain at the x of partial discharge source position based on tested power cable3(t) and reflection partial discharge
Signal v4(t), the signal energy p (x) that power cable in testing time length assumes partial discharge source position x is obtained:
P (x)=∫ (| v3(t)|+|v4(t)|)2dt；
S35 repeats step S32 to step S34, using 0.5 meter of the first spacing, changes hypothesis office on tested power cable
Source position x is put, step S32 to step S34 is repeated, so that it is traversed the entire tested power cable of length l, seeks each hypothesis
The energy p (x) of partial discharge source position, and make mutually it will be assumed that entire power cable position-energy of partial discharge source position and energy relationship
Measure spectrogram；For the tested electricity in part within the scope of 5 meters before and after position very big in entire power cable position-energy spectrogram of acquisition
Power cable is changed on locally tested power cable using 0.01 meter of the second spacing and assumes partial discharge source position x, repeated step S32
To step S34, the part for choosing its traversal is tested power cable, seeks each energy p (x) for assuming partial discharge source position, and
It makes mutually it will be assumed that partial discharge source position and the part of energy relationship are tested power cable position-energy spectrogram, then part is tested electric
Energy maximum point in power cable position-energy spectrogram is the partial discharge position of power cable.
Power cable position-energy spectrogram that the present embodiment obtains is as shown in figure 8, as can be seen from the figure in tested electric power
There are an energy maximum at cable 249m, exactly tested power cable partial discharge position at this.
Those of ordinary skill in the art will understand that embodiment here be to help reader understand it is of the invention
Principle, it should be understood that protection scope of the present invention is not limited to such specific embodiments and embodiments.This field it is common
Technical staff disclosed the technical disclosures can make the various various other tools for not departing from essence of the invention according to the present invention
Body variations and combinations, these variations and combinations are still within the scope of the present invention.
Claims (4)
1. a kind of power cable partial discharge positioning method based on time reversal, it is characterised in that the following steps are included:
The test of S1 power cable transmission function, finds out the attenuation constant α (ω) and phase constant β (ω) of power cable
Transmission function e is carried out with the reference power cable of batch with model to tested power cable or with tested power cable-(α(ω)+jβ(ω))lTest, l be tested power cable or with reference to power cable length,And the biography obtained according to test
Delivery function finds out the attenuation constant α (ω) and phase constant β (ω) of power cable；
S2 is tested the test of power cable shelf depreciation
Shelf depreciation test is carried out to tested power cable, and obtained local discharge signal s (t) is subjected to time reversal, is obtained
It separates, obtains to time reversal local discharge signal s (T-t), then by the local discharge signal s (T-t) after time reversal
Incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t)；Wherein t is the time, and T is any normal no more than testing time length
Number；
S3 is tested the positioning of power cable shelf depreciation
S31 utilizes Fast Fourier Transform (FFT) by incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t) frequency domain is transformed to, is obtained
The incident Partial discharge signal V of frequency domain1(ω) and reflection Partial discharge signal V2(ω)；
S32 assumes a partial discharge source position x on tested power cable, utilizes the obtained attenuation constant α (ω) of step S1 and phase
Bit constant β (ω) is to incident Partial discharge signal V1(ω) and reflection Partial discharge signal V2(ω) is modified, and obtains revised incident office
Discharge signal V3(ω) and reflection Partial discharge signal V4(ω):
V3(ω)=V1(ω)e-(α(ω)+jβ(ω))x
V4(ω)=V2(ω)e-(α(ω)+jβ(ω))(2l-x)
Wherein, V3(ω)、V4(ω) is respectively incident Partial discharge signal v1(t) and reflection Partial discharge signal v2(t) x and 2l-x is propagated respectively
Frequency-domain expression apart from rear signal；L is tested power cable length；e-(α(ω)+jβ(ω))xFor incoming signal v1(t) propagate x away from
Transmission function expression formula from after, e-(α(ω)+jβ(ω))(2l-x)To reflect signal v2(t) the transmission function table after propagating 2l-x distance
Up to formula；
S33 utilizes Fast Fourier Transform Inverse by revised incident Partial discharge signal V3(ω) and reflection Partial discharge signal V4(ω) becomes
Time domain is changed to, the incident Partial discharge signal v that power cable assumes time domain at the x of partial discharge source position is obtained3(t) and reflection Partial discharge signal v4
(t)；
S34 assumes the incident Partial discharge signal v of time domain at the x of partial discharge source position based on tested power cable3(t) and reflection Partial discharge signal
v4(t), the signal energy p (x) that power cable in testing time length assumes partial discharge source position x is obtained:
P (x)=∫ (| v3(t)|+|v4(t)|)2dt；
S35 repeats step S32 to step S34, and the entire power cable that traversal length is l obtains entire power cable position-energy
Spectrogram is measured, then the energy maximum point in spectrogram is the partial discharge position of power cable.
2. the power cable partial discharge positioning method based on time reversal according to claim 1, it is characterised in that use
Time domain method or frequency domain method carry out transmission function test to tested power cable.
3. the power cable partial discharge positioning method according to claim 1 or claim 2 based on time reversal, it is characterised in that step
In rapid S35, tested power cable is chosen using equidistant mode and assumes partial discharge source position.
4. the power cable partial discharge positioning method based on time reversal according to claim 3, it is characterised in that step
In S35, entire power cable is traversed using the first spacing and chooses tested power cable hypothesis partial discharge source position, repeats step S32
To step S34, entire power cable position-energy spectrogram is obtained；Using the second spacing less than the first spacing in spectrogram pole
It is tested the tested power cable of selection in 1%~5% range of power cable length before and after big position and assumes partial discharge source position, repeats to walk
Rapid S32 to step S34 obtains the position-energy being tested before and after very big position in 1%~5% range of power cable length in spectrogram
Spectrogram is measured, energy maximum point is the partial discharge position of power cable.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710657493.XA CN107390100B (en) | 2017-08-03 | 2017-08-03 | A kind of power cable partial discharge positioning method based on time reversal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710657493.XA CN107390100B (en) | 2017-08-03 | 2017-08-03 | A kind of power cable partial discharge positioning method based on time reversal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107390100A CN107390100A (en) | 2017-11-24 |
| CN107390100B true CN107390100B (en) | 2019-07-19 |
Family
ID=60343102
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710657493.XA Active CN107390100B (en) | 2017-08-03 | 2017-08-03 | A kind of power cable partial discharge positioning method based on time reversal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107390100B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3879284A1 (en) * | 2020-03-13 | 2021-09-15 | General Electric Technology GmbH | Fault localization system and associated fault localization method |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108089103B (en) * | 2017-12-11 | 2019-07-12 | 西南交通大学 | Vehicle-mounted cable power-frequency voltage positive half period shelf depreciation start-phase evaluation method |
| CN107907805B (en) * | 2017-12-11 | 2019-09-06 | 西南交通大学 | The vehicle-mounted cable low temperature shelf depreciation start-phase evaluation method of power-frequency voltage negative half period |
| CN108983032B (en) * | 2018-09-17 | 2021-02-12 | 华北电力大学（保定） | Power cable open-circuit fault positioning method based on time reversal |
| CN110488169A (en) * | 2019-09-24 | 2019-11-22 | 黄志彭 | A kind of high voltage power cable partial discharge positioning method |
| CN111157845B (en) * | 2020-01-07 | 2021-10-26 | 国网辽宁省电力有限公司朝阳供电公司 | Time domain inversion technology-based fault positioning method suitable for power distribution network |
| CN112083298B (en) * | 2020-09-08 | 2021-07-06 | 四川大学 | Insulating paperboard partial discharge failure probability characterization method based on double risk competition |
| CN111929553B (en) * | 2020-10-19 | 2020-12-22 | 四川大学 | Partial discharge positioning method based on phase velocity frequency-varying characteristics |
| CN112505411A (en) * | 2020-10-30 | 2021-03-16 | 山东科技大学 | Method for identifying cable transfer function on line |
| CN112505483A (en) * | 2020-11-25 | 2021-03-16 | 山东科技大学 | Cable partial discharge positioning method based on online transfer function |
| CN112557851A (en) * | 2020-12-07 | 2021-03-26 | 国网天津市电力公司电力科学研究院 | Power cable partial discharge on-line positioning system and method based on transfer function |
| CN113281613B (en) * | 2021-05-18 | 2023-03-31 | 国网江苏省电力有限公司无锡供电分公司 | Power cable transmission parameter acquisition method |
| CN113820568B (en) * | 2021-07-14 | 2024-04-02 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Partial discharge positioning method and device for cable, computer equipment and storage medium |
| CN114689997B (en) * | 2022-05-06 | 2023-09-26 | 国网江苏省电力有限公司宿迁供电分公司 | Distribution network cable fault identification and positioning method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2577457A1 (en) * | 2006-10-27 | 2008-04-27 | Artenga Inc. | Microbubble medical devices |
| CN101873692A (en) * | 2010-06-23 | 2010-10-27 | 电子科技大学 | Wireless sensor network node positioning method based on time reversal |
| CN103149513A (en) * | 2013-02-26 | 2013-06-12 | 中国科学院电工研究所 | Positioning method and device for reestablishing local discharge ultrasonic source of transformer |
| CN104639491A (en) * | 2015-01-29 | 2015-05-20 | 武汉铭为电子科技有限公司 | Time reversal orthogonal FDM (Frequency-Division Multiplexing) wireless communication method based on SFME (Solvent Free Microwave Extraction) |
| CN105093083A (en) * | 2015-08-31 | 2015-11-25 | 广州供电局有限公司 | Apparatus and method for locating cable partial discharge signals |
| CN106405341A (en) * | 2016-07-15 | 2017-02-15 | 国网甘肃省电力公司电力科学研究院 | Method used for reducing partial discharge ultrahigh frequency signal |
-
2017
- 2017-08-03 CN CN201710657493.XA patent/CN107390100B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2577457A1 (en) * | 2006-10-27 | 2008-04-27 | Artenga Inc. | Microbubble medical devices |
| CN101873692A (en) * | 2010-06-23 | 2010-10-27 | 电子科技大学 | Wireless sensor network node positioning method based on time reversal |
| CN103149513A (en) * | 2013-02-26 | 2013-06-12 | 中国科学院电工研究所 | Positioning method and device for reestablishing local discharge ultrasonic source of transformer |
| CN104639491A (en) * | 2015-01-29 | 2015-05-20 | 武汉铭为电子科技有限公司 | Time reversal orthogonal FDM (Frequency-Division Multiplexing) wireless communication method based on SFME (Solvent Free Microwave Extraction) |
| CN105093083A (en) * | 2015-08-31 | 2015-11-25 | 广州供电局有限公司 | Apparatus and method for locating cable partial discharge signals |
| CN106405341A (en) * | 2016-07-15 | 2017-02-15 | 国网甘肃省电力公司电力科学研究院 | Method used for reducing partial discharge ultrahigh frequency signal |
Non-Patent Citations (2)
| Title |
|---|
| On the use of the Time Reversal of Electromagnetic fields to locate lightning discharges;Gaspard Lugrin等;《2012 International Conference on Lightning Protection》;20120615;第1-4页 |
| 基于时间反演技术的电力电缆局部放电定位方法;谢敏;《中国电机工程学报》;20180605;第38卷(第11期);第3402-3409页 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3879284A1 (en) * | 2020-03-13 | 2021-09-15 | General Electric Technology GmbH | Fault localization system and associated fault localization method |
| WO2021180637A1 (en) * | 2020-03-13 | 2021-09-16 | General Electric Technology Gmbh | Fault localization system and associated fault localization method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107390100A (en) | 2017-11-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107390100B (en) | A kind of power cable partial discharge positioning method based on time reversal | |
| Papazyan et al. | Calibration for time domain propagation constant measurements on power cables | |
| Sheng et al. | Partial discharge pulse propagation in power cable and partial discharge monitoring system | |
| Koltunowicz et al. | Synchronous multi-channel PD measurements | |
| Song et al. | Detection and location of multiple wiring faults via time–frequency-domain reflectometry | |
| Wagenaars et al. | Measurement of transmission line parameters of three-core power cables with common earth screen | |
| Wang et al. | Electromagnetic time reversal applied to fault location: On the properties of back-injected signals | |
| Hashmi et al. | Determining wave propagation characteristics of MV XLPE power cable using time domain reflectometry technique | |
| Shi et al. | Application of iterative deconvolution for wire fault location via reflectometry | |
| CN113820568B (en) | Partial discharge positioning method and device for cable, computer equipment and storage medium | |
| Shi et al. | Detection and location of single cable fault by impedance spectroscopy | |
| Gautschi et al. | Calibration of UHF sensors for GIS: Comparison of different methods and testing of a calibration system based on a conical antenna | |
| Clark et al. | Partial discharge pulse propagation, localisation and measurements in medium voltage power cables | |
| Thayoob et al. | Analysis of wave propagation in Time Domain Reflectometry circuit simulation model | |
| Lee et al. | Impedance change localization for live underground cable using time-frequency domain reflectometry | |
| Shi et al. | A new method of locating the single wire fault | |
| Chen et al. | Comparison of location accuracy between frequency domain reflectometry and line resonance analysis for power cables | |
| Ragusa et al. | Practical Evaluation of Electromagnetic Time Reversal to Locate Partial Discharges on Power Networks in the Presence of Noise | |
| Mier et al. | Methods for Partial Discharge Calibration in Gas-Insulated Substations for HVDC Power Grids and Charge Evaluation Uncertainty | |
| Yao et al. | Cable fault location and signal separation based on continuous wavelet transform | |
| Tao et al. | Design of aircraft cable fault diagnose and location system based on aircraft airworthiness requirement | |
| Hashmi et al. | Comparing wave propagation characteristics of MV XLPE cable and covered-conductor overhead line using time domain reflectometry technique | |
| CN106053945B (en) | A kind of S function Smith chart analytical equipment and method in short-term | |
| Burghardt et al. | System Identification of a Branched 50 Ohm Network by Transient Excitation | |
| Zhao et al. | Research on Propagation Characteristics of Partial Discharge Signal in Cable Joint Based on CST |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhou Kai Inventor after: Huang Yonglu Inventor after: Xie Min Inventor after: He Min Inventor after: Zhao Shilin Inventor after: Zhang Fuzhong Inventor after: Ran Li Inventor before: Zhou Kai Inventor before: Xie Min Inventor before: He Min Inventor before: Zhao Shilin Inventor before: Zhang Fuzhong Inventor before: Huang Yonglu Inventor before: Ran Li |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |