郭海新

(昆侖能源湖北黃岡液化天然氣有限公司，湖北黃岡 438000)

在液化天然氣(LNG)儲(chǔ)罐使用的過程中，由于介質(zhì)中含有一定的腐蝕性物質(zhì)及水分，腐蝕性物質(zhì)溶于水后會(huì)導(dǎo)致儲(chǔ)罐腐蝕，進(jìn)而影響儲(chǔ)罐的使用安全，大大縮短儲(chǔ)罐的使用壽命[1-2]。此外，腐蝕產(chǎn)物還會(huì)直接進(jìn)入到LNG之中，對(duì)LNG的質(zhì)量產(chǎn)生嚴(yán)重的影響。如能及時(shí)掌握儲(chǔ)罐的腐蝕狀況，就可以提前采取預(yù)防措施，保障儲(chǔ)罐的運(yùn)行安全[3-4]。

目前，隨著LNG行業(yè)的飛速發(fā)展，對(duì)LNG儲(chǔ)罐安全問題的研究較多。單彤文[5]針對(duì)儲(chǔ)罐運(yùn)行安全，對(duì)其研究進(jìn)展及發(fā)展趨勢(shì)進(jìn)行了研究。結(jié)果表明，為了提高LNG儲(chǔ)罐使用的安全性，需要將理論計(jì)算方法與監(jiān)測(cè)措施相結(jié)合，及時(shí)了解儲(chǔ)罐內(nèi)的腐蝕狀況，同時(shí)，還需要采取較精細(xì)的管理措施。劉旭霞等[6]通過進(jìn)行腐蝕試驗(yàn)的方式對(duì)LNG儲(chǔ)罐試壓過程中的腐蝕問題進(jìn)行合理全面的研究。結(jié)果表明，為了保障儲(chǔ)罐運(yùn)行的安全性，需要采取合理的陰極保護(hù)措施，同時(shí)，需要根據(jù)儲(chǔ)罐材質(zhì)的不同，對(duì)保護(hù)電位進(jìn)行合理的選擇，若保護(hù)電位不合理，不但會(huì)出現(xiàn)嚴(yán)重的腐蝕問題，還可能會(huì)出現(xiàn)氫脆風(fēng)險(xiǎn)。黃剛?cè)A等[7]對(duì)LNG儲(chǔ)罐腐蝕的影響因素進(jìn)行了分析和研究。結(jié)果表明，硫化氫、氧氣以及水是引起儲(chǔ)罐腐蝕問題的主要因素，大量的腐蝕性物質(zhì)會(huì)在儲(chǔ)罐的頂部位置沉積，進(jìn)而使其頂部位置的腐蝕速率加快，為了保障儲(chǔ)罐的運(yùn)行安全，必須采取合理的防腐措施。

通過采用不同試驗(yàn)方式對(duì)儲(chǔ)罐的腐蝕情況進(jìn)行研究，將智能算法引入該領(lǐng)域，可及時(shí)了解LNG儲(chǔ)罐的腐蝕狀況。因此，筆者采用GWO-LSSVM組合算法，并將該類型的組合算法應(yīng)用到儲(chǔ)罐的腐蝕速率預(yù)測(cè)中，對(duì)儲(chǔ)罐的腐蝕狀況進(jìn)行了分析研究。

1 LNG儲(chǔ)罐腐蝕理論分析

在LNG儲(chǔ)罐的介質(zhì)中，含有大量的腐蝕性物質(zhì)以及水分，儲(chǔ)罐長(zhǎng)時(shí)間使用后，出現(xiàn)腐蝕問題的概率相對(duì)較大，若儲(chǔ)罐的腐蝕問題較為嚴(yán)重，不但可能會(huì)出現(xiàn)各種類型的安全風(fēng)險(xiǎn)，還可能會(huì)大大降低儲(chǔ)罐的使用壽命[8]。此外，儲(chǔ)罐使用過程中防腐處理問題較易被忽略，儲(chǔ)罐被腐蝕后，大量的腐蝕產(chǎn)物無法排除罐外，會(huì)在儲(chǔ)罐內(nèi)大量聚集，并與儲(chǔ)存介質(zhì)相互接觸，對(duì)介質(zhì)的質(zhì)量產(chǎn)生嚴(yán)重影響。通過對(duì)LNG儲(chǔ)罐腐蝕問題進(jìn)行調(diào)研后發(fā)現(xiàn)，儲(chǔ)罐內(nèi)部腐蝕問題的類型相對(duì)較多，主要包括磨損腐蝕、點(diǎn)腐蝕及均勻腐蝕等。從化學(xué)角度進(jìn)行分析，石化領(lǐng)域的腐蝕問題主要可以分為2種類型，分別為化學(xué)腐蝕以及電化學(xué)腐蝕，其中，化學(xué)腐蝕主要出現(xiàn)在特定條件下，LNG儲(chǔ)罐內(nèi)的腐蝕問題主要為電化學(xué)腐蝕。在儲(chǔ)罐正式使用前，一般均會(huì)對(duì)其進(jìn)行合理的防腐處理，采取一定的防腐措施，但是受到施工工藝以及儲(chǔ)罐操作等因素的影響，使用一段時(shí)間后，防腐措施較易失效，儲(chǔ)罐出現(xiàn)腐蝕問題的概率升高[9]。

