三峡水库蓄水后仙女山断裂北段数值模拟及其孕震机理的研究.docVIP

三峡水库蓄水后仙女山断裂北段数值模拟及其孕震机理的研究 　　摘要：基于线弹性理论有限元方法，建立三峡库首区的三维有限元模型，模拟计算了库首区尤其是仙女山断裂北段在蓄水后的全位移场和构造应力场变化。结果表明：三峡水库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日，仙女山断裂北段附近分别发生了MS4.1、MS4.2和MS4.5地震，引发了社会各界对于三峡水库是否诱发了地震的讨论。 　　笔者参加了2014年3月30日秭归MS4.5地震的应急调查工作，该次地震震中烈度为Ⅴ度，有感范围较大。在郭家坝镇头道河村、郭家坝村等地出现Ⅵ度烈度异常区，且造成了一定程度的房屋破坏和财产损失。根据地震监测记录波形分析，认为此次地震为构造型水库地震，与仙女山微地块的活动有关。 　　值得注意的是，蓄水后三峡库区的地震活动多发生在仙女山断裂和九畹溪断裂的围限区附近，属于这两条断裂的北端部。断层端部构造应力场及孕震机理一直是国内外研究的热点。断层的末端扩展受阻，往往会在断层的端部产生应力集中而引发地震（Obara et al，2009；孔礼健等，2009）。而三峡水库蓄水后的形变观测结果显示，仙女山断裂北段地块的倾斜场发生了改变，由蓄水前的正断层运动转为蓄水后的逆断运动性质（张燕等，2005；邢灿飞等，2006）。因此，综合考虑蓄水效应，研究蓄水后仙女山断裂北段构造应力场及其孕震机理具有重要意义。 　　本文应用三维有限元方法模拟计算了三峡水库蓄水对仙女山断裂北段构造应力场的影响，并对其孕震机理进行了分析和讨论，以期对三峡水库诱发地震机理作进一步探索。 　　1 地震活动概况 　　三峡水库自2003年6月蓄水以来，仙女山断裂北段附近共发生了41次MS≥2.0地震，其中MS4.0～4.9地震3次，MS3.0～3.9地震3次，MS2.0～2.9地震35次，且大部分位于长江南岸，呈密集面状分布特征，如图1所示，地震目录取自湖北省地震局台网中心。 　　车用太等（2009）认为，2008年11月22日胡家坪MS4.1地震是在水库水体荷载与库水下渗的共同作用下沿仙女山断裂发生的构造型水库诱发地震；吴建超等（2012）认为胡家坪MS4.1地震发生在等效应力高值区与低值区的梯度带上，为等效应力变化最为剧烈、最易产生应变能积累和应力集中的区域；陈俊华等（2013）从震群活动时序特征、震源体大小、震源机制及地震波形的频谱特征等方面进行研究，认为2012年10月31日屈原镇MS3.2地震为水库触发型构造地震；徐长朋（2012）对三峡库区119次小震的震源机制解进行研究，认为仙女山断裂北段地震震源机制以走滑为主。 　　上述研究多针对蓄水所诱发的地震，未对蓄水后仙女山断裂上的位移场和深部应力场变化进行综合研究。因此，有必要通过对蓄水后仙女山断裂北段的孕震动力学环境进行数值模拟，从上地壳深部探究其孕震机理。 　　2 研究区三维有限元模型 　　本文构建的三维有限元模型范围为110km×85km，深度取至上地壳底部埋深12km处。魏贵春等（2013）利用双差定位法对三峡数字遥测地震台网记录到的秭归附近地区73个MS≥1.0的地震进行了重新定位，结果显示，绝大多数地震事件的震源深度均在7km左右。这一深度刚好位于本文三维有限元模型的中部，计算过程中避免了受边界条件的影响，结果相对准确。因此，本文构建的三维有限元模型是合理与可靠的。图2为考虑地形地貌、断裂和构造层建立的三维模型，模型以（110°00′00″N，31°56′35″E）为坐标原点，X轴正方向为SE33°，Y轴正方向为NE33°，z轴垂直于水平面，向上为正。 　　模型按上地壳岩石层的物性参数可划分为沉积盖层和基底构造层。其中，沉积盖层厚约7km，主要为花岗岩、灰岩和砂岩；基底构造层层厚约5km，主要为花岗岩。依据陈蜀俊（2005）研究成果，断裂带的参数分别取为：弹性模量E=1.07×1010Pa，泊松比v=0.24，密度ρ=2570kg?m-3。有限元网格划分以