朱成杰 王鳴雁 吳章玉

【摘? ?要】? ?針對目前鋰離子電池RUL的預(yù)測中用來表達退化特征的健康因子（HI）存在信息不足的缺陷，提出一種鋰離子電池RUL預(yù)測模型。首先提取出多個可以反映電池退化特征的參數(shù);其次用主成分分析法（PCA）對提取的多個參數(shù)進行去冗余降維并充分包含退化信息;再引入集合經(jīng)驗?zāi)B(tài)分解（EEMD），對HI進行多尺度分解;最后基于長短期記憶網(wǎng)絡(luò)（LSTM）構(gòu)建預(yù)測模型。實驗證明該預(yù)測模型在預(yù)測鋰電池的RUL時有較高的精度、適應(yīng)性較強。

【關(guān)鍵詞】? ?鋰離子電池;長短期記憶網(wǎng)絡(luò);集合經(jīng)驗?zāi)B(tài)分解;主成分分析法

Indirect Prediction of RUL for Li-ion Batteries Based

on Principal Component Analysis

Zhu Chengjie， Wang Mingyan*， Wu Zhangyu

（Anhui University of Science and Technology， Huainan 232001，China）

【Abstract】? ? A lithium-ion battery RUL prediction model is proposed to address the effect of insufficient information on the HI used to express the degradation characteristics in the current prediction of lithium-ion battery RUL. Firstly， several parameters that can reflect the degradation characteristics of the battery are extracted. Second， these multiple parameters were de-redundantly dimensionalized using principal component analysis （PCA）， and the degradation information is fully included. Then， the ensemble empirical modal decomposition （EEMD） is introduced to perform multi-scale decomposition of HI. Finally， the prediction model is constructed based on the LSTM network. Experiments show that the prediction model has high accuracy and strong adaptability in predicting RUL of lithium-ion batteries.

【Key words】? ? ?lithium ion battery; long short term memory; ensemble empirical mode decomposition; principal component analysis

〔中圖分類號〕 TM912? ? ? ? ? ? ? ? ? ? ? ? ?   ? ? ? 〔文獻標識碼〕? A ? ? ? ? ? ? ?〔文章編號〕 1674 - 3229（2022）01- 0030 - 07

0? ? ?引言

鋰離子電池作為目前新能源領(lǐng)域使用最廣泛的能源動力，具有環(huán)境污染小、輸出功率大、循環(huán)使用壽命長[1]、安全和節(jié)省成本等優(yōu)點。因為環(huán)境溫度、充放電電流、電壓等參數(shù)的影響，鋰離子電池容易隨著時間的推移而退化[2]，就有可能會影響使用鋰離子電池作為動力的新能源設(shè)備的正常運轉(zhuǎn)，帶來損失。對鋰離子電池的剩余有用壽命（RUL）進行預(yù)測可以幫助我們有效地了解電池的狀態(tài)及時做好維修更換準備。數(shù)據(jù)驅(qū)動的方法是目前本領(lǐng)域預(yù)測鋰離子電池使用壽命的最常用方法。該方法不用考慮電池內(nèi)部發(fā)生的復(fù)雜化學(xué)反應(yīng)，通過對歷史數(shù)據(jù)進行分析搭建算法模型，從而得到預(yù)測結(jié)果，有更高的泛用性。

數(shù)據(jù)驅(qū)動的方法需要大量的歷史監(jiān)測數(shù)據(jù)，并且預(yù)測精度受其特征參數(shù)影響，因此選取合適的健康因子（Health Indicator， HI）來表征退化特征就變得尤為重要。文獻[3]采用了等壓降放電時間作為間接HI，文獻[4]采用了多個退化參數(shù)融合建立新的HI方式，都取得了一定的成果。本文提出了一種提取多個退化參數(shù)用主成分分析（PCA）對多個參數(shù)降維處理形成新的HI，然后用集合經(jīng)驗?zāi)B(tài)分解（EEMD）的方法將新的HI信號分解為多個尺度，最后通過長短期記憶網(wǎng)絡(luò)（LSTM）對每個信號建立預(yù)測模型，計算鋰離子電池的RUL。

1? ? ?提取退化參數(shù)

實驗數(shù)據(jù)來源于NASA PCoE公布的鋰離子電池數(shù)據(jù)集中BatteryAgingARC-FY08Q4測試組中的05、06、07和18四組電池。四組電池的容量退化情況如圖1所示。NASA的電池數(shù)據(jù)以電池額定容量下降30%為閾值。

本文采用恒流充電時間、放電功率、平均充電電壓升、初始電壓驟降值、等壓降放電時間、放電電壓均方根等6個能代表電池退化狀態(tài)的參數(shù)[5-6]，以B5電池為例進行分析。觀察圖2中B5電池充電電流曲線發(fā)現(xiàn)，鋰離子電池在恒流充電階段，隨著循環(huán)周期的增加，恒流充電時間逐漸減小。

將每個周期的恒流充電時間組成序列用來表征電池的退化特性。隨著鋰離子電池的使用，電池的性能逐漸退化，單位時間內(nèi)的放電量也會隨著變?nèi)?，?dǎo)致放電功率減小。圖3是B5電池充電電壓曲線，從圖中可以發(fā)現(xiàn)隨著電池循環(huán)的增加，充電到達截止電壓的時間越短。平均充電電壓升序列的計算公式如式（1）所示。取1000S到1500S內(nèi)的[m]個數(shù)據(jù)，充電截止電壓為4.2V。

[ACVRi=j=1mVE-Vjm] （1）

式中：[m]為1000S到1500S內(nèi)的電壓個數(shù)，[i]為當前循環(huán)數(shù)，[VE]為充電截至電壓，[Vj]為第[j]個電壓，則平均充電電壓升序列為：

[ACVR=ACVR1，ACVR2，…，ACVRn] ? （2）

圖4為電池放電電壓曲線，隨著放電循環(huán)不斷增加，電池電壓驟降的幅值也在變大。文獻[7]中的研究發(fā)現(xiàn)單個放電周期中電壓從高值下降到低值所需的時間就是放電時間的差值，隨著電池循環(huán)的增加，差值會減小。

[Ti=TiVH-TiVL，i=1，2，3，…，m] ? （3）

（3）式中：[Ti]為第[i]個周期的等壓降放電時間差，[TiVH]為鋰離子電池從放電開始電壓下降到達選定的高電壓時所用的時間，[TiVL]為鋰離子電池從放電開始電壓下降到達選定的低電壓時所用的時間，[m]為電池循環(huán)周期數(shù)。本文取電池從3.8V降至3.5V時對應(yīng)的T1和T2，T1與T2的差值會隨著放電循環(huán)周期的增加不斷縮小。

放電電壓均方根如（4）式，均方根可以表示樣本序列的變化趨勢。

[Xj=1ni=1nVi2] （4）

式中：[n]為電壓樣本個數(shù)，[Vi]為第[i]個樣本，[Xj]為第[j]個循環(huán)。

2? ? ?基于PCA的鋰離子電池健康因子構(gòu)建

為了驗證上述提取的參數(shù)能否正確地表達鋰離子電池的退化特性，計算6個退化參數(shù)的秩相關(guān)系數(shù)，計算公式如式（5）所示。

[RS=i=1n（xi-x）（yi-y）i=1n（xi-x）2i=1n（yi-y）2] （5）

式中：[n]為序列中變量的個數(shù)，[i]為變量在序列中的位置，[x]和[y]為平均值。[RS]取值范圍為（-1，1），值為-1或者1時相關(guān)性最強。計算結(jié)果見表1，可以看出鋰離子電池的各個退化參數(shù)之間存在相關(guān)性，這使得信息之間存有一定的冗余。本文使用主成分分析法來解決這一問題。

表1中C為電池的容量，[Tec]為電池恒流充電時間，P為放電功率，ACVR為平均充電電壓升，ΔV初始電壓驟降值，ΔT為等壓降放電時間差，RMSV為放電電壓均方根。

在使用PCA算法之前，必須對每個退化參數(shù)進行KMO檢驗和Bartlett's檢驗。如果KMO值大于0.7，Bartlett's檢驗的顯著性小于0.01，意味著變量之間存在一定的相關(guān)性，適用于PCA算法。對B5電池提取的各個退化參數(shù)進行檢驗，結(jié)果如表2所示。

分析表2，提取的退化參數(shù)KMO值為0.869，Sig.為0.000。兩個檢驗結(jié)果說明選取的退化參數(shù)適用于PCA算法。

PCA算法可以用來減少和消除數(shù)據(jù)之間的相關(guān)性獲得包含多數(shù)信息的線性不相關(guān)的主成分[8]。PCA算法構(gòu)建健康因子的主要步驟如下[9-10]：

①將提取的退化參數(shù)組成矩陣并對其標準化。設(shè)矩陣X如式（6），其中i=1，2，…，n（n為提取的退化參數(shù)個數(shù)）;j=1，2，…，m（m為退化參數(shù)中電池的指標數(shù)）。

[X=x11x12…x1jx21x22…x2j??…?xi1xi2…xij] （6）

②計算X的協(xié)方差矩陣∑，如式（7）。

[∑=1nXXT] （7）

③計算∑的特征值和特征向量，并將其降序排列。

④計算特征根的貢獻率[Ei]和累計貢獻率E，計算公式如式（8）和式（9）。

[Ei=λii=1Sλi] （8）

[E=i=1mλii=1Sλi] （9）

按上述步驟計算，結(jié)果如表3所示。

表3中第一個主成分的貢獻率已經(jīng)達到98%，基本包含了大多數(shù)的信息。為了驗證由PCA算法構(gòu)建的HI是否能夠表征鋰離子電池的退化狀況，計算其秩相關(guān)系數(shù)如表4所示。PCA算法構(gòu)建的HI與原本的容量的秩相關(guān)系數(shù)均在0.97以上，相關(guān)性較強。

圖5為PCA提取HI與電池實際容量的關(guān)系圖。由圖5可知PCA構(gòu)建的HI與實際容量之間有著很高的相關(guān)性，綜合表4和圖5所述，PCA算法構(gòu)建的HI與原本的容量相關(guān)性較強，因此該HI能表征鋰離子電池退化狀態(tài)，可以用于后續(xù)的RUL預(yù)測。

3? ? ?集合經(jīng)驗?zāi)B(tài)分解算法

由于鋰離子電池的正常使用工況中，存在長時間的擱置階段，該階段存在著容量回升效應(yīng)，導(dǎo)致退化數(shù)據(jù)出現(xiàn)局部波動。用集合經(jīng)驗?zāi)B(tài)分解（Ensemble Empirical Mode Decomposition， EEMD）算法對新構(gòu)建的HI進行分解可以大大降低容量回升引起的局部波動對算法預(yù)測性能的影響。EEMD是Wu和Huang[11]等人針對經(jīng)驗?zāi)B(tài)分解（EMD）提出的一種改進算法。EMD分解的本征模態(tài)函數(shù)（Intrinsic Modal Functions， IMF）包含不同時間尺度的特征，導(dǎo)致模態(tài)混疊[12]。EEMD在原信號中加入了高斯白噪聲信號使整個頻段的極值點均勻分布。EEMD的分解步驟[13]如下：

①在原始信號[S（t）]中加入均值為零的隨機高斯白噪聲[Ni（t）]，式（10）：

[Si（t）=S（t）+Ni（t）] （10）

②對信號進行EMD分解，得到IMF分量[Cij]和殘余分量[Ri（t）]，式（11）。IMF分量滿足的條件：

a. 函數(shù)的局部極值點和過零點的數(shù)量相同或者最多相差一個。

b. 上包絡(luò)線和下包絡(luò)線均值為零。

[Si（t）=i=1nCij+Ri（t）] （11）

③由于高斯白噪聲均值為0，對[i]組分量求均值最終得到EEMD分解的IMF為式（12）：

[Ci=1Ni=1NCij（t）] （12）

將HI的信號進行EEMD分解，結(jié)果如圖6所示。原始信號共分成了5個子信號。其中IMF1-IMF3突變性較強，波形復(fù)雜，在原始信號中占比較小，有明顯的多尺度特征。而IMF4-IMF5信號較平緩，變化幅度不大，在原始信號中占比較大，呈現(xiàn)整個信號的大體趨勢。

4? ? ?長短期記憶網(wǎng)絡(luò)（LSTM）

長短期記憶網(wǎng)絡(luò)（Long Short-Term Memory，LSTM）是基于循環(huán)神經(jīng)網(wǎng)絡(luò)改變而來的。LSTM有著長期記憶的特性[14]，能夠很好地處理時間序列的問題。而鋰離子電池的數(shù)據(jù)有強烈的時間序列性，LSTM在處理這類數(shù)據(jù)的預(yù)測上有著一定的優(yōu)越性。圖7為LSTM結(jié)構(gòu)示意圖，其采用控制門的機制，主要由輸入門、遺忘門和輸出門組成[15]，主要通過這3個門控來控制信息的傳遞與更新，在t時刻的信息更新過程如下。

①通過遺忘門來決定舊細胞傳遞過來的信息遺忘多少，由式（13）求得，其中σ為Sigmoid激活函數(shù)，[Wf]為遺忘門的權(quán)重，[bf]為遺忘門的偏置。

[ft=σ（Wf?[Ht-1，Xt]+bf）] （13）

②計算輸入門[it]的值和細胞更新信息的候選值[Ct]。如式（14）和式（15），其中[Wi]為輸入門的權(quán)重，[bi]為輸入門的偏置，tanh為雙曲正切函數(shù)。

[it=σ（Wi?[Ht-1，Xt]+bi）] （14）

[Ct=tanh（WC?[Ht-1，Xt]+bC）] （15）

③更新細胞狀態(tài)，計算記憶細胞[Ct]的值，見式（16）。

[Ct=ft?Ct-1+it?Ct] （16）

④輸出門確定輸出哪些信息，如式（17）和式（18）。其中[Wo]為輸出門的權(quán)重，[bo]為輸出門的偏置。

[Ot=σ（Wo?[Ht-1，Xt]+bo）] （17）

[Ht=Ot?tanh（Ct）] （18）

5? ? ?實驗的設(shè)計與分析

5.1? ?預(yù)測框架

基于主成分分析的鋰離子電池RUL間接預(yù)測框架如圖8所示，主要分為4個步驟。

①PCA重構(gòu)健康因子。從鋰離子電池的充放電數(shù)據(jù)中提取所需的退化參數(shù)，并對該退化參數(shù)進行評估，再用PCA算法對提取的參數(shù)進行去冗余降維處理。

②構(gòu)建健康因子預(yù)測模型。對新構(gòu)建的HI進行EEMD分解，對每個序列構(gòu)建預(yù)測模型，并將預(yù)測結(jié)果累加組合。

③構(gòu)建RUL預(yù)測模型。將一組的電池全部容量數(shù)據(jù)作為模型的輸入，PCA構(gòu)建的HI作為輸出，訓(xùn)練LSTM網(wǎng)絡(luò)。再將②中得到的預(yù)測數(shù)據(jù)作為輸入得到預(yù)測結(jié)果。

④計算鋰離子電池RUL。為了定量地對模型精度進行分析，本文采用鋰離子電池的RUL誤差以及均方根誤差（RMSE）來作為評價標準，如式（19）和式（20）所示。

[RULerror=RULact-RULpre] （19）

[RMSE=1ni=1nxi-xi2] （20）

5.2? ?結(jié)果與分析

按照上述步驟將NASA數(shù)據(jù)中的B5、B6、B7和B18四組電池數(shù)據(jù)帶入模型，結(jié)果如圖9所示。

圖9中EOL為壽命終止（End of life， EOL），B5、B6和B7號電池組實驗充放電循環(huán)次數(shù)為168，取預(yù)測起點為第102循環(huán)。B18電池組實驗充放電循環(huán)次數(shù)為132，取預(yù)測起點為第82循環(huán)。4組電池的PCA-EEMD-LSTM預(yù)測模型的預(yù)測結(jié)果接近真實的容量值，取得了較好的預(yù)測結(jié)果。

其中B7號電池組在實驗中原始容量曲線最低點為1.40045Ah，高于實驗選取的額定容量70% （1.4Ah）的失效閾值，所以B7號電池并未到達閾值條件。表5中[RULe]為式（19）中的[RULerror];“-”為未到達閾值點。分析表5可知，該方法預(yù)測的RUL誤差均在2個以內(nèi)，RMSE較前兩者均有較大提升，B7由于自身數(shù)據(jù)問題并未到達閾值點，但是RMSE較前兩者依舊有提升。

6? ? ?結(jié)語

本文提出了一種基于主成分分析的間接預(yù)測鋰離子電池RUL的模型，用于預(yù)測鋰離子電池的RUL，并用實驗驗證了模型的適用性和精確性。由于鋰離子電池容量檢測的復(fù)雜性，選用易于監(jiān)測的退化參數(shù)來表征鋰離子電池的退化特征，并用PCA算法去除各退化參數(shù)之間的冗余。構(gòu)建了PCA-EEMD-LSTM預(yù)測模型，用新構(gòu)建的HI作為輸入，得到鋰離子電池的RUL。為了驗證模型的準確性，實驗選取了4組鋰離子電池分別用LSTM、PCA-LSTM和PCA-EEMD-LSTM模型進行分析。實驗結(jié)果驗證了該模型能夠有效地預(yù)測鋰離子電池的RUL，且適用性較廣、預(yù)測精度較高。

[參考文獻]

[1] 蔡艷平，陳萬，蘇延召，等.鋰離子電池剩余壽命預(yù)測方法綜述[J].電源技術(shù)，2021，45（5）：678-682.

[2] 桂長清.溫度對LiFePO_4鋰離子動力電池的影響[J].電池，2011，41（2）：88-91.

[3] Datong Liu. Satellite Lithium-Ion Battery Remaining Cycle Life Prediction with Novel Indirect Health Indicator Extraction[J]. Energies， 2013， 6（8） ： 3654-3668.

[4] Yang Zhang，Bo Guo. Online Capacity Estimation of Lithium-Ion Batteries Based on Novel Feature Extraction and AdaptiveMulti-Kernel Relevance Vector Machine[J]. Energies， 2015， 8（11） ： 12439-12457.

[5] 史永勝，施夢琢，丁恩松，等.基于多退化特征的鋰離子電池剩余壽命預(yù)測[J].電源技術(shù)，2020，44（6）：836-840.

[6] 龐曉瓊，王竹晴，曾建潮，等.基于PCA-NARX的鋰離子電池剩余使用壽命預(yù)測[J].北京理工大學(xué)學(xué)報，2019，39（4）：406-412.

[7] 龐景月，馬云彤，劉大同，等.鋰離子電池剩余壽命間接預(yù)測方法[J].中國科技論文，2014，9（1）：28-36.

[8] 吳曉冬，李小軍.基于PCA和SVM的鋰電池電極涂布質(zhì)量預(yù)測[J].工業(yè)技術(shù)創(chuàng)新，2019，6（3）：56-59.

[9] 宋哲，高建平，潘龍帥，等.基于主成分分析和改進支持向量機的鋰離子電池健康狀態(tài)預(yù)測[J].汽車技術(shù)，2020（11）：21-27.

[10] Fengtao Wang，Bei Wang，Bosen Dun，et al. Remaining life prediction of rolling bearing based on PCA and improved logistic regression model[J]. Journal of Vibroengineering，2016，18（8）：5192-5203.

[11] ZHAOHUA WU，NORDEN E. HUANG. ENSEMBLE EMPIRICAL MODE DECOMPOSITION： A NOISE-ASSISTED DATA ANALYSIS METHOD[J]. Advances in Adaptive Data Analysis，2009，1（1）.

[12] 雷炳銀，王子馳，蘇雨晴，等.基于EEMD-CS-LSSVM的短期負荷預(yù)測方法研究[J].電力系統(tǒng)及其自動化學(xué)報，2021，33（9）：117-122.

[13] 張飛，張志偉，萬樂斐，等.基于EEMD-GRNN方法的光伏電站短期出力預(yù)測[J].太陽能學(xué)報，2020，41（12）：103-109.

[14] 李瑞津，劉斌，張學(xué)敏，等.基于改進LSTM的變電站鉛酸電池壽命預(yù)測[J].電池，2020，50（6）：560-564.

[15] 吉祥飛，李明東，陶衛(wèi)國，等.基于多通道LSTM-CNN模型的Twitter情感分析[J].廊坊師范學(xué)院學(xué)報（自然科學(xué)版），2019，19（2）：21-24+37.