摘要：用三维流变非连续变形与有限元相结合（DDA＋FEM）的方法，在青藏高原及邻近地区三维构造块体相互制约的大背景中，考虑了龙门山断裂带东西两侧地势、地壳厚度和分层的明显变化，及断裂带东侧四川盆地及鄂尔多斯块体坚硬地壳阻挡的影响，通过用GPS资料做位移速率边界约束和震源机制约束，计算得到研究区的速度场和应力场与该地区GPS测量结果和震源机制分布结果基本一致．在此基础上，模拟计算现今构造块体边界断层上表征剪应力及法向应力等综合影响的危险度分布．结果表明，上、中地壳层危险度分布中危险度较高的地段多数与近几十年来发生的七级以上大震区域基本一致．包括2008年汶川8．0级等大震的发震断层．通过分别对龙门山断裂带东西两侧的两种不同构造格局进行试算表明，龙门山断裂带东西两侧地势、地壳厚度、分层与物性明显变化对汶川大震的孕育发生均起了关键性作用．计算得到的应变率强度分布图可见，高原东部整个边缘地带均接近应变率强度的陡变带．其中以龙门山断裂带上的陡变最为明显，西侧应变率强度是东侧的近4倍，而且断裂带东侧应变率强度等值线衰减比西侧快．反映了汶川大震逆冲型发震断层地区独特的特征．此外，由计算得到的应变能密度分布图可见，龙门山断裂带在上、中地壳层中均位于宽度相同、其走向与龙门山断裂带走向一致的高应变能密度带中，在上地壳层这个带的东西两侧则是应变能密度较低的地区，而在中地壳层，其强度在断裂带东侧逐渐向东衰减，西侧应变能密度高，而东侧应变能密度较低．表明在印度板块强烈推挤作用和高原各构造块体相互制约及龙门山断裂带东西两侧特殊构造环境中，高原地壳物质向东水平运动，受到龙门山断裂带东侧介质刚性强度较大的�

关键词：汶川地震 孕震机理 青藏高原 龙门山断裂带 非连续变形分析方法 

单位：中国科学院地质与地球物理研究所地球深部研究重点实验室 北京100029 中国地震局地球物理研究所 北京100081 成都理工大学地质灾害防治与地质环境保护国家重点实验室 成都610059