2018年1月23日阿拉斯加Mw7.9地震发震节面讨论

2018年1月23日美国阿拉斯加湾科迪亚克岛东南280 km发生M_W 7.9级走滑地震,地震震中位于太平洋板块.震源机制解显示地震破裂面是走向东西(左旋)或走向南北(右旋)的高倾角断层,美国地震调查局(USGS)和中国地震局(CEA)分别发布了走向东西和南北的破裂模型.随着同震GPS数据的发布、余震数据的积累和对区域地质构造的了解,提供了更多资料来约束发震面,我们通过建立全球横向非均匀并行椭球型地球弹性模型,计算此次阿拉斯加湾大地震产生的同震形变、库仑应力变化以讨论上述破裂模型的差别.初步计算结果表明:USGS模型计算的最大地表水平位移约为4.5 m,CEA模型计算的最大水平位移约为3 m, USGS模型计算结果和GPS观测吻合更好;采用最优破裂面投影,计算USGS模型和CEA模型引起的库仑应力变化,分别有92.3%和72.7%的余震落在库仑应力大于0.01 MPa区域;在地质构造上,USGS给出的阿拉斯加地震节面与海底断层走向一致,但余震分布呈现出共轭条带分布,指示了地震破裂过程的复杂性.从静态应力转移角度,本次地震引起了科迪亚克破裂段(Kodiak Segment)的库仑应力增加最大0.01 MPa,加速了未来地震发生的可能性.     

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断裂系参数计算库仑应力变化时,发现新西兰北岛和南岛震中以南区域的库仑应力均增加,可触发部分余震的发生.     

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

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

基于库仑应力变化分析2015年尼泊尔MS8.1地震对中国大陆的影响

本文采用分层黏弹性介质模型, 模拟了2015年4月25日尼泊尔M_S8.1地震产生的同震和震后地表位移场, 计算了尼泊尔大地震引起的青藏高原及其周缘主要断裂上的同震和震后库仑应力变化。 地表位移场结果显示, 此次尼泊尔8.1级地震对中国大陆的影响区域主要是拉萨地块和羌塘地块, 对拉萨块体的影响主要表现为水平向南朝喜马拉雅构造带的汇聚作用, 垂直同震位移以下降为主, 震后以上升为主。 静态库仑破裂应力变化的计算结果显示, 尼泊尔大地震对青藏块体中南部的拉张性断层影响最为显著, 其中, 使尼泊尔地震北部的拉张断层的库仑应力显著增加, 个别断层库仑应力增加量超过0.01MPa, 而使其两侧的拉张断层库仑应力明显降低; 对青藏块体中部的走滑断裂则以正影响为主; 另外, 对南北地震带主要以负影响为主, 但量值微小。     

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,有可能加速强震的发生.     

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量级的区域集中在震中,但近场大部分余震分布在库仑应力减小区域,有效摩擦系数变化和区域构造应力场的耦合作用可能是其原因;利用改进的库仑应力变化计算方法和最优破裂方向计算得出的结果显示库仑应力触发理论可较好地解释余震分布.     

球体和黏弹性球体位错理论在尼泊尔Ms8.1地震中的应用与分析

本文利用球体位错理论和黏弹性球体位错理论,基于USGS发布的尼泊尔地震断层滑动模型,分别计算该地震造成的同震水平和垂直位移场,两种位错理论计算结果高度吻合,且都与实测GPS同震位移在空间分布和量级上具有较好的一致性,表明了两种位错理论的可靠性和实用性.为了更好的比较和分析这两种位错模型,分别模拟尼泊尔地震同震水平和垂直位移场,两种位错模型模拟结果均验证了尼泊尔地震主要以逆冲滑动为主,该次地震造成的水平位移较大,地震造成的南北方向上的水平位移最突出,且集中在加德满都附近区域;但模拟结果也存在差异,在近场两组结果水平位移和垂直位移差异占同震信号的不足3%,在远场两组结果的水平位移差异占信号的5%～9%左右,而垂直位移差异占信号的比例普遍在10%以上,显示出地球的黏滞性对近场水平同震位移和垂直同震位移影响较小,对远场水平位移影响有限,但是对垂直位移影响较大,即表明在计算远场同震位移时应该考虑地球黏滞性的影响.     

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

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

昆仑山口西M_S 8.1地震同震影响场的数值模拟

根据野外调查得到的地震破裂以及同震位错分布,利用三维有限元模型研究了昆仑山口西MS8.1地震的同震位移场和同震应力场的分布,讨论了震后地震活动与同震应力场变化的关系。模拟的MS8.1地震同震位移场显示,东昆仑断裂带南北两侧的同震位移衰减程度存在差异,断裂带南部的同震位移衰减大于北部同震位移的衰减;同震应力场研究结果则表明,MS8.1地震的同震最大剪切应力主要集中在东昆仑断裂带的周边地区,最大剪切应力等值线沿断裂带密集分布,形成了高梯度带。MS8.1地震的发生使青藏活动地块区的地壳应力发生了剧烈的变化,可能导致研究区的地壳应力场得到松弛。库伦破裂应力的研究表明,昆仑山地震对其后发生的5次强余震有一定的触发作用,昆仑山口西MS8.1地震可能也触发了2003年的德令哈6.6级地震     

2015尼泊尔M_S8.1地震中等余震震源机制研究

本研究利用西藏台网记录的波形数据,采用gCAP方法反演了2015年4月25日尼泊尔MS8.1大震5次中等余震(5.0≤MS≤6.5)及西藏定日MS5.9地震震源机制解.结果显示,6次地震包含2个正断、2个走滑及2个逆冲型地震.其中2个正断型地震位于主震的东北方向,即发震断层的上盘,表明该区域受到主震同震位移的影响,表现出应力拉张的变化特征;2个走滑型地震在主震破裂的东南方向上,说明随着破裂往东南方向延伸,余震的走滑分量增强;另外2个逆冲型地震位于5月12日MS7.5强余震区域,与MS7.5地震的滑移状态一致,可能与主震同震位移引起该区域处于应力挤压状态密切相关.这些结果表明,尼泊尔MS8.1主震发生后,由于同震位移的影响,不同区域处于不同的应力状态,从而使中等余震表现出不同的震源类型.     

2015年尼泊尔M_S8.1地震的地壳重力均衡背景与地表形变响应特征

对2015年尼泊尔MS8.1地震的地壳均衡背景及其引起的地表形变特征进行了研究,结果表明:(1)尼泊尔MS8.1地震震中以南的印度板块岩石圈有效弹性厚度大约为9km,加载主要来自地幔;地震以北的拉萨地块岩石圈有效弹性厚度大约为2km,加载主要来自地表.(2)尼泊尔MS8.1地震震中以南地区的地壳均衡异常大约为-100mGal(10-5 m·s-2),但其北部的地壳均衡异常则为300~400mGal,尼泊尔MS8.1地震发生在地壳均衡负异常向正异常过渡的高梯度带上.(3)尼泊尔MS8.1地震使震中周围地区的地壳整体向南运动,最大水平位移超过1.5m,分布在震中东南.震中以北的同震垂向位移总体为负值,最大下降幅度超过0.5m,同震重力变化总体为正值,最大超过60μGal(10-8 m·s-2);震中以南的垂向位移总体为正值,最大升幅超过0.7m,同震重力变化总体为负值,最大降幅超过-120μGal.(4)尼泊尔MS8.1地震使"世界屋脊"喜马拉雅山脉产生沉降,最大同震降幅超过120mm,震后松弛效应将使"世界屋脊"持续缓慢下降.该强震使世界最高峰珠穆朗玛峰降低了2~3mm,有可能被GPS、InSAR等现代大地测量工具检测到.     

芦山7级地震的同震位移估计和震源滑动模型反演尝试

