Characterization and Identification of Coal and Carbonaceous Shale Using Visible and Near-Infrared Reflectance Spectroscopy

Because of the high organic carbon concentration in carbonaceous shale, a large proportion of carbonaceous shales are often misclassified into coals using visible and near-infrared (VIS-NIR) reflectance spectroscopy in the field of coal-gangue identification of hyperspectral remote sensing of coal mine. In order to study spectral characterization of coal and carbonaceous shale, three bituminite samples and three carbonaceous shales were collected from a coal mine of China, and their spectral reflectance curves were obtained by a field spectrometer in the wavelength range of 350–2500 nm. Only one carbonaceous shale could be easily identified from the three bituminite samples according to obvious absorption valleys near 1400 nm, 1900 nm, and 2200 nm of its reflectance curve while the other two carbonaceous shales have similar reflectance curves to the three bituminite samples. +e effect of carbon concentration on reflectance curve was simulated by the mixed powder of ultralow ash bituminite and clay in 0.5mm grain size under various mixing ratios. It was found that absorption valleys near 1400 nm, 1900 nm, and 2200 nm of the mixed powder become not obvious when the bituminite content is more than 30%. In order to establish an effective identification method of coal and carbonaceous shale, 250 other samples collected from the same coal mine were divided into 150 training samples and 100 prediction samples. Principal component analysis (PCA) and Gauss radial basis kernel principal component analysis (GRB-KPCA) were employed to extract principal components (PCs) of continuum removed (CR) spectra of the training samples in eight selected wavelength regions which are related to the main mineral and organic compositions. Two support vector machine(SVM-) based models PCA-SVM and GRB-KPCA-SVMwere established. +e results showed that the GRB-KPCA-SVM model had better identification accuracies of 94% and 92% for powder and nature block prediction samples, respectively.


Introduction
Coal is the main energy source in China and still the irreplaceable main energy in the foreseeable future.Carbonaceous shale which is sandwiched in coal seam and was formed during coal formation is mainly composed of minerals and organic carbon.During the coal production, lots of carbonaceous shales which are the main part of gangues are often cut together with coals and transported to the ground.ese gangues which occupy a large amount of land on the ground can cause damages to the environment, and the ones which are mixed in coals increase the ineffective transports [1,2].
Due to the huge coal demand and coal production capacity requirement, lots of advanced spectral technologies have been employed in coal mine investigation, coal-gangue identification, and coal property determination [3][4][5][6].Reflectance spectrum is easily acquired using cheap instrument and suitable for online analysis in the visible and nearinfrared (VIS-NIR) region [7].However, in coal mine areas, most carbonaceous shales have similar black color to coals.Meanwhile, due to the high carbon concentration in carbonaceous shale, there is no obviously different spectral absorption feature between coals and many carbonaceous shales, and they may have similar spectral characterization in the VIS-NIR region.erefore, a large proportion of carbonaceous shales are often misclassified into coals using VIS-NIR reflectance spectroscopy in the field of coal-gangue identification of hyperspectral remote sensing of coal mine [8][9][10].
Under the same measurement conditions, the difference between reflectance spectra of coal and carbonaceous shale mainly depends on the material compositions.Coal which is mainly composed of complicated organic matters including aliphatic structures, aromatic structures and oxygencontaining functional group structures often contains many minerals while carbonaceous shale which is mainly composed of clay minerals and quartz often contains many organic matters [11].
e spectral properties of coal and carbonaceous shale had been studied by some scholars in the VIS-NIR region [12,13].However, the effect of carbon concentration on the VIS-NIR reflectance spectra of coal and carbonaceous shale had seldom been investigated.Some scholars have focused on various experiments, algorithms, and models to improve the spectral identification accuracies of coals and gangues including a large number of carbonaceous shales in coal mine areas, such as Mao's normalized difference coal index (NDCI) model [8], Le's improved distribution maps model, and Song's classification method based on combination of reflectance spectra in the VIS-NIR region and spectral emissivity in the thermal infrared (TIR) region [9,10].However, all the abovementioned scholars have not accurately discriminated coals and carbonaceous shales by the reflectance spectra in the VIS-NIR region, and carbonaceous shale is the main type of gangue that reduces the classification accuracies of coals and gangues.For instance, in Song's experiment, based on the reflectance spectra in the VIS-NIR region, all carbonaceous shales were misclassified into coals in all samples including 12 coals and 9 carbonaceous shales [10].erefore, an effective identification method based on the reflectance spectra of coal and carbonaceous shale in the VIS-NIR region is very necessary to be studied.
In this paper, spectral reflectance characterization (350-2500 nm) of typical bituminite and carbonaceous shale samples of China was studied.
e effect of carbon concentration on spectral reflectance curves of coal and carbonaceous shale in the VIS-NIR region was investigated through an experiment.An effective coal-carbonaceous shale identification method based on VIS-NIR reflectance spectroscopy was established.

Spectra and Compositions Acquisition of Typical Samples.
Figure 1 shows six typical samples including three bituminite samples and three carbonaceous shales collected from Malan coal mine area in Shanxi province of China.It is obvious that the appearances of carbonaceous shales and bituminite samples in Figure 1 are similar.
Reflection spectra of rough or mat solid surfaces of remote samples are closer to reflection spectra of powder samples in larger grain size (0.25-1.6 mm) than smaller grain size to determine the compositions of the surfaces [14].For powder samples of sedimentary rock, larger particle size (0.25-1.2 mm) displays lower reflectance, but shows more obvious spectral features in their spectral bands [12].
e ancillary files of sedimentary rocks in ASTER spectral library version 2.0 of JPL pointed out that reflectance spectra of fresh (rough) solid surfaces of sedimentary or metamorphic rocks could be approximately simulated by spectra of the very coarse particulate (0.5-1.5 mm) samples [15].Some scholars had used reflectance spectra of powdered carbonate rock (a type of sedimentary rock) samples with grain size fractions between 0.125 mm and 0.5 mm and known chemical compositions to simulate reflectance spectra of fresh solid surfaces of block carbonate rock samples [16].Coal and carbonaceous shale belong to sedimentary rock.erefore, in this study, the six typical samples were crushed, air-dried, and sieved to achieve size fraction of 0.5 mm.ese powder samples are divided into three fractions which were used for spectra measurement, X-ray fluorescence (XRF) analysis, and proximate analysis, respectively.
e sample powder in 0.5 mm grain size used for spectra measurement was placed in a culture dish, and the powder surface was smoothed with the culture dish cover.en the powder surface was measured by an ASD FieldSpec field spectrometer (Analytical Spectral Devices Inc., USA) in a dark laboratory.e spectral range is 350-2500 nm, and the spectral resolution is 3 nm for 350-1000 nm, 8.5 nm for 1000-1800 nm, and 6.5 nm for 1800-2500 nm.As is shown in Figure 2, each sample was horizontal, and the fiber optic cable with field of view (FOV) of 25 °of the spectrometer was 15 cm vertical above the surface of the sample.With the distance of 15 cm to the sample, the recorded spectrum was an average spectral reflectance of the surface area of 35 cm 2 , which equals a circle of 6.5 cm in diameter.A halogen lamp at a distance of 15 cm above the surface of the sample was the light source, and the incidence angle was 45 °.Ten spectral reflectance curves of each sample were recorded in the first direction, then the sample was rotated 90 °in the horizontal plane and measured to record the second 10 curves and so on, obtaining 40 curves of the sample in all four directions.Mean of the 40 curves of each sample was calculated, and a total of 6 spectral curves of the samples in Figure 1 were obtained with a spectral dimension of 2151 in the wavelength range of 350-2500 nm.
rough XRF by an S8 Tiger spectrometer (Bruker Inc., Germany) and proximate analysis by a VMF10/6 volatile muffle furnace (Carbolite Inc., United Kingdom), an AAF12/18 ash muffle furnace (Carbolite Inc., United Kingdom), and a BS124S electronic analytical balance (Sartorius Inc., Germany), the contents of major mineral elements, volatile, and ash of these six samples were acquired.Categories (coal or rock) of these samples were determined by their ash yields according to ISO 11760 that ash yield of coal is less than or equal to 50% comparing with rock.

Spectra Acquisition of Bituminite and Clay Mixtures.
Ultralow ash bituminite powder (2% ash yield, 8% volatile) with a grain size of 0.5 mm was mixed evenly with clay powder (62% kaolinite, 19% quartz, 10% hydromuscovite, 4% montmorillonite, 3% limonite, and 1.2% illite) with a grain size of 0.5 mm under various mixing ratios to simulate the effect of carbon concentration on reflectance 2 Journal of Spectroscopy curve.e weight mixing ratios of bituminite to clay were 0 : 100, 5 : 95, 10 : 90, . .., 95 : 5, and 100 : 0 varying at the interval of 5%, and a total of 21 mixed powder samples were made.Reflectance spectra of these 21 samples were obtained under the same condition as the sample in Figure 2.

Spectra and Categories Acquisition of Modeling Samples.
Two hundred other bituminite and carbonaceous shale samples collected from Malan coal mine were randomly divided into 150 training samples and 50 prediction samples in order to establish an effective identification model of coal and carbonaceous shale.All these 200 modeling samples were crushed, air-dried, and sieved to achieve a size fraction of 0.5 mm.Each powder sample was divided into two fractions which were used for spectra measurement and proximate analysis, respectively.e training set included 56 bituminite samples and 94 carbonaceous shales, and the prediction set included 22 bituminite samples and 28 carbonaceous shales according to ISO 11760 and the ash yields acquired through proximate analysis of all these 200 samples.Reflectance spectra of these 200 powder samples in 0.5 mm grain size were obtained under the same condition as the sample in Figure 2. e preprocessing method of continuum removal (CR) was employed to remove background noise of the reflectance curves and isolate particular absorption features [17].

Spectra and Categories Acquisition of Block Samples.
Fifty extra nature block samples collected from the same coal mine as above 200 samples were air-dried and measured using the same field spectrometer and halogen lamp in the same dark laboratory.Reflectance spectra of these 50 block samples were used to verify the practicability of the identification model established by the training samples discussed above.As is shown in Figure 3, each block sample was horizontal, and the fiber optic cable of the spectrometer was vertical above the top solid surface of the sample.e recorded spectrum was the average spectral reflectance of the circular area of about 35 cm 2 .As the sample in Figure 2, mean of 40 reflectance curves in all four directions of each sample was acquired in the wavelength range of 350-2500 nm. e spectral reflectance curves of top solid surfaces of the 50 block samples were preprocessed by CR.
After reflectance spectra acquisition, the small sample pieces were cut from the measured circular area of the top solid surface of each sample and crushed into powder for proximate analysis.
e 50 block samples included 19 bituminite samples and 31 carbonaceous shales according to ISO 11760 and the ash yields acquired through proximate analysis of these 50 cut samples.

Kernel Principal Component Analysis.
Principal component analysis (PCA) is used for spectral feature extraction, and good results can be achieved when dealing with linear problems [18][19][20].However, this method has some limitations when dealing with many spectral analyses, which are often not simple linear problems.Kernel principal component analysis (KPCA) is based on the PCA [21][22][23].With the kernel function introduced, the data are mapped into a high-dimensional space, and then, the nonlinear features are extracted by the PCA method, thus the unsatisfactory results of PCA under nonlinear data distribution are improved.
e procedures of KPCA and PCA are similar, but in KPCA, the kernel function is used instead of the original data [22].For the linearly inseparable spectral data set X l×p (X l×p ⊂ R p ) of the samples, where l is the number of samples and p is the spectral dimension, and they could be mapped into a higher dimension space F through nonlinear transformation Φ.
Φ(X l×p ) is defined as and kernel function matrix is defined as Eigenvector V and eigenvalue λ of C satisfy the relationship CV � λ V and α � (α 1 , α 2 , . .., α l ) exists, making 4 Journal of Spectroscopy Based on Equation ( 5), the relationship Kα � lλ α can be concluded and α can be calculated.Based on Equation ( 5) and α, the nonlinear features of any inputted spectral data x in the higher dimension space F can be expressed by the following equation: e above process is derived when the following equation, is assumed, but in actual spectral data, K ′ is used instead of K by where In this study, Gauss radial basis (GRB) kernel function k(x i , x j ) is chosen and expressed by where x i and x j are the spectral data of the i-th and j-th sample and σ is the width parameter.
2.6.Support Vector Machine.Support vector machine (SVM) is a machine-learning algorithm based on statistic learning theory (SLT) established by Vapnik et al. in 1990s, in which structural risk minimization (SRM) criterion is used to reduce the upper bound of the model generalization error while minimizing the sample error in order to improve the generalization ability of the model [24].SVM can solve the problem of pattern recognition in small samples and nonlinear and high-dimensional data space [25,26].When dealing with nonlinear problems, the support vector classification (SVC) maps data in original space into highdimensional space and constructs an optimal classification plane in the high-dimensional space so that the distances of all samples from the plane are maximized.Radial basis function (RBF) is chosen as the kernel function of SVM in this study.

