用于可再充电电池容量的实时估计的系统和方法与流程

用于可再充电电池容量的实时估计的系统和方法
1.本申请是申请日为2016年8月2日、国际申请号为pct/us2016/045084、国家申请号为201680082627.3、发明名称为“用于可再充电电池容量的实时估计的系统和方法”的发明专利申请的分案申请。
技术领域
2.本公开总体上涉及电池和电池模块的领域。具体来说，本公开涉及估计可再充电电池的实时容量。


背景技术：

3.本节旨在向读者介绍可能涉及下述的本公开的各个方面的本领域中的各个方面。据信，此论述有助于为读者提供背景信息，以便于更好地理解本发明的各个方面。因此，应理解，这些表述应当从这个角度来阅读，而不是作为对现有技术的承认。
4.使用一个或多个电池系统来为车辆提供全部或一部分动力的车辆能够被称为xev，其中术语“xev”在本文中被定义成包括所有以下车辆或者其任何变型或组合：其将电力用作其车辆动力的全部或一部分。例如，xev包括使用电力作为全部动力的电动车辆(ev)。如本领域中的技术人员应了解，也被视作xev的混合动力电动车辆(hev)将内燃机推进系统与电池供能的电动推进系统(例如48伏(v)或130v系统)相组合。术语hev可以包括混合动力电动车辆的任何变型。例如，全混合动力系统(fhev)可以使用一个或多个电动机、仅使用内燃机，或者使用这两者向车辆提供动力和其他电力。相比之下，轻度混合动力系统(hmev)在车辆怠速时停用内燃机，并且使用电池系统来继续向空气调节单元、无线电或者其他电子设备供电，以及在需要推进时重新起动发动机。轻度混合动力系统还可以应用一定程度的动力辅助，例如在加速期间，以作为对内燃机的补充。轻度混合动力通常在96v到130v，并且通过皮带或曲柄集成式起动器发电机回收制动能量。此外，微型混合动力电动车辆(mhev)也使用类似于轻度混合动力的“启
‑
停”系统，但微型混合动力系统可以向内燃机提供动力辅助或不提供动力辅助，并且在60v以下的的电压下操作。为了本论述的目的，应认识到，mhev通常在技术上不使用直接提供到曲轴或传动系的电力作为车辆的任何部分的动力，但是mhev仍然可以视作xev，因为它在车辆怠速(其中内燃机停用)时的确使用电力来作为对车辆动力需求的补充，并且通过集成式起动器发电机来回收制动能量。此外，插电式电动车辆(pev)是能够从外部电源，例如壁式插座充电的任何车辆，并且储存在可再充电电池组中的能量驱动或者帮助驱动车轮。pev是ev的一个子类，包括全电动或电池电动车辆(bev)、插电式混合动力电动车辆(phev)以及混合动力电动车辆和传统内燃机车辆的改装型电动车辆。
5.相对于仅使用内燃机和传统电力系统(通常是由铅酸电池供能的12v系统)的更传统的气动车辆而言，上述xev可提供若干优势。例如，相对于传统内燃车辆，xev可以产生更少的不良排放产物，并且可以呈现更高的燃料效率，并且在一些情况下，此类xev可以完全不使用汽油，如同某些类型的ev或pev一样。
6.随着技术的不断发展，需要为这种车辆的电池模块提供改进的状态指示器。例如，xev使用的电力可以由可再充电电池提供。在可再充电电池工作时，可能难以准确地描绘可再充电电池的充电状态或容量。本公开总体涉及估计在可再充电电池和/或xev的工作期间可再充电电池的实时参数。


技术实现要素：

7.下文概述了本文中所公开的特定实施例。应了解，这些方面仅向读者提供关于这些特定实施例的概述，并且这些方面并不旨在限制本公开的范围。实际上，本公开可以涵盖下文未描述的多个方面。
8.本公开涉及一种电池系统，其包括耦接到电气系统的电池。电池系统还包括电耦接到电池的电池控制模块。电池控制模块监测电池的至少一个监测参数，且电池控制模块使用两个线性回归模型递归地计算所述电池的第一容量估计值。两个线性回归模型基于至少一个等效电路模型、至少一个监测参数和卡尔曼滤波器。
9.本公开还涉及一种确定耦接到电气系统的可再充电电池的实时容量的方法。该方法包括通过耦接到可再充电电池的一个或多个传感器，利用电池控制模块监测电池工作期间可再充电电池的至少一个监测参数。该方法还包括使用两个线性回归模型，通过电池控制模块递归地计算可再充电电池的第一实时容量。两个线性回归模型基于至少一个等效电路模型、至少一个监测参数和卡尔曼滤波器。
10.本公开还涉及一种用于车辆的能量存储部件。能量存储部件包括壳体、第一端子和第二端子、以及设置在壳体中的可再充电电池。可再充电电池耦接到第一端子和第二端子。能量存储部件还包括电池控制模块，其监测所述能量存储部件的至少一个监测参数。电池控制模块还使用基于至少一个等效电路模型、至少一个监测参数和卡尔曼滤波器的两个线性回归模型，递归地计算能量储存部件的第一容量估计值。另外，电池控制模块至少部分地基于能量存储部件的两个弛豫开路电压测量值和能量存储部件的库仑计数，递归地计算能量存储部件的第二容量估计值。
附图说明
11.通过阅读以下详细描述并参考附图，可以更好地理解本公开的各个方面，其中：
12.图1是根据本方法的一个实施例的具有对车辆全部或部分电力作出贡献的电池系统的车辆(xev)的立体图；
13.图2是根据本方法的一个实施例的混合动力电动车辆(hev)形式的图1的xev的剖面示意图；
14.图3是根据本方法的一个实施例的图1的xev的电池系统的示意图；
15.图4是根据本方法的一个实施例的图1的xev的能量存储部件的1
‑
rc等效电路模型；
16.图5是根据本方法的一个实施例的图1的xev的能量存储部件的开路电压(ocv)和充电状态(soc)之间关系的图表；
17.图6a和6b是描述根据本方法的实施例的用于计算能量存储部件参数和确定能量存储部件的参数收敛性的方法的过程流程图；
18.图7a和7b是根据本方法的一个实施例的用于存储数据的环形缓冲器，该数据是使用图6a和6b的过程流程图计算出的；
19.图8a和8b是描述根据本方法的一个实施例的使用两个线性回归模型计算能量存储部件的容量的方法的过程流程图；
20.图9a和9b是描述根据本方法的一个实施例的使用两个弛豫开路电压测量值和库仑计数来计算能量存储部件的容量的方法的过程流程图；
21.图10是描述根据一实施例的用于对能量存储部件的容量估计进行方向验证的方法的过程流程图；以及
22.图11是根据一实施例来描绘图10的过程流程图的实施例的图表。
具体实施方式
23.下文将描述一个或多个具体实施例。为提供对这些实施例的简要描述，本说明书中并未描述实际实施方式的所有特征。应了解到，在任何此类实际实施方式的开发中，如在任何工程或设计项目中，必须做出许多专门针对实施方式的决策，以实现开发人员的具体目标，例如符合系统相关和业务相关的约束条件，这些约束条件可能会随着实施方式的不同而有所不同。此外，应了解到，此类开发工作可能较为复杂并且耗时，但是尽管如此，这也应当是受益于本发明的普通技术人员在设计、制作和制造中的常规任务。
24.本文描述的电池系统可用于向各种类型的电动车辆(xev)和其他高电压能量存储/消耗应用(例如，电网电力存储系统)提供电力。这种电池系统可以包括一个或多个电池模块，每个电池模块具有多个电池单元(例如，锂离子(li
‑
ion)电化学电池单元)，其布置并电互连以提供有益于对例如xev的一个或多个部件供电的特定电压和/或电流。作为另一个例子，根据本实施例的电池模块可以与固定电力系统(例如，非汽车系统)结合或向固定电力系统提供电力。
25.基于优于传统燃气动力车辆的优点，通常生产传统气动车辆的制造商可能希望在其车辆生产线内使用改进的车辆技术(例如，再生制动技术)。通常，这些制造商可以利用其传统车辆平台之一作为起点。由于传统的气动车辆设计为使用12伏电池系统，因此可以使用12伏锂离子电池来补充12伏铅酸电池。更具体地，12伏锂离子电池可用于更有效地捕获再生制动期间产生的电能，并随后供应电能以为车辆的电气系统提供供电。
26.随着车辆技术的进步，高电压电气设备也可能包含在车辆的电气系统中。例如，锂离子电池可以向轻度混合动力电动车辆中的电动机供应电能。通常，这些高电压电气设备使用大于12伏的电压，例如高达48伏的电压。因此，在一些实施例中，可以使用dc
‑
dc转换器来对12伏锂离子电池的输出电压进行升压，以向高电压装置供电。另外或可替代地，48伏锂离子电池可用于补充12伏铅酸电池。更具体地，48伏锂离子电池可用于更有效地捕获再生制动期间产生的电能，并随后提供电能以为高电压装置供电。
27.因此，关于是使用12伏锂离子电池还是48伏锂离子电池的设计选择可直接取决于特定车辆中包括的电气装置。然而，尽管电压特性可能不同，但12伏锂离子电池和48伏锂离子电池的操作原理大体是相似的。更具体地，如上所述，两者都可用于在再生制动期间捕获电能并随后供应电能以为车辆中的电气设备供电。
28.因此，为了简化以下讨论，将结合具有12伏锂离子电池和12伏铅酸电池的电池系
统来描述本技术。然而，本领域普通技术人员能够使本技术适用于其他电池系统，例如具有48伏锂离子电池和12伏铅酸电池的电池系统。
29.本公开涉及电池和电池模块。更具体地，本公开涉及估计可再充电电池的实时参数。具体实施例涉及可用于车辆环境(例如，混合动力电动车辆)以及其他能量存储/消耗应用(例如，用于电网的能量存储)的锂离子电池单元。
30.考虑到前述内容，本公开描述了用于估计可再充电电池的实时参数的技术。传统上，可再充电电池的额定容量可以通过在室温下使用恒定放电速率将可再充电电池从充满电状态完全放电来确定。这种方法对于车辆环境中使用的可再充电电池可能是不切实际的。因此，可使用可再充电电池的部分放电来估计可再充电电池的容量。因此，本公开中描述的电池控制单元可以使用下面详细描述的系统和方法实时或接近实时地估计可再充电电池的容量和其他参数。
31.为了帮助说明，图1是可利用再生制动系统的车辆10的一实施例的立体图。尽管以下讨论涉及具有再生制动系统的车辆，但是本文描述的技术可适用于利用电池捕获/存储电能的其他车辆，其可包括电动和气动车辆。
32.如上所述，期望电池系统12在很大程度上与传统的车辆设计兼容。因此，电池系统12可以布置在车辆10中本应容纳传统电池系统的位置。例如，如图所示，车辆10可包括电池系统12，其所在的位置类似于典型的内燃机车辆的铅酸电池的位置(例如，在车辆10的引擎盖下)。此外，如下面将更详细描述的，电池系统12可以定位成便于管理电池系统12的温度。例如，在一些实施例中，将电池系统12定位在车辆10的引擎盖下方可以使空气管道能够使空气流流过电池系统12上方并冷却电池系统12。
33.图2中描述了电池系统12的更详细视图。如图所示，电池系统12包括能量存储部件14，其耦接到点火系统16、交流发电机18、车辆控制台20，并且可选地耦接到电动机22。通常，能量存储部件14可以捕获/存储在车辆10中产生的电能并输出电能以为车辆10中的电气设备供电。
34.换句话说，电池系统12可以向车辆电气系统的部件提供电力，这些部件可以包括散热器冷却风扇、气温控制系统、电动助力转向系统、活动悬架系统、自动泊车系统、电动油泵、电动增压器/涡轮增压器、电动水泵、加热挡风玻璃/除雾器、车窗升降电机、梳妆灯、胎压监测系统、天窗电机控制器、电动座椅、警报系统、信息娱乐系统、导航特征、车道偏离警告系统、电动驻车制动器、外部灯或者它们的任何组合。说明性地，在图示的实施例中，能量存储部件14向车辆控制台20、车辆内的显示器21和点火系统16供电，所述点火系统可用于起动(例如，经曲柄起动)内燃机24。
35.此外，能量存储部件14可以捕获交流发电机18和/或电动机22产生的电能。在一些实施例中，交流发电机18可以在内燃机24运行时产生电能。具体来说，交流发电机18可以将内燃机24旋转所产生的机械能转换成电能。另外或可替代地，当车辆10包括电动机22时，电动机22可以通过将车辆10运动(例如，车轮旋转)所产生的机械能转换成电能来产生电能。因此，在一些实施例中，能量存储部件14可以捕获交流发电机18和/或电动机22在再生制动期间产生的电能。因此，交流发电机18和/或电动机22在本文中被统称为再生制动系统。
36.为便于捕获和供应电能，能量存储部件14可以经由总线26电耦接到车辆的电气系统。例如，总线26可以使得能量存储部件14能够接收交流发电机18和/或电动机22产生的电
能。此外，总线26可以使得能量存储部件14向点火系统16和/或车辆控制台20输出电能。因此，当使用12伏电池系统12时，总线26可以承载通常在8
‑
18伏之间的电功率。
37.另外，如图所示，能量存储部件14可包括多个电池模块。例如，在图示的实施例中，能量存储部件14包括根据本实施例的铅酸(例如，第一)电池模块28、以及锂离子(例如，第二)电池模块30，其中每个电池模块28、30包括一个或多个电池单元。在其他实施例中，能量存储部件14可包括任意数量的电池模块。另外，尽管第一电池模块28和第二电池模块30在图示中为彼此相邻，但是它们可以定位在车辆周围的不同区域中。例如，第二电池模块30可以定位在车辆10的内部或周围，而第一电池模块28可以定位在车辆10的引擎盖下方。
38.在一些实施例中，能量存储部件14可包括多个电池模块以利用多种不同的电池化学组成。例如，第一电池模块28可以使用铅酸电池化学组成，第二电池模块30可以使用锂离子电池化学组成。在这样的实施例中，电池系统12的性能可以得到改善，因为锂离子电池化学组成通常具有比铅酸电池化学组成更高的库仑效率和/或更高的功率充电接受率(例如，更高的最大充电电流或充电电压)。这样，可以改善电池系统12的捕获、存储和/或分配效率。
39.为了便于控制电能的捕获和存储，电池系统12可以另外包括控制模块32。更具体地，控制模块32可以控制电池系统12中的部件的操作，例如能量存储部件14、交流发电机18和/或电动机22内的继电器(例如，开关)的操作。例如，控制模块32可以调节由每个电池模块28或30捕获/提供的电能量(例如，降低电池系统12的额定值和重新确定电池系统12的额定值)，在电池模块28和30之间执行负载平衡，确定每个电池模块28或30的充电状态，确定每个电池模块28或30的温度，确定任一电池模块28和30的预测温度轨迹，确定电池模块28或30的预测寿命，确定任一电池模块28或30的燃料经济性贡献，控制交流发电机18和/或电动机22输出的电压或电流的大小等。
40.因此，控制模块(例如，单元)32可以包括一个或多个处理器34和一个或多个存储器36。更具体地，一个或多个处理器34可以包括一个或多个专用集成电路(asic)、一个或多个现场可编程门阵列(fpga)、一个或多个通用处理器、或其任意组合。通常，处理器34可以执行与本文描述的过程有关的计算机可读指令。另外，处理器34可以是定点处理器或浮点处理器。
41.另外，一个或多个存储器36可以包括易失性存储器，例如随机存取存储器(ram)，和/或非易失性存储器，例如只读存储器(rom)、光盘驱动器、硬盘驱动器或固态驱动器。在一些实施例中，控制模块32可包括车辆控制单元(vcu)的一些部分和/或单独的电池控制模块。另外，如图所示，控制模块32可以与能量存储部件14分离，例如作为独立模块。在其他实施例中，电池管理系统36可以包括在能量存储部件14内。
42.在某些实施例中，控制模块32或处理器34可以从布置在能量存储部件14内和/或周围的各种传感器38接收数据。传感器38可以包括用于测量关于电池模块28或30的电流、电压、温度等的各种传感器。在从传感器38接收数据之后，处理器34可以将原始数据转换成电池模块28和30的参数估值。这样，处理器34可以将原始日期呈现为数据，该数据可以向车辆10的操作者提供与电池系统12的操作有关的有价值信息，并且可以在显示器21上显示与电池系统12的操作有关的信息。显示器21可以显示由设备10产生的各种图像，例如用于操作系统的gui或图像数据(包括静止图像和视频数据)。显示器21可以是任何合适类型的显
示器，例如液晶显示器(lcd)、等离子显示器或有机发光二极管(oled)显示器。另外，显示器21可以包括触敏元件，其可以为控制模块32的调整参数或者由处理器34处理的数据提供输入。
43.能量存储部件14可以具有与典型的铅酸电池的尺寸相当的尺寸，以避免为了容纳电池系统12而不得不对车辆10设计作出改变。例如，能量存储部件14可以具有与h6电池类似的尺寸，其尺寸可以约为13.9英寸x6.8英寸x7.5英寸。如图所示，能量存储部件14可以包括在单个连续壳体内。在其他实施例中，能量存储部件14可包括耦接在一起的多个壳体(例如，包括第一电池28的第一壳体和包括第二电池30的第二壳体)。在又一些其他实施例中，如上所述，能量存储部件14可包括位于车辆10的引擎盖下方的第一电池模块28，并且第二电池模块30可位于车辆10的内部。
44.更具体地，图3示出了电池系统12的部件的示意图。如上面在图2的讨论中所提到的，控制模块32可以调节由每个电池模块28或30捕获/提供的电能量(例如，降低电池系统12的额定值和重新确定电池系统12的额定值)，在电池模块28和30之间执行负载平衡，确定每个电池模块28或30的充电状态，确定每个电池模块28或30的温度，确定电池模块28或30的预测温度轨迹，确定电池模块28或30的预测寿命，确定电池模块28或30的燃料经济性贡献，控制交流发电机18和/或电动机22输出的电压或电流的大小等。特别地，控制模块32可以基于由传感器38测量的电池参数来测量充电状态(soc)和/或健康状态(soh)。
45.在一些实施例中，能量存储部件14可包括单个锂离子电池或串联耦接的多个锂离子电池。另外，预期其他可再充电电池化学组成。能量存储部件14可以将存储的能量释放到负载40，负载40可以包括点火系统16、车辆控制台20、显示器21、电动机22以及车辆10的任何其他电气部件。当能量存储部件14将存储的能量释放到负载40时，交流发电机18和/或电动机22可以向能量存储部件14提供能量以补充先前释放到负载40的存储的能量。传感器38可以测量能量存储部件14的电池参数，并且传感器38可以将测量值传输到控制模块32。能量存储部件14的电池参数可以包括端子电压测量值、端子电流测量值和电池温度测量值。控制模块32处理测量的电池参数，如下面详细描述的，以估计能量存储部件14的soc、与能量存储部件14的等效电路模型相关联的两个电阻、以及与等效电路模型相关联的电容。如下面关于图4所讨论的，等效电路模型可以是1
‑
rc等效电路模型。
46.此外，可以理解的是，本文描述的系统和方法可以用于改变能量存储部件14的化学组成。例如，能量存储部件14的soc、电阻和电容可以代表单电池单元的锂离子电池或多电池单元的锂离子电池、单电池单元的铅酸电池或多电池单元的铅酸电池、其某些组合(例如，与铅酸电池并联电耦接的锂离子电池)，或任何其他单电池单元电池化学组成或多电池单元电池化学组成。此外，在具有并联电耦接的多个电池化学组成的能量存储部件14中，soc、电阻和电容可以代表整个能量存储部件14，或者可以为多种电池化学组成中的每一种电池化学组成计算soc、电阻和电容。
47.图4示出了能量存储部件14的1
‑
rc等效电路模型42。1
‑
rc等效电路模型42将电池参数(例如，开路电压(ocv)44、电阻46和48以及电容50)与由传感器38测量的测得参数(例如，端子电压52、端子电流和电池温度)相关联。另外，1
‑
rc等效电路模型42提供了在车辆10运行期间实时估计ocv的机制。使用其他方法来测量ocv，例如库仑计数方法，当电池系统12处于休眠状态达到很长一段时间时，ocv可以是可测量的。也就是说，在电池系统12已经休
眠了一个或多个小时之后，电池系统12的ocv可以是可测量的。因此，通过使用1
‑
rc等效电路模型42，不再使用休眠时段，并且可以在电池系统12在负载40下运行的同时估计ocv。
48.在1
‑
rc等效电路模型42中，电阻46(即，r0)表示能量存储部件14的电流路径的欧姆电阻，电阻48(例如，r1)表示能量存储部件14的电荷转移电阻，电容50(例如c1)表示能量存储部件14的双层电容。1
‑
rc等效电路模型42被视为归属于单个电阻器
‑
电容器对(例如，电阻48和电容50)的1
‑
rc等效电路模型。使用1
‑
rc等效电路模型42使得能够在车辆10的实时驾驶状态期间确定ocv 44。在1
‑
rc等效电路模型42中，电阻46和48以及电容50通常可以是能量存储部件14的非时变参数。或者，可以用于确定能量存储部件14的充电状态的ocv 44通常可以是能量存储部件14的时变参数。即，当能量存储部件14在一段时间充电和放电时，ocv 44将随时间增加和减少。
49.ocv 44、电阻46和48以及电容50的准确估计可以有益于控制能量存储部件14以获得更长的电池寿命和提高混合动力电动车辆的燃料效率。例如，图5示出了图表54，其提供了能量存储部件14的ocv 44与充电状态(soc)之间的关系。soc沿着图表54的横坐标56以百分比显示。另外，ocv 44沿图表54的纵坐标58以电压显示。曲线60表示能量存储部件14的ocv 44与soc之间的关系。例如，曲线60可用作查找表以提供能量存储部件14的准确soc表示。当基于测量的电池参数确定1
‑
rc等效电路模型42的ocv 44时，ocv 44的值可以与相应的soc百分比匹配。soc百分比可以在车辆10的运行期间向车辆10的操作者实时地提供能量存储部件14的剩余电池寿命的准确指示。此外，可以理解，ocv 44随着soc的变化而变化。例如，当能量存储部件14的soc增加或减少时，ocv 44不会稳定在某一电压下。
50.返回图4的讨论。1
‑
rc等效电路模型42可以最初以离散时间导出，以将ocv 44、电阻46和48以及电容50的估计值(即，电池参数估计值)与能量存储部件14的测得的端子电压52和测得的电流相关联。卡尔曼滤波器用于这种关联，以根据测得的端子电压52和测得的电流确定电池参数估计值。使用卡尔曼滤波器，控制模块32可以在有限地依赖于预定义的电池参数的情况下实时更新电池参数估计值。
51.对于任何任意电流源i，可以通过杜哈梅尔叠加定理计算能量存储部件14的电压52(例如，v)：
[0052][0053]
其中ξ是积分的虚变量。等式1右侧的前两个项(即，v
oc
和ir0)给出能量存储部件14的欧姆描述，因为电压52与ocv44减去欧姆降ir0有关。此外，等式1右侧的第三项对应于叠加积分，通过该叠加积分，过往电流对ocv44的影响超过对能量存储部件14予以表征的平均soc的变化的一阶效果。由于指数加权函数，较旧的电流
‑
电位数据点的影响指数级别地小于最近的数据点。
[0054]
可以在两个任意时间步t
k
‑1和t
k
评估等式1，以得到：
[0055][0056]
和
[0057][0058]
可以假设以固定的时间间隔测量电池电流和电压，例如：
[0059][0060]
因此，可以组合等式2和3以得到：
[0061][0062]
如果i(ξ)近似为l
k
‑1的阶跃电流，则可以减去等式5以得到：
[0063]
v
k
＝av
k
‑1 (1
‑
a)v
oc
‑
i
k
r0‑
i
k
‑1[(1
‑
a)r1‑
ar0]
ꢀꢀꢀ
(6)
[0064]
其中另外，如果通过(i
k
‑1 i
k
)/2的阶跃电流对i(ξ)求近似，则可以减去等式5以得到：
[0065][0066]
此外，如果通过逐段线性方程(i
k
‑1 (i
k
‑
i
k
‑1))/((t
k
‑
t
k
‑1)x(ξ
‑
t
k
‑
t
))对i(ξ)求近似，可以减去等式5以得到：
[0067][0068]
由于i(ξ)对于所有时间间隔的连续性质，等式8可以比等式6和7更精确地计算电池电压，等式6和7是基于不连续的阶跃电流导出的。为方便起见，等式8可以改写为：
[0069]
v
k
＝θ1‑
i
k
θ3‑
i
k
‑1θ2 θ4v
k
‑1ꢀꢀꢀ
(9)
[0070]
其中
[0071]
θ1＝(1
‑
θ4)v
oc
ꢀꢀꢀ
(10)；
[0072][0073]
和
[0074][0075]
在本实施例中，等式9用作电池模型，以使用卡尔曼滤波方法从测得的电压和电流数据同时递归地估计四个参数，θ1‑
θ4。物理参数，例如v
oc
(例如，ocv 44)、r0(例如，电阻46)、r1(例如，电阻48)、c1(例如，电容50)和时间常数τ＝r1c1是根据以下等式从θ1‑
θ4中提取的：
[0076]
(16)
[0077][0078][0079]
和
[0080][0081]
可以使用双模型过程来实现标准卡尔曼滤波方法。为了估计这四个参数θ1‑
θ4，使用标准卡尔曼滤波方法，状态转移模型可以描述为：
[0082]
θ(k 1)＝θ(k) r(k)
ꢀꢀꢀ
(18)
[0083]
其中r(k)是过程噪声向量，θ(k)是状态向量，其中θ1‑
θ4是状态向量的四个参数。与等式18和其他卡尔曼滤波方法相反，下面描述的备用卡尔曼滤波方法没有明确地计算状态转移。相反，可以使用以下等式描述状态测量模式：
[0084]
v(k)＝φ(k)
′
θ(k) w(k)
ꢀꢀꢀ
(19)
[0085]
其中w(k)是测量噪声，φ(k)是如下等式描述的回归向量：
[0086]
φ(k)＝[1
‑
i(k
‑
1)
‑
i(k)v(k
‑
1)]
′ꢀꢀꢀ
(20)
[0087]
此外，通过使用替代卡尔曼滤波方法，使用以下等式计算用于soc估计的卡尔曼增益k(k)：
[0088][0089]
因此，状态转换模型可以更新如下：
[0090]
θ(k)＝θ(k
‑
1) k(k)[y
‑
φ
′
θ(k
‑
1)]
ꢀꢀꢀ
(22)
[0091]
其中y是在时间步k测得的电池电压52。此外，协方差矩阵p(k)通过以下等式更新：
[0092]
p(k)＝[i
‑
k(k)φ(k)
′
]p(k
‑
1) r
ꢀꢀꢀ
(23)
[0093]
在替代卡尔曼滤波方法中，θ1被识别为快时变参数，θ2‑
θ4被识别为慢时变参数。为了使用替代卡尔曼滤波方法实现此特征，使用以下等式定义过程噪声矩阵r：
[0094][0095]
其中r是小值(例如，r可以设置为大约0.001)，用于使用随机游走概念控制θ1随时间的可变性。这使卡尔曼滤波器能够在不使用状态转移模型的情况下估计可再充电电池的时变参数，其通常与通过电流积分或库仑计数计算soc变化有关。另外，通过将过程噪声矩阵r(k)的对角线上的三个其他值建立为零，不变参数或慢时变参数θ2‑
θ4通常保持不变。
[0096]
协方差矩阵p(k)的计算扮演的角色是提高估计电池参数的数值效率和准确度。由于在计算机中的截断和舍入误差，替代卡尔曼滤波方法可能导致p(k)的对称性和正定性丧失，或者可能导致发散。因此，在替代卡尔曼滤波方法中，实现用于计算p(k)的u
‑
d因数分解方法。u
‑
d因数分解方法改善了p(k)的正定性和对称性，这可以导致获得高估计精度和鲁棒性。
[0097]
使用u
‑
d因数分解方法实现卡尔曼滤波器函数可以包括计算以下等式：
[0098]
f＝u
′
φ，g＝df，α0＝1
ꢀꢀꢀ
(25)
[0099]
在计算等式25时，对于从j＝1到j＝m的值的范围，计算以下等式：
[0100][0101]
另外，对于从j＝1到j＝i
‑
1的值的范围，计算以下等式：
[0102][0103]
此外，r的值通过以下等式加到d(1,1)：
[0104]
d(1，1)＝d(1，1) r
ꢀꢀꢀ
(28)
[0105]
增益也使用以下等式计算出：
[0106]
k＝b/α
ꢀꢀꢀ
(29)
[0107]
在完成等式25
‑
29的计算后，使用以下等式计算新估计值：
[0108]
θ＝θ k(y
‑
φ
′
θ)
ꢀꢀꢀ
(30)
[0109]
现在转向图6a和6b，方法70的流程图示出了使用上述替代卡尔曼滤波方法确定能量存储部件14的电池参数的方法。最初，在方框72，由控制单元32执行初始化。在初始化期间，可以将参数θ1‑
θ4的任意初始值分配给状态向量θ。参数的初始值可以是任意的，因为替
代卡尔曼滤波方法使用起始点来确定参数的最终值，并且替代卡尔曼滤波方法不需要与参数θ1‑
θ4的实际值接近的起始点。因此，可以在初始化期间将状态向量初始化为[0 0 0 0]’。另外，在初始化期间，u，可以被初始化为单位对角矩阵(即，具有分配给对角元素的1值的矩阵)，可以将d初始化为具有分配给每个对角元素的大值(例如，1000)的对角矩阵。当来自传感器38的能量存储部件14的新电压和电流测量值对控制模块32变得可用时，该初始化可以使方法70开始训练测量模型，其由上面的等式9表示。此外，初始化可以包括设置相对误差容限(rtol)和比例因子的值。rtol可以表示平滑度和平坦度标准，比例因子可以是电池电流的因子，以避免在定点微处理器上可能遇到的溢出数学错误。
[0110]
随后，在方框74，当来自传感器38新数据在控制模块32处变得可用时，可以更新数据计数器。例如，以下等式可以表示数据计数器：
[0111]
k＝k 1
ꢀꢀꢀ
(31)
[0112]
其中k是当前数据计数。此外，在方框76处，将在当前时间步和先前时间步测得的数据值分配给回归向量φ(k)。测得的数据值可以包括测得的电池电压和电流值。
[0113]
在方框78处，执行卡尔曼滤波器函数。卡尔曼滤波器函数可以包括上述等式25
‑
30。在方框80处，通过执行卡尔曼滤波器函数，可以从θ1‑
θ4的值中提取电池参数(例如，ocv 44、电阻46和48、电容50和时间常数τ)的值。此外，在方框82处，使用从θ1‑
θ4的值中提取的ocv 44的值，可以基于存储在控制模块32的存储器36中的ocv至soc查找表来确定能量存储部件14的soc。ocv至soc查找表通常可以基于类似于图5中所示的曲线60的曲线。
[0114]
方框84
‑
96涉及监测电阻46和48的估计值的收敛。通过监测电阻46和48的收敛，控制模块32可以确认电阻46和48的估计值是非时变的。因此，在方框84处，递归地确定电阻46和48的样本数量l的移动样本窗口的平均值和方差。具体地，以下等式可用于确定电阻46和48的平均值和方差：
[0115]
如果k<＝l，那么
[0116][0117]
否则
[0118][0119]
其中μ
k
和σ
k2
是在时间步k评估的电阻r
j
的样本均值和方差。如等式32和33所示，为计算样本均值和方差而开发的递归公式不涉及每次数据复制，只有新数据和最旧数据涉及平均值和方差计算。因此，控制模块32有效地计算等式32和33。
[0120]
随后，在方框86，控制模块32计算每个电阻46和48的方差和平均值平方的比值，并且控制模块32确定两个电阻46和48中的哪一个具有更大的比值。此外，一旦确定了哪个比值更大，控制模块32将更大的比值与相对误差容限(rtol)的平方进行比较。该比较用作平滑度和平坦度测试。平滑度和平坦度测试可用于在系统的学习时间段之后确定所估计的参数是否不再随时间变化。因此，为了判断电阻46和48的收敛，平滑度和平坦度测试需要连续n次地合格。n的值可以是5个循环、10个循环、15个循环、或者可以可靠地指示电阻46和48的收敛的任何其他数量的循环。此外，如果对于方法70使用1秒的采样频率，则方框86的平滑度和平坦度测试的每个循环可以运行0.1纳秒的总时间段。
[0121]
随后，在方框88，如果较大的比值小于rtol的平方，则通过将1加到前一循环计数器值来更新循环计数器(例如，icheck)。相反，如果较大的比值大于rtol的平方，则在方框90，将循环计数器重置为零，并且方法70返回方框72处的初始化。如果在方框88处更新循环计数器(即，如果较大的比值小于rtol的平方)，则在方框92处确定循环计数器是否超过n的值。如果循环计数器未超过n的值，则该方法返回方框72处的初始化。
[0122]
然而，如果循环计数器已经超过n的值，则在方框94，将收敛电阻46和48的平均值和温度值存储到控制模块32的存储器36中。通过根据温度测量电阻46和48的值，控制模块32可以观察到电阻46和48如何随时间降低。测量电阻46和48的降低可以使得能够从能量存储部件14的欧姆电阻增长的角度表征能量存储部件14的健康状态(soh)。一旦存储了收敛电阻46和48的平均值以及温度值，则在方框96处将循环计数器重置为零，并且方法70返回方框72处的初始化。可以重复方法70，直到控制模块32提供停止操作的指示。
[0123]
现在转到图7a和7b，具有样本数量l的电阻46、48的样本平均值和方差的递归计算结果(例如，如图6a的方框84所实现的)可以被存储在具有l个存储位置102的环形缓冲器100中。例如，对于10的样本数量l，环形缓冲器100可以包括10个存储位置102，如图所示。当计算包括电阻46和48的样本平均值和方差的新数据104时，新数据104可以存储在存储位置102之一中。在图7a中，空的存储位置102按照记录每个样本的时间顺序用新数据104填充。或者，如图7b所示，新数据104以如下方式存储在环形缓冲器100中最旧数据值的存储位置
102中(例如，r
j
(k
‑
9))：使得新数据104存储在环形缓冲器100内并且从环形缓冲器100中移除旧数据106(例如，r
j
(k
‑
9))。因此，只有l个存储位置102可用于存储样本数量l的样本平均值和方差，并且旧数据106从环形缓冲区100移除。
[0124]
在实时估计能量存储部件14的容量方面也可能有益处。能量存储部件14的容量可以被称为能量存储部件14的健康状态(soh)。能量存储部件14的soh可以指示能量存储部件的额定容量的变化。下面详细讨论的是用于估计能量存储部件14的容量的两种互补方法。结合图8a和图8b讨论的第一种方法使用卡尔曼滤波器和等效电路电池模型提供实时电池电流和电压的线性回归。关于图9a和图9b讨论的第二种方法涉及监测两个开路电压弛豫事件，并从两个开路电压弛豫事件计算soh。控制模块32可以执行两种方法中的每一种以估计能量存储部件的容量。
[0125]
现在转向图8a和8b，其示出了用于计算能量存储部件14的soh的方法120。方法120基于在每个时间步串行运行的两个线性回归模型。因为方法120依赖于线性回归模型，所以方法120可以提供用于估计soh的数值稳定性，因为即使在任意选择参数的初始值时线性回归模型的参数也将收敛。任意选择参数的初始值可以提供优于扩展卡尔曼滤波器(ekf)模型的独有优势，扩展卡尔曼滤波器(ekf)模型由于ekf模型的非线性通常依赖于对可再充电电池的参数的准确初始猜测。另外，方法120可以提供比ekf模型更高程度的对测量噪声和数据缺陷的容限。
[0126]
可以使用能量存储部件14的部分放电来使用以下等式确定能量存储部件14的soh(即，容量)：
[0127][0128]
其中q是以安培小时(ah)为单位的能量存储部件14的容量，i是以安培为单位的电流，soc是当前soc值，并且soc0是初始soc值。如上所述，当可获得开路电压的估计值或已知值时，可以从类似于图5的曲线60的查找表获得soc。等式34也可以重写为以下等式：
[0129]
soc＝soc0 w∫
0t
idt
ꢀꢀꢀ
(35)其中w等于100/q且是能量存储部件14的累积ah吞吐量。等式35可以用作方法120中的电池容量估计的主要等式。soc0和w都是非时变参数，如果soc值是已知的，它们是可以估计出的。如上面在图6a和6b讨论中所讨论的那样，可以使用卡尔曼滤波技术估计soc值。因此，在方法120中，使用两个线性回归模型。也就是说，方框122至132可以确定能量存储部件14的soc值，方框134
‑
156可以使用作为输入的soc和ah吞吐量以及收敛确定来提供soc0和w的估计。因此，可以根据w的估计值容易地计算电池容量。
[0130]
在方框122处，方法120被初始化。在初始化期间，可以为状态向量θ将任意初始值分配给θ1‑
θ4。参数的初始值可以是任意的，因为替代卡尔曼滤波方法使用起始点来确定参数的最终值，但是替代卡尔曼滤波方法不依赖于接近参数θ1‑
θ4的实际值的起始点。因此，状态向量θ可以在初始化事件期间被初始化为[0 0 0 0]’。另外，在初始化期间，u，可以被初始化为单位对角矩阵(即，具有分配给对角元素的1值的矩阵)，d可以被初始化为具有分配给每个对角元素的大值(例如，1000)的对角矩阵。当来自传感器38的能量存储部件14的新
电压和电流测量值对于控制模块32变得可用时，该初始化可以使方法70能够开始训练测量模型，其由上面的等式9表示。此外，初始化可以包括设置相对误差容限(rtol)和比例因子的值。rtol可以表示平滑度和平坦度标准，比例因子可以是电池电流的因子，以避免在定点微处理器上可能遇到的溢出数学错误。可以对用于soc估计的线性回归模型和用于soc0和w的线性回归模型这两者执行方框122处的初始化。
[0131]
在方框124处，当来自传感器38的新数据在控制模块32处变得可用时，可以更新数据计数器。例如，等式31可以用作当前数据计数的表示。此外，当数据可用时，还可以使用以下等式更新安培小时吞吐量：
[0132]
q＝q δq
k
ꢀꢀꢀ
(36)
[0133]
其中q是安培
‑
小时吞吐量。此外，在方框126处，将当前时间步和先前时间步测得的数据值分配给回归向量φ(k)。测得的数据值可以包括测得的电池电压和电流值。
[0134]
在方框128处，执行卡尔曼滤波器函数。卡尔曼滤波器函数可以包括上述等式25
‑
30。通过执行卡尔曼滤波器函数，除了电阻46和48、电容50和时间常数τ之外，还可以在方框130处从θ1‑
θ4的值中提取ocv 44的值。另外，在方框132处，使用从θ1‑
θ4的值中提取的ocv 44的值，可以基于存储在控制模块32的存储器36中的ocv至soc查找表来确定能量存储部件14的soc。ocv至soc查找表通常可以基于类似于图5中所示的曲线60的曲线。
[0135]
随后，在方框134处，控制模块32确定估计的soc是否可用作输入以继续进行能量存储部件14的后续容量估计。例如，控制模块32可以基于是否满足电池端子温度要求(例如，能量存储部件的端子的温度小于约25摄氏度)以及是否已经超过最小soc学习时间段(例如，电流测量和初始条件之间的时间间隔大于100秒)进行确定。如果不满足这些条件中的任何一个，则在方框136处，将电流和电压测量值重置为先前值，并且方法120在方框124的数据计数器更新处重新开始。
[0136]
如果soc可用于容量估计，则在方框138处，跟踪最小和最大soc值以用于计算最大soc摆动。此外，在方框140处将安培
‑
小时吞吐量分配给第二回归向量。第二回归向量由φi表示并且表示以下等式：
[0137]
φ1(k)＝[1q/scale 1]
′ꢀꢀꢀ
(37)
[0138]
其中q是安培小时吞吐量而标度l是在定点微处理器上计算时为了避免溢出数学错误的估计容量的比例因子。
[0139]
一旦将安培
‑
小时吞吐量分配给第二回归向量，在方框142执行卡尔曼滤波器函数。执行卡尔曼滤波器函数可以更新协方差矩阵u1和d1以及用于容量估计的参数向量θ1，其与用于soc估计的参数向量θ一起被初始化为[0 0 0 0]'。此外，用于容量估计的卡尔曼增益k(k)由以下等式表示：
[0140][0141]
并且协方差矩阵p1(k)通过以下等式更新：
[0142]
p1(k)＝[i
‑
k1(k)φ1(k)
′
]p1(k
‑
1)
ꢀꢀꢀ
(39).
[0143]
当与soc估计的等式23比较时，等式39不包括过程噪声矩阵r，因为soc0和w都是非时变参数。
[0144]
随后，在方框144，可以从参数向量θ1中提取估计容量值q。例如，估计容量值q可以由以下等式表示：
[0145]
q＝
‑
100/θ12ꢀꢀꢀ
(40)
[0146]
其中θ12是来自参数向量θ1的第二参数。在提取估计容量值q之后，可以使用方框146
‑
156来监视容量估计的收敛。
[0147]
因此，在方框146处，递归地确定估计容量值q的具有样本数量l的移动样本窗口的平均值和方差，并将其存储在类似于上述环形缓冲器100的环形缓冲器中。在确定平均值和方差之后，在方框148处，将方差与均方值平方的比值与rtol的平方进行比较。该比较用作平滑度和平坦度测试。平滑度和平坦度测试可用于确定估计参数(例如，估计容量值q)在系统的学习时间段之后是否不再随时间改变。因此，为了判断估计容量值q的收敛，平滑度和平坦度测试要连续n次地合格。n的值可以是5个循环、10个循环、15个循环、或者可以可靠地指示估计容量值q的收敛的任何其他循环数量。此外，如果对于方法120使用1秒的采样频率，则方框148的平滑度和平坦度测试的每个循环可以运行0.1纳秒的总时间段。
[0148]
随后，在方框150，如果该比值小于rtol的平方，则通过将1加到前一循环计数器值来更新循环计数器(例如，icheck)。相反，如果该比值大于rtol的平方，则在方框136处，将电流和电压测量值重置为先前值，并且方法120在方框124的数据计数器更新处重新开始。如果循环计数器在方框150处被更新(即，较大的比值小于rtol的平方)，在方框152处确定循环计数器是否已超过n的值以及soc摆动(最大soc减去最小soc)是否大于soc0的20％。在方框152处用于soc摆动的阈值也可以是soc0的另一百分比，这取决于具体的能量存储部件14。例如，百分比可以是25％或高达50％，或者百分比可以是soc0的10％或15％。如果循环计数器没有超过n的值和/或soc摆动不够大，则方法120返回到方框136以重置电压和电流值。
[0149]
如果循环计数器已经超过n的值并且soc摆动足够，则在方框154，将提取的估计容量值q保存到存储器36。在估计容量值q被保存到存储器36之后，循环计数器被重置为零，并且方法120返回到方框136。另外，方法120可以如上所述重复，直到控制模块32接收到停止估计能量存储部件14的估计容量值q的指示。
[0150]
图9a和9b示出了用于通过监测两个有效开路电压(ocv)弛豫事件并且在两个有效ocv弛豫事件之间随时间对电流求积分来计算能量存储部件14的soh的方法160。可以在较长的休眠时段之后(例如，在车辆熄火超过2小时之后)启动车辆10时立即测量两个有效ocv弛豫事件中的第一个，可以在能量存储部件14已经弛豫一段较长的时间段时(例如，在车辆10熄火并且已经过预定量的时间之后)测量两个有效ocv弛豫事件中的第二个。另外，两个有效ocv弛豫事件之间的时间可以被限制为限制当前偏移误差的累积并且提高soh测量的准确度。
[0151]
在方框162处，方法160被初始化。在初始化期间，设置开始安培小时吞吐量∑q，的时间计数器t1和t2以及二进制ocv有效性状态iocv。如果在控制模块32醒来的时刻(例如，当车辆10起动时)能量存储部件14的ocv完全弛豫，则在初始化期间可以将iocv的值设置为1。否则，在初始化期间可以将iocv的值设置为0。在室温下停放车辆一段较长的时间段(例如，大于一小时)之后，ocv可以完全弛豫，其中具有小于弛豫电流阈值l
relax
的可忽略的电流汲取值。确认ocv是否完全弛豫的一种方法是监测ocv随时间的变化率，因为ocv可能随时间渐
近地变化。在室温下较长时间段之后，ocv可以随时间非常缓慢地变化，并且对于方法160可以使用这种ocv以足够的精度表示弛豫的ocv。
[0152]
在初始化之后，在方框164处，当来自传感器38的新数据在控制模块32处变得可用时，计数器被更新。例如，等式31可以用作当前数据计数的表示。一旦更新了计数器，在方框166处，从能量存储部件14的传感器38读取数据。数据可包括电池电流、电压和安培
‑
小时吞吐量的阶跃变化。因此，随着安培
‑
小时吞吐量的阶跃变化，在方框168，使用以下等式更新累积的安培
‑
小时吞吐量：
[0153]
∑q＝∑q δq
k
ꢀꢀꢀ
(41)
[0154]
其中∑q表示累积安培小时吞吐量，而δq
k
表示安培小时吞吐量的阶跃变化。
[0155]
在方框170处，控制模块32确定电池电流是否满足预定的弛豫电流阈值l
relax
。可以通过在当前时间步测量电池电流的绝对值并将该绝对值与i
relax
进行比较来完成该确定。如果电池电流的绝对值大于i
relax
，则电池没有充分弛豫。此时，可以监视时间步，并且在方框172处iocv值可以保持为0或者设置为0。因此，方法160可以在方框164的计数器更新处再次开始。
[0156]
或者，如果电池电流的绝对值小于i
relax
，则可以在方框174处确定iocv有效性状态。即，如果iocv的值是1，则能量存储部件14已经处于弛豫状态足够的时间量，并且在方框178处可以将ocv设置为当前电压读数v
k
。或者，如果iocv是0，则在方框176，确定是否已过去足够时间量以至于使能量存储部件14充分弛豫。例如，可以将从初始化到当前时间步的时间量与弛豫阈值时间t
relax
进行比较。如上所述，在标准操作条件下，t
relax
可设定为1小时、2小时或更长时间。然而，t
relax
可以基于诸如电池温度的外部条件而增加或减少。通常，电池的温度越高，t
relax
时间可能越短。如果尚未满足弛豫阈值时间t
relax
，则方法160返回到方框164的计数器更新以再次开始方法160。
[0157]
或者，如果已经满足弛豫阈值时间t
relax
，能量存储部件14已经处于弛豫状态足够的时间量，则可以在方框178处将ocv设置为当前电压读数v
k
。根据由当前电压读数v
k
确定的ocv值，可以从类似于图5的曲线60的ocv对soc查找表确定当前soc。
[0158]
随后，在方框182处，确定计数器值k是否大于样本数量l。如果计数器值k大于l，则通过在方框184将该值加1来更新模数变量m的值。模数变量m表示已满足ocv标准的累积次数。或者，如果计数器值k不大于样本数量l，则在方框186处将l的值设置为计数器值k加1。在更新模数变量m和/或样本数量l的值之后，在方框188执行模数运算，以在控制模块32的存储器36内的两个存储器存储位置之间交替记录当前soc值。因此，方法160将交替地在方框190将当前soc值存储为soc1，并在方框192将当前soc值存储为soc2。此外，时间和累积安培
‑
小时吞吐量也可以在方框190和192处存储在存储器中。
[0159]
在方框194，确定在为soc1存储的值和为soc2存储的值之间是否存在足够的soc摆动。足够的soc摆动可以被描述为soc
min
，并且soc
min
可以表示两个存储的soc值之间的阈值百分比差异。例如，soc
min
可以是5％摆动、10％摆动、15％摆动或为了能量存储部件14的精确容量估计而在两个soc值之间建立足够的差异的另一百分比的摆动。如果摆动不够，则方法160可以返回到方框164处的计数器更新。
[0160]
如果soc摆动足够，则在方框196，控制模块32可以确定总电流积分时间是否小于最大允许时间。例如，如果总电流积分时间大于最大允许时间，则电流偏移误差可能以不期
望的方式影响容量估计。因此，可能希望将总电流积分时间限制为小于约50小时。如果已经超过最大允许时间，则方法160可以返回到方框164的计数器更新。
[0161]
或者，如果尚未超过最大允许时间，则可以在方框198计算能量存储部件14的估计电池容量q。为了计算估计电池容量q，可以使用以下等式：
[0162][0163]
其中，用于计算的值是从方框190和192中存储的值获得的。随后，在方框200，将估计电池容量q存储到存储器36中。此外，在存储估计电池容量q时，方法160可以返回到方框164的计数器更新，并且方法160可以递归地操作，直到向控制模块32提供停止方法160的指示。
[0164]
上述两种容量估计方法120和160在实践中彼此互补。方法120可以在两个ocv弛豫事件之间不使用soc摆动，而可以在车辆10运行时使用实时soc摆动。因此，由于提前启停和混合动力电动车辆的设计原理，方法120可以特别适用于提前启停和混合动力电动车辆，以在存在动能的过剩供应时最大化收获能量的充电率并且在存在峰值功耗需求时最大化放电率以改善燃料经济性。因此，能量存储部件14可以在大约50％的soc处建立待机(休眠)模式。因此，能量存储部件14通常在两次ocv测量之间具有小的间隔。或者，由于方法160的简单性、准确性和低实施成本，当车辆10在典型操作期间经历若干休眠时段时，方法160可适用以增强能量存储部件14的容量估计值的鲁棒性和准确性。
[0165]
偶尔验证能量存储部件14的容量估计结果可能是有益的。因此，图10是描述能量存储部件14的估计电池容量q的验证过程的流程图210。在方框212处，控制模块32确定能量存储部件14的容量估计是否有效。例如，如果在一段较长的时间量(例如，自上次有效容量估计结果以来超过两个月)尚未有来自方法120和160的任何有效容量估计结果，如果最近计算的容量估计值超过变化阈值(例如，最近计算的容量估计值比之前计算的容量估计值大5％、大10％、大15％、或者大更多)，或者如果有效容量估计值的准确度不在可接受的范围内，则容量估计值可能是无效的。当满足这些标准中的任何一个时，控制模块32可以初始化流程图210的验证过程。
[0166]
在确定有效容量估计值的准确度是否不在可接受的范围内时，可以使用soc/ocv测量值的误差以及经由呈现当前测量精度的电流积分来计算容量估计值的误差。例如，卡尔曼滤波方法可以包括误差估计值或最大准确度估计值(例如，大约3％)，并且如本文所述的收敛标准的评估可以对这种误差估计进行精炼。在对误差估计进行精炼之后，基于在电池系统12的观察条件下的预期老化的与时间或能量吞吐量相关的误差增加可以应用于误差估计。因此，如果在一定精度内确定容量估计值，则误差稳定地增加，直到有机会以较低误差达到新估计值为止，这将重置误差估计值。
[0167]
如果确定容量估计值有效，则可以在方框214更新当前容量估计值。在方框214处更新的容量可以在方框216处用作电流积分过程期间的能量存储组件14的实际容量，该电流积分过程使用电流容量估计(例如，使用上面的等式35)的。使用能量存储部件的实际容量的电流积分过程的结果是在方框218处计算出的在一积分时间段上的soc(1)的值。
[0168]
或者，如果确定容量估计值无效(例如，平均误差太高，并且重置不会提高准确
度)，则可以在方框220处，在与方框216的电流积分过程并行的电流积分过程中使用候选容量。候选容量可以比方框216的电流积分过程中使用的实际容量小5％。使用候选容量的电流积分过程的结果是在方框222处计算出的在一积分时间段上的soc(2)的值。此外，在方框216和220处执行电流积分过程时，可以使用能量存储部件14的测量参数224。测量参数224可以包括系统/传感器规格、能量存储部件14的温度、电流和电压。
[0169]
为了有效地使用流程图210的验证过程，可以使用并行电流积分过程开始时的准确初始soc。因此，测量参数224可以为电池控制模块32提供值以计算准确的初始soc。例如，如果车辆10在方框216和220处的并行电流积分过程之前退出能量存储部件14的较长的休眠时段，则可以从测得的参数224测量能量存储部件14的开路电压44，如上面详细讨论的那样。使用开路电压44，存储在存储器36中的开路电压
‑
充电状态查找表可以由控制模块32查阅，以计算准确的初始soc值。还可以理解，在并行电流积分过程开始时的初始soc与评估变化容量估计值的时间点处的最终soc之间的显著soc增量(例如，soc摆动)对准确评估容量估计值的有效性可能是有益的。该值可由δsoc
δq
表示。
[0170]
另外，为了与并行电流积分过程产生的soc值进行比较，可以在评估容量估计值有变化的时间点在方框226处使用卡尔曼滤波器。例如，使用上面讨论的与6a和6b有关的方法70，能量存储部件14的soc估计值(soc(3))可以在方框228由控制模块32计算。可选地或另外地，如果能量存储部件14处于休眠状态，则可以在方框230从能量存储部件14获得ocv 44。在方框232由控制模块32对ocv至soc查找表的咨询导致方框234处对soc(4)的计算。
[0171]
在计算出soc(1)
‑
soc(3)和/或soc(4)之后，可以在方框236处由控制模块32执行soc值的方向比较。例如，可以将在评估容量估计值有变化的时间点处的soc(3)和soc(4)的值与在相同时间点的soc(1)和soc(2)的值进行比较。该比较可以提供关于能量存储部件14的新容量估计值的许多细节。如果soc(3)和/或soc(4)的值落入soc(1)和soc(2)的值的范围内或者soc(3)和/或soc(4)的值相对于soc(1)处于与soc(2)的值相同的方向上，则控制模块32可以将从soc(3)和/或soc(4)计算的估计容量值视为有效的。在这种情况下，控制模块32可以计算此时间点的容量，并在方框214处将实际容量值更新为新计算出的估计容量。
[0172]
或者，如果soc(3)和/或soc(4)的值不落在soc(1)和soc(2)的值的范围内，则控制模块32可以在方框236认为来自soc(3)和/或soc(4)的估计容量值无效。在这种情况下，流程图210的过程可以在不同的初始soc值处重新开始，或者可以在稍后的时间计算soc(3)和/或soc(4)以确定基于soc(3)和/或soc(4)的容量值是否有效。
[0173]
转到图11中，图表240示出了图10的流程图210中详述的过程。纵坐标242以百分比表示能量存储部件14的soc。横坐标244表示测量soc的时间。线246表示使用能量存储部件的实际容量时能量存储部件14的电流积分，线248表示使用能量存储部件14的候选容量时能量存储部件14的并行电流积分，其如图所示，比能量存储部件14的实际容量低5％。
[0174]
如图所示，并行电流积分过程在时间t1以准确的初始soc值250开始。可以通过开路电压(ocv)测量值以及ocv测量值与ocv至soc查找表的比较来获得准确的初始soc值250。另外，可以执行并行电流积分过程直到时间t2。时间t2可以表示并行电流积分过程的δsoc
q
(例如，充电状态摆动)足以准确验证能量存储部件14的估计容量的时间。例如，充电状态摆动可能是用来准确验证估计容量的大约20％的摆动。在其他情况下，充电状态摆动可以是用来准确验证估计容量的大约10％、15％或高达25％或更大的摆动。
[0175]
随着并行积分过程接近时间t2，绘制了例如使用卡尔曼滤波方法计算的能量存储部件的若干soc测量值252和254。深色soc测量值252表示对于容量估计更新无效的soc测量值。例如，与并行电流积分线246和248的soc摆动(例如，δsoc
δq
)相比，线246和soc测量值252之间的soc变化(例如，δsoc
k
)处于相反的方向。因此，如果由控制模块32计算深色soc测量值252的值，则控制模块32将等待随后的有效soc测量值以更新容量估计值。
[0176]
或者，浅色soc测量值254表示对于容量估计值更新有效的soc测量值。例如，线246和soc测量值254之间的soc变化(例如，δsoc
k
)与并行电流积分线246和248的soc摆动(例如，δsoc
δq
)方向相同。因此，如果由控制模块32计算浅色soc测量值254的值，则控制模块32将在计算浅色soc测量值254时更新容量估计值。
[0177]
所公开的实施例中的一个或多个可以单独或组合地提供一种或多种技术效果，包括确定电池时变变量和非时变变量、确定电池的充电状态、确定电池的健康状态、以及验证电池的健康状态。说明书中的技术效果和技术问题是示例性的而非限制性的。应当注意，说明书中描述的实施例可以具有其他技术效果并且可以解决其他技术问题。
[0178]
虽然仅示出和描述了某些特征和实施例，但是本领域技术人员可以想到许多修改和变化(例如，各种元件的大小、尺寸、结构、形状和比例的变化，参数的值(例如，温度、压力等)，安装布置，材料的使用，颜色，方向等)，而不实质上脱离所公开的主题的新颖教导和优点。根据替代实施例，可以改变或重新排序任何过程或方法步骤的顺序或次序。因此，应理解，所附权利要求旨在涵盖本公开的真实精神内的所有这些修改和变化。此外，为了提供示例性实施例的简明描述，可能未描述实际实施方式的所有特征。应当理解，在任何这样的实际实施方式的开发中，如在任何工程或设计项目中，可以做出许多实施方式的具体决定。这样的开发努力可能是复杂且耗时的，但是对于受益于本公开的普通技术人员来说，却是设计、生产和制造的常规任务，而无需过多的实验。

技术特征：
1.一种电池系统，其包括：电池，配置为耦接到电气系统；和电池控制模块，配置为电耦接到所述电池，其中所述电池控制模块配置为：监测所述电池的至少一个监测参数；和使用基于至少一个等效电路模型、所述至少一个监测参数和卡尔曼滤波器的两个线性回归模型，递归地计算所述电池的第一容量估计值。2.根据权利要求1所述的电池系统，其中所述两个线性回归模型中的第一线性回归模型配置为计算所述电池的充电状态估计值，其中所述两个线性回归模型中的第二线性回归模型配置为至少部分地基于所述充电状态估计值来计算所述第一容量估计值。3.根据权利要求2所述的电池系统，其中基于所述电池的估计开路电压，通过开路电压
‑
充电状态查找表来计算所述第一线性回归模型的所述充电状态估计值。4.根据权利要求1所述的电池系统，其中所述电池的至少一个监测参数包括端子电压、电流、温度或其任意组合。5.根据权利要求1所述的电池系统，其中所述等效电路模型包括所述电池的1
‑
rc等效电路模型。6.根据权利要求1所述的电池系统，其中所述电池控制模块配置为使用所述电池的两个弛豫开路电压测量值和所述电池的库仑计数来递归地计算所述电池的第二容量估计值。7.根据权利要求6所述的电池系统，其中所述电池控制模块配置为在所述电池超过一小时不活动之后接收所述两个弛豫开路电压测量值。8.根据权利要求6所述的电池系统，其中所述电池控制模块配置为相隔5小时内接收所述两个弛豫开路电压测量值。9.根据权利要求6所述的电池系统，其中所述电池控制模块配置为基于所述两个弛豫开路电压测量值计算两个弛豫充电状态值，其中所述第二容量估计值至少部分地基于所述两个弛豫充电状态值。10.根据权利要求1所述的电池系统，其中所述第一容量估计值存储在环形缓冲器中，并且所述电池控制模块配置为至少部分地基于在所述环形缓冲器的时间上较早的位置中存储的所述第一容量估计值的先前值来确定所述第一容量估计值的收敛。11.根据权利要求1所述的电池系统，其中所述电池控制模块被配置为对所述电池的所述至少一个监测参数使用比例因子，以能够在定点微处理器上递归计算所述电池的第一容量估计值。12.根据权利要求1所述的电池系统，其中所述电池控制模块配置为测试所述电池的第一容量估计值的收敛。13.一种确定配置为耦接到电气系统的可再充电电池的实时容量的方法，其包括：通过耦接到所述可再充电电池的一个或多个传感器，利用电池控制模块监测所述电池工作期间所述可再充电电池的至少一个监测参数；和使用基于至少一个等效电路模型、所述至少一个监测参数和卡尔曼滤波器的两个线性回归模型，通过所述电池控制模块递归地计算所述可再充电电池的第一实时容量。14.根据权利要求13所述的方法，其中所述两个线性回归模型是串行执行的。15.根据权利要求13所述的方法，其中所述两个线性回归模型中的第一回归模型包括
计算所述可再充电电池的充电状态估计值，其中所述两个线性回归模型的第二线性回归模型包括至少部分地基于所述充电状态估计值计算所述第一实时容量。16.根据权利要求13所述的方法，包括至少部分地基于所述可再充电电池的两个弛豫开路电压测量值和所述可再充电电池的库仑计数来递归地计算所述可再充电电池的第二实时容量。17.根据权利要求16所述的方法，其中在所述可再充电电池不活跃一个或多个小时之后，所述两个弛豫开路电压测量值被所述电池控制模块使用。18.一种用于车辆的能量存储部件，其包括：壳体；第一端子和第二端子；可再充电电池，设置在所述壳体内，并配置成与所述第一端子和所述第二端子耦接；电池控制模块，配置为：监测所述能量存储部件的至少一个监测参数；使用基于至少一个等效电路模型、所述至少一个监测参数和卡尔曼滤波器的两个线性回归模型，递归地计算所述能量储存部件的第一容量估计值；以及至少部分地基于所述能量存储部件的两个弛豫开路电压测量值和所述能量存储部件的库仑计数，递归地计算所述能量存储部件的第二容量估计值。19.根据权利要求18所述的能量存储部件，其中所述至少一个监测参数包括端子电压、电流、温度或其任意组合。20.根据权利要求18所述的能量存储部件，其中所述等效电路模型包括1
‑
rc等效电路模型。21.根据权利要求18所述的能量储存部件，其中所述两个线性回归模型中的第一回归模型包括计算所述能量储存部件的第一充电状态估计值，其中所述两个线性回归模型中的第二线性回归模型包括至少部分地基于所述第一充电状态估计值计算所述第一容量估计值。22.根据权利要求18所述的能量储存部件，其中所述电池控制模块配置为基于所述两个弛豫开路电压测量值计算第二充电状态估计值和第三充电状态估计值，其中所述第二容量估计值至少部分地基于所述第二充电状态估计值和第三充电状态估计值。23.根据权利要求22所述的能量储存部件，其中所述电池控制模块配置为确定所述能量储存部件的第一容量估计值的收敛。
技术总结
一种电池系统包括耦接到电气系统的电池。电池系统还包括电耦接到电池的电池控制模块。电池控制模块监测电池的至少一个监测参数，且电池控制模块使用基于至少一个等效电路模型、至少一个监测参数和卡尔曼滤波器的两个线性回归模型，递归地计算电池的第一容量估计值。递归地计算电池的第一容量估计值。递归地计算电池的第一容量估计值。


技术研发人员：郭庆之 克里斯汀
受保护的技术使用者：柯锐世先进解决方案有限责任公司
技术研发日：2016.08.02
技术公布日：2021/6/25