包從望,江偉,劉永志,張彩紅
(六盤(pán)水師范學(xué)院礦業(yè)與土木工程學(xué)院,貴州六盤(pán)水 553000)
采煤機(jī)作為煤礦生產(chǎn)的核心開(kāi)采設(shè)備,是礦井高效安全開(kāi)采的前提,一旦采煤機(jī)發(fā)生故障,將導(dǎo)致礦井停產(chǎn)甚至發(fā)生重大事故。因此,研究深部煤層中,復(fù)雜環(huán)境下,采煤機(jī)的運(yùn)行狀態(tài)感知技術(shù)與故障智能診斷是煤炭智能開(kāi)采、綠色開(kāi)采亟需解決的關(guān)鍵問(wèn)題之一。滾動(dòng)軸承作為采煤機(jī)的主要傳動(dòng)支撐件,是采煤機(jī)穩(wěn)定運(yùn)行的“重要關(guān)節(jié)”,其故障率占采煤機(jī)故障的61%。因此高效、準(zhǔn)確地實(shí)現(xiàn)深部煤層下采煤機(jī)軸承的故障辨識(shí)是采煤機(jī)故障診斷的研究重點(diǎn)。而采煤機(jī)常在強(qiáng)噪聲環(huán)境下運(yùn)行,致使軸承的故障信號(hào)呈現(xiàn)非平穩(wěn)、非線性且故障信息微弱的特點(diǎn)。為此,研究變工況、強(qiáng)噪聲環(huán)境下采煤機(jī)的軸承故障診斷具有重要意義。
卷積神經(jīng)網(wǎng)絡(luò)作為深度學(xué)習(xí)的一種,能從復(fù)雜的數(shù)據(jù)中提取性能良好的故障特征,目前已成為主流的特征提取方法之一。彭雪瑩等基于圖卷積網(wǎng)絡(luò),實(shí)現(xiàn)軸承在不同工況下的服役故障診斷。LIANG等基于特征融合與并行卷積神經(jīng)網(wǎng)絡(luò)組建了故障診斷模型,解決了小數(shù)據(jù)集時(shí)神經(jīng)網(wǎng)絡(luò)過(guò)擬合的問(wèn)題,實(shí)現(xiàn)滾動(dòng)軸承在小數(shù)據(jù)集時(shí)的故障診斷。WANG等從軸承的振動(dòng)信號(hào)和聲音信號(hào)中提取故障特征,并在一維卷積神經(jīng)網(wǎng)絡(luò)中實(shí)現(xiàn)特征融合,實(shí)現(xiàn)了融合振動(dòng)傳感器和聲音傳感器的軸承故障診斷方法。吳晨芳等將改進(jìn)的LeNet-5模型融入到卷積神經(jīng)網(wǎng)絡(luò),根據(jù)信號(hào)特點(diǎn)自動(dòng)調(diào)整卷積網(wǎng)絡(luò)的輸入尺寸,通過(guò)卷積網(wǎng)絡(luò)自適應(yīng)提取有限轉(zhuǎn)速情形下故障特征,從而實(shí)現(xiàn)故障診斷。
卷積神經(jīng)網(wǎng)絡(luò)在提取故障過(guò)程中需要輸入大量帶標(biāo)簽的數(shù)據(jù)樣本,而采煤機(jī)長(zhǎng)期工作在變工況環(huán)境中,其服役軸承的故障數(shù)據(jù)量極少,且受工作載荷、環(huán)境噪聲等的干擾,導(dǎo)致不同工況下采集到的故障訓(xùn)練信號(hào)也存在差異,因此直接利用卷積網(wǎng)絡(luò)對(duì)采煤機(jī)進(jìn)行故障診斷是行不通的。為此,文中引入故障模型的遷移理論,實(shí)現(xiàn)采煤機(jī)故障信號(hào)少數(shù)據(jù)的故障診斷。遷移學(xué)習(xí)能夠?qū)崿F(xiàn)源域故障模型遷移至目標(biāo)域,目標(biāo)域在已建模型的基礎(chǔ)上只需要稍加微調(diào)即可達(dá)到較好的學(xué)習(xí)效果。WANG等提出了一種基于可變形的CNN-DLSTM滾動(dòng)軸承故障診斷遷移學(xué)習(xí)方法,利用小數(shù)據(jù)集故障樣本對(duì)預(yù)訓(xùn)練模型進(jìn)行微調(diào),從而實(shí)現(xiàn)了不同工況下小數(shù)據(jù)樣本的軸承故障診斷。QIAN等提出一種基于卷積自編碼器的新型深度遷移學(xué)習(xí)網(wǎng)絡(luò),該網(wǎng)絡(luò)無(wú)需進(jìn)行數(shù)據(jù)標(biāo)記,即可實(shí)現(xiàn)變載荷和強(qiáng)噪聲干擾下的齒輪箱遷移故障診斷。ZOU等提出一種基于深度卷積域?qū)咕W(wǎng)絡(luò)的遷移學(xué)習(xí)方法,增加了特征的聚合性和源域不同類別特征間的邊界,從而實(shí)現(xiàn)變工況下軸承故障診斷。CAO等基于深度遷移學(xué)習(xí)算法,以偽分類最大平均差異作為優(yōu)化目標(biāo),實(shí)現(xiàn)領(lǐng)域自適應(yīng)智能故障診斷,并驗(yàn)證了診斷模型的穩(wěn)定性。ZHENG等以其他類型滾動(dòng)軸承的故障數(shù)據(jù)集作為輔助,提出卷積神經(jīng)網(wǎng)絡(luò)下的遷移學(xué)習(xí)網(wǎng)絡(luò),實(shí)現(xiàn)了不同滾動(dòng)軸承之間的故障遷移診斷,具有準(zhǔn)確率高、泛化能力強(qiáng)的特點(diǎn)。
由于采煤機(jī)服役軸承故障數(shù)據(jù)匱乏,深度學(xué)習(xí)網(wǎng)絡(luò)無(wú)法實(shí)現(xiàn)充分訓(xùn)練,也無(wú)法直接用于采煤機(jī)故障診斷,本文作者基于卷積神經(jīng)網(wǎng)絡(luò),設(shè)計(jì)了變工況下的采煤機(jī)軸承故障診斷遷移學(xué)習(xí)方法。針對(duì)采煤機(jī)實(shí)際故障樣本缺乏的問(wèn)題,通過(guò)充足的人工模擬故障,在卷積神經(jīng)網(wǎng)絡(luò)中學(xué)習(xí)故障知識(shí),利用采煤機(jī)故障診斷實(shí)驗(yàn)臺(tái)的數(shù)據(jù)對(duì)卷積網(wǎng)絡(luò)進(jìn)行微調(diào),并將學(xué)習(xí)后的故障模型遷移到采煤機(jī)故障診斷中,實(shí)現(xiàn)人工模擬故障下的故障診斷模型到現(xiàn)場(chǎng)采煤機(jī)故障診斷的遷移。
卷積神經(jīng)網(wǎng)絡(luò)已廣泛應(yīng)用于圖像處理中,基本結(jié)構(gòu)包括輸入層、卷積層、池化層、全連接層和輸出層,如圖1所示。通過(guò)卷積層和池化層計(jì)算,直接提取數(shù)據(jù)中的深層特征,免去了原始方法中人工提取幾何特征的繁瑣操作。
圖1 卷積神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)模型
卷積神經(jīng)網(wǎng)絡(luò)的經(jīng)典操作方式是將圖像進(jìn)行卷積計(jì)算和池化壓縮,通過(guò)卷積核將圖像中的特征以局部區(qū)域的形式提取出來(lái),其計(jì)算過(guò)程如式(1)所示
(1)
池化層壓縮卷積運(yùn)算的結(jié)果,將區(qū)域數(shù)據(jù)集內(nèi)的典型特征進(jìn)一步提取出來(lái),常用的方法有均值池化和最大池化函數(shù)。為凸顯采煤機(jī)故障時(shí)頻特征,采用最大池化的方法池化數(shù)據(jù)集。
網(wǎng)絡(luò)的最終輸出層將池化像素結(jié)果按順序一一排列構(gòu)成一維特征向量,形成全連接層,組成一維特征向量。
卷積神經(jīng)網(wǎng)絡(luò)有較好的圖像特征提取能力,但主要針對(duì)的對(duì)象是圖像,而采煤機(jī)的故障信號(hào)振動(dòng)信號(hào)為一維信號(hào),為充分利用卷積神經(jīng)網(wǎng)絡(luò)的特征提取能力,通常需將一維振動(dòng)信號(hào)以圖像的形式輸入到卷積神經(jīng)網(wǎng)絡(luò)中。常見(jiàn)的轉(zhuǎn)換方式有:直接時(shí)域振動(dòng)信號(hào),按時(shí)間順序排列振動(dòng)信號(hào)幅值形成圖像,快速傅里葉變換為頻譜圖,小波變換等。為充分保留振動(dòng)信號(hào)中的時(shí)頻信息,提高原始信號(hào)轉(zhuǎn)換為圖像過(guò)程中的特征信息量,結(jié)合小波變換對(duì)時(shí)頻信號(hào)的分析優(yōu)點(diǎn),采用小波變換方式預(yù)處理原始振動(dòng)信號(hào),獲得源域和目標(biāo)域樣本。當(dāng)振動(dòng)信號(hào)為()時(shí),小波變換可計(jì)算為
(2)
式中:為尺度因子;為平移因子;(·)為小波基函數(shù)。對(duì)連續(xù)采集到的軸承振動(dòng)信號(hào),其周期變化信號(hào)的特征是相似的,因此選擇cmor小波基。
小波變換后的時(shí)頻圖包含了各類故障信息,為進(jìn)一步凸顯故障特征,采用圖像增強(qiáng)手段豐富故障信息量,以提高故障識(shí)別效果,圖像預(yù)處理過(guò)程中采用簡(jiǎn)單有效的直方圖均值化,通過(guò)改變頻譜圖中的灰度值分布而改變各類故障信息之間的對(duì)比度。
假設(shè)頻譜圖的像素為,其灰度值范圍為[0,-1],則灰度直方圖可表述為
()=
(3)
式中:()為第個(gè)灰度級(jí)的直方圖;為第個(gè)灰度級(jí)對(duì)應(yīng)的個(gè)數(shù)。
直方圖均值化變換如式(4)所示,式中:round(·)為取整運(yùn)算,為灰度分布函數(shù),其計(jì)算為式(5)
=round[(-1),]
(4)
(5)
遷移學(xué)習(xí)算法屬于機(jī)器學(xué)習(xí),遷移方法可分為樣本、特征、模型、關(guān)系等方面的遷移。遷移過(guò)程是將不同領(lǐng)域之間的知識(shí)進(jìn)行相互轉(zhuǎn)化,即源域下的模型遷移到目標(biāo)域中,從而實(shí)現(xiàn)實(shí)驗(yàn)?zāi)M到現(xiàn)場(chǎng)設(shè)備之間的故障遷移。因采煤機(jī)故障診斷時(shí)的結(jié)構(gòu)參數(shù)未知,因此可將它作為目標(biāo)域,遷移學(xué)習(xí)的示意如圖2所示。
圖2 遷移學(xué)習(xí)結(jié)構(gòu)簡(jiǎn)圖
為度量從源域到目標(biāo)域遷移的效果,以最大均值差異(Maximum Mean Discrepancy,MMD)作為優(yōu)化準(zhǔn)則,觀測(cè)源域和目標(biāo)域數(shù)據(jù)在再生核希爾伯特空間的空間分布(Reproducing Kernel Hilbert Space, RKHS),可表示為
(6)
式中:={,, …,}、={,,…,}分別表示源域和目標(biāo)域的數(shù)據(jù)0,為源域數(shù),為目標(biāo)域數(shù);(·)為核函數(shù),通過(guò)核函數(shù)將源域和目標(biāo)域表述為RKHS。因高斯核函數(shù)對(duì)映射空間的核函數(shù)沒(méi)有要求,因此,文中選擇高斯核函數(shù)實(shí)現(xiàn)特征空間的映射,如式(7)所示:
(,)=e(-||-||)(2)
(7)
式中:為控制核函數(shù)徑向范圍的帶寬。
針對(duì)采煤機(jī)故障樣本不易獲取的特點(diǎn),將遷移學(xué)習(xí)的故障診斷方法引入到采煤機(jī)軸承故障診斷,即將大量的、容易獲取的源域數(shù)據(jù)樣本訓(xùn)練后的模型用于采煤機(jī)現(xiàn)場(chǎng)故障診斷中。如圖3所示,首先對(duì)源域數(shù)據(jù)樣本進(jìn)行信號(hào)預(yù)處理,包括數(shù)據(jù)歸一化、小波時(shí)頻變換、圖像尺寸自適應(yīng)和圖像增強(qiáng),得到二維時(shí)頻圖集樣本;然后,將圖集樣本輸入到卷積神經(jīng)網(wǎng)絡(luò)實(shí)現(xiàn)特征提取,得到特征空間,并將該特征模型遷移至源域中;最后,以采煤機(jī)故障診斷實(shí)驗(yàn)臺(tái)的數(shù)據(jù)作為目標(biāo)域,通過(guò)與源域相同的信號(hào)預(yù)處理得到二維時(shí)頻圖集,以該圖集作為遷移模型的樣本輸入,并通過(guò)降低學(xué)習(xí)率,以MMD作為優(yōu)化指標(biāo)實(shí)現(xiàn)權(quán)值微調(diào),從而實(shí)現(xiàn)故障診斷模型的遷移學(xué)習(xí)。最終通過(guò)SVM分類器驗(yàn)證模型的故障識(shí)別效果。
圖3 變工況下采煤機(jī)故障診斷模型
為使診斷模型有更好的遷移效果,將目標(biāo)域數(shù)據(jù)經(jīng)信號(hào)預(yù)處理后輸入到遷移模型中,并通過(guò)模型提取目標(biāo)域的特征空間。計(jì)算目標(biāo)域特征和源域特征之間的MMD值,根據(jù)經(jīng)驗(yàn),當(dāng)MMD值大于0.01時(shí),需重新初始化卷積層。因此,當(dāng)MMD值大于0.01時(shí),將模型的學(xué)習(xí)率降低為原學(xué)習(xí)率的1/2,從而實(shí)現(xiàn)模型權(quán)值的微調(diào)。與其他遷移算法不同,文中所提方法并未鎖定遷移后的模型,而是以MMD值作為優(yōu)化指標(biāo)對(duì)權(quán)值進(jìn)行微調(diào)。
文中的源域數(shù)據(jù)來(lái)源于凱斯西儲(chǔ)大學(xué)(Case Western Reserve University,CWRU)的軸承故障數(shù)據(jù)集,實(shí)驗(yàn)臺(tái)包含電機(jī)、故障軸承、加載及數(shù)據(jù)采集。選取6203型深溝球軸承的正常、內(nèi)圈、外圈和滾子故障4種運(yùn)行狀態(tài)作為源域數(shù)據(jù)集。各類運(yùn)行狀態(tài)在變工況下組成數(shù)據(jù)集,如表1所示,除了健康狀態(tài)以外,其余3種故障直徑分別為0.177 8 mm和0.533 4 mm,從數(shù)據(jù)集中交錯(cuò)滑動(dòng)取樣,得到每類工況下的480個(gè)樣本,每個(gè)樣本點(diǎn)數(shù)均為1 200個(gè)。當(dāng)采樣頻率為12 kHz、電機(jī)加載功率為2.2 kW、轉(zhuǎn)速為1 730 r/min、故障尺寸為0.177 8 mm時(shí)4種軸承的振動(dòng)信號(hào)如圖4所示。
表1 軸承振動(dòng)數(shù)據(jù)集
圖4 源域樣本振動(dòng)時(shí)域信號(hào)
圖5所示為井下采煤機(jī)開(kāi)采現(xiàn)場(chǎng),當(dāng)割到煤和巖石時(shí)軸承受到不同程度的振動(dòng),且截割力矩也有所不同。因采煤機(jī)現(xiàn)場(chǎng)故障數(shù)據(jù)樣本較少且不易獲取,針對(duì)采煤機(jī)開(kāi)采時(shí)的運(yùn)行特點(diǎn),搭建如圖6所示的實(shí)驗(yàn)臺(tái)模擬采煤機(jī)軸承故障信息,實(shí)驗(yàn)臺(tái)包含變頻調(diào)速電機(jī)、速度控制箱、磁力加載器及信號(hào)采集系統(tǒng)等,通過(guò)調(diào)節(jié)加載器的加載力矩模擬采煤機(jī)割煤時(shí)的隨機(jī)工況。軸承選用NU202,并利用激光在內(nèi)圈、外圈和滾動(dòng)體上加工直徑0.4 mm、深度為0.5 mm的故障,設(shè)置不同的電機(jī)轉(zhuǎn)速,目標(biāo)域的組成如表2所示,其中數(shù)據(jù)集A-C中有40個(gè)樣本為訓(xùn)練集微調(diào),20個(gè)樣本為測(cè)試集,而數(shù)據(jù)集D中的所有數(shù)據(jù)集不參與訓(xùn)練,只有20個(gè)樣本作為測(cè)試集。
圖5 采煤機(jī)截割部運(yùn)行環(huán)境
圖6 采煤機(jī)故障診斷實(shí)驗(yàn)臺(tái)
表2 目標(biāo)域數(shù)據(jù)集
選用cmor小波基,對(duì)時(shí)域振動(dòng)信號(hào)進(jìn)行小波變換,并將頻域圖像歸一化為600像素×600像素的樣本,如圖7所示。圖中可看出:4種軸承狀態(tài)存在較大的區(qū)別,因此可將其作為故障模型的輸入樣本,按照該處理方式將其余樣本經(jīng)小波變換得到圖像樣本。
圖7 軸承時(shí)頻圖
為進(jìn)一步增強(qiáng)樣本的故障特征,利用直方圖均衡化原理對(duì)時(shí)頻圖作圖像增強(qiáng)。如圖8所示,與原始時(shí)域圖相比,增強(qiáng)后的圖像有更明顯的故障特征,有利于模型訓(xùn)練與故障模式識(shí)別。
圖8 軸承時(shí)頻增強(qiáng)圖
基于TensorFlow平臺(tái)調(diào)用卷積神經(jīng)網(wǎng)絡(luò)工具箱,利用Python語(yǔ)言編程,在i7-10750H型處理器下實(shí)現(xiàn)模型的驗(yàn)證。選取源域中的每種運(yùn)行樣本10個(gè)作為測(cè)試集,即共有測(cè)試集160個(gè),所有測(cè)試樣本均隨機(jī)抽取,選用3層卷積層和4層池化層作為網(wǎng)絡(luò)模型,最大訓(xùn)練次數(shù)為60次,初始學(xué)習(xí)率設(shè)為0.01。利用源域數(shù)據(jù)集完成模型的訓(xùn)練,將訓(xùn)練后的模型遷移至采煤機(jī)故障診斷實(shí)驗(yàn)臺(tái)的目標(biāo)域,通過(guò)目標(biāo)域數(shù)據(jù)集對(duì)訓(xùn)練完后的模型以0.000 5的學(xué)習(xí)率進(jìn)行微調(diào)。
為驗(yàn)證變工況下采煤機(jī)故障診斷的遷移效果,對(duì)比分析該方法與以下方法的差異。方法一:未對(duì)故障振動(dòng)信號(hào)進(jìn)行預(yù)處理,直接將一維振動(dòng)信號(hào)輸入到卷積神經(jīng)網(wǎng)絡(luò),以同樣的模型遷移方法研究遷移效果;方法二:振動(dòng)信號(hào)經(jīng)小波變換變?yōu)闀r(shí)頻圖,但未經(jīng)圖像增強(qiáng)直接輸入到卷積神經(jīng)網(wǎng)絡(luò)中完成訓(xùn)練、測(cè)試與遷移;方法三:訓(xùn)練及測(cè)試樣本均經(jīng)過(guò)小波時(shí)域變換與圖像增強(qiáng)處理,僅鎖定遷移后的低層權(quán)值,未對(duì)所有權(quán)值進(jìn)行微調(diào);方法四:源域數(shù)據(jù)僅取自1 750 r/min轉(zhuǎn)速下、電機(jī)加載功率為1.5 kW的工況,數(shù)據(jù)集并非變工況數(shù)據(jù),其余與所提方法的處理方式一樣;方法五:對(duì)樣本數(shù)據(jù)集進(jìn)行時(shí)域變換和圖像增強(qiáng)處理實(shí)現(xiàn)網(wǎng)絡(luò)的訓(xùn)練和模型遷移,并以MMD作為優(yōu)化指標(biāo),利用目標(biāo)域的樣本對(duì)遷移后的模型進(jìn)行權(quán)值微調(diào),即所提方法。
為直觀體現(xiàn)所提方法的優(yōu)越性,比較每種方法在遷移學(xué)習(xí)權(quán)值微調(diào)過(guò)程中的學(xué)習(xí)精度,如圖9所示。結(jié)果顯示:方法一遷移有較大的波動(dòng),且平穩(wěn)時(shí)的精度不到75%;在方法二中,當(dāng)步數(shù)到達(dá)60步時(shí)的精度依然存在較大波動(dòng),且最大精度不足80%;方法三的精度較穩(wěn)定,但值始終未達(dá)到80%;方法四相比其他方法有較高的精度,且當(dāng)訓(xùn)練步數(shù)達(dá)50步以后偏向穩(wěn)定,精度可達(dá)92%,但遷移精度存在一定波動(dòng);方法五到第15步的時(shí)候精度已達(dá)95%,且精度較穩(wěn)定,最終可達(dá)100%的遷移精度。
圖9 模型的遷移精度對(duì)比
為進(jìn)一步說(shuō)明所提方法的有效性,分別以數(shù)據(jù)集A-D中的20個(gè)樣本作為驗(yàn)證數(shù)據(jù),驗(yàn)證幾種方法的遷移效果,其遷移識(shí)別精度如表3所示。方法一~方法四的精度均低于方法五,直接將二維振動(dòng)信號(hào)輸入到診斷模型中的識(shí)別率最低,相比方法四與方法五,當(dāng)源域訓(xùn)練中融入變工況數(shù)據(jù)時(shí),遷移效果更好,微調(diào)訓(xùn)練過(guò)的目標(biāo)域樣本集A、B、C的識(shí)別率明顯高于未訓(xùn)練過(guò)的樣本集D。所提方法可有效地實(shí)現(xiàn)采煤機(jī)軸承故障診斷的遷移學(xué)習(xí)。
表3 方法對(duì)比
針對(duì)采煤機(jī)故障診斷過(guò)程中故障樣本數(shù)不足、難以充分訓(xùn)練卷積神經(jīng)網(wǎng)絡(luò)的問(wèn)題,提出變工況下采煤機(jī)故障診斷的遷移學(xué)習(xí)方法,以CWRU數(shù)據(jù)集作為源域,采煤機(jī)故障診斷實(shí)驗(yàn)臺(tái)的故障數(shù)據(jù)作為目標(biāo)域,有效地實(shí)現(xiàn)利用少量樣本微調(diào)權(quán)值以獲得遷移診斷模型。通過(guò)研究得到以下結(jié)論:
(1)以MMD作為優(yōu)化指標(biāo),利用源域數(shù)據(jù)實(shí)現(xiàn)遷移模型的權(quán)值微調(diào),既保留了源域中的模型特征,同時(shí)又融入了目標(biāo)域的權(quán)值特征,由此提高故障診斷識(shí)別效果;
(2)通過(guò)小波變換將一維振動(dòng)信號(hào)變?yōu)闀r(shí)域波形圖,同時(shí)保留了診斷信號(hào)的時(shí)域和頻域特征,并利用圖像增強(qiáng)的方法,強(qiáng)化時(shí)域波形圖中的故障特征,對(duì)比發(fā)現(xiàn):經(jīng)圖像增強(qiáng)后的識(shí)別效果明顯優(yōu)于未經(jīng)增強(qiáng)的;
(3)遷移模型微調(diào)過(guò)程中,僅鎖定低層網(wǎng)絡(luò)的權(quán)值,不對(duì)各層權(quán)值進(jìn)行微調(diào)的診斷效果明顯不如文中所提微調(diào)各層權(quán)值的方法;
(4)源域中融入不同工況數(shù)據(jù)為源域訓(xùn)練集的識(shí)別效果比單獨(dú)用一種工況數(shù)據(jù)訓(xùn)練故障模型的遷移效果好。