Next Article in Journal
Environmental Application of Biogenic Magnetite Nanoparticles to Remediate Chromium(III/VI)-Contaminated Water
Next Article in Special Issue
The Role of Mineral Matter in Concentrating Uranium and Thorium in Coal and Combustion Residues from Power Plants in Poland
Previous Article in Journal
He-Ar Isotopes and Trace Gas Compositions of Fluid Inclusions in Massive Sulphides from the Yushui Copper-Polymetallic Deposit, South China: Metallogenic Implications
Previous Article in Special Issue
Nano-Scale Rare Earth Distribution in Fly Ash Derived from the Combustion of the Fire Clay Coal, Kentucky
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Threshold Value Determination Using Machine Learning Algorithms for Ba Interference with Eu in Coal and Coal Combustion Products by ICP-MS

1
College of Geoscience and Survey Engineering, China University of Mining and Technology, Beijing 100083, China
2
Department of Computing, Hong Kong Polytechnic University, Hung Hom, Kowloon, HKSAR, Hong Kong, China
*
Author to whom correspondence should be addressed.
Minerals 2019, 9(5), 259; https://doi.org/10.3390/min9050259
Submission received: 15 March 2019 / Revised: 24 April 2019 / Accepted: 24 April 2019 / Published: 29 April 2019

Abstract

:
Ba-based ion interference with Eu in coal and coal combustion products during quadrupole-based inductively coupled plasma mass spectrometry procedures is problematic. Thus, this paper proposes machine-learning-based prediction models for determination of the threshold value of Ba interference with Eu, which can be used to predict such interference in coal. The models are trained for Eu, Ba, Ba/Eu, and Ba interference with Eu. Under different user-defined parameters, different prediction models based on the corresponding model tree can be applied to Ba interference with Eu. We experimentally show the effectiveness of these different prediction models and find that, when the Ba/Eu value is less than 2950, the Ba-Eu interference prediction model is y = 0.18419411 + 0.00050737 × x ,   0 < x < 2950 . Further, when the Ba/Eu value is between 2950 and 189,523, the Ba-Eu interference prediction model of y   =   0.293982186   +   0.00000181729975   ×   x ,   2950   <   x   <   189 , 523 yields the best result. Based on the optimal model, a threshold value of 363 is proposed; i.e., when the Ba/Eu value is less than 363, Ba interference with Eu can be neglected during Eu data interpretation. Comparison of this threshold value with a value proposed in earlier works reveals that the proposed prediction model better determines the threshold value for Ba interference with Eu.

1. Introduction

Rare earth elements and yttrium (REY, or REE if Y is excluded) in coal and coal combustion products (CCPs), e.g., fly and bottom ash, have attracted much attention in recent years, not only because of the high international demand for these technologically important elements, but also because of the restrictions on export from China [1,2]. Seredin and Dai [3] and Dai et al. [4] have shown that coal has high potential as a REY source, given that the average concentration of REY oxides (REO) in world coal ash is 485 μg/g, which is half the cut-off grade of REO in CCPs (1000 μg/g). In some cases, CCPs contain >1000 μg/g REO; thus, they could constitute an economically viable source for REY extraction. Previous investigations have shown that some coals from China [5,6,7], Russia [3,8], and the USA contain high concentrations of REY [9,10,11], comparable to or even higher than those of conventional REY deposits [3]. Other studies concerning REY resources [12,13], modes of REY occurrence in coal and CCPs [14,15], and extraction technology [16] have also suggested the great potential of coal as REY source.
The REY (including Eu) concentration in coal and CCPs can be determined via several methods, including X-ray fluorescence spectrometry (XRF) [17,18], instrumental neutron activation analysis (INAA) [19,20,21], laser-ablation inductively-coupled plasma mass spectroscopy (La-ICP-MS) [22], sector-field inductively-coupled plasma mass spectroscopy ICP-MS [23], laser-induced break-down spectrometry [24,25,26], and quadrupole-based ICP-MS (ICP-Q-MS) [4,27]. Among these techniques, ICP-Q-MS has mostly been used for determination of REY concentrations in coal and CCPs. This is because of this method’s ability to rapidly, precisely, and accurately measure REY content at very low detection limits in liquid and solid samples, with relatively simple spectra and a wide linearity range [14,27,28,29]. However, in some cases, Eu concentrations in coal and CCPs cannot be accurately measured by the ICP-Q-MS. This difficulty is attributed to overlapping ion interference (such as those from M+, MO+, and/or MOH+ ions) with Eu, which affects interpretation of the ICP-Q-MS spectra [30,31,32,33,34,35]. Because Ba concentrations in coal and CCPs are generally much higher than those of Eu [4,36], the most significant types of interference encountered for Eu in coal and CCPs are due to 135Ba16O, 134Ba16OH, 137Ba16O, and/or 136Ba16OH interfering with 151Eu and/or 153Eu. For example, the ratio of Ba vs. Eu for world coals is as high as ~882 (with an average concentration of 150 μg/g Ba and 0.47 μg/g Eu) [36]. Therefore, as noted by Dai et al. [4], Eu content values in coal and CCPs determined via ICP-Q-MS should be treated with great caution.
Recently, Yan et al. [27] described a reliable analytical method to avoid Ba-based ion interference with Eu in coal, CCPs, and sedimentary rocks during ICP-Q-MS procedures, which was based on AG50W-X8 cation exchangeable chromatography. Yan et al. [27] provided an estimated Ba/Eu threshold value of 1000 and demonstrated that the determined Eu concentration should not be used if the Ba/Eu ratio exceeds this threshold value and if no effective action has been taken to avoid Ba interference with Eu. However, the Ba/Eu threshold value estimated by Yan et al. [27] is an approximation. In a different context, Loges et al. [35] also suggested an experience-based threshold value, i.e., Ba/Eu > 1000, to exclude interference of Ba with Eu in ICP-Q-MS analyses of Eu content values; however, this is not an accurate estimation either. Determination of an accurate Ba/Eu threshold value is important not only to determine if coal can potentially be used as a REY source, but also to deduce the geological setting of coal formation using Eu as a reliable indicator.
In this paper, we propose a threshold value at which Ba interference with Eu has a meaningful effect on ICP-Q-MS results, which is determined via prediction models created using machine learning algorithms. All Eu and Ba training data employed in this study are from Yan et al. [27]. Yan’s work [27] has significantly diminished the inference of 137Ba16O, 136Ba17O, 135Ba18O, and 134Ba18OH on 153Eu in related samples. They used a Bio-Rad AG50W-x8 cation exchange resin to effectively separate Ba from digested solutions of related coal and CCP samples. The results as presented by Yan et al. [27] showed that the determination of Ba and Eu in the National Institute of Standards Technology (NIST) standard references of coal and fly ash samples using this method is quite reliable. Also, Yan et al. [27] showed that the detection limits for Ba and Eu are very low, 0.030 μg/L and 0.006 μg/L, respectively, and the determination coefficient of their calibration curves is >0.9999. Unlike other ICP-Q-MS data which were obtained based on non-separation between Ba and Eu in solutions digested from solid samples (e.g., coal samples in the U.S. Geological Survey’s WoCQI database, Palmer et al. [37], and in other numerous published papers, for example but not limited to references [38,39,40,41,42,43,44,45,46]), the data by Yan et al. [27] provided a good opportunity for determining the threshold value using machine learning algorithms for Ba interference with Eu in coal and coal combustion products by ICP-Q-MS. However, there have been some studies to assess the interference of 135Ba16O, 134Ba16OH, 137Ba16O, and/or 136Ba16OH on 151Eu and/or 153Eu, e.g., determining the yield of Ba-based oxide and hydroxide ions using a single-element solution of Ba (e.g., 500 ng/mL Ba in BaCl2 solution) and compare the yield of potential interfering irons (e.g., 135Ba16O) with the Eu ions (151Eu) that has the same mass number in a single-element solution of Eu [47,48,49]. A study by Dulski [48] shows that that 1000 ng/g Ba could cause 0.22 ng/g Eu. Another approach to evaluating the degree of the interference is to analyze the correlation relation between Ba and Eu in the related samples, i.e., a linear Ba-Eu correlation indicating distinct interference of Ba with Eu [50,51,52,53,54,55,56].
Classification and regression are two typical algorithms in machine learning [57], with the difference between them being that their target variables are discrete and continuous, respectively [58]. Here, we employ a model tree [57] based on linear regression and a regression tree to construct prediction models for this interference. Based on analysis of Ba, Eu, and Ba/Eu, the target variables for Ba interference with Eu are found to be continuous. Thus, we adopt the regression method for prediction of Ba interference with Eu. Empirically, we find that it is difficult to precisely construct a global prediction model using linear regression because of the element data complexity. To overcome this problem, the models are adjusted to split the element data according to many partitions. In the case of the regression tree, classification and regression tree (CART) algorithms [59] are then applied.

2. Proposed Machine Learning Models for Prediction of Ba Interference with Eu

In this study, machine learning algorithms were used to develop models of Ba interference with Eu in coal. Note that, when a pair of element datasets (e.g., ( x 1 Ba , , x i Ba ) and ( x 1 Eu , , x i Eu ) for Ba and Eu, respectively) is established, the interference between them is difficult to determine. Effective Ba-Eu interference prediction depends on various factors including the element concentration, element interference, and samples.
Here, three machine learning models were used to predict Ba interference with Eu, incorporating linear regression, regression trees, and model trees. The problem of threshold value identification can be represented as the problem of constructing a prediction model between a Ba interference with Eu dataset ( y 1 prediction Ba / Eu , , y i prediction Ba / Eu ) and Ba, Eu, and Ba/Eu ratio datasets ( x 1 Ba , , x i Ba ) , ( x 1 Eu , , x i Eu ) and ( x 1 Ba / x 1 Eu , , x i Ba / x i Eu ) respectively.

2.1. Linear Regression Model for Prediction of Ba Interference with Eu

2.1.1. Linear Regression Model

The ( x 1 Ba , , x i Ba ) and ( x 1 Eu , , x i Eu ) datasets used in this study were based on a number of coal and ash samples. Furthermore, based on experiments by Yan et al. [27], we calculated ( x 1 Ba / x 1 Eu , , x i Ba / x i Eu ) . The linear regression prediction model for Ba-Eu interference is:
y i prediction Ba / Eu = ( x i Ba / x i Eu ) T w i prediction .
Specifically, the target variable of Ba-Eu interference y i prediction Ba / Eu is related to the different element concentrations, i.e., ( x 1 Ba , , x i Ba ) , ( x 1 Eu , , x i Eu ) and ( x 1 Ba / x 1 Eu , , x i Ba / x i Eu ) . Linear regression for Ba-Eu interference predicts target interference values. In the model given in Equation (1), the vector w i prediction is the regression weight. Regression is used to find w i prediction and hence, the Ba-Eu interference values are predicted.

2.1.2. Ba-Eu Interference Prediction Error

The error is defined as the difference between the actual Ba interference with Eu y i Ba / Eu and y i prediction Ba / Eu ; i.e.,:
i = 1 m ( ( x i Ba / x i Eu ) T w i prediction y i Ba / Eu ) 2 .

2.1.3. Machine Learning Process for Ba-Eu Interference Prediction

Training: First, all input concentrations of Ba, Eu, and the Ba/Eu ratio values (i.e., ( x 1 Ba , , x i Ba ) , ( x 1 Eu , , x i Eu ) and ( x 1 Ba / x 1 Eu , , x i Ba / x i Eu ) ), and the interference values of Ba with Eu ( y 1 Ba / Eu , , y i Ba / Eu ) are entered. All the input training element data values are prepared and converted into matrices. The interference error of the above elements can also be expressed in matrix notation, as:
w i prediction = ( ( x i Ba / x i Eu ) T ( x i Ba / x i Eu ) ) 1 ( x i Ba / x i Eu ) T y i Ba / Eu
This equation is solved using the ordinary least squares method. Hence, w i prediction is predicted according to the best estimate based on the training element data values.
Prediction: Based on the training process described above, the prediction for Ba-Eu interference can be estimated from the formula:
y i prediction Ba / Eu = ( x i Ba / x i Eu ) T w i prediction .

2.2. Regression Tree Model for Prediction of Ba Interference with Eu

Regression Tree Model

The elements have many features, i.e., x i Ba , x i Eu , x i Ba / x i Eu and y i Ba / Eu , and samples C i Ba and C i Eu . The linear regression model cannot achieve good fitting of ( x 1 Ba , , x i Ba ) , ( x 1 Eu , , x i Eu ) and ( x 1 Ba / x 1 Eu , , x i Ba / x i Eu ) with ( y 1 prediction Ba / Eu , , y i prediction Ba / Eu ) , as shown in Figure 1; thus, a nonlinear model emerges. The nonlinear model partitions the element concentrations and Ba-Eu interference. Every partition can be constructed with linear regression models. Note that in Figure 1, the x-axis is the ratio of Ba vs. Eu in digested solutions derived from solid samples before Ba is separated from Eu in the solutions; and y-axis is the Eu concentration contributed from Ba ions.
The steps of the CART algorithm for constructing the regression tree model for prediction of Ba interference with Eu are as follows. The extracted feature is Ba/Eu, i.e., ( x 1 Ba / x 1 Eu ) , , ( x i Ba / x i Eu ) . The extracted feature value is the interference of Ba with Eu, i.e., y 1 Ba / Eu , , y i Ba / Eu . For every ( x 1 Ba / x 1 Eu ) , , ( x i Ba / x i Eu ) , binary splits are executed to yield two parts ( x 1 Ba / x 1 Eu ) , , ( x j Ba / x j Eu ) < s , ( x j Ba / x j Eu ) , , ( x i Ba / x i Eu ) > s . Then, for every feature in the two different parts, y 1 prediction Ba / Eu , , y i prediction Ba / Eu is calculated based on Equation (5). The process surveys every feature and value to find the best split that minimizes the error:
min { min x i Ba / Eu s ( y i prediction Ba / Eu c 1 ) 2 + min ( x i Ba / Eu > s ( y i prediction Ba / Eu c 2 ) 2 ) } ; c 1 = 1 N x i Ba / Eu s y i Ba / Eu ,   c 2 = 1 M x i Ba / Eu > s y i Ba / Eu ,
where N and M are the feature value numbers in the different parts.
Based on the binary split process described above, for every feature ( x 1 Ba / x 1 Eu ) , , ( x i Ba / x i Eu ) , if the feature value y 1 prediction Ba / Eu , , y i prediction Ba / Eu is greater than the best split value, we traverse the left side of the regression tree, i.e., the left subtree T r e e l e f t . If the feature value is lower than the best split value, we traverse the right side of the regression tree, i.e., T r e e r i g h t . For T r e e l e f t and T r e e r i g h t , we survey every feature and feature value to find the best split until the minimum error is achieved (cf. Equation (5)). The binary split process is recursive for many iterations until the feature cannot be split; then, its feature value is the leaf node. Hence, the Ba interference with the Eu regression trees can be determined.

2.3. Model Tree

