高學(xué)金, 劉爽爽, 高慧慧

(1. 北京工業(yè)大學(xué) 信息學(xué)部, 北京 100124; 2. 數(shù)字社區(qū)教育部工程研究中心, 北京 100124;3. 城市軌道交通北京實(shí)驗(yàn)室, 北京 100124; 4. 計(jì)算智能與智能系統(tǒng)北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100124)

1 前 言

發(fā)酵過(guò)程是現(xiàn)代流程工業(yè)中常見(jiàn)的一種生產(chǎn)方式[1]，通常表現(xiàn)出動(dòng)態(tài)性、非線性等固有特性，難以建立模型在線監(jiān)測(cè)。因此，亟需建立合適的監(jiān)測(cè)模型實(shí)現(xiàn)過(guò)程安全穩(wěn)定運(yùn)行?；跀?shù)據(jù)驅(qū)動(dòng)的多元統(tǒng)計(jì)方法在工業(yè)過(guò)程領(lǐng)域取得廣泛應(yīng)用，典型方法有主成分分析(principal component analysis，PCA)[2-3]、獨(dú)立分量分析(independent component analysis， ICA)[4-5]、時(shí)間片典型相關(guān)分析(canonical correlation analysis，CCA)[6]，但傳統(tǒng)多元統(tǒng)計(jì)方法模型較為簡(jiǎn)單，不能自適應(yīng)提取數(shù)據(jù)深層特征，在當(dāng)今的大數(shù)據(jù)背景下，面對(duì)數(shù)據(jù)量大、數(shù)據(jù)復(fù)雜多變等問(wèn)題，監(jiān)測(cè)能力遠(yuǎn)遠(yuǎn)不足。

深度學(xué)習(xí)[7-9]在大數(shù)據(jù)特征提取上具有其他算法不可比擬的優(yōu)勢(shì)，通過(guò)逐層的非線性特征提取能深入挖掘數(shù)據(jù)內(nèi)在信息，表現(xiàn)出強(qiáng)大的信息學(xué)習(xí)能力和表達(dá)能力，并且對(duì)時(shí)序相關(guān)的序列數(shù)據(jù)同樣具有較強(qiáng)的適用性。近年來(lái)深度學(xué)習(xí)方法因優(yōu)越性能被廣泛用于工業(yè)過(guò)程，如深層置信網(wǎng)絡(luò)[10]、循環(huán)神經(jīng)網(wǎng)絡(luò)[11]、堆疊去噪自編碼[12]和卷積神經(jīng)網(wǎng)絡(luò)(convolution neural network，CNN)[13]，尤其是CNN 優(yōu)越的特征提取能力越來(lái)越受到科研人員的青睞。Lee 等[14]提出故障檢測(cè)和分類(lèi)卷積神經(jīng)網(wǎng)絡(luò)方法，應(yīng)用于半導(dǎo)體制造過(guò)程。Zhang 等[15]提出一種第一層寬卷積核的深度卷積神經(jīng)網(wǎng)絡(luò)方法，解決CNN 應(yīng)用于故障診斷準(zhǔn)確率不高的問(wèn)題。Jiang 等[16]針對(duì)齒輪箱振動(dòng)信號(hào)固有的多尺度特性，提出一種多尺度卷積神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)，實(shí)現(xiàn)齒輪箱振動(dòng)信號(hào)的多尺度特征提取和分類(lèi)。但多數(shù)學(xué)者直接處理原數(shù)據(jù)訓(xùn)練模型用于軸承、齒輪等故障診斷。

