王志会,屈鸿伟,张晶涛,王东杰,邓骏杰
(空军预警学院,湖北 武汉 430019)
随着红外探测技术的不断发展,通过对目标和背景的红外辐射特性分析,已经成为有效探测与识别目标的重要方法。尤其对于夜间作战,依据红外辐射特性可清晰掌握目标和战场情况[1]。大气透过率对空中目标红外特性的传输产生重要的影响,而红外探测系统接受的信号强度主要由目标红外辐射强度和大气透过率共同决定。尤其对军事领域中的飞机、导弹等红外目标实现精确地跟踪、警戒和搜索,必须考虑大气透过率的影响,而且大气透过率往往为决定性因素[2]。
关于大气透过率的获取,目前常用的主要的方法有三种:数值模拟计算[2-4]、软件计算(如LOWTRAN、MODTRAN)和实地测量[2]。其中,数值模拟方法主要结合经验公式进行计算,在计算过程中,利用测得的大气水蒸气含量[5],可以快速地计算出不同地区的红外波段大气透过率。路远等人利用实验数据计算求得大气中水蒸气和CO2的光谱吸收系数,并通过数值模拟的方法计算出海平面水平路径传输情况下的大气吸收透过率[6]。张亮、毛峡等人充分考虑温度和气压随海拔高度变化对大气吸收透过率的影响,建立水平传输及倾斜传输两种情况大气吸收透过率的计算模型[7-8]。实际中,在低层大气环境下,大气密度随海拔高度而变化,其对水蒸气及CO2的吸收透过率也会产生重要的影响,不容忽视[8]。基于温度和气压对大气吸收透过率的影响,本文增加大气密度因素的影响,建立水蒸气、CO2红外吸收透过率计算数值模型,分析了水平路程传输及倾斜路程传输两种情况下的大气吸收透过率,并将计算结果与MODTRAN5软件计算结果进行比较,充分验证了该模型的有效性。
实际情况下,大气中的分子对红外辐射的吸收往往为选择性吸收,而这些吸收峰与分子的共振频率是对应的。分子振动和自旋共振频率都处于红外区域,而且主要分布在2~15 μm。大气中,水蒸气是吸收红外辐射的主要因素,且主要集中在2~3 km大气层下;而大气中的CO2虽然只占大气气体的0.03 %~0.05 %,但它却是红外大气吸收的另一重要因素。此外,在低层大气环境中,氧化氮、一氧化碳等气体产生的吸收比较小,通常可以忽略不计[9]。图1为大气垂直分层情况。
图1 大气垂直分层示意图
水蒸气对红外辐射的吸收与其传输的路程中水分子的总量有直接关系。对于大气一段路程中水汽的含量通常用可凝结水毫米数来表示,它可以根据大气温度和相对湿度计算得到。即水蒸气在大气中的含量ω(mm)为:
ω=f·Hr·R
(1)
相应的水蒸气吸收透过率为:
τH2O(λ)=exp(-u(λ)·ω)
(2)
式中,f(mm/km)为环境温度T(℃)时的饱和空气中水蒸气质量;Hr为相对湿度;R(km)为海平面水平传输距离;u(λ)(km-1)为水蒸气单色红外线性吸收系数。实验数据表明,在红外辐射通过水蒸气含量相等的路程情况下,大气透过率基本相等[10]。由此可以得出水蒸气大气吸收透过率为:
(3)
式中,u1(λ)(km-1)为大气温度为5 ℃,相对湿度为100 %时海平面水平路径情况下水蒸气的单色红外线性吸收系数;f为大气温度为T时饱和空气中的水蒸气质量;f0(mm/km)为大气温度5 ℃时饱和空气中的水蒸气质量,通过查询文献[5]可得到相应数据。
由于水蒸气对红外辐射的吸收随温度与气压变化相应改变,因此当高度变化明显时需要进行相应的修正。在计算过程中,由于吸收能力受温度的影响表较小(≤4 %),可以忽略不计,所以往往只需要考虑气压的影响。假设ωe(mm)表示辐射传输路程中依据吸收能力折算到大气近地层水蒸气的等效含量,而ωH(mm)表示H(km)高度下水蒸气含量,则有:
ωe=ωHβH2O
(4)
式中,βH2O为高度修正系数,可由经验公式得到:
βH2O=exp(-0.0654H)
(5)
考虑大气密度随高度而降低,大气湿度随高度的变化满足经验公式[9]:
Ha,H=Ha,0exp(-0.45H)
(6)
式中,Ha,H(mg/L)为高度为H处的绝对湿度;Ha,0(mg/L)为近地处或海平面处的绝对湿度。
综合水蒸气的吸收能力和水蒸气量随高度的变化,可将高度H的辐射沿着水平传输路程的水蒸气有效含量表示为:
ωe=f·Hr·R·exp(-0.5154H)
(7)
在倾斜传输路程中,水蒸气的有效含量可通过以下方法计算得到。在高度为h(km)处的大气单元层dr(km)(图2)中的水蒸气含量为:
图2 辐射传输中(倾斜路线)吸收物质层有效厚度
dωe=f·Hr·exp(-0.5154h)dr
(8)
式中,h=Rsinθ;θ(°)为观测高度角。
忽略地面曲率影响,对式(8)进行积分可求得从高度H1(km)到H2(km)的倾斜路程中水蒸气含量为:
(9)
将公式(7)、(9)分别代入公式(3)中,可得不同高度H水平路程传输的水蒸气吸收透过率为:
(10)
倾斜路程传输的水蒸气吸收透过率为:
(11)
λ1~λ2(μm)波段内的水蒸气平均吸收透过率为:
(12)
除了水蒸气,CO2是大气中吸收红外辐射最多的因素。在低层大气中,CO2含量相对稳定,因此由CO2吸收所产生的辐射衰减近似认为与气象条件无关,吸收透过率仅仅和辐射传输的距离有关,即:
τCO2(λ)=exp(-∂(λ)·R)
(13)
式中,∂(λ)(km-1)为CO2单色红外线性吸收系数。大气温度5 ℃,相对湿度100 %时海平面水平路径CO2的红外单色线性吸收系数可通过文献[6]查得。
考虑高度对CO2红外吸收能力的影响,对CO2传输距离进行高度修正,得到以下经验公式:
Re=RH·e-0.19H
(14)
式中,H(km)为红外辐射传输的高度;RH(km)为高度H上红外辐射传输的距离,Re(km)为等效成近地层之后传输的有效距离。
由于空气压强随海拔高度成指数变化:
PH=P0·e-0.123H
(15)
式中,P0(Pa)为海平面上的气压;PH(Pa)为高度H处的气压。
因此,综合考虑CO2的吸收能力及质量随海拔高度的变化,可得到折算成近地层水平路程传输的等效距离:
Re=RH·e-0.19H·e-0.123H=RH·e-0.313H
(16)
在倾斜路程传输的情况下,类似得到折算成近地层水平路程传输的等效距离为:
(17)
将式(16)、(17)分别代入式(13)中,得到高度H情况下CO2水平路程传输的吸收透过率为:
τCO2(λ)=exp(-∂(λ)·RH·e-0.313H)
(18)
倾斜路程传输的吸收透过率为:
(19)
λ1~λ2波段内的CO2平均吸收透过率为:
(20)
综合考虑水蒸气、CO2的吸收,可得低层大气吸收透过率:
τz(λ)=τH2O(λ)·τCO2(λ)
(21)
在λ1~λ2波段内的低层大气平均吸收透过率为:
(22)
考虑空中军事目标红外辐射主要分布在3~5 μm和8~12 μm波段[7,9,11-12],因此实验重点对这两个波段的大气吸收透过率进行了仿真计算,进行MODTRAN5软件计算时,选用1976年美国标准大气模式。
设置水平路程传输距离为1 km,近地层大气相对湿度为46 %,大气温度为15 ℃(由文献[5]可得该温度对应的水蒸气质量),探测系统高度分为0.2 km、3 km、8 km三种情况。仿真结果如图3、图4(左图为水蒸气吸收透过率,右图为CO2吸收透过率)所示。图3为不同海拔高度3~5 μm波段仿真计算结果,图4为不同海拔高度8~12 μm波段仿真计算结果。
图3 3~5 μm波段不同海拔高度水蒸气和CO2吸收透过率
图4 8~12 μm波段不同海拔高度水蒸气和CO2吸收透过率
由图3、图4可以看出,水平路程传输情况下,随着海拔高度的增加,水蒸气及CO2的吸收透过率在不断增大,且通过对大气密度随高度变化影响的修正,吸收透过率得到了明显提升,其中水蒸气的提升幅度大于CO2,这是因为水蒸气质量受海拔高度的影响大于CO2。
综合水蒸气及CO2的吸收,计算不同海拔高度大气平均吸收透过率,结果如表1所示。
表1 不同海拔高度大气平均吸收透过率计算结果
通过表1比较不同海拔高度情况下修正前后与MODTRAN5软件计算的结果,利用公式(23)计算出相对误差曲线如图5所示:
(23)
图5 不同海拔高度大气平均吸收透过率修正前后相对误差计算结果
由图5可以看出,水平传输路程情况下,在海拔高度为0.2 km时,大气密度受高度的影响比较小,所以修正前后的相对误差都比较小,在3~5 μm波段均小于2 %,其中出现修正后的误差大于修正前的情况是因为MODTRAN5软件本身在计算该波段吸收透过率时存在一定误差,尤其是针对CO24.3un吸收带时,计算误差较大[13]。但总体来看,随着海拔高度的不断增加,当高度达到8 km时,修正后的大气平均吸收透过率计算结果明显优于修正前,该现象在8~12 μm波段更为清晰。
设置H1为0,H2分别为1 km、2 km、5 km,观测高度角θ为90°,即垂直路程传输,近地层大气相对湿度为46 %,大气温度为15 ℃,仿真结果如下图6、图7(左图为水蒸气吸收透过率,右图为CO2吸收透过率)所示。图6为不同传输距离3~5 μm波段的仿真计算结果,图7为不同传输距离8~12 μm波段的仿真计算结果。
图6 3~5 μm波段不同垂直传输距离水蒸气和CO2吸收透过率
图7 8~12 μm波段不同垂直传输距离水蒸气和CO2吸收透过率
通过观察图6和图7可以发现,在垂直路程传输情况,随着传输距离的增加,水蒸气及CO2的吸收透过率在不断减小,而修正前后的差别也越来越大,从而表明在垂直路程传输情况下,随着传输距离的增加,受大气密度的影响越来越大,主要是因为传输距离的增加使得大气密度对吸收透过率造成的影响积累的越来越大。
综合水蒸气及CO2的吸收,计算不同垂直传输距离大气平均吸收透过率,结果如表2所示。
表2 不同垂直传输距离大气平均吸收透过率计算结果
通过表2比较不同传输距离情况下修正前后与MODTRAN5软件计算的结果,利用公式(23)计算出相对误差曲线如图8所示:
图8 不同垂直传输距离大气平均吸收透过率修正前后相对误差计算结果
同样,由图8可以看出,在垂直路程传输情况下,当传输距离为1 km时,修正前后的计算误差都比较小,小于5 %,其中在3~5 μm波段也出现修正前优于修正后的情况。但随着垂直传输距离的增加,修正后的计算结果明显优于修正前,当传输距离为5 km时,修正前的相对误差超过了25 %。因此,在远距离传输过程中,计算大气吸收透过率必须考虑大气密度的影响。
总之,在低层大气环境中,无论是计算不同海拔高度水平传输情况下的大气吸收透过率,还是计算不同传输距离倾斜路程传输情况下的大气吸收透过率,大气密度都会造成一定的影响,且随着海拔高度及传输距离的增加,影响会越来越大,需要进行修正。
红外辐射在通过大气时会受到大气吸收、散射、反射和漫射的影响。其中,大气吸收作为导致低层红外辐射衰减的重要因素,对其进行分析计算是研究红外大气传输的关键,也是有效进行红外目标预警的基础。本文针对低层大气中水蒸气和CO2两个吸收红外辐射的主要因素,在考虑分子吸收能力受高度影响的基础上,增加大气密度随高度变化的影响,得到修正后的低层大气红外吸收透过率计算模型。该模型在获得大气温度、相对湿度的情况下,便可以计算出任意海拔高度水平路程传输及任意观测角倾斜路径传输的低层大气吸收透过率。最终,通过实验验证了该模型的有效性。由于大气影响因素分析中引用了经验公式,所以一定会影响到计算结果的精度,但对于精度要求不高的工程计算,基本上可满足要求,具有较好的参考价值。