王順利，謝 非，陳 蕾，張 麗，王 露

（1.西南科技大學(xué) 信息工程學(xué)院/檢測研究所，綿陽 621010；2.西南科技大學(xué) 信息工程學(xué)院，綿陽 621010）

0 引言

當(dāng)航空鋰電池處于放電情況下，電池的SOC是電池的一個重要參數(shù)，能夠清楚的表示電池的當(dāng)前剩余電量，可以為電池管理系統(tǒng)提供控制策略依據(jù)。那么能否準(zhǔn)確的估計電池的SOC就顯得至關(guān)重要。

國外在SOC估算方面的研究相對較成熟，美國的GregoyL.Plett考慮到電池參數(shù)存在的非線性特性，采用擴(kuò)展卡爾曼濾波算法，用于SOC估算。他還提出了“sigma”點卡爾曼濾波算法，能夠更加準(zhǔn)確地預(yù)測鋰電池組的SOC值。德國研究人員認(rèn)為電池管理系統(tǒng)的功能應(yīng)該包括均衡充電、估算電池荷電狀態(tài)以及健康狀況，并根據(jù)SOC和SOH以及電池溫度來控制電池的放電電流等。

在國內(nèi)，也有大量高校、企業(yè)和科研院所對SOC估算方法進(jìn)行了相關(guān)的設(shè)計研究。

鄧濤等人結(jié)合傳統(tǒng)的安時積分法和開路電壓法，充分考慮了電池充放電效率和電池使用溫度等因素，提出了一種帶有補償措施的基于安時積分法的SOC估算方法，使得SOC估算值的偏差在3%以內(nèi)，該估算方法精確度較高，為電池SOC的在線實時估算和檢測提供了重要參考。

本文研究主從無跡卡爾曼濾波（UKF）對航空鋰電池的SOC進(jìn)行估算，在simulink上建立航空鋰電池的等效電路模型，進(jìn)行仿真實驗。仿真實驗表明該算法SOC估算精度誤差穩(wěn)定在百分之五左右。

1 理論與方法研究

1.1 卡爾曼濾波技術(shù)研究

卡爾曼濾波是一種遞推線性最小方差估計的運算方法，該法利用系統(tǒng)觀測量從中提獲取所需要估算的系統(tǒng)狀態(tài)向量。該方法引入系統(tǒng)狀態(tài)向量和系統(tǒng)觀測向量理論，可以建立狀態(tài)方程和觀測方程，是一種最優(yōu)化自回歸數(shù)據(jù)處理算法?？柭鼮V波法用反饋控制的方法估計過程狀態(tài)。因此卡爾曼濾波法可以分為兩個部分：狀態(tài)方程和測量更新方程。

在式(1)中，k為離散時間，X（k+1）和X（k）分別為系統(tǒng)在k+1和k時刻的狀態(tài)，W(k)為輸入白噪聲。

在式(2)中，k為離散時間，Y（k）為對應(yīng)時刻的觀測狀態(tài)，X（k）為系統(tǒng)k時刻的狀態(tài)，V（k）為觀測噪聲。

由于航空鋰電池具有很強的非線性，普通的線性卡爾曼濾波算法并不能很好的保證精度。無跡卡爾曼濾波算法對非線性的系統(tǒng)具有良好的濾波效果。因此我們選擇了無跡卡爾曼濾波算法進(jìn)行處理。

1.2 對Thevenin模型進(jìn)行研究

根據(jù)對Thevenin模型的分析得到下面的式子：

同時又因為E(t)為電源電動勢，它與SOC存在著非線性函數(shù)的關(guān)系式為E(t)=f(SOC(t))，再結(jié)合安時積分法

選擇SOC為系統(tǒng)狀態(tài)量，端電壓V（t）作為觀測量，i(t)作為系統(tǒng)輸入量，但由于我們這是個離散的系統(tǒng)，所以經(jīng)過離散化處理后我們得到以下關(guān)系：

式(6)是系統(tǒng)狀態(tài)方程，S（k+1）和S（k）分別代表k+1以及k時刻系統(tǒng)的SOC值，Q是鋰電池的容量，ΔT是單位時間內(nèi)（采樣時間），i（k）系統(tǒng)電流大小，V（k）為噪聲。

式(7)是系統(tǒng)測量方程，V（k）為k時刻電源電動勢的值，f（SOC(k)）為電壓與SOC之間的非線性關(guān)系，R1為歐姆內(nèi)阻，i（k）為k時刻測得的電流值，Uc（k）為k時刻開路電壓值，V（k）為噪聲。

1.3 無跡卡爾曼算法研究

對于一個非線性系統(tǒng)來說，它的狀態(tài)方程和測量方程分別為：

k為當(dāng)前所處時刻，f（XK-1，UK）為非線性系統(tǒng)狀態(tài)轉(zhuǎn)移方程，g（Xk-1）為非線性測量方程，Xk為狀態(tài)變量，Uk為已知的輸入，Yk為測量信號；Wk叫過程噪聲，Vk測量噪聲。我們假定Wk和Vk是不相關(guān)的均值均為零高斯白噪音，其協(xié)方差分別為Qw，Rv。具體過程如下：

1）初始值計算：

2）建立Sigma點：

3）狀態(tài)更新方程如下：

式子中的k/k-1為基于k-1時刻對k時刻的估計值。

4）測量更新方程如下：

5）重復(fù)上面四個步驟。

由上述公式可知，只要給定初始條件X0和P0，根據(jù)k-1時刻的狀態(tài)值、k時刻獲得的輸入值和觀測值，就可以估算k時刻的狀態(tài)最優(yōu)估計值Xk。

2 系統(tǒng)設(shè)計與實驗分析

2.1 電源等效電路模型的選擇

鋰電池等效電路模型的選擇對航空鋰電池的SOC估算系統(tǒng)來說十分重要。工程上比較常用而且易于實現(xiàn)的電池等效電路模型有Rint模型，RC模型，Thevenin模型和PNGV模型。Rint模型是一種理想狀態(tài)的簡單模型，U1來表示電池的開路電壓，R1來表示電池的內(nèi)阻，但在實際中電池的電壓值、電流值都是隨時變動的，不是靜止不變的，故該模型誤差較大，不適合做高精度模型應(yīng)用。

圖1 Rint模型

RC模型認(rèn)為電池可由電容Cc和Cb，電阻R1、R2、Rc等效組成。Cc為一個取值非常的電容描述電池儲存的電量，Rc和Cc串聯(lián)作為儲能單元的計劃電阻；Cb表示電池極化電容；R1、R2表示電池內(nèi)阻。

圖2 RC模型

Thevenin模型，該模型參數(shù)描述如下：Uc表示電池的開路電壓，R1表示電池歐姆電阻，R2表示電池極化電阻，C表示電池極化電容。

圖3 Thevenin模型

PNGV模型，該模型在Thevenin模型基礎(chǔ)上添加一個電容Cb，Cb表示負(fù)載電流的變化而產(chǎn)生的開路電壓變化也可用1/OVC表示，Ub表示電容Cb兩端的電壓?？芍狿NGV模型參數(shù)數(shù)學(xué)關(guān)系式：

在此，對選擇電池等效電路模型需要考慮兼顧以下三點要求：

1）模型最終是為了精確估算電池SOC值，因此需要盡可能準(zhǔn)確的用各參數(shù)反映電池內(nèi)部電壓、電流、內(nèi)阻等特點；

2）電池需要適應(yīng)動態(tài)充放電環(huán)境中，需要電池模型能準(zhǔn)確的計算出電池在動態(tài)、變電流充放電電流的情況下，精確計算電池剩余電量；

3）電池模型結(jié)構(gòu)盡量精簡易行，在不影響估算電池SOC精度前提下盡量精簡模型階數(shù)。

我們綜合考慮了本實驗的精度要求及模型的復(fù)雜程度，最終選擇了Thevenin模型。Thevenin模型的關(guān)系式如下：

2.2 實驗分析

圖4 soc-voc關(guān)系圖

為了驗證基于Thevenin模型的UKF算法的有效性。根據(jù)對航空鋰電池的實驗數(shù)據(jù)計算得到了電源電動勢與SOC的關(guān)系。

圖5為卡爾曼濾波系統(tǒng)框圖，描述了卡爾曼濾波算法的過程。

圖6為航空鋰電池電池部分的模型圖。

由圖7分析可知，卡爾曼濾波（紅色）得到的soc值非常接近于真實的soc值。達(dá)到誤差在5%之內(nèi)的目標(biāo)。

3 結(jié)束語

圖5 卡爾曼濾波系統(tǒng)框圖

圖6 電池模型

圖7 卡爾曼濾波波形圖

我們討論了關(guān)于航空鋰電池的等效電路模型及關(guān)于鋰電池的SOC的常用幾種估算方法，較為詳細(xì)的介紹了線性卡爾曼濾波及無跡卡爾曼濾波。動力電池等效模型對SOC估算具有重要意義。在建模方面，此處采用的Thevenin模型考慮了極化等因素對估算的影響，在算法上，在UKF中將比例修正法加入對稱采樣中 ，避免了局部效應(yīng)問題；無跡卡爾曼在蓄電池的SOC估算系統(tǒng)中比擴(kuò)展卡爾曼更容易實現(xiàn)，無跡卡爾曼并能實現(xiàn)更高的狀態(tài)估計精度?？梢灶A(yù)見，基于適合的電池等效模型，無跡卡爾曼在其他各類蓄電池的估計方面也有著廣闊的應(yīng)用空間，因此進(jìn)一步實現(xiàn)基于的估算方法的工程化是很有必要的。

[1]黃小平,王巖.卡爾曼濾波原理及應(yīng)用[M].電子工業(yè)出版社,2015.

[2]何志超,楊耕,盧蘭光.基于恒流外特性和SOC的電池直流內(nèi)阻測試方法[J].清華大學(xué)學(xué)報,2015;55(5):532-537.

[3]Nejad S, Gladwin DT, Stone DAA systematic review of lumpedparameter equivalent circuit models for real-time estimation of lithium-ion battery states[J].Journal of Power Sources,2016,316:183-196.

[4]Barai A, Chouchelamane GH, Guo Y.A study on the impact of lithium-ion cell relaxation on electrochemical impedance spectroscopy[J].Journal of Power Sources,2015,280:74-80.

[5]Sun BX, Jiang JC, Zheng FD. Practical state of health estimation of power batteries based on Delphi method and grey relational grade analysis[J].Journal of Power Sources,2015,282:146-157.

[6]Yuan SF, Wu HJ, Ma XR. Stability analysis for li-ion battery model parameters and state of charge estimation by measurement uncertainty consideration[J].Energies,2015,8(8):7729-7751.

[7]Bazinski SJ, Wang X. Experimental study on the influence of temperature and state-of-charge on the thermophysical properties of an LFP pouch cell[J].Journal of Power Sources,2015,293:283-291.

[8]Marongiu A, Roscher M, Sauer DU. Influence of the vehicleto-grid strategy on the aging behavior of lithium battery electric vehicles[J].Applied Energy,2015,137:899-912.