2.7.
e Proposed Procedure for Establishment of Coal-Carbonaceous Shale Prediction Models.Principal components (PCs) were obtained through PCA and KPCA of CR spectra of the training samples in eight selected wavelength regions which are related to the main mineral and organic compositions.With the PCs combined with the categories (coal and rock) as input parameters, the SVM model was trained.
e optimal width parameter σ in formula (10) of Gauss radial basis kernel principal component analysis (GRB-KPCA) was determined by cross validation.For each width parameter σ from 1 to 100 in step length of 1, the maximum average prediction accuracy (PA) of validation set of the training set was obtained through 5-fold cross validation in SVM.At the same time, based on the maximum average PA, the optimal penalty factor c and RBF variance parameter c of SVM for each σ were obtained.en the optimal width parameter σ was determined by the maximum value of all maximum average PAs of all width parameters σ.GRB-KPCA-SVM model was established according to the abovementioned procedure.Figure 4 shows flowcharts of procedures of the two prediction models: PCA-SVM and GRB-KPCA-SVM.
Prediction accuracies (PAs) of categories of the 50 powder and 50 nature block prediction samples were employed to assess the prediction models.Establishment of prediction models and assessment of the models were all performed in Matlab R2013b software (MathWorks Inc., USA).

Spectral Reflectance Curves and Compositions of Typical
Samples.Spectral reflectance curves of the six representative samples in Figure 1 in the wavelength range of 350-2500 nm are shown in Figure 5. Eight spectral feature bands are marked using No. 1-8 in Figure 5, and interactive processes of the related major compositions are shown in Table 1.In Figure 5, absorption valleys of 1400 nm (5) (5 is the spectral feature band number, and so on), 1900 nm (6), and 2200 nm (7) bands of the reflectance curve of carbonaceous shale (f ) are more obvious than those of the other five samples.Carbonaceous shale (f ) could be easily identified from the three bituminite samples according to the above three obvious absorption valleys, while the other two carbonaceous shales have similar reflectance curves to the three bituminite samples.e reflectance curves of carbonaceous shale (d) and bituminite (c) have similar obvious absorption valleys of 700 nm (2) and 870 nm (3) bands.e frequent absorption valleys of 400-550 nm (1) and 1000-1100 nm (4) bands appear in all reflectance curves, and the typical strong absorption bands from 2300 nm to 2500 nm (8) resulting from organic matters appear in all reflectance curves, too.ere is another obvious different feature between reflectance curve of carbonaceous shale (f ) and that of the other five samples that the overall shape of the reflectance curve of carbonaceous shale (f ) is straight while the other five are obviously concave.
Table 2 shows the contents of major mineral elements (SiO 2 , Al 2 O 3 , and Fe 2 O 3 ) through XRF analysis and the contents of volatile and ash through proximate analysis of these six samples in Figure 1.Matrix of XRF refers to C, H, O, and N in organic matters and matrix of carbonaceous shale (d) and (e) (47.80% and 41.25%) is much more than that of carbonaceous shale (f ) (20.24%). is may be the main reason why reflectance curves of carbonaceous shale (d) and (e) are similar to that of the three bituminite samples [12,31,32].
e obvious absorption valleys of 1400 nm, 1900 nm, and 2200 nm bands of the reflectance curve of carbonaceous shale (f ) depends on the low matrix content (20.24%) and clay and quartz minerals from which the high ash yield (78.94%) originates [27,28].e top two highest ferriferous mineral element contents (3.39% and 3.14% Fe 2 O 3 ) in Table 2 are the main cause of obvious absorption valleys of 700 nm and 870 nm bands of carbonaceous shale (d) and bituminite (c) [27,28].
According to the main compositions in Table 2 and the major related interactive processes in Table 1, the eight selected wavelength regions in Table 1 can reflect the main composition information of a bituminite or carbonaceous shale sample.From the waveforms of reflectance curves in Figure 5 and composition contents in Table 2, it can easily be deduced that matrix content and ash yield which are correlated to organic matters and minerals of a sample are the key factors to determine the waveform of reflectance curve of the sample.In other words, the effect of carbon concentration on the VIS-NIR reflectance spectra of coal and carbonaceous shale should be systemically investigated.

Spectral Reflectance Characterization of Bituminite and
Clay Mixtures.Figure 6 shows spectral reflectance curves of the 21 mixed powder samples of dry ultralow ash bituminite and clay in 0.5 mm grain size in the wavelength range of 350-2500 nm. e bituminite content increases from 0% to 100% at the interval of 5%.e eight spectral feature bands of Table 1 are also marked in Figure 6. e pure clay powder has the highest spectral reflectance curve.However, spectral reflectance curve of the mixed powder of 5% bituminite   6 Journal of Spectroscopy content dramatically falls comparing with the curve of pure clay powder especially in the near-infrared (NIR) region (780-2500 nm).When bituminite content is more than 5%, with increasing bituminite content, spectral reflectance curve of the mixed powder falls at a smaller interval.With increasing bituminite content, the depths of absorption valleys of 1400 nm, 1900 nm, and 2200 nm bands of the reflectance curve all decrease.When the bituminite content is more than 30%, absorption valleys of 1400 nm, 1900 nm, and 2200 nm bands of the reflectance curve become not obvious, and the overall shape of the reflectance curve changes from convex to concave.e overall shapes of the curves are more dense, similar, and horizontal when the bituminite contents are more than 40%.e pure bituminite powder has the lowest spectral reflectance curve, and the overall shape of the curve is close to the horizontal line.e experimental results of mixed powder of bituminite and clay explain the spectral reflectance curves of the six representative samples in Figure 5.It could be deduced that it is difficult to identify carbonaceous shale from bituminite according to the waveform of spectral reflectance curve of the sample in the VIS-NIR region when carbon concentration of the carbonaceous shale is more than about 30%.However, the material composition information contained  in the reflectance spectra in the eight selected wavelength regions in Table 1 might provide reference for identifying carbonaceous shale with high carbon concentration from bituminite.And an effective identification method based on the reflectance spectra of coal and carbonaceous shale in the VIS-NIR region should be studied.

Spectra of Modeling Samples and Establishment of
Identification Models.In order to establish an effective identification model of coal and carbonaceous shale, spectral reflectance curves of 200 powder samples in 0.5 mm grain size including 150 training samples and 50 prediction samples in the wavelength range of 350-2500 nm were obtained and shown in Figure 7(a).Figure 7(b) shows the 200 spectral curves preprocessed by CR.It is obvious that less reflectance curves (12 samples) with overall convex shapes could be easily identified from others according to the CR method in Figure 7(b).But it is difficult to identify most of the samples because of many overlappings, crossovers, and similarities among these spectral curves.Figure 7(b) indicates that the absorption features are amplified.CR spectra of the training set in the eight selected wavelength regions in Table 1 were used to establish identification models.PCA was first employed to eliminate irrelevant information and reduce dimensionality of preprocessed spectra in the eight regions.Table 3   e first ten PCs of the training set through PCA combining with the categories (0 stands for coal and 1 stands for rock) were taken as input variables to train the SVM model.5-fold cross validation was employed to determine the optimal penalty parameter c and radial basis function variance parameter c. e search ranges of c and c were both from 2 −10 to 2 10 , and step length of the exponent was 0.5.e process of determining the optimal parameter c and c is shown in Figure 8(b).With c � 2 1.5 and c � 2 −2 as optimal parameters, the maximum average PA of categories of the validation set is 90.67%.
KPCA which has more abilities to deal with the nonlinear data distribution might get a better result when employed to reduce dimensionality of preprocessed spectra of the training set in the eight regions.When the width parameter σ in Formula (10) of GRB-KPCA increases from 1 to 100 in step length of 1, Figure 9 shows the variation trend of the maximum average PA of categories of the validation set through GRB-KPCA-SVM.As is shown in Figure 9, the optimal width parameter σ is 63 for GRB-KPCA-SVM.Table 3 lists the cumulative contribution rates of the first ten PCs of CR spectra of the training set in the eight regions through GRB-KPCA (σ � 63).e cumulative contribution rate of the first three PCs   10(a).Figure 10(a) indicates that the 3D spatial cluster of coal and rock (carbonaceous shale) balls is better than that of PCA, and the mixed part of balls decreases.When the width parameter σ in Formula (10) of GRB-KPCA is 63, the first ten PCs of the training set through GRB-KPCA combining with the categories (0 stands for coal and 1 stands for rock) were taken as input variables to train the SVM model.5-fold cross validation was employed to determine the optimal penalty parameter c and radial basis function variance parameter c. e search ranges of c and c were also both from 2 −10 to 2 10 and step length of the exponent was also 0.5.e process of determining the optimal parameter c and c is shown in Figure 10(b).e optimal parameters are c � 2 4 and c � 2 0 based on GRB-KPCA-SVM, and the maximum average PA of categories of the validation set is 98%.
From the results of cross validation of the validation set, the dimension reduction method using KPCA is better than that of PCA for CR spectra of the training set in the eight selected wavelength regions in Table 1.Two coal-carbonaceous shale prediction models, PCA-SVM and GRB-KPCA-SVM, were established according to above processes.

Identification of Powder Prediction
Samples.CR spectra of the 50 prediction samples in 0.5 mm grain size in the eight selected wavelength regions in Table 1 were input into the above PCA-SVM (c � 2 1.5 , and c � 2 −2 ) and GRB-KPCA-SVM (σ � 63, c � 2 4 , and c � 2 0 ) models, and the results predicted by the two models are shown in Figure 11.12 indicates that spectral curves of block samples are more similar to spectral curves of powder samples in 0.5 mm grain size in Figure 7.Likewise, it is difficult to identify most of these block samples using the observed spectral curves in Figure 12.

Spectra and Identification of Block
In order to further verify the practicability of the above two trained models, categories of the 50 extra nature block samples were predicted.CR spectra of the 50 block samples in the eight selected wavelength regions in Table 1 were input into the above PCA-SVM (c � 2 1.5 , and c � 2 −2 ) and GRB-KPCA-SVM (σ � 63, c � 2 4 , and c � 2 0 ) models.
e identification results of categories of the 50 block samples by the two models are shown in Figure 13.  4. e 150 training samples were also predicted by the two models and are listed in Table 4.
According to the results in Table 4, it can be seen that the categories of coal and carbonaceous shale samples based on GRB-KPCA-SVM have higher PAs comparing with PAs of PCA-SVM.
e PA of 92% for nature block samples is acceptable, and GRB-KPCA-SVM is a more practicable and ideal coal-carbonaceous shale identification method.And the GRB-KPCA-SVM model provides a valid reference for coal-carbonaceous shale identification using VIS-NIR reflectance spectroscopy.

Conclusions
In order to solve the problem of misclassification between coal and carbonaceous shale using VIS-NIR reflectance spectroscopy in the field of coal-gangue identification of hyperspectral remote sensing of coal mine, spectral reflectance characterization of six typical bituminite, and carbonaceous shale samples of China in the wavelength range of 350-2500 nm were studied.In the six samples, two carbonaceous shales without obvious absorption valleys have similar reflectance curves to the three bituminite samples.rough simulation experiment of the effect of bituminite content on reflectance curve of mixed powder of ultralow ash bituminite and clay, it was found that absorption valleys near 1400 nm, 1900 nm, and 2200 nm of the mixed powder become not obvious when bituminite content is more than 30%.It could be deduced that when carbon concentration of the carbonaceous shale is more than about 30%, it is difficult to identify carbonaceous shale from bituminite according to the waveform of spectral reflectance curve of the sample in the VIS-NIR region.Eight wavelength regions including 400-550 nm, 650-750 nm, 820-920 nm, 1000-1100 nm, 1350-1450 nm, 1850-1950 nm, 2150-2250 nm, and 2300-2500 nm which are related to the main mineral and organic compositions of the samples were selected.CR spectra of 150 other coal and carbonaceous shale samples of China in the eight  Journal of Spectroscopy selected wavelength regions were used to establish two coal-carbonaceous shale prediction models: PCA-SVM and GRB-KPCA-SVM.e two models were verified through predicting categories of 50 powder samples and 50 nature block samples.
e GRB-KPCA-SVM model with prediction accuracies of 94% and 92% for powder and nature block prediction samples, respectively, is a more practicable and ideal coal-carbonaceous shale identification method.e basic principles and models in this paper provide references for coal-carbonaceous shale identification using VIS-NIR reflectance spectroscopy in coal mines.

Figure 2 :Figure 3 :
Figure 2: Reflectance spectra acquisition of powder sample in 0.5 mm grain size using a field spectrometer and a halogen lamp light source.
of all maximum average PAs of all with parameters σ e best σ, c, γ and trained GRB-KPCA-SVM model (b)

Figure 5 :
Figure 5: Spectral reflectance curves of the six representative samples in Figure 1 in the wavelength range of 350-2500 nm.

Figure 6 :
Figure 6: Spectral reflectance curves of 21 mixed powder samples of dry bituminite and clay in 0.5 mm grain size under different mixing ratios in the wavelength range of 350-2500 nm.
lists cumulative contribution rates of the first ten PCs of CR spectra of the training set in the eight regions through PCA. e cumulative contribution rate of the first three PCs is 90.0403%, and the 3D spatial distribution of the first three PCs of the 150 training samples through PCA is exhibited in Figure 8(a).

Figure 8 (
a) indicates that a large part of rock (carbonaceous shale) balls are mixed in coal balls.

Figure 7 :
Figure 7: Spectral reflectance curves (a) and continuum-removed spectral curves (b) of 200 modeling samples in the wavelength range of 350-2500 nm.

Figure 9 :Figure 8 :
Figure 9: Process of determining the optimal parameter σ of GRB-KPCA: the variation trend of the maximum average PA of categories of the validation set through GRB-KPCA-SVM with the change of σ.

3. 4 . 3 .
Performance of Prediction Models.Identification accuracies of categories of the 50 powder prediction samples and the 50 extra block samples are listed in Table

Figure 10 :
Figure 10: (a) 3D spatial distribution of the first three PCs of CR spectra of the 150 training samples in the eight regions in Table 1 through GRB-KPCA (σ � 63); (b) process of searching for the optimal penalty parameter c and radial basis function variance parameter c of GRB-KPCA-SVM model through cross validation with the first ten PCs as input variables.

Figure 12 :Figure 13 :
Figure 12: Spectral reflectance curves (a) and continuum-removed spectral curves (b) of 50 nature block samples in the wavelength range of 350-2500 nm.

Table 2 :
X-ray fluorescence (XRF) and proximate analysis of samples in Figure1.

Table 3 :
Cumulative contribution rates of the first ten PCs of continuum-removed spectra of the 150 training samples in the eight selected wavelength range in Table1.

Table 4 :
Identification accuracies of categories of all bituminite and carbonaceous shale samples collected from Malan coal mine area in Shanxi province of China by the prediction models.