2012, Vol. 55
2. 广西地震局龙滩台网中心,广西 天峨 547300;
3. 广西天峨县地震局,广西 天峨 547300
2. Longtan Network Center, Guangxi Earthquake Administration, Tian'e, Guangxi 547300,China;
3. Gaungxi Tian'e Earthquake Administration, Tian'e, Guangxi 547300,China
近几十年来，随着国内外大型水库的兴建，出现了很多水库蓄水诱发地震的例子.其中1962年3月20日我国新丰江水库地震达到M6.2 级，1967 年12月10日印度的Koyna地震达到M6.3级[1-6].较大水库诱发地震的发生会破坏水库大坝本身结构，从而威胁着库区附近人民生命财产安全和水库工程安全，需要加强对水库诱发地震的研究.
水库诱发地震的发生与库区的先存地质构造、断裂发育、岩性分布、水文地质构造、库水水位变化等密切有关[7-14].水库诱发地震发生的物理效应基于库伦定理，当岩石空隙压力增大或者摩擦系数、内聚力减小时，就会使得岩石强度减小，导致岩石失稳、破裂，从而引发地震.当库区蓄水后，水库蓄水水位迅速增加，会导致库区正应力和孔隙水压力瞬时增加.库水水位增加，库容量增大，水体自重增加，库水会顺着库区破碎带向库区深部渗透，使得库区深部岩石摩擦系数减小，特别是在碳酸盐岩石中这种现象更明显，这是由于碳酸盐岩在水和二氧化碳的作用下能发生化学反应后发生化学侵蚀，生成可溶的碳酸氢钙.当库水顺着先存断裂破碎带扩散到地下一定深度，水的渗入使得碳酸盐岩发生化学侵蚀，使得岩石内聚力减小或者岩石的摩擦系数减小，致使岩石失稳、破裂，从而产生地震[3, 15-19].
当前地球物理探测仪器和技术发展很快，采用先进的地球物理观测技术探明库区先存断裂带体系、深部构造环境，在水库加卸载过程中观测库水对深部介质属性的变化，进而研究库区诱发地震的可能分布区域、水库诱发地震类型等，是当前对水库诱发地震研究的重要发展方向之一.1998 年在印度6.3级Koyna地震区进行过地震探测工作，发现在地震区5~10km 左右为低速层，大量库区诱发地震多分布在该深度区域，说明深部低速层的存在与地震分布相关[20-22].库水加载会影响库水在地下介质里的渗透变化，水的渗透对地下介质的电导率属性有极大的影响，大地电磁方法对深部介质电导率变化反映最灵敏、分辨力最高[23-24]，广泛应用于地震区深部孕震背景、断裂带内部结构探测研究[25-32].
广西龙滩库区是目前中国国内在建的仅次于长江三峡电站及金沙江溪洛渡电站的特大型水电工程.龙滩水库自2006年10月开始蓄水，到2008 年蓄水水位达到最高蓄水深度375 米.库区蓄水前仅有少量的地震活动，随着水库蓄水，库区地震频繁发生，至今龙滩遥测地震台网记录到上千次地震，已发生MS4.0级以上地震3次[33-34].2006 年至2007 年度发生的水库诱发地震精定位结果显示，龙滩地震主要分布在水库蓄水后先淹没的5 个深水区，明显成丛于5个区域[35].
鉴于广西龙滩水库工程的安全，在科技部支撑项目的支持下，利用电磁方法在广西龙滩库区开展库水加卸载动态过程电磁观测，获得库区深部地下电阻率随库水加卸载及渗透过程的变化特征，研究其与水库加卸载及渗透过程之间的关系.本文将介绍2009和2010年2个年度在龙滩库区库水加卸载过程中4个测点5 期次电磁定点观测结果，结合库区区域深部电性结构探测成果，对库区地下介质电性与库水加卸载进程中库水渗透的关系、水库地震孕震环境和地震类型等进行分析研究.
2 龙滩库区地质构造和电磁定点、电磁探测剖面位置广西龙滩水库横跨广西天峨县和贵州罗甸县，地处云贵高原南缘向桂西北山区与丘陵过渡的斜坡地带.库区主要由南盘江、蒙江、曹渡河、布柳河和红水河上游组成.区内主要发育4 条北东向和2 条北北东向断裂，从西往东依次为罗甸-罗苏断裂(F1)、八茂-百塘断裂(F2)、拉浪-达良断裂(F3)、龙凤-八腊断裂(F4)、天峨-金谷断裂(F5)、拉色-辉马断裂(F6)[36].研究区主要出露的地层有泥盆系(D)、石碳系(C)和二迭系(P)，三迭系(T).泥盆系地层为粉砂岩、石英砂岩、页岩碎屑岩和碳酸盐岩建造特征.石碳系(C)主要是灰岩等碳酸盐岩建造为主.二迭系分上、下统，下统(P1)为浅灰、灰色中厚状灰岩夹白云质灰岩等碳酸盐岩建造为主，上统(P2)为硅质灰岩、硅质页岩夹砂岩、页岩等碎屑岩和碳酸盐岩建造为特征.三迭系主要为砂岩、页岩、泥质灰岩等碎屑建造和碳酸盐岩建造.根据各地层的主要岩性性质，可知研究区内的泥盆系(D)、二迭系上统(P2)和三迭系(T)地层含水性差，石炭系(C)和二迭系下统(P1)地层易含水[36-37].
图 1中也放置了2006年至2010年发生的水库诱发地震精定位结果[35]，结果显示龙滩水库地震主要分布在水库蓄水后先淹没的5 个深水区，明显成丛于5个区域，即罗妥(丛Ⅰ)、八茂(丛Ⅱ)、坝首(丛Ⅲ)、拉浪(丛Ⅳ)和布柳河(丛Ⅴ).其中罗妥(丛Ⅰ)和八茂(丛Ⅱ)两个地震丛集区沿八茂-百塘断裂(F2)分布，并分别位于南盘江和蒙江水域；坝首(丛Ⅲ)、拉浪(丛Ⅳ)两个地震丛集区位于拉浪-达良断裂(F3)南北两端，分别分布在红水河和曹渡河水域，布柳河(丛Ⅴ)地震丛分布在布柳河水域.
|
图 1 广西龙滩库区地形（上）、地质构造（下）和电磁测量点位置图 (1):南盘江；2):蒙江；（3):曹渡河；（4):布柳河；（5):红水河；F1:罗甸一罗苏断裂；F2:八茂一百塘断裂；F3 :拉浪一达良断裂；F4 :龙凤一八腊断裂；F5 :天峨一金谷断裂；F6 :拉色一辉马断裂. Fig. 1 Map showing topography (upper) , tectonics (lower) and magnetotelluric sites in the Longtan reservoir (1) Nanpan River ； (2) MengRiver; (3) Caodu River ； (4) Buliu River ； (5) Hongshu River ； F1 ： Luodian-Luosu Fault ； F2 ： Bamao-Baitang Fault ；F3 ： Lalang-Daliang Fault ； F4 ： Longfeng-Bala Fault ； F5 ： Tian；e-Jinggu Fault ； F6 ： Lase-Huima Fault. |
两条大地电磁探测剖面基本垂直构造走向为北西-南东向布置，两条电磁剖面基本平行.其中A剖面西北起于贵州境内的立亭，向东南经过八茂、东南止于广西的龙腊，穿过八茂(丛Ⅱ)和拉浪(丛Ⅳ)地震丛集区.B 剖面西北起于贵州罗甸县境内的罗苏以西，向东南经过罗妥、拉当，东南止于广西的板别，穿过罗妥(丛Ⅰ)和坝首(丛Ⅲ)地震丛集区.
电磁定点位置位于库区中上游的坪上(定点1)，库区中部的中良坪测点(定点2)、库首区域的交必屯(定点3)和向阳(定点4)4 个地点作为电磁观测定点.坪上(定点1)位于电磁探测A 剖面的A07点上，交必屯(定点3)位于电磁探测剖面的B08 号点上.坪上(点1)和中良坪(点2)测点位于二迭(P)和三迭系(T)地层交接位置，向阳(点4)位于三迭系上统(T2)地层上，交比屯(点3)测点位于二迭系上统(P2)地层上(图 1).
3 库区库水加卸载进程中定点电阻率观测结果和分析 3.1 5期次电磁定点观测时段图 2给出了2006年10月1日到2010年10月31日的水库库水水位(上)和2.0 级以上地震时序图(下).图中可见龙滩水库未蓄水前地震较少，随着水库蓄水后水位上升，库区地震频繁发生.地震多发时段出现在库水高水位的2007年7-8月、2008年6-12月和2010年10月期间.
|
图 2 龙滩库区2006年10月1日到2010年10月31日库水水位变化(上)､ML>2.0地震时序图(下)和5期次电磁测量时段(上图中粗黑线) Fig. 2 Water level change curve (upper)，and time series of earthquakes (ML〉2. 0) (lower) from October 1, 2006 to October 31, 2010 in the Longtan reservoir (In the upper figure, vertical dashed lines denote 5 times of magnetotelluric measurements) |
为了在库水位最高、最低水位时段实施电磁定点测量，以获得每一年水库加卸载过程中水位变化最大时段的电磁观测数据，我们分析了2006 年10月1日水库蓄水以来的水位变化规律，并密切跟踪2009年到2010年11月期间研究区的降雨、库区水位变化趋势，最终在水位最佳时段进行了5 期次电磁观测工作.5期次电磁测量具体时段(图 2中黑色粗黑线)分别为(1):2009年4月23日到5月9日，水位约342m，(2):2009年9月4日到9月18日，水位约359m，(3):2009年12月27日到2010年1月6日，水位约350 m，(4):2010 年5 月27 日到5月6日，水位约331 m；(5)2010 年10 月1 日到10月16日，水位约372m(见图 2).
3.2 野外观测、处理和结果分析5期次电磁定点野外测量使用了3套加拿大凤凰公司生产的MTU-5A 大地电磁仪器.野外观测5分量，测量相互正交的2个的电场分量(Hx和Hy)和3个磁场分量(Hx、Hy和Hz).由于龙滩库区山势险峻，又是水库库水新淹没区，野外观测场地极难选择，4 个电磁观测点电极布置只能依照实际测量场地布置.具体观测参数见表 1.
|
表 1 4个电磁观测点的布置参数 Table 1 Layout parameters of 4 MT sites |
为了得到高信噪比的资料，每个测点每次测量时野外记录时间都在3 天以上.对所有资料全部采用Robust[38]的处理方法进行了处理.Robust数据处理方法根据观测误差的剩余功率谱的大小对数据进行加权，注重未被干扰的数据，降低突变点数据的权，使它对大地电磁阻抗估算值的影响降低到最小，从而改善受电磁噪声干扰的大地电磁测深资料.
对处理得到的谱数据利用大地电磁数据管理系统“MTDATAMNG"[39]计算得到了各测点5 期次的视电阻率和相位数据.图 3给出了交比屯(JB)(点3)、坪上(PS)(点1)、中良坪(ZL)(点2)、向阳(XY)(点4)测点5期次观测的视电阻率和阻抗相位曲线图，5期次测量结果分别用1、2、3、4、5数字来表示.全部4个测点测量得到的阻抗相位曲线在5期次测量过程中变化不大，且在低频段的观测结果信噪比小，在本文中只重点分析视电阻率测量结果.另外交比屯测点实地测量环境是一南北方向长、东西方向窄的紧靠水库区域的地点，南北方向的电极距长度78m，东西方向的电极距长度较短仅为15 m，重点关注南北方向的视电阻率测量结果.
|
图 3 交比屯(JB)、坪上（PS)、中良坪(ZL)和向阳（XY)测点5期次观测的视电阻率和阻抗相位曲线 Fig. 3 The apparent resistivity and impedance phase curves of five times of MT measurements in two directions at sites of Jiaobitun(JB)，Pingshang(PS)，Zhongliangping(ZL)and Xiangyang(XY) |
(1) 视电阻率曲线形态反应出地下深部电性结构特征:图 3中明显可见交比屯(点3)测量得到的视电阻率曲线形态、数值和5 期次测量结果的变化都与坪上(点1)、中良坪(点2)和向阳(点4)测点的明显不同.交必屯测点测量得到南北方向视电阻率曲线数值自高频到低频均在100Ωm以上，高频首枝到频率1 Hz左右视电阻率值约为200 Ωm，1Hz以下视电阻率值减小，到0.1 Hz 左右达到100Ωm，随后随频率再减小视电阻率数值增大.而坪上(点1)、中良坪(点2)和向阳(点4)测点2个方向的视电阻率曲线形态值基本一致，从高频到低频表现为高、低、低的特征.坪上和中良坪测点上的视电阻率曲线数值更一致，南北方向的视电阻率曲线首枝数值为200Ωm 左右，随着频率减小视电阻率曲线数值减小，到0.01 Hz附近减小到约1Ωm；东西方向的视电阻率曲线首枝数值略大于南北方向，接近1000Ωm，随着频率减小视电阻率数值也减小，到0.01Hz附近为几欧姆米，到更低频段，视电阻率值又为增大的趋势.向阳测点南北方向视电阻率曲线首枝数值为10Ωm 左右，随着频率减小视电阻率数值减小，到0.01Hz附近减小到1Ωm 以下，东西方向视电阻率曲线首枝数值略大于南北方向，接近100Ωm，随着频率减小视电阻率数值也减小，到0.01Hz附近约为1Ωm.说明交比屯测点下的深部结构在一定深度范围都表现为高电阻特性，在坪上(点1)、中良坪(点2)和向阳(点4)测点下的深部电性结构为高、低、高三层结构.
(2) 首枝视电阻率反应了出露地层的电性差异:4个测点的视电阻率曲线首枝视电阻率值不尽相同，交比屯测点上最大，坪上和中良坪测点次之，向阳测点上最小.反映出4 个测点上出露地层的电性差异，交比屯测点位于二迭系上统(P1)上，坪上测点靠近二迭系下统(P2)地层，中良坪点位于二迭和三迭系地层交接位置，向阳点则位于三迭系下统(T2)地层出露的腹地.二迭系上统(P2)地层主要是硅质灰岩、硅质页岩夹砂岩等碎屑岩为主，三迭系上统(T2)以砂页岩和泥质灰岩等为主，说明二迭系上统(P1)地层的电阻率数值较三迭系下统(T2)地层的稍大.
(3) 5期次观测的视电阻率数值变化趋势:为了对比4个测点上5期次观测的视电阻率数值与库水水位变化的关系，绘制了4个测点5期次观测的11个频率(80 Hz、10 Hz、3 Hz、1.13 Hz、0.56 Hz、0.28Hz、0.14Hz、0.07Hz、0.035Hz、0.0176Hz、0.0088Hz)2个方向的视电阻率值和水位、时间变化图(见图 4).图 4中各测点2个方向上80Hz、10Hz、3Hz3个频率的视电阻率值用黑色实线绘制，其余7个频率的视电阻率值用红色实线绘制.
|
图 4 交比屯(JB)、坪上（PS)、中良坪(ZL)和向阳（XY)测点5期次观测的11个频率的视电阻率值随水位变化图 (1)水位 342 m;(2)水位 359 m;(3)水位 350 m;(4)水位 331m;(5)水位 372 m. Fig. 4 The apparent resistivity values of 11 frequency of five times of MT measurements at sites of Jiaobitun(JB), Pingshang(PS) ,Zhongliangping(ZL) and Xiangyang(XY) changing with water level Frequency:1:80 Hz;2:10 Hz;3:3 Hz;4:1.13 Hz;5:0.56 Hz;6:0.28 Hz;7:0.14 Hz;8:0.07 Hz;9:0.035 Hz;A:0.0176 Hz;B:0.0088 Hz. (1) Water level 342 m;(2) Water level 359 m;(3) Water level 350 m;(4) Water level 331 m;(5) Water level 372 m. |
图 4中可清楚看出4个测点5次测量的视电阻率曲线数值在频率1Hz以上的高频都基本一致.交比屯(JB)测点在0.1Hz到较低频段，第1、2、3、和4次测量结果一致，水位最高时的第5 次测量的视电阻率数值小于其他4 次测量结果.坪上(PS)、中良坪(ZL)和向阳(XY)测点从频率1Hz开始到0.001Hz的视电阻率数值有明显变化，而且变化趋势一致，即2009年9 月(2)和2010 年10 月(5)水位最高时段也是水库库容最大、水体自重最大时段测量的视电阻率数值最低，2009年12月水位下降阶段(3)测量的视电阻率数值最高，在2009年5月(1)和2010年5月(4)最低水位时段测量的视电阻率数值居中.
4 结论和讨论 4.1 库水加卸载及渗透过程中的深部电阻率变化2009和2010年度龙滩库区库水加卸载过程中在交比屯、坪上、中良坪和向阳4个测点5期次电磁观测结果显示，在交比屯测点测量得到的视电阻率曲线和数值以及5期次观测结果视电阻率数值变化趋势与坪上、中良坪和向阳测点上测量结果明显不同.4个测点上5期次观测到视电阻率值变化特征，可以从龙滩库区区域深部电性结构得到进一步说明.
图 5为2009年实施的跨过龙滩库区的两条大地电磁测深剖面(A 和B剖面)自地表到深度20km的二维电性结构图像[40].A 剖面上的A07号测点为坪上电磁监测点，B 剖面上B08 测点为交比屯电磁监测点.根据地质构造图上断裂带的地表位置[36]，在图 5中标识出区内6条断裂带位置以及跨过库区河流的位置.同时也放置了2006-2010年5月龙滩库区地震精定位结果[35].在A 剖面放置的地震个数是垂直A 剖面两侧约0.08°范围的地震个数，在B剖面放置的地震是沿B 剖面西南侧0.1°范围的地震个数.在图 5中还放置了2010年9月18日发生的MS4.3级地震震源位置.
|
图 5 龙滩库区A剖面（上）和B剖面（下）二维电性结构和地震精定位分布图 F1:罗甸一望谟断裂；F2 :里落一白塘断裂；F3 :拉浪一达良断裂；F4 :龙凤一八腊断裂；F5 :天峨一金谷断裂；F6 :拉色一辉马断裂. Fig. 5 2-D geo-electrical structure model obtained by the NLCG(b) along the A and B profiles and the distribution of relocated earthquakes in the Longtan reservoir. F1: Luodian-Wangmu Fault； F2: Liluo-Baitang Fault； F3 : Lalang-Daliang Fault ；F4: Longfeng-Bala Fault； F5 : Tian；e-JingguFault； F6 : Lase-Huima Fault. |
图 5中清楚看出交比屯测点位于一上宽下窄似“铆钉"状的高阻体上方，高阻体内部电性结构较完整，向下延伸到约10km.根据区域地质资料[36]，推测该高阻体对应于由硅质灰岩、硅质页岩夹砂岩等组成的二迭系上统(P2)地层，其岩性不易含有地下水也不易受库水水位的影响，因而在该区域的交比屯测点仅在最高水位时段测量得视电阻率值在频率0.1Hz到较低频段，第1、2、3、和4 次测量结果一致，水位最高时的第5 次测量的视电阻率数值略小于其他4次测量结果.坪上测点处于自地表到地下20km 左右为高-低-高的近似3层的电性结构区域，在地表到深度约1km 是上千欧姆米的高电阻层，其下电阻率逐渐减小，低电阻层可以延伸到地下7~9km 左右，在深度约5km 深度附近电阻率值小，在9km 以下电阻率值又增大.根据地层岩性等资料[36-37]，推测坪上测点下浅表到深度1km 左右的高电阻对应三迭系下统(T1)和二迭系上统(P2)地层，其下低阻层对应于由易富水的灰岩组成的二迭系下统(P1)和石炭系(C)地层，再往深部的高阻层对应于含水性差的泥盆纪(D)地层.由于深度1km以下是灰岩组成的二迭系下统(P1)和石炭系(C)地层，易含地下水而易受水位影响，使得在坪上测点5期次观测得到的视电阻率数值在频率1 Hz左右到低频段发生变化，表现出在高水位时段测量的视电阻率数值小、在低水位时段测量的视电阻率数值稍大.中良坪和向阳测点上观测的5 期次视电阻率值也有相似变化特征.
深部电性结构和龙滩库区4个地点在库水加卸载及渗透过程中5期次测量的视电阻率数值和变化趋势说明龙滩库区库水的渗透作用可影响到地下深度7km 左右范围，不同岩性和不同深部电性结构的区域表现明显不同.当水库水位增加、水库库容面积增大，使得水体的自重加大，库水渗透过程受水体自重应力主要控制[41]，从而影响地下水在地下介质中的渗透和连通作用，从而是地下介质的电性发生变化.
4.2 龙滩库区深部孕震结构和地震类型图 1中显示库区地震丛Ⅰ(罗妥)、丛Ⅱ(八茂)、丛Ⅲ(坝首)、丛Ⅳ(拉浪)4 个地震丛的震源区垂直投影到地表都位于库区蓄水后水域淹没区的河流拐弯或者汇聚处，亦位于二迭和三迭系地层衔接位置，地震丛Ⅰ、丛Ⅱ 和丛Ⅳ 分别位于里落-白塘断裂(F2)和拉浪-达良断裂(F3)附近，地震丛Ⅲ位于龙凤-八腊断裂F4的西北侧.深部电性结构显示这4个地震丛集区的震源位置位于高、低电阻边界，其中地震丛Ⅰ、丛Ⅱ和丛Ⅳ 的震源位于低阻到高阻过渡边界.基于龙滩库区4个地震丛集区附近的水域、断裂、岩性、深部电性结构特征、地震丛集区震源结构特征、库水加卸载及渗透过程的深部电阻率观测结果、库区速度研究[42-43]、库水渗透及水库诱发地震机制研究[3, 17-18]等，对龙滩库区4 个地震丛集区地震类型和成因分析如下:
罗妥地震丛(丛Ⅰ)、八茂地震丛(丛Ⅱ)和拉浪地震丛(丛Ⅳ):3 个地震丛集区的震源区深部电性特征相似，都位于从浅部低阻层底部(约4km)到深部高阻层的过渡边界附近，位于由灰岩组成的具有低阻特性的石炭系地层的下部.低阻层对应于低VP、QP 和QS 值区[42-43].低阻层的形成推测是由于水库蓄水后库水随断裂破碎带渗透到地下易含水的主要是灰岩组成的二迭系下统(P2)和石炭系(C)地层中，易形成岩溶水体[3].岩溶水体易受到水库水位和水体自重变化的影响，导致在拉浪地震丛(丛Ⅳ)上坪上测点5 期次观测到的视电阻率数值在频率1Hz左右以下发生变化.岩溶含水体的渗透性表现为非均一性和各向异性，易引起孔隙压力在特定方向和特定部位的增高，从而诱发地震.推测这3个地震丛集区发生的地震主要是以岩溶诱发地震为主.岩溶诱发地震一般不大于4 级，在八茂地震丛(丛Ⅱ)和拉浪地震丛(丛Ⅳ)发生的地震至今没有大于4.0级地震.在罗妥地震丛(丛Ⅰ)于2007年3月17日发生了MS4.0 地震和2010 年9 月18 日发生了MS4.4地震，根据地震震级和数量统计，在罗妥地震丛(丛Ⅰ)发生的2~4级地震数远小于其他3个地震丛，另外从图 5中可见，罗妥地震丛(丛Ⅰ)上伏的低阻体的范围大于八茂地震丛(丛Ⅱ)和拉浪地震丛(丛Ⅳ)，推测在罗妥地震丛(丛Ⅰ)上伏的岩溶水体范围大，再加上2~4 级地震少，从而在该地震丛区发生了2个较大地震.
坝首地震丛(丛Ⅲ):处于库区水域最宽、蓄水水位最高的坝首区域，其震源区位于具有一定规模的上宽下窄似“铆钉状"高阻体下部，该高阻体具有一定深度(大约10km 左右)和宽度，内部电性结构较完整，不易含有地下水和不易受水库水位变化影响，因而在该区域交比屯测点上5期次电磁定点观测的视电阻率数值仅在水位最高时段发生变化.在坝首地震丛(丛Ⅲ)于2007 年7 月17 日发生了MS4.5地震，在2006-2009年期间3~4级地震有3个.推测该地震丛集区发生的地震是由于水库蓄水前，在构造应力场作用下高电阻脆性岩体内聚集了高变形能，在水库蓄水后水体压力增大和库水渗透作用下引起了脆性断裂引发的地震.
2009年和2010年两个年度在广西龙滩库区库水加卸载过程中的4个地点进行的5期次电磁观测以及库区区域深部电性结构探测研究，揭示了龙滩水库区区域以及4个地震丛集区的震源区深部电性结构特征和深部孕震环境，观测到库水加卸载及渗透过程中深部介质电性的变化，说明库水在加卸载及渗透过程中对地下介质有明显影响，库区上游库水渗透影响深度可达7km 左右，在坝首区域可达10km.电磁探测结果说明库水渗透对龙滩水库区地震孕育和发生起着重要、直接的作用.
致谢野外数据采集得到江汉石油管理局物探公司711电法队大力协助，赵翠萍研究员提供了龙滩库区地震精定位数据，蒋海昆、单新建研究员在本专题实施过程中给以了具体指导，广西地震局、广西天峨县地震局给予了大力支持，在此一并感谢.
| [1] | Chander R. Can dams and reservoirs cause earthquakes? Triggering of earthquakes. Resonance , 1999, 11: 1-10. |
| [2] | Gupta H K. A review of recent studies of triggered earthquakes by artificial water reservoirs with special emphasis on earthquakes in Koyna, India. Earth Science Reviews , 2002, 58(3-4): 279-310. DOI:10.1016/S0012-8252(02)00063-6 |
| [3] | Chen L Y, Talwan P. Reservoir-induced Seismicity in China. Pure Aappl. Geophys. , 1998, 153(1): 133-149. DOI:10.1007/s000240050188 |
| [4] | Simpson D W, Negmatullaev S K. Induced seismicity at Nurek reservoir, Tadjikistan, USSR. Bulletin of the Seismological Society of America , 1981, 71(5): 1561-1586. |
| [5] | 胡毓良. 诱发地震. 国际地震动态 , 1994(6): 37–38. Hu Y L. The reservoir induced earthquakes. Recent Developments in Word Seismology (in Chinese) , 1994(6): 37-38. |
| [6] | 易立新, 车用太, 王广才. 水库诱发地震研究的历史、现状与发展趋势. 华南地震 , 2003, 23(1): 28–37. Yi L X, Che Y T, Wang G C. Retrospection and prospect of the research on reservoir induced seismicity. Southern China Earthquake (in Chinese) , 2003, 23(1): 28-37. |
| [7] | Rajendran K. Sensitivity of a seismically active reservoir to low-amplitude fluctuations: observations from Lake Jocassee, South Carolina. Pure Appl. Geophys. , 1995, 145(1): 87-95. DOI:10.1007/BF00879485 |
| [8] | Hu P, Hu Y L. Advances in reservoir-induced seismicity research in China. Tecronophysics , 1992, 209(1-4): 331-337. DOI:10.1016/0040-1951(92)90043-6 |
| [9] | Hu Y L, Liu Z Y, Yang Q Y, et al. Induced seismicity at Wujiangdu Reservoir, China: A case induced in the karst area. Pure Appl. Geophys. , 1996, 147(2): 400-418. |
| [10] | Leith W, Simpson D W, Alvarez W. Structure and permeability: Geologic controls on induced seismicity at Nurek reservoir, Tadjikistan, USSR. Geology , 1981, 9(10): 440-444. DOI:10.1130/0091-7613(1981)9<440:SAPGCO>2.0.CO;2 |
| [11] | Talwani P. On the nature of reservoir-induced seismicity. Pure Appl. Geophys. , 1997, 150(3-4): 473-492. DOI:10.1007/s000240050089 |
| [12] | 杨主恩, 林传勇, 高振寰, 等. 东江水库诱发地震的地震地质背景研究. 地震地质 , 1995, 17(3): 241–252. Yang Z E, Ling C Y, Gao Z H, et al. Studies on seismogeological background for Dongjiang reservoir-Induced earthquake. Seismology and Geology (in Chinese) , 1995, 17(3): 241-252. |
| [13] | 林松建, 黄昭, 连玉平, 等. 水口水库诱发地震特征与发震断裂的探讨. 大地测量与地球动力学 , 2007, 27(1): 40–44. Lin S J, Huang Z, Lian Y P, et al. On characteristics and seismogenic faults of induced earthquake of ShuiKou reservoir in Fujian Province. Journal of Geodesy and Geodynamics (in Chinese) , 2007, 27(1): 40-44. |
| [14] | 陈蜀俊, 姚运生, 曾佐勋. 三峡水库蓄水对库区孕震环境及潜在震源影响研究. 大地测量与地球动力学 , 2005, 25(3): 116–120. Chen S J, Yao Y S, Zeng Z X. Effect of impounding on seismogenic environment and potential focus in head region of Three Gorge Reservoir. Journal of Geodesy and Geodynamics (in Chinese) , 2005, 25(3): 116-120. |
| [15] | Talwan P, Acree S. Pore pressure diffusion and the mechanism of reservoir-induced seismicity. Pure Appl. Geophys. , 1984, 85(122): 947-965. |
| [16] | Talwani P. Seismogenic properties of the crust inferred from recent studies of reservoir-induced seismicity-Application to Koyna. Current Science , 2000, 79(9): 1327-1333. |
| [17] | 雷兴林, 马胜利, 闻学泽, 等. 地表水体对断层应力与地震时空分布影响的综合分析-以紫坪铺水库为例. 地震地质 , 2008, 30(94): 1046–1064. Lei X L, Ma S L, Wen X Z, et al. Integrated analysis of stress and regional seismicity by surface loading-a case study of Zipingpu reservoir. Seismology and Geology (in Chinese) , 2008, 30(94): 1046-1064. |
| [18] | 周斌, 薛世峰, 邓志辉, 等. 水库诱发地震时空演化与库水加卸载及渗透过程的关系-以紫坪铺水库为例. 地球物理学报 , 2010, 53(11): 2651–2670. Zhou B, Xue S F, Deng Z H, et al. Relationship between the evolution of reservoir-induced seismicity in space-time and the process of reservoir water body load-unloading and water infiltration-a case study of Zipingpu reservoir. Chinese J. Geophys. (in Chinese) , 2010, 53(11): 2651-2670. |
| [19] | Nascimento A D, Lunn R J, Cowie P A. Modelling heterogeneous hydraulic properties of faults using constraints from reservoir-induced seismicity. Journal of Geophysical Research , 2005, 110: B09201. DOI:10.1029/2004JB003398 |
| [20] | Talwani P. Seismotectonics of the Koyna-Warna area, India. Pure Appl. Geophys. , 1997, 150(3-4): 511-550. DOI:10.1007/s000240050091 |
| [21] | Mandal P, Rastogi B K, Sarma C S P. Source parameters of Koyna earthquakes. India. Bull. Seismol. Soc. Am. , 1998, 88(3): 833-842. |
| [22] | Krishna V G, Kaila K L, Reddy P R. Synthetic seismogram modeling of crustal seismic record sections from the Koyna DSS profiles in western India. Monogr. IUGG , 1989, 51(6): 143-157. |
| [23] | 陈乐寿, 王光锷. 大地电磁测深法. 北京: 地质出版社, 1990 . Chen L S, Wang G E. Method of Magetotelluric Sounding (in Chinese). Beijing: Geological Publishing House, 1990 . |
| [24] | 徐常芳. 壳内流体演化及地震成因(三). 地震学报 , 1997, 19(2): 139–144. Xu C F. Crustal fluid evolution and causes of earthquakes (3). Acta Seismologica Sinica (in Chinese) , 1997, 19(2): 139-144. |
| [25] | 赵国泽, 刘国栋, 詹艳, 等. 张北-尚义地震区及其邻区地壳上地幔结构. 地震地质 , 1998, 20(2): 155–162. Zhao G Z, Liu G D, Zhan Y, et al. The structure of the crust and upper mantle for the Zhangbei-Shangyi seismic area and its neighbouring region. Seismology and Geology (in Chinese) , 1998, 20(2): 155-162. |
| [26] | 詹艳, 赵国泽, 王继军, 等. 1927年古浪地震区及其两侧地块的深部电性结构. 地球物理学报 , 2008, 51(2): 220–229. Zhen Y, Zhao G Z, Wang J J, et al. Deep electric structure beneath the epicentre of the 1927 Gulang M8 earthquake and its adjacent areas from magnetotelluric sounding. Chinese J. Geophys. (in Chinese) , 2008, 51(2): 220-229. |
| [27] | 詹艳, 赵国泽, 陈小斌, 等. 宁夏海原大震区西安州-韦州剖面大地电磁探测与研究. 地球物理学报 , 2004, 47(2): 305–313. Zhan Y, Zhao G Z, Chen X B, et al. Crustal structure from magnetotelluric profiling in the Haiyuan earthauake area, Ningxia Hui Autonomous region, China. Chinese J. Geophys. (in Chinese) , 2004, 47(2): 305-313. DOI:10.1002/cjg2.v47.2 |
| [28] | Ogawa Y M, Mishina T, Goto H, et al. Magnetotelluric imaging of fluids in intraplate earthquake zones, NE Japan back arc. Geophys. Res. Lett. , 2001, 28(19): 3741-3744. DOI:10.1029/2001GL013269 |
| [29] | Unsworth M, Bedrosian P A. On the geoelectric structure of major strike-slip faults and shear zones. Earth Planets Space , 2004, 56: 1177-1184. DOI:10.1186/BF03353337 |
| [30] | Tank S B, Honkura Y, Ogawa Y, et al. Resistivity structure in the western part of the fault rupture zone associated with the 1999. Izmit earthquake and its seismogenic implication.. Earth Planets Space , 2003, 55: 437-442. DOI:10.1186/BF03351777 |
| [31] | Ingham M. High resolution electrical image of fault zone. Physics of the Earth and Planetary Interiors , 2005, 150(1-3): 93-105. DOI:10.1016/j.pepi.2004.08.017 |
| [32] | Philip E, Wannamaker T, Grant Caldwell, et al. Fault zone fluids and seismicity in compressional and extensional environments inferred from electrical conductivity: the New Zealand Southern Alps and U. S. Great Basin. Earth Planets Space , 2004, 56: 1171-1176. DOI:10.1186/BF03353336 |
| [33] | 史水平, 于学庆, 龙政强, 等. 龙滩库区蓄水后地震性质分析. 地震地磁观测与研究 , 2010, 31(3): 40–45. Shi S P, Yu X Q, Long Z Q, et al. Analysis on the nature of earthquake in Longtan area after the reservoir storage water. Seismological and Geomagnetic Observation and Research (in Chinese) , 2010, 31(3): 40-45. |
| [34] | 广西龙滩地震台网地震目录月报. 2006-2010. Monthly Earthquake Catalog from Guangxi Longtan Seismic Network. 2006-2010. |
| [35] | 陈翰林, 赵翠萍, 修济刚, 等. 龙滩水库地震精定位及活动特征研究. 地球物理学报 , 2009, 52(8): 2035–2043. Chen H L, Zhao C P, Xiu J G, et al. Study on precise relocation of Longtan reservoir earthquakes and its seismic activity. Chinese J. Geophys. (in Chinese) , 2009, 52(8): 2035-2043. |
| [36] | 广西壮族自治区地质局. 南丹幅地质图. 北京: 地质出版社, 1968 . Guangxi Bureau of Geology. Nandan Geological Map (in Chinese). Beijing: Geological Publishing House, 1968 . |
| [37] | 向宏发, 周庆. 广西红水河龙滩水电站工程场地地震动参数复核. 中国地震局地质研究所工作报告,2006. Xiang H F, Zhou Q. Review report of ground motion parameters of red river Longtan Hydropower Station in Guangxi. Institute of Geology, China Seismological Bureau (in Chinese),2006. |
| [38] | Egbert G D, Booker J R. Robust estimation of geomagnetic transfer functions. Geophys. J. R. Astr. Soc. , 1986, 87(1): 173-194. DOI:10.1111/gji.1986.87.issue-1 |
| [39] | 陈小斌, 赵国泽, 詹艳. MT资料处理与解释的Windows可视化集成系统. 石油地球物理勘探 , 2004, 39(S 1): 11–16. Chen X B, Zhao G Z, Zhan Y. A visual intergrated windows system for MT data process and interpretation. Oil Geophysical Prospecting (in Chinese) , 2004, 39(S 1): 11-16. |
| [40] | 詹艳, 王立凤, 王继军, 等. 广西龙滩库区深部孕震结构大地电磁探测研究. 地球物理学报 , 2012, 55(4): 1400–1410. Zhan Y, Wang L F, Wang J J, et al. Electromagnetic survey of the seismogenic structures beneath the Longtan reservoir in Guangxi province. Chinese J. Geophys. (in Chinese) , 2012, 55(4): 1400-1410. |
| [41] | 万力, 蒋小伟, 王旭升. 含水层的一种普遍规律: 渗透系数随深度衰减. 高校地质学报 , 2010, 16(1): 7–12. Wan L, Jiang X W, Wang X S. A common regularity of aquifers: The decay in hydraulic conductivity with depth. Geological Journal of China Universities (in Chinese) , 2010, 16(1): 7-12. |
| [42] | Zhou L Q, Zhao C P, Zheng X, et al. Inferring water infiltration in the Longtan reservoir area by three-dimensional attenuation tomography. Geophys. J. Int. , 2011, 186(3): 1045-1063. DOI:10.1111/j.1365-246X.2011.05124.x |
| [43] | Zhou L Q, Zhao C P, Chen Z L, et al. Three-dimensional VP and VP/VS structure in the Longtan reservoir area by local earthquake tomography. Pure Appl. Geophys , 2011, 169(1-2): 123-139. |