川滇地区构造块体系统稳定性影响的三维DDA+FEM方法数值模拟 3-D numerical simulation on influence of 1997 Mani earthquake occurrence to stability of tectonic blocks system in Qingzang and Chuandian zone using DDA+FEM method.pdfVIP

第51卷第5期 地 球 物 理 学 报 V01．51，No．5 2008年9月 CHINESE OFGEOPHYSICS Sep．，2008 JOURNAL 陈祖安，林邦慧，白武明．1997年玛尼地震对青藏川滇地区构造块体系统稳定性影响的三维DDA+FEM方法数值模拟．地 球物理学报，2008，51(5)：1422～1430 Chen B W numerical oninfluenceof1997 occurrencetO of ZA。LinH．BaiM．3-D simulation Mani stability earthquake tectonicblocks in andChuandianzone DDA+FEMmethod．Chinese Chinese)，2008， systemQingzang using J．Geophys．(in 51(5)：1422～1430 1997年玛尼地震对青藏川滇地区构造块体系统 稳定性影响的三维DDA+FEM方法数值模拟 陈祖安h3，林邦慧2，白武明1 1中国科学院地质与地球物理研究所地球深部研究重点实验室，北京 一 2中国地震局地球物理研究所，北京100081 3地质灾害防治与地质环境保护国家重点实验室成都理工大学，成都 一 擒 要 本文用三维非连续变形与有限元相结合(DDA+FEM)的方法，在青藏川滇地区三维构造块体相互制约的 大背景巾，通过用GPS资料做位移速率边界约束和震源机制约束，计算得到研究区的初始位移场和应力场与该地 区GPS测量结果和震源机制分布结果基本一致．在此基础上进一步数值模拟1997年玛尼7．9级大震的发生过程， 研究大震引起研究区各块体边界断层应力状态变化的特征．(1)发震断层两侧发生左旋走滑错动，最大水平位错大 约7m；(2)深部位错面上位错分布与用地震波资料震源反演的结果类似；(3)最大差应力变化等值线图与由星载 D-INSAR技术获取的地表形变场图像相似；(4)地表垂直位移表明地震断层面略向北逆冲．计算模拟得到了玛尼地 震发生引起青藏川滇地区构造块体系统各边界断层上库仑破裂应力变化的分布，表明玛尼大震的发生除了使其发 震断层的两端库仑破裂应力增大，应力进一步集中外，位于上地壳层上东昆仑断裂中段的2001年昆仑山8．1大震 (H一11 kin)发震断层段的库仑破裂应力增加约2MPa，位于中地壳层上喀拉昆仑断裂带中的2008年改则6．9级 地震(H=30km)发震断层段的库仑破裂应力也增加约0．7MPa，可见这两个已接近破裂强度地段的失稳对发生 大震起了一定促进作用．研究结果也表明：作者发展的三维DDA+FEM方法能有效地用于大震活动与各构造块体 相互作用关系的研究． 关键词玛尼大震，青藏川镇构造块体，非连续变形分析方法，数值模拟 中图分类号P541．P315 文章编号0001—5733(2008)05-1422—09 simulationoninfluenceof1997Mani 3-Dnumerical earthquake occurrenceto oftectonicblocks in stability systemQingzang method andChuandianzone DDA+FEM