留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

北斗三号与超高频GNSS同震形变监测:以2021年青海玛多Mw 7.4地震为例

柴海山 陈克杰 魏国光 方荣新 邹蓉 祝会忠

柴海山, 陈克杰, 魏国光, 方荣新, 邹蓉, 祝会忠. 北斗三号与超高频GNSS同震形变监测:以2021年青海玛多Mw 7.4地震为例[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 946-954. doi: 10.13203/j.whugis20220140
引用本文: 柴海山, 陈克杰, 魏国光, 方荣新, 邹蓉, 祝会忠. 北斗三号与超高频GNSS同震形变监测:以2021年青海玛多Mw 7.4地震为例[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 946-954. doi: 10.13203/j.whugis20220140
CHAI Haishan, CHEN Kejie, WEI Guoguang, FANG Rongxin, ZOU Rong, ZHU Huizhong. Coseismic Deformation Monitoring Using BDS-3 and Ultra-High Rate GNSS: A Case Study of the 2021 Maduo Mw 7.4 Earthquake[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 946-954. doi: 10.13203/j.whugis20220140
Citation: CHAI Haishan, CHEN Kejie, WEI Guoguang, FANG Rongxin, ZOU Rong, ZHU Huizhong. Coseismic Deformation Monitoring Using BDS-3 and Ultra-High Rate GNSS: A Case Study of the 2021 Maduo Mw 7.4 Earthquake[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 946-954. doi: 10.13203/j.whugis20220140

北斗三号与超高频GNSS同震形变监测:以2021年青海玛多Mw 7.4地震为例

doi: 10.13203/j.whugis20220140
基金项目: 

国家自然科学基金 42074024

详细信息

Coseismic Deformation Monitoring Using BDS-3 and Ultra-High Rate GNSS: A Case Study of the 2021 Maduo Mw 7.4 Earthquake

Funds: 

The National Natural Science Foundation of China 42074024

More Information
  • 摘要: 随着北斗三号完成全球组网以及硬件技术发展,地壳形变监测网络的全球导航卫星系统(global navigation satellite system, GNSS)接收机正朝着多系统、超高频(≥50 Hz)更新升级。基于陆态网络和中国青海省连续运行参考站的GNSS观测数据,获取2021年青海玛多Mw 7.4地震的同震位移波,并进行了震前精度统计和功率谱密度分析。结果表明,北斗三号在本次地震中的位移拾取精度比全球定位系统(global positioning system,GPS)更优;采样率提高到50 Hz对削弱位移时间序列混叠效应作用不显著,在10 Hz左右达到饱和;对50 Hz同震位移波差分获取到的速度和加速度序列噪声较大,难以利用。综合应用场景和成本考虑,GNSS地震监测的采样率最高设置在5~10 Hz即可。
  • 图  1  2021年玛多Mw 7.4地震的测站分布及同震位移

    Figure  1.  Distribution of GNSS Stations and Coseismic Displacement of the 2021 Maduo Mw 7.4 Earthquake

    图  2  SCGZ测站时间序列和各测站平均功率谱密度

    Figure  2.  Time Series at SCGZ Station and Average Power Spectrum Density of All Stations

    图  3  QSHZ测站的时间序列和功率谱密度

    Figure  3.  Time Series and Power Spectrum Density at QSHZ Station

    图  4  QHMQ站与DAW63站的同震位移、速度和加速度

    Figure  4.  Coseismic Deformation, Velocity and Acceleration at QHMQ and DAW63 Stations

    图  5  各采样率的时间序列功率谱密度

    Figure  5.  Power Spectrum Density of Time Series with Different Sampling Rates

    表  1  GNSS数据信息

    Table  1.   Information of GNSS Data

    数据来源 测站个数 GNSS采样率/Hz 接收机型号 卫星系统
    陆态网络参考站 7 1,50 Trimble NetR9 GRC2,GREC2
    青海省QSHZ站 1 1 中海达VNet8 GREC2+3
    下载: 导出CSV

    表  2  PPP处理策略和参数设置

    Table  2.   Processing Strategy and Parameters Setting of PPP

    策略与参数 具体描述
    组合方式 无电离层组合
    截止高度角/(°) 7
    观测值定权 高度角加权,伪距、载波先验误差:0.6 m、0.01周
    系统定权 G∶R∶E∶C2∶C3=1.5∶1∶1.5∶1∶1.5
    估计方式 最小二乘
    天线改正 IGS14.atx
    周跳探测 Turboedit
    地球自转参数 IGS EOP
    对流层估计 Saastamoinen模型+分段线性估计
    精密星历/钟差 CODE:30 s钟差,15 min轨道(GREC2+3)
    下载: 导出CSV

    表  3  不同系统时间序列的标准差均值/mm

    Table  3.   Average Standard Deviation of Time Series of Diffferent Systems /mm

    系统 E方向 N方向 U方向
    G 9.4 12.4 39.1
    GRC2 8.1 7.1 17.0
    GREC2 7.8 6.5 16.8
    下载: 导出CSV

    表  4  QSHZ测站不同系统时间序列标准差/mm

    Table  4.   Standard Deviation of Time Series of Diffferent Systems at QSHZ Station /mm

    系统 E方向 N方向 U方向
    G 9.2 9.5 43.6
    E 7.5 10.4 19.6
    C2+3 5.3 5.4 18.1
    GREC2 4.8 5.4 16.2
    GREC2+3 3.4 4.7 12.4
    下载: 导出CSV
  • [1] Bock Y, Prawirodirdjo L, Melbourne T I. Detection of Arbitrarily Large Dynamic Ground Motions with a Dense High-Rate GPS Network[J]. Geophysical Research Letters, 2004, 31(6): 1-4
    [2] Trifunac M D, Todorovska M I. Evolution of Accelerographs, Data Processing, Strong Motion Arrays and Amplitude and Spatial Resolution in Recording Strong Earthquake Motion[J]. Soil Dynamics and Earthquake Engineering, 2001, 21 (6) : 537-555 doi:  10.1016/S0267-7261(01)00013-6
    [3] Bock Y, Melgar D, Crowell B W. Real-Time Strong-Motion Broadband Displacements from Collocated GPS and Accelerometers[J]. Bulletin of the Seismological Society of America, 2011, 101(6): 2904-2925 doi:  10.1785/0120110007
    [4] Chen K J, Ge M R, Babeyko A, et al. Retrieving Real-Time Co-Seismic Displacements Using GPS/ GLONASS: A Preliminary Report from the Septem ber 2015 Mw 8. 3 Illapel Earthquake in Chile[J]. Geophysical Journal International, 2016, 206(2): 941-953 doi:  10.1093/gji/ggw190
    [5] Li X X, Ge M R, Zhang X H, et al. Real-Time High-Rate Co-Seismic Displacement from Ambiguity-Fixed Precise Point Positioning: Application to Earthquake Early Warning[J]. Geophysical Research Letters, 2013, 40(2): 295-300 doi:  10.1002/grl.50138
    [6] Avallone A, Marzario M, Cirella A, et al. Very High Rate(10 Hz)GPS Seismology for ModerateMagnitude Earthquakes: The Case of the Mw 6.3 L'Aquila(Central Italy)Event[J]. Journal of Geophysical Research, 2011, 116(B2): B02305
    [7] Geng T, Xie X, Fang R X, et al. Real-Time Capture of Seismic Waves Using High-Rate Multi-GNSS Observations: Application to the 2015 Mw 7. 8 Nepal Earthquake[J]. Geophysical Research Letters, 2016, 43(1): 161-167 doi:  10.1002/2015GL067044
    [8] Larson K M. GPS Seismology[J]. Journal of Geodesy, 2009, 83(3/4): 227-233
    [9] Yang Y X, Mao Y, Sun B J. Basic Performance and Future Developments of BeiDou Global Navigation Satellite System[J]. Satellite Navigation, 2020, 1: 1 doi:  10.1186/s43020-019-0006-0
    [10] Wu W W, Guo F, Zheng J Z. Analysis of Galileo Signal-in-Space Range Error and Positioning Performance During 2015-2018[J]. Satellite Navigation, 2020, 1: 6 doi:  10.1186/s43020-019-0005-1
    [11] Geng J H, Jiang P, Liu J N. Integrating GPS with GLONASS for High-Rate Seismogeodesy[J]. Geophysical Research Letters, 2017, 44(7): 3139-3146 doi:  10.1002/2017GL072808
    [12] Li X X, Zheng K, Li X, et al. Real-Time Capturing of Seismic Waveforms Using High-Rate BDS, GPS and GLONASS Observations: The 2017 Mw 6. 5 Jiuzhaigou Earthquake in China[J]. GPS Solutions, 2018, 23(1): 1-12
    [13] Xu P L, Shu Y M, Niu X J, et al. High-Rate Multi-GNSS Attitude Determination: Experiments, Comparisons with Inertial Measurement Units and Applications of GNSS Rotational Seismology to the 2011 Tohoku Mw 9. 0 Earthquake[J]. Measurement Science and Technology, 2019, 30(2): 024003
    [14] Geng J, Pan Y, Li X, et al. Noise Characteristics of High-Rate Multi-GNSS for Subdaily Crustal Deformation Monitoring[J]. Journal of Geophysical Research: Solid Earth, 2018, 123(2): 1987-2002 doi:  10.1002/2018JB015527
    [15] Fang R X, Lv H H, Shu Y M, et al. Improved Performance of GNSS Precise Point Positioning for High-Rate Seismogeodesy with Recent BDS-3 and Galileo[J]. Advances in Space Research, 2021, 68 (8): 3255-3267 doi:  10.1016/j.asr.2021.06.012
    [16] Genrich J F, Bock Y. Instantaneous Geodetic Positioning with 10-50 Hz GPS Measurements: Noise Characteristics and Implications for Monitoring Networks[J]. Journal of Geophysical Research: Solid Earth, 2006, 111(B3): 1-12
    [17] Yi T H, Li H N, Gu M. Experimental Assessment of High-Rate GPS Receivers for Deformation Monitoring of Bridge[J]. Measurement, 2013, 46(1): 420-432 doi:  10.1016/j.measurement.2012.07.018
    [18] Xu P L, Shi C, Fang R X, et al. High-Rate Precise Point Positioning(PPP)to Measure Seismic Wave Motions: An Experimental Comparison of GPS PPP with Inertial Measurement Units[J]. Journal of Geodesy, 2013, 87(4): 361-372 doi:  10.1007/s00190-012-0606-z
    [19] 方荣新. 高采样率GPS数据非差精密处理方法及其在地震学中的应用研究[D]. 武汉: 武汉大学, 2010

    Fang Rongxin. High-Rate GPS Data Non-Difference Precise Processing and Its Application on Seismology [D]. Wuhan: Wuhan University, 2010
    [20] Smalley R. High-Rate GPS: How High do We Need to Go?[J]. Seismological Research Letters, 2009, 80(6): 1054-1061 doi:  10.1785/gssrl.80.6.1054
    [21] 刘刚, 聂兆生, 贾治革, 等. 高频GNSS地震监测中频率混叠现象的仿真实验研究[J]. 地震研究, 2017, 40 (2): 186-192 doi:  10.3969/j.issn.1000-0666.2017.02.003

    Liu Gang, Nie Zhaosheng, Jia Zhige, et al. Simulation Experiment on the Aliasing of High Rate GNSS Used for Earthquake Monitoring[J]. Journal of Seismological Research, 2017, 40(2): 186-192 doi:  10.3969/j.issn.1000-0666.2017.02.003
    [22] Li M, Huang D F, Yan L, et al. Quasi-Real Time Determination of 2013 Lushan Mw 6. 6 Earthquake Epicenter, Trigger Time and Magnitude Using 50 Hz GPS Observations[C]// China Satellite Navigation Conference(CSNC), Nanjing, China, 2014
    [23] Xu P L, Du F, Shu Y M, et al. Regularized Reconstruction of Peak Ground Velocity and Acceleration from very High-Rate GNSS Precise Point Positioning with Applications to the 2013 Lushan Mw 6. 6 Earthquake[J]. Journal of Geodesy, 2021, 95(1): 1-22 doi:  10.1007/s00190-020-01448-7
    [24] 李志才, 丁开华, 张鹏, 等. GNSS观测的2021年青海玛多地震(Mw 7.4)同震形变及其滑动分布[J]. 武汉大学学报·信息科学版, 2021, 46(10): 1489-1497 doi:  10.13203/j.whugis20210301

    Li Zhicai, Ding Kaihua, Zhang Peng, et al. Coseismic Deformation and Slip Distribution of 2021 Mw 7.4 MadoiEarthquake from GNSS Observation[J]. Geomatics and Information Science of Wuhan University, 2021, 46(10): 1489-1497 doi:  10.13203/j.whugis20210301
    [25] 姜卫平, 许才军, 李志伟, 等. 利用空间观测技术研究青海玛多7.4级地震孕育发生变形时空特征[J]. 地球物理学报, 2022, 65(2): 495-508 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX202202006.htm

    Jiang Weiping, Xu Caijun, Li Zhiwei, et al. Using Space Observation Techniques to Study Temporal and Spatial Characteristics of Seismogenic Process, Occurrence and Deformation of the Qinghai MadoiMw 7.4 Earthquake[J]. Chinese Journal of Geophysics, 2022, 65(2): 495-508 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX202202006.htm
    [26] 王迪晋, 王东振, 赵斌, 等. 2021年青海玛多Mw 7.4地震GNSS同震形变场及其断层滑动分布[J]. 地球物理学报, 2022, 65(2): 537-551 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX202202018.htm

    Wang Dijin, Wang Dongzhen, Zhao Bin, et al. 2021 Qinghai Madoi Mw 7.4 Earthquake Coseismic Deformation Field and Fault-Slip Distribution Using GNSS Observations[J]. Chinese Journal of Geophysics, 2022, 65(2): 537-551 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX202202018.htm
    [27] 余鹏飞, 熊维, 陈威, 等. 基于GNSS和InSAR约束的2021年玛多Ms 7.4地震同震滑动分布及应用[J]. 地球物理学报, 2022, 65(2): 509-522 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX202202007.htm

    Yu Pengfei, Xiong Wei, Chen Wei, et al. Slip Model of the 2021 Ms 7.4 MadoiEarthquake Constrained by GNSS and InSAR Coseismic Deformation[J]. Chinese Journal of Geophysics, 2022, 65(2): 509-522 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX202202007.htm
    [28] Lou Y D, Zhang W X, Shi C, et al. High-Rate(1-Hz and 50-Hz)GPS Seismology: Application to the 2013 Mw 6. 6 Lushan Earthquake[J]. Journal of Asian Earth Sciences, 2014, 79: 426-431 doi:  10.1016/j.jseaes.2013.10.016
    [29] Shi C, Zhao Q L, Geng J H, et al. Recent Development of PANDA Software in GNSS Data Processing [C]// International Conference on Earth Observation Data Processing and Analysis (ICEODPA), Wuhan, China, 2008
    [30] 耿涛, 赵齐乐, 刘经南, 等. 基于PANDA软件的实时精密单点定位研究[J]. 武汉大学学报·信息科学版, 2007, 32(4): 312-315 http://ch.whu.edu.cn/article/id/1877

    Geng Tao, Zhao Qile, Liu Jingnan, et al. RealTime Precise Point Positioning Based on PANDA Software[J]. Geomatics and Information Science of Wuhan University, 2007, 32(4): 312-315 http://ch.whu.edu.cn/article/id/1877
    [31] 方荣新, 施闯, 宋伟伟, 等. 实时GNSS地震仪系统实现及精度分析[J]. 地球物理学报, 2013, 56(2): 450-458 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201302008.htm

    Fang Rongxin, Shi Chuang, Song Weiwei, et al. Real-Time GNSS Seismometer and Its Accuracy [J]. Chinese Journal of Geophysics, 2013, 56(2): 450-458 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201302008.htm
    [32] 方荣新, 施闯, 陈克杰, 等. GPS地震仪: PANDA软件测试结果与验证[J]. 武汉大学学报· 信息科学版, 2011, 36(4): 453-456 http://ch.whu.edu.cn/article/id/522

    Fang Rongxin, Shi Chuang, Chen Kejie, et al. GPS Seismometer: PANDA Software Testing Results and Validation[J]. Geomatics and Information Science of Wuhan University, 2011, 36(4): 453-456 http://ch.whu.edu.cn/article/id/522
    [33] 张小红, 郭斐, 李星星. IGS卫星钟差产品采样间隔对PPP精度的影响[J]. 武汉大学学报·信息科学版, 2010, 35(2): 152-155 http://ch.whu.edu.cn/article/id/860

    Zhang Xiaohong, Guo Fei, Li Xingxing. Impact of Sample Rate of IGS Satellite Clock on Precise Point Positioning[J]. Geomatics and Information Science of Wuhan University, 2010, 35(2): 152-155 http://ch.whu.edu.cn/article/id/860
    [34] Wang R, Schurr B, Milkereit C, et al. An Improved Automatic Scheme for Empirical Baseline Correction of Digital Strong-Motion Records[J]. Bulletin of the Seismological Society of America, 2011, 101(5): 2029-2044 doi:  10.1785/0120110039
  • [1] 李艳艳, 殷海涛, 韩林桥.  利用MSMPCA去噪的CEEMD方法监测高频GNSS同震形变 . 武汉大学学报 ● 信息科学版, 2022, 47(3): 352-360. doi: 10.13203/j.whugis20190356
    [2] 曾添, 隋立芬, 贾小林, 肖国锐, 戴卿, 甘雨.  小型化LEO星座与BDS-3全星座联合定轨仿真 . 武汉大学学报 ● 信息科学版, 2022, 47(1): 61-68. doi: 10.13203/j.whugis20190426
    [3] 余鹏飞, 陈威, 乔学军, 赵斌, 李刚, 熊维.  基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究 . 武汉大学学报 ● 信息科学版, 2022, 47(6): 898-906. doi: 10.13203/j.whugis20220093
    [4] 江颖, 刘子维, 张晓彤, 张丽娜, 韦进.  2021年青海玛多Mw 7.4地震前后b值的变化特征研究 . 武汉大学学报 ● 信息科学版, 2022, 47(6): 907-915. doi: 10.13203/j.whugis20220071
    [5] 李志才, 丁开华, 张鹏, 温扬茂, 赵利江, 陈建峰.  GNSS观测的2021年青海玛多地震(Mw 7.4)同震形变及其滑动分布 . 武汉大学学报 ● 信息科学版, 2021, 46(10): 1489-1497. doi: 10.13203/j.whugis20210301
    [6] 蔡洪亮, 孟轶男, 耿涛, 谢新.  北斗三号卫星星地星间联合精密定轨初步结果 . 武汉大学学报 ● 信息科学版, 2020, 45(10): 1493-1500. doi: 10.13203/j.whugis20180499
    [7] 李水平, 陈刚, 何平, 丁开华, 陈云锅, 王琪.  2015年尼泊尔地震同震滑动及震后余滑的三角位错模型反演 . 武汉大学学报 ● 信息科学版, 2019, 44(12): 1787-1796. doi: 10.13203/j.whugis20180128
    [8] 周仁宇, 胡志刚, 苏牡丹, 李军正, 李鹏博, 赵齐乐.  北斗全球系统广播电离层模型性能初步评估 . 武汉大学学报 ● 信息科学版, 2019, 44(10): 1457-1464. doi: 10.13203/j.whugis20170428
    [9] 苏利娜, 张勇.  GPS时间序列中同震和震后形变的自动识别和估计 . 武汉大学学报 ● 信息科学版, 2018, 43(10): 1504-1510. doi: 10.13203/j.whugis20170016
    [10] 许才军, 邓长勇, 周力璇.  利用方差分量估计的地震同震滑动分布反演 . 武汉大学学报 ● 信息科学版, 2016, 41(1): 37-44. doi: 10.13203/j.whugis20150500
    [11] 杨国华, 朱爽, 梁洪宝, 杨博.  芦山Ms 7.0级地震震前及同震地表形变 . 武汉大学学报 ● 信息科学版, 2015, 40(1): 121-127.
    [12] 沈飞, 李建成, 郭斐.  单历元PPP分析日本Mw9.0地震引起的我国东部沿海同震地表形变 . 武汉大学学报 ● 信息科学版, 2012, 37(11): 1345-1347.
    [13] 王永哲, 李志伟, 朱建军, 胡俊.  融合多平台DInSAR数据解算拉奎拉地震三维同震形变场 . 武汉大学学报 ● 信息科学版, 2012, 37(7): 859-863.
    [14] 许才军, 何平, 温扬茂, 张磊.  日本2011 Tohoku-Oki Mw 9.0级地震的同震形变及其滑动分布反演:GPS和InSAR约束 . 武汉大学学报 ● 信息科学版, 2012, 37(12): 1387-1391.
    [15] 刁法启, 熊熊, 李军, 郑勇.  利用GPS观测反演Tokachi-Oki地震(日本)及强余震的同震滑移分布 . 武汉大学学报 ● 信息科学版, 2010, 35(3): 270-273.
    [16] 邢乐林, 李建成, 李辉, 孙文科.  Sumatra-Adaman大地震同震和震后形变的GRACE卫星检测 . 武汉大学学报 ● 信息科学版, 2009, 34(9): 1080-1084.
    [17] 彭颖, 许才军, 刘洋.  联合地震位错模型和InSAR数据构建2017年九寨沟Mw6.5地震同震三维形变场 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20200289
    [18] 阚昊宇, 胡志刚, 吕逸飞, 谢新, 周仁宇, 赵齐乐.  利用不同时间同步体制钟差评估北斗三号星载原子钟性能 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20210286
    [19] 宋伟伟, 赵新科, 楼益栋, 孙唯彬, 赵正予.  北斗三号PPP-B2b服务性能评估 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20200686
    [20] 王西龙, 许小龙, 赵齐乐.  北斗三号系统信号质量分析及轨道精度验证 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20200577
  • 加载中
图(5) / 表(4)
计量
  • 文章访问数:  187
  • HTML全文浏览量:  58
  • PDF下载量:  64
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-17
  • 刊出日期:  2022-06-05

北斗三号与超高频GNSS同震形变监测:以2021年青海玛多Mw 7.4地震为例

doi: 10.13203/j.whugis20220140
    基金项目:

    国家自然科学基金 42074024

    作者简介:

    柴海山,硕士生,主要研究方向为GNSS数据处理及其应用。12032535@mail.sustech.edu.cn

    通讯作者: 陈克杰,博士,助理教授。chenkj@sustech.edu.cn
  • 中图分类号: P237

摘要: 随着北斗三号完成全球组网以及硬件技术发展,地壳形变监测网络的全球导航卫星系统(global navigation satellite system, GNSS)接收机正朝着多系统、超高频(≥50 Hz)更新升级。基于陆态网络和中国青海省连续运行参考站的GNSS观测数据,获取2021年青海玛多Mw 7.4地震的同震位移波,并进行了震前精度统计和功率谱密度分析。结果表明,北斗三号在本次地震中的位移拾取精度比全球定位系统(global positioning system,GPS)更优;采样率提高到50 Hz对削弱位移时间序列混叠效应作用不显著,在10 Hz左右达到饱和;对50 Hz同震位移波差分获取到的速度和加速度序列噪声较大,难以利用。综合应用场景和成本考虑,GNSS地震监测的采样率最高设置在5~10 Hz即可。

English Abstract

柴海山, 陈克杰, 魏国光, 方荣新, 邹蓉, 祝会忠. 北斗三号与超高频GNSS同震形变监测:以2021年青海玛多Mw 7.4地震为例[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 946-954. doi: 10.13203/j.whugis20220140
引用本文: 柴海山, 陈克杰, 魏国光, 方荣新, 邹蓉, 祝会忠. 北斗三号与超高频GNSS同震形变监测:以2021年青海玛多Mw 7.4地震为例[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 946-954. doi: 10.13203/j.whugis20220140
CHAI Haishan, CHEN Kejie, WEI Guoguang, FANG Rongxin, ZOU Rong, ZHU Huizhong. Coseismic Deformation Monitoring Using BDS-3 and Ultra-High Rate GNSS: A Case Study of the 2021 Maduo Mw 7.4 Earthquake[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 946-954. doi: 10.13203/j.whugis20220140
Citation: CHAI Haishan, CHEN Kejie, WEI Guoguang, FANG Rongxin, ZOU Rong, ZHU Huizhong. Coseismic Deformation Monitoring Using BDS-3 and Ultra-High Rate GNSS: A Case Study of the 2021 Maduo Mw 7.4 Earthquake[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 946-954. doi: 10.13203/j.whugis20220140
  • 实时高频(≥1 Hz)的全球导航卫星系统(global navigation satellite system,GNSS)可以克服地震仪饱和效应和强震仪基线偏差[1-2],在同震形变监测领域得到广泛应用[3-7],因此GNSS地震学这一交叉学科应运而生[8]。由于全球定位系统(global positioning system,GPS)建成时间最早、使用最广泛、表现最成熟,GNSS地震学在一段时间内多特指GPS。随着GLONASS、北斗三号、Galileo的建设发展[9-10],GNSS多系统组合在同震形变监测中精度更高、稳定性更好、震源参数反演更可靠[11-14]。2020年6月,北斗三号完成全球组网,文献[15]使用振动平台进行了一系列模拟实验,证实了北斗三号和Galileo对地震位移拾取的积极作用,但北斗三号在天然地震中的应用报道尚不多见。

    随着硬件和储存技术的发展,当前GNSS接收机采样率已越来越接近地震仪和强震仪,可达到50~100 Hz[16-18](本文称之为超高频),中国大陆构造环境监测网络(简称陆态网络)的GNSS测站也已普遍升级,可提供50 Hz的观测记录[19]。但是目前仍有许多国家的GNSS同震形变观测网以1 Hz、5 Hz为主流,如美国的板块边界观测网络(plate boundary observatory,PBO),仅部分站可达10 Hz[20]。虽然提高采样率有利于削弱位移时间序列的混叠效应[20-21],也存在获得超高采样率的速度、加速度序列的可能性[22-23],使GNSS接收机实现类似强震仪、地震仪的作用,但同时也极大地增加了数据储存、管理、计算的成本。

    北京时间2021-05-22T02:04:11青海玛多发生了Mw 7.4地震(https://news.ceic.ac.cn/),其同震过程被附近的陆态网络和当地的连续运行参考站(continuously operating reference stations,CORS)记录[24-27]。本文以此次地震为例,研究多系统组合同震位移拾取的精度,分析北斗三号的地表位移拾取性能,并且针对GNSS地震监测采样率设置存在的一些争议[82028],测试超高频数据获取同震位移、速度、加速度的性能,分析是否有必要将采样率提升至50 Hz。

    • 本文选取了震中附近7个陆态网络参考站(平均距离约280 km),获取2021年玛多Mw 7.4地震的1 Hz同震位移进行分析,如图 1所示。其中G表示GPS,R表示GLONASS,E表示Galileo,C2表示北斗二号,C3表示北斗三号,C2+3表示北斗二号和北斗三号组合,其他字母组合含义类似。图 1(b)中的位移波从上到下按照测站到震中的距离排序,均可见明显同震信号。现阶段陆态网络中绝大部分参考站支持GRC2三系统,部分参考站支持GREC2四系统。2021年玛多地震发生时,陆态网络附近参考站尚不支持北斗三号,为了开展相关研究,本文选取了中国青海省CORS站中的QSHZ站观测数据,该站可观测包括北斗三号在内的四大卫星导航系统,其观测数据详细的信息见表 1。本文所涉及的GNSS数据均采用PANDA软件[29-31]进行计算,其精密单点定位(precise point positioning,PPP)功能可以处理高达50 Hz的数据[32],处理策略和参数设置见表 2

      图  1  2021年玛多Mw 7.4地震的测站分布及同震位移

      Figure 1.  Distribution of GNSS Stations and Coseismic Displacement of the 2021 Maduo Mw 7.4 Earthquake

      表 1  GNSS数据信息

      Table 1.  Information of GNSS Data

      数据来源 测站个数 GNSS采样率/Hz 接收机型号 卫星系统
      陆态网络参考站 7 1,50 Trimble NetR9 GRC2,GREC2
      青海省QSHZ站 1 1 中海达VNet8 GREC2+3

      表 2  PPP处理策略和参数设置

      Table 2.  Processing Strategy and Parameters Setting of PPP

      策略与参数 具体描述
      组合方式 无电离层组合
      截止高度角/(°) 7
      观测值定权 高度角加权,伪距、载波先验误差:0.6 m、0.01周
      系统定权 G∶R∶E∶C2∶C3=1.5∶1∶1.5∶1∶1.5
      估计方式 最小二乘
      天线改正 IGS14.atx
      周跳探测 Turboedit
      地球自转参数 IGS EOP
      对流层估计 Saastamoinen模型+分段线性估计
      精密星历/钟差 CODE:30 s钟差,15 min轨道(GREC2+3)

      在处理陆态网络和青海省QSHZ站1 Hz数据时,考虑到卫星钟差在30 s范围内高频抖动并不剧烈,对于1 s采样率的GNSS数据,使用5 s间隔钟差与30 s间隔钟差获得的精度基本一致[33],故选择欧洲定轨中心(Center for Orbit Determination in Europe,CODE)提供的30 s钟差和15 min精密星历产品。

      通过PPP获取位移时间序列后,对其进行一阶离散差分便可获得速度序列,再对速度差分可获得加速度序列,计算式为:

      vi=si+1-si-12t ]]>
      ai=vi+1-vi-12t ]]>

      式中,sva分别表示位移、速度、加速度;i表示历元;t表示时间间隔。

    • 为了研究超高频GNSS差分获取速度和加速度的能力,本文选择QHMQ站及其附近2 km内的一个强震仪测站DAW63进行比较,DAW63与QHMQ具有高度相关性,可提供采样率为200 Hz的加速度序列。通过对原始加速度序列积分可获得速度序列,但由于基线偏差的存在,此时的速度序列包含线性漂移,以4阶巴特沃斯高通滤波的方式去除速度序列中频率在0.02 Hz以下的低频信号,再次积分即可得到地表形变。以上数据处理过程与经验基线自动校正方法[34]有一定区别,会对静态同震位移造成一定的削弱,但考虑到本文重点关注高频同震位移波且QHMQ站静态同震位移较小,上述处理方法是合理的。

    • 为了分析陆态网络测站1 Hz数据同震位移拾取精度,计算2021-05-21T13:00—18:00共计5 h震前时间序列与功率谱密度,出于篇幅考虑,图 2仅展示SCGZ单站的时间序列和各测站的平均功率谱密度。

      图  2  SCGZ测站时间序列和各测站平均功率谱密度

      Figure 2.  Time Series at SCGZ Station and Average Power Spectrum Density of All Stations

      图 2(a)可以直观看出,SCGZ单站的多系统组合(GRC2)相较于单系统G能够在一定程度上削弱低频噪声和跳变,使地表形变拾取更加稳健可靠。由图 2(b)可以看出,各测站的U方向功率谱密度较大,在任意频段均要高于E、N方向。随着系统数目的增加,各频段噪声普遍明显减弱,GREC2相比GRC2效果更好。单系统G在E方向噪声与GRC2差异不大,甚至部分频段噪声更低,但GREC2仍然表现最好。高频段出现个别峰值的原因在于位移时间序列存在一定的跳变,而这些跳变可能是由精密钟差的高频钟跳所导致。

      计算不同系统时间序列的标准差并作站间平均,结果见表 3。单系统G在E、N方向的标准差均值约为1 cm,U方向约为4 cm,主要是因为受卫星几何构型和多路径效应等影响,U方向误差较大。GRC2三系统组合相比于单系统G,在E、N、U方向上的精度分别提升了14%、42%、56%,GREC2四系统组合相比于单系统G,在E、N、U方向上的精度分别提升了17%、47%、57%,印证了多系统组合的优越性。

      表 3  不同系统时间序列的标准差均值/mm

      Table 3.  Average Standard Deviation of Time Series of Diffferent Systems /mm

      系统 E方向 N方向 U方向
      G 9.4 12.4 39.1
      GRC2 8.1 7.1 17.0
      GREC2 7.8 6.5 16.8
    • 北斗二号虽然在亚太地区卫星数与GPS相当,但其包含若干颗地球同步轨道(geosynchronous earth orbit,GEO)卫星,几何构型较差,对定位精度提升有限,单独采用北斗二号并不能取得良好的同震位移拾取效果。随着北斗三号的全面建成,这一情况得到了极大改善。

      对青海省QSHZ站进行PPP解算,截取与§2.1相同的时段进行时间序列分析,并计算了不同系统组合时间序列的标准差,结果见表 4。QSHZ站震前时间序列的功率谱密度如图 3所示。

      表 4  QSHZ测站不同系统时间序列标准差/mm

      Table 4.  Standard Deviation of Time Series of Diffferent Systems at QSHZ Station /mm

      系统 E方向 N方向 U方向
      G 9.2 9.5 43.6
      E 7.5 10.4 19.6
      C2+3 5.3 5.4 18.1
      GREC2 4.8 5.4 16.2
      GREC2+3 3.4 4.7 12.4

      图  3  QSHZ测站的时间序列和功率谱密度

      Figure 3.  Time Series and Power Spectrum Density at QSHZ Station

      图 3(a)可以看出,附加北斗三号的四系统组合精度最高,相较于单GPS定位结果,其在E、N、U方向上的精度均有显著提高。单GPS结果受到几何构型、多路径效应等影响,难以单独提供稳定的时间序列,尤其是在N方向末段和U方向上结果较差,此部分误差和SCGZ测站较为相似,说明误差的起因更多在于卫星几何构型较差。就单系统而言,随着北斗三号的建成,北斗系统不仅可以单独实现同震位移监测,而且在单系统中效果最佳,其原因在于地震发生在亚太地区,北斗可观测卫星数更多,几何结构更好,且北斗三号搭载的原子钟更稳定,随着北斗三号系统的加入,同震形变的噪声能进一步得到削弱,提高了拾取精度和可靠性;从此时段结果来看,Galileo高频噪声较大,但低频噪声要低于GPS。由此可知,单系统受制于卫星几何构型等因素,其时间序列不如多系统稳健,采用多系统组合的方式对同震位移拾取有重要意义。

      图 3(b)可以看出,QSHZ测站震前时间序列的功率谱密度的基本特征与陆态网络平均功率谱密度一致。单系统之间相比,北斗在各频段的噪声均为最低;Galileo次之,且在0.1 Hz频段及以上噪声略大;GPS低频噪声水平最大,不够稳定,对应了图 3(a)中N方向末段和U方向的结果,四系统组合的噪声水平相比于单系统GPS在任意频段均有大幅削弱。以上结果表明,北斗不仅能单独实现高精度定位,且相比于GPS和Galileo噪声更小;在加入北斗三号后,GREC2+3的噪声相对于GREC2(代表现阶段陆态网络的普遍配置)在各个频段均有明显降低,时间序列稳定性进一步上升,说明北斗三号可以对同震位移拾取做出重要贡献。

    • 以QHMQ和DAW63为例,设置一系列采样率(1 Hz、2 Hz、5 Hz、10 Hz、25 Hz、50 Hz),获得相应的同震位移、速度和加速度序列。考虑到2021年青海玛多Mw 7.4地震为走滑型地震,以E、N方向的位移为主导,且GNSS在U方向的精度较差,本文重点研究E方向和N方向,GNSS与强震仪结果比较如图 4所示。为保证超高频处理结果的可靠性,将钟差间隔提升至5 s,选用武汉大学提供的5 s钟差以及15 min轨道计算,GRC2三系统参与运算,其余策略和表 2一致。

      图  4  QHMQ站与DAW63站的同震位移、速度和加速度

      Figure 4.  Coseismic Deformation, Velocity and Acceleration at QHMQ and DAW63 Stations

      图 4进行分析可得如下结论:

      1)就位移拾取而言,GNSS和强震仪在E、N方向上的位移拾取保持高度一致性;混叠效应随采样率提升而削弱,在5~10 Hz处达到饱和;N方向相比于E方向振幅更大,混叠效应整体较弱,50 Hz和1 Hz相比结果无明显改善。

      2)就速度拾取而言,差分噪声随着采样率的提升而增强,但同时混叠效应减弱;GNSS在N方向的2~5 Hz与强震仪保持了较好的一致性,但在E方向效果一般,原因在于N方向比E方向振幅更大,故速度序列噪声整体要更小;速度差分效果随振幅增大而改善,这表明GNSS更适合于近场以及大型地震监测。

      3)就加速度拾取而言,通过二次差分获得的结果始终未取得良好效果;在10 Hz时达到了强震仪的量级,但符合程度较差;随着采样率提升,噪声水平持续增强,有效信号被掩盖。

      对于其他震级的地震是否能得出类似结论有待进一步研究,值得注意的是,受限于定位精度,GNSS对中小型地震同震位移的拾取能力较弱,需要缩短测站与震中距离以实现更好的观测效果。

      为了进一步研究超高频数据对时间序列存在的混叠效应的削弱,选择3个与震中距离不同且具有代表性的测站QHMD(38 km)、QHMQ(175 km)和QHGE(366 km)进行功率谱密度分析,分析时段为2021-05-21T17:00—20:00(包含同震形变),结果如图 5所示。由图 5可以看出,随着采样率提升,高频段的混叠效应逐渐被削弱,在采样率达到10 Hz时饱和;继续提高采样率并不能进一步削弱混叠效应,这一特点在E、N方向上一致,在距震源远近不同的测站也表现一致;同震形变反映在0.01~0.1 Hz处的峰值,N方向较E方向形变更大,故功率谱密度的峰值更高。

      图  5  各采样率的时间序列功率谱密度

      Figure 5.  Power Spectrum Density of Time Series with Different Sampling Rates

      综上所述,50 Hz超高频数据未能显著削弱时间序列存在的混叠效应,其效果与5 Hz、10 Hz接近;在振幅达到分米级时,通过差分的方式可以获得较高精度的速度序列;二次差分获取加速度较为困难,难以取得良好的效果。

    • 本文以2021年青海玛多Mw 7.4地震为例,研究了GNSS多系统组合,尤其是北斗三号对地表位移拾取精度的提升,还分析了50 Hz超高频数据拾取位移、速度、加速度的能力,得到如下结论:

      1)多系统组合对于削弱低频多路径噪声、提高位移拾取稳定性有显著作用。融入了北斗三号的GREC2+3四系统组合效果最优,在E、N、U方向上的精度分别达到了3.4 mm、4.7 mm、12.4 mm,相较于GREC2的组合方式分别提高了29%、13%、23%。

      2)对2021年玛多地震设置1~50 Hz内多个采样率,位移序列混叠效应的削弱在5~10 Hz时达到饱和;采样率在5 Hz时差分速度序列效果最佳;在1~50 Hz,通过二次差分获取加速度的效果均不理想;超高频GNSS在同震位移监测中的作用需要进一步挖掘。

      3)将采样率提升到50 Hz,增大了数据运算、储存、管理的成本(数据量约20 GB/站/天),既未能显著削弱混叠效应,也难以提供有效的超高频速度、加速度序列。5~10 Hz已经能满足大部分场景的需求,加密基准站或升级原有接收机使其可观测北斗三号卫星,将会取得更好的同震形变监测效果。

参考文献 (34)

目录

    /

    返回文章
    返回