留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Thermal radiation stray light integration method of infrared camera in geostationary orbit

Lu Qiang

陆强. 地球同步轨道红外相机的热辐射杂散光集成法[J]. 红外与激光工程, 2020, 49(5): 20190457. doi: 10.3788/IRLA20190457
引用本文: 陆强. 地球同步轨道红外相机的热辐射杂散光集成法[J]. 红外与激光工程, 2020, 49(5): 20190457. doi: 10.3788/IRLA20190457
Lu Qiang. Thermal radiation stray light integration method of infrared camera in geostationary orbit[J]. Infrared and Laser Engineering, 2020, 49(5): 20190457. doi: 10.3788/IRLA20190457
Citation: Lu Qiang. Thermal radiation stray light integration method of infrared camera in geostationary orbit[J]. Infrared and Laser Engineering, 2020, 49(5): 20190457. doi: 10.3788/IRLA20190457

地球同步轨道红外相机的热辐射杂散光集成法

doi: 10.3788/IRLA20190457
详细信息
  • 中图分类号: O439;TN215

Thermal radiation stray light integration method of infrared camera in geostationary orbit

  • 摘要: 地球同步轨道上的三轴稳定卫星所处的空间环境复杂,仪器背景辐射变化较大。传统杂散光分析法无法模拟非均匀温度场,且无法实时计算仪器背景,仿真误差较大。提出用热辐射杂散光集成法来分析在轨红外相机的仪器背景辐射,通过更趋于在轨真实温度的且具有温度梯度的实时温度场,结合辐射传递因子,计算探测器上的仪器背景辐射以及相机的信杂比。将热辐射杂散光集成法、传统杂光分析法计算信杂比与在轨实测信杂比进行对比,热辐射杂散光集成法误差小于17%,而传统杂光分析法误差达114%。表明热辐射杂散光集成法的仿真结果更趋近于在轨实际情况,仿真效率和仿真精度更高。
  • Figure  1.  Schematic diagram of light path

    Figure  2.  Calculating process of instrument background radiation stray light for TRSLIM

    Figure  3.  Thermal network model

    Figure  4.  External heat flow into the hood

    Figure  5.  Temperature curve of primary mirror in one day

    Figure  6.  Temperature gradient curve of primary mirror in one day

    Figure  7.  Radiation curve of instrument background radiation produced by primary mirror on detector in one day

    Figure  8.  Total instrument background radiation curve generated by optical-mechanical components on detector in one day

    Figure  9.  SCR curve formed by instrument background radiation in one day

    Figure  10.  Variation curve of SCR at band 3 in one day

    Table  1.   Radiation calculation of primary mirror at different time in infrared band

    TimeTemperature/°CEmissivityM/(W/m2)
    12:00−1.30.040.049 0
    20:00−5.60.040.040 5
    6:002.30.040.057 2
    下载: 导出CSV

    Table  2.   Energy transfer factor Bij of nodes of primary mirror to detector

    Cond_idNode 1Node 2Area*e*BijBijBji
    1COMET.9284COMET.67902.728 026e-050.018 5320.012 89
    2COMET.9284COMET.67915.640 484e-070.003 455 9210.002 322 2
    3COMET.9284COMET.73081.156 538e-060.013 559 960.001 076
    4COMET.9284COMET.73091.401 487e-060.007 363 9610.003 127 9
    ……
    N-3COMET.9310COMET.73096.087 829E-085.18E-050.001 187
    N-2COMET.9310COMET.73101.103 783E-060.009 397 840.004 149
    N-1COMET.9310COMET.73116.306 837E-070.005 369 770.001 819
    NCOMET.9310COMET.73122.469 422E-060.002 102 5180.007 124
    下载: 导出CSV

    Table  3.   Radiation calculation of nodes of primary mirror on detector at 12:00

    Cond_idNode 1Node 2BijEj/(W/m2)E/(W/m2)
    1COMET.9284COMET.67790.005 853 20.000 7263.712E-04
    2COMET.9284COMET.67800.000 345 592 10.000 135
    3COMET.9284COMET.72700.003 559 960.000 532
    4COMET.9284COMET.72710.001 363 9610.000 289
    ……
    N-3COMET.9310COMET.72715.18E-052.03E-06
    N-2COMET.9310COMET.72720.009 397 840.000 368
    N-1COMET.9310COMET.72730.000 536 9770.000 21
    NCOMET.9310COMET.72740.002 102 5188.24E-05
    下载: 导出CSV

    Table  4.   Radiation calculation of nodes of primary mirror on detector at 20:00

    Cond_idNode 1Node 2BijEj/(W/m2)E/(W/m2)
    1COMET.9284COMET.67790.005 853 20.000 5951.820E-04
    2COMET.9284COMET.67800.000 345 592 10.000 111
    3COMET.9284COMET.72700.003 559 960.000 435
    4COMET.9284COMET.72710.001 363 9610.000 236
    ……
    N-3COMET.9310COMET.72715.18E-051.66E-06
    N-2COMET.9310COMET.72720.009 397 840.000 302
    N-1COMET.9310COMET.72730.000 536 9770.000 172
    NCOMET.9310COMET.72740.002 102 5186.75E-05
    下载: 导出CSV

    Table  5.   Radiation calculation of nodes of primary mirror on detector at 6:00

    Cond_idNode 1Node 2BijEj//(W/m2)E/(W/m2)
    1COMET.9284COMET.67790.005 853 20.000 8545.855E-04
    2COMET.9284COMET.67800.000 345 592 10.000 159
    3COMET.9284COMET.72700.003 559 960.000 625
    4COMET.9284COMET.72710.001 363 9610.000 339
    ……
    N-3COMET.9310COMET.72715.18E-052.39E-06
    N-2COMET.9310COMET.72720.009 397 840.000 433
    N-1COMET.72730.000 536 9770.000 248
    NCOMET.9310COMET.72740.002 102 5189.69E-05
    下载: 导出CSV

    Table  6.   Comparison of SCR from on-orbit measured, TRSLIM and TSLAM

    TimeOn-orbit test valuesTRSLIMTSLAM
    Simulation valueErrors from on-orbit test valuesSimulation valueErrors from on-orbit test values
    3:000.0590.061.69%0.11289.83%
    8:000.0640.07415.63%0.13103.13%
    12:000.0780.09116.67%0.14484.62%
    15:000.0890.10315.73%0.189112.36%
    18:000.0950.10914.74%0.204114.74%
    22:000.0890.09810.11%0.179101.12%
    下载: 导出CSV
  • [1] Zhao Fei, Wang Sen. Study of stray light for the Xinglong 1-meter optical telescope [J]. Astronomical Research & Technology, 2010, 7(2): 158−167.
    [2] Huang Qiang, Niu Xinhua, Shen Xuemin. Stray radiation analysis caused by interior heat radiation in infrared optical systems [J]. Infrared Technology, 2006, 28(6): 348−352.
    [3] Niu Jinxing, Zhou Renkui, Liu Zhaohui. Analysis of stray light caused by thermal radiation of infrared detection system [J]. Acta Optica Sinica, 2010, 30(8): 2267−2271. doi:  10.3788/AOS20103008.2267
    [4] Pei Linlin, Xiang Libin, Lv Qunbo. Research on the stray light suppression for super-resolution optical system [J]. Acta Photonica Sinca, 2017, 46(11): 1122002. doi:  10.3788/gzxb20174611.1122002
    [5] Wang Weili, Lyu Yong, Li Xiaoying, et al. Stray light analysis on infrared channel of dual-band imaging system [J]. Journal of Applied Optics, 2018, 39(2): 262−267.
    [6] Lu Haiping, Yuan Xiangyan, Zhang Kaiyuan. Infrared background radiation removing design and simulation of Antarctic survey telescope [J]. Acta Optica Sinica, 2014, 34(11): 1122002. doi:  10.3788/AOS201434.1122002
    [7] Zhao Yuchen, Xu Yanjun, Sha Wei. Stray light analyze and suppress of the Space-Borne infrared optical system [J]. Chinese Journal of Lasers, 2015, 42(10): 1016001. doi:  10.3788/CJL201542.1016001
    [8] Yang Jiaqiang, Li Ronggang, Peng Qingqing. Research of stray light in the IR system based on Lighttools and Matlab [J]. Laser & Infrared, 2014, 44(8): 888−891.
    [9] Zurmehly G. Thermal and structural analysis of the GOES scan mirror’s on-orbit performance[C]//SPIE, 1991, 1532: 170-176.
    [10] Lightsey Paul A, Wei Zongying. Stray light performance for the James Webb Space Telescope[C]//SPIE, 2014, 9143: 91433P.
    [11] Gerhard S, David C, Gerard P. Advancements in Integrated Structural/Thermal/Optical (STOP) analysis of optical systems[C]//SPIE, 2007, 6675: 66750D.
    [12] Stoeckel Gerhard P, Doyle Keith B. Integrated optomechanical analysis and testing software development at MIT Lincoln Laboratory[C]// SPIE, 2013,8840:88400E.
    [13] Min Guirong. Satellite Thermal Control Technology [M]. Beijing: Aerospace Press, 1991.
  • [1] 厉夫兵, 游啟, 冷俊敏, 杨林昊.  漫反射与非漫反射表面间的辐射传递系数快速计算 . 红外与激光工程, 2024, 53(3): 20230611-1-20230611-11. doi: 10.3788/IRLA20230611
    [2] Fu Yanjun, Cai Xiaoqi, Zhong Kejun, Ma Baiheng, Yan Zhanjun.  Method for phase-height mapping calibration based on fringe projection profilometry . 红外与激光工程, 2022, 51(4): 20210403-1-20210403-9. doi: 10.3788/IRLA20210403
    [3] Fan Yu, Yuan Qianqian, Jiang Haitao.  Fabrication of low Mg content MgxZn1-xO nanowires ultraviolet photosensors via chemical vapour deposition method . 红外与激光工程, 2021, 50(9): 20200448-1-20200448-7. doi: 10.3788/IRLA20200448
    [4] Yang Sen.  Hybrid PSO-AMLS-based method for data fitting in the calibration of the infrared radiometer . 红外与激光工程, 2021, 50(8): 20200471-1-20200471-9. doi: 10.3788/IRLA20200471
    [5] Liang Yan, Zhang Qingdong, Zhao Ning, Li Chuanmiao.  Indoor location method based on UWB and inertial navigation fusion . 红外与激光工程, 2021, 50(9): 20200484-1-20200484-14. doi: 10.3788/IRLA20200484
    [6] Li Ning, Wei Daozhi, Zhang Dongyang, Yao Liangfu.  Information fusion method of ranging-imaging guidance integrated fuze . 红外与激光工程, 2021, 50(11): 20210039-1-20210039-8. doi: 10.3788/IRLA20210039
    [7] Xu Baoxiang, Xiong Zhi, Huang Jixun, Yu Haicheng.  Research on method of improving magnetic field adaptability of high-precision IFOG . 红外与激光工程, 2021, 50(4): 20200239-1-20200239-8. doi: 10.3788/IRLA20200239
    [8] Guo Ming, Zhang Yongxiang, Zhang Wenying, Li Hong.  Thermal damage of monocrystalline silicon irradiated by long pulse laser . 红外与激光工程, 2020, 49(3): 0305002-0305002-9. doi: 10.3788/IRLA202049.0305002
    [9] Guo Ming, Zhang Yongxiang, Zhang Wenying, Li Hong.  Thermal damage of monocrystalline silicon irradiated by long pulse laser . 红外与激光工程, 2020, 49(3): 1-9.
    [10] Yang Xu, Jiang Pengfei, Wu Long, Xu Lu, Zhang Jianlong, Hu Haili, Liu Yuehao, Zhang Yong.  Underwater Fourier single pixel imaging based on water degradation function compensation method . 红外与激光工程, 2020, 49(11): 20200281-1-20200281-12. doi: 10.3788/IRLA20200281
    [11] Niu Chaojun, Wang Xiaobin, Lu Fang, Han Xiang’e.  Validity of beam propagation characteristics through oceanic turbulence simulated by phase screen method . 红外与激光工程, 2020, 49(7): 20190452-1-20190452-7. doi: 10.3788/IRLA20190452
    [12] Liu Lu, Xi Dongdong, Cheng Lei, Wang Yuwei, Cai Bolin, Zhou Huiyu.  Enhanced Gray-code method for three-dimensional shape measurement . 红外与激光工程, 2020, 49(11): 20200314-1-20200314-8. doi: 10.3788/IRLA20200314
    [13] Fang Zeyuan, Yin Lu, Yan Mingjian, Han Zhigang, Shen Hua, Zhu Rihong.  Study on signal light transmission efficiency enhancement of backward pump-signal combiners in high-power fiber lasers . 红外与激光工程, 2020, 49(10): 20200014-1-20200014-10. doi: 10.3788/IRLA20200014
    [14] 曾瀚林, 孟祥勇, 钱惟贤.  高斯差分滤波图像融合方法 . 红外与激光工程, 2020, 49(S1): 20200091-20200091. doi: 10.3788/IRLA20200091
    [15] 陈鹏, 赵继广, 宋一铄, 王燊.  气溶胶环境下FMCW与脉冲激光探测性能对比 . 红外与激光工程, 2020, 49(6): 20190399-1-20190399-10. doi: 10.3788/IRLA20190399
    [16] Chengmingyue Li.  Optical phase conjugation (OPC) for focusing light through/inside biological tissue . 红外与激光工程, 2019, 48(7): 702001-0702001(9). doi: 10.3788/IRLA201948.0702001
    [17] Ji Yi.  Visible light optical coherence tomography in biomedical imaging . 红外与激光工程, 2019, 48(9): 902001-0902001(9). doi: 10.3788/IRLA201948.0902001
    [18] Jitao Zhang, Mingguang Tuo, Min Liang, Wei-Ren Ng, Michael E. Gehm, Hao Xin.  Terahertz radiation of a butterfly-shaped photoconductive antenna(invited) . 红外与激光工程, 2019, 48(4): 402001-0402001(9). doi: 10.3788/IRLA201948.0402001
  • 加载中
计量
  • 文章访问数:  985
  • HTML全文浏览量:  640
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-12
  • 修回日期:  2020-01-15
  • 网络出版日期:  2020-05-08
  • 刊出日期:  2020-05-27

Thermal radiation stray light integration method of infrared camera in geostationary orbit

doi: 10.3788/IRLA20190457
  • 中图分类号: O439;TN215

摘要: 地球同步轨道上的三轴稳定卫星所处的空间环境复杂,仪器背景辐射变化较大。传统杂散光分析法无法模拟非均匀温度场,且无法实时计算仪器背景,仿真误差较大。提出用热辐射杂散光集成法来分析在轨红外相机的仪器背景辐射,通过更趋于在轨真实温度的且具有温度梯度的实时温度场,结合辐射传递因子,计算探测器上的仪器背景辐射以及相机的信杂比。将热辐射杂散光集成法、传统杂光分析法计算信杂比与在轨实测信杂比进行对比,热辐射杂散光集成法误差小于17%,而传统杂光分析法误差达114%。表明热辐射杂散光集成法的仿真结果更趋近于在轨实际情况,仿真效率和仿真精度更高。

English Abstract

    • For the three-axis stabilized satellite in geosynchronous orbit, its space environment temperature is complex. The satellite temperature is affected by the sun, which not only changes with the season, but also varies greatly at different time of the same day. The temperature inside the camera determines the background radiation of the instrument and affects the dynamic range of the camera. Therefore, the simulation and calculation of the background radiation of the instrument becomes an important factor affecting the accuracy and quantification of the infrared camera in geosynchronous orbit.

      For the analysis of stray light in instrument background, the traditional stray light analysis method (TSLAM) is based on Tracepro, Lighttools and other software, which greatly simplifies the temperature field of each component. Usually, the radiation analysis of uniform temperature field can not simulate the actual non-uniform temperature field in orbit, nor can the real-time calculation of the temperature field changing in one day being carried out[1-8], The simulation error is usually greater than 100%. However, because of the non-uniformity of the distribution of the external light source and the internal heat source, the temperature field of the camera components varies with time. Especially for geosynchronous orbit space camera, the temperature of each optical component varies dramatically with time in one day. The pointing mirror of the Geostationary Operational Environment Satellite system (GOES) is exposed to the outside, the maximum temperature difference can reach 7 K, and the temperature gradient on the surface of GOES’s pointing mirror will change with the sun's relative orientation[9]. The temperature difference in different regions of the primary mirror of the James Webb Space Telescope (JWST) is as high as 14 K. If the average temperature field is used in stray light simulation, the error is large. Therefore, the method of subregional simulation is adopted[10].

      It can be seen that the temperature field of geosynchronous orbit varies greatly with time. The background radiation of the instrument caused by the temperature field has become an important factor affecting the imaging quality of infrared cameras. Real-time analysis of the influence of stray light from instrument background radiation on camera detection ability caused by non-uniform temperature field in one day is a necessary step in optical system design. In this paper, the thermal radiation stray light integration method (TRSLIM) is proposed to analyze the instrument background radiation of on-orbit infrared camera. The temperature field of each optical component is directly simulated by thermal radiation stray light software. The temperature field is real-time and has temperature gradient, which tends to be more realistic. The background radiation of the instrument on the detector can be calculated by combining the temperature field of each node of the opto-mechanical element with the radiation transfer factor. The signal-to-clutter ratio (SCR) of the camera can be calculated by combining the radiation characteristics of the target.

    • The optical model of the space camera is an off-axis three-mirror structure. The optical path diagram is shown in Fig.1. After the parallel light is reflected by the pointing mirror, the primary mirror, the secondary mirror and the third mirror are illuminated successively through the aperture stop. The field stop is set at the convergence point between the secondary mirror and the third mirror. The back light path adopts the way of splitting by splitters, and achieves the splitting of channel 1, channel 2 and channel 3 by setting splitter 1 and splitter 2.

      Figure 1.  Schematic diagram of light path

      For channel 3 (infrared channel), the radiation flux of the point target received by the camera on the focal plane device is as follows:

      $${{{P}}_0} = {{J}} \cdot \frac{{{{\pi D}}_0^2}}{{4{{{I}}^2}}} \cdot {{{\tau }}_{{a}}} \cdot {{{\tau }}_0} \cdot {{Enc}}$$ (1)

      The irradiance of the point target received on the focal plane is as follows:

      $${{{H}}_0} = {{J}} \cdot \frac{{{{\pi D}}_0^2}}{{4{{{I}}^2}{{{a}}^2}}} \cdot {{{\tau }}_{{a}}} \cdot {{{\tau }}_0} \cdot {{Enc}}$$ (2)

      Among them, J is the radiation intensity of the point target; I is the detection distance of the point target; τa is the atmospheric transmittance in infrared band, the value is 0.65; τ0 is the optical transmittance in the channel 3 of the camera, the value is 0.48; Enc is the concentration of diffraction energy in one pixel of the channel 3 of the camera, the value is 70%; D0 is the optical aperture of the optical system, and a is the pixel size. Finally, the irradiance of the point target on the detector of channel 3 is calculated to be 1.887×10−3 W/m2.

    • In recent years, the commercial software Structural/Thermal/Optics Performance (STOP) for optical-mechanical thermal integration analysis has been gradually integrated into the engineering application[11-12]. It can complete the simulation analysis of the image quality change of the camera under the influence of temperature field change and stress release on the orbit, but it can not complete the simulation analysis of the stray light change of the camera when it works on the orbit. In terms of the nature of light, both heat and light are electromagnetic waves, which are only the differences of manifestations. Thermal simulation analysis software and stray light simulation analysis software have common features. The most prominent point is that the two analysis software adopts MCM, a relatively mature method.

      At present, TracePro and Lighttools are mainly used to analyze instrument background radiation. These software uses mature Monte-Carlo Method (MCM) for ray tracing. However, the stray light software itself can not analyze the orbital parameters of the space camera and the periodic variation of the temperature field of the instrument, so the non-real-time and uniform temperature field is often used for approximate analysis. Thermal analysis software such as Thermal Desktop, TMG and so on, also uses MCM calculation method, but can obtain real-time instrument temperature field with temperature gradient. The radiation transfer factor is introduced into the thermal software, which includes the energy exchange between the opto-mechanical elements and the focal plane, which is concerned by stray light analysis. The radiation transfer factor remains unchanged in the process of solving the energy transfer problem. Finally, the instrument background radiation can be calculated by the instrument temperature field and the radiation transfer factor.

      Therefore, it is feasible to use thermal analysis software for stray light analysis.

      The radiation transfer factor Bij is defined as the percentage of the energy directly projected from surface 1 to surface 2 to the radiation energy of surface 1. The definition can be written as:[13]

      $${{{B}}_{ij}} = \frac{{{Q_{12}}}}{{{Q_1}}}$$ (3)

      Among them, Q1 is the effective radiation of surface 1 and Q12 is the projection radiation of surface 1 to surface 2.

      If the number of volume Vi and panel Sj in the thermal radiation system is N and M, the energy equation of volume Vi is

      $$4{\rm{\sigma }}{{\rm{\kappa }}_i}{V_i}T_i^4 =\mathop \sum \nolimits_{j = 1}^N 4{\rm{\sigma }}{{\rm{\kappa }}_j}{V_j}T_j^4{B_{ji}} + \mathop \sum \nolimits_{k = 1}^M {S_k}{\rm{\sigma }}{{\rm{\varepsilon }}_k}T_k^4{B_{ki}}$$ (4)

      The energy equation of panel Si is

      $${S_i}{{\rm{\varepsilon }}_i}{\rm{\sigma }}T_i^4 = \mathop \sum \nolimits_{j = 1}^N 4{\rm{\sigma }}{{\rm{\kappa }}_j}{V_j}T_j^4{B_{ji}} + \mathop \sum \nolimits_{k = 1}^M {S_k}{\rm{\sigma }}{{\rm{\varepsilon }}_k}T_k^4{B_{ki}}$$ (5)

      The radiation transfer factor emitted by volume Vi and absorbed by panel Sj is as follows:

      $${{B}}\left( {{{{V}}_i}{{,}}{{{S}}_j}} \right) = \frac{{{N_j}\left( {{{{S}}_j}} \right)}}{{{N_i}\left( {{{{V}}_i}} \right)}}$$ (6)

      Among them, Ni(Vi) is the total number of energy beams emitted by volume Vi and Nj(Sj) is the total number of energy beams absorbed by panel Sj.

    • The method of calculating instrument background radiation for thermal radiation stray light integrated analysis is based on the orbital model and the radiation characteristics of the camera itself. The thermal analysis of the camera is carried out, and the non-uniform temperature field of each optical component of the camera is obtained. For the element of interest, the node temperature and the specific emissivity of the element with time are taken and calculated according to Planck's law of radiation. The radiation illumination of the element node with time varying in the spectrum of interest is as follows:

      $${M_{i\lambda }} = \varepsilon \frac{{2\pi h{c^2}}}{{{\lambda ^5}}}\frac{1}{{{e^{ch/\lambda k{T_i}}} - 1}}$$ (7)

      Through the thermal radiation module of the thermal radiation stray light integrated analysis software, the radiation transfer factor Bij between each element node and detector node is calculated, and then the instrument background irradiance of each element on the detector which changes with time is calculated. The calculation formula is as follows:

      $${{E}} = \frac{{\mathop \sum \nolimits_j \mathop \sum \nolimits_i \mathop \smallint \nolimits_{{\lambda _1}}^{{\lambda _2}} {M_{i\lambda }}{d_\lambda } \cdot {B_{{{ij}}}}}}{{\mathop \sum \nolimits^ j}}$$ (8)

      Then, the total instrument background irradiance on the detector is obtained by adding the instrument background irradiance of each component. According to the calculation results of the irradiance of the point target on the detector in the last section, the SCR caused by the instrument background radiation with time is calculated. The calculation process is shown in Fig.2.

      Figure 2.  Calculating process of instrument background radiation stray light for TRSLIM

    • Because of the complex thermal characteristics of space camera in orbit, it is necessary to simplify the geometric model, the external heat flow and the thermal characteristics of nodes before building the model. The simplification is mainly based on the technical status and heat transfer characteristics of the space camera. In the process of building thermal model, the partition of nodes and grids is also an important step. If the main thermal characteristics of space camera can be embodied, the number of nodes will be reduced as much as possible. If the optical-mechanical element is isothermal, it can be regarded as a node. If the temperature gradient of the optical-mechanical element is large, the node should be subdivided. In addition, the space camera is divided into several different components, including the main frame component, the pointing mechanism component, the main optical component, the rear light path component and the refrigerator Dewar component. The thermal contact relationship between the components is defined separately and connected through the grid assembly relationship. The final thermal analysis model of space camera is as shown in Fig.3:

      Figure 3.  Thermal network model

    • The orbit of the space camera is geosynchronous. Figure 4 shows the external heat flow of the camera in orbit around 4 hours of midnight under five working conditions of + 8.8 °, - 8.8 °, spring and autumn equinox, winter solstice and summer solstice. It can be seen from the picture that the external heat flow into the camera is different around midnight in different seasons. At + 8.8 °, the external heat flow into the camera is the most, at the winter solstice the least, and at the spring and autumn equinoxes there is a 72 min shadow period. Because of the most external heat flow and the most drastic change of camera temperature field in +8.8 ° working condition, we choose the solar altitude angle +8.8 ° working condition to calculate the temperature field inside the camera, and then analyses the influence of the corresponding instrument background radiation.

      Figure 4.  External heat flow into the hood

    • The TRSLIM is carried out for the infrared channel of channel 3. Firstly, the background irradiance of each optical-mechanical element should be calculated.

      The temperature change curve and temperature gradient change curve of each optical-mechanical element in one day were obtained by thermal model and orbit parameters. Take the primary mirror as an example, Fig.5 and Fig.6 are the temperature change curve and the temperature gradient curve of the primary mirror in one day, respectively. The temperature of the primary mirror at 12:00, 20:00 and 6:00 is used to calculate the radiation in the infrared band. Table 1 calculates the average irradiance of the primary mirror M in the required infrared band. Table 2 calculates the radiation transfer factor Bij between the primary mirror nodes and the infrared detector nodes by TRSLIM, in which the primary mirror nodes are Node1 and the infrared detector nodes are Node2. Based on the actual irradiance and Bij of the primary mirror nodes in the infrared band, the average background radiation of the primary mirror on the infrared detector at 12:00, 20:00 and 6:00 can be calculated, as shown in Tab.3-Tab.5. The irradiance of the primary mirror on the infrared detector in one day can be calculated by calculating the irradiance of the three points of time mentioned above, as shown in Fig.7. Use this method, the background irradiance of other optical-mechanical elements can be calculated. The total background irradiance of the instrument can be obtained by adding the background irradiance of all optical-mechanical elements on the infrared detector, as shown in Fig.8. According to the irradiance of the point target on channel 3 detector in section 2, which is 1.887×10−3 W/m2, the SCR curve formed by instrument background radiation of channel 3 is calculated, as shown in Fig.9.

      Table 1.  Radiation calculation of primary mirror at different time in infrared band

      TimeTemperature/°CEmissivityM/(W/m2)
      12:00−1.30.040.049 0
      20:00−5.60.040.040 5
      6:002.30.040.057 2

      Table 2.  Energy transfer factor Bij of nodes of primary mirror to detector

      Cond_idNode 1Node 2Area*e*BijBijBji
      1COMET.9284COMET.67902.728 026e-050.018 5320.012 89
      2COMET.9284COMET.67915.640 484e-070.003 455 9210.002 322 2
      3COMET.9284COMET.73081.156 538e-060.013 559 960.001 076
      4COMET.9284COMET.73091.401 487e-060.007 363 9610.003 127 9
      ……
      N-3COMET.9310COMET.73096.087 829E-085.18E-050.001 187
      N-2COMET.9310COMET.73101.103 783E-060.009 397 840.004 149
      N-1COMET.9310COMET.73116.306 837E-070.005 369 770.001 819
      NCOMET.9310COMET.73122.469 422E-060.002 102 5180.007 124

      Table 3.  Radiation calculation of nodes of primary mirror on detector at 12:00

      Cond_idNode 1Node 2BijEj/(W/m2)E/(W/m2)
      1COMET.9284COMET.67790.005 853 20.000 7263.712E-04
      2COMET.9284COMET.67800.000 345 592 10.000 135
      3COMET.9284COMET.72700.003 559 960.000 532
      4COMET.9284COMET.72710.001 363 9610.000 289
      ……
      N-3COMET.9310COMET.72715.18E-052.03E-06
      N-2COMET.9310COMET.72720.009 397 840.000 368
      N-1COMET.9310COMET.72730.000 536 9770.000 21
      NCOMET.9310COMET.72740.002 102 5188.24E-05

      Table 4.  Radiation calculation of nodes of primary mirror on detector at 20:00

      Cond_idNode 1Node 2BijEj/(W/m2)E/(W/m2)
      1COMET.9284COMET.67790.005 853 20.000 5951.820E-04
      2COMET.9284COMET.67800.000 345 592 10.000 111
      3COMET.9284COMET.72700.003 559 960.000 435
      4COMET.9284COMET.72710.001 363 9610.000 236
      ……
      N-3COMET.9310COMET.72715.18E-051.66E-06
      N-2COMET.9310COMET.72720.009 397 840.000 302
      N-1COMET.9310COMET.72730.000 536 9770.000 172
      NCOMET.9310COMET.72740.002 102 5186.75E-05

      Figure 5.  Temperature curve of primary mirror in one day

      Figure 6.  Temperature gradient curve of primary mirror in one day

      Table 5.  Radiation calculation of nodes of primary mirror on detector at 6:00

      Cond_idNode 1Node 2BijEj//(W/m2)E/(W/m2)
      1COMET.9284COMET.67790.005 853 20.000 8545.855E-04
      2COMET.9284COMET.67800.000 345 592 10.000 159
      3COMET.9284COMET.72700.003 559 960.000 625
      4COMET.9284COMET.72710.001 363 9610.000 339
      ……
      N-3COMET.9310COMET.72715.18E-052.39E-06
      N-2COMET.9310COMET.72720.009 397 840.000 433
      N-1COMET.72730.000 536 9770.000 248
      NCOMET.9310COMET.72740.002 102 5189.69E-05

      Figure 7.  Radiation curve of instrument background radiation produced by primary mirror on detector in one day

      Figure 8.  Total instrument background radiation curve generated by optical-mechanical components on detector in one day

      Figure 9.  SCR curve formed by instrument background radiation in one day

    • When the satellite attitude is stable to the ground, the SCR varies with the temperature field of the instrument as shown in Fig.10 according to the image data of channel 3 collected by the camera on April 12, 2017 (corresponding to the solar altitude angle of + 8.8 °). The variation range of SCR varies from 0.059 to 0.096 with time as the x-axis.

      Figure 10.  Variation curve of SCR at band 3 in one day

      The instrument background radiation of Channel 3 is simulated by TRSLIM and TSLAM respectively. Table 6 shows the simulated SCR and the on-orbit SCR. 3:00, 8:00, 12:00, 15:00, 18:00 and 22:00 are selected for comparison, and the TRSLIM, the TSLAM and on-orbit test are compared. The error of TRSLIM and on-orbit measurement is less than 17%, and that of TSLAM and on-orbit measurement method is 114%. It can be seen that the simulation accuracy of TRSLIM is much higher than that of TSLAM.

      In TSLAM, the temperature field is several characteristic temperature conditions provided by the thermal control personnel, and then these temperature fields are manually assigned to the optical and mechanical components in TracePro for background radiation simulation. It is impossible to calculate the changing temperature field in real time, and the simulation time takes one or two weeks. If TRSLIM is used, the temperature field of each optical machine element is the real-time temperature field on-orbit simulated directly by the thermal radiation stray light integrated software. Combined with the radiation factor, the instrument background radiation can be calculated directly. The simulation time is only half a day, which greatly improves the efficiency.

      It can be seen that TRSLIM has two advantages over TSLAM. First, the real-time calculation of the SCR can be obtained by the temperature field changing with time. The background radiation analysis and thermal analysis is integrated into one process. Secondly, the non-uniform temperature field in orbit can be simulated. The temperature field is closer to the actual situation in orbit, and the simulation accuracy is improved greatly.

      Table 6.  Comparison of SCR from on-orbit measured, TRSLIM and TSLAM

      TimeOn-orbit test valuesTRSLIMTSLAM
      Simulation valueErrors from on-orbit test valuesSimulation valueErrors from on-orbit test values
      3:000.0590.061.69%0.11289.83%
      8:000.0640.07415.63%0.13103.13%
      12:000.0780.09116.67%0.14484.62%
      15:000.0890.10315.73%0.189112.36%
      18:000.0950.10914.74%0.204114.74%
      22:000.0890.09810.11%0.179101.12%
    • In this paper, the TRSLIM is used to analyze the instrument background radiation of on-orbit infrared camera. The approximate temperature and uniform temperature field in the TSLAM are replaced by the real temperature field which tends to be on-orbit and has a temperature gradient. Combined with the radiation transfer factor, the background radiation of the instrument and the SCR of the camera on the detector is calculated. Comparing the SCR calculated by TRSLIM and TSLAM, the error of TRSLIM is less than 17%. Compared with the TSLAM, the TRSLIM is real-time and intuitive, the calculation results are closer to the actual situation in orbit, and has higher simulation efficiency and accuracy.

WeChat 关注分享

返回顶部

目录

    /

    返回文章
    返回