2013年岷县漳县6.6级地震孕育

张希 张四新 秦姗兰 贾鹏 李瑞莎

摘要：基于西秦岭构造区区域水准、长水准剖面观测资料，结合GPS速度场约束和跨断层短水准观测数据，借助三维负位错联合反演，分析2013年岷县漳县6.6级地震前相关断裂的中长期闭锁背景，讨论断裂活动的中—短临动态演变以及震后变化。结果表明：1999—2006年，岷县漳县地震前的震区所在块体即西秦岭北缘断裂与临潭—宕昌断裂之间水平面最大剪应变率、挤压性质的面应变率、垂向梯度率都较大，与构造和动力学背景一致的三维差异运动明显，有利于能量积累和孕育强震。包含漳县地区的西秦岭北缘断裂中段、临潭—宕昌断裂中东段即发震构造的主断裂段闭锁明显，挤压逆断兼左旋，显示中长期能量积累背景。跨西秦岭北缘断裂漳县段至临潭—宕昌断裂中东段岷县地区的长水准剖面显示震前的2008—2011年，西秦岭北缘断裂具有与GPS观测一致的差异运动增强特性;震前数天，位于这两条闭锁段或其端部的四店、毛羽沟短水准场地出现显著短临突变异常;震后均显示垂向差异运动呈减弱—恢复趋势。

关键词：三维负位错联合反演;岷县漳县地震;断裂闭锁;断裂活动动态变化

中图分类号：P315.725文献标识码：A文章编号：1000-0666（2019）01-0040-09

0引言

西秦岭构造区位于甘肃省的东南部、青藏高原东北缘，东昆仑断裂带东段及相连的塔藏断裂以北，从北至南由西秦岭北缘、临潭—宕昌、礼县—罗家堡、光盖山—迭山北麓和南麓、成县盆地北缘、迭部—白龙江、康县北等多条断裂构成。2013年7月22日岷县、漳县间（34.5°N，104.2°E）发生6.6级强震，打破了该区1936年6级与6.0级地震后长达77年的6级以上地震平静期。此次地震发震构造为临潭—宕昌断裂中东段的次级断层，震源深度约20km，发震断裂以逆冲为主、兼左旋走滑性质。临潭—宕昌断裂是东昆仑断裂与西秦岭北缘断裂之间的构造转换区，青藏高原东北部东昆仑断裂的向北挤压和向东运动、西秦岭北缘断裂走滑及向南北两侧逆冲“花状构造”，是临潭—宕昌断裂带构造应力集中、继而发生强震的主要动力因素（何文贵等，2013;郑文俊等，2013）。

2013年岷县漳县6.6级地震发生前，GPS连续站资料无显著异常反映（朱爽，周伟，2015）。而对1999—2007年、2009—2011年GPS区域站速度场资料的剖面计算和应变场分析表明，西秦岭构造区存在水平挤压高应变积累（王双绪等，2013），2008年汶川8.0级大震及震后调整，还使岷县漳县地震孕震区的左旋运动增强，从而加速了能量积聚（杨博等，2014）。负位错反演显示，西秦岭北缘断裂中—西段1999—2007年相对闭锁（张希等，2014a;赵静等，2016），2009—2013年包含漳县地区的断裂中段的闭锁程度增强，但GPS站点在西秦岭构造这一小范围区域仅数十个，且点距较大，针对逆冲为主的岷县漳县地震错动，更缺乏近场和垂向观测控制，反映的闭锁特性不够全面、精细，对临潭—宕昌发震构造更是没有做闭锁程度、深度等的相关反演。王双绪等（2013）综合1999—2006年西秦岭构造区水准观测所得垂直形变场图像以及GPS、重力资料总结了该次地震的孕育背景，区域水准观测从20世纪70年代始，整网复测2006年截止，沿测线点间距较小（1～5km左右），可作为三维负位错联合反演的垂向位移控制。另外，西秦岭构造区还有十几处跨断层短水准场地，始测于20世纪80年代，常规复测每年3期，测线数百米长、场地跨度不超过1km，以监视近场构造活动为目的，岷县漳县地震前短期阶段部分场地观测到突跳性前兆异常（王双绪等，2013;杨晓东等，2014;张希等，2014b）。

Matsuura等（1986）的负位错模型认为，受现今地壳运动驱动下的各个活动地块间的相对运动在地块边界处有可能受到阻碍，从而导致应力应变积聚。即视块体边界区域的地表位移为刚性块体的运动减去边界上部（由若干断裂段构成，每条断裂段用1985年Okada弹性均匀介质半空间的矩形位错模拟）对块体相对运动的锁定或部分锁定在地表产生的位移。这样，利用地表位移观测、块体划分与负位错模型，就可以通过反演获取与强震孕育有关的断裂闭锁信息。因而被许多专家应用于GPS或水准资料反演，并在模型中增加块体应变、倾斜参量或拓展至分层非均匀位错模型（伍吉仓，许才军，2002;张希等，2004，2005，2011;赵静等，2016）。GPS与区域水准有三维互补的观测意义，但负位错联合反演的研究却很少，仅张希等（2007）开展过研究，且所用资料、研究区与本文不同，更没有结合跨断层短测线即近断层观测数据。

地震是地壳、上地幔物質在构造力作用下长期变形的结果，地壳形变是反映强震孕育—发生过程的直接证据。本文针对以逆冲为主的2013年岷县漳县6.6级地震，基于西秦岭构造区区域水准观测资料，结合GPS速度场和跨断层短水准观测约束，首次尝试3种手段资料的负位错联合反演，获取孕震过程中断裂闭锁的定性定量分析结果。进而联合跨西秦岭北缘断裂、临潭—宕昌断裂的长水准剖面测线在20世纪80年代至2014年多时段速率变化数据，以及相关的跨断层短水准场地观测结果，分析岷县漳县地震前跨断裂垂直形变的动态演变以及震后变化。

1地壳运动的三维负位错联合模型

1.1资料说明与模型构建

将西秦岭构造区作为研究区，西秦岭北缘断裂、临潭—宕昌断裂作为反演对象，根据断裂（张培震等，2013）与测点分布划分块体与反演断裂段，见图1。其中，GPS资料为1999—2007年42个站点的水平运动速度场结果;区域水准资料为王双绪等（2013）所用1999—2006年垂直形变速率值，测点135个，时间仅差1年，可以构建三维联合负位错模型，获得公共时段，即1999—2006年的联合反演结果。另外，西秦岭构造区17处跨断层短水准场地中有5处场地（口子门、四店、柳家沟、奈子沟、毛羽沟）的跨断层情况、场地条件与资料连续性好，正好位于这两条待反演断裂带上，且距区域水准测线很近。1999—2006年5个场地年变速率分别为-0.13，0.04，-0.34，0.10，-0.02mm/a（负值、正值分别表示逆断、正断特性），对口子门和四店场地、奈子沟和毛羽沟场地年变速率分别求均值，结果为-0.05mm/a，0.04mm/a，加上柳家沟场地-0.34mm/a的年变速率，分别作为图1中CD、VW、EH段反演时的垂向蠕滑量约束。

F1：马衔山断裂;F2：会宁—义岗断裂;F3：拉脊山断裂;F4：青海南山—循化南山断裂;F5：西秦岭北缘断裂;F6：临潭—宕昌断裂;F7：礼县—罗家堡断裂;F8：光盖山—迭山北麓断裂;F9：光

1.2地壳运动的三维负位错联合模型

针对西秦岭构造区42个GPS站点、135个水准测点的年速率值（形变场等值线与速度矢量分布见图3a、图4），负位错模型由图1所示4个块体和5条断裂段组成。设（xi，yi）为GPS或水准测点经纬度经高斯投影所得平面坐标，位于第j（j=1，2，3，4）个块体内，Vj0x，Vj0y，Vj0z为该块体质心（xj0，yj0）的正东向、正北向及垂向平移量，εjx，εjxy，εjy，ωj，tjx，tjy为在原Matsuura负位错模型基础上借助Taylor展式增加的块体（相对质心）水平面均匀应变（x向正应变、剪应变、y向正应变、旋转）及倾斜分量（因只在曲面上观测，故仅设倾斜量），那么该测点的水平（GPS站点）与垂直（水准测点）运动观测方程可表示如下：

式中：等式左边第1，2项分别为观测与改正值;Δx=xi-xj0，Δy=yi-yj0;右边求和项的具体函数解析式见Okada（1985）的研究，Uk1，Uk2，δk，Wk分别为第k个断裂段的锁定走滑量、倾滑量、倾角、宽度。图2为Okada弹性均匀介质半空间矩形位错模型，并假设断裂延伸至地表。

