汶川8.0级地震前GPS与跨断层资料反映的运动与变形演化特征

为了研究与总结2008年5月12日汶川8.0级地震前GPS与跨断层资料反映的龙门山断裂带及其周边地区的运动、构造变形、应变积累演化过程,以及汶川地震临震阶段可能的物理机制,本文综合1999~2007期GPS速度场、1999~2008年大尺度GPS基线时间序列、1985~2008年跨断层短水准等资料进行了相关分析与讨论。结果表明:(1)GPS速度剖面结果显示,宽达500km的川西高原在震前有明显的连续变形,而四川盆地一侧和跨龙门山断裂带基本没有变形趋势,表明震前川西高原在持续不断地为已经处于闭锁状态的龙门山断裂带提供能量积累。(2)GPS应变率结果显示,震前龙门山断裂带中北段的NW侧EW向挤压变形明显,变形幅度从远离断裂带较大到靠近断裂带逐渐减小,而断裂带变形微弱;龙门山断裂带西南段周边形成了显著的EW向挤压应变集中区,应变积累速率明显大于中北段。(3)断层闭锁程度反演结果显示,除了汶川地震的震源位置闭锁相对较弱,且西南段有大概20km宽度断层在12~22.5km深度为蠕滑状态以外,震前整条龙门山断裂基本处于强闭锁状态。(4)大尺度GPS基线结果显示,跨南北地震带区域的NE向基线从2005年开始普遍出现压缩转折,反映NE向地壳缩短的相对运动增强。(5)跨断层短水准场地结果显示,震前年均垂直变化速率和形变累积率很低,表明断层近场垂向活动很弱、闭锁较强。通过以上分析认为,在相对小尺度的地壳变形中,震前龙门山断裂带深浅部均处于强闭锁状态,断裂带水平与垂直变形都很微弱,这可能经历了一个缓慢的过程,而且越是临近地震的发生,微弱变形的范围可能越大;在相对大尺度的地壳变形中,震前龙门山断裂带西侧的巴颜喀拉块体东部地区经历了地壳缓慢且持续的缩短挤压变形,为龙门山断裂带应变积累持续提供了动力支持。     

近年来龙门山断裂GPS剖面变形与应变积累分析

利用1999~2012年4期GPS速度场数据,对龙门山断裂天全—映秀段及安县—广元段的GPS剖面变形与应变积累作了探讨。在剖面应变积累对地震触发成因及程度方面,将芦山MS7.0与汶川MS8.0地震进行了对比分析。结果表明:(1)芦山地震之前,该地区压性速率远高于1999~2007年,它是汶川大震后的应力调整过程;(2)汶川MS8.0、芦山MS7.0两次强震间隔如此之短,不排除是汶川地震之后龙门山断裂西南段应力场高值压应力的调整过程对此次芦山地震起到触发作用的可能性。     

龙门山断裂带西南段闭锁与变形特征

为研究龙门山断裂带西南段断层闭锁程度与变形状态,并分析该区域地震危险性,利用1999—2007、2009—2011、2011—2013和2013—2015年共4期GPS速度场,采用DEFNODE负位错反演程序估算了该断裂的闭锁程度和滑动亏损空间动态分布,并结合1990—2017年跨断层水准资料分析了断层的三维运动变形特征。结果表明:(1)龙门山断裂带西南段在汶川地震前后一直处于较强的闭锁状态,且汶川地震使西南段应变积累速度加快,加速了芦山地震的孕育进程;芦山地震后西南段闭锁程度并没有明显减弱,芦山地震对西南段的应变能释放是局部的和有限的。(2)龙门山断裂带4期垂直断层滑动亏损速率均为挤压亏损速率,汶川地震后西南段亏损速率明显增大,而后2期西南段亏损速率逐渐减小,目前依然明显高于汶川地震前量值。(3)汶川地震前跨断层水准结果显示龙门山断裂带西南段处于完全闭锁状态,汶川地震后多个场地的年均垂直变化速率明显增大,并随时间呈逐渐衰减状态,目前已经恢复至汶川地震前正常变化速率,因此西南段仍然处于闭锁状态。综合GPS反演结果和跨断层水准结果分析认为,目前龙门山断裂带西南段在大部分段落处于强闭锁状态下依然有发生大地震的可能性。     

基于青衣江变形河流阶地研究龙门山断裂带南段的构造活动性

龙门山断裂带南段第四纪沉积差,断层出露不明显,晚第四纪构造活动性资料零星。为了提高对龙门山断裂带南段构造活动性的认识,探索芦山地震的发震构造,文中在分析龙门山断裂带南段的地貌以及构造演化的基础上,对跨盐井-五龙断裂、大川-双石断裂和芦山盆地的青衣江不同段的6级河流阶地进行了差分GPS连续测量和细致研究,结合对高分辨率航拍影像的地质解译,得到了龙门山断裂带南段青衣江各段的河流阶地横剖面,通过不同河段河流阶地的对比分析,建立了龙门山断裂带南段青衣江河流阶地纵剖面。通过对河流阶地的变形分析,发现龙门山断裂带南段晚第四纪以来,盐井-五龙断裂的平均垂向断错速率为0.6~1.2mm/a,大川-双石断裂没有明显的垂向活动,芦山地震的发震断层控制的山前褶皱最新活动。结合龙门山断裂带南段的地壳深部结构资料和芦山地震的精定位余震资料等,认为芦山地震的发震构造不是大川-双石断裂,而是龙门山断裂带南段的山前盲逆断层和反冲断层。     

利用水准测量数据分析芦山M_S7.0地震前龙门山断裂带南段应变积累

利用2010年和2013年两期一等水准测量数据分析了龙门山断裂带南段芦山MS7.0地震前的应变积累。结果表明:(1)汶川地震的发生明显加速了该区域应变积累的过程,龙门山断裂带南段主要断层区域垂直形变速率为6~9 mm/a,表明汶川地震的发生加速了此次芦山地震的孕育过程;靠近鲜水河断裂带区域垂直形变速率为1~5 mm/a,低于汶川地震前的隆升速度。(2)芦山地震并未释放该区域长期积累的应变能,龙门山断裂带南段仍然具有发生破坏性地震的可能。     

龙门山断裂带中南段的一个破裂空段——芦山地震的震后效应

2013年4月20日,在青藏高原东缘的龙门山断裂带的南段发生了芦山7.0级地震.通过地壳剪切波分裂的分析,揭示了龙门山断裂带域的地壳主压应力方向及其与断裂之间的关联.芦山地震余震的定位结果显示,芦山地震的破裂与汶川地震的破裂没有贯通,在芦山与汶川之间形成了一个"破裂空段".本文建议,应关注芦山地震与汶川地震之间"破裂空段"的应力变化.采用不依赖于天然地震记录的井下人工源观测技术是一种有意义的科学尝试.     

2008年汶川8.0级地震前横跨龙门山断裂带的震间形变

利用区域GPS和水准测量资料,结合地震构造背景的分析,本文研究2008年汶川8.0级地震前横跨龙门山断裂带地区的震间地壳形变,探讨引起这种形变的活动构造与动力学模式,并由此认识汶川地震的孕育与成因机制.主要结果表明：1997～2007年期间,自龙门山断裂带中段朝北西约230 km的地带内存在垂直于断裂的水平缩短变形、以及平行于断裂的水平右旋剪切变形,缩短率为1.3×10^-8/a （即：0.013 mm/km/a）,角变形速率为2.6×10^-8/a;同一地带在1975～1997年期间还表现出垂直上隆变形,上隆速率在龙门山前山断裂与中央断裂之间仅0.6 mm/a,而至龙门山后山断裂及其以西达2～3 mm/a.这些反映了在汶川地震之前至少10～30余年,龙门山断裂带中段的前山与中央断裂业已闭锁、并伴有应变积累.造成这种形变的主要原因是：以壳内的低速层为“解耦”带,巴颜喀拉地块上地壳朝南东的水平运动在四川盆地西缘受到华南地块的阻挡、转换成龙门山断裂带中段的逆冲运动;由于该断裂段的震间闭锁,致使西侧的巴颜喀拉地块的上地壳发生横向缩短以及平行断裂的右旋剪切变形.然而,龙门山断裂带北段在1997～2007年期间除了有大约0.9 mm/a的右旋剪切变形外,横向的缩短变形极微弱,这可能与该断裂段西侧的岷江、虎牙、龙日坝等断裂带吸收了巴颜喀拉地块朝东水平运动的大部分有关.另外,汶川地震前,横跨龙门山断裂带中段与北段的地壳形变特征的差异,与汶川地震时能量释放的空间分布吻合.     

