考慮不同充電策略的鋰電池健康狀態(tài)區(qū)間估計(jì)

張孝遠(yuǎn), 張金浩, 楊立新

(河南工業(yè)大學(xué) 電氣工程學(xué)院,鄭州 450001)

在碳達(dá)峰、碳中和目標(biāo)背景下,電動(dòng)汽車(chē)得到大力發(fā)展.鋰離子電池作為電動(dòng)汽車(chē)的核心部件之一,使用量呈現(xiàn)井噴式增長(zhǎng).與此同時(shí),由鋰離子電池健康狀態(tài)(State of Health, SOH)劣化引發(fā)的安全問(wèn)題日益凸顯,在一定程度上限制了電動(dòng)汽車(chē)行業(yè)的發(fā)展.SOH通過(guò)百分比的形式量化電池的老化情況,是評(píng)價(jià)電池老化程度的綜合評(píng)價(jià)標(biāo)準(zhǔn)[1].合理精確的SOH估計(jì)有助掌握電池的衰減情況,對(duì)電池安全使用、維護(hù)管理和經(jīng)濟(jì)性評(píng)估都有十分重要的意義[2-3].

SOH無(wú)法通過(guò)傳感器直接測(cè)量,現(xiàn)有SOH評(píng)估方法主要分為模型法和數(shù)據(jù)驅(qū)動(dòng)法[4].模型法基于電化學(xué)機(jī)理建立物理等效電路模型來(lái)描述鋰離子電池SOH的變化,模型參數(shù)可依據(jù)電池的不同衰退情況而變化,但鋰電池內(nèi)部化學(xué)變化是一個(gè)動(dòng)態(tài)、非線性的復(fù)雜電化學(xué)系統(tǒng),建立準(zhǔn)確的電化學(xué)模型非常困難[5].相較之下,數(shù)據(jù)驅(qū)動(dòng)法無(wú)需考慮電池內(nèi)部復(fù)雜的電化學(xué)機(jī)理,僅基于大量充放電數(shù)據(jù)提取能夠表征電池SOH的特征參數(shù),即可建立特征參數(shù)與SOH之間的定量映射關(guān)系.因此,基于數(shù)據(jù)驅(qū)動(dòng)的鋰電池SOH估計(jì)得到眾多研究者的青睞[6-7].

數(shù)據(jù)驅(qū)動(dòng)法一般包含特征參數(shù)選取和SOH回歸估計(jì)兩大步驟.提取出能夠充分表征SOH的特征參數(shù)是近年來(lái)的一個(gè)研究熱點(diǎn),文獻(xiàn)[8]中從充放電曲線中獲取增量容量(Incremental Capacity, IC)曲線,研究了IC曲線隨循環(huán)次數(shù)的變化趨勢(shì),并驗(yàn)證IC峰值、IC峰值位置等特征與SOH之間的相關(guān)性,但不同型號(hào)電池的IC峰值數(shù)量并不一致,所用特征并不適用于所用電池.文獻(xiàn)[9]中進(jìn)一步研究不同型號(hào)電池包含多個(gè)IC峰值的情況.但在IC曲線的提取過(guò)程中必須設(shè)計(jì)復(fù)雜的濾波算法來(lái)平滑原始電壓測(cè)量值,同時(shí),IC峰值位置會(huì)隨著電池老化而逐漸偏移;并且要獲得所有老化周期下的IC特性,需要在較大的電壓范圍內(nèi)進(jìn)行采樣,十分耗時(shí).針對(duì)這個(gè)問(wèn)題,文獻(xiàn)[10]中從放電電壓曲線及溫度曲線中獲取放電平均電壓和平均溫度作為外部特征,估算SOH.但在實(shí)際應(yīng)用中,電池放電深度不一致導(dǎo)致平均電壓變化且溫度受環(huán)境影響較大.為此,文獻(xiàn)[11]中選擇從恒壓充電曲線中獲取平均電流等作為外部特征,并用等效電路模型表征恒壓曲線中的電池老化現(xiàn)象.上述鋰電池特征參數(shù)提取的最新研究較大程度上推動(dòng)了基于數(shù)據(jù)驅(qū)動(dòng)的鋰電池SOH估計(jì)的發(fā)展.但上述研究有一個(gè)共性不足——都是基于恒流恒壓充電策略下進(jìn)行的鋰電池充放電實(shí)驗(yàn),所用數(shù)據(jù)均為單一充電策略下所得數(shù)據(jù),未考慮不同充電策略或鋰電池實(shí)際放電情況等問(wèn)題,在實(shí)際應(yīng)用中缺乏普適性.研究高效準(zhǔn)確的適用不同充電策略的鋰電池特征參數(shù)選擇方法是迫切的現(xiàn)實(shí)需求.

獲得與SOH高度相關(guān)的特征參數(shù)之后,需要建立所獲取特征參數(shù)與SOH之間的映射關(guān)系,其本質(zhì)是一類(lèi)多元回歸問(wèn)題.諸多機(jī)器學(xué)習(xí)方法可以用于解決這類(lèi)問(wèn)題.反向傳播(BP)神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)、高斯回歸、決策樹(shù)等方法相繼被用來(lái)建立SOH的回歸估計(jì)模型[4].如文獻(xiàn)[12]中搭建BP神經(jīng)網(wǎng)絡(luò)對(duì)鋰電池SOH進(jìn)行在線估計(jì);文獻(xiàn)[3]中采用粒子群優(yōu)化的支持向量回歸模型實(shí)時(shí)估計(jì)不同老化狀態(tài)下的電池SOH;文獻(xiàn)[13]中基于高斯過(guò)程回歸(Gaussian Process Regression, GPR)建立SOH估算模型,取得較精確的估算精度;文獻(xiàn)[10]中搭建布谷鳥(niǎo)搜索支持向量模型實(shí)現(xiàn)鋰離子電池的SOH快速估計(jì).結(jié)合前述特征參數(shù)選取方法,上述回歸估計(jì)模型在特定條件下能夠較為精確地估計(jì)SOH值,但上述研究所采用的SOH估計(jì)模型均為點(diǎn)估計(jì)方法.實(shí)際工作環(huán)境中鋰電池剩余容量存在容量回升效應(yīng)[14],隨著循環(huán)次數(shù)的推移,鋰電池剩余容量會(huì)有一定幅度的波動(dòng)甚至上升,導(dǎo)致在某一時(shí)刻的SOH估計(jì)結(jié)果無(wú)法充分揭示SOH的退化趨勢(shì).同時(shí),在實(shí)際運(yùn)行中電動(dòng)汽車(chē)工作環(huán)境復(fù)雜、負(fù)載波動(dòng)較大,導(dǎo)致鋰電池退化模式不固定,SOH估計(jì)結(jié)果具有一定的不確定性[15],傳統(tǒng)點(diǎn)估計(jì)方法無(wú)法對(duì)鋰電池退化過(guò)程進(jìn)行合理量化.

針對(duì)鋰電池SOH估計(jì)的不確定性量化估計(jì)問(wèn)題,近年來(lái)有研究人員提出兩類(lèi)方法來(lái)構(gòu)建電池SOH的預(yù)測(cè)區(qū)間.第一類(lèi)為理論方法,基于以下假設(shè):預(yù)測(cè)誤差遵循具有零均值的確定分布,通常為正態(tài)分布或拉普拉斯分布[16].文獻(xiàn)[17]中開(kāi)發(fā)了一個(gè)基于平均熵和相關(guān)向量機(jī)的多步預(yù)測(cè)模型,并將其應(yīng)用于電池的SOH評(píng)估和剩余壽命預(yù)測(cè),該方法基于高斯分布構(gòu)造95%的預(yù)測(cè)結(jié)果置信區(qū)間.但實(shí)測(cè)的鋰電池充放電數(shù)據(jù)涉及復(fù)雜過(guò)程和噪聲干擾,很難確保滿(mǎn)足高斯假設(shè).第二類(lèi)為神經(jīng)網(wǎng)絡(luò)方法,該類(lèi)方法無(wú)需考慮誤差分布,因此受到研究者青睞.如文獻(xiàn)[18]中構(gòu)建了一個(gè)改進(jìn)的上下邊界估計(jì)神經(jīng)網(wǎng)絡(luò)對(duì)鋰電池SOH進(jìn)行區(qū)間估計(jì).但此算法需要大量數(shù)據(jù)集進(jìn)行訓(xùn)練,而鋰電池全生命期,即SOH從100%退化到80%,一般循環(huán)次數(shù)不超過(guò) 2 000,屬于典型的小樣本問(wèn)題,這類(lèi)方法效果有限.

