CN106772067B - 多时间尺度估计动力电池荷电状态与健康状态的方法 - Google Patents
多时间尺度估计动力电池荷电状态与健康状态的方法 Download PDFInfo
- Publication number
- CN106772067B CN106772067B CN201611077392.7A CN201611077392A CN106772067B CN 106772067 B CN106772067 B CN 106772067B CN 201611077392 A CN201611077392 A CN 201611077392A CN 106772067 B CN106772067 B CN 106772067B
- Authority
- CN
- China
- Prior art keywords
- state
- parameter
- value
- power battery
- estimation
- 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
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/3644—Constructional arrangements
- G01R31/3648—Constructional arrangements comprising digital calculation means, e.g. for performing an algorithm
-
- 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/005—Testing of electric installations on transport means
- G01R31/006—Testing of electric installations on transport means on road vehicles, e.g. automobiles or trucks
- G01R31/007—Testing of electric installations on transport means on road vehicles, e.g. automobiles or trucks using microprocessors or computers
-
- 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/392—Determining battery ageing or deterioration, e.g. state of health
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Secondary Cells (AREA)
Abstract
本申请涉及利用IAPF滤波估计动力电池荷电状态与健康状态的方法,首先进行准备工作:建立容量‑SOC‑OCV三维响应面,然后利用多时间尺度IAPF滤波算法,使用宏观时间尺度获得动力电池***参数估计值、使用微观时间尺度估计***状态,最后提出荷电状态SOC和健康状态SOH的表征量,评估所述动力电池的荷电状态与健康状态,形成基于多时间尺度的动力电池参数和状态的联合估计方法,实现动力电池荷电状态和可用容量在不确定性应用环境中的精确联合估计,使得估计结果在全工作周期内更加稳定可靠且实现减小BMS计算量。
Description
技术领域:
本发明涉及动力电池管理技术领域,尤其是车载动力电池***参数辨识、荷电状态和健康状态估计领域。
背景技术:
在前专利申请CN201610375853.4首次提出了联合估计动力电池的荷电状态SOC和参数的联合估计方法,分别利用HF算法进行在线参数辨识,利用UKF算法和前述在线参数辨识的结果进行在线SOC估计,实现估计方法针对***的实时参数更新和状态估计。
然而,上述方法中,参数根据采集数据实时更新,但是参数是一个缓慢变化的过程,而状态可以认为是一个快速实时变化的过程,所以采用HF-UKF联合估计方法虽然可以根据电池的工作情况实时更新电池的参数以便于获取更准确的荷电状态SOC,但是实时更新缓慢变化的参数对于车载电池管理***(以下简称BMS)而言计算量比较大,而更合理的方法是采用不同时间尺度来更新电池的参数和荷电状态。此外,HF-UKF联合估计方法无法获取电池的健康状态。
基于此,本发明对上述方法进行改进,利用Improved Adeptive Particle Filter(以下简称IAPF)算法来进行动力电池***荷电状态与健康状态的联合估计。
本发明针对动力电池***状态量的快速时变特性与参数量的缓慢时变特性,采用IAPF算法在微观时间尺度估计动力电池的SOC,利用微观时间尺度估计结果,采用IAPF算法在宏观时间尺度估计动力电池的模型参数与可用容量,形成基于多时间尺度的动力电池SOC和容量的联合估计方法,实现动力电池SOC和容量在不确定性应用环境中的精确联合估计。相对于HF-UKF联合估计电池参数与状态方法,此方法减少了计算量,并可以在线获取电池的健康状态。
发明内容:
本发明的一种动力电池状态和参数估计的多时间尺度IAPF滤波方法,包括:
每个微观采样点进行状态估计,每隔L个微观采样点为一个宏观采样点并进行参数估计;进而实现在微观和宏观两个时间尺度的动力电池估计。
初始化:初始化所述滤波算法的宏观参数观测器和微观状态观测器的初始参数;
步骤①:状态估计:利用当前微观采样点k下***的电流值、上个微观采样点k-1的状态粒子集和上个宏观采样点l-1的参数估计值进行状态粒子集值预估,然后计算状态权重值,最后对状态权重值进行归一化处理,得到当前微观采样点k下的状态估计值;
步骤②:参数估计:通过更新后的当前***开路电压计算参数粒子集值预估、参数权重值,然后对参数权重值进行归一化处理,得到当前宏观采样点l下的参数估计值。
优选地,所述初始参数包括状态噪声和参数噪声;
步骤③:判断k+1是否超过该算法最大计算次数,超过,则停止所述多时间尺度IAPF滤波算法算法,如果未超过,则进行步骤④;
步骤④:通过计算状态噪声和参数噪声的协方差,更新所述状态噪声和所述参数噪声;
步骤⑤:针对状态粒子和参数粒子退化现象,进行重新采样生成新的状态粒子集和参数粒子集;
经过上述五步之后,把k+1作为新的当前微观采样点状态,判断k+1是否能被L整除,如果能,则进行步骤①;否则直接输出参数估计值。
优选地,状态为极化电压和荷电状态,参数包括容量、极化电容、极化电阻和欧姆内阻。
优选地,步骤②中利用最近一次宏观采样点的参数估计值、当前微观采样点k 状 态估计值和根据容量建立的容量-SOC-OCV(开路电压)三维响应面更新所述当前***开路电压,
优选地,所述步骤⑤中,利用随机重采样方法生成新的粒子集。
本发明还涉及一种联合估计动力电池***荷电状态与健康状态的方法,其特征在于:
首先,建立动力电池的容量-状态-开路电压三维响应面;
其次,数据在线获取,实时采集动力电池单体或动力电池组的电压值和电流值;
然后,利用在前所述的IAPF滤波算法,获得的当前所述状态估计值与所述参数估计值;
最后,在线荷电状态与健康状态提取,利用所述状态估计值与参数估计值,估计动力电池***荷电状态与健康状态。
本发明所提出的动力电池荷电状态和健康状态估计方法与传统方法相比具有以下优势:
(1)IAPF算法更适用于动力电池中的锂离子电池这种强非线性模型,估计精度更高,IAPF算法相比传统的APF算法可以有效提高SOC估计的收敛速度和鲁棒性;
(2)多时间尺度IAPF算法可以在SOC初值均不准确的情况下均能快收敛到真值,即实现了电池最大容量未知时,SOC的准确估计,解决了传统SOC估计算法以最大可用容量已知为前提而无法成功运用到实车上的难题;
(3)电池的容量与内阻均为衡量电池健康状态(SOH)的重要指标,因而上述联合估计算法在一定程度上实现了SOC与SOH的联合估计;
(4)容量-SOC-OCV三维响应面能够更好的实现电池在不同使用环境下的SOC和容量的联合估计;
(5)减少了车载电池管理***(BMS)的计算量。
附图说明:
图1动力电池荷电状态与健康状态的估计方法
图2多时间尺度IAPF算法;
图3动力电池Thevenin等效电路模型;
图4容量-SOC-OCV三维响应面;
图5UDDS工况电流和SOC参考值。其中:(a)、电流;(b)、SOC参考值;
图6多时间尺度IAPF算法端电压和SOC估计结果。其中:(a)、端电压预测值与测量值对比;(b)、端电压预测误差;(c)、SOC估计值与参考值对比;(d)、SOC估计误差;
图7多时间尺度IAPF算法容量估计结果。其中(a)、容量估计值与参考值对比;(b)、容量估计误差;
具体实施方式:
本发明说设计的动力电池***包括动力电池单体或者成组后的动力电池***。
本发明所述的一种基于多时间尺度的动力电池***荷电状态(以下简称SOC)与健康状态(以下简称SOH)的联合估计方法如附图1所示。
本发明***荷电状态指代时时变化的***指标,即电池SOC。而***的健康状态指相对于荷电状态变化而言的较慢的***指标,如电池容量和电池模型参数,其在一次完全充放电过程中几乎没有变化。本发明下述***非特殊说明皆指代动力电池***,***荷电状态优选地对应电池***的SOC。***健康状态对应电池***参数或容量,优选地对应电池***SOH或最大可用容量。
该联合估计方法包括以下四方面:容量-SOC-OCV(开路电压)三维响应面的建立、数据在线获取、多时间尺度IAPF算法以及在线SOC与SOH提取。下面分别对上述四个方面就行详细叙述:
准备工作:容量-SOC-OCV三维响应面的建立
作为SOC估计算法的修正曲线的SOC-OCV曲线通常在某一特定工况或温度条件下获得,但电池在不同的使用环境(不同温度和不同老化程度)下,该曲线会发生较为明显的变化,进而导致在不同使用环境下的SOC估计精度降低。本发明将温度、老化程度等因素对该曲线的影响直接反映到电池容量的差异之上,利用容量、SOC与OCV三者的关系作为容量与SOC联合估计算法的修正曲面。具体过程如下:
在不同电池容量(即温度、老化程度变化时)下进行开路电压试验,以获取不同电池容量下的SOC与OCV对应关系,采用组合模型(如式(1)所示)分别对不同容量下的SOC与OCV关系进行拟合,从而得到各个不同容量下的α0,α1,…,α6参数值,最后采用二次函数(如式(2)所示)对参数α0,α1,…,α6与容量的关系进行拟合,至此完成容量-SOC-OCV三维响应面的建立。
Uoc(Ca,z)=α0+α1z+α2z2+α3z3+α4/z+α5ln(z)+α6ln(1-z) (1)
Ca为电池容量;
z为电池SOC;
Uoc(Ca,z)表示开路电压OCV,其表示为电池容量与SOC的函数;
α0,α1,…,α6为组合模型的系数;
上标T表示矩阵的转置;
Λ为7×3常数矩阵。
本发明使用Thevenin动力电池等效电路模型为例来阐述该动力电池SOC与SOH联合估计方法。上述Thevenin电池模型也可以采用现有技术中的其他模型替代:电化学模型、状态方程模型和其他结构的等效电路模型等。图3为Thevenin动力电池等效电路模型,该模型由电压源、欧姆内阻、以及RC网络三部分组成。其相应数学模型如式(3)所示。
Up为极化电压,为其导数;
Cp为极化电容;
Rp为极化电阻;
iL为输入电流;
Ut为端电压;
Uoc为开路电压;
R0为欧姆内阻。
动力电池SOC的计算方程为:
z0表示SOC的初值;
Ca为动力电池最大可用容量(下文简称为容量),同时电池最大可用容量是表征电池健康状态(SOH)的重要参数,即相同使用条件下,电池最大可用容量越小,电池衰退越明显,电池健康状态(SOH)越差。
在式(3)所示的动力电池***的数学模型的基础上,由于采样时间是离散的,因此基于IAPF算法,本发明建立了适用于动力电池***的IAPF算法的非线性离散***,所述IAPF算法的非线性离散***如式(5)所示:
x表示***的状态向量,xk=[Up,k zk]T,设xk服从一阶马尔科夫过程,即tk时刻的概率仅与tk-1时刻的概率相关;
θl表示***的参数向量,θl=[R0,l Rp,l Cp,l Ca,l]T;
y表示***的测量向量,yk=Ut,k;
u表示***的输入向量,uk=iL,k;
下标k表示tk时刻***采样时间点,同时也代表了状态估计的时间尺度,即在每个采样时间点下均进行一次状态估计。微观时间尺度,即所述状态估计的时间尺度;
下标l表示参数估计的时间尺度,其数值等于k除以L的商(L为时间尺度转换限值),即每隔L个采样时间点进行一次参数辨识,且每次参数辨识结果被用来估计tl×L时刻之后的L个时刻下的状态值。宏观时间尺度,即所述参数估计的时间尺度;
f(xk-1,θl,uk-1)表示模型的状态函数;
g(xk,θl,uk)表示模型的观测函数;
wk-1和ρl-1分别为***状态噪声和参数噪声,vk为测量噪声,在IAPF算法之中,所述***状态噪声、参数噪声和测量噪声被设计为随机且未知的,突破了传统滤波算法状态噪声、参数噪声和测量噪声为白噪声这一假设,因而与实际生产结合更加紧密。在本发明中,我们假设所有噪声均服从正态高斯分布的白噪声,因此***状态噪声和参数噪声的均值皆为0,更新状态噪声和参数噪声的协方差即能够实现更新***状态噪声和参数噪声。
Δt表示时间尺度k的单位时间间隔;
η(iL,k-1)表示充放电效率。
至此,已完成动力电池非线性离线***中各相关参数的定义。
1、数据在线获取
当电动汽车运行时,动力电池***中的BMS能够实时采集动力电池单体或动力电池组的电压、电流等信息,并储存在相应的存储器,为下面的多时间尺度IAPF算法提供实时信息输入,所述信息输入包括tk时刻***的测量值yk=Ut,k,tk时刻***的输入信息uk=iL,k。其中iL,k为控制电流;Ut,k为端电压。
2、多时间尺度IAPF算法
本发明使用多时间尺度IAPF算法来实现动力电池参数与状态联合估计。
下面对该算法具体过程进行描述:
算法的初始化:分别设置宏观参数观测器IAPFθ和微观状态观测器IAPFx的初始参数值。包括:
p(x0)、N、p(θ0)、M、w0、ρ0、v0 (6)
p(x0)为已知的***状态向量的先验概率密度函数,由初始的所述***状态向量的先验概率密度函数p(x0)产生k=0时刻***状态向量粒子集每个粒子的权重值均为1/N;
N为基于p(x0)随机产生的初始状态向量粒子数,由用户选择参数N作为在计算量和估计精度之间的权衡;
p(θ0)为已知的***参数向量的先验概率密度函数,由初始的所述***参数向量的先验概率密度函数p(θ0)产生k=0时刻***参数向量粒子集每个粒子的权重值均为1/M;
M为基于p(θ0)随机产生的初始参数向量粒子数,由用户选择参数M作为在计算量和估计精度之间的权衡;
w0和ρ0分别为初始的***状态噪声和参数噪声,v0为初始的测量噪声,用户可以自行根据经验设定这些噪声的初值,因为IAPF算法可以在噪声未知或初值不准确的情况下自动更新噪声以达到收敛效果;
当采样时间k∈{1,2,...,∞}时,基于电流、电压等信息的不断输入,计算:
步骤①:基于微观时间尺度的状态观测器IAPFx的状态估计
利用微观采样点k下的***的电流值、最近一次微观采样点k-1下的状态向量粒子集和最近一次容量估计宏观采样点l-1下的参数向量估计值,进行微观采样点k下的状态向量粒子集值预估,然后进行重要性采样即计算状态向量粒子权重值,最后对***状态向量权重值进行归一化处理,进而得到微观采样点k下的状态向量估计值。
状态向量粒子集值预估:
重要性采样即状态向量粒子权重值:
状态向量粒子权重值归一化:
状态向量估计值:tk-1时刻即微观采样点k-1;
为第tk时刻的状态向量粒子集,为tk-1时刻状态向量粒子集,为tk-1时刻***参数向量估计值,uk-1为tk-1时刻***的输入信息,其为已知量;
w1,k i为第i个粒子在tk时刻的状态权重值;
w11,k i为归一化后的第i个粒子在tk时刻的状态权重值;
yk和yk i分别为***真实量测值和***量测估计值,***量测估计值是把第i个粒子的状态向量代入量测方程得出的值;
R1为***状态噪声wk的协方差值;
为tk时刻状态向量估计值;
得到上述值后,进行状态向量粒子退化现象的判断。通过公式(7、8、10)计算状态向量粒子的权重Neff1,如果小于权重阈值(Nth1),则判断发生状态向量粒子退化现象。
步骤②:基于微观时间尺度的状态观测器IAPFθ的参数更新
利用最近一次容量估计宏观采样点的***参数向量估计值、当前微观采样点tk时刻的状态向量估计值和容量-SOC-OCV三维响应面,更新当前***开路电压得到开路电压OCVk。
进行宏观采样点l下的参数向量粒子集值预估,然后进行重要性采样即参数向量粒子权重值计算,之后计算并归一化参数向量权重值,进而可以计算出宏观采样点l的参数向量估计值。
参数向量粒子集值预估:
重要性采样即参数向量粒子权重值计算:
参数向量粒子权重值归一化:
参数向量估计值: 为第tk时刻的参数向量的粒子集,为tk-1时刻参数向量的粒子集,为步骤②中输出的状态向量估计值,uk-1为tk-1时刻***的输入信息,其为已知量;w2,k j为第j个粒子在tk时刻的参数权重值;
w21,k j为归一化后的第j个粒子在tk时刻的参数权重值;
yk和yk j分别为***真实量测值和***量测估计值,***量测估计值是把第j个粒子的参数向量代入量测方程得出的值;
R2为***参数噪声vk的协方差值;
为tk时刻参数向量估计值;
得到上述值后,进行参数向量粒子退化现象的判断。计算参数向量粒子的权重Neff2,如果小于权重阈值Nth2,则判断发生参数向量粒子退化现象。
步骤③:判断k+1是否存在,如不存在,则停止算法,如果存在,则继续进行步骤④更新状态噪声和参数噪声。
换而言之,步骤③判断k+1是否小于等于设定的IAPE最大计算次数或最大微观采样次数,若小于等于IAPE最大计算次数或最大微观采样次数,则继续进行步骤④,此外停止IAPF算法。此步骤实现IAPF算法在微观采样次数达到设定的IAPE最大计算次数或最大微观采样次数时停止。
步骤④:状态噪声协方差和参数噪声协方差更新-σx,k,σθ,l,为下一时刻的状态更新和参数更新提供***状态噪声和参数噪声。
状态向量新息矩阵:
状态噪声协方差更新:
参数向量新息矩阵:
参数噪声协方差更新:
ex,k为第tk时刻的状态向量的信息矩阵;
σx,k为第tk时刻的状态噪声的协方差值,其上下限参见公式(16);
eθ,l为第tk时刻的参数向量的信息矩阵;
σθ,k为第tk时刻的参数噪声的协方差值,其上下限参见公式(18);
步骤⑤:在完成状态噪声协方差和参数噪声协方差更新后,若产生状态向量粒子退化现象,重新采样生成新的***状态向量粒子集现有技术中有较多成熟的重新采样方法,本发明采用随机重新采样方法生成新的粒子集,其中去除低权值的粒子,复制高权值的粒子,也就是重新采样生成更新后的***状态向量粒子集同时所有***状态向量粒子的权重值都为1/N。
若产生向量粒子退化现象,重新采样生成新的***参数向量粒子集现有技术中有较多成熟的重新采样方法,本发明采用随机重新采样方法生成新的粒子集,其中去除低权值的粒子,复制高权值的粒子,也就是重新采样生成更新后的***参数向量粒子集同时所有***参数向量粒子的权重值都为1/M。
经过上述五步之后,获得了tk时刻下***参数向量估计值与***状态向量估计值之后需要把k+1作为新的微观采样点,判断k+1是否能被L整除,如果能,则进行步骤①;否则直接输出当前时刻的参数,优选的输出参数状态向量估计值;
3、SOC与SOH提取
基于上述多时间尺度H∞滤波算法,得到实时的电池参数与状态通过式(19)提取出状态量sk,参数量Ca,l、R0,l与Rp,l。
Ca,l表示tl×L时刻下更新的电池的容量值;R0,l与Rp,l分别表示tl×L时刻下更新的电池欧姆内阻与极化内组值。
式中,状态量sk即为***状态估计值修正是实时的荷电状态(SOC);参数量Ca,l、R0,l和Rp,l与***参数估计值修正相关,则能直接实时地反映电池的健康状态(SOH)。
在本发明中,电池的荷电状态SOC通过微观时间尺度实时更新,而健康状态我们选择用电池的容量来衡量,因为对于锂电池来说,电池容量越小,反应电池老化越严重,同时意味着电池健康状态(SOH)越差,在此算法过程中,电池容量精度较高,可以以此作为主要的SOH衡量参数;同时,考虑到电池内阻估计精度未得到充分检验,因而仅将其作为SOH的辅助衡量参数。
下面本发明通过选用某一型号镍钴锰NMC三元锂离子电池为例进行试验,得到基于多时间尺度的荷电状态和健康状态的估计值。
镍锰钴NMC三元锂离子电池的额定容量为2.5Ah,充放电截止电压分别为4.15V、2.75V。准备试验包括三个固定温度点(10℃、25℃、40℃)下的基础容量、开路电压、UDDS循环工况三项试验,以及室温条件下的基础容量与DST循环工况试验。建立容量-SOC-OCV三维响应面如图4所示。
通过上述多时间尺度IAPF算法来实现SOC与容量的联合估计。具体过程为:
首先,完成联合估计算法程序的调试。基于相应的UDDS试验数据,完成上述基于多时间尺度IAPF算法的SOC与容量联合估计算法程序的调试。
然后将室温下UDDS试验数据直接调入上述调试好的联合估计算法程序中,为了节省计算时间,我们选取SOC范围为80%到20%。将算法中SOC初值设置为70%、容量初值设置为2.1Ah(准确初值为2.096Ah),这里取尺度转换限值L=1s(微观时间尺度)和L=50s(宏观时间尺度)完成动力电池SOC与容量估计.
图5为UDDS工况的电流图和SOC图,其中的SOC是通过安时积分法计算的,作为多时间尺度IAPF算法SOC估计结果的参考值。
多时间尺度IAPF算法估计结果分别如图6、图7所示。图6为多时间尺度IAPF算法相应的电压和SOC结果,图6(a)为端电压的实测值和估计值对比图,端电压的误差如图6(b)所示。图6(c)为SOC参考值和估计值的对比图,SOC估计误差如图6(d)所示。从图6可以看出端电压绝对误差最大值小于0.05V,SOC除了在不精确的初值开始时,误差始终小于1.5%。图7为容量参考值和估计值的对比图,从图中可以看出容量估计值始终在参考值上下1%区间波动。从图6和图7可以看出,即使在SOC初值不准确的情况下,使用多时间尺度IAPF算法估计的SOC与容量精度都非常高,这是因为IAPF算法比传统的APF在估计过程中能够生产更合适的噪声协方差。
从上述分析得出,本发明所提出的多时间尺度IAPF算法进行***荷电状态与健康状态联合估计方法与传统方法相比具有以下优势:
(1)IAPF算法更适用于锂离子电池这种强非线性模型,估计精度更高,IAPF算法相比传统的APF算法可以有效提高SOC估计的收敛速度和鲁棒性;
(2)时间尺度IAPF算法可以在SOC初值均不准确的情况下均能快收敛到真值,即实现了电池最大容量未知时,SOC的准确估计,解决了传统SOC估计算法以最大可用容量已知为前提而无法成功运用到实车上的难题;
(3)容量与内阻都是衡量电池健康状态(SOH)的重要指标,因而上述联合估计算法在一定程度上实现了SOC与SOH的联合估计;
(4)容量-SOC-OCV三维响应面能够更好的实现电池在不同使用环境下的SOC和容量的联合估计。
Claims (9)
1.一种动力电池状态和参数估计的多时间尺度IAPF滤波方法,其特征在于
每个微观采样点进行状态估计,每隔L个微观采样点为一个宏观采样点并进行参数估计;
初始化:初始化所述IAPF滤波的宏观参数观测器和微观状态观测器的初始参数;
步骤①:状态估计:利用当前微观采样点k下动力电池的电流值、上个微观采样点k-1下的状态粒子集和上个宏观采样点l-1的参数估计值进行状态粒子集值预估,然后计算状态权重值,最后对状态权重值进行归一化处理,得到当前微观采样点k下的状态估计值;
步骤②:参数估计:通过更新后的动力电池开路电压计算参数粒子集值预估、参数权重值,然后对参数权重值进行归一化处理,得到当前宏观采样点l下的参数估计值。
2.如权利要求1所述的方法,其特征在于:还包括:
所述初始参数包括状态噪声和参数噪声。
3.如权利要求2所述的方法,其特征在于:还包括:
步骤③:判断k+1是否超过最大计算次数,超过,则停止估计,如果未超过,则进行步骤④;
步骤④:通过计算状态噪声和参数噪声的协方差,更新所述状态噪声和所述参数噪声。
4.如权利要求3所述的方法,其特征在于:还包括:
步骤⑤:针对状态粒子和参数粒子的退化现象,进行重新采样更新状态粒子集和参数粒子集;
经过上述五步之后,把k+1作为新的当前微观采样点状态,判断k+1是否能被L整除,如果能,则进行步骤①;否则直接输出参数估计值。
5.如权利要求1-4任意一项所述的方法,其特征在于:状态为极化电压或荷电状态,参数包括容量、极化电容、极化电阻或欧姆内阻。
6.如权利要求1-4任意一项所述的方法,其特征在于:步骤②中利用上个宏观采样点下的参数估计值、当前微观采样点k下的状态估计值和根据容量建立的动力电池的容量-荷电状态-开路电压三维响应面更新所述动力电池开路电压。
7.如权利要求4所述的方法,其特征在于:所述步骤⑤中,利用随机重采样方法更新状态粒子集和参数粒子集。
8.一种估计动力电池荷电状态与健康状态的方法,其特征在于:
首先,建立动力电池的容量-荷电状态-开路电压三维响应面;
其次,数据在线获取,实时采集动力电池单体或动力电池组的电压值和电流值;
然后,利用如权利要求1-7任一项所述的多时间尺度IAPF滤波方法,获得动力电池所述状态估计值与所述参数估计值;
最后,在线荷电状态与健康状态提取,利用所述状态估计值与参数估计值,估计动力电池荷电状态与健康状态。
9.如权利要求8所述的方法,其特征在于:所述状态是所述动力电池的荷电状态;利用参数估计值中的容量估计值来评估健康状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611077392.7A CN106772067B (zh) | 2016-11-30 | 2016-11-30 | 多时间尺度估计动力电池荷电状态与健康状态的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611077392.7A CN106772067B (zh) | 2016-11-30 | 2016-11-30 | 多时间尺度估计动力电池荷电状态与健康状态的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106772067A CN106772067A (zh) | 2017-05-31 |
CN106772067B true CN106772067B (zh) | 2018-12-28 |
Family
ID=58898739
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611077392.7A Active CN106772067B (zh) | 2016-11-30 | 2016-11-30 | 多时间尺度估计动力电池荷电状态与健康状态的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772067B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108414947B (zh) * | 2018-06-06 | 2021-05-28 | 哈尔滨工业大学 | 一种基于多时间尺度的空间锂离子电池状态联合估计方法 |
CN110018678B (zh) * | 2019-02-14 | 2021-05-14 | 电子科技大学 | 一种网联汽车控制***的故障诊断方法 |
CN110221221A (zh) * | 2019-04-24 | 2019-09-10 | 吉林大学 | 锂离子电池荷电状态和健康状态联合估计方法 |
CN110275118B (zh) * | 2019-06-27 | 2021-06-22 | 金龙联合汽车工业(苏州)有限公司 | 一种功率型动力电池健康状态估算方法 |
CN110888058B (zh) * | 2019-12-02 | 2021-12-31 | 西安科技大学 | 一种基于动力电池soc和soh联合估计的算法 |
CN112433104A (zh) * | 2020-10-30 | 2021-03-02 | 国网山东省电力公司青岛供电公司 | 一种智能变电站二次电压核相装置及方法 |
CN113030752B (zh) * | 2021-04-12 | 2024-03-29 | 安徽理工大学 | 一种基于变遗忘因子在线参数辨识和soc联合估计方法 |
CN113504481B (zh) * | 2021-06-07 | 2022-07-26 | 华南理工大学 | 一种电动汽车动力电池多状态联合估计***及方法 |
CN114089076A (zh) * | 2021-11-29 | 2022-02-25 | 国网北京市电力公司 | 一种动力电池的模拟电路 |
CN117686936B (zh) * | 2024-02-01 | 2024-05-03 | 杭州高特电子设备股份有限公司 | 电池***剩余使用寿命的在线检测方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103995464A (zh) * | 2014-05-26 | 2014-08-20 | 北京理工大学 | 一种估计电动车辆的动力***的参数和状态的方法 |
CN105319515A (zh) * | 2015-11-18 | 2016-02-10 | 吉林大学 | 锂离子电池荷电状态和健康状态联合估算方法 |
CN105842627A (zh) * | 2016-02-01 | 2016-08-10 | 北京理工大学 | 基于数据模型融合的动力电池容量和荷电状态的估计方法 |
-
2016
- 2016-11-30 CN CN201611077392.7A patent/CN106772067B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103995464A (zh) * | 2014-05-26 | 2014-08-20 | 北京理工大学 | 一种估计电动车辆的动力***的参数和状态的方法 |
CN105319515A (zh) * | 2015-11-18 | 2016-02-10 | 吉林大学 | 锂离子电池荷电状态和健康状态联合估算方法 |
CN105842627A (zh) * | 2016-02-01 | 2016-08-10 | 北京理工大学 | 基于数据模型融合的动力电池容量和荷电状态的估计方法 |
Non-Patent Citations (1)
Title |
---|
Adaptive Unscented Particle Filter Based on Predicted Residual;Hua-jian Wang 等;《2011 6th IEEE Joint International Information Technology and Artificial Intelligence Conference》;20111231;第2卷;第181-184页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106772067A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772067B (zh) | 多时间尺度估计动力电池荷电状态与健康状态的方法 | |
CN106291381B (zh) | 一种联合估计动力电池***荷电状态与健康状态的方法 | |
CN108196200B (zh) | 一种锂电池健康和荷电状态的联合模拟评估方法 | |
CN110261779B (zh) | 一种三元锂电池荷电状态与健康状态在线协同估计方法 | |
CN106443474B (zh) | 一种动力电池***寿命衰退特征快速识别的方法 | |
CN105807226B (zh) | 基于等效电路模型的锂离子电池soc预测方法及装置 | |
CN110361652A (zh) | 一种基于模型参数优化的卡尔曼滤波锂电池soc估计方法 | |
Seo et al. | Innovative lumped-battery model for state of charge estimation of lithium-ion batteries under various ambient temperatures | |
CN105510829B (zh) | 一种新型锂离子动力电池soc估计方法 | |
CN106093783A (zh) | 卡尔曼滤波与数据驱动融合的电池soc估计方法 | |
CN109870659A (zh) | 应用滑窗寻优策略的锂离子电池健康状态估算方法 | |
CN109164392A (zh) | 一种动力电池的soc估算方法 | |
CN110687462B (zh) | 一种动力电池soc与容量全生命周期联合估计方法 | |
CN106093782A (zh) | 动态建模的最小二乘支持向量机soc估计方法 | |
CN112557906A (zh) | 一种动力电池全生命周期内soc和容量在线联合估计方法 | |
CN106772104A (zh) | 一种动力电池soc值估算方法 | |
CN106291393A (zh) | 一种用于在线识别电池模型参数的方法 | |
CN109239602A (zh) | 一种动力电池的欧姆内阻的估算方法 | |
CN109490790B (zh) | 采用补偿脉冲法的锂动力电池功率特性测试方法及装置 | |
CN108829911A (zh) | 一种开路电压与soc函数关系优化方法 | |
CN109901072A (zh) | 基于历史数据与实验室测试数据的退役电池参数检测方法 | |
CN110196395A (zh) | 蓄电池soc估算方法 | |
CN106443496A (zh) | 一种带改进型噪声估计器的电池荷电状态估计方法 | |
CN115158076A (zh) | 计量误差评估方法、装置及计算机可读存储介质 | |
Chen et al. | A novel sliding mode observer for state of charge estimation of EV lithium batteries |
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 |