青藏高原东缘川西地区的现今构造变形、应变分配与深部动力过程

青藏高原东缘的川西地区是中国大陆南北地震带的中段,构造变形复杂,断裂活动强烈,控制着一系列历史强震的发生,2008年5月12日汶川8．0级地震就发生在南北地震带中段的龙门山断裂带．川西地区构造变形图像、运动特性和深部驱动机制的研究,不仅对于理解青藏高原东边界的动力过程具有重要意义,还有助于认识南北带未来强震危险性．研究表明,北西西走向的鲜水河断裂是中国大陆内部活动最强烈的断裂带之一,左旋走滑速率达9-11mm/a,其左旋走滑运动在东部的石绵一带被分解到了安宁河、大凉山和则木河等断裂带之上,仍以左旋走滑为主;再向南,左旋剪切走滑运动沿小江断裂带发生,并以8-9mm／a的走滑速率进入云南,跨过著名的红河断裂,延伸入缅甸境内．龙门山断裂则以地壳缩短和右旋走滑为特征,但其总体滑动速率只有约3mm/a．GPS观测结果也证实了这种运动方式和应变分解图像．青藏高原东边界的构造变形可能是一种近似连续的旋转弧形构造,以左旋剪切为主要特征,只是在向北东方向突出的弧顶部位出现少量的挤压逆冲构造．造成这种弧形左旋走滑运动的深部驱动机制可能是中下地壳的流动．因为川西高原的中下地壳流变强度比正常地壳软弱,并且其与临近四川盆地和华南地块的地壳厚度差达20—30km、地貌高度差达3000—4000m,从而形成能够驱动中下地壳软弱层发生流动的横向压力差,从底部拖曳着被断裂切割的上部脆性地壳,发生以左旋走滑为主的变形和运动,并导致应变在不同的断裂上积累和释放,形成强烈地震．这种由北西西走向左旋剪切转换为近南北向左旋剪切兼挤压的变形模式在中国大陆具有普遍意义．     

2008年汶川M_S8.0地震后GPS资料反映的龙门山断裂带南段地壳水平变形动态特征分析

本文对龙门山断裂带南段GPS位移时间序列进行了欧拉基准转换、最小二乘配置插值和滤波分析等数据后处理,并利用块体刚性-弹塑性运动-应变模型计算了龙门山南段各区域应变参数时间序列。结合各站点位移时间序列、速度场及各区域应变参数序列,分析发现2008年汶川MS8.0地震对龙门山南段的影响由西北至东南由大变小,引起最大变形速度达背景速度6倍,增强了隐伏断层对两侧地表影响的差异,使南段的地震危险性由北至南增加。根据龙门山南段的历史地震及构造背景,对各区域地震危险性进行了比较,发现东南侧危险性最大,历史应变的积累及汶川地震的加速影响导致了芦山地震的提前发生和破裂断层的左旋错动。汶川地震对较远的鲜水河断裂南段也有持续两年的应力微释放影响,但其远不及川滇块体的应力压缩作用。     

芦山地震断层的滑动分布与汶川地震断层的关系

本文利用2013年芦山M_S7.0级地震同震GPS数据反演了芦山断层几何与断层滑动分布,结果表明:芦山地震发震断层具有南陡北缓、上陡下缓的特征,低倾角的区域位于发震断层北段且靠近映秀断层的一侧;滑动分布模型的最大滑动量为0.82m,其深度为13.67km与小震发生集中平均深度12.5km接近.我们选取1998—2014年龙门山断裂带区域地壳形变观测数据,拟合获得了龙门山断裂带走向方向上的速度分量,发现在汶川M_S8.0地震与芦山M_S7.0地震之间宽度约30km破裂空区,龙门山断裂带西南段与东北段的形变分量以破裂空区为界方向相反.断裂带东北段(汶川地震主要发震断层)的形变分量方向与断层右旋走滑运动方向一致,而在断裂带西南段(芦山地震发震断层)的形变分量方向与断层左旋走滑运动方向一致.芦山地震走滑方向与汶川地震走滑方向相反是因为该断裂带构造运动在特有几何构造下受青藏高原东南向挤压,遇龙门山中段岩石圈楔状构造的阻挡,在汶川M_S8.0地震与芦山M_S7.0地震间的地震空区,形成了构造运动向其两侧分流的结果.     

汶川地震前龙门山断裂带闭锁程度和滑动亏损分布研究

利用1999—2007期GPS水平速度场数据,采用Defnode负位错反演程序估算了龙门山断裂在汶川地震前的闭锁程度和滑动亏损分布,结合龙门山断裂带附近地表水平应变率场结果,综合分析了震前地壳变形特征.反演结果表明,震前龙门山断裂中北段处于完全闭锁状态,闭锁深度达到21 km(闭锁比例0.99)左右,垂直断层方向的挤压滑动亏损速率约为2.2 mm/a,平行断层方向的右旋滑动亏损速率约为4.6 mm/a.龙门山断裂南段只有地表以下12 km闭锁程度较高(闭锁比例0.99),垂直断层方向滑动亏损速率约为1.4 mm/a,平行断层方向滑动亏损速率约为4.6 mm/a;在12~16 km处闭锁比例约为0.83,垂直断层方向滑动亏损速率约为1.2 mm/a,平行断层方向滑动亏损速率约为3.8 mm/a;在16~21 km处闭锁比例约为0.75,垂直断层方向滑动亏损速率约为1.1 mm/a,平行断层方向滑动亏损速率约为3.5 mm/a.在21~24 km处整条断裂均逐步转变为蠕滑.上述反演结果与区域应变计算获得的龙门山断裂带中北段整体应变积累速率较低、南段应变积累速率较高相一致,均表明中北段闭锁程度高、南段闭锁程度稍低,该特征可以较好地解释汶川地震时从震中向北东向单向破裂现象.     

汶川8.0级地震前区域地壳运动与变形动态过程

利用GPS等观测资料研究了2008年汶川8.0级地震前的区域地壳运动与构造变形、应变积累,以及大区域地壳运动微动态变化过程,结合同震位移场分布等分析、讨论了汶川地震前近10年区域地壳变形的表象所反映的大震孕震最后阶段的物理过程.结果表明,发生汶川8.0级地震的龙门山断裂带1999～2004年相对稳定的华南地块无明显的相对运动,2004～2007年有一定程度的相对运动显示,以右旋活动为主,年速率为1.6mm/a.龙门山断裂带西北侧的巴颜喀拉地块东部为右旋剪切为主兼有推挤的大尺度缓慢变形状态,右旋扭动变形率为2.1±0.2(10^-8/a),地壳总体缩短变形率为-0.7±0.1(10^-8/a).由GPS基准站资料计算的基线时间序列显示,2005年以来大尺度北东向地壳缩短的相对运动明显增强,青藏块体相对华南地块的北东向运动明显增强.对汶川大地震前应变积累的特殊性等问题进行了初步讨论,分析认为,汶川大地震的发生是四川盆地西缘的龙门山断裂带受到其西侧巴颜喀拉地块推挤导致大尺度、长时期、缓慢的地壳应变积累的结果.在发震前不仅龙门山断裂带为显著闭锁状态,且孕震区域可能存在地壳弹性变形趋于极限后难以发生变形的相持阶段,大尺度地壳运动与动力环境的动态变化对汶川大地震的发生有促进影响.     

汶川地震前龙门山及其周缘断裂形变运动与应力累积的数值模拟

为揭示汶川地震前龙门山及其周缘断裂的形变与应力累积状态,文中构建了包含龙门山、龙日坝、岷江和虎牙4条断裂的三维黏弹性有限元模型,以1999—2004年GPS结果为约束,模拟了龙门山断裂带及其周缘区域的形变运动。得到以下结论:1)平行于龙门山断裂带的速度分量主要被龙日坝断裂吸收,垂直于龙门山断裂带的速度分量主要被其自身吸收;岷江和虎牙断裂对龙门山断裂带北段起到一定的屏障作用,导致其北段压缩量明显低于南段。2)沿龙门山断裂带由SW向NE方向延伸,主压应力与断层走向的夹角由接近垂直逐步转至约45°;断层南段挤压、剪切应力累积速率高,且压应力大于剪应力,北段应力累积速率低,压应力与剪应力接近。这与龙门山断裂带SW段中小地震频发、地震活动强烈,NE段偶有小震、地震活动微弱相吻合;也与汶川M_S8.0地震逆冲兼具右旋走滑、芦山M_S7.0地震逆冲破裂的方式相一致。3)假设发生震级、类型相同的地震所需应力积累量相同,那么研究区内岷江断裂、龙门山断裂南段和虎牙断裂破裂以逆冲运动为主,3条断裂的地震复发周期依次变长;龙日坝断裂北段和龙门山断裂北段以逆冲兼具右旋走滑为主,前者地震复发周期短于后者;龙日坝断裂南段则以纯右旋走滑为主,地震复发周期有可能最短。     

芦山地震前龙门山断裂带闭锁程度 与变形动态特征研究

利用1999-2007和2009-2011年中国大陆GPS水平速度场数据, 采用DEFNODE(反演计算弹性岩石圈块体旋转、 应变和块体边界断层闭锁或同震滑动的Fortran程序)负位错反演程序估算了芦山地震前龙门山断裂带的三维闭锁程度, 并结合剖面结果分析了断层深浅部变形特征. GPS反演结果表明, 1999-2007年, 龙门山断裂中北段(闭锁比例为0.99)处于强闭锁(本文将闭锁比例大于0.97的称为强闭锁)状态; 龙门山断裂南段地表以下深度16 km内为强闭锁, 深度16-21 km处闭锁比例降低为0.62, 深度21-24 km处整条断裂逐渐转变为蠕滑状态. 2009-2011年, 即汶川地震后, 龙门山断裂中北段处于震后蠕滑状态; 龙门山断裂南段深度16-21 km处闭锁比例降低为0.45, 其它位置闭锁程度保持不变. GPS剖面结果显示, 2009-2011年, 即汶川地震后, 龙门山断裂中北段为逆冲兼右旋走滑运动; 而南段断层不能自由滑动、 变形宽度较大. 综合分析认为, 汶川地震时, 龙门山断裂南段并没有发生破裂, 一直处于较强的闭锁状态, 汶川地震的发生又加速了芦山地震的孕育进程; 由于龙门山断裂带南段的闭锁深度较中北段浅, 因此芦山地震较汶川地震强度低、 震级小、 破裂范围窄.      

汶川8.0级地震前区域地壳运动与应变场动态特征

利用中国地壳运动观测网络GPS资料, 通过获取水平相对运动、 水平应变场分布变化等, 研究了汶川8.0级地震前的区域水平运动与应变率场变化, 以及大尺度地壳运动动态特征等。 研究结果表明, 发生汶川8.0级地震的龙门山断裂带, 由于受到其西侧巴颜喀拉地块向东运动的构造动力作用, 处于缓慢的应变积累状态。 在汶川地震前龙门山断裂带相对华南地块的差异运动小于GPS观测误差。 川滇地区应变场图像显示2004-2007年面应变率负值最高区出现在汶川8.0级地震震中区及其附近, 可能反应了局部挤压增强。 GPS基准站资料反映的大尺度北东向地壳缩短的相对运动增强, 也形成了促进龙门山断裂带发生大破裂的区域构造动力增强的背景条件。     

龙门山断裂带中段及邻区构造变形特征及应力场演化的数值模拟研究

根据龙门山断裂带中段及邻区的主要构造特征,同时考虑地壳速度结构和构造运动的深浅部差异,建立该地区的二维接触有限元模型,数值模拟了长期构造过程中断裂带位移和应力的时间演化过程。结果显示:①汶川—茂县断裂、映秀—北川断裂、灌县—安县断裂的逆冲速率与实际观测结果基本一致,映秀—北川断裂逆冲性质最强,巴颜喀拉块体和断裂带区域呈水平缩短与增厚抬升趋势,四川盆地保持相对稳定;②断裂带附近岩体大小主应力比随深度增加而减小,总体上与实际测试结果相符,且主应力轴方向和倾角与利用震源机制解反演得到的构造应力场特征亦较为接近;③龙门山断裂带中段区域构造应力场呈低态稳定—应力累积—临界稳定的阶段式演化过程,与该区所处的应力累积—地震发生—应力再累积的地震活动特征相呼应。此外,仅根据地表调查和震前GPS观测所揭示的断裂活动速率评价活动断裂区的地震危险程度具有一定局限性,需结合其深部动力学过程和区域应力状态进行综合分析。     

2008年汶川MS 8．0地震后GPS资料反映的龙门山断裂带南段地壳水平变形动态特征分析

本文对龙门山断裂带南段GPS位移时间序列进行了欧拉基准转换、最小二乘配置插值和滤波分析等数据后处理,并利用块体刚性一弹塑性运动一应变模型计算了龙门山南段各区域应变参数时间序列。结合各站点位移时间序列、速度场及各区域应变参数序列,分析发现2008年汶川M 8．0地震对龙门山南段的影响由西北至东南由大变小,引起最大变形速度达背景速度6倍,增强了隐伏断层对两侧地表影响的差异,使南段的地震危险性由北至南增加。根据龙门山南段的历史地震及构造背景,对各区域地震危险性进行了比较,发现东南侧危险性最大,历史应变的积累及汶川地震的加速影响导致了芦山地震的提前发生和破裂断层的左旋错动。汶川地震对较远的鲜水河断裂南段也有持续两年的应力微释放影响,但其远不及川滇块体的应力压缩作用。     

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

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

地壳流变结构控制作用下的龙门山断裂带地震发生机理

青藏高原东缘低地形变速率的龙门山断裂带上相继发生了2008汶川M_w7.9级地震和2013芦山M_w6.6级地震.地震勘探与震源定位结果揭示了龙门山区域地震空间分布特征：纵向上,龙门山断裂带这两次地震主震均发生在龙门山断裂带上地壳的底部（14~19 km）,绝大部分余震均发生在上地壳范围（5~25 km）,而在其中、下地壳深度范围内鲜见余震发生;横向上,地震（M_w>3）在龙门山断裂带青藏高原一侧密集分布且曾有大震发生,而四川盆地地震稀少（M_w>3）.为探讨龙门山断裂带地震发生机理,并解释以上龙门山区域地震空间分布特征,本文建立了龙门山断裂带西南段跨芦山地震震中区域的四种不同流变结构的龙门山断裂带三维岩石圈模型,以地表GPS观测资料为约束边界条件,数值模拟龙门山断裂带岩石圈在数千年以上长期匀速构造挤压作用下的应力积累特征,探讨了地壳分层流变性质对地壳应力积累的影响,分析了该区域地震空间分布与构造应力积累速率的关系.计算结果表明：该区域在数千年的应力积累过程中,脆性上地壳中应力表现近于恒定值的线性增长趋势,龙门山断裂带上地壳底部出现应力集中积累现象,这一应力集中现象可以解释龙门山断裂带汶川地震与芦山地震主震的发生,及其大部分余震在脆性上地壳中的触发;青藏高原一侧上地壳应力积累速率远远高于四川盆地的应力积累速率,这一应力积累分布现象可以解释龙门山区域青藏高原一侧地震密集而四川盆地地震稀少的地震空间分布特征;通过比较不同流变结构模型中的应力积累状态,认为导致这一应力积累空间分布状态的重要控制因素在于青藏高原中、下地壳较低的黏滞系数与四川盆地中、下地壳较高的黏滞系数的差异.在柔性的中、下地壳内,应力增长近于指数形式,稳定状态之后其应力增长速率近于零,构造应力积累难以达到岩石破裂强度,因而鲜见地震发生.地壳各层位的应力增长率差异与地震成层分布的现象共同揭示了龙门山区域岩石圈分层流变结构：脆性上地壳、韧性中、下地壳（青藏高原一侧较弱,四川盆地一侧较强）、韧性岩石圈上地幔.     

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

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