1. 引言
在航天发射任务中,要求测量设备尤其是新研或改造后在正式投入任务使用前和使用一段时间后,必须进行严格的精度鉴定。现行的设备精度鉴定通常采用传统的飞机校飞模式,由校飞组织单位牵头协调信标合作单位、待校飞单位、数据处理单位等,协调飞机架次,设计飞行航线,收集测量数据和飞机航迹数据,处理测量数据和航迹数据最终给出待鉴定设备的精度鉴定结果。其中飞机架次还受环境气候等因素制约,严重影响设备的精度鉴定效率。2019年,我国天平一号卫星系统服役并投入使用,天平一号卫星系统在空间建立了目标散射特性基准和目标运动状态基准系统[1],为空间目标探测设备的鉴定与标定提供服务,为空间目标探测技术的研究与实验提供支持。
但从目前的使用情况来看,天平一号卫星系统的应用还不够充分,多数设备单纯用于目标跟踪训练,基于其轨道数据的精度鉴定应用较少,未能有效挖掘大量的设备测量数据和卫星轨道数据的应用效能;并且缺乏应用层面统一的跟踪精度鉴定系统,也未直接有效地提供测量数据服务与计算处理服务等功能,导致未能充分发挥该卫星系统的使用效益。因此亟需建设一个模式统一,功能集成度高,响应速度快,使用过程简单方便的测量设备精度鉴定系统,以解决传统精度鉴定工作中存在的问题,为所有测量设备用户提供统一的、快速的、简便的精度鉴定服务。
2. 功能架构与流程设计
2.1. 功能架构设计
本文设计并实现航天测量设备的快速、简便、统一的精度鉴定系统,以网络化服务模式收集各设备的跟踪测量数据,运用网络爬虫技术[2]自动获取天平一号的轨道数据,进行数据处理、分析、存储及显示,以线上服务模式为各设备用户提供任务规划、精度鉴定、精度优化、精度监测等服务,大幅度提升设备精度鉴定的效率,同时也为用户提供长效的设备精度变化监测,使其更好地掌握设备的技术性能指标。各种服务的主要功能包括:
(1) 标校星跟踪任务规划
获取标校星预报轨道数据,根据各测量站的地理分布,自动化分析计算各设备的可视情况,同时给出时间规划和空间规划,以完成针对单站跟踪或多站联合跟踪开展跟踪任务规划、引导弹道生成等。
(2) 跟踪精度鉴定
基于各测量设备跟踪标校星的测量数据与标校星精密轨道数据的设备跟踪精度鉴定方法,完成各测量设备的距离精度、方位精度、俯仰精度以及速度精度的鉴定,为用户提供以网络化服务模式的统一的精度鉴定结果。
(3) 测量数据精度优化
通过对测量数据与精密轨道数据的综合分析,研究设备标定数据质量评估方法、测量设备测量数据的误差修正算法,对设备进行测量数据误差修正分析,给出距离零值、角度零值、重力下垂、大盘不水平、机电偏差、方位俯仰不正交度等参数的估算值,以为用户提供数据精度优化的建议。
(4) 设备跟踪精度监测
通过对历次设备的精度鉴定结果数据存储至数据库,积累成精度变化监测的源数据,完成对设备精度的监测。主要包含距离、方位角、俯仰角及速度四个测元的总误差、随机误差和系统误差。
系统的总体功能架构如图1所示。
Figure 1. Overall functional architecture
图1. 总体功能架构图
2.2. 流程设计
2.2.1. 任务规划
任务规划是自动获取天平一号卫星预报轨道数据,完成针对单站跟踪或多站联合跟踪开展鉴定任务的时间和空间规划、引导弹道生成等。主要通过网络爬虫技术定时下载预报轨道数据,根据各用户注册时提供的站址坐标进行坐标转换,再进行可视性判断,默认为俯仰2˚的判断标准,之后给出所有用户近三日内的所有可跟踪任务的时间规划视图,并提供各可跟踪时段的空间(方位、俯仰角度)覆盖范围,便于用户根据覆盖性需求开展空间规划。同时生成对应规划时段的引导弹道文件,以供用户下载使用。任务规划流程如图2所示。
Figure 2. Task planning process
图2. 任务规划流程
2.2.2. 精度鉴定
精度鉴定[3] [4]以运行于空间近地轨道的天平一号标校卫星为基准目标。其基本原理为:被鉴定测量设备跟踪测量卫星目标,获取测量数据;同时,获取该卫星对应于测量弧段的精密轨道数据,将测量数据与卫星精确轨道数据进行对比分析,将统计数据与被鉴定设备战技术指标进行比对,从而对该测量设备精度作出评估。设备测量精度反映测量结果与真值的一致程度,一般由系统差和随机差综合表示。
由于受测量误差、系统动态滞后、大气折射[5]等因素影响,测量数据在做精度鉴定前,用户需要先对数据进行动态滞后修正、系统误差修正及电波折射修正等。同时,在与轨道数据进行比对前,还要做数据对齐处理,包括数据合法性校验、平滑滤波、秒时间处理、丢帧处理等,以解决测量数据与卫星轨道数据采样时刻不一致、时间格式错误、跨零跳点、丢帧等问题。在对数据进行前期对齐处理后,运用卫星定轨数据作为比较标准值,对设备测量数据进行精度分析,最终得到测量设备的跟踪精度。精度鉴定流程如图3所示。
Figure 3. Accuracy appraisal process
图3. 精度鉴定流程
2.2.3. 精度优化
精度优化是通过对测量数据与定轨数据的综合分析,运用误差修正模型,对设备进行测量数据误差参数估算,优化数据精度。在做精度优化分析前,除了不做系统误差修正外,其他数据处理过程与精度鉴定过程相同。之后运用误差估算模型结合定轨数据进行误差参数估算,以脉冲雷达为例,估算的参数包括:方位零值、俯仰零值、距离零值、大盘不水平方位角、大盘不水平度、方位俯仰不正交度、重力下垂、方位俯仰机电偏差。精度优化流程如图4所示。同时,运用设备标定数据质量评估方法,结合精度鉴定功能可对估算的误差参数进行闭环验证,检验其精度优化的质量、效果。
Figure 4. Precision optimization process
图4. 精度优化流程
2.2.4. 精度监测
Figure 5. Precision monitoring process
图5. 精度监测流程
精度监测是通过对历次设备的精度鉴定结果数据存储至数据库,积累成精度变化监测的源数据,完成对设备精度的监测。用户通过选定需要查看的数据范围,主要包括跟踪的圈次、时间等,调取精度鉴定结果数据,包含距离、方位角、俯仰角及速度四个测元的总误差、随机误差和系统误差,将所有数据以原始数据和精度变化监视曲线的方式直观地展示给用户,并提供曲线图片和详细数据的下载服务。精度监测流程如图5所示。
3. 关键算法模型
在数据处理过程中,涉及到了动态滞后修正、系统误差修正及电波折射修正,同时还要做数据对齐处理,包括数据合法性校验、平滑滤波、秒时间处理、丢帧处理等,其中有部分处理算法均为常见且较为通用的算法,不再赘述,主要对精度分析与误差估计算法进行介绍。
3.1. 精度计算算法
精度评定的主要指标有系统误差和随机误差,以及总误差。以方位角A为例,各误差的计算如下:
测量数据与天平一号定轨轨道数据比对作差,按公式(1)得到两测量数据之差:
(1)
式中:
为测量设备跟踪天平一号标校卫星得到的方位角度测量数据;
为天平一号的定轨轨道方位角度数据,作为标准比对值;
为区段序号;
为某区段中的采样点序号。
规定:n为400个点,每个区段为20秒400个点。如果在一段内有效数据小于200时,此段就不做处理。
(1) 总误差
由公式(2)计算第j区段的方位角度测量误差的总误差:
(2)
而m个区段测量的方位角度总误差按公式(3)计算:
(3)
(2) 随机误差
第j个区段的测量数据方位角度随机误差按公式(4)的变量差分法计算求解:
(4)
取差分阶数P = 3,步长t = 4,n = 400个点,则:
(5)
而m个区段的测量数据的方位角度随机误差按公式(10)计算:
(6)
(3) 系统误差
第j个区段的测量数据的方位角度系统误差按公式(7)计算:
(7)
而m个区段的测量数据的方位角度系统误差按公式(8)计算:
(8)
3.2. 误差估算模型
系统误差修正项目包括天线座水平误差、方位轴和俯仰轴垂直误差、光机轴平行误差、动态滞后误差、天线重力变形误差等。根据标校得到的各项误差数据,对设备的外测数据实施方位、俯仰、距离进行系统误差的修正[6] [7],为了方便对误差分量估计[8] [9],对方位、俯仰的系统误差修正模型做一个简化处理:
(9)
(10)
(11)
式中:Am为雷达方位测量数值(或码盘提供数值) (˚);Em为雷达俯仰测量数值(或码盘提供数值) (˚);A0为方位零值(˚);E0为俯仰零值(˚);
为方位轴铅锤度(大盘水平度) (˚);Amax为不水平方位角(˚);
为方位轴与俯仰轴的垂直度(˚);
为天线重力下垂(˚);Rt为信标延迟(μs);R0为距离零值(m),其中,Vaz为方位机电不匹配角,Vel为俯仰机电不匹配角。
由式(9)得到
(12)
由式(10)得到
(13)
令
(14)
(15)
(16)
(17)
(18)
(19)
则有
(20)
同样,由式(11)得到
(21)
令
(22)
(23)
(24)
(25)
(26)
则有
(27)
其中,
,
,
分别表示雷达修正值与原始测量值的距离、方位角、仰角残差。
为了估计误差分量,测量设备跟踪天平一号卫星得到测量数据,并用卫星精轨数据代替系统误差修正后的数据,则上节中的残差表示设备测量数据与卫星精轨数据的残差。
根据最小二乘准则可以估算系统误差分量矩阵:
(28)
(29)
(30)
根据式(15)~(19)及式(23)~(26),可以求得系统误差分量。
需要说明的是,算法中对仰角E0,Vel进行了合并处理,只给出了E0 + Vel的估计值,不会影响系统误差修正的结果。
4. 算法验证与系统实现
4.1. 算法验证
4.1.1. 精度计算算法应用验证
对某段测量数据进行精度分析验证,数据经过了合理性检验、跨零跳点修正、系统误差修正和电波折射修正的处理,再与对应时段的天平一号定轨轨道数据进行精度分析。区段数据点数为400,共计算了14段数据,各区段数据的测元总误差、随机误差及系统误差如表1所示,该测量数据的总的误差结果如表2所示。
Table 1. Results of segmented accuracy analysis
表1. 分段精度分析结果
误差类型 |
总误差 |
随机误差 |
系统误差 |
区段 序号 |
|
|
|
|
|
|
|
|
|
|
|
|
1. |
0.014705 |
0.186243 |
37.07957 |
0.241776 |
0.00274 |
0.004133 |
0.834339 |
0.066143 |
0.014447 |
0.186197 |
37.07018 |
0.232553 |
2. |
0.014458 |
0.138293 |
34.30857 |
0.331057 |
0.002896 |
0.002622 |
0.652093 |
0.063379 |
0.014165 |
0.138269 |
34.30237 |
0.324934 |
3. |
0.01403 |
0.101797 |
36.42415 |
0.48121 |
0.002599 |
0.002226 |
0.949732 |
0.109613 |
0.013787 |
0.101772 |
36.41177 |
0.46856 |
4. |
0.015074 |
0.06723 |
29.20641 |
0.327708 |
0.002017 |
0.002511 |
1.009057 |
0.002794 |
0.014938 |
0.067183 |
29.18897 |
0.327696 |
5. |
0.015722 |
0.053362 |
30.85997 |
0.3769 |
0.002309 |
0.002125 |
1.146447 |
0.014036 |
0.015552 |
0.05332 |
30.83867 |
0.376638 |
6. |
0.015477 |
0.050035 |
23.77126 |
0.372591 |
0.001725 |
0.002364 |
0.857013 |
0.029163 |
0.015381 |
0.049979 |
23.7558 |
0.371448 |
7. |
0.014881 |
0.046627 |
24.06435 |
0.407267 |
0.001465 |
0.001589 |
0.94963 |
0.017936 |
0.014808 |
0.0466 |
24.04561 |
0.406872 |
8. |
0.018282 |
0.039326 |
19.34072 |
0.496378 |
0.001274 |
0.001096 |
0.752982 |
0.056239 |
0.018238 |
0.03931 |
19.32605 |
0.493182 |
9. |
0.021681 |
0.032801 |
18.91825 |
0.598222 |
0.00099 |
0.001136 |
0.929057 |
0.059913 |
0.021658 |
0.032782 |
18.89542 |
0.595214 |
10. |
0.029286 |
0.020643 |
13.73426 |
0.694877 |
0.001063 |
0.001232 |
0.21893 |
0.002969 |
0.029267 |
0.020606 |
13.73251 |
0.694871 |
11. |
0.033772 |
0.007586 |
14.19678 |
24.73917 |
0.001444 |
0.001227 |
0.258025 |
0.71336 |
0.033742 |
0.007486 |
14.19443 |
24.72888 |
12. |
0.017654 |
0.006266 |
11.8203 |
63.62971 |
0.001545 |
0.001251 |
0.405145 |
0.004067 |
0.017587 |
0.00614 |
11.81335 |
63.62971 |
13. |
0.012169 |
0.050297 |
7.60993 |
63.35626 |
0.001445 |
0.001393 |
0.243327 |
0.008804 |
0.012083 |
0.050277 |
7.606039 |
63.35626 |
14. |
0.00735 |
0.035498 |
7.592265 |
63.17761 |
0.001478 |
0.001705 |
0.792096 |
0.003412 |
0.0072 |
0.035457 |
7.550833 |
63.17761 |
Table 2. Accuracy analysis results
表2. 精度分析结果
误差类型 |
总误差 |
随机误差 |
系统误差 |
|
|
|
|
|
|
|
|
|
|
|
|
数值 |
0.017467 |
0.059715 |
22.0662 |
15.65934 |
0.001785 |
0.001901 |
0.714134 |
0.082273 |
0.017347 |
0.05967 |
22.05229 |
15.65603 |
4.1.2. 估算算法应用验证
为了确保误差参数解算的精度,要求尽量多的数据应用于参数解算,其中测量数据应覆盖方位0˚~360˚空域,仰角在保精度范围内:以20˚~65˚的范围为佳,采用8组满足上述要求的跟踪测量数据进行精度优化分析。
所有测量数据经过了合理性检验、跨零跳点修正及电波折射修正的处理,再以该跟踪时刻对应的天平一号定轨轨道数据进行误差分量估算,得到简化模型中的各误差分量结果如表3所示。
Table 3. Estimation results of error components
表3. 误差分量估算结果
分量名称 |
估算值 |
分量名称 |
估算值 |
方位零值A0 |
−0.017635 |
方位俯仰不正交角δm |
0.00212 |
俯仰零值E0 |
−0.209912 |
重力下垂ρg |
−0.0037 |
距离零值R0 |
15.00322 |
方位机电偏差Vaz |
−0.065 |
大盘不水平方位角Amax |
65.035 |
俯仰机电偏差Vel |
0 |
大盘不水平度θm |
0.010336 |
|
|
备注:仰角E0,Vel进行了合并处理,故给出Vel为零。
将估算的结果装填进修正模型,对未进行系统误差修正的测量数据进行系统误差修正后,一并与运用常规标定的参数进行修正后的测量值同时与天平一号定轨轨道数据进行作差比对,各分量的残差进行分析,结果表明,运用估算的误差分量进行修正后,测量结果基本与常规标定的参数修正基本一致。
再对修正后的数据进行精度分析,结果如表4所示。对比表2的数据结果,可以看出修正后的数据精度也与人工标定的参数修正结果基本一致。
Table 4. Results of accuracy analysis of measurement data after conventional calibration and parameter correction for estimation
表4. 常规标定与估算参数修正后的测量数据精度分析结果
误差类型 |
总误差 |
随机误差 |
系统误差 |
修正类型 |
|
|
|
|
|
|
|
|
|
常规参数 |
0.014705 |
0.186243 |
37.07957 |
0.00274 |
0.004133 |
0.834339 |
0.014447 |
0.186197 |
37.07018 |
估算参数 |
0.013955 |
0.206532 |
32.43043 |
0.00253 |
0.005641 |
1.021215 |
0.011425 |
0.200981 |
31.409215 |
4.2. 系统功能实现
4.2.1. 任务规划功能实现
任务规划功能中,用户可选择需要跟踪的卫星类别,以查看跟踪规划计划。实现所有用户近三日的任务规划情况视图,如图6所示。同时计划视图页面支持放大、缩小、平移等功能,放大后,可以将光标置于弧段前后时显示时间点,便于用户获知精准的时间,而点击空间规划按键,在各时间规划的图条后方给出了空间范围。用户还可以根据需要点击某设备的标签按钮,取消查看该设备规划情况,选取需要查看的特定用户的规划情况,如图7所示。在需要确定实施跟踪任务计划时,系统实现了当前用户所有规划弧段的引导弹道下载,供用户开展鉴定任务或跟踪训练任务的引导文件使用。
Figure 6. Task planning view
图6. 任务规划视图
Figure 7. Enlarged user planning view
图7. 放大后的用户规划视图
4.2.2. 精度鉴定功能实现
在精度鉴定页面,有提供用户选择的卫星类型、轨道产品类型及电波折射修正(大气修正),以便用户根据自己的跟踪情况及对比轨道情况选择用于精度鉴定的卫星类型及轨道产品。卫星类型的选择依据用户所跟踪的实际卫星而定,卫星类型分为TP-1A星和TP-1B星;轨道产品选择则根据用户对精度鉴定的比对标准数据的精度要求而定,分为:准实时定轨轨道产品、超快速定轨轨道产品、快速定轨轨道产品及最终定轨轨道产品4类,主要区别为定轨精度不同。
上传经过修正后的测量数据,根据弹道文件绘制的R、A、E图像,同时系统开始自动获取轨道数据,并对两组数据进行处理后,给用户提供A、E、R、V分别的总误差、随机误差、系统误差的精度鉴定结果及数据下载。如图8所示。
Figure 8. Analysis of total error in accuracy appraisal
图8. 精度鉴定的总误差分析
4.3.3. 精度优化功能实现
同样,在精度优化功能实现上,也具备卫星类型及轨道产品类型的选择,文件上传支持单个及多个测量数据文件同时上传。完成上传文件后,系统对文件进行处理,并根据处理流程进行误差参数估算。并将估算结果显示于用户,如图9所示。估算参数包括:方位零值、俯仰零值、距离零值、大盘不水平方位角、大盘不水平度、方位俯仰不正交度、重力下垂、方位机电偏差、俯仰机电偏差。
Figure 9. The revised analysis results
图9. 修正分析结果显示
4.3.4. 精度监测功能实现
在精度监测功能实现上,主要通过曲线和数据直接展示用户选定范围内的精度鉴定结果,如图10所示。同时,可以对绘制的曲线图片及表格中的精度结果数据进行下载,以供用户保存使用。
Figure 10. Accuracy monitoring effect
图10. 精度监测效果
5. 结束语
本文主要设计并实现了为执行航天发射任务的测量设备开展线上服务模式的精度鉴定的系统,包括单套设备的或多套设备的鉴定任务规划、测量设备的精度鉴定、精度优化以及精度监测。规范了设备精度鉴定的各个环节,避免了繁琐的沟通协调工作,省去了校飞飞机的飞行任务,节省了大量的人力物力财力;实现了实时快速规划出单套设备或多套设备联合的鉴定任务,并提供引导弹道数据下载;实现了随时随地由单个设备即可组织本设备的精度鉴定工作,并快速给出真实有效的精度鉴定结果,大幅提升了精度鉴定的工作效率;实现了设备测量数据精度优化建议,有效提高测量设备的跟踪精度;实现了通过积累各设备精度鉴定结果数据,长期有效地监测设备的精度变化,为各测量设备更好地掌握其精度指标提供数据支撑。
NOTES
*通讯作者。