将GPS水平与区域水准资料作为地表位移约束，通过贝叶斯最小二乘法（张希等，2004，2005，2007，2014a）可以反演块体及其边界断裂段的各项参数（由于地质资料已定断裂轨迹和块体边界，断裂长度、走向、中点坐标不作为参数反演）。反演时先利用只有块体运动与应变、倾斜参量的地壳运动模型，根据最小二乘法则确定等式右边前9个参数的初值及其先验误差，再通过这几项初值、断裂附近站点观测值反映的块体间相对运动及边界断裂活动状况、表1所示断裂活动信息等确定断层参数初值。实际计算中还要扣除正好位于断裂及其延伸线上的测点，这样的反演过程稳定，不会出现奇异问题。

由于块体相对运动在其边界导致的应变能积聚由负位錯描述，可利用年均负位错地震矩Mk0=μLkWk（Uk1）2+（Uk2）2来反映断裂段积累能量的大小（μ为地壳刚性系数，取4×1010N/m2），由于真实的初始状态未知，得到的是观测时段的年均增量;并用锁定位错总量（水平面走滑与法向张压、垂向三分量的平方和求根）与断裂两侧块体间差异运动三分量的平方和求根之百分比来度量闭锁程度。由块体水平面的正应变、剪应变及垂向倾斜量，还可以推算水平面的最大剪应变、面应变率以及梯度值（tjx）2+（tjy）2，以反映块体内部三维运动差异程度。

1.3模型特点与反演可靠性分析

首先，位错模型能够在三维空间描述断层的几何与运动学特征，由其发展而来的负位错模型则有反映断裂闭锁状态的特点和优势;考虑到本文研究区有数个块体和断裂段，还要综合多手段资料，故采用相对简化的Okada位错模型来反映每条断裂段的总体活动特征与能量积累状况。

其次，对GPS与区域水准观测分别进行Globk与拟稳平差后，中误差均为1mm/a量级，相差不大，考虑到中误差量值大小与平差方法甚至处理时参数设定等有关，反演前去掉误差较大的测点，再取所用观测值（GPS与区域水准）等权，能更充分、客观地利用两种手段观测资料来反映地壳三维运动状况。笔者也尝试对两种资料设权比参数，结果显示：若降低GPS观测值的权，本来点数就少的GPS资料，其参与联合反演的意义势必削弱;若增加GPS观测值的权，其拟合效果还不如等权。

反演都有多解性，本文仅反演每条断裂段的锁定走滑量、倾滑量、倾角、宽度4个参数。图3a、图4为实际观测所得垂直形变速率场（最小二乘配置拟合所得等值线分布）和水平运动矢量分布。左右旋、正逆断特性由断裂两侧实际观测反映的差异运动所定;断层倾向、倾角根据表1由程序加以控制（反演倾角不超出表1给出的倾角范围）;约束断层深度（断层宽度乘以倾角的正弦函数）小于25km（认为锁定作用主要发生在上地壳），对断层宽度也有一定约束。将跨断层短水准观测所得年变速率作为反演的垂向蠕滑量约束，即要求待定倾滑锁定参量在垂向方向的投影小于或等于块体相对垂向变化减去同时段跨断层短水准变化量，即前述CD，EH，VW段垂向活动量分别为-0.05mm/a，-0.34mm/a，0.04mm/a，所选场地跨度虽不到1km，但确定跨断层，且资料连续性好。虽不能完全排除反演多解性，至少可以获得既符合观测实际，也符合地质信息的相对合理、可靠的解，有一定意义。

图3、图4为负位错模型理论速率场与实际观测所得速率场的对比结果。尽管本文所用的Okada位错模型相对简化，但由于块体水平变形包含平移、应变及旋转多个参数，因此图4显示模型理论值与实际观测值较为接近。垂向方向只有平移和倾斜参量，实际观测值和理论值存在差异，但图3显示变化趋势基本一致。负位错反演的拟合中误差为0.85mm/a，小于GPS站点年速率值的平均中误差1.28mm/a和水准测点年速率值的平均中误差1.47mm/a，左右旋、正逆断特性与表1所示地质资料一致，结果相对合理和可靠。

（1）由表2可见，西秦岭构造区各块体均表现为水平面挤压特征，与青藏高原东北缘近NE向挤压的动力背景一致，最大剪应变率、面应变率和垂向梯度率最高的都是发生岷县漳县6.6级地震的块体3，即西秦岭北缘断裂与临潭—宕昌断裂、礼县—罗家堡断裂所围区域，揭示相对显著的差异运动特征，有利于应变积累和孕育强震（江在森等，2001）。（2）由表3可见，能量积累年增量（用负位错地震矩度量）和闭锁程度（80%以上）较大的是发震构造临潭—宕昌断裂中东段VW段及其东北侧的西秦岭北缘断裂中段CD段;西秦岭北缘断裂西段AB段也反映一定程度闭锁（闭锁程度约50%）。（3）表3显示VW段闭锁深度18.5km，与地质调查认为的震源深度（约20km）接近。

综合分析，岷县漳县地震前的1999—2006年，GPS、区域水准和跨断层短水准3种地形变手段资料的负位错联合反演显示，震区两侧的西秦岭北缘断裂中段至岷县—宕昌断裂段逆断挤压兼左旋的闭锁最为显著，反映中长期能量积累背景。赵静等（2016）利用单一GPS资料对西秦岭北缘断裂进行了负位错反演，显示1999—2007年西秦岭北缘断裂中—西段相对闭锁，闭锁程度0.6～1.0、深度不小于20km，与本文反演结果相比，闭锁深度较深，而闭锁特性、位置和程度与本文大体一致。

3长水准剖面、跨断层短水准观测反映的地震前后垂直形变尽管西秦岭构造区整网水准复测于2006年截止，但跨前述负位错反演揭示的西秦岭北缘断裂中段和临潭—宕昌断裂中东段闭锁段却有一条测线（图4中蓝点标注），距岷县漳县地震震中最近，且2007，2008，2011和2014年实施过复测，有十几个测点，2006—2008年变化不大，但2008—2011年变幅明显增大，且2008年观测时间还在汶川8级地震后。本文根据该测线上各测点距图1中西秦岭北缘断裂CD段的垂向距离排序，获得了1987—1999年，1987—2008年，1987—2011年，1987—2014年的长时间尺度的垂直形变速率变化曲线（图5a）和1987—1999年、1999—2008年、2008—2011年、2011—2014年4个主要时段的垂直形变速率动态变化（图5b）。

然而图4、5所示剖面在西秦岭北缘断裂以北和岷县—宕昌断裂段以南点位太少，无法进行负位错反演。但从图5b还是可以看出，1987—1999年、1999—2008年速率曲线在-1.6～0.6mm/a的较小范围内波动，图5a显示的1987—2008年长时间段也如此，与前述负位错反演所得闭锁特性相符;然而2008—2011年，西秦岭北缘断裂南侧20km以外加速下降，西秦岭北缘、临潭—宕昌断裂两侧差异运动明显加快，最大下降速率达6.7mm/a。尽管由于西秦岭北缘断裂的SW向倾斜和临潭—宕昌断裂的NE向倾斜，导致图5b显示出2008—2011年西秦岭北缘断裂正断加速、临潭—宕昌断裂逆断加速的特性，但这与杨博等（2014）认为的汶川地震后2009—2011年岷县漳县地震的孕震区左旋运动增强、能量积累加速的推测并不矛盾。图5b显示2011—2014年（震后一年）反向抬升、包含同震变形，图5a所示1987—2014年最长时间尺度变化则揭示断裂两侧差异运动减小，有恢复以往变化形态的趋势。陈长云和贺建明（2016）根据GPS资料分析认为西秦岭北缘断裂漳县段2009—2013年较1999—2007年滑动速率和主张应变增大，鉴于岷县漳县地震发生在2013年7月底，2009—2013年4年的变化不仅含同震变形，更应包含震前变化。因此，图5b显示的2008—2011年差异运动增强，及西秦岭北缘断裂局部段落显示正断特性也有了佐证。考虑到图5测线长仅60km，不足以反映断裂活动全貌，相对断裂单侧二三百千米长的GPS剖面（王双绪等，2013），测点密度和空间尺度差异较大，且水准测线仅穿过断裂局部段落，由于断裂展布的复杂性和可能的分段活动差异，出现局部段落正断变化也属正常;加之水准观测實施于地表，受地形地貌影响，左旋运动增强也可能引起垂向方向的差异变化，综合分析可以得出汶川地震后至2011年西秦岭北缘断裂两侧三维差异运动增大的结论。

毛羽沟、四店场地是对岷县漳县地震有较好短临前兆反映的跨断层短水准场地，震前4天这2个场地有显著正断突跳（四店场地高达5.5mm），于震后转折恢复。两处场地分别位于1999—2006年明显闭锁的发震构造南端和西秦岭北缘断裂中段（图1）。图6所示观测曲线反映的是断层下盘测点相对上盘测点的垂向变化，曲线向上为正断特性。毛羽沟场地变化形态与杨晓东（2013）、王双绪（2013）所用曲线形态相反，原因是2017年《全国跨断层场地优化改造》项目对该场地地质调查后认为以往的场地图将断裂倾向标反了。需要说明的是，跨断层短水准场地因测线和观测周期都短，多呈振荡波动状态，两处场地起测30多年来，逆断下降（尽管速度缓慢）的时段稍长，即逆断特征略占优势，与地质资料相符，周边永登5.8级、景泰5.9级地震前观测到显著的变幅加剧或尖点突跳变化。岷县漳县6.6级地震前数天出现的正断突跳与逆冲为主的地震错动不一致，可能反映临震前断层解锁甚至小范围区域的亚失稳事件或其产生的扰动影响（这种与大区域GPS挤压特性相反的正断变化也不排除与断裂展布复杂性的可能关联）。

4结论

本文基于2013年岷县漳县6.6级地震前西秦岭构造区1999—2007年GPS观测、1999—2006年区域水准观测所得水平、垂直三维速度场资料，结合跨断层观测数据，首次尝试3种手段资料的负位错联合反演，获取孕震过程的能量积累定性定量分析结果;进而联合跨西秦岭北缘断裂、临潭—宕昌断裂的长水准剖面测线20世纪80年代至2014年多时段速率变化数据，以及相关的跨断层短水准场地观测结果，分析岷县漳县地震前跨断裂垂直形变的中—短临动态演变及震后变化。

以上分析结果表明：岷县漳县6.6级地震前的1999—2006年（GPS、区域水准资料公共时间段），震区所在块体即西秦岭北缘断裂与临潭—宕昌断裂之间水平面最大剪应变率、挤压性质的面应变率、垂向梯度率都较大，与构造和动力学背景一致的三维差异运动明显，有利于能量积累和孕育强震;而西秦岭北缘断裂中—西段尤其是包含漳县地区的断裂中段、临潭—宕昌断裂中东段，即发震构造的主断裂段闭锁明显，具有挤压逆断兼左旋特性，显示中长期能量积累背景。而汶川地震后至岷县漳县地震前2年的2008—2011年，跨西秦岭北缘断裂中段漳县地区和临潭—宕昌断裂中东段岷县地区的长水准剖面显示断裂两侧差异运动增强，与GPS资料一致;震前数天，位于西秦岭北缘断裂中段闭锁段的四店短水准场地、临潭—宕昌断裂闭锁段南端的毛羽沟短水准场地甚至观测到最高达5.5mm的短临突变异常;震后跨断裂长水准剖面、短水准场地呈现的垂向差异运动均大幅减弱，呈恢复趋势。

参考文献：

陈长云，贺建明.2016.西秦岭北缘断裂分段运动变形特征分析[J].大地测量与地球动力学，36（9）：784-788.

陈长云，任金卫，孟国杰，等.2012.巴颜喀拉块体北东缘主要断裂现今活动性分析[J].大地测量与地球动力学，32（3）：27-30.

何文贵，郑文俊，王爱国，等.2013.临潭—宕昌断裂新活动特征与岷县漳县MS6.6地震关系研究[J].地震工程学报，35（4）：751-760.

江在森，丁平，王双绪，等.2001.中国西部大地形变监测与地震预测[M].北京：地震出版社.

王双绪，蒋锋云，李宁，等.2013.岷县漳县6.6级地震前区域地壳运动变形背景与断层形变异常特征[J].地震工程学报，35（3）：503-512.

伍吉仓，许才军.2002.利用GPS资料反演华北块体运动的负位错模型参数[J].武汉大学学报（信息科学版），27（4）：352-357.

杨博，朱爽，杨国华，等.2014.岷县漳县MS6.6地震前区域形变场的特征[J].地震，34（2）：115-123.

杨晓东，李宁，刘立炜，等.2014.2013年甘肅岷县漳县MS6.6地震跨断层形变异常研究[J].地震研究，37（4）：578-587.

张波，何文贵，袁道阳，等.2012.西秦岭北缘断裂带西端晚第四纪活动特征及其西延问题[J].地震，32（1）：136-143.

张培震，邓起东，张竹琪，等.2013.中国大陆的活动断裂、地震灾害及其动力过程[J].中国科学：地球科学，43（10）：1607-1620.

张希，崔笃信，郝明，等.2014a.全国主要构造区近期应变积累动态演化研究[J].地震工程学报，36（1）：61-68.

张希，江在森，王琪，等.2005.青藏块体东北缘弹性块体边界负位错反演与强震地点预测[J].地震学报，27（6）：620-629.

张希，江在森，王双绪，等.2007.青藏块体东北缘GPS与水准资料的三维负位错联合反演[J].国际地震动态，343（7）：61-66.

张希，李心然，唐红涛.2014b.青藏块体东北缘断层形变综合指标异常与岷县漳县6.6级等强震关系[J].地震研究，37（4）：572-577.

张希，王庆良，唐红涛，等.2011.汶川地震孕震背景与同震变化的铲形断层位错模拟[J].地球学报，32（2）：189-194.

张希，张四新，王双绪.2004.昆仑山口西8.1级地震前后地壳垂直运动的负位错模型[J].地震研究，27（2）：153-158.

赵静，牛安福，李强，等.2016.陇西块体周边断层闭锁程度与滑动亏损特征研究[J].地震研究，39（3）：352-358.

郑文俊，袁道阳，何文贵，等.2013.甘肃东南地区构造活动与2013年岷县—漳县MS6.6级地震孕震机制[J].地球物理学报，2013，56（12）：4058-4071.

朱爽，周伟.2015.甘肃岷县漳县6.6级地震前后区域地壳形变分析[J].地震工程学报，37（3）：731-738.

AbstractBasedonobservationaldataoftheregionallevelingandlongdistancelevelingprofilesattheWesternQinlingMountainstectonicarea，combinedwiththeGPSvelocityfieldconstraintsandthecrossfaultshortlevelingobservationaldata，weanalyzedthemiddleandlongtermlockingbackgroundoftherelatedfracturesbeforethe2013MinxianZhangxianMS6.6earthquake，anddiscussedthemiddleshorttermdynamicandaftershockvariation，aidedbythefirstassociatedthreedimensionalnegativedislocationinversion.Theresultsshowthat，themaximumshearstrainrate，compressionalsurfacestrainrateandverticalgradientrateoftheblockbetweentheWestQinlingMountainsnorthedgefaultandLintanDangchangfaultwerelargerduring1999and2007beforetheearthquake.Thethreedimensionaldifferencemovement，whichwasconsistentwiththetectonicanddynamicbackground，wasobvious，whichwasbeneficialtotheaccumulationofenergyandtheoccurrenceofstrongearthquake.ThemiddlesegmentoftheWestQinlingMountainsnorthedgefaultcontainingtheZhangxianarea，andthemiddleeasternsegmentofLintanDangchangfaultwhichisthemainfractureofseismictectonic，displayedobviouslockingfeature，andthecompressivereverseandleftlateralstrike，showedthebackgroundofthemediumandlongtermenergyaccumulation.Moreover，thelongdistancelevelingprofilefromtheZhangxiansegmentoftheWestQinlingMountainsnorthedgefaulttotheMinxianareaofthemiddleeasternsegmentoftheLintanDangchangfault，showedthatthedifferencemovementoftheWestQinlingMountainsnorthedgefaultenhancedintwoyearsbeforethisearthquake，wasconsistentwiththatoftheGPSobservations.Severaldaysbeforetheearthquake，SidianandMaoyugousiteslocatedatthesetwofaultsreflectedlargeshortimpendingsuddenlychangeabnormity;aftertheearthquake，theverticaldifferentialmovementdecreasedandrecovered.

Keywords：joint3Dnegativedislocationinversion;theMinxianZhangxianearthquake;lockingoffracture;dynamicvariationoffractureactivities