1 引言
2022年3月17日甘肃张掖市肃南县5.1级地震发生后，韩立波、蒋长胜、郭祥云、张辉、中国地震台网中心等采用不同方法和资料得到的该地震的多个震源机制解。这些震源机制有一定的离散度，为地震动力学分析或其他应用带来抉择的困难。这些结果都是震源错动方式的一种测量，因此可以按照多种测量结果给出一个中心值供以后的地震发生背景、地震应力触发、地壳应力场分析以及地震前应力方向改变的地震前兆研究。这里整理了各个机构给出的震源机制解，求出了与所有测定的震源机制的差别平方和最小的一个解作为中心震源机制解，并计算了该地震在周围地表产生的位移场和应变场。
2 多个震源机制的中心解的确定
根据国外机构发布的该地震震源机制结果整理得到表1。我们分别以各个震源机制为初始解得到的中心震源机制给出的标准差（表1第5列）大体一致（在小数点4位后有一定涨落），表明采用这种方法得到的解是稳定的。尽管如此，本研究将各个机构测定的震源机制分别作为初始解，比较得到标准差最小的解作为最终结果。发现将韩立波、蒋长胜得到的震源机制作为初始解得到的震源机制的标准差最小。本研究以此（节面I走向343.03°，倾角75.75°，滑动角159.43°，节面II走向78.31°，倾角70.09°，滑动角15.18°）作为最终结果，P轴走向31.47°，倾伏角3.83°,不确定范围分别为18.00~45.00°和-6.21~13.91°;T轴的走向299.72°，倾伏角24.52°,不确定范围分别为286.25~313.25°和15.13~34.12°;B轴的走向129.77°倾伏角65.15°,不确定范围分别为99.67~161.12°和56.16~76.87°。得到的中心震源机制的最小空间旋转角见表1第6列。所得到的中心震源机制及其不确定性绘于图1。从表和图中可以看出，该地震震源机制解距中心解的空间旋转角最大达8.31°，最小空间旋转角为1.85°。
图1甘肃张掖市肃南县5.1级地震的中心震源机制解(a)及空间三维辐射花样(b)
(a)中的黑色弧线表示中心震源机制的两个节面，绿色弧线覆盖区域为其不确定范围；红色、蓝色和黄色的点表示中心震源机制解的P轴、T轴和B轴，其周围对应颜色的封闭曲线表示其不确定性范围；绿点和黑点表示各个机构得到的震源机制的P轴和T轴的投影；紫色弧线表示各个机构和作者得到的震源机制节面。(b)中的压缩区域和膨胀区域分别用蓝色和红色表示。
3 震源机制中心解的空间表示
地震的震源机制通常采用震源球的辐射花样来表示。在此我们分别用震源球的上半球和下半球的动画来表示（图2），在震源球上采用P波辐射的振幅相对大小的颜色来填充。为立体表现震源球的辐射花样，我们使震源球水平旋转，并且在各个不同象限采用P波辐射相对振幅大小和方向绘制在震源球上，这样可更为直观地表现震源的空间辐射，见图3。
图2甘肃张掖市肃南县5.1级地震下半球、上半球
图3甘肃张掖市肃南县5.1级地震的中心震源机制解的空间辐射花样表示
图中颜色为P波辐射花样在震源球面上的填充，红色表示向外，蓝色表示向内，从红到蓝表示P波辐射花样由向外最大逐渐过渡到向内最大。箭头表示P波辐射花样振幅的相对大小。
4 该地震在周围产生的位移场与应变场
本次地震震中位于东经97.66°，北纬39.02°，震源深度约为9公里。统计本次地震50公里范围内，1970年以来历史地震M3.0级及以上地震共发生30次，其中3.0~3.9级地震21次，4.0~4.9级地震6次，5.0~5.9级地震3次，最大为1975年01月04日青海哈拉湖北5.3级地震，距本次震中约24公里。100公里范围内M3.0级及以上地震共发生112次，其中3.0~3.9级地震91次，4.0~4.9级地震15次，5.0~5.9级地震6次，最大为2002年12月14日甘肃玉门5.9级地震，距本次震中约92公里。
图4 研究区域地形图
为了更好的了解甘肃张掖市肃南县5.1级地震对周围地区的影响，我们基于弹性半空间模型Okada（1992），选择中心震源机制解节面Ⅱ作为发震断层面，根据Wells和Coppersmith（1994）给出的统计公式估计了断层面大小和滑动量（长度：3.9084km，宽度：4.1400km，滑动量：10.3129cm）计算了甘肃张掖市肃南县5.1级地震在周围地区产生的地表同震位移场、面应变、体应变、北向应变、东向应变和北东向应变（图5~7）。
图5 甘肃张掖市肃南县5.1级地震产生的同震位移场
图中白色方框代表断层在地表的投影，箭头代表此次地震所产生的水平位移，颜色代表垂直位移，上升为正
从图5可以看出，此次地震周围地区的位移场具有以下初步特征：从地表的水平位移场来看，震中东北，西南物质涌入震中，而东南的物质明显向外涌出，在发震断层附近则呈现明显的走滑机制。与水平位移场相对应，垂直位移场在震中东南部表现为隆升，而西北部表现为沉降。地震为走滑型地震，水平位移大于垂直位移，以上特征符合走滑型地震所产生位移场的特征。
图6 甘肃张掖市肃南县5.1级地震产生的各分量应变
（a）体应变（b）北向应变（c）东向应变（d）北东向应变
图中底色代表应变大小，拉张为正，单位为1e-9.
从图6可以看出此次地震产生的体应变,在震中东南侧附近呈现伸张，而在震中西南，东北两侧呈现压缩，北向应变、东向应变、北东向应变均呈现明显八花瓣样式分布。北向应变在震中东北侧和东南侧和西侧压缩，东侧伸张。东向应变在震中附近东侧伸张，东侧较远处压缩，震中西侧压缩，震中北侧压缩较小。北东向应变在震中西南和东北两侧压缩，西南侧压缩大于东北侧压缩，在震中东南侧伸张。
图7 甘肃张掖市肃南县5.1级地震产生的水平主应变和面应变场
图中黑色箭头和白色箭头分别表示水平主压应变和水平主张应变（单位为1e-9）;底色表示水平面应力，拉张为正
从图7可以看出面应变的分布与体应变形态大体类似，表现为面应变为体应变的继承。在面应变为负（蓝色）区域内，物质的运动方向（图5）与主压应变大体一致，而在面应变为正（红色）区域内，物质的运动方向（图5）与主张应变大体一致。
本文所用震源机制解来自韩立波、蒋长胜、郭祥云、张辉、中国地震台网中心特此致谢。图件由基于MATLAB的FM3Dplot软件绘制。
