吴建超　陈蜀俊　陈俊华　蔡永建

摘要：基于线弹性理论有限元方法，建立三峡库首区的三维有限元模型，模拟计算了库首区尤其是仙女山断裂北段在蓄水后的全位移场和构造应力场变化。结果表明：三峡水库175m水位蓄水，造成仙女山断裂和九畹溪断裂北段处于全位移场归州低值区和新滩高值区的梯度带上：同时，仙女山断裂和九畹溪断裂北段深部3组主应力值较蓄水前陡增。这样的地段全位移差值变化明显，3组主应力值增高，剪应力和剪应变能易于集中，容易加速仙女山断裂和九畹溪断裂北端失稳而诱发水库地震。

关键词：三峡水库；仙女山断裂；数值模拟；孕震机理；水库地震

中图分类号：P315.7 文献标识码：A 文章编号：1000-0666（2016）02-0218-06

0 引言

三峡水库自2003年6月蓄水以来，在库区周边发生了数千次地震，绝大多数为M≤2.0非构造成因的诱发地震（车用太等，2009）。然而，伴随着三峡水库175m试验性蓄水进程，2008年11月22日、2014年3月27日和3月30日，仙女山断裂北段附近分别发生了M_S4.1、M_S4.2和M_S4.5地震，引发了社会各界对于三峡水库是否诱发了地震的讨论。

笔者参加了2014年3月30日秭归M_S4.5地震的应急调查工作，该次地震震中烈度为Ⅴ度，有感范围较大。在郭家坝镇头道河村、郭家坝村等地出现Ⅵ度烈度异常区，且造成了一定程度的房屋破坏和财产损失。根据地震监测记录波形分析，认为此次地震为构造型水库地震，与仙女山微地块的活动有关。

值得注意的是，蓄水后三峡库区的地震活动多发生在仙女山断裂和九畹溪断裂的围限区附近，属于这两条断裂的北端部。……