杨登科

(铁道第三勘察设计院集团有限公司,天津 300142)

欧洲地区IGS基准站坐标时间序列各分量方向最优噪声模型变化探讨

杨登科

(铁道第三勘察设计院集团有限公司,天津 300142)

针对现有坐标时间序列研究中各分量方向间噪声特性对比研究较少,本文选取欧洲地区58个IGS基准站11年的坐标时间序列,采用不同的噪声模型组合对它们进行噪声分析并就各分量方向间的噪声特性进行对比。结果表明,各测站在NEU分量表现出不同的噪声特性,因此分析时间序列变化特征的物理解释时需对各分量区别对待;垂直分量的有色噪声、振幅等明显大于水平方向,与一般认为的垂直分量精度低于水平分量相一致。

坐标时间序列;最优噪声模型;分量方向;极大似然估计;噪声特性

0 引 言

20多年累积的IGS基准站坐标时间序列序列,已成为大地测量学和地球动力学研究(地壳形变监测、全球板块动力学研究、区域构造运动等)的主要数据之一。针对GPS坐标时间序列建立适当的噪声模型,可实现形变信号和噪声的有效分离,其结果可精化速度场,为板块运动引起的线性构造形变以及各种非线性形变的分离提供准确的基础数据[1]。

当前,国内外学者对GPS坐标时间序列各分量方向的噪声特性进行了研究。如文献[2]中李昭对中国区域11个IGS基准站的坐标时间序列数据进行了噪声特性分析,并指出中国区域IGS基准站各分量的噪声特性存在多样性;文献[3]中蒋志浩等利用CGCS2000坐标系下的29个CORS站坐标时间序列进行了噪声分析,得出我国CORS站的东西分量比南北分量具有较大的闪烁噪声,并提出这主要与我国周边IGS站框架点分布有关;文献[4]中李靖分析了“陆态网络”基准站坐标时间序列的噪声特性,指出在水平方向,E分量比N分量具有较大的随机漫步噪声和闪烁噪声。由于对各分量方向间的噪声特性缺乏系统详尽的对比分析,因此采用不同的噪声模型组合对欧洲地区58个IGS基准站坐标时间序列数据进行噪声分析,并就各分量方向间的噪声特性进行对比是极为必要的。

本文的主要目的是采用多种类的噪声模型组合分析欧洲地区58个IGS基准站坐标时间序列的随机特性,并对各分量方向内的噪声特性进行对比分析,从而得出有意义的结论。

1 数据及分析策略

1.1数据介绍

使用SOPAC提供的欧洲地区58个IGS基准站2001.0014-2011.9986的坐标时间序列数据,地理位置如图1所示(由于站数较多,未标注IGS基准站站名)。

1.2数据处理策略

能有效确定GPS坐标时间序列噪声特性的方法通常有两种:极大似然估计和功率谱分析,分别从时间域和频率域对GPS坐标时间序列进行噪声分析。极大似然估计可避开频谱分析需要均匀采样、依赖频谱平均的局限性,被认为是目前最准确的噪声分析方法。本文选取WN、FN+WN、RWN+WN、FN+RWN+WN、PL+WN、FOGM+WN、FOGM+RWN+WN共七个噪声模型,采用CATS软件对IGS基准站各分量方向进行噪声分析[5]。

1.3最优噪声模型评价准则

选用Langbein提出的保守准则判断不同模型的优劣,具体的最优噪声模型评价准则可参见文献[1]。

2 GPS坐标时间序列最优噪声模型的建立

2.1谱指数的求取

谱指数计算是一种有效判断噪声模型的方法,利用CATS软件求取各测站在不同分量方向的谱指数,由于站数较多,本文仅列出5个测站的谱指数(下同),结果如表1所示。

表1 各分量方向的谱指数

从表1可以看出,各分量方向的谱指数均不相同,且均不等于0,表明各分量方向均包含有色噪声。对欧洲地区58个IGS基准站各分量方向谱指数的平均值进行了统计,如图2所示。

图2 N、E、U三方向谱指数对比

由图2可见,N、E、U三方向谱指数平均值差别不大,N方向谱指数平均值最大。

2.2GPS时间序列最优模型的建立

采用1.2节的数据处理策略和1.3节的最优噪声模型评价准则,对选取的欧洲地区58个IGS基准站坐标时间序列各分量方向进行上述7种组合噪声模型分析,得到各分量方向的最优噪声模型[6-8],结果如表2所示。

表2 各分量方向的最优噪声模型

由表2可知,各测站在不同分量方向的最优噪声模型是可变的,时间序列在一个分量方向求得的最优噪声模型仅能代表此方向,因此在对时间序列进行噪声模型等的研究时,需指出相应方向。对欧洲地区58个IGS基准站各分量方向的最优噪声模型进行统计,结果如图3所示。

图3 N、E、U三方向最优噪声模型对比

由图3可知,N、E、U三方向最优噪声模型主要为FN+WN和PL+WN,其余噪声模型组合所占比例都较小;各测站的最优噪声模型是可变的,且N、E、U三方向表现出不同的最优噪声模型变化特征。

2.3FN+RW+WN噪声分量对比

对FN+RW+WN中各噪声分量在每个方向所占比值进行了统计,如图4所示。

图4 N、E、U三方向噪声分量对比

由上图可知,闪烁噪声为噪声的主要成分,在N、E、U三方向所占比例均达到了50%;垂直分量有色噪声(FN和RW均为有色噪声)明显大于水平分量,这与一般认为的垂直分量精度低于水平分量相一致。

2.4振幅的求取

图5 N、E、U三方向振幅对比

由图5可知,N、E方向振幅大致相同,且远小于垂直方向。

3 结束语

随着GPS的快速发展,针对GPS坐标时间序列建立适当的噪声模型,实现形变信号与噪声的有效分离,是GPS数据处理领域的热点问题。本文利用欧洲地区58个IGS基准站坐标时间序列数据,采用多种类的噪声模型组合对其进行噪声分析并对比各分量方向间的噪声特性。结果表明:各测站在N、E、U分量方向表现出不同的噪声特性,因此分析时间序列变化特征的物理解释时需对各分量区别对待;垂直分量的有色噪声、振幅等明显大于水平方向,这与一般认为的垂直分量精度低于水平分量相一致。

[1] 杨登科,安向东,黄广利,等. 不同GPS时间序列跨度对噪声模型建立的影响[J]. 测绘科学,2016(5):33-37.

[2] 李昭. GPS坐标时间序列的非线性变化研究[D]. 武汉:武汉大学,2012:1-178.

[3] 蒋志浩,张鹏,秘金钟,等. 顾及有色噪声影响的CGCS2000下我国CORS站速度估计[J]. 测绘学报,2010,39(4):355-363.

[4] 李靖. “陆态网络”基准站坐标时间序列变化特性分析[D]. 郑州:解放军信息工程大学,2013.

[5] 杨登科,邓连生,安向东,等. IGS基准站坐标时间序列最优噪声模型变化探讨[J]. 测绘地理信息,2016(1):7-10.

[6] 张彦芬,毕刚,陈华,等. 山西CORS坐标时间序列噪声特征分析[J]. 大地测量与地球动力学,2012,32(4):61-66.

[7] 田云锋. GPS位置时间序列中的中长期误差研究[D]. 北京:中国地震局地质研究所,2011:1-124.

[8] 苏利娜,丁晓光,张彦芬,等. 陕西连续GPS基准站坐标时间序列分析[J]. 大地测量与地球动力学,2014,34(5):106-110.

TheAnalysisontheBestNoiseModelChangeoftheComponentDirectionofIGSReferenceStationCoordinateTimeSeriesinEurope

YANGDengke

(The3rdRailwaySurveyamp;DesignInstitute,Tianjin300142,China)

According to the fact that there is on enough research on the comparison of noise characteristics among different components, this paper focusing on the coordinate time series of 58 reference stations in Europe for 11 years, noise analysis has been implemented using different noise model combination and the noise characteristics of each component direction are compared. From these analysis, the noise characteristics in the NEU component of each station is variable, so it is necessary to treat each component differently when analysing the physical interpretation of time series variation characteristics; the coloured noise and amplitude of the vertical component are obviously larger than the horizontal direction, which is consistent with the generally considered the accuracy of the vertical component lower than the horizontal component.

Position time series; optimal noise model; component direction; maximum likelihood estimation; noise characteristics

10.13442/j.gnss.1008-9268.2017.05.004

P228.4

A

1008-9268(2017)05-0021-04

2017-05-30

联系人: 杨登科 E-mail: 624892261@qq.com

杨登科(1990-),男,河南商丘人,硕士，助理工程师,主要从事GPS坐标时间序列分析及应用、变形监测等研究。