CN112816876B - 一种用于可充电电池的低温电池剩余电量估算方法及装置 - Google Patents

一种用于可充电电池的低温电池剩余电量估算方法及装置 Download PDF

Info

Publication number
CN112816876B
CN112816876B CN202011581305.8A CN202011581305A CN112816876B CN 112816876 B CN112816876 B CN 112816876B CN 202011581305 A CN202011581305 A CN 202011581305A CN 112816876 B CN112816876 B CN 112816876B
Authority
CN
China
Prior art keywords
battery
state
discrete time
temperature
low
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
Application number
CN202011581305.8A
Other languages
English (en)
Other versions
CN112816876A (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.)
Hunan Aerospace Jiecheng Electronic Equipment Co ltd
Original Assignee
Hunan Aerospace Jiecheng Electronic Equipment Co ltd
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 Hunan Aerospace Jiecheng Electronic Equipment Co ltd filed Critical Hunan Aerospace Jiecheng Electronic Equipment Co ltd
Priority to CN202011581305.8A priority Critical patent/CN112816876B/zh
Publication of CN112816876A publication Critical patent/CN112816876A/zh
Application granted granted Critical
Publication of CN112816876B publication Critical patent/CN112816876B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/382Arrangements for monitoring battery or accumulator variables, e.g. SoC
    • 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/389Measuring internal impedance, internal conductance or related variables

Landscapes

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

Abstract

本发明提供一种用于可充电电池的低温电池剩余电量估算方法及装置,所述装置包括电量估算单元、温度监控单元、加热控制单元、加热单元、电池、电压电流监控单元,所述电量估算单元建立算法模型,结合温度监控单元对数据进行拟合,提升低温下电池电量估算准确度。在估算中,广泛采用的是利用电池组电阻分量的阻抗电流积分法,通过对电池的建模,首先,选择电池的状态参数。在此基础上,通过输入电流和输出电压得到SOC的状态变量OCV(Zk),建立了电池的状态方程,可以对其剩余电量进行测量。然后利用扩展卡尔曼滤波算法利用测量电压y值估计电池的状态x,即通过扩展卡尔曼滤波算法估计电池的SOC,并使状态估计观测值
Figure DDA0002865102180000011
的中误差最小。

Description

一种用于可充电电池的低温电池剩余电量估算方法及装置
技术领域
本发明属于低温可充电电池剩余电量估算技术领域,具体涉及一种用于可充电电池的低温电池剩余电量估算方法及装置。
背景技术
近年来,充电锂电池应用范围迅速扩大,但锂电池制造工艺决定电池可放电电量与温度成直接关系,为保证电池低温下具备足够容量,通常会对电池进行加热,同时电池在放电过程中也会产生自温升,实际情况是低温放电初期,电池处于放电状态,但电池剩余电量(SOC)增加,此现象较影响用户体验。
发明内容
本发明针对上述缺陷,提供一种借助机器算法,建立数据模型,充分考虑电池自放电的热耗及加热功耗,通过预计方法实现电池低温下SOC预估,并避免放电状态下电池剩余电量SOC增加的可充电电池的低温电池剩余电量估算方法及装置。
本发明提供如下技术方案:一种用于可充电电池的低温电池剩余电量估算方法,包括以下步骤:
S1:采用温度监控单元采集每个离散时间k的温度T并传递给电量估算单元,采用电压电流监控单元采集每个时间步的所述电池的输出电压yk、输入电流ik、所述可充电电池充电状态下的电池内部电阻R+、所述可充电电池放电状态下的电池内部电阻R-和迟滞现象带来的电阻M并传递给所述电量估算单元;
S2:所述电量估算单元构建带有低通滤波器Rs的状态空间方程的开路电压电池模型,利用所述状态空间方程和状态变量对所述可充电电池于离散时间k时的状态x进行估计,构建电池建模参数θ,得到所述低通滤波器Rs于离散时间k+1时的滤波器状态向量fk+1和于离散时间k时的滤波器输出电压向量
Figure BDA0002865102160000011
S3:所述电量估算单元向加热控制单元发出指令,使所述加热控制单元控制所述加热单元将环境温度加热至不同的温度T下,所述温度T从25℃起始,所述温度监控单元和所述电压电流监控单元再次采集监控所述温度T下的各项参数并传递给所述电量估算单元,所述电量估算单元根据所述不同的温度T及其各项参数构建的不同电池建模参数θ3、θ2、θ1、θ0再次构建根据温度进行电磁建模参数优化后的滤波器输出电压向量y(T):
y(T)=θ3T32T21T10
S4:采用离散时间k时的所述S3步骤得到的输出电压向量值y(k)估计电池于离散时间k+1时的电池状态x(k+1);
S5:构建离散时间k初始值为0时的电池状态x(k)的期望值
Figure BDA0002865102160000021
模型和协方差矩阵模型
Figure BDA0002865102160000022
并进行逐渐递增的离散之间k的运算;
S6:使用扩展卡尔曼滤波算法电池状态x的估计值
Figure BDA0002865102160000023
进行计算,得到的误差最小的估计值
Figure BDA0002865102160000024
为低温电池剩余电量SOC值:
Figure BDA0002865102160000025
进一步地,所述步骤S2包括以下步骤:
S21:构建带有低通滤波器Rs的状态空间方程的开路电压电池模型:
V电池=Vocv+Rik+Vp1+Vp2+Vh
其中,所述Vocv为开路电压,所述R为电池内阻,所述Vh为迟滞现象带来的电压,所述Vp1为第一电阻C1的电压,所述Vp2为第二电阻C2的电压;
S22:利用所述S21步骤得到的状态空间方程和状态变量对所述可充电电池的状态x=[fk+1 hk+1 Zk+1]T进行估计,得到所述可充电电池的电池剩余电量Zk对应的开路电压OCV(Zk),所述状态变量包括输入电流ik和输出电压yk
Figure BDA0002865102160000031
OCV(Zk)=yk-Rik-Mhk-filt(ik);
其中,所述fk为离散时间k时的所述低通滤波器状态向量,所述hk为离散时间k时的迟滞现象导致的电流,所述Zk为所述为离散时间k时的所述可充电电池的电池剩余电量;所述yk为模型于离散时间k的输出电压,所述filt(ik)为所述低通滤波器接收电池充电/放电产生的电流,所述M为迟滞现象带来的电阻,所述Cn为第n电阻的电池标称容量,所述Δt为监测离散时间间隔,所述Bf为所述低通滤波器的输入增益矩阵;
S23:根据代表所述滤波器极点的α1、所述低通滤波器输出常数g1,所述电池的磁滞率常数γ、第n电阻的电池标称容量Cn、所述可充电电池充电状态下的电池内部电阻R+、所述可充电电池放电状态下的电池内部电阻R-和迟滞现象带来的电阻M的向量得到电池建模参数θ,电池建模参数θ表示为θ=[α1,g1,γ,Cn,R+,R-,M];
S24:构建所述低通滤波器的线性状态空间方程:
fk+1=Affk+Bfik
Figure BDA0002865102160000032
其中,所述fk为离散时间k处的滤波器状态向量,所述ik为输入电流向量,所述
Figure BDA0002865102160000033
为滤波器输出电压向量,所述Af为所述低通滤波器的状态转移矩阵,所述Bf为所述低通滤波器的输入增益矩阵,所述G是滤波器的输出矩阵。
进一步地,所述F(ik)被定义为:
Figure BDA0002865102160000034
进一步地,所述S4步骤包括以下步骤:
x(k+1)=f(x(k),u(k))+wk
y(k)=g(x(k),u(k))+vk
Figure BDA0002865102160000041
Figure BDA0002865102160000042
其中,所述wk为服从正态分布的估算过程噪声,所述vk为服从正态分布的估算测量噪声,所述f(x(k),u(k))为所述电池于离散时间k的预测状态
Figure BDA0002865102160000043
的泰勒级数展开式,所述g(x(k),u(k))为所述电池于离散时间k的卡尔曼估算期望状态
Figure BDA0002865102160000044
的泰勒级数展开式,所述
Figure BDA0002865102160000045
为所述电池于离散时间k的预测系数矩阵,所述
Figure BDA0002865102160000046
为所述电池于离散时间k的卡尔曼估算期望系数矩阵。
进一步地,所述S5步骤包括以下步骤:
S51:当初始值k=0时,计算所述电池于离散时间k的卡尔曼估算期望状态
Figure BDA0002865102160000047
和所述协方差矩阵
Figure BDA0002865102160000048
Figure BDA0002865102160000049
Figure BDA00028651021600000410
其中,所述
Figure BDA00028651021600000411
是所述电池于离散时间k的状态xk的期望卡尔曼估算值,所述
Figure BDA00028651021600000412
是所述电池于于初始值0时的状态x0的期望卡尔曼估算值;
S52:进入k值运算,所述k=1,2,……,所述电池状态估计时间更新公式:
Figure BDA00028651021600000413
其中,所述
Figure BDA00028651021600000414
为所述电池于离散时间k-1时的状态xk-1的预测值,所述uk-1为所述电量估算单元(1)于离散时间k-1时的控制输入电流值;
S53:所述电池状态预测误差协方差
Figure BDA00028651021600000415
计算:
Figure BDA00028651021600000416
其中,所述
Figure BDA00028651021600000417
为所述电池于离散时间k-1时的状态xk-1的卡尔曼估算期望系数矩阵,所述
Figure BDA00028651021600000418
为所述
Figure BDA00028651021600000419
的转置矩阵,所述
Figure BDA00028651021600000420
为所述电池于离散时间k-1时的卡尔曼估计误差协方差矩阵,所述Sw为所述估算过程噪声矩阵;
S54:所述电池状态卡尔曼增益Lk计算:
Figure BDA0002865102160000051
其中,所述
Figure BDA0002865102160000052
为所述电池于离散时间k的预测系数矩阵,所述
Figure BDA0002865102160000053
为所述
Figure BDA0002865102160000054
的转置矩阵,所述Sv为估算测量噪声矩阵;
S55:所述电池状态估算计算:
Figure BDA0002865102160000055
其中,所述
Figure BDA0002865102160000056
为所述电池于离散时间k的预测状态,所述
Figure BDA0002865102160000057
为所述电池状态测量余量;
S56:所述电池状态误差协方差计算:
Figure BDA0002865102160000058
以进一步更新所述电池于离散时刻k时的卡尔曼估算协方差矩阵
Figure BDA0002865102160000059
本发明还提供采用上述可充电电池的低温电池剩余电量估算方法的估算装置,包括电量估算单元、温度监控单元、加热控制单元、加热单元、电池、电压电流监控单元,所述电量估算单元构建带有低通滤波器Rs的状态空间方程的开路电压电池模型,基于所述开路电压电池模型,所述电量估算单元向加热控制单元发出指令,使所述加热控制单元控制所述加热单元将环境温度加热至不同的温度T下,所述温度T从25℃起始,所述温度监控单元和所述电压电流监控单元再次采集监控所述温度T下的各项参数并传递给所述电量估算单元,结合所述温度监控单元对数据进行拟合,采用离散时间k时的得到的输出电压向量值y(k)估计电池于离散时间k+1时的电池状态x(k+1);构建离散时间k初始值为0时的电池状态x(k)的期望值
Figure BDA00028651021600000510
模型和协方差矩阵模型
Figure BDA00028651021600000511
并进行逐渐递增的离散之间k的运算;使用扩展卡尔曼滤波算法电池状态x的估计值
Figure BDA00028651021600000512
进行计算,得到的误差最小的估计值
Figure BDA00028651021600000513
为低温电池剩余电量SOC值,提升低温下电池电量估算准确度。
进一步地,所述温度监控单元采用热敏电阻NTC,通过上拉电阻进行分压,输出至电量估算单元。
进一步地,所述加热控制单元通过IO电路实现加热单元控制;所述加热单元通过加热膜,实现对电池单元加热。
进一步地,所述电池选用低温18650锂电池电芯,通过串并方式满足电池容量需求。
进一步地,电压电流监控单元选用多节电池的电池组监视器,实现对电池单体电压、总电压及电流的监控。
本发明的有益效果为:
1、本发明提供的估算方法考虑了温度特性的可充电电池的低温放电后电池荷电状态(SOC)估计算法。将可充电电池的非线性特性以数学形式表示为微分方程,通过电池建模,用状态空间方程表示,提取可充电电池的特性参数。采用扩展卡尔曼滤波(EKF)算法作为状态观测器,对电池状态进行估计,进而估算出可充电电池的剩余电量。
2、本发明提供的方法及装置借助机器算法,建立数据模型,充分考虑电池自放电的热耗及加热功耗,通过预计方法实现电池低温下SOC预估,并避免放电状态下电池剩余电量SOC增加。
3、本发明提供的装置中电量估算单元通过多次试验数据,对不同环境温度下电池温升速率建立模型,并根据加热单元的功耗数据以及环境温度,确定设备热平衡点,依据电池的温度-电量曲线、电压电量信息以及放电前的剩余电量(SOC),实现对电池剩余电量(SOC)精确预估。
附图说明
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中:
图1为本发明实施例1提供的方法构建的带有低通滤波器Rs的状态空间方程的开路电压电池模型示意图;
图2为本发明实施例2提供的可充电电池的低温电池剩余电量估算装置结构示意图。
具体实施例方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
本实施例提供一种用于可充电电池的低温电池剩余电量估算方法,包括以下步骤:
S1:采用温度监控单元2采集每个离散时间k的温度T并传递给电量估算单元1,采用电压电流监控单元6采集每个时间步的所述电池的输出电压yk、输入电流ik、所述可充电电池充电状态下的电池内部电阻R+、所述可充电电池放电状态下的电池内部电阻R-和迟滞现象带来的电阻M并传递给所述电量估算单元1;
S2:所述电量估算单元1构建带有低通滤波器Rs的状态空间方程的开路电压电池模型,利用所述状态空间方程和状态变量对所述可充电电池于离散时间k时的状态x进行估计,构建电池建模参数θ,得到所述低通滤波器Rs于离散时间k+1时的滤波器状态向量fk+1和于离散时间k时的滤波器输出电压向量
Figure BDA0002865102160000071
S3:所述电量估算单元1向加热控制单元3发出指令,使所述加热控制单元3控制所述加热单元4将环境温度加热至不同的温度T下,所述温度T从25℃起始,所述温度监控单元2和所述电压电流监控单元6再次采集监控所述温度T下的各项参数并传递给所述电量估算单元1,如表1所示,所述电量估算单元1根据所述不同的温度T及其各项参数构建的不同电池建模参数θ3、θ2、θ1、θ0再次构建根据温度进行电磁建模参数优化后的滤波器输出电压向量y(T):
y(T)=θ3T32T21T10
S4:采用离散时间k时的所述S3步骤得到的输出电压向量值y(k)估计电池于离散时间k+1时的电池状态x(k+1);
S5:构建离散时间k初始值为0时的电池状态x(k)的期望值
Figure BDA0002865102160000081
模型和协方差矩阵模型
Figure BDA0002865102160000082
并进行逐渐递增的离散之间k的运算;
S6:使用扩展卡尔曼滤波算法电池状态x的估计值
Figure BDA0002865102160000083
进行计算,得到的误差最小的估计值
Figure BDA0002865102160000084
为低温电池剩余电量SOC值:
Figure BDA0002865102160000085
表1
Figure BDA0002865102160000086
其中,步骤S2包括以下步骤:
S21:如图1所示,构建带有低通滤波器Rs的状态空间方程的开路电压电池模型:
V电池=Vocv+Rik+Vp1+Vp2+Vh
其中,所述Vocv为开路电压,所述R为电池内阻,所述Vh为迟滞现象带来的电压,所述Vp1为第一电阻C1的电压,所述Vp2为第二电阻C2的电压;
S22:利用所述S21步骤得到的状态空间方程和状态变量对所述可充电电池的状态x=[fk+1 hk+1 Zk+1]T进行估计,得到所述可充电电池的电池剩余电量Zk对应的开路电压OCV(Zk),所述状态变量包括输入电流ik和输出电压yk
Figure BDA0002865102160000091
OCV(Zk)=yk-Rik-Mhk-filt(ik);
其中,所述fk为离散时间k时的所述低通滤波器状态向量,所述hk为离散时间k时的迟滞现象导致的电流,所述Zk为所述为离散时间k时的所述可充电电池的电池剩余电量;所述yk为模型于离散时间k的输出电压,所述filt(ik)为所述低通滤波器接收电池充电/放电产生的电流,所述M为迟滞现象带来的电阻,所述Cn为第n电阻的电池标称容量,所述Δt为监测离散时间间隔,所述Bf为所述低通滤波器的输入增益矩阵;
F(ik)被定义为:
Figure BDA0002865102160000092
在该方程中,低通滤波器filt(ik)的函数使输出电压yk在经过一段静止时间后收敛到OCV(Zk),并且在恒流充放电时,输出yk收敛到包括OCV(Zk)在内的OCV(Zk)+hk+Rik
S23:由于电解质浓度、电池标称容量(Cn)、电池内部电阻估计值对于充电(R+)和放电(R-)等条件的不同,在极化和放电期间出现的滞后(hk)不同,因此需要根据代表所述滤波器极点的α1、所述低通滤波器输出常数g1,所述电池的磁滞率常数γ、第n电阻的电池标称容量Cn、所述可充电电池充电状态下的电池内部电阻R+、所述可充电电池放电状态下的电池内部电阻R-和迟滞现象带来的电阻M的向量得到电池建模参数θ,可表示为θ=[α1,g1,γ,Cn,R+,R-,M];
S24:构建所述低通滤波器的线性状态空间方程:
fk+1=Affk+Bfik
Figure BDA0002865102160000101
其中,所述fk为离散时间k处的滤波器状态向量,所述ik为输入电流向量,所述
Figure BDA0002865102160000102
为滤波器输出电压向量,所述Af为所述低通滤波器的状态转移矩阵,所述Bf为所述低通滤波器的输入增益矩阵,所述G是滤波器的输出矩阵。
S4步骤包括以下步骤:
x(k+1)=f(x(k),u(k))+wk
y(k)=g(x(k),u(k))+vk
Figure BDA0002865102160000103
Figure BDA0002865102160000104
其中,所述wk为服从正态分布的估算过程噪声,所述vk为服从正态分布的估算测量噪声,所述f(x(k),u(k))为所述电池于离散时间k的预测状态
Figure BDA0002865102160000105
的泰勒级数展开式,所述g(x(k),u(k))为所述电池于离散时间k的卡尔曼估算期望状态
Figure BDA0002865102160000106
的泰勒级数展开式,所述
Figure BDA0002865102160000107
为所述电池于离散时间k的预测系数矩阵,所述
Figure BDA0002865102160000108
为所述电池于离散时间k的卡尔曼估算期望系数矩阵。
S5步骤包括以下步骤:
S51:当初始值k=0时,计算所述电池于离散时间k的卡尔曼估算期望状态
Figure BDA0002865102160000109
和所述协方差矩阵
Figure BDA00028651021600001010
Figure BDA0002865102160000111
Figure BDA0002865102160000112
其中,所述
Figure BDA0002865102160000113
是所述电池于离散时间k的状态xk的期望卡尔曼估算值,所述
Figure BDA0002865102160000114
是所述电池于于初始值0时的状态x0的期望卡尔曼估算值;
S52:进入k值运算,所述k=1,2,……,所述电池状态估计时间更新公式:
Figure BDA0002865102160000115
其中,所述
Figure BDA0002865102160000116
为所述电池于离散时间k-1时的状态xk-1的预测值,所述uk-1为所述电量估算单元(1)于离散时间k-1时的控制输入电流值;
S53:所述电池状态预测误差协方差
Figure BDA0002865102160000117
计算:
Figure BDA0002865102160000118
其中,所述
Figure BDA0002865102160000119
为所述电池于离散时间k-1时的状态xk-1的卡尔曼估算期望系数矩阵,所述
Figure BDA00028651021600001110
为所述
Figure BDA00028651021600001111
的转置矩阵,所述
Figure BDA00028651021600001112
为所述电池于离散时间k-1时的卡尔曼估计误差协方差矩阵,所述Sw为所述估算过程噪声矩阵;
S54:所述电池状态卡尔曼增益Lk计算:
Figure BDA00028651021600001113
其中,所述
Figure BDA00028651021600001114
为所述电池于离散时间k的预测系数矩阵,所述
Figure BDA00028651021600001115
为所述
Figure BDA00028651021600001116
的转置矩阵,所述Sv为估算测量噪声矩阵;
S55:所述电池状态估算计算:
Figure BDA00028651021600001117
其中,所述
Figure BDA00028651021600001118
为所述电池于离散时间k的预测状态,所述
Figure BDA00028651021600001119
为所述电池状态测量余量;
S56:所述电池状态误差协方差计算:
Figure BDA00028651021600001120
以进一步更新所述电池于离散时刻k时的卡尔曼估算协方差矩阵
Figure BDA00028651021600001121
在估算中,广泛采用的是利用电池组电阻分量的阻抗电流积分法。本文通过对电池的建模,建立了电池的状态方程,然后利用状态观测器来估计电池的剩余容量。首先,选择电池的状态参数。在此基础上,通过输入电流ik和输出电压yk得到SOC的状态变量OCV(Zk)。系统由一个状态变量组成,通过输入和输出,可以对其进行测量。扩展卡尔曼滤波算法利用测量电压y值估计电池的状态x,即通过EKF估计电池的SOC,并使状态估计观测值
Figure BDA0002865102160000121
的中误差最小。
实施例2
如图2所示,为利用实施例1的一种用于可充电电池的低温电池剩余电量估算方法的估算装置,包括电量估算单元1、温度监控单元2、加热控制单元3、加热单元4、电池5、电压电流监控单元6,所述电量估算单元构建带有低通滤波器Rs的状态空间方程的开路电压电池模型,基于所述开路电压电池模型,所述电量估算单元向加热控制单元发出指令,使所述加热控制单元控制所述加热单元将环境温度加热至不同的温度T下,所述温度T从25℃起始,所述温度监控单元2和所述电压电流监控单元6再次采集监控所述温度T下的各项参数并传递给所述电量估算单元,结合所述温度监控单元对数据进行拟合,采用离散时间k时的得到的输出电压向量值y(k)估计电池于离散时间k+1时的电池状态x(k+1);构建离散时间k初始值为0时的电池状态x(k)的期望值
Figure BDA0002865102160000122
模型和协方差矩阵模型
Figure BDA0002865102160000123
并进行逐渐递增的离散之间k的运算;使用扩展卡尔曼滤波算法电池状态x的估计值
Figure BDA0002865102160000124
进行计算,得到的误差最小的估计值
Figure BDA0002865102160000125
为低温电池剩余电量SOC值,
本实施例的电量估算单元采用STM32F429,基于试验数据建立计算模型,模型分为只放电不加热、边放电边加热、先加热后放电,同时分析确认不同模型热平衡点,并依据电量温度曲线以及电压电流信息,进行数据拟合,通过机器算法实现电池剩余电量预计。STM32F429处理器内部集成内存(SRAM)、非遗失性存储器(ROM),单芯片可构成完整的、小巧的系统,同时支持正常工作模式(Normal)、低功耗模式(Stanby),可通过外部电路唤醒。
温度监控单元2采用热敏电阻NTCALUG01T103F,其电阻阻值根据温度进行相应变化,通过上拉电阻进行分压,输出至电量估算单元,其可靠性非常高。
加热控制单元3采用大电流固态继电器ASSR-1530,可对外提供60V/1A电流,驱动接触器吸合,通过IO电路实现加热单元,4控制,对加热膜进行供电;加热单元4选用电加热膜实现,加热膜由电阻丝、绝缘包覆层、引出导线和接插件组成,绝缘包覆层背面背胶保证安装可靠性,实现对电池单元5加热。
电池5选用国产宽温18650锂电池电芯SW18650-22SP,单节电芯4.2V/2200mAh,,通过串并方式满足电池容量需求,可在低温下稳定可靠放电。
电压电流监控单元6选用凌特公司第3代多节电池的电池组监视器LTC6804,可测量多达12个串接电池的电压,并具备低于1.2mV的总测量误差,实现对电池单体电压、总电压及电流的监控。
本发明采用电子器件,借助机器算法,建立数据模型,充分考虑电池自放电的热耗及加热功耗,通过预计方法实现电池低温下SOC预估,并避免放电状态下电池剩余电量(SOC)增加。
电量估算单元1通过多次试验数据,对不同环境温度下电池温升速率建立模型,并根据加热单元4的功耗数据以及环境温度,确定设备热平衡点,依据电池的温度-电量曲线、电压电量信息以及放电前的剩余电量SOC,实现对电池剩余电量SOC精确预估。
虽然已经参考优选实施例对本发明进行了描述,但在不脱离本发明的范围的情况下,可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是,只要不存在结构冲突,各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例,而是包括落入权利要求的范围内的所有技术方案。

Claims (6)

1.一种可充电电池的低温电池剩余电量估算方法,其特征在于,包括以下步骤:
S1:采用温度监控单元(2)采集每个离散时间k的温度T并传递给电量估算单元(1),采用电压电流监控单元(6)采集每个时间步的所述电池的输出电压yk、输入电流ik、所述可充电电池充电状态下的电池内部电阻R+、所述可充电电池放电状态下的电池内部电阻R-和迟滞现象带来的电阻M并传递给所述电量估算单元(1);
S2:所述电量估算单元(1)构建带有低通滤波器Rs的状态空间方程的开路电压电池模型,利用所述状态空间方程和状态变量对所述可充电电池于离散时间k时的状态x进行估计,构建电池建模参数θ,得到所述低通滤波器Rs于离散时间k+1时的滤波器状态向量fk+1和于离散时间k时的滤波器输出电压向量
Figure FDA0003300994290000011
S3:所述电量估算单元(1)向加热控制单元(3)发出指令,使所述加热控制单元(3)控制加热单元(4)将环境温度加热至不同的温度T下,所述温度T从25℃起始,所述温度监控单元(2)和所述电压电流监控单元(6)再次采集监控所述温度T下的各项参数并传递给所述电量估算单元(1),所述电量估算单元(1)根据所述不同的温度T及其各项参数构建的不同电池建模参数θ3、θ2、θ1、θ0再次构建根据温度进行电磁建模参数优化后的滤波器输出电压向量y(T):
y(T)=θ3T32T21T10
S4:采用离散时间k时的所述S3步骤得到的输出电压向量值y(k)估计电池于离散时间k+1时的电池状态x(k+1);
S5:构建离散时间k初始值为0时的电池状态x(k)的期望值
Figure FDA0003300994290000016
模型和协方差矩阵模型
Figure FDA0003300994290000013
并进行逐渐递增的离散时间k的运算;
S6:使用扩展卡尔曼滤波算法电池状态x的估计值
Figure FDA0003300994290000014
进行计算,得到的误差最小的估计值
Figure FDA0003300994290000015
为低温电池剩余电量SOC值:
Figure FDA0003300994290000021
所述步骤S2包括以下步骤:
S21:构建带有低通滤波器Rs的状态空间方程的开路电压电池模型:
V电池=Vocv+Rik+Vp1+Vp2+Vh
其中,所述Vocv为开路电压,所述R为电池内阻,所述Vh为迟滞现象带来的电压,所述Vp1为第一电阻C1的电压,所述Vp2为第二电阻C2的电压;
S22:利用所述S21步骤得到的状态空间方程和状态变量对所述可充电电池的状态x=[fk+1 hk+1 Zk+1]T进行估计,得到所述可充电电池的电池剩余电量Zk对应的开路电压OCV(Zk),所述状态变量包括输入电流ik和输出电压yk
Figure FDA0003300994290000022
OCV(Zk)=yk-Rik-Mhk-filt(ik);
其中,所述fk为离散时间k时的所述低通滤波器状态向量,所述hk为离散时间k时的迟滞现象导致的电流,所述Zk为所述为离散时间k时的所述可充电电池的电池剩余电量;所述yk为模型于离散时间k的输出电压,所述filt(ik)为所述低通滤波器接收电池充电/放电产生的电流,所述M为迟滞现象带来的电阻,所述Cn为第n电阻的电池标称容量,所述Δt为监测离散时间间隔,所述Bf为所述低通滤波器的输入增益矩阵;
S23:根据代表所述低通滤波器极点的α1、所述低通滤波器输出常数g1,所述电池的磁滞率常数γ、第n电阻的电池标称容量Cn、所述可充电电池充电状态下的电池内部电阻R+、所述可充电电池放电状态下的电池内部电阻R-和迟滞现象带来的电阻M的向量得到电池建模参数θ,所述电池建模参数θ表示为θ=[α1,g1,γ,Cn,R+,R-,M];
S24:构建所述低通滤波器的线性状态空间方程:
fk+1=Affk+Bfik
Figure FDA0003300994290000031
其中,所述fk为离散时间k处的滤波器状态向量,所述ik为输入电流向量,所述
Figure FDA0003300994290000032
为滤波器输出电压向量,所述Af为所述低通滤波器的状态转移矩阵,所述Bf为所述低通滤波器的输入增益矩阵,所述G是滤波器的输出矩阵;
所述S5步骤包括以下步骤:
S51:当初始值k=0时,计算所述电池于离散时间k的卡尔曼估算期望状态
Figure FDA0003300994290000033
和所述协方差矩阵
Figure FDA0003300994290000034
Figure FDA0003300994290000035
Figure FDA0003300994290000036
其中,所述
Figure FDA0003300994290000037
是所述电池于离散时间k的状态xk的期望卡尔曼估算值,所述
Figure FDA0003300994290000038
是所述电池于于初始值0时的状态x0的期望卡尔曼估算值;
S52:进入k值运算,所述k=1,2,......,所述电池状态估计时间更新公式:
Figure FDA0003300994290000039
其中,所述
Figure FDA00033009942900000310
为所述电池于离散时间k-1时的状态xk-1的预测值,所述uk-1为所述电量估算单元(1)于离散时间k-l时的控制输入电流值;
S53:所述电池状态预测误差协方差
Figure FDA00033009942900000321
计算:
Figure FDA00033009942900000312
其中,所述
Figure FDA00033009942900000313
为所述电池于离散时间k-1时的状态xk-1的卡尔曼估算期望系数矩阵,所述
Figure FDA00033009942900000314
为所述
Figure FDA00033009942900000315
的转置矩阵,所述
Figure FDA00033009942900000316
为所述电池于离散时间k-1时的卡尔曼估计误差协方差矩阵,所述Sw为所述估算过程噪声矩阵;
S54:电池状态卡尔曼增益Lk计算:
Figure FDA00033009942900000317
其中,所述
Figure FDA00033009942900000318
为所述电池于离散时间k的预测系数矩阵,所述
Figure FDA00033009942900000319
为所述
Figure FDA00033009942900000320
的转置矩阵,所述Sv为估算测量噪声矩阵;
S55:电池状态估算计算:
Figure FDA0003300994290000041
其中,所述
Figure FDA0003300994290000042
为所述电池于离散时间k的预测状态,所述
Figure FDA0003300994290000043
为所述电池状态测量余量;
S56:电池状态误差协方差计算:
Figure FDA0003300994290000044
以进一步更新所述电池于离散时刻k时的卡尔曼估算协方差矩阵
Figure FDA0003300994290000045
所述F(ik)被定义为:
Figure FDA0003300994290000046
所述S4步骤包括以下步骤:
x(k+1)=f(x(k),u(k))+wk
y(k)=g(x(k),u(k))+vk
Figure FDA0003300994290000047
Figure FDA0003300994290000048
其中,所述wk为服从正态分布的估算过程噪声,所述vk为服从正态分布的估算测量噪声,所述f(x(k),u(k))为所述电池于离散时间k的预测状态
Figure FDA0003300994290000049
的泰勒级数展开式,所述g(x(k),u(k))为所述电池于离散时间k的卡尔曼估算期望状态
Figure FDA00033009942900000410
的泰勒级数展开式,所述
Figure FDA00033009942900000411
为所述电池于离散时间k的预测系数矩阵,所述
Figure FDA00033009942900000412
为所述电池于离散时间k的卡尔曼估算期望系数矩阵。
2.应用权利要求1所述的低温电池剩余电量估算方法的一种用于可充电电池的低温电池剩余电量估算装置,其特征在于,包括电量估算单元(1)、温度监控单元(2)、加热控制单元(3)、加热单元(4)、电池(5)、电压电流监控单元(6),所述电量估算单元(1)构建带有低通滤波器Rs的状态空间方程的开路电压电池模型,基于所述开路电压电池模型,所述电量估算单元(1)向加热控制单元(3)发出指令,使所述加热控制单元(3)控制所述加热单元(4)将环境温度加热至不同的温度T下,所述温度T从25℃起始,所述温度监控单元(2)和所述电压电流监控单元(6)再次采集监控所述温度T下的各项参数并传递给所述电量估算单元(1),结合所述温度监控单元(2)对数据进行拟合,采用离散时间k时的得到的输出电压向量值y(k估计电池于离散时间k+1时的电池状态x(k+1;构建离散时间k初始值为0时的电池状态x(k的期望值
Figure FDA0003300994290000052
模型和协方差矩阵模型
Figure FDA0003300994290000051
并进行逐渐递增的离散之间k的运算;使用扩展卡尔曼滤波算法电池状态x的估计值
Figure FDA0003300994290000053
进行计算,得到的误差最小的估计值
Figure FDA0003300994290000054
为低温电池剩余电量SOC值,提升低温下电池电量估算准确度。
3.根据权利要求2所述的一种用于可充电电池的低温电池剩余电量估算装置,其特征在于,所述温度监控单元(2)采用热敏电阻NTC,通过上拉电阻进行分压,输出至电量估算单元。
4.根据权利要求2所述的一种用于可充电电池的低温电池剩余电量估算装置,其特征在于,所述加热控制单元(3)通过IO电路实现加热单元(4)控制;所述加热单元(4)通过加热膜,实现对电池单元(5)加热。
5.根据权利要求2所述的一种用于可充电电池的低温电池剩余电量估算装置,其特征在于,所述电池(5)选用低温18650锂电池电芯,通过串并方式满足电池容量需求。
6.根据权利要求2所述的一种用于可充电电池的低温电池剩余电量估算装置,其特征在于,电压电流监控单元(6)选用多节电池的电池组监视器,实现对电池单体电压、总电压及电流的监控。
CN202011581305.8A 2020-12-28 2020-12-28 一种用于可充电电池的低温电池剩余电量估算方法及装置 Active CN112816876B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011581305.8A CN112816876B (zh) 2020-12-28 2020-12-28 一种用于可充电电池的低温电池剩余电量估算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011581305.8A CN112816876B (zh) 2020-12-28 2020-12-28 一种用于可充电电池的低温电池剩余电量估算方法及装置

Publications (2)

Publication Number Publication Date
CN112816876A CN112816876A (zh) 2021-05-18
CN112816876B true CN112816876B (zh) 2021-12-07

Family

ID=75854222

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011581305.8A Active CN112816876B (zh) 2020-12-28 2020-12-28 一种用于可充电电池的低温电池剩余电量估算方法及装置

Country Status (1)

Country Link
CN (1) CN112816876B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114497818B (zh) * 2022-01-26 2023-10-24 中南大学 一种低温环境下锂离子电池的内部快速加热方法
CN114454745B (zh) * 2022-03-01 2024-03-08 雅迪科技集团有限公司 一种锂电池低温充电系统、方法及电动二轮车
CN115563704B (zh) * 2022-09-23 2023-12-08 四川新能源汽车创新中心有限公司 电池状态预测模型的优化方法、容量预测方法及相关装置

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103616647A (zh) * 2013-12-09 2014-03-05 天津大学 一种用于电动汽车电池管理系统的电池剩余电量估计方法
EP2985618A1 (en) * 2005-12-27 2016-02-17 Toyota Jidosha Kabushiki Kaisha Charged state estimating device and charged state estimating method of secondary battery
CN105699910A (zh) * 2016-04-21 2016-06-22 中国计量大学 一种锂电池剩余电量在线估计方法
CN106814329A (zh) * 2016-12-30 2017-06-09 深圳市麦澜创新科技有限公司 一种基于双卡尔曼滤波算法的电池soc在线估计方法
CN107741569A (zh) * 2017-11-16 2018-02-27 温州大学 一种基于分段扩展卡尔曼滤波的锂电池荷电状态的估算方法
WO2018162021A1 (en) * 2017-03-06 2018-09-13 Volvo Truck Corporation A battery cell state of charge estimation method and a battery state monitoring system
CN109669134A (zh) * 2019-02-27 2019-04-23 浙江科技学院 一种基于卡尔曼滤波法的soc的估算方法
CN110395141A (zh) * 2019-06-27 2019-11-01 武汉理工大学 基于自适应卡尔曼滤波法的动力锂电池soc估算方法
CN110441702A (zh) * 2019-07-31 2019-11-12 湘潭大学 一种用扩展卡尔曼滤波估算锂离子电池充电电量的方法
CN111537894A (zh) * 2020-05-29 2020-08-14 合肥工业大学 一种用于估计锂电池soc和sop的方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2985618A1 (en) * 2005-12-27 2016-02-17 Toyota Jidosha Kabushiki Kaisha Charged state estimating device and charged state estimating method of secondary battery
CN103616647A (zh) * 2013-12-09 2014-03-05 天津大学 一种用于电动汽车电池管理系统的电池剩余电量估计方法
CN105699910A (zh) * 2016-04-21 2016-06-22 中国计量大学 一种锂电池剩余电量在线估计方法
CN106814329A (zh) * 2016-12-30 2017-06-09 深圳市麦澜创新科技有限公司 一种基于双卡尔曼滤波算法的电池soc在线估计方法
WO2018162021A1 (en) * 2017-03-06 2018-09-13 Volvo Truck Corporation A battery cell state of charge estimation method and a battery state monitoring system
CN107741569A (zh) * 2017-11-16 2018-02-27 温州大学 一种基于分段扩展卡尔曼滤波的锂电池荷电状态的估算方法
CN109669134A (zh) * 2019-02-27 2019-04-23 浙江科技学院 一种基于卡尔曼滤波法的soc的估算方法
CN110395141A (zh) * 2019-06-27 2019-11-01 武汉理工大学 基于自适应卡尔曼滤波法的动力锂电池soc估算方法
CN110441702A (zh) * 2019-07-31 2019-11-12 湘潭大学 一种用扩展卡尔曼滤波估算锂离子电池充电电量的方法
CN111537894A (zh) * 2020-05-29 2020-08-14 合肥工业大学 一种用于估计锂电池soc和sop的方法

Also Published As

Publication number Publication date
CN112816876A (zh) 2021-05-18

Similar Documents

Publication Publication Date Title
CN112816876B (zh) 一种用于可充电电池的低温电池剩余电量估算方法及装置
CN112098851B (zh) 智能电池与其荷电状态在线估计方法及应用
Qays et al. Recent progress and future trends on the state of charge estimation methods to improve battery-storage efficiency: A review
Gholizadeh et al. Estimation of state of charge, unknown nonlinearities, and state of health of a lithium-ion battery based on a comprehensive unobservable model
CN109633454B (zh) 一种实现锂离子电池等效温度在线估计的方法
CN112444749B (zh) 一种基于温度修正模型的锂电池荷电状态联合估计方法
CN111914503B (zh) 一种锂离子电池功率输入电热耦合模型建立方法
CN106443473A (zh) 一种动力锂离子电池组soc估算方法
CN106887877B (zh) 一种基于电池能量状态估计的电池组主动均衡控制系统
CN104849675A (zh) 锂离子电池电化学和热耦合模型的获取方法
Li et al. A new parameter estimation algorithm for an electrical analogue battery model
CN113253140A (zh) 电池健康状态在线估算方法
EP0505333B1 (en) Estimating the charge of batteries
CN110816366A (zh) 适用于单体电池内部的温度估算方法、系统、介质及设备
CN109752660B (zh) 一种无电流传感器的电池荷电状态估计方法
Lu et al. Modeling discharge characteristics for predicting battery remaining life
CN108169687A (zh) 一种基于云平台的蓄电池soc估算方法
CN108829911A (zh) 一种开路电压与soc函数关系优化方法
CN114740357A (zh) 一种并联电池组支路电流、荷电状态和功率状态的联合估计方法
CN110632520A (zh) 一种动力电池soc的估算装置及其估算方法
CN114779107A (zh) 一种考虑温度影响的锂离子电池soc估计方法
CN114928133A (zh) 一种用于锂离子电池的充电控制方法、系统、设备和介质
CN117293975B (zh) 一种锂电池的充放电调节方法、装置、设备及存储介质
CN113325327B (zh) 基于内阻测试的动力电池瞬态产热率的测算方法
CN114089193A (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