1引言
根据中国地震台网测定,2024年3月05日10时07分,青海省杂多县发生5.3级地震,地震发生后,中国地震台网中心、中国地震局地球物理研究所和国外机构等采用不同方法和资料得到的该地震的多个震源机制解。这些震源机制有一定的离散度,为地震动力学分析或其他应用带来抉择的困难。这些结果都是震源错动方式的一种测量,因此可以按照多种测量结果给出一个中心值供以后的地震发生背景、地震应力触发、地壳应力场分析以及地震前应力方向改变的地震前兆研究。这里整理了各个机构给出的震源机制解,求出了与所有测定的震源机制的差别平方最小的一个解作为中心震源机制解,供分析余震发展趋势参考。
2 多个震源机制的中心解的确定
根据国内外机构网站和多位作者发布的该地震震源机制结果整理得到表1。我们分别以各个震源机制为初始解得到的中心震源机制给出的标准差(表1第4列)大体一致(在小数点3位后有一定涨落),表明采用这种方法得到的解是稳定的。尽管如此,本研究将各个机构测定的震源机制分别作为初始解,比较得到标准差最小的解作为最终结果。发现将王勤彩得到的震源机制作为初始解得到的震源机制的标准差最小。本研究以此(节面I走向336.88°,倾角77.57°,滑动角-168.59°,节面II走向244.39°,倾角78.86°,滑动角-12.68°)作为最终结果,P轴走向200.50°,倾伏角16.79°,不确定范围分别为181.18~ 219.18°和3.68~ 30.30°;B轴的走向23.72°,倾伏角73.19°,不确定范围分别为-51.46~ 99.70°和67.78~ 89.90°;T轴的走向290.77°,倾伏角0.89°,不确定范围分别为271.45~ 309.45°和-12.87~ 14.27°。得到的中心震源机制和各个机构测定震源机制的最小空间旋转角见表1第5列。所得到的中心震源机制及其不确定性绘于图1。从表和图中可以看出,该地震震源机制解距中心解的空间旋转角最大达31.52°,最小空间旋转角为6.67°。这些数据表明不同机构和作者得到的震源机制解较为集中。
图1 青海省杂多县发生5.3级地震的中心震源机制解(a)及空间三维辐射花样(b)
(a)中的黑色弧线表示中心震源机制的两个节面,绿色弧线覆盖区域为其不确定范围;红色、蓝色和黄色的点表示中心震源机制解的P轴、T轴和B轴,其周围对应颜色的封闭曲线表示其不确定性范围;绿点和黑点表示各个机构得到的震源机制的P轴和T轴的投影;紫色弧线表示各个机构和作者得到的震源机制节面。(b)中的压缩区域和膨胀区域分别用蓝色和红色表示
3 震源机制中心解的空间表示
地震的震源机制通常采用震源球的辐射花样来表示。在此我们分别用震源球的上半球和下半球的动画来表示(图2),在震源球上采用P波辐射的振幅相对大小的颜色来填充。为立体表现震源球的辐射花样,我们使震源球水平旋转,并且在各个不同象限采用P波辐射相对振幅大小和方向绘制在震源球上,这样可更为直观地表现震源的空间辐射,见图3。
图2 青海省杂多县发生5.3级地震的中心震源机制解的下半球、上半球表示
图3 青海省杂多县发生5.3级的中心震源机制解的空间辐射花样表示
图中颜色为P波辐射花样在震源球面上的填充,红色表示向外,蓝色表示向内,从红到蓝表示P波辐射花样由向外最大逐渐过渡到向内最大。箭头表示P波辐射花样振幅的相对大小。
4 青海玉树州杂多县5.3级地震在周围产生的位移场与应变场
本次地震周边5公里内的村庄有扎钦,50公里内无乡镇驻地分布,100公里内的乡镇有索加乡、永曲乡、查旦乡。震中距西藏聂荣县171公里、距西藏安多县187公里、距西藏索县196公里、距西藏巴青县203公里、距杂多县225公里,距西藏那曲市245公里,距拉萨市467公里。
图4 研究区域地形图
为了了解青海玉树州杂多县5.3级地震对周围地区的影响,我们基于弹性半空间模型Okada(1992),根据Wells和Coppersmith(1994)给出的统计公式估计了断层面大小和滑动量(长度:5.2km,宽度:4.69km,滑动量:13.66cm),根据中国地震台网中心提供的深度10km,选择震源机制中心解节面Ⅰ。计算了青海玉树州杂多县5.3级地震在周围地区产生的地表同震位移场、面应变、体应变、北向应变、东向应变和北东向应变(图5~7)。
图5 青海玉树州杂多县5.3级地震产生的同震位移场
图中白色方框代表断层在地表的投影,箭头代表此次地震所产生的水平位移,颜色代表垂直位移,上升为正
从图5可以看出,此次地震近处的位移场具有以下初步特征:从地表的水平位移场来看,震中东南侧和西北侧的物质向外涌出,而东北侧的物质涌入震中,发震断层附近呈现明显的走滑机制。与水平位移场相对应,垂直位移场在震中东南侧和西北侧表现为隆升,而震中东北侧表现为沉降。从震源机制解判断地震为走滑型地震。
图6 青海玉树州杂多县5.3级地震产生的各分量应变
(a)体应变(b)北向应变(c)东向应变(d)北东向应变
图中底色代表应变大小,拉张为正,单位为1e-9.
从图6可以看出此次地震在地表产生的体应变,在震中附近东北侧呈现压缩,而在震中及东南侧和西北侧呈现拉张。北向应变在震中北侧压缩,在北侧较远处、东侧和西侧拉张。东向应变在震中北侧压缩,在震中东侧和西侧拉张。北东向应变在震中附近东北侧呈现压缩,在东北侧较远处、东南侧和西北侧区域呈现拉张。
图7 青海玉树州杂多县5.3级地震产生的水平主应变和面应变场
图中黑色箭头和白色箭头分别表示水平主压应变和水平主张应变(单位为1e-9);底色表示水平面应力,拉张为正
从图7可以看出面应变的分布与体应变形态大体类似,表现为面应变为体应变的继承。在离开断层区域,面应变为负(蓝色)区域内,物质的运动方向(图5)与主压应变大体一致,而面应变为正(红色)区域内,物质的运动方向(图5)与主张应变大体一致。
致谢:本文所用震源机制解来自中国地震局地球物理研究所郭祥云,GFZ,罗钧和赵翠萍,李启雷,张喆,中国地震台网中心,韩立波和蒋长胜,王勤彩,GCMT,特此致谢。图件由基于MATLAB的FM3Dplot软件绘制。如若发现有错误,请反馈给我们,我们将及时修改,QQ:1978250163