从龙门山地质地貌分段性探讨芦山地震与汶川地震的关系*

在龙门山中、北段发生汶川Ms8.0级特大地震5年之后,2013年4月20日在龙门山南段发生了芦山Ms7.0级地震。基于芦山地震基本特征及其所处的龙门山断裂带的构造特征、地貌水系特征、重力异常分布特征,分析芦山地震发生的区域地质背景,并探讨2013年芦山地震与2008年汶川地震之间的关系。初步获得以下认识：①芦山地震并不是汶川地震的余震,它们是两次独立的地震,芦山地震与汶川地震具有一定的关联性,汶川地震可能促进了芦山地震的发生;②下地壳流向上仰冲可能是芦山地震和汶川地震共同的成因机制,这一成因机制可以很好地解释汶川地震与芦山地震之间的空区;③龙门山南段在山前发育多条断裂和褶皱带,使得这一地区由北西向南东的应力在山前多个断裂和褶皱带得到一定的释放,因此不具备发生类似汶川地震这样特大地震的构造条件。     

汶川大地震·龙门山

5·12汶川大地震的根本原因是板块活动的结果.进入新生代以来,南方的印度板块一直向北方推移,经过强烈挤压,使原来的古地中海东延部分逐渐消失,然后又使喜马拉雅山脉隆起.     

基于水准数据的芦山7.0级地震震间期和同震位移场特征

利用1985~2010年跨龙门山断裂带南段的精密水准观测,并消除2008年汶川地震同震变形的影响,获取了龙门山断裂带南段震间期的地壳垂直形变场.结果表明,双石-大川断裂(龙门山前山断裂)和新开店断裂(大邑断裂南段)之间的芦山地块处于正、负垂直形变梯度带相交的区域,说明该区域在芦山地震前25年内已闭锁.根据2010~2013年穿过震区的水准观测数据,同时消除震前近3年的地壳垂直变形,得到了芦山地震同震垂向位移.同震位移场揭示距离震中最近的达多35水准点最大同震位移达到198.4 mm(相对马雅165基),西北盘(上盘)相对东南盘(下盘)向上运动,显示发震断层以逆冲错动为主.对比芦山地震垂向同震变形和震间期变形,可发现该地震符合弹性回跳学说,芦山地块震间期积累的弹性应变能(位移亏损)在震时得到弹性回跳释放.     

汶川-映秀M_S8.0地震的地球物理场与动力过程

汶川-映秀M_S8.0大地震发生在当今构造、地震远不如相邻地带的龙门山断裂系上,震间和震后在地表形成了一条长300 km,宽近80 km,发生了70000多次余震的长廊地带,造成了严重的地表破坏和大量人员伤亡,而震前却未见到可视为与大地震发生与发展有因果关系的征兆现象,可是大地震就在这里发生了.多年的研究表明,这次M_S8.0大地震的发生乃是由于印度洋板块与欧亚板块碰撞、挤压,且东构造结向北"顶进"插入青藏高原东北缘力系作用下导致的该区构造、地震均强烈活动,在此背景下:(1)高原腹地壳、幔物质向东运移受阻、应力集中,深部物质重新分异、调整,物质与能量进行强烈交换;(2)震源深处介质属性与结构变异、破裂,且导致重力场的高度不均衡;(3)川西高原,相对低速壳、幔物质在坚硬的四川盆地壳、幔物质阻隔下,沿龙门山陡峭主断裂面相上逆冲,并在龙门山3条不同角度西倾断裂向下收敛处强烈碰撞;(4)震源深处汇聚断裂的形成构成了汶川-映秀M_S8.0大地震的发震断裂.     

汶川地震对芦山地震孕震环境影响的数值模拟

利用有限元数值模拟方法,以汶川地震和芦山地震的野外地质调查、同震位移和汶川地震前形变观测结果为约束,分析研究了高海拔地形蓄积的重力位能作用下,中地壳低速层的存在以及中央断裂和前山断裂同时破裂的条件下汶川地震的发生对龙门山南段孕震环境的影响.模拟计算结果显示,在一定的模型壳幔介质属性及空间结构下,龙门山断裂带西侧低速层的存在、青藏高原向四川盆地过渡带的地形特征以及汶川地震时中央断裂和前山断裂同时破裂的条件,是控制汶川地震对芦山地震孕震环境影响程度的重要因素.其中,地形和低速层对汶川地震引发的芦山地震孕震环境库仑破裂应力变化的影响最为显著,在所选剖面的北东段,地形特征的影响更为明显,而汶川地震双破裂面的影响则较小.     

地球是否进入地震活跃期

四川芦山地震勾起了国人对2008年汶川大地震的回忆,再想到近年来印度尼西亚、智利、日本等国多次发生8～9级大地震的情况,人们不禁会想:地球进入地震活跃期了吗?美国地质勘探局地球物理学家杰西卡.特纳说,尽管近年来世界各地发生了多次震级较大的地震,但"我们只能说地球在运动,目前没有明显证据显示地震增加"。同属美国地质勘探局的专家伊丽莎白.科克伦提供的数据显示,全球平均每     

21世纪以来青藏高原大震前重力变化

青藏高原自21世纪以来连续发生了多次大震:2001年昆仑山口西8.1级、2008年新疆于田7.3级和四川汶川8.0级、2010年青海玉树7.1级、2013年四川芦山7.0级地震、2014年新疆于田7.3级和2017年四川九寨沟7.0级大震.这些地震前,中国地震局在青藏高原开展过多期流动重力观测,并观测到震中附近可靠的重力随时间的变化.本研究综合利用地面绝对重力、相对重力资料,通过对多种重力观测资料的整体处理分析,研究了青藏高原区域重力场变化及其与大震发生的关系.结果表明:(1)大震易发生在与构造活动块体边界有关联的重力变化正、负异常区过渡的高梯度带上,重力变化等值线的拐弯部位,构造活动块体边界是物质变迁和构造变形差异运动强烈的地带,易产生剧烈重力变化,积累应力应变而孕育地震.(2) 2001年昆仑山口西8.1级地震发生在重力变化高梯度带与东昆仑断裂带交汇附近, 2008年汶川Ms8.0和2013年芦山Ms7.0地震发生在龙门山断裂带的重力变化高梯度带上, 2008和2014年两次于田Ms7.3地震发生在与康西瓦断裂走向基本一致的重力变化高梯度带零值线及梯度带的拐弯部位,2010年玉树7.1级地震发生在重力变化高梯度带与甘孜-玉树断裂上, 2017年九寨沟7.0级地震发生在具有显著重力变化的塔藏断裂和岷江断裂构造活动断裂带附近.(3)重力资料对青藏高原2001年以来发生的7次大震均有较好反映,地震前震中区及其附近观测到明显的区域性重力异常及重力变化高梯度带,可能是地震孕育过程中观测到的重力前兆信息.根据重力资料显示的异常变化,对2008年于田Ms7.3和汶川Ms8.0、2013年芦山Ms7.0、2014年于田Ms7.3和2017年九寨沟Ms7.0地震均进行了较准确的中期预测,尤其是震中位置的判定.     

汶川大地震·成都

5.12汶川大地震发生后,人们十分关心一个问题:位于龙门山附近的成都市区距离震中比北川、青川近得多,为什么几乎毫发未伤? 这首先要从发震构造说起.成都市区在成都平原东缘,虽然平原内也有两条小断层,但它们和龙门山断裂带没有直接连通的关系,不会因为龙门山中发生地震而直接受影响.     

东喜马拉雅构造结上新世以来快速抬升的裂变径迹证据

研究区位于东喜马拉雅构造结的东北侧,近东西向喜马拉雅山与北东向念青唐古拉山、近南北向横断山的交汇地段.这里是印度板块与欧亚板块碰撞的前缘部位,壳层的挤压变形十分强烈,无论是地质边界,还是地形都围绕此构造结发生了最急剧的转折,塑造了世界上著名的高山峡谷地形,其平均切割深度达5000m左右,是探讨青藏高原隆升历史最有利的典型地区.研究区地质上属于冈底斯岛弧的东延部分,其南部主要由一套角闪岩相的变质岩组成,可能是岛弧或拉萨地块的基底,北部主要为古生代的沉积盖层,在基底和盖层中有大量花岗岩     

浅谈桥梁抗震设计的问题及其对策研究

中国是一个地震多发的国家,历史上唐山大地震、汶川大地震等强烈的地震灾害给中国人民造成了无可估量的损失.为了对抗自然灾害,尽可能地避免地震灾害造成的问题,文章针对桥梁抗震设计中的问题与对策进行了总结,旨在为桥梁抗震设计工作贡献一份力量.     

高频GPS测定的汶川Ms8.0级地震震时近场地表变形过程

以龙门山断裂带下盘的四川GPS连续观测网高频(1 Hz)GPS观测资料为基础, 得出了2008年5月12日发生在四川汶川的Ms8.0级大地震的震时近场地表形变过程. 结果显示: 地震在近场产生的最大形变量明显大于震后位移, 破裂带北段各站水平分量震时先向震中方向运动, 后转折垂直于破裂带方向运动, 南段各站水平分量形变相对较小且基本为可恢复性变形. 各站垂向均先下沉, 然后呈周期性上下起伏波动. 将高频GPS所得位移与强震仪所得结果进行对比分析, 发现震动初期两者具有较好的一致性, 但中后期虽然相位基本同步, 但其振幅存在10 cm左右的差异, 其具体原因还需进一步研究论证. 此次记录到的汶川大地震近场地表变形过程, 可为进一步研究地表破裂过程和地震波的传播方式提供非常有价值的基础资料.     

青藏高原东缘地震各向异性、应力及汶川地震影响

使用国家地震台网固定台站和区域流动地震台站资料,通过对最长达18年记录开展的剪切波分裂分析,得到青藏高原东缘构造域(包括龙门山断裂带)的上地壳各向异性空间分布和时间变化特征.通过快剪切波(快波)偏振和慢剪切波(慢波)时间延迟参数,获得了快波偏振的分区分布特征和汶川地震前后不同分区统计特征的时间变化.剪切波分裂参数受到应力场和断裂构造的影响,块体边界或断裂附近的快波偏振特征比块体内部更为复杂,导致一些分区显示两个优势方向.在汶川MS8.0地震前后,可以观察到一些台站或区域的剪切波分裂参数的变化,既有快波偏振方向的改变,也有慢波时间延迟的变化.地震前后有明显变化的台站,基本符合两种情况:一是震中距较近;二是处于大的断裂或构造边界附近.研究发现,汶川地震前后,龙门山断裂带区域的慢波时间延迟的降幅显著大于周边区域;龙门山断裂带北段的慢波时间延迟降幅,大于龙门山断裂带中段,大于龙门山断裂带南端与鲜水河断裂、安宁河断裂的交汇区,这反映了地壳应力及介质物性状态的变化.通过分析剪切波分裂参数,推断出龙门山断裂带域及周边区域的主压应力分布.本文同时还分析了芦山地震前后剪切波分裂特征的变化.研究认为,利用地震剪切波特性监测应力变化,进而可应用于地震应力预测研究.需要注意是,不同的构造部位对应力变化有不同的响应.     

鄂尔多斯地块地壳上地幔速度结构及构造意义

基于鄂尔多斯地块及周缘地区稠密分布的地震台站所记录的远震波形资料,采用FMTOMO走时层析成像方法获得了深至400 km的高分辨率地壳上地幔P波速度结构,并分析和探讨了研究区的深部构造与火山、强震活动及克拉通破坏等相关的地球动力学问题.地壳上地幔速度结构横向变化和纵向差异的构造含义如下:鄂尔多斯地块下方正异常速度的巨厚岩石圈是鄂尔多斯地块缺乏克拉通大规模破坏的证据;大同火山群下方深至400km的柱状低速异常表明大同火山的岩浆作用起源于地幔过渡带,由太平洋板块俯冲作用引起熔融的玄武质岩浆上涌造成;鄂尔多斯地块与四川盆地之间呈东西向分布的上地幔低速异常带,可能是青藏高原东北缘软流圈流出的通道;鄂尔多斯周缘地区深至上地幔的低速异常分布表明,青藏高原北东向推挤和太平洋板块西向俯冲的共同作用,促使鄂尔多斯地块整体的抬升和周缘拉张断陷带的形成,为强震的孕育和发生提供了良好的应力环境,导致鄂尔多斯地块周缘地区地震频繁发生.     

玉树M_S7.1级地震地表破裂与历史大地震

野外调查表明,青海玉树MS7.1级地震较清晰的同震地表破裂带由3条主破裂左阶组成,走向310°～320°,总长约31km,左旋走滑性质.另在隆宝镇东侧一带见有长约2km的雁列式张裂缝带,如以该点为破裂带的北端点,则破裂带总长约51km.地表破裂带由一系列挤压鼓包与张裂缝相间排列或雁列式裂缝组成,实测最大水平位错约1.8m.地表破裂带沿甘孜-玉树断裂展布,显示该断裂是此次地震的发震构造.甘孜-玉树断裂历史上记载过多次7级左右地震,古地震遗迹明显,具有短周期的大地震重复特征.玉树地震的孕育机制与汶川地震一样,都是青藏高原东扩、地块边界应力积累和释放的结果,不同的是玉树地震为巴颜喀拉地块与川滇块体向东不均匀挤出产生的左旋走滑型地震.     

休闲时光

青藏高原成因有新解青藏高原隆起之谜困扰国际地学界多年,如今成都民间学者、四川地球起源演化研究所所长杨槐对此作出了新的解释.被称为"世界屋脊"的青藏高原的成因,目前国内外普遍采用板块说观点,认为是印度大陆北漂,碰撞欧亚大陆所致.杨槐以大量新资料论证印度洋膨胀的事实,继而合理建立起相应力学模型及该模型在青藏地区的具体应用,从而揭示     

青藏高原东边缘地壳“管流”层的电磁探测证据

通过对青藏高原东边缘及其附近地区石棉-乐山剖面大地电磁资料的研究发现, 青藏高原东边缘带和四川地块的电性结构有明显差别, 东部的四川地块地壳总体电阻率大, 西部的青藏高原东边缘带地壳总体电阻率小. 西部地壳分为上中下3层, 中地壳为厚约10~15 km的低阻层, 电阻率最小达3~10 Ωm, 推测它含有较低黏滞度的部分熔融和(或)含盐流体, 易于变形和流动, 是青藏高原东边缘带向东南方向挤出作用下形成的“管流”层. 它使上地壳和下地壳解耦, 上地壳高阻的脆性层以左旋走滑和逆冲的断层运动为主, 地表抬升, 地震主要发生在上地壳内. 低阻层把鲜水河-安宁河深大断裂带截成上下两段. 地壳上、中、下层各自厚度的横向变化, 综合产生了青藏高原东边缘带地壳厚度西部厚、东部薄、地形西部高、东部低的过渡带.     

青藏高原GPS位移绕喜马拉雅东构造结顺时针旋转成因的数值模拟

GPS测量结果表明青藏高原水平位移场在其东部绕喜马拉雅东构造结呈顺时针转动, 其地球动力学成因难以在二维模型中模拟. 由于青藏高原地壳巨厚、下地壳温度高, 岩石层流变性质可能会呈现脆性的上地壳、柔性的下地壳和黏滞系数较大的岩石圈上地幔, 在岩石圈下又存在黏滞系数较低的软流层的复杂分层流变结构. 在上述研究成果的基础上, 利用三维球坐标系下的Maxwell体有限元模型模拟了青藏高原在印度板块推挤下的变形. 模拟结果显示, 柔软下地壳的存在使整个青藏高原在印度板块的推挤下表现为整体抬升, 高原南缘和喜马拉雅东构造结抬升迅速, 而由于周边地块下地壳相对较硬而封闭、仅高原东南部存在高温柔软的通道, 青藏高原软的下地壳的存在还使得高原整体隆起达到一定高度后, 下地壳和软流层的物质向东、东南流动, 并拖曳上地壳作类似运动, 形成绕喜马拉雅东构造结的顺时针转动.     

GPS观测的2015年尼泊尔M_S 8.1级地震震前应变积累及同震变形特征

2015年4月25日,在印度板块与欧亚板块交界区的喜马拉雅地震带上发生了尼泊尔MS 8.1级大地震.震前GPS速度场和应变率场显示,喜马拉雅地震带整体表现为15.94±1.82 mm/a的压缩特征,同时还具有分段活动特征.此次地震发生在速度场顺时针旋转和逆时针旋转的分界带,该处最大主压应变率的量值在喜马拉雅地震带并非最大.GPS观测的同震位移场揭示了尼泊尔MS 8.1级地震引起的地壳变形特征,分别有9和6个测站观测到明显的水平向和垂向同震位移,其水平分量的运动方向整体表现为南向运动,位于震中东南侧的3个测站垂直分量表现为上升,其余测站为下降.中国境内距离震中最近的5个测站的垂向同震位移显示,此次地震造成珠穆朗玛峰的沉降量约为4mm.依据GPS观测到的同震位移场,利用非负最小二乘方法反演震源断层面上的滑动分布.反演结果表明最大滑动量为6.84 m,滑动量较大的区域分布在加德满都附近及其以北区域的下方,这可能是造成加德满都地区具有较大破坏的原因之一,该滑动分布模型能够很好地解释GPS观测到的同震位移.利用此滑动分布模型计算的地震矩为8.21×1020 N m,对应的矩震级为MW7.9.     

滇西南地区地壳速度结构的区域波形反演

利用中国地震局地壳应力研究所于2010~2011年在滇西南地区布设的宽频带流动地震台所记录的三分量高质量波形资料,采用小生境遗传算法通过波形反演,获得了研究区P波速度模型.结果显示,滇西南地区地壳厚度为38~43 km,上地幔顶部平均速度为7.9 km/s,地壳速度存在明显的横向不均匀性.思茅地块的上地壳和下地壳厚度均呈现出明显的由南向北、自东向西增厚趋势,说明该地块地壳增厚由上地壳和下地壳共同增厚形成.腾冲地块中地壳和下地壳的平均厚度及速度均小于保山地块和思茅地块,思茅地块部分路径上的下地壳和上地幔速度均较低,可能与印度板块东向俯冲引起地幔热物质上涌等动力过程密切相关.这些结果对于认识青藏高原形成与演化动力过程提供重要地震学约束.     

2004年印度尼西亚9级大地震前的潜热通量异常

研究了2004年12月26日印度尼西亚苏门答腊岛西海岸发生的9级大地震前后地球表面潜热通量的时空演化过程, 发现如下几个现象: (1) 在印尼9级大地震前, 在震中及其附近地区, 出现了显著的潜热通量异常; (2) 地震前最大的潜热通量异常出现于缅甸小板块中段俯冲带上, 这里正是地震破裂带的中部, 也是余震集中的区域; (3) 异常发展具有由弱到强、由零散到集中的过程, 异常开始出现在缅甸小板块张性东边界及其东侧海盆, 然后向压性西边界(即海沟俯冲带)和印度洋发展, 地震后潜热通量异常消失. 震源系统是一个耗散系统, 地震前应力增强使系统与外界的物质和能量交换加激, 提高了海面和大气之间水热交换的速率, 这可能是导致潜热通量异常的原因.