基于连续－离散单元法的汶川地震动力学二维自发破裂全周期模拟研究

首先利用地震地质调查、地震剖面探测结果对汶川地震的发震断层几何形态及周边块体介质参数进行约束,建立了3个包含有不同主断层的龙门山断裂带二维黏弹性有限元模型;然后基于连续-离散单元法对汶川地震的孕育与发生进行了全周期模拟,包括震间应力积累阶段、同震应力释放阶段及震后恢复调整阶段。模型对比研究结果表明,龙门山断裂带中由于断层几何结构及断层面剪切应力与正应力的综合作用,使得最先达到断层破裂准则阈值的灌县-江油断层率先破裂,初始破裂在深度15~20km处,破裂进一步传导至映秀-北川断层;同时由于灌县-江油断层与映秀-北川断层的破裂导致系统内应力大幅度降低,使得灌县-江油断层不具备产生同震破裂的应力条件。然而,尽管汶川-茂汶断层上的应力水平与汶川地震前相比有所减弱,但仍积累了较高的震间应力,可能预示其在汶川地震后具有较高的地震危险性。模拟结果还表明,铲型逆冲断层系统的深部断层破裂能够推动浅部高倾角区域发生被动破裂,因而其发震断层的倾角上限很可能比传统认识的大。