摘要：2008年5月12日汶川地震突发在现今并不活动的龙门山断裂带上，该地震发生的动力学机制问题引起广泛关注．文中利用黏弹性接触问题的有限元方法，考虑重力作用，对青藏高原东缘的应力场空间分布及其随时间的演化进行了数值模拟，结果显示应力在空间由分散分布逐渐向龙门山及周边地区转移集中．基于前人的研究成果及计算分析，初步认为汶川地震孕育发生的动力学过程如下：青藏高原的物质东流在向东运动过程中由于受到稳定的四川盆地的阻挡，一部分东流物质在川西地区囤积，造成龙门山隆升；高角度（50°～70°）、犁状的龙门山断层面上的正应力随着川西高原向东运动而不断增大，导致该断层的闭锁性逐步加强，并且分布在断层附近的变质杂岩为存贮高密度弹性应变能提供物质保障．但另一方面随着青藏高原较柔软的下地壳物质的不断向东运动，囤积的东流物质对龙门山断裂带上盘的推挤作用会不断加强，从而导致断裂带上剪应力越来越大；当剪应力超过摩擦强度时，断层解锁产生滑动，发生地震．模拟结果还表明龙门山断层面上的摩擦系数较高，断裂带上地震的平均复发周期约为3163年，这与其他资料结果有一致性．

关键词：黏弹性有限元 物质东流 动力学机制 汶川地震 

单位：中国地震局地壳应力研究所 北京100085 中国地震局地质研究所地震动力学国家重点实验室 北京100029