針對(duì)當(dāng)前有關(guān)鋰電池SOH估計(jì)方面的研究在特征參數(shù)選取普適性和SOH估計(jì)不確定性量化兩方面的不足,提出一種適用于多種充電策略的鋰電池SOH不確定量化估計(jì)方法.該方法充分考慮當(dāng)前常用充電策略下充放電曲線的變化,從電壓、電流及IC曲線中提取出多個(gè)外部特征,再通過(guò)交叉驗(yàn)證方法自動(dòng)選取出適用于不同充電策略的個(gè)性化最優(yōu)特征組合.首次提出使用支持向量分位數(shù)回歸(Support Vector Quantile Regression, SVQR)模型進(jìn)行鋰電池SOH區(qū)間估計(jì).所采用的SVQR方法利用分位數(shù)回歸(Quantile Regression, QR)模型解決了區(qū)間分布問(wèn)題,同時(shí),支持向量回歸(Support Vector Regression, SVR)模型在小樣本充放電數(shù)據(jù)下具有較好的泛化性能.將所提方法在常用的3種不同充電策略的鋰電池?cái)?shù)據(jù)集上進(jìn)行驗(yàn)證,并與QR、分位數(shù)回歸神經(jīng)網(wǎng)絡(luò)(Quantile Regression Neural Network,QRNN)及GPR進(jìn)行對(duì)比,驗(yàn)證了所提方法的優(yōu)越性.

1 所提方法的總體框架

所提方法的總體流程如圖1所示,分為離線訓(xùn)練和在線評(píng)估兩個(gè)階段.在離線訓(xùn)練階段,首先從循環(huán)電壓曲線、電流曲線及IC曲線中提取多個(gè)通用特征參數(shù),然后采用交叉驗(yàn)證方法自動(dòng)選擇與SOH相關(guān)性最好的最優(yōu)特征參數(shù)組合,并將其輸入SVQR模型進(jìn)行訓(xùn)練.在線階段將充放電數(shù)據(jù)輸入訓(xùn)練好的SVQR模型,實(shí)現(xiàn)電池SOH的在線區(qū)間估計(jì),從而給出鋰電池使用安全指導(dǎo)建議.

圖1 所提SOH區(qū)間估計(jì)方法流程Fig.1 Flow chart of the proposed SOH interval estimation method

2 考慮不同充電策略的鋰電池SOH 區(qū)間估計(jì)方法

2.1 鋰電池常用充電策略及其特點(diǎn)

鋰電池充電速度是制約電動(dòng)汽車(chē)進(jìn)一步發(fā)展的主要瓶頸之一[19].人們希望提高充電速度,但充電速度的提高會(huì)造成溫度的快速升高和電池壽命的快速衰退.因此,探索能夠兼顧充電快速性、溫度控制與合理壽命預(yù)期的充電策略一直是一個(gè)研究熱點(diǎn)[20].目前鋰電池充電策略已由最初的恒壓充電、恒流充電,發(fā)展到恒壓恒流充電結(jié)合的多階段充電[21].結(jié)合當(dāng)前鋰電池充電策略發(fā)展的現(xiàn)狀,重點(diǎn)考慮如下3類(lèi)常用的充電策略.第一類(lèi)為恒流恒壓(Constant Current Constant Voltage, CC-CV)充電.CC-CV充電第一階段以恒定電流充電;當(dāng)電壓達(dá)到預(yù)定值時(shí)轉(zhuǎn)入第二階段進(jìn)行恒壓充電,此時(shí)電流逐漸減小;當(dāng)充電電流下降到0時(shí),蓄電池完全充滿(mǎn).該充電策略可有效避免電池極化效應(yīng),但充電時(shí)間較長(zhǎng);第二類(lèi)為恒流(Constant Current, CC)充電.CC充電采用恒定電流對(duì)電池進(jìn)行充電.由于電流一般較大,電池會(huì)快速達(dá)到截止電壓,充電時(shí)間較短,但控制不得當(dāng)則會(huì)增大電芯的極化;第三類(lèi)為兩階段快速充電(Two Step Fast Charge, 2-step).該充電策略歷經(jīng)兩個(gè)充電階段.充電過(guò)程可表達(dá)為“C1(Q1)-C2”,其中C1和C2分別代表兩個(gè)充電階段,Q1代表充電電流切換時(shí)電池的荷電狀態(tài)(State of Charge, SOC),第一個(gè)階段采用恒流充電至SOC為Q1結(jié)束;第二個(gè)階段采用CC-CV方式直至電池充滿(mǎn).上述3種充電策略的電壓變化曲線如圖2所示.

圖2 不同充電策略下電壓變化Fig.2 Voltage change of different charging strategies

2.2 SOH特征選擇

鋰電池充放電曲線中蘊(yùn)含豐富的電池老化信息.圖3為某單體電池不同循環(huán)次數(shù)下電壓、電流變化曲線.可見(jiàn),隨著循環(huán)次數(shù)增多,充放電電壓曲線和電流曲線也隨之變化,充放電曲線的形狀與電池的老化情況高度相關(guān).如何從充放電曲線中提取出與SOH高度敏感的特征參數(shù)是SOH估計(jì)的關(guān)鍵.但不同充電策略所對(duì)應(yīng)的敏感特征參數(shù)可能不同,很難找到普適于多個(gè)充電策略的特征參數(shù).當(dāng)前有關(guān)鋰電池SOH估計(jì)方法大多針對(duì)某一特定充電策略,沒(méi)有考慮充電策略的不同及其進(jìn)一步發(fā)展.為此,提出一種更泛化的特征參數(shù)提取流程.首先,充分分析不同充電策略的單體電池在不同循環(huán)次數(shù)下電壓、電流、IC曲線變化特性,并參考文獻(xiàn)[2,6,15]中的結(jié)論,從IC曲線提取IC曲線峰值(PIC)、IC曲線峰值對(duì)應(yīng)電壓(UIC)、IC曲線峰左斜率(kL)、IC曲線峰右斜率(kR) 4個(gè)特征,如圖4(a)所示;從電流曲線提取CV充電下降到給定電流所用時(shí)間(tD)和CV充電過(guò)程最后階段的平均電流(IA),如圖4(b)所示;從電壓曲線提取CC充電上升到給定電壓所用時(shí)間(tU)和單次循環(huán)后靜置時(shí)間內(nèi)的電壓回升數(shù)值(UV),如圖4(c)所示.然后,為了適應(yīng)不同充電策略,采用交叉驗(yàn)證方法自動(dòng)選擇針對(duì)特定充電策略的個(gè)性化最優(yōu)特征參數(shù).交叉驗(yàn)證方法在上述所提取的通用特征集中遍歷所有可能的特征參數(shù)組合, 然后對(duì)各個(gè)特征組合的估計(jì)結(jié)果進(jìn)行排序以選出最優(yōu)特征參數(shù)組合.

圖3 某單體電池充放電曲線在不同循環(huán)周期下的變化情況Fig.3 Change of voltage and current in charging and discharging of a single battery of different cycles

圖4 電流、電壓和IC曲線的提取特征Fig.4 Extracted features from current, voltage, and IC curve

2.3 分位數(shù)回歸

i=1, 2, …,t;τ=τ1,τ2, …,τM

(1)

式中:xi=[xi,1xi,2…xi,d]為包含d個(gè)分量的獨(dú)立變量,xi,j是xi的第j個(gè)分量,j=1, 2, …,d;yi為對(duì)應(yīng)xi的因變量;qτ(xi)為因變量yi在自變量xi下的第τ分位數(shù),對(duì)于任意分位數(shù)0<τ<1,β(τ)=[β1(τ)β2(τ) …βd(τ)]′為τ分位數(shù)下的回歸系數(shù).β(τ)的估計(jì)值通過(guò)最小化檢驗(yàn)函數(shù)ρτ(σ)獲得.

(2)

檢驗(yàn)函數(shù)ρτ(σ)為一分段線性函數(shù),定義為

ρτ(σ)=σ(τ-I(σ))

(3)

式中:I(σ)為一個(gè)指示函數(shù),定義為

(4)

2.4 支持向量分位數(shù)回歸

(5)

(6)

s.t.yi-wTφ(xi)-b≤ξi

Shim等[25]采用半?yún)?shù)方法通過(guò)構(gòu)造拉格朗日函數(shù)的方法對(duì)式(6)優(yōu)化問(wèn)題進(jìn)行求解,即

(7)

(8)

式中:xs、xt為任意變量;σK為核函數(shù)的寬度.采用廣義近似交叉驗(yàn)證準(zhǔn)則[26]獲得最優(yōu)超參數(shù)組合,即

(9)

3 實(shí)驗(yàn)設(shè)計(jì)

3.1 實(shí)驗(yàn)設(shè)置

采用3個(gè)具有不同充電策略的公開(kāi)數(shù)據(jù)集來(lái)驗(yàn)證所提方法的有效性.分別為馬里蘭大學(xué)高級(jí)生命周期工程中心公開(kāi)數(shù)據(jù)集(簡(jiǎn)稱(chēng)CALCE)[27]、牛津電池老化數(shù)據(jù)集(簡(jiǎn)稱(chēng)Oxford)[28]和麻省理工-斯坦福-豐田研究中心電池?cái)?shù)據(jù)集(簡(jiǎn)稱(chēng)TRI)[29].3個(gè)數(shù)據(jù)集的基本信息如表1所示,更為詳細(xì)的信息可參閱文獻(xiàn)[27-29].圖5分別以每個(gè)數(shù)據(jù)集中的1個(gè)電池為例,展示3個(gè)數(shù)據(jù)集各自電池容量隨循環(huán)次數(shù)的變化情況.可見(jiàn)盡管電池容量都隨循環(huán)次數(shù)增加而下降,但不同數(shù)據(jù)集的變化情況不同.

表1 3個(gè)數(shù)據(jù)集的基本信息Tab.1 Basic information of three datasets

圖5 CALCE、Oxford、TRI數(shù)據(jù)集電池容量變化Fig.5 Battery capacity change of CALCE, Oxford and TRI datasets

在每個(gè)數(shù)據(jù)集中選用1個(gè)放電程度較深的單體電池充放電數(shù)據(jù)作為訓(xùn)練集,在其余單體電池中任選3個(gè)作為測(cè)試集來(lái)驗(yàn)證所提方法的性能.所選訓(xùn)練集和測(cè)試集的情況如表2所示.

表2 訓(xùn)練集和測(cè)試集情況Tab.2 Basic information of training set and test set

3.2 參與比較的方法

所提方法的優(yōu)勢(shì)在于兩個(gè)方面:①可以適用于不同的充電策略;②采用的SVQR方法可以在小樣本數(shù)據(jù)集上更好地進(jìn)行SOH區(qū)間估計(jì).第一個(gè)優(yōu)勢(shì)顯而易見(jiàn),如文獻(xiàn)[13-14]的方法都針對(duì)某一特定充電策略,在遇到不同充電策略時(shí)所提取的特征參數(shù)不再適用.為此,在實(shí)驗(yàn)部分重點(diǎn)論證所采用的SVQR相比于其他區(qū)間預(yù)測(cè)方法的優(yōu)越性.首先對(duì)比SVQR與QR、QRNN,QR把分位數(shù)的概念融入普通的線性回歸中,本質(zhì)上是一種線性方法;QRNN將分位數(shù)概念與神經(jīng)網(wǎng)絡(luò)結(jié)合以解決非線性回歸問(wèn)題.此外,最近GPR被報(bào)道在處理鋰電池SOH區(qū)間估計(jì)方面具有優(yōu)秀性能[15],因此還將SVQR方法與GPR方法進(jìn)行對(duì)比.

3.3 評(píng)價(jià)指標(biāo)

為全面驗(yàn)證所提方法的性能,采用如下3個(gè)指標(biāo)來(lái)評(píng)價(jià)SOH區(qū)間估計(jì)結(jié)果的優(yōu)劣:平均區(qū)間分?jǐn)?shù)(Average Interval Score, AIS)[30]、平均預(yù)測(cè)區(qū)間中心偏差(Mean Prediction Interval Center Deviation, MPICD)[31]以及平均絕對(duì)百分比誤差(Mean Absolute Percentage Error,MAPE)[32].上述指標(biāo)包含對(duì)估計(jì)結(jié)果區(qū)間寬度、區(qū)間覆蓋率以及估計(jì)偏差等信息的衡量,各自定義如下:

