测绘学报 ›› 2019, Vol. 48 ›› Issue (2): 191-197.doi: 10.11947/j.AGCS.2019.20180261

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

动态GNSS监测瞬态无震蠕滑信号提取

侯争1, 郭增长1,2, 杜久升1   

  1. 1. 河南理工大学, 河南 焦作 454003;
    2. 河南测绘职业学院, 河南 郑州 451464
  • 收稿日期:2018-06-05 修回日期:2018-11-05 出版日期:2019-02-20 发布日期:2019-03-02
  • 通讯作者: 郭增长 E-mail:guohpu@163.com
  • 作者简介:侯争(1982-),男,博士生,研究方向为大地测量数据处理和地球动力学。E-mail:houzheng.com@163.com
  • 基金资助:

    国家自然科学基金(41774041);教育部人文社会科学青年基金(17YJCZH041);煤炭工业协会科学技术研究指导性计划(MTKJ2016-212);河南省自然科学基金(182300410111)

Extraction of transient aseismic creep signal monitored by dynamic GNSS

HOU Zheng1, GUO Zengzhang1,2, DU Jiusheng1   

  1. 1. Henan Polytechnic University, Jiaozuo 454003, China;
    2. Henan College of Surveying and Mapping, Zhengzhou 451464, China
  • Received:2018-06-05 Revised:2018-11-05 Online:2019-02-20 Published:2019-03-02
  • Supported by:

    The National Natural Science Foundation of China (No. 41774041);The Humanities and Social Science Project of the Ministry of Education (No. 17YJCZH041);The Guiding Program for National Coal Association Science and Technology Research (No. MTKJ2016-212);The National Natural Science Foundation of Henan Province of China (No. 182300410111)

摘要:

断层的瞬态无震蠕滑常诱发震级高、破坏性强的地震。针对其滑动速度缓慢,难以探测的问题,本文基于GNSS连续坐标时间序列的异常波动提出一种断层瞬态无震蠕滑信息的自动探测方法。首先利用独立成分分析进行时空滤波,提高坐标时间序列的信噪比;然后计算坐标时间序列波动的相对强度指数及峰度值;最后通过累积分布函数将其转换为蠕滑信号概率,进而探测断层蠕滑事件。本文模拟500 d GNSS地表位移时间序列进行仿真试验,其中包含25 d瞬态蠕滑信号。试验结果表明,当信号强度至少与噪声水平相当时可有效探测出断层的蠕滑信息。计算Akutan地区连续3年的GNSS数据后探测到一个蠕滑信号,推断其可能为火山岩强烈运动引起的无震蠕滑。通过对四川省陆态网18个测站7年的GNSS数据处理后发现了4个异常信号,结果表明这些信号可能与震前应力积累造成的地表异常移动和震后断层余滑有密切关系。

关键词: 无震蠕滑, 独立分量分析, 相对强度指数, 峰度值, 累积分布函数, 慢滑移

Abstract:

Transient aseismic creep of the fault frequently induces earthquakes with high magnitude and destructiveness. In allusion to the challenges of slow sliding speed and difficult to detect, an automatic detection method of transient aseismic creep information for the fault was proposed, which based on the abnormal fluctuations of GNSS continuous coordinate time series. First, independent component analysis was used to improve the signal-to-noise ratio (SNR); then the relative strength index and kurtosis value of the fluctuations of the coordinate time series were calculated; finally a creep signal probability was converted through the cumulative distribution function, so as to the fault creep event was detected. In this paper, a 500-day GNSS surface displacement time series was simulated that included a 25-day transient creep signal. The experimental results show that the creep information of the fault can be effectively detected when the signal strength was at least equivalent to the noise level. After calculating the GNSS data for three consecutive years in the Akutan Zone, onecreep signalwas detected and it might be an aseismic creep signal related to the strong movement of the volcanic. In accordance with the processing results of seven-year GNSS data from 18 stations of the Crustal Movement Observation Network of China in Sichuan province, four abnormal signals was found. The analysis results indicated that these signals may be closely related to the abnormal displacement caused by stress accumulation and fault creep caused by the earthquake.

Key words: aseismic creep, independent component analysis, relative strength index, kurtosis value, cumulative distribution function, slow slip

中图分类号: