CN107831412A - Local discharge of electrical equipment localization method and system, storage medium and terminal device - Google Patents
Local discharge of electrical equipment localization method and system, storage medium and terminal device Download PDFInfo
- Publication number
- CN107831412A CN107831412A CN201710942874.2A CN201710942874A CN107831412A CN 107831412 A CN107831412 A CN 107831412A CN 201710942874 A CN201710942874 A CN 201710942874A CN 107831412 A CN107831412 A CN 107831412A
- Authority
- CN
- China
- Prior art keywords
- kurtosis
- local discharge
- frequency
- harmonics
- electrical equipment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
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/1209—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 using acoustic measurements
-
- 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/08—Locating faults in cables, transmission lines, or networks
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Testing Relating To Insulation (AREA)
- Filters And Equalizers (AREA)
Abstract
The present invention provides a kind of harmonics cut-off frequency of the power frequency of the sound pressure signal of local discharge of electrical equipment localization method and system, storage medium and terminal device including obtaining electrical equipment；Using the more stopband comb filter of linear phase shift to being that the sound pressure signal of electrical equipment is filtered, to obtain the residual signals obtained after the harmonics for filtering out power frequency；The kurtosis of the residual signals is calculated, and sets kurtosis threshold value；If kurtosis is more than the kurtosis threshold value, within the period corresponding to the kurtosis shelf depreciation occurs for judgement；The local discharge signal of phase property is included based on the shelf depreciation period judged and the generation of local discharge signal model；Realize that shelf depreciation positions according to the phase property of local discharge signal.The local discharge of electrical equipment localization method and system of the present invention, storage medium and terminal device can realize the accurate fast positioning of transformer station's shelf depreciation.
Description
Technical field
The present invention relates to the technical field of signal transacting, more particularly to a kind of local discharge of electrical equipment localization method and
System, storage medium and terminal device.
Background technology
Either transmission line of electricity due to contamination or corrosion caused by manufacturing process or longtime running, can cause electric field for transformer station
The change of intensity and shelf depreciation occurs.Phenomena such as shelf depreciation is with acousto-optic heat and chemical reaction, can produce electric flux
The aging of simultaneously acceleration components is lost, influences the stability of power supply unit work.Therefore, partial discharge monitoring turns into electrical equipment
The important component of on-line monitoring.However, because local discharge of electrical equipment noise has instantaneity, and in failure early stage
Caused shelf depreciation noise signal is more faint, while test site is there is substantial amounts of harmonics noise and other interference,
The accuracy of noise source monitoring and local breakdown location is had a strong impact on.
Shelf depreciation can determine its orientation according to multichannel acoustic pressure signal phase feature.But in implementation process, collection
Obtained sound pressure signal is frequently subjected to the interference such as transformer body harmonics noise, fan noise, and these noises substantially exceed
Noise caused by electric discharge.Therefore, harmonics noise, quick separating discharge noise and the fidelity for keeping its phase are removed, it is right
It is accurately positioned and makes great sense in the shelf depreciation based on phase property.Traditional signal processing of partial discharge method has FFT
Threshold filter method, adaptive digital filtering method, Empirical mode decomposition, wavelet analysis method etc..These methods are deposited to some extent
Threshold value select difficulty, modal overlap, slow and computationally intensive convergence rate the shortcomings of.In recent years, it is special according to signal priori
Sign construction atom function, local discharge signal can be separated using the method for Its Sparse Decomposition, but its most distinct issues is to calculate
Measure huge.
The content of the invention
In view of the above the shortcomings that prior art, it is an object of the invention to provide a kind of local discharge of electrical equipment to determine
Position method and system, storage medium and terminal device, by the spectrum analysis to transformer station's noise, using the more stopbands of linear phase shift
Comb filter filters out harmonics noise, then the period that doubtful shelf depreciation occurs is determined by calculating the kurtosis of residue signal,
Shelf depreciation noise signal is extracted finally by Its Sparse Decomposition, so as to realize the accurate fast positioning of transformer station's shelf depreciation.
In order to achieve the above objects and other related objects, it is local to provide a kind of electrical equipment based on phase property by the present invention
Breakdown location method, comprises the following steps：Obtain the harmonics cut-off frequency of the power frequency of the sound pressure signal of electrical equipment；Using linear
The more stopband comb filter of phase shift obtain to being that the sound pressure signal of electrical equipment is filtered after filtering out the harmonics of power frequency with acquisition
Residual signals；The harmonics centre frequency of the stopband of more stopband comb filter is less than the harmonics cutoff frequency including frequency
Rate and amplitude exceed all harmonics frequencies of the cutoff threshold；The kurtosis of the residual signals is calculated, and sets kurtosis threshold value；
If kurtosis is more than the kurtosis threshold value, within the period corresponding to the kurtosis shelf depreciation occurs for judgement；Based on being judged
The shelf depreciation period and the generation of local discharge signal model include the local discharge signal of phase property；According to shelf depreciation
The phase property of signal come realize shelf depreciation position.
In one embodiment of the invention, the harmonics cut-off frequency of power frequency of the sound pressure signal of electrical equipment is obtained including following
Step：
FFT is carried out to the sound pressure signal of electrical equipment, obtains noise spectrum feature；
The spectrum signature of the harmonics of power frequency in the noise spectrum, highest harmonics amplitude is determined, and hundred will be preset
Divide the highest harmonics amplitude of ratio as cutoff threshold；
If the amplitude more than or equal to the harmonics of a certain frequency is respectively less than cutoff threshold, the frequency is set as harmonics cutoff frequency
Rate.
In one embodiment of the invention, more stopband comb filter are FIR filter.
In one embodiment of the invention, the kurtosis of the residual signals is calculated, and sets kurtosis threshold value；If kurtosis is more than institute
Kurtosis threshold value is stated, then judges that shelf depreciation occurs within the period corresponding to the kurtosis to be comprised the following steps：
It is determined that sliding length of window；
Calculate the kurtosis of the residual signalsWherein i represents the sequence number of kurtosis value, and W is
Slide length of window, SrResidual signals are represented, j represents the sequence number of residual signals in window；
Kurtosis threshold value is set, if kurtosis is more than the kurtosis threshold value, judgement is sent out within the period corresponding to the kurtosis
Raw shelf depreciation.
In one embodiment of the invention, by the upward rounding of sampling number of the sound pressure signal in power frequency period to even number
Value as sliding length of window.
In one embodiment of the invention, the kurtosis threshold value is Represent the high and steep of the residual signals respectively with σ
The average and variance of degree.
In one embodiment of the invention, based on the shelf depreciation period judged and local discharge signal model generation bag
Local discharge signal containing phase property comprises the following steps：
Construct local discharge signal model
And local discharge noise is normalized；Wherein γjIt is the waveform parameter collection of discretization,α tables
Show bandwidth factor, at the time of τ represents that shelf depreciation occurs, f0Partial discharge pulse's centre frequency is represented,Represent initial phase；A
And f0Its span is determined by shelf depreciation electric discharge priori, τ determines according to the period for judging generation shelf depreciation
Its span,It is the numerical parameter for needing optimizing to solve；
For the parameter combination of different local discharge signal models, the shelf depreciation noise signal of normalized is calculated
With the inner product of residual signals, the maximum local discharge signal of inner product is obtained, and obtain the phase information of the local discharge signal.
Correspondingly, the present invention also provides a kind of local discharge of electrical equipment alignment system based on phase property, including obtains
Modulus block, filtration module, determination module, generation module and locating module；
Acquisition module is used for the harmonics cut-off frequency for obtaining the power frequency of the sound pressure signal of electrical equipment；
Filtration module is used for using the more stopband comb filter of linear phase shift to being that the sound pressure signal of electrical equipment is filtered
Ripple, to obtain the residual signals obtained after the harmonics for filtering out power frequency；The harmonics center of the stopband of more stopband comb filter
Frequency includes all harmonics frequencies that frequency exceedes the cutoff threshold less than the harmonics cut-off frequency and amplitude；
Determination module is used for the kurtosis for calculating the residual signals, and sets kurtosis threshold value；If kurtosis is more than the kurtosis
Threshold value, then judge shelf depreciation occurs within the period corresponding to the kurtosis；
Generation module is used to include phase based on the shelf depreciation period judged and the generation of local discharge signal model
The local discharge signal of feature；
Locating module is used to realize that shelf depreciation positions according to the phase property of local discharge signal.
In addition, the present invention also provides a kind of computer-readable recording medium, computer program is stored thereon with, the program quilt
The above-mentioned local discharge of electrical equipment localization method based on phase property is realized during computing device.
Finally, the present invention also provides a kind of terminal device, including：Processor and memory；The memory is used to store
Computer program, the processor is used for the computer program for performing the memory storage, so that the terminal device performs
The above-mentioned local discharge of electrical equipment alignment system method based on phase property.
As described above, the local discharge of electrical equipment localization method and system of the present invention, storage medium and terminal device, tool
There is following beneficial effect：
(1) by filtering out harmonics noise, the signal to noise ratio of local discharge signal is improved, improves and judges whether that part occurs
The accuracy of electric discharge；
(2) by calculating residual signals kurtosis, it can quickly and accurately estimate the doubtful period discharged, subtract significantly
The small amount of calculation of Its Sparse Decomposition algorithm；
(3) filtered by more stopband comb filter with linear phase characteristic, maintain in residual signals and locally put
The phase property of electric signal, to realize that the shelf depreciation noise source based on phase positions；
(4) compared to prior art, not only solving speed is fast, and can keep the feature invariant of each channel phases, from
And realize the discharge noise source based on phase property and be accurately positioned.
Brief description of the drawings
Fig. 1 is shown as the local discharge of electrical equipment localization method based on phase property of the present invention in an embodiment
Flow chart；
Fig. 2 is shown as the local discharge of electrical equipment localization method based on phase property of the present invention in an embodiment
Noise time-domain diagram and spectrogram；
Fig. 3 is shown as the local discharge of electrical equipment localization method based on phase property of the present invention in an embodiment
The curve map of more stopband comb filter amplitude-frequencies and phase frequency；
Fig. 4 is shown as the local discharge of electrical equipment localization method based on phase property of the present invention in an embodiment
Sound pressure signal, residual noise, the curve map of kurtosis and local discharge signal；
Fig. 5 is shown as the local discharge of electrical equipment alignment system based on phase property of the present invention in an embodiment
Structural representation；
Fig. 6 is shown as structural representation of the terminal device of the present invention in an embodiment.
Component label instructions
51 acquisition modules
52 filtration modules
53 determination modules
54 generation modules
55 locating modules
6 terminal devices
61 processors
62 memories
Embodiment
Illustrate embodiments of the present invention below by way of specific instantiation, those skilled in the art can be by this specification
Disclosed content understands other advantages and effect of the present invention easily.The present invention can also pass through specific realities different in addition
The mode of applying is embodied or practiced, the various details in this specification can also be based on different viewpoints with application, without departing from
Various modifications or alterations are carried out under the spirit of the present invention.It should be noted that in the case where not conflicting, following examples and implementation
Feature in example can be mutually combined.
It should be noted that the diagram provided in following examples only illustrates the basic structure of the present invention in a schematic way
Think, only show the component relevant with the present invention in schema then rather than according to component count, shape and the size during actual implement
Draw, kenel, quantity and the ratio of each component can be a kind of random change during its actual implementation, and its assembly layout kenel
It is likely more complexity.
The local discharge of electrical equipment localization method and system based on phase property, the storage medium and terminal of the present invention is set
It is standby to identify, separate shelf depreciation noise, and based on carrying out multichannel shelf depreciation noise linearity phase shift, substantially increase and ask
Speed is solved, and passes through the mutual phase difference of the multi channel signals separated and remains unchanged, so as to realize based on phase property
Discharge noise source is accurately positioned.
As shown in figure 1, in an embodiment, the local discharge of electrical equipment localization method of the invention based on phase property
Comprise the following steps：
Step S1, the harmonics cut-off frequency of the power frequency of the sound pressure signal of electrical equipment is obtained.
Specifically, spectrum analysis is carried out by the sound pressure signal of the electrical equipment to being gathered to obtain the frequency spectrum of noise spy
Sign, then determine the harmonics cut-off frequency of power frequency.Wherein power frequency refers to generating, transmission of electricity, power transformation and the controller switching equipment of power system
And the rated frequency that industry uses with civil electrical equipment, unit hertz Hz.For example, power frequency is 50Hz used by China,
Power frequency is 60Hz used by some countries.
Preferably, the electrical equipment includes transformer.
In an embodiment, the harmonics cut-off frequency for obtaining the power frequency of the sound pressure signal of electrical equipment comprises the following steps：
11) FFT is carried out to the sound pressure signal of electrical equipment, obtains noise spectrum feature.
For transformer, noise signal includes transformer body harmonics noise, discharge noise, fan noise and other
Noise.Wherein, transformer body harmonics noise concentrates on power frequency and the harmonics composition using power frequency as fundamental frequency；Discharge noise is in short-term
Impulsive noise；Fan noise is low frequency periodic noise；Other noises include ambient noise and random noise.
12) spectrum signature of the harmonics of the power frequency in the noise spectrum, highest harmonics amplitude is determined, and will be default
The highest harmonics amplitude of percentage is as cutoff threshold.
Preferably, preset percentage 10%.
If 13) amplitude more than or equal to the harmonics of a certain frequency is respectively less than cutoff threshold, sets the frequency and end as harmonics
Frequency.
Step S2, the sound pressure signal of electrical equipment is filtered using linear phase shift more stopband comb filter, to obtain
Take the residual signals obtained after the harmonics for filtering out power frequency；The harmonics centre frequency of the stopband of more stopband comb filter includes
Frequency exceedes all harmonics frequencies of the cutoff threshold less than the harmonics cut-off frequency and amplitude.
In one embodiment of the invention, the sound pressure signal of electrical equipment is entered using linear phase shift more stopband comb filter
Row filtering, comprised the following steps with obtaining the residual signals obtained after the harmonics for filtering out power frequency：
21) more stopband comb filter are generated, the stopband center frequency of more stopband comb filter is low including frequency
Exceed all harmonics frequencies of the cutoff threshold in the harmonics cut-off frequency and amplitude.
Specifically, according to the spectrogram of noise and harmonics cut-off frequency, determine frequency less than the harmonics cut-off frequency and
Amplitude exceedes stopband center frequency of all harmonics frequencies of the cutoff threshold as more stopband comb filter.It is preferred that
Ground, the stopband center frequency of more stopband comb filter is all harmonics frequencies from power frequency to harmonics cut-off frequency.
Preferably, more stopband comb filter are FIR filter.(Finite Impulse Response, have FIR
Limit for length's unit impulse response) wave filter, also known as non-recursive type wave filter are members most basic in digital information processing system
Part, it can be while any amplitude versus frequency characte be ensured with strict linear phase-frequency characteristic, while its unit sample respo is
It is time-limited, thus wave filter is stable system.Thus it is guaranteed that the linear phase shift of more stopband comb filter is special
Property.
22) sound pressure signal of electrical equipment is filtered using more stopband comb filter, work is filtered out to obtain
The residual signals obtained after the harmonics of frequency.
Step S3, the kurtosis of the residual signals is calculated, and sets kurtosis threshold value；If kurtosis is more than the kurtosis threshold value,
Then judge shelf depreciation occurs within the period corresponding to the kurtosis.
In one embodiment of the invention, the kurtosis of the residual signals is calculated, and sets kurtosis threshold value；If kurtosis is more than institute
Kurtosis threshold value is stated, then judges that shelf depreciation occurs within the period corresponding to the kurtosis to be comprised the following steps：
31) sliding length of window is determined.
Specifically, the upward rounding of sampling number using sound pressure signal in power frequency period to even number value as sliding window
Length.
32) kurtosis of the residual signals is calculatedWherein i represents the sequence number of kurtosis value, W
To slide length of window, SrResidual signals are represented, j represents the sequence number of residual signals in window.
33) kurtosis threshold value is set, if kurtosis is more than the kurtosis threshold value, is judged within the period corresponding to the kurtosis
Generation shelf depreciation.
Specifically, the kurtosis threshold value is Represent average and the side of the kurtosis of the residual signals respectively with σ
Difference.
Step S4, phase property is included based on the shelf depreciation period judged and the generation of local discharge signal model
Local discharge signal.
In one embodiment of the invention, based on the shelf depreciation period judged and local discharge signal model generation bag
Local discharge signal containing phase property comprises the following steps：
41) local discharge signal model is constructed
And local discharge noise is normalized；Wherein γjIt is the waveform parameter collection of discretization,α
Bandwidth factor is represented, at the time of τ represents that shelf depreciation occurs, f0Partial discharge pulse's centre frequency is represented,Represent initial phase.
Specifically, τ ∈ { τi, α ∈ { αk,A and f0Can be discharged priori by shelf depreciation
Knowledge determines its span, and τ determines its span according to the period for judging generation shelf depreciation,It is to need optimizing
The numerical parameter of solution.Preferably, α is a positive constant.
42) parameter combination of different local discharge signal models is directed to, the shelf depreciation for calculating normalized is made an uproar letter
Inner product number with residual signals, the maximum local discharge signal of inner product is obtained, and obtain the phase information of the local discharge signal.
Step S5, realize that shelf depreciation positions according to the phase property of local discharge signal.
Specifically, for multichannel sound pressure signal, each passage is estimated according to the phase property of local discharge signal
Delay, the space orientation of shelf depreciation noise source is realized in conjunction with space coordinates.
The local discharge of electrical equipment based on phase property of the present invention is expanded on further below by specific embodiment
Localization method.Electrical equipment is transformer station in this embodiment, then transformer station's partial discharge positioning method comprises the following steps：
1) the harmonics cut-off frequency of the power frequency of the sound pressure signal of transformer station is obtained.
Specifically, FFT is carried out to the multichannel sound pressure signal of transformer station, obtains noise spectrum feature；Such as Fig. 2 institutes
Show, cutoff threshold is used as using frequency as the 10% of the amplitude of 400Hz highest harmonics；If more than or equal to the harmonics of a certain frequency
Amplitude is respectively less than cutoff threshold, sets the frequency as harmonics cut-off frequency.
Wherein, transformer station's noise mainly includes transformer body noise, shelf depreciation noise, fan noise and other made an uproar
Sound.This bulk noise of transformer or reactor concentrates on 50Hz and the harmonics composition using 50Hz as fundamental frequency；Shelf depreciation noise is
Short term impulsive noise, main frequency range are audible sound section 20Hz~20KHz；Fan noise is low frequency periodic signals, other noises
Including white noise and ambient noise etc..
2) using the more stopband comb filter of linear phase shift to being that the sound pressure signal of electrical equipment is filtered, to obtain filter
Except the residual signals obtained after the harmonics of power frequency；The harmonics centre frequency of the stopband of more stopband comb filter includes frequency
Less than all harmonics frequencies that the harmonics cut-off frequency and amplitude exceed the cutoff threshold.
As shown in figure 3, identified harmonics cut-off frequency and spectrum curve are foundation, more stopband comb filter are determined
Stopband center frequency is [100,250,300,350,400,450,550].
3) kurtosis of the residual signals is calculated, and sets kurtosis threshold value；If kurtosis is more than the kurtosis threshold value, judge
Shelf depreciation occurs within the period corresponding to the kurtosis.
As shown in Figure 4, it is first determined sliding window width when kurtosis calculates；The kurtosis of residual signals is calculated again, and is painted
Kurtosis curve processed；The average and variance of calculating kurtosis finally find the abnormal data of kurtosis to obtain kurtosis threshold value at 0.2s
Section, you can shelf depreciation occurs herein for judgement.
4) part of phase property is included based on the shelf depreciation period judged and the generation of local discharge signal model
Discharge signal.
For the local discharge signal that inner product is maximum, the time delay of the local discharge signal is obtained, phase, frequency, decay refer to
The information such as number.
5) realize that shelf depreciation positions according to the phase property of local discharge signal.
Specifically, the locus of local discharge signal is determined by multichannel phase information
As shown in figure 5, in an embodiment, the local discharge of electrical equipment alignment system of the invention based on phase property
Including acquisition module 51, filtration module 52, determination module 53, generation module 54 and locating module 55.
Acquisition module 51 is used for the harmonics cut-off frequency for obtaining the power frequency of the sound pressure signal of electrical equipment.
Specifically, spectrum analysis is carried out by the sound pressure signal of the electrical equipment to being gathered to obtain the frequency spectrum of noise spy
Sign, then determine the harmonics cut-off frequency of power frequency.Wherein power frequency refers to generating, transmission of electricity, power transformation and the controller switching equipment of power system
And the rated frequency that industry uses with civil electrical equipment, unit hertz Hz.For example, power frequency is 50Hz used by China,
Power frequency is 60Hz used by some countries.
Preferably, the electrical equipment includes transformer.
In an embodiment, the harmonics cut-off frequency that acquisition module 51 obtains the power frequency of the sound pressure signal of electrical equipment performs
Following steps：
A) FFT is carried out to the sound pressure signal of electrical equipment, obtains noise spectrum feature.
For transformer, noise signal includes transformer body harmonics noise, discharge noise, fan noise and other
Noise.Wherein, transformer body harmonics noise concentrates on power frequency and the harmonics composition using power frequency as fundamental frequency；Discharge noise is in short-term
Impulsive noise；Fan noise is low frequency periodic noise；Other noises include ambient noise and random noise.
B) spectrum signature of the harmonics of the power frequency in the noise spectrum, highest harmonics amplitude is determined, and will be default
The highest harmonics amplitude of percentage is as cutoff threshold.
Preferably, preset percentage 10%.
If c) amplitude more than or equal to the harmonics of a certain frequency is respectively less than cutoff threshold, the frequency is set as harmonics cutoff frequency
Rate.
Filtration module 52 is connected with acquisition module 51, for using the more stopband comb filter of linear phase shift to electrical equipment
Sound pressure signal be filtered, filter out the residual signals obtained after the harmonics of power frequency to obtain；More stopband comb filter
The harmonics centre frequency of stopband include frequency and exceed all of the cutoff threshold less than the harmonics cut-off frequency and amplitude
Harmonics frequency.
In one embodiment of the invention, filtration module 52 is using the more stopband comb filter of linear phase shift to electrical equipment
Sound pressure signal is filtered, and following steps are performed to obtain the residual signals bag obtained after the harmonics for filtering out power frequency：
A) more stopband comb filter are generated, the stopband center frequency of more stopband comb filter is less than including frequency
The harmonics cut-off frequency and amplitude exceed all harmonics frequencies of the cutoff threshold.
Specifically, according to the spectrogram of noise and harmonics cut-off frequency, determine frequency less than the harmonics cut-off frequency and
Amplitude exceedes stopband center frequency of all harmonics frequencies of the cutoff threshold as more stopband comb filter.It is preferred that
Ground, the stopband center frequency of more stopband comb filter is all harmonics frequencies from power frequency to harmonics cut-off frequency.
Preferably, more stopband comb filter are FIR filter.(Finite Impulse Response, have FIR
Limit for length's unit impulse response) wave filter, also known as non-recursive type wave filter are members most basic in digital information processing system
Part, it can be while any amplitude versus frequency characte be ensured with strict linear phase-frequency characteristic, while its unit sample respo is
It is time-limited, thus wave filter is stable system.Thus it is guaranteed that the linear phase shift of more stopband comb filter is special
Property.
B) sound pressure signal of electrical equipment is filtered using more stopband comb filter, power frequency is filtered out to obtain
Harmonics after obtained residual signals.
Determination module 53 is connected with filtration module 52, for calculating the kurtosis of the residual signals, and sets kurtosis threshold value；
If kurtosis is more than the kurtosis threshold value, within the period corresponding to the kurtosis shelf depreciation occurs for judgement.
In one embodiment of the invention, determination module 53 calculates the kurtosis of the residual signals, and sets kurtosis threshold value；If
Kurtosis is more than the kurtosis threshold value, then judges that shelf depreciation occurs within the period corresponding to the kurtosis performs following steps：
A) sliding length of window is determined.
Specifically, the upward rounding of sampling number using sound pressure signal in power frequency period to even number value as sliding window
Length.
B) kurtosis of the residual signals is calculatedWherein i represents the sequence number of kurtosis value, W
To slide length of window, SrResidual signals are represented, j represents the sequence number of residual signals in window.
C) kurtosis threshold value is set, if kurtosis is more than the kurtosis threshold value, is judged within the period corresponding to the kurtosis
Generation shelf depreciation.
Specifically, the kurtosis threshold value is Represent average and the side of the kurtosis of the residual signals respectively with σ
Difference.
Generation module 54 is connected with determination module 53 and filtration module 52, for based on the shelf depreciation period judged
The local discharge signal of phase property is included with the generation of local discharge signal model.
In one embodiment of the invention, generation module 54 is based on the shelf depreciation period judged and local discharge signal
Local discharge signal of the model generation comprising phase property performs following steps：
A) local discharge signal model is constructed
And local discharge noise is normalized；Wherein γjIt is the waveform parameter collection of discretization,α tables
Show bandwidth factor, at the time of τ represents that shelf depreciation occurs, f0Partial discharge pulse's centre frequency is represented,Represent initial phase.
Specifically, τ ∈ { τi, α ∈ { αk,A and f0Can be discharged priori by shelf depreciation
Knowledge determines its span, and τ determines its span according to the period for judging generation shelf depreciation,It is to need optimizing
The numerical parameter of solution.Preferably, α is a positive constant.
B) parameter combination of different local discharge signal models is directed to, calculates the shelf depreciation noise letter of normalized
Inner product number with residual signals, the maximum local discharge signal of inner product is obtained, and obtain the phase information of the local discharge signal.
Locating module 55 is connected with generation module 54, and local put is realized for the phase property according to local discharge signal
Electricity positioning.
Specifically, for multichannel sound pressure signal, each passage is estimated according to the phase property of local discharge signal
Delay, the space orientation of shelf depreciation noise source is realized in conjunction with space coordinates.
It should be noted that it should be understood that the modules of system above division be only a kind of division of logic function,
Can completely or partially it be integrated on a physical entity when actually realizing, can also be physically separate.And these modules can be with
All realized in the form of software is called by treatment element；All it can also realize in the form of hardware；Can also part mould
Block calls the form of software to realize by treatment element, and part of module is realized by the form of hardware.For example, x modules can be
The treatment element individually set up, it can also be integrated in some chip of said apparatus and realize, in addition it is also possible to program generation
The form of code is stored in the memory of said apparatus, is called by some treatment element of said apparatus and is performed above x moulds
The function of block.The realization of other modules is similar therewith.In addition these modules can completely or partially integrate, can also be only
It is vertical to realize.Treatment element described here can be a kind of integrated circuit, have the disposal ability of signal.In implementation process,
Each step of the above method or more modules can pass through the integrated logic circuit or soft of the hardware in processor elements
The instruction of part form is completed.
For example, the above module can be arranged to implement one or more integrated circuits of above method, such as：
One or more specific integrated circuits (ApplicationSpecificIntegratedCircuit, abbreviation ASIC), or, one
Or multi-microprocessor (digitalsingnalprocessor, abbreviation DSP), or, one or more field-programmable gate array
Arrange (FieldProgrammableGateArray, abbreviation FPGA) etc..For another example, some module is dispatched by treatment element more than
When the form of program code is realized, the treatment element can be general processor, such as central processing unit
(CentralProcessingUnit, abbreviation CPU) or it is other can be with the processor of caller code.For another example, these modules can
To integrate, realized in the form of on-chip system (system-on-a-chip, abbreviation SOC).
Computer program is stored with the storage medium of the present invention, the program is realized when being executed by processor above-mentioned is based on phase
The local discharge of electrical equipment localization method of position feature.One of ordinary skill in the art will appreciate that：Realize that above-mentioned each method is real
Applying all or part of step of example can be completed by the related hardware of computer program.Foregoing computer program can be deposited
It is stored in a computer-readable recording medium.The program upon execution, execution the step of including above-mentioned each method embodiment；It is and preceding
The storage medium stated includes：ROM, RAM, magnetic disc or CD etc. are various can be with the medium of store program codes.
As shown in fig. 6, in an embodiment, terminal device 6 of the invention includes processor 61 and memory 62.
The memory 62 is used to store computer program, and the processor 61 is used to perform what the memory 62 stored
Computer program, so that the terminal device performs the above-mentioned local discharge of electrical equipment localization method based on phase property.
Preferably, the processor 61 can be general processor, including central processing unit
(CentralProcessingUnit, abbreviation CPU), network processing unit (NetworkProcessor, abbreviation NP) etc.；Can be with
It is digital signal processor (DigitalSignalProcessing, abbreviation DSP), application specific integrated circuit
(ApplicationSpecificIntegratedCircuit, abbreviation ASIC), field programmable gate array (Field-
ProgrammableGateArray, abbreviation FPGA) either other PLDs, discrete gate or transistor logic device
Part, discrete hardware components.
In summary, local discharge of electrical equipment localization method and system of the invention, storage medium and terminal device lead to
Cross and filter out harmonics noise, improve the signal to noise ratio of local discharge signal, improve the accuracy for judging whether that shelf depreciation occurs；
By calculating residual signals kurtosis, it can quickly and accurately estimate the doubtful period discharged, substantially reduce Its Sparse Decomposition
The amount of calculation of algorithm；Filtered, maintained local in residual signals by more stopband comb filter with linear phase characteristic
The phase property of discharge signal, to realize that the shelf depreciation noise source based on phase positions；Compared to prior art, not only solve
Speed is fast, and can keep the feature invariant of each channel phases, so as to realize that the discharge noise source based on phase property is accurate
Positioning.So the present invention effectively overcomes various shortcoming of the prior art and has high industrial utilization.
The above-described embodiments merely illustrate the principles and effects of the present invention, not for the limitation present invention.It is any ripe
Know the personage of this technology all can carry out modifications and changes under the spirit and scope without prejudice to the present invention to above-described embodiment.Cause
This, those of ordinary skill in the art is complete without departing from disclosed spirit and institute under technological thought such as
Into all equivalent modifications or change, should by the present invention claim be covered.
Claims (10)
- A kind of 1. local discharge of electrical equipment localization method based on phase property, it is characterised in that：Comprise the following stepsObtain the harmonics cut-off frequency of the power frequency of the sound pressure signal of electrical equipment；Using the more stopband comb filter of linear phase shift to being that the sound pressure signal of electrical equipment is filtered, power frequency is filtered out to obtain Harmonics after obtained residual signals；The harmonics centre frequency of the stopband of more stopband comb filter is less than institute including frequency State harmonics cut-off frequency and amplitude exceedes all harmonics frequencies of the cutoff threshold；The kurtosis of the residual signals is calculated, and sets kurtosis threshold value；If kurtosis is more than the kurtosis threshold value, judge described Shelf depreciation occurs in period corresponding to kurtosis；Based on the shelf depreciation period judged and the generation of local discharge signal model, the shelf depreciation comprising phase property is believed Number；Realize that shelf depreciation positions according to the phase property of local discharge signal.
- 2. the local discharge of electrical equipment localization method according to claim 1 based on phase property, it is characterised in that：Obtain The harmonics cut-off frequency of the power frequency of the sound pressure signal of electrical equipment is taken to comprise the following steps：FFT is carried out to the sound pressure signal of electrical equipment, obtains noise spectrum feature；The spectrum signature of the harmonics of power frequency in the noise spectrum, determines highest harmonics amplitude, and by preset percentage Highest harmonics amplitude as cutoff threshold；If the amplitude more than or equal to the harmonics of a certain frequency is respectively less than cutoff threshold, the frequency is set as harmonics cut-off frequency.
- 3. the local discharge of electrical equipment localization method according to claim 1 based on phase property, it is characterised in that：Institute It is FIR filter to state more stopband comb filter.
- 4. the local discharge of electrical equipment localization method according to claim 1 based on phase property, it is characterised in that：Meter The kurtosis of the residual signals is calculated, and sets kurtosis threshold value；If kurtosis is more than the kurtosis threshold value, judge in the kurtosis pair Shelf depreciation occurs in the period answered to comprise the following steps：It is determined that sliding length of window；Calculate the kurtosis of the residual signalsWherein i represents the sequence number of kurtosis value, and W is sliding Length of window, SrResidual signals are represented, j represents the sequence number of residual signals in window；Kurtosis threshold value is set, if kurtosis is more than the kurtosis threshold value, judges the generation office within the period corresponding to the kurtosis Discharge in portion.
- 5. the local discharge of electrical equipment localization method according to claim 4 based on phase property, it is characterised in that：Will The value of the upward rounding of sampling number of the sound pressure signal in power frequency period to even number is used as sliding length of window.
- 6. the local discharge of electrical equipment localization method according to claim 4 based on phase property, it is characterised in that：Institute Stating kurtosis threshold value is Represent the average and variance of the kurtosis of the residual signals respectively with σ.
- 7. the local discharge of electrical equipment localization method according to claim 1 based on phase property, it is characterised in that：Base In the shelf depreciation period judged and the generation of local discharge signal model, the local discharge signal comprising phase property includes Following steps：Construct local discharge signal modelAnd Local discharge noise is normalized；Wherein γjIt is the waveform parameter collection of discretization,α is represented Bandwidth factor, at the time of τ represents that shelf depreciation occurs, f0Partial discharge pulse's centre frequency is represented,Represent initial phase；A and f0Its span is determined by shelf depreciation electric discharge priori, τ determines it according to the period for judging generation shelf depreciation Span,It is the numerical parameter for needing optimizing to solve；For the parameter combination of different local discharge signal models, calculate the shelf depreciation noise signal of normalized with it is residual The inner product of difference signal, the maximum local discharge signal of inner product is obtained, and obtain the phase information of the local discharge signal.
- A kind of 8. local discharge of electrical equipment alignment system based on phase property, it is characterised in that：Including acquisition module, filtering Module, determination module, generation module and locating module；Acquisition module is used for the harmonics cut-off frequency for obtaining the power frequency of the sound pressure signal of electrical equipment；Filtration module is used for using the more stopband comb filter of linear phase shift to being that the sound pressure signal of electrical equipment is filtered, with Obtain the residual signals obtained after the harmonics for filtering out power frequency；The harmonics centre frequency bag of the stopband of more stopband comb filter Include all harmonics frequencies that frequency exceedes the cutoff threshold less than the harmonics cut-off frequency and amplitude；Determination module is used for the kurtosis for calculating the residual signals, and sets kurtosis threshold value；If kurtosis is more than the kurtosis threshold value, Then judge shelf depreciation occurs within the period corresponding to the kurtosis；Generation module is used to include phase property based on the shelf depreciation period judged and the generation of local discharge signal model Local discharge signal；Locating module is used to realize that shelf depreciation positions according to the phase property of local discharge signal.
- 9. a kind of computer-readable recording medium, is stored thereon with computer program, it is characterised in that the program is held by processor The local discharge of electrical equipment localization method based on phase property any one of claim 1 to 7 is realized during row.
- A kind of 10. terminal device, it is characterised in that including：Processor and memory；The memory is used to store computer program, and the processor is used for the computer journey for performing the memory storage Sequence, so that the terminal device performs the electrical equipment based on phase property as any one of claim 1 to 7 and locally put Electric localization method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710942874.2A CN107831412B (en) | 2017-10-11 | 2017-10-11 | Electrical equipment partial discharge positioning method and system, storage medium and terminal equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710942874.2A CN107831412B (en) | 2017-10-11 | 2017-10-11 | Electrical equipment partial discharge positioning method and system, storage medium and terminal equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107831412A true CN107831412A (en) | 2018-03-23 |
| CN107831412B CN107831412B (en) | 2020-02-21 |
Family
ID=61648002
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710942874.2A Active CN107831412B (en) | 2017-10-11 | 2017-10-11 | Electrical equipment partial discharge positioning method and system, storage medium and terminal equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107831412B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112150474A (en) * | 2020-10-12 | 2020-12-29 | 山东省科学院海洋仪器仪表研究所 | Underwater bubble image feature segmentation and extraction method |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20080076328A (en) * | 2007-02-15 | 2008-08-20 | 주식회사 효성 | Electric power appliance remote part discharge diagnosis system and method to use frequency/status spectra |
| CN104316846A (en) * | 2014-08-28 | 2015-01-28 | 国家电网公司 | Identification method, apparatus, and system of partial discharger mode of electrical equipment |
| CN105606966A (en) * | 2015-12-21 | 2016-05-25 | 安徽理工大学 | Partial discharge pattern recognition method based on mixed neural network algorithm |
| CN106501690A (en) * | 2016-10-28 | 2017-03-15 | 国网辽宁省电力有限公司辽阳供电公司 | A kind of XLPE power cables partial discharge diagnostic method and system |
| CN206161731U (en) * | 2016-10-31 | 2017-05-10 | 国家电网公司 | On -line monitoring is put with digital power frequency phase sensor in office of gas -insulated transformer substation |
| CN107036709A (en) * | 2017-04-11 | 2017-08-11 | 国网上海市电力公司经济技术研究院 | A kind of transformer station's noise matching separation method |
| CN107167716A (en) * | 2017-07-11 | 2017-09-15 | 国网福建省电力有限公司泉州供电公司 | A kind of shelf depreciation default kind identification method and device |
-
2017
- 2017-10-11 CN CN201710942874.2A patent/CN107831412B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20080076328A (en) * | 2007-02-15 | 2008-08-20 | 주식회사 효성 | Electric power appliance remote part discharge diagnosis system and method to use frequency/status spectra |
| CN104316846A (en) * | 2014-08-28 | 2015-01-28 | 国家电网公司 | Identification method, apparatus, and system of partial discharger mode of electrical equipment |
| CN105606966A (en) * | 2015-12-21 | 2016-05-25 | 安徽理工大学 | Partial discharge pattern recognition method based on mixed neural network algorithm |
| CN106501690A (en) * | 2016-10-28 | 2017-03-15 | 国网辽宁省电力有限公司辽阳供电公司 | A kind of XLPE power cables partial discharge diagnostic method and system |
| CN206161731U (en) * | 2016-10-31 | 2017-05-10 | 国家电网公司 | On -line monitoring is put with digital power frequency phase sensor in office of gas -insulated transformer substation |
| CN107036709A (en) * | 2017-04-11 | 2017-08-11 | 国网上海市电力公司经济技术研究院 | A kind of transformer station's noise matching separation method |
| CN107167716A (en) * | 2017-07-11 | 2017-09-15 | 国网福建省电力有限公司泉州供电公司 | A kind of shelf depreciation default kind identification method and device |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112150474A (en) * | 2020-10-12 | 2020-12-29 | 山东省科学院海洋仪器仪表研究所 | Underwater bubble image feature segmentation and extraction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107831412B (en) | 2020-02-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105640545B (en) | Fetal electrocardiosignal extraction method and device | |
| EP3163574B1 (en) | Method and apparatus for detecting abnormal frame | |
| CN110967599A (en) | Electric energy quality disturbance detection and positioning algorithm | |
| CN105508276B (en) | A kind of City Substation main transformer fan noise evaluation method | |
| CN102792173B (en) | Instrument and method for detecting partial electrical discharges | |
| US20190036333A1 (en) | Method and Device for Online Limit Early-Warning to Shunt Capacitor Bank | |
| CN105928701B (en) | Effective IMF decision method during EMD based on correlation analysis | |
| CN105486938B (en) | A kind of substation's mixed noise separation method | |
| CN111160146B (en) | Hydroelectric generating set state monitoring signal digital filtering method, device and system based on time-frequency conversion | |
| CN105046025B (en) | A kind of method being respectively separated in nuclear magnetic resonance multiphase flow measurement | |
| CN111239565B (en) | Oil-filled casing partial discharge pulse signal processing method and system based on layered denoising model | |
| CN107831412A (en) | Local discharge of electrical equipment localization method and system, storage medium and terminal device | |
| CN103197143A (en) | Harmonic and inter-harmonic detection method based on Hanning-window FFT algorithm and traversal filtering | |
| Feizifar et al. | Application of continuous wavelet transform for fault location in combined overhead line and cable distribution networks | |
| Silaghi et al. | Voice command of an industrial robot in a noisy environment | |
| CN111948454B (en) | Synchronous phasor measurement method and system for resisting attenuated direct current component interference | |
| CN103007431B (en) | Device and method for shockable signal detection | |
| CN111289800A (en) | Small-resistance vibration monitoring method based on generalized regression neural network | |
| JP6343771B2 (en) | Head related transfer function modeling apparatus, method and program thereof | |
| CN115856429A (en) | Current harmonic detection method, system and storage medium | |
| CN114034375B (en) | Ultra-high voltage transmission line noise measurement system and method | |
| CN106352974B (en) | A kind of digital sound level meter pulse weighted method and device | |
| CN107765074A (en) | The detection method and detecting system of traction convertor harmonic current signal intensity | |
| CN103929143A (en) | Method for acquiring amplitude-frequency response characteristic of filter | |
| CN117894341B (en) | Method, device, equipment and medium for identifying abnormal rhythm noise of engine |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |