王德順,薛金花,馬林康
(1.中國(guó)電力科學(xué)研究院有限公司,江蘇南京 210003;2.東南大學(xué) 自動(dòng)化學(xué)院,江蘇 南京 210096;3.北京交通大學(xué)國(guó)家能源主動(dòng)配電網(wǎng)技術(shù)研發(fā)中心,北京 100044)
儲(chǔ)能電池是電化學(xué)儲(chǔ)能電站的核心,儲(chǔ)能電站建設(shè)成本主要源于儲(chǔ)能電池系統(tǒng)。鋰離子電池在循環(huán)過程中的性能衰退情況受循環(huán)倍率、環(huán)境溫度、循環(huán)充放電區(qū)間等因素的影響[1]。電池健康狀態(tài)的劣化原因主要有電池內(nèi)阻增大、電池內(nèi)部不可逆的鋰離子損失、電池的正負(fù)極材料損失等[2]。目前對(duì)電池性能衰退的研究還是主要集中在電動(dòng)汽車應(yīng)用場(chǎng)景,鮮有針對(duì)儲(chǔ)能電站工況的電池衰退研究,對(duì)儲(chǔ)能系統(tǒng)電池組的衰退及不一致性演變研究也比較少。
柳楊等[3]基于電池電壓曲線隨老化會(huì)發(fā)生平移、伸縮變換的假設(shè),采用DTW(動(dòng)態(tài)時(shí)間扭曲算法)對(duì)電池組內(nèi)部單體的容量差異進(jìn)行估計(jì),用來描述電池組內(nèi)單體容量的不一致性。李建林等[4]認(rèn)為電池能量狀態(tài)主要通過電池SOC和電池的SOH來表征,并總結(jié)了常用的電池SOC和SOH估計(jì)方法。董漢成等[5]針對(duì)鋰離子電池SOH估計(jì)和剩余可用壽命(RUL)預(yù)測(cè)的問題,應(yīng)用支持向量回歸機(jī)(SVR-PF)對(duì)電池的健康狀態(tài)進(jìn)行估計(jì)同時(shí)實(shí)現(xiàn)電池RUL 的預(yù)測(cè)。張偉等[6]針對(duì)電池容量估計(jì)的問題,提出以平均充電壓升作為電池健康指標(biāo),經(jīng)Box-Cox變換后通過線性回歸模型估算電池容量。鄒涵宇[7]通過電池充電曲線的相似性對(duì)電池模組的一致性進(jìn)行辨識(shí)和分析,通過控制變量研究不同的衰退因素對(duì)電池組不一致性的影響,并利用相似性矩陣進(jìn)行電池分選。在電池的衰退機(jī)理研究方面,王永紅等[8]研究了三元電池衰退的主要原因,包括晶體結(jié)構(gòu)的改變和相變、活性材料的損失、電解質(zhì)的分解和消耗、可脫嵌鋰離子的損耗以及固體電解質(zhì)界面的形成。
上述方法大多是對(duì)電池單體進(jìn)行容量估算或壽命預(yù)測(cè),缺乏對(duì)電池組一致性的研究。本文以某型號(hào)磷酸鐵鋰電池模組為研究對(duì)象,設(shè)計(jì)典型工況模擬實(shí)驗(yàn)及電池性能測(cè)試來探究電池模組及單體衰退過程。通過電池循環(huán)前后容量增量曲線變化分析單體電池老化主要因素,同時(shí)提出一種串聯(lián)電池單體容量估計(jì)方法,進(jìn)一步研究了電池組隨工況循環(huán)SOC一致性演變過程,為電池組的維護(hù)提供依據(jù)。
實(shí)驗(yàn)采用的磷酸鐵鋰電池組結(jié)構(gòu)為1 并8 串,電池基本參數(shù)如表1 所示。
表1 電池組基本參數(shù)
在電池組進(jìn)行循環(huán)測(cè)試前,首先進(jìn)行電池性能測(cè)試,具體包括模組容量測(cè)試,模組內(nèi)阻測(cè)試,模組小電流測(cè)試以及單體電池容量測(cè)試。其中,小電流測(cè)試采用0.1C(15 A)恒流充放電以獲取電池組當(dāng)前最大可用容量。電池組容量測(cè)試多采用多階段降電流充電,1/3C恒流放電的充放電制式,以電池組放電的容量作為電池組的實(shí)際容量。其中,1#模組容量為159 Ah,2#模組容量為169 Ah。
模組內(nèi)阻測(cè)試采用復(fù)合脈沖功率特性(HPPC)測(cè)試方法,因?yàn)榱姿徼F鋰電池在50%SOC附近的內(nèi)阻比較穩(wěn)定[9],所有本文提到的內(nèi)阻均為電池模組SOC為50%時(shí)對(duì)應(yīng)的單體電池的0.1 s 的歐姆內(nèi)阻,兩組電池共計(jì)16 只單體的歐姆內(nèi)阻辨識(shí)結(jié)果如圖1 所示。
圖1 電池單體內(nèi)阻分布
由圖1 可以發(fā)現(xiàn),在測(cè)試的兩個(gè)電池模組中單體之間內(nèi)阻差異較大,大部分單體內(nèi)阻分布在0.2~0.3 mΩ 之間,少部分單體內(nèi)阻較大,分布在0.4~0.5 mΩ 的范圍,單體最大內(nèi)阻達(dá)到了最小內(nèi)阻的2倍及以上,且1#模組內(nèi)阻一致性相對(duì)較差。
電池單體的小電流測(cè)試采用0.1C,對(duì)模組內(nèi)每一只單體電池全區(qū)間充放電2~3 次,以最后一次放電容量作為單體電池的最大可用容量,兩個(gè)模組內(nèi)各單體容量測(cè)試結(jié)果如圖2所示。兩個(gè)模組內(nèi)單體容量差異明顯,其中1#模組單體容量均值177 Ah,極差為34 Ah,2#模組單體容量均值也為177 Ah,極差為26 Ah,綜合來看1#模組容量一致性劣于2#模組。
圖2 電池單體容量分布
電池組容量由電池組內(nèi)單體容量及單體SOC共同決定,文獻(xiàn)[10]指出,串聯(lián)電池組容量可以表示為電池組中單體電池最小的放電容量和最小的充電容量之和,即:
式中:SOCi為電池組內(nèi)第i個(gè)單體的SOC;Ci為第i個(gè)單體容量;SOCj為電池組內(nèi)第j個(gè)單體的SOC;Cj為第j個(gè)單體容量;Cpack為電池模組容量。
為了研究電池組在能量型工況下的衰退特性,本研究設(shè)計(jì)的模擬工況類型為典型儲(chǔ)能電站峰谷調(diào)節(jié)工況,考慮到電池的安全和使用壽命,設(shè)定電池循環(huán)SOC區(qū)間為10%~90%,充電倍率為3/8C,放電倍率為1/2C(每次性能測(cè)試后對(duì)模組重新進(jìn)行容量標(biāo)定,保證SOC循環(huán)區(qū)間不發(fā)生變化)。
電池組每100 次工況循環(huán)進(jìn)行一次性能測(cè)試,每500 次模組循環(huán)進(jìn)行一次小電流充放電測(cè)試。由于1#模組在第600~700 次循環(huán)過程中發(fā)生了較為嚴(yán)重的鼓包現(xiàn)象,且電池容量保持率已小于80%,在電池組完成700 次循環(huán)測(cè)試后僅進(jìn)行一次性能測(cè)試,不再進(jìn)行工況循環(huán)測(cè)試。
將1#模組和2#模組按照1.2 節(jié)中的電池循環(huán)測(cè)試流程開展工況循環(huán)測(cè)試,測(cè)試設(shè)備采用Arbin BT2000 電池組測(cè)試系統(tǒng),電池組前700 次循環(huán)對(duì)應(yīng)的電池組容量衰退曲線如圖3所示。
圖3 電池組工況循環(huán)容量衰退情況
由電池模組的容量衰退曲線可知,兩組電池在循環(huán)初期出現(xiàn)了較快的容量衰減,隨后衰退速度變緩,衰退速率基本保持不變,并且一致性較差的模組衰退速度要快于一致性較好的電池模組。在500 次后的循環(huán)測(cè)試中,1#模組容量衰退速度基本不變,2#模組衰退速度變得更慢。1#模組經(jīng)歷700次循環(huán)后容量保持率為75%,對(duì)應(yīng)的模組容量為120 Ah,2#模組經(jīng)歷700 次循環(huán)后其容量保持率為84%,對(duì)應(yīng)的模組容量為142 Ah。
從電池容量衰退過程來看,電池組在相同的工況下衰退速度受電池組初始容量、內(nèi)阻分布一致性和循環(huán)過程中單體衰退差異的綜合因素影響。在電池工況循環(huán)測(cè)試之前,電池組進(jìn)行了均衡(SOC高端對(duì)齊),因此,電池組內(nèi)單體容量及其衰退不一致可能是電池組衰退出現(xiàn)差異的主要原因。
容量增量分析法(ICA)作為一種原位非破壞的分析方法,常用于對(duì)鋰離子電池的老化分析和診斷[11]。IC(容量增量)曲線可以將電池原始的V-Q曲線上的電壓平臺(tái)轉(zhuǎn)換為易于分析和識(shí)別的氧化還原峰。磷酸鐵鋰電池充電V-Q曲線和IC 曲線對(duì)應(yīng)關(guān)系如圖4 所示,磷酸鐵鋰電池負(fù)極在充電過程中會(huì)發(fā)生五個(gè)相變過程,對(duì)應(yīng)電池V-Q曲線上的五個(gè)電壓平臺(tái),即IC 曲線上的五個(gè)容量增量峰,分別為⑤峰~①峰。
圖4 磷酸鐵鋰電池IC和V-Q 曲線關(guān)系圖
在第500 次和第700 次循環(huán)結(jié)束后對(duì)電池單體進(jìn)行了小電流充放電測(cè)試,通過數(shù)值的方法繪制出磷酸鐵鋰電池的IC曲線,電池循環(huán)前后IC 曲線變化如圖5 所示,可以發(fā)現(xiàn)幾乎所有的電池IC 曲線的變化都表現(xiàn)出類似的特征。電池IC 峰上最為明顯的是①峰的衰退,其次為②峰衰退,其余峰衰退不明顯,大多數(shù)電池IC 峰位置出現(xiàn)向高電壓方向的偏移。除此之外,電池循環(huán)前后電池IC 曲線的高端(V>3.4 V)和低端(V<3.2 V)基本重合,即電池高端的充電電量基本不變,近似為一定值。對(duì)于磷酸鐵鋰電池而言,電池不同的衰退模式對(duì)應(yīng)著不同IC 變化特征[12]。在峰谷調(diào)節(jié)工況下,由于電池①峰衰退最為明顯,可以推測(cè)該批次電池衰退的主要方式為鋰離子損失或電池正極材料損失;②峰雖然也出現(xiàn)明顯衰退,但與①峰相比其減小比例較小,即電池發(fā)生了少量的負(fù)極材料損失。同時(shí)電池IC 峰的位置出現(xiàn)了偏移,這是由于電池在循環(huán)過程中,內(nèi)阻增大導(dǎo)致充電過程極化程度增加,從而導(dǎo)致IC 峰和谷的位置向高電壓方向發(fā)生偏移。
圖5 循環(huán)前后電池單體IC曲線變化規(guī)律
由于模組內(nèi)部電池單體老化存在差異,其衰退速度必定不同,為了研究單體電池容量變化規(guī)律,有必要對(duì)單體電池的容量進(jìn)行估計(jì)。由于電池組最大可用容量是由部分單體決定的,除了組內(nèi)一只或者兩只單體會(huì)達(dá)到充電或者放電截止電壓,其他電池均不能實(shí)現(xiàn)滿充滿放,因此,如何基于電池模組測(cè)試數(shù)據(jù)來估計(jì)內(nèi)部單體容量成為要解決的主要問題。在電池模組小電流測(cè)試過程中,大部分電池只有部分區(qū)間充放電段,相應(yīng)地也只能得到部分IC 曲線,但電池衰退主要表現(xiàn)在V-Q曲線上平臺(tái)區(qū)的縮短,即IC 峰面積的減小,因此可以基于模組測(cè)試數(shù)據(jù)和單體電池初始測(cè)試數(shù)據(jù)實(shí)現(xiàn)電池單體容量估計(jì)。
由電池老化分析結(jié)果可知,電池高端部分(一般取3.4 V以上部分)充電容量近似為一定值,且電池低電壓區(qū)間內(nèi)⑤峰和④~③峰面積近似不變,因此可以利用數(shù)值方法擬合電池高端充電曲線,得到電池高端部分V-Q函數(shù)關(guān)系:
式中:Q為電池高端部分待充電電量;V為電池的充電電壓。當(dāng)電池模組進(jìn)行小電流充放電測(cè)試之后,可以得到某只電池組單體電池在電池組充滿電時(shí)的電壓記為Vt,進(jìn)而根據(jù)該單體末端充電容量擬合函數(shù)[式(2)]得到該只電池的末端充電容量Qhigh,其計(jì)算方法如式(3)所示:
式中:Vcut_off為單體電池充電截止電壓。
通過上述方法,可以得到電池高端部分完整的充電電量。對(duì)于單體的低端部分,同樣可以采用類似的方法估計(jì)低端電量,最終,電池的總?cè)萘烤涂梢员硎緸椋?/p>
式中:Qmodule5_1為模組充電數(shù)據(jù)中各個(gè)單體平臺(tái)區(qū)充電電量(特征峰面積之和);Qlow為電池低端部分電池充電電量。
通過上述方法,可以在獲取某只單體電池起始完整的V-Q曲線后,僅利用模組的性能測(cè)試數(shù)據(jù),就可以對(duì)電池組內(nèi)單體電池容量進(jìn)行估計(jì),節(jié)省測(cè)試時(shí)間,實(shí)現(xiàn)每一只單體電池的容量快速估計(jì)。
利用上述容量估計(jì)方法,每進(jìn)行一次模組測(cè)試數(shù)據(jù)就可以對(duì)單體電池容量進(jìn)行估計(jì),由于電池組循環(huán)過程中僅進(jìn)行了兩次單體性能測(cè)試,取第二次單體性能測(cè)試數(shù)據(jù)為驗(yàn)證數(shù)據(jù),得到兩個(gè)模組第500 次循環(huán)結(jié)束后的單體電池容量估計(jì)結(jié)果及估計(jì)誤差,估計(jì)結(jié)果如圖6 所示。
由圖6 可以看出,1#模組中7 號(hào)單體電池容量估算誤差較大,相對(duì)誤差達(dá)到為5.9%。除此之外,其余單體電池容量估計(jì)誤差均在4%以內(nèi),表明該方法可以實(shí)現(xiàn)單體電池容量快速和可靠估計(jì),且該方法不依賴于大量歷史數(shù)據(jù)的訓(xùn)練,在獲取單體電池初始數(shù)據(jù)和模組充電數(shù)據(jù)后就可以實(shí)現(xiàn)單體電池容量估計(jì)。
圖6 模組內(nèi)單體電池容量估計(jì)結(jié)果
電池組在工況循環(huán)過程中單體電池的SOC會(huì)發(fā)生不同程度變化,進(jìn)而影響電池組實(shí)際最大可用容量,甚至使得電池組實(shí)際最大可用容量小于模組中最小容量單體,電池容量不能得到有效利用。因此,有必要對(duì)電池組內(nèi)單體電池的SOC變化趨勢(shì)進(jìn)行分析,在獲取某只電池單體容量估計(jì)值之后,其高低端SOC可以由式(5)獲?。?/p>
式中:SOC_H為電池高端SOC值。由上述單體SOC估算方法,可以獲得每100 次循環(huán)后的單體電池高端SOC估計(jì)值,進(jìn)而得到前500 次循環(huán)模組內(nèi)部各單體高端SOC隨模組循環(huán)過程中的變化情況,如圖7 所示。
圖7 模組內(nèi)電池單體高端SOC變化趨勢(shì)
由圖7 可知,兩個(gè)模組內(nèi)大部分單體電池在500 次循環(huán)之后單體SOC仍保持在較高水平。然而,1#模組的一號(hào)單體電池及2#模組五號(hào)單體電池SOC出現(xiàn)了較為明顯的降低,其中2#模組五號(hào)單體經(jīng)過500 次模組循環(huán)之后高端SOC減小為88%,電池組內(nèi)部的SOC極差達(dá)到了11.7%。通過之前的分析可知,2#模組五號(hào)單體電池及1#模組一號(hào)單體電池在模組中具有最高的單體容量,與模組內(nèi)部其他容量較低的電池相比,其衰退速度也較慢,同時(shí)電池間的充放電效率也存在一定差異,二者相互疊加導(dǎo)致電池組內(nèi)單體SOC變化的差異。
為了探究磷酸鐵鋰電池組在儲(chǔ)能系統(tǒng)應(yīng)用中的衰退規(guī)律,開展了電池組模擬工況循環(huán)老化研究。選取兩個(gè)磷酸鐵鋰電池模組,設(shè)計(jì)了模組電池峰谷調(diào)節(jié)工況的循環(huán)老化實(shí)驗(yàn)。700 次循環(huán)之后電池組容量保持率分別為初始容量的75%和84%,且一致性較差的模組容量衰退較快。通過容量增量分析法對(duì)單體電池循環(huán)前后進(jìn)行老化分析,發(fā)現(xiàn)電池容量衰退的主要原因是鋰離子損失和正極活性材料損失。為了辨別單體電池容量衰退差異,基于電池模組和單體性能測(cè)試數(shù)據(jù),提出了一種串聯(lián)電池單體的容量估計(jì)方法,并進(jìn)行了驗(yàn)證,大部分電池單體容量估計(jì)誤差在4%以內(nèi)?;陔姵貑误w容量估計(jì)結(jié)果,進(jìn)一步對(duì)單體電池的SOC變化進(jìn)行分析,結(jié)果表明模組內(nèi)容量較高的電池SOC變化相對(duì)于其他容量較低的單體電池變化更明顯,這種現(xiàn)象出現(xiàn)的原因是單體電池衰退不一致和單體電池充放電效率差異的共同作用。