基于改進粒子濾波算法的動力鋰離子電池荷電狀態(tài)估計

劉 淑 杰， 郝 昆 昆， 王 永， 鄧 威 威

( 大連理工大學 機械工程學院, 遼寧 大連 116024 )

0 引 言

當前環(huán)境污染和能源危機問題日益突顯，氣候變暖和溫室效應對全球環(huán)境也造成了巨大的影響，促進新型能源和儲能系統(tǒng)的研究是應對環(huán)境和能源危機一個較好的解決方案．在眾多儲能裝置中，鋰離子電池由于其應用廣泛，而且作為新能源電動汽車中的儲能裝置，被看作具有良好發(fā)展前景的一種儲能方式．與傳統(tǒng)類型的鉛酸、鎳氫電池相比較，鋰離子電池具有高能量密度和功率密度，循環(huán)使用壽命長以及自放電率極低、污染小等優(yōu)點[1]．電動汽車推廣使用對減少溫室氣體排放有著重要意義，電池包作為電動汽車的重要供能裝置，在汽車使用過程中的可靠性和安全性應予以保證，電池管理系統(tǒng)(battery management system，BMS)可以實現(xiàn)對電池包信息監(jiān)測和狀態(tài)預測功能，防止電池包過充電過放電行為發(fā)生，保證電池使用過程中性能安全可靠，使其最佳能效得到發(fā)揮并延長使用壽命[2]．

電池包在工作過程中可觀測的物理量僅有電壓、電流和溫度，而電池的一些關鍵指標如荷電狀態(tài)(state of charge，SOC)、功率狀態(tài)和健康狀態(tài)等信息則無法通過直接測量獲得，需要借助可觀測量和相應算法或模型來了解電池包的真實狀態(tài)．SOC基本定義是電池工作過程中剩余容量與電池額定容量百分比，是表征電池基本狀態(tài)信息的關鍵指標，其意義等同于燃油汽車中的燃油表[3]．大量文獻研究電池SOC估計方法，其大致可以分為安時法、開路電壓法、模型法、阻抗法以及基于模糊邏輯和機器學習的聯(lián)合估計方法．安時法雖然簡單易實施，但是需要準確獲取電池的初始容量和充放電電流，當充放電時間較長時，累計誤差對最終估計結(jié)果造成的影響便難以忽略，此方法需要進行定期校正，因此安時法通常與開路電壓法結(jié)合使用[4-5]．基于模型類的估計方法可分為電路模型法與電化學模型法[6-7]，電路模型法可以分為SOC直接估計方法與基于自適應濾波器及觀測器的改進估計方法[8-12]，直接進行SOC估計是一種開環(huán)方法，而基于濾波器和觀測器類的SOC估計方法是閉環(huán)方法，閉環(huán)方法能修正狀態(tài)估計過程中存在的偏差，對模型準確性以及測量信號精度要求不是非常嚴格；另一種模型方法是電化學模型方法，該方法能從本質(zhì)上反映電池內(nèi)部電化學行為與電池SOC及溫度之間的關系，無須電路模型方法中要借助查表方法來獲取不同SOC和溫度間的關系，但電化學模型方法的缺點是模型中含有大量偏微分方程，使算法計算量和模型復雜度要遠大于電路模型法．另一類常用的方法為模糊邏輯和機器學習算法，如模糊邏輯算法、神經(jīng)網(wǎng)絡算法、基于模糊邏輯的神經(jīng)網(wǎng)絡算法和支持向量機等機器學習類算法，這類方法的估計精度與模型訓練數(shù)據(jù)量及數(shù)據(jù)的準確性有較大關系[13-16]．

電動汽車中電池包SOC估計常采用等效電路模型方法[17-18]，等效電路模型根據(jù)電學原理用電子元器件近似模擬電池中的電化學物理量[19-20]，建立相應的電路方程，推導出電路方程與所估計狀態(tài)量間的數(shù)學關系，是一種數(shù)學模型且計算量遠小于電化學模型．借助相應的電路模型，將SOC整合到電路方程中，根據(jù)可測量的電學量估計OCV(open circuit voltage，OCV)，從而建立SOC與電池開路電壓OCV間關系．利用OCV對電池SOC進行估計，可靠的SOC-OCV映射關系和準確的OCV估算是電池SOC估計時需要考慮的兩個關鍵問題．Plett[20]應用擴展卡爾曼濾波(extended Kalman filter，EKF)方法提出了一種狀態(tài)和模型參數(shù)雙估計器，優(yōu)化的參數(shù)包括OCV、歐姆內(nèi)阻、庫侖效率和電池的容量等；He等[21]在Plett研究基礎上將遞歸最小二乘法與EKF結(jié)合應用于估計器中，對電池狀態(tài)和模型參數(shù)進行估計；El Mejdoubi等[22]假定OCV是一個隨時間緩慢變化的信號，建立電池動態(tài)數(shù)學模型，用EKF方法和李雅普諾夫觀測器進行OCV預測，但是此方法中模型觀測器參數(shù)很難調(diào)整，而EKF方法僅適用于高斯噪聲環(huán)境；Dang等[23-24]應用雙神經(jīng)網(wǎng)絡擴散模型基于OCV進行SOC估計，一個神經(jīng)網(wǎng)絡用于辨識電池模型參數(shù)，另一個神經(jīng)網(wǎng)絡用于建立電池OCV和SOC的函數(shù)映射關系，之后又基于電池電學特性在文獻[24]中采用受控自回歸滑動平均模型估計電池OCV，運用BP神經(jīng)網(wǎng)絡對SOC與OCV間的關系進行映射，但是此類方法在模型訓練階段需要大量數(shù)據(jù)；Chen等[25]采用粒子濾波(particle filter，PF)算法對一階等效電路模型中OCV進行估計，根據(jù)實驗中建立的SOC-OCV關系估計電池的SOC，雖然PF算法可用于求解非線性和非高斯噪聲條件下系統(tǒng)估計問題，但在迭代求解過程中存在粒子退化現(xiàn)象，當粒子退化嚴重時，最終的估算結(jié)果會產(chǎn)生較大誤差．針對PF算法進行SOC估計時存在的問題，本文分別采用改進的擴展粒子濾波(extended particle filter，EPF)和無跡粒子濾波(unscented particle filter，UPF)算法對電池SOC進行估計，首先建立電池模型并運用帶遺忘因子的遞歸最小二乘方法對模型中參數(shù)進行實時估計，將辨識后的參數(shù)代入系統(tǒng)方程中，分別以EKF和UKF作為EPF和UPF中的建議密度函數(shù)，對電池進行SOC估計．

1 電池建模和參數(shù)估計

1.1 電池建模

電路模型通常將一個或多個并聯(lián)的電容-電阻(RC)網(wǎng)絡進行串聯(lián)，模擬電壓的緩慢變化過程，串聯(lián)歐姆電阻模擬電壓的快速變化過程，Zhang等[26]對RC個數(shù)與模型計算復雜度和估計誤差進行了研究，結(jié)果表明隨RC個數(shù)增加估算的精度會提高，但同時算法的運算時間也會增加．本文綜合考慮模型精度與算法復雜度，采用一階RC電路模型作為電池模型[27]，如圖1所示．

圖1電路模型中包含一個歐姆內(nèi)阻和一個RC網(wǎng)絡，UOCV為電池開路電壓；Rin為電池歐姆內(nèi)阻，用于模擬充放電過程中輸出端電壓快速變化過程；Rp、Cp代表電池內(nèi)部的極化內(nèi)阻和極化電容，模擬工作過程中電池兩端電壓緩慢變化過程；Ut為鋰離子電池的輸出電壓；URin為電池歐姆內(nèi)阻兩端有電流通過時產(chǎn)生的電壓；Up為電池內(nèi)部并聯(lián)Rp、Cp網(wǎng)絡兩端的電壓，并規(guī)定放電時電流符號為負，充電時為正．根據(jù)電路學基本原理可得到電路模型的數(shù)學表達式如下：

(1)

Ut(t)=UOCV(t)-(URin(t)+Up(t))

(2)

根據(jù)SOC定義，其表達式可寫為

(3)

式中：η為電池庫侖效率，Qn代表電池額定容量，SOC(t0)表示t0時刻電池初始SOC值，I(t)表示電池工作過程中電流．

方程(1)進行整理后可得

(4)

方程(2)、(4)進行拉氏變換得到

Ut(s)-UOCV(s)=-(UR(s)+Up(s))=

(5)

(6)

令U(s)=I(s)，Y(s)=Ut(s)-UOCV(s)，得系統(tǒng)傳遞函數(shù)表達式：

(7)

(8)

式(8)中，T為離散后系統(tǒng)采樣步長，本研究中采樣間隔為1 s．

根據(jù)離散后的傳遞函數(shù)表達式(8)，可得到輸出電壓Ut(k)、開路電壓UOCV(k)和電路電流I(t)間的數(shù)學關系式：

E(k)=a1E(k-1)+a2I(k)+a3I(k-1)

(9)

其中

(10)

方程(9)中，E(k)表示電池輸出電壓Ut(k)與開路電壓UOCV(k)間電壓差，UOCV(k)的獲取則是通過實驗建立的SOC-OCV曲線，如圖2通過查表法獲得的；方程(10)是方程(9)中各系數(shù)的數(shù)學表達式，其余各參數(shù)含義與方程(1)、(2)中相同．

1.2 模型參數(shù)估計

如圖1中所示的一階Thevenin電路模型，該模型需要辨識的參數(shù)為Rin、Rp、Cp．室溫環(huán)境(25±1) ℃下，采用動態(tài)應力測試(dynamic stress test，DST)方法[28]對單體電池進行實驗室條件下動態(tài)充放電測試，使用變化電流激勵電池，用于模擬電池動態(tài)工作特性，實驗中采用的電池信息見表1，采用的DST放電策略見表2，在該放電策略下測得的電池兩端電壓、電流數(shù)據(jù)如圖3所示．

表2 DST工況測試策略

表1 實驗電池詳細參數(shù)

本文中采用小電流充放電方法對電池SOC與OCV特性進行研究，此方法一定程度上可以降低極化現(xiàn)象對電池端電壓的影響，還可以提高實驗效率[29]．根據(jù)文獻[30]中的21點數(shù)據(jù)擬合方法，選取實驗中測得的21個SOC-OCV離散數(shù)據(jù)點，采用9階多項式對SOC-OCV間關系式進行擬合：

UOCV(k)=p1SOC9(k)+…+p9SOC(k)+p10

(11)

SOC-OCV擬合曲線結(jié)果如圖2所示，擬合后方程(11)各系數(shù)如表3所示．

表3 SOC-OCV多項式擬合系數(shù)

表4為SOC-OCV曲線擬合效果優(yōu)劣的評定參數(shù)，誤差平方和越接近0，表明模型的擬合效果越好；決定系數(shù)取值范圍為[0 1]，越接近1說明模型對數(shù)據(jù)的擬合效果越好；校正后決定系數(shù)對模型評價的效果與決定系數(shù)相同，取值范圍也是[0 1]，該值越大模型的擬合效果越好；擬合標準差評價效果與誤差平方和相同，通常也是越小越好．從以上評價指標數(shù)值可以看出，方程(11)對SOC與OCV的關系有良好的擬合效果．

表4 SOC-OCV擬合曲線評定指標

為保證模型SOC估算時初值可靠，同時避免電池模型參數(shù)估計過程中舊數(shù)據(jù)累計造成下一次參數(shù)估計無法進行有效的修正情況，本文中模型參數(shù)辨識采用帶遺忘因子的最小二乘方法．此方法能夠通過引入遺忘因子減弱歷史數(shù)據(jù)累計對下次更新時參數(shù)權重的影響，利用新觀測數(shù)據(jù)對舊估計值進行修正．在進行參數(shù)估計時，將輸入數(shù)據(jù)定義為u(k)，待辨識參數(shù)定義為θ(k)，y(k)定義為輸出參數(shù)，利用帶遺忘因子的最小二乘方法與方程(9)、(10)對電路模型中參數(shù)Rin、Rp和Cp進行辨識，辨識結(jié)果如圖4所示．

動態(tài)參數(shù)辨識后需要對提出模型的精度進行驗證，將圖4中辨識后模型參數(shù)代入方程(2)中，得到模型估計下的輸出電壓值如圖5(a)所示，并將得到的估計電壓與DST工況下的電壓進行對比，得到兩者的誤差如圖5(b)所示．

從圖5(a)中可以看出僅在4C放電或2C充電時，模型的估計電壓誤差最大，但此類電流在一個循環(huán)周期中作用時間非常短，其余時刻均能實現(xiàn)對實際電壓良好估計，電壓的平均估計誤差僅為0.001 6 V，因此可對本文SOC估計方法驗證．

2 基于改進粒子濾波算法的SOC估計

卡爾曼濾波(Kalman filter，KF)是一種線性系統(tǒng)和高斯噪聲條件下運用最小均方誤差原則進行系統(tǒng)動態(tài)估計的方法，實際應用中多數(shù)系統(tǒng)為非線性系統(tǒng)，常用處理方式是利用線性化方法將非線性問題轉(zhuǎn)化為近似線性問題進行次優(yōu)求解運算．EKF是一種常用的非線性處理方法，當系統(tǒng)維數(shù)較高時，系統(tǒng)的觀測噪聲協(xié)方差陣和過程噪聲協(xié)方差陣容易出現(xiàn)非正定，會使濾波器發(fā)散；另一種近似非線性濾波方法是UKF，由于該方法非線性函數(shù)轉(zhuǎn)化過程中未忽略高階項，得到的估計結(jié)果精度要高于EKF，但是實際應用中的多數(shù)非線性系統(tǒng)，其過程噪聲和觀測噪聲并不服從高斯分布，所以KF及其擴展方法在非高斯分布的系統(tǒng)中并不適用．PF算法在求解非線性特性與非高斯條件下系統(tǒng)狀態(tài)估計問題時較為常用，它通過對系統(tǒng)概率密度函數(shù)采樣得到的集合進行預測與更新實現(xiàn)狀態(tài)估計．PF算法中的重要密度函數(shù)是先驗信息概率密度函數(shù)，如果測量精度沒有嚴格要求時，可以得到滿意的結(jié)果，但PF算法由于沒有考慮當前測量值，使得以先驗密度函數(shù)采樣獲得的樣本集合和以實際后驗概率密度采樣的粒子集合兩者之間產(chǎn)生嚴重偏差，同時重要性采樣中的權重方差也會隨時間變大，導致樣本中粒子的權重發(fā)生退化，在幾輪迭代后粒子集中部分粒子權重可能退化至很小，不僅浪費計算資源，同時采樣后的粒子集也不能真實表示粒子集合后驗概率密度的分布情況．

粒子退化是PF算法存在的一個嚴重問題，可以采用一些方法對其進行改進，常用的方法有增加采樣粒子數(shù)目、重采樣和選用合理建議密度函數(shù)．重采樣方法能夠減少粒子集合中權重較低的粒子，關注權重較高的粒子，然而這種方法會降低粒子的多樣性；而增加粒子數(shù)目抑制權重退化又會使算法計算時間大大增加．選擇合理的建議密度分布函數(shù)，通過重要性重采樣技術得到一組樣本點集“覆蓋”真實粒子狀態(tài)，便能保證采樣的有效性和結(jié)果的可靠性．本文從建議密度函數(shù)角度改進PF算法，分別以EKF、UKF作為改進濾波算法EPF和UPF的建議密度函數(shù)，與PF算法中以先驗信息密度函數(shù)作為重要密度函數(shù)不同，EKF算法能夠根據(jù)非線性問題的一階泰勒展開形式，結(jié)合最新量測值對后驗分布進行估計，并將近似后驗密度函數(shù)作為建議密度函數(shù)實現(xiàn)粒子更新；UKF相比EKF，在處理非線性問題時采用無跡變換處理一步預測中的均值和方差，結(jié)合量測值采用一系列樣本近似逼近后驗概率密度函數(shù)，避開了EKF方法中非線性系統(tǒng)一階泰勒展開而忽略高階項引起的誤差問題，并以得到的后驗概率密度函數(shù)作為建議密度函數(shù)對粒子進行更新，對電池SOC進行估計．

本研究中采用的EPF和UPF算法執(zhí)行過程如下．

2.1 EPF算法執(zhí)行步驟

EPF算法執(zhí)行步驟如下：

(1)初始化，從先驗分布中抽取初始化狀態(tài)：

(12)

(2)用EKF更新粒子：

(13)

(3)更新產(chǎn)生粒子：

(14)

(4)計算粒子權值：

(15)

(5)歸一化粒子權值：

(16)

(6)狀態(tài)估計：

(17)

(7)重采樣：

(18)

2.2 UPF算法執(zhí)行步驟

UPF算法執(zhí)行步驟如下：

(2)k=1，2，…

①重要性采樣：i=1，2，…，N，使用UKF算法更新粒子．

時間更新：

(19)

測量更新：

(20)

③輸出融合結(jié)果：

(21)

系統(tǒng)離散處理后，得到電池模型的狀態(tài)方程和觀測方程，如式(22)和(23)：

狀態(tài)方程

xk+1=f(xk，uk)+wk

(22)

觀測方程

yk+1=h(xk+1，uk+1)+vk+1

(23)

式(22)、(23)中，xk表示系統(tǒng)狀態(tài)向量，uk表示輸入向量，yk表示觀測向量．f(·)是狀態(tài)方程映射函數(shù)，h(·)是觀測方程映射函數(shù)，wk和vk分別模擬系統(tǒng)過程噪聲與系統(tǒng)觀測噪聲．

當SOC(t0)已知時，SOC(t)可由方程(3)計算得到，根據(jù)方程(1)、(2)，系統(tǒng)狀態(tài)方程和相應的觀測方程可寫為式(24)～(27)：

(24)

Yk+1=(1 -1 -I(k+1))·

(25)

(26)

(27)

式中：Δt是數(shù)據(jù)采樣間隔；Qn為電池的額定容量；wSOCk、wUp,k、wRin，k為過程噪聲，vk+1為觀測噪聲，過程噪聲協(xié)方差為(10-410-510-5)，觀測噪聲協(xié)方差為10-4；其余模型參數(shù)定義與前述相同．

3 實驗結(jié)果分析

為驗證研究中所用兩種改進后粒子濾波算法的有效性，本文采用第1章中的DST工況數(shù)據(jù)對電池的SOC進行估計，實驗中電池電流、電壓數(shù)據(jù)如圖3所示．根據(jù)電池模型狀態(tài)方程和觀測方程，將辨識后的參數(shù)代入模型中，按照EPF、UPF算法流程執(zhí)行，可以得到DST工況下，基于EKF、UKF、PF、EPF和UPF共5種算法的電池SOC估計值，如圖6所示．

從圖6曲線可以看出，SOC估計所用的5種算法中，PF算法及其改進的EPF和UPF算法下SOC估計值波動程度要小于EKF和UKF算法，對參考SOC的追蹤能力好于KF類算法．圖7為5種濾波算法下的SOC與參考SOC兩者間偏差D對比圖，由圖7可知，EKF、UKF和PF濾波算法中，PF算法估計偏差曲線波動最小，EKF算法估計偏差曲線波動最大，而UKF算法估計效果介于EKF和PF算法之間．PF、EPF和UPF算法SOC估計偏差對比結(jié)果表明，改進后EPF與UPF算法估計的結(jié)果要好于PF算法下的估計結(jié)果．以UKF為建議密度函數(shù)的UPF算法SOC估計偏差明顯小于以EKF為建議密度函數(shù)的EPF算法的估計偏差，并且UPF方法的偏差波動范圍較小，穩(wěn)定性好于EPF．圖8是5種濾波算法估計誤差e和估計誤差統(tǒng)計特征量f分析，從圖8(a)中可以看出SOC的估計誤差與圖7中偏差的分析結(jié)論相似，EPF和UPF兩種改進后的PF算法在進行估計時穩(wěn)定性高于EKF、UKF和PF 3種傳統(tǒng)濾波算法，從兩種改進后的算法估計誤差曲線看出，UPF算法估計誤差和誤差曲線波動程度要小于EPF算法．圖8(b)是5種算法SOC估計誤差結(jié)果的統(tǒng)計特征量，分別是SOC估計誤差的均值、標準差和最大估計誤差，采用EKF、UKF和PF算法估計情況下，平均誤差分別為3.26%、2.19%和1.36%，標準差為2.14、1.68和1.35，表明PF算法的估計效果優(yōu)于EKF和UKF；而PF、EPF和UPF 3種算法下估計誤差結(jié)果的均值、標準差分別為1.36%、1.09%、0.71%和1.35、1.25、0.79，EPF和UPF的估計精度相對于PF算法分別提高19.8%、47.8%，說明以UKF為建議密度函數(shù)進行重采樣的UPF算法SOC估計效果優(yōu)于以EKF為建議密度函數(shù)進行重采樣的EPF估計效果，兩者與PF算法相比在一定程度上抑制了粒子退化現(xiàn)象．

4 結(jié) 語

本文選擇Thevenin一階等效電路模型作為電池研究模型，提出了基于EPF、UPF算法的動力電池SOC估計方法．在實驗室DST數(shù)據(jù)條件下，采用帶遺忘因子的最小二乘方法辨識電池模型中未知參數(shù)，將辨識后的參數(shù)代入模型狀態(tài)方程和觀測方程，結(jié)合改進后的EPF和UPF算法，對電池的SOC進行估計．結(jié)果表明，以UKF為建議密度函數(shù)進行重采樣后的UPF算法得到的SOC估計誤差均值為0.71%，標準差為0.79，低于以EKF為建議密度函數(shù)重采樣的EPF算法下的均值1.09%和標準差1.25；與PF算法相比，UPF算法通過選擇合理的建議密度函數(shù)抑制了PF算法中粒子權重退化的現(xiàn)象，提高了SOC估計的精度，比以EKF為建議密度函數(shù)的EPF方法在SOC估計中更有意義．