趙又群　周曉鳳　劉英杰

南京航空航天大學(xué),南京,210016

基于擴(kuò)展卡爾曼粒子濾波算法的鋰電池SOC估計(jì)

趙又群周曉鳳劉英杰

南京航空航天大學(xué),南京,210016

鋰電池荷電狀態(tài)用來描述電池剩余電量的多少，進(jìn)而反映電動(dòng)汽車的續(xù)駛里程，是電池管理系統(tǒng)中的核心參數(shù)。電池循環(huán)次數(shù)、瞬間大電流以及溫度等因素都會(huì)使電池特性發(fā)生變化,使用擴(kuò)展卡爾曼濾波算法對(duì)電池荷電狀態(tài)進(jìn)行估計(jì)，會(huì)有較大的誤差甚至導(dǎo)致算法不收斂。為了有效地抑制發(fā)散以及噪聲的影響,基于鋰電池混合噪聲模型，應(yīng)用擴(kuò)展卡爾曼粒子濾波算法對(duì)鋰電池荷電狀態(tài)和電流漂移噪聲進(jìn)行同步估計(jì)。最后根據(jù)充放電試驗(yàn)數(shù)據(jù)進(jìn)行仿真分析，結(jié)果證明了該算法的優(yōu)越性。

鋰電池；荷電狀態(tài)；混合噪聲模型；擴(kuò)展卡爾曼粒子濾波

0　引言

電動(dòng)汽車以環(huán)保、節(jié)能及輕便等特點(diǎn)被人們所青睞，而作為關(guān)鍵技術(shù)部件的電池，其性能直接影響到整車性能的好壞。鋰離子電池具有電壓高、比能量高、充放電壽命長(zhǎng)、無記憶效應(yīng)、無污染、自放電率低、工作溫度范圍寬和安全可靠等優(yōu)點(diǎn)，是目前電動(dòng)汽車首選的核心動(dòng)力源[1]。它具有高度的非線性，且受溫度、充放電倍率以及容量等因素影響，其動(dòng)態(tài)過程無法用單一的模型來描述。因此，如何利用電池的可測(cè)參數(shù)來實(shí)現(xiàn)當(dāng)前電池剩余電量準(zhǔn)確估計(jì)，一直以來是電動(dòng)汽車電池管理系統(tǒng)的核心問題和急需解決的技術(shù)難點(diǎn)，也是一項(xiàng)重要而富有挑戰(zhàn)性的任務(wù)[2]。只有準(zhǔn)確估計(jì)荷電狀態(tài)(SOC)才能起到優(yōu)化電池性能、提高電池安全性和延長(zhǎng)電池使用壽命等作用。

目前,工程上最常用的算法是A·h計(jì)量法[3-5],它是一種開環(huán)估計(jì)方法，簡(jiǎn)單易實(shí)現(xiàn),但累計(jì)誤差大,估計(jì)精度低。針對(duì)誤差引起的算法發(fā)散問題,文獻(xiàn)[6]提出基于電池狀態(tài)空間模型的卡爾曼濾波算法,實(shí)現(xiàn)了電池SOC的閉環(huán)估計(jì)。文獻(xiàn)[7-11]均采用擴(kuò)展卡爾曼濾波算法對(duì)電池SOC進(jìn)行估計(jì)，使得極化效應(yīng)的影響大大減弱，提高了電池SOC計(jì)算的精度。為減小鋰電池非線性引起的誤差，也有研究人員將粒子濾波算法應(yīng)用到電池SOC估計(jì)中。文獻(xiàn)[12]為了較好地解決粒子濾波算法中粒子的退化問題，根據(jù)遺傳算法的原理，提出了遺傳粒子濾波算法并給出其估計(jì)電池SOC的實(shí)現(xiàn)方法，通過實(shí)例驗(yàn)證，并與標(biāo)準(zhǔn)粒子濾波算法進(jìn)行對(duì)比分析，結(jié)果顯示遺傳粒子濾波算法對(duì)電池SOC的估計(jì)更加準(zhǔn)確，精度更高。

電池實(shí)際工作過程中，溫度、循環(huán)次數(shù)、瞬間大電流等因素使電池的極化內(nèi)阻、可使用容量等存在較大差異，在估計(jì)電池SOC時(shí)往往存在較大誤差甚至發(fā)散。為抑制系統(tǒng)非線性和發(fā)散，單純地使用粒子濾波算法在粒子數(shù)目比較少時(shí)可能會(huì)出現(xiàn)粒子匱乏的現(xiàn)象，而且算法對(duì)模型依賴性強(qiáng)并且存在抖動(dòng)。因此為了消除算法抖動(dòng)和粒子退化的現(xiàn)象，本文基于鋰電池混合噪聲模型，應(yīng)用擴(kuò)展卡爾曼粒子濾波(EKPF)算法對(duì)鋰電池SOC進(jìn)行估計(jì)，并將其與擴(kuò)展卡爾曼濾波(EKF)算法和粒子濾波(PF)算法進(jìn)行比較分析。

1　鋰電池混合噪聲模型

電池SOC是不能直接測(cè)量的，必須通過溫度、電流、電壓等能測(cè)量的參數(shù)進(jìn)行估算，因此，建立一個(gè)準(zhǔn)確的電池模型對(duì)精確估算電池SOC十分重要。一個(gè)理想的電池模型應(yīng)該能夠較好地反映電池的各種特性,而且階數(shù)不能太高,便于計(jì)算,易于工程實(shí)現(xiàn)。本文借鑒文獻(xiàn)[2]提出的結(jié)合Shepherd模型、Unnewehr模型和Nernst模型的復(fù)合模型,考慮電流漂移噪聲的干擾,建立混合噪聲模型,具體數(shù)學(xué)表達(dá)式如下:

K3lnxk+K4ln(1-xk)

(1)

式中，yk為電池工作電壓;ireal為真實(shí)放電電流；R為電池內(nèi)阻,本文考慮充放電內(nèi)阻的不同,充電時(shí)R=Rc,ireal為負(fù),放電時(shí)R=Rd,ireal為正;K0、K1、K2、K3、K4為模型參數(shù);xk為k時(shí)刻的瞬時(shí)荷電狀態(tài)。

為了表征溫度、充放電倍率、容量等對(duì)電池SOC估計(jì)的影響，對(duì)SOC作如下定義：

(2)

κ=|i(t)|/CN

式中,δSOC為電池的SOC值;η為充放電效率;T為溫度;κ為充放電倍率;CN為常溫下電池的標(biāo)稱容量；C(T，κ，t)為在不同溫度、不同充放電倍率下的電池可用容量。

為了抑制電流漂移噪聲的干擾，將電池荷電狀態(tài)及電流漂移噪聲一起作為狀態(tài)變量進(jìn)行同步估計(jì)[13]。狀態(tài)方程中的電流應(yīng)當(dāng)是除去電流漂移的電流真實(shí)值,即

(3)

放電電流的測(cè)量不可避免存在漂移和噪聲。將放電電流的測(cè)量值i分成兩個(gè)部分，即電流真實(shí)值ireal和電流噪聲值is,則

i=ireal+is

(4)

將式(4)代入式(3)得

(5)

將式(5)作離散化處理得

(6)

式中,ik-1為k-1時(shí)刻的測(cè)量電流,is,k-1為k-1時(shí)刻的電流噪聲值;Δt為時(shí)間間隔。

該混合噪聲模型能有效抑制電流漂移的干擾,簡(jiǎn)單易操作,計(jì)算量不大，易于工程實(shí)現(xiàn)，而且模型中的參數(shù)易于辨識(shí),只要獲得N組電流、電壓和真實(shí)SOC數(shù)據(jù),就可以由最小二乘法辨識(shí)得到模型參數(shù)。

2　模型參數(shù)辨識(shí)

系統(tǒng)辨識(shí)是在已知輸入和輸出數(shù)據(jù)的基礎(chǔ)上，在某種準(zhǔn)則條件下，估計(jì)出系統(tǒng)模型的未知參數(shù)的。目前常用的辨識(shí)方法有最小二乘法、梯度校正法、極大似然法等。最小二乘法由于原理簡(jiǎn)單、收斂較快、易于理解、易于編程實(shí)現(xiàn)等特點(diǎn)，在系統(tǒng)參數(shù)估計(jì)中應(yīng)用相當(dāng)廣泛[14]。

對(duì)于電動(dòng)汽車來說，被控對(duì)象通常都可以不斷提供新的輸入輸出數(shù)據(jù)，而且還希望利用這些新的信息來改善估計(jì)精度,因此需要在線辨識(shí)電池模型的參數(shù)。本文采用遞推最小二乘法對(duì)鋰電池復(fù)合模型進(jìn)行參數(shù)辨識(shí)，遞推公式為

(7)

(8)

P(k)=(I-K(k)φT(k))P(k-1)

(9)

遞推最小二乘法進(jìn)行模型參數(shù)辨識(shí)一般步驟[14]如下：

(2)采樣當(dāng)前輸出y(k)和輸入u(k);

(4)k→k+1,返回步驟(2),繼續(xù)循環(huán)。

本文對(duì)額定電壓3.2 V、額定容量20 A·h的鋰離子電池的充放電試驗(yàn)數(shù)據(jù),采用遞推最小二乘法對(duì)建立的混合噪聲模型進(jìn)行參數(shù)辨識(shí)，辨識(shí)結(jié)果如表1所示。

表1　鋰離子電池復(fù)合模型參數(shù)表

3　基于EKPF算法的SOC估計(jì)

EKPF算法是解決非線性系統(tǒng)狀態(tài)變量估計(jì)問題最常用的方法,粒子濾波擺脫了解決非線性濾波問題時(shí)隨機(jī)變量必須滿足高斯分布的制約條件,能夠應(yīng)用在任意非線性非高斯隨機(jī)系統(tǒng)中,在許多領(lǐng)域取得了成功應(yīng)用。粒子濾波器中密度分布的建議值一般決定了粒子濾波器應(yīng)用的性能。為了避免粒子匱乏和算法抖動(dòng)，本文以擴(kuò)展卡爾曼濾波器作為建議分布，實(shí)現(xiàn)基于擴(kuò)展卡爾曼粒子濾波算法對(duì)鋰電池SOC的估計(jì)。

首先建立鋰電池的離散狀態(tài)空間模型。鋰電池SOC作為狀態(tài)變量,負(fù)載電壓作為觀測(cè)變量,將上文建立的復(fù)合模型的狀態(tài)空間方程進(jìn)行離散化。則狀態(tài)方程：

(10)

觀測(cè)方程：

K3ln(xk)+K4ln(1-xk)+vk

(11)

其中,wk為系統(tǒng)噪聲,vk為觀測(cè)噪聲,它們均為零均值的高斯白噪聲,wk～N(0,Q),vk～N(0,R),Δt為離散系統(tǒng)采樣周期,模型參數(shù)K0、R、K1、K2、K3、K4已由上文第2節(jié)的方法辨識(shí)得到,其值詳見表1。

根據(jù)鋰離子電池狀態(tài)空間方程,利用EKPF實(shí)現(xiàn)電池SOC估計(jì)的算法流程如圖1所示。

圖1　EKPF算法流程圖

具體算法步驟如下:

(12)

式中,g(·)為已知電流和荷電狀態(tài)情況下，由混合噪聲模型計(jì)算得到的電壓值。

然后,進(jìn)行歸一化處理：

(13)

(14)

若Ne f f

(15)

(4)循環(huán)迭代k→k+1,轉(zhuǎn)步驟(2)循環(huán)迭代。

4　仿真試驗(yàn)及結(jié)果分析

為了驗(yàn)證EKPF算法對(duì)鋰電池SOC估計(jì)的精確性,對(duì)3.2 V/20 A·h的鋰電池進(jìn)行放電試驗(yàn)。首先將電池組充滿電(δSOC(0)=1),然后在恒溫條件下對(duì)電池組放電,放電電流波形如圖2所示,整個(gè)放電時(shí)間為105s,采樣間隔Δt為10 s。

圖2　放電電流波形

試驗(yàn)中,粒子數(shù)N設(shè)定為50,系統(tǒng)噪聲和觀測(cè)噪聲的方差分別取2×10-14、25×10-6，分別采用EKF、PF和EKPF三種算法對(duì)電池SOC進(jìn)行估計(jì),估計(jì)結(jié)果圖如圖3所示。

圖3　三種不同算法的SOC估計(jì)結(jié)果

圖3中的真實(shí)SOC值是通過放電試驗(yàn)法得到的。對(duì)電池進(jìn)行持續(xù)放電，測(cè)出在一定的很小的時(shí)間間隔后釋放的電量，再根據(jù)SOC的定義便可計(jì)算得到該時(shí)刻的真實(shí)SOC值。為了定量地比較EKF、PF和EKPF的性能，定義試驗(yàn)的均方根誤差為

(16)

為了考察算法的時(shí)間復(fù)雜度,定義第k次試驗(yàn)算法的耗時(shí)時(shí)間為tk,則算法的單次平均消耗時(shí)間為

(17)

通過計(jì)算,EKF、PF和EKPF三種不同電池SOC估計(jì)算法的均方根誤差和單次平均消耗時(shí)間的值如表2所示。

表2　EKF、PF和EKPF估計(jì)SOC性能比較

從圖3和表2的仿真試驗(yàn)結(jié)果可以看出,基于鋰電池混合噪聲模型的EKPF算法精度較高,既說明了算法的優(yōu)越性,也表明了模型的準(zhǔn)確性。單獨(dú)使用PF算法時(shí)抖動(dòng)很明顯，而EKPF算法則比較平滑，在均方根誤差的指標(biāo)上也優(yōu)于EKF算法和PF算法，但由于算法復(fù)雜，計(jì)算量相對(duì)較大，故單次平均消耗時(shí)間相對(duì)較長(zhǎng)。

5　結(jié)束語

外界溫度、循環(huán)次數(shù)、瞬間大電流等因素的影響，會(huì)使電池的極化內(nèi)阻、可使用容量存在較大差異，在估計(jì)電池SOC時(shí)往往存在較大誤差甚至發(fā)散。為抑制電流漂移的干擾、系統(tǒng)非線性、算法發(fā)散以及算法本身的抖動(dòng)，本文基于鋰電池混合噪聲模型，應(yīng)用EKPF算法對(duì)鋰電池SOC進(jìn)行估計(jì)。仿真試驗(yàn)結(jié)果表明,EKPF算法能有效地估計(jì)鋰電池荷電狀態(tài)。

[1]郭炳焜,徐徽,王先友,等.鋰離子電池[M].長(zhǎng)沙:中南大學(xué)出版社,2002.

[2]Cai C H,Du D,Liu Z Y.Battery State-of-charge(SOC) Estimation Using Adaptive Neuro-fuzzy Inference System(ANFIS)[C]//IEEE International Conference on Fuzzy Systems. St. Louis,MO,United States,2003:1068-1073.

[3]林成濤,王軍平,陳全世,等.用改進(jìn)的安時(shí)計(jì)量法估計(jì)電動(dòng)汽車動(dòng)力電池SOC[J].清華大學(xué)學(xué)報(bào),2006,46(2):247-251.

Lin Chengtao,Wang Junping,Chen Quanshi,et al.Improved Ah Counting Method for State of Charge Estimation of Electric Vehicle Batteries[J].Journal of Tsinghua University(Science and Technology),2006,46(2):247-251.

[4]Mark E.Adaptive State of Charge Algorithm for Nickel Metal Hydride Batteries Including Hysteresis Phenomena[J].Journal of Power Source,2004,126(1):236-249.

[5]曾潔,卜凡濤.基于多項(xiàng)式回歸算法的鋰電池SOC估測(cè)[J].大連交通大學(xué)學(xué)報(bào),2011,32(4):70-74.

Zeng Jie,Bu FanTao.Study of SOC Estimation and Measurement of Li-Ion Battery Based on Polynomial Regression Algorithm[J].Journal of Dalian Jiaotong University,2011,32(4):70- 74.

[6]Plett G L.Extended Kalman Filtering for Battery Management Systems of LiPB-based HEV Battery Packs[J].Journal of Power Sources,2004,134(2):277-292.

[7]范波,田曉輝,馬建偉.基于EKF的動(dòng)力鋰電池SOC狀態(tài)預(yù)測(cè)[J].電源技術(shù),2010,34(8):797-799.Fan Bo,Tian Xiaohui,Ma Jianwei.Ekf-Based Estimation of Lithium-Ion Traction-Battery SOC[J].Chinese Journal of Power Sources,2010,34(8):797-799.[8]夏超英,張術(shù),孫宏濤.基于推廣卡爾曼濾波算法的SOC估算策略[J].電源技術(shù),2007,31(5):414-417.

Xia Chaoying,Zhang Shu,Sun Hongtao.A Strategy of Estimating State of Charge Based on Extended Kalman Filter[J].Chinese Journal of Power Sources,2007,31(5):414-417.

[9]李超,劉忠慶.基于參數(shù)與狀態(tài)聯(lián)合估算的SOC預(yù)測(cè)法[J].電源技術(shù),2012,36(3):362-365.Li Chao,Liu Zhongqing. SOC Estimate Based on Joint Determination of Parameters and States[J].Chinese Journal of Power Sources,2012,36(3):362-365.[10]戴海峰,魏學(xué)哲,孫澤昌.基于擴(kuò)展卡爾曼濾波算法的燃料電池車用鋰離子動(dòng)力電池荷電狀態(tài)估計(jì)[J].機(jī)械工程學(xué)報(bào),2007,43(2):92-95.

Dai Haifeng,Wei Xuezhe,Sun Zechang,Estimate State of Charge of Power Lithium-Ion Batteries Used on Fuel Cell Hybrid Vehicle with Method Based on Extended Kalman Filtering[J].Chinese Journal of Mechanical Engineering,2007,143(2):92-95.

[11]Liye Wang,Lifang Wang,Chenglin Liao.Research on Improved EKF Algorithm Applied on Estimate EV Battery SOC[C]//2010 Adia-Pacific Power and Engergy Engineering Conference.Chengdu,China,2010:1-4.

[12]張潔.基于粒子濾波算法的電動(dòng)汽車剩余電量動(dòng)態(tài)估計(jì)研究[D].北京:北京交通大學(xué),2012.

[13]劉新天.電源管理系統(tǒng)設(shè)計(jì)及參數(shù)估計(jì)策略研究[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué),2011.

[14]龐中華,崔紅.系統(tǒng)辨識(shí)與自適應(yīng)控制MATLAB仿真[M].北京:北京航空航天大學(xué)出版社,2009.

(編輯郭偉)

SOC Estimation for Li-Ion Battery Based on Extended Kalman Particle Filter

Zhao YouqunZhou XiaofengLiu Yingjie

College of Energy and Power Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing,210016

As the key parameter for power battery management,the SOC of Li-ion battery described the residual capacity,and indicated the remainder driving range of electric vehicles.The cycles,instantaneous high current,abnormal temperatures and other factors would change cell characteristics,which might introduce larger errors even divergence over time if the extended Kalman filter algorithm were applied to the SOC estimation.To suppress the divergence and noise,this paper proposed a method based on EKPF algorithm to realize accurate SOC and the current drift estimation on the Li-ion battery mixed noise model.Finally,the superiority of this method was validated by simulation results.

Li-ion battery;state-of-charge(SOC);mixed noise model;extended Kalman particle filter(EKPF)

2013-03-28

2014-07-28

國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)資助項(xiàng)目(2011AA11A210,2011AA11A220)

U463< class="emphasis_italic">DOI

：10.3969/j.issn.1004-132X.2015.03.019

趙又群,男,1968年生。南京航空航天大學(xué)能源動(dòng)力學(xué)院教授、博士研究生導(dǎo)師。研究方向?yàn)檐囕v系統(tǒng)動(dòng)力學(xué)。發(fā)表論文120余篇。周曉鳳,女,1988年生。南京航空航天大學(xué)能源動(dòng)力學(xué)院車輛工程系碩士研究生。劉英杰,男,1982年生。南京航空航天大學(xué)能源動(dòng)力學(xué)院車輛工程系博士研究生。