# Applying Softcomputing for Copper Recovery in Leaching Process

Authors: Claudio Leiva; Víctor Flores; Felipe Salgado; Diego Poblete; Claudio Acuña
Journal: Scientific Programming
 (2017)

---

## Abstract

The mining industry of the last few decades recognizes that it is more profitable to simulate model using historical data and available mining process knowledge rather than draw conclusions regarding future mine exploitation based on certain conditions. The variability of the composition of copper leach piles makes it unlikely to obtain high precision simulations using traditional statistical methods; however the same data collection favors the use of softcomputing techniques to enhance the accuracy of copper recovery via leaching by way of prediction models. In this paper, a predictive modeling contrasting is made; a linear model, a quadratic model, a cubic model, and a model based on the use of an artificial neural network (ANN) are presented. The model entries were obtained from operation data and data of piloting in columns. The ANN was constructed with 9 input variables, 6 hidden layers, and a neuron in the output layer corresponding to copper leaching prediction. The validation of the models was performed with real information and these results were used by a mining company in northern Chile to improve copper mining processes.


---

## Body

## 1. Introduction

Due to the complexity of extraction mining worldwide, computer models are becoming an essential tool for reducing production costs [1]. Computational techniques have recently been incorporated into the copper industry to improve both the process and the results obtained through leaching, thus achieving tangible results such as higher levels of production, identification of low grade ore, and reduction of production costs [2].Copper mining in Chile is the most profitable industry and contributes to approximately 12% of world copper production [3, 4]. In recent years, the mining industry in Chile has begun to utilize and show an increasing interest in Artificial Intelligence techniques in processes such as copper production predictions, where the traditional methods tend to lack certainty and knowledge, and one of these techniques is softcomputing [5].Softcomputing is the branch of Artificial Intelligence that groups paradigms and techniques working with incomplete and imprecise information in critical processes, in order for the company to obtain useful solutions for tasks such as prediction or discovery of information or knowledge [4–6]. A definition of softcomputing is“an emerging approach to computing which parallels that remarkable ability of the human mind to reason and learn in an environment of uncertainty and imprecision”[7].This paper describes the work done with ANN to generate behavior predictions of the pile in the copper leaching domain. The use of ANNs has been used in combination with other techniques commonly used as a predictor in the mining industry, such as linear regression, with the aim of contrasting and enhancing the quality of results. An ANN can be defined as a set of computational units (neurons) that are highly interconnected. Each neuron is also called node which represents the biological neuron, and the connection among them represents the biological neuronal network [8]. The basic element of data processing in a neural network is the perceptron, a simple element that with an input vector provides a single output, for which it uses an activation function. An extended and well-accepted way to classify ANN is according to architecture and learning. Architecture refers to the arrangement and connections between neurons, number of input and output variables, and number of hidden layers, while learning classification refers to training the network by way of patterns, and the patterns are performed iteratively until the restriction Od-Oi≤δ is satisfied, Od being the desired output, Oi the output of the nth iteration, and δ a small numerical value [9, 10].This paper describes the copper recovery prediction process for an extraction mine in northern Chile, using both a statistical model (methods traditionally used in the copper leaching industry) and a neural model known as multilayer perceptron with backpropagation [11, 12]. The multilayer perceptron model is considered the most widely used neural network due to its efficiency and ease of understanding and interpretation of both the network and its results. The parameters of entry to the ANN were selected with respect to the objectives and conditions of the company. Specifically, this paper reflects the processes and results obtained in the SCM Franke Company, of the international group KGHM International, which has been operating in Chile since 2009.


## 2. Background and Related Works

The process of applying different techniques to adjust the weights related to each of the input variables of neurons in an ANN is known as learning [8, 13]. For this process, which is also known as network training, several techniques can be followed but the most common ones are as follows: supervised and unsupervised training [14]. The supervised training consists of using a series of pairs (Ii,Oi) where Iirepresents an input data vector and Oi represents the desired output vector. The training consists of using an algorithm to make the values of Ii as close as possible to the vector Oi, while in the unsupervised training only the vector Ii is available and the training algorithm attempts to find the hidden structures between these data to adjust the weights of each neuron in the network [15]. In both cases, the learning rules act by modifying the weights to achieve the goal, for which there are several techniques as well. One of the most used techniques is error correction via adjustment of weights that essentially consists of obtaining a value δ very close to 0.Currently, the applications that include ANN are in all areas of science and engineering, that is, energy production estimation systems [16]. In engineering in particular, the increase in ANN usage in recent decades is significant, being applied to tasks such as prediction and regression. For example, in [17] a comparison between ANN and regression analysis is described, in terms of notation and implementation of both paradigms, and also highlights the advantages in terms of implementation. In [4, 18] works are presented related to metal-mechanical industry that compares the performance of the Bayesian networks with ANNs in the prediction of surface quality in high-speed machining processes. The most frequently used techniques for this process are ANN [13] and linear or multiple regression methods [14]. In [19], ANNs are used to predict the influence of biological and nonbiological parameters on the precipitation of ferric iron through a bioleaching process. In [20] a process of construction of a three-layer backpropagation ANN to predict the concentration of heavy metals (Zn, Ni, Cd, and Pb) as waste in a zinc leach cell is described.In [21] an ANN used to predict the copper production process is presented as well. Moreover, in the work recently reported in [22] an ANN is used to predict copper flotation rates under different operating conditions; different dosages of chemical reagents used in the process, feed rate, and granulometry are used for this process. A three-layer backpropagation neural network (input layer, two hidden layers, and one output), topology (9-10-10-3), was used, and the quality of the prediction in the testing process was 93%.


## 3. Materials and Methods

The SCM Franke Company uses three industrial processes widely known in the copper industry to produce metallic copper via hydrometallurgy: leaching in dynamic cells, solvent extraction, and electrowinning. The goal in these processes is to achieve the highest copper production by saving resources and having the lowest possible environmental impact. In this sense, the company considered carrying out simulations to predict the copper production through leaching processes. The objective was to produce an efficient model (with an accuracy greater than 95%) for the estimation of copper leaching recovery based on historical data.A previous analysis by the SCM Franke Company conducted to determine which of these processes is the least controlled (number of influencing factors and homogeneity of irrigation) found it to be the leaching process in dynamic stacks. Based on this previous analysis, the parameters to be used in the prediction models and the desired quality level (95%) were identified and employed as the desired adjustment value in the simulation to test the results of the models.The study summarized the prediction of copper extraction in dynamic stacks, applying mathematical models (statistical model) that require complete information to give precise results and an ANN, which is considered to have advantages over the treatment of incomplete information in terms of generating predictions in industrial production.Pile leaching copper is a percolation process that operates above ground. The procedure is illustrated in Figure1. The oxide copper ore is piled up on leach pads which have an inclination of 3° approx. and a rubber lining that seals the ground beneath. These heaps are 3-meter-high and 72 meters in length with a base area of 2,880 square meters, covering 120,000 tons of ore.Figure 1
Process flow diagram of the pile leaching process.A sprinkler system is installed on top of each heap allowing diluted sulfuric acid to be fed uniformly over the ore. While the solution percolates through the pile, copper is leached out of the ore. The pregnant leaching solution is collected by a drainage system below the pile and is led through a collection ditch into a pond.


## 4. Methodology

The research consists of two stages, the first called “copper recovery modeling” in which models are generated with each of the aforementioned methods and the second stage called “Evaluation” where the results are compared and the quantity and quality of information are obtained.Both stages are beneficial for the SCM Frank Company in order to understand the leaching process with the characteristics of the ore and to be able to predict the recovery of copper in dynamic stacks, focused on achieving a less than 5% margin of error in estimation. These two stages are possible due to the previous analysis and selection of variables and parameters representing the inputs to the prediction methods; Figure2 graphically shows this work scheme.Figure 2
Work scheme for copper recovery prediction in dynamic stacks.


### 4.1. Statistical Models

In order to perform the modeling, the most recent literature in leaching process was considered, taking into account the variables that affect the recovery of copper. In addition, the historical operational results data and process-related pilot testing were taken into account, which we call pilot data. Pilot tests consist of testing columns with strict control measures regarding irrigation rates, acid concentrations in irrigation solutions, and operating cycles, conditions that vary depending on the test to be performed. A case study was generated with a database of approximately 30,000 pieces of data.


#### 4.1.1. Data Collection

For the purpose of this study, the historical database corresponding with industrial and piloting performance from SCM Franke Company was used. The data stored in the database correspond to both plant data operation (weighted values and accumulated daily) and pilot data.The operating plant data were obtained at a frequency of 4 hours for 1 year. For some periods, the irrigation was stopped on some cells or modules in service; due to these periods, inconsistent results were disregarded and corresponding data was not considered for the data collection process; some “noise” in the system and useless information were disregarded as well. Regarding the pilot data, the method of information collection was the same used for the data of operation.The next step was to identify the parameters that affect the copper recovery and to form a robust database with this information for use in the preparation and evaluation of process models using different techniques in order to determine one that meets the plant requirements (adjustment greater than 95%) with respect to plant and pilot operating data. The parameters considered in the statistical model are detailed in Table1.Table 1
Entry parameters to the statistical model.
Notation
Name and description
Optimum values
x
1
Monoclass granulometry (refers to the mineral)
Between 11,5 mm and 15 mm
x
2
Irrigation rates
Between 14 (l×h)/m2 and 6 (l×h)/m2
x
3
Total acid added
Between 0,5 g/l and 100 g/l
x
4
Pile high
Between 1 m and 5 m
x
5
Total copper grade
Between 0,5% and 2%
x
6
CO3 grade
Between 0,5% and 10%
x
7
Leaching ratio
Less than 15 m3/TMS in pilot and Less than 8 m3/TMS in plant
x
8
Operation days
Between 90 and 120 days
x
9
Soluble copper grade stacked
70% of the total copper gradeA search was performed to find the combination of suitable variables to produce the most successful model for the desired response. The system response was defined as the percentage of copper extraction. To validate the model, the response variable was monitored. The model was monitored in order to obtain the best fit for the operating conditions and variations. When model misalignment was below acceptable (much lower than the desired 95%), the model was readjusted and remade considering the operating ranges that were initially not considered; parameters that were initially not considered were now carried out.To select the subset of variables that ensure a practical model, coefficients were generated using the Minitab software. In detail, the calculated coefficients were as follows: Vars,R-cuad, adjusted R-cuad, Cp of Mallows, and S (standard dispersion). A preliminary statistical model was made considering the Press factor, used to avoid overadjustment. Additionally, for this model the following restrictions were considered: (1) that the adjusted values of S, R-cuad (R2), and R-squared were the highest possible and (2) that the value of the Press factor was the highest grade possible.Originating from the preliminary model and utilizing the Minitab tools, three statistical models were generated: linear adjustment model, quadratic adjustment model, and cubic adjustment model. Table2 summarizes each of these models and the adjustments results according to the value R2.Table 2
Comparative ofR2 values.
Model
R
2 value
Lineal model
69,8%
Quadratic model (standard deviation = 9,4)
89,9%
Cubic model (standard deviation = 8,3)
92,3%
ANN model
97,9%For the lineal model (1),(1)Y=∑i=19αi∗χi+C.For the quadratic and cubic models (2)(2)Y=∑i=19αi∗χi+∑j=19∑i=19bi∗χi∗χj+C,Y=∑i=19ai∗χi+∑j=19∑i=19bi∗χi∗χj+∑k=19∑j=19∑i=19ci∗χi∗χj∗χk+C,where(i)
χi, χj, and χk represent the input value of the variable;(ii)
ai, bi, and ci represent the weight of the input;(iii)
C represents a constant


### 4.2. ANN Model Base

To perform the RNA modeling, the MATLAB program was applied. For the network programming, the following parameters were used: 9 income variables (see Table1), hidden layers, 60% of the information for the network, 40% of the information for network validation, 0.2 of error as cycle-error, or maximum of 500 iterations. The data for the training and validation was randomly selected.In the variables selection the following was considered: the experience accumulated by operation experts of the SCM Franke Company, the results of the piloting, bibliography, and the results observed with the statistical models described above. In the ANN configuration, nine input variables were used (see Table1) for a multilayer perceptron model with backpropagation algorithm.The output layer was formed by a corresponding neuron to the prediction of copper recovery in dynamic piles, equally focusing on achieving a less than 5% error in the prediction estimation. The mathematical expression used in each neuron to obtain the output value is indicated in (3) where (xi) represents the input value to the neuron and the variable (wij) represents the weight associated with the neuron. The task of initializing the weights associated with each input variable was performed randomly in a range normalized between (-1,1) and using the MAPMINMAX function of MATLAB (a function used to normalize the input and output values of the neurons), and threshold values that were used in the activation functions (sigmoid function) were initialized in a similar way (these data were not reflected in this document because they are subject to confidentiality agreements with the company). The mathematical expression of the sigmoid function is indicated in (3)y=∑j=19Lj∑i=99wiji∗xi,where(i)
xi represents the ith input value to the neuron;(ii)
wiji represents the ith weight associated with the neuron i of the layer j;(iii)
Lj represents the constant of the jth layer.


## 4.1. Statistical Models

In order to perform the modeling, the most recent literature in leaching process was considered, taking into account the variables that affect the recovery of copper. In addition, the historical operational results data and process-related pilot testing were taken into account, which we call pilot data. Pilot tests consist of testing columns with strict control measures regarding irrigation rates, acid concentrations in irrigation solutions, and operating cycles, conditions that vary depending on the test to be performed. A case study was generated with a database of approximately 30,000 pieces of data.


### 4.1.1. Data Collection

For the purpose of this study, the historical database corresponding with industrial and piloting performance from SCM Franke Company was used. The data stored in the database correspond to both plant data operation (weighted values and accumulated daily) and pilot data.The operating plant data were obtained at a frequency of 4 hours for 1 year. For some periods, the irrigation was stopped on some cells or modules in service; due to these periods, inconsistent results were disregarded and corresponding data was not considered for the data collection process; some “noise” in the system and useless information were disregarded as well. Regarding the pilot data, the method of information collection was the same used for the data of operation.The next step was to identify the parameters that affect the copper recovery and to form a robust database with this information for use in the preparation and evaluation of process models using different techniques in order to determine one that meets the plant requirements (adjustment greater than 95%) with respect to plant and pilot operating data. The parameters considered in the statistical model are detailed in Table1.Table 1
Entry parameters to the statistical model.
Notation
Name and description
Optimum values
x
1
Monoclass granulometry (refers to the mineral)
Between 11,5 mm and 15 mm
x
2
Irrigation rates
Between 14 (l×h)/m2 and 6 (l×h)/m2
x
3
Total acid added
Between 0,5 g/l and 100 g/l
x
4
Pile high
Between 1 m and 5 m
x
5
Total copper grade
Between 0,5% and 2%
x
6
CO3 grade
Between 0,5% and 10%
x
7
Leaching ratio
Less than 15 m3/TMS in pilot and Less than 8 m3/TMS in plant
x
8
Operation days
Between 90 and 120 days
x
9
Soluble copper grade stacked
70% of the total copper gradeA search was performed to find the combination of suitable variables to produce the most successful model for the desired response. The system response was defined as the percentage of copper extraction. To validate the model, the response variable was monitored. The model was monitored in order to obtain the best fit for the operating conditions and variations. When model misalignment was below acceptable (much lower than the desired 95%), the model was readjusted and remade considering the operating ranges that were initially not considered; parameters that were initially not considered were now carried out.To select the subset of variables that ensure a practical model, coefficients were generated using the Minitab software. In detail, the calculated coefficients were as follows: Vars,R-cuad, adjusted R-cuad, Cp of Mallows, and S (standard dispersion). A preliminary statistical model was made considering the Press factor, used to avoid overadjustment. Additionally, for this model the following restrictions were considered: (1) that the adjusted values of S, R-cuad (R2), and R-squared were the highest possible and (2) that the value of the Press factor was the highest grade possible.Originating from the preliminary model and utilizing the Minitab tools, three statistical models were generated: linear adjustment model, quadratic adjustment model, and cubic adjustment model. Table2 summarizes each of these models and the adjustments results according to the value R2.Table 2
Comparative ofR2 values.
Model
R
2 value
Lineal model
69,8%
Quadratic model (standard deviation = 9,4)
89,9%
Cubic model (standard deviation = 8,3)
92,3%
ANN model
97,9%For the lineal model (1),(1)Y=∑i=19αi∗χi+C.For the quadratic and cubic models (2)(2)Y=∑i=19αi∗χi+∑j=19∑i=19bi∗χi∗χj+C,Y=∑i=19ai∗χi+∑j=19∑i=19bi∗χi∗χj+∑k=19∑j=19∑i=19ci∗χi∗χj∗χk+C,where(i)
χi, χj, and χk represent the input value of the variable;(ii)
ai, bi, and ci represent the weight of the input;(iii)
C represents a constant


## 4.1.1. Data Collection

For the purpose of this study, the historical database corresponding with industrial and piloting performance from SCM Franke Company was used. The data stored in the database correspond to both plant data operation (weighted values and accumulated daily) and pilot data.The operating plant data were obtained at a frequency of 4 hours for 1 year. For some periods, the irrigation was stopped on some cells or modules in service; due to these periods, inconsistent results were disregarded and corresponding data was not considered for the data collection process; some “noise” in the system and useless information were disregarded as well. Regarding the pilot data, the method of information collection was the same used for the data of operation.The next step was to identify the parameters that affect the copper recovery and to form a robust database with this information for use in the preparation and evaluation of process models using different techniques in order to determine one that meets the plant requirements (adjustment greater than 95%) with respect to plant and pilot operating data. The parameters considered in the statistical model are detailed in Table1.Table 1
Entry parameters to the statistical model.
Notation
Name and description
Optimum values
x
1
Monoclass granulometry (refers to the mineral)
Between 11,5 mm and 15 mm
x
2
Irrigation rates
Between 14 (l×h)/m2 and 6 (l×h)/m2
x
3
Total acid added
Between 0,5 g/l and 100 g/l
x
4
Pile high
Between 1 m and 5 m
x
5
Total copper grade
Between 0,5% and 2%
x
6
CO3 grade
Between 0,5% and 10%
x
7
Leaching ratio
Less than 15 m3/TMS in pilot and Less than 8 m3/TMS in plant
x
8
Operation days
Between 90 and 120 days
x
9
Soluble copper grade stacked
70% of the total copper gradeA search was performed to find the combination of suitable variables to produce the most successful model for the desired response. The system response was defined as the percentage of copper extraction. To validate the model, the response variable was monitored. The model was monitored in order to obtain the best fit for the operating conditions and variations. When model misalignment was below acceptable (much lower than the desired 95%), the model was readjusted and remade considering the operating ranges that were initially not considered; parameters that were initially not considered were now carried out.To select the subset of variables that ensure a practical model, coefficients were generated using the Minitab software. In detail, the calculated coefficients were as follows: Vars,R-cuad, adjusted R-cuad, Cp of Mallows, and S (standard dispersion). A preliminary statistical model was made considering the Press factor, used to avoid overadjustment. Additionally, for this model the following restrictions were considered: (1) that the adjusted values of S, R-cuad (R2), and R-squared were the highest possible and (2) that the value of the Press factor was the highest grade possible.Originating from the preliminary model and utilizing the Minitab tools, three statistical models were generated: linear adjustment model, quadratic adjustment model, and cubic adjustment model. Table2 summarizes each of these models and the adjustments results according to the value R2.Table 2
Comparative ofR2 values.
Model
R
2 value
Lineal model
69,8%
Quadratic model (standard deviation = 9,4)
89,9%
Cubic model (standard deviation = 8,3)
92,3%
ANN model
97,9%For the lineal model (1),(1)Y=∑i=19αi∗χi+C.For the quadratic and cubic models (2)(2)Y=∑i=19αi∗χi+∑j=19∑i=19bi∗χi∗χj+C,Y=∑i=19ai∗χi+∑j=19∑i=19bi∗χi∗χj+∑k=19∑j=19∑i=19ci∗χi∗χj∗χk+C,where(i)
χi, χj, and χk represent the input value of the variable;(ii)
ai, bi, and ci represent the weight of the input;(iii)
C represents a constant


## 4.2. ANN Model Base

To perform the RNA modeling, the MATLAB program was applied. For the network programming, the following parameters were used: 9 income variables (see Table1), hidden layers, 60% of the information for the network, 40% of the information for network validation, 0.2 of error as cycle-error, or maximum of 500 iterations. The data for the training and validation was randomly selected.In the variables selection the following was considered: the experience accumulated by operation experts of the SCM Franke Company, the results of the piloting, bibliography, and the results observed with the statistical models described above. In the ANN configuration, nine input variables were used (see Table1) for a multilayer perceptron model with backpropagation algorithm.The output layer was formed by a corresponding neuron to the prediction of copper recovery in dynamic piles, equally focusing on achieving a less than 5% error in the prediction estimation. The mathematical expression used in each neuron to obtain the output value is indicated in (3) where (xi) represents the input value to the neuron and the variable (wij) represents the weight associated with the neuron. The task of initializing the weights associated with each input variable was performed randomly in a range normalized between (-1,1) and using the MAPMINMAX function of MATLAB (a function used to normalize the input and output values of the neurons), and threshold values that were used in the activation functions (sigmoid function) were initialized in a similar way (these data were not reflected in this document because they are subject to confidentiality agreements with the company). The mathematical expression of the sigmoid function is indicated in (3)y=∑j=19Lj∑i=99wiji∗xi,where(i)
xi represents the ith input value to the neuron;(ii)
wiji represents the ith weight associated with the neuron i of the layer j;(iii)
Lj represents the constant of the jth layer.


## 5. Results

Table2 shows the R2 values obtained by both, algorithmic models and ANN model; as seen in the table, the second-best value of R2 corresponds to the cubic fit model. None of these values exceed the desired minimum value (95%); however, the values obtained in the simulation, using cubic algorithmic models, are close to what was obtained.The best value corresponding to the ANN model during the learning phase and validation of the ANN, a 97,9% of adjustment, was obtained. Figure3 shows the adjustment performed with the ANN.Figure 3
ANN validation model.


## 6. Discussion

The result obtained through linear modeling was far from the required (89.9% v/s 95%), which indicates that the capacity to make the needed adjustment far exceeds what can be delivered by this type of modeling. The adjustment results of the quadratic model were approximately 22% better than those obtained previously and theR2 value was very similar to R2 (adjusted), which means the quadratic model is better able to generate more realistic data projections. With respect to the cubic model, the observed adjustment was close to what was obtained using the quadratic model, which indicates that the use of higher-grade models does not present a significant improvement in relation to the model adjustment, due to the limitations given by the techniques used and the high variability presented by the data used.The use of ANN for system modeling led to an improvement from 92.3% to 97.9%, compared to the cubic model. It was observed that the use of ANN achieved the objective of adjusting to fit, which was due to the complex network that this process uses to adjust the system response to the existing parameters.To achieve these results through the use of ANN, the use of a program such as MATLAB or another advanced calculation program is required. Due to the high complexity of such a modeling technique, it requires a database with an abundance of information so as to give the program the information necessary to evaluate as many cases as possible. On the other hand, ANN has the disadvantage of not being easily detected compared to statistical models, so this could present a problem for new users.


## 7. Conclusion

SCM Franke has carried out simulations on other copper production processes and knows from experience that the prediction is much cheaper than the experimental work, but, in this case, our study describes a prediction model that represents the behavior of the leaching plant, considering the variables initially defined with adjustment results higher than 95% obtained.Accordingly, it was determined that ANN was the best model for SCM Franke mining leach plant, due to the high variability of the existing plant results and tests. The adjustment obtained was 97.9%, which is higher than the adjustment of 95% initially requested.This study has served to obtain a comparison between prediction models and it can be intuited by the precision of the adjustments that the neural model has potential for use in future copper production prediction process. Furthermore, experience was gained in defining the model of an ANN, which can be used in future process simulations related to the improvement of copper attainment via softcomputing techniques in the SCM Franke Company or in similar companies.



---

*Source: 6459582-2017-12-05.xml*