過(guò)程擾動(dòng)實(shí)際上發(fā)生在不同的時(shí)間尺度上，即測(cè)量過(guò)程數(shù)據(jù)具有時(shí)序相關(guān)性[17]。數(shù)據(jù)的動(dòng)態(tài)特性對(duì)于統(tǒng)計(jì)量的統(tǒng)計(jì)特性影響很大，進(jìn)而影響過(guò)程的監(jiān)測(cè)結(jié)果[18]。Zhang 等[19-21]研究討論了此問(wèn)題，但研究中仍然假設(shè)不同時(shí)間段提取的潛在變量特征為獨(dú)立狀態(tài)，與變量先前值和其他變量之間不相關(guān)。Zhao 等[22-23]提出慢特性分析挖掘批處理過(guò)程時(shí)變動(dòng)態(tài)信息，依賴(lài)于復(fù)雜的核計(jì)算，并且不能自動(dòng)學(xué)習(xí)內(nèi)在時(shí)變特征。

鑒此，本文提出一種多變量深度卷積神經(jīng)網(wǎng)絡(luò)(multivariable deep convolution neural network，MDCNN)故障監(jiān)測(cè)方法。利用深度卷積神經(jīng)網(wǎng)絡(luò)提取分解變換后數(shù)據(jù)的高層次特征。該方法對(duì)具有動(dòng)態(tài)性的多個(gè)過(guò)程變量進(jìn)行時(shí)頻分析，更完整地表征數(shù)據(jù)特征，使模型監(jiān)測(cè)準(zhǔn)確性大大提高。

2 LeNet-5 卷積神經(jīng)網(wǎng)絡(luò)

卷積神經(jīng)網(wǎng)絡(luò)最早于20 世紀(jì) 80 年代末提出，隨后幾年卷積神經(jīng)網(wǎng)絡(luò)性能持續(xù)提升[24]。LeNet-5[25]是典型的 CNN 模型，包含輸入層，2 個(gè)交替的卷積和池化層，2 個(gè)全連接層，用于分類(lèi)的輸出層。LeNet-5 卷積網(wǎng)絡(luò)模型如圖1 所示。

圖1 卷積神經(jīng)網(wǎng)絡(luò)模型Fig.1 Schematic diagram of the convolutional neural network model

3 基于MDCNN 的故障監(jiān)測(cè)

3.1 數(shù)據(jù)預(yù)處理

發(fā)酵過(guò)程數(shù)據(jù)為三維數(shù)據(jù)結(jié)構(gòu)，還有變量間和時(shí)刻間耦合相關(guān)特性，需將三維數(shù)據(jù)X(I×J×K)按變量方向展開(kāi)為二維矩陣X(KI×J)，其中，I、J、K 分別為批次數(shù)、過(guò)程變量、采樣時(shí)刻。

滑動(dòng)窗在很多領(lǐng)域用來(lái)分析動(dòng)態(tài)信號(hào)。采用滑窗對(duì)測(cè)量的所有變量在采樣時(shí)間上進(jìn)行滑動(dòng)截取，如圖2所示。對(duì)一個(gè)批次，令滑動(dòng)窗的大小為k×j，k 為截取的采樣時(shí)刻數(shù)，j 為所有測(cè)量變量數(shù)，滑動(dòng)步長(zhǎng)為1，可以獲得所有變量每個(gè)采樣時(shí)刻的時(shí)序相關(guān)性，直至獲得所有批次數(shù)據(jù)。

圖2 按變量展開(kāi)數(shù)據(jù)滑窗獲取Fig.2 Schematic diagram of expand data window acquisition by variables

3.2 時(shí)頻圖獲取

當(dāng)工業(yè)過(guò)程發(fā)生故障時(shí)，信號(hào)的幅值、相位和頻率都會(huì)發(fā)生變化，故障監(jiān)測(cè)實(shí)質(zhì)上也是信號(hào)分析處理的過(guò)程。希爾伯特-黃變換(Hilbert-Huang transform，HHT)[26]是性能較好的時(shí)頻分析方法。它包含2部分：經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition，EMD)和希爾伯特變換。HHT 方法在機(jī)械等方面對(duì)單變量進(jìn)行分析，本文針對(duì)發(fā)酵過(guò)程多批次的多個(gè)變量進(jìn)行時(shí)頻分析，具體步驟如下：

