元亚菲，张以文，孔庆珠，李朋朋，李志远，张雨阳

（1.中国矿业大学 环境与测绘学院，江苏 徐州 221116；2.水发规划设计有限公司，山东 济南 250100）

0 引 言

近几十年来，高精度、高时间分辨率的GNSS区域网建设取得了长足发展，为研究板块运动机制、建立和维护地球参考框架、揭示大地构造形变过程等提供了重要的数据支撑。大量研究结果表明，由于GNSS站易受外部环境、自身观测技术和地球构造运动等影响，GNSS站坐标时间序列存在明显的有色噪声［1-2］，且不同方向上的有色噪声存在一定相关性［3］。因此，研究GNSS站坐标时间序列的最优噪声模型，确定不同方向噪声的耦合关系，可以获取更为准确的运动特征，进而更加合理、全面地进行地球物理现象研究。针对GNSS站坐标时间序列的最优噪声模型和运动特征，国内外众多学者进行了深入分析，得到了一系列有益成果。如：A.Mao等［4］对全球23个GPS站三年的坐标序列进行分析，认为白噪声（white noise，WN）和闪烁噪声（flicker noise，FN）的模型组合可以作为全球GPS站噪声模型的最佳选择；S.D.P.Williams等［5］利用极大似然估计法（MLE）对全球954个GPS站坐标序列进行分析，证明了WN+FN能更好地表示坐标序列的噪声模型；明锋等［6］利用STL法提取了中国区域IGS站的长时间变化趋势，并利用MLE法估计了各站点的最优噪声模型；马俊等［7］采用最小范数二次无偏估计法（MINQUE）分析了中国区域20个IGS站的噪声分布情况；He X等［8］采用MLE法估计了全球110个测站的最优噪声模型，并分析低频噪声对趋势估值的长期影响。

已有GNSS站坐标时间序列研究中，噪声模型估计方法主要为极大似然估计法和方差协方差分量估计法（VCE）。……