三峡库区巴东地震前后视应力变化特征
戴苗1,2,3, 吴桂桔1,*, 刘坚1, 申重阳1, 孙少安1, 申学林1,2,3, 魏贵春1,2,3
1 中国地震局地震研究所, 地震大地测量重点实验室, 武汉 430071
2 中国地震局地震研究所, 地震预警湖北省重点实验室, 武汉 430071
3 湖北省地震局, 武汉 430071
*通讯作者: 吴桂桔, 副研究员, E-mail:wugjsky@126.com

〔作者简介〕 戴苗,男, 1982年生, 2008年于中国地质大学获固体地球物理学硕士学位, 工程师, 主要从事测震学理论、 方法和应用研究, E-mail:theocean123@163.com

摘要

三峡库区地震成因复杂, 为了研究该区域地震的成因及发震规律, 文中利用湖北省区域台网、 三峡台网提供的数字波形资料, 计算2010年1月1日至2015年12月31日三峡库区394个 ML≥2.0地震的震源参数, 获得了三峡库区视应力的时空分布图像, 分析了三峡库区研究时间段显著地震事件前后的视应力时空演化特征。研究表明: 1)巴东 ML5.5、 秭归 ML4.7以及 ML5.1地震发生前后, 不同震级分段计算结果显示, 震前震中区附近新华-水田坝断裂和高桥断裂同时出现高值, 高值区域分布呈现协同化程度高, 而震后零散的现象, 表明该区域在震前积累了较高的应力。2)研究区域视应力在地震前后经历了一个明显的 “上升-下降”过程, 在应力上升初期发生了巴东 ML5.5、 秭归 ML4.7以及 ML5.1 三次显著地震。3)研究区域视应力深度剖面显示, 不同深度视应力与震级大小呈现正相关性, 震后研究区内小震属于低应变释放背景下的低应力释放, 未新增明显视应力异常集中区域, 震后巴东高桥断裂、 周家山-牛口断裂以及秭归仙女山断裂短期内发生破坏性地震的可能性不大。

关键词: 视应力; 时空特征; 三峡库区; 巴东地震
中图分类号:P315.72 文献标志码:A 文章编号:0253-4967(2017)04-0837-16
APPARENT STRESS VARIATION CHARACTERISTICS BEFORE AND AFTER BADONG EARTHQUAKE IN THE THREE GORGES RESERVOIR AREA
DAI Miao1,2,3, WU Gui-ju1, LIU Jian1, SHEN Chong-yang1, SUN Shao-an1, SHEN Xue-lin1,2,3, WEI Gui-chun1,2,3
1 Key Laboratory of Earthquake Geodesy, Institute of Seismology, China Earthquake Administration, Wuhan 430071, China
2 Hubei Key Laboratory of Earthquake Early Warning, Institute of Seismology, China Earthquake Administration, Wuhan 430071, China
3 Earthquake Administration of Hubei Province, Wuhan 430071, China
Abstract

The causes of earthquakes in the Three Gorges reservoir area are complex. In order to study the cause of earthquakes happening in the region, we calculated the source parameters of 394 ML≥2.0 earthquakes occurring in Three Gorges reservoir area based on waveform data observed by the regional seismic network of Hubei Province and Three Gorges, obtained the apparent stress spatial-temporal variation map of Three Gorges reservoir area, and analyzed apparent stress spatial-temporal variation characteristics before and after the main earthquakes in the Three Gorges reservoir area. The results show that: 1)Before the Badong ML5.5, Zigui ML4.7 and Zigui ML5.1 earthquake, high apparent stress of earthquakes with different magnitudes is concentrated in Xinhua-Shuitianba Fault and Gaoqiao Fault. The distribution of high value area shows the high degree of synergism before the earthquake and the scattering after the earthquake, which indicates that the area accumulated a high stress before the earthquake, and the fault was in a locked state; 2)In the study area, the apparent stress before and after the earthquake showed significant rise in the first and then decline, the earthquake occurred in the process of rising; 3)Apparent stress depth profiles show that apparent stress at different depths has a positive correlation with the size of magnitude of earthquake, and the phenomenon of “small magnitude and strong apparent stress” did not appear. The small earthquakes occurring after the major earthquakes in the study area belong to low strain release under the background of low stress release, and there are no new apparent stress anomaly concentration areas appearing, this indicates that the Badong-Gaoqiao Fault, Zhoujiashan-Niukou Fault and Zigui Xiannushan Fault have been effective in releasing after the Badong earthquake and Zigui earthquakes and the probability of destructive earthquake is small on these faults.

Keyword: apparent stress; spatial-temporal feature; Three Gorges; Badong earthquake
0 引言

在震源介质为均匀弹性的脆性性质假定下, 地震视应力可由地震波信息获得, 是1个与区域平均应力呈正比的物理量, 与地震辐射能量和标量地震矩有关(Choy et al., 1995; 吴忠良, 2001; 吴忠良等, 2002)。在传统的地震学里, 地震辐射能量和标量地震矩只能利用震级通过经验公式换算得到, 误差较大且不能相互独立。目前, 随着中国观测技术水平的不断提升, 特别是宽频带数字地震计的使用, 使得利用数字波形资料反演地震震源谱的拐角频率、 低频水平、 地震应力降、 标量地震矩以及震源半径等震源参数成为可能, 利用这些参数可以获取相对独立的地震辐射能量, 进而获取地震视应力的大小。考虑到地震震级对视应力计算结果的影响, 具体在利用视应力研究区域应力场变化时, 为尽可能地保证视应力能真实反映区域应力场的变化, 采用震级分段研究是经国内实践验证的最直接也是最有效的手段和方法(李艳娥等, 2015; 陈丽娟等, 2015; 岳晓媛等, 2016)。近年来, 王鹏等(2014, 2015)尝试利用视应力值和近震震级的拟合关系, 将计算得到的视应力值减去近震震级拟合应力值, 得到视应力差值(即差视应力), 以期能排除震级对视应力的影响, 最终去除视应力的 “ 线性趋势” , 但由于能量的测定和地震矩的测定都存在一定误差, 该方法又高度依赖统计获取的近似拟合关系, 2种效应叠加在一起, 获得的结果误差可能会更大。

地震前的应力集中, 可通过视应力间接反映, 视应力可作为地壳应力水平的下限估计。研究表明视应力在中强地震前会有升高的趋势, 可作为震后趋势快速判定的依据(李艳娥等, 2012a, b; 陈丽娟等, 2015)。云南地区中小地震的视应力值在云南地区的部分地震序列中有很好的前兆特征表现, 在部分地区, 地震序列中只要发生了视应力值超过1MPa的中小地震, 其后就一定会有中强震发生; 如果序列中中强地震发生后再没有发生视应力值超过1MPa的中小地震, 则该序列就不会有强余震发生(刘红桂等, 2006)。2006年河北文安5.1级地震前2a开始震中区附近出现了明显的视应力升高过程(陈学忠等, 2007)。龙门山断裂带上小震频次高值区视应力在2003-2005年间, 曾出现过2次较为明显的上升-下降变化过程, 汶川地震临震前, 视应力出现突然上升(李艳娥等, 2012b)。

2013年12月16日, 湖北巴东县发生ML5.5地震, 随后不到4个月在相距不到50km的秭归再次发生ML4.9、 ML5.1地震, 这在湖北有历史地震记录以来是相当罕见的。3次地震震中均位于秭归向斜SW侧, 即黄陵背斜之西翼马蹄形盆地的西南边缘, 其区域主压应力场方向为NW-SE向。巴东地震震源机制解表明, 此次巴东地震的主张应力轴走向及主压应力轴走向和震源错动性质与蓄水之前巴东区域的构造应力场背景有近似性, 本次主震的最大主压应力轴和最大主张应力轴走向与蓄水前几乎一致, 部分学者认为巴东地震活动与震中区域NE向断裂的区域构造应力调整有关, 属于构造地震(陈俊华等, 2014)。有限元数值模拟巴东地震前库区等效应力场的变化结果显示, 巴东地震震中位于蓄水引起的巴东与牛口等效应力差值区的梯度带上, 故部分学者指出巴东地震与蓄水相关, 震源区附近发育大规模条带岩溶, 认为该地震是区域构造应力场叠加蓄水引起的局部构造应力状态调整触发的塌陷型水库地震(陈蜀俊等, 2014)。此外, 对于此次地震的发震构造, 由于湖北巴东地震及较大余震震源机制解显示, 近EW向节面为巴东地震发震断层, 部分学者认为巴东地震发震构造可能为EW向大坪断裂(赵凌云等, 2014)。与此同时, 部分学者根据巴东地震区域地震构造环境和地震烈度分布等结论分析认为, 巴东高桥断裂、 马鹿池断裂及周家山断裂围限的楔形微块体构成巴东地震的发震单元(雷东宁等, 2014), 该区域浅表构造和深部构造非一一对应关系, 故地表找不到相应的发震构造。

巴东地震是三峡2003年蓄水以来库区最强地震, 开展巴东地震前后三峡库区地震震源参数及视应力的相关研究, 对深入了解巴东地震成因、 发震构造以及获取三峡库区近期后续地震趋势具有重要意义。

1 地震构造背景及资料的选取

