从库仑破裂应力和余震分布角度探讨汶川地震和芦山地震的关系

本文以龙门山及周边地区为研究对象,考虑区域地质构造差异、主要活动断裂带、地表附加重力影响,建立能反映地表起伏和岩石圈分层结构的龙门山地区三维粘弹性有限元模型。以地壳水平运动速率观测值为约束条件重建研究区现今构造背景应力场,在此基础上分别模拟了汶川地震和芦山地震的发生机理。通过分析同震库仑破裂应力变化与余震空间分布的关系,探讨了2次地震主震对余震的触发作用以及汶川地震对芦山地震的影响。研究表明,汶川地震和芦山地震的余震大部分由其主震触发,汶川地震对芦山地震的余震有约6.78%的触发作用。汶川地震的同震库仑破裂应力在芦山地震主震位置的增加值约为0.016MPa,如果龙门山断裂带南段库仑破裂应力年累积速率按照0.4×10-3-0.6×10-3MPa·a-1计算,汶川地震使芦山地震提前了约27-40年。计算还表明汶川地震和芦山地震的发生使鲜水河断裂带南段和虎牙断裂的库仑破裂应力增加,这些断裂带在未来发生地震的可能性增加。     

2013年芦山M_S7.0地震产生的静态库仑应力变化及其对余震空间分布的影响

2013年4月20日芦山MS7.0地震发生在龙门山断裂带的西南段,距2008年汶川MS8.0地震仅约85km,时间上仅相隔5年.首先计算了汶川地震的静态库仑应力变化对本次芦山地震的影响,得出芦山地震是由汶川地震触发造成的(库仑应力上升了0.012 MPa);进一步计算了芦山地震与汶川地震这两次大地震共同产生的静态库仑应力变化.结果表明,芦山地震的余震受前面两次大地震的共同影响,而不仅仅是芦山地震单独作用的结果,超过85%的余震发生在两次地震共同产生的静态库仑应力变化增大的地方,而芦山地震本身触发不了本次的余震序列(仅48.7%的余震位于主震所产生的应力加载区).此外,计算结果表明芦山地震本身对周边断层影响较小,仅龙门山断裂带的东北段受到一定的加载作用;而由于汶川地震的作用,安宁河断裂、大凉山断裂、马尔康断裂、岷江断裂和虎牙断裂呈卸载趋势,仅鲜水河断裂东南段和龙门山断裂中段受到加载作用,这均会加速断层上新地震的发生.     

龙门山地区历史强震数值模拟

考虑区域地质构造差异、主要活动断裂分布特征和地表附加重力影响,建立反映龙门山地区地表起伏和岩石圈分层的三维粘弹性有限元模型。以GPS为约束重建研究区现今构造应力场,依次模拟龙门山地区1900年以来发生的5次M_S 7.0以上地震,从库仑应力和等效应力角度,分析应力场演化对强震的影响以及强震间的相互作用关系。研究结果表明:从库仑应力角度,有3次地震对后续地震有促进作用,其中汶川地震对芦山地震有触发作用;从等效应力角度,有4次地震对后续地震的发生具有加速作用。     

基于库仑应力变化分析巴颜喀拉地块东端的强震相互关系

本研究基于分层黏弹介质模型,考虑同震位错效应和震后黏滞松弛效应,分析巴颜喀拉地块东端1976年松潘地震序列、2008年汶川8.0级地震、2013年芦山7.0级地震和2017年九寨沟7.0级地震等多次大地震的可能存在的触发关系,计算大地震引起的周边各主要断裂的库仑应力变化.结果显示,1976年松潘地震序列各次地震间关系密切,存在明显的相继触发作用;综合考虑同震和震后效应,汶川8.0级地震对同属于龙门山断裂带的芦山7.0级地震有触发作用,且震后效应影响不可忽略;1976年地震序列,特别是1976年8月16日7.2级地震促进了2017年8月8日九寨沟7.0级地震的发生;汶川地震对九寨沟地震的影响研究中,采用不同的汶川地震同震位错模型,计算结果有差异.综合考虑多次大地震对周边断裂带的影响,龙门山断裂带南段、鲜水河断裂带中南段、平武—青川断裂北段、灌县—安县断裂北段、文县断裂的累积库仑应力增加显著,巴颜喀拉地块东端的东昆仑断裂带东段、迭部—白龙江断裂带西段以及金沙江断裂带库仑应力亦有所增加.综合考虑各重要断裂带已有的大地震危险性分析结果和库仑应力变化计算结果,龙门山断裂带南段、鲜水河断裂带中南段、东昆仑断裂带玛沁—玛曲段和金沙江断裂带的发震紧迫性有所增强,需引起关注.     

2008年汶川地震同震-震后应力演化及其对2017年九寨沟M_S 7.0地震的影响

2017年8月8日四川九寨沟发生了MS7.0地震,距离汶川地震震中约为140km.地震的发生会调整区域应力分布和影响区域地震活动性变化.那么,汶川地震的发生对此次九寨沟地震有何影响?区域应力如何演化?为回答这些问题,本研究基于高性能并行有限元方法和网格自适应加密技术,采用含地形起伏的黏弹非均匀椭球型地球模型,分别计算了汶川地震引起的同震-震后应力时空动态演化及其对此次九寨沟地震的影响.同时,考虑同震静态应力和震后黏弹性应力调整,计算了汶川、芦山和九寨沟地震的发生对周边断层应力积累的影响.计算结果显示汶川地震对九寨沟发震断层的同震库仑应力加载约为0.008 MPa,震后9年的黏弹性应力加载约为0.012~0.016 MPa,可能会使九寨沟地震提前发生.此外,汶川、芦山和九寨沟地震的发生引起龙门山断裂映秀以南段、东昆仑断裂、龙日坝断裂东段和岷江断裂北段及鲜水河断裂康定段应力加载大于0.01 MPa,增加了这些断裂带发生地震的危险性.     

汶川地区区域构造应力场特征及汶川MS8.0地震对周围主要断层面的影响

为剖析2008年汶川M_S8.0地震对后期地震的影响及发震区域构造应力场特征,首先利用2008年5月12日—2013年4月19日汶川地震及其邻区的1660条震源机制解,同时采用同一地震多个震源机制中心解的方法筛去重复地震事件的震源机制解,最终获得911个震源机制解。其次,通过网格搜索法分段反演出区域构造应力场。结果显示:东北区主要受WNW-ESE向的挤压,西南区受W-E向的挤压,中区受WSW-ENE向的挤压。西南到东北主压应力轴方向有所变化,这可能与龙门山地区受到来自印度板块北北东方向的俯冲推挤、四川盆地的阻挡和巴颜喀拉块体东南向挤压的联合作用有关。然后,基于USGS给出的汶川M_S8.0地震的破裂模型,计算出该地震对附近强震的触发关系,结果表明,本次汶川地震对同属龙门山断裂西南端的芦山地震触发作用明显,对位于东昆仑断裂上的玛多地震也有一定的触发作用。最后,计算汶川M_S8.0地震对周围断层的同震库仑应力变化,发现本次地震造成龙门山断裂南北两端、秦岭南缘断裂、鲜水河断裂东南端、东昆仑断裂、白玉断裂的库仑破裂应力增加,龙门山断裂南北两端和秦岭南缘断裂增加最为明显,对分析地震危险性有一定的参考意义。      

MS8.0汶川地震对MS7.0芦山地震成核失稳的影响

依据速率-状态依赖性摩擦本构关系,并结合中下地壳和上地幔的黏弹性松弛效应,以震后库仑破裂应力变化和同震动态库仑破裂应力变化的计算为基础,模拟了两种应力变化对芦山地震断层的失稳发震时间的影响,研究了2008年M_S8.0汶川地震与2013年M_S7.0芦山地震之间的触发关系.计算得到,汶川地震在芦山地震震源断层面上产生的动态应力变化的峰值为0.127 MPa;此外,经过近5年中下地壳和上地幔的黏弹性松弛效应,芦山地震震源断层面上受到的震后应力变化值为0.025 MPa.结果表明,芦山地震的震源断层在应力积累逐渐接近临界状态的某一特定时期内,受到了汶川地震产生的动态应力变化、静态应力变化以及黏弹性松弛效应造成应力变化的共同触发作用,且动态应力的延迟触发作用可能更为显著.最后对芦山地震之后研究区域的应力变化场进行了初步探讨.     

2008年汶川Ms8.0地震在2013年芦山Ms7.0地震和2014年康定Ms6.3地震破裂区引起的库仑破裂应力

2008年5月12日四川汶川发生M_S8.0地震发生之后, 先后于2013年4月20日和2014年11月22日分别在汶川M_S8.0地震震中西南约80 km的芦山县和约178 km的康定县发生了M_S7.0和M_S6.3地震。 芦山地震位于龙门山前主边界处, 康定地震位于鲜水河断裂带上。 芦山地震发生前有几位作者先后计算了汶川M_S8.0地震引起的库仑破裂应力, 分析了其对周围断层的影响。 本文对这些研究结果进行了简要回顾, 并根据芦山地震和康定地震的实际发震断层面参数, 计算了汶川M_S8.0地震在芦山地震和康定地震震源深度处的水平面上以及其发震断层面上产生的库仑破裂应力, 还给出了其震中处库仑破裂应力随深度的变化。 结果表明, 汶川M_S8.0地震发生使芦山地震震源断层面上有利于其错动发震的应力大面积增加而使康定地震震源断层面上有利于其错动发震的应力大面积减小。 在芦山地震初始破裂点处引起的库仑破裂应力达0.245 MPa, 在康定地震初始破裂点处引起的库仑破裂应力为-0.00063 MPa。 因此, 汶川M_S8.0地震的发生对芦山地震具有明显的促进作用, 而对康定地震的作用不明显。 在目标断层面参数已知的情况下, 根据库仑破裂应力在目标断层面上的分布, 可能为未来地震发生的地点提供线索, 进而对地震发生的危险性进行预测。 若一次地震的发生使目标断层面上有利于其错动发震的应力大面积显著增加, 这种情况下库仑破裂应力对未来地震具有预测意义。     

利用三维孔隙弹性模型探讨紫坪铺水库对汶川地震的影响

本文基于三维孔隙弹性理论,建立了紫坪铺水库及周边地区的有限元模型.根据紫坪铺水库开始蓄水到汶川地震发震时刻的水位变化情况,计算了整个区域的孔隙压力和库仑应力.详细讨论了断层及周围地层的弹性模量和扩散系数对计算结果的影响.计算结果表明:从弹性角度看,断层的弹性模量对汶川地震震源处的库仑应力影响很小;震源处的库仑应力随着断层和周围地层的扩散系数增大而增大.当给定弹性模量和扩散系数代表性值的时候,计算结果表明在汶川地震发震时刻,震源处的库仑应力变化量为+1 kPa左右,这表明紫坪铺水库使得汶川地震发震断层更加危险.是否这个量级的库仑应力就能够触发汶川地震还需要进一步探讨.通过分析库区周边小震的分布,发现小震分布区域均是库仑应力增加的地区,因此紫坪铺水库周边的小震应该与紫坪铺水库蓄水有直接关系.     

1900年以来巴颜喀拉块体应力演化与周缘强震关系的数值模拟研究

以巴颜喀拉块体为研究对象,考虑区域地质构造差异,主要活动断裂带,活动块体和边界断裂带的划分结果,引入深部三维速度结构,建立能反映地表起伏和岩石圈分层结构的青藏高原地区三维黏弹性有限元模型.以地壳水平运动速率观测值和最大主压应力方向测量值为约束条件重建研究区现今构造背景应力场.在此基础上模拟了自1900年以来巴颜喀拉块体周缘的7级以上强震序列,从库仑破裂应力角度研究了应力演化与强震的关系、强震之间的相互作用关系以及长期构造加载对强震的影响.研究结果表明,巴颜喀拉块体周缘强震的发生可能与震源区总应力的增加有关.2008年汶川地震导致龙门山断裂带南段应力增加,表明汶川地震对2013年芦山地震有促进作用.鲜水河断裂带上的7级以上强震序列对发生在邻近龙门山断裂带上的2008年汶川地震和2013年芦山地震有延迟作用.     

鲜水河断裂带断层间相互作用的触震与缓震效应

探讨了断层间相互作用产生的同震库仑应力改变及对地震的触发与延缓效应，并以鲜水河断裂带不同断裂段时间上连续发生的4次MS6.0以上地震为例，计算和分析了每次地震发生后，在周围其它断裂上产生的同震库仑应力改变及其对后续地震的触发，以及1973年炉霍MS7.6地震和1981年道孚MS6.9地震发生后，在其周围最易破坏失稳的微破裂上产生的同震库仑应力改变及对余震活动的影响．在其它条件保持不变的情况下，将这4次地震的累积触震与缓震效应加以定量考虑，对鲜水河断裂带各断裂段的地震潜势进行了重新计算，并与已有预测结果进行对比分析，检验和评估了鲜水河断裂带断层间相互作用触震与缓震效应的重要性．结果表明:鲜水河断裂带每次地震都发生于受其先前发生的地震影响而产生同震库仑应力增加的断裂段上，不同断层间相互作用的触震和缓震效应导致地震复发概率的改变可高达30.5％以上， 主震后的余震大多发生于同震库仑应力增加较高的微破裂上．      

汶川地震和芦山地震的孕震机理及震前中长期地震危险性研究

运用非连续变形分析法与三维有限元法相结合的方法,以GPS资料作为位移速率和震源机制的约束条件,通过数值模拟研究了青藏高原及其东侧邻区构造地块的运动、变形、相互作用及其与近30年来发生于该区的大地震之间的关系。研究中引入了以应力与摩擦强度的比值定义的断层“失稳危险度”,通过数值模拟计算得到了研究区地壳块体边界断层的失稳危险度分布。结果表明,失稳危险度高的地段与近期该区发生的M_S≥7.0地震所在的位置基本一致,其中龙门山断裂带上包括汶川和芦山大地震的发震断层均为失稳危险度最高值地区。计算得到的应变率强度分布图显示,青藏高原东部边缘整条地带均为应变率强度的陡变带,特别是以龙门山断裂带上的陡变最为明显,其西侧应变率强度为东侧的近4倍,而且,这个带位于宽度相同、走向与龙门山断裂带走向相一致的高应变能密度带中,表明这两次大地震前,作为其发震断层的龙门山断裂带已积累了相当高的应变能,失稳危险度高,处于力学上的不稳定状态。模拟计算得到在上地壳层中,2001年昆仑山口西M_S8.1地震引起汶川、芦山地震发震断层的库仑破裂应力增加约0.016 MPa,相当于龙门山断裂带约两年的应力积累,也就是说,使汶川、芦山地震发震断层的失稳破裂提前了约两年。 此外,关于2008年汶川M_S8.0地震的模拟计算表明,汶川地震的发生也使包括芦山地震发震断层的龙门山断裂带西南段和东昆仑断裂带东南端的库仑破裂应力增大,应变能积累增强,这说明汶川M_S8.0地震的发生对已处于失稳危险度较高状态的2013年芦山地震和2017年九寨沟地震发震断层的提前失稳破裂起到了促进作用。      

汶川地震前后b值的时空变化及构造意义

地震b值与地应力具有密切的关系,即高应力对应低b值,研究汶川地震区b值的时空变化可以获得地震前后的应力演化过程.本文通过震源位置和速度模型联合反演的方法和双差定位法对龙门山地区汶川地震(M_s8.0)前后的地震进行了重定位工作,在此基础上精确计算了汶川地震前后不同时间段内b值大小的变化趋势及b值空间分布.结果显示:(1)b值在深度上存在明显的分层性,15 km以上b值较高,15 km以下b值较低,反映上部岩石为脆性,下部岩石为韧性的特征;(2)区域平均b值在汶川地震之后出现下降,而后有所恢复并逐渐增加,芦山地震后又降低,反映了龙门山断裂带的同震变化和震后愈合过程;(3)分时段b值空间分布结果显示,在芦山地震之后,汶川地震破裂区b值下降,即应力增加,这意味着汶川芦山地震之间的空段存在韧性变形,且芦山地震的发生促进了汶川地震破裂区的愈合;(4)映秀附近彭灌杂岩区在震前呈现低b值异常,反映该区域存在较大的差应力,可能是震后最大位错发生在该区域的根本原因.     

NUMERICAL SIMULATION OF DEFORMATION MOVEMENT AND STRESS ACCUMULATION IN LONGMENSHAN AND ITS ADJACENT FAULTS BEFORE WENCHUAN EARTHQUAKE

In order to reveal the deformation and cumulative stress state in Longmenshan and its adjacent faults before Wenchuan earthquake,a 3D viscoelastic finite element model,which includes Longmenshan,Longriba,Minjiang and Huya faults is built in this paper.Using the GPS measurement results of 1999-2004 as the boundary constraints,the deformation and movement of Longmenshan fault zone and its adjacent zones before Wenchuan earthquake are simulated.The conclusions are drawn in this paper as follows:First,velocity component parallel to Longmenshan Fault is mainly absorbed by Longriba Fault and velocity component perpendicular to the Longmenshan Fault is mainly absorbed by itself.Because of the barrier effect of Minjiang and Huya faults on the north section of Longmenshan Fault,the compression rate in the northern part of Longmenshan Fault is lower than that in the southern part.Second,extending from SW to NE direction along Longmenshan Fault,the angle between the main compressive stress and the direction of the fault changes gradually from the nearly vertical to 45 degrees. Compressive stress and shear stress accumulation rate is high in southwest segment of Longmenshan Fault and compressive stress is greater;the stress accumulation rate is low and the compressive stress is close to shear stress in the northeast segment of the fault.This is coincident with the fact that small and medium-sized earthquakes occurred frequently and seismic activity is strong in the southwest of the fault,and that there are only occasional small earthquakes and the seismic activity is weak in the northeast of the fault.It is also coincident with the rupture type of thrust and right-lateral strike-slip of the Wenchuan earthquake and thrust of the Lushan earthquake.Third,assuming that the same type and magnitude of earthquake requires the same amount of stress accumulation,the rupture of Minjiang Fault,the southern segment of Longmenshan Fault and the Huya Fault are mainly of thrust movement and the earthquake recurrence period of the three faults increases gradually.In the northern segment of Longriba Fault and Longmenshan Fault,earthquake rupture is of thrusting and right-lateral strike-slip. The earthquake recurrence period of former is shorter than the latter.In the southern segment of Longriba Fault,earthquake rupture is purely of right-lateral strike-slip,it is possible that the earthquake recurrence period on the fault is the shortest in the study region.     

2008年汶川地震前震源区的异常应力状态

我们收集了紫坪铺水库区域台网记录到的2004年8月16日至2008年5月12日汶川地震前发生在震中区及附近区域的486个小震震源机制资料,运用碎裂分析法对这些震源机制数据进行分析,获得了汶川地震前3年多时间内震源区应力水平随时间的演化过程.结果显示,强震发生前震源及附近区域约在2007年6月发生了应力突变,出现了2个应力异常高值区,其形成过程同时伴随着龙门山中央断裂的应力水平降低,最终强震发生在高应力区和低应力区之间的应力梯度带上.汶川地震前应力状态的演化过程,类似于实验研究中出现的非均匀断层失稳前的地震成核现象,符合地震成核的宏观表现.巴颜喀拉块体的南边界-鲜水河断裂的短基线测量也显示出,在同一时间段鲜水河断裂南段由左旋滑动转为异常的右旋错动,反映了汶川强震前整个巴颜喀拉块体有向东加速运动的异常过程.     

东昆仑活动断裂带大地震之间的黏弹性应力触发研究

对青藏高原北部东昆仑破裂带大地震之间的应力转移和断层相互作用进行研究. 考虑1937年以来沿此破裂带发生的5个M≥7的地震：1937年M7.5花石峡地震，1963年MS7.1都兰地震，1973年MS7.3玛尼地震，1997年MW7.5玛尼地震和2001年MW7.8可可西里地震，模拟了黏弹性成层介质中地震断层错动产生的应力演化过程，并计算了在后续地震破裂面上产生的库仑破裂应力变化. 结果表明，前面4个地震均造成2001年可可西里地震断层面上库仑破裂应力的增加，并且中地壳和下地壳的黏弹性松弛效应使得库仑破裂应力场随着时间的推移而逐渐加强. 在计算过程中定量估计了可可西里地震发生时前面4个地震同震形变和黏弹性松弛导致可可西里地震破裂面上库仑破裂应力变化之间的比值，发现前3个地震由黏弹性松弛造成的变化远远大于同震形变所造成的变化. 可可西里地震之后应力场的模拟表明东昆仑断层中段的东大滩－西大滩断层段（位于可可西里地震破裂以东及都兰地震以西）的库仑破裂应力显著增加，变化值达0.05～0.1 MPa，预示这一地区地震危险性的增加.     

龙门山断裂带上“破裂空段”的发震概率研究

本文采用分层粘弹性介质模型计算了汶川地震对芦山震中产生的库仑应力加载的影响,进而结合Dieterich(1994)提出的速率状态摩擦定律给出芦山附近区域6级地震累积发震概率随时间的变化。结果显示,2013年芦山7.0级地震时其累积发震概率达18%,说明汶川地震产生的应力扰动加速了芦山地震的发生。本文还计算了汶川、芦山2次地震对其间"破裂空段"处产生的累积库仑应力扰动的影响,结合背景地震发生率,给出了"破裂空段"处6级地震累积发震概率变化。虽然计算结果可能受到大邑地震、介质模型参数的选取和背景地震发生概率等因素影响而存在一定误差,但"破裂空段"在2次强震应力加载下累积发震概率是不断增大的,因此我们认为"破裂空段"处发生中强地震的紧迫性不断增强。     

四川芦山7.0级地震及其与汶川8.0级地震的关系

2013年4月20日在四川省雅安市芦山县发生M7.0级地震.根据四川省台网资料和收集的国内外相关资料,我们分析了芦山地震的基本参数、余震分布、序列衰减等特征.结果表明:芦山地震位于龙门山断裂南段,其震源力学机制显示为纯逆冲性质,与龙门山断裂构造特征相符合;芦山地震的余震较丰富,震后15天震区已发生7800多次余震,其中,5级以上余震4次,最大余震是4月21日17时5分芦山、邛崃交界M5.4级地震;余震分布形成的图形显示其长轴走向与龙门山断裂构造走向一致,余震分布显示密集区长轴约40 km,短轴约20 km.与汶川M8.0级地震在震源力学机制、破裂过程、余震空间展布以及地表破裂等对比分析后表明:芦山地震与汶川地震的震源错动类型、破裂过程、地表破裂以及余震活动等特征存在明显差异;芦山地震与汶川地震震中位置相距90 km,两次地震的余震密集区相距50 km;汶川8.0级地震造成龙门山断裂中北段较充分破裂,芦山7.0级地震则展布于龙门山断裂南段且破裂尺度有限;两者有发震构造上的联系,但两次地震是相对独立的地震事件.     

2013年芦山MS7.0地震和2008年汶川MS8.0地震井水位同震变化的比较分析

以芦山、汶川地震震中为中心,选择了8个井点,对2个地震引起的井水位同震变化进行了比较,重点分析了3口同震变化反向的井点,即重庆的荣昌华江井、北碚柳荫井和四川的泸沽湖井。为使对比研究结果可靠,进一步收集了2011年3月11日日本9.0级地震、2012年4月11日苏门答腊北部海域8.6级地震引起的井水位同震变化资料。对汶川、芦山地震引起的同震体应变的计算结果显示：汶川地震时,荣昌华江井、北碚柳荫井和泸沽湖井水位同震变化与井点位置处的同震体应变一致,即位于同震体应变压缩区的井水位上升,位于体应变膨胀区的井水位下降;芦山地震时,3口井的水位同震变化与同震体应变不一致,表现出与日本、苏门答腊地震时相似的同震变化特征。     

COULOMB STRESS EVOLUTION AND SEISMIC HAZARD ALONG THE EASTERN BOUNDARY OF THE SICHUAN-YUNNAN BLOCK

Using a more realistic model of multi-layered viscoelastic media, and considering the effects of the coseismic dislocation and the postseismic viscoelastic relaxation caused by the 34 great earthquakes occurring along the eastern boundary of the Sichuan-Yunnan block since 1480 and the interseismic stress accumulation caused by the tectonic loading generated by plate motions which were modeled by introducing "virtual negative displacements" along the major fault segment in the region under study, we calculated the evolution of the Coulomb stress change in each fault plane of 18 major fault segments along the eastern boundary caused by the coseismic, postseismic and interseismic effects. We studied the interactions of the Xianshuihe, Anninghe, Zemuhe and Xiaojiang fault zones on the eastern boundary of the Sichuan-Yunnan block. By evaluating if the previous earthquake could bring another earthquake closer to or farther from failure, we analyzed the interactions of the earthquakes which occurred in the different segments in the same fault zone, or in the different fault zones respectively. And further based on the calculation results of the Coulomb stress change on the fault planes, we analyzed the seismic hazard of each fault segment.The results show that the previous earthquake may trigger another earthquake which can occur in the same fault zone or in the different fault zone. And the calculation results on the evolution of the cumulative Coulomb stress change in the each fault segment show that, the Coulomb stress increases significantly in the middle section and the Moxi segment of the Xianshuihe fault zone, the Mianning-Xichang segment of the Anninghe fault zone, the Qiaojia-dongchuan segment and the Jianshui segment of the Xiaojiang fault zone, and the seismic hazard in these fault segments is worthy paying attention to.