2015年智利Illapel M_W8.3地震同震效应及其对南美大陆地震危险性影响

基于高性能并行有限元方法,建立含地表地形和Moho面起伏的大规模非均匀椭球地球模型,计算了2015年9月17日智利Illapel M_W8.3地震同震效应,同时讨论2010年智利Maule M_W8.8地震和2014年智利Iquique M_W8.1地震对2015年智利Illapel地震的加载作用,进一步分析这3个地震共同作用对南美大陆及周围断层的地震危险性影响.结果表明:对此次智利Illapel地震CEA(中国地震局地球物理研究所)滑动模型计算的地表最大水平位移约为3 m,USGS(美国地调局)模型的最大水平位移约6 m.特大地震影响范围广,南美大陆接近一半区域同震水平位移量级达到0.5 mm;受3次地震共同作用,智利2010年Maule地震南部俯冲带、2015年Illapel地震和2014年Iquique地震之间俯冲带,以及2014年Iquique地震以北俯冲带,受到10 kPa东西向压应力加载;2010年Maule地震和2014年Iquique地震共同加速了2015年Illapel地震滑动;1960年Valdivia M_W9.5地震破裂区、2014年Iquique地震和1995年Antofagasta地震之间区域以及2010年Maule地震和2015年Illapel地震的未破裂区重合处库仑应力变化均超过10 kPa,有可能加速强震的发生.     

基于InSAR同震形变观测反演2010年新西兰南岛Mw7.1 Darfield地震同震破裂分布

2010年9月4日新西兰南岛Canterbury平原发生了M_w7.1地震,震源深度约为10 km.本次地震发生在一条震前不为人所知的断层上.我们利用覆盖整个震区的合成孔径雷达(SAR)观测资料,通过干涉处理分析获得雷达视线向(LOS)同震形变场;以此资料为约束反演了断层的几何参数以及同震破裂分布.结果显示,该地震造成四条相对独立断层的破裂.大部分的地震矩释放发生在Greendale断层(编号1-4),其错动以右旋走滑为主,最大破裂约为8.5 m.其它三条断层中,经过震源的逆冲断层最大破裂为5.1 m (编号6),位于Greendale断层以西的逆冲断层最大破裂为3.5 m (编号5),位于Greendale断层北面的走滑断层最大破裂为1.9 m(编号7).反演的Greendale断层地表滑动与地质调查得到的地表破裂在形态和数值上均吻合较好.本次地震释放的地震矩为5.0×10¹⁹N·m,矩震级为7.1.板块边界带形变场分析表明,Darfield地震的发生受边界带应变分配在该地区残留构造应力场控制,其复杂性体现了区域构造应力场的特点.地震对其周围地区的应力场影响较大,库仑应力增加区与余震分布有一定对应关系,并在2011年Christchurch 6.3级地震发震断层区域造成约0.1bar的库仑应力增加,对此地震有一定的触发作用.     

句容16井远大震同震效应的初步研究

通过分析江苏句容16井2001～2007年间几次远大震的同震效应现象,发现苏16井动水位和水温资料对震级大于Ms7．5,震中距800～5000km不等的远大震,有较明显的同震效应反映,且动水位和水温的同震阶变总是上升,幅值随震级的增大而增大,随震中距的增大而减小。本文对苏16井动水位和水温同震效应的同步变化机理进行了初步探讨。     

地方震及近震地下水同震震后效应研究

同震效应、震后效应是地下水异常变化的显著特征。通过多年预报研究、实践, 依据1982年以来的观测资料, 分析同震效应、震后效应的机理、特征, 提出震后效应的单井、密集、条带、迁移等空间分布, 可对老震区余震、外围区及井位附近牵动性后继地震地点预报有指示意义, 亦可对中期、短期预报及后继地震震型判定提供线索。     

北京塔院井数字化观测水温的同震效应研究

基于北京塔院井数字化水温观测资料,分析了远震引起的水温同震效应,注意到塔院井水温同震变化总是具有下降－上升－恢复的过程,不受地震方位和震源机制的影响;水温同震下降幅度随震级的增大而增大,随震中距的增大而减小, 三者之间有较好的关系;震后水温后效恢复上升幅值受水位动态影响. 最后,对塔院井水温同震效应机理进行了探讨,初步研究结果显示,井孔中的水体受震荡激发而加速对流与掺混是导致水温先下降的主要原因：当井水受到地震波的作用时,对流加速,井内深部较热的水体上涌, 而浅部较冷的水体下沉,水温探头将先观测到温度下降现象; 震后水震波逐渐平息,探头附近井水温逐步恢复上升.     

首次直接观测到与理论预测一致的同震静态应力偏量变化——2016年4月7日山西原平ML4.7地震的钻孔L应变观测

地震会引起同震位移和应力场变化,造成现存各断层上的库仑应力变化量的增加或减小并影响后续地震风险,对于地震预测具有重要意义.现有的大量地震实例中,理论计算的同震位移场与GPS、InSAR等实测的同震位移吻合,间接证实了基于均匀理想弹性假定下可以使用半无限空间Okada公式计算与位移微分相关的同震应力变化.然而,由于地壳并...     

昌马MS7.6地震三维同震形变场数值模拟及库仑应力触发研究

1932年12月25日,甘肃玉门市昌马镇附近发生M_S7.6地震(震中：39.70°N,96.70°E),此次地震是继1920年海原M_S8.5地震、1927年古浪M_S8.0地震之后青藏高原东北缘祁连山—河西走廊地震带内的又一强震。这组大地震活动无论在时间序列上还是在空间序列上都极其罕见,引起了众多国内外学者的广泛关注,普遍认为是一次陆内地震活动沿相关断裂带由东向西迁移的典型样例。基于PSGRN/PSCMP程序,以前人地表地质调查相关结果为约束,利用弹性位错理论获取昌马地震断层破裂模型,以此为基础,基于黏弹性半空间分层模型模拟计算昌马地震产生的地表同震三维形变场。通过研究单次地震与地震叠加的库仑应力变化对后续地震的触发关系,分别对海原地震、古浪地震以及它们的叠加对昌马地震的触发作用进行研究。结果表明：海原地震、古浪地震以及它们的叠加均未超过可能触发地震的经验阈值,故认为昌马地震可能不受海原地震和古浪地震所造成的库仑应力变化的影响。此研究可为利用实际形变解释昌马地震的孕震过程研究提供理论依据。     

天津井水位同震效应分析

通过整理天津张道口、静海、宝坻等多口观测井的数字化水位的观测资料,并针对全球近年来发生的8级以上地震分析了井水位同震变化的特征。结果表明,天津地区这几口观测井对远、大地震水位同震效应显著,地震中水震波变化为振荡型与震荡阶变型。     

新疆地区地下水位同震及震后效应的研究

简要回顾了目前国内同震、震后效应研究的现状,阐述了新疆的新04、10、21井同震、震后效应的形态及特征,分析并总结了新04、10、21井的同震、震后效应特征可能具有的前兆特征。     

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

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

2010智利莫尔M_W8.8地震对其破裂区及附近地区地震活动性的影响

强震的发生与地壳中应力增加有关。库仑破裂应力被用于研究由于某次强震发生所产生的应力变化。2010年智利莫尔M_W8.8地震发生之后,其破裂区以南地区没有发生M≥6.0地震,而在其以北地区于2015年在智利发生了伊亚佩尔M_W8.3级地震。在这两个地区地震活动性悬殊。本文利用库仑破裂应力研究了2010年智利莫尔M_W8.8地震在这两个地区产生的应力变化。结果表明,2010年智利莫尔M_W8.8地震在其破裂区以北地区产生的应力变化为正值,有利于断层错动发震,而在其破裂区以南地区产生的应力变化为负值,具有抑制断层错动发震的作用。     

2016年3月2日苏门答腊MS7.8地震同震位移和应力场数值模拟研究

苏门答腊地区发育多条左旋走滑性质断层,地震活动活跃.2006年3月2日该区西南海域发生了M_S7.8大地震.大地震的发生常常会引发区域位移场和应力场及周围断层应力状态发生变化.本文建立全球PREM有限元地球模型,据已有的断层滑动模型计算了此次苏门答腊地震引发的同震位移和应力及库仑应力变化,并进一步讨论了此次地震对周围断层的影响以及区域构造应力场对库仑应力变化计算的影响.初步结果表明此次苏门答腊M_S7.8地震造成较大的南北向水平位移且集中在探测者破裂区（Investigator Fracture Zone）,最大水平位移量约6.74 m,断层倾角接近垂直,下盘向北运动而上盘向南,进一步表明M_S7.8地震为典型的左旋走滑为主的地震,发生海啸的可能性较低;库仑应力变化达MPa量级的区域集中在震中,但近场大部分余震分布在库仑应力减小区域,有效摩擦系数变化和区域构造应力场的耦合作用可能是其原因;利用改进的库仑应力变化计算方法和最优破裂方向计算得出的结果显示库仑应力触发理论可较好地解释余震分布.

     

历史大地震断层滑动模型的建立及其对同震数值计算的影响——以1920年宁夏海原MS8.5大地震为例

大地震的发生会引起区域位移场和应力场发生变化,进而改变区域内及临近断层的应力状态和地震活动性.目前,研究学者可据已有的断层滑动模型来计算分析大地震同震应力变化,同时采用库仑应力触发理论来进一步分析震后余震分布和断层危险性.然而,历史上曾经发生过不少大地震,例如,1920年的海原M_S8.5大地震,是全球范围内少见的特大地震之一.局限于无确切的地震台站地震波等资料,前人在研究历史地震的影响时往往给出一些简单的断层滑动模型,将断层面上错动量视为均匀分布.为更准确地了解历史地震对后续地震的影响,基于前人研究和一般地震滑动形态分布规律及地震反射剖面等资料,以海原M_S8.5大地震为例,探讨了如何建立海原大地震断层滑动模型,并分别搭建了简单断层滑动模型和复杂断层滑动模型的全球同震横向不均匀并行椭球型地球模型.通过对海原M_S8.5地震的同震位移场和应力场的计算,发现采用复杂断层滑动模型比简单断层滑动模型地表位错分布更切合实际.同时,进一步计算和分析了此次大地震对青藏高原东北缘近100年历史地震和周围断层的应力触发作用,得出断层滑动模型对同震计算结果的影响集中在发震断层附近而对远场影响较小.

     

2008年汶川地震的破裂模型及引起的余震活动分析

用两种发震构造的破裂位移分布模型对2008年汶川地震造成的区域主要断裂的库仑应力变化进行了模拟.通过对这两种破裂位移模型模拟结果的对比分析,表明破裂位移模型的逆冲兼右旋滑动的位移模式及南北两段分别有两个滑动量高值区这样的特点控制了区域库仑应力变化的总体分布特征.汶川地震的余震分布和计算结果还显示:所获得的区域主要断裂的库仑应力变化分布和正应力变化分布同时控制了该地震发生后余震的分布,库仑应力变化增强同时断层面受到主震断层位错拉张的区域和余震的分布强烈吻合,而库仑应力变化增强但是相应部位断层面受到主震的挤压作用却没有余震发生,显示了强震后区域断层受到的正应力变化和库仑应力变化均对余震的发生起到了重要作用,而不仅仅是库仑应力变化增强是触发余震的机制;两种模型计算结果的细微差异显示了倾角向下变缓的震源破裂模型能更好地解释余震在破裂断层北部的分布特点.对昆仑山断裂、鲜水河断裂等周边强震多发断层的库仑应力变化及对未来地震活动的可能影响进行了初步评价.     

2016年11月13日新西兰凯库拉(Kaikoura)MW7.8地震同震位移和应力数值分析

2016年11月13日新西兰南岛北端凯库拉（Kaikoura）发生了M_W7.8大地震,造成了强烈的地表变形并引发大面积滑坡和海啸的发生.基于美国地质调查局（USGS）断层滑动模型,建立全球同震横向不均匀并行椭球型地球模型,计算了此次新西兰凯库拉大地震产生的同震形变和应力及库仑应力变化.初步计算结果表明：新西兰凯库拉M_W7.8地震造成断层上盘东北向抬升,下盘西南俯冲;引起发震区域同震位移较大,从凯库拉到坎贝拉（Campbell）以及首都惠灵顿（Wellington）整体上东北向抬升,最大同震水平位移1.2 m,垂直位移1.1 m.此次大地震释放了发震断层上积累的压应力,但增加了发震断层两端的挤压力;同时,同震剪应力变化增加了NE-SW向断层发生右旋滑动的危险性;采用此次地震发震断层参数计算得出的最大库仑应力变化增加区域集中在发震断层两端,可达到MPa量级.当分别采用新西兰北岛Awatere断裂系和南岛Wellington断裂系参数计算库仑应力变化时,发现新西兰北岛和南岛震中以南区域的库仑应力均增加,可触发部分余震的发生.

     

2015年4月25日尼泊尔M_S8.1大地震的同震效应

应用有限单元方法,计算了2015年尼泊尔MS8.1大地震发生产生的同震变形和应力变化.计算中考虑地球为球体以确保远场应力场变化得到可靠结果,采用PREM模型的地球分层模型,考虑了中国地震局(CEA)和美国地质调查局(USGS)各自提供的断层滑动模型.结果表明:尼泊尔MS8.1地震是一个比较典型的低角度逆冲地震,水平位移和应力降较大;地震造成南北方向上的水平位移最突出,且集中在首都加德满都附近区域.USGS断层滑动模型地表最大位移量达到3.5m,CEA滑动模型最大为1.2m;东西向和垂直方向上的同震位移相对较小;同震位移量级在0.1m的影响区域可达300km;地震造成尼泊尔地区最大库仑应力变化可达到MPa量级,地震危险性依然较大.此次MS8.1地震对我国西藏地区有一定影响,特别是雅鲁藏布江地区和拉萨块体南北走向的正断层,库仑应力变化为正,量级可达数千帕乃至十余千帕,应该注意该区被诱发中强震的可能性.     

智利M_W8.3地震与M_W8.8地震的震源过程及其引起的库仑应力特征

2015年9月17日6时54分32秒(北京时间)智利中部伊拉佩尔附近(震中31.57°S,71.67°W)发生了一次M_w8.3大地震,在此次地震震中以南约500 km处的马乌莱地区曾于2010年2月27日14时34分11秒发生过一次M_w8.8强震(震中36.12°S,72.90°W),两次地震余震分布区之间有约75 km的地震空区.本文利用远场体波与面波波形,基于有限断层模型,反演了这两次地震的震源破裂过程.结果显示这两次地震均为逆冲型大地震,2015年伊拉佩尔M_w8.3地震的平均滑动角度为107°,平均滑动量为2.43 m,平均破裂速度为1.82 km·s~(-1),标量地震矩为3.28×10~(21)Nm,95%的标量地震矩在104 s内得到了释放.最大滑动量约8 m,位于沿走向75 km,深度8 km处.2010年马乌莱M_w8.8地震的平均滑动角度为109°,平均滑动量为4.95 m,平均破裂速度1.90 km·s~(-1),标量地震矩为1.86×10~(22)Nm,95%的标量地震矩在121 s内得到了释放.最大滑动量约12.5 m,位于沿走向100 km,深度21 km处.2015年伊拉佩尔M_w8.3地震浅部更大的滑动量应该是其引起了较大海啸的一个原因.基于破裂滑动分布,我们计算了这两次地震引起的周边俯冲带上静态库仑应力变化,结果显示两次地震均显著增加了周边俯冲带上的库仑应力,2010年马乌莱地震使得2015.年伊拉佩尔地震震源区附近的库仑应力增加了(0.01~0.15)×10~5Pa,从应力积累的角度看,2010年马乌莱地震有利于2015年伊拉佩尔地震的发生,对后者的发生起到了促进作用.     

2001年昆仑山口西MS8.1地震对周围断层的应力影响数值分析

大地震的发生往往会引起周围区域形变场和应力场变化,且对临近断层上的应力状态也有影响.2001年11月4日,昆仑山口西发生了半个世纪以来中国最大的M_S8.1级地震.本文基于已有的滑动模型,建立了三维含地形高程的横向不均匀性椭球型地球有限元模型,采用等效体力方法,分析了此次M_S8.1地震产生的全球同震位移和应力场变化.与解析方法相比,该模型考虑了地形、Moho面起伏和地球介质横向不均匀性;与一般的有限元数值模拟相比,该模型考虑了地球曲率和椭率,合理地规避了有限块体模型假定边界位移为零所引入的误差.计算得出同震位移与GPS观测数据可以很好地吻合.据库仑破裂应力准则和震源参数,计算得出昆仑山口西M_S8.1地震的发生造成了汶川、芦山、改则和当雄地震的发震断层上库仑应力增加,对这些地震的发生起促进作用;而造成玉树和德令哈地震发震断层上的库仑应力变化为负值,在一定程度上抑制了这些断层的地震活动性.此外,计算结果显示地球地形高程、介质非均匀性和椭率对昆仑山口西M_S8.1地震同震变化计算有一定的影响,其中地形和椭率造成的同震位移场相对误差约10%.