研究范围属三峡库区地震重点监视区域, 区域构造主要以NNE、 NW向为主, 如NNE向新华-水田坝断裂(F3)、 九畹溪断裂(F5)、 周家山断裂(F2)、 高桥断裂(F1), NW走向的断裂主要有仙女山断裂(F4)、 天阳坪断裂(F7)等(图1)。其中, 新华-水田坝断裂新华段属基底断裂, 早期具左旋走滑性质; 水田坝段属盖层断裂, 断层面W倾, 倾角约70° , W盘向S运动, 断层具张扭性质, 第四纪以来活动性不明显。仙女山断裂走向340° ~355° , 早期表现为逆冲性质, 断面向E或W陡倾; 深部资料显示该断裂未有越过长江的迹象; 断裂最新活动时间为早、 中更新世(雷东宁等, 2014); 以右旋走滑兼正断倾滑性质为主; 2008年11月26日, 该断裂附近发生了三峡蓄水后首次4级以上地震, 震级达ML4.6。天阳坪断裂走向320° , 倾向SW, 倾角30° ~70° ; 西接仙女山断裂, 西段倾角较陡, 剖面表现为逆冲叠瓦状冲断组合; 第四纪以来活动性较弱。周家山-牛口断裂发育于震中东部, 总体走向20° , 倾向300° , 倾角约70° ; 此断裂为周家山向斜内部的1条横断裂, 表现为左行平移逆冲活动的特点, 喜山期以来为张性正断性质; 其最新活动年代测定为晚更新世早期(华卫等, 2010)。此次巴东地震震中位置正好位于高桥断裂与周家山-牛口断裂之间。

本文主要选用湖北区域台网、 三峡台网记录的2010年1月1日至2015年12月31日共计394个ML≥ 2.0地震的数字地震波形资料; 该时间段内ML2.0以上地震月频次相对较为稳定, 频次较高, 且台网观测站分布均匀, 选用台站及地震事件时空分布见图1, 2。选用台站仪器型号主要为BBVS-60、 FSS-3M和少量JC-V104-3D, 观测频带主要为0.017~40Hz、 0.5~70Hz, 个别为1~40Hz, 采样率均为100Hz。在近震源条件下, 计算时主要选用震中距在25~150km以内的台站, 选取信噪比较高的波形数据, 并对这些数据进行去倾和仪器响应校正处理。

图1 研究区地震活动性、 测震台站分布及构造分布图(2010-01-01-2015-12-30)
F1高桥断裂; F2周家山断裂; F3新华-水田坝断裂; F4仙女山断裂; F5九畹溪断裂; F6亩田湾断裂; F7天阳坪断裂
Fig. 1 Distribution map of earthquakes, seismic stations and faults in the Three Gorges reservoir area(2010-01-01-2015-12-30).

图2 三峡库区 2010年1月1日至2015年12月30日地震活动性参数图Fig. 2 Seismicity parameters in the Three Gorges reservoir area during 2010-01-01-2015-12-30.

2 计算方法及过程

首先, 将波形数据进行傅里叶变换计算观测谱。然后, 从观测谱中扣除传播路径效应和场地效应等因素的影响, 得到地震震源位移谱; 对于中小地震, 震源位移谱符合Brune圆盘模型, 震源位移谱 Ω(f)可表示为(李艳娥等, 2012a, b; 陈丽娟等, 2015)

Ω(f)=Ω01+ffc2(1)

式(1)中, Ω0为离开震源质心单位距离处的震源位移谱低频水平; fc为震源位移谱的拐角频率。给定 Ω0fc即可确定震源谱。可采用以下方法确定 Ω0fc

本文采用的数据为三分向速度记录, 各个分向可以得到相应的速度谱, 总的速度谱是各个分量速度谱的矢量合成, 即

V(f)=V(f)Z2+V(f)NS2+V(f)EW2(2)

速度谱和位移谱 Ω(f)的关系为

V(f)=2πfΩ(f)(3)

所以可得到位移谱:

Ω(f)=V(f)Z2+V(f)NS2+V(f)EW22πf(4)

在实际计算中, 选取谱形状较好的台站进行计算, 根据实际情况选取位移谱的低频段水平部分, 低频范围假设为f1~f2, 同时选取高频段的衰减起始频率f3 。在低频段(f1~f2), Ω (f)≈ Ω 0, 则位移谱平方的积分为

SD=f1f2Ω(f)2df=Ω02f2-f1(5)

Ω0=SDf2-f1(6)

对于实际位移谱的低频段水平部分, 用式(5)进行积分, 根据式(6)可得到位移谱零频极限值 Ω0得到 Ω0后, 可以通过如下方法确定拐角频率 fc, f1~f3频段范围, 将得到的 Ω0代入式(1)可得到这个频段内的理论值Ω ti, 而同频段内的实际观察谱为Ω i, 二者之间的距离可由式(7)计算:

RΩ=i=1N(Ωi-Ωti)2(7)

RΩ随频率变化, 当其最小时, 即说明理论值与实际观测值之间的偏差最小; 因此, 将 RΩ为最小值时的 f值作为拐角频率 fc

中小地震的地震矩可利用式(8)计算:

M0=4πρυ3Ω0dR(8)

式(8)中, ρ为地壳密度; υ为S波速度; d为震源距; R为辐射因子。

对于1条数字化的速度波形地震记录, 在进行仪器响应校正、 自由表面效应校正、 介质衰减校正后, 其地震辐射能量可表示为

Es=4πρβ·20V(f)2df(9)

具体在计算的过程中, 考虑到中小地震含有较多的高频成分, 有限的频带宽度可能会导致地震辐射能量低估, 所以在计算过程中还需考虑低频和高频补偿。在得到地震矩 M0和地震能量 Es后, 即可计算地震视应力:

σa=μEsM0(10)

式(10)中, μ为剪切模量。具体对某次地震时, 一般会有多个台站记录到波形资料, 可根据上述方法分别求得每个台站的震源力学参数, 再由各台站的值求平均得到多台平均值。

采用上述方法, 我们计算了研究区域2010年以来近6a内ML2.0以上地震的拐角频率、 地震矩、 地震能量以及视应力等震源动力学参数, 表1中给出了研究区域部分地震(ML≥ 2.7地震)的震源动力学参数计算结果。其中, 研究区域ML2.0以上地震视应力平均值为0.05MPa, 这和华东地区视应力场平均值0.08MPa(黄显良等, 2005)以及江西地区平均视应力0.04MPa(汤兰荣等, 2007)非常接近。

表1 研究区域ML≥ 2.7地震的震源动力学参数计算结果 Table1 Seismic source dynamic parameters of ML≥ 2.7 earthquakes in the study area
3 三峡库区地震震源参数特征

地震震源参数与震级的关系国内外已有相关研究, 但地震震源参数和震级的相关性存在着地区差异, 图3给出了三峡库区ML2.0以上地震事件各震源参数与震级的统计关系。

图3 研究区震源动力学参数与震级之间的关系Fig. 3 Relationship between seismic source dynamic parameters and magnitude.

3.1 地震矩与震级的关系

图3a为研究区域地震矩与震级之间的定标关系图。从计算结果来看, 三峡库区地震矩的对数和震级基本呈现线性关系, 两者统计关系为 logM0=0.85×ML+13.6表明该地区随着震级的增大, 地震矩也相应增大。定标关系的斜率可反映研究区域介质的非均匀度(钟羽云等, 2004)。三峡地区地震矩与震级定标关系的斜率为0.85, 低于中国平均水平(易桂喜等, 2011; 阮祥等, 2011; 王鹏等, 2014, 2015; 李艳娥等, 2015), 表明该地区介质的均匀度较好。

3.2 地震辐射能量、 视应力与震级的关系

图3b和图3c分别为三峡地区地震辐射能量、 视应力与震级之间的定标关系图, 根据计算结果, 它们的定标关系分别为 logEs=1.392×ML+8.213, logσa=0.548×ML-3.004三峡库区地震的辐射能量以及视应力均随震级的增大而增大。表明震级越大, 地震辐射的能量越多, 视应力越高; 如果个体地震与视应力出现非关联性, 可以考察其视应力是否存在异常特征。地震视应力随震级增大而增大意味着地震越大, 单位地震矩所释放的地震能量越多, 反映出大地震是比小地震具有更高效率的地震能量辐射体, 大地震和小地震可能在破裂速度等方面有较大差异。

3.3 地震拐角频率与震级的关系

图3d为三峡地区地震拐角频率与震级之间的定标关系图, 定标关系为logfc=-0.08× ML+0.725。统计结果显示, 中小地震的拐角频率较大, 大地震往往拥有较小的拐角频率。此外, 不同类型地震的 P波和S波拐角频率也会有差别(张丽芬等, 2013), 统计结果显示的三峡地区拐角频率与震级之间的定标关系斜率为-0.08, 计算结果低于中国其他非水库区域的计算结果(易桂喜等, 2011; 阮祥等, 2011; 王鹏等, 2014; 李艳娥等, 2015; 王鹏等, 2015), 这可能和三峡库区ML2.0~2.5震级段的地震中含有部分水库诱发的塌陷型地震有关。

4 研究区域视应力的时空变化特征
4.1 地震前后视应力的空间分布

2013年12月16日湖北巴东发生ML5.5地震, 2014年3月27日和3月30日, 湖北秭归先后发生ML4.7和ML5.1地震, 4个月内, 接连发生3次ML4.5以上地震。 以2013年12月16日湖北巴东ML5.5地震和2014年3月30日湖北秭归ML5.1地震为时间节点, 将时间段分成地震前、 地震后及地震时3个时间段, 为消除余震对震后计算结果的影响, 具体计算时震后的时间起始时间从2014年4月10日开始计算。以0.05° × 0.05° 的空间窗, 对视应力进行空间扫描, 可以得到研究区域视应力的空间分布(图4)。

图4 研究区巴东、 秭归地震前后ML≥ 2.0及ML≥ 2.5地震视应力空间分布Fig. 4 The spatial distribution of apparent stress with ML≥ 2.0 and ML≥ 2.5 earthquakes in the study area before and after the Badong earthquake and Zigui earthquake.