1) 首先利用 EMD 方法對(duì)滑窗獲得數(shù)據(jù)中每個(gè)變量分解為有限個(gè)本征模態(tài)函數(shù)(intrinisic mode function，IMF)，所分解出來(lái)的各IMF 分量包含原信號(hào)的一段時(shí)間尺度的局部特征信號(hào)；公式為

如圖3 所示為EMD 分解圖和多個(gè)變量m 個(gè)時(shí)刻的時(shí)頻圖。

圖3 EMD 分解圖和多變量 m 個(gè)時(shí)刻的時(shí)頻圖Fig.3 EMD decomposition diagram and time-frequency diagram of multivariable at m-times

3.3 MDCNN 網(wǎng)絡(luò)訓(xùn)練

針對(duì)發(fā)酵過(guò)程特有的動(dòng)態(tài)性和非線性，本文提出MDCNN 的網(wǎng)絡(luò)如圖4 所示。

圖4 MDCNN 網(wǎng)絡(luò)Fig.4 Schematic diagram of mDCNN network

CNN 輸入層每個(gè)輸入分布并不一致，批量歸一化(BN)[27]層旨在減少內(nèi)部協(xié)方差的轉(zhuǎn)移，加速深層神經(jīng)網(wǎng)絡(luò)的訓(xùn)練過(guò)程。在本節(jié)中，基于LeNet-5 模型加入BN 層來(lái)提高模型訓(xùn)練精度，在第2 層卷積層之后加入BN 層。批量歸一化操作為

在大型網(wǎng)絡(luò)訓(xùn)練中容易產(chǎn)生過(guò)擬合，dropout 通過(guò)在訓(xùn)練階段的迭代中省略一些節(jié)點(diǎn)來(lái)避免過(guò)擬合現(xiàn)象，在實(shí)驗(yàn)中設(shè)置dropout 率ρ=0.5 在全連接層使用，卷積層沒(méi)有很多參數(shù)，不設(shè)置dropout。

MDCNN 監(jiān)測(cè)方法的結(jié)構(gòu)決定著特征提取的完整與否和模型訓(xùn)練時(shí)間；參數(shù)、深度和神經(jīng)元個(gè)數(shù)直接影響實(shí)驗(yàn)結(jié)果的精度。實(shí)驗(yàn)多種MDCNN 網(wǎng)絡(luò)，其模型結(jié)構(gòu)如表1 所示。網(wǎng)絡(luò)的深度實(shí)驗(yàn)在表1中從1 層卷積到4 層卷積層，更深層次的網(wǎng)絡(luò)監(jiān)測(cè)時(shí)間會(huì)大大提高，耗時(shí)較多，但監(jiān)測(cè)準(zhǔn)確性提高不明顯。所以在其合適的范圍內(nèi)設(shè)計(jì)實(shí)驗(yàn)，表1 中列出7 個(gè)模型進(jìn)行對(duì)比分析，具有最好故障監(jiān)測(cè)性能的網(wǎng)絡(luò)作為實(shí)驗(yàn)?zāi)Ｐ?。?duì)于網(wǎng)絡(luò)參數(shù)的確定，以發(fā)酵過(guò)程的數(shù)據(jù)特性和訓(xùn)練模型的精確度與耗時(shí)長(zhǎng)短為依據(jù)研究網(wǎng)絡(luò)的參數(shù)大小選擇。參數(shù)主要包括每個(gè)卷積層的濾波器數(shù)量和FC 層的輸出長(zhǎng)度。

