# A Fundamental Wave Amplitude Prediction Algorithm Based on Fuzzy Neural Network for Harmonic Elimination of Electric Arc Furnace Current

Authors: Wanjun Lei; Yanxia Wang; Lu Wang; Hui Cao
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

Electric arc furnace (EAF) causes the harmonics to impact on the supply network greatly and harmonic elimination is a very important research work for the power quality associated with EAF. In the paper, a fundamental wave amplitude prediction algorithm based on fuzzy neural network for harmonic elimination of EAF current is proposed. The proposed algorithm uses the learning ability of the neural network to refine Takagi-Sugeno type fuzzy rules and the inputs are the average of the current measured value in different time intervals. To verify the effectiveness of the proposed algorithm, some experiments are performed to compare the proposed algorithm with the back-propagation neural networks, and the field data collected at an EAF are used in the experiments. Moreover, the measured amplitudes of fundamental waves of field data are obtained by the sliding-window-based discrete Fourier transform on the field data. The experiments results show that the proposed algorithm has higher precision. The real curves also verify that the amplitude of fundamental wave current could be predicted accurately and the harmonic elimination of EAF would be realized based on the proposed algorithm.


---

## Body

## 1. Introduction

Electric arc furnace (EAF) smelts the raw metal by generating arcs and has been widely used in the steel industry. For the characteristics of electric arc and meltdown processes, EAF causes some undesirable disturbing effects, such as harmonics and flicker, which would have high impact on the supply network [1]. Therefore, harmonic elimination is of important theoretical significance and practical motivation for the power quality associated with EAF.Traditionally, most researches focus on the steady-state modeling for harmonic elimination and some improved models of EAF are presented recently. The cubic spline interpolation is proposed to model the voltage-current characteristic of EAF [2], and a new EAF-specific model based on field measurements of instantaneous furnace voltages and currents is presented [3]. Moreover, a mathematical model of the 3-phase AC of EAF is proposed in a more macroscopic level of modeling [4]. However, due to the instability of smelting process and the randomness of the arc generation, establishing a precise model of the EAF is extremely difficult. Moreover, because some physical responses of EAF cannot be obtained from the theoretical study easily, the assumptions of parameters would simplify the EAF model and the analysis results would be affected. Artificial neural network (ANN) learns the relationship between output and input by adjusting the interconnections of neurons and has been adopted for harmonic elimination [5–11]. In particular, an ANN-based inductances on-line estimation of the swinging power cables is discussed for efficient power control of EAF [12]. Nevertheless, opaqueness is major shortcoming of ANN. Hence, to deal with this problem, the fuzzy neural network (FNN) is proposed [13]. FNN is a combined system with ANN and fuzzy logic technique. The fuzzy concepts could improve the transparency for understanding the inner working of ANN, and for fuzzy logic technique, the fuzzy rules have clear semantic meanings; namely, FNN has the advantages of learning optimization ability of ANN and the human-like thinking of fuzzy logic technique. In addition, the parameter of FNN could be determined easily. Adaptive neurofuzzy inference system (ANFIS) is one kind of FNN and a first-order Takagi-Sugeno fuzzy system. Since ANFIS could excellently estimate a system with uncertainty, a new method based on ANFIS is proposed to evaluate the slag quality in EAF using power quality indices [14].This paper proposes a fundamental wave amplitude prediction algorithm based on FNN for harmonic elimination of EAF current. For the proposed algorithm, ANFIS is used to predict the real-time amplitude of fundamental wave current of EAF, which could provide a guarantee for the realization of harmonic elimination. The inputs are the average of the current measured value in different time intervals. To verify the effectiveness of the proposed algorithm, some experiments are performed on the field data collected at an EAF. In the experiments, the proposed algorithm is compared with the back-propagation neural networks (BPNN) for estimating the fundamental wave amplitude. Moreover, the measured values of fundamental wave amplitude of field data are obtained by the sliding-window-based discrete Fourier transform (SW-DFT) on the field data. The paper is organized as follows. The power supply system model and the equivalent circuit of EAF are described in Section2. In Section 3, the proposed algorithm is presented in detail. In Section 4, the experiments results are discussed to verify the effectiveness of the proposed algorithm. Finally, Section 5 concludes the paper.


## 2. Electric Arc Furnace

The EAF power supply system model is shown in Figure1. T
1 is a supply transformer and T
2 is a furnace transformer. Z
S is the grid impedance and Z
C is the busbar impedance of EAF. According to the state variables of EAF, such as the arc length, the current value, and the temperature, the controller could regulate the stalls of output voltage of the furnace transformer and the electrode position. The equivalent circuit is shown in Figure 2. R
s and L
s are equivalent to the grid impedance. The transformer with a ratio of T
p:1 is equivalent to T
1 and T
2. R
c and L
c are equivalent to the busbar impedance. V
s is the power source. Time-varying resistance R
E and V
b are the arc equivalent model.Figure 1
Diagram of EAF power supply system model.Figure 2
Equivalent circuit of EAF.


## 3. The Proposed Algorithm

According to the characteristics of EAF, a fundamental wave amplitude prediction algorithm is presented in the paper and the architecture of the proposed algorithm is shown in Figure3. By predicting the fundamental wave amplitude, the nonfundamental components of EAF current are separated out and used as the command signal to control the power electronic converter to generate the harmonic current, which could counteract the harmonic of EAF, so the grid current only contains the fundamental wave. The output of the proposed algorithm, i
f
(
t
+
n
+
1
), is the prediction amplitude of fundamental wave current and n is the number of inputs. The inputs of the proposed algorithm are the average of the current measured value in different time intervals; for example, i
m
(
t
n
) is the nth inputs and could be obtained by the following equation:(1)
i
m
t
n
=
1
k
∑
j
=
0
k
i
m
t
+
n
-
1
k
+
1
+
j
,where k is a fixed integer and equals the number of samplings in time intervals.Figure 3
The proposed algorithm for harmonic elimination of EAF current.For the proposed algorithm, ANFIS is adopted for predicting the amplitude of fundamental wave current of EAF, and the architecture is shown in Figure2. The algorithm has five layers. In the same layer, the nodes are with the same activation function. To facilitate the implementation of the proposed algorithm, each input uses the same linguistic term, and the linguistic terms set is {
A
1
,
A
2
,
…
,
A
m
}, where m is the number of fuzzy linguistic terms. Since the membership functions are determined automatically by training process, the value of m may not be set beforehand. In layer 1, each node function is the membership values of each input with respect to its linguistic term A
L, L
∈
1,2
,
…
,
m, and the Gaussian function is used as the membership functions. In layer 2, each node is a circle node labeled ∏ and plays the role of a simple multiplier. There are m
n nodes in this layer and the output of each node represents the firing strength of the rule. In layer 3, each node is a circle node labeled N and calculates the ratio of the activation level of the ith rule to the total of all activation level. In layer 4, there are m
n nodes in the layer and each node calculates the contribution of the overall output; namely, it is simply the product of the normalized firing strength and the function of consequent, which is a first-order polynomial. Layer 5 has only one node labeled ∑ to indicate that it performs the function of a simple summer. The output of the node is final output, i
f
(
t
+
n
+
1
), which is the prediction amplitude of fundamental wave current and is given by(2)
i
f
t
+
n
+
1
=
∑
q
=
1
m
n
w
i
∑
p
=
1
m
n
w
p
·
∑
c
=
1
m
a
c
q
x
c
+
b
q
,where a
c
q and b
q are the consequent parameters.For the proposed algorithm, a hybrid learning process is adopted to tune all the modifiable parameters to make the output match the measured value. The hybrid learning process adjusts the consequent parameters in a forward pass and the antecedent parameters in a backward pass [15]. This is a linear combination of the modifiable parameters. For this observation, we can divide the parameter set into two sets. Consider(3)
P
=
P
1
⊕
P
2
,where P is the set of total parameters. P
1 is the set of antecedent parameters. P
2 is the set of consequent parameters. ⊕ is the sum operation.For the forward path, the least square method is used to identify the consequent parameters. Now, for a given set of values of elements ofP
1, we can obtain the following matrix equation based on the train data:(4)
B
X
=
Y
,where B contains the unknown parameters in P
2, X is the input matrix, and Y is the output.A least squares estimate is sought to minimizeB
X
-
Y
2, which is a problem for linear regression. Thus, the calculation for the parameters set B is:(5)
B
∗
=
X
T
X
-
1
X
T
Y
.In the backward path, the error signals propagate backward. The antecedent parameters are updated by the descent method, through minimizing the error between the output of the network and that of real system, which is a quadratic function:(6)
E
Ψ
=
1
2
∑
d
=
1
N
y
d
-
i
f
d
2
,where Ψ represents the training data set and N is the number of the training data.In a recursive manner, with respect toΨ, the update of the parameters of membership function can be written as(7)
c
A
L
z
+
1
=
c
A
L
z
-
η
∂
E
∂
c
A
L
σ
A
L
z
+
1
=
σ
A
L
z
-
η
∂
E
∂
σ
A
L
,where c
A
L and σ
A
L are the centers of the function and the widths of the membership function. η is the learning rate and η
>
0.Therefore, the amplitude of fundamental wave current could be predicted correctly by the proposed algorithm and the harmonic elimination of EAF would be realized. In the next section, the experiments results will show the effectiveness of the proposed algorithm.


## 4. Experiments Results

In the section, some experiments are performed to evaluate the effectiveness of the proposed algorithm. The experiments focus on comparing the proposed algorithm with BPNN for estimating the fundamental wave amplitude and the field data are used in the experiments. The field data are collected at an EAF at the second side of furnace transformer and the sampling rate is 50 kHz (1000 amples/per cycle for the ideal case off = 50 Hz). To facilitate analysis, the values of field data are normalized. For the experiments, the 4000 samples from field data are divided into the training set and the test set based on the Monte Carlo cross-validation at the ratio that equaled 4 : 1. The root mean squares error of prediction (RMSEP) is used to estimate the predictive ability of the proposed algorithm and BPNN. For the proposed algorithm, the training epoch number is 100, the training error goal is zero, the initial step size is 0.01, the step size decrease rate is 0.7, and the step size increase rate is 1.3. For BPNN, the number of training epochs is 100, the learning rate is 0.01, the number of hidden neurons is 10, and the training error goal is zero. Moreover, the measured amplitudes of fundamental waves of field data are obtained by SW-DFT on the field data. Because the inputs are the average of measured value in different time intervals, the number of inputs may not be fixed and could be set according to the analysis demand. Hence, the experiments are performed with different numbers of inputs and the results are discussed in the following.For the number of inputs is 10, the RMSEP values of 10 times experiments are shown in Table1. Since assessing the sampling bias, the division of field data and the calculation process would be implemented repeatedly. The average RMSEP values of the proposed algorithm and BPNN are 4.8764E − 05 and 1.8450E − 04, respectively. The proposed algorithm is more effective and improves the prediction precision by 73.57% than BPNN. According to the average RMSEP value, the measured value curve and predictions curves of the proposed algorithm and BPNN for other 2000 samples from field data are shown in Figure 5, and the models of the proposed algorithm and BPNN are based on the number 2 experiment results. It also shows that the proposed algorithm has higher effectiveness.Table 1
RMSEP values of 10 times experiments for 10 inputs.
Number
The proposed algorithm
BPNN
1
5.73
E
-
05
2.6432
E
-
04
2
4.54
E
-
05
1.7707
E
-
04
3
5.81
E
-
05
1.2414
E
-
04
4
4.94
E
-
05
1.5082
E
-
04
5
4.08
E
-
05
1.6229
E
-
04
6
3.26
E
-
05
1.4230
E
-
04
7
4.34
E
-
05
1.7170
E
-
04
8
5.53
E
-
05
1.8885
E
-
04
9
3.47
E
-
05
2.4323
E
-
04
10
7.08
E
-
05
2.2032
E
-
04
The average RMSEP
4.8764
E
-
05
1.8450
E
-
04For the number of inputs is 20, the RMSEP values of 10 times experiments are shown in Table2. The average RMSEP values of the proposed algorithm and BPNN are 3.92E − 06 and 6.77E − 05, respectively. The proposed algorithm is more effective and improves the prediction precision by 94.21% than BPNN. According to the average RMSEP value, the measured value curve and predictions curves of the proposed algorithm and BPNN for the 2000 samples used in Figure 4 are shown in Figure 6, and the models of the proposed algorithm and BPNN are based on the number 9 experiment results. The predictive ability of BPNN is decreased with the number of inputs being increased and that of the proposed algorithm is still better. It also shows that the proposed algorithm has higher effectiveness.Table 2
RMSEP values of 10 times experiments for 20 inputs.
Number
The proposed algorithm
BPNN
1
3.92
E
-
06
6.77
E
-
05
2
3.73
E
-
06
7.23
E
-
05
3
3.72
E
-
06
3.73
E
-
05
4
3.60
E
-
06
4.12
E
-
05
5
3.63
E
-
06
9.32
E
-
05
6
3.45
E
-
06
2.59
E
-
05
7
3.45
E
-
06
2.17
E
-
05
8
3.50
E
-
06
9.61
E
-
05
9
3.68
E
-
06
4.65
E
-
05
10
3.47
E
-
06
9.52
E
-
05
The average RMSEP
3.92
E
-
06
6.77
E
-
05Figure 4
The structure of ANFIS used in the proposed algorithm.Figure 5
Experiments results for 10 inputs.Figure 6
Experiments results for 20 inputs.Therefore, the experiments results verify that the effectiveness of the proposed algorithm is higher; namely, the amplitude of fundamental wave current could be predicted accurately, and the harmonic elimination of EAF would be realized based on the proposed algorithm.


## 5. Conclusions

The paper proposes a fundamental wave amplitude prediction algorithm based on FNN for harmonic elimination of EAF current. The proposed model has some advantages as follows. First, it has higher prediction capability for fundamental wave amplitude. Second, it is almost not affected by the number of inputs changed. Third, it uses the learning ability of the neural network to implement and refine Takagi-Sugeno type fuzzy rules to predict the fundamental wave amplitude. Fourth, by the hybrid learning process, all the best parameters of membership functions are obtained. The experiments results also verify the effectiveness of the proposed algorithm. Since the train time may be affected by the number of inputs, in the future research work, we will use some feature extraction approaches to further improve the training performance of the proposed algorithm.



---

*Source: 268470-2015-10-04.xml*