基于改进粒子滤波的锂离子电池健康状态估计

徐超 李立伟 杨玉新 （

【摘要】为了准确估计电池当前健康状态（ SOH），提出一种基于改进粒子滤波算法的电池SOH在线评估方法。首先针对传统萤火虫优化算法的不足，提出了一种改进萤火虫算法替代传统粒子滤波的重采样，然后从锂离子电池工作时的可测参数中提取在线健康指标（ HI），建立HI与SOH之间的映射模型，并将其应用于状态空间模型的观测。利用马里兰大学先进寿命周期工程中心（ CALCE）公布的试验测量数据进行验证，结果表明，该方法对具有非线性和非高斯特性的锂离子电池降解过程的状态估计具有良好的适应性。

徐超1 李立伟2 杨玉新3

1. 266071；2. 266071；3. （ 青岛大学，电气工程学院，青岛 青岛大学，威海创新研究院，青岛 青岛大学，图书馆， 266071）青岛

主题词：锂离子电池 SOH估计 粒子滤波算法 萤火虫算法 健康指标TM912 A 10.19620/j.cnki.1000-3703.20191421中图分类号： 文献标识码： DOI:

State of Health Estimation of Lithium-Ion Battery Based on Improved Particle Filter

Xu Chao1, Li Liwei2, Yang Yuxin3 1. School of Electrical Engineering, Qingdao University, Qingdao 266071; 2. Weihai Innovation Institute, Qingdao （ University, Qingdao 266071; 3. Library, Qingdao University, Qingdao 266071） Abstract The State Of Health (SOH) estimation is one of the key functions of lithium- ion Battery Management【 】System (BMS). In order to accurately estimate the current health state of battery, this paper proposes an improved firefly algorithm and replaces the resampling of the traditional particle filter in view of the shortcomings of the traditional firefly algorithm. And then the measurable parameters of lithium- ion battery are extracted on line health index Hi, the mapping model between Health Index (HI) and battery SOH is established and applied to the observation of state space model. An on- line evaluation method of battery SOH based on improved particle filter algorithm is proposed. Finally, based on a set of experimental data of battery capacity published by the University of Maryland Center for Advanced Life Cycle Engineering (CALCE), it is proved that the method has good adaptability to the degradation process of lithium- ion batteries with nonlinear and non-Gaussian characteristics.

Key words: Lithium- ion battery, SOH estimation, Particle filter algorithm, Firefly algorithm, Health indicators 1 前言

锂离子电池具有自放电率低、能量密度高等显著优点[1]，在交通、航天、国防等领域得到了广泛应用[2]。但在反复充、放电过程中，锂离子电池容量会下降，性能会降低，严重时甚至会发生灾难性事件[3- 4]。因此，准确估

State Of Health，SOH）

计锂离子电池的健康状态（ 对于电池的安全可靠运行具有重要意义。

SOH

电池容量是电池 的直接指标，在许多估计方

Chen

法中都有应用。 等人[5]提出了一种基于欧姆内阻

SOH

与容量衰减关系的电池 估计方法。与此同时，基

Health Index，HI）

于可测电池参数的健康指标（ 在线构建也得到了越来越多的研究[6- 8]。此外还有基于数据驱动的预测方法，但是需要与电池电容和可测量的物理参数建立映射关系，对训练数据集具有高度依赖性。电池在复杂条件下工作时，该方法可能会降低精度。因此，基于模型的预测方法越来越受到人们的重视，该方法通过提取表征电池动态老化和失效过程的内部参数来建立物理模型。由于模型参数随电池退化而变化，因此可以采用滤波方法递归更新模型参数[9-14]。

Kalman Filter，KF）

然而，卡尔曼滤波（ 存在一些固有的缺陷，不适应非线性的模型，而且要求系统噪声和测量噪声必须满足高斯分布，但锂电池退化过程并不符

Particle Filter，PF）

合这些要求。相比之下，粒子滤波（

SOH

算法在处理非高斯和非线性电池 方面具有优势[15]，

PF

因而受到许多学者的关注和重视。然而， 算法存在粒子多样性退化的问题，为解决这一问题，可以采取的

2

主要策略有 种：一是改进重采样技术[16-17]，虽然重采样技术可以在一定程度上解决粒子数量不足的问题，但是会导致粒子的多样性缺失；二是选择更合理的粒子建议分布[18]。近年来，随着生物智能算法的出现和发展，采用智能群体优化算法优化粒子滤波来解决粒子退化贫化问题具有非常好的应用价值[19- 21]。因此，本文选用萤火虫算法优化粒子滤波，并对传统萤火虫算法加以改进，模拟萤火虫优化思想的吸引行为和移动行为代替粒子滤波重采样来解决粒子贫化问题。

2 基本算法2.1 粒子滤波算法

粒子滤波是一种基于蒙特卡洛仿真的近似贝叶斯滤波的概率统计算法[22]。许多非线性滤波问题可以用动态空间模型来描述，状态方程和观测方程分别为：

式中， x 、y 分别为 k时刻的状态值和观测值； f 为状k k k态转移函数； h 为观测函数； u 为过程噪声； v 为观测k k k噪声。

设初始的状态概率密度为P( x0|y0)= P( x0) ，则得到

( k- 1) )

时刻状态变量xk后验概率密度分布P( x |y1: 1 的状k k态预测方程为：

2.2 萤火虫算法

通过对热带萤火虫的发光机制和运动行为进行模

Yang

拟， 提出了一种新型的自然启发式群体智能优化算法，即萤火虫算法（ Firefly Algorithm，FA）[ 23]。真实的萤火虫发光自然地呈现出一种离散的闪烁模式，而萤火虫算法假设它们总是在发光[24]。

2.2.1 荧光亮度荧光亮度取决于萤火虫所在位置的适应度，并与其近似成正比。亮度越高，表示该个体的吸引力越强，其他个体向其移动的概率越大，并且存在亮度随着距离而弱化的过程，定义荧光亮度参数为： - γr2 ij

式中， I0 为 r= 0时的光强，即最大光强； γ为介质的吸ij

1；

收系数，通常取 r 为萤火虫i和 j 之间的欧几里得ij距离：

式中， x 和x 分别为萤火虫i和 j所处的空间位置； D为i j维度； x 和x 分别为萤火虫i和 j 所在的第d维分量。id jd 2.2.2

吸引度萤火虫的吸引度随荧光亮度的变化而变化，定义吸引度参数为：

式中， β0为最亮萤火虫的吸引度。

2.2.3 位置移动萤火虫在受到比自身更亮的个体吸引时，会向更亮的个体移动，位置的更新表示为：

度； α( R- 0.5)为随机项，其中α为步长因子， R ∈[0,1] 为随机数。

3 改进算法3.1 改进萤火虫算法

FA

是一种参数少、稳定性强，且操作方便的算法，

但是在迭代后期会出现随萤火虫之间距离逐渐缩小，其相对吸引度增大的现象，从而导致萤火虫算法局部搜索能力变弱，并且有可能在峰值附近反复振荡，使得算法的计算速度和求解精度下降。与此同时，萤火虫的个体多样性变差，易出现早熟收敛现象。针对传统萤火虫算

FA

法的不足之处，本文对 进行改进。

3.1.1

采用动态搜索步长

FA

在标准 中， α是固定的，未考虑随着迭代次数增加，萤火虫之间距离逐渐减小的特征。因此，本文采用随迭代变化自适应控制的随机步长因子α来减小萤火虫的随机运动。α在每次迭代中进行动态调谐：

1

式中， tmax 为最大迭代次数。

3.1.2

引入惯性权重因子

13）

虽然式（ 最后的带有特定系数的随机项可以在一定程度上避免局部振荡的延续，但是可能由于迭代次

FA

数过多，精度无法得到满足。因此，为了提高 的局部搜索能力，将惯性权重引入位置更新公式：

式中， ωt为惯性权重因子，一般取0~1范围内的常数。

3.1.3

改进移动策略每一只萤火虫都在逐步向更亮的萤火虫移动，因为每只萤火虫的亮度都需要相互比较，而且每次比较都伴随着一次移动，这就会使得跟踪时间过长。为了克服这一问题，以所有比萤火虫i亮的位置坐标的平均值作为代表点，萤火虫i只向这一代表点移动而不必向其他较亮萤火虫移动，这将大幅缩短跟踪所需时间。最终改进的位置移动公式为：

位置的平均坐标； Z为较亮萤火虫的数量。

3.2 改进萤火虫算法优化粒子滤波步骤

a. 设定所需参数，包括粒子数N、最大吸引度β0、步长因子α、传播媒介的吸收系数γ、最大迭代次数tmax和算法迭代终止阈值ϵ。

b. 粒子初始化。根据先验分布 P( x0) 得到{ xi 0, i= 1,2,…, N}，其中 xi 服从重要性密度函数：

0

c. 采用改进萤火虫算法来优化粒子，模拟萤火虫优化思想的吸引行为和移动行为代替粒子滤波重采样。首先求出 x ，并获取粒子吸引度，然后根据式16）

（ 更新粒子位置，迭代拟合最小方差估计。粒子吸引度计算公式为：

d.

判断拟合终止。当荧光亮度大于设置的终止阈值ϵ时，停止算法迭代，否则继续迭代至最大迭代次数；当算法符合设定的ϵ时，说明粒子已收敛至粒子真实值，或者达到最大迭代次数，此时停止优化；否则转入

步骤c。

e. 权重补偿及更新。萤火虫算法与粒子滤波结

合的核心思想是，为了得到更准确的预测值，对PF中的粒子逐个进行迭代寻优，粒子会向后验概率密度值高的区域聚集。但是萤火虫算法会导致粒子滤波分布存在贝叶斯滤波理论丢失的情况，因此，在粒子位置更新的同时对权重进行补偿及更新：

式中， si为 k时刻的粒子i。k f. 20）计算各新粒子的权重wi ，采用式（ 对新粒子k权重进行归一化：

g. 21），将归一化后的权重代入式（ 计算滤波状态估计值并输出：

4 仿真及结果分析4.1 锂电池数据集

Center for

本文采用先进寿命周期工程中心（

Advanced Life Cycle Engineering，CALCE）

公布的利用美

Arbin BT2000

国 锂电池试验系统对锂电池进行连续的充、放电试验获得的一组电池容量试验测量数据。试验

1.1 A·h

在室温条件下进行，所用锂电池的额定容量为 。

0.675 A 4.2 V

首先，以恒定电流 充电，当电压升高到 时

1.350 A

充电结束。然后，以恒定电流 放电，放电速率为

1 C， 2.7 V

当电压降低到 时放电结束。最后，对锂电池

0.1~5.0 Hz

进行电化学阻抗谱扫描，频率扫描范围 。

4.2 锂离子电池在线SOH估计方法

本文通过提取电池实际运行过程中可以直接监测HI，到的放电电压序列和时间序列构造一个在线的 对

HI SOH

构造的 与电池 之间的相关性进行评估和映射，并将此映射关系应用到系统状态空间方程中以实现电

SOH

池 估计。等放电电压差时间间隔潜在退化模型与电池容量退化模型类似，取等放电电压采样的最大值Vmax和最小值 Vmin分别为4V 3.5 V，

和 一定充、放电周期内，提取等

HI

放电电压时间间隔的在线 为：

| | 1,2,…, 22）

H=t -t , i= k （

i V max V min

式中， tVmax和 tVmin分别为电压采样的上限和下限时间。

HI

在线 序列可以表示为：

HI ={H1, H2,…, H} 23）

（ k

HI SOH HI

在线 可以作为电池 的间接指标，而在线

SOH

与电池 之间的相关分析和转换关系也具有一定的

[25] SOH HI

参考价值。由文献 可知，电池 与在线 之间的转换关系可以描述为： i)+

O = β0 + β1H + β2 ln( 24）

H ε （ i i式中， O 为电池健康状态； β0 为常数； β1 和 β2 为转换i关系系数； ε为误差。这种近似方法应用线性基础展开，是广义线性回归模型的扩展之一，这种转换关系的优点是可以添加平滑函数来描述两个变量之间的线性或非线性关系。

4.2.1

退化模型和状态转移方程锂离子电池反复充、放电造成的退化是一个复杂的电化学反应过程，因此很难建立合适的模型来表征整个降解过程，但是不难看出，电池退化的总体趋势是指数

Xing

衰减的。而且 等人[26]的研究表明，双指数退化模型能够很好地描述电池的退化过程：

exp( exp( dk) 25）

O = a· bk) + c· （ k

式中， a、c为内阻抗； bd、 为退化速率。

1

电池双指数模型所确定的函数如图 所示。由于原数据集样本量较大，先对数据样本进行约减预处理，

5

原数据样本每隔 次循环抽取一个新数据点，构建缩减

R= 0.995 5，

后的数据样本。拟合相关系数 说明双指数模型表征电池降解过程的能力很强。

选择退化模型参数作为系统状态空间变量，状态转移方程可建立为：

; ; ; ] 26）

X =[a b cd （ k k k k k a) ì a =a +υ σ

1 1, υ~ N( 0, k k- a, k - a b)

b =b +υ σ

1 1, υ~ N( 0,

27） í k k- b, k - b c) （

c =c +υ σ

1 1, υ~ N( 0, k k- c, k - c d)

d =d +υ σ

1 1, υ~ N( 0, î k k- d, k - d

N( 0, σm)(m= a,b,c,d)

式中， 为均值和标准差都为零的高斯噪声。

4.2.2

测量方程测量方程描述了系统状态变量与观测信息之间的

HI

函数关系，选择在线作为状态空间方程的测量值。

HI SOH

考虑到电池退化模型以及在线 与电池 之间的映射关系，测量方程为：

exp( ) exp( k)

O =a· bk + · d

k k k N( 0, k 28） ) μ)

H = g(O + μ , μ ~ σ （ k k k k

HI SOH

式中， g()∙为在线 与电池 之间的转换关系。

4.2.3 SOH

状态更新和电池 估计

SOH

在第 k 个循环周期，每个粒子对电池 的估计为：

exp( ) exp( k) 29）

Oi = ai · bi k + ci · di （ k k k k SOH

则估计的电池 表示为： Σ N

30）

O= wi Oi （ k k k 1 i=

4.3 仿真设计

CALCE A12

首先，对 数据集中 系列电池进行在线

HI HI SOH

提取，并对构建的 与 进行相关性分析。然

HI SOH

后，对测量方程建立在线 与电池 之间的映射函数，对状态转移方程建立退化模型。最后，利用改进算

SOH

法对电池 估计的状态参数进行优化。为便于与电

HI

池数据对比，将在线 除以某一基准值，所构造的在线

HI SOH 2 2 HI

和真实 如图 所示，由图 可知，在线 与电池

SOH SOH

有很强的相关性，说明电池 的变化可以通过

HI

提出的 表示。

24） SOH HI 3

由式（ 和电池 获得的映射在线 如图 所示。利用最小二乘法计算系数 β0 、β1 和 β2 。最大误0.045 7，

差为 表明利用映射关系可以建立系统状态空间方程。

4.4 评价指标

为了分析改进算法的估计性能，对粒子数相同的

PF

算法进行分析和比较。评价指标包括均方根误差Root- Mean- Square Error，RMSE） ERMSE Mean （ 、最大误差（ Error，ME） EME Maximum Relative Error，、最大相对误差（

MRE） EMRE Average Error，AE） EAE：和平均误差（

式中， Ok和 Ok分别为第k个周期的估计值和真实值； est true

L为迭代长度。

4.5 仿真结果与讨论

PF

传统 算法与改进萤火虫算法优化的粒子滤波

Improved Firefly Algorithm- Particle Filter，IFA- PF）

（ 算

A12 SOH 4

法对 系列电池的 估计结果如图 所示，量化性

1

能比较结果如表 所示。4 1

从图 和表 可以看出：经过滤波的估计值与真实

0.05 10%

值之间的误差都在 以内，最大相对误差都在 以

PF SOH

内，具有较高的准确性，说明 算法在电池 估计中具有良好的性能，适用于锂离子电池等非线性、非高斯

HI IFA-PF

系统；结合退化模型和在线 的 算法最大相对

6%

误差可减小到 以内，表明改进算法估计结果具有较

PF

高的精度，优于传统的 算法。

5 结束语

本文提出了一种基于改进粒子滤波算法的锂离子

SOH

电池 在线估计方法。首先从放电电压和放电时间

SOH

方面构建了在线测量的健康指标，建立了其与电池之间的映射关系，并将其应用于系统状态空间方程。通过建立双指数退化模型，采用改进算法实时调整模型参

SOH IFA- PF

数估计电池 。仿真结果表明，基于 算法的

SOH

电池 估计精度高，能够有效地适应具有非线性和非高斯特性的锂离子电池。然而，当电池在动态条件下工作时，情况将更为复杂，测量精度可能降低，因此在未来将对这一具有挑战性的问题进行研究。

参考文献

[1] , , .严旭明 冯光炷 黄雪 锂离子电池负极材料的制备及应[J]. , 2019, 47(7): 22-25.用进展 化工新型材料[2] HE W, WILLIARD N, OSTERMAN M, et al. Prognostics of Lithium- Ion Batteries Based on Dempster- Shafer Theory and the Bayesian Monte Carlo Method[J]. Journal of Power Sources, 2011, 196(23): 10314-10321. [3] HENDRICKS C, WILLIARD N, MATHEW S, et al. A Failure Modes, Mechanisms, and Effects Analysis (FMMEA) of Lithium- Ion Batteries[J]. Journal of Power Sources, 2015, 297: 113-120.