汶川地震区地壳上地幔三维S波速度结构初步研究

2008年5月12日汶川MW7.9地震发生在龙门山断裂带。龙门山断裂带及其邻域的地壳上地幔三维速度结构的研究对于理解汶川大地震的动力学背景具有重要的意义。2006年10月至2009年10月,在国家重大基础研究项目(973)的支持下,中国地震局地质研究所地震动力学国家重点实验室在川西地区(26°~32°N,100°~105°E)布设了由297台宽频带数字地震仪组成的流动观测台阵(简称川西台阵)。根据川西台阵记录的环境噪声和远震波形数据,利用噪声成像技术和接收函数方法,我们研究了川西地区(29°~32°N,100°~105°E)地壳上地幔100km深度范围内的三维S波速度结构。本文得到的结果为研究川西高原和四川盆地的地壳结构提供了新的高分辨率观测证据。我们的结果表明:1)观测台阵覆盖的川滇地块、松潘-甘孜地块和四川盆地的地壳上地幔S波速度结构具有显着差异,龙门山断裂和鲜水河断裂带,作为地块间的边界断裂带,对两侧地壳结构具有明显的控制作用。2)观测台阵覆盖区域的地壳厚度存在明显差异,川滇地块的地壳厚度为60~64km,松潘-甘孜地块的地壳厚度为52~56km,四川盆地前陆的地壳厚度为46~52km,沿龙门山断裂带松潘-甘孜地块和四川盆地形成镶嵌结构,汶川地震震中处南北两侧的壳幔边界存在约6km的断错。3)四川盆地前陆低速特征表明相应区域存在厚度8~10km的沉积盖层,松潘-甘孜地块和川滇地块的中下地壳具有大面积分布的S波低速区,松潘-甘孜地块地壳平均泊松比高达0.29~0.31,汶川地震余震绝大多数分布在低速区上方的高速介质区域内,而四川盆地的中下地壳呈现整体性的高速特征,以汶川地震的震中为界,龙门山断裂带北段和南段的S波速度结构显示了明显的速度分段特征,其北段的S波速度总体上高于南段。4)本文给出的研究区地壳三维S波速度结构表明,川西高原中下地壳较为软弱,而四川盆地中下地壳的强度应明显高于松潘-甘孜地块,意味着四川盆地坚硬中下地壳可以阻挡松潘-甘孜地块向东的逃逸;另一方面,川西高原和川滇地块的中下地壳虽然均存在大面积的S波低速区,但松潘-甘孜地块内的地壳速度结构相对来说较为复杂,并形成了高、低速相间的结构特征,表明在四川盆地的阻挡作用下,该地块形成了折皱变形的结构。5)与S波低速区相应,松潘-甘孜地块和川滇地块中下地壳应处于部分熔融的状态,这对该区域存在中下地壳通道流(Channelflow)的推断是一个支持;但是,松潘-甘孜地块内是否存在中下地壳通道流仍有待进一步的深入研究。6)接收函数方位各向异性的偏振分析表明,以汶川地震震中为界,龙门山断裂西南侧处于挤压状态,而其东北侧的主压应力方向与断层走向大体平行,推断先存应力场可能驱动了汶川地震逆冲破裂之后沿龙门山断裂向北东方向的走滑破裂。     

基于2001年MW7.8可可西里地震震后形变模拟研究藏北地区岩石圈流变学结构

