胡 晨,金 翼,崔邴晗,杜春雨?
(1.中國電力科學研究院有限公司新能源與儲能運行控制國家重點實驗室,北京 100192;2.哈爾濱工業大學化工與化學學院,特種化學電源研究所,黑龍江哈爾濱 150001)
人們針對閥控式鉛酸(VRLA)電池健康狀態(SOH)在線估算開展的研究,所用的電化學模型主要可分為電化學阻抗模型、經驗模型及數據驅動模型[1]。電化學阻抗模型通過對電池進行阻抗測試,獲取隨電池老化而變化的內部關鍵參數(如歐姆阻抗、電荷傳遞阻抗和雙電層電容等),來反映電池內部電化學反應過程。梅成林等[2]引入小波神經網絡對電池壽命進行估算,預測平均誤差為1.49%,但實驗采集的樣本量較少,泛化性可能較差。經驗模型通過對電池測試數據進行擬合并合理外推,獲得電池的SOH[3],但測試數據無法反映多工況條件下的狀態,不能應用于動態工況場景。數據驅動模型將電池視為黑箱,無須了解內部機理,利用智能算法實現對SOH的預測。Z.Chen等[4]利用隨機森林算法,通過電壓與電流平均值求出電池SOH,誤差低于2%;M.Talha等[5]利用人工神經網絡算法估算電池的開路電壓,通過等效電路模型預測電池荷電狀態(SOC),再通過SOH與SOC的關系計算出SOH,SOC的估算誤差為1%~2%,但由于SOH是通過SOC來估算的,誤差會進一步增大,造成精度降低。數據驅動模型具有較高的靈活性與估算精度,且計算量小,便于實際應用,但不對電池的衰減老化特性進行分析,僅依靠優化數據驅動模型算法,會導致數據驅動模型泛化能力差、精度無法滿足要求。
簡單的電壓電流數據無法完全反映鉛酸電池的衰減信息,可將對衰減敏感性高的電壓增量作為特征參數。鉛酸電池具有特定的衰減規律,且深度神經網絡(DNN)模型精度高,只需經過訓練便可使用,同時可在線學習,便于實際應用。本文作者將電壓增量引入基于深度學習的DNN模型中,提出以電壓增量為特征參數的DNN模型。將反映電池內部衰減信息的宏觀特征引入到DNN模型中,增強模型的泛化性,利用DNN模型的計算能力,增強SOH的估算精度,以實現優勢互補。
SOH可以反映出電池的可用容量、內阻的變化與老化信息,是描述電池衰減程度的指標。通過電池容量或內阻的變化,反映電池SOH的計算公式,分別見式(1)、(2)。

式(1)、(2)中:SOH為電池的SOH;Q為電池的容量;R為電池的阻抗;下標now、aged和full分別代表相應物理量的當前值、衰減到壽命終點處的值和額定值。
鉛酸電池在實際應用中,常常無法全充全放,僅以某一放電深度(DOD)進行放電和充電,上述兩種SOH計算方法在不進行容量標定時,難以計算以DOD為條件的SOH,因此,提出一種基于電池放電截止電壓估算SOH的方法:

式(3)中:U為電池的端電壓;下標DOD為放電到某DOD時的值;U’DOD為首次循環放電到某DOD時的端電壓。
選取電壓、電流、SOC及電池的微分電壓作為模型的特征值,微分電壓描述的是一段時間(t)或SOC(SOC)變化對應的電壓(U)增加或減少量(ΔU/Δt,ΔU/ΔSOC)。 通過對鉛酸電池的充放電數據進行微分,可得到電池的電壓增量曲線,以便監測隨著電池衰減,充放電曲線的微小變化[6]。電壓增量曲線與充放電曲線相比,對工況改變和電池衰減具有更高的敏感性,因此構造以電壓增量作為特征參數的DNN模型,可以提高SOH的估算精度。
DNN由輸入層、隱藏層及輸出層構成,網絡結構見圖1。

圖1 DNN的結構圖Fig.1 Schematic of deep neural network(DNN)
所構建的DNN模型,包括4層隱藏層。特征值(U、I、SOC、dU/dt、dU/dSOC)通過組合排列,形成訓練集X;SOH估算公式得到對應X條件下的電池SOH,作為訓練標簽Y;每個神經元之間的權重參數矩陣為w;神經元之間的偏置矩陣為b;訓練樣本數為m;選取的特征值數目為n;隱藏層神經元數目分別為n1、n2、n3和n4;隱藏層中神經元所使用的激活函數為σ;z為隱藏單元中計算結果。用ReLu激活函數對z進行運算,以增強神經網絡對數據的非線性處理及擬合能力,得到估算值a,a為上一層激活函數值。通過計算a與真實值y之間的差值,得到損失值L。DNN模型由前向傳播算法與反向傳播算法構成,首先建立前向傳播算法數學模型:

均方誤差(MSE)對異常點的魯棒性差,平均絕對誤差(MAE)在最優點附近梯度較大,難以收斂,Huber損失函數二階微分不一定可導。相對而言,雙曲余弦對數損失函數對于較小的x近似等于MSE,對于較大的x近似等于MAE,對異常點的魯棒性強,梯度隨L的減小而減小,且二階導數處處可微。得到損失值后,可通過反向傳播算法對權重值w及偏置值b進行更新,建立反向傳播數學模型:

式(7)對所有樣本的Loss進行求和后,得到代價函數J,從輸出層開始反向求導;式(8)為第4層隱藏層權重值w的更新梯度;式(9)為第4層隱藏層偏置值b的更新梯度。以此類推,得到所有層中權重值w及偏置值b的更新梯度,之后對w及b進行更新,輸出層得到新預測值a,模型通過最小化L值以得到最佳的w、b及a,即完成訓練。
對于包含多個隱藏層的深度學習問題,在訓練過程中,隱藏層的數據分布由權重值w及偏置值b的變化而不斷變化,最終將數據分布引入到非線性激活函數取值區間,造成在反向傳播過程中較前隱藏層的梯度較小甚至消失,訓練速度慢。對隱藏層的輸出進行歸一化,就能將輸出值分布在非線性激活函數對輸入比較敏感的區域,得到較大的反向傳播更新梯度,避免產生梯度消失及局部收斂的問題,加快神經網絡的訓練速度。
建立的DNN模型包含1層輸入層、4層隱藏層、5層BN層和1層輸出層。輸入特征值矩陣X后,經過前向傳播算法得到損失值L,通過批標準化(BN)層調控輸出值的數據分布,再經過反向傳播算法更新權重值及偏置值,以最小化L,最終得到鉛酸電池SOH估算的最佳精度。
實驗平臺由8-GCMH型鉛酸電池(江蘇產)、5 V/5 A充放電測試系統(武漢產)和中央控制電腦組成。在25℃的各種工況條件下對電池進行充放電,并實時記錄工作電壓、工作電流、SOC及測試時間等數據。實驗電池的額定容量為20 Ah,充放電電壓為1.75~2.40 V,定義電池容量為額定容量的80%達到壽命終點。測試電池編號及工況列于表1。

表1 電池編號及電池測試工況Table 1 Battery numbers and test conditions
1號、2號和3號電池的老化測試為全充全放:先將電池放電到1.75 V,靜置1 h后,以表1中倍率恒流充電至充電截止電壓2.4 V,之后恒壓充電至截止電流0.84 A。再次靜置1 h后,放電至1.75 V,獲取電池可用容量。記錄初始電池容量為Q1,第i次循環電池容量為Qi。4號、5號電池為不同DOD的老化測試:先將電池放電到1.75 V,靜置1 h后,以標準的恒流恒壓充電法充電到2.40 V,再次靜置1 h后,通過控制放電容量將電池放電至不同的DOD,記錄初始截止電壓U1,第i次放電截止電壓Ui,將電池不斷循環,最終的放電截止電壓為3號電池達到壽命終點時的工作電壓。
1號電池測試工況下的電壓電流采樣曲線見圖2。

圖2 1號電池測試工況下的電壓電流采樣曲線Fig.2 Voltage and current sampling curves of No.1 battery under test conditions
為驗證以對電池衰減高度敏感的電壓增量作為特征參數的DNN模型,對鉛酸電池SOH估算精度的優越性,先對測試的鉛酸電池數據進行處理,得到包含電池電壓電流數據的特征值矩陣X1及包含電壓增量信息的特征值矩陣X2,構建對特征值矩陣X1及X2進行非線性擬合的DNN模型。
對以電壓增量為特征參數的DNN模型得到的不同循環倍率下使用容量計算的SOH進行驗證,結果如圖3所示,估算誤差(E)如圖4所示。

圖3 以電壓增量為特征參數的DNN模型得到的以容量為計算方法、循環倍率為影響因素的電池SOH估算值Fig.3 Estimated value of battery SOH obtained by DNNmodelwith voltage increment as the characteristic parameter,capacity as the calculation method and cycle rate as the influencing factor

圖4 DNN模型得到的電池SOH估算誤差 Fig.4 Estimation error of battery SOH obtained by DNNmodel
從圖3可知,以不同倍率進行充放電循環,鉛酸電池SOH與循環次數近似成線性關系,說明在循環過程中,容量的衰減率基本不變;當循環倍率升高時,衰減速度減慢、壽命延長。
從圖4可知,不同循環倍率條件下的鉛酸電池在全循環壽命周期內的SOH估算精度頗高,1號電池在全周期循環中誤差均小于1.0%,估算精度最高;2號電池的最大誤差來源于循環末期,約為3.0%,其他循環誤差在1.0%左右;3號電池的最大誤差同樣位于循環末期,最大誤差約為1.2%。
使用電壓、電流數據的數據驅動模型對不同循環倍率下使用容量計算得到的SOH進行驗證,結果如圖5所示,估算誤差如圖6所示。

圖5 使用電壓、電流數據的數據驅動模型得到的以容量為計算方法、循環倍率為影響因素的電池SOH估算值Fig.5 Estimated value of battery SOH obtained by data-drivenmodel of voltage and currentdatawith capacity as the calculationmethod and cycle rate as the influencing factor

圖6 數據驅動模型得到的電池SOH估算誤差 Fig.6 Estimation error of battery SOH obtained by data-driven model
從圖5可知,SOH估算值與真實值之間的誤差相比于圖3要大,說明特征值的選取對SOH的估算精度有很大影響。從圖6可知,使用電壓、電流數據的數據驅動模型對不同循環倍率下的鉛酸電池SOH估算精度,相比于以對電池衰減高度敏感的電壓增量作為特征參數的DNN模型有一定的差距,最大誤差均在5.0%左右,約為DNN模型的5倍。
對以電壓增量為特征參數的DNN模型得到的不同循環倍率下使用放電截止電壓計算的SOH進行驗證,結果見圖7。

圖7 以電壓增量為特征參數的DNN模型得到的以放電截止電壓為計算方法、DOD為影響因素的電池SOH估算值Fig.7 Estimated value ofbattery SOH obtained by DNNmodelwith voltage incrementas the characteristic parameter,discharge cut-off voltage as the calculationmethod and DOD as the influencing factor
從圖7可知,在DOD為80%時,與全充放結果相同,SOH和循環次數近似呈線性關系。當DOD減小到60%時,SOH與循環次數的關系發生了變化,即隨著循環的進行,電池的衰減速度逐漸加快。將放電截止電壓用于SOH的計算,電池的SOH估算擬合誤差(EMAE)均小于0.5%,相比于以容量計算的SOH,誤差得到減小,說明使用放電截止電壓對SOH進行計算,可對模型進行優化,提高性能。
對使用電壓、電流數據的數據驅動模型得到的不同DOD條件下使用放電截止電壓計算的SOH進行驗證,結果見圖8。

圖8 使用電壓、電流數據的數據驅動模型得到的以放電截止電壓為計算方法、DOD為影響因素的電池SOH估算值Fig.8 Estimated value of battery SOH obtained by data-drivenmodel of voltage and current data with discharge cut-off voltage as the calculation method and DOD as the influencing factor
從圖8可知,相比于以循環倍率為影響因素的3只電池,當以DOD為影響因素時,模型誤差增大。以80%DOD進行充放電,電池SOH估算的最大誤差在25.0%左右,以60%DOD進行充放電,最大誤差在40.0%左右,相比于以對電池衰減高度敏感的電壓增量作為特征參數的DNN模型,模型誤差增大約80倍。因此,DNN模型可對模型進行優化,減少模型誤差,提高估算精度,增強模型性能。
使用均方根誤差RMSE、MSE及MAE[4]對5只電池的誤差進行分析。電池誤差分析結果見表2。

表2 電池誤差分析結果Table 2 Analysis results of battery error
從表2可知,80%DOD、0.2C倍率循環下,電池的誤差相對較高,原因是在電池循環終點處,電池SOH的估算值與真實值相差較大。盡管如此,4號電池的估算誤差仍在1.0%內,說明模型對于鉛酸電池SOH預測的準確性較高。使用對電池衰減高度敏感的電壓增量作為特征參數的DNN模型,與應用電壓、電流的數據驅動模型相比,前者平均誤差在0.2%左右,后者平均誤差在1.4%左右。使用對電池衰減高度敏感的電壓增量作為特征參數的DNN模型誤差為后者的1/7。相比于數據驅動模型,DNN模型具有更少的異常點,且異常點造成的誤差(5.0%)遠小于數據驅動模型(40.0%)。說明與數據驅動模型相比,使用對電池衰減高度敏感的電壓增量作為特征參數的DNN模型具有更強的魯棒性。
針對鉛酸電池SOH在線估算問題,本文作者提出使用對電池衰減高度敏感的電壓增量作為特征參數的DNN模型的鉛酸電池SOH估算方法,選取對電池衰減老化更加敏感的電壓增量參數作為模型特征值,并進行全壽命周期不同倍率(0.10C、0.15C和0.20C),不同 DOD(60%、80%和100%)的鉛酸電池老化實驗對模型精度進行驗證。得到以下結論:
鉛酸電池在進行全充放循環時,SOH與循環次數成線性關系,隨著循環倍率的增大,電池循環壽命延長。
使用DNN模型以循環倍率作為影響因素估算電池的SOH,模型誤差減少至0.2%;以DOD作為影響因素估算電池的SOH,模型誤差減少至5.0%,達到了模型的優化;模型減少了異常點的數目,增強了魯棒性;相比于使用電壓,電流的數據驅動模型,誤差減少為原來的1/7。
使用DNN模型以循環倍率及DOD作為影響因素估算電池的SOH,仍存在一些問題。如在電池循環壽命終期,容量衰減速度較循環前中期加快,導致模型在此階段誤差增大。需要將鉛酸電池在其他測試條件下的循環數據引入模型,增強模型對于循環壽命終期衰減速度加快的適應性,增加模型的估算精度,減少異常點的數目。
對于下一步研究,應適當增加鉛酸電池測試工況,增強模型泛化性,尋找容量最優的工況條件。