一個(gè)樣本為采用滑窗截取大小為 32×10 的二維輸入數(shù)據(jù)序列，例如模型5，包含2 個(gè)卷積層(C)，2 個(gè)池化層(P)，1 個(gè)全連接層(FC)和 1 個(gè)輸出層(O)。在 2 個(gè)卷積層中，內(nèi)核大小為5×5，步幅為 1。第 1 個(gè)卷積層包含64 個(gè)濾波器，第2 個(gè)卷積層包含32 個(gè)濾波器；BN 層加在第2 個(gè)卷積層之后進(jìn)行映射；2 個(gè)最大池化層位于第1 和第2 卷積層之后，2×2 為在2 個(gè)池化層中的內(nèi)核大??；FC 層的輸入是一維向量，大小為(20×8×32)；20 為FC 層的輸出行，“dropout”用于該層。輸出層具有“softmax”函數(shù)的樣本矩陣類(lèi)，模型的輸出長(zhǎng)度為2 的向量，表示監(jiān)測(cè)結(jié)果的二分類(lèi)。

表1 用于發(fā)酵過(guò)程故障監(jiān)測(cè)的MDCNN 模型結(jié)構(gòu)設(shè)計(jì)Table 1 Design of MDCNN network structure for fault monitoring of fermentation processes

3.4 基于MDCNN 的發(fā)酵過(guò)程故障監(jiān)測(cè)

發(fā)酵過(guò)程具有周期性，在監(jiān)測(cè)過(guò)程中，樣本分為正常批次和故障批次，在建立模型時(shí)分為離線和在線階段。訓(xùn)練階段樣本既包含正常批次也包含故障批次，經(jīng)過(guò)數(shù)據(jù)預(yù)處理變換之后的數(shù)據(jù)作用于監(jiān)測(cè)過(guò)程樣本?；贛DCNN 故障監(jiān)測(cè)方法如圖5 所示。

圖5 基于MDCNN的故障診斷方法框架Fig.5 Fault diagnosis framework based on MDCNN

離線階段：

步驟1：從發(fā)酵過(guò)程數(shù)據(jù)庫(kù)中獲取歷史數(shù)據(jù)；

步驟2：將三維數(shù)據(jù)變量方法展開(kāi)為二維數(shù)據(jù)，采用滑動(dòng)窗采集m×n 的數(shù)據(jù)矩陣，進(jìn)行時(shí)頻分析，獲得所需樣本，并用相應(yīng)的類(lèi)別標(biāo)記，包括“正常”及“故障”類(lèi)別標(biāo)簽；

步驟3：將包括其相應(yīng)標(biāo)簽的樣本分為訓(xùn)練集和測(cè)試集；

步驟4：設(shè)計(jì)網(wǎng)絡(luò)模型，輸入訓(xùn)練樣本進(jìn)行模型訓(xùn)練；

步驟5：輸入測(cè)試樣本進(jìn)行模型測(cè)試；

步驟 6：如果測(cè)試結(jié)果精確度較高，輸出測(cè)試結(jié)果，用于在線故障監(jiān)測(cè)；否則，重新設(shè)計(jì)網(wǎng)絡(luò)進(jìn)行訓(xùn)練。

在線監(jiān)測(cè)：

步驟1：從發(fā)酵過(guò)程獲取新批次的m 時(shí)刻數(shù)據(jù)；

步驟2：對(duì)獲得的變量數(shù)據(jù)采用離線模型數(shù)據(jù)處理方式獲得所需樣本；

步驟3：將在線獲得樣本輸入模型測(cè)試其監(jiān)測(cè)結(jié)果，分為有故障和無(wú)故障類(lèi)型；

步驟4：監(jiān)測(cè)結(jié)束之后，若測(cè)試結(jié)果準(zhǔn)確度較高，在線數(shù)據(jù)可作為歷史數(shù)據(jù)訓(xùn)練模型。

4 實(shí)驗(yàn)驗(yàn)證

4.1 青霉素發(fā)酵過(guò)程仿真驗(yàn)證

青霉素是人類(lèi)大規(guī)模用于臨床的一種抗生素[28]，其發(fā)酵過(guò)程是一種典型的生化反應(yīng)過(guò)程。復(fù)雜工業(yè)過(guò)程中一個(gè)變量的波動(dòng)可能會(huì)引起其他多個(gè)變量的波動(dòng)，采用單變量監(jiān)控會(huì)導(dǎo)致各監(jiān)測(cè)結(jié)果之間難以綜合考慮，影響監(jiān)測(cè)的質(zhì)量。所以實(shí)驗(yàn)選取發(fā)酵過(guò)程的 10 個(gè)變量(x1，x2，……x10)進(jìn)行監(jiān)測(cè)，實(shí)現(xiàn)對(duì)過(guò)程運(yùn)行在線監(jiān)測(cè)。所選變量如表2 所示。

青霉素發(fā)酵過(guò)程每個(gè)批次的反應(yīng)時(shí)間為400 h，采樣間隔為0.5 h。實(shí)驗(yàn)數(shù)據(jù)集使用Pensim 2.0 平臺(tái)仿真得到的5 批正常工況數(shù)據(jù)和5 批故障數(shù)據(jù)作為訓(xùn)練樣本，3 批不同于訓(xùn)練樣本的故障數(shù)據(jù)作為測(cè)試樣 本 ， 即 三 維 數(shù) 據(jù) X(10×10×800)、Y(3×10×800)，實(shí)驗(yàn)進(jìn)行二分類(lèi)，測(cè)試結(jié)果分為有故障和無(wú)故障。發(fā)酵過(guò)程故障大致分為 3 大類(lèi)：傳感器故障、過(guò)程故障(包括誤操作)、機(jī)械與電氣故障。目前仿真的故障主要集中在過(guò)程故障和機(jī)械故障上。Pensim 可以設(shè)定的故障變量和故障擾動(dòng)類(lèi)型如表3 所示。

表2 青霉素發(fā)酵過(guò)程變量Table 2 Process variables of penicillin fermentation

表3 青霉素實(shí)驗(yàn)故障類(lèi)型Table 3 Fault types of penicillin fermentation

在完成MDCNN 網(wǎng)絡(luò)結(jié)構(gòu)構(gòu)建之后，為了顯示故障監(jiān)測(cè)的實(shí)驗(yàn)結(jié)果，定義實(shí)驗(yàn)的精確度計(jì)算。所有訓(xùn)練數(shù)據(jù)集都訓(xùn)練一次，稱(chēng)為epoch，批次大小設(shè)為256。訓(xùn)練階段，1 個(gè)批次樣本通過(guò)網(wǎng)絡(luò)訓(xùn)練一次，包含一次前向傳播和一次后向傳播，每迭代一次權(quán)重更新一次；測(cè)試階段，1 個(gè)批次樣本通過(guò)網(wǎng)絡(luò)測(cè)試一次包含一次前向傳播，一次epoch 需要完成的批次數(shù)量等于迭代次數(shù)。對(duì)于不同的數(shù)據(jù)集，數(shù)據(jù)的多樣性會(huì)影響合適的 epoch 數(shù)量，隨著epoch 的增加，訓(xùn)練和測(cè)試的樣本集呈現(xiàn)不同精確度，如圖6 所示。

當(dāng)網(wǎng)絡(luò)訓(xùn)練次數(shù) epoch>100 時(shí)，精確度為飽和狀態(tài)，每訓(xùn)練一次耗時(shí)將增加，實(shí)驗(yàn)設(shè)為 epoch=90，訓(xùn)練樣本精確度為98.05%，測(cè)試樣本精確度為 95.2%。在達(dá)到訓(xùn)練的快速性的同時(shí)保證精確度。

表4 列出了表1 中模型的測(cè)試平均精確度和訓(xùn)練時(shí)間。模型6 訓(xùn)練和測(cè)試平均精確度最高(98.05%、95.2%)，訓(xùn)練時(shí)間30 s×90=45 min。模型1 和2 淺層的結(jié)構(gòu)具有較低的故障監(jiān)測(cè)精確度，相對(duì)于模型4，在卷積層之后加入BN 層。雖然訓(xùn)練時(shí)間相對(duì)增加，但測(cè)試精確度卻大幅度提高。模型 5 訓(xùn)練精確度與模型6 相近，但測(cè)試精確度低。綜合比較，選擇模型6 為下面實(shí)驗(yàn)?zāi)Ｐ?。通過(guò)實(shí)驗(yàn)，將模型 6 的參數(shù)設(shè)置為{C(64)、C(64)、C(32)、FC(20)}。

圖6 迭代過(guò)程的訓(xùn)練精度和測(cè)試精度Fig.6 Training accuracy and testing accuracy of the iteration processes

圖7 多個(gè)模型監(jiān)測(cè)結(jié)果Fig.7 Monitoring results of multiple models

表4 模型平均監(jiān)測(cè)精度Table 4 Average monitoring accuracy of different models

在工業(yè)過(guò)程故障診斷研究中，Wu 等[29]提出的基于 DCNN 的模型用于田納西伊斯曼(tennessee eastman, TE)過(guò)程取得了很好的診斷結(jié)果。為了說(shuō)明本文方法的性能，將實(shí)驗(yàn)?zāi)Ｐ秃推渌Ｐ陀糜诎l(fā)酵過(guò)程數(shù)據(jù)的結(jié)果進(jìn)行比較，如圖7 所示。對(duì)于訓(xùn)練數(shù)據(jù)集，基于支持向量機(jī)(support vector machine，SVM)、BP 神經(jīng)網(wǎng)絡(luò)的模型，訓(xùn)練精確度在90%以下?；贒CNN 的模型在訓(xùn)練數(shù)據(jù)集的故障監(jiān)測(cè)表現(xiàn)出了較高準(zhǔn)確度，而本文提出MDCNN 模型精確度高達(dá)98.05%，比DCNN 模型精確度高2.91%，比標(biāo)準(zhǔn)CNN模型精確度高出4.99%。對(duì)于測(cè)試數(shù)據(jù)集，本文方法平均精確度達(dá)到95%以上，表現(xiàn)出優(yōu)異的監(jiān)測(cè)性能。

如圖8 所示為發(fā)酵過(guò)程二分類(lèi)可視化結(jié)果，圖(a)～(d)分別為SVM，CNN，DCNN，MDCNN 監(jiān)測(cè)結(jié)果。圖中框圖表示為期望(標(biāo)簽)測(cè)試集分類(lèi)結(jié)果，星號(hào)表示實(shí)際測(cè)試集分類(lèi)結(jié)果，橫坐標(biāo)為樣本個(gè)數(shù)，如果框圖和星號(hào)重疊度越高表明監(jiān)測(cè)的精確度越高。從圖中可以看出本文方法基本沒(méi)有錯(cuò)分和漏分情況，其測(cè)試精確度高達(dá)95%以上，說(shuō)明將滑動(dòng)窗和時(shí)頻分析綜合考慮之后采用深度卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行模型訓(xùn)練，測(cè)試精確度更高，更貼近實(shí)際生產(chǎn)過(guò)程。

圖8 青霉素發(fā)酵過(guò)程監(jiān)測(cè)可視化結(jié)果Fig.8 Visual monitoring results of penicillin fermentation processes

4.2 大腸桿菌發(fā)酵過(guò)程實(shí)驗(yàn)研究

用于大腸桿菌發(fā)酵過(guò)程故障監(jiān)測(cè)的實(shí)驗(yàn)數(shù)據(jù)來(lái)源于工廠，實(shí)際的生產(chǎn)環(huán)境往往存在諸多不確定因素，對(duì)于算法研究更具有現(xiàn)實(shí)意義[30]。

