1 研究背景

2021年5月22日02时04分青海省果洛藏族自治州玛多县发生M_S 7.4地震，震中位于（98.34°E，34.59°N），震源深度17 km。震后发生多次余震，属于典型的主震—余震型地震。该地震序列并不是沿着平直断层发生的，东南段出现了沿走向的东向偏转，这种转向处的应力场相对于破裂断层的应力场有哪些改变，破裂性质是否发生变化，这是该区域地震动力学应该考虑的重要问题。本研究根据该地震序列的震源机制求解该地区的非均匀应力场，并给出动力学解释。

2 资料和方法

搜集前人提供的玛多地震序列震源区的150条震源机制数据，采用万永革（2019）提出的同一地震多个震源机制中心解确定算法进而得到震源机制中心解。然后按照沿断层自东南向西北的原则将玛多震源区震源机制解事件位置进行排列，接着由东南至西北为事件分区，由此可以获得断层区不同区域的应力状态。之后采用滑动窗口法对震源机制进行分区，沿断层走向将研究区域事件分成8组，每个窗口含有30—31个震源机制，前移一个窗口更新13个震源机制解数据。最后采用网格搜索法求解每个区域的地壳应力场。

3 研究结果

反演结果展示在图 1中，沿破裂区走向观察，破裂区的（a）区域，R = 0.9主要以挤压应力为主，中间轴与张应力轴表现为拉张作用，以至于中间轴与张应力轴不易区分开来；破裂区的（b）区域R值减少为0.7，相比较与（a）区域中间轴与张应力轴逐渐区分开来，此时张轴的倾伏角也突然增大，向垂直向偏转；破裂区（c）区域R值突然减小为0.3，主要以拉张应力为主，中间轴与张应力轴表现为压缩作用，且张轴倾伏角继续增大，表现为近乎垂直向；破裂区（d）区域R值为0.6，该区域的张应力轴与中间应力轴又变得不能区分开来，且张轴倾伏角急剧减小，数值和（a）区域相近；破裂区（e）区域R值为0.4，主要以拉张应力为主，张轴倾伏角轻微增加；破裂区（f）区域R值为0.5，表现为拉张应力和挤压应力平衡的状态，中间轴和两轴可以很好的区分开来，且张轴倾伏角轻微增加；破裂区（g）区域R值为0.7，中间轴和张应力轴表现为拉张作用，且变得不易区分开来，张轴倾伏角也与（f）区域相近；破裂区（h）区域R值为0.1，中间轴与压应力轴表现为压缩作用，两轴不易区分开来，且张轴倾伏角降低许多。这些分区的应力场结果大部分与整体的应力场结果相差较小，是整体应力场的范围内进行的较小的扰动，表明了整体应力场和分区非均匀应力场的一致性。

将分区反演的非均匀应力场总结为：破裂区的应力场的拉张轴均呈现不同程度的张轴倾伏角，区域（a）（b）（d）（e）（f）（g）（h）的应力场反演结果分别表现为带有逆冲分量的走滑型，而区域（c）的应力场反演结果表现为逆断型，这与整体应力场所表现的类型不符，而（c）区域正位于地震序列走向东向偏转的位置。

4 结束语

地震序列东南段走向偏转处（c）区域的张应力轴近乎垂直，其他部分的应力场的张应力轴也均呈现不同程度的张轴倾伏角，表现为兼有逆冲分量的走滑型。破裂区（c）的应力场结果与整体应力场具有较大区别且对地震序列走向产生了较大影响（图 2），由此推测该区域的明显应力场变化可能是由于青藏高原在印度板块的北向推挤下，巴颜喀拉块体内部昆仑山口—江错断裂破裂区物质向NW向移动时，受到了东侧相对高速体的阻挡，使得断层东侧局部受到了挤压变形而产生逆冲型应力场。

区域（a）和区域（d）的地震序列走向分布相较于其他区域逆时针旋转了大约30°。断层东西两侧的泊松比差异明显，东侧的泊松比明显高于西侧，因此东部物质相对西部而言的抗形变能力更差，更容易发生形变，在东侧受到挤压产生逆冲应力类型后，由于东侧的泊松比更高，更易发生形变，所以地震继续向东侧破裂，且破裂强度更大，以至于向东延伸横穿玛多—甘德断裂，具有“裁弯取直”的特点，受到局部应力调整由此产生的新的震源破裂机制，新的震源破裂机制应力场表现为兼逆冲分量的走滑型。

文中图件采用GMT6软件和MATLAB软件绘制而成，中国科学技术大学张海江研究小组为本研究提供地球岩石圈速度结构模型电子文件，特此感谢！