The model tree for prediction of Ba interference with Eu is based on the linear regression and regression tree models described above. The steps of the CART algorithm for this model tree are similar to those for the regression tree.
The feature and feature values extracted here are ( x 1 Ba / x 1 Eu ) , , ( x i Ba / x i Eu ) and y 1 Ba / Eu , , y i Ba / Eu , respectively. For every ( x 1 Ba / x 1 Eu ) , , ( x i Ba / x i Eu ) , we execute binary splits that yield two parts ( x 1 Ba / x 1 Eu ) , , ( x j Ba / x j Eu ) < s , ( x j Ba / x j Eu ) , , ( x i Ba / x i Eu ) > s . Then, for every feature in the two different parts, we calculate y 1 prediction Ba / Eu , , y i prediction Ba / Eu based on Equation (6). We repeat the process for every feature and every value to find the best split that minimizes the error; i.e.,:
min { min x i Ba / Eu s ( ( x i Ba / x i Eu ) T w i prediction c 1 ) 2 + min ( x i Ba / Eu > s ( ( x i Ba / x i Eu ) T w i prediction c 2 ) 2 ) } c 1 = 1 N x i Ba / Eu s y i Ba / Eu ;       c 2 = 1 M x i Ba / Eu > s y i Ba / Eu y i prediction Ba / Eu = ( x i Ba / x i Eu ) T w i prediction .
Based on the binary split process above, for every feature ( x 1 Ba / x 1 Eu ) , , ( x i Ba / x i Eu ) , if the feature value { ( x 1 Ba / x 1 Eu ) T w 1 prediction , , ( x i Ba / x i Eu ) T w i prediction } is greater than the best split value, we traverse T r e e l e f t . If the feature value is lower than the best split value, we traverse T r e e r i g h t . For T r e e l e f t and T r e e r i g h t , we survey every feature and feature value to find the best split until the minimum error is achieved (cf. Equation (6)). The binary split process is recursive for many iterations until the feature cannot be split; then, its feature value is the leaf node. Hence, model trees for Ba interference with Eu can be formed. The difference between the regression tree and model tree is that the leaf nodes of the regression tree are constant sets with y 1 prediction Ba / Eu , , y i prediction Ba / Eu , but the leaf nodes of the model tree are linear model sets with { ( x 1 Ba / x 1 Eu ) T w 1 prediction , , ( x i Ba / x i Eu ) T w i prediction } .

2.4. Machine Learning Process for Ba-Eu Interference Prediction

Based on the constructed regression tree and model tree for Ba-Eu interference prediction, the proposed machine learning process is implemented as follows.
Training: All ( x 1 Ba , , x i Ba ) , ( x 1 Eu , , x i Eu ) , ( x 1 Ba / x 1 Eu , , x i Ba / x i Eu ) , and ( y 1 Ba / Eu , , y i Ba / Eu ) are entered.
Regression tree for prediction: After regression tree training, we perform binary splits to obtain P parts recursively, and obtain a prediction of the Ba-Eu interference y i prediction Ba / Eu =   c p I ( x i Ba / Eu P ) based on determination of the best split that minimizes the error:
P min { min x i Ba / Eu s ( y i prediction Ba / Eu c 1 ) 2 + min ( x i Ba / Eu > s ( y i prediction Ba / Eu c 2 ) 2 ) } , c 1 = 1 N x i Ba / Eu s y i Ba / Eu ,   c 2 = 1 M x i Ba / Eu > s y i Ba / Eu .
Model tree for prediction: After the Ba-Eu interference model tree training, we execute binary splits to obtain P parts recursively, and obtain y i prediction Ba / Eu =   c p I ( x i Ba / Eu P ) based on determination of the best split that minimizes the error:
p min { min x i Ba / Eu s ( ( x i Ba / x i Eu ) T w i prediction c 1 ) 2 + min ( x i Ba / Eu > s ( ( x i Ba / x i Eu ) T w i prediction c 2 ) 2 ) } , y i prediction Ba / Eu = ( x i Ba / x i Eu ) T w i prediction .
From the above analysis, the model tree for prediction of Ba interference with Eu is selected.

3. Results and Performance Evaluation

3.1. Simulation Setup

To implement the Ba-Eu interference prediction models and calculate the Ba/Eu threshold value, the Python programming language [60] was used. The prediction models for the Ba interference with Eu were constructed as follows:
(1) All relevant element data were collected, as detailed in Table 1 and Table 2.
(2) All input element data were prepared. Note that all Ba, Eu and Ba/Eu concentrations were prepared in a standard Python format list.
(3) The element data were analyzed. Note that all input element data could be analyzed for feature selection. These features included x i Ba , x i Eu , x i Ba / x i Eu , y i Ba / Eu , and C i Ba , C i Eu .
(4) The algorithm was trained. To achieve the target variable, i.e., the Ba interference with Eu, and the threshold points of the element, i.e., Ba/Eu, we implemented the model tree based on linear regression and the regression tree in Python.
(5) The algorithm was tested. That is, the performance of the interference prediction model obtained in the above step was tested.

3.2. Model Tree for Prediction of Ba Interference with Eu

The input element data sets were prepared as detailed in Table 3. For execution of the model tree for element interference prediction in Python, two variables were necessary: tolS and tolN, the tolerance of the Ba-Eu interference error reduction and the minimum Ba related to Eu data instances in a split, respectively. Note that the model tree was sensitive to the tolS and tolN settings. and different settings yielded different prediction models.
We performed model tree experiments for element interference prediction by inputting rare earth element data sets. All prediction models are detailed in Table 4. For (tolS,tolN) = (0,1) and (0,2), the prediction model for Ba interference with Eu is shown in Figure 2a. There are six split values: (1) when the Ba/Eu value is greater than 40,726.5, the Ba-Eu interference prediction model is a linear regression, where y = 0.389451044 + 0.00000121646954 × x , x > 40 , 726.5 ; (2) when the Ba/Eu value is less than 40,726.5 and greater than 17,344.79, the interference prediction model is a linear regression, with y = 0.759759680 0.0000125166582 × x , 17 , 344.79 < x < 40 , 726.5 ; (3) when the Ba/Eu value is less than 17,344.79 and greater than 2950, the interference prediction model is a linear regression, where y = 0.0364849362 + 0.0000332367781 × x , 2950 < x < 17 , 344.79 ; (4) when the Ba/Eu value is less than 2950 and greater than 2247.06, the interference prediction model is a linear regression, with y = 1.29 + 0.0002 × x , 2247.06 < x < 2950 ; (5) when the Ba/Eu value is less than 2247.06 and greater than 10.88, the interference prediction model is a linear regression, with y = 0.0131201121 + 0.0000291815483 × x , 10.88 < x < 2247.06 ; and (6) when the Ba/Eu value is less than 10.88, the interference prediction model is a linear regression, where y = 0.06066667 + 0.00833333 × x , x < 10.88 .
For (tolS, tolN) = (0,3), the model trees for prediction of Ba interference with Eu are shown in Figure 2b, having four split values: (1) when the Ba/Eu value is greater than 27,144.6, the interference prediction model is a linear regression, with y = 0.259761293 + 0.00000193096467 × x , x > 27 , 144.6 ; (2) when the Ba/Eu value is greater than 2950 and less than 27,144.6, the interference prediction model is a linear regression, with y = 0.0966488103 + 0.0000156280402 × x , 2950.0 < x < 27 , 144.6 ; (3) when the Ba/Eu value is less than 2950 and greater than 1938.06, the interference prediction model is a linear regression, where y = 5.88574265 + 0.00265247237 × x , 1938.06 < x < 2950.0 ; and (4) when the Ba/Eu value is less than 1938.06, the interference prediction model is a linear regression, where y = 0.0169820389 + 0.0000320448913 × x , x < 1938.06 .
For (tolS, tolN) = (0,4), (0,5), (0,6), (1,4), (1,5), (1,6), (2,4), (2,5), (2,6), the model trees for interference prediction of Ba on Eu are shown in Figure 2c. (1) When the Ba/Eu value is less than 2950, the Ba interference with Eu prediction model is a linear regression, with y = 0.18419411 + 0.00050737 × x , 0 < x < 2950 ; (2) when the Ba/Eu value is greater than 2950 and less than 189,523, the interference prediction model is a linear regression, where y = 0.293982186 + 0.00000181729975 × x , 2950 < x < 189 , 523 .
For (tolS, tolN) = (0,7,…,∞), (1,7,…,∞), (2,7,…,∞), (3,…,∞,1, …,∞), the model trees for prediction of Ba interference with Eu are shown in Figure 2d.
For (tolS, tolN) = (1,1), (1,2), (2,1), and (2,2), the model trees for prediction of Ba interference with Eu are shown in Figure 2e. (1) When the Ba/Eu value is greater than 2950, the interference prediction model is a linear regression, with y = 0.293982186 + 0.00000181729975 × x , x > 2950 ; (2) when the Ba/Eu value is less than 2950 and greater than 2247.06, the interference prediction model is a linear regression, where y = 1.29 + 0.0002 × x , 2247.06 < x < 2950 ; and (3) when the Ba/Eu value is less than 2247.06, the interference prediction model is a linear regression, with y = 0.017931963 + 0.0000267651624 × x , x < 2247.06 .
For (tolS, tolN) = (1,3), (2,3), the model trees for prediction of Ba interference with Eu are shown in Figure 2f. (1) When the Ba/Eu value is greater than 2950, the interference prediction model is a linear regression, with y = 0.293982186 + 0.00000181729975 × x , x > 2950 ; (2) when the Ba/Eu value is less than 2950 and greater than 1938.06, the interference prediction model is a linear regression, where y = 5.88574265 + 0.00265247237 × x , 1938.06 < x < 2950 ; and (3) when the Ba/Eu value is less than 1938.06, the interference prediction model is a linear regression, with y = 0.0169820389 + 0.0000320448913 × x , x < 1938.06 .

3.3. Results

For all the model trees of the Ba-Eu interference prediction model illustrated in Figure 2, the training data sets of the Ba/Eu ratio and the Ba interference with Eu were scattered, as shown in Figure 1. The points (2900,1.87) and (2950,1.88) were outliers from the other Ba/Eu ratio and Ba interference with Eu data points.
All prediction models with (tolS, tolN) = (0,1), (0,2), (0,3), (1,1), (1,2), (2,1), (2,2), (1,3), and (2,3) contained outlier points of (2950,1.88) and (2900,1.87). The prediction models with (tolS,tolN) = (0,7,…,∞), (1,7,…,∞), (2,7, …,∞), (3,…,∞,1,…,∞) yielded lower prediction accuracy; thus, these models were imprecise.
The optimal values of (tolS, tolN) for the prediction model of Ba interference with Eu were found to be (0,4), (0,5), (0,6), (1,4), (1,5), (1,6), (2,4), (2,5), and (2,6). When the Ba/Eu value was less than 2950, a linear regression was obtained for the interference prediction model, where y = 0.18419411 + 0.00050737 × x , 0 < x < 2950 . Further, when the Ba/Eu value was greater than 2950 and less than 189,523, the interference prediction model was found to be a linear regression, with y = 0.293982186 + 0.00000181729975 × x , 2950 < x < 189 , 523 . From the optimal models, a threshold point value of 363.0370538 could be determined. Note that, when the Ba/Eu value is 363.0370538, it is not necessary to consider the Ba interference with Eu; thus, the Eu values can be interpreted from the data for the investigated samples.

3.4. Performance Evaluation

To verify the threshold value for Ba interference with Eu proposed in this paper, a wide dataset of Ba/Eu values covering 2–361 through 379–938 to 1042–3305 from previously published literature was used (Table 5 and Table 6) [62,64,65]. The data for the testing were selected from Dai et al. [62,64] and Duan et al. [65], because these data points were all obtained via ICP-Q-MS. Thus, the Ba concentrations were expected to interfere with the Eu concentrations in the samples if the Ba/Eu values exceeded the threshold value, either at 1000 (as proposed in previous works) or at 363 (as proposed in this study). A total of 41 coal bench samples from a boehmite-rich 36.37-m-thick Pennsylvanian coal seam in Inner Mongolia, northern China, were considered, which were reported by Dai et al. [64]. A total of 60 coal bench samples from three Ge-rich Neogene coals from Lincang, Yunnan Province, southwestern China, were considered, which were reported by Dai et al. [62]. Further, a total of 27 coal bench samples from Reshuihe, Zhenxiong, Yunnan Province, China, were considered, which were reported by Duan et al. [65]. The test datasets presented in Table 5 and Table 6 could be classified into three groups: Ba/Eu < 363, Ba/Eu = 363–1000, and Ba/Eu > 1000. We compared the threshold value of 363 determined from our proposed model with the value of 1000 proposed by others (e.g., [4,27,35]). Based on Table 5, the correlation coefficient of Ba and Eu was 0.1326 and 0.9659 when Ba/Eu was < 363 and >1000, respectively. When Ba/Eu varied from 363 to 1000, the Ba and Eu correlation coefficient remained as high as 0.9545, as illustrated in Figure 3A. Based on Table 6, the correlation coefficient of Ba and Eu was 0.231 and 0.9318 when Ba/Eu was < 363 and >1000, respectively. When Ba/Eu varied from 363 to 1000, the Ba and Eu correlation coefficient remained as high as 0.9317, as illustrated in Figure 3B. The distinctively different correlation coefficients for Ba and Eu in the different Ba/Eu ranges show that the threshold value of 363 is more accurate for determination of Ba interference with Eu than the previously proposed value of 1000.

4. Conclusions

In conclusion, to determine the threshold value of Ba interference with Eu in the context of ICP-Q-MS data analysis, three machine learning techniques—namely, the linear regression, regression tree, and model tree methods—were used to construct prediction models of Ba interference with Eu in coal and coal-related samples. The CART algorithm was applied to the tree regression. To apply the models for prediction of Ba interference with Eu, all related data, including that on Ba, Eu, Ba/Eu, and Eu interference, were collected and prepared. A Ba-Eu interference linear regression model, regression tree, and model tree were implemented in Python for prediction. The results showed that the model tree is far superior to the regression tree for determination of Ba/Eu threshold points. The extracted feature was Ba/Eu and the extracted feature value was the interference of Ba with Eu. From all obtained prediction models, an optimal threshold point value of 363 was determined. This indicates that, when the Ba/Eu value is <363, the Ba interference with Eu can be neglected; thus, the Eu concentrations in samples can be determined based on ICP-Q-MS data. Based on the results of simulations in which the threshold value of 363 proposed in this study and that of 1000 proposed in other works (e.g., [4,27,35]) were compared, the former is more accurate for determining whether Ba interferes with Eu in investigated samples. In the future, we will use deep learning techniques [66,67,68] to determine the threshold value of Ba interference with Eu.

Author Contributions

Conceptualization, N.X. and Q.L.; methodology, N.X. and Q.L.; validation, N.X. and Q.L.; formal analysis, N.X. and Q.L.; Writing—Original Draft preparation, N.X.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 61772320, and 111 Projects, grant number B17042.

Acknowledgments

We thank Xudong Mao of the City University of Hong Kong for help with the manuscript revision.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barakos, G.; Mischo, H.; Gutzmer, J. An outlook on the rare earth elements mining industry. 2016. Available online: https://www.ausimmbulletin.com/feature/an-outlook-on-the-rare-earth-elements-mining-industry/ (accessed on 1 November 2018).
  2. Hower, J.C.; Granite, E.J.; Mayfield, D.B.; Lewis, A.S.; Finkelman, R.B. Notes on contributions to the science of rare earth element enrichment in coal and coal combustion by-products. Minerals 2016, 6, 32. [Google Scholar] [CrossRef]
  3. Seredin, V.V.; Dai, S. Coal deposits as potential alternative sources for lanthanides and yttrium. Int. J. Coal Geol. 2012, 94, 67–93. [Google Scholar] [CrossRef]
  4. Dai, S.; Graham, I.T.; Ward, C.R. A review of anomalous rare earth elements and yttrium in coal. Int. J. Coal Geol. 2016, 159, 82–95. [Google Scholar] [CrossRef]
  5. Dai, S.; Ward, C.R.; Graham, I.T.; French, D.; Hower, J.C.; Zhao, L.; Wang, X. Altered volcanic ashes in coal and coal-bearing sequences: A review of their nature and significance. Earth Sci. Rev. 2017, 175, 44–74. [Google Scholar] [CrossRef]
  6. Dai, S.; Yan, X.; Ward, C.R.; Hower, J.C.; Zhao, L.; Wang, X.; Zhao, L.; Ren, D.; Finkelman, R.B. Valuable elements in Chinese coals: A review. Int. Geol. Rev. 2018, 60, 590–620. [Google Scholar]
  7. Dai, S.; Finkelman, R.B. Coal as a promising source of critical elements: Progress and future prospects. Int. J. Coal Geol. 2018, 186, 155–164. [Google Scholar] [CrossRef]
  8. Seredin, V.V. REE-bearing coals from Russian Far East deposits. Int. J. Coal Geol. 1996, 30, 101–129. [Google Scholar] [CrossRef]
  9. Hower, J.C.; Ruppert, L.F.; Eble, C.F. Lanthanide, yttrium, and zirconium anomalies in the fire clay coal bed, Eastern Kentucky. Int. J. Coal Geol. 1999, 39, 141–153. [Google Scholar] [CrossRef]
  10. Hower, J.C.; Eble, C.F.; Dai, S.; Belkin, H.E. Distribution of rare earth elements in eastern Kentucky coals: Indicators of multiple modes of enrichment? Int. J. Coal Geol. 2016, 160–161, 73–81. [Google Scholar] [CrossRef]
  11. Hower, J.C.; Dai, S. Petrology and chemistry of sized Pennsylvania anthracite, with emphasis on the distribution of rare earth elements. Fuel 2016, 185, 305–315. [Google Scholar] [CrossRef]
  12. Lin, R.; Soong, Y.; Granite, E.J. Evaluation of trace elements in U.S. coals using the USGS COALQUAL database version 3.0. Part I: Rare earth elements and yttrium (REY). Int. J Coal Geol. 2018, 192, 1–13. [Google Scholar] [CrossRef]
  13. Taggart, R.K.; Hower, J.C.; Dwyer, G.S.; Hsu-Kim, H. Trends in the rare earth element content of U.S.-based coal combustion fly ashes. Environ. Sci. Technol. 2016, 50, 5919–5926. [Google Scholar] [CrossRef]
  14. Kolker, A.; Scott, C.; Hower, J.C.; Vazquez, J.A.; Lopano, C.L.; Dai, S. Distribution of rare earth elements in coal combustion fly ash, determined by SHRIMP-RG ion microprobe. Int. J. Coal Geol. 2017, 184, 1–10. [Google Scholar] [CrossRef]
  15. Laudal, D.A.; Benson, S.A.; Addleman, R.S.; Palo, D. Leaching behavior of rare earth elements in Fort Union lignite coals of North America. Int. J Coal Geol. 2018, 191, 112–124. [Google Scholar] [CrossRef]
  16. Zhang, W.; Honaker, R.Q. Rare earth elements recovery using staged precipitation from a leachate generated from coarse coal refuse. Int. J Coal Geol. 2018, 195, 189–199. [Google Scholar] [CrossRef]
  17. Vito, I.E.D.; Olsina, R.A.; Masi, A.N. Enrichment method for trace amounts of rare earth elements using chemofiltration and XRF determination. Fresen. J. Anal. Chem. 2000, 368, 392–396. [Google Scholar] [CrossRef]
  18. Smoliński, A.; Stempin, M.; Howaniec, N. Determination of rare earth elements in combustion ashes from selected Polish coal mines by wavelength dispersive X-ray fluorescence spectrometry. Spectrochim. Acta B 2016, 116, 63–74. [Google Scholar] [CrossRef]
  19. Ohde, S. Determination of rare earth elements in carbonatites from the kangankunde mine, Malawi by instrumental neutron activation analysis. J. Radioanal. Nucl. Chem. 2003, 257, 433–435. [Google Scholar] [CrossRef]
  20. Kumar, K.; Saion, E.; Halimah, M.K.; Yap, C.K.; Hamzah, M.S. Rare earth element (REE) in surface mangrove sediment by instrumental neutron activation analysis. J. Radioanal. Nucl. Chem. 2014, 301, 667–676. [Google Scholar] [CrossRef]
  21. El-Taher, A. Rare earth elements content in geological samples from the Eastern Desert, Egypt, determined by instrumental neutron activation analysis. Appl. Radiat. Isot. 2010, 68, 1859–1863. [Google Scholar] [CrossRef]
  22. Cook, N.J.; Ciobanu, C.L.; O’Rielly, D.; Wilson, R.; Das, K.; Wade, B. Mineral chemistry of rare earth element (REE) mineralization, Browns Ranges, Western Australia. Lithos 2013, 172–173, 192–213. [Google Scholar] [CrossRef]
  23. Thompson, R.L.; Bank, T.; Roth, E.; Granite, E. Resolution of rare earth element interferences in fossil energy by-product samples using sector-field ICP-MS. Fuel 2016, 185, 94–101. [Google Scholar] [CrossRef]
  24. Gaft, M.; Dvir, E.; Modiano, H.; Schone, U. Laser induced breakdown spectroscopy machine for online ash analyses in coal. Spectrochim Acta B 2008, 63, 1177–1182. [Google Scholar] [CrossRef]
  25. Stankova, A.; Gilon, N.; Dutruch, L.; Kanicky, V. A simple LIBS method for fast quantitative analysis of fly ashes. Fuel 2010, 89, 3468–3874. [Google Scholar] [CrossRef]
  26. Haider, A.F.M.Y.; Rony, M.A.; Lubna, R.S.; Abedin, K.M. Detection of multiple elements in coal samples from Bangladesh by laser-induced breakdown spectroscopy. Opt. Laser Technol. 2011, 43, 1405–1410. [Google Scholar] [CrossRef]
  27. Yan, X.; Dai, S.; Graham, I.T.; He, X.; Shan, K.; Liu, X. Determination of Eu concentrations in coal, fly ash and sedimentary rocks using a cation exchange resin and inductively coupled plasma mass spectrometry (ICP-MS). Int. J. Coal Geol. 2018, 191, 152–156. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Jiang, Z.; He, M.; Hu, B. Determination of trace rare earth elements in coal fly ash and atmospheric particulates by electrothermal vaporization inductively coupled plasma mass spectrometry with slurry sampling. Environ. Pollut. 2007, 148, 459–467. [Google Scholar] [CrossRef]
  29. Ardini, F.; Soggia, F.; Rugi, F.; Udisti, R.; Grotti, M. Comparison of inductively coupled plasma spectrometry techniques for the direct determination of rare earth elements in digests from geological samples. Anal. Chim. Acta 2010, 678, 18–25. [Google Scholar] [CrossRef]
  30. Gray, A.L.; Williams, J.G. Oxide and doubly charged ion response of a commercial inductively coupled plasma mass spectrometry instrument. J. Anal. At. Spectrom. 1987, 2, 51–82. [Google Scholar] [CrossRef]
  31. Jarvis, K.E.; Gray, A.L.; McCurdy, E. Avoidance of spectral interference on europium in inductively coupled plasma mass spectrometry by sensitive measurement of the doubly charged ion. J. Anal. At. Spectrom. 1989, 4, 743–747. [Google Scholar] [CrossRef]
  32. Raut, N.M.; Huang, L.S.; Aggarwal, S.K.; Lin, K.C. Determination of lanthanides in rock samples by inductively coupled plasma mass spectrometry using thorium as oxide and hydroxide correction standard. Spectrochim. Acta B 2003, 58, 809–822. [Google Scholar] [CrossRef]
  33. Raut, N.M.; Huang, L.S.; Lin, K.C.; Aggarwal, S.K. Uncertainty propagation through correction methodology for the determination of rare earth elements by quadrupole based inductively coupled plasma mass spectrometry. Anal. Chim. Acta 2005, 530, 91–103. [Google Scholar] [CrossRef]
  34. Zawisza, B.; Pytlakowska, K.; Feist, B.; Polowniak, M.; Kita, A.; Sitko, R. Determination of rare earth elements by spectroscopic techniques: A review. J. Anal. Atom. Spectrom. 2011, 26, 2373–2390. [Google Scholar] [CrossRef]
  35. Loges, A.; Wagner, T.; Barth, M.; Bau, M.; Göb, S.; Markl, G. Negative Ce anomalies in Mn oxides: The role of Ce4+ mobility during water-mineral interaction. Geochim. Cosmochim. Acta. 2012, 86, 296–317. [Google Scholar] [CrossRef]
  36. Ketris, M.P.; Yudovich, Y.E. Estimations of Clarkes for Carbonaceous biolithes: World average for trace element contents in black shales and coals. Int. J. Coal Geol. 2009, 78, 135–148. [Google Scholar] [CrossRef]
  37. Palmer, C.A.; Oman, C.L.; Park, A.J.; Luppens, J.A. The U.S. Geological Survey Coal Quality (COALQUAL) Database Version 3.0. Data Series; The U.S. Geological Survey: Reston, VA, USA, 2015; p. 57.
  38. Finkelman, R.B.; Palmer, C.A.; Wang, P. Quantification of the modes of occurrence of 42 elements in coal. Int. J. Coal Geol. 2018, 185, 138–160. [Google Scholar] [CrossRef]
  39. Dai, S.F.; Ren, D.Y.; Chou, C.-L.; Li, S.S.; Jiang, Y.F. Mineralogy and geochemistry of the No. 6 Coal (Pennsylvanian) in the Junger Coalfield, Ordos Basin, China. Int. J. Coal Geol. 2006, 66, 253–270. [Google Scholar] [CrossRef]
  40. Dai, S.F.; Wang, X.B.; Zhou, Y.P.; Hower, J.C.; Li, D.H.; Chen, W.M.; Zhu, X.W.; Zou, J.H. Chemical and mineralogical compositions of silicic, mafic, and alkali tonsteins in the Late Permian coals from the Songzao Coalfield, Chongqing, Southwest China. Chem. Geol. 2011, 282, 29–44. [Google Scholar] [CrossRef]
  41. Liu, J.J.; Yang, Z.; Yan, X.Y.; Ji, D.; Yang, Y.P.; Hu, L.C. Modes of occurrence of highly-elevated trace elements in superhigh-organic-sulfur coals. Fuel 2015, 156, 190–197. [Google Scholar] [CrossRef]
  42. Johnston, M.N.; Hower, J.C.; Dai, S.; Wang, P.; Xie, P.; Liu, J. Petrology and Geochemistry of the Harlan, Kellioka, and Darby Coals from the Louellen 7.5-Minute Quadrangle, Harlan County, Kentucky. Minerals 2015, 5, 894–918. [Google Scholar] [CrossRef]
  43. Dai, S.; Ren, D.; Zhou, Y.; Chou, C.-L.; Wang, X.; Zhao, L.; Zhu, X. Mineralogy and geochemistry of a superhigh-organic-sulfur coal, Yanshan Coalfield, Yunnan, China: Evidence for a volcanic ash component and influence by submarine exhalation. Chem. Geol. 2008, 255, 182–194. [Google Scholar] [CrossRef]
  44. Dai, S.; Xie, P.; Jia, S.; Ward, C.R.; Hower, H.C.; Yan, X.; French, D. Enrichment of U-Re-V-Cr-Se and rare earth elements in the Late Permian coals of the Moxinpo Coalfield, Chongqing, China: Genetic implications from geochemical and mineralogical data. Ore Geol. Rev. 2017, 80, 1–17. [Google Scholar] [CrossRef]
  45. Hower, J.C.; Eble, C.F.; O’Keefe, J.M.K.; Dai, S.; Wang, P.; Xie, P.; Liu, J.; Ward, C.R.; French, D. Petrology, Palynology, and Geochemistry of Gray Hawk Coal (Early Pennsylvanian, Langsettian) in Eastern Kentucky, USA. Minerals 2015, 5, 592–622. [Google Scholar] [CrossRef]
  46. Liu, J.; Ward, C.R.; Graham, I.T.; French, D.; Dai, S.; Song, X. Modes of occurrence of non-mineral inorganic elements in lignites from the Mile Basin, Yunnan Province, China. Fuel 2018, 222, 146–155. [Google Scholar] [CrossRef]
  47. Chang, C.; Fu, Q.; Wang, X. Linear Correlation of Ba and Eu Contents by Hydrothermal Activities: A Case Study in the Hetang Formation, South China. Geofluids 2019. [Google Scholar] [CrossRef]
  48. Dulski, P. Interferences of oxide, hydroxide and chloride analyte species in the determination of rare earth elements in geological samples by inductively coupled plasma-mass spectrometry. Fresenius J. Anal. Chem. 1994, 350, 194–203. [Google Scholar] [CrossRef]
  49. Smirnova, E.V.; Fedorova, I.N.; Sandimirova, G.P.; Petrov, L.L.; Balbekina, N.G.; Lozhkin, V.I. Determination of rare earth elements in black shales by inductively coupled plasma mass spectrometry. Spectrochim. Acta B 2003, 58, 329–340. [Google Scholar] [CrossRef]
  50. Wang, J.; Chen, D.; Wang, D.; Yan, D.; Zhou, X.; Wang, Q. Petrology and geochemistry of chert on the marginal zone of Yangtze platform, western Hunan, South China, during the Ediacaran-Cambrian transition. Sedimentology 2012, 59, 809–829. [Google Scholar] [CrossRef]
  51. Jiang, S.Y.; Zhao, H.X.; Chen, Y.Q.; Yang, T.; Yang, J.H.; Ling, H.F. Trace and rare earth element geochemistry of phosphate nodules from the lower Cambrian black shale sequence in the Mufu Mountain of Nanjing, Jiangsu province, China. Chem. Geol. 2007, 244, 584–604. [Google Scholar] [CrossRef]
  52. Shields, G.; Stille, P. Diagenetic constraints on the use of cerium anomalies as palaeoseawater redox proxies: An isotopic and REE study of Cambrian phosphorites. Chem. Geol. 2001, 175, 29–48. [Google Scholar] [CrossRef]
  53. Dai, S.; Guo, W.; Nechaev, V.P.; French, D.; Ward, C.R.; Spiro, B.F.; Finkelman, R.B. Modes of occurrence and origin of mineral matter in the Palaeogene coal (No. 19-2) from the Hunchun Coalfield, Jilin Province, China. Int. J. Coal Geol. 2018, 189, 94–110. [Google Scholar] [CrossRef]
  54. Dai, S.; Ji, D.; Ward, C.R.; French, D.; Hower, J.C.; Yan, X.; Wei, Q. Mississippian anthracites in Guangxi Province, southern China: Petrological, mineralogical, and rare earth element evidence for high-temperature solutions. Int. J. Coal Geol. 2018, 197, 84–114. [Google Scholar] [CrossRef]
  55. Dai, S.; Xie, P.; Ward, C.R.; Yan, X.; Guo, W.; French, D.; Graham, I.T. Anomalies of rare metals in Lopingian super-high-organic-sulfur coals from the Yishan Coalfield, Guangxi, China. Ore Geol. Rev. 2017, 88, 235–250. [Google Scholar] [CrossRef]
  56. Wang, Z.; Dai, S.; Zou, J.; French, D.; Graham, I.T. Rare earth elements and yttrium in coal ash from the Luzhou power plant in Sichuan, Southwest China: Concentration, characterization and optimized extraction. Int. J. Coal Geol. 2019, 203, 1–14. [Google Scholar] [CrossRef]
  57. Harrington, P. Machine Learning in Action; Manning Publications: Shelter Island, NY, USA, 2012; p. 384. [Google Scholar]
  58. Han, J.; Pei, J.; Kamber, M. Data Mining: Concepts and Techniques; Morgan Kaufmann: Burlington, MA, USA, 2011; p. 744. [Google Scholar]
  59. Wu, X.; Kumar, V.; Ross Quinlan, J.; Ghosh, J.; Yang, Q.; Motoda, H.; McLachlan, G.J.; Ng, A.; Liu, B.; Yu, P.S.; et al. Top 10 algorithms in data mining. Knowl. Inf. Syst. 2008, 14, 1–37. [Google Scholar] [CrossRef]
  60. McKinney, W. Python for Data Analysis; O’Reilly Media: Sebastopol, CA, USA, 2017; p. 511. [Google Scholar]
  61. Dai, S.; Wang, X.; Seredin, V.V.; Hower, J.C.; Ward, C.R.; O’Keefe, J.M.K.; Huang, W.; Li, T.; Li, X.; Liu, H.; et al. Petrology, mineralogy, and geochemistry of the Ge-rich coal from the Wulantuga Ge ore deposit, Inner Mongolia, China: New data and genetic implications. Int. J. Coal Geol. 2012, 90–91, 72–99. [Google Scholar] [CrossRef]
  62. Dai, S.; Wang, P.; Ward, C.R.; Tang, Y.; Song, X.; Jiang, J.; Hower, J.C.; Li, T.; Seredin, V.V.; Wagner, N.J.; et al. Elemental and mineralogical anomalies in the coal-hosted Ge ore deposit of Lincang, Yunnan, southwestern China: Key role of N2–CO2-mixed hydrothermal solutions. Int. J. Coal Geol. 2015, 152, 19–46. [Google Scholar] [CrossRef]
  63. Wei, Q.; Rimmer, S.M.; Dai, S. Distribution of trace elements in fractions after micronization and density-gradient centrifugation of high-Ge coals from the Wulantuga and Lincang Ge ore deposits, China. Energy Fuels 2017, 31, 11818–11837. [Google Scholar] [CrossRef]
  64. Dai, S.; Li, D.; Chou, C.L.; Zhao, L.; Zhang, Y.; Ren, D.; Ma, Y.; Sun, Y. Mineralogy and geochemistry of boehmite-rich coals: New insights from the Haerwusu Surface Mine, Jungar Coalfield, Inner Mongolia, China. Int. J. Coal Geol. 2008, 74, 185–202. [Google Scholar] [CrossRef]
  65. Duan, P.; Wang, W.; Liu, X.; Sang, S.; Ma, M.; Zhang, W. Differentiation of rare earth elements and yttrium in different size and density fractions of the Reshuihe coal, Yunnan Province, China. Int. J. Coal Geol. 2019, 207, 1–11. [Google Scholar] [CrossRef]
  66. Li, L.; Ota, K.; Dong, M.; Borjigin, W. Eyes in the dark: Distributed scene understanding for disaster management. IEEE Trans. Parallel Distrib. Syst. 2017, 28, 3458–3471. [Google Scholar] [CrossRef]
  67. Li, L.; Ota, K.; Dong, M. DeepNFV: A Light-weight Framework for Intelligent Edge Network Functions Virtualization. IEEE Network 2019, 33, 136–141. [Google Scholar] [CrossRef]
  68. Ota, K.; Dao, M.S.; Mezaris, V.; De Natale, F.G. Deep Learning for Mobile Multimedia: A Survey. ACM Trans. Multimedia Comput. Commun. Appl. (TOMM) 2017, 13, 21. [Google Scholar] [CrossRef]
Figure 1. Training data (from Yan et al. [27]).
Figure 1. Training data (from Yan et al. [27]).
Minerals 09 00259 g001
Figure 2. Prediction models with variables (tolS, tolN) of (a) (0,1), (0,2); (b) (0,3); (c) (0,4), (0,5), (0,6), (1,4), (1,5), (1,6), (2,4), (2,5), (2,6); (d) (0,7,…,∞), (1,7,…,∞), (2,7,…,∞), (3,…,∞, 1,…,∞); (e) (1,1), (1,2), (2,1), (2,2); (f) (1,3) (2,3).
Figure 2. Prediction models with variables (tolS, tolN) of (a) (0,1), (0,2); (b) (0,3); (c) (0,4), (0,5), (0,6), (1,4), (1,5), (1,6), (2,4), (2,5), (2,6); (d) (0,7,…,∞), (1,7,…,∞), (2,7,…,∞), (3,…,∞, 1,…,∞); (e) (1,1), (1,2), (2,1), (2,2); (f) (1,3) (2,3).
Minerals 09 00259 g002
Figure 3. Performance comparison with test data of (A) Table 5; (B) Table 6.
Figure 3. Performance comparison with test data of (A) Table 5; (B) Table 6.
Minerals 09 00259 g003
Table 1. Description for Ba/Eu in related samples used by Yan et al. [27].
Table 1. Description for Ba/Eu in related samples used by Yan et al. [27].
SampleBa/EuTypeDescription
SRM2682b2247CoalNational Institute of Standards and Technology (NIST) standard reference samples
SRM2685b292Bituminous Coal
SRM26902900Fly ash
SRM26912950
WLTG C6-218,598Low-rank CoalNo. 6 coal of Wulantuga Deposit (Shengli Coalfield, Inner Mongolia [61])
ZJ-4-63813Low-rank CoalNo. 4 coal of Zhoujing Mine, Baise Coalfield, Guangxi Province
ZJ-5-122083No. 5 coal from Zhoujing Mine, Baise Coalfield, Guangxi Province
X1-1R202,200Carbonate metasomatitesDazhai Mine, Lincang Ge ore deposit, Yunnan Province [62,63]
X1-2R42,236
Z2-15F51,027Quartz-carbonate metasomatites
Z2-16F33,816
LL5-K3-813.18Semi-anthraciteNo. K3 coal from the La-Lang 5 Mine, Yishan Coalfield, Guangxi Province [55]
LL5-K3-1310.69
Table 2. Certified (Cer), Before separation (BS) and after separation (Steps 1–4) values of Eu and Ba (μg/g) in coal and coal-related samples [27] (BDL: below detection limit).
Table 2. Certified (Cer), Before separation (BS) and after separation (Steps 1–4) values of Eu and Ba (μg/g) in coal and coal-related samples [27] (BDL: below detection limit).
ElementsSRM2690 (Ba/Eu = 2900)X1-1R (Ba/Eu = 202,200)X1-2R (Ba/Eu = 42,236)
CerBSStep 1Step 2Step 3Step 4BSStep 1Step 2Step 3Step 4BSStep 1Step 2Step 3Step 4
153Eu2.004.01BDL0.001.872.000.52BDL0.000.620.010.25BDLBDL0.250.02
137Ba5800.0063901.060.345852.23112.101895.230.450.072022.0022.03814.534.020.09844.6712.66
ElementsSRM2691 (Ba/Eu = 2950)Z2-15F (Ba/Eu = 51,027)Z2-16F (Ba/Eu = 33,816)
CerBSStep 1Step 2Step 3Step 4BSStep 1Step 2Step 3Step 4BSStep 1Step 2Step 3Step 4
153Eu2.004.08BDLBDL1.881.930.44BDLBDL0.450.030.38BDLBDL0.420.05
137Ba5900.006109.000.880.526392.00193.001493.230.790.031530.8932.791357.232.850.261690.8922.93
ElementsSRM2682b (Ba/Eu = 2247)ZJ-4-6 (Ba/Eu = 3813)ZJ-5-12 (Ba/Eu = 2083)
CerBSStep 1Step 2Step 3Step 4BSStep 1Step 2Step 3Step 4BSStep 1Step 2Step 3Step 4
153Eu0.170.23BDLBDL0.070.160.18BDL0.000.080.100.24BDLBDL0.080.17
137Ba382.00368.77BDL0.54407.210.83350.47BDL0.42381.32BDL329.471.800.60353.99BDL
ElementsSRM2685b (Ba/Eu = 292)WTGC6-2 (Ba/Eu = 18,598)LL5-K3-8 (Ba/Eu = 13.18)
CerBSStep 1Step 2Step 3Step 4BSStep 1Step 2Step 3Step 4BSStep 1Step 2Step 3Step 4
153Eu0.360.33BDLBDL0.020.340.61BDLBDL0.540.142.59BDLBDL0.032.41
137Ba105.0097.60BDL0.40113.10BDL2428.270.400.402603.77BDL26.233.430.6831.773.72
ElementsLL5-K3-13(Ba/Eu = 10.69)
BSStep 1Step 2Step 3Step 4
153Eu2.23BDL0.01BDL2.24
137Ba19.00BDLBDL23.942.15
Table 3. Training data [27].
Table 3. Training data [27].
SampleGroup NoBa/EuBa Interference with Eu
SRM2682b18.480.01
SRM2685b210.880.03
SRM26903291.670.02
SRM269141938.060.08
WLTGC6-252247.060.07
ZJ-4-6629001.87
ZJ-5-12729501.88
X1-1R83504.70.08
X1-2R917,344.790.54
Z2-15F1027,144.60.42
Z2-16F1140,726.50.25
LL5-K3-81249,774.330.45
LL5-K3-1313189,5230.62
Table 4. Prediction models for Ba interference with Eu based on model tree.
Table 4. Prediction models for Ba interference with Eu based on model tree.
Variables (tolS, tolN)Prediction Models for Ba Interference with Eu Based on Model Tree
(0,1), (0,2) y = 0.389451044 + 0.00000121646954 × x , x > 40 , 726.5 ;
y = 0.759759680 0.0000125166582 × x , 17 , 344.79 < x < 40 , 726.5 ;
y = 0.0364849362 + 0.0000332367781 × x , 2950 < x < 17 , 344.79 ;
y = 1.29 + 0.0002 × x , 2247.06 < x < 2950 ;
y = 0.0131201121 + 0.0000291815483 × x , 10.88 < x < 2247.06 ;
y = 0.06066667 + 0.00833333 × x , x < 10.88
(0,3) y = 0.259761293 + 0.00000193096467 × x , x > 27 , 144.6
y = 0.0966488103 + 0.0000156280402 × x , 2950.0 < x < 27 , 144.6
y = 5.88574265 + 0.00265247237 × x , 1938.06 < x < 2950.0
y = 0.0169820389 + 0.0000320448913 × x , x < 1938.06
(0,4), (0,5), (0,6), (1,4), (1,5), (1,6), (2,4), (2,5), (2,6) y = 0.18419411 + 0.00050737 × x , 0 < x < 2950
y = 0.293982186 + 0.00000181729975 × x , 2950 < x < 189 , 523
(0,7,…,∞), (1,7,…,∞), (2,7,…,∞), (3,…,∞, 1,…,∞) y = 0.471722528 + 5.54453477 × 10 7 × x
(1,1), (1,2), (2,1), (2,2) y = 0.293982186 + 0.00000181729975 × x , x > 2950
y = 1.29 + 0.0002 × x , 2247.06 < x < 2950
y = 0.017931963 + 0.0000267651624 × x , x < 2247.06
(1,3) (2,3) y = 0.293982186 + 0.00000181729975 × x , x > 2950
y =   5.88574265 + 0.00265247237 × x , 1938.06 < x < 2950
y = 0.0169820389 + 0.0000320448913 × x , x < 1938.06
Table 5. Ba, Eu, and Ba/Eu data for test (Data from Dai et al. [62,64]).
Table 5. Ba, Eu, and Ba/Eu data for test (Data from Dai et al. [62,64]).
SampleBaEuBa/Eu
S3-1R4340.61711.4754098
S3-2R3151.63193.2515337
S3-481.30.12677.5
S3-587.30.14623.5714286
S3-685.90.16536.875
S3-793.00.10930
S3-81130.091255.555556
S3-9F2130.59361.0169492
S3-10F5360.79678.4810127
S3-11F5610.79710.1265823
WA-S390.60.13696.9230769
Z2-1R4660.60776.6666667
Z2-21030.14735.7142857
Z2-31110.101110
Z2-4P4480.172635.294118
Z2-5P94.60.071351.428571
Z2-5LP1560.121300
Z2-6P2850.51558.8235294
Z2-71010.091122.222222
Z2-897.00.091077.777778
Z2-979.40.11721.8181818
Z2-101710.31551.6129032
Z2-11P2130.26819.2307692
Z2-121260.091400
Z2-131070.13823.0769231
Z2-142260.211076.190476
Z2-15F13980.433251.162791
Z2-16F13050.413182.926829
WA-Z21220.13938.4615385
X1-1R18180.553305.454545
X1-2R7390.252956
X1-3R8430.342479.411765
X1-484.20.13647.6923077
X1-577.70.12647.5
X1-686.10.12717.5
X1-767.80.12565
X1-899.30.11902.7272727
X1-91460.141042.857143
X1-1087.00.16543.75
X1-1182.50.14589.2857143
X1-1275.70.22344.0909091
X1-1386.90.13668.4615385
X1-1470.60.10706
X1-151430.28510.7142857
X1-16F2761.39198.5611511
X1-17F2471.16212.9310345
X1-18F5820.80727.5
WA-X1910.14650
1418-1580.18322.2222222
1418-2540.13415.3846154
1418-31480.39379.4871795
H-15250.550
H-16280.835
H-17240.4652.17391304
H-18160.2176.19047619
H-19230.4254.76190476
H-20260.22118.1818182
H-21180.3946.15384615
H-22210.3953.84615385
H-22-23-P250.14178.5714286
H-23180.3847.36842105
H-24200.3066.66666667
H-24-25-P240.17141.1764706
H-25340.937.77777778
H-26160.532
H-27240.4948.97959184
H-28300.933.33333333
H-29280.931.11111111
H-B1170.3253.125
H-B2281.125.45454545
H-B31181.769.41176471
WG-11760.21838.0952381
CS-1680.16425
1104/150.60.1506
H-T10291.01029
H-1250.550
H-1-2-P453.214.0625
H-4741.452.85714286
H-5320.935.55555556
H-5-6-P1571.151.81818182
H-5-6-P2250.462.5
H-6350.938.88888889
H-7461.141.81818182
H-8300.650
H-8-9-P230.09255.5555556
H-9151.69.375
H-10280.740
H-11220.4450
H-12260.4261.9047619
H-13210.4151.2195122
H-14220.544
S3-481.30.12677.5
1418-4730.08912.5
Table 6. Ba, Eu, and Ba/Eu data for test (Data from Duan et al. [65]).
Table 6. Ba, Eu, and Ba/Eu data for test (Data from Duan et al. [65]).
SampleBaEuBa/Eu
141.370.4396.2093023
237.070.36102.972222
357.960.41141.365854
4102.690.44233.386364
5142.940.44324.863636
6223.670.33678
734.810.3891.6052632
830.590.3880.5
933.970.3791.8108108
1052.040.42123.904762
11110.630.41269.829268
12158.850.44361.022727
13240.50.4601
1432.90.3984.3589744
1545.710.26175.807692
1667.950.32212.34375
17132.770.38349.394737
18231.420.35661
19364.430.42868
20546.40.62881
2125.160.2793.1851852
2228.980.3974.3076923
23749.070.461628.41304
241304.880.582249.7931
251648.570.722289.68056
262296.270.942442.84043
273086.021.591940.89308
2828.980.3974.30769231
1–7, Size (mm): 6–13, Density (kg/L): <1.4, 1.4–1.5, 1.5–1.6, 1.6–1.7, 1.7–1.8, >1.8, Feed coal; 8–14, Size (mm): 3–6 mm, Density (kg/L): <1.4, 1.4–1.5, 1.5–1.6, 1.6–1.7, 1.7–1.8, >1.8, Feed coal; 15–21, Size (mm): 0.5–3 mm; Density (kg/L): <1.4, 1.4–1.5, 1.5–1.6, 1.6–1.7, 1.7–1.8, >1.8, Feed coal; 22–28, Size (mm): <0.5 mm, Density (kg/L): <1.4, 1.4–1.5, 1.5–1.6, 1.6–1.8, >1.8, Feed coal.

Share and Cite

MDPI and ACS Style

Xu, N.; Li, Q. Threshold Value Determination Using Machine Learning Algorithms for Ba Interference with Eu in Coal and Coal Combustion Products by ICP-MS. Minerals 2019, 9, 259. https://doi.org/10.3390/min9050259

AMA Style

Xu N, Li Q. Threshold Value Determination Using Machine Learning Algorithms for Ba Interference with Eu in Coal and Coal Combustion Products by ICP-MS. Minerals. 2019; 9(5):259. https://doi.org/10.3390/min9050259

Chicago/Turabian Style

Xu, Na, and Qing Li. 2019. "Threshold Value Determination Using Machine Learning Algorithms for Ba Interference with Eu in Coal and Coal Combustion Products by ICP-MS" Minerals 9, no. 5: 259. https://doi.org/10.3390/min9050259

APA Style

Xu, N., & Li, Q. (2019). Threshold Value Determination Using Machine Learning Algorithms for Ba Interference with Eu in Coal and Coal Combustion Products by ICP-MS. Minerals, 9(5), 259. https://doi.org/10.3390/min9050259

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop