CN112757961B - 锂电池soc估算方法、主动均衡控制方法及系统 - Google Patents

锂电池soc估算方法、主动均衡控制方法及系统 Download PDF

Info

Publication number
CN112757961B
CN112757961B CN202011544477.8A CN202011544477A CN112757961B CN 112757961 B CN112757961 B CN 112757961B CN 202011544477 A CN202011544477 A CN 202011544477A CN 112757961 B CN112757961 B CN 112757961B
Authority
CN
China
Prior art keywords
soc
battery
equalization
value
battery pack
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
CN202011544477.8A
Other languages
English (en)
Other versions
CN112757961A (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.)
Nanjing Institute of Technology
Original Assignee
Nanjing Institute of Technology
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 Nanjing Institute of Technology filed Critical Nanjing Institute of Technology
Priority to CN202011544477.8A priority Critical patent/CN112757961B/zh
Publication of CN112757961A publication Critical patent/CN112757961A/zh
Application granted granted Critical
Publication of CN112757961B publication Critical patent/CN112757961B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
    • B60L58/00Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles
    • B60L58/10Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling batteries
    • B60L58/12Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling batteries responding to state of charge [SoC]
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
    • B60L58/00Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles
    • B60L58/10Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling batteries
    • B60L58/18Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling batteries of two or more battery modules
    • B60L58/22Balancing the charge of battery modules
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
    • B60L2240/00Control parameters of input or output; Target parameters
    • B60L2240/40Drive Train control parameters
    • B60L2240/54Drive Train control parameters related to batteries
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/60Other road transportation technologies with climate change mitigation effect
    • Y02T10/70Energy storage systems for electromobility, e.g. batteries

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Power Engineering (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)

Abstract

本发明公开了一种锂电池SOC估算方法、主动均衡控制方法及系统,提高了电池组的工作效率与使用寿命。本发明所述的锂电池SOC估算方法采用平滑滤波算法对强跟踪平方根容积滤波算法进行改进,解决了现有SOC算法存在的计算精度低、平滑性差的问题,同时减轻了计算机运算的负担,提高了计算效率;针对多组串联电池组充电过程与电池间电流均衡过程,提出优化的均衡控制方法,并利用改进的SOC估算法解算出的SOC值作为控制判据,进一步提高了均衡控制的精准度,增强了锂电池管理系统控制效果。

Description

锂电池SOC估算方法、主动均衡控制方法及系统
技术领域
本发明涉及一种锂电池充放电管理方法,尤其涉及一种锂电池SOC估算方法、主动均衡控制方法及系统。
背景技术
随着电动汽车的飞速发展,锂电池因在能量密度、功率特性和使用寿命的优势成为电动汽车的理想动力源。为保障动力电池正常工作,电池管理系统需对动力电池的充放电过程进行控制,进而提高电池组的工作效率与使用寿命。
锂电池荷电状态(SOC,State ofCharge)可以从本质反应锂电池间的不一致性,提高SOC估算的精度有利于增强电池管理系统的控制效果。近年来随着SOC估算技术的快速发展,以SOC为控制判据的均衡策略也得到长足的发展。强跟踪平方根容积滤波算法结合了强跟踪滤波器与平方根容积滤波算法,解决了复杂模型的滤波问题,是一个成熟有效的算法,但仍存在计算精度较低、解算平滑性差的问题。
目前,锂电池管理系统(BMS,BatteryManagementSystem)发展日益成熟,控制效果得到改善。电池均衡管理技术作为其核心技术,能够有效提高充放电过程均衡性,近年在被动均衡与主动均衡两种工作模式上均取得突破。被动均衡电路结构简单,可靠性高,整体造价便宜,效率与均衡速度较低;主动均衡电路结构复杂,效率与均衡速度高,但其工作需要结合相应的控制方法,进而对控制器的要求更高。但是,现有的锂电池管理系统仍然存在充放电过程均衡性较差的问题。
发明内容
发明目的:针对以上问题,本发明提出一种锂电池SOC估算方法、主动均衡控制方法及系统,能够解决现有SOC算法存在的计算精度低、平滑性差的问题,并进一步增强电池管理系统的控制效果,实现对锂动力电池的有效管理,提高电池组的工作效率与使用寿命。
技术方案:本发明所采用的技术方案是一种锂电池SOC估算方法,采用平滑滤波算法对强跟踪平方根容积滤波算法进行改进,在计算锂电池状态的平滑估算值时引入平滑增益,采用固定区间平滑滤波模式,滤波和平滑过程同时迭代进行,实现实时SOC估算。
所述的在计算锂电池状态的平滑估算值时引入平滑增益,平滑估算值
Figure BDA0002852956120000011
计算式为:
Figure BDA0002852956120000012
Figure BDA0002852956120000021
式中,
Figure BDA0002852956120000022
为k时刻下锂电池状态的平滑估算值,
Figure BDA0002852956120000023
为k+1时刻下锂电池状态的平滑初始值;
Figure BDA0002852956120000024
为k时刻下锂电池状态的估计值,
Figure BDA0002852956120000025
为k+1时刻下锂电池状态的预测状态值,
Figure BDA0002852956120000026
为k+1时刻下锂电池状态的估计值;其中平滑增益按下式计算:
Figure BDA0002852956120000027
式中,Ks k为k时刻的平滑增益,Ps xx为互协方差矩阵,Sk+1|k为带渐消因子的预测误差协方差阵特征平方根。
基于上述改进的SOC估算方法,本发明还提出一种锂电池主动均衡控制方法,包括以下步骤:
步骤1、利用改进的强跟踪平方根容积滤波算法对电池SOC进行实时估算,所述改进的强跟踪平方根容积滤波算法是采用平滑滤波算法对强跟踪平方根容积滤波算法进行改进,在计算锂电池状态的平滑估算值时引入平滑增益,采用固定区间平滑滤波模式,滤波和平滑过程同时迭代进行,实现实时SOC估算;
步骤2、针对多组串联电池组的充电过程,以步骤1所得的SOC估算值为均衡变量,控制各个串联电池组的充电电流;同时,针对串联电池组内充放电过程,采用考虑测量噪声的SOC-电压混合均衡控制方法控制串联电池组内各单体电池间的均衡电路的工作状态。
其中,步骤(2)所述的控制各个串联电池组的充电电流,包括以下步骤:
(21)确定第n组的初始充电电流In 0
Figure BDA0002852956120000028
式中,In 0为第n组的初始充电电流,I1 0为第1组初始充电电流,SOCaim为充电预期,取值为0.85~0.95,
Figure BDA0002852956120000029
为初始充电时第n组电池串联电池组中SOC估算值的中位数,初始值最小的串联电池组记为第1组,
Figure BDA00028529561200000210
为初始充电时第1组串联电池组的SOC估算值的中位数;所述的第1组初始充电电流I1 0计算式为:
Figure BDA00028529561200000211
式中,α为充电接受比,η为充放电效率,CN为电池标称容量。
(22)确定第n组在t时刻的协调充电电流In t
Figure BDA00028529561200000212
式中,In t为第n组在t时刻的协调充电电流,
Figure BDA0002852956120000031
为第n组电池SOC估算值在t时刻的中位数。
步骤(2)中所述的考虑测量噪声的SOC-电压混合均衡控制方法,包括以下步骤:
(31)根据电池组内SOC估算值的众数SOCmode选择采用均衡模式;所述的根据电池组内SOC估算值的众数SOCmode选择采用均衡模式,包括电压均衡模式和SOC均衡模式,选择条件为:
Figure BDA0002852956120000032
(32)计算状态判断函数F;
所述的计算第i次测量时的状态判断函数值Fi,在电压均衡模式下,计算式为:
Figure BDA0002852956120000033
式中,βset为电压均衡使能阈值,
Figure BDA0002852956120000034
其中Vmax i为第i次测量时电池组内单体电池最高电压,
Figure BDA0002852956120000035
为第i次测量时电池组的平均电压值,βi-1为第i-1次最高电压与第i-1次电压平均值之差;
在SOC均衡模式下,计算第i次测量时状态判断函数值Fi的计算式为:
Figure BDA0002852956120000036
式中,δset为SOC均衡使能阈值;δi为第i次测量时电池组SOC估算值方差,δi-1为第i-1次测量时电池组SOC估算值的方差。
(33)计算判断函数ρ,计算式为:
Figure BDA0002852956120000037
式中,ρi为第i次测量时判断函数值,用于控制均衡电路的开关,ρi+1为第i+1次测量时判断函数值;Fi为第i次测量时状态判断函数值;Fi+1为第i+1次测量时状态判断函数值;判断函数值ρi的下限为0,当计算值为负数时记为0;
(34)判断均衡电路开关:当判断函数值ρi增长到设定阈值时,均衡电路开启,判断函数值ρi重置为0,当均衡过程结束后返回步骤(31)重新计算。
本发明提出一种锂电池主动均衡控制系统,包括处理器和均衡电路,处理器执行上述锂电池主动均衡控制方法中的步骤,用以控制串联电池组内各单体电池间的均衡电路的工作状态。
有益效果:相比于现有技术,本发明具有以下优点:(1)在强跟踪平方根容积滤波算法的基础上,结合平滑滤波算法对强跟踪平方根容积滤波算法进行改进,引入平滑增益计算平滑估算值,解决了现有SOC算法存在的计算精度低、平滑性差的问题,同时减轻了计算机运算的负担,提高了计算效率;(2)针对多组串联电池组的充电过程,结合最优充电理论,以SOC为均衡变量形成多组串联电池组最优协调充电电流控制策略,提高了均衡过程中的充电速率;(3)针对充放电过程中组内电池间剩余容量不一致导致的电流均衡问题,引入判断函数ρ,形成一种考虑测量噪声的SOC-电压混合均衡策略,减少测量噪声产生的均衡电路误动作,提高了均衡电路的工作效率。上述控制策略结合改进SOC估计值作为控制判据,进一步提高了控制精度。本发明实现了SOC实时高精度计算,并以此为基础建立了一种新型均衡控制方法,进而提高电池组的工作效率与使用寿命。
附图说明
图1是本发明所述的锂电池SOC估算方法的流程图;
图2是Thevenin等效电路模型图;
图3是本发明所述的考虑测量噪声的SOC-电压混合均衡策略的流程图。
具体实施方式
下面结合附图和实施例对本发明的技术方案作进一步的说明。
本发明所述的锂电池SOC估算方法,利用改进的强跟踪平方根容积滤波算法对电池SOC进行实时估算。该方法采用平滑滤波算法对强跟踪平方根容积滤波算法进行改进,引入平滑增益Ks k,采用固定区间平滑滤波模式,滤波和平滑过程同时迭代进行,实现SOC实时估算。图1展示了改进后的算法流程,包括以下具体步骤:
步骤1.1,如图2所示,结合Thevenin等效电路模型具体讲述如何将电池SOC估算与滤波算法结合的过程。这里只是展示了一种等效电路的状态变量、输入变量、输出变量。针对其他的模型有不同的具体表现形式,但是主体思想是相类似的。采用Thevenin等效电路模型,建立锂电池离散的状态方程:
Figure BDA0002852956120000041
以及锂电池离散的量测方程:
Figure BDA0002852956120000051
式中,Up,k为k时刻的极化电压;Up,k-1为k-1时刻的极化电压;SOCk为k时刻锂电池的荷电状态;SOCk-1为k-1时刻锂电池的荷电状态;Ut,k为k时刻的电池端电压;Ik-1为k-1时刻的负载电流;Ik为k时刻的负载电流;Uoc,k为k时刻的开路电压,由SOC-OCV曲线确定;Rp为极化内阻;R0为欧姆内阻;Δt为采样间隔;τ为时间常数,数值等于RpCp,Cp为极化电容;CN为电池额定容量;w1,k-1、w2,k-1是k-1时刻的一些不可测量的随机输入量对系统状态变量的干扰,为系统噪声;vk为k时刻电池端电压的量测噪声。
系统噪声与量测噪声为互不相关的零均值高斯白噪声,有噪声统计特性:
Figure BDA0002852956120000052
Figure BDA0002852956120000053
Figure BDA0002852956120000054
式中,wk为k时刻的系统噪声矩阵;cov为协方差算子;Qk为k时刻的系统噪声协方差阵;Rk为k时刻的量测噪声协方差阵。
根据Thevenin等效电路模型,结合式(1)和式(2)得到输入滤波算法中的状态方程:
Figure BDA0002852956120000055
与量测方程:
Figure BDA0002852956120000056
式中,
Figure BDA0002852956120000057
为k时刻的状态变量;
Figure BDA0002852956120000058
为1-1时刻的状态变量;uk-1为k-1时刻的系统输入;uk为k时刻的系统输入;
Figure BDA0002852956120000059
为k时刻的系统输出;f(·)为系统的状态函数;h(·)为系统的量测函数。各参数的物理含义由式(8)-(12)显示:
Figure BDA00028529561200000510
Figure BDA00028529561200000511
uk-1=[Ik-1] (10)
uk=[Ik] (11)
Figure BDA0002852956120000061
步骤1.2,参数初始化,输入初值包括初始状态x0|0、初始预测误差协方差阵P0|0和初始噪声协方差阵Q0和R0,设定遗忘因子ρ和弱化因子η。
步骤1.3,对P0|0进行乔列斯基分解,得到预测误差协方差阵特征平方根初值S0|0
S0|0=chol(P0|0)T (13)
式中,S0|0为预测误差协方差阵特征平方根初值;chol(·)为乔列斯基分解。
步骤1.4,估计预测状态值
Figure BDA0002852956120000062
计算k-1时刻的状态容积点,
Figure BDA0002852956120000063
式中,Xi,k-1|k-1为k-1时刻
Figure BDA0002852956120000064
的状态容积点;Sk-1|k-1为k-1时刻预测误差协方差矩阵的特征平方根,当k=1时,Sk-1|k-1=S0|0;ξi为状态容积点向量;
Figure BDA0002852956120000065
为k-1时刻的状态估计值,当k=1时,
Figure BDA0002852956120000066
Figure BDA0002852956120000067
式中,X* i,k|k-1为状态容积点Xi,k-1|k-1上的预测状态值;uk-1为k-1时刻的输入向量;f(·)为系统的状态函数。
Figure BDA0002852956120000068
式中,
Figure BDA0002852956120000069
为k时刻的预测状态值;m为状态容积点的个数。
步骤1.5,估计带渐消因子的预测误差协方差阵特征平方根Sk|k-1
Figure BDA00028529561200000610
式中,Qk-1为k-1时刻的系统噪声协方差阵,当k=1时,Qk-1=Q0;SQ,k-1与ST Q,k-1为Qk-1分解得到,ST Q,k-1为SQ,k-1的转置。
Figure BDA00028529561200000611
其中,χ* k|k-1表示为:
Figure BDA00028529561200000612
式中,qr(·)为QR分解,Qd和Rd分别为矩阵经QR分解得到的正交阵与上三角阵;λk-1为k-1时刻的强跟踪渐消因子,当k=1时λ0=1。
Sk|k-1=Rd(1:n,:)T (20)
式中,Sk|k-1为带渐消因子的预测误差协方差阵特征平方根;Rd(1:n,:)T表示取矩阵Rd的前n行和前n列转置后赋值给Sk|k-1;n为矩阵χ* k|k-1的行数。
步骤1.6,估计预测量测值
Figure BDA0002852956120000071
计算更新后的状态容积点及容积点上的预测量测值,
Figure BDA0002852956120000072
Zi,k|k-1=h(Xi,k|k-1,uk) (22)
式中,Xi,k|k-1为更新后的状态容积点;ξi为状态容积点向量;Zi,k|k-1为更新后状态容积点上的预测量测值;h(·)为系统的量测函数;uk为k时刻的输入向量。
Figure BDA0002852956120000073
式中,
Figure BDA0002852956120000074
为k时刻预测量测值。
步骤1.7,估计新息协方差特征平方根Szz,k|k-1和互协方差阵Pxz,k|k-1
Figure BDA0002852956120000075
式中,Rk为k时刻量测噪声协方差阵;SR,k与ST R,k为Rk分解得到,ST R,k为SR,k的转置。
Figure BDA0002852956120000076
其中,Z* k|k-1表示为:
Figure BDA0002852956120000077
式中,Qd和Rd分别为矩阵经QR分解得到的正交阵与上三角阵。
Szz,k|k-1=Rd(1:l,:)T (27)
Figure BDA0002852956120000078
式中,Szz,k|k-1为新息协方差特征平方根;Rd(1:l,:)T表示取矩阵Rd的前l行和前l列转置后赋值给Szz,k|k-1;l为矩阵Z* k|k-1的行数;Pxz,k|k-1为互协方差阵。
步骤1.8,更新强跟踪渐消因子λk
Figure BDA0002852956120000079
其中,λk计算为:
Figure BDA0002852956120000081
其中,Nk和Mk分别计算为:
Figure BDA0002852956120000082
Figure BDA0002852956120000083
式中,λk为k时刻的强跟踪渐消因子,在初始时设置为1,即λ0=1;tr(·)为求迹运算;ST k|k-1为Sk|k-1的转置;ST zz,k|k-1为Szz,k|k-1的转置;η为弱化因子,取η>0,作用是为避免有可能造成的过调节,一般取η=2。
其中,Vk为实际输出的残差序列的协方差阵,计算为:
Figure BDA0002852956120000084
Figure BDA0002852956120000085
式中,εk为k时刻输出残差序列;zk为k时刻的实际量测值;ρ为遗忘因子,一般取ρ=0.95。
步骤1.9,计算滤波增益Kk
Figure BDA0002852956120000086
式中,Kk为k时刻的滤波增益。
步骤1.10,估计状态估计值
Figure BDA0002852956120000087
Figure BDA0002852956120000088
式中,
Figure BDA0002852956120000089
为k时刻的状态估计值。
步骤1.11,估计预测误差协方差阵的特征平方根Sk|k
Figure BDA00028529561200000810
Sk|k=Rd(1:n,:)T (38)
其中,χk|k-1表示为:
Figure BDA00028529561200000811
式中,Qd和Rd分别为矩阵经QR分解得到的正交阵与上三角阵;Sk|k为k时刻预测误差协方差阵的特征平方根;Rd(1:n,:)T表示取矩阵Rd的前n行和前n列转置后赋值给Sk|k;n为矩阵χk|k-1的行数。
步骤1.12,为平滑初始值赋值:
Figure BDA0002852956120000091
式中,xs k为k时刻的平滑初始值,用于计算k-1时刻的平滑估算值。
并参考步骤1.4,估计预测状态值
Figure BDA0002852956120000092
Figure BDA0002852956120000093
Figure BDA0002852956120000094
Figure BDA0002852956120000095
式中,Xi,k|k为k时刻
Figure BDA0002852956120000096
的状态容积点;Sk|k为k时刻预测误差协方差矩阵的特征平方根;ξi为状态容积点向量;X* i,k+1|k为状态容积点Xi,k|k上的预测状态值;uk为k时刻的输入向量;
Figure BDA0002852956120000097
为k+1时刻预测状态值。
步骤1.13,估计互协方差矩阵Ps xx和带渐消因子的预测误差协方差阵特征平方根Sk+1|k。Ps xx的计算式为:
Figure BDA0002852956120000098
其中,χk|k和χs k+1|k表示为:
Figure BDA0002852956120000099
Figure BDA00028529561200000910
式中,Ps xx为互协方差矩阵。
参考步骤1.5,估计带渐消因子的预测误差协方差阵特征平方根Sk+1|k
Figure BDA00028529561200000911
Figure BDA00028529561200000912
Figure BDA00028529561200000913
Sk+1|k=Rd(1:n,:)T (50)
式中,Qk为k时刻的系统噪声协方差阵;SQ,k与ST Q,k为Qk分解得到,ST Q,k为SQ,k的转置;Qd和Rd分别为矩阵经QR分解得到的正交阵与上三角阵;Rd(1:n,:)T表示取矩阵Rd的前n行和前n列转置后赋值给Sk+1|k;n为矩阵χ* k+1|k的行数;λk为k时刻的强跟踪渐消因子;Sk+1|k为更新的带渐消因子的预测误差协方差阵特征平方根。
步骤1.14,计算平滑增益Ks k
Figure BDA0002852956120000101
式中,Ks k为k时刻的平滑增益,用于计算k时刻的平滑估算值。
步骤1.15,计算平滑估算值xs k|k
Figure BDA0002852956120000102
Figure BDA0002852956120000103
式中,xs k|k为k时刻下的平滑估算值,为k时刻状态变量的最终估计值;xs k+1为k+1时刻的平滑初始值;
Figure BDA0002852956120000104
为k+1时刻的状态估计值。由上述步骤可以看出,平滑解算过程蕴藏在强跟踪平方根容积滤波过程中,并且系统对于上一时刻存储的值较少,在提高计算精度、增强解算平滑性的同时减轻了计算机运算的负担,也保障了实时性。
步骤2、针对电池组串联充电过程,提出以SOC为均衡变量的多组串联电池组最优协调充电策略。同时,针对串联电池组内充放电过程,采用考虑测量噪声的SOC-电压混合均衡控制方法控制串联电池组内各单体电池间的均衡电路的工作状态。
步骤2中,针对电池组串联充电过程,结合最优充电电流理论与安时积分公式,建立以SOC为均衡变量的多组串联电池组最优协调充电电流策略。
在锂电池充电时,最优充电电流理论表示为:
I=I0e-αt (54)
式中,I为充电电流;I0为初始充电电流;α为充电接受比;t为充电时间。
根据安时积分公式,
Figure BDA0002852956120000105
结合式(54),最终化简为:
Figure BDA0002852956120000106
式中,SOC(t)为充电t时刻电池的荷电状态;SOC(0)为充电初始时电池的荷电状态;η为充放电效率;CN为电池标称容量。
可以考虑,当充电过程结束时,充电电流I近似为0,而此时的电池荷电状态为1,可以表示为:
Figure BDA0002852956120000111
结合式(56)和式(57),得到基于SOC的最优充电电流模型:
Figure BDA0002852956120000112
传统最优充电电流理论的均衡变量为时间,但考虑到充电时电池间容量均衡过程的影响,以及电池容量衰退对充电时间的影响,按照时间调整充电电流可能会造成电池的过充或无法充满就退出充电,而SOC能本质反应锂电池间的不一致性及电池本身的状态,基于步骤1的改进SOC估算方法,结合式(58),提出以SOC为均衡变量的多组串联电池组最优协调充电策略,包括以下具体步骤:
步骤2.1,确定单体电池参数,包括充电接受比α,充放电效率η和电池标称容量CN
步骤2.2,对各电池进行SOC估算,记录初始充电时第n组串联电池组中SOC估算值的中位数为SOCm,n 0,初始值最小的组记为第1组,即第1组串联电池组初始充电时SOC估算值的中位数记为SOCm,1 0,上标m表示该组SOC估算值的中位数。
步骤2.3,引入充电预期SOCaim,确定第1组初始充电电流I1 0
Figure BDA0002852956120000113
式中,SOCaim为充电预期,通常选取在0.85-0.95之间,当充电达到充电预期时,可转为涓流充电对电池组进行充电,在提高充电效率的同时也对电池进行保护,延长电池的使用寿命。
步骤2.4,为保障各电池组达到充电预期的时间相同,减少电池组间均衡电路的动作,结合安时积分公式,确定第n组的初始充电电流In 0
Figure BDA0002852956120000114
式中,In 0为第n组的初始充电电流。
步骤2.5,确定各串联电池组充电策略:
Figure BDA0002852956120000121
式中,In t为第n组在t时刻协调充电电流;SOCm,n t为第n组电池SOC估算值在t时刻的中位数。
步骤2中,针对充放电过程中组内电池间剩余容量不一致导致的电流均衡过程,提出一种考虑测量噪声的SOC-电压混合均衡策略。
将流过每一个电池的电流Ic分为两部分,一部分为流经该串联电池组的充放电电流,记为I0,这部分电流为主电路流经各电池的电流,幅值较大并且每个电池都相等;另一部分为促使各电池容量变化速率一致,调节各电池容量的均衡电流,记为I1,该电流的流通通道为辅助电路中的均衡电路,幅值较小且每个电池流过的电流不相同。均衡控制方法是控制均衡电路的开启与关闭的,在必要时为均衡电流提供通路,进而降低电池间的不一致性,在完成均衡后关闭,减少损耗。总体而言,Ic=I0+I1。本发明所述的电流均衡过程,即是指涉及到均衡电流的部分。
针对充放电过程中组内电池间剩余容量不一致导致的电流均衡过程,引入判断函数ρ对均衡电路的开关进行控制,提出一种考虑测量噪声的SOC-电压混合均衡策略。
观察电池的充放电曲线时,发现SOC在20%到80%间电池电压的一致性较好,采用电压均衡效果较差,此时应当引入SOC均衡;当SOC大于80%或SOC小于20%时,电压的一致性较差,并考虑到倘如充电末期进行涓流充电,SOC估计精度降低,此时应当采用电压均衡。图3展示了SOC-电压混合均衡策略的流程,包括以下具体步骤:
步骤(31),采集电池组内各电池的电压,并对各电池SOC进行估算。根据电池组内SOC估算值的众数SOCmode选择采用均衡模式:
Figure BDA0002852956120000122
步骤(32),计算状态判断函数F:
在电压均衡过程中,记第i次测量时,电池组内单体电池最高电压记为Vmax i,电池组的平均电压值为
Figure BDA0002852956120000123
记最高电压减去电压平均值为βi
Figure BDA0002852956120000124
则第i次测量时状态判断函数值Fi为:
Figure BDA0002852956120000131
式中,βset为电压均衡使能阈值;βi-1为第i-1次最高电压与第i-1次电压平均值之差。
在SOC均衡过程中,记第i次测量时,各电池SOC估算值记为SOCj i,电池组SOC估算值的方差记为δi
Figure BDA0002852956120000132
Figure BDA0002852956120000133
式中,n为电池组内单体电池的个数。
则第i次测量时状态判断函数值Fi为:
Figure BDA0002852956120000134
式中,δset为SOC均衡使能阈值;δi-1为第i-1次测量时电池组SOC估算值的方差。
当Fi为-1时,说明组内电池间的差异性正在变大,电池间进行均衡的必要性增加;Fi为1时,说明组内电池间的差异性正在变小,电池间进行均衡的必要性减小;Fi为0时,系统无法判断电池间的差异性变化情况,处于观察状态。状态判断函数值F的引入有利于减少因测量噪声产生的误动作,增强了控制系统的鲁棒性。
步骤(33),计算判断函数ρ:
Figure BDA0002852956120000135
式中,ρi为第i次测量时判断函数值,用于控制均衡电路的开关,ρi+1为第i+1次测量时判断函数值;Fi+1为第i+1次测量时状态判断函数值;其中当Fi×Fi+1<0时,说明状态判断函数输出出错,为避免均衡电路误动作,判断函数应重新计数。ρi的下限为0,若计算值出现负数,仍然记ρi为0。
步骤(34),判断均衡电路开关:
均衡电路的开启条件为:
ρi=ρset (69)
式中,ρset为均衡电路开启阈值,取值应当适中。当判断函数值达到ρset时均衡电路开启,ρi记为0,当均衡过程结束时再重新计数。均衡电路是现有电路,按照能量的转移方式可分为电容式、电感式、变压器式、DC-DC变换器式以及多电平变换器式等,典型的电路有开关电容均衡电路、LC谐振均衡电路等。均衡过程是指,自均衡电路开启,至均衡电路完成均衡,即均衡电路中电流为0时结束这一整个过程,在一次充放电过程中往往会有多次均衡过程。
引入判断函数ρ,结合考虑测量噪声的SOC-电压混合均衡策略,能够有效减少测量噪声产生的均衡电路误动作,以提高均衡电路的工作效率。

Claims (7)

1.一种锂电池主动均衡控制方法,其特征在于,包括以下步骤:
步骤1、对电池组内各电池SOC进行估算;
步骤2、针对多组串联电池组的充电过程,以步骤1所得的SOC估算值为均衡变量,控制各个串联电池组的充电电流;同时,针对串联电池组内充放电过程,采用考虑测量噪声的SOC-电压混合均衡控制方法控制串联电池组内各单体电池间的均衡电路的工作状态;
其中,步骤(2)所述的控制各个串联电池组的充电电流,包括以下步骤:
(21)确定第n组的初始充电电流In 0
Figure FDA0003567479920000011
式中,In 0为第n组的初始充电电流,I1 0为第1组初始充电电流,SOCaim为充电预期,取值为0.85~0.95,
Figure FDA0003567479920000012
为初始充电时第n组串联电池组中SOC估算值的中位数,初始值最小的串联电池组记为第1组,
Figure FDA0003567479920000013
为初始充电时第1组串联电池组的SOC估算值的中位数;
(22)确定第n组在t时刻的协调充电电流In t
Figure FDA0003567479920000014
式中,In t为第n组在t时刻的协调充电电流,
Figure FDA0003567479920000015
为第n组电池SOC估算值在t时刻的中位数。
2.根据权利要求1所述的锂电池主动均衡控制方法,其特征在于:步骤(21)中所述的第1组初始充电电流I1 0计算式为:
Figure FDA0003567479920000016
式中,α为充电接受比,η为充放电效率,CN为电池标称容量。
3.根据权利要求1所述的锂电池主动均衡控制方法,其特征在于,步骤(2)所述的考虑测量噪声的SOC-电压混合均衡控制方法,包括以下步骤:
(31)根据电池组内SOC估算值的众数SOCmode选择均衡模式;
(32)计算状态判断函数F;
(33)计算判断函数ρ,计算式为:
Figure FDA0003567479920000021
式中,ρi为第i次测量时判断函数值,用于控制均衡电路的开关,ρi+1为第i+1次测量时判断函数值;Fi为第i次测量时状态判断函数值;Fi+1为第i+1次测量时状态判断函数值;判断函数值ρi的下限为0,当计算值为负数时记为0;
(34)判断均衡电路开关:当判断函数值ρi增长到设定阈值时,均衡电路开启,判断函数值ρi重置为0,当均衡过程结束后返回步骤(31)重新计算。
4.根据权利要求3所述的锂电池主动均衡控制方法,其特征在于:步骤(31)中所述的根据电池组内SOC估算值的众数SOCmode选择均衡模式,包括电压均衡模式和SOC均衡模式,选择条件为:
Figure FDA0003567479920000022
5.根据权利要求3所述的锂电池主动均衡控制方法,其特征在于:步骤(32)中所述的状态判断函数F,在电压均衡模式下,计算式为:
Figure FDA0003567479920000023
式中,Fi为第i次测量时的状态判断函数值,βset为电压均衡使能阈值,
Figure FDA0003567479920000024
其中Vmax i为第i次测量时电池组内单体电池最高电压,
Figure FDA0003567479920000025
为第i次测量时电池组的平均电压值,βi-1为第i-1次最高电压与第i-1次电压平均值之差;
在SOC均衡模式下,计算第i次测量时状态判断函数值Fi的计算式为:
Figure FDA0003567479920000026
式中,δset为SOC均衡使能阈值;δi为第i次测量时电池组SOC估算值方差,δi-1为第i-1次测量时电池组SOC估算值的方差。
6.根据权利要求1所述的锂电池主动均衡控制方法,其特征在于,步骤1中所述的对电池组内电池进行SOC估算,是利用改进的强跟踪平方根容积滤波算法对电池SOC进行实时估算,所述改进的强跟踪平方根容积滤波算法是采用平滑滤波算法对强跟踪平方根容积滤波算法进行改进,在计算锂电池状态的平滑估算值时引入平滑增益,采用固定区间平滑滤波模式,滤波和平滑过程同时迭代进行,实现实时SOC估算。
7.一种锂电池主动均衡控制系统,包括处理器和均衡电路,其特征在于:处理器执行权利要求1-6任一项所述的锂电池主动均衡控制方法中的步骤,用以控制串联电池组内各单体电池间的均衡电路的工作状态。
CN202011544477.8A 2020-12-23 2020-12-23 锂电池soc估算方法、主动均衡控制方法及系统 Active CN112757961B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011544477.8A CN112757961B (zh) 2020-12-23 2020-12-23 锂电池soc估算方法、主动均衡控制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011544477.8A CN112757961B (zh) 2020-12-23 2020-12-23 锂电池soc估算方法、主动均衡控制方法及系统

Publications (2)

Publication Number Publication Date
CN112757961A CN112757961A (zh) 2021-05-07
CN112757961B true CN112757961B (zh) 2022-06-03

Family

ID=75695485

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011544477.8A Active CN112757961B (zh) 2020-12-23 2020-12-23 锂电池soc估算方法、主动均衡控制方法及系统

Country Status (1)

Country Link
CN (1) CN112757961B (zh)

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FI123467B (fi) * 2011-07-08 2013-05-31 Europ Batteries Oy Menetelmä ja järjestelmä litium-ioni-kennoston varaustilan hallinnoimiseksi
US9077181B2 (en) * 2013-01-11 2015-07-07 GM Global Technology Operations LLC Battery section balancing methods and systems
DE102015202939A1 (de) * 2015-02-18 2016-08-18 Robert Bosch Gmbh Vorrichtung und Verfahren zum Ausgleichen des Ladezustands von Batteriezellen sowie Batteriemodul, Batterie, Batteriesystem, Fahrzeug, Computerprogramm und Computerprogrammprodukt
CN105406530A (zh) * 2015-11-27 2016-03-16 中科动力(福建)新能源汽车有限公司 一种锂电池单体电量的均衡方法及均衡电路
CN108761340A (zh) * 2018-05-31 2018-11-06 天津工业大学 基于噪声干扰的强跟踪容积卡尔曼滤波的电池估算方法
CN109038707B (zh) * 2018-07-05 2020-05-19 华中科技大学 一种电池组分段混合均衡控制方法
CN109017381B (zh) * 2018-07-31 2021-09-24 电子科技大学 一种动力电池组复合均衡控制方法
CN109617169A (zh) * 2018-12-20 2019-04-12 西安交通大学 一种基于模糊控制的锂离子电池均衡方法
CN110824363B (zh) * 2019-10-21 2022-01-11 江苏大学 一种基于改进ckf的锂电池soc和soe联合估算方法
CN111707953A (zh) * 2019-11-24 2020-09-25 华南理工大学 一种基于后向平滑滤波框架的锂电池soc在线估计方法
CN111509806A (zh) * 2020-04-29 2020-08-07 广东电网有限责任公司东莞供电局 电池均衡管理方法、装置、设备及存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《动力锂电池均衡充电技术研究》;杨红华;《工程科技Ⅱ辑》;20121231(第12期);全文 *
《电动汽车锂离子电池组均衡充电技术研究》;肖峰;《工程科技Ⅱ辑》;20120131(第1期);全文 *

Also Published As

Publication number Publication date
CN112757961A (zh) 2021-05-07

Similar Documents

Publication Publication Date Title
CN107368619B (zh) 一种扩展卡尔曼滤波soc估算方法
CN110303945B (zh) 一种蓄电池组电量自适应优化平衡控制方法
CN109061506A (zh) 基于神经网络优化ekf的锂离子动力电池soc估计方法
CN109633456B (zh) 一种基于分段电压识别法的动力锂电池组soc估算方法
CN109031147B (zh) 一种磷酸铁锂电池组的soc估算方法
CN112557906A (zh) 一种动力电池全生命周期内soc和容量在线联合估计方法
CN112649747A (zh) 一种分数阶扩展卡尔曼的锂电池soc估算方法
CN113777510A (zh) 一种锂电池荷电状态估计方法及装置
CN109346787B (zh) 一种电动汽车动力电池自适应优化充电方法
CN108872865B (zh) 一种防滤波发散的锂电池soc估算方法
CN115598541B (zh) 基于遗忘因子自适应反馈修正的电池能量状态评估方法
CN114217234B (zh) 一种基于ide-asrckf的锂离子电池参数辨识与soc估计方法
CN113807039A (zh) 一种串联型电池系统的功率状态预测方法
CN113777501A (zh) 一种电池模组的soh的估算方法
CN112580289A (zh) 一种混合电容器功率状态在线估计方法及系统
CN115906613A (zh) 一种基于粒子群和自适应模糊神经网络的电池组均衡方法
CN114814591A (zh) 一种锂电池soe估算方法、装置及系统
CN112946480B (zh) 一种提高soc估计实时性的锂电池电路模型简化方法
CN111044924A (zh) 一种全工况电池剩余容量确定方法及系统
CN110232432B (zh) 一种基于人工生命模型的锂电池组soc预测方法
CN112757961B (zh) 锂电池soc估算方法、主动均衡控制方法及系统
CN111293378B (zh) 一种提高锂离子电池循环寿命的充电方法
TWI527337B (zh) The method of determining the charging current of multi - stage lithium battery using Taguchi orthogonal table and fuzzy algorithm
CN116718920A (zh) 基于rnn优化扩展卡尔曼滤波的锂电池soc估算方法
CN116559673A (zh) 一种基于igwo-ekf算法的锂电池soc估算方法

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