To control the response of the gas absorber, the range of linguistic values of the output of each feedback fuzzy controller is tuned between (– 3) and (3). Comparison between the output response with fuzzy controller (when the number of linguistic values of the controller input – output pair is three) and without controller is illustrated in Fig.7 and Fig.8. In Fig.7, the controller is tuned to interfere the natural decay of the system. In Fig. 8, the fuzzy controller is adjusted to improve the response of the gas absorber. It is noted that the response of controlled system has less overshoot, less steady state error and faster compared

**Figure 8.** The controller is tuned to improve the response of the system.

## *2.1.5. Discussion*

This section has presented a new and simple fuzzy controller for a gas absorber system to enhance the response of the output. The simulation results have shown that the controller guarantees well-damped behaviour of the controlled gas absorber system.

New Areas in Fuzzy Application 395

and providing qualitative "rules of thumb" (qualitative control rules in terms of vague and

A concerted effort has been made to formally reduce the size of the fuzzy rule base to make fuzzy control attractive to interconnected systems. Two of the difficulties with the design of

The properties that a fuzzy membership function is used to characterize are usually fuzzy. Therefore, we may use different membership functions to characterize the same description. Conceptually, there are two approaches to determine a membership function. The first approach is to use the knowledge of human experts. Usually this approach can only give a rough formula of the membership function; fine-tuning is required. In the second approach, data are collected from various sensors to determine the member ship functions. Specifically, the structures of the membership functions are specified first, then fine-tuning of the membership function parameters should be implemented based on the collected data

In this section, we contribute to the further development of intelligent control techniques of interconnected systems. It provides a new approach to fuzzy control design for interconnected system. The approach consists of two stages: In the first stage, a group of local state estimator is constructed to generate the data base of input-output pairs. In the second stage, an array of feedback fuzzy controllers is designed and implemented to ensure the asymptotic satiability of the interconnected system. Simulation studies on a large-scale system with unstable eigenvalus are carried out to illustrate the features and capability of

In the sequel, the terms large-scale and interconnected are used interchangeable. The term large scale system (LSS) does not have a unique established meaning, but it covers systems that possess several particular feature, such as multiple subsystem, [14,17] multiple control, multiple objectives, decentralized and/or hierarchical information structures. Any LSS includes many variables but their control is faced by a well-know fact [16] that the states are

Many authors have considered the state estimation of large-scale systems in input decentralized fashion. Here we summarize one convenient algorithm [15]. Let the state

N

i ii ii ij j x t Ax t Bu t G x (14)

i ii y t C x t , i, j 1, 2, N (15)

fuzzy sentences).

[8].

this new approach.

any fuzzy control system are:

The choice of fuzzy rules.

The shape of the membership functions.

*2.2.2. State estimation of interconnected systems*

model of the ith subsystem described by

not always available for measurement and state must be estimated.

## **2.2. Large scale fuzzy controller**

In this section, we shall develop a new approach to the control of interconnected system using fuzzy system theory. The approach shall be based on incorporating a group of local estimators on the system level to generate the input-output database. An array of feedback fuzzy controllers shall then be designed to ensure the asymptotic stability of the closed loop system. The developed technique shall be applied to an unstable large-scale system and extensive simulation studies shall be carried out to illustrate the potential of this new approach..

## *2.2.1. Background*

In control engineering research, problems of decentralized control and stabilization of interconnected systems are receiving considerable interest in recent years [14,15] where most of the effort is focused on dealing with the interaction patterns. It is concluded that a systematic approach to deal with the problems of interconnected systems is twofold: first is to base the analysis and design effort on the subsystem level using conventional control methods and second is to deal with interactions effectively. These methods are facilitated, in general, by virtue of several mathematical tools including linearization, delay approximation, decomposition and model reduction. This constitutes the so-called modelbased control system approach for which we have seen numerous techniques [16]. Most of the available results have so far overlooked the operational knowledge of the interconnected system under consideration. In [17], a knowledge-based control system approach has been suggested to deal with the analysis and design problems of interconnected systems by incorporating both the simplest available model as well as the best available knowledge about the system. For single physical systems, one of the earlier efforts along this direction has been on the development of an expert learning system [18-19]. An alternative approach has been on integrating elements of discrete event systems with differential equations [20]. A practically-supported third approach has been through the use of fuzzy logic control by successfully applying fuzzy sets and systems theory [21].

For interconnected systems, the foregoing approach motivates the research into intelligent control by combining techniques of control and systems theory with those from artificial intelligence. The main focus should be on integrating a knowledge base, an approximate (humanlike) reasoning and/or a learning process within a hierarchical structure.

Fuzzy logic controllers [23-25] are generally considered applicable to plants that are mathematically poorly understood (there is no acceptable mathematical model for the plant) and where experienced human operators are available for satisfactorily controlling the plant and providing qualitative "rules of thumb" (qualitative control rules in terms of vague and fuzzy sentences).

A concerted effort has been made to formally reduce the size of the fuzzy rule base to make fuzzy control attractive to interconnected systems. Two of the difficulties with the design of any fuzzy control system are:

394 Fuzzy Controllers – Recent Advances in Theory and Applications

successfully applying fuzzy sets and systems theory [21].

**2.2. Large scale fuzzy controller**

This section has presented a new and simple fuzzy controller for a gas absorber system to enhance the response of the output. The simulation results have shown that the controller

In this section, we shall develop a new approach to the control of interconnected system using fuzzy system theory. The approach shall be based on incorporating a group of local estimators on the system level to generate the input-output database. An array of feedback fuzzy controllers shall then be designed to ensure the asymptotic stability of the closed loop system. The developed technique shall be applied to an unstable large-scale system and extensive simulation studies shall be carried out to illustrate the potential of this new

In control engineering research, problems of decentralized control and stabilization of interconnected systems are receiving considerable interest in recent years [14,15] where most of the effort is focused on dealing with the interaction patterns. It is concluded that a systematic approach to deal with the problems of interconnected systems is twofold: first is to base the analysis and design effort on the subsystem level using conventional control methods and second is to deal with interactions effectively. These methods are facilitated, in general, by virtue of several mathematical tools including linearization, delay approximation, decomposition and model reduction. This constitutes the so-called modelbased control system approach for which we have seen numerous techniques [16]. Most of the available results have so far overlooked the operational knowledge of the interconnected system under consideration. In [17], a knowledge-based control system approach has been suggested to deal with the analysis and design problems of interconnected systems by incorporating both the simplest available model as well as the best available knowledge about the system. For single physical systems, one of the earlier efforts along this direction has been on the development of an expert learning system [18-19]. An alternative approach has been on integrating elements of discrete event systems with differential equations [20]. A practically-supported third approach has been through the use of fuzzy logic control by

For interconnected systems, the foregoing approach motivates the research into intelligent control by combining techniques of control and systems theory with those from artificial intelligence. The main focus should be on integrating a knowledge base, an approximate

Fuzzy logic controllers [23-25] are generally considered applicable to plants that are mathematically poorly understood (there is no acceptable mathematical model for the plant) and where experienced human operators are available for satisfactorily controlling the plant

(humanlike) reasoning and/or a learning process within a hierarchical structure.

guarantees well-damped behaviour of the controlled gas absorber system.

*2.1.5. Discussion*

approach..

*2.2.1. Background*

The properties that a fuzzy membership function is used to characterize are usually fuzzy. Therefore, we may use different membership functions to characterize the same description.

Conceptually, there are two approaches to determine a membership function. The first approach is to use the knowledge of human experts. Usually this approach can only give a rough formula of the membership function; fine-tuning is required. In the second approach, data are collected from various sensors to determine the member ship functions. Specifically, the structures of the membership functions are specified first, then fine-tuning of the membership function parameters should be implemented based on the collected data [8].

In this section, we contribute to the further development of intelligent control techniques of interconnected systems. It provides a new approach to fuzzy control design for interconnected system. The approach consists of two stages: In the first stage, a group of local state estimator is constructed to generate the data base of input-output pairs. In the second stage, an array of feedback fuzzy controllers is designed and implemented to ensure the asymptotic satiability of the interconnected system. Simulation studies on a large-scale system with unstable eigenvalus are carried out to illustrate the features and capability of this new approach.

## *2.2.2. State estimation of interconnected systems*

In the sequel, the terms large-scale and interconnected are used interchangeable. The term large scale system (LSS) does not have a unique established meaning, but it covers systems that possess several particular feature, such as multiple subsystem, [14,17] multiple control, multiple objectives, decentralized and/or hierarchical information structures. Any LSS includes many variables but their control is faced by a well-know fact [16] that the states are not always available for measurement and state must be estimated.

Many authors have considered the state estimation of large-scale systems in input decentralized fashion. Here we summarize one convenient algorithm [15]. Let the state model of the ith subsystem described by

$$\mathbf{x}\_{\mathbf{i}}\left(\mathbf{t}\right) = \mathbf{A}\_{\mathbf{i}}\mathbf{x}\_{\mathbf{i}}\left(\mathbf{t}\right) + \mathbf{B}\_{\mathbf{i}}\mathbf{u}\_{\mathbf{i}}\left(\mathbf{t}\right) + \boldsymbol{\Sigma}^{\mathbf{N}}\mathbf{G}\_{\mathbf{i}\parallel}\mathbf{x}\_{\mathbf{j}} \tag{14}$$

$$\mathbf{y}\_i(\mathbf{t}) = \mathbb{C}\_i \mathbf{x}\_i(\mathbf{t}), \text{i.j.} = \mathbf{1}, \text{2. } \dots \text{N} \tag{15}$$

Where all vectors and matrices are appropriately defined and gi(.) is the interaction function between the ith subsystem and the rest of the system. It is considered that (Ci, Ai) is completely observable for i = 1, 2, …….. N.

The following algorithm finds the optimal states of a large-scale system based on decentralized estimation and control [17]:

*Algorithm 1:*

*Step 1:*

Read the matrices Ai, Bi and select Qi 0 and Ri 0 as weighted matrix.

*Step 2:*

Solve the following 2N algebraic Riccati equations for Hi, Ki

$$\text{Hil}(\text{A}^{\text{T}}\_{\text{i}} + a\text{I}\_{\text{i}}) + (\text{A}\_{\text{i}} + a\text{I}\_{\text{i}})\text{H}\_{\text{i}} - \text{H}\_{\text{i}}\text{D}\_{\text{i}}\text{H}\_{\text{i}} + \text{Q}\_{\text{i}} = \text{0} \tag{16}$$

New Areas in Fuzzy Application 397

(20)

11 12 13 14 15 21 22 23 24 25 31 32 33 34 35 41 42 43 44 45 51 52 53 54 55

>

25 0.1 0.5

(21)

(22)

0.4 0.2 1 0.5 0.1 r1 r1 r1 r1 r1 A=

r2 r2 r2 r2 r2 0.1 0 0 0.2 0 r2 r2 r2 r2 r2 0 0.2 1 0 0.4 r2 r2 r2 r2 r2 0.6 0.1 0.25 2 0 r2 r2 r2 r2 r2 0.4 0.2 1 0.5 0.1

1100000000<PHONE_NUMBER>

0000011000<PHONE_NUMBER>

Which is considered to be composed of two-coupled subsystems; each of order 5. The coupling parameters are r1jk and r2jk where j and k take values of 1,2,3,4 and 5. In the sequel,

> A11 G12(r1) B1 x= x+ v G21(r2) A22 B2

For a typical values [4] of r115=-0.1, r124=0.1, r142=0.2, r222=0.1,r242=0.15,r251=0.11 and allothe values of coupling parameters are zeros, we examined the stability of the system by computing the eigenvalues of matrix A. They are {-1.0915, -1.0641, 0.477 + j0.0206, 0.477 – j0.00206, 0.022 + j0.0544, 0.022 – j0.0544, -1.8709 + j0.1713, -1.8709 – j0.1713, -1.9306 + j0.1413, -

(23)

r2 r2 r2 r2 r2 1.5 0.3 0.

B=

11 12 13 14 15

21 22 23 24 25 31 32 33 34 35 41 42 43 44 45 51 52 53 54 55

C=

we refer to the structure of the interconnected system model as:

Where G12(r1) and G21(r2) are the coupling mtrices.

$$\text{Ki}(\text{A}^{\text{T}}\_{\text{i}} + a\text{I}\_{\text{i}}) + (\text{A}\_{\text{i}} + a\text{I}\_{\text{i}})\text{K}\_{\text{i}} - \text{K}\_{\text{i}}\text{S}\_{\text{i}}\text{K}\_{\text{i}} + \text{Q}\_{\text{i}} = \text{0} \tag{17}$$

Where Di = CTiCi , Si = Bi R-1i BT

*Step 3:*

Integrate the following set of N simultaneous equation for ei(t), I = 1, 2 …. N, using the initial condition ei(0) = x i(0)

$$
\begin{pmatrix}
\dot{\mathbf{e}}\_{1} \\
\vdots \\
\dot{\mathbf{e}}\_{N}
\end{pmatrix} = \begin{pmatrix}
\mathbf{A}\_{1}\cdot\mathbf{S}\_{1}\mathbf{K}\_{1} & \cdots & \mathbf{G}\_{1N} \\
\vdots & \ddots & \\
\mathbf{G}\_{N\_{1}} & & \mathbf{A}\_{N}\cdot\mathbf{H}\_{N}\mathbf{D}\_{N}
\end{pmatrix} \begin{pmatrix}
\mathbf{e}\_{1} \\
\vdots \\
\mathbf{e}\_{N}
\end{pmatrix} + \begin{pmatrix}
\mathbf{B}\_{1}\mathbf{v}\_{1} \\
\vdots \\
\mathbf{B}\_{N}\mathbf{v}\_{N}
\end{pmatrix} \tag{18}
$$

*Step 4:*

Integrate the following set of n simultaneous equations for x 1(t), i = 1, 2, ….. N

$$
\begin{pmatrix}
\dot{\mathbf{x}}\_{1} \\
\vdots \\
\dot{\mathbf{x}}\_{N}
\end{pmatrix} = \begin{pmatrix}
\mathbf{A}\_{1}\cdot\mathbf{S}\_{1}\mathbf{K}\_{1} & \cdots & \mathbf{G}\_{1N} & \mathbf{S}\_{1}\mathbf{K}\_{1}\cdots & 0 \\
\vdots & & \ddots & \\
\mathbf{G}\_{N\_{1}} & & \mathbf{A}\_{N}\cdot\mathbf{S}\_{N}\mathbf{I}\_{N} & 0 & \mathbf{S}\_{N}\mathbf{K}\_{N}
\end{pmatrix} \begin{pmatrix}
\mathbf{x}\_{1} \\
\vdots \\
\mathbf{x}\_{N}
\end{pmatrix} + \begin{pmatrix}
\mathbf{B}\_{1}\mathbf{v}\_{1} \\
\vdots \\
\mathbf{B}\_{N}\mathbf{v}\_{N}
\end{pmatrix} \tag{19}
$$

ẋ N GN1 AN - SNIN 0 SNKN x N BNvN

*Step 5:*

Generate the input-output pairs {vi, ŷi = ci ˆ i}.

#### *2.2.3. Interconnected system*

Assume the following interconnected system of order 10 [17]:

11 12 13 14 15 21 22 23 24 25 31 32 33 34 35 41 42 43 44 45 51 52 53 54 55 11 12 13 14 15 -1.5 -0.3 -0.25 0.1 0.5 r1 r1 r1 r1 r1 0.1 0 0 0.2 0 r1 r1 r1 r1 r1 0 0.2 1 0 0.4 r1 r1 r1 r1 r1 0.6 0.1 0.25 2 0 r1 r1 r1 r1 r1 0.4 0.2 1 0.5 0.1 r1 r1 r1 r1 r1 A= r2 r2 r2 r2 r2 1.5 0.3 0. 21 22 23 24 25 31 32 33 34 35 41 42 43 44 45 51 52 53 54 55 25 0.1 0.5 r2 r2 r2 r2 r2 0.1 0 0 0.2 0 r2 r2 r2 r2 r2 0 0.2 1 0 0.4 r2 r2 r2 r2 r2 0.6 0.1 0.25 2 0 r2 r2 r2 r2 r2 0.4 0.2 1 0.5 0.1 (20)

$$\mathbf{B} = \begin{pmatrix} 0 & 0 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 1 \end{pmatrix} \tag{21}$$

$$\mathbf{C} = \begin{pmatrix} 1 & 1 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 1 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 1 & 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 & 1 & 0 & 1 \end{pmatrix} \tag{22}$$

Which is considered to be composed of two-coupled subsystems; each of order 5. The coupling parameters are r1jk and r2jk where j and k take values of 1,2,3,4 and 5. In the sequel, we refer to the structure of the interconnected system model as:

$$
\dot{\mathbf{x}} = \begin{pmatrix}
\text{A11} & \cdots & \text{G12(r1)} \\
\vdots & \ddots & \\
\text{G21(r2)} & & \text{A22}
\end{pmatrix} \mathbf{x} + \begin{pmatrix}
\text{B1} \\
\vdots \\
\text{B2}
\end{pmatrix} \mathbf{v} \\
\tag{23}
$$

Where G12(r1) and G21(r2) are the coupling mtrices.

396 Fuzzy Controllers – Recent Advances in Theory and Applications

completely observable for i = 1, 2, …….. N.

decentralized estimation and control [17]:

Where Di = CTiCi , Si = Bi R-1i BT

condition ei(0) = x i(0)

*Algorithm 1:*

*Step 1:*

*Step 2:*

*Step 3:*

*Step 4:*

*Step 5:*

Where all vectors and matrices are appropriately defined and gi(.) is the interaction function between the ith subsystem and the rest of the system. It is considered that (Ci, Ai) is

The following algorithm finds the optimal states of a large-scale system based on

T Hi(A I ) (A I )H H D H Q 0 i i i i i iii i

T Ki(A I ) (A I )K K S K Q 0 i i i i i ii i i

Integrate the following set of N simultaneous equation for ei(t), I = 1, 2 ….

Solve the following 2N algebraic Riccati equations for Hi, Ki

1

Integrate the following set of n simultaneous equations for x 1(t), i = 1, 2, ….. N

1

Assume the following interconnected system of order 10 [17]:

ẋ N GN1 AN - SNIN 0 SNKN x N BNvN

*2.2.3. Interconnected system*

Generate the input-output pairs {vi, ŷi = ci ˆ i}.

For a typical values [4] of r115=-0.1, r124=0.1, r142=0.2, r222=0.1,r242=0.15,r251=0.11 and allothe values of coupling parameters are zeros, we examined the stability of the system by computing the eigenvalues of matrix A. They are {-1.0915, -1.0641, 0.477 + j0.0206, 0.477 – j0.00206, 0.022 + j0.0544, 0.022 – j0.0544, -1.8709 + j0.1713, -1.8709 – j0.1713, -1.9306 + j0.1413, -

1.9306 – j0.1413}, and it is quite clear that there are four eigenvalues lie in the open right half of the complex plane and thus the interconnected system is unstable. Further, it is easy to check that the interconnected system is both controllable and observable.

New Areas in Fuzzy Application 399

Ni+1 at zero.

*2.2.5. Design of an array of fuzzy controller*

*Step 1:*

*Step 2:*

Mi

*Step 3:*

*Step 4:*

Ni+2, ---, Mi

we use Ni fuzzy set M1i, ---, Mi

and yb-L of the fuzzy set ki

design the fuzzy controller.

The following 2Ni+1 rules are considered

**Figure 11.** Block diagram of the proposed fuzzy feedback controller array.

of the respective subsystem outputs, where i = 1, 2, 3, 4.

In order to build each fuzzy controller, the following steps are implemented:

The range of the inputs to each fuzzy controller [i, i] are driven from the estimated value

2N+1 fuzzy set MLi in [i, i] that are normal, consistent and complete with triangular membership functions [24], are defined for each controller, where L = 1, 2--, 2Ni+1. That is

2Ni+1 to cover the positive internal (0, ], and the center of fuzzy set Mi

ai i bi i i IF y is M or y is M then u is K

Where L = 1, 2,---, 2Ni+1, and ai, bi are the input to the fuzzy controller i, and the center yaiL

Product inference engine, singleton fuzzyfier, and center average defuzzifier are selected to

L are chosen such that

ai bi i

y and y 0 for1= N +1

LL L

0 for1=1, ,N

0 for1= N +2, ,2N +1

i

(24)

i i

Ni to cover the negative internal [I, 0), the other Ni fuzzy sets

We are going to treat the interconnected system at hand as being composed of two identical and coupled subsystems. The control system to be designed is such that each subsystem has its own fuzzy negative feedback controller which its input being the output of the respective subsystem (Fig. 11). Each subsystem fuzzy controller is constructed using two fuzzy systems.

## *2.2.4. Estimation of the system state variables and outputs*

A Matlab program is written to implement the computational Algorithm (1) of section 1.2.2 on the interconnected system. Different positive and negative step input are applied to estimate the outputs. The results of two cases are illustrated in Fig. 9 and Fig. 10 It is observed that the outputs tend to track conveniently the input signals.

**Figure 9.** Simulation Results for case 1.

**Figure 10.** Simulation Results for case 2.

## *2.2.5. Design of an array of fuzzy controller*

We are going to treat the interconnected system at hand as being composed of two identical and coupled subsystems. The control system to be designed is such that each subsystem has its own fuzzy negative feedback controller which its input being the output of the respective subsystem (Fig. 11). Each subsystem fuzzy controller is constructed using two fuzzy systems.

**Figure 11.** Block diagram of the proposed fuzzy feedback controller array.

In order to build each fuzzy controller, the following steps are implemented:

*Step 1:*

398 Fuzzy Controllers – Recent Advances in Theory and Applications

**Figure 9.** Simulation Results for case 1.

**Figure 10.** Simulation Results for case 2.

1.9306 – j0.1413}, and it is quite clear that there are four eigenvalues lie in the open right half of the complex plane and thus the interconnected system is unstable. Further, it is easy to

A Matlab program is written to implement the computational Algorithm (1) of section 1.2.2 on the interconnected system. Different positive and negative step input are applied to estimate the outputs. The results of two cases are illustrated in Fig. 9 and Fig. 10 It is

check that the interconnected system is both controllable and observable.

observed that the outputs tend to track conveniently the input signals.

*2.2.4. Estimation of the system state variables and outputs*

The range of the inputs to each fuzzy controller [i, i] are driven from the estimated value of the respective subsystem outputs, where i = 1, 2, 3, 4.

#### *Step 2:*

2N+1 fuzzy set MLi in [i, i] that are normal, consistent and complete with triangular membership functions [24], are defined for each controller, where L = 1, 2--, 2Ni+1. That is we use Ni fuzzy set M1i, ---, Mi Ni to cover the negative internal [I, 0), the other Ni fuzzy sets Mi Ni+2, ---, Mi 2Ni+1 to cover the positive internal (0, ], and the center of fuzzy set Mi Ni+1 at zero.

*Step 3:*

The following 2Ni+1 rules are considered

$$\text{IF } \mathbf{y}\_{\text{ai}} \text{ is } \mathbf{M}\_{\mathbf{i}}^{\text{L}} \text{ or } \mathbf{y}\_{\text{bi}} \text{ is } \mathbf{M}\_{\mathbf{i}}^{\text{L}} \text{ then } \mathbf{u} \text{ is } \mathbf{K}\_{\mathbf{i}}^{\text{L}}.$$

Where L = 1, 2,---, 2Ni+1, and ai, bi are the input to the fuzzy controller i, and the center yaiL and yb-L of the fuzzy set ki L are chosen such that

$$\mathbf{y\_{ai}^{-\mathrm{L}}} \text{ and } \mathbf{y\_{bi}^{-\mathrm{L}}} \begin{cases} \leq -0 & \text{for} \, \mathbf{1} = \mathbf{1}, \cdots, \mathbf{N\_i} \\ = 0 & \text{for} \, \mathbf{1} = \mathbf{N\_i} + \mathbf{1} \\ \geq 0 & \text{for} \, \mathbf{1} = \mathbf{N\_i} + \mathbf{2}, \cdots, \mathbf{2N\_i} + \mathbf{1} \end{cases} \tag{24}$$

*Step 4:*

Product inference engine, singleton fuzzyfier, and center average defuzzifier are selected to design the fuzzy controller.

## *2.2.6. Simulation results*

The behaviour of the interconnected system outputs after implementing the fuzzy controllers with unity step function input are shown in Fig 12 and Fig: 13. It is clearly evident that the system becomes asymptotically stable by using the negative fuzzy feedback controller array.

New Areas in Fuzzy Application 401

The following figures illustrate the above test cases:

**Figure 14.** Case 2 Outputs y1 against y2

**Figure 15.** Case 2 Outputs y3 against y4

**Figure 16.** Case 3 Outputs y1 against y2

**Figure 17.** Case 3 Outputs y3 against y4

**Figure 18.** Case 4 Outputs y1 against y2

**Figure 12.** Outputs y1 against y2 y2

**Figure 13.** Outputs y3 against y4 y4

## *2.2.7. Performance of the proposed fuzzy feedback controller array*

Now, we examine the effect of coupling matrices on the performance of fuzzy controlled interconnected system. Five additional cases with deferent coupling ranks are implemented. Fine tuning of membership functions was required to adjust their ranges. The following table summarizes the test cases:

**Table 1.** Results summary for 6 test cases.

The following figures illustrate the above test cases:

400 Fuzzy Controllers – Recent Advances in Theory and Applications

The behaviour of the interconnected system outputs after implementing the fuzzy controllers with unity step function input are shown in Fig 12 and Fig: 13. It is clearly evident that the system becomes asymptotically stable by using the negative fuzzy feedback

*2.2.6. Simulation results*

**Figure 12.** Outputs y1 against y2 y2

**Figure 13.** Outputs y3 against y4 y4

table summarizes the test cases:

**2Norm**

**Table 1.** Results summary for 6 test cases.

**Case No. A11,A2**

*2.2.7. Performance of the proposed fuzzy feedback controller array*

**G12 Norm**

**G12 Sparsty**

Now, we examine the effect of coupling matrices on the performance of fuzzy controlled interconnected system. Five additional cases with deferent coupling ranks are implemented. Fine tuning of membership functions was required to adjust their ranges. The following

> **G21 Sparsty**

1 (Fig 12, 13) 2.2529 3/25 0.2 3/25 1.8028 Unstable Stable 2 (Fig. 14, 15) 2.2529 12/25 0.4712 3/25 0.1803 Unstable Stable 3 (Fig. 16, 17) 2.2529 3/25 .2 12/25 0.5341 Unstable Stable 4(Fig. 18, 19) 2.2529 1 3.0361 3/25 0.1803 Unstable Stable 5 (Fig. 20, 21) 2.2529 3/25 .2 1 3.0364 Unstable Stable 6 (Fig.22, 23) 2.2529 1 3.0361 1 3.0417 Unstable Stable

**G21 Norm**

**System Stability without controller**

**System Stability with controller**

controller array.

**Figure 15.** Case 2 Outputs y3 against y4

**Figure 16.** Case 3 Outputs y1 against y2

**Figure 18.** Case 4 Outputs y1 against y2

New Areas in Fuzzy Application 403

This section has developed a new fuzzy control design approach to interconnected system. It has been shown the approach consists of two stages: In stage 1, a group of local state estimator has been constructed to generate the input-output database. Then an array of feedback controllers has been designed and implemented to guarantee the overall asymptotically system stability. Extensive simulation studies have been performed to

This section presents the use of fuzzy logic technique to control the reactive power of a load and hence improve the power factor. A shunt compensator is used, which consists of a reactor in series with a phase controlled Thyristor bridge in parallel with a capacitor. The control composed of two independent fuzzy controllers, the Fuzzy Grouse Controller (FGC) and the Fuzzy Fine Controller (FFC). These fuzzy controllers are used to control the firing angle of the Thyristor Bridge until the source power factor reaches a desired value. Simulations for three different practical study cases are presented and the results show how

Power factor, nowadays, is an important issue. The over increasing utilization of power electronics in all kinds of industry applications and the severe standards requirements are pushing the research toward new solutions to keep the industrial power factor within

A mathematical formulation for the optimal reactive power control is discussed in [27]. Also, the optimized Fuzzy logic and digital PID controllers for a single phase power factor correction converter used in online UPS are demonstrated in [28, 29]. Parameters such as input membership functions, output membership functions, inference rules of fuzzy logic controller and proportional gain, integral gain and derivative PID controller are selected and optimized by genetic algorithms. In additional to that, the applications of a hybrid converter are implemented in [30]. The hybrid converter is basically a converter bridge with two GTOs. A control strategy based on learning is proposed. The learning structure is coded into Fuzzy conditional rules to train a neural network in manipulating the converter variables.

The common method of correction is by means of using static capacitors, whether connected in series or parallel. These are installed as a single unit or as a bank, to regulate the voltage and the reactive power flow at the point of connection. In shunt compensation arrangement, a reactor is connected in parallel with conventional capacitor compensation. The shunt reactor current can be varied via a phase-controlled thyristor bridge connected in series with the shunt reactor [31]. Changing the thyristor firing angle varies the amount. of the current flowing through the reactor. Thus, this thyristor controlled reactor acts as a variable reactor

*2.2.8. Discussion*

*2.3.1. Background*

certain ranges [26].

support the developed design approach.

the designed controller is fast and accurate.

**2.3. Power factor correction**

**Figure 19.** Case 4 Outputs y3 against y4

**Figure 20.** Case 5 Outputs y1 against y2

**Figure 21.** Case 5 Outputs y3 against y4

**Figure 22.** ase 6 Outputs y1 against y2

**Figure 23.** Case 6 Outputs y3 against y4

## *2.2.8. Discussion*

402 Fuzzy Controllers – Recent Advances in Theory and Applications

**Figure 19.** Case 4 Outputs y3 against y4

**Figure 20.** Case 5 Outputs y1 against y2

**Figure 21.** Case 5 Outputs y3 against y4

**Figure 22.** ase 6 Outputs y1 against y2

**Figure 23.** Case 6 Outputs y3 against y4

This section has developed a new fuzzy control design approach to interconnected system. It has been shown the approach consists of two stages: In stage 1, a group of local state estimator has been constructed to generate the input-output database. Then an array of feedback controllers has been designed and implemented to guarantee the overall asymptotically system stability. Extensive simulation studies have been performed to support the developed design approach.

## **2.3. Power factor correction**

This section presents the use of fuzzy logic technique to control the reactive power of a load and hence improve the power factor. A shunt compensator is used, which consists of a reactor in series with a phase controlled Thyristor bridge in parallel with a capacitor. The control composed of two independent fuzzy controllers, the Fuzzy Grouse Controller (FGC) and the Fuzzy Fine Controller (FFC). These fuzzy controllers are used to control the firing angle of the Thyristor Bridge until the source power factor reaches a desired value. Simulations for three different practical study cases are presented and the results show how the designed controller is fast and accurate.

### *2.3.1. Background*

Power factor, nowadays, is an important issue. The over increasing utilization of power electronics in all kinds of industry applications and the severe standards requirements are pushing the research toward new solutions to keep the industrial power factor within certain ranges [26].

A mathematical formulation for the optimal reactive power control is discussed in [27]. Also, the optimized Fuzzy logic and digital PID controllers for a single phase power factor correction converter used in online UPS are demonstrated in [28, 29]. Parameters such as input membership functions, output membership functions, inference rules of fuzzy logic controller and proportional gain, integral gain and derivative PID controller are selected and optimized by genetic algorithms. In additional to that, the applications of a hybrid converter are implemented in [30]. The hybrid converter is basically a converter bridge with two GTOs. A control strategy based on learning is proposed. The learning structure is coded into Fuzzy conditional rules to train a neural network in manipulating the converter variables.

The common method of correction is by means of using static capacitors, whether connected in series or parallel. These are installed as a single unit or as a bank, to regulate the voltage and the reactive power flow at the point of connection. In shunt compensation arrangement, a reactor is connected in parallel with conventional capacitor compensation. The shunt reactor current can be varied via a phase-controlled thyristor bridge connected in series with the shunt reactor [31]. Changing the thyristor firing angle varies the amount. of the current flowing through the reactor. Thus, this thyristor controlled reactor acts as a variable reactor [32]. By varying the firing angle, the total reactive power of the system can be controlled and hence the power factor of the system is improved. The use of fuzzy logic to derive a practical control scheme for a boost rectifier with reactive power factor correction was applied [33]. The control action is primarily derived from a set of linguistic rules used to generate a slow-varying DC signal to determine the PWM ramp function. The proposed technique uses lesser sensing elements than the classical rectifier. A new FACTS controller known as the Bootstrap variable inductance can emulate a variable positive and negative inductance [34]. The bootstrap variable inductance has a variety of FACTS applications such as series compensation of lines, fault current limiting, reactive-power control and load power factor improvement.

New Areas in Fuzzy Application 405

**Figure 24.** A detailed block diagram of the load, the fuzzy controllers, the source and the control

The following procedures describe the steps of designing the power factor controller:

The sizes of the 'Inductance' and the capacitance bank are selected such that their maximum available reactive power (in VAR) is equal to the maximum load reactive power (MLQ). Since the full source voltage is continuously applied on the capacitance bank (assuming the voltage drop across the short cable is negligible), then capacitance value 'C' of bank can be

> 2 source MLQ Farad

However, this is not the case for the inductance since it is connected in series with a full wave controller. The existence of such controller will limit the available maximum reactive consumed by the inductance depending on the firing angle action of the thyristor bridge. Thus, the effect of the single-phase full-wave circuit is considered when determining the size

As listed in [35], the general formulas for the r.m.s value of current (*Iload*) and voltage (*Vload*) across a load, comprises of inductance in series with resistance, controlled by single-phase

r( )(α / ω t) 1/2 l 2 1 I [ {sin( ) sin( ) } ] *source*

2 V 1 V [ {sin(ω t θ) sin(α θ)e }dω t] z π

*t e dt*

(26)

(27)

β r( )(α / ω t) source 1/2 l

*f* (25)

2 V

Where Vsource is "source" r.m.s. voltage in Volt and *f* is "source" frequency in Hz.

*C*

full-wave circuit. These formulas are given as follows:

π

*V*

*z*

scheme

*<IP_ADDRESS>. Elements sizing*

determined as follows:

of the inductance.

*Load*

Load

Here in this section, we shall focuses on the use of fuzzy logic sets to control the supply power factor. Unlike the conventional capacitive approach for power factor improvement in ac power system, the proposed control scheme has the advantage to avoid complexities associated with the non-linear mathematical modelling of switching converters. The proposed fuzzy logic controlling scheme consists of two controllers. The first controller (FGC) is designed to give the nearest desired value of the firing angle required to compensate for the source reactive power. However, the output correction of this controller is not efficiently accurate and hence, another correction step is needed. Thus, the second controller (FFC) checks the value of the source power factor and improves it above a pre-set desired value. The discussion includes the following:

## *2.3.2. Fuzzy power factor controller*

Figure 24 illustrates the block diagram for a single-phase variable load, with variable lag power factor, supplied by sinusoidal AC power source. Capacitor bank in parallel with inductance, controlled by single-phase full-wave circuit, are connected in parallel to the load in order to govern the total reactive power of the circuit. Fuzzy controller is designed to tune the firing angle of the single-phase full-wave circuit in order to adjust the voltage applied across the Inductance. By this way, the total source reactive power can be minimized to improve the source power factor. For three-phase circuits, one controller is dedicated for each phase. Here, we considered single-phase circuit for simplicity.

The structure of the controller contains two independent fuzzy controllers: Fuzzy Grouse Controller (FGC) and Fuzzy Fine Controller (FFC). FGC input is the load reactive power, the output of this controller gives the nearest value of the desired firing angle, which required to minimize the source power factor. FFC input is the source power factor. The output of FFC corrects the firing angle of the single-phase full-wave controller until the source power factor reaches or exceeds the pre-setted desired value.

**Figure 24.** A detailed block diagram of the load, the fuzzy controllers, the source and the control scheme

The following procedures describe the steps of designing the power factor controller:

#### *<IP_ADDRESS>. Elements sizing*

404 Fuzzy Controllers – Recent Advances in Theory and Applications

desired value. The discussion includes the following:

Description of the design steps of the power factor controller.

each phase. Here, we considered single-phase circuit for simplicity.

factor reaches or exceeds the pre-setted desired value.

Simulation of the proposed technique by testing it for three different study cases.

Figure 24 illustrates the block diagram for a single-phase variable load, with variable lag power factor, supplied by sinusoidal AC power source. Capacitor bank in parallel with inductance, controlled by single-phase full-wave circuit, are connected in parallel to the load in order to govern the total reactive power of the circuit. Fuzzy controller is designed to tune the firing angle of the single-phase full-wave circuit in order to adjust the voltage applied across the Inductance. By this way, the total source reactive power can be minimized to improve the source power factor. For three-phase circuits, one controller is dedicated for

The structure of the controller contains two independent fuzzy controllers: Fuzzy Grouse Controller (FGC) and Fuzzy Fine Controller (FFC). FGC input is the load reactive power, the output of this controller gives the nearest value of the desired firing angle, which required to minimize the source power factor. FFC input is the source power factor. The output of FFC corrects the firing angle of the single-phase full-wave controller until the source power

Illustration of the proposed control scheme.

*2.3.2. Fuzzy power factor controller*

power factor improvement.

[32]. By varying the firing angle, the total reactive power of the system can be controlled and hence the power factor of the system is improved. The use of fuzzy logic to derive a practical control scheme for a boost rectifier with reactive power factor correction was applied [33]. The control action is primarily derived from a set of linguistic rules used to generate a slow-varying DC signal to determine the PWM ramp function. The proposed technique uses lesser sensing elements than the classical rectifier. A new FACTS controller known as the Bootstrap variable inductance can emulate a variable positive and negative inductance [34]. The bootstrap variable inductance has a variety of FACTS applications such as series compensation of lines, fault current limiting, reactive-power control and load

Here in this section, we shall focuses on the use of fuzzy logic sets to control the supply power factor. Unlike the conventional capacitive approach for power factor improvement in ac power system, the proposed control scheme has the advantage to avoid complexities associated with the non-linear mathematical modelling of switching converters. The proposed fuzzy logic controlling scheme consists of two controllers. The first controller (FGC) is designed to give the nearest desired value of the firing angle required to compensate for the source reactive power. However, the output correction of this controller is not efficiently accurate and hence, another correction step is needed. Thus, the second controller (FFC) checks the value of the source power factor and improves it above a pre-set

The sizes of the 'Inductance' and the capacitance bank are selected such that their maximum available reactive power (in VAR) is equal to the maximum load reactive power (MLQ). Since the full source voltage is continuously applied on the capacitance bank (assuming the voltage drop across the short cable is negligible), then capacitance value 'C' of bank can be determined as follows:

$$C = \frac{\text{MLQ}}{2 \times \pi \times f \times \text{V}\_{\text{source}}^2} \text{ Farad} \tag{25}$$

Where Vsource is "source" r.m.s. voltage in Volt and *f* is "source" frequency in Hz.

However, this is not the case for the inductance since it is connected in series with a full wave controller. The existence of such controller will limit the available maximum reactive consumed by the inductance depending on the firing angle action of the thyristor bridge. Thus, the effect of the single-phase full-wave circuit is considered when determining the size of the inductance.

As listed in [35], the general formulas for the r.m.s value of current (*Iload*) and voltage (*Vload*) across a load, comprises of inductance in series with resistance, controlled by single-phase full-wave circuit. These formulas are given as follows:

$$\mathbf{I}\_{Load} = \frac{\sqrt{2 \times} V\_{source}}{z} \mathbf{\frac{1}{\pi}} \Big| \frac{1}{\pi} \Big| \{ \sin(\alpha \, t - \theta) - \sin(\alpha - \theta)e^{\binom{N}{1}(\alpha \,/\omega - \mathbf{t})} \} \text{dof } t \Big|^{1/2} \tag{26}$$

$$\mathbf{V}\_{\text{Load}} = \frac{\sqrt{2 \times} V\_{\text{source}}}{\mathbf{z}} \mathbf{l} \frac{1}{\pi} \Big| \frac{\beta}{\pi} \Big| \{ \sin(\omega \cdot \mathbf{t} - \Theta) - \sin(\alpha - \Theta) \mathbf{e}^{\int\_{\gamma}^{\zeta} (\alpha \cdot (\omega - \mathbf{t})}) d\omega \, \mathbf{t} \}^{1/2} \tag{27}$$

Where:

 =2f radian/second = Firing angle = Extinction angle (cut-off angle) = tan-1 (l/r) l = Load Inductance r = Load Resistance t = Time z = Load impedance

Since the conducting angle = - cannot exceed, the firing angle may not be less than and the control range of the firing angle is:

$$
\pi \succeq \alpha \ge \theta \tag{28}
$$

New Areas in Fuzzy Application 407

Select an initial value for the inductance as L = 1 *(*2*C)*

source VAR will determine the desired value of the inductance.

 If the source VAR equals zero*,* then stop the trials. If not, increase slightly the value of the inductance until the source VAR reaches the zero value. Then, this value of the

After using the Simulink model shown in Figure 25, the inductance value is obtained from the iteration. Then, the same model is used again to find the ratio of (QL/MLQ) for a range of () starting from /2 to degree where (QL) is the Inductance reactive power at certain value of (). This ratio (QL/MLQ) verses the firing angle () is then plotted. Typical curve is shown

After that, the resultant nonlinear curve is divided to N sections where each nonlinear section is approximated by the nearest linear section. Thus, N Gross Fuzzy controllers (FGC) are built, one FGC for each section, such that the input of each controller is equal to 1- (QLoad/(MLQ), where (QLoad) is the load reactive power. The output of the FGC is the firing angle (), which results in (QL) approximately equal to MLQ minus QLoad. The accuracy of

In order to design each FGC, let the full range of the input fuzzy membership function FGCMF(in)a for each controller is set to the (QL/MLQ) limits of the respective linearzed section, and the range of the output fuzzy membership function FGCMF(out)b for each controller is set to the () limits of the respective linearzed section. For example, in Case 3 (as will be explained later into details), the first controller is designed to take an action in case of (QL/MLQ) ratio reaches a value between rl=0.95 and 1and, the output firing angle of this controller shall take a value between 90 degree and 1 = 100 degree. However, the second controller is designed to take an action in case of (QL/MLQ) ratio reaches a value between r2 = 0.175 and rl = 0.950 and, the output firing angle of this controller shall take a

The resultant output of each fuzzy controller can be obtained based on the respective

linearzed section using the following functions for fuzzy implication process:

Run the model and record the source VAR.

*<IP_ADDRESS>. Fuzzy Grouse Controller (FGC) design*

**Figure 26.** Typical curve for (QL/MLQ) verses firing angle ()

the resultant (QL) depends on the curve linearization.

value between 1 = 100 degree and 2 = 105degree, and so on.

Triangle type for Membership functions.

Mamdany engine.

 Product for And. Max for Or.

in Figure 26.

For maximum available reactive power consumed by pure inductance MARP where r is ignored in the above equation 24 & 25), the maximum conducting angle is considered. Therefore, the value of the inductance can be calculated from the following equation assuming that =/2 and = (neglecting the impedance of the short Cable and assuming ideal thyristors):

$$\text{MARP} = \left(V\_{load}\right) \times \left(I\_{load}\right) \tag{29}$$

This MARP value must be equal to the maximum load reactive power MLQ for complete compensation. However, the two equations listed above are nonlinear and difficult to solve. Another simple procedure is needed when determining the size of the inductance L.

**Figure 25.** Simulink circuit used to determine and plot the ratio (QL/MLQ) verses firing angle ()in the range from /2 to degree.

Thus, the value of the inductance can be found by a practical and fast method using Simulink tool [26]. Figure 25, illustrates the proposed model that been used. The model consists of an ac power source connected to a reactive power compensator controlled by a thyristor circuit. The procedure to find the value of the inductance L can be summarized by the following steps:

Since the conducting angle = - cannot exceed, the firing angle may not be less than

For maximum available reactive power consumed by pure inductance MARP where r is ignored in the above equation 24 & 25), the maximum conducting angle is considered. Therefore, the value of the inductance can be calculated from the following equation assuming that =/2 and = (neglecting the impedance of the short Cable and assuming

This MARP value must be equal to the maximum load reactive power MLQ for complete compensation. However, the two equations listed above are nonlinear and difficult to solve.

**Figure 25.** Simulink circuit used to determine and plot the ratio (QL/MLQ) verses firing angle ()in the

Thus, the value of the inductance can be found by a practical and fast method using Simulink tool [26]. Figure 25, illustrates the proposed model that been used. The model consists of an ac power source connected to a reactive power compensator controlled by a thyristor circuit. The procedure to find the value of the inductance L can be summarized by

Another simple procedure is needed when determining the size of the inductance L.

παθ (28)

MARP *V I load load* (29)

Where:

= tan-1 (l/r)

t = Time

 =2f radian/second = Firing angle

l = Load Inductance r = Load Resistance

z = Load impedance

ideal thyristors):

range from /2 to degree.

the following steps:

= Extinction angle (cut-off angle)

and the control range of the firing angle is:

 If the source VAR equals zero*,* then stop the trials. If not, increase slightly the value of the inductance until the source VAR reaches the zero value. Then, this value of the source VAR will determine the desired value of the inductance.

#### *<IP_ADDRESS>. Fuzzy Grouse Controller (FGC) design*

After using the Simulink model shown in Figure 25, the inductance value is obtained from the iteration. Then, the same model is used again to find the ratio of (QL/MLQ) for a range of () starting from /2 to degree where (QL) is the Inductance reactive power at certain value of (). This ratio (QL/MLQ) verses the firing angle () is then plotted. Typical curve is shown in Figure 26.

**Figure 26.** Typical curve for (QL/MLQ) verses firing angle ()

After that, the resultant nonlinear curve is divided to N sections where each nonlinear section is approximated by the nearest linear section. Thus, N Gross Fuzzy controllers (FGC) are built, one FGC for each section, such that the input of each controller is equal to 1- (QLoad/(MLQ), where (QLoad) is the load reactive power. The output of the FGC is the firing angle (), which results in (QL) approximately equal to MLQ minus QLoad. The accuracy of the resultant (QL) depends on the curve linearization.

In order to design each FGC, let the full range of the input fuzzy membership function FGCMF(in)a for each controller is set to the (QL/MLQ) limits of the respective linearzed section, and the range of the output fuzzy membership function FGCMF(out)b for each controller is set to the () limits of the respective linearzed section. For example, in Case 3 (as will be explained later into details), the first controller is designed to take an action in case of (QL/MLQ) ratio reaches a value between rl=0.95 and 1and, the output firing angle of this controller shall take a value between 90 degree and 1 = 100 degree. However, the second controller is designed to take an action in case of (QL/MLQ) ratio reaches a value between r2 = 0.175 and rl = 0.950 and, the output firing angle of this controller shall take a value between 1 = 100 degree and 2 = 105degree, and so on.

The resultant output of each fuzzy controller can be obtained based on the respective linearzed section using the following functions for fuzzy implication process:

where,

a = 1,2, ... F (fuzzy membership function number) b = F-a+l(fuzzy membership function number) a and b are fuzzy membership function numbers F is the number of the fuzzy membership function.

Simulink automatic switching system SS 1 as shown in Figure 24 is designed to check the active range of (QL/MLQ) in order to select the proper FGC based on that value of (QL/MLQ) ratio.

New Areas in Fuzzy Application 409

*<IP_ADDRESS>. Discrete control signal design*

controlled.

*2.3.3. Case study*

Simulink Sum Block is used to add the output of GFC to the output of FFC. Sampler system shown in Figure 24 is designed to convert the resultant analog control signal to discrete signal. The sampling time is selected to be greater than the system time constant. The discrete signal is connected to the input for the synchronized pulse generator to control the thyristor firing angle. Accordingly, the network var and the source power factor are

The Simulink circuit shown in Figure 24 is used as a base to study three 'test' cases. These cases are assigned to check the capability of the controller to operate within a considerable

 **Case 1 Case 2 Case 3 Source Voltage (Volt)<PHONE_NUMBER>**

Q1(kVAR) 3.197 142.857 142.857

Load Stages Q2 (kVAR) 12.789 571.429 428.571

F (Number of fuzzy membership functions) 7 7 7

**Table 2.** Circuit data and parameters for the test cases

Stage 1 P1 (kW) 2 40 40

Stage 2 P2 (kW) 12 200 240

Stage 3 P3(kW) 32 1200 1240 Q3 (kVAR) 22.38 1000 1000 Stage 4 P4 (kW) 42 2240 2280 Q4 (kVAR) 0 0 0 MLQ (MVAR) 0.02238 1 1

k1 1 1 1 k2 2.8 0.95 0.5 k3 0.44 1 2.3 k4 5.7 4.5 5.3 r1 0.64 0.52 0.95 r2 0.28 0.21 0.2 r3 0.03 0.02 0.04

ALPHA1 (DEGREE) 99 102 100 ALPHA2 (DEGREE) 124 113.5 105 ALPHA3 (DEGREE) 128 125.5 128 N (Number of sections) 4 4 4

variation of power factor values at different loading and voltage level.

## *<IP_ADDRESS>. Fuzzy Fine Controller (FFC) design*

Since the FGC output is not accurate due to the linearization process, FFC is designed to tune the firing g angle in order to achieve the desired power factor.

Two fuzzy control1ers are used along with automatic switching system to select the proper controller based on the power factor type (Lead or Lag), which is determined from the source VAR sign (+ or -) as shown in Figure 24.

The input of each FFC controller is the source power factor (PF) and the output is the corrective firing angle (M), which is required to fine tune the FGC output. Since the power factor value varies between zero and one, then the range of each FPC input membership function FFCMF(in)a is [0, 1].

The output membership function of each FFC is FFCMF(out-Lag)a for lagging input and FFCMF(out-Lead)a for leading input and, the magnitude of the firing angle range for the first linearzed section (1) is considered as a base to scale the FFC output as given hereinafter with the fact that any other section can be selected as the base. Also, [0, ∆1] and [-∆1 , 0] are assigned to FFCMF(out-Lag)b and FFCMF(out-Lead)b respectively. In addition to that, N multiplier factors (kN = ∆N*,/* 1*)* are used to scale the FFC output in order to match the magnitude of the firing angle range of the respective linearzed section (∆N). By this method, less number of FFC's are used.

Another Simulink automatic switching system SS2 as shown in Figure 24 synchronized with SS 1 is designed to check the active range of (QL/MLQ) in order to select the proper multiplier factor (KN). The fuzzy implication process functions used for FGC design are used for each FFC design as well, but with the following fuzzy rules:

$$\text{IF (PF)} \text{is } \text{FFCMF (in)}^{\text{a}} \text{THEN}(a\_{\Lambda}) \text{ is } \text{FFCMF (out } \text{Lag} \text{)}^{\text{a}} \tag{30}$$

$$\text{IF (PF)} \text{is } \text{FFCMF (in)}^{\text{a}} \text{THEN}(a\_{\Lambda}) \text{ is } \text{FFCMF (out } \text{Lead} \text{)}^{\text{a}} \tag{31}$$

### *<IP_ADDRESS>. Discrete control signal design*

Simulink Sum Block is used to add the output of GFC to the output of FFC. Sampler system shown in Figure 24 is designed to convert the resultant analog control signal to discrete signal. The sampling time is selected to be greater than the system time constant. The discrete signal is connected to the input for the synchronized pulse generator to control the thyristor firing angle. Accordingly, the network var and the source power factor are controlled.

## *2.3.3. Case study*

408 Fuzzy Controllers – Recent Advances in Theory and Applications

 Largest of Maximum for Defuzzification. Fuzzy rule : IF is MF(in)a THEN () is MF(out)b

a = 1,2, ... F (fuzzy membership function number) b = F-a+l(fuzzy membership function number) a and b are fuzzy membership function numbers F is the number of the fuzzy membership function.

*<IP_ADDRESS>. Fuzzy Fine Controller (FFC) design*

source VAR sign (+ or -) as shown in Figure 24.

function FFCMF(in)a is [0, 1].

method, less number of FFC's are used.

for each FFC design as well, but with the following fuzzy rules:

Simulink automatic switching system SS 1 as shown in Figure 24 is designed to check the active range of (QL/MLQ) in order to select the proper FGC based on that value of (QL/MLQ)

Since the FGC output is not accurate due to the linearization process, FFC is designed to

Two fuzzy control1ers are used along with automatic switching system to select the proper controller based on the power factor type (Lead or Lag), which is determined from the

The input of each FFC controller is the source power factor (PF) and the output is the corrective firing angle (M), which is required to fine tune the FGC output. Since the power factor value varies between zero and one, then the range of each FPC input membership

The output membership function of each FFC is FFCMF(out-Lag)a for lagging input and FFCMF(out-Lead)a for leading input and, the magnitude of the firing angle range for the first linearzed section (1) is considered as a base to scale the FFC output as given hereinafter with the fact that any other section can be selected as the base. Also, [0, ∆1] and [-∆1 , 0] are assigned to FFCMF(out-Lag)b and FFCMF(out-Lead)b respectively. In addition to that, N multiplier factors (kN = ∆N*,/* 1*)* are used to scale the FFC output in order to match the magnitude of the firing angle range of the respective linearzed section (∆N). By this

Another Simulink automatic switching system SS2 as shown in Figure 24 synchronized with SS 1 is designed to check the active range of (QL/MLQ) in order to select the proper multiplier factor (KN). The fuzzy implication process functions used for FGC design are used

> a a IF PF is FFCMF in THEN( ) is FFCMF out Lag

 a a IF PF is FFCMF in THEN( ) is FFCMF out Lead

(30)

(31)

tune the firing g angle in order to achieve the desired power factor.

Proportional for Aggregation.

where,

ratio.

The Simulink circuit shown in Figure 24 is used as a base to study three 'test' cases. These cases are assigned to check the capability of the controller to operate within a considerable variation of power factor values at different loading and voltage level.

**Table 2.** Circuit data and parameters for the test cases

## *2.3.4. Test cases data*

Each test case consists of four load stages. The load stages are selected such that the power factor varies from 0.3 to 1.0 lag. However, for practical cases, the power factor varies between 0.6 to 0.8 and thus the proposed wide range of power factor tested here is to demonstrate the capability of the designed controller. Table 2 summaries the circuit parameters for three voltage levels 120. 480 and 4160 Volts respectively. Figures 27-29 illustrate the linearization process for each case.

New Areas in Fuzzy Application 411

**Figure 29.** Linearization results for case 3

response time of the controller and its accuracy.

The controller is adjusted to correct the power factor of the test cases to a value greater than a desired value of 0.97. This value is the pre-set value and it can be any chosen practical value. Figures 30 -32 illustrate the results for cases 1-3 respectively. These figures show the variation of the load active and reactive power with respect to time for each case. The response of the controller during the test period represented by the firing angle is also shown. In addition to that, source and load power factor values are plotted to check the

*2.3.5. Results*

**Figure 27.** Linearization results for case1.

**Figure 28.** Linearization results for case 2

**Figure 29.** Linearization results for case 3

#### *2.3.5. Results*

410 Fuzzy Controllers – Recent Advances in Theory and Applications

illustrate the linearization process for each case.

**Figure 27.** Linearization results for case1.

**Figure 28.** Linearization results for case 2

Each test case consists of four load stages. The load stages are selected such that the power factor varies from 0.3 to 1.0 lag. However, for practical cases, the power factor varies between 0.6 to 0.8 and thus the proposed wide range of power factor tested here is to demonstrate the capability of the designed controller. Table 2 summaries the circuit parameters for three voltage levels 120. 480 and 4160 Volts respectively. Figures 27-29

*2.3.4. Test cases data*

The controller is adjusted to correct the power factor of the test cases to a value greater than a desired value of 0.97. This value is the pre-set value and it can be any chosen practical value. Figures 30 -32 illustrate the results for cases 1-3 respectively. These figures show the variation of the load active and reactive power with respect to time for each case. The response of the controller during the test period represented by the firing angle is also shown. In addition to that, source and load power factor values are plotted to check the response time of the controller and its accuracy.

New Areas in Fuzzy Application 413

Seconds

**Figure 31.** Results of test case 2.

**Figure 30.** Results of test case 1.

**Figure 31.** Results of test case 2.

**Figure 30.** Results of test case 1.

Seconds

New Areas in Fuzzy Application 415

The test results for the three cases show clearly how efficient is the controller. Even when the load reactive power is very small at both high and low power factor, the controller was successful in reaching an accurate level. During the stage where the load power factor is greater than 0.97 and, hence no need for capacitor compensation, the controller will check the source power factor at the beginning of that stage and if it drops below 0. 97, it will take an action in order to eliminate the compensation added in the previous stage. That is why the controller took an action as shown in Figure 30 for the fourth load stage of case no. 1 where the source power factor drops below 0.97. However, if the source power factor stays above the pre-set power factor value of 0.97 during the load stage where the load power factor is greater than 0.97, then no action will be taken as shown in Figure 32 for the fourth load stage of case no. 3. The time required for the controller to improve the power factor in all three cases is relatively short compared with practical applications. In real cases, the power factor does remain unchanged for relatively longer time. The maximum time for power factor correction was 0.35 second recorded in test case no.1. Overall, the graphs show that the controller works satisfactory under different load conditions and when there is no

As mentioned before, power factor correction is really an important issue. The designed controller presented in this section shows an efficient, fast and accurate technique in reactive power control. As seen from the overall structure of the controller, it is applicable for lagging power factor loads. Practically, this is almost true but not always where at rare occasions the power factor of the total load is leading not lagging. This will bring the attention towards generalizing the presented controller such that it will work for both cases. Several issues are also need to be considered in the future such as the dynamics of motors. As known that most connected loads are motors which really necessitate testing this controller under these circumstances. From the test results, it was seen that the speed of the controller depends on the system time constant and hence, a time delay is needed to assure that thc dynamics of the motors reach its equilibrium. Other issues related to Thyristors such as the harmonics are a1so need to be taken care by describing a harmonics filter. In addition to that, protection devices such as relays need to be checked during the controller action. Finally, the work presented was based on a single phase and it can be extended for

Most of the more advanced prediction techniques can be subdivided into two separate tasks. In a rst step, the modelling step, the algorithm uses a set of training data to identify a model of a process, from which the training data could have been obtained. In a second step, the simulation step, the algorithm uses the previously identied model to make predictions outside the training data set. The modelling algorithm can either attempt to identify the true structure of the system, from which the training data were obtained, or it can content itself with identifying any process able to explain the training data set. In the former case, we talk

*2.3.6. Discussion*

need for capacitor compensation.

three phase system.

**3. Trending and prediction**

**Figure 32.** Results of test case 3.

### *2.3.6. Discussion*

414 Fuzzy Controllers – Recent Advances in Theory and Applications

**Figure 32.** Results of test case 3.

Seconds

The test results for the three cases show clearly how efficient is the controller. Even when the load reactive power is very small at both high and low power factor, the controller was successful in reaching an accurate level. During the stage where the load power factor is greater than 0.97 and, hence no need for capacitor compensation, the controller will check the source power factor at the beginning of that stage and if it drops below 0. 97, it will take an action in order to eliminate the compensation added in the previous stage. That is why the controller took an action as shown in Figure 30 for the fourth load stage of case no. 1 where the source power factor drops below 0.97. However, if the source power factor stays above the pre-set power factor value of 0.97 during the load stage where the load power factor is greater than 0.97, then no action will be taken as shown in Figure 32 for the fourth load stage of case no. 3. The time required for the controller to improve the power factor in all three cases is relatively short compared with practical applications. In real cases, the power factor does remain unchanged for relatively longer time. The maximum time for power factor correction was 0.35 second recorded in test case no.1. Overall, the graphs show that the controller works satisfactory under different load conditions and when there is no need for capacitor compensation.

As mentioned before, power factor correction is really an important issue. The designed controller presented in this section shows an efficient, fast and accurate technique in reactive power control. As seen from the overall structure of the controller, it is applicable for lagging power factor loads. Practically, this is almost true but not always where at rare occasions the power factor of the total load is leading not lagging. This will bring the attention towards generalizing the presented controller such that it will work for both cases. Several issues are also need to be considered in the future such as the dynamics of motors. As known that most connected loads are motors which really necessitate testing this controller under these circumstances. From the test results, it was seen that the speed of the controller depends on the system time constant and hence, a time delay is needed to assure that thc dynamics of the motors reach its equilibrium. Other issues related to Thyristors such as the harmonics are a1so need to be taken care by describing a harmonics filter. In addition to that, protection devices such as relays need to be checked during the controller action. Finally, the work presented was based on a single phase and it can be extended for three phase system.

## **3. Trending and prediction**

Most of the more advanced prediction techniques can be subdivided into two separate tasks. In a rst step, the modelling step, the algorithm uses a set of training data to identify a model of a process, from which the training data could have been obtained. In a second step, the simulation step, the algorithm uses the previously identied model to make predictions outside the training data set. The modelling algorithm can either attempt to identify the true structure of the system, from which the training data were obtained, or it can content itself with identifying any process able to explain the training data set. In the former case, we talk

about a deep model, whereas models in the latter category are referred to as shallow models.

New Areas in Fuzzy Application 417

(34)

A = {x U|

Plausibility: The point y\* should represent B` from an intuitive point of view.

Continuity: A small Change in B` should not result in a large change in y\*.

represents the fuzzy set B`.

*3.1.1. Fuzzy C-Means Clustering Algorithm(FCM)*

to fuzzy membership μij, and cluster centroid Vi.

number of data points and f is a fuzziness index (greater than 1)

a. Initialize memberships μij of Xj belonging to cluster i such that

V =

d. Repeat steps 2 and 3 until the value of Jm is no longer decreasing.

μij , but different choices of initial μij might lead to local minima.

The FCM algorithm is executed by the following steps:

b. Compute the fuzzy centroid Vi from i=1 to i=c using

c. Update the fuzzy memberships μij using

Computational simplicity.

where X = [X1,….,Xn]t

*3.1.2. Cluster validity*

μA x }

c. Defuzzification is defined as a mapping from fuzzy set B` in V R to crisp point y\* V. Conceptually, the task of the defuzzification is to specify a point in V that best

The following three criteria should be considered in choosing the defuzification method:

The fuzzy c-means (FCM) clustering algorithm is the fuzzy equivalent of the nearest hard clustering algorithm [43,44], which minimizes the following objective function with respect

c n f 2

m ij j j i m ij j

(μ ) X

j i ij ( 2/(m 1)) j i j

The FCM always converges to strict local minimum of Jm starting from an initial guess of

The quality of a clustering is indicated by how closely the data points are associated to the cluster centers, and it is the membership functions, which measure the level of association or

(||(X ,V )||) μ (||(X ,V )||)

(μ )

( 2/(m 1))

m ij j i I j J ( μ ) || X ,V || (35)

ij j (μ ) 1 (36)

(37)

(38)

is a vector representing the data, c is the number of clusters, n is the

The identied model can be either a quantitative or a qualitative model. A quantitative model operates on the measurement data directly, whereas a qualitative model rst discretizes the measurement data, and then reasons about the discrete classes only. Also the model can be either a parametric model or a nonparametric model. A parametric model maps the knowledge contained in the training data set onto a set of model parameters. During the simulation phase, the training data are no longer needed, since the information contained in them is now stored in the parameter values. A non-parametric model only classies the training data during the modelling phase, and refers back to these classied training data during the simulation phase.

Fuzzy logic are used now a day in many application for diagnostic, prediction forecast and understanding the behaviour of very nonlinear systems such as marketing, electrical load forecast, work load analysis, technical analysis etc…

In this chapter Section 2.1, we shall introduce an important algorithm for classifying "clustering" the data based on fuzzy logic. Then two new fuzzy trending and prediction application shall be discussed. In Section 2.2 Accident rates Estimation Modelling Based on Human Factors shall be introduces, and in the next section 2.3, Fault Location in Distribution Networks shall be introduced.

## **3.1. Clustering algorithm and validity criteria**

Clustering attempts to assess the relationships among patterns of the data set by organizing the patterns into groups or clusters such that patterns within a cluster are more similar to each other than are patterns belonging to different clusters. Many algorithms for hard and fuzzy clustering have been developed to accomplish this[42]. An intimately related important issue is the cluster validity, which deals with the significance of the structure imposed by a clustering method [43].

For fuzzy sets, the following definitions are recalled from [36]:

Such that

$$\sum\_{i} \mu\_{i\mathbf{j}} = 1, \forall\_{\mathbf{j}'} 1 \le \mathbf{j} \le \mathbf{n}.\tag{33}$$

An α-cut of fuzzy set A is a crisp set Aα that contains all the elements in U that have membership value in A greater thanα, that is

$$\mathbf{A}\_{\boldsymbol{\alpha}} = \{ \mathbf{x} \in \mathbf{U} \, | \, \mu\_{\mathcal{A}} \left( \mathbf{x} \right) \ge a \} \tag{34}$$

c. Defuzzification is defined as a mapping from fuzzy set B` in V R to crisp point y\* V. Conceptually, the task of the defuzzification is to specify a point in V that best represents the fuzzy set B`.

The following three criteria should be considered in choosing the defuzification method:

416 Fuzzy Controllers – Recent Advances in Theory and Applications

training data during the simulation phase.

Distribution Networks shall be introduced.

imposed by a clustering method [43].

Such that

that takes values in the interval [0,1].

the set Mfc = {U Vcn ij [0,1]} (32)

membership value in A greater thanα, that is

forecast, work load analysis, technical analysis etc…

**3.1. Clustering algorithm and validity criteria**

For fuzzy sets, the following definitions are recalled from [36]:

models.

about a deep model, whereas models in the latter category are referred to as shallow

The identied model can be either a quantitative or a qualitative model. A quantitative model operates on the measurement data directly, whereas a qualitative model rst discretizes the measurement data, and then reasons about the discrete classes only. Also the model can be either a parametric model or a nonparametric model. A parametric model maps the knowledge contained in the training data set onto a set of model parameters. During the simulation phase, the training data are no longer needed, since the information contained in them is now stored in the parameter values. A non-parametric model only classies the training data during the modelling phase, and refers back to these classied

Fuzzy logic are used now a day in many application for diagnostic, prediction forecast and understanding the behaviour of very nonlinear systems such as marketing, electrical load

In this chapter Section 2.1, we shall introduce an important algorithm for classifying "clustering" the data based on fuzzy logic. Then two new fuzzy trending and prediction application shall be discussed. In Section 2.2 Accident rates Estimation Modelling Based on Human Factors shall be introduces, and in the next section 2.3, Fault Location in

Clustering attempts to assess the relationships among patterns of the data set by organizing the patterns into groups or clusters such that patterns within a cluster are more similar to each other than are patterns belonging to different clusters. Many algorithms for hard and fuzzy clustering have been developed to accomplish this[42]. An intimately related important issue is the cluster validity, which deals with the significance of the structure

a. A fuzzy set in a universe discourse U is characterized by a membership function A(x)

b. Let X={x1,,xn} be any set, Vcn be the set of real cn matrices U=[ij], c, i, j be integer numbers with 2 c n , 1 i c and 1 j n .Then the fuzzy partition matrix for X is

An α-cut of fuzzy set A is a crisp set Aα that contains all the elements in U that have

μij j 1, ,1 j n. *i* (33)

Continuity: A small Change in B` should not result in a large change in y\*.

#### *3.1.1. Fuzzy C-Means Clustering Algorithm(FCM)*

The fuzzy c-means (FCM) clustering algorithm is the fuzzy equivalent of the nearest hard clustering algorithm [43,44], which minimizes the following objective function with respect to fuzzy membership μij, and cluster centroid Vi.

$$J\_{\mathbf{m}} = \sum\_{\mathbf{l}} \, ^c\_{\mathbf{l}} \sum\_{\mathbf{j}} ^n (\mu\_{\mathbf{ij}})^\mathbf{f} \mid \mathbb{I} \left( \mathbf{X}\_{\mathbf{j}'} \mathbf{V}\_{\mathbf{i}} \right) \mid \mathbb{I}^2 \tag{35}$$

where X = [X1,….,Xn]t is a vector representing the data, c is the number of clusters, n is the number of data points and f is a fuzziness index (greater than 1)

The FCM algorithm is executed by the following steps:

a. Initialize memberships μij of Xj belonging to cluster i such that

$$\sum\_{\mathbf{j}} (\mu\_{\mathbf{i}\mathbf{j}}) = 1 \tag{36}$$

b. Compute the fuzzy centroid Vi from i=1 to i=c using

$$\mathbf{V}\_{\mathbf{i}} = \frac{\sum\_{\mathbf{j}} (\mu\_{\tilde{\mathbf{i}}\tilde{\mathbf{j}}})^{\mathbf{m}} \times \mathbf{X}\_{\mathbf{j}}}{\sum\_{\mathbf{j}} (\mu\_{\tilde{\mathbf{i}}\tilde{\mathbf{j}}})^{\mathbf{m}}} \tag{37}$$

c. Update the fuzzy memberships μij using

$$\mu\_{\rm ij} = \frac{\left( \mid \mid (\mathbf{X}\_{\mathbf{j}}, \mathbf{V}\_{\mathbf{i}}) \mid \mid \right)^{\left( -2/(\mathfrak{m}-1) \right)}}{\sum\_{\mathbf{j}} \left( \mid \mid (\mathbf{X}\_{\mathbf{j}}, \mathbf{V}\_{\mathbf{i}}) \mid \mid \right)^{\left( -2/(\mathfrak{m}-1) \right)}} \tag{38}$$

d. Repeat steps 2 and 3 until the value of Jm is no longer decreasing.

The FCM always converges to strict local minimum of Jm starting from an initial guess of μij , but different choices of initial μij might lead to local minima.

#### *3.1.2. Cluster validity*

The quality of a clustering is indicated by how closely the data points are associated to the cluster centers, and it is the membership functions, which measure the level of association or

classification. If the value of one of the membership is significantly larger than the others for a particular data point, then that point is identified as being a part of the subset of the data represented by the corresponding cluster center. But, each data point has c memberships; so, it is desirable to summarize the information contained in the memberships by a single number, which indicates how well the data point is classified by the clustering. This can be done in a variety of ways; for example, for the data point Xj with memberships {μij, ...,μcj}, one could use any of the following:

$$\text{Index}\mathbf{1} = \sum\_{\mathbf{i}} (\mu\_{\mathbf{i}\mathbf{j}})^2 \tag{39}$$

New Areas in Fuzzy Application 419

Policies. Line management is encouraged to conduct regular inspections of the workplace and employees are trained to behave safely and are given the appropriate protective equipment. The impact of such initiatives could be seen in the overall downward trend in accident statistics from 1990 to 1998/99. After 2000, accident statistics started rising in many UK industrial sectors. In the Quarry Industry, for example, there has been a 60% rise in the

Another study by Stanford University [39] indicated that risk taking is often a normal part of human psychology. We sometimes drive too fast or take chances we should not. Risk which is taken on the job site, however, can be fatal. This study found that many workers believe that taking unnecessary risks is an accepted part of the job process. This risk acceptance attitude leads to carelessness and accidents. The results of the Stanford study show that workers who are likely to have lost-time accidents share similar characteristics. These workers have a negative attitude toward doing their jobs safely, and they accept unnecessary risk and, therefore, do not work safely. Taking unnecessary risks and adopting a poor safety attitude simply makes workers more prone to accident occurrences. The conclusion of the Stanford study clearly supports the contention that employee actions and attitudes can affect the number and type of workplace accidents. Employers can and do address this attitude of risk taking through safety education, safety rules, and training programs. However, no employer can supervise each employee every minute of the work

In [40], the paper focuses on the development and representation of linguistic variables to model risk levels subjectively. These variables are then quantified using fuzzy set theory. In this paper the development of two safety evaluation frameworks, using fuzzy logic approaches for maritime engineering safety based decision support in the concept design stage are presented. An example is used to illustrate and compare the proposed approaches. The paper also suggests that future risk analysis in maritime engineering applications may

The field of fuzzy systems has been making rapid progress over the past decade [36]. There

a. The problem is too complicated for precise description to be obtained, therefore approximation, or fuzziness, must be introduced in order to be a reasonable and net

b. As we move into the information era, knowledge is becoming increasingly important and the need for a theory to formulate human knowledge in a systematic manner

But as a general principle, a good engineering theory should be capable of making use of all the available information effectively. For many practical systems, important information comes from two sources: one source is from human experts who describe their knowledge about the system in natural languages; the other is sensory measurements or mathematical

take full advantages of fuzzy logic approaches to complement existing ones.

are two kinds of justification for fuzzy systems to be used to achieve our objective:

number of fatalities.

day.

traceable model.

becomes the norm not the exception.

models that are derived from to physical laws.

$$\text{Index2} = \sum\_{\text{i}} \mu\_{\text{i}\text{j}} \log(\mu\_{\text{i}\text{j}}) \tag{40}$$

$$\text{Index3=max}\_{i}(\mu\_{ij})\tag{41}$$

$$\mathsf{Index}4 = \min\_{\mathsf{i}} \left( \mu\_{\mathsf{i}\mathsf{j}} \right) / \mathsf{max}\_{\mathsf{i}} (\mu\_{\mathsf{i}\mathsf{j}}) \tag{42}$$

In fact, theses four indices of these are used as measure of the quality of clustering and are the basis for the *validity functional*, *partition coefficient*, *classification entropy, and proportion exponent*, respectively.

To illustrate the use of validity functional, we shall focus on the partition coefficient technique because of its simplicity. It is based on using Sj =Σi (μij)2 as a measure of how well the jth data point has been classified. This is a reasonable indicator because the closer a data point is to a cluster center, the closer Sj is to 1, the maximum value it could have. Conversely, the further away the kth point is from all the cluster centers the closer the value of Sj is to 1/c, the minimum possible value. The partition coefficient is then the average over the data set of the Sj's. In particular, for a data set X={x1,…,xi} and a specific choice of c and m one obtains the output of fuzzy c-means and computes the partition coefficient (PC) by PC=Σj(Σi(μij)2/n. The closer this value is to one the better the data are classified. So, in theory, one computes PC for the outputs of a variety of values of c and m selects the best clustering as the one corresponding to the highest partition coefficient [44,45].

## **3.2. Accident rates estimation modeling based on human factors using fuzzy cmean clustering techniques**

Several individual books [37] and projects shed light on worker accident causation. One study on the Bonneville Dam project, reported that seven times the number of work accidents that had occurred on this project were due to unsafe employee actions rather than to unsafe site conditions. In addition, this study found that the negative attitude of the workers toward safety was a major factor in accident occurrence.

In [38]. Many organizations spend a lot of time and effort trying to improve safety. As well as addressing technical and hardware issues, many conduct safety management system audits to discover deviations from the performance standards set in their Health & Safety Policies. Line management is encouraged to conduct regular inspections of the workplace and employees are trained to behave safely and are given the appropriate protective equipment. The impact of such initiatives could be seen in the overall downward trend in accident statistics from 1990 to 1998/99. After 2000, accident statistics started rising in many UK industrial sectors. In the Quarry Industry, for example, there has been a 60% rise in the number of fatalities.

418 Fuzzy Controllers – Recent Advances in Theory and Applications

corresponding to the highest partition coefficient [44,45].

workers toward safety was a major factor in accident occurrence.

one could use any of the following:

*exponent*, respectively.

**mean clustering techniques**

classification. If the value of one of the membership is significantly larger than the others for a particular data point, then that point is identified as being a part of the subset of the data represented by the corresponding cluster center. But, each data point has c memberships; so, it is desirable to summarize the information contained in the memberships by a single number, which indicates how well the data point is classified by the clustering. This can be done in a variety of ways; for example, for the data point Xj with memberships {μij, ...,μcj},

In fact, theses four indices of these are used as measure of the quality of clustering and are the basis for the *validity functional*, *partition coefficient*, *classification entropy, and proportion*

To illustrate the use of validity functional, we shall focus on the partition coefficient technique because of its simplicity. It is based on using Sj =Σi (μij)2 as a measure of how well the jth data point has been classified. This is a reasonable indicator because the closer a data point is to a cluster center, the closer Sj is to 1, the maximum value it could have. Conversely, the further away the kth point is from all the cluster centers the closer the value of Sj is to 1/c, the minimum possible value. The partition coefficient is then the average over the data set of the Sj's. In particular, for a data set X={x1,…,xi} and a specific choice of c and m one obtains the output of fuzzy c-means and computes the partition coefficient (PC) by PC=Σj(Σi(μij)2/n. The closer this value is to one the better the data are classified. So, in theory, one computes PC for the outputs of a variety of values of c and m selects the best clustering as the one

**3.2. Accident rates estimation modeling based on human factors using fuzzy c-**

Several individual books [37] and projects shed light on worker accident causation. One study on the Bonneville Dam project, reported that seven times the number of work accidents that had occurred on this project were due to unsafe employee actions rather than to unsafe site conditions. In addition, this study found that the negative attitude of the

In [38]. Many organizations spend a lot of time and effort trying to improve safety. As well as addressing technical and hardware issues, many conduct safety management system audits to discover deviations from the performance standards set in their Health & Safety

2

ij i Index1= ( μ ) (39)

ij ij i Index2 = μ log(μ ) (40)

i ij Index3= max (μ ) (41)

i ij i ij Index4= min (μ ) / max (μ ) (42)

Another study by Stanford University [39] indicated that risk taking is often a normal part of human psychology. We sometimes drive too fast or take chances we should not. Risk which is taken on the job site, however, can be fatal. This study found that many workers believe that taking unnecessary risks is an accepted part of the job process. This risk acceptance attitude leads to carelessness and accidents. The results of the Stanford study show that workers who are likely to have lost-time accidents share similar characteristics. These workers have a negative attitude toward doing their jobs safely, and they accept unnecessary risk and, therefore, do not work safely. Taking unnecessary risks and adopting a poor safety attitude simply makes workers more prone to accident occurrences. The conclusion of the Stanford study clearly supports the contention that employee actions and attitudes can affect the number and type of workplace accidents. Employers can and do address this attitude of risk taking through safety education, safety rules, and training programs. However, no employer can supervise each employee every minute of the work day.

In [40], the paper focuses on the development and representation of linguistic variables to model risk levels subjectively. These variables are then quantified using fuzzy set theory. In this paper the development of two safety evaluation frameworks, using fuzzy logic approaches for maritime engineering safety based decision support in the concept design stage are presented. An example is used to illustrate and compare the proposed approaches. The paper also suggests that future risk analysis in maritime engineering applications may take full advantages of fuzzy logic approaches to complement existing ones.

The field of fuzzy systems has been making rapid progress over the past decade [36]. There are two kinds of justification for fuzzy systems to be used to achieve our objective:

But as a general principle, a good engineering theory should be capable of making use of all the available information effectively. For many practical systems, important information comes from two sources: one source is from human experts who describe their knowledge about the system in natural languages; the other is sensory measurements or mathematical models that are derived from to physical laws.

An important task, therefore, is to combine these two types of information into system designs. Therefore, the key question is how to transform human knowledge base into a mathematical formula or model. Essentially, what a fuzzy system does is to perform this transformation in a systematic way.

New Areas in Fuzzy Application 421

*<IP_ADDRESS>. Questionnaire design*

same field made the same judgments.

linguistic evaluations.

work.

*<IP_ADDRESS>. The response*

*<IP_ADDRESS>. Limitations*

only to one sex.

In the design of questionnaire (Appendix) we have selected some certain features of human nature that, we believe, have a great potential on the accident causation in the local market [41]. The first page of the questionnaire concentrated on the personal information of the workers: i.e. 'height', 'weight', 'optical status', 'hearing ability', 'general health', 'education' and 'adherence to safety rules'. In this part we used some linguistic evaluations like 'high', 'low', 'fair', 'good', 'medium' etc., and in some others we have used numerical evaluations like in height, weight, as well as education. Since the objective here is to create a fuzzy model, high accuracy is not important and we considered that the respondents from the

The second page was designed to concentrate on the work information: 'overtime work', 'experience', 'work nature', 'work type', 'hazardous level', 'needs for safety-gears', 'work location' and 'level of boredom'. Again, we have used some linguistic evaluations as well as numerical evaluations. We also concentrated on the managerial factors: 'salary received on time', 'level of training' 'importance' and level of importance placed on safety', with only

The third page was a mix of both external factors: 'noise', 'live with family' and 'communication 'language', and accident history focusing on the number of accidents the laborer has faced during his work in the local construction market, which was the most important data that we needed to develop our model. The severity of accidents has not been

The cases obtained for this study were collected from three different construction companies selected to represent the local market. We have tried to select the cases from different ranks of the workforce, from higher levels to the lower levels to be able to study the different accident level cases that serve the purpose as well as adding versatility and diversity to this

From the original cases that we collected on 95 people, we included only 76 cases and excluded 19 cases, which were incomplete. This has produced a very high response rate that reached 82.1%, which is relatively high, especially as the questionnaire is lengthy and a little bit complicated. More cases could have been obtained. However, since the aim of this study was to develop a model for accidents, the number of cases is not set a priority. Rather, the cases are accumulated and the algorithms are stopped when a cluster validity criterion is

Limitations in this study should be noted. One of the limitations is that we did not include any specific information related to the accident consequences. Another limitation is that this study was made only on males and no female cases were studied, which makes it specific

taken into account since it is not considered a factor influencing the accident rate.

satisfied thereby yielding the optimum number of clusters.

In this part of Section 2, we attempt to use a completely different approach to analyze accidents. A model shall be developed for data collected from an accident rate questionnaire filled-in by laborers working for a reputable construction company. This questionnaire was designed to include information about human factors, as well as other factors such as work type, managerial factors, training, physical factors and the historical accident rate for each labor during his period of employment in this particular construction company, and his experience during his career life time. The collected data shall be split into a training set for model construction and a test for model verification. The training information shall be classified into a number of groups or clusters, the centroids of these clusters were subsequently used to generate a set of rules to develop a fuzzy engine, which can then predict and forecast the rate of accidents. The test cases shall be used to verify and validate the developed model. Discussion on the results obtained from using fuzzy logic techniques shall be carried out.

## *3.2.1. Data organization*

Construction sites are very dynamic and complex by nature, creating the potential for hazards that change constantly. So, what was safe yesterday may no longer be safe today. Thus, safety precautions should be followed and controlled. Unsafe working conditions and accidents are usually warning signs that something is wrong and has to be rectified.

Different government authorities measure safety at construction sites [41], however coordination and sharing of information with each other is still lacking. In addition, the data available on construction site accidents are neither accurate nor complete, due to the absence of a reliable accident reporting and recording system. Incomplete records are due to the poor accident investigation that may be a result of:

For these reasons we endeavored to avoid the normal way of doing the job, and instead, we focused on the laborer, himself, and his accident rates during his years of work experience as an expert source of data. We have tried to design the questionnaire in such a way as to serve our purpose of analyzing the data, and selected the interview method to get the maximum precise data possible.

## *<IP_ADDRESS>. Questionnaire design*

420 Fuzzy Controllers – Recent Advances in Theory and Applications

transformation in a systematic way.

shall be carried out.

*3.2.1. Data organization*

 Narrow interpretation Judgmental behavior.

maximum precise data possible.

poor accident investigation that may be a result of:

 Reluctance of reporters to assert authority. Inexperienced and untrained investigators

 Incomplete or erroneous conclusions. Delays in accident investigations.

An important task, therefore, is to combine these two types of information into system designs. Therefore, the key question is how to transform human knowledge base into a mathematical formula or model. Essentially, what a fuzzy system does is to perform this

In this part of Section 2, we attempt to use a completely different approach to analyze accidents. A model shall be developed for data collected from an accident rate questionnaire filled-in by laborers working for a reputable construction company. This questionnaire was designed to include information about human factors, as well as other factors such as work type, managerial factors, training, physical factors and the historical accident rate for each labor during his period of employment in this particular construction company, and his experience during his career life time. The collected data shall be split into a training set for model construction and a test for model verification. The training information shall be classified into a number of groups or clusters, the centroids of these clusters were subsequently used to generate a set of rules to develop a fuzzy engine, which can then predict and forecast the rate of accidents. The test cases shall be used to verify and validate the developed model. Discussion on the results obtained from using fuzzy logic techniques

Construction sites are very dynamic and complex by nature, creating the potential for hazards that change constantly. So, what was safe yesterday may no longer be safe today. Thus, safety precautions should be followed and controlled. Unsafe working conditions and

Different government authorities measure safety at construction sites [41], however coordination and sharing of information with each other is still lacking. In addition, the data available on construction site accidents are neither accurate nor complete, due to the absence of a reliable accident reporting and recording system. Incomplete records are due to the

For these reasons we endeavored to avoid the normal way of doing the job, and instead, we focused on the laborer, himself, and his accident rates during his years of work experience as an expert source of data. We have tried to design the questionnaire in such a way as to serve our purpose of analyzing the data, and selected the interview method to get the

accidents are usually warning signs that something is wrong and has to be rectified.

In the design of questionnaire (Appendix) we have selected some certain features of human nature that, we believe, have a great potential on the accident causation in the local market [41]. The first page of the questionnaire concentrated on the personal information of the workers: i.e. 'height', 'weight', 'optical status', 'hearing ability', 'general health', 'education' and 'adherence to safety rules'. In this part we used some linguistic evaluations like 'high', 'low', 'fair', 'good', 'medium' etc., and in some others we have used numerical evaluations like in height, weight, as well as education. Since the objective here is to create a fuzzy model, high accuracy is not important and we considered that the respondents from the same field made the same judgments.

The second page was designed to concentrate on the work information: 'overtime work', 'experience', 'work nature', 'work type', 'hazardous level', 'needs for safety-gears', 'work location' and 'level of boredom'. Again, we have used some linguistic evaluations as well as numerical evaluations. We also concentrated on the managerial factors: 'salary received on time', 'level of training' 'importance' and level of importance placed on safety', with only linguistic evaluations.

The third page was a mix of both external factors: 'noise', 'live with family' and 'communication 'language', and accident history focusing on the number of accidents the laborer has faced during his work in the local construction market, which was the most important data that we needed to develop our model. The severity of accidents has not been taken into account since it is not considered a factor influencing the accident rate.

The cases obtained for this study were collected from three different construction companies selected to represent the local market. We have tried to select the cases from different ranks of the workforce, from higher levels to the lower levels to be able to study the different accident level cases that serve the purpose as well as adding versatility and diversity to this work.

#### *<IP_ADDRESS>. The response*

From the original cases that we collected on 95 people, we included only 76 cases and excluded 19 cases, which were incomplete. This has produced a very high response rate that reached 82.1%, which is relatively high, especially as the questionnaire is lengthy and a little bit complicated. More cases could have been obtained. However, since the aim of this study was to develop a model for accidents, the number of cases is not set a priority. Rather, the cases are accumulated and the algorithms are stopped when a cluster validity criterion is satisfied thereby yielding the optimum number of clusters.

#### *<IP_ADDRESS>. Limitations*

Limitations in this study should be noted. One of the limitations is that we did not include any specific information related to the accident consequences. Another limitation is that this study was made only on males and no female cases were studied, which makes it specific only to one sex.

#### *<IP_ADDRESS>. The feature matrix*

The feature matrix (FM) is the most important part of our work in this section , since by using this matrix we have been able to convert the linguistic variables into numerical variables. Thus, we can deal easily with practical cases and reduce the required operations of processing the output.

New Areas in Fuzzy Application 423

In order to obtain the models, the following steps have been implemented:

MATLAB fuzzy toolbox has been used to implement the above three steps [46].

between the number of clusters and the corresponding error where:

**Figure 33.** Relation between number of clusters and the corresponding error

column.

the normalized data.

number of clusters.

*<IP_ADDRESS>. Scaling of data*

the following formula:

Xi = Vector element Xmin = Min. Vector element Xmax = Max. Vector element

Where:

*<IP_ADDRESS>. Clustering results and discussions*

The results can be summarized as follows:

 Each column in the feature matrix is normalized by dividing all the numbers in this column by the maximum number of the absolute values of all the numbers in the said

Cluster validity study is implemented to determine the optimum number of clusters for

FCM technique is implemented to determine the centroids matrix of the selected

a. The optimum number of cluster (twelve) is determined by implementing cluster validity technique. The result is illustrated in Figure 33 which gives the relation

b. Table 9 (Appendix) illustrates the centroide matrix for the optimum number of clusters (twelve), where FW(j) stands for feature weight described in Table 10 (Appendix).

In order to reduce the error in the estimation of the membership functions ranges (e.g. the estimated range of FW5 was from 1 to 5, while from centroide matrix the range for the same feature is found to be from 3.750 to 4.9953) each feature vector element is scaled according to

X X X /X X scaled i min max min (43)

The columns of the FM matrix represent the feature variables which we obtained from the questionnaire. The rows of the matrix represent the different cases of laborers that we selected for interviews. Thus, for each case in the matrix we mapped the linguistic meanings into numbers according to weights we have proposed. For example, in case labelled (S1), the feature weight 1 (FW1) represents the rate of accidents per year of experience, which is the actual representation of the accident rates. FW2 represents the ratio between weight and height (specific weight). FW3 represents optical status and FW4 represents hearing ability. These feature weights are scaled on a scale of five from 1 to 5, to represent the linguistic variables. Therefore, 1 means `very bad`, 2 means `bad`, 3 means `medium`, 4 means `good` and 5 means `very good`. All the other variables are dealt in the same way until the matrix was generated. A sample of the feature matrix is shown in Table3.

**Table 3.** Sample of feature matrix illustrates the weight of some features for the first 10 cases

## *3.2.2. Modelling*

In this stage, FCM techniques will be implemented on the feature matrix after normalizing the data, based on column maximum values for ease and as being more indicative, then deciding the optimum number of clusters, by applying cluster validity techniques. The centroides for these optimum clusters are considered perfect models represent the feature matrix.

In order to obtain the models, the following steps have been implemented:

MATLAB fuzzy toolbox has been used to implement the above three steps [46].

#### *<IP_ADDRESS>. Clustering results and discussions*

422 Fuzzy Controllers – Recent Advances in Theory and Applications

The feature matrix (FM) is the most important part of our work in this section , since by using this matrix we have been able to convert the linguistic variables into numerical variables. Thus, we can deal easily with practical cases and reduce the required operations

The columns of the FM matrix represent the feature variables which we obtained from the questionnaire. The rows of the matrix represent the different cases of laborers that we selected for interviews. Thus, for each case in the matrix we mapped the linguistic meanings into numbers according to weights we have proposed. For example, in case labelled (S1), the feature weight 1 (FW1) represents the rate of accidents per year of experience, which is the actual representation of the accident rates. FW2 represents the ratio between weight and height (specific weight). FW3 represents optical status and FW4 represents hearing ability. These feature weights are scaled on a scale of five from 1 to 5, to represent the linguistic variables. Therefore, 1 means `very bad`, 2 means `bad`, 3 means `medium`, 4 means `good` and 5 means `very good`. All the other variables are dealt in the same way until the matrix

was generated. A sample of the feature matrix is shown in Table3.

**Weight/ Height**

**Table 3.** Sample of feature matrix illustrates the weight of some features for the first 10 cases

In this stage, FCM techniques will be implemented on the feature matrix after normalizing the data, based on column maximum values for ease and as being more indicative, then deciding the optimum number of clusters, by applying cluster validity techniques. The centroides for these optimum clusters are considered perfect models represent the feature

**Optical status**

S1 1/12 71/160 6/18 5 5 4 5 2 S2 5/12 77/ 170 6/60 5 4 3 5 3 S3 1/21 90/175 6/6 4 5 5 5 5 S4 0/8 54/165 6/60 4 5 2 5 20 S5 6/3.5 68/187 6/36 3 4 3 5 0 S6 10/11 85/177 6/6 4 4 3 5 10 S7 2/19 76/173 6/60 4 5 5 5 14 S8 18/25 72/170 6/18 3 4 4 4 4 S9 45/14 80/176 6/6 4 4 3 4 8 S10 3/20 81/174 6/6 4 4 5 4 1

**Hearing Ability**

**General Health**

**Adherence**

**to Safety**

**Education**

**Overtime work**

*<IP_ADDRESS>. The feature matrix*

of processing the output.

**Feature**

**Accidents/**

**experience**

**Case No**

*3.2.2. Modelling*

matrix.

The results can be summarized as follows:

a. The optimum number of cluster (twelve) is determined by implementing cluster validity technique. The result is illustrated in Figure 33 which gives the relation between the number of clusters and the corresponding error where:

**Figure 33.** Relation between number of clusters and the corresponding error

b. Table 9 (Appendix) illustrates the centroide matrix for the optimum number of clusters (twelve), where FW(j) stands for feature weight described in Table 10 (Appendix).

#### *<IP_ADDRESS>. Scaling of data*

In order to reduce the error in the estimation of the membership functions ranges (e.g. the estimated range of FW5 was from 1 to 5, while from centroide matrix the range for the same feature is found to be from 3.750 to 4.9953) each feature vector element is scaled according to the following formula:

$$\mathbf{X}\_{\text{scaled}} = \mathbf{X}\_{\text{i}} - \mathbf{X}\_{\text{min}} \left( \mathbf{X}\_{\text{max}} - \mathbf{X}\_{\text{min}} \right) \tag{43}$$

Where:

Xi = Vector element Xmin = Min. Vector element Xmax = Max. Vector element

### *<IP_ADDRESS>. Model development*

By comparing each row of the centroids matrix with the contents of the questionnaire after scaling, one can infer the structure of the respective model, for example maximum and minimum rate of accident as illustrated in Table 11 (Appendix). From structure, following features can be extracted:

New Areas in Fuzzy Application 425

**Figure 34.** Fuzzy Accident Prediction Flowchart

results are shown as follows in Table 4:

**Table 4.** Standard deviation results for the test cases

results are given in hereinafter.

*<IP_ADDRESS>. Relevant results*

accident rate.

Mamdani inference engine with centroid defuzzification and proportional aggregation is used in the construction of the fuzzy system. MATLAB Fuzzy Toolbox has been used to implement the required fuzzy prediction system. Eight cases have been tested and the

a. Eight additional random cases have been chosen from a local construction company. Features, actual accident rate and output of fuzzy prediction system for each test case are given in Table 12 (Appendix). The standard deviation between the actual accident rate and the predicted accident rate for each test case is calculated then the average standard deviation is obtained in order to determine the validity of the model. The

> **Case Number Standard Deviation** Case#1 0.3536 Case#2 0.2828 Case#3 0.4243 Case#4 0.7071 Case#5 01414 Case#6 0.1414 Case#7 0.2121 Case#8 0.0318 Average 0.2868

b. In (Appendix), Figure 37 shows the output of the prediction system for a case that fires all the linguistic values at the middle. It is observed from this result that the laborers with `average` personal information, `average` work condition, `average` managerial condition and `average` external effects are exposed to `average` annual

## *3.2.3. Accident rate prediction*

Now, fuzzy logic techniques will be implemented using the models obtained in Section 2.2.2, as perfect fuzzy rules, to predict the accident rate for any laborer who works in the construction field. The following flow chart (Figure 34) describes the fuzzy accident prediction system:

The models obtained from the fuzzy c -means the clustering process has been considered as very good and suitable fuzzy rules that govern the relation between the laborers in construction field and the expected annual rate of accident.

The beauty of using this type of clustering is not only to achieve the required models, but also these models are fuzzy, and can be geared in the fuzzy engine.

In order to fuzzify the model variables, a suitable number of Gauss functions (linguistic variable) is selected for each linguistic value so that any rule must fire all the linguistic values and the rules are given as follows [36]:

$$\text{IF } \left( \text{FW1 is } \text{mf}\_{\text{a}}^{-1} \right) \text{ and } \left( \text{FW2 is } \text{mf}\_{\text{b}}^{-2} \right) \text{ and } \dots \left( \text{FW2} \text{ is } \text{mf}\_{\text{v}}^{-2} \right) \text{ THEN } \left( \text{FW2} \text{ is } \text{mf}\_{\text{w}}^{-2} \right) \tag{44}$$

Where,

FW(1to22) : input linguistic variable FW23 : output linguistic variable mf a,b ..w : semantic rule a,b ..w: integer number from 1 to 5

**Figure 34.** Fuzzy Accident Prediction Flowchart

Mamdani inference engine with centroid defuzzification and proportional aggregation is used in the construction of the fuzzy system. MATLAB Fuzzy Toolbox has been used to implement the required fuzzy prediction system. Eight cases have been tested and the results are given in hereinafter.

#### *<IP_ADDRESS>. Relevant results*

424 Fuzzy Controllers – Recent Advances in Theory and Applications

By comparing each row of the centroids matrix with the contents of the questionnaire after scaling, one can infer the structure of the respective model, for example maximum and minimum rate of accident as illustrated in Table 11 (Appendix). From structure, following

b. The expected range of accident rate in company varies from 0.1581 to 2.8894 accidents

c. By comparing the above two extreme models, one can easily discover that receiving

d. The highest rate of accidents occurs to non-local workers (Live without their family

e. The twelve models obtained above shall be utilized later as fuzzy rules representing the this local construction company to predict the rate of accident for any person working

Now, fuzzy logic techniques will be implemented using the models obtained in Section 2.2.2, as perfect fuzzy rules, to predict the accident rate for any laborer who works in the construction field. The following flow chart (Figure 34) describes the fuzzy accident

The models obtained from the fuzzy c -means the clustering process has been considered as very good and suitable fuzzy rules that govern the relation between the laborers in

The beauty of using this type of clustering is not only to achieve the required models, but

In order to fuzzify the model variables, a suitable number of Gauss functions (linguistic variable) is selected for each linguistic value so that any rule must fire all the linguistic

 1 2 22 23 ab v w IF FW1 is mf and FW2 is mf and . . . FW22 is mf THEN FW23 is mf (44)

a. Labours in this construction company can be classified into 12 models.

salary on time has no effective impact on accident rate.

construction field and the expected annual rate of accident.

values and the rules are given as follows [36]:

FW(1to22) : input linguistic variable FW23 : output linguistic variable

a,b ..w: integer number from 1 to 5

mf a,b ..w : semantic rule

also these models are fuzzy, and can be geared in the fuzzy engine.

*<IP_ADDRESS>. Model development*

features can be extracted:

most of the time).

on construction field.

*3.2.3. Accident rate prediction*

per year.

prediction system:

Where,

a. Eight additional random cases have been chosen from a local construction company. Features, actual accident rate and output of fuzzy prediction system for each test case are given in Table 12 (Appendix). The standard deviation between the actual accident rate and the predicted accident rate for each test case is calculated then the average standard deviation is obtained in order to determine the validity of the model. The results are shown as follows in Table 4:

b. In (Appendix), Figure 37 shows the output of the prediction system for a case that fires all the linguistic values at the middle. It is observed from this result that the laborers with `average` personal information, `average` work condition, `average` managerial condition and `average` external effects are exposed to `average` annual accident rate.

c. Figure 38(Appendix) correlates the laborers general health and specific weight with their annual accident rate considering all other factors are `average`. It is clear that the annual accident rate increases with the significant increase of the specific weight and significant poorness of the general health.

New Areas in Fuzzy Application 427

**3.3. Fault location in distribution networks using fuzzy c-mean clustering**

In last section of this chapter we shall studies an existing 13.8 kilovolt distribution network which, serves an oil production field spread over an area of approximately sixty kilometers square, in order to locate any fault that may occur anywhere in the network using fuzzy c-

In addition, we shall introduce several methods for normalizing data and selecting the optimum number of clusters in order to classify data. Results and conclusion shall be also

A joint venture oil company possesses two production areas, Area1 and Area2. For each area a power distribution system is provided. The power supply required for the two fields is provided from 130 MVA capacity power generation plant located in Area1. Twenty MW is transmitted to Area2 via 35KM long over head transmission line (OHTL) on wooden

Area2 distribution system consists of three radial overhead transmission lines, 8-10 KM long each, serve submersible oil pumps and other loads scattered in the field (60-kilometer square). These overhead transmission lines have neither differential relays nor sectionalizing fuses. Programmable logic Controller (PLC) is connected to the incomers, outgoings and generators auxiliaries at Area2 power station. This PLC records the running

Due to the aging of the system, remote area problems and harsh desert weather, repeated faults are experienced in the grid. Because of unavailability of differential relays and/or sectionalizing fuses, it is very difficult and long time is consumed to locate any fault in this network. The problem is reflecting passively on the oil productivity of this important area.

In order to measure the features of the faults at 176 nodes of the network; load flow study is implemented to determine the respective power loss for each short circuit case and also short circuit study is carried out to determine the feature vector for each short circuit case.

The results, obtained from the load flow study and the short circuit study, shall be used to form the network *feature matrix,* which will be clustered and analyzed hereinafter. The

Fuzzy clustering technique shall be used to classify the possible fault locations, which can be near to any node in the network or near to a chosen set of nodes based on the operator experience, into groups. The optimum number of groups (clusters) is computed using validity clustering technique. The remaining 13 cases are used as test cases. Euclidean

parameters that selected to build the feature matrix are shown in Table 5.

given to show the feasibility for the using the fuzzy logic to locate the fault location.

poles. Area2 power system also contains two 3.16 MVA stand by generators.

and tripping information for all bus-bar compartments.

**techniques**

mean classification techniques.

*3.3.1. Network description*

*3.3.2. Feature Matrix*

## *3.2.4. Discussion*

The main objective of in this example to generate a human model which reflects the interaction between human factors in addition to other factors such as managerial factors, accident information, and work information, using fuzzy clustering techniques. Secondly, to predict annual rate of accident for any sample of workers by applying fuzzy logic techniques. Some of the important results obtained can be summarized as follows:

## **3.3. Fault location in distribution networks using fuzzy c-mean clustering techniques**

In last section of this chapter we shall studies an existing 13.8 kilovolt distribution network which, serves an oil production field spread over an area of approximately sixty kilometers square, in order to locate any fault that may occur anywhere in the network using fuzzy cmean classification techniques.

In addition, we shall introduce several methods for normalizing data and selecting the optimum number of clusters in order to classify data. Results and conclusion shall be also given to show the feasibility for the using the fuzzy logic to locate the fault location.

## *3.3.1. Network description*

426 Fuzzy Controllers – Recent Advances in Theory and Applications

significant poorness of the general health.

accidents) can be done via the fuzzy engine.

accident of this particular labourer.

suitable workers for any particular task

particular company.

hearing ability.

*3.2.4. Discussion*

c. Figure 38(Appendix) correlates the laborers general health and specific weight with their annual accident rate considering all other factors are `average`. It is clear that the annual accident rate increases with the significant increase of the specific weight and

d. In Figure 39 (Appendix), the laborers educational level and safety adherence are plotted against their annual accident rate considering all other factors are `average`. It is noticed from this illustration that the laborers with the two educational level extremes are

e. Figure 40 (Appendix) correlates the laborers optical status and hearing ability with their annual accident rate considering all other factors to be `average`. It is clear that the annual accident rate increases with the significant poorness of the optical status and

f. From the last plot (Figure 41 - Appendix), it is clear that if the company is not keen on the level of safety, the annual accident rate will be increased, considering all other factors are `average`. In addition, it is noticeable from the figure that the tasks that are

The main objective of in this example to generate a human model which reflects the interaction between human factors in addition to other factors such as managerial factors, accident information, and work information, using fuzzy clustering techniques. Secondly, to predict annual rate of accident for any sample of workers by applying fuzzy logic

a. Fuzzy clustering techniques can be used to build a model that characterizes the different features of workers in local construction field against their rate of accidents. b. The model obtained from clustering is considered as rules to be used in a fuzzy logic engine to predict the rate of accidents for any worker in the construction field. c. For any specific case, 231 correlations (between any two features and the rate of

d. Optimum training to improve the safety attitude for certain laborer with minimum cost can be estimated by analyzing the correlation between the level of training and rate of

e. By analyzing the correlation between level of safety importance and rate of accidents for the workers in a particular company, the limits of the effective safety improvement

f. A similar technique can be applied to a particular company to predict the rate of accident in order to estimate the insurance rate for the people who work in this

g. Using the accident rate fuzzy prediction techniques companies can select the most

can be predicted in order to evaluate the investment in this direction.

exposed to accidents more than the `average` educational level laborers.

considered to be dangerous and need safety-gear are a source of accidents.

techniques. Some of the important results obtained can be summarized as follows:

A joint venture oil company possesses two production areas, Area1 and Area2. For each area a power distribution system is provided. The power supply required for the two fields is provided from 130 MVA capacity power generation plant located in Area1. Twenty MW is transmitted to Area2 via 35KM long over head transmission line (OHTL) on wooden poles. Area2 power system also contains two 3.16 MVA stand by generators.

Area2 distribution system consists of three radial overhead transmission lines, 8-10 KM long each, serve submersible oil pumps and other loads scattered in the field (60-kilometer square). These overhead transmission lines have neither differential relays nor sectionalizing fuses. Programmable logic Controller (PLC) is connected to the incomers, outgoings and generators auxiliaries at Area2 power station. This PLC records the running and tripping information for all bus-bar compartments.

Due to the aging of the system, remote area problems and harsh desert weather, repeated faults are experienced in the grid. Because of unavailability of differential relays and/or sectionalizing fuses, it is very difficult and long time is consumed to locate any fault in this network. The problem is reflecting passively on the oil productivity of this important area.

## *3.3.2. Feature Matrix*

In order to measure the features of the faults at 176 nodes of the network; load flow study is implemented to determine the respective power loss for each short circuit case and also short circuit study is carried out to determine the feature vector for each short circuit case.

The results, obtained from the load flow study and the short circuit study, shall be used to form the network *feature matrix,* which will be clustered and analyzed hereinafter. The parameters that selected to build the feature matrix are shown in Table 5.

Fuzzy clustering technique shall be used to classify the possible fault locations, which can be near to any node in the network or near to a chosen set of nodes based on the operator experience, into groups. The optimum number of groups (clusters) is computed using validity clustering technique. The remaining 13 cases are used as test cases. Euclidean

distance technique is implemented to find out the group of nodes, which the fault may be found near to, for each test case.

New Areas in Fuzzy Application 429

d. The norms between each test data and the full data in the chosen cluster is examined, accordingly the corresponding cluster for each test data is decided based on the

e. Alpha-cut defuzzification is used with alpha equal to 90% of the average of norms

f. The nearest node to the fault is checked in order to determine whether it is included in

Matlab program is written to implement the above six steps and the results are analyzed

a. The optimum number of cluster (twenty five) is determined by implementing cluster validity technique discussed in 2.1. The result is illustrated in Fig 35 which give the

b. Table 6 shows the effort saved to locate 13 fault cases. From this table, it can be noticed

Effort Saving = 1-Ratio between the number of possible locations to the total

 Percentage of successful trials (Ratio of the cases of including the nearest node in the possible locations to the total number of testing cases\*100%) = 92%. (46) Average effort saving (Summation of the Effort Saving percentages divided by the

It is important to notice that fuzzy cluster technique failed to locate the fault of test data number 4 due to the lack of information near to this location. However, it is expected that

for more available information the performance of this technique will be improved.

relation between the number of clusters and the corresponding error where:

**Figure 35.** The relation between number of clusters and the corresponding error

number of nodes\*100%. (45)

total number of the cases) =75% (47)

Also, we can define the following terms:

minimum norm obtained from the said examination.

between the cluster center and its data.

the possible locations or not.

*<IP_ADDRESS>. Results and discussion*

and summarized as follows:

that:

**Table 5.** Summary of the parameters that are selected to build the feature matrix.

**Assumptions:** The following assumptions are considered:

Assumption (a) is valid since the transmission lines are short [47]. For assumption (b), the same work can be repeated for all other types of failures.

## *3.3.3. Fault location using column maximum normalization*