测绘学报 ›› 2023, Vol. 52 ›› Issue (8): 1245-1254.doi: 10.11947/j.AGCS.2023.20220675

• 大地测量学与导航 • 上一篇    下一篇

GNSS连续观测站构造运动变化特征自适应提取方法

苏小宁1, 石睿娟1, 鲍庆华1, 朱庆2, 孟国杰3, 闫浩文1   

  1. 1. 兰州交通大学测绘与地理信息学院, 甘肃 兰州 730070;
    2. 西南交通大学地球科学与环境工程学院, 四川 成都 611756;
    3. 中国地震局地震预测研究所, 北京 100036
  • 收稿日期:2022-11-28 修回日期:2023-03-28 发布日期:2023-09-07
  • 通讯作者: 朱庆 E-mail:zhuq66@263.net
  • 作者简介:苏小宁(1986-),男,博士,副教授,博士生导师,主要从事地壳变形特征与变形机理研究。E-mail:suxiaoning_666@mail.lzjtu.cn
  • 基金资助:
    国家自然科学基金(42174003);甘肃省科技计划资助(22JR5RA315;21JR7RA318)

Self-adaptive extraction method of tectonic movement change recorded by GNSS continuous observations

SU Xiaoning1, SHI Ruijuan1, BAO Qinghua1, ZHU Qing2, MENG Guojie3, YAN Haowen1   

  1. 1. Faculty of Geomatics, Lanzhou Jiaotong University, Lanzhou 730070, China;
    2. Faculty of Geosciences and Environmental Engineering, Southwest Jiaotong University, Chengdu 611756, China;
    3. Institute of Earthquake Forecasting, China Earthquake Administration, Beijing 100036, China
  • Received:2022-11-28 Revised:2023-03-28 Published:2023-09-07
  • Supported by:
    The National Natural Science Foundation of China (No. 42174003); Gansu Province Science and Technology Plan (Nos. 22JR5RA315;21JR7RA318)

摘要: 随着GNSS连续观测站数据积累时长和空间密度的不断增加,坐标时间序列呈现出随时间变化的构造运动变化信息,因此基于现有函数模型对坐标时间序列建模已不具备普遍适用性。本文提出一种GNSS坐标时间序列中构造运动变化特征的自适应提取方法,以获取构造运动随时间变化特征。首先,构建表征GNSS坐标时间序列中构造运动随时间变化的改进函数模型,引入基于贝叶斯框架的马尔可夫链·蒙特卡洛算法自适应求解模型参数最优解及其误差;然后,试验模拟不同噪声水平和不同分段线性速率差值的坐标时间序列,分析函数模型中线性速率和时间节点参数的精度,结果表明在当前GNSS单日解水平分量定位精度下,水平分量线性速率精度为0.1 mm/a,可分辨出0.4 mm/a及以上随时间变化的差异性构造运动,时间节点的精度与分段线性速率差异值有关,差异值为0.4 mm/a时其精度为1.0 a;最后,利用青藏高原东北隅GNSS连续站观测数据对本文方法有效性进行验证,并在青藏高原东北隅最前沿区域提取出始于2017年东北向运动减缓的趋势性构造运动变化特征。

关键词: GNSS时间序列, 传统函数模型, 改进函数模型, 构造运动变化, 自适应提取

Abstract: With the increase of spatial density and time span of GNSS observations, more and more abundant information including the time-dependent tectonic deformation has been presented in GNSS coordinate time series, which makes the traditional function model of GNSS coordinate timeseries fitting no longer have universal applicability. To overcome this problem, we propose an adaptive algorithm to extract the characteristics of tectonic deformation in GNSS coordinate timeseries. Firstly, a modified functional model of coordinate time series is constructed to characterize the variation characteristics of tectonic motion of GNSS stations and the adaptive algorithm of Bayesian framework is introduced to solve the time node of linear rate change as well as the posterior probability density function of parameters to obtain the optimal solution of parameters. Secondly, the coordinate time series of different noise levels and different piecewise linear rate differences are simulated. The accuracy of the linear rate and time nodes are analyzed. The results show that with the current GNSS daily positioning accuracy for the horizontal component, the calculated linear rate accuracy is 0.1 mm/a, the tectonic motion of more than 0.4 mm/a can be distinguished. The resolution of the tectonic motion change time node is related to the difference value of the tectonic motion, and its accuracy is 1.0 year when the difference value is 0.4 mm/a. is 1.0 years. Finally, the effectiveness of the method is verified by using the GNSS continuous observations on the northeastern corner of the Tibetan Plateau. At the forefront of the northeastern corner of the Tibetan Plateau, the regional tectonic kinematic variation characteristics that occurred in 2017 and the obstruction of northeastward movement are extracted.

Key words: GNSS coordinate timeseries, traditional functional model, modified functional model, tectonic movement change, self-adaptive extraction

中图分类号: