九寨沟MS 7.0地震静态库仑应力触发及对附近断层影响分析

朱航 屈勇

四川省地震局，成都市人民南路3段29号 610041

0 引言

2017年8月8日九寨沟MS7.0地震发生在巴颜喀拉块体东边界。根据先前的认识，巴颜喀拉块体位于青藏活动地块区中北部，其北、南边界分别被具有左旋走滑性质的东昆仑断裂带、甘孜-玉树断裂带和鲜水河断裂带所包围，东边界则由龙门山断裂带、岷江-虎牙断裂带等共同构成（张培震等，2003）。自1997年西藏玛尼MS7.9地震发生以来，巴颜喀拉块体边界共发生8次7级以上地震（含2次8级地震），其中，3次发生在四川境内，即2008年汶川MS8.0、2013年芦山MS7.0和2017年九寨沟MS7.0等地震；这3次地震都发生在巴颜喀拉块体东边界。此外，巴颜喀拉块体东边界1973～1976年曾发生4次MS≥6.5地震，即1973年8月11日松潘黄龙MS6.5以及1976年8月 16～23日松潘MS7.2、MS6.7、MS7.2地震。上述7次地震共同构成了1973～2017年巴颜喀拉块体东边界强震序列。

为了了解巴颜喀拉块体东边界强震序列的发生规律和相关物理机制，我们曾对其应力触发过程进行了系列跟踪研究。朱航等（2009）计算和研究了1973～1976年四川松潘4次强震组成的序列引起的静态库仑应力变化图像，分析了由该序列各次事件引起的近场应力变化及其与后续强震间的关系，同时分析了该序列引起的远场应力变化与随后25年区域中、强地震活动间的关系。结果显示，1973年8月11日松潘黄龙6.5级地震导致虎牙断裂带中段库仑应力显著增加，并触发了1976年8月16日的7.2级地震；此后，又沿断层向南相继触发了1976年8月22日6.7级地震和8月23日7.2级地震。此外，在该强震序列发生后的25年间，在距该序列发震断层中部约200km范围内，6次5.0～6.6级地震均发生在由该序列引起的远场、微量的库仑应力变化增加区中。汶川8.0级地震发生后，Zhu等（2010）发现，2008年汶川8.0级地震的震中位于1973～1976年四川松潘强震序列引起的库仑应力触发区内，说明松潘强震序列的发生导致了龙门山断裂带中、南段应力水平的增强，从而可能在一定程度上使得汶川8.0级地震提前发生；计算松潘强震序列和汶川8.0级强震共同引起的库仑应力变化，结果显示，汶川8.0级地震触发了2008年5月25日龙门山断裂带北东段上的青川6.4级地震；龙门山断裂带南段也处于应力触发区内；与汶川地震发震断层平行的茂县-汶川断裂（龙门山断裂带后山断裂）则处于应力影区内，发生强震的危险性应相对减小。值得说明的是，2013年芦山7.0级地震发生在该研究提出的龙门山断裂带南段应力触发区内。屈勇等（2017）计算和研究了1973～2014年发生在巴颜喀拉块体东-南边界断裂带上的、由9次MS≥6.3地震组成的强震序列引起的库仑应力变化图像，分析序列中先发地震破裂对后发地震的应力触发作用，结果表明：①1973～1976年在巴颜喀拉块体东边界的虎牙断裂带上发生的4次强震存在显著的应力触发关系，同时改变了龙门山断裂带中南段及周围区域的应力状态，相继触发了汶川8.0级、芦山7.0级地震；②南边界鲜水河断裂带自1973年以来发生的3次强震同样存在应力触发关系，2014年康定6.3级地震是巴颜喀拉块体东、南边界强震活动导致的库仑应力变化共同触发的结果。

2017年九寨沟7.0级地震发生在1973年松潘黄龙6.5级地震震中以北的邻近区域，它与先前发生在巴颜喀拉块体东边界的各次强震间的应力触发关系，以及它向邻近断层的应力转移情况，是本文所关注的问题。本文将计算有关的静态库仑应力变化，探讨九寨沟7.0级地震的应力触发机制，并指出由九寨沟7.0级地震发生所引起的库仑应力增加而可能被触发强震的断层。

1 研究方法

库仑破裂应力是指断层面上存在着一定程度的应力积累，当积累应力大小超过断层面的滑动阻抗强度时，断层面将会滑动和破裂。关于库仑破裂准则以及地震的应力触发理论已有大量论述（King et al，1994；Stein，1999；Pollitz et al，2002；Deng et al，1997；Harris，1998；Okada，1992），库仑破裂准则确定了岩体中先存断层重新发生摩擦滑动或者“破裂”的应力条件（Scholz，1990；Jaeger et al，1979）。断层上的库仑应力若发生改变，将影响该断层重新滑动或破裂的失稳条件；超过0.01MPa的微小库仑应力变化就会影响断层的应力分布状态，从而增加或缩短断层发生潜在破裂所需的时间（King et al，1994；Stein，1999）。研究断层趋于破裂的过程可使用库仑应力变化Δσf进行数值描述

式中，Δτs为断层面的剪切应力变化（沿滑动方向为正）；Δσn为正应力变化（张应力为正），μ′为断层面的摩擦系数。

当一次强震发生时，发震断层的应力将在一定程度上得以释放，同时，地震时的同震位错会引起附近区域断层的静态应力变化，即应力转移。当所引起的相邻断层段上的库仑应力变化为正值时，该断层段重新发生粘滑的可能性增大，或者会很快发生另一次地震，这种情况即地震的应力触发作用；当库仑应力变化为负值时，该断层段重新发生粘滑的可能性降低。一般称库仑应力变化的正值区域为应力触发区，负值区域为应力影区（King et al，1994）。引起应力变化的地震和破裂面称为“源地震”和“源断层”，区域中的其它先存断层称为“接受断层”。根据Okada（1992）提出的均匀各向同性弹性半空间模型，可由源断层的同震位错计算接受断层的库仑应力变化量。本文采用King等（1994）推导的库仑应力变化公式以及发展的相关算法，并使用Coulomb3.3软件分别计算1973～2017年间巴颜喀拉块体东边界强震序列引起的静态库仑应力变化，分析可能存在的触发作用，并进一步研究该强震序列对周边地区中、强地震活动的影响。

根据以上库仑应力原理可知，由Coulomb3.3软件计算的接受断层上的静态库仑应力变化量取决于源断层的几何特征（断层走向、倾角、滑动角）与源地震的同震位错，以及接受断层的几何特征及断层摩擦系数等（King et al，1994）。

2 发震构造与位错模型

受印度洋板块北向推挤作用的影响，青藏活动地块区整体朝S-SE方向运动，位于其中的巴颜喀拉块体受华南地块阻挡，在东边界带上产生了强烈的构造运动和变形。巴颜喀拉块体北边界主要由东昆仑断裂带构成，南边界的南东段主要由鲜水河、甘孜-玉树断裂带构成，均表现为左旋走滑运动。本文以巴颜喀拉块体东边界及邻近地区为主体，取29.5°～35.0°N、101.5°～106.5°E的矩形范围作为研究区。区内1973～2017年强震序列分别发生在巴颜喀拉块体东边界的NW向树正断裂、NNW向虎牙断裂带以及NE向龙门山断裂带上，其中，树正断裂表现为左旋走滑，而虎牙断裂带则表现为左旋走滑兼逆冲的运动方式，虎牙断裂带向南延伸逐渐过渡到由3条主干断裂组成的龙门山逆冲推覆断裂带（图1）。

图1显示了1973～2017年巴颜喀拉块体东边界强震序列的震中位置及发震构造，并标注了各次地震的震源机制解。由图1可见，2008年汶川MS8.0、2013年芦山MS7.0等地震均发生在龙门山断裂带中、南段，震源机制为逆冲型；1976年8月16～23日2次松潘MS7.2地震发生在NNW向的虎牙断裂带，震源机制为左旋走滑兼逆冲型；1976年8月22日松潘MS6.7地震发生在NE向叶塘断裂，该断裂为逆冲型。2017年九寨沟MS7.0地震发生在NW向树正断裂，震源机制为左旋走滑型。

