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

第 12 卷 第 1 期 震灾防御技术 Vol. 12 ，No. 1 2017 年 3 月 Technology for Earthquake Disaster Prevention Mar. ，2017 刘盼，李平恩，廖力，2017 ．从库仑破裂应力和余震分布角度探讨汶川地震和芦山地震的关系．震灾防御技术，12 （1）：40—55． doi ：10.11899/zzf 从库仑破裂应力和余震分布角度探讨 汶川地震和芦山地震的关系1 刘 盼 李平恩 廖 力 （中国地震局地球物理研究所，北京 100081 ） 摘要 本文以龙门山及周边地区为研究对象，考虑区域地质构造差异、主要活动断裂带、地表附加重力影 响，建立能反映地表起伏和岩石圈分层结构的龙门山地区三维粘弹性有限元模型。以地壳水平运动速率观 测值为约束条件重建研究区现今构造背景应力场，在此基础上分别模拟了汶川地震和芦山地震的发生机理。 通过分析同震库仑破裂应力变化与余震空间分布的关系，探讨了 2 次地震主震对余震的触发作用以及汶川 地震对芦山地震的影响。研究表明，汶川地震和芦山地震的余震大部分由其主震触发，汶川地震对芦山地 震的余震有约 6.78% 的触发作用。汶川地震的同震库仑破裂应力在芦山地震主震位置的增加值约为 0.016MPa，如果龙门山断裂带南段库仑破裂应力年累积速率按照 0.4×10−3—0.6×10−3MPa ⋅a−1 计算，汶川地 震使芦山地震提前了约 27—40 年。计算还表明汶川地震和芦山地震的发生使鲜水河断裂带南段和虎牙断裂 的库仑破裂应力增加，这些断裂带在未来发生地震的可能性增加。 关键词：龙门山断裂带 库仑破裂应力 余震分布 汶川地震 芦山地震 有限元方法 引言 龙门山断裂带位于青藏高原东缘和四川盆地交汇处。在不到 5 年时间内，龙门山断裂带的 中段和南段相继发生了2008 年 5 月 12 日汶川M 8.0 地震和2013 年4 月20 日芦山M 7.0 地震， S S 给人民的生命财产造成了巨大损失。地震发生后，国内外学者从不同角度对这 2 次地震进行了 深入研究。其中，库仑破裂应力作为一个重要的指标，在分析地震之间的相互关系、探讨震后 断层应力水平的变化以及估算地震的发震周期方面具有独特的意义（Harris ，1998；Parsons 等， 2008 ；陈连旺等，2008 ；王辉等，2008 ；邵志刚等，2010 ；徐晶等，2013 ；李玉江等，2013a ； 2013b ；2014 ；Wang 等，2014 ；Jia 等，2014 ）。因此，许多学者从库仑破裂应力角度研究地震 的影响。Parsons 等（2008 ）采用地震波反演得到的均匀弹性有限断层模型，通过计算后发现 汶川地震引起的同震应力变化使龙门山断裂带南段库仑破裂应力增加了 0.1MPa，该断裂带正 是芦山地震发生的位置。Toda 等（2008 ）和万永革等（2009 ）采用弹性位错模型分别计算了汶 川地震后周围断层上库仑破裂应力的变化，认为龙门山断裂北部和最南端、鲜水河断裂带南段、 1 基金项目 地震行业科研专项（201408014 ）和中国地震局地球物理研究所基本科研业务费（DQJB14C06 ）共同资助。 ［收稿日期］ 2016-04-25 ［作者简介］刘盼，男，生于 1987 年。硕士研究生。目前从事地球动力学数值模拟方面的研究。E-mail ：liupan_0214@163.com ［通讯作者］ 李平恩，男，生于 1977 年。副研究员。目前从事地球动力学数值模拟方面的研究。E-mail ：pingen2000@163.com 1 期 刘盼等：从库仑破裂应力和余震分布角度探讨汶川地震和芦山地震的关系 41 岷江断裂的地震危险性增加。在此基础上，单斌等（2009 ）进一步考虑了地壳自重及分层结构。 董培育等（2013 ）计算了芦山地震后周边断层的库仑应力变化。单斌等（2013