研究者對(duì)LNG儲(chǔ)罐腐蝕問題開展的研究時(shí)間相對(duì)較早，也取得了較多的研究成果。20世紀(jì)50年代，美國(guó)石油化工協(xié)會(huì)就開始組織對(duì)LNG儲(chǔ)罐的腐蝕問題的研究，研究結(jié)果發(fā)現(xiàn)在LNG儲(chǔ)罐內(nèi)有游離態(tài)水存在的前提下，若H2S的分壓相對(duì)較高，超過了0.000 35 MPa，則H2S對(duì)腐蝕的影響就會(huì)加重。LNG儲(chǔ)罐內(nèi)的環(huán)境狀況以及儲(chǔ)罐的材料類型也是影響腐蝕問題的關(guān)鍵因素，研究結(jié)果表明，若儲(chǔ)罐內(nèi)環(huán)境條件不同，其腐蝕速率存在較大差別；若儲(chǔ)罐內(nèi)LNG的質(zhì)量分?jǐn)?shù)相對(duì)較高，則儲(chǔ)罐的使用周期可以達(dá)到20~25年，若儲(chǔ)罐內(nèi)的LNG質(zhì)量分?jǐn)?shù)相對(duì)較低，儲(chǔ)罐使用5~10年后，就會(huì)出現(xiàn)較嚴(yán)重的腐蝕問題。

LNG儲(chǔ)罐的腐蝕主要可以分為3個(gè)方面：分別是罐底的腐蝕、罐壁的腐蝕以及罐頂?shù)母g，其中罐底的腐蝕問題相對(duì)較嚴(yán)重，主要是由于該位置的腐蝕情況較為復(fù)雜，部分腐蝕性物質(zhì)會(huì)產(chǎn)生直接的化學(xué)腐蝕，同時(shí)，腐蝕性產(chǎn)物可能會(huì)與某些酸性成分相互作用，進(jìn)而產(chǎn)生電化學(xué)腐蝕。罐頂處會(huì)聚集大量的水蒸氣以及氧氣，產(chǎn)生電化學(xué)腐蝕。罐底位置主要以濃差腐蝕為主，腐蝕程度相對(duì)較輕[10]。

2 算法理論模型

2.1 灰狼優(yōu)化算法

灰狼優(yōu)化算法(GWO算法)是一種智能優(yōu)化算法，使用該類型的算法可以在全局之內(nèi)尋找最優(yōu)解，其原理主要是通過對(duì)灰狼群的覓食行為進(jìn)行合理的模擬，最終達(dá)到目標(biāo)最優(yōu)化，該種類型算法在使用的過程中尋優(yōu)的精度相對(duì)較高，其收斂的速度相對(duì)較快[11-14]。

2.2 最小二乘支持向量機(jī)算法

最小二乘支持向量機(jī)(LSSVM)算法是一種智能預(yù)測(cè)方法，其屬于支持向量機(jī)(SVM)算法的一種改進(jìn)模式，在使用過程中，可以降低數(shù)據(jù)樣本的訓(xùn)練復(fù)雜性，提升運(yùn)算速度，在進(jìn)行機(jī)器學(xué)習(xí)的過程中具有較強(qiáng)的廣泛性及先進(jìn)性。在使用LSSVM算法的過程中，核函數(shù)的類型相對(duì)較多，其中RBF核函數(shù)的應(yīng)用相對(duì)較廣，因此，筆者選用RBF核函數(shù)，該類型核函數(shù)中含有g(shù)am以及sig2參數(shù)，這2種參數(shù)的設(shè)定會(huì)對(duì)最終的預(yù)測(cè)結(jié)果產(chǎn)生較大影響。

2.3 GWO-LSSVM組合模型建模

該研究建立的GWO-LSSVM組合算法實(shí)質(zhì)上是使用GWO算法對(duì)LSSVM算法中的部分參數(shù)進(jìn)行尋優(yōu)，即通過使用GWO算法確定LSSVM算法中g(shù)am及sig2參數(shù)的數(shù)值，利用優(yōu)化后的模型，對(duì)LNG儲(chǔ)罐罐底、罐壁以及罐頂位置的腐蝕速率分別進(jìn)行預(yù)測(cè)。應(yīng)用組合模型的步驟為：①分別輸入LNG儲(chǔ)罐罐底、罐壁以及罐頂3個(gè)位置2007—2016年的腐蝕速率數(shù)據(jù)；②對(duì)灰狼的數(shù)量、尋優(yōu)過程中迭代的次數(shù)、數(shù)據(jù)的維度、灰狼的位置以及尋優(yōu)參數(shù)的范圍進(jìn)行設(shè)定，同時(shí)，根據(jù)狼群所處的區(qū)間，對(duì)每個(gè)灰狼的位置進(jìn)行隨時(shí)調(diào)整；③使用gam以及sig2參數(shù)的初始化處理和LNG儲(chǔ)罐罐底、罐壁以及罐頂3個(gè)位置2007—2016年的腐蝕速率數(shù)據(jù)對(duì)LSSVM模型進(jìn)行訓(xùn)練；④對(duì)適應(yīng)度數(shù)值進(jìn)行計(jì)算，將LSSVM模型運(yùn)行過程中的擬合值與測(cè)量值進(jìn)行比較，將兩者之間的誤差作為適應(yīng)度的函數(shù)；⑤使用GWO算法進(jìn)行全局性質(zhì)的尋優(yōu)，對(duì)灰狼所處的位置進(jìn)行更新，直到滿足算法終止的條件為止，適應(yīng)度最小時(shí)，a等級(jí)灰狼所處的位置即為gam和sig2的數(shù)值；⑥將gam和sig2的最佳數(shù)值輸入到LSSVM模型中，對(duì)2017—2020年LNG儲(chǔ)罐罐底、罐壁以及罐頂3個(gè)位置腐蝕速率進(jìn)行預(yù)測(cè)，并將預(yù)測(cè)結(jié)果與實(shí)際情況進(jìn)行對(duì)比，計(jì)算平均相對(duì)誤差。

3 實(shí)際應(yīng)用研究

3.1 研究對(duì)象

某液化天然氣企業(yè)內(nèi)的LNG儲(chǔ)罐已經(jīng)連續(xù)使用14年，每年均會(huì)對(duì)罐底、罐壁以及罐頂位置的腐蝕速率分別進(jìn)行檢測(cè)，其檢測(cè)結(jié)果見表1。

表1 LNG儲(chǔ)罐腐蝕速率統(tǒng)計(jì)

由表1可見：儲(chǔ)罐3個(gè)位置的腐蝕速率均隨著時(shí)間的延長(zhǎng)而增大。

3.2 GWO-LSSVM模型預(yù)測(cè)

使用GWO-LSSVM組合模型對(duì)儲(chǔ)罐罐底、罐壁以及罐頂進(jìn)行腐蝕速率預(yù)測(cè)，分別將2007—2016年的腐蝕速率數(shù)據(jù)作為整個(gè)模型的學(xué)習(xí)數(shù)據(jù)集，對(duì)2017—2020年的腐蝕速率進(jìn)行全面的預(yù)測(cè)，預(yù)測(cè)過程中使用了MATLAB軟件。組合模型中GWO優(yōu)化算法中的灰狼數(shù)量設(shè)定為50，最大的迭代次數(shù)為20，參數(shù)取值的界限設(shè)定為gam=[10,105]以及sig2=[10-6,106]，所選取得適應(yīng)度函數(shù)為學(xué)習(xí)數(shù)據(jù)集擬合過程中的均方誤差。使用GWO算法對(duì)LSSVM算法中的參數(shù)進(jìn)行優(yōu)化后得到gam=30 621.942 361，sig2=2.745 623，使用最佳的參數(shù)對(duì)LSSVM模型進(jìn)行設(shè)定，并對(duì)2007—2016年罐頂、罐壁以及罐底3個(gè)位置的腐蝕速率進(jìn)行擬合，擬合結(jié)果見圖1。對(duì)3個(gè)位置2017—2020年的腐蝕速率進(jìn)行預(yù)測(cè)，預(yù)測(cè)結(jié)果見圖2。

圖1 2007—2016年GWO-LSSVM擬合結(jié)果

圖2 2017—2020年GWO-LSSVM預(yù)測(cè)結(jié)果

由圖1可見：使用組合模型對(duì)2007—2016年3個(gè)位置腐蝕速率進(jìn)行擬合的過程中出現(xiàn)了過擬合的現(xiàn)象，這主要由于使用GWO算法已經(jīng)找到了LSSVM算法中的最佳參數(shù)，導(dǎo)致在對(duì)訓(xùn)練數(shù)據(jù)進(jìn)行擬合的過程中出現(xiàn)了過擬合問題。

由圖2可見：由于在進(jìn)行數(shù)據(jù)訓(xùn)練的過程中出現(xiàn)了過擬合問題，在使用訓(xùn)練數(shù)據(jù)對(duì)2017—2020年罐頂、罐壁以及罐底等位置進(jìn)行腐蝕速率預(yù)測(cè)的過程中，3個(gè)位置腐蝕速率的預(yù)測(cè)結(jié)果幾乎為直線，無法真實(shí)反映腐蝕速率的波動(dòng)情況，預(yù)測(cè)的結(jié)果的準(zhǔn)確性相對(duì)較差。

3.3 經(jīng)過交叉驗(yàn)證的GWO-LSSVM模型預(yù)測(cè)

為了防止數(shù)據(jù)擬合過程中出現(xiàn)過擬合問題，將所有的數(shù)據(jù)集合分為3個(gè)部分，其中，2007—2013年的腐蝕速率數(shù)據(jù)作為學(xué)習(xí)數(shù)據(jù)集，2014—2016年的腐蝕速率數(shù)據(jù)作為交叉驗(yàn)證數(shù)據(jù)集，2017—2020年的腐蝕速率數(shù)據(jù)作為預(yù)測(cè)數(shù)據(jù)集。在進(jìn)行模型訓(xùn)練的過程中，GWO-LSSVM組合算法的相關(guān)參數(shù)設(shè)定保持不變，將訓(xùn)練數(shù)據(jù)集和預(yù)測(cè)數(shù)據(jù)集兩者產(chǎn)生的均方誤差加和，并將其作為適應(yīng)度函數(shù)，經(jīng)過處理以后，組合模型的擬合結(jié)果見圖3，預(yù)測(cè)結(jié)果見圖4。

圖3 經(jīng)交叉驗(yàn)證后2007—2016年GWO-LSSVM擬合結(jié)果

由圖3可見：在經(jīng)過交叉驗(yàn)證處理后，2007—2016年3個(gè)位置腐蝕速率的擬合結(jié)果相對(duì)較好，擬合結(jié)果的變化趨勢(shì)與真實(shí)的腐蝕速率基本一致，過擬合問題已經(jīng)解決，因此，可以使用交叉驗(yàn)證處理以后的GWO-LSSVM組合對(duì)LNG儲(chǔ)罐腐蝕速率進(jìn)行預(yù)測(cè)。

由圖4可見：在使用交叉驗(yàn)證后的GWOLSSVM組合模型對(duì)儲(chǔ)罐腐蝕速率進(jìn)行預(yù)測(cè)的過程中，預(yù)測(cè)結(jié)果的變化趨勢(shì)與真實(shí)腐蝕速率的變化情況相同，預(yù)測(cè)準(zhǔn)確度相對(duì)較好，其中罐底位置處的預(yù)測(cè)結(jié)果與實(shí)際腐蝕速率最為接近。

圖4 經(jīng)交叉驗(yàn)證后2017—2020年GWO-LSSVM預(yù)測(cè)結(jié)果

3.4 自回歸積分滑動(dòng)平均模型(ARIMA)預(yù)測(cè)

