基于改進(jìn)GWOLSTM的船舶主機(jī)性能預(yù)測(cè)模型

石彪　王海燕　焦品博

摘要：為提高船舶主機(jī)性能預(yù)測(cè)的精度，更好地幫助輪機(jī)員制定主機(jī)的維修保養(yǎng)計(jì)劃，提出一種利用改進(jìn)灰狼優(yōu)化算法（grey wolf optimizer， GWO）對(duì)長(zhǎng)短期記憶網(wǎng)絡(luò)（long shortterm memory network， LSTM）進(jìn)行優(yōu)化的船舶主機(jī)性能預(yù)測(cè)模型，簡(jiǎn)稱為改進(jìn)的GWOLSTM模型?；贛etropolis接受準(zhǔn)則的思想，在傳統(tǒng)GWO中引入一種隨機(jī)搜索機(jī)制來(lái)解決傳統(tǒng)GWO前期收斂速度慢和后期易陷入局部最優(yōu)的問(wèn)題。分別建立單步預(yù)測(cè)模型和多步預(yù)測(cè)模型。與GWOLSTM和傳統(tǒng)LSTM的預(yù)測(cè)結(jié)果進(jìn)行對(duì)比，結(jié)果表明，改進(jìn)GWOLSTM單步預(yù)測(cè)的均方根誤差分別降低了3236%和5038%，多步預(yù)測(cè)的均方根誤差分別降低了2616%和3557%。

關(guān)鍵詞：? Metropolis接受準(zhǔn)則; 長(zhǎng)短期記憶網(wǎng)絡(luò)（LSTM）; 灰狼優(yōu)化算法（GWO）; 性能預(yù)測(cè)

中圖分類號(hào)：? U6641211文獻(xiàn)標(biāo)志碼：? A

Prediction model of marine main engine performance

based on improved GWOLSTM

Abstract： In order to improve the prediction accuracy of marine main engine performance and help engineers make better maintenance plan of main engines， a prediction model of marine main engine performance is proposed， where the improved grey wolf optimizer （GWO） is used to optimize the long shortterm memory network （LSTM）. The model is called a modified GWOLSTM for short. Based on the idea of Metropolis acceptance criterion， a random search mechanism is introduced to the traditional GWO to solve the problem of slow convergence speed in the early stage and being easy to fall into the local optimum in the later stage of the traditional GWO. The onestep prediction model and the multistep prediction model are established， respectively. Compared with the GWOLSTM and the traditional LSTM， the results show that， the root mean square error of the onestep prediction of the improved GWOLSTM is reduced by 3236% and 5038%， respectively， and the root mean square error of the multistep prediction of the improved GWOLSTM is reduced by 2616% and 3557%， respectively.

Key words： Metropolis acceptance criterion; long shortterm memory network （LSTM）; grey wolf optimizer （GWO）; performance prediction

引言

船舶主機(jī)作為船舶的心臟，長(zhǎng)期處于復(fù)雜多變的工作環(huán)境中，運(yùn)行時(shí)間越長(zhǎng)其工作性能就越差。掌握主機(jī)性能變化趨勢(shì)，有助于輪機(jī)員制訂維修保養(yǎng)計(jì)劃[1]，不僅可以預(yù)防故障的發(fā)生，還能減少停機(jī)、停航時(shí)間，節(jié)約成本。

船舶主機(jī)性能受多種因素影響，對(duì)船舶主機(jī)性能變化的預(yù)測(cè)是一種復(fù)雜且不穩(wěn)定的非線性問(wèn)題。目前，主機(jī)性能趨勢(shì)的預(yù)測(cè)方法由基于模型驅(qū)動(dòng)逐漸轉(zhuǎn)變?yōu)榛跀?shù)據(jù)驅(qū)動(dòng)。基于模型驅(qū)動(dòng)的方法過(guò)于復(fù)雜，且泛化性能較差，很難處理非線性問(wèn)題。在基于數(shù)據(jù)驅(qū)動(dòng)的方法中，循環(huán)神經(jīng)網(wǎng)絡(luò)（recurrent neural network，RNN）因其具有特殊的網(wǎng)絡(luò)結(jié)構(gòu)而適合處理這種非線性問(wèn)題，并且由于RNN還引入了時(shí)序的概念，它能更好地分析主機(jī)性能變化趨勢(shì)。1997年HOCHREITER等[2]對(duì)RNN的結(jié)構(gòu)做了進(jìn)一步的改進(jìn)，提出了長(zhǎng)短期記憶網(wǎng)絡(luò)（long shortterm memory，LSTM），不僅有效地緩解了RNN中梯度爆炸和梯度消失的問(wèn)題，而且解決了RNN的長(zhǎng)時(shí)依賴的問(wèn)題。

近年來(lái)，LSTM在船舶設(shè)備性能、壽命等方面的評(píng)估中得到了廣泛的應(yīng)用。王子銘等[3]使用LSTM建立了船舶主機(jī)油耗實(shí)時(shí)預(yù)測(cè)模型，結(jié)果表明該模型預(yù)測(cè)結(jié)果可靠、誤差較小。易文海等[4]使用LSTM建立波浪與船舶的因果關(guān)系模型，對(duì)隨機(jī)橫浪中船舶的橫搖運(yùn)動(dòng)進(jìn)行短期預(yù)測(cè)，結(jié)果表明此方法可行。QIAO等[5]使用LSTM實(shí)時(shí)預(yù)測(cè)動(dòng)態(tài)系泊線響應(yīng)，結(jié)果表明該方法預(yù)測(cè)精度較高。SHI等[6]利用LSTM對(duì)船舶柴油機(jī)的整體性能做了預(yù)測(cè)，結(jié)果表明該方法可以用于船舶性能的短期預(yù)測(cè)和長(zhǎng)期預(yù)測(cè)。

LSTM在訓(xùn)練數(shù)據(jù)時(shí)普遍使用隨時(shí)間反向傳播（back propogation through time，BPTT）算法，但該算法因復(fù)雜度較高、學(xué)習(xí)效率低而效率低下。灰狼優(yōu)化算法（grey wolf optimizer，GWO）作為啟發(fā)式算法的一種，相比于其他啟發(fā)式算法具有結(jié)構(gòu)簡(jiǎn)單、容易實(shí)現(xiàn)和需要人為調(diào)整的超參數(shù)較少的優(yōu)點(diǎn)[78]，適合用于本文的研究?；诖耍疚奶岢鲇肎WO代替BPTT算法尋優(yōu)的思路，針對(duì)傳統(tǒng)GWO前期收斂速度慢和后期易陷入局部最優(yōu)的缺點(diǎn)，利用Metropolis接受準(zhǔn)則思想提出隨機(jī)搜索機(jī)制，對(duì)GWO的結(jié)構(gòu)進(jìn)行改進(jìn)。使用改進(jìn)后的GWO來(lái)優(yōu)化LSTM，提高船舶主機(jī)性能預(yù)測(cè)精度。

1主機(jī)性能指標(biāo)

1.1監(jiān)測(cè)參數(shù)與數(shù)據(jù)來(lái)源

本文研究數(shù)據(jù)采集于上海海事大學(xué)自動(dòng)化機(jī)艙實(shí)驗(yàn)室，采集對(duì)象是MAN B&W 6S35MEB9型號(hào)的主機(jī)，該主機(jī)的主要技術(shù)參數(shù)見(jiàn)表1。

綜合信息系統(tǒng)通過(guò)傳感器收集監(jiān)測(cè)參數(shù)，參數(shù)經(jīng)過(guò)處理后保存至數(shù)據(jù)庫(kù)。傳感器采樣間隔時(shí)間為1 min。柴油機(jī)的轉(zhuǎn)速和功率信號(hào)通過(guò)磁電式測(cè)試儀采集，溫度信號(hào)由pt100溫度傳感器采集，壓力信號(hào)由壓力變送器采集。

1.2性能參數(shù)選擇

影響船舶主機(jī)性能的熱力參數(shù)眾多，合理選擇熱力參數(shù)是預(yù)測(cè)船舶主機(jī)性能的前提。本文的主機(jī)監(jiān)測(cè)數(shù)據(jù)均是在實(shí)驗(yàn)室固定環(huán)境下采集的，因此環(huán)境溫度、環(huán)境壓力、環(huán)境濕度等外界參數(shù)可以視為固定不變。本文選取的主機(jī)性能參數(shù)見(jiàn)表2。

表2本文選取的主機(jī)性能參數(shù)性能參數(shù)選擇原因功率、轉(zhuǎn)速反映主機(jī)的動(dòng)力性[910]掃氣壓力、掃氣溫度反映主機(jī)的換氣性能冷卻滑油溫度、滑油入口溫度、

本文參數(shù)的選擇基于兩個(gè)方面的原因：一是便于采集，不需要加裝額外的傳感器;二是輪機(jī)員也依靠這些參數(shù)來(lái)判斷主機(jī)的性能狀況。

1.3性能指標(biāo)

采集該主機(jī)運(yùn)行1745 h的性能數(shù)據(jù)作為數(shù)據(jù)源，對(duì)數(shù)據(jù)處理后，選取50%負(fù)荷、轉(zhuǎn)速穩(wěn)定在114 r/min下的數(shù)據(jù)作為觀測(cè)樣本，觀測(cè)樣本共1 009個(gè)，每個(gè)樣本均含有20個(gè)性能數(shù)據(jù)。以主機(jī)出廠試車數(shù)據(jù)中50%負(fù)荷試車工況（見(jiàn)表3）下的100組試車樣本數(shù)據(jù)作為初始性能樣本，計(jì)算觀測(cè)樣本到初始性能樣本的馬氏距離。馬氏距離可以描述兩個(gè)樣本之間的偏離程度。柴油機(jī)的性能參數(shù)具有量綱不同、變化不一和相關(guān)性強(qiáng)的特點(diǎn)，通過(guò)直接觀測(cè)多個(gè)性能參數(shù)判斷柴油機(jī)的工作狀態(tài)容易產(chǎn)生混亂，特別是在某些性能參數(shù)的變化與柴油機(jī)的整體性能變化不一致甚至趨勢(shì)相反時(shí)。馬氏距離的優(yōu)勢(shì)在于：（1）馬氏距離的計(jì)算中包括樣本總體協(xié)方差的計(jì)算，它考慮了樣本的總體分布，消除了各變量之間的相關(guān)性;（2）馬氏距離不受數(shù)據(jù)量綱的影響，即兩點(diǎn)之間的馬氏距離與原始數(shù)據(jù)的測(cè)量單位無(wú)關(guān)。因此，本文用馬氏距離來(lái)表示柴油機(jī)性能的退化程度。馬氏距離越大，說(shuō)明主機(jī)性能退化得越嚴(yán)重。馬氏距離數(shù)學(xué)表達(dá)式如下：（1）式中：DM為觀測(cè)樣本到初始樣本的馬氏距離;x為觀測(cè)樣本矩陣;u為初始樣本的均值矩陣;S-1為初始樣本的協(xié)方差矩陣。

為了較直觀地描述主機(jī)性能的退化過(guò)程，將馬氏距離歸一化成性能指標(biāo)。性能指標(biāo)值越小，說(shuō)明主機(jī)性能退化越嚴(yán)重，其數(shù)學(xué)表達(dá)式如下：（2）式中：α為調(diào)節(jié)因子。α的計(jì)算式為（3）式中：P0為主機(jī)工作初期的性能指標(biāo)值;DM，0為主機(jī)工作初期性能的馬氏距離。

根據(jù)式（1）計(jì)算觀測(cè)樣本到初始樣本的馬氏距離。從圖1可以看出，隨著運(yùn)行時(shí)間的增加，馬氏距離越來(lái)越大。在0號(hào)樣本處，DM為529×105，說(shuō)明此時(shí)主機(jī)的整體性能已經(jīng)遠(yuǎn)遠(yuǎn)偏離初始性能，且處于退化加劇期。柴油機(jī)性能退化過(guò)程主要有兩個(gè)階段：一是緩慢退化期，此階段柴油機(jī)性能接近其出廠性能，馬氏距離較小;二是退化加劇期，此階段柴油機(jī)性能退化速度加劇，馬氏距離較大。

根據(jù)式（3）確定式（2）中的α值，再根據(jù)式（2）將不同時(shí)刻的DM歸一化成P。當(dāng)α=1.2時(shí)，P均勻分布在0與1之間。由圖2可以看出，主機(jī)性能曲線整體上呈下降趨勢(shì)，表示主機(jī)性能逐漸退化。

將前80%的觀測(cè)樣本作為訓(xùn)練樣本進(jìn)行模型訓(xùn)練和參數(shù)調(diào)整，后20%的觀測(cè)樣本作為測(cè)試樣本進(jìn)行模型評(píng)估和預(yù)測(cè)結(jié)果分析。

