引言

前人对于动态和静态应力变化的研究表明， 同震过程导致了断层的相互作用和地震的触发． 大震后区域静态库仑应力变化常常被用于解释区域地震活动性速率的变化(Reasenberg,Simpson, 1992; Simpson,Reasenberg, 1994; Harris,Simpson, 1998; Toda et al, 1998)、 主震断层外余震的发生(Dieterich, 1972, 1994; Rybicki, 1973; Das, Scholz, 1981; Stein et al, 1992, 1994; Dieterich,Kilgore, 1996; Gross,Kisslinger, 1997; Gross, Bürgmann, 1998)以及即将失稳断层的地震发生概率的变化(Toda et al, 1998)； 而动态模型则被用于研究由动态应力变化引起的断层间相互作用和地震的触发(Harris,Day, 1993; Gomberg et al, 1997, 1998)． 1992年美国Landers地震为我们提供了一些小范围内断层相互作用和远距离地震活动性触发的实例． King等(1994)计算得到Landers地震对后续Big Bear M_S6.5地震产生的库仑破裂应力增量约为0.2 MPa， 认为是Landers地震加速了Big Bear地震的发生． Spudich等(1995)还讨论了关于Big Bear余震是由1992年Landers地震触发的问题， 并指出无论用动态触发还是静态触发，都很难解释这3.5个小时的时间延迟．Gomberg等(1998)解释了摩擦失稳模型中近距离(静态和动态)触发和远距离(动态)触发的主要特征， 并讨论了静态和动态加载是怎样改变了断层的失稳时间．Harris(1998) 则详细地总结了应力触发、 应力影区及其对地震危险性的影响．

关于汶川地震应力触发的相关研究受到了国内外专家学者的重视． 汶川地震震中位于龙门山断裂的中段， 震源深度14 km(中国地震台网中心， 2013a)， 破裂沿着北东走滑、 西部倾斜的龙门山逆冲断层带向前扩展． 该地震引起地壳深部的岩石破裂长达300多千米(陈运泰, 2008; 王卫民等, 2008)． Toda等(2008)利用几条主要活动断裂插值得到了汶川地震产生的应力变化分布． 万永革等(2009)和单斌等(2009)计算了汶川地震导致的周围断层应力变化， 解朝娣等(2010)在其基础上探讨了汶川地震引起的周边断裂应力变化后地震活动性随时间变化的情况． Parsons等(2008)还初步探讨了汶川地震产生的应力变化和对周围断层的影响．

汶川地震发生近5年后， 2013年4月20日8点2分46秒， 我国四川省雅安市芦山县(30.3°N， 103.0°E)发生了M_S7.0地震， 震区发生大面积的滑坡和建筑物损毁． 地震定位和震源机制结果显示该地震是发生在龙门山断裂带上的又一次逆冲型地震， 也是汶川地震后发生在龙门山断裂带上最强的一次地震． 芦山M_S7.0地震发生在龙门山断裂带的南端， 在汶川地震发生时这一地段没有发生破裂滑动． “芦山地震是否为汶川地震的余震”在学界引起了一定的争论． 陈运泰等(2013)认为芦山地震是汶川地震的最强余震； 王卫民等(2013)指出“芦山地震在宏观上可视为汶川地震一次‘迟到’的强余震”； 而刘杰等(2013)认为两次地震的余震区存在约45 km的间隔， 芦山M_S7.0地震不是汶川地震的余震．

然而， 无论持有哪种观点， 各位专家都不否认芦山M_S7.0与汶川M_S8.0地震有密切的联系． 汶川地震发生后， 雅安芦山地区的库仑应力增加的事实是肯定的(Parsons et al, 2008)． 那么汶川地震究竟是怎样触发芦山地震？ 是否还会有类似芦山地震的强余震在周围地区发生？ 这才是研究汶川地震后区域库仑应力变化的最主要目的和最核心的问题．

本文利用含(滑移)速率和状态的摩擦定律(Dieterich, 1994)， 结合汶川地震前后的地震活动性水平， 深入阐述了芦山M_S7.0地震的静态应力触发过程， 从地震发生概率的变化上理解芦山地震的发震根源， 并计算周边的鲜水河断层和熊坡断层的发震概率变化， 从而为区域地震活动性以及危险性评价提供重要的参考．

1.   含(滑移)速率和状态的摩擦定律

含(滑移)速率和状态的摩擦本构关系为我们提供了一个对断层属性进行复杂的定量实验观测的基本框架(Dieterich, 1979， 1981， 1994; Ruina, 1983)． 最简单地表述含(滑移)速率和状态的断层应力是由Ruina, 1983)基于Dieterich(1994))的理论提出的． 如果将其推广为多态变量(Rice,Gu, 1983)， 则可以写作

This page contains the following errors:

Below is a rendering of the page up to the first error.

图  1  弹簧-滑块模型示意图， 由此得到的滑块运动方程为：τ =-k(δ-δ₀).其中τ为剪切摩擦应力， k为弹簧的有效弹性系数，δ和δ₀分别为滑动位移和初始滑动位移

Figure  1.  Spring-slide block model. The equation of motion for the slider is τ =-k(δ-δ₀)， where k， δ， δ₀，τ are the spring constant， the displacement of the block，the displacement of the right-hand end of the spring， and the shear friction, respectively

下载: 全尺寸图片 幻灯片

This page contains the following errors:

Below is a rendering of the page up to the first error.

进一步， 假设第n次子事件的失稳时间t_f=nΔt， 其中Δt是第n-1次事件到第n次事件的时间间隔， 那么可以得到第n次子事件的滑动速率

将式(4)代入式(5)， 可得地震(事件)次数

This page contains the following errors:

Below is a rendering of the page up to the first error.

This page contains the following errors:

Below is a rendering of the page up to the first error.

由于地震活动性与区域应力场有关， 那么区域应力场的改变必将导致相应的地震危险性概率的变化． 假定地震序列在时间t内， 对震级M以上的发震概率服从泊松模型(Dieterich,Kilgore， 1996)

那么， 利用式(8)的概率分布， 结合式(7)可以计算在应力Δτ的扰动下造成大于M级地震的发震概率．

2.   发震概率计算

截至2013年2月28日， 汶川地震已有M>5.0的余震128次， M>4.0的余震841次， M>2.0的余震2 174次(中国地震台网中心, 2013b)， 而且余震还在持续． 图2显示了汶川地震前(1970-01-01—2008-05-11)后(2008-05-12—2010-05-11)的地震分布及M>2.0地震的年发生率． 震后汶川周边地区的应力状态发生了明显地变化(Parsons et al, 2008; Toda et al, 2008; 万永革等, 2009; 单斌等, 2009; 解朝娣等, 2010)， 从而导致了地震活动性分布的显著不同．

图  2  汶川地震前(a)、 后(b)的M>2.0地震分布及地震年发生率网格大小为0.05°×0.05°， 右侧色标显示的年发生率为网格内M>2.0的年发生率. 红色圆圈代表M>5.0地震， 黑色圆圈代表2.0＜M＜5.0地震， 灰色实线代表断层位置

Figure  2.  Distribution of M>2.0 earthquakes and seismicity rate of Longmenshan fault before (a) and after (b) Wenchuan earthquake Red and black circles represent M>5.0 and 2.0＜M＜5.0 earthquakes, respectively.The seismicity on colorbar on the right is in the grid of 0.05°×0.05°. The gray line is the fault

下载: 全尺寸图片 幻灯片

假设汶川地区的余震持续时间为85年(申文豪等, 2013)， 大于设定震级M的地震在该地区的年发生率为0.01次， A=0.01， σ_n=10 MPa (Dieterich，Kilgore， 1996)． 在不同应力扰动的情况下， 根据式(6)和(7)计算所得发生大于设定震级M的概率如图3所示． 从图3可见， 大应力扰动可导致地震发生率的快速上升， 因此在受到大应力扰动的地区地震发生的时间会大幅提前．

图  3  不同应力扰动下发生大于设定震级M的概率变化

Figure  3.  The probability of earthquake occurrence of magnitude >M on different shear stress disturbance

下载: 全尺寸图片 幻灯片

对于雅安地区(30.0°—30.7°N， 102.5°—103.5°E)， 汶川地震后库仑应力增加了0.1 MPa (Parsons et al, 2008)． 从1970年1月1日—2008年1月1日， 共发生了23次M>4.0地震， 其中仅发生一次M>6.0的地震． 那么根据背景场地震活动性速率公式r_≥M=N/Δt(其中N为发生大于震级M的地震次数， Δt为时间间隔)可以求出r_≥6.0=0.026 3次/年． 依然假设A=0.01， σ_n=10 MPa (Dieterich,Kilgore， 1996)， 那么在汶川地震后雅安地区的地震活动性概率变化如图4所示． 从图4可以看出， 发震断层的余震持续 时间对于发震

图  4  雅安地区发生M>6.0地震的概率变化不同颜色曲线代表了不同的余震持续时间, 黑色箭头所指位置为现今雅安地区对应的发震概率

Figure  4.  The probability of earthquake occurrence of M>6.0 in Ya'an area Different colors give difference aftershock durations, the black arrow shows the earthquake probability in Ya'an area now

下载: 全尺寸图片 幻灯片

概率的影响并不大， 而受到扰动之后的几年内发震概率的增幅较快． 芦山地震发生在汶川地震后近5年， 而由于受汶川地震的影响， 雅安地区此时发生M6.0地震的概率已经增加到25%—30%．

事实上， 汶川地震之后不仅雅安地区的库仑应力增加， 其周边的鲜水河断层和熊坡断层也存在库仑应力的明显增加(Parsons et al, 2008), 这些断层仍然是今后需要防范的重点地区． 图5显示了汶川地震后鲜水河断层不同震级的发震概率变化． 根据中国地震台网中心给出的地震目录， 从汶川地震后到2013年2月28日， 鲜水河断层已发生M>4.0地震10次， 其中2010年4月27日发生一次M5.0地震， 此时M>5.0的发震概率已超过80%， 而至今未发生M>6.0的地震， 而M>6.0的发震概率已约为60%， 所以鲜水河断层附近将成为M>6.0地震的重点防范地区. 鲜水河断层未来发生M>6.0地震概率由表1给出．

图  5  汶川地震后鲜水河断层发生地震的概率变化实箭头所指位置为现今鲜水河断层对应的M>6.0的发震概率， 虚箭头所指的位置为2010年4月27日发生的M5.0地震对应的发震概率

Figure  5.  The probability of earthquake occurrence of Xianshuihe fault after Wenchuan earthquake The black arrow shows the probabilities of earthquake occurrence of M>6.0 on Xianshuihe fault at present, the dashed arrow points out the probabilities of earthquake occurrence when the M5.0 earthquake happened on April 27, 2010

下载: 全尺寸图片 幻灯片

表  1  鲜水河断层未来发生M>6.0的地震概率

Table  1.  The probabilities of earthquake occurrence of M>6.0 on Xianshuihe fault

下载: 导出CSV 

| 显示表格

图6显示了汶川地震后熊坡断层不同震级的发震概率变化． 根据中国地震台网给出的地震目录， 从汶川地震后到2013年2月28日， 尚未发生M>4.0地震， 此时M>4.0的发震概率已接近90%，虽然M>5.0的发震概率仅为26%， 但是该地区将来仍旧有发生中强地震的危险性. 熊坡断层未来发生M>4.0地震概率由表2给出．

图  6  汶川地震后熊坡断层发生地震的概率变化黑色虚线所指的位置为现今对应的发震概率

Figure  6.  The probabilities of earthquake occurrence of Xiongpo fault after Wenchuan earthquake The dashed line shows the probabilities of earthquake occurrence at present

下载: 全尺寸图片 幻灯片

表  2  熊坡断层未来发生M＞4.0地震概率

Table  2.  The probabilities of earthquake occurrence of M＞4.0 on Xiongpo fault

下载: 导出CSV 

| 显示表格

3.   讨论与结论

芦山M_S7.0地震的发生重新引起了对汶川地震的热议． 本文利用含(滑移)速率和状态的摩擦定律(Dieterich, 1994)， 结合汶川地震前后的地震活动性水平， 从地震发生概率 的变化上解释了芦山地震的发震根源; 计算得到了芦山地震发生时， 雅安地区发生M>6.0 地震的概率已经增加了25%—30%． 此外， 本文对于库仑应力增加的鲜水河断层和熊坡断层也进行了发震概率的计算． 根据中国地震台网中心给出的地震目录， 从M>4.0地震的发生情况来看， 计算结果与实际情况基本符合． 鲜水河断层从汶川地震后至今未发生M>6.0地震， 而M>6.0的发震概率已约为60%; 熊坡断层自汶川地震以来尚未发生M>4.0地震， 此时M>4.0的发震概率已接近90%． 所以， 我们认为鲜水河断层附近将成为M>6.0地震的重点防范地区， 熊坡地区将来仍旧有发生中强地震的危险性．

然而， 在计算过程中存在很多不确定性． 首先， 参数的选取可能带来一定误差， 比如A取0.01， σ_n取10 MPa (Dieterich,Kilgore, 1996)等， 与实际不完全一致； 其次， 库仑应力状态的计算过程依赖模型的选取， 不同研究者可能给出不同的应力变化值(Parsons et al, 2008; Toda et al, 2008; 单斌等, 2009; 万永革等, 2009; 解朝娣等, 2010). 本文仅选用Parsons等(2008)结果进行计算也会带来一定的误差； 最后， 地震目录特别是小震目录的缺失使得对背景场地震活动性的估计不够准确． 所以， 应用条件概率模型会更具可靠性． 但是由于缺乏发震断层的滑移历史、 发震的准周期及上次地震至今的时间， 条件概率模型在目前的研究中， 尤其对我国西南地区来讲， 仍然难以应用．