Using Orthogonal Transformation Method to Select Independent Baselines and Independent Double-difference Ambiguities
-
摘要: 介绍了利用Gram-Schmidt正交化算法构建独立向量集的原理,并利用关联矩阵和向量权的思想构建独立基线向量集和独立双差模糊度向量集,讨论了三种构建双差模糊度向量的特点,给出了选择独立基线和独立双差模糊度的正交化算法。通过一个全球GPS网定轨算例,给出了计算结果,验证了方法的有效性。
-
关键词:
- Gram-Schmidt正交化 /
- 独立基线 /
- 独立双差模糊度 /
- 关联矩阵
Abstract: The principle of constructing independent vector sets using Gram-Schmidt orthogonal transformation is proposed. The method of setting independent baseline and independent double-difference ambiguity data sets refers to the idea of vertex incidence matrix and vector weight. The three ways to build double-difference ambiguity data sets are discussed. The computing process of selecting independent baselines and independent double-difference ambiguities is provided. An example of orbit determination using global GPS track nets is tested. The result shows that the orthogonal transform algorithm is effective. -
近年来，随着各类GPS监测网络的建立，尤其是中国地壳运动观测网络(简称“网络工程”)和中国大陆构造环境监测网络(简称“陆态网络”)的投入运行以来，产生了大量的GPS连续和流动观测数据。依据这些数据得到的GPS速度场及其衍生产品在地壳运动监测、地球动力学研究等领域发挥了重要作用[1-7]。在上述研究中，GPS流动水平速度场应用较多，而GPS流动垂向速度场受非构造运动影响较大且缺乏足够的观测数据来消除这些影响，使得GPS流动观测解算的垂向速度场精度不高，应用较少。
目前，我国部分区域(如京津唐、川滇地区)GPS连续观测网络已初具规模，如何利用这些较高密度的连续观测网络提高同一区域的流动观测垂向速率精度，是挖掘GPS垂向速度场价值的当务之急，具有重要的意义。鉴于此，本文提出了利用GPS连续观测修正流动观测垂向速率的方法，并选用云南地区陆态网络GPS连续和流动观测数据进行了验证。
1 GPS连续观测年周期信号提取与应用
GPS垂向时间序列存在较为显著的季节性运动，且年周期运动最为突出[8-9]，因此，本文只考虑年周期运动，则有：
(1) 式中，Y为垂向时间序列；a为常数项；V0为速率；A和B为年周期运动正弦波和余弦波振幅；C为年周期运动总振幅；φ为年周期运动初相；ε为残差；t为时间。
由于网络工程和陆态网络采用的流动观测模式为每隔2~3 a观测一期，每期连续观测4 d，这种观测模式难以获取准确的年周期运动信息，因此,求解时只能顾及线性运动：
(2) 在用最小二乘理论求解参数时，如果忽略部分参数，则估计参数的平差值会包含未估计参数的影响[8]，即按式(2)求得的速率V1包含了年周期运动的影响，与V0存在一定偏差。如果能够获取每个流动站的年周期运动振幅并在时间序列Y中扣除年周期运动得到Y′，再根据Y′只顾及线性运动求解速率V2，则V2应该能够扣除年周期运动的影响，等于(或者接近于)V0：
(3) 由于一些区域内GPS时间序列中的非构造噪声存在区域相关性[10-11]，因此，可以利用同一个区域的连续站年周期运动信息通过空间加权内插得到流动站的年周期运动信息。对于每个流动站，年周期运动振幅A0和B0分别取其周边200 km内所有连续站的年周期运动振幅A和B的加权平均(权取距离d的平方倒数)：
(4) 2 连续站数据检核
2.1 数据概述
为了检核修正方法，选取陆态网络GPS连续站密度较高的云南地区进行试验。陆态网络在该地区布设了27个连续站，其中KMIN和XIAG站从1999年开始观测，其他25个测站是从2010年陆续开始观测。本文选取这27个测站2010-01-01以后的观测数据，至2013年底这些测站观测时间均超过960 d，能够较为准确地获取该地区垂向速率和年周期运动信息。利用GAMIT/GLOBK[12]和QOCA[8]软件解算各连续站观测数据，得到各连续站2010~2013年的时间序列，然后按式(1)分别求解各连续站垂向运动速率和年周期运动振幅(表 1和图 1)。
![]()
图 1 云南地区GPS连续站年周期运动 (2010~2013年，参考点2010.0)Figure 1. Annual Terms of Vertical Position Variations for GPS Fiducial Stations(from 2010 to 2013, the Reference Time as 2010.0) in Yunnan Area表 1 云南地区GPS连续站垂向运动信息(2010~2013年)Table 1. Vertical Motion for GPS Fiducial Station in Yunnan Area from 2010 to 2013点名 速率 /(mm\5a-1) 年周期 观测天数 振幅/mm 初相/(°) KMIN 2.93 9.45 17.61 1 305 XIAG 0.17 9.5 357.3 1 365 YNCX* 1.52 8.55 4.11 1 108 YNDC* 2.21 9.44 9.89 1 050 YNGM -0.24 6.35 350.08 988 YNHZ 1.21 8.29 8.56 973 YNJD 1.04 10.97 8.31 1126 YNJP -1.9 8.88 17.38 1 134 YNLA -0.68 14.31 9.05 1 147 YNLC* 0.6 11.49 5.49 1 008 YNLJ* 1.39 10.56 18.84 1 115 YNMH -1.02 11.6 7.86 1 150 YNMJ -0.3 9.89 0.15 1 129 YNML 1.07 7 12.25 1 149 YNMZ* -3.17 10.65 24.8 1 163 YNRL -0.64 11.8 10.57 1 087 YNSD* 2.18 12.27 11.65 975 YNSM* -1.54 10.94 5.44 1 162 YNTC 0.84 12.73 6.87 1 149 YNTH* 1.16 8.27 17.31 1 164 YNWS 3.17 8.47 21.88 968 YNXP 2.15 7.95 6.63 1 215 YNYA 2.38 9.68 14.57 1 006 YNYL 1.95 12.85 16.35 1 163 YNYM* 1.45 10.22 17.14 1 121 YNYS 2.11 11.11 12.93 1 160 YNZD 0.86 8.97 22.98 980 注: 标*的点为检核站，未标注的点为基准站。 表 1是云南地区GPS连续站2010~2013年垂向运动信息，从中可以看出该地区各GPS连续站年周期运动较为显著，年周期运动振幅较大(7~14 mm)，垂向运动速率相对较小(4 mm/a以内)，因此，测站年周期运动会对垂向运动速率的求解产生较大的影响，在该地区有必要对流动观测垂向运动进行修正。
图 1是云南地区GPS连续站2010~2013年年周期运动，时间参考点为2010.0，箭头长度表示振幅,从正东逆时针旋转的方位角表示初相位方向，箭头指向东、南、西、北方向依次表示极值发生在0.25、0.50、0.75、0.0 年。结合表 1和图 1可以看出云南地区GPS连续站的年周期运动空间一致性较好，各连续站年周期运动振幅和初相均比较一致，满足本方法的使用条件。
2.2 检核结果及分析
在云南地区27个连续站中选取18个测站(图 1中蓝色点)记为基准站，9个测站(表 1标*点和图 1中黑色点)记为检核站，用以检核修正方法效果。首先根据每个连续站垂向时间序列Y按式(1)用最小二乘方法求解速率V0和年周期运动振幅A和B；然后根据每个检核站垂向时间序列Y按式(2)用最小二乘方法求解速率V1；再次，每个检核站按式(4)根据基准站的年周期运动振幅插值求得年周期运动振幅A0和B0，并在时间序列Y中扣除依据A0和B0计算得到的年周期运动，得到时间序列Y′，然后按式(3)用最小二乘方法求解速率V2。比较V1和V2与V0的差异(见表 2)，如果V2与V0的差异比V1与V0的差异小，则说明修正方法有效。
表 2 检核站解算结果Table 2. Results of Check Sites检核站 年周期运动振幅/mm 速率结果/(mm\5a-1) A B A0 B0 V0 V1 V2 V1-V0 V2-V0 YNCX 8.53 0.61 9.55 1.66 1.52±0.12 2.04±0.18 1.55±0.12 0.51 0.03 YNDC 9.30 1.62 8.25 1.39 2.21±0.19 1.21±0.25 2.09±0.19 -1.01 -0.12 YNLC 11.44 1.10 8.95 0.23 0.6±0.15 1.21±0.26 0.71±0.16 0.61 0.11 YNLJ 9.99 3.41 10.43 2.38 1.39±0.11 1.29±0.2 1.24±0.1 -0.10 -0.15 YNMZ 9.66 4.47 8.17 2.36 -3.17±0.16 -3.01±0.23 -3.22±0.16 0.16 -0.06 YNSD 12.02 2.48 11.00 1.28 2.18±0.17 -0.22±0.29 1.86±0.17 -2.40 -0.32 YNSM 10.89 1.04 10.81 1.01 -1.54±0.12 -1.14±0.21 -1.54±0.12 0.40 0.00 YNTH 7.90 2.46 8.09 1.63 1.16±0.11 1.34±0.17 1.1±0.12 0.18 -0.05 YNYM 9.76 3.01 9.36 2.10 1.45±0.13 1.31±0.21 1.34±0.13 -0.14 -0.11 从表 2可以看出，不顾及年周期运动求得的速率(V1)与顾及年周期运动求得的速率(V0)差别较大，最大值达到-2.4 mm/a，说明在年周期运动较为显著的地区，在求解线性速率时必须考虑年周期运动。由于该区域年周期运动空间一致性较好，9个检核站内插得到的年周期运动振幅比较接近于实际值(差异在3 mm以内)。由于内插得到的检核站年周期运动振幅精度较高，使得本方法对垂向速率的修正效果较好，各检核站经过年周期运动修正后得到的速率(V2)比未经过修正的速率(V1)更为接近V0，9个检核站V2与V0的差异量值在0.32 mm/a以内，且V2的误差比V1小，接近V0的误差。以上结果说明，在垂向时间序列年周期运动显著且空间一致性较好的区域，通过本方法能够较好地消除年周期运动对垂向速率的影响，提高垂向速率的解算精度。
3 流动站修正结果及分析
陆态网络在云南地区布设了157个GPS流动站，分别于2011年3~8月和2013年4~9月进行了一期观测，每期观测时段长度为4 d。利用GAMIT/GLOBK和QOCA软件解算了这两期流动观测数据(解算模式与§2中GPS连续站时间序列的解算模式保持一致，即GAMIT解算参数、QOCA解算中的框架点的选取等设置一致)，得到了各流动站2011~2013年的时间序列，按式(2)求得了各流动站未扣除年周期运动的2011~2013年垂向速率(图 2(a))。利用§2计算得到的云南地区27个连续站年周期运动振幅，按式(4)内插得到各流动站年周期运动振幅，然后对各流动站的时间序列扣除年周期运动，再按式(3)求得扣除年周期运动的2011~2013年垂向速率(图 2(b))。
![]()
图 2 云南地区GPS垂向速度场(流动观测，2011~2013年)Figure 2. GPS Vertical Velocity Fields in Yunnan Area (Mobile Observations,from 2011 to 2013)由图 2可知，未进行年周期改正的垂向速度场(图 2(a))和改正后的垂向速度场(图 2(b))整体趋势类似，即大体以北纬24°为界，北纬24°以北的区域以隆升为主，北纬24°以南以沉降为主，大体与由GPS连续站测定的云南地区垂向速度场(图 3)特征一致。改正后的垂向速度场比未改正的速度场在空间一致性上有所改善，且更加接近GPS连续站测定的垂向速度场和水准观测结果得到的该区地壳垂直运动特征[13-14]，如图 2中两个红色椭圆框内区域。未改正的速度场显得较为混乱，方向上有升也有降，量值有大有小，改正后的速度场基本上都是上升，且量值空间分布更为均匀。图 2中绿色椭圆框区域的速度场经过改正后更接近郝明等的文章中图 7结果[14]。总体而言，本方法对GPS流动垂向速度场的改进取得了较好的效果。
![]()
图 3 云南地区GPS垂向速度场(连续观测，2010~2013年)Figure 3. GPS Vertical Velocity Fields in Yunnan Area (Continuous Observations,from 2010 to 2013)4 结 语
本文提出了利用连续观测站年周期运动信息通过空间插值获取流动观测站年周期运动信息、并在流动观测站时间序列中加以扣除、进而求解垂向速率的方法，并利用云南地区陆态网络2010~2013年连续观测数据和2011~2013年流动观测数据进行了验证，结果表明在年周期运动空间一致性较好的区域(云南地区)能够取得较好的效果。
本方法的前提是区域内年周期运动空间一致性较好，只有满足这个条件，才能利用连续站年周期运动信息内插得到较为准确的流动站年周期运动振幅，才能得到较为准确的垂向速率。由于本文的主要目的是验证利用连续站信息修正流动观测垂向速率的方法有效性，因此，只考虑了季节性运动中最主要的年周期运动，但本方法也适用于其他季节性运动(如半年周期运动)，只需在式(1)和式(3)增加相关参数求取并加以修正。
我国国土辽阔，各区域情况不尽相同，可先处理区域内连续观测数据，分析区域垂向季节性运动机制。如果与趋势(线性)运动相比，季节性运动振幅较小，对垂向速率的影响可忽略不计，则无需进行季节性改正；如果区域内季节性运动影响较大，且有多种较为显著的周期运动(年周期、半年周期)，则可以根据分析结果选取对垂向速率影响较大的季节项加以修正。为了进一步提高方法的准确性，还将进行空间插值方法的择优选取、在GPS连续站垂向季节性运动信息求取前扣除一些能够定量计算的非构造形变(如大气、积雪等负荷效应造成的非构造形变)[15-18])等分析研究。
致谢: 感谢GFZ葛茂荣研究员在算法上给予的指导。感谢IGS提供的基准站数据。 -
表 1 独立双差模糊度选取情况
Table 1 Selection Statistics of Independent Double-difference Ambiguity
基线号 双差模糊度数 独立双差模糊度数 独立双差模糊度数占比 1 278 45 0.162 2 248 47 0.190 3 288 46 0.160 4 231 48 0.208 5 277 45 0.163 6 319 46 0.144 7 271 50 0.185 8 327 50 0.153 9 457 63 0.138 10 360 53 0.147 … … … … 总计 31 485 5 635 0.179 
下载: 导出CSV
-
[1] Ge Maorong, Gendt G, Dick G, et al. Improving Carrier-Phase Ambiguity Resolution in Global GPS Network Solution[J]. Journal of Geodesy , 2005(4):103-110
[2] 施闯. 大规模、高精度GPS网平差处理与分析理论及其应用[D]. 武汉:武汉大学,1999 Shi Chuang. The Theory and Application of the Large Scale and High Precision GPS Network Adjustment and Analysis[D]. Wuhan, Wuhan University,1999
[3] 王磊,李盼,吕翠仙.关联矩阵法在独立基线及独立双差模糊度选择中的应用[J].武汉大学学报·信息科学版,2010,35(6): 715-718 Wang Lei, Li Pan, Lv Cuixan. Application of Vertex Incidence Matrix Algorithm to Independent Baselines Search and Independent Double-Difference Ambiguity Selection[J]. Geomatics and Information Science of Wuhan University , 2010,35(6): 715-718
[4] 薛定宇. 科学运算语言Matlab5.3程序设计与应用[M].北京:清华大学出版社,2000 Xue Dingyu. Scientific Computing Language Matlab5.3 Program Design and Application[M]. Beijing: Tsinghua University Press, 2000
[5] Dach R, Hugentobler U, Fridez P, et al.Bernese GPS Software Version 5.0[R].Astronomical Institute, University of Bern, 2007
[6] NAPEOS, Mathematical Models and Algorithms[R].European Space Operation Centre,2009
[7] 龚冬保, 魏战线. 大学数学教程(第2卷):线性代数与空间解析几何[M].西安:西安交通大学出版社,2000 Gong Dongbao, Wei Zhanxian. University Mathematics Tutorial (Second Volume): Linear Algebra and Space Analytic Geometry[M]. Xi'an: Xi'an Jiao Tong University Press, 2000
计量
- 文章访问数: 1257
- HTML全文浏览量: 64
- PDF下载量: 329