利用自动经验基线校正方法,分析2013年4月20日芦山M_S7.0级地震13个近场强震动台的观测资料,以估算同震位移场分布,并据此反演了震源滑动模型.经与GPS结果比较,两种不同方法给出的芦山7级地震的水平近场同震位移场幅度都不超过cm级,均显示为典型的逆冲型地震（兼有少量左旋走滑错动）.强震最大水平和垂直永久位移分别为4.9 cm和4.4 cm,分别出现在51YAM台和51QLY 台.两种资料反演的震源滑动模型虽显示多事件特征,但主要滑动均集中在第一次事件,即初始滑动点两侧的走向长约30 km、倾向长约25 km的相对集中的较小范围内,强震和GPS模型的最大滑动量分别为1.14 m和1.09 m,较为一致.其余子事件滑动量小且分布零散,不能排除其数值效应的因素.反演矩震级均在M_w6.7左右,地表破裂应该不明显.文章还讨论了目前在我国利用近场强震动记录估计M_w6~7级地震同震位移场存在的困难和问题,为今后类似工作提供参考.     

基于无人机倾斜航空摄影三维点云测量同震倾滑变形研究

2021年5月22日青海玛多发生M_S7.4地震,震源断层错动在地表形成了长达160 km的同震地表破裂。可靠的地震地表破裂带参数是研究震源断层活动机制和评价地震危险性的重要基础。采用无人机倾斜摄影测量技术可以获得高精度的点云数据并产出DOM和DEM数据。通过跨破裂带的地形测量,获取了玛多M_S7.4地震同震地表变形的垂直位移、水平缩短量和水平拉张量等参数。测量结果显示,玛多M_S7.4地震发震断层在不同破裂段具有不同性质和大小的倾滑分量,其中具有压扭性质的野马滩观测点断层垂直位移为0.69～1.01 m,倾向水平缩短量为0.17～0.41 m,倾滑位移为0.71～1.09 m;具有张扭性质的朗玛加合日段断层垂直位移为0.34～0.54 m,倾向水平拉张量为1.99～2.08 m。     

2015年尼泊尔M_S8.1强震对中国大陆静态应力触发影响的初探

本研究基于地震危险性分析中地震重复原则和构造外推原则,以及一段时期内活动断裂上发生地震机制是相似或相近的假设,利用中国已发生的震级为MS5.0以上的地震震源机制解,将其中的一个节面视为实际地质断层参数和静态应力触发计算中的接收断层面,计算了2015年4月25日尼泊尔MS8.1地震对中国大陆地区的静态应力触发情况.研究结果表明:2015年4月25日尼泊尔MS8.1地震,触发了其后续的2次强余震的发生;本次地震对中国大陆产生的应力变化量值很小,其产生的应力加载主要集中在其邻近的西藏和新疆地区的部分断层上,而其余地区的断层则主要受到了应力卸载作用.     

基于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的库仑应力增加,对此地震有一定的触发作用.     

2010智利Maule特大地震的同震效应

2010年智利Maule大地震是人类有观测记录以来第六大地震,基于高性能并行有限元方法,建立含地表地形和Moho面起伏的大规模非均匀椭球地球模型,计算了该地震同震效应,并根据库仑应力变化分析周围断层地震活动性和主震对余震的触发关系.结果表明:对于2010年Maule地震,地球曲率和层状结构对同震水平位移影响不可忽略,如果不采用球面分层模型,在低纬度远场可能会造成位移和应变量值均被低估,在高纬度远场可能会造成位移被夸大而应变被低估.如果不考虑介质横向不均匀和地形,同震应力降被显著低估,近场同震东西向应力变化误差最高95.4%,南北向应力变化误差最高90.8%.Illapel地区库仑应力变化约10 kPa,相当于15年构造积累,Maule地震加速了2015年Illapel地震发生.以最优破裂方向投影得到库仑应力变化,经统计70.9%余震落在库仑应力变化量值大于10 kPa区域.本次特大地震触发了Pichilemu附近M_w6.9和M_w7.0两次正断型余震.