# A Fault Diagnosis Method of Power Systems Based on Gray System Theory

Authors: Huang Darong; Tang Jianping; Zhao Ling
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

To provide some decision-making suggestions for fault diagnosis in power systems, a new model for identifying fault component is constructed by using Gray theory. Firstly, the basic concepts of Gray theory are introduced and explained in detail. And then the recognition algorithm of the power supply interrupted districts and the assignment principle of fault state vectors are depicted according to the working principle of protective relays (PRs) and circuit breakers (CBs). Secondly, based on the concept of the Gray correlation degree, the fault information explanation degree model is constructed and the judging method of malfunction and rejection for PRs and CBs is established. Meanwhile, to achieve the goal of the fault diagnosis, the fault diagnosis procedure that determined which components malfunction is designed for power systems. Finally, some simple experiments have already verified that the proposed method and model are effective and reasonable and the trend of further research is analyzed and summarized.


---

## Body

## 1. Introduction

Fault diagnosis of power systems is a method, which uses the information collected from protective relays and circuit breaker to recognize fault component and malfunctioned or tripped PRs and CBs. It is generally known that the fault component recognition is the crucial problem in engineering application [1]. In recent year, many fault diagnosis methods of power systems are proposed, including expert system [2], artificial neural network [3], optimization technology [4], rough set theory [5, 6], Petri net [7], and Bayesian network [6]. In order to meet certain preconditions in the existing diagnosis methods, some assumptions are made for them. But in some cases these assumptions may be contrary to reality and may even cause error diagnosis results. For example, when the information of protective device is incomplete, the fault diagnosis’ conclusion may be incorrect. So the crucial difficulty of fault diagnosis is how to guarantee the corrective and effective diagnosis results for power systems with incomplete information [8]. In other words, if we want to get the corrective and effective diagnosis results, we have to know how to obtain the unknown information by utilizing the known information under the condition of incomplete information [9–11].Fortunately, the Gray system theory can simulate and determine the unknown information according to the known information of systems. In addition, using the Gray systems theory to construct the fault diagnosis model involves a smaller sample and little information. The process is easier to be operated. However, the anomalies contrary to qualitative analysis will not be produced in the process of Gray correlation analysis. Hence, under the condition of information incomplete and not explicit, the Gray system theory has its unique superiority; that is, it is an effective method in case of incomplete and uncertain information.Based on the above, application of the Gray theory to the fault diagnosis of power systems is proposed in this paper. The layout of the rest of the paper is arranged as follows. in Section2 the Gray system theory and Gray correlation analysis were introduced briefly. Section 3 described the fault diagnosis method and model of power systems based on Gray system theory and designed the diagnostic procedure for fault components. In Section 4 we have discussed the simulative results by experiments. Finally, Section 5 concluded this paper with inferences and directions for future work.


## 2. Basis Introduction of Gray System Theory and Gray Correlation Analysis

### 2.1. Basal Principle of Gray System Theory

Gray system theory, which is established by Chinese scholar Professor Deng Ju-Long in 1982, has been used to research uncertain problem with lack of data and information. In general speaking, some information is known and other pieces of information are unknown in uncertain data systems. The main idea of Gray system theory is to describe correctly the systems’ evolution law and to monitor effectively their running behavior by extracting valuable information from the known information. Over forty years later, the structure system of Gray system theory has basically taken shape. The theoretical basis includes Gray Matrix, Gray Algebra, and Gray Equation. And the Gray model (GM) may also implement the analysis, evaluation, prediction, and control decision of uncertain data systems by using spatial association rule and sequence generation method. For future understanding of the advantages of Gray system theory, six basic principles are introduced as follows.


#### 2.1.1. Differential Information Principle

There is the difference among differential information. In other words, the difference is the information. For example, if there are two differential objects or systems, everyone has unique information that is not similar to another piece of information. One of the most basic pieces of information which the human society perceives is that the world comes from the difference between matter and matter.


#### 2.1.2. Nonuniqueness of Solutions

The solution is nonunique under the condition of information being incomplete and vague. The uncertainty of systems leads to the existence of uncertain information and then causes nonuniqueness of solutions.


#### 2.1.3. Smallest Information Principle

The basic idea of Gray system theory is to utilize the smallest information achieving from known knowledge data to accomplish a given task. The smallest information getting by researching the uncertain problem with the small samples and poor information is the fundamental basis of making a distinction between Gray area and no-Gray.


#### 2.1.4. Cognitive Foundation Principle

The Foundation of Cognition Is Information. The accurate and complete cognition is determined according to definite and precise knowledge; the uncertain and incomplete knowledge may also lead to vague cognition. Correspondingly, if there is no information of systems, the cognition of systems is not also completed. Thus the cognition should be studied based on information.


#### 2.1.5. Innovation Priority Principle

The function of new information is more important than old information for cognitive behaviors, because the new information directly reflects the current states of the system and mainly influences the future trend of the development.


#### 2.1.6. Gray Indestructibility Principle

Notably, the incomplete and uncertain information is very much strongly entrenched in real systems. As new information is continuously generated in real engineering, the cognition conclusion is improved gradually and the level of cognition will tend to rationality and correctness. Accordingly, the Gray systems do not disappear.Based on the above principles, the theoretical model implementing a given task may be constructed according to small sample and poor information. However, the running states and trend of systems are determined by many factors in real application. So we need to discern the primary factors and lesser factors. And fortunately, the Gray correlation analysis, which is the important part of Gray system theory, can judge the connection in accordance with the approximation degree of the two-dimensional curve in time domain and frequencies domain. The higher the similarity is, the greater the correlation is. For further analysis, the Gray relational axioms are introduced firstly.


### 2.2. Gray Correlation Axiom

Suppose that the behavior sequence of the system isX
0
=
(
x
0
(
1
)
,
x
0
(
2
)
,
…
,
x
0
(
n
)
); thus the corresponding factor series is as follows:(1)
X
1
=
x
1
1
,
x
1
2
,
…
,
x
1
n
,
W
i
⋮
X
i
=
x
i
1
,
x
i
2
,
…
,
x
i
n
,
W
i
⋮
X
m
=
x
m
1
,
x
m
2
,
…
,
x
m
n
.For a given real r
(
x
0
(
k
)
,
x
i
(
k
)
), if r
(
X
0
,
X
i
)
=
(
1
/
n
)
∑
k
=
1
n
r
(
x
0
(
k
)
,
x
i
(
k
)
) satisfies(1)
normalization:0
<
r
(
X
0
,
X
i
)
≤
1, and X
0
=
X
i
⇒
r
(
X
0
,
X
i
)
=
1;
(2)
integrity: letX
i
,
X
j
∈
X
=
{
X
s
∣
s
=
0,1
,
…
,
m
,
m
≥
2
}; if i
≠
j, then r
X
i
,
X
j
≠
r
X
j
,
X
i
.
(3)
even symmetry: ifX
=
{
X
i
,
X
j
}, X
=
{
X
i
,
X
j
}
⇔
r
(
X
i
,
X
j
)
=
r
(
X
j
,
X
i
);
(4)
accessibility: the value of the Euclid distancex
0
(
k
)
-
x
i
(
k
) is inversely proportional to the value of the given real r
(
x
0
(
k
)
,
x
i
(
k
)
).Thus r
(
X
0
,
X
i
)
=
(
1
/
n
)
∑
k
=
1
n
r
(
x
0
(
k
)
,
x
i
(
k
)
) is called as Gray relational degree, where r
(
X
0
,
X
i
) represents the correlation coefficient of each pair of variables X
i and X
j. Four conditions (1–4) are also regarded as Gray relational four axioms.Notice that four conditions just do positively mean these things. Normalization illustrates that there is a correlation between two arbitrary behavior sequences. Integrity describes that the Gray relational degree is influenced by the external environment. If the outside environment is changed, the relational degree is also varied. So the symmetry principle is not necessarily true. Meanwhile, even symmetry represents that the symmetry principle is true while the set of factors contains just two factors. Accessibility may constrain the relational quantization.


### 2.3. Gray Correlation Analysis

To get the computing formula of Gray related degree, the distance measure between vectorsx
0
(
k
) and x
i
(
k
) is defined as follows:(2)
Δ
0
i
k
=
x
0
k
-
x
i
k
.And suppose that(3)
Δ
m
a
x
=
max
i
⁡
max
k
⁡
Δ
0
i
k
,
Δ
m
i
n
=
min
i
⁡
min
k
⁡
Δ
0
i
k
.According to the formulas (2)-(3), the correlation coefficient between vectors x
0
(
k
) and vector x
i
(
k
) is defined as follows:(4)
r
x
0
k
,
x
i
k
=
Δ
m
i
n
+
ξ
Δ
m
a
x
Δ
m
i
n
+
ξ
Δ
m
a
x
;that is,(5)
r
x
0
k
,
x
i
k
=
min
i
⁡
min
k
⁡
x
0
k
-
x
i
k
+
ξ
max
i
⁡
max
k
⁡
x
0
k
-
x
i
k
x
0
k
-
x
i
k
+
ξ
max
i
⁡
max
k
⁡
x
0
k
-
x
i
k
,where ξ is called the resolution coefficients and the values of ξ are usually restricted to a certain range (
0,1
).Notice that the discriminatory power varies depending on the different correlation coefficients: the smaller theξ is, the higher the differences between two correlation coefficients are and the stronger the discriminatory power is.Letr
x
0
(
k
)
,
x
i
(
k
)
=
r
0
i
(
k
), and then we define(6)
r
X
0
,
X
i
=
1
n
∑
k
=
1
n
r
x
0
k
,
x
i
k
=
1
n
∑
k
=
1
n
r
0
i
k
,and r
X
0
,
X
i is considered as the Gray correlation degree between reference sequence X
0 and compare sequence X
i. Obviously, r
X
0
,
X
i satisfy Gray relational four axioms (1–4).


### 2.4. Computing Algorithm of Gray Correlation Degree

By the definition of Gray correlation degree, the computational steps of Gray correlation are made as follows.(1) Collect the evaluation data on the evaluation index system; then the sequences of data may be stated in matrix form as follows:(7)
X
0
,
X
1
,
…
,
X
m
=
x
0
1
x
1
1
⋯
x
m
1
x
0
2
x
1
2
⋯
x
m
2
⋮
⋮
⋯
⋮
x
0
n
x
1
n
⋯
x
m
n
,where n represents the number of indexes. And X
i
=
(
x
i
(
1
)
,
x
i
(
2
)
,
…
,
x
i
(
n
)
), i
=
1,2
,
…
,
m.(2) Apply dimensionless method to the original data sequences; let the dimensionless model be(8)
x
i
′
k
=
x
i
k
1
/
n
∑
k
=
1
n
x
i
k
,
i
=
0,1
,
…
,
m
;
k
=
1,2
,
…
,
n
.Thus the new data sequences processed by dimensionless model may be rewritten as(9)
X
0
,
X
1
,
…
,
X
m
=
x
0
′
1
x
1
′
1
⋯
x
m
′
1
x
0
′
2
x
1
′
2
⋯
x
m
′
2
⋮
⋮
⋯
⋮
x
0
′
n
x
1
′
n
⋯
x
m
′
n
.(3) Define the reference sequenceX
0
′. The reference sequence consists of the most optimal value or the worst value of every index. That is, X
0
′
=
(
x
0
′
(
1
)
,
x
0
′
(
2
)
,
…
,
x
0
′
(
n
)
). Accordingly, the rest of the data is as compare sequence.(4) Compute the distance measure between the corresponding elements of the reference sequenceX
0
′
(
k
) and compare sequence X
i
′
(
k
); that is, Δ
0
i
k
=
∥
x
0
′
(
k
)
-
x
i
′
(
k
)
∥.(5) CalculateΔ
m
a
x
,
Δ
m
i
n using the formula (3); that is,(10)
Δ
m
a
x
=
max
i
⁡
max
k
⁡
x
0
′
k
-
x
i
′
k
,
Δ
m
i
n
=
min
i
⁡
min
k
⁡
x
0
′
k
-
x
i
′
k
.(6) Compute the correlation coefficient between vectorx
0
′
(
k
) and vector x
i
′
(
k
) by formulas (4) and (5); that is, (11)
r
x
0
′
k
,
x
i
′
k
=
Δ
m
i
n
+
ξ
Δ
m
a
x
Δ
0
i
k
+
ξ
Δ
m
a
x
.Notice that ξ is the resolution coefficients and the values of ξ are usually restricted to a certain range (
0,1
). To keep things simple, put ξ
=
0.5 in this paper.(7) Compute the correlation degree between reference sequenceX
0 and compare sequence X
i by formula (6), and then obtain the evaluation conclusion by comparing the size of correlation degree.


## 2.1. Basal Principle of Gray System Theory

Gray system theory, which is established by Chinese scholar Professor Deng Ju-Long in 1982, has been used to research uncertain problem with lack of data and information. In general speaking, some information is known and other pieces of information are unknown in uncertain data systems. The main idea of Gray system theory is to describe correctly the systems’ evolution law and to monitor effectively their running behavior by extracting valuable information from the known information. Over forty years later, the structure system of Gray system theory has basically taken shape. The theoretical basis includes Gray Matrix, Gray Algebra, and Gray Equation. And the Gray model (GM) may also implement the analysis, evaluation, prediction, and control decision of uncertain data systems by using spatial association rule and sequence generation method. For future understanding of the advantages of Gray system theory, six basic principles are introduced as follows.


### 2.1.1. Differential Information Principle

There is the difference among differential information. In other words, the difference is the information. For example, if there are two differential objects or systems, everyone has unique information that is not similar to another piece of information. One of the most basic pieces of information which the human society perceives is that the world comes from the difference between matter and matter.


### 2.1.2. Nonuniqueness of Solutions

The solution is nonunique under the condition of information being incomplete and vague. The uncertainty of systems leads to the existence of uncertain information and then causes nonuniqueness of solutions.


### 2.1.3. Smallest Information Principle

The basic idea of Gray system theory is to utilize the smallest information achieving from known knowledge data to accomplish a given task. The smallest information getting by researching the uncertain problem with the small samples and poor information is the fundamental basis of making a distinction between Gray area and no-Gray.


### 2.1.4. Cognitive Foundation Principle

The Foundation of Cognition Is Information. The accurate and complete cognition is determined according to definite and precise knowledge; the uncertain and incomplete knowledge may also lead to vague cognition. Correspondingly, if there is no information of systems, the cognition of systems is not also completed. Thus the cognition should be studied based on information.


### 2.1.5. Innovation Priority Principle

The function of new information is more important than old information for cognitive behaviors, because the new information directly reflects the current states of the system and mainly influences the future trend of the development.


### 2.1.6. Gray Indestructibility Principle

Notably, the incomplete and uncertain information is very much strongly entrenched in real systems. As new information is continuously generated in real engineering, the cognition conclusion is improved gradually and the level of cognition will tend to rationality and correctness. Accordingly, the Gray systems do not disappear.Based on the above principles, the theoretical model implementing a given task may be constructed according to small sample and poor information. However, the running states and trend of systems are determined by many factors in real application. So we need to discern the primary factors and lesser factors. And fortunately, the Gray correlation analysis, which is the important part of Gray system theory, can judge the connection in accordance with the approximation degree of the two-dimensional curve in time domain and frequencies domain. The higher the similarity is, the greater the correlation is. For further analysis, the Gray relational axioms are introduced firstly.


## 2.1.1. Differential Information Principle

There is the difference among differential information. In other words, the difference is the information. For example, if there are two differential objects or systems, everyone has unique information that is not similar to another piece of information. One of the most basic pieces of information which the human society perceives is that the world comes from the difference between matter and matter.


## 2.1.2. Nonuniqueness of Solutions

The solution is nonunique under the condition of information being incomplete and vague. The uncertainty of systems leads to the existence of uncertain information and then causes nonuniqueness of solutions.


## 2.1.3. Smallest Information Principle

The basic idea of Gray system theory is to utilize the smallest information achieving from known knowledge data to accomplish a given task. The smallest information getting by researching the uncertain problem with the small samples and poor information is the fundamental basis of making a distinction between Gray area and no-Gray.


## 2.1.4. Cognitive Foundation Principle

The Foundation of Cognition Is Information. The accurate and complete cognition is determined according to definite and precise knowledge; the uncertain and incomplete knowledge may also lead to vague cognition. Correspondingly, if there is no information of systems, the cognition of systems is not also completed. Thus the cognition should be studied based on information.


## 2.1.5. Innovation Priority Principle

The function of new information is more important than old information for cognitive behaviors, because the new information directly reflects the current states of the system and mainly influences the future trend of the development.


## 2.1.6. Gray Indestructibility Principle

Notably, the incomplete and uncertain information is very much strongly entrenched in real systems. As new information is continuously generated in real engineering, the cognition conclusion is improved gradually and the level of cognition will tend to rationality and correctness. Accordingly, the Gray systems do not disappear.Based on the above principles, the theoretical model implementing a given task may be constructed according to small sample and poor information. However, the running states and trend of systems are determined by many factors in real application. So we need to discern the primary factors and lesser factors. And fortunately, the Gray correlation analysis, which is the important part of Gray system theory, can judge the connection in accordance with the approximation degree of the two-dimensional curve in time domain and frequencies domain. The higher the similarity is, the greater the correlation is. For further analysis, the Gray relational axioms are introduced firstly.


## 2.2. Gray Correlation Axiom

Suppose that the behavior sequence of the system isX
0
=
(
x
0
(
1
)
,
x
0
(
2
)
,
…
,
x
0
(
n
)
); thus the corresponding factor series is as follows:(1)
X
1
=
x
1
1
,
x
1
2
,
…
,
x
1
n
,
W
i
⋮
X
i
=
x
i
1
,
x
i
2
,
…
,
x
i
n
,
W
i
⋮
X
m
=
x
m
1
,
x
m
2
,
…
,
x
m
n
.For a given real r
(
x
0
(
k
)
,
x
i
(
k
)
), if r
(
X
0
,
X
i
)
=
(
1
/
n
)
∑
k
=
1
n
r
(
x
0
(
k
)
,
x
i
(
k
)
) satisfies(1)
normalization:0
<
r
(
X
0
,
X
i
)
≤
1, and X
0
=
X
i
⇒
r
(
X
0
,
X
i
)
=
1;
(2)
integrity: letX
i
,
X
j
∈
X
=
{
X
s
∣
s
=
0,1
,
…
,
m
,
m
≥
2
}; if i
≠
j, then r
X
i
,
X
j
≠
r
X
j
,
X
i
.
(3)
even symmetry: ifX
=
{
X
i
,
X
j
}, X
=
{
X
i
,
X
j
}
⇔
r
(
X
i
,
X
j
)
=
r
(
X
j
,
X
i
);
(4)
accessibility: the value of the Euclid distancex
0
(
k
)
-
x
i
(
k
) is inversely proportional to the value of the given real r
(
x
0
(
k
)
,
x
i
(
k
)
).Thus r
(
X
0
,
X
i
)
=
(
1
/
n
)
∑
k
=
1
n
r
(
x
0
(
k
)
,
x
i
(
k
)
) is called as Gray relational degree, where r
(
X
0
,
X
i
) represents the correlation coefficient of each pair of variables X
i and X
j. Four conditions (1–4) are also regarded as Gray relational four axioms.Notice that four conditions just do positively mean these things. Normalization illustrates that there is a correlation between two arbitrary behavior sequences. Integrity describes that the Gray relational degree is influenced by the external environment. If the outside environment is changed, the relational degree is also varied. So the symmetry principle is not necessarily true. Meanwhile, even symmetry represents that the symmetry principle is true while the set of factors contains just two factors. Accessibility may constrain the relational quantization.


## 2.3. Gray Correlation Analysis

To get the computing formula of Gray related degree, the distance measure between vectorsx
0
(
k
) and x
i
(
k
) is defined as follows:(2)
Δ
0
i
k
=
x
0
k
-
x
i
k
.And suppose that(3)
Δ
m
a
x
=
max
i
⁡
max
k
⁡
Δ
0
i
k
,
Δ
m
i
n
=
min
i
⁡
min
k
⁡
Δ
0
i
k
.According to the formulas (2)-(3), the correlation coefficient between vectors x
0
(
k
) and vector x
i
(
k
) is defined as follows:(4)
r
x
0
k
,
x
i
k
=
Δ
m
i
n
+
ξ
Δ
m
a
x
Δ
m
i
n
+
ξ
Δ
m
a
x
;that is,(5)
r
x
0
k
,
x
i
k
=
min
i
⁡
min
k
⁡
x
0
k
-
x
i
k
+
ξ
max
i
⁡
max
k
⁡
x
0
k
-
x
i
k
x
0
k
-
x
i
k
+
ξ
max
i
⁡
max
k
⁡
x
0
k
-
x
i
k
,where ξ is called the resolution coefficients and the values of ξ are usually restricted to a certain range (
0,1
).Notice that the discriminatory power varies depending on the different correlation coefficients: the smaller theξ is, the higher the differences between two correlation coefficients are and the stronger the discriminatory power is.Letr
x
0
(
k
)
,
x
i
(
k
)
=
r
0
i
(
k
), and then we define(6)
r
X
0
,
X
i
=
1
n
∑
k
=
1
n
r
x
0
k
,
x
i
k
=
1
n
∑
k
=
1
n
r
0
i
k
,and r
X
0
,
X
i is considered as the Gray correlation degree between reference sequence X
0 and compare sequence X
i. Obviously, r
X
0
,
X
i satisfy Gray relational four axioms (1–4).


## 2.4. Computing Algorithm of Gray Correlation Degree

By the definition of Gray correlation degree, the computational steps of Gray correlation are made as follows.(1) Collect the evaluation data on the evaluation index system; then the sequences of data may be stated in matrix form as follows:(7)
X
0
,
X
1
,
…
,
X
m
=
x
0
1
x
1
1
⋯
x
m
1
x
0
2
x
1
2
⋯
x
m
2
⋮
⋮
⋯
⋮
x
0
n
x
1
n
⋯
x
m
n
,where n represents the number of indexes. And X
i
=
(
x
i
(
1
)
,
x
i
(
2
)
,
…
,
x
i
(
n
)
), i
=
1,2
,
…
,
m.(2) Apply dimensionless method to the original data sequences; let the dimensionless model be(8)
x
i
′
k
=
x
i
k
1
/
n
∑
k
=
1
n
x
i
k
,
i
=
0,1
,
…
,
m
;
k
=
1,2
,
…
,
n
.Thus the new data sequences processed by dimensionless model may be rewritten as(9)
X
0
,
X
1
,
…
,
X
m
=
x
0
′
1
x
1
′
1
⋯
x
m
′
1
x
0
′
2
x
1
′
2
⋯
x
m
′
2
⋮
⋮
⋯
⋮
x
0
′
n
x
1
′
n
⋯
x
m
′
n
.(3) Define the reference sequenceX
0
′. The reference sequence consists of the most optimal value or the worst value of every index. That is, X
0
′
=
(
x
0
′
(
1
)
,
x
0
′
(
2
)
,
…
,
x
0
′
(
n
)
). Accordingly, the rest of the data is as compare sequence.(4) Compute the distance measure between the corresponding elements of the reference sequenceX
0
′
(
k
) and compare sequence X
i
′
(
k
); that is, Δ
0
i
k
=
∥
x
0
′
(
k
)
-
x
i
′
(
k
)
∥.(5) CalculateΔ
m
a
x
,
Δ
m
i
n using the formula (3); that is,(10)
Δ
m
a
x
=
max
i
⁡
max
k
⁡
x
0
′
k
-
x
i
′
k
,
Δ
m
i
n
=
min
i
⁡
min
k
⁡
x
0
′
k
-
x
i
′
k
.(6) Compute the correlation coefficient between vectorx
0
′
(
k
) and vector x
i
′
(
k
) by formulas (4) and (5); that is, (11)
r
x
0
′
k
,
x
i
′
k
=
Δ
m
i
n
+
ξ
Δ
m
a
x
Δ
0
i
k
+
ξ
Δ
m
a
x
.Notice that ξ is the resolution coefficients and the values of ξ are usually restricted to a certain range (
0,1
). To keep things simple, put ξ
=
0.5 in this paper.(7) Compute the correlation degree between reference sequenceX
0 and compare sequence X
i by formula (6), and then obtain the evaluation conclusion by comparing the size of correlation degree.


## 3. Model and Algorithm of Fault Diagnosis of Power Systems Based on Gray System Theory

### 3.1. Overall Description of Fault Diagnosis of Power Systems

Firstly, when some relevant components of power systems malfunction, the corresponding protective relays and circuit breakers will work actively and the fault components will be disconnected with the power supplier. So the power supply interrupted districts will be formed. And then analyzing and eliminating breakdown will be also made in the power supply interrupted districts. Secondly, in the blackout area, the working principle of protective relays and circuit breakers will be analyzed in detail, and the information of the protective relays and circuit breakers, which are in fault state, may be deduced. And based on what is mentioned above, we may establish the state vector of the fault modes and quarantined modes of the components and give the assignment principle of the value for each member of the vector. Finally, the correlation degree between the reference component and compare component will be calculated and sorted. The discriminant criterion is that the component with maximum correlation value is as the fault component. Moreover, the numbers of fault components increase constantly. When the information explanation degree of fault components reaches a specified threshold and the numbers of protective relays and circuit breakers which fail to operate or refused operation do not increase, the number of fault components is determined. According to the order of correlation degree, those components which have same account of maloperation and refused operation are regarded as the fault component. So the ultimate goal of the fault diagnosis is also accomplished.Obviously, some basic concepts should be explained such as the power supply interrupted districts, fault information explanation degree, and judging rule of maloperation and refused operation. Next, these related concepts and model will be explained step by step.


### 3.2. Quick Recognition for the Power Supply Interrupted Districts

The formation mechanism of the power supply interrupted districts is established on the difference of the topological structures of power systems before and after fault occurrence. When the fault happens, the protection action will be implemented to trip the relevant circuit breakers and the fault components can be isolated from the systems for avoiding the expansion of accident. On the basis of the real time information of the circuit breaker, the topologies of systems before and after fault occurrence are recognized by the real time topology analysis on power systems. And the fault components can automatically form some passive networks according to the difference of the information of the topological structures. So these passive networks are regarded as the power supply interrupted districts. Thus, the recognition of fault components may be limited to the blackout area.The specific steps [12] are shown as follows.(1)
Based on topology analysis for the normal power systems, the corresponding equivalent power or the generator are remarked as the active nodes.
(2)
The topology analysis is made for the faulty system again, and several subsystems are achieved at the same time according to the different information of the topological structures.
(3)
The nodes of each subsystem are searched for one by one, and every component connected with each node is active or is not judged. If the component is active, the subsystem is in the normal state, and the search ends; otherwise, if all components in a subsystem are searched for and no active components are found, the subsystem is regarded as the power supply interrupted district.


### 3.3. Construction for Fault State Vectors and Assignment Principle for Vector Elements

The state vectors are established by means of the protective relays and circuit breakers. If there aren circuit breakers and m protective relays, the form of state vectors can be expressed as(12)
F
k
=
R
,
C
=
r
1
,
r
2
,
…
,
r
m
,
c
1
,
c
2
,
…
,
c
m
,
k
=
1,2
,
…
,
n
,where R is the state value of protective relays and C is the state value of circuit breakers. The assignment principle of vector elements is as follows.(1)
If the protective relay run in active mode, thenr
i
=
1, i
=
1,2
,
…
,
m; otherwise r
i
=
0.
(2)
Analogously, if the circuit breaker is kept in off state, thenc
i
=
1, i
=
1,2
,
…
,
n; otherwise c
i
=
0.
(3)
When the state of the element possesses two modes simultaneously, the information of the protective relays conflicts with the information of the circuit breakers. That is to say, the value of state is both 1 and 0. In that case, the value ofR and C is given as 0.5; that is, r
i
=
0.5, i
=
1,2
,
…
,
m; c
i
=
0.5, i
=
1,2
,
…
,
n.


### 3.4. Modeling the Fault Information Explanation Degree

Fault information explanation degree (FIED) is the matching degree of quarantined state vector and fault state vector of multicomponent system. In the other words, it is used for explaining the behavior of protective relays and circuit breakers when the fault of multicomponent systems happens.LetF
A
,
F
B
,
…
,
F
X, whose sizes of dimension are all 1
×
(
m
+
n
), be the fault state vector, respectively, to single-fault components A
,
B
,
…
,
X. F
0 is quarantined state vector with 1
×
(
m
+
n
) dimension.Therefore, the fault information explanation degree (FIED) of single component can be expressed by the following formula:(13)
F
I
E
D
=
t
h
e
n
u
m
b
e
r
o
f
e
l
e
m
e
n
t
s
w
h
o
s
e
v
a
l
u
e
s
a
r
e
a
l
l
F
A
+
F
0
1
-
1
W
W
W
W
i
g
r
e
a
t
e
r
t
h
a
n
1
a
m
o
n
g
F
A
+
F
0
W
W
W
W
×
t
h
e
n
u
m
b
e
r
o
f
e
l
e
m
e
n
t
s
a
l
l
e
q
u
a
l
t
o
1
W
W
W
W
W
w
h
o
s
e
v
a
l
u
e
s
a
r
e
a
l
l
e
q
u
a
l
t
o
1
-
1
×
100
%
.If there are more than one faulty component in power system, FIED can be presented as follows:(14)
FIED
=
t
h
e
n
u
m
b
e
r
o
f
e
l
e
m
e
n
t
s
w
h
o
s
e
v
a
l
u
e
s
a
r
e
F
A
+
F
0
e
q
u
a
l
t
o
1
-
1
W
W
i
W
i
i
W
a
l
l
g
r
e
a
t
e
r
t
h
a
n
1
W
W
i
W
i
i
W
i
n
F
A
+
F
0
∨
F
B
+
F
0
∨
F
C
+
F
0
∨
⋯
W
W
W
i
i
W
×
t
h
e
n
u
m
b
e
r
o
f
e
l
e
m
e
n
t
s
W
W
W
W
i
i
W
w
h
o
s
e
v
a
l
u
e
s
a
r
e
a
l
l
e
q
u
a
l
t
o
1
-
1
×
100
%
,where the symbol “∨” denotes that the greater value is selected for the values of corresponding elements.To locate out the fault components in systems, we need to compare the size of FIED. Once FIED satisfiesF
I
E
D
=
100
%, the fault components may be determined by formulas (13) and (14).


### 3.5. Judging Method of Malfunction and Rejection for Protective Relays and Circuit Breakers

According to the assumption and analysis above, the judging rule is designed as follows.(1) If the value of element of the vector (F
A
-
F
0) is negative, then the corresponding protection of protective relays or action of circuit breakers is regarded as malfunction for single-component fault. Otherwise, the protection of protective relays or trip of circuit breakers is regarded as rejecting action.(2) If there is more than one fault component in the system, the judging rule is as follows. The protection of protective relays or action of circuit breakers is regarded as malfunction while the value of the vectorF
A
∨
F
B
∨
F
C
∨
⋯
-
F
0 is negative; otherwise the action is regarded as rejecting action.Notice that when failures happen, the main protection of PRs is first implemented to trip the circuit breakers. If the main protective relays do not work, then the backup protective relays begin to work to trip the circuit breakers. Therefore, the situation of the backup PRs deduced by (1) and (2) was the rejecting action; it may be that the action of PRs has already tripped the action of CBs. So we concern only the main protective relays and the circuit breakers and then give priority to main PRs when the malfunction and rejecting action are considered in this paper.


### 3.6. Fault Diagnosis Model of Single Component

The problem of recognition for single-fault component can be solved by seeking for the fault assumptions, which can best explain the alarming information and can be expressed by the following maximization problem:(15)
R
i
=
1
n
∑
k
=
1
n
min
i
⁡
min
k
⁡
F
0
k
-
F
i
k
+
ξ
max
i
⁡
max
k
⁡
F
0
k
-
F
i
k
-
1
W
W
W
W
i
W
+
ξ
max
i
⁡
max
k
⁡
F
0
k
-
F
i
k
W
W
W
W
W
×
F
0
k
-
F
i
k
max
k
⁡
F
0
k
-
F
i
k
W
W
W
W
W
i
W
+
ξ
max
i
⁡
max
k
⁡
F
0
k
-
F
i
k
-
1
,where i stands for the total number of the single-fault, and its value is determined by adding the number of components and the number of lines; n denotes the number of vector elements, and its value is the number of protective relays and circuit breakers; ξ is on behalf of the resolution coefficient; F
0
(
k
) indicates the quarantined faulty state vector; F
i
(
k
) means the fault state vector of single-line or single component and also is called the single-fault state vector; R
i represents the similarity degree between the quarantined fault state vector and single-fault state vector and also is called correlation degree.Notice that if the state vector of single-fault componentA
,
B
,
…
,
X is F
A
,
F
B
,
…
,
F
X, respectively, then the dimension of these state vectors is 1
∗
(
m
+
n
) and the dimension of the quarantined fault state vector F
0 is also 1
∗
(
m
+
n
). Meanwhile, we may calculate the correlation degree and sort the calculated results. Based on above discussion, the recognition rule of fault component is given as follows.Rule 1.
If the correlation degree between the reference state vectorF
0 and the compare state vector F
i
(
i
=
A
,
B
,
…
,
X
) is the largest, the corresponding part is fault component or fault line.Rule 2.
Let the threshold of fault diagnosis be the mean of correlation degree. If the correlation degree between the reference state vectorF
0 and new compare state vector F
N is more than the threshold, the part is regarded as the fault component or fault line.Therefore, the possibility that fault happens may be determined for every single component according to Rules1 and 2.


### 3.7. Diagnostic Procedure for Fault Components

In the actual project, there is more than one fault component when faults occur. But the diagnostic procedure may be designed on the basis of the fault diagnosis model of single component. So the specific steps of fault diagnosis are shown as follows.


#### 3.7.1. Data Preprocessing

The data, which contains the information of protective relays and circuit breakers and the topology of the power system, is written in the system.


#### 3.7.2. Recognition for the Power Supply Interrupted Districts

The blackout area of power systems is recognized, and the tripping principle of protective relays and circuit breaker in blackout area is listed.


#### 3.7.3. Recognition of State Information

Under single-fault state, the state information of protective relays and circuit breaker is obtained in accordance with the working principle of protective relays and circuit breakers.


#### 3.7.4. Construction for Fault State Vectors

The state vector of fault mode for each component is given by the assignment principle of vector elements, and then the quarantined fault state vector is loaded.


#### 3.7.5. Calculating and Sorting the Correlation Degree

The correlation coefficients are calculated, respectively, and then these results are sorted. If the correlation degree between quarantined state vector and reference state vector is the highest, the probability that component malfunctions is also the largest.


#### 3.7.6. Computing Fault Information Explanation Degree

FIED is calculated by formulas (13) and (14) and is sorted. If FIED of some fault components satisfies F
I
E
D
=
100
%, go to next step.


#### 3.7.7. Fault Diagnosis Analysis

According to Rules1 and 2, the fault component can be recognized and diagnosed. And then judge the malfunction and rejection for protective relays and circuit breakers.The specific flowchart for fault diagnosis is shown as in Figure1.Figure 1
The flowchart for fault diagnosis.


## 3.1. Overall Description of Fault Diagnosis of Power Systems

Firstly, when some relevant components of power systems malfunction, the corresponding protective relays and circuit breakers will work actively and the fault components will be disconnected with the power supplier. So the power supply interrupted districts will be formed. And then analyzing and eliminating breakdown will be also made in the power supply interrupted districts. Secondly, in the blackout area, the working principle of protective relays and circuit breakers will be analyzed in detail, and the information of the protective relays and circuit breakers, which are in fault state, may be deduced. And based on what is mentioned above, we may establish the state vector of the fault modes and quarantined modes of the components and give the assignment principle of the value for each member of the vector. Finally, the correlation degree between the reference component and compare component will be calculated and sorted. The discriminant criterion is that the component with maximum correlation value is as the fault component. Moreover, the numbers of fault components increase constantly. When the information explanation degree of fault components reaches a specified threshold and the numbers of protective relays and circuit breakers which fail to operate or refused operation do not increase, the number of fault components is determined. According to the order of correlation degree, those components which have same account of maloperation and refused operation are regarded as the fault component. So the ultimate goal of the fault diagnosis is also accomplished.Obviously, some basic concepts should be explained such as the power supply interrupted districts, fault information explanation degree, and judging rule of maloperation and refused operation. Next, these related concepts and model will be explained step by step.


## 3.2. Quick Recognition for the Power Supply Interrupted Districts

The formation mechanism of the power supply interrupted districts is established on the difference of the topological structures of power systems before and after fault occurrence. When the fault happens, the protection action will be implemented to trip the relevant circuit breakers and the fault components can be isolated from the systems for avoiding the expansion of accident. On the basis of the real time information of the circuit breaker, the topologies of systems before and after fault occurrence are recognized by the real time topology analysis on power systems. And the fault components can automatically form some passive networks according to the difference of the information of the topological structures. So these passive networks are regarded as the power supply interrupted districts. Thus, the recognition of fault components may be limited to the blackout area.The specific steps [12] are shown as follows.(1)
Based on topology analysis for the normal power systems, the corresponding equivalent power or the generator are remarked as the active nodes.
(2)
The topology analysis is made for the faulty system again, and several subsystems are achieved at the same time according to the different information of the topological structures.
(3)
The nodes of each subsystem are searched for one by one, and every component connected with each node is active or is not judged. If the component is active, the subsystem is in the normal state, and the search ends; otherwise, if all components in a subsystem are searched for and no active components are found, the subsystem is regarded as the power supply interrupted district.


## 3.3. Construction for Fault State Vectors and Assignment Principle for Vector Elements

The state vectors are established by means of the protective relays and circuit breakers. If there aren circuit breakers and m protective relays, the form of state vectors can be expressed as(12)
F
k
=
R
,
C
=
r
1
,
r
2
,
…
,
r
m
,
c
1
,
c
2
,
…
,
c
m
,
k
=
1,2
,
…
,
n
,where R is the state value of protective relays and C is the state value of circuit breakers. The assignment principle of vector elements is as follows.(1)
If the protective relay run in active mode, thenr
i
=
1, i
=
1,2
,
…
,
m; otherwise r
i
=
0.
(2)
Analogously, if the circuit breaker is kept in off state, thenc
i
=
1, i
=
1,2
,
…
,
n; otherwise c
i
=
0.
(3)
When the state of the element possesses two modes simultaneously, the information of the protective relays conflicts with the information of the circuit breakers. That is to say, the value of state is both 1 and 0. In that case, the value ofR and C is given as 0.5; that is, r
i
=
0.5, i
=
1,2
,
…
,
m; c
i
=
0.5, i
=
1,2
,
…
,
n.


## 3.4. Modeling the Fault Information Explanation Degree

Fault information explanation degree (FIED) is the matching degree of quarantined state vector and fault state vector of multicomponent system. In the other words, it is used for explaining the behavior of protective relays and circuit breakers when the fault of multicomponent systems happens.LetF
A
,
F
B
,
…
,
F
X, whose sizes of dimension are all 1
×
(
m
+
n
), be the fault state vector, respectively, to single-fault components A
,
B
,
…
,
X. F
0 is quarantined state vector with 1
×
(
m
+
n
) dimension.Therefore, the fault information explanation degree (FIED) of single component can be expressed by the following formula:(13)
F
I
E
D
=
t
h
e
n
u
m
b
e
r
o
f
e
l
e
m
e
n
t
s
w
h
o
s
e
v
a
l
u
e
s
a
r
e
a
l
l
F
A
+
F
0
1
-
1
W
W
W
W
i
g
r
e
a
t
e
r
t
h
a
n
1
a
m
o
n
g
F
A
+
F
0
W
W
W
W
×
t
h
e
n
u
m
b
e
r
o
f
e
l
e
m
e
n
t
s
a
l
l
e
q
u
a
l
t
o
1
W
W
W
W
W
w
h
o
s
e
v
a
l
u
e
s
a
r
e
a
l
l
e
q
u
a
l
t
o
1
-
1
×
100
%
.If there are more than one faulty component in power system, FIED can be presented as follows:(14)
FIED
=
t
h
e
n
u
m
b
e
r
o
f
e
l
e
m
e
n
t
s
w
h
o
s
e
v
a
l
u
e
s
a
r
e
F
A
+
F
0
e
q
u
a
l
t
o
1
-
1
W
W
i
W
i
i
W
a
l
l
g
r
e
a
t
e
r
t
h
a
n
1
W
W
i
W
i
i
W
i
n
F
A
+
F
0
∨
F
B
+
F
0
∨
F
C
+
F
0
∨
⋯
W
W
W
i
i
W
×
t
h
e
n
u
m
b
e
r
o
f
e
l
e
m
e
n
t
s
W
W
W
W
i
i
W
w
h
o
s
e
v
a
l
u
e
s
a
r
e
a
l
l
e
q
u
a
l
t
o
1
-
1
×
100
%
,where the symbol “∨” denotes that the greater value is selected for the values of corresponding elements.To locate out the fault components in systems, we need to compare the size of FIED. Once FIED satisfiesF
I
E
D
=
100
%, the fault components may be determined by formulas (13) and (14).


## 3.5. Judging Method of Malfunction and Rejection for Protective Relays and Circuit Breakers

According to the assumption and analysis above, the judging rule is designed as follows.(1) If the value of element of the vector (F
A
-
F
0) is negative, then the corresponding protection of protective relays or action of circuit breakers is regarded as malfunction for single-component fault. Otherwise, the protection of protective relays or trip of circuit breakers is regarded as rejecting action.(2) If there is more than one fault component in the system, the judging rule is as follows. The protection of protective relays or action of circuit breakers is regarded as malfunction while the value of the vectorF
A
∨
F
B
∨
F
C
∨
⋯
-
F
0 is negative; otherwise the action is regarded as rejecting action.Notice that when failures happen, the main protection of PRs is first implemented to trip the circuit breakers. If the main protective relays do not work, then the backup protective relays begin to work to trip the circuit breakers. Therefore, the situation of the backup PRs deduced by (1) and (2) was the rejecting action; it may be that the action of PRs has already tripped the action of CBs. So we concern only the main protective relays and the circuit breakers and then give priority to main PRs when the malfunction and rejecting action are considered in this paper.


## 3.6. Fault Diagnosis Model of Single Component

The problem of recognition for single-fault component can be solved by seeking for the fault assumptions, which can best explain the alarming information and can be expressed by the following maximization problem:(15)
R
i
=
1
n
∑
k
=
1
n
min
i
⁡
min
k
⁡
F
0
k
-
F
i
k
+
ξ
max
i
⁡
max
k
⁡
F
0
k
-
F
i
k
-
1
W
W
W
W
i
W
+
ξ
max
i
⁡
max
k
⁡
F
0
k
-
F
i
k
W
W
W
W
W
×
F
0
k
-
F
i
k
max
k
⁡
F
0
k
-
F
i
k
W
W
W
W
W
i
W
+
ξ
max
i
⁡
max
k
⁡
F
0
k
-
F
i
k
-
1
,where i stands for the total number of the single-fault, and its value is determined by adding the number of components and the number of lines; n denotes the number of vector elements, and its value is the number of protective relays and circuit breakers; ξ is on behalf of the resolution coefficient; F
0
(
k
) indicates the quarantined faulty state vector; F
i
(
k
) means the fault state vector of single-line or single component and also is called the single-fault state vector; R
i represents the similarity degree between the quarantined fault state vector and single-fault state vector and also is called correlation degree.Notice that if the state vector of single-fault componentA
,
B
,
…
,
X is F
A
,
F
B
,
…
,
F
X, respectively, then the dimension of these state vectors is 1
∗
(
m
+
n
) and the dimension of the quarantined fault state vector F
0 is also 1
∗
(
m
+
n
). Meanwhile, we may calculate the correlation degree and sort the calculated results. Based on above discussion, the recognition rule of fault component is given as follows.Rule 1.
If the correlation degree between the reference state vectorF
0 and the compare state vector F
i
(
i
=
A
,
B
,
…
,
X
) is the largest, the corresponding part is fault component or fault line.Rule 2.
Let the threshold of fault diagnosis be the mean of correlation degree. If the correlation degree between the reference state vectorF
0 and new compare state vector F
N is more than the threshold, the part is regarded as the fault component or fault line.Therefore, the possibility that fault happens may be determined for every single component according to Rules1 and 2.


## 3.7. Diagnostic Procedure for Fault Components

In the actual project, there is more than one fault component when faults occur. But the diagnostic procedure may be designed on the basis of the fault diagnosis model of single component. So the specific steps of fault diagnosis are shown as follows.


### 3.7.1. Data Preprocessing

The data, which contains the information of protective relays and circuit breakers and the topology of the power system, is written in the system.


### 3.7.2. Recognition for the Power Supply Interrupted Districts

The blackout area of power systems is recognized, and the tripping principle of protective relays and circuit breaker in blackout area is listed.


### 3.7.3. Recognition of State Information

Under single-fault state, the state information of protective relays and circuit breaker is obtained in accordance with the working principle of protective relays and circuit breakers.


### 3.7.4. Construction for Fault State Vectors

The state vector of fault mode for each component is given by the assignment principle of vector elements, and then the quarantined fault state vector is loaded.


### 3.7.5. Calculating and Sorting the Correlation Degree

The correlation coefficients are calculated, respectively, and then these results are sorted. If the correlation degree between quarantined state vector and reference state vector is the highest, the probability that component malfunctions is also the largest.


### 3.7.6. Computing Fault Information Explanation Degree

FIED is calculated by formulas (13) and (14) and is sorted. If FIED of some fault components satisfies F
I
E
D
=
100
%, go to next step.


### 3.7.7. Fault Diagnosis Analysis

According to Rules1 and 2, the fault component can be recognized and diagnosed. And then judge the malfunction and rejection for protective relays and circuit breakers.The specific flowchart for fault diagnosis is shown as in Figure1.Figure 1
The flowchart for fault diagnosis.


## 3.7.1. Data Preprocessing

The data, which contains the information of protective relays and circuit breakers and the topology of the power system, is written in the system.


## 3.7.2. Recognition for the Power Supply Interrupted Districts

The blackout area of power systems is recognized, and the tripping principle of protective relays and circuit breaker in blackout area is listed.


## 3.7.3. Recognition of State Information

Under single-fault state, the state information of protective relays and circuit breaker is obtained in accordance with the working principle of protective relays and circuit breakers.


## 3.7.4. Construction for Fault State Vectors

The state vector of fault mode for each component is given by the assignment principle of vector elements, and then the quarantined fault state vector is loaded.


## 3.7.5. Calculating and Sorting the Correlation Degree

The correlation coefficients are calculated, respectively, and then these results are sorted. If the correlation degree between quarantined state vector and reference state vector is the highest, the probability that component malfunctions is also the largest.


## 3.7.6. Computing Fault Information Explanation Degree

FIED is calculated by formulas (13) and (14) and is sorted. If FIED of some fault components satisfies F
I
E
D
=
100
%, go to next step.


## 3.7.7. Fault Diagnosis Analysis

According to Rules1 and 2, the fault component can be recognized and diagnosed. And then judge the malfunction and rejection for protective relays and circuit breakers.The specific flowchart for fault diagnosis is shown as in Figure1.Figure 1
The flowchart for fault diagnosis.


## 4. Experiment and Its Results Analysis

Firstly, regarding the simple circuit structure shown in Figure2 as the power supply interrupted districts, the following sections will explain and demonstrate in detail how to use the working principle of protection action to identify the state vectors of fault components.Figure 2
Framework of simple circuit structure.Obviously, the simple system shown in Figure2 consists of 4 elements (i.e., A, B, L1, and L2) and 5 circuit breakers (i.e., CB1, CB2, CB3, CB4, and CB5). So there are 14 protection actions, respectively, that is, Am, Bm, L1Am, L1Bm, L2Bm, L2Cm, L1Ap, L1Bp, L2Bp, L2Cp, L1As, L1Bs, L2Bs, and L2Cs. To make things simple, these protection actions are represented symbolically by r
1
,
r
2
,
…
,
r
14. Note that, in this example, A and B are Bus, L represents Line, the subscript m denotes the primary protection, s stands for the first backup protection, and p is on behalf of the second backup protection. Through analysis of the circuit structure, the working principle of the primary protection and the backup protection is listed in Table 1.Table 1
Working principle of protective relays and circuit breakers.
Number
Name
Corresponding action of protective relays and circuit breakers
r
1
Am
Fault: A; action: Am; tripping: CB1 and CB2
r
2
Bm
Fault: B; action: Bm; tripping: CB3 and CB4
r
3
L1Am
Fault: L1; action: L1Am; tripping: CB2
r
4
L1Bm
Fault: L1; action: L1Bm; tripping: CB3
r
5
L2Bm
Fault: L2; action: L2Bm; tripping: CB4
r
6
L2Cm
Fault: L2; action: L2Cm; tripping: CB5
r
7
L1Ap
Fault: L1; unaction: L1Am; action: L1Ap; tripping: CB2
r
8
L1Bp
Fault: L1; unaction: L1Bm; action: L1Bp; tripping: CB3
r
9
L2Bp
Fault: L2; unaction: L2Bm; action: L2Bp; tripping: CB4
r
10
L2Cp
Fault: L2; unaction: L2Cm; action: L2Cp; tripping: CB5
r
11
L1As
Fault: B; untripping: CB3; action: L1As; tripping: CB2; or fault: L2; untripping: CB3 and CB4; action: L1As; tripping: CB2
r
12
L1Bs
Fault: A; untripping: CB2; action: L1Bs; tripping: CB3
r
13
L2Bs
Fault: C; untripping: CB5; action: L2Bs; tripping: CB4
r
14
L2Cs
Fault: B; untripping: CB4; action: L2Cs; tripping: CB5;or fault: L1; untripping: CB3 and CB4; action: L2Cs; tripping: CB5According to the working principle of primary protection and backup protection shown in Table1, the state information of protective relays and circuit breakers is deduced as shown in Table 2.Table 2
State information of protective relays and circuit breakers.
Fault nodes
The state information for protective relays and circuit break
A
r
1 = 1, CB1 = 1 and CB2 = 1 or r
12 = 1, CB2 = 0 and CB3 = 1
B
r
2 = 1, CB3 = 1 and CB4 = 1 or r
11 = 1, CB3 = 0 and CB2 = 1 or r
14 = 1, CB4 = 0 and CB5 = 1
L1
r
3 = 1 and CB2 = 1 or r
4 = 1 and CB3 = 1 or r
7 = 1, r
3 = 0, and CB2 = 1 or r
8 = 1, r
4 = 0, and CB3 = 1 or r
14 = 1, CB3 = 0, CB4 = 0, and CB5 = 1
L2
r
5 = 1 and CB4 = 1 or r
6 = 1 and CB5 = 1 or r
9 = 1, r
5 = 0, and CB4 = 1 or r
10 = 1, r
6 = 0, and CB5 = 1 or r
11 = 1, CB3 = 0, CB4 = 0, and CB2 = 1Letc
i
(
i
=
1,2
,
…
,
5
) denote the action of five circuit breakers. Thus, by combining with the state information shown in Table 2, the fault state vector of every component may be determined using the assignment principle introduced in Section 3.3. Therefore, the fault state vector of every element is shown in Table 3.Table 3
State vectors of fault component.
Fault node
The state vectors(
r
1
,
r
2
,
…
,
r
14,c
1
,
c
2
,
…
,
c
5
)
Simple forms
A
(1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0.5, 1, 0, 0)
(r
1, r
12, c
1, c
2/2, c
3)
B
(0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0.5, 0.5, 1)
(r
2, r
11, r
14, c
2, c
3/2, c
4/2, c
5)
L1
(0, 0, 0.5, 0.5, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0.5, 0, 1)
(r
3/2, r
4/2, r
7, r
8, r
14, c
2, c
3/2, c
5)
L2
(0, 0, 0, 0, 0.5, 0.5, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0.5, 1)
(r
5/2, r
6/2, r
9, r
10, r
11, c
2, c
4/2, c
5)If a quarantined fault state vector obtained by alarming information is(
r
1
,
r
5
,
r
6
,
c
1
,
c
2
,
c
4
,
c
5
) and it needs to determine which node is in fault, we need to compare the size of correlation coefficient between quarantined fault state and reference state vector. According the data in Table 3, the correlation degree may be computed by formula (15); that is,(16)
r
A
=
0.8333
,
r
B
=
0.8247
,
r
L
1
=
0.7807
,
r
L
2
=
0.8509
.The sorting result is(17)
r
L
2
>
r
A
>
r
B
>
r
L
1
.So we can infer that the probability that the Line L2 malfunctions is maximum, and then the number of fault elements (FEN) is 1; that is, there is only one fault component L2 in systems. In this case, we may get(18)
F
I
E
D
=
5
÷
7
×
100
%
=
71.4
%
.Obviously, the fault component cannot be completely determined. Through analysis of the information data, when there are two fault components L2 and A, FIED is(19)
F
I
E
D
=
7
÷
7
×
100
%
=
100
%
.Therefore, the quarantined state vector indicates that two components L2 and A are failure.Meanwhile,(20)
F
L
2
∨
F
A
=
r
1
,
r
5
2
,
r
6
2
,
r
9
,
r
10
,
r
11
,
r
12
,
c
1
,
c
2
,
c
3
,
c
4
2
,
c
5
,
(21)
F
L
2
∨
F
A
-
F
0
=
-
r
5
2
,
-
r
6
2
,
r
9
,
r
10
,
r
11
,
r
12
,
c
3
,
-
c
4
.By the judging rule of malfunction and rejection, the primary protection which corresponds to r
5 and r
6 and the action c
4 which corresponds to the circuit breaker CB4 is regarded as malfunction. Note that the alarming information is (
r
1
,
r
5
,
r
6
,
c
1
,
c
2
,
c
4
,
c
5
), and the action c
3 shown in formula (21) is a rejection. In other words, the information which corresponds to the action of the circuit breaker CB3 has lost. Namely, the more accurate diagnosis results can be obtained by the algorithm even if some information has lost in systems.The above computing process has primarily illustrated the calculated steps of fault diagnosis algorithm based on Gray systems theory. And the results indicate that the algorithm is effective and reasonable. To further explain and verify the rationality and effectiveness of the method mentioned in this paper, we take the classic system structure of local power relay protection shown by [9], which contains 28 components, 84 protective relays, and 40 circuit breakers, to analyze the diagnosis procedure. The structure of the local power systems is shown in Figure 3.Figure 3
Classic system structure of the local power relay protection.Where A and B denote Bus, L is Line, T represents Transformer, and CB is circuit breaker.As can be seen in Figure3, 28 components are A1, A2, A3, A4, T1, T2, T3, T4, T5, T6, T7, T8, B1, B2, B3, B4, B5, B6, B7, B8, L1, L2, L3, L4, L5, L6, L7, and L8; 84 protections are consisting of 36 primary protections and 48 backup protections, where 36 primary protections are, respectively, A1m, B1m, B2m, B3m, B4m, B5m, B6m, B7m, B8m, L1Sm, L1Rm, L2Sm, L2Rm, L3Sm, L3Rm, L4Sm, L4Rm, L5Sm, L5Rm, L6Sm, L6Rm, L7Sm, L7Rm, L8Sm, and L8Rm; 48 backup protections are as follows: T1p, T2p, T3p, T4p, T5p, T6p, T7p, T8p, T1s, T2s, T3s, T4s, T5s, T6s, T7s, T8s, L1Sp, L1Rp, L2Sp, L2Rp, L3Sp, L3Rp, L4Sp, L4Rp, L5Sp, L5Rp, L6Sp, L6Rp, L7Sp, L7Rp, L8Sp, L8Rp, L1Ss, L1Rs, L2Ss, L2Rs, L3Ss, L3Rs, L4Ss, L4Rs, L5Ss, L5Rs, L6Ss, L6Rs, L7Ss, L7Rs, L8Ss, and L8Rs.And the meaning of these symbols is, respectively, as follows: S is Sending End of Line and R is Receiving End,m is primary protection, p represents first backup protection, and s is second backup protection. The working principle of these protections may be seen in the Resources [11, 13] section.By comparing the method proposed in this paper with the methods in [13, 14], we are testing out 4 most complicated circumstances in case of incomplete information and complete information to verify the effectiveness of the method. The diagnosis results are listed as in Table 4.Table 4
Part of the test results.
Test sequences
Relay action information
Missing information
Diagnosis result
Judgment for malfunction and rejection
Diagnosis result in [11]
Diagnosis result in [12]
1
Behavior: T5s, T6s; tripping:   C21, C22, C23, C24, and C25
unexisting
A3
Malfunction: A3m; rejection (tripping): C21, C24, and C25
A3
A3
2
Behavior: B1m, L2Rs, L4Rs, C4; tripping: C5, C7, C9, C12, and C27
unexisting
B1
Rejection (tripping): C6
B1
B1
3
Behavior: B1m, L4Rs; tripping: C4, C7, C9, C12, and C27
L2Rs, C5
B1
Rejection (tripping): C5 and C6
B1
4
Behavior: B1m, T1m, T2m, L1sm, L2RP; tripping: C2, C3, C4, C5, C7, C9, C11, and C28
T7s, C6
L1, B1, T1, T2
Rejection (tripping): L1Rm, C6; malfunction: C28
L1, B1, T1, T2As is shown in Table4, in case of complete information, the diagnosis result for test sequence 1 and sequence 2 is identical to the results in [11, 12]; in case of incomplete information, the diagnosis result for test sequence 3 and test sequence 4 is identical to the results in [11], but it is not identical to the result in [12]. This is because there are some lost information in test sequence 3 and test sequence 4.Therefore, to further test this algorithm proposed in this paper, we apply the presented algorithm in this paper and the Bayesian algorithm in [12] to diagnose the power supply in intelligent traffic systems (ITS). The detecting information of the power supply in ITS contains the PT break-phase, CT signal, volts d.c. and current, synchronous signal, temperature of sensors, and the Airborne capacity. In real application, the fitting curve between using the two methods is as in Figure 4.Figure 4
(a) The fitting curve of test sequence (0–780) using the two methods; (b) The fitting curve of test sequence (0–301) using the two methods; (c) The fitting curve of test sequence (301–500) using the two methods (Note: + indicates the test sequence using Bayesian algorithm in [12]; ○ indicates the test sequence using Gray theory algorithm).
(a)
(b)
(c)Figure4(a) showed the whole effect of diagnosis procedure. Correspondingly, Figures 4(b) and 4(c) described the part effect of diagnosis algorithm. This indicates that the error between using the two methods is very small. This indicates that the algorithm presented in this paper is reasonable.However, the running stability of real system needs to be ensured for any algorithm, whether or not the algorithm is effective and reasonable. That is to say, the stability of the diagnosis procedure should to be analyzed. For verifying the rationality of the designed algorithm, the anticipant precision threshold is set as 0.01 and every running test time is set as 1000 epochs. The total running performance is shown as in Figure5.Figure 5
The total running performance chart of fault diagnosis for the power supply in ITS.The running effect of diagnosis system indicated that the abnormal phenomenon of the security running took place at 20th second, and then the parameter of systems needs to be adjusted to ensure the systems normal running. Figure5 displays that the abnormal phenomena of the designed diagnosis systems seem to change small in warning, and yet the running performance is very stable during the whole running process. To find out why the danger happens, we take the running performance chart from 0th second to 40th second as in Figure 6.Figure 6
The part running performance chart of fault diagnosis for the power supply in ITS.Figure6 displays that the interior change of diagnosis system was quite rapid. Simultaneously, there exists a glacis from the initial warning to normal running. As a result, the running effect is in accord with the practical situation. This is because certain period of time is used to adjust parameter of the power supply for ensuring the final stable running of the algorithm.


## 5. Conclusion

Based on the theory of Gray system, a fault diagnosis method in power systems is proposed with respect to the diagnosis problem of incomplete information for protective relays devices in power system, and the analysis of overall architecture and the diagnosis process are conducted. Some simple experiments show that the method presented in this paper is effective and can also carry out the same fault diagnosis task that another traditional diagnosis method completed. Through comparing these diagnosis results, we know that the method may diagnose the fault in case of incomplete and complete information. In addition, through the practical engineering application in ITS, the stability of diagnosis procedure was also analyzed. The simulation result demonstrated that the running stability of the diagnosis procedure is practical and consistent with the operation requirements.However, owing to the fact that information data comes from different parts of power supply systems in ITS, the running situation of diagnosis systems was influenced. Thereby, the investigation about how to logically distribute the information of power supply device in ITS remains an interesting area for further research. Meanwhile, the method and algorithm presented in this paper can only be used to diagnose and determine which component malfunctions after the abnormal phenomenon of power systems has occurred. But in practical application of power supply devices in ITS, the users always hope that the diagnosis systems can forecast the faults before the fault states of components will happen. So, how to fuse the incomplete and complete information of protective relays and circuit breakers to forecast the faults by GM model, which provides the maintaining decision for the components and lines, is also very important problem in power systems.



---

*Source: 971257-2015-06-08.xml*