地球物理学报  2017, Vol. 60 Issue (5): 1735-1745   PDF    
川滇活动块体中-北部主要活动断裂带现今应力状态的分区特征
吴微微, 吴朋, 魏娅玲, 孙玮     
四川省地震局, 成都 610000
摘要: 恢复2009年1月1日—2015年12月31日间川滇活动块体中-北部1012次2.0≤ML≤5.0地震的震源谱,计算标量地震矩M0、震源尺度r和应力降Δσ等震源参数并拟合各参数间的定标关系,基于区域地震构造背景、活动断裂展布以及地震活动的成丛分布将研究区划分成四个统计单元,分别讨论各构造单元的应力分布特征、地震应力降随地点位置的变化以及应力-应变加载作用与区域变形的动力学过程的关联.结果显示:中小地震释放的应力降Δσ在0.1~10 MPa;标量地震矩M0与近震震级ML呈现较好线性关系(lgM0=0.92ML+10.46);应力降与地震大小的关系与Nuttli的增加应力降(ISD)模型比较吻合(lgΔσ=0.31 lgM0-3.92).震源应力降结果显示:① 金沙江断裂端部为低应力区,断裂单元整体滑动速率较高、强震活动极少,不具备强震发生的应力高度集中条件;在3条次级断裂构成的条带断裂结构中,理塘断裂上的应力-应变加载作用自北西向南东逐渐减弱,相对闭锁的北西段较其他部位更易积累应变. ② 鲜水河断裂带的地震应力降以康定为界南低北高,南段(康定—石棉)短期内难以积累较高应变,北段(甘孜—康定)应力水平较高,已发生的中强地震尚未能填充地震矩释放的亏空区,段落局部仍有较高的应力积累. ③ 安宁河—则木河断裂上高应力降地震事件集中,该单元的应变积累强、应力水平最高,地震危险性大. ④ 丽江—小金河断裂上不同震级地震的应力降特征并不相同,推测与当地复杂的构造背景有关,具体原因尚需深入探讨;木里地区应力背景较低,可能受当地构造环境的影响.研究表明,地震应力降随地点位置而系统变化,高应力降地震事件多发生在断裂与断裂的交汇部位,而断裂无闭锁条件、断裂以蠕滑为主且断面松弛、断裂端部为高温或破碎塑性变形带时,多以低应力降地震事件为主;与通常所认为的“走滑断裂不易积累应力”相反的是,鲜水河断裂带、安宁河—则木河断裂带均表现出较高的应力水平,其原因一方面可能是因为已发生的中强地震无论数量还是强度都尚不足以释放已经积累的能量,另一方面也许是在区域变形的复杂动力学过程中,当构造单元间阻碍断层运动和协助积累应力的作用占主导时,相同震级的地震会释放更多的应力.
关键词: 川滇活动地块      震源参数      定标关系      应力降      活动断裂     
Regional characteristics of stress state of main seismic active faults in mid-northern part of Sichuan-Yunnan block
WU Wei-Wei, WU Peng, WEI Ya-Ling, SUN Wei     
Sichuan Earthquake Administration, Chengdu 610000, China
Abstract: We restore the seismic source spectrums of 1012 earthquakes (2.0≤ML≤5.0) in the mid-northern part of Sichuan-Yunnan seismic block between January 1, 2009 and December 31, 2015, then calculate the source parameters (e.g., seismic moments M0, focal scale r and stress drop Δσ) and fit the calibration relationship between these parameters. Based on the regional seismic tectonic background, the distribution of active faults and seismicity, the study area is divided into four statistical units. For each unit the stress distribution characteristics, change of stress drop with location, correlation between the stress-strain loading and the dynamic process of regional deformation, are discussed respectively. The results show that seismic moments M0 are consistent with the magnitude-moment relation lgM0=0.92ML+10.46. The relationship between stress drop and magnitude is consistent with the result gained by Nuttli that intraplate earthquake follows the ISD model, with a statistical relationship lgΔσ=0.31 lgM0-3.92. Seismic source stress drop results show the following, (1) The stress at the end of the Jinshajiang fault is low, the overall sliding rate of the fault unit is high, and strong earthquake activity is very rare. In the fault belt consisting of three secondary faults, stress-strain loading deceases gradually from northwest to southeast along Litang fault, the northwest section which is relatively locked is more likely to accumulate strain than southeast section. (2) Stress drop of Xianshuihe fault zone is divided by Kangding, the southern section is low and northern section is high. Southern section (Kangding-Shimian) is difficult to accumulate higher strain in the short term, but in northern section (Garzê-Kangding), moderate and strong earthquakes have not filled the gaps of seismic moment release, there is still a high stress accumulation in partial section. (3) High stress-drop events were concentrated on Anninghe-Zemuhe fault zone, strain accumulation of this unit is strong, and stress level is the highest, earthquake risk is high. (4) On Lijiang-Xiaojinhe fault zone, stress drop characteristics of different magnitude earthquakes are not the same, which is related to complex tectonic setting, the specific reasons still need to be discussed deeply. Stress background in the Muli area is low and may be affected by the local tectonic environment. The study also shows that, (1) Stress drops display a systematic change with different faults and locations, high stress-drop events occur mostly on the fault intersection area. Faults without locking condition and mainly creeping, are mainly characterized by low stress drop. (2) Contrasting to what is commonly thought that "strike-slip faults are not easy to accumulate stress", Xianshuihe fault zone and Anninghe-Zemuhe fault zone all exhibit high stress levels, which may be due to that the magnitude and intensity of medium-strong earthquakes are not enough to release the accumulated energy. On the other hand, in the complex dynamics of regional deformation, when the tectonic unit blocking fault movement and its contribution to accumulation of stress play a key role, the earthquake of same magnitude will release higher stress drop.
Key words: Sichuan-Yunnan active faulted-block      Source parameter      Scaling relation      Stress drop      Active faults     
1 引言

大地震的震源特征,多反映较大区域范围的地震构造背景、震源受力方式和破裂过程结果,中小地震则更多地反映小区域范围应力状态的非均匀性变化.研究震源参数是在地震事件波形数据反演震源谱的基础上进行的,拟合地震震源谱与理论震源谱,即可得到相应的震源谱参数,进而估算标量地震矩M0、震源尺度r和应力降Δσ等震源参数.特别地,通过分析“应力降”参数的时空演变过程,可以帮助判断活动断裂现今应力应变的状态,区分出正处于相对高应力积累的段落,为探索强震发生的空间及时间的前兆信息提供重要依据 (Abercrombie,1995Allmann and Sheare, 2009).

目前,开展此项研究最多的是美国加州圣安德烈斯断层及其周围区域,Shearer等 (2006)利用加州中小地震震源参数的计算结果研究了南加州地震的定标率关系、应力降的时空变化过程及其与强震的相关性.Allmann和Sheare (2007)的研究指出2004年12月MS6.0地震前,震源区应力降显著高于断层上的其他地区,而在MS6.0地震发生后,震源区应力降明显下降.Hardebeck和Allegra (2009)认为,断层上的高应力降分布区域代表着这里的介质更强或者承受更高的外加剪应力,高应力降集中分布的区域也许是中强以上地震的潜在震源成核区.“十五”以来,由中国地震局监测预报司、中国地震局预测研究所牵头,与区域数字地震台网联合开展震源动力学参数研究,获得了中国大陆13个构造区域上中小地震震源特征的空间分布特征和时间变化过程,极大丰富了大陆块体内部地震震源定标率关系研究的样本.华卫等 (2009)的研究中,汶川地震序列中1070次ML≥3.0地震的地震矩与震级之间表现出较好线性关系,应力降和视应力的大小与震级大小有关,汶川余震序列地震应力降总体上随时间衰减,预示着主震发生后整个余震区应力降呈逐渐衰减的状态.赵翠萍等 (2011)认为中国大陆大部分地区的中小地震释放的应力降随震级增大而增大,不遵循常数应力降或分段模型,其中新疆天山中东段地区、青藏块体东北缘及龙门山断裂带东北缘是高应力降地震集中的区域.

川滇活动块体中-北部主要活动断裂包括金沙江断裂带及理塘断裂、鲜水河断裂带、安宁河—则木河断裂带以及丽江—小金河断裂带.作为川滇断块的一部分,研究区 (26°N—33°N,99°E—104°E) 仅过去300多年就至少发生过12次≥7级地震,并且区域主干、次级断裂的地震活动特征存在较大差异,7级以上强震多发生在鲜水河断裂上,安宁河断裂自1536年后、则木河断裂自1850年后、金沙江断裂自1870年后、理塘断裂自1948年后均未发生过≥7级的地震 (国家地震局震害防御司,1995中国地震局震害防御司,1999).历史震例表明长期的强震空段可能具有发生强震甚至大地震的活动构造背景,因此川滇活动块体中-北部中-长期大地震潜势及其现今应力状态变化已经成为关注重点.

为此,我们收集整理了四川、云南区域地震台网的数字波形资料 (国家测震台网数据备份中心,2007郑秀芬等,2009) 和中国地震台网中心的震相数据,自0.1~20 Hz的SH观测波形数据中扣除仪器、噪声、几何扩散、传播路径的介质衰减、台站场地响应等影响,恢复2009年1月1日—2015年12月31日间川滇活动块体中-北部1012次2.0≤ML≤5.0地震的震源谱;采用ω平方模型拟合震源谱,利用Brune的圆盘模型 (Brune,1970) 计算标量地震矩M0、震源尺度r和应力降Δσ等震源参数,拟合各参数间的定标关系;结合研究区地质构造背景、活动断裂展布以及地震活动在边界带的成丛分布,将研究区划分成四个基本统计单元 (金沙江断裂带及理塘断裂单元、鲜水河断裂带单元、安宁河—则木河断裂带单元和丽江—小金河断裂带单元),分别讨论各构造单元上的应力分布特征、地震应力降随地点位置的变化以及应力-应变加载作用与区域变形的动力学过程的关联,从而为进一步利用震源参数开展强震孕育的时空演化过程提供参考 (图 1).

图 1 川滇活动块体中-北部主要活动断裂、统计单元和震中分布 绿色虚线标示统计单元:Z1金沙江断裂带及理塘断裂单元,Z2鲜水河断裂带单元,Z3安宁河—则木河断裂带单元,Z4丽江—小金河断裂带单元;圆代表震中,灰色圆代表研究时段中的地震,蓝色圆代表计算应力降的地震;黑色实线代表主要断裂:F1巴塘断裂,F2理塘断裂,F3金沙江断裂带,F4鲜水河断裂带,F5安宁河断裂带,F6则木河断裂带,F7大凉山断裂带,F8丽江—小金河断裂,F9小江断裂;左下角为区域索引图,三角形代表区域地震台网宽频带固定台站;断层数据闻学泽等 (2003)的研究结果;地震目录来源:孙成民 (2010). Fig. 1 Distribution of main active faults, statistic units and earthquakes in the mid-northern part of Sichuan-Yunnan block Green dotted lines represent four statistic units; Circles represent earthquake epicenters; Black solid lines represent the major faults: F1 Batang fault, F2 Litang fault, F3 Jinshajiang fault, F4 Xianshuihe fault, F5 Anninghe fault, F6 Zemuhe fault, F7 Daliangshan fault, F8 Lijiang-Xiaojinhe fault, F9 Xiaojiang fault; Index map is shown on the down-left side, triangles are broadband seismic stations. Fault data from research results of Wen et al. (2003); Earthquake catalogs from Sun (2010).
2 计算结果 2.1 介质品质因子和台站场地响应

计算中小地震震源参数前,需先计算研究区的介质品质因子Q(f) 和所用地震台站的场地响应.选取互相衔接的三段地震波几何衰减模型,采用Atkinson和Mereu方法 (Atkinson and Mereu,1992) 计算介质品质因子Q(f),参照Moya方法 (Moya et al., 2000) 联合多台多地震数据求得地震台站的场地响应 (黄玉龙等,2003刘杰等,2003).由于四川地区介质横向不均匀性强烈,地震波穿过不同地块时会表现出不同衰减特征,因此将研究区分为川西高原 (图 1中Z1、Z2单元) 和攀枝花—西昌地区 (图 1中Z3、Z4单元,下文统一简称“攀西地区”),分别采用各自单元内的地震和台站计算介质衰减特征和地震台站场地响应.计算每个地震震源参数时,根据所用地震台站的位置,选择对应Q(f) 作为观测记录扣除非弹性衰减时采用的Q值.

拟合后的衰减模型分别是Q川西(f)=136.6f0.5813Q攀西(f)=101.9f0.6663,1~10 Hz频段上川西高原介质品质因子较攀西地区略高,10 Hz以上的高频段两者几乎一致.场地响应方面,川西台站场地增益在0.83~6.01倍,除REG、RTA地震台外,大部分台站增益明显;攀西台站场地增益在1.13~6.93倍,MGU、PGE和SMK地震台场地增益较大,YYU、YYC地震台场地放大倍数在1~2之间平稳变化,LGH地震台场地响应随频率增大而快速衰减 (吴微微等,2016).

2.2 定标关系

研究共使用地震波形数据万余条,按照一个地震至少被3个以上台站记录到、一个台站至少记录到3个以上地震及信噪比>1.5倍、震源谱拟合误差<0.3的要求,最终恢复2009年1月1日—2015年12月31日间川滇活动块体中-北部1012次2.0≤ML≤5.0地震的震源谱,获得这些地震的标量地震矩M0、震源尺度r和应力降Δσ等震源参数,拟合各震源参数之间的定标关系.

1012次地震的地震矩M0在1012~1016N·m (1.6≤MW≤4.6),拐角频率fc在0~20 Hz,大部分地震释放的应力降集中在0.1~10 MPa (图 2);矩震级MW2.0~4.8地震的震源尺度在80~1000 m间 (图 3a),面波震级MS与矩震级MW存在线性偏差 (图 3b),单对数坐标下标量地震矩M0与近震震级ML的线性关系 (图 3c) 和其他地震学家的结果基本一致 (Randall,1973Chen and Bai, 1992华卫等,2009赵翠萍等,2011):

(1)

图 2 拐角频率fc与标量地震矩M0、矩震级MW的关系 红线分别对应于应力降Δσ为0.1 MPa、1 MPa、10 MPa和100 MPa时M0fc的关系曲线 (ΔσM00.25) Fig. 2 Map showing the relationship between fc and M0, MW
图 3 矩震级MW与震源尺度 (a)、MS震级与矩震级MW(b)、近震震级ML与标量地震矩M0(c) 的关系 Fig. 3 MW versus source scale (a), MS versus MW(b) and ML versus M0(c)

中小地震释放的应力降与标量地震矩M0的定标关系拟合为

(2)

地震释放的应力降随着地震强度增大而增加 (图 4),与Nuttli的板内地震为增加应力降 (ISD) 模型 (Nuttli,1983) 的结果吻合.

图 4 应力降与标量地震矩M0(a)、应力降与矩震级MW(b) 的关系 Fig. 4 Stress drop versus M0(a), stress drop versus MW(b)
2.3 应力降计算结果

川滇活动块体中-北部438次2.0≤ML≤2.9地震平均应力降为1.70 MPa,554次3.0≤ML≤3.9地震平均应力降为2.86 MPa,20次4.0≤ML≤5.0地震平均应力降为4.52 MPa (下文将不同震级范围对应的平均应力降统一简称为“Δσ均值”).

考虑川滇活动块体中-北部地质构造背景、主干及次级活动断裂展布以及地震活动在边界带的成丛分布 (闻学泽,2000闻学泽和易桂喜,2003徐锡伟等, 2003a, 2005aWen et al., 2008),将研究区划分成四个基本统计单元分别讨论 (图 5):Z1金沙江断裂带及理塘断裂单元 (图 5a),Z2鲜水河断裂带单元 (图 5b),Z3安宁河—则木河断裂带单元 (图 5c),Z4丽江—小金河断裂带单元 (图 5d).

图 5 应力降空间分布图 (a) Z1金沙江断裂带及理塘断裂单元;(b) Z2鲜水河断裂带单元;(c) Z3安宁河—则木河断裂带单元;(d) Z4丽江—小金河断裂带单元;圆代表地震震中,不同圆的颜色随应力降增加由红变蓝;其他图例同图 1. Fig. 5 Spatial distribution of stress drop (a) Z1 Jinshajiang fault and Batang fault zone; (b) Z2 Xianshuihe fault zone; (c) Z3 Anninghe-Zemuhe fault zone; (d) Z4 Lijiang-Xiaojinhe fault zone. Circles represent earthquakes, circle color varies with the stress drop increases from red to blue. The legend same as Fig. 1.

比较各统计单元的震源应力降结果和相同震级范围的均值 (表 1),可以看出应力降大小在不同震级段的整体结果基本相似,应力降表现出随地点位置的系统变化,小区域上的应力分布并不均匀.

表 1 各统计单元的应力降计算结果 Table 1 Regional characteristics of stress drop in the mid-northern part of Sichuan-Yunnan block
3 讨论 3.1 金沙江断裂带及理塘断裂单元

研究时段中,金沙江断裂及理塘断裂单元2.0≤ML≤2.9地震应力降平均值为1.76 MPa,与1.70 MPa的Δσ均值相近;3.0≤ML≤3.9地震平均值为1.86 MPa,小于2.86 MPa的Δσ均值 (表 1图 5a);应力降相对较大的地震事件多分布在理塘断裂的北西段.整体上看,金沙江断裂及理塘断裂单元以频度不高的小震活动为主,目前处于中偏低应力状态.

由数字卫星影像和野外考察可知,金沙江断裂带由多条向东凸起的近SN向弧形次级逆断裂和NNE向巴塘主干断裂共同组成.据可查资料记载,弧形次级断层上几乎没有地震活动,巴塘断裂也仅发生过1870年MS71/2地震和1989年MS6.7震群.从活动构造上看,弧形次级逆断裂是由众多小断层和断层组合的复杂破碎带,应力不易集中,而巴塘主干断裂同时兼有右旋走滑和倾滑 (逆) 运动分量,全新世滑动速率并不低 (水平向6~7 mm·a-1,垂直 (逆) 速率>2~3 mm·a-1)(徐锡伟等,2005a).历史震例中,当断裂端部为低应力区 (高温、破碎的塑性变形带等)、断裂变形速率较高 (无闭锁条件) 时,断裂不具备强震发生的应力高度集中条件,会出现高滑动速率断裂的无强震活动现象 (强祖基和张立人,1983).结合我们的计算结果,认为金沙江断裂单元符合上述特征,地震危险性小.

理塘断裂是左旋走滑断裂,由三条次级断裂组成,地震破裂具有从北西向南东单方面时空迁移的特征,北西段同震错动会触发或加速相邻段上地震发生 (徐锡伟等,2005b).2009—2010年GNSS观测资料表明,理塘断裂周围地区存在显著的、具有一定空间规模的面应变象限分布,最大正应变和最大剪切应变均有闭锁迹象 (杨国华等,2012).应力降空间分布图像中 (图 5a),理塘断裂北西段地震应力降明显高于南东段,小震活动频度也较高.认为在3条次级断裂构成的条带断裂结构中,理塘断裂上的应力-应变加载作用沿自北西向南东逐渐减弱,相对闭锁的北西段较断裂其他部分更易积累应变,地震危险性也更大.

3.2 鲜水河断裂带单元

鲜水河断裂带是晚第四纪强烈活动的大型左旋走滑断裂,北与甘孜—玉树断裂斜列,南与安宁河断裂在康定一带相接,是中国大陆地壳运动变形最强烈的断裂带之一.据已知强震和大地震破裂背景资料 (Wen et al., 2008孙成民,2010) 以及小震活动时空分布显示:以康定为界,北侧甘孜—炉霍—道孚—康定以北 (图 5b中Ⅰ段) 地震频度高、强度大;康定以南—石棉以北 (图 5b中Ⅱ段) 小震、弱震活动频繁而集中.中小地震的应力降大小,以康定为界,呈南低北高的空间分布特征 (图 5b).

甘孜—康定段 (图 5b中Ⅰ段),2.0≤ML≤2.9地震应力降平均值为2.78 MPa,3.0≤ML≤3.9地震平均值为4.02 MPa,均高于Δσ均值 (1.70 MPa和2.86 MPa).显然,计算结果与通常所认为的“鲜水河走滑断裂断层面间的断层泥使应力降大大降低 (Engelder et al., 1975周蕙兰,1985)”并不吻合.

我们推测,在左旋走滑的鲜水河断裂上集中了较多的高应力降地震事件,可能是因为长久以来甘孜—康定段积累了巨大的能量,虽然研究时段中炉霍、道孚、康定相继发生过4次MS≥5.0地震,最近100多年也曾先后发生4次MS≥7.0强震 (Wen et al., 2008),但新近有强震的应变释放段与新近无强震段空间交替 (闻学泽,2000),在应力状态变化复杂的构造部位还存在着地震矩释放的亏空区,并且迄今已发生的中强地震无论数量还是强度都未能填充这些亏空,导致断裂的局部位置仍有较高的应力积累.甚至,随着时间推移,考虑强震的同震、震后以及震间效应的综合影响,个别敏感单元的应力水平可能会不同程度地增加.

康定—石棉段 (图 5b中Ⅱ段)2.0≤ML≤2.9地震应力降平均值为1.01 MPa,3.0≤ML≤3.9地震平均值为1.18 MPa,4.0≤ML≤5.0地震平均值为2.81 MPa,均低于Δσ均值 (1.70 MPa、2.86 MPa、4.52 MPa).该段落以相对集中的弱震活动为主,同时兼有左旋走滑和倾滑 (逆) 运动分量,水平和垂直滑动速率均较高 (9.6±1.7 mm·a-1、3.2±0.7 mm·a-1)(徐锡伟等,2003b),滑动断面松弛,处于低应力状态,短期内难以重新积累较高应力.

3.3 安宁河—则木河断裂带单元

安宁河—则木河断裂带单元上地震活动成丛分布,空间非均匀性强.2.0≤ML≤2.9地震应力降平均值为3.02 MPa,3.0≤ML≤3.9地震平均值为4.57 MPa,4.0≤ML≤5.0地震平均值为6.38 MPa (表 2),均显著高于Δσ均值 (1.70 MPa、2.86 MPa、4.52 MPa),是本文四个统计单元中应力状态最高的区域 (图 5c).

表 2 安宁河—则木河断裂带单元4.0≤ML≤5.0地震的基本参数 Table 2 4.0≤ML≤5.0 earthquake parameters of Anninghe-Zemuhe fault zone

单元北段 (图 5c中Ⅰ段) 是鲜水河断裂、安宁河断裂以及丽江—小金河断裂等多条活动断裂的交汇部位.以石棉为界,以北的康定—石棉段 (图 5b中Ⅱ段),弱震活动密集,应力水平较低;石棉南侧至冕宁,地震活动频度高,中小地震平均震级偏大,研究时段集中发生了5次ML>4.0地震 (表 2).震源应力降结果显示,相同震级的地震应力降表现出沿鲜水河断裂向安宁河断裂逐渐增大的趋势,段落整体应力水平中-偏高.

单元中段 (图 5c中Ⅱ段) 是安宁河断裂南段、则木河断裂北段以及大凉山断裂的交汇区,其中冕宁—西昌段和西昌—普格段分别是“安宁河地震空区”和“则木河空段”(闻学泽等,2008).研究时段中,随着普格、宁南、昭觉几次ML>4.0地震的发生以及大凉山断裂带小震活动的增强,空区面积明显收缩,这种现象被认为是安宁河断裂高度闭锁、应力应变积累在区域地震活动性上的响应,具有中期强震危险背景 (闻学泽,2003闻学泽等,2008).空间b值和地壳水平运动速度场同样表明,安宁河断裂呈现显著左旋剪应变积累,是川滇地区目前应变积累的主要地段 (江在森等,2005易桂喜等,2008谭毅培等,2015)

单元南段 (图 5c中Ⅲ段) 是则木河断裂南段和小江断裂北段的交汇区,集中了较多的高应力降地震事件,其中会东ML4.7和ML4.9地震的应力降值高达12.01 MPa和11.58 MPa (表 2).地震活动方面,2002年以来,沿则木河断裂—小江断裂带在较短时间内形成近南北向4级地震活动条带 (闻学泽,2003),随着占主导作用的大应力降事件——会东地震的发生,近年来该区地震活动明显增强.蒋长胜和庄建仓 (2010)的研究表明,1995年至今巧家地区b值明显下降,Δb负值变化持续时间较长.

综上所述,无论是震源参数研究获得的应力降空间分布,还是地震地质、地震活动性参数的研究结果,不同时间尺度、不同研究方法都揭示出安宁河—则木河断裂带单元现今正处于相对较高的应力水平,地震危险性大.

由此引出的问题是,为什么震源应力降在石棉一带会出现沿鲜水河断裂向安宁河断裂逐渐增大的趋势?应力-应变加载作用与活动断裂的动力学特征之间是否存在某种关联?我们注意到,鲜水河断裂带、安宁河—则木河断裂带、丽江—小金河断裂均以左旋走滑运动为主,走滑断层的错动方向是近水平的,断层扩展方向也是近水平的,因此走滑断层的错动可以促进断层沿走向的扩展 (马瑾等,2013),与此同时,应力-应变加载作用也沿断层走向进行传递和调整.在石棉一带,鲜水河断裂的左旋走滑运动被分解到安宁河、大凉山和则木河等断裂带上 (张培震,2008),在这个动力学过程中,安宁河断裂旁侧NE和NW向次级断层与近SN向安宁河断裂交汇形成“构造楔形体”,楔形断块尖端强力挤入并错断安宁河主断裂带,使主断裂的走滑活动受阻形成高度的应力集中 (黄润秋,2005),导致相同震级的地震所释放的应力降明显偏高.

此外,在安宁河—则木河断裂带单元上,高应力降地震事件多发生在断裂与断裂的交汇部位.在区域变形过程中,震源区和外围区里不同构造背景的主干和次级活动断裂局部交叉、错切,相互作用形成复杂的动力学系统,当阻碍断层运动和协助积累应力的作用占主导时,常导致同等震级下、高应力降地震事件的发生.

3.4 丽江—小金河断裂带单元

该单元的活动断裂系统以NE向丽江—小金河断裂带为主以及与之平行的次级断裂,此外还分布着木里弧形断裂、中甸—大具断裂和宁蒗断裂等.这些主干和次级活动断裂局部交叉、错切,导致该单元无论历史与现今都有频繁的地震活动.其中,丽江—小金河断裂切割川滇菱形块体并将其划分为川西北次级块体和滇中次级块体2个部分,宁蒗—盐源 (图 5d中Ⅰ段) 和木里地区 (图 5d中Ⅱ段) 正处于南、北两个次级块体的交界地带.考虑到木里地处锦屏水电站库区,其震源参数特征与构造地震并不相同,将其单独讨论.

研究时段中,宁蒗—盐源 (图 5d中Ⅰ段) 小震活动频繁,2.0≤ML≤2.9地震应力降平均值为1.96 MPa,大于1.70 MPa的Δσ均值;3.0≤ML≤3.9地震平均值为2.79 MPa,与2.86 MPa的Δσ均值相近;4.0≤ML≤5.0地震平均值为3.14 MPa,小于4.52 MPa的Δσ均值.显然该地区应力降大小在不同震级段上各有差异,推测与当地复杂的构造背景有关,具体原因尚需深入探讨.

2013年7月—12月木里地区集中发生了大量微震和有感小震群活动,蔡一川等 (2015)重新定位的震源分布显示,小震群的发生与木里地处锦屏水电站库区有关.统计结果显示,2.0≤ML≤2.9地震应力降平均值为1.07 MPa,3.0≤ML≤3.9地震平均值为1.35 MPa,4.0≤ML≤5.0地震平均值为3.22 MPa,均明显低于Δσ均值 (1.70 MPa、2.86 MPa、4.52 MPa).这种低应力现象被认为是,由于水库蓄水造成的地下介质孔隙压力增大或水的润滑作用,导致了较低构造应力下的地震发生 (Fehler and Phillips, 1991Abercrombie and Leary, 1993).

4 结论

我们利用四川、云南区域地震台网的数字波形资料和中国地震台网中心的震相数据,恢复2009—2015年川滇活动块体中-北部1012次2.0≤ML≤5.0地震的震源谱,计算了标量地震矩M0、震源尺度r和应力降Δσ等震源参数,并拟合它们的定标关系.结果显示:川滇活动块体中-北部中小地震释放的应力降在0.1~10 MPa;标量地震矩M0与近震震级ML呈现较好的线性关系,统计关系为lgM0=0.92ML+10.46;地震应力降与震级之间的关系和Nuttli的板内地震为增加应力降 (ISD) 模型较吻合,应力降与标量地震矩M0的定标关系拟合为lgΔσ=0.31lgM0-3.92.

震源应力降结果显示:四个构造单元的应力分布并不均匀;地震应力降随地点位置而系统变化,高应力降地震事件多发生在断裂与断裂的交汇部位,而断裂无闭锁条件、断裂以蠕滑为主且断面松弛、断裂端部为高温或破碎塑性变形带时,多以低应力降地震事件为主;与通常所认为的“走滑断裂不易积累应力”相反的是,鲜水河断裂带、安宁河—则木河断裂带均表现出较高的应力状态,其原因一方面可能是因为段落长久以来积累了巨大的能量,已发生的中强地震尚未能填充地震矩释放的亏空区,段落局部仍有较高的应力积累,另一方面可能是在区域变形的复杂动力学过程中,当阻碍断层运动和协助积累应力的作用占主导时,导致同等震级下、高应力降地震事件的发生.

致谢

文章承蒙杨建思研究员、闻学泽研究员、赵翠萍研究员、郑斯华研究员指点,衷心感谢;华卫研究员、王勤彩研究员指导相关计算;中国地震局地球物理研究所国家测震台网数据备份中心 (doi: 10.11998/SeisDmc/SN)、中国地震台网中心、四川、云南地震台网为本研究提供地震波形数据;两位匿名审稿人提出宝贵的修改意见,在此一并致谢.

参考文献
Abercrombie R E, Leary P. 1993. Source parameters of small earthquakes recorded at 2.5 km depth, Cajon Pass Southern California: implications for earthquake scaling. Geophysical Research Letters, 20(14): 1511-1514. DOI:10.1029/93GL00367
Abercrombie R E. 1995. Earthquake source scaling relationships from-1 to 5 ML using seismograms recorded at 2.5 km depth. Journal of Geophysical Research: Solid Earth, 100(B12): 24015-24036. DOI:10.1029/95JB02397
Allmann B P, Shearer P M. 2007. Spatial and temporal stress drop variations in small earthquakes near Parkfield, California. Journal of Geophysical Research: Solid Earth, 112(B4): B04305. DOI:10.1029/2006JB004395
Allmann B P, Shearer P M. 2009. Global variations of stress drop for moderate to large earthquakes. Journal of Geophysical Research: Solid Earth, 114(B1): B01310. DOI:10.1029/2008JB005821
Atkinson G M, Mereu R F. 1992. The shape of ground motion attenuation curves in southeastern Canada. Bull. Seism. Soc. Am., 82(5): 2014-2031.
Brune J N. 1970. Tectonic stress and the spectra of seismic shear waves from earthquakes. Journal of Geophysical Research, 75(26): 4997-5009. DOI:10.1029/JB075i026p04997
Cai Y C, Dai S G, Wei Y L. 2015. The active characteristics of small earthquake swarms in the Muli area of Sichuan province. Technology for Earthquake Disaster Prevention, 10(4): 933-940. DOI:10.11899/zzfy20150411
Chen P S, Bai T X. 1992. Quantification relations between the source parameters. Acta Seismologica Sinica, 5(3): 435-446. DOI:10.1007/BF02650539
Department of Earthquake Disaster Prevention, State Seismological Bureau. The Catalogue of Chinese Historical Strong Earthquakes. Beijing: Seismological Press, 1995: 1-514.
Department of Earthquake Disaster Prevention, China Seismological Bureau. The Catalogue of Modern Earth-quakes of China. Beijing: Chinese Science and Technology Press, 1999: 1-637.
Engelder J T, Logan J M, Handin J. 1975. The sliding characteristics of sandstone on quartz fault-gouge. Pure Appl. Geophys., 113(1): 69-86. DOI:10.1007/BF01592900
Fehler M C, Phillips S W. 1991. Simultaneous inversion for Q and source parameters of microearthquakes accompanying hydraulic fracturing in granitic rock. Bull. Seism. Soc. Am., 81(2): 553-575.
Hardebeck J L, Allegra A. 2009. Earthquake stress drops and inferred fault strength on the Hayward fault, East San Francisco Bay, California. Bull. Seism. Soc. Am., 99(3): 1801-1814. DOI:10.1785/0120080242
Hua W, Chen Z L, Zheng S H. 2009. A study on segmentation characteristics of aftershock source parameters of Wenchuan M8.0 earthquake in 2008. Chinese J. Geophys., 52(2): 365-371.
Huang R Q. 2005. Main characteristics of high rock slopes in southwestern China and their dynamic evolution. Advances in Earth Science, 20(3): 292-297.
Huang Y L, Zheng S H, Liu J, et al. 2003. Attenuation of ground motion and site response in Guangdong region. Chinese J. Geophys., 46(1): 54-61.
Jiang C S, Zhuang J C. 2010. Evaluation of background seismicity and potential source zones of strong earthquakes in the Sichuan-Yunan region based on the space-time ETAS model. Chinese J. Geophys., 53(2): 305-317. DOI:10.3969/j.issn.0001-5733.2010.02.008
Jiang Z S, Niu A F, Wang M, et al. 2005. Quantitative analysis for tectonic deformation on active rupture zones. Acta Seismologica Sinica, 27(6): 610-619.
Liu J, Zheng S H, Huang Y L. 2003. The inversion of non-elasticity coefficient, source parameters, site response using genetic algorithms. Acta Seismologica Sinica, 25(2): 211-218.
Ma J, Liu P X, Liu Y Z. 2013. Features of seismogenic process of the Longmenshan fault zone derived from analysis on the temporal-spatial evolution of earthquakes. Seismology and Geology, 35(3): 461-471.
Moya A, Aguirre J, Irikura K. 2000. Inversion of source parameters and site effects from strong ground motion records using genetic algorithms. Bull. Seism. Soc. Am., 90(4): 977-992. DOI:10.1785/0119990007
Nuttli O W. 1983. Average seismic source-parameter relations for mid-plate earthquakes. Bull. Seism. Soc. Am., 73(2): 519-535.
Qiang Z J, Zhang L R. 1983. Regionalization of quaternary active fractures and seismicity in China. Acta Geologica Sinica(4): 358-368.
Randall M J. 1973. The spectral theory of seismic sources. Bull. Seism. Soc. Am., 63(3): 1133-1144.
Shearer P M, Prieto G A, Hauksson E. 2006. Comprehensive analysis of earthquake source spectra in southern California. Journal of Geophysical Research: Solid Earth, 111(B6). DOI:10.1029/2005JB003979
State Earthquake Information Service-Data Management Centre. 2007. Waveform data of China National Seismic Network (in Chinese). Institute of Geophysics, China Earthquake Administration, doi: 10.11998/SeisDmc/SN, http://www.seisdmc.ac.cn.
Sun C M. Earthquake Records of Sichuan Province. Chengdu: Sichuan People's Publishing House, 2010.
Tan Y P, Xu F K, Cao J Q, et al. 2015. Missing earthquakes detection and temporal decay characteristics analysis of the aftershock sequence of the 2014 Ludian, Yunan, MS6.5 earthquake. Progress in Geophysics, 30(4): 1592-1596. DOI:10.6038/pg20150413
Wen X Z. 2000. Character of rupture segmentation of the Xianshuihe-Anninghe-Zemuhe Fault Zone, Western Sichuan. Seismology and Geology, 22(3): 239-249.
Wen X Z. 2003. Basic situation and work proposal of the Western Sichuan Strengthened Monitoring Region on Earthquake. Earthquake Research in Sichuan(2): 1-6.
Wen X Z, Yi G X. 2003. Re-zoning of statistic units of seismicity in Sichuan-Yunnan region. Journal of Seismological Research(z1): 1-9.
Wen X Z, Ma S L, Xu X W, et al. 2008. Historical pattern and behavior of earthquake ruptures along the eastern boundary of the Sichuan-Yunnan faulted-block, southwestern China. Phys. Earth. Planet. Inter., 168(1-2): 16-36. DOI:10.1016/j.pepi.2008.04.013
Wen X Z, Fan J, Yi G X, et al. 2008. A seismic gap on the Anninghe fault in western Sichuan, China. Science in China Series D: Earth Sciences, 51(10): 1375-1387. DOI:10.1007/s11430-008-0114-4
Wu W W, Su J R, Wei Y L, et al. 2016. Discussion on attenuation characteristics, site response and magnitude determination in Sichuan. Seismology and Geology, 38(4): 1005-1018.
Xu X W, Cheng G L, Yu G H, et al. 2003a. Tectonic and paleomagnetic evidence for the clockwise rotation of the Sichuan-Yunnan rhombic block. Seismology and Geology, 25(1): 61-70.
Xu X W, Wen X Z, Zheng R Z, et al. 2003b. Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region, China. Science in China Series D: Earth Sciences, 46(S2): 210-220.
Xu X W, Zhang P Z, Wen X Z, et al. 2005a. Features of active tectonics and recurrence behaviors of strong earthquakes in the western Sichuan province and its adjacent regions. Seismology and Geology, 27(3): 446-461.
Xu X W, Wen X Z, Yu G H, et al. 2005b. Average slip rate, earthquake rupturing segmentation and recurrence behavior on the Litang fault zone, western Sichuan province, China. Science in China Series D: Earth Sciences, 48(8): 1183-1196. DOI:10.1360/04yd0072
Yang G H, Yang B, Zhan W, et al. 2012. Spatial distribution of regional horizontal deformations before and after the Yushu and Wenchuan earthquakes. Earthquake, 32(2): 40-51.
Yi G X, Wen X Z, Su Y J. 2008. Study on the potential strong-earthquake risk for the eastern boundary of the Sichuan-Yunnan active faulted-block, China. Chinese J. Geophys., 51(6): 1719-1725.
Zhang P Z. 2008. The tectonic deformation, strain distribution and deep dynamic processes in the eastern margin of the Qinghai-Tibetan plateau. Scientia Sinica Terrae, 38(9): 1041-1056.
Zhao C P, Chen Z L, Hua W, et al. 2011. Study on source parameters of small to moderate earthquakes in the main seismic active regions, China mainland. Chinese J. Geophys., 54(6): 1478-1489. DOI:10.3969/j.issn.0001-5733.2011.06.007
Zheng X F, Ouyang B, Zhang D N, et al. 2009. Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake. Chinese J. Geophys., 52(5): 1412-1417. DOI:10.3969/j.issn.0001-5733.2009.05.031
Zhou H L. 1985. Some characteristics of source processes of large shallow strike-slip earthquakes. Acta Geophysica Sinica, 28(6): 579-587.
蔡一川, 戴仕贵, 魏娅玲. 2015. 四川木里小震群活动特征研究. 震灾防御技术, 10(4): 933–940. DOI:10.11899/zzfy20150411
国家测震台网数据备份中心. 2007. 国家测震台网地震波形数据. 中国地震局地球物理研究所, doi: 10.11998/SeisDmc/SN, notenoalianjie.
国家地震局震害防御司. 中国历史强震目录. 北京: 地震出版社, 1995: 1-514.
华卫, 陈章立, 郑斯华. 2009. 2008年汶川8.0级地震序列震源参数分段特征的研究. 地球物理学报, 52(2): 365–371.
黄润秋. 2005. 中国西南岩石高边坡的主要特征及其演化. 地球科学进展, 20(3): 292–297.
黄玉龙, 郑斯华, 刘杰, 等. 2003. 广东地区地震动衰减和场地响应的研究. 地球物理学报, 46(1): 54–61.
蒋长胜, 庄建仓. 2010. 基于时-空ETAS模型给出的川滇地区背景地震活动和强震潜在危险区. 地球物理学报, 53(2): 305–317. DOI:10.3969/j.issn.0001-5733.2010.02.008
江在森, 牛安福, 王敏, 等. 2005. 活动断裂带构造变形定量分析. 地震学报, 27(6): 610–619.
刘杰, 郑斯华, 黄玉龙. 2003. 利用遗传算法反演非弹性衰减系数、震源参数和场地响应. 地震学报, 25(2): 211–218.
马瑾, 刘培洵, 刘远征. 2013. 地震活动时空演化中看到的龙门山断裂带地震孕育的几个现象. 地震地质, 35(3): 461–471.
强祖基, 张立人. 1983. 中国第四纪活动断裂分区与地震活动性. 地质学报(4): 358–368.
孙成民. 四川地震全纪录. 成都: 四川人民出版社, 2010.
谭毅培, 徐甫坤, 曹井泉, 等. 2015. 2014年云南鲁甸MS6.5地震早期余震序列遗漏地震检测与衰减特征分析. 地球物理学进展, 30(4): 1592–1596. DOI:10.6038/pg20150413
闻学泽. 2000. 四川西部鲜水河—安宁河—则木河断裂带的地震破裂分段特征. 地震地质, 22(3): 239–249.
闻学泽. 2003. 川西地震强化监视区的基本情况与工作设想. 四川地震(2): 1–6.
闻学泽, 易桂喜. 2003. 川滇地区地震活动统计单元的新划分. 地震研究(z1): 1–9. DOI:10.3969/j.issn.1000-0666.2003.z1.002
闻学泽, 范军, 易桂喜, 等. 2008. 川西安宁河断裂上的地震空区. 中国科学D辑:地球科学, 38(7): 797–807.
吴微微, 苏金蓉, 魏娅玲, 等. 2016. 四川地区介质衰减、场地响应与震级测定的讨论. 地震地质, 38(4): 1005–1018.
徐锡伟, 程国良, 于贵华, 等. 2003a. 川滇菱形块体顺时针转动的构造学与古地磁学证据. 地震地质, 25(1): 61–70.
徐锡伟, 闻学泽, 郑荣章, 等. 2003b. 川滇地区活动块体最新构造变动样式及其动力来源. 中国科学D辑:地球科学, 33(增刊): 151–162.
徐锡伟, 张培震, 闻学泽, 等. 2005a. 川西及其邻近地区活动构造基本特征与强震复发模型. 地震地质, 27(3): 446–461.
徐锡伟, 闻学泽, 于贵华, 等. 2005b. 川西理塘断裂带平均滑动速率、地震破裂分段与复发特征. 中国科学D辑:地球科学, 35(6): 540–551.
杨国华, 杨博, 占伟, 等. 2012. 玉树和汶川地震前后区域水平形变的空间分布. 地震, 32(2): 40–51.
易桂喜, 闻学泽, 苏有锦. 2008. 川滇活动地块东边界强震危险性研究. 地球物理学报, 51(6): 1719–1725.
张培震. 2008. 青藏高原东缘川西地区的现今构造变形、应变分配与深部动力过程. 中国科学D辑:地球科学, 38(9): 1041–1056.
赵翠萍, 陈章立, 华卫, 等. 2011. 中国大陆主要地震活动区中小地震震源参数研究. 地球物理学报, 54(6): 1478–1489. DOI:10.3969/j.issn.0001-5733.2011.06.007
郑秀芬, 欧阳飚, 张东宁, 等. 2009. "国家数字测震台网数据备份中心"技术系统建设及其对汶川大地震研究的数据支撑. 地球物理学报, 52(5): 1412–1417. DOI:10.3969/j.issn.0001-5733.2009.05.031
中国地震局震害防御司. 中国近代地震目录. 北京: 中国科学技术出版社, 1999: 1-637.
周蕙兰. 1985. 浅源走滑大震震源过程的某些特征. 地球物理学报, 28(6): 579–587.