(1) AIS指標(biāo).AIS指標(biāo)定義為所有區(qū)間估計(jì)點(diǎn)的區(qū)間分?jǐn)?shù)均值,可綜合衡量區(qū)間覆蓋率和區(qū)間寬度,其值越大越優(yōu),按下式計(jì)算.

(10)

式中:α為置信區(qū)間寬度;Ntest為測(cè)試集電池循環(huán)次數(shù);S(α)(xi)為第i個(gè)估計(jì)區(qū)間的區(qū)間分?jǐn)?shù),定義為

S(α)(xi)=

(11)

(2) MPICD指標(biāo).MPICD表示預(yù)測(cè)區(qū)間中心位置與實(shí)際值的偏差,該值越小越優(yōu),按下式計(jì)算.

(12)

(3) MAPE指標(biāo).MAPE按下式計(jì)算,該指標(biāo)可衡量預(yù)測(cè)區(qū)間的寬度信息,狹窄的區(qū)間比寬大的區(qū)間更好.

(13)

式中:Li、Ui分別為預(yù)測(cè)下限和上限.

4 實(shí)驗(yàn)結(jié)果與分析

經(jīng)過(guò)特征提取和交叉驗(yàn)證特征自動(dòng)優(yōu)選過(guò)程,最終對(duì)CALCE數(shù)據(jù)集選擇IC曲線峰值、IC曲線峰值對(duì)應(yīng)電壓、IC曲線左斜率、CC充電上升到給定電壓所用時(shí)間、CV充電下降到給定電流所用時(shí)間5個(gè)特征;對(duì)Oxford數(shù)據(jù)集選擇IC曲線峰值、IC曲線左斜率、IC曲線右斜率3個(gè)特征;對(duì)TRI數(shù)據(jù)集選擇IC曲線峰值、IC曲線右斜率、CV充電下降到給定電流所用時(shí)間3個(gè)特征.為了實(shí)驗(yàn)對(duì)比的公平性,其余3種方法在評(píng)估時(shí)針對(duì)不同的數(shù)據(jù)集采用與SVQR方法同樣的特征參數(shù).表3～5分別展示所提方法與QR、QRNN和GPR方法分別在CALCE、Oxford和TRI數(shù)據(jù)集中各自任選的3個(gè)單體電池上的SOH評(píng)估結(jié)果,其中用粗體標(biāo)出的為最優(yōu)結(jié)果.圖6～8分別展示4種方法在3個(gè)數(shù)據(jù)集中CS35、Cell7、b0c32這3個(gè)電池上的SOH估計(jì)效果.

首先,從圖6～8和表3～5可以看出,在所有數(shù)據(jù)集上,對(duì)于所有的評(píng)價(jià)指標(biāo)QRNN與QR結(jié)果比較接近,但都遠(yuǎn)劣于SVQR;QRNN略?xún)?yōu)于QR,這印證了QR本質(zhì)上是線性回歸方法,其在處理鋰電池SOH非線性估計(jì)上性能不足.而QRNN性能遠(yuǎn)劣于SVQR,這充分印證了在本文數(shù)據(jù)集的量級(jí)上,SVQR處理小樣本非線性問(wèn)題時(shí)具有更好的泛化能力.其次,在所有實(shí)驗(yàn)中SVQR方法的結(jié)果與GPR接近,但在多數(shù)電池上多數(shù)指標(biāo)都優(yōu)于GPR方法.尤其在MAPE指標(biāo)方面,SVQR方法在所有數(shù)據(jù)集上全優(yōu)于GPR.AIS指標(biāo)方面,SVQR方法在Oxford和TRI數(shù)據(jù)集上均取得最好結(jié)果,僅在CALCE數(shù)據(jù)集上,劣于GPR方法;MPICD指標(biāo)方面,SVQR方法在CALCE和Oxford數(shù)據(jù)集上除Cell7電池外,均明顯優(yōu)于GPR,僅在TRI數(shù)據(jù)集上略劣于GPR.文獻(xiàn)[15]報(bào)道,GRP方法是目前處理鋰電池SOH估計(jì)的最優(yōu)方法.但實(shí)際上,GPR需要假設(shè)數(shù)據(jù)服從高斯分布,而實(shí)測(cè)的鋰電池充放電數(shù)據(jù)涉及復(fù)雜過(guò)程和噪聲干擾,很難確保滿(mǎn)足高斯假設(shè).本文所提SVQR則無(wú)需此假設(shè),這可以解釋上述實(shí)驗(yàn)中SVQR總體優(yōu)于GPR的原因.另外值得指出的是,實(shí)驗(yàn)采用同一充電協(xié)議下的某一放電深度較深電池的全生命期循環(huán)充放電數(shù)據(jù)來(lái)訓(xùn)練SVQR模型,用同數(shù)據(jù)集中其他電池來(lái)測(cè)試,這與其他多數(shù)研究采用同一電池的部分充放電數(shù)據(jù)來(lái)預(yù)測(cè)其剩余循環(huán)的退化情況相比更有實(shí)際意義.上述結(jié)論充分證明本文方法在SOH估計(jì)方面的優(yōu)勢(shì).

圖6 4種方法在CALCE數(shù)據(jù)集CS35電池的SOH區(qū)間估計(jì)對(duì)比結(jié)果Fig.6 Comparison of four methods for estimating SOH interval of CS35 battery in CALCE dataset

圖7 4種方法在Oxford數(shù)據(jù)集Cell7電池的SOH區(qū)間估計(jì)對(duì)比結(jié)果Fig.7 Comparison of four methods for SOH interval estimation of Cell7 battery in Oxford dataset

圖8 4種方法在TRI數(shù)據(jù)集b0c32電池的SOH區(qū)間估計(jì)對(duì)比結(jié)果Fig.8 Comparison of four methods in estimating SOH interval of b0c32 battery in TRI dataset

表3 CALCE數(shù)據(jù)集3個(gè)單體電池的SOH估計(jì)對(duì)比結(jié)果Tab.3 Comparison of SOH estimation of three single cells in CALCE dataset 10-3

5 結(jié)論

提出一種考慮不同充電策略的鋰電池剩余容量SVQR區(qū)間估計(jì)方法.該方法在離線階段,從鋰電池充放電電壓、電流及容量增量曲線中獲取多個(gè)特征,并采用交叉驗(yàn)證方法自動(dòng)選擇最優(yōu)特征組合.選取放電程度較深的單體電池循環(huán)充放電數(shù)據(jù)訓(xùn)練SVQR模型.在SOH在線估計(jì)階段,采用和訓(xùn)練模型相同的特征作為模型輸入,對(duì)服役中的鋰電池進(jìn)行SOH區(qū)間估計(jì).最后在3個(gè)采用不同充電策略的數(shù)據(jù)集驗(yàn)證所提方法,并與QR、QRNN和GPR方法進(jìn)行對(duì)比,可得出以下結(jié)論:

(1) 在不同充電策略下,利用相同外部特征對(duì)于鋰電池SOH的估計(jì)結(jié)果也不同.采用兩步特征提取策略,即首先針對(duì)不同充電策略提取足夠多特征,再選取針對(duì)不同電池及充電策略的個(gè)性化特征參數(shù),可以顯著提高電池SOH估計(jì)的精度,且具有良好的普適性.

(2) 神經(jīng)網(wǎng)絡(luò)類(lèi)方法的訓(xùn)練需要較大數(shù)據(jù)集樣本,但在樣本量有限的鋰電池充放電數(shù)據(jù)集中神經(jīng)網(wǎng)絡(luò)類(lèi)算法估計(jì)效果有限.而SVQR利用SVR算法適用于小樣本數(shù)據(jù)集的特點(diǎn),即使是在小樣本條件下也可以較好地在置信區(qū)間內(nèi)估計(jì)鋰電池SOH.

(3) 相較于文獻(xiàn)報(bào)告的當(dāng)前最優(yōu)的GPR方法,所提SVQR方法無(wú)需原始數(shù)據(jù)服從高斯分布,可在一定置信區(qū)間內(nèi)實(shí)現(xiàn)較好的SOH估計(jì)結(jié)果,具有更好的適應(yīng)性和魯棒性.