摘要
了解地震发生的动力学机制是研究地震发震原因的关键，而数值模拟的方法是高速、有效的手段。2011年3月11日日本东北部宫城县发生9.0级大地震，文中以该次大地震所在的日本俯冲带为研究对象，通过使用黏弹性有限元数值模拟，并引用接触对，建立了研究区二维数值模型，模拟俯冲带与上覆板片之间的滑动、黏滞到再滑动的过程，亦即断层失稳发生地震的过程。模拟结果显示，随着太平洋板块不断俯冲，在俯冲带上自发出现了断层闭锁、解锁到再闭锁的黏滑过程，且这种过程呈现一定的准周期性，大事件主要集中分布在20～30 km的深度范围内。根据俯冲带可能在俯冲过程中角度的变化，建立了不同的模型，进行模拟对比研究，结果表明，俯冲带的几何形态，以及俯冲角度变化所在的不同深度，对模拟的结果有不同的影响。