吳華偉，張遠(yuǎn)進(jìn)，葉從進(jìn)

（1湖北文理學(xué)院純電動汽車動力系統(tǒng)設(shè)計與測試湖北省重點(diǎn)實驗室，湖北襄陽441053；2湖北文理學(xué)院汽車與交通工程學(xué)院，湖北襄陽441053）

電動汽車具有節(jié)能環(huán)保、噪音低、易保養(yǎng)、維護(hù)方便等優(yōu)點(diǎn)。鋰電池作為電動汽車的主要能源部件，對車輛的加速爬坡、剎車制動起到關(guān)鍵性作用。電池荷電狀態(tài)（state of charge，SOC）是反映鋰電池能量的重要參數(shù)，準(zhǔn)確地估算SOC 能延長電池的使用壽命，避免電池的過充/放電，是保障電動汽車安全行駛的前提。SOC 屬于電池內(nèi)部特性參數(shù)無法直接測量，而且SOC 同電池電壓、電流、溫度等參數(shù)間表現(xiàn)出強(qiáng)烈的非線性關(guān)系。因此，如何提高SOC 估算的精度和準(zhǔn)確性是當(dāng)前電動汽車領(lǐng)域亟待解決的問題。

當(dāng)前，電動汽車常用的SOC 估算方法有：安時積分法[1-2]、神經(jīng)網(wǎng)絡(luò)法[3-5]、卡爾曼濾波法[6-8]等。文獻(xiàn)[9]為了使荷電狀態(tài)量SOC 估計精確，提出以遺傳算法優(yōu)化最小二乘支持向量機(jī)的方法對電池的SOC進(jìn)行估算的模型。為解決BP 網(wǎng)絡(luò)收斂速度慢和容易陷入局部最小值的問題，文獻(xiàn)[10]提出了一種蓄電池荷電狀態(tài)估算的遺傳算法和BP 網(wǎng)絡(luò)相結(jié)合方法。文獻(xiàn)[11]以無跡卡爾曼濾波（UKF）算法為基礎(chǔ)，采用自適應(yīng)無跡卡爾曼濾波（AUKF）算法估算電動汽車電池SOC。文獻(xiàn)[12]將基于遺傳算法的徑向基函數(shù)（GA-RBF）神經(jīng)網(wǎng)絡(luò)應(yīng)用于磷酸鐵鋰電池SOC 的估算，提高了算法的自適應(yīng)能力。文獻(xiàn)[13]提出一種基于擴(kuò)展Kalman 濾波（EKF）的算法對SOC 進(jìn)行估算，實驗表明，該方法提高了SOC 估算精度。

上述方法分別從不同角度（電池估算模型、測量參數(shù)、控制策略等）對動力電池SOC 估算進(jìn)行了改進(jìn)，而本文主要著眼于鋰電池SOC 估算影響因素的關(guān)聯(lián)性及復(fù)雜性。針對BP 神經(jīng)網(wǎng)絡(luò)算法的非線性處理能力和螢火蟲算法的全局尋優(yōu)能力，并充分考慮鋰電池端電壓和放電電流對算法的影響，本文提出了一種基于螢火蟲算法（firefly algorithm，F(xiàn)A）和BP 神經(jīng)網(wǎng)絡(luò)算法相結(jié)合的方法應(yīng)用于鋰電池SOC 估算。充分發(fā)揮二者的優(yōu)勢，基于FA 對BP神經(jīng)網(wǎng)絡(luò)的權(quán)值和閾值進(jìn)行優(yōu)化，然后利用BP 神經(jīng)網(wǎng)絡(luò)找到其最優(yōu)解。以磷酸鐵鋰電池為測試對象，在ARBIN 公司生產(chǎn)的EVTS 電動車動力電池測試系統(tǒng)裝置上進(jìn)行測試，收集鋰電池的各項性能參數(shù)。仿真實驗結(jié)果表明，與現(xiàn)有的BP 神經(jīng)網(wǎng)絡(luò)估算方法相比，基于FA-BP 神經(jīng)網(wǎng)絡(luò)的鋰電池SOC 估算方法準(zhǔn)確度高，具備很好的實用性。

1 基于FA-BP神經(jīng)網(wǎng)絡(luò)的SOC估算模型

1.1 螢火蟲算法的基本原理

螢火蟲算法[14]于2008年由劍橋大學(xué)Xin-She YANG 首次在文獻(xiàn)中提出，作為群集智能優(yōu)化算法領(lǐng)域的最新算法，該算法具有更好的尋優(yōu)效果，且更易于工程實現(xiàn)。近些年專家學(xué)者的研究表明，螢火蟲算法在全局尋優(yōu)能力上比遺傳算法（genetic algorithm，GA）及粒子群算法（particle swarm optimtzation，PSO）等更有效，成功率更高[15-17]。螢火蟲算法是基于生態(tài)環(huán)境下的螢火蟲通過發(fā)光相互吸引，進(jìn)而實現(xiàn)位置的更迭的自然現(xiàn)象演變而來。

1.2 螢火蟲算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)

BP神經(jīng)網(wǎng)絡(luò)具有很強(qiáng)的并行處理能力以及自適應(yīng)性，因此非常適用于解決復(fù)雜的非線性問題。然而，BP 神經(jīng)網(wǎng)絡(luò)本質(zhì)上是一個具有多極點(diǎn)的非線性函數(shù)，因此在運(yùn)行過程中BP 神經(jīng)網(wǎng)絡(luò)極有可能陷入這些局部極小值。為提高BP 神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)能力及識別精度，基于螢火蟲算法對BP 神經(jīng)網(wǎng)絡(luò)的部分參數(shù)進(jìn)行優(yōu)化。具體步驟如下。

（1）編碼設(shè)置

BP 神經(jīng)網(wǎng)絡(luò)由輸入層N、隱含層H 和輸出層O組成。其中，網(wǎng)絡(luò)中的權(quán)值和閾值可以表示為ωnh、ωHO和θH、θO。在FA-BP 神經(jīng)網(wǎng)絡(luò)訓(xùn)練時，螢火蟲個體采用實數(shù)編碼的方式。因此，F(xiàn)A-BP 神經(jīng)網(wǎng)絡(luò)亟待優(yōu)化的參數(shù)L可表示為

螢火蟲種群間位置向量公式為

（2）參數(shù)初始化

構(gòu)建BP 神經(jīng)網(wǎng)絡(luò)的網(wǎng)絡(luò)參數(shù)，包括訓(xùn)練次數(shù)、初始權(quán)值、訓(xùn)練誤差和初始閾值等。設(shè)置螢火蟲算法的初始值，包括螢火蟲數(shù)目n、隨機(jī)步長α、最大吸引度因子β0和最大迭代次數(shù)T等。算法開始前，在可行區(qū)域內(nèi)隨機(jī)初始所有個體位置。

（3）適應(yīng)度函數(shù)計算

計算FA-BP 神經(jīng)網(wǎng)絡(luò)的目標(biāo)函數(shù)，也就是計算螢火蟲個體的適應(yīng)度值。FA-BP 神經(jīng)網(wǎng)絡(luò)訓(xùn)練過程中，每一個螢火蟲即代表著一個神經(jīng)網(wǎng)絡(luò)。FA-BP神經(jīng)網(wǎng)絡(luò)的訓(xùn)練精度由均方誤差函數(shù)E確定，即

式中，E為FA-BP 神經(jīng)網(wǎng)絡(luò)的均方誤差；Ti為第i個節(jié)點(diǎn)的輸出值；Ki為第i個節(jié)點(diǎn)的目標(biāo)值；N為輸出的個數(shù)。

螢火蟲個體的適應(yīng)度值，即目標(biāo)函數(shù)f為

（4）確定螢火蟲個體間吸引度和相對亮度

種群之間相互的吸引力β的數(shù)學(xué)公式可表示為

式中，r為任意兩個螢火蟲間的距離；β0為r=0的螢火蟲最大吸引力；γ為光照強(qiáng)度的吸收系數(shù)。

任意兩螢火蟲間滿足如下數(shù)學(xué)公式

其中rij為螢火蟲i和螢火蟲j之間的距離，d為維度。

（5）位置更新

根據(jù)螢火蟲個體間的適應(yīng)度值，進(jìn)行位置移動更新。

螢火蟲i的發(fā)光亮度遠(yuǎn)遠(yuǎn)高于螢火蟲j的發(fā)光亮度，則螢火蟲j的位置可做如下更新：

式中，α為（0，1）之間的隨機(jī)數(shù)；Xj為螢火蟲j更新后的位置。

在FA-BP 神經(jīng)網(wǎng)絡(luò)的訓(xùn)練過程中，為了避免螢火蟲間的不規(guī)則運(yùn)動，種群間的移動范圍擴(kuò)大導(dǎo)致螢火蟲間的有效距離增加，進(jìn)而造成種群尋優(yōu)轉(zhuǎn)變?yōu)殡S機(jī)移動的結(jié)果。本文在螢火蟲位置更迭的公式中添加吸引力下限βmin，提高算法的搜索能力，即

（6）條件判斷

迭代次數(shù)達(dá)到最大或訓(xùn)練誤差收斂，結(jié)束，否則返回上一步。

（7）確定權(quán)值與閾值

將螢火蟲最優(yōu)位置向量映射給BP 神經(jīng)網(wǎng)絡(luò)，作為BP 神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值。

（8）BP 神經(jīng)網(wǎng)絡(luò)訓(xùn)練

訓(xùn)練結(jié)束，得到適應(yīng)度值最大個體。將測試樣本輸入，進(jìn)行電動汽車SOC 估算。

2 電池測試系統(tǒng)及數(shù)據(jù)選取

2.1 電池測試系統(tǒng)

以磷酸鐵鋰電池為測試對象，在ARBIN 公司生產(chǎn)的EVTS 電動車動力電池測試系統(tǒng)裝置上進(jìn)行測試，由與其連接的上位機(jī)保存及讀取數(shù)據(jù)。EVTS電動車動力電池測試系統(tǒng)是大功率自動電池測試系統(tǒng)，專門用于電動車電池或混合電動車電池的研究測試，可以采集電池的各項性能參數(shù)（電池的電壓、電流、SOC 等）。動力電池測試系統(tǒng)如圖1 所示，磷酸鐵鋰電池技術(shù)參數(shù)如表1 所示。

圖1 電池測試系統(tǒng)Fig.1 Battery testing system

表1 磷酸鐵鋰電池技術(shù)參數(shù)Table 1 Technical parameters of lithium iron phosphate battery

2.2 鋰電池數(shù)據(jù)庫的建立

根據(jù)EVTS 電動車動力電池測試系統(tǒng)進(jìn)行不同恒流放電模式的鋰電池測試，保持測試的環(huán)境溫度為25 ℃。采取3 種不同的放電倍率（1 C、2 C、3 C）進(jìn)行恒流放電，在測試過程中設(shè)備會自動記錄鋰電池各個參數(shù)的變化，采樣頻率為1 Hz。訓(xùn)練與測試樣本數(shù)據(jù)如表2 所示。

表2 訓(xùn)練與測試樣本數(shù)據(jù)Table 2 Training and testing sample data

3 仿真與結(jié)果分析

本文采用螢火蟲算法和BP 神經(jīng)網(wǎng)絡(luò)相結(jié)合的方法進(jìn)行電動汽車鋰電池SOC 估算。螢火蟲算法參數(shù)設(shè)置：最大吸引度因子β=0.08，最大迭代次數(shù)T=100，螢火蟲算法的種群大小n=50，初始隨機(jī)步長α=0.03，光吸收系數(shù)γ=0.6。BP 神經(jīng)網(wǎng)絡(luò)算法參數(shù)設(shè)置：訓(xùn)練次數(shù)N=1000，學(xué)習(xí)率η=0.1，網(wǎng)絡(luò)誤差E=0.00001。

為了測試網(wǎng)絡(luò)的性能，將表2 中的測試樣本輸入網(wǎng)絡(luò)，其結(jié)果由圖2 所示。由圖可見，電池SOC網(wǎng)絡(luò)估算結(jié)果和實驗測量結(jié)果非常吻合，其絕對誤差平均值為1.01%，表明FA-BP 神經(jīng)網(wǎng)絡(luò)具有較高的準(zhǔn)確性。

圖2 鋰電池SOC 網(wǎng)絡(luò)估算結(jié)果與實驗測量結(jié)果Fig.2 Estimated and measured results of SOC network for lithium batteries

為了進(jìn)一步驗證FA-BP 神經(jīng)網(wǎng)絡(luò)的實用性和泛化能力，將訓(xùn)練好的網(wǎng)絡(luò)模型編程嵌入到EVTS電動車動力電池測試系統(tǒng)中，進(jìn)行不同放電倍率下的放電測試。圖3 為隨機(jī)選取的放電倍率為0.5 C情況下磷酸鐵鋰電池放電電壓與SOC 變化曲線。其實際測試SOC 值與FA-BP 神經(jīng)網(wǎng)絡(luò)估算SOC 值的誤差在1.54%左右，而與BP 神經(jīng)網(wǎng)絡(luò)估算SOC 值的誤差范圍在2.78%左右，表明基于FA-BP 的動力電池SOC 估算具有良好的估算效果。

圖3 鋰電池放電電壓與SOC 關(guān)系圖Fig.3 Relation between discharge voltage and SOC

4 結(jié) 論

（1）本文建立了基于FA-BP 神經(jīng)網(wǎng)絡(luò)的鋰電池SOC 估算模型，并應(yīng)用于磷酸鐵鋰電池SOC 的估算中。仿真結(jié)果表明，基于FA-BP 神經(jīng)網(wǎng)絡(luò)的SOC估算方法可行，比BP 神經(jīng)網(wǎng)絡(luò)估算的結(jié)果更加精確，實現(xiàn)了電動汽車SOC 的最優(yōu)估算。

（2）由于實驗條件的影響，本文僅考慮了放電電流和工作電壓對電動汽車SOC 估算模型的影響，人為排除了電池壽命、室內(nèi)溫度等因素。下一步的工作是考慮多種因素下SOC 估算模型的建立，實現(xiàn)復(fù)雜工況下的電動汽車SOC 估算。