2017年MW6.5九寨沟地震激发的同震库仑应力变化及其对周边断层的影响

2017年8月8日21时19分在我国四川省北部阿坝州九寨沟县发生了M_W6.5左旋走滑型地震.该地震发生在青藏高原巴颜喀拉块体东北缘,东昆仑断裂南东段的塔藏断裂、岷江断裂和虎牙断裂的交汇地带.包括此次地震,近年来在巴颜喀拉块体周缘已发生了九次6级以上强震,表明巴颜喀拉块体周缘主要活动断裂上的应力水平仍处于不断调整之中.本文采用库仑应力模型研究2017年M_W6.5九寨沟地震激发的库仑应力变化、该地震与周边地震的应力触发关系以及强震对周边主要活动断裂的应力扰动.强震序列包括周边区域1536-1975年M≥6历史强震和1976-2017年的M_W≥ 6 gCMT地震目录中的强震,共计32个.研究结果表明：（1）2017年M_W6.5九寨沟地震激发的同震库仑应力变化仅在局部范围内超过0.1×10⁵Pa,且75%的余震（~12.7天）受到该地震明显的同震应力触发作用,而其余25%的余震落在应力影区,采用最优破裂面可以进一步提高同震库仑应力变化与余震分布的空间相关性;（2）2008年M_W7.9汶川地震对2017年M_W6.5九寨沟地震的发生有一定的促进作用,在后者震源处激发的同震库仑应力变化为（0.026~0.263）×10⁵Pa,震后黏弹性库仑应力变化为（0.010~0.032）×10⁵Pa.该库仑应力的变化范围取决于汶川地震源断层参数和九寨沟地震接收断层参数.2013年M_W6.6芦山地震对九寨沟地震的发生几乎没有影响（< 0.001×10⁵Pa）;（3）1654年M8.0甘肃天水南地震对九寨沟地震的发生有明显的促进作用,在九寨沟地震震源处激发的同震库仑应力变化为（0.410~1.266）×10⁵Pa,震后库仑应力变化为（0.147~0.490）×10⁵Pa.1879年M8.0甘肃武都地震可能有比1654年M8.0甘肃天水南地震更强的应力触发作用,但也有可能对九寨沟地震的发生起到抑制作用.在选取的8个九寨沟地震接收断层面上,其中6个接收断层面上该地震所激发的同震库仑应力变化为（0.913~2.364）×10⁵Pa,2个接收断层面上该地震所激发的同震库仑应力变化为（－1.326~－0.454 ）×10⁵Pa;在4个接收断层面上震后库仑应力变化为（0.094~1.072）×10⁵Pa,在另外4个接收断层面上震后库仑应力变化为（－1.593~－0.106）×10⁵Pa.1933年四川叠溪地震对九寨沟地震的发生影响较弱,其所激发的同震库仑应力变化为（0.015~0.080）×10⁵Pa,震后库仑应力变化为（－0.029~0.025）×10⁵Pa;（4）九寨沟地震仅在其附近的岷江断裂北段、塔藏断裂和虎牙断裂南段造成较明显的同震库仑应力变化,其分别为0.09×10⁵Pa、（0.14~2.03）×10⁵Pa和0.25×10⁵Pa.而进一步顾及其余31个强震的库仑应力作用则发现,同震库仑应力增加非常显著的主要活动断裂分段为：岷江断裂北段南侧和岷江断裂南段的库仑应力变化分别升高5.6×10⁵Pa和9.8×10⁵Pa.鲜水河断裂北段南侧库仑应力升高23.0×10⁵Pa,鲜水河断裂南段道孚-康定段的北侧库仑应力升高9.0×10⁵Pa,而最南端库仑应力升高3.0×10⁵Pa;龙门山断裂带中段的北侧库仑应力变化为（6.1~7.4）×10⁵Pa,中段库仑应力增加（2.1~11.5）×10⁵Pa;西秦岭北缘断裂东段库仑应力变化为4.4×10⁵Pa;龙日坝断裂北段最北侧的库仑应力变化为2.0×10⁵Pa;小金河断裂北段库仑应力变化为1.7×10⁵Pa;安宁河断裂北段库仑应力变化为1.6×10⁵Pa;（5）由于下地壳和上地幔的黏弹性松弛作用,所有强震在九寨沟地震震后20年造成的黏弹性库仑应力变化在鲜水河断裂、龙门山断裂中段、塔藏断裂以及秦岭南缘断裂西段比较显著,其分别为：（1.0~3.0）×10⁵Pa、2.8×10⁵Pa、（2.3~2.7）×10⁵Pa和0.9×10⁵Pa.但总体上黏弹性库仑应力变化没有改变各断裂上的同震库仑应力变化空间分布.总的库仑应力变化在鲜水河断裂北段南侧和南段的道孚至康定段北侧、龙门山断裂中段北侧、岷江断裂南段和北段南侧、虎牙断裂、塔藏断裂以及西秦岭北缘东段很显著（均超过4×10⁵Pa）.由于库仑应力明显升高可能预示着地震潜在危险性增强,因此这些断裂分段可能将来需要重点加以关注.     

基于库仑应力变化分析巴颜喀拉地块东端的强震相互关系

本研究基于分层黏弹介质模型,考虑同震位错效应和震后黏滞松弛效应,分析巴颜喀拉地块东端1976年松潘地震序列、2008年汶川8.0级地震、2013年芦山7.0级地震和2017年九寨沟7.0级地震等多次大地震的可能存在的触发关系,计算大地震引起的周边各主要断裂的库仑应力变化.结果显示,1976年松潘地震序列各次地震间关系密切,存在明显的相继触发作用;综合考虑同震和震后效应,汶川8.0级地震对同属于龙门山断裂带的芦山7.0级地震有触发作用,且震后效应影响不可忽略;1976年地震序列,特别是1976年8月16日7.2级地震促进了2017年8月8日九寨沟7.0级地震的发生;汶川地震对九寨沟地震的影响研究中,采用不同的汶川地震同震位错模型,计算结果有差异.综合考虑多次大地震对周边断裂带的影响,龙门山断裂带南段、鲜水河断裂带中南段、平武-青川断裂北段、灌县-安县断裂北段、文县断裂的累积库仑应力增加显著,巴颜喀拉地块东端的东昆仑断裂带东段、迭部-白龙江断裂带西段以及金沙江断裂带库仑应力亦有所增加.综合考虑各重要断裂带已有的大地震危险性分析结果和库仑应力变化计算结果,龙门山断裂带南段、鲜水河断裂带中南段、东昆仑断裂带玛沁-玛曲段和金沙江断裂带的发震紧迫性有所增强,需引起关注.

     

九寨沟M_S7.0地震对巴颜喀拉块体东北缘活动断裂影响的有限元模拟

在详细调研地震地质资料的基础上,构建了巴颜喀拉地块东北缘三维有限元模型。以九寨沟M_S7.0地震同震位错为荷载,模拟计算了九寨沟地震的发生对巴颜喀拉块体东北缘主要活动断裂加卸载效应的影响。模拟结果显示,九寨沟地震的发生对龙日坝断裂、虎牙断裂、青川-平武断裂西段、迭部-白龙江断裂西段和东段、临潭-宕昌断裂东段,以及处于甘青川交界危险区内的东昆仑断裂东段、塔藏断裂西段,处于六盘山南-西秦岭东危险区的西秦岭北缘断裂东段表现为库仑应力加载;对岷江断裂、塔藏断裂东段库仑应力卸载效应显著。     

松潘-甘孜块体东北端强震间相互作用及地震危险性研究

松潘-甘孜块体位于中国大陆西南部、南北地震带的中段,其东段与扬子块体相接,拥有多条活动断裂带,是青藏高原北部的主要构造单元.该地区地震活动性强烈,历史上曾发生过多次灾难性地震.本文基于地震触发原理和黏弹松弛分层地壳模型,计算了松潘-甘孜块体东北端历史强震之间应力传输和相互作用的过程.模型结果显示,受之前地震导致的库仑应力场变化的影响,1879年武都地震和1976年8月23日松潘M7.2级地震震中库仑应力积累提升,将促进这些地震提前发生;1933年M7.5叠溪地震和1973年M6.5松潘地震震中库仑应力降低,前续地震的影响可能使得这两次地震的发震时间推迟;在研究历史地震对1960年漳腊M6.7级地震、1976年8月16日M7.2级和1976年8月22日M6.7级松潘地震的作用时,有效摩擦系数的取值十分重要,当有效摩擦系数取0.8时,前续地震导致的应力场变化将促进以上三次地震的发生.松潘-甘孜块体东北端的强震活动有效地增强了西秦岭北缘断裂、东昆仑断裂玛沁-玛曲段、鲜水河断裂康定-道孚段和岷江断裂中段上的库仑应力积累,将提升这些断裂今后发生地震的概率;有效降低了龙日坝断裂上库仑应力的积累,降低了该断层上发生地震的概率.松潘-甘孜块体的地震活动降低了汶川地震震中位置的库仑破裂应力,但提升了破裂面东北段的应力积累,有助于汶川地震向东北端破裂.     

