不同背景下高超声速飞行器红外可探测性分析

于晓杰，郑永超*，郭崇岭，董士奎，杨 霄

(1.北京空间机电研究所 先进光学遥感技术北京市重点实验室，北京 100094；2. 哈尔滨工业大学 能源科学与工程学院， 哈尔滨 150001)

引 言

临近空间高超声速飞行器一般飞行在20km～100km范围内，飞行速度大于马赫数Ma=5，能够实现高速飞行、远程打击、快速突防等作战目标，具有重要的军事价值，其中以第2代猎鹰高超声速飞行器(falcon hypersonic technology vehicle 2，HTV-2)为代表[1]。该类目标飞行速度非常快，雷达散射截面非常小，使得利用雷达探测手段探测目标的难度很大。传统可见光探测高超飞行器、隐身飞机需要较高分辨率，不能用于夜间探测，且覆盖范围小，不适用于高超目标的侦察探测[2-3]。随着高超技术的发展，现有的雷达、可见光等探测手段难以有效探测该类空天目标，亟待研究发展有效探测高超目标的新技术及手段。研究发现，高超飞行器以高超声速在大气层内飞行时，由于周围空气强烈的压缩和摩擦作用，会在本体周围形成高温绕流场，此时高温绕流场内的多组分空气流会产生极强的热辐射[4-6]。这种热辐射现象为利用红外探测系统对高超飞行器的探测提供了可能途径。红外探测系统[7-9]对目标进行探测时是通过识别目标和背景的红外辐射信号来实现的，其中目标与背景对比度是一个重要的研究方向。

国内外针对巡航导弹、弹道导弹等目标的红外可探测性进行了大量的理论研究和仿真实验验证，但在高超目标红外可探测性方面的探讨比较少。WANG等人[10]分析了高超飞行器的红外辐射特征，讨论了实现该类目标红外探测预警的思路和数学模型。ZHOU等人[11]分析研究了地基红外探测器对乘波体类高超飞行器的探测性能，建立了一种新的红外系统作用距离计算模型。YANG等人[12]分析并建立了飞艇的红外探测系统对“乘波者”类高超飞行器的探测距离模型。YANG等人[13]系统分析了地基平台、浮空器、天基卫星对高超声速滑翔飞行器的红外探测能力，得出不同探测波段及不同探测平台对目标的最大作用距离。上述研究是针对单一背景进行分析，没有考虑复杂背景情况下背景辐射对高超目标红外可探测性的影响。

本文中以类HTV-2飞行器为研究对象，对飞行弹道上典型飞行条件下的红外辐射特性进行了分析和计算。通过计算卫星正俯视探测时不同飞行工况、不同背景条件下的高超目标与地球/大气背景在 1μm～14μm波段内的对比度变化规律，给出了针对类HTV-2目标滑翔段的最佳探测谱段范围。

1 地球/大气背景和类HTV-2目标的红外辐射特性计算

1.1 地球/大气背景特性的计算

类HTV-2高超目标滑翔段探测背景主要是地面和大气系统，背景的复杂多样性直接影响系统探测性能。地球/大气背景辐射包括由地球大气系统自身发射的辐射和由该系统反射和散射太阳的辐射。由于计算地球/大气背景辐射所要求的光谱分辨率不高，并且需要考虑太阳的散射以及云的影响，因此,可选择MODTRN软件中的相关模型对不同地球/大气背景的红外辐射特性进行计算。计算初始条件如下：美国1976年标准大气，地表为海洋/森林，气溶胶模式为海洋/乡村，太阳天顶角、方位角都设置为0°，云背景为晴天无积云/有积云，考虑大气多次散射。

下面用MODTRN软件计算了高超目标所在位置到探测器处不同波长λ的大气透过率，如图1所示。目标在两个典型高度30km,50km处多种条件下的背景亮度Lb如图2所示。

从图2中可以看出,差异较明显的谱段为1.10μm～2.40μm，该波段是短波红外窗口。两端受大气中水汽与二氧化碳吸收作用限制，且水汽在1.80μm处有一个吸收带，所以本窗口可分为1.50μm～1.75μm和2.10μm～2.40μm两个小窗口。在地球辐射中，0.30μm～3.00μm波段范围内主要反射太阳辐射，而地球本身及表面物体辐射极弱，可忽略不计。

Fig.1 Relationship between atmospheric transmittance and wavelength at various altitudes

Fig.2 Infrared radiation of geo-atmospheric backgrounds

a—with and without cumulus in the ocean situation b—with and without cumulus in the forest situation

2.40μm～5.00μm波段是中红外窗口。由于二氧化碳在4.30μm处吸收作用很强，所以又可分为3.40μm～4.20μm和4.60μm～5.00μm两个窗口。有积云时，海洋和森林地表下的背景辐射亮度基本一致。这主要是因为地球/大气背景红外辐射主要是云层反射太阳辐射，故海洋、森林两种背景下的红外辐射相对差异较小；无积云时，3.00μm～4.20μm范围内，森林地表相对于海洋而言，背景辐射亮度高。主要是由于对于海洋而言，在该谱段内其反射率极低，吸收率较高，而在0.70μm以上的红外区，植被有较高反射率，所以无积云时，森林相比于海洋，红外辐射更强。

8.00μm～14.00μm是热红外窗口，是热红外谱段，主要是地物本身的热辐射。在这个窗口主要是用来获得地面目标的热辐射信息，能有效地探测地面常温物体。这一波段范围内主要是由于地球本身的热辐射明显大于表面物体的辐射，所以在6.00μm～14.00μm处不同背景下的辐射亮度曲线是基本一致的。另外由于二氧化碳在9.40μm处有弱吸收带，故在9.40μm附近处会有波谷存在，与仿真图中吻合。

1.2 类HTV-2目标红外辐射特性计算

类HTV-2飞行器处于滑翔段时，头部附近来流受到强烈压缩，在头部周围形成激波层，激波层内的来流动能转化为内能，使得温度急剧上升，对飞行器形成气动加热效应。由于与周围环境相比，高超声速飞行器表面温度明显要高，故可只考虑目标自身辐射，且将目标视为灰体，发射率为常数。

计算目标温度场主要采用流热耦合的方式。通过3维软件Pro/e对目标建模，然后用 ICEM软件对其进行网格划分。使用商业软件CFD对目标蒙皮外表面的气动加热进行计算，求得蒙皮温度场[13]如图3所示。

Fig.3 Surface temperature of HTV-2 type vehicle

从图中可以看出，30km高度、马赫数Ma=7和50km高度、马赫数Ma=17两种工况对应球头部分的最高温度分别为1500K和2400K。该计算结果与相关文献[14]中报道的HTV-2本体表面所能承受的温度极限2200K相接近。根据得到的温度场分布，利用普朗克公式对目标本体辐射进行计算。

定义目标观测方向如图4所示。以目标为坐标原点，设置目标的球头方向为x轴的正半轴方向，水平方向为y轴，z轴的正半轴为竖直向上方向，θ为观测方向的天顶角，φ为观测方向的圆周角。

采用普朗克公式[15]计算目标本体的红外辐射。普朗克定律引入了微观粒子能量不连续假设，确定了黑体辐射强度与其绝对温度和波长之间的关系,其数学表达式为：

Fig.4 Diagram of the observed angle

(1)

式中,λ为波长；Eλ,b为黑体辐射光谱力，单位W/(m2·μm)；h为普朗克常数，h=6.624×10-34J·s；k为玻尔兹曼常数，k=(1.380662±0.000044)×10-23J/K;c为光速，c=(2.99792458±0.000000012)×108m/s；c1为第一辐射常数，c1=(3.7415±0.0003)×10-16W·m2；c2为第二辐射常数，c2=(1.43879±0.00019)×104μm·K。

目标的光谱辐射强度[16]定义为在一定方向上，单位投影面积内的光谱辐射总和(W/(μm·sr))，计算公式为：

Iλ=LλAt

(2)

式中，At为目标的可见面积。

设目标的光谱发射率为ελ，用朗伯源模型计算光谱辐射亮度Lλ(W/(μm·sr·m2))，计算公式为：

(3)

式中，Eλ为目标光谱辐射出射度。

设目标表面由N个表面发射率为εi的微元面组成，发射率取为0.9，沿观测方向上每个微元的可见面积为Ai(m2)，则目标沿观测方向上的光谱辐射强度Iλ(W/(μm·sr))由(1)式～(3)式可知：

(4)