實(shí)驗(yàn)選用來(lái)源于北京某生物制藥廠的大腸桿菌發(fā)酵數(shù)據(jù)。發(fā)酵罐容量15 L，發(fā)酵時(shí)間為6～7 h，采樣間隔為5 min。由于實(shí)際數(shù)據(jù)批次之間的差異，本實(shí)驗(yàn)采取批次之間的等長(zhǎng)78 時(shí)刻樣本數(shù)據(jù)，選取了7 個(gè)主要過(guò)程變量，變量名如表5 所示。

為驗(yàn)證所研究監(jiān)測(cè)方法的有效性，用于模型訓(xùn)練和測(cè)試的數(shù)據(jù)共15 批，包含7 批正常數(shù)據(jù)和8 批故障數(shù)據(jù)，其中，6 批正常數(shù)據(jù)和6 批故障數(shù)據(jù)用于訓(xùn)練；1 批正常數(shù)據(jù)和2 批故障數(shù)據(jù)用于測(cè)試。測(cè)試所用故障類(lèi)型如表 6 所示，分別模擬大腸桿菌發(fā)酵過(guò)程中可能發(fā)生的攪拌電機(jī)故障和蒸汽管道故障。用于實(shí)驗(yàn)的MDCNN的網(wǎng)絡(luò)結(jié)構(gòu)為模型6。

表7 是作用于實(shí)際大腸桿菌數(shù)據(jù)的模型監(jiān)測(cè)結(jié)果，從結(jié)果中可以看出，在測(cè)試實(shí)際生產(chǎn)數(shù)據(jù)中，與其他模型相比，本文提出的模型取得了較好的結(jié)果。平均監(jiān)測(cè)精度高達(dá)93.42%，優(yōu)于其他所有方法。其他模型只有DCNN 達(dá)到90% 以上，基于SVM 和BP 神經(jīng)網(wǎng)絡(luò)模型監(jiān)測(cè)效果較差，足以證明本文提出方法具有顯著的實(shí)際監(jiān)測(cè)性能。

表5 大腸桿菌發(fā)酵過(guò)程變量Table 5 Variables of E.coli fermentation processes

表6 大腸桿菌實(shí)驗(yàn)故障類(lèi)型Table 6 Fault types of E.coli fermentation

表7 實(shí)驗(yàn)?zāi)Ｐ团c其他模型監(jiān)測(cè)結(jié)果Table 7 Monitoring results of current model and other models

如圖9 所示為大腸桿菌發(fā)酵過(guò)程監(jiān)測(cè)可視化結(jié)果，圖(a)～(d)為SVM，CNN，DCNN，MDCNN 模型監(jiān)測(cè)結(jié)果。從圖中可以看出本文方法分類(lèi)準(zhǔn)確度最高，具有較好的監(jiān)測(cè)結(jié)果。

圖9 大腸桿菌發(fā)酵過(guò)程監(jiān)測(cè)可視化結(jié)果Fig.9 Visual monitoring results of Escherichia coli fermentation processes

5 結(jié) 語(yǔ)

本文提出一種基于多變量深度卷積神經(jīng)網(wǎng)絡(luò)的故障監(jiān)測(cè)方法。該方法針對(duì)具有動(dòng)態(tài)性和非線性的發(fā)酵過(guò)程，利用希爾伯特-黃變換提取多變量動(dòng)態(tài)數(shù)據(jù)的初級(jí)特征，更好地表征了數(shù)據(jù)時(shí)頻域信息；然后輸入卷積神經(jīng)網(wǎng)絡(luò)自動(dòng)學(xué)習(xí)數(shù)據(jù)的高層特征，提高了分類(lèi)準(zhǔn)確度和監(jiān)測(cè)性能。本文在青霉素發(fā)酵過(guò)程和大腸桿菌生產(chǎn)數(shù)據(jù)上分別進(jìn)行了實(shí)驗(yàn)驗(yàn)證。結(jié)果證明本文所提方法具有更高的監(jiān)測(cè)精度，在實(shí)際工業(yè)生產(chǎn)中有較強(qiáng)的應(yīng)用價(jià)值。