-
该节首先介绍美国防太空架构基本情况,再详述天基跟踪层星座组成及主要设计参数。
-
国防太空架构在航天发展局2019年7月1日发布的信息征询书中首先被提出,旨在构建一种灵活、弹性、敏捷的系统架构,以应对高超声速武器等新型威胁目标,进一步保持美在太空领域的优势。根据信息征询书,国防太空架构包括7个功能层[1],如图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 accuracy Obit velocity accuracy 10 m (3σ) 0.02 m/s (3σ) TT & C frequency Data rate Bit error rate Ka (30-31 GHz) 1 Mbps 1e-7 Data downlink frequency Data rate Ka (20.2-21.2 GHz) At least 500 Mbps (250 GB/d) V-band Crosslink frequency Transmission rate 40-75 GHz 1 Gb/s 第0阶段8颗卫星将部署在两个轨道面内,每个轨道面4颗,为确保两颗卫星形成立体跟踪能力,星间最优相位约为30°,且最大应不超过60°。0批次8颗跟踪层卫星在轨主要任务:一是验证新的低轨高灵敏度红外相机对高超声速武器的探测能力;二是验证星星、星地链路传输能力;三是验证利用传输层卫星,通过战术数据链传输目标跟踪信息的能力;四是验证不同供应商提供卫星之间的互操作性;五是验证星上处理、星间通信以及提升星间通信带宽算法的能力;六是开发全球目标跟踪的作战概念;七是评估未来全功能系统对传感器数量、类型、星地数据处理以及通信等的需求。
-
基于前述跟踪层卫星主要参数,重点对其星座覆盖性能、相机探测能力、目标跟踪能力等建模,研究分析其星座核心能力。
-
跟踪层系统对地覆盖性能与星座轨道高度、星间相位、卫星数量、相机瞬时探测视场、总探测视场等相关[3],如图2所示。
假设跟踪层卫星相机为矩形视场,其视场张角称为瞬时探测视场角,记为
$\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所示。
根据目标视轴指向的定义:
$$ {{\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颗卫星,但考虑低轨卫星轨道运动特定,实际所需卫星数量更多。跟踪层卫星设计采用总探测视场覆盖方式,通过整星机动实现对地可见范围内任意位置目标的跟踪。由于
$\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 parameters Satellite number Coverage rate for at least1 valid access Coverage rate for at least2 valid access 7/5/2/180 35 99.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所示。
根据跟踪层卫星载荷在1.5 microflick[2]背景探测能力要求,并结合1000 km轨道高度、450 μrad相机角分辨率等参数推断,跟踪层卫星载荷探测灵敏度优于700 W/sr;但是受背景杂波影响,单帧图像信杂比较低,跟踪层卫星必需通过高帧频探测降低信息处理算法对信杂比的需求或研究低信杂比条件下目标检测与跟踪算法。
-
在表2跟踪层卫星星座参数实现对地表高度20~100 km高度空间双重覆盖条件下,结合公式对跟踪层卫星星座全球目标跟踪性能进行了分析,如图6所示。
从图6仿真结果可以推断,跟踪层卫星完成组网后,在全球范围内在对高超声速武器最优跟踪精度将达到200 m左右,较现有高轨天基红外系统(SBIRS)有了极大的改善。
Performance modeling and simulation of space based tracking layer constellation
-
摘要: 随着高超声速武器技术的快速发展和实战化部署,美国有针对性地提出了建设国防太空架构跟踪层计划的设想。首先,介绍了跟踪层卫星主要设计性能参数;其次,建立了跟踪层星座覆盖性能分析模型、探测能力分析模型和跟踪性能分析模型;最后,通过仿真计算对跟踪层系统性能进行了分析,并推断了系统全星座最小卫星数量、探测灵敏度和全球目标最优跟踪精度等核心能力。分析结果对相关系统信息处理算法研究、载荷与星座方案设计具有重要参考价值。Abstract: With the rapid development of technology and practical deployment of hypersonic weapon, the United States has carried out a new plan of building Tracking Layer of National Defense Space Architecture (NDSA). Firstly, the main design parameters of tracking layer were introduced; Secondly, the analysis model of tracking layer constellation coverage performance, detection capability and tracking performance were established; Finally, the performance of tracking layer was analyzed through simulation and calculation, and the core capabilities such as the minimum satellites number of the whole constellation, detection sensitivity and optimal global target tracking accuracy were deduced. The analysis results have important reference value for the research of data processing algorithm and the design of similar payload and constellation.
-
表 1 跟踪层卫星其他主要参数
Table 1. Other parameters of tracking layer satellite
Obit position accuracy Obit velocity accuracy 10 m (3σ) 0.02 m/s (3σ) TT & C frequency Data rate Bit error rate Ka (30-31 GHz) 1 Mbps 1e-7 Data downlink frequency Data rate Ka (20.2-21.2 GHz) At least 500 Mbps (250 GB/d) V-band Crosslink frequency Transmission rate 40-75 GHz 1 Gb/s 表 2 覆盖性能与星座参数关系表
Table 2. Relationship between coverage performance and constellation parameters
Constellation parameters Satellite number Coverage rate for at least1 valid access Coverage rate for at least2 valid access 7/5/2/180 35 99.3% 52.48% 14/7/2/180 98 100% 99.46% -
[1] Hu Yini, Zhong Jiangshan, Wei Chenxi, et al. Analysis of US next generation space architecture [J]. Spacecraft Engineering, 2021, 30(2): 108-117. (in Chinese) doi: 10.3969/j.issn.1673-8748.2021.02.015 [2] SDA. Technical requirements of tracking layer to WFOV[EB/OL]. [2020-05-11]. http://beta.sam.gov/. [3] Wang Bo, Xu Dan, Chen Yanjun, et al. Analysis of the coverage of LEO constellation based infrared staring sensors [J]. Journal of Hunan University (Nature Sciences), 2009, 36(10): 68-74. (in Chinese) [4] Xue Yonghong, An Wei, Zhang Tao, et al. Real-time line of sight calibration algorithm for high earth orbit infrared scanning sensor [J]. Acta Optica Sinica, 2013, 33(6): 0628003. (in Chinese) doi: 10.3788/AOS201333.0628003 [5] Liu Tao, Yang Degui, Jiang Weidong, et al. IR imaging modeling and simulation of ballistic target group outside atmosphere [J]. Infrared and Laser Engineering, 2011, 40(4): 589-594. (in Chinese) doi: 10.3969/j.issn.1007-2276.2011.04.003 [6] Zhang Tao, Xue Yonghong, Ding Yonghe, et al. Analysis of imaging characteristics of background in elliptical orbit infrared surveillance platform [J]. Infrared and Laser Engineering, 2016, 45(3): 0326001. (in Chinese) doi: 10.3788/irla201645.0326001 [7] 孙建奇, 方小平. 红外物理[M]. 西安: 西安电子科技大学出版社, 2004. Sun Jianqi, Fang Xiaoping. Infrared Physics[M]. Xi'an: Xidian University Press, 2004. (in Chinese) [8] 盛卫东. 天基光学传感器网络目标跟踪技术研究[D]. 长沙: 国防科学技术大学, 2011. Sheng Weidong. Research on target tracking technologies for space-based optical sensor network[D]. Changsha: National University of Defense Technology, 2011. (in Chinese) [9] Yang Xiao, Niu Qinglin, He Zhihong, et al. Analysis of infrared radiation characteristics and detectability of HTV-2-like hypersonic gliding aircrafts [J]. Acta Optica Sinica, 2017, 37(12): 1204001. (in Chinese) doi: 10.3788/AOS201737.1204001 [10] Wang Yingying, He Ping, Meng Changliang, et al. Detection distance computation of the space borne infrared detector for hypersonic flight vehicles [J]. Infrared and Laser Engineering, 2019, 48(7): 0704003. (in Chinese) doi: 10.3788/IRLA201948.0704003 [11] Yu Xiaoya, Liu Lituo, Li Rui, et al. Measurements of absolute radiative emissions for supersonic reentry [J]. Chinese Journal of Optics, 2020, 13(1): 87-94. (in Chinese) doi: 10.3788/co.20201301.0087