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年历史地震和周围断层的应力触发作用,得出断层滑动模型对同震计算结果的影响集中在发震断层附近而对远场影响较小.     

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

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

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地震的发生对芦山地震具有明显的促进作用, 而对康定地震的作用不明显。 在目标断层面参数已知的情况下, 根据库仑破裂应力在目标断层面上的分布, 可能为未来地震发生的地点提供线索, 进而对地震发生的危险性进行预测。 若一次地震的发生使目标断层面上有利于其错动发震的应力大面积显著增加, 这种情况下库仑破裂应力对未来地震具有预测意义。     

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

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

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

地球自转减速对2008年汶川MS8.0地震的作用

本文将地球看成是一个质量分布均匀的球对称弹性球体, 给出了地球以匀角速率旋转情况下产生的应力场的解析表达式, 计算了地球匀速自转时在2008年汶川M_S8.0地震发震断层面上产生的应力, 分析了地球自转速率变化对汶川地震发震断层的影响。 结果表明, 地球自转减速有利于汶川地震发生。 如果地球自转减速对汶川地震的发生起到了触发作用, 那么在震前震区或许会出现中小地震发生受到地球自转减速触发的现象。 分析表明, 从2006年起汶川地震震中附近地区M_L≥2.7地震受到地球自转减速的触发, 这从另一个侧面说明了地球自转减速对汶川地震的发生具有促进作用。     

强地震前后重力观测中异常变化现象的研究

分析了2001年昆仑山口西M_S8.1地震和2004年印尼M_S8.9地震前后部分重力观测的变化, 对这些采样间隔在1分钟以下的重力数据进行了高通滤波和频谱分析。 结果表明, 在昆仑山口西M_S8.1地震前, 乌鲁木齐和库尔勒观测台在二个时段分别同步出现了异常变化。 一次是在2001年11月3~5日, 另一次在2001年11月11~13日; 这二次异常都表现为颤动信号的出现, 而且后一次的强度要大于前一次; 颤动信号的周期主要集中在几十分钟。 印尼M_S8.9地震前, 这二个台也记录到类似的异常信息, 但颤动的幅度很小; 同时武汉观测台的超导重力仪也记录到了颤动信号, 该信号一直持续到地震的发生。 这些信息对于识别、 认识慢地震, 进而探索与强地震发生之间的关系具有重要意义。     

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地震对我国西藏地区有一定影响,特别是雅鲁藏布江地区和拉萨块体南北走向的正断层,库仑应力变化为正,量级可达数千帕乃至十余千帕,应该注意该区被诱发中强震的可能性.     

2022年芦山MS 6.1、马尔康MS 6.0地震井水位、水温同震响应特征

基于全国地震地下流体台网数据,分析了芦山M_S 6.1、马尔康M_S 6.0地震引起的地下流体井水位、水温同震响应特征。结果表明,对于芦山M_S 6.1地震水位同震响应观测井数量较多,以上升变化为主,同震变化幅度较大;而对于马尔康M_S 6.0地震水位同震响应观测井数量较少,以振荡为主,变化幅度较小;2次地震引起的水位同震响应能力均强于水温测项,水温记录同震响应的前提是同井能记录到水位同震变化;2组地震引起的同震响应特征差异主要与震源机制解、台站分布密度、同震响应机理不同等有关。     

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

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

2017年九寨沟地震所受历史地震黏弹性库仑应力作用及其后续对周边断层地震危险性的影响

首先对2017年九寨沟M_S7.0地震周边断裂活动和历史地震特征进行了阐述;然后利用黏弹性地壳模型,计算了1933年叠溪地震、1976年松潘震群和2008年汶川地震对2017年九寨沟地震的同震和震后库仑应力作用.该结果显示1933年叠溪地震对九寨沟地震具有延缓作用,而1976年松潘震群和2008年汶川地震对九寨沟地震的黏弹性库仑应力作用为正;随着下地壳和上地幔黏弹性物质的持续作用,前述几次地震总的黏弹性库仑应力在九寨沟地震破裂中心点处负的库仑应力逐渐减弱,而在破裂北段这些库仑应力逐渐转为正值,并促进了九寨沟地震的发生.本文也计算了九寨沟地震后对周边断层的库仑影响,并将此影响值转换为对断层能量积累的影响时间上,结果显示塔藏断裂带西段和中段在内的多条断裂带受到黏弹性库仑应力影响时间值超过10年.将库仑应力影响时间值加入到部分已知离逝时间的断层段上,也得到了这些断层段的未来30年特征地震发生概率.最终结果认为玛沁断裂带、玛曲断裂带、哈南—稻畦子断裂中段和西段等断层段的强震危险性需要重点关注.     

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主震发生后,由于同震位移的影响,不同区域处于不同的应力状态,从而使中等余震表现出不同的震源类型.     

天山北部强震活动对区域应力扰动的分析以2012年伊犁M6.6地震和2017年精河M6.6地震为例

天山地区为典型地震活跃区, 为定量分析该构造活跃区强震对周边构造变形和地震活动的影响, 本文基于地震位错理论和岩石圈分层模型计算了天山北部近期发生的2012年伊犁和2017年精河两次M6.6地震对周围地壳形变和应力的影响。 计算结果显示伊犁地震和精河M6.6地震引起震中附近地表同震位移达数厘米, 地表同震应变量级约为10^-7; 对比天山北部地区年平均构造形变特征, M6.6强震释放了震中附近近十年的构造主压应变积累; 地震引起震中附近(80 km内)同震库仑应力变化大于1 kPa, 而距震中较远区域活动断层上库仑应力变化微弱。 结合天山北部现今地壳变形特征及区域地震分布, 初步推测两次M6.6地震的发生对震后余震有显著的触发作用, 而对区域后续微震活动的影响微弱。     

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两次正断型余震.     

九寨沟M_s7.0地震的InSAR观测及发震构造分析

2017年8月8日四川省九寨沟县发生M_s7.0地震.本文基于Sentinel-1 SAR影像,利用InSAR技术获取了此次地震的同震形变场,反演获得同震滑动分布,计算了同震位错对余震分布和周边断层的静态库仑应力变化,并对发震构造进行了分析讨论.结果表明:①InSAR同震形变场显示,九寨沟地震造成地表形变最大量级约为20 cm(雷达视线方向),同震形变存在非对称性分布特征.②同震位错以左旋走滑为主,主要发生在4~16 km深度,最大滑动量约为77 cm,位于9 km深处.反演得到的矩震级为Mw6.46.同震错动未破裂到地表.③大部分余震发生在库仑应力增加区.此次地震增加了震中周边地区一些断裂的库仑应力,如东昆仑断裂带东段、龙日坝断裂、虎牙断裂等.④东昆仑断裂东段的未来地震危险性值得关注.⑤九寨沟地震的发震断层为树正断裂,可能是虎牙断裂的北西延伸隐伏部分,此次地震是巴颜喀拉块体南东向运动受到华南块体的强烈阻挡过程中发生的一次典型构造事件.