蓄電池荷電狀態(tài)閉環(huán)動態(tài)估算模型

目前，由于鉛酸蓄電池的經濟性和技術成熟性，使其成為豐要的儲能設備。為了達到優(yōu)化蓄電池電力系統(tǒng)效率的目的，對蓄電池容量的實時監(jiān)控必不可少。而由于蓄電池的非線性特性，反映其容量的關鍵參數荷電狀態(tài)(SOC)，作為電池的內特性不可能直接進行測量。SOC數值只能使用工作電壓、電流等直接測量得到的外特性參數估算獲得。

本文使用最優(yōu)估計理論建立蓄電池的動態(tài)工作模型，實現蓄電池SOC的實時估算。該動態(tài)模型被劃分為兩個部分：第一部分是蓄電池數學解析描述，即對蓄電池工作特性的開環(huán)描述；第二部分是加入動態(tài)過程的描述，實現蓄電池工作特性的閉環(huán)描述。對于蓄電池的解析模型，較為通用的方式是建立描述輸入輸出之間關系的數學模型，通過實驗來確定模型的某些參數，或者模型內部的某些狀態(tài)量。然而，僅僅使用開環(huán)描述模型得到動態(tài)輸出與實際的動態(tài)情況常常存在偏差，這種誤差主要歸咎于測量過程中的異常偏差。當這種誤差出現時，只有閉環(huán)描述模型才能根據這些誤差對模型進行調整。本文使用基于電化學理論的安時模型實現電池數學解析描述，而動態(tài)過程描述則使用帶有自矯正能力的擴展卡爾曼濾波算法。

1 基于電化學的安時模型

普通的安時計量法使用下式估算蓄電池的SOC。

式中：s(0)為初始時刻的蓄電池SOC數值，若從充滿開始放電，其值可以設為1；s(t)為t時刻的SOC實時值；Q為蓄電池的標稱容量；η為庫侖因子。通過調整庫侖因子可以滿足不同放電電流下的SOC計算。實際應用中，庫侖因子多通過試驗確定為常數或是關于放電電流I的函數。但是，蓄電池的標稱容量不等于實際容量，且實際容量在使用中也會衰減。同時，確定庫侖因子過程中產生的誤差，也會影響到安時估算的精度。為了對上述問題進行改進，提高安時法 SOC估算的精度。本文使用電化學理論，構造新的基于安時法的SOC估算模型。

1.1 電解液活性物質濃度損失函數

蓄電池內部電解液所含有的活性物質，其濃度損失百分比可以表示為：

式中：C*為初始濃度；C(t)為電解液中t時刻活性物質的濃度；時間t的取值范圍[0，L]，L為放電總時間。

當使用蓄電池一維的電化學模型，根據電化動力學理論，最終可以得到電解液活性物質濃度損失百分比函數：

式中：v為反應中電子的數目；F為法拉利常數；A為電極的面積；D為擴散系數。

1.2 電化學安時模型

由于電解液的活性物質濃度和電池的SOC成正比的關系，設比例系數為M，可以直接得出電池t時刻的SOC解析表達式：

若考慮電流值為I的恒流放電過程，放電截止時ρ(L)1，則可以得到以下等式：

對于給定的恒流放電集合{I*，*=1，2，…，n}，可以使用最小二乘法得到最優(yōu)的α、β參數，其中：

得到模型參數之后，為方便模型的實際應用，使用積分的矩形近似方法改寫(4)式，用以獲得離散時間上的近似遞推模型，在間隔周期△t足夠小的情況下，遞推模型可以寫為：

式中：sk表示k時刻的電池SOC的實時值；Ik表示k時刻的電池電流。對比式(1)的標準安時估算模型，可以發(fā)現α等于電池的標稱容量Q，庫侖因子則由β 和放電時間k△t決定。從電化學角度分析，表達式(7)的括號中的第二項表示蓄電池中無法使用的總電量，當β數值增加的時候，第二項趨向于零。因此，較大的β數值意味著蓄電池可以被看作理想儲能元件，所有充電電量都可以完全通過放電過程釋放。這是因為大的β數值表明更快的擴散效應，蓄電池電解液中的活性物質可以更快的到達電極的表面。反之，小的β數值表明蓄電池儲能損失大，大量的充電電量無法在放電過程中釋放。

2 擴展卡爾曼濾波閉環(huán)估算模型

改進安時模型能夠較好地反應電池的動態(tài)特性，但這種蓄電池SOC計算方式只是一種開路的估算方式，存在著傳統(tǒng)安時計量法的缺點，即對電流測量中的測量偏差十分敏感，某一個時刻出現的測量偏差，可以影響到該時刻后所有的SOC估算值。如果將估算模型構造成閉環(huán)反饋的模式，則可以自動修正電流測量中的偏差，給出正確的SOC估算值。在(7)式遞推模型的基礎上，可以使用卡爾曼濾波器方法構造出具有閉環(huán)特性的電池SOC估算模型。

首先將(7)式作為蓄電池SOC估算系統(tǒng)的狀態(tài)方程，蓄電池SOC為狀態(tài)量，蓄電池的工作電流作為系統(tǒng)的輸入。然后，利用蓄電池的工作電壓構造系統(tǒng)的觀測方程。

蓄電池負載電壓與當前時刻蓄電池的開路電壓(Vcc)之間的關系是：

式中：R為蓄電池內阻。又由于Vcc和內阻都與其SOC有著直接的關系，故可以使用關于sk的函數，得到卡爾曼濾波算法中的觀測方程：

式中：uk表示k時刻的電池端電壓，則(7)式和(9)式組成了蓄電池SOC估算的卡爾曼濾波系統(tǒng)。確定(9)式的具體過程將在實驗部分詳細分析。