胡长德,张海波,王朝晖,王林旭,李子杨,刘婕
(1.航天工程大学 士官学校,北京 102200;2.航天工程大学 职业教育中心,北京 101416)
直至20 世纪70 年代中期,反导防御系统一直采用雷达探测识别目标。随着反雷达识别技术不断提高,采用雷达从大量干扰物中识别出弹头变得愈来愈困难。而因为弹头、弹体、诱饵等目标在高速飞行中的温度变化及其热辐射有着很大差别[1-3],运用红外技术进行目标识别就具有一定优势。在高度深冷空间,轻诱饵由于本身热容量小,其表面温度将下降很快,而弹头的温度变化却很小,因此可以利用彼此“冷热”的不同来区分真假目标。目前,红外辐射特性测量设备越来越多地被运用到了导弹突防性能鉴定试验中。这些设备的操作维护一般由士官来完成。因此在士官院校的日常教学训练中,就必须配备红外辐射特性测量设备,用以学习相应的结构原理,训练相应的操作技能。而现在用于武器鉴定试验的红外辐射特性测量设备往往价格昂贵,耗资均在千万元以上,在士官院校很难实现按需配备。且用于执行任务的红外辐射特性测量设备的集成化程度高,操作与任务结合紧密,一些核心器件本身的指标,如红外相机的响应率、迟滞、重复性等,无法在现有测量设备中检测获得。而这些指标参数对于学员更好地学习理解设备组成原理,掌握操作技能,进而获取更精确的任务数据具有重要作用[4-6]。因此,研制价格相对低廉、具备实践教学功能的教学系统就显得尤为必要。
本研究提出了一种红外辐射特性测量教学系统,采用非制冷普通红外相机、黑体以及红外辐射特性测量软件。该系统价格不到20 万元/套,仅为实装设备售价的百分之一。除了具备实际设备所具有的辐射特性测量、标定及非均匀性校正等操作科目外,所提出的教学系统还增加了诸如响应率、迟滞、重复性等红外辐射特性测量系统中红外相机器件相关指标测试的实践科目。操作科目总计近20 个,对于该系统的结构、原理学习大有帮助。
结合红外辐射特性测量的设备任务及岗位要求,除设备和岗位的操作测试外,本教学系统增加了红外探测器这一核心器件的指标探测,同时增加典型材料红外辐射特性数据库,以帮助进行目标分析与识别。
红外辐射测量系统的组成如图1 所示,由非制冷型红外探测器、光学镜头、黑体和红外辐射测量软件组成。该系统可以根据获取的红外图像数据,反演出目标的红外辐射亮度和亮温。首先由光学镜头收集来自目标的红外辐射,经红外探测器转换为电信号,再经红外探测器接口电路送至数字处理板进行处理,然后通过主机显示并由辐射测量处理软件完成相应处理。具体实物图见图2。
图1 红外辐射特性测量系统原理图Fig.1 Principle of the infrared radiation characteristic measurement system
图2 红外辐射特性测量系统实物图Fig.2 Photo of the infrared radiation characteristic measurement system
辐射测量处理软件由计算机和对测量图像数据进行处理分析的各种软件组成,具有目标辐射功率、辐射强度、辐射出射度、红外探测器性能等指标的测试功能。辐射测量处理软件均基于Windows 平台,运用VC 语言编程。在软件设计的过程中重点考虑了可重用性和可扩展性。可重用性是指所设计的代码不作修改或稍加改动就可多次重复使用;可扩展性是指随着红外辐射测量技术的不断发展,其功能可不断加以完善。
从功能划分,红外辐射特性测量系统软件主要包括辐射测量控制软件模块和数据综合处理软件模块。辐射测量控制软件模块:完成辐射测量系统的各项控制管理功能;提供标定过程中所需的软件支撑;定义红外辐射标定单元的软件接口;将标定数据结果供其他模块使用。数据综合处理软件模块:精确计算各项辐射测量、反演及材料分析识别结果,并形成分析报告。
由图3 可知,系统软件的工作流程如下:
图3 红外辐射特性测量系统软件逻辑组成框图Fig.3 Diagram of the infrared radiation characteristic measurement system with interpolation
1)辐射测量控制软件模块
辐射测量控制模块主要包括系统设置、系统自检、定标和测量控制、数据显示、数据记录和回放等功能,并在用户信息窗口显示系统当前工作状态(正常或异常)、自检结果,以及定标和测量控制、数据接收等功能模块的工作信息。
红外辐射测量数据显示功能主要完成目标红外图像、温度分布、波段辐射亮度分布等测量数据的处理和分析。实时显示计算的X轴和Y轴的温度曲线;实时计算和显示最高温度、最低温度、平均温度、任意点温度值及其坐标;计算矩形区域内/外温度的最高值、最低值、平均值;三维方式显示图像;放大图像中任意区域。红外探测器指标测试包括响应率测量、迟滞测量、重复性测量、非均匀性测量及校正等。
2)数据综合处理软件模块
数据综合处理模块主要对红外探测器和辐射标定采集的图像和数据进行预处理、综合计算分析,反演计算目标和背景在当时气象条件下的真实红外辐射特性(温度分布、温度变化率测量),并将结果信息存储,建立靶场数据库,并根据数据库中已有的材料辐射特性进行目标材料的分析识别。
红外跟踪与测量分系统的主要功能是对目标的红外辐射量进行测量,因此就需要对系统的辐射响应进行精确标定[5]。其具体标定方法如下:
在实验室采用均匀辐射圆形孔状的标准黑体源作为均匀红外辐射目标,该目标通过辐射计光学系统成像在辐射计探测器红外焦平面上。平行光管是无挡光的离轴光学结构,其口径大于辐射计望远镜口径,黑体孔的像大于红外焦平面一个单元的面积。标定时采用中心点标定。
在系统的工作动态范围中取一系列不同的黑体辐射温度作为红外辐射标定点。直接改变黑体辐射温度,可从被测系统的辐射响应信号输出中得到一系列测量值,该响应值经模数转换后用灰度级表示。由此可得到被测系统辐射响应值同黑体目标辐射量的对应曲线[7-8]。
根据普朗克黑体辐射定律,黑体在 λ1~λ2波段中的辐射发射量为:
式中:C1=3.741 5×10-12W·cm2;C2=1.438 8 cm·K。
对充满红外焦平面一个像元的黑体目标,入射在该敏感单元上的辐射功率为:
式中:τ0为光学系统透过率;A0为望远镜通光面积;ω为系统瞬时视场。
辐射计焦平面单元尺寸为 12 μm×12 μm,光学系统焦距为48.3cm,瞬时视场立体角为ω=4μ2=6.86×10-3Sr。那么黑体在系统入射孔径上的辐射照度H与黑体辐射发射量W的关系为:
定标时,将黑体设定在不同的温度T,计算出8~14 μm 波段的黑体辐射发射量W,进一步得到系统入射孔径上的辐射照度值。当系统输出信噪比为1 时即为系统的噪声等效辐射照度(即系统的灵敏度)。辐射计用H进行标定,得到灰度级与辐射照度的曲线。
各灰度等级G对应的目标辐射照度H可按下式计算:
式中:G0为显示灰级值(满度为256 灰级);W为黑体辐射发射量;H为辐射照度。
当得到辐射计输出的灰度级后,先根据目标所在区域和波段进行非均匀性校正,而后通过定标曲线得到照度值H,即可计算出实际目标的辐射强度[9-12]。
为了更好地理解辐射特性测量设备的工作原理,掌握辐射特性测量设备的操作技能,在研究和设计该系统时增加了红外探测器指标测试操作科目。这些指标参数可以通过手动测试计算,也可以通过软件自动计算获得。下面主要阐述响应率、迟滞这两个指标测试操作科目的设计实现。
1)红外探测器响应率测量
红外探测器的响应率表征的是探测器对于目标红外辐射的响应程度。即探测器电输出信号的基本量与输入功率基本量的比值,单位是V/W。
式中:Uλ为来自目标的均值波长为λ的红外辐射经探测器转换输出的电压值;Pλ为来自目标的均值波长为λ的红外辐射功率。对于已知温度的黑体,其辐射功率也是已知的。探测器的口径也是一定的,比黑体要小,在进行实际操作时,红外探测器与黑体的距离很近,可以忽略掉大气衰减,实际进入探测器能量与红外探测器的口径成正比[13-15],假设探测器口径和黑体面积的比值为α,Uλ为红外探测器实际输出。
对于该操作科目,检测时设置一系列黑体温度(N个),分别为T1、T2、T3、···、TN,分别记录对应的红外辐射功率Pλ1、Pλ2、Pλ3、···、PλN,红外探测器输出电压值Uλ1、Uλ2、Uλ3、···、UλN。可以得到响应率为:
2)红外探测器迟滞测试
迟滞是指传感器在正(输入量增大)反(输入量减小)行程中输出不一致的程度。该指标说明,红外探测器在测量目标时,对于升温目标与降温目标,同样的温度变化,其输出结果不一致的程度[16-17]。
对于该科目,检测时设置一系列正(升温)行程黑体温度(N个),分别为T1、T2、T3、···、TN,分别记录对应红外探测器输出灰度值H1、H2、H3、···、HN。然后开始降温,设置同样的反(降温)行程黑体温度(N个),分别为TN、TN-1、TN-2、···、T1,分别记录对应的红外探测器输出灰度值不难看出,最大灰度值为HN,即。则迟滞 γ为:
式中:i为自然数,1≤i≤N;为正反行程中同温度值时的最大差值。
红外辐射特性测量系统实验如图4 所示,由黑体辐射源、光学镜头、红外相机、辐射测量软件组成。根据任务及教学训练要求,主要阐述分析标定、目标红外辐射特性、红外探测器迟滞实验结果。
图4 红外辐射特性测量系统实验组成Fig.4 Experimental composition of the infrared radiation characteristic measurement system
选择自动定标,自动定标参数起始温度设置为30 ℃,截止温度设置为200 ℃,步长为5 ℃。重复定标5 次,获得定标系数如表1 所示。
表1 定标数据Tab.1 Calibration data
由表1 数据可得,该系统的红外辐射响应和温度变化基本呈线性关系,在30~200 ℃的范围内的响应系数为136.4。
完成标定后,在红外辐射特性测量反演软件中导入标定实验的实验结果,系数K和偏移量B,对实测目标红外辐射特性进行测量,获取数据如表2 所示。
表2 目标红外辐射特性实验数据Tab.2 Experimental data of the infrared radiation characteristics of the target
由表2 数据可得,该系统可以进行目标温度和目标温度分布测量,具备目标特性分析功能,能够完成该方面的实践教学。
按照3.2 节中分析得出的迟滞计算模型式(7),测试红外探测器对黑体温度在30~100 ℃响应迟滞情况。从30 ℃开始,每间隔5 ℃记录1 次,依次获得黑体温度正(增大)行程时的探测器响应值。然后从100 ℃开始记,依次记录黑体温度反(减小)行程时的探测器响应值。所得数据如表3。
表3 红外探测器迟滞实验数据Tab.3 Experimental data of the infrared detector hysteresis
由表3 的数据可得,该系统的迟滞为0.031,能够帮助完成红外辐射特性测量系统迟滞指标测试实践教学,有助于学员进一步理解红外辐射特性测量系统。
提出了一种红外辐射特性测量系统。该系统结合红外辐射特性测量设备任务及设备结构原理,设计了近20 个操作科目,阐述了红外辐射特性测量系统标定原理,设计实现了红外辐射特性测量教学系统的软硬件,研究并确立了红外探测器指标测试操作科目的计算模型。最后进行了一系列的实验验证,结果表明该系统能满足教学训练要求。所设计的系统对于士官院校日常教学中红外辐射特性测量系统原理、指标参数的学习大有帮助。