楊 恩，王世博，葛世榮，向 陽(yáng)

(中國(guó)礦業(yè)大學(xué) 機(jī)電工程學(xué)院，江蘇 徐州 221116)

由于煤炭形成沉積過(guò)程的影響，煤層頂板與煤層交界處多為富含碳質(zhì)物質(zhì)的碳質(zhì)頁(yè)巖或黑色頁(yè)巖[1-2]，常使得頂板巖石外觀顏色上呈現(xiàn)出類似于相鄰煤層煤炭的黑色[3-4]。由于黑色的頂板巖石與煤層煤炭外觀色澤相近，在采煤機(jī)滾筒截割煤層的工作過(guò)程中，往往難以準(zhǔn)確區(qū)分煤炭與周圍巖石，導(dǎo)致煤巖界面處黑色的巖石被滾筒同時(shí)截割下落混入煤中，此現(xiàn)象是導(dǎo)致開采煤炭中矸石增多的主要因素之一，也是無(wú)人化開采煤巖識(shí)別研究中需要解決的主要問(wèn)題[5-7]。因此研究能夠準(zhǔn)確區(qū)分煤炭和黑色頂板巖石的方法至關(guān)重要。中國(guó)國(guó)家標(biāo)準(zhǔn)GB/T 5751《中國(guó)煤炭分類》和國(guó)際標(biāo)準(zhǔn)ISO 11760 Classification of coals對(duì)煤炭給出了相同的定義標(biāo)準(zhǔn)，即煤炭是主要由植物遺體經(jīng)煤化作用轉(zhuǎn)化而成的富含碳的固體可燃有機(jī)沉積巖，含有一定量的礦物質(zhì)，其灰分產(chǎn)率小于或等于50%。據(jù)此，如果能夠準(zhǔn)確、快速地得到原位固體煤炭與周圍黑色巖石的灰分產(chǎn)率，即可判定煤巖類型，進(jìn)而控制采煤機(jī)滾筒的截割位置。

近紅外反射光譜技術(shù)是一種高信噪比、實(shí)時(shí)、低廉、針對(duì)原位物質(zhì)的光譜分析技術(shù)，不僅在煤礦和巖礦遙感等定性識(shí)別探測(cè)領(lǐng)域得到了應(yīng)用[8-11]，而且已在煤炭、礦物、土壤等成分定量分析檢測(cè)領(lǐng)域中得到了較多應(yīng)用[12-14]，其中灰分產(chǎn)率分析是基于近紅外反射吸收光譜的煤質(zhì)成分定量分析研究中的主要內(nèi)容，已有國(guó)內(nèi)外眾多學(xué)者建立了多種基于近紅外反射吸收光譜的煤炭灰分產(chǎn)率預(yù)測(cè)方法。ANDRES等[15-16]基于主成分分析法(PCA)和偏最小二乘法(PLS)建立了煤炭灰分指標(biāo)與其近紅外光譜之間的定量數(shù)學(xué)關(guān)系模型，此模型對(duì)未知煤樣灰分具有較高的預(yù)測(cè)精度和較快的預(yù)測(cè)速度。KIM等[17]在煤炭的近紅外光譜中選擇了1 680,1 942,2 100,2 180,2 300 nm等信息含量最高的波長(zhǎng)點(diǎn)，建立了灰分的多元回歸分析(MRA)預(yù)測(cè)模型，在煤炭灰分的在線分析中取得了有效的結(jié)果。趙凱等[18]采用主成分分析、自組織映射神經(jīng)網(wǎng)絡(luò)-模糊C均值聚類(SOM-FCM)雙層聚類方法有效地優(yōu)化煤炭樣本，降低了煤炭近紅外光譜數(shù)據(jù)維數(shù)，并搭建了基于GA-BP神經(jīng)網(wǎng)絡(luò)的煤炭灰分預(yù)測(cè)模型，有效提高了模型學(xué)習(xí)精度。雷萌[19]針對(duì)基于機(jī)器學(xué)習(xí)的煤炭灰分近紅外光譜分析，研究了多種提高建模光譜數(shù)據(jù)質(zhì)量和預(yù)測(cè)模型性能的學(xué)習(xí)算法。基于近紅外反射吸收光譜分析的煤炭灰分預(yù)測(cè)的研究為煤礦井下煤巖灰分產(chǎn)率的原位獲取提供了啟示。

然而以上文獻(xiàn)中的各種煤炭灰分預(yù)測(cè)模型，在煤炭光譜采集時(shí)，大部分試樣為粉碎后的均質(zhì)粉末試樣，并且近紅外光譜儀光纖探頭均為靠近粉末試樣獲取煤炭光譜。在常見煤礦井下現(xiàn)場(chǎng)工作面煤壁、原位煤炭和頂板巖石均為塊狀試樣，同時(shí)采煤機(jī)機(jī)身離煤層頂端較遠(yuǎn)，頂板與煤層交界巖石離采煤機(jī)機(jī)身距離常見為1.5～3.0 m，鮮有相關(guān)文獻(xiàn)對(duì)于實(shí)時(shí)獲取原位塊狀煤炭和頂板黑色巖石灰分產(chǎn)率的研究。然而地物反射光譜成分定量探測(cè)和定量遙感已取得了長(zhǎng)足的發(fā)展[20-21]，也為煤礦井下煤巖灰分等成分定量探測(cè)提供了參考和思路。獲取原位煤巖的反射光譜并研究煤巖感知識(shí)別方法，為降低含矸量、煤巖識(shí)別等研究提供依據(jù)。

筆者收集了同一井下綜采工作面的碳質(zhì)頁(yè)巖和煙煤塊狀試樣75個(gè)，在實(shí)驗(yàn)室利用近紅外光譜儀和光纖準(zhǔn)直鏡采集了距離塊狀煤巖試樣表面3 m處的近紅外波段(1 000～2 500 nm)反射光譜?；诿簬r灰分產(chǎn)率及其反射光譜建立了兩類支持向量機(jī)(SVM)煤巖預(yù)測(cè)模型。

1 煤巖試樣及其反射光譜采集

從山西馬蘭煤礦同一井下綜采工作面現(xiàn)場(chǎng)收集了頂板與煤層交界處的碳質(zhì)頁(yè)巖以及煙煤大塊狀試樣75個(gè)，試樣收集后放入自封袋保存，其中煤樣42個(gè)，巖樣33個(gè)，試樣外觀均呈黑色，較為相似。

如圖1所示，在實(shí)驗(yàn)室采集每塊試樣較平整表面的反射光譜。為減小被測(cè)煤巖物質(zhì)雙向反射特性的影響[20,22-23]，使用4個(gè)50 W鹵鎢聚光燈按上下左右90°間隔排列，均以45°入射角照射所選平整表面中心，形成直徑約10 cm，光照度約20 000 lux的圓形光斑。使用荷蘭Avantes AvaSpec近紅外光譜儀進(jìn)行光譜采集，其光譜波長(zhǎng)為1 000～2 500 nm，光譜分辨率為8.9 nm，采集每條光譜曲線耗時(shí)約為1 ms，同時(shí)采用PTFE材料白板進(jìn)行反射參照標(biāo)定。使用直徑2.5 cm準(zhǔn)直鏡進(jìn)行光束準(zhǔn)直，準(zhǔn)直鏡連接于Y型石英光纖合并端，分支端的一端連接光譜儀進(jìn)行光譜采集，另一端連接650 nm激光光源以進(jìn)行目標(biāo)試樣對(duì)準(zhǔn)。準(zhǔn)直鏡固定于三維旋轉(zhuǎn)臺(tái)以進(jìn)行探測(cè)目標(biāo)對(duì)準(zhǔn)調(diào)節(jié)，調(diào)節(jié)旋轉(zhuǎn)臺(tái)保持準(zhǔn)直鏡軸線垂直對(duì)準(zhǔn)試樣表面的光斑中心，準(zhǔn)直鏡與光斑中心保持距離l=3 m。準(zhǔn)直鏡視場(chǎng)角(FOV)θ=0.12°，鏡頭直徑d=2.5 cm，如圖1(b)的實(shí)驗(yàn)原理示意所示，光譜儀所采集光譜為準(zhǔn)直鏡視場(chǎng)角所形成圓臺(tái)在試樣表面底面圓區(qū)域的平均反射光譜，底面圓的直徑D、面積A分別約為

(1)

(2)

計(jì)算機(jī)通過(guò)USB3.0與光譜儀連接，實(shí)時(shí)顯示并可同時(shí)采集反射光譜曲線。每個(gè)試樣的探測(cè)區(qū)域共采集10條反射光譜曲線，進(jìn)行斷點(diǎn)修正和去噪，然后計(jì)算其均值作為該試樣表面的反射光譜。

圖1 塊狀煤巖試樣反射光譜采集實(shí)驗(yàn)Fig.1 Reflectance spectra acquisition experiment of block coal or rock sample

在每塊試樣進(jìn)行反射光譜采集時(shí)，標(biāo)記光譜采集表面的7.5 cm2圓形區(qū)域。在所有試樣完成反射光譜采集后，將每塊試樣的圓形光譜采集區(qū)域打磨下表面至0.5 cm深度內(nèi)的試樣，碾成粉末，將粉末放至灰化馬弗爐進(jìn)行灰分產(chǎn)率分析，獲得每塊試樣反射光譜采集區(qū)域的灰分產(chǎn)率。

2 煤巖灰分產(chǎn)率與其反射光譜關(guān)系

2.1 煤巖光譜反射率曲線

為有效去除由于儀器性能以及環(huán)境等因素影響所產(chǎn)生的光譜曲線噪聲，筆者采用Savitzky-Golay(SG)卷積平滑算法對(duì)實(shí)驗(yàn)采集光譜數(shù)據(jù)進(jìn)行濾波去噪處理[24-25]，其核心求解算法為

(3)

圖2為75個(gè)塊狀煤巖試樣表面在近紅外波段(1 000～2 500 nm)的經(jīng)斷點(diǎn)修正、濾波去噪、取均值后的光譜反射率曲線。受黏土礦物、石英等礦物影響，反射光譜曲線在近紅外波段明顯地表現(xiàn)出水平至正向斜率傾斜的整體變化趨勢(shì)[4]，小部分試樣光譜曲線受原位試樣中水分影響，1 400和1 900 nm波段附近出現(xiàn)吸收特征。由圖2可以看出一大部分煤巖試樣的反射光譜波形相似并且重疊交叉，從波形難以有效區(qū)分。

圖2 75個(gè)煤巖試樣的光譜反射率曲線Fig.2 Reflectance spectra curves of 75 coal and rock samples

為了研究能夠根據(jù)原位塊狀煤巖的反射光譜有效區(qū)分煤巖的方法并建立識(shí)別模型，從75條試樣光譜曲線中隨機(jī)選擇50條為訓(xùn)練樣本，其中煤樣28個(gè)，巖樣22個(gè)，其余25條光譜曲線為測(cè)試樣本，其中煤樣14個(gè)，巖樣11個(gè)。

2.2 煤巖反射光譜預(yù)處理

為有效增強(qiáng)光譜反射吸收特征強(qiáng)度，利用一階微分(FD)、二階微分(SD)、連續(xù)統(tǒng)去除(CR)、標(biāo)準(zhǔn)正態(tài)變量變換(SNV)4種方法預(yù)處理圖2中所有75條煤巖反射光譜曲線[14,24]。4種預(yù)處理方法的核心求解算法如下:

(1)一階微分

(4)

其中，xi,λ+Δλ1，xi,λ-Δλ1為光譜向量xi,λ中的元素，即編號(hào)為i的樣本在波長(zhǎng)λ+Δλ1，λ-Δλ1處的光譜值;Δλ1為一階微分的波長(zhǎng)間隔;x′i,λ為波長(zhǎng)點(diǎn)λ的一階微分光譜反射率。本文選擇Δλ=3 nm。

(2)二階微分

(5)

其中，x′i,λ+Δλ2，x′i,λ-Δλ2分別指編號(hào)為i的樣本在波長(zhǎng)λ+Δλ2，λ-Δλ2處的光譜值(光譜反射率)對(duì)波長(zhǎng)的一階微分;Δλ2為二階微分的波長(zhǎng)間隔;x″i,λ為波長(zhǎng)點(diǎn)λ的二階微分光譜反射率。

(3)連續(xù)統(tǒng)去除

(6)

其中，xi,λ為光譜向量xi,λ中的元素，即編號(hào)為i的樣本在波長(zhǎng)λ處的光譜值;zi,λ為波長(zhǎng)點(diǎn)λ的包絡(luò)線值;cri,λ為波長(zhǎng)點(diǎn)λ的連續(xù)統(tǒng)去除光譜反射率。

(4)標(biāo)準(zhǔn)正態(tài)變量變換

(7)

xi,λ,A=-lgxi,λ

(8)

(9)

結(jié)合文獻(xiàn)[4]，為了進(jìn)行光譜預(yù)處理分析說(shuō)明，在圖2中選擇了灰分產(chǎn)率在50%附近的5個(gè)煤樣和5個(gè)巖樣的光譜曲線為代表，如圖3所示，其中，煤1～煤5的灰分產(chǎn)率分別為39%,40%,48%,41%,43%，巖1～巖5的灰分產(chǎn)率分別為60%,57%,59%,52%,53%。由圖3可知，這10條反射光譜曲線的整體輪廓波形較為相似。圖4分別為此10條反射光譜曲線經(jīng)一階微分、二階微分、連續(xù)統(tǒng)去除、標(biāo)準(zhǔn)正態(tài)變量變換預(yù)處理后的光譜曲線。由圖4可知，經(jīng)4種算法預(yù)處理后，光譜曲線整體波形發(fā)生變化，反射吸收特征有所增強(qiáng)，但各個(gè)試樣的光譜曲線依舊波形趨勢(shì)相似并且重疊交叉，難以有效區(qū)分。然而，因物質(zhì)光譜的反射吸收特征是其物質(zhì)成分的“指紋”，含有豐富的物質(zhì)成分信息，反射吸收特征的增強(qiáng)更利于煤巖光譜與其物質(zhì)成分關(guān)系的研究。

圖3 代表性煤巖試樣的光譜反射率曲線Fig.3 Spectral reflectance curves of representative coal and rock samples

圖4 代表性煤巖試樣光譜反射率曲線經(jīng)4種方法預(yù)處理后的光譜曲線Fig.4 Spectral curves of representative coal and rock samples whose spectral reflectance curves had been preprocessed by four methods

2.3 煤巖灰分產(chǎn)率與其反射光譜的相關(guān)性

煤巖灰分主要來(lái)源于其黏土礦物、石英等礦物成分，礦物和有機(jī)物成分影響其反射光譜吸收特征，物質(zhì)成分含量與光譜值的相關(guān)性分析是一種有效的獲取特征波長(zhǎng)點(diǎn)的方法[24,26]。為此選取前述50個(gè)訓(xùn)練樣本，對(duì)其灰分產(chǎn)率與去噪后反射光譜以及前述4種算法預(yù)處理后光譜各波長(zhǎng)點(diǎn)處的光譜值進(jìn)行相關(guān)性分析，其中波長(zhǎng)λ處的相關(guān)系數(shù)rλ為

(10)

式中,X′i,λ為波長(zhǎng)λ處的光譜值;yi為第i個(gè)樣本的灰分產(chǎn)率。

此50個(gè)煤巖樣本的灰分產(chǎn)率與去噪后反射光譜以及4種算法預(yù)處理后光譜在各波長(zhǎng)點(diǎn)處的相關(guān)系數(shù)如圖5所示，圖中同時(shí)標(biāo)出了每種相關(guān)系數(shù)曲線的最大值。由圖5可知，連續(xù)統(tǒng)去除預(yù)處理方法的最大相關(guān)系數(shù)出現(xiàn)在波長(zhǎng)1 698 nm處，為0.777，大于其他4種相關(guān)系數(shù)曲線的最大值，其中二階微分預(yù)處理方法的最大相關(guān)系數(shù)最小。對(duì)連續(xù)統(tǒng)去除1 698 nm處進(jìn)行線性、指數(shù)、冪函數(shù)回歸分析，如圖6所示。

由圖6中的各類回歸方法可知，灰分產(chǎn)率回歸分析的決定系數(shù)均較低，因此常見線性和非線性回歸函數(shù)均不能有效進(jìn)行灰分預(yù)測(cè)。然而，連續(xù)統(tǒng)去除預(yù)處理方法最大相關(guān)系數(shù)波長(zhǎng)點(diǎn)1 698 nm位于1 700 nm光譜帶附近，1 700 nm為煤巖中主要有機(jī)基團(tuán)倍頻和合頻的典型譜帶[27]，如芳香結(jié)構(gòu)中C—H伸縮振動(dòng)3 030 cm-1與脂肪側(cè)鏈中CH2對(duì)稱伸縮振動(dòng)2 850 cm-1的合頻[28]。本實(shí)驗(yàn)煤巖主要有機(jī)成分為芳香環(huán)狀結(jié)構(gòu)化合物、脂肪側(cè)鏈結(jié)構(gòu)化合物[29-30]，煤巖灰分產(chǎn)率間接取決于這些主要有機(jī)成分的含量，與此同時(shí)，連續(xù)統(tǒng)去除預(yù)處理方法1 698 nm波長(zhǎng)點(diǎn)附近左右區(qū)間相關(guān)系數(shù)曲線連續(xù)，為最大相關(guān)系數(shù)曲線段，因此考慮采用1 698 nm波長(zhǎng)點(diǎn)附近左右區(qū)間連續(xù)統(tǒng)去除光譜建立原位煤巖灰分產(chǎn)率預(yù)測(cè)模型。本文選取的連續(xù)統(tǒng)去除光譜波長(zhǎng)點(diǎn)區(qū)間為1 693～1 703 nm，所含波長(zhǎng)點(diǎn)對(duì)應(yīng)相關(guān)系數(shù)見表1。

3 基于反射光譜的煤巖預(yù)測(cè)模型

3.1 煤巖預(yù)測(cè)模型的建立

根據(jù)前述對(duì)50個(gè)煤巖樣本主要有機(jī)成分的特征譜帶、灰分產(chǎn)率-光譜相關(guān)波長(zhǎng)點(diǎn)的分析，借鑒定量遙感和化學(xué)計(jì)量學(xué)的思路[21,24,31]，選取表1中的11個(gè)波長(zhǎng)點(diǎn)處連續(xù)統(tǒng)去除預(yù)處理光譜為訓(xùn)練數(shù)據(jù)，建立煤巖灰分產(chǎn)率機(jī)器學(xué)習(xí)預(yù)測(cè)模型。與其他機(jī)器學(xué)習(xí)算法相比，支持向量機(jī)(SVM)算法更適合于本實(shí)驗(yàn)的小訓(xùn)練樣本、多維度成分因素、非線性關(guān)系問(wèn)題[21]，為此本文建立煤巖灰分產(chǎn)率支持向量回歸(SVR)預(yù)測(cè)模型。同時(shí)以表1中11個(gè)波長(zhǎng)點(diǎn)處連續(xù)統(tǒng)去除預(yù)處理光譜為訓(xùn)練數(shù)據(jù)，50個(gè)試樣的煤巖類型為預(yù)測(cè)目標(biāo)數(shù)據(jù)，建立煤巖類型支持向量分類(SVC)預(yù)測(cè)模型，以對(duì)比回歸和分類兩類支持向量機(jī)方法對(duì)基于反射光譜原位煤巖識(shí)別的有效性。

圖5 50個(gè)訓(xùn)練集煤巖試樣灰分產(chǎn)率與反射光譜及預(yù)處理后光譜的相關(guān)系數(shù)曲線Fig.5 Correlation coefficient curves between ash yields and reflectance spectra and preprocessed spectra of 50 training coal and rock samples

圖6 圖5(d)中1 698 nm處灰分產(chǎn)率與光譜值的回歸分析Fig.6 Regression analyses between ash yields and spectral values at 1 698 nm inFig.5(d)

表1 50個(gè)訓(xùn)練集煤巖試樣灰分產(chǎn)率與其連續(xù)統(tǒng)去除預(yù)處理光譜的最大相關(guān)系數(shù)區(qū)間
Table 1 Interval of maximum correlation coefficientsbetween ash yields and continuum removal-preprocessedspectra of 50 training coal and rock samples

波長(zhǎng)點(diǎn)/nm相關(guān)系數(shù)1 6930.7701 6940.7711 6950.7731 6960.7731 6970.7741 6980.7771 6990.7761 7000.7721 7010.7691 7020.7691 7030.768

SVR基本原理為將光譜數(shù)據(jù)映射到高維空間并建立最優(yōu)界面，使得所有數(shù)據(jù)到此界面的距離最小，此時(shí)最優(yōu)界面回歸函數(shù)為

(11)

(12)

(13)

其最優(yōu)界面求解過(guò)程與SVR類似，亦借助Lagrange變換法[32]，各參數(shù)意義與SVR類似，其中,yic為第i個(gè)訓(xùn)練樣本的類型代號(hào)，本文使用0代表煤，1代表巖。

對(duì)于SVR和SVC煤巖預(yù)測(cè)模型的建立，本文采用5折交叉驗(yàn)證求解最優(yōu)γ和c，兩個(gè)參數(shù)的搜索范圍均取2-10～210，指數(shù)步長(zhǎng)均取0.5，由50個(gè)煤巖訓(xùn)練試樣的表1中11個(gè)波長(zhǎng)點(diǎn)處連續(xù)統(tǒng)去除預(yù)處理光譜值及其灰分產(chǎn)率、煤巖類型為訓(xùn)練數(shù)據(jù)，求解最優(yōu)γ和c的過(guò)程如圖7所示。由圖7(a)可知，煤巖灰分產(chǎn)率SVR預(yù)測(cè)模型訓(xùn)練中，煤巖試樣5折交叉驗(yàn)證的驗(yàn)證集煤巖灰分產(chǎn)率平均最小均方差(MSE)為0.002 182，此時(shí)最優(yōu)解為c=24.5，γ=2-8。由圖7(b)可知，煤巖類型SVC預(yù)測(cè)模型訓(xùn)練中，煤巖試樣5折交叉驗(yàn)證的驗(yàn)證集煤巖類型平均最大預(yù)測(cè)精度(PA)為94%，此時(shí)最優(yōu)解為c=24，γ=2-4.5。此時(shí)，由前述50個(gè)煤巖樣本為訓(xùn)練集，建立了煤巖灰分產(chǎn)率支持向量回歸(SVR)預(yù)測(cè)模型和煤巖類型支持向量分類(SVC)預(yù)測(cè)模型。

3.2 煤巖預(yù)測(cè)模型的預(yù)測(cè)效果評(píng)估

為進(jìn)一步驗(yàn)證以上所建立的SVR和SVC兩種煤巖預(yù)測(cè)模型的有效性，取前述25個(gè)測(cè)試煤巖試樣連續(xù)統(tǒng)去除預(yù)處理后光譜在表1中11個(gè)波長(zhǎng)點(diǎn)處的光譜值，用上述兩種模型進(jìn)行灰分產(chǎn)率和煤巖類型預(yù)測(cè)，預(yù)測(cè)結(jié)果如圖8所示。表2列出了SVR模型對(duì)25個(gè)測(cè)試煤巖試樣灰分產(chǎn)率預(yù)測(cè)的均方根誤差(RMSE)、決定系數(shù)(R2)，以及兩類SVM方法對(duì)測(cè)試試樣的預(yù)測(cè)精度(PA)，3個(gè)評(píng)價(jià)參數(shù)的計(jì)算方法如下

(14)

(15)

圖7 交叉驗(yàn)證求解兩類支持向量機(jī)預(yù)測(cè)模型的最優(yōu)參數(shù)Fig.7 Processes of searching for the optimal parameters of two types of SVM prediction models by cross validations

圖8 兩類支持向量機(jī)預(yù)測(cè)模型對(duì)測(cè)試煤巖樣本的預(yù)測(cè)結(jié)果Fig.8 Prediction results of test coal and rock samples by two types of SVM prediction models

式中，yi,yip分別為第i個(gè)測(cè)試試樣的實(shí)測(cè)、預(yù)測(cè)灰分產(chǎn)率;N為測(cè)試試樣類型被正確預(yù)測(cè)的個(gè)數(shù)。

同時(shí)，使用兩類SVM預(yù)測(cè)模型對(duì)50個(gè)訓(xùn)練樣本的同類光譜數(shù)據(jù)進(jìn)行了預(yù)測(cè)，預(yù)測(cè)結(jié)果數(shù)據(jù)同時(shí)也列于表2中。

表2 兩類支持向量機(jī)預(yù)測(cè)模型對(duì)煤巖樣本的預(yù)測(cè)結(jié)果數(shù)據(jù)
Table 2 Prediction result data of coal and rock samples bytwo types of SVM prediction models

預(yù)測(cè)模型測(cè)試試樣RMSER2PA/%訓(xùn)練試樣RMSER2PA/%SVR0.050.88920.030.9298SVC——96——98

由表2可知，測(cè)試試樣的預(yù)測(cè)精度均達(dá)到90%以上，訓(xùn)練試樣均達(dá)到95%以上，SVC類型預(yù)測(cè)方法略高于SVR灰分產(chǎn)率回歸預(yù)測(cè)方法。由圖8(a)可知，測(cè)試試樣預(yù)測(cè)灰分產(chǎn)率與實(shí)測(cè)灰分產(chǎn)率整體趨勢(shì)是一致的，同時(shí)SVR灰分產(chǎn)率回歸RMSE為5%，故本文SVR模型對(duì)天然原位煤巖試樣固體表面灰分產(chǎn)率的預(yù)測(cè)較為精確。

本實(shí)驗(yàn)在仿真模型預(yù)測(cè)測(cè)試試樣時(shí)，還測(cè)試了兩種預(yù)測(cè)模型對(duì)每一個(gè)測(cè)試集試樣光譜的預(yù)測(cè)時(shí)間，均在0.1 s以內(nèi)，加上每個(gè)試樣單條光譜曲線采集時(shí)間約為1 ms，即每個(gè)試樣光譜采集和類型預(yù)測(cè)總耗時(shí)不超過(guò)0.1 s。因此，以本實(shí)驗(yàn)中50個(gè)煤巖訓(xùn)練樣本在表1中11個(gè)波長(zhǎng)點(diǎn)處連續(xù)統(tǒng)去除預(yù)處理光譜值，及其灰分產(chǎn)率、煤巖類型建立的SVR,SVC兩種原位煤巖感知識(shí)別模型對(duì)同一煤礦試樣能夠有效識(shí)別，并且能夠?qū)υ嚇宇愋蛯?shí)時(shí)快速預(yù)測(cè)。

此外，考慮到光纖準(zhǔn)直鏡探頭將來(lái)應(yīng)用于煤礦井下時(shí)應(yīng)具有防爆功能，針對(duì)防爆玻璃窗對(duì)近紅外反射信號(hào)的衰減問(wèn)題，可采取增加煤巖探測(cè)區(qū)域光照度的措施。本實(shí)驗(yàn)煤巖表面探測(cè)區(qū)域光照度只有20 000 lux，GOETZ等[13]通過(guò)設(shè)計(jì)專用聚光探頭提高了光照度，其在指定波段光照度為40倍日光照度，在100 000 lux以上，實(shí)現(xiàn)了傳送帶上矸石中礦物含量的在線精確測(cè)定。同時(shí)，煤巖探測(cè)區(qū)域光照度的增加不僅可增加探測(cè)距離，而且還可增加反射光譜的信噪比。

4 結(jié) 論

(1)依據(jù)常見頂板高度，距煙煤和碳質(zhì)頁(yè)巖原位塊狀試樣表面3 m處通過(guò)光纖準(zhǔn)直鏡采集的近紅外波段(1 000～2 500 nm)反射光譜，通過(guò)SG去噪，經(jīng)FD,SD,CR,SNV四種方法預(yù)處理后，與試樣灰分產(chǎn)率的最大相關(guān)系數(shù)出現(xiàn)在1 698 nm波長(zhǎng)點(diǎn)處，由連續(xù)統(tǒng)去除(CR)預(yù)處理方法獲得，且其波長(zhǎng)點(diǎn)位于與煤巖主要有機(jī)成分有關(guān)的1 700 nm光譜帶附近。

(2)煤巖訓(xùn)練樣本的近紅外反射光譜在經(jīng)CR預(yù)處理后，以1 693～1 703 nm波長(zhǎng)區(qū)間的11個(gè)波長(zhǎng)點(diǎn)處預(yù)處理光譜值，及訓(xùn)練樣本灰分產(chǎn)率、煤巖類型建立了SVR煤巖灰分產(chǎn)率預(yù)測(cè)模型、SVC煤巖類型預(yù)測(cè)模型，兩種原位煤巖感知識(shí)別模型對(duì)同一煤礦試樣識(shí)別率均達(dá)到90%以上，并且能夠?qū)υ嚇宇愋蛯?shí)時(shí)快速預(yù)測(cè)。

(3)本實(shí)驗(yàn)所用原位煤巖感知預(yù)測(cè)方法，為煤礦井下采煤工作面中煤層與頂板間煤巖界面的遠(yuǎn)距離、準(zhǔn)確、快速、低成本識(shí)別或界定提供了參考。