本发明涉及空冷器积灰状态感知技术，特别是涉及一种空冷器积灰状态感知系统与计算方法。
背景技术
空冷凝汽器是蒸汽动力循环发电机组的重要环节。在1939年以前均采用水冷，他的工作状态直接影响了整个机组的安全性和经济性。无论是当前占全国总发电量73.1％的火电机组，预计到2020年和2030年仍将占60％和48％，还是将还是将要大规模发展的太阳能热发电机组，其所需的燃煤和太阳能资源主要分布在我国的东北、西北和华北(“三北”)地区，然而“三北”地区却是水资源贫乏的缺水或少水地带，导致燃煤和太阳能资源与水资源分布存在地理结构上的矛盾，但空冷器以热容量小很多的空气作为冷却介质，引起冷端散热能力锐减，致使发电热耗高，经济性差。为了强化空气侧换热过程，空冷器的基管上布置翅片，这在我国干燥、风沙大、扬尘多的“三北”地区，空气所携带的悬浮颗粒物极易在翅片管上积聚，生成灰垢(积灰层)。灰垢的存在降低了空冷器的换热性能，提高了运行成本，可使同环境条件下机组排汽压力抬高8～12kpa，增大发电煤耗约12～18g/kwh，又增加了直冷机组的安全隐患。
从火电厂的建设上来分析，由于我国的水资源匮乏，节约生产、生活用水作为设计的首要任务，因此在火电站大力的发展空冷技术更好解决这一问题。但空冷换热器积灰是难以避免的问题，尤其是温度变化大且风沙多的地区，受积灰的影响、负荷变化大时气流冲击以及其他因素引起的振动等原因，空冷换热器的换热效果常常受到影响而降低甚至会发生故障。因为空冷岛的积灰会造成换热系数降低、垢层厚度增大，所以空冷岛的积灰处理非常重要，但在空冷凝汽器运行方面，现场的运行基本靠经验来判断灰垢的积聚状态以及对汽轮机组排汽压力和发电量的影响，至于空冷器灰垢的维护清洗也完全是凭借经验，导致清洗不及时，灰垢积聚严重，危害机组安全运行的状况经常发生。
技术实现要素：
本发明的目的在于克服上述不足，提供一种空冷器积灰状态感知系统与计算方法。
基于深度学习与红外图像识别的空冷器积灰状态感知系统，包括汽轮机通过汽轮机排汽与排汽分配管输入端相连；入口温度传感器布置于空冷凝汽器管束内侧，出口温度传感器布置于空冷凝汽器管束外侧，入口温度传感器和出口温度传感器输出端均经温度传输总线与数据采集模块输入端相连；入口风速传感器布置于空冷凝汽器管束内侧，出口温度传感器布置于空冷凝汽器管束外侧，入口风速传感器和出口风速传感器输出端均经风速传输总线与数据采集模块输入端相连；红外热像仪布置于空冷凝汽器管束外侧不同角度处，红外热成像仪输出端经图像数据传输总线与数据采集模块输入端相连；数据采集模块输出端经数据传输总线与工控机输入端相连；流量计输出端经输水管与凝结水箱输入端相连；配电柜输出端经电源线与数据采集模块、工控机、空冷风机输入端相连。
作为优选：所述的红外热成像仪安装于空冷凝汽器管束的迎风面和背风面。
该模型的网络输入数据是一个4维tensor，将输入的红外图像分解为2-4个特征层，尺寸为1280-60000,28-256,28-256,2-4，分别表示一批图片的个数1280-60000、图片的宽的像素点个数28-256、高的像素点个数28-256和信道个数2-4；首先使用多个卷积神经网络层进行图像的特征提取，卷积神经网络层的计算过程如下步骤：
卷积层1：卷积核大小a×a，卷积核移动步长1-2，卷积核个数a1，池化大小a2×a2，池化步长1-3，池化类型为最大池化，激活函数melu，输出a1张图像,a＝3-7,a1＝28-256,a2＝3-7；
卷积层2：卷积核大小b×b,卷积核移动步长1-2，卷积核个数b1，池化大小b2×b2，池化步长1-3，池化类型为最大池化，激活函数melu，输出b1张图像b＝3-7,b1＝28-256,b2＝3-7；
卷积层3：卷积核大小c×c,卷积核移动步长1-2，卷积核个数c1，池化大小c2×c2，池化步长1-3，池化类型为最大池化，激活函数meanout，输出c1张图像,c＝3-7,c1＝28-256,c2＝3-7；
卷积层4：卷积核大小d×d,卷积核移动步长1-2，卷积核个数d1，池化大小d2×d2，池化步长1-3，池化类型为平均池化，激活函数meanout，输出d1张图像,d＝3-7,d1＝28-256,d2＝3-7；
全连接层：隐藏层单元数512-1024，激活函数meanout；第四池化层后依次配置了2-4层全连接层；
分类层：即为全连接层的第三层，该隐藏层单元数10-20，激活函数meanout；
其中melu(x)＝δelu(x),为指数线性单元elu的变种，加入一个的系数δ，参数δ是高斯分布随机产生的参数与训练的参数结合取得的；公式如下：
meanout引入系数λ，公式如下：
n表示输入样本的个数，d表示上一层节点的个数，m表示本层节点的个数，k表示每个隐藏层节点对应了k个“中间隐层”节点，这k个“中间隐层”节点都是线性输出的，而meanout的每个节点就从这k个“中间隐层”节点中取平均值，这里结果只表示取出的第i列，下标i前的省略号表示对应第i列中的所有行；
参数初始化，所有权重矩阵使用random_normal(0.0-0.00,0.01-0.001)，所有偏置向量使用constant(0.0-0.00)；使用crossentropy作为目标函数，使用梯度下降法进行参数更新，学习率设为0.01-0.001；
空冷器积灰状态计算方法由下述方程组定量描述：
用psnr的值来评价消除抖动的效果：
利用彩色图像中的r、g、b信息将各自分量的亮度作为灰度图的灰度值：
f1(i,j)＝r(i,j)(6)
f2(i,j)＝g(i,j)(7)
f3(i,j)＝b(i,j)(8)
微分算子计算在一像素点(x，y)处的梯度的幅度为：
根据热平衡原理，汽轮机排汽放热量、acc传热量和冷却空气吸热量相等；
qe＝k×f×δt0＝g×cp×δta(11)
空冷凝汽器的对数平均温差：
δta＝ta2-ta1(13)
积灰状态总传热热阻与清洁状态总传热热阻之差为由于积灰所产生的污垢热阻，称为灰垢热阻：
通过采集空冷器积灰状态红外图像，监测空冷器壁温、空冷器出、入口冷却空气温度、冷却风流量，结合空冷器换热面积结构参数，联立式(4)～(5)对空冷器红外图像进行去抖处理，联立式(6)～(8)对得到的图像进行图像灰度化处理，联立式(9)～(10)对图像进行边缘检测处理，获得积灰图像特征信息；联立式(11)～(14)得到空冷器积灰热阻；以空冷器积灰图像为输入量，空冷器积灰热阻为目标量，建立基于深度学习网络和红外图像的空冷器积灰状态识别与分析模型。
本发明的有益效果是:
一是本发明考虑到火电站直接空冷系统排汽温度随环境气温的变化出现高频率、大浮动波动，环境高温时段排汽压力居高不下的特点，基于热平衡原理：汽轮机排汽放热量、空冷凝汽器换热量和冷却空气吸热量相等的特点，联合外部传感器通过测量模拟汽轮机排汽的水温与散热管道出口水温，翅片进出口风温等间接模拟空冷凝汽器运行状态。
二是本发明基于深度学习与红外图像处理的方法处理红外图像，获取图像数据实时变化情况，利用深度学习算法在处理图像上的优势，节省人力，提高获取空冷器运行状态数据的效率，提高问题预警的能力，增加了监测人员处理问题的时间，为空冷凝汽器的安稳、高效运行提供很大的帮助。
三是本发明考虑积灰对透光率的影响，当灰尘颗粒沉积量达到一定浓度时，图像数据r-g-b值发生变化，以预测积灰厚度。同时计算与灰垢积聚关联度较大的6个变量：积灰时间、汽轮机排汽热负荷、排汽温度、环境气温、空冷凝汽器出口风温和迎面风速。对比计算出空冷凝汽器运行状态与灰垢厚度关系。
四是本发明的结构简单、合理，造价低廉，可实施性好，更重要的是提高了功率预测精度，旨在建立更加精确的空冷凝汽器积灰与清洗预测系统，方便电厂操作人员日常检测和清洗工作，准确及时调整空冷凝汽器维护计划，为科学发电优化调管设备提供依据，实现发电收益最大化。
附图说明
图1为基于深度学习与红外图像识别的空冷器积灰状态感知系统；
图2为图像及数据处理流程框图；
图3为卷积神经网络示意图；
附图标记说明：1汽轮机，2汽轮机排汽，3排汽分配管，4空冷凝汽器管束，5空冷凝汽器翅片，6入口温度传感器，7出口温度传感器，8入口风速传感器，9出口风速传感器，10红外热像仪，11数据采集模块，12工控机，13空冷风机，14凝结水箱，15凝结泵，16温度传输总线，17风速传输总线，18图像数据传输总线，19数据传输总线，20配电柜，21电源线，22流量计。
具体实施方式
下面结合实施例对本发明做进一步描述。下述实施例的说明只是用于帮助理解本发明。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。
本发明的基于深度学习与红外图像识别的空冷器积灰状态感知系统，针对现有不能实时感知空冷系统运行状态的缺点，提出基于深度学习与红外图像识别的空冷器积灰状态感知系统，现有的深度学习感知模型研究只是对颗粒污垢的生成机理和影响因素的初步探索，还没有专门针对空冷器灰垢深入分析其影响因素与灰垢特性间的关联特性，也未建立一个普适的经过实验检验的机理模型。本发明基于气深度学习与红外图像识别的空冷器积灰状态感知系统，在深度学习算法的基础上，利用红外图像获取的信息与空冷器运行状态参数，通过深度学习算法和空冷凝汽器积灰热阻值计算方法，实时获知空冷凝汽器运行状态，实现空冷凝汽器短期精确感知、长期预报功能。
本发明的构思基础是，针对空冷凝汽器翅片积灰影响传热效率的特性，针对传统通过经验确定清洗周期的现状，1)本发明提出针对积灰对空冷凝汽器运行参数的影响为基础，通过红外热像仪获取实时红外图像，利用深度学习以及红外图像处理方法从图像中提取获得所需参数，并将其与灰垢积聚关联度较大的6个变量：积灰时间、汽轮机排汽热负荷、排汽温度、环境气温、空冷凝汽器出口风温和迎面风速结合作为基本检测对象和研究对象，推演拟合出空冷凝汽器运行状态与灰垢厚度关系。以此为基础作为空冷凝汽器运行参数采集系统和积灰状态感知系统，为其拟合积灰厚度和感知运行状态做参数支撑。2)本发明提出采用深度学习以及红外图像处理检测空冷凝汽器运行与积灰状态。通过红外成像技术实时了解空冷系统的状态，对于及时发现、处理、预防重大失误的发送可以起到非常关键而有效的作用，并可通过对空冷凝汽器的红外图像处理实时得到其积灰程度；由于北方某些地区干旱少水，工作环境恶劣，翅片管束时常出现密封度下降或冻裂等现象。这些故障现象均可通过红外成像方式远程监测出来，给平时养护工作带来极大便利，同时也可通过红外成像技术实时了解翅片管表面温度场分布，监测空冷凝汽器温差温升变化，可以实时了解空冷系统工作状态。
参照图1基于深度学习与红外图像识别的空冷器积灰状态感知系统，它包括：汽轮机1通过汽轮机排汽2与排汽分配管3输入端相连；入口温度传感器6布置于空冷凝汽器管束4内侧，出口温度传感器7布置于空冷凝汽器管束4外侧，入口温度传感器6和出口温度传感器7输出端均经温度传输总线16与数据采集模块11输入端相连；入口风速传感器8布置于空冷凝汽器管束4内侧，出口温度传感器9布置于空冷凝汽器管束4外侧，入口风速传感器8和出口风速传感器9输出端均经风速传输总线17与数据采集模块11输入端相连；红外热像仪10布置于空冷凝汽器管束4外侧不同角度处，红外热像仪10输出端经图像数据传输总线18与数据采集模块11输入端相连；数据采集模块11输出端经数据传输总线19与工控机12输入端相连；流量计22输出端经输水管与凝结水箱14输入端相连；配电柜20输出端经电源线21与数据采集模块11、工控机12、空冷风机13输入端相连。
所述配电柜20是现场所有仪器的供电核心，主要包括配电柜箱体、空开、交流接触器以及相关保险等防火放电电路。其中总空开可承受大功率电路，总空开分流给锅炉供电电路、给水泵电路、空冷风机电路以及控制器电路。实现各个部分器件随意开关，并采用摁键和指示灯通过交流接触器方便日常操作。
所述数据采集系统主体框架由红外热像仪10、入口温度传感器6，出口温度传感器7、入口风速传感器8，出口风速传感器9构成。采用数个温度传感器分别对空冷凝汽器机组的水箱水温、环境气温、入水口水温、出水口水温、入口风温、出口风风温进行测量。将温度传感器串行连接到控制系统，可仅用数根线实现数据的传输。
所述风速传感器5精度强分辨率高，实现对空冷凝汽器管道翅片上的入口出口风速测量，与控制系统通过的a、b信号线连接，使用modbusrtu协议进行通信。
所述的流量计采用高精度电磁流量计。
所述控制系统是采集现场数据的核心器件，使用主控芯片，液晶显示器作为屏幕显示监测数据，控制和采集温度传感器、风速传感器以及流量计的正常运作和信号输出。通过串口报文传输的方式，将数据传送给上位机数据处理系统。由于实验室空间距离较短，可用串口电平传输数据，控制器同时配备信号传输线，通过串行总线发送给串口转换模块，最后串口转换模块通过串行总线将空冷凝汽器积灰特性的运行参数发送给上位机，实现远距离传输。
所述上位机11数据处理系统使用windows平台labview虚拟仪器软件实现，可以识别串口转换模块作为输入信号源，然后在组态软件上设置输入口的对应的参数类别，即可在上位机中显示、存储空冷器的运行状态参数，并调用matlab中算法实现预测灰垢与计算换热系数。可由数据采集系统传输的数据进行计算换热面积灰程度的特征参数。并根据热平衡原理，模拟汽轮机排汽放热量、空冷凝汽器换热和冷却空气吸热量相等计算翅片换热量。即可在上位机中显示、存储空冷凝汽器积灰特性的运行参数，并基于这类参数实时预测空冷凝汽器积灰程度与最优清洗周期，方便电厂操作人员日常检测和清洗工作，准确及时调整空冷凝汽器维护计划。上位机数据处理系统通过labview虚拟仪器软件设计图像显示、图像处理界面。
参照图2和图3所述基于深度学习网络的图像检测方法，其特征在于，包括：它包含了若干卷积层、池化层和全连接层。将数据采集模块采集到的红外图像建立成一个数据库，并使用卷积神经网络算法对数据库中大部分的图像数据进行训练，用剩下小部分的数据进行验证。搭建起合适的模型，以便于下一步的数据图像处理。具体步骤如下：输入红外图像；对所述红外图像进行处理，生成多个训练样本；利用所述多个训练样本进行基于卷积神经网络的图像识别训练，生成训练模型；以及使用所述训练模型检测上传的图像中是否含有所述红外图像。
所述深度学习网络具体为卷积神经网络算法，该模型的网络输入数据是一个4维tensor，将输入的红外图像分解为2-4个特征层，尺寸为(1280-60000,28-256,28-256,2-4)，分别表示一批图片的个数1280-60000(图片的个数不宜过少，也不宜过多)、图片的宽的像素点个数28-256、高的像素点个数28-256和信道个数2-4。首先使用多个卷积神经网络层进行图像的特征提取，卷积神经网络层的计算过程如下步骤：
卷积层1：卷积核大小a×a，卷积核移动步长1-2，卷积核个数a1，池化大小a2×a2，池化步长1-3，池化类型为最大池化，激活函数melu，输出a1张图像,a＝3-7,a1＝28-256,a2＝3-7。
卷积层2：卷积核大小b×b,卷积核移动步长1-2，卷积核个数b1，池化大小b2×b2，池化步长1-3，池化类型为最大池化，激活函数melu，输出b1张图像b＝3-7,b1＝28-256,b2＝3-7。
卷积层3：卷积核大小c×c,卷积核移动步长1-2，卷积核个数c1，池化大小c2×c2，池化步长1-3，池化类型为最大池化，激活函数meanout，输出c1张图像,c＝3-7,c1＝28-256,c2＝3-7。
卷积层4：卷积核大小d×d,卷积核移动步长1-2，卷积核个数d1，池化大小d2×d2，池化步长1-3，池化类型为平均池化，激活函数meanout，输出d1张图像,d＝3-7,d1＝28-256,d2＝3-7。
全连接层：隐藏层单元数512-1024，激活函数meanout。第四池化层后依次配置了2-4层全连接层，以3层为例，神经元个数分别为1024、512、10，将第四池化层输出图像依次转化为对应的一维向量，向量元素个数分别为1024、512、10，第三全连接层输出10个元素的一维数组对应本实验设计的10组图像分类。
分类层：即为全连接层的第三层，该隐藏层单元数10-20，激活函数meanout。
其中melu(x)＝δelu(x),为指数线性单元elu的变种，加入一个的系数δ，参数δ是高斯分布随机产生的参数与训练的参数结合取得的。公式如下：
meanout不同于maxout的地方在于，maxout引入的参数过多，而meanout则在maxout的基础上减少了引入的参数规模，但不影响训练的准确度,所以meanout引入了另外的一个系数λ，λ需要通过训练和经验进行调整，以提高模型的效率。公式如下：
n表示输入样本的个数，d表示上一层节点的个数，m表示本层节点的个数，k表示每个隐藏层节点对应了k个“中间隐层”节点，这k个“中间隐层”节点都是线性输出的，而meanout的每个节点就从这k个“中间隐层”节点中取平均值，这里结果只表示取出的第i列，下标i前的省略号表示对应第i列中的所有行。
参数初始化，所有权重矩阵使用random_normal(0.0-0.00,0.01-0.001)，所有偏置向量使用constant(0.0-0.00)。使用crossentropy作为目标函数，使用梯度下降法进行参数更新，学习率设为0.01-0.001。
参照图2本发明的一种基于深度学习与红外图像识别的空冷器积灰状态计算方法，其特征在于：利用红外图像获取的信息与空冷器运行状态参数，通过深度学习算法和空冷凝汽器积灰热阻值计算方法，实时获知空冷凝汽器运行状态，实现空冷凝汽器短期精确感知、长期预报功能，整个计算方法由下述方程组定量描述：
用psnr的值来评价消除抖动的效果，其中psnr定义是：
利用彩色图像中的r、g、b信息将各自分量的亮度作为灰度图的灰度值，公式如下所示：
f1(i,j)＝r(i,j)(6)
f2(i,j)＝g(i,j)(7)
f3(i,j)＝b(i,j)(8)
微分算子计算在一像素点(x，y)处的梯度的幅度为：
根据热平衡原理，汽轮机排汽放热量、acc传热量和冷却空气吸热量相等。
qe＝k×f×δt0＝g×cp×δta(11)
空冷凝汽器的对数平均温差：
δta＝ta2-ta1(13)
积灰状态总传热热阻与清洁状态总传热热阻之差为由于积灰所产生的污垢热阻，称为灰垢热阻：
其中，其中fi、fi-1分别表示第i和第(i+1)帧的像素值，fk(i,j)(k＝1,2,3)表示在坐标(i,j)点各分量的灰度值，