謝 旭,蒲嫻怡,2,畢貴紅,王 凱,高 晗,3
(1.昆明理工大學(xué)電力工程學(xué)院,云南昆明 650500;2.云南電網(wǎng)有限責(zé)任公司玉溪供電局,云南玉溪 653100;3.云南電網(wǎng)有限責(zé)任公司昆明供電局,云南昆明 650011)
在鋰電池使用過程中,由于內(nèi)部老化現(xiàn)象,其提供能量和存儲能量的能力會降低[1]。為了進一步優(yōu)化鋰電池的性能,延長其生命周期,對電池健康狀態(tài)(SOH)進行監(jiān)控是至關(guān)重要的[2]。SOH通常指鋰電池相對于初始值存儲能量的能力,一般使用容量或電阻等參數(shù)來對SOH進行定量[3]。當(dāng)使用容量作為電池SOH指標(biāo)時,SOH即為當(dāng)前循環(huán)的電池容量與壽命開始時的額定容量或電池制造商提供的初始容量的比值[4]。當(dāng)前,一些常見的SOH評估方法包括:直接測量法,一般直接測量電池的關(guān)鍵參數(shù)例如內(nèi)阻[5-6]、阻抗[7]、開路電壓[8]和充放電電流[9],再使用公式計算SOH,這類方法要滿足嚴(yán)格的外部條件,且精度極度依賴于測量技術(shù),魯棒性較差,只適用于實驗室條件;基于模型方法,如電化學(xué)模型[10-11]和等效電路模型[12],這類方法普適性較差,針對不同電池型號和工作狀態(tài)需要不斷調(diào)整模型,因此精度難以保障,且模型復(fù)雜、計算成本較大;數(shù)據(jù)驅(qū)動方法,通過構(gòu)建模型捕獲各種電池參數(shù)中的復(fù)雜非線性關(guān)系,獲得較好的估算結(jié)果,因此近些年來被廣泛應(yīng)用[13]。
基于數(shù)據(jù)驅(qū)動算法通常使用溫度、放電深度、充電速率和平均放電電壓、平均放電溫度、容量及容量差等作為特征向量對SOH進行估計[14-16]。但這些指標(biāo)不能直接表征電池的性能退化程度,對于在線檢測存在一定困難,有些甚至不能夠通過在線測量獲取。為此,龐等[17]提出使用電池循環(huán)對應(yīng)周期的等壓降放電時間作為表征向量,建立SOH與等壓降放電時間的關(guān)系模型。
然而,鋰離子電池是一種復(fù)雜的電化學(xué)器件,它在不同的內(nèi)外條件下具有獨特的非線性運行狀態(tài)[18],隨著循環(huán)進行,其容量及等壓降放電時間的測量值會受多種因素影響,導(dǎo)致數(shù)據(jù)存在再生現(xiàn)象及噪聲波動。目前,應(yīng)對這一問題的有效方法為結(jié)合模式分解算法和機器學(xué)習(xí)的混合模型。胡等[19]提出使用經(jīng)驗?zāi)B(tài)分解(EEMD)對獲取的鋰電池能量數(shù)據(jù)進行特征提取,并結(jié)合深度學(xué)習(xí)得到了良好的實驗結(jié)果。曾等[20]提出了一種變分模態(tài)分解可變模式分解(VMD)和高斯過程回歸的融合方法,通過分解將原始數(shù)據(jù)的影響因素分開,精度較高。殷等[21-22]提出了二層分解技術(shù),并將該技術(shù)應(yīng)用于新能源技術(shù)、電力負(fù)荷等領(lǐng)域,通過二層分解技術(shù)處理后的數(shù)據(jù)能夠具備更為明顯的特征性和規(guī)律性,有助于模型性能的提升。
以上提出的一系列傳統(tǒng)的機器學(xué)習(xí)方法中,其模型建立容易受數(shù)據(jù)集影響,其網(wǎng)絡(luò)參數(shù)需要人工依賴經(jīng)驗設(shè)定,且不能很好地反映長時序數(shù)據(jù)間的關(guān)系,精度難以保證。吳等[23]使用螢火蟲算法對BP 神經(jīng)網(wǎng)絡(luò)進行參數(shù)優(yōu)化,實驗效果良好。肖等[24]和王等[15]使用長短期記憶(LSTM)神經(jīng)網(wǎng)絡(luò)對SOH進行估算,實驗證明了該算法的有效性。
結(jié)合上述問題,本文提出了一種基于二層分解技術(shù)的鋰離子電池SOH估算方法。從容量的角度定義SOH,選擇等壓降放電時間作為健康因子,該健康因子能夠?qū)崟r獲取,且能較好地反映容量數(shù)據(jù)存在的再生現(xiàn)象。通過改進快速集合經(jīng)驗?zāi)B(tài)(FEEMD)與相關(guān)性分析(CA)將等壓降放電時間序列進行多尺度分解,分離出波動和主趨勢數(shù)據(jù),然后使用VMD 將波動趨勢數(shù)據(jù)分解成一系列模態(tài)分量。二層分解技術(shù)進一步將波動數(shù)據(jù)中的特征信息提取出來,降低噪聲及再生現(xiàn)象影響的同時,能夠?qū)⒂行У臄?shù)據(jù)用作模型訓(xùn)練的輸入,減少了數(shù)據(jù)量的冗余。采用粒子群算法(PSO)的LSTM 模型建立等壓降放電時間波動趨勢分量、主趨勢數(shù)據(jù)與電池容量間的關(guān)系,PSO 能夠通過計算獲得最優(yōu)網(wǎng)絡(luò)參數(shù),進一步提高了模型的估算精度。最終,使用NASA 鋰離子電池數(shù)據(jù)集驗證了該方法的有效性,并與現(xiàn)有方法進行了對比實驗。
鋰離子電池SOH一般通過容量或內(nèi)阻衡量[25-26],然而這兩個參數(shù)難以直接在線獲取。選擇可以在線測量的等壓降放電時間作為健康因子,能夠反映電池的局部再生現(xiàn)象,非常適合作為鋰離子電池SOH估算的表征量。以NASA 鋰離子電池數(shù)據(jù)集5 號電池為例,如圖1 所示,選取放電區(qū)間為3.7~4 V,得到等壓降放電時間曲線,等壓降放電時間曲線與5號電池的實際容量退化趨勢相似。所以,選擇將容量作為輸出,從容量的角度定義SOH并實現(xiàn)估算[27]。

圖1 NASA鋰離子電池數(shù)據(jù)集5號電池等壓降放電時間和容量變化曲線
本文提出了一種基于二層分解技術(shù)的粒子群優(yōu)化LSTM神經(jīng)網(wǎng)絡(luò)鋰電池SOH估算方法,模型實驗方法如圖2 所示,包括數(shù)據(jù)處理、估算模型及優(yōu)化、SOH估算三個模塊。

圖2 FEEMD-VMD-PSO-LSTM估算方法
數(shù)據(jù)處理部分使用了FEEMD 和VMD 結(jié)合的二層分解方法,對等壓降放電時間D(t)進行多尺度分解,以降低各分量數(shù)據(jù)的復(fù)雜性,提高模型估算精度。FEEMD[28]相較于EEMD[29]優(yōu)化了停止判據(jù)準(zhǔn)則,對信號的處理速度更快,能夠滿足數(shù)據(jù)處理的實時性需求。VMD[30]將信號分解轉(zhuǎn)化為對變分約束模型求最優(yōu)解的過程,能夠自適應(yīng)將信號分解成一系列固定模式分量。首先,利用FEEMD 對D(t)進行一次分解,得到一組本征模態(tài)分量(IMF);其次,利用CA 將IMF重構(gòu),獲得等壓降放電時間序列的全局退化趨勢和波動趨勢;最后,使用VMD 對一次分解獲得的復(fù)雜波動趨勢再次分解,被二次分解后的各波動分量規(guī)律性增強,復(fù)雜性降低,可進一步提高各分量的估算精度。
考慮到數(shù)據(jù)的長時序特點,使用LSTM 神經(jīng)網(wǎng)絡(luò)作為估算模型,并加入了PSO 算法以進一步提高模型的估算精度。LSTM 神經(jīng)網(wǎng)絡(luò)解決了循環(huán)神經(jīng)網(wǎng)絡(luò)(RNN)存在的梯度消失、梯度爆炸等問題,且對長時間依賴問題,有著良好的計算效果,它能夠篩選歷史網(wǎng)絡(luò)單元的有用信息,捕捉數(shù)據(jù)序列中更有意義、更清晰的信息;PSO 算法能夠通過模型訓(xùn)練自動尋找LSTM 神經(jīng)網(wǎng)絡(luò)中的最優(yōu)超參數(shù),避免了傳統(tǒng)人工憑借經(jīng)驗選擇超參數(shù)的隨機性、不確定性。
用容量來表征SOH,使用全局退化趨勢數(shù)據(jù)和波動趨勢分量經(jīng)VMD 分解后得到的二層本征模態(tài)分量組成輸入數(shù)據(jù)集,對應(yīng)的電池容量組成輸出數(shù)據(jù)集,利用PSO 優(yōu)化獲得的PSO-LSTM 模型對電池容量進行估算。
實驗數(shù)據(jù)來自于美國航空航天局卓越故障預(yù)測研究中心(NASA PCoE)。使用4 塊額定容量2 Ah 的18650 型鋰離子電池(B5、B6、B7、B18)在室溫24 ℃下進行充電、放電和阻抗測量實驗,并記錄監(jiān)測數(shù)據(jù)。本文使用B5 和B6 兩塊電池數(shù)據(jù)進行實驗。
以NASA 鋰離子電池數(shù)據(jù)集5 號電池為例,某一次放電循環(huán)Ti(i=1,2···T)對應(yīng)的PSO 優(yōu)化的LSTM 估算方法輸入變量如表1 所示。輸出值是5 號電池放電循環(huán)Ti的對應(yīng)容量值。