2預(yù)測(cè)模型

本文首先對(duì)傳統(tǒng)GWO的結(jié)構(gòu)進(jìn)行改進(jìn)，然后在改進(jìn)GWO的基礎(chǔ)上建立改進(jìn)GWOLSTM船舶主機(jī)性能預(yù)測(cè)模型。

2.1LSTM

LSTM與RNN最大的不同在于，LSTM神經(jīng)元中具有3個(gè)“控制門”和2種狀態(tài)，見(jiàn)圖3。在這種結(jié)構(gòu)下，LSTM的細(xì)胞狀態(tài)可以隨時(shí)間進(jìn)行信息傳遞，并通過(guò)“門”增刪信息，使得LSTM具有記憶功能[13]。

首先，遺忘門進(jìn)行信息篩選，剔除掉部分信息，公式如下：（4）式中：ft為t時(shí)刻神經(jīng)元的遺忘門矩陣;Wf為遺忘門權(quán)重矩陣;Xt為t時(shí)刻神經(jīng)元的輸入矩陣;ht-1為t-1時(shí)刻神經(jīng)元的隱藏狀態(tài)輸出矩陣;bf為遺忘門偏置矩陣。

其次，輸入門決定細(xì)胞狀態(tài)的哪些信息被保留，表達(dá)式如下：（5）

（6）

（7）式中：It為t時(shí)刻神經(jīng)元的輸入門矩陣;C～t為t時(shí)刻神經(jīng)元的細(xì)胞候選狀態(tài)矩陣;Ct為t時(shí)刻神經(jīng)元的細(xì)胞狀態(tài)矩陣;WI為輸入門權(quán)重矩陣;bI為輸入門偏置矩陣;Wc為細(xì)胞候選狀態(tài)權(quán)重矩陣;bc為細(xì)胞候選狀態(tài)偏置矩陣;⊙為矩陣的點(diǎn)乘。

最后，輸出門決定當(dāng)前時(shí)刻神經(jīng)元的輸出，表達(dá)式如下：（8）

（9）式中：ot為t時(shí)刻神經(jīng)元的輸出門矩陣;Wo為輸出門權(quán)重矩陣;bo為輸出門偏置矩陣;ht為t時(shí)刻神經(jīng)元的隱藏狀態(tài)輸出矩陣。

2.2GWO

2014年，MIRJALILI等[14]提出一種新的群體智能優(yōu)化算法——GWO。GWO通過(guò)模擬狼群中的等級(jí)制度，將狼群中的狼分為由高到低4個(gè)等級(jí)：Alpha、Beta、Delta和Omega。Alpha狼是狼群中唯一的領(lǐng)導(dǎo)者，領(lǐng)導(dǎo)著種群捕獵。捕獵分為3個(gè)階段：

（1）包圍。在Alpha狼的帶領(lǐng)下，狼群先包圍獵物，數(shù)學(xué)模型如下：（10）

（11）式中：下標(biāo)i為迭代次數(shù)，i=0，1，…，n;D為灰狼與獵物之間的距離;Xp，i為迭代次數(shù)為i時(shí)的獵物位置向量;Xi為迭代次數(shù)為i時(shí)的灰狼位置向量;Xi+1為迭代次數(shù)為i+1時(shí)的灰狼位置向量;A和C為系數(shù)向量。A和C表達(dá)式為（12）

（13）式中：r1和r2為元素取值在區(qū)間[0，1]內(nèi)的隨機(jī)向量;a為控制參數(shù)，隨迭代次數(shù)的增加從2線性減少至0。

（2）獵捕。狼群包圍獵物后，Alpha狼、Beta狼、Delta狼離獵物的距離最近，狼群將在Alpha狼、Beta狼和Delta狼的帶領(lǐng)下向獵物逼近，Omega狼的位置更新如下：（14）

（15）

（16）

（17）

（18）

（19）

（20）式中：Dα、Dβ、Dδ分別為當(dāng)前Omega狼與Alpha狼、Beta狼、Delta狼之間的距離;Xα、Xβ、Xδ分別為Alpha狼、Beta狼、Delta狼的位置向量。

（3）進(jìn)攻。該階段的目標(biāo)是捕獲獵物，即獲得最優(yōu)解。GWO通過(guò)減小a的值模擬逼近獵物的過(guò)程。當(dāng)a逐漸減小時(shí)，A的元素取值也在[-a，a]區(qū)間內(nèi)。當(dāng)A<1時(shí)，狼群能夠攻擊到獵物;當(dāng)A>1時(shí)狼群會(huì)向四周散去，導(dǎo)致狼群失去最優(yōu)位置，這也是GWO容易陷入局部最優(yōu)的原因[15]。

2.3改進(jìn)GWO

GWO的缺點(diǎn)是前期收斂速度慢和后期易陷入局部最優(yōu)[16]。本文借鑒模擬退火的Metropolis接受準(zhǔn)則思想，對(duì)GWO的不足進(jìn)行改進(jìn)。Metropolis接受準(zhǔn)則是指以一定的概率接受新狀態(tài)，從而使算法有能力跳出局部極值和避免過(guò)早收斂。接受新狀態(tài)的概率為

（21）

式中：p（ON）為從狀態(tài)O轉(zhuǎn)移為狀態(tài)N的概率;f（N）和f（O）分別為狀態(tài)O和狀態(tài)N的目標(biāo)函數(shù)值;c為控制參數(shù)。

基于Metropolis接受準(zhǔn)則思想，本文提出隨機(jī)搜索機(jī)制，即在灰狼種群圍捕獵物的過(guò)程中，賦予每只灰狼在自身一定范圍內(nèi)的隨機(jī)搜索能力，當(dāng)灰狼到達(dá)新位置時(shí)，對(duì)灰狼的新位置進(jìn)行適應(yīng)度值計(jì)算，然后與舊位置的適應(yīng)度值進(jìn)行比較。若新位置的適應(yīng)度值小于舊位置的適應(yīng)度值，灰狼將留在新位置;否則，灰狼返回舊位置。數(shù)學(xué)模型如下：（22）

（23）式中：f（Xi）為灰狼原位置的適應(yīng)度值;f（Xi，new）為灰狼新位置的適應(yīng)度值;vmax和vmin分別為解空間的上限和下限;r為隨機(jī)向量，其元素取值在[-1，1]區(qū)間內(nèi);w為隨機(jī)搜索因子。w由下式進(jìn)行更新：

（24）

式中：inow和imax分別為當(dāng)前迭代次數(shù)和最大迭代次數(shù);wmax為最大搜索因子，設(shè)置為0.01;wmin為最小搜索因子，設(shè)置為0.001。

圖4搜索因子隨迭代次數(shù)的變化w的值隨迭代次數(shù)的增加非線性遞減，見(jiàn)圖4。當(dāng)w的值較大時(shí)，灰狼隨機(jī)搜索的范圍較大，可以更好地尋找全局最優(yōu)解;當(dāng)w的值較小時(shí)，灰狼進(jìn)行隨機(jī)搜索的范圍變小，能精確地找到局部最優(yōu)解。

2.4改進(jìn)GWOLSTM船舶主機(jī)性能預(yù)測(cè)模型

2.4.1模型優(yōu)化

使用改進(jìn)GWO來(lái)優(yōu)化LSTM，其優(yōu)化參數(shù)為L(zhǎng)STM中的權(quán)重和偏置。使用向量J來(lái)表示LSTM所有的權(quán)重和偏置參數(shù)，公式如下：（25）式中：J=（J1，J2，…，JK）;fre（·）為可以將矩陣轉(zhuǎn)換成列向量的函數(shù);W為L(zhǎng)STM的各權(quán)重矩陣;b為L(zhǎng)STM的各偏置矩陣。

優(yōu)化目標(biāo)是使測(cè)試樣本的預(yù)測(cè)值與實(shí)際值之間的差值最小。目標(biāo)函數(shù)如下：min f（Y′，Y）

（26）式中：Y′為預(yù)測(cè)值;Y為實(shí)際值。

本文使用均方根誤差（root mean square error，RMSE）函數(shù)作為3種智能算法優(yōu)化的目標(biāo)函數(shù)，數(shù)學(xué)表達(dá)式為f（Y′，Y）=1nni=1（Y′i-Yi）2（27）式中：n為樣本數(shù)。

2.4.2模型預(yù)測(cè)流程

模型預(yù)測(cè)流程見(jiàn)圖5，具體步驟如下：