青藏高原岩石圈的流变学结构和形变机制是地学界长期争论的重大科学问题.2001年发生在东昆仑断裂带的M_W7.8可可西里地震造成青藏高原北部地区岩石圈构造应力场的很大改变,引起下地壳与上地幔的快速弛豫形变,从而为研究这一问题提供了难得的机会.本研究采用该区域的GPS震后观测,反演这一地区岩石圈的流变学参数并探讨其形变机制.反演所采用的数据来自45个GPS观测点,其中包括一个中国地壳运动观测网络的基准站,数据最长时间跨度达6.4年.大地震震后形变场主要来源于地壳、上地幔的黏弹性松弛与断层面上的震后余滑,因此本研究同时反演介质的黏滞系数和断层的震后余滑.考虑到东昆仑断层南侧的巴颜喀拉-羌塘地区与北侧的柴达木盆地地区具有明显不同的地壳结构,断层南北两侧采用不同的Burgers体流变学结构,其下地壳-上地幔的短期和长期黏滞系数采用网格搜索法获得;断层震后余滑反演则同时施加近似正比于库仑应力的约束.最终结果显示：东昆仑断层北侧柴达木盆地地区下地壳-上地幔短期和长期黏滞系数分别为5×10¹⁸ Pa·s和1.5×10²⁰ Pa·s;东昆仑断层南侧巴颜喀拉-羌塘地区下地壳-上地幔短期和长期黏滞系数分别为1.5×10¹⁸ Pa·s和1.5×10¹⁹ Pa·s.这一结果表明：巴颜喀拉-羌塘地区下地壳-上地幔黏滞系数显著低于柴达木盆地,意味着巴颜喀拉-羌塘地区下地壳可能存在部分熔融,其地壳形变模式更趋近于连续形变,而柴达木盆地形变模式更趋近于块体运动.研究区下地壳长期黏滞系数比下地壳流模型所主张的黏滞系数高2~3个数量级,表明下地壳流在本地区可能不存在.     

岩石圈结构对大地震震后形变的影响——以1976年唐山大地震和2001年昆仑山大地震为例

大地震发生之后通常会诱发一系列的余震序列,对比1976年M_S7.8唐山大地震和2001年M_S8.1昆仑山大地震周边区域的地震事件可以看出,唐山大地震余震活动时间要明显长于昆仑山大地震余震活动时间.余震序列往往与震后形变密切相关,而影响震后形变的因素不仅与地震发震断层和震级有关,同时与岩石圈的结构有关.考虑到唐山大地震的发震区华北地块和昆仑山大地震的发震区青藏高原有着较大的岩石圈结构差异,本文采用PSGRN/PSCMP软件计算了岩石圈分层模型的大地震同震和震后形变,分析了地壳弹性模量、弹性厚度以及黏滞性系数对同震和震后形变的影响,进而讨论了影响唐山地震和昆仑山地震余震序列差异的原因.计算结果显示,震后形变会在黏弹性效应的作用下逐渐调整,震后形变的持续时间与地壳弹性模量、地壳弹性厚度和下地壳黏滞性系数有关.上地壳和下地壳弹性模量越大,震后形变达到稳定值的时间越短,弹性模量对震后形变稳定值影响很小.地壳弹性厚度越大,震后形变达到稳定值的时间越短,当断层面底端深度小于地壳弹性厚度时,地壳弹性厚度的增加会引起震后形变稳定值的减小;下地壳厚度对震后形变达到稳定值的时间和稳定值基本无影响.下地壳黏滞性系数越大,震后形变达到稳定值的时间越长,反之亦然.结合唐山地震区的华北地块和昆仑山地震的青藏高原深部结构发现,两者之间的上地壳弹性模型差别不大,唐山地震区地壳弹性厚度略大于昆仑山地震区,但昆仑山地震区下地壳黏滞性系数明显低于唐山地震区.这些因素均决定了昆仑山地震的震后形变持续时间短（余震时间序列短）而唐山地震的震后形变持续时间长（余震时间序列长）.由此可见,岩石圈结构差异可能是导致唐山地震和昆仑山地震余震序列差异的主要因素之一.     

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

利用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年, 即汶川地震后, 龙门山断裂中北段为逆冲兼右旋走滑运动; 而南段断层不能自由滑动、 变形宽度较大. 综合分析认为, 汶川地震时, 龙门山断裂南段并没有发生破裂, 一直处于较强的闭锁状态, 汶川地震的发生又加速了芦山地震的孕育进程; 由于龙门山断裂带南段的闭锁深度较中北段浅, 因此芦山地震较汶川地震强度低、 震级小、 破裂范围窄.      

从震后形变探讨青藏高原下地壳黏滞系数

青藏高原下地壳黏滞系数究竟是多少,已成为深入定量研究中的突出问题,它的数值量级将极大地影响定量模拟的结果．为了获得青藏高原下地壳的黏滞系数,从3条途径对该参数进行了计算：一是基于对青藏高原深部温度状态改进了的估计,用流变定律和GPS求得的应变速率对高原北部下地壳黏滞系数重新进行了估算,获得的昆仑山地区中地壳等效黏滞系数为10^20～10^22Pa．S量级,下地壳等效黏滞系数在10^19～10^21 Pa．s之间．二是用3种流变模型对2001年昆仑山Ms8．1级地震震后跨断层GPS站点记录到的震后变形进行了模拟,得到的下地壳黏滞系数为10^17 Pa．s量级．三是用黏弹性模型对炉霍地震后的跨断层形变曲线进行了拟合,得到的下地壳黏滞系数为10伸Pa．s量级．前人研究等效黏滞系数时忽视了等效黏滞系数与应变速率存在的非线性关系,本研究结果协调地解释了实验室实验、大地震后较短时期的变形和大地震后较长时间变形下,其等效黏滞系数存在差异的问题．     

青藏高原东缘的地壳流及动力过程

黏滞性地壳流对地壳及上地幔变形作用及动力机制,是大陆新生代造山带的一个重要研究内容.青藏高原中下地壳存在部分熔融或含水物质的黏滞性流体,已为一系列地球物理及岩石学研究所证实.为研究青藏高原东缘地壳流的动力作用,本文用密集的被动源宽频带地震台的观测数据,反演了地壳上地幔精细速度结构和泊松比.研究表明,川西及滇西北高原的中地壳内普遍存在低速层,而高泊松比的地壳只分布在川西北地区.位于中地壳的黏滞性地壳流从青藏高原腹地羌塘高原流出,自北西向南东流入青藏高原东缘.这些黏滞性地壳流带动了上地壳块体水平移动,当它们受到刚强的四川盆地及华南地块阻挡时将发生分层作用,地壳流将分为二或更多分支不同方向的分流,向上的一支地壳流将对上地壳产生挤压,引起地面隆升,向下的一支地壳流将使莫霍面下沉加厚下地壳·黏滞性地壳流的运动在地壳中产生应变破裂发生强烈地震活动,地震的空间分布与震源机制也受到地壳流动力作用控制.     

GPS观测揭示的芦山M_S7.0地震前龙门山断裂带南段变形演化特征

横跨龙门山断裂带南段的连续GPS测网记录到了2013年4月20日芦山MS7.0地震孕育过程相关的地壳变形信息,为研究此次地震前孕震区地壳变形动态演化过程提供重要的基础资料.研究表明,汶川地震的发生导致茂县-汶川断裂南段及以东地区挤压应变和左旋剪切应变加载.GPS跨单条断裂的基线平均缩短速率约为1~2 mm/a,跨越整个断裂带的基线平均缩短速率约为8~10 mm/a,且均表现出随芦山地震临近年均缩短速率逐渐减小的特征;多站组合的应变参数时序结果显示,龙门山断裂带南段主压应变率自西向东逐渐减小,主压应变方向为N30°~45°W近似垂直于断裂带;北川-映秀断裂以东地区以挤压变形为主兼有明显的左旋剪切变形,且面应变和第一剪应变随着芦山地震的临近应变率逐渐减小;北川-映秀断裂以西则表现为在时间进程上逐渐增强的右旋剪切变形.区域GPS变形场结果显示汶川震后龙门山断裂带南段挤压应变积累速率显著大于震前,且茂县-汶川断裂以东地区表现出左旋剪切应变积累特征.综合分析认为,汶川地震后巴颜喀拉块体东向运动加速,运动速度自西向东递减,致使在汶川地震中未破裂的龙门山断裂带南段的挤压应变积累水平进一步增强.     

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

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

川西地区台阵环境噪声瑞利波相速度层析成像

2006年中国地震局地质研究所地震动力学国家重点实验室在川西地区(26°N~32°N,100°E~105°E)布设了由297台宽频带数字地震仪组成的流动观测台阵.利用该密集台阵29°N以北156个台站2007年1~12月份的地震环境噪声记录和互相关技术,我们得到了所有台站对的面波经验格林函数和瑞利波相速度频散曲线,并进一步反演得到了观测台阵下方2~35 s周期的瑞利波相速度分布图像．本文结果表明,观测台阵覆盖的川滇地块、松潘-甘孜地块和四川盆地的地壳速度结构存在显著差异,具体表现为：(1)短周期（2~8 s）相速度分布与地表构造特征相吻合,作为川滇地块、松潘-甘孜地块和四川盆地之间的边界断裂,龙门山断裂带和鲜水河断裂带对上述三个地块上地壳的速度结构具有明显的控制作用,四川盆地前陆低速特征表明相应区域存在较厚的（约10 km）沉积盖层;(2)中周期（12~18 s）相速度分布表明,川滇地块和松潘-甘孜地块中上地壳速度结构存在明显的不均匀横向变化,并形成了尺度不同且高、低速相间的分块结构,而四川盆地中地壳整体上已经表现出相对高速;(3)长周期（25~35 s）相速度分布表明,松潘-甘孜地块,特别是川滇地块中下地壳表现为广泛的明显低速异常,意味着它们的中下地壳相对软弱,而四川盆地的中下地壳呈现整体性的相对高速,意味着四川盆地具有相对坚硬的中下地壳,并且以汶川地震的震中为界,龙门山断裂带的地壳结构显示了北段为高速异常,南段为低速异常的分段特征．     

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

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

阿坝—遂宁宽角地震剖面重建藏东缘龙门山地区地壳速度结构

龙门山断裂带位于青藏高原东缘,在中生代和晚新生代经历强烈的构造变形,急剧抬升,是研究青藏高原隆升和扩展动力学过程的重要窗口.本文利用起伏地形下的高精度成像方法,对"阿坝一龙门山一遂宁"宽角反射/折射地震数据重新处理,通过走时反演重建研究区地壳速度结构.剖面自西向东跨越松潘一甘孜块体、龙门山断裂带和四川盆地,不同块体速度结构表现了显著的差异.松潘甘孜块体地表复理石沉积层内有高速岩体侵入,低速层低界面起伏不平反映了该区的逆冲推覆构造.中下地壳速度横向上连续变化,平均速度较低(约6.26 km·s~(-1)).四川盆地沉积层西厚东薄,并在西侧出现与挤压和剥蚀作用相关的压扭形态.中下地壳西薄东厚,平均速度较高(约6.39 km.s~(-1)).龙门山断裂带是地壳速度和厚度的陡变带,Moho面自西向东抬升约13 km.在整个剖面上Moho面表现为韧性挠曲,中下地壳横向上连续变化,推测古扬子块体已到达松潘甘孜块体下方.松潘甘孜块体下方中下地壳韧性变形,并在底部拖曳着被断裂切割的脆性上地壳,应力在不同断裂上积累和释放,诱发大量地震.     

汶川Ms8．0地震：地壳上地幔S波速度结构的初步研究

2008年5月12日我国四川省汶川地区发生了震惊世界的Ms8.0地震.历史上,同类地震在大陆内部极为罕见.该地震深部构造背景的研究对理解其成因极为重要.本文利用中国地震局地质研究所地震动力学国家重点实验室在川西地区布设的大规模密集流动宽频带地震台阵记录的远震P波波形数据和接收函数非线性反演方法,得到了沿北纬31°线的19个台站下方120 km深度范围内的S泼速度结构及台站下方地壳的平均泊松比.该观测剖面穿越了主震区,总长度约为420 km.我们的结果揭示了川滇地块、松潘-甘孜地块和四川盆地三个不同地块构造差异.上述三个地块的地壳结构特征口J以概括为:(1)四川盆地前陆壳幔界面向西侧倾斜井有较为明显的横向变肜,地壳厚度存在46～52 km的横向变化,中下地壳S波速度存在横向变化,地壳平均泊松比值较高(0.28～0.31),但在龙门山断裂带附近,显示了坚硬地壳的特征,地壳平均泊松比仪为0.2;(2)松潘-甘孜地块地壳厚度由西侧靠近鲜水河断裂的60 km,向东减薄为52 km,在14～50 km深度范围内存在S波速度2.75～3.15 km/s的楔状低速区,其厚度由西侧的～30 km向东逐渐减薄为～15 km,相应区域的地壳平均泊松比高达0.29～0.31;(3)鲜水河断裂西侧,川滇地块地壳结构相对简单,地壳厚度为58 km,并在26 km深度存在约10 km厚度的高速层,地壳内平均泊松比约为0.25;(4)汶川大震区在12～23 km深度上具有近乎4.0 km/s的S波高速结构,而其下方的地壳为低速结构,地壳平均泊松比0.31～0.32,汶川大震的余震序列主要分布在高速介质区域内.本文的结果表明松潘-甘孜地块的地壳相对软弱;而且并不仔在四川盆地向西侧的俯冲,我们认为在青藏高原东向挤压的长期作用下,四川盆地强硬地壳的阻挡作用可导致松潘-甘孜地块内部蓄积很大的应变能量以及上、下地壳存壳内低速层顶部边界的解耦,在龙门山断裂带附近形成上地壳的铲形逆冲推覆.汶川大地震及其邻近区域所具有的坚硬上地壳和四川盆地的阻挡作用为低应变率下的高慢度应力积累创造了必要条件,而松潘-甘孜地块长期变形积累的高应变能构成了孕育汶川大地震的动力来源.     

汶川MS8.0地震：地壳上地幔S波速度结构的初步研究

2008年5月12日我国四川省汶川地区发生了震惊世界的M_S8.0地震.历史上,同类地震在大陆内部极为罕见.该地震深部构造背景的研究对理解其成因极为重要.本文利用中国地震局地质研究所地震动力学国家重点实验室在川西地区布设的大规模密集流动宽频带地震台阵记录的远震P波波形数据和接收函数非线性反演方法,得到了沿北纬31°线的19个台站下方120 km深度范围内的S波速度结构及台站下方地壳的平均泊松比.该观测剖面穿越了主震区,总长度约为420 km. 我们的结果揭示了川滇地块、松潘-甘孜地块和四川盆地三个不同地块构造差异.上述三个地块的地壳结构特征可以概括为：（1）四川盆地前陆壳幔界面向西侧倾斜并有较为明显的横向变形,地壳厚度存在46～52 km的横向变化,中下地壳S波速度存在横向变化,地壳平均泊松比值较高（0.28～0.31）,但在龙门山断裂带附近,显示了坚硬地壳的特征,地壳平均泊松比仅为0.2;（2）松潘-甘孜地块地壳厚度由西侧靠近鲜水河断裂的60 km,向东减薄为52 km,在14～50 km深度范围内存在S波速度2.75～3.15 km/s的楔状低速区,其厚度由西侧的～30 km向东逐渐减薄为～15 km,相应区域的地壳平均泊松比高达0.29～0.31; （3）鲜水河断裂西侧,川滇地块地壳结构相对简单,地壳厚度为58 km,并在26 km深度存在约10 km厚度的高速层,地壳内平均泊松比约为0.25;（4）汶川大震区在12～23 km深度上具有近乎4.0 km/s的S波高速结构,而其下方的地壳为低速结构,地壳平均泊松比0.31～0.32,汶川大震的余震序列主要分布在高速介质区域内. 本文的结果表明松潘-甘孜地块的地壳相对软弱;而且并不存在四川盆地向西侧的俯冲.我们认为在青藏高原东向挤压的长期作用下,四川盆地强硬地壳的阻挡作用可导致松潘-甘孜地块内部蓄积很大的应变能量以及上、下地壳在壳内低速层顶部边界的解耦,在龙门山断裂带附近形成上地壳的铲形逆冲推覆.汶川大地震及其邻近区域所具有的坚硬上地壳和四川盆地的阻挡作用为低应变率下的高强度应力积累创造了必要条件,而松潘-甘孜地块长期变形积累的高应变能构成了孕育汶川大地震的动力来源.     

汶川地震震后短期青藏高原东北缘地壳变形特征成因讨论

2007—2009年时段GPS数据资料显示青藏高原东北缘地区地壳运动较前期(1999-2007年时段)发生明显变化.本文分别对汶川地震同震影响、震后粘弹性松弛影响和震后余滑影响进行了分析.结果表明,同震对该区域地壳变形特征的形成存在一定影响,而震后粘弹性松弛和余滑产生的影响均可忽略不计.分析认为,对于中国大陆这种多块体系统而言,现有的震后变形机制并不能充分解释震后短期某些区域地壳运动变形特征的形成,汶川地震震后短期青藏高原东北缘的地壳变形特征的形成可能主要与震后块体运动调整有关.     

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

本研究基于分层黏弹介质模型,考虑同震位错效应和震后黏滞松弛效应,分析巴颜喀拉地块东端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级地震的发生;汶川地震对九寨沟地震的影响研究中,采用不同的汶川地震同震位错模型,计算结果有差异.综合考虑多次大地震对周边断裂带的影响,龙门山断裂带南段、鲜水河断裂带中南段、平武—青川断裂北段、灌县—安县断裂北段、文县断裂的累积库仑应力增加显著,巴颜喀拉地块东端的东昆仑断裂带东段、迭部—白龙江断裂带西段以及金沙江断裂带库仑应力亦有所增加.综合考虑各重要断裂带已有的大地震危险性分析结果和库仑应力变化计算结果,龙门山断裂带南段、鲜水河断裂带中南段、东昆仑断裂带玛沁—玛曲段和金沙江断裂带的发震紧迫性有所增强,需引起关注.     

昆仑山Ms 8.1级地震震后变形场数值模拟与成因机理探讨

2001年11月14日,在青海和新疆交界处发生了昆仑山M_s8.1级强烈地震,GPS后观测显示,此次地震震后形变不仅在断裂南北两侧存在很大的差异,而且在短时间调整后断裂南北两侧表现为同向运动.本文以观测的地震形变为约束,通过有限元数值模拟分析昆仑山地震震后形变的物理机制.建立有关的有限元虚功方程,通过有限元数值方法模拟震后形变,从理论上分析介质的非均匀性、黏滞性松弛、流体调整对震后形变的影响.采用网格搜索确定昆仑断裂南北两侧下地壳的黏滞系数分别为5.0×10¹⁷Pa·s, 9.0×10¹⁸Pa·s左右,正是这十余倍的差异引起了断裂两侧震后形变的非对称性和同向运动,这一差异既是长期地质作用的结果,又是现代地球动力学环境的决定因素之一.通过数值模拟定性讨论了断裂北侧地表形变在震后短期内的调整,对于靠近断裂附近的测点可能是黏弹性松弛和孔隙流体调整共同作用的结果,所以在分析短期震后形变时综合考虑黏弹松弛和孔隙流体调整是很有必要的.     

金川—芦山—犍为剖面重力异常和地壳密度结构特征

重力剖面金川—芦山—犍穿越芦山震区,近垂直于龙门山断裂带南段,长约300km,测点距平均2.5km,采用高精度绝对重力控制下的相对重力联测与同址GPS三维坐标测量,获得了沿剖面的自由空气异常和布格重力异常,并对布格重力异常进行了剩余密度相关成像和密度分层结构正反演研究.结果表明,芦山地震所在的龙门山断裂带南段存在垂直断裂走向的宽广的巨型重力梯级带,重力变化达252×10-5 m·s-2以上(龙泉山以西),反映出四川盆地与松潘—甘孜地块地壳厚度陡变(约14.5km)性质;四川盆地与松潘—甘孜地块过渡区(龙门山断裂带与新津—成都—德阳断裂之间)存在(30~50)×10-5 m·s-2的剩余异常"凹陷",可能与上地壳低密度体、山前剥蚀与松散堆积和推覆体前缘较为破碎有关;剩余密度相关成像显示地壳密度呈现分段性特征,在芦山地震位置出现高低密度变化;地壳呈现三层结构,四川盆地上、中、下地壳底界面平缓,反映其稳定阻挡作用,而松潘—甘孜块体上、中、下地壳底界面明显往盆地逐步抬升,反映出青藏高原往东的强烈挤压作用;松潘—甘孜块体往东推覆变形主要集中在上地壳范围内,推覆深度随离龙门山断裂带愈近而越浅.本文通过对密度分布及结构特征的研究,分析了芦山地震及龙门山地区地壳构造背景和当前活动性的深部动力环境特征.     

大地测量约束下的阿尔泰山岩石圈流变结构

本文根据大地测量数据得到过去50年左右阿尔泰山富蕴区域形变场，破裂带中段的相对位移最大，平均速率达6mm/a，总体上表现为沿断裂的走滑运动特征．该形变场可以用1931年富蕴8级地震的震后黏弹性松弛模型进行模拟，反演得到下地壳黏滞系数为1.6×10¹⁹～7.9×10¹⁹Pa·s，上地幔黏滞系数为16×1018～63×1019Pa·s，与华北、Nevada等地区利用震后变形资料推算的黏度基本一致．根据该地区最佳黏弹性分层模型，最近五十年由于岩石圈下部应力松弛引起的地震破裂带两侧最大水平速率约为4mm/a． 我们的研究表明：大陆7～8级大震在几十年后仍可能有可观的地表变形，GPS监测得到的现今变形场可能包含震后变形成分．     

基于模拟退火法的青藏高原东缘噪声层析成像研究

青藏高原是全球造山带研究的热点地区,此前在青藏高原开展的三维层析成像研究大多基于线性反演方法.本文利用青藏高原东缘及邻区布设的127个宽频带固定地震台站记录的连续波形资料,首先通过噪声互相关提取了3～50sRayleigh波群速度频散曲线并反演得到群速度分布,再进一步采用模拟退火法反演了研究区的三维S波速度及泊松比结构.结果显示:(1)松潘—甘孜地块的中下地壳低速异常主要分布在龙日坝断裂带、鲜水河断裂带、龙门山断裂带和岷山隆起所围限的区域,而该区域的中下地壳仅具有中等泊松比值,推测松潘—甘孜地块中下地壳的低速物质可能是青藏高原与扬子块体长期相互作用产生的塑性低速滑脱层;上地壳脆性物质在板块作用下沿中地壳低速滑脱层顶界面发生逆冲增厚,造成龙门山的持续抬升和地形起伏,并在构造边界带形成了应变积累和应力集中;而龙门山断裂带的上地壳低速软弱物质为地壳发生破裂提供了有利条件,从而在某种程度上促进了汶川地震和芦山地震的发生.(2)岷山隆起一带中下地壳的高泊松比异常呈"凸起"形态,结合前人研究发现的较高热流和岩石快速抬升现象,推测岷山隆起一带可能存在岩石圈的拆沉,导致地幔热物质上涌而形成下地壳高泊松比物质.(3)川滇地块的北部和南部具有不同的S波速度和泊松比分布特征.30km深度下川滇地块北部具有明显的低速异常,而该深度下并不具有明显的高泊松比值特征;此外剖面成像结果也显示川滇地块内的低速异常与高泊松比的分布不一致,因此川滇地块的研究结果不支持下地壳流模型.综合其他地震学证据,本文认为川滇地块的变形模式为上地壳纯剪切增厚,块体变形主要受块体内部的走滑断裂及活动边界断裂控制.     

青藏高原东缘及四川盆地的壳幔导电性结构研究

自从2008年M_S8.0级汶川大地震发生以来,青藏高原东缘便成为地质与地球物理研究的热点区域.该区域的龙门山断裂带标志着青藏高原东缘与四川盆地的边界.汶川地震即发生于龙门山断裂带内的映秀-北川断裂上.该地区现有的研究工作多集中于青藏高原东缘及四川盆地的西部,对四川盆地东部构造情况的研究目前较少.在SinoProbe项目的资助下,完成了一条跨越青藏高原东缘及整个四川盆地的大地电磁测深剖面.该剖面自西北始于青藏高原内部的松潘-甘孜地块,向东南延伸穿过龙门山断裂带、四川盆地内部及四川盆地东部的华蓥山断裂,最终止于重庆东南的川东滑脱褶皱带附近.维性分析表明剖面数据整体二维性较好,通过二维反演得到了最终的电性结构模型.该模型表明,从电性结构上看,沿剖面可分为三个主要的电性结构单元,分别为:浅部高阻、中下地壳低阻的松潘-甘孜地块,浅部低阻、中下地壳相对高阻的四川盆地,以及华蓥山以东整体为高阻特征的扬子克拉通地块.龙门山断裂带在电性结构上表现为倾角较缓、北西倾向的逆冲低阻体,反映了青藏高原东缘相对四川盆地的推覆作用.其在地下向青藏高原内部延伸,深度约为20 km左右.在标志逆冲推覆滑脱面的低阻层下存在一电性梯度带,表征着低阻的青藏高原中下地壳与高阻的扬子地壳之间的电性转换.位于四川盆地东边界的华蓥山断裂在电性结构上表现为一倾向为南东向的低阻体插入高阻的扬子克拉通结晶基底,切割深度约为30 km左右.这一结构反映出华蓥山向西的推覆作用.在电性结构模型的基础上,进一步讨论了青藏高原东缘的壳内物质流、青藏块体与扬子块体的深部关系以及青藏高原东部的隆升机制等构造问题.