GNSS观测的2021年青海玛多地震(Mw 7.4)同震形变及其滑动分布

李志才, 丁开华, 张鹏, 温扬茂, 赵利江, 陈建峰

李志才, 丁开华, 张鹏, 温扬茂, 赵利江, 陈建峰. GNSS观测的2021年青海玛多地震(Mw 7.4)同震形变及其滑动分布[J]. 武汉大学学报 ( 信息科学版), 2021, 46(10): 1489-1497. DOI: 10.13203/j.whugis20210301
引用本文: 李志才, 丁开华, 张鹏, 温扬茂, 赵利江, 陈建峰. GNSS观测的2021年青海玛多地震(Mw 7.4)同震形变及其滑动分布[J]. 武汉大学学报 ( 信息科学版), 2021, 46(10): 1489-1497. DOI: 10.13203/j.whugis20210301
LI Zhicai, DING Kaihua, ZHANG Peng, WEN Yangmao, ZHAO Lijiang, CHEN Jianfeng. Coseismic Deformation and Slip Distribution of 2021 Mw 7.4 Madoi Earthquake from GNSS Observation[J]. Geomatics and Information Science of Wuhan University, 2021, 46(10): 1489-1497. DOI: 10.13203/j.whugis20210301
Citation: LI Zhicai, DING Kaihua, ZHANG Peng, WEN Yangmao, ZHAO Lijiang, CHEN Jianfeng. Coseismic Deformation and Slip Distribution of 2021 Mw 7.4 Madoi Earthquake from GNSS Observation[J]. Geomatics and Information Science of Wuhan University, 2021, 46(10): 1489-1497. DOI: 10.13203/j.whugis20210301

GNSS观测的2021年青海玛多地震(Mw 7.4)同震形变及其滑动分布

基金项目: 

国家自然科学基金 41774011

国家自然科学基金 41874011

国家自然科学基金 41804038

国家重点研发计划 2016YFB0501405

详细信息
    作者简介:

    李志才,博士,正高级工程师,主要从事卫星导航定位及大地测量反演研究。zcli@ngcc.cn

    通讯作者:

    张鹏,博士,正高级工程师。zhangpeng@ngcc.cn

  • 中图分类号: P227

Coseismic Deformation and Slip Distribution of 2021 Mw 7.4 Madoi Earthquake from GNSS Observation

Funds: 

The National Natural Science Foundation of China 41774011

The National Natural Science Foundation of China 41874011

The National Natural Science Foundation of China 41804038

the National Key Research and Development Program of China 2016YFB0501405

More Information
    Author Bio:

    LI Zhicai, PhD,senior engineer, specializes in GNSS data processing and geodesy inversion. E-mail: zcli@ngcc.cn

    Corresponding author:

    ZHANG Peng, PhD, senior engineer, E-mail: zhangpeng@ngcc.cn

  • 摘要: 北京时间2020年5月22日2时,中国青海省果洛州玛多县发生Mw 7.4地震。收集震中附近9个连续运行卫星定位基准站(continuously operating reference stations,CORS)观测数据,基于模糊度解算的精密单点定位(precise point positioning with ambiguity resolution,PPP-AR)技术处理了2 h时段的高频(1 Hz采样)全球导航卫星系统(global navigation satellite system,GNSS)数据,快速确定本次地震的近场同震位移,其中水平方向最大约为0.6 m。结合远场12个CORS震前、震后各3 d低频(30 s采样)数据的非差PPP解算同震形变场,基于弹性位错模型反演了玛多地震断层几何参数和滑动分布。结果显示:玛多地震属于典型走滑事件,发震断层走向278.49°,倾角为64.38°,滑动角为-10.90°,破裂长度约为138.72 km,宽度为4.82 km;滑动量超过3 m的滑移主要集中在东部小于18 km深度的区域,最大破裂可达4.2 m。反演的地震矩为1.85×1020 Nm,相当于矩震级7.45,比美国地质调查局利用地震波的反演结果略大。
    Abstract:
      Objectives  At 2:00 on May 22, 2020 (UTC+8), an Mw 7.4 earthquake struck Madoi County, Qinghai Province, China. The earthquake occurred in the unattended area of the Qinghai Tibet Plateau, with less injuries and building damage, but caused obvious surface rupture.
      Methods  9 high frequency continuously operating reference stations(CORS) data with 1 Hz sampling rate were collected near the epicenter, and the co-seismic deformation was derived from high frequency observation data with 1 hour before and after the earthquake based on precise point positioning with ambiguity resolution(PPP-AR). The low frequency data with 30 s sampling of 12 CORS around Madoi County were collected before and after the earthquake within 3 days to solve for the permeant co-seismic deformation based on the non-differential precise point positioning(PPP).
      Results  The comparison between these two kinds of co-seismic deformation illustrates that the co-seismic deformation from the high frequency (1 Hz) data is slightly smaller than that of the low frequency (30 s) data. The maximum observed global navigation satellite system(GNSS) co-seismic displacement is about 0.6 m.Further, the fault geometry and slip distribution of the earthquake were inverted based on the permanent co-seismic deformation.
      Conclusions  The inversion results indicate that the Madoi earthquake is a typical left-lateral strike slip event (rake angle: -10.90°), with a strike of 278.49°and dip angle of 64.38°. The rupture length is about 138.72 km. The slip distribution reveals that more than 3 m slip is mainly concentrated in the eastern asperity with the depth of less than 18 km with the maximum slip of 4.2 m. The geodetic seismic moment is about 1.85×1020 Nm, equivalent to Mw 7.45, which is slightly larger than the result of United States Geological Survey(USGS) from seismic waves.
  • 据中国地震台网正式测定,北京时间2021年5月22日2时4分,中国青海省果洛州玛多县发生7.4级地震,震源深度17 km,震中位于(98.34°E,34.59°N)[1]。据报道,此次地震可能产生1 m左右的地表破裂。地震导致的剧烈地壳变动信号除被远场地震台以动态波形方式完整记录外,还能被连续运行卫星定位基准站(continuously ope-rating reference stations,CORS)的全球导航卫星系统(global navigation satellite system,GNSS)接收机可靠监测,通过数据传输网络实时传到数据中心,可快速解算确定地震形变。同时,高频(1 Hz)GNSS数据还可以用于快速确定地震发震位置、滑动大小等[2-5]

    随着CORS地震监测技术的不断发展,其在大震应急响应中的作用日渐突出,例如2011年日本东北(Mw 9.0)、2013年中国四川芦山(Ms 7.0)、2013年尼泊尔(Mw 7.8)、2013年美国Craig(Mw 7.5)、2015年智利(Mw 8.3)、2016年新西兰(Mw 7.8)、2018年美国阿拉斯加(Mw 7.9)以及2019年美国加州(Mw 7.1)地震[2-10],密集CORS可有效地确定地震同震形变场,快速开展震源机制确定,也为此后基于多源数据的联合反演提供了参考[11-15]

    中国在十三五期间建设了大量CORS,其中大部分站点为基岩站,采用混凝土基座确保站址稳固[16]。针对此次青海玛多地震,地方自然资源部门、中国地震局及中国移动通讯公司CORS获得了有效的地震观测数据,本文收集了这一区域发震当天的高频(采样率为1 Hz)观测数据,以及地震前后各3 d的低频(采样率30 s)观测数据,经过高精度处理获得该区域的同震形变场信息;结合地震地质资料,进行发震断层参数及滑动分布模型反演,试图探索该区域的地震破裂机制,为今后CORS大震监测与应急响应工作提供借鉴。

    玛多地震发生于青藏高原内部的巴颜喀拉块体,又称松潘-甘孜地块,南北宽约500 km,东西长约2 000 km,总体呈北西-南东走向(见图 1),巴颜喀拉块体被柴达木、羌塘和华南块体所围限,是青藏高原北部的重要次级活动块体[17-18],其北以东昆仑走滑断裂带为界,甘孜-玉树-鲜水河走滑断裂带、龙门山逆冲断裂带为其南、东边界[19]。除块体边缘的深大活动断裂带外,其内部还发育着一系列重要的次级活动断裂,如达日断裂带、白龙江断裂带、虎牙和岷江断裂带以及龙日坝断裂带等[20]。受块体间相互作用、断裂带内部变形及深部中下地壳向东流转等共同作用,巴颜喀拉块体表现出复杂的变形样式与强烈的地震活动性[21-25]。例如东昆仑断裂带自1900年以来共发生5次7级以上大震[26],形成近千千米的地表破裂带(其中仅2001年昆仑山口西Ms 8.1地震的地表破裂就长达426 km),留下东-西大滩和玛沁-玛曲两段地震空区[27]。南甘孜-玉树-鲜水河断裂带也是大震频发,20世纪发生多次强震(如1973年7.6级地震),最新一次是2010年Ms 7.1玉树地震。而龙门山断裂带近来强震活动尤其显著,分别在2008、2010、2017年发生Ms 8.0汶川、Ms 7.0芦山和Ms 7.0九寨沟地震。基于活动构造、震源机制解、GPS站速度、地震破裂展布以及历史地震资料,文献[28]认为巴颜喀拉块体主要边界断裂处于应变加速释放期,发生大震的可能性仍较大。此次地震发生于巴颜喀拉块体内部玛多-甘德断裂[29],此前对该断裂活动构造研究较少,因而从各个方面研究玛多地震具有重大的科学意义。

    图  1  青海玛多地震构造背景
    Figure  1.  Tectonic Settings of Mado Earthquake, Qinghai

    在青海玛多地震震中及周边区域分布有多方建设的GNSS永久性观测站,包括中国大陆构造环境监测网络10个基准站(含玛多县城1站),青海省自然资源厅7个省级CORS站,中国移动通讯公司于2020年建成的4个CORS站(具体分布见图 1)。本文利用其中9站地震当日1 Hz采样率的高频数据,通过模糊度解算的精密单点定位(precise point positioning with ambiguity resolution,PPP-AR)技术计算地震前后的单历元动态解[30];此外,还收集21站地震前后各3 d低频(30 s采样)数据,采用非差PPP技术进行精密处理,获得更稳定的同震形变结果[31]

    本文首先采用武汉大学Pride PPP-AR软件对震中附近200 km范围9站(分别为陆态网络:QHMD;青海省CORS网:MADU、HSHX、KANQ、JDUO;中国移动CORS网:QHAE、QHAG、QHAH、QHAJ)2 h同震观测进行测站动态坐标解算。所用软件基于PPP-AR技术,可实现高频数据的单历元解[30],解算利用武汉大学发布的预报星历,各类数表文件基于最新物理模型,处理结果如图 2所示。需要说明的是,HSHX和QHAH站位于花石峡镇,QHMD、MADU和QHAE 3个站点均位于玛多县城。因而从图 2中可以看出,HSHX和QHAH站变形发生时间和发展趋势在南北和东西方向基本一致,QHMD和QHAE站变形发生时间和发展趋势基本一致,MADU站在震时20 s内,GPS观测卫星仅2颗,无法解算出可靠结果,但基本变形趋势也与QHMD和QHAE站相同。

    图  2  高频GPS观测的单历元动态解
    Figure  2.  Single Epoch Processing Kinematic Solution from High Frequency GPS Observation

    考虑到地震发生时刻高频单历元结果的跳跃性较大,如HSHX、QHMD、QHAE等站东向坐标分量急剧变化,因此本文直接剔除了发震时刻的高动态解,仅保留发震前后较平滑的处理单历元结果分别计算地震前后坐标平均值,通过差分获得可靠同震形变(见表 1图 3)。从表 1中可看出,高频GNSS观测的同震形变在水平方向上误差较小,可靠性较高;垂向误差较大,可靠性略差。QHAE和QHMD的同震水平形变基本一致,略小于MADU;QHAJ变形水平分量接近0.6 m(其中东向变形-0.475 m,北向变形-0.315 m),在所有CORS站中最大,移动方向朝西偏南,垂向变形仅有4.6 cm。

    表  1  基于高频GPS观测的同震形变
    Table  1.  Co-Seismic Deformation form 1 Hz GPS Data
    站代码 纬度/(°) 经度/(°) 东向形变/m 北向形变/m 垂向形变/m 东向误差/m 北向误差/m 垂向误差/m
    HSHX 35.10 98.86 -0.068 0.009 -0.012 0.009 0.010 0.019
    JDUO 34.49 97.98 0.258 -0.058 -0.017 0.010 0.011 0.018
    KANQ 34.91 97.57 -0.159 0.210 -0.021 0.009 0.010 0.018
    MADU 34.92 98.22 -0.325 0.005 0.125 0.020 0.014 0.031
    QHAE 34.92 98.21 -0.239 0.082 0.025 0.009 0.008 0.018
    QHAG 35.11 98.87 -0.066 0.008 0.012 0.009 0.010 0.020
    QHAH 34.33 97.93 0.143 -0.022 0.066 0.012 0.015 0.029
    QHAJ 34.62 99.13 -0.475 -0.315 0.046 0.010 0.011 0.021
    QHMD 34.92 98.21 -0.237 0.079 0.040 0.011 0.015 0.025
    下载: 导出CSV 
    | 显示表格
    图  3  基于低频和高频GNSS观测的同震形变
    Figure  3.  Co-Seismic Deformation from 1 Hz and 30 s Sampling GPS Data

    基于瑞士伯尔尼大学开发的Bernese(V5.2)软件[31],对地震发生前后各3 d 21站的低频GPS观测数据,采用PPP非差技术进行了精密处理。为提高计算效率,实际计算时用300 s间隔数据,卫星截止高度角设为10°,卫星轨道、钟差、地球自转参数、全球电离层数值改正模型均采用欧洲定轨中心(Center for Orbit Determination in Europe,CODE)产出的事后精密数据产品,对流层改正采用Saastamoinen模型。采用最新的各类地球物理模型及其他辅助文件,如大地水准面模型文件EGM2008、章动模型文件IAU2000R06.NUT、IERS常量文件IERS2010XY.SUB、海潮模型FES2004、地球固体潮模型TIDE2000、行星星历等文件。通过精密处理获得了玛多地震附近21个CORS站点的地震前后的依据ITRF2014参考框架的单天解。对地震前后各3 d的单天解3个坐标分量分别取平均,再从震后减去震前平均值,即得高精度的同震形变[13],具体结果参见图 3表 2

    表  2  基于低频GPS观测的同震形变
    Table  2.  Co-Seismic Deformation from 30 s Sampling GPS Data
    站代码 纬度/(°) 经度/(°) 东向形变/m 北向形变/m 垂向形变/m 东向误差/m 北向误差/m 垂向误差/m
    HSHX 35.10 98.86 -0.074 0.011 0.011 0.004 0.004 0.011
    JDUO 34.49 97.98 0.281 -0.058 0.008 0.005 0.006 0.008
    KANQ 34.91 97.57 -0.167 0.204 -0.037 0.018 0.019 0.033
    MADU 34.92 98.22 -0.24 0.088 -0.019 0.005 0.006 0.006
    QHAE 34.92 98.21 -0.237 0.087 -0.002 0.011 0.007 0.009
    QHAG 35.11 98.87 -0.065 0.008 0.002 0.013 0.005 0.003
    QHAH 34.33 97.93 0.15 -0.007 0.008 0.006 0.005 0.004
    QHAJ 34.62 99.13 -0.493 -0.317 0.005 0.009 0.004 0.004
    QHBM 32.93 100.74 0.006 -0.011 -0.008 0.004 0.005 0.009
    QHDL 36.30 98.10 -0.003 0.008 0.002 0.005 0.005 0.006
    QHGE 36.15 94.77 0.000 -0.002 0.003 0.009 0.004 0.003
    QHMD 34.92 98.21 -0.238 0.087 -0.008 0.006 0.006 0.013
    QHME 37.47 101.40 0.005 -0.008 0.004 0.01 0.006 0.005
    QHMQ 34.48 100.25 -0.010 -0.012 0.009 0.004 0.005 0.009
    QHMY 38.48 90.80 0.005 -0.005 0.008 0.005 0.005 0.008
    QHQI 38.09 100.37 -0.009 -0.004 0.021 0.005 0.002 0.018
    QHTT 34.22 92.44 0.002 -0.003 0.010 0.002 0.004 0.014
    QSHE 33.81 97.14 -0.017 0.016 -0.069 0.025 0.017 0.008
    WENQ 35.41 99.44 -0.026 -0.025 0.01 0.002 0.005 0.012
    XNIN 36.60 101.77 -0.001 -0.006 0.002 0.006 0.006 0.011
    XRID 36.05 97.51 -0.008 0.016 0.00 0.006 0.004 0.006
    下载: 导出CSV 
    | 显示表格

    表 2可看出,低频GPS观测的同震形变在水平方向上误差较小,且相对高频数据,垂向误差下降了近一倍。QHMD、MADU和QHAE 3个站点变形在水平量级上保持着良好的一直性,均显示往西约0.24 m、往北约0.09 m的位移,观测精度较高频结果有所提高。低频GPS观测通常可采用更长观测时段,自然较瞬时、高频同震观测可靠性高。同样,在全部CORS站中,QHAJ点受此次地震破裂影响最大,水平分量往西为0.493 m,往南为0.317 m,垂向有5 mm的隆升,与高频结果基本一致。图 3同震形变结果显示玛多地震是一个典型的左旋走滑断层,断层上盘往西北方向运动,下盘向东南方向运动,发震断层破裂方式与青藏高原向东运移一致。

    基于玛多地震的低频GPS同震形变数据,本文对发震断层几何参数及其滑动分布进行反演。反演过程包括两个步骤:首先,采用文献[32]方法进行断层位置、走向、倾角、长度与宽度等参数的非线性反演,反演中采用单纯形算法进行参数搜索,各参数初始值综合考虑中国地震台网中心(China Earthquake Networks Center,CENC)、全球地震矩张量计划(Global Centroid Moment Tensor,GCMT)和美国地质调查局(United States Geological Survey,USGS)等3家机构公布的震源参数[133-34],搜索空间依据初始值误差设置;然后,在参数搜索范围内随机生成模型,并依据这些随机模型正演出地表GNSS站点的同震形变模拟值,进一步与观测值求差以评估所选参数的误差。本文共测试了100组模型,反演了对应的几何参数值及其误差,最终确定的震源几何参数见表 3。从表 3中可以看出反演的断层大致与USGS和GCMT的结果一致,但本文确定的断层位置在经度方向上较GCMT结果略偏西,断层走向为278.49°,倾角为64.38°,滑动角为-10.90°,断层破裂长度为138.72 km,宽度为4.82 km,对应的破裂上表面深度为4.34 km。由于采用了均一滑动模型,反演的断层深度为破裂面上表面深度,因而与地震波反演的震源深度(10~17 km)有较大差异。

    表  3  震源几何参数反演结果
    Table  3.  Focal Mechanism Solutions of the Madu Earthquake
    模型 纬度/(°) 经度/(°) 走向角/(°) 倾角/(°) 滑动角/(°) 长度/km 宽度/km 深度/km
    CENC 34.59 98.246 17
    USGS 34.613 98.246 92 67 -40 10
    200 53 -151
    GCMT 34.65 98.46 13 81 -173 12
    282 83 -9
    本文方法 34.62±0.02 98.38±0.02 278.49±0.83 64.38±12.88 -10.90±5.03 138.72±1.29 4.82±2.53 4.34±2.28
    下载: 导出CSV 
    | 显示表格

    在确定断层几何参数后,本文再对发震断层的滑动分布进行反演,不同于前面步骤中采用均一滑动模型进行反演,反演中各子断层的滑动量是可变的[10]。采用表 3反演的断层位置、走向角、倾角等几何信息,将断层沿走向延展至180 km左右,沿宽度至30 km,进一步将断层划分为62×10个子断层,每个子断层的尺度约为3 km×3 km,以反演各子断层的滑动量,最终的断层滑动分布如图 4所示。

    图  4  断层滑移分布及GNSS同震形变拟合情况
    Figure  4.  Fault Slip Distribution Model and Its Fitting to the GPS Data

    图 4显示断层的破裂沿走向存在两个凹凸体,分别处于GCMT震中的两侧,西边区域小且滑动量较小,约为1~2 m,地表破裂长度约45 km;东部滑动量超过3 m的滑移主要集中在小于18 km深度的区域,最大破裂可达4.2 m,地表破裂长度可达95 km。地震矩为1.85×1020 Nm,相当于矩震级7.45,相较于USGS利用地震波反演的地震矩张量1.31×1020 Nm,矩震级Mw 7.34,本文反演结果略大。

    图 4还展示了低频GNSS计算的同震形变值、模型反演值以及两者残差,残差的均值为1.9 cm,最大值5.1 cm(靠近震中的KANQ站)。相比于同震形变平均值11.1 cm来说,仍有约20%的信号未被同震滑动分布模型所解释。考虑本文仅采用相对简单的断层几何模型,比如断层的走向假设为单一值,实际情况可能会复杂得多。此外,观测到的低频同震形变要比模拟值大,而本文同震形变计算采用了震前、震后各3 d位置的平均值作差求得,由于在震后早期可能存在的震后余滑也被计算为同震滑移,因此残差值加大。同震形变的高频结果要小于静态结果(见图 3)说明了这一点。此外,观测的同震形变可能包含有早期震后余滑,导致地震矩要比USGS等机构给定的大。

    因此,要获得更为准确可靠的断层滑移分布信息,有待加入周边的GPS大地控制点同震位移,以及合成孔径雷达干涉测量(interferometric synthetic aperture radar,InSAR)资料获得的同震位移,同时考虑更为复杂的断层几何结构等。但是,本文基于目前所能搜集到的连续站资料得到的同震位移,反演的同震滑移分布可以较好地反映出同震滑移的大体特征,这对震后早期及时快速了解该地震的总体特征具有重要的参考意义。

    本文针对青海玛多7.4级地震,在震后快速收集了震区附近9个高频GNSS震时观测数据以及远场12站地震前后各3 d的低频GNSS数据,获得高精度同震形变场;并基于位错模式约束发震断层的几何参数和滑动发布。结果表明:

    1) 利用PPP-AR快速解算高动态时间序列可获得较高精度的动态地震位移,可快速确定地震形变区域、大小。CORS监测到最大水平位移近0.6 m,清晰揭示断层左旋走滑活动特征。

    2) 破裂模型显示发震断层走向为278.49°,倾角为64.38°,滑动角为-10.90°,长度为138.72 km,宽度为4.82 km,与利用地震波反演结果基本一致;反演最大滑动约4.2 m,主要集中在发震断层东段小于18 km深度区域,总矩张量1.85×1020 Nm,对应矩震级7.45。

    3) 此次地震的监测实践显示,现有CORS网具备快速确定区域大震(7级以上地震)震源机制和震级的能力。

    致谢: 本文GNSS观测数据来自中国大陆构造环境监测网络、青海省自然资源厅和中国移动公司建设的卫星导航定位基准站,在此表示感谢。
  • 图  1   青海玛多地震构造背景

    Figure  1.   Tectonic Settings of Mado Earthquake, Qinghai

    图  2   高频GPS观测的单历元动态解

    Figure  2.   Single Epoch Processing Kinematic Solution from High Frequency GPS Observation

    图  3   基于低频和高频GNSS观测的同震形变

    Figure  3.   Co-Seismic Deformation from 1 Hz and 30 s Sampling GPS Data

    图  4   断层滑移分布及GNSS同震形变拟合情况

    Figure  4.   Fault Slip Distribution Model and Its Fitting to the GPS Data

    表  1   基于高频GPS观测的同震形变

    Table  1   Co-Seismic Deformation form 1 Hz GPS Data

    站代码 纬度/(°) 经度/(°) 东向形变/m 北向形变/m 垂向形变/m 东向误差/m 北向误差/m 垂向误差/m
    HSHX 35.10 98.86 -0.068 0.009 -0.012 0.009 0.010 0.019
    JDUO 34.49 97.98 0.258 -0.058 -0.017 0.010 0.011 0.018
    KANQ 34.91 97.57 -0.159 0.210 -0.021 0.009 0.010 0.018
    MADU 34.92 98.22 -0.325 0.005 0.125 0.020 0.014 0.031
    QHAE 34.92 98.21 -0.239 0.082 0.025 0.009 0.008 0.018
    QHAG 35.11 98.87 -0.066 0.008 0.012 0.009 0.010 0.020
    QHAH 34.33 97.93 0.143 -0.022 0.066 0.012 0.015 0.029
    QHAJ 34.62 99.13 -0.475 -0.315 0.046 0.010 0.011 0.021
    QHMD 34.92 98.21 -0.237 0.079 0.040 0.011 0.015 0.025
    下载: 导出CSV

    表  2   基于低频GPS观测的同震形变

    Table  2   Co-Seismic Deformation from 30 s Sampling GPS Data

    站代码 纬度/(°) 经度/(°) 东向形变/m 北向形变/m 垂向形变/m 东向误差/m 北向误差/m 垂向误差/m
    HSHX 35.10 98.86 -0.074 0.011 0.011 0.004 0.004 0.011
    JDUO 34.49 97.98 0.281 -0.058 0.008 0.005 0.006 0.008
    KANQ 34.91 97.57 -0.167 0.204 -0.037 0.018 0.019 0.033
    MADU 34.92 98.22 -0.24 0.088 -0.019 0.005 0.006 0.006
    QHAE 34.92 98.21 -0.237 0.087 -0.002 0.011 0.007 0.009
    QHAG 35.11 98.87 -0.065 0.008 0.002 0.013 0.005 0.003
    QHAH 34.33 97.93 0.15 -0.007 0.008 0.006 0.005 0.004
    QHAJ 34.62 99.13 -0.493 -0.317 0.005 0.009 0.004 0.004
    QHBM 32.93 100.74 0.006 -0.011 -0.008 0.004 0.005 0.009
    QHDL 36.30 98.10 -0.003 0.008 0.002 0.005 0.005 0.006
    QHGE 36.15 94.77 0.000 -0.002 0.003 0.009 0.004 0.003
    QHMD 34.92 98.21 -0.238 0.087 -0.008 0.006 0.006 0.013
    QHME 37.47 101.40 0.005 -0.008 0.004 0.01 0.006 0.005
    QHMQ 34.48 100.25 -0.010 -0.012 0.009 0.004 0.005 0.009
    QHMY 38.48 90.80 0.005 -0.005 0.008 0.005 0.005 0.008
    QHQI 38.09 100.37 -0.009 -0.004 0.021 0.005 0.002 0.018
    QHTT 34.22 92.44 0.002 -0.003 0.010 0.002 0.004 0.014
    QSHE 33.81 97.14 -0.017 0.016 -0.069 0.025 0.017 0.008
    WENQ 35.41 99.44 -0.026 -0.025 0.01 0.002 0.005 0.012
    XNIN 36.60 101.77 -0.001 -0.006 0.002 0.006 0.006 0.011
    XRID 36.05 97.51 -0.008 0.016 0.00 0.006 0.004 0.006
    下载: 导出CSV

    表  3   震源几何参数反演结果

    Table  3   Focal Mechanism Solutions of the Madu Earthquake

    模型 纬度/(°) 经度/(°) 走向角/(°) 倾角/(°) 滑动角/(°) 长度/km 宽度/km 深度/km
    CENC 34.59 98.246 17
    USGS 34.613 98.246 92 67 -40 10
    200 53 -151
    GCMT 34.65 98.46 13 81 -173 12
    282 83 -9
    本文方法 34.62±0.02 98.38±0.02 278.49±0.83 64.38±12.88 -10.90±5.03 138.72±1.29 4.82±2.53 4.34±2.28
    下载: 导出CSV
  • [1] 中国地震台网中心. 青海果洛州玛多县7.4级地震[EB/OL]. (2021-06-02). https://news.ceic.ac.cn/CC20210522020411.html

    China Earthquake Network Center. M 7. 4-Madoi, Guoluo, Qinghai, China[EB/OL]. (2021-06-02). https://news.ceic.ac.(cn/CC20210522020411.html

    [2]

    Fang R, Shi C, Song W, et al. Determination of Earthquake Magnitude Using GPS Displacement Waveforms from Realtime Precise Point Positioning[J]. Geophysical Journal International, 2014, 196(1): 461-472 doi: 10.1093/gji/ggt378

    [3]

    Branzanti M, Colosimo G, Crespi M, et al. GPS Near-Real-Time Coseismic Displacements for the Great Tohoku-Oki Earthquake[J]. IEEE Geoscience Remote Sensing Letter, 2013, 10(2): 372-376 doi: 10.1109/LGRS.2012.2207704

    [4]

    Chen K, Ge M, Babeyko A, et al. Retrieving Real-Time Coseismic 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] 李志才, 张鹏, 蒋志浩, 等. 基于单历元方法研究汶川8.0级地震对中国大陆的同震变形[J]. 测绘科学, 2009, 34(5): 47-48

    Li Zhicai, Zhang Peng, Jiang Zhihao, et al. Research on Chinese Continent Coseismic Deformation Effect due to Wenchuan Ms 8.0 Earthquake Based on Single Epoch Method[J]. Science of Surveying and Mapping, 2009, 34(5): 47-48

    [6]

    Jiang Z, Wang M, Wang Y, et al. GPS Constrained Coseismic Source and Slip Distribution of the 2013 Mw 6. 6 Lushan, China, Earthquake and Its Tectonic Implications[J]. Geophysical Research Letters, 2014, 41: 407-413 doi: 10.1002/2013GL058812

    [7]

    Wen Y, Guo Z, Xu C, et al. Coseismic and Postseismic Deformation Associated with the 2018 Mw 7. 9 Kodiak, Alaska, Earthquake from Low-Rate and High-Rate GPS Observations[J]. Bulletin of Seismological Society of America, 2019, 109(3): 908-918 doi: 10.1785/0120180246

    [8]

    Ding K, Freymueller J T, Wang Q, et al. Coseismic and Early Postseismic Deformation of the 5 January 2013 Mw 7. 5 Craig Earthquake from Static and Kinematic GPS Solutions[J]. Bulletin of Seismological Society of America, 2015, 105(2B): 1 153-1 164 doi: 10.1785/0120140172

    [9]

    Zheng K, Zhang X, Li X, et al. Capturing Coseismic Displacement in Real Time with Mixed Singleand Dual-Frequency Receivers: Application to the 2018 Mw 7. 9 Alaska Earthquake[J]. GPS Solutions, 2018, DOI: 10.1007/s10291-018-0794-y

    [10]

    Zhong S, Xu C, Yi L, et al. Focal Mechanisms of the 2016 Central Italy Earthquake Sequence Inferred from High-Rate GPS and Broadband Seismic Wave Forms[J]. Remote Sensing, 2018, 10(4): 1-18 http://www.onacademic.com/detail/journal_1000040544025910_1fdd.html

    [11] 李志才, 张鹏, 温扬茂, 等. 基于GPS和海底基准点观测形变反演日本大地震(Mw 9. 0)同震断层滑动分布[J]. 武汉大学学报∙信息科学版, 2013, 38(1): 40-43 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH201301010.htm

    Li Zhicai, Zhang Peng, Wen Yangmao, et al. Coseismic Slip and Rupture of the 2011 Mw 9. 0 Tohoku Earthquake from GPS and Sea-Floor Point Observations[J]. Geomatics and Information Science of Wuhan University, 2013, 38(1): 40-43 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH201301010.htm

    [12] 李志才, 张鹏, 金双根, 等. 基于GPS观测数据的汶川地震断层形变反演分析[J]. 测绘学报, 2009, 38(2): 108-113 doi: 10.3321/j.issn:1001-1595.2009.02.003

    Li Zhicai, Zhang Peng, Jin Shuanggen, et al. Wen- chuan Earthquake Deformation Fault Inversion and Analysis Based on GPS Observations[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(2): 108-113 doi: 10.3321/j.issn:1001-1595.2009.02.003

    [13]

    Li Z, Wen Y, Zhang P, et al. Joint Inversion of GPS, Leveling, and InSAR Data for the 2013 Lushan(China)Earthquake and Its Seismic Hazard Implications[J]. Remote Sensing, 2020, 12(4): 715 doi: 10.3390/rs12040715

    [14]

    Zhang G, Hetland E, Shan X, et al. Triggered Slip on a Back Reverse Fault in the Mw 6. 8 2013 Lushan, China Earthquake Revealed by Joint Inversion of Local Strong Motion Accelerograms and Geodetic Measurements[J]. Tectonophysics, 2016, 672/ 673: 24-33

    [15]

    Huang Y, Qiao X, Freymueller J T, et al. Fault Geometry and Slip Distribution of the 2013 Mw 6. 6 Lushan Earthquake in China Constrained by GPS, InSAR, Leveling, and Strong Motion Data[J]. Journal of Geophysical Research Solid Earth, 2019, 124: 7 341-7 353 doi: 10.1029/2019JB017451

    [16] 陈明, 武军郦, 李志才. 国家卫星导航定位基准站建设[J]. 地理信息世界, 2018, 25(1): 42-46 doi: 10.3969/j.issn.1672-1586.2018.01.009

    Chen Ming, Wu Junli, Li Zhicai. Construction of the National GNSS Continuous Operation Reference Stations[J]. Geomatics World, 2018, 25(1): 42-46 doi: 10.3969/j.issn.1672-1586.2018.01.009

    [17] 邓起东, 张培震, 冉勇康, 等. 中国活动构造基本特征[J]. 中国科学(D辑): 地球科学, 2002, 32(12): 1 020-1 030 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200212006.htm

    Deng Qidong, Zhang Peizhen, Ran Yongkang, et al. Basic Characteristics of Chinese Active Tectonics [J]. Science in China(Series D): Earth Sciences, 2002, 32(12): 1 020-1 030 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200212006.htm

    [18] 张培震, 邓起东, 张国民, 等. 中国大陆的强震活动与活动地块[J]. 中国科学(D辑): 地球科学, 2003, 33(S): 12-20 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK2003S1001.htm

    Zhang Peizhen, Deng Qidong, Zhang Guomin, et al. Active Tectonic Blocks and Strong Earthquake Activities in the Continent of China[J]. Science in China(Series D): Earth Sciences, 2003, 33(S): 12-20 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK2003S1001.htm

    [19] 邓起东, 高翔, 陈桂华, 等. 青藏高原昆仑-汶川地震系列与巴颜喀喇断块的最新活动[J]. 地学前缘, 2010, 17(5): 163-178 https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201005017.htm

    Deng Qidong, Gao Xiang, Chen Guihua, et al. Recent Tectonic Activity of Bayankala Fault-Block and the Kunlun-Wenchuan Earthquakes Series of the Tibetan Plateau[J]. Earth Science Frontiers, 2010, 17(5): 163-178 https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201005017.htm

    [20] 徐锡伟, 闻学泽, 陈桂华, 等. 巴颜喀拉地块东部龙日坝断裂带的发现及其大地构造意义[J]. 中国科学(D辑): 地球科学, 2008, 38(5): 529-542 doi: 10.3321/j.issn:1006-9267.2008.05.001

    Xu Xiwei, Wen Xueze, Chen Guihua, et al. Discovery of Longriba Fault Zone in Eastern Bayankala Block, China and Its Tectonic Implication[J]. Science in China(Series D): Earth Sciences, 2008, 38(5): 529-542 doi: 10.3321/j.issn:1006-9267.2008.05.001

    [21]

    Gan W, Zhang P, Shen Z K, et al. Present-Day Crustal Motion Within the Tibetan Plateau Inferred from GPS Measurements[J]. Journal of Geophysical Research Solid Earth, 2007, 112(B8): 582-596 http://perso-sdt.univ-brest.fr/~jacdev/pdf/gan07_Tibet_GPS.pdf

    [22] 王椿镛, 楼海, 吕智勇, 等. 青藏高原东部地壳上地幔S波速度结构: 下地壳流的深部环境[J]. 中国科学(D辑): 地球科学, 2008, 38(1): 22-32 doi: 10.3321/j.issn:1006-9267.2008.01.003

    Wang Chunyong, Lou Hai, Lü Zhiyong, et al. S-wave Crustal and Upper Mantle's Velocity Structure in the Eastern Tibetan Plateau-Deep Environment of Lower Crustal Flow[J]. Science in China: Series D(Earth Sciences), 2008, 38(1): 22-32 doi: 10.3321/j.issn:1006-9267.2008.01.003

    [23]

    Bai D, Unsworth M J, Meju M J, et al. Crustal Deformation of the Eastern Tibetan Plateau Revealed by Magnetotelluric Imaging[J]. Nature Geoscience, 2010, 3(5): 358-362 doi: 10.1038/ngeo830

    [24] 李冲, 李建成, 黄瑞金, 等. 青藏高原东部地壳物质流变模型及汶川地震机理探讨[J]. 武汉大学学报∙信息科学版, 2015, 40(6): 810-815 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH201506020.htm

    Li Chong, Li Jiancheng, Huang Ruijin, et al. Discussion of Crustal Flow Beneath the Eastern Tibetan Plateau and Mechanism of the Wenchuan Earthquake[J]. Geomatics and Information Science of Wuhan University, 2015, 40(6): 810-815 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH201506020.htm

    [25] 中国地震局震害防御司. 中国近代地震目录(公元1912— 1990年)[M]. 北京: 地震出版社, 1999

    Department of Earthquake Disaster Prevention, State Seismological Bureau. The Catologure of Chinese Historical Strong Earthquakes 1912—1990)[M]. Beijing: Seismological Press, 1999

    [26]

    Xu X, Chen W, Ma W, et al. Surface Rupture of the Kunlunshan Earthquake(Ms 8. 1), Northern Tibetan Plateau, China[J]. Seismological Research Letter, 2002, 73(6): 884-892 doi: 10.1785/gssrl.73.6.884

    [27] 闻学泽, 杜方, 张培震, 等. 巴颜喀拉块体北和东边界大地震序列的关联性与2008年汶川地震[J]. 地球物理学报, 2011, 54(3): 706-716 doi: 10.3969/j.issn.0001-5733.2011.03.010

    Wen Xueze, Du Fang, Zhang Peizhen, et al. Correlation of Major Earthquake Sequences on the Northern and Eastern Boundaries of the Bayan-Har Block, and Its Relation to the 2008 Wenchuan Earthquake[J]. Chinese Journal of Geophysics, 2011, 54(3): 706-716 doi: 10.3969/j.issn.0001-5733.2011.03.010

    [28] 高翔, 邓起东. 巴颜喀喇断块边界断裂强震活动分析[J]. 地质学报, 2013, 87(1): 9-19 doi: 10.3969/j.issn.0001-5717.2013.01.002

    Gao Xiang, Deng Qidong. Activity Analysis of Large Earthquake in Boundary Faults Around the Bayankaka Faulting Block[J]. Acta Geologica Sinica, 2013, 87(1): 9-19 doi: 10.3969/j.issn.0001-5717.2013.01.002

    [29] 詹艳, 梁明剑, 孙翔羽, 等. 2021年5月22日青海玛多Ms 7. 4级地震深部环境及发震构造模式[J]. 地球物理学报, 2021, 64(7): 2 232-2 252 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX202107002.htm

    Zhan Yan, Liang Mingjian, Sun Xiangyu, et al. Deep Structure and Seismogenic Pattern of the 2021, 5. 22 Madoi(Qinghai)Ms 7. 4 Earthquake[J]. Chinese Journal of Geophysics, 2021, 64(7): 2 232-2 252 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX202107002.htm

    [30]

    Geng J, Chen X, Pan Y, et al. PRIDE PPP-AR: An Open-Source Software for GPS PPP Ambiguity Resolution[J]. GPS Solutions, 2019, 23: 91 doi: 10.1007/s10291-019-0888-1

    [31]

    Dach R, Lutz S, Walser P, et al. Bernese GNSS Software Version 5. 2. User Manual[R]. Switzerland: Bern Open Publishing, 2015

    [32]

    Clarke P J, Paradissis D, Briole P, et al. Geodetic Investigation of the 13 May 1995 Kozani-Grevena (Greece) Earthquake[J]. Geophysical Research Letters, 1997, 24(6): 707-710 http://www.onacademic.com/detail/journal_1000035777530210_0e63.html

    [33]

    USGS. M 7. 3-Southern Qinghai, China[EB/ OL]. [2021-06-02]. https://earthquake.usgs.gov/earthquakes/eventpage/us7000e54r/executive

    [34]

    GCMT. Global CMT Catalog Search[EB/OL]. [2021-06-02]. https://www.globalcmt.org/CMT-search.html

  • 期刊类型引用(33)

    1. 孟思晨,孟秋,陈启志,胡才博. 2021年5月22日青海玛多M_S7.4地震同震变形的数值模拟及其动力学启示. 中国科学院大学学报. 2024(01): 70-80 . 百度学术
    2. 李思辉,刘超,张开坤,王向阳,段伟,孙景强. 基于时间卷积模型的GNSS多路径误差削弱. 江西科学. 2024(01): 142-147 . 百度学术
    3. 余娜,王建军,张晓清,杨晓霞. 临震微波动现象回溯性研究——以青海地区为例. 大地测量与地球动力学. 2024(04): 336-340 . 百度学术
    4. 高志钰,李彦川,单新建,黄星. 基于模拟高频GNSS数据的安宁河断裂强震滑动分布快速反演研究. 武汉大学学报(信息科学版). 2024(09): 1527-1537 . 百度学术
    5. 赵忠海,岳昊,关沧海. 基于PRIDE PPP-AR的精密单点定位精度分析. 测绘与空间地理信息. 2024(09): 37-39+43 . 百度学术
    6. 吴凯,朱守彪. 2021年玛多M_W7.4地震超剪切破裂的动力学机制. 大地测量与地球动力学. 2024(12): 1264-1270 . 百度学术
    7. 张彩红,谭凯,鲁小飞,李琦,李承涛,李志才. 玛多M_W7.3地震同震滑动分布的三维有限元模拟分析. 大地测量与地球动力学. 2023(01): 42-45+60 . 百度学术
    8. 王运生,许强,赵伟华,吉锋,杨莹辉,黄健,崔圣华,董秀军,李为乐. 青海玛多7.4级地震发震断裂继承性及新生性探讨. 灾害学. 2023(01): 57-62+71 . 百度学术
    9. 王利,许豪,舒宝,李新瑞,田云青. 滑坡灾害GNSS监测技术研究进展与展望. 导航定位与授时. 2023(01): 12-26 . 百度学术
    10. 王浩浩,郝明,庄文泉. GNSS实时卫星钟差估计在地震监测中的应用. 导航定位与授时. 2023(03): 108-116 . 百度学术
    11. 谭明明,樊亚灵,吕柯欣,蒿赏,徐涛,蔡国田,张杰,李志才. 不同卫星导航系统组合监测的玛多地震同震形变与精度分析. 测绘通报. 2023(07): 80-84 . 百度学术
    12. 王帅鹏,徐克科,王小怡. 2021年玛多M_W7.4地震强震加载及其对周围地区应力扰动的影响. 地震学报. 2023(05): 875-891 . 百度学术
    13. Shuai Wang,Chuang Song,ShanShan Li,Xing Li. Resolving co-and early post-seismic slip variations of the 2021 M_W 7.4 Madoi earthquake in east Bayan Har block with a block-wide distributed deformation mode from satellite synthetic aperture radar data. Earth and Planetary Physics. 2022(01): 108-122 . 必应学术
    14. 余鹏飞,熊维,陈威,乔学军,王迪晋,刘刚,赵斌,聂兆生,李瑜,赵利江,张怀. 基于GNSS和InSAR约束的2021年玛多M_S7.4地震同震滑动分布及应用. 地球物理学报. 2022(02): 509-522 . 百度学术
    15. 王阅兵,李瑜,蔡毅,蒋连江,师宏波,江在森,甘卫军. GNSS观测的2021年5月22日玛多M_S7.4地震同震位移及其约束反演的滑动破裂分布. 地球物理学报. 2022(02): 523-536 . 百度学术
    16. 王迪晋,王东振,赵斌,李瑜,赵利江,王阅兵,聂兆生,乔学军,王琪. 2021年青海玛多M_W7.4地震GNSS同震形变场及其断层滑动分布. 地球物理学报. 2022(02): 537-551 . 百度学术
    17. 刘博研,解孟雨,史保平. 青海玛多M_S7.4地震对周边活动断裂的库仑应力加载及发震概率增量的计算. 地球物理学报. 2022(02): 563-579 . 百度学术
    18. 张岚,李琦,唐河,孙文科. 2021年青海玛多M_W7.5地震的同震变形及断层滑动模型研究. 地球物理学报. 2022(03): 1044-1056 . 百度学术
    19. Yashan FENG,Xiong XIONG,Bin SHAN,Chengli LIU. Coulomb stress changes due to the 2021 M_S7.4 Maduo Earthquake and expected seismicity rate changes in the surroundings. Science China(Earth Sciences). 2022(04): 675-686 . 必应学术
    20. Jihong LIU,Jun HU,Zhiwei LI,Zhangfeng MA,Lixin WU,Weiping JIANG,Guangcai FENG,Jianjun ZHU. Complete three-dimensional coseismic displacements due to the 2021 Maduo earthquake in Qinghai Province, China from Sentinel-1 and ALOS-2 SAR images. Science China(Life Sciences). 2022(04): 687-697 . 必应学术
    21. Yashan FENG,Xiong XIONG,Bin SHAN,Chengli LIU. Coulomb stress changes due to the 2021 M_S7.4 Maduo Earthquake and expected seismicity rate changes in the surroundings. Science China(Life Sciences). 2022(04): 675-686 . 必应学术
    22. 刘计洪,胡俊,李志伟,马张烽,吴立新,姜卫平,冯光财,朱建军. 融合哨兵1号和ALOS-2数据的2021年青海玛多地震高精度三维同震形变场研究. 中国科学:地球科学. 2022(05): 882-892 . 百度学术
    23. 冯雅杉,熊熊,单斌,刘成利. 2021年玛多M_S7.4地震导致的周边地区库仑应力加载及地震活动性变化. 中国科学:地球科学. 2022(06): 1100-1112 . 百度学术
    24. 江颖,刘子维,张晓彤,张丽娜,韦进. 2021年青海玛多Mw 7.4地震前后b值的变化特征研究. 武汉大学学报(信息科学版). 2022(06): 907-915 . 百度学术
    25. 柴海山,陈克杰,魏国光,方荣新,邹蓉,祝会忠. 北斗三号与超高频GNSS同震形变监测:以2021年青海玛多Mw 7.4地震为例. 武汉大学学报(信息科学版). 2022(06): 946-954 . 百度学术
    26. Yanbo Zhang,Yueren Xu,Wenqiao Li,Runchao Liu,Ruoyu Mu,Jiayi Li,Da Zhang,Haofeng Li,Qinjian Tian. Coseismic surface rupture characteristics and earthquake damage analysis of the eastern end of the 2021 M_S7.4 Madoi(Qinghai) earthquake. Earthquake Research Advances. 2022(02): 49-58 . 必应学术
    27. 岳冲,孙玺皓,牛安福,吉平,李晓帆,赵静. 2021年青海玛多M_S7.4地震同震应力场与形变台站同震阶变分布关系探讨. 中国地震. 2022(02): 189-198 . 百度学术
    28. 何福秀,杨洋,饶太国,陈聪,刘江,张澜,罗凌,姚贤良,张翼. 四川区域GNSS网络观测到玛多7.4级地震引起的远场同震位移. 四川地震. 2022(02): 9-13+20 . 百度学术
    29. 吴立新,卢菁琛,毛文飞,胡俊,周子龙,李志伟,齐源,姚汝冰. 基于断层倾角分段变化的玛多地震发震断层构造应力场演化数值模拟分析. 地球物理学报. 2022(10): 3844-3857 . 百度学术
    30. 熊仕昭,刘刚,林牧,王迪晋,熊维,赵利江,赵斌,乔学军. 基于GPS的2021年玛多M_W7.4地震同震静、动态形变场. 大地测量与地球动力学. 2022(10): 1041-1046 . 百度学术
    31. 王守文,季灵运,朱良玉,刘传金. 基于InSAR技术监测2021年青海玛多M_w7.4级地震同震形变场与断层滑动. 地球科学与环境学报. 2022(06): 1016-1026 . 百度学术
    32. 杨君妍,孙文科,洪顺英,苑争一,李瑜,陈伟,孟国杰. 2021年青海玛多7.4级地震的同震变形分析. 地球物理学报. 2021(08): 2671-2683 . 百度学术
    33. 曹泽林,陶夏新,陶正如. 2021年玛多7.4级地震近断裂三分量地震动场合成. 世界地震工程. 2021(04): 1-11 . 百度学术

    其他类型引用(15)

图(4)  /  表(3)
计量
  • 文章访问数:  2018
  • HTML全文浏览量:  925
  • PDF下载量:  287
  • 被引次数: 48
出版历程
  • 收稿日期:  2021-06-07
  • 发布日期:  2021-10-04

目录

/

返回文章
返回