芦山地震完全库仑应力变化与应力场SEM数值模拟

摘要

地震库仑应力是研究地震应力触发的重要物理参数.关于地震应力触发的研究一般分为静态应力触发和动态应力触发,前者研究的是地震的断层位错产生的静态库伦应力变化对后续地震的触发作用,后者研究的是地震破裂扩展与地震波通过时产生的动态库伦应力变化对后续地震的触发作用.由于两者的触发性质不同:静态应力是通过持久性地改变施加在断层上的应力来触发地震,而动态应力是通过较强的瞬态应力冲击来触发后续地震,Kilb等提出了包含动态和静态库仑破裂应力的完全库仑破裂应力概念,并强调了动态库仑破裂应力的触发作用.rn 首先进行空间网格剖分;然后求出控制方程，暨根据所需求解的非分方程的弱解形式进行空间离散，再进行总体刚度矩阵合成。Maday等将拉格朗日插值函数引入谱元法中，并与Gauss-Lobatto-Legendre(GLL)积分相结合来合成对角质量矩阵，在时间域可采用显式的时间差分算法，有效提高计算效率。因此，将构造应力场的模拟考虑成粘弹性有限元问题，通过CMT反演获得芦山地震的余震矩张量解，然后根据应力主轴的空间约束，分别建立了空间四面体、六面体和八面体网格，最后利用GLL积分法求解粘弹性力学的控制方程。