2008年汶川大地震单侧破裂过程的动力学机制研究

朱守彪，袁杰

1 中国地震局地壳应力研究所，北京 1000852 中国科学院计算地球动力学重点实验室，北京 1000493 中国地震局地球物理研究所，北京 100081



2008年汶川大地震单侧破裂过程的动力学机制研究

朱守彪1,2，袁杰1,3

1 中国地震局地壳应力研究所，北京 1000852 中国科学院计算地球动力学重点实验室，北京 1000493 中国地震局地球物理研究所，北京 100081

2008年汶川大地震的破裂过程极其发杂，向东北方向的破裂距离长达300 km，而向西南方向的破裂长度很小，呈现出单侧破裂的主要特征.尽管汶川地震破裂呈单侧传播的现象引起许多地震学家的关注，但其物理机制至今还不是十分清楚.本文利用有限单元计算方法，模拟了汶川地震的破裂过程.模型中根据龙门山断裂带两侧(东南侧为四川盆地，西北侧为川西高原)实际的地震波速度来确定模型的介质物性参数，利用目前观测的应力环境来选定初始应力条件.模拟结果表明：破裂在汶川地震的震中处成核后，先向断层两侧自发传播，但向东北方向的传播距离明显大于向西南方向；断层面上的正应力在东北方向(破裂的正方向)随着传播距离的增大而不断减小，位错速率随着破裂的传播距离而越来越大，其脉冲变得越来越尖锐，即产生了Weertman脉冲.研究结果显示：由于这种脉冲的出现，破裂在正方向上(东北方向)能够自己放大、自己愈合、自行维持，摩擦热极小，所以破裂能够沿着东北方向一直传播，直到应力场方位发生变化，不利于破裂时才最后终止.但在西南方向，破裂过程中断层面上的正应力增大，阻碍破裂继续扩展.最后就出现了汶川地震中破裂朝东北方向单侧优势传播的基本格局.模拟结果还表明：若断层面两侧介质均匀，则破裂向两侧是对称传播，且破裂距离很短，因此这种情况无法产生像汶川大地震那样的特大地震.因此，文中的模拟结果表明龙门山断裂带两侧的物性差异是造成汶川大地震单侧传播的决定性因素.断层两侧物性差异(bimaterial contrast)影响断层破裂过程的研究对于深入认识地震动力学过程、地震灾害预测及评估等有重要的科学意义.

破裂动力学过程；断层两侧物性差异；Weetman脉冲；单侧破裂；汶川地震；有限单元

1 引言

2008年发生的汶川大地震，断层破裂过程非常发杂，地表所能观测到的破裂长度达300余千米.地震发生后，国内外很多学者陆续利用多种方法和不同资料(如：地震波形、地表破裂、强地面运动、GPS和InSAR资料等)对汶川地震的破裂过程的运动学模型进行了研究(Ji and Hayes,2008；王卫民等，2008；张勇等，2008；杜海林等，2009；Shen et al.，2009；Xu et al.,2009；赵翠萍等，2009；Nakamura et al.，2010；Zhang and Ge，2010；Wang et al.，2011；Fielding et al.,2013;Pang et al.,2014;Zhang et al.,2014;Zhang and Wang,2015)，取得了多种反演结果.尽管不同作者给出的结果不尽相同，但描述的汶川地震破裂的运动学特征大体如下：

破裂首先在位于映秀—北川断裂带上的震源处起始，在最先的开始阶段，表现为双侧的破裂形式，在东北及西南的两个方向同时产生破裂，但向西南方向的破裂很快就终止.而在断层的东北方向，破裂继续前行，跨越了多个断层阶区，破裂经过的断层面形状复杂(如弯曲、拐折、交叉等)，破裂最后在青川断裂带上终止.汶川大地震破裂的总长度大约为300 km,破裂持续时间大约100 s;破裂速度在不同的地方差异很大(但基本上都小于剪切波速度)，静态位错分布也很不均匀，断层的倾角由南向北越来越大，位错方向由近纯逆冲逐渐过渡到走滑类型.

尽管运动学模型告诉我们断层破裂的长度、破裂速度、位错大小等很多物理量，也为认识汶川地震及抗震救灾发挥了重要作用.但是，是什么因素造成了这样的破裂格局？或者说，汶川地震破裂过程中，什么物理量起着决定性作用？这些问题至今依然没有给出回答.实际上，这些问题在物理上是动力学问题，只有通过对汶川地震的破裂过程进行动力学分析，才有可能加以解决.汶川地震发生后，其震源破裂的动力学过程倍受人们的关注，国内外不少专家对此进行了探索，获得了一些认识.

Duan(2010)利用三维有限元方法对汶川地震的破裂动力学过程进行了研究.模型中，汶川地震的破裂面假定为一条平直断层面，其倾角为33°；介质材料均匀，初始应力场不均匀并随深度变化.通过研究发现，随着初始应力场方向的改变，断层的破裂速度会发生很大的变化，初始应力场的方向还可以导致截然不同的位错格局及地表运动形态.特别是当初始场的指向与断层走向之间的夹角处在最不利于断层破裂的方位时，破裂会自然地停止.汶川地震在断层东北端的终止以及两处破坏最强的区域，也许就是初始场的方向在起着作用.Wen等(2012) 利用显示的三维有限元方法，所用模型与Duan(2010)的模型基本相同，设定了不均匀的初始应力条件，模拟结果表明初始应力场决定了破裂速度、位错的空间分布等主要破裂特征.

Zhu和Zhang(2010，2013)利用改进的有限单元模型，分两步(即先反演，后正演)分别模拟了汶川地震震间变形及同震破裂过程.第一步利用大地测量资料(水准测量及GPS观测)作约束，反演震间的最佳数值模型，获得汶川地震前的应力、应变状态；然后据此最佳模型，进行正演计算模拟断层的自发破裂过程.计算结果显示：下部缓倾角断层上的破裂对上部陡倾角上的错动有触发作用，破裂过程中，断层面上的剪切应力和法向应力都有利于断层的错动，这可能也是高倾角逆冲断层上发生大地震的原因之一.

虽然上述学者对汶川地震破裂动力学机制研究方面作了不少探索，取得一定的进展，但汶川地震破裂的最基本特征为单侧破裂(或者认为是不对称的双侧破裂)，对这个问题至今未见有文献报道.此外，我们注意到龙门山断裂带两侧介质物理性质差异很大(Lei and Zhang,2009;Zhu and Zhang,2013;Liu et al.,2014)，而断层两侧介质物理性质差异(bimaterial contrast)会严重影响断层的破裂过程和破裂特征(Weertman,1980;Adam,1995;Andrews and Ben-Zion,1997;Ben-Zion,2001;Ranjith and Rice,2001;Ben-Zion and Shi,2005;Ben-Zion et al.,2007;Rubin and Rice,2007;Ampuero and Ben-Zion,2008;Duan,2008;Ma and Beroza,2008;Olsen-Kettle et al.,2008;Brietzke et al.,2009;DeDontney et al.,2011;Langer et al.,2012;Scholz,2014;Yuan and Zhu,2016).所以本文将利用有限元方法，考虑龙门山断裂带两侧介质物性的差异，模拟汶川地震的破裂过程，试图定量给出汶川地震为什么会出现单侧破裂的动力学机制.

2 有限元模型

为研究汶川地震单侧破裂的物理机制，首先必须对实际复杂问题进行抽象与简化.

汶川地震中，断层破裂规模巨大，断层几何形状复杂多变，既有平面也有曲面，断层面倾角自南向北不断发生变化，断层在空间上还存在不连续分布等多种特征.此外，地震中断层的位错主要表现为由逆冲向走滑性质过渡：从震源处的几乎纯逆冲，然后不断地由南向北过渡到几乎纯走滑类型.研究中，为抓住主要矛盾，重点突出汶川地震单侧破裂过程的物理本质，研究中将实际的三维地质体简化为二维模型(通过震源的水平面)，将复杂的断层几何简化为平面，复杂的断层破裂过程简化为简单的走滑类型.模拟中，同时还假定汶川大地震之前龙门山断裂带及其附近区域处在均匀的应力环境里，这样文中研究的破裂问题就简化为标准的滑开型(Ⅱ-型)破裂模式.计算时，首先让破裂在汶川地震的震源处成核，然后让破裂自由地向东南—西北两个方向自发地传播、演化.

2.1 基本控制方程

对于二维模型，当模型四周利用无限单元来吸收地震波时，破裂动力学过程的能量表达形式为(Hibbitt et al.,2006)

(1)

此外假定断层面上的法向位移为连续，即

(2)

断层面上的质点所处的运动状态表达式为

(3)

其中σn法向接触应力，τ为剪切应力,μ为摩擦系数.

利用有限元方法对于上述 (1)—(3)的方程进行求解,即可得到模型中位移、速度、加速度，应力，应变以及能量等物理量.

2.2 有限元时间积分方法及计算时间步长

2.3 有限元模型

图1 有限元模型图中坐标原点处的粗实线为汶川地震的震中(成核区).龙门山断裂带的走向沿X轴正方向(东北方向).Fig.1 Sketch of the finite element modelThe thicksolid line at the origin of coordinates stands for the epicenter of the 2008 Wenchuan earthquake.The strike of the Longmen Shan fault is in the direction of X-axis,namely northeast.

根据汶川地震的震中位置以及龙门山断裂带的实际构造，研究中构造了二维有限元模型，如图1所示.模型空间几何尺度为500 km×500 km.为减少地震波辐射对破裂行为的影响，模型四周施加了吸收边界(利用ABAQUS中的无限单元来实施).地震破裂的成核位置位于模型的中心(图中粗黑线段).有限元网格全部为三角形单元，单元最小边长为100 m,单元最大边长为500 m.破裂带长度为500 km,利用ABAQUS中接触单元来表示.由于初始应力场影响破裂行为，所以模拟中根据对龙门山断裂带及附近地区地应力的观测结果(Wu et al.,2009)来选择初始应力场.

