基于GPS的2021年玛多MW7.4地震同震静、动态形变场

2022-10-10 01:53熊仕昭王迪晋熊赵利江乔学军
大地测量与地球动力学 2022年10期
关键词:波速测站断裂带

熊仕昭 刘 刚,2 林 牧 王迪晋熊 维 赵利江 赵 斌 乔学军

1 中国地震局地震大地测量重点实验室,武汉市洪山侧路40号,430071 2 防灾科技学院,河北省三河市学院街465号,065201 3 青海省基础测绘院,西宁市黄河路13号,810001

据中国地震台网测定,2021-05-22 02:04青海省果洛藏族自治州玛多县发生MW7.4地震,震中位于34.61°N、98.37°E,震源深度约为17 km,是2008年汶川地震以来中国大陆发生的震级最大的地震。截至2021-05-24 10:00共记录到余震1 872次,其中5.0~5.9级余震1次,4.0~4.9级余震10次,3.0~3.9级余震30次。

此次地震发生于青藏高原中北部地区的巴颜喀拉块体。该块体东西长2 200 km,南北较窄(最小宽度80 km),沿北西走向,是新生代以来非常具有代表性的侧向挤出活动块体之一[1-2]。巴颜喀拉块体东西南北边界处分别有龙门山逆冲断裂带、阿什库勒-平顶山正滑断裂带和硝尔库勒走滑-正滑断裂带、玛尼-玉树-鲜水河走滑断裂带、东昆仑走滑断裂带,平均海拔5 000 km左右,是目前青藏高原强震活动的主体区域,在整个青藏高原中地势最为平坦[3]。自1997年玛尼7.5级地震以来,中国大陆所有7级以上的地震均发生于此区域,包括2001年昆仑山口西(MS8.1)、2008年于田(MS7.3)、2008年汶川(MS8.0)、2010年玉树(MS7.1)、2013年芦山(MS7.0)、2014年于田(MS7.3)、2017年九寨沟(MS7.0)等地震,且目前强震活动仍在继续[4-5]。但巴颜喀拉块体内部相对稳定,仅在1947年发生过一次7.7级达日地震[6]。本次地震震中位于东昆仑断裂带南部活动速度相对较低的昆仑山口-江错断裂带的东南段(图1)。该断裂带强震复发周期偏长[7-8],且此前对该区域的构造活动研究较少,因此从不同角度分析玛多地震具有重要的研究价值。

本文通过解算青海省CORS网记录到的2021年玛多地震的同震GPS数据,得到同震永久位移、动态形变波形,并据此对此次地震的同震形变场特征进行分析。

1 静态数据处理与同震永久形变场

收集玛多地震震中附近300 km范围内43个青海省CORS网站点地震前后共5 d的数据,采样率为30 s。为确保参考框架的统一与自洽,保证同震位移解算结果的精度与可靠性,加入周边地区15个IGS台站的数据,得到由58个站点组成的一个观测网。以24 h为一个测段,使用GAMIT/GLOBK 10.71软件对原始观测数据进行基线解算获得单日松弛解,该软件利用双差观测值可以较好地消除GPS卫星钟差和接收机钟差,输出的单日解文件包括定位和定轨数据信息,是求解玛多地震同震位移的可靠观测值。解算策略中,采用最新的地球自转参数、太阳星历、月亮星历、章动参数及天线相位中心参数文件等,并加入电离层模型改正。采用分段线性插值法,每2 h估算 一次对流层天顶延迟参数,同时估算大气的延迟水平梯度。为避免系统误差,并保证和IGS精密轨道产品间较好的自洽性,采用全球海潮模型FES2004修正海潮荷载引起的站点地壳形变,并考虑其导致的地球质心变化。在网平差处理阶段,将上一步解算的松弛约束区域解与SOPAC(Scripps Orbital and Permanent Array Center)生成的全球解相结合,并通过选定全球范围内分布的15个IGS参考站,采用平移、旋转和缩放的七参数相似变换将得到的单日解坐标转换到ITRF2014框架下。因数据质量较好且完整,故对于震前只采用4 d的数据,而采用震后1 d的数据是为了规避较大的余震对处理结果的影响。对于这些连续站,将震后1 d的坐标值与震前4 d的坐标平均值作差,所得即为此次玛多地震的同震位移量。水平位移形变场见图1,三维形变数值结果见表1。

误差椭圆置信区间为70%图1 玛多MW7.4地震GPS水平方向同震形变场Fig.1 GPS horizontal coseismic deformation field of Maduo MW7.4 earthquake

表1 基于静态GPS解算的同震位移

GPS观测数据在水平方向精度较高,误差为3~5 mm,垂向误差约为其2~3倍。从图1看出,KANQ、MADU、JDUO这3个站距离震中最近,位移量也较大,水平向均接近0.3 m。说明此次地震在玛多县城及其附近区域产生显著的水平同震位移,该结果与李志才等[9]低频GPS数据解算的形变特征一致。结合紧靠断裂面地段的地形图可以判断,位于断层南侧的JDUO等站发生朝南东向的位移,位于断层北侧的KANQ、MADU、HSHX等站发生朝北西向的位移,符合此次玛多地震的左旋走滑机制,与断层所在的巴颜喀拉块体具有向东运动的特征[3]相符。在发震断层的东端,可以看到BUDR和GAND站除向东运动外还有向南运动的趋势,可能是走滑错动的端部拉张变形所致。距离震中150 km以内的站点均能反映出cm级的位移;从距离震中200 km开始,如LAJA、YUSH等站,同震位移快速衰弱,只能反映出mm级的位移;对于距离震中300 km的测站,如YANH、GERM、JZHA等,不能反映出明显的同震位移,可认为这些远场测站受到此次同震的影响偏小。综上,本次玛多MW7.4强震的同震影响集中在震中距300 km的区域范围内。

2 动态数据处理与同震动态波形

垂直虚线为发震时刻,倾斜虚线为本次地震的走时曲线图2 震时地表运动图Fig.2 Earth surface movement during earthquake

表2 测站地震波速信息

从图2看出,测站记录到了清晰的同震位移波形。GCMT给出的发震时刻为UTC 18:04:29,在震后11 s JDUO站首次产生明显的峰谷变化,随后各个测站的位移波形显示出地震波由震中向外传播的过程,并随着震中距的增加,位移量下降,表明传播过程中伴随着地震波能量的衰减。从图3看出,JDUO和KANQ站的震时三维运动轨迹显示出一定的逆进椭圆样式,同时伴随着水平向振荡,表明记录到了层状介质中激发的Love波和Rayleigh波。

图3 JDUO和KANQ站三维运动轨迹Fig.3 Three-dimensional motion trajectory at JDUO and KANQ stations

确定测站波形的初动时刻和震中距后,可计算出地震波传播速度,结果见表2,表中瞬时位移时刻以UTC 18:00(即动态解算起始时间)为起点,采用误差阈值法来确定每个测站的地震波初至时刻[10],其与震源发震时刻之差即为地震的传播时间。方荣新[11]假定地震到达各个测站的波速相等,基于GPS观测数据对汶川地震震中进行反演,计算得到的平均地震波速为3.939 km/s。本文得到的地震波速约2.8~4.9 km/s,与地震学估计的一维速度模型中S波速度基本一致[2]。距离震中较远的BUDR和GAND站的波速大于其他近场站,可能是破裂向东传播的方向性或者较高的破裂速度导致的。上述特征表明,地震波能量传播不仅与震中距离有关,同时也会受到地壳介质分布及震源破裂的动态效应、方向效应的影响。本文得到的S波速仅是在有限的数据中获取的估算值,当数据量足够时,依据本文的解算精度可以进一步反演出可靠的震中信息、发震时刻的波速及其随传播路径的时间变化。

Gutenberg[12]给出远场面波的峰值位移和震中距与震级之间的关系为:

M=lg(A)+1.66lg(Δ)+2.0

(1)

式中,A是以μm为单位的水平峰值位移,Δ是以(°)为单位的震中距。

Fang等[13]使用3个地震实例证明了该公式对于含有近场站的GPS位移仍然是适用的。本文基于该公式获取15个震幅明显的站点的震时地表峰值位移与震中距之间的关系,结果见图4,其中,图4(a)表示峰值位移随震中距的衰减,图4(b)表示在10的幂次方尺度下峰值位移与震中距之间呈斜率为-1.66的线性关系,蓝色直线的值即为拟合震级。依据上述公式,本文估计的震级区间在6.8~7.6,拟合平均震级为7.35±0.15。为了研究该经验公式估算玛多地震可靠震级所需要的时间,本文根据站点解算数据作为时间的函数,得到图5。结果表明,地震发生后33 s内可实现震级的估算,时间主要取决于站点位置与地震波的速度。由于该估算震级的大小受数据解算质量、测站地质环境等的影响,因此存在一定的波动性,但总体能在70 s后稳定于一个大致的范围内。

图4 峰值位移-震中距函数图Fig.4 Function diagram of peak displacements and epicentral distance

图5 震级-震后时间函数图Fig.5 Function diagram of estimated magnitudeand time after earthquake onset

3 结 语

本文基于青海省CORS网连续GPS观测数据,利用GAMIT/GLOBK软件对玛多地震同震的静态数据及动态数据进行解算,分别获取地表的同震位移场与震时地表运动状态,结果表明:

1)位于断层北侧的近场站朝北西向运动,位于断层南侧的近场站朝南东向运动,符合左旋走滑机制,震中150 km以内的站点均能反映出cm级的同震位移。

2)近场动态GPS时序显示出明显的地震波传播特征,估计得到的地震波速与该区位的S波速度相当,断层面东端站点估计的波速明显大于其他站点,可能与破裂方向性或者破裂速度有关。

3)由震级经验公式估计得到震级为7.35。若进行实时估算,在震后70 s内可得到稳定的震级。

致谢:GPS数据来源于青海省CORS网,玛多地震震源机制解来源于GCMT,余震数据来源于中国地震台网,数据处理软件为GAMIT,在此一并表示感谢。

猜你喜欢
波速测站断裂带
花岗岩物理参数与纵波波速的关系分析
波速球帮你增强核心、协调性和平衡感(下)
冷冻断裂带储层预测研究
WiFi室内定位测站布设优化的DOP数值分析
受载岩体破坏全过程声波响应特征及工程意义
海洋潮汐负荷对精密单点定位的影响研究
依兰—伊通断裂带黑龙江段构造运动特征
利用背景噪声研究福建金钟库区地壳介质波速变化
浅析董事会断裂带