為了進(jìn)一步驗(yàn)證該研究中采用模型的先進(jìn)性，使用ARIMA算法進(jìn)行對(duì)比驗(yàn)證。通過對(duì)3個(gè)位置腐蝕速率的變化情況進(jìn)行分析可以發(fā)現(xiàn)，其數(shù)據(jù)序列十分不平穩(wěn)，腐蝕速率逐年提升，通過對(duì)其自相關(guān)系數(shù)以及偏自相關(guān)系數(shù)進(jìn)行計(jì)算后發(fā)現(xiàn)，腐蝕速率數(shù)據(jù)確實(shí)為非平穩(wěn)的數(shù)據(jù)序列。對(duì)3個(gè)位置腐蝕速率進(jìn)行二階差分處理，處理后的時(shí)序圖見圖5。

由圖5可見：對(duì)罐頂、罐壁及罐底3個(gè)位置多年的腐蝕速率數(shù)據(jù)進(jìn)行二階差分處理，結(jié)果表明其序列相對(duì)較為平穩(wěn)，處理后的數(shù)據(jù)處于上下波動(dòng)狀態(tài)，可以滿足使用ARIMA模型進(jìn)行預(yù)測(cè)的要求。因此，可以構(gòu)建ARIMA模型對(duì)3個(gè)位置的腐蝕速率進(jìn)行預(yù)測(cè)。建立ARIMA模型前，需要計(jì)算3個(gè)位置腐蝕速率的自相關(guān)系數(shù)及偏自相關(guān)系數(shù)，計(jì)算結(jié)果見表2。

圖5 罐頂、罐壁及罐底位置的二階差分時(shí)序圖

由表2可見：罐頂?shù)淖韵嚓P(guān)系數(shù)在第2年位置處最??；罐壁的自相關(guān)系數(shù)在第8年位置處最小，第2年位置處次之；罐底的自相關(guān)系數(shù)在第9年位置處最小，第2年位置處也相對(duì)較小。綜合考慮，將ARIMA模型的自相關(guān)系數(shù)選定為2。按照偏自相關(guān)系數(shù)相對(duì)較大的原則，將ARIMA模型的偏自相關(guān)系數(shù)定為1，最終確定模型為ARIMA(0,2,1)。使用該模型對(duì)2007—2016年3個(gè)位置的腐蝕速率進(jìn)行擬合，擬合結(jié)果見圖6，對(duì)2017—2020年3個(gè)位置的腐蝕速率進(jìn)行預(yù)測(cè)，預(yù)測(cè)結(jié)果見圖7。

表2 二階差分后的自相關(guān)系數(shù)和偏自相關(guān)系數(shù)

由圖6可見：使用該模型對(duì)2007—2016年3個(gè)位置腐蝕速率的擬合效果相對(duì)較好，擬合結(jié)果與實(shí)際腐蝕速率的變化趨勢(shì)基本一致，沒有出現(xiàn)過擬合或者欠擬合的問題，因此，可以使用該種模型對(duì)儲(chǔ)罐3個(gè)位置處的腐蝕速率進(jìn)行預(yù)測(cè)。

圖6 2007—2016年ARIMA模型擬合結(jié)果

由圖7可見：使用ARIMA(0,2,1)對(duì)儲(chǔ)罐3個(gè)位置的腐蝕速率預(yù)測(cè)的結(jié)果與實(shí)際腐蝕速率的變化趨勢(shì)一致，其中，罐底位置處的ARIMA(0,2,1)模型預(yù)測(cè)結(jié)果與實(shí)際腐蝕速率之間的誤差相對(duì)較小。

圖7 2017—2020年ARIMA模型預(yù)測(cè)結(jié)果

僅依靠腐蝕速率的變化趨勢(shì)圖無法對(duì)ARIMA(0,2,1)以及GWO-LSSVM的先進(jìn)性進(jìn)行對(duì)比，因此，需要對(duì)2種模型在各個(gè)位置的預(yù)測(cè)誤差進(jìn)行評(píng)價(jià)。ARIMA(0,2,1)和GWO-LSSVM模型的預(yù)測(cè)誤差對(duì)比情況見表3。

表3 罐頂、罐壁及罐底位置腐蝕速率預(yù)測(cè)結(jié)果對(duì)比

由表3可見：使用交叉驗(yàn)證后的GWOLSSVM模型對(duì)罐頂腐蝕速率預(yù)測(cè)的最大相對(duì)誤差為4.69%，平均相對(duì)誤差為3.10%；對(duì)罐壁腐蝕速率預(yù)測(cè)的最大相對(duì)誤差為6.12%，平均相對(duì)誤差為4.30%；對(duì)罐底腐蝕速率進(jìn)行預(yù)測(cè)的最大相對(duì)誤差為2.90%，平均相對(duì)誤差為2.07%。使用ARIMA(0,2,1)模型對(duì)罐頂位置腐蝕速率預(yù)測(cè)的最大相對(duì)誤差為8.77%，平均相對(duì)誤差為7.17%；對(duì)罐壁位置腐蝕速率預(yù)測(cè)的最大相對(duì)誤差為10.53%，平均相對(duì)誤差為7.16%；對(duì)罐底位置腐蝕速率預(yù)測(cè)的最大相對(duì)誤差為4.69%，平均相對(duì)誤差為3.13%。

將2種模型的預(yù)測(cè)結(jié)果進(jìn)行對(duì)比分析，結(jié)果表明：GWO-LSSVM模型預(yù)測(cè)結(jié)果的誤差波動(dòng)相對(duì)較小，罐頂、罐壁以及罐底3個(gè)位置預(yù)測(cè)結(jié)果的平均相對(duì)誤差均小于ARIMA模型，表明GWO-LSSVM模型相對(duì)較為先進(jìn)，可使用該模型對(duì)儲(chǔ)罐3個(gè)位置的腐蝕速率進(jìn)行預(yù)測(cè)。

4 結(jié)論

針對(duì)LNG儲(chǔ)罐罐頂、罐壁以及罐底3個(gè)位置的腐蝕問題，對(duì)其腐蝕原因進(jìn)行了分析，提出基于GWO-LSSVM組合模型的預(yù)測(cè)算法，通過構(gòu)建組合模型的方式，對(duì)3個(gè)位置2007—2016年的腐蝕速率數(shù)據(jù)進(jìn)行了分析和研究，并對(duì)3個(gè)位置2017—2020年的腐蝕速率進(jìn)行了預(yù)測(cè)，同時(shí)與其他模型計(jì)算結(jié)果進(jìn)行了對(duì)比，得到了以下結(jié)論。

1)LNG儲(chǔ)罐的腐蝕主要可以分為罐底腐蝕、罐壁腐蝕及罐頂腐蝕，其中罐底腐蝕較為嚴(yán)重，主要是由于該位置影響腐蝕的因素較為復(fù)雜；罐頂位置會(huì)聚集大量的水蒸氣以及氧氣，產(chǎn)生電化學(xué)腐蝕；罐底主要以濃差腐蝕為主。

2)單獨(dú)使用GWO-LSSVM模型對(duì)2007—2016年3個(gè)位置腐蝕速率進(jìn)行擬合會(huì)出現(xiàn)過擬合現(xiàn)象，這主要是由于使用GWO算法已經(jīng)找到了LSSVM算法中的最佳參數(shù)，同時(shí)，受過擬合影響，3個(gè)位置腐蝕速率的預(yù)測(cè)結(jié)果幾乎為直線，無法反應(yīng)真實(shí)腐蝕速率的波動(dòng)情況，預(yù)測(cè)的結(jié)果相對(duì)較差。

3)使用經(jīng)過交叉驗(yàn)證后的GWO-LSSVM模型對(duì)數(shù)據(jù)進(jìn)行訓(xùn)練，其擬合結(jié)果相對(duì)較好，預(yù)測(cè)結(jié)果的變化趨勢(shì)與真實(shí)腐蝕速率的變化情況相同，預(yù)測(cè)效果相對(duì)較好。罐頂腐蝕速率預(yù)測(cè)的最大誤差為4.69%，平均相對(duì)誤差為3.10%；罐壁腐蝕速率預(yù)測(cè)的最大誤差為6.12%，平均相對(duì)誤差為4.30%；罐底腐蝕速率預(yù)測(cè)的最大誤差為2.90%，平均相對(duì)誤差為2.07%。3個(gè)位置的預(yù)測(cè)結(jié)果均好于ARIMA模型，表明采用GWO-LSSVM模型對(duì)儲(chǔ)罐3個(gè)位置的腐蝕速率進(jìn)行預(yù)測(cè)具有較強(qiáng)的先進(jìn)性。