俯冲带深部应国场的二维弹性有限元数值模拟.pdfVIP

维普资讯 第 24卷 第 3期 地 震 学 报 ol24 Na 3 2002年 5月 (277～284) ACTA SEISM ()LOGICA SINICA M av．2002 文章编号 ：0253—3782(2002)03—0277—08 俯冲带深部应力场的二维弹性 有限元数值模拟 毛兴华 刘亚静 叶国扬 宁杰远 (中国北京 100871北京大学地球 物理 学系) 摘要 基于一些简化 的数值模型，根据弹性本构关系 ．用平面应变有限元方法计算 了相变 引 起的体积变化、板块 内部温度差、密度异常及边界力产生的应力场分布情况．数值模拟结果 显示 ，热应力能够解释俯冲带深源地震的应力场方 向特征 ，但解释不 了深源地震 的深度分布 特征 ；有亚稳态橄榄石存在时密度异常所产生 的应力场特征与地震观测结果所显示 的应 力场 特征有所偏离 ；虽然橄榄石 尖 晶石相变体积变化所产生的应力在橄榄石 尖晶石相变过渡 区附近有最大值 ，其数值远远超过温度差和密度异常产生 的应力场 的最大剪应力数值 ，但在 相变界面附近的区域 ，主压应力方 向垂直于相变界面的方 向，与地震观测结果所显示 的主压 应力的方 向不一致．所以，不能用弹性模拟得 出的俯冲带温度变化所产生 的热应力、俯冲带 密度变化所产生 的应力、相变体积变化所产生的应力来对深源地震进行简单化 的解释 ． 关键词 俯冲带 应力场 数值模拟 弹性 深源地震 中图分类号 ：P315．727 文献标识码 ：A 引言 俯冲带深源地震的形成机制是地球科学界十分关注的问题．在认识深源地震形成机制 方面的一个重要进展是 ，Sung和 Burns(1976a，b)基于相变动力学计算 ，提 出在低温 的俯 冲板块核心，橄榄石将 以亚稳态的形式存在于橄榄石一尖晶石稳态相变界面以下 的深度． 这一推测为 20世纪 90年代以来 的一系列相变动力学实验计算结果所证实(Rubie，Ross， 1994)．Iidaka和Suetsugu(1992)利用 日本 台网的观测资料 ，研究了发生在伊豆 小笠原俯 冲带震源深度为 3O0～500km地震的走时，发现当设伊豆一小笠原俯冲带在 4O0～500km 存在～5Yoo的低速异常时，能更好地拟合观测资料．宁杰远和臧绍先 (1999，2001)的数值 模拟结果也显示俯冲带有可能有这样 的亚稳态橄榄石楔存在． 人们进行 了大量实验室实验 (Rubie，1993)．尤其值得一提的是 ，人们根据实验室模拟 结果提 出了与相变相联系的断层 的概念模型 (Kirby，1987；Green，Burnley，1989；Kirbyet 国家科委 “九五”攀登计划预选项E1(95专一05)和高等学校骨干教师资助计划项 目(GG 170—10001—1022)共同资 助． 2002—01—15收到初稿 ，200205—22收到修改稿并决定采用． 维普资讯 地 震 学 报 a1．，199l；Burnleyeta1．，1991)．与地震观测资料所显示的深源地震以剪切破裂为主的图 象 (Kawakatsu，1991)相吻合 ，受到了大家的关注．但是 ，这种模型无法解释不存在亚稳态 橄榄楔区域 的俯冲带 的深源地震活动．而且 ，实验室模拟受时间和空间尺度的限制 ，基于 实验室模拟结果提出的概念模型需要用数值模拟 的方法来进行检验． 俯冲带应力场的数值模拟是认识俯冲带深源地震形成机制必不可少 的工作．Goto等 (1983，1987)对俯冲带的深部应力场进行 了二维粘弹性有 限元数值模拟，初步给出与深源 地震震源机制解所反映的应力场方 向相一致的俯冲带深部应力场计算结果．Bina(1997)和 Yoshioka等 (1997)分别基于橄榄石 尖晶石相变动力学计算结果，用二维弹性模型对俯 冲板块 的应力场进行了数值模拟 ，认为俯冲板块与周围地幔的密度差所产生 的应力能够解 释深源地震震源机制解所反 映的应力场方 向特征，也与地震应力 降 (Kikuchi，Fukao，