CN102662132A - Method and system for locating partial discharge of electrical equipment - Google Patents

Method and system for locating partial discharge of electrical equipment Download PDF

Info

Publication number
CN102662132A
CN102662132A CN2012101732313A CN201210173231A CN102662132A CN 102662132 A CN102662132 A CN 102662132A CN 2012101732313 A CN2012101732313 A CN 2012101732313A CN 201210173231 A CN201210173231 A CN 201210173231A CN 102662132 A CN102662132 A CN 102662132A
Authority
CN
China
Prior art keywords
signal
array
algorithm
electrical equipment
matrix
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
Application number
CN2012101732313A
Other languages
Chinese (zh)
Other versions
CN102662132B (en
Inventor
谢庆
律方成
程述一
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
State Grid Economic and Technological Research Institute
Original Assignee
North China Electric Power University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN201210173231.3A priority Critical patent/CN102662132B/en
Publication of CN102662132A publication Critical patent/CN102662132A/en
Application granted granted Critical
Publication of CN102662132B publication Critical patent/CN102662132B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

A method for locating partial discharge of electrical equipment is used for detecting and locating partial discharge of the electrical equipment. The technical scheme includes the steps: firstly, using three ultrasonic array sensors to receive partial discharge ultrasonic signals at different positions of the electrical equipment and forming an array model; then, converting broadband signals received by the three ultrasonic array sensors into narrow-band signals by applying a TLS (transport layer security)-based broadband focus algorithm, and obtaining azimuth information of partial discharge sources according to a narrow-band direction-finding algorithm combining the phase matching principle and cloud optimization search; and finally, determining specific positions of electrical equipment partial discharge by means of a partial discharge source location algorithm based on global optimization search of a modified genetic algorithm. The invention simultaneously provides an improved system for locating partial discharge. Compared with a traditional method for locating partial discharge, the method has the advantages of small calculation amount, high direction-finding precision, accuracy in location, high convergence rate and the like, and the system for locating partial discharge is simple in structure, low in cost and high in location precision.

Description

A kind of local discharge of electrical equipment localization method and system
Technical field
The present invention relates to a kind of method and detection system that is used to detect the inside electric appliance partial discharge position, belong to the detection technique field.
Background technology
The detection of shelf depreciation is with accurately the location is significant for the latency insulation defect of timely discovery inside electric appliance, trouble-saving generation.
The detection method of shelf depreciation has multiple, and supercritical ultrasonics technology is used widely because of having advantages such as principle is simple, the location is convenient.But because ultrasound wave is comparatively complicated in the inner propagation of transformer; Common sonac is difficult to detect ultrasonic signal in some cases; Therefore adopt supersonic array sensor to replace single sonac in the correlative study, utilize array signal process technique to realize the accurate location of shelf depreciation with strong jamming inhibition ability and high-space resolution ability.
Research shows that the array signal localization method is a kind of effective means that present local discharge of electrical equipment detects.But not enough below existing supersonic array localization method exists:
(1) deficiency of shelf depreciation supersonic array localization method aspect.
Traditional array signal direction finding algorithm generally is applicable to the direction finding estimation of narrow band signal, and office's discharge signal of electrical equipment is a kind of typical broadband signal, and existing focusing algorithm operand is bigger, and precision is not high, is not suitable for the processing of electrical equipment office discharge signal;
Traditional M USIC direction finding algorithm the convergence speed is slow, and under low signal-to-noise ratio and condition of small sample, can not differentiate the closer signal of space length;
The traditional geometry localization method like simple, the convenience of calculation of independent positioning method principle in direction finding Cross Location Method and the two platform common vertical lines, can be found out the geometric space position of Partial Discharge Sources to a certain extent.But still there are problems such as precision is not enough, error amplification in these methods.
(2) deficiency of shelf depreciation supersonic array signal antinoise method aspect.
Environment of living in is extremely complicated during electric equipment operation, has strong ground unrest on every side, because the type local-discharge ultrasonic signal is very faint, usually is buried in the strong background noise, makes the precision of Partial Discharge Detection descend greatly even causes and detect failure.Traditional denoising method has the good denoising effect for single channel, single array element ultrasonic signal; But; For electrical equipment type local-discharge ultrasonic array localization method, lack the ripe denoising method of a cover at present, if use traditional denoising method; Then the phase information of array signal is lost seriously, and this just office put the critical data of location.
(3) deficiency of shelf depreciation supersonic array Sensor Design aspect.
Used supersonic array sensor mostly adopts planar square array and linear array sensor in the existing localization method, and the planar square sensor array can acquisition office be put two dimension position angle, source and angle of pitch information, but because element number of array is many, cost is higher; The linear array Fundamentals of Sensors are simple, when guaranteeing certain direction finding precision, can reduce redundant array element, greatly reduce cost, but can only the office of obtaining put the one dimension azimuth information in source, ignored the spatial feature of signal.
Summary of the invention
The objective of the invention is to overcome prior art deficiency, a kind of new local discharge of electrical equipment localization method is provided, realize the accurate location of shelf depreciation, simultaneously, the present invention also improves the supersonic array sensor.
Problem according to the invention realizes with following technical proposals:
A kind of local discharge of electrical equipment localization method, it at first utilizes three supersonic array sensors to receive the shelf depreciation ultrasonic signal at the diverse location of electrical equipment, and forms Array Model; The wideband focusing algorithm of using then based on total least square method (TLS) converts the broadband signal that three supersonic array sensors are received to narrow band signal; The azimuth information in source is put in the arrowband direction finding algorithm acquisition office that adopts phase matching principle and cloud optimization searching to combine again, adopts the office of global optimization search to put the particular location that the source location algorithm is confirmed local discharge of electrical equipment at last.
Above-mentioned local discharge of electrical equipment localization method, it may further comprise the steps:
A, utilize three supersonic array sensors to receive the shelf depreciation ultrasonic signal at the diverse location of electrical equipment;
B, the wideband focusing algorithm of using based on TLS convert the broadband signal that three supersonic array sensors are received to narrow band signal:
Suppose that the wideband array signal model is:
Figure 47486DEST_PATH_IMAGE001
In the formula,
Figure 981944DEST_PATH_IMAGE002
It is the reception data vector of array signal;
Figure 240887DEST_PATH_IMAGE003
The direction matrix, whole broadband be divided into into JIndividual arrowband;
Figure 526375DEST_PATH_IMAGE004
Data vector for original signal;
Figure 263386DEST_PATH_IMAGE005
Be that the j section receives data noise spatial data matrix.
The basic step of conversion is:
The initial value of
Figure 317930DEST_PATH_IMAGE006
estimated signal, and selected reference frequency
Figure 747774DEST_PATH_IMAGE007
;
Figure 5712DEST_PATH_IMAGE008
Utilize the initial value of signal, construct the array flow pattern of each Frequency point
Figure 546414DEST_PATH_IMAGE009
, wherein,
Figure 721044DEST_PATH_IMAGE010
Be jSection receives the array of data flow pattern,
Figure 56210DEST_PATH_IMAGE011
Be the focusing transform matrix;
Figure 316290DEST_PATH_IMAGE012
constructs the focussing matrix of each Frequency point:
Figure 395105DEST_PATH_IMAGE013
; Wherein,
Figure 424241DEST_PATH_IMAGE014
is the focussing matrix of each Frequency point;
Figure 930308DEST_PATH_IMAGE015
, the matrix that left singular vector that
Figure 926952DEST_PATH_IMAGE016
and is
Figure 693100DEST_PATH_IMAGE018
respectively and right singular vector are formed;
Figure 635648DEST_PATH_IMAGE019
utilizes a series of focussing matrix array reception data to carry out focusing transform, obtains the data covariance matrix of single-frequency point.
The azimuth information in source is put in c, the arrowband direction finding algorithm acquisition office that adopts phase matching principle and cloud optimization searching to combine:
Suppose to arrive initial point OAcoustical signal do s( t), and the signal amplitude of arrival primitive is identical, n( t) be white noise, the reception signal of primitive is respectively so:
Figure 604741DEST_PATH_IMAGE020
?,
In the formula
Figure 25358DEST_PATH_IMAGE021
,
Figure 29086DEST_PATH_IMAGE022
Be array element distance, Be the azimuth information of signal source,
Figure 84078DEST_PATH_IMAGE024
Be the velocity of sound, NBe the length of signal data section, signal frequency-domain was expressed as after following formula was done the conversion of Fu's formula:
Figure 308386DEST_PATH_IMAGE025
Figure 166620DEST_PATH_IMAGE026
is respectively the phase place of signal, wanted signal and the noise of reception in the formula, and its matrix form is:
Figure 450971DEST_PATH_IMAGE027
?,
Adopt the corresponding angle of maximal value of cloud optimization searching algorithm search DOA estimation criterion
Figure 394656DEST_PATH_IMAGE028
; The azimuth information of can the office of obtaining putting the source; Be the position angle and the angle of pitch of direction finding line
Wherein, Tr (.) is the diagonal element sum of matrix,
Figure 157076DEST_PATH_IMAGE029
Figure 869817DEST_PATH_IMAGE030
The particular location that the source location algorithm is confirmed local discharge of electrical equipment is put in d, the office that adopts global optimization to search for:
At first set up the equation that certain spatial is put each bar direction finding linear distance sum according to the home position of position angle, the angle of pitch and the sensor array of the direction finding line of three supersonic array sensors; Then global search its apart from the sum minimum value, this is worth the locus that the source is put in the pairing spatial point office of being.
Above-mentioned local discharge of electrical equipment localization method before the broadband signal that the supersonic array sensor is received focuses on, should adopt the algorithm based on quick independent component analysis (FastICA) that signal is carried out denoising, and concrete steps are:
1. play a game and put the supersonic array signal and carry out centralization (going average) and whitening pretreatment:
Centralization: each group observation data
Figure 325069DEST_PATH_IMAGE031
thus deduct its average to make it become an average be 0 variable;
Albefaction: the covariance matrix to mixing observation signal
Figure 5318DEST_PATH_IMAGE032
carries out characteristic value decomposition; Obtain feature matrix
Figure 571429DEST_PATH_IMAGE033
and be the diagonal matrix
Figure 138676DEST_PATH_IMAGE034
of diagonal element with the eigenwert; Then the signal after the albefaction is ; Wherein
Figure 683107DEST_PATH_IMAGE036
is the albefaction matrix.
2. the blind separation of type local-discharge ultrasonic array signal
Expectation function based on the FastICA algorithm of plural number is:
Figure 318488DEST_PATH_IMAGE037
Figure 474663DEST_PATH_IMAGE038
is separating vector in the formula, several kinds of functions below function
Figure 553608DEST_PATH_IMAGE039
can be selected:
Figure 162444DEST_PATH_IMAGE040
Figure 867095DEST_PATH_IMAGE041
;
Figure 877776DEST_PATH_IMAGE042
all is approximately 0.1; The maximum value of calculating target function, the two step formulas that obtain plural FastICA are:
Figure 473023DEST_PATH_IMAGE044
In the formula;
Figure 450206DEST_PATH_IMAGE045
respectively representative function
Figure 580973DEST_PATH_IMAGE046
once reach second derivative
When having a plurality of independent source signal; For fear of the same source signal of repeated isolation; Before
Figure 500256DEST_PATH_IMAGE047
individual vector is separated;
Figure 83685DEST_PATH_IMAGE048
for iteration generation each time, carry out following operation:
Figure 849832DEST_PATH_IMAGE050
In the formula, the separation matrix that
Figure 690749DEST_PATH_IMAGE051
expression is made up of preceding
Figure 761473DEST_PATH_IMAGE052
individual separating vector;
3. the extraction of type local-discharge ultrasonic signal:
Calculate waveform similarity coefficient (
Figure 346039DEST_PATH_IMAGE053
):
Figure 467709DEST_PATH_IMAGE054
Figure 948369DEST_PATH_IMAGE055
representes that respectively the Bureau of Standards puts the type local-discharge ultrasonic signal of ultrasonic signal and extraction, extracts the maximum signal of absolute value wherein as the type local-discharge ultrasonic signal;
4. the array flow pattern of type local-discharge ultrasonic signal is estimated:
Being estimated as of array flow pattern:
Figure 771969DEST_PATH_IMAGE056
;
5. the reconstruct of type local-discharge ultrasonic array signal:
Be the type local-discharge ultrasonic signal if confirmed
Figure 629066DEST_PATH_IMAGE047
road output signal ; With the isolated component output matrix after other each road signal zero setting is
Figure 506073DEST_PATH_IMAGE058
, the array signal
Figure 82547DEST_PATH_IMAGE059
after then recovering.
The local discharge of electrical equipment positioning system a kind of as above-mentioned localization method uses; Comprise supersonic array sensor, multi-channel synchronous data acquisition device and computing machine in the formation; Said supersonic array sensor comprises Stainless Steel Shell and a plurality of piezoelectricity array elements that are installed in the Stainless Steel Shell; Said a plurality of piezoelectricity array element is arranged in cruciform; Their signal receiving end is towards the signal receiving port of Stainless Steel Shell, and the signal output part of a plurality of piezoelectricity array elements connects the different input ends of multi-channel synchronous data acquisition device respectively through the array element lead-in wire; The signal output part of said multi-channel synchronous data acquisition device connects computing machine.
Above-mentioned local discharge of electrical equipment positioning system, the signal receiving port of said Stainless Steel Shell is provided with the stainless steel cuticula, between Stainless Steel Shell and piezoelectricity array element, is provided with the silicon rubber filling thing.
Above-mentioned local discharge of electrical equipment positioning system is provided with the acoustic damping piece between piezoelectricity array element and array element lead-in wire.
Above-mentioned local discharge of electrical equipment positioning system, said piezoelectricity array element is provided with nine.
The present invention adopts the wideband focusing algorithm based on TLS that the type local-discharge ultrasonic broadband signal is focused on; The azimuth information in source is put in the arrowband direction finding algorithm acquisition office that adopts phase matching principle and cloud optimization searching to combine, and adopts three platform direction findings and correction genetic algorithm to confirm the position of shelf depreciation.Compare with traditional partial discharge positioning method, the present invention has advantages such as operand is little, direction finding precision is high, accurate positioning, fast convergence rate.
The piezoelectricity array element of the used supersonic array sensor of the present invention is arranged in cruciform; Can record the deflection and the angle of pitch of type local-discharge ultrasonic ripple signal simultaneously; But element number of array obviously reduces than rectangle plane type sensor array; So just when guaranteeing position success rate and bearing accuracy, simplify the structure of detection system greatly, reduced the cost of detection system.
Description of drawings
Below in conjunction with accompanying drawing the present invention is further specified.
Fig. 1 is nine yuan of cruciform supersonic arrays of the instant-plugging sensor vertical view among the present invention;
Fig. 2 is nine yuan of cruciform supersonic arrays of the instant-plugging sensor sectional view among the present invention;
Fig. 3 among the present invention based on the shelf depreciation positioning system synoptic diagram of supersonic array sensor;
Fig. 4 is the process flow diagram of FastICA algorithm among the present invention;
Fig. 5 among the present invention based on the process flow diagram of the searching algorithm of CTOA;
Fig. 6 is three platform direction finding positioning principle synoptic diagram among the present invention;
Fig. 7 is for revising the process flow diagram of genetic algorithm among the present invention;
Fig. 8 is the noisy type local-discharge ultrasonic array signal of 9 a passages time domain waveform;
Fig. 9 is the direction finding spectrogram of original signal;
Figure 10 is the contour map of original signal;
Figure 11 is the type local-discharge ultrasonic signal after extracting;
Figure 12 is that 9 passage restructuring arrays receive time domain plethysmographic signal figure;
Figure 13 is the direction finding spectrogram behind the signal denoising;
Figure 14 is the contour map behind the signal denoising.
Each label is among the figure: 1, cable; 2, gland nut; 3, stainless steel cover; 4, array element lead-in wire; 5, shape of the mouth as one speaks circle; 6, acoustic damping piece; 7, piezoelectricity array element; 8, Stainless Steel Shell; 9, silicon rubber filling thing; 10, stainless steel cuticula; 11, pressurization and coupling device; 12, bushing; 13, discharging model and stationary installation; 14, electrical equipment fuel tank; 15, supersonic array sensor; 16, multi-channel synchronous data acquisition device; 17, computing machine.
The text of each symbol is:
Figure 477757DEST_PATH_IMAGE060
is the frequency j, focusing matrix;
Figure 818695DEST_PATH_IMAGE061
direction matrix;
Figure 906736DEST_PATH_IMAGE062
The data for the signal vector; is the j noise spatial data segment received data matrix;
Figure 638249DEST_PATH_IMAGE064
is the reference frequency;
Figure 838286DEST_PATH_IMAGE065
is
Figure 97229DEST_PATH_IMAGE066
left singular vector matrix composed;
Figure 117138DEST_PATH_IMAGE067
is
Figure 119729DEST_PATH_IMAGE066
right singular vector matrix composed;
Figure 925005DEST_PATH_IMAGE024
is the speed of sound; Tr (.) is the diagonal elements of the matrix, and; FastICA is fast independent component analysis;
Figure 89270DEST_PATH_IMAGE068
is the standard PD ultrasonic signals; is the extraction of PD ultrasound signal.
Embodiment
The present invention aims to provide a kind of supersonic array sensor and novel array signal process technique of utilizing and realizes the method and system of local discharge of electrical equipment location.
The present invention has designed a kind of shelf depreciation positioning system based on the supersonic array sensor, uses it for to receive the shelf depreciation ultrasonic signal, and forms Array Model; Proposed a kind of then based on the type local-discharge ultrasonic array signal Denoising Algorithm of quick independent component analysis and the overall evaluation system of denoising effect; On this basis, proposed a cover complete comprise focusings of type local-discharge ultrasonic array signal wideband, arrowband direction finding, sterically defined method, for accurately assessing the local discharge of electrical equipment position effective reference is provided.
For realizing above-mentioned purpose, concrete technical scheme is following:
(1) a kind of shelf depreciation positioning system based on the supersonic array sensor.
1. nine yuan of cruciform supersonic arrays of instant-plugging Sensor Design scheme
Like accompanying drawing 1, shown in 2, nine yuan of cruciforms of arrangement mode of supersonic array sensor 15 among the present invention, structure, material are identical.Piezoelectricity array element 7 center distance are 6mm, the half-wavelength of corresponding common type local-discharge ultrasonic ripple signal.Each piezoelectricity array element 7 cross section is a disc, and radius 2mm, array element thickness are 12mm, and each right cylinder array element afterbody installs acoustic damping piece 6 additional, and thickness is 4mm, and uniform silicon rubber filling thing 9 is smeared in the periphery, and thickness is 0.5mm.
Supersonic array sensor among the present invention can be according to actual conditions; Array element is inserted in the center distance hole such as corresponding; And be coupled directly to the electrical equipment outer wall; Insert multi-channel synchronous data acquisition device 16 through 9 cables, can obtain the type local-discharge ultrasonic signal of coupling and form Array Model, for realizing establishing good basis based on the type local-discharge ultrasonic detection of supersonic array sensor.
2. based on the shelf depreciation positioning system of supersonic array sensor
Shelf depreciation positioning system based on the supersonic array sensor comprises: supersonic array sensor 15, multi-channel synchronous data acquisition device 16, computing machine 17 and connecting line thereof etc., its structural representation is shown in accompanying drawing 3.
Wherein, multi-channel synchronous data acquisition device 16 is 9 passages, can realize functions such as data acquisition, Signal Pretreatment (comprising filtering, the amplification of supersonic array signal), data high-speed transmission.Wherein the gain that can realize of each passage is 40dB, and filter range is 20~250kHz, and SF is 256kHz~10MHz, and sampled data transfers to computing machine 17 by USB interface.
(2) one covers are based on the partial discharge positioning method of array signal process technique.
1. based on the type local-discharge ultrasonic array signal Denoising Algorithm of quick independent component analysis and the comprehensive evaluation of denoising effect
(1) based on the type local-discharge ultrasonic array signal Denoising Algorithm of quick independent component analysis
1. the pre-service of type local-discharge ultrasonic array signal
Be shortcut calculation, improve travelling speed,, comprise centralization (going average) and albefaction two parts carrying out to carry out pre-service to mixed signal before the FastICA algorithm.
Centralization is the most basic pre-service, just each group observation data
Figure 402757DEST_PATH_IMAGE070
thus deduct its average to make it become an average be 0 variable.
To mixing the albefaction of observation signal ; Be exactly through certain linear transformation ; Make the matrix after the albefaction correlation matrix satisfies ; Through removing the correlativity between each component of signal, make that second-order statistics is independent between the component of the signal after the albefaction.Concrete grammar is that the covariance matrix to observation signal
Figure 529850DEST_PATH_IMAGE075
carries out characteristic value decomposition; Obtain feature matrix
Figure 35918DEST_PATH_IMAGE033
and be the diagonal matrix
Figure 783294DEST_PATH_IMAGE034
of diagonal element with the eigenwert, the albefaction matrix is
Figure 665800DEST_PATH_IMAGE076
.
2. the blind separation of type local-discharge ultrasonic array signal
Can know by the central limit theorem in the theory of probability; Under certain condition; The distribution of a plurality of independent random variable sums more is tending towards Gaussian distribution, the mixed signal of promptly forming
Figure 549442DEST_PATH_IMAGE077
by several independent source signals than any source signal
Figure 226411DEST_PATH_IMAGE078
all more near Gaussian distribution.Therefore; Can be with non-Gauss's property of separating the output signal
Figure 195504DEST_PATH_IMAGE079
after mixing as optimization aim; When non-Gauss's property of a certain output signal
Figure 881700DEST_PATH_IMAGE080
reaches maximum; Explain that it can not be mixed signal
Figure 370582DEST_PATH_IMAGE081
more; Just approach certain source signal
Figure 15190DEST_PATH_IMAGE082
more, promptly maximized non-Gauss's property will produce an estimation to the independent source signal.Characteristic by kurtosis can be known; The kurtosis of Gaussian random variable equals zero; The kurtosis of owing gaussian distributed random variable is a negative value, the kurtosis of this distribution variables of superelevation be on the occasion of, therefore adopting the non-Gauss's property tolerance based on the conduct of kurtosis criterion is the good criterion of ICA.
Expectation function based on the FastICA algorithm of plural number does
Figure 940420DEST_PATH_IMAGE083
(1)
(2)
Figure 22963DEST_PATH_IMAGE038
is separating vector in the formula; Choose and the choosing of real number field of function
Figure 838472DEST_PATH_IMAGE039
are distinguishing, below several kinds of functions be proved to be and be applicable to complex field.
Figure 765846DEST_PATH_IMAGE085
Figure 793844DEST_PATH_IMAGE041
,
Figure 506586DEST_PATH_IMAGE042
all is approximately 0.1.The maximum value of calculating target function, the two step formulas that obtain the quick ICA of plural number are:
Figure 492996DEST_PATH_IMAGE086
In the formula,
Figure 923977DEST_PATH_IMAGE087
respectively representative function
Figure 490088DEST_PATH_IMAGE046
once reach second derivative.When there being a plurality of independent source signals; For fear of the same source signal of repeated isolation, need to produce a plural separation matrix (forming) observation signal that is mixed by a plurality of independent source signals is separated by all separating vectors.Before
Figure 434221DEST_PATH_IMAGE047
individual vector is separated;
Figure 86920DEST_PATH_IMAGE048
for iteration generation each time, carry out following operation:
Figure 456721DEST_PATH_IMAGE088
(8)
Figure 144054DEST_PATH_IMAGE089
(9)
In the formula, the separation matrix that
Figure 941109DEST_PATH_IMAGE051
expression is made up of preceding
Figure 81103DEST_PATH_IMAGE090
individual separating vector.
3. the extraction of type local-discharge ultrasonic signal
Because the amplitude after FastICA separates and the uncertainty of ordering; The present invention adopts waveform similarity coefficient (
Figure 769443DEST_PATH_IMAGE053
) to extract the type local-discharge ultrasonic signal
Figure 780124DEST_PATH_IMAGE091
(10)
Figure 279238DEST_PATH_IMAGE092
representes that respectively the Bureau of Standards puts the type local-discharge ultrasonic signal of ultrasonic signal and extraction.The value of
Figure 375370DEST_PATH_IMAGE053
is between-1 to 1, and waveform was reverse before and after-1 expression separated; 0 represents two waveform quadratures; 1 representative is identical.The absolute value of
Figure 618133DEST_PATH_IMAGE053
is more near 1, the waveform that is near the mark more of the type local-discharge ultrasonic waveform after just explanation separates.So choose the foundation that the type local-discharge ultrasonic signal is extracted in the wherein maximum conduct of absolute value.
4. the array flow pattern of type local-discharge ultrasonic signal is estimated
Because in the array location model; It is decision accurate positioning whether important parameter that source signal arrives the phase information that sensor produced; Therefore need not to consider the variation of array flow pattern matrix on ratio scale estimated, what we more were concerned about is its entrained phase information.Through the FastICA algorithm; Obtain separation matrix
Figure 748900DEST_PATH_IMAGE051
; Albefaction matrix
Figure 418916DEST_PATH_IMAGE072
is known again; The simple operation that carries out matrix can obtain the estimation of array flow pattern, i.e.
Figure 736764DEST_PATH_IMAGE093
.
5. the reconstruct of type local-discharge ultrasonic array signal
Carry out the direction finding Position Research of array signal, need reconstruct array signal.Be the type local-discharge ultrasonic signal if confirmed
Figure 533950DEST_PATH_IMAGE047
road output signal
Figure 519224DEST_PATH_IMAGE057
; With the isolated component output matrix after other each road signal zero setting is
Figure 360141DEST_PATH_IMAGE058
, the array signal
Figure 165286DEST_PATH_IMAGE094
after then recovering.The process flow diagram of above-mentioned algorithm is shown in accompanying drawing 4.
(2) comprehensive evaluation of type local-discharge ultrasonic array signal Denoising Algorithm
Because the needs of the characteristics of array received signal and array location; Adopt ICA to carry out having obtained estimation again after its denoising to source signal and array flow pattern; Thereby also reconstruct the array received signal; Denoising appraisement system proposed by the invention is based on the theoretical research of array location, has adopted the relevant evaluating index parameter respectively to the estimated result of above three aspects.
1. based on the evaluation of source signal waveform similarity property
Above-mentioned waveform similarity coefficient can be estimated the similarity degree of denoising front and back waveform on the whole; But to the concussion situation of waveform, the evaluation of
Figure 749851DEST_PATH_IMAGE053
is incomplete.
Figure 120789DEST_PATH_IMAGE095
characterized the similarity degree of two wave form varies trend; Weighed the concussion situation of waveform to a certain extent;
Figure 601449DEST_PATH_IMAGE095
explains that more near 1 the concussion of waveform is similar more.
Figure 939896DEST_PATH_IMAGE096
(11)
Figure 531414DEST_PATH_IMAGE097
is the slope sum by the infinitesimal straight line of a waveform under ascendant trend for the ascendant trend parameter; Compare with the slope sum of the infinitely small straight line of another waveform under ascendant trend and to obtain, promptly
Figure 22438DEST_PATH_IMAGE098
(12)
Wherein, is original waveform;
Figure 984895DEST_PATH_IMAGE100
is the waveform after the denoising;
Figure 380104DEST_PATH_IMAGE099
all greater than
Figure 725635DEST_PATH_IMAGE101
,
Figure 548097DEST_PATH_IMAGE100
is all greater than
Figure 97022DEST_PATH_IMAGE102
.
Figure 295922DEST_PATH_IMAGE103
is the slope sum by the infinitesimal straight line of a waveform under downtrending for the downtrending parameter; Compare with the slope sum of the infinitely small straight line of another waveform under downtrending and to obtain, promptly
(13)
Wherein,
Figure 489323DEST_PATH_IMAGE105
.
2. the evaluation that distorts based on array flow pattern phase differential
The array positioning principle is to arrive " phase differential " that produced between each array element according to the type local-discharge ultrasonic signal to carry out DOA estimation, and these phase informations are embodied in the phase angle information of each element in its formed array flow pattern.So, after obtaining denoising, in the type local-discharge ultrasonic signal, also to realize the estimation of array flow pattern, the restructuring array signal just can position research.The phase angle difference of so-called array flow pattern is meant the phase angle difference between the array flow pattern component of array flow pattern component and reference array element of (except that reference array element) other array elements.The phase difference value of gathering denoising front and back array flow pattern compares, and weighs this Denoising Algorithm indirectly for keeping " phase differential " superiority that information had.
Define a vector
Figure 774811DEST_PATH_IMAGE106
Figure 777402DEST_PATH_IMAGE107
(14)
The middle component
Figure 18262DEST_PATH_IMAGE110
of array flow pattern
Figure 996210DEST_PATH_IMAGE109
of theoretical value and the phase angle difference between the reference component are got in expression in the formula; In like manner;
Figure 558965DEST_PATH_IMAGE111
array flow pattern
Figure 468015DEST_PATH_IMAGE112
for estimating,
Figure 68761DEST_PATH_IMAGE113
is sensor array element number.
Figure 63261DEST_PATH_IMAGE106
shown that directly the phase angle difference of two array flow patterns changes number percent situation shared in the standard phase angle difference; And the distortion situation of " phase differential " information before and after the explanation that can be indirect, denoising.Its value is more approaching zero, and it is little to explain that phase angle difference changes, and promptly " phase differential " distortion is little.
3. based on the evaluation of the array signal signal to noise ratio (S/N ratio) of reconstruct
The mixed signal that each array element of sensor receives is that the type local-discharge ultrasonic signal mixes with the unknown of noise signal, if in the mixed signal of each array element, the composition of noise is excessive, will directly influence the accuracy of type local-discharge ultrasonic location.Therefore, the signal to noise ratio (S/N ratio) of each mixed signal before and after the denoising as an independently denoising metric parameter, is estimated the superiority of denoising performance.
Input signal-to-noise ratio:
Figure 407655DEST_PATH_IMAGE114
(15)
Output signal-to-noise ratio:
Figure 171212DEST_PATH_IMAGE115
(16)
Figure 942859DEST_PATH_IMAGE116
; representes sampling number;
Figure 57893DEST_PATH_IMAGE118
representes the array received signal that mixed by pure type local-discharge ultrasonic signal respectively, the array received signal of reconstruct after the denoising and the array received signal that is mixed by type local-discharge ultrasonic signal and noise.
Based on the type local-discharge ultrasonic array signal wideband of array signal process technique focus on, arrowband direction finding and space orientation algorithm
(1) based on the broadband type local-discharge ultrasonic array signal focusing algorithm of TLS (Total Least Square) algorithm
The wideband array signal processing algorithm mainly is divided into two big types:
One type is based on irrelevant Signal Processing method (ISM).The main thought of ISM is with the narrow band data of broadband signal data decomposition to the nonoverlapping bands; On this basis; Each frequency band is carried out the subspace of narrow band signal and handle, estimate, obtain the final angle estimated result through certain combination again these initial estimation thereby obtain angle.ISM class algorithm operation quantity is big, and focusing effect is general.
Two types of disposal routes (CSM) that are based on coherent signal in addition.The basic thought of CSM is to focus on some reference frequency points to the signal space on the not overlapping Frequency point in the broadband signal frequency ranges of data, obtains the data covariance of reference frequency point, uses the method for narrow band signal processing again and carries out the DOA estimation.The CSM method compares with the ISM method that estimated accuracy is higher, operand is less relatively.
If the wideband array signal model is:
Figure 941536DEST_PATH_IMAGE119
(17)
In the formula,
Figure 884084DEST_PATH_IMAGE120
It is the reception data vector of array signal;
Figure 853177DEST_PATH_IMAGE121
The direction matrix, whole broadband be divided into into JIndividual arrowband;
Figure 539373DEST_PATH_IMAGE122
Data vector for original signal;
Figure 277522DEST_PATH_IMAGE123
Be that the j section receives data noise spatial data matrix.
Our purpose is through a focusing transform matrix
Figure 390972DEST_PATH_IMAGE124
, make corresponding to JThe data conversion of individual different frequency range becomes the data of a centre frequency.The basic thought of CSM class broadband focusing algorithm is exactly through focussing matrix the data of each Frequency point to be become the data of the existing Frequency point of reference, thereby forms correlation matrix, and its key is the selection of focussing matrix.
The present invention selects the broadband signal focusing algorithm based on total least square (TLS) for use.Compare with other algorithms, its superior performance, operand is little and noise inhibiting ability is strong.
Under ideal conditions, the array flow pattern after certain Frequency point focuses on should equate with the array flow pattern of reference frequency point fully, that is:
Figure 96628DEST_PATH_IMAGE125
(18)
The following formula brief note is for .Suppose that focussing matrix satisfies formula:
Figure 179171DEST_PATH_IMAGE127
(19)
As focussing matrix
Figure 197942DEST_PATH_IMAGE128
is carried out svd
Figure 141628DEST_PATH_IMAGE129
, then:
(20)
So
(21)
Suppose definition
Figure 619511DEST_PATH_IMAGE132
; Then focussing matrix satisfies
Figure 253754DEST_PATH_IMAGE133
, and formula (21) can be reduced to:
Figure 351023DEST_PATH_IMAGE134
(22)
By above-mentioned relation, on the definition total least square meaning cost function:
Figure 652692DEST_PATH_IMAGE135
(23)
In the formula,
Figure 810003DEST_PATH_IMAGE136
and
Figure 931543DEST_PATH_IMAGE137
is the perturbation matrix on the total least square meaning.
Find the solution formula (23), can get focussing matrix
Figure 612929DEST_PATH_IMAGE138
(24)
Wherein, the matrix that left singular vector that
Figure 769104DEST_PATH_IMAGE139
is
Figure 97317DEST_PATH_IMAGE140
respectively and right singular vector are formed.
The basic step of this algorithm is:
The initial value of
Figure 706153DEST_PATH_IMAGE006
estimated signal, and selected reference frequency point;
Figure 145225DEST_PATH_IMAGE008
utilizes initial value, constructs the array flow pattern of each Frequency point;
Figure 155906DEST_PATH_IMAGE012
utilizes formula (24) to construct the focussing matrix of each Frequency point;
Figure 655020DEST_PATH_IMAGE019
utilizes a series of focussing matrix array reception data to carry out focusing transform, obtains the data covariance matrix of single-frequency point;
Figure 33043DEST_PATH_IMAGE141
utilizes Estimation of Spatial Spectrum method estimated signal incident direction.
(2) the type local-discharge ultrasonic array signal direction-finding method that adopts phase matching principle and cloud optimization searching to combine
Suppose to arrive initial point OAcoustical signal do s( t), and the signal amplitude of arrival primitive is identical. n( t) be white noise.The reception signal of primitive is respectively so
Figure 744647DEST_PATH_IMAGE142
(25)
In the formula
Figure 140994DEST_PATH_IMAGE143
Be array element distance,
Figure 811009DEST_PATH_IMAGE144
Be the azimuth information of signal source,
Figure 394437DEST_PATH_IMAGE024
Be the velocity of sound, NLength for the signal data section.Signal frequency-domain was expressed as after formula (25) was done the conversion of Fu's formula
Figure 440891DEST_PATH_IMAGE145
(26)
Wherein is respectively the phase place of signal, wanted signal and the noise of reception.
According to the signal phase matching principle (its principle is for to carry out phase compensation to the output of each array element, the phase differential that aspect causes if the phase place of compensation equals, then the direction of arrival of each array element is identical, just the space phase of signal is mated.), when having only the signal space phase place identical, just can find the solution signal by formula (26) SReal part and imaginary part.Therefore use
Figure 74988DEST_PATH_IMAGE147
respectively, multiplier (26) both sides get
Figure 145713DEST_PATH_IMAGE148
?(27)
Its matrix form is:
Figure 730278DEST_PATH_IMAGE149
(28)
Wherein
Figure 570058DEST_PATH_IMAGE150
(29)
Figure 581876DEST_PATH_IMAGE151
(30)
Figure 139896DEST_PATH_IMAGE152
(31)
Figure 262573DEST_PATH_IMAGE153
(32)
Figure 238751DEST_PATH_IMAGE154
(33)
At last; Through the corresponding angle of DOA estimation criterion
Figure 155891DEST_PATH_IMAGE155
search
Figure 466787DEST_PATH_IMAGE156
maximal value, can the office of obtaining put the azimuth information in source.Wherein Tr (.) is defined as matrix trace, promptly is defined as the diagonal element sum of matrix.
For avoiding traditional method of exhaustion maximizing arithmetic speed shortcoming slowly, the present invention adopts the cloud optimization searching algorithm, and (Cloud Theory Optimization Algorithm CTOA) carries out the search of optimum solution.Cloud model has the characteristics of " have uncertainty in confirming, change again in stable " when expressing knowledge, demonstrated fully " spore " ultimate principle.For cloud model
Figure 861996DEST_PATH_IMAGE157
; The excellent characteristic of
Figure 207527DEST_PATH_IMAGE158
representative " parent " individual inheritance; It is " filial generation " succession to " parent "; The uncertainty of
Figure 29989DEST_PATH_IMAGE158
and
Figure 828181DEST_PATH_IMAGE159
expression genetic process has showed " variation " characteristic in the spore process.Theoretical according to cloud; Normal state cloud operator can be accomplished the conversion of concept space to numerical space; This conversion is to confirm with accurate on the one hand, because each water dust of numerical space all is the once quantification realization of qualitativing concept, all is the representative of this qualitativing concept to a certain extent; On the other hand, this conversion is at random and blurs that conversion each time obtains different water dust set again, any one representative during same qualitativing concept can be gathered with water dust, and also different water dusts is represented definite degree difference of this notion.
Therefore; Can fully use for reference the evolution thought of genetic algorithm " survival of the fittest; the survival of the fittest "; The individual excellent characteristic of representative " parent " with ; Use
Figure 476386DEST_PATH_IMAGE161
and the degree of
Figure 469750DEST_PATH_IMAGE159
control of heredity and variation; Accomplish the conversion of concept space with normal state cloud operator to numerical space; Produce population; Carry out genetic manipulation, realizing plays a game puts the search of source space coordinate.
Realize that cloud optimization calculationization needs several steps below the process: 1) select the elite individual; 2) confirm the excellent individual vector; 3) the initialization population that evenly distributes; 4) select; 5) evolve; 6) variation.Its process flow diagram is shown in accompanying drawing 5.
(3) put the source location algorithm based on the office of global optimization search
At first utilize 3 groups of supersonic array sensors to carry out Partial Discharge Detection at the diverse location of electrical equipment; Obtain the position angle and the angle of pitch of the direction finding line of every group pattern sensor; And set up the equation that certain spatial is put each bar direction finding linear distance sum based on the home position of sensor array, look the locus that the source is put in its spatial point be office corresponding apart from the sum minimum value.Asking for of minimum value is based upon based on the global search of revising genetic algorithm.
This method has taken into full account the non-intersect actual conditions of antarafacial of direction finding line in experimental study and the engineering survey and the influence of various errors, and effectively the bearing accuracy in source is put in raising office.Concrete method is derived as follows.
Coupling place of if
Figure 755237DEST_PATH_IMAGE162
group pattern sensor (establishing each sensor array) coordinate is
Figure 492249DEST_PATH_IMAGE163
,
Figure 546793DEST_PATH_IMAGE164
.The deflection and the angle of pitch of the sensor array that then can obtain being coupled through the direction finding algorithm are respectively .On the basis of the position of direction finding angle and sensor array, can confirm the straight-line equation at direction finding line place.
If the direction vector of this line
Figure 749421DEST_PATH_IMAGE166
is
Figure 24545DEST_PATH_IMAGE167
; Suppose the local discharge of electrical equipment position again for
Figure 949907DEST_PATH_IMAGE168
, having no under the situation of error the office of intersecting to put source for many direction of arrival straight lines of diverse location.But in experimental study and the engineering survey because the influence of error;
Figure 545153DEST_PATH_IMAGE162
bar straight line is a different surface beeline, can the office of obtaining put minimum and the function of source to the straight line at each bar direction finding line place thus.It is:
Figure 623967DEST_PATH_IMAGE169
(34)
Within the boundary space condition, carry out the global optimizing search; Seeking certain point
Figure 653103DEST_PATH_IMAGE170
makes the space put the minimum apart from sum of each direction finding line; Promptly
Figure 424750DEST_PATH_IMAGE171
value is minimum, promptly looks this point and puts the source space coordinate for office.Principle is shown in accompanying drawing 6.
Shown in accompanying drawing 6; If the direction vector of the straight line at the direction finding line of sensor array
Figure 906547DEST_PATH_IMAGE172
place is
Figure 569479DEST_PATH_IMAGE173
, and its position angle, the angle of pitch are
Figure 453121DEST_PATH_IMAGE174
.The direction vector of XOY plane normal is
Figure 395669DEST_PATH_IMAGE175
, and the expression formula of direction vector is again:
Figure 364762DEST_PATH_IMAGE176
(35)
The straight-line equation that therefore can get the direction finding line place of sensor array
Figure 785379DEST_PATH_IMAGE172
is:
Figure 789107DEST_PATH_IMAGE177
(36)
The coupling position of sensor array
Figure 902557DEST_PATH_IMAGE178
; Shelf depreciation power supply locus is , get final product thus shelf depreciation power supply to the straight line at each direction finding line place following apart from expression formula:
(37)
Wherein
Figure 926642DEST_PATH_IMAGE181
.Therefore, the straight line that this moment, the shelf depreciation power supply belonged to each direction finding line apart from sum be:
Figure 210992DEST_PATH_IMAGE182
(38)
Can estimate the volume coordinate of the shelf depreciation power supply of minimum
Figure 154678DEST_PATH_IMAGE183
and this moment through global search.
The present invention adopts a kind of global search orientation problem that genetic algorithm solves the local discharge of electrical equipment power supply of revising.Compare with traditional genetic algorithm, this algorithm can effectively be avoided local optimal searching and precocious convergence, and fast convergence rate, and its novelty mainly contains following three aspects: the one, and correction algorithm is taked encoding mechanism technology and the method that produces initial population at random; The 2nd, correction algorithm is in order to stablize diversity individual in the population and to enlarge the search volume, the adaptive algorithm that algorithm use variation probability changes with fitness automatically; The 3rd, correction algorithm is deleted relatively poor solution space in evolution; And by deletion group size breeding better solutions space; This method has strengthened the density of excellent solution space and defect individual thereof, makes evolutionary process avoid local optimal searching with local precocious, and has improved the search capability and the speed of convergence of algorithm.Its basic step is following:
coding: the genotype string structure data that the feasible solution tables of data of solution space are shown as hereditary space;
Figure 895418DEST_PATH_IMAGE008
generates initial population: produce N original string structured data at random; Each string structure data is called body one by one, and this N string structure data are the primary iteration data;
Figure 85091DEST_PATH_IMAGE012
assessment detects fitness value;
selects: from current colony, select good individuality by the adaptable individual big principle of probability that produces one or more offsprings, make it have an opportunity to be breeding descendants of future generation as parent; Selection course has embodied the thought of revising genetic algorithm, has realized Darwinian survival of the fittest principle;
Figure 331450DEST_PATH_IMAGE141
intersects: topmost genetic manipulation is interlace operation in the correction genetic algorithm.Can obtain inheriting the offspring individual of the characteristic of its former generation's individuality through it, intersection has embodied the thought of message exchange;
Figure 164277DEST_PATH_IMAGE184
variation: variation is in colony, to select body one by one to change the value of certain displacement in the string structure data with certain probabilistic ground at random, and making a variation provides chance for new individual generation.
Figure 524851DEST_PATH_IMAGE185
is optimum to be preserved: for the enhancement algorithms stability and convergence; In new algorithm, make the excellent individual of some in the parent directly get into follow-on operation; Can prevent like this to duplicate, intersect or make a variation in the destruction of excellent individual, make optimizing process global convergence on probabilistic.Revise the genetic algorithm flow process shown in accompanying drawing 7.
Variable GEN is current evolutionary generation in the process flow diagram, and N is a population scale, and M is the maximum algebraically that algorithm is carried out.In revising genetic algorithm, group size (N) is generally got 20-100; Stop evolutionary generation (M) and get 100-500; Crossover probability (
Figure 443128DEST_PATH_IMAGE186
) is got 0.4-0.99%; Variation probability (
Figure 812930DEST_PATH_IMAGE187
) is got 0.1%-1%.
Through above-mentioned steps, can realize the accurate location of local discharge of electrical equipment.
3. advantage of the present invention
(1) nine yuan of cruciform supersonic arrays of the instant-plugging sensor that designs among the present invention; Have the control of wave beam flexibly, high signal gain, the higher spatial resolving power except that having the planar square sensor array; Also can reduce redundant array element to a certain extent, reduce cost;
(2) the type local-discharge ultrasonic array positioning system adopts existing typing instrument and equipment among the present invention, on market, can directly buy, and can realize the online detection to electrical equipment, and method is simple, and can not cause equipment under test impaired, and positional accuracy is high.
(3) employing can effectively focus on the type local-discharge ultrasonic broadband signal based on the broadband signal focusing algorithm of TLS in the localization method of the present invention, and operand is little.
(4) the type local-discharge ultrasonic array signal direction-finding method that adopts phase matching principle and cloud optimization searching to combine in the localization method among the present invention, direction finding precision is high, fast convergence rate.
(5) employing is put the source location method based on three platform direction findings with the office of revising genetic algorithm in the localization method of the present invention, and bearing accuracy is high, fast convergence rate.
4, checking of the present invention
(1), verifies based on the focusing of type local-discharge ultrasonic array signal wideband, arrowband direction finding and the space-location method of array signal process technique
1. obtain one group of measured signal based on this system
Local discharge of electrical equipment positioning system shown in the accompanying drawing 3 array element distance of nine yuan of cruciform supersonic array sensors (select 5mm) is placed the laboratory of abundant simulated field noisy environment, and the volume coordinate of at first setting discharge source is (30,120; 40) cm; The sensor array volume coordinate is (50,0,0) cm; The theoretical position angle and the angle of pitch that can calculate the relative sensor array of this position discharge source thus are (99.5 °, 18.2 °); Setting image data length then is 8000, and 256 times of signal amplification factor, SF are 10MHz, and filter range 60kHz to 300kHz, triggering mode are external trigger.To record 9 noisy type local-discharge ultrasonic array signals of passage as shown in Figure 8 in experiment, transverse axis express time wherein, and unit is ms, and the longitudinal axis is represented voltage magnitude, and unit is V.
2. focus on and the direction finding result
Adopt broadband signal focusing algorithm proposed by the invention and narrow band signal direction finding algorithm; Can be (105.6 ° by the direction finding result that the local discharge of electrical equipment supersonic array signal that above-mentioned experiment records is obtained this position; 23.3 °), with the error at theoretical direction finding angle be (6.1 °, 5.1 °).Direction finding spectrogram and contour map are shown in Fig. 9 and accompanying drawing 10.
3. locate
The present invention is based on three platform direction finding principles and revise the genetic algorithm source of putting of playing a game and position, on above-mentioned direction finding result's basis, is (0,50 at coordinate respectively; 0) two position repeated experiments processes of cm, (0,100,0) cm, its direction finding result is respectively (71.8 °; 33.5 °), (39.8 °, 39.8 °), with theoretical value (66.8 °, 27.7 °), (33.7 °; 48.0 °) error be (5.0 °, 5.8 °), (6.1 °, 8.2 °).
According to aforementioned theory, adopt based on the office of revising genetic algorithm and put the source location method, can the office of obtaining putting the source space coordinate is (23.5,126.8,46.3) cm, positioning error is 9.8cm.
4. error analysis
Be the correctness of checking partial discharge positioning method that the present invention puies forward and system, other 3 groups of positioning results are as shown in table 1.
Can find out through table 1; The shelf depreciation array localization method and the system that utilize the present invention to carry position shelf depreciation; The positioning error of experiment has certain engineering practicability all about 11 cm several times, but because experiment is in the noisy environment of simulated field, to carry out; The type local-discharge ultrasonic signal that collects is aliasing in noise signal, makes the precision of Partial Discharge Detection descend greatly.
Table 1 three platform direction finding positioning experiment results
Figure 234684DEST_PATH_IMAGE188
(2), verify based on the type local-discharge ultrasonic array signal denoising method of quick independent component analysis
1. carry out denoising, reconstruction signal to above-mentioned reception signal
The signal that obtains after adopting algorithm that the present invention proposes to the signal denoising that collects among Fig. 8 is shown in figure 11, transverse axis express time wherein, and unit is ms, and the longitudinal axis is represented voltage magnitude, and unit is V.
Array received signal by after the estimation reconstruct acquisition denoising of Figure 11 type local-discharge ultrasonic signal and array flow pattern is shown in figure 12, and wherein transverse axis is the time, and unit is ms, and the longitudinal axis is a voltage magnitude, and unit is V.Compare with Fig. 8, the signal to noise ratio (S/N ratio) of each channel signal significantly improves.
2. the focusing of signal and direction finding result after the denoising
Shelf depreciation supersonic array signal after adopting broadband signal focusing algorithm proposed by the invention and narrow band signal direction finding algorithm to denoising reconstruct is handled; The direction finding result that can try to achieve this position is (102.6 °, 21.6 °), with the error at theoretical direction finding angle be (3.1 °; 3.4 °); And the error at the direction finding angle of original signal and theoretical direction finding angle is (6.1 °, 5.1 °), and the spatial information (si) of supersonic array signal has obtained preserving greatly after the visible denoising.The direction finding spectrogram of signal and contour map such as Figure 13 and shown in Figure 14 after the denoising.Compare with Figure 10 with Fig. 9, its spectrum peak is more sharp-pointed, secondary lobe is littler, directive property is better.
3. locate
Employing is (0,50,0) cm, (0,100 to coordinate respectively based on the type local-discharge ultrasonic array signal Denoising Algorithm of quick independent component analysis; 0) the supersonic array signal of gathering on two of cm positions carries out denoising, then it is carried out direction finding, and the result is respectively (70.3 °; 27.5 °), (35.4 °, 46.6 °), with (66.8 ° of theoretical values; 27.7 °), the error of (33.7 °, 48.0 °) be (3.5 °, 0.2 °), (1.7 °; 1.4 °), to compare with the direction finding result of original signal, precision obviously improves.
According to aforementioned theory, adopt based on the office of revising genetic algorithm and put the source location method, can the office of obtaining putting the source space coordinate is (26.8,122.5,41.7) cm, positioning error is 4.4cm.
4. error analysis
The detection position of source and sensor array is put in change office, and the array signal before and after the denoising is carried out three direction finding Position Research, and the result is as shown in table 2.
Can find out that from table 2 result the FastICA algorithm that utilizes the present invention to propose carries out the absolute error of positioning result after the denoising in 6cm, and without the absolute error of the positioning result of denoising about 11cm.Hence one can see that, thereby the noise that can effectively remove in the type local-discharge ultrasonic array signal based on the Denoising Algorithm of FastICA improves bearing accuracy, it is thus clear that the necessity and the feasibility of Denoising Algorithm proposed by the invention.
Table 2 three platform direction finding positioning experiment results
Figure 31739DEST_PATH_IMAGE189

Claims (7)

1. a local discharge of electrical equipment localization method is characterized in that, it at first utilizes three supersonic array sensors to receive the shelf depreciation ultrasonic signal at the diverse location of electrical equipment, and forms Array Model; The wideband focusing algorithm of using then based on TLS converts the broadband signal that three supersonic array sensors are received to narrow band signal; The azimuth information in source is put in the arrowband direction finding algorithm acquisition office that adopts phase matching principle and cloud optimization searching to combine again, adopts the office of global optimization search to put the particular location that the source location algorithm is confirmed local discharge of electrical equipment at last.
2. according to the said local discharge of electrical equipment localization method of claim 1, it is characterized in that it may further comprise the steps:
A, utilize three supersonic array sensors to receive the shelf depreciation ultrasonic signal at the diverse location of electrical equipment;
B, the wideband focusing algorithm of using based on TLS convert the broadband signal that three supersonic array sensors are received to narrow band signal:
Suppose that the wideband array signal model is:
Figure 653090DEST_PATH_IMAGE001
;
In the formula,
Figure 425874DEST_PATH_IMAGE002
It is the reception data vector of array signal;
Figure 700997DEST_PATH_IMAGE003
The direction matrix, whole broadband be divided into into JIndividual arrowband;
Figure 875627DEST_PATH_IMAGE004
Data vector for original signal;
Figure 210793DEST_PATH_IMAGE005
Be that the j section receives data noise spatial data matrix,
The basic step of conversion is:
The initial value of
Figure 470873DEST_PATH_IMAGE006
estimated signal, and selected reference frequency ;
Figure 329556DEST_PATH_IMAGE008
utilizes the initial value of signal, constructs the array flow pattern
Figure 101203DEST_PATH_IMAGE009
of each Frequency point;
constructs the focussing matrix of each Frequency point:
Figure 465505DEST_PATH_IMAGE011
; Wherein,
Figure 349147DEST_PATH_IMAGE012
, the matrix that left singular vector that
Figure 291696DEST_PATH_IMAGE013
and
Figure 260789DEST_PATH_IMAGE014
is respectively and right singular vector are formed;
Figure 934401DEST_PATH_IMAGE016
utilizes a series of focussing matrix array reception data to carry out focusing transform, obtains the data covariance matrix of single-frequency point;
The azimuth information in source is put in c, the arrowband direction finding algorithm acquisition office that adopts phase matching principle and cloud optimization searching to combine:
Suppose to arrive initial point OAcoustical signal do
Figure 47851DEST_PATH_IMAGE017
, and the signal amplitude of arrival primitive is identical, n( t) be white noise, the reception signal of primitive is respectively so:
Figure 504240DEST_PATH_IMAGE018
In the formula, ,
Figure 321203DEST_PATH_IMAGE020
Be array element distance,
Figure 136713DEST_PATH_IMAGE021
Be the azimuth information of signal source,
Figure 549239DEST_PATH_IMAGE022
Be the velocity of sound, NBe the length of signal data section, signal frequency-domain was expressed as after following formula was done the conversion of Fu's formula:
Figure 859129DEST_PATH_IMAGE023
Its matrix form is:
Figure 40712DEST_PATH_IMAGE024
Adopt the corresponding angle of maximal value of cloud optimization searching algorithm search DOA estimation criterion
Figure 761543DEST_PATH_IMAGE025
; The azimuth information of can the office of obtaining putting the source; Be the position angle and the angle of pitch of direction finding line
Wherein, Tr (.) is the diagonal element sum of matrix,
Figure 661366DEST_PATH_IMAGE026
D, employing are put the particular location that the source location algorithm is confirmed local discharge of electrical equipment based on the office of revising genetic algorithm global optimization search:
At first set up the equation that certain spatial is put each bar direction finding linear distance sum according to the home position of position angle, the angle of pitch and the sensor array of the direction finding line of three supersonic array sensors; Adopt then based on revise the genetic algorithm global search its apart from the sum minimum value, this is worth the locus that the source is put in the pairing spatial point office of being.
3. according to claim 1 or 2 said local discharge of electrical equipment localization methods; It is characterized in that; Before the broadband signal that the supersonic array sensor is received focuses on, should adopt algorithm that signal is carried out denoising based on quick independent component analysis (FastICA), concrete steps are:
1. play a game and put the supersonic array signal and carry out centralization (going average) and whitening pretreatment:
Centralization: each group observation data
Figure 758635DEST_PATH_IMAGE027
thus deduct its average to make it become an average be 0 variable;
Albefaction: the covariance matrix to mixing observation signal
Figure 60303DEST_PATH_IMAGE028
carries out characteristic value decomposition; Obtain feature matrix
Figure 952036DEST_PATH_IMAGE029
and be the diagonal matrix of diagonal element with the eigenwert; Then the signal after the albefaction is
Figure 489383DEST_PATH_IMAGE031
, wherein
Figure 379978DEST_PATH_IMAGE032
;
2. the blind separation of type local-discharge ultrasonic array signal,
Expectation function based on the FastICA algorithm of plural number is:
Figure 317027DEST_PATH_IMAGE034
Figure 21678DEST_PATH_IMAGE035
is separating vector in the formula, several kinds of functions below function can be selected:
Figure 627606DEST_PATH_IMAGE038
Figure 621101DEST_PATH_IMAGE039
;
Figure 156304DEST_PATH_IMAGE041
all is approximately 0.1; The maximum value of calculating target function, the two step formulas that obtain plural FastICA are:
Figure 520607DEST_PATH_IMAGE043
In the formula; respectively representative function
Figure 346797DEST_PATH_IMAGE045
once reach second derivative
When having a plurality of independent source signal; For fear of the same source signal of repeated isolation; Before
Figure 417521DEST_PATH_IMAGE046
individual vector is separated;
Figure 251354DEST_PATH_IMAGE047
for iteration generation each time, carry out following operation:
Figure 622292DEST_PATH_IMAGE048
In the formula, the separation matrix that expression is made up of preceding
Figure 926552DEST_PATH_IMAGE050
individual separating vector;
3. the extraction of type local-discharge ultrasonic signal:
Calculate waveform similarity coefficient (
Figure 783649DEST_PATH_IMAGE051
):
representes that respectively the Bureau of Standards puts the type local-discharge ultrasonic signal of ultrasonic signal and extraction, extracts the maximum signal of absolute value wherein as the type local-discharge ultrasonic signal;
4. the array flow pattern of type local-discharge ultrasonic signal is estimated:
Being estimated as of array flow pattern:
Figure 987863DEST_PATH_IMAGE054
;
5. the reconstruct of type local-discharge ultrasonic array signal:
Be the type local-discharge ultrasonic signal if confirmed
Figure 914231DEST_PATH_IMAGE046
road output signal ; With the isolated component output matrix after other each road signal zero setting is
Figure 82224DEST_PATH_IMAGE056
, the array signal
Figure 349257DEST_PATH_IMAGE057
after then recovering.
4. local discharge of electrical equipment positioning system; It is characterized in that; It comprises supersonic array sensor (15), multi-channel synchronous data acquisition device (16) and computing machine (17); Said supersonic array sensor (15) comprises Stainless Steel Shell (8) and is installed in a plurality of piezoelectricity array elements (7) in the Stainless Steel Shell (8); Said a plurality of piezoelectricity array elements (7) are arranged in cruciform, and their signal receiving end is towards the signal receiving port of Stainless Steel Shell (8), and the signal output part of a plurality of piezoelectricity array elements (7) connects the different input ends of multi-channel synchronous data acquisition device (16) respectively through array element lead-in wire (4); The signal output part of said multi-channel synchronous data acquisition device (16) connects computing machine (17).
5. according to the said local discharge of electrical equipment positioning system of claim 4; It is characterized in that; The signal receiving port of said Stainless Steel Shell (8) is provided with stainless steel cuticula (10), between Stainless Steel Shell (8) and piezoelectricity array element (7), is provided with silicon rubber filling thing (9).
6. according to claim 4 or 5 said local discharge of electrical equipment positioning systems, it is characterized in that, between piezoelectricity array element (7) and array element lead-in wire (4), be provided with acoustic damping piece (6).
7. according to the said local discharge of electrical equipment positioning system of claim 6, it is characterized in that said piezoelectricity array element (7) is provided with nine.
CN201210173231.3A 2012-05-30 2012-05-30 A kind of local discharge of electrical equipment localization method and system Expired - Fee Related CN102662132B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210173231.3A CN102662132B (en) 2012-05-30 2012-05-30 A kind of local discharge of electrical equipment localization method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210173231.3A CN102662132B (en) 2012-05-30 2012-05-30 A kind of local discharge of electrical equipment localization method and system