结果显示, 对于ML≥ 2.0的地震, 巴东和秭归地震前视应力高值沿新华-水田坝断裂(F3)及高桥断裂(F1)走向展布, 协同化程度较高, ML≥ 2.5地震视应力主要集中在新华-水田坝断裂(F3)带上。震后, 无论是ML≥ 2.0地震, 还是ML≥ 2.5以上地震的视应力高值分布均较为分散和无序。2013年12月16日湖北巴东ML5.5地震发生在ML≥ 2.0地震视应力震前高桥断裂(F1)高值异常体的边缘, 而2014年3月27日和3月30日湖北秭归ML4.7和ML5.1地震均发生在新华-水田坝断裂(F3)异常条带连线的南端附近。巴东和秭归地震后, ML≥ 2.0地震和ML≥ 2.5地震计算所得视应力结果均显示, 受巴东地震余震的影响, 巴东地震震源体周边出现了非常显著的视应力高值区。

为了消除震级跨度对视应力计算结果的影响, 综合考虑湖北省不同强度地震频次的实际情况, 我们将研究区域地震震级按照ML2.0~2.4、 ML2.5~3.5分段, 通过不同震级段视应力的空间分布研究三峡地区巴东和秭归3次显著地震事件前震源体附近视应力的高值分布, 以期最终获得具有指示意义的震级段视应力高值异常。分段计算的结果如图5所示。

巴东地震前, ML2.0~2.4、 ML2.5~3.5 两个震级段的视应力高值分布在新华-水田坝断裂(F3)和巴东高桥断裂(F1)上。震后, ML2.0~2.4震级段高桥断裂(F1)高值异常消失, 新华-水田坝断裂(F1)和周家山断裂(F2)出现新的高值区; 受余震的影响, ML2.5~3.5震级段视应力在巴东地震震源体附近的高桥断裂(F1)中段以及周家山断裂(F2)出现大面积高值。

4.2 研究区域地震视应力随时间的变化

为剔除主震本身以及震级差异性对计算结果的影响, 我们将震级分为ML2.0~2.4、 ML2.5~3.5两个震级段进行研究, 由图6可知, ML2.5~3.5 震级段地震视应力平均值要大于ML2.0~2.4震级段, 视应力随着震级的增加而增加, 这与前人的研究结果是一致的, 也间接证实了本文计算视应力方法的可行性以及震级分段研究的必要性。

图5 研究区域巴东、 秭归地震前后ML2.0~2.4和ML2.5~3.5震级范围内地震视应力空间分布Fig. 5 The spatial distribution of apparent stress with ML2.0~2.4 and ML2.5~3.5 earthquakes in the study area before and after the Badong earthquake and Zigui earthquake.

由图6a所示, 研究区域内ML2.0~2.4地震视应力2011年之前一直较为稳定, 2011年1月份开始下降, 2014年11月底降至最低, 然后开始缓慢抬升, 上升过程中先后发生2013年12月16日湖北巴东ML5.5地震, 2014年3月27日、 3月30日湖北秭归ML4.7、 ML5.1地震, 震后研究区域视应力开始缓慢下降, 2015年6月左右恢复到震前水平。 ML2.5~3.5地震视应力计算结果受研究区域地震个数的影响, 曲线的采样率虽然不高, 但从曲线的整体趋势仍可以看出2012年之前研究区域视应力一直较为稳定, 2012年1月份开始下降, 2013年6月份左右降至最低, 然后开始抬升, 同样也是在上升的过程中先后发生巴东地震和秭归地震, 震后2015年下半年度视应力恢复到震前的平均水平。研究区域ML2.5~3.5震级段地震视应力震前开始抬升的拐点比ML2.0~2.4时间段计算结果提前(图6b), 这可能是受该时间段ML2.5~3.5震级段地震数量少导致采样率低所致; 即便如此, 从图像的整体变化特征仍能看出ML2.0~2.4、 ML2.5~3.5两个震级段地震视应力随时间的变化同步性较好。研究区域视应力经历了明显的 “ 上升-下降” 变化过程。这种变化过程说明视应力在研究区域内巴东和秭归地震前有明显的升高异常变化。

图6 研究区域不同震级段视应力随时间的变化Fig. 6 The changes of apparent stress with different magnitudes in the study area.

图7 研究区域视应力沿新华-水田坝断裂深度剖面分布图(2010-01-01-2015-12-31)Fig. 7 The depth profiles of apparent stress along Xinhua-Shuitianba Fault in the area(2010-01-01-2015-12-31).

图8 研究区域视应力沿长江深度剖面分布图(2010-01-01-2015-12-31)Fig. 8 Depth profiles of apparent stress along the Yangtze River in the study area(2010-01-01-2015-12-31).

4.3 研究区域视应力深度剖面分布特征

基于研究区域ML2.5以上地震精定位结果, 将NE向新华-水田坝断裂(F3, 图1)以及长江附近50km范围内地震分别投影到沿断裂和长江的剖面上分析视应力的深度变化特征, 以新华-水田坝断裂和长江为横坐标, 震源深度为纵坐标, 色标为视应力的大小, 地震资料选用时间段从2010年1月1日至2015年12月31日。从图7, 8可以看出, 剖面在两河口、 水田坝以及信陵镇和郭家坝下方出现视应力的高值集中区域, 这跟研究区内秭归双震和巴东地震视应力相关。此外, 值得关注的是, 图7显示巴东及秭归地震震源体附近或周边出现了地震视应力高值副中心, 经核实主要和巴东地震以及秭归地震各自的几次ML3.0以上的余震有关。震后, 2条深度剖面视应力的高值和震级基本呈现正相关, 属于低应变释放背景下的低应力释放, 震后未出现 “ 小震级大视应力” 的现象, 未新增视应力高值异常区域。

4.4 研究区水位涨落期间视应力变化特征

为研究三峡蓄水对研究区域视应力的影响, 选取低水位期和高水位期2个时间段进行对比研究。低水位与高水位期视应力计算结果分别为2010-2015年6a低水位期(6月1日至8月31日)76个和高水位期(10月1日至次年1月31日)127个地震视应力的累计结果。低水位运行期坝前水位维持在145m左右, 高水位运行期坝前水位为165~175m。

计算结果显示: 低水位运行期间, 兴山县境内视应力局部出现高值, 整个库区其他区域视应力分布较为均匀; 高水位运行期间, 巴东高桥断裂沿线出现高值; 值得关注的是, 巴东、 秭归长江沿线库岸10km范围内视应力分布和水位高低无明显关系。

图9 研究区域水位涨落期间ML2.0~3.0震级范围内地震视应力空间分布Fig. 9 The spatial distribution of apparent stress with ML2.0~3.0 earthquakes in the study area before and after the impoundment of the Three Gorges Reservoir.

5 结论

三峡库区地震成因较为复杂, 既有天然构造地震, 也有因蓄水导致的塌陷、 边坡失稳等水库地震, 同时也包含蓄水因素诱发的构造地震。通过对三峡库区2010年1月1日至2015年12月31日ML2.0以上地震震源参数的计算, 得到了三峡库区ML2.0以上地震的地震矩、 辐射能量、 拐角频率、 视应力等震源参数, 三峡库区主要震源参数和震级的统计结果以及三峡库区视应力的时空分布总结得出以下结论:

(1)三峡地区拐角频率的对数和震级的定标关系斜率为-0.08, 该值低于中国非水库地区的计算结果, 这与三峡库区ML2.0~2.5震级段地震含有部分塌陷地震事件有关。

(2)三峡库区视应力空间分布显示, 巴东地震以及秭归地震震前, 三峡库区不同震级段地震视应力在巴东和秭归地震震源体附近或周边出现高值集中现象, 高值分别沿高桥断裂带和新华-水田坝断裂分布且协同化程度震前较高, 震后分布较为零散, 表明该区域在震前积累了较高的应力。巴东地震和秭归地震震中分别位于高桥断裂高值异常体边缘和新华-水田坝断裂异常条带延长线附近, 这可能与震源体在应力的长期连续作用下, 弹性应变不断积累, 导致构造的突然失稳破裂有关。考虑到地震孕育的复杂性, 单纯的依赖视应力的空间分布, 很难确定具体的发震构造。

(3)三峡库区视应力的时间序列显示, 从2013年底至2015年初, 研究区域应力经历了1个 “ 上升-下降” 过程, 在应力上升初期发生了巴东ML5.5、 秭归ML4.7以及ML5.1等3次显著地震。

(4)一般情况下, 震级相当时, 视应力值相对高的震源区应力水平高; 视应力值相等时, 震级较低的震源区应力水平高。研究区域视应力的深度投影显示, 巴东及秭归地震震后, 该区域中小地震的视应力和震级呈现较好的正相关性, 未出现 “ 小地震大视应力” 的现象, 未新增明显视应力异常集中区域, 目前三峡库区仍处于巴东MS5.1地震震后的调整阶段, 综合考虑三峡库区历史地震背景, 短期内三峡库区发生破坏性地震的可能性不大。

(5)三峡水位不同区间视应力对比分析结果显示, 三峡水位涨落对三峡库区长江沿岸视应力分布影响有限。表明巴东、 秭归地震前库区视应力空间异常和三峡蓄水关联不大, 巴东、 秭归地震的主要成因为区域构造应力调整, 巴东、 秭归地震均为构造地震。

致谢 中国地震局地球物理勘探中心姚家骏老师对本文提供了全程指导, 审稿专家提出的修改意见对稿件质量的提升帮助很大, 在此表示衷心感谢。

The authors have declared that no competing interests exist.

参考文献
[1] 陈俊华, 陈正松, 蒋玲霞, . 2014. 巴东MS5. 1地震成因研究[J]. 大地测量与地球动力学, 34(3): 10-14.
CHEN Jun-hua, CHEN Zheng-song, JIANF Ling-xia, et al. 2014. On causes of Badong MS5. 1 earthquake[J]. Journal of Geodesy and Geodynamics, 34(3): 10-14(in Chinese). [本文引用:1]
[2] 陈丽娟, 李艳娥, 杨立明, . 2015. 2013年7月22日岷县MS6. 6地震前后视应力研究[J]. 地震, 35(4): 30-42.
CHEN Li-juan, LI Yan-e, YANG Li-ming, et al. 2015. Temporal-spatial variations of apparent stress before and after the 22 July 2013 Minxian MS6. 6 earthquake[J]. Earthquake, 35(4): 30-42(in Chinese). [本文引用:3]
[3] 陈蜀俊, 姚运生, 吴建超, . 2014. 巴东MS5. 1地震: 一种新的水库地震类型[J]. 大地测量与地球动力学, 34(3): 1-5.
CHEN Shu-jun, YAO Yun-sheng, WU Jian-chao, et al. 2014. Badong MS5. 1 earthquake: A new type of reservoir earthquake[J]. Journal of Geodesy and Geodynamics, 34(3): 1-5(in Chinese). [本文引用:1]
[4] 陈学忠, 李艳娥. 2007. 2006年7月4日河北文安MS5. 1地震前震中周围地区小震视应力随时间的变化[J]. 中国地震, 23(4): 327-336.
CHEN Xue-zhong, LI Yan-e. 2007. Temporal changes in apparent stresses of small earthquakes around the epicentral area prier to the July 4, 2006 Wen'an, Hebei earthquake(MS=5. 1)[J]. Earthquake Research in China, 23(4): 327-336(in Chinese). [本文引用:1]
[5] 华卫, 陈章立, 郑斯华, . 2010. 三峡水库地区震源参数特征研究[J]. 地震地质, 32(4): 533-542. doi: 103969/j.issn.0253-4967.201004. 001.
HUA Wei, CHEN Zhang-li, ZHENG Si-hua, et al. 2010. A study on characteristics of source parameters in the Three Gorges reservoir area[J]. Seismology and Geology, 32(4): 533-542(in Chinese). [本文引用:1]
[6] 黄显良, 刘东旺, 沈小七, . 2005. 华东地区视应力场的研究[J]. 地震地磁观测与研究, 26(6): 29-32.
HUANG Xian-liang, LIU Dong-wang, SHEN Xiao-qi, et al. 2005. The study on the apparent stress field in East China[J]. Seismological and Geomagnetic Observation and Research, 26(6): 29-32(in Chinese). [本文引用:1]
[7] 雷东宁, 陈俊华, 张丽芬, . 2014. 湖北巴东MS5. 1地震构造背景分析[J]. 大地测量与地球动力学, 34(3): 6-9, 14.
LEI Dong-ning, CHEN Jun-hua, ZHANG Li-fen, et al. 2014. Analysis of seismogenic background of Badong MS5. 1 earthquake[J]. Journal of Geodesy and Geodynamics, 34(3): 6-9, 14(in Chinese). [本文引用:2]
[8] 李艳娥, 陈丽娟, 王生文, . 2015. 山东地区地震视应力时空变化特征研究[J]. 地震, 35(2): 80-90.
LI Yan-e, CHEN Li-juan, WANG Sheng-wen, et al. 2015. Temporal and spatial variations of apparent stress in Shand ong Province, China[J]. Earthquake, 35(2): 80-90(in Chinese). [本文引用:3]
[9] 李艳娥, 陈学忠, 付虹. 2012 a. 2007年云南泞洱MS6. 4地震前滇西南地区震源动力学参数时空变化特征[J]. 地震, 32(1): 28-39.
LI Yan-e, CHEN Xue-zhong, FU Hong. 2012a. Temporal and spatial variation of focal dynamic parameters in southwest Yunnan before the 2007 MS6. 4 Ning'er earthquake[J]. Earthquake, 32(1): 28-39(in Chinese). [本文引用:2]
[10] 李艳娥, 陈学忠, 王恒信. 2012 b. 汶川8. 0级地震前四川地区地震视应力时空变化特征[J]. 地震, 32(4): 113-122.
LI Yan-e, CHEN Xue-zhong, WANG Heng-xin. 2012b. Temporal and spatial variation of apparent stress in Sichuan area before the MS8. 0 Wenchuan earthquake[J]. Earthquake, 32(4): 113-122(in Chinese). [本文引用:2]
[11] 刘红桂, 刘杰, 丁页岭, . 2006. 地震视应力在云南地震序列中的前兆特征[J]. 地震学报, 28(5): 462-471.
LIU Hong-gui, LIU Jie, DING Ye-ling, et al. 2006. Precursory speciallties of apparent stresses in Yunnan earthquake series[J]. Acta Seismologica Sinica, 28(5): 462-471(in Chinese). [本文引用:1]
[12] 刘红桂, 王培玲, 杨彩霞, . 2007. 地震视应力在地震预测中的应用[J]. 地震学报, 29(4): 437-445.
LIU Hong-gui, WANG Pei-ling, YANG Cai-xia, et al. 2007. Application of apparent stress in earthquake prediction[J]. Acta Seismologica Sinica, 29(4): 437-445(in Chinese). [本文引用:1]
[13] 阮祥, 程万正, 乔惠珍, . 2011. 安宁河-则木河地震带震源参数及应力状态的研究[J]. 西北地震学报, 33(1): 46-50.
RUAN Xiang, CHENG Wan-zheng, QIAO Hui-zhen, et al. 2011. Research on source parameters and stress state in Anninghe-Zemuhe earthquake belt[J]. Northwestern Seismological Journal, 33(1): 46-50(in Chinese). [本文引用:2]
[14] 汤兰荣, 杨锐, 杨雅琼, . 2007. 九江-瑞昌地震序列的视应力变化研究[J]. 华南地震, 27(4): 19-27.
TANG Lan-rong, YANG Rui, YANG Ya-qiong, et al. 2007. A study on the variations in apparent stress of Jiujiang-Ruichang seismic sequence[J]. South China Journal of Seismology, 27(4): 19-27(in Chinese). [本文引用:1]
[15] 王鹏, 郑建常. 2014. 鲁东地区视应力时空变化特征分析[J]. 地震, 34(4): 70-77.
WANG Peng, ZHENG Jian-chang. 2014. Temporal and spatial variation of apparent stress in eastern Shand ong Province[J]. Earthquake, 34(4): 70-77(in Chinese). [本文引用:2]
[16] 王鹏, 郑建常, 刘希强, . 2015. 郯庐断裂带山东段震源参数及应力状态[J]. 地震地质, 37(4): 966-981. doi: 103969/j.issn.0253-4967.2015.04.003.
WANG Peng, ZHENG Jian-chang, LIU Xi-qiang, et al. 2015. Research of source parameters and stress state in Shand ong segment of Tanlu fault zone[J]. Seismology and Geology, 37(4): 966-981(in Chinese). [本文引用:2]
[17] 吴忠良. 2001. 由宽频带辐射能量目录和地震矩目录给出的视应力及其地震学意义[J]. 中国地震, 17(1): 8-15.
WU Zhong-liang. 2001. Apparent stress obtained from broadband radiated energy catalogue and seismic moment catalogue and its seismological significance[J]. Earthquake Research in China, 17(1): 8-15(in Chinese). [本文引用:1]
[18] 吴忠良, 黄静, 林碧苍. 2002. 中国西部地震视应力的空间分布[J]. 地震学报, 24(3): 293-301.
WU Zhong-liang, HUANG Jing, LIN Bi-cang. 2002. Distribution of apparent stress in western China[J]. Acta Seismologica Sinica, 24(3): 293-301(in Chinese). [本文引用:1]
[19] 易桂喜, 闻学泽, 辛华, . 2011. 2008年汶川MS8. 0地震前龙门山-岷山构造带的地震活动性参数与地震视应力分布[J]. 地球物理学报, 54(6): 1490-1500.
YI Gui-xi, WEN Xue-ze, XIN Hua, et al. 2011. Distributions of seismicity parameters and seismic apparent stresses on the Longmenshan-Minshan tectonic zone before the 2008 MS8. 0 Wenchuan earthquake[J]. Chinese Journal of Geophysics, 54(6): 1490-1500(in Chinese). [本文引用:2]
[20] 岳晓媛, 武安绪, 冯刚, . 2016. 首都圈中强震前后视应力时空演化过程研究[J]. 地震, 36(2): 119-131.
YUE Xiao-yuan, WU An-xu, FENG Gang, et al. 2016. Temporal-spatial evolution of apparent stresses before moderate-strong earthquakes in the capital circle area of China[J]. Earthquake, 36(2): 119-131(in Chinese). [本文引用:1]
[21] 张丽芬, 姚运生, 李井冈, . 2013. 三峡库区构造和塌陷地震的拐角频率特征[J]. 大地测量与地球动力学, 33(2): 27-30, 40.
ZHANG Li-fen, YAO Yun-sheng, LI Jing-gang, et al. 2013. Corner frequency characteristic of tectonic earthquakes and collapsed ones in Three Gorges region[J]. Journal of Geodesy and Geodynamics, 33(2): 27-30, 40(in Chinese). [本文引用:1]
[22] 赵凌云, 张辉, 陈俊华, . 2014. 基于CAP方法的2013湖北巴东5. 1级地震震源机制与发震构造研究[J]. 地震工程学报, 36(4): 1014-1018.
ZHAO Ling-yun, ZHANG Hui, CHEN Jun-hua, et al. 2014. Focal mechanism solutions and seismogenic structure of the Badong MS5. 1 earthquake of 2013 using the CAP method[J]. China Earthquake Engineering Journal, 36(4): 1014-1018(in Chinese). [本文引用:1]
[23] 钟羽云, 朱新运, 张震峰. 2004. 不同类型地震的地震矩: 震级标度关系研究[J]. 西北地震学报, 26(1): 57-61.
ZHONG Yu-yun, ZHU Xin-yun, ZHANG Zhen-feng. 2004. Study on relations between seismic moment and magnitude for various types of earthquake sequence[J]. Northwestern Seismological Journal, 26(1): 57-61(in Chinese). [本文引用:1]
[24] Choy G L, Boatwright J L. 1995. Global patterns of radiated seismic energy and apparent stress[J]. Journal of Geophysical Research, 100(B9): 18205-18228. [本文引用:1]