大气气溶胶红外散射透过率计算研究

吴辉阳，王泽洋，黄兴军，张 艳，石孝楠，张云飞

大气气溶胶红外散射透过率计算研究

吴辉阳1，王泽洋1，黄兴军2，张 艳2，石孝楠1，张云飞1

（1. 北京航空航天大学 航空科学与工程学院，北京 100191；2. 北京机电工程研究所，北京 100074）

红外辐射在大气中传播时大气气溶胶散射是能量衰减的原因之一。通过考虑气溶胶密度随高度变化，结合大气能见度参数，建立了红外线在水平均匀传播和斜程传播下的大气气溶胶散射透过率计算方法。水平均匀传输下，在中长波波段对大气气溶胶散射透过率采用常规积分求均值方法计算，与采用波长中值的工程计算公式结果对比，表明工程计算公式具有足够的工程精度；对于斜程传播情况，通过对高度积分得到红外线大气气溶胶散射透过率的斜程工程计算公式。计算分析了探测器高度、波段和大气能见度对大气气溶胶散射透过率的影响。本文对于建立自主的精确大气透过率计算模型、计算机载探测器对红外目标的下视探测性能以及深入了解大气气溶胶散射透过率的影响因素具有重要意义。

红外辐射；气溶胶密度；大气散射透过率；能见度

0 引言

在研究红外探测器的探测性能以及现代飞行器的红外隐身性能时，通常依据红外作用距离方程[1-2]，这要考虑3个方面的因素：①红外飞行器的红外辐射特性，与目标本身的红外辐射特性有关；②红外探测器的性能参数，与探测器有关，可以采用所谓红外探测器系统特征常数来描述[3]；③大气衰减和背景辐射干扰，它会降低目标的辐射强度和干扰探测器的探测性能，与周围环境的大气条件有关。

大气中的分子和粒子对红外线辐射强度的减弱作用定义为红外线在大气中传播一定距离后的辐射能与入射辐射能的比值，用大气透过率表示。

大气衰减作用有3个方面[4]：①大气中某些气体分子的吸收；②大气分子、气溶胶的散射；③气象条件（雨、雪）造成的衰减。大气光谱透过率可表示为：

a()＝吸()×散()×气象() (1)

式中：吸()、散()、气象()分别为大气吸收、大气散射和受气象条件影响的大气气象光谱透过率。

大气中的水蒸气（H2O）、二氧化碳（CO2）、臭氧（O3）以及其他微量分子如CO、NO等具有极性的分子对于红外辐射在大气中的传输特性具有吸收作用。其中，水蒸气和二氧化碳因在大气中的含量较高且稳定，因而对于红外线具有最强的吸收作用。因此大气吸收透过率可写为：

大气散射衰减是大气分子散射、云和雾以及气溶胶等颗粒散射造成的红外辐射衰减。对于红外波段来说，由于大气分子的尺度远小于红外线波长，研究其散射吸收是没有意义的，因此在研究大气对红外辐射的散射衰减时，只考虑颗粒散射的影响。

国内外针对于大气吸收作用影响的研究有大量的著作和文献[5-9]，对于大气气溶胶散射的影响也有许多研究。张建奇[4]给出了气溶胶密度随高度的变化规律，方义强[10]等在对天空大气背景的红外辐射的理论分析及分层数值计算中，给出了考虑高度影响的大气气溶胶散射透过率计算公式。易亚星[11]在分析影响目标红外辐射亮度的因素时，计算了大气气溶胶散射波段透过率随路径的变化，两者大致呈线性递减关系，文中公式虽然考虑了高度影响，但并没有给出计算结果及高度影响的规律。魏合理[12]等研制了一款通用辐射大气传输计算软件CART，考虑了大气吸收、大气散射、高度等对大气透过率的影响，并与国外的计算软件进行了对比。国内外研发的大气背景红外辐射的多种计算软件，如LOWTRAN、MODTRAN、FAS-CODE、CART等是建立在大量观测数据的基础上，有较高的可信度，但这些计算模型对于用户来说是不可见的，不利于问题的分析研究以及实际应用的借鉴和利用。总之，各种研究考虑红外辐射在斜程传播中高度变化影响的研究较少，模型不透明，有哪些影响参数，影响程度如何也不甚清晰。

现代红外探测与制导武器发展迅速，而且预警机及战斗机搭载红外探测器的情况也日趋普遍。因此在分析地面红外探测器探测中高空目标以及机载红外探测器下视探测红外目标时，必须考虑高度对于大气吸收透过率以及大气气溶胶散射透过率的影响。

本文建立了红外辐射水平均匀传输及斜程传输情况下的大气气溶胶散射透过率计算方法，并对高度、波段和能见度的影响进行了综合计算分析。

1 大气气溶胶散射透过率影响因素计算

大气中的各种颗粒对红外线的散射作用称为粒子散射或米氏（Mie）散射，其影响除了由云和雾这些不确定的气象因素所造成外，更为常见的是由大气中广泛存在的气溶胶颗粒所造成。为区别于衡量大气吸收衰减作用的大气吸收透射率，将大气颗粒散射所引起的衰减作用通过大气气溶胶散射透射率来衡量。本文只研究气溶胶颗粒的影响。

1.1 水平均匀传输下大气散射透过率计算

在红外波段，大气气溶胶散射对辐射造成的衰减可用散射透过率表示[4]。当一束单色辐射在均匀介质中传播距离后，由于散射作用将使辐射按指数规律衰减，因此在仅含散射物质（即无吸收物质）的均匀大气中传播时，通过距离的纯散射的介质透过率为：

式中：(0)、()分别为经过距离前、后的光谱辐射功率，()为散射系数[13]：

()＝p()2(4)

式中：为每立方厘米内的微粒数（气溶胶浓度）；()为散射效率因子，对于单粒径均匀粒子而言，它为散射截面和粒子几何截面之比，与气溶胶颗粒的半径、介质的折射率、波长有关；为散射粒子的半径。由于散射系数随波长变化，因此透射率也是波长的函数。

气溶胶是指悬浮在气体中的小粒子，尺度范围为10－3～10mm。在近地面大气中气溶胶的浓度约为每100～1000个/cm3，随高度呈指数递减。气溶胶浓度随高度变化一般表示为拟合公式：

式中：(0)为海平面的气溶胶浓度；为大气层高度；*为标高，它与地面能见度大小有关，当能见度为2～6km时，变化范围为0.8～1km，当能见度为6～25km时，变化范围为1～1.4 km，这时可取*=1.2km[14]。综合(3)～(5)式，某一海拔高度处薄层大气中，红外辐射传播距离的大气气溶胶散射透过率为：

式中：(0,)是海平面传播距离的大气气溶胶散射透过率。

气象学上常采用气象能见度的方法来表征不同气溶胶浓度对大气散射的影响，该方法具有足够的精度且方便快捷[15]。此时海平面处散射系数与能见度关系为：

式中：0为大气吸收谱线的波长，通常取0.55mm或0.61mm；为经验系数。当大气能见度特别好（≥60km）时，＝1.6；中等能见度时，＝1.3；能见度很差（≤6km）时，＝0.5851/3[13]。因此，在高度处红外辐射传播距离时的大气气溶胶散射光谱透过率为：

对波长积分可得1～2波段的大气气溶胶散射透过率：

因此可以计算出长波和中波波段、几种典型能见度情况下，高度对大气气溶胶散射透过率与距离关系曲线的影响（图1～图3）。

图1 海拔高度对大气气溶胶散射透过率影响（两波段，=23km）

图2 海拔高度对大气气溶胶散射透过率影响（两波段，=10km）

可见，随高度增大，大气气溶胶散射透过率将会增大。如果不考虑高度对大气气溶胶散射透过率的影响，则透过率计算值将会偏小，目标的红外辐射强度将会被低估。另外，其他情况相同时，能见度越小，高度对大气气溶胶散射透过率影响越大，而8～14mm长波大气气溶胶散射透过率比3～5mm短波都要高一些。

式(9)的积分对于工程计算应用来说较为麻烦，故在此对其进行简化。

对于中等能见度情况取＝1.3，由拉格朗日中值定理，高度处红外辐射传播距离为的大气气溶胶散射波段透过率为：

式中：为波段系数，经数值积分法计算，8～14mm波段取中值波长11mm时，＝0.080，3～5mm波段取中值波长4mm时，＝0.30。

因为气溶胶散射透过率随波长变化较慢，如图4所示，因此在衰减不太大的情况下，采用中值法的计算结果与精确算法的非常接近，误差很小。采用8～14mm波段的中值波长11mm及对应波段系数计算的大气气溶胶散射波段透过率曲线如图5所示。

1.2 斜程传输下大气散射透过率计算

在高度处红外辐射传播距离为d的大气散射波段透过率为：

因为d＝dh/sin，故可转化为：

图4 大气气溶胶散射透过率与波长关系（=1km，=10km，=10km）

图5 海拔高度对大气气溶胶散射透过率影响（=11mm，=10km）

探测器到目标之间总的大气散射波段透过率为：

式中：a为探测器高度；t为目标高度；为视线仰角，即红外线与水平面的夹角。

按上式可以计算出长波和中波波段、几种典型能见度情况下，视线仰角对大气散射透过率与距离关系曲线的影响（图6，图7）。这里探测器位于海平面高度。

图6 视线仰角对大气气溶胶散射透过率影响（两波段，=23km）

图7 视线仰角对大气气溶胶散射透过率影响（两波段，=10km）

可见，当视线仰角不变时，随着传播距离的增大，大气散射透过率逐渐减小。当传播距离一定时，随视线仰角增大，大气散射透过率将会增大。

1.3 大气散射透过率计算工程方法

对于水平均匀传输和斜程传输，大气气溶胶散射透过率工程计算的公式分别为式(10)和式(13)。

利用这两个计算公式，计算得到探测器在海平面高度，传输距离20km时，不同波段和能见度情况下，大气散射波段透过率随仰角的变化（图8）。

可见，在波段和能见度一定情况下，仰角越大，红外辐射传输所处高度越高，则大气气溶胶散射波段透过率越大。海平面情况（图中0°视线仰角）的透过率是最小的，在仰角不为0时的大气气溶胶散射波段透过率都比海平面的要高，仰角越大，透过率越大，而且越接近于1。与水平均匀传输情况类似，能见度越低，波长越短，仰角对大气气溶胶散射波段透过率的影响越大。

2 结论

1）通过理论分析与工程近似推导，得到了红外线水平均匀传输情况和斜程传输情况的大气气溶胶散射透过率的工程计算公式，考虑了红外波段、大气能见度、海拔高度、视线仰角、传输距离等因素的影响；

2）对于大气气溶胶散射波段透过率，各种情况下，该值随海拔高度而增大（当海拔高度达到5km以上时，该值接近于1），随视线仰角增加而增大，随红外波长增加而增大，随能见度提高而增大；

3）本文结果完善了大气透过率的计算方法，可用于计算地面红外探测器探测空中目标以及机载红外探测器下视探测地空和地面红外目标的探测性能，对于飞行器的红外隐身性能计算与评估具有重要理论意义与工程应用价值。

[1] Richard D Hudson, Jr.[M]. John Wiley & Sons Inc., 1969.

[2] 牟达, 王建立, 陈涛. 凝视型红外搜索跟踪系统对高速飞机作用距离的分析[J]. 光学技术,2007, 33(3): 420-423.

MOU Da, WANG Jianli, CHEN Tao. Analysis of the range of high-speed aircraft by staring infrared search and tracking system [J]., 2007, 33(3): 420-423.

[3] 赖德雄, 张云飞. 飞行器红外隐身性能评估方法研究[J]. 红外技术, 2014, 36(5): 398-403.

LAI Dexiong, ZHANG Yunfei. Research on the evaluation method of aircraft infrared stealth performance[J]., 2014, 36(5): 398-403.

[4] 张建奇, 方小平. 红外物理[M]. 西安: 西安电子科技大学出版社, 2007:143.

ZHANG Jianqi, FANG Xiaoping.[M]. Xi'an: Xidian University Press, 2007: 143.

[5] 彭德权, 周成平, 丁明跃. 海空背景长波红外大气透过率的仿真计算方法[J]. 红外与激光工程, 2001, 30(2): 112-117.

PENG Dequan, ZHOU Chengping, DING Mingyue. Simulation calculation method of atmospheric transmittance of long-wave infrared in sea and sky background[J]., 2001, 30(2): 112-117.

[6] 王迪, 李承芳, 熊飞. 大气红外辐射传输的简便算法与MODTRAN的比较[J]. 光学技术, 2006, 32(8):293-295.

WANG Di, LI Chengfang, XIONG Fei. Comparison of a simple algorithm for atmospheric infrared radiation transfer and MODTRAN [J]., 2006, 32(8):293-295.

[7] 毛峡, 董旭阳, 刁伟鹤. 巡航导弹红外辐射及大气衰减计算模型[J]. 北京航空航天大学学报, 2008, 34(8): 952-955.

MAO Xia, DONG Xuyang, DIAO Weihe. Calculation model of cruise missile infrared radiation and atmospheric attenuation[J]., 2008, 34(8): 952-955.

[8] 周国辉, 刘湘伟, 徐记伟. 一种计算红外辐射大气透过率的数学模型[J]. 红外技术, 2008, 30(6): 331-334.

ZHOU Guohui, LIU Xiangwei, XU Jiwei. A mathematical model for calculating the atmospheric transmittance of infrared radiation[J]., 2008, 30(6): 331-334.

[9] 吴晗平. 红外搜索系统[M]. 北京: 国防工业出版社, 2013. WU Hanping.[M]. Beijing: National Defense Industry Press, 2013.

[10] 方义强, 樊祥, 程正东, 等. 天空大气背景红外辐射建模与计算[J]. 激光与红外, 2013, 43(8): 896-902.

FANG Yiqiang, FAN Xiang, CHENG Zhengdong, et al. Modeling and calculation of sky atmospheric background infrared radiation[J].d, 2013, 43(8): 896-902.

[11] 易亚星, 姚梅, 吴军辉, 等. 影响红外目标探测亮度的因素[J]. 红外与激光工程, 2014, 43(1): 13-18.

YI Yaxing, YAO Mei, WU Junhui, et al. Factors affecting infrared target detection brightness[J]., 2014, 43(1): 13-18.

[12] 魏合理, 陈秀红, 饶瑞中. 通用大气辐射传输软件CART介绍[J]. 大气与环境光学学报, 2007, 2(6): 446-450.

WEI Heli, CHEN Xiuhong, RAO Ruizhong. Introduction to the general atmospheric radiation transfer software CART[J]., 2007, 2(6): 446-450.

[13] 潘皓. 云背景红外辐射计算及图像仿真[D]. 西安: 西安电子科技大学, 2006.

PAN Hao. Infrared Radiation Calculation and Image Simulation in Cloud Background[D]. Xi'an: Xidian University, 2006.

[14] 杨风暴.红外物理与技术[M]. 北京: 电子工业出版社, 2014.

YANG Fengbao.[M]. Beijing: Electronic Industry Press, 2014.

[15] Paul W Kruse, Laurence D McGlauchlin, Richmond B McQuistan.[M]. New York: John Wiley, 1962.

Calculation of Infrared Scattering Transmittance of Aerosol

WU Huiyang1，WANG Zeyang1，HUANG Xinjun2，ZHANG Yan2，SHI Xiaonan1，ZHANG Yunfei1

(1.,,100191,;2.,ing 100074,)

Atmospheric scattering is one of the factors influencing the energy attenuation of infrared radiation in the atmosphere. Considering the variation in the aerosol densitywith altitude, the calculation method of atmospheric scattering transmittance is developed for the horizontal homogeneous transmission and oblique transmission integrating the atmospheric visibility parameters. In the case of horizontal homogeneous transmission, the average transmission method of atmospheric aerosol scattering transmission in the mid-long wave band is calculated by performing conventional integration. Compared with the result of the engineering calculation formula using the median wavelength, the engineering calculation formula is shown to have sufficient engineering accuracy. For the case of slope propagation, a slope engineering calculation formula for infrared atmospheric aerosol scattering transmittance was obtained by integrating the altitude. The effects of the altitude, waveband, and visibility on the atmospheric scattering transmittance are calculated and analyzed. This study has great significance for the establishment of an independent accurate atmospheric transmittance calculation model, the downward-looking detection performance of computer-borne detectors for infrared targets, and provides a good understanding of the factors influencing atmospheric aerosol scattering transmittance.

infrared radiation, aerosol density, atmospheric scattering transmittance, visibility

TN012

A

1001-8891(2022)02-0151-05

2020-03-30；

2020-06-04.

吴辉阳（1996-），男，河南周口人，硕士研究生，主要从事红外隐身和大气红外辐射研究。E-mail：wuhuiyang@buaa.edu.cn。

张云飞（1964-），男，湖南长沙人，教授，研究方向为飞机总体设计、隐身技术、无人机设计。E-mail：cloud_zhang@buaa.edu.cn。