摘要：采用2008年6月至2013年10月糯扎渡水库地震台网记录的6300个地震，对水库蓄水前后的地震活动进行分析，发现库区蓄水后，地震频度与水位同步增长，地震活动显著增强，表明库区有水库诱发地震发生。蓄水前后的地震活动图像显示，澜沧江上景谷、双江交界处的小震活动增加最为显著，采用Hypodd双差相对定位方法对该震群进行了重新定位，定位后地震在空间上更为丛集，地震震源深度从0～20 km都有分布。该震群发生在窝拖寨断裂上，该断裂历史上曾发生过5.3级地震，因此分析认为，该震群的发生时间虽与水位相关，但可能是触发了构造地震，所以震源深度显现的是构造活动的特征。
关键词：糯扎渡水库；地震活动性特征；水库地震；Hypdd双差定位方法
中图分类号：P315.7文献标识码：A文章编号：1000-0666（2015）02-0189-07
0 引言
自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×lO8 m3，库长214 km，水库同水与上游人朝山水电站衔接，电站下游与景洪水电站相连，是澜沧江中下流梯级规划巾工程规模、调节库容、装机容量和发电量最人的电站，于2011年11月19日开始蓄水。糯扎渡电站水库地震监测台网于2008年6月建成，并开始进行地震监测的工作。鉴于水库工程的重要性，研究水库地震事件的活动特征是很有必要的。本文利用糯扎渡水库台网的基础地震资料对糯扎渡水库蓄水前后的区域地震活动性特征进行对比研究，为研究水库诱发地震提供科学依据。
1 糯扎渡水库地质构造环境
糯扎渡水电站工程场地位于耿马一澜沧地震带和思茅一普洱地震区2个二级地震带（区）之间，区域外围控制性的主下断裂有3组近南北向的澜沧江断裂，北西向的营盘山断裂、普文断裂、木戛一谦迈断裂、澜沧一勐遮断裂（耿马澜沧断裂南段）、勐满勐海断裂和北东向的南汀河断裂、澜沧一孟连断裂等。库区构造较发育，有近南北向的澜沧江断裂，北西向的帮东断裂、谦六断裂、平掌寨断裂、白马山断裂、酒房断裂、李子箐断裂、肖塘断层、麻栗坪断裂、肖塘箐断裂和北东向的窝拖寨断层（图1）。库区主库、各支库的部分库段分布在人范围的碳酸盐岩、块状花岗岩和块状石英闪长岩地区；库区断裂多为逆断或局部逆断性质，犬部分的断裂穿越主库或支库，库水淹没断层的现象较普遍。
根据区域地质构造资料，糯扎渡水电站西部属上芒播褶皱束，有近南北向的小青乃向斜、东西向的龙潭向斜，东部为近南北向的上芒播背斜。在水库的库首段的东侧出露海西晚期混合花岗岩，西侧出露下元冉界角闪岩、黑云岩混合岩和花岗闪长质混合岩。库区出露上三叠统花岗岩、砾岩及岩屑砂岩；在外围出露中侏罗统粉沙质泥岩、钙质泥岩。由于库区及周边复杂的地质构造，历史上就是地震多发区，如图1所示。李志祥等（2008）曾采用概率计算方法得出糯扎渡水库蓄水后将发生水库诱发地震，可能的最大震级为5.5级。
2 糯扎渡水库蓄水前后地震活动分析
2.1 资料选取
水库地震的选取范围，目前尚无统一标准，一般按照天然地震活动与诱发地震的空间分布特征来考虑。李永莉等（2012）选取的范围为库区20 km内，也有研究选取库区50 km范围，如毛玉平等（2008）和王绍晋等（2005）对云南漫湾水库和人朝山水库的研究；李家荣（1999）将在水域外20 km范围内引起的异常地震定义为水库诱发地震活动。欧作畿（2005）指出多数水库地震具有以下特点：震中多集中在库坝附近，一般集中分布在水库岸边几千米至十几千米范围内。震源深度较浅，一般在地表下10 km以内，个别达20 km，以4～7 km居多。陈章立小组的最新研究结果认为水库地震具有“双十”特征。本文选取从糯扎渡电站库首段到库尾段的澜沧江水域外10 km范围作为研究范围，如图2巾黑圈所示，时间段为2008年6月至2013年10月。同时考虑到定位的精度，只选取了使用4个以上台站定位的地震共6300个，其中糯扎渡水库蓄水前2008年6月至2011年11月18日773个MT≥-1.1地震，蓄水后2011年11月19日至2013年10月31日5 527个ML≥-1.1地震作为研究样本。
2.2 地震活动的空间分布
图2为研究区在蓄水前和蓄水后的地震震中分布情况。由图2a可看出水库蓄水前的三年多时间里的地震分布呈现为零散的自然分布，并且地震数目较少，ML≥2.0地震也很少，只有5个。但和图2b对比后可看出，蓄水后的近两年时问内地震明显增多，尤其是库坝区和库区巾段，地震活动增加明显，蓄水前研究区内地震有773个，而蓄水后有5527个。不仪弱震人量增加，而且ML≥2.0地震也有显著的增加，如图3b所示。 2.3 地震频度和强度与水位分析
根据2008年6月至2013年10月31日糯扎渡水库台网的地震资料，无论是蓄水前还是蓄水后，研究区的地震活动以弱震为主，特别是ML0.0～0.9级占绝对优势。蓄水前的观测资料表明，研究区地震频度的月平均值为18.6次，而蓄水后的月平均值为239.9次，是蓄水前的13倍。由图3可看出蓄水前的月频度普遍小于50，而蓄水后地震开始活跃。
图3为蓄水前后月频度与水位的关系图，糯扎渡水库设计蓄水位高程为812 m，从2011年11月19日开始蓄水后，水位一直都在上升，至2013年10月至最高水位811.92 m。糯扎渡水库从2011年11月19日至2012年7月，水库水位变化很大，从610.33 m上升至765 m左右，由图3可看出，蓄水初期2011年12月至2012年3月地震活动平静，与蓄水前无差别。随着水位的不断上升地震活动活跃度不断加强，地震活动频次不断增加，尤其到了2012年6月地震频次大幅上升，2012年7月最为活跃，月地震频次达到503。至2012年8月水库水位上升变缓，地震活动也相应地变弱，2012年9月至2013年2月水位变化都不大，但也在缓慢上升，由765 m上升至774 m，地震活动变化也比较活跃。2013年3～7月地震活动呈逐渐下降的趋势，相对应的水位值在774 m左右。2013年8月水位开始明显上升，至2013年10月水位由775 m上升至812 m，地震活动频次也明显活跃。由此可见，研究区的地震活动频次与水库水位变化具有一定相关性。
图4为ML≥2.0地震的M-t与水位的关系图。由图可看出蓄水前的MT≥2.0地震很少，只有5个，分布也无规律，最大为发生在库区中段的ML3.1地震。蓄水后发生的最人地震为景谷ML 4.5地震，发震位置在糯扎渡水库库区上游，这是糯扎渡水库地震台网建成后所发生的最火地震。蓄水后ML≥2.0地震有79个，由图2b可看出蓄水后的ML≥2.0地震人多分布于库坝区及库区上游。如图4所示，蓄水后ML≥2.0地震开始增多，2012年6月至2013年4月ML≥2.0地震尤为密集，此时间段又可分为3个阶段，即2012年6月至2012年8月上旬，2012年9月中旬至2013年1月上旬，2013年1月下旬至2013年4月。第一阶段糯扎渡水库水位由 725.5 m上升至767 m，第二阶段水位由767 m上升至775 m，第三阶段水位由775 m下降至771 m，又上升至775 m。此3个阶段的间隔时间里都没有ML≥ 2.0地震，水位都保持一个数值上下变化。
综合上述资料分析，糯扎渡水库蓄水后，库区的地震活动明显增强，频度和强度都有显著增加。且地震频度和强度于2012年6月都明显提高，虽然糯扎渡水库是从2011年11月19日就开始蓄水，但水位到达一定高度后、蓄水量达到一定量级后，地震活动与水位才开始相关。地震频度与强度的明显变化都伴随着水位的变化，虽然地震并不是都发生在水位变化剧烈的时期，但是水位都有变化，所以地震频度及强度的变化和水库水位有一定的对应关系。
3 景谷双江小震群
2012年6月1日至9月17日在云南省景谷县和双江县的交界处持续发生小地震，随后在2012年9月18日发生了景谷M4.2地震，将此时段的地震统称为景谷双江小地震群，该震群的震中分布见图5a。由于定位结果中的经纬度只保留到小数点后两位，故在图中呈现为整齐的分散的棋盘形式。
地震事件的精确定位是研究水库地震发生环境及机理的基础，而深度的精确测定是揭示特定属性的水库体系在库水的动态加载过程中库区下部介质和应力场发生响应的重要“云图”（陈翰林等，2009）。故选择Hypodd双差定位方法对该震群进行了重定位，为了保障定位结果的可靠性，所选择的数据为4个以上台站定位的地震事件共991个。重定位后的震中分布如图5b所示，与初始定位结果相比，震源位置在重定位后空间分布上更为集中，更突出其丛集性，而重定位前的结果则比较分散。重定位前后都可看出此震群主要分布在窝拖寨断裂上，由图1可看出窝拖寨位于糯扎渡水库库区范围内，据澜沧江糯扎渡水电站地震安全性评价报告中指出窝拖寨断裂现今构造活动不明显，但沿断层的地震活动有密集成带的现象，并伴有中等强度的地震发生。在云南地区历史地震目录中，自云南台网建成以后，沿窝拖寨断裂的确有密集成带的地震活动发生，并伴有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）。图7b与图3b相比较，可以发现虽然图3b为整个研究区ML≥2.0地震的频度，但是两幅图中所示ML≥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居多，随时问无太火变化。
由区域地质构造资料可知，窝拖寨断裂走向北东，倾向南东，属于逆断层，由图5a可看出震群震中分布区域至断裂仪2～5 km，主要分布在断裂的南边。所以该震群与窝拖寨断裂有一定的关系。所以重定位后的震源深度，受蓄水影响，震源深度在0～10 km范围内的分布较多，但在10～20 km范围内也有分布，所以这部分地震的出现与窝拖寨断裂有关。
4 结论与讨论
本文研究结果表明，糯扎渡水库蓄水前后的研究区的地震活动特征有明显差异，具体表现如下：
（1）云南省糯扎渡水库在蓄水后，库区地震大量增加，ML≥2.0的地震也大幅增加。且地震频度及强度的变化和水库水位的变化有一定的对应关系，表明蓄水后库区有水库诱发地震发生。
（2）对蓄水后地震活动表现较强的景谷双江小震群采用Hypodd双差定位方法进行重定位，重定位后的结果相对于人工定位的结果在空间分布上更加集中，更具有震群的丛集分布特征。且重定位后的震源深度分布为0～20 km。其中，震源深度8～9 km分布最多，可能是受蓄水影响较大，结合窝拖寨断裂的产状特征，震源深度在10～20 km范围内的地震可能与断裂有关。通过查阅云南地区历史地震目录，自云南台网建成以后，沿窝拖寨断裂就有密集成带的地震活动发生，并伴有中等地震发生，此次震群也具有这些特征。此外小震群的最大震级和ML≥2.0地震频度都未超过该地区的本底。综合以上分析认为，该震群的发生时间虽与水位相关，但可能是触发了构造地震，所以震源深度显现的也是构造活动的特征。
