CN108828401B - A kind of method and device of failure line selection - Google Patents
A kind of method and device of failure line selection Download PDFInfo
- Publication number
- CN108828401B CN108828401B CN201810362236.8A CN201810362236A CN108828401B CN 108828401 B CN108828401 B CN 108828401B CN 201810362236 A CN201810362236 A CN 201810362236A CN 108828401 B CN108828401 B CN 108828401B
- Authority
- CN
- China
- Prior art keywords
- modal components
- phase space
- route
- zero
- sequence current
- 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.)
- Expired - Fee Related
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/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
The embodiment of the present invention provides a kind of method and device of failure line selection, which comprises obtains the zero-sequence current of each route；The zero-sequence current is decomposed respectively, to obtain the corresponding each modal components of each route；Phase space reconfiguration is carried out according to parameter preset to first modal components in each modal components, to obtain phase space reconfiguration vector；The similarity of the phase space reconfiguration vector with identical phase point number two-by-two is sought, and constructs phase space similarity matrix；The integrated correlation coefficient that the corresponding All other routes of each route are sought according to the phase space similarity matrix, using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.Described device executes the above method.Method and device provided in an embodiment of the present invention can be improved the accuracy of faulty line identification.
Description
Technical field
The present embodiments relate to technical field of electric system protection, and in particular to a kind of method and dress of failure line selection
It sets.
Background technique
As specific gravity of the cable run in middle voltage distribution networks is continuously increased, faulty line is accurately identified, power distribution network is pacified
Full stable operation is of great significance.
Currently, in the distribution system of neutral by arc extinction coil grounding mainly using selection method have steady state information method,
Transient information method, external signal method, multi-information fusion method etc..Common steady state information method, such as harmonic component method, zero sequence admittance method
Deng being influenced by factors such as the methods of operation, route selection accuracy rate is low.Change when external signal method includes additional distinctive signal and failure
Two methods of neutral grounding mode are influenced greatly by ground resistance size, position of failure point etc., it may be necessary to which configuration is additionally set
It is standby, it is less economical.Multi-information fusion method integrates multifactor evaluation method and intelligent algorithm mostly, and adaptability is high, but detection process
It is complicated.Traditional transient information method includes first half-wave method, transient energy method etc., and electromagnetism interference is poor, and in certain complexity
Under fault condition, route selection effect is undesirable.In order to solve the deficiency of traditional transient information method, wavelet analysis, S variation, experience
The signal analysis methods such as mode decomposition (EMD) largely use in selection method, greatly improve the accuracy rate of route selection.
Empirical mode decomposition is the process that tranquilization processing is carried out to signal, there is good time-frequency characteristic and adaptivity,
Excellent noise reduction effect, and the bulk information that the high fdrequency component obtained includes original signal is separated, a certain of zero-sequence current can be extracted
Characteristic quantity is as route selection foundation.But EMD is influenced etc. to lack in decomposable process there are modal overlap, end effect, by sample frequency
Point, under strong jamming or the original waveform discomposing effect containing interval signal, pulse signal is undesirable.
Therefore, drawbacks described above how is avoided, improves the accuracy of faulty line identification, becoming need solve the problems, such as.
Summary of the invention
In view of the problems of the existing technology, the embodiment of the present invention provides a kind of method and device of failure line selection.
In a first aspect, the embodiment of the present invention provides a kind of method of failure line selection, which comprises
Obtain the zero-sequence current of each route；
The zero-sequence current is decomposed respectively, to obtain the corresponding each modal components of each route；
Phase space reconfiguration is carried out according to parameter preset to first modal components in each modal components, to obtain phase
Space Reconstruction vector；
The similarity of the phase space reconfiguration vector with identical phase point number two-by-two is sought, and constructs phase space similarity moment
Battle array；
The integrated correlation coefficient of the corresponding All other routes of each route is sought according to the phase space similarity matrix,
Using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.
Second aspect, the embodiment of the present invention provide a kind of device of failure line selection, and described device includes:
Acquiring unit, for obtaining the zero-sequence current of each route；
Decomposition unit, for being decomposed respectively to the zero-sequence current, to obtain the corresponding each mould of each route
State component；
Reconfiguration unit, for carrying out phase space according to parameter preset to first modal components in each modal components
Reconstruct, to obtain phase space reconfiguration vector；
Construction unit for seeking the similarity of the phase space reconfiguration vector with identical phase point number two-by-two, and constructs
Phase space similarity matrix；
Unit is sought, for seeking the corresponding All other routes of each route according to the phase space similarity matrix
Integrated correlation coefficient, using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.
The third aspect, the embodiment of the present invention provide a kind of electronic equipment, comprising: processor, memory and bus, wherein
The processor and the memory complete mutual communication by the bus；
The memory is stored with the program instruction that can be executed by the processor, and the processor calls described program to refer to
Order is able to carry out following method:
Obtain the zero-sequence current of each route；
The zero-sequence current is decomposed respectively, to obtain the corresponding each modal components of each route；
Phase space reconfiguration is carried out according to parameter preset to first modal components in each modal components, to obtain phase
Space Reconstruction vector；
The similarity of the phase space reconfiguration vector with identical phase point number two-by-two is sought, and constructs phase space similarity moment
Battle array；
The integrated correlation coefficient of the corresponding All other routes of each route is sought according to the phase space similarity matrix,
Using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.
Fourth aspect, the embodiment of the present invention provide a kind of non-transient computer readable storage medium, comprising:
The non-transient computer readable storage medium stores computer instruction, and the computer instruction makes the computer
Execute following method:
Obtain the zero-sequence current of each route；
The zero-sequence current is decomposed respectively, to obtain the corresponding each modal components of each route；
Phase space reconfiguration is carried out according to parameter preset to first modal components in each modal components, to obtain phase
Space Reconstruction vector；
The similarity of the phase space reconfiguration vector with identical phase point number two-by-two is sought, and constructs phase space similarity moment
Battle array；
The integrated correlation coefficient of the corresponding All other routes of each route is sought according to the phase space similarity matrix,
Using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.
The method and device of failure line selection provided in an embodiment of the present invention carries out phase space by the modal components to route
Reconstruct, and integrated correlation coefficient is obtained by constructing, seeking phase space similarity matrix, by the smallest integrated correlation coefficient of numerical value
Corresponding route can be improved the accuracy of faulty line identification as faulty line.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is this hair
Bright some embodiments for those of ordinary skill in the art without creative efforts, can be with root
Other attached drawings are obtained according to these attached drawings.
Fig. 1 is the method flow schematic diagram of failure line selection of the embodiment of the present invention；
Fig. 2 is the specific flow chart that VMD of the embodiment of the present invention is decomposed；
Fig. 3 is the specific flow chart of phase space reconfiguration of the embodiment of the present invention；
Fig. 4 is the detailed process figure that integrated correlation coefficient of the embodiment of the present invention calculates；
Fig. 5 is the apparatus structure schematic diagram of failure line selection of the embodiment of the present invention；
Fig. 6 is electronic equipment entity structure schematic diagram provided in an embodiment of the present invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
Fig. 1 is the method flow schematic diagram of failure line selection of the embodiment of the present invention, as shown in Figure 1, the embodiment of the present invention provides
A kind of failure line selection method, comprising the following steps:
S101: the zero-sequence current of each route is obtained.
Specifically, device obtains the zero-sequence current of each route.Each route in the embodiment of the present invention, which can be, to be based on
What the resonator system of variation mode decomposition (VMD) was chosen, VMD is substantially one group of adaptive wiener filter, the onrecurrent of use
Method will not transmission error, have good noise robustness, the advantage of EMD can not only be retained, and can make up well EMD's
It is insufficient.The zero-sequence current of each route can be obtained by way of half waveform after extracting failure, which is
Transient zero-sequence current.Uniform sampling frequency intercepts failure the latter cycle data, for the zero-sequence current of any one route
(assuming that having N route), it is assumed that the data sampled are fi(t)。
S102: respectively decomposing the zero-sequence current, to obtain the corresponding each modal components of each route.
Specifically, device respectively decomposes the zero-sequence current, to obtain the corresponding each mode of each route
Component.VMD decomposition, f are carried out to each sampled datai(t) it decomposes to obtain each modal components by VMD:
{u1,u2,...uk}
Fig. 2 is the specific flow chart that VMD of the embodiment of the present invention is decomposed, as shown in Fig. 2, specifically can be in the following way:
First to fi(t) the variational problem model of belt restraining is established, these modal components have the limited band of centre frequency
Width, and estimate the sum of bandwidth minimum:
In formula: { ukIt is each modal components；{ωkIt is each modal components centre frequency；δ (t) is impulse function；It is to t
Partial derivative is sought, argument Lagrange function is introduced, converts non-binding variational problem for restrictive variational problem:
In formula: α is secondary penalty factor, and λ is Lagrangian.It is sought using alternating direction Multiplier Algorithm (ADMM)
The optimal solution of constraint variation model can be obtained in the saddle point for stating Lagrangian.
VMD is decomposed it needs to be determined that modal components number (optimal each modal components number) and secondary penalty factor two ginsengs
Number.Secondary penalty factor can directly choose the numerical values recited of above-mentioned sample frequency.The determination of optimal each modal components number
Step can be such that using db10 wavelet packet generating function, and 5 layers of wavelet decomposition of zero sequence current signal progress are obtained 32 son frequencies
Band；Calculate each frequency band energy and energy ratio；Selected modal components number allows upper limit kmax, it is more than 1/k with frequency band energy ratiomax
Sub-band number primarily determine modal components number；If there is mode centre frequency and finally overlaps phenomenon in VMD decomposition,
Modal components number is gradually reduced, optimal modal components number is finally obtained.
The Wave data of zero-sequence current involved in the method for the embodiment of the present invention answers sample frequency having the same, if
Sample frequency is different, needs to carry out multi-sampling rate signal processing to waveform: reducing sample frequency by extracting data；Pass through interpolation
Improve sample frequency；Any fraction time sampling frequency conversion is realized with the cascade method of interpolation by extracting.
S103: carrying out phase space reconfiguration according to parameter preset to first modal components in each modal components, with
Obtain phase space reconfiguration vector.
Specifically, device carries out phase space weight according to parameter preset to first modal components in each modal components
Structure, to obtain phase space reconfiguration vector.The parameter preset includes d dimension, time delays τ.Fig. 3 is phase space of the embodiment of the present invention
The specific flow chart of reconstruct, as shown in figure 3, phase space reconfiguration can be in the following way:
Take first modal components u1={ x1,x2,...,xnCarry out d dimension, the phase space reconfiguration that time delays are τ, it obtains
Phase space indicate are as follows:
X={ X1,X2,...,Xm}
Wherein, Xi={ xi,xi+τ,...,xi+(d-1)τIt is each phase point；
M=n- (d-1) τ is mutually to count；
Dimension and time delays parameter can be acquired by C-C algorithm, seek the average value that all n dimensions are calculated and when
Between the average value that is delayed respectively as the specific value of dimension and time delays parameter, (the unified dimension reconstructed in corresponding diagram 3 is joined
Several and time delays parameter).
S104: the similarity of the phase space reconfiguration vector with identical phase point number two-by-two is sought, and constructs Phase Space
Like degree matrix.
Specifically, device seeks the similarity of the phase space reconfiguration vector with identical phase point number two-by-two, and construct phase
Space similarity matrix.Fig. 4 is the detailed process figure that integrated correlation coefficient of the embodiment of the present invention calculates, as shown in figure 4, firstly,
Phase space similarity is calculated using normalized crosscorrelation method, for two sample phase space X=with identical phase point number
{X1,X2,...,XmAnd Y={ Y1,Y2,...,Ym, then phase space similarity are as follows:
For N route, the phase space similarity of route zero-sequence current two-by-two is calculated, constructs phase space similarity matrix,
Construct phase space similarity matrix are as follows:
M=(Pij)NN
The phase space similarity matrix is a symmetrical matrix, and diagonal element is all 1, the element marked as (i, j)
Indicate the similarity of the phase space two-by-two of i-th line road and j-th strip route.
S105: the comprehensive related of the corresponding All other routes of each route is sought according to the phase space similarity matrix
Coefficient, using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.
Specifically, device seeks the comprehensive of the corresponding All other routes of each route according to the phase space similarity matrix
Related coefficient is closed, using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.
Integrated correlation coefficient of the every route relative to All other routes is sought according to phase space similarity matrix:
The integrated correlation coefficient is to take every data line, and rower i therein indicates i-th line road, removal 1 and minimum
Value, takes integrated correlation coefficient of the average value of remaining N-2 number as i-th line road.
Since the size of the transient zero-sequence current of sound circuit is only related with each line mutual-ground capacitor size, change frequency base
This is identical, and the variation of faulty line frequency relatively lags behind, and zero-sequence current amplitude is also big compared with sound circuit, it may be assumed that faulty line zero sequence
Period, amplitude and the phase of electric current are larger with sound circuit difference.So the smallest route of integrated correlation coefficient is failure
Route, all elements of phase space similarity matrix are both greater than a setting valve, and maximum integrated correlation coefficient and minimum are comprehensive
Phase relation difference is less than setting valve and is then determined as bus-bar fault.
The method of the embodiment of the present invention has the advantage that as follows:
1, for different circuit types (including overhead transmission line, cable run, aerial cable joint line), different faults item
Single-phase fault under part (positional distance bus distance occurs for failure, initial phase angle occurs for failure, ground resistance etc. occurs for failure), all
It can correct route selection.
2, being decomposed using VMD has very strong noise immunity.
3, has the advantages that adaptivity route selection.
The method of failure line selection provided in an embodiment of the present invention carries out phase space reconfiguration by the modal components to route,
And integrated correlation coefficient is obtained by constructing, seeking phase space similarity matrix, the smallest integrated correlation coefficient of numerical value is corresponding
Route as faulty line, can be improved the accuracy of faulty line identification.
On the basis of the above embodiments, described that the zero-sequence current is decomposed respectively, to obtain each route point
Not corresponding each modal components, comprising:
Belt restraining Variation Model is established to the zero-sequence current, each modal components have the limited band of centre frequency
Width, and estimate that the sum of bandwidth minimum, the belt restraining Variation Model are as follows:
Wherein, { ukIt is each modal components；{ωkIt is each modal components centre frequency；δ (t) is impulse function；It is to t
Seek partial derivative；F (t) is the zero-sequence current.
Specifically, device establishes belt restraining Variation Model to the zero-sequence current, each modal components have center frequency
The finite bandwidth of rate, and estimate that the sum of bandwidth minimum, the belt restraining Variation Model are as follows:
Wherein, { ukIt is each modal components；{ωkIt is each modal components centre frequency；δ (t) is impulse function；It is to t
Seek partial derivative；F (t) is the zero-sequence current.It can refer to above-described embodiment, repeat no more.
Following argument Lagrange function is introduced, converts non-binding variational problem for restrictive variational problem:
Wherein, α is secondary penalty factor, λ is Lagrangian；The numbers of optimal each modal components and the α according to
Preset rules, which calculate, to be obtained.
Specifically, device introduces following argument Lagrange function, non-binding change is converted by restrictive variational problem
Divide problem:
Wherein, α is secondary penalty factor, λ is Lagrangian；The numbers of optimal each modal components and the α according to
Preset rules, which calculate, to be obtained.It can refer to above-described embodiment, repeat no more.
The saddle point of the argument Lagrange function is sought using alternating direction Multiplier Algorithm, it is described wait constrain change to obtain
The optimal solution of sub-model, and using the optimal solution as each modal components.
Specifically, device seeks the saddle point of the argument Lagrange function using alternating direction Multiplier Algorithm, to obtain
The optimal solution to constraint variation model, and using the optimal solution as each modal components.It can refer to above-described embodiment,
It repeats no more.
The method of failure line selection provided in an embodiment of the present invention, by solving the belt restraining Variation Model established and converted,
Using solving result as each modal components, each modal components can be quickly got.
On the basis of the above embodiments, the zero-sequence current is collected according to identical sample frequency；Correspondingly,
It is described to include: according to the number and the α of preset rules calculating acquisition optimal each modal components
Using the numerical value of the identical sample frequency as the numerical value of the α.
Specifically, device is using the numerical value of the identical sample frequency as the numerical value of the α.It can refer to above-mentioned implementation
Example, repeats no more.
The determination of the number of optimal each modal components, comprising:
Using db10 wavelet packet generating function, zero sequence current signal is subjected to 5 layers of wavelet decomposition, obtains 32 sub-bands.
Specifically, device uses db10 wavelet packet generating function, zero sequence current signal is subjected to 5 layers of wavelet decomposition, obtains 32
A sub-band.It can refer to above-described embodiment, repeat no more.
Each sub-band energy is calculated, and sub-band energy ratio is calculated according to the sub-band energy.
Specifically, device calculates each sub-band energy, and sub-band energy ratio is calculated according to the sub-band energy.It can
Referring to above-described embodiment, repeat no more.
Selected modal components number allows upper limit kmax, it is more than 1/k with the sub-band energy ratiomaxSub-band number at the beginning of
Step determines the number of each modal components.
Specifically, device, which selectes modal components number, allows upper limit kmax, it is more than 1/k with the sub-band energy ratiomax's
Sub-band number primarily determines the number of each modal components.It can refer to above-described embodiment, repeat no more.
If decomposing the mode centre frequency occurred to overlap, the number of each modal components is gradually reduced, so that
Mode centre frequency does not overlap, and the number of each modal components at this time is determined as to the number of optimal modal component.
If being gradually reduced described each specifically, device judgement knows that the mode centre frequency decomposed and occurred overlaps
The number of each modal components at this time is determined as optimal by the number of modal components so that mode centre frequency does not overlap
The number of modal components.It can refer to above-described embodiment, repeat no more.
The method of failure line selection provided in an embodiment of the present invention, by reasonably determine optimal each modal components number and
α can guarantee going on smoothly for this method.
On the basis of the above embodiments, the method also includes:
If the sample frequency is different, it is adjusted in the following way:
The sample frequency is reduced by extracting data；Or, improving the sample frequency by interpolated data.
If being adjusted in the following way specifically, device judgement knows that the sample frequency is different:
The sample frequency is reduced by extracting data；Or, improving the sample frequency by interpolated data.It can refer to
Embodiment is stated, is repeated no more.
The method of failure line selection provided in an embodiment of the present invention ensure that zero-sequence current is root by making corresponding adjustment
It is collected according to identical sample frequency, it ensure that being normally carried out for this method.
On the basis of the above embodiments, the parameter preset includes d dimension, time delays τ；Correspondingly, described to described
First modal components in each modal components carry out phase space reconfiguration according to parameter preset, to obtain phase space reconfiguration vector,
Include:
Take first modal components u1={ x1,x2,...,xnCarry out d dimension, time delays τ phase space reconfiguration, obtain
Phase space reconfiguration vector are as follows:
X={ X1,X2,...,Xm}
Wherein, Xi={ xi,xi+τ,…,xi+(d-1)τ } it is each phase point；M=n- (d-1) τ is mutually to count.
Specifically, device takes first modal components u1={ x1,x2,…,xnCarry out d dimension, the phase space of time delays τ
Reconstruct, obtained phase space reconfiguration vector are as follows:
X={ X1,X2,...,Xm}
Wherein, Xi={ xi,xi+τ,...,xi+(d-1)τ } it is each phase point；M=n- (d-1) τ is mutually to count.It can refer to above-mentioned
Embodiment repeats no more.
The method of failure line selection provided in an embodiment of the present invention can effectively and reasonably carry out each modal components mutually empty
Between reconstruct, be further able to improve faulty line identification accuracy.
On the basis of the above embodiments, the determination of the d dimension and the specific value of the time delays τ, comprising:
Using the average value of the average value for all n dimensions being calculated using C-C algorithm and time delays as d
The specific value of the specific value of dimension and the time delays τ.
Specifically, device is by the average value of all n dimensions being calculated using C-C algorithm and being averaged for time delays
It is worth the specific value respectively as the d specific value tieed up and the time delays τ.It can refer to above-described embodiment, repeat no more.
The method of failure line selection provided in an embodiment of the present invention, by rationally determining the specific value of parameter preset, into one
Step can be improved the accuracy of faulty line identification.
On the basis of the above embodiments, described each route is sought according to the phase space similarity matrix to respectively correspond
All other routes integrated correlation coefficient, comprising:
The integrated correlation coefficient is calculated according to the following formula:
Wherein, ρiFor sum that the integrated correlation coefficient, N are route, PijFor in the phase space similarity matrix
Element.
Specifically, device calculates the integrated correlation coefficient according to the following formula:
Wherein, ρiFor sum that the integrated correlation coefficient, N are route, PijFor in the phase space similarity matrix
Element.It can refer to above-described embodiment, repeat no more.
The method of failure line selection provided in an embodiment of the present invention calculates integrated correlation coefficient by specific formula, effectively obtains
Integrated correlation coefficient has been taken, has been further able to improve the accuracy of faulty line identification.
Fig. 5 is the apparatus structure schematic diagram of failure line selection of the embodiment of the present invention, as shown in figure 5, the embodiment of the present invention provides
A kind of device of failure line selection, including acquiring unit 501, decomposition unit 502, reconfiguration unit 503, construction unit 504 and ask
Take unit 505, in which:
Acquiring unit 501 is used to obtain the zero-sequence current of each route；Decomposition unit 502 is for respectively to the zero sequence electricity
Stream is decomposed, to obtain the corresponding each modal components of each route；Reconfiguration unit 503 is used for each modal components
In first modal components according to parameter preset carry out phase space reconfiguration, to obtain phase space reconfiguration vector；Construction unit 504
For seeking the similarity of the phase space reconfiguration vector with identical phase point number two-by-two, and construct phase space similarity matrix；
Unit 505 is sought for seeking the synthesis phase of the corresponding All other routes of each route according to the phase space similarity matrix
Relationship number, using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.
Specifically, acquiring unit 501 is used to obtain the zero-sequence current of each route；Decomposition unit 502 is used for respectively to institute
It states zero-sequence current to be decomposed, to obtain the corresponding each modal components of each route；Reconfiguration unit 503 is used for described each
First modal components in modal components carry out phase space reconfiguration according to parameter preset, to obtain phase space reconfiguration vector；Structure
Unit 504 is built for seeking the similarity of the phase space reconfiguration vector with identical phase point number two-by-two, and constructs Phase Space
Like degree matrix；Unit 505 is sought for seeking other corresponding lines of each route according to the phase space similarity matrix
The integrated correlation coefficient on road, using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.
The device of failure line selection provided in an embodiment of the present invention carries out phase space reconfiguration by the modal components to route,
And integrated correlation coefficient is obtained by constructing, seeking phase space similarity matrix, the smallest integrated correlation coefficient of numerical value is corresponding
Route as faulty line, can be improved the accuracy of faulty line identification.
The device of failure line selection provided in an embodiment of the present invention specifically can be used for executing the place of above-mentioned each method embodiment
Process is managed, details are not described herein for function, is referred to the detailed description of above method embodiment.
Fig. 6 is electronic equipment entity structure schematic diagram provided in an embodiment of the present invention, as shown in fig. 6, the electronic equipment
It include: processor (processor) 601, memory (memory) 602 and bus 603；
Wherein, the processor 601, memory 602 complete mutual communication by bus 603；
The processor 601 is used to call the program instruction in the memory 602, to execute above-mentioned each method embodiment
Provided method, for example, obtain the zero-sequence current of each route；The zero-sequence current is decomposed respectively, to obtain
Take the corresponding each modal components of each route；To first modal components in each modal components according to parameter preset
Phase space reconfiguration is carried out, to obtain phase space reconfiguration vector；Seek phase space reconfiguration vector two-by-two with identical phase point number
Similarity, and construct phase space similarity matrix；Each route is sought according to the phase space similarity matrix to respectively correspond
All other routes integrated correlation coefficient, using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.
The present embodiment discloses a kind of computer program product, and the computer program product includes being stored in non-transient calculating
Computer program on machine readable storage medium storing program for executing, the computer program include program instruction, when described program instruction is calculated
When machine executes, computer is able to carry out method provided by above-mentioned each method embodiment, for example, obtains the zero of each route
Sequence electric current；The zero-sequence current is decomposed respectively, to obtain the corresponding each modal components of each route；To described each
First modal components in modal components carry out phase space reconfiguration according to parameter preset, to obtain phase space reconfiguration vector；It asks
The similarity of the phase space reconfiguration vector with identical phase point number two-by-two is taken, and constructs phase space similarity matrix；According to institute
The integrated correlation coefficient that phase space similarity matrix seeks the corresponding All other routes of each route is stated, numerical value is the smallest comprehensive
The corresponding route of related coefficient is closed as faulty line.
The present embodiment provides a kind of non-transient computer readable storage medium, the non-transient computer readable storage medium
Computer instruction is stored, the computer instruction makes the computer execute method provided by above-mentioned each method embodiment, example
It such as include: the zero-sequence current for obtaining each route；The zero-sequence current is decomposed respectively, it is right respectively to obtain each route
Each modal components answered；Phase space reconfiguration is carried out according to parameter preset to first modal components in each modal components,
To obtain phase space reconfiguration vector；The similarity of the phase space reconfiguration vector with identical phase point number two-by-two is sought, and is constructed
Phase space similarity matrix；The synthesis of the corresponding All other routes of each route is sought according to the phase space similarity matrix
Related coefficient, using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line.
Those of ordinary skill in the art will appreciate that: realize that all or part of the steps of above method embodiment can pass through
The relevant hardware of program instruction is completed, and program above-mentioned can be stored in a computer readable storage medium, the program
When being executed, step including the steps of the foregoing method embodiments is executed；And storage medium above-mentioned includes: ROM, RAM, magnetic disk or light
The various media that can store program code such as disk.
The embodiments such as electronic equipment described above are only schematical, wherein it is described as illustrated by the separation member
Unit may or may not be physically separated, and component shown as a unit may or may not be object
Manage unit, it can it is in one place, or may be distributed over multiple network units.It can select according to the actual needs
Some or all of the modules therein is selected to achieve the purpose of the solution of this embodiment.Those of ordinary skill in the art are not paying wound
In the case where the labour for the property made, it can understand and implement.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
It realizes by means of software and necessary general hardware platform, naturally it is also possible to pass through hardware.Based on this understanding, on
Stating technical solution, substantially the part that contributes to existing technology can be embodied in the form of software products in other words, should
Computer software product may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, CD, including several fingers
It enables and using so that a computer equipment (can be personal computer, server or the network equipment etc.) executes each implementation
Method described in certain parts of example or embodiment.
Finally, it should be noted that the above various embodiments is only to illustrate the technical solution of the embodiment of the present invention, rather than it is right
It is limited；Although the embodiment of the present invention is described in detail referring to foregoing embodiments, the ordinary skill of this field
Personnel are it is understood that it is still possible to modify the technical solutions described in the foregoing embodiments, or to part
Or all technical features are equivalently replaced；And these are modified or replaceed, it does not separate the essence of the corresponding technical solution
The range of each embodiment technical solution of the embodiment of the present invention.
Claims (8)
1. a kind of method of failure line selection characterized by comprising
Obtain the zero-sequence current of each route；
The zero-sequence current is decomposed respectively, to obtain the corresponding each modal components of each route；
Phase space reconfiguration is carried out according to parameter preset to first modal components in each modal components, to obtain phase space
Reconstruct vector；
The similarity of the phase space reconfiguration vector with identical phase point number two-by-two is sought, and constructs phase space similarity matrix；
The integrated correlation coefficient that the corresponding All other routes of each route are sought according to the phase space similarity matrix, will count
It is worth the corresponding route of the smallest integrated correlation coefficient as faulty line；
It is described that the zero-sequence current is decomposed respectively, to obtain the corresponding each modal components of each route, comprising:
Belt restraining Variation Model is established to the zero-sequence current, each modal components have the finite bandwidth of centre frequency, and
Estimate that the sum of bandwidth minimum, the belt restraining Variation Model are as follows:
Wherein, { ukIt is each modal components；{ωkIt is each modal components centre frequency；δ (t) is impulse function；It is that inclined is asked to t
Derivative；F (t) is the zero-sequence current；
Following argument Lagrange function is introduced, converts non-binding variational problem for restrictive variational problem:
Wherein, α is secondary penalty factor, λ is Lagrangian；The number of optimal each modal components and the α are according to default
Rule, which calculates, to be obtained；
The saddle point of the argument Lagrange function is sought using alternating direction Multiplier Algorithm, it is described to constraint variation mould to obtain
The optimal solution of type, and using the optimal solution as each modal components；
The zero-sequence current is collected according to identical sample frequency；Correspondingly, described calculated according to preset rules obtains
The number of optimal each modal components and the α include:
Using the numerical value of the identical sample frequency as the numerical value of the α；
The determination of the number of optimal each modal components, comprising:
Using db10 wavelet packet generating function, zero sequence current signal is subjected to 5 layers of wavelet decomposition, obtains 32 sub-bands；
Each sub-band energy is calculated, and sub-band energy ratio is calculated according to the sub-band energy；
Selected modal components number allows upper limit kmax, it is more than 1/k with the sub-band energy ratiomaxSub-band number it is tentatively true
The number of fixed each modal components；
If decomposing the mode centre frequency occurred to overlap, the number of each modal components is gradually reduced, so that mode
Centre frequency does not overlap, and the number of each modal components at this time is determined as to the number of optimal modal component.
2. the method according to claim 1, wherein the method also includes:
If the sample frequency is different, it is adjusted in the following way:
The sample frequency is reduced by extracting data；
Or,
The sample frequency is improved by interpolated data.
3. method according to claim 1 or 2, which is characterized in that the parameter preset includes d dimension, time delays τ；Phase
It answers, first modal components in each modal components carry out phase space reconfiguration according to parameter preset, to obtain
Phase space reconfiguration vector, comprising:
Take first modal components u1={ x1,x2,...,xnThe phase space reconfiguration of d dimension, time delays τ is carried out, obtained phase is empty
Between reconstruct vector are as follows:
X={ X1,X2,...,Xm}
Wherein, Xi={ xi,xi+τ,...,xi+(d-1)τIt is each phase point；M=n- (d-1) τ is mutually to count.
4. according to the method described in claim 3, it is characterized in that, the specific value of d dimension and the time delays τ
It determines, comprising:
The average value of all n dimensions and the average value of time delays that are calculated using C-C algorithm are tieed up as d
The specific value of specific value and the time delays τ.
5. method according to claim 1 or 2, which is characterized in that described to be sought according to the phase space similarity matrix
The integrated correlation coefficient of the corresponding All other routes of each route, comprising:
The integrated correlation coefficient is calculated according to the following formula:
Wherein, ρiFor sum that the integrated correlation coefficient, N are route, PijFor the element in the phase space similarity matrix.
6. a kind of device of failure line selection characterized by comprising
Acquiring unit, for obtaining the zero-sequence current of each route；
Decomposition unit, for being decomposed respectively to the zero-sequence current, to obtain the corresponding each mode point of each route
Amount；
Reconfiguration unit, for carrying out phase space weight according to parameter preset to first modal components in each modal components
Structure, to obtain phase space reconfiguration vector；
Construction unit for seeking the similarity of the phase space reconfiguration vector with identical phase point number two-by-two, and constructs mutually empty
Between similarity matrix；
Unit is sought, for seeking the synthesis of the corresponding All other routes of each route according to the phase space similarity matrix
Related coefficient, using the corresponding route of the smallest integrated correlation coefficient of numerical value as faulty line；
The decomposition unit is specifically used for:
Belt restraining Variation Model is established to the zero-sequence current, each modal components have the finite bandwidth of centre frequency, and
Estimate that the sum of bandwidth minimum, the belt restraining Variation Model are as follows:
Wherein, { ukIt is each modal components；{ωkIt is each modal components centre frequency；δ (t) is impulse function；It is that inclined is asked to t
Derivative；F (t) is the zero-sequence current；
Following argument Lagrange function is introduced, converts non-binding variational problem for restrictive variational problem:
Wherein, α is secondary penalty factor, λ is Lagrangian；The number of optimal each modal components and the α are according to default
Rule, which calculates, to be obtained；
The saddle point of the argument Lagrange function is sought using alternating direction Multiplier Algorithm, it is described to constraint variation mould to obtain
The optimal solution of type, and using the optimal solution as each modal components；
The zero-sequence current is collected according to identical sample frequency；Correspondingly, the decomposition unit is also used to:
Using the numerical value of the identical sample frequency as the numerical value of the α；
The determination of the number of optimal each modal components, comprising:
Using db10 wavelet packet generating function, zero sequence current signal is subjected to 5 layers of wavelet decomposition, obtains 32 sub-bands；
Each sub-band energy is calculated, and sub-band energy ratio is calculated according to the sub-band energy；
Selected modal components number allows upper limit kmax, it is more than 1/k with the sub-band energy ratiomaxSub-band number it is tentatively true
The number of fixed each modal components；
If decomposing the mode centre frequency occurred to overlap, the number of each modal components is gradually reduced, so that mode
Centre frequency does not overlap, and the number of each modal components at this time is determined as to the number of optimal modal component.
7. a kind of electronic equipment characterized by comprising processor, memory and bus, wherein
The processor and the memory complete mutual communication by the bus；
The memory is stored with the program instruction that can be executed by the processor, and the processor calls described program to instruct energy
Enough execute method as claimed in claim 1 to 5.
8. a kind of non-transient computer readable storage medium, which is characterized in that the non-transient computer readable storage medium is deposited
Computer instruction is stored up, the computer instruction makes the computer execute method as claimed in claim 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810362236.8A CN108828401B (en) | 2018-04-20 | 2018-04-20 | A kind of method and device of failure line selection |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810362236.8A CN108828401B (en) | 2018-04-20 | 2018-04-20 | A kind of method and device of failure line selection |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108828401A CN108828401A (en) | 2018-11-16 |
| CN108828401B true CN108828401B (en) | 2019-11-15 |
Family
ID=64154521
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810362236.8A Expired - Fee Related CN108828401B (en) | 2018-04-20 | 2018-04-20 | A kind of method and device of failure line selection |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108828401B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110261723B (en) * | 2019-06-05 | 2021-05-18 | 河南理工大学 | Low-current grounding line selection method based on coefficient of variation and high-order cumulant |
| CN110161125A (en) * | 2019-06-17 | 2019-08-23 | 哈尔滨工业大学 | The Aeroengine Smart monitoring method combined based on acceleration with sound emission cognition technology |
| CN111157843B (en) * | 2020-01-06 | 2022-04-12 | 长沙理工大学 | Power distribution network line selection method based on time-frequency domain traveling wave information |
| CN111413588B (en) * | 2020-03-31 | 2022-03-11 | 陕西省地方电力（集团）有限公司咸阳供电分公司 | Power distribution network single-phase earth fault line selection method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104849050A (en) * | 2015-06-02 | 2015-08-19 | 安徽工业大学 | Rolling bearing fault diagnosis method based on composite multi-scale permutation entropies |
| CN104865499A (en) * | 2015-05-11 | 2015-08-26 | 昆明理工大学 | Super-high voltage direct-current power transmission line region internal and external fault identification method |
| CN106375157A (en) * | 2016-10-31 | 2017-02-01 | 华侨大学 | Phase-space-reconstruction-based network flow correlation method |
| CN107491832A (en) * | 2017-07-12 | 2017-12-19 | 国网上海市电力公司 | Energy quality steady-state index prediction method based on chaology |
-
2018
- 2018-04-20 CN CN201810362236.8A patent/CN108828401B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104865499A (en) * | 2015-05-11 | 2015-08-26 | 昆明理工大学 | Super-high voltage direct-current power transmission line region internal and external fault identification method |
| CN104849050A (en) * | 2015-06-02 | 2015-08-19 | 安徽工业大学 | Rolling bearing fault diagnosis method based on composite multi-scale permutation entropies |
| CN106375157A (en) * | 2016-10-31 | 2017-02-01 | 华侨大学 | Phase-space-reconstruction-based network flow correlation method |
| CN107491832A (en) * | 2017-07-12 | 2017-12-19 | 国网上海市电力公司 | Energy quality steady-state index prediction method based on chaology |
Non-Patent Citations (3)
| Title |
|---|
| 基于EEMD和关联维数的小电流接地故障选线方法;彭泓等;《测控技术》;20170130;第36卷(第1期);124-128 * |
| 基于改进变分模态分解和Hilbert变换的变压器局部放电信号特征提取及分类;朱永利等;《电工技术学报》;20170925;第32卷(第9期);221-235 * |
| 谐振接地电网故障选线相关分析法;束洪春等;《电力自动化设备》;20080930;第28卷(第9期);6-9,48 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108828401A (en) | 2018-11-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108828401B (en) | A kind of method and device of failure line selection | |
| CN103344875B (en) | Classification line selection method for single-phase earth fault of resonance earthing system | |
| CN102854437B (en) | Fault line selection method of low current grounding system using time-frequency atom decomposition theory | |
| Patel | A new FDOST entropy based intelligent digital relaying for detection, classification and localization of faults on the hybrid transmission line | |
| Aslan et al. | Artificial neural-network-based fault location for power distribution lines using the frequency spectra of fault data | |
| US10312041B2 (en) | Frequency measurement for electric power delivery system | |
| CN112051480A (en) | Neural network power distribution network fault diagnosis method and system based on variational modal decomposition | |
| CN109100568A (en) | Power grid harmonic wave general measurement analysis method, device, equipment and medium | |
| CN115480203A (en) | Online quantitative evaluation method and system for error state of current transformer | |
| Fan et al. | Transmission line fault location using deep learning techniques | |
| CN110967594A (en) | Method and device for positioning faults of power distribution network containing inverter type distributed power supply | |
| Arghandeh et al. | Phasor-based assessment for harmonic sources in distribution networks | |
| Mishra et al. | Intrinsic time decomposition based fault location scheme for unified power flow controller compensated transmission line | |
| CN102570490A (en) | Intelligent charge-discharge control method for electric vehicle | |
| CN110687390A (en) | Low-current fusion line selection method based on extension theory | |
| Morsi et al. | On the application of wavelet transform for symmetrical components computations in the presence of stationary and non-stationary power quality disturbances | |
| CN106372800A (en) | Method and device for acquiring subway power supply quality | |
| Jose et al. | Smart grid data compression of power quality events using wavelet transform | |
| Patel | Superimposed components of Lissajous pattern based feature extraction for classification and localization of transmission line faults | |
| CN106548252A (en) | Based on the distribution network structure constructing system for improving K means clusters | |
| Paolone et al. | An automatic system to locate phase-to-ground faults in medium voltage cable networks based on the wavelet analysis of high-frequency signals | |
| CN109116142A (en) | The anti-interference emulation test method of track circuit receiver and system | |
| CN114062832A (en) | Method and system for identifying short-circuit fault type of power distribution network | |
| CN115144696A (en) | Fault line selection method, device, equipment and medium for low-current grounding system | |
| Sudha et al. | A comparison between different approaches for fault classification in transmission lines |
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 | ||
| 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: 20191115 |