青藏高原东北缘地震时空迁移的有限元数值模拟

地震在大陆内部断层系统中的时空迁移和丛集的基本力学机制一直是地球科学家关注的重要问题.青藏高原东北缘地震活动频繁,其地震时空迁移和地震丛集现象显著,是研究这个问题的重要区域.我们建立了一个三维黏弹塑性有限元模型,模拟了青藏高原东北缘主要活动断层系统的地震循环和地震时空迁移;计算了断层系统的应力演化;并探讨了断层之间的相互作用及地震时空迁移和地震丛集的原因.模拟结果显示断层之间的相互作用通过增加或降低断层上的库仑应力,加速或延缓了地震发生,使得区域地震可以在短时间内集中发生,从而形成地震丛集;另外,区域经过多个地震循环的长期演化,一些孕震断层上的应力状态恰好都达到屈服的临界状态附近,从而也可以导致这些断层上的地震在短期内集中发生,因此产生地震丛集和地震迁移.我们发现当区域经历地震丛集之后,该区域的应力大大释放,区域进入地震平静期;随着构造加载的持续,区域应力逐渐恢复,为下一次地震丛集或地震序列累积应力和能量;上述过程可以重复发生.因此地震丛集期与平静期交替出现.我们还统计了各个断层的大地震相互迁移的模拟结果,结果显示青藏高原东北缘下一次大地震有很大的概率会发生在海原断层上.

青藏高原东缘中下地壳流与地壳变形

地壳缩短导致青藏隆升造山是普遍的认识.然而,在青藏东部,越来越多的观测数据和研究支持了中下地壳流与隆升造山的关系.目前,地壳缩短造山机制和中下地壳流造山机制仍然处于争论之中.本文建立了二维黏弹塑性有限元模型,模拟了龙门山断层带的多个地震循环的应变与变形,探讨了无与有中下地壳流情况下,地壳地表的位移、速度与变形的分布和演化;以及有中下地壳流情况下,不同流动范围、速度与黏度对模型结果的影响;并结合地形变观测数据的约束,推测了青藏东缘中下地壳流的流动状态.模拟结果显示,通过对比有和无中下地壳流的模拟结果,发现青藏东部震间的地表垂向速度在变形样式及数值上存在较大差异,即存在地壳流的地表垂向抬升速率显著大于无地壳流;震间在龙门山断层西侧附近产生的垂向凸状隆起随中下地壳流的速度、黏度及通道长度的变化而变化.此外,本文研究结果对青藏其他地区可能存在的地壳流的研究也具有一定的参考意义.