詹杰  郭瑞鹏  魏合理 黄宏华 徐青山 饶瑞中

（安徽光学精密机械研究所 合肥 230031）

摘  要:本文详细叙述了利用Bear-Langley-Bouguer定律测量整层大气透过率的原理。从硬件实现方法上分别介绍了晴朗天气下白天利用太阳直接辐射测量太阳的连续谱和晚上利用ICCD测量恒星分波段辐射值,通过定标实时计算整层大气透过率。文中概述了从可见到近红外波段白天透过率测量的硬件实现;晚上可根据待测天区的需要,任意选择该天区内的0等到6等恒星，实时测量该天区方向的大气透过率。测量结果和相关仪器设备测量结果对比，证明该大气透过率的测量技术是可行的。

关键词:大气透过率，增强CCD(ICCD)，辐射，定标，整层大气。

1.引  言

在地基测量空间目标的光学特性研究中，由于受到大气的影响，在同样照明条件下，大气的传输特性不同，导致测量得到的目标光学特性有不同的表观结果。大气传输的影响是空间目标光学特性定量测量所必须扣除的因素。大气光谱透过率是反映大气辐射传输性质的重要参数，在大气辐射、地球资源遥感、空气质量监测等问题的研究中具有重要意义[1,2]。大气是一个非常复杂的系统，在不同的地区、环境、季节、时间等条件下，其差别是非常巨大的，大气的变化导致光在大气中的传输有显著不同。美国于上个世纪70年代开发了低分辨率大气透过率计算软件LOWTRAN，以及九十年代的升级版本MODTRAN，这些程序根据不同地理纬度和季节建立了不同的大气模型，考虑了各种影响大气的不同因素，是目前最完整的一种计算大气透过率和背景辐射的软件。然而我国幅员辽阔，不同地区的地理条件差异非常明显，大气条件同样千差万别，仅仅依靠统计的、平均的理论模式难以对千变万化的大气现象进行精确的定量描述，因此实时测量显得尤为重要。大气透过率实时测量可以实时测量出具体地基、具体时域、具体天区内的大气透过率，由此可以建立可靠的整层大气透过率测量系统，并掌握大气传输特性测量分析的关键技术和方法。可对其它地基测量的光度数据进行大气传输修正，保证观测数据的唯一性和可靠性；通过对利用大气透过率谱和光学传感器的响应特性分析，为选择优化探测光谱波段、提高光学传感器的探测性能提供设计参考。

大气透过率是一个多因子函数，它与探测路径长度、视角、辐射波长、气压、气温和大气成分等有关。大气透过率衰减由大气分子和悬浮在大气中的气溶胶颗粒的吸收和散射造成的。大气的吸收衰减主要是由于大气中水汽H2O、CO2和O3的强烈辐射吸收造成[3]。大气的散射衰减主要是由于大气分子和大气中悬浮的气溶胶粒子如雾和雨滴对光波的散射造成的，在可见光波段（0.38μm -0.77μm）是一个大气窗口，大气吸收在该窗口很少，辐射在该波段的衰减主要是由于大气分子和气溶胶粒子散射造成，因此，准确、实时测量和监测可见光波段的大气中传输衰减，对研究大气的光学特性、测量气溶胶浓度以及辐射测量的大气修正是十分必要的。测量数据可在遥感技术、气象等关系国民经济建设的其它领域也有广泛的应用。

大气透过率的测量方式通常是通过太阳辐射计对太阳直射光谱进行测量[4,5]，本文的研究内容是在可见光到近红外波段在同一个系统中实现白天和晚上透过率的测量，白天利用光栅和线阵CCD测量太阳的连续辐射，晚上利用望远镜从选定星表中记录的恒星光谱和增强CCD（ICCD）非常直观地分波段测量恒星的弱辐射，利用长法计算大气透过率。

2．大气透过率测量的定标方法

在可见光和近红外波段，经过整层大气到达地表的太阳光或星光（以下简称目标）光谱强度可以表示为：

                             (1)

其中，为该目标的大气外界的星光光谱强度，（假定目标的大气外光谱已知），为总的垂直方向大气消光光学厚度（大气衰减），为大气质量，为天顶角，。大气衰减包括大气分子消光和气溶胶消光。

仪器接收星光强度，转化成电信号，其电压输出为：

                        (2)

定标的目的：

1）, 得到仪器的响应常数；

：                                   (3)

其中，已知，:定标拟合直线在Y轴的截距值B, ,知道、就可以计算出。不同的对应不同的。

2）, 对晚上恒星测量，由每颗星的天外谱得到该颗星的无大气衰减的

仪器输出信号：                                           (4)

这样，晚上由任何时刻测量的星光强度就可以得到穿过整层的垂直透过率或对应恒星方向的斜程大气透过率：

垂直透过率：                          (5)

目标方向的斜程大气透过率：                                 (6)

3．据恒星测量结果推算其附近空间目标处的整层大气透过率：

对天体目标的辐射测量进行线性拟合，为太阳或恒星的天顶角，为整层垂直大气光学厚度，则天顶角为的被测空间目标对应的斜程大气光学厚度为，相应地斜程大气透过率为。

研制该仪器的目的是为了计算地面沿卫星方向的大气透过率，为了计算，采取的方法如下。

天顶方向的大气透过率为，光学厚度为，

测量目标（太阳、恒星）的大气透过率为光学厚度

待测目标的大气透过率为,光学厚度

则：，   ∴

     =                       

4．白天测量

设计自动跟踪太阳系统，设计光学镜筒和后接光学系统、截止滤光片、光栅、线阵CCD及驱动系统、数据采集及处理；光栅前放置截止滤光片以消除光栅在宽光谱范围内的一级和二级（或多级）光谱重叠现象。采用宽光谱范围的平场凹面光栅（光路图如图1所示）。CCD采用宽响应范围，高灵敏度的线阵CCD。

光栅采用平场凹面光栅，其优点是成像光场平场，减小CCD测量的误差，且相当于平面光栅和会聚透镜的作用。选用法国JY公司的平场凹面光栅，波长范围300nm－1170nm, 根据光栅参数可计算[6]：

线色散：DL=28μm/nm(f=210mm);谱线的半角宽度： =2.9nm;

LA                                 入射狭缝

                                            N

     光栅                          LB1                                    

                              LB2                λ1

                                              线阵CCD

                                            λ2

                                    H

                          图1，太阳连续光谱测量光路图

线阵CCD采用美国PerkinElmer公司生产的专门用于光谱测量的R1210LGQ,像元2.5mmХ25μm峰值相应在650nm处，光谱范围300nm－1000nm,动态范围达到70000。仪器定标后测量出以下光谱辐射值如图2所示，通过计算对应的透过率如图3所示。

在可见到近红外波段，主要吸收气体是氧气和水汽。水汽主要集中在对流层下部，其含量与气温及天气形势有密切关系。在0.83μm和0.94μm附近的吸收峰主要是由水汽造成的；而均匀混合气体在可见光波段的吸收主要是氧气造成的，在0.68μm、0.73μm和0.76μm附近就是氧气的吸收峰。可见光波段的0.685μm至0.705μm，0.73μm至0.755μm，红外波段的0.83μm至0.μm，0.99μm至1.05μm，大气斜程透过率数值出现极大值，是比较理想的大气窗口。

图2，实际测量的太阳辐射值                        图3，实时大气透过率曲线

4.晚上测量

不同于白天太阳的测量,晚上以恒星为测量对象。测量恒星辐射有两个特点:一是供选择测量的目标多,可根据所选择的星表选择某一天区的0等到6等恒星作为被测目标,这样测量该天区的大气透过率值就不是通过换算得来的,而是实测得到的,准确度较高；另一个特点是恒星属于弱辐射,不同的星等辐射亮度不一样,因此不能象白天太阳辐射测量那样通过分光谱连续测量,实际测量方法是通过望远镜接受恒星辐射,在目镜处加多个宽带滤光片,光聚焦后入射到增强CCD靶面上,用CCD直观探测，这样可以消除因望远镜视场小而有选择性的去掉视场内的其他恒星的干扰辐射,如果用光电倍增管来采集信号不可能消除其它恒星的影响.

望远镜选用美国MEADE公司14″史密特－卡塞格林折反式天文望远镜14″LX200GPS-SMT（实物图见图4），既有大口径采光特点又有反射后折射到焦点后成像的高质量和高分辨率，主镜直径370mm，净口径356mm,焦长3556mm,焦比f10，分辨率0.32弧秒，极限目视星等18.5，望远镜视场： 0.65°;能自动驱动望远镜到数据库（145000个太空目标）中的任何一个太空，能GPS自动定位（经度、纬度、时间信息）和实现全自动跟踪；可观测陆地和天空任何目标;PC电脑通过RS232接口控制望远镜;微聚焦器调整目镜焦距。带8X50导星镜。跟踪时可先用它来粗定目标;运行速度分九等可调;可生成漫游表实现多目标跟踪，可增加观测点到用户定义观测点列表中。

ICCD的选取: 考虑望远镜接收面直径为356mm,设入射光λ=700nm，焦距f=3560mm, 衍射极限半径R=1.22λf/D=8μm.考虑3倍衍射极限时光斑直径约为50μm。考虑望远镜次镜影响，估计有效接收光面为320mm,在望远镜焦面附近光通亮放大倍数估计为(320/0.05)2=4x107倍，考虑大气吸收、镜面效率、接收面偏离焦面等其它影响也有107的放大倍数，望远镜接收光通亮为3Х10-13W/m2的入射光在接收面处估计为3X10-6W/m2,在555nm处换算为Lx,大约为2X10-3Lx,当滤光片带宽选50nm时，这个值会稍高一些。这是选择CCD最低照度的主要依据，根据计算选用北方夜视公司生产的ICCD—1XZ18/18WS-7。可以通过电压调整增益值,也可通过脉冲宽度控制曝光时间，试验中控制增益是用XICOR公司可编程电阻器X9C103，它是数字控制的微调电阻器，包含有99个电阻单元的电阻阵列，用它调整增益控制电压值，每档精度为0.25V,通过生产厂家定标得到每一控制电压值所对应的增益值。试验仪器外形图如图5所示。

星表选择依巴谷星表，该卫星是欧洲航天局于19年8月8日发射升空，专门用于天文测量的卫星，所测恒星数据受大气影响很小，在星表中记录了全天118 218个天体非常精确的赤经、赤纬、光谱等数据，其位置中值精度大约为1毫角秒，因此以依巴谷星表作为数据源，奠定了恒星位置精度的基础。

图4，望远镜实物图                         图5，试验系统外形图

恒星跟踪：

光学系统选取适当数量和波长的宽带滤光片组，考虑恒星辐射特点滤光片采中心波长分别为300、350、400、450、500、550、600、650、700、750nm，带宽为50nm可基本覆盖整个可见光波段。滤光片由单片机控制定位和自动快速切换。通过计算可以得到大气透过率。

     测量结果：

参考文献

1  路远，凌永顺，红外辐射大气透射比的简易计算，红外技术，2003，25（5）：45~49

2  Kambezidis H D,Vera D-P,Adamopoulos A D,Radiative transfer.Ⅰ.Atmospheric transmission monitoring  with modeling and  ground-based  multispectral  measurements, App.Opt.,

1997,36(27): 6976 ~  6982

3  范伟. 基于大气背景光学性质的目标探测波段选择方法研究:[博士论文].中国科学院安徽光学精密机械研究所，2005，42~49

4  詹杰，谭锟，便携式自动太阳辐射计，量子电子学报，2001，Vol.18,：551～555

5  詹杰，谭锟，新型多功能太阳辐射计用于大气光学厚度测量，安徽气象，1999：76～79

6  赵凯华，钟锡华，光学，下册，第1 版，北京：北京大学出版社，2001：18～21下载本文