CN115598535A - Battery state of charge estimation method and system considering state of health - Google Patents
Battery state of charge estimation method and system considering state of health Download PDFInfo
- Publication number
- CN115598535A CN115598535A CN202211214705.4A CN202211214705A CN115598535A CN 115598535 A CN115598535 A CN 115598535A CN 202211214705 A CN202211214705 A CN 202211214705A CN 115598535 A CN115598535 A CN 115598535A
- Authority
- CN
- China
- Prior art keywords
- battery
- state
- charge
- health
- open circuit
- 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
Links
Images
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/367—Software therefor, e.g. for battery testing using modelling or look-up tables
-
- 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/378—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC] specially adapted for the type of battery or accumulator
-
- 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/385—Arrangements for measuring battery or accumulator variables
- G01R31/387—Determining ampere-hour charge capacity or SoC
- G01R31/388—Determining ampere-hour charge capacity or SoC involving voltage measurements
-
- 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)
- Secondary Cells (AREA)
- Tests Of Electric Status Of Batteries (AREA)
Abstract
本发明公开了一种考虑健康状态的电池荷电状态估计方法及系统,属于电池荷电状态估计技术领域,该方法包括步骤:获取电池充放电数据,构建数据集;建立相应电池模型;基于数据集和电池模型,估计电池的开路电压;将估计的开路电压和电池健康状态作为神经网络模型的输入,初步估计得到SOC;结合电流信息通过滤波算法得到SOC的最终估计值。本发明结合电池模型方法和数据驱动方法,实现对不同健康状态下电池荷电状态的准确估计。
The invention discloses a method and system for estimating the battery state of charge considering the state of health, and belongs to the technical field of battery state of charge estimation. Set the battery model to estimate the open circuit voltage of the battery; use the estimated open circuit voltage and battery health state as the input of the neural network model, and obtain the SOC through preliminary estimation; combine the current information to obtain the final estimated value of SOC through the filtering algorithm. The invention combines a battery model method and a data-driven method to realize accurate estimation of the state of charge of the battery under different health states.
Description
技术领域technical field
本发明属于电池荷电状态估计技术领域,具体涉及一种考虑健康状态的电池荷电状态估计方法及系统。The invention belongs to the technical field of battery state-of-charge estimation, and in particular relates to a battery state-of-charge estimation method and system considering the state of health.
背景技术Background technique
电池荷电状态(State of charge,SOC)是用于表征电池短期电量变化情况的重要参数,一般定义为当前剩余电量与当前最大可用容量之比,SOC估计是电池管理系统的核心功能之一。准确的SOC估计可以提高电池能量利用率,防止电池发生过充电、过放电,保障电池在使用过程中的安全性和可靠性,延长电池使用寿命。然而,受电池本体非线性、时变性等特征和复杂工况环境的影响,难以获得高精度的SOC估计结果。The state of charge (SOC) of the battery is an important parameter used to characterize the battery's short-term power changes. It is generally defined as the ratio of the current remaining power to the current maximum available capacity. SOC estimation is one of the core functions of the battery management system. Accurate SOC estimation can improve battery energy utilization, prevent battery overcharge and overdischarge, ensure battery safety and reliability during use, and prolong battery life. However, it is difficult to obtain high-precision SOC estimation results due to the nonlinear and time-varying characteristics of the battery body and the complex working environment.
常用的电池SOC估计方法主要包括:安时积分法、开路电压法、等效模型法和数据驱动法等。其中,安时积分法的计算简洁高效,但存在明显的误差累计效应;开路电压法需要对电池进行超过半小时的静置,难以应用于实际;等效模型法通常采用电池模型加滤波算法观测器实现SOC估计,能够在估计精度和运算复杂度之间实现较好的平衡,但适用性和泛化能力有限;数据驱动法通常以电流电压信息作为输入建立机器学习或深度学习模型直接估计得到SOC,估计精度较高,但该方法具有开环发散风险,运算量较大难以在线运行,且建立的模型为黑箱模型,缺乏可解释性,同时数据驱动法估计SOC的可解释性尚未得到充分研究。Commonly used battery SOC estimation methods mainly include: ampere-hour integration method, open circuit voltage method, equivalent model method, and data-driven method. Among them, the calculation of the ampere-hour integration method is simple and efficient, but there is an obvious error accumulation effect; the open circuit voltage method requires the battery to stand still for more than half an hour, which is difficult to apply in practice; the equivalent model method usually uses a battery model plus a filtering algorithm to observe SOC estimation can be achieved by a device, which can achieve a good balance between estimation accuracy and computational complexity, but its applicability and generalization ability are limited; data-driven methods usually use current and voltage information as input to establish machine learning or deep learning models to directly estimate SOC, the estimation accuracy is high, but this method has the risk of open-loop divergence, the amount of calculation is large and it is difficult to run online, and the established model is a black box model, which lacks interpretability. At the same time, the interpretability of the data-driven method to estimate SOC has not been fully obtained. Research.
上述估计方法大都没有考虑出现电池容量衰减的情况下如何继续保证SOC精度,即随着电池健康状态的改变,SOC估计误差逐渐增大。Most of the above estimation methods do not consider how to continue to ensure the SOC accuracy in the case of battery capacity decay, that is, as the battery health status changes, the SOC estimation error gradually increases.
发明内容Contents of the invention
为了克服上述现有技术的缺点,本发明的目的在于提供一种考虑健康状态的电池荷电状态估计方法及系统,能够有效解决随着健康状态改变,SOC估计出现的误差增大、数据驱动法开环发散风险以及可解释性差的技术问题。In order to overcome the shortcomings of the above-mentioned prior art, the object of the present invention is to provide a method and system for estimating the battery state of charge considering the state of health, which can effectively solve the problem of increasing errors in SOC estimation as the state of health changes, and data-driven methods. Open-loop divergence risks and technical issues with poor explainability.
为了达到上述目的,本发明采用以下技术方案予以实现:In order to achieve the above object, the present invention adopts the following technical solutions to achieve:
本发明公开了一种考虑健康状态(SOH)的电池荷电状态(SOC)估计方法,包括以下步骤:The invention discloses a battery state of charge (SOC) estimation method considering the state of health (SOH), comprising the following steps:
步骤1:利用公开数据或通过充放电设备,获取电池充放电数据及其对应的健康状态,构建数据集;Step 1: Use public data or charge and discharge equipment to obtain battery charge and discharge data and their corresponding health status, and build a data set;
步骤2:建立相应的电池模型;Step 2: Establish the corresponding battery model;
步骤3:基于步骤1的数据集和步骤2的电池模型估计电池的开路电压;Step 3: Estimate the open circuit voltage of the battery based on the data set in step 1 and the battery model in
步骤4:构建以开路电压(OCV)和健康状态作为输入、电池荷电状态作为输出的神经网络模型;Step 4: Construct a neural network model with open circuit voltage (OCV) and state of health as input and battery state of charge as output;
步骤5:结合步骤1的数据集中的电流信息和步骤3的神经网络模型输出的电池荷电状态估计初值,通过滤波算法进行融合,得到最终电池荷电状态估计值。Step 5: Combining the current information in the data set in step 1 and the estimated initial value of the battery state of charge output by the neural network model in step 3, perform fusion through a filtering algorithm to obtain the final estimated value of the battery state of charge.
优选地,步骤1中是利用公开数据或通过充放电设备获取电池充放电数据及对应的健康状态。Preferably, in step 1, the battery charging and discharging data and corresponding health status are obtained by using public data or through charging and discharging equipment.
优选地,充放电设备为电池检测系统或充电桩。Preferably, the charging and discharging device is a battery detection system or a charging pile.
优选地,所述电池为铅酸电池、镍氢电池、镍铬电池或锂离子电池。Preferably, the battery is a lead-acid battery, a nickel-metal hydride battery, a nickel-chromium battery or a lithium-ion battery.
优选地,步骤2中,构建电池模型,包括以下步骤:Preferably, in
步骤21、基于电池的一阶RC等效电路模型建立状态空间方程:Step 21, establish a state space equation based on the first-order RC equivalent circuit model of the battery:
其中,E0为开路电压,R1为欧姆内阻,R2为极化内阻,C2为极化电容,U2为极化电压,I为电流,U0为电池端电压;Among them, E 0 is the open circuit voltage, R 1 is the ohmic internal resistance, R 2 is the polarization internal resistance, C 2 is the polarization capacitance, U 2 is the polarization voltage, I is the current, and U 0 is the battery terminal voltage;
步骤22、计算状态空间方程的传递函数,并增加零阶保持器环节,计算传递函数G(s):Step 22. Calculate the transfer function of the state space equation, and add a zero-order retainer link to calculate the transfer function G(s):
其中,T为采样间隔,H(s)为零阶保持器的传递函数,s为一个变量符号;e为一个无理数;Among them, T is the sampling interval, H(s) is the transfer function of the zero-order holder, s is a variable symbol; e is an irrational number;
步骤23、对传递函数G(s)做反拉氏变换,并进行离散化处理:Step 23, perform an inverse Laplace transform on the transfer function G(s), and perform discretization processing:
求得k时刻的差分方程:Find the difference equation at time k:
步骤24、结合公式(1),计算得到k时刻的开路电压:Step 24, combined with formula (1), calculate the open circuit voltage at time k:
其中, in,
步骤25、令Vk=E0(k)-U0(k),建立电池的受控自回归滑动平均模型,即CARMA模型:Step 25. Set V k =E 0(k) -U 0(k) to establish a controlled autoregressive moving average model of the battery, namely the CARMA model:
其中,θ1=a,θ2=[R1+R2(1/a-1)],θ3=aR0,q-1为一步平移算子,ξk为系统噪声。Wherein, θ 1 =a, θ 2 =[R 1 +R 2 (1/a-1)], θ 3 =aR 0 , q -1 is a one-step translation operator, and ξ k is system noise.
优选地,步骤3中,采用带遗忘因子的递推最小二乘法对电池CARMA模型参数进行在线估计,包括以下步骤:Preferably, in step 3, the battery CARMA model parameters are estimated online using the recursive least squares method with a forgetting factor, including the following steps:
步骤31、由于开路电压在极短时间内变化量极小,故假设在一个采样间隔中开路电压不变,即E0(k)=E0(k-1),依据电池CARMA模型,得到:Step 31. Since the open circuit voltage changes very little in a very short time, it is assumed that the open circuit voltage does not change in a sampling interval, that is, E 0(k) =E 0(k-1) , and according to the battery CARMA model, it is obtained:
U0(k)=[U0(k-1) Ik -Ik-1 1][θ1k θ2k θ3k θ4k]T (26)U 0(k) =[U 0(k-1) I k -I k-1 1][θ 1k θ 2k θ 3k θ 4k ] T (26)
其中,θ4k=(1-a)E0(k),由此求得开路电压E0(k)的估计值:Among them, θ 4k =(1-a)E 0(k) , thus obtain the estimated value of the open circuit voltage E 0(k) :
步骤32、采用带遗忘因子的递推最小二乘法进行开路电压的在线估计,初始化估计值:Step 32. Use the recursive least squares method with forgetting factor to estimate the open circuit voltage online, and initialize the estimated value:
其中,x为待估计参数,x0为待估计参数的初始值,P0为误差的协方差矩阵初始值,E为单位矩阵;Among them, x is the parameter to be estimated, x 0 is the initial value of the parameter to be estimated, P 0 is the initial value of the covariance matrix of the error, and E is the identity matrix;
步骤33、进行迭代求解,并依据公式(9)估计得到开路电压,迭代过程如下:Step 33, perform an iterative solution, and estimate the open circuit voltage according to formula (9), the iterative process is as follows:
ek=yk-Hkxk-1 (29)e k =y k -H k x k-1 (29)
xk=xk-1+Kkek (31)x k =x k-1 +K k e k (31)
其中,k均指代k时刻,yk为电池端电压,ek为yk的预计误差,Hk=[U0(k-1) Ik -Ik-11],xk=[θ1k θ2k θ3k θ4k]T为待估计参数阵,Kk为增益系数,Pk为协方差矩,λ为遗忘因子。Among them, k refers to time k, y k is the battery terminal voltage, e k is the expected error of y k , H k = [U 0(k-1) I k -I k-1 1], x k = [ θ 1k θ 2k θ 3k θ 4k ] T is the parameter array to be estimated, K k is the gain coefficient, P k is the covariance moment, and λ is the forgetting factor.
优选地,步骤4中,构建神经网络模型的具体方法,包括:Preferably, in step 4, the concrete method of constructing neural network model includes:
构建包含输入层、隐含层和输出层的BP神经网络模型,其输入特征序列为x=[OCVSOH]T,即电池的开路电压和健康状态,输出序列为y=SOC,即电池的荷电状态,隐含层神经元个数设置为64,激活函数为ReLU函数,整个网络对于任一样本的输出为:Construct a BP neural network model including an input layer, a hidden layer and an output layer. The input feature sequence is x=[OCVSOH] T , which is the open circuit voltage and state of health of the battery, and the output sequence is y=SOC, which is the charge of the battery state, the number of neurons in the hidden layer is set to 64, and the activation function is the ReLU function. The output of the entire network for any sample is:
其中,βi为隐含层第i个神经元与输出层之间的连接权值,g(x)为激活函数,ai为输入层与隐含层连接权重矩阵的第i行,xj为第j个样本的输入矩阵,bi为输入层与隐含层连接偏置矩阵的第i行,bi i为输出层偏置值,P为样本个数。Among them, β i is the connection weight between the i-th neuron in the hidden layer and the output layer, g(x) is the activation function, a i is the i-th row of the connection weight matrix between the input layer and the hidden layer, x j is the input matrix of the jth sample, b i is the ith row of the bias matrix connecting the input layer and the hidden layer, b i i is the output layer bias value, and P is the number of samples.
优选地,步骤5中,通过滤波算法进行融合,具体包括:Preferably, in step 5, fusion is performed through a filtering algorithm, specifically including:
步骤51、计算电流与电池荷电状态之间的关系:Step 51, calculate the relationship between current and battery state of charge:
其中,Q为电池当前最大可用容量,T为采样间隔,I为电流;Among them, Q is the current maximum available capacity of the battery, T is the sampling interval, and I is the current;
步骤52、结合BP神经网络估计的SOC结果和公式(16),得到状态空间方程的递推式:Step 52, in conjunction with the SOC result estimated by the BP neural network and formula (16), obtain the recursive formula of the state space equation:
其中,SOCk为k时刻的SOC估计值,SOCBP,k为通过BP神经网络模型得到的k时刻SOC估计值,Ik为k时刻的电流;Wherein, SOC k is the SOC estimated value of k moment, SOC BP,k is the SOC estimated value of k moment obtained by BP neural network model, I k is the electric current of k moment;
步骤53、化简求得系统的状态方程和观测方程,并采用卡尔曼滤波算法,进行时间更新及测量更新,得到最终电池荷电状态估计值:Step 53. Simplify and obtain the state equation and observation equation of the system, and use the Kalman filter algorithm to perform time update and measurement update to obtain the final estimated value of the battery state of charge:
其中,xk为系统状态变量,即SOC,yk为输出量,即BP神经网络模型估计得到的SOC,uk为输入量,A为状态矩阵,B为输入矩阵,C为输出矩阵,ωk为过程噪声,vk为观测噪声。Among them, x k is the system state variable, that is, SOC, y k is the output, that is, the SOC estimated by the BP neural network model, u k is the input, A is the state matrix, B is the input matrix, C is the output matrix, ω k is the process noise, and v k is the observation noise.
本发明还公开了一种考虑健康状态的电池荷电状态估计系统,包括:The present invention also discloses a battery state of charge estimation system considering the state of health, including:
数据集构建模块,用于获取电池充放电数据及对应的健康状态;The data set building block is used to obtain battery charge and discharge data and corresponding health status;
电池模型构建模块,用于构建电池模型;Battery model building blocks for building battery models;
开路电压估计模块,用于根据构建的数据集和电池模型估计电池的开路电压;The open circuit voltage estimation module is used to estimate the open circuit voltage of the battery according to the constructed data set and battery model;
神经网络模型构建模块,用于构建以开路电压和充放电数据健康状态作为输入、以电池荷电状态作为输出的神经网络模型;The neural network model building block is used to construct a neural network model that takes the open circuit voltage and the state of health of the charge and discharge data as input, and takes the state of charge of the battery as the output;
电池荷电状态估计模块,用于结合数据集中的电流信息和神经网络模型输出的电池荷电状态估计初值,通过滤波算法进行融合,得到最终电池荷电状态估计值。The battery state of charge estimation module is used to combine the current information in the data set and the estimated initial value of the battery state of charge output by the neural network model, and fuse them through a filtering algorithm to obtain the final estimated value of the battery state of charge.
与现有技术相比,本发明具有以下有益效果:Compared with the prior art, the present invention has the following beneficial effects:
本发明所提出的估计方法通过电压、电流数据和电池模型估计电池的开路电压,并将其作为神经网络的输入量,赋予神经网络领域知识,相较于传统的电压、电流输入量,本方法具有一定的可解释性。本发明采用滤波算法融合估计结果,可以避免数据驱动法带来的开环发散风险,使得方法具有较高的鲁棒性。本发明考虑电池的健康状态,将健康状态同样作为神经网络的输入量,实现在不同的健康状态下均能准确估计电池SOC,对于提升电池管理系统的性能具有重要意义。The estimation method proposed by the present invention estimates the open circuit voltage of the battery through voltage, current data and battery model, and uses it as the input of the neural network to endow the neural network with domain knowledge. Compared with the traditional voltage and current input, this method have a certain degree of interpretability. The invention adopts a filter algorithm to fuse estimation results, which can avoid the open-loop divergence risk brought by the data-driven method, so that the method has higher robustness. The invention considers the health state of the battery, and uses the health state as the input of the neural network to realize accurate estimation of the battery SOC under different health states, which is of great significance for improving the performance of the battery management system.
附图说明Description of drawings
图1为本发明提出的方法流程图。Fig. 1 is a flow chart of the method proposed by the present invention.
图2为本发明实施例中一阶RC等效电路模型。FIG. 2 is a first-order RC equivalent circuit model in an embodiment of the present invention.
图3为本发明实施例构建的BP神经网络结构图。Fig. 3 is a structural diagram of the BP neural network constructed by the embodiment of the present invention.
图4为本发明实施例中电池SOH为97.5%时的荷电状态估计结果。Fig. 4 is the estimation result of the state of charge of the battery when the SOH of the battery is 97.5% in the embodiment of the present invention.
图5为本发明实施例中电池SOH为73.2%时的荷电状态估计结果。Fig. 5 is the estimation result of the state of charge of the battery when the SOH of the battery is 73.2% in the embodiment of the present invention.
具体实施方式detailed description
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。In order to enable those skilled in the art to better understand the solutions of the present invention, the following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only It is an embodiment of a part of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts shall fall within the protection scope of the present invention.
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。It should be noted that the terms "first" and "second" in the description and claims of the present invention and the above drawings are used to distinguish similar objects, but not necessarily used to describe a specific sequence or sequence. It is to be understood that the data so used are interchangeable under appropriate circumstances such that the embodiments of the invention described herein can be practiced in sequences other than those illustrated or described herein. Furthermore, the terms "comprising" and "having", as well as any variations thereof, are intended to cover a non-exclusive inclusion, for example, a process, method, system, product or device comprising a sequence of steps or elements is not necessarily limited to the expressly listed Those steps or elements may instead include other steps or elements not explicitly listed or inherent to the process, method, product or apparatus.
下面结合附图对本发明做进一步详细描述:The present invention is described in further detail below in conjunction with accompanying drawing:
参见图1,本发明提供的一种考虑健康状态的电池荷电状态估计方法,包括以下步骤:Referring to Fig. 1, a method for estimating the state of charge of a battery considering the state of health provided by the present invention includes the following steps:
第一步,利用公开数据或通过充放电设备(电池检测系统、充电桩等),获取电池充放电数据及其对应的健康状态(SOH),构建数据集;The first step is to use public data or charge and discharge equipment (battery detection system, charging pile, etc.) to obtain battery charge and discharge data and its corresponding state of health (SOH), and build a data set;
第二步,建立相应的电池模型。具体建立过程,包括以下步骤:The second step is to establish the corresponding battery model. The specific establishment process includes the following steps:
步骤21、基于电池的一阶RC等效电路模型建立状态空间方程:Step 21, establish a state space equation based on the first-order RC equivalent circuit model of the battery:
其中,E0为开路电压,R1为欧姆内阻,R2为极化内阻,C2为极化电容,U2为极化电压,I为电流,U0为电池端电压;Among them, E 0 is the open circuit voltage, R 1 is the ohmic internal resistance, R 2 is the polarization internal resistance, C 2 is the polarization capacitance, U 2 is the polarization voltage, I is the current, and U 0 is the battery terminal voltage;
步骤22、计算状态空间方程的传递函数,并增加零阶保持器环节,计算传递函数G(s):Step 22. Calculate the transfer function of the state space equation, and add a zero-order retainer link to calculate the transfer function G(s):
其中,T为采样间隔,H(s)为零阶保持器的传递函数,s为一个变量符号;e为一个无理数;Among them, T is the sampling interval, H(s) is the transfer function of the zero-order holder, s is a variable symbol; e is an irrational number;
步骤23、对传递函数G(s)做反拉氏变换,并进行离散化处理:Step 23, perform an inverse Laplace transform on the transfer function G(s), and perform discretization processing:
求得k时刻的差分方程:Find the difference equation at time k:
步骤24、结合公式(1),计算得到k时刻的开路电压:Step 24, combined with formula (1), calculate the open circuit voltage at time k:
其中, in,
步骤25、令Vk=E0(k)-U0(k),建立电池的受控自回归滑动平均模型,即CARMA模型:Step 25. Set V k =E 0(k) -U 0(k) to establish a controlled autoregressive moving average model of the battery, namely the CARMA model:
其中,θ1=a,θ2=[R1+R2(1/a-1)],θ3=aR0,q-1为一步平移算子,ξk为系统噪声。Wherein, θ 1 =a, θ 2 =[R 1 +R 2 (1/a-1)], θ 3 =aR 0 , q -1 is a one-step translation operator, and ξ k is system noise.
第三步,基于数据集和电池模型,估计电池的开路电压(OCV)。具体采用带遗忘因子的递推最小二乘法对电池CARMA模型参数进行在线估计,包括以下步骤:In the third step, based on the dataset and the battery model, the open circuit voltage (OCV) of the battery is estimated. Specifically, the recursive least squares method with forgetting factor is used to estimate the parameters of the battery CARMA model online, including the following steps:
步骤31、由于开路电压在极短时间内变化量极小,故假设在一个采样间隔中开路电压不变,即E0(k)=E0(k-1),依据电池CARMA模型,得到:Step 31. Since the open circuit voltage changes very little in a very short time, it is assumed that the open circuit voltage does not change in a sampling interval, that is, E 0(k) =E 0(k-1) , and according to the battery CARMA model, it is obtained:
U0(k)=[U0(k-1) Ik -Ik-1 1][θ1k θ2k θ3k θ4k]T (44)U 0(k) =[U 0(k-1) I k -I k-1 1][θ 1k θ 2k θ 3k θ 4k ] T (44)
其中,θ4k=(1-a)E0(k),由此求得开路电压E0(k)的估计值:Among them, θ 4k =(1-a)E 0(k) , thus obtain the estimated value of the open circuit voltage E 0(k) :
步骤32、采用带遗忘因子的递推最小二乘法进行开路电压的在线估计,初始化估计值:Step 32. Use the recursive least squares method with forgetting factor to estimate the open circuit voltage online, and initialize the estimated value:
其中,x为待估计参数,x0为待估计参数的初始值,P0为误差的协方差矩阵初始值,E为单位矩阵;Among them, x is the parameter to be estimated, x 0 is the initial value of the parameter to be estimated, P 0 is the initial value of the covariance matrix of the error, and E is the identity matrix;
步骤33、进行迭代求解,并依据公式(9)估计得到开路电压,迭代过程如下:Step 33, perform an iterative solution, and estimate the open circuit voltage according to formula (9), the iterative process is as follows:
ek=yk-Hkxk-1 (47)e k =y k -H k x k-1 (47)
xk=xk-1+Kkek (49)x k =x k-1 +K k e k (49)
其中,k均指代k时刻,yk为电池端电压,ek为yk的预计误差,Hk=[U0(k-1) Ik -Ik-11],xk=[θ1k θ2k θ3k θ4k]T为待估计参数阵,Kk为增益系数,Pk为协方差矩,λ为遗忘因子。Among them, k refers to time k, y k is the battery terminal voltage, e k is the expected error of y k , H k = [U 0(k-1) I k -I k-1 1], x k = [ θ 1k θ 2k θ 3k θ 4k ] T is the parameter array to be estimated, K k is the gain coefficient, P k is the covariance moment, and λ is the forgetting factor.
第四步,构建以开路电压和健康状态作为输入,SOC作为输出的神经网络模型。构建神经网络模型的具体方法,具体包括:The fourth step is to build a neural network model with open circuit voltage and health status as input and SOC as output. The specific method of constructing the neural network model includes:
构建包含输入层、隐含层和输出层的BP神经网络模型,其输入特征序列为x=[OCVSOH]T,即电池的开路电压和健康状态,输出序列为y=SOC,即电池的荷电状态,隐含层神经元个数设置为64,激活函数为ReLU函数,整个网络对于任一样本的输出为:Construct a BP neural network model including an input layer, a hidden layer and an output layer. The input feature sequence is x=[OCVSOH] T , which is the open circuit voltage and state of health of the battery, and the output sequence is y=SOC, which is the charge of the battery state, the number of neurons in the hidden layer is set to 64, and the activation function is the ReLU function. The output of the entire network for any sample is:
其中,βi为隐含层第i个神经元与输出层之间的连接权值,g(x)为激活函数,ai为输入层与隐含层连接权重矩阵的第i行,xj为第j个样本的输入矩阵,bi为输入层与隐含层连接偏置矩阵的第i行,bi i为输出层偏置值,P为样本个数。Among them, β i is the connection weight between the i-th neuron in the hidden layer and the output layer, g(x) is the activation function, a i is the i-th row of the connection weight matrix between the input layer and the hidden layer, x j is the input matrix of the jth sample, b i is the ith row of the bias matrix connecting the input layer and the hidden layer, b i i is the output layer bias value, and P is the number of samples.
第五步,结合数据集中的电流信息和神经网络模型输出的SOC估计初值,通过滤波算法进行融合,得到最终准确的SOC估计值。The fifth step is to combine the current information in the data set with the estimated initial value of SOC output by the neural network model, and fuse it through a filtering algorithm to obtain the final and accurate estimated value of SOC.
通过滤波算法进行融合,具体包括:Fusion is performed through filtering algorithms, including:
步骤51、计算电流与电池荷电状态之间的关系:Step 51, calculate the relationship between current and battery state of charge:
其中,Q为电池当前最大可用容量,T为采样间隔,I为电流;Among them, Q is the current maximum available capacity of the battery, T is the sampling interval, and I is the current;
步骤52、结合BP神经网络估计的SOC结果和公式(16),得到状态空间方程的递推式:Step 52, in conjunction with the SOC result estimated by the BP neural network and formula (16), obtain the recursive formula of the state space equation:
其中,SOCk为k时刻的SOC估计值,SOCBP,k为通过BP神经网络模型得到的k时刻SOC估计值,Ik为k时刻的电流;Wherein, SOC k is the SOC estimated value of k moment, SOC BP,k is the SOC estimated value of k moment obtained by BP neural network model, I k is the electric current of k moment;
步骤53、化简求得系统的状态方程和观测方程,并采用卡尔曼滤波算法,进行时间更新及测量更新,得到最终电池荷电状态估计值:Step 53. Simplify and obtain the state equation and observation equation of the system, and use the Kalman filter algorithm to perform time update and measurement update to obtain the final estimated value of the battery state of charge:
其中,xk为系统状态变量,即SOC,yk为输出量,即BP神经网络模型估计得到的SOC,uk为输入量,A为状态矩阵,B为输入矩阵,C为输出矩阵,ωk为过程噪声,vk为观测噪声。Among them, x k is the system state variable, that is, SOC, y k is the output, that is, the SOC estimated by the BP neural network model, u k is the input, A is the state matrix, B is the input matrix, C is the output matrix, ω k is the process noise, and v k is the observation noise.
具体应用实施例:Specific application examples:
本实施例中的测试对象为标称容量10Ah的多氟多三元软包锂电池(DFDPSP1265132-10Ah),对其进行脉冲特性测试以及不同健康状态下的城市道路工况(UDDS)和新欧洲驾驶周期(NEDC)放电测试,获取测试过程中的电压、电流、时间及其对应的SOH,构建数据集。其中,采用不同老化程度下的UDDS工况放电数据进行建模,采用不同老化程度下的NEDC工况放电数据进行验证。The test object in this embodiment is a polyfluorine multi-ternary soft-pack lithium battery (DFDPSP1265132-10Ah) with a nominal capacity of 10Ah, and it is tested for pulse characteristics and urban road conditions (UDDS) and new European conditions under different health states. Driving cycle (NEDC) discharge test, obtain the voltage, current, time and their corresponding SOH during the test, and build a data set. Among them, the discharge data of UDDS operating conditions under different aging degrees are used for modeling, and the discharge data of NEDC operating conditions under different aging degrees are used for verification.
基于电池的一阶RC等效电路模型建立微分方程数学模型,参见图2,其中E0为开路电压,R1为欧姆内阻,R2为极化内阻,C2为极化电容,U2为极化电压,I为电流,U0为电池端电压。通过对数学模型进行离散与差分处理,建立电池的受控自回归滑动平均模型(CARMA模型),如下式所示:Based on the first-order RC equivalent circuit model of the battery, a differential equation mathematical model is established, see Figure 2 , where E0 is the open circuit voltage, R1 is the ohmic internal resistance, R2 is the polarization internal resistance, C2 is the polarization capacitance, U 2 is the polarization voltage, I is the current, and U 0 is the battery terminal voltage. By performing discrete and differential processing on the mathematical model, a controlled autoregressive moving average model (CARMA model) of the battery is established, as shown in the following formula:
式中:θ1为θ2为[R1+R2(1/a-1)];θ3为aR0;q-1为一步平移算子;ξk为系统噪声。In the formula: θ 1 is θ 2 is [R 1 +R 2 (1/a-1)]; θ 3 is aR 0 ; q -1 is a one-step translation operator; ξ k is system noise.
为避免“数据饱和”问题,在递推最小二乘法中引入遗忘因子,兼顾程序的跟踪性能和稳定后的误差需求。通过带遗忘因子的递推最小二乘法(FFRLS法)对电池CARMA模型参数进行在线估计,从而得到开路电压估计值。其中递推最小二乘法中的遗忘因子设置为0.9,协方差矩阵初始值P0为单位矩阵。In order to avoid the "data saturation" problem, the forgetting factor is introduced into the recursive least squares method, taking into account the tracking performance of the program and the error requirements after stabilization. The parameters of the battery CARMA model are estimated online by the recursive least squares method with forgetting factor (FFRLS method), so as to obtain the estimated value of the open circuit voltage. The forgetting factor in the recursive least squares method is set to 0.9, and the initial value P 0 of the covariance matrix is the identity matrix.
采用BP神经网络模型为主体实现基于OCV和SOH的SOC估计,即由电池的SOH和第三步得到的OCV作为BP神经网络的输入特征,电池SOC作为输出。BP神经网络模型分为输入层、隐含层和输出层三层结构,其结构参见图3。BP神经网络模型的隐含层神经元个数设置为64,激活函数为ReLU函数。The BP neural network model is used as the main body to realize the SOC estimation based on OCV and SOH, that is, the SOH of the battery and the OCV obtained in the third step are used as the input characteristics of the BP neural network, and the battery SOC is used as the output. The BP neural network model is divided into a three-layer structure of input layer, hidden layer and output layer, and its structure is shown in Figure 3. The number of neurons in the hidden layer of the BP neural network model is set to 64, and the activation function is the ReLU function.
由第四步可以初步实现SOC估计,但估计精度相对较低且波动较大,同时考虑到数据驱动法的开环发散风险,因此在第五步结合电流信息通过卡尔曼滤波算法进一步提升估计精度。采用最大绝对误差(MaxAE)和均方根误差(RMSE)作为估计误差量化指标。The fourth step can preliminarily realize the SOC estimation, but the estimation accuracy is relatively low and fluctuates greatly. At the same time, considering the open-loop divergence risk of the data-driven method, the fifth step combines the current information with the Kalman filter algorithm to further improve the estimation accuracy. . The maximum absolute error (MaxAE) and the root mean square error (RMSE) are used as the quantitative indicators of the estimated error.
本实施例中实现了不同SOH下的准确SOC估计,并与基于一阶RC模型的双自适应衰减卡尔曼滤波(DAFEKF)进行对比:本发明所提出的方法在电池SOH从97.7%到73.2%之间的范围内MaxAE均不超过1.5%,且RMSE在0.5%左右,而作为对比的DAFEKF方法随着SOH的降低,SOC估计误差逐渐增大,表明本发明所提出方法的优异性。In this embodiment, accurate SOC estimation under different SOHs is realized, and compared with the double-adaptive attenuation Kalman filter (DAFEKF) based on the first-order RC model: the method proposed by the present invention can improve the battery SOH from 97.7% to 73.2% In the range between, the MaxAE is not more than 1.5%, and the RMSE is about 0.5%, while the DAFEKF method as a comparison, with the decrease of SOH, the SOC estimation error gradually increases, indicating the superiority of the method proposed by the present invention.
图4和图5详细展示了两个不同SOH下的SOC估计结果与对比。在SOH为97.5%时,采用本发明提出的方法估计得到的SOC与采用DAFEKF方法得到的SOC误差接近,误差均不超过2%,都能较为准确的估计得到SOC值。但当电池SOH达到73.2%时,DAFEKF方法的估计误差已经超过6%,而本发明提出的SOC估计方法依然能够保证最大绝对误差不超过1.5%,且可以克服常见滤波类方法中的初值依赖问题,在计算初期也无明显的收敛过程。由此证明所提出SOC估计方法的估计精度在不同电池健康状态下都能得到保证。Figure 4 and Figure 5 show in detail the SOC estimation results and comparisons under two different SOHs. When the SOH is 97.5%, the error of the SOC estimated by the method proposed by the present invention is close to the SOC obtained by the DAFEKF method, and the error is no more than 2%, and the SOC value can be estimated more accurately. However, when the battery SOH reaches 73.2%, the estimation error of the DAFEKF method has exceeded 6%, while the SOC estimation method proposed by the present invention can still ensure that the maximum absolute error does not exceed 1.5%, and can overcome the initial value dependence in common filtering methods problem, there is no obvious convergence process at the initial stage of calculation. This proves that the estimation accuracy of the proposed SOC estimation method can be guaranteed under different battery health states.
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。The above content is only to illustrate the technical ideas of the present invention, and cannot limit the protection scope of the present invention. Any changes made on the basis of the technical solutions according to the technical ideas proposed in the present invention shall fall within the scope of the claims of the present invention. within the scope of protection.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211214705.4A CN115598535B (en) | 2022-09-30 | 2022-09-30 | A battery state of charge estimation method and system considering health status |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211214705.4A CN115598535B (en) | 2022-09-30 | 2022-09-30 | A battery state of charge estimation method and system considering health status |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115598535A true CN115598535A (en) | 2023-01-13 |
CN115598535B CN115598535B (en) | 2025-05-09 |
Family
ID=84844240
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211214705.4A Active CN115598535B (en) | 2022-09-30 | 2022-09-30 | A battery state of charge estimation method and system considering health status |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115598535B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116068449A (en) * | 2023-03-23 | 2023-05-05 | 中国科学技术大学 | Lithium battery evaluation system and method for physical model constraint deep neural network |
CN117686937A (en) * | 2024-02-02 | 2024-03-12 | 河南科技学院 | A method for estimating the health status of individual cells in a battery system |
CN119029350A (en) * | 2024-08-14 | 2024-11-26 | 东莞市广劲动力科技有限公司 | Battery control method and system |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301511A (en) * | 2015-11-12 | 2016-02-03 | 华晨汽车集团控股有限公司 | Method for utilizing joint extended Kalman filter to pre-estimate battery state of charge |
CN105842633A (en) * | 2016-05-30 | 2016-08-10 | 广西大学 | Method for estimating SOC (State of Charge) of lithium ion battery based on gray extended Kalman filtering algorithm |
US20190036356A1 (en) * | 2017-07-31 | 2019-01-31 | Robert Bosch Gmbh | Method and System for Estimating Battery Open Cell Voltage, State of Charge, and State of Health During Operation of the Battery |
-
2022
- 2022-09-30 CN CN202211214705.4A patent/CN115598535B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301511A (en) * | 2015-11-12 | 2016-02-03 | 华晨汽车集团控股有限公司 | Method for utilizing joint extended Kalman filter to pre-estimate battery state of charge |
CN105842633A (en) * | 2016-05-30 | 2016-08-10 | 广西大学 | Method for estimating SOC (State of Charge) of lithium ion battery based on gray extended Kalman filtering algorithm |
US20190036356A1 (en) * | 2017-07-31 | 2019-01-31 | Robert Bosch Gmbh | Method and System for Estimating Battery Open Cell Voltage, State of Charge, and State of Health During Operation of the Battery |
Non-Patent Citations (2)
Title |
---|
宁博 等: "采用等效电路的参数自适应电池模型及电池荷电状态估计方法", 西安交通大学学报, vol. 49, no. 10, 31 October 2015 (2015-10-31), pages 67 - 71 * |
黄玉莎 等: "基于CARMA模型的动力电池荷电状态估计", 控制与信息技术, no. 2, 31 December 2020 (2020-12-31), pages 36 - 40 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116068449A (en) * | 2023-03-23 | 2023-05-05 | 中国科学技术大学 | Lithium battery evaluation system and method for physical model constraint deep neural network |
CN117686937A (en) * | 2024-02-02 | 2024-03-12 | 河南科技学院 | A method for estimating the health status of individual cells in a battery system |
CN117686937B (en) * | 2024-02-02 | 2024-04-12 | 河南科技学院 | Method for estimating health state of single battery in battery system |
CN119029350A (en) * | 2024-08-14 | 2024-11-26 | 东莞市广劲动力科技有限公司 | Battery control method and system |
CN119029350B (en) * | 2024-08-14 | 2025-04-04 | 东莞市广劲动力科技有限公司 | Battery control method and system |
Also Published As
Publication number | Publication date |
---|---|
CN115598535B (en) | 2025-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106026260B (en) | A kind of series winding battery pack SOC estimation method with equalizing circuit | |
CN103197251B (en) | An identification method for the second-order RC equivalent model of power lithium battery | |
CN111581904A (en) | Lithium battery SOC and SOH collaborative estimation method considering influence of cycle number | |
CN115598535A (en) | Battery state of charge estimation method and system considering state of health | |
CN110261779A (en) | A kind of ternary lithium battery charge state cooperates with estimation method with health status online | |
CN113777510A (en) | Method and device for estimating state of charge of lithium battery | |
CN107576919A (en) | Power battery charged state estimating system and method based on ARMAX models | |
CN109061506A (en) | Lithium-ion-power cell SOC estimation method based on Neural Network Optimization EKF | |
CN110275113A (en) | A method for estimating the state of charge of a lithium battery | |
CN112083331A (en) | A fusion method to improve the estimation accuracy of lithium-ion battery state of charge | |
CN108445402A (en) | A kind of lithium-ion-power cell state-of-charge method of estimation and system | |
CN103323781B (en) | Power battery pack on-line parameter detection system and SOC method of estimation | |
CN106250576A (en) | A kind of modeling method of lithium battery model based on motional impedance | |
CN107037374A (en) | A kind of SOC and SOH combined estimation methods of lithium ion battery | |
CN112345939B (en) | Lithium ion battery model parameter identification method based on continuous impulse response | |
CN110554324A (en) | SOC and SOH joint estimation method | |
CN113128672B (en) | Lithium ion battery pack SOH estimation method based on transfer learning algorithm | |
CN106772067A (en) | The method that Multiple Time Scales IAPF filters estimated driving force battery charge state and health status | |
CN111400935B (en) | Lithium battery three-order model parameter identification method and system based on likelihood function | |
CN109917299B (en) | Three-layer filtering estimation method for state of charge of lithium battery | |
CN112946481A (en) | Based on federation H∞Filtering sliding-mode observer lithium ion battery SOC estimation method and battery management system | |
CN116224074A (en) | A method, device, and storage medium for estimating the state of charge of a soft-pack lithium-ion battery | |
CN113740735A (en) | Method for estimating SOC of lithium ion battery | |
CN114814591A (en) | A lithium battery SOE estimation method, device and system | |
CN115327389A (en) | SOC estimation method of lithium battery based on double Kalman filter improved by genetic algorithm |
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 |