具体实施方式

参考结合附图详细描述的实施例，本公开的优点和特征以及实现它们的方法将变得明显。然而，应当理解的是，本公开不限于下面呈现的实施例，而是可以以各种不同的形式实施，并且包括本公开的精神和范围中包括的所有变型、等同物和替代物。下面给出的实施例被提供以完成发明的公开，并且被提供以充分地告知本公开所属领域的普通技术人员发明的范围。在本公开的描述中，相关领域的某些详细解释在被认为其会不必要地模糊发明的本质时被省略。

本说明书中使用的术语仅用于描述特定实施例，并且不旨在限制本公开。以单数使用的表述涵盖复数的表述，除非其在上下文中具有明显不同的含义。在本说明书中，将理解的是，诸如“包括”或“具有”等的术语旨在指示说明书中公开的特征、数量、步骤、动作、组件、部分或其组合的存在，并且不旨在排除可以存在或可以添加一个或更多个其他特征、数量、步骤、动作、组件、部分或其组合的可能性。虽然“第一”、“第二”等这样的术语可以用于描述各种组件，但是这样的组件不必限于上述术语。上述术语仅用于将一个组件与另一组件区分开。

在下文中，下面将参照附图更详细地描述本公开的实施例。相同或对应的那些组件被赋予相同的附图标记，而与图号无关，并且省略了冗余的解释。

图1是根据实施例的用于执行用于检测异常电池单体的方法的电池系统的示意性配置图。

参照图1，电池系统100可以包括电池110、电压测量单元120、电流测量单元130、微处理器140和存储单元150。

作为用于存储电力的单元的电池110可以包括在第一端子101与第二端子102之间彼此电连接的多个电池单体111。多个电池单体111可以串联连接、并联连接或者以串联和并联的组合连接。每个电池单体111可以理想地放电和充电相同大小的电流，并且可以具有相同大小的容量。然而，实际上，电池单体111的内部状态可以不同，因此，它们可以具有不同的内部电阻和电动势。此外，每个电池单体111可以具有不同的G参数值和H参数值。

每个电池单体111可以包括可充电的二次电池。例如，电池单体111可以包括镍-镉电池、铅酸电池、镍金属氢化物(NiMH)电池、锂离子电池或锂聚合物电池。可以根据电池110所需的电容、输出电压和输出电流来确定电池单体111的数量。

图1中示出了一个电池110。然而，多个电池110可以并联和/或串联连接，并且可以通过第一端子101和第二端子102连接到负载和/或充电装置。尽管未在图1中示出，但是电池110连接到负载和/或充电装置并且在使用中。也就是说，电池110可以将电流放电到负载或从充电装置充入电力。

电压测量单元120可以在使用中连接到每个电池单体111的两个端子，并且可以测量每个电池单体111的电压且可以周期性地产生每个电池单体111的电压值。因为当电池单体111并联连接时，并联连接的电池单体111具有相同的电压，所以电压测量单元120可以测量一个电池单体111的电压并且可以产生剩余的电池单体111的相同电压值。

例如，电压测量单元120可以在预设时间周期Δt测量每个电池单体111的两个端子处的电压。当前或最近测量的电压值被称为当前电压值，并且由V(t)表示。在时间周期Δt之前测量的电压值被称为紧前电压值，并且由V(t-1)表示。时间周期Δt可以是例如1秒。然而，这是示例性的，时间周期Δt可以被设定为另一时间，例如，0.1秒、0.5秒、2秒、5秒、10秒等。可以根据连接到电池系统100的电动系统适当地设定时间周期Δt。

每个电池单体111的电压值可以根据电池单体111的内部状态而彼此不同。第j电池单体111的当前电压值和紧前电压值分别由V_j(t)和V_j(t-1)表示。这里，j表示电池单体111的序数，可以大于或等于1并且可以小于或等于电池110中包括的电池单体111的总数量。

电流测量单元130可以测量使用中的每个电池单体111的电流，并且可以周期性地产生每个电池单体111的电流值。因为当电池单体111串联连接时，串联连接的电池单体111具有相同的电流，所以电流测量单元130可以测量一个电池单体111的电流，并且可以产生剩余的电池单体111的相同电流值。

电压测量单元120和电流测量单元130可以彼此同步，并且可以同时测量每个电池单体111的电压和电流。电流测量单元130可以在时间周期Δt测量每个电池单体111的电流。由电流测量单元130测量的电流值可以在电流是充电电流时显示为正(+)，并且可以在电流是放电电流时显示为负(-)。当前或最近测量的电流值被称为当前电流值并且由I(t)表示，在时间周期Δt之前测量的电流值被称为紧前电流值并且由I(t-1)表示。

根据电池单体111的内部状态，每个电池单体111的电流值可以彼此不同。第j电池单体111的当前电流值和紧前电流值分别由I_j(t)和I_j(t-1)表示。

微处理器140可以从由电压测量单元120提供的每个电池单体111的电压值和由电流测量单元130提供的每个电池单体111的电流值实时更新G参数值和H参数值，G参数值和H参数值是指示每个电池单体111的当前状态的G参数和H参数中的每个的数值。这里，G参数是指示使用中的电池单体111的电压相对于电流变化的灵敏度的参数，H参数是指示由使用中的电池单体110中的局部平衡电位分布和电阻分布确定的有效电位的参数。第j电池单体111的G参数和H参数分别由G_j和H_j表示。第j电池单体111的G参数值和H参数值由G_j(t)和H_j(t)表示。

微处理器140可以使用自适应滤波器，以便从每个电池单体111的电压值和电流值产生每个电池单体的G参数值和H参数值。自适应滤波器可以是使用递归最小二乘(RLS)的滤波器或者使用加权最小二乘(WLS)的滤波器。在下文中，将详细描述其中微处理器140使用利用RLS的自适应滤波器的实施例。

微处理器140可以通过使用实时产生的每个电池单体111的G参数值和H参数值来确定电池单体111之中的异常劣化的电池单体或异常减少容量的电池单体。可以通过使用自适应滤波器从每个电池单体111的电压值和电流值实时产生每个电池单体111的G参数值和H参数值。异常劣化的电池单体和异常减少容量的电池单体可以被称为异常电池单体。

根据实施例，微处理器140可以从每个电池单体111的G参数值计算代表G参数值。代表G参数值可以是G参数值的平均值或中值。代表G参数值可以被确定为可以代表G参数的不同值。

微处理器140可以从每个电池单体111的H参数值计算代表H参数值。代表H参数值可以是H参数值的平均值或中值。代表H参数值可以被确定为可以代表H参数的不同值。

微处理器140可以通过使用每个电池单体111的G参数值、每个电池单体111的H参数值、代表G参数值和代表H参数值来确定每个电池单体111是异常劣化的电池单体还是异常减少容量的电池单体。

根据实施例，因为微处理器140仅使用对四则算术运算的简单运算来检测异常电池单体，所以微处理器140可以包括在电池组的电池管理系统(BMS)中。根据另一实施例，根据本公开的用于检测异常电池单体的方法可以由电动车辆的BMS中的微控制器或电子控制单元(ECU)执行。根据另一实施例，根据本公开的用于检测异常电池单体的方法可以由能量存储系统的集成控制器执行。根据另一实施例，用于检测异常电池单体的方法可以由通过通信连接到电池系统或能量存储系统的服务器的处理器执行。

存储单元150可以存储执行根据本实施例的用于检测异常电池单体的方法所需的命令和数据。因为，在根据本实施例的方法中，可以基于在时间周期Δt产生的每个电池单体111的电压值和电流值而产生每个电池单体111的G参数值和H参数值，并且基于每个电池单体111的G参数值和H参数值检测异常电池单体，所以每个电池单体111的当前电压值、当前电流值以及紧前值可以存储在存储单元150中，并且其他的电压数据和电流数据可以不存储在存储单元150中。也就是说，不需要将大量的电压数据和电流数据存储在存储单元150中。

根据根据本实施例的用于检测异常电池单体的方法，存储单元150可以存储包括每个电池单体111的G参数和H参数的状态向量的紧前值和当前值以及RLS方法所需的协方差矩阵的紧前值和当前值。此外，存储单元150可以存储代表G参数值和代表H参数值。因此，因为不需要将大量命令和数据存储在存储单元150中，所以可以将存储单元150实现为具有小尺寸的存储器。例如，存储单元150可以被实现为微处理器140内的存储器。

代表电池的当前状态的重要估算值包括荷电状态(SOC)、健康状态(SOH)和功率能力状态(SOP)。在根据相关领域的用于估算电池的内部状态的方法中，预先实验地找到可测量变量(电流I、电压V和温度T)与待估算的内部状态之间的相关性，然后将相关性制成表的形式，并且使用制成表的形式的相关性来估算电池的内部状态，诸如电池单体的内部电阻、内部电容和开路电压(OCV)。这些用于估算电池状态的方法在成本、可扩展性和适应性方面有不足。

根据根据相关领域的用于估算电池的内部状态的方法，需要通过实验预先收集数据，并且基于所收集的数据来估算电池的内部状态。因为内部状态受各种变量影响，所以为了精确估算内部状态，需要通过配置电池的荷电状态(SOC)、电池的电流I、电压V和温度T的各种组合来收集数据。因此，追求的精度越高，需要的数据越多，因此，时间和财务成本也增加。另外，内部状态需要通过对超过由实验测量的变量范围的值进行外推来估算。在这种情况下，估算的可靠性会是一个问题。另外，当通过实验估算的内部状态不准确时，估算的内部状态的效用降低。预先实验地找到相关性并使用相关性来估算状态变量的方法也具有类似的问题。

为了提取电池组的状态信息，需要关于电池组的实验数据。然而，因为收集关于电池组的数据需要高成本，所以通常替代地使用从电池单体测量的数据来估算电池组的状态信息。然而，该方法基于以下假设：从电池单体测量的变量之间的关系是对更大系统的可扩展机制。因此，如果没有通过附加验证预先验证系统之间的机制的可扩展性，则通过该方法估算的电池组的状态信息的可靠性是有疑问的，并且会存在准确性的问题。

根据相关领域的内部状态估算的适应性是降低的。电池随时间具有不同的特性。这被称为老化效应。可测量变量(诸如电池的电压V、电流I和温度T)与内部状态之间的相关性也随时间而变化。然而，如果通过实验预先完成指示可测量变量与待估算的内部状态之间的相关性的表，则难以反映老化效应。为了反映老化效应，需要预先执行关于老化效应的各种实验，或者需要预先预测老化效应以在实验模型中反映。然而，前者具有成本增加的问题，后者具有如果没有预先准确地预测老化效应则降低状态估算的可靠性的风险。电池单体的设计以及时间的变化也会导致电池单体特性的变化。然而，如果这些变化未被预先预测并且未反映在实验模型上，则当存在诸如在电池单体的制造中使用的活性物质或形状的设计变化时，不能使用现有的实验模型。

除了实验模型之外，还可以使用考虑电池中的电化学现象和热行为的电化学模型(例如，Newman模型)来估算电池内部状态。通过使用该模型，不能直接测量的物理状态可以用作电池管理系统的控制参数。然而，这种模型的缺点在于它需要过多的时间和成本来开发，并且计算所需的资源太大。因此，它在电池管理系统中没有广泛使用。

为了克服根据相关领域的方法的问题，本公开提出了一种通过使用G参数和H参数来区分正常电池单体与异常电池单体的方法，G参数和H参数是指示电池单体的当前状态的参数。根据本公开的用于检测异常电池单体的方法可以相对简单地实现以在电池管理系统(BMS)中执行，并且可以在没有附加驱动条件的情况下具有高精度。

G参数是指示电池单体的端子电压对施加到使用中的电池单体的电流的变化的灵敏度的状态量，并且具有电阻的单位。H参数是在使用期间由电池单体内的局部平衡电位分布和电阻分布确定的有效电位。可以使用理论模型将电池单体的G参数和H参数量化为电池材料性质与设计变量之间的显式相关性。在下文中，将描述电池单体的G参数和H参数。

在电池单体中，可以假设电压V和电流i具有与V＝f(i；x,p)相同的关系。这里，x是指示电池单体的内部状态的物理量，并且p是参数。

函数f是非线性隐式函数。当函数f可以被划分为快速变化量g和缓慢变化量h时，上述关系为V＝g(i；x,p)+h(i；x,p)。

如果假设存在相对于电流i缓慢变化的函数G(i；x,p)＝dg/di，则上述关系可以表示为V＝G(i；x,p)i+H(i；x,p)。

在上述关系中，dG/di和dH/di具有非常小的值。换句话说，当满足上述假设时，G和H是相对于电流i缓慢变化的函数。因此，表示电压V与电流i之间的非线性关系的函数f可以表示为类似于上述关系的准线性关系。

这里，G被称为G参数，H被称为H参数。当电流i是充电/放电电流并且Ueq是电池单体的平衡电位时，可以使用G参数G和H参数H将放电过电压计算为Ueq-V＝-G·i+(Ueq-H)。

这里，是由电池单体产生的使电流流过端子的过电压，并且包括反应动态极化量以及电子和离子电阻极化量。(Ueq-H)是由电池单体的局部热力学平衡状态与整个系统的平衡状态的偏差引起的过电压。也就是说，(Ueq-H)是由电池单体内部的热力学不均匀性引起的低效率。当电池单体的内部系统达到热力学平衡状态时，H参数H变得等于平衡电位Ueq。

在根据实施例的用于检测异常电池单体的方法中，通过使用例如RLS方法从由每个电池单体测量的电压值和电流值直接提取每个电池单体的G参数和H参数，并且使用每个电池单体的G参数G和H参数H检测异常电池单体。

图2是示出根据实施例的用于检测异常电池单体的方法的流程图。

参照图2连同图1，微处理器140可以通过使用自适应滤波器来执行图2中所示的用于检测异常电池单体的方法。

微处理器140可以通过使用电压测量单元120和电流测量单元130来检测电连接在第一端子101与第二端子102之间并且在使用中的多个电池单体111中的每个的电压和电流。微处理器140可以周期性地产生每个电池单体111的电压值V和电流值I(S110)。可以产生第j电池单体(例如，Cell_j)的电压值V_j和电流值I_j。第j电池单体(例如，Cell_j)可以指电池110中包括的电池单体111之中的一个。

微处理器140可以通过使用存储在存储单元150中的自适应滤波器实时地从每个电池单体(例如，Cell_j)的电压值V_j和电流值I_j更新G参数值G_j(t)和H参数值H_j(t)，G参数值G_j(t)和H参数值H_j(t)是指示每个电池单体(例如，Cell_j)的当前状态的G参数和H参数中的每个的数值(S120)。下面将详细描述实时更新G参数值G_j(t)和H参数值H_j(t)的操作(S120)。

微处理器140可以计算代表电池单体111的G参数值G(t)的代表G参数值G_ref，并且可以计算代表电池单体111的H参数值H(t)的代表H参数值H_ref(S130)。根据实施例，代表G参数值G_ref可以是G参数值G(t)的平均值，并且代表H参数值H_ref可以是H参数值H(t)的平均值。根据另一实施例，代表G参数值G_ref可以是G参数值G(t)的中值，并且代表H参数值H_ref可以是H参数值H(t)的中值。

微处理器140可以基于G参数值G_j(t)、H参数值H_j(t)、代表G参数值G_ref和代表H参数值H_ref来确定每个电池单体(例如，Cell_j)是否是异常劣化的电池单体。例如，微处理器140可以相对于每个电池单体(例如，Cell_j)通过从每个电池单体(例如，Cell_j)的G参数值G_j(t)中减去代表G参数值G_ref来计算每个电池单体(例如，Cell_j)的G参数偏差值G_j(t)-G_ref，可以对于每个电池单体(例如，Cell_j)通过从每个电池单体(例如，Cell_j)的H参数值H_j(t)中减去代表H参数值H_ref来计算每个电池单体(例如，Cell_j)的H参数偏差值H_j(t)-H_ref，并且可以确定每个电池单体(例如，Cell_j)的G参数偏差值G_j(t)-G_ref是否大于预设G参数设定值(G_th)，或者确定每个电池单体(例如，Cell_j)的H参数偏差值H_j(t)-H_ref是否大于预设H参数设定值H_th(S140)。

当每个电池单体(例如，Cell_j)的G参数偏差值G_j(t)-G_ref大于预设G参数设定值G_th或者每个电池单体(例如，Cell_j)的H参数偏差值H_j(t)-H_ref大于预设H参数设定值H_th时，微处理器140可以将对应的电池单体(例如，Cell_j)确定为异常劣化的电池单体(S150)。例如，G参数设定值G_th可以被预设为约0.5，H参数设定值H_th可以被预设为约0.01。然而，这些值是示例性的，可以根据电池单体111的特性将G参数设定值G_th和H参数设定值H_th设定为不同值。

如果每个电池单体(例如，Cell_j)的G参数偏差值G_j(t)-G_ref不大于预设G参数设定值G_th并且每个电池单体(例如，Cell_j)的H参数偏差值H_j(t)-H_ref不大于预设H参数设定值H_th，则微处理器140可以通过对电池单体(例如，Cell_k)执行操作S140的确定来确定电池单体Cell_k是否是异常劣化的电池单体。

如果所有的电池单体111的G参数偏差值G(t)-G_ref都不大于预设G参数设定值G_th并且所有的电池单体111的H参数偏差值H(t)-H_ref都不大于预设H参数设定值H_th，则微处理器140可以进行到操作S110以在预设时间周期Δt之后新产生每个电池单体111的电压值和电流值。

当微处理器140在操作S150中将电池单体(例如，Cell_j)确定为异常劣化的电池单体时，微处理器140可以确定电池单体(例如，Cell_j)的G参数偏差值G_j(t)-G_ref是否大于预设G参数设定值G_th，并且确定电池单体(例如，Cell_j)的H参数偏差值H_j(t)-H_ref是否大于预设H参数设定值H_th(S160)。

当电池单体(例如，Cell_j)的G参数偏差值G_j(t)-G_ref大于预设G参数设定值G_th并且电池单体(例如，Cell_j)的H参数偏差值H_j(t)-H_ref大于预设H参数设定值H_th时，微处理器140可以确定对应的电池单体Cell_j是异常减少容量的电池单体(S170)。

当在操作S160中确定电池单体(例如，Cell_j)的G参数偏差值G_j(t)-G_ref不大于预设G参数设定值G_th或者电池单体(例如，Cell_j)的H参数偏差值H_j(t)-H_ref不大于预设H参数设定值H_th时，以及当在操作S170中确定电池单体Cell_j为异常减少容量的电池单体时，微处理器140可以进行到操作S110以在预设时间周期Δt之后新产生每个电池单体111的电压值和电流值。

根据实施例，微处理器140可以通过使用存储在存储单元150中的使用RLS方法的滤波器来执行图2的操作。

图3是示出图2的用于检测异常电池单体的方法的操作S110至S120的流程图。

参照图1至图3，当微处理器140使用RLS滤波器时，微处理器140可以初始化每个电池单体111的状态向量Θ＝[G；H]和协方差矩阵P(S101)。例如，可以存在与第j电池单体Cell_j对应的状态向量Θ_j＝[G_j；H_j]和协方差矩阵P_j，并且在操作S101中，微处理器140可以初始化第j电池单体(例如，Cell_j)的状态向量Θ_j＝[G_j；H_j]和协方差矩阵P_j。

图2的操作S110可以包括图3的操作S111、S112和S113。在操作S111中，微处理器140可以产生每个电池单体111的当前电压值V(t)和当前电流值I(t)。在时间周期Δt之前产生的每个电池单体111的当前电压值V(t)和当前电流值I(t)可以是每个电池单体111的紧前电压值V(t-1)和紧前电流值I(t-1)。例如，微处理器140可以产生第j电池单体(例如，Cell_j)的紧前电压值V_j(t-1)和紧前电流值I_j(t-1)，并且可以在时间周期Δt之后产生第j电池单体(例如，Cell_j)的当前电压值V_j(t)和当前电流值I_j(t)。

在操作S112中，微处理器140可以计算对于每个电池单体111的当前电流值I(t)与紧前电流值I(t-1)之间的电流差ΔI，例如，对于第j电池单体(例如，Cell_j)的当前电流值I_j(t)与紧前电流值I_j(t-1)之间的电流差ΔI_j。

在操作S113中，微处理器140可以将每个电池单体111的电流差ΔI与预设范围进行比较。可以通过使用下限值ΔImin和上限值ΔImax来定义预设范围。例如，微处理器140可以确定第j电池单体(例如，Cell_j)的电流差ΔI_j是否在下限值ΔImin与上限值ΔImax之间。

当第j电池单体(例如，Cell_j)的电流差ΔI_j在下限值ΔImin与上限值ΔImax之间时，在操作S121中，可以更新第j电池单体(例如，Cell_j)的状态向量Θ_j和协方差矩阵P_j。因为第j电池单体(例如，Cell_j)的状态向量Θ_j包括第j电池单体(例如，Cell_j)的G参数值G_j(t)和H参数值H_j(t)，所以当在操作S121中更新状态向量Θ_j时，可以更新G参数值G_j(t)和H参数值H_j(t)(S120)。

随后，微处理器140可以执行图2的操作S130。

根据本实施例，可以使用状态向量和协方差矩阵P(t)。在下文中，将描述实时更新电池单体111之中的一个电池单体(例如，Cell_j)的G参数值G_j(t)和H参数值H_j(t)的方法。因为G参数值G_k(t)和H参数值H_k(t)对于其他电池单体(例如，Cell_j)以相同的方式计算，所以省略了分别标识电池单体111的序数(例如，j或k)。

状态向量可以包括G参数和H参数并且可以被定义为

协方差矩阵P(t)可以包括第一值P₁(t)和第二值P₂(t)，并且可以被定义为

可以对每个时间周期Δt产生电池单体111的电压值V(t)和电流值I(t)，并且可以通过使用递归方法对每个时间周期Δt更新状态向量和协方差矩阵P(t)。也可以对每个时间周期Δt更新G参数和H参数

在操作S101中，微处理器140可以初始化状态向量和协方差矩阵P(t)。例如，状态向量可以被初始化为并且协方差矩阵P(t)可以被初始化为在本实施例中，状态向量的元素和协方差矩阵P(t)的元素都被初始化为1，但是这是示例性的，并且可以被初始化为其他值。

例如，微处理器140可以通过使用电压测量单元120和电流测量单元130来测量每个电池单体111的电压和电流，并且可以对每个时间周期Δt周期性地产生每个电池单体111的电压值和电流值(S120)。微处理器140可以产生每个电池单体111的紧前电压值V(t-1)和紧前电流值I(t-1)，并且可以在时间周期Δt之后产生每个电池单体111的当前电压值V(t)和当前电流值I(t)。

微处理器140可以通过使用RLS方法通过使用对每个时间周期Δt产生的当前电压值V(t)和当前电流值I(t)来对每个时间周期Δt更新G参数和H参数

微处理器140可以基于当前电流值I(t)和状态向量的紧前值来计算电池单体111的当前电压估算状态向量的紧前值可以基于在时间周期Δt之前的紧前电压值V(t-1)和紧前电流值I(t-1)来计算，并且可以包括G参数的紧前值和H参数的紧前值作为通过将H参数的紧前值与当前电流值I(t)和G参数的紧前值的乘积相加而获得的值的电池单体111的当前电压估算可以被计算为

微处理器140可以基于协方差矩阵的当前电流值I(t)和紧前值P(t-1)来更新增益矩阵L(t)。根据协方差矩阵P(t)的定义，协方差矩阵的紧前值P(t-1)可以基于在时间周期Δt之前的紧前电压值V(t-1)和紧前电流值I(t-1)来计算，并且可以包括协方差矩阵的第一紧前值P₁(t-1)和协方差矩阵的第二紧前值P₂(t-1)。增益矩阵L(t)对应于每个电池单体111。

增益矩阵L(t)用于更新状态向量和协方差矩阵P(t)。增益矩阵L(t)可以包括增益矩阵的第一值L₁(t)和增益矩阵的第二值L₂(t)，并且可以如下计算。

这里，λ₁是第一遗忘因子，并且与G参数相关。λ₂是第二遗忘因子，并且与H参数相关。第一遗忘因子λ₁和第二遗忘因子λ₂分别是用于指示在G参数和H参数的计算中过去的电压值和过去的电流值对G参数的当前值和H参数的当前值的影响的值。第一遗忘因子λ₁和第二遗忘因子λ₂可以随着它们更接近1而在长时间段内具有对G参数的当前值和H参数的当前值的影响，并且可以随着它们更接近0而在短时间段内具有对G参数的当前值和H参数的当前值的影响。

根据实施例，第一遗忘因子λ₁和第二遗忘因子λ₂可以是0.9以上且1以下。根据另一实施例，第一遗忘因子λ₁可以被设定为大于或等于第二遗忘因子λ₂的值。例如，第一遗忘因子λ₁可以被设定为0.99999，并且第二遗忘因子λ₂可以被设定为0.95。这些设定值可以根据电池单体111的特性而改变。

本公开的发明人发现，在对特定电池单体执行的实验中，当第一遗忘因子λ₁和第二遗忘因子λ₂分别为0.99999和0.95时，获得了高可靠性结果。然而，上述值是示例性的，也可以根据电池单体111的特性设定为不同的值。例如，第一遗忘因子λ₁可以被设定为0.9999，并且第二遗忘因子λ₂也可以被设定为0.98。

根据另一实施例，第一遗忘因子λ₁和第二遗忘因子λ₂中的全部也可以被设定为1。在这种情况下，可以看出未应用第一遗忘因子λ₁和第二遗忘因子λ₂。

微处理器140可以基于当前电流值I(t)、增益矩阵L(t)和协方差矩阵的紧前值P(t-1)来更新协方差矩阵P(t)。协方差矩阵P(t)可以如下计算。

微处理器140可以将当前电压值V(t)与上述计算的当前电压估算之间的电压误差e(t)计算为

微处理器140可以基于状态向量的紧前值增益矩阵的当前值L(t)和电压误差e(t)来更新状态向量可以更新状态向量使得可以产生G参数的当前值和H参数的当前值

状态向量的当前值可以如下计算为通过将增益矩阵的当前值L(t)和电压误差e(t)的乘积与状态向量的紧前值相加而获得的值。

用于递归地表达状态向量的上述等式可以如下导出。

首先，可以如下定义应用第一遗忘因子λ₁和第二遗忘因子λ₂的损失函数ε。

这里，V(i)是第i电压值，I(i)是第i电流值。V(t)和I(t)分别是当前电压值和当前电流值，并且V(t-1)和I(t-1)可以分别是紧前电压值和紧前电流值。

G(i)和H(i)分别是第i G参数和第i H参数的实数值，并且和分别是G参数的当前值估算和H参数的当前值估算。

当分别相对于和微分损失函数ε的结果为0时，相对于和的损失函数ε被最小化。

其中相对于微分损失函数ε的结果为0的可以如下获得。

重新排列上述等式，如下。

其中相对于微分损失函数ε的结果为0的可以如下获得。

重新排列上述等式，如下。

对于实时估算，上述获得的和使用状态向量以递归形式总结如下。

如上所述，因为当前电压估算被定义为并且电压误差e(t)被定义为所以状态向量可以被表达如下。

这里，分别如上所述，增益矩阵L(t)和协方差矩阵P(t)可以如下计算。

因为根据本实施例的用于检测异常电池单体的方法使用递归方法，所以存储单元150可以存储当前电压值V(t)、当前电流值I(t)、状态向量和协方差矩阵P(t)。根据另一实施例，存储单元150还可以存储第一遗忘因子λ₁和第二遗忘因子λ₂。存储单元150不需要存储过去的电压值和过去的电流值中的全部。

因为图2中所示的用于检测异常电池单体的方法使用递归方法，所以计算非常简单，并且可以通过使用具有几kB的小尺寸的存储单元150来执行计算。此外，因为每当接收到电压值和电流值时新更新状态向量和协方差矩阵P(t)，所以每个电池单体111的电压变化和电流变化可以基本上实时地反映在G参数和H参数中。

微处理器140可以对测量每个电池单体111的电压值V(t)和电流值I(t)的每个时间周期Δt重复执行操作S120至S170。因此，当电池单体111之中发生异常劣化的电池单体和异常减少容量的电池单体时，微处理器140可以实时检测到它。

根据图3的实施例，G参数和H参数可以不对每个时间周期Δt更新。

如上所述，微处理器140可以通过将电流差ΔI与其中将电流差ΔI定义为下限值ΔImin和上限值ΔImax的特定范围进行比较来确定电流差ΔI是否包括在特定范围内(S113)。可以根据电池系统100的充电/放电电流模式来设定下限值ΔImin和上限值ΔImax。下限值ΔImin和上限值ΔImax可以基于电池系统100的充电/放电速率来设定。例如，下限值ΔImin可以被设定为0.2C，并且上限值ΔImax可以被设定为1C。

当电流差ΔI太小或太大时，由于递归方法，电压值与电流值之间的感测误差导致G参数和H参数的大偏差。当电流值ΔI包括在特定范围内时，该方法可以进行到操作S121，使得可以基于当前电压值V(t)和当前电流值I(t)来更新状态向量和协方差矩阵P(t)。然而，当电流差ΔI超过特定范围时，该方法可以进行到操作S111，使得可以不更新状态向量和协方差矩阵P(t)。也就是说，当电流差ΔI超过特定范围时，当前电压值V(t)和当前电流值I(t)不会影响状态向量和协方差矩阵P(t)。

只有当电流差ΔI包括在特定范围内时，才会更新G参数和H参数使得G参数和H参数的可靠性可以保持为高的。

下面将参照图4至图6描述测量或估算正常电池单体、异常劣化的电池单体和异常减少容量的电池单体的参数的结果。

图4示出了在左侧的在包括正常电池单体的电池组中测量的电池组电流(a)和电池组电压(b)以及针对电池组估算的G参数值(c)和H参数值(d)以及在右侧的在包括在电池组中的每个正常电池单体中测量的电池单体电压(e)和电池单体SOC(f)以及针对电池单体估算的G参数值(g)和H参数值(h)。

图4的曲线图a-h使用20Ah的能量存储系统(ESS)的电流/电压实际使用的操作模式数据。使用其中四个电池单体串联连接的结构，并且时间周期Δt为1秒。流过正常电池单体的电流与电池组的电流相同。图4的曲线图a-h示出了特定参数值根据时间(秒)的变化。

图5示出了在用图4的电池组中的异常劣化的电池单体替换一个正常电池单体之后测量或估算一个正常电池单体的结果。图5的曲线图a示出了正常电池单体和异常劣化的电池单体的G参数值，并且图5的曲线图b示出了正常电池单体和异常劣化的电池单体的电池单体电压。

通过曲线图b的电池单体电压不能区分正常电池单体和异常劣化的电池单体，但是通过使用曲线图a的G参数值清楚地区分正常电池单体和异常劣化的电池单体。可以看出，异常劣化的电池单体的G参数值大于正常电池单体的G参数值。这样，根据本公开的实施例，将特定电池单体的G参数值与其他电池单体的G参数值进行比较，使得可以清楚地看出特定电池单体是否是异常劣化的电池单体。

当电池单体劣化时，电池单体的内部电阻增加。因为G参数表示电压相对于电池单体的电流变化的灵敏度，所以可以直接反映电池单体的内部电阻。因此，基于G参数值确定异常劣化的电池单体的方法可以具有可靠性。

图6示出了在用图4的电池组中的异常减少容量的电池单体替换一个正常电池单体之后测量或估算一个正常电池单体的结果。图6的曲线图a示出了正常电池单体和异常减少容量的电池单体的电池单体电压，并且图6的曲线图b示出了正常电池单体和异常减少容量的电池单体的电池单体SOC，图6的曲线图c示出了正常电池单体和异常减少容量的电池单体的G参数值，并且图6的曲线图d示出了正常电池单体和异常减少容量的电池单体的H参数值。

参照曲线图c和d，可以看出，不仅异常减少容量的电池单体的G参数值大于正常电池单体的G参数值，而且异常减少容量的电池单体的H参数值也大于正常电池单体的H参数值。因此，根据本公开的实施例，将特定电池单体的G参数值和H参数值分别与其他电池单体的G参数值和H参数值进行比较，使得可以清楚地看出特定电池单体是否是异常减少容量的电池单体。

本公开的精神不限于上述实施例并且不应被限定，并且权利要求以及等同于这些权利要求或从这些权利要求改变的所有范围都落入本公开的范围内。