1 引言
根据中国地震台网测定,2026年10月22日22时33分,西藏那曲发生7.0级地震。震中位于东经93.79°,北纬32.4°,震源深度约为10公里。地震发生后,中国地震台网中心、中国地震局地球物理研究所张喆、北京大学张勇、中国科学院青藏高原研究所王卫民和国外机构GCMT、GFZ、NEIC、CPPT、AUST等采用不同方法和资料得到的该地震的多个震源机制解。这些震源机制有一定的离散度,为地震动力学分析或其他应用带来抉择的困难。这些结果都是震源错动方式的一种测量,因此可以按照多种测量结果给出一个中心值供以后的地震发生背景、地震应力触发、地壳应力场分析以及地震前应力方向改变的地震前兆研究。这里整理了各个机构给出的震源机制解,求出了与所有测定的震源机制的差别平方和最小的一个解作为中心震源机制解,并计算了该地震破裂对周围产生的地表变形场,供分析余震发展趋势参考。
2 多个震源机制中心解的确定
根据国内外机构网站和多位作者发布的该地震震源机制结果整理得到表1。我们分别以各个震源机制为初始解得到的中心震源机制给出的标准差(表1第5列)大体一致(在小数点3位后有一定涨落),表明采用这种方法得到的解是稳定的。尽管如此,本研究将各个机构测定的震源机制分别作为初始解,比较得到标准差最小的解作为最终结果。发现以AUST得到的震源机制作为初始解得到的震源机制的标准差最小。本研究以此(节面Ⅰ走向18.78°,倾角51.94°,滑动角-87.86°,节面Ⅱ走向195.31°,倾角38.11°,滑动角-92.73°)作为最终结果,P轴走向301.09°,倾伏角82.88°,不确定范围分别为242.48~454.97°和77.28~95.64°;B轴的走向197.46°,倾伏角1.69°,不确定范围分别为184.33~210.33°和-8.39~12.00°;T轴的走向107.26°,倾伏角6.92°,不确定范围分别为94.13~120.13°和-5.95~20.04°。得到的中心震源机制和各个机构测定震源机制的最小空间旋转角见表1第6列。所得到的中心震源机制及其不确定性绘于图1。从表和图中可以看出,该地震震源机制解距中心解的空间旋转角最大达24.95°,最小空间旋转角为2.93°。这些数据表明不同机构和作者得到的震源机制解较为集中。
表1 不同机构给出西藏那曲发生7.0级地震震源机制及得到的中心震源机制解的标准差


图1 西藏那曲7.0级地震的震源机制中心解地震的中心震源机制解(a)及空间三维辐射花样(b)
(a)中的黑色弧线表示中心震源机制的两个节面,绿色弧线覆盖区域为其不确定范围:红色、蓝色和黄色的点表示中心震源机制解的P轴、轴和B轴,其周围对应颜色的封闭曲线表示其不确定性范围;绿点和黑点表示各个机构得到的震源机制的P轴和T轴的投影;紫色弧线表示各个机构和作者得到的震源机制节面。(b)中的压缩区域和膨胀区域分别用红色和蓝色表示。
3 震源机制中心解的空间表示
地震的震源机制通常采用震源球的辐射花样来表示。在此我们分别用震源球的上半球和下半球的动画来表示(图2),在震源球上采用P波辐射的振幅相对大小的颜色来填充。为立体表现震源球的辐射花样,我们使震源球水平旋转,并且在各个不同象限采用P波辐射相对振幅大小和方向绘制在震源球上,这样可更为直观地表现震源的空间辐射,见图3。


图2 西藏那曲7.0级地震的中心震源机制解的上半球、下半球表示

图3西藏那曲7.0级地震的中心震源机制解的空间辐射花样表示
图中颜色为P波辐射花样在震源球面上的填充,红色表示向外,蓝色表示向内,从红到蓝表示P波辐射花样由向外最大逐渐过渡到向内最大。箭头表示P波辐射花样振幅的相对大小。
4 西藏那曲7.0级地震在周围产生的位移场与应变场
本次西藏那曲7.0级地震发生在青藏高原上,为印度板块向欧亚板块碰撞导致青藏高原隆升和地震最为活跃的地区之一。高地形的重力势能导致高原东西向扩展,产生了大量正断和走滑地震。西藏那曲7.0级地震震中以东北仅约2km发育着全新世活动的巴青-类乌齐断裂。自1900年以来50km范围内M3.0级及以上地震共发生38次,其中7.0~7.9级地震1次,其中6.0~6.9级地震1次,其中5.0~5.9级地震4次,4.0~4.9级地震7次,3.0~3.9级地震25次,最大为本次地震。100km范围内M3.0级及以上地震共发生136次,其中7.0~7.9级地震1次,其中6.0~6.9级地震1次,其中5.0~5.9级地震20次,4.0~4.9级地震27次,3.0~3.9级地震87次,最大为本次地震。从震源机制解,判断地震为正断层型地震。

图4 研究区域构造背景
为了了解西藏那曲7.0级地震对周围地区的影响,我们根据USGS网站(https://earthquake.usgs.gov/earthquakes)提供的破裂模型,计算了西藏那曲7.0级地震在周围地区产生的地表同震位移场、面应变、体应变、北向应变、东向应变和北东向应变(图5~7)。

图5 西藏那曲7.0级地震产生的同震位移场
图中箭头代表此次地震所产生的水平位移,颜色代表垂直位移,上升为正。
从图5可以看出,此次地震近处的位移场具有以下初步特征:从地表的水平位移场来看,震中西侧与东侧的物质向外涌出,而南侧与远端北侧的物质涌入震中。与水平位移场相对应,垂直位移场在震中远端西北侧与北侧表现为隆升,而震中及震中南侧、西南侧表现为沉降。

图6 西藏那曲7.0级地震产生的各分量应变
(a)体应变(b)北向应变(c)东向应变(d)北东向应变
图中底色代表应变大小,拉张为正,单位为10-9。
从图6可以看出此次地震在地表产生的体应变,在震中及东西两侧呈现压缩,在震中西南侧及远端北侧呈现拉张。北向应变在震中周围呈现八花瓣压缩、拉张相间的模式,东向应变呈现东西向侧挤压,南侧、北侧和西南侧远方呈现拉张,北东向应变东南侧、西南侧和北侧远处呈现拉张,东北侧和西侧呈现压缩。

图7 西藏那曲7.0级地震产生的水平主应变和面应变场
图中黑色箭头和白色箭头分别表示水平主压应变和水平主张应变(单位为10-9);底色表示水平面应力,拉张为正。
从图7可以看出面应变的分布与体应变形态大体类似,表现为面应变为体应变的继承。在离开断层区域,面应变为负(蓝色)区域内,物质的运动方向(图5)与主压应变大体一致,而面应变为正(红色)区域内,物质的运动方向(图5)与主张应变大体一致。
致谢:本文所用震源机制解来自中国地震台网中心、中国地震局地球物理研究所张喆、北京大学张勇、中国科学院青藏高原研究所王卫民和国外机构GCMT、GFZ、NEIC、CPPT、AUST,特此致谢。震源机制动画由研究小组基于MATLAB开发的FM3Dplot软件绘制。如若发现有错误,请反馈给我们,我们将及时修改,联系方式:13573902798@163.com