Publications (2)

Publication Number Publication Date
CN102662132A true CN102662132A (en) 2012-09-12
CN102662132B CN102662132B (en) 2015-09-09

Family

ID=46771661

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210173231.3A Expired - Fee Related CN102662132B (en) 2012-05-30 2012-05-30 A kind of local discharge of electrical equipment localization method and system

Country Status (1)

Country Link
CN (1) CN102662132B (en)

Cited By (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103529364A (en) * 2013-09-20 2014-01-22 华北电力大学(保定) In-oil multi-local discharge source positioning method based on ultrasonic thinned array sensor
CN103760479A (en) * 2014-01-10 2014-04-30 三峡大学 Portable corona positioning instrument
CN104251961A (en) * 2013-06-27 2014-12-31 国家电网公司 Partial discharging signal source positioning device and system of gas-insulation completely-enclosed combined electrical appliance
CN104808124A (en) * 2015-04-30 2015-07-29 国家电网公司 Acoustic-electric combined positioning method of switch cabinet local discharge fault positions
CN105719003A (en) * 2016-01-13 2016-06-29 中国南方电网有限责任公司超高压输电公司 Quantum genetic algorithm-based converter transformer partial-discharge ultrasonic location method
CN105911497A (en) * 2016-06-15 2016-08-31 国网北京市电力公司 Correction method and device for partial discharge positioning error
CN105929312A (en) * 2016-04-28 2016-09-07 华北电力大学(保定) Electrical equipment partial discharge source direction finding apparatus
CN106019098A (en) * 2016-06-15 2016-10-12 国网北京市电力公司 Partial discharge positioning method, device and system
CN106556783A (en) * 2016-12-05 2017-04-05 西安交通大学 A kind of shelf depreciation direction-finding method in transformer station based on superfrequency phased array principle
CN107271860A (en) * 2017-06-09 2017-10-20 北京兴迪仪器有限责任公司 Partial discharge pulse waveform similarity recognition method
CN107677941A (en) * 2017-09-30 2018-02-09 北京华电智成电气设备有限公司 A kind of high-frequency noise synthesizes synchronous common mode noise-reduction method
CN107741554A (en) * 2017-09-30 2018-02-27 北京华电智成电气设备有限公司 A kind of spatial noise synthesizes synchronization combining noise-reduction method
CN107942212A (en) * 2017-11-17 2018-04-20 国网天津市电力公司 A kind of substation's partial discharge positioning method without blur estimation based on spatial spectrum
CN107966206A (en) * 2017-11-07 2018-04-27 华北电力大学(保定) A kind of array signal denoising method of fiber ultrasonic sensor array
EP3234981A4 (en) * 2014-12-17 2018-07-11 Lockheed Martin Corporation Sensor array packaging solution
CN108445447A (en) * 2018-02-27 2018-08-24 国家电网有限公司 A kind of station domain space Mutual coupling system of substation's discharge source
CN109116202A (en) * 2018-10-17 2019-01-01 江苏方天电力技术有限公司 A kind of PVDF ultrasound array sensor and its array signal denoising method
CN109444691A (en) * 2018-11-22 2019-03-08 天津大学 Multiple physical field based on ultrasonic method couples silicon rubber partial discharge detecting system
CN109829416A (en) * 2019-01-28 2019-05-31 国网山东省电力公司日照供电公司 A kind of method and system detecting ultrahigh frequency partial discharge signal from random noise
CN110470956A (en) * 2019-08-05 2019-11-19 上海电机学院 A kind of power equipment shelf depreciation ultrasound locating method
CN110865287A (en) * 2019-11-27 2020-03-06 国网重庆市电力公司电力科学研究院 GIS partial discharge positioning device
CN111060871A (en) * 2019-12-26 2020-04-24 南京长峰航天电子科技有限公司 Five-element array positioning method and device based on improved genetic algorithm
CN111289861A (en) * 2020-03-26 2020-06-16 云南电网有限责任公司电力科学研究院 Method for detecting position of partial discharge source
CN111428596A (en) * 2020-03-16 2020-07-17 重庆邮电大学 Grinding sound signal detection method based on three sound pickups
CN111665422A (en) * 2020-06-08 2020-09-15 郑州精铖电力设备有限公司 FPGA-based microphone array non-invasive type broadband sound wave real-time imaging detection system
CN112255507A (en) * 2020-09-10 2021-01-22 深圳供电局有限公司 Partial discharge positioning method and device, computer equipment and storage medium
CN113361607A (en) * 2021-06-08 2021-09-07 云南电网有限责任公司电力科学研究院 Medium-voltage distribution network line problem analysis method and device
CN114280436A (en) * 2021-12-24 2022-04-05 中国科学院电工研究所 F-P ultrasonic sensor array implantation device for monitoring partial discharge of power equipment
CN114966547A (en) * 2022-05-18 2022-08-30 珠海视熙科技有限公司 Compensation method, system and device for improving sound source positioning precision
CN117347797A (en) * 2023-09-28 2024-01-05 中国长江电力股份有限公司 Ultrasonic array detection device and method for partial discharge

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101702000A (en) * 2009-11-26 2010-05-05 华北电力大学(保定) Phased ultrasound locating method and system of partial discharge of transformer
CN102183710A (en) * 2011-03-02 2011-09-14 华北电力大学(保定) Method and system for identifying numbers of partial discharge sources of transformers

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101702000A (en) * 2009-11-26 2010-05-05 华北电力大学(保定) Phased ultrasound locating method and system of partial discharge of transformer
CN102183710A (en) * 2011-03-02 2011-09-14 华北电力大学(保定) Method and system for identifying numbers of partial discharge sources of transformers

Non-Patent Citations (11)

* Cited by examiner, † Cited by third party
Title
SHAHROKH VALAEE ET AL.: "Localization of Wideband Signals Using Least-Squares and Total Least-Squares Approaches", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING,1999》 *
乔梁 等: "扩展的快速独立分量分析算法", 《哈尔滨工程大学学报》 *
孙进才 等: "利用信号相位匹配原理的声源方位估计实验研究术", 《声学学报》 *
孙进才 等: "平面阵的方位估计算法及其误差分析", 《声学学报》 *
李楠 等: "基于SGA与LS-SVM的变压器油中溶解气体含量预测", 《电力设备监测与故障诊断》 *
肖卉 等: "基于信号相位匹配原理的二维方位估计方法", 《信号处理》 *
谢庆 等: "基于多平台测向及全局搜索的局部放电超声阵列定位方法", 《电工技术学报》 *
谢庆 等: "超声相控阵与宽带阵列信号处理相结合的油中局部放电定位方法", 《中国电机工程学报》 *
谢庆: "基于空间谱估计的变压器局放超声阵列定位方法研究", 《中国优秀博士学位论文全文数据库工程科技II辑》 *
贾银洁 等: "基于FastICA的混合音频信号盲分离", 《信息与电子工程》 *
顾江 等: "基于独立分量分析的声发射信号去噪方法", 《江南大学学报(自然科学版)》 *

Cited By (43)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104251961A (en) * 2013-06-27 2014-12-31 国家电网公司 Partial discharging signal source positioning device and system of gas-insulation completely-enclosed combined electrical appliance
CN103529364B (en) * 2013-09-20 2016-05-25 华北电力大学(保定) Many Partial Discharge Sources localization method in oil based on ultrasonic thinned array sensor
CN103529364A (en) * 2013-09-20 2014-01-22 华北电力大学(保定) In-oil multi-local discharge source positioning method based on ultrasonic thinned array sensor
CN103760479A (en) * 2014-01-10 2014-04-30 三峡大学 Portable corona positioning instrument
CN103760479B (en) * 2014-01-10 2016-11-02 三峡大学 Portable corona position finder
EP3234981A4 (en) * 2014-12-17 2018-07-11 Lockheed Martin Corporation Sensor array packaging solution
CN104808124A (en) * 2015-04-30 2015-07-29 国家电网公司 Acoustic-electric combined positioning method of switch cabinet local discharge fault positions
CN105719003A (en) * 2016-01-13 2016-06-29 中国南方电网有限责任公司超高压输电公司 Quantum genetic algorithm-based converter transformer partial-discharge ultrasonic location method
CN105929312B (en) * 2016-04-28 2019-03-26 华北电力大学(保定) A kind of electrical equipment partial discharge source direction-finding device
CN105929312A (en) * 2016-04-28 2016-09-07 华北电力大学(保定) Electrical equipment partial discharge source direction finding apparatus
CN106019098A (en) * 2016-06-15 2016-10-12 国网北京市电力公司 Partial discharge positioning method, device and system
CN105911497A (en) * 2016-06-15 2016-08-31 国网北京市电力公司 Correction method and device for partial discharge positioning error
CN106556783B (en) * 2016-12-05 2019-07-19 西安交通大学 Shelf depreciation direction-finding method based on superfrequency phased array principle in a kind of substation
CN106556783A (en) * 2016-12-05 2017-04-05 西安交通大学 A kind of shelf depreciation direction-finding method in transformer station based on superfrequency phased array principle
CN107271860A (en) * 2017-06-09 2017-10-20 北京兴迪仪器有限责任公司 Partial discharge pulse waveform similarity recognition method
CN107271860B (en) * 2017-06-09 2019-09-24 北京兴迪仪器有限责任公司 Partial discharge pulse waveform similarity recognition method
CN107677941A (en) * 2017-09-30 2018-02-09 北京华电智成电气设备有限公司 A kind of high-frequency noise synthesizes synchronous common mode noise-reduction method
CN107741554A (en) * 2017-09-30 2018-02-27 北京华电智成电气设备有限公司 A kind of spatial noise synthesizes synchronization combining noise-reduction method
CN107677941B (en) * 2017-09-30 2019-12-13 北京华电智成电气设备有限公司 high-frequency noise synthesis synchronous common-mode noise reduction method
CN107741554B (en) * 2017-09-30 2019-12-13 北京华电智成电气设备有限公司 Synchronous combined noise reduction method for airspace noise synthesis
CN107966206B (en) * 2017-11-07 2019-01-04 华北电力大学(保定) A kind of array signal denoising method of fiber ultrasonic sensor array
CN107966206A (en) * 2017-11-07 2018-04-27 华北电力大学(保定) A kind of array signal denoising method of fiber ultrasonic sensor array
CN107942212A (en) * 2017-11-17 2018-04-20 国网天津市电力公司 A kind of substation's partial discharge positioning method without blur estimation based on spatial spectrum
CN108445447A (en) * 2018-02-27 2018-08-24 国家电网有限公司 A kind of station domain space Mutual coupling system of substation's discharge source
CN108445447B (en) * 2018-02-27 2021-09-28 国家电网有限公司 Substation area space direction of arrival estimation system of substation discharge source
CN109116202A (en) * 2018-10-17 2019-01-01 江苏方天电力技术有限公司 A kind of PVDF ultrasound array sensor and its array signal denoising method
CN109444691A (en) * 2018-11-22 2019-03-08 天津大学 Multiple physical field based on ultrasonic method couples silicon rubber partial discharge detecting system
CN109829416A (en) * 2019-01-28 2019-05-31 国网山东省电力公司日照供电公司 A kind of method and system detecting ultrahigh frequency partial discharge signal from random noise
CN110470956A (en) * 2019-08-05 2019-11-19 上海电机学院 A kind of power equipment shelf depreciation ultrasound locating method
CN110865287A (en) * 2019-11-27 2020-03-06 国网重庆市电力公司电力科学研究院 GIS partial discharge positioning device
CN111060871A (en) * 2019-12-26 2020-04-24 南京长峰航天电子科技有限公司 Five-element array positioning method and device based on improved genetic algorithm
CN111428596A (en) * 2020-03-16 2020-07-17 重庆邮电大学 Grinding sound signal detection method based on three sound pickups
CN111289861A (en) * 2020-03-26 2020-06-16 云南电网有限责任公司电力科学研究院 Method for detecting position of partial discharge source
CN111289861B (en) * 2020-03-26 2022-01-25 云南电网有限责任公司电力科学研究院 Method for detecting position of partial discharge source
CN111665422A (en) * 2020-06-08 2020-09-15 郑州精铖电力设备有限公司 FPGA-based microphone array non-invasive type broadband sound wave real-time imaging detection system
CN112255507A (en) * 2020-09-10 2021-01-22 深圳供电局有限公司 Partial discharge positioning method and device, computer equipment and storage medium
CN112255507B (en) * 2020-09-10 2022-07-22 深圳供电局有限公司 Partial discharge positioning method and device, computer equipment and storage medium
CN113361607A (en) * 2021-06-08 2021-09-07 云南电网有限责任公司电力科学研究院 Medium-voltage distribution network line problem analysis method and device
CN113361607B (en) * 2021-06-08 2023-01-20 云南电网有限责任公司电力科学研究院 Medium-voltage distribution network line problem analysis method and device
CN114280436A (en) * 2021-12-24 2022-04-05 中国科学院电工研究所 F-P ultrasonic sensor array implantation device for monitoring partial discharge of power equipment
CN114280436B (en) * 2021-12-24 2024-01-19 中国科学院电工研究所 F-P ultrasonic sensor array implantation device for monitoring partial discharge of power equipment
CN114966547A (en) * 2022-05-18 2022-08-30 珠海视熙科技有限公司 Compensation method, system and device for improving sound source positioning precision
CN117347797A (en) * 2023-09-28 2024-01-05 中国长江电力股份有限公司 Ultrasonic array detection device and method for partial discharge

Also Published As

Publication number Publication date
CN102662132B (en) 2015-09-09

Similar Documents

Publication Publication Date Title
CN102662132A (en) Method and system for locating partial discharge of electrical equipment
CN102866385B (en) Multi-sound-source locating method based on spherical microphone array
CN108318862B (en) Sound source positioning method based on neural network
CN109782231B (en) End-to-end sound source positioning method and system based on multi-task learning
CN101482610A (en) Self-adapting special interference restraint technology for phased array radar
CN102183710A (en) Method and system for identifying numbers of partial discharge sources of transformers
CN101702000A (en) Phased ultrasound locating method and system of partial discharge of transformer
CN103278810B (en) Method for extracting dimension characteristics of underwater target based on space similarity
CN103529364A (en) In-oil multi-local discharge source positioning method based on ultrasonic thinned array sensor
CN106680815A (en) Tensor-sparse-representation-based MIMO radar's imaging method
CN104502904B (en) Torpedo homing beam sharpening method
CN103235286A (en) High-precision locating method for electric noise sources
CN103217211A (en) Substation noise source distribution measuring method based on synthetic aperture principle
CN103616661A (en) Robust far-field narrowband signal source number estimation method
Jiang et al. Raypath separation with high resolution processing
CN103983946A (en) Method for processing singles of multiple measuring channels in sound source localization process
CN110133578B (en) Seabed reflection sound ray incident angle estimation method based on semi-cylindrical volume array
CN103048641A (en) Vector array steady focusing processing method with combined constraint optimization form
CN110208736B (en) Non-circular signal uniform array direction-of-arrival angle estimation method based on fourth-order cumulant
Yan et al. Two‐Step Root‐MUSIC for Direction of Arrival Estimation without EVD/SVD Computation
CN107635181A (en) A kind of multiple access based on channel study perceives the feedback optimized method in source
CN111352075B (en) Underwater multi-sound-source positioning method and system based on deep learning
CN113030849A (en) Near-field source positioning method based on self-encoder and parallel network
CN102183755A (en) Novel high-resolution orientation-estimating method based on Cauchy Gaussian model
CN111190167B (en) Target positioning method of bionic sonar robot

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: STATE POWER ECONOMIC RESEARCH INSTITUTE

Effective date: 20150813

C41 Transfer of patent application or patent right or utility model
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Xie Qing

Inventor after: Liu Dan

Inventor after: Tao Junhan

Inventor after: Li Tong

Inventor after: Liu Xuying

Inventor after: Cheng Shuyi

Inventor after: Formula of law

Inventor before: Xie Qing

Inventor before: Formula of law

Inventor before: Cheng Shuyi

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: XIE QING LV FANGCHENG CHENG SHUYI TO: XIE QING LIU DAN TAO JUNHAN LI TONG LIU XUYING CHENG SHUYI LV FANGCHENG

TA01 Transfer of patent application right

Effective date of registration: 20150813

Address after: 071003 Hebei province Baoding Yonghua No. 619 North Street

Applicant after: NORTH CHINA ELECTRIC POWER University (BAODING)

Applicant after: STATE POWER ECONOMIC Research Institute

Address before: 071003 Hebei province Baoding Yonghua No. 619 North Street

Applicant before: North China Electric Power University (Baoding)

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150909