CN105842627A - 基于数据模型融合的动力电池容量和荷电状态的估计方法 - Google Patents

基于数据模型融合的动力电池容量和荷电状态的估计方法 Download PDF

Info

Publication number
CN105842627A
CN105842627A CN201610069076.9A CN201610069076A CN105842627A CN 105842627 A CN105842627 A CN 105842627A CN 201610069076 A CN201610069076 A CN 201610069076A CN 105842627 A CN105842627 A CN 105842627A
Authority
CN
China
Prior art keywords
electrokinetic cell
charge
described electrokinetic
state
max
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.)
Granted
Application number
CN201610069076.9A
Other languages
English (en)
Other versions
CN105842627B (zh
Inventor
熊瑞
杨瑞鑫
何洪文
孙逢春
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201610069076.9A priority Critical patent/CN105842627B/zh
Publication of CN105842627A publication Critical patent/CN105842627A/zh
Application granted granted Critical
Publication of CN105842627B publication Critical patent/CN105842627B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/385Arrangements for measuring battery or accumulator variables
    • G01R31/387Determining ampere-hour charge capacity or SoC

Abstract

本发明涉及动力电池参数和状态的估计领域,尤其涉及一种基于数据模型融合的动力电池容量和荷电状态的估计方法。为提高电动车辆动力电池荷电状态的估计精度,本发明提出一种基于数据模型融合的动力电池容量和荷电状态的估计方法,对动力电池进行老化实验,得出由动力电池在不同老化程度下的开路电压、荷电状态和可用容量形成的三维响应面关系,建立动力电池的等效电路模型,输入动力电池的实测端电压和充放电电流到动力电池的模型方程中,辨识出动力电池的等效电路模型的模型参数以及动力电池的最大可用容量和荷电状态初始值,估计出动力电池的荷电状态。本发明方法可对不同老化程度的动力电池的容量和荷电状态进行估计,实用性强,估计精度较高。

Description

基于数据模型融合的动力电池容量和荷电状态的估计方法
技术领域
本发明涉及动力电池参数和状态的估计领域,尤其涉及一种基于数据模型融合对电动车辆的车载动力电池的容量和荷电状态进行估计的方法。
背景技术
以电动车辆为主的新能源汽车是我国战略新兴产业之一,也是“中国制造2025”重点推进的领域之一。动力电池是电动车辆的核心部件,因此,动力电池的参数及状态的估计精准度直接影响电动车辆应用的安全性和可靠性。
近年来,动力电池的荷电状态(state of charge,简称SOC)的估计一直是学术界和工业界研究的热点,并且已经取得了一定的成就。但是,目前,对动力电池的荷电状态SOC进行估计时,多是采用确定的最大可用容量对动力电池的荷电状态进行估计,而在实际使用过程中,动力电池的最大可用容量会随着动力电池的老化而逐渐退化,所以,利用确定的最大可用容量对动力电池的荷电状态SOC进行估计,难以保证估计结果的合理性和可靠性,进而易因荷电状态SOC估计不准确而导致动力电池过充、过放以及滥用,甚至发生安全故障。
发明内容
为提高电池车辆的动力电池的荷电状态的估计精度低,本发明提出一种基于数据模型融合的动力电池容量和荷电状态的估计方法,该方法包括如下步骤:
步骤1,建立由所述动力电池在不同的老化程度下的开路电压、荷电状态和可用容量形成的三维响应面
在相同温度下,对所述动力电池进行老化实验,采集所述动力电池在不同老化程度下的开路电压UOC以及与该开路电压UOC相对应的荷电状态z和可用容量,建立由所述动力电池在不同老化程度程度下的开路电压UOC以及与该开路电压UOC相对应的荷电状态z和可用容量形成的三维响应面,且UOC=f(z,Qmax),
其中,
f为所述动力电池的开路电压函数,
Qmax为所述动力电池的最大可用容量;
步骤2,建立所述动力电池的等效电路模型;
步骤3,估计所述动力电池的最大可用容量Qmax和荷电状态z
首先,根据所述动力电池的等效电路模型建立所述动力电池的模型方程;
然后,在所述动力电池充放电过程中,在一段连续时间段T内对所述动力电池的实测端电压U和充放电电流I进行采样,采样时间间隔为Δt;
接着,将z=g(z0,Qmax)和UOC=f(z,Qmax)带入到所述动力电池的模型方程中,并将采样得到的实测端电压U和充放电电流I输入到所述动力电池的模型方程中,辨识出所述动力电池的等效电路模型的模型参数以及所述动力电池的最大可用容量Qmax和荷电状态初始值z0
其中,g为所述动力电池的荷电状态函数;
最后,根据所述动力电池的最大可用容量Qmax和荷电状态初始值z0估计出所述动力电池的荷电状态z。
本发明方法对动力电池的最大可用容量和荷电状态进行离线估计时,利用由动力电池在不同老化程度下的开路电压以及与该开路电压相对应的荷电状态、可用容 量形成的三维响应面和动力电池的等效电路模型配合,并在对动力电池的等效电路模型的模型参数进行辨识时,直接辨识得出动力电池的最大可用容量和荷电状态初始值,并根据动力电池的最大可用容量和荷电状态初始值估计出动力电池的荷电状态,计算简单,计算量小,且估计精度较高。
优选地,在所述步骤1中,采集所述动力电池的可用容量、开路电压UOC及荷电状态z时的采样时间间隔Δt'为定值,便于采集数据,且采集到的动力电池的可用容量及荷电状态z呈等量递减,便于绘制三维响应面。
优选地,在对所述动力电池进行老化实验时,所述动力电池所处老化程度为循环0次、失效状态和失效的中间状态,对多种老化程度下的动力电池进行老化实验,提高制得的三维响应面的准确性。
优选地,在所述步骤2中,所建立的所述动力电池的等效电路模型为具有一状态滞后的一阶RC网络等效电路模型,且所述模型方程为:
其中,
UD,k为所述动力电池在充放电过程中k时刻的极化电压,
UD,k+1为所述动力电池在充放电过程中k+1时刻的极化电压,
τD为所述动力电池的等效电路模型中的RC网络的时间常数,且τD=RDCD
RD为所述动力电池的等效电路模型中的RC网络的极化电阻,
CD为所述动力电池的等效电路模型中的RC网络的极化电容,
κ为衰减因子,
Ik为所述动力电池在充放电过程中k时刻的充放电电流,
hk为所述动力电池在充放电过程中k时刻的滞后电压,
hk+1为所述动力电池在充放电过程中k+1时刻的滞后电压,
H为所述动力电池的滞后电压的最大值,
Uk为所述动力电池在充放电过程中k时刻的实测端电压,
Ro为所述动力电池的欧姆内阻。
采用具有一状态滞后的RC网络等效电路作为动力电池的等效电路,将动力电池在充放电过程中因充放电电流的大小变化引起的滞后电压引入到动力电池的等效电路中,建立的动力电池的等效电路模型更为准确。
优选地,在所述步骤3中,
其中,t为所述动力电池进行充放电时的充放电时间。
优选地,在所述步骤3中,采用遗传算法对由所述动力电池的等效电路模型的模型参数以及所述动力电池的最大可用容量Qmax和荷电状态初始值z0形成的优化变量进行辨识,其中,
为所述动力电池充电时的欧姆内阻,
为所述动力电池放电时的欧姆内阻,
H+为所述动力电池充电时的滞后电压,
H-为所述动力电池放电时的滞后电压。
优选地,在进行辨识时,设定目标函数
其中,UL为所述动力电池的估计端电压,
n为采样得到的所述动力电池的实测端电压U和充放电电流I数据的长度,且
设定约束条件为
其中,
Cmin为所述动力电池的最大可用容量Qmax的最小值,Cmax为所述动力电池的最大可用容量Qmax的最大值,
为所述动力电池的充电时的欧姆内阻的最小值,为所述动力电池充电时的欧姆内阻的最大值,
为所述动力电池放电时的欧姆内阻的最小值,为所述动力电池放电时的欧姆内阻的最大值,
κmin为所述衰减因子κ的最小值,κmax为所述衰减因子κ的最大值,
为所述动力电池充电时的滞后电压H+的最小值,为所述动力电池充电时的滞后电压H+的最大值,
为所述动力电池放电时的滞后电压H-的最小值,为所述动力电池放电时的滞后电压H-的最大值,
RD min为所述动力电池的RC网络的极化电阻RD的最小值,RD max为所述动力电池的RC网络的极化电阻RD的最大值,
τD min为所述动力电池的RC网络的时间常数τD的最小值,τD max为所述动力电池的RC网络的时间常数τD的最大值。
优选地,采用安时积分法根据所述动力电池的最大可用容量Qmax和荷电状态初始值z0估计出所述动力电池的荷电状态z,估计步骤简单、方便。
使用本发明方法对动力电池的最大可用容量和荷电状态进行估计时,采用由动力电池在不同老化程度下的开路电压、荷电状态及可用容量形成的三维响应面和动力电池的等效电路模型配合,通过参数辨识得出动力电池的最大可用容量和荷电状态初始值,并根据动力电池的最大可用容量和荷电状态初始值估计出动力电池的荷电状态,计算简单,计算量小,且估计得出的最大可用容量的相对误差在3%以内,荷电状态初始值的相对误差在2%以内,荷电状态的相对误差在3%以内,估计精度高于现行的误差标准5%,估计精度高。另外,采用本发明方法可对处于不同老化程度下的动力电池的最大可用容量和荷电状态进行估计,具有普适性,实用性强。
附图说明
图1为本发明估计动力电池的最大可用容量和荷电状态初始值的流程图;
图2为动力电池具有一状态滞后的一阶RC网络等效电路图;
图3为由试验用动力电池的开路电压与相应的荷电状态、最大可用容量形成的三维响应面图;
图4为循环0次的试验用动力电池在UDDS工况下的充放电电流随时间变化的曲线图;
图5为循环0次的试验用动力电池在UDDS工况下的实测端电压随时间变化的曲线图;
图6为从图5中截取的1000s-4600s时间段内的实测端电压随时间变化的曲线图;
图7为循环0次的试验用动力电池在图6所示对应时间段内的估计端电压随时间变化的曲线图;
图8为循环0次的试验用动力电池在图6所示对应时间段内的荷电状态的真实值和估计值随时间变化的曲线;
图9为循环200次的试验用动力电池在UDDS工况下的充放电电流随时间变化的曲线图;
图10为循环200次的试验用动力电池在UDDS工况下的实测端电压随时间变化的曲线图;
图11为从图9中截取的1000s-4600s时间段内的实测端电压随时间变化的曲线图;
图12为循环200次的试验用动力电池在图11所示对应时间段内的估计端电压随时间变化的曲线图;
图13为循环200次的试验用动力电池在图11所示对应时间段内的荷电状态的真实值和估计值随时间变化的曲线。
具体实施方式
下面结合图1-13对本发明基于数据模型融合的动力电池容量和荷电状态的估计方法进行详细说明。
如图1所示,对电动车辆的动力电池进行老化实验,得出该动力电池在不同的老化程度下的开路电压、荷电状态及最大可用容量之间的对应关系,建立由动力电池在不同老化程度下的开路电压以及与该开路电压相对应的荷电状态和可用容量形成的三维响应曲面。建立动力电池的等效电路模型,并根据设定的目标函数、约束条件辨识出动力电池的等效电路模型的模型参数、动力电池的荷电状态初始值z0及最大可用容量Qmax,从而利用辨识得到的动力电池的最大可用容量Qmax和荷电状态初始值z0估计出动力电池的荷电状态z。
具体步骤如下:
步骤1,建立由动力电池在不同老化程度下的开路电压和与该开路电压相对应的荷电状态、可用容量形成的三维响应面
在相同温度下,对电动车辆的动力电池进行老化实验,并对该动力电池处于不同老化程度的开路电压UOC以及与该开路电压UOC相对应的荷电状态z和可用容量Q进行采集,且采样时间间隔为Δt'为定值。动力电池的老化程度可表示为循环0次,失效状态和失效的中间状态,其中,循环0次的动力电池指的是新电池;处于失效状态的动力电池指的是当前状态的最大可用容量低于或等于其额定容量的80%的动力电池;处于失效的中间状态的动力电池指的是当前状态的最大可用容量大约为且额定容量的90%的动力电池。在进行老化实验时,动力电池的放电电流I的大小保持不变,每经过一个采样时间间隔Δt'后,动力电池中的可用容量Q从采样起始时刻对应的最大可用容量Qmax开始呈等量递减。也就是说,在对该动力电池进行老化实验时,采集得到的该动力电池的荷电状态z从采样起始时刻对应的荷电状态初始值z0开始呈等量递减。根据采样数据建立由动力电池在不同老化程度的开路电压UOC以及与该开路电压UOC相对应的荷电状态z和可用容量Q形成的三维响应面,且UOC=f(z,Qmax),其中,f为动力电池的开路电压函数。
步骤2,建立动力电池的等效电路模型
由于动力电池在进行充电和放电的过程中,存在一定的迟滞现象,故选用如图2所示的具有一状态滞后的一阶RC网络等效电路作为动力电池的等效电路模型。该等效电路除引入用来描述在充放电电流变化的时候引起的电压滞后现象的滞后电压h外,还包括电压源-OCV、欧姆内阻-Ro和RC网络三部分,其中,电压源-OCV为动力电池的开路电压UOC;欧姆内阻-Ro表示动力电池中电极材料、电解液、隔膜电阻及其他零件的接触电阻;RC网络使用极化内阻RD和极化电容CD来描述动力电池的 动态特性,该动态特性包括动力电池的极化特性和扩散效应,RC网络的时间常数τD=RDCD
步骤3,估计出动力电池的最大可用容量和荷电状态
根据动力电池的等效电路模型建立动力电池的模型方程
其中,
UD,k为动力电池在充放电过程中k时刻的极化电压,
UD,k+1为动力电池在充放电过程中k+1时刻的极化电压,
κ为衰减因子,用来调整滞后电压h的大小,
Ik为动力电池在充放电过程中k时刻的充放电电流,
hk为动力电池在充放电过程中k时刻的滞后电压,
hk+1为动力电池在充放电过程中k+1时刻的滞后电压,
H为动力电池的滞后电压的最大值,
Uk为动力电池在充放电过程中k时刻的实测端电压。
在动力电池充放电过程中,在一段连续时间段T内对动力电池的实测端电压V和充放电电流I进行采样,采样时间间隔为Δt。
由于动力电池的荷电状态z可表示为z=g(z0,Qmax),
其中,g为动力电池的荷电状态函数,
故可将z=g(z0,Qmax)和UOC=f(z,Qmax)带入到动力电池的模型方程中,并将动力电池的最大可用容量Qmax和荷电状态SOC的初始值z0及该动力电池的等效电路模型的模型参数共同设置为优化变量χ,且其中,
为动力电池充电时的欧姆内阻,
为动力电池放电时的欧姆内阻,
H+为动力电池充电时的滞后电压,
H-为动力电池放电时的滞后电压。
优选地,其中,t为动力电池充放电时间。
将采样得到的动力电池的实测端电压U和充放电电流I输入到动力电池的模型方程中,并采用遗传算法对优化变量χ进行辨识,从而辨识得出动力电池的等效电路模型的模型参数以及动力电池的最大可用容量Qmax和荷电状态初始值z0,进而根据动力电池的最大可用容量Qmax和荷电状态初始值z0估计出动力电池在充放电过程中的荷电状态z。在进行辨识时,设定目标函数F使动力电池的估计端电压与该动力电池的实测端电压之间的误差的平方和最小,即
其中,
UL为动力电池的估计端电压,
n为采样得到的动力电池的实测端电压U和充放电电流I数据的长度,且
设定约束条件为
其中,
Cmin为动力电池的最大可用容量Qmax的最小值,Cmax为动力电池的最大可用容量Qmax的最大值,
为动力电池的充电时的欧姆内阻的最小值,为动力电池充电时的欧姆内阻的最大值,
为动力电池放电时的欧姆内阻的最小值,为动力电池放电时的欧姆内阻的最大值,
κmin为衰减因子κ的最小值,κmax为衰减因子κ的最大值,
为动力电池充电时的滞后电压H+的最小值,为动力电池充电时的滞后电压H+的最大值,
为动力电池放电时的滞后电压H-的最小值,为动力电池放电时的滞后电压H-的最大值,
RD min为动力电池的RC网络的极化电阻RD的最小值,RD max为动力电池的RC网络的极化电阻RD的最大值,
τD min为动力电池的RC网络的时间常数τD的最小值,τD max为动力电池的RC网络的时间常数τD的最大值。
约束条件中的相关参数的取值范围是根据动力电池本身的特性以及动力电池生产厂家提供的技术参数确定的。比如,动力电池的最大可用容量Qmax的取值范围根据生产厂家提供的技术参数确定的,通常不允许动力电池的最大可用容量降低至其标称容量Q标称的80%;动力电池的荷电状态的最佳工作区间为0.3~1,故在约束条件中,动力电池荷电状态初始值的取值范围可确定为0.1~1;动力电池的欧姆电阻Ro的取值范围根据该动力电池本身的特性以及生产厂家提供的技术参数计算确定。
在对优化变量χ进行辨识时,也可以采用其他方法进行辨识,比如粒子群优化算法。
在根据动力电池的最大可用容量Qmax和荷电状态初始值z0对动力电池在充放电过程中的荷电状态z进行估计时,可采用安时积分方法对动力电池的荷电状态z进行估计,
下面以标称容量Q标称为25Ah、上截止电压为4.1伏特(V)、下截止电压为3.0V的三元聚合物锂电池作为试验对象,验证本发明在估计电动车辆上的动力电池的最大可用容量及荷电状态时存在的优势。
首先,对待试验的动力电池进行老化实验,实验温度为25℃。实验过程中,分别对该动力电池在不同老化程度下的开路电压及与该开路电压相对应的荷电状态和可用容量进行采集,且采集过程中,采样时间间隔Δt'为定值,并根据采集得到的数据绘制三维响应面图。具体地,在动力电池的老化程度分别为循环0次、循环100次、循环200次和循环300次时对该动力电池的开路电压、荷电状态及可用容量进行采集,且采样时间间隔Δt'为1s,采集到的部分数据如表1所示,并根据采集得到 的数据绘制形成如图3所示的三维响应面图。
表1
接着,分别以老化程度分别为循环次数为0次和循环次数为200次的动力电池为例,说明采用本发明方法对动力电池的最大可用容量和荷电状态进行估计时的估计效果,估计采用的约束条件中,且
Cmin=80%Q标称,Cmax=110%Q标称
R min + = 10 - 4 Ω , R max + = 10 - 2 Ω ;
R min - = 10 - 4 Ω , R max - = 10 - 2 Ω ;
κmin=-0.1,κmax=0.1;
H m i n + = - 0.1 V , H m a x + = 0.1 V ;
H min - = - 0.1 V , H max - = 0.1 V ;
RD min=10-4Ω,RD max=10-2Ω;
τD min=0,τD max=300。
Eg1.老化程度为循环0次的动力电池
对该动力电池进行城市道路循环(Urban Dynamometer Driving Schedule,简称UDDS)工况测试,且环境温度为25℃,在测试过程中,电动车辆的电池能量管理系统记录的该动力电池的放电电流I和实测端电压U随时间变化的曲线分别如图4和5所示。从中截取时间点范围为1000s-4600s,即一个小时的时间段内的放电电流及端电压数据,且采用时间间隔Δt为1s,且实测端电压U随时间变化的曲线如图6所示。将采样得出的放电电流I和实测端电压U输入到该动力电池具有一状态滞后的一阶RC网络等效电路模型的模型方程中,利用遗传算法辨识得出的动力电池的最大可用容量Qmax和荷电状态初始值z0如表2所示,并得到该动力电池的估计端电压如图7所示。
此外,通过实验测试得出,在测试时间点为1000s即采样起始时刻,该动力电池的可用容量为27.48Ah,荷电状态为85.5%,即该动力电池的最大可用容量Qmax的实测值为27.48Ah,荷电状态初始值z0的实测值为85.5%。
表2
采用安时积分法根据辨识得出的动力电池的最大可用容量和荷电状态初始值对动力电池在放电过程中的荷电状态进行估计,估计结果随该动力电池的放电时间变化的曲线如图8中的虚线所示。而根据电动车辆的电池能量管理系统记录的该动力电池的荷电状态即动力电池的荷电状态的实测值在相应时间段内随时间变化的曲线如图8中的实线所示。经计算可知,在估计时间段内,该动力电池的荷电状态的估计值相对于该动力电池的荷电状态的实测值的最大误差为0.6%。
Eg2.老化程度为循环200次的动力电池
对该动力电池进行UDDS工况测试,且环境温度为25℃,在测试过程中,电动车辆的电池能量管理系统记录的该动力电池的放电电流I和实测端电压U随时间变化的曲线分别如图9和10所示。从中截取时间点范围为1000s-4600s,即一个小时的时间段内的放电电流及端电压数据,且采用时间间隔Δt为1s,且实测端电压U随时间变化的曲线如图11所示。将采样得出的放电电流I和实测端电压U输入到该动力电池具有一状态滞后的一阶RC网络等效电路模型的模型方程中,利用遗传算法辨识得出的动力电池的最大可用容量Qmax和荷电状态初始值z0如表3所示,并得到该动力电池估计端电压如图12所示。
此外,通过实验测试得出,在测试时间电为1000s即采样起始时刻,该动力电池的可用容量为25.21Ah,荷电状态为93.4%,即该动力电池的最大可用容量Qmax的实测值为25.21Ah,荷电状态初始值z0的实测值为93.4%。
表3
采用安时积分法根据辨识得出的动力电池的最大可用容量和荷电状态初始值对动力电池在放电过程中的荷电状态进行估计,估计结果随该动力电池的放电时间变化的曲线如图13中的虚线所示。而根据电动车辆的电池能量管理系统记录的该动力电池的荷电状态即动力电池的荷电状态的实测值在相应时间段内随时间变化的曲线如图13中的实线所示。经计算可知,在估计时间段内,该动力电池的荷电状态的估计值相对于该动力电池的荷电状态的实测值的最大误差为2.3%。
综上可见,本发明所提出的估计动力电池的最大可用容量Qmax和荷电状态初始值z0的方法具有如下优势:
1、该估计方法估计得到的动力电池的最大可用容量Qmax的相对误差在3%以内,荷电状态初始值z0的相对误差在2%以内,荷电状态z的相对误差在3%以内,估计精度较高;
2、该估计方法可对处于不同老化程度下的动力电池的最大可用容量Qmax和荷电状态初始值z0进行估计,且估计精度较高,具有一定的普适性,实用性强。

Claims (8)

1.一种基于数据模型融合的动力电池容量和荷电状态的估计方法,其特征在于,该方法包括如下步骤:
步骤1,建立由所述动力电池在不同的老化程度下的开路电压、荷电状态和可用容量形成的三维响应面
在相同温度下,对所述动力电池进行老化实验,采集所述动力电池在不同老化程度下的开路电压UOC以及与该开路电压UOC相对应的荷电状态z和可用容量,建立由所述动力电池在不同老化程度程度下的开路电压UOC以及与该开路电压UOC相对应的荷电状态z和可用容量形成的三维响应面,且UOC=f(z,Qmax),
其中,
f为所述动力电池的开路电压函数,
Qmax为所述动力电池的最大可用容量;
步骤2,建立所述动力电池的等效电路模型;
步骤3,估计所述动力电池的最大可用容量Qmax和荷电状态z
首先,根据所述动力电池的等效电路模型建立所述动力电池的模型方程;
然后,在所述动力电池充放电过程中,在一段连续时间段T内对所述动力电池的实测端电压U和充放电电流I进行采样,采样时间间隔为Δt;
接着,将z=g(z0,Qmax)和UOC=f(z,Qmax)带入到所述动力电池的模型方程中,并将采样得到的实测端电压U和充放电电流I输入到所述动力电池的模型方程中,辨识出所述动力电池的等效电路模型的模型参数以及所述动力电池的最大可用容量Qmax和荷电状态初始值z0
其中,g为所述动力电池的荷电状态函数;
最后,根据所述动力电池的最大可用容量Qmax和荷电状态初始值z0估计出所述动力电池的荷电状态z。
2.根据权利要求1所述的基于数据模型融合的动力电池容量和荷电状态的估计方法,其特征在于,在所述步骤1中,采集所述动力电池的可用容量、开路电压UOC及荷电状态z时的采样时间间隔Δt'为定值。
3.根据权利要求2所述的基于数据模型融合的动力电池容量和荷电状态的估计方法,其特征在于,在对所述动力电池进行老化实验时,所述动力电池所处老化程度为循环0次、失效状态和失效的中间状态。
4.根据权利要求1-3中任意一项所述基于数据模型融合的动力电池容量和荷电状态的估计方法,其特征在于,在所述步骤2中,所建立的所述动力电池的等效电路模型为具有一状态滞后的一阶RC网络等效电路模型,且所述模型方程为: U D , k + 1 = exp ( - Δ t / τ D ) U D , k + R D [ 1 - exp ( - | κ I k Δ t | ) ] I k h k + 1 = exp ( - | κI k Δ t | ) h k + [ 1 - exp ( - | κI k Δ t | ) ] H U k = f ( z , Q max ) - R o I k - U D , k + h k ,
其中,
UD,k为所述动力电池在充放电过程中k时刻的极化电压,
UD,k+1为所述动力电池在充放电过程中k+1时刻的极化电压,
τD为所述动力电池的等效电路模型中的RC网络的时间常数,且τD=RDCD
RD为所述动力电池的等效电路模型中的RC网络的极化电阻,
CD为所述动力电池的等效电路模型中的RC网络的极化电容,
κ为衰减因子,
Ik为所述动力电池在充放电过程中k时刻的充放电电流,
hk为所述动力电池在充放电过程中k时刻的滞后电压,
hk+1为所述动力电池在充放电过程中k+1时刻的滞后电压,
H为所述动力电池的滞后电压的最大值,
Uk为所述动力电池在充放电过程中k时刻的实测端电压,
Ro为所述动力电池的欧姆内阻。
5.根据权利要求4所述的基于数据模型融合的动力电池容量和荷电状态的估计方法,其特征在于,在所述步骤3中,
其中,t为所述动力电池进行充放电时的充放电时间。
6.根据权利要求4所述的基于数据模型融合的动力电池容量和荷电状态的估计方法,其特征在于,在所述步骤3中,采用遗传算法对由所述动力电池的等效电路模型的模型参数以及所述动力电池的最大可用容量Qmax和荷电状态初始值z0形成的优化变量进行辨识,其中,
为所述动力电池充电时的欧姆内阻,
为所述动力电池放电时的欧姆内阻,
H+为所述动力电池充电时的滞后电压,
H-为所述动力电池放电时的滞后电压。
7.根据权利要求6所述的基于数据模型融合的动力电池容量和荷电状态的估计方法,其特征在于,在进行辨识时,设定目标函数
其中,UL为所述动力电池的估计端电压,
n为采样得到的所述动力电池的实测端电压U和充放电电流I数据的长度且 n = T Δ t + 1 ;
设定约束条件为 s . t . C m i n ≤ Q m a x ≤ C m a x 0.1 ≤ z 0 ≤ 1 R m i n + ≤ R o + ≤ R m a x + R min - ≤ R o - R max - κ m i n ≤ κ ≤ κ m a x H min + ≤ H + ≤ H m a x + H min - ≤ H - ≤ H max - R D m i n ≤ R D ≤ R D m a x τ D m i n ≤ τ D ≤ τ D m a x ,
其中,
Cmin为所述动力电池的最大可用容量Qmax的最小值,Cmax为所述动力电池的最大可用容量Qmax的最大值,
为所述动力电池的充电时的欧姆内阻的最小值,为所述动力电池充电时的欧姆内阻的最大值,
为所述动力电池放电时的欧姆内阻的最小值,为所述动力电池放电时的欧姆内阻的最大值,
κmin为所述衰减因子κ的最小值,κmax为所述衰减因子κ的最大值,
为所述动力电池充电时的滞后电压H+的最小值,为所述动力电池充电时的滞后电压H+的最大值,
为所述动力电池放电时的滞后电压H-的最小值,为所述动力电池放电时的滞后电压H-的最大值,
RDmin为所述动力电池的RC网络的极化电阻RD的最小值,RDmax为所述动力电池的RC网络的极化电阻RD的最大值,
τDmin为所述动力电池的RC网络的时间常数τD的最小值,τDmax为所述动力电池的RC网络的时间常数τD的最大值。
8.根据权利要求7所述的基于数据模型融合的动力电池容量和荷电状态的估计方法,其特征在于,采用安时积分法根据所述动力电池的最大可用容量Qmax和荷电状态初始值z0估计出所述动力电池的荷电状态z。
CN201610069076.9A 2016-02-01 2016-02-01 基于数据模型融合的动力电池容量和荷电状态的估计方法 Active CN105842627B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610069076.9A CN105842627B (zh) 2016-02-01 2016-02-01 基于数据模型融合的动力电池容量和荷电状态的估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610069076.9A CN105842627B (zh) 2016-02-01 2016-02-01 基于数据模型融合的动力电池容量和荷电状态的估计方法

Publications (2)

Publication Number Publication Date
CN105842627A true CN105842627A (zh) 2016-08-10
CN105842627B CN105842627B (zh) 2018-06-01

Family

ID=56586812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610069076.9A Active CN105842627B (zh) 2016-02-01 2016-02-01 基于数据模型融合的动力电池容量和荷电状态的估计方法

Country Status (1)

Country Link
CN (1) CN105842627B (zh)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106291381A (zh) * 2016-08-16 2017-01-04 北京理工大学 一种联合估计动力电池系统荷电状态与健康状态的方法
CN106324521A (zh) * 2016-09-05 2017-01-11 北京理工大学 一种联合估计动力电池系统参数与荷电状态的方法
CN106772067A (zh) * 2016-11-30 2017-05-31 北京理工大学 多时间尺度iapf滤波估计动力电池荷电状态与健康状态的方法
CN106842045A (zh) * 2017-01-20 2017-06-13 北京理工大学 一种基于自适应权重方法的电池多模型融合建模方法和电池管理系统
CN107037374A (zh) * 2017-05-04 2017-08-11 广东电网有限责任公司电力科学研究院 一种锂离子电池的soc和soh联合估计方法
EP3306332A1 (en) * 2016-10-09 2018-04-11 Optimum Battery Co., Ltd. System and method for estimating parameters of a battery model using partical swarm optimisation
CN108061858A (zh) * 2017-10-30 2018-05-22 东华大学 基于欧姆内阻的锂电池soc评估方法
CN108414937A (zh) * 2017-12-08 2018-08-17 国网北京市电力公司 充电电池荷电状态确定方法及装置
CN108736083A (zh) * 2017-04-25 2018-11-02 通用汽车环球科技运作有限责任公司 用于评估含有表现出电压滞后的材料的电池的方法和设备
CN109239605A (zh) * 2018-11-01 2019-01-18 西南交通大学 一种磷酸铁锂动力电池soc估计方法
CN109425831A (zh) * 2017-08-21 2019-03-05 加百裕工业股份有限公司 估测电池Qmax的方法
CN109950661A (zh) * 2019-03-25 2019-06-28 清华大学 一种动力电池组的内外部同时加热的装置及方法
CN110208704A (zh) * 2019-04-29 2019-09-06 北京航空航天大学 一种基于电压滞后效应的锂电池建模方法和系统
CN110226097A (zh) * 2017-03-03 2019-09-10 康奈可关精株式会社 充电率推定装置以及充电率推定方法
CN110333448A (zh) * 2018-03-30 2019-10-15 比亚迪股份有限公司 电动汽车及动力电池的能量状态soe计算方法、装置
CN110632521A (zh) * 2019-10-23 2019-12-31 北京理工大学 一种锂离子电池容量的融合估计方法
CN110687462A (zh) * 2019-11-04 2020-01-14 北京理工大学 一种动力电池soc与容量全生命周期联合估计方法
CN112526883A (zh) * 2020-11-30 2021-03-19 北京航空航天大学 一种基于智能网联信息的整车能量管理方法
CN112904211A (zh) * 2021-01-04 2021-06-04 北京空间飞行器总体设计部 一种深空探测用锂离子电池组剩余容量估计方法
CN113391216A (zh) * 2021-06-29 2021-09-14 珠海创芯科技有限公司 基于响应曲面设计的铅酸电池可用容量估算方法和装置
CN115494401A (zh) * 2022-11-14 2022-12-20 湖北工业大学 一种基于信息融合的动力电池云端数据清洗方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608542A (zh) * 2012-04-10 2012-07-25 吉林大学 动力电池荷电状态估计方法
CN104535932A (zh) * 2014-12-20 2015-04-22 吉林大学 一种锂离子电池荷电状态估计方法
CN104569835A (zh) * 2014-12-16 2015-04-29 北京理工大学 一种估计电动汽车的动力电池的荷电状态的方法
CN105259511A (zh) * 2015-11-02 2016-01-20 国网黑龙江省电力有限公司伊春供电公司 基于蓄电池运行状态归算的荷电状态估算方法
US20160018468A1 (en) * 2014-07-21 2016-01-21 Richtek Technology Corporation Method of estimating the state of charge of a battery and system thereof
CN105277898A (zh) * 2015-10-27 2016-01-27 浙江大学 一种电池荷电状态的检测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608542A (zh) * 2012-04-10 2012-07-25 吉林大学 动力电池荷电状态估计方法
US20160018468A1 (en) * 2014-07-21 2016-01-21 Richtek Technology Corporation Method of estimating the state of charge of a battery and system thereof
CN104569835A (zh) * 2014-12-16 2015-04-29 北京理工大学 一种估计电动汽车的动力电池的荷电状态的方法
CN104535932A (zh) * 2014-12-20 2015-04-22 吉林大学 一种锂离子电池荷电状态估计方法
CN105277898A (zh) * 2015-10-27 2016-01-27 浙江大学 一种电池荷电状态的检测方法
CN105259511A (zh) * 2015-11-02 2016-01-20 国网黑龙江省电力有限公司伊春供电公司 基于蓄电池运行状态归算的荷电状态估算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
AMIR VASEBI ET AL: "A novel combined battery model for state-of-charge estimation in lead-acid batteries based on extended Kalman filter for hybrid electric vehicle applications", 《JOURNAL OF POWER SOURCES》 *
熊瑞: "基于数据模型融合的电动车辆动力电池组状态估计研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106291381A (zh) * 2016-08-16 2017-01-04 北京理工大学 一种联合估计动力电池系统荷电状态与健康状态的方法
CN106291381B (zh) * 2016-08-16 2018-09-11 北京理工大学 一种联合估计动力电池系统荷电状态与健康状态的方法
CN106324521A (zh) * 2016-09-05 2017-01-11 北京理工大学 一种联合估计动力电池系统参数与荷电状态的方法
CN106324521B (zh) * 2016-09-05 2018-09-11 北京理工大学 一种联合估计动力电池系统参数与荷电状态的方法
EP3306332A1 (en) * 2016-10-09 2018-04-11 Optimum Battery Co., Ltd. System and method for estimating parameters of a battery model using partical swarm optimisation
CN106772067B (zh) * 2016-11-30 2018-12-28 北京理工大学 多时间尺度估计动力电池荷电状态与健康状态的方法
CN106772067A (zh) * 2016-11-30 2017-05-31 北京理工大学 多时间尺度iapf滤波估计动力电池荷电状态与健康状态的方法
CN106842045A (zh) * 2017-01-20 2017-06-13 北京理工大学 一种基于自适应权重方法的电池多模型融合建模方法和电池管理系统
CN106842045B (zh) * 2017-01-20 2023-04-28 北京理工大学 一种基于自适应权重方法的电池多模型融合建模方法和电池管理系统
CN110226097B (zh) * 2017-03-03 2022-07-05 康奈可关精株式会社 观测器增益的设定方法
CN110226097A (zh) * 2017-03-03 2019-09-10 康奈可关精株式会社 充电率推定装置以及充电率推定方法
CN108736083B (zh) * 2017-04-25 2021-06-15 通用汽车环球科技运作有限责任公司 用于评估含有表现出电压滞后的材料的电池的方法和设备
US10901042B2 (en) 2017-04-25 2021-01-26 GM Global Technology Operations LLC Method and apparatus for evaluating battery cells containing materials that exhibit voltage hysteresis
CN108736083A (zh) * 2017-04-25 2018-11-02 通用汽车环球科技运作有限责任公司 用于评估含有表现出电压滞后的材料的电池的方法和设备
CN107037374A (zh) * 2017-05-04 2017-08-11 广东电网有限责任公司电力科学研究院 一种锂离子电池的soc和soh联合估计方法
CN107037374B (zh) * 2017-05-04 2019-10-08 广东电网有限责任公司电力科学研究院 一种锂离子电池的soc和soh联合估计方法
CN109425831A (zh) * 2017-08-21 2019-03-05 加百裕工业股份有限公司 估测电池Qmax的方法
CN108061858A (zh) * 2017-10-30 2018-05-22 东华大学 基于欧姆内阻的锂电池soc评估方法
CN108414937A (zh) * 2017-12-08 2018-08-17 国网北京市电力公司 充电电池荷电状态确定方法及装置
CN110333448B (zh) * 2018-03-30 2021-02-23 比亚迪股份有限公司 电动汽车及动力电池的能量状态soe计算方法、装置
CN110333448A (zh) * 2018-03-30 2019-10-15 比亚迪股份有限公司 电动汽车及动力电池的能量状态soe计算方法、装置
CN109239605A (zh) * 2018-11-01 2019-01-18 西南交通大学 一种磷酸铁锂动力电池soc估计方法
CN109239605B (zh) * 2018-11-01 2019-09-27 西南交通大学 一种磷酸铁锂动力电池soc估计方法
CN109950661A (zh) * 2019-03-25 2019-06-28 清华大学 一种动力电池组的内外部同时加热的装置及方法
CN109950661B (zh) * 2019-03-25 2022-04-12 清华大学 一种动力电池组的内外部同时加热的装置及方法
CN110208704A (zh) * 2019-04-29 2019-09-06 北京航空航天大学 一种基于电压滞后效应的锂电池建模方法和系统
CN110632521A (zh) * 2019-10-23 2019-12-31 北京理工大学 一种锂离子电池容量的融合估计方法
CN110687462A (zh) * 2019-11-04 2020-01-14 北京理工大学 一种动力电池soc与容量全生命周期联合估计方法
CN110687462B (zh) * 2019-11-04 2020-09-04 北京理工大学 一种动力电池soc与容量全生命周期联合估计方法
CN112526883A (zh) * 2020-11-30 2021-03-19 北京航空航天大学 一种基于智能网联信息的整车能量管理方法
CN112526883B (zh) * 2020-11-30 2021-09-24 北京航空航天大学 一种基于智能网联信息的整车能量管理方法
CN112904211B (zh) * 2021-01-04 2022-10-28 北京空间飞行器总体设计部 一种深空探测用锂离子电池组剩余容量估计方法
CN112904211A (zh) * 2021-01-04 2021-06-04 北京空间飞行器总体设计部 一种深空探测用锂离子电池组剩余容量估计方法
CN113391216A (zh) * 2021-06-29 2021-09-14 珠海创芯科技有限公司 基于响应曲面设计的铅酸电池可用容量估算方法和装置
CN115494401A (zh) * 2022-11-14 2022-12-20 湖北工业大学 一种基于信息融合的动力电池云端数据清洗方法
CN115494401B (zh) * 2022-11-14 2023-03-10 湖北工业大学 一种基于信息融合的动力电池云端数据清洗方法

Also Published As

Publication number Publication date
CN105842627B (zh) 2018-06-01

Similar Documents

Publication Publication Date Title
CN105842627A (zh) 基于数据模型融合的动力电池容量和荷电状态的估计方法
Wang et al. State of charge estimation for LiFePO4 battery via dual extended kalman filter and charging voltage curve
CN106443474B (zh) 一种动力电池系统寿命衰退特征快速识别的方法
CN107991623B (zh) 一种考虑温度和老化程度的电池安时积分soc估计方法
CN104569835B (zh) 一种估计电动汽车的动力电池的荷电状态的方法
CN103344917B (zh) 一种锂电池循环寿命快速测试方法
Hu et al. Online estimation of an electric vehicle lithium-ion battery using recursive least squares with forgetting
CN106716158B (zh) 电池荷电状态估算方法和装置
CN108663620B (zh) 一种动力电池组荷电状态估算方法及系统
CN103020445B (zh) 一种电动车车载磷酸铁锂电池的soc与soh预测方法
CN105425153B (zh) 一种估计电动车辆的动力电池的荷电状态的方法
CN106291378B (zh) 一种电动汽车动力电池soh的测算方法
Coleman et al. State-of-charge determination from EMF voltage estimation: Using impedance, terminal voltage, and current for lead-acid and lithium-ion batteries
CN104977537B (zh) 电池soc的确定方法及使用该方法的电池管理系统
Liu et al. A comparative study of equivalent circuit models and enhanced equivalent circuit models of lithium-ion batteries with different model structures
CN104535932A (zh) 一种锂离子电池荷电状态估计方法
CN109870659A (zh) 应用滑窗寻优策略的锂离子电池健康状态估算方法
CN105425154B (zh) 一种估计电动汽车的动力电池组的荷电状态的方法
CN104678316A (zh) 锂离子电池荷电状态估算方法和装置
CN108732508B (zh) 一种锂离子电池容量的实时估计方法
CN102565710A (zh) 用于估计蓄电池健康状态的方法和装置
US9555718B2 (en) Estimation and compensation of battery measurement and asynchronization biases
CN102468521A (zh) 用于估计蓄电池健康状态的方法和装置
CN102937704A (zh) 一种动力电池rc等效模型的辨识方法
CN104122447A (zh) 一种电动汽车动力电池组直流阻抗的在线估算方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant