地气光辐射对空间目标成像特性影响分析

2021-09-11张路平陈盛胡谋法

北京航空航天大学学报 2021年8期

张路平陈盛胡谋法

(国防科技大学电子科学学院, 长沙 410073)

天基光电测量技术在航天器捕获追踪、交汇测量等方面有着广泛的应用。空间目标探测过程中,光电探测系统除了接收源自于太阳辐照的目标辐射能量外,还会接收地气光辐射。地气光辐射是指以地球表面为下限、大气上界为上限的地球-大气系统向宇宙空间释放的各类辐射的总称,主要以反射的太阳光为主[1]。以地气光辐射为背景的地球临边背景噪声严重影响空间弱小目标的成像,研究如何消减地气光辐射对空间目标成像影响成为当前空间弱小目标探测识别亟需解决的问题。文献[2]建立了简化的航天器地气光分布计算模型,分析了地气光在航天器视场中的分布、随季节更替的变化。文献[3]建立了地表反射光模型,提出将地气光作为辅助照明光源照射空间目标的表面照度计算方法。文献[4]分析了深空背景影响下,空间目标可见光相机探测的原理。文献[5]建立了简易的地气光对星敏感器辐照度影响模型。文献[6]提出了一种利用姿态调整调节机械装置规避地气光的空间碎片探测方法。上述方法均未系统地分析地气光辐射对空间目标成像的影响,未考虑地气光辐射无法通过调节机械装置进行完全抑制时,如何通过软件算法消减其对目标成像的影响。分析地气光辐射对空间目标成像特性的影响成为开展地气光背景抑制研究的基础。

本文以地球同步轨道(GEO)卫星上的天基空间光电探测系统为探测平台,将空间目标及地气光辐射放到同一场景,采用微元法进行辐射建模,通过卫星工具包(STK)设计近地轨道(LEO)及高椭圆轨道(HEO)2 个不同轨道上的目标运动场景,分析角度及距离参数变化对不同轨道空间目标及地气光背景等效星等的影响,以及不同观测时刻地气光辐射对空间目标信噪比的影响,为进行空间目标探测识别提供理论计算依据,同时也为减小、规避地气光辐射影响,在空间目标探测时间段的选择上提供时间维度参考。

1 空间目标成像及地气光辐射建模

1.1 空间目标成像可见性分析

在分析地气光辐射对空间目标成像特性的影响时,首先,考虑太阳、地球及空间探测平台三者之间的空间位置关系,确保空间探测平台与空间目标之间的视线不应该受到地球的遮挡,如图1所示。目标1 位于探测视场内,目标2 位于地球遮挡临界点,目标探测会受到影响,目标3 位于地球阴影区,完全被遮挡。其次,要求探测器探测空间目标时避开太阳直射,否则会使视场处于高亮背景,对探测器造成损坏。最后,从能量角度要求空间目标辐射信号强度满足探测器的最低探测阈值,确保目标可被探测[7]。

图1 地球遮挡示意图Fig.1 Schematic diagram of earth occlusion

1.2 空间目标辐射建模

将可见光探测器看成一个积分光子计数器,空间目标的辐射能量用光子数来表示。空间目标辐射能量主要来自于太阳光辐照。根据普朗克方程可计算出太阳光在可见光光谱范围内的辐射强度Is(单位W/sr)为

式中:α为太阳光与空间目标、探测器之间形成的探测相角;ρt为空间目标的光谱反射率;At为空间目标的等效截面积。

将探测器采集到的空间目标信号能量转化为光子数,得到进入探测器的空间目标信号光子流量密度Φt(单位photons·m-2·s-1)为

式中:Ad为探测器光学系统的有效入瞳面积;τ为光学系统透过率;Q为探测器对信号光的量子效率;t为信号曝光时间;Np为空间目标成像像素数。

1.3 地气光辐射模型

地气光辐射对空间目标成像的影响是通过增加地气光背景噪声从而改变空间目标信噪比实现的[8]。由于地球自身辐射很小,地气光辐射主要考虑来自于地球反射的太阳光。假定太阳到地球大气层表面的距离为Rse,太阳对地球大气层表面的辐照度Ese(单位W/m2)为

将地球看作是一个朗伯球体,探测器对轨道空间目标及地球表面地气光辐射的观测光路如图2所示。图中:β为太阳到大气层表面面元的连线与面元区域法线方向的夹角,即光照角,θ为面光源对探测器投射光方向与面元法向的夹角,即出射角,ζ为面光源对探测器投射光方向与探测器光轴方向的夹角,即离轴角。

图2 空间目标与地气光辐射观测示意图Fig.2 Schematic diagram of space object and earth-atmosphere radiation observation