图5中分别给出了目标在30km高度、马赫数Ma=7和50km高度，马赫数Ma=17两个工况，x轴、y轴和z轴3个特定观测方向下的1μm～14μm波段内的光谱辐射强度计算结果。在相同高度下，俯视观测即z轴观测方向(θ=0°,φ=0°)时，目标的光谱辐射强度最大，前视即x轴观测方向(θ=90°,φ=0°)时辐射强度最小。这是因为探测角度的不同，导致飞行器的可见表面积和对应的温度分布也有所不同。俯视观测光谱强度较高，是其可见表面积明显大于另外两个角度的可见表面积。30km高度、马赫数Ma=7工况下辐射光谱的峰值波长约在3.00μm附近，而50km工况峰值向左移动约1μm，这是因为壁面光谱辐射峰值与表面温度相关，且温度升高峰值波长变短，这也符合维恩位移定律的描述。观测角度相同时，50km高度、马赫数Ma=17工况下的光谱辐射强度要高于30km高度、马赫数Ma=7工况。

Fig.5 Spectral radiant intensity with different observed anglesa—30km，Ma=7 b—50km，Ma=17

目标在30km、50km两个典型飞行高度时，改变观测方向的天顶角(0°～180°)和圆周角(0°～360°)，计算得出3μm～5μm，8μm～12μm波段内的总辐射强度大小如图6～图9所示。相同条件下，3μm～5μm波段内总的辐射强度比8μm～12μm波段内的大。目标在50km高度下各个方向的辐射强度都要高于30km，这主要是由于50km高度、马赫数Ma=17工况下，气动加热效应更强，导致表面温度更高。

Fig.6 Total radiation intensity of 3μm～5μm in 30km，Ma=7

Fig.7 Total radiation intensity of 8μm～12μm in 30km，Ma=7

Fig.8 Total radiation intensityof 3μm～5μm in 50km，Ma=17

Fig.9 Total radiation intensity of 8μm～12μm in 50km，Ma=17

2 目标背景对比度计算及分析

目标背景对比度是表征探测系统性能的一个重要指标，其值越大，越容易实现对目标的探测。目标与背景对比度受多种因素影响，其与大气透过率有关，不同大气背景、不同高度的大气透过率不同，产生的背景辐射也是不同的；与目标及探测器的相对位置也有关系。在点目标探测系统中，目标背景对比度如下[17-19]：

(5)

式中，It表示目标红外辐射强度，τa表示大气透过率，Lb表示背景辐射亮度，At表示目标沿视线方向投影面积，Ap表示瞬时视场在高超目标所在位置投影面积。

以低轨道探测器为例，其轨道高度为1000km，分辨率为100m，目标位于卫星正下方海拔30km，50km高度。利用(5)式计算目标背景对比度光谱变化规律，如图10所示。

Fig.10 Target-background contrast

a—without cumulus in the ocean situation b—with cumulus in the ocean situation c—with out cumulus in the forest situation d—with cumulus in the forest situation

从图10中可以看出,在2.50μm～3.00μm谱段范围内的目标背景对比度较大，在1.80μm～2.00μm,3.00μm～3.50μm及4.30μm～4.50μm处的目标背景对比度也有一定的峰值，但与2.50μm～3.00μm相比量级太小，可以忽略不计。为了更精确的分析，下面计算了2.50μm～3.00μm处的目标与背景对比度，如图11所示。

Fig.11 Target-background contrast in the vicinity of 2.7μm

a—without cumulus in the ocean situation b—with cumulus in the ocean situation c—without cumulus in the forest situation d—with cumulus in the forest situation

从图10、图11的光谱变化中可分析出：(1)类HTV-2高超目标在30km高度、马赫数Ma=7和50km高度、马赫数Ma=17两种工况下，2.70μm波段附近目标背景对比度先增大后有所下降；50km高度、马赫数Ma=17工况下的目标背景对比度远大于30km高度、马赫数Ma=7工况下的；不同背景下的30km高度、马赫数Ma=7工况的目标背景对比度基本一致，这说明目标在这种工况下，背景对目标背景对比度的影响很小;(2)当探测背景无积云时，海洋地表下的目标背景对比度相比于森林地表情况下是基本一致的；有积云时，也符合上述情况,这说明地表类型对目标背景对比度的影响不大;(3)在不同工况和任一背景下，目标背景对比度在2.65μm～2.85μm这个谱段范围内都较大，所以不考虑其它影响因素时，可将其作为设计红外探测系统时的参考。

3 结 论

计算了高超目标在两种典型飞行工况下的红外辐射特性、不同背景的辐射强度变化。以此为基础，计算获得了多种情况下的目标背景对比度的变化曲线。通过对仿真结果的分析，得出2.65μm～2.85μm这个谱段范围更适合卫星探测。由于本文中只是从目标背景对比度单一角度开展卫星最佳谱段的定量研究，所以在接下来的工作中要综合考虑信噪比、探测系统灵敏度等探测指标的影响，从多角度深入开展研究。