2010年4月14日玉树MS7.1地震对余震的触发研究

研究了2010年4月14日青海玉树7.1级地震产生的静态库仑破裂应力变化对余震及区域地震活动空间分布的影响.在考虑震源区附近区域构造应力场的基础上,由震源区构造应力和主震破裂产生的应力叠加,计算得到最容易破裂的余震断层面,进而计算玉树主震在上述余震断层面上产生的库仑破裂静态应力变化.结果表明,库仑应力变化图像与余震分布较为吻合,说明玉树主震对大部分余震有触发作用.利用相同方法计算了玉树主震对周边更大范围内地震的应力触发影响,发现大部分区域地震活动空间分布与库仑应力变化图像基本一致,中小地震大多发生在正向触发区域,但量值较小.对余震断层面与主震完全一致的情形进行了同样的计算,与前述方法计算结果的对比显示,考虑区域应力场的计算方式所得库仑破裂应力变化图像,能够更好地与余震及区域中小地震空间分布图像相吻合.     

库仑应力计算及应用过程中若干问题的讨论——以汶川地震为例

大地震发生后,估计后续地震发展趋势是人们关心的问题.目前,学术界经常利用大地震造成的库仑应力变化研讨大震对后续地震的影响,但库仑应力的计算和应用中尚有一些问题被忽视而有待探讨.本文修正了传统的库仑应力计算中沿地震破裂面滑动方向计算剪应力变化的近似方法,考虑震后主应力方向可能改变对剪应力变化量计算的影响,对改进方法和传统方法获得的结果进行了对比分析.本文基于龙门山地区背景应力场测量资料,根据不同研究者反演的地震破裂模型,计算了汶川地震造成的静态库仑应力变化,考察了不同地震破裂模型下库仑应力分布差异.在断层面附近,使用改进方法与传统方法计算的库仑应力分布差异相当大;利用多个地震破裂模型计算大地震库仑应力判断大震后余震发展趋势时,要注意破裂模型不确定性对危险性估计的影响.     

云南澜沧-耿马Ms7.6地震的完全库仑破裂应力变化与后续地震的动态、静态应力触发

计算和研究了1988年云南澜沧-耿马Ms7.6地震产生的完全库仑破裂应力变化的时空演化图像, 对继澜沧地震13 min后发生的耿马M7.2地震的应力触发问题进行了分析, 同时也对澜沧-耿马地震后24 d内发生的Ms5.0~6.9后续强余震的应力触发问题进行了探讨. 结果显示, 澜沧地震断层破裂产生的完全库仑破裂应力变化空间分布图像具有很强的非对称性, 正值的动态和静态库仑破裂应力变化区域均具有与强余震分布位置吻合较好的现象. 耿马Ms7.2地震受到了澜沧Ms7.6地震产生动态和静态库仑破裂应力的触发作用; 绝大多数后续强余震受到了动、静态库仑破裂应力的综合触发作用.     

震源参数因素对玉树地震库仑破裂应力图像及其与余震活动统计关系的影响

2010年4月14日青海省玉树发生MS7.1强震, 根据国内外3个机构公布的5套玉树地震震源参数,计算了玉树地震产生的库仑破裂应力变化,研究了震源参数的差异对库仑破裂应力变化及其与余震活动关系的影响． 结果表明,由于不同来源震源参数存在差异,计算得到的库仑破裂应力变化图像也存在不同程度差异,从而影响库仑破裂应力变化与余震活动的关系．在震源断层参数中,断层滑动角和位错量的影响大于其它参数．此外还发现对余震统计时间的长短也会影响库仑破裂应力变化与余震活动关系的统计结果．      

2010年MS7.1级玉树地震同震库仑应力变化以及对2011年MS5.2级囊谦地震的影响

2010年4月14日青海省玉树藏族自治州发生M_S7.1级地震.和传统的板内地震相比,玉树M_S7.1级地震的余震具有数量少、震级大的特点.研究玉树地震主震与余震之间的关系,对于我们了解余震的发震机理具有十分重要的参考价值.本文利用弹性位错理论和分层岩石圈模型,计算玉树地震引起的同震及震后黏弹松弛应力场变化,讨论M_S7.1级玉树地震对余震分布的影响以及与2011年囊谦M_S5.2级地震之间的触发关系.结果显示,玉树地震导致了四处明显的库仑应力增强的扇区,2010年4月13日至6月17日的870次M_L>1.0级余震主要分布于主震破裂面附近区域以及破裂面东北端的应力增强扇区.分析玉树地震对余震分布的影响时,有效摩擦系数以及计算深度的选取对计算结果的影响较小,是否考虑区域构造应力场的影响较大.考虑区域构造应力场时,占总数86.7%的余震位于库仑应力增强区,地震应力触发理论较好地解释了余震的分布.选取囊谦地震震源机制解的两个节面作为库仑应力计算中的接收断层参数,并且考虑不同黏滞系数下的玉树地震同震及震后黏弹松弛效应,模型计算结果均表明囊谦地震位于玉树地震所导致应力影区,仅依靠地震的静态、震后黏弹松弛应力触发理论,无法解释囊谦地震的发生,说明该次地震可能是一次独立的事件.     

2006年珊溪水库ML 4.6震群静态库仑应力变化及其对后续地震的影响

浙江省珊溪水库2006年2月4日至11日发生8个M_L 4以上地震,最大地震为M_L 4.6地震,截至2009年12月10日,共发生1 063次小地震。为探讨8个M_L 4以上地震主发震断层的库仑应力变化特征以及库仑应力变化对后续地震发生的促进作用,基于库仑破裂准则,计算珊溪水库该M_L 4.6震群的静态库仑应力变化,讨论不同有效摩擦系数对库仑应力变化的影响,并依据Dieterich速率—状态摩擦定律计算并讨论其静态库伦应力变化对后续地震活动性的影响。结果表明,珊溪水库8次M_L 4.0以上地震的静态库仑应力变化及库水下渗导致的静摩擦力变化,对后续地震的发生起到一定促进作用,约76.3%的地震发生在库仑应力变化为正的区域;有效摩擦系数的取值对库仑应力场与余震分布之间的关系影响不大;计算F_11-1、F_11-2及F_11-3分支断层中部的库仑应力变化,发现库仑应力在F_11-1、F_11-2     

北天山西段中强地震应力触发作用初步研究

计算了新疆北天山西段1944年以来部分5级以上地震产生的静态库仑破裂应力变化(ΔCFS),分析该区先发生的中强震产生的累加库仑应力变化以及单个中强地震对后续地震的触发作用。研究中依据静态库仑破裂应力变化的空间分布确定后续地震活动的可能增强区范围,应用累积贝尼奥夫应变研究增强区内中等地震活动的加速矩释放特征,定量论证中强震产生的库仑应力变化对后续地震产生的影响的可靠性。结果表明,北天山西段多数中强震产生的ΔCFS有利于其邻区后续中强震的发生,其中正ΔCFS区是后续中强地震的可能发生地点,该区中强震产生的库仑应力变化作用的时空影响范围分别为几天至几十年1、00 km范围内,这为未来地震趋势判定提供了依据。     

俯冲带上特大地震静态库仑应力变化对后续余震触发效果的研究

地震静态触发研究在全世界范围内广泛开展,并取得显著成效;但是否所有的大地震,地震静态触发都有很好的效果,对此目前还不甚清楚.本文通过计算最近发生在俯冲带上的三次特大地震(2011年日本东北地震(M_w=9.1)、2010年智利地震(M_w=8.8)与2004年苏门答腊—安达曼地震(M_w=9.0))所产生的静态库仑应力变化,考察主震库仑应力变化对其后续余震空间分布的影响,从而研究俯冲带上特大地震对地震的触发效果.计算结果显示:对于2011年日本东北地震,仅有47%的后续余震发生在库仑应力增加的区域;2010年智利地震也只有47.6%的余震位置处于库仑应力变化的正值区;2004年苏门答腊地震触发了49.8%的后续余震.文中通过进一步改变模型参数(如:采用不同的有效摩擦系数,使用不同作者给出的震源模型等)进行计算,结果表明这三大地震对后续余震的触发比例仍然不高(最好情况下,触发比例也不超过60%).但对于板内地震(如:2008年汶川地震,1999年集集地震),主震对后续余震的触发比例超过85%.由此可以推知,对于俯冲带上的特大地震,地震静态触发效果不显著.因此,对于俯冲带上大地震的触发问题,还需深入研究.     

以2008年汶川8.0级地震为例分析地震触发的静态库仑应力计算的不确定性

人们常用地震产生的静态库仑应力变化解释地震的空间迁移和余震的空间分布规律,但是,在计算静态库仑应力场时,室内试验数据和观测资料往往无法对计算所涉及的参数进行足够的约束,而参数取值的不确定性会直接影响到静态库仑应力的计算结果。不同研究组得到的汶川地震产生的静态库仑应力变化结果并不一样,为研究不同计算模型和参数的选取对静态库仑应力计算的影响,我们通过建立模型并依次改变参数,分析比较了各个参数对计算结果的影响。结果表明:重力、接收断层位置和走向等的影响较小,其他参数均会对静态库仑应力变化量值和符号造成不同程度的影响,其影响在地震破裂面周边尤其明显,计算时参数敏感性分析不可忽略。     

2008年新疆于田MS7.3地震对后续 地震的完全库仑应力触发作用

在离散波数法基础上计算2008年3月21日新疆于田M_S7.3地震造成的近场区域完全库仑应力变化, 分析该变化对余震发生所产生的影响, 得到了此次地震在2008年10月5日乌恰M_S6.8地震震中处所产生的动态库仑破裂应力变化. 计算结果表明, 该地震近场区域库仑应力变化图像演化大概持续了60 s, 库仑应力变化对余震的触发率达到90%以上, 其中动态库仑应力变化图像更好地解释了余震的分布. 余震震中处的完全库仑应力变化计算结果表明, 其动态库仑应力变化远远大于静态库仑应力变化. 于田M_S7.3地震在乌恰M_S6.8地震震中处造成的最大动态库仑应力变化为0.12 MPa, 说明后者可能受到了于田M_S7.3地震的动态应力触发作用, 但不显著;而静态库仑应力则对其影响很小.      

2010年4月14日玉树Ms7.1地震对余震的触发研究

This paper calculates the static stress changes generated by the Yushu M_S 7. 1 earthquake in Qinghai Province. On the basis of regional stress,we take account of the static stress change triggered by the Yushu M_S 7. 1 earthquake to find the optimally oriented fault planes,then calculate the Coulomb stress change on the optimally oriented fault plane. The results indicate that most of the aftershocks are triggered by the mainshock. The image of Coulomb stress changes is also in accord with regional earthquakes ( M_L ≥3. 0 ) distribution,but the value is lower than 0. 01MPa. In addition,this paper calculates the Coulomb stress changes in the case that the aftershock fault plane is the same as the main shock. Through comparison,we find that the image of Coulomb stress changes obtained using the "optimally oriented fault"approach is more consistent with the distribution of Yushu aftershocks and regional earthquakes.     

2021年云南漾濞MS6.4地震序列静态应力触发研究

为深度剖析2021年云南漾濞M_S6.4地震对周围断层的影响以及地震序列间的静态库仑应力影响关系,文章基于主震的破裂模型和Okada给出位移对空间偏导数的解析式,首先计算了主震在周围断层上产生的静态库仑破裂应力,结果表明维西—乔后断裂带中段、澜沧江断裂带北端、红河断裂带北端以及怒江断裂带中段库仑应力均有千帕量级的增加。其次,计算了此次地震对周围地区产生的水平应力场及位移场,发现震中东西两侧物质向外流出,南北两侧向震中汇聚;震中南北两侧沉降,东西两侧隆升;产生的应力场呈EW向挤压,NS向拉张,在一定程度上抵消了该区域背景构造应力场。最后计算了前震-主震-余震序列间的静态库仑应力影响,结果表明前震产生的静态库仑应力促进了主震的发生;在2 km、13 km和18.5 km深度附近,触发的余震(前震和主震产生的库仑应力变化为正)比例很高,但在7 km深度处(同震破裂模型中滑移量最大)大部分余震分布在库仑应力负值区(应力影区),考虑到该深度余震与主震震源机制相差较大,因此通过模拟最易错动的断层面作为余震接受断层面,从而计算出最大静态库仑破裂应力,发现应力影区的余震仍有被触发的...     

以2008年汶川8.0级地震为例分析地震触发的静态库仑应力计算的不确定性

Static Coulomb stress change induced by earthquake slip is frequently used to explain earthquake activities and aftershock distribution.However,some parameters for the Coulomb stress calculation are unable to be well constrained from laboratory experiments and field observations.Different parameters may directly affect the pattern of static Coulomb stress.The static Coulomb stress changes induced by the Wenchuan earthquake calculated by six research groups are not consistent with each other.To investigate how the parameters affect the calculation results,we change the parameters in turn through modeling and compare the results of different calculation parameters.We find that gravity,position and strike of receiver faults have little influence on coseismic Coulomb stress calculations,but other parameters can change the value and sign of the results in various degrees especially around the earthquake rupture plane.Therefore the uncertainty analysis of static Coulomb stress change induced by earthquake should be taken into consideration in the earthquake hazard analysis.     

"Stress triggering" between different rupture events in several earthquakes

IntroductionMoreandmoreobservationssuggestthatmostintermediateandstrongearthquakeshavecomplexsourceprocesses.MostintermediateandstrongearthquakesoccurredincontinentalregionofChinaincludetwoormoreruptureevents.Thesubeventshavetimeintervalfromseveralseconds(suchasthe1990Gonghe,Qinghaiearthquakeandthe1996Lijiang,Yunnanearthquake),severalminutes(suchasthe1976Tangshanearthquake)toseveraldays(suchasthe1966Xingtaiearthquake).Tostudytherelationshipbetweendifferentruptureeventsoccurredsubsequentlyis…     

2010年4月玉树M_S7.3地震序列的断层结构

利用双差定位方法对玉树地震序列2010年4月14日至10月31日间发生的ML≥1.0地震进行双差定位,得到1545个地震的重定位结果.综合分析地震双差定位结果和玉树地震序列中强地震震源机制解,发现玉树MS7.3地震发震构造由北西向和北东东向两条相交断层组成,主震发生在北西走向的甘孜—玉树断裂带上,5月29日的MS5.9余震序列发生在北东东走向的一条隐伏断裂上,两条断裂均接近直立.甘孜—玉树断裂是羌塘地块和巴彦喀拉地块的构造边界,由于羌塘地块和巴颜喀拉地块的差异运动使甘孜—玉树断裂强耦合段应力高度积累,在应变能超过岩石强度时破裂失稳发生了MS7.3地震.主震断层的左旋滑动导致北东东向断层的正应力减小,库伦应力增加,45天后触发了MS5.9余震序列的活动.     

数值模拟静态应力扰动下的断层失稳:结果分析兼与Dieterich模型和Coulomb模型的对比

基于大量岩石力学实验,Dieterich和其他研究者(Dieterich,1978;Ruina,1983)首先提出了描述岩石摩擦过程的速率-状态摩擦定律(R-S摩擦定律).如今R-S摩擦定律已成为研究地震成核等地震演化机制的有效手段.Dieterich(1992,1994)最早提出了描述受静态剪应力扰动后断层失稳时间提前或推后的余震触发机制的解析模型.现在Dieterich模型已经成为解释余震随时间衰减规律的Omori定律等地震观测现象的有力工具.与之相对应,广泛使用的Coulomb应力失稳模型也可以给出断层受到静态剪应力扰动后,断层失稳时间的提前和推后量.Dieterich模型和Coulomb应力失稳模型基于不同的物理方法,所以在进行地震危险性评估时,二者均有各自的局限性.本文利用R-S摩擦定律控制的1-D弹簧-滑块模型,模拟计算了理论地震循环以及在不同静态剪切应力扰动下,失稳时间的提前和推后量的变化情况,然后将计算得到的时间提前和推后量分别与Dieterich模型和Coulomb模型的相应计算结果进行了定量化对比和差异性分析,并给出了相应的解释.数值模拟结果显示,对于R-S摩擦定律在参数不变的条件下,断层模型失稳时间的提前和推后量的大小强烈依赖于静态剪应力扰动的大小和作用时间,而且绝对值相同的正、负向静态剪应力扰动造成的失稳时间的提前和推后量的变化情况并不完全一致.在震后松弛/滑移阶段和闭锁阶段,时间提前和推后量是常数,且随静态剪应力扰动绝对值的增大而增大,两者的比值接近于1.0,这与Dieterich模型和Coulomb模型的结果是一致的,相应的差值小于两模型结果的10%.而在自加速阶段,模拟计算结果则存在与Dieterich模型和Coulomb模型结果不同的特征.首先,在自加速阶段模拟计算结果均偏离Coulomb模型,而且时间提前和推后量的比值小于1.0,相异于Coulomb模型的论断.不过当受到正向静态剪应力扰动后,Dieterich模型的结果和模拟计算结果是一致的,最大相差量不超过Dieterich模型结果的7%,可接近0.对于负向静态剪应力扰动,当其绝对值较小时,Dieterich模型的结果很接近模拟计算结果,相差量不超过该结果的14%.但对于绝对值较大的情况,模拟计算结果远大于Dieterich模型的结果,最大可达Dieterich模型结果的35倍,这是由于负向静态剪应力扰动后使得Ω=δ·θ/dc1的条件不再成立,进而使得Dieterich模型不再成立.总的来说,与模拟计算相比Dieterich模型可以很好地描述1-D断层受扰后失稳时间提前和推后量的变化情况,并且可以体现出正、负静态剪切应力扰动后失稳时间提前量和推后量变化的差异性,而Coulomb模型则不能完整地给出受到静态剪应力扰动后断层失稳时间提前或推后的估计值.