2022年墨西哥米却肯州7.6级地震序列的触发关系及发震动力学探讨.pdf

书书书第 卷第 期（）年 月中 国 地 震 关兆萱，万永革，黄少华，等 年墨西哥米却肯州 级地震序列的触发关系及发震动力学探讨 中国地震，（）：年墨西哥米却肯州 级地震序列的触发关系及发震动力学探讨关兆萱）万永革，）黄少华）李枭）冯淦）防灾科技学院，河北三河 ）河北省地震动力学重点实验室，河北三河 摘要 年 月 日墨西哥米却肯州发生 地震，该地震位于北美板块与科科斯板块交汇处。为进一步研究本次地震对周围断层和后续地震的影响，剖析该地震的孕震背景和发生条件，首先使用包括 在内的多个国外地震机构得到的墨西哥 地震以及后续的两次 地震的震源机制解，分别计算出震源机制中心解，然后通过计算主震产生库伦破裂应力变化来研究本次地震对后续两次地震的触发作用，最后收集了 年 月 日至 年 月 日发生在本次地震附近的 的 条震源机制解数据，反演该地区的构造应力场，并模拟在局部构造及其作用下产生的各种断层形状及相对正应力和剪应力的分布情况。研究结果表明：震源机制中心解表明主震和 级余震为逆冲型地震，而 级余震的震源机制为正断型地震；本次地震对 级和 级余震的库伦破裂力变化均超过 的阈值，表明这两次强余震可能是在主震的触发下发生的；震源区应力场的主压应力轴为 向，主张应力轴近乎垂直。本地震序列的断层破裂近乎沿应力场的最大剪应力平面发生，最大限度地释放了构造运动积累的应力。关键词：墨西哥米却肯 地震震源机制中心解构造应力场库伦破裂力变化 文章编号 （）中图分类号 文献标识码 收稿日期 修定日期 项目类别 国家自然科学基金（、）、中央高校科研业务费专项（）、河北省地震动力学重点实验室开放基金项目（）、河北省高等学校自然科学技术项目（）共同资助 作者简介 关兆萱，男，年生，硕士研究生，主要从事构造应力场、地震应力触发方面的研究工作。：万永革，通讯作者，男，年生，研究员，主要从事构造应力场、地震应力触发等方面的研究。：引言根据美国地质调查局（，简称 ）测定，世界标准时间（，简称 ）年 月 日 时 分，墨西哥西海岸米却肯州发生 地震，震中位置为 ，震源深度约为 。月 日 时 分（）在主震以西约 处（，）发生 地震，月 日 时 期关兆萱等：年墨西哥米却肯州 级地震序列的触发关系及发震动力学探讨 分（）在主震东南约 处（，）发生 地震（图 ）。认识该地震序列的触发关系可以从应力转移角度理解地震发生的机理，此外该地震序列发生与当地应力场的关系也是地球科学所关心的问题。注：震源机制采用下半球施密特投影的海滩球表示，带有颜色的卦限为压缩区，红色为逆断型，蓝色为正断型，青色为正走滑型；图（）给出了研究区域的位置，其中红色线框表示研究区域的位置，黑色线条表示为板块边界。图 墨西哥米却肯地震序列周围的震源机制分布墨西哥位于北美板块，其地震活动主要受控于科科斯板块、里维拉板块以及北美板块的碰撞、挤压与俯冲结果，其中主要受科科斯板块俯冲北美板块的作用，该俯冲带形成的海底构造为墨西哥太平洋沿岸近海的中美海沟，因而在墨西哥西海岸线 范围内发生一系列强震，这些地震以逆冲型地震为主，震源深度大部分在 左右，破裂长度一般在 （张洪由，）。本次地震序列发生在北美板块与科科斯板块交界处，为板间型地震（，），附近有科科斯板块、太平洋板块、北美板块交汇分布（图 ），这种特殊的地理位置和复杂的地质构造使得该地区成为世界主要的地震高发区之一（，；孙克虬，），曾发生过多次超过 级的大地震，如 年 月 日墨西哥西海岸外太平洋底 级强震和 年 月 日墨西哥中部莫雷洛斯州 级地震，通常这些地震都发生在墨西哥西海岸的俯冲带上，并具有相似的发震机制（，）。前人对科科斯板块俯冲带做了大量的研究，等（）根据地震震源位置和震源机制研究区域的应力分布，从而确定科科斯板块俯冲的形状；等（）提出了一个全新的俯冲带几何模型，根据二维剖面详细描述了科科斯俯冲带三维几何模型的构造；等（）对科科斯板块的移动速度进行研究，结果表明板块以 的速度向北美板块俯冲移动；等（）对科 中国地震 卷科斯板块俯冲速度的不均匀性进行研究，发现这种不均匀性可能与奥罗斯科（）断裂带和特万特佩克（）岭的阻碍作用有关。但前人的研究工作主要针对该地区的俯冲带结构，而在构造应力场方面的研究却很少。为了解该地震序列的地震破裂机制，本研究根据多家机构提供的该地震序列的震源机制解分别计算出主震及两次较大余震的震源机制中心解，然后根据主震同震破裂模型计算在余震断层面和滑动方向上产生的库伦破裂应力变化。为进一步研究本次地震孕育和发生的动力学过程，收集了全球矩心矩张量（，简称 ）提供的历史上 次 地震的震源机制解数据（图 ），反演该地区的构造应力场。并且，研究了所确定的构造应力场投影到地震断层面上的相对剪应力和正应力，进而讨论该地震序列和应力场之间的关系。震源机制中心解的确定地震发生后，等国外机构很快获得了 地震的多个震源机制解。由于使用的方法和数据不同，这些震源机制具有相当大的离散度，这些结果都是震源错动方式的一种测量，因此可以根据这些测量结果给出中心值，用于对未来地震发生背景、地震应力触发、地壳应力场分析以及地震前应力方向改变的地震前兆进行研究。本研究首先搜集多家机构采用不同的方法和数据得到的该地震序列的震源机制解，为避免选择不同震源机制的困难，给出一个相对可靠的震源机制解，采用万永革（）提出的“同一地震多个震源机制中心解的确定”算法进行处理，震源机制中心解可以在多个震源机制中找到一个与其他震源机制空间旋转角平方和最小的震源机制，得出“最优解”，为进一步反演应力场提供了可靠的震源机制解。根据国外 、和 六个机构网站发布的主震震源机制结果整理得到表 。本研究分别以各个震源机制为初始解得到中心震源机制，给出大体一致（在小数点 位后有一定涨落）的标准差（表 第 列），表明采用这种方法得到的解是稳定的。尽管如此，本研究将各个机构测定的震源机制分别作为初始解，将得到标准差最小的解作为最终结果，发现将 得到的震源机制作为初始解得到的震源机制的标准差最小，本研究以此作为最终结果（表 第 行）。得到的中心震源机制和各个机构测定震源机制的最小空间旋转角见表 第 列并将所得到的中心震源机制及其不确定性绘于图 。从表和图中可以看出，该地震震源机制解距中心解的空间旋转角最大达 ，最小空间旋转角为 。这些数据表明不同机构得到的震源机制解较为集中。本文还收集了 和 等机构提供的 月 日 地震和 月 日 地震的震源机制解并分别计算其震源机制中心解。采用 、和 机构的震源机制解，得到 月 日 地震震源机制中心解见表 ，所得震源机制解距中心解的空间旋转角最大达 ，最小空间旋转角为 。采用了 、和 机构的震源机制解，得到 月 日 地震的震源机制中心解（表 ），所得震源机制解距中心解的空间旋转角最大达 ，最小空间旋转角为 。主震和 月 日 地震均为逆冲型地震，月 日 地震为正断型地震。期关兆萱等：年墨西哥米却肯州 级地震序列的触发关系及发震动力学探讨表 不同机构给出的墨西哥米却肯州 地震震源机制及作为初始解得到的震源机制中心解和标准差序号震源机制走向，倾角，滑动角（）机构作为初始解得到的震源机制解走向，倾角，滑动角（）作为初始解得到的标准差 （）以得到的中心解与其他震源机制的最小空间旋转角（），美国地质调查局 ，全球矩心矩张量组 ，法属波利尼西亚探测与地理物理实验室 ，德国地球科学研究中心 ，巴黎地球物理研究所 ，法国蔚蓝海岸台 ，表 墨西哥米却肯州地震序列震源机制中心解及 、轴走向和倾伏角震级震源机制中心解走向，倾角，滑动角（）轴走向（）轴倾伏角（）轴走向（）轴倾伏角（）轴走向（）轴倾伏角（），（）（）（）（）（）（），（）（）（）（）（）（），（）（）（）（）（）（）注：括号内容为置信范围。注：（）中的黑色弧线表示中心震源机制的两个节面，绿色弧线覆盖区域为其不确定范围；红色、蓝色和黄色的点表示中心震源机制解的 轴、轴和 轴，其周围对应颜色的封闭曲线表示其不确定性范围；绿点和黑点表示各个机构得到的震源机制的 轴和 轴的投影；紫色弧线表示各个机构得到的震源机制节面；（）中的压缩区域和膨胀区域分别用蓝色和红色表示（万永革等，）。图 墨西哥米却肯州 地震的中心震源机制解（）及空间三维辐射花样（）中国地震 卷 墨西哥 地震对后续地震的触发作用为讨论墨西哥米却肯州 月 日 地震对 月 日 和 月 日 地震是否存在触发关系，本文收集 提供的墨西哥 地震的断层破裂模型，该模型结果基于矩心矩张量节面（走向 ，倾角 ）计算，地震矩释放为 （），该破裂面长 ，宽 ，最大破裂量为 。地震产生的库伦破裂应力变化 表示为 （）其中，为视摩擦系数，取值范围在 之间，通用值 ；为地震产生的剪切应力变化量，与接收应力的断层滑动方向一致时取正值，当 与接收应力的断层滑动方向相反时取负值；为法向正应力变化量，拉张为正。这样就可以计算出此断层的库仑应力变化。选取前文所做的 和 地震的中心震源机制解中与俯冲带走向大体平行的节面（节面）作为接收断层面，根据 测定结果设置 地震的深度为 ，地震的深度为 。同时仿照前人的研究（，；，、），视摩擦系数取 ，使用 等（、）提出的 计算方法，当 为正值时，说明该次地震对研究区域断层的破裂具有促进作用；当 为负值时，说明该次地震对研究区域断层的破裂具有抑制作用。同时研究（，、；，）表明，库伦应力变化大于 （静态应力触发阈值）时会对后续地震的时空分布造成影响。按此方法计算出 月 日 地震在 日 地震断层面的滑动方向上产生的剪切应力变化为 ，正应力变化 ，库伦破裂应力变化为 ，超过 的阈值（图 （）。再计算出 日 地震断层面（节面）的滑动方向上产生的剪 切 应 力 变 化 ，正 应 力 变 化 ，库 伦 破 裂 应 力 为 ，超 过 的阈值（图 （）。墨西哥米却肯州 地震对 和 地震可能具有促进作用。主震产生的应力张量在 余震破裂面的滑动方向上产生的库伦破裂应力变化在主震破裂面附近均为正值，在断层面两端均为负值，而在 余震破裂面的滑动方向上产生的库伦破裂应力变化在主震破裂面附近均为负值，在破裂面两端均为正值，这是由于 为正断型地震，与 地震的破裂面几何形状大体接近，但滑动方向相反所致。这提醒我们，在计算库伦破裂应力变化研究地震之间的相互作用时，不仅要考虑两个地震位置的差异，还要考虑源地震（产生应力张量的地震）和接收断层形状和滑动方向的相互关系。地震序列周围应力场与地震序列的关系 应力场反演为进一步研究本次地震的孕震背景和发生条件，搜集了 提供的主震震源区及其邻 ：期关兆萱等：年墨西哥米却肯州 级地震序列的触发关系及发震动力学探讨图 墨西哥米却肯州 地震在 地震（）和 地震（）断层面和滑动方向上的库伦破裂应力变化表 墨西哥 地震附近 年的 次震源机制解序号日期（年 月 日）经度（）维度（）深度 走向（）倾角（）滑动角（）震级类型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 正走滑型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 逆断型 正走滑型 中国地震 卷区（，）范围内，年 月 日至 年 月 日共 次 地震的震源机制解数据（表 ）。根据震源机制类型划分标准（万永革，），逆断型地震有 个，占总体数量的 ；正断型地震有 个，占总体数量的 。本研究采用 等（）提出的网格搜索法对应力场进行反演，网格搜索法是在应力场参数的模型空间中找到应力模型与实际地震数据间平均残差最小的最佳应力模型。相较于其他应力场反演方法，网格搜索法能够设置更加精细化的应力场格值参数来搜索构造应力场参数的最优解（李枭等，），从而得到更为准确的结果，并给出求解参数在一定置信度下的置信范围。采用上述方法进行求解得到应力场的最优解参数为主压应力轴 的走向为 ，倾伏角为 ，不确定范围分别为 和 ；中间轴 的走向为 ，倾伏角为 ，不确定范围分别 和 ；主张力轴的走向为 ，倾伏角为 ，不确定范围为 和 。应力比 ，表明该区域主要以挤压应力为主，中间轴与张应力轴表现为拉张作用，张应力轴表现为接近垂直的状态（图 ）。反演的构造应力场结果表明，在科科斯板块 向俯冲作用下，该地呈现较强的 向挤压应力，俯冲过程导致了墨西哥沿岸附近区域隆升，造就了山脉隆起。注：（）应力场反演结果及数据拟合情况的下半球施密特投影图，图中蓝色大箭头表示主张力轴的水平最佳方向，红色大箭头表示为主压力轴的水平最优方向；蓝色小箭头表示为“可能断层面”的观测滑动方向，红色小箭头表示为“可能断层面”的理论滑动方向；绿色断层面表示置信水平为 的应力场最大剪应力平面；黄色箭头表示该平面的滑动方向；黑色弧线表示为“可能断层面”的施密特投影；（）墨西哥地震地区应力状态三维示意图，红色为压缩区，蓝色为膨胀区。图 墨西哥米却肯地震周围应力场反演结果（）及三维示意图（）应力场与地震序列的关系万永革（）对挤压应力体系和拉张应力体系下所产生震源机制及其剪应力和正应力的表现进行了模拟。为分析本文反演的应力场与本次地震序列的关系，基于该方法，得到墨西哥西海岸地区构造应力场在不同几何形态的断层上产生的相对剪应力、正应力（采用最大剪应力（最大主应力与最小主应力之差的一半）进行归一化）和震源机制类型，并将其作用于墨西哥 地震及两次余震的震源机制中心解的断层面上产生的相对剪应力与正应力，如图 所示。期关兆萱等：年墨西哥米却肯州 级地震序列的触发关系及发震动力学探讨图 墨西哥米却肯州西海岸构造应力场作用下不同几何形态断层面上的相对剪应力（）和相对正应力（）墨西哥米却肯 地震震源机制中心解节面 相对剪应力和相对正应力分别为 （最大为 ）和 ，剪应力滑动角为 ，与观测滑动角相差 ；节面相对剪应力和相对正应力分别为 和 ，剪应力滑动角为 ，与观测滑动角相差 。墨西哥 地震震源机制中心解节面 相对剪应力和相对正应力分别为 和 ，剪应力滑动角为 ，与观测滑动角相差 ；节面相对剪应力和相对正应力分别为 和 ，剪应力滑动角为 ，与观测滑动角相差 。墨西哥 地震震源机制中心解节面 相对剪应力和相对正应力分别为 和 ，剪应力滑动角为 ，与观测滑动角相差 ；节面 相对剪应力和相对正应力分别为 和 中国地震 卷 ，剪应力滑动角为 ，与观测滑动角相差 。很明显，这三次地震的震源机制中心解的节面基本都处于剪应力较大区域，地震和 地震的剪应力滑动角与对应节面的观测滑动角比较接近，但 地震的剪应力滑动角与对应节面的观测滑动角相差很大（几乎呈反方向滑动），而这三个地震的断层面几何形态基本相似（地震断层面与 地震断层面的法向夹角为 、地震断层面与 地震断层面的法向夹角为 ，图 中这三个地震相距很近也说明了这一问题），这表明 地震不同于 和 地震。俯冲带发生正断型地震的可能原因一直处于探索之中，如 等（）提出墨西哥中部之下的俯冲科科斯板片导致的热应力是一种可能的原因，然而如果按照这种解释，正断型地震应该在一个较大区域是集中的。墨西哥俯冲带是地震多发区域，但没有看到统计学意义的正断型地震较为集中区域。根据 和 地震的断层面和发震时间的相近性及断层滑动几乎完全相反的情况，推测该正断型地震的发震原因可能是：科科斯板块向北美板块俯冲，因而北美板块边界形成了巨大的压力导致了断层突然滑动，引发 逆冲型地震，但由于突然释放的逆冲分量过大，因此需要一定程度上的正断分量来缓解过多的逆冲分量。另外根据地震的发生时间，主震发生 后发生了 地震，而又间隔了 才发生了 地震，这就表明主震发生很短时间就对应力状态进行了调整，而后续间隔了较长时间才发生同种性质的地震。此外，地震距主震破裂非常近，或者就发生在破裂面上，且断层面的几何形状非常接近，根据这些情况推测，地震的正断层震源机制可能是由于主震逆冲过度造成的。本次地震对该区域的整体地球动力学研究没有影响，只是短暂时间、非常局部、震级较小的一次扰动事件。在墨西哥米却肯州西海岸地区构造应力条件下，三个地震的震源机制中心解的节面震前基本都处于最大剪应力加载状态，但 地震的滑动方向与加载应力方向完全相反。由此可见，和 地震都是在科科斯板块 向俯冲产生逆断型构造应力背景作用下，基本沿应力场的剪应力最大的断层面进行能量释放（冯淦等，），而 地震是一次与背景构造应力场差别较大的地震事件。结论本文结合不同国外机构给出的墨西哥米却肯地震序列的震源机制数据结果，得到 主震及 和 余震的震源机制中心解，计算主震产生的库伦破裂力对两次较大余震的触发作用。通过反演墨西哥米却肯地区的构造应力场并研究本次地震发生的地球动力学背景，主要得出以下结论：（）地震震源机制中心解的结果为：节面 走向 ，倾角 ，滑动角 ；节面走向 ，倾角 ，滑动角 。地震震源机制中心解的结果为：节面 走向 ，倾角 ，滑动角 ；节面走向 ，倾角 ，滑动角 。地震震源机制中心解的结果为：节面 走向 ，倾角 ，滑动角 ；节面走向 ，倾角 ，滑动角 。结果表明 和 地震为逆冲型地震，地震为正断型地震。期关兆萱等：年墨西哥米却肯州 级地震序列的触发关系及发震动力学探讨（）利用 提供的破裂模型探讨这三次地震的关系，发现 地震对 和 地震 可 能 具 有 促 进 作 用，库 伦 破 裂 应 力 分 别 为 和 ，均 超 过 的阈值，该结果为研究本次地震对震源周围的影响及以后地震活动趋势提供了基础资料。（）应力场反演结果显示墨西哥米却肯地震序列震源区总体上受到 向和 向的挤压力、垂直向的拉张力，整体表现为逆冲型应力类型，与在科科斯板块 向的俯冲作用下墨西哥西海岸沿岸发生了一系列逆冲型地震的地震活动背景相符合。在墨西哥米却肯州西海岸地区构造应力条件下，三次地震中，每个地震的两个节面震前基本都处于最大剪应力加载状态，由此可见，此次墨西哥米却肯地震序列是在科科斯板块 向俯冲产生逆冲型构造应力场作用下，经过一定时间的应力积累，沿应力场最大剪应力平面进行的正常能量释放。墨西哥地震活动基本上由于科科斯板块俯冲导致，这与中国台湾受到太平洋板块和菲律宾板块的俯冲有相似特点，本研究可为中国台湾地震活动分析提供借鉴。致谢：本研究所用的破裂模型来自 ，绘图采用 软件（，）和 软件绘制，特此致谢。参考文献冯淦，万永革，许鑫，等 年青海玛多 地震对周围地区的应力影响 地球物理学报，（）：李枭，万永革，许鑫，等 地震断层带破碎非均匀程度研究 以新马德里地震带 断层为例 地球物理学报，（）：孙克虬 墨西哥地震与板块运动 海洋地质动态，（）：万永革，盛书中，许雅儒，等 不同应力状态和摩擦系数对综合 波辐射花样影响的模拟研究 地球物理学报，（）：万永革 同一地震多个震源机制中心解的确定 地球物理学报，（）：万永革 震源机制与应力体系关系模拟研究 地球物理学报，（）：万永革（）震源机制水平应变花面应变的地震震源机制分类方法及序列震源机制总体特征分析 地球科学，：张洪由 墨西哥的地震、火山监测及防灾减灾工作 国际地震动态，（）：，：，（）：，：，（）：，：，（）：，（）：，：，（）：，（）：，（）：中国地震 卷 ，：，（）：，（）：，：，（）：，（）：，（）：，（）：，：，（）：，（）：，（）：期关兆萱等：年墨西哥米却肯州 级地震序列的触发关系及发震动力学探讨 ），），），），），），（），：，：；