松潘-甘孜块体东北端强震间相互作用及地震危险性研究

松潘-甘孜块体位于中国大陆西南部、南北地震带的中段,其东段与扬子块体相接,拥有多条活动断裂带,是青藏高原北部的主要构造单元.该地区地震活动性强烈,历史上曾发生过多次灾难性地震.本文基于地震触发原理和黏弹松弛分层地壳模型,计算了松潘-甘孜块体东北端历史强震之间应力传输和相互作用的过程.模型结果显示,受之前地震导致的库仑应力场变化的影响,1879年武都地震和1976年8月23日松潘M7.2级地震震中库仑应力积累提升,将促进这些地震提前发生;1933年M7.5叠溪地震和1973年M6.5松潘地震震中库仑应力降低,前续地震的影响可能使得这两次地震的发震时间推迟;在研究历史地震对1960年漳腊M6.7级地震、1976年8月16日M7.2级和1976年8月22日M6.7级松潘地震的作用时,有效摩擦系数的取值十分重要,当有效摩擦系数取0.8时,前续地震导致的应力场变化将促进以上三次地震的发生.松潘-甘孜块体东北端的强震活动有效地增强了西秦岭北缘断裂、东昆仑断裂玛沁-玛曲段、鲜水河断裂康定-道孚段和岷江断裂中段上的库仑应力积累,将提升这些断裂今后发生地震的概率;有效降低了龙日坝断裂上库仑应力的积累,降低了该断层上发生地震的概率.松潘-甘孜块体的地震活动降低了汶川地震震中位置的库仑破裂应力,但提升了破裂面东北段的应力积累,有助于汶川地震向东北端破裂.     

Earthquake stress interaction in the northeastern Songpan-Garz(e) block and its implication for earthquake hazard

松潘—甘孜块体位于中国大陆西南部、南北地震带的中段,其东段与扬子块体相接,拥有多条活动断裂带,是青藏高原北部的主要构造单元.该地区地震活动性强烈,历史上曾发生过多次灾难性地震.本文基于地震触发原理和黏弹松弛分层地壳模型,计算了松潘—甘孜块体东北端历史强震之间应力传输和相互作用的过程.模型结果显示,受之前地震导致的库仑应力场变化的影响,1879年武都地震和1976年8月23日松潘M7.2级地震震中库仑应力积累提升,将促进这些地震提前发生;1933年M7.5叠溪地震和1973年M6.5松潘地震震中库仑应力降低,前续地震的影响可能使得这两次地震的发震时间推迟;在研究历史地震对1960年漳腊M6.7级地震、1976年8月16日M7.2级和1976年8月22日M6.7级松潘地震的作用时,有效摩擦系数的取值十分重要,当有效摩擦系数取0.8时,前续地震导致的应力场变化将促进以上三次地震的发生.松潘—甘孜块体东北端的强震活动有效地增强了西秦岭北缘断裂、东昆仑断裂玛沁—玛曲段、鲜水河断裂康定—道孚段和岷江断裂中段上的库仑应力积累,将提升这些断裂今后发生地震的概率;有效降低了龙日坝断裂上库仑应力的积累,降低了该断层上发生地震的概率.松潘—甘孜块体的地震活动降低了汶川地震震中位置的库仑破裂应力,但提升了破裂面东北段的应力积累,有助于汶川地震向东北端破裂.     

东昆仑断裂滑动速率变化及其对2017年九寨沟地震的应力加载

