余震触发的动力学机制和其理论模型.pdfVIP

·272· 中国地球物理2011 (8)地球内部结构及其动力学 余震触发的动力学机制及其理论模型 mechanismofaftershock anditstheoreticalmodel Dynamic triggering 胡才博蔡永恩 HuCaiboCai Yongkang 北京大学地球与空间科学学院 100871 大地震后余震的动力学机制及其时空分布是人们十分关心的问题。基于均匀、弹性半无限空间的地震 位错模型，对于走滑型大地震而言，计算得到的平行于主震断层面上的库仑应力变化(ACFS)在主震断 层两侧都小于零，无法解释主震断层两侧的余震分布。这主要是因为该模型不能考虑介质的不均匀性以及 区域应力场的影响，该模型本质上是一个运动学模型。 实验观测到的岩石样品破裂前的损伤特征，本文从动力学的角度出发，提出了一种研究余震触发动力学机 制的理论模型，认为余震是在主震后区域应力场作用下，由于研究区域内介质的损伤导致的。理论模型将 研究区域分为两个部分：主震震源区及其外部区域。主震震源区由主震断层带、裂隙以及原岩组成，介质 不均匀性较大。本文利用有限元方法研究区域应力场和地震触发因子Il捌的演化，并且假设地震触发因子大 于或者等于l的有限元网格将发生余震。研究区域的边界允许施加位移、应力或者两者混合的弹簧支撑。 因此，该模型是一个动力学模型，可以用来研究在区域应力场作用下，在充分考虑研究区域介质的不均匀 性，特别是主震震源区的不均匀性的情况下，大地震后余震的动力学机制及其时空分布特征。 为了验证本文提出的理论模型，建立了一个的二维平面应变弹性有限元模型，尺寸为1000km×1000km。 位于主震震源区的中央。每个裂隙的长度不等，长10～20km，厚200m，在主震震源区内不规则分布。主震震 的杨氏模量和泊松比分别为95GPa和0．25。有限元模型西边界和南边界法向位移约束为零，东边界和北边界 分别施加1 与此同时，裂隙的剪切模量从29．2GPa降为1GPa，其他区域的材料与主震之前相同。 模拟结果表明： (1)震前地震触发因子在断层带及其端点附近分布不均匀；震后地震触发因子大于或者等于l的区域， 不仅位于主震断层带的两个端点附近的四个叶形区域，而且也存在主震断层带两侧的附近区域。这可以用 来解释走滑型大地震后断层两侧存在大量的余震分布。这是本文提出的余震触发动力学理论模型与常见的 基于地震位错理论的库仑应力变化方法的主要不同之处。 (2)发现在主震震源区内部，最大主应力和正东方向的夹角可以在一30。到30。左右变化，这要比地震 位错得到的主盈利方向变化差不多大一个量级。这可以用来解释一些余震的震源机制为何与主震不同。 (3)在震源区外部，该夹角随离开震源区的距离迅速衰减为零，逐渐与远场(或边界)的主应力方向 一致。这就是远离震源区，余震的震源机制基本不变的原岗。这个结果与实际观测得到的余震震源机制的 空间分布较为一致。 金(20090460130)资助。 参考文献 a1．Anewfiniteelementmodelin andcontinuousevolutionofstress [1】HuC，et studying triggering earthquake Ser．D-Earth 1430-009-0082-3． Sci，2009，doi：10．1007／sI 【2】Hu,C．et of ina crust anewfiniteelement a1．Studyearthquaketriggeringheterogeneoususing Research Letters，2009，80：5．doi：10．1785／gssrl．80．5．799．