论文部分内容阅读
地球定向参数(EOP)是实现天球参考架与地球参考架相互转换的参数，本文在EOP高频变化的VLBI资料检测方法、实测资料归算和海潮激发方面进行了有益的尝试。主要工作与结论如下： (1) 对地球定向参数高频变化的研究现状进行了调研，包括基于VLBI和GPS实测资料对EOP高频变化的测定以及对高频变化的机理解释等内容。主要结论为，虽然关于周日和亚周日地球自转高频变化的研究取得了较大进展，但也仍然存在很多亟待进一步分析的问题。例如，不同作者基于实测资料得到的EOP高频变化潮汐项估值存在差异的原因是什么，如何从实测资料中得到更为可靠的EOP高频变化解算序列，岁差章动模型的缺陷对求解EOP高频变化有何影响等等。对EOP高频变化的求解和机理解释仍需要进一步深化。 (2) 在CALC／SOVLE软件系统的相关分析和模型改进方面进行了以下4个方面的工作。 (a) 介绍了稀疏矩阵的四种常见形式以及稀疏矩阵技术在天测与测地VLBI数据处理中的应用。推演了天测与测地VLBI数据综合解算中所用稀疏矩阵形式下待估参数求解和协方差矩阵估算的算法。通过对是否采用稀疏矩阵技术时方程求解(乘法和加法)运算对数的估算和比较，表明普通最小二乘方法的运算对数约为参数总数的3次方，而采用稀疏矩阵技术时的运算对数近似与参数总数成线性关系，从而能够在现代空间大地测量技术的大样本数据处理中显著缩短计算时间、降低资源需求。 (b) 对上海佘山站参与的天测与测地VLBI实验抽样选取27次，分别进行了单次解算。分析了解算参数随剩余钟行为和大气效应分段拟合长度的变化，结果表明，选取不同分段拟合长度时，站坐标解算结果和时延残差加权均方根存在差异，最大分别至厘米级和数十皮秒，因而分段拟合长度不是随意选取的。分段拟合时段长度不宜过长，否则钟和大气剩余效应的短周期变化不能被很好地模型化。同时，为了保证待估参数解算时有足够的自由度，拟合时段不宜过短，否则将导致法方程近于或出现奇异，也达不到较好控制噪声的效果。因而分段拟合长度存在某一合理取值范围。由于各次实验以及同一次实验中的不同台站其钟和大气条件各异，因而有必要对每次实验以及每次实验中的各观测台站分别分析寻找合适的分段拟合长度。这在实际操作中显然是相当烦琐的工作。一般而言，在剩余钟行为拟合长度(CALC／SOLVE)缺省值60min情况下，剩余大气效应分段拟合长度介于10min至40min为宜，在剩余大气效应分段拟合长度缺省值20min情况下，剩余钟行为分段拟合长度介于20min至100min为宜。 (c) 上海天文台运用CALC／SOLVE系统进行全球全时段天测与测地VLBI观测资料综合归算工作，前后已持续10余年的时间。其天球参考架、地固参考架和地球定向参数序列归算结果，曾多次被IERS、IVS等相关国际学术组织采纳为综合结果的输入解，特别是2000年的地固参考