表1 方法的輸入變量
以5 號電池為例,F(xiàn)EEMD 分解等壓降放電時間序列得到6 個分量和剩余量,相關(guān)性分析計算結(jié)果如表2 所示,關(guān)聯(lián)閾值Th為0.396。IMF5、IMF6 分量與剩余量疊加組成全局退化趨勢,剩余分量疊加作為波動趨勢,如圖3 所示。VMD 分解波動趨勢,得到若干模態(tài)分量,如圖4 所示。VMD 將波動趨勢的局部特征分別體現(xiàn)在各分量當(dāng)中,主退化趨勢呈現(xiàn)很好地單調(diào)下降趨勢。

表2 FEEMD 本征模態(tài)函數(shù)的相關(guān)系數(shù)

圖3 全局退化趨勢和波動趨勢

圖4 波動趨勢VMD分解結(jié)果
分別將兩塊電池數(shù)據(jù)中的前80 次放電循環(huán)數(shù)據(jù)作為訓(xùn)練數(shù)據(jù)集,后88 次放電循環(huán)數(shù)據(jù)作為測試數(shù)據(jù)集建立PSOLSTM 估算方法。
誤差評價函數(shù)使用平均絕對值百分比誤差(MAPE)和均方根誤差(RMSE),如式(1)~(2)所示:

式中:Ci表示第i次放電循環(huán)對應(yīng)容量的真實值;表示第i次放電循環(huán)對應(yīng)容量的估算值;N表示放電循環(huán)次數(shù)。

表3 為PSO-LSTM 估算方法的網(wǎng)絡(luò)參數(shù)設(shè)置表。
使用以上數(shù)據(jù)和網(wǎng)絡(luò)設(shè)置參數(shù)構(gòu)建SOH估算方法。為驗證所提模型對SOH估算的精確性,建立了FEEMD-LSTM、FEEMD-VMD-LSTM、FEEMD-VMD-PSO-LSTM 共3 種模型進行對比實驗,結(jié)果如圖5 所示,文中提出方法得到的容量曲線能夠準(zhǔn)確表現(xiàn)真實容量退化中存在的再生現(xiàn)象,這說明經(jīng)過二層分解后的數(shù)據(jù)更能表現(xiàn)局部特征,這有利于神經(jīng)網(wǎng)絡(luò)更加有針對性地進行訓(xùn)練,相較于單獨使用FEEMD 方法,二層分解技術(shù)有效抑制了模態(tài)混疊的出現(xiàn),避免了由于再生現(xiàn)象或模態(tài)混疊現(xiàn)象對估算模型的精度造成不良影響。

圖5 5號、6號電池容量估算結(jié)果
3 種模型的SOH估算精度如表4 所示,與其他模型相比,F(xiàn)EEMD-VMD-PSO-LSTM 模型SOH估算誤差最小。以5 號電池為例,比較FEEMD-LSTM 和FEEMD-VMD-LSTM 兩模型誤差精度,后者的MAPE降低明顯;比較FEEMD-VMDLSTM 和FEEMD-VMD-PSO-LSTM 兩模型誤差精度,加入粒子群優(yōu)化的模型精度提升明顯,說明使用PSO 進行網(wǎng)絡(luò)參數(shù)優(yōu)化相較于人工經(jīng)驗設(shè)置參數(shù)要更具科學(xué)化,能夠有效提升估算精度。

表4 5 號和6 號電池的估算精度
本文選擇等壓降放電時間作為鋰離子電池健康因子,提出了一種基于二層分解技術(shù)的粒子群優(yōu)化長短期記憶循環(huán)神經(jīng)網(wǎng)絡(luò)的估算方法。
本文提出的方法選擇了與容量數(shù)據(jù)特征擬合較好的等壓降放電時間構(gòu)建與SOH的關(guān)系,使用二層分解技術(shù),將等壓降放電時間數(shù)據(jù)中的全局退化趨勢和波動趨勢很好地分離開,降低使用單一分解方法作為特征輸入量的隨機性、波動性,通過均方根誤差評價參數(shù)的對比,能夠體現(xiàn)出FEEMDVMD-LSTM 相較于FEEMD-LSTM 方法的優(yōu)勢,5 號和6 號電池誤差分別降低了2%和9.7%。
LSTM 針對時序數(shù)據(jù)表現(xiàn)良好,能夠有效捕捉數(shù)據(jù)中的時序特征,使用PSO 進行LSTM 神經(jīng)網(wǎng)絡(luò)超參數(shù)優(yōu)化有利于提高估算精度,且提升較為明顯,通過均方根誤差評價參數(shù)的對比,體現(xiàn)出FEEMD-VMD-PSO-LSTM 相較于FEEMDVMD-LSTM 方法的優(yōu)勢,5 號和6 號電池誤差分別降低了3.3%和1%。
總的來看,F(xiàn)EEMD-VMD-PSO-LSTM 方法相較于FEEMD-LSTM 方法,5 號和6 號電池均方根誤差分別降低了5.3%和10.5%,這說明文中提出的估算方法在多樣本電池數(shù)據(jù)中也能夠保證較好的精度,能夠為電池管理系統(tǒng)提供一定的決策依據(jù)。