欢迎访问《地震》,

地震 ›› 2021, Vol. 41 ›› Issue (4): 57-67.doi: 10.12196/j.issn.1000-3274.2021.04.004

• • 上一篇    下一篇

基于GNSS数据的三种应变计算方法在青藏高原东南缘地区的应用

金阳1, 金红林1, 高原1, 黄星2   

  1. 1.中国地震局地震预测研究所, 地震预测重点实验室, 北京 100036;
    2.中国地震局地质研究所, 北京 100029
  • 收稿日期:2021-04-29 修回日期:2021-06-27 出版日期:2021-10-31 发布日期:2022-03-25
  • 通讯作者: 金红林,副研究员。E-mail:jhl_1968@126.com
  • 作者简介:金阳(1996-),男,江苏南通人,在读硕士研究生,主要从事地壳形变研究。
  • 基金资助:
    国家重点研发计划(2018YFC1503704); 自然科学基金重点项目 (41730212); 中国地震局地震预测研究所重点项目(2017IES0102)

Application of Three Strain Calculation Methods Based on GNSS Data in the Southeastern Margin of Tibet Plateau

JIN Yang1, JIN Hong-lin1, GAO Yuan1, HUANG Xing2   

  1. 1. Insititute of Earthquake Forcasting, Key Laboratory of Earthquake Prediction, Beijing 100036, China;
    2. Institute of Geology, China Earthquake Administration, Beijing 100029, China
  • Received:2021-04-29 Revised:2021-06-27 Online:2021-10-31 Published:2022-03-25

摘要: 利用GNSS速度场求解应变率的方法存在多元化特征, 可以区分为数学方法和物理方法。 本文选取求解应变率的三种数学方法: 最小二乘配置法、 多尺度球面小波法、 Gussian加权内插法, 对1999—2019年青藏高原东南缘(97°E~106°E, 21°N~30°N)GNSS速度场以及均匀GNSS速度场进行解算, 探讨以上三种方法在该区域的适用性特征。 通过分析由两组不同数据分布密度获得的主应变率(最大/最小)、 最大剪切应变率、 面应变率结果, 发现利用多尺度球面小波法获得的应变参数在青藏高原东南缘区域有更高的空间分辨率表现。 在数据分布密度高时可以详细反映小区域内空间变形特征, 有助于详细确定断裂变形特征。 本研究结果仅代表在青藏高原东南缘研究区域结果, 对于其他区域计算方法的选取仍需考虑其构造变形背景、 测点分布密度等条件。

关键词: GNSS数据, 应变率, 最小二乘配置法, 多尺度球面小波法, Gussian加权内插法

Abstract: The methods of using the GNSS velocity field to solve the strain rate have multiple characteristics, which can be divided into mathematical methods and physical methods. This paper selects three mathematical methods to solve the strain rate: Least Squares Collocation Method, Multi-scale Spherical Wavelet Method, Gussian Weighted Interpolation Method, to resolve the GNSS velocity field and uniform GNSS velocity field on the southeastern margin of Tibet Plateau from 1999 to 2019 (97°E~106°E, 21°N~30°N), and the applicability characteristics of the above three methods in this area is discussed. By analyzing the results of principal strain rate (maximum/minimum), maximum shear strain rate, and surface strain rate obtained from two sets of different data distributions, we found that the strain parameters obtained by multi-scale spherical wavelet method have higher resolution in the southeastern margin of Tibet Plateau. In the high-density data distribution, the spatial deformation characteristics in a small area can be reflected in detail, which is helpful to determine the fault deformation characteristics. The result only represents the studied area on the southeastern margin of Tibet Plateau, and the selection of calculation methods for other areas still needs to consider the conditions of the structural deformation background and the distribution density of measurement points.

Key words: GNSS data, Strain rate, Least squares collocation method, Multi-scale spherical wavelet method, Gussian weighted interpolation method

中图分类号: