空间滤波对国家cors坐标时间序列的影响.docxVIP

空间滤波对国家cors坐标时间序列的影响

1cors站坐标时间序列噪声定

国家导航卫星系统的连续运营指导站（cois）是“国家空间数据基础设施”中最重要的组成部分。这是中国制定和维护高精度和动态地心坐标模式的唯一途径。国家CORS站运动速度估值和误差估计是维持CGCS2000(中国大地坐标系统2000)框架精确性、动态性和稳定性的重要基础。目前在估计国内CORS站运动速度参数时,通常由CORS站日解坐标值组成CORS站原始坐标时间序列,在假定原始坐标时间序列误差仅具有白噪声性质的前提下,采用线性拟合方法得到CORS站运动速度估值和误差估计。

随着CORS网坐标时间序列资料的不断积累,CORS站坐标时间序列中所包含的噪声的性质及其定量估计问题也逐渐引起了人们的重视,国内外一些学者采用了不同噪声模型初步分析了CORS站坐标时间序列噪声性质,研究表明CORS站坐标时间序列中同时包含有白噪声和有色噪声的影响,而且白噪声不是噪声的主要成分。但是上述学者的研究只是分析给出了CORS网坐标时间序列的短期成果和局部区域内的有色噪声的数值,没有定量分析在CGCS2000框架下有色噪声对我国CORS站速度估值和误差估计的整体性影响。本文首先采用线性拟合的方式去除国家CORS网坐标时间序列中的地壳构造形变信息,并采用主成分空间滤波方法(principalcomponentanalysis,PCA,简称空间滤波)提取CGCS2000框架下国家CORS网坐标时间序列中公共误差(commonmodeerrors,CME)的时空特性,提高CORS站本地信号(坐标时间序列)的信噪比,并比较分析空间滤波前后CORS站成果的精确性。其次对空间滤波后的国家CORS站本地周期性信号进行三角函数拟合,精确地提取国家CORS站本地周期性变化信息和本地噪声信息(坐标残差时间序列),采用功率谱方法分析公共误差序列的频谱特性和国家CORS站坐标残差时间序列的噪声性质,在此基础上采用白噪声、闪烁噪声和随机漫步噪声组合的误差模型,利用最大似然估计方法定量估计了坐标残差时间序列中的有色噪声分量。最后采用加权最小二乘法评定了顾及不同噪声影响的CGCS2000框架下的国家CORS网年速度估值和实际精度。

2glodk软件平差处理

对国家CORS网的数据处理采用三个步骤。第一步,利用GAMIT软件获得国家CORS网点坐标和卫星轨道的单日松弛解;第二步,利用GLOBK软件将获得的国家CORS网单日松弛解结果和全球IGS连续运行站的松弛解合并,得到一个包含国内CORS网和全球IGS(internationalGPSservice)站的单日松弛解结果;第三步,利用GLOBK软件平差处理上一步合并后的单日松弛解,并选取全球IGS站在ITRF1997框架的坐标和速度值作为约束(选取方法和结果见参考文献),采用相似变换方法求解出国家CORS网成果在CGCS2000坐标框架中的单日解坐标值。采用上述方法对近十年CORS站数据进行处理,获得了国家CORS站29站的坐标时间序列,CORS站坐标时间序列水平分量的重复性(加权均方根误差)优于5mm,垂直分量的重复性优于12mm,达到了较高的精度。

3时间序列中ti1

国家CORS站本地坐标时间序列拟合函数模型为

y(ti)=a+bti+csin(2πti)+dcos(2πti)+esin(4πti)+fcos(4πti)+j0∑j=1ΟjΗ(ti-Τj)+vi(1)y(ti)=a+bti+csin(2πti)+dcos(2πti)+esin(4πti)+fcos(4πti)+∑j=1j0OjH(ti?Tj)+vi(1)

式中,ti是以年为单位的时间序列中解算历元;a为常数项;b为线性速度;c和d组合表示全年性周期运动;e和f组合表示半年性周期运动;Oj为在Tj时刻CORS站的位置突变,H(ti-Tj)表示阶梯函数,j0∑j=1ΟjΗ(ti-Τj)∑j=1j0OjH(ti?Tj)表示发生在历元Tj处大小为Oj的j0个偏移常量;vi代表残差序列。

为了消除国家CORS网坐标时间序列中的点位突变项对速度估值的影响(如2001年11月14日昆仑山口西8.1级地震对DLHA站的影响,2008年5月12日汶川8.0级地震对泸州站的影响等),通过建立阶梯函数(公式(1)中H),将点位突变项与点位速度同时估计。

3.1cors站的公共误差

我国CORS站在CGCS2000框架下的原始坐标时间序列具有高度的相关性,即在不同站点坐标时间序列中包含时空相关的公共误差,公共误差是国家CORS网数据处理的主要误差源,公共误差的存在使得我国CORS站的