青海玛多7．4级地震GNSS同震水平位移的快速获取

第４３卷　第５期２０２１年１０月
地　震　地　质
ＳＥＩＳＭＯＬＯＧＹＡＮＤＧＥＯＬＯＧＹ
Ｖｏｌ．４３，Ｎｏ．５Ｏｃｔ．，２０２１
ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５３－４９６７．２０２１．０５．００３
李经纬，陈长云，占伟，等．２０２１．青海玛多７．４级地震ＧＮＳＳ同震水平位移的快速获取［Ｊ］．地震地质，４３（５）：１０７３—１０８４．ＬＩＪｉｎｇ ｗｅｉ
，ＣＨＥＮＣｈａｎｇ ｙｕｎ，ＺＨＡＮＷｅｉ，ｅｔａｌ．２０２１．ＲｅｓｅａｒｃｈｏｎｆａｓｔａｃｑｕｉｓｉｔｉｏｎｏｆＧＮＳＳｃｏｓｅｉｓｍｉｃｈｏｒｉｚｏｎｔａｌｄｉｓｐｌａｃｅｍｅｎｔｏｆＭａｄｕｏＭＳ７ ４ｅａｒｔｈｑｕａｋｅｉｎＱｉｎｇｈａｉＰｒｏｖｉｎｃｅ［Ｊ］．
ＳｅｉｓｍｏｌｏｇｙａｎｄＧｅｏｌｏｇｙ，４３（５）：１０７３—１０８４．青海玛多７．４级地震ＧＮＳＳ同震
水平位移的快速获取
李经纬１，２）　陈长云１）
　占　伟１）　武艳强１
）１）中国地震局第一监测中心，天津　３００１８０
２
）中国科学院精密测量科学与技术创新研究院，武汉　４３００７１摘　要 ２０２１年５月２１日青海玛多７ ４级地震导致玛多县及其周边地区受灾严重。

为了快速获取Ｇ
ＮＳＳ同震位移场服务于震情快速研判，文中采用多种方法对玛多地震的同震位移开展了针对性研究。

采用１Ｈｚ高频动态差分、震后５ｈ静态解算及多天静态解算３种方法获取了同震位移，经对比分析发现ＧＡＭＩＴ／ＧＬＯＢＫ的动态差分ＴＲＡＣＫ模块在利用高频数据获取同震形变时适用于近场形变量较大的测站，对远场形变量小的测站则适用性不强；震后５ｈ静态解算获得的同震位移与多天静态解算结果相比时效性更强，两者的数值相差１～４ｍｍ，所反映的位移大小、方向及影响范围都保持较强的一致性。

综合ＧＮＳＳ资料分析表明，本次地震的发震断层为昆仑山口－江错断裂，且震源机制以左旋走滑为主，距离震中区约４０ｋｍ的ＱＨＭＤ站（青海玛多）的Ｗ向位移约为２４ｃｍ，Ｎ向位移约为８ｃｍ。

从近场震后５ｈ静态解和多天静态解的结果分析可知震后变形显著，仍需持续关注。

关键词 青海玛多７ ４级地震　ＧＮＳＳ　同震位移中图分类号：Ｐ３１５ ２
文献标识码：Ａ
文章编号：０２５３－４９６７（２０２１）０５－１０７３－１２
〔收稿日期〕　２０２１－０５－３１收稿，２０２１－０８－１１改回。

〔基金项目〕　国家自然科学基金（４１９７４０１１，１１８０３０３２）、国家重点研发计划项目（２０１８ＹＦＣ１５０３６００）和中国地
震局震情跟踪定向工作任务（２０２１０１０２２２５）共同资助。

通讯作者：陈长云，男，１９８１年生，高级工程师，主要研究方向为活动构造与现今地壳变形，Ｅ－ｍａｉｌ：ｃｃｙ＿６６６＠１６３ ｃｏｍ。

０　引言
据中国地震台网正式测定，２０２１年５月２２日２时４分，在青海省果洛藏族自治州玛多县（３４ ５９°Ｎ，９８ ３４°Ｅ）发生７ ４级地震，震源深度１７ｋｍ，此后发生了多次４～５级余震（图１ｂ）。

经专家学者现场考察及地震序列的精定位初步确定，此次地震发生昆仑山口－江错断裂（李智敏等，２０２１；潘家伟等，２０２１；王未来等，２０２１）（图１）。

昆仑山口－江错断裂位于巴颜喀拉块体北边界———东昆仑断裂带以南，向Ｅ与东昆仑断裂衔接，向Ｗ延伸至甘德附近，全长＞５００ｋｍ，位于玛多－甘德断裂与甘德南缘断裂之间（邓起东等，２００７；徐锡伟等，２０１６；李智敏等，２０２１）。

该断裂所在地区发育一系列大型活动断层，包括西藏大沟－昌马河断裂、玛多－甘
地　震　地　质４３卷
①ｈｔｔｐｓ：∥ｅａｒｔｈｑｕａｋｅ．ｕｓｇｓ．ｇｏｖ／ｅａｒｔｈｑｕａｋｅｓ。

②ｆｔｐ：∥６０．３０．７７．１９。

图１　玛多７ ４级地震的构造背景及余震分布
Ｆｉｇ．１　ＴｈｅｔｅｃｔｏｎｉｃｓｅｔｔｉｎｇｏｆｔｈｅＭａｄｕｏＭＳ７ ４ｅａｒｔｈｑｕａｋｅａｎｄｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆａｆｔｅｒｓｈｏｃｋｓ．
ａ青海玛多７ ４级地震震中区历史地震、活动断裂及ＧＮＳＳ速度场；ｂ青海玛多７ ４级地震的主震及余震分布。

断层数据来自于 邓起东等（２００７）和徐锡伟等（２０１６），同震和余震数据来自于ＵＳＧＳ①，ＧＮＳＳ速度场结果由中国地震局第一监测中心计算②
德断裂、甘德南缘断裂及达日断裂等。

第四纪以来该地区地震活动性显著，近百年来发生了多次大地震。

例如，此次地震震中Ｎ
Ｗ的东昆仑断裂在１９３７年和１９６３年发生了２次７级以上地震、ＳＥ的达日断裂中段在１９４７年发生了Ｍ７ ７地震（中国地震局震害防御司，１９９９）。

该地区的ＧＮＳＳ在ＩＴＲＦ参考框架下呈现Ｎ偏Ｅ的运动趋势，ＥＷ向运动分量由西向东逐渐增大，由南向北逐渐减小。

研究表明，北部东昆仑断裂带的走滑速率自西向东由１３～１４ｍｍ／ａ逐步减弱至３ｍｍ／ａ，其中托索湖段约为５ｍｍ／ａ（高翔等，２０１３；梁明剑等，２０２０）。

Ｚｈｕ等（２０２１）通过２００５—２０１０年的ＧＮＳＳ观测数据给出昆仑山口－江错的滑动速率约为（１ ２±０ ８）ｍｍ／ａ。

虽然该断层的运动速率小于附近的东昆仑托索湖段运动速率，但这次地震的发生证明该地区的构造形变发生在整个区域而不是少数几条主要断层上（张培震等，２
００３）。

４７０１
５期李经纬等：青海玛多７．４级地震ＧＮＳＳ同震水平位移的快速获取同震近、远场位移场是认识发震断裂运动学、动力学特征必不可少的资料（Ｂａｎｅｒｊｅｅｅｔａｌ．，２００５；Ｇａｎｅｔａｌ．，２００７；Ｗｕｅｔａｌ．，２０１６；Ｗａｎｇｅｔａｌ．，２０２０），快速产出该资料对于震后快速判定未来地震趋势、研究区域地壳变形特征及地震动力学特征具有极其重要的意义。

快速获取ＧＮＳＳ测站的坐标常采用动态数据解算的方法，可在发震后利用１５ｍｉｎ的高频ＧＮＳＳ观测数据获取，常用的软件包括各类精密单点定位（ＰＰＰ）软件及ＧＡＭＩＴ／ＧＬＯＢＫ的ＴＲＡＣＫ动态差分处理模块①（
Ｐｓｉｍｏｕｌｉｓｅｔａｌ．，２０１６）。

由于受到星历、钟差等产品精度的限制，ＰＰＰ的精度常为ｃｍ级（Ｇｕｏｅｔａｌ．，２０１５；Ｓｕｅｔａｌ．，２０２０）。

ＴＲＡＣＫ模块采用双差卡尔曼滤波方法可消除部分误差，但其主要适用于中短距离的基线，定位精度也在ｃｍ级（殷海涛等，２０１２；Ｔｉｒｙａｋｉｏｇｌｕｅｔａｌ．，２０１７）。

除了采用ＧＮＳＳ高频数据动态定位外，可利用４～８ｈ观测数据获取ｍｍ级精度的定位结果，这在测绘工程的高等级控制网观测中已有众多案例（周命端等，２０１５）。

因此，可利用地震前、后数ｈ的静态观测数据采用双差的静态解算模式快速获取同震位移结果。

本文利用震中所在的巴颜喀拉块体内部及周边１８个连续ＧＮＳＳ站点的数据（其中７个测站记录了１Ｈｚ高频数据），采用震前、后１５ｍｉｎ１Ｈｚ高频数据动态解算，震后５ｈ静态解算和多天静态解算３种方法获取了同震位移并进行对比分析，探索和测试了采用震前、震后数ｈ的观测资料静态解算获取形变资料的可行性，该测试成果快速揭示了青海玛多７ ４级地震产生的近、远场同震位移，为灾情评估和地震趋势的快速判定提供了基础资料，有助于加深对青海玛多７ ４级地震发生机理的认识。

１　ＧＮＳＳ数据及处理方法
地震发生后，中国地震局台网中心及时发布了周边地区的“陆态网络”ＧＮＳＳ连续站的原始观测数据。

本文挑选了震中附近的１８个ＧＮＳＳ连续观测站的数据。

观测数据可分为２类，一类是采样间隔为３０ｓ的观测数据，包含同震当天、震前３ｄ和震后１ｄ的数据；另一类是其中７个测站记录的地震当天１Ｈｚ的高频观测数据。

针对这２类数据，分别采用静态解算和动态解算方法进行计算。

在进行静态解算时，采用Ｇ
ＡＭＩＴ／ＧＬＯＢＫ网解模式。

为满足网解需求，本文选取了１８个中国大陆周边的ＩＧＳ站作为框架参考站。

首先利用ＧＡＭＩＴ／ＧＬＯＢＫ软件获取ＧＮＳＳ基线松弛解，然后以ＩＧＳ站点作为框架点采用ＧＬＯＢＫ进行约束平差获得ＩＴＲＦ１４参考框架下的测站坐标①。

这种双差模式的优势在于能够在双差组合方程中消除钟差、电离层延迟等带来的误差影响。

在动态解算中，本文采用ＧＡＭＩＴ／ＧＬＯＢＫ的ＴＲＡＣＫ模块进行计算，其可利用震前、震后约１０ｍｉｎ的高频数据获取地震发生瞬间的动态变化过程。

在解算中首先需要找一个基准点，然后进行载波相位动态差分解算，这种方法获取结果的速度较快，但基线＞１００ｋｍ时的解算精度为ｃｍ级（范士杰等，２０１３）。

无论是进行静态解算还是动态解算，都是为了快速而准确地获取地震前后的坐标，从而获取同震位移。

在获取震前、震后坐标时，根据不同的解算模式及所用数据的时长，将其划分为３种方法：１）利用１Ｈｚ高频观测数据进行动态差分定位获取逐历元解，然后分别在震前、震后取平均作差获得最终同震位移，这种方法可以在获取高频数据后快速产出结果；２
）仅采用震后数ｈ的观测值进行静态解算，解算前需要将单天数据中的震前、震后数据分割为独立的文件
５
７０１①
ＨｅｒｒｉｎｇＴ，ＫｉｎｇＲ，ＭｃｃｌｕｓｋｙＳ，等，２０１８，ＧＡＭＩＴ／ＧＬＯＢＫｒｅｆｅｒｅｎｃｅｍａｎｕａｌｓ，Ｒｅｌｅａｓｅ１０ ７，ＭＩＴＴｅｃｈｎｉｃａｌＲｅｐｏｒｔｓ。

地　震　地　质４３卷
进行静态解算，针对此次青海玛多７ ４级地震，本文将同震当天的数据去除ＵＴＣ１８：００—１８：３０时刻的数据后分割为２个观测文件，同时为了验证５ｈ结果的可靠性，对震前１８ｈ的文件也进行了分割处理；３
）利用多天的数据获取精确解，本文选取了震前４ｄ和震后１ｄ获取的稳定结果，由于该地区在ＩＴＲＦ参考框架下，变化最大的方向量值≤５０ｍｍ／ａ，即单天变化约为０ １ｍｍ，因此认为在震前３～４ｄ测站是固定的，不存在较大的线性速率问题，在获取震前坐标时直接采用４
ｄ的平均值。

图２　利用高频ＧＮＳＳ数据获取的同震位移及分布（参考站为ＧＳＭＡ）
Ｆｉｇ．２　Ｃｏｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔａｎｄｄｉｓｔｒｉｂｕｔｉｏｎｇａｉｎｅｄｆｒｏｍｈｉｇｈ ｆｒｅｑｕｅｎｃｙＧＮＳＳｄａｔａ（ｒｅｆｅｒｅｎｃｅｓｔａｔｉｏｎＧＳＭＡ）．
ａ利用高频ＧＮＳＳ数据通过动态差分获取的同震位移（相对于ＧＳＭＡ站）分布；ｂＧＳＭＡ－ＱＨＭＤ基线的长度变化；ｃＧＳＭＡ－ＱＨＢＭ基线的长度变化。

图ｂ、ｃ中的灰色虚线表示玛多地震发生时刻所在历元，蓝色实线表示基线坐标
差值序列，红色实线表示震前、震后各自的中位值线
２　结果与讨论
２ １　１Ｈｚ高频数据的动态差分方法结果
依据ＴＲＡＣＫ模块的解算特点，其获取的结果为２点之间基线长度的动态变化值，因此对应地选取ＵＴＣ１７：４５—１８：１５时段的数据进行动态差分解算，卫星星历为ＩＧＳ提供的快速星历（
ｉｇｒ）。

图２ａ展示了ＧＳＭＡ－ＱＨＭＱ测站的结果，基线向Ｗ和向Ｓ分别发生了近１ｃｍ和近２ｃｍ的位移。

图２ｂ展示了ＸＺＣＤ－ＱＨＭＤ测站的结果，从结果中可以明显地看出地震发生过程中在近断层处的位移过程，基线向Ｅ发生了近２２ｃｍ的位移，向Ｎ发生了近６ｃｍ的位移。

图２ｃ展示了ＧＳＭＡ－ＱＨＢＭ测站结果，基线显示向Ｅ发生了近１ ４ｃｍ位移。

利用此结果结合地质构造背景可以基本判断此次地震发生于左旋走滑型断层上。

从图２中我们发现，由于Ｔ
ＲＡＣＫ模块采用卡尔曼滤波方法，解算中历元间值存在关联，６７０１
５期李经纬等：青海玛多７．４级地震ＧＮＳＳ同震水平位移的快速获取在地震瞬时历元附近的序列会抖动，作为离震中较远的基准站ＧＳＭＡ受到地震影响的时间晚于其他站，因此在基线时间序列中存在２次抖动的特征。

由于ＧＳＭＡ站与ＱＨＢＭ站距离较近，２次抖动存在重叠，同时基线Ｎ向位移量值较小，在图２ｃ中显示不明显。

在解算精度方面，经统计基线解算收敛后的整体解算精度为１～２ｃｍ，除上述几个测站外其他基线的变化量值较小，离散度已经超过了变化量值，较难采用平均值或中值的办法来获取准确的位移量。

除此之外，从地震发生时刻到地震波传播到测站位置时存在时间差，由于数据离散程度较大，当位移量较小时较难从数据的角度确定准确的到达时刻。

因此，Ｔ
ＲＡＣＫ模块受制于精度的影响，对于本次解算中位移量较大的ＱＨＭＤ、ＱＨＭＱ站可以采用这种方法合理地估算出位移量，而对于其他测站则较难获取合理的位移量值，但从位移的趋势与分布可以准确地呈现断层两侧的运动现象，在快速揭示破裂性质方面仍有着重要作用。

２ ２　震后５ｈ数据的静态解算方法结果
此次地震发生在５月２１日ＵＴＣ１８：０４：１３时，在震后第２天（北京时间５月２２日）获取了相关数据，地震当天的数据中包括震后５ｈ的观测值，满足ＧＮＳＳ静态双差处理的时长需求。

卫星星历采用Ｉ
ＧＳ中心发布的快速星历（ｉｇｒ）。

采用ＧＦＺＲＮＸ软件①对当天的数据进行分时段提取，提取时段分为第１时段（ＵＴＣ００：００—１８：００）、第２时段（ＵＴＣ００：００—０９：００）、第３时段（
ＵＴＣ０９：００—１８：００）及第４时段（ＵＴＣ１８：３０—２４：００）。

