一种锂电池SOH预测方法及装置与流程

一种锂电池soh预测方法及装置
技术领域
1.本发明涉及锂电池健康管理技术领域，特别涉及一种锂电池soh预测方法及装置。


背景技术：

2.锂离子电池是电动汽车、电子消费产品甚至是航天器等的核心电源，因此，锂电池电池的可靠性和安全性是实际应用过程中的一个极为关键的问题。锂离子电池的充放电过程会 导致电池内部产生化学反映，使得锂离子不断流失，电池性能不断下降，这种锂电池的老化过程是不可逆的，这种情况会影响电气设备的正常运行，如电动汽车行驶里程变短，还可能带来严重的后果，例如近年来发生的电池爆炸、电动汽车自燃以及一些发电厂的电池储能箱爆炸，因此，锂电池的健康状态（state of health）成为了一项受人关注的指标。soh是一个表征当前电池相对于新电池存储电能的能力，以百分比的形式表示电池从寿命开始到寿命结束期间所处的状态，用来定量描述当前的性能状态。为了能够尽量避免锂电池soh的衰退对整体造成的影响，对soh的预测工作成为了健康管理（phm）领域的热点课题。soh的预测能够帮助相关技术人员掌握电池老化的影响因素，为电池的使用和维护提供理论指导。
3.目前，针对soh的预测方法大致上分为两类，基于模型的方法和数据驱动的方法。模型的方法运用的较多的是建立电化学模型和等效电路模型，但前者建模过程非常复杂，后者虽然结构简单，但是耗时且需要获得大量不同的数据集，且预测效果较差。数据驱动法包含机器学习、人工智能的算法等。由于电池衰退过程具有一定的规律，选择不同的算法的效果也有很大的区别。普通的数据驱动预测算法如lstm、svm、bp相比传统的如线性回归，对于非线性关系的数据，能够更好的捕捉到数据之间的特性，预测效果有所提高，但依然存在许多问题，lstm由于网络的输出部分是通过随机分配权重的方式，容易丢失掉关键的信息，并且容易受到数据噪声和冗余特征的影响。bp神经网络的泛化能力差，且在训练过程中容易陷入局部最优的情况，布谷鸟算法随让能够避免网络陷入局部最优，但是在输入训练集数量下降的时候，预测模型的精度会受到影响，因此鲁棒性较差。svm在数据集较大时，会耗费大量的机器内存，且运算时间比较长。由上分析，已有的数据驱动的预测算法在网络稳定、抵抗数据噪声和冗余特征影响这几方面的效果较差。容易影响网络的预测精度。在soh预测工作上，现有的方法普遍没有考虑影响电池容量衰退的影响因子，而是直接将soh输入到网络中，实际上，影响soh衰退的因素由很多，这些影响因子或多或少都会对soh的衰退趋势造成影响。


技术实现要素：

4.为解决上述现有技术中所存在的传统soh预测方法所暴露出来的缺点，如稳定性差、泛化能力差、预测精度低等问题，本发明提供了一种锂电池soh预测方法及装置，能够进行预测精度较高、稳定性较强的锂电池soh预测。
5.为了实现上述技术目的，本发明提供了如下技术方案：一种锂电池soh预测方法，
包括：获取锂电池的样本数据，通过灰色关联度分析方法对样本数据进行筛选，得到相关联数据；获取锂电池soh的原始数据，通过集合经验模态分解对原始数据进行分解，生成soh分量训练数据及标签数据；构建深度学习模型，其中所述深度学习模型为基于注意力机制的cnn-bilstm模型；通过相关联数据、soh分量训练数据及标签数据对深度学习模型进行训练，生成优化模型，通过优化模型对实测数据进行预测，生成锂电池soh预测结果以实现锂电池soh预测。
6.可选的，所述样本数据包括soc数据、测量电压、测量电流、工作温度。
7.可选的，对样本数据进行筛选的过程包括：将样本数据与对应soh中对应的故障特征序列分别进行差值计算，统计所述差值计算结果中的最大差值及最小差值，基于最大差值及最小差值计算得到相关系数，并对相关系数进行均值计算，基于均值计算结果对样本数据进行筛选，得到相关联数据。
8.可选的，对原始数据进行分解的过程包括：将白噪声添加到原始数据中，得到原始信号，统计所述原始信号中的局部极值点，通过三次样条差值方法对局部极值点进行处理，得到包络线；计算包络线的均值，对包络线的均值及原始信号进行计算，得到中间信号，对中间信号进行判断，基于判断结果，将所述中间信号作为分量信号，根据原始信号及分量信号计算得到剩余分量，对剩余分量重复计算及判断过程，直到最后一个残余信号为单调函数，得到若干个分量信号，基于若干个分量信号及残余信号，得到分解后的信号，分解后的信号包括soh分量训练数据及标签数据，其中soh分量训练数据为标签数据的上一时刻数据。
9.可选的，构建深度学习模型的过程包括：构建bilstm网络，在所述bilstm网络输入端添加cnn网络，生成cnn-bilstm模型，并在所述bilstm网络的输出端添加注意力机制，并在注意力机制后添加dense层进行映射，生成基于注意力机制的cnn-bilstm模型。
10.可选的，所述cnn-bilstm模型的计算过程为：
11.其中，x
t
表示t时刻的输入，wf、wr、wz表示cnn-bilstm网络遗忘门、输出门、输入门对应的内置权重，ef、ez和er表示网络遗忘门、输入门、输出门对应的偏置。g(t)、j(t)、c
t
、r(t)分别表示的是遗忘门、输入门、记忆储存、输出门，y(t)表示网络的输出，tanh代表双曲线正切函数，kf、kz、kr表示遗忘门、输入门、输出门对应的输出权重值。
12.可选的，所述注意力机制中，通过softmax函数计算注意力机制中的注意力权重。
13.可选的，对深度学习模型进行训练之前还包括：分别对相关联数据进行预处理，其中预处理过程包括归一化处理。
14.为了更好的实现上述技术目的，本发明还提供了一种锂电池soh预测装置，包括：样本获取模块用于获取锂电池的样本数据，通过灰色关联度分析方法对样本数据进行筛选，得到相关联数据；标签获取模块用于获取锂电池soh的原始数据，通过集合经验模态分解对原始数据进行分解，生成标签数据；构建模块用于构建深度学习模型，其中所述深度学习模型为基于注意力机制的cnn-bilstm模型；优化测量模块用于通过相关联数据及标签数据对深度学习模型进行训练，生成优化模型，通过优化模型对实测数据进行预测，生成锂电池soh预测结果以实现锂电池soh预测。
15.本发明具有如下技术效果：本方案相较于cnn-bilstm，在网络框架中添加了注意力机制，能够代替传统bilstm神经元随机分配权重以计算输出值的方式，通过计算结果与该点的误差大小合理分配权重，更能保留数据之前的关键特征。此外，把控了网络的输入部分，利用gra-eemd的方法合理的控制网络的输入数据，避免冗余特征和数据之间存在的噪声对网络的影响，从而整体提高网络的稳定性和泛化能力。
附图说明
16.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
17.图1为本发明实施例提供的方法流程示意图；图2为本发明实施例提供的方案效果示意图。
具体实施方式
18.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
19.如图1所述，本发明提供了一种锂电池soh预测方法及装置，方法所涉及相关技术：注意力机制结合cnn-bilstm模型，此外，在数据处理阶段利用灰色关联度分析（gra）和集合经验模态分解（eemd）的方法对数据进行筛选和降噪处理。其方法步骤包括：步骤1.特征提取：锂电池循环测试数据集中除了容量，普遍包含多种实验过程种采集的其他相关数据，如soc（state of charge）、测量电压、测量电流、工作温度，这些数据称之为相关数据。
soh的计算方法为soh=当前容量/标称容量，在计算出soh后，传统方案普遍是将相关数据包括soh一同输入到网络中，与传统方法不同，本专利采用了gra（gray relation analysis）方法筛选出与soh关联度较好的非冗余特征输入到网络中。
20.灰色关联分析是将soh序列和相关序列通过合理的方法转化成几何形状并通过对比，判断其相似程度来确定相关序列与soh序列之间的关联程度。表1为关联度分析表，由表1可知，当关联系数小于0.6时，相关序列与故障特征属于失调状态，如果此类失调数据作为输入被传输到网络中，会影响模型的精度和稳定性，因此，应该选择关联系数大于0.6的相关序列。
21.表1
22.步骤方法如下：1）将相关序列视为如下矩阵：
23.其中，m代表每个相关序列的数据个数。将所有的相关序列转换成如式（2）所示
24.其中n表示的是相关序列的个数。
25.2) 确定soh序列ri0:
26.3）依次计算每个相关序列与故障特征序列对应位置元素的差值：
27.4）根据3），确定各序列分别与故障特征序列的最大差值和最小差值，然后再从这些差值里选取最值，既确定两级最大差和两级最小差，分别为和。
28.5）根据式5计算每个相关序列和故障特征序列对应的位置的相关系数，为分辨系数，取0.5。
[0029][0030]
6）计算5）中得到的每个相关序列与故障特征序列相关系数，利用式5求取均值，得到的结果既为相关序列与故障特征序列的关联度，由表1可知，关联度越高，证明相关性越强，协同度越高。
[0031][0032]
利用此方法，可以合理的确定网络的输入参数，保证网络的稳定性不受冗余特征的影响。
[0033]
步骤2：数据分解：由于soh数据之间存在数据噪声的情况，此类情况会影响网络的整体性能。本专利使用了集合经验模态分解(ensemble empirical mode decomposition，eemd)将高频和不规则分量从soh中分离出来，可以减少不同数据之间的噪声对模型精度的影响。除此之外，eemd相比emd，过程信号极值点分更加均匀，可以更加有效的抑制由间歇性高频分量等因素造成的模态混叠。步骤方法如下：1）将正态分布的白噪声加到soh中，并且作为一个整体，将其视为原始信号，记为v（t）。
[0034]
2）找出原始信号v(t)中全部局部极值点，并且利用三次样条差值的方法，将局部极大值点和局部极小值点连成上包络线、下包络线，分别为emax(t)，emin(t)，且所有数据点都在上、下包络线构成的范围内部。
[0035]
3）计算上、下包络线的均值
[0036]
4）将原始信号v(t)减去n1(t)，得到中间信号c(t)，既：
[0037]
5）判断信号c(t)是否满足固有模态函数imf的两个条件，分别是：在整个过程中，极值点个数与过零点个数相等或最多相差1；任意时刻内，上、下包络线相对于时间轴局部对称。
[0038]
若满足，则将中间信号c1(t)视为一个imf分量，若不满足，则以该信号为基础，重新循环1）到步骤4），直到得到第一个imf分量。
[0039]
6）将原始信号v(t)减去第一个imf分量c1(t)，得到剩余分量r(t)，既：
[0040]
将剩余分量r(t)重复步骤1到步骤4既可得到第二个imf分量c2(t)，在通过步骤5，如此分解下去，直到最后一个残余信号d(t)成为单调函数无法分解后，既完成了对故障特征的分解工作，分解后的原始信号v(t)如式8所示，所有imf分量与残余信号之和为原始信号v(t)。
[0041][0042]
经eemd分解后得到多个分量，将单个分量与其他特征输入到网络中，能够降低数据噪声对网络的影响，还可以提高网络的泛化能力。
[0043]
步骤3：数据预处理：经步骤一筛选后的非冗余特征数据的按时间顺序将前70%视为训练集，后30%视为测试集。训练集用于网络的训练，测试集用于评估网络的预测精度。为了能够输出数据之间的量纲差异，需要对数据进行预处理，具体的，对除soh数据外获取的数据进行归一化处理，归一化公式如下所示。
[0044][0045]
其中，x’表示样本数据，x
max
表示样本数据的最大值，x
min
表示样本数据的最小值。
[0046]
步骤4：构建cnn-bilstm网络结构lstm网络结构是由循环神经网络演变而来，是为了针对循环神经网络在训练过程中容易出现容易出现的梯度消失或梯度爆炸的情况。lstm与其他神经网络不同的是，它内部包含3个门控模块，分别是遗忘门、输入门和输出门，这三个门控模块的作用分别是：遗忘门：遗忘门可以通过将上一个输入c
t-1
乘以遗忘权重来选择性的忘记。遗忘权重是将h
t-1
和x
t
堆叠后乘以遗忘门的权重wf，再通过σ函数激活，得到的值称之为“遗忘权重”。σ函数的表达式如式9所示。
[0047][0048]
输入门：输入门能够对输入选择性的记忆，将h
t-1
和x
t
堆叠后乘以输入门的权重wi然后使用σ函数激活，得到的权重称之为“输入权重”。
[0049]
输出门：输出门可以通过将当前信息乘以输出权重来决定当前信息哪些可以当成当前状态作为输出，输出权重是将h
t-1
和x
t
堆叠后乘以输出门的权重wo，使用σ函数激活，得到的值称之为“输出权重”。
[0050]
w则是用于循环神经网络中的运算，经过循环神经网络计算得到的权重是通过输入信息乘以w后通过激活函数tanh激活得到的权重。
[0051]
输入信息经过lstm神经元的处理后，得到的结果如式11所示。
[0052][0053]
若将h
t
进行softmax操作，则得到真正意义上的输出既y(t)，否则，它就将被作为下一个神经元的输入。
[0054]
lstm是一个正向的过程，它只考虑上文的信息，如果将输入的序列反转，在重新按照lstm的方式重新计算一遍，最后得到的结果与正向计算得到的结果简单的叠加，考虑正向、反向lstm便是bilstm，通过bilstm更能够捕捉双向数据的变化趋势。
[0055]
cnn主要由卷积层和池化层构成，在bilstm网络中添加卷积神经网络，卷积层能够
提高网络的非线性特征提取能力，池化层通过过滤器的结果max操作，负责压缩提取的特征并且生成更具有特征信息的信号，max可以保留显著的信息，经cnn处理后的结果作为bilstm的输入。经过cnn-bilstm处理后，得到的结果如下公式所示：
[0056]
式中，x
t
表示t时刻的输入，wf、wr、wz表示cnn-bilstm网络遗忘门、输出门、输入门对应的内置权重，ef、ez和er表示网络遗忘门、输入门、输出门对应的偏置。g(t)、j(t)、c
t
、r(t)分别表示的是遗忘门、输入门、记忆储存、输出门，y(t)表示网络的输出，tanh代表双曲线正切函数，kf、kz、kr表示遗忘门、输入门、输出门对应的输出权重值。网络的数据流是bilstm网络通过读取上一时刻的输出y(t)和记忆存储c
t
，通过上述的公式得到g(t)、j(t)、r(t)，随后得到下一时刻的输出和记忆存储。
[0057]
步骤5：网络的输出添加注意力机制。
[0058]
注意力机制被提出的初衷就是为了模仿人脑的资源分配机制，这种机制能够让人将注意力集中在需要关注的地方，从而获得更多的有用信息，忽略无用信息。注意力机制是通过合理分配权重的方法，忽略无用信息，将关注点聚焦在关键部位，从而提高模型的预测精度。在所提出的网络中添加注意力机制，可以有效的避免bilstm在针对长序列预测工作中出现的有效信息丢失的情况。
[0059]
具体的，在bilstm的输出部分添加注意力机制，以此取代bilstm神经元随机分配权值的方式。计算神经元内部隐含层h
i’和前面的随机初始化的四个隐含层hi的相似度ei，公式如式11所示，下一步，使用softmax函数计算出注意力权重，得到注意力权重w后，与神经元的输出加权整合，通过dense层映射，得到预测值output。softmax公式如式12所示。
[0060][0061]
步骤6：模型训练与测试。
[0062]
经gra-eemd处理后，将t-1时刻的soh分量和非冗余特征作为训练数据，t时刻的soh分量作为目标标签，利用时间滑窗制造出训练数据与标签对。标签值即为实际值，用于展示训练过程中预测值与实际值的误差。在cnn-bilstm网络的输出中添加注意力机制，将训练集作为网络的输入，经过训练得到一个预训练模型，随后，将测试集输入到预训练模型中，得到的结果为分量1的预测值，随后，将soh分量2替换分量1，在重复上述训练-测试步骤，最后得到了每个分量的预测值，最后，将多个分量的预测值叠加即为soh预测结果即测试结果。构建预训练模型的流程如图1所示。
[0063]
本方案可用于在线预测或离线分析。在线预测需要先进行离线分析，先通过离线分析电池包的终身数据，经过gra确定非冗余参数。在线预测过程中，根据离线分析得到的经gra处理的结果，选择满足要求的参数作为网络的输入。根据预训练模型的建立时长，如预训练模型建立需要4分钟，那么选择输送到网络中的时间片段为五分钟，既每五分钟，数据输入到网络中，由于预训练模型建立后无需重新建立，因此，可以在未来五分钟的数据传入到网络之前得到预测结果。综上所述，根据预训练模型的建立时长，延长一定的时间确定每个输入到网络中的时间片段，确保下次数据传输之前得到预测结果，实现在线应用。
[0064]
根据本方案阐述的步骤，采用nasa锂电池公开数据集作为验证，利用该数据集中的容量参数计算soh，计算公式如下。该数据集包含多个参数，根据方案的步骤，先利用gra分析各参数与soh的关联度，最后发现只有工作温度和工作电压与soh的关联度超过0.6，其它参数都属于失调状态，因此，将这两类参数作为网络的输入。随后将工作温度、工作电压和soh划分成训练集和测试集，前70%为训练集，后30%为测试集。将工作温度和工作电压进行归一化处理后，将soh进行eemd分解，得到10个分量和1个残余信号。将soh序列替换成分量1序列，三类参数输入到添加了注意力机制的cnn-bilstm网络中，通过训练集对网络进行训练，经过训练得到预训练模型后，将测试集输入到预训练模型中，得到分量1的预测结果。后续将分量1替换成分量2，重复上述步骤，得到分量2的预测结果。重复该步骤，直到10个分量和残余信号都完成了上述步骤，最终得到了10个分量和残余信号的预测值，将这些预测值叠加即为soh预测结果即测试集所生成结果。如图2所示，可以发现，本方案相较于cnn-bilstm，更能够追踪soh的变化曲线。为了能够直观的看出模型之间的差异，利用公式15.16和17计算预测值与实际值的误差。式中，xi表示实际值，yi表示预测值，误差越小代表预测精度越高。计算误差如表2所示。表2为误差指标表，由表2所示，本方案在三个误差指标分别为0.0832、0.0991和4.67%上，相比cnn-lstm的误差都更低，因此本方案的预测精度相较于cnn-lstm更高。
[0065]
在实际运用过程中，除循环次数外，其他参数如温度，工作电压等参数都具有一定的周期性，既重复性，只有循环次数是基于上次的次数逐渐叠加，因此，在已有的数据中，未来时刻的除循环次数以外的参数与训练集中的数据相同。以nasa数据集为例，选择训练集的工作电压和工作温度作为未来该电池的这两类参数，循环次数基于之前的基础上叠加，由此构建了未来时刻的众参数即实测数据，将其导入模型中即可得到未来时刻的soh。
[0066][0067][0068]
表2
[0069]
实施例二为了更好的实现上述技术目的，本发明还提供了一种锂电池soh预测装置，包括：样本获取模块用于获取锂电池的样本数据，通过灰色关联度分析方法对样本数据进行筛选，得到相关联数据；标签获取模块用于获取锂电池soh的原始数据，通过集合经验模态分解对原始数据进行分解，生成标签数据；构建模块用于构建深度学习模型，其中所述深度学习模型为基于注意力机制的cnn-bilstm模型；优化测量模块用于通过相关联数据及标签数据对深度学习模型进行训练，生成优化模型，通过优化模型对实测数据进行预测，生成锂电池soh预测结果以实现锂电池soh预测。本装置与上述方法相对应，此处不做赘述。
[0070]
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。