2.4 摩擦本构关系及正则化

摩擦本构关系决定着破裂的形式及其演化特征，根据前人对断层两侧物性不同发生破裂情况的模拟经验(Rubin and Ampuero,2007;Ampuero and Ben-Zion,2008;Ma and Beroza,2008;Finzi and Langer,2012;Gabriel et al.,2012；Langer et al.,2012)，研究中选择速度-状态相依的摩擦本构关系，具体数学表达式为

(4)

其中μs为静摩擦系数，V为滑移速率，Vc为特征速度，αf及βf表示演化效果的常数，θ为状态参量，tc表示特征时间尺度.

通常情况下，当断层两侧介质物性不同时，断层自发破裂会出现所谓的Adams不稳定现象(Adams,1995，1998;Harris and Day,1997;Simõnes and Martins,1998；Cochard and Rice,2000;Ben-Zion,2001;Ranjith and Rice,2001)，这种现象在物理上表现为破裂产生的脉冲随着传播距离变得更加尖锐并且发散.在数学上，该系统对微扰的响应是病态的(ill-posedness)，并且数值计算结果随着网格的变化而变化(Ben-Zion and Huang,2002).为此，要得到稳定的解，需要进行正则化处理.根据前人数值模拟的结果(Cochard and Rice,2000;Ranjith and Rice,2001;Rubin and Ampuero,2007;Ampuero and Ben-Zion,2008；Brietzke et al.,2009)，计算时通过减缓断层面上的剪切应力随正应力变化的过程来实现正则化，从而提高计算的稳定性.实际数值模拟中，断层面上的剪切应力随时间的变化率为(DeDontney et al.,2011)

(5)

其中σ、τ、μf分别为正应力、剪切应力以及摩擦系数；t*为正则化时间尺度，取值为：4Δx/Cs1(其中，Δx为滑移面上的单元尺度，Cs1为断层两侧中软弱介质的S波速度)，理论上t*应远大于计算的时间步长，通常在实际计算时大一个数量级即可.具体数值见表1.

2.5 地震成核过程

在数值模拟中，要使破裂过程启动，通常利用减小摩擦系数来实现.具体计算中是将成核区内的静态摩擦系数突然减小，表达式为

(6)

破裂在成核区一旦形成，在其外的区域里，破裂行为由摩擦本构关系及应力环境来决定.破裂可以自发地向断层两侧传播、演化，直至最后终止.

3 有限元模拟结果

为模拟2008年汶川地震的破裂过程，我们根据前人对龙门山地区的地震层析成像、以及其他地球物理反演结果等资料，建立具体的有限元模型.模型中龙门山断裂带及附近地区介质的地震波速度根据地震层析成像结果(以汶川地震震源深度处)来选取(Lei and Zhao,2009;Liu et al.,2014)，数值见表1.这样龙门山断裂带两侧(川西高原与四川盆地)介质的地震波速度相差大约为20%.模拟中关于摩擦行为的基本参数参考前人的模拟结果来选定(Rubin and Ampuero,2007;Ampuero and Ben-Zion,2008;Ma and Beroza,2008;Finzi and Langer,2012).

图2为模拟给出的汶川地震震中(成核中心)两侧断层面上的滑移量分布(为清楚起见，该图只展示了局部区域，以下同).由图可知，破裂为明显的不对称分布，向北东方向的破裂距离明显大于向西南方向的距离，位错量也是在东北方向的大于向西南向的.因此许多学者认为汶川地震为单侧破裂类型(实际为不对称的双侧破裂).破裂向东北方向传播距离远，而向西南方向传播的距离很小，仅10余千米.

图2 震中两侧断层面上的位错在不同时刻的空间分布Fig.2 Spatial distribution of fault displacements on both sides of the epicenter at different times

参量物理性质参数(单位)数值四川盆地S波速度Vs1(m·s-1)4120四川盆地P波速度Vp1(m·s-1)7150川西高原S波速度Vs2(m·s-1)3380川西高原P波速度Vp2(m·s-1)5850初始剪应力τ0(MPa)69初始正应力σ0(MPa)200介质密度ρ(kg·m-3)2900静摩擦系数μs0.6直接演化系数αf0.01演化系数βf0.40特征速度Vc(m·s-1)0.7特征时间tc(s)0.3速率弱化特征长度Dc(m)0.21成核区大小Lnucl(km)1正则化时间尺度t*(s)1.2×10-3

此外，图2还显示，断层面上任一点的位错持续时间很短，远远小于整个断层破裂的持续时间，因此该破裂为脉冲型(pulse-like)破裂(袁杰和朱守彪，2014a).另外，我们还可以从图2中大体知道破裂的传播速度为2.7 km·s-1,因此属于亚剪切破裂(破裂速度小于剪切波速度).该破裂速度与张勇等(2008)以及Xu等(2009)通过反演给出的运动学震源破裂模型的破裂速度大体相当.

根据图1所示的应力边界条件以及川西高原所在一侧是物质相对较柔软的格局，我们知道本研究中破裂的正方向应该是朝龙门山断裂带的东北方向，即图中的x轴正向；那么破裂的负方向就是x轴的反向(Andrews and Ben-Zion,1997;Ben-Zion and Andrews,1998;Ben-Zion,2001;Ranjith and Rice,2001;Ben-Zion and Shi,2005;Shi and Ben-Zion,2006;Ben-Zion et al.,2007;Rubin and Ampuero,2007;Ampuero and Ben-Zion,2008).对于双层不同性质的断层，破裂的正向应为破裂的优势取向，即破裂在东北方向应有较好的破裂效果，图2的结果正是说明了这一点.

图3为滑移速率在断层上的不同位置的分布.由图可见，在破裂的正方向上(东北向)，速度脉冲的幅值随着破裂传播距离的增大而越来越大，说明破裂具有自己放大的能力，这样地震破裂就能持续很长距离而不会衰减；此外，位错速率脉冲的宽度(最大振幅一半处的脉冲宽度)在正方向上是越来越小、越来越尖锐，说明破裂能量向低波长方向转移.但在破裂的负方向，脉冲变小，最后终止.图3还显示，每个脉冲起始后，在很短的时间里就达到最大值，然后迅速降低为0(意味着断层很快就愈合).可见，这种脉冲型破裂具有自己愈合的能力，这正是脉冲型破裂(pulse-like)的典型特征.图3所示的破裂脉冲就是通常所称谓的Weertman脉冲，它能自己愈合、自行维持、自行放大，对地震破裂过程的影响非常深远(Weertman,1980;Andrews and Ben-Zion,1997;Ben-Zion,2001;Shi and Ben-Zion,2006;Ma and Beroza,2008;Ampuero and Ben-Zion,2008;Brietzke et al.,2009;DeDontney et al.,2011;Finzi and Langer,2012;Langer et al.,2012;Scholz,2014;Yuan and Zhu,2016),一直受到国际地震学界的广泛重视.由于破裂过程中产生的Weertman 脉冲具有自己放大、自行维持的属性，所以这种破裂不易于停止，这就是为什么汶川地震的震级是如此巨大的重要因素.

图3 断层面上滑移速率沿着断层上不同位置时的分布(时间间隔为1.5 s) Fig.3 Distribution of slip rates at different locations along the fault at regular time interval (1.5 s)

图4 断层面上的正应力及剪切应力在不同时刻的分布(a) 正应力；(b) 剪切应力.Fig.4 Distribution of normal and shear stresses on the surface of the fault at different times(a) Normal stress;(b) Shear stress.

图5 当模型中介质性质均匀时，破裂传播过程中断层面上的正应力及剪切应力在不同时刻的分布(a) 正应力；(b) 剪切应力.Fig.5 Distribution of normal and shear stresses on the surface of the fault at different times in rupture propagation when the material in the model is not uniform(a) Normal stress;(b) Shear stress.

图4给出了断层面上的正应力及剪切应力在不同时刻的分布.图中清楚地显示，断层上的正应力是随着时间变化的.但在模拟时模型施加的初始应力场为均匀场，即断层面上的正应力应该是处处相等，不会发生变化.图4结果表明，随着破裂传播距离的增大，张性正应力变化会越大越大，特别是当断层很长，破裂传播的距离很远时，可以让断层处在拉张的应力状态，这样破裂就会变得更加容易，摩擦产生的热几乎为0，断层的双侧介质效应会愈加明显.所以震级俞大，Weertman脉冲的影响就会愈加显著(Ranjith and Rice,2001;Ben-Zion and Huang，2002).

为了搞清楚断层面上的正应力变化的机理，我们又构建了另外一个模型.该模型中断层两侧介质的物理性质完全一样，即整个模型介质均匀(其他模型参量均保持不变).这时，模拟给出的破裂传播过程中断层面上的应力分布如图5所示.通过与图4的比较可以发现，当模型中介质均匀时，破裂传播过程中断层面上的正应力不随时间及位置的变化而变化.

实际上，早在1980年，对于断层两侧介质性质不同的破裂传播问题，Weertman(1980)就通过数学解析推导，给出了断层面上的剪切应力变化及正应力变化的数学表达式为

(7)

(8)

Andrews和Ben-Zion(1997)将介质在断层物性较软弱一侧的运动形象地比作为地毯在地板上的移动.当地毯在地面上移动时则会出现褶皱(wrinkle)，介质瞬时从接触面离开，因此在褶皱传播方向上就产生了一个净的剪切位移，导致正应力变小(Andrews and Ben-Zion,1997;Ben-Zion and Andrews,1998;Cochard and Rice,2000;Ranjith and Rice,2001;Ben-Zion and Huang,2002;Ben-Zion and Shi,2005).

事实上，正是由于龙门山断裂带两侧介质物性的差异，导致了破裂在东北方向的传播过程中断层面上的正应力减小，特别是当传播距离越来越远时，甚至会出现完全的拉张正应力状态(Rubin and Ampuero,2007)，这非常有利于破裂的传播(破裂的正方向)，使得破裂能够自己维持、自己愈合、自行放大，破裂过程中摩擦产生的热能极小，导致破裂传播的距离很远，最终导致地震的震级很大(Brietzke et al.,2009).

