王玉坤, 姜久春, 張彩萍, 張維戈, 馬澤宇
(北京交通大學(xué)電氣工程學(xué)院,北京 100044)
鋰離子電池是新型綠色環(huán)保蓄電池,與其他傳統(tǒng)蓄電池相比,它具有比能量高、大電流放電能力強(qiáng)、循環(huán)壽命長等特性,這些特性決定了其在新能源汽車、儲能電源及小型分布式電站方面將具有良好的應(yīng)用前景[1]。2008年北京奧運(yùn)會上投入了50輛純電動大客車,實(shí)現(xiàn)了電動汽車的規(guī)?;瘧?yīng)用,做到了奧林匹克公園核心區(qū)內(nèi)的“零”排放。經(jīng)過3年使用后,電池容量有不同程度衰退。根據(jù)國家電動汽車電池使用標(biāo)準(zhǔn),電動汽車用動力電池有效容量不得低于額定容量的80%。為降低整車成本,當(dāng)電池有效容量低于80%以后,可以將其用于儲能系統(tǒng)或一些小型的儲能裝置中進(jìn)行二次利用。
電池測試中,常用的容量測試方法有恒流恒壓測試方法[2],每塊電池約8小時(shí)充放電測試以得到準(zhǔn)確的容量值;另有一種根據(jù)電池開路電壓與剩余容量的對應(yīng)關(guān)系,容量增量與SOC增量之比的計(jì)算方法[3],同樣需要對電池進(jìn)行充放電實(shí)驗(yàn),但測得容量值準(zhǔn)確性沒有驗(yàn)證。車用淘汰電池?cái)?shù)目巨大,逐一檢測電池容量耗時(shí)、耗力、耗材,快速估算淘汰電池容量成為梯次利用的關(guān)鍵。
本文首先對一整車2008年奧運(yùn)會大巴車用淘汰電池容量和內(nèi)阻的分布特性進(jìn)行分析,研究表明兩者均符合近似正態(tài)分布且沒有線性相關(guān)性,基于不同電池不同采樣時(shí)間直流內(nèi)阻所建立的極化內(nèi)阻均不相同,建立不同采樣時(shí)間直流內(nèi)阻作為輸入,容量作為輸出的非線性BP神經(jīng)網(wǎng)絡(luò)模型。并抽取30個電池模塊不同采樣時(shí)間直流內(nèi)阻作為輸入進(jìn)行仿真,仿真結(jié)果與傳統(tǒng)容量測試結(jié)果基本吻合,實(shí)現(xiàn)了梯次利用電池容量的快速估計(jì),并且驗(yàn)證了該方法的準(zhǔn)確性。
本文以一整車10箱100節(jié)北京奧運(yùn)會純電動大客車淘汰的錳酸鋰能量型電池模塊為測試對象。單體電池額定容量90 Ah,電池模塊為4個單體電池并聯(lián),額定容量為360 Ah,電池模組由8個或者16個模塊串聯(lián)組成。電池模塊測試采用拜特電池充放電測試設(shè)備,電池模組測試采用Arbin充放電測試設(shè)備。實(shí)驗(yàn)測試環(huán)境溫度為(20±5)℃。
電池模塊容量采用恒流恒壓測試方法[2],為了得到準(zhǔn)確的容量值,連續(xù)測試5個循環(huán)以上,當(dāng)連續(xù)3次測試容量的極差與其平均值之比在3%以內(nèi),停止實(shí)驗(yàn),并將3次測試容量的平均值作為電池實(shí)際容量。
直流內(nèi)阻的測量采用復(fù)合脈沖電流法[4],通過對電池輸入電流階躍信號(△I)并測定對應(yīng)的電壓變化值(△U),利用歐姆定律得到直流內(nèi)阻R=△U/△I。測試方法是將電池組放電至50%荷電狀態(tài),靜置1 h,使電池極化退去,狀態(tài)趨于穩(wěn)定。電池首先以1/3C放電10 s,再以1/3C充電持續(xù)10 s,電壓響應(yīng)如圖1所示。在第30 s電流信號由-1/3C變?yōu)?/3C時(shí)電壓增量△U,以此計(jì)算直流內(nèi)阻。由于電池已經(jīng)處于老化狀態(tài),為了避免過充過放,加在電池上的電流不宜過大。這種方法電流由負(fù)值直接變?yōu)檎?,相?dāng)于電流幅值加倍,電壓響應(yīng)也同時(shí)加倍,減小設(shè)備采集誤差及極化對直流內(nèi)阻測量的影響。
圖1 復(fù)合脈沖電流及電壓響應(yīng)
由于電池在車中所處的位置不同,溫度、震動程度、連接情況的差異使電池容量、內(nèi)阻的衰退呈現(xiàn)出不一致性。梯次利用電池模塊容量和內(nèi)阻分布及擬合統(tǒng)計(jì)分布曲線分別如圖2、3所示。新電池模塊容量為360 Ah,極差小于1 Ah;1 s直流內(nèi)阻為0.325mΩ,極差小于0.01mΩ。使用三年后電池模塊容量均值200 Ah,容量衰退率約45%,電池容量極差30 Ah;1 s直流內(nèi)阻均值0.42mΩ,內(nèi)阻增加率約30%,極差0.216mΩ。容量衰退率分布較為集中,內(nèi)阻增加率分布較為分散。較新電池而言,梯次利用電池容量、內(nèi)阻一致性明顯下降。
對100只電池模塊的容量和內(nèi)阻數(shù)據(jù)做非參數(shù)檢驗(yàn)(K-S檢驗(yàn)及S-W檢驗(yàn)),檢驗(yàn)其是否服從正態(tài)分布[5]。當(dāng)樣本含量n≤2 000時(shí),結(jié)果以S-W檢驗(yàn)為準(zhǔn),以K-S檢驗(yàn)為參考,兩種檢驗(yàn)方式結(jié)果均表明容量和內(nèi)阻都服從正態(tài)分布,近似概率密度函數(shù)分別如式(1)和式(2)所示,容量近似服從N(200,34)分布,內(nèi)阻近似服從 N(0.42,0.003 85)的分布。
圖2 容量分布圖
圖3 內(nèi)阻分布圖
通過散點(diǎn)圖粗略地給出變量相關(guān)關(guān)系的關(guān)聯(lián)類型與程度,再通過統(tǒng)計(jì)分析,求出描述變量間相關(guān)程度與變化方向的量數(shù),即相關(guān)系數(shù)。相關(guān)系數(shù)絕對值越接近1,兩個變量之間的相關(guān)關(guān)系越密切。
已驗(yàn)證電池容量衰退率和內(nèi)阻增加率均服從正態(tài)分布,滿足Pearson相關(guān)系數(shù)的使用條件,故采用Pearson相關(guān)系數(shù)衡量兩者的相關(guān)性程度。取采樣時(shí)間1 s計(jì)算內(nèi)阻增加率,分析得出容量衰退率和內(nèi)阻增加率的Pearson相關(guān)系數(shù)為0.445,可知在置信度(雙側(cè))為0.01時(shí),容量衰退率和內(nèi)阻衰退率具有一定的相關(guān)性。電池容量衰退率與內(nèi)阻增加率散點(diǎn)圖如圖4所示,容量衰退率與內(nèi)阻增加率點(diǎn)均勻散落在圖中,同一容量衰退率對應(yīng)多個內(nèi)阻增加率,同一內(nèi)阻增加率對應(yīng)多個容量衰退率,兩者并沒有明顯的趨勢線,增勢不明顯。這說明容量的衰退會伴隨一定的內(nèi)阻增加,在電池的使用過程中內(nèi)阻增加和容量衰退是同時(shí)發(fā)生的,由于電池所處環(huán)境不同,內(nèi)阻增加和容量衰退速率不一致。取容量衰退率在0.45到0.48區(qū)間再次計(jì)算容量衰退率和內(nèi)阻增加率的Pearson相關(guān)系數(shù)只有0.066,容量衰退率和內(nèi)阻增加率幾乎是完全不相關(guān)的,所以不能通過內(nèi)阻和容量的線性關(guān)系估算電池容量。
圖4 容量衰退率與內(nèi)阻增加率散點(diǎn)圖
人工神經(jīng)網(wǎng)絡(luò)(ANN),簡稱神經(jīng)網(wǎng)絡(luò),它由大量人工神經(jīng)元互聯(lián)組成,通過調(diào)整內(nèi)部大量節(jié)點(diǎn)相互連接的關(guān)系,模仿腦神經(jīng)網(wǎng)絡(luò)行為,進(jìn)行分布式并行處理信息。通用的神經(jīng)元模型如圖5所示。x=(x1,x2,…,xn)T為神經(jīng)元輸入,w=(w1,w2,…,wn)T為可調(diào)的輸入權(quán)值,θ為偏移信號,用于建模神經(jīng)元的興奮閾值。U(·)和f(·)分別表示神經(jīng)元的基函數(shù)和激活函數(shù)?;瘮?shù)U(·)是一個多輸入單輸出函數(shù)u=u(x,w,θ);激活函數(shù)的一般作用是對基函數(shù)輸出u進(jìn)行“擠壓”:y=f(u),即通過非線性函數(shù)f(·)將u變換指定范圍內(nèi)。神經(jīng)元之間的連接可以有任意形式,最常見的結(jié)構(gòu)是前向神經(jīng)網(wǎng)絡(luò)和反饋神經(jīng)網(wǎng)絡(luò)[6]。
圖5 通用的神經(jīng)元模型
前向神經(jīng)網(wǎng)絡(luò)中的BP網(wǎng)絡(luò)由于其較強(qiáng)的適應(yīng)性和學(xué)習(xí)功能廣泛應(yīng)用于實(shí)際中。BP網(wǎng)絡(luò)是一種多層前饋神經(jīng)網(wǎng)絡(luò),由一個輸入層、一個輸出層和若干個隱含層所組成。位于同一層的單元之間沒有連接,各層的單元只能向高層的單元輸出激活信號。BP算法有兩部分組成:信息的正向傳遞與誤差的反向傳播。在正向傳播過程中,輸入信息從輸入經(jīng)隱含層逐層計(jì)算傳向輸出層,每一層神經(jīng)元的狀態(tài)只影響下一層神經(jīng)元的狀態(tài)。如果在輸出層沒有得到期望的輸出,則計(jì)算輸出層的誤差變化值,然后轉(zhuǎn)向反向傳播,通過網(wǎng)絡(luò)將誤差信號沿原來的連接通路反傳回來修改各層神經(jīng)元的權(quán)值,直至達(dá)到期望目標(biāo)。
錳酸鋰電池在0.1 s左右建立歐姆內(nèi)阻[5],極化內(nèi)阻隨時(shí)間開始建立,取1、5、10 s采樣時(shí)間表征電池不同時(shí)間極化建立程度。實(shí)驗(yàn)表明,0.1、1、5、10 s四個采樣時(shí)間的直流內(nèi)阻均服從正態(tài)分布,故采用Pearson相關(guān)系數(shù)來衡量容量與內(nèi)阻的相關(guān)性程度。不同采樣時(shí)間對應(yīng)的容量和直流內(nèi)阻的Pearson相關(guān)系數(shù)如表1所示。
表1 容量與內(nèi)阻的Pea rson相關(guān)系數(shù)
由表1可以看出,不同采樣時(shí)間內(nèi)阻與容量的Pearson相關(guān)系數(shù)并不相同。這說明不同程度衰退的電池在1、5、10 s建立的極化內(nèi)阻均不相同,不同采樣時(shí)間直流內(nèi)阻可以表征電池不同內(nèi)阻特性。采樣時(shí)間為5 s時(shí),容量和內(nèi)阻的負(fù)相關(guān)性最強(qiáng)。為全面表征電池內(nèi)阻-容量非線性關(guān)系,采用0.1、1、5、10 s四個采樣時(shí)間的直流內(nèi)阻作為神經(jīng)網(wǎng)絡(luò)模型的輸入?yún)?shù)。
電池模塊的直流內(nèi)阻隨SOC變化曲線如圖6所示。當(dāng)電池荷電狀態(tài)處于20%~80%區(qū)間內(nèi),不同采樣時(shí)間直流內(nèi)阻隨SOC變化較小[3],所以在建立神經(jīng)網(wǎng)絡(luò)模型或者利用模型進(jìn)行大規(guī)模電池容量快速估計(jì)時(shí),只需要電池的荷電狀態(tài)不處于兩段狀態(tài),即可將電池模塊串聯(lián)成組,利用復(fù)合脈沖電流法進(jìn)行內(nèi)阻測試。
圖6 不同采樣時(shí)間直流內(nèi)阻隨SOC變化曲線
由Kolmogorov定理和BP定理[7]可知,對于一個三層BP神經(jīng)網(wǎng)絡(luò),只要隱含層節(jié)點(diǎn)數(shù)足夠多,就具有逼近任意復(fù)雜的非線性映射的能力。通過上文中容量和內(nèi)阻測試及其特性研究,本文選擇已測試的70個電池模塊的采樣時(shí)間為0.1、1、5、10 s的內(nèi)阻值作為模型輸入,容量值作為模型輸出。學(xué)習(xí)規(guī)則采用有導(dǎo)制最速下降BP算法,權(quán)值和閾值的修正是在所有樣本輸入之后,計(jì)算其總的誤差后進(jìn)行的,在樣本數(shù)比較多時(shí),這種批處理方式比分別處理方式的收斂速度快;隱含層中的神經(jīng)元均采用S型傳輸函數(shù);由于模型用于函數(shù)逼近,所以輸出層神經(jīng)元采用purelin線性傳輸函數(shù)。首先建立3層神經(jīng)網(wǎng)絡(luò),即1個隱含層的神經(jīng)網(wǎng)絡(luò),輸入層節(jié)點(diǎn)數(shù)為4,輸出層節(jié)點(diǎn)數(shù)為1,根據(jù)經(jīng)驗(yàn)公式(3)計(jì)算[8],預(yù)設(shè)隱含層節(jié)點(diǎn)數(shù)為13。
式中:n為輸入層節(jié)點(diǎn)數(shù);m為輸出層節(jié)點(diǎn)數(shù);a為l~10之間的常數(shù)。
理論上增加模型隱含層層數(shù)及隱含層節(jié)點(diǎn)數(shù)提高預(yù)測精度,研究表明對比本模型增加隱含層層數(shù)及節(jié)點(diǎn)數(shù)并沒有提高模型精度,最終確定1個隱含層,14個隱含層節(jié)點(diǎn)數(shù)模型預(yù)測精度最高,建立模型如圖7所示。
圖7 BP神經(jīng)網(wǎng)絡(luò)模型
模型精度如圖8所示,實(shí)際容量與估計(jì)容量平均誤差1.2 Ah(實(shí)際容量的0.6%),最大誤差7 Ah(實(shí)際容量的3.5%),模型精度較高。
圖8 模型估計(jì)值與實(shí)驗(yàn)值對比
應(yīng)用上述建立的BP神經(jīng)網(wǎng)絡(luò),取具有相同的行駛工況的30個電池模塊進(jìn)行預(yù)測驗(yàn)證。估計(jì)容量與實(shí)際測量容量如圖9所示。
圖9 模型估計(jì)圖
實(shí)際容量與估計(jì)容量誤差在6 Ah(實(shí)際容量的3%)以內(nèi),占抽取電池總數(shù)的80%,最大誤差9 Ah(實(shí)際容量的4.5%),模型預(yù)測精度較高。對比傳統(tǒng)逐一檢測大量淘汰電池容量的方法,神經(jīng)網(wǎng)絡(luò)法只需要檢測一車電池的容量和內(nèi)阻值建立BP網(wǎng)絡(luò)模型,通過脈沖法測試其他車電池的不同采樣時(shí)間直流內(nèi)阻,基于R-Q(內(nèi)阻-容量)神經(jīng)網(wǎng)絡(luò)模型可以實(shí)現(xiàn)快速準(zhǔn)確估計(jì)電池的容量。
目前北京市有大量2008年奧運(yùn)會電動汽車淘汰的鋰離子電池。本論文以此為應(yīng)用背景,研究這些淘汰電池的快速容量估計(jì)方法,對有效挖掘電池剩余價(jià)值、提高成組后電池組的使用壽命,大幅降低電動汽車總體擁有成本有著重要意義。本文得出如下結(jié)論:
(1)容量衰退率和內(nèi)阻增加率不存在線性相關(guān)關(guān)系;
(2)以行駛工況的電動汽車淘汰電池中任意一車?yán)肂P神經(jīng)網(wǎng)絡(luò)算法,以20%~80%荷電狀態(tài)的0.1、1、5、10 s采樣時(shí)間直流內(nèi)阻作為輸入?yún)?shù),容量作為輸出信號,建立非線性估計(jì)模型,仿真表明,電池容量估計(jì)精度在5%以內(nèi)。因此,其他車淘汰電池均可利用該車電池所建立的神經(jīng)網(wǎng)絡(luò)模型,整箱檢測20%~80%荷電狀態(tài)時(shí)0.1、1、5、10 s采樣時(shí)間直流內(nèi)阻估計(jì)電池容量,實(shí)現(xiàn)車用淘汰電池容量快速準(zhǔn)確估計(jì)。
[1]沈曉彥,黃鐘琪,周建新,等.鋰電池在風(fēng)光發(fā)電儲能系統(tǒng)中的應(yīng)用分析[J].電源技術(shù),2011,35(5):602-604.
[2]國家發(fā)展和改革委員會.電動汽車用鋰離子蓄電池測試規(guī)范[S].北京:國家發(fā)展和改革委員會,2006.
[3]文鋒.純電動汽車用鋰離子電池組管理技術(shù)基礎(chǔ)問題研究[D].北京:北京交通大學(xué),2010.
[4]SCHWEIGER HG,OBEIDIO.Comparison of severalmethods for determ ining the internal resistance of lithium ion cells[J].Sensors,2010(2):5604-5625.
[5]魯妍.梯次利用鋰離子動力電池特性及仿真研究[D].北京:北京交通大學(xué),2012.
[6]張德豐.MATLAB神經(jīng)網(wǎng)絡(luò)應(yīng)用設(shè)計(jì)[M].北京:機(jī)械工業(yè)出版社,2012.
[7]焦淑華,夏冰,徐海靜,等.BP神經(jīng)網(wǎng)絡(luò)預(yù)測的MATLAB實(shí)現(xiàn)[J].哈爾濱金融高等專科學(xué)校學(xué)報(bào),2009,97(1):55-56.
[8]葉斌,雷燕.關(guān)于BP網(wǎng)中隱含層層數(shù)及其節(jié)點(diǎn)數(shù)選取方法淺析[J].商丘職業(yè)技術(shù)學(xué)院學(xué)報(bào),2004,15(3):52-54.
[9]林志強(qiáng),孫水華.利用在線內(nèi)阻測試方法捕捉蓄電池組中的異常個體[J].貴州大學(xué)學(xué)報(bào),2010,27(6):74-80.
[10]郭宏榆,姜久春,王吉松,等.功率型鋰離子動力電池的內(nèi)阻特性[J].北京交通大學(xué)學(xué)報(bào),2011,35(5):119-123.
[11]李剛,謝永成,李光升,等.基于自適應(yīng)神經(jīng)網(wǎng)絡(luò)模糊推理系統(tǒng)的蓄電池SOH預(yù)測[J].技術(shù)與方法,2011,30(22):82-87.
[12]何洪文,余曉江.電動車輛動力電池的性能評價(jià)[J].吉林大學(xué)學(xué)報(bào):工學(xué)版,2006,36(5):659-663.
[13]劉險(xiǎn)峰,鄒積巖,李立偉.大容量蓄電池組的連接方式[J].吉林大學(xué)學(xué)報(bào):工學(xué)版,2007,37(3):672-674.
[14]高飛,李建玲,趙淑紅,等.鋰動力電池性能衰退規(guī)律研究進(jìn)展[J].中國稀土學(xué)報(bào),2010,28:819-823.
[15]邱綱,陳勇.電動汽車用動力電池組SOC的神經(jīng)網(wǎng)絡(luò)估計(jì)[J].遼寧工程技術(shù)大學(xué)學(xué)報(bào),2006,25(2):230-233.