CN105759169B - A kind of method for diagnosing faults based on WAMS times section information and topology information - Google Patents
A kind of method for diagnosing faults based on WAMS times section information and topology information Download PDFInfo
- Publication number
- CN105759169B CN105759169B CN201610127868.7A CN201610127868A CN105759169B CN 105759169 B CN105759169 B CN 105759169B CN 201610127868 A CN201610127868 A CN 201610127868A CN 105759169 B CN105759169 B CN 105759169B
- Authority
- CN
- China
- Prior art keywords
- phase
- fault
- measuring point
- mode vector
- amplitude
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 239000013598 vector Substances 0.000 claims abstract description 98
- 238000003745 diagnosis Methods 0.000 claims abstract description 65
- 230000002123 temporal effect Effects 0.000 claims abstract description 24
- 230000000977 initiatory effect Effects 0.000 claims abstract description 6
- 238000005259 measurement Methods 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 abstract description 7
- 230000005540 biological transmission Effects 0.000 description 16
- 230000002159 abnormal effect Effects 0.000 description 6
- 230000005611 electricity Effects 0.000 description 5
- 238000003909 pattern recognition Methods 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000013024 troubleshooting Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
- Testing Electric Properties And Detecting Electric Faults (AREA)
Abstract
The invention discloses a kind of method for diagnosing faults based on WAMS times section information and topology information.This method includes:(1) three-phase voltage of each measuring point and the amplitude of electric current and phase-angle data in interconnected network are obtained by WAMS;(2) it is calculated using the amplitude data of the three-phase voltage of acquisition and electric current according to fault diagnosis start-up criterion;(3) real-time characteristic amount extraction is carried out to the amplitude data of measuring point three-phase voltage and electric current for meeting fault initiating criterion, forms real-time mode vector;(4) real-time mode vector and the reference mode vector set of preset known fault type that measuring point i is extracted are matched, preliminary fail result is provided according to matching result;(5) fault diagnosis result for meeting the measuring point of start-up criterion is integrated according to temporal information and power network topology information, obtains final fault diagnosis result.For the present invention by WAMS data for electric network failure diagnosis, diagnosis speed is fast, as a result accurately, can improve power system security stable operation level.
Description
Technical field
The invention belongs to electric network failure diagnosis technical fields, and WAMS (wide area measurement systems are based on more particularly, to one kind
System, Wide-Area Measurement System) time section information and topology information method for diagnosing faults.
Background technology
With the high speed development of Chinese national economy, power grid scale constantly expands, electric network composition even more complex, current electricity
Force system is faced with increasing challenge.Electric network fault caused by exceedingly odious weather and equipment fault is steady to the safety of power grid
Fixed operation causes greatly to threaten.When an error occurs, if dispatching of power netwoks operations staff can obtain fault type, event in time
The information such as Downtime contribute to it to make correct judgement accordingly, take appropriate processing mode to avoid the further of accident
Diffusion.The Wide Area Measurement System (WAMS) gradually risen in recent years provides optional data source for management and running personnel.
WAMS is to acquire information for base cell with phasor measurement unit (Phasor Measurement Unit, PMU), by communication
System is uploaded to control centre, realizes the operational monitoring to electric system.PMU utilizes global positioning system (Global
Positioning System, GPS) clock synchronize the characteristics of, obtain various states of each measuring point under same time coordinate
Amount.At this stage Wide Area Measurement System monitoring data are applied to satisfy the requirements in electric network failure diagnosis.
The existing method for diagnosing faults based on WAMS data for having put into actual motion is mainly based upon pattern-recognition
Thinking.In a broad sense, it is present in the things of observable in time and space, if we can distinguish their whether identical or phases
Seemingly, " pattern " can be referred to as.And thering is pattern to compare, be registrated with oneself object observing, the process for judging its generic is exactly
Pattern-recognition is (quoted from document " mode identification method general introduction ", electronic design engineering, 2012,20 (19):48-51).But initial
Method for diagnosing faults based on pattern-recognition does not consider characteristic of the variable quantity in time scale (quoted from document " based on WAMS's
Grid disturbance recognition methods ", Automation of Electric Systems, 2006,30 (5):24-28), cause the identification to single accident occur compared with
More wrong report and the problem of fail to report.In entitled " a kind of data grid fault identification method decomposed based on WAMS sequential " (specially
Profit takes into account timing information application No. is in patent application document 201210507381.3), to original fault diagnosis
Method makes improvement.But one side this method is only applicable to single failure situation, and it is short can not to successfully manage external area error in area
The time interior complex faults such as appearance simultaneously;Another aspect this method can simultaneously be perceived for same failure due to multiple measuring points
The failure and there is the fault diagnosis result of redundancy, be easy to cause when faults frequent occurs fault diagnosis result disorder, so as to
Limit the application level of the fault diagnosis functions based on WAMS data.
Invention content
In view of the drawbacks of the prior art, the purpose of the present invention is to provide one kind to be based on WAMS times section information and topology
The method for diagnosing faults of information, it is intended to solve it is existing based on the method for diagnosing faults that WAMS sequential is decomposed there are the scope of application compared with
The problem of small and fault diagnosis result is disorderly.
The present invention provides a kind of method for diagnosing faults based on WAMS times section information and topology information, this method packets
It includes:(1) three-phase voltage of each measuring point, the amplitude of three-phase current and phase-angle data in interconnected network are obtained by WAMS, wherein,
Interconnected network specifically refers to the interconnected network equipped with Wide Area Measurement System, and i installing synchronized phasor is shared in interconnected network
The measuring point of measuring unit;(2) using step (1) obtain three-phase voltage, three-phase current amplitude data according to fault diagnosis
Start-up criterion is calculated, if meeting start-up criterion, performs step (3), otherwise return to step (1);(3) it opens meeting failure
The measuring point three-phase voltage of dynamic criterion, the amplitude data of three-phase current carry out real-time characteristic amount extraction, form real-time mode vector;
(4) by the real-time mode vector that measuring point i is extracted with it is multiple in the reference mode vector set of preset known fault type
Reference mode vector is matched respectively, if successful match, illustrates that measuring point i has occurred and the reference mode vector set pair
The failure answered, and according to matching result provide with the corresponding fault diagnosis result of reference mode vector set, if it fails to match,
Then illustrate that measuring point i does not have failure;If the measuring point diagnosis for all meeting start-up criterion is completed, (5) are entered step, otherwise
Return to step (4) diagnosis residue meets the measurement point data of start-up criterion;(5) according to temporal information and power network topology information to step
Suddenly the fault diagnosis result of different measuring points that (4) obtain is integrated, and obtains final fault diagnosis result.This method has original
The characteristics of reason is clear and definite, diagnosis speed is fast and accuracy is high, can quickly carry out electric network failure diagnosis, be operation of power networks dispatcher
The auxiliary reference of troubleshooting and decision is provided.
In order to achieve the above object, the present invention adopts the following technical scheme that:
A kind of method for diagnosing faults based on WAMS times section information and topology information, this method include the following steps:
(1) three-phase voltage of each measuring point, the amplitude of three-phase current and phase-angle data in interconnected network are obtained by WAMS;
Wherein, interconnected network specifically refers to the interconnected network equipped with Wide Area Measurement System, and shares i installing in interconnected network and synchronize
The measuring point of phasor measurement unit.
(2) using step (1) obtain three-phase voltage, three-phase current amplitude data according to fault diagnosis start-up criterion into
Row calculates, if meeting start-up criterion, performs step (3), otherwise return to step (1).
(3) to meeting the measuring point three-phase voltage of fault initiating criterion, the amplitude data progress real-time characteristic of three-phase current
Amount extraction forms real-time mode vector.
(4) the real-time mode vector for extracting measuring point i and the reference mode vector set of preset known fault type
In multiple reference mode vectors matched respectively, if successful match, illustrate measuring point i have occurred with the reference mode to
The corresponding failure of duration set, and according to matching result provide with the corresponding fault diagnosis result of reference mode vector set, if
It fails to match, then illustrates that measuring point i does not have failure;If the measuring point diagnosis for all meeting start-up criterion is completed, into step
Suddenly (5), otherwise return to step (4) diagnosis residue meet the measurement point data of start-up criterion.Wherein, preset known fault type
Reference mode vector set obtained by following methods:The amplitude amount of three-phase voltage, three-phase current to known fault type
Discontinuity surface decomposition when measured data carries out (three-phase voltage that will mutually in the same time, three-phase current amplitude metric data combinations of values
Into one group of data), establish the benchmark fault mode vector set comprising temporal information corresponding to specific fault type, Mei Geji
It closes comprising multiple benchmark fault mode vectors corresponding with specific fault type.
(5) fault diagnosis result of different measuring points obtained according to temporal information and power network topology information to step (4)
It is integrated, obtains final fault diagnosis result.
Wherein, the calculating formula of the fault diagnosis start-up criterion described in step (2) is
WithIn formula:ΔVφi=Vφi-VφNi(φ=A, B, C), VφiIt is current
The φ phase voltage amplitudes of i-th of PMU measuring point, VφNiFor the φ phase rated voltage amplitudes of i-th of PMU measuring point, Δ VsetFor voltage
Start-up criterion threshold value;ΔIφi=Iφi-IφNi(φ=A, B, C), Iφiφ phase current width for current i-th of PMU measuring points
Value, IφNiFor the φ phase nominal current magnitudes of i-th of PMU measuring point, Δ IsetFor current start criterion threshold value.ΔVsetWith Δ Iset
0.10 times of specified amplitude when what is chosen is normal operation before failure.
Wherein, the reference mode vector set described in step (4) is shaped like P(j)=[P1 (j )P2 (j)...Pi (j)] (i=1,
2 ..., m, m represent reference mode vector number, m=4, j=1,2 ..., f, f are chosen in the present invention and represents reference mode vector
Gather number, corresponding to fault type number), each benchmark fault mode vector P thereini (j)Comprising 7 characteristic quantities, shaped like Pi (j)=[UAi (j),UBi (j),UCi (j),IAi (j),IBi (j),ICi (j),ti (j)]T.In formula, UAi (j),UBi (j),UCi (j)Respectively characterize PMU
The symbolic variable of A, B, C three-phase voltage amplitude size of measuring point, 0 represents that normal voltage is horizontal, and 1 represents abnormal voltage levels;
IAi (j),IBi (j),ICi (j)The respectively symbolic variable of A, B, C three-phase current amplitude size of characterization PMU measuring points, 0 represents normal
Load current is horizontal, and 1 represents that abnormal current is horizontal;ti (j)For temporal characteristics amount, time during GPS pairs of current time is represented.Its
In:0:Phase voltage, the amplitude of phase current are represented between 0.8-1.2 times of load value, i.e., normal voltage (electric current) is horizontal;1:Generation
Table phase voltage, the amplitude of phase current are less than 0.8 times of load value, i.e. abnormal voltage (electric current) level;Here load value refers to event
Amplitude before barrier.
Wherein, the fault diagnosis result described in step (5), including fault moment, fault type and failure evolution process.
The present invention by real-time phasor data collected to phasor measurement unit in interconnected network carry out real-time characteristic to
The extraction of amount, and with multiple reference mode vectors in the reference mode vector set of preset known fault type according to the time
Sequence is matched, and preliminary fault diagnosis result is obtained, in conjunction with temporal information and power network topology information to different measuring points
Preliminary fault diagnosis result is integrated, and obtains final simple accident classification diagnostic result.Multiple reference mode vectors
Using ensure that this method can diagnose complex fault;Reference mode vector set pre-set with temporal information and
The illustrative simplicity for considering the accuracy that ensure that this method and diagnostic result of power network topology information.
Description of the drawings
Fig. 1 is the flow chart of the method for diagnosing faults the present invention is based on WAMS times section information and topology information.
Fig. 2 is the specific steps flow chart for the implementation steps (4) that the present invention is directed to transmission line single-phase fault.
Fig. 3 shows that troubles inside the sample space and external area error are in short-term in the transmission line single-phase fault that electric system actually occurs
The interior three-phase voltage and three-phase current amplitude waveform of transmission line side PMU records occurred in succession.Wherein, scheme (a), (b),
(c), (d), (e), (f) represent the A phase voltages amplitude of the transmission line of electricity side, B phase voltages amplitude, C phase voltages amplitude, A respectively
The waveform of phase current magnitude, B phase current magnitudes and C phase current magnitudes.
Fig. 4 shows that troubles inside the sample space and external area error are in short-term in the transmission line single-phase fault that electric system actually occurs
The interior three-phase voltage and three-phase current amplitude waveform of transmission line opposite side PMU records occurred in succession.Wherein, scheme (a),
(b), (c), (d), (e), (f) represent the A phase voltages amplitude, B phase voltages amplitude, C phase voltages of the transmission line of electricity opposite side respectively
Amplitude, A phase current magnitudes, B phase current magnitudes and C phase current magnitudes waveform.
Specific embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to the accompanying drawings and embodiments, it is right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
The present invention provides a kind of method for diagnosing faults based on WAMS times section information and topology information, this method packets
It includes:(1) three-phase voltage of each measuring point, the amplitude of three-phase current and phase-angle data in interconnected network are obtained by WAMS, wherein,
Interconnected network specifically refers to the interconnected network equipped with Wide Area Measurement System, and i installing synchronized phasor is shared in interconnected network
The measuring point of measuring unit;(2) using step (1) obtain three-phase voltage, three-phase current amplitude data according to fault diagnosis
Start-up criterion is calculated, if meeting start-up criterion, performs step (3), otherwise return to step (1);(3) it opens meeting failure
The measuring point three-phase voltage of dynamic criterion, the amplitude data of three-phase current carry out real-time characteristic amount extraction, form real-time mode vector;
(4) by the real-time mode vector that measuring point i is extracted with it is multiple in the reference mode vector set of preset known fault type
Reference mode vector is matched respectively, if successful match, illustrates that measuring point i has occurred and the reference mode vector set pair
The failure answered, and according to matching result provide with the corresponding fault diagnosis result of reference mode vector set, if it fails to match,
Then illustrate that measuring point i does not have failure;If the measuring point diagnosis for all meeting start-up criterion is completed, (5) are entered step, otherwise
Return to step (4) diagnosis residue meets the measurement point data of start-up criterion;(5) according to temporal information and power network topology information to step
Suddenly the fault diagnosis result of different measuring points that (4) obtain is integrated, and obtains final fault diagnosis result.
The present invention is specifically using method comprising the following steps:
(1) three-phase voltage of each measuring point, the amplitude of three-phase current and phase-angle data in interconnected network are obtained by WAMS;
Wherein, interconnected network specifically refers to the interconnected network equipped with Wide Area Measurement System, and shares i installing in interconnected network and synchronize
The measuring point of phasor measurement unit.
It should be noted that the three-phase voltages that measure of PMU, current data are phasors, i.e. plural number of time domain sinusoidal signal etc.
Valency representation, plural mould correspond to the amplitude of sinusoidal signal, and argument (polar form) is corresponding to the phase angle of sinusoidal signal.
(2) using step (1) obtain three-phase voltage, three-phase current amplitude data according to fault diagnosis start-up criterion into
Row calculates, if meeting start-up criterion, performs step (3), otherwise return to step (1);
It should be noted that the calculating formula of fault diagnosis start-up criterion isWithIn formula:ΔVφi=Vφi-VφNi(φ=A, B, C), VφiIt is current
The φ phase voltage amplitudes of i PMU measuring point, VφNiFor the φ phase rated voltage amplitudes of i-th of PMU measuring point, Δ VsetIt is opened for voltage
Dynamic criterion threshold value;ΔIφi=Iφi-IφNi(φ=A, B, C), IφiFor the φ phase current magnitudes of current i-th of PMU measuring points,
IφNiFor the φ phase nominal current magnitudes of i-th of PMU measuring point, Δ IsetFor current start criterion threshold value.ΔVsetWith Δ IsetChoosing
0.10 times of specified amplitude when what is taken is normal operation before failure.The criterion can effectively find three-phase voltage and three-phase current width
The variation of value ensures the high sensitivity of the method for diagnosing faults.
When meeting at least one in above-mentioned two calculating formula, then step (3) is performed;Otherwise return to step (1).
(3) to meeting the measuring point three-phase voltage of fault initiating criterion, the amplitude data progress real-time characteristic of three-phase current
Amount extraction forms real-time mode vector.
It is carried it should be noted that only measuring point data to the PMU for meeting fault initiating criterion here and carrying out real-time characteristic amount
It takes.It is assumed that PMU measuring points k meets fault initiating criterion, the GPS moment is t at this time0, the real-time characteristic of measuring point k is extracted immediately
Amount forms real-time mode vector Xk=[UAk,UBk,UCk,IAk,IBk,ICk,tk]T.In formula, UAk,UBk,UCkRespectively characterize PMU amounts
The symbolic variable of A, B, C three-phase voltage amplitude size of measuring point k, 0 represents that normal voltage is horizontal, and 1 represents abnormal voltage levels;
IAk,IBk,ICkThe respectively symbolic variable of A, B, C three-phase current amplitude size of characterization PMU measuring points k, 0 represents normal duty
Levels of current, 1 represents that abnormal current is horizontal;tkFor temporal characteristics amount, time during GPS pairs of current time is represented.Wherein:
0:Phase voltage, the amplitude of phase current are represented between 0.8-1.2 times of load value, i.e. normal voltage (electric current) water
It is flat;
1:Phase voltage, the amplitude of phase current are represented less than 0.8 times of load value, i.e. abnormal voltage (electric current) level;
Here load value refers to the amplitude before failure.
Real-time mode vector X can be obtained according to above rulek, and t at this timek=t0.For same PMU measuring points k,
In t=t0+ Δ t to t=t0+tdPeriod calculates the real-time mode vector at corresponding moment every a time interval Δ t.Wherein,
Δ t be sampling time interval, Δ t=40 milliseconds here;tdFor preset trouble duration, td=20 seconds here.By upper
Calculating is stated, can obtain one group of real-time mode vector
(4) the real-time mode vector for extracting measuring point i and the reference mode vector set of preset known fault type
In multiple reference mode vectors matched respectively, if successful match, illustrate measuring point i have occurred with the reference mode to
The corresponding failure of duration set, and according to matching result provide with the corresponding fault diagnosis result of reference mode vector set, if
It fails to match, then illustrates that measuring point i does not have failure;If the measuring point diagnosis for all meeting start-up criterion is completed, into step
Suddenly (5), otherwise return to step (4) diagnosis residue meet the measurement point data of start-up criterion.Wherein, preset known fault type
Reference mode vector set obtained by following methods:The amplitude amount of three-phase voltage, three-phase current to known fault type
Discontinuity surface decomposition when measured data carries out (three-phase voltage that will mutually in the same time, three-phase current amplitude metric data combinations of values
Into one group of data), establish the benchmark fault mode vector set comprising temporal information corresponding to specific fault type, Mei Geji
It closes comprising multiple benchmark fault mode vectors corresponding with specific fault type.
It should be noted that step (4) includes following specific implementation step:
A) vectorial first element of real-time mode is chosenSuccessively with the first of f real-time reference pattern vector set
A element P1 (j)(j=1,2 ..., f) is matched, if with x-th of vector P1 (x)Successful match, then being chosen in subsequent match should
Real-time reference pattern vector set P where vectorial(x)Remaining reference mode vector matched, and will real-time mode at this time
The last one characteristic quantity (temporal characteristics amount) of vector is assigned to P1 (x)Temporal characteristics amount t1 (x), enter step b), if can not
Match, then choose next real-time mode vector, repeat step a), until completing the last one real-time mode vector
Match;
B) choose next real-time mode vector, and with real-time reference pattern vector set P(x)Second element P2 (x)
It is matched, if successful match, the last one characteristic quantity (temporal characteristics amount) of the vector of real-time mode at this time is assigned to P2 (x)'s
Temporal characteristics amount t2 (x), step b) is repeated, if it fails to match, calculates t2 (x)-t1 (x)And with preset time parameter Tset1 (x)Compare
Size, if t2 (x)-t1 (x)>Tset1 (x)It then enters step c), if t2 (x)-t1 (x)≤Tset1 (x)Then terminate this calculating, output and institute
Corresponding fault diagnosis result (the benchmark failure mould of matched basic schema vector total number in matched reference mode vector set
Formula vector set and corresponding fault diagnosis result are ready in advance), return to step (1);
C) choose next real-time mode vector, and with real-time reference pattern vector set P(x)Third element P3 (x)
It is matched, if successful match, the last one characteristic quantity (temporal characteristics amount) of the vector of real-time mode at this time is assigned to P3 (x)'s
Temporal characteristics amount t3 (x), step c) is repeated, if it fails to match, calculates t3 (x)-t2 (x)And with preset time parameter Tset2 (x)Compare
Size, if t3 (x)-t2 (x)>Tset2 (x)It then enters step d), if t3 (x)-t2 (x)≤Tset2 (x)Then terminate this calculating, output and institute
Corresponding fault diagnosis result (the benchmark failure mould of matched basic schema vector total number in matched reference mode vector set
Formula vector set and corresponding fault diagnosis result are ready in advance), return to step (1);
D) choose next real-time mode vector, and with real-time reference pattern vector set P(x)The 4th element P4 (x)
It is matched, if successful match, the last one characteristic quantity (temporal characteristics amount) of the vector of real-time mode at this time is assigned to P4 (x)'s
Temporal characteristics amount t4 (x), step d) is repeated, if it fails to match, calculates t4 (x)-t3 (x)And with preset time parameter Tset3 (x)Compare
Size, if t4 (x)-t3 (x)>Tset3 (x)It then enters step e), if t4 (x)-t3 (x)≤Tset3 (x)Then terminate this calculating, output and institute
Corresponding fault diagnosis result (the benchmark failure mould of matched basic schema vector total number in matched reference mode vector set
Formula vector set and corresponding fault diagnosis result are ready in advance), return to step (1);
If real-time mode vector e) is traversedWhole elements, then terminate
Epicycle calculates, and output failure corresponding with basic schema vector total number matched in matched reference mode vector set is examined
Disconnected result (benchmark fault mode vector set and corresponding fault diagnosis result are ready in advance) is simultaneously extracted next full
The real-time mode vector of the measuring point of start-up criterion described in sufficient step (3) repeats step (4), otherwise returns to previous step.
Wherein, time parameter Tset1 (x)、Tset2 (x)、Tset3 (x)Be according to the preset numerical value of known fault type,
Number is (m-1) a (m represents reference mode vector number).
Wherein, step b) to step d) increases with the reference mode vector number m in corresponding reference mode vector set
Subtract, m=4 in this patent, therefore shared (m-1=3) 3 circulation steps (i.e. step b) to step d)).
It should be noted that the reference mode vector number m in the reference mode vector set of known fault type can root
Factually border fault type is increased and decreased, and can also be increased and decreased with the corresponding fault diagnosis result of reference mode vector set.
(5) fault diagnosis result of different measuring points obtained according to temporal information and power network topology information to step (4)
It is integrated, obtains final fault diagnosis result.
It should be noted that the integration work of the fault diagnosis result of different measuring points is carried out there are two foundation, first, not
With the time of failure gap size of failure, second is that the Topology connection situation of the corresponding PMU measuring points of different faults.If two
The fault diagnosis result of event of failure meets:Time of failure interval is less than preset time interval Tint, and two failure things
Two event of failure are then integrated into an event of failure by the corresponding PMU measuring points of part there are during direct Topology connection;Otherwise
Think that two event of failure can not be integrated.All event of failure are integrated one by one, obtain the failure of event of failure to the end
Diagnostic result.
Wherein, TintFor the preset minimum interval for distinguishing two event of failure, value is by actual accidents data
Arrangement obtains, and T is chosen in this patentint=5s.
Be below apply the present invention in transmission line single-phase fault troubles inside the sample space and in the external area error short time in succession
The case study on implementation of generation.Transmission line single-phase fault is directed to, corresponding reference mode vector set and time parameter are shown in Table 1 institute
Show, corresponding fault diagnosis output the results are shown in Table 2.
Table 1
Table 2
Fig. 2 show the specific steps flow chart for the implementation steps (5) for being directed to transmission line single-phase fault.Fig. 3 is shown
Troubles inside the sample space and occur in succession in the external area error short time in the transmission line single-phase fault that electric system actually occurs
The three-phase voltage and three-phase current amplitude waveform of transmission line side PMU records.Fig. 4 shows that electric system actually occurs primary
Troubles inside the sample space and the transmission line opposite side PMU occurred in succession in the external area error short time are recorded in transmission line single-phase fault
Three-phase voltage and three-phase current amplitude waveform.Apply the present invention on the actual measurement PMU data, the data and circuit can be obtained
The reference mode vector sets match of the B phase faults in single-phase fault is transmitted, corresponding benchmark has been marked in attached drawing 3 and attached drawing 4
The match condition of pattern vector.Table 2 is corresponded to according to matched result, the fault diagnosis result that can obtain attached drawing 3 is:The event of B phases
Barrier, fault moment:13:34:56.960 unsuccessful reclosing, Reclosing Instant:13:34:57.960;The fault diagnosis result of Fig. 4
For:B phase faults, fault moment:13:34:56.120 unsuccessful reclosing, Reclosing Instant:13:34:57.760.Two PMU amounts
In the fault diagnosis result that measuring point obtains, 0.84s is divided between fault moment, less than Tint=5s, and two PMU measuring points measure
Be same transmission line of electricity both ends, there are the connection of direct electric topology, therefore the two event of failure can be integrated into:B phases
Failure, fault moment:13:34:56.960 unsuccessful reclosing, Reclosing Instant:13:34:57.960.
Similar, other failures that transmission line of electricity, power plant and DC converter station occur can also be used and this implementation
The similar method of example carries out fault diagnosis.
The present invention by real-time phasor data collected to phasor measurement unit in interconnected network carry out real-time characteristic to
The extraction of amount, and with multiple reference mode vectors in the reference mode vector set of preset known fault type according to the time
Sequence is matched, and preliminary fault diagnosis result is obtained, in conjunction with temporal information and power network topology information to different measuring points
Preliminary fault diagnosis result is integrated, and obtains final simple accident classification diagnostic result.Multiple reference mode vectors
Using ensure that this method can diagnose complex fault;Reference mode vector set pre-set with temporal information and
The illustrative simplicity for considering the accuracy that ensure that this method and diagnostic result of power network topology information.
As it will be easily appreciated by one skilled in the art that the foregoing is merely illustrative of the preferred embodiments of the present invention, not to
The limitation present invention, all any modification, equivalent and improvement made all within the spirits and principles of the present invention etc., should all include
Within protection scope of the present invention.
Claims (5)
1. a kind of method for diagnosing faults based on WAMS times section information and topology information, which is characterized in that including following step
Suddenly:
(1) amplitude of three-phase voltage of each measuring point and phase-angle data, the width of three-phase current in interconnected network are obtained by WAMS
Value and phase-angle data;Wherein, interconnected network refers to the interconnected network equipped with Wide Area Measurement System, and in the interconnected network altogether
There is the measuring point of n installing synchronous phasor measurement unit;
(2) to the amplitude data of three-phase voltage of each measurement point, the amplitude data of three-phase current according to fault diagnosis start-up criterion
It is calculated, if meeting start-up criterion, performs step (3), otherwise return to step (1);
(3) amplitude data of the three-phase voltage of measuring point, the amplitude data of three-phase current that meet fault initiating criterion are carried out
Real-time characteristic amount is extracted, and forms real-time mode vector;
(4) in the reference mode vector set of the real-time mode vector for extracting measuring point i and preset known fault type
Multiple reference mode vectors are matched respectively, if successful match, measuring point i has occurred and the reference mode vector set pair
The failure answered, and provided and the corresponding fault diagnosis result of reference mode vector set according to matching result;If it fails to match,
Then measuring point i does not have failure;Enter step (5) after the measuring point for all meeting start-up criterion diagnoses completion;
(5) fault diagnosis result of different measuring points that step (4) obtains is carried out according to temporal information and power network topology information
It integrates, obtains final fault diagnosis result.
2. method for diagnosing faults as described in claim 1, which is characterized in that the fault diagnosis startup described in step (2) is sentenced
According to forWith
Wherein, Δ ViThe virtual value of three-phase voltage deviation for current i-th of PMU measuring points, Δ IiFor current i-th of PMU measuring points three
The virtual value of phase current deviation, Δ Vφi=Vφi-VφNi(φ=A, B, C), Vφiφ phase voltage width for current i-th of PMU measuring points
Value, VφNiFor the φ phase rated voltage amplitudes of i-th of PMU measuring point, Δ VsetFor voltage starting criterion threshold value;ΔIφi=
Iφi-IφNi(φ=A, B, C), IφiFor the φ phase current magnitudes of current i-th of PMU measuring points, IφNiFor i-th PMU measuring point
φ phase nominal current magnitudes, Δ IsetFor current start criterion threshold value, Δ VsetWith Δ IsetThe normal operation before being failure chosen
When 0.10 times of specified amplitude.
3. method for diagnosing faults as claimed in claim 1 or 2, which is characterized in that in step (4), the reference mode to
Duration set P(j)=[P1 (j) P2 (j) ... Pi (j)], i=1,2 ..., m, m represent reference mode vector number, j=1,2 ...,
F, f represent reference mode vector set number, corresponding to fault type number, each benchmark fault mode vector Pi (j)Include 7
Characteristic quantity Pi (j)=[UAi (j),UBi (j),UCi (j),IAi (j),IBi (j),ICi (j),ti (j)]T;UAi (j),UBi (j),UCi (j)Respectively characterize
The symbolic variable of A, B, C three-phase voltage amplitude size of PMU measuring points;IAi (j),IBi (j),ICi (j)Respectively characterize PMU measuring points
A, B, C three-phase current amplitude size symbolic variable;ti (j)For temporal characteristics amount, time during GPS pairs of current time is represented.
4. method for diagnosing faults as claimed in claim 1 or 2, which is characterized in that described preset known in step (4)
The reference mode vector set of fault type is obtained by following methods:
Discontinuity surface when the amplitude metric data of amplitude metric data, three-phase current to the three-phase voltage of known fault type carries out
It decomposes, establishes the benchmark fault mode vector set comprising temporal information corresponding to specific fault type, each set includes
Multiple benchmark fault mode vectors corresponding with specific fault type.
5. method for diagnosing faults as claimed in claim 1 or 2, which is characterized in that in step (5), the final failure is examined
Disconnected result includes fault moment, fault type and failure evolution process.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610127868.7A CN105759169B (en) | 2016-03-07 | 2016-03-07 | A kind of method for diagnosing faults based on WAMS times section information and topology information |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610127868.7A CN105759169B (en) | 2016-03-07 | 2016-03-07 | A kind of method for diagnosing faults based on WAMS times section information and topology information |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105759169A CN105759169A (en) | 2016-07-13 |
CN105759169B true CN105759169B (en) | 2018-07-06 |
Family
ID=56331814
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610127868.7A Active CN105759169B (en) | 2016-03-07 | 2016-03-07 | A kind of method for diagnosing faults based on WAMS times section information and topology information |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105759169B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112766426B (en) * | 2021-04-06 | 2021-09-07 | 中国铁道科学研究院集团有限公司通信信号研究所 | Fault type identification method based on time constraint |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2249226C2 (en) * | 2003-02-04 | 2005-03-27 | Бодров Павел Александрович | Method of measuring remoteness of single-phase short circuit in three-phase power line |
CN101661069A (en) * | 2009-09-25 | 2010-03-03 | 北京四方继保自动化股份有限公司 | Dynamic process real-time estimation method of weak observable non-PMU measuring point not depending on state matrix |
CN101750562A (en) * | 2010-01-13 | 2010-06-23 | 湖北省电力公司 | Non-PMU measure point dynamic process estimation method based on flow equation sensitiveness analysis |
CN103760400A (en) * | 2014-01-08 | 2014-04-30 | 东北电力大学 | Voltage stabilization on-line monitoring method based on wide area measuring system |
CN105044554A (en) * | 2015-06-30 | 2015-11-11 | 安徽博诺思信息科技有限公司 | Power grid fault detection and diagnosis method |
-
2016
- 2016-03-07 CN CN201610127868.7A patent/CN105759169B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2249226C2 (en) * | 2003-02-04 | 2005-03-27 | Бодров Павел Александрович | Method of measuring remoteness of single-phase short circuit in three-phase power line |
CN101661069A (en) * | 2009-09-25 | 2010-03-03 | 北京四方继保自动化股份有限公司 | Dynamic process real-time estimation method of weak observable non-PMU measuring point not depending on state matrix |
CN101750562A (en) * | 2010-01-13 | 2010-06-23 | 湖北省电力公司 | Non-PMU measure point dynamic process estimation method based on flow equation sensitiveness analysis |
CN103760400A (en) * | 2014-01-08 | 2014-04-30 | 东北电力大学 | Voltage stabilization on-line monitoring method based on wide area measuring system |
CN105044554A (en) * | 2015-06-30 | 2015-11-11 | 安徽博诺思信息科技有限公司 | Power grid fault detection and diagnosis method |
Also Published As
Publication number | Publication date |
---|---|
CN105759169A (en) | 2016-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kezunovic | Smart fault location for smart grids | |
Almas et al. | RT-HIL implementation of the hybrid synchrophasor and GOOSE-based passive islanding schemes | |
Xu et al. | Synchronized waveforms–a frontier of data-based power system and apparatus monitoring, protection, and control | |
CN109582004B (en) | Intelligent detection platform | |
Bastos et al. | Synchrowaveform measurement units and applications | |
CN103618385B (en) | One puies forward high-precision state estimation data correction system and method | |
CN110783946A (en) | Method for locating phase faults in microgrid | |
CN109283407B (en) | Voltage loop monitoring system based on total station data contrastive analysis | |
Guo et al. | A statistical process control approach for automatic anti-islanding detection using synchrophasors | |
Rinaldi et al. | A testing framework for the monitoring and performance analysis of distributed energy systems | |
Zamora-Cárdenas et al. | Hybrid state estimator considering SCADA and synchronized phasor measurements in VSC-HVDC transmission links | |
CN103905045B (en) | Alternating-current analog-digital conversion device real-time response detection method based on physical layer sampled value | |
CN105759169B (en) | A kind of method for diagnosing faults based on WAMS times section information and topology information | |
Naumann et al. | Experience with PMUs in industrial distribution networks | |
CN108693442A (en) | The detection method and system of electric network fault | |
Mondal et al. | Application of PMU and SCADA data for estimation of source of forced oscillation | |
Gomathi et al. | Optimal location of PMUs for complete observability of power system network | |
CN109768619A (en) | Monitoring information automatic debugging system and adjustment method | |
Tlusty et al. | The monitoring of power system events on transmission and distribution level by the use of phasor measurements units (PMU) | |
CN105067949B (en) | Method for synchronizing time and system based on circuit both end voltage phase contrast | |
Chandra et al. | PMU based real time power system state estimation using ePHASORsim | |
Wang et al. | Real time fault monitoring and diagnosis method for power grid monitoring and its application | |
CN114287093A (en) | Method and system for estimating and monitoring power distribution network conditions | |
CN111813825A (en) | Distribution transformer outlet power failure event missing report automatic detection method and system | |
Blas et al. | Micro-Synchrophasors for Event Detection in an Electric Power System Simulator |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |