三峡库区范家坪滑坡地表形变InSAR监测与综合分析

位于湖北省秭归县的范家坪滑坡是长江三峡库区干流上的大型岩质滑坡之一。阐述了高分辨率合成孔径雷达干涉测量技术(InSAR)监测滑坡地表形变的工作方法与技术体系,采用22景3m空间分辨率的TerraSAR-X数据,辅以人工反射体布设和GPS测量,对范家坪滑坡变形进行监测,发现滑坡处在缓慢匀速变形状态,其谭家河滑坡体的形变比木鱼包滑坡体更为强烈,形变最大处的雷达视线向形变速率达到300mm/a。通过综合分析滑坡区2012年大气降雨和长江水位资料,发现年度内范家坪滑坡变形受水位变化和大气降雨影响微弱。     

准噶尔盆地南缘古牧地背斜多期构造变形特征

古牧地背斜位于博格达山前阜康断裂带的北缘。为探讨博格达山隆升在沉积盆地内的构造响应特征,作者结合野外地表露头观测、地震反射特征和磷灰石裂变径迹分析,认为古牧地背斜垂向上可划分为深浅两大构造层:深层侏罗系西山窑组(J_2x)以下地层具有双重构造的特征,活动时间为150~130和100~80 Ma;浅层侏罗系西山窑组及以上地层具有断层突破型断层传播褶皱的特征,活动时间约为30 Ma。以上研究成果一方面反映出博格达山晚侏罗世以来多期构造活动的特点,另一方面反映古牧地背斜深层双重构造具有形成时间早、埋深相对浅的特点,是准噶尔盆地南缘有利的油气勘探领域。     

孟加拉盆地东部褶皱带差异构造变形及成因

基于野外地质调查和地震资料分析,证实了孟加拉盆地东部褶皱带发育底部和中部两套滑脱层。底部滑脱层在整个褶皱带均发育,而中部滑脱层仅发育于北段—中段的陆上和陆架区,在中段的陆坡-洋盆区及南段则不发育。中部滑脱层的发育受古陆坡富泥沉积控制。北段—中段主要发育滑脱褶皱及被断层突破的滑脱褶皱,而南段主要发育泥底辟和泥火山构造。孟加拉盆地东部褶皱带滑脱层系及构造样式的差异变形特征主要受区域挤压作用、沉积作用以及地层流体超压等方面存在着差异的共同影响。     

基于GPS观测研究2010年青海玉树MS7.1地震震后地壳形变特征及其机制

2010年4月14日青海玉树M_S7.1地震发生在青藏高原东南部甘孜-玉树地震带,在震后7~10天内,我们快速建立了由15个GPS测站组成的跨地震破裂带观测剖面,包括1个连续站,3个半连续站和11个流动站,对所有站进行了240多天的观测,获取了该次地震的震后形变时空特征.采用欧拉矢量和位错模型解算了背景速度场,并从GPS观测的形变场中扣除该分量.采用分层黏弹性位错模型计算余震引起的地表形变,结果表明余震对部分测站的位移造成不可忽视的影响.采用对数模型拟合位移时间序列,表明特征衰减时间为6.7±1.2天.利用最速下降法反演震后余滑时空分布,反演结果表明震后断层活动以左旋滑动为主,断层南盘具有少量的抬升.在空间分布上,余滑主要位于同震破裂区的两侧,西北侧的余滑几乎达到地表,而东南区的余滑基本在同震破裂区的下方,余滑最大的区域位于结古镇东南下方10~20 km的深度范围.随着震后离逝时间的增加,2个余滑区在空间上保持不变,余滑区的面积逐渐扩大.余滑的矩释放为（1.5~5.1）×10¹⁸Nm,相当于1个M_W6.1~6.4地震释放的能量.分层岩石圈黏弹性模型计算的地壳孔隙弹性反弹形变与地表观测值相差较大,不能解释观测到的震后变形.采用麦克斯维尔流变体模型计算下地壳和上地幔松弛引起的地表形变,显示出其对地表形变的贡献较小.GPS观测得到的震后形变所具有的快速衰减特征,以及余滑模型能够较好地拟合GPS地表形变,表明2010年玉树M_S7.1地震后早期阶段的地壳形变主要是由余滑机制决定的.     

南北构造带北段上地幔各向异性特征

对布设在南北构造带北段的中国地震科学探测台阵项目二期674个宽频带流动台站和鄂尔多斯台阵21个宽频带流动台站记录的远震XKS(SKS、SKKS和PKS)波形资料作偏振分析,采用最小切向能量的网格搜索法和"叠加"分析方法求得每一个台站的XKS波的快波偏振方向和快、慢波的时间延迟,并结合该区域出版的122个固定台站的分裂结果,获得了南北构造带北段上地幔各向异性图像.快波方向分布显示青藏高原东北缘、阿拉善块体和鄂尔多斯块体西缘的快波方向主要表现为NW—SE方向,秦岭造山带的快波方向为近E—W方向,鄂尔多斯块体内部的快波方向在北部为近N—S方向,南部表现为近E—W方向.时间延迟分布来看,鄂尔多斯块体的时间延迟不仅明显小于其周缘地区,而且小于其他构造单元,特别是在高原东北缘、阿拉善块体和鄂尔多斯块体的交汇地区的时间延迟很大,反映了构造稳定单元的时间延迟小于构造活跃单元.通过比较快波方向的横波分裂测量值与地表变形场模拟的预测值,并结合研究区地质构造和岩石圈结构特征分析表明,在青藏高原东北缘、阿拉善块体和鄂尔多斯块体西缘各向异性主要由岩石圈变形引起,地表变形与地幔变形一致,地壳耦合于地幔,是一种垂直连贯变形模式;秦岭造山带的各向异性不仅来自于岩石圈,而且其岩石圈板块驱动的软流圈地幔流作用不可忽视;鄂尔多斯块体内部深浅变形不一致,具有弱的各向异性、厚的岩石圈和构造稳定的特征,我们认为其各向异性可能保留了古老克拉通的"化石"各向异性.     

用跨断层形变资料分析鲜水河断裂西北段的运动特征

本文利用主成分分析、震间位错反演和小波技术分析了鲜水河断裂西北段的跨断层形变资料(1986—2013).结果表明:以左旋走滑为主的断层长期运动为跨断层资料的主要信息,且符合负指数函数的运动规律,随着断层深度的增加,滑动量逐渐减弱.从炉霍段、道孚段到乾宁段,断层闭锁程度逐渐增强.2001年昆仑山M8.1地震发生后,鲜水河断裂西北段地壳浅层(地表以下30km)的左旋走滑明显减弱,到2013年的累积减弱量为3~13mm.5·12汶川M8.0地震发生前,断层滑动出现周期4~5年的增强信号.4·20芦山M7.0地震发生前,断层滑动首先出现低频信号(4~5年周期)增强,随着地震发生临近,强信号频率逐渐升高,直到出现周期为1年的强信号.     

青藏高原地壳密度变形带及构造分区

将区域重力场多尺度刻痕分析用于提取青藏高原地壳变形带的信息,可了解高原内地壳变形带从浅到深的变化和平面分布特征,并对青藏高原主要地体的空间分布定位,为岩石圈研究提供地表地质难以取得的新信息.多尺度脊形化系数的图像刻划不同深度平面上的地壳变形带.青藏高原地壳变形带从上到下由细密逐渐变为粗稀型,而且细密型变形区分布的范围逐渐缩小,到下地壳完全消失.从这种情况可以推测,以垂直地面方向上看,地壳变形带应该是树形的,下地壳粗稀型的变形带为树的主干,而中地壳粗稀型的变形带为树的分枝,上地壳的变形带为树枝的小枝杈.上地壳细密型变形分布区反映了与中新生代地壳缩短变形区的范围,下地壳清晰连续的变形带反映了青藏高原的构造骨架.多尺度边界刻痕系数的图像刻画不同深度平面上的地体边界,下地壳的刻痕边界系数与密度剧烈变化带位置吻合;因此,由多尺度刻痕分析划分地体时同时取得地体密度信息.青藏高原内密度较高的地体包括喜马拉雅地体、克什米亚地体、察隅河地体、柴达木地体、巴颜喀拉地体和羌塘地体.柴达木地体、巴颜喀拉地体和羌塘地体是青藏高原中有壳根的核,而密度最高的克什米亚和察隅河地体在大陆碰撞时不易碎裂,对东西两个构造结的形成起了关键作用.     

基于形变观测分析2011年日本9.0级地震与断层运动间关系

2011年3月11日日本发生9.0级地震,本文以此次地震的震间、同震和震后形变观测为约束,依据不同时段断层运动空间分布特征分析日本海沟地区强震与断层运动间关系.震间日本海沟地区,断层运动闭锁线深度约为60km,闭锁线以上从深到浅依次为断层运动强闭锁段、无震滑移段和弱闭锁段.由同震位错反演结果,2011年日本9.0级地震同震存在深浅两个滑移极值区,同震较浅的滑移极值区(同震位错量10~50m,深度小于30km)震间为断层弱闭锁段;同震较深的滑移极值区(同震位错量10~20m,深度在40km左右)震间为断层强闭锁段;而在两者之间的过渡带同震位错相对较小,震间断层运动表现为无震滑移.震后初期断层运动主要分布在在闭锁线以上的同震较深滑移极值区,而同震较浅的滑移极值区能量释放比较彻底,断层震后余滑量相对较小.依据本文同震和震间断层运动反演结果,震间强闭锁段积累10m同震位错需要100多年时间,与该区域历史上7级地震活动复发周期相当;震间弱闭锁段积累30~50m同震位错约需要300~600年时间,与相关研究给出的日本海沟9级左右地震复发周期比较一致.在实际孕震能力判定的工作中,由于不同性质的断层段在同震过程中会表现更多的组合形式,断层发震能力判定结果存在更多的不确定性,但利用区域形变观测等资料给出震间断层运动特征的研究工作对于断层强震发震能力的判定具有非常重要的实际意义.     

基于水准剖面的山西断陷带现今活动性研究

根据跨山西断陷带的水准剖面观测资料,计算其垂直形变量和垂直形变速率,研究其现今垂直活动性。2006年以来的垂直形变揭示,山西断陷带东北段延怀盆地和蔚县盆地呈掀斜上升活动,北段阳原盆地和大同盆地以下降活动为主,晋中盆地为上下波动,临汾盆地呈现上升活动。2013年山西断陷带活动出现转折性变化,蔚县盆地大幅度下降,大同盆地、代县-定襄盆地、晋中盆地和临汾盆地均呈现不同程度的上升。     

CRInSAR与PSInSAR融合探测地表线性变形速率

提出利用CRInSAR与PSInSAR融合解算地表变形参数的算法。将CR点上获取的线性变形速率与高程改正结果作为PS基线网络的起算数据,依据最小二乘原理求解PS点上待求线性变形速率与高程改正值的最优解。这一算法可融合PSInSAR与CRInSAR两种算法的优势,起到很好的互补作用。实验显示,这一算法获取的线性变形速率精度可达±0.37 mm/a,高程改正值精度可达±0.5 m,证实该算法在实际地表变形监测中具备可行性。