为了解东昆仑断裂活动对2017年8月8日九寨沟M_S7.0地震的影响,本文选取1999-2007年、2013-2017年GPS速度场作为约束,基于块体-位错模型反演计算东昆仑断裂两个时间段的块体运动速率、断裂滑动速率和滑动亏损率,并进一步研究青藏高原东缘最大剪应变率场和九寨沟震区的震间库仑应力累积速率.结果显示,东昆仑断裂中西段左旋走滑速率较高,东段走滑速率较低,自西向东逐步递减,存在明显的梯度.在两个时间段,阿坝块体刚性运动的方向顺时针偏转0.2°,运动速率由12.22 mm·a^-1增大到15.96 mm·a^-1;东昆仑断裂左旋走滑速率升高,其中西段较为明显（升高约1.2±0.3 mm·a^-1）;东昆仑断裂东段闭锁深度和闭锁程度增加;2013-2017年,东昆仑断裂滑动引起的九寨沟震区库仑应力累积速率是1999-2007年的3倍,最大剪应变率也明显升高.因此本文认为：2008年汶川地震和2013年芦山地震后,龙门山断裂部分解锁,阿坝地块活动性增强,东昆仑断裂滑动速率增大,导致九寨沟震区库仑应力加载速率增加,加速了九寨沟地震的孕育过程.

     

2008年汶川MS8.0地震对周边断层地震活动的影响

为分析2008年5月12日四川汶川M_S8.0级地震对周边断层地震活动的影响,本文首先基于Burgers体黏滞松弛模型计算汶川M_S8.0级地震引起的库仑应力动态演化,分析认为2008年汶川M_S8.0级地震在周边断层上引起的库仑应力显著增加的主要有四个断层段,分别为鲜水河断裂道孚-康定段、东昆仑断裂东段玛曲段、青川断裂和龙门山断裂南段.而且震后4年内黏滞松弛引起的库仑应力变化量可能与同震变化相当,相当于再发生一次汶川地震所造成的影响,因此震后效应在分析强震影响时不应忽略.本文基于强震引起的库仑应力变化动态演化,结合背景地震发生率、由Dieterich(1994)模型给出地震发生概率,结合相关构造地质、历史地震、余震活动等方面资料的综合分析认为,上述4个断裂段地震危险性由高到低依次为鲜水河断裂道孚-康定段、龙门山断裂南段、东昆仑断裂东段玛曲段和青川断裂.     

2008年新疆于田7.3级地震对周围断层的影响及其正断层机制的区域构造解释

计算了2008年3月21日新疆于田Ms7.3地震导致周围断层的库仑破裂应力变化. 结果表明,库仑应力在贡嘎错断裂的西南段和东北段、康西瓦断裂东段和西段、普鲁断裂西段、龙木错-邦达错西段、阿尔金断裂和西昆仑山前逆冲断裂增加分别达6.9×104 Pa、4.4×103 Pa、8.0×105 Pa、2.8×103 Pa、2.3×104 Pa、5.7×103 Pa、1.6×103 Pa和1.1×103 Pa,而使得贡嘎错断裂中段、康西瓦断裂中段、龙木错-邦达错东段、普鲁断裂中段和喀喇昆仑北段断层上的库伦破裂应力减小分别达3.9×105 Pa、5.8×105 Pa、1.5×104 Pa、1.5×104 Pa和1.1×103 Pa. 根据康西瓦的断层滑动速率估计该地震使得中段特征地震延迟1~127年,使得东段特征地震提前2~138年. 计算还得到2005年10月8日巴基斯坦Mw7.6地震在于田地震破裂面和滑动方向上产生的库仑破裂应力增加为10~20 Pa,虽然有利于于田地震的发生,但数值太小,难以看到明显的触发作用. 分析了该地震及其周围已发生地震出现正断层机制的原因. 认为在印度板块向欧亚板块北北东向推挤作用下,相对塔里木块体,青藏高原北缘的柴达木块体沿阿尔金断裂东向运动,而兴都库什块体沿喀喇昆仑断裂整体北西向运动. 地震发生区位于此两块体之间,在双边动力相互作用下导致了东西向拉张的形变模式.     

九寨沟M_S7.0地震静态库仑应力触发及对附近断层影响分析

计算和研究了1973~2013年巴颜喀拉块体东边界发生的6次M_S≥6.5地震引起的静态库仑应力变化图像,认为先前发生的6次地震对2017年8月8日九寨沟M_S7.0地震具有触发作用。九寨沟M_S7.0地震发生后,对震中附近主要断裂应力转移的计算结果显示,在巴颜喀拉块体东边界7次M_S≥6.5地震的共同作用下,东昆仑断裂带东段及龙日坝断裂带库仑应力水平升高,具有触发作用;岷江断裂带处于库仑应力影区内,具有减震效应。     

1997年以来巴颜喀拉块体周缘强震之间的黏弹性触发研究

2021年5月22日在青海玛多发生M_W7.4地震,为了探究玛多地震的不同滑动模型对周围地区及断层应力的加卸载作用,本文首先以GNSS数据为约束,结合中国地震局地质研究所公布的玛多地震同震滑动模型（模型A）断层面几何结构反演获得同震滑动模型（模型C）,再分别利用模型A、模型B（USGS）、模型C计算玛多地震对周围地区及断层的应力加卸载作用。结果显示：① 模型C矩震级为M_W7.46,最大滑动量为3.39 m,主体破裂位于0—10 km深度范围,整体破裂东侧大于西侧,滑动分布相对于模型A也更加均匀平滑,反演效果较好;② 不同模型计算的应力分布基本相同,沿破裂段同震库仑应力加载区域面积随着深度的增加而增加,且在发震断裂带西端、东端分别各有两处明显的库仑应力加载区域,在昆仑山口—江错断裂东、西段、甘孜—玉树断裂、东昆仑断裂东段、玛多—甘德断裂、清水河断裂中、西段、达日断裂西段均产生了明显的应力加载,但模型B计算结果有所差异,昆仑山口—江错断裂中段处于应力卸载状态,震后10年断层应力状态变化不大,但清水河断裂东段在震后应力调整中卸载作用较为明显,地震危险性降低;③ 为了探究强震对玛多地震的影响,本文分别计算了2008年汶川地震以后巴颜喀拉地块M≥7.0强震同震及震后效应对玛多地震的应力加卸载,结果表明所有强震均对玛多地震有应力加载作用,但累积库仑应力并未超过触发阈值。

     

青藏高原东北缘主要断裂带的库仑应力演化及其强震危险性

自1920年海原发生M8.5地震以来,青藏高原东北缘接连发生了1927年古浪M8.0地震、1932年昌马M7.6地震等一系列大地震,使其进入了强震活动的丛集期。为了探究青藏高原东北缘这一系列地震间的相互作用及区域地震危险性,建立青藏高原东北缘的三维Maxwell黏弹性有限元模型,模拟了区域自1920年以来17次M6.7以上地震的同震及震后库仑应力演化。结果显示：研究区自1920年海原M8.5大地震之后,后续的16次地震中,有13次地震发生在库仑应力变化为正的区域,说明了地震间的相互作用可能是导致区域地震丛集的主要原因之一。系列地震发生后,阿尔金断裂、柴达木盆地断裂西段、东昆仑断裂中段、鄂拉山断裂北段、共和盆地断裂南段、日月山断裂南段、庄浪河断裂、礼县—罗家堡断裂、成县盆地断裂西段、文县断裂西段、龙首山断裂南段、六盘山断裂东段、西秦岭北缘断裂东段、海原断裂西段和祁连断裂东段位于库仑应力变化为正的区域,且大部分断裂或断裂段的累积库仑应力变化超过了0.01 MPa,它们未来的地震危险性较高。     

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

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

基于GNSS的青藏高原东北缘地壳运动场及强震趋势研究

利用“中国大陆构造环境监测网络”GNSS数据研究1998—2018年青藏高原东北缘排除同震影响等干扰后的速度场、主应变率场、最大剪切应变率场、面应变场等的变化,活动断裂滑动速率变化、跨活动断裂基线变化等。将研究区域内的二级块体再分区,获得各次级块体内部的应变率变化;获取研究区域地壳运动场的趋势性、动态特征。研究结果显示,阿尔金断裂带中东段、祁连块体和柴达木块体交界、巴颜喀拉块体与羌塘块体交界、祁连块体南边界中段、海原—六盘山断裂带和西秦岭北缘断裂带西段的逆冲运动,祁连块体北边界西段、庄浪河断裂的左旋走滑运动,祁连块体北边界东段、西秦岭北缘断裂带东段的左旋逆走滑运动,都属于造成一定程度地壳变形的持续性局部应变增强活动。阿尔金断裂带东段、东昆仑断裂带中西段、祁连块体北边界、庄浪河断裂北段、海原断裂南段、六盘山断裂北段、西秦岭北缘断裂带东段可能存在闭锁,未来十年可能发生M_S6.0以上地震。     

巴颜喀拉块体周缘强震间应力作用与丛集活动特征初步分析

巴颜喀拉块体是近年来青藏高原内部板内强震高发地带,对于该块体强震发生时空模式的研究可为大范围板内地震的深入研究提供基础;而断层间的相互作用在强震发生模式中的意义也需要深入研究。文中使用了黏弹性库仑应力模型,计算了1893年以来巴颜喀拉块体边缘和内部7级左右及以上历史地震破裂之间的黏弹性库仑应力相互影响,并使用强震所在断层的滑动速率将这一应力影响值转换为影响时间,并在历史地震发生时间中减去该影响值,比较历史地震影响值对于丛集现象的影响情况,从而分析历史地震间的相互作用对于巴颜喀拉块体内部历史地震丛集特征的影响。从结果看,巴颜喀拉块体的强震自1893年以来存在着前期相隔16a左右的准周期活动,然后经过相对较长时间平静后大约在1997年后发生丛集现象,目前巴颜喀拉块体仍然位于这一强震丛集活动中,仍存在着发生强震的危险;计算出的未来30a强震发生条件概率显示,危险程度较高的断层包括玛沁断裂、玛曲断裂、阿万仓断裂、塔藏断裂罗叉段、鲜水河断裂带磨西段、当江断裂,其他破裂段也存在发生MS7.0左右强震的危险。     

构造转换区强震间相互作用机制的数值模拟研究: 以2008—2020年4次MS＞6.0于田地震为例

2008年于田M_S7.3地震后, 阿尔金断裂带西南段的构造转换区内接连发生了2012年M_S6.2、2014年M_S7.3和2020年M_S6.4等强震.研究强震间相互作用对于理解强震的发震机制具有重要意义.本文通过构建阿尔金断裂西南段及邻区三维黏弹性有限元模型, 以强震同震滑动模型为约束, 综合同震位错与震后黏弹性松弛效应的影响, 研究4次强震间相互作用以及强震发生所引起的区域内断裂库仑应力变化.计算结果表明: 考虑同震和震后效应的共同作用, 后续3次地震发生前震源处库仑应力变化分别为-1.89×10⁴ Pa、3.63×10⁴ Pa和1.04×10⁴ Pa, 表明2008年于田地震延迟了2012年于田地震的发生, 2008和2012年于田地震共同触发了2014年于田地震的发生, 而前3次地震同样触发了2020年于田地震的发生.从不同构造部位强震间相互作用的结果来看, 阿尔金断裂西南段尾端张性区内发生的正断型强震, 对区域内北东向断层上的左旋走滑型强震有加载效应, 造成强震的提前发生, 同时也会延迟尾端张性区内近南北向的断层上正断型强震的发生.综合分析认为, 阿尔金断裂和郭扎措断裂西段应力显著增加(＞1.0×10⁵ Pa), 未来发生强震的可能性很高, 需要重点关注.

     

2021年青海玛多MS7.4地震对周围地区的应力影响