4 讨论与结论

本文运用完全动力学方法模拟了断层的破裂过程，显示了龙门山断裂带两侧介质物理性质差异是导致汶川地震单侧破裂的决定性因素.为研究介质物性差异对地震破裂过程的影响，下面再展示一下上文提到的将模型中断层两侧材料性质视为均匀的模型中位错在空间的分布结果.图6是当模型中介质均匀时，位错在空间的分布图.图中显示，当断层两侧介质性质完全相同时(其他模型参量不变)，破裂确实是向两侧传播，并且在空间上是对称分布.但传播的距离非常之短，总共只有16 km,因此在这种情况下，产生不了像汶川地震(Ms8.0)那样的特大地震.由此可见，断层两侧物性的差异对于破裂过程的影响是多么巨大.通过图6与图2的比较，我们清楚地认识到龙门山地震带两侧介质物理性质的差异导致了汶川大地震的单侧传播过程(或不对称的双侧破裂).

在此之上，文中还专门针对断层两侧介质物性差异程度对破裂传播过程的影响进行了模拟分析.当断层两侧介质物性差异分别为5%、10%以及15%时，断层破裂都会出现不对称的双侧破裂行为，但不对称性随着物性差异程度的增大而不断增大(Yuan and Zhu,2016).

龙门山断裂带及周边地区的应力状态应该是不均匀的.但Duan(2010)以及Wen等(2012)的研究结果表明，初始应力场的不均匀性会严重影响破裂的静态位错、应力降以及破裂速度的分布，所以文中所用的均匀初始场不会影响对汶川地震单侧破裂传播过程这一物理本质的认识.此外，文中的简单模型导致断层破裂很难终止，之所以汶川地震最后终止于青川断裂带，Duan(2010)认为是初始应力场的方位改变造成了破裂的终止.所以初始应力场的非均匀性，不影响我们对汶川地震破裂过程不对称的数值模拟.但是，需要特别指出的是，汶川地震的向南终止绝对不是初始应力场的方位发生变化造成的，因为在汶川地震震中附近一直向其西南方向很大的范围里，应力场方位几乎是相同的(Wu et al.,2009).

此外，本文在模拟过程中，将龙门山断裂带的逆冲兼走滑性质的断层简化为纯走滑断层，将实际的三维断层结构简化为二维.这种简化的合理性主要是：首先从地图上看，汶川地震的破裂是由南向北传播，所以从整体上看该地震可以等价为一个走滑型事件.其次，在地震的起始阶段，破裂表现为逆冲型，但随着时间的推移，走滑特征越来越明显，特别是到破裂的后期，几乎就是纯走滑性质.所以，文中抓住汶川地震单侧破裂这个根本性的科学问题，将其简化处理，这不会从本质上影响我们对这个问题的认识.当然，今后可以利用三维模型，对破裂过程进行更为符合实际的模拟，对文中的结果进行进一步的验证.顺便指出，震源破裂过程不对称性的运动学研究成果非常丰富(高原等，1997；Gao et al.,2000)，这些结果对于深入研究震源破裂的动力学过程提供非常重要的参考.

图6 当断层两侧介质均匀时，位错随在不同时刻的空间分布图Fig.6 Spatial distribution of displacements on both sides of the fault from at different times when the material on both sides of the fault is uniform

根据文中的计算与分析，得出如下初步结论：

利用有限单元计算方法，可以模拟断层的自发破裂动力学过程.尽管对于断层两侧介质性质不同的破裂问题，模拟中会出现Adams不稳定现象，但通过正则化可以得到稳定的数值解.

龙门山断裂带的西北侧(青藏高原)介质柔软，而东南侧(四川盆地)介质刚硬.由于断层两侧介质物性的不同，导致了汶川地震破裂的不对称性空间分布.

模拟结果表明，汶川地震在破裂过程中，断层面上的正应力在东北方向(破裂的正方向)随着传播距离的增大而不断减小，位错速率随着破裂的传播距离而越来越大，其脉冲变得越来越尖锐，即产生了能够自己放大、自己愈合、自行维持的Weertman脉冲，所以汶川地震在东北方向破裂的距离很大.但在西南方向，破裂过程中断层面上的正应力在不断增大，不利于破裂继续扩展.所以，在龙门山断裂带及附近区域的应力环境不变时，断裂带两侧的物性差异造成了汶川大地震中发生了单侧传播现象.研究断层两侧物性差异影响断层的破裂行为，对于深入认识地震动力学过程、分析未来地震的震级大小、预测地震灾害等有重要的科学意义.

致谢 两位审稿专家提出了十分宝贵的建议，在此表示衷心的感谢！

Adams G G.1995.Self-excited oscillations of two elastic half-spaces sliding with a constant coefficient of friction.Journal of Applied Mechanics,62(4):867-872.