同时假定ＵＴＣ００：００—１８：００时段的结果为真，以此为基准经处理后获得的每个测站的震前水平方向坐标差值见表１，第４段与第１时段的差值表示ＧＮＳＳ的同震位移结果，如图３所示。

表１中，第１时段和第２时段结果与第４时段坐标结果差异的中位值为０ ３ｍｍ。

除了青海的ＱＨＹＳ站和西藏的ＸＺＣＤ站及震中距约为４０ｋｍ的ＱＨＭＤ站（青海玛多）在震前的结果差异＞２ｍｍ外，其他测站的差异均＜２ｍｍ。

经统计，第２、３时段水平方向解算精度的中位值为５ ４ｍｍ，第４时段为６ ７ｍｍ，第１时段为２ ２ｍｍ。

结果表明，观测时长缩短会降低解算中的内符合精度，虽然表１中展现的结果差异大部分＜２ｍｍ，但第４时段解算的内符合精度约为６ ７ｍｍ，同时由于在震后５ｈ受到诸多余震的影响，其结果也可能存在偏差。

我们认为将位移量值小于第３时段内符合精度值的结果作为参考值是可行的，但作为最终结果并不可靠。

图３中的蓝色箭头展示了５ｈ静态解算获得的位移场，结果表现出明显的左旋走滑特性，这与该地区的地质构造背景运动性质一致。

在距震中２００～３００ｋｍ的ＱＨＤＬ、ＱＨＭＱ站发生了约１０ｍｍ位移，在ＱＨＢＭ、ＱＨＹＳ发生了约５ｍｍ的位移，超过该范围的测站的位移＜５ｍｍ。

综上所述，分析认为采用震后数ｈ获取的双差静态结果作为快速产出的同震位移场资料在近场基本可靠，在远场位移量值较小的测站位移可作为初步的参考值。

２ ３　多天数据的静态解算方法结果
为了获取更为准确的同震形变结果，本文对震前４ｄ及震后１ｄ的数据进行了分析处理，在处理前去除了地震当天数据文件中震后时段（ＵＴＣ１８：３０—２４：００）的数据。

卫星星历仍然采用ＩＧＳ提供的快速星历（ｉｇｒ），同时添加高阶电离层改正。

经解算统计获取的多天单日解基线ＮＲＭＳ＜０ ２５，约束平差后的测站坐标解算精度为１～２ｍｍ，获取的ＧＮＳＳ同震位移结果如图３
７
７０１①
ＮｉｓｃｈａｎＴ，２０１６，ＧＦＺＲＮＸ ＲＩＮＥＸＧＮＳＳＤａｔａＣｏｎｖｅｒｓｉｏｎａｎｄＭａｎｉｐｕｌａｔｉｏｎＴｏｏｌｂｏｘ，ＧＦＺＤａｔａＳｅｒｖｉｃｅｓ，ｈｔｔｐ：∥ｄｘ．ｄｏｉ．
ｏｒｇ／１０．５８８０／ＧＦＺ．１．１．２０１６．００２。

地　震　地　质４３卷
表１　利用不同时间段的数据获取的坐标值差异（以ＵＴＣ００：００—１８：００为基准）（单位：ｍｍ）
Ｔａｂｌｅ１　Ｔｈｅｄｉｆｆｅｒｅｎｃｅｏｆｃｏｏｒｄｉｎａｔｅｖａｌｕｅｓｏｂｔａｉｎｅｄｆｒｏｍｄａｔａｏｆｄｉｆｆｅｒｅｎｔｔｉｍｅｐｅｒｉｏｄｓ
（ＵＴＣ００：００—１８：００ａｓｔｈｅｂｅｎｃｈｍａｒｋ）（ｕｎｉｔ：ｍｍ）
测站名称ＵＴＣ００：００—０９：００ＵＴＣ０９：００—１８：００
ＮＥＮＥＤＬＨＡ－０．１０．１０．５０．０ＧＳＪＴ０．３－０．１－０．３－０．３ＧＳＬＺ－０．７－１．２１．４０．３ＧＳＭＡ－０．１－０．６－０．５１．０ＧＳＭＸ１．０－０．２－０．４０．８ＱＨＢＭ－１．８－１．１１．４１．０ＱＨＤＬ－０．６０．３０．３０．３ＱＨＧＥ０．１－１．６－０．２０．４ＱＨＭＤ－１．２２．２２．２－０．８ＱＨＭＥ－０．７－０．５１．１０．４ＱＨＭＱ０．６０．１－０．５０．７ＱＨＹＳ４．６－０．４２．１－１．６ＳＣＧＺ０．４－１．９－１．５２．５ＳＣＬＨ１．９－１．７－０．２０．７ＳＣＭＸ０．７０．７－０．８－１．０ＳＣＳＰ０．４－１．０－０．４１．２ＳＣＸＪ－０．２０．７０．２０．０ＸＮＩＮ－１．０－０．７０．８０．５ＸＺＣＤ
２．３
－０．５
－２．５
－０．２
中的绿色箭头所示。

图３结果显示，多天静态解算与震后５ｈ静态解算获取的同震位移结果基本保持一致，从表２发现两者的数值相差１～４ｍｍ。

从数值和地理分布分析，近场测站的ＥＷ向和向ＳＮ向位移都呈现继续在原有方向上增加的趋势，且离震中越近这种增量越大，如ＱＨＭＤ站增加了约４ｍｍ位移、ＱＨＭＱ增加了约１ ５ｍｍ位移。

造成这种现象主要是由于在断层面瞬间发生破裂后断层仍将保持快速滑动而使得形变累积增大，这种现象也是早期震后形变的重要特征（Ｌｉｕ Ｚｅｎｇｅｔａｌ．，２０２０；Ｌｉｅｔａｌ．，２０２１）。

当我们使用震后１ｄ的单日解作为震后测站坐标时，其代表震后１ｄ的平均坐标。

因此，如果考虑震后的持续变形过程，其坐标值必将大于震后５ｈ的静态解结果，因此存在这种差值是合理的，在分析ｄｍ级的同震位移时可以忽略。

但如果作为震后形变分析，表明震后十几ｈ内近场发生了显著变形，例如ＱＨＭＤ站出现了Ｗ向３ ７ｍｍ、Ｎ向３ ３ｍｍ的位移，因此该地区后续的震情仍需持续关注。

除此之外，较远区域测站的多天静态结果与５ｈ静态解结果相比量值更小。

造成这种现象的原因主要是由于震后５ｈ数据受到观测时长及余震的影响，其解算结果的内符合精度较低，多天静态解算的精度能够提高１～２ｍｍ，在位移量值较小的情况下这种改正非常显著。

另外，
８７０１
５期李经纬等：青海玛多７．４级地震ＧＮＳＳ同震水平位移的快速获取
图３　青海玛多７ ４级地震ＧＮＳＳ震后５ｈ解和多天单日解获取的同震位移结果
Ｆｉｇ．３　ＣｏｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔｏｂｔａｉｎｅｄｒｅｓｐｅｃｔｉｖｅｌｙｆｒｏｍｔｈｅＧＮＳＳｓｔａｔｉｏｎ ｓｆｉｖｅ ｈｏｕｒａｎｄｍｕｌｔｉ ｄａｙｓｏｌｕｔｉｏｎｓ
ｏｆｔｈｅＭ
Ｓ
７ ４ＭａｄｕｏｅａｒｔｈｑｕａｋｅｉｎＱｉｎｇｈａｉＰｒｏｖｉｎｃｅ．
箭头表示测站的同震位移，红色和蓝色箭头表示震后５ｈ静态解结果，金色和绿色箭头表示多天单日解结果
表２　利用震后５ｈ数据与多天数据通过静态解方法获取的同震位移量及差值（单位：ｍｍ）Ｔａｂｌｅ２　Ｃｏｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔａｎｄｄｉｆｆｅｒｅｎｃｅｏｂｔａｉｎｅｄｂｙｐｏｓｔｓｅｉｓｍｉｃ５ ｈｏｕｒ
ａｎｄｍｕｌｔｉ ｄａｙｓｔａｔｉｃｓｏｌｕｔｉｏｎｓ，ｒｅｓｐｅｃｔｉｖｅｌｙ（ｕｎｉｔ：ｍｍ）
测站名称
震后５ｈ多天差值
ＮＥＮＥＮＥ
ＤＬＨＡ２．９－５．７３．２－４．６－０．３－１．１ＧＳＪＴ－２．３－２．７－０．８－１．９－１．５－０．８ＧＳＬＺ－３．３－３．２－１．２－２．３－２．０－０．９ＧＳＭＡ－４．１３．０－２．１－１．０－２．０４．０ＧＳＭＸ－１．２－２．５－０．６－０．９－０．６－１．６ＱＨＢＭ－５．８４．５－６．４４．３０．７０．１ＱＨＤＬ１０．７－５．５１１．１－７．２－０．４１．８ＱＨＧＥ０．６１．８２．５－２．６－１．９４．５ＱＨＭＤ８５．８－２３５．３８９．１－２３９．０－３．３３．７ＱＨＭＥ－１．９－２．１－３．４－０．１１．５－１．９ＱＨＭＱ－６．５－１０．５－７．６－１１．３１．２０．８ＱＨＹＳ６．１５．１６．２５．９０．０－０．８ＳＣＧＺ－６．７２．３－４．２－０．３－２．６２．６ＳＣＬＨ－３．８２．９－４．２１．００．４１．９ＳＣＭＸ－２．６－２．７－１．００．１－１．５－２．８ＳＣＳＰ－２．３－２．３－０．２－０．３－２．１－２．０ＳＣＸＪ－４．７－０．４－２．４１．３－２．３－１．７ＸＮＩＮ－１．６－３．０－３．２－４．０１．６１．０ＸＺＣＤ－０．６－１．０－１．０－０．３０．４－０．７
９７
０１
地　震　地　质４３卷
青海玉树测站（ＱＨＹＳ）的运动速率呈现与其他测站运动不一致的现象，该测站的解算精度较其他测站低，具体还需后续与测站环境等因素结合进行讨论。

２ ４　同震形变特征
根据本文所获得的ＧＮＳＳ位移场结果，结合华俊等（２０２１）的ＩｎＳＡＲ结果、王未来等（２０２１）的余震精定位结果及李智敏等（２０２１）的野外考察结果，综合分析认为此次同震地表形变的影响范围主要集中在断裂带两侧９７°～１０１°Ｅ的范围内。

距离震中约４０ｋｍ的ＱＨＭＤ测站的Ｗ向位移约为２４０ｍｍ，Ｎ向位移约为９０ｍｍ，运动方向与断层的走向保持一致。

在此方向延伸处的ＱＨＧＥ站（格尔木）存在同向的微小位移。

震中北侧的ＱＨＤＬ和ＤＬＨＡ站呈现ＮＷ向的位移，但与ＱＨＭＤ站的运动方向相比明显更偏于Ｎ向。

位于青海玛曲的ＱＨＭＱ测站的Ｓ向运动约为８ｍｍ，Ｗ向运动约为１１ｍｍ，Ｓ向运动的比例相对于同区域ＱＨＭＤ、ＱＨＤＬ站的位移方向存在明显差异，造成这种现象可能是ＱＨＭＱ站位于破碎带上。

位于断层另一侧的ＱＨＢＭ站（班玛），其Ｓ向位移约为６ｍｍ，Ｅ向位移约为４ｍｍ，呈现ＳＥ向的运动位移。

同时甘孜、炉霍等地区受到较小的波动，ＳＣＧＺ、ＳＣＬＨ和ＳＣＸＪ测站的Ｅ向位移不明显，呈现Ｓ向位移特征，位移量约为４ｍｍ。

整体而言，断裂北侧呈现ＮＷ向位移，南侧呈现ＳＥ向位移，在远离断层处垂直于断层方向存在背离断层面的运动趋势。

３　结论
本文对通过１Ｈｚ高频数据动态差分、震后５ｈ静态解算及多天静态解算３种方法获取的青海玛多７
４级地震地表同震位移开展了对比分析，得到以下结论：（１）ＴＲＡＣＫ模块具有快速解算的能力，其高频动态差分结果在近场位移量较大的站点上同震位移显著，但在远场受测站距离所限解算精度为ｃｍ级，较难显示＜１ｃｍ的形变量。

虽然ＴＲＡＣＫ的结果不精确，但可以作为快速获取近场测站位移量及初步判定断层运动特征的手段；
（２）震后５ｈ静态解算获得的同震位移与多天静态解算结果相比，两者的位移大小、方向以及影响范围都保持较强的一致性。

多天静态解的同震结果与震后５ｈ静态解的结果相差１～４ｍｍ。

因此，考虑震后５ｈ数据的时效性，认为震后５ｈ静态解算方法的结果作为同震位移快速结果是可信的。

（３）综合多方面结果开展分析，可初步判断此次地震是以左旋走滑为主的地震。

目前测站中离震中最近的约４０ｋｍ的ＱＨＭＤ站的最大位移约为２４０ｍｍ。

地震的影响范围波及了震中玛多县周边区域，包括北部西宁、都兰、德令哈地区及南部斑马、甘孜地区。

由于多天静态解包含了震后１ｄ的形变，其形变方向与震后５ｈ形变保持方向一致，且量值增加，这种现象与震后余滑较吻合，因此认为近场的震后形变现象显著，后期仍需持续关注。

致谢　中国地震局台网中心提供了ＧＮＳＳ高频数据，在此表示衷心感谢！
参
考
文
献
邓起东，冉勇康，杨晓平，等．２００７．中国活动构造图［ＣＭ］．北京：地震出版社．
ＤＥＮＧＱｉ ｄｏｎｇ，ＲＡＮＹｏｎｇ ｋａｎｇ，ＹＡＮＧＸｉａｏ ｐｉｎｇ，ｅｔａｌ．２００７．ＡｃｔｉｖｅＴｅｃｔｏｎｉｃＭａｐｏｆＣｈｉｎａ［ＣＭ］．Ｓｅｉｓｍｏｌｏｇｉｃａｌ
０８０１
５期李经纬等：青海玛多７．４级地震ＧＮＳＳ同震水平位移的快速获取Ｐｒｅｓｓ，Ｂｅｉｊｉｎｇ（ｉｎＣｈｉｎｅｓｅ）．
范士杰，刘焱雄，张健，等．２０１３．ＰＰＰ与ＧＡＭＩＴ／ＴＲＡＣＫ在地震监测中的应用［Ｊ］．测绘科学，３８（２）：１８４—
１８６．
ＦＡＮＳｈｉ ｊｉｅ，ＬＩＵＹａｎ ｘｉｏｎｇ，ＺＨＡＮＧＪｉａｎ，ｅｔａｌ．２０１３．ＡｐｐｌｉｃａｔｉｏｎｏｆＰＰＰａｎｄＧＡＭＩＴ／ＴＲＡＣＫｍｏｄｕｌｅｉｎ
ｅａｒｔｈｑｕａｋｅｍｏｎｉｔｏｒｉｎｇ［Ｊ］．ＳｃｉｅｎｃｅｏｆＳｕｒｖｅｙｉｎｇａｎｄＭａｐｐｉｎｇ，３８（２）：１８４—１８６（ｉｎＣｈｉｎｅｓｅ）．高翔，邓起东．２０１３．巴颜喀喇断块边界断裂强震活动分析［Ｊ］．地质学报，８７（１）：９—１９．
ＧＡＯＸｉａｎｇ
，ＤＥＮＧＱｉ ｄｏｎｇ．２０１３．ＡｎａｌｙｓｉｓｏｆｌａｒｇｅｅａｒｔｈｑｕａｋｅｓｉｎｂｏｕｎｄａｒｙｆａｕｌｔｓａｒｏｕｎｄＢａｙａｎｋａｌａｆａｕｌｔ ｂｌｏｃｋ［Ｊ］．ＡｃｔａＧｅｏｌｏｇｉｃａＳｉｎｉｃａ，８７（１）：９—１９（ｉｎＣｈｉｎｅｓｅ）．
华俊，赵德政，单新建，等．２０２１．２０２１年青海玛多ＭＷ７ ３地震ＩｎＳＡＲ的同震形变场、断层滑动分布及其对周
边区域的应力扰动［Ｊ］．地震地质，４３（３）：６７７—６９１．ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５３－４９６７．２０２１．０３．０１３．ＨＵＡＪｕｎ，ＺＨＡＯＤｅ ｚｈｅｎｇ，ＳＨＡＮＸｉｎ ｊｉａｎ，ｅｔａｌ．２０２１．Ｃｏｓｅｉｓｍｉｃｄｅｆｏｒｍａｔｉｏｎｆｉｅｌｄ，ｓｌｉｐｄｉｓｔｒｉｂｕｔｉｏｎａｎｄＣｏｕｌｏｍｂ
ｓｔｒｅｓｓｄｉｓｔｕｒｂａｎｃｅｏｆｔｈｅ２０２１ＭＷ７ ３ＭａｄｕｏｅａｒｔｈｑｕａｋｅｕｓｉｎｇＳｅｎｔｉｎｅｌ－１ＩｎＳＡＲｏｂｓｅｒｖａｔｉｏｎｓ［Ｊ］．ＳｅｉｓｍｏｌｏｇｙａｎｄＧｅｏｌｏｇｙ，４３（３）：６７７—６９１（ｉｎＣｈｉｎｅｓｅ）．
李智敏，李文巧，李涛，等．２０２１．２０２１年５月２２日青海玛多ＭＳ７ ４地震的发震构造和地表破裂初步调查
［
Ｊ］．地震地质，４３（３）：７２２—７３７．ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５３－４９６７．２０２１．０３．１６．ＬＩＺｈｉ ｍｉｎ，ＬＩＷｅｎ ｑｉａｏ，ＬＩＴａｏ，ｅｔａｌ．２０２１．ＳｅｉｓｍｏｇｅｎｉｃｆａｕｌｔａｎｄｃｏｓｅｉｓｍｉｃｓｕｒｆａｃｅｄｅｆｏｒｍａｔｉｏｎｏｆｔｈｅＭａｄｕｏ
ＭＳ７ ４ｅａｒｔｈｑｕａｋｅｉｎＱｉｎｇｈａｉ，Ｃｈｉｎａ：Ａｑｕｉｃｋｒｅｐｏｒｔ［Ｊ］．ＳｅｉｓｍｏｌｏｇｙａｎｄＧｅｏｌｏｇｙ，４３（３）：７２２—７３７（ｉｎＣｈｉｎｅｓｅ）．
梁明剑，杨耀，杜方，等．２０２０．青海达日断裂中段晚第四纪活动性与１９４７年Ｍ７
地震地表破裂带再研究［Ｊ］．地震地质，４２（３）：７０３—７１４．ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５３－４９６７．２０２０．０３．０１１．
ＬＩＡＮＧＭｉｎｇ ｊｉａｎ
，ＹＡＮＧＹａｏ，ＤＵＦａｎｇ，ｅｔａｌ．２０２０．ＬａｔｅＱｕａｔｅｒｎａｒｙａｃｔｉｖｉｔｙｏｆｔｈｅｃｅｎｔｒａｌｓｅｇｍｅｎｔｏｆｔｈｅＤａｒｉＦａｕｌｔａｎｄｒｅｓｔｕｄｙｏｆｔｈｅｓｕｒｆａｃｅｒｕｐｔｕｒｅｚｏｎｅｏｆｔｈｅ１９４７Ｍ７ Ｄａｒｉｅａｒｔｈｑｕａｋｅ．ＱｉｎｇｈａｉＰｒｏｖｉｎｃｅ［Ｊ］．ＳｅｉｓｍｏｌｏｇｙａｎｄＧｅｏｌｏｇｙ，４２（３）：７０３—７１４（ｉｎＣｈｉｎｅｓｅ）．
潘家伟，白明坤，李超，等．２０２１ ２０２１年５月２２日青海玛多ＭＳ７ ４地震地表破裂带及发震构造［Ｊ］．地质学
报，９５（６）：１６５５—１６７０．
ＰＡＮＪｉａ ｗｅｉ，ＢＡＩＭｉｎｇ ｋｕｎ，ＬＩＣｈａｏ，ｅｔａｌ．２０２１．Ｃｏｓｅｉｓｍｉｃｓｕｒｆａｃｅｒｕｐｔｕｒｅａｎｄｓｅｉｓｍｏｇｅｎｉｃｓｔｒｕｃｔｕｒｅｏｆｔｈｅ２０２１
－０５－２２Ｍａｄｕｏ（Ｑｉｎｇｈａｉ）ＭＳ７ ４ｅａｒｔｈｑｕａｋｅ［Ｊ］．ＡｃｔａＧｅｏｌｏｇｉｃａＳｉｎｉｃａ，９５（６）：１６５５—１６７０（ｉｎＣｈｉｎｅｓｅ）．王未来，房立华，吴建平，等．２０２１．２０２１年青海玛多ＭＳ７ ４地震序列精定位研究［Ｊ］．中国科学（Ｄ辑），
５１（７）：１—１０．
ＷＡＮＧＷｅｉ ｌａｉ，ＦＡＮＧＬｉ ｈｕａ，ＷＵＪｉａｎ ｐｉｎｇ，ｅｔａｌ．２０２１．Ａｆｔｅｒｓｈｏｃｋｓｅｑｕｅｎｃｅｒｅｌｏｃａｔｉｏｎｏｆｔｈｅ２０２１ＭＳ７ ４Ｍａｄｕｏ
ｅａｒｔｈｑｕａｋｅ，Ｑｉｎｇｈａｉ，Ｃｈｉｎａ［Ｊ］．ＳｃｉｅｎｃｅｉｎＣｈｉｎａ（ＳｅｒＤ），５１（７）：１—１０（ｉｎＣｈｉｎｅｓｅ）．徐锡伟，韩竹军，杨晓平，等．２０１６．中国及邻近地区地震构造图［ＣＭ］．北京：地震出版社．
ＸＵＸｉ ｗｅｉ
，ＨＡＮＺｈｕ ｊｕｎ，ＹＡＮＧＸｉａｏ ｐｉｎｇ，ｅｔａｌ．２０１６．ＳｅｉｓｍｏｔｅｃｔｏｎｉｃＭａｐｏｆＣｈｉｎａａｎｄＡｄｊａｃｅｎｔＡｒｅａｓ［ＣＭ］．ＳｅｉｓｍｏｌｏｇｉｃａｌＰｒｅｓｓ，Ｂｅｉｊｉｎｇ（ｉｎＣｈｉｎｅｓｅ）．
殷海涛，肖根如，张磊，等．２０１２．ＴＲＡＣＫ高频ＧＰＳ定位中震时参考站的选取方法［Ｊ］．大地测量与地球动力
学，３
２（４）：１５—１９．ＹＩＮＨａｉ ｔａｏ
，ＸＩＡＯＧｅｎ ｒｕ，ＺＨＡＮＧＬｅｉ，ｅｔａｌ．２０１２．ＳｔｕｄｙｏｎｍｅｔｈｏｄｆｏｒｓｅｌｅｃｔｉｎｇｒｅｆｅｒｅｎｃｅｓｔａｔｉｏｎｉｎｈｉｇｈｒａｔｅＧＰＳｐｏｓｉｔｉｏｎｉｎｇｕｓｉｎｇｔｒａｃｋｄｕｒｉｎｇｅａｒｔｈｑｕａｋｅ［Ｊ］．ＪｏｕｒｎａｌｏｆＧｅｏｄｅｓｙａｎＧｅｏｄｙｎａｍｉｃｓ，３２（４）：１５—１９（ｉｎＣｈｉｎｅｓｅ
）．张培震，邓起东，张国民，等．２００３．中国大陆的强震活动与活动地块［Ｊ］．中国科学（Ｄ辑），３３（Ｓ１）：１２—２０．ＺＨＡＮＧＰｅｉ ｚｈｅｎ，ＤＥＮＧＱｉ ｄｏｎｇ，ＺＨＡＮＧＧｕｏ ｍｉｎ，ｅｔａｌ．２００３．Ｓｔｒｏｎｇｅａｒｔｈｑｕａｋｅａｃｔｉｖｉｔｙａｎｄａｃｔｉｖｅｂｌｏｃｋｓｉｎ
ｍａｉｎｌａｎｄＣｈｉｎａ［
Ｊ］．ＳｃｉｅｎｃｅｉｎＣｈｉｎａ（ＳｅｒＤ），３３（Ｓ１）：１２—２０（ｉｎＣｈｉｎｅｓｅ）．１
８０１
地　震　地　质４３卷
中国地震局震害防御司．１９９９．中国近代地震目录（公元１９１２—１９９０年ＭＳ≥４ ７）［Ｚ］．北京：中国科学技术出
版社．
ＤｅｐａｒｔｍｅｎｔｏｆＥａｒｔｈｑｕａｋｅＤｉｓａｓｔｅｒＰｒｅｖｅｎｔｉｏｎ，ＣｈｉｎａＥａｒｔｈｑｕａｋｅＡｄｍｉｎｉｓｔｒａｔｉｏｎ．１９９９．ＣａｔａｌｏｇｕｅｏｆＲｅｃｅｎｔＣｈｉｎｅｓｅ
Ｅａｒｔｈｑｕａｋｅｓ（１９１２—１９９０ＡＤ，ＭＳ≥４ ７）［Ｚ］．ＣｈｉｎａＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙＰｒｅｓｓ，Ｂｅｉｊｉｎｇ（ｉｎＣｈｉｎｅｓｅ）．周命端，覃发超，杨玉忠，等．２０１５．利用ＣＧＯ分析处理ＧＰＳ工程控制网的研究与实践［Ｊ］．测绘通报，（Ｓ１）：
８４—８７．
ＺＨＯＵＭｉｎｇ ｄｕａｎ，ＱＩＮＦａ ｃｈａｏ，ＹＡＮＧＹｕ ｚｈｏｎｇ，ｅｔａｌ．２０１５．ＲｅｓｅａｒｃｈｏｎｕｓｉｎｇＣＧＯｔｏａｎａｌｙｚｅａｎｄｐｒｏｃｅｓｓｔｈｅ
ＧＰＳｅｎｇｉｎｅｅｒｉｎｇｃｏｎｔｒｏｌｎｅｔｗｏｒｋ［Ｊ］．ＢｕｌｌｅｔｉｎｏｆＳｕｒｖｅｙｉｎｇａｎｄＭａｐｐｉｎｇ，（Ｓ１）：８４—８７（ｉｎＣｈｉｎｅｓｅ）．ＢａｎｅｒｊｅｅＰ，ＰｏｌｌｉｔｚＦ，ＢüｒｇｍａｎｎＲ．２００５．ＴｈｅｓｉｚｅａｎｄｄｕｒａｔｉｏｎｏｆｔｈｅＳｕｍａｔｒａ Ａｎｄａｍａｎｅａｒｔｈｑｕａｋｅｆｒｏｍｆａｒ ｆｉｅｌｄ
ｓｔａｔｉｃｏｆｆｓｅｔｓ
［Ｊ］．Ｓｃｉｅｎｃｅ，３０８（５７２９）：１７６９—１７７２．ＧａｎＷＪ，ＺｈａｎｇＰＺ，ＳｈｅｎＺＫ，ｅｔａｌ．２００７．Ｐｒｅｓｅｎｔ ｄａｙｃｒｕｓｔａｌｍｏｔｉｏｎｗｉｔｈｉｎｔｈｅＴｉｂｅｔａｎｐｌａｔｅａｕｉｎｆｅｒｒｅｄｆｒｏｍＧＰＳ
ｍｅａｓｕｒｅｍｅｎｔｓ
［Ｊ］．ＪｏｕｒｎａｌｏｆＧｅｏｐｈｙｓｉｃａｌＲｅｓｅａｒｃｈ：ＳｏｌｉｄＥａｒｔｈ，１１２（Ｂ８）：Ｂ０８４１６．ＧｕｏＢＦ，ＺｈａｎｇＸＨ，ＲｅｎＸ，ｅｔａｌ．２０１５．Ｈｉｇｈ ｐｒｅｃｉｓｉｏｎｃｏｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔｅｓｔｉｍａｔｉｏｎｗｉｔｈａｓｉｎｇｌｅ ｆｒｅｑｕｅｎｃｙ
ＧＰＳｒｅｃｅｉｖｅｒ［Ｊ］．ＧｅｏｐｈｙｓｉｃａｌＪｏｕｒｎａｌＩｎｔｅｒｎａｔｉｏｎａｌ，２０２（１）：６１２—６２３．
ＬｉＪＷ，Ｚ
ｈａｎＷ，ＧｕｏＢＦ，ｅｔａｌ．２０２１．ＣｏｍｂｉｎａｔｉｏｎｏｆｔｈｅＬｅｖｅｎｂｅｒｇ ＭａｒｑｕａｒｄｔａｎｄｄｉｆｆｅｒｅｎｔｉａｌｅｖｏｌｕｔｉｏｎａｌｇｏｒｉｔｈｍｓｆｏｒｔｈｅｆｉｔｔｉｎｇｏｆｐｏｓｔｓｅｉｓｍｉｃＧＰＳｔｉｍｅｓｅｒｉｅｓ［Ｊ］．ＡｃｔａＧｅｏｐｈｙｓｉｃａ，６９（２）：４０５—４１４．
Ｌｉｕ ＺｅｎｇＪ
，ＺｈａｎｇＺ，ＲｏｌｌｉｎｓＣ，ｅｔａｌ．２０２０．Ｐｏｓｔｓｅｉｓｍｉｃｄｅｆｏｒｍａｔｉｏｎｆｏｌｌｏｗｉｎｇｔｈｅ２０１５ＭＷ７ ８Ｇｏｒｋｈａ（Ｎｅｐａｌ）ｅａｒｔｈｑｕａｋｅ：ＮｅｗＧＰＳｄａｔａ，ｋｉｎｅｍａｔｉｃａｎｄｄｙｎａｍｉｃｍｏｄｅｌｓ，ａｎｄｔｈｅｒｏｌｅｓｏｆａｆｔｅｒｓｌｉｐａｎｄｖｉｓｃｏｅｌａｓｔｉｃｒｅｌａｘａｔｉｏｎ［Ｊ］．ＪｏｕｒｎａｌｏｆＧｅｏｐｈｙｓｉｃａｌＲｅｓｅａｒｃｈ：ＳｏｌｉｄＥａｒｔｈ，１２５（９）：ｅ２０２０ＪＢ０１９８５２．
ＰｓｉｍｏｕｌｉｓＰ，ＨｏｕｌｉéＮ，ＭｉｃｈｅｌＣ，ｅｔａｌ．２０１６．Ａｎａｌｙｓｉｓｏｆ１ ＨｚＧＰＳｄａｔａｆｏｒｔｈｅｅｓｔｉｍａｔｉｏｎｏｆｌｏｎｇ ｐｅｒｉｏｄｓｕｒｆａｃｅ
ｍｏｔｉｏｎｏｆＴｏｈｏｋｕ ＯｋｉＭＷ９ ０２０１１ｅａｒｔｈｑｕａｋｅ［Ｃ］．３ｒｄＪｏｉｎｔＩｎｔｅｒｎａｔｉｏｎａｌＳｙｍｐｏｓｉｕｍｏｎＤｅｆｏｒｍａｔｉｏｎＭｏｎｉｔｏｒｉｎｇ
，Ｖｉｅｎｎａ，Ａｕｓｔｒｉａ．ＳｕＫ，ＪｉｎＳＧ．２０２０．Ｒｅａｌ ｔｉｍｅｓｅｉｓｍｉｃｗａｖｅｆｏｒｍｓｅｓｔｉｍａｔｉｏｎｏｆｔｈｅ２０１９ＭＷ＝６ ４ａｎｄＭＷ＝７ １Ｃａｌｉｆｏｒｎｉａ
ｅａｒｔｈｑｕａｋｅｓｗｉｔｈｈｉｇｈ ｒａｔｅｍｕｌｔｉ ＧＮＳＳｏｂｓｅｒｖａｔｉｏｎｓ
［Ｊ］．ＩＥＥＥＡｃｃｅｓｓ，８（１）：８５４１１—８５４２０．ＴｉｒｙａｋｉｏｇｌｕＩ，ＹｉｇｉｔＣＯ，ＹａｖａｓｏｇｌｕＨ，ｅｔａｌ．２０１７．Ｔｈｅｄｅｔｅｒｍｉｎａｔｉｏｎｏｆｉｎｔｅｒｓｅｉｓｍｉｃ，ｃｏｓｅｉｓｍｉｃａｎｄｐｏｓｔｓｅｉｓｍｉｃ
ｄｅｆｏｒｍａｔｉｏｎｓｃａｕｓｅｄｂｙｔｈｅＧ ｋ ｅａｄａ Ｓａｍｏｔｈｒａｋｉｅａｒｔｈｑｕａｋｅ（２０１４，ＭＷ６ ９）ｂａｓｅｄｏｎＧＮＳＳｄａｔａ［Ｊ］．ＪｏｕｒｎａｌｏｆＡｆｒｉｃａｎＥａｒｔｈＳｃｉｅｎｃｅｓ，１３３：８６—９４．
ＷａｎｇＭ，Ｓ
ｈｅｎＺＫ．２０２０．Ｐｒｅｓｅｎｔ ｄａｙｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｏｆｃｏｎｔｉｎｅｎｔａｌＣｈｉｎａｄｅｒｉｖｅｄｆｒｏｍＧＰＳａｎｄｉｔｓｔｅｃｔｏｎｉｃｉｍｐｌｉｃａｔｉｏｎｓ［Ｊ］．ＪｏｕｒｎａｌｏｆＧｅｏｐｈｙｓｉｃａｌＲｅｓｅａｒｃｈ：ＳｏｌｉｄＥａｒｔｈ，１２５（２）：ｅ２０１９ＪＢ０１８７７４．
ＷｕＹＱ，ＪｉａｎｇＺＳ，ＬｉａｎｇＨＢ，ｅｔａｌ．２０１６．Ｃｏｓｅｉｓｍｉｃｄｅｆｏｒｍａｔｉｏｎｓｏｆｔｈｅ２０１５ＭＷ７ ８Ｇｏｒｋｈａｅａｒｔｈｑｕａｋｅａｎｄ
ｉｎｔｅｒｓｅｉｓｍｉｃｓｔｒａｉｎａｃｃｕｍｕｌａｔｉｏｎｉｎｔｈｅＨｉｍａｌａｙａｎｔｅｃｔｏｎｉｃｂｅｌｔａｎｄＴｉｂｅｔａｎｐｌａｔｅａｕ［Ｊ］．Ｔｅｃｔｏｎｏｐｈｙｓｉｃｓ，６７０：１４４—１５４．
ＺｈｕＹＧ
，ＤｉａｏＦＱ，ＦｕＹＣ，ｅｔａｌ．２０２１．Ｓｌｉｐｒａｔｅｏｆｔｈｅｓｅｉｓｍｏｇｅｎｉｃｆａｕｌｔｏｆｔｈｅ２０２１ＭａｄｕｏｅａｒｔｈｑｕａｋｅｉｎｗｅｓｔｅｒｎＣｈｉｎａｉｎｆｅｒｒｅｄｆｒｏｍＧＰＳｏｂｓｅｒｖａｔｉｏｎｓ［Ｊ］．ＳｃｉｅｎｃｅＣｈｉｎａＥａｒｔｈＳｃｉｅｎｃｅｓ，６４（８）：１３６３—１３７０．
２８０１
５期李经纬等：青海玛多７．４级地震ＧＮＳＳ同震水平位移的快速获取ＲＥＳＥＡＲＣＨＯＮＦＡＳＴＡＣＱＵＩＳＩＴＩＯＮＯＦＧＮＳＳＣＯＳＥＩＳＭＩＣ
ＨＯＲＩＺＯＮＴＡＬＤＩＳＰＬＡＣＥＭＥＮＴＯＦＭＡＤＵＯＭＳ７．４
ＥＡＲＴＨＱＵＡＫＥＩＮＱＩＮＧＨＡＩＰＲＯＶＩＮＣＥ
ＬＩＪｉｎｇ ｗｅｉ１，２）　ＣＨＥＮＣｈａｎｇ ｙｕｎ１）　ＺＨＡＮＷｅｉ１）　ＷＵＹａｎ ｑｉａｎｇ
１）１）ＴｈｅＦｉｒｓｔＭｏｎｉｔｏｒｉｎｇａｎｄＡｐｐｌｉｃａｔｉｏｎＣｅｎｔｅｒ，ＣｈｉｎａＥａｒｔｈｑｕａｋｅＡｄｍｉｎｉｓｔｒａｔｉｏｎ，Ｔｉａｎｊｉｎ　３００１８０，Ｃｈｉｎａ
２）ＩｎｎｏｖａｔｉｏｎＡｃａｄｅｍｙｆｏｒＰｒｅｃｉｓｉｏｎＭｅａｓｕｒｅｍｅｎｔＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ，
ＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ，Ｗｕｈａｎ　４３００７１，Ｃｈｉｎａ
Ａｂｓｔｒａｃｔ
ＴｈｅＭａｙ２１，２０２１，ＭａｄｕｏＭＳ７ ４ｅａｒｔｈｑｕａｋｅｉｎＱｉｎｇｈａｉＰｒｏｖｉｎｃｅｃａｕｓｅｄｓｅｒｉｏｕｓｄｉｓａｓｔｅｒｓｉｎ
ＭａｄｕｏＣｏｕｎｔｙａｎｄｉｔｓｓｕｒｒｏｕｎｄｉｎｇａｒｅａｓ．ＴｈｅＧＮＳＳｃｏ ｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔｆｉｅｌｄｄａｔａｃａｎｐｌａｙａｋｅｙｒｏｌｅｉｎｑｕｉｃｋｌｙｄｅｔｅｒｍｉｎｉｎｇｔｈｅｉｎｆｌｕｅｎｃｅｒａｎｇｅｏｆｔｈｅｅａｒｔｈｑｕａｋｅａｎｄｓｅｒｖｉｎｇｆｏｒｔｈｅｒａｐｉｄ
ｉｎｖｅｓｔｉｇａｔｉｏｎ．Ａｆｔｅｒｔｈｅｅａｒｔｈｑｕａｋｅ
，ｗｅｉｍｍｅｄｉａｔｅｌｙｃｏｌｌｅｃｔｅｄｔｈｅｄａｔａｏｆ１８ＧＮＳＳｓｔａｔｉｏｎｓｓｕｒ ｒｏｕｎｄｉｎｇｔｈｅｅｐｉｃｅｎｔｅｒ，ｉｎｃｌｕｄｉｎｇ７ｓｔａｔｉｏｎｓｔｈａｔｒｅｃｏｒｄｅｄ１Ｈｚｈｉｇｈ ｆｒｅｑｕｅｎｃｙｏｂｓｅｒｖａｔｉｏｎｄａｔａ．ＶａｒｉｏｕｓｄａｔａｗｅｒｅｕｓｅｄｔｏｒａｐｉｄｌｙｏｂｔａｉｎｔｈｅＧＮＳＳｃｏ ｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔｓ，ｓｕｃｈａｓ，ｔｈｅ１５ ｍｉｎｕｔｅ
ｈｉｇｈ ｆｒｅｑｕｅｎｃｙｄａｔａ
，５ｈｏｕｒｓａｆｔｅｒｅａｒｔｈｑｕａｋｅａｎｄｍｕｌｔｉ ｄａｙｄｉｓｐｌａｃｅｍｅｎｔｄａｔａ．Ｉｎｔｈｉｓｐａｐｅｒ，ｗｅｕｓｅｄｔｈｒｅｅｍｅｔｈｏｄｓｔｏｏｂｔａｉｎｔｈｅｃｏ ｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔ
，ｉｎｃｌｕｄｉｎｇｔｈｅｄｙｎａｍｉｃｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄｆｏｒ１ＨｚｆｒｅｑｕｅｎｃｙｄａｔａｂｙＧＡＭＩＴ／ＧＬＯＢＫＴｒａｃｋｍｏｄｕｌｅ，ａｎｄｔｈｅｓｔａｔｉｃｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄｆｏｒｔｈｅｐｏｓｔ ｓｅｉｓｍｉｃ５ ｈｏｕｒｄａｔａａｎｄｆｏｒｔｈｅｐｒｅ ａｎｄｐｏｓｔ ｓｅｉｓｍｉｃｍｕｌｔｉ ｄａｙｄａｔａｂｙＧＡＭＩＴ／ＧＬＯＢＫ．Ｔｈｅ
ｒｅｓｕｌｔｓａｒｅｓｈｏｗｎａｓｆｏｌｌｏｗｓ
：（１）Ｔｈｅｄｙｎａｍｉｃｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄｆｏｒ１ＨｚｆｒｅｑｕｅｎｃｙｄａｔａｂｙＧＡＭＩＴ／ＧＬＯＢＫＴｒａｃｋｍｏｄｕｌｅｈａｓａｂｉｌｉｔｙｔｏｑｕｉｃｋｌｙｐｒｏｃｅｓｓｔｈｅｄａｔａａｎｄａｃｑｕｉｒｅｔｈｅｃｏ ｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔ．Ｗｈｅｎｕｓｉｎｇｔｈｅｈｉｇｈ ｆｒｅｑｕｅｎｃｙｄａｔａｔｏｏｂｔａｉｎｃｏ ｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔｂｙＴｒａｃｋｍｏｄｕｌｅ，ｉｔｉｓｓｕｉｔａｂｌｅｆｏｒｔｈｅｎｅａｒｆｉｅｌｄ
ｓｔａｔｉｏｎｓｗｈｉｃｈｈａｖｅａｌａｒｇｅｖａｌｕｅｏｆｃｏ ｓｅｉｓｍｉｃｄｅｆｏｒｍａｔｉｏｎ．Ｈｏｗｅｖｅｒ
，ｉｎｔｈｅｆａｒｆｉｅｌｄ，ｔｈｅａｃｃｕｒａｃｙｏｆｔｈｅｓｏｌｕｔｉｏｎｉｓａｔｃｍｌｅｖｅｌｒｅｓｔｒｉｃｔｅｄｂｙｔｈｅｄｉｓｔａｎｃｅｏｆｓｔａｔｉｏｎｓ．Ｉｎａｄｄｉｔｉｏｎ
，ｔｈｅｒｅｓｕｌｔｏｆｔｈｅＴｒａｃｋｉｓｉｎｆｌｕｅｎｃｅｄｂｙｔｈｅｓｔａｂｉｌｉｔｙｏｆｒｅｆｅｒｅｎｃｅｓｔａｔｉｏｎ．ＡｌｔｈｏｕｇｈｔｈｅｒｅｓｕｌｔｓｏｂｔａｉｎｅｄｂｙＴｒａｃｋａｒｅｎｏｔａｃｃｕｒａｔｅ，ｉｔｃａｎｂｅｕｓｅｄａｓａｍｅｔｈｏｄｔｏｑｕｉｃｋｌｙｊｕｄｇｅｔｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓａｎｄａｍｏｕｎｔｏｆｃｏｓｅｉｓｍｉｃｓｕｒｆａｃｅｍｏｔｉｏｎ．
（２）Ｃｏｍｐａｒｉｎｇｔｈｅｒｅｓｕｌｔｓｏｂｔａｉｎｅｄｆｒｏｍｔｈｅｐｏｓｔ ｓｅｉｓｍｉｃ５ ｈｏｕｒｄａｔａａｎｄｔｈｅｐｒｅ ａｎｄｐｏｓｔ ｓｅｉｓｍｉｃｍｕｌｔｉ ｄａｙｄａｔａ，ｔｈｅＧＮＳＳｓｔａｔｉｏｎｓ ｄｉｓｐｌａｃｅｍｅｎｔｓｈａｖｅｇｏｏｄｃｏｎｓｉｓｔｅｎｃｙｉｎｔｈｅｍａｇｎｉｔｕｄｅ，ｄｉｒｅｃｔｉｏｎａｎｄｉｎｆｌｕｅｎｃｅｒａｎｇｅ，ｅｓｐｅｃｉａｌｌｙｉｎｔｈｅｎｅａｒｆｉｅｌｄ．Ｔｈｅｄｉｆｆｅｒｅｎｃｅｏｆｔｈｅｒｅｓｕｌｔｓｂｙｔｈｅｔｗｏｍｅｔｈｏｄｓｉｓｆｒｏｍ１ｍｍｔｏ４ｍｍ．ＣｏｎｓｉｄｅｒｉｎｇｔｈｅｐｒｏｃｅｓｓｉｎｇａｃｃｕｒａｃｙｏｆｔｈｅＧＡＭＩＴ／ＧＬＯＢＫ，ｔｈｅｖａｌｕｅ
ｏｆｔｈｅｄｉｆｆｅｒｅｎｃｅｉｓｎｏｔｕｎｒｅａｓｏｎａｂｌｙｈｉｇｈ．Ｗｈｅｎｔｈｅｄｉｓｐｌａｃｅｍｅｎｔｖａｌｕｅｉｓｓｍａｌｌ
，ｉｔｉｓｄｉｆｆｉｃｕｌｔｔｏｏｂｔａｉｎａｃｃｕｒａｔｅｒｅｓｕｌｔｓ．Ｉｎａｄｄｉｔｉｏｎ，ｔｈｅｄｉｒｅｃｔｉｏｎｏｆｔｈｅｐｒｅ ａｎｄｐｏｓｔ ｓｅｉｓｍｉｃｍｕｌｔｉ ｄａｙｒｅｓｕｌｔｉｓｃｏｎｓｉｓｔｅｎｔｗｉｔｈｔｈａｔｆｒｏｍｔｈｅｐｏｓｔ ｓｅｉｓｍｉｃ５ ｈｏｕｒｄａｔａ，ａｎｄｔｈｅｖａｌｕｅｉｎｃｒｅａｓｅｄ．Ｉｆｗｅｒｅｇａｒｄｔｈｅｒｅｓｕｌｔｏｆｔｈｅｐｒｅ ａｎｄｐｏｓｔ ｓｅｉｓｍｉｃｍｕｌｔｉ ｄａｙｄａｔａａｓｔｈｅｒｅｓｕｌｔｏｆｏｎｅｄａｙｄａｔａａｆｔｅｒｔｈｅｅａｒｔｈｑｕａｋｅｗｈｉｃｈｉｓ３８０１
地　震　地　质４３卷ｉｎｃｌｕｄｅｄｉｎｔｈｅｐｏｓｔ ｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔ，ｔｈｉｓｐｈｅｎｏｍｅｎｏｎｃｏｉｎｃｉｄｅｓｗｉｔｈｔｈｅａｆｔｅｒｓｌｉｐｄｅｆｏｒｍａｔｉｏｎ，ａｎｄｔｈｅｄｉｆｆｅｒｅｎｃｅｍａｙｂｅｃａｕｓｅｄｂｙｔｈｅａｆｔｅｒｓｌｉｐ，ｅｓｐｅｃｉａｌｌｙｉｎｔｈｅｎｅａｒｆｉｅｌｄ．Ａｌｔｈｏｕｇｈｔｈｅｄｉｆｆｅｒｅｎｃｅｅｘｉｓｔｓ，ｔａｋｉｎｇｉｎｔｏａｃｃｏｕｎｔｔｈｅｔｉｍｅｌｉｎｅｓｓａｎｄｔｈｅｏｖｅｒａｌｌｃｏｎｓｉｓｔｅｎｃｙ，ｗｅｂｅｌｉｅｖｅｔｈａｔｕｓｉｎｇｔｈｅｐｏｓｔｓｅｉｓｍｉｃ５ ｈｏｕｒｄａｔａｔｏｑｕｉｃｋｌｙｏｂｔａｉｎｔｈｅｃｏ ｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔｉｓｃｒｅｄｉｂｌｅｉｎａｎｅｍｅｒｇｅｎｃｙ．
（３）Ｂａｓｅｄｏｎｔｈｅａｎａｌｙｓｉｓｏｆｖａｒｉｏｕｓｒｅｓｕｌｔｓ，ｉｔｉｓｐｒｅｌｉｍｉｎａｒｉｌｙｊｕｄｇｅｄｔｈａｔｔｈｅＭａｄｕｏｅａｒｔｈｑｕａｋｅ
ｉｓｄｏｍｉｎａｔｅｄｂｙｌｅｆｔ ｈａｎｄｅｄｓｔｒｉｋｅ ｓｌｉｐ．ＴｈｅｍａｘｉｍｕｍｄｉｓｐｌａｃｅｍｅｎｔａｔｔｈｅｓｔａｔｉｏｎＱＨＭＤ
，ｗｈｉｃｈｉｓａｂｏｕｔ４０ｋｍｆｒｏｍｔｈｅｅｐｉｃｅｎｔｅｒ，ｉｓａｂｏｕｔ２４ｃｍｔｏｔｈｅｗｅｓｔａｎｄ８ｃｍｔｏｔｈｅｎｏｒｔｈ．ＴｈｅｅａｒｔｈｑｕａｋｅａｆｆｅｃｔｅｄｔｈｅａｒｅａａｒｏｕｎｄｅｐｉｃｅｎｔｅｒｉｎｃｌｕｄｉｎｇＭａｄｕｏ，Ｘｉｎｉｎｇ，Ｄｕｌａｎ，Ｄｅｌｉｎｇｈａｉｎｔｈｅｎｏｒｔｈ，ａｎｄＺｅｂｒａａｎｄＧａｎｚｉａｒｅａｓｉｎｔｈｅｓｏｕｔｈ．Ｆｒｏｍｔｈｅｃｏｍｐａｒｉｓｏｎｏｆｔｈｅｒｅｓｕｌｔｓｏｆｔｈｅｓｔａｔｉｃｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄ
ｆｏｒｔｈｅ５ｈｏｕｒｓａｎｄｍｕｌｔｉ ｄａｙｄａｔａ
，ｉｔｉｓｂｅｌｉｅｖｅｄｔｈａｔｔｈｅｐｏｓｔ ｓｅｉｓｍｉｃｄｅｆｏｒｍａｔｉｏｎｔａｋｉｎｇｐｌａｃｅｉｎｔｈｅｎｅａｒｆｉｅｌｄｉｓｓｉｇｎｉｆｉｃａｎｔ，ａｎｄｃｏｎｔｉｎｕｏｕｓａｔｔｅｎｔｉｏｎｉｓｒｅｑｕｉｒｅｄｉｎｔｈｅｌａｔｅｒｓｔａｇｅｓ．
Ｋｅｙｗｏｒｄｓ　ｔｈｅＭａｄｕｏＭＳ７ ４ｅａｒｔｈｑｕａｋｅｉｎＱｉｎｇｈａｉＰｒｏｖｉｎｃｅ，ＧＮＳＳ，ｃｏｓｅｉｓｍｉｃｄｉｓｐｌａｃｅｍｅｎｔ
〔作者简介〕　李经纬，男，１９８９年生，２０１６年于中国地质大学（武汉）获测绘科学与技术硕士学位，工程师，
主要从事ＧＰＳ数据处理与地壳形变分析研究，Ｅ－ｍａｉｌ：ｌｉ ｊｉｎｇｗｅｉ＠ｃｕｇ．ｅｄｕ．ｃｎ。

４８０１。