MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2015/268470 268470 Research Article A Fundamental Wave Amplitude Prediction Algorithm Based on Fuzzy Neural Network for Harmonic Elimination of Electric Arc Furnace Current Lei Wanjun http://orcid.org/0000-0001-8830-046X Wang Yanxia Wang Lu Cao Hui Guo William State Key Laboratory of Electrical Insulation and Power Equipment School of Electrical Engineering Xi’an Jiaotong University Xi’an, Shaanxi 710049 China xjtu.edu.cn 2015 4102015 2015 11 04 2015 25 05 2015 4102015 2015 Copyright © 2015 Wanjun Lei et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

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 . 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 , and a new EAF-specific model based on field measurements of instantaneous furnace voltages and currents is presented . Moreover, a mathematical model of the 3-phase AC of EAF is proposed in a more macroscopic level of modeling . 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 . In particular, an ANN-based inductances on-line estimation of the swinging power cables is discussed for efficient power control of EAF . Nevertheless, opaqueness is major shortcoming of ANN. Hence, to deal with this problem, the fuzzy neural network (FNN) is proposed . 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 .

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 Section 2. 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 Figure 1. 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.

Diagram of EAF power supply system model.

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 Figure 3. 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 n th 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.

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 Figure 2. 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 i th 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 . 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 of P 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 minimize B 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 of f = 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 Table 1. 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.

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 - 04

For the number of inputs is 20, the RMSEP values of 10 times experiments are shown in Table 2. 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.

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 - 05

The structure of ANFIS used in the proposed algorithm.

Experiments results for 10 inputs.

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.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (51207126), the Delta Science & Technology Educational Development Program (DREG2011009), and the Fundamental Research Funds for the Central University.

Rahmat A. H. Mahdi T. E. Adaptive filter design based on the LMS algorithm for delay elimination in TCR-FC compensators ISA Transactions 2011 50 2 142 149 10.1016/j.isatra.2010.12.002 2-s2.0-79952452193 Liu Y.-J. Chang G. W. Hong R.-C. Curve-fitting-based method for modeling voltage-current characteristic of an ac electric arc furnace Electric Power Systems Research 2010 80 5 572 581 10.1016/j.epsr.2009.10.015 2-s2.0-77149153063 Göl M. Salor Ö. Alboyací B. Mutluer B. Çadírcí I. Ermis M. A new field-data-based EAF model for power quality studies IEEE Transactions on Industry Applications 2010 46 3 1230 1242 10.1109/TIA.2010.2046280 2-s2.0-77952711527 Logar V. Dovžan D. Škrjanc I. Mathematical modeling and experimental validation of an electric arc furnace ISIJ International 2011 51 3 382 391 10.2355/isijinternational.51.382 2-s2.0-79960007640 Hooshmand R. A. Noohi S. Neural network-based algorithm for harmonic analysis of industrial systems using asynchronous simulations in time and frequency domains IEEJ Transactions on Electrical and Electronic Engineering 2010 5 6 688 694 10.1002/tee.20593 2-s2.0-79251520682 Chang G. W. Chen C.-I. Teng Y.-F. Radial-basis-function-based neural network for harmonic detection IEEE Transactions on Industrial Electronics 2010 57 6 2171 2179 10.1109/tie.2009.2034681 2-s2.0-77952283808 Han Y. Xu L. Khan M. M. Yao G. Zhou L.-D. Chen C. Study on a novel approach to active power filter control using neural network-based harmonic identification scheme Electrical Engineering 2010 91 6 313 325 10.1007/s00202-009-0143-9 2-s2.0-77951025792 Temurtas H. Temurtas F. An application of neural networks for harmonic coefficients and relative phase shifts detection Expert Systems with Applications 2011 38 4 3446 3450 10.1016/j.eswa.2010.08.131 2-s2.0-78650693545 Ketabi A. Sadeghkhani I. Feuillet R. Using artificial neural network to analyze harmonic overvoltages during power system restoration European Transactions on Electrical Power 2011 21 7 1941 1953 10.1002/etep.523 2-s2.0-80055045795 Filho F. Maia H. Z. Mateus T. H. A. Ozpineci B. Tolbert L. M. Pinto J. O. P. Adaptive selective harmonic minimization based on ANNs for cascade multilevel inverters with varying DC sources IEEE Transactions on Industrial Electronics 2013 60 5 1955 1962 10.1109/TIE.2012.2224072 2-s2.0-84873334154 Guellal A. Larbes C. Bendib D. Hassaine L. Malek A. FPGA based on-line artificial neural network selective harmonic elimination pwm technique International Journal of Electrical Power & Energy Systems 2015 68 33 43 10.1016/j.ijepes.2014.11.030 Khoshkhoo H. Sadeghi S. H. H. Moini R. Talebi H. A. An efficient power control scheme for electric arc furnaces using online estimation of flexible cable inductance Computers and Mathematics with Applications 2011 62 12 4391 4401 10.1016/j.camwa.2011.10.009 2-s2.0-82255175592 Jang J.-S. R. ANFIS: adaptive-network-based fuzzy inference system IEEE Transactions on Systems, Man and Cybernetics 1993 23 3 665 685 10.1109/21.256541 2-s2.0-0027601884 Dehkordi B. M. Moallem M. Parsapour A. Predicting foaming slag quality in electric arc furnace using power quality indices and fuzzy method IEEE Transactions on Instrumentation and Measurement 2011 60 12 3845 3852 10.1109/tim.2010.2085190 2-s2.0-81255134141 Lee K. C. Gardner P. Adaptive neuro-fuzzy inference system (ANFIS) digital predistorter for RF power amplifier linearization IEEE Transactions on Vehicular Technology 2006 55 1 43 51 10.1109/tvt.2005.861171 2-s2.0-32144439918