一种具有温度补偿功能的锂电池电荷状态估算方法与流程

本发明属于锂电池充放电技术领域，具体涉及一种具有温度补偿功能的锂电池电荷状态估算方法。

背景技术：

随着煤、石油等化石能源的逐渐枯竭以及人们对环保问题的不断重视，以动力电池作为主要能源供给的电动汽车以零污染、节能效率高的优点受到人们越来越高的重视。单体锂电池产品额定电压为3.2v和4.2v两种，如若想满足动力电池所需的电压及容量，需要将大量的单体电池进行串并联成组使用。动力电池作为电动汽车最核心的零部件之一，需要其运行状态进行实时精准的检测。电荷状态soc(stateofcharge)作为锂电池最关键的技术参数之一，在实际情况下是无法进行直接测量的，因此，人们对锂电池soc的估算方法进行了深入的研究。

目前，主流的锂电池soc估算方法有安培积分法、开路电压法、卡尔曼滤波算法及扩展卡尔曼滤波算法等。安培分法对锂电池的soc初始值依赖性较大，若soc初始值有误差，则会导致锂电池soc的估算不准确。开路电压法较为简单，只需要将锂电池充分静置后查表便可以得到锂电池soc，但静置时间一般在2小时以上，所以该方法所消耗的时间较长，不易推广使用。扩展卡尔曼滤波算法对锂电池soc初始值依赖性不大，但在计算时结果容易发散导致估算不准确。

并且，无论哪一种soc估算算法，温度对锂电池soc的估算精度都有着很大的影响，所以在对锂电池soc进行估算时，有必要引入温度补偿功能来提高其估算精确性。

技术实现要素：

为了解决目前锂电池电荷状态估算算法对初始值的依赖性、估算结果容易发散并且未考虑温度对锂电池电荷状态估算的精度的影响等相关问题，提出了一种具有温度补偿功能的锂电池电荷状态估算方法。

为实现上述技术目的，所采用的技术方案是：一种具有温度补偿功能的锂电池电荷状态估算方法，其特征在于，包括以下步骤：

步骤1、建立各单体锂电池二阶rc网络等效电路模型；

步骤2、建立单体锂电池温度模型；

dee＝m*cp*dtr(2)

qloss＝qconv+qcond(4)

qconv＝hconv1sarea(tr-tair)+hconv2sarea(ty-tz)(5)

其中，k表示时间步长；ee:电池内部能量；qgen(k)：电池内部产热速率；m：电池质量，m>0；cp：电池比热容，取值为130-880j/kg/k；r0(soc,tb)：电池内阻；i：工作电流；与分别表示时间步长k处电池内部电化学极化电压与浓差极化电压；r1(soct,tb)k：表示时间步长k处电化学极化内阻；r2(soct,tb)k：表示时间步长k处浓差极化内阻；dtcell：单体电池随时间的温度变化；hconv1：空气与锂电池之间对流换热系数，其取值为5-10w/(m²*k)；hconv2：锂电池与锂电池之间对流换热系数，取值为5-10w/(m²*k)；sarea：热交换面积，取值大于0；kt：材料导热系数，取值为100-300w/(m*k)；a：垂直于热流方向的面积，取值大于0；d:层间距离，取值大于0；tr：第r节单体电池温度；ty：单体电池y的温度；tz：单体电池z的温度；联立以上公式求解，即可得出各单体电池温度；

步骤3、建立单体锂电池二阶rc网络等效电路的离散状态空间模型；

步骤4、以步骤一的各单体锂电池二阶rc网络等效电路模型、步骤二的单体锂电池温度模型以及步骤三的单体锂电池二阶rc网络等效电路的离散状态空间模型为基础，利用无迹卡尔曼滤波算法估算锂电池的电荷状态。

单体锂电池二阶rc网络等效电路的离散状态空间模型为：

其中，ts为采样时间，wk为过程噪声，sock表示时间步长k处锂电池的荷电状态，cq表示锂电池额定容量，c1(sock,tb)k，c2(sock,tb)k分别表示时间步长k处的锂电池电化学极化电容和浓差极化电容，表示时间步长k处测量输出电压，em(sock,tb)k表示时间步长k处开路电压，vk为测量噪声，i表示工作电流。

利用无迹卡尔曼滤波算法估算锂电池的电荷状态的具体方法是：

步骤4.1、使用状态初始值x[0]和状态估计误差协方差p初始化滤波器

是状态估计，表示使用在时间步长0,1,2,…,kb的测量值对在时间步长ka的状态估计，这里其中，soc0∈[0,1],分别表示锂电池电荷状态、浓差极化电压电化学极化电压的初始估计协方差，

步骤4.2、对每一个时间步长k，使用测量数据y[k]更新状态估计和状态估计误差协方差：

4.2.a、在时间步长k处选择ε点

其中，c＝α²(m+ζ)，α取值为[0,1]，ζ取值为0，m取值为3；

4.2.b、根据方程式(8)计算每个ε点的预测测量值

其中，um[k]表示时间步长k处方程式(8)的输入；

4.2.c、将每个ε点的预测测量值结合起来得到时间步长k处的预测测量值

4.2.d、估计步骤4.2.c得到的时间步长k处的预测测量值协方差

其中，r[k]为时间步长k处的测量噪声协方差矩阵，取值为[0,1]，β取值为2；

4.2.e、估计与之间的互协方差

4.2.f、得到在时间步长k处所估计的状态变量值及状态估计误差协方差

其中，为卡尔曼增益矩阵，

步骤4.3、预测下一时间步长的状态变量值及状态估计误差协方差

4.3.a、选择时间步长k处的ε点

4.3.b、根据方程式(7)计算每一个ε点的所预测的状态变量值

其中，us[k]表示表示时间步长k处方程式(7)的输入；

4.3.c、结合每一个ε点的所预测的状态变量值得到步长k+1处的预测状态量值

4.3.d、计算步骤4.3.c得到的步长k+1处的预测状态量值协方差

其中，为过程噪声协方差矩阵，(max(|dsock|))²∈[0，1]，

本发明的有益效果是：

1、在锂电池二阶rc网络等效电路模型的基础上，考虑了温度对锂电池内部参数(开路电压、内阻、极化电压)等的影响并建立了相应的模型，提高了锂电池模型的精度。

2、利用无极卡尔曼滤波算法并结合考虑温度影响的二阶rc网络等效模型实现对锂电池电荷状态的估计，电荷状态估算精度更高。

3、该方法容易实现，将相应算法编写为c语言并固化到单片机中即可实现，该算法对锂电池组中给单体电池电荷状态的估计可以无限制扩展。

附图说明

图1为锂电池电荷状态估算方法整体框图；

图2为锂电池二阶rc网络等效电路模型图；

图3为锂电池动态脉冲放电曲线图；

图4为锂电池二阶rc网络等效电路模型参数辨识曲线及相对误差图；

图5为锂电池二阶rc网络等效电路模型参数辨识曲线及相对误差局部放大图；

图6为锂电池电荷状态估算方法simulink模型图；

图7为锂电池组各单体电池温度模型图；

图8为无迹卡尔曼滤波估计锂电池电荷状态方法图；

图9为各单体电池电荷状态实测与估计曲线图；

图10为各单体电池两种电荷状态估算方法与测量值的相对误差曲线图；

图11为各单体电池不同温度模型电荷状态实测与估计曲线图；

图12为各单体电池不同温度模型电荷状态估算方法与测量值的相对误差曲线。

具体实施方式

下面结合具体实施例对本发明作进一步说明，以使本领域的技术人员可以更好的理解本发明并能予以实施，但所举实施例不作为对本发明的限定。

一种具有温度补偿功能的锂电池电荷状态估算方法，包括以下步骤：

步骤1、建立各单体锂电池二阶rc网络等效电路模型。

步骤1.1、在5℃、20℃、40℃恒温条件下对实验用锂电池进行动态脉冲放电实验，其中，20℃条件下动态脉冲放电曲线如图3所示。

步骤1.2、根据实验测得的脉冲放电数据，利用最小二乘法对锂电池二阶rc网络等效电路模型进行参数辨识，辨识结果及相对误差如图4所示，平均相对误差3.3196mv，局部放大图如图5所示。

步骤2、结合锂电池相关物理属性(质量、体积、比热容)，根据热力学第一定律、傅里叶定律及牛顿冷却定律建立单体电池温度模型。

dee＝m*cp*dtr(2)

qloss＝qconv+qcond(4)

qconv＝hconv1sarea(tr-tair)+hconv2sarea(ty-tz)(5)

其中，k表示时间步长，ee:电池内部能量；qgen(k)：电池内部产热速率；m：电池质量，取值1kg；cp：电池比热容，取值为810.5j/kg/k；i：工作电流；r0(soc,tb)k：表示时间步长k处电池内阻；与分别表示时间步长k处电池内部电化学极化电压与浓差极化电压；r1(soct,tb)k：表示时间步长k处电化学极化内阻；r2(soct,tb)k：表示时间步长k处浓差极化内阻；dtcell：单体电池随时间的温度变化；hconv1：空气与锂电池对流换热系数，其取值为10w/(m²*k)；hconv2：锂电池与锂电池对流换热系数，取值为5w/(m²*k)；sarea：热交换面积，取值为0.102m²；kt：材料导热系数，取值为200w/(m*k)；a：垂直于热流方向的面积，取值为1e-3m²；d:层间距离(材料厚度)，取值为0.1m；tr：第r节单体电池温度；ty：单体电池y的温度；tz：单体电池z的温度；联立以上公式求解，即可得出各单体电池温度。

步骤3、建立单体锂电池二阶rc网络等效电路的离散状态空间模型

其中，ts为采样时间，wk为过程噪声，sock表示时间步长k处锂电池的荷电状态，cq表示锂电池额定容量，c1(sock,tb)k，c2(sock,tb)k分别表示时间步长k处的锂电池电化学极化电容和浓差极化电容，表示时间步长k处测量输出电压，em(sock,tb)k表示时间步长k处开路电压，vk为测量噪声，i表示工作电流。

步骤4、利用无迹卡尔曼滤波算法(简称ukf)估算锂电池soc，具体步骤为：步骤4.1、使用状态初始值x[0]和状态估计误差协方差p初始化滤波器：

是状态估计，表示使用在时间步长0,1,2,…,kb的测量值对在时间步长ka的状态估计，这里

步骤4.2、对每一个时间步长k，使用测量数据y[k]更新状态估计和状态估计误差协方差：

4.2.a、在时间步长k处选择ε点

其中，c＝α²(m+ζ)，取决于状态m的数量及参数α、ζ，这里α取值为1，ζ取值为0，m取值为3。

4.2.b、根据方程式(8)计算每个ε点的预测测量值

这里，um[k]表示时间步长k处方程式(8)的输入，um[k]的输入包括单体电池的温度tr和工作电流i。

4.2.c、将每个ε点的预测的测量值结合起来得到时间步长k处的预测测量值

4.2.d、估计步骤4.2.c得到的时间步长k处的预测测量值协方差

这里，r[k]为时间步长k处的测量噪声协方差矩阵，取值为1e-3，β取值为2。

4.2.e、估计与之间的互协方差

4.2.f、得到在时间步长k处所估计的状态变量值及状态估计误差协方差

其中，为卡尔曼增益矩阵。

步骤4.3、预测下一时间步长的状态变量值及状态估计误差协方差

4.3.a、选择时间步长k处的ε点

4.3.b、根据方程式(7)计算每一个ε点的所预测的状态变量值

这里，us[k]表示表示时间步长k处方程式(7)的输入，us[k]的输入包括单体电池的温度tr，工作电流i和额定容量cq。

4.3.c、结合每一个ε点的所预测的状态变量值得到步长k+1处的预测状态量值

4.3.d、计算步骤4.3.c得到的步长k+1处的预测状态量值协方差

这里，

步骤5、在matlab/simulink中对上述模型及算法进行建模及仿真分析建立整体锂电池soc估算算法simulink模型如图6所示，整体模型由锂电池组模型、锂电池组温度模型及无迹卡尔曼滤波估算soc算法三部分组成。分别如图7、图8所示。

锂电池组温度模型经过计算锂电池与环境之间对流换热、锂电池与锂电池之间对流换热、锂电池与锂电池之间传导传热，得到锂电池的温度tr(其中，r＝1,2,…,8,表示第r节锂电池)，得到单体电池r的温度之后，结合此时锂电池的输出电压，利用无迹卡尔曼滤波算法估计出锂电池此时所对应的soc。

为了验证该算法的精度，选取八节额定容量30ah额定电压3.7v锂电池串联成组作为实验对象，采用上述温度模型的ukf估算锂电池soc和未采用上述温度模型的ukf估算锂电池soc仿真及实验结果如图9和图10所示。

从图9可以看出，八节串联锂电池在考虑温度影响后的ukf对soc的估计与未考虑温度影响的ukf对soc的估计相比，与实测值更为接近。

从图10可以看出，考虑温度影响的ukf对各单体电池的soc估计比未考虑温度影响的ukf对各单体电池的soc估计相对误差更低，相对误差在1％左右。未考虑温度影响的ukf对各单体电池的soc估计相对误差在2％左右，相对精度提高了50％。同时，为了更进一步验证该算法的精度，将上述温度模型与其他温度模型做了比较，阅览相关文献，也有构建类似温度补偿模型，通过研究，了解到类似温度补偿模型是在安培积分法的基础上加入温度修正因子以表示温度对锂电池soc估计的影响，这种方法可以是考虑了锂电池内部电阻产热，但是未考虑锂电池与锂电池之间的对流传热与传导传热。因此将采用本文所构建的温度模型与只考虑锂电池内部电阻产热所构建的模型做了对比，对比结果如图11和图12所示。

从图11可以看出，采用本文所构建的温度模型对各单体电池soc的预测值与只考虑锂电池内部电阻产热所构建的模型对各单体电池soc的预测值相比，与真实值更为接近。

从图12可以看出，采用本文所构建的温度模型对各单体电池soc的预测值与只考虑锂电池内部电阻产热所构建的模型对各单体电池soc的预测值相比，与真实值的相对误差更低，相对误差在1％左右，只考虑锂电池内部电阻产热所构建的模型对各单体电池soc的预测相对误差在1.5％左右，相对精度提高了33％。

因此所提出的考虑温度补偿的ukf估算锂电池soc方法与未考虑温度补偿的ukf估算锂电池soc方法和只考虑锂电池内部电阻产热所构建的温度补偿的ukf估算锂电池soc方法相比，相对精度分别提高了50％和33％。该方法更有助于提高锂电池soc的估计精度。