将地球表面看成是一个等效朗伯球体,当受到太阳光照射时,由于漫反射,该等效朗伯球体表面又将反射的太阳光辐射到达探测器,为便于分析,只考虑能够进入探测器视场的地气光影响。

采用微元法分析地气光辐射模型。太阳光照射到地球表面的某一区域内,将该区域分成尽可能多的微小面元,该微小面元作为一个面光源产生的辐射强度为Ie(单位W/sr),若地球表面的反射率为ρe,则Ie可表示为

该微小面元又将太阳光辐射到探测器。假定地气光辐射进入探测器遮光罩入口平面的地球表面区域总面积为Ae,到探测器之间的距离为Red,则在探测器表面产生的地气光辐照度Eed(单位为W/m2)

式中:E0=2.96 ×10-8W/m2为0 星等空间目标在地球大气层外产生的辐照度。

地气光背景噪声在探测器上产生的光子数Ne为

式中:Φe=5 ×1010/2.512me为单位立体角地气光背景噪声光子流量密度;Ap为探测器单个像元的角面积。

将式(9)、式(10)代入式(11)可得地气光背景噪声光子数为

1.4 空间目标信噪比分析

式(5)、式(12)、式(15)为计算和分析地气光辐射影响下的空间目标信噪比的理论基础。

2 地气光辐射对空间目标信噪比的影响

2.1 STK 场景设计

为分析地气光辐射对空间目标信噪比的影响,通过STK 对天基探测平台与卫星目标进行场景设计[15]。探测平台为GEO 卫星上搭载的可见光探测器,探测目标为逆时针运动的LEO 卫星与HEO 卫星。假设卫星目标的等效球体半径为2 m,空间目标的光谱反射率ρt为0.33。部分探测器参数参照了美国“空间中段监视”试验卫星上搭载的天基可见光相机(SBV)参数及文献[4]中仿真设计的探测器参数,详细参数设置如表1所示。

表1 探测系统参数Table 1 Parameters of observation system

通过STK 场景数据报告可得到GEO、HEO、LEO 轨道卫星及太阳在J2000 地球惯性坐标系下的实时空间坐标与运动速度。 STK 软件仿真场景如图3 所示。

图3 STK 场景Fig.3 STK scene

2.2 地气光辐射对LEO 及HEO 轨道目标信噪比的影响

在GEO 探测平台布置2 个探测器,采用主动方式追踪LEO 轨道目标与HEO 轨道目标(即目标离轴角ζ=0°)。为分析地气光辐射对空间目标信噪比的影响,根据太阳、地球、不同轨道上卫星及探测平台之间的位置关系选择探测时间段,确保地气光辐射进入探测器目标视场内,即空间目标及地气光辐射均进入探测器视场。以地气光进入和离开探测器目标视场的时刻为观测开始及结束时段,称该时间段内的观测为一次完整观测。

本文仅分析一个GEO 轨道周期内(23.93 h)地气光辐射对空间目标信噪比的影响。 HEO 和LEO 轨道周期分别为11. 96 h 和1. 64 h,根据GEO 特性,在一个GEO 轨道周期内,理论上HEO轨道目标至多有1 次完整观测,LEO 轨道目标至多有7 次完整观测。场景仿真时选择的观测日期为2019 年3 月21 日,根据目标可观测条件及地气光辐射进入探测器目标视场时间的综合要求,HEO 轨道目标有1 次完整观测,LEO 轨道目标有4 次完整观测,表2、表3 分别给出了2 个轨道目标的完整观测时间段。

表2 地气光辐射进入HEO 轨道目标探测视场时间段Table 2 Period when earth-atmosphere radiation enters the HEO orbit target detection field of view

表3 地气光辐射进入LEO 轨道目标探测视场时间段Table 3 Period when earth-atmosphere radiation enters the LEO orbit target detection field of view

2.2.1 HEO 及LEO 轨道目标的第1 次完整观测

影响空间目标信噪比的因素包括与空间目标光子数相关的α、Rtd、Rst,以及与地气光光子数相关的β、θ、Red、Rse。表4 分析了第1 次完整观测时间段内上述参数的统计值。

表4 距离和角度参数变化Table 4 Variation of distance and angle parameters

采样间隔为1 s,第1 次完整观测时间段内空间目标等效星等及相关参数变化如图4 所示。

综上所述,微信群护理人员培训手段的应用,是近几年来临床上新兴的一种护理理论培训措施,其主要是通过对需要培训的护理人员在工作后的业余时间,配合完成各项护理理论知识的培训,确保护理工作有序开展,也是一种培养优秀护理人才的好的培训方法。同时在培训过程中不断总结经验,探讨改进措施,完善培训方法,努力提高培训质量。

图4 空间目标等效星等及相关参数变化Fig.4 Variation of space target equivalent magnitude and related parameter

第1 次完整观测时间段内,HEO 轨道目标由远地点向近地点运动,远离探测器及太阳,Rst与Rtd单调递增,太阳、探测平台及空间目标之间形成的探测相角α单调递减,变化范围为49° ～55°;LEO 轨道目标从地球北极上空经赤道向南极上空运动,Rst与Rtd先减后增,α单调递减,取值范围约为25° ～34°。

根据地气光辐射模型得到观测时段内地气光背景等效星等及相关参数变化,如图5 所示。

图5 地气光背景等效星等及相关参数变化Fig.5 Variation of earth-atmosphere background equivalent magnitude and related parameter

HEO 及LEO 轨道目标从靠近北极上空经赤道再向南极运动过程中,地气光辐射进入空间目标探测视场,Rse与Red先减后增,地气光光照角β、出射角θ先减后增。

得到空间目标及地气光背景等效星等后,根据空间目标信噪比模型,可得到第1 次完整观测时间段内LEO 及HEO 轨道目标的信噪比变化,如图6 所示。

图6 空间目标信噪比变化Fig.6 Variation of space target SNR

在地气光辐射进入和离开空间目标探测视场时,空间目标信噪比变化较大,HEO 轨道目标信噪比在前200 s 观测内从7.76 下降到1.75,随后稳定在1.5,地气光辐射离开空间目标探测视场前200 s 内,空间目标信噪比由1.47 上升到6.94;LEO 轨道目标信噪比在地气光出入视场时刻变化幅度相对HEO 轨道目标要小。

2.2.2 LEO 轨道目标的4 次完整观测

LEO 轨道目标的4 次完整观测时段内,空间目标与地气光背景等效星等及信噪比变化如图7所示。

图7 不同观测时段空间目标及地气光背景特性变化Fig.7 Characteristic variation of target and earth-atmosphere background in different observation periods

LEO 轨道目标等效星等均呈现先减小后增大趋势。地气光辐射第3 次进入LEO 轨道目标探测视场后的200 s 内,空间目标信噪比从13.19下降到2.13,随后稳定在2.0 左右,地气光辐射离开空间目标探测视场的前200 s 内,LEO 轨道目标信噪比从2.05 升到19.11,整个观测时段内信噪比变化幅度约为18。

地气光背景等效星等在地气光进出空间目标探测视场时刻变化较大,其等效星等比空间目标等效星等低,即地气光背景噪声强度比空间目标信号强度大,空间目标信噪比主要受地气光辐射影响。地气光辐射越强,空间目标信噪比越低,反之越高。在不同探测时段,得到的空间目标及地气光背景等效星等相差不大。

2.3 模型验证

2.3.1 空间目标等效星等验证

陈荣利等[9]利用普森公式将空间目标照度转化为星等数来描述,其球体模型如下:

式中:mt为空间目标等效星等数;Rtd以100 km 为单位;D为空间目标等效直径;σ为位相角。

利用该模型计算得到的结果与紫金山天文台及国家天文台兴隆站观测的实测数据基本吻合,仅存在0.5 个星等的差异。经分析,其星等误差是由仿真模型表面光学特性假设、仿真模型形状假设与真实空间目标之间的差异引起的。虽然本文未能通过实测数据验证空间目标等效星等模型的正确性与合理性,但可通过与基于普森公式构建的模型进行对比,从而间接验证本文模型的正确性与合理性。

将本文建模所用的空间目标反射率ρt、空间目标等效球体直径(4 m)、空间目标探测距离Rtd、位相角σ等参数代入普森公式,计算出基于该模型的空间目标等效星等,同时与本文模型得到的空间目标等效星等进行对比,对比结果如图8及表5 所示。

表5 两种模型空间目标等效星等数据对比Table 5 Comparison of space target equivalent magnitude between two models

图8 两种模型得到的空间目标等效星等变化Fig.8 Equivalent magnitude variation of space objects with two models

图8 及表5 对比结果表明,由于空间目标探测距离与位相角参数取值的有效位数误差,本文模型计算得到的HEO 轨道目标与LEO 轨道目标等效星等与基于普森公式计算得到的等效星等仅有0.1 个星等的差异。 2 种模型计算结果表明,本文模型计算得到的空间目标等效星等与真实观测数据星等吻合,构建的空间目标等效星等模型正确合理。

2.3.2 地气光辐射模型验证

原育凯等[2]通过建立坐标系,对不同季节的地气辐射变化进行了分析。肖相国等[3]将地气光辐射作为空间目标的辅助照明光源进行研究。这些研究成果都涉及地气光辐射建模问题,本文借鉴上述成果研究思路,采用微元法对地气光的辐射进行建模分析,从而得到地气光背景等效星等计算方式,并进一步推导出空间目标在探测器中的成像信噪比变化,上述研究成果为本文地气光辐射模型构建提供了理论支撑,确保本文构建模型的正确性及合理性。