一种基于健康状态的锂离子电池荷电状态估计方法

1.本发明涉及锂电池技术领域，特别涉及一种基于健康状态等效电路模型的锂离子电池荷电状态估计方法。


背景技术：

2.随着能源危机和环境危机的日趋严峻，世界各国越来越重视新能源的发展。其中，锂离子电池作为环保型储能器件，具有功率密度高、能量密度高和自放电率低等优势，已被广泛应用于新能源汽车汽车和智能电网等领域。在电池管理系统中，电池荷电状态的估计是至关重要的，对于保障系统安全有着重要作用。然而，模型参数受噪声干扰，电池寿命受内外部因素影响以及电池充放电的循环使用，使得电池荷电状态的估计仍存在许多挑战。专利【一种基于噪声自适应粒子滤波的电池荷电状态估计方法和装置】采用基于噪声的自适应粒子滤波实现电池的荷电状态估计，但是其没有考虑电池的老化状态即容量的衰退。专利【一种基于自适应模型的锂电池荷电状态估计装置及方法】虽然考虑了电池的老化状态，但是其采用拓展卡尔曼滤波实现电池的荷电状态估计，拓展卡尔曼滤波噪声为高斯噪声，由于锂离子电池具有较强的非线性关系且系统噪声为非高斯噪，其处理效果不是很理想。专利【锂离子电池荷电状态和健康状态联合估计方法】采用双滑膜观测器实现电池荷电状态和健康状态的联合估计，其考虑了容量的衰减，但是荷电状态和健康状态变化的时间尺度不同，健康状态的变化比较缓慢，荷电状态处于实时的变化过程，不适应用于实际的成产。


技术实现要素：

3.本发明的目的是提供一种基于健康状态的锂离子电池荷电状态估计的装置及方法，能够根据容量衰减程度不同、环境温度不同和充放电倍率不同，提高锂离子电池荷电状态估计的精度。
4.本发明所提供的技术方案是：
5.一种基于健康状态的锂离子电池荷电状态估计方法，包括：
6.步骤s1，对锂电池进行循环充放电实验，提取容量、恒流充电时间、恒压充电时间、放电压差与放电时间的比值作为锂电池老化的有效特征，建立电池剩余使用寿命预测模型，将提取的老化特征组成特征向量作为预测模型的输入；
7.步骤s2，建立锂离子电池的等效电路模型，采取小电流充放电的方式，对电池的开路电压-荷电状态曲线进行测定，再对电池等效电路模型进行参数辨识；
8.步骤s3，在25℃室温下，对锂电池进行循环充放电实验，并对充放电电流进行积分，获得每个充放电循环下的电池实际容量，并绘制实际可用容量与循环次数的关系曲线；将步骤s1中预测得到电池的剩余使用寿命替换掉循环次数，通过拟合得到电池实际容量与剩余使用寿命的函数关系式；
9.步骤s4，在不同的温度下对电池进行彻底的恒流放电；建立了电池工作温度对电
池容量的补偿方程；
10.步骤s5，不同倍率下对电池充放电实验，建立电池不同充放电倍率下对电池容量的补偿方程；
11.步骤s6，根据步骤s3、s4和s5获得电池最终的实际可用容量，然后根据步骤s2中等效电路模型建立的状态空间方程，采用粒子滤波的方法实现电池荷电状态的估计。
12.一种基于健康状态的锂离子电池荷电状态估计装置，包括电源供电模块，ac-dc转换模块、dc-dc转换模块、控制模块，模拟量采集模块和大功率电阻；所述电源供电模块，用于提供各模块所需的电源；所述ac-dc模块，用于将常用220v交流电转换为直流电；所述dc-dc模块：实现锂电池的放电操作，将电池中的电量向功率电阻输出所述控制模块，用于控制锂离子电池的充放电电流；所述数据采集模块，用于采集锂离子电池端电压、充电电流以及工作温度；所述大功率电阻，用作负载。
13.本发明的技术效果在于：
14.(1)本发明通过锂离子电池剩余使用寿命的预测映射为电池的实际基准容量，考虑了电池的健康状态(即老化状态)，同时将工作温度和充放电倍率作为影响因子，实现了电池实际可用容量的补偿。
15.(2)本发明采用二阶等效电路模型，同时采用粒子滤波的方法降低了过程噪声和测量噪声的干扰，提高了锂离子电池的荷电状态估计的准确性。
附图说明
16.图1为本发明中锂离子电池荷电状态估计装置的框图。
17.图2为本发明中锂离子电池荷电状态估计方法的流程图。
18.图3为本发明中锂离子电池二阶等效电路模型示意图。
19.图4为本发明中剩余使用寿命预测框架图。
20.图5为本发明中粒子滤波算法的流程图。
具体实施方式
21.下面参照附图对本发明进行更全面的描述，在附图中，相同的标识表示相同或相似的组件或元素。
22.图1表示本发明提供的荷电状态估计的装置示意图，包括电源供电模块，ac-dc转换模块，dc-dc转换模块，控制模块，模拟量采集模块，大功率电阻，betaq电池循环老化测试仪，恒温箱，温度、电压和电流传感器。
23.其中，模拟量采集模块与ac-dc转换模块、锂电池组、dc-dc转换模块和控制板连接，控制板与ac-dc转换模块和dc-dc转换模块连接，ac-dc转换模块与锂电池组、模拟量采集模块和控制版连接，dc-dc转换模块与锂电池组、模拟量采集模块、控制板和大功率电阻连接；
24.电源供电模块用于将外部输入电压转换为控制模块，模拟量采集模块以及数据采集模块所需的电源电压。
25.ac-dc模块用于将常用220v交流电转换为直流电，实现对电池恒流恒压的充电方式。
26.dc-dc模块实现锂电池的放电操作，将电池中的电量向功率电阻输出。
27.控制模块用于向ac-dc转换模块和dc-dc转换模块输出pwm脉冲信号，控制锂离子电池的充放电电流；同时还用于荷电状态估计中方法的处理和计算。
28.所述大功率电阻模块作为负载，实现锂电池能量的消耗。
29.各种传感器模块用于锂离子电池端电压，充电电流以及工作温度等数据的采集。
30.恒温箱用于控制锂离子电池充放电的工作温度。
31.betaq电池循环老化测试仪用于锂离子电池循环充放电实验，并提供数据信号供各传感器采集。
32.如图2所示，本发明提供的一种基于健康状态的锂离子电池荷电状态估计的方法包括如下步骤：
33.步骤s1，利用betaq电池循环老化测试仪对电池进行循环充放电，提取每次循环充放电的老化特征容量、恒流充电时间、恒压充电时间、放电压差与放电时间的比值，然后采用基于邻域分析的方法进行特征增强，建立电池剩余使用寿命预测模型。
34.图4为本发明的剩余使用寿命预测框架图。锂离子电池得剩余使用寿命预测主要分为两部。第一步离线训练，首先需要电池循环充放电采集数据，然后对数据进行特征得提取和筛选，再采用邻域分析法对特征增强，最后作为训练数据集采用相关向量机的方法对数据进行训练，得到锂离子电池剩余使用寿命预测模型。第二部分为在线运行，对数据的处理与离线训练类似，最后将处理得到的测试数据作为预测模型的输入，最终得到电池的剩余使用寿命。
35.本发明优选8个电池进行实验，先通过一段时间的循环充放电实验获得训练数据集，将训练数据集作为输入，剩余使用寿命作为输出，经过训练后，得到了电池剩余使用寿命的预测模型。然后，可以采用该模型进行在线的剩余使用寿命预测。
36.步骤s2，首先需要对电池的开路电压-荷电状态曲线进行测定，然后采用带遗忘因子的递推最小二乘实现电池模型的参数辨识。
37.先让电池以恒流恒压的方式充电，直至充电电流下降至电池的下限截止电流。将电池充满电结束时的荷电状态定义为100％；
38.将锂离子电池静置5小时，实现电池的退极化。然后以0.25c的放电倍率进行放电，放电次数为20次，直至电池的电压下降至其下限截止电压3v。将放电结束时的电池荷电状态定义为0％。
39.在每次脉冲放电结束时，将电池置于开路状态5个小时，同样实现电池的退极化，然后对电池的终端电压进行测量并定义为放电的开路电压。另外，将电池充电池荷电状态为100％的状态，获得电池充电的开路电压。
40.对获得的放电开路电压和充电的开路电压取平均值，将平均值作为实际的开路电压。
41.获得的开路电压-荷电状态的曲线是非线性的，若直接对该曲线进行线性化，将会增加荷电状态估计的难度。本技术中将荷电状态-开路电压曲线分解为一组线性段共9段并进行简化，每个线段的形式如下所示：
42.ocv＝c0 c1soc
43.其中ocv为开路电压，soc为荷电状态，c0和c1为常数系数。测定的开路电压-荷电状
态曲线如图2所示，以查表的方式从该曲线中得到电池开路电压信息。
44.图3为本发明的锂离子电池二阶等效电路模型示意图，主要由等效内阻rs，极化电容c
p1
和极化电阻r
p1
，浓差极化电容c
p2
和浓差极化电阻r
p2
，开路电压v
oc
、端电压v
t
和充电电流i组成。
45.根据等效电路模型，可得其等效电路表达式为：
[0046][0047]
通过上述获得的开路电压-荷电状态函数关系以及建立的等效电路模型，接下来采用带遗忘因子的递推最小二乘实现等效电路模型的参数辨识，带遗忘因子的递推最小二乘法公式如下所示：
[0048][0049]
其中，θ表示待辨识参数，k表示增益矩阵，p表示协方差矩阵，为测量值，λ表示遗忘因子。将电池等效电路模型公式进行一系列的频域以及差分形式的转换为带遗忘因子的递推最小二乘法公式形式：
[0050][0051]
具体步骤如下：
[0052]
设置初始化的参数，设定荷电状态soc(0),开路电压v
oc
(0)和时间间隔δt的初始值；
[0053]
实时的计算荷电状态soc(k)和开路电压v
oc
(k)；
[0054]
设置协方差矩阵p0、待辨识参数和遗忘因子λ的初始值；
[0055]
通过电池等效电路模型转换的带遗忘因子的递推最小二乘法公式求出k1、k2、k3、k4和k5。
[0056]
求出等效电路模型参数rs、r
p1
、r
p2
、c
p1
和c
p2
。
[0057]
步骤s3，在25℃室温下，根据实验平台对电池进行循环充放电实验，通过高精度电流传感器采集电流，并对电流进行积分，获得每个充放电循环下的电池实际容量，并绘制实际可用容量与循环次数的关系曲线。然后，将步骤s1中预测得到电池的剩余使用寿命替换掉循环次数，通过拟合获得电池实际容量与剩余使用寿命的函数关系式。
[0058]cn
(x)＝ax3 bx2 cx d
[0059]
其中，cn表示电池的实际容量，a、b、c和d为常数系数，x表示电池的剩余使用寿命。
[0060]
步骤s4，在不同的温度下对电池进行彻底的恒流放电。首先在25℃的环境下，采用
1c的放电倍率放电，同样通过电流传感器获得放电电流，并对电流进行积分，得到此时25℃工作环境下的电池实际可用容量。然后，再不同的温度环境下，获得不同的实际可用容量，并与25℃环境下的实际可用容量作比值。最后建立了电池工作温度对电池容量的补偿方程。
[0061]
ξ(x)＝a1x3 b1x2 c1x d
[0062]
其中，ξ表示温度对电池实际容量的补偿，a、b、c和d为常数系数，x表示电池的工作温度。
[0063]
步骤s5，不同倍率下对电池充放电实验。与步骤s4类似，首先获得1c放电倍率下的电池实际可用容量，然后获得不得放电倍率下的实际可用容量，因此得到关于容量的比值。最后建立电池不同充放电倍率下对电池容量的补偿方程。
[0064]
η(x)＝a2x3 b2x2 c2x d2[0065]
其中，η表示充放电倍率对电池实际容量的补偿，a、b、c和d为常数系数，x表示电池的充放电倍率。
[0066]
步骤s6，将步骤s3中获得电池的老化状态与实际可用容量间的定量关系得到电池的实际可用容量，并以此为基准，将步骤s4和s5中获得的电池工作温度和充放电倍率与实际可用容量的定量补偿系数对实际可用容量进一步补偿和校正，得到最终的电池实际可用容量。
[0067]
qr＝cn(rul)
·
ξ(t)
·
η(c)
[0068]
锂离子电池荷电状态可用如下公式所示：
[0069][0070]
其中soc(t0)为初始的荷电状态，cn为锂离子电池的实际可用容量，i(t)为电池的充放电电流。将该公式转换为离散形式为：
[0071][0072]
其中t为离散周期，k为采样点。
[0073]
然后根据等效电路模型，等效电路模型表达式为：
[0074][0075]
其中，u1和u2分别为极化电阻r
p1
和r
p2
上的电压。让τ成为rc电路的时间常数，即τ＝rc。将上述公式结合，通过离散化得到了锂离子电池的状态方程和观测方程如下：
[0076][0077]vk
＝v
oc(soc,k)-u
p1,k-u
p2,k-reik vk[0078]
估计中的状态变量为：
[0079]
xk＝[u
p1,k
,u
p2,k
,sock]
[0080]
假设荷电状态估计过程中的过程噪声和测量噪声服从高斯分布，可得：
[0081][0082]
g(xk,uk,σk)＝v
oc(soc,k)-u
p1,k-u
p2,k-reik[0083]
根据步骤s3中将筛选并增强得到的老化特征值输入构建的锂离子电池剩余使用寿命模型，得到电池的剩余使用寿命rul，并将该值带入容量基准公式中，获得电池的基准容量为cn,然后将温度和充放电倍率对电池容量的补偿系数分别带入中，通过计算获得最终的电池实际可用容量。最后采用粒子滤波的方法实现电池荷电状态的估计。
[0084]
图5为本发明所述的粒子滤波算法的流程图，实现电池荷电状态的在线估计。
[0085]
(1)首先，初始化，k＝0，设置噪声方差值、有效阈值；初始的极化电压值和荷电状态值；粒子的数量等。初始化粒子，从初始密度分布函数中随机抽取n个状态粒子：x
0i
(i＝1,2,...n)；初始化各个粒子的权重：w
0i
＝1/n(i＝1,2,...n)
[0086]
(2)然后，更新状态的时间：k＝k-1；并利用状态方程公式(11)预测k时刻的值
[0087][0088]
(3)其次，利用观测方程(12)计算每个粒子的权值，计算公式为：
[0089][0090]
(4)再次，对k时刻的状态进行估计：
[0091][0092]
(5)判断是否需要进行重采样，若需要重采样，需将粒子的权值重新赋值为1/n。
[0093]
(6)判断循环是否结束，结束则实现了荷电状态的估计，否则返回步骤(2)。