楊 恩，王世博，葛世榮，向 陽

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

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

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

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

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

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

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

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

(1)

(2)

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

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

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

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

2.1 煤巖光譜反射率曲線

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

(3)

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

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

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

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

為有效增強光譜反射吸收特征強度，利用一階微分(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,λ中的元素，即編號為i的樣本在波長λ+Δλ1，λ-Δλ1處的光譜值;Δλ1為一階微分的波長間隔;x′i,λ為波長點λ的一階微分光譜反射率。本文選擇Δλ=3 nm。

(2)二階微分

(5)

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

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

(6)

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

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

(7)

xi,λ,A=-lgxi,λ

(8)

(9)

結(jié)合文獻[4]，為了進行光譜預(yù)處理分析說明，在圖2中選擇了灰分產(chǎn)率在50%附近的5個煤樣和5個巖樣的光譜曲線為代表，如圖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ā)生變化，反射吸收特征有所增強，但各個試樣的光譜曲線依舊波形趨勢相似并且重疊交叉，難以有效區(qū)分。然而，因物質(zhì)光譜的反射吸收特征是其物質(zhì)成分的“指紋”，含有豐富的物質(zhì)成分信息，反射吸收特征的增強更利于煤巖光譜與其物質(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)性

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

(10)

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

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

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

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

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

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

圖5 50個訓(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個訓(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

波長點/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ù)到此界面的距離最小，此時最優(yōu)界面回歸函數(shù)為

(11)

(12)

(13)

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

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

3.2 煤巖預(yù)測模型的預(yù)測效果評估

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

(14)

(15)

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

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

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

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

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

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

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

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

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

4 結(jié) 論

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

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

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