地震发生过程的有限单元法模拟——以苏门答腊俯冲带上的大地震为例 Simulation of earthquake processes by finite element method The case of megathrust earthquakes on the Sumatra subduction zone.pdfVIP

第5l卷第2期 地球物理学报 V01．51．No．2 2008年3月 CHINESE GEOPHYSICS JOURNALOF Mat．，2008 朱守彪，邢会林，谢富仁等．地震发生过程的有限单元法模拟——以苏门答腊俯冲带上的大地震为例．地球物理学报，2008。 51(2)：460一468 ZhuS H F a1．Simulationof caseof the L，XieR，et finiteelementmethod：The on B．Xing eaahquakeprocessesby megathrustearthquakes Sumatra zone．Chinese subducfion Chinese)，2008．51(2)：460—468 J．Geophys．(in 地震发生过程的有限单元法模拟 以苏门答腊俯冲带上的大地震为例 朱守彪1’2”，邢会林2，谢富仁1，石耀霖3 1中国地震局地壳应力研究所，北京 100085 2Earth Science of Systems ComputationalCentre(ESSCC)，The 4072。Australia UniversityQueensland，St．Lucia．QLD 3中国科学院研究生院计算地球动力学重点实验室，北京 100049 摘要地震预报由目前的经验预报走向物理预报，数值模拟地震过程是其中的关键．文中应用统一的数学公式 时间积分方法，基于R最小策略，控制时间步长以保持力学状态变化稳定，从而保证有限元计算过程平稳、收敛．以 2004年发生过M。=9．3特大地震的苏门答腊俯冲带为例，模拟了俯冲带上俯冲板片与上伏板块之间的闭锁、解锁、 滑动到再闭锁这一准周期性过程，即地震的孕育、发生过程．计算结果表明，俯冲带上具有较大尺度、介质均匀、摩 擦系数相同的区域是产生大规模、大幅度整体突然滑动(即大地震)的条件；模拟的苏门答腊俯冲带上的大地震在 时间上有准周期性，空间上有迁移特征，破裂由深部向浅部进行；此外，俯冲带的几何特征对大地震的震源位置有 很大的影响． 关键词 有限单元，摩擦接触，地震，苏门答腊俯冲带 文章编号0001—5733(2008102—0460一09中图分类号P315 Simulationof finiteelement earthquakeprocessesby method： Thecaseof ontheSumatra megathrustearthquakes subductionzone ZHU Shou．Bia01·2，XINGHui．Lin2，XIEFu．Renl，SHIYao．Lin3 1 Institute of Crustal功wm讧，ChinaganhquakeAdministration，miji,皤100085，China 2