面向多元锂离子电池数据的全生命周期演化趋势预测方法

1.本发明涉及锂离子电池领域，特别是涉及一种面向多元锂离子电池数据的全生命周期演化趋势预测方法。


背景技术：

2.锂离子电池作为最成熟的储能技术，是电动车辆和储能系统(ess)最核心的组成部分，因而面向蓄能产品的故障预测和健康管理(prognostics and health management,phm)对于运行维护至关重要。从电池全生命周期管理的角度来讲，锂离子电池的电化学机理在使用或搁置过程会逐渐退化，使系统性能的供电能力和安全特性等管理指标也表现为衰减趋势，为降低电池运行成本、实现及时维护和故障防控，基于电池管理系统(bmss)对健康状态(soh)进行在线监控是很有必要的。从电池管理的角度来讲，故障预测的核心是针对电池系统运行欠佳或失效的异常进行刨析，健康管理的核心是面向容量和阻抗等指标的健康状态(state of health,soh)估计。
3.目前，基于模型驱动的健康状态预测占多数。该类方法是根据有穷小规模数据的背景分布或领域经验设定数学假设，通过假设建模最优拟合度的映射函数，并基于函数运算值实现预测。基于模型驱动的健康状态预测方法更易于调整模型参数、适应性更强，针对相似机理或同等工况的模型具备普适性。然而，现有驱动模式的时序预测方法存在难以挖掘现在和未来间细粒度映射关系、难以提升智能预测服务准确性、鲁棒性和可解释性的技术问题。


技术实现要素：

4.针对现有技术中的上述不足，本发明提供的一种面向多元锂离子电池数据的全生命周期演化趋势预测方法解决了现有预测方法难以挖掘现在和未来间细粒度映射关系、难以提升智能预测服务准确性、鲁棒性和可解释性的技术问题。
5.为了达到上述发明目的，本发明采用的技术方案为：一种面向多元锂离子电池数据的全生命周期演化趋势预测方法，所述方法包括以下步骤：
6.s1：针对产业链采集的锂离子电池数据集，基于相邻样本数为6的knn优化算法和调整平滑参数的vondrak滤波算法实现缺失值填补和序列平滑降噪的预处理操作，并根据训练、评估、测试比率为6:2:2进行集合划分；
7.s2：根据经预处理操作的锂离子电池数据集td，利用特征再抽象算法，基于最小绝对偏差和非局部季节滤波改进robuststl算法，并基于时序数据的加法概念模型对辨识性特征进行低耦合的属性分解获得低阶特征分量；
8.s3：基于低阶特征分量，设置噪音剔除评价指标和关联强度评估算法分别判定低阶特征分量的白噪声高频占比以及seasonality分量、trend分量与数据曲线的关联强度，实现锂离子电池的属性降维和关联增强的特征提取；
9.s4：基于属性降维和关联增强的特征提取，利用矩阵点积运算构建时间戳强约束
的时序预测训练集合；
10.s5：基于时序预测训练集合，利用基于自注意加权算法的lstm底层逻辑设计lbtsa深度模型，通过深度模型的链接权重和学习参数实现网络建模，根据集合划分设置机理属性和强约束时间戳作为深度模型的控制输入，获取当前时刻的soh状态和未来时间区间的soh状态估计值，实现锂离子电池的全生命周期演化趋势的预测。
11.上述方案的有益效果是：本发明基于数据驱动的大数据智能分析框架提出了针对多元长周期锂离子电池序列的soh状态预测算法，是以fra再抽象序列算法为核心的特征提取和以lbtsa深度模型为主体的数据建模。解决了现有预测方法难以挖掘现在和未来间细粒度映射关系、难以提升智能预测服务准确性、鲁棒性和可解释性的技术问题。
12.进一步地，s1中利用knn优化算法，根据欧氏距离确定距离缺失样本点接近的6个样本，并通过计算6个样本观测点的加权平均值估计缺失数据。
13.上述进一步方案的有益效果是：序列缺失数值将激化数据稀疏性敏感问题，价值挖掘会因空值现象而陷入混乱，导致不可靠信息的输出，导致时序预测表现的不确定性因素更加突出，蕴含的确定性因素更难把握，采用上述技术方案，有效解决序列缺失问题。
14.进一步地，s1中vondrak滤波算法q公式为
[0015][0016]
其中，f为逼真拟合度，λ为正常值，s为粗糙光滑度，n为数据规模，i为数据规模取值，pi为观测点的数值权重，yi为观测数值，为观测值的平均值，δ为差分符号，δ3为三阶差分；
[0017]
设定无量纲正数ε为
[0018][0019]
将无量纲正数ε定义为平滑因子，在观测数据的绝对逼真和绝对拟合之间起到平滑作用。
[0020]
上述进一步方案的有益效果是：针对产业链数据集存在噪音干扰及离群点密集的预处理问题，本发明通过调参优化了vondrak滤波机制。其中，vondrak滤波是面向航天高精度序列数据的处理算法，可以剔除高频噪音并最大程度保留价值区间，通过定义平滑因子，在观测数据的绝对逼真和绝对拟合之间起到平滑作用。
[0021]
进一步地，s2中改进的robuststl算法通过设置内循环与外循环实现鲁棒局部加权回归的权重调节、趋势拟合及周期分量的系统运算，使经改进的robuststl算法分解的特征分量对异常波动和趋势突变具有强的鲁棒性；
[0022]
所述改进的robuststl算法包括基于白噪音平滑的residual更新阶段，基于序列差分的trend分解阶段和基于领域覆盖的seasonality分解阶段。
[0023]
上述进一步方案的有益效果是：通过对robuststl算法进行改进，提升算法的鲁棒性，利于方案的实施。
[0024]
进一步地，基于白噪音平滑的residual更新阶段包括以下公式：
[0025]
根据序列加法概念模型将residual分量描述为
[0026]
residuali＝ai+ni[0027]
其中，residuali为由白噪声ni、尖峰ai或低谷构成的残差集合；
[0028]
利用双边滤波算法剔除不稳定因素以保证residual分量和trend分量的稳健性，所述双边滤波算法公式如下：
[0029][0030]
其中，td'i为去噪样本，j为长度为2h+1的滤波窗口，为时刻i的观测点邻接区间为j的滤波权重；
[0031][0032]
其中，σ为无量纲的归一化相关参数，e为指数函数，和为算法预设的平滑度控制参数，tdj为邻接区间为j的观测值，tdi为时刻i的观测值；
[0033]
则更新后的一阶residual分量为
[0034][0035]
其中，residual'i为更新后的一阶residual分量，为趋势trendi中白噪音的错误或干扰因素。
[0036]
上述进一步方案的有益效果是：通过上述公式，用于实现基于白噪音平滑的residual更新阶段，residual表现为采集工况中因信号错误或噪声污染掺杂的不稳定因素，通过利用双边滤波技术剔除不稳定因素以保证trend分量和seasonality分量的稳健性。
[0037]
进一步地，基于序列差分的trend分解阶段包括以下公式：
[0038]
通过设定基于seasonality分量和residual分量的季节性差分值趋于极小值的假设条件减轻季节性特征在序列成分中的影响因子，季节性差分gi公式为
[0039][0040]
其中，t为循环周期，代表差分操作，seasonalityi为seasonality分量；
[0041]
基于季节性差分gi构建加权目标函数恢复一阶加权目标函数公式为
[0042][0043]
其中为季节性二阶差分值，λ1和λ2为目标函数的权衡参数，n为观测序列的样本点数量；
[0044]
通过加权目标函数，得到trend分量公式为
[0045][0046]
其中，为trend分量的优化输出，上标～为优化输出的标识，为首次
迭代输出的趋势分量，x为fra算法的迭代次数；
[0047]
更新后的二阶residual分量residual”i
为
[0048][0049]
其中，trendi为trend分量；
[0050]
基于lad损失回归函数分解trend分量后，去趋势集合td”i
表示为
[0051][0052]
上述进一步方案的有益效果是：通过上述公式用于实现基于序列差分的trend分解阶段。
[0053]
进一步地，基于领域覆盖的seasonality分解阶段中采用非局部季节性过滤算法公式为
[0054][0055][0056]
ω＝{(i',j)|(i'＝t-k
×
t,j＝i'
±
h)}k＝1,2,
…
,k；h＝1,2,
…
,h
[0057]
其中，为长度为2h+1的滤波窗口的邻接区间为(i
′
，j)的滤波权重，ω为邻接区间的集合表示，tdj″
为邻接区间(i
′
，j)时的去趋势集合，t、k、k、h和h均为正实数。
[0058]
上述进一步方案的有益效果是：通过上述公式用于实现基于领域覆盖的seasonality分解阶段，为确保seasonality的估算精度，本方案采用非局部季节性过滤算法。
[0059]
进一步地，s3中包括以下公式
[0060]
面向序列白噪音分解纯度的评估阶段，针对特征分量的二阶抽象运算函数tsr表示为
[0061]
tsr＝fra(trend，seasonality，residual|θ，μ，td，ψ)
[0062]
其中，fra为再抽象算法，θ为预设定的提取纯度，μ为白噪音分解纯度，ψ为fra的算法参数，td为观测值；
[0063]
判定白噪音分解占比μ
′n的公式为
[0064][0065]
其中，n为实验序列长度，μ
′n为第n次算法迭代时residual分量相较第n-1次residual分量的收敛性，residual
i,n
为第n次residual分量，residual
i,n-1
为第n-1次residual分量；
[0066]
通过比较μ
′n与预设值，评估白噪声的高频占比；
[0067]
面向trend分量和seasonality分量提取纯度的评估阶段，利用pearson相关系数
衡量特征分量和评估参数定距变量间的线性关系，通过健康指标soh判定特征分量提取纯度的公式为
[0068][0069][0070]
其中，和为不同序列在第n次迭代时的总体相关评估强度，sohi为第i个时间戳的健康指标值，为trend的平均值，为soh的平均值。
[0071]
上述进一步方案的有益效果是：通过上述技术方案，获得白噪声高频占比以及seasonality分量、trend分量和数据曲线的关联强度，实现序列演化规律和波动频率的无监督学习。
[0072]
进一步地，s4中通过搭建矩阵运算逻辑，构建时间戳强约束的时序预测训练集合，包括以下公式：
[0073]
根据季节趋势分量，通过矩阵点积运算构建一阶训练集tsf'为
[0074][0075][0076]
其中，tsr为特征分量，p为2
×
3阶矩阵，上标t为矩阵的转置，
⊙
为同或运算符；
[0077]
锂离子电池数据集合中，基于序列周期cycle和采集频率af构建的时间戳time位于维度为13的第12列，基于实际电容c构建的soh评估参数位于第13列，自注意力机制在信息加权时要设定时间特性的数学假设，其中，一对一映射的时间戳time标签列能够表现时间依赖性，soh标签设定为建模的目标参数，在训练阶段用于深度模型的负反馈优化，测试阶段用于预测性能评估，提取参数的矩阵点积公式为
[0078][0079]
其中，q为2
×
13阶矩阵；
[0080]
基于一阶训练集tsf'，通过矩阵点积运算构建二阶训练集tsf”为
[0081][0082]
上述进一步方案的有益效果是：通过上述技术方案，利用矩阵点积运算构建时序预测训练集合，点积(多重复制)注意力函数，其计算规则与传统self-attention的运算规则一致，在实践中，点积注意力的运行时效性表现很优秀，速度更快，空间效率更高。
[0083]
进一步地，s5中设计lbtsa深度模型利用嵌套的建模架构机制，融合self-attention机制的序列加权及lstm神经网络的远距离信息继承，实现针对lstm单元门控机制和seq2seq架构的优化，所述嵌套的建模架构机制，基于lstm神经网络的seq2seq架构实现，包括以下分步骤：
[0084]
s5-1：将lstm神经网络的概念化模型拆分为编码器和解码器，利用lstm encoder对特征训练集tsf进行序列编码，利用门控输出针对时间戳time、trend和seasonality列进行计算，隐藏状态针对soh列进行计算，得到全局时间步的序列输出和最后时间步的隐藏状态，包括以下公式：
[0085]
采用点积注意力加权算法，并构建了基于点积注意力加权算法的lstm解码器模型lbtsa，所述点积注意力加权算法构建query，key，value矩阵向量为
[0086][0087]
其中，query、keyi和valuei为抽象向量，hd为隐藏状态，qi为第i个观测值的特征向量，ωq、ωk和ωv为神经网络的参数；
[0088]
s5-2：利用self-attention算法建立序列输出和隐藏状态间的自注意机制，基于序列输出对隐藏状态关联度的加权操作，得到融合自注意力机制的序列矩阵，包括以下公式：
[0089]
基于抽象向量query，key，value构建注意力得分矩阵g为
[0090][0091]
其中，dk为key的维度；
[0092]
为缩放因子，防止内积数值过大影响lbtsa深度模型的无监督自适应学习机制；
[0093]
针对各时间步特征向量构建的注意力权重矩阵a为
[0094][0095]
其中，gi为时间步为i时注意力得分矩阵的列向量，exp为针对注意力得分矩阵的运行以e为底的指数运算；
[0096]
lbtsa深度模型的self-attention算法构建加权序列z2(q,k,v)为
[0097][0098]
s5-3：利用lstm decoder对序列矩阵进行解码，得到预测序列。
[0099]
上述进一步方案的有益效果是：通过上述技术方案，解决了结构化长周期序列预测的数据稀疏性敏感问题，实现了长序预测的性能优化。
附图说明
[0100]
图1为一种面向多元锂离子电池数据的全生命周期演化趋势预测方法流程图。
[0101]
图2为特征再抽象算法执行框图。
[0102]
图3为lbtsa深度模型的原理图。
具体实施方式
[0103]
下面结合附图和具体实施例对本发明做进一步说明。
[0104]
如图1所示，一种面向多元锂离子电池数据的全生命周期演化趋势预测方法，其特征在于，所述方法包括以下步骤：
[0105]
s1：针对产业链采集的锂离子电池数据集，基于相邻样本数为6的knn优化算法和调整平滑参数的vondrak滤波算法实现缺失值填补和序列平滑降噪的预处理操作，并根据训练、评估、测试比率为6:2:2进行集合划分；
[0106]
s2：根据经预处理操作的锂离子电池数据集td，利用特征再抽象算法，基于最小绝对偏差和非局部季节滤波改进robuststl算法，并基于时序数据的加法概念模型对辨识性特征进行低耦合的属性分解获得低阶特征分量；
[0107]
s3：基于低阶特征分量，设置噪音剔除评价指标和关联强度评估算法分别判定低阶特征分量的白噪声高频占比以及seasonality分量、trend分量与数据曲线的关联强度，实现锂离子电池的属性降维和关联增强的特征提取；
[0108]
s4：基于属性降维和关联增强的特征提取，利用矩阵点积运算构建时间戳强约束的时序预测训练集合；
[0109]
s5：基于时序预测训练集合，利用基于自注意加权算法的lstm底层逻辑设计lbtsa深度模型，通过深度模型的链接权重和学习参数实现网络建模，根据集合划分设置机理属性和强约束时间戳作为深度模型的控制输入，获取当前时刻的soh状态和未来时间区间的soh状态估计值，实现锂离子电池的全生命周期演化趋势的预测。
[0110]
s1中利用knn优化算法，根据欧氏距离确定距离缺失样本点接近的6个样本，并通过计算6个样本观测点的加权平均值估计缺失数据。
[0111]
s1中vondrak滤波算法q公式为
[0112][0113]
其中，f为逼真拟合度，λ为正常值，s为粗糙光滑度，n为数据规模，i为数据规模取值，pi为观测点的数值权重，yi为观测数值，为观测值的平均值，δ为差分符号，δ3为三阶差分；
[0114]
设定无量纲正数ε为
[0115][0116]
将无量纲正数ε定义为平滑因子，在观测数据的绝对逼真和绝对拟合之间起到平滑作用。
[0117]
f表示逼真拟合度，是指滤波曲线和待平滑样本曲线的偏离程度，f的设定数值越小意味拟合程度越理想，当f＝0时说明滤波曲线完全重合观测序列，称为“绝对逼真”；s表示粗糙光滑度，是指滤波曲线和待平滑样本曲线的拟合程度，s的设定数值越小意味曲线光滑程度越明显，当s＝0时说明滤波平滑值位于穿越观测点的直线上，称为“绝对平滑”。
[0118]
s2中改进的robuststl算法通过设置内循环与外循环实现鲁棒局部加权回归的权重调节、趋势拟合及周期分量的系统运算，使经改进的robuststl算法分解的特征分量对异常波动和趋势突变具有强的鲁棒性；
[0119]
所述改进的robuststl算法包括基于白噪音平滑的residual更新阶段，基于序列差分的trend分解阶段和基于领域覆盖的seasonality分解阶段。
[0120]
residual是由通俗意义上的白噪音和尖峰或低谷构成的残差集合，表现为采集工况中因信号错误或噪声污染掺杂的不稳定因素。本发明利用双边滤波技术剔除不稳定因素以保证trend分量和seasonality分量的稳健性，基于相似的邻接观测值平滑序列数据，完整保留趋势trend的突变、尖峰或低谷。
[0121]
基于白噪音平滑的residual更新阶段包括以下公式：
[0122]
根据序列加法概念模型将residual分量描述为
[0123]
residuali＝ai+ni[0124]
其中，residuali为由白噪声ni、尖峰ai或低谷构成的残差集合；
[0125]
利用双边滤波算法剔除不稳定因素以保证residual分量和trend分量的稳健性，所述双边滤波算法公式如下：
[0126][0127]
其中，td'i为去噪样本，j为长度为2h+1的滤波窗口，为时刻i的观测点邻接区间为j的滤波权重；
[0128][0129]
其中，σ为无量纲的归一化相关参数，e为指数函数，和为算法预设的平滑度控制参数，tdj为邻接区间为j的观测值，tdi为时刻i的观测值；
[0130]
则更新后的一阶residual分量为
[0131][0132]
其中，residual'i为更新后的一阶residual分量，为趋势trendi中白噪音的错误或干扰因素。
[0133]
基于序列差分的trend分解阶段包括以下公式：
[0134]
通过设定基于seasonality分量和residual分量的季节性差分值趋于极小值的假设条件减轻季节性特征在序列成分中的影响因子，季节性差分gi公式为
[0135][0136]
其中，t为循环周期，代表差分操作，seasonalityi为seasonality分量；
[0137]
基于季节性差分gi构建加权目标函数恢复一阶加权目标函数公式为
[0138][0139]
其中为季节性二阶差分值，λ1和λ2为目标函数的权衡参数，n为观测序列的样本点数量；
[0140]
通过加权目标函数，得到trend分量公式为
[0141][0142]
其中，为trend分量的优化输出，上标～为优化输出的标识，为首次迭代输出的趋势分量，x为fra算法的迭代次数；
[0143]
更新后的二阶residual分量residual”i
为
[0144][0145]
其中，trendi为trend分量；
[0146]
基于lad损失回归函数分解trend分量后，去趋势集合td”i
表示为
[0147][0148]
在本发明的一个实施例中，设计目标函数通过最小化lad损失回归函数的经验误差(加权目标函数公式的第一项)恢复一阶季节差分项，使对序列的不稳定因素具备鲁棒性。加权目标函数公式的第二项表示trend分量的一阶差分算子约束，尽管上述公式设定假设，但依旧表现出突发性的水平波动。加权目标函数公式的第三项表示trend分量的二阶差分算子约束，当序列趋势保持平滑且分段线性，呈现稀疏性的假设条件时，trend分量可以捕获突发性移动和渐变的序列表现。
[0149]
本发明建模seasonality为time
″i的权重线性组合，time
″i是去趋势集合td
″i基于邻域的时间戳。time
″i和td
″i的权值比重既依赖于锂离子电池数据集在time标签维度上的聚类中心，又与time
″i和td
″i的相似程度关联。即non-local seasonality算法基于time
″i季节相似程度为td
″
i-kt
邻接序列的td
″
i-kt+h
点赋予权重，使季节性检索算法高效的识别相似点，解决了季节性偏移问题，针对序列异常值提升鲁棒性。
[0150]
基于领域覆盖的seasonality分解阶段中采用非局部季节性过滤算法公式为
[0151][0152][0153]
ω＝{(i',j)|(i'＝t-k
×
t,j＝i'
±
h)}k＝1,2,
…
,k；h＝1,2,
…
,h
[0154]
其中，为长度为2h+1的滤波窗口的邻接区间为(i
′
，j)的滤波权重，ω为邻接区间的集合表示，td
″j为邻接区间(i
′
，j)时的去趋势集合，t、k、k、h和h均为正实数。
[0155]
s3中包括以下公式
[0156]
面向序列白噪音分解纯度的评估阶段，针对特征分量的二阶抽象运算函数tsr表示为
[0157]
tsr＝fra(trend，seasonality，residual|θ，μ，td，ψ)
[0158]
其中，fra为再抽象算法，θ为预设定的提取纯度，μ为白噪音分解纯度，ψ为fra的算法参数，td为观测值；
[0159]
判定白噪音分解占比μ
′n的公式为
[0160][0161]
其中，n为实验序列长度，μ
′n为第n次算法迭代时residual分量相较第n-1次residual分量的收敛性，residual
i,n
为第n次residual分量，residual
i,n-1
为第n-1次residual分量；
[0162]
通过比较μ
′n与预设值，评估白噪声的高频占比；
[0163]
面向trend分量和seasonality分量提取纯度的评估阶段，利用pearson相关系数衡量特征分量和评估参数定距变量间的线性关系，通过健康指标soh判定特征分量提取纯度的公式为
[0164][0165][0166]
其中，和为不同序列在第n次迭代时的总体相关评估强度，sohi为第i个时间戳的健康指标值，为trend的平均值，为soh的平均值。
[0167]
pearson相关性系数是一个无量纲的统计指标，形式化的描述为x和y两个变量的协方差与标准差乘积的比值，取值范围为-1≤ρ≤1，当ρ≤0时为负相关，当ρ≥0时为正相关，且|ρ|的数值越接近1，表示序列变量关联程度越紧密。依据pearson相关性系数对运算结果的评估，当|ρ|∈[0.8,1]区间时，即证实tsf向量和预期目标参数之间具备极强的相关性，能够通过相关强度匹配深度学习算法的复杂程度。
[0168]
s4中通过搭建矩阵运算逻辑，构建时间戳强约束的时序预测训练集合，包括以下公式：
[0169]
经所述s2和s3中提取的n
×
3特征训练集tsr，其列标签依次对应trend分量、
seasonality分量和residual分量，根据无量纲特征的序列关联性可知，trend分量、seasonality分量是通过波动方向和幅度，以及波动的周期拐点反映序列集合的演变特征，residual分量和集合的关联因素只包含序列白噪音，所以构建时序预测的训练集合仅需要季节趋势分量；
[0170]
根据季节趋势分量，通过矩阵点积运算构建一阶训练集tsf'为
[0171][0172][0173]
其中，tsr为特征分量，p为2
×
3阶矩阵，上标t为矩阵的转置，
⊙
为同或运算符；
[0174]
锂离子电池数据集合中，基于序列周期cycle和采集频率af构建的时间戳time位于维度为13的第12列，基于实际电容c构建的soh评估参数位于第13列，自注意力机制在信息加权时要设定时间特性的数学假设，其中，一对一映射的时间戳time标签列能够表现时间依赖性，soh标签设定为建模的目标参数，在训练阶段用于深度模型的负反馈优化，测试阶段用于预测性能评估，提取参数的矩阵点积公式为
[0175][0176]
其中，q为2
×
13阶矩阵；
[0177]
基于一阶训练集tsf'，通过矩阵点积运算构建二阶训练集tsf”为
[0178][0179]
s5中设计lbtsa深度模型利用嵌套的建模架构机制，融合self-attention机制的序列加权及lstm神经网络的远距离信息继承，实现针对lstm单元门控机制和seq2seq架构的优化，所述嵌套的建模架构机制，基于lstm神经网络的seq2seq架构实现，包括以下分步骤：
[0180]
s5-1：将lstm神经网络的概念化模型拆分为编码器和解码器，利用lstm encoder对特征训练集tsf进行序列编码，利用门控输出针对时间戳time、trend和seasonality列进行计算，隐藏状态针对soh列进行计算，得到全局时间步的序列输出和最后时间步的隐藏状态，包括以下公式：
[0181]
采用点积注意力加权算法，并构建了基于点积注意力加权算法的lstm解码器模型lbtsa，所述点积注意力加权算法构建query，key，value矩阵向量为
[0182][0183]
其中，query、keyi和valuei为抽象向量，hd为隐藏状态，qi为第i个观测值的特征向量，ωq、ωk和ωv为神经网络的参数；
[0184]
s5-2：利用self-attention算法建立序列输出和隐藏状态间的自注意机制，基于
序列输出对隐藏状态关联度的加权操作，得到融合自注意力机制的序列矩阵，包括以下公式：
[0185]
基于抽象向量query，key，value构建注意力得分矩阵g为
[0186][0187]
其中，dk为key的维度；
[0188]
为缩放因子，防止内积数值过大影响lbtsa深度模型的无监督自适应学习机制；
[0189]
针对各时间步特征向量构建的注意力权重矩阵a为
[0190][0191]
其中，gi为时间步为i时注意力得分矩阵的列向量，exp为针对注意力得分矩阵的运行以e为底的指数运算；
[0192]
lbtsa深度模型的self-attention算法构建加权序列z2(q,k,v)为
[0193][0194]
s5-3：利用lstm decoder对序列矩阵进行解码，得到预测序列。
[0195]
本发明提出的lbtsa深度模型采用的注意力函数为点积(多重复制)注意力函数，其计算规则与传统self-attention的运算规则一致，缩放比例因子采用在实践中，点积注意力的运行时效性表现很优秀，速度更快，空间效率更高。
[0196]
在本发明的一个实施例中，基于数据驱动的大数据智能分析框架提出了针对多元长周期锂离子电池序列的soh状态预测算法，是以fra再抽象序列算法为核心的特征提取和以lbtsa深度模型为主体的数据建模。本发明针对产业链采集的工业数据，设置训练、测试及评估的比例为6：2：2，设置机理属性和强约束时间戳作为状态预测模型的输入控制，以获取当前时刻的soh状态和未来时间区间的soh状态估计值，实现面向锂离子电池的故障预测和健康管理的智能化控制。本发明根据锂离子电池安全状态风险评估白皮书预设了soh标准，记作soh_sta。面向锂离子电池故障预测和风险评估的管理系统是基于当前时刻的soh状态估计值和soh_sta的判定结果，当前健康状态满足soh_sta时，说明锂离子电池出现热短路、化学自燃和高温自燃的风险较低，反之亦然。面向锂离子电池健康管理和剩余使用寿命(rul)评估的管理系统是基于未来时间区间的soh状态估计值和soh_sta扼值的判定结果，健康状态的演化曲线趋近soh_sta的程度反映了rul的极值。
[0197]
本发明提出了基于robuststl优化分解和双机制关联评估的特征再抽象(fra)算法，如图2所示，可以异步抽取高精度、弱噪音及低耦合的涨幅趋势、季节波动和残差分量，并利用白噪音分解纯度判定和pearson关联强度判定迭代更新提取分量，解决了结构化高维序列预测的数据稀疏性敏感问题，实现了由多元标签属性序列向二元特征向量的关联性降维映射，在数据基础上保障整体预测性能的准确性和鲁棒性；提出了基于transformer架构下自注意力机制优化的lbtsa深度模型，如图3所示，可以利用self-attention机制的信息加权特性强调序列的价值片段并针对遗忘门的底层逻辑嵌入权重决策机制，结合lstm单
元掌握长距离依赖、选择性遗忘的信息继承能力，解决了结构化长周期序列预测的数据稀疏性敏感问题，实现了长序预测的性能优化；本发明基于“数据+模型”的数据驱动架构提出面向多元锂离子电池数据的全生命周期演化趋势预测方法，通过设置序列分解和关联评估以捕获低维强关联的可辨识特征，实现全广域、深层次的信息挖掘；通过优化lstm网络底层的信息继承逻辑以构建lbtsa深度模型，实现全周期、高性能的未知推测，进一步提升状态预测方法挖掘细粒度映射关系的能力以及预测服务的准确性、鲁棒性及可解释性。
[0198]
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在发明的保护范围内。

技术特征：
1.一种面向多元锂离子电池数据的全生命周期演化趋势预测方法，其特征在于，所述方法包括以下步骤：s1：针对产业链采集的锂离子电池数据集，基于相邻样本数为6的knn优化算法和调整平滑参数的vondrak滤波算法实现缺失值填补和序列平滑降噪的预处理操作，并根据训练、评估、测试比率为6:2:2进行集合划分；s2：根据经预处理操作的锂离子电池数据集td，利用特征再抽象算法，基于最小绝对偏差和非局部季节滤波改进robuststl算法，并基于时序数据的加法概念模型对辨识性特征进行低耦合的属性分解获得低阶特征分量；s3：基于低阶特征分量，设置噪音剔除评价指标和关联强度评估算法分别判定低阶特征分量的白噪声高频占比以及seasonality分量、trend分量与数据曲线的关联强度，实现锂离子电池的属性降维和关联增强的特征提取；s4：基于属性降维和关联增强的特征提取，利用矩阵点积运算构建时间戳强约束的时序预测训练集合；s5：基于时序预测训练集合，利用基于自注意加权算法的lstm底层逻辑设计lbtsa深度模型，通过深度模型的链接权重和学习参数实现网络建模，根据集合划分设置机理属性和强约束时间戳作为深度模型的控制输入，获取当前时刻的soh状态和未来时间区间的soh状态估计值，实现锂离子电池的全生命周期演化趋势的预测。2.根据权利要求1所述的面向多元锂离子电池数据的全生命周期演化趋势预测方法，其特征在于，所述s1中利用knn优化算法，根据欧氏距离确定距离缺失样本点接近的6个样本，并通过计算6个样本观测点的加权平均值估计缺失数据。3.根据权利要求2所述的面向多元锂离子电池数据的全生命周期演化趋势预测方法，其特征在于，所述s1中vondrak滤波算法q公式为其中，f为逼真拟合度，λ为正常值，s为粗糙光滑度，n为数据规模，i为数据规模取值，p
i
为观测点的数值权重，y
i
为观测数值，为观测值的平均值，δ为差分符号，δ3为三阶差分；设定无量纲正数ε为将无量纲正数ε定义为平滑因子，在观测数据的绝对逼真和绝对拟合之间起到平滑作用。4.根据权利要求3所述的面向多元锂离子电池数据的全生命周期演化趋势预测方法，其特征在于，所述s2中改进的robuststl算法通过设置内循环与外循环实现鲁棒局部加权回归的权重调节、趋势拟合及周期分量的系统运算，使经改进的robuststl算法分解的特征分量对异常波动和趋势突变具有强的鲁棒性；所述改进的robuststl算法包括基于白噪音平滑的residual更新阶段，基于序列差分的trend分解阶段和基于领域覆盖的seasonality分解阶段。5.根据权利要求4所述的面向多元锂离子电池数据的全生命周期演化趋势预测方法，
其特征在于，所述基于白噪音平滑的residual更新阶段包括以下公式：根据序列加法概念模型将residual分量描述为residual
i
＝a
i
+n
i
其中，residual
i
为由白噪声n
i
、尖峰a
i
或低谷构成的残差集合；利用双边滤波算法剔除不稳定因素以保证residual分量和trend分量的稳健性，所述双边滤波算法公式如下：其中，td
i
'为去噪样本，j为长度为2h+1的滤波窗口，为时刻i的观测点邻接区间为j的滤波权重；其中，σ为无量纲的归一化相关参数，e为指数函数，和为算法预设的平滑度控制参数，td
j
为邻接区间为j的观测值，td
i
为时刻i的观测值；则更新后的一阶residual分量为其中，residual
′
i
为更新后的一阶residual分量，为趋势trend
i
中白噪音的错误或干扰因素。6.根据权利要求5所述的面向多元锂离子电池数据的全生命周期演化趋势预测方法，其特征在于，所述基于序列差分的trend分解阶段包括以下公式：通过设定基于seasonality分量和residual分量的季节性差分值趋于极小值的假设条件减轻季节性特征在序列成分中的影响因子，季节性差分g
i
公式为其中，t为循环周期，代表差分操作，seasonality
i
为seasonality分量；基于季节性差分g
i
构建加权目标函数恢复一阶加权目标函数公式为其中为季节性二阶差分值，λ1和λ2为目标函数的权衡参数，n为观测序列的样本点数量；通过加权目标函数，得到trend分量公式为其中，为trend分量的优化输出，上标～为优化输出的标识，为首次迭代输出的趋势分量，x为fra算法的迭代次数；
更新后的二阶residual分量residual
″
i
为其中，trend
i
为trend分量；基于lad损失回归函数分解trend分量后，去趋势集合td”i
表示为7.根据权利要求6所述的面向多元锂离子电池数据的全生命周期演化趋势预测方法，其特征在于，所述基于领域覆盖的seasonality分解阶段中采用非局部季节性过滤算法公式为公式为ω＝{(i',j)|(i'＝t-k
×
t,j＝i'
±
h)}k＝1,2,
…
,k；h＝1,2,
…
,h其中，为长度为2h+1的滤波窗口的邻接区间为(i
′
，j)的滤波权重，ω为邻接区间的集合表示，td
″
j
为邻接区间(i
′
，j)的去趋势集合，t、k、k、h和h均为正实数。8.根据权利要求7所述的面向多元锂离子电池数据的全生命周期演化趋势预测方法，其特征在于，所述s3中包括以下公式面向序列白噪音分解纯度的评估阶段，针对特征分量的二阶抽象运算函数tsr表示为tsr＝fra(trend，seasonality，residual|θ，μ，td，ψ)其中，fra为再抽象算法，θ为预设定的提取纯度，μ为白噪音分解纯度，ψ为fra的算法参数，td为观测值；判定白噪音分解占比μ
′
n
的公式为其中，n为实验序列长度，μ
′
n
为第n次算法迭代时residual分量相较第n-1次residual分量的收敛性，residual
i，n
为第n次residual分量，residual
i，n-1
为第n-1次residual分量；通过比较μ
′
n
与预设值，评估白噪声的高频占比；面向trend分量和seasonality分量提取纯度的评估阶段，利用pearson相关系数衡量特征分量和评估参数定距变量间的线性关系，通过健康指标soh判定特征分量提取纯度的公式为
其中，和为不同序列在第n次迭代时的总体相关评估强度，soh
i
为第i个时间戳的健康指标值，为trend的平均值，为soh的平均值。9.根据权利要求8所述的面向多元锂离子电池数据的全生命周期演化趋势预测方法，其特征在于，所述s4中通过搭建矩阵运算逻辑，构建时间戳强约束的时序预测训练集合，包括以下公式：根据季节趋势分量，通过矩阵点积运算构建一阶训练集tsf'为根据季节趋势分量，通过矩阵点积运算构建一阶训练集tsf'为其中，tsr为特征分量，p为2
×
3阶矩阵，上标t为矩阵的转置，
⊙
为同或运算符；锂离子电池数据集合中，基于序列周期cycle和采集频率af构建的时间戳time位于维度为13的第12列，基于实际电容c构建的soh评估参数位于第13列，自注意力机制在信息加权时要设定时间特性的数学假设，其中，一对一映射的时间戳time标签列能够表现时间依赖性，soh标签设定为建模的目标参数，在训练阶段用于深度模型的负反馈优化，测试阶段用于预测性能评估，提取参数的矩阵点积公式为其中，q为2
×
13阶矩阵；基于一阶训练集tsf'，通过矩阵点积运算构建二阶训练集tsf”为10.根据权利要求9所述的面向多元锂离子电池数据的全生命周期演化趋势预测方法，其特征在于，所述s5中设计lbtsa深度模型利用嵌套的建模架构机制，融合self-attention机制的序列加权及lstm神经网络的远距离信息继承，实现针对lstm单元门控机制和seq2seq架构的优化，所述嵌套的建模架构机制，基于lstm神经网络的seq2seq架构实现，包括以下分步骤：s5-1：将lstm神经网络的概念化模型拆分为编码器和解码器，利用lstm encoder对特征训练集tsf进行序列编码，利用门控输出针对时间戳time、trend和seasonality列进行计算，隐藏状态针对soh列进行计算，得到全局时间步的序列输出和最后时间步的隐藏状态，
包括以下公式：采用点积注意力加权算法，并构建了基于点积注意力加权算法的lstm解码器模型lbtsa，所述点积注意力加权算法构建query，key，value矩阵向量为其中，query、key
i
和value
i
为抽象向量，h
d
为隐藏状态，q
i
为第i个观测值的特征向量，ω
q
、ω
k
和ωv为神经网络的参数；s5-2：利用self-attention算法建立序列输出和隐藏状态间的自注意机制，基于序列输出对隐藏状态关联度的加权操作，得到融合自注意力机制的序列矩阵，包括以下公式：基于抽象向量query，key，value构建注意力得分矩阵g为其中，d
k
为key的维度；为缩放因子，防止内积数值过大影响lbtsa深度模型的无监督自适应学习机制；针对各时间步特征向量构建的注意力权重矩阵a为其中，g
i
为时间步为i时注意力得分矩阵的列向量，exp为针对注意力得分矩阵的运行以e为底的指数运算；lbtsa深度模型的self-attention算法构建加权序列z2(q,k,v)为s5-3：利用lstm decoder对序列矩阵进行解码，得到预测序列。

技术总结
本发明公开了一种面向多元锂离子电池数据的全生命周期演化趋势预测方法，根据时序数据的无量纲演化特性对产业链采集的锂离子电池数据利用FRA特征再抽象算法，制作时间戳强约束、关联特性增强的待检测低阶特征集合，利用自注意力机制表达序列型数据高价值区间，利用权值矩阵优化LSTM神经网络的底层门控逻辑设计LBTSA深度模型，利用预训练的链接权重和学习参数实现网络建模，通过LBTSA深度模型挖掘特征集合与预设指标间的映射关系，解决产业链针对多元全周期时序预测提出的性能需求问题。题。题。


技术研发人员：周建涛 王昊 王燕 倪磊 任思宇
受保护的技术使用者：内蒙古大学
技术研发日：2023.05.23
技术公布日：2023/8/21