关于1973年松潘黄龙MS6.5地震的发震构造尚存疑问。据《中国震例（1966－1975）》（程式等，1988）记载：“1973年8月11日松潘黄龙MS6.5地震震中附近没有发现明显的断裂构造，因此，使用等震线走向、余震分布和震源机制解节面的走向判断，这次地震的构造位置位于NNW向虎牙断裂向北的延伸线上。”由图1可见，1973年松潘黄龙MS6.5地震与2017年九寨沟MS7.0地震的震中都位于树正断裂的西侧，它们的震源机制非常相似，均为走滑型，其NW向节面走向与树正断裂的走向一致。据此，我们可以认定，这2次地震都发生在NW走向的树正断裂，均为左旋剪切破裂。

为了获取计算所需的各次地震事件的震中位置、发震断层产状及同震位错量，我们在综合多种资料的基础上，构建利用Coulomb3.3软件计算所需的断层模型。表1为巴颜喀拉块体东边界强震序列7次事件发震断层产状及同震位错参数，其数据来源见表注。表1中1973年松潘黄龙MS6.5、2017年8月8日九寨沟MS7.0等地震的同震位错量由以下闻学泽（1995）建立的走滑型地震断层的经验公式估计得出

图1 巴颜喀拉块体东边界活动断裂及1973年以来MS≥6.5地震分布

表1 巴颜喀拉块体东边界强震序列7次事件发震断层产状及同震位错参数

式中，MS为面波震级；u×l为同震平均位错u与破裂长度l的乘积。

3 库仑应力变化图像与分析

计算之前，将研究区内1973年以前的库仑应力变化量归零。计算中，地壳的剪切模量取3.3×1010Pa，泊松比取0.25，断层面摩擦系数μ′取 0.4；该取值是 King等（1994）根据不同地区计算的库仑应力变化得出的经验值。计算深度为地下10km。

图2 1973～2013年发生的6次MS≥6.5地震共同引起的以树正断裂为接受断层的库仑应力变化图像

3.1 巴颜喀拉块体东边界先发地震对九寨沟MS 7.0地震的触发作用

图2是巴颜喀拉块体东边界1973～2013年发生的6次MS≥6.5地震共同引起的以树正断裂为接受断层的库仑破裂应力变化图像。从表1给出的九寨沟地震发震断层产状可知，断层运动为左旋走滑，断层走向NW，高倾角，断层面接近直立。由图2可见，在巴颜喀拉块体东边界大部分段落已经破裂的情况下，树正断裂北段的库仑应力水平显著升高，而九寨沟MS7.0地震即发生在应力触发区内，因此认为，巴颜喀拉块体东边界先发地震对九寨沟地震具有明显的触发作用。

3.2 巴颜喀拉块体东边界强震序列对邻近断裂的应力转移

九寨沟MS7.0地震是巴颜喀拉块体边界断裂带自1997年以来发生的第8次MS≥7.0地震，至此，巴颜喀拉块体东、西边界大部分已经破裂。那么，巴颜喀拉块体东边界强震序列的发生对未来强震形势有何影响？

九寨沟MS7.0地震震中位于巴颜喀拉块体东北角的树正断裂北段，其南面的巴颜喀拉块体东边界断层段多已破裂，积累的应变能得以释放，因此，不具备未来10年乃至数十年内强震复发的条件。本文主要关注九寨沟地震震中的西、北面3条断裂带的近场库仑应力变化情况，具体如下。

3.2.1 巴颜喀拉块体东边界强震序列引起的岷江断裂带库仑应力变化

岷江断裂带南起松潘以南，向北与东昆仑断裂带东段相交，大致沿岷江河谷总体呈SN向延伸，长约120km。岷江断裂带曾发生过2次较强地震，分别是1748年漳腊北面的6.5级、1938年松潘6.0级地震。钱洪等（1995）认为，岷江断裂带活动强度不大，至少没有发生过造成断错地貌的强烈破坏性地震。

岷江断裂带以逆冲为主，兼有一定左旋运动分量，断层倾向W，倾角约80°（邓起东等，1994）。图3（a）是九寨沟地震发生后，巴颜喀拉块体东边界1973～2017年发生的7次MS≥6.5地震共同引起的以岷江断裂带为接受断层的库仑破裂应力变化图像。由图3（a）可见，岷江断裂带处于应力影区范围内，即邻近的树正断裂破裂错动，对岷江断裂带而言，具有减震效应。

3.2.2 巴颜喀拉块体东边界强震序列引起的东昆仑断裂带库仑应力变化

东昆仑断裂带是巴颜喀拉块体北边界断裂系统的一条巨型左旋走滑活动断裂，自有仪器记录的1902年以来，沿该断裂带先后发生了7次MS6.9～8.1的大震，震级最大的为2001年11月14日在西段发生的昆仑山口西8.1级地震；而在断裂带东段上存在一处地震空区——玛曲空区（付俊东等，2012），因此，该段具有发生强震的潜在危险。

东昆仑断裂带东段呈帚状散开的几何特征。总体来看，东昆仑断裂带东段以左旋走滑运动为主，兼有一定逆冲分量，断裂大体走向 300°，高倾角，断层面倾向 SW（李陈侠等，2009）。图3（b）是九寨沟地震发生后，巴颜喀拉块体东边界1973～2017年发生的7次MS≥6.5地震共同引起的以东昆仑断裂带东段为接受断层的库仑应力变化图像。由图3（b）可见，该断裂处于应力触发区范围内，即东边界强震序列的发生，对东昆仑断裂带东段，特别是对玛曲空区具有强震触发作用。

3.2.3 巴颜喀拉块体东边界强震序列引起的龙日坝断裂带库仑应力变化

龙日坝断裂带位于龙门山断裂带NW侧约200km的巴颜喀拉块体东部，其北东段由走向 54°N±5°E、相距约30km的龙日曲、毛尔盖2条分支断层组成，分别表现出向SE方向逆冲和右旋走滑的运动方式。以龙日坝断裂带为界，巴颜喀拉地块分为西部阿坝和东部龙门山2个次级块体；龙门山次级块体的整体缩短和隆升反映出从龙门山断裂带到龙日坝断裂带是巴颜喀拉地块SE向运移过程中由于受到华南地块的强烈阻挡而形成的后展式推覆构造系统，龙日坝断裂带正是这一系统中晚第四纪新生的活动断裂带（徐锡伟等，2008）。

图3 1973～2017年发生的7次MS≥6.5地震共同引起的库仑应力变化图像

图3（c）是九寨沟地震发生后，巴颜喀拉块体东边界1973～2017年发生的7次MS≥6.5地震共同引起的以龙日坝断裂带为接受断层的库仑破裂应力变化图像。由图3（c）可见，该断裂带处于应力触发区范围内，即东边界强震序列的发生有可能触发龙日坝断裂破裂。

4 结论与讨论

（1）2017年8月8日九寨沟7.0级地震是巴颜喀拉块体东边界自1973年以来发生的第7次MS≥6.5地震，利用库仑应力计算结果分析表明，先前发生的6次地震对九寨沟地震具有触发作用。

（2）对九寨沟7.0级地震震中附近主要断裂应力转移的计算结果显示，在巴颜喀拉块体东边界7次MS≥6.5地震的共同作用下，东昆仑断裂带东段及龙日坝断裂带库仑应力水平升高，这可能导致这些断裂已积累的应变能提前释放，即具有触发作用；岷江断裂带处于库仑应力影区内，具有减震效应。

（3）发生在巴颜喀拉块体东边界的强震是否会短期内触发北、南边界的地震，这是我们所关注的问题。2014年11月22日南边界的南端发生了康定6.3级地震，其与2013年4月20日芦山7.0级地震的时间间隔仅有1.6年，研究认为，这是由于巴颜喀拉块体东、南边界强震活动导致的库仑应力变化共同触发的结果（屈勇等，2017）。九寨沟7.0级地震发生在巴颜喀拉块体东边界的北端，其相邻的北边界东昆仑断裂带东段目前不仅位于应力增加区内，而且还存在玛曲地震空区，所以，应特别关注该段未来地震危险性。