一种锂离子电池单体热失控核心温度估计方法

摘要

本发明公开了一种锂离子电池单体热失控核心温度估计方法。所述方法包括：计算锂离子电池单体热失控过程中的产热功率；根据热失控过程中锂离子电池表面温度及热失控过程中的产热功率建立热失控集总参数热模型；基于建立的模型利用递推最小二乘法在线辨识模型参数，同时结合扩展卡尔曼滤波算法估计热失控过程中的核心温度。该方法采用电池表面温度测量与模型相结合的方式，可实现在锂离子电池发生热失控时，实时估计电池核心温度，为热失控抑制方案优化提供指导。

说明书

技术领域

本发明涉及一种锂离子电池单体热失控核心温度估计方法，属于锂离子电池热安全技术领域。

背景技术

世界范围内环境和能源问题愈发突出。中国采取了很多节能减排举措，其中就包括政策引导新能源汽车的发展。新能源汽车在我国乃至全球都处于高速发展时期，锂离子电池因其高能量密度、长寿命、技术相对成熟等诸多优点，成为新能源汽车最有竞争力的动力源。用户对于新能源汽车续航里程的担心，促使电池科技工作者向高能量密度、大容量动力电池方向努力，而能量密度的提升是以牺牲安全性为代价的，锂离子电池系统安全问题是目前制约新能源汽车普及的重要问题之一。对于目前使用广泛的锂离子动力电池来说，热失控是其安全问题的核心科学问题。

“热失控”是指电池内部出现放热连锁反应引起电池温升速率急剧变化的过热现象。锂离子电池单体的能量密度高一旦发生短路，其内部的化学能将以热能方式迅速释放，电池温度急剧上升最高可至上千摄氏度。

当前针对锂离子电池热失控的研究实验中多采用在电池单体外壳上固定热电偶的形式记录热失控过程中的温度变化，但锂离子电池热失控过程反应剧烈，产热功率较高，极易产生较大的内外温度梯度，而在实际工况中，很难通过测量的手段得到热失控过程中锂离子电池内部的核心温度，因此迫切需要提出一种估计方法用于获得锂离子电池的核心温度，从而为后续热失控抑制方案的研究提供指导。

发明内容

鉴于上述现有技术的不足之处，本发明的目的在于提供一种锂离子电池单体热失控核心温度估计方法，通过将物理建模与算法估计相结合的手段实现热失控过程中电池单体核心温度的准确获取。

本发明提出的一种锂离子电池单体热失控核心温度估计方法，具体包括如下步骤：

步骤一：计算电池单体热失控过程中的产热功率

根据锂离子电池单体热失控过程的电化学及内短路反应，计算随反应时间变化的产热功率；

步骤二：建立热失控集总参数热模型

根据锂离子电池热失控表面温度及热失控过程的产热功率，结合锂离子电池结构特性，建立锂离子电池单体热失控集总参数热模型；

步骤三：基于模型估计热失控过程中的核心温度

根据所述热失控集总参数热模型，利用递推最小二乘法在线辨识模型参数；用状态空间方程对热失控集总参数热模型进行描述，以核心温度T

作为本发明的优选方案，所述步骤一具体包括:

根据所述锂离子电池单体热失控过程的电化学反应计算反应的产热功率Q

Q

根据所述锂离子电池单体热失控过程的内短路过程计算电能释放产热功率Q

根据所述锂离子电池单体热失控过程中的化学反应产热功率Q

Q＝Q

其中，Q

作为本发明的优选方案，所述步骤二具体包括：

假设热失控电池单体的产热功率是从电池核心产生，传递路线为电池核心-电池表面-电池所处环境，且该路线中的产热功率为均匀传递，结合锂离子电池结构特性可建立热失控集总参数热模型：

其中，C

作为本发明的优选方案，所述步骤三具体包括：

对式(4)进行拉普拉斯变换，整理后得锂离子电池热失控集总参数热模型的传递函数式：

将电池热失控的总产热功率Q作为输入，电池内外温差ΔT作为输出，将其离散化可得：

ΔT(k)＝aΔT(k-1)+bΔT(k-2)+cQ(k-1)+dQ(k-2) (6)

其中，ΔT为锂离子电池内外温差，a、b、c、d为系统待辨识参数，k为采样点，t为采样时间，θ(k)为系统待辨识的参数矩阵，h(k)为输入量和输出量的历史集合矩阵，根据递推最小二乘法的递推公式(8)实现模型参数的在线辨识与更新；

其中，

用状态空间方程对热失控集总参数热模型进行描述，以核心温度T

将状态空间方程转化为规范形式：

其中，所述

根据扩展卡尔曼滤波算法递推公式(12)，利用热失控集总参数热模型输出的表面温度与通过实际测量获得的热失控过程电池表面温度的差值不断修正系统状态变量，从而实现热失控过程中锂离子电池核心温度的准确估计

其中，S(k)为实际测量获得的电池表面温度，

与现有技术相比，本发明所具有的有益效果有：

(1)本发明采用了电池表面温度测量与模型相结合的方式估计热失控过程的电池核心温度，克服了现有热失控核心温度测量技术中需要对电池结构进行改造以安装内部温度传感器从而影响电池性能的问题，为实际工况下电池热失控事故核心温度的获取奠定了基础。

(2)本发明采用了递推最小二乘法对热失控集总参数热模型的参数进行在线辨识，克服了实际工况下电池表面与环境之间换热条件实时变化的问题，为实际工况下电池热失控事故核心温度的获取提供了技术支持。

(3)本发明采用了扩展卡尔曼滤波算法对热失控核心温度进行估计，在热失控过程温度发生突变时能够快速跟随并收敛，保证方法的估计精度。

附图说明

图1是本发明提供的锂离子电池单体热失控核心温度估计方法的流程图。

图2是锂离子电池单体热失控过程中总产热功率变化曲线。

图3是本发明提供的圆柱形锂离子电池单体热失控集总参数热模型。

图4是锂离子电池单体热失控过程核心及表面测量温度对比曲线。

图5是锂离子电池单体热失控过程核心估计温度与核心测量温度对比曲线。

图6是锂离子电池单体热失控核心温度估计方法的绝对误差与相对误差。

具体实施方式

下面结合附图并举实施例，对本发明进行详细描述。

本发明涉及一种锂离子电池单体热失控核心温度估计方法，流程详见图1，本实施例以圆柱型锂离子电池热失控过程的核心温度估计为例进行说明：

步骤一：计算圆柱型锂离子电池单体热失控过程中的产热功率，根据所述锂离子电池单体热失控过程的电化学反应计算反应的产热功率Q

Q

Q＝Q

其中，Q

步骤二：建立热失控集总参数热模型，根据锂离子电池热失控表面温度及热失控过程的总产热功率，结合圆柱型锂离子电池结构特性，假设热失控电池单体的产热功率是从电池核心产生，传递路线为电池核心-电池表面-电池所处环境，且该路线中的产热功率为均匀传递，建立锂离子电池单体热失控集总参数热模型，如图3所示，模型表达式为：

其中，C

步骤三：基于模型估计热失控过程中的核心温度，首先对所建热失控集总参数热模型进行在线参数辨识，对式(4)进行拉普拉斯变换，整理后可得锂离子电池热失控集总参数热模型的传递函数式：

将电池热失控的总产热功率Q作为输入，电池内外温差ΔT作为输出，将其离散化可得：

ΔT(k)＝aΔT(k-1)+bΔT(k-2)+cQ(k-1)+dQ(k-2) (6)

其中，ΔT为锂离子电池内外温差，a、b、c、d为系统待辨识参数，k为采样点，t为采样时间选用0.1s，θ(k)为系统待辨识的参数矩阵，h(k)为输入量和输出量的历史集合矩阵，根据递推最小二乘法的递推公式(8)，取得每一次新的ΔT(k)后，就在前一次参数估计结果

其中，

然后，用状态空间方程对热失控集总参数热模型进行描述，以核心温度T

将状态空间方程转化为规范形式：

其中，所述

根据扩展卡尔曼滤波算法，首先设定

其中，S(k)为实际测量获得的电池表面温度，

估计结果：

如图4所示，为了对本发明的方法进行验证，本发明通过对圆柱型锂离子电池进行改造，分别在其几何中心处及表面中心处固定热电偶，引发电池热失控获得锂离子电池单体从自产热直至热失控最高温度过程的电池核心测量温度及表面测量温度，其中核心测量温度仅作为方法验证参考值(本发明方法实际应用中，无需对锂离子电池进行改造，无需布置内部热电偶)，表面测量温度作为方法输入值。

如图5所示，本发明提出的锂离子电池单体热失控核心温度估计方法获得的核心估计温度与实验测得的核心测量温度非常接近，仅在电池发生喷阀时及热失控极速升温时有所偏差但迅速收敛；通过图6可以看出热失控极速升温过程中绝对误差最大为25.6℃，相对误差最大为5.8％，符合热失控瞬间突变的特性，热失控过程最高温度绝对误差仅为8.9℃，相对误差1.1％，误差均在可接受范围内。

综上，本发明提出的锂离子电池单体热失控核心温度估计方法，在不破坏电池结构的前提下，可以很好的反映出锂离子电池单体热失控过程核心温度变化情况，同时对热失控过程最高温度的准确估计可为热失控抑制方案优化提供指导。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，依然可以有各种修改。凡在本发明的的精神和原则之内所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。