Infrared Temperature Measurement Calibration and Temperature Field Reconstruction of Aero-Engine Blade Surface
-
摘要:
随着我国航空业的蓬勃发展，航空发动机技术取得了长足的进步。作为航空发动机的关键组件，航空发动机叶片的研制关乎全局。对其温度进行实时、精确的测量监控有助于推动航空发动机叶片相关技术的突破。本文提出了一种基于红外测温技术的航空发动机叶片温度校准算法和温度场三维重构策略。该工作利用多层感知机网络的函数拟合能力，建立了对目标点精确温度与其周围温度分布的函数关系，并将红外测温的误差控制在1.24℃。在此基础上，采用投影法和三维点云法向量估计的创新策略，实现了二维红外图像到三维空间位置映射。该方法能从单视角红外图像得到温度场三维重构，重构过程简单，摆脱了对多视角图像的依赖。
Abstract:With the vigorous development of the aviation industry in China, significant progress has been made in aeroengine technology. As key components of aeroengines, the development of aeroengine blades is crucial. Real-time and accurate monitoring of the blade temperature will help promote breakthroughs in related technologies for aeroengine blades. This study proposes a temperature calibration algorithm and three-dimensional reconstruction strategy for the temperature field based on infrared temperature measurements for aeroengine blade temperature monitoring. Based on the function fitting ability of the multilayer perceptron network, the functional relationship between the precise temperature of the target point and its surrounding temperature distribution is established, and the error of infrared temperature measurement is controlled within 1.24 ℃. On this basis, through an innovative projection method and 3D point cloud normal vector estimation, the position mapping from 2D infrared image to 3D space has been achieved. We achieved a three-dimensional temperature field reconstruction via single-view infrared images through a simple process without dependence on multiview images.
-
0. 引言
在化工、冶金、火电等工业项目以及一些制冷相关系统的民用项目中，都存在着需要对项目内部温度进行控制的环节，以保证项目内部的温度符合工业或民用需求。此外，为了节约水资源，项目中常采用循环水冷却系统作为制冷系统。这种循环冷却水系统在设计时于系统中换热器、管路、冷却塔布水器以及水泵处都留有2%~5%的压力余能（水头）[1-2]，系统中这些压力余能累积就有相当可观的剩余能量。因此，在此理论基础上，多位工程师成功研制出纯水动力冷却塔以达到节约能源的目的。
然而，由于水动力冷却塔是通过循环水系统中的余压进行工作所以相较于电动力冷却塔而言工作情况更复杂，而且，运行过程中水轮机整机性能也会一直受到富余水头变化情况的影响。因此为了避免由水循环系统中富余水头变化带来的冷却塔制冷功率变化，可在纯水动力冷却塔的基础上进行结构优化，设计出水电混合动力冷却塔，同时基于热学成像原理设计出适用于水电混合动力冷却塔的温控系统，以保证在任何情况下冷却塔都能有效地完成工作。在设计冷却塔时，需构建冷却塔热力计算模型进行热力计算。本文将根据逆流式冷却塔热交换的基本原理以及水电混合动力冷却塔的结构特性，基于Merkel方程构建数学模型进行热力计算。同时，为了能够正确地识别冷却塔所在的瞬时工况并判断其是否需要接入电机以进行动力补偿工作，使冷却塔能够完成冷却目标，本文也基于冷却塔的工作环境进行了用于冷却塔的红外热像温控系统设计。并通过设计的温控系统进行实验对所建立热学计算模型进行验证。
1. 冷却塔热力计算模型
逆流式水电混合动力冷却塔内，温度较高的循环水在经由进水管导入水轮机做功后由尾水管导出至洒水机构，并经过喷头由上至下经过喷洒区、填料区以及雨区[3]。其中，在喷洒区和雨区中，水以水滴的形式通过，在中间的薄膜式填料区内，水以水流的形式通过。冷却塔内空气的方向则与水流方向相反，自下而上，依次在雨区、填料区和喷洒区内部与水进行热质交换，整个热力交换过程如图 1所示。
![]()
图 1 冷却塔热力交换示意图Figure 1. Heat exchange of cooling tower schematic diagram1 Motor; 2 Fan; 3 Water turbine; 4 Sprayer; 5 Packing为了便于计算，将图 1中的冷却塔内部热力交换过程进行简化，建立冷却塔热力计算模型示意图，如图 2所示。
![]()
图 2 冷却塔热力计算模型示意图Figure 2. Thermal calculation model of cooling tower schematic diagram如图 2所示，循环水系统中冷却水自上而下运动，进口和出口处的水温和水质量流率分别为tw、tw－dtw和mw、mw－dmw，空气自下而上运动，空气焓、空气温度、空气含湿量以及空气质量流率在进口和出口处分别为h、h+dh，ta、ta+dta，φ、φ+dφ以及ma、ma+d ma。
由此根据质量守恒定律，可基于循环水量的变化等于空气含湿量的变化的原理建立空气含湿量的变化方程：
$$ {\rm{d}}\varphi = {\rm{d}}{m_{\rm{w}}} = {\beta _{{\rm{xv}}}}\left( {{\varphi ^{\prime \prime }} - \varphi } \right){\rm{d}}V $$ (1) 式中：φ为初始空气含湿量；φ"为饱和空气含湿量，kg[水]/kg[空气]；βxv为填料容积散质系数，kg[水]/m3；dV为填料微元体积，m3；mw为水质量流率，kg/s。
在冷却塔内热质交换过程中，空气温度升高，湿度变大，依靠风机旋转产生的压力差驱动源源不断地流入塔内。由于冷却塔内水和空气之间相对温差较小，因此可将辐射传热忽略，两者之间只考虑接触传热和蒸发传热的存在，其传热过程主要发生于填料区，故可从填料体积微元dV内的质量和能量守恒方程建立冷却塔的基本热力计算理论模型[4]，即Merkel模型。
在冷却塔填料体积微元dV内，热量交换方程可分为由水汽化传递的热量方程(2)和循环水降温传递的热量方程(3)两部分：
$${\rm{d}}{h_1} = \gamma {\beta _{{\rm{xv}}}}\left( {\varphi '' - \varphi } \right){\rm{d}}V$$ (2) $${\rm{d}}{h_2} = \alpha \left( {t - \theta } \right){\rm{d}}V$$ (3) 式中：γ为水的汽化热，kJ/kg[水]；v为水的散热系数，kJ/(m3·℃)；t为当前水温，℃；θ为当前湿空气温度，℃。
综上式(2)、式(3)，可得出冷却塔内部交换的总热量为：
$${\rm{d}}h = \gamma {\beta _{{\rm{xv}}}}\left( {\varphi '' - \varphi } \right){\rm{d}}V + \alpha \left( {t - \theta } \right){\rm{d}}V$$ (4) 同时，根据质量与能量守恒定律可得循环水在热交换过程中散发的总热量：
$${\rm{d}}h{\rm{ = }}C{p_{\rm{w}}}{m_{\rm{w}}}{\rm{d}}t = \,\gamma {\beta _{{\rm{xv}}}}\left( {\varphi '' - \varphi } \right){\rm{d}}V + \alpha \left( {t - \theta } \right){\rm{d}}V$$ (5) 式中：Cpw为水的比热容，kJ/(kg·℃)；dt为填料体积微元dV中循环水的温降，℃。
式(5)为逆流式冷却塔内部热交换的基本方程。
在式(5)的基础上，根据Merkel方程热力计算模型[5]引入焓概念的理论基础进行简化，获得Merkel数（冷却数）Me：
$${\rm{Me}} = \frac{{{\beta _{{\rm{xv}}}}V}}{{{m_{\rm{w}}}}} = \int_{{t_2}}^{{t_1}} {\frac{{C{p_{\rm{w}}}{\rm{d}}t}}{{h'' - h}}} $$ (6) 式中：t1、t2分别表示冷却塔内循环水的进塔和出塔水温；h"为水温t℃时饱和湿空气的比焓；h为湿空气温度θ℃时饱和湿空气的比焓。
最终，根据上文建立的冷却塔热力计算模型可获得冷却塔出口空气温度θ2：
$${\theta _2}{\rm{ = }}{t_1} - \frac{{\int_{{t_2}}^{{t_1}} {{{\rm{e}}^{{\rm{Me}}/\lambda }}} {\rm{d}}t + {t_2} - {\theta _1}}}{{{{\rm{e}}^{{\rm{Me}}/\lambda }}}}$$ (7) 式中：θ1为进口空气温度，℃；λ为气水比。
2. 运用红外热成像技术建立冷却塔温控系统
2.1 冷却塔温控系统的构建
红外热成像技术是将不可见的红外辐射转化为可见图像的技术，利用这一技术研制的装置统称为红外热成像装置或红外热像仪。红外热像仪的主要特点为：①为非接触式测量，能够检测运动目标、微小目标及带电目标温度[6]；②测温范围广：高于－273.15℃绝对零度的物体都会向外产生红外辐射能，所以理论上讲，红外辐射测温是没有上限的，具有宽广的测温范围[7]；③温度分辨率高，可准确区分较小温差，获取精确温度值；④数据输出方式主要为灰度或彩色云图，目标各点温度值读取方便；⑤可进行数据存储和计算机处理，便于自动化设计。
根据上文所建冷却塔热力计算模型，设计出如图 3所示的逆流式冷却塔热力性能监测模型，将实时测量获得的环境参数、运行参数、电参数代入监测计算模型，求解方程结果，实现冷却塔工作状态的动态测量。
![]()
图 3 逆流式冷却塔热力性能监测模型Figure 3. Thermal performance monitoring model of countercurrent cooling tower在设计冷却塔温控系统时，可利用红外热成像技术进行冷却塔出口空气温度的采集，并将所获取数据反馈至控制器中进行分析，判断冷却塔所需风量是否满足需求，最终通过控制器对风机进行转速调节，形成温控系统闭环，温控系统示意图如图 4所示。
![]()
图 4 冷却塔温控系统示意图Figure 4. Schematic diagram of temperature control system of cooling tower1 The cooling tower; 2 Thermal infrared imager; 3 The controller; 4 Environmental parameter sensor2.2 冷却塔温控系统温控精度误差分析
红外热成像仪通过光学系统采集被测物体产生的红外辐射，经传感器将光信号转换为电信号，再通过放大电路及补偿电路处理后，并最终由系统终端输出结果，其系统结构如图 5所示。
由于红外热成像仪为非接触式测量仪器，其测量精度较低，根据其技术原理及系统结构可以推断出，红外热成像仪的测量误差主要为光学系统造成的误差。光学系统误差产生的因素较复杂：①发射率，被测物体的材料形状、表面粗糙度及凹凸度均可影响物体的发射率，而红外热成像仪能从物体上接受的辐射能量大小与其发射率成正比，因此，当红外热成像仪对形状怪异以及表面质量较差的物体进行测温时测量精度较低，测量误差较大。②被测物体尺寸，当被测物体小于等于红外热成像仪的测温视场时，仪器会受到测量区域以外的环境影响，从而产生测量误差。③被测物体与红外成像仪之间的距离，当热成像仪与被测物体之间的距离过大时，热成像仪的分辨率越低，测量精度越差。④环境因素，由于红外热成像仪时非接触式测量，因此被测物体与热成像仪之间存在其它物质，这些物质不仅会和被测物体产生热交换，也会吸收被测物体产生的红外线辐射能量并将其转变成其它形式的能量，从而产生测量误差。
综上因素所产生的红外热成像仪中光学系统的测量误差最终会经过传感器变为电信号后，并通过热成像仪内部的信号处理电路进行信号处理后由终端输出至冷却塔温控系统中的控制器中影响整个温控系统的温控精度。
由图 6可以得知，由红外热成像仪的终端输出的出口空气温度是冷却塔温控系统运行的重要依据之一，因此红外热成像仪的测量精度误差是影响冷却塔温控系统精度的重要原因。相较之下，温控系统中的其它传感器均为接触式传感器，测量精度高，对温控系统的温控精度误差影响小。最终，根据红外热成像仪测量精度F1、出口水温温度传感器测量精度F2、环境参数传感器测量精度F3以及安全系数S设定冷却塔温控系统温控精度F：
$$ F = \left\{ {1 - \left[ {\left( {1 - {F_1}} \right) \times \left( {1 - {F_2}} \right) \times \left( {1 - {F_3}} \right)} \right]} \right\} \times S \times 100\% $$ (8) 3. 利用手持式红外测温仪进行实验与分析
为了进一步验证冷却塔热力计算模型的准确性，本应该利用红外热成像技术对冷却塔工作过程进行温度监控，根据红外特性判别冷却塔出口处空气温度，如图 7所示。
![]()
图 7 红外热像温控系统监测图Figure 7. Monitoring diagram of infrared imaging temperature control system然而由于红外温控系统结构复杂，建造周期长，故将实验中的探测过程简化，特采用手持式红外测温仪进行真机的现场实验。
1）实验装置
实验装置如图 8所示。
2）实验装置参数
冷却塔参数：流量Q＝100 t/h，工作水头H＝10 m，塔高h＝3.1 m，水轮机效率η＝0.85，风机额定转速n＝360 r/min，风机额定功率P＝5 kW，风机叶片直径D＝1.8 m。
手持式红外测温仪参数：测温范围：－32℃~330℃，测量精确度±2%，分辨率：0.1℃。
3）实验步骤
① 通过温度计获得实验当地空气干球温度29.6℃，湿球温度23.3℃，相对湿度70.9%。
② 预设冷却塔进塔水温，根据冷却塔热力计算模型以及文献[8]获得理论出塔空气温度。
③ 根据预设进塔水温进行真机实验，获得实际出塔空气温度并将其与理论出塔空气温度进行对比分析。
④ 现场实验数据及分析
现场实验数据与理论数据对比分析如图 9所示。
从图中可以明显地看出现场实验数据与理论计算模型得出的数据之间存在着一一对应关系，成功验证了上文所建立的冷却塔热力计算模型的正确性，同时也能证明手持式红外测温仪的参数能够满足实验要求。
4. 结论
根据水电混合动力冷却塔的结构特性，基于热质平衡方程可成功建立适用于冷却塔温控系统的热力计算模型。由于实验条件难以实现，故将实验过程简化，使用手持式红外测温仪进行真机的现场实验。虽然手持式红外测温仪的分辨率较低，实验结果与理论数据之间有所偏差，但从获得现场实验数据和理论计算数据的对比图中可以看出，数据之间存在对应关系，可成功验证冷却塔热力计算模型的正确性以及该模型应用于冷却塔温控系统设计的可行性。未来实验条件允许的情况下，可采用红外热像技术进行真机实验，在提高设备分辨率的前提下获得更加准确的实验数据进行实验论证。
-
表 1 样本平均温度误差E/mm
Table 1 Average temperature error of test samples
Test sample Average temperature error of sample/℃ Average value/℃ 1 1.0770 1.24 2 1.8330 3 0.9736 4 0.9483 5 1.3708 
下载: 导出CSV
-
[1] 胡承波. 中国航空制造业企业技术创新长效机制研究[D]. 武汉: 武汉理工大学, 2011. HU Chengbo. The Research on the Long-term Mechanism of Enterprise Technological Innovation in Aviation Manufacturing of China[D]. Wuhan: Wuhan University of Technology, 2011.
[2] 王魁汉, 李友, 王柏忠. 温度测量技术的最新动态及特殊与实用测温技术[J]. 自动化仪表, 2001(8): 3-9. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDYB200108000.htm WANG Kuihan, LI You, WANG Baizhong. The latest development status of the technology of the temperature measurement and the unique and practical measurement technology[J]. Process Automation Instrumentation, 2001(8): 3-9. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDYB200108000.htm
[3] YU Y B, CHOW W K. Review on an advanced high-temperature measurement technology: the optical fiber thermometry[J]. Journal of Thermodynamics, 2009(1): 823482.
[4] 张志学, 刘忠奎, 张玉新, 等. 航空发动机壁温测量方法综述[C]// 2015航空试验测试技术学术交流会论文集, 2015: 139-142. ZHANG Zhixue, LIU Zhongkui, ZHANG Yuxin, et al. Survey on the wall temperature measurement of aeroengine[C]//Proceedings of the 2015 Aviation Test Technology Academic Exchange Conference, 2015: 139-142.
[5] 唐文彬. 某型涡桨发动机燃烧室温度场热电偶动态响应规律的技术研究[D]. 南京: 南京航空航天大学, 2013. TANG Wenbin. Research on Dynamic Response Rule of Thermocouple for One Type of Turboprop Engines Combustor[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2013.
[6] 黄名海, 臧树升, 葛冰, 等. 热风洞中涡轮叶片温度场红外热像测量方法[J]. 航空动力学报, 2014, 29(11): 2679-2683. https://www.cnki.com.cn/Article/CJFDTOTAL-HKDI201411019.htm HUANG Minghai, ZANG Shusheng, GE Bing, et al. Method of infrared thermography measurement for temperature field of turbine vane in hot wind tunnel[J]. Journal of Aerospace Power, 2014, 29(11): 2679-2683. https://www.cnki.com.cn/Article/CJFDTOTAL-HKDI201411019.htm
[7] CHRISTENSEN L, CELESTINA R, SPERLING S, et al. Infrared temperature measurements of the blade tip for a turbine operating at corrected engine conditions[J]. Turbomach, 2021, 143(10): 101005. DOI: 10.1115/1.4050675
[8] 肖冬杰, 李其申, 周翠岭. 基于稀疏表示的自适应图像融合方法研究[J]. 计算机应用与软件, 2014, 31(3): 203-206. https://www.cnki.com.cn/Article/CJFDTOTAL-JYRJ201403054.htm XIAO Dongjie, LI Qishen, ZHOU Cuiling. Study on sparse representation based adaptive image fusion methods[J]. Computer Applications and Software, 2014, 31(3): 203-206. https://www.cnki.com.cn/Article/CJFDTOTAL-JYRJ201403054.htm
[9] 刘兆栋. 基于稀疏表示理论的图像去噪与融合算法研究[D]. 重庆: 重庆大学, 2016. LIU Zhaodong. Research on Image Denoising and Fusion Based on Sparse Representation[D]. Chongqing: Chongqing University, 2016.
[10] 胡建文. 基于多尺度滤波和稀疏表示的图像融合方法研究[D]. 长沙: 湖南大学, 2013. HU Jianwen. The Research on Image Fusion Based on Multiscale Filter and Sparse Representation[D]. Changsha: Hunan University, 2013.
[11] 肖传民, 亓琳, 史泽林. 多尺度双边滤波及其在图像分割中的应用[J]. 信息与控制, 2009, 38(2): 229-233. DOI: 10.3969/j.issn.1002-0411.2009.02.018 XIAO Chuanmin, QI Lin, SHI Zelin. Multi-scale bilateral filtering with application to image segmentation[J]. Information and Control, 2009, 38(2): 229-233. DOI: 10.3969/j.issn.1002-0411.2009.02.018
[12] 王华伟. 基于红外热成像的温度场测量关键技术研究[D]. 西安: 中国科学院研究生院(西安光学精密机械研究所), 2013. WANG Huawei. Research on the Key Technologies of Temperature Field Measurement Based on Thermal Infrared Imager[D]. Xi'an: Xi'an Institute of Optics & Precision Mechnics, Chinese Academy of Sciences, 2013.
[13] 李云红. 基于红外热像仪的温度测量技术及其应用研究[D]. 哈尔滨: 哈尔滨工业大学, 2010. LI Yunhong. Research on Temperature Measurement Technology and Application Based on Infrared Thermal Imager[D]. Harbin: Harbin Institute of Technology, 2010.
[14] 程胜. 人体三维远红外成像及其测温的研究[D]. 上海: 上海交通大学, 2009. CHENG Sheng. Three-Dimensional Far-infrared Thermal Imaging and Its Temperature Measurement for Medical Applications[D]. Shanghai: Shanghai Jiao Tong University, 2009.
[15] 李勇. 基于红外测温的炉内温度场重构方法初步研究[D]. 北京: 华北电力大学, 2010. LI Yong. Research on Reconstruction Algorithm of Temperature Field Based on Infrared Radiation Measurement in Boiler[D]. Beijing: North China Electric Power University, 2010.
-
期刊类型引用(1)
1. 李裴，刘慧，李满堂，朱莉，王换换. 基于STC单片机的智能化病房温控系统设计与实现. 计算机测量与控制. 2022(03): 126-132 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 71
- HTML全文浏览量: 10
- PDF下载量: 24
- 被引次数: 1
