The content of the invention
The object of the present invention is to make an uproar to effectively suppress the signal obtained in high voltage electric equipment partial discharge monitoring
Sound, the invention discloses a kind of Weighted Threshold wavelet de-noising method for partial discharge on-line monitoring.
The technical scheme is that the present invention, which is first according to similarity criterion, selectes a kind of mother wavelet, to including noise
Local discharge signal carry out n-layer wavelet decomposition;Then to signal noise variances sigma2Estimated, adopted with reference to coefficient of wavelet decomposition
Estimate threshold value with corresponding threshold strategies, the comprehensive uniform threshold of each layer coefficients selection, the Birge-Massart thresholds obtained to decomposition
Value, penalty threshold values and the weighting of multi thresholds strategy, to coefficient acting threshold process;The coefficient after processing is finally passed through into small echo
Rebuild and recover original signal, that is, obtain the time domain plethysmographic signal after denoising.
Denoising method of the present invention comprises the following steps:
1) a local discharge signal f (n) is s (n) after noise pollution, its basic noise model can represent
For:
S (n)=f (n)+σ e (n) (1)
Wherein e (n) is noise, and σ is noise intensity, defines noise w (n)=σ e (n).In the simplest situations can be false
If e (n) is white Gaussian noise, σ=1.
To local discharge signal, the related coefficient γ that is determined with formula (2) assesses discharge pulse signal and small echo atom
Similitude.
In formula, xiFor partial discharge pulse's signal data sequence,For xiAverage;yiFor the data sequence of small echo atom,For yiAverage.The similitude of the bigger expression discharge signals of γ and selected small echo is higher.
Local discharge signal is precipitous with rising edge, and low amplitude value, the duration is short, and decay fast characteristic.With shelf depreciation
The characteristic of signal is compared, and there is higher similitude in compact schemes orthogonal wavelet Daubechies (db) wavelet functions race.In db small echos
In family of functions, db5 small echo supports length is 9, vanishing moment 5 so that the energy after conversion is more concentrated, and is provided simultaneously with preferable side
Boundary's characteristic.Decomposition order is also larger on noise reduction influence, and decomposing thinner low frequency part will preferably be suppressed, but be increased
Decomposition order needs more processing times, considers factors above, selects 5 layers of db5 wavelet decompositions to local discharge signal
Carry out noise reduction.
2) uniform threshold is estimated.Donoho and Johnstone proposes the concept of uniform threshold, and proves in all pairs
In angular estimation, factor 2logeN is optimal, and uniform threshold calculation formula is as follows:
In formula, σ2For the noise variance (σ is standard deviation) estimated from signal, N is wavelet details coefficient lengths at different levels, by
The formula calculates gradient thresholds at different levels.
It is zero that if M, which is P average, variance isIndependent Gaussian stochastic variable absolute value middle position, then can prove:
Ε{M}≈0.6745σ0 (4)
The M mathematic expectaions that Ε { M } is, by ignoring partial discharge pulse's signal influence of itself, using detail coefficients
Middle position M estimates the variance of noise W:
Bringing formula (5) into formula (3) can obtain being calculated as below the equivalence formula of uniform threshold:
3) Birge-Massart threshold estimations.Birge-Massart threshold strategies include two empirical coefficients M and α, meter
It is as follows to calculate rule:
A Decomposition order j specified is given, to j+1 and higher, all coefficients retain.To i-th layer (1≤i≤
J), the n of maximum absolute value is retainediA coefficient, niDetermined by following formula:
ni=M (j+2-i)α (7)
M and α is empirical coefficient in formula, and M=L (1), wavelet coefficient after namely the 1st layer decomposition of L (1) are taken under default condition
Length, under normal circumstances, M should meet L (1)≤M≤2L (1);α takes α=3 in the case of noise reduction.
4) penalty threshold estimations.Penalty threshold values are derived from from Donoho-Johnstone methods, similar to formula
(6) uniform threshold provided, factor log2N is replaced, and computational methods see below formula:
The definition of M and N is identical with formula (6) in formula.
5) multi thresholds strategy weights, and obtains improved gradient threshold Tw(j):
In formula, i identifies for threshold strategies, and j is wavelet systems several levels, and for lev to divide Decomposition order, N is total threshold strategies number,
Ti(j) the j-th stage threshold value determined for i-th kind of threshold strategies, wi(j) it is corresponding weight.To definite series j=j0For, power
Value coefficient wi(j0) meet following constraint formula:
The present invention is weighted using 5 grades of wavelet decompositions and 3 kinds of threshold strategies, therefore lev=5, N=3.
6) coefficient after processing is recovered into original signal by wavelet reconstruction, that is, obtains the time domain plethysmographic signal after denoising.
The invention has the advantages that uniform threshold de-noising signal is too smooth, some letters of signal in itself are lost
Breath, does not meet similarity criterion.Burr is more after penalty threshold deniosings, however it remains obvious noise signal.By table 1
Threshold values at different levels can be seen that uniform threshold is larger, eliminate part useful signal, add denoising risk.Penalty threshold values are firm
Well in contrast, threshold value is smaller, and noise reduction is undesirable.Birge-Massart threshold deniosing effects are preferable, but to signal amplitude
Have a certain impact.Utilize Tw(j) noise reduction is carried out to original signal, it is found that Weighted Threshold significantly improves Birge-Massart
The signal amplitude attenuation problem that threshold band is come, while ensure that preferable signal local characteristic.
The present invention is suitable for local inside the high voltage electric equipments such as high-tension switch cabinet, transformer, power cable, GIS device
The denoising of discharge signal.
Embodiment:
The present embodiment scene local discharge signal comes from certain 500kV transformer online monitoring system.As shown in Figure 1, this reality
Apply example and line monitor signal denoising is carried out by Fig. 1 flows to local discharge signal.
As shown in Figure 2.Analyze field data to find, the present embodiment local discharge signal includes substantial amounts of interference signal, has
Even be submerged in interference.
Yn is the shelf depreciation electric discharge original signal of a power frequency period in Fig. 2, and other is coefficient of wavelet decomposition, wherein A5
For the 5th grade of approximation coefficient, Di (1≤i≤5) is i-stage detail coefficients.When D2 shows wavelet decomposition to the second layer, interference is
Substantially it is stripped.
The threshold value of 3 kinds of threshold value (uniform threshold, Birge-Massart threshold values, penalty threshold values) policy calculations more than
As shown in table 1, Fig. 3 shows the noise reduction of each threshold strategies, and wherein yn is original signal, and yd1~yd4 is de-noising signal.
Uniform threshold noise reduction is more satisfactory, but signal local feature is lost seriously;Penalty threshold values are generally less than normal, after noise reduction still
Disturbed in the presence of part;Birge-Massart threshold values are larger to the less impulse attenuation of amplitude, and part positive-negative polarity impulse attenuation is extremely
Unipolarity (preceding 5 pulses of yd3 in Fig. 3).
1 on-site signal of table is layered wavelet threshold
According to the characteristics of each threshold strategies, Weighted Threshold T is calculated by formula (10)w(j) it is shown in Table 1.Utilize Tw(j) to original
Beginning signal carries out noise reduction, and obtained de-noising signal is as shown in yd4 in Fig. 3.Compare preceding 5 pulses of yd3 and yd4 signals, can be with
It was found that Weighted Threshold significantly improves the signal amplitude attenuation problem that Birge-Massart threshold bands are come, while ensure that preferably
Signal local characteristic.