自1931年希腊的马拉松水库4.7级地震后提出水库诱发地震问题后,20世纪60年代发生过几个具有较大影响的水库地震,分别为:1962年中国新丰江水库6.1级地震及1967年印度的柯依那水库6.5级地震,使得科学家们将水库地震问题作为一个科学问题进行较系统的研究,许多学者做了大量的探索性研究工作,使这一学科得到了迅速发展(Gupta,Rastogi,1976; Gupta,1985; 丁原章等,1989)。但水库诱发地震是一个复杂的问题,受众多条件因素影响,对水库诱发地震的孕震过程及其机理认识一直存在较大争议,主要原因是缺乏对水库地震过程及其有关现象的观测研究,对水库诱发地震机理,诱发和发展过程研究尚处于推测阶段(毛玉平等,2007)。不过随着数字地震观测技术诞生和应用,水库地震的研究也朝着纵深方向推进。Roeloffs(1988)指出引起地震活动的孔压增加与湖水位变化的频度与幅度有关,Simpson等(1988)将大型水库注水后的地震响应类型分为延迟响应和快速响应两种,Pradeep和Acree(1985)研究了这两种响应类型的地震活动特征的差异。

糯扎渡水电站位于云南省普洱市思茅区(左岸)与澜沧县(右岸)交界处的澜沧江南北向峡谷地段,设计最大坝高262.5 m,设计正常蓄水位高程812 m,水库总库容237.03×10⁸ m³,库长214 km,水库回水与上游大朝山水电站衔接,电站下游与景洪水电站相连,是澜沧江中下流梯级规划中工程规模、调节库容、装机容量和发电量最大的电站,于2011年11月19日开始蓄水。糯扎渡电站水库地震监测台网于2008年6月建成,并开始进行地震监测的工作。鉴于水库工程的重要性,研究水库地震事件的活动特征是很有必要的。本文利用糯扎渡水库台网的基础地震资料对糯扎渡水库蓄水前后的区域地震活动性特征进行对比研究,为研究水库诱发地震提供科学依据。

糯扎渡水电站工程场地位于耿马—澜沧地震带和思茅—普洱地震区2个二级地震带(区)之间,区域外围控制性的主干断裂有3组近南北向的澜沧江断裂,北西向的营盘山断裂、普文断裂、木戛—谦迈断裂、澜沧—勐遮断裂(耿马—澜沧断裂南段)、勐满—勐海断裂和北东向的南汀河断裂、澜沧—孟连断裂等。库区构造较发育,有近南北向的澜沧江断裂,北西向的帮东断裂、谦六断裂、平掌寨断裂、白马山断裂、酒房断裂、李子箐断裂、肖塘断层、麻栗坪断裂、肖塘箐断裂和北东向的窝拖寨断层(图1)。库区主库、各支库的部分库段分布在大范围的碳酸盐岩、块状花岗岩和块状石英闪长岩地区; 库区断裂多为逆断或局部逆断性质,大部分的断裂穿越主库或支库,库水淹没断层的现象较普遍 中国地震局地质研究所,糯扎渡水电站地震工程研究队.1993.澜沧江糯扎渡水电站地震安全性评价报告.。

根据区域地质构造资料 云南省地质局.1977.区域地质调查报告(思茅幅). 云南省地质局.1977.区域地质调查报告(景洪幅).,糯扎渡水电站西部属上芒播褶皱束,有近南北向的小青乃向斜、东西向的龙潭向斜,东部为近南北向的上芒播背斜。在水库的库首段的东侧出露海西晚期混合花岗岩,西侧出露下元古界角闪岩、黑云岩混合岩和花岗闪长质混合岩。库区出露上三叠统花岗岩、砾岩及岩屑砂岩; 在外围出露中侏罗统粉沙质泥岩、钙质泥岩。由于库区及周边复杂的地质构造,历史上就是地震多发区,如图1所示。李志祥等(2008)曾采用概率计算方法得出糯扎渡水库蓄水后将发生水库诱发地震,可能的最大震级为5.5级。

2012年6月1日至9月17日在云南省景谷县和双江县的交界处持续发生小地震,随后在2012年9月18日发生了景谷M4.2地震,将此时段的地震统称为景谷双江小地震群,该震群的震中分布见图5a。由于定位结果中的经纬度只保留到小数点后两位,故在图中呈现为整齐的分散的棋盘形式。

地震事件的精确定位是研究水库地震发生环境及机理的基础,而深度的精确测定是揭示特定属性的水库体系在库水的动态加载过程中库区下部介质和应力场发生响应的重要“云图”(陈翰林等,2009)。故选择Hypodd双差定位方法对该震群进行了重定位,为了保障定位结果的可靠性,所选择的数据为4个以上台站定位的地震事件共991个。重定位后的震中分布如图5b所示,与初始定位结果相比,震源位置在重定位后空间分布上更为集中,更突出其丛集性,而重定位前的结果则比较分散。重定位前后都可看出此震群主要分布在窝拖寨断裂上,由图1可看出窝拖寨位于糯扎渡水库库区范围内,据澜沧江糯扎渡水电站地震安全性评价报告 中国地震局地质研究所,糯扎渡水电站地震工程研究队.1993.澜沧江糯扎渡水电站地震安全性评价报告.中指出窝拖寨断裂现今构造活

图5 利用Hypodd双差定位方法重定位前(a)后(b)的景谷—双江小地震群震中分布图
Fig.5 The epicenter distribution map of the Jinggu-Shuangjiang small earthquake swarm before(a)and after(b)relocation by using the Hypodd double difference location method

动不明显,但沿断层的地震活动有密集成带的现象,并伴有中等强度的地震发生。在云南地区历史地震目录中,自云南台网建成以后,沿窝拖寨断裂的确有密集成带的地震活动发生,并伴有M4.0～5.0中等地震发生。图7为选取1986～2008年窝拖寨断裂附近(23°～23.5°N,99.8°～100.2°E)的地震事件的M-t图及月频度图,从图中可以看出1990年、1992～1994年均有密集的地震活动,且都有一个M4.0～5.0中等地震相随。所以此次小震群与之前所查阅的历史地震有着相似的特征,同样是沿着窝拖寨断裂成丛集出现,并伴有一个M4.2地震。但是如图7所示,此次震群的最大震级并未超过该地区的本底(历史地震目录中最大的地震为M5.3)。图7 b与图3b相比较,可以发现虽然图3b为整个研究区M_L≥2.0地震的频度,但是两幅图中所示M_L≥2.0地震发生的频度最大值相差不多,所以月频度也并未超过该地区的本底。

图6a为重定位前的地震震源深度,从图中可以看出深度主要分布在5～10 km,且大部分分布在10 km处,这是因为人工定位软件的设定为如果震源深度≤5 km,则自动跳至10 km,至2012年10月才进行了修改,深度不再自动给出,所以图中所示序号800以后的地震深度主要分布在0～5 km。由图6b可看出,重定位后的震源深度从0～20 km范围内都有分布,其中8～9 km居多,随时间无太大变化。

由区域地质构造资料 云南省地质局.1977.区域地质调查报告(双江幅).可知,窝拖寨断裂走向

图6 重定位前(a)、后(b)的深度分布图
Fig.6 The depth distribution before(a)and after(b)the relocation

图7 窝拖寨断裂及其附近历史地震M-t图(a)与月频度图(b)
Fig.7 M-t(a)and monthly frequency diagram(b)of historical earthquake on the Wotuozhai Fault and its surrounding area

北东,倾向南东,属于逆断层,由图5a可看出震群震中分布区域至断裂仅2～5 km,主要分布在断裂的南边。所以该震群与窝拖寨断裂有一定的关系。所以重定位后的震源深度,受蓄水影响,震源深度在0～10 km范围内的分布较多,但在10～20 km范围内也有分布,所以这部分地震的出现与窝拖寨断裂有关。

本文研究结果表明,糯扎渡水库蓄水前后的研究区的地震活动特征有明显差异,具体表现如下:

(1)云南省糯扎渡水库在蓄水后,库区地震大量增加,M_L≥2.0的地震也大幅增加。且地震频度及强度的变化和水库水位的变化有一定的对应关系,表明蓄水后库区有水库诱发地震发生。

(2)对蓄水后地震活动表现较强的景谷双江小震群采用Hypodd双差定位方法进行重定位,重定位后的结果相对于人工定位的结果在空间分布上更加集中,更具有震群的丛集分布特征。且重定位后的震源深度分布为0～20 km。其中,震源深度8～9 km分布最多,可能是受蓄水影响较大,结合窝拖寨断裂的产状特征,震源深度在10～20 km范围内的地震可能与断裂有关。通过查阅云南地区历史地震目录,自云南台网建成以后,沿窝拖寨断裂就有密集成带的地震活动发生,并伴有中等地震发生,此次震群也具有这些特征。此外小震群的最大震级和M_L≥2.0地震频度都未超过该地区的本底。综合以上分析认为,该震群的发生时间虽与水位相关,但可能是触发了构造地震,所以震源深度显现的也是构造活动的特征。