MS8.0汶川地震对MS7.0芦山地震成核失稳的影响

依据速率-状态依赖性摩擦本构关系,并结合中下地壳和上地幔的黏弹性松弛效应,以震后库仑破裂应力变化和同震动态库仑破裂应力变化的计算为基础,模拟了两种应力变化对芦山地震断层的失稳发震时间的影响,研究了2008年M_S8.0汶川地震与2013年M_S7.0芦山地震之间的触发关系.计算得到,汶川地震在芦山地震震源断层面上产生的动态应力变化的峰值为0.127 MPa;此外,经过近5年中下地壳和上地幔的黏弹性松弛效应,芦山地震震源断层面上受到的震后应力变化值为0.025 MPa.结果表明,芦山地震的震源断层在应力积累逐渐接近临界状态的某一特定时期内,受到了汶川地震产生的动态应力变化、静态应力变化以及黏弹性松弛效应造成应力变化的共同触发作用,且动态应力的延迟触发作用可能更为显著.最后对芦山地震之后研究区域的应力变化场进行了初步探讨.     

从应力演化角度讨论龙门山区域7级强震序列对芦山地震的影响

计算包括同震静态库仑破裂应力变化以及震后黏滞松弛引起的应力变化等,可以更好地解释余震分布、地震序列等地震观测结果。在芦山地震之前,从1900年以来龙门山区域发生了4次7级以上地震,分析这几次地震的同震应力变化以及震后黏弹性松弛对芦山地震的产生的影响,芦山地震以后区域断裂带上的应力伴随强震如何演化,芦山地震与汶川地震的破裂空段呈现何种应力状态,探讨这些问题可能会为了解芦山地震震源处震前的应力状态及该区域未来地震风险评估提供一定的依据。文中采用有限元数值模拟方法,根据地质构造、速度、密度结构深部反演结果以及GPS及应力观测资料等,建立龙门山地区三维黏弹性有限元模型进行研究。模拟结果显示:龙门山断裂带南段及鲜水河断裂带的库仑应力年变化速率在研究区域中相对更高,这与研究区域的地震活动性一致。芦山地震的前4次地震,除叠溪地震外,康定、松潘、汶川等3次地震在芦山地震震中位置产生的同震库仑破裂应力变化大于0,表明这3次地震可能促进了芦山地震的发生,汶川地震的同震库仑破裂应力超过了0.01MPa,同震触发效应十分显著。震间的黏弹性松弛对芦山震源处起加载作用,从1900年以来这种持续的加载作用也超过0.01MPa,因此在模拟应力演化的时候,介质的黏弹性松弛效应不能被忽略。从库仑破裂应力的角度计算龙门山区域断层的应力演化,可以发现龙门山断裂带上汶川地震和芦山地震破裂的空段,在芦山地震之后仍然属于相对应力水平较高的区域。     

汶川地震后中下地壳及上地幔的粘弹性效应引起的应力变化与芦山地震的发生机制

采用多层粘弹性模型计算了2008年5月12日汶川Mw7.9地震对周围地区尤其是龙门山断裂带南段的影响,结合地震活动性分析,探讨了2013年4月20日芦山Mw6.6地震的发生机制,并对沿龙门山断裂带芦山地震与汶川地震之间的地震空段进行了分析.计算结果表明,由于下部地壳及上地幔的粘弹性效应,芦山地震震源处对应的汶川地震同震库仑应力变化(△CFS)为-0.1bar以下量级,其后随时间逐渐增加,在芦山地震前增加到了0.4bar(有效摩擦系数为0.4)或0.6bar(有效摩擦系数为0.2)以上.这表明,芦山地震的发生与汶川地震的非弹性触发密切相关.芦山地震与汶川地震之间的地震破裂空段存在1970年大邑M6.2地震的库仑应力阴影区(下降区),汶川地震及芦山地震未能使阴影区发生根本变化,因此,如没有其他不可知因素,该空段短期内整体发生破裂,引发大地震的可能性不大.     

汶川8.0级地震和芦山7.0级地震对周边断层的影响

采用分层粘弹性介质模型,研究龙门山断裂带发生的两次大地震(2008年5月12日汶川8.0级和2013年4月20日芦山7.0级地震)的同震位错和震后粘滞松弛效应引起的周边断裂带上的库仑应力变化。结果表明,汶川8.0级地震同震和震后效应均引起了芦山7.0级地震震中位置的库仑应力增加,有利于芦山地震的发生;汶川和芦山地震的同震和震后效应引起的周边各断层上的累积库仑应力变化分布显示,这两次大地震引起了其间的破裂空段的库仑应力显著增加;此外,累积库仑应力显著升高的断裂还有:鲜水河断裂带中段、灌县一安县断裂北段、文县断裂的一.部分、平武一青川断裂北段、略阳一勉县断裂。累积库仑应力显著降低的断裂有:马尔康断裂、岷江断裂南端、虎牙断裂、龙泉山断裂北端、北川一映秀断裂北段。     

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

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

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

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

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

从库仑破裂应力和余震分布角度探讨汶川地震和芦山地震的关系

本文以龙门山及周边地区为研究对象,考虑区域地质构造差异、主要活动断裂带、地表附加重力影响,建立能反映地表起伏和岩石圈分层结构的龙门山地区三维粘弹性有限元模型。以地壳水平运动速率观测值为约束条件重建研究区现今构造背景应力场,在此基础上分别模拟了汶川地震和芦山地震的发生机理。通过分析同震库仑破裂应力变化与余震空间分布的关系,探讨了2次地震主震对余震的触发作用以及汶川地震对芦山地震的影响。研究表明,汶川地震和芦山地震的余震大部分由其主震触发,汶川地震对芦山地震的余震有约6.78%的触发作用。汶川地震的同震库仑破裂应力在芦山地震主震位置的增加值约为0.016MPa,如果龙门山断裂带南段库仑破裂应力年累积速率按照0.4×10-3-0.6×10-3MPa·a-1计算,汶川地震使芦山地震提前了约27-40年。计算还表明汶川地震和芦山地震的发生使鲜水河断裂带南段和虎牙断裂的库仑破裂应力增加,这些断裂带在未来发生地震的可能性增加。     

强震后地表变形的动力学机制研究—以1999年台湾集集地震为例

强震后地表变形的动力学机制是地球动力学研究的重要方面,现在普遍认为震后变形主要由断层的震后余滑或由介质的黏弹性松弛所至。1999年台湾集集地震GPS观测系统记录到了空前的资料,为研究震后变形的动力学机制提供了难得的机会。本研究认为集集地震后地表变形由震后断层余滑、下地壳／上地幔的黏弹性松弛、震源区介质的破裂、孔隙弹性回跳、地下流体的运移、介质孔隙度及孔隙压的变化等多种因素共同影响决定。为抓住重点,研究中将介质的破裂、地下流体的运移和孔隙弹性回跳等因素等效为震源区介质的物性变化。文中运用黏弹性有限单元模型（麦克斯威尔体）、利用GPS观测的时间序列资料对震后余滑、地壳／地幔黏度以及等效的震源区介质物性变化进行了反演。反演模型给出了震后余滑的分布及变化特征,反演结果初步显示台湾地区的下地壳／上地幔的黏度分别为2.7×10^18,4.2×10^20Pa·s。此外,反演结果还给出每种影响因素对地表变形的贡献大小,在集集地震后的450d时间里,断层的震后余滑引起的地表变形占总变形的44.6％,下地壳／上地幔的黏性松弛占34.7％,等效的震源区介质的物性变化占20.7％。     

芦山Ms7.0地震产生的静态库仑应力影响研究

2013年4月20日在四川芦山发生了Ms7.0地震,震源深度17 km,发震断层位于龙门山断裂带南段的前山断裂附近,为逆冲断层机制,芦山地震震中位置位于2008年汶川Ms8.0地震西南方向85 km处,两个地震的余震分布区中间存在一个约20~30 km的破裂空段。首先计算汶川地震的静态库伦应力,明确汶川地震对芦山地震的触发影响,芦山地震的发震位置处于汶川地震后的应力加载区域;进一步计算芦山地震对周围的静态库仑应力,结果显示,芦山地震对二者之间破裂空段的影响为应力加载;最后计算汶川地震和芦山地震周围的静态库仑应力的综合影响,二者对川滇交界东部的应力加载作用是明确的,对中间存在的20~30 km的空段也存在一定的应力加载影响。     

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%的余震位于库仑应力增强区,地震应力触发理论较好地解释了余震的分布.选取囊谦地震震源机制解的两个节面作为库仑应力计算中的接收断层参数,并且考虑不同黏滞系数下的玉树地震同震及震后黏弹松弛效应,模型计算结果均表明囊谦地震位于玉树地震所导致应力影区,仅依靠地震的静态、震后黏弹松弛应力触发理论,无法解释囊谦地震的发生,说明该次地震可能是一次独立的事件.     

汶川地震对龙泉山断裂地震发生概率的影响研究

通过构建汶川Ms8.0级地震震源断层模型和龙泉山断裂模型,计算龙泉山断裂在汶川地震后的库仑应力变化,结果显示:汶川地震使龙泉山断裂北部(A区)库仑应力增大0.46～0.58 bars,中部(B区)库仑应力增加0.12～0.17 bars,南部(C区)基本上无库仑应力变化.汶川地震后至2010年1月龙泉山断裂沿线的Ms≥...     

2015年尼泊尔M_W7.9地震对青藏高原活动断裂同震、震后应力影响

2015年尼泊尔MW7.9地震重烈度区从震中向东延伸,致灾范围包括尼泊尔、印度北部、巴基斯坦、孟加拉和中国藏南地区,其应力调整对邻区和周边活动断裂可能产生重要影响.本文基于地震应力触发理论,采用岩石圈地壳分层黏弹性位错模型,计算了尼泊尔MW7.9地震引起的周边断裂,特别是青藏高原活动断裂的同震和震后库仑应力变化.结果显示,尼泊尔地震同震效应引起大部分震区库仑应力升高,余震主要分布在最大同震滑动等值线外部库仑应力升高区域;少量余震靠近最大滑动量区域,可能该区域积累的地震能量在主震期间没有完全释放.尼泊尔地震同震库仑应力对青藏高原,特别是中尼边境区域活动断裂有一定影响.亚东—谷露地堑南段、北喜马拉雅断裂西段、当惹雍错—定日断裂和甲岗—定结断裂同震库仑应力升高,其中当惹雍错—定日断裂南端,北喜马拉雅断裂西段同震库仑应力变化峰值超过0.01 MPa;帕龙错断裂、班公错断裂、改则—洞措断裂库仑应力降低,其地震发生概率有所降低.震后应力影响方面,未来40年内黏弹性松弛作用导致北喜马拉雅断裂、改则—洞措断裂和喀喇昆仑断裂整体应力卸载;藏南一系列正断层震后应力持续上升,其中帕龙错断裂南段受到震后黏弹性库仑应力影响,由应力阴影区逐渐转化为应力增强区,当惹雍错—定日断裂南段应力进一步加强,震后40年其南端应力变化峰值达到0.1345 MPa,亚东—谷露断裂南段应力亦持续增强.藏南正断层的地震活动性值得进一步关注.     

东昆仑活动断裂带大地震之间的黏弹性应力触发研究

对青藏高原北部东昆仑破裂带大地震之间的应力转移和断层相互作用进行研究. 考虑1937年以来沿此破裂带发生的5个M≥7的地震：1937年M7.5花石峡地震，1963年MS7.1都兰地震，1973年MS7.3玛尼地震，1997年MW7.5玛尼地震和2001年MW7.8可可西里地震，模拟了黏弹性成层介质中地震断层错动产生的应力演化过程，并计算了在后续地震破裂面上产生的库仑破裂应力变化. 结果表明，前面4个地震均造成2001年可可西里地震断层面上库仑破裂应力的增加，并且中地壳和下地壳的黏弹性松弛效应使得库仑破裂应力场随着时间的推移而逐渐加强. 在计算过程中定量估计了可可西里地震发生时前面4个地震同震形变和黏弹性松弛导致可可西里地震破裂面上库仑破裂应力变化之间的比值，发现前3个地震由黏弹性松弛造成的变化远远大于同震形变所造成的变化. 可可西里地震之后应力场的模拟表明东昆仑断层中段的东大滩－西大滩断层段（位于可可西里地震破裂以东及都兰地震以西）的库仑破裂应力显著增加，变化值达0.05～0.1 MPa，预示这一地区地震危险性的增加.     

汶川地震和芦山地震的孕震机理及震前中长期地震危险性研究

运用非连续变形分析法与三维有限元法相结合的方法,以GPS资料作为位移速率和震源机制的约束条件,通过数值模拟研究了青藏高原及其东侧邻区构造地块的运动、变形、相互作用及其与近30年来发生于该区的大地震之间的关系。研究中引入了以应力与摩擦强度的比值定义的断层“失稳危险度”,通过数值模拟计算得到了研究区地壳块体边界断层的失稳危险度分布。结果表明,失稳危险度高的地段与近期该区发生的M_S≥7.0地震所在的位置基本一致,其中龙门山断裂带上包括汶川和芦山大地震的发震断层均为失稳危险度最高值地区。计算得到的应变率强度分布图显示,青藏高原东部边缘整条地带均为应变率强度的陡变带,特别是以龙门山断裂带上的陡变最为明显,其西侧应变率强度为东侧的近4倍,而且,这个带位于宽度相同、走向与龙门山断裂带走向相一致的高应变能密度带中,表明这两次大地震前,作为其发震断层的龙门山断裂带已积累了相当高的应变能,失稳危险度高,处于力学上的不稳定状态。模拟计算得到在上地壳层中,2001年昆仑山口西M_S8.1地震引起汶川、芦山地震发震断层的库仑破裂应力增加约0.016 MPa,相当于龙门山断裂带约两年的应力积累,也就是说,使汶川、芦山地震发震断层的失稳破裂提前了约两年。 此外,关于2008年汶川M_S8.0地震的模拟计算表明,汶川地震的发生也使包括芦山地震发震断层的龙门山断裂带西南段和东昆仑断裂带东南端的库仑破裂应力增大,应变能积累增强,这说明汶川M_S8.0地震的发生对已处于失稳危险度较高状态的2013年芦山地震和2017年九寨沟地震发震断层的提前失稳破裂起到了促进作用。      

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

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

芦山M_s7.0地震和汶川M_s8.0地震对川滇地区的库仑应力影响

2013年4月20日,龙门山断裂带南段发生了芦山M_s7.0地震,震前的静态库仑应力计算表明此处位于2008年汶川M_s8.0地震静态库仑应力显著增强地区,芦山地震的发生对区域应力进行了释放和调整。利用川滇地区的活动断裂数据和地震震源机制解,构建了研究区的0.1°×0.1°的网格化接收断层矩阵。计算了自重多层粘弹性介质中,不同模型参数下,芦山地震和汶川地震对川滇地区的静态库仑应力影响。结果发现,芦山地震与汶川地震之间的破裂空段、东昆仑断裂东段玛曲段、塔藏段以及鲜水河断裂中南段道孚-乾宁段静态库仑应力增加明显,地震危险性较高。     

2017年米林6.9级地震与1950年察隅8.6级地震的关系及两次地震对周边活动断层的影响

2017年11月18日,我国西藏林芝市米林县发生M6.9地震.其东南220 km距离处,1950年发生了察隅8.6级大地震.二者同处喜马拉雅东构造结大拐弯处,周边分布多条活动断裂带.察隅地震作为20世纪我国内陆最大的地震,它的发生对周边断层活动性的影响目前研究尚少.本文基于分层半无限空间黏弹性地球模型,计算了察隅地震对周围活动断裂带和米林地震的影响,同时利用有限断层破裂模型计算了米林地震对周边活动断层产生的同震库仑应力变化,并分析了两次地震的关系,讨论了下地壳及以下介质的黏滞系数和断层有效摩擦系数对计算结果的影响.结果表明,察隅地震影响范围广、强度大且持续时间长,对喜马拉雅东构造结周边活动断层均产生了较大的影响,库仑应力变化达到数MPa量级.在下地壳黏滞系数为1.0×10²⁰ Pa·s、摩擦系数为0.4情况下,按照青藏高原主要活动断裂每年由正常构造运动应力累积1~4 kPa计算,察隅地震的发生使米林地震提前了相当于2.79~11.15年的时间.米林6.9级地震影响范围和强度有限,只对震源附近的活动断层产生一定影响,库仑应力变化最大为数十kPa量级,对远处的断层影响较小.下地壳及以下介质的黏滞系数在震后黏滞松弛过程中影响逐渐明显,断层有效摩擦系数对同震库仑应力计算影响较大.     

印尼苏门答腊8.7级大震对中国陆区的影响

发生在2004年12月26日的印尼苏门答腊8．7级地震，是全球有地震台网以来100多年中记录到的少数几次特大地震之一。初步汇集了在中国陆区观测到的地震活动、地壳形变、地下流体等多学科震时和震后效应的异常变化，包括地震活动性（特别是云南地区）的显著增强；地壳形变观测，尤其是钻孔应变观测记录到地壳应力-应变的震时和震后显著变化，以及地下水温度、化学成分、特别是水位的突出变化。同时，还从库仑破裂应力触发、动态应力触发以及下地壳和上地幔流动变形机制等方面，对这些大震效应作了一些成因机理讨论。     

唐山地震序列应力触发的粘弹性力学模型研究

根据前人对唐山地震破裂分布、地壳波速和粘性结构的研究,考虑局部应力场、孔隙流体压力和断层附近软介质的影响,计算了唐山地震产生的,投影到后续大余震断层面和滑动方向上的库仑破裂应力变化。结果表明,随后发生的滦县地震和宁河地震均发生在唐山地震产生的库仑破裂应力变化为正的区域。为研究唐山地震、滦县地震和宁河地震对后续小震的触发作用,根据前人对该地区构造应力场和地震破裂分布的研究,假定构造应力量值为10 MPa,求得了震源附近各处可能的小震震源机制。将上述3次地震产生的应力变化投影到可能的小震破裂面和滑动方向上,发现唐山地震、滦县地震和宁河地震产生的正库仑破裂应力变化的ldquo;蝴蝶rdquo;形分布与后续小震发生的空间分布具有较好的一致性,95%的余震发生在库仑破裂应力变化增加的区域,说明唐山地震序列中前面的大震对后续小震的发生起到了调制作用。该研究结果对大震后余震的危险性快速评估具有一定意义。如果大地震发生后能够快速确定详细的破裂分布和震源区域详细断层及滑动特性资料,本文方法可用来预测未来大余震的发震位置。