为了探讨 2021 年青海玛多MS7.4地震对后续地震的影响,基于 2021 年青海玛多MS7.4 地震的破裂模型及均匀弹性半空间模型,本文计算了该地震在震中附近主要断层面上产生的同震库仑应力变化,结果表明库仑破裂应力加载区主要集中在东昆仑断裂东段西部、玛多—甘德断裂中段、昆仑山口—江错断裂西段,而库仑破裂应力影区主要集中在东昆仑断裂中段东部、昆仑山口—江错断裂东段、玛多—甘德断裂西段、达日断裂.其次,本文计算了该地震对周围地区造成的水平面应力变化及位移场.位移场的水平分量表明震中西南和东北两侧物质向震中汇聚,而震中东南和西北两侧物质向外流出,从位移场的垂直分量来看,震中西南和东北两侧表现出明显的隆升,而在震中东南和西北两侧表现出明显的沉降.水平面应力在震中的北东-南西两侧(东昆仑—柴达木断块东侧和巴颜喀拉断块西侧)增加(拉张),而在震中的北西-南东两侧(东昆仑—柴达木断块西部和巴颜喀拉断块东部)降低(压缩),在昆仑山口—江错断裂段附近,水平面最小主应力近 NW向,水平面最大主应力近 SW向,在一定程度上抵消了该区域构造应力场(SW向挤压,NW向拉张),结合位移场说明该次地震是在构造应力场作用下的一次正常应力释放.从整体来看,水平面最小主应力和水平面最大主应力的方向分布类似于磁场线的分布,并且两者互相垂直.     

2001年昆仑山地震同震—震后效应与2021年玛多地震关系探讨

位于青藏高原中北部的巴颜喀拉地块是我国西部近年来的主体地震活动区,一系列M_W7.0以上强震均发生在该次级块体周边,而其北边界东昆仑断裂带是一条长达2000 km、规模最大、活动性最强的深大断裂带.2001年在东昆仑断裂带中段发生了M_W7.8昆仑山地震,2021年5月在其震中东南部大约450 km处巴颜喀拉块体内部一次级断裂上发生了M_W7.3玛多地震.玛多地震对人们以往认为强震更可能发生在巴颜喀拉块体边界断裂上的认识提出挑战,但是也为研究巴颜喀拉块体边界断裂与块体内部次级断裂活动关系、地震触发关系带来机遇.本文利用前期基于2001年昆仑山地震后积累的大量InSAR数据获得的震后大范围形变场时空演化图像和库仑应力变化模型,探讨昆仑山地震与玛多地震的关系.InSAR震后观测结果显示：昆仑山地震后沿东昆仑断裂带出现了长达500 km的大范围南北不对称震后形变场,其中南盘形变宽度和量级均明显大于北盘,南盘形变宽度达到250 km,断层近场相对平均形变速率达到＞20 mm·a^-1,而且南盘向南衰减梯度小,整体衰减缓慢,意味着震后形变对巴颜喀拉块体形成持续东向加载作用,并将分摊到块体内部的一系列次级断裂上,应力加载增加次级断层的地震危险性.2015—2020年InSAR震间应变率场则显示次级断裂——昆仑山口—江错断裂呈高剪切应变率特征.本文计算了昆仑山地震同震破裂和震后形变引起的玛多震区多条SE向次级断裂的累积库仑应力变化,结果显示昆仑山地震同震和震后形变对玛多地震发震断裂(昆仑山口—江错断裂)形成了一定的应力加载.本文认为昆仑山地震同震和长时间尺度震后形变加速了巴颜喀拉块体的东向运动,而断层本身运动学性质和区域应力扰动共同影响了玛多地震的发生.

     

巴颜喀拉块体东南缘应力演化及鲁甸地震自发破裂模拟研究

巴颜喀拉块体活动导致的昆仑—汶川地震系列,是目前中国大陆唯一一组7级以上地震序列,未来几年该块体仍是发生7级地震的主体地区。巴颜喀拉块体东缘从走滑运动向逆冲运动转换的构造动力学性质使得活动块体两个相邻边界带存在大地震序列的关联性。研究表明近期逐渐加速的块体东边界强震序列是北边界强震活动的响应,巴颜喀拉块体南边界的鲜水河断裂带同样以左旋走滑的方式向东边界挤压,但巴颜喀拉块体东南缘地震呈现以下两个特点,一是强震大都发生在巴颜喀拉块体东边界断裂带,东边界强震序列呈现加速状态;二是南边界鲜水河断裂带在1973年炉霍7.6级地震后,已有40多年没有发生7级以上地震,鲜有6.5级以上强震,在此次地震序列中表现出一种大震缺失的状态。因此,巴颜喀拉块体东南缘主要断裂带的地震危险性值得关注和研究。同震库仑破裂应力及震后粘弹性松弛在解释余震序列分布,分析强震序列相互作用关系及未来地震危险性分析等方面发挥了重要的作用。比如汶川地震之后,一些研究者利用弹性位错模型计算了汶川地震的同震库仑破裂应力变化,并对龙门山区域未来的地震危险性进行了估计。一些研究者分别计算了汶川地震和芦山地震后周边断层的同震库仑应力变化,以及汶川地震和芦山地震引起的同震和震后粘弹性松弛应力场变化,得到了汶川地震可能有效地促进了芦山地震的发生的结论。在应力演化过程中,震间构造应力加载对强震危险性的作用也非常重要,目前比较流行的计算这部分应力积累的做法是采用"负位错"的理论进行计算,这一方法虽然简洁,但假设了断层是完全闭锁的,可能会造成一定误差。本文拟以巴颜喀拉块体区域构造活动作为三维有限元模型的边界约束,综合研究区域介质速度、密度结构最新反演结果计算介质的弹性材料参数,从初始应力场出发模拟在构造应力和重力作用下的现今背景应力场,计算强震序列引起的区域内断裂带震前、震间、震后的应力演化,获得更为准确可靠的应力场结果。利用速率-状态关系计算库仑应力变化对研究区域背景地震发生率的影响,利用泊松概率过程计算研究区域在汶川地震和芦山地震震后未来10年6级以上强震的发震概率及研究区域主要断裂带的6级以上强震发震概率。除以上述准动力学的方法研究巴颜喀拉块体东南缘应力演化及地震危险性,本文还以鲁甸地震为例,研究鲁甸地震从成核区开始破裂,然后扩展直至终止的自发破裂动力学过程,探索影响地震自发破裂位错的主要因素。研究内容:(1)建立巴颜喀拉块体东南缘三维粘弹性有限元模型。(2)模拟巴颜喀拉块体东南缘主要断裂带上1900年以来7级以上强震序列的应力演化过程,并对强震序列对芦山地震的触发作用进行了研究,对比库仑应力场及等效应力场演化结果与背景地震活动的关系。(3)应用速率-状态关系研究库仑破裂应力对汶川地震、芦山地震震后研究区域的地震发生率的影响,应用泊松过程给出研究区域主要断层未来的强震发生概率。(4)应用可以描述断层摩擦机制的FAULTS有限元软件对芦山地震断层错动方式,初始破裂点进行了模拟,研究了汶川地震后,断层强度的变化对芦山地震的影响。(5)应用曲线有限差分方法建立鲁甸地震非平面断层自发破裂模型,模拟鲁甸地震在包谷垴—小河断裂上的自发破裂过程,分析影响鲁甸地震自发破裂过程的主要因素。本文的创新点:(1)以三维速度、密度反演结果建立能反映更真实地下介质的三维粘弹性有限元模型;(2)从初始应力场出发,进行数十万年的应力场演化模拟直至获得与现今构造应力场更为吻合的背景应力场用于地震序列模拟;(3)采用泊松概率模型计算地震发生概率时直接基于三维粘弹性有限元模型得到的库仑应力结果,充分考虑了同震、震间和震后效应的影响。同时,在考虑背景地震发生率时采用了更全面完整的地震目录,这可能使得对未来地震危险性的估计更准;(4)对鲁甸地震自发破裂过程进行了动力学模拟,将震间应力积累过程与同震地震破裂过程的模拟结合起来。     

汶川8.0级地震的余震触发作用和对断层的应力加卸载作用

采用可变滑动震源模型,计算和研究了汶川M_S 8.0地震产生的应力变化及对后续6次强余震的应力触发作用,并定量分析了汶川地震对附近活动断裂带的库仑应力加卸载作用。结果显示:汶川地震产生的库仑破裂应力变化分布图案复杂性强,主要由发震断裂复杂的空间展布形态以及震源断层活动具有逆冲兼右旋走滑两种性质决定。汶川地震产生的库仑破裂应力变化对多数强余震存在一定的触发作用。汶川地震的发生对附近的断裂带有不同程度的影响,其中青川断裂中段和北段、岷江断裂南段和灌县-江油断裂南段,主要是以强烈的库仑破裂应力加载作用占主导地位,有利于强余震序列的孕育、加速以至于发生;东昆仑断裂带和鲜水河断裂带受到了轻微的应力加载作用;龙泉山断裂带、华莹山断裂带和西秦岭北缘断裂带主要受到了轻微的卸载作用。