青藏高原东南缘构造应力场重构

采用有限元方法，针对青藏高原东南缘建立更细致、更精确的三维有限元弹性模型。选取9种不同的应力边界条件，分别进行优化分析后处理，将对应台站形变模拟值与GPS实测值进行误差分析，最终选取最佳方案作为古构造应力场。结果表明，青藏高原东南缘4 Ma BP的古应力场主要起源于中国大陆周围板块的相互作用，特别是印度板块NNE向强烈碰撞作用，成为中国大陆尤其是西南部青藏高原地区构造应力场最主要的动力来源，控制各个块体相互作用的方式和运动格局。青藏高原东南缘古应力场主要包括几个力源：西北部青藏高原侧向挤压造成的WE向应力约105 MPa；西南部直接来自于印度板块的NE-WS向应力约70 MPa；南部NS向作用力33 MPa；东南部扬子块体侧向NW-SE阻挡力56 MPa；北东部受扬子块体强烈EW向阻挡力90 MPa。这些力源共同作用于青藏高原东南缘，形成现今复杂应力场。
     

青藏高原东南缘壳幔力学耦合及其动力学意义

现代大陆动力学及构造学研究认为，壳内软弱带控制着地壳与地幔间相互作用的强弱，制约岩石圈层间的耦合关系，同时也影响地块的相互作用方式。这种耦合作用可以通过比较不同深度的形变和应力的相关特征来分析。青藏高原东南缘20～35 km的深度广泛存在低速带，有可能在现今地震及构造活动区内由壳内部分熔融承担了潜在的解耦功能。基于1999～2015年GPS坐标时间序列数据得到速度场，结合前人给出的SKS剪切波分裂数据，分析该区水平方向形变耦合特征；同时搜集莫霍面和岩石圈底面形变数据，给出垂向上的形变耦合特征。利用EGM2008重力模型的11～36阶球谐系数计算青藏高原东南缘岩石圈底面地幔对流应力场；搜集并整理2000年至今的1 131个震源机制解数据，采用区域应力张量阻尼反演法得到该区多震层应力场；进一步分析两种应力场的相关性。在ABAQUS软件中建立青藏高原东南缘岩石圈三维粘弹性有限元模型，同时加入上地壳的边界位移约束条件和底面地幔对流拖曳力加载条件，分别对水平方向和垂直方向的应力应变作分析，并给出强震剖面结果，探讨壳幔耦合强震孕育机制，据此分析青藏高原东南缘各个断层的地震危险性。研究认为，应力的不均匀发展是应变能积累的关键，同时垂向上圈层间的耦合关系极大地影响块体间的作用方式；壳内软弱带广泛发育的块体应力出现分层现象，边缘区域横向纵向应力变化明显，与之接触的上伏地壳地震危险性较高。     

青藏高原东南缘壳幔耦合强震孕育模型研究

结合青藏高原东南缘软弱带的分布及地震活动性的深度空间特征，提出强震孕育的壳幔耦合模型。利用中国国家数字地震台网和四川云南区域数字地震台网及布设在青藏高原东南缘的地震台站所记录的远震SKS波形资料，通过剪切波分裂得出上地幔各向异性结果数据，描绘出青藏高原东南缘的上地幔各向异性图像；再分别同该地区GPS观测的速度场数据、区域构造应力场数据和13 km深度的介质密度数据进行相关性处理与分析，并针对本文提出的孕震模型进行一定程度的验证。对形变场和应力场耦合特征的研究表明，各向异性的快慢波时间延迟、快波偏振方向总体分别与GPS 测量的地壳运动速度场方向变化和水平最大应力方向相关性较大，壳幔各圈层相互作用及软流圈的流动所展示的地震各向异性都与应力传递机制息息相关。通过相关分析可以得到壳内软弱带（低速高导体）在壳幔耦合运动研究中的关键作用，而该类耦合运动与地震孕育息息相关。     

青藏高原东南缘地幔对流与发震层应力场耦合关系分析

利用EGM2008重力模型的11~36阶球谐系数计算青藏高原东南缘岩石圈底面地幔对流应力场，收集整理青藏高原东南缘2000年至今的1 131个震源机制解数据，采用区域应力张量阻尼反演法得到该区芦山地震前和现今的发震层应力场，进一步分析2种应力场的相关性，探讨不同区域的力学耦合情况与强震发生的关系。结果表明：1)青藏高原东南缘大部分位于耦合与解耦的中间地带，耦合区域基本按块体分布，东侧华南块体强耦合，西北部藏北块体、巴颜喀拉块体和西南部滇西南块体部分耦合，解耦主要发生在松潘-甘孜块体附近，以龙门山断裂带连接强耦合的华南块体；2)分析孕震原因，提出青藏高原东南缘地幔对流应力场与发震层应力场耦合程度强弱交界处为地震危险性较高的区域，从现今耦合关系来看，龙门山断裂带仍处于耦合强度变化梯度非常大的区域，具备孕震储能条件，地震危险性较高。其他危险区域大致有：岷江断裂带附近、鲜水河-安宁河-则木河断裂带附近、红河断裂带及南汀河断裂带附近。