一种基于红外图像储能电池堆三维温度可视化方法与流程

本发明属于储能电池堆热管理领域，为电池堆的热管理提供一种三维温度可视化的方法。

背景技术：

近年来，锂离子电池因具有功率/能量密度高、循环寿命长、自放电率低等优点，受到学术界和产业界的广泛关注，已应用于电子消费品、电动汽车、分布式储能、大规模储能等不同场景。然而受相关领域的安全性能升级需求影响，作为这些领域动力系统的核心部件，锂离子电池的安全性受到了行业内外的广泛关注。特别是在大规模储能应用领域，当锂离子电池发生热失控，引发起火、爆炸等事故时，整个储能电站将毁于一旦，而且对电站周边环境、公众的安全与财产产生一定的负面效应。国内外近期发生多起锂离子电池储能电站火灾事故，2018年7月2日，韩国一风力发电园区内ess储能设备发生重大火灾事故，造成706m2规模电池建筑和3500块以上锂电池全部烧毁。

锂离子电池储能电站的安全问题是需要警钟长鸣的重大课题。随着锂离子电池新材料的研发、电池制作技术的创新以及众多科研机构和企业的参与，锂离子电池的性能正日益提高，单体安全性能也得到极大提高。但由于大规模储能系统单体电池容量更大，电池簇单体数量更多，电池簇并联数量更大，电池堆电流更大，电池簇充放电深度更深，电池簇运行一致性和寿命要求更为严格，在使用过程中极易出现局部热失控现象，存在巨大的安全隐患。

在日常生活、工业制造等众多领域，温度的测量与控制时刻都在进行，因此温度测量应用极其广泛，在温度测量技术上国内外众多的研究人员也开展了大量的研究工作。当前，比较常用的温度测量方法主要分为两大类:接触式测量和非接触式测量。接触式测量应用最为广泛的是热电阻、热电偶测温，两者都是点式测温，通过与被测物体单点接触进行测温，只能反映物体某一点的温度。红外诊断技术属于非接触故障，能够实现不停机、不断电进行检测故障，能够及时准确地发现热故障，具有检测灵活方便、安全、检测范围广等优点。近几年的计算机技术和微电子技术的迅速发展，推进了红外诊断技术的发展。如今红外诊断技术己经比较成熟，检测具有较高的准确度，检测灵活方便，加上先进的图像处理技术和科学的诊断算法，可更加智能化对电池堆热故障进行定性和定位分析。

技术实现要素：

为了克服现有的分布式光伏电站储能系统热故障检测的准确率较低、检测麻烦的不足，本发明提供一种具有较高的准确度、检测灵活方便的基于红外图像的电池堆三维立体温度场重构方法，该方法结合电池堆表面温度分布情况、特定点温度以及电池阻抗，快速重构电池堆三维温度场。

本发明解决其技术问题所采用的技术方案是：

一种储能电池堆三维温度可视化方法，所述方法包括步骤如下：

1)红外摄像辐射定标

以面源黑体为基准，采集红外热像仪在不同黑体温度下的灰度图像，然后利用图像的灰度和黑体实际温度进行拟合，再辅以环境辐射修正等措施，建立温度与灰度映射模型temval＝m(grayval)，其中temval表示映射的温度值，grayval表示红外图像灰度值，m()表示拟合的非线性方程；

2)获取电池堆表面温度

首先对获取的红外图像预处理，通过小波阈值函数滤除图像的噪声；接着与可见光图像配准，结合图像分割与形态学技术分割出红外电池堆表面；最后结合电池堆节点分布将电池堆表面进行子单元划分，在各个表面子单元内连续插值，对表面子单元内所有插值点积分并除以面积得到相应灰度值，辅以辐射理论对其矫正，根据温度与灰度映射模型temval＝m(grayval)得到电池堆表面温度；

3)基于反距离权重插值的三维温度可视化模型

首先基于反距离权重插值原理，以电池堆表面温度插值得到电池堆三维温度temv1(x,y,z)；接着根据n个温度传感器的温度误差插值得到电池堆三维温度误差temv2(x,y,z)，该分量与temv1(x,y,z)叠加，即得到修正的电池堆三维温度模型temv3(x,y,z)；最后，预留的温度传感器的输出用于评价电池堆三维温度模型temv3(x,y,z)的精度；

4)依据阻抗修正立体子单元内部温度

首先以所用电池型号为基准，采集其不同内芯温度、不同激励频率以及不同电荷状态下的电池阻抗，建立电池阻抗与温度映射模型imp＝g(f,tin,soc)，其中imp表示电池阻抗，f表示激励频率，tin表示电池内芯温度，soc表示电荷状态，imp＝g()表示拟合的非线性方程；接着根据建立的阻抗与温度映射模型估算电池堆立体子单元内部温度最后以获取的立体子单元内部温度修正三维温度模型。

进一步，所述步骤2)中，电池堆表面子单元以积分方式进行灰度计算，而不是采取简单的均值，以其中一个表面为例，即有：

式(1)中y，z表示坐标，y1，y1表示y轴起点和终点，z1，z2表示z轴起点和终点，fg1(y,z)表示主视灰度分布情况，f(y,z)表示插值点灰度。

再进一步，所述步骤3)中，基于反距离权重插值的三维温度可视化步骤为：

3.1)根据步骤2)得到的电池堆表面温度分布反距离权重插值电池堆三维温度temv1(x,y,z)，即有：

temv1(x,y,z)＝λ1tout1(y,z)+λ2tout2(x,z)+λ3tout3(x,y)(2)

式(2)、(3)中x，y，z表示坐标，temv1(x,y,z)表示立体子单元(x,y,z)的温度，tout1(y,z)，tout1(x,z)，tout1(x,y)表示不同表面子单元的温度，d1，d2，d3分别表示立体子单元距离(y,z)，(x,z)，(x,y)表面的垂直距离，λ1，λ2，λ3分别表示不同表面子单元对其立体子单元的权重，p表示指数值，默认值为－2；

3.2)以n个采集点的温度误差进行反距离权重插值得到电池堆三维温度误差temv2(x,y,z)，即有

式(4)、(5)temv2(x,y,z)表示立体子单元(x,y,z)处的温度误差，δti(x,y,z)表示采集点i的温度误差；n表示采集点的数量，λi表示采集点i的权重，di表示插值点与采集点i处的直线距离；

3.3)用步骤3.2)计算得到的temv2(x,y,z)补偿步骤3.1)计算得到的temv1(x,y,z),即有修正的电池堆三维温度模型temv3(x,y,z)；

3.4)预留的温度传感器的输出对步骤3.3)重构的三维温度可视化模型以均方误差mse为评价准则进行评价。

更进一步，所述步骤4)中，在soc未知的情况下，电池阻抗与温度映射模型imp＝g(f,tin,soc)通过对soc平均来构造阻抗，以使模型独立于这些影响，即有：

式(6)中n1表示soc的取值数量，n2表示相等soc下测量的次数，vj表示零均值高斯噪声，g(f,tin)表示在soc未知情况下，阻抗与温度映射模型；

所述步骤4)中，为确保电池阻抗模型相对于温度的灵敏度，需要确定激励频率f或多个频率fi,i∈{1,...,n}，从而得出估算电池温度的估算即有：

g1(fi,tin,zi)＝re(g(fi,tin)-zi)(8)

g2(fi,tin,zi)＝im(g(fi,tin)-zi)(9)

式(7)、(8)、(9)中n表示测量过程中激励频率f的数量，zi表示测量阻抗，表示在笛卡尔坐标系下的超参数，

本发明的有益效果为：结合电池堆表面温度分布情况、特定点温度以及电池阻抗，快速重构电池堆三维温度场；具有较高的准确度、检测灵活方便。

附图说明

图1是一种基于红外图像储能电池堆三维温度可视化方法流程图。

具体实施方式

结合附图对本发明的实施例作详细说明：本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。

参照图1，一种基于红外图像储能电池堆三维温度可视化方法，所述方法包括步骤如下：

1)红外摄像辐射定标

以面源黑体为基准，采集红外热像仪在不同黑体温度下的灰度图像，然后利用图像的灰度和黑体实际温度进行拟合，再辅以环境辐射修正等措施，建立温度与灰度映射模型temval＝m(grayval)，其中temval表示映射的温度值，grayval表示红外图像灰度值，m()表示拟合的非线性方程；

2)获取电池堆表面温度

