王 杰

(1.瓦斯災(zāi)害監(jiān)控與應(yīng)急技術(shù)國家重點(diǎn)實(shí)驗(yàn)室，重慶 400037；2.中煤科工集團(tuán)重慶研究院有限公司，重慶 400039)

目前，煤礦瓦斯災(zāi)害治理的根本途徑是煤層瓦斯抽采。其中，巷道法預(yù)抽鄰近層及卸壓層瓦斯是瓦斯治理中有效的區(qū)域性措施，該方法需在底(頂)板巷中向突出煤層施工大量的鉆孔，通過鉆孔工程的擾動使煤體部分卸壓[1-3]。根據(jù)瓦斯抽采鉆孔的施工要求，底抽巷穿層鉆孔需進(jìn)入煤層底板后再穿出煤層頂板不少于0.5 m。而現(xiàn)階段的全液壓回轉(zhuǎn)鉆機(jī)，在施工過程中無法準(zhǔn)確地判斷鉆頭層位，只能依靠人的經(jīng)驗(yàn)知識來判斷鉆孔是否到達(dá)設(shè)計(jì)層位，一般根據(jù)鉆孔泥漿顏色、排渣特征、鉆機(jī)壓力變化情況等信息來進(jìn)行巖層或煤層識別，該方法受經(jīng)驗(yàn)積累程度、信息滯后、干擾因素多的影響，判斷結(jié)果往往具有一定的偏差、延遲和誤差。若不能及時識別煤層的頂、底板，則可能出現(xiàn)鉆孔不到位，導(dǎo)致抽采空白帶的出現(xiàn)。

在工程實(shí)踐中，基于鉆進(jìn)參數(shù)的地層信息識別方法廣泛用于地層評價和界面識別[4-5]。在煤礦巷道圍巖分類方面，提出了一種基于鉆進(jìn)參數(shù)的核模糊Cmeans 聚類(Kernel Fuzzy C-means，KFCM)算法的含煤地層巖性模糊識別方法，主要用于巷道圍巖特征識別[6]。研發(fā)的煤礦巷道液壓錨桿鉆機(jī)隨鉆參數(shù)系統(tǒng)可通過計(jì)算破巖比能對巖石的力學(xué)性能進(jìn)行識別[7]。在煤礦井下自動鉆進(jìn)方面，提出了一種基于地層識別的自動鉆進(jìn)控制方法，系統(tǒng)通過比功法自動識別3 種地層并匹配相應(yīng)的鉆進(jìn)參數(shù)[8]。在地質(zhì)鉆進(jìn)方面，采用數(shù)字式鉆孔過程監(jiān)測(Drilling Process Monitor，DPM)系統(tǒng)對香港地區(qū)復(fù)雜風(fēng)化花崗巖地層的界面進(jìn)行了成功識別[9]。在鑿巖掘進(jìn)方面，提出一種基于PCA-BP（Principal Component Analysis-Back Propagation）神經(jīng)網(wǎng)絡(luò)隨鉆參數(shù)巖性智能感知方法，通過鑿巖鉆機(jī)試驗(yàn)臺采集的6 種敏感鉆進(jìn)參數(shù)對泥巖、砂巖和軟弱層進(jìn)行了識別[10]。在石油鉆井方面，將3 種機(jī)器學(xué)習(xí)模型應(yīng)用于基于鉆進(jìn)參數(shù)的地層巖性和頂部實(shí)時預(yù)測，研究成果表明優(yōu)化后的人工神經(jīng)網(wǎng)絡(luò)模型(Artificial Neural Network，ANN)的精度性能優(yōu)于自適應(yīng)神經(jīng)模糊推理系統(tǒng)(Adaptive Neuro-Fuzzy Inference System，ANFIS)和函數(shù)神經(jīng)網(wǎng)絡(luò)(Functional Neural Network，F(xiàn)NN)[11]。上述研究表明了基于鉆進(jìn)參數(shù)進(jìn)行地層識別的可行性，但多數(shù)應(yīng)用于巷道圍巖分類、地質(zhì)鉆進(jìn)、鑿巖掘進(jìn)和油氣井等工程領(lǐng)域，在瓦斯抽采鉆孔工程中的相關(guān)文獻(xiàn)較少。

由于煤礦井下環(huán)境的特殊性和鉆孔類型的多樣性，鉆進(jìn)參數(shù)的獲取和分析難度較大。因此，在煤礦井下鉆進(jìn)參數(shù)采集系統(tǒng)設(shè)計(jì)和煤巖界面識別方法中開展相關(guān)研究工作，并將機(jī)器學(xué)習(xí)算法應(yīng)用于煤巖分類中，獲得鉆進(jìn)參數(shù)和地層之間的分類關(guān)系，建立基于鉆進(jìn)參數(shù)的煤巖界面識別模型，以期實(shí)現(xiàn)煤巖界面的自動識別，為鉆孔層位判定提供一套技術(shù)裝備。

1 基本原理

1.1 鉆進(jìn)參數(shù)與破碎比功

對于回轉(zhuǎn)切削的成孔方式，鉆進(jìn)過程包括回轉(zhuǎn)和推進(jìn)兩個主要做功部分，鉆進(jìn)參數(shù)主要有鉆桿轉(zhuǎn)速、回轉(zhuǎn)扭矩、推進(jìn)力(鉆壓)和推進(jìn)速度(鉆速)4 個基本參數(shù)。研究表明某些巖石的物理性質(zhì)與鉆進(jìn)參數(shù)之間存在一定的響應(yīng)關(guān)系，在理想環(huán)境下，鉆進(jìn)參數(shù)的差異可明顯區(qū)分不同的模擬地層，可較準(zhǔn)確地對地層進(jìn)行界面識別[12-14]。

在實(shí)際應(yīng)用中，由于地質(zhì)條件的復(fù)雜性及鉆孔過程的隨機(jī)性，鉆機(jī)的工作參數(shù)是動態(tài)變化的，難以用準(zhǔn)確的閾值范圍來區(qū)分不同地層類型。而能量指標(biāo)綜合考慮了各種因素的影響大小，有效的避免了上述不利因素。在施工設(shè)備和技術(shù)水平固定的條件下，破碎單位體積巖石所需的能量與巖石物理力學(xué)性質(zhì)之間存在一定的相關(guān)性，利用能量參數(shù)可較為精確的表征鉆進(jìn)過程中物理參數(shù)與地層信息之間的響應(yīng)關(guān)系[13]。在各種能量指標(biāo)中，破碎比功是近年來被許多科學(xué)家認(rèn)可的比較合理的一種方法，即通過破碎單位體積巖石消耗的能量評價巖石破碎的難易程度。其計(jì)算公式[15-16]為：

式中：E為破碎比功，MPa；Fd為推進(jìn)軸向力，N；Az為鉆孔橫截面積，m2；n為鉆桿轉(zhuǎn)速，r/min；Tr為鉆機(jī)回轉(zhuǎn)扭矩，N·m；vd為推進(jìn)速度，m/min。

從式(1)可以看出，當(dāng)鉆孔橫截面積確定后，破碎比功主要由4 個基本參數(shù)決定，是鉆進(jìn)參數(shù)的一個派生參數(shù)。因此，只需要在鉆孔過程中采集鉆進(jìn)參數(shù)，就可得到頂板巖層和煤層的破碎比功值，近似反映當(dāng)前鉆進(jìn)地層的強(qiáng)度大小，從而識別鉆進(jìn)地層。

1.2 支持向量機(jī)分類原理

支持向量機(jī)(SVM)是機(jī)器學(xué)習(xí)中的一種分類方法，本質(zhì)上是一種按監(jiān)督學(xué)習(xí)方式對數(shù)據(jù)進(jìn)行二元分類的廣義線性分類器。其基本思想是：對于線性可分的數(shù)據(jù)集，在原空間尋找兩類樣本的最優(yōu)分類超平面，該超平面不僅能正確分類每一個樣本，且使得每一類樣本中距離超平面最近的數(shù)據(jù)到超平面的距離盡可能遠(yuǎn)。當(dāng)線性不可分時，通過使用非線性映射將低維度輸入空間的樣本映射到高維度空間使其變?yōu)榫€性可分，從而在該映射空間中尋找最優(yōu)分類超平面[17-18]。其分類原理如圖1 所示。

圖1 支持向量機(jī)分類原理Fig.1 Classification principle of SVM

圖1 中，實(shí)心和空心圓點(diǎn)分別代表兩類二維數(shù)據(jù)，分割它們的任一超平面可用線性方程描述為：

式中：w為法向量，決定了超平面的方向；x為二維數(shù)據(jù)構(gòu)成的樣本向量；b為位移量，決定了超平面與原點(diǎn)的距離，w和b這兩組參數(shù)可通過訓(xùn)練得到。

對于訓(xùn)練樣本(xi，yi)，定義yi=1 時為正樣本，yi=-1 時為負(fù)樣本，則可用以下公式表達(dá)：

距離超平面最近的幾個樣本點(diǎn)使得yi(wTxi+b)=1成立，它們被稱為“支持向量”，支持向量到超平面的距離之和被稱為“間隔”。由解析幾何中的距離公式可得到間隔d的表達(dá)為：

則優(yōu)化目標(biāo)為求間隔d最大值時對應(yīng)w和b，顯然，最大化間隔等于最小化||w||，也等價于最小化||w||2，加上約束條件后，求解的優(yōu)化問題可以寫成：

1.3 煤巖界面識別流程

煤巖界面的識別工作包括煤巖分類預(yù)測和界面位置確定兩方面需求。整個過程的識別流程如圖2 所示。其中，界面位置的確定由鉆機(jī)的工況識別系統(tǒng)實(shí)現(xiàn)鉆桿數(shù)量的自動記錄，通過鉆機(jī)液壓系統(tǒng)的各種條件判斷和邏輯規(guī)則識別進(jìn)、退桿工況，在進(jìn)桿時增加鉆桿數(shù)量，實(shí)現(xiàn)鉆孔深度的記錄，具體原理在本文中不再詳細(xì)贅述。

圖2 基于鉆進(jìn)參數(shù)的煤巖界面識別流程Fig.2 Coal-rock interface recognition process based on drilling parameters

煤巖分類預(yù)測采用機(jī)器學(xué)習(xí)的方法進(jìn)行自動識別，基于鉆機(jī)工作時壓力、扭矩和相關(guān)能量的動力傳遞過程，在鉆機(jī)每根鉆桿的鉆進(jìn)過程中，對鉆機(jī)的4 個鉆進(jìn)參數(shù)和破碎比功開展隨鉆檢測與自動采集。在此基礎(chǔ)上，對鉆進(jìn)數(shù)據(jù)進(jìn)行整理和標(biāo)記工作，構(gòu)建基于鉆進(jìn)參數(shù)的樣本數(shù)據(jù)庫和基于支持向量機(jī)的煤巖分類識別模型。通過模型訓(xùn)練和參數(shù)調(diào)整不斷提高預(yù)測準(zhǔn)確率和泛化能力，直到獲得理想的輸出為止，最后再根據(jù)此最優(yōu)的模型訓(xùn)練結(jié)果去識別未知的鉆進(jìn)參數(shù)。

2 煤巖界面識別系統(tǒng)設(shè)計(jì)

2.1 系統(tǒng)總體方案設(shè)計(jì)

為了實(shí)現(xiàn)煤巖界面識別過程的自動化，設(shè)計(jì)了基于鉆進(jìn)參數(shù)的鉆孔煤巖界面識別系統(tǒng)，整套系統(tǒng)由數(shù)據(jù)感知層、采集層和分析層組成。其中，數(shù)據(jù)感知層主要用于測量和變送鉆機(jī)參數(shù)，是該系統(tǒng)的基礎(chǔ)部分，由各類物理量傳感器和相應(yīng)的電氣線路及液壓管路組成。數(shù)據(jù)采集層主要用于鉆進(jìn)參數(shù)的采集與計(jì)算、鉆進(jìn)工況識別、鉆孔深度記錄和破碎比功曲線繪制，是整個系統(tǒng)的核心部分，由數(shù)據(jù)采集儀、工控計(jì)算機(jī)和相應(yīng)的數(shù)據(jù)采集程序組成。數(shù)據(jù)分析層主要用于樣本數(shù)據(jù)收集、數(shù)據(jù)學(xué)習(xí)、煤巖識別和結(jié)果輸出，是該系統(tǒng)的人機(jī)交互部分，由專用軟件和輸入輸出設(shè)備組成。整個系統(tǒng)的結(jié)構(gòu)框架如圖3 所示。

圖3 基于鉆進(jìn)參數(shù)的煤巖界面識別系統(tǒng)架構(gòu)Fig.3 Architecture of coal-rock interface recognition system based on drilling parameters

2.2 鉆機(jī)數(shù)據(jù)采集系統(tǒng)設(shè)計(jì)

鉆機(jī)數(shù)據(jù)采集系統(tǒng)主要用于采集鉆機(jī)液壓系統(tǒng)工作參數(shù)，同時完成數(shù)據(jù)計(jì)算、顯示和記錄等功能[19]。該系統(tǒng)由礦用隔爆兼本質(zhì)安全型電控箱、礦用本質(zhì)安全型傳感器、礦用防爆手機(jī)和配套附件等組成，其中數(shù)據(jù)采集器、工控計(jì)算機(jī)和各類電氣輔件均安裝在礦用隔爆兼本質(zhì)安全型電控箱內(nèi)，從而滿足煤礦井下防爆環(huán)境的使用要求；數(shù)據(jù)采集器通過安全隔離電路完成各類傳感器的模擬量和頻率量信號采集，再交由工控計(jì)算機(jī)進(jìn)行數(shù)據(jù)處理和顯示，所有的人機(jī)交互操作通過礦用防爆手機(jī)完成，其電氣原理如圖4 所示。

圖4 鉆機(jī)數(shù)據(jù)采集系統(tǒng)電氣原理Fig.4 Electrical principle of drilling rig data acquisition system

2.3 煤巖界面識別軟件功能設(shè)計(jì)

煤巖界面識別軟件是基于鉆進(jìn)參數(shù)的數(shù)據(jù)分析類軟件，運(yùn)行在工控計(jì)算機(jī)的操作系統(tǒng)上。因此，在功能上應(yīng)具備數(shù)據(jù)采集、鉆進(jìn)參數(shù)計(jì)算和數(shù)據(jù)分析與輸出三大功能。其中，數(shù)據(jù)采集模塊負(fù)責(zé)與數(shù)據(jù)采集器通信獲取傳感器輸出信號，并實(shí)時顯示傳感器當(dāng)前數(shù)值。鉆進(jìn)參數(shù)計(jì)算模塊主要是讀取配置參數(shù)、根據(jù)傳感器數(shù)值實(shí)時計(jì)算各類鉆進(jìn)參數(shù)。數(shù)據(jù)分析與輸出模塊可完成樣本標(biāo)記、模型訓(xùn)練、煤巖識別結(jié)果輸出等功能。整個軟件的工作流程如圖5 所示，工作界面如圖6 所示。

圖5 煤巖界面識別軟件工作流程Fig.5 Workflow of coal-rock interface recognition software

圖6 煤巖界面識別軟件工作界面Fig.6 User interface of coal-rock interface recognition software

3 現(xiàn)場試鉆參數(shù)采集

3.1 試驗(yàn)地點(diǎn)及鉆孔設(shè)計(jì)

煤巖界面識別系統(tǒng)的現(xiàn)場應(yīng)用在河南省鶴壁中泰礦業(yè)有限公司實(shí)施，鉆孔施工巷道為3116 上底抽巷，施工位置距二煤層底板約40 m，煤層走向132°～270°、傾向42°～160°、傾角4°～18°，平均煤層厚度8.1 m。預(yù)測瓦斯含量在6～8 m3/t，瓦斯壓力0.60～0.74 MPa，煤的堅(jiān)固性系數(shù)為0.8～1.5。

本次試驗(yàn)共設(shè)計(jì)4 個鉆孔，1-3 號鉆孔為試驗(yàn)鉆孔，4 號鉆孔為測試鉆孔，各鉆孔的設(shè)計(jì)參數(shù)見表1。

表1 鉆孔設(shè)計(jì)參數(shù)Table 1 Design parameters of boreholes

3.2 鉆機(jī)數(shù)據(jù)采集系統(tǒng)安裝

鉆進(jìn)參數(shù)的獲取采用間接測量法，即通過液壓系統(tǒng)相關(guān)回路的流量和壓力參數(shù)來計(jì)算對應(yīng)的鉆進(jìn)參數(shù)。因此，鉆孔施工前，需將對應(yīng)的傳感器安裝在鉆機(jī)的液壓系統(tǒng)中，使用的傳感器類型、數(shù)量及安裝位置見表2。

表2 傳感器參數(shù)及安裝位置Table 2 Sensor parameters and positions

安裝完畢后，為了確保傳感器均能正常工作，進(jìn)行了數(shù)據(jù)采集測試，如圖7 所示。

圖7 傳感器數(shù)據(jù)采集測試Fig.7 Sensor data acquicition testing

3.3 鉆進(jìn)參數(shù)采集

現(xiàn)場施工采用ZYWL-4000 型礦用全液壓鉆機(jī)，鉆孔過程中，控制鉆機(jī)轉(zhuǎn)速約65 r/min，鉆壓約14 MPa，所有鉆進(jìn)參數(shù)由鉆機(jī)數(shù)據(jù)采集系統(tǒng)按設(shè)定的采樣頻率(1 Hz)自動記錄，再用其平均值作為當(dāng)前鉆桿的鉆進(jìn)參數(shù)，試驗(yàn)現(xiàn)場如圖8 所示。

圖8 煤巖界面識別系統(tǒng)試驗(yàn)現(xiàn)場Fig.8 Experimental site of coal-rock interface identification

在1-4 號鉆孔的每個鉆進(jìn)階段，鉆機(jī)數(shù)據(jù)采集系統(tǒng)實(shí)時顯示傳感器采集參數(shù)和鉆孔深度，鉆進(jìn)參數(shù)中的轉(zhuǎn)速和鉆壓兩個參數(shù)由鉆機(jī)液壓系統(tǒng)控制，基本處于穩(wěn)定狀態(tài)；回轉(zhuǎn)扭矩、推進(jìn)速度和破碎比功3 個參數(shù)隨地層和孔深的變化發(fā)生改變。鉆孔結(jié)束后繪制出這3 個參數(shù)隨鉆孔深度的變化規(guī)律，如圖9 所示。

圖9 鉆進(jìn)參數(shù)隨鉆孔深度的變化規(guī)律Fig.9 Variation law of drilling parameters with drilling depth

從圖9 中可以看出以下規(guī)律。

(1)回轉(zhuǎn)扭矩、推進(jìn)速度、破碎比功3 個鉆進(jìn)參數(shù)的變化過程和地層之間具有明顯的對應(yīng)關(guān)系，這些參數(shù)在煤巖界面處均會出現(xiàn)明顯的漲落，與煤巖分類結(jié)果之間具有較大的相關(guān)性。

(2)由參數(shù)變化過程可以判斷出：1 號鉆孔的實(shí)際鉆孔深度為28 m，其中19～27 m 為煤層段，2—4 號鉆孔的深度和煤層段可同理得出，該結(jié)果與現(xiàn)場有經(jīng)驗(yàn)的司鉆人員判斷一致。

(3)整個鉆進(jìn)過程中，除開孔時的數(shù)據(jù)外，巖層中鉆進(jìn)的回轉(zhuǎn)扭矩在1 200～2 800 N·m；煤層中鉆進(jìn)時基本在1 200 N·m，但2 號鉆孔的回轉(zhuǎn)扭矩存在較大波動。

(4)在穩(wěn)壓鉆進(jìn)模式下，推進(jìn)速度會產(chǎn)生明顯變化，在巖層中推進(jìn)速度基本在0.3～1.0 m/min，在煤層中該值則為1.2～2.0 m/min，差異明顯。

(5)在煤層中鉆進(jìn)的破碎比功基本處于30～60 MPa；在巖層中多數(shù)大于90 MPa，具有較強(qiáng)的區(qū)分性，但在巖石中開孔時的數(shù)值和煤層段存在重疊，可能出現(xiàn)誤判。

4 模型訓(xùn)練與測試

4.1 支持向量機(jī)分類算法及訓(xùn)練步驟

采用Python 機(jī)器學(xué)習(xí)庫Scikit-Learn 中的支持向量機(jī)算法(SVM)進(jìn)行鉆進(jìn)參數(shù)的煤巖分類學(xué)習(xí)，相關(guān)步驟[18]如下。

(1)對樣本數(shù)據(jù)集進(jìn)行分類標(biāo)記和標(biāo)簽的數(shù)字化處理，用不同的數(shù)字表示各種分類。

(2)數(shù)據(jù)集分割：將樣本數(shù)據(jù)集分割為訓(xùn)練集和測試集，訓(xùn)練集用于模型訓(xùn)練和參數(shù)調(diào)整；測試集用于模型準(zhǔn)確度評價。

(3)對數(shù)據(jù)集進(jìn)行標(biāo)準(zhǔn)化或歸一化處理。

(4)選擇SVM 算法的核函數(shù)和對應(yīng)的超參數(shù)。

(5)通過訓(xùn)練集獲得最優(yōu)超參數(shù)下的分類模型。

(6)利用測試集測試模型的分類準(zhǔn)確度，觀察預(yù)測結(jié)果。

具體流程如圖10 所示。

圖10 支持向量機(jī)(SVM)模型訓(xùn)練步驟Fig.10 Model training steps of SVM

4.2 樣本數(shù)據(jù)集的建立與分割

根據(jù)鉆進(jìn)參數(shù)與地層信息之間的響應(yīng)過程，將3 個與煤巖分類結(jié)果具有較大相關(guān)性的鉆進(jìn)參數(shù)作為模型訓(xùn)練的特征參數(shù)。確定特征參數(shù)后將1-4 號鉆孔的數(shù)據(jù)按順序進(jìn)行煤巖分類標(biāo)記和標(biāo)簽數(shù)字化，用數(shù)字“0”和“1”分別作為巖石和煤層分類的數(shù)據(jù)標(biāo)簽，得到用于模型訓(xùn)練的樣本數(shù)據(jù)集，見表3。

表3 樣本數(shù)據(jù)集Table 3 Sample data sets

根據(jù)模型訓(xùn)練需要，再將數(shù)據(jù)集分割為訓(xùn)練集和測試集，各數(shù)據(jù)集的樣本信息見表4，訓(xùn)練集和測試集的數(shù)量比例約為3∶1。

表4 數(shù)據(jù)集分類信息Table 4 Classification information on data sets

4.3 模型參數(shù)搜尋

模型參數(shù)搜尋主要是選擇SVM 算法的核函數(shù)和對應(yīng)的超參數(shù)，常用的核函數(shù)有線性核函數(shù)，徑向基核函數(shù)(高斯核函數(shù))和多項(xiàng)式核函數(shù)等，其中線性核函數(shù)主要用于線性可分的數(shù)據(jù)集，徑向基是SVM 默認(rèn)的核函數(shù)，多項(xiàng)式核函數(shù)可以將低維的輸入空間映射到高維的特征空間。每種核函數(shù)都有其對應(yīng)的調(diào)整參數(shù)，其中有兩個重要的參數(shù)，即懲罰系數(shù)C和影響系數(shù)γ；C為分類誤差的寬容度，數(shù)值越大則對誤差的容忍程序越??；γ定義了單個樣本對于整個超平面的影響大小，數(shù)值越大支持向量越少，越小則支持向量越多[20-21]。

模型參數(shù)的確定通過訓(xùn)練集獲得，為了解訓(xùn)練集各樣本點(diǎn)的分布規(guī)律，首先繪制了數(shù)據(jù)散點(diǎn)圖，用黑色圓點(diǎn)表示巖石樣本點(diǎn)，灰色圓點(diǎn)表示煤層樣本點(diǎn)，則分布規(guī)律如圖11 所示。該圖可以表示3 個維度的信息：x軸為樣本點(diǎn)的推進(jìn)速度參數(shù)，y軸為回轉(zhuǎn)扭矩參數(shù)，圓點(diǎn)的面積由每個樣本點(diǎn)的鉆進(jìn)參數(shù)計(jì)算的破碎比功參數(shù)確定，面積越大則代表該樣本點(diǎn)消耗的能量越高。

圖11 訓(xùn)練集數(shù)據(jù)的散點(diǎn)分布Fig.11 Scatter plots of the training set data

從圖11 中可以看出，巖石和煤層樣本點(diǎn)之間具有明顯的分界線，說明特征參數(shù)和分類結(jié)果之間為線性可分的關(guān)系，因此，使用SVM 中的線性核函數(shù)即可滿足分類要求。

線性核函數(shù)可以調(diào)整的超參數(shù)只有懲罰系數(shù)C一個參數(shù)，采用網(wǎng)格搜索函數(shù)(GridSearchCV)在指定的參數(shù)范圍內(nèi)按步長依次調(diào)整參數(shù)并進(jìn)行比較，直到從所有的參數(shù)中找到在驗(yàn)證集上精度最高的參數(shù)。設(shè)置搜索范圍0.1～10.0，步長50，最終獲得C的最佳參數(shù)為0.1。

4.4 模型精度測試

對于小樣本數(shù)據(jù)，為了最大程度地驗(yàn)證模型的精確度和泛化能力，在模型訓(xùn)練時常采用交叉驗(yàn)證法劃分訓(xùn)練集和驗(yàn)證集。考慮到訓(xùn)練集分類數(shù)據(jù)的不平衡性，采用分層K折交叉驗(yàn)證法，首先將訓(xùn)練集樣本劃分成K個大小相等的樣本子集，并確保每一個子集中的樣本類別分布與整個數(shù)據(jù)集的類別分布一致；依次遍歷這K個子集，每次把當(dāng)前子集作為驗(yàn)證集，其余K-1 個子集作為新的訓(xùn)練集，進(jìn)行模型的訓(xùn)練和評估；最后把K次評估指標(biāo)的平均值作為最終的評估指標(biāo)。

本次訓(xùn)練進(jìn)行了K=10 的交叉驗(yàn)證實(shí)驗(yàn)，具體分組結(jié)果如圖12 所示。

圖12 訓(xùn)練集10 折交叉驗(yàn)證分組結(jié)果Fig.12 Cross validation results of training set (K=10)

圖中橫坐標(biāo)代表樣本號0～88，縱坐標(biāo)代表10 次分組情況，用分組1～10 表示。其中，黑色部分為劃分的訓(xùn)練集樣本號，綠色部分為驗(yàn)證集樣本號。分組確定后，再用同一測試集對所有分組進(jìn)行模型評估，得到最終的測試結(jié)果見表5。

表5 10 折交叉驗(yàn)證模型精度Table 5 Model accuracy of cross validation model (K=10)

從表中可以看出，只有分組7 對應(yīng)的模型驗(yàn)證集精度為0.888 9，其他9 組模型的驗(yàn)證集精度全部為1.0，驗(yàn)證集的平均精度為0.988 9；測試集的平均精度為1.0。

分組7 雖在驗(yàn)證集中出現(xiàn)了分類誤差，但在測試集中分類全部正確。通過數(shù)據(jù)定位，找到分組7 中驗(yàn)證集分類錯誤的樣本點(diǎn)序號為60，該樣本為2 號鉆孔第31 m 時的鉆進(jìn)數(shù)據(jù)，即該鉆孔穿出煤層再次進(jìn)入巖層中的第一根鉆桿的鉆進(jìn)參數(shù)，該數(shù)據(jù)的分類錯誤主要是由于煤巖界面處的樣本點(diǎn)數(shù)量太少導(dǎo)致的。

5 結(jié)論

a.開發(fā)的基于鉆進(jìn)參數(shù)的煤巖界面識別系統(tǒng)，初步解決了底抽巷穿層鉆孔煤巖界面識別缺少技術(shù)手段的問題。通過分析鉆進(jìn)參數(shù)與煤巖地層之間的響應(yīng)關(guān)系，找到了與地層變化相關(guān)性較大的3 個特征參數(shù)：回轉(zhuǎn)扭矩、推進(jìn)速度和破碎比功。其中，破碎比功是綜合考慮鉆進(jìn)參數(shù)影響的能量指標(biāo)，在煤巖識別方面，較其他兩個特征參數(shù)具有更強(qiáng)的參考意義，但在鉆孔開孔段的適用性較差。

b.在煤巖分類方面，采用的支持向量機(jī)(SVM)分類算法在3 個特征值的二元分類中具有非常高的準(zhǔn)確性，特別適合小樣本數(shù)據(jù)集的訓(xùn)練，使用線性核函數(shù)的分類模型在測試集中取得了100%的準(zhǔn)確度，該結(jié)果也驗(yàn)證了特征參數(shù)和分類結(jié)果之間是線性可分的。

c.經(jīng)現(xiàn)場試驗(yàn)驗(yàn)證，該煤巖界面識別系統(tǒng)可準(zhǔn)確采集和記錄鉆機(jī)的鉆進(jìn)參數(shù)，并通過特征參數(shù)分析獲得當(dāng)前鉆進(jìn)地層的煤巖分類結(jié)果，從而指導(dǎo)司鉆人員進(jìn)行鉆孔層位判定，確保鉆孔達(dá)標(biāo)，避免抽采空白帶的產(chǎn)生。

d.現(xiàn)有煤巖分類模型的學(xué)習(xí)結(jié)果僅對應(yīng)用煤礦的識別效果最佳。當(dāng)在新鉆場應(yīng)用(地質(zhì)條件改變)或改變鉆機(jī)轉(zhuǎn)速和鉆壓時，為保證識別的準(zhǔn)確度，應(yīng)增加2～3 個鉆孔的鉆進(jìn)參數(shù)進(jìn)行數(shù)據(jù)學(xué)習(xí)和模型更新。后續(xù)可開展半監(jiān)督或無監(jiān)督機(jī)器學(xué)習(xí)算法在煤巖分類識別中的應(yīng)用，從而減少人的參與程度。