CN107209227A - 电池组的电池单元的充电状态的自动估计方法 - Google Patents
电池组的电池单元的充电状态的自动估计方法 Download PDFInfo
- Publication number
- CN107209227A CN107209227A CN201580069203.9A CN201580069203A CN107209227A CN 107209227 A CN107209227 A CN 107209227A CN 201580069203 A CN201580069203 A CN 201580069203A CN 107209227 A CN107209227 A CN 107209227A
- Authority
- CN
- China
- Prior art keywords
- battery unit
- moment
- mrow
- capacity
- predicted value
- 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
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/396—Acquisition or processing of data for testing or for monitoring individual cells or groups of cells within a battery
-
- 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/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/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
- G01R31/3842—Arrangements for monitoring battery or accumulator variables, e.g. SoC combining voltage and current 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/389—Measuring internal impedance, internal conductance or related variables
Landscapes
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Computer Hardware Design (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Secondary Cells (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Tests Of Electric Status Of Batteries (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Power Engineering (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
Abstract
本发明涉及一种用于估计电池单元的充电状态的方法,其包括对所述电池单元的容量预测值Cn,k3进行修正,其中,所述修正包括:在时刻k3获取(176)由以下关系式限定的可测量物理量zk3:公式(I)其中k是最接近时刻k3的时刻,并且N是大于或等于一的整数,‑使用以下观测模型来计算(176)所述可测量物理量zk3的预测值zk3:公式(II)其中:‑SOCk‑N是该电池单元在更早时刻k‑N的充电状态,‑Cn,k3是就在修正之前在时刻k3计算的该电池单元容量的预测值,以及‑根据所获取的物理量zk3和所述预测值zk3之间的差来修正(178)容量预测值Cn,k3。
Description
本发明涉及一种电池组的电池单元的充电状态的自动估计方法以及用于实施该方法的存储介质和电池组管理***。本发明还涉及包括该电池组管理***的机动车辆。
已知的用于估计电池组的电池单元的充电状态的方法包括:
a)在时刻k获取电池单元的端子之间的电压的测量值yk和对该电池单元进行充电或放电的电流的测量强度ik;以及
b)根据测量值yk、测量强度ik和电池单元的容量Cn,k3来估计该电池单元在时刻k的充电状态SOCk,该容量Cn,k3以Ah表示,其代表在时刻k3能够存储在电池单元中的最大电能量,所述时刻k3是最接近时刻k的已经估计该容量Cn,k3的时刻;以及
c)估计时刻k3的容量Cn,k3,所述时刻k3不如时刻k频繁,所述估计包括:
·使用将容量Cn,k3与同一电池单元在前一时刻k3-1的容量Cn,k3-1相关联的状态模型来计算该容量Cn,k3的预测值,然后
·根据对电池单元进行充电或放电的电流在最接近时刻k3的时刻的测量强度来修正容量Cn,k3的预测值。
例如,这样的方法在以下文献的部分3的第4.2段中进行描述:L.Plett等:“Extended Kalman filtering for battery management systems of LiPB-based HEVbattery packs”,Journal of Power Sources,2004年,第252-292页。以下,该文献用缩写“Plett 2004”来指代。
还由US2012/3168010A1和US2012/133369A1已知现有技术。
本发明旨在改进Plett 2004部分4.2中描述的方法,以便在降低实施该方法所需的计算功率的同时提高电池单元的充电状态的估计精度。因此,其主题之一是根据权利要求1所述的估计方法。
在所要求保护的方法中,用于计算预测值的观测模型是容量Cn,k3的线性函数。因此,对预测值的计算不带有与关于值Cn,k3的观测模型的线性化相关的误差。相反,Plett 2004的部分4.2.2中描述的观测模型是非线性的,并且因此必须关于容量Cn,k3进行线性化。这种线性化必然导致在对所修正的容量Cn,k3的估计中的误差。此外,观测模型的这种线性化操作消耗计算功率。
这种自动估计方法的实施例可以包括从属权利要求的特征中的一个或更多个。
所述自动估计方法的这些实施例另外具有以下优点:
-响应越过阈值而自动触发容量Cn,k3的估计的事实允许执行该估计的频率自动适应该电池单元的实际使用。这特别使得可以避免在不必要时对容量Cn,k3进行估计,并且因此限制执行该方法所需的计算功率而不降低对该容量的估计的精度。此外,只在充电状态阈值或电压阈值或充电流量阈值已被越过时才触发对容量Cn,k3的估计,使得可以在电池单元的每次充电和放电循环中仅进行一次估计,这足以精确地估计该容量。
-使用在对容量Cn,k3的两次相继估计之间形成的强度ik的N次测量值来估计该容量,允许提高该估计的精度。
-使用限定电池单元的几乎完整的放电循环的第一和第二预设阈值,允许在考虑该电池单元的这个几乎完整的放电循环期间所测量的所有值的同时对容量Cn,k3进行估计。这提高了容量Cn,k3的估计精度。
-在Plett 2004中,根据量dk3的测量值中的噪声的协方差来修正预测值Cn,k3,该预测值Cn,k3未被测量。在这些情况下,很难确定所述量dk3的测量噪声的协方差。相反,在所要求保护的方法中,使用的是强度ik的测量值中的噪声的协方差。强度ik的测量噪声的这种协方差易于通过在用于测量强度ik的电流计上进行的测试或由电流计制造商传送的数据来实验地确定。这因此使得实施所要求保护的方法的实施比在Plett 2004中描述的方法的实施容易得多。
本发明的另一主题是数据存储介质,该数据存储介质包括指令,所述指令用于在这些指令被电子处理器执行时执行以上自动估计方法。
本发明的另一主题是根据权利要求11所述的用于管理电池组的***。
最后,本发明的另一主题是包括所要求保护的电池组管理***的机动车辆。
细读仅作为非限制性示例给出的并参照附图做出的以下说明,将更好地理解本发明,在附图中:
-图1是配备有电气电池组的机动车辆的局部示意图,
-图2是图1的车辆的电池组的电池单元的电气模型的示意图,
-图3是用于估计图1的车辆的电池组的电池单元的充电状态的估计器的布置的示意图,
-图4至9示出图3的估计器所使用的不同状态模型和观测模型的等式,
-图10是用于借助于图3的估计器来估计电池单元的充电状态的方法的流程图,
-图11是图1的车辆的电池组的充电状态的确定方法的流程图,
-图12是用于对电池组的不同电池单元的充电状态估计的刷新时刻进行规划的方法的流程图,
-图13是示出借助于图12的方法规划的不同刷新时刻的时序图,
-图14是被用于估计图1的车辆的电池组的电池单元的充电状态的另一估计器的布置的示意图,
-图15和16分别示出图14的估计器所使用的状态模型和观测模型,
-图17是用于借助于图14的估计器来估计电池单元的充电状态的方法的流程图,
-图18是用于预测电池组的电池单元的内阻和电容的另一可行状态模型的图示。
在这些图中,相同的附图标记被用于指示相同的元件。在本说明书的其余部分中,本领域的技术人员所熟知的特征和功能未详细地描述。
在本说明书中,“计算功率”指电子计算器所执行的操作次数。因此,减少计算功率意味着减少为了获得相同结果或相同性质的结果所要执行的操作次数。
图1示出电牵引式机动车辆2,其更广泛地称为“电动车辆”。电动车辆是众所周知的,并且仅呈现了理解本说明书的其余内容所需的结构性元件。车辆2包括:
-电动马达4,其能够转动地驱动动力轮6以使得车辆2沿道路8行驶,以及
-电池组10,其给马达4提供电能。
电池组10包括用于电连接的两个端子12、14和在这些端子12和14之间电连接的多个电气电池单元。端子12和14连接到待供电的电负载。在此,这些端子因此特别地连接到电动马达4。
为了简化图1,仅示出四个电气电池单元18至21。这些电气电池单元通常被分组成多个级并且这些级串联地连接在端子12和14之间。在此,仅示出两个级。第一级包括电池单元18和19,而第二级包括电池单元20和21。每个级包括并联连接的多个分支。级的每个分支包括一个电气电池单元或串联的多个电气电池单元。在此,第一级包括两个分支,并且每个分支包括单个电气电池单元。在图1中示出的示例中,第二级与第一级在结构上是相同的。
在此,电池组10的所有电池单元除制造公差以外在结构上是相同的。因此,现在仅更详细地描述电池单元18。
电池单元18包括将其电连接到其他电池单元以及电池组10的端子12和14的两个电连接端子30、32。电池单元18还无自由度地机械固定到电池组10的其他电池单元,以形成通常所称的电池单元“包”。电池单元18能够在不被使用时储存电能。这种储存的电能进而被用于给马达4供电,这使得电池单元18放电。替代地,电池单元18也可以接收为其充电的电能。电池单元完全放电随后进行其完全重新充电构成所谓的充电/放电循环,或在下文中被简单地称作“电池单元循环”。
电池单元18是已知类型的电池单元,例如是LiPB(锂离子聚合物电池组)或其他电池单元。
电池单元18由初始标称容量Cn ini、初始内阻ROini、电流强度Imax、最大电压Umax、最小电压Umin和函数OCV(SOCk)表征。容量Cn ini是电池单元18的初始容量。电池单元的容量代表能够被储存在该电池单元中的最大电能量。该容量用Ah表达。随着电池单元18老化,即随着充电和放电循环的数量增加,电池单元的容量减小。电池单元18在时刻k的标称容量应记为Cn,k。
初始内阻ROini是电池单元18在其开始老化之前的内阻值。电池单元的内阻是电气电池单元的大部分电气模型中所具有的物理量。当电池单元老化时,内阻通常增大。电池单元18在时刻k的内阻记为ROk。
Imax是电池单元18在不会损坏的情况下所能够传送的最大电流强度。
Umax是可以持续存在于电池单元的端子30和32之间而不损坏该电池单元的最大电压。
电压Umin是电池单元18完全放电时端子30和32之间的最小电压。在下文中,我们将认为Imax、Umax、Umin是不随着时间过程而变化的恒定物理量。
OCV(SOCk)是根据电池单元18的充电状态SOCk来返回电池单元18空载电压的预定函数。空载电压是在电池单元18已经与任何电负载电绝缘两个小时之后在端子30和32之间能够测得的电压。
电池单元18在时刻k的充电状态被记为SOCk。充电状态代表电池单元18的填充率。其在储存在电池单元18中的电能量等于其容量Cn,k时等于100%。其在储存在电池单元18中的电能量为零,即不再能从电池单元18提取电能以给电负载供电时等于0%。
参数Cn ini、ROini、Imax、Umax、Umin和函数OCV(SOCk)是电池单元的已知参数。例如,这些参数由电池单元的制造商提供,或根据在该电池单元上实施的测量来实验性地确定。
电池组10还对于每个电池单元包括:
-测量该电池单元的端子之间的电压的电压计,以及
-测量该电池单元的充电或放电电流强度的电流计。
为了简化图1,仅示出电池单元18的电压计34和电流计36。
与电池单元18的所述各个参数不同,电池单元18的充电状态SOCk不能够被测量。因此需要对其进行估计。为此,车辆2包括用于管理电池组10的***40,或BMS(电池组管理***)。该***40的功能特别在于确定电池组10的充电状态以及该电池组的健康状态。为了确定该充电状态和该健康状态,***40能够估计电池组10的每个电池单元的充电状态和健康状态。电池单元的健康状态代表该电池单元的老化进程状态。在此,电池单元在时刻k的健康状态被记为SOHk。在下文中,该健康状态由比值Cn,k/Cn ini来测量。为了计算电池单元的健康状态,***40因此还能够估计该电池单元在当前时刻k的容量Cn,k。
为了执行这些各个估计,***40电连接到电池组10的每个电压计和每个电流计,以获取每个电池单元的端子之间的电压和电流强度的测量值。
在此,***40包括存储器42和能够执行保存在存储器42中的指令的可编程电子计算器44。为此,存储器42包含执行图10至12和/或图17的方法所必需的指令。该存储器42还包含执行这些方法所必需的不同参数的初始值。***40的结构因此与已知电池组管理***的结构相同或相似,并且不再更详细地说明。
图2示出电池单元18的电气模型50。该模型称为“一阶Thévenin模型”或“集中参数模型”。该模型包括从端子32开始直至端子30相继串联地连接的:
-空载电压OCV(SOCk)生成器52,
-并联RC电路54,以及
-内阻56,在下文中其在时刻k被称作“内阻ROk”。
电路54包括电容器CD,其与值为RD的电阻器并联连接。在下文中,我们将认为模型50的这两个参数CD和RD是已知的并且是随时间恒定的。电路54的端子处的在时刻k的电压记为VD,K。电池单元18的端子30和32之间的电压在时刻k的值记为yk,并且电池单元18的充电或放电电流强度在相同时刻记为ik。
图3示出在***40中实施的、用于估计电池单元18的充电状态和健康状态的估计器的布置的第一实施例。每个估计器都以由计算器执行的估计算法的形式来实施。由此,下文既会说“执行估计器”也会说“执行估计算法”。在该第一实施例中,***40包括用于基于电压的测量值yk和测量强度ik来估计充电状态SOCk和电压VD,K的估计器60。估计器60在此以卡尔曼滤波器的形式来实施。该估计器因此使用状态模型62(图4)和观测模型64(图5)。在这些图4和5中,通过使用此前限定的符号来表示这些模型的等式。符号ROk2和Cn,k3分别代表电池单元18分别在时刻k2和k3的内阻和容量。这些时刻k2和k3会在下文中限定。此外,在模型62中,xk指代时刻k的状态向量[SOCk,VD,K]T。在本说明书中,符号“T”指代数学转置运算。乘法运算用运算符“·”或“*”表示。
在下文中,我们会假设起始时间对应于时刻k的零值。在这些条件下,当前时刻k等于kTe,其中Te是用于电池组10的电流计和电压计的测量的采样周期。由此,Te是分隔通过***40获取电压和电流强度的任何两个相继时刻k和k-1的时间周期。周期Te通常是从0.1s到10s的常数。在此,周期Te等于1s±20%。例如,Te等于一秒。
在模型62中,wk是状态噪声向量。在此,噪声wk是居中的(centered)高斯白噪声。该噪声代表所用模型的不确定性。噪声wk在时刻k的协方差矩阵记为Qk。其由以下关系式限定:Qk=E(wk*wk T),其中E(…)是数学期望函数。模型62也以Xk+1=Fkxk+Bkik+wk的形式写出,其中:
-Fk是时刻k的状态转换矩阵,
-Bk是时刻k的控制向量。
模型62尤其允许根据之前的充电状态SOCk来预测时刻k+1的充电状态SOCk+1。
模型64允许根据充电状态SOCk、电压VD,K和测量强度ik来预测时刻k的电压值yk。在该模型中,vk是居中的测量高斯白噪声。在下文中,噪声vk在时刻k的协方差矩阵记为Rk。在此描述的具体情况中,该矩阵Rk是具有单列和单行的矩阵。其由关系式Rk=E(vk*vk T)来限定。该噪声vk独立于噪声wk和初始状态向量x0。
要指出的是,模型64是非线性的,这是因为函数OCV(SOCk)一般是非线性的。因此,估计器60实施卡尔曼滤波器的扩展版本,更常称为EKF(扩展卡尔曼滤波器)在该扩展版本中,通过使模型64在向量xk的邻域中线性化而以yk=Hkxk+ROk2.ik+vk的线性观测模型的形式结束。通常,模型64在向量xk的邻域中展开为泰勒级数。然后,忽略自第二阶起的导数的影响。在此,矩阵Hk因此等于函数OCV在充电状态SOCk的邻域中的一阶导数。模型64的这种线性化通常对于充电状态SOCk的每个新值实现。
为了能够估计充电状态SOCk+1,估计器60需要获悉容量Cn,k3和内阻ROk2。电池单元18的容量和内阻随着其老化而变化。为了考虑这种老化,电池单元18的容量和内阻分别在时刻k3和k2来估计。在此,估计器66根据测量值yk2、测量强度ik2和充电状态SOCk2来估计内阻ROk2。另一估计器68根据强度ik3和充电状态SOCk3来估计容量Cn,k3。
电池单元18的内阻和容量比其充电状态变化得更慢。由此,为了限制估计电池单元的充电状态所要求的计算功率而不因此损害该估计的精确性,估计器66和68比估计器60执行得较不频繁。在下文中,估计器66和68的执行时刻分别记为k2和k3,以与时刻k区分。在此,时刻k2的集合和时刻k3的集合是时刻k的集合的子集。因此,在两个相继时刻k2和k2-1之间以及两个相继时刻k3和k3-1之间经过了多个周期Te和多个时刻k。
这些估计器66和68各自也以卡尔曼滤波器的形式来实施。估计器66使用状态模型70(图6)和观测模型72(图7)。在这些模型中,噪声w2,k2和v2,k2是居中的高斯白噪声。噪声w2,k2和v2,k2的协方差在下文中分别记为Q2,k2和R2,k2。观测模型72允许我们对能够直接测量的物理量uk2的值进行预测。物理量uk2在此是最后N个测量值yk的和。其由以下关系式限定:
N是严格大于一的整数,其如下所述地来进行计数。在上述关系式中以及在模型72中,时刻k等于时刻k2。
模型72不仅考虑到在时刻k=k2测量的充电状态SOCk、电压VD,K和测度ik,而且还考虑到时刻k2与k2-1之间估计器60的前N个估计和前N个测量强度。对时刻k2与k2-1之间的中间测量和估计的考虑使得可以提高对内阻ROk2的估计的精确性。
估计器68使用状态模型74(图8)和观测模型76。在模型74和76中,噪声w3,k3和v3,k3是居中的高斯白噪声。噪声w3,k3和v3,k3的协方差在下文中分别记为Q3,k3和R3,k3。要注意的是,模型76是线性模型,使得可以对于估计器68使用简单的卡尔曼滤波器而不是扩展卡尔曼滤波器。
观测模型76允许我们估计可直接测量的物理量zk3。物理量zk3在此是最后N个测量强度ik的和。该物理量由以下关系式限定:
在上述关系式中并在模型76中,时刻k等于时刻k3。该物理量zk3不仅考虑到时刻k3之前的时刻k-1测量的强度ik-1,还考虑到时刻k3与k3-1之间测量的前N个强度。在此,N是严格大于一的整数,其按照以下将进一步说明的方式来计数。其不一定等于在模型72中引入的N。考虑时刻k3与k3-1之间的中间测量和估计使得可以提高容量Cn,k3的估计的精确性。
现在将借助于图10的方法并且在估计电池单元18的充电状态的具体情况中描述估计器60、66和68的运行。
该方法起始于调节执行估计器60、66和68所必需的各个协方差矩阵的阶段100。更准确地说,在操作102期间,借助于以下关系式来自动调节估计器60的协方差矩阵Qk和Rk:Qk=[N0G0,k(N0)]-1和Rk=I,其中:
-N0是严格大于1的预定整数,
-I是单位矩阵,以及
-G0,k(N0)由以下关系式来限定:
N0是一般在***40设计期间选择的,并且然后一次性设定。一般,N0小于100。例如,N0为5到15。在此,N0被选择为等于10。
如将在下文中看见的,使用前述关系式显著地简化对矩阵Q0和R0的调节以及对矩阵Qk的Rk调节。实际上,唯一要选择的参数是整数N0的值。
在操作104期间,协方差Q2,0和R2,0也被调节。例如,协方差Q2,0被选择为等于[(β*ROini)/(3*NC eol*NS)]2,其中:
-β是被选择为大于或等于0.3或0.5、并且优选地大于0.8、并且一般小于三的常数,
-NC eol是在电池单元18达到其寿命结束之前的、电池单元18的充电和放电循环的预测数量,以及
-NS是电池单元18的每个充电和放电循环估计内阻的次数。
以除以100的百分比表达的常数β代表初始内阻值ROini和其在寿命结束时的值之间的差。通常,β由使用者设定,或实验性地测得。NC eol是可以实验性地测得或根据电池单元18的制造商的数据获得的循环数。NS通过由计算器44实施的充电状态估计方法来设定。在该实施例中,如将在下文中看见的,每个循环仅估计内阻一次。因此,NS取值为等于1。
作为说明,协方差R2,0被选择为等于(2εmUmax/300)2,其中εm是以百分比表达的电压计34最大误差。
在下文中,协方差Q2,k2和R2,k2被视为恒定的且分别等于Q2,0和R2,0。
在操作106期间,协方差Q3,0和R3,0被调节。例如,协方差Q3,0被取值为等于[γ*Cn ini/(3*NC eol*NS)]2,其中,按除以100的百分比表达的γ代表容量Cn ini与电池单元18的寿命结束时的容量之间的差。γ是由使用者在0.05到0.8之间、优选为0.05到0.3之间选择的常数。在此,γ=0.2。
协方差R3,0例如被选择为等于[2*εim*Imax/300]2,其中εim是按百分比表达的电流计36的最大误差。
在下文中,协方差Q3,k3和R3,k3被视为常数,并且分别等于Q3,0和R3,0。
一旦已调节了协方差矩阵,就可以开始对电池单元18的充电状态进行估计。
在阶段110期间,在每个时刻k,电压计34和电流计36分别对值yk和强度ik进行测量,并且这些测量立即被***40获取并保存在存储器42中。阶段110在每个时刻k被重复。
并行地,估计器60执行估计电池单元18在时刻k的充电状态的阶段114。
为此,在步骤116期间,估计器60计算时刻k电池单元18的相应充电状态和电路54端子处的电压VD的预测值和预测值VD,k/k-1。在于此使用的符号中,指数k/k-1指示该预测值是通过仅考虑时刻0到k-1之间实现的测量而做出的由此谈及先验预测。指数k/k指示在时刻k的预测值考虑了时刻0到k之间实现的所有测量。由此说后验预测值。预测值和VD,k/k-1借助于模型62根据测量强度ik-1和容量Cn,k3来计算。要指出的是,在模型62中,无论k如何,状态转换矩阵Fk-1都是恒定的,并且因此不需要在每个时刻k进行重新估值。
在步骤117期间,估计器60还计算用于状态向量xk的估计误差的协方差矩阵的预测值Pk/k-1。通常,这借助于以下关系式实现:
Pk/k-1=Fk-1Pk-1/k-1Fk-1 T+Qk-1
这些不同的矩阵Fk-1、Pk-1/k-1和Qk-1此前已限定。
然后,在步骤118期间,估计器60通过使得模型64围绕预测值和VD,k/k-1线性化来构建矩阵Hk。
在步骤120期间,协方差矩阵Qk和Rk自动更新。在此,为此,步骤120与操作102相同,这次考虑到在步骤118期间构建的矩阵Hk。
在此之后,在步骤122期间,估计器60根据测量值yk与通过模型64预测的值k之间的差来修正预测值和VD,k/k-1。所述差被称为“更新量(innovation)”。该步骤122通常包括:
-计算预测值的操作124,然后
-修正预测值和VD,k/k-1以及矩阵Pk/k-1以获得经修正的预测值VD,k/k和Pk/k的操作126。
在操作124期间,预测值借助于模型64来计算,在该模型中,充电状态的值被取为等于并且电压VD,k的值被取为等于VD,k/k-1。测量值yk和其预测值之间的差在下文中被记为Ek。
存在用于基于更新量Ek来修正先验估计值和VD,k/k-1的许多方法。例如,在操作126期间,借助于卡尔曼增益Kk来修正这些估计。增益Kk由关系式Kk=Pk/k-1HT k(HkPk/k-1HT k+Rk)-1给出然后,借助以下关系式来修正先验预测值:xk/k=xk/k-1+KkEk。
借助于以下关系式来修正矩阵Pk/k-1:
Pk/k=Pk/k-1-KkHkPk/k-1。
步骤116至122在需要实现电池单元18充电状态的新估计时的每个时刻k被重复。在每次新的重复期间,通过在阶段114的先前重复期间对于电池单元18获得的值来初始化状态向量xk-1。
并行地,在步骤130期间,计算器44将强度的每个新测量值ik与预定电流阈值SHi比较。只要测量强度不越过该阈值SHi,就禁止估计器66的执行。相反地,一旦测量强度ik越过该阈值SHi,则立即执行估计器66。阈值SHi一般大于Imax/2,并且有利地大于0.8*Imax或0.9*Imax。
估计器66执行估计时刻k2的内阻ROk2的阶段140。在此,时刻k2等于强度ik越过阈值SHi的时刻k。
为此,在步骤142期间,估计器66根据模型70来计算内阻的先验预测值
然后,在步骤144期间,估计器66计算内阻估计误差的协方差矩阵的预测值P2,k2/k2-1。例如,该预测值借助于以下关系式来计算:P2,k2/k2-1=P2,k2-1/k2-1+Q2,0。要指出的是,在此,模型72是状态变量的线性函数。因此不需要将其在预测值的邻域中线性化以获得矩阵H2,k2。在此,该矩阵H2,k2等于-N。
在步骤148期间,估计器66根据测量的物理量uk2与该同一物理量的预测值之间的差来修正预测值在此,N是被严格选择为一、并且优选地大于10或30的预定常数。量uk2由估计器66随着值yk而测量和获取。
更准确地说,在操作150期间,计算器44获取测得量uk2,并计算预测值量uk2的获取是通过对测量值yk的最后N个测量值求和实现的。预测值是借助于模型72来计算的。在该模型72中,值ROk2取为等于之前计算的ROk2/k2-1。
然后,在操作152期间,估计器66根据更新量Ek2来修正预测值更新量Ek2等于测得量uk2与预测量之间的差。例如,在操作152期间,使用与在操作126期间实施的方法相同的方法。由此,在此不更详细地说明该操作152。然后,在估计器60随后的执行期间代替先前的估计ROk2-1/k2-1而使用新的估计ROk2/k2。
仅在测量强度ik高的时候才触发估计器66的执行允许提高估计内阻的精确性,同时减小实施该方法所需的计算功率。实际上,在强度ik更高时电流计的测量精确性更高。
同样与阶段110和114并行地,所述方法包括步骤160,在该步骤期间,在每个时刻k,估计值SOCk与预定的上阈值SHSOC进行比较。如果估计值SOCk下降到该阈值SHSOC之下,那么所述方法立即继续进行步骤162和164。否则,在下一时刻k重复步骤160。通常,阈值SHSOC为90%到100%。
在步骤162期间,计算器44通过将计数器初始化为零而开始,然后从该步骤开始起在每次测量新的强度ik时递增1。而且,在每个时刻k,测量强度ik和在相同时刻生成的估计值SOCk以与该时刻k关联的方式被保存在数据库中。
与步骤162并行地,在步骤164期间,计算器44将每个新的估计值SOCk与预定阈值SLSOC进行比较。阈值SLSOC例如为0%到10%。只要估计值SOCk保持高于该阈值SLSOC,那么就在下一时刻k重复步骤162。否则,一旦电池单元18的估计值SOCk下降到该阈值SLSOC以下,那么计算器44就立即触发估计器68的执行,并停止递增计数器。由此,只要未越过该阈值SLSOC,就禁止估计器68的执行。
在阶段166期间,估计器68估计时刻k3的容量Cn,k3。时刻k3因此等于触发估计器68的执行的时刻。
对于阶段140也相同,假定在每个时刻k未执行估计器68,则时刻k3-1不对应于时刻k-1。相反,时刻k3和k3-1分隔大于或等于NTe的时间间隔,其中N是在步骤162期间所计数的数量。
估计器68的卡尔曼滤波器的参数通过在阶段166的时刻k3-1的之前重复结束时获得的、这些参数的之前的值来初始化。
阶段166包括:
-在步骤170期间,借助于模型74来计算预测值Cn,k3/k3-1,
-在步骤172期间,计算容量的估计误差的协方差矩阵的预测值P3,k3/k3-1,以及
-在步骤174期间,修正预测值Cn,k3/k3-1和P3,k3/k3-1。
在步骤172和174期间,可观测性矩阵H3,k3等于[(SOCk-SOCk-N)]*3600/(NTe)。N在此是从估计的充电状态下降到阈值SHSOC之下的时刻到估计的充电状态已经下降到阈值SLSOC之下的时刻之间经过的时刻k的数量。值N等于在步骤162期间所计数的值。
步骤174包括获取测量的物理量zk3并计算量zk3的预测值的操作176。量zk3的获取在此包括计算时刻k-1与k-N之间测量的最后N个强度的和。预测值根据模型76来获得。
然后,在操作178期间,估计器68根据测得量zk3与预测量之间的差来修正预测值Cn,k3/k3-1,以获得容量Cn,k3/k3的后验估计值。该修正例如像在步骤126期间描述地那样实现。
然后,容量Cn,k3/k3被发送给估计器60,该估计器将其用于估计电池单元18在随后的时刻的充电状态。
仅在电池单元18大部分放电之后才触发估计器68的执行允许提高估计精确性,同时减小实施该方法所需的计算功率。
在阶段166结束时,在步骤180期间,计算器借助于以下公式来计算时刻k3的健康状态SOHk3:SOHk3=Cn,k3/Cn ini。
图11示出一种确定电池组10的充电状态的方法。在时刻k,根据该电池组的每个电池单元的充电状态来确定电池组10的充电状态。例如,这是以如下方式来实现的。在步骤190期间,计算器44确定该电池组的每个级的充电状态,这是通过将该级的每个电池单元的充电状态加和来实现的。
然后,在步骤192期间,电池组的充电状态被取值为等于在步骤190期间确定的级充电状态中的最小值。
如由图11的方法所示,对电池组在每个时刻k的充电状态的确定要求仅具有每个电池单元在时刻k的充电状态估计值。第一解决方案因此包括通过在每个时刻k执行阶段114来并行地对于每个电池单元执行图10的估计方法。然而,为了限制所要求的计算功率而不损害对于电池组确定的充电状态的精确性,还可以像关于图12的方法所述的那样规划电池单元充电状态估计的执行。
图12的方法是在这样的简化情况中描述的:在该简化情况中,仅使用三个优先级,其分别被称作高优先级、中等优先级和低优先级。而且,假设其优先级高的电池单元的充电状态需要在每个时刻k并且因此以频率fe来估计。其优先级中等的电池单元的充电状态仅需要须以三分之一的频率并且因此以频率fe/3来估计。最后,其优先级低的电池单元的充电状态需要以十分之一的频率并且因此以频率fe/10来估计。在该示例中,对于高和中等优先级,存在提前知悉的有限位置数。换句话说,被赋予高优先级的电池单元的数量被限制为提前预定的最大数量。对于被赋予中等优先级的电池单元的数量也是同样的。
为了规划需要刷新每个电池单元充电状态估计的时刻,计算器通过在步骤198期间给每个电池单元赋予优先级而开始。
步骤198起始于操作200,在该操作期间,***40获取每个电池单元端子之间的电压的测量值yk。
然后,在步骤202期间,如果测量值yk高于上阈值SHy或相反地低于下阈值SLy,则计算器44给该电池单元赋予高优先级,只要该级别还有足够的位置即可。阈值SHy大于或等于0.9*Umax,并且优选地大于0.95*Umax。至于阈值SLy,其大于或等于Umin,并小于1.1*Umin或1.05*Umin。重要之处在于频繁地刷新其电压接近Umax或相反地接近Umin的电池单元的充电状态的估计。实际上,在这样的情况中对电池单元充电状态的估计的误差可能会导致损害该电池单元的电气和机械特性。
然后,对于其他电池单元,在操作204期间,计算器44计算当前测量值yk与之前的值yk-X之间的电压差,其中X是大于或等于一的并且一般小于5或10的预定整数。在此,X=1。
在操作205期间,计算器44识别出相似电池单元。电池单元如果在相同时刻k具有相同的电压差和相同的测量值yk则被视为是“相似”。为此,在步骤205期间,计算器44将一个电池单元的电压差和测量值yk与其他电池单元在相同时刻的电压差和测量值yk进行比较,以在这些其他电池单元中识别出该电池单元的相似电池单元。该电池单元的标识符和该电池单元或被识别为其相似电池单元的标识符则合并为集合,该集合被保存在存储器42中。上述比较例如对于电池组10的、其标识符还未被包含到所保存的相似电池单元集合之一中的每个电池单元来实施。在下文中,仅给相似电池单元的每个集合中的电池单元之一赋予优先级。由此,操作206和随后的步骤208和210仅对于不具有任何相似电池单元的电池单元和对于相似电池单元的每个集合中的一个电池单元来实施。
在操作206期间,计算器将电池单元按在操作204期间计算的差的绝对值降序地分类。然后,计算器赋予该分类的第一电池单元与高优先级相关的剩余位置。然后,计算器将与中等优先级相关的剩余位置赋予该分类中随后的电池单元。最后,计算器给该分类最后的电池单元赋予低优先级。
一旦给每个电池单元都已经被赋予优先级,在步骤208期间,计算器44根据电池单元的优先级规划电池单元充电状态估计的刷新时刻。步骤208被实现以便遵循用于与每个优先级相关联的估计的刷新频率。为此,例如,计算器44首先预留出需要刷新高优先级的电池单元的估计的时刻。然后,考虑到已预留的刷新时刻,计算器预留出需要刷新中等优先级的电池单元的充电状态估计的时刻。最后,计算器以同样的方式对待被赋予低优先级的电池单元。
为了说明此,假设高优先级被赋予电池单元18、中等优先级被赋予电池单元19和20,并且低优先级被赋予电池单元21。而且,设定在周期期间,计算器最多执行两次图10的方法的阶段114。用这些假设所获得的结果在图13中示出。在该图中,时刻k至k+11在横坐标上被示出。在这些时刻k中的每一个之上,两个格子表示计算器44可以在每个时刻k执行两次图10方法的阶段114。这些格子的每一个中都指示出对其执行阶段114的电池单元的编号。当该格子中没有出现任何编号时,这意味着图10的方法没有执行,并且因此节省的计算功率可以用于其他目的(例如执行估计器66和68)。
最后,在步骤210期间,对于赋予了优先级的每个电池单元,计算器44在为该电池单元规划的时刻执行阶段114。在这些规划时刻之外,计算器禁止对该电池单元完整地执行阶段114。同样地,阶段114对于没有被赋予优先级的相似电池单元的执行也受到禁止。
并行地,在步骤212期间,对于没有被赋予优先级的每个相似电池单元,该电池单元的充电状态的估计值被取值为等于在步骤210期间对于该电池单元的相似电池单元所计算的最后的估计值。由此,仅对于相似电池单元中的一个执行阶段114。这使得可以减小确定电池组充电状态所要求的计算功率而不损害所述确定的精确性。
可选地,与步骤210并行地,在每个时刻k,计算器44还执行对在步骤210期间在该时刻k没有被处理的每个电池单元的充电状态进行预测的步骤214。步骤214包括对于在相同时刻没有执行完整估计阶段114的所有电池单元仅执行预测步骤116,而不执行修正步骤122。实际上,预测步骤116比步骤122消耗少得多的计算功率,并且因此可以例如在每个时刻k执行。由此,当步骤214被实施时,在每个时刻k都具有对电池组的每个电池单元的充电状态的新的估计值。
步骤198和208按照规则间隔重复,以更新赋予这些电池单元中每个的优先级并因此更新这些电池单元的充电状态估计的刷新频率。电池单元充电状态估计的该刷新时刻规划方法使得可以限制所需计算功率而不因此损害对于电池组确定的充电状态的精确性。实际上,图12的方法利用的是电压差小的电池单元是少量放电或充电并且因此其充电状态不快速地变化的电池单元。因此能够以较低的频率来估计这些电池单元的充电状态而不损害对于电池组而确定的充电状态的精确性。
在执行图10和11的方法的过程中,每次电池单元在给定时刻的充电状态SOCk需要被用于计算时,充电状态SOCk都被取值为等于对于该电池单元估计或预测的最后的充电状态。换句话说,认为充电状态在其被估计或预测的两个相继时刻之间保持恒定。
同样要指出的是,每次计算器44对于电池单元执行估计阶段114时,该计算器都从在该阶段对于同一电池单元的之前执行结束时获得的数值来检索该执行所必需的信息。这尤其是例如对于状态变量的情况。然而,要指出的是,之前的执行时刻不一定是时刻k-1,而根据赋予该电池单元的优先级可以是时刻k-3或k-10。
电池单元充电状态估计方法的许多其他实施例是可行的。例如,图14示出另一估计器布置。该另一布置与图3的布置相同,只是估计器66和68被单个估计器230代替。估计器230同时估计电池单元18的容量和内阻。估计器230比估计器60更不频繁地执行。在此,将执行估计器230的时刻记为k4,并且因此将在时刻k4估计的容量和内阻记为Cn,k4和ROk4。时刻k4的集合是时刻k的子集。
估计器230同时估计容量Cn,k4和内阻ROk4。该估计器230实施使用状态模型232(图15)和观测模型234(图16)的卡尔曼滤波器。
现在将参照图17的方法并且在电池单元18的具体情况中说明该估计器230的运行。图17的该方法与图10的方法相同,只是步骤130至174由步骤240、242、244和估计容量和内阻的阶段246来代替。
在步骤240期间,计算器44在每个时刻k将测量值yk与上阈值SHy2进行比较。通常,该阈值SHy2大于或等于0.8*Umax或0.9*Umax。仅在测量值yk下降到该阈值SHy2之下时才执行步骤242和244。
在步骤242期间,计算器44以初始化计数器为零开始,然后在每个新的时刻k使该计时器递增1。而且,在这些时刻k中的每一个,测量强度ik、值yk、充电状态SOCk和估计的电压VD,K以与该时刻k关联的方式被保存在数据库中。
与步骤242并行地,在步骤244期间,计算器44在每个时刻k将新的测量值yk与低电压阈值SLy2进行比较。该阈值SLy2小于或等于1.2*Umin或1.1*Umin,并且大于或等于Umin。
一旦测量值yk下降到阈值SLy2之下,计数器在步骤242期间的递增就停止,并且估计器230的执行就被触发。另一方面,只要测量值yk保持大于该阈值SLy2,就禁止估计器230的执行。
估计器230执行阶段246。如前所述的,要指出的是,时刻k4和k4-1分隔大于或等于NTe的时间间隔,其中N是在步骤242期间增加的计数器值。估计器230的运行能够从此前对于估计器66和68所述的运行中推导出。因此不再在此更详细说明。
其他电气模型和因此其他状态模型可以被用于估计电池单元18的充电状态。例如,在一个简化变型中,电路54被省去。另一方面,更为复杂的电气模型可包含彼此串联地电连接的多个并联RC电路。电池单元18的状态模型则因此必须被修改以对应于电池单元的该新电气模型。然而,所有在上文中所述的内容都轻易地适用这样的经修改的状态模型。对于经修改的状态模型的示例,读者可参考专利申请WO2006057468。
模型50的参数RD和CD也可以被估计而不是被视为是预定的恒定参数。为此,这两个参数RD和CD例如被引入到状态向量xk中,该状态向量由此变成[SOCk,VD,k,RD,k和CD,k]T。例如,状态模型被修改为包含以下两个等式:RD,k+1=RD,k和CD,k+1=CD,k。
状态向量xk也可以由电池单元的温度来补充,以在估计该电池单元的充电状态的同时估计该温度。
电池单元也可以配备有比如温度传感器的附加传感器。在该情况中,对观测模型进行修改以考虑这些附加的测量物理量。读者就更改的观测模型示例可以参考专利申请WO2006057468。
用于模型化电气电池单元的其他可行电气模型还在Plett 2004第2部分的第3.3章中说明。
协方差矩阵Rk和Qk的连续自动调节可以以不同的方式来实现。例如,可应用在以下文献中描述的所谓的“协方差匹配”的方法:Mehra,R.K:“On the identification ofvariances and adaptative Kalman Filtering”,Automatic Control,IEEE Transactionon,第15卷No.2,第175-184页,1970年4月。例如像在步骤102期间描述地那样,该方法在对矩阵R0和Q0进行初始设定之后应用。
在另一变型中,矩阵Q0、R0、Qk和Rk不是像参照操作102和120所述的那样进行调节。例如,通过实施常规方法来调节这些矩阵。在简化的情况中,这些矩阵是恒定的。例如,矩阵R0则利用由传感器制造商提供的数据或基于在这些传感器上实施的测试来调节,而矩阵Q0则通过连续的测试来调节。
修正预测值的步骤122或178可以以不同的方式来实施。例如,在优选方法中,对充电状态的和电压VD,k的预测值的修正通过最小化包括以下两个项的二次成本函数J来实现:
-与测量值的预测误差相关的一项,和
-与状态向量的估计误差相关的另一项。
该方法在以下书籍的第10.5.2章中详细地说明:Y.Bar-Shalom等:“EstimationWith Applications to Tracking and Navigation”,Theory Algorithms and Software,Wiley Inter-science,2001年。
在另一变型中,估计器60不以卡尔曼滤波器的形式来实施。例如,充电状态通过以RII(无限脉冲响应)滤波器的形式模拟其在时间上的变化来估计,该滤波器的系数由RLS(递推最小二乘估计)方法来估计。
其他状态模型可被用于估计电池单元的内阻和容量。例如,模型232可以由在图18中示出的模型250来代替。在模型250中,α、β和γ是其值根据电池单元的制造商的数据或实验性测得的数据来获得的常数。通常:
-α等于1±30%或10%,
-β也等于1±30%或10%,以及
-γ通常为0.1到0.5。例如,γ等于0.2±30%或10%。
在模型250中,Nc k等于在时刻k之前实现的电池单元的充电/放电循环数。该循环数例如通过计数电池单元的充电状态下降到上阈值SHSOC之下、然后下降到下阈值SLSOC之下的次数来测量。wa d,k是居中的高斯白噪声。γ是按除以100的百分比表达的、电池单元初始容量Cn ini与其寿命结束时的容量之间的差。该模型考虑到以下事实:
-内阻随着电池单元老化而增大,以及
-电池单元的容量随着电池单元老化而减小。
类似地,状态模型70可以由以下状态模型来代替:ROk2+1=(α+βNC k2/NC EOL)ROk2+w2,k2,其中该模型的各个符号之前已经说明。
状态模型74可以由以下状态模型来代替:Cn,k3+1=(1-γNC k3/NC EOL)Cn,k3+v3,k3,其中该模型的各个符号之前已经说明。
根据估计器68所使用的观测模型,量zk3可以以不同的方式来计算。例如,量zk3等于时刻k与k-N+1之间测得的最后N个强度的和。在该情况中,当N等于1时,zk3=ik3。
以上对于协方差矩阵Qk和Rk所述的内容也可以适用于对估计器68和230的协方差矩阵进行初始化。
作为变型,估计器68不以卡尔曼滤波器的形式来实施。例如,容量通过以RII(无限脉冲响应)滤波器的形式模拟其在时间上的变化来估计,该滤波器的系数由RLS(递推最小二乘估计)方法来估计。
图10和17的方法可以通过将N取值为等于预定常数而受到简化。在该情况中,N不被计数,并且步骤160、162、240和242可以被省去。例如,N被选择为等于一,或相反地被选择为严格大于1或5或10。
在另一变型中,在时刻k3与k3-1之间的每个时刻k,仅执行计算预测值Cn,k的步骤170,而不执行修正该预测值的步骤174。由此,获得电池单元在这些时刻k中每一个处的容量的新预测值,同时限制所需的计算功率。类似地,在时刻k4与k4-1之间的每个时刻k,仅执行计算容量的和内阻的预测值的步骤,而不执行修正这些预测值的步骤。由此,在这些变型中,在每个时刻k预测电池单元的容量,但该预测值仅在时刻k3或k4才被修正。该容量的估计算法因此在时刻k3与k3-1之间或k4与k4-1之间仅部分地被执行,并且仅在时刻k3或k4才被完整地执行。
在时刻k3与k3-1之间或时刻k4与k4-1之间的每个时刻k,容量可以通过执行第一算法来估计,然后在时刻k3或k4,容量通过执行不同于第一算法并且需要更大的计算功率的第二算法来估计。第一和第二算法不一定如前所述地分别对应于卡尔曼滤波器的单个步骤170和阶段166或246。这些算法也可以是两个完全不同的估计算法。
估计电池单元的容量的步骤166或246可以如参照图10所述地响应越过充电状态阈值而被触发,或如参照图17所述地响应越过电压阈值而被触发。这些步骤166和246也可以响应越过电流输出量阈值而被触发。为此,自电池单元的电压或充电状态下降到预定上阈值之下时的时刻开始,在每个时刻k,计算器44借助于以下关系式来计算电流输出量QCk:QCk=QCk-1+ikTe。一旦QCk穿过上阈值SHQ,则执行阶段166或246。另一方面,只要量QCk保持大于阈值SHQ,就禁止阶段166或246的执行。作为变型,量QCk也可以在包含最后N个时刻k的滑动窗口上计算,其中N为预定常数。
在另一实施例中,省略响应越过阈值而触发对容量和/或内阻的估计。例如,这些估计以规则的间隔触发。如果可用的计算功率足以在每个时刻k估计该容量和该内阻,那么该规则间隔等于Te。
图12的方法的众多其他实施例是可行的。例如,可以省略操作205。在该情况中,不识别相似电池单元,并且还省略步骤212。
操作202可用以不同的方式来实施。例如,仅使用上阈值和下阈值中的一个。也可以省略操作202。
优先级的数量可以是任意的,并且可以大于至少两个或三个。用于给电池单元赋予优先级的其他方法是可行的。例如,电池单元的优先级可以借助于公式来计算,所述公式使该电池单元的优先级关联到其电压差和其电压。在该后一情况中,省略比较操作。
所述用于根据电池单元的优先级使刷新时刻与电池单元关联的方法仅是一个示例。用于根据任务优先级对这些任务进行排序的任何其他已知方法可以适用于在此所述的对用于估计电池单元充电状态的刷新时刻进行排序的情况。
可以省略参照图12所述的用于估计每个电池单元充电状态的刷新时刻的规划。例如,当在每个时刻k估计每个电池单元的充电状态所需的计算功率可用时情况将会如此。
作为变型,计算器44包括多个可编程子计算器,其每个都能够并行地对于相应的电池单元执行图10或17的估计方法。
电池单元的健康状态也可以借助于以下关系式来计算:SOHK=ROK/ROini。
电池组10可以由任何类型的电池组(例如:铅电池组、超大容量电池组或燃料电池)来代替。在该情况中,估计器60的状态模型和/或观测模型可以可选地适用于考虑电池组的技术。
以上描述的内容也适用于混合动力车辆,即其动力轮的驱动由电动马达和热内燃机同时或交替提供的车辆。车辆2也可以是卡车、摩托车或三轮车,并且一般而言为能够通过借助于由电池组供电的电动马达来驱动动力轮而移动的任何车辆。例如,其可以是起重机。
电池组10可以借助于可以将其连接到输电干线的电源插座来进行再充电。电池组10也可以由热内燃机来再充电。
Claims (12)
1.一种电池组的电池单元的充电状态的自动估计方法,所述方法包括:
a)在时刻k获取(110)所述电池单元的端子之间的电压的测量值yk和对所述电池单元进行充电或放电的电流的测量强度ik;以及
b)估计(166;246)所述电池单元在时刻k3的容量Cn,k3,所述容量Cn,k3以Ah表示,其代表所述时刻k3能够存储在所述电池单元中的最大电能量,所述时刻k3是最接近所述时刻k的已经估计所述容量Cn,k3的时刻,并且所述时刻k3不如所述时刻k频繁,所述估计包括:
使用将所述容量Cn,k3与同一电池单元在前一时刻k3-1的容量Cn,k3-1相关联的状态模型来计算(170)所述容量Cn,k3的预测值,然后
根据对所述电池单元进行充电或放电的电流在最接近所述时刻k3的时刻的测量强度来修正(174)所述容量Cn,k3的预测值;以及
c)根据所述测量值yk、所述测量强度ik和所述电池单元的容量Cn,k3来估计(114)所述电池单元在所述时刻k的充电状态SOCk;
其特征在于,对所述容量Cn,k3的预测值的修正包括:
-在所述时刻k3获取(176)由以下关系式限定的可测量物理量zk3:
<mrow>
<msub>
<mi>z</mi>
<mrow>
<mi>k</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>k</mi>
<mo>-</mo>
<mi>N</mi>
</mrow>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>i</mi>
<mi>m</mi>
</msub>
</mrow>
其中k是最接近所述时刻k3的时刻,并且N是大于或等于1的整数,当N等于一时,zk3等于ik-1;
-使用以下观测模型来计算(176)所述可测量物理量zk3的预测值
<mrow>
<msub>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>SOC</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<mi>SOC</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>NT</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mn>.3600.</mn>
<msub>
<mi>C</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>k</mi>
<mn>3</mn>
</mrow>
</msub>
</mrow>
其中:
-k是最接近所述时刻k3的时刻,
-SOCk-N是所述电池单元在先前时刻k-N的充电状态,以及
-Cn,k3是所述观测模型中所述电池单元在所述时刻k3的容量的预测值,即,就在进行修正之前计算的预测值;以及
-根据所获取的物理量zk3和计算自所述观测模型的预测值之间的差来修正(178)所述容量Cn,k3的预测值。
2.根据权利要求1所述的方法,其中:
-如果参数的值没有越过第一预设阈值,则禁止(164;244)对所述容量Cn,k3的预测值进行修正,从所述测量值yk、所估计的充电状态SOCk以及在时刻k和先前时刻之间由所述电池单元所传递的电流量Qk中选择所述参数,并且
-响应以下事实而触发(164;244)对所述容量Cn,k3的预测值的修正:当所述参数是所述测量值yk或所述充电状态SOCk的估计值时,所述参数已经下降到低于所述第一预设阈值,或者当所述参数是通过的电流量Qk时,所述参数已经升高到高于所述第一预设阈值。
3.根据权利要求2所述的方法,其中,所述方法包括对在所述测量值yk或所估计的充电状态SOCk下降到低于第二预设阈值的时刻t1与所述测量值yk或所估计的充电状态SOCk下降到低于所述第一预设阈值的随后时刻t2之间的步骤a)的迭代次数进行计数(162;242),然后,将所述容量Cn,k3的预测值的修正(178)中的整数N选择为等于在所述时刻t1和t2之间所计数的迭代次数并且严格大于一。
4.根据权利要求3所述的方法,其中,所述第一预设阈值对应于所述电池单元的低于20%的充电状态,并且所述第二预设阈值对应于所述电池单元的高于90%的充电状态。
5.根据上述权利要求中任一项所述的方法,其中,所述方法包括限定(106)所述强度ik的测量中的噪声的协方差,然后根据所述强度ik的测量中的噪声的所述协方差来修正(178)所述容量Cn,k3的预测值。
6.根据权利要求5所述的方法,其中,根据所获取的物理量zk3和所计算的预测值之间的差,使用以下关系式来修正所述容量Cn,k3的预测值:Cn,k3/k3=Cn,k3/k3-1+Kk3Ek3,其中:
-Cn,k3/k3是所述容量Cn,k3的修正的预测值,
-Cn,k3/k3-1是所述容量Cn,k3的修正之前的预测值,
-Ek3是所获取的物理量zk3和所计算的预测值之间的差,并且
-Kk3是根据所述强度ik的测量中的噪声的协方差针对所述时刻k3计算的卡尔曼增益。
7.根据上述权利要求中任一项所述的方法,其中,在对所述容量Cn,k3的预测值的计算中所使用的状态模型如下:Cn,k3=Cn,k3-1,并且所述方法包括:
-根据由所使用的状态模型引入的建模误差的噪声的协方差Q3,0来计算(172)在所述容量Cn,k3的估计值中的误差的协方差P3,k3的预测值;以及
-使用以下关系式来自动限定(106)所述协方差Q3,0:Qd,k3=[(γCn ini)/(3Nc eolNs)]2,其中:
-γ是预设常数,其值严格小于一并且大于零,
-Cn ini是所述电池单元的最大标称容量,
-Nc eol是在所述电池单元到达其寿命终点之前所述电池单元的放电循环的次数,所述次数是预设的,以及
-Ns是所述电池单元的一次放电循环中步骤c)的迭代次数。
8.根据上述权利要求中任一项所述的方法,其中,在所述时刻k3和k3-1之间执行至少N次步骤a)和b),其中N是大于三的整数。
9.根据上述权利要求中任一项所述的方法,其中,所述方法还包括使用以下关系式来计算(180)所述电池单元的健康状态SOHk3:SOHk3=Cn,k3/Cn ini,其中:
-Cn,k3是在最近时刻k3估计的所述电池单元的容量,以及
-Cn ini是所述电池单元的最大标称容量。
10.一种数据存储介质(42),其特征在于,其包括指令,所述指令用于在这些指令被电子处理器执行时执行如上述权利要求中任一项所述的估计方法。
11.一种用于管理配备有至少一个电池单元的电池组的***,所述***包括电子处理器(44),其被编程为:
a)在时刻k获取所述电池单元端子的端子之间的电压的测量值yk和对所述电池单元进行充电或放电的电流的测量强度ik;以及
b)估计所述电池单元在时刻k3的容量Cn,k3,所述容量Cn,k3以Ah表示,其代表所述时刻k3能够存储在所述电池单元中的最大电能量,所述时刻k3是最接近所述时刻k的已经估计所述容量Cn,k3的时刻,并且所述时刻k3不如所述时刻k频繁,所述估计包括:
使用将所述容量Cn,k3与同一电池单元在先前时刻k3-1的容量Cn,k3-1相关联的状态模型来计算所述容量Cn,k3的预测值,然后
根据对所述电池单元进行充电或放电的电流在最接近所述时刻k3的时刻的测量强度来修正所述容量Cn,k3的预测值;以及
c)根据所述测量值yk、所述测量强度ik和所述电池单元的容量Cn,k3来估计所述电池单元在所述时刻k的充电状态SOCk;
其特征在于,所述电子处理器(44)还被编程为通过执行以下步骤来修正所述容量Cn,k3的预测值:
-在所述时刻k3获取由以下关系式限定的可测量物理量zk3:
<mrow>
<msub>
<mi>z</mi>
<mrow>
<mi>k</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>k</mi>
<mo>-</mo>
<mi>N</mi>
</mrow>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>i</mi>
<mi>m</mi>
</msub>
</mrow>
其中k是最接近所述时刻k3的时刻,并且N是大于或等于1的整数,当N等于一时,zk3等于ik-1;
-使用以下观测模型来计算所述可测量物理量zk3的预测值
<mrow>
<msub>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>SOC</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<mi>SOC</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>NT</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mn>.3600.</mn>
<msub>
<mi>C</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>k</mi>
<mn>3</mn>
</mrow>
</msub>
</mrow>
其中:
-k是最接近所述时刻k3的时刻,
-SOCk-N是所述电池单元在先前时刻k-N的充电状态,以及
-Cn,k3是所述观测模型中所述电池单元在所述时刻k3的容量的预测值,即,就在进行修正之前计算的预测值;以及
-根据所获取的物理量zk3和计算自所述观测模型的预测值之间的差来修正所述容量Cn,k3的预测值。
12.一种机动车辆,包括:
-至少一个驱动轮(6);
-电动马达(4),其能够驱动所述驱动轮转动以使得所述机动车辆移动;
-电池组(10),其包括能够储存电能并交替地释放电能以给所述电动马达供电的至少一个电池单元(18-21),所述电池单元包括两个端子(30、32),所述电池单元通过所述两个端子电连接到所述电动马达;
-电压计(34),其电连接在所述电池单元的端子之间以测量这些端子之间的电压;
-电流计(36),其与电气电池单元串联连接以测量对所述电池单元进行充电或放电的电流的强度;以及
-用于管理所述电池组的***(40),所述***连接到所述电压计和所述电流计,并且包括能够根据所述电压计和电流计的测量值来估计所述电池组的电池单元的充电状态的可编程电子处理器(44),
其特征在于,用于管理所述电池组的所述***(40)是根据权利要求11所述的***。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1461616 | 2014-11-28 | ||
FR1461616A FR3029297B1 (fr) | 2014-11-28 | 2014-11-28 | Procede automatique d'estimation de l'etat de charge d'une cellule d'une batterie |
PCT/FR2015/053240 WO2016083754A1 (fr) | 2014-11-28 | 2015-11-26 | Procede automatique d'estimation de l'etat de charge d'une cellule d'une batterie |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107209227A true CN107209227A (zh) | 2017-09-26 |
CN107209227B CN107209227B (zh) | 2020-03-17 |
Family
ID=52392120
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580069203.9A Active CN107209227B (zh) | 2014-11-28 | 2015-11-26 | 电池组的电池单元的充电状态的自动估计方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US10379171B2 (zh) |
EP (1) | EP3224635B1 (zh) |
JP (1) | JP6831325B2 (zh) |
KR (1) | KR102534688B1 (zh) |
CN (1) | CN107209227B (zh) |
FR (1) | FR3029297B1 (zh) |
WO (1) | WO2016083754A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114325395A (zh) * | 2020-09-30 | 2022-04-12 | 北京昇科能源科技有限责任公司 | 一种电池状态的确定方法及装置 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3029315B1 (fr) * | 2014-11-28 | 2016-12-09 | Renault Sa | Procede automatique d'estimation de la capacite d'une cellule d'une batterie |
FR3029299B1 (fr) * | 2014-11-28 | 2016-12-09 | Renault Sa | Procede automatique de determination de l'etat de charge d'une batterie |
CN106970332B (zh) * | 2017-05-04 | 2019-09-03 | 成都雅骏新能源汽车科技股份有限公司 | 基于回馈自适应的充电soc修正方法 |
CN109444747B (zh) * | 2018-11-09 | 2021-04-20 | 北京经纬恒润科技股份有限公司 | 一种对象的信息检测方法及装置 |
US11069926B1 (en) * | 2019-02-14 | 2021-07-20 | Vcritonc Alpha, Inc. | Controlling ongoing battery system usage via parametric linear approximation |
US20220229114A1 (en) * | 2019-06-12 | 2022-07-21 | Volvo Truck Corporation | A method for estimating a battery state |
WO2020249203A1 (en) * | 2019-06-12 | 2020-12-17 | Volvo Truck Corporation | A method for estimating a battery state |
CN110632520A (zh) * | 2019-10-17 | 2019-12-31 | 西安工业大学 | 一种动力电池soc的估算装置及其估算方法 |
CN113884884B (zh) * | 2021-10-21 | 2022-07-26 | 山东大学 | 一种基于相关性的动力电池组故障诊断方法及*** |
FR3144868A1 (fr) | 2023-01-05 | 2024-07-12 | Entroview | Procédé automatique d'estimation de l'état de charge d'une cellule d'une batterie |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101116003A (zh) * | 2005-01-27 | 2008-01-30 | 松下电动车辆能源股份有限公司 | 二次电池充电/放电电量的估计方法和装置,二次电池极化电压的估计方法和装置,以及二次电池剩余容量的估计方法和装置 |
CN101598769A (zh) * | 2009-06-29 | 2009-12-09 | 杭州电子科技大学 | 一种基于采样点卡尔曼滤波的电池剩余电量估计方法 |
US20120133369A1 (en) * | 2010-11-30 | 2012-05-31 | GM Global Technology Operations LLC | Algorithm for determining the capacity of a battery while in service |
US20120316810A1 (en) * | 2011-06-08 | 2012-12-13 | GM Global Technology Operations LLC | Battery limit calibration based on battery life and performance optimization |
CN103852727A (zh) * | 2014-02-14 | 2014-06-11 | 清华大学深圳研究生院 | 动力电池荷电状态在线估算方法和装置 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0843504A (ja) * | 1994-07-27 | 1996-02-16 | Nippondenso Co Ltd | 電気自動車用電池残存容量計 |
EP1512189A2 (en) * | 2002-03-29 | 2005-03-09 | Estco Battery Management, Inc. | Fuel cell operating control system |
US8103485B2 (en) * | 2004-11-11 | 2012-01-24 | Lg Chem, Ltd. | State and parameter estimation for an electrochemical cell |
US7315789B2 (en) * | 2004-11-23 | 2008-01-01 | Lg Chem, Ltd. | Method and system for battery parameter estimation |
BRPI0419140B8 (pt) | 2004-11-29 | 2023-01-10 | Lg Chemical Ltd | Método e aparato para estimar os valores descritivos da presente condição operacional de uma bateria |
KR100804698B1 (ko) * | 2006-06-26 | 2008-02-18 | 삼성에스디아이 주식회사 | 배터리 soc 추정 방법 및 이를 이용하는 배터리 관리시스템 및 구동 방법 |
KR100962856B1 (ko) * | 2008-04-03 | 2010-06-09 | 현대자동차주식회사 | 배터리의 잔존용량 추정 방법 |
KR101671998B1 (ko) * | 2010-04-02 | 2016-11-04 | 에스케이이노베이션 주식회사 | 배터리의 용량 열화 상태 측정 장치 및 방법 |
JP5535968B2 (ja) * | 2011-03-08 | 2014-07-02 | 三菱重工業株式会社 | 充電率推定装置、充電率推定方法、及びプログラム |
US20120316800A1 (en) | 2011-06-07 | 2012-12-13 | Gregory David Shteinhauz | System for predicting vehicle vibration or acoustic response |
JP5595361B2 (ja) | 2011-09-27 | 2014-09-24 | プライムアースEvエナジー株式会社 | 二次電池の充電状態推定装置 |
US20140236511A1 (en) * | 2011-09-30 | 2014-08-21 | Kpit Cummins Infosystems Ltd. | System and method for determining state of charge of a battery |
FR3029299B1 (fr) * | 2014-11-28 | 2016-12-09 | Renault Sa | Procede automatique de determination de l'etat de charge d'une batterie |
FR3029315B1 (fr) * | 2014-11-28 | 2016-12-09 | Renault Sa | Procede automatique d'estimation de la capacite d'une cellule d'une batterie |
FR3029296B1 (fr) * | 2014-11-28 | 2016-12-30 | Renault Sa | Procede automatique d'estimation de l'etat de charge d'une cellule d'une batterie |
FR3029298B1 (fr) * | 2014-11-28 | 2016-12-30 | Renault Sa | Procede automatique d'estimation de l'etat de charge d'une cellule d'une batterie |
-
2014
- 2014-11-28 FR FR1461616A patent/FR3029297B1/fr not_active Expired - Fee Related
-
2015
- 2015-11-26 CN CN201580069203.9A patent/CN107209227B/zh active Active
- 2015-11-26 US US15/531,182 patent/US10379171B2/en active Active
- 2015-11-26 KR KR1020177017126A patent/KR102534688B1/ko active IP Right Grant
- 2015-11-26 EP EP15808756.9A patent/EP3224635B1/fr active Active
- 2015-11-26 JP JP2017528422A patent/JP6831325B2/ja active Active
- 2015-11-26 WO PCT/FR2015/053240 patent/WO2016083754A1/fr active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101116003A (zh) * | 2005-01-27 | 2008-01-30 | 松下电动车辆能源股份有限公司 | 二次电池充电/放电电量的估计方法和装置,二次电池极化电压的估计方法和装置,以及二次电池剩余容量的估计方法和装置 |
CN101598769A (zh) * | 2009-06-29 | 2009-12-09 | 杭州电子科技大学 | 一种基于采样点卡尔曼滤波的电池剩余电量估计方法 |
US20120133369A1 (en) * | 2010-11-30 | 2012-05-31 | GM Global Technology Operations LLC | Algorithm for determining the capacity of a battery while in service |
US20120316810A1 (en) * | 2011-06-08 | 2012-12-13 | GM Global Technology Operations LLC | Battery limit calibration based on battery life and performance optimization |
CN103852727A (zh) * | 2014-02-14 | 2014-06-11 | 清华大学深圳研究生院 | 动力电池荷电状态在线估算方法和装置 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114325395A (zh) * | 2020-09-30 | 2022-04-12 | 北京昇科能源科技有限责任公司 | 一种电池状态的确定方法及装置 |
CN114325395B (zh) * | 2020-09-30 | 2024-01-26 | 北京昇科能源科技有限责任公司 | 一种电池状态的确定方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
KR20170092589A (ko) | 2017-08-11 |
FR3029297B1 (fr) | 2016-12-30 |
EP3224635B1 (fr) | 2019-05-01 |
KR102534688B1 (ko) | 2023-05-19 |
CN107209227B (zh) | 2020-03-17 |
WO2016083754A1 (fr) | 2016-06-02 |
EP3224635A1 (fr) | 2017-10-04 |
FR3029297A1 (fr) | 2016-06-03 |
US20180031639A1 (en) | 2018-02-01 |
JP2017538929A (ja) | 2017-12-28 |
US10379171B2 (en) | 2019-08-13 |
JP6831325B2 (ja) | 2021-02-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107209227A (zh) | 电池组的电池单元的充电状态的自动估计方法 | |
CN107110915A (zh) | 电池组的电池单元的充电状态的自动估计方法 | |
CN107110916A (zh) | 电池组的电池单元的充电状态的自动估计方法 | |
CN107003360A (zh) | 电池组的充电状态的自动确定方法 | |
US10838010B2 (en) | Automatic method for estimating the capacitance of a cell of a battery | |
CN107690585B (zh) | 用于确定锂硫电池组的健康状况和充电状态的方法和装置 | |
CN103270668B (zh) | 为电池快速充电的设备和方法 | |
JP2023518778A (ja) | 電池充電状態を決定する方法及び装置、電池管理システム | |
CN104267261B (zh) | 基于分数阶联合卡尔曼滤波的二次电池简化阻抗谱模型参数在线估计方法 | |
JP2017538931A5 (zh) | ||
CN105319515A (zh) | 锂离子电池荷电状态和健康状态联合估算方法 | |
JP6711981B2 (ja) | バッテリのパラメータ推定装置 | |
CN112666475A (zh) | 一种电动车的电池组的荷电状态估算方法 | |
CN113466726B (zh) | 确定电池等效电路模型的参数的方法、装置、存储介质及电子设备 | |
Maithani et al. | Advanced State Estimation Methods for Lithium-Ion Battery Cell-A Comprehensive Review | |
CN105487013A (zh) | 用于确定能量存储设备的运行工况的方法和*** | |
CN113945854A (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 |