2017年11月23日重庆武隆 MS5.0地震序列重定位及发震断层分析
李翠平, 唐茂云, 郭卫英*, 黄世源, 王小龙, 高见
重庆市地震局, 重庆 401147
*通讯作者: 郭卫英, 女, 1969年生, 高级工程师, 研究方向为地震综合预报与大形势预测, E-mail: weiying9564@163.com

〔作者简介〕 李翠平, 女, 1989年生, 2015年于中国地震局地质研究所获构造地貌学专业硕士学位, 助理工程师, 主要从事地震活动性与地震危险性研究, 电话: 023-67086630, E-mail: wulicuiping@126.com

摘要

武隆 MS5.0地震是发生在七曜山-金佛山及彭水2条基底断裂所夹持的武隆凹陷束中的一次中强震, 震源区位于芙蓉江断层、 文复断层和马武断层的交会部位。文中采用CAP波形反演方法, 计算了武隆 MS5.0地震的震源机制解与震源矩心深度, 同时利用双差定位方法对武隆地震及其余震序列进行重定位, 分析了本次地震的发震断层。结果显示, 武隆 MS5.0地震震源机制解节面Ⅰ的走向为24°、 倾角为16°、 滑动角为-108°, 节面Ⅱ的走向为223°、 倾角为75°、 滑动角为-85°; 重定位后的余震序列主要分布在主震的SW向, 深度集中于5~8km。分析表明, 余震分布特征与震源机制解节面Ⅱ一致, 且符合文复断层走向SW、 倾向NW、 倾角60°~73°的产状特征。由此推测此次武隆 MS5.0地震的发震构造为文复断层, 并初步探讨了其可能的发震模式。

关键词: 武隆地震; 震源机制; 双差定位; 发震断层
中图分类号:P315.2 文献标志码:A 文章编号:0253-4967(2019)03-0603-16
RELOCATION OF THE 23 NOVEMBER 2017 WULONG MS5.0 EARTHQUAKE SEQUENCE AND ANALYSIS OF ITS SEISMOGENIC FAULT
LI Cui-ping, TANG Mao-yun, GUO Wei-ying, HUANG Shi-yuan, WANG Xiao-long, GAO Jian
Chongqing Earthquake Agency, Chongqing 401147, China
Abstract

The Wulong MS5.0 earthquake on 23 November 2017, located in the Wolong sap between Wenfu, Furong and Mawu faults, is the biggest instrumentally recorded earthquake in the southeastern Chongqing. It occurred unexpectedly in a weak earthquake background with no knowledge of dramatically active faults. The complete earthquake sequences offered a significant source information example for focal mechanism solution, seismotectonics and seismogenic mechanism, which is helpful for the estimation of potential seismic sources and level of the future seismic risk in the region. In this study, we firstly calculated the focal mechanism solutions of the main shock using CAP waveform inversion method and then relocated the main shock and aftershocks by the method of double-difference algorithm. Secondly, we determined the seismogenic fault responsible for the MS5.0 Wulong earthquake based on these calculated results. Finally, we explored the seismogenic mechanism of the Wulong earthquake and future potential seismic risk level of the region.
The results show the parameters of the focal mechanism solution, which are: strike24°, dip 16°, and rake -108° for the nodal plane Ⅰ, and strike223°, dip 75°, and rake -85° for the nodal plane Ⅱ. The calculations are supported by the results of different agencies and other methods. Additionally, the relocated results show that the Wulong MS5.0 earthquake sequence is within a rectangular strip with 4.7km in length and 2.4km in width, which is approximately consistent with the scales by empirical relationship ofWells and Coppersmith(1994). Most of the relocated aftershocks are distributed in the southwest of the mainshock. The NW-SE cross sections show that the predominant focal depth is 5~8km. The earthquake sequences suggest the occurrence features of the fault that dips northwest with dip angle of 63° by the least square method, which is largely consistent with nodal planeⅡof the focal mechanism solution. Coincidentally, the field outcrop survey results show that the Wenfu Fault is a normal fault striking southwest and dipping 60°~73° by previous studies. According to the above data, we infer that the Wenfu Fault is the seismogenic structure responsible for Wulong MS5.0 earthquake.
We also propose two preliminary genetic mechanisms of “local stress adjustment” and “fluid activation effect”. The “local stress adjustment” model is that several strong earthquakes in Sichuan, such as M8.0 Wenchuan earthquake, M7.0 Luzhou earthquake and M7.0 Jiuzhaigou earthquake, have changed the stress regime of the eastern margin of the Sichuan Basin by stress transference. Within the changed stress regime, a minor local stress adjustment has the possibility of making a notable earthquake event. In contract, the “fluid activation effect” model is mainly supported by the three evidences as follows: 1)the maximum principle stress axial azimuth is against the regional stress field, which reflects NWW-SEE direction thrusting type; 2)the Wujiang River crosscuts the pre-existing Wenfu normal fault and offers the fluid source; and 3)fractures along the Wenfu Fault formed by karst dissolution offer the important fluid flow channels.

Keyword: Wulong earthquake; focal mechanism; double-difference relocation; seismogenic fault
0 引言

