CN108318823B - 一种基于噪声跟踪的锂电池荷电状态估算方法 - Google Patents
一种基于噪声跟踪的锂电池荷电状态估算方法 Download PDFInfo
- Publication number
- CN108318823B CN108318823B CN201711462522.3A CN201711462522A CN108318823B CN 108318823 B CN108318823 B CN 108318823B CN 201711462522 A CN201711462522 A CN 201711462522A CN 108318823 B CN108318823 B CN 108318823B
- Authority
- CN
- China
- Prior art keywords
- estimation
- noise
- soc
- model
- state
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/385—Arrangements for measuring battery or accumulator variables
- G01R31/387—Determining ampere-hour charge capacity or SoC
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Secondary Cells (AREA)
- Tests Of Electric Status Of Batteries (AREA)
Abstract
本发明涉及一种基于噪声跟踪的锂电池荷电状态估算方法,包括以下步骤:1)离线模型构建,所述离线模型包括开路电压模型和等效电路模型;2)在线估计,利用所述离线模型建立SOC在线估计模型,基于噪声跟踪实现SOC估计。所述在线估计具体包括以下步骤:201)基于电流积分公式和所述离线模型建立非线性状态空间方程;202)结合滚动时域估计策略,建立增广非线性状态空间方程和SOC在线估计模型;203)根据检测电压和电流,利用SOC在线估计模型实现过程噪声估计、测量噪声估计和SOC估计。与现有技术相比,本发明具有通过跟踪过程噪声可减少电流积分法中的电流测量误差,具有锂电池SOC估计精确性和可靠性高等优点。
Description
技术领域
本发明涉及电池管理系统,尤其是涉及一种基于噪声跟踪的锂电池荷电状态估算方法。
背景技术
因具有能量密度大、输出功率高、充放电寿命长等优点,锂离子电池已广泛应用于便携电子设备、电动汽车、家庭储能及空间技术等新兴技术领域。电池荷电状态(State ofCharge,即SOC)估计作为锂电池管理系统的核心功能之一,对于提高电池使用率、延长电池使用寿命、提高电池使用安全性提高至关重要。目前电池管理系统SOC估算中主要采用电流积分法。该方法SOC估计精度主要受限于初始SOC估计误差和电流测量误差两方面。由于电流积分法缺乏消除初始误差的反馈机制,且无法对电流测量噪声进行及时跟踪和校正,因此精度较低,无法完全满足实际需求。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于噪声跟踪的锂电池荷电状态估算方法。
本发明的目的可以通过以下技术方案来实现:
一种基于噪声跟踪的锂电池荷电状态估算方法,包括以下步骤:
1)离线模型构建,所述离线模型包括开路电压模型和等效电路模型;
2)在线估计,利用电流积分公式和所述离线模型建立SOC在线估计模型,基于噪声跟踪实现SOC估计。
所述开路电压模型为开路电压与SOC的函数关系,具体表示为:
其中,VOC为开路电压,SOC为电池荷电状态,j为多项式函数中第j阶,β1j为多项式系数,M为多项式总阶数,下标k为采样时刻。
所述等效电路模型为电路参数与SOC的函数关系,所述电路参数包括开路电压、欧姆内阻以及一阶或多阶RC网络的电阻和电容。
所述等效电路模型为一阶或多阶RC网络的等效电路模型。
所述一阶或多阶RC网络的等效电路模型中,欧姆内阻、RC网络的电阻和电容与SOC的函数关系为:
其中,n为RC网络阶数,R0为欧姆内阻,Rn为第n阶RC网络上的极化电阻,Cn为第n阶RC网络上的等效电容,SOC为电池荷电状态,j为多项式函数中第j阶,β2j、β2n+1,j及β2n+2,j为多项式系数,M为多项式总阶数,下标k为采样时刻。
所述在线估计具体包括以下步骤:
201)基于电流积分公式和所述离线模型建立非线性状态空间方程;
202)结合滚动时域估计策略,建立增广非线性状态空间方程和基于噪声追踪的SOC在线估计模型;
203)根据检测电压和电流,利用所述SOC在线估计模型实现过程噪声估计、测量噪声估计和SOC估计。
所述非线性状态空间方程包括:
状态方程:xk+1=F(xk,uk)+wk
观测方程:yk=h(xk,uk)+vk
其中,状态矢量x=[SOC,V1,…,Vn]T,输入变量u=I,观测变量y=Vb,Vb为电池电压,w和v分别为过程噪声和测量噪声,两者相互独立且均为高斯白噪声,其协方差分别为Qw和R。I为负载电流,Δt为采样周期,VOC为开路电压,C为电池容量,SOC为电池荷电状态,n为RC网络阶数,R0为欧姆内阻,Rn为第n阶RC网络上的极化电阻,Cn为第n阶RC网络上的等效电容,Vn为第n阶RC网络上的电压,τn=RnCn为第n阶RC网络的时间常数,下标k为采样时刻。
步骤202)中,将过程噪声作为状态变量,建立增广非线性状态空间方程,且所述增广非线性状态空间方程的状态方程中F(xk,uk)替换为F(zk,uk),
其中,z=[SOC,V1,…,Vn,w0,…,wn]T为增广状态向量,过程噪声相应转化为γ=[w,θ]T,与测量噪声相互独立且为高斯白噪声,其协方差为Q。
所述SOC在线估计模型表示为:
步骤203)具体包括:
231)初始化;
232)在T时刻对所述SOC在线估计模型进行求解,获得当前状态估计值、过程噪声估计和测量噪声估计值;
233)根据状态方程计算获得T时刻的SOC;
234)更新误差方差矩阵;
235)令T=T+1,构造新的测量数据集yT,返回步骤232)。
所述估计误差方协方差P的更新公式为:
与现有技术相比,本发明具有以下有益效果:
1、本发明方法通过跟踪过程噪声可减少电流积分法中的电流测量误差,通过跟踪测量噪声估计可提高回馈电压对SOC估计值的修正能力,从而保证工业检测环境下,锂电池SOC估计值的精确性和可靠性,最终提升电池管理系统整体性能。
2、本发明建立的离线模型包括开路电压模型和等效电路模型,准确性高,为在线估计提供基础。
附图说明
图1为本发明方法的原理示意图;
图2为本发明实施方式中SOC估算装置的结构;
图3为本发明实施方式中电流激励和电压响应波形图;
图4为本发明实施方式中锂电池等效电路模型图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,本发明提供一种基于噪声跟踪的锂电池荷电状态估算方法,包括以下步骤:1)离线模型构建,所述离线模型包括开路电压模型和等效电路模型;2)在线估计,利用所述离线模型建立SOC在线估计模型,基于噪声跟踪实现SOC估计。在线估计过程具体为:201)基于电流积分公式和所述离线模型建立非线性状态空间方程;202)结合滚动时域估计策略,建立增广非线性状态空间方程和基于噪声跟踪的SOC在线估计模型;203)根据检测电压和电流,利用所述SOC在线估计模型实现过程噪声估计、测量噪声估计和SOC估计。该方法通过跟踪过程噪声可减少电流积分法中的电流测量误差,通过跟踪测量噪声估计可提高回馈电压对SOC估计值的修正能力,从而保证工业检测环境下,锂电池SOC估计值的精确性和可靠性,最终提升电池管理系统整体性能。
上述方法可应用于锂电池管理系统,进行锂电池储能设备的荷电状态估计。在本发明具体实施方式中,锂电池SOC估算装置的结构如图2所示,包括微控制器100、存储器102、电流和电压表104、SOC估算器106。微控制器100总体上控制SOC估算装置,电流和电压表。存储器102用于存储控制器所执行程序。电流和电压表104根据控制器100的控制来测量电流和电压。SOC估算器106根据控制器100的控制估算SOC,以及向控制器100提供估算结果。SOC估算器的建立包括离线模型构建和在线算法应用。
上述基于噪声跟踪的锂电池荷电状态估算方法的具体过程如下:
步骤S11,对所述电池进行开路电压实验,建立开路电压模型。所述电池先以恒流恒压方式充电至截止电压,并静置一定时间;再以1C倍率电流对电池持续放电至特定SOC后,静置1小时。整个过程中,电池端电压及负载电流均以1Hz采样频率同步采集。根据各静置点的SOC及其所对应的开路电压测量值,建立开路电压与SOC的函数关系。本发明具体实施中采用10阶多项式形式表示开路电压与SOC的函数关系:
式中,所需辨识的参数为多项式系数β1i,所需SOC根据电流积分法计算:
式中,SOC(0)为电池初始SOC值,C为电池容量,I为负载电流。结合关系式(1)和(2),采用最小二乘法对β1j进行参数辨识,确定开路电压与SOC的函数关系。
步骤S12,对所述电池进行充放电特性测试,并基于所测数据建立等效电路模型。所述电池先以恒流恒压方式充电至截止电压,并静置一定时间;再对所述电池实施特定充放电测试工况。整个过程中,电池端电压及负载电流均以1Hz采样频率同步采集。本发明具体实施中采用HPPC标准测试工况作为充放电测试工况,如图3所示。等效电路模型包括三部分:开路电压VOC,欧姆内阻R0以及一阶或多阶RC网络,其中RC网络由极化电阻和等效电容组成,开路电压VOC由步骤S11确定。
本发明具体实施中采用一阶等效电路模型,如图4所示,其中Vb为电池电压,I为负载电流。该等效电路模型符合以下电压电流关系:
Vb=VOC-V1-IR0(4)
对于一个采样周期Δt内,关系式(3)和(4)的离散化形式表示可得表示为:
Vb,k=VOC(SOCk)-V1,k-IkR0,k (6)
其中,时间常数τ1=R1C1。在本发明具体实施中,Δt为1s。
在本发明具体实施中,式(5)和(6)中电路参数R0、R1和C1与SOC的函数关系采用6阶多项式形式表示:
其中,所需辨识参数为多项式系数β2j、β3j及β4j。在辨识过程中,基于式(5)和式(6),采用最小二乘法对图2中的电压响应曲线进行拟合,获得β2j、β3j及β4j。至此,离线模型构建完成。
步骤S21,基于电流积分公式和离线模型建立非线性状态空间方程。该状态空间方程可表示为:
状态方程:xk+1=F(xk,uk)+wk (10)
观测方程:yk=h(xk,uk)+vk (11)
满足约束条件:
xk∈X,wk∈W,vk∈V (12)
xk∈[xL,xU] (13)
在本发明中,定义状态矢量为xk=[SOCk,V1,k]T,输入变量为uk=Ik,观测变量为yk=Vb,k。wk和vk分别表示过程噪声和测量噪声,相互独立且均为高斯白噪声,两者协方差分别为Qwk和Rk。状态方程和观测方程中的非线性函数F(xk,uk)和h(xk,uk)分别为:
h(xk,uk)=VOC(SOCk)-V1,k-IkR0,k (15)
式中,VOC、R0、R1和C1与SOC的函数关系由离线模型所得。
步骤S22,为同步实施噪声跟踪及状态估计,将过程噪声转化为状态变量,即:wk+1=wk+θk,θk为协方差Qθk的高斯白噪声,建立增广非线性状态空间方程。定义z=[SOC,V1,w0,w1]T为增广状态向量,式(10)相应转化为:
zk+1=F(zk,uk)+γk (16)
其中,过程噪声γk=[wk,θk]T,与测量噪声相互独立且为高斯白噪声,其协方差为Qk。
其中,z=[SOC,V1,w0,w1]T为增广状态向量。
结合非线性状态空间方程和滚动时域估计策略,建立SOC在线估计模型,整定各算法参数。假设系统初始状态为z0,对于k时刻,所有测量数据为干扰序列为并且z0的先验估计值满足均值为协方差为P0的正态分布。若固定数据时域为L,那么在T时刻,状态估计问题可等价为如下滚动时域优化问题:
满足约束条件(7)-(9),(11),(15)-(17)
0≤SOCk≤1 (19)
上述模型中,为到达代价,采用近似替代。参数R为过程噪声协方差,反映电流积分过程中的电流测量误差和离线模型误差。参数Q为测量噪声协方差,反映电压测量误差。参数P为估计误差协方差,反映对初始估计的信心。通过求解该问题,可得到当前时刻SOC估计值、过程噪声估计值和测量噪声估计值。在本发明具体实施中,L综合估计精度和精算时间进行优选。Q,R根据各状态变量数量级进行整定。
式中P采用下式更新:
式中各项定义如下:
至此,优化目标(16)可表示为:
步骤S23,根据检测电压和电流,利用在线估计模型,实现过程噪声估计、测量噪声估计和SOC估计。在采用非线性滚动时域法进行估算时,包括以下五个步骤:
4、根据式(20)计算下一时刻误差方差矩阵PT-L;
5、在T+1时刻,测量yT,构造新的测量数据集,返回步骤2。
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。
Claims (8)
1.一种基于噪声跟踪的锂电池荷电状态估算方法,其特征在于,包括以下步骤:
1)离线模型构建,所述离线模型包括开路电压模型和等效电路模型;
2)在线估计,利用所述离线模型建立SOC在线估计模型,基于噪声跟踪实现SOC估计,所述在线估计具体包括以下步骤:
201)基于电流积分公式和所述离线模型建立非线性状态空间方程;
202)结合滚动时域估计策略,建立增广非线性状态空间方程和基于噪声跟踪的SOC在线估计模型;
203)根据检测电压和电流,利用所述SOC在线估计模型实现过程噪声估计、测量噪声估计和SOC估计;
步骤203)具体包括:
231)初始化;
232)在T时刻对所述SOC在线估计模型进行求解,获得当前状态估计值、过程噪声估计和测量噪声估计值;
233)根据状态方程计算获得T时刻的SOC和噪声;
234)更新估计误差协方差;
235)令T=T+1,构造新的测量数据集yT,返回步骤232)。
3.根据权利要求1所述的基于噪声跟踪的锂电池荷电状态估算方法,其特征在于,所述等效电路模型为电路参数与SOC的函数关系,所述电路参数包括开路电压、欧姆内阻以及一阶或多阶RC网络的电阻和电容。
5.根据权利要求1所述的基于噪声跟踪的锂电池荷电状态估算方法,其特征在于,所述非线性状态空间方程包括:
状态方程:xk+1=F(xk,uk)+wk
观测方程:yk=h(xk,uk)+vk
其中,状态矢量x=[SOC,V1,…,Vn]T,输入变量u=I,观测变量y=Vb,Vb为电池电压;w和v分别为过程噪声和测量噪声,两者相互独立且均为高斯白噪声,其协方差分别为Qw和R;I为负载电流,Δt为采样周期,VOC为开路电压,C为电池容量,SOC为电池荷电状态,n为RC网络阶数,R0为欧姆内阻,Rn为第n阶RC网络上的极化电阻,Cn为第n阶RC网络上的等效电容,Vn为第n阶RC网络上的电压,τn=RnCn为第n阶RC网络的时间常数,下标k为采样时刻。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711462522.3A CN108318823B (zh) | 2017-12-28 | 2017-12-28 | 一种基于噪声跟踪的锂电池荷电状态估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711462522.3A CN108318823B (zh) | 2017-12-28 | 2017-12-28 | 一种基于噪声跟踪的锂电池荷电状态估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108318823A CN108318823A (zh) | 2018-07-24 |
CN108318823B true CN108318823B (zh) | 2020-06-02 |
Family
ID=62892634
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711462522.3A Active CN108318823B (zh) | 2017-12-28 | 2017-12-28 | 一种基于噪声跟踪的锂电池荷电状态估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108318823B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
PL3839532T3 (pl) * | 2019-02-15 | 2023-02-13 | Lg Energy Solution, Ltd. | Urządzenie i sposób szacowania soc |
CN112379270B (zh) * | 2020-11-13 | 2024-01-30 | 哈尔滨工业大学 | 一种电动汽车动力电池荷电状态滚动时域估计方法 |
CN114252797B (zh) * | 2021-12-17 | 2023-03-10 | 华中科技大学 | 一种基于不确定性估计的锂电池剩余使用寿命预测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5694088B2 (ja) * | 2011-08-23 | 2015-04-01 | トヨタ自動車株式会社 | 二次電池の劣化管理システム |
CN106154168A (zh) * | 2016-04-01 | 2016-11-23 | 储盈新能源科技(上海)有限公司 | 数据驱动的动力电池荷电状态估计方法 |
CN106842060A (zh) * | 2017-03-08 | 2017-06-13 | 深圳市海云图新能源有限公司 | 一种基于动态参数的动力电池soc估算方法及系统 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102289557B (zh) * | 2011-05-17 | 2013-08-07 | 杭州电子科技大学 | 一种电池模型参数与剩余电量联合异步在线估计方法 |
CN104007395B (zh) * | 2014-06-11 | 2016-08-24 | 北京交通大学 | 锂离子电池荷电状态与参数自适应联合估计方法 |
EP2963434B1 (en) * | 2014-06-30 | 2021-08-11 | Foundation Of Soongsil University-Industry Cooperation | Battery state estimation method and system using dual extended kalman filter, and recording medium for performing the method |
CN105572596B (zh) * | 2016-01-20 | 2019-02-19 | 上海交通大学 | 锂电池soc估算方法及系统 |
CN105607009B (zh) * | 2016-02-01 | 2018-05-01 | 深圳大学 | 一种基于动态参数模型的动力电池soc估计方法和系统 |
CN106324521B (zh) * | 2016-09-05 | 2018-09-11 | 北京理工大学 | 一种联合估计动力电池系统参数与荷电状态的方法 |
CN106646265A (zh) * | 2017-01-22 | 2017-05-10 | 华南理工大学 | 一种锂电池soc估计方法 |
CN107037374B (zh) * | 2017-05-04 | 2019-10-08 | 广东电网有限责任公司电力科学研究院 | 一种锂离子电池的soc和soh联合估计方法 |
-
2017
- 2017-12-28 CN CN201711462522.3A patent/CN108318823B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5694088B2 (ja) * | 2011-08-23 | 2015-04-01 | トヨタ自動車株式会社 | 二次電池の劣化管理システム |
CN106154168A (zh) * | 2016-04-01 | 2016-11-23 | 储盈新能源科技(上海)有限公司 | 数据驱动的动力电池荷电状态估计方法 |
CN106842060A (zh) * | 2017-03-08 | 2017-06-13 | 深圳市海云图新能源有限公司 | 一种基于动态参数的动力电池soc估算方法及系统 |
Non-Patent Citations (2)
Title |
---|
Online Parameter Identification and State of Charge Estimation of Lithium-Ion Batteries Based on Forgetting Factor Recursive Least Squares and Nonlinear Kalman Filter;Bizhong Xia 等;《Energies》;20171221;全文 * |
基于模型的锂离子电池SOC估计研究;郭凯;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20140315(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108318823A (zh) | 2018-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108008320B (zh) | 一种锂离子电池荷电状态与模型参数自适应联合估计方法 | |
Chen et al. | Online state of charge estimation of Li-ion battery based on an improved unscented Kalman filter approach | |
CN110488194B (zh) | 一种基于电化学阻抗模型的锂电池soc估算方法及其系统 | |
Wang et al. | An adaptive working state iterative calculation method of the power battery by using the improved Kalman filtering algorithm and considering the relaxation effect | |
CN105116343B (zh) | 最小二乘支持向量机的动力电池电荷状态估计方法及系统 | |
CN105607009B (zh) | 一种基于动态参数模型的动力电池soc估计方法和系统 | |
CN104267261B (zh) | 基于分数阶联合卡尔曼滤波的二次电池简化阻抗谱模型参数在线估计方法 | |
CN106842060A (zh) | 一种基于动态参数的动力电池soc估算方法及系统 | |
CN109061506A (zh) | 基于神经网络优化ekf的锂离子动力电池soc估计方法 | |
CN103344921B (zh) | 锂离子动力电池健康状态评估系统及方法 | |
CN108318823B (zh) | 一种基于噪声跟踪的锂电池荷电状态估算方法 | |
JP2016529478A (ja) | バッテリに蓄積された電荷のモニタリング | |
CN103616647A (zh) | 一种用于电动汽车电池管理系统的电池剩余电量估计方法 | |
CN112989690B (zh) | 一种混合动力汽车锂电池多时间尺度荷电状态估计方法 | |
US9891287B2 (en) | Temperature-compensated state of charge estimation for rechargeable batteries | |
Wang et al. | Adaptive state-of-charge estimation method for an aeronautical lithium-ion battery pack based on a reduced particle-unscented kalman filter | |
CN109085505A (zh) | 一种动力电池充放电状态估算方法 | |
CN112305426B (zh) | 一种多约束条件下的锂离子电池功率状态估计系统 | |
CN110716146A (zh) | 一种动力电池开路电压的估计方法 | |
CN113011012A (zh) | 一种基于Box-Cox变化的储能电池剩余寿命预测方法 | |
CN111060822B (zh) | 一种基于模型切换及融合的荷电状态估计方法 | |
CN113466725B (zh) | 确定电池荷电状态的方法及装置,存储介质及电子设备 | |
Wang et al. | Lithium-ion battery security guaranteeing method study based on the state of charge estimation | |
CN109782182B (zh) | 一种串联电池组能量状态在线估计方法和装置 | |
CN114740385A (zh) | 一种自适应的锂离子电池荷电状态估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |