CN115598540A - 一种考虑宽温度的锂电池参数辨识及状态联合估计方法 - Google Patents

一种考虑宽温度的锂电池参数辨识及状态联合估计方法 Download PDF

Info

Publication number
CN115598540A
CN115598540A CN202211309229.4A CN202211309229A CN115598540A CN 115598540 A CN115598540 A CN 115598540A CN 202211309229 A CN202211309229 A CN 202211309229A CN 115598540 A CN115598540 A CN 115598540A
Authority
CN
China
Prior art keywords
battery
state
soc
temperature
parameter
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.)
Pending
Application number
CN202211309229.4A
Other languages
English (en)
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.)
Yancheng Institute of Technology
Harbin Institute of Technology Weihai
Original Assignee
Yancheng Institute of Technology
Harbin Institute of Technology Weihai
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 Yancheng Institute of Technology, Harbin Institute of Technology Weihai filed Critical Yancheng Institute of Technology
Priority to CN202211309229.4A priority Critical patent/CN115598540A/zh
Publication of CN115598540A publication Critical patent/CN115598540A/zh
Pending legal-status Critical Current

Links

Images

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/367Software therefor, e.g. for battery testing using modelling or look-up tables
    • 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/396Acquisition or processing of data for testing or for monitoring individual cells or groups of cells within a battery

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Tests Of Electric Status Of Batteries (AREA)

Abstract

本发明公布了一种考虑宽温度的锂电池参数辨识及状态联合估计方法,所述方法为:由电池充放电实验得到电池的开路电压和SOC,拟合得到UOCV‑SOC曲线,从而建立电池一阶RC等效电路模型,获得空间状态方程;依次通过基于改进递推最小二乘法的参数辨识模块确定T‑α关系式、基于改进无迹卡尔曼滤波的状态估计模块确定T‑μ关系式;获取电池实时温度,经温度判别器决定给定温度是否变化,再根据T‑α、T‑μ关系式分别得到参数辨识最佳遗忘因子α、状态估计最佳遗忘因子μ,将α反馈给模型参数辨识模块,以更新模型参数,更新的模型参数再与μ一起反馈给状态估计模块,以进一步更新电池状态估计值,即得到更准确的SOC值。

Description

一种考虑宽温度的锂电池参数辨识及状态联合估计方法
技术领域
本发明涉及动力电池管理领域,特别是涉及动力电池系统参数与荷电状态高精度估计。
背景技术
当今的汽车正在向电动化、智能化、网联化、共享化发展,电动汽车无疑是研究的热点领域,电动汽车的安全性也备受人们关注。而电动汽车的安全问题的重要一环是动力电池系统。在动力电池系统中,电池的荷电状态(SOC)估计是其中的一个重要问题。本发明主要针对的是电池的参数辨识和荷电状态SOC的联合估计,其中电池模型参数对应电池内阻R0、极化内阻Rp、极化电容Cp和开路电压UOCV
目前,对SOC估计常采用基于锂电池模型的滤波类方法,该方法的核心思路是在现有的电池模型的基础上,结合当前状态控制理论技术,更高精度、更广泛地估计SOC。在估计方法上,卡尔曼滤波(Kalman filtering)算法通过构建线性系统状态方程,对比系统输入输出数据,对系统进行最优估计,是一种目前广泛应用的滤波方法。但现有技术大多是在某一个动态工况、温度等条件下仿真和实验验证,对于动力电池全工作温度环境下的SOC估计结果往往偏离客观情况,这是由于不同温度和老化等工作条件下,电池的参数变化很大,而电池模型的参数和状态具有耦合性的缘故。因此,有研究学者提出了多种模型参数与SOC的联合估计方法。其中参数估计常采用递推最小二乘法,SOC仍采用卡尔曼类滤波器的方法。为了解决数据饱和问题,有学者分别在参数辨识和状态估计过程分别引入遗忘因子α和μ,即带有遗忘因子的递推最小二乘法(Recursive Least Square with Forgetting Factor,FFRLS)和带衰减记忆的卡尔曼滤波类算法(或称为带遗忘因子的卡尔曼滤波类算法)。
在其他领域,对于FFRLS,已有人提出建立关于变遗忘因子公式实时更新α,但锂电池易受温度影响,而以往建立的遗忘因子递推经验公式未考虑电池温度T变化因素,即使实时更新递推遗忘因子α,得到的α也并非是相应温度T下的最优或次优值,且实时递推计算量较大。对于带衰减记忆卡尔曼滤波器中的衰减因子(也可称为遗忘因子)μ通常通过经验设定一恒定值而应用于各温度T下,也并非是各温度T下的最优或次优值。因此,目前尚缺少一种适用于宽温度下锂电池参数辨识和状态联合估算过程的遗忘因子高效确定方法。
发明内容
本发明解决的问题是在于提供一种考虑宽温度的锂电池参数及状态联合高精度估计方法。电池参数辨识阶段使用FFRLS,即采用基于改进递推最小二乘法的参数辨识模块来实现电池模型参数辨识,状态估计阶段使用带遗忘因子的无迹卡尔曼滤波UKF,即采用基于改进的无迹卡尔曼滤波状态估计模块来实现电池荷电状态估计。
为实现上述目的,本发明提供如下技术方案:
S1、通过锂电池充放电实验获得电池实验数据,包括电池开路电压UOCV、端电压UL、荷电状态SOC、电流I、温度T,由电池UOCV、SOC实验数据经拟合得到UOCV-SOC关系式,建立电池一阶RC等效电路模型,获得电池空间状态方程;
S2、根据电池充放实验获得的电池实验数据,结合电池空间状态方程,经基于改进递推最小二乘法的参数辨识模块确定温度-参数遗忘因子(T-α)关系式;
S3、根据由基于改进递推最小二乘法的参数辨识模块得到的电池模型参数(包括电池内阻R0、极化内阻Rp、极化电容Cp和开路电压UOCV),以及根据电池充放实验获得的电池实验数据,再结合电池空间状态方程,经基于改进无迹卡尔曼滤波的状态估计模块确定温度-状态遗忘因子(T-μ)关系式;
S4、获取电池实时温度,经温度判别器决定给定温度是否变化,再根据T-α、T-μ关系式分别得到最佳参数遗忘因子α、最佳状态遗忘因子μ,将α反馈给参数辨识模块,以更新电池模型参数,更新的电池模型参数再与μ一起反馈给状态估计模块,以进一步更新电池状态估计值,即得到更准确的SOC预测值。
在步骤S1中,在锂电池常用工作温度范围内(-10℃~50℃),每次以10℃为间隔分别进行:以电池厂家推荐的充放电倍率进行的电池容量实验、开路电压实验、动态工况实验(UDDS工况),在动态工况下测量电池电流i、电压UP和温度T;结合容量实验和开路电压实验,获得每个采样点的电池UOCV和SOC,采用下式拟合得到UOCV-SOC关系式:
UOCV=k0+k1SOC+k2SOC2+k3SOC3+k4log(SOC)+k5log(1-SOC)
式中,k0~k5为待确定的参数。
在步骤S1中,根据基尔霍夫定律得到锂电池一阶RC等效电路模型为
UL=UOCV-IR0-UP
Figure BDA0003907284500000021
式中,UL为电池端电压,UOCV为电池开路电压,UP为极化电压,I是电池电流,R0为电池内阻,Rp、Cp分别为极化内阻和极化电容;
电池模型经过离散化处理后得到电池系统的状态空间方程为
Figure BDA0003907284500000031
[UL,k]=UOCV(SOCk)-Up,k-R0Ik+vk
式中,下标k表示第k时刻的数据,Ts为实验采集数据时间,wk与vk是均值为零、协方差为Q和R的高斯白噪声,η表示库伦效率,QN表示电池容量。
在步骤S2中,产生T-a关系式步骤如下:
1)对电池模型进行线性化处理:
UL,k=UOCV,k-a1UOCV,k-1+a1UL,k-1+a2Ik+a3Ik-1
式中,a1、a2、a3是与电池模型参数有关的系数,具体如下式:
Figure BDA0003907284500000032
由式(1)得到相应的参数矩阵θ和数据矩阵
Figure BDA0003907284500000033
Figure BDA0003907284500000034
2)初始化UOCV、R0、Rp和Cp四个需要辨识的电池模型参数以及协方差矩阵P,并通过式(2)转换为mk、a1、a2、a3,即得到参数矩阵初值
Figure BDA0003907284500000035
和协方差矩阵P0,同时给定参数遗忘因子α的初始值α0
3)代入动态工况的电池电流和电压数据,更新电池端电压测量值yk
4)参数及误差协方差矩阵更新:
Figure BDA0003907284500000036
Figure BDA0003907284500000037
Figure BDA0003907284500000038
式中,
Figure BDA0003907284500000039
为上一时刻参数矩阵预测值,
Figure BDA00039072845000000310
为本时刻电池端电压的预测值,kk为增益系数,Pk为协方差矩阵;
5)将参数矩阵θk逆变换得到需要的电池模型参数:
Figure BDA0003907284500000041
6)由辨识出的电池模型参数经式(1)计算电池端电压预测值
Figure BDA0003907284500000042
并与测量值yk比较计算出平均绝对误差(MAE)与均方根误差(RMSE):
Figure BDA0003907284500000043
Figure BDA0003907284500000044
其中
Figure BDA0003907284500000045
用以评价辨识结果的好坏,N为自然数;
7)确定当前温度下最佳参数遗忘因子:由于α是一个大于1且不易过大的数值,首先令α在1.00~1.10之间按等间距0.001来逐一取值,将每个α取值代入步骤4)~步骤6),然后根据MAE和RMSE来找出辨识效果最好的参数遗忘因子α,作为当前温度Tk下的最佳参数遗忘因子αk
8)使用不同温度Tj下的动态工况数据,通过对初始参数及其对应的协方差矩阵进行更新。重复上述步骤3)~步骤7),同样得到该温度Tj下对应的最佳参数遗忘因子αj,j为大于的自然数,最后使用七次多项式拟合得到T-α关系式:
α=k0+k1T+k2T2+k3T2+k5T5+k6T6+k7T7
其中k0~k7是拟合系数。
在步骤S3中,产生T-μ关系式步骤如下:
1)设定SOC初始值为1,状态变量数为2个,状态变量x=[SOC,Up],Up为极化电压,观测变量设定为电池模型参数R0、Rp、Cp和UOCV
2)对SOC、Up以及相应的过程噪声方差Q、测量噪声方差R进行初始化;
3)将SOC与Up进行sigma化,并计算其权重ω:
①均值和协方差初始化,具体公式如下:
Figure BDA0003907284500000046
Figure BDA0003907284500000047
Figure BDA0003907284500000051
是状态变量均值的初值,P0是协方差的初值;
②产生2n+1个sigma点Xi,计算相应的权重ω:
Figure BDA0003907284500000052
式中,λ为尺度调节系数,修改其大小用来降低预测误差,采样点xi与x具有近似的高斯分布,其中
Figure BDA0003907284500000053
表示矩阵
Figure BDA0003907284500000054
的第i列;
Figure BDA0003907284500000055
式中,n表示状态变量维度,λ=ε2(n+m)-n,ε为常数,取值范围为10-4≤ε<1;β是状态分布参数;ωm、ωc分别为方差及均值的权重;
4)计算卡尔曼增益:
①计算这些sigma点经过非线性变换的结果,具体公式如下:
Figure BDA0003907284500000056
Figure BDA0003907284500000057
式中,
Figure BDA0003907284500000058
是由k-1时刻状态变量经非线性状态方程f(·)后预测得到的k时刻状态变量,
Figure BDA0003907284500000059
是更新后的状态变量预测值;
②状态变量协方差更新,具体公式如下:
Figure BDA00039072845000000510
式中,Qk为k时刻的过程噪声方差;
③观测变量更新,具体公式如下:
Figure BDA00039072845000000511
Figure BDA00039072845000000512
式中,
Figure BDA0003907284500000061
是状态变量经测量方程h(·)后预测得到新的观测变量预测值,
Figure BDA0003907284500000062
是更新后的观测变量预测值;
④误差协方差更新,在观测噪声协方差矩阵中引入了遗忘因子μ,具体公式如下:
Figure BDA0003907284500000063
Rk+1=μ×Rk
Figure BDA0003907284500000064
式中Rk是k时刻的观测噪声方差,μ为状态遗忘因子;
⑤卡尔曼增益更新,具体公式如下:
Figure BDA0003907284500000065
式中,kk为卡尔曼增益,表征模型预测误差与量测误差的比重;
5)根据卡尔曼增益计算出系统状态和协方差更新:
Figure BDA0003907284500000066
Px,k/k=Px,k/k-1-kkPyy,kkk T
式中,
Figure BDA0003907284500000067
是k时刻的估计值,Px,k/k是其误差协方差;
6)由步骤5)中得到电池SOC预测值,再与荷电状态参考值SOCr比较,分别得到SOC平均绝对误差MAE及SOC均方根误差RMSE,
Figure BDA0003907284500000068
Figure BDA0003907284500000069
其中,SOCerrori=SOCg-SOCr,SOCg为SOC预测值;
7)确定当前温度下最佳状态遗忘因子:首先令状态遗忘因子μ在1.00~1.30以0.01等间距逐一取值,将每个μ取值代入步骤4)~步骤6),然后根据SOCMAE和SOCRMSE来找出辨识效果最好的状态遗忘因子α,作为当前温度Tk下的最佳参数遗忘因子μk
8)重复使用-10℃~50℃的数据,重复上述步骤4)~步骤7),同样得到该温度Tj下对应的最佳参数遗忘因子μj,使用七次多项式拟合得到T-μ关系式:
μ=k0+k1T+k2T2+k3T3+k5T5+k6T6+k7T7
式中,k0~k7是拟合系数。
在步骤S4中,所述的温度判别器设计如下:令温度差ΔT=Tk-Tk-1,Tk为本时刻电池实时温度,Tk-1为上一时刻电池实时温度,若ΔT大于0.5℃,则在Tk温度下进行遗忘因子的更新,否则在Tk-1温度下进行遗忘因子的更新。
本发明提出的考虑宽温度影响的锂电池参数辨识和状态联合估计方法,与传统方法相比具有以下优势:采用代有随温度变化的遗忘因子算法进行电池模型参数辨识和状态估计,能够适应电池工作环境温度的变化,提高全工作温度环境下SOC估计精度。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用来解释本发明,并不构成对本发明的限制。
在附图中:
图1为本发明一阶RC等效电路模型图;
图2为本发明方法流程图;
图3为本发明25℃、DST工况在线参数辨识结果图;
图4为本发明不同遗忘因子下SOC估计及误差分析图。
具体实施方式
下面结合附图对本发明所提供的考虑宽温度的锂电池参数辨识和状态联合估计方法进行详尽的阐述。
在本发明的一个优选实施例中,基于一款未经使用过的锂电池作为研究对象,所用电池型号为4.1V/2Ah,动态工况采用DST工况。认为其初始容量为100%。
1.建模及实验
①为了改善内阻模型的不足之处,考虑到锂电池的极化作用,在原有电路基础上并联上一个RC模块,构成了一阶RC等效电路模型。相较于内阻模型而言,一阶RC等效电路模型能更好的展现电池内部动态变化,其具体结构如附图1所示。其中UOCV表示电池开路电压,这里不把它看作恒定不变的值,R0表示电池的欧姆内阻,Rp和Cp分别表示电池的极化电阻和极化电容。欧姆内阻R0能够在短时间内变化,从而体现电压随充放电时的变化,在加入了RC电路以后,RC模块能够表示电池充电或者放电结束以后电池内部电压的逐渐下降或者上升。一阶RC等效电路模型在结构上同样简单,但是却拥有很高的精度,因此应用也较为广泛。
②一阶RC等效电路模型需要辨识欧姆内阻R0、极化内阻Rp、极化电容Cp和开路电压UOCV。根据基尔霍夫定律,可以得到如下关系式:
Figure BDA0003907284500000081
Figure BDA0003907284500000082
结合式(1)和式(2),经过离散化处理,得到电池系统的状态空间方程如式(3)和式(4):
Figure BDA0003907284500000083
[UL,k]=UOCV(SOCk)-Up,k-R0Ik+vk (4)
式中,下标k表示第k时刻的数据,Ts为实验采集数据时间,wk与vk是均值为零、协方差为Q和R的高斯白噪声,η表示库伦效率,QN表示电池容量。
③在上述三个温度下进行:电池容量实验、OCV实验、动态工况实验DST工况。在本发明中采用单向脉冲加间隔静置法来获取开路电压。首先用单向脉冲对电池充至电压上限,全程保持温度不变,使电池满电,然后静置1小时,使电池处于电位平衡状态;接着以1A的电流给电池放电12分钟,此时认为电池电路中的SOC达到初始值的90%,之后再静置1小时。重复以上步骤即可获得SOC从100%到0%的开路电压值,25℃的开路电压实验结果如附图3所示。
结合容量实验和OCV实验获得每个采样点的开路电压和SOC,使用下式(5)拟合得到UOCV-SOC关系式。
UOCV=k0+k1SOC+k2SOC2+k3SOC3+k4 log(SOC)+k5 log(1-SOC) (5)
2.由基于改进递推最小二乘法的参数辨识模块确定温度-参数遗忘因子(T-α)关系式
①将上述状态空间方程进行转换得到
UL,k=UOCV,k-a1UOCV,k-1+a1UL,k-1+a2Ik+a3Ik-1 (6)
其中的a1,a2,a3具体如式(7):
Figure BDA0003907284500000084
相应的参数矩阵θ和数据矩阵
Figure BDA0003907284500000085
为:
Figure BDA0003907284500000086
式中的mk采用开路电压的递推算法求得。由上面的公式,把原本需要辨识的四个参数转换为mk、a1、a2、a3。在每次辨识出参数后,需要通过计算转换成需要的参数。具体转化公式(9)如下:
Figure BDA0003907284500000091
②初始化UOCV、R0、Rp和Cp四个要辨识的参数以及协方差矩阵,同时给出遗忘因子的初始值。以25℃情况为例,
Figure BDA0003907284500000092
α=0.95。
③获取离线的动态工况的电流和电压数据,通过动态工况实验得到的电流及电压数据用如下递推最小二乘法公式重复进行参数辨识,一直递推到没有新的数据。
Figure BDA0003907284500000093
Figure BDA0003907284500000094
Figure BDA0003907284500000095
得到辨识参数随时间的变化,如附图3所示,用辨识参数通过公式(5)计算端电压预测值,与测量值比较计算平均绝对误差与均方根误差,用以评价辨识结果的好坏。
Figure BDA0003907284500000096
Figure BDA0003907284500000097
④在1.00~1.10之间以0.001的间距改变遗忘因子α,找出辨识效果最好的遗忘因子α作为该温度下的最佳遗忘因子。在标准25℃工况下,通过调试参数,最终调节遗忘因子α=1.031时,有较好的参数辨识效果。图3为本发明25℃、DST工况在线辨识结果图
⑤使用不同温度下的动态工况数据,通过对初始参数及对应的协方差矩阵进行调节重复上述实验过程。导入10℃、DST工况的电流及电压数据,最终调节至遗忘因子α=1.053时,有较好的参数辨识效果,同理,导入40℃、DST工况的电流及电压数据,最终调节至遗忘因子α=1.004时,有较好的参数辨识效果。由此我们得出结论:在运用带有遗忘因子α的递推最小二乘法进行参数辨识时,遗忘因子α变化较小,并随着温度的升高应适当降低α。使用多项式拟合得到T-α关系式:
α=-0.00001T2-0.00108T+1.0648 (15)
3.由基于改进无迹卡尔曼滤波的状态估计模块确定温度-状态遗忘因子(T-μ)关系式
①设定SOC初始值为1,状态变量数为2个,状态变量x=[SOC,Up],观测变量设定为电池模型参数R0、Rp、Cp和UOCV
②进行SOC、Up,以及相应的噪声Q、R的初始化。在本例中,设定初始化的噪声R=0.0001,Up=0。初始化
Figure BDA0003907284500000101
③接着按照无迹变换计算步骤,首先将SOC与Up进行sigma化并计算其权重,同时计算出相应的卡尔曼增益,并更新状态空间方程以及协方差矩阵。
S1.均值和协方差初始化。
Figure BDA0003907284500000102
Figure BDA0003907284500000103
S2.产生2n+1个sigma点,计算相应的权重。
Figure BDA0003907284500000104
Figure BDA0003907284500000105
其中,n=2表示状态变量维度,ε=0.5为离散程度系数,ε在实际选取中应尽量小一些,越小则越接近于总样本水平;β=2表示验前分布系数,当处于高斯分布时,取β=2时最优;m=2是保证m+n≠0的辅助尺度因子;λ是尺度调节系数:
λ=ε2(n+m)-n (20)
S3.计算这些sigma点经过非线性变换的结果,具体公式如下:
Figure BDA0003907284500000106
Figure BDA0003907284500000107
S4.状态变量协方差更新具体公式如下:
Figure BDA0003907284500000108
S5.观测变量更新,具体公式如下:
Figure BDA0003907284500000109
Figure BDA0003907284500000111
S6.误差协方差更新,在观测噪声协方差矩阵中引入了遗忘因子μ,具体公式如下:
Figure BDA0003907284500000112
Rk+1=μ×Rk (27)
Figure BDA0003907284500000113
S7.卡尔曼增益更新,具体公式如下:
Figure BDA0003907284500000114
S8.状态更新和协方差更新,具体公式如下:
Figure BDA0003907284500000115
Figure BDA0003907284500000116
执行完毕后将估计出的端电压值与实际端电压值进行对比,得到的误差对误差矩阵进行更新,一直重复循环直到锂电池放电结束。
④选取25℃、DST工况作为标准工况,由经验可得该遗忘因子μ是一个比1大的值,但不会过大,从μ=1.00~1.30以0.01等间距取值估计锂电池SOC,当遗忘因子μ取值较大时,SOC在后半段估计时有较大的误差,当μ=1.08时,整体SOC估计误差控制在2%内,因有较好的估计效果,因此,把μ=1.08作为25度时SOC状态估计最终取值。图4为本发明不同遗忘因子下状态估计及误差分析图。
⑤相应的改变DST况的温度,获得不同温度下的SOC状态估计效果及误差分析图以探究温度对遗忘因子μ的影响。在改变温度时应相应改变初始状态协方差矩阵及过程噪声、观测噪声参数值,最终调节10度时μ=1.11,40度时μ=1.13时有较好的状态估计效果,使用多项式拟合得到T-μ关系式:
μ=1.1744T2-0.0082T+1.1744 (32)
4.在线实验调用实现精确SOC估计
在给定的温度下,如15℃,通过拟合关系式计算出该温度下最佳的遗忘因子α=1.046,μ=1.091确定遗忘因子后重复上述参数辨识和状态估计算法实现更加精确的SOC估计效果。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。

Claims (6)

1.一种考虑宽温度的锂电池参数辨识及状态联合估计方法,所述方法步骤如下:
S1、通过锂电池充放电实验获得电池实验数据,包括电池开路电压UOCV、端电压UL、荷电状态SOC、电流I、温度T,由电池UOCV、SOC实验数据经拟合得到UOCV-SOC关系式,建立电池一阶RC等效电路模型,获得电池空间状态方程;
S2、根据电池充放实验获得的电池实验数据,结合电池空间状态方程,经基于改进递推最小二乘法的参数辨识模块确定温度-参数遗忘因子(T-α)关系式;
S3、根据由基于改进递推最小二乘法的参数辨识模块得到的电池模型参数(包括电池内阻R0、极化内阻Rp、极化电容Cp和开路电压UOCV),以及根据电池充放实验获得的电池实验数据,再结合电池空间状态方程,经基于改进无迹卡尔曼滤波的状态估计模块确定温度-状态遗忘因子(T-μ)关系式;
S4、获取电池实时温度,经温度判别器决定给定温度是否变化,再根据T-α、T-μ关系式分别得到最佳参数遗忘因子α、最佳状态遗忘因子μ,将α反馈给参数辨识模块,以更新电池模型参数,更新的电池模型参数再与μ一起反馈给状态估计模块,以进一步更新电池状态估计值,即得到更准确的SOC预测值。
2.根据权利要求1所述的一种考虑宽温度的锂电池参数辨识及状态联合估计方法,其特征在于,所述的T-α关系式产生步骤如下:
1)对电池模型进行线性化处理得:
UL,k=UOCV,k-a1UOCV,k-1+a1UL,k-1+a2Ik+a3Ik-1 (1)
式中,a1、a2、a3是与电池模型参数有关的系数,具体计算如下:
Figure FDA0003907284490000011
由式(1)得到相应的参数矩阵θ和数据矩阵
Figure FDA0003907284490000012
Figure FDA0003907284490000013
2)初始化UOCV、R0、Rp和Cp四个需要辨识的电池模型参数以及协方差矩阵P,并通过式(2)转换为mk、a1、a2、a3,即得到参数矩阵初值
Figure FDA0003907284490000014
和协方差矩阵P0,同时给定参数遗忘因子α的初始值α0
3)代入动态工况的电池电流和电压数据,更新电池端电压测量值yk
4)参数及误差协方差矩阵更新:
Figure FDA0003907284490000015
Figure FDA0003907284490000016
Figure FDA0003907284490000021
式中,
Figure FDA0003907284490000022
为上一时刻参数矩阵预测值,
Figure FDA0003907284490000023
为本时刻电池端电压的预测值,kk为增益系数,Pk为协方差矩阵;
5)将参数矩阵θk逆变换得到需要的电池模型参数:
Figure FDA0003907284490000024
6)由辨识出的电池模型参数经式(1)计算电池端电压预测值
Figure FDA0003907284490000025
并与测量值yk比较计算出平均绝对误差(MAE)与均方根误差(RMSE):
Figure FDA0003907284490000026
Figure FDA0003907284490000027
其中
Figure FDA0003907284490000028
用以评价辨识结果的好坏;
7)确定当前温度下最佳参数遗忘因子:由于α是一个大于1且不易过大的数值,首先令α在1.00~1.10之间按等间距0.001来逐一取值,将每个α取值代入步骤4)~步骤6),然后根据MAE和RMSE来找出辨识效果最好的参数遗忘因子α,作为当前温度Tk下的最佳参数遗忘因子αk
8)使用不同温度Tj下的动态工况数据,通过对初始参数及其对应的协方差矩阵进行更新。重复上述步骤3)~步骤7),同样得到该温度Tj下对应的最佳参数遗忘因子αj,j为大于的自然数,最后使用七次多项式拟合得到T-α关系式:
α=k0+k1T+k2T2+k3T3+k5T5+k6T6+k7T7 (10)
其中k0~k7是拟合系数。
3.根据权利要求1所述的一种考虑宽温度的锂电池参数辨识及状态联合估计方法,
其特征在于,所述的T-μ关系式产生步骤如下:
1)设定SOC初始值为1,状态变量数为2个,状态变量x=[SOC,Up],Up为极化电压,观测变量设定为电池模型参数R0、Rp、Cp和UOCV
2)对SOC、Up以及相应的过程噪声方差Q、测量噪声方差R进行初始化;
3)将SOC与Up进行sigma化,并计算其权重ω:
第一步,均值和协方差初始化:
Figure FDA0003907284490000029
Figure FDA00039072844900000210
式中,
Figure FDA00039072844900000211
是状态变量均值的初值,P0是协方差的初值;
第二步,产生2n+1个sigma点Xi,计算相应的权重ω:
Figure FDA0003907284490000031
式中,λ为尺度调节系数,修改其大小用来降低预测误差,采样点xi与x具有近似的高斯分布,
Figure FDA0003907284490000032
表示矩阵
Figure FDA0003907284490000033
的第i列;
Figure FDA0003907284490000034
式中,n表示状态变量维度,λ=ε2(n+m)-n,ε为常数,取值范围为10-4≤ε<1;β是状态分布参数;ωm、ωc分别为方差及均值的权重;
4)计算卡尔曼增益:
第一步,状态变量更新,具体公式如下:
Figure FDA0003907284490000035
Figure FDA0003907284490000036
Figure FDA0003907284490000037
是由k-1时刻状态变量经非线性状态方程f(·)后预测得到的k时刻状态变量,
Figure FDA0003907284490000038
是更新后的状态变量预测值;
第二步,状态变量协方差Px,k/k-1更新,具体公式如下:
Figure FDA0003907284490000039
式中,Qk为k时刻的过程噪声方差;
第三步,观测变量更新,具体公式如下:
Figure FDA00039072844900000310
Figure FDA00039072844900000311
式中,
Figure FDA00039072844900000312
是状态变量经测量方程h(·)后预测得到新的观测变量预测值,
Figure FDA00039072844900000313
是更新后的观测变量预测值;
第四步,误差协方差更新,在观测噪声协方差矩阵中引入了状态遗忘因子μ,具体公式如下:
Figure FDA00039072844900000314
Rk+1=μ×Rk (21)
Figure FDA0003907284490000041
式中Rk是k时刻的观测噪声方差,
第五步,卡尔曼增益更新,具体公式如下:
Figure FDA0003907284490000042
kk为卡尔曼增益,表征模型预测误差与量测误差的比重;
5)根据卡尔曼增益更新系统状态和协方差:
Figure FDA0003907284490000043
Figure FDA0003907284490000044
式中,
Figure FDA0003907284490000045
是k时刻的估计值,Px,k/k是其误差协方差;
6)由步骤5)中得到电池SOC预测值,再与荷电状态参考值SOCr比较,分别得到SOC平均绝对误差MAE及SOC均方根误差RMSE;
Figure FDA0003907284490000046
Figure FDA0003907284490000047
其中SOCerrori=SOCg-SOCr,SOCg为SOC预测值;
7)确定当前温度下最佳状态遗忘因子:首先令状态遗忘因子μ在1.00~1.30以0.01等间距逐一取值,将每个μ取值代入步骤4)~步骤6),然后根据SOCMAE和SOCRMSE来找出辨识效果最好的状态遗忘因子α,作为当前温度Tk下的最佳参数遗忘因子μk
8)重复使用-10℃~50℃的数据,重复上述步骤4)~步骤7),同样得到该温度Tj下对应的最佳参数遗忘因子μj,使用七次多项式拟合得到T-μ关系式:
μ=k0+k1T+k2T2+k3T3+k5T5+k6T6+k7T7 (28)
式中,k0~k7是拟合系数。
4.根据权利要求1所述的一种考虑宽温度的锂电池参数辨识及状态联合估计方法,其特征在于,所述的温度判别器设计如下:令温度差ΔT=Tk-Tk-1,Tk为本时刻电池实时温度,Tk-1为上一时刻电池实时温度,若ΔT大于0.5℃,则在Tk温度下进行遗忘因子的更新,否则在Tk-1温度下进行遗忘因子的更新。
5.根据权利要求1所述的一种考虑宽温度的锂电池参数辨识及状态联合估计方法,其特征在于,所述的UOCV-SOC关系式产生如下:在锂电池常用工作温度范围内(-10℃~50℃),每次以10℃为间隔分别进行:以电池厂家推荐的充放电倍率进行的电池容量实验、开路电压实验、动态工况实验(DST工况),在动态工况下测量电池电流i、电压UP和温度T;结合容量实验和开路电压实验,获得每个采样点的电池UOCV和SOC,采用下式拟合得到UOCV-SOC关系式:
UOCV=k0+k1SOC+k2SOC2+k3SOC3+k4log(SOC)+k5log(1-SOC) (29)
式中,k0~k5为待确定的参数。
6.根据权利要求1所述的一种考虑宽温度的锂电池参数辨识及状态联合估计方法,其特征在于,所述的电池空间状态方程如下:
根据基尔霍夫定律得到锂电池一阶RC等效电路模型为
UL=UOCV-IR0-UP (30)
Figure FDA0003907284490000051
式中,UL为电池端电压,UOCV为电池开路电压,UP为极化电压,I是电池电流,R0为电池内阻,Rp、Cp分别为极化内阻和极化电容;
电池模型经过离散化处理后得到电池系统的状态空间方程为
Figure FDA0003907284490000052
[UL,k]=UOCV(SOCk)-Up,k-R0Ik+vk (33)
式中,下标k表示第k时刻的数据,Ts为实验采集数据时间,wk与vk是均值为零、协方差为Q和R的高斯白噪声,η表示库伦效率,QN表示电池容量。
CN202211309229.4A 2022-10-25 2022-10-25 一种考虑宽温度的锂电池参数辨识及状态联合估计方法 Pending CN115598540A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211309229.4A CN115598540A (zh) 2022-10-25 2022-10-25 一种考虑宽温度的锂电池参数辨识及状态联合估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211309229.4A CN115598540A (zh) 2022-10-25 2022-10-25 一种考虑宽温度的锂电池参数辨识及状态联合估计方法

Publications (1)

Publication Number Publication Date
CN115598540A true CN115598540A (zh) 2023-01-13

Family

ID=84848158

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211309229.4A Pending CN115598540A (zh) 2022-10-25 2022-10-25 一种考虑宽温度的锂电池参数辨识及状态联合估计方法

Country Status (1)

Country Link
CN (1) CN115598540A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116540104A (zh) * 2023-05-22 2023-08-04 车百中汽科技(北京)有限公司 一种考虑温度影响的锂离子电池健康度确定方法及系统
CN117074966A (zh) * 2023-10-19 2023-11-17 东方电子股份有限公司 基于衰减记忆卡尔曼滤波器的电池soc估计方法及系统

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116540104A (zh) * 2023-05-22 2023-08-04 车百中汽科技(北京)有限公司 一种考虑温度影响的锂离子电池健康度确定方法及系统
CN116540104B (zh) * 2023-05-22 2024-03-22 车百中汽科技(北京)有限公司 一种考虑温度影响的锂离子电池健康度确定方法及系统
CN117074966A (zh) * 2023-10-19 2023-11-17 东方电子股份有限公司 基于衰减记忆卡尔曼滤波器的电池soc估计方法及系统

Similar Documents

Publication Publication Date Title
CN113466723B (zh) 确定电池荷电状态的方法及装置,电池管理系统
CN111142036B (zh) 基于容量增量分析的锂离子电池在线快速容量估计方法
CN110398691B (zh) 基于改进自适应双无迹卡尔曼滤波器的锂离子动力电池SoC估计方法
CN110488194B (zh) 一种基于电化学阻抗模型的锂电池soc估算方法及其系统
CN111505506A (zh) 一种多尺度卡尔曼滤波与无迹卡尔曼滤波融合的电池soc估算方法
CN115598540A (zh) 一种考虑宽温度的锂电池参数辨识及状态联合估计方法
He et al. State-of-charge estimation of lithium ion batteries based on adaptive iterative extended Kalman filter
CN112946499A (zh) 基于机器学习的锂电池健康状态及荷电状态联合估算方法
CN109669132B (zh) 一种基于变分贝叶斯滤波的电池荷电状态估计方法
CN109839596B (zh) 基于ud分解的自适应扩展卡尔曼滤波的soc估算方法
CN111060824B (zh) 一种基于模型与算法分步融合的荷电状态估计方法
CN113156321B (zh) 一种锂离子电池荷电状态soc的估算方法
CN114740386B (zh) 一种基于健康状态的锂离子电池荷电状态估计方法
CN112989690B (zh) 一种混合动力汽车锂电池多时间尺度荷电状态估计方法
CN109917299B (zh) 一种锂电池荷电状态的三层滤波估算方法
CN109633470B (zh) 基于ekf-gpr和日常片段数据的电池实时全充时间的估算方法
CN112269133B (zh) 一种基于预充电路模型参数识别的soc估计方法
CN112098849A (zh) 基于求积分卡尔曼滤波的锂电池剩余电量估计方法
Yan et al. Electric vehicle battery soc estimation based on gnl model adaptive kalman filter
CN113466725B (zh) 确定电池荷电状态的方法及装置,存储介质及电子设备
CN115656848A (zh) 一种基于容量修正的锂电池soc估算方法
CN110412472B (zh) 一种基于正态伽马滤波的电池荷电状态估计方法
CN113495214B (zh) 一种基于温度变化模型的超级电容荷电状态估计方法
CN113805062B (zh) 锂电池等效电路模型参数在线鲁棒自适应辨识方法
CN114861545A (zh) 基于rnn神经网络与多参数约束的锂电池sop在线估算方法

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