一种锂离子电池荷电状态估计方法技术

一种锂离子电池荷电状态估计方法，本发明专利技术涉及锂离子电池荷电状态估计方法。本发明专利技术是为了解决现有开路电压法估计电池荷电状态耗时长及不适于在线实时应用的问题。具体是按照以下步骤进行的：步骤一：开路电压-荷电状态数据库的获取；步骤二：开路电压拟合公式的获取；步骤三：获取电池终止工作后端电压随时间变化的U(t)数据；步骤四：根据步骤二和步骤三得到公式(7)中的γ值；步骤五：根据步骤四获得的γ值和步骤一获得的开路电压-荷电状态数据库，得到电池的荷电状态。本发明专利技术应用于锂离子电池荷电状态估计领域。

【技术实现步骤摘要】

本专利技术涉及锂离子电池荷电状态估计方法。
技术介绍
锂离子电池荷电状态(StateofCharge，SOC)经常被用来表征电池的剩余可用电量，定义为电池剩余电量与电池额定容量的百分比。准确及时地估计电池荷电状态对电池的安全、高效利用具有重大意义。电池荷电状态不能直接测量，必须经由与荷电状态相关的电池外特性参数(如电压、电流、内阻等)间接获取。现有的锂离子电池荷电状态估计方法中，完全放电实验法虽然估计精度高，但不可作为电池实时在线估计方法，只能作为实验室估计方法；计时电流法(安时法)操作简单，易于实现，但初始荷电状态值选取不当以及电流测量的不精确使得电池荷电状态估计的误差较大，且这一误差会逐渐累积，导致估计误差过大；阻抗法需要电池长时间静置达到稳定状态，且测试仪器成本高，荷电状态与阻抗的非线性关系选择不当会引入较大误差，无法满足在线应用；最新出现的模型法与估计算法(如卡尔曼滤波算法、粒子滤波算法等)结合的方法中，多数无法考虑电池在充放电过程中内部复杂的电化学反应机理，采取线性化近似处理电池充放电的非线性特性，引入了较大误差，且大量的计算也对微处理器提出了较高的要求，增加了估计成本；开路电压法的原理是根据开路电压与荷电状态的关系曲线，测定电池的开路电压，再根据查表法等读取电池荷电状态值，其估计精度较高，但传统的开路电压法需要将电池静置数小时达到稳定状态，耗时过长，不适于实时应用。
技术实现思路
本专利技术是为了解决现有开路电压法估计电池荷电状态耗时长及不适于在线实时应用的问题，而提出的一种锂离子电池荷电状态估计方法。一种锂离子电池荷电状态估计方法按以下步骤实现：步骤一：开路电压-荷电状态数据库的获取；步骤二：开路电压拟合公式的获取：U(t)＝γ-a·tb-c·tdln(t)(7)其中所述U(t)为电池非平衡态的端电压，γ为电池的稳态开路电压，a,b,c,d为拟合参数，t为时间；步骤三：获取电池终止工作后端电压随时间变化的数据；步骤四：根据步骤二和步骤三得到公式(7)中的γ值；步骤五：根据步骤四获得的γ值和步骤一获得的开路电压-荷电状态数据库，得到电池的荷电状态。专利技术效果：本专利技术从电池工作的电化学基本原理出发，获得的开路电压估计方法具有普适通用、快速准确、可实时应用的特性，易于植入嵌入式电池管理系统估计荷电状态，只需几分钟之内的实验数据即可拟合获得精度达到99％荷电状态估计值，耗时长度至少缩小了约两个数量级，大大改善了现有开路电压法估计荷电状态需要耗时数小时的缺点，实用性较强，能够满足现有设备的需求。且开路电压的估计公式对由任何电极材料组成的电池，在0-100％的荷电状态范围内、在保证电池安全操作的温度范围及放电倍率条件下均适用。附图说明图1为锂离子电池开路电压-荷电状态的关系曲线图；图2为锂离子电池拓扑结构示意图；Li为金属锂，M为过渡金属，O为氧，C6为石墨，e-为电子，下标x代表氧原子数；图3为电流中断后电池端电压随时间的变化曲线图；图4为拟合300s范围内的实验数据获得的拟合结果与实验结果对比图；图5为拟合600s范围内的实验数据获得的拟合结果与实验结果对比图。具体实施方式具体实施方式一：一种锂离子电池荷电状态估计方法包括以下步骤：步骤一：开路电压-荷电状态数据库的获取；步骤二：开路电压拟合公式的获取：U(t)＝γ-a·tb-c·tdln(t)(7)其中所述U(t)为电池非平衡态的端电压，γ为电池的稳态开路电压，a,b,c,d为拟合参数，t为时间；步骤三：获取电池终止工作后端电压随时间变化的U(t)数据；步骤四：根据步骤二和步骤三得到公式(7)中的γ值；步骤五：根据步骤四获得的γ值和步骤一获得的开路电压-荷电状态数据库，得到电池的荷电状态。即根据获得的γ值，查表开路电压-荷电状态数据库，即可获得电池荷电状态估计值。具体实施方式二：本实施方式与具体实施方式一不同的是：所述步骤一中开路电压-荷电状态数据库的获取的具体过程为：使电池在预设温度条件下静置1h-3h后，用C/100-C/10的倍率全充放循环两次，记录电池端电压和充放电容量值，对容量值进行归一化处理，获得电压-荷电状态数据，由于倍率非常小，所以充放电过程中的极化可以认为近似相等，由公式(1)获得电池开路电压公式(2)，即将充电电压和放电电压取平均，获得的电压值作为电池平衡态的开路电压值，最终获得开路电压-荷电状态数据库如图1所示；其中所述预设温度为-30℃-60℃；其中所述Vc为充电电压，Vdis为放电电压，V0为稳态电压，I为电流，Ri为阻抗。具体实施方式三：本实施方式与具体实施方式一或二不同的是：所述步骤二中稳态开路电压拟合公式的获取过程为：(1)如图2所示，设计锂离子电池，所述电池由正极、负极、隔膜和电解液组成；(2)采用有限元法将锂离子电池分割为s段，s取值为∞；视每一段内电池的各个物理化学性质相同，中断外置电路时，电池终止工作，电池内部没有净电流，电池端电压的缓慢恢复过程主要受温度、带电粒子浓度及扩散等因素影响；(3)在相邻两段之间设计电极反应，电极反应电势符合Nernst方程，即公式(3)；由于电池内部各个界面上的带电粒子浓度在不断变化，界面之间的离子浓度梯度也变随时间不断变化，所以认为活度αi(t)是时间的函数；其中αi为粒子活度，αi(t)为时间的函数，Ei(t)为电极反应电势，Ei'为形式电势，R为气体常数，T为温度，F法拉第常数，n为电极反应电子转移数目，ci(t)为第i段粒子浓度；电池的端电压等于电池两端电极反应电势的差值；(4)根据对数函数的性质，电池终止放电后，电池端电压的弛豫恢复过程符合公式(5)；U(t)＝γ-δ(t)-β(t)lnλ(t)(5)δ(t)、β(t)和λ(t)是时间t的函数(δ(t)、β(t)和λ(t)为中间变量)；(5)设置公式(5)中δ(t)、β(t)和λ(t)得到公式(6)；将公式(6)代入公式(5)获得电池终止放电后开路电压拟合公式(7)U(t)＝γ-a·tb-c·tdln(t)(7)。具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：所述步骤三中获取电池终止工作后随时间变化的U(t)数据的具体过程为：电池终止工作后，外电路电流为0，电压的弛豫过程表现为电压逐渐变大，直至达到稳态，采用电池测试设备获取电池端电压随时间的变化数据。具体实施方式五：本实施方式与具体实施方式一至四之一不同的是：所述步骤四中得到公式(7)中的γ值的具体过程为：采用最小二乘法辨识方法，当公式(8)中的E达到最小值时拟合获得公式(7)中的γ值，即为电池的稳态开路电压值；拟合过程为将步骤三获取的U(t)数据，输入数据处理器对数据进行拟合，若拟合相关系数R≥99％则停止拟合，获得电池稳态时的开路电压值，否则本文档来自技高网...

【技术保护点】
一种锂离子电池荷电状态估计方法，其特征在于，所述状态估计方法包括以下步骤：步骤一：开路电压‑荷电状态数据库的获取；步骤二：开路电压拟合公式的获取：U(t)＝γ‑a·tb‑c·td ln(t)        (7)其中所述U(t)为电池非平衡态的端电压，γ为电池的稳态开路电压，a,b,c,d为拟合参数，t为时间；步骤三：获取电池终止工作后端电压随时间变化的U(t)数据；步骤四：根据步骤二和步骤三得到公式(7)中的γ值；步骤五：根据步骤四获得的γ值和步骤一获得的开路电压‑荷电状态数据库，得到电池的荷电状态。

【技术特征摘要】
1.一种锂离子电池荷电状态估计方法，其特征在于，所述状态估计方法包括以下步骤：
步骤一：开路电压-荷电状态数据库的获取；
步骤二：开路电压拟合公式的获取：
U(t)＝γ-a·tb-c·tdln(t)(7)
其中所述U(t)为电池非平衡态的端电压，γ为电池的稳态开路电压，a,b,c,d为拟合参
数，t为时间；
步骤三：获取电池终止工作后端电压随时间变化的U(t)数据；
步骤四：根据步骤二和步骤三得到公式(7)中的γ值；
步骤五：根据步骤四获得的γ值和步骤一获得的开路电压-荷电状态数据库，得到电池
的荷电状态。
2.根据权利要求1所述的一种锂离子电池荷电状态估计方法，其特征在于所述步骤一
中开路电压-荷电状态数据库的获取的具体过程为：
使电池在预设温度条件下静置1h-3h后，用C/100-C/10的倍率全充放循环两次，记录
电池端电压和充放电容量值，对容量值进行归一化处理，获得电压-荷电状态数据，由公式
(1)获得电池开路电压公式(2)，即将充电电压和放电电压取平均，获得的电压值作为电
池平衡态的开路电压值，最终获得开路电压-荷电状态数据库；
其中所述预设温度为-30℃～60℃；
Vc=V0+IΣRiVdis=V0-IΣRi---(1)]]>V0=Vc+Vdis2---(2)]]>其中所述Vc为充电电压，Vdis为放电电压，V0为稳态电压，I为电流，Ri为阻抗。
3.根据权利要求2所述的一种锂离子电池荷电状态估计方法，其特征在于所述步骤二
中稳态开路电压拟合公式的获取过程为：
(1)设计锂离子电池，所述电池由正极、负极、隔膜和电解液组成；
(2)采用有限元法将锂离子电池分割为s段，s取值为∞；
(3)在相邻两段之间设计电极反应，电极反应电势符合Nernst方程，即公式(3)；
Ei(t)=Ei′-RTnFlnciαi(t)(t)ci+1αi+1(t)(t),i=1,2,3...,s-1---(3)]]>其中αi为粒子活度，αi(t)为时间的函数，Ei(t)为电极反应电势，Ei'为形式电势，R为
气体常数，T为温度，F法拉第常数，n为电极反应电子转移数目，ci(t)为第i段粒子浓度；
电池的端电压等于电池两端电极...

【专利技术属性】
技术研发人员：尹鸽平，杨杰，高云智，杜春雨，申斌，崔瑛志，左朋建，程新群，马玉林，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江;23