陆态网络GNSS基准站阶跃偏移量的精确估计.pdfVIP

第 9卷 增刊 震灾防御技术 Vo1．9，Suppl 2014年 l0月 TechnologyforEarthquakeDisasterPrevention Oct．，2014 梁洪宝，朱爽，王友，顾焕杰 ，2014．陆态网络GNSS基准站阶跃偏移量的精确估计．震灾防御技术，9(增刊)：622—631 doi：10．11899／zzfy2014s108 陆态网络 GNSS基准站阶跃 偏移量的精确估计’ 梁洪宝 朱 爽 王 友 顾焕杰 (中国地震局第一监测中心，天津 300180) 摘要 通过对发生阶跃的GNSS基准站坐标序列谱指数的求解，表明了其具有有色噪声的特性，并对 三分量分别建立了最优噪声模型，以此对 GNSS基准站坐标的阶跃偏移量进行了精确估计，探明了地 震和仪器更换对测站产生的精细影响。同时利用求解的阶跃偏移量对坐标序列进行修正，对获取连续 性变化的测站坐标序列和获取精确的测站三维运动参数具有重要意义。 关键词：陆态网络 GNSS基站站 阶跃 偏移量 引言 在我国，2007年国家重大科技基础设施中国大陆构造环境监测网络 (简称 “陆态网络”) 在 1997年国家重大科学工程中国地壳运动观测网络 (简称 “网络工程”)的基础上，建设了260 个GNSS基准站，2011年全面建设完成并投入运行 (李强等，2012；牛之俊等，2002)。陆态 网络运行过程 中，GNSS测站经常因远场大地震引起的同震位移、仪器或天线更换、观测墩维 护或改建和一些未知等原因，使得GNSS坐标位置存在阶跃偏移，如图1所示。其中甘肃岷县 (GsMx)站时问序列发生阶跃是由2013年 7月22日甘肃岷县6．6级地震引起的同震位移；而 山东荣成 (SDRC)站时间序列发生两处阶跃偏移，分别是2011年 3月 11日日本9．0级特大 地震引起的同震位移和2013年 11月3日因仪器天线升级引起的阶跃偏移。因此，由包含阶跃 偏移的时间序列得到的测站运动参数是不正确的，必须予以剔除。通常情况下，假定GNSS坐 标时间序列是纯 白色噪声，采用数值拟合方法来求解阶跃偏移量，但一系列研究成果证明， GNSS坐标时间序列具有有色噪声特性。研究成果表明，GNSS坐标时间序列的噪声更接近 白 噪声和闪烁噪声的叠加 (Nikolaidis，2002a；黄立人等，2006~2007)，因此，仅基于白色噪声 模型下的阶跃偏移量求解是不准确的，甚至可能导致一个数量级的过高估计。为此，本文选取 陆态网 1月间坐标位置明显发生阶跃偏移的52个测站，测站分布如图2所示， 对这 52个测站坐标三分量分别进行噪声分析，并对三分量分别建立相应的最优噪声模型，在 1基金资助 中国地震局 “三结合”课题 (143301)；中国地震局震情跟踪课题 (2014020205) [收稿 日期] 2014．05—27 [作者简介] 梁洪宝，男，生于 1984年。助理工程师。主要从事GNSS区域监测、数据处理与形变分析。E—mail：lhbl31421 @126．c0m 624 震灾防御技术 9卷 1．2 粗差剔除 根据 GNSS测站坐标时间序列的特点，粗差的剔除通常可用 “拉依达准则”(常被称为 3 倍中误差准则)。考虑到时间序列中可能存在周期性和线性变化，本文采用相邻两个观测值之 差来剔除粗差，具体做法如下 (黄立人等，2012)： 设序列中第 i个观测值为 ，对应的误差估计为m，相邻的后一个观测值为 小 对应 的误差估计为m川，但第一个观测值本身不是粗差，可用下式进行判别： IHl=1Hi+1一Hl3×√( + 2+『1) (1) 如果上式成立，则第 i+1个观测值为粗差，予以剔除，这样剔除了第 i+1个观测值后原 来的第 i+2个观测值就变成了第 i+1个观测值，以此类推，并继续进行后续序列中粗差的剔 除。在上述算法的基