0.   引言

变电站电力设备长期处于工况运行状态，易发生各类故障而威胁电网安全稳定运行。研究显示^[1]，在电气设备故障中超过半数均伴有设备异常发热现象。

20世纪中期国外的一些研究机构与电力企业尝试使用热像仪记录设备温度进行故障检测^[2]。发展至今，国内外的技术人员通过对运行设备定期热拍照，由检测者观察设备热图像的温度分布并判定设备中存在的异常温升，确定热故障类型^[3]，该方法可有效检测出电力设备的异常温升。加上其无接触检测的优点，适应电力系统要求，因而在电力系统故障检测中得到推广运用。然而当前电力设备热图谱普遍采用人工观察，判断较主观^[4-5]，且其准确程度对检测人员经验依赖大，准确度波动大且不够高效，并且此类传统的人工图像特征分析方法较难提炼异常发热区域特征与分布。近年来计算机技术和人工智能技术飞速发展，研究者们提出了一些具有自动化和智能化功能的设备热像图故障识别方法，以提高判定与故障处理效率与准确性。对于自动化算法来说，如何利用电力设备的各类特点敏锐地发现热像图中的异常温升区域是热故障识别的关键。胡世征等人通过检测设备对应测量点间温差与其中较热点温升的比值，实现相对温差判断^[6]。郭英军等提出使用红外摄像仪在线监测与计算变压器铁心的温升变化情况^[7]。王如意等根据设备热图谱在异常温升处的色差进行图像分割而提取故障点^[8]。孙怡等人提出利用不同电气设备出现异常温升的经验数值进行分类阈值判定的方法^[9]。王淼等人通过分析导线热像图诊断输电线路^[10]。冯振新等采用局部区域聚类的方法对热图谱温升区进行提取^[11]。而现有研究较多关注于故障点程度的分割、计算和分析工作，而运用三相设备自有特征进行故障区域自动查找定位的研究却较少涉及。

在目前现状下，本文提出一种三相自搜寻比较的电力设备热像图自动诊断算法，变电站三相设备或并列运行的同类设备通常具有完全相同的外观特点，其电气运行状态相同或相近（严重三相不平衡和短路故障情况除外），因而正常运行时三相设备一般具有相同的热状态，对应部位的发热状态一致，这是热检测的重要判据。当其中某相设备发生局部热故障时，正常相设备可成为其余两相设备发热异常判定的良好参照物。通过计算机自动搜寻算法提取并利用三相对应的热状态信息对比，可自动高效地发现热故障并准确判定其发生部位，进行诊断分析。

1.   三相自搜寻比较算法

三相自搜寻比较法的思路是以热点区域为追踪目标，对三相设备热像图进行计算机处理，自动化地得到可直接相比较分析的三相图片。其方法实施流程如图 1所示。

图  1  三相自搜寻比较法执行流程图

Figure  1.  Three-phase self - search comparison execution flowchart

下载: 全尺寸图片 幻灯片

为实现自动热区域搜寻和比较，本路线中先利用数学形态学的方法对图像中视图相连的三相设备进行分割处理，进一步将每相提取为单独图片，另外由于取景角度和杂物的影响，需对单相图像进行图像调整使其保持相同图像大小、角度，通过这样的自动化处理就得到了可进行对应区域映射对比的单相图片。接着为了最优程度地进行对应处比较，采用数据模块化方法，将处理完成的单相图像按照同等规格进行区块划分，由于先前的单相图同等化处理和同规格的分区，每相图片相同坐标的区块温度很大程度上反应了三相设备相同部位的温度水平。对每区温度取平均值，并对三相设备全区块进行温度分布统计，可很好地反应出各相的温度水平和温升情况，基于此相同的数据模型，可以便捷地实现设备温升判定和热数据信息分析，完成自动温升搜寻和故障比较识别。

2.   三相设备热像图的提取分割

2.1   热像图形态学处理

实现三相设备或同类设备之间的热像图对比，需要将图片中三相设备按相进行分离。由于图像中存在连接体、环境异物影响，以及拍摄角度不规范导致设备局部重叠等原因，致使照片中的设备间常有连接物存在，增加了三相设备主体形域的提取分割难度，如图 2所示。

图  2  断路器三相设备连接图及其二值化

Figure  2.  Circuit breaker three-phase equipment connection diagram and its binarization

下载: 全尺寸图片 幻灯片

图 2为电力断路器三相设备热像图及其提取的二值化图像。热像图中最亮处温度达68.1℃，最暗处为－0.2℃。从图中可看出，对图像二值化后，这些干扰一般体现为设备主体形状间的纤细连接，为去除这些连接干扰，可以采用对三相设备热像图的二值图像进行形态学开操作运算。

形态学是由动植物形态研究发展而来。数学形态学的处理可实现从图像中提取边界、骨架等表达描绘图形的图像分量，常用操作有腐蚀膨胀，开闭运算等。

形态学上称类似的图形间的纤细连接为不可靠连接，可以使用形态腐蚀膨胀等操作将其切断去除。其中，在二维图像Z²（即一个二维整数空间）以一个一定像素大小的结构元B对图形对象A进行空间卷积，$\hat B$为B的补集，B对A腐蚀定义为：

而膨胀则是：

可得出，腐蚀操作可消除运算对象接触的所有边界点，使其边界向内收缩，该过程能够消除目标对象附近小而无意义的物体（可干扰体）^[2]，而膨胀操作则相反，其使目标区域范围“变大”。将目标区域接触的背景点合并到该目标物中，使目标边界向外部扩张^[2]，可以用来填补目标区域中某些空洞以及消除包含在目标区域中的小颗粒噪声。

在此基础上，利用腐蚀和膨胀操作可以衍生出开启运算：

即对图形对象进行先腐蚀后膨胀的操作。

以及闭合运算：

即对图形对象进行先膨胀后腐蚀的操作。它们对图形操作具有不同的效果，都能部分平滑对象轮廓，开操作断开较窄的狭颈并消除细的突出物，闭操作则一般弥合较窄的间断和细长的沟壑，消除小的孔洞，填补轮廓线中的断裂等^[12]。

结合三相设备热像图及其得到的二值图像的特点，可对其分别进行合适的开闭类操作，实现对不可靠连接干扰的处理。

如图 3所示，可看出开运算对短小纤细处等不可靠连接进行了“切断”处理，并去除了周边的小干扰体，而闭运算弥合了狭窄的目标间缝隙，使得目标扩大连接在一起。

图  3  开运算和闭运算操作

Figure  3.  Opening and closing operations

下载: 全尺寸图片 幻灯片

针对去除不可靠连接，分离三相设备主体形态的目的，可使用合适程度的开运算操作，现采用合适大小的圆状结构元，对图 4所示三相断路器设备的热像图进行开运算操作。

图  4  三相设备开运算

Figure  4.  Opening operation for three-phase equipment

下载: 全尺寸图片 幻灯片

可以看出，开运算操作对三相设备热像图二值图像中不可靠连接具有较好的“切断”作用。且开运算运用像素半径不同的结构元执行，将得到有区别的执行效果，结构元尺寸减小，图像将由从连接处“猛烈截断”到保留了更多区域细节，应根据图像情况调整合适的结构元大小。

2.2   三相设备主体的标定与提取

经过形态处理后还需进一步将三相设备主体提取出来，观察形态处理后的二值图像，可以看出白色目标区域主要是设备主体区域以及无用的杂散对象，在图像中，把这些具有相同像素值且相互邻接的一片像素区域作为连通区域或连通分量，在二值图像中即前景白色的数个区域。为进一步提取设备区域，可将图像中的连通区域分配不同数字标签并标定。根据其分布特性，将每个白色（像素值1）的连通分量像素标为同一值，不同连通分量分别标定1, 2, 3…等标签，形成新的标签矩阵，并将其用不同颜色表示，如图 5所示。

图  5  连通分量的标定

Figure  5.  Calibration of connected components

下载: 全尺寸图片 幻灯片

根据分析与标签结果，图像中标定的区域主要是设备主体区域以及附带设备、背景或未完全去除的线管等。依据开运算操作时处理尺度不同，标定区域的大小与数量也不同，而在所有的标签区域中，需要提取三相设备主体所在的连通分量。

根据设备热像图的拍摄特点分析，设备本体是拍摄的主要目标，对合格的设备热像图来说，电气设备主体应占据照片主要区域。利用这个特点，三相设备热像图中，标定连通区域中像素数量最大的3个区域即是三相设备的主体区域，即设备目标区域。

如图 6所示，搜集所有标定区域的像素数与标签号，找到具有前三多像素数的区域标签，并将其标记到新的矩阵图像中并进行颜色标定。可以看出，图像中标定了三相设备主体所在的最大的3个连通区域，去除了其他干扰区域的影响。

图  6  三相主体区域提取

Figure  6.  Three phase main area extraction

下载: 全尺寸图片 幻灯片

因二值图像目标区域像素值为1，而背景为0，进一步可将新的三相区域二值矩阵与原设备热像图执行对应元素相乘，如图 7所示，能够得到去除所有干扰与背景元素的三相设备热像对比图。

图  7  三相设备主体目标热图像

Figure  7.  Thermal image of main target of three-phase equipment

下载: 全尺寸图片 幻灯片

2.3   设备图像分相及调整

得到三相设备主体的热图像后，为实现每相设备之间的热状态对比，应将其每相分离至单张图片上，并裁切调整至一致。由于每相所在的连通区域已标定不同标签，能够确定标签及所在位置，因此可为每相区域创建相同大小图像，将三相区域分别标定，即能得到分离的三相设备热图像。

将每相设备图像如图 8单独分离后，由于其大小有差异以及背景位置不同，为进行对应位置的比对，应将其裁切并调整为相同图片大小与方向。利用每相二值图像记录的设备区域位置裁切，并可通过图像分辨率调整将其调整为等大图像，如图 9所示。

图  8  设备热图像多相分离

Figure  8.  Device thermal image multiphase separation

下载: 全尺寸图片 幻灯片

图  9  单相设备热像图调整

Figure  9.  Thermal image adjustment of single-phase equipment

下载: 全尺寸图片 幻灯片

3.   三相设备区块对比与温升区域判定

采用同类比较法，判定设备某区域是否出现异常温升，在完成了设备提取和图像调整后，三相热像图应是大小一致的图像。热像图中由颜色亮暗代表温度高低，由比色条色彩标定，由此将热像图矩阵转变为对应像素的温度分布矩阵，即可进行温度对比。可将彩色热像图的RGB像素值转换为反应亮暗程度的灰度图像：

由热像图比色条中给出的图像温度区间[T_min, T_max]，将灰度图像素值矩阵映射至该温度区间，即：

由此可得相同大小设备温度分布矩阵。而根据电气设备的故障发热特点，某处的故障可能会造成该处一片区域发热而产生异常温升，因此，可将每相的温度分布矩阵划分为数个区块，每个区块温度取其均值，以反映每处的温度分布水平。如图 10是实验设备三相的分区块均温分布表。

图  10  三相分区温度分布表（℃）

Figure  10.  Three-phase partition temperature distribution table (℃)

下载: 全尺寸图片 幻灯片

根据各个区块的平均温度值分布情况，还可对各相不同温度范围内的区块数量进行统计，更直观地分析其温度分布情况，如图 11将示例三相断路器的每相的温度分布的区块数量以直方图统计。由图可看出，排除背景的低温区块，该设备温度区块主要分布在30℃~40℃温度范围内。而在更高的温度层级，B相、C相有少量区块的温度跨过了一般范围直接达到60℃~70℃温度范围，相对大多数区块所在温度范围属于较明显的温升，而A相则相对正常，无少量区块相较有明显的温度跃升。

图  11  每相各温度区块数量分布直方图

Figure  11.  Histogram of quantity distribution of each temperature block in each phase

下载: 全尺寸图片 幻灯片

除讨论不同温度范围的区块数量外，各区块相对位置的温度分布情况是判定温升区域的重点。为了能够清晰地看出每相设备的区块温度分布情况，使用常规的红外图像比色条对温度分布表中每区块温度进行映射，如图 12所示。

图  12  热像图比色条

Figure  12.  Thermal image colorimetric bar

下载: 全尺寸图片 幻灯片

将所有区块按温度均值映射回红外图像标记色。如图 13可见，可清晰获取感知每区块的均值温度，而按红外图像温度颜色标记的区块图可清晰地体现区块的温度水平及每相对应区块的温差。为更清晰体现温度比对水平，可进一步将红外比色区块图制为柱状图，如图 14。

图  13  三相设备区块对比

Figure  13.  Three-phase equipment block comparison

下载: 全尺寸图片 幻灯片

图  14  三相设备温度对比柱状图

Figure  14.  Three phase equipment temperature comparison bar chart

下载: 全尺寸图片 幻灯片

至此，每相设备分区块的温度对比即已清晰明了，只要根据设定的阈值标准，即可判定异常温升的存在，在热成像电力设备故障检测领域中，即可根据同类比较，如同组三相设备、同相设备之间及同类设备之间对应部位的温差进行比较分析，这类判定方法中的相对温差法是一种良好的根据周边温度差值进行温升判定的检测方法^[6]，利用相对温差δ_t即同类设备对应测量点间的温差与其中较高温点温升的比值：

式中：τ₁和T₁为发热点的温升和温度；τ₂和T₂为正常相对应点的温升和温度；T₀为环境温度参照体的温度。以示例三相断路器图谱分析，设定异常温升参考值为20 K，进行区块温差比较计算发现，B相区块中部和C相顶部连接处出现异常温度高值，计算相对四周的δ_t过大、超限，结果符合图 15的温度区块数分布图的统计分析结果。由此能判定B相和C相均存在区域发生了超限温升热故障。

图  15  异常温升区块判定

Figure  15.  Determination of abnormal temperature rise block

下载: 全尺寸图片 幻灯片

如图 15所示，准确地分析出设备存在的发热区域，是快速识别诊断热故障的关键。运用该方法可有效实现三相设备异常温升区域的自动查找。图 16和图 17分别为用该方法实现两组不同角度拍摄的三相电流互感器设备发热故障搜寻效果图。

图  16  三相自搜寻比较法发热查找

Figure  16.  Three-phase self-search comparative method heat search

下载: 全尺寸图片 幻灯片

图  17  三相自搜寻比较法发热查找

Figure  17.  Three-phase self-search comparative method heat search

下载: 全尺寸图片 幻灯片

基于此，三相自搜寻比较法按照图 1本文方法流程对电力设备热像图进行内部比较判断，自动分析识别，从而能更有效地识别设备温升与热故障。

该方法适用于具有成套相同设备、多相并列运行等易于进行相间温度比较的电力设备或者相内温升明显易确定相对温差的设备。若针对拍摄角度随意的多相设备热图谱，其记录的各相热信息缺失，较难准确地对各部位的温升情况做出故障判定，另外对于形态较为独特或不连续的设备，其判定能力可能受限。因此本方法对热像图拍摄质量有一定要求，后续研究可通过积累大量热故障判定实例来提高多相判定的准确性和对各种热图谱的兼容性，以提升方法智能化水平。

4.   结语

以电力设备红外热成像检测中同类比较法与相对温差法为基础，本文提出了一种基于电气设备热信息的三相自搜寻比较检测算法。对电网常见的三相设备、同类设备的检测热像图，本算法针对三相设备热像图的故障显示特征，基于图像处理和温差分析实施对热故障的自动搜寻与处理，获得三相温升比较对象，并进行分区块温度分布统计与自动温升比较分析，实现以同类比较法和相对温差法为分析基础的电力设备热故障高效自动化检测。与传统的设备热像图人工检测分析方法相比，本方法能够实现设备热故障的计算机自动搜寻检测，相对人工检测具有更高的效率和准确度，且能够实现三相设备热故障的全自动化搜寻检测。不足之处是对热像图拍摄角度、质量有一定要求，工程实际中需要对红外检测作业人员的拍摄角度进行必要规范，形成符合要求的三相设备热像图。