据中国地震台网中心正式测定, 北京时间2017年11月23日17时43分, 重庆武隆区发生MS5.0地震, 震中位于29.39° N, 107.99° E, 震源深度为10km。据现场调查, 本次地震最高烈度为Ⅵ 度, 等震线长轴走向NNE, 长轴长20.9km, 短轴长11.8km(① http://www.cqdzj.gov.cn/content.php?id=6343& page=11。)。根据当地政府的灾害调查结果, 地震共造成武隆、 彭水两地10人受伤, 直接经济损失约4 079万元。此次地震是有仪器记录以来发生在重庆渝东南地区的最大地震, 历史上该地区仅发生过2次5级以上地震, 最强和最近的地震应追溯到距此次地震震中90km外的1856年黔江咸丰 6$\frac{1}{4}$级地震。在构造背景上, 该区域除了存在3条年代较老的深大基底断裂外, 还发育一系列中等规模断续分布的断层, 长度多为几十km, 年代学研究表明这些断层大多为早— 中更新世断层(丁仁杰等, 2004)。此次地震没有形成地表破裂带, 且震源区位于基岩山区, 从地质方面去研究发震构造存在一定的困难。另外, 现今GNSS数据、 中小震综合机制解反演、 现代地壳形变与地应力测量资料表明该区域处于以NWW-SEE向水平主压应力为主的现代构造应力场环境中(唐荣昌等, 1993; 刘玉亮, 2009; 陈涛等, 2018)。而本次MS5.0地震发生后, 不同机构给出的震源机制均显示该地震为一次NE走向的正断型张性破裂事件(表1), 这与区域压性应力场环境(丁仁杰等, 2004; 李峰等, 2013)不吻合。

表1 武隆MS5.0地震震源机制解 Table1 Focal mechanism solutions of the Wulong MS5.0 earthquake

武隆MS5.0地震在弱地震活动以及与其震源机制不一致的应力场背景下突然发生, 给该地区地震风险认识带来了新的挑战, 其发震构造及成因机制成为众研究者共同关注的问题。尽管此次地震震级不大, 造成的破坏也不严重, 但完整的地震序列为确定该地震的发震构造提供了重要的震源参数等信息, 为进一步研究该地区的发震断层、 地震成因机制以及渝东南地区未来地震形势的研判提供了参考信息。本文基于区域地震台网记录的武隆地震序列地震波形资料和震相到时数据, 采用CAP(Cut and Paste)方法反演了武隆MS5.0地震的震源机制和震源矩心深度; 同时, 利用双差定位方法对主震及其余震序列重新进行精确定位, 获得该序列较为准确的震源参数和余震活动的时空展布特征; 最后结合武隆MS5.0地震的区域构造环境和烈度分布等, 探讨了该地震序列可能的发震断层和成因机制。

1 区域地震构造环境

武隆MS5.0地震发生在七曜山-金佛山和彭水2条基底断裂所夹持的武隆凹陷束中, 位于芙蓉江逆断层与文复、 马武正断层交会的部位(图1)。该区域发育有一系列形成于燕山期的NE向、 NNE向高陡背斜带及断裂带(谭淋耘等, 2015; 方志雄等, 2016)(表2), 全新世区域构造活动较弱, 但继承性的局部褶皱、 断裂活动仍在继续, 局部表现为垂直或掀斜振荡, 这些活动的褶皱、 断裂也成为该地区孕震和发震的主要构造(丁仁杰等, 2004)。

图 1 研究区地震台站、 历史震中分布图(a)和区域断裂图(b)
震中附近的黄色椭圆为地震烈度Ⅵ 度区, 断裂信息见表2
Fig. 1 Seismic stations and historical earthquakes(a)and faults(b)in the study area.

表2 武隆震区及周边断裂性质 Table2 The fault properties in and around the Wulong MS5.0 earthquake area

本次地震震中西侧主要以压性或压扭性断裂为主, 多表现为发育于背斜两翼的逆冲推覆断裂, 如老场断裂(王赞军等, 2016)。而七曜山-金佛山基底断裂表现为高倾角、 NW倾向的正断层, 该断裂带自晚更新世以来没有明显活动, 仅在与长寿-遵义基底断裂交会处曾发生1854年南川 5$\frac{1}{2}$级地震。震区东侧断裂主要以发育在背斜轴部的正断层为主, 且断面倾向均为NW, 包括彭水基底断裂在地表显示的郁山正断层以及文复断层(四川省地质局107地质队, 1975, 黔江幅1:20万区域地质调查报告。)。彭水基底断裂性质相对复杂, 其顶部倾向SE, 下部倾向NW, 是一条总体倾向NW、 具有特殊地表(顶部)正断、 实质为逆冲性质特征的弧形断裂。该断裂的形成时代仍不清楚, 仅郁山正断层与NNW向横断层交切部位曾发生1855年彭水 4$\frac{3}{4}$级地震。郁山、 文复断层以及震中距10km范围内一系列倾向NW的次级正断层, 都是在燕山期晚期受NW-SE向挤压应力的逆时针方向张扭性应力作用而形成的(谢忱等, 2013)。除此之外, 震中北侧马武断层为倾向SE的正断性质, 但该断裂活动较弱, 仅在邻近地区发生了2013年7月18日石柱MS4.3地震。震中南侧发育一系列由NNE走向的三会冲、 芙蓉江及火石垭等逆断层组成的叠瓦式构造, 在第四纪可能有过活动, 这些老断层在外界环境影响下可能活化并重新活动, 如2003年武隆江口ML3.5地震与震中附近水库蓄水导致的芙蓉江断层活化有着重要关联(丁仁杰等, 2004; 齐文浩等, 2006)。

2 武隆MS5.0地震震源机制解

本文采用目前广泛应用的CAP波形反演方法(Zhao et al., 1994)获取武隆MS5.0主震的震源机制解。该方法将地震波形分为体波(Pnl)和面波2部分, 分别对截取的体波、 面波进行带通滤波, 计算理论波形与观测波形之间的误差函数, 然后利用网格搜索获取给定参数空间中误差函数值最小的最佳解, 解决了由于地震波速度模型不准确导致的不同震相到时差异产生的拟合偏差, 其不仅能有效减少对速度模型的依赖性(Tan et al., 2006; 龙锋等, 2010; Wang et al., 2016; 易桂喜等, 2016, 2017a, b), 还能降低对台站数量和方位角分布的要求, 提高反演的信噪比和精度。

本文采用王小龙等(2013)基于地震背景噪声层析成像技术反演得到的武隆地区精细地壳速度模型(波速比为1.79)(表3), 并基于台站方位角均匀分布的原则, 选择位于震中距250km范围内的重庆市、 湖北省和贵州省区域地震台网的9个台站(见图 1中蓝色三角形)高信噪比宽频带波形资料, 反演MS5.0主震的震源机制。计算中, 首先对挑选的三分量地震波形进行去均值和去除仪器响应处理, 旋转到大圆路径得到径向、 切向和垂向的速度记录, 再分成体波和面波, 并分别用0.05~0.2Hz和0.08~0.1Hz的频率范围进行带通滤波; 最后计算各台站的格林函数, 以5° 为网格步长搜索不同震源深度的最佳双力偶解。图 2显示了MS5.0主震的反演方差和震源机制解随深度的变化, 可见误差最小的最佳拟合震源深度为10km。图 3给出了此次地震在最佳拟合深度10km处的观测波形和理论波形拟合图, 绝大部分分量的观测波形和理论波形相关系数 > 80% , 二者具有较好的拟合关系。

表3 武隆地区速度模型 Table3 The seismic velocity model used in the location

图 2 武隆MS5.0地震震源机制解反演方差随不同深度的变化图Fig. 2 Variation of fitting error with depth during the focal mechanism inversion of the Wulong MS5.0 earthquake.

图 3 武隆MS5.0地震观测波形和理论波形拟合图
黑色实线为观测地震图, 红色实线为理论地震图; 波形左侧台站上方的数字为震中距(km), 下方的2行数字分别为理论地震图相对观测地震图的位移时间(s)和相关系数
Fig. 3 The fitting between observed and synthetic waveforms of the Wulong MS5.0 earthquake.

最终反演得到的主震震源机制解结果为: 节面Ⅰ 走向24° 、 倾角16° 、 滑动角-108° , 节面Ⅱ 走向223° 、 倾角75° 、 滑动角-85° , P轴方位140° 、 仰角60° , T轴方位309° 、 仰角30° 。结果显示此次地震主要是因NW-SE向的拉张作用而发生的。

此外, 我们还利用传统的P波初动方法, 选取震中距250km范围内具有清晰P波初动且方位角分布较均匀的23个台站求解主震震源机制。 计算结果显示, 节面Ⅰ 走向37° 、 倾角20° 、 滑动角-93° , 节面Ⅱ 走向220° 、 倾角70° 、 滑动角-89° , 与CAP波形反演结果基本一致。上述基于不同方法及不同机构给出的主震震源机制解(表1)均显示本次MS5.0地震为正断型。

3 地震序列重定位

据重庆地震台网资料统计, 震后1个月余震分布在NE向长轴约16km、 短轴约7km的区域内, 其深度优势分布范围非常模糊(图4)。本文选取了观测报告中地震序列震级在ML0.5以上, 至少有4个台站记录的201个地震事件, 采用双差定位方法(Waldhauser et al., 2000)进行重新定位。该方法采用震相走时差作为输入数据, 能够减小地壳速度模型的不确定性带来的部分误差, 适用于小台网地震观测中的大地震及其余震序列的重定位(杨智娴等, 2003; Waldhauser et al., 2008; 王小龙等, 2012; 张广伟等, 2014; 王清东等, 2015)。武隆地震发生在重庆地区地震监测能力相对较高的区域, 震中距100km范围内包含近20个测震台(图1), 密集分布的台站很好地包围了震源区, 使得该地区监测能力可达ML0.5, 同时近台资料也能提高地震的定位精度(张天中等, 2007; 黄媛等, 2008)。

图 4 武隆MS5.0地震序列初始空间分布图
a 震中及其邻区断层与余震分布图; b 平行于主震震源机制走向的震源深度剖面AB, 剖面比例为2:1; c 垂直于主震震源机制走向的震源深度剖面CD, 剖面比例为2:1。 黑色五角星表示主震。 WF 文复断层; HSY 火石垭断层
Fig. 4 Original spatial distribution of the MS5.0 Wulong earthquake sequence.

本文定位采用的速度模型见表3。在进行双差定位时, 选取震相对的数目≥ 8、 震源间距< 5km的地震走时进行配对组合, 采用LSQR反演算法进行5轮迭代, 共迭代25次, 其间不断调整迭代参数, 使得表征双差方程条件数的CND在40~60之间(一般为40~80比较合理)。迭代结束时双差数据的标准差(RMS)从最初的181ms减小为51ms, 平均震源改变量(DX, DY, DZ, DT)分别从530m、 662m、 9 274m和126ms减小为9m、 9m、 42m和1ms, 在接受误差范围之内。被定位到空气中的个数(AQ)也从4个减少到结束时的0个, 说明迭代过程趋于稳定及合理。

本文最终获得了武隆地震序列163个地震的重新定位结果, 重定位率达到84%。EW、 SN方向及深度上的平均定位误差分别为0.16km、 0.16km和0.26km。重定位结果显示, 武隆MS5.0主震震中位置为29.416° N, 107.983° E, 震源深度10.5km, 震中位置与中国地震台网正式目录的结果(29.39° N, 107.99° E)较为接近; 主震震源深度与CAP方法反演得到的主震矩心深度(10km)一致。另外, 序列中ML3.1最大余震的重定位震中位置为29.389° N, 107.957° E, 震源深度5km, 与中国地震台网正式目录的结果(29.386° N, 107.964° E, 震源深度为5km)几乎完全相同。

图 5展示了重新定位后的武隆地震序列空间分布特征。重定位后大部分余震分布在主震的SW侧, 震源深度分布范围较为集中, 优势深度为5~8km(图6), 且82%的余震位于主震的上部, 显示出武隆地震破裂向上扩展的过程。震源深度剖面表现为NW倾向、 宽约3km的余震条带, 通过最小二乘拟合得到的倾角约63° , 与主震震源机制解的节面Ⅱ 基本一致, 推测该余震条带对应发震断层, 大致勾勒出断层的产状(图5c)。余震主要分布在断层面的上盘, 即NW盘, 下盘余震较少, 这也符合正断型地震余震分布的特点。

图 5 重定位地震分布图
剖面和符号意义同图2, 虚线矩形框表示深度剖面的余震密集条带
Fig. 5 The distribution map of relocated earthquakes.

图 6 震源深度对比统计直方图Fig. 6 The statistical comparison of histograms of focal depth of original epicenter(a)and relocations(b).

为避免由于余震分布随时间的扩展或触发型余震对确定主震破裂尺度的干扰, 通常采用主震后前3d的早期余震分布估计破裂的尺度(龙锋等, 2006)。震后3d内的余震主要为SW向单侧破裂模式, 震源深度分布较离散, 余震多集中在主震附近, 余震密集区长约4.7km、 宽约2.4km, 平行于附近断裂的走向(图7)。根据正断层的矩震级与破裂长度和宽度的经验公式(Wells et al., 1994), 可以得出武隆MS5.0地震的理论破裂长2.3km、 宽2.6km, 与重定位后早期余震空间展布范围接近。震后4~6d余震呈SE向扩展模式, 且震源多集中于浅部, 主震深度以下几乎没有余震分布。

图 7 主震后不同时间段地震空间分布
剖面和符号意义同图2
Fig. 7 The spatial distribution of relocated earthquakes for different periods since the main shock.

4 讨论与结论
4.1 武隆MS5.0地震发震断层的判定

除一些特殊条件外, 中小地震一般很难形成地表破裂(唐茂云等, 2015)。尽管传统的槽探和钻探方法较为可靠, 但需要合适的地质条件和大量的经费支持, 因此存在一定的局限性。此外, 重庆武隆地区为基岩山区, 二叠纪、 三叠纪灰岩大量出露, 断层露头往往弥散性隐匿于山体裂缝与解理面中, 断续出露或不出露, 同时缺乏可靠的年代限制, 调查其长期活动性难度较大。利用地震序列的空间分布特征、 震源机制解以及烈度分布等推断发震断层(易桂喜等, 2015, 2017a, b)是目前较为有效的一种方法。

本文的重定位结果显示, 武隆MS5.0主震震中位于文复断层地表露头西侧约5km处, 震源深度10.5km。震后3d内余震震群主要分布于主震的SW侧, 表明在水平方向上地震序列向SW呈单侧破裂, 并成1条平行且靠近文复断层的密集条带, 趋势性较为明显。在深度剖面图中, 余震序列表现为倾向NW、 倾角为63° 的条带, 优势深度为5~8km, 与主震震源机制解节面Ⅱ 相吻合, 由此推测本次MS5.0地震的发震断层为一条走向223° 、 倾向NW且倾角为63° ~75° 的正断层, 且地震序列应主要位于地表可见断层的NW盘。烈度Ⅶ 度区长轴呈NNE向展布, 约$\frac{3}{4}$的主体面积位于主震的NW侧, 与正断层的烈度在断层倾向方向衰减较慢的认识一致(闻学泽, 2014), 指示发震断层倾向NW, 以正断活动为主。据地表露头考察获知, 文复断层是一条走向NE, 倾向NW, 倾角为60° ~73° 的正断层(四川省地质局107地质队, 1975, 黔江幅1:20万区域地质调查报告。) (丁仁杰等, 2004), 且武隆地震序列均位于其NW侧。综合上述因素, 推测文复断层为本次武隆MS5.0地震序列的发震构造(图8)。

图 8 武隆震区地质构造(a)和发震构造剖面(b)Fig. 8 Geological structure of the Wulong MS5.0 earthquake source area(a)and the seismic-tectonic profile(b).

4.2 武隆地震成因机制

GNSS数据显示重庆地区地壳整体向SE运动, 周边震源机制和小震综合面解显示主压应力场P轴的优势方位为NWW-SEE向, 倾角一般< 20° , 近水平(李峰等, 2013), 说明该地区整体应力状态仍以挤压为主。而本文的主震震源机制解结果显示该地震是以NW-SE向拉张作用下形成的正断层作用为主, 究其原因可能有2种解释:

(1)受外部应力扰动作用

区域强震不仅可以改变周边构造应力场的方向(张致伟等, 2015), 同时还可能通过应力传递的形式触发或延迟周围地震(Parsons et al., 2008; 万永革等, 2015; Han et al., 2017), 如2010年四川遂宁与重庆潼南交界的M5.0地震就可能为2008年汶川地震之后川东南应力调整的产物(张永久等, 2011)。武隆地区处于四川盆地东边界的渝东南地区, 属于弱地震背景, 表明该地区长期处于应变积累状态。 2008年汶川M8.0、 2013年芦山M7.0和2017年九寨沟M7.0地震不仅使得四川盆地周缘活动断裂上的应力发生变化(汪建军等, 2017; 徐晶等, 2017), 也可能对四川盆地东边界的渝东南地区的应力场产生扰动。Hill(2012)指出正断层在远程应力扰动作用下比其它断层类型更容易受到触发作用, 武隆地震则可能是在区域内多次强震扰动下, 因局部应力调整而发生的一次显著性事件。

(2)外部流体作用

流体作用对断层的影响与区域应力场的张性区或挤压区无关(马文涛等, 2013), 其对断层的影响主要通过2种方式, 一是流体的荷载作用, 二是孔隙压力改变。武隆地震震中位置与乌江水系最近的距离仅为2km(图3), 位于上游彭水电站与下游银盘电站水位段之间(共13级电站)。据现场调查, 银盘电站最大坝高78.5m, 库容3.2× 108m3, 电站正常蓄水位215m, 起蓄水位203m; 蓄水段自2011年开始蓄水, 历年水位变化不大, 一直在210~215m小幅波动, 经计算可知, 与地质体内原有应力状态相比, 该段蓄水荷载的影响很小。因此, 流体直接的荷载作用似乎不能解释此次地震的发震成因。对于孔隙压力的改变, Pandey等(2003)研究指出, 在流体扩散过程中, 约2~3bar的孔隙压力改变就足以诱发地震, 其中良好的渗透通道和局部构造条件是其最主要的2大条件。本次地震震中附近普遍发育二叠纪和三叠纪石灰岩地层, 沿文复、 芙蓉江等断层形成了一系列近乎垂直的 “ 天坑” 、 溶洞和地下暗河(夏凯生等, 2010)。而溶洞、 地下暗河等岩溶地貌是水体渗透扩散的有利通道(易立新等, 2000)。文复断层为张性断裂, 流体沿着这些通道向深部渗透, 在结构面上产生孔隙压力变化, 使得有效正压力减小, 岩石抗剪强度降低, 从而活化文复断层, 继承性地诱发此次地震。

4.3 结论及认识

(1)基于2017年11月23日重庆市武隆MS5.0地震的区域构造环境、 主震震源机制、 余震重定位结果和烈度分布特征等, 推测此次武隆地震的发震构造为文复断层。

(2)此次武隆MS5.0地震的震源机制解与区域构造应力场方向不一致, 经初步研究有2种成因解释: 一是区域内多次强震应力扰动导致局部区域应力场发生了偏转或调整; 二是震中区域岩溶地貌发育, 遍布溶洞和地下暗河, 加之震源附近乌江水系发达, 流体的深部渗透作用可能导致局部断层的活化。

致谢 Waldhauser为本研究提供了双差定位程序; 重庆地震台网中心提供了观测报告; 审稿人为本文提供了宝贵的意见和建议。在此一并表示感谢!

参考文献
[1] 陈涛, 张锐, 黄世源. 2018. 三峡重庆地区形变场及构造应力场研究[J]. 大地测量与地球动力学, 38(2): 129-132. doi: 1014075/j. jgg. 2018. 02. 004.
CHEN Tao, ZHANG Rui, HUANG Shi-yuan. 2018. Study on deformation field and tectonic stress field in Chongqing area of the Three Gorges[J]. Journal of Geodesy and Geodynamics, 38(2): 129-132(in Chinese). [本文引用:1]
[2] 丁仁杰, 李克昌. 2004. 重庆地震研究 [M]. 北京: 地震出版社.
DING Ren-jie, LI Ke-chang. 2004. Research of Earthquakes in Chongqing [M]. Seismological Press, Beijing(in Chinese). [本文引用:5]
[3] 方志雄, 何希鹏. 2016. 渝东南武隆向斜常压页岩气形成与演化[J]. 石油与天然气地质, 37(6): 819-827.
FANG Zhi-xiong, HE Xi-peng. 2016. Formation and evolution of normal pressure shale gas reservoir in Wulong syncline, southeast Chongqing, China[J]. Oil & Gas Geology, 37(6): 819-827(in Chinese). [本文引用:1]
[4] 黄媛, 吴建平, 张天中, . 2008. 汶川8. 0级大地震及其余震序列重定位研究[J]. 中国科学(D辑), 38(10): 1242-1249.
HUANG Yuan, WU Jian-ping, ZHANG Tian-zhong, et al. 2008. Relocation of the MS8. 0 Wenchuan earthquake and its aftershock sequence[J]. Science in China(Ser D), 51(12): 1703-1711. [本文引用:1]
[5] 李峰, 张效亮, 刘华国. 2013. 重庆市主要构造地震危险性评价[J]. 地震地质, 35(3): 518-531. doi: 103969/j. issn. 0253-4967. 2013. 03. 006.
LI Feng, ZHANG Xiao-liang, LIU Hua-guo. 2013. Seismic hazard analysis of the main faults in Chongqing urban area[J]. Seismology and Geology, 35(3): 518-531(in Chinese). [本文引用:2]
[6] 刘玉亮. 2009. 重庆黔江地区第四纪构造活动与地震地质研究 [D]. 广州: 中山大学.
LIU Yu-liang. 2009. Study on Quaternary tectonic movements and seismogeology in Qianjiang area of Chongqing Municipality [D]. Sun Yat-Sen University, Guangzhou(in Chinese). [本文引用:1]
[7] 龙锋, 闻学泽, 徐锡伟. 2006. 华北地区地震活断层的震级-破裂长度、 破裂面积的经验关系[J]. 地震地质, 28(4): 511-535.
LONG Feng, WEN Xue-ze, XU Xi-wei. 2006. Empirical relationships between magnitude and rupture length, and rupture area, for seismogenic active faults in North China[J]. Seismology and Geology, 28(4): 511-535(in Chinese). [本文引用:1]
[8] 龙锋, 张永久, 闻学泽, . 2010. 2008年8月30日攀枝花-会理6. 1级地震序列 ML≥4. 0事件的震源机制解[J]. 地球物理学报, 53(12): 2852-2860.
LONG Feng, ZHANG Yong-jiu, WEN Xue-ze, et al. 2010. Focal mechanism solutions of ML≥4. 0 events in the MS6. 1 Panzhihua-Huili earthquake sequence of Aug. 30, 2008[J]. Chinese Journal of Geophysics, 53(12): 2852-2860(in Chinese). [本文引用:1]
[9] 马文涛, 蔺永, 苑京立, . 2013. 水库诱发地震的震例比较与分析[J]. 地震地质, 35(4): 914-929. doi: 103969/j. issn. 0253-4967. 2013. 04. 020.
MA Wen-tao, LIN Yong, YUAN Jing-li, et al. 2013. Comparison and analysis on the basic features of reservoir-induced seismicity[J]. Seismology and Geology, 35(4): 914-929(in Chinese). [本文引用:1]
[10] 齐文浩, 薄景山, 李贞德. 2006. 江口水库诱发地震的有限元分析[J]. 地震工程与工程振动, 26(6): 7-10.
QI Wen-hao, BO Jing-shan, LI Zhen-de. 2006. Finite element analysis of Jiangkou reservoir induced earthquake[J]. Earthquakes Engineering and Engineering Vibration, 26(6): 7-10(in Chinese). [本文引用:1]
[11] 谭淋耘, 徐铫, 李大华, . 2015. 渝东南地区五峰组—龙马溪组页岩气成藏地质条件与有利区预测[J]. 地质学报, 89(7): 1308-1317.
TAN Lin-yun, XU Yao, LI Da-hua, et al. 2015. Geological condition of shale gas accumulation and favorable area prediction for the Wufeng-Longmaxi Formation in southeastern Chongqing[J]. Acta Geologica Sinica, 89(7): 1308-1317(in Chinese). [本文引用:1]
[12] 唐茂云, 刘静, 邵延秀, . 2015. 中小震级事件产生地表破裂的震例分析[J]. 地震地质, 37(4): 1193-1214. doi: 103969/j. issn. 0253-4967. 2015. 04. 020.
TANG Mao-yun, LIU-ZENG Jing, SHAO Yan-xiu, et al. 2015. Analysis about the minimum magnitude earthquake associated with surface ruptures[J]. Seismology and Geology, 37(4): 1193-1214(in Chinese). [本文引用:1]
[13] 唐荣昌, 韩渭滨. 1993. 四川活动断裂与地震 [M]. 北京: 地质出版社.
TANG Rong-chang, HAN Wei-bin. 1993. Research of Active Faults and Earthquakes in Sichuan [M]. Geological Publishing House, Beijing(in Chinese). [本文引用:1]
[14] 万永革, 盛书中, 李祥, . 2015. 2015年尼泊尔强震序列对中国大陆的应力影响[J]. 地球物理学报, 58(11): 4277-4286.
WAN Yong-ge, SHENG Shu-zhong, LI Xiang, et al. 2015. Stress influence of the 2015 Nepal earthquake sequence on Chinese mainland [J]. Chinese Journal of Geophysics, 58(11): 4277-4286(in Chinese). [本文引用:1]
[15] 汪建军, 许才军. 2017. 2017年 MW6. 5九寨沟地震激发的同震库仑应力变化及其对周边断层的影响[J]. 地球物理学报, 60(11): 4398-4420.
WANG Jian-jun, XU Cai-jun. 2017. Coseismic Coulomb stress changes associated with the 2017 MW6. 5 Jiuzhaigou earthquake(China)and its impacts on surrounding major faults[J]. Chinese Journal of Geophysics, 60(11): 4398-4420(in Chinese). [本文引用:1]
[16] 王清东, 朱良保, 苏有锦, . 2015. 2012年9月7日彝良地震及余震序列双差定位研究[J]. 地球物理学报, 58(9): 3205-3221. doi: 106038/cjg20150916.
WANG Qing-dong, ZHU Liang-bao, SU You-jin, et al. 2015. Double-difference relocation of the 7 September 2012 Yiliang earthquake and its aftershock sequence[J]. Chinese Journal of Geophysics, 58(9): 3205-3221(in Chinese). [本文引用:1]
[17] 王小龙, 马胜利, 郭志, . 2013. 利用地震背景噪声成像技术反演三峡库区及邻近地区地壳剪切波速度结构[J]. 地球物理学报, 56(12): 4113-4124. doi: 106038/cjg20131216.
WANG Xiao-long, MA Sheng-li, GUO Zhi, et al. 2013. S-wave velocity of the crust in the Three Gorges reservoir and the adjacent region inverted from seismic ambient noise tomography[J]. Chinese Journal of Geophysics, 56(12): 4113-4124(in Chinese). [本文引用:1]
[18] 王小龙, 马胜利, 雷兴林, . 2012. 重庆荣昌诱发地震区精细速度结构及2010年 ML5. 1地震序列精确定位[J]. 地震地质, 34(2): 348-358. doi: 103969/j. issn. 0253-4967. 2012. 02. 013.
WANG Xiao-long, MA Sheng-li, LEI Xing-lin, et al. 2012. Fine velocity structure and relocation of the 2010 ML5. 1 earthquake sequence in Rongchang gas field[J]. Seismology and Geology, 34(2): 348-358(in Chinese). [本文引用:1]
[19] 王赞军, 王宏超, 何宏林, . 2016. 七曜山-金佛山断裂带武隆土坎段最后活动时代研判厘定[J]. 四川地震, (1): 25-29.
WANG Zan-jun, WANG Hong-chao, HE Hong-lin, et al. 2016. Dating of the last active age for Wulong ridge section of Qiyaoshan-Jinfoshan faults[J]. Earthquake Research in Sichuan, (1): 25-29(in Chinese). [本文引用:1]
[20] 闻学泽. 2014. 内蒙古1979年五原6. 0级和1996年包头6. 4级地震的震源区构造[J]. 地震地质, 36(3): 586-597. doi: 103969/j. issn. 0253-4967. 2014. 03. 004.
WEN Xue-ze. 2014. Structures of source regions of the 1979 MS6. 0 Wuyuan earthquake and the 1996 MS6. 4 Baotou earthquake in Inner Mongolia, China[J]. Seismology and Geology, 36(3): 586-597(in Chinese). [本文引用:1]
[21] 夏凯生, 袁道先, 谢世友, . 2010. 乌江下游岩溶地貌形态特征初探: 以重庆武隆及其邻近地区为例[J]. 中国岩溶, 29(2): 196-204.
XIA Kai-sheng, YUAN Dao-xian, XIE Shi-you, et al. 2010. Preliminary study on morphology features of karst in the lower reaches of Wujiang River: A case in Wulong, Chongqing and the surrounding area[J]. Carsologica Sinica, 29(2): 196-204(in Chinese). [本文引用:1]
[22] 谢忱, 张金川, 李玉喜, . 2013. 渝东南渝科1井下寒武统富有机质页岩发育特征与含气量[J]. 石油与天然气地质, 34(1): 11-15.
XIE Chen, ZHANG Jin-chuan, LI Yu-xi, et al. 2013. Characteristics and gas content of the Lower Cambrian dark shale in Well Yuke 1, Southeast Chongqing[J]. Oil & Gas Geology, 34(1): 11-15(in Chinese). [本文引用:1]
[23] 徐晶, 邵志刚, 刘静, . 2017. 基于库仑应力变化分析巴颜喀拉地块东端的强震相互关系[J]. 地球物理学报, 60(10): 4056-4068.
XU Jing, SHAO Zhi-gang, LIU-ZENG Jing, et al. 2017. Analysis of interaction between great earthquakes in the eastern Bayan Har block based on changes of Coulomb stress[J]. Chinese Journal of Geophysics, 60(10): 4056-4068(in Chinese). [本文引用:1]
[24] 杨智娴, 陈运泰, 郑月军, . 2003. 双差地震定位法在中国中西部地区地震精确定位中的应用[J]. 中国科学(D辑), 33(S1): 129-134.
YANG Zhi-xian, CHEN Yun-tai, ZHENG Yue-jun, et al. 2003. The application of double difference earthquake location method in precise positioning of earthquakes in the central and western regions in China[J]. Science in China(Ser D), 33(S1): 129-134(in Chinese). [本文引用:1]
[25] 易桂喜, 龙锋, 梁明剑, . 2017a. 2017年8月8日九寨沟 M7. 0地震及余震震源机制解与发震构造分析[J]. 地球物理学报, 60(10): 4083-4097. doi: 106038/cjg20171033.
YI Gui-xi, LONG Feng, LIANG Ming-jian, et al. 2017a. Focal mechanism solution and seismogenic structure of the 8 August 2017 M7. 0 Jiuzhaigou earthquake and its aftershocks, northern Sichuan[J]. Chinese Journal of Geophysics, 60(10): 4083-4097(in Chinese). [本文引用:2]
[26] 易桂喜, 龙锋, 梁明剑, . 2017b. 2016年9月23日四川理塘 M4. 9和 M5. 1地震发震构造分析[J]. 地震地质, 39(5): 949-963. doi: 103969/j. issn. 0253-4967. 2017. 05. 006.
YI Gui-xi, LONG Feng, LIANG Ming-jian, et al. 2017b. Seismogenic structure of the M4. 9 and M5. 1 Litang earthquakes on 23 Sept 2016, southwestern China[J]. Seismology and Geology, 39(5): 949-963(in Chinese). [本文引用:2]
[27] 易桂喜, 龙锋, 闻学泽, . 2015. 2014年11月22日康定 M6. 3地震序列发震构造分析[J]. 地球物理学报, 58(4): 1205-1219. doi: 106038/cjg20150410.
YI Gui-xi, LONG Feng, WEN Xue-ze, et al. 2015. Seismogenic structure of the M6. 3 kangding earthquake sequence on 22 Nov. 2014, southwestern china[J]. Chinese Journal of Geophysics, 58(4): 1205-1219(in Chinese). [本文引用:1]
[28] 易桂喜, 龙锋, 赵敏, . 2016. 2014年10月1日越西 M5. 0地震震源机制与发震构造分析[J]. 地震地质, 38(4): 1124-1136. doi: 103969/j. issn. 0253-4967. 2016. 04. 025.
YI Gui-xi, LONG Feng, ZHAO Min, et al. 2016. Focal mechanism and seismogenic structure of the M5. 0 Yuexi earthquake on 1 Oct. 2014, southwestern China[J]. Seismology and Geology, 38(4): 1124-1136(in Chinese). [本文引用:1]
[29] 易立新, 车用太. 2000. 水库诱发地震及其水文地质条件和诱震机理[J]. 中国地质灾害与防治学报, 11(2): 46-50.
YI Li-xin, CHE Yong-tai. 2000. Reservoir-induced earthquake and its hydrogeological conditions and induced mechanism[J]. The Chinese Journal of Geological Hazard and Control, 11(2): 46-50(in Chinese). [本文引用:1]
[30] 张广伟, 雷建设, 梁姗姗, . 2014. 2014年8月3日云南鲁甸 MS6. 5地震序列重定位与震源机制研究[J]. 地球物理学报, 57(9): 3018-3027. doi: 106038/cjg20140926.
ZHANG Guang-wei, LEI Jian-she, LIANG Shan-shan, et al. 2014. Relocations and focal mechanism solutions of the 3 August 2014 Ludian, Yunnan MS6. 5 earthquake sequence[J]. Chinese Journal of Geophysics, 57(9): 3018-3027(in Chinese). [本文引用:1]
[31] 张天中, 武巴特尔, 黄媛, . 2007. 近台资料对近震相对定位算法的影响[J]. 地球物理学报, 50(4): 1123-1130.
ZHANG Tian-zhong, WU Bater, HUANG Yuan, et al. 2007. Effect of the data recorded at nearby station on earthquake relative location[J]. Chinese Journal of Geophysics, 50(4): 1123-1130(in Chinese). [本文引用:1]
[32] 张永久, 龙锋, 苏金蓉. 2011. 2010年1月31日遂宁5. 0级地震波形特征及其成因浅析[J]. 地震研究, 34(2): 119-125.
ZHANG Yong-jiu, LONG Feng, SU Jin-rong. 2011. Preliminary analysis of the waveform feature and the causes for the MS5. 0 Suining earthquake on Jan. 31, 2010[J]. Journal of Seismological Research, 34(2): 119-125(in Chinese). [本文引用:1]
[33] 张致伟, 阮祥, 王晓山, . 2015. 汶川、 芦山地震前后四川地区应力场时空演化[J]. 地震, 35(4): 136-146.
ZHANG Zhi-wei, RUAN Xiang, WANG Xiao-shan, et al. 2015. Spatial-temporal evolution of stress fields in Sichuan area before and after the 2008 Wenchuan and the 2013 Lushan earthquake[J]. Earthquake, 35(4): 136-146(in Chinese). [本文引用:1]
[34] Han L B, Peng Z G, Johnson C W, et al. 2017. Shallow microearthquakes near Chongqing, China triggered by the Rayleigh waves of the 2015 M7. 8 Gorkha, Nepal earthquake[J]. Earth and Planetary Science Letters, 479: 231-240. [本文引用:1]
[35] Hill D P. 2012. Dynamic stresses, Coulomb failure, and remote triggering—Corrected[J]. Bulletin of the Seismological Society of America, 102(6): 2313-2336. doi:10.1785/0120120085. [本文引用:1]
[36] Pand ey A P, Chadha R K. 2003. Surface loading and triggered earthquakes in the Koyna-Warna region, western India[J]. Physics of the Earth and Planetary Interiors, 139(3-4): 207-223. [本文引用:1]
[37] Parsons T, Ji C, Kirby E, et al. 2008. Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan Basin[J]. Nature, 454(7203): 509-510. [本文引用:1]
[38] Tan Y, Zhu L P, Helmberger D V, et al. 2006. Locating and modeling regional earthquakes with two stations[J]. Journal of Geophysical Research: Solid Earth, 111: B01306. doi:10.9/2005JB003775. [本文引用:1]
[39] Waldhauser F, Ellsworth W L. 2000. A double-difference earthquake location algorithm: Method and application to the Northern Hayward Fault, California[J]. Bulletin of the Seismological Society of America, 90(90): 1353-1368. [本文引用:1]
[40] Waldhauser F, Schaff D P. 2008. Large-scale relocation of two decades of Northern California seismicity using cross-correlation and double-difference methods[J]. Translated World Seismology, 113(B8): 4177-4183. [本文引用:1]
[41] Wang X L, Lei X L, Gou X B, et al. 2016. Focal mechanism solution of the February 23, 2014, ML3. 8 Rongchang earthquake with the CAP(Cut and Paste)method[J]. Earthquake Research in China, 30(2): 233-239. [本文引用:1]
[42] Wells D L, Coppersmith K J. 1994. New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement[J]. Bulletin of the Seismological Society of America, 84(4): 974-1002. [本文引用:1]
[43] Zhao L S, Helmberger D V. 1994. Source estimation from broadband regional seismograms[J]. Bulletin of the Seismological Society of America, 84(1): 91-104. [本文引用:1]