步驟1根據(jù)權(quán)重和偏置取值范圍隨機(jī)初始化灰狼的位置向量，數(shù)學(xué)表達(dá)式為X0=（r（1），r（2），…，r（K））（28）式中：X0為當(dāng)前灰狼的初始位置向量;r（1），r（2），…，r（K）均為0與1之間的隨機(jī)數(shù)。

步驟2將灰狼的初始位置向量賦值給LSTM中的權(quán)重和偏置（數(shù)學(xué)表達(dá)式如式（28）），利用訓(xùn)練樣本對(duì)模型進(jìn)行訓(xùn)練，以RMSE為適應(yīng)度函數(shù)，將RMSE值中最低的3只灰狼分別作為Alpha狼、Beta狼和Delta狼。（29）步驟3更新灰狼的位置，并計(jì)算各灰狼的適應(yīng)度值，選出新的Alpha狼、Beta狼和Delta狼。

步驟4判斷是否達(dá)到最大迭代次數(shù)，若未達(dá)到則返回步驟3，否則將Alpha狼的位置向量賦值給LSTM中的權(quán)重和偏置，數(shù)學(xué)表達(dá)式為（30）式中：Xα，max為經(jīng)歷了最大次數(shù)的迭代后Alpha狼的位置向量。

步驟5輸入測(cè)試樣本進(jìn)行預(yù)測(cè)。

3實(shí)驗(yàn)與結(jié)果

基于.NET 45.2平臺(tái)使用C#編程語(yǔ)言分別建立LSTM、GWOLSTM、改進(jìn)GWOLSTM的單步預(yù)測(cè)模型和多步預(yù)測(cè)模型，并利用MATLAB繪制預(yù)測(cè)結(jié)果。單步預(yù)測(cè)是指以固定長(zhǎng)度的歷史序列預(yù)測(cè)下一步的序列值;多步預(yù)測(cè)是指根據(jù)固定長(zhǎng)度的歷史序列預(yù)測(cè)未來(lái)連續(xù)多個(gè)時(shí)間步的值。使用第1.3節(jié)的樣本分別訓(xùn)練預(yù)測(cè)3次，然后對(duì)預(yù)測(cè)結(jié)果進(jìn)行對(duì)比分析。

3.1模型評(píng)價(jià)標(biāo)準(zhǔn)

通過(guò)輸入測(cè)試樣本得到預(yù)測(cè)結(jié)果的RMSE和平均絕對(duì)值誤差（mean absolute error，MAE）作為模型的評(píng)價(jià)指標(biāo)。預(yù)測(cè)結(jié)果的RMSE和MAE越低，表示模型預(yù)測(cè)值與實(shí)際值的偏差越小，預(yù)測(cè)精度越高，模型的預(yù)測(cè)效果越好。MAE的數(shù)學(xué)表達(dá)式為（31）

3.2模型參數(shù)設(shè)置

該改進(jìn)GWOLSTM模型由輸入層、1層LSTM層和輸出層組成，輸入維數(shù)為1，輸出維數(shù)為1，灰狼種群數(shù)為30，迭代次數(shù)為500。隱藏層節(jié)點(diǎn)數(shù)一般依靠經(jīng)驗(yàn)設(shè)定，先根據(jù)輸入維數(shù)和輸出維數(shù)確定隱藏層節(jié)點(diǎn)數(shù)范圍，然后根據(jù)實(shí)驗(yàn)選定單步預(yù)測(cè)模型的隱藏層節(jié)點(diǎn)數(shù)為5，輸入層節(jié)點(diǎn)數(shù)為5，輸出層節(jié)點(diǎn)數(shù)為1。

3.3單步預(yù)測(cè)

改進(jìn)GWOLSTM模型、GWOLSTM模型和LSTM模型的預(yù)測(cè)結(jié)果見(jiàn)圖6。各模型評(píng)價(jià)指標(biāo)對(duì)比見(jiàn)表4和5。由圖6、表4和表5可以看出，在3組實(shí)驗(yàn)中，改進(jìn)GWOLSTM模型的預(yù)測(cè)值曲線與真實(shí)值曲線最為吻合，其RMSE和MAE均比其他模型的低。改進(jìn)GWOLSTM模型的RMSE值比LSTM模型、GWOLSTM模型的分別低了5038%、3236%。這說(shuō)明改進(jìn)GWOLSTM模型在船舶主機(jī)性能單步預(yù)測(cè)方面，預(yù)測(cè)效果最好，預(yù)測(cè)精度最高。

3.4多步預(yù)測(cè)

多步預(yù)測(cè)主要有兩種策略：遞歸多步策略（又叫迭代多步策略）和直接多步策略[17]。遞歸多步策略先訓(xùn)練出一個(gè)單步模型，使用該模型預(yù)測(cè)下一個(gè)時(shí)間步的值，然后將該預(yù)測(cè)值作為輸入預(yù)測(cè)下一個(gè)時(shí)間步的值，逐漸循環(huán)，直到預(yù)測(cè)到第b步為止。該方法使用預(yù)測(cè)值代替真實(shí)值，因此該方法會(huì)累積預(yù)測(cè)誤差，導(dǎo)致預(yù)測(cè)結(jié)果的偏差較大，且該模型的性能會(huì)隨著時(shí)間的增加而急劇下降。直接多步策略是直接訓(xùn)練出一個(gè)b步預(yù)測(cè)模型，該模型使用a-h至a時(shí)間步的時(shí)間序列作為輸入，預(yù)測(cè)a+1至a+b時(shí)間步的值。因此本文采取直接多步策略進(jìn)行預(yù)測(cè)。先訓(xùn)練一個(gè)五步預(yù)測(cè)模型，每個(gè)時(shí)間步的長(zhǎng)度為1 min;該模型的隱藏層節(jié)點(diǎn)數(shù)為10，輸入層節(jié)點(diǎn)數(shù)為5，輸出層節(jié)點(diǎn)數(shù)為5。

根據(jù)預(yù)測(cè)結(jié)果，計(jì)算這3種模型的評(píng)價(jià)指標(biāo)，結(jié)果見(jiàn)表6和7。從表6、表7和圖7（各模型多步預(yù)測(cè)結(jié)果比較）可以看出，3組數(shù)據(jù)中，改進(jìn)GWOLSTM模型的預(yù)測(cè)值與實(shí)際值最為貼近，其RMSE和MAE也比其他模型的低。改進(jìn)GWOLSTM模型預(yù)測(cè)結(jié)果的RMSE值比LSTM模型、GWOLSTM模型的分別低了3557%、2616%。這說(shuō)明改進(jìn)GWOLSTM模型在船舶主機(jī)性能多步預(yù)測(cè)方面，仍具有最好的預(yù)測(cè)精度和預(yù)測(cè)效果。

4結(jié)論

為提高對(duì)船舶主機(jī)性能的預(yù)測(cè)精度，在傳統(tǒng)灰狼優(yōu)化算法（GWO）中增加一種隨機(jī)搜索機(jī)制，用改進(jìn)的GWO代替隨時(shí)間反向傳播（BPTT）算法優(yōu)化長(zhǎng)短期記憶網(wǎng)絡(luò)（LSTM）。分別用改進(jìn)GWOLSTM模型、GWOLSTM模型和LSTM模型對(duì)船舶主機(jī)性

能進(jìn)行預(yù)測(cè)，結(jié)果表明：改進(jìn)GWOLSTM模型相對(duì)于GWOLSTM模型和LSTM模型，單步預(yù)測(cè)的RMSE值分別低了3236%和5038%，MAE值分別低了3513%和522%;多步預(yù)測(cè)的RMSE值分別低了2616%和3557%，MAE值分別降低了3049%和4301%。參考文獻(xiàn)：

[1]劉伊凡. 營(yíng)運(yùn)船舶主機(jī)性能退化和健康管理方法研究[D]. 大連： 大連海事大學(xué)， 2018.

[2]HOCHREITER S， SCHMIDHUBER J. Long shortterm memory[J]. Neural Computation， 1997， 9（8）： 17351780. DOI： 101162/neco.19979.81735.

[3]王子銘， 陳順懷. 基于LSTM神經(jīng)網(wǎng)絡(luò)的船舶主機(jī)油耗實(shí)時(shí)預(yù)報(bào)[J]. 武漢理工大學(xué)學(xué)報(bào)（交通科學(xué)與工程版）， 2020， 44（5）： 923927. DOI： 103963/j.issn.20953844202005030.

[4]易文海， 高志亮. 基于LSTM神經(jīng)網(wǎng)絡(luò)的隨機(jī)橫浪中船舶橫搖運(yùn)動(dòng)極短期預(yù)報(bào)[J/OL]. 武漢理工大學(xué)學(xué)報(bào)（交通科學(xué)與工程版）， 2021， 45（6）： 11131117. DOI： 103963/j.issn.20953844202106020.

[5]QIAO D S， LI P， MA G， et al. Realtime prediction of dynamic mooring lines responses with LSTM neural network model[J]. Ocean Engineering， 2021， 219： 108368. DOI： 101016/j.oceaneng.2020108368.

[6]SHI B J， SHI H J， WANG H Y. Performance prediction of marine diesel engine based on long shortterm memory network[C]∥Proceedings of the 2nd International Conference on Artificial Intelligence and Computer Science. Journal of Physics： Conference Series， 2020： 1631（1）： 012135.

[7]HU P， PAN J S， CHU S C. Improved binary grey wolf optimizer and its application for feature selection[J]. KnowledgeBased Systems， 2020， 195： 105746. DOI： 101016/j.knosys.2020105746.

[8]張新明， 姜云， 劉尚旺， 等. 灰狼與郊狼混合優(yōu)化算法及其聚類優(yōu)化[J/OL]. （20200326）[20210518]. 自動(dòng)化學(xué)報(bào)： 117. DOI： 1016383/j.aas.c190617.

[9]劉世元， 杜潤(rùn)生， 楊叔子. 利用轉(zhuǎn)速波動(dòng)信號(hào)診斷內(nèi)燃機(jī)失火故障的研究（1）： 診斷模型方法[J]. 內(nèi)燃機(jī)學(xué)報(bào)， 2000， 18（3）： 315319. DOI： 1016236/j.cnki.nrjxb.200003020.

[10]孫秀成， 譚親明， 曾向明， 等. 基于油耗率的船用低速柴油機(jī)燃燒故障診斷[J]. 中國(guó)航海， 2018， 41（4）： 2327.

[11]張博， 張萍， 郭旭， 等. 基于試驗(yàn)設(shè)計(jì)遺傳算法的船用柴油機(jī)冷卻系統(tǒng)多目標(biāo)優(yōu)化[J]. 推進(jìn)技術(shù)， 2020， 41（11）： 25182529. DOI： 1013675/j.cnki.tjjs.200317.

[12]夏洪永. 船舶主機(jī)冷卻水出口溫度控制系統(tǒng)[J]. 艦船科學(xué)技術(shù)， 2018， 40（12）： 121123. DOI： 103404/j.issn.1672764920186A.041.

[13]LIU J Q， PAN C L， LEI F， et al. Fault prediction of bearings based on LSTM and statistical process analysis[J]. Reliability Engineering & System Safety， 2021， 214： 107646. DOI： 101016/j.ress.2021107646.

[14]MIRJALILI S， MIRJALILI S M， LEWIS A. Grey wolf optimizer[J]. Advances in Engineering Software， 2014， 69： 4661. DOI： 101016/j.advengsoft.201312007.

[15]張曉鳳， 王秀英. 灰狼優(yōu)化算法研究綜述[J]. 計(jì)算機(jī)科學(xué)， 2019， 46（3）： 3038. DOI： 1011896/j.issn.1002137X.201903004.

[16]龍文， 蔡紹洪， 焦建軍， 等. 一種改進(jìn)的灰狼優(yōu)化算法[J]. 電子學(xué)報(bào)， 2019， 47（1）： 169175. DOI： 103969/j.issn.03722112201901022.

[17]TAIEB S B， ATIYA A F. A bias and variance analysis for multistepahead time series forecasting[J]. IEEE Transactions on Neural Networks and Learning Systems， 2016， 27（1）： 6276. DOI： 101109/TNNLS.20152411629.

（編輯賈裙平）

收稿日期： 20210511修回日期： 20210926

基金項(xiàng)目： 上海市科技計(jì)劃（20DZ2252300）

作者簡(jiǎn)介： 石彪（1997—），男，山東菏澤人，碩士研究生，研究方向?yàn)檩啓C(jī)自動(dòng)化、船舶主機(jī)性能預(yù)測(cè)，（Email）18800205042@163.com;

王海燕（1976—），男，河北平山人，副教授，博士，研究方向?yàn)檩啓C(jī)自動(dòng)化、船舶動(dòng)力裝置仿真，（Email）wanghaiyan@shmtu.edu.cn