Adams G G.1998.Steady sliding of two elastic half-spaces with friction reduction due to interface stick-slip.Journal of applied Mechanics,65(2):470-475.

Ampuero J P,Ben-Zion Y.2008.Cracks,pulses and macroscopic asymmetry of dynamic rupture on a bimaterial interface with velocity-weakening friction.Geophys.J.Int.,173(2):674-692.

Andrews D J,Ben-Zion Y.1997.Wrinkle-like slip pulse on a fault between different materials.Journal of Geophysical Research,102(B1):553-571.Ben-Zion Y.2001.Dynamic ruptures in recent models of earthquake faults.Journal of the Mechanics and Physics of Solids,49(9):2209-2244.

Ben-Zion Y,Huang Y Q.2002.Dynamic rupture on an interface between a compliant fault zone layer and a stiffer surrounding solid.Journal of Geophysical Research:Solid Earth (1978—2012),107(B2):ESE 6-1—ESE 6-13.

Ben-Zion Y,Peng Z,Zhao P,et al.2007.Variations of the velocity contrast and rupture properties of M6 earthquakes along the parkfield section of the san andreas fault.∥American Geophysical Union,Fall Meeting 2007.AGU.

Ben-Zion Y,Shi Z Q.2005.Dynamic rupture on a material interface with spontaneous generation of plastic strain in the bulk.Earth and Planetary Science Letters,236(1-2):486-496.

Brietzke G B,Cochard A,Igel H.2007.Dynamic rupture along bimaterial interfaces in 3D.Geophysical Research Letters,34(11):L11305.

Brietzke G B,Cochard A,Igel H.2009.Importance of bimaterial interfaces for earthquake dynamics and strong ground motion.Geophysical Journal International,178(2):921-938.

Cochard A,Rice J R.2000.Fault rupture between dissimilar materials:Ill-posedness,regularization,and slip-pulse response.J.Geophys.Res.,105(B11):25891-25907.

DeDontney N,Rice J R,Dmowska R.2012.Finite element modeling of branched ruptures including off-fault plasticity.Bulletin of the Seismological Society of America,102(2):541-562.

DeDontney N,Templeton-Barrett E L,Rice J R,et al.2011.Influence of plastic deformation on bimaterial fault rupture directivity.Journal of Geophysical Research:Solid Earth,116(B10):B10312.

Dor O,Rockwell T K,Ben-Zion Y.2006.Geological observations of damage asymmetry in the structure of the San Jacinto,San Andreas and Punchbowl faults in Southern California:a possible indicator for preferred rupture propagation direction.Pure and Applied Geophysics,163(2-3):301-349.

Dor O,Yildirim C,Rockwell T K,et al.2008.Geological and geomorphologic asymmetry across the rupture zones of the 1943 and 1944 earthquakes on the North Anatolian Fault:possible signals for preferred earthquake propagation direction.Geophysical Journal International,173(2):483-504.

Du H L,Xu L S,Chen Y T.2009.Rupture process of the 2008 great Wenchuan earthquake from the analysis of the Alaska-array data.Chinese J.Geophys.(in Chinese),52(2):372-378.

Duan B C.2008.Asymmetric off-fault damage generated by bilateral ruptures along a bimaterial interface.Geophysical Research Letters,35(14):L14306.

Duan B C.2010.Role of initial stress rotations in rupture dynamics and ground motion:A case study with implications for the Wenchuan earthquake.J.Geophys.Res.,115(B5):B05301,doi:10.1029/2009JB006750.

Fielding E J,Sladen A,Li Z H,et al.2013.Kinematic fault slip evolution source models of the 2008 M7.9 Wenchuan earthquake in China from SAR interferometry,GPS and teleseismic analysis and implications for Longmen Shan tectonics.Geophysical Journal International,194(2):1138-1166.

Finzi Y,Langer S.2012.Predicting rupture arrests,rupture jumps and cascading earthquakes.J.Geophys.Res.,117:B12303,doi:10.1029/2012JB009544.

Gabriel A A,Ampuero J P,Dalguer L A,et al.2012.The transition of dynamic rupture styles in elastic media under velocity-weakening friction.J.Geophys.Res.,117:B09311,doi:10.1029/2012JB009468.

Gao Y,Wu Z,Liu Z,et al.2000.Seismic source characteristics of nine strong earthquakes from 1988 to 1990 and earthquake activity since 1970 in the Sichuan-Qinghai-Xizang (Tibet) zone of China.Pure and Applied Geophysics,157(9):1423-1443.

Gao Y,Zhou H L,Liu Z.1997.Analysis on rupture feature of unilateral rupture source.Journal of the Graduate School of the Chinese Academy of Science (in Chinese),14(1):39-42.

Hibbitt H D,Karlsson B I,Sorensen P.2006.ABAQUS theory manual,version 6.3.Pawtucket,Rhode Island,USA.

Ji C,Hayes G.2008.Preliminary result of the May 12,2008 Mw7.9 eastern Sichuan,China earthquake.US Geol.Surv.

Langer S,Olsen-Kettle L,Weatherley D.2012.Identification of supershear transition mechanisms due to material contrast at bimaterial faults.Geophysical Journal International,190(2):1169-1180.

Lei J S,Zhao D P.2009.Structural heterogeneity of the Longmenshan fault zone and the mechanism of the 2008 Wenchuan earthquake (MS8.0).Geochem.Geophys.Geosyst.,10(10):Q10010,doi:10.1029/2009GC002590.

Liu Q Y,van der Hilst R D,Li Y,et al.2014.Eastward expansion of the Tibetan Plateau by crustal flow and strain partitioning across faults.Nature Geoscience,7(5):361-365.

Ma S,Beroza G C.2008.Rupture dynamics on a bimaterial interface for dipping faults.Bulletin of the Seismological Society of America,98(4):1642-1658.

Nakamura T,Tsuboi S,Kaneda Y,et al.2010.Rupture process of the 2008 Wenchuan,China earthquake inferred from teleseismic waveform inversion and forward modeling of broadband seismic waves.Tectonophysics,491(1-4):72-84,doi:10.1016/j.tecto.2009.09.020.

Olsen-Kettle L M,Weatherley D,Saez E,et al.2008.Analysis of slip-weakening frictional laws with static restrengthening and their implications on the scaling,asymmetry,and mode of dynamic rupture on homogeneous and bimaterial interfaces.Journal of Geophysical Research:Solid Earth (1978—2012),113(B8).

Ranjith K,Rice J R.2001.Slip dynamics at an interface between dissimilar materials.Journal of the Mechanics and Physics of Solids,49(2):341-361.

Rubin A M,Ampuero J P.2007.Aftershock asymmetry on a bimaterial interface.J.Geophys.Res.,112(B5).

Scholz C H.2014.The rupture mode of the shallow large-slip surge of the Tohoku-Oki earthquake.Bulletin of the Seismological Society of America,104(5):2627-2631.

Shen Z K,Sun J B,Zhang P Z,et al.2009.Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake.Nature Geoscience,2(10):718-724,doi:10.1038/ngeo636.

Shi Z Q,Ben-Zion Y.2006.Dynamic rupture on a bimaterial interface governed by slip-weakening friction.Geophysical Journal International,165(2):469-484.

Simõnes F M F,Martins J A C.1998.Instability and ill-posedness in some friction problems.Int.J.Eng.Sci.,36(11):1265-1293.

Wang W M,Zhao L F,Li J,et al.2008.Rupture process of the MS8.0 wenchuan earthquake of Sichuan,China.Chinese J.Geophys.(in Chinese),51(5):1403-1410.

Weertman J.1980.Unstable slippage across a fault that separates elastic media of different elastic constants.Journal of Geophysical Research:Solid Earth (1978—2012),85(B3):1455-1461.

Wen Y Y,Oglesby D,Duan B,et al.2012.Dynamic rupture simulation of the 2008 MW7.9 Wenchuan earthquake with heterogeneous initial stress.Bulletin of the Seismological Society of America,102(4):1892-1898.

Wu M L,Zhang Y Q,Liao C T,et al.2009.Preliminary results of in-situ stress measurements along the Longmenshan fault zone after the Wenchuan MS8.0 earthquake.Acta Geologica Sinica (English Edition),83(4):746-753.

Xu Y,Koper K D,Sufri O,et al.2009.Rupture imaging of the MW7.9 12 May 2008 Wenchuan earthquake from back projection of teleseismic P waves.Geochem.Geophys.Geosyst.,10:Q04006,doi:10.1029/2008GC002335.

Yuan J,Zhu S B.2014a.FEM simulation of the dynamic processes of fault spontaneous rupture.Chinese J.Geophys.(in Chinese),57(1):138-156,doi:10.6038/cjg20140113.

Yuan J,Zhu S B.2014b.Effects of stepover on rupture propagation.Chinese J.Geophys.(in Chinese),57(5):1510-1521,doi:10.6038/cjg20140515.

Yuan J,Zhu S B.2016.Distributions of strong ground motion due to dynamic ruptures across a bimaterial fault:Implications for seismic hazard analyses.J.Asian Ear.Sci.,in review.

Zaliapin I,Yehuda Ben-Zion Y.2011.Asymmetric distribution of aftershocks on large faults in California.Geophys.J.Int.,185(3):1288-1304.

Zhang H,Ge Z X.2010.Tracking the rupture of the 2008 Wenchuan earthquake by using the relative back-projection method.Bull.Seismol.Soc.Am.,100(5B):2551-2560,doi:10.1785/0120090243.

Zhang Y,Feng W P,Xu L S,et al.2009.Spatio-temporal rupture process of the 2008 great Wenchuan earthquake.Science in China Series D:Earth Sciences,52(2):145-154.

Zhang Y,Wang R J.2015.Geodetic inversion for source mechanism variations with application to the 2008 MW7.9 Wenchuan earthquake.Geophysical Journal International,200(3):1627-1635.

Zhang Y,Wang R J,Zschau J,et al.2014.Automatic imaging of earthquake rupture processes by iterative deconvolution and stacking of high-rate GPS and strong motion seismograms.Journal of Geophysical Research:Solid Earth,119(7):5633-5650.

Zhao C P,Chen Z L,Zhou L Q,et al.2010.Rupture process of the Wenchuan 8.0 earthquake of Sichuan China:The segmentation feature.Chinese Science Bulletin,55(3):284-292,doi:10.1007/s11434-009-0425-7.

Zhu S B.2013.Numerical simulation of dynamic mechanisms of the 2008 Wenchuan MS8.0 earthquake:implications for earthquake prediction.Natural Hazards,69(2):1261-1279.

Zhu S B,Zhang P Z.2010.Numeric modeling of the strain accumulation and release of the 2008 Wenchuan,Sichuan,China,Earthquake.Bulletin of the Seismological Society of America,100(5B):2825-2839.

Zhu S B,Zhang P Z.2013.FEM simulation of interseismic and coseismic deformation associated with the 2008 Wenchuan Earthquake.Tectonophysics,584:64-80.

附中文参考文献

杜海林,许力生,陈运泰.2009.利用阿拉斯加台阵资料分析2008年汶川大地震的破裂过程.地球物理学报,52(2):372-378.

高原,周蕙兰,刘振.1997.单侧破裂震源的破裂特征分析.中国科学院研究生院学报,14(1):39-42.

王卫民,赵连锋,李娟等.2008.四川汶川8.0级地震震源过程.地球物理学报,51(5):1403-1410.

袁杰,朱守彪.2014a.断层自发破裂动力过程的有限单元法模拟.地球物理学报,57(1):138-156,doi:10.6038/cjg20140113.

袁杰,朱守彪.2014b.断层阶区对震源破裂传播过程的控制作用研究.地球物理学报,57(5):1510-1521,doi:10.6038/cjg20140515.张勇,冯万鹏,许力生等.2008.2008年汶川大地震的时空破裂过程.中国科学:D辑,38(10):1186-1194.

赵翠萍,陈章立,周连庆等.2009.汶川MW8.0级地震震源破裂过程研究:分段特征.科学通报,54(22):3475-3482.

(本文编辑 张正峰)

Mechanisms for the fault rupture of the 2008 Wenchuan earthquake (MS=8.0) with predominately unilateral propagation

ZHU Shou-Biao1,2,YUAN Jie1,3

1 Institute of Crustal Dynamics,China Earthquake Administration,Beijing 100085,China2 Key Lab of Computational Geodynamics,Chinese Academy of Sciences,Beijing 100049,China3 Institute of Geophysics,China Earthquake Administration,Beijing 100081,China

The rupture processes of the 2008 Wenchuan earthquake are very complex.The rupture propagated northeastward a large distance as long as 300 km,while the rupture length is short in the southwest direction.Although a great deal of attention has been paid to why the rupture propagated preponderantly northeastward during this event,the physical mechanism remains unclear.By means of a finite element model in which bimaterial contrast across the Longmen Shan fault is taken into account,this study suggests that the nucleation process is initiated at the center of the fault,ruptures then spread out outward spontaneously in both directions.The simulation results show that the different materials between both sides of the fault lead to tensile changes of normal stresses on the fault,which enhances the propagation of the rupture.We find that bimaterial mechanism is important for the earthquake rupture,and the mode II rupture evolves with propagation distance along a bimaterial interface to a unilateral wrinkle-like Weertman pulse in the direction of slip on the more compliant side of the fault,namely in the positive direction (i.e,the northeast direction).This Weertman pulse can be self-amplified,self-sustained and self-healing,which gives rise to little frictional energy and long propagation distance.This may be the reason why the 2008 Wenchuan earthquake is a unilateral fault rupture and has a high seismic magnitude.In addition,the modelling results suggest that the rupture distance would be much smaller if the material in the model is homogeneous,in which no huge earthquakes can occur like the 2008 Wenchan event.Hence,the study of the bimaterial effect will play an important role in understanding fault dynamics,seismic hazard estimation and prediction.

Dynamic rupture process;Bimaterial contrast;Weertman pulse;Unilateral fault rupture;Wenchuan earthquake;Finite element

朱守彪，袁杰.2016.2008年汶川大地震单侧破裂过程的动力学机制研究.地球物理学报，59(11):4063-4074,

10.6038/cjg20161111.

Zhu S B,Yuan J.2016.Mechanisms for the fault rupture of the 2008 Wenchuan earthquake (MS=8.0) with predominately unilateral propagation.Chinese J.Geophys.(in Chinese),59(11):4063-4074,doi:10.6038/cjg20161111.

国家自然科学基金项目(41574041)、北京市自然科学基金项目(8152034)以及广西科技攻关项目(桂科攻1377002)共同资助.

朱守彪，男，研究员，理学博士，博士生导师，现主要从事地球动力学及地震活动性研究.E-mail:zhusb@pku.edu.cn;zhushoubiao@gmail.com

10.6038/cjg20161111

P315，P542

2016-04-05，2016-07-20收修定稿