留言板

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

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

天基跟踪层星座性能建模与仿真

薛永宏 樊士伟 乔凯 张磊

薛永宏, 樊士伟, 乔凯, 张磊. 天基跟踪层星座性能建模与仿真[J]. 红外与激光工程, 2021, 50(11): 20210067. doi: 10.3788/IRLA20210067
引用本文: 薛永宏, 樊士伟, 乔凯, 张磊. 天基跟踪层星座性能建模与仿真[J]. 红外与激光工程, 2021, 50(11): 20210067. doi: 10.3788/IRLA20210067
Xue Yonghong, Fan Shiwei, Qiao Kai, Zhang Lei. Performance modeling and simulation of space based tracking layer constellation[J]. Infrared and Laser Engineering, 2021, 50(11): 20210067. doi: 10.3788/IRLA20210067
Citation: Xue Yonghong, Fan Shiwei, Qiao Kai, Zhang Lei. Performance modeling and simulation of space based tracking layer constellation[J]. Infrared and Laser Engineering, 2021, 50(11): 20210067. doi: 10.3788/IRLA20210067

天基跟踪层星座性能建模与仿真

doi: 10.3788/IRLA20210067
详细信息
    作者简介:

    薛永宏,男,助理研究员,博士,主要研究方向为天基信息获取与处理

  • 中图分类号: TJ86; V47

Performance modeling and simulation of space based tracking layer constellation

图(6) / 表(2)
计量
  • 文章访问数:  319
  • HTML全文浏览量:  114
  • PDF下载量:  64
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-26
  • 修回日期:  2021-04-28
  • 网络出版日期:  2021-12-03
  • 刊出日期:  2021-11-30

天基跟踪层星座性能建模与仿真

doi: 10.3788/IRLA20210067
    作者简介:

    薛永宏,男,助理研究员,博士,主要研究方向为天基信息获取与处理

  • 中图分类号: TJ86; V47

摘要: 随着高超声速武器技术的快速发展和实战化部署,美国有针对性地提出了建设国防太空架构跟踪层计划的设想。首先,介绍了跟踪层卫星主要设计性能参数;其次,建立了跟踪层星座覆盖性能分析模型、探测能力分析模型和跟踪性能分析模型;最后,通过仿真计算对跟踪层系统性能进行了分析,并推断了系统全星座最小卫星数量、探测灵敏度和全球目标最优跟踪精度等核心能力。分析结果对相关系统信息处理算法研究、载荷与星座方案设计具有重要参考价值。

English Abstract

    • 天基跟踪层是美国防太空架构(National Defense Space Architecture, NDSA)的重要组成部分,也是美下一代持续红外系统(Next Generation Overhead Persistent Infrared System, NG-OPIR)的重要组成部分,主要承担高超声速武器等先进导弹威胁目标(Advanced Missile Threats)探测、预警、跟踪和指示等任务。

      2020年10月5日,美航天发展局(Space Development Agency, SDA)向SpaceX和L3Harris两家公司分别授予价值1.935亿美元和1.49亿美元的合同,用于天基跟踪层0批次8颗卫星研制,其中每家公司负责研制4颗卫星,计划2022年9月起的6个月内完成发射部署。

      文中以美天基跟踪层为研究背景,重点对其对地覆盖性能、探测性能、跟踪性能等进行仿真分析,并推断了跟踪层系统性能。论文主要结构安排如下:第一部分介绍天基跟踪层主要任务及设计参数;第二部分建立覆盖性能、探测性能、跟踪性能分析模型;第三部分对天基跟踪层性能进行仿真分析,并给出相关结论。

    • 该节首先介绍美国防太空架构基本情况,再详述天基跟踪层星座组成及主要设计参数。

    • 国防太空架构在航天发展局2019年7月1日发布的信息征询书中首先被提出,旨在构建一种灵活、弹性、敏捷的系统架构,以应对高超声速武器等新型威胁目标,进一步保持美在太空领域的优势。根据信息征询书,国防太空架构包括7个功能层[1],如图1所示。

      图  1  国防太空架构[1]

      Figure 1.  National Defense Space Architecture[1]

      (1)传输层:为全球范围内的作战人员或平台提供稳定、可靠、低时延的数据和通信连接。

      (2)战斗管理层:为战役行动中时敏杀伤链闭合提供任务分派、指挥控制及数据分发等功能。

      (3)跟踪层:对全球高超声速武器等先进导弹威胁目标进行预警、跟踪、指示等。

      (4)监视层:对高时敏目标发射前征候进行全天时(24×7)、全天候监视,支撑目标指示。

      (5)导航层:在GPS拒止环境下,提供定位、导航、授时功能。

      (6)威慑层:抑制地月空间的敌对行动。

      (7)支撑层:确保地面和发射部分对太空架构的快速响应支持。

      2020年10月,航天发展局公布了NDSA第0批次采购计划,共包括28颗卫星;其中,20颗为传输层卫星,8颗为跟踪层卫星。0批次系统包括六大验证目标:一是验证低延迟数据传输系统;二是演示从外部天基传感器向作战人员传送数据;三是演示有限的C3作战管理功能;四是从传输层网络向作战人员传输综合广播系统数据(Integrated Broadcast System);五是通过网络实时存储、中继和传输16链路数据;六是不依赖于GPS系统的通用时空基准。第1批次计划2024年底完成,卫星数量将增加到100~150颗,由多家不同的单位承担研制任务。

    • 卫星部署在1000 km轨道高度、80°~100°倾角的圆轨道,采用三轴稳定控制方式,平台具备不小于0.3(°)/s机动指向能力,能够确保相机在100 s内指向任意区域并稳定探测;卫星寿命不少于4年。

      每颗卫星装载一台宽视场红外相机,相机装配红外短波谱段(Short Wave Infrared, SWIR)和红外短波透过谱段(SWIR See-to-Ground, SWSTG),相机视场50°×50°,像元角分辨率450 μrad,任意视场地面空间分辨率不大于1.5 km;卫星利用恒星进行光轴校正,光轴校正精度优于0.1个像元;星上具备数据存储和处理能力,能够存储不小于50 min原始探测数据,星上处理算法参数如门限、阈值等可在轨调整。每颗卫星装载4条激光链路,其中两个用于同轨道面卫星建链、两个用于异轨道面卫星建链或星地链路);装载一个V频段星间链路,作为激光链路的备份;星间链路采用与传输层卫星相同的路由协议以及基于IP的安全策略;星上时统精度优于100 ns。跟踪层卫星其他主要参数见表1[2]

      表 1  跟踪层卫星其他主要参数

      Table 1.  Other parameters of tracking layer satellite

      Obit position accuracyObit velocity accuracy
      10 m (3σ)0.02 m/s (3σ)
      TT & C frequencyData rateBit error rate
      Ka (30-31 GHz)1 Mbps1e-7
      Data downlink frequencyData rate
      Ka (20.2-21.2 GHz)At least 500 Mbps (250 GB/d)
      V-band Crosslink frequencyTransmission rate
      40-75 GHz1 Gb/s

      第0阶段8颗卫星将部署在两个轨道面内,每个轨道面4颗,为确保两颗卫星形成立体跟踪能力,星间最优相位约为30°,且最大应不超过60°。0批次8颗跟踪层卫星在轨主要任务:一是验证新的低轨高灵敏度红外相机对高超声速武器的探测能力;二是验证星星、星地链路传输能力;三是验证利用传输层卫星,通过战术数据链传输目标跟踪信息的能力;四是验证不同供应商提供卫星之间的互操作性;五是验证星上处理、星间通信以及提升星间通信带宽算法的能力;六是开发全球目标跟踪的作战概念;七是评估未来全功能系统对传感器数量、类型、星地数据处理以及通信等的需求。

    • 基于前述跟踪层卫星主要参数,重点对其星座覆盖性能、相机探测能力、目标跟踪能力等建模,研究分析其星座核心能力。

    • 跟踪层系统对地覆盖性能与星座轨道高度、星间相位、卫星数量、相机瞬时探测视场、总探测视场等相关[3],如图2所示。

      图  2  卫星姿态与传感器指向关系示意图

      Figure 2.  Sketch map of satellite attitude and sensor point

      假设跟踪层卫星相机为矩形视场,其视场张角称为瞬时探测视场角,记为$\alpha $;通过卫星姿态机动或转台转动可扩展相机探测范围,相机对地可见探测范围张角称为总探测视场角,记为$\;\beta $;另假设卫星轨道高度为$H$,地球半径为$R$,则可得对地覆盖率与卫星数量间的关系:

      $$ {N_{S}=}\left\{ {\begin{array}{*{20}{l}} {\dfrac{{4\pi {R^2}}}{{{{\left( {\alpha H} \right)}^2}}}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\rm Instantaneous\;field\;of\;view\;(IFOV)}}\\ {\dfrac{{4\pi {R^2}}}{{\pi {{\left( {\beta H/2} \right)}^2}}}\;\;\;\;\;\;\;\;\;{\rm Total\;field\;of\;view\;(TFOV)}} \end{array}} \right. $$ (1)
    • 跟踪层卫星探测性能分析需首先建立红外相机成像模型,具体包括几何投影模型[4-5]和信号传递模型[6]两部分。其中,几何投影模型主要描述目标与背景在像平面的成像位置关系,论文不再赘述;信号传递模型主要描述像平面上目标与背景能量幅度关系,直接与系统探测性能相关。

      根据参考文献[6-7],相机获取目标与背景能量分别可表示为:

      $$ \left\{ \begin{array}{l} Q_{{\rm{foc}}}^{\left( {\rm{B}} \right)} = E_{{\rm{foc}}}^{\left( {\rm{B}} \right)} \cdot {t_{d}} \cdot d_{{\rm{foc}}}^2 \\ Q_{\rm foc}^{\left( {\rm T} \right)} = E_{\rm foc}^{\left( {\rm T} \right)} \cdot {t_{d}} \cdot \displaystyle\int_{{y_c} - \frac{{{d_{\rm foc}}}}{2}}^{{y_c} + \frac{{{d_{\rm foc}}}}{2}} {\displaystyle\int_{{x_c} - \frac{{ - {d_{\rm foc}}}}{2}}^{{x_c} + \frac{{{d_{\rm foc}}}}{2}} {{p_{({x_0},{y_0})}}(x,y){\rm d}x{\rm d}y} } \\ \end{array} \right. $$ (2)

      式中:${t_{d}}$为相机积分时间;${d_{{\rm{foc}}}}$为单个像元尺寸;$({x_c},{y_c})$为目标投影在像平面的位置;${p_{({x_0},{y_0})}}(x,y)$为相机点扩散函数,反映相机能量集中度,记为$\;\rho $$E_{\rm foc}^{\left( {\rm T} \right)}$$E_{{\rm{foc}}}^{\left( {\rm{B}} \right)}$分别为目标与背景在像平面的辐射照度。

      $$ \left\{ \begin{array}{l} E_{{\rm{foc}}}^{\left( {\rm{B}} \right)} = \displaystyle\int {{L_\lambda }} {\rm d}\lambda \cdot \dfrac{{f_{{\rm{IFOV}}}^2 \cdot {\tau _0} \cdot \pi \cdot {D^2} \cdot \cos \left( {{\theta _{\rm{0}}}} \right)}}{{4d_{{\rm{foc}}}^2}} \\ E_{\rm foc}^{\left( {\rm T} \right)} = \dfrac{{\displaystyle\int {\tau _a^\lambda \cdot {I_\lambda }{\rm d}\lambda } }}{{{l^2}}} \cdot \dfrac{{{\tau _0} \cdot \pi \cdot {D^2} \cdot \cos \left( {{\theta _1}} \right)}}{{4d_{{\rm{foc}}}^2}} \\ \end{array} \right. $$ (3)

      式中:${L_\lambda }$为大气背景光谱辐射亮度;${I_\lambda }$为目标光谱辐射强度;$l$为探测距离;$\tau _a^\lambda $为大气透过率;${f_{{\rm{IFOV}}}}$为相机单个像元的角分辨率;${\theta _{\rm{0}}}$${\theta _1}$分别为背景与目标入射光线与相机口面法线的夹角;${\tau _0}$为相机成像效率;$D$为相机光学口径尺寸。

      采用目标与背景杂波辐射对比度[6],即信杂比表述相机探测性能,利用公式(2)、(3)可得:

      $$ CR = \frac{{Q_{\rm foc}^{\left( {\rm T} \right)}}}{{Q_{{\rm{foc}}}^{\left( {\rm{B}} \right)}}} = \frac{{\rho \cdot \displaystyle\int {\tau _a^\lambda \cdot {I_\lambda }{\rm d}\lambda } }}{{{l^2} \cdot f_{{\rm{IFOV}}}^2 \cdot \displaystyle\int {{L_\lambda }} {\rm d}\lambda }} $$ (4)

      由此可知,探测性能主要与相机探测谱段、像元角分辨率、探测距离、目标与背景光谱辐射强度/亮度、大气透过率等因素相关;其中,目标与背景光谱辐射强度/亮度、大气透过率为客观因素,不会随着相机或人为因素改变而改变。

    • 跟踪层卫星对高超声速目标跟踪属于仅测角跟踪系统,必须利用两颗或者多颗卫星的观测数据实现对目标的三维定位跟踪[8],如图3所示。

      图  3  光学交叉定位原理图

      Figure 3.  Schematic diagram of optical cross location

      根据目标视轴指向的定义:

      $$ {{\boldsymbol{v}}_T} = \frac{{{\boldsymbol{r}} - {{\boldsymbol{r}}_S}}}{{\left\| {{\boldsymbol{r}} - {{\boldsymbol{r}}_S}} \right\|}} $$ (5)

      式中:${\boldsymbol{r}} = {\left( {x,y,{\textit{z}}} \right)^{\rm T}}$${{\boldsymbol{r}}_S} = {\left( {{x_S},{y_S},{{\textit{z}}_S}} \right)^{\rm T}}$分别表示目标和卫星的位置矢量。对上式求一阶微分可得:

      $$ {\rm d}{{\boldsymbol{v}}_T} = {\boldsymbol{H}}\left( {{\rm d}{\boldsymbol{r}} - {\rm d}{{\boldsymbol{r}}_S}} \right) $$ (6)

      进一步改写可得到:

      $$ {\boldsymbol{H}}{\rm d}{\boldsymbol{r}} = {\rm d}{{\boldsymbol{v}}_T} + {\boldsymbol{H}}{\rm d}{{\boldsymbol{r}}_S} $$ (7)

      可知,目标跟踪定位性能主要与卫星位置误差和目标视轴指向误差相关。式中${\boldsymbol{H}}$反映了目标视轴指向误差对目标位置误差、卫星位置误差的灵敏程度,可表示为:

      $$\begin{split}& {\boldsymbol{H}} =\\& \frac{{\left( {\begin{array}{*{20}{c}} {{{\left( {y - {y_S}} \right)}^2} + {{\left( {{\textit{z}} - {{\textit{z}}_S}} \right)}^2}}&{ - \left( {x - {x_S}} \right)\left( {y - {y_S}} \right)}&{ - \left( {x - {x_S}} \right)\left( {{\textit{z}} - {{\textit{z}}_S}} \right)} \\ { - \left( {x - {x_S}} \right)\left( {y - {y_S}} \right)}&{{{\left( {x - {x_S}} \right)}^2} + {{\left( {{\textit{z}} - {{\textit{z}}_S}} \right)}^2}}&{ - \left( {y - {y_S}} \right)\left( {{\textit{z}} - {{\textit{z}}_S}} \right)} \\ { - \left( {x - {x_S}} \right)\left( {{\textit{z}} - {{\textit{z}}_S}} \right)}&{ - \left( {y - {y_S}} \right)\left( {{\textit{z}} - {{\textit{z}}_S}} \right)}&{{{\left( {x - {x_S}} \right)}^2} + {{\left( {y - {y_S}} \right)}^2}} \end{array}} \right)}}{{{{\left\| {{\boldsymbol{r}} - {{\boldsymbol{r}}_{\rm{S}}}} \right\|}^3}}} \end{split}$$ (8)

      由于$\left| {\boldsymbol{H}} \right| \equiv 0$${\boldsymbol{H}}$不可逆;因此单星观测条件下无法求解目标跟踪性能,当且仅当卫星数量${N_S} \geqslant 2$时,才能对目标跟踪性能进行分析。进一步构建目标跟踪性能分析模型为:

      $$ \left( {\begin{array}{*{20}{c}} {{{\boldsymbol{H}}_1}} \\ \vdots \\ {{{\boldsymbol{H}}_{{N_S}}}} \end{array}} \right){\rm d}{\boldsymbol{r}} = \left( {\begin{array}{*{20}{c}} {{\rm d}{{\boldsymbol{v}}_{T,1}}} \\ \vdots \\ {{\rm d}{{\boldsymbol{v}}_{T,{N_S}}}} \end{array}} \right) + {{{\rm{\varLambda}} }}\left( {\begin{array}{*{20}{c}} {{\rm d}{{\boldsymbol{r}}_{S,1}}} \\ \vdots \\ {{\rm d}{{\boldsymbol{r}}_{S,{N_S}}}} \end{array}} \right) $$ (9)

      其中,${\boldsymbol{\varLambda }} = diag\left( {{{\boldsymbol{H}}_1}, \cdots ,{{\boldsymbol{H}}_{{N_S}}}} \right)$;进一步求解可得:

      $$ {\rm d}{\boldsymbol{r}} = {\left( {{{\boldsymbol{\Sigma }}_H}} \right)^{ - 1}}{{\boldsymbol{C}}_H}\left( {\begin{array}{*{20}{c}} {{\rm d}{{\boldsymbol{v}}_{T,1}}} \\ \vdots \\ {{\rm d}{{\boldsymbol{v}}_{T,{N_S}}}} \end{array}} \right) + {\left( {{{\boldsymbol{\Sigma }}_H}} \right)^{ - 1}}{{\boldsymbol{B}}_H}\left( {\begin{array}{*{20}{c}} {{\rm d}{{\boldsymbol{r}}_{S,1}}} \\ \vdots \\ {{\rm d}{{\boldsymbol{r}}_{S,{N_S}}}} \end{array}} \right) $$ (10)

      其中,

      $$ \begin{array}{l} {{\boldsymbol{\Sigma }}_H} = \displaystyle\sum\limits_{i = 1}^{{N_S}} {{\boldsymbol{H}}_i^{\rm T}{{\boldsymbol{H}}_i}} \\ {{\boldsymbol{C}}_H} = \left( {{\boldsymbol{H}}_1^{\rm T}, \cdots ,{\boldsymbol{H}}_{{N_S}}^{\rm T}} \right) \\ {{\boldsymbol{B}}_H} = \left( {{\boldsymbol{H}}_1^{\rm T}{{\boldsymbol{H}}_1}, \cdots ,{\boldsymbol{H}}_{{N_S}}^{\rm T}{{\boldsymbol{H}}_{{N_S}}}} \right) \\ \end{array} $$ (11)
    • 根据跟踪层卫星设计参数,对系统覆盖性能、探测能力及跟踪性能进行仿真计算。

    • 根据跟踪层系统参数计算可知,利用瞬时探测视场视线对地100%一重覆盖,至少需671颗卫星,如图4所示;利用卫星总探测视场($\;\beta {\rm{ = }}120°$)实现对地100%一重覆盖理论上仅需15颗卫星,但考虑低轨卫星轨道运动特定,实际所需卫星数量更多。

      图  4  瞬时视场与一重覆盖卫星数量关系

      Figure 4.  Relationship between FOV and satellite number under one coverage

      跟踪层卫星设计采用总探测视场覆盖方式,通过整星机动实现对地可见范围内任意位置目标的跟踪。由于$\alpha {\rm{ = }}50° \ll \beta {\rm{ = }}120°$,故而可推断跟踪层卫星不具备目标搜索发现能力,必需在其他系统(如高轨SBIRS卫星)引导信息支持下,实现对高超声速武器的连续跟踪探测。

      基于覆盖性能计算软件工具,通过构建场景进一步分析跟踪层星座参数、卫星数量与覆盖性能关系见表2。为实现对高超声速武器的立体跟踪定位,跟踪层卫星必需实现对高超声速武器飞行空间,即距地表20~100 km高度[9]实现双重覆盖,由此推断,跟踪层卫星总数量需不少于98颗,考虑一定的功能备份,卫星总数量不应少于100颗。

      表 2  覆盖性能与星座参数关系表

      Table 2.  Relationship between coverage performance and constellation parameters

      Constellation parametersSatellite numberCoverage rate for at least1 valid accessCoverage rate for at least2 valid access
      7/5/2/1803599.3%52.48%
      14/7/2/180 98 100% 99.46%
    • 根据参考文献[9-11],高超声速武器辐射能量主要集中在2~5 μm红外短波谱段,与跟踪层卫星谱段设计相匹配。其中,高超声速武器在红外短波谱段的光谱辐射强度在10000 W·sr−1·µm−1以上[9-10],背景在红外短波谱段的光谱辐射亮度约为0.01W·sr−1·m−2·µm−1[9-10]。假设大气综合吸收率为0.5,则结合式可得到探测距离、相机角分辨率与目标背景辐射对比度CR间的关系,如图5所示。

      图  5  探测距离与角分辨率与辐射对比度关系图

      Figure 5.  Relationship between range, IFOV and CR

      根据跟踪层卫星载荷在1.5 microflick[2]背景探测能力要求,并结合1000 km轨道高度、450 μrad相机角分辨率等参数推断,跟踪层卫星载荷探测灵敏度优于700 W/sr;但是受背景杂波影响,单帧图像信杂比较低,跟踪层卫星必需通过高帧频探测降低信息处理算法对信杂比的需求或研究低信杂比条件下目标检测与跟踪算法。

    • 表2跟踪层卫星星座参数实现对地表高度20~100 km高度空间双重覆盖条件下,结合公式对跟踪层卫星星座全球目标跟踪性能进行了分析,如图6所示。

      图  6  跟踪层星座全球GDOP图

      Figure 6.  Global GDOP map of tracking layer

      图6仿真结果可以推断,跟踪层卫星完成组网后,在全球范围内在对高超声速武器最优跟踪精度将达到200 m左右,较现有高轨天基红外系统(SBIRS)有了极大的改善。

    • 跟踪层星座设计与目标、背景特性等密切相关,通过对跟踪层卫星覆盖性能、探测性能、跟踪性能进行建模与仿真,推断了跟踪层星座卫星数量、探测能力边界、最优跟踪性能等核心指标参数,主要启示有:一是在目标探测跟踪方面,从载荷探测灵敏度、跟踪性能看,跟踪层卫星对高超声速武器具有较好的探测效果,但受制于卫星姿态机动跟踪模式,星下点附近对高超目标跟踪能力较弱;二是在目标识别方面,跟踪层卫星主要采用红外短波谱段,仅依赖红外辐射特性难以有效进行目标分类判别,但综合利用目标飞行高度、速度等特征信息,系统具备高超声速武器分类判别能力;三是在信息处理算法方面,跟踪层卫星探测图像中目标背景辐射对比度较低,需研究低信噪比信息处理算法确保对目标的可靠检测处理。

参考文献 (11)

目录

    /

    返回文章
    返回