首先对获取的红外图像预处理，通过小波阈值函数滤除图像的噪声；接着与可见光图像配准，结合图像分割与形态学技术分割出红外电池堆表面；最后结合电池堆节点分布将电池堆表面进行子单元划分，在各个表面子单元内连续插值，对表面子单元内所有插值点积分并除以面积得到相应灰度值，辅以辐射理论对其矫正，根据温度与灰度映射模型temval＝m(grayval)得到电池堆表面温度；

3)基于反距离权重插值的三维温度可视化模型

首先基于反距离权重插值原理，以电池堆表面温度插值得到电池堆三维温度temv1(x,y,z)；接着根据n个温度传感器的温度误差插值得到电池堆三维温度误差temv2(x,y,z)，该分量与temv1(x,y,z)叠加，即得到修正的电池堆三维温度模型temv3(x,y,z)；最后，预留的温度传感器的输出用于评价电池堆三维温度模型temv3(x,y,z)的精度；

4)依据阻抗修正立体子单元内部温度

首先以所用电池型号为基准，采集其不同内芯温度、不同激励频率以及不同电荷状态下的电池阻抗，建立电池阻抗与温度映射模型imp＝g(f,tin,soc)，其中imp表示电池阻抗，f表示激励频率，tin表示电池内芯温度，soc表示电荷状态，imp＝g()表示拟合的非线性方程；接着根据建立的阻抗与温度映射模型估算电池堆立体子单元内部温度最后以获取的立体子单元内部温度修正三维温度模型。

进一步，所述步骤2)中，电池堆表面子单元以积分方式进行灰度计算，而不是采取简单的均值，以其中一个表面为例，即有：

式(1)中y，z表示坐标，y1，y1表示y轴起点和终点，z1，z2表示z轴起点和终点，fg1(y,z)表示主视灰度分布情况，f(y,z)表示插值点灰度。

再进一步，所述步骤3)中，基于反距离权重插值的三维温度可视化步骤为：

3.1)根据步骤2)得到的电池堆表面温度分布反距离权重插值电池堆三维温度temv1(x,y,z)，即有：

temv1(x,y,z)＝λ1tout1(y,z)+λ2tout2(x,z)+λ3tout3(x,y)(2)

式(2)、(3)中x，y，z表示坐标，temv1(x,y,z)表示立体子单元(x,y,z)的温度，tout1(y,z)，tout1(x,z)，tout1(x,y)表示不同表面子单元的温度，d1，d2，d3分别表示立体子单元距离(y,z)，(x,z)，(x,y)表面的垂直距离，λ1，λ2，λ3分别表示不同表面子单元对其立体子单元的权重，p表示指数值，默认值为－2；

3.2)以n个采集点的温度误差进行反距离权重插值得到电池堆三维温度误差temv2(x,y,z)，即有

式(4)、(5)temv2(x,y,z)表示立体子单元(x,y,z)处的温度误差，δti(x,y,z)表示采集点i的温度误差；n表示采集点的数量，λi表示采集点i的权重，di表示插值点与采集点i处的直线距离；

3.3)用步骤3.2)计算得到的temv2(x,y,z)补偿步骤3.1)计算得到的temv1(x,y,z),即有修正的电池堆三维温度模型temv3(x,y,z)；

3.4)预留的温度传感器的输出对步骤3.3)重构的三维温度可视化模型以均方误差mse为评价准则进行评价。

更进一步，所述步骤4)中，在soc未知的情况下，电池阻抗与温度映射模型imp＝g(f,tin,soc)通过对soc平均来构造阻抗，以使模型独立于这些影响，即有：

式(6)中n1表示soc的取值数量，n2表示相等soc下测量的次数，vj表示零均值高斯噪声，g(f,tin)表示在soc未知情况下，阻抗与温度映射模型；

所述步骤4)中，为确保电池阻抗模型相对于温度的灵敏度，需要确定激励频率f或多个频率fi,i∈{1,...,n}，从而得出估算电池温度的估算即有：

g1(fi,tin,zi)＝re(g(fi,tin)-zi)(8)

g2(fi,tin,zi)＝im(g(fi,tin)-zi)(9)

式(7)、(8)、(9)中n表示测量过程中激励频率f的数量，zi表示测量阻抗，表示在笛卡尔坐标系下的超参数，

最后，还需要注意的是，以上列举的仅是本发明的一个具体实施例。显然，本发明不限于以上实施例，还可以有许多变形。本领域的普通技术人员能从本发明公开的内容直接导出或联想到的所有变形，均应认为是本发明的保护范围。