CN103487723A - Electric system fault diagnosis method and system - Google Patents

Electric system fault diagnosis method and system Download PDF

Info

Publication number
CN103487723A
CN103487723A CN201310388672.XA CN201310388672A CN103487723A CN 103487723 A CN103487723 A CN 103487723A CN 201310388672 A CN201310388672 A CN 201310388672A CN 103487723 A CN103487723 A CN 103487723A
Authority
CN
China
Prior art keywords
fault
protection
fault diagnosis
power system
warning information
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201310388672.XA
Other languages
Chinese (zh)
Other versions
CN103487723B (en
Inventor
陈亦平
周华锋
李矛
赵旋宇
熊卫斌
文福拴
吴文可
李晓露
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Alstom Electric Power Network Technique Center Co Ltd
Zhejiang University ZJU
China Southern Power Grid Co Ltd
Original Assignee
Alstom Electric Power Network Technique Center Co Ltd
Zhejiang University ZJU
China Southern Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Alstom Electric Power Network Technique Center Co Ltd, Zhejiang University ZJU, China Southern Power Grid Co Ltd filed Critical Alstom Electric Power Network Technique Center Co Ltd
Priority to CN201310388672.XA priority Critical patent/CN103487723B/en
Publication of CN103487723A publication Critical patent/CN103487723A/en
Application granted granted Critical
Publication of CN103487723B publication Critical patent/CN103487723B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention provides an electric system fault diagnosis method and system. The method includes the following steps of receiving alarm information generated after electric system faults appear, determining a fault area according to the alarm information, conducting modeling on fault components according to the logic relation between each fault component in the fault area and the motion of a corresponding protective breaker, setting up a weighting fuzzy time sequence Petri network fault diagnosis model of the fault components, and setting up a corresponding matrix according to reasoning analysis of the weighting fuzzy time sequence Petri network fault diagnosis model to conduct reasoning calculation to diagnose the fault components, wherein the alarm information includes time sequence information. According to the electric system fault diagnosis method and system, the fault tolerance and the accuracy of the fault diagnosis in an electric system are effectively improved.

Description

Power system failure diagnostic method and system
Technical field
The present invention relates to the power system security processing technology field, particularly relate to a kind of power system failure diagnostic method and a kind of power system fault diagnosis.
Background technology
Power system failure diagnostic is exactly that the warning information produced after utilizing fault to occur is located fault element quickly and accurately, for the dispatcher provides decision support, thereby rapidly removing faults, recovery system, to normal operating condition, improves the operation stability of system.Done a large amount of research work both at home and abroad on the power system failure diagnostic direction, proposed the different fault diagnosis method, wherein simple, the physical significance of model clearly Petri net method for diagnosing faults be widely applied.
Traditional application Petri net carries out in the method for power system failure diagnostic, for protection and isolating switch, the uncertain problems such as tripping or malfunction may occur, and improves the fault-tolerance of fault diagnosis by setting up the Fuzzy Petri Net inference pattern.Yet said method only utilizes the protection that receives and the action message of isolating switch, when complex fault occurring and follow protection and isolating switch tripping or malfunction and warning information self distorted, may cause obtaining correct diagnostic result.
Summary of the invention
Based on this, the invention provides a kind of power system failure diagnostic method and system, can improve the accuracy of fault diagnosis.
For achieving the above object, the present invention adopts following technical scheme:
A kind of power system failure diagnostic method comprises the following steps:
Receive after electric power system fault occurs the warning information produced, described warning information comprises time sequence information;
Determine fault zone according to described warning information;
To each fault element modeling, set up the Weighted Fuzzy temporal Petri nets fault diagnosis model of fault element according to the logical relation between each fault element in described fault zone and corresponding protection, isolating switch action;
According to the rational analysis of described Weighted Fuzzy temporal Petri nets fault diagnosis model, construct corresponding matrix and carry out the reasoning computing, diagnose the element that is out of order.
A kind of power system fault diagnosis comprises:
The warning information receiver module, for receiving after electric power system fault occurs the warning information produced, described warning information comprises time sequence information;
The fault zone determination module, for determining fault zone according to described warning information;
Model building module, to each fault element modeling, set up the Weighted Fuzzy temporal Petri nets fault diagnosis model of fault element for the logical relation between each fault element according to described fault zone and corresponding protection, isolating switch action;
The reasoning computing module, for the rational analysis according to described Weighted Fuzzy temporal Petri nets fault diagnosis model, construct corresponding matrix and carry out the reasoning computing, diagnoses the element that is out of order.
By above scheme, can be found out; a kind of power system failure diagnostic method of the present invention and system; the time delay binding feature of existence between protection and isolating switch action and the possibility of protection and isolating switch malfunction and tripping have been considered; propose a kind of electric system Weighted Fuzzy temporal Petri nets fault diagnosis model that can take into account this time delay constraint on the basis of existing fault diagnosis model, and carried out fault diagnosis according to this Weighted Fuzzy temporal Petri nets fault diagnosis model.Because fault diagnosis model of the present invention can be processed time delay restricted problem, therefore effectively improved fault-tolerance and the accuracy of fault diagnosis in the electric system.
The accompanying drawing explanation
The schematic flow sheet that Fig. 1 is a kind of power system failure diagnostic method in the embodiment of the present invention;
Fig. 2 is the fault diagnosis model schematic flow sheet based on the Weighted Fuzzy temporal Petri nets in the embodiment of the present invention;
The Weighted Fuzzy temporal Petri nets fault diagnosis model schematic diagram that Fig. 3 is embodiment of the present invention median generatrix;
The Weighted Fuzzy temporal Petri nets fault diagnosis model schematic diagram that Fig. 4 is circuit in the embodiment of the present invention;
Fig. 5 is IEEE New England 10 machine 39 node system schematic diagram in the embodiment of the present invention;
The structural representation that Fig. 6 is a kind of power system fault diagnosis in the embodiment of the present invention.
Embodiment
Shown in Figure 1, a kind of power system failure diagnostic method comprises the following steps:
Step S101, receive after electric power system fault occurs the warning information produced, and described warning information comprises time sequence information, then enters step S102.
Step S102, determine fault zone according to described warning information, then enters step S103.
As an embodiment preferably, the process of described definite fault zone specifically can comprise: adopt breadth-first search to determine described fault zone.
Step S103, to each fault element modeling, set up the Weighted Fuzzy temporal Petri nets fault diagnosis model of fault element according to the logical relation between each fault element in described fault zone and corresponding protection, isolating switch action, then enters step S104.
Step S104, according to the rational analysis of described Weighted Fuzzy temporal Petri nets fault diagnosis model, construct corresponding matrix and carry out the reasoning computing, diagnoses the element that is out of order.
As an embodiment preferably, after element is out of order in diagnosis, can also comprise the steps:
Step S105, carry out backward reasoning to the fault element of diagnosing out, draws malfunction and the tripping situation of protection and isolating switch.
1), determine power supply interrupted district (being fault zone) in another embodiment, fault diagnosis model of the present invention can be divided into following four layers:; 2), set up the element fault diagnostic model; 3), the reasoning based on matrix computations; 4), protection and isolating switch action evaluation.Detailed process is shown in Figure 2:
1, the isolating switch displacement information provided according to network topology structure and data acquisition and surveillance (supervisory control and data acquisition, SCADA), adopt BFS (Breadth First Search) to determine fault zone.If fault zone only comprises an element, this element is fault element; If comprise two or more elements in fault zone, each element is wherein set up respectively to Weighted Fuzzy temporal Petri nets network model and carry out the failure judgement element;
2., according to the logical relation between fault element and corresponding protection, isolating switch action, build the Weighted Fuzzy temporal Petri nets diagnostic model of element, then merged according to each interelement topological relation;
3, the Weighted Fuzzy temporal Petri nets diagnostic model based on developed, realize rational analysis by matrix operation, the final diagnosis element that is out of order;
4, the fault element of diagnosing out is taked to backward reasoning, with malfunction and the tripping situation of judgement protection and isolating switch.
Below troubleshooting step of the present invention is elaborated.
One, the mathematical description of Weighted Fuzzy temporal Petri nets:
1), definition Weighted Fuzzy temporal Petri nets is 11 tuples:
S WFTPN={P,T,I,O,A cc,ΔT min,ΔT max,U,T h,W,M} (1)
In formula: P={p 1, p 2..., p nby storehouse, collected; T={t 1, t 2..., t mit is the transition collection; For characterizing inference rule; I:P → T is the mapping of reflection transition that storehouse is arrived; I=[δ ij] be n * m matrix; Work as p it jinput (have p ito t jdirected arc) time δ ij=1; Otherwise δ ij=0; O:T → P reflection be transitted towards storehouse mapping; O=[γ ij] be m * n matrix; Work as p jt ioutput (have t ito p jdirected arc) time γ ij=1; Otherwise γ ij=0; A cc=[a ij] be n * n matrix; Characterize general storehouse arrive the purpose storehouse path; Work as p istorehouse institute path through p jthe time a ij=1; Otherwise a ij=0; Δ T min=[Δ τ 1min, Δ τ 2min..., Δ τ nmin] be storehouse with the minimum time delay constraint of preposition transition; Δ T max=[Δ τ 1max, Δ τ 2max,, Δ τ nmax] be storehouse with the constraint of the maximum delay of preposition transition; If Δ τ min=Δ τ max=0; Transition moment activation; U=[μ 1, μ 2..., μ m] for the degree of confidence vector of transition; If for any j, μ is arranged j=1; Model is the simple Petri net that does not contain fuzzy variable; T h=[λ 1, λ 2..., λ m] for the igniting threshold vector of transition; W=diag (w 1, w 2..., w n) for inputting the weight matrix of arc; The reflection precondition is to regular influence degree; The event type that its value characterizes to storehouse is relevant; M=[α (p 1), α (p 2) ..., α (p n)] be storehouse institute degree of confidence vector; α (p i) p of library representation institute idegree of confidence.
The Petri net can pass through the node P={p of storehouse institute 1, p 2..., p n, transition node T={t 1, t 2..., t mand directed arc mean.Directed arc comprises input arc and output arc two classes, and the input arc is by storehouse transition pointed, and the output arc points to storehouse institute by transition.Wherein, storehouse circle used " zero " means, vertical line for transition " | " expression.The structure of Weighted Fuzzy Petri Net and the maximum difference of simple Petri net be to consider the input arc weights, storehouse degree of confidence, the degree of confidence of transition, the transition process of probable value etc.In addition, the probable value of the event of institute's value representation storehouse, Weighted Fuzzy Petri Net storehouse representative, even transition are activated, input magazine is worth and remains unchanged, library representation institute degree of confidence still, and can be along with transition produce storehouse institute.Consider that on the basis of Weighted Fuzzy Petri Net model the time delay between the institute of storehouse retrains, and is Weighted Fuzzy temporal Petri nets model.This model can be processed time delay restricted problem, can improve fault-tolerance and the accuracy of model.
2), the rational analysis of transition temporal constraint, forward direction collection and the backward collection of transition t ∈ T are defined as respectively .t={p| δ pt=1}; t .={ p| γ tp=1}.
Component library action delay constraint be to being connected to same transition t both sides Liang Ge storehouse institutes (successive relay trip) .t and t .the time delay interval
Figure BDA00003750885700051
with
Figure BDA00003750885700052
carry out interconnection constraint, i.e. Δ T t=[Δ τ tmin, Δ τ tmax], be embodied in these two component libraries the connection transition on.
Forward inference: known
Figure BDA00003750885700053
with Δ T t=[Δ τ tmin, Δ τ tmax], can obtain component library institute .the follow-up component library of t action delay constraint:
T t . = T t . + Δ T t = [ τ t min . , τ t max . ] + [ Δτ t min , Δτ t max ] = [ τ t min . + Δτ t min . , τ t max . + Δ τ t max ] - - - ( 2 )
Backward reasoning: known
Figure BDA00003750885700055
with Δ T t=[Δ τ tmin, Δ τ tmax], can obtain the t of component library institute .forerunner's component library action delay constraint:
T t . = T t . - Δ T t = [ τ t . min , τ t . max ] - [ Δ τ t min , Δτ t max ] = [ τ t . min - Δτ t max , τ t . max - Δτ t min ] - - - ( 3 )
3), for Weighted Fuzzy Petri Net, its fuzzy reasoning depends on fuzzy rule set R={R 1, R 2..., R m}:
R i:Ifd jThend k(CF=μ i)(i=1,2,…,m) (4)
In formula: d jand d kfor the proposition containing fuzzy variable, its confidence alpha (d j) and α (d k) value in interval [0,1], in the Petri net, be presented as storehouse degree of confidence; μ i∈ [0,1] shows regular R ithe degree of confidence value, it is presented as the degree of confidence of transition in Petri net.
In the reasoning process of Weighted Fuzzy temporal Petri nets, need to consider the time delay constraint simultaneously.Suppose that A, B and C are m * n matrix, and D is m * q matrix, E is q * n matrix; Definition: (1) addition operator ⊕: C=A ⊕ B, c ij=max (a ij, b ij); (2) comparison operator
Figure BDA000037508857000512
: work as a ij>=b ijthe time c ij=1, otherwise c ij=0; (3) directly take advantage of c of operator e:C=AeB ij=a ijb ij; (4) multiplication operator
Figure BDA000037508857000510
:
Figure BDA000037508857000511
Figure BDA00003750885700057
(5) matrix multiplication g:C=DgE,
Figure BDA00003750885700058
The reasoning process of Weighted Fuzzy temporal Petri nets is for obtaining a stable network state, and the value of storehouse institute degree of confidence matrix M is no longer carried out with iteration and the state that changes.At first retrain alarm is screened by time delay.Suppose that the k time iteration obtains the degree of confidence matrix M k, obtain the degree of confidence matrix M the k+1 time k+1reasoning process as follows:
Step1: definition storehouse institute path matrix is A cc.According to temporal constraint, the minimum cumulative delay constraint matrix Σ Δ T that storehouse institute path is corresponding is obtained in rational analysis min, and cumulative maximum time delay constraint matrix Σ Δ T max, be all n * 1 vector:
ΣΔT min=(A ccg(ΔT min) T) T (5)
ΣΔT max=(A ccg(ΔT maxT) T (6)
Step2: the sequential consistance of check warning information, screen on this basis alarm.Definition F=[f 1, f 2..., f n] be the time delay constraint qualification vector of protection and isolating switch, element wherein characterizes time delay that whether corresponding storehouse meets that model bank is concentrated constraint, f i=1 and 0 (i=1; 2; N) mean respectively to meet and do not meet constraint.Delayed data set Δ T by protection and isolating switch in warning information mesminwith Δ T mesmaxcompare with the time delay constraint, shown in (7), can obtain:
Figure BDA00003750885700061
Step3: given original state M 0, make k=0;
Step4: the weights that calculate the input arc:
W inarc=WgI (8)
Step5: the synthetic input confidence level of calculating transition:
E k=M kgW inarc (9)
Step6: synthetic input confidence level and the threshold value of transition are compared; Be met the transition set of activation condition:
Figure BDA00003750885700062
Step7: the synthetic input confidence level that calculating can make transition activate:
H k=E keG k (11)
Step8: the M that calculates the k+1 time reckoning of storehouse institute k+1:
Figure BDA00003750885700063
Step9: if M k+1=M k, the degree of confidence matrix of Petri net is stable, otherwise makes k=k+1, returns to Step4.
Two, Weighted Fuzzy temporal Petri nets fault diagnosis model, after the element generation permanent fault in electric system, relevant protection and isolating switch can move, and finally form one or more power supply interrupted districts, and fault element is certainly among power supply interrupted district.The fault diagnosis model below developed is exactly for element included in power supply interrupted district.
The IEEE39 node system of take is set up the Weighted Fuzzy temporal Petri nets fault diagnosis model of bus and circuit as example.Suppose that every circuit two ends have configured main protection, nearly back-up protection and back-up protection far away.Back-up protection on the main protection of having supposed every bus bar configuration, the circuit that is connected.Mean respectively bus, circuit, main protection, nearly back-up protection and back-up protection far away with B, L, m, p and s.For example, CB (4)-14mean the isolating switch of circuit L4-14 near bus B4 side, R (4)-14mmean the main protection of circuit L4-14 near bus B4 side, the rest may be inferred.
1), the protection sequence of movement is followed successively by main protection, nearly back-up protection and back-up protection far away.Defining them is respectively with respect to the time delay of fault moment
Figure BDA00003750885700071
with
Figure BDA00003750885700072
the isolating switch that defines all kinds of protection correspondences is respectively with respect to the time delay of operating time of protection Δ T mc ∈ [ Δτ mc min , Δτ mc max ] , Δ T pc ∈ [ Δτ pc min , Δτ pc max ] With Δ T sc ∈ [ Δτ sc min , Δτ sc max ] .
2), basic ideas that bus is set up to fault diagnosis model are: the model of its each closure main protection of model, back-up protection, then build bus resultant fault diagnostic model.Similarly, the basic ideas of circuit being set up to fault diagnosis model are: all kinds of protection models at its two ends of model, then build circuit resultant fault diagnostic model.
The bus B14 of take in Fig. 5 is example, and it has configured main protection R b14mand the back-up protection R far away of connection line (4)-14s, R (13)-14sand R 14-(15) s, obtain the Weighted Fuzzy temporal Petri nets fault diagnosis unified model of this bus as shown in Figure 3.Take circuit L4-14 as example, and it has configured main protection R (4)-14mand R 4-(14) m, nearly back-up protection R (4)-14pand R 4-(14) p, back-up protection R far away (3)-4s, R 4-(5) s, R 14-(15) sand R (13)-14s, obtain the Weighted Fuzzy temporal Petri nets fault diagnosis unified model of this circuit as shown in Figure 4.
3) parameter determination method, after having built fault diagnosis model, next step just is based on the reasoning process of matrix operation, and it is consistent with the fuzzy reasoning based on fuzzy rule.Wherein, relevant parameters is determined as follows:
1. the arc weights are inputted in transition
For each end of circuit, consider main protection, nearly back-up protection and the back-up protection far away Different Effects to circuit, with and with the coordinating of respective circuit breakers.Like this, input arcs for protection storehouse institute, isolating switch storehouse institute to two of transition and give respectively weights, as shown in table 1.
The protection that table 1 is given and isolating switch weights
Figure BDA00003750885700076
2. the degree of confidence initial value is given
To protection and the corresponding storehouse of isolating switch alarm institute, given its degree of confidence initial value, as shown in table 2.
The degree of confidence of the protection that table 2 is given and isolating switch alarm
Figure BDA00003750885700081
Consider that the warning information received may be mistake or incomplete, for the protection in alarm and isolating switch degree of confidence all are not set to 0.2.From the fault-tolerance angle, the degree of confidence that transition are set is 0.95, and the transition threshold value is 0.2, and transition storehouse institute degree of confidence is 1, and transition transition threshold value is 0, and transition transition time delay is 0.
3. time delay constraint
Definition main protection, nearly back-up protection, back-up protection far away, and the time delay of corresponding isolating switch constraint is respectively Δ T mr∈ [10,40], Δ T pr∈ [300,500], Δ T sr∈ [600,1100], Δ T mc∈ [40,60], Δ T pc∈ [20,40], Δ T sc∈ [40,100], unit is ms.
The fault diagnosis model developed has been considered between element fault and protection action, has been protected the time delay between action and circuit breaker trip to retrain; and by model can find out isolating switch and the protection storehouse be all by transition directly arrive bus and line library, like this isolating switch storehouse Δ T minwith Δ T maxtime delay constraint constantly should occur with respect to fault for what forward inference was tried to achieve in assignment.
Three, the adaptability of fault diagnosis model and fault-tolerance:
1), the adaptability in the network topology change situation.For example, when the existing fault diagnosis model based on Petri net changes (have newly-increased circuit or existing line out of service) in NETWORK STRUCTURE PRESERVING POWER SYSTEM, often need to rebuild the Petri pessimistic concurrency control, workload is large.And the Weighted Fuzzy temporal Petri nets model in the present invention; setting up bus main protection, back-up protection; and, after the submodel of power line main protection, nearly back-up protection, back-up protection far away, then carry out simple the fusion and form integrated diagnosis model, simple for structurely understand.When having newly-increased bus or existing bus out of service, for the model built, only need newly-increased or delete the model of associated bus, because model structure and the reasoning process of each bus is similarly, so just the energy rapid pin be made correction to network topology change.For circuit model; situation and bus model are similar; but consider that circuit may also be responsible for the standby of associated bus or circuit/far back-up protection simultaneously; when having newly-increased circuit or existing line out of service; not only need to increase or delete the model of this circuit; also need standby to respective bus bars and circuit/back-up protection far away to be adjusted, increase or delete corresponding protection submodel.The fault diagnosis model developed can adapt to network topology change fast, and counting yield is high.
2), the fault-tolerance in the incorrect situation of warning information.The fault diagnosis model of the temporal based Fuzzy Petri Net of setting up in the embodiment of the present invention is by analyzing the degree of confidence of candidate's fault element; can take into account protection and the malfunction of isolating switch/tripping situation; and warning information is not right-on situation, there is fault-tolerance preferably.In addition, according to configuration relation and the action logic of element, protection, isolating switch, adopt backward reasoning can differentiate the malfunction of protection and isolating switch/tripping situation.
Below; for having protection and likely malfunction or tripping of isolating switch in fault diagnosis; in the alarm upload procedure, also may occur uploading not in time, the uncertain factor such as distortion or loss, the present invention be take the IEEE39 node system shown in Fig. 5 and is verified the fault diagnosis model developed as example.Suppose to receive the warning information with sequential: protection R b14m(20ms), R (4)-14s(750ms) and R 14-(15) s(371ms) action; Isolating switch CB (14)-15(73ms), CB (14)-13(81ms) and CB (4)-14(87ms) tripping.
1), fault zone identification.Warning information based on received, determine fault zone by BFS (Breadth First Search), and possible fault element is B14 and L4-14.
2), to each element modeling in fault zone.Set up bus B14 as shown in Figure 3, Figure 4 and the Weighted Fuzzy temporal Petri nets fault diagnosis model of circuit L4-14.
3), the rational analysis based on matrix operation.According to the rational analysis of Weighted Fuzzy temporal Petri nets, construct corresponding matrix and carry out the reasoning computing:
1. adopt matrix operation to carry out rational analysis to the Weighted Fuzzy temporal Petri nets model of bus B14:
At first according to the time delay constraint, the alarm of receiving is screened, storehouse institute path matrix, minimum cumulative delay constrained vector and cumulative maximum time delay constrained vector that storehouse institute path is corresponding are respectively:
(annotate: the R of storehouse institute b14mstorehouse institute path only need count t 1, t 3and t 5in one get final product because three transition are of equal value.)
ΣΔT min=(A ccg(ΔT min) T) T=[640,640,640,50,50,50,10,600,600,600,0,0,0,0]
ΣΔT max=(A ccg(ΔT max) T) T=[1200,1200,1200,100,100,100,40,1100,1100,1100,0,0,0,0]
Afterwards, the consistance of check sequential.By the protection in warning information and isolating switch delayed data set Δ T mesminwith Δ T mesmaxwith the time delay constraint, compare, the protection and the isolating switch that are met constraint are judged vector (f i=0 means that sequential is inconsistent):
F = [ f B 14 m , f R ( 4 ) - 14 s , f R 14 - ( 15 ) s , f CB ( 14 ) - 15 , f CB ( 14 ) - 13 , f CB ( 4 ) - 14 ]
Figure BDA00003750885700103
Like this, the effective warning information after time delay constraint screening is: protection R b14m(20ms) and R (4)-14s(750ms) action; Isolating switch CB (14)-15(73ms), CB (14)-13(81ms) and CB (4)-14(87ms) action tripping operation.
The storehouse of Petri pessimistic concurrency control collects:
p={CB (4)-14,CB (13)-14,CB 14-(15),CB 4-(14),CB (14)-13,CB (14)-15,R B14m,R (4)-14s,R (13)-14s,R 14-(15)s,L4-14,L13-14,L14-15,B14}
The transition collection is:
T={t 1,t 2,t 3,t 4,t 5,t 6,t 7}
The input matrix of transition:
Figure BDA00003750885700111
The output matrix of transition:
Figure BDA00003750885700112
Transition degree of confidence vector:
U=[0.95,0.95,0.95,0.95,0.95,0.95,1]
Transition activation threshold vector:
T h=[0.2,0.2,0.2,0.2,0.2,0.2,0]
Input arc weight vector:
W=diag(0.50,0.50,0.50,0.40,0.40,0.40,0.60,0.50,0.50,0.50,0.33,0.33,0.33,0)
Given storehouse institute original state is:
M 0=[0.65,0.20,0.20,0.20,0.85,0.85,0.90,0.70,0.20,0.20,0.00,0.00,0.00,0.00]
Adopting the matrix reasoning process can obtain the degree of confidence that bus B14 breaks down is 0.771.
2. adopt matrix operation to carry out rational analysis to the Weighted Fuzzy temporal Petri nets model of circuit L4-14:
At first according to the time delay constraint, the alarm of receiving is screened, storehouse institute path matrix, minimum cumulative delay constrained vector and cumulative maximum time delay constrained vector that storehouse institute path is corresponding are respectively:
Figure BDA00003750885700121
ΣΔT min=(A ccg(ΔT min) T) T=[640,640,640,640,320,320,50,50,10,10,300,300,600,600,600,600,0,0,0,0,0]ΣΔT max=(A ccg(ΔT max) T) T=[1200,1200,1200,1200,540,540,100,100,40,40,500,500,1100,1100,1100,1100,0,0,0,0,0]
Afterwards, the consistance of check sequential.By the protection in warning information and isolating switch delayed data set Δ T mesminwith Δ T mesmaxwith the time delay constraint, compare, the protection and the isolating switch that are met constraint are judged vector (f i=0 means that sequential is inconsistent):
F = [ f B 14 m , f R ( 4 ) - 14 s , f R 14 - ( 15 ) s , f CB ( 14 ) - 15 , f CB ( 14 ) - 13 , f CB ( 4 ) - 14 ]
Figure BDA00003750885700123
Like this, the effective warning information after time delay constraint screening is: protection R b14m(20ms) and R (4)-14s(750ms) action; Isolating switch CB (14)-15(73ms), CB (14)-13(81ms) and CB (4)-14(87ms) action tripping operation.
The storehouse of Petri pessimistic concurrency control collects:
p={CB 4-(5),CB (3)-4,CB 14-(15),CB (13)-14,CB′ (4)-14,CB′ 4-(14),CB (4)-14,CB 4-(14),R (4)-14m,R 4-(14)m,
R (4)-14p,R 4-(14)p,R 4-(5)s,R (3)-4s,R 14-(15)s,R (13)-14s,p 1,p 2,L(4)-14,L4-(14),L4-14}
The transition collection is:
T={t 1,t 2,t 3,t 4,t 5,t 6,t 7,t 8,t 9,t 10,t 11}
The input matrix of transition:
The output matrix of transition:
Figure BDA00003750885700132
The time delay constrained vector:
ΔT min=[640,640,640,640,320,320,50,50,10,10,300,300,600,600,600,600,0,0,0,0,0]
Δ T max=[1200,1200,1200,1200,540,540,100,100,40,40,500,500,1100,1100,1100,1100,0,0,0,0,0] transition degree of confidence vector:
U=[0.95,0.95,0.95,0.95,0.95,0.95,0.95,0.95,1,1,1]
Transition activation threshold vector:
T h=[0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0,0,0]
Input arc weight vector:
W=diag(0.50,0.50,0.50,0.50,0.45,0.45,0.40,0.40,0.60,0.60,0.55,0.55,0.50,0.50,0.50,0.50,1,1,0.50,0.50,0
Adopting the matrix reasoning process can obtain the degree of confidence that circuit L4-14 breaks down is 0.219.
4), protection and isolating switch action evaluation
After above-mentioned steps is judged fault element, carry out now backward reasoning.Known after bus B14 breaks down, main protection R b14maction, trigger isolating switch CB (14)-15; , CB (14)-13and CB 4-(14)action tripping operation, and isolating switch CB wherein 4-(14)tripping, therefore corresponding back-up protection R (4)-14saction tripping isolating switch CB (4)-14thereby, isolated fault.Like this; Evaluation result is isolating switch CB 4-(14)tripping.
For the fault diagnosis model that checking develops better, in the embodiment of the present invention, also the various faults scene is tested.Table 3 has been listed the diagnostic result to the partial fault scene.Numerical results shows, method proposed by the invention can be processed protection and isolating switch has incorrect operation and the wrong complex fault situation of warning information.
Table 3 numerical testing
Figure BDA00003750885700141
Figure BDA00003750885700151
Identical with above-mentioned a kind of power system failure diagnostic method, the present invention also provides a kind of power system fault diagnosis, as shown in Figure 6, comprising:
Warning information receiver module 101, for receiving after electric power system fault occurs the warning information produced, described warning information comprises time sequence information;
Fault zone determination module 102, for determining fault zone according to described warning information;
Model building module 103, to each fault element modeling, set up the Weighted Fuzzy temporal Petri nets fault diagnosis model of fault element for the logical relation between each fault element according to described fault zone and corresponding protection, isolating switch action;
Reasoning computing module 104, for the rational analysis according to described Weighted Fuzzy temporal Petri nets fault diagnosis model, construct corresponding matrix and carry out the reasoning computing, diagnoses the element that is out of order.
As an embodiment preferably, described power system fault diagnosis can also comprise:
The action evaluation module, for after element is out of order in diagnosis, carry out backward reasoning to the fault element of diagnosing out, draws malfunction and the tripping situation of protection and isolating switch.
As an embodiment preferably, described fault zone determination module can adopt breadth-first search to determine described fault zone.
Other technical characterictic of above-mentioned a kind of power system fault diagnosis is identical with a kind of power system failure diagnostic method of the present invention, and it will not go into details herein.
By above scheme, can find out; a kind of power system failure diagnostic method of the present invention and system; the time delay binding feature of existence between protection and isolating switch action and the possibility of protection and isolating switch malfunction and tripping have been considered; propose a kind of electric system Weighted Fuzzy temporal Petri nets fault diagnosis model that can take into account this time delay constraint on the basis of existing fault diagnosis model, and carried out fault diagnosis according to this Weighted Fuzzy temporal Petri nets fault diagnosis model.Because fault diagnosis model of the present invention can be processed time delay restricted problem, therefore effectively improved fault-tolerance and the accuracy of fault diagnosis in the electric system.
The above embodiment has only expressed several embodiment of the present invention, and it describes comparatively concrete and detailed, but can not therefore be interpreted as the restriction to the scope of the claims of the present invention.It should be pointed out that for the person of ordinary skill of the art, without departing from the inventive concept of the premise, can also make some distortion and improvement, these all belong to protection scope of the present invention.Therefore, the protection domain of patent of the present invention should be as the criterion with claims.

Claims (6)

1. a power system failure diagnostic method, is characterized in that, comprises the following steps:
Receive after electric power system fault occurs the warning information produced, described warning information comprises time sequence information;
Determine fault zone according to described warning information;
To each fault element modeling, set up the Weighted Fuzzy temporal Petri nets fault diagnosis model of fault element according to the logical relation between each fault element in described fault zone and corresponding protection, isolating switch action;
According to the rational analysis of described Weighted Fuzzy temporal Petri nets fault diagnosis model, construct corresponding matrix and carry out the reasoning computing, diagnose the element that is out of order.
2. power system failure diagnostic method according to claim 1, is characterized in that, after element is out of order in diagnosis, also comprises step:
The fault element of diagnosing out is carried out to backward reasoning, draw malfunction and the tripping situation of protection and isolating switch.
3. power system failure diagnostic method according to claim 1 and 2, is characterized in that, the process of described definite fault zone comprises: adopt breadth-first search to determine described fault zone.
4. a power system fault diagnosis, is characterized in that, comprising:
The warning information receiver module, for receiving after electric power system fault occurs the warning information produced, described warning information comprises time sequence information;
The fault zone determination module, for determining fault zone according to described warning information;
Model building module, to each fault element modeling, set up the Weighted Fuzzy temporal Petri nets fault diagnosis model of fault element for the logical relation between each fault element according to described fault zone and corresponding protection, isolating switch action;
The reasoning computing module, for the rational analysis according to described Weighted Fuzzy temporal Petri nets fault diagnosis model, construct corresponding matrix and carry out the reasoning computing, diagnoses the element that is out of order.
5. power system fault diagnosis according to claim 4, is characterized in that, also comprises:
The action evaluation module, for after element is out of order in diagnosis, carry out backward reasoning to the fault element of diagnosing out, draws malfunction and the tripping situation of protection and isolating switch.
6. according to the described power system fault diagnosis of claim 4 or 5, it is characterized in that, described fault zone determination module adopts breadth-first search to determine described fault zone.
CN201310388672.XA 2013-08-30 2013-08-30 Fault diagnosis method of electric power system and system Expired - Fee Related CN103487723B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310388672.XA CN103487723B (en) 2013-08-30 2013-08-30 Fault diagnosis method of electric power system and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310388672.XA CN103487723B (en) 2013-08-30 2013-08-30 Fault diagnosis method of electric power system and system

Publications (2)

Publication Number Publication Date
CN103487723A true CN103487723A (en) 2014-01-01
CN103487723B CN103487723B (en) 2017-09-15

Family

ID=49828098

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310388672.XA Expired - Fee Related CN103487723B (en) 2013-08-30 2013-08-30 Fault diagnosis method of electric power system and system

Country Status (1)

Country Link
CN (1) CN103487723B (en)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103901320A (en) * 2014-03-11 2014-07-02 国家电网公司 Method for diagnosing power system fault considering multi-source data
CN104020396A (en) * 2014-06-12 2014-09-03 国家电网公司 Grid fault diagnosis method based on data mining
CN104065510A (en) * 2014-06-27 2014-09-24 北京思特奇信息技术股份有限公司 PetriNet-based system operation and maintenance monitoring method and PetriNet-based system operation and maintenance monitoring system
CN104297637A (en) * 2014-10-31 2015-01-21 国家电网公司 Power system fault diagnosis method comprehensively using electricity amount and timing sequence information
CN105005644A (en) * 2015-06-29 2015-10-28 山东科技大学 Method for detecting fault of three-phase asynchronous motor
CN105548815A (en) * 2016-01-14 2016-05-04 山东科技大学 Petri net power grid fault detection method based maximum likelihood decoding
CN106353690A (en) * 2016-09-20 2017-01-25 上海海事大学 Method for diagnosing lithium battery faults by Petri net
CN106771853A (en) * 2016-11-29 2017-05-31 华中科技大学 A kind of intelligent power network method for diagnosing faults based on event origin
CN106841910A (en) * 2016-12-20 2017-06-13 国网辽宁省电力有限公司沈阳供电公司 Imitative electromagnetism algorithm is melted into the Fault Diagnosis Method for Distribution Networks of timing ambiguity Petri network
CN108957230A (en) * 2018-06-28 2018-12-07 武汉科技大学 A kind of Fault Diagnosis Method for Distribution Networks and system based on layering timing sequence adaptive fuzzy H net
CN109459659A (en) * 2018-11-14 2019-03-12 国网黑龙江省电力有限公司信息通信公司 It is a kind of for predicting the system and method for electric power system fault
CN109729073A (en) * 2018-11-30 2019-05-07 全球能源互联网研究院有限公司 Network Abnormal recognition methods and system in a kind of electric network information physical system
CN109884473A (en) * 2019-03-29 2019-06-14 国网山东省电力公司费县供电公司 A kind of electric power overhaul system and method
CN110174465A (en) * 2019-05-27 2019-08-27 武汉科技大学 A kind of small internal leakage device of on-line monitoring hydraulic cylinder
CN110188421A (en) * 2019-05-16 2019-08-30 南京理工大学 Power grid cascading fault modeling emulation mode based on colored stochastic Petri net
CN110488156A (en) * 2019-08-30 2019-11-22 西南交通大学 A kind of high-speed railway distribution substation fault component diagnosis method
CN110927509A (en) * 2019-10-14 2020-03-27 山东科技大学 Fault diagnosis method for power distribution network
CN111337792A (en) * 2020-04-14 2020-06-26 上海海事大学 Power system fault diagnosis method based on improved Petri network
CN112180195A (en) * 2020-09-30 2021-01-05 国网福建省电力有限公司厦门供电公司 Method for acquiring all-state information and diagnosing faults of garden power supply and utilization equipment
CN113295969A (en) * 2021-06-30 2021-08-24 广东电网有限责任公司 Method, system, terminal and storage medium for rapidly positioning power distribution network fault
CN114252727A (en) * 2021-11-12 2022-03-29 国网浙江省电力有限公司杭州供电公司 Power grid fault rapid diagnosis method based on artificial intelligence technology

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6681374B1 (en) * 1999-06-09 2004-01-20 Lucent Technologies Inc. Hit-or-jump method and system for embedded testing
CN101726688A (en) * 2009-09-17 2010-06-09 浙江大学 Method for diagnosing multi-data-source information fusion-based power system fault
US7944843B2 (en) * 2004-02-13 2011-05-17 Alcatel-Lucent Usa Inc. Path based network management method and apparatus for data communication networks
CN102663500A (en) * 2012-03-27 2012-09-12 许继集团有限公司 Fuzzy Petri net transformer station fault diagnosis method based on time membership analysis
CN102680817A (en) * 2012-04-28 2012-09-19 辽宁省电力有限公司朝阳供电公司 Transformer fault diagnosis method based on fuzzy Petri net
CN103001328A (en) * 2012-11-19 2013-03-27 山东大学 Fault diagnosis and assessment method of intelligent substation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6681374B1 (en) * 1999-06-09 2004-01-20 Lucent Technologies Inc. Hit-or-jump method and system for embedded testing
US7944843B2 (en) * 2004-02-13 2011-05-17 Alcatel-Lucent Usa Inc. Path based network management method and apparatus for data communication networks
CN101726688A (en) * 2009-09-17 2010-06-09 浙江大学 Method for diagnosing multi-data-source information fusion-based power system fault
CN102663500A (en) * 2012-03-27 2012-09-12 许继集团有限公司 Fuzzy Petri net transformer station fault diagnosis method based on time membership analysis
CN102680817A (en) * 2012-04-28 2012-09-19 辽宁省电力有限公司朝阳供电公司 Transformer fault diagnosis method based on fuzzy Petri net
CN103001328A (en) * 2012-11-19 2013-03-27 山东大学 Fault diagnosis and assessment method of intelligent substation

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
文福拴等: "电力***故障诊断与不可观测保护状态识别的解析模型", 《华南理工大学学报(自然科学版)》 *
杨健维等: "基于方向性加权模糊Petri网的电网故障诊断方法", 《中国电机工程学报》 *
杨健维等: "基于时序模糊Petri网的电力***故障诊断", 《电力***自动化》 *

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103901320A (en) * 2014-03-11 2014-07-02 国家电网公司 Method for diagnosing power system fault considering multi-source data
CN104020396A (en) * 2014-06-12 2014-09-03 国家电网公司 Grid fault diagnosis method based on data mining
CN104065510A (en) * 2014-06-27 2014-09-24 北京思特奇信息技术股份有限公司 PetriNet-based system operation and maintenance monitoring method and PetriNet-based system operation and maintenance monitoring system
CN104297637A (en) * 2014-10-31 2015-01-21 国家电网公司 Power system fault diagnosis method comprehensively using electricity amount and timing sequence information
CN105005644B (en) * 2015-06-29 2019-03-01 山东科技大学 A method of detection threephase asynchronous failure
CN105005644A (en) * 2015-06-29 2015-10-28 山东科技大学 Method for detecting fault of three-phase asynchronous motor
CN105548815A (en) * 2016-01-14 2016-05-04 山东科技大学 Petri net power grid fault detection method based maximum likelihood decoding
CN106353690A (en) * 2016-09-20 2017-01-25 上海海事大学 Method for diagnosing lithium battery faults by Petri net
CN106353690B (en) * 2016-09-20 2019-04-02 上海海事大学 Utilize the method for Petri network diagnosis lithium battery failure
CN106771853A (en) * 2016-11-29 2017-05-31 华中科技大学 A kind of intelligent power network method for diagnosing faults based on event origin
CN106841910B (en) * 2016-12-20 2019-05-24 国网辽宁省电力有限公司沈阳供电公司 Imitative electromagnetism algorithm is melted into the Fault Diagnosis Method for Distribution Networks of timing ambiguity Petri network
CN106841910A (en) * 2016-12-20 2017-06-13 国网辽宁省电力有限公司沈阳供电公司 Imitative electromagnetism algorithm is melted into the Fault Diagnosis Method for Distribution Networks of timing ambiguity Petri network
CN108957230A (en) * 2018-06-28 2018-12-07 武汉科技大学 A kind of Fault Diagnosis Method for Distribution Networks and system based on layering timing sequence adaptive fuzzy H net
CN108957230B (en) * 2018-06-28 2024-04-23 武汉科技大学 Power distribution network fault diagnosis method and system based on hierarchical timing sequence self-adaptive fuzzy H network
CN109459659A (en) * 2018-11-14 2019-03-12 国网黑龙江省电力有限公司信息通信公司 It is a kind of for predicting the system and method for electric power system fault
CN109729073A (en) * 2018-11-30 2019-05-07 全球能源互联网研究院有限公司 Network Abnormal recognition methods and system in a kind of electric network information physical system
CN109729073B (en) * 2018-11-30 2023-03-10 全球能源互联网研究院有限公司 Network anomaly identification method and system in power grid information physical system
CN109884473A (en) * 2019-03-29 2019-06-14 国网山东省电力公司费县供电公司 A kind of electric power overhaul system and method
CN110188421B (en) * 2019-05-16 2022-11-04 南京理工大学 Power grid cascading failure modeling simulation method based on colored stochastic Petri network
CN110188421A (en) * 2019-05-16 2019-08-30 南京理工大学 Power grid cascading fault modeling emulation mode based on colored stochastic Petri net
CN110174465A (en) * 2019-05-27 2019-08-27 武汉科技大学 A kind of small internal leakage device of on-line monitoring hydraulic cylinder
CN110488156A (en) * 2019-08-30 2019-11-22 西南交通大学 A kind of high-speed railway distribution substation fault component diagnosis method
CN110927509A (en) * 2019-10-14 2020-03-27 山东科技大学 Fault diagnosis method for power distribution network
CN111337792A (en) * 2020-04-14 2020-06-26 上海海事大学 Power system fault diagnosis method based on improved Petri network
CN112180195A (en) * 2020-09-30 2021-01-05 国网福建省电力有限公司厦门供电公司 Method for acquiring all-state information and diagnosing faults of garden power supply and utilization equipment
CN112180195B (en) * 2020-09-30 2023-12-12 国网福建省电力有限公司厦门供电公司 Method for collecting full-state information and diagnosing faults of electric equipment supplied by park
CN113295969B (en) * 2021-06-30 2022-05-27 广东电网有限责任公司 Method, system, terminal and storage medium for rapidly positioning power distribution network fault
CN113295969A (en) * 2021-06-30 2021-08-24 广东电网有限责任公司 Method, system, terminal and storage medium for rapidly positioning power distribution network fault
CN114252727A (en) * 2021-11-12 2022-03-29 国网浙江省电力有限公司杭州供电公司 Power grid fault rapid diagnosis method based on artificial intelligence technology

Also Published As

Publication number Publication date
CN103487723B (en) 2017-09-15

Similar Documents

Publication Publication Date Title
CN103487723A (en) Electric system fault diagnosis method and system
CN103308824B (en) Power system fault diagnostic method based on probability Petri net
CN103001328B (en) Fault diagnosis and assessment method of intelligent substation
CN106226658B (en) Power grid fault diagnosis method based on multi-data fusion
Zhao et al. Power system fault diagnosis based on history driven differential evolution and stochastic time domain simulation
CN106771848B (en) A kind of electric network failure diagnosis method based on the fusion of topological primitive information
CN107394773A (en) Consider the distribution information physical system reliability estimation method of troubleshooting overall process
CN103901320A (en) Method for diagnosing power system fault considering multi-source data
CN105183952B (en) A kind of power transmission network method for diagnosing faults based on separation time Fuzzy Petri Net
CN101251835A (en) Method for evaluating reliability of +/- 800kV current exchanging station main wire connection
CN104281899A (en) Novel fault diagnosis method based on information fusion
CN105678337B (en) Information fusion method in intelligent substation fault diagnosis
CN110045207A (en) A kind of complex fault diagnostic method based on power grid architecture and multisource data fusion
CN105512011A (en) Electronic device testability modeling evaluation method
CN112052573A (en) Distributed fault identification method based on finite time observer
CN103197168A (en) Fault diagnosis control method in power system based on event set causal chain
Daigle A qualitative event-based approach to fault diagnosis of hybrid systems
CN106771853A (en) A kind of intelligent power network method for diagnosing faults based on event origin
CN103311960A (en) Method for dividing coherent stability regions of strong-connection power grids
CN114252727A (en) Power grid fault rapid diagnosis method based on artificial intelligence technology
Naseem et al. Probability estimation for the fault detection and isolation of pmu-based transmission line system of smart grid
CN106226657B (en) Power failure area identification method for power grid fault diagnosis
Reshmila et al. Diagnosis of power system failures using observer based discrete event system
Ye et al. Diagnosability analysis of discrete event systems with autonomous components
Meynen et al. Robust fault detection and isolation for distributed and decentralized systems

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170915

Termination date: 20180830

CF01 Termination of patent right due to non-payment of annual fee