# An ELM Based Online Soft Sensing Approach for Alumina Concentration Detection

Authors: Sen Zhang; Xi Chen; Yixin Yin
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

The concentration of alumina in the electrolyte is of great significance during the production of aluminum; it may affect the stability of aluminum reduction cell and the current efficiency. However, the concentration of alumina is hard to be detected online because of the special circumstance in the aluminum reduction cell. At present, there is lack of fast and accurate soft sensing methods for alumina concentration and existing methods can not meet the needs for online measurement. In this paper, a novel soft sensing method based on a modified extreme learning machine (MELM) for online measurement of the alumina concentration is proposed. The modified ELM algorithm is based on the enhanced random search which is called incremental extreme learning machine in some references. It randomly chooses the input weights and analytically determines the output weights without manual intervention. The simulation results show that the approach can give more accurate estimations of alumina concentration with faster learning speed compared with other methods such as BP and SVM.


---

## Body

## 1. Introduction

In the industrial aluminum reduction cells, the stability of the alumina concentration is the key issue to maintain high efficiency during the production of aluminum. It is easy to lead to the occurrence of the so-called “anode effect” if the alumina content in the electrolyte becomes too low (e.g., 1%–1.5%). When it occurs, cell voltage rises abruptly to 30v–50v, which directly affects the energy balance of aluminum reduction cells. It is important to avoid the high alumina concentration because alumina-rich sludge will accumulate at the base of the cell and the cell operation can be seriously disrupted [1]. Therefore, how to detect the distribution of alumina concentration in the cell in real time is the key problem to control the production of aluminum. The site environment of aluminum reduction cell, which has the characteristics of large current, strong magnetic field, high temperature, high humidity, and a lot of dust, is very severe and complex. Aluminum reduction cell is a severe nonlinear, multi-input multioutput, slow time-variety, long time-delay, and coupled process. Therefore, it is a challenge to seek an online soft sensing method for alumina concentration. Numerous investigations have been carried out on the soft measurement method for alumina concentration by the researchers. A prediction model based on wavelet neural network was proposed by Li et al. [2]. The prediction method based on linear regression and orthogonal transform is applied to improve the accuracy of the alumina concentration forecast by Lin et al. [3]. Yan and Liang proposed a predictive model of aluminum reduction cell based on LS-SVM [4]. Li et al. [5] proposed a new fuzzy expert control method based on smart identification, multicontrol mode, and decision-making mechanism to achieve the alumina concentration prediction and real time control. The GM(1,1) model is introduced into the aluminum concentration estimate by Zhang et al. [6]. However, the computational burden of the above nonlinear predictive models is still large when the dimension of input variable increases; the learning speed and accuracy of networks are in general far slower and can not meet the requirement of real time detection.Figure1 shows the experimental environment in the ZunYi aluminium electrolysis factory in China. In industries, there often exist some crucial variables that can not be directly measured online due to the fact that no sensor is available in the certain complex environment or due to its high cost. Soft sensing technology is a way to solve the problems. The soft sensing technology is widely used and becomes one of the important developing directions of surveying and processing control area. In practice, the application and research of this technology have gotten more extensions and many related technologies based on it have emerged [7–10].Figure 1
Experimental environment.The main thought of the soft sensing technology is to build one model which uses variables that could be directly online measured as input and thus the variables to be estimated as output, to use many kinds of complex calculating and evaluation and to get the values of detecting variables by computer software [11–14]. So the main problem of soft sensing technology is how to build the relation model between detecting variables and other easy getting variables. By now, there are many methods to build the models, and many methods are in the trend of intersecting and mixing together. In all the methods, the artificial intelligence method is used often, such as the method based on model identification, the method based on artificial neural network, or the method based on fuzzy set theory. Some professors and scholars have done some research about the technology, such as global asymptotic stability of neural networks with multiple time-varying delays and fuzzy model-based robust networked control for a class of nonlinear systems proposed by Zhang et al. [15, 16], weighted least squares support vector machines, robustness and sparse approximation, proposed by Suykens et al. [17], Artificial Neural Network in process engineering proposed by Willis et al. [18], and the modified extreme learning machine method by Cao et al. [19]. Soft sensing technology is applied to various fields; Wang and Ren proposed soft sensing method for wastewater treatment based on BP neural network [20]. Rolling bearing fault detection based on the Teager Energy Operator and Elman Neural Network was proposed by Liu et al. [21].So far, the widely used method of determining alumina concentration in the industrial factory is to use the spectrometer to analyze the sampled electrolyte. It is an offline method. The improved ELM method is applied to build up the soft sensing online detection approach of alumina concentration in this paper. This ELM algorithm tends to provide the best generalization performance at extremely fast learning speed. The experimental result shows that the new method can produce the best generalization performance and learn much faster than the traditional prediction models.The following sections are organized as follows. Section2 shows the theory of the extreme learning machine. Section 3 proposes the modified ELM algorithm. Data selection and data preprocessing are shown in Section 4. Section 5 gives the simulation results. Section 6 summarizes the conclusions.


## 2. The Theory of Extreme Learning Machine

Huang proposed extreme learning machine (ELM) algorithm; ELM was originally proposed for the single-hidden-layer feedforward neural networks (SLFNs) and then it extends to the generalized SLFNs. ELM can randomly generate the input weights and the bias of hidden nodes. It uses the theory of least squares to get the output weights. The learning speed of ELM can be thousands of times faster than traditional feedforward network learning algorithms like backpropagation (BP) algorithm while obtaining better generalization performance and smaller training error [22, 23]. Figure 2 shows the single-hidden-layer feedforward network (SLFN) architecture.Figure 2
SLFN network architecture.Considering there areN arbitrary distinct samples (Xi,Ti), when Xi=[xi1,xi2,…,xin]T∈Rn and Ti=[ti1,ti2,…,tim]T∈Rm, standard SLFNs with L hidden neurons and activation function h(x) are mathematically modeled as(1)∑i=1L‍βihWi·Xj+bi=Oj,j=1,2,…,N,where Wi=[wi1,wi2,…,win]T is the weight vector connecting the ith hidden neuron and the input neurons, βi=[βi1,βi2,…,βim]T is the weight vector connecting the ith hidden neuron and the output neurons, and bi is the threshold of the ith hidden neuron.The fact that standard SLFNs withL hidden neurons with activation function h(x) can approximate these N samples with zero error means that ∑j=1LOj-Tj=0; there exist βi, Wi, and bi such that(2)∑i=1L‍βihWi·Xj+bi=Tj,j=1,2,…,N.The aboveN equations can be written compactly as (3)Hβ=T,where(4)HW1,…,WL,b1,…,bL,X1,…,XN=hW1·X1+b1⋯hWL·X1+bL⋮⋱⋮hW1·XN+b1⋯hWL·XN+bLN×L,β=β1T⋮βLTL×m,T=t1T⋮tLTN×m.H is called the hidden layer output matrix of the neural network. ELM is to minimize the training error as well as the norm of the output weight. Minimize Hβ-T2 and β.The minimal norm least squares method instead of the standard optimization method was used in the original implementation of ELM:(5)β^=HΓT,where HΓ is the Moore-Penrose generalized inverse of H. Several methods can be used to calculate the HΓ; these methods may include orthogonal projection, orthogonalization method, iterative method, and singular value decomposition (SVD), and ELM algorithm makes use of SVD, where HΓ=(HTH)-1HT, so (6)β^=HTH-1HTT.The algorithm ELM can be summarized as Three-Step Learning Model. Given a training set ℵ={Xi,Ti∣Xi∈Rn,Ti∈Rm,i=1,2,…,N}, an activation function h(x), and the hidden neuron number L, we have the following steps.Step 1.
Assign arbitrary input weighsWi and bias of hidden layer nodes bi,i=1,2,…,L.Step 2.
Calculate the hidden layer output matrixH.Step 3.
Calculate the output weightsβ^.


## 3. The Modified ELM Algorithm

### 3.1. Ridge Regression Based ELM Algorithm

Ridge Regression, which was proposed by Horel and Kennard in1970, proposes the idea to adding a small positive number on the main diagonal of the design matrix. Considering ELM algorithm, the estimation of β^* can be obtained by employing the following formula:(7)β^*=HTH+kI-1HTT,where k is the ridge parameter. Superficially, it is possible to get the inversion of design matrix HTH. The following is the new results of (8) and (9):(8)β^*=∑i=1N~‍HiTTi∑i=1N~‍Hi2+k=∑i=1N~‍HiTHiβi+εi∑i=1N~‍Hi2+k=∑i=1N~‍Hi2βi∑i=1N~‍Hi2+k+∑i=1N~‍HiTεi∑i=1N~‍Hi2+k,1Eβ^*=∑i=1N~‍Hi2∑i=1N~‍Hi2+kβ≠β,which denotes that Ridge Regression is biased.Consider(9)2Vβ^*=σ2∑i=1N~Hi2‍∑i=1N~Hi2+k2=σ2∑i=1N~λi‍∑i=1N~λi+k2<Vβ^,which denotes that Ridge Regression makes the estimation more stable.Consider(10)3MSEβ^*=1N~×σ2∑i=1N~Hi2+β2k2‍∑i=1N~Hi2+k2=σ2∑i=1N~λi‍∑i=1N~λi+k2+β2k2∑i=1N~λi+k2=1N~γ1k+γ2k.Hoerl and Kennard had proved that Ridge Regression has less mean square error than the ordinary regression under the proper ridge parameter. It is as follows:(11)dγ1(k)dk=-2σ2∑i=1N~‍λi∑i=1N~λi+k3,dγ2(k)dk=2β2k∑i=1N~‍λi∑i=1N~λi+k3.So, whenk→0+, (12)limk→0+dγ1(k)dk=-2σ2∑i=1N~λi2<0,limk→0+dγ2(k)dk=0.From the above equations,MSE(β^*) is an increasing function of k when k∈(0,δ), where dγ1(k)/dkk=δ=dγ2(k)/dkk=δ. Therefore, the selection of parameter k is essential to the performance of Ridge Regression. In our ER-ELM algorithm, a method to determine the ridge parameter proposed by Huang [24] is used:(13)k=σ^2β^,where β^ is the ordinary ELM algorithm estimation and σ^2=∑i=1N(yi-β^xi)2/v,v=n-1.


### 3.2. To Select the Number of Hidden Nodes Based on the Improved ELM

The Error Minimized Extreme Learning Machine algorithm starts from a small size of ELM hidden layer and adds random hidden node (nodes) to the hidden layer, while the output weights are updated incrementally.Suppose a SLFN,H1=H(w1,…,wl0,b1,…,bl0,k1), denotes the hidden layer output matrix with l0 hidden nodes and ridge parameter k1 calculated by (19). Considering the poor performance due to lower number of hidden nodes, additional l hidden nodes are added to the SLFN. A new hidden layer output matrix H2 is composed of H1 and other l extra hidden nodes as(14)H2=H1,H,where H=G(wl0+1,bl0+1,x1,k2)⋯G(wl1,bl1,x1,k2)⋮⋱⋮G(wl0+1,bl0+1,xN,k2)⋯G(wl1,bl1,xN,k2)N×(l1-l0),l1-l0=l, and k2 is ridge parameter of HTH.Huang had provedE(H2)=min⁡H2β2-T≤E(H1)=min⁡H1β1-T, where E(H) denotes the output error function of SLFNs. We set a stopping criterion for the iterative algorithm as follows:(15)EHn+1-EHn<ε,where ε>0 is called the target error.Consider(16)H2+=H2TH2+K-1H2T=H1THTH1H+K100K2-1H1THT,where K=kI.In order to facilitate the calculation, the inversion is substituted as follows:(17)Φ=φ11φ12φ21φ22=H1THTH1H+K100K2-1,where(18)φ11=H1TH1+K1-1+H1TH1+K1-1H1THR-1HTH1H1TH1+K1-1,φ12=-H1TH1+K1-1H1THR-1,φ21=-R-1HTH1H1TH1+K1-1,φ22=R-1,R=HTH+K2-HTH1H1TH1+K1-1H1TH=HTH+K2-HTH1H1+H.So(19)H2+=UD=φ11H1T+φ12HTφ21H1T+φ22HT,(20)D=R-1HT-R-1HTH1H1TH1+K1-1H1T=R-1HT-R-1HTH1H1+=HTH+K2-HTH1H1+H-1HTI-H1H1+=HTI-H1H1+H+K2-1HTI-H1H1+=HTMH+K2-1HTM,where M=I-H1H1+.Similarly, we getU=H1+-H1+HD.Now, a new hidden layer output matrix is obtained which has less output error. Then, we can update the output weight matrix based on the new hidden layer output matrix.


## 3.1. Ridge Regression Based ELM Algorithm

Ridge Regression, which was proposed by Horel and Kennard in1970, proposes the idea to adding a small positive number on the main diagonal of the design matrix. Considering ELM algorithm, the estimation of β^* can be obtained by employing the following formula:(7)β^*=HTH+kI-1HTT,where k is the ridge parameter. Superficially, it is possible to get the inversion of design matrix HTH. The following is the new results of (8) and (9):(8)β^*=∑i=1N~‍HiTTi∑i=1N~‍Hi2+k=∑i=1N~‍HiTHiβi+εi∑i=1N~‍Hi2+k=∑i=1N~‍Hi2βi∑i=1N~‍Hi2+k+∑i=1N~‍HiTεi∑i=1N~‍Hi2+k,1Eβ^*=∑i=1N~‍Hi2∑i=1N~‍Hi2+kβ≠β,which denotes that Ridge Regression is biased.Consider(9)2Vβ^*=σ2∑i=1N~Hi2‍∑i=1N~Hi2+k2=σ2∑i=1N~λi‍∑i=1N~λi+k2<Vβ^,which denotes that Ridge Regression makes the estimation more stable.Consider(10)3MSEβ^*=1N~×σ2∑i=1N~Hi2+β2k2‍∑i=1N~Hi2+k2=σ2∑i=1N~λi‍∑i=1N~λi+k2+β2k2∑i=1N~λi+k2=1N~γ1k+γ2k.Hoerl and Kennard had proved that Ridge Regression has less mean square error than the ordinary regression under the proper ridge parameter. It is as follows:(11)dγ1(k)dk=-2σ2∑i=1N~‍λi∑i=1N~λi+k3,dγ2(k)dk=2β2k∑i=1N~‍λi∑i=1N~λi+k3.So, whenk→0+, (12)limk→0+dγ1(k)dk=-2σ2∑i=1N~λi2<0,limk→0+dγ2(k)dk=0.From the above equations,MSE(β^*) is an increasing function of k when k∈(0,δ), where dγ1(k)/dkk=δ=dγ2(k)/dkk=δ. Therefore, the selection of parameter k is essential to the performance of Ridge Regression. In our ER-ELM algorithm, a method to determine the ridge parameter proposed by Huang [24] is used:(13)k=σ^2β^,where β^ is the ordinary ELM algorithm estimation and σ^2=∑i=1N(yi-β^xi)2/v,v=n-1.


## 3.2. To Select the Number of Hidden Nodes Based on the Improved ELM

The Error Minimized Extreme Learning Machine algorithm starts from a small size of ELM hidden layer and adds random hidden node (nodes) to the hidden layer, while the output weights are updated incrementally.Suppose a SLFN,H1=H(w1,…,wl0,b1,…,bl0,k1), denotes the hidden layer output matrix with l0 hidden nodes and ridge parameter k1 calculated by (19). Considering the poor performance due to lower number of hidden nodes, additional l hidden nodes are added to the SLFN. A new hidden layer output matrix H2 is composed of H1 and other l extra hidden nodes as(14)H2=H1,H,where H=G(wl0+1,bl0+1,x1,k2)⋯G(wl1,bl1,x1,k2)⋮⋱⋮G(wl0+1,bl0+1,xN,k2)⋯G(wl1,bl1,xN,k2)N×(l1-l0),l1-l0=l, and k2 is ridge parameter of HTH.Huang had provedE(H2)=min⁡H2β2-T≤E(H1)=min⁡H1β1-T, where E(H) denotes the output error function of SLFNs. We set a stopping criterion for the iterative algorithm as follows:(15)EHn+1-EHn<ε,where ε>0 is called the target error.Consider(16)H2+=H2TH2+K-1H2T=H1THTH1H+K100K2-1H1THT,where K=kI.In order to facilitate the calculation, the inversion is substituted as follows:(17)Φ=φ11φ12φ21φ22=H1THTH1H+K100K2-1,where(18)φ11=H1TH1+K1-1+H1TH1+K1-1H1THR-1HTH1H1TH1+K1-1,φ12=-H1TH1+K1-1H1THR-1,φ21=-R-1HTH1H1TH1+K1-1,φ22=R-1,R=HTH+K2-HTH1H1TH1+K1-1H1TH=HTH+K2-HTH1H1+H.So(19)H2+=UD=φ11H1T+φ12HTφ21H1T+φ22HT,(20)D=R-1HT-R-1HTH1H1TH1+K1-1H1T=R-1HT-R-1HTH1H1+=HTH+K2-HTH1H1+H-1HTI-H1H1+=HTI-H1H1+H+K2-1HTI-H1H1+=HTMH+K2-1HTM,where M=I-H1H1+.Similarly, we getU=H1+-H1+HD.Now, a new hidden layer output matrix is obtained which has less output error. Then, we can update the output weight matrix based on the new hidden layer output matrix.


## 4. Data Selection and Preprocessing

In the aluminum production, through referring to relative documents and soliciting experts opinion, there are many factors that affect the alumina concentration, such as alumina feeding speed, cell voltage, series current, current of anode rod, voltage between anode rod and cathode bar, and bath temperature [25]. Existing predictive method of alumina concentration is to get an average alumina content in the aluminum reduction cell; the distribution of alumina concentration in the cell is not known yet. At present, people want to know the distribution situation of alumina content in the cell in order to know the process of production better and control the production of aluminum more accurately. Our experimental aluminum reduction cell has 24 anode rods; we can get the current signal of 24 anode rods, voltage signal between anode rods and cathode bars, and alumina concentration signal below the 24 anode rods. Through experiment and analysis, we find that the voltage signal between anode rod and cathode bar can reflect the alumina concentration below the anode rod more accurately. So we can use the voltage signal between anode rod and cathode bar in the aluminum reduction cell to reflect the distribution of alumina concentration. To use the ELM algorithm, we decide to select voltage between anode rod and cathode bar variable as model input according to experiments, and output variable is alumina concentration in the electrolyte below the anode rod.It is necessary to do the data denoising preprocessing before using the algorithm, because the process of data collection may introduce noise, where data collection is affected by interference and influence of all kinds of noise signal. At present, there are two denoising methods, the traditional filtering method and the wavelet denoising method. Traditional denoising method is based on Fourier analysis and can always be used in the environment where signal and noise are very small. While wavelet analysis is known as the “microscope” of mathematical analysis, it is a time-frequency analysis method of signal, with the characteristics of multiresolution analysis [26]. Different signals may choose the different denoising methods. Through experiment and comparision, we decide to use the wavelet denoising method in view of the voltage signal between anode rod and cathode bar.


## 5. Simulative Results of Alumina Concentration

In this paper, we select voltage between anode rod and cathode bar as model input and alumina concentration in the electrolyte below the anode rod as model output.First, ELM algorithm is applied to build soft sensing method model of alumina concentration. In general, active functions play an important role in computing of neural networks. Widely used active functions in the ELM algorithm aresig, sin, hardlim, and radbas. Through comparison, we find that the sig function has more outstanding performance than sin, hardlim, and radbas. So we apply sig function as active function in the ELM algorithm. In order to make comparison, we also use the BP and SVM algorithm to build soft sensing models of alumina concentration [27, 28]. The BP parameters are chosen as 17 hidden layer neurons and 1 output layer neuron and transfer function of hidden layer is tan-sigmoid and transfer function of output layer is linear.Data in this paper came from a 350 kA prebaked aluminum reduction cell in the ZunYi aluminium electrolysis factory; Figure3 shows experimental aluminum reduction cell and anode rod. We got the data of voltage between anode rod and cathode bar through the voltage measurement instruments. At the same time, we obtained the electrolyte below the experimental anode rod in the experimental aluminum reduction cell, and then we used the X fluorescence spectrometer in the laboratory to analyze the alumina concentration. Figure 4 shows the sampled electrolyte. We got the experimental data through this method. We select 100 pairs as samples to construct sample set (Xi,Ti), where Xi∈R denotes voltage between anode rod and cathode bar and Ti∈R denotes alumina concentration obtained at the same time with the voltage siginal.Figure 3
The voltage between number 20 anode rod and cathode bar and the alumina concentration in the electrolyte below the anode rod.Figure 4
The voltage between number 21 anode rod and cathode bar and the alumina concentration in the electrolyte below the anode rod.The ELM, BP, and SVM soft sensing models are trained by the same sample set. All the data are preprocessed before training and simulating in the algorithm. For the purpose of comparing the three models quantitatively, we substitute tenXi of new sample (out of the sample set) into ELM, BP, and SVM models, respectively, to calculate corresponding T^i; simulation results of ELM model are shown in Figure 5; the alumina concentration results of three models are shown in Table 1. To see the results more clearly, the actual alumina concentration values and simulation results are displayed under the same axis which is shown in Figure 6. There are five performance indicators to measure the quality of algorithm: Training Time, Testing Time, Training RMSE (root mean square error), Testing RMSE, and Average Relative Error (ARE), where (21)ARE=110∑i=110‍Ti-T^iTi×100%.The simulation results of five performance indicators are shown in Table 2.Table 1
Simulation alumina concentration results of three models.
1
2
3
4
5
6
7
8
9
10
Actual value (%)
1.94
2.37
2.03
1.95
2.16
2.02
1.99
1.94
2.64
2.94
Output of ELM (%)
1.8164
2.0641
2.0739
2.1544
1.9423
1.9065
1.9166
2.047
2.6761
2.6443
Output of BP (%)
1.7551
2.015
1.9953
2.4558
1.9594
1.9047
1.9510
2.0183
3.0956
2.9993
Output of SVM (%)
1.9176
2.003
2.006
2.0917
1.8828
1.8834
1.8979
1.8833
2.301
2.3042Table 2
The comparison of performance indicators.
Training Time (s)
Testing Time (s)
Training RMSE
Testing RMSE
ARE (%)
ELM
0
0
0.0926
0.1784
6.78
BP
28.3281
0.0156
0.1698
0.2626
9.17
SVM
0.0156
0
0.2635
0.2797
8.83Figure 5
Experimental aluminum reduction cell and anode rod.Figure 6
Sampled electrolyte.From Figure7, the horizontal axis denotes the values of voltage between anode rod and cathode bar and the vertical axis denotes the alumina concentration below the experimental anode rode. The simulation results of ELM approach the sample data in a certain range. As is shown in Table 1 and Figure 8, the values of output of ELM model are closer to the actual values of alumina concentration. From Table 2, Training Time of ELM model is 0 s whose Training RMSE is 0.0926% and Testing RMSE is 0.1784%, while Training Time of BP model is 28.3281 s whose Training RMSE is 0.1698% and Testing RMSE is 0.2626%. In the SVM model, Training Time is 0.0156 s and Training RMSE and Testing RMSE are 0.2635% and 0.2797%. In terms of Average Relative Error (ARE) performance indicator, ELM has the smallest ARE (6.78%) in all of the three algorithms. It is clear that SVM model has faster learning speed and less Average Relative Error than BP model, but SVM model is slower and less accurate than ELM model. So the soft sensing measurement model based on ELM algorithm has better performance than BP model and SVM model.Figure 7
Simulation result of ELM model.Figure 8
Actual alumina concentration values and simulation results.


## 6. Conclusions

Alumina concentration is very important in the aluminum electrolysis, which may affect the performance of aluminum reduction cell. It is difficult to measure the alumina content due to the complicated environment of aluminum reduction cell. This paper proposes a novel soft sensing method of alumina concentration in the electrolyte based on extreme learning machine (ELM) and builds the relationship between alumina concentration and voltage between anode rod and cathode bar. Through the simulation results and comparison of BP model and SVM model, we can see the validity and advantage of the proposed method. This method is able to effectively achieve rapid and reliable estimation of alumina concentration in a relatively short time.



---

*Source: 268132-2015-05-27.xml*