Take into account the cable local discharge signal processing and localization method of DISCHARGE PULSES EXTRACTION and denoising
Technical field
Power cable Partial Discharge Detecting Technology of the present invention field is related to partial discharge signal wave analysis processing technique,
It is fixed more particularly to a kind of signal processing of partial discharge method for taking into account DISCHARGE PULSES EXTRACTION and signal denoising and power cable shelf depreciation
Position method.
Background technique
Shelf depreciation (partial discharge, PD) (abbreviation partial discharge) is that power cable insulation the performance deteriorated occurs
One of form, while being also one of the main reason for power cable insulation further deteriorates and leads to failure of insulation.By to electricity
Power cable carries out Partial Discharge Detection and analysis, can diagnose for power cable insulation and provide reliable judging basis with monitoring.
Traditional Partial discharge signal processing method mainly includes two major classes: 1) feature based on phase distribution mode statistical spectrogram
Extracting method, this method count the local discharge signal (abbreviation Partial discharge signal) of a period of time interior, several power frequency period first, draw
Shelf depreciation amplitude is with the X-Y scheme of phase change, i.e. shelf depreciation phase distribution (phase resolved partial out
Discharge, PRPD) figure, in order to more intuitively reflect local discharge characteristic, then from PRPD extract electric discharge operating frequency phase φ,
Discharge capacity q and discharge time n constitutes two dimension or three-dimensional spectrum (shelf depreciation fingerprint chromatogram), and by fingerprint chromatogram, extracts and use
In the characteristic quantity of characterization shelf depreciation spectrogram, such as degree of skewness Sk, kurtosis ku, degree of asymmetry cc etc.;2) based on the spy of pulse
Extracting method is levied, this method usually utilizes partial discharge pulse to extract the characteristic quantity for characterizing local discharge signal feature, such as
Equivalent time width T, equivalent bandwidth F etc., different type Partial discharge signal can effectively be identified by extracting different characteristic quantity and
Effective assessment to power cable insulation state can be achieved.However, in above two Partial discharge signal processing method, with partial discharge arteries and veins
It purges with and is taken as to obtain the basis of characteristic quantity, but in general, partial discharge DISCHARGE PULSES EXTRACTION is built upon Partial discharge signal denoising base
On plinth, i.e., " first denoising, rear extraction ", therefore denoising result directly affects the validity of DISCHARGE PULSES EXTRACTION.
Application No. is CN201611094903.6, a kind of entitled " pulse suitable for high-frequency current Partial Discharge Detection
The Chinese invention patent of extracting method ", provide it is a kind of based on wavelet decomposition-reconstruct partial discharge denoising method, then using from
Dynamic threshold calculations extract pulse burst.Due to live Partial discharge signal complexity, only by the peak values of coefficient of wavelet decomposition at different levels with
Missing inspection, erroneous detection can inevitably occur for the method that the ratio between virtual value determines signal type.Partial discharge pulse is chosen by automatic threshold
Method is only applicable to the higher situation of signal-to-noise ratio, and when partial discharge pulse biggish there are amplitude difference, the lesser part of amplitude
Pulse is easier to be identified as noise.
Summary of the invention
Extract that validity is poor, is influenced etc. to ask by denoising result for the local pulse currently based on " first denoising, then extract "
Topic, the present invention provides a kind of signal processing of partial discharge methods for taking into account DISCHARGE PULSES EXTRACTION and signal denoising, can effectively extract
Local pulse, even if still there is stronger white noise rejection ability in low signal ratio.
Another object of the present invention is intended to provide a kind of power cable partial discharge positioning method.
The present invention is primarily based on random in short-term using the local discharge signal data processing policy of " first extracting, rear denoising "
Singular value decomposition and in short-term singular value decomposition are realized to the quick of pulse, accurate positioning and are extracted, and singular value is further utilized
Decomposition-reconstruct method is denoised, and local signal white noise is made to be effectively suppressed.
The DISCHARGE PULSES EXTRACTION and the signal processing of partial discharge method of signal denoising provided by the invention taken into account includes following step
It is rapid:
(1) Partial discharge signal is filtered, obtains Partial discharge signal sequence s (n), n=1,2,3 ..., N0, N0It is whole
A Partial discharge signal sequence length;
(2) it is calculated by the white noise standard deviation η of Partial discharge signal sequence s (n) and determines the unusual threshold value of Partial discharge signal global optimum
τ*;
It (3) is T using lengthWTime window intercept Partial discharge signal section y (n) from the partial signal sequence s (n), n=1,
2,3 ..., N, N are Partial discharge signal sequence length in time window;Matrix is constructed using the Partial discharge signal in the signal segment, and to structure
The matrix built carries out singular value decomposition, and the maximum singular value σ of Partial discharge signal in this section is calculatedmax;By σmaxWith step (2)
In the unusual threshold tau * of global optimum that is calculated compare, if σmax≥τ*, which is recorded as pulse and is searched
Rope initial time ti, i indicate extract umber of pulse, enter step (4);If σmax<τ*, time window is along the sliding of the direction of sequence ends
One step first sets step-length, intercepts local signal section again using time window, rebuilds matrix, calculates σmax, with τ*Compare,
Circulation is until σmax≥τ*;
(4) the Partial discharge signal data in time window are rebuild into matrix, and singular value point is carried out to the matrix of building
Solution, is calculated the unusual estimation threshold value of local optimum of Partial discharge signal in this sectionBy optimal unusual estimation threshold valueWith most
Big singular value δmaxCarry out following compare:
If1. time window utilizes time window weight along the sliding second setting step-length that moves a step in the direction that sequence originates
New interception local signal section, rebuilds matrix, calculates δmax, withCompare, circulation untilThe time window is risen
Moment beginning is recorded as the starting t of partial discharge DISCHARGE PULSES EXTRACTIONs;2. time window is back to pulse search initial time ti, and by the time
The sliding third setting step-length that moves a step in the direction of windowsill sequence ends, is intercepted local signal section again using time window, rebuild
Matrix calculates δmax, withCompare, circulation untilThe time window end time is recorded as partial discharge DISCHARGE PULSES EXTRACTION
Finish time te;Extract ts-teThe Partial discharge signal data at moment are i-th of partial discharge pulse signal;
IfThe time window initial time is recorded as to the starting t of partial discharge DISCHARGE PULSES EXTRACTIONs, and by time window edge
The sliding third setting step-length that moves a step in the direction of sequence ends, is intercepted local signal section again using time window, rebuilds square
Battle array calculates δmax, withCompare, circulation untilThe time window end time is recorded as partial discharge DISCHARGE PULSES EXTRACTION
Finish time te;Extract ts-teThe Partial discharge signal data at moment are i-th of partial discharge pulse signal;
(5) singular value decomposition-reconstruct is carried out to i-th of partial discharge pulse signal that step (4) is extracted, after obtaining denoising
I-th of partial discharge pulse signal;
(6) after i-th of partial discharge pulse signal after being denoised, return step (3)-(5) are searched for and are extracted next
Non- partial discharge pulse train is all set to zero, completion office until searching original Partial discharge signal end by partial discharge pulse signal
Put DISCHARGE PULSES EXTRACTION and signal denoising.
The above-mentioned signal processing of partial discharge method for taking into account DISCHARGE PULSES EXTRACTION and signal denoising in step (1), utilizes high pass
Filter is filtered Partial discharge signal, to filter out low-frequency interference signal.The cutoff frequency range of the high-pass filter
Are as follows: 50kHz~100kHz.The Partial discharge signal handled by high pass filter filters further reaches acquisition dress through coaxial cable
It sets to obtain Partial discharge signal sequence s (n), the sample rate of the acquisition device is at least 50MS/s, to guarantee that single partial discharge pulse has
More data point is convenient for subsequent analysis.
The above-mentioned signal processing of partial discharge method for taking into account DISCHARGE PULSES EXTRACTION and signal denoising determines noise by the following method
Standard deviation and the calculating unusual threshold value of global optimum:
(21) determine the white noise standard deviation η of Partial discharge signal sequence s (n): using wavelet transformation (such as db4, db6 or
Db8 small echo etc.) wavelet decomposition is carried out to Partial discharge signal sequence s (n), by first layer wavelet details coefficient w1,j1 × j matrix is constructed,
J=N0/ 2, withFor approximate evaluation white noise standard deviation, wherein median refers to matrix element w1,j
Absolute value take median;
(22) the unusual threshold tau of Partial discharge signal global optimum is calculated as follows*,In formulaβ=L0/K0, K0=N0-L0+1;L0Value range be usually N0/20
~N0/2。
The above-mentioned signal processing of partial discharge method for taking into account DISCHARGE PULSES EXTRACTION and signal denoising is small using db4, db6 or db8
Wave etc. carries out wavelet decomposition to Partial discharge signal sequence s (n).The present invention carries out several layers of wavelet decompositions to Partial discharge signal using small echo
There is no particular/special requirement, regardless of carrying out several layers of wavelet decompositions, the present invention only takes first layer wavelet details coefficient.
The above-mentioned signal processing of partial discharge method for taking into account DISCHARGE PULSES EXTRACTION and signal denoising determines maximum odd by the following method
Different value and partial discharge pulse search POS INT point:
It (31) is T using lengthWTime window intercept Partial discharge signal segment y (n) from the partial signal sequence s (n), n=
1,2,3 ..., N, N are Partial discharge signal sequence length in time window;
(32) Hankel matrix A is constructedL×K:
In formula, K=N-L+1, L=N/3;
(33) calculating matrix Y=AX, X are the random matrix X=randn (K, 2) of K × 2, and by Y orthonormal,
Obtain the rank orthogonal basis of L × 2 W1;
(34) calculating matrix B=W1 T·A;
(35) according to B=U Σ VTBy matrix B singular value decomposition, in formula, U and V are respectively that 2 × 2 peacekeeping K × K dimension is orthogonal
Matrix;It is diagonal matrix, diagonal element is that the singular value of matrix B arranges in descending order, σ1It is as maximum odd
Different value σmax;
(36) by σmaxWith the unusual threshold tau of global optimum being calculated in step (2)*Compare, if σmax≥τ*, when by this
Between window initial time be recorded as pulse search initial time ti, enter step (4);I indicates to extract umber of pulse;If σmax<τ*, when
Between windowsill sequence ends the sliding first setting step-length that moves a step in direction, intercept the return of local signal section again using time window and hold
Row (31)-(36) rebuilds matrix, calculates σmax, with τ*Compare, circulation is until σmax≥τ*。
The above-mentioned signal processing of partial discharge method for taking into account DISCHARGE PULSES EXTRACTION and signal denoising determines i-th by the following method
Partial discharge pulse signal:
(41) the Partial discharge signal data in time window are rebuild into Hankel matrix Z:
In formula, K=N-L+1, L=N/3;
(42) according to Z=C Σ ' DTThe Hankel matrix Z of step (41) building is subjected to singular value decomposition, Z=C Σ ' DT,
C and D is respectively L × L peacekeeping K × K dimension orthogonal matrix in formula;Σ '=diag (δ1,δ2,…,δp) (p=min (L, K)) be pair
Angular moment battle array, diagonal element are that the singular value of matrix Z arranges in descending order;Amount of orientation δ=(δ1,δ2,...,δm), m is non-in Σ '
The number and m≤p of neutral element are the order of matrix Z;δ1For maximum singular value δmax;
(43) withFor the unusual estimation threshold value of local optimum of data in data frame, ω in formula
(β)≈0.56β3-0.95β2+ 1.82 β+1.43, β=L/K, median, which refers to, takes median to the absolute value of matrix delta;
(44) by optimal unusual estimation threshold valueWith maximum singular value δmaxCompare,
If1. time window utilizes time window weight along the sliding second setting step-length that moves a step in the direction that sequence originates
New interception local signal section, which returns, executes (41)-(44), rebuilds matrix, calculates δmax, withCompare, circulation untilThe time window initial time is recorded as to the starting t of partial discharge DISCHARGE PULSES EXTRACTIONs;2. time window is back to pulse to search
Rope initial time ti, and time window is cut along the sliding third setting step-length that moves a step in the direction of sequence ends using time window again
It takes local signal section to return and executes (41)-(44), rebuild matrix, calculate δmax, withCompare, circulation untilThe time window end time is recorded as to the finish time t of partial discharge DISCHARGE PULSES EXTRACTIONe;Extract ts-teThe partial discharge at moment
Signal data is i-th of partial discharge pulse signal;
IfThe time window initial time is recorded as to the starting t of partial discharge DISCHARGE PULSES EXTRACTIONs, and by time window edge
The sliding third setting step-length that moves a step in the direction of sequence ends, is intercepted local signal section again using time window and returns to execution (41)-
(44), matrix is rebuild, δ is calculatedmax, withCompare, circulation untilThe time window end time is recorded as
The finish time t of partial discharge DISCHARGE PULSES EXTRACTIONe;Extract ts-teThe Partial discharge signal data at moment are i-th of partial discharge pulse signal.
The above-mentioned signal processing of partial discharge method for taking into account DISCHARGE PULSES EXTRACTION and signal denoising, carries out singular value by the following method
Decomposition-reconstruct:
(51) according to ts-teI-th of partial discharge pulse data signal at moment constructs Hankel matrix Z':
In formula, K'=N'-L'+1, L'=N'/3;N' is ts-teThe Partial discharge signal sequence length at moment;
(52) according to Z '=C ' Σ " D 'TThe Hankel matrix Z' progress singular value decomposition that step (51) are constructed, Z '=
C′Σ″D′T, C' and D' is respectively L' × L' peacekeeping K' × K' dimension orthogonal matrix in formula;Σ "=diag (δ '1,δ'2,…, δ'p')
(p'=min (L', K')) is diagonal matrix, and diagonal element is that the singular value of matrix Z' arranges in descending order;Amount of orientation δ '=
(δ′1,δ′2,...,δ′m′), the number and m'≤p' of nonzero element in m' Σ " are the order of matrix Z';
(53) withFor the unusual estimation threshold value of local optimum of data in data frame, ω in formula
(β′)≈0.56β′3-0.95β′2+ 1.82 β '+1.43, β '=L'/K', median refer to matrix delta ' absolute value take middle position
Number;
(54) by singular value matrix Σ "=diag (δ '1,δ'2,…,δ'p') in be lower thanSingular value be set as zero, obtain
Singular value matrix Σ " ' after to denoising;According to Z "=C ' Σ " ' D 'TThe Hankel matrix Z " of signal after being denoised;
(55) for the matrix Z " obtained after reconstruct, to t in the way of being averageds-teThe Partial discharge signal at moment carries out
Reconstruct
I-th of partial discharge pulse signal after being denoised in turn and its position in entire original signal sequence.
The above-mentioned signal processing of partial discharge method for taking into account DISCHARGE PULSES EXTRACTION and signal denoising, Partial discharge signal extracts in order to balance
Validity and extraction efficiency, T in step (31)WTake 1~2 μ s, first setting N/4~N/2 data of step-length in step (36)
Point, step (42) second set step-length as 2~5 data points, and third sets step-length as 2~5 data points.
The above-mentioned signal processing of partial discharge method for taking into account DISCHARGE PULSES EXTRACTION and signal denoising, the partial discharge pulse signal extracted
Be primarily present following three kinds of situations: (1) the partial discharge pulse signal that white noise is effectively suppressed, (2) there are the obvious higher-orders of oscillation
The partial discharge pulse signal of noise, (3) denoising result show as the noise signal of the higher-order of oscillation.In order to further suppress partial discharge arteries and veins
It rushes the higher-order of oscillation noise in signal, after step (6), 1 layer scattering is further utilized to extracted partial discharge pulse signal
Wavelet transformation carries out threshold denoising to detail coefficients, which is the conventional means that this field has disclosed, referring to Shen
It please be number for partial discharge signal denoising method and Zhang Bo based on wavelet adaptive threshold disclosed in CN201310473326.1
Deng disclosed Wavelet Denoising Method to research (Zhang Bo, Liu Chengguo, Xu Zhong, the soughing of the wind in forest trees for improving GIS ultrasound local discharge signal discrimination
Research [J] the electrical engineering journal of Wavelet Denoising Method to raising GIS ultrasound local discharge signal discrimination, 2017,12 (11):
41-45).After singular value decomposition-discrete wavelet joint denoising, the high fdrequency component in signal extracted effectively is pressed down
System;In addition, for the noise signal for being mistaken for partial discharge pulse signal, since high fdrequency component has obtained obvious inhibition, after denoising its
Local optimum singular value threshold value when maximum singular value will be significantly less than denoising, so reject only include high-frequency noise signal
Sequence;To the partial discharge pulse containing high-frequency noise, high-frequency noise is further suppressed;It is most of for partial discharge pulse signal
Energy is retained, and the maximum singular value variation of denoising front and back less, is achieved in and the further of partial discharge pulse is extracted and gone
It makes an uproar.
Invention further provides a kind of power cable partial discharge positioning methods, at the local discharge signal
Reason method extracts local pulse signal, extract two partial discharge pulse signal peak values for being no more than setting time interval and to it is corresponding when
Carve t1,t2, foundationDetermine position of the partial discharge source on power cable, d is that partial discharge source is tested apart from power cable
The distance at end, l are power cable total length, and v is pulse propagation velocity.
Compared with prior art, the signal processing of partial discharge provided by the invention for taking into account DISCHARGE PULSES EXTRACTION and signal denoising
Method has the advantages that following very prominent and advantageous effects:
1, the present invention is to Partial discharge signal, using the data processing order of " first extracting, rear to denoise ", first with singular value point
Solution, optimal unusual threshold value is compared with maximum singular value, realizes quick, accurate positioning to partial discharge pulse, and to partial discharge
Pulse extracts, and recycles singular value reconstruct further to denoise the local pulse signal of extraction, effectively avoids low amplitude value office
Pulse signal is put directly to be filtered out.
2, the present invention is in terms of the DISCHARGE PULSES EXTRACTION to Partial discharge signal, based on random singular value decomposition and global optimum in short-term
Unusual threshold value realizes the quick positioning of partial discharge pulse signal search starting point;Simultaneously using singular value decomposition in short-term and based on the time
The method of the unusual threshold estimation of the local optimum of data carries out porch search, it can be achieved that partial discharge pulse signal originates in window
With the accurate positioning of final position, effective extraction to partial discharge pulse signal is completed.
3, the present invention passes through the singular value decomposition-of the unusual threshold estimation of local optimum in terms of the denoising to Partial discharge signal
Reconstruct, even if being still able to achieve effective inhibition of Partial discharge signal white noise under the conditions of compared with low signal-to-noise ratio.
4, the present invention also has good expansion, combines 1 layer scattering wavelet transformation, carries out threshold value to detail coefficients and goes
It makes an uproar, the High-frequency Interference in local signal after removal white noise can be further suppressed.
5, the present invention is by may be implemented electricity to the precise positioning of partial discharge pulse signal on power cable and effectively extraction
The accurate positioning in power cable partial discharge source (i.e. shelf depreciation point), to ensure that power cable operation provides reliable basis.
Detailed description of the invention
Fig. 1 is the signal processing of partial discharge method flow schematic diagram that the present invention takes into account DISCHARGE PULSES EXTRACTION and signal denoising.
Fig. 2 is the Partial discharge signal drawn in the embodiment of the present invention 1;Wherein (a) is original Partial discharge signal under power frequency, (b) is
Partial discharge signal after denoising, (c) the partial discharge pulse signal to extract.
Fig. 3 is the Partial discharge signal drawn in the embodiment of the present invention 2;Wherein (a) is original Partial discharge signal under oscillation wave, (b)
For the Partial discharge signal after denoising.
Fig. 4 is that the partial discharge source drawn in embodiment 2 positions schematic diagram.
Specific embodiment
Provide the embodiment of the present invention below with reference to attached drawing, and by embodiment to technical solution of the present invention carry out into
Clear, the complete explanation of one step.Obviously, the embodiment is only a part of the embodiments of the present invention, rather than whole realities
Apply example.Based on the content of present invention, those of ordinary skill in the art obtained institute without making creative work
There are other embodiments, belongs to the range that the present invention is protected.
Embodiment 1
The research object that the present embodiment is directed to is the cold of preset longitudinal tool marks defect (long 100mm, width 0.2mm, depth 1mm)
Contracting cable termination, cable model YJV22-8.7/15 will design defective cable and test partial discharge under 15kV power-frequency voltage
Signal, and carry out DISCHARGE PULSES EXTRACTION and denoising.
The signal processing of partial discharge method provided in this embodiment for taking into account DISCHARGE PULSES EXTRACTION and signal denoising, such as Fig. 1 institute
Show, comprising the following steps:
S1 is filtered Partial discharge signal, obtains Partial discharge signal sequence s (n), n=1,2,3 ..., N0, N0It is entire
Partial discharge signal sequence length.
The present embodiment uses High Frequency Current Sensor to obtain the Partial discharge signal on cable, the Partial discharge signal warp of acquisition first
High-pass filter is exported to acquisition device.The cutoff frequency range for the high-pass filter that the present embodiment uses are as follows: 50kHz~
100kHz.The acquisition device model Tektronix TDS7104 that the present embodiment uses, sample rate are set as 50MS/s.From adopting
The original Partial discharge signal sequence s (n) of 2ms duration is chosen in the original Partial discharge signal sequence that acquisition means obtain (such as Fig. 2 (a) institute
Show) carry out partial discharge DISCHARGE PULSES EXTRACTION and denoising, n=1,2,3 ..., N0, N0=105For entire Partial discharge signal sequence length.
S2 is calculated by the white noise standard deviation η of Partial discharge signal sequence s (n) and is determined the unusual threshold value of Partial discharge signal global optimum
τ*:
S21 determines the white noise standard deviation of Partial discharge signal sequence s (n): being played a game using small echo (morther wavelet selection ' db8 ')
Discharge signal sequence s (n) carries out wavelet decomposition, takes first layer wavelet details coefficient w1,j, and utilize first layer wavelet details coefficient
Construct 1 × j matrix, j=N0/ 2, withFor approximate evaluation white noise standard deviation, wherein median refers to
To matrix element w1,jAbsolute value take median;
S22 calculates the unusual threshold tau of Partial discharge signal global optimum as follows*,In formulaβ=L0/K0, K0=N0-L0+1;L0=N0/ 10, N0=105。
S3 obtains maximum singular value using random singular value decomposition, determines partial discharge pulse search POS INT point:
S31 is T using lengthWThe time window of (1.5 μ s) intercepts Partial discharge signal segment y from partial signal sequence s (n)
(n), n=1,2,3 ..., N, N=76 are Partial discharge signal sequence length in time window;
S32 constructs Hankel matrix AL×K:
In formula, L=25, K=N-L+1=52;
S33 calculating matrix Y=AX, X are the rank random matrix of K × 2 X=constructed using randn function in matlab
Randn (52,2), and matrix Y orthonormal is obtained into the rank orthogonal basis of L × 2 W using orth function in matlab1;
S34 calculating matrix B=W1 T·A;
S35 is according to B=U Σ VTBy matrix B singular value decomposition, in formula, U and V are respectively that 2 × 2 peacekeepings 52 × 52 dimension is orthogonal
Matrix;It is diagonal matrix, diagonal element is that the singular value of matrix B arranges in descending order, σ1It is as maximum odd
Different value σmax;
S36 is by σmaxWith the unusual threshold tau of global optimum being calculated in step S2*Compare, if σmax≥τ*, by the time
Window initial time is recorded as pulse search initial time ti, enter step S4;I indicates to extract umber of pulse, every to have extracted an arteries and veins
Signal is rushed, umber of pulse is extracted and is increased by 1, carry out the extraction of next pulse signal;If σmax<τ*, time window is along sequence ends
Direction cunning move a step the first setting step-length (38 data points), intercepted again using time window local signal section return execute
S31-S36 rebuilds matrix, calculates σmax, compared with τ *, circulation is until σmax≥τ*。
S4 executes side search, obtains i-th of partial discharge pulse signal:
Partial discharge signal data in time window are constructed Hankel matrix Z by S41:
In formula, L=25, K=N-L+1=52;
S42 is according to Z=C Σ ' DTThe Hankel matrix Z that step S41 is constructed carries out singular value decomposition, Z=C Σ ' DT, formula
Middle C and D is respectively 25 × 25 peacekeepings 52 × 52 dimension orthogonal matrix;Σ '=diag (δ1,δ2,…, δp) (p=25) be to angular moment
Battle array, diagonal element are that the singular value of matrix Z arranges in descending order;Amount of orientation δ=(δ1,δ2,...,δm), m is non-zero entry in Σ '
The number and m≤p of element are the order of matrix Z;δ1For maximum singular value δmax;
S43 withFor the unusual estimation threshold value of local optimum of data in data frame, ω in formula
(β)≈0.56β3-0.95β2+ 1.82 β+1.43, β=L/K, median, which refers to, takes median to the absolute value of matrix delta;
S44 is by optimal unusual estimation threshold valueWith maximum singular value δmaxCompare:
IfNeed to complete following two step: 1. time window is along sliding second setting that moves a step in the direction that sequence originates
Step-length (3 data points) is intercepted local signal section again using time window and returns to execution S41-S44, rebuild matrix, count
Calculate δmax, withCompare, circulation untilThe time window initial time is recorded as to the starting t of partial discharge DISCHARGE PULSES EXTRACTIONs;
2. time window is back to pulse search initial time ti, and time window is set along the sliding third that moves a step in the direction of sequence ends
Fixed step size (3 data points) is intercepted local signal section again using time window and returns to execution S41-S44, rebuilds matrix,
Calculate δmax, withCompare, circulation untilThe time window end time is recorded as to the end of partial discharge DISCHARGE PULSES EXTRACTION
Moment te;Extract ts-teThe Partial discharge signal data at moment are i-th partial discharge pulse signal;
IfTime window default is searched for since the initial position of Partial discharge signal sequence, which is originated
Moment is recorded as the starting t of partial discharge DISCHARGE PULSES EXTRACTIONs, and by time window along the sliding third setting step that moves a step in the direction of sequence ends
Long (3 data points) is intercepted local signal section again using time window and returns to execution S41-S44, rebuild matrix, calculate
δmax, withCompare, circulation untilThe time window end time is recorded as to the finish time of partial discharge DISCHARGE PULSES EXTRACTION
te;Extract ts-teThe Partial discharge signal data at moment are i-th of partial discharge pulse signal.
S5 carries out singular value decomposition-reconstruct to i-th of partial discharge pulse signal that step S4 is extracted, after being denoised
I-th of partial discharge pulse signal:
S51 is according to ts-teI-th of partial discharge pulse data signal at moment constructs Hankel matrix Z':
In formula, K'=N'-L'+1, L'=N'/3;N' is ts-teThe Partial discharge signal sequence length at moment;
S52 is according to Z '=C ' Σ " D 'TThe Hankel matrix Z' that step S51 is constructed carries out singular value decomposition, Z '=C '
Σ″D′T, C' and D' is respectively L' × L' peacekeeping K' × K' dimension orthogonal matrix in formula;Σ "=diag (δ '1,δ'2,…, δ'p')
(p'=min (L', K')) is diagonal matrix, and diagonal element is that the singular value of matrix Z' arranges in descending order;Amount of orientation δ '=
(δ′1,δ′2,...,δ′m′), the number and m'≤p' of nonzero element in m' Σ " are the order of matrix Z';
S53 withFor the unusual estimation threshold value of local optimum of data in data frame, ω in formula
(β′)≈0.56β′3-0.95β′2+ 1.82 β '+1.43, β '=L'/K', median refer to matrix delta ' absolute value take middle position
Number;
S54 is by singular value matrix Σ "=diag (δ '1,δ'2,…,δ'p') in be lower thanSingular value be set as zero, obtain
Singular value matrix Σ " ' after denoising;The Hankel matrix Z " of signal after being denoised according to Z "=C ' Σ " ' D ' T;
S55 is for obtained matrix Z " after reconstruct, to t in the way of being averageds-teThe Partial discharge signal at moment carries out weight
Structure
I-th of partial discharge pulse signal after being denoised in turn and its position in entire original signal sequence are such as schemed
Shown in 2 (c).
S6 denoised after i-th of partial discharge pulse signal after return step S3-S5 search for and extract next partial discharge
Pulse (i.e. i+1 partial discharge pulse signal), until original Partial discharge signal end is searched, non-partial discharge pulse train is whole
It is set as zero, completes partial discharge DISCHARGE PULSES EXTRACTION and signal denoising.
S7 carries out a layer scattering Wavelet Denoising Method to the partial discharge pulse signal of said extracted, filters out High-frequency Interference, finally obtains
Partial discharge signal such as Fig. 2 (b) shown in, it can be seen from the figure that the white noise in original Partial discharge signal has obtained effective inhibition,
Low amplitude value partial discharge pulse signal also accurately is extracted out simultaneously.
Embodiment 2
The research object that the present embodiment is directed to is length 498m 10kV XLPE cable (away from test lead 249m transition joint
Locate default semi-conductive layer and overlap bad defect).Defective cable will be designed and test Partial discharge signal under oscillation wave, and carried out
Defect location is realized in DISCHARGE PULSES EXTRACTION and denoising.
The signal processing of partial discharge method provided in this embodiment for taking into account DISCHARGE PULSES EXTRACTION and signal denoising, such as Fig. 1 institute
Show, comprising the following steps:
S1 is filtered Partial discharge signal, obtains Partial discharge signal sequence s (n), n=1,2,3 ..., N0, N0It is entire
Partial discharge signal sequence length.
The present embodiment is first under oscillation wave condition (ceiling voltage is no more than 20kV, frequency 67Hz), using high-frequency electrical
Flow sensor obtains the Partial discharge signal on cable, and the Partial discharge signal of acquisition is exported through high-pass filter to acquisition device.This implementation
The cutoff frequency range for the high-pass filter that example uses are as follows: 50kHz~100kHz.The acquisition device model that the present embodiment uses
For Tektronix TDS7104, sample rate 50MS/s.6ms is chosen from the original Partial discharge signal sequence that acquisition device obtains
Original Partial discharge signal sequence s (n) (shown in such as Fig. 3 (a)) the progress partial discharge DISCHARGE PULSES EXTRACTION of duration and denoising, n=1,2,
3,…,N0, N0=3 × 105For entire Partial discharge signal sequence length.
S2 is calculated by the white noise standard deviation η of Partial discharge signal sequence s (n) and is determined the unusual threshold value of Partial discharge signal global optimum
τ*:
S21 determines the white noise standard deviation of Partial discharge signal sequence s (n): being played a game using small echo (morther wavelet selection ' db8 ')
Discharge signal sequence s (n) carries out wavelet decomposition, takes first layer wavelet details coefficient w1,j, and utilize first layer wavelet details coefficient
Construct 1 × j matrix, j=N0/ 2, withFor approximate evaluation white noise standard deviation, wherein median refers to
To matrix element w1,jAbsolute value take median;
S22 calculates the unusual threshold tau of Partial discharge signal global optimum as follows*,In formulaβ=L0/K0, K0=N0-L0+1;L0=N0/ 10, N0=3 × 105。
S3 obtains maximum singular value using random singular value decomposition, determines partial discharge pulse search POS INT point:
S31 is T using lengthWThe time window of (1.5 μ s) intercepts Partial discharge signal segment y from partial signal sequence s (n)
(n), n=1,2,3 ..., N, N=76 are Partial discharge signal sequence length in time window;
S32 constructs Hankel matrix AL×K:
In formula, L=25, K=N-L+1=52;
S33 calculating matrix Y=AX, X are the rank random matrix of K × 2 X=constructed using randn function in matlab
Randn (52,2), and matrix Y orthonormal is obtained into the rank orthogonal basis of L × 2 W using orth function in matlab1;
S34 calculating matrix B=W1 T·A;
S35 is according to B=U Σ VTBy matrix B singular value decomposition, in formula, U and V are respectively that 2 × 2 peacekeepings 52 × 52 dimension is orthogonal
Matrix;It is diagonal matrix, diagonal element is that the singular value of matrix B arranges in descending order, σ1It is as maximum odd
Different value σmax;
S36 is by σmaxWith the unusual threshold tau of global optimum being calculated in step S2*Compare, if σmax≥τ*, by the time
Window initial time is recorded as pulse search initial time ti, enter step S4;I indicates to extract umber of pulse, every to have extracted an arteries and veins
Signal is rushed, umber of pulse is extracted and is increased by 1, carry out the extraction of next pulse signal;If σmax<τ*, time window is along sequence ends
Direction cunning move a step the first setting step-length (38 data points), intercepted again using time window local signal section return execute
S31-S36 rebuilds matrix, calculates σmax, with τ*Compare, circulation is until σmax≥τ*。
S4 executes side search, obtains i-th of partial discharge pulse signal:
Partial discharge signal data in time window are constructed Hankel matrix Z by S41:
In formula, L=25, K=N-L+1=52;
S42 is according to Z=C Σ ' DTThe Hankel matrix Z that step S41 is constructed carries out singular value decomposition, Z=C Σ ' DT, formula
Middle C and D is respectively 25 × 25 peacekeepings 52 × 52 dimension orthogonal matrix;Σ '=diag (δ1,δ2,…, δp) (p=25) be to angular moment
Battle array, diagonal element are that the singular value of matrix Z arranges in descending order;Amount of orientation δ=(δ1,δ2,...,δm), m is non-zero entry in Σ '
The number and m≤p of element are the order of matrix Z;δ1For maximum singular value δmax;
S43 withFor the unusual estimation threshold value of local optimum of data in data frame, ω in formula
(β)≈0.56β3-0.95β2+ 1.82 β+1.43, β=L/K, median, which refers to, takes median to the absolute value of matrix delta;
S44 is by optimal unusual estimation threshold valueWith maximum singular value δmaxCompare:
IfNeed to complete following two step: 1. time window is along sliding second setting that moves a step in the direction that sequence originates
Step-length (3 data points) is intercepted local signal section again using time window and returns to execution S41-S44, rebuild matrix, count
Calculate δmax, withCompare, circulation untilThe time window initial time is recorded as to the starting t of partial discharge DISCHARGE PULSES EXTRACTIONs;
2. time window is back to pulse search initial time ti, and time window is set along the sliding third that moves a step in the direction of sequence ends
Fixed step size (3 data points) is intercepted local signal section again using time window and returns to execution S41-S44, rebuilds matrix,
Calculate δmax, withCompare, circulation untilThe time window end time is recorded as to the end of partial discharge DISCHARGE PULSES EXTRACTION
Moment te;Extract ts-teThe Partial discharge signal data at moment are i-th partial discharge pulse signal;
IfTime window default is searched for since the initial position of Partial discharge signal sequence, which is risen
Moment beginning is recorded as the starting t of partial discharge DISCHARGE PULSES EXTRACTIONs, and by time window along the sliding third setting that moves a step in the direction of sequence ends
Step-length (3 data points) is intercepted local signal section again using time window and returns to execution S41-S44, rebuild matrix, count
Calculate δmax, withCompare, circulation untilAt the end of the time window end time is recorded as partial discharge DISCHARGE PULSES EXTRACTION
Carve te;Extract ts-teThe Partial discharge signal data at moment are i-th of partial discharge pulse signal.
S5 carries out singular value decomposition-reconstruct to i-th of partial discharge pulse signal that step S4 is extracted, after being denoised
I-th of partial discharge pulse signal:
S51 is according to ts-teI-th of partial discharge pulse data signal at moment constructs Hankel matrix Z':
In formula, K'=N'-L'+1, L'=N'/3;N' is ts-teThe Partial discharge signal sequence length at moment;
S52 is according to Z '=C ' Σ " D 'TThe Hankel matrix Z' that step S51 is constructed carries out singular value decomposition, Z '=C '
Σ″D′T, C' and D' is respectively L' × L' peacekeeping K' × K' dimension orthogonal matrix in formula;Σ "=diag (δ '1,δ'2,…, δ'p')
(p'=min (L', K')) is diagonal matrix, and diagonal element is that the singular value of matrix Z' arranges in descending order;Amount of orientation δ '=
(δ′1,δ′2,...,δ′m′), the number and m'≤p' of nonzero element in m' Σ " are the order of matrix Z';
S53 withFor the unusual estimation threshold value of local optimum of data in data frame, ω in formula
(β′)≈0.56β′3-0.95β′2+ 1.82 β '+1.43, β '=L'/K', median refer to matrix delta ' absolute value take middle position
Number;
S54 is by singular value matrix Σ "=diag (δ '1,δ'2,…,δ'p') in be lower thanSingular value be set as zero, obtain
Singular value matrix Σ " ' after denoising;According to Z "=C ' Σ " ' D 'TThe Hankel matrix Z " of signal after being denoised;
S55 is for obtained matrix Z " after reconstruct, to t in the way of being averageds-teThe Partial discharge signal at moment carries out weight
Structure
I-th of partial discharge pulse signal after being denoised in turn and its position in entire original signal sequence.
S6 denoised after i-th of partial discharge pulse signal after return step S3-S5 search for and extract next partial discharge
Pulse (i.e. i+1 partial discharge pulse signal), until original Partial discharge signal end is searched, non-partial discharge pulse train is whole
It is set as zero, completes partial discharge DISCHARGE PULSES EXTRACTION and signal denoising.
S7 carries out a layer scattering Wavelet Denoising Method to the pulse signal of putting of said extracted, filters out High-frequency Interference, finally obtained
Shown in Partial discharge signal such as Fig. 3 (b), extracted partial discharge pulse signal enlarged drawing is corresponded to for dashed box in Partial discharge signal figure below.
It can be seen from the figure that the white noise in original Partial discharge signal has obtained effective inhibition, while low amplitude value partial discharge pulse signal
Accurately it is extracted out.
To the Partial discharge signal after denoising obtained in step S7, extraction time interval is no more than two partial discharge pulses of 6 μ s
Signal peak and corresponding moment t1,t2.By calibration, pulse propagation velocity is v=163.23m/ μ s, is calculate by the following formula partial discharge
Source is away from test lead distanceL is total cable length.Repeatedly being calculated can obtain after being averaged, partial discharge source be located at away from
From (as shown in Figure 4) at test lead 247.17m, it coincide substantially with initial preset defective locations (away from test lead 249m).Therefore,
The signal processing of partial discharge method provided through the invention can be accurately positioned and effectively extract partial discharge pulse signal, in turn
The accurate position for obtaining partial discharge source.