刘自凤 龙锋 彭关灵 赵小艳

摘要：利用1965年1月—2018年11月云南地震台网提供的月报目录和1999年1月—2018年10月中小地震震源机制解，对滇西北地区主要断裂的b值和震源机制一致性参数进行了空间扫描，结合历史强震空区，分析了滇西北地区的强震危险性。结果表明：①2000年1月以来程海断裂、龙蟠—乔后断裂及维西—乔后断裂南段b值显著降低，而中甸断裂地震空区两端、程海断裂及邻区则为低b值（b<07）主体区域。②震源机制一致性较好的区域分布于程海断裂中部、龙蟠—乔后断裂及维西—乔后断裂南段及邻区。分析认为，低b值和震源机制一致性好可能是区域应力水平增强所致，据此判断程海断裂、龙蟠—乔后断裂及维西—乔后断裂南段及邻区可能是滇西北地区未来发生强震的潜在危险区。

关键词：b值;震源机制一致性;应力张量方差;强震危险性;滇西北

中图分类号：P31572   文献标识码：A   文章编号：1000-0666（2019）03-0330-08

0 引言

强震活动往往发生在活动断裂带上具有高应力积累的凹凸体段或者闭锁段（Aki，1984;Wiemer，Wyss，1997;Wyss et al，2000）。地震一旦成核，破裂将会在高应力环境下扩展，并穿过凹凸体及其附近的非凹凸体的断裂段（Wyss et al，2004）。古登堡-里克特所提出的震级-频度关系（Gutenber，Richter，1944）中的b值与介质的非均匀度、有效剪应力等参数有关，并与应力呈反比关系，这一现象早已在岩石破裂实验（Scholz，1968）、与流体抽取有关的地震活动（Wyss，1973）以及地下矿山岩石破裂（Urbancic et al，1992）中被观测到。特定岩体内构造应力的大小与b值成反比，因此低b值可作为评判高有效剪应力的一个指标（Wyss et al，2000）。目前，利用b值的空间分布揭示和推断一个区域内活动断裂带不同段落的相对应力水平，已成为分析不同断裂段现今活动习性与强震（大震）潜在危险地段的重要手段之一（Wyss et al，2000;Wyss，Stefansson，2006;钱晓东，秦嘉政，2008;易桂喜等，2008，2010，2013;邵延秀等，2015;吴萍萍等，2015;毛燕等，2016;韩晓明等，2016;张广伟，2016）。

震源机制一致性参数（陈颙，1978;刁桂苓等，2004;王俊国，刁桂苓，2005;程万正等，2006;泽仁志玛等，2010）在震情跟踪工作中已经成为判断地震危险性的一个有效判据，基于震源机制解反演应力场的时空变化捕捉强震前兆信息的方法也被广泛运用。

滇西北地区处于青藏高原东南缘、川滇菱形块体西南部，受印度板块与欧亚板块碰撞作用（Dewey，Sun，1988;Rowley，1996）的影响，新构造活动显著，发育了众多受近南北向断裂控制的地堑或半地堑盆地，构成了“Z”字形的滇西北裂陷带（吴中海等，2008;黄小巾等，2014;罗睿洁等，2015）。该区域历史上M≥60地震频发，且大多发生在主要断裂及其附近区域，但M7专项工作组（2012）研究表明，在M≥60地震主体活动断裂上仍存在3个地震空区，分别为中甸断裂地震空区、小金河断裂地震空区和宾川地震空区。1996年2月3日丽江M70地震后，滇西北地区M≥60地震活动较弱，仅发生2组6级地震，最大震级MS62，可见该区自1996年丽江地震后强震进入弱活动状态，地震空区持续发展。2013年3月以来，滇西北地區ML≥40地震显著增强，且主要集中分布在维西—乔后断裂、程海断裂及其附近地区。2013年3月—2018年11月共计发生ML≥40地震20次，其中40～49级16次、50级以上4次，最大为2013年3月3日洱源55级地震。因此，探讨滇西北地区的强震危险性具有重要意义。本文结合滇西北地区历史强震空区，计算b值与震源机制一致性参数，分析该区域主要断裂和历史强震空区现今应力状态和强震危险性。

1 资料选取与预处理

11 地震目录及最小完整性震级

本文利用云南地震台网提供的1965年1月—2018年11月月报目录计算滇西北地区（250°～285°N，990°～1015°）b值。为避免统计结果 受资料选取的影响，余震的删除、小震群去丛和地震目录的完整性分析非常重要。本文在利用k-k法初步去除余震的基础上，人工进行了二次核对，最终获得0级以上地震94 148次。再采用EMR方法（基于G-R关系基础上的完整性震级范围法）计算研究区地震目录最小完整性震级MC值。由图1a可见，MC值随时间呈减小趋势，这与地震台网的扩展、地震分析方法的不断改进和监测能力的不断提高有关;2000年前后MC值大幅度降低，这是因为该时段云南地区台站实现了从模拟记录向数字化记录的转变（皇甫岗，李忠华，2010）。为了减弱MC的大幅变化和研究区内发生的1996年2月3日丽江M70地震对b值计算结果的影响，分2个时段分别计算b值，即1965年1月—1995年12月（Ⅰ时段）、2000年1月—2018年11月（Ⅱ时段）。图1b，c显示，MC在Ⅰ时段为22±003，在Ⅱ时段为13±003。

12 震源机制解

本文根据2002年立项的国家重大基础研究前期研究专项“新一代地震参数目录测定”、专著（徐彦，2013）、2016年立项的地震监测预报专项和中国地震局滇西地震预报实验场小孔径台网给出的相关结果整理了研究区1999年1月—2018年10月绝大部分中小地震震源机制参数。为了排除序列的余震对震源机制一致性参数的影响，剔除了余震的震源机制解，最终获得1 158次ML≥25地震的震源机制解，其中25～29级地震456次、30～39级地震576次、40～49级地震107次、50～59级地震15次、60～69级地震4次。

2 研究方法

21 b值计算

古登堡-里克特提出的震级与对应频度的对数线性关系为：

按照一定间隔将研究区进行网格化处理，采用最大似然法计算出ML≥MC地震的各统计单元的b值;再根据每一单元的b值得到相应区域的b值;最后基于应力与b值呈反比关系（Wyss et al，2000）分析滇西北主要断裂带和M7专项工作组（2012）提出的历史强震空区的地震危险性。易桂喜等（2008）研究指出，b>10的区域表明区域内断层具有以频繁小震活动为主要特征的蠕滑性质，未来发生强震的可能性相对较小，反之，则区域整体应力水平相对较高，而Wyss等（2004）所获得的凹凸体b<07，因此一般b<07的区域为低b值异常区。

22 震源机制一致性计算

已有研究提出多种表征震源机制一致性的参数，如单个地震震源机制解的P，T，B轴和构造应力场的3个正交应力主轴在三维空间的夹角之和（刁桂苓等，2004;赵英萍等，2004）、单个地震的滑动矢量与在平均剪切应力作用下产生的滑动矢量之间的夹角misfit角（Michael，1984，1987）、应力张量方差（Michael et al，1990;Lu et al，1997）和体波谱振幅相关系数（朱航等，2006;段云歌等，2018）等。本文采用Michael等（1990）反演研究区的應力张量方差，并用此表征震源机制一致性。Michael等（1990）的方法是按照剪应力方向与断层滑动方向相一致的原则提出的将非线性问题转化为线性反演应力张量的算法。计算过程中，由于无法区分震源机制解的断层面和辅助面，因此在反演应力张量时将2个节面等同看待，通过节面与假定应力张量的拟合函数关系选取其中拟合误差较小的节面作为震源机制断层面。基于研究区每个网格节点及其周围一定范围内的多个震源机制解，用每个震源机制解给激发断层运动的应力张量1个限制，通过若干个解的逼近则可以得到作用于该节点的应力张量解的方向，在应力张量结果的基础上求取应力张量方差和misfit角。

Michael等（1990）认为misfit角大小或应力张量方差高低均能表征震源机制一致性程度，其低值反映震源机制一致性好，应力场是均匀的;高值则反映震源机制一致性差，应力场是非均匀的。Lu等（1997）研究指出当应力张量方差大于02时，说明该区域的应力场在时空分布上都是非均匀的，震源机制一致性程度低，可解释为研究区区域应力水平不高;当应力张量方差小于01时，说明震源机制一致性程度高，且可以用统一的应力张量解释在该区域观测到的震源机制解，也可以理解为研究区域统一受控于较高的区域应力场。通过反演震源机制一致性参数的空间分布特征，可确定相对高应力区，再对相对高应力区的一致性参数进行时间扫描，根据变化形态可分析区域内的强震危险程度和紧迫性。综上，misfit角或应力张量方差的时空变化特征可以用来分析研究区未来强震的可能发生地点和时间。

3 结果分析

31 b值空间变化特征

本文研究区构造复杂，展布有龙蟠—乔后断裂、维西—乔后断裂、中甸断裂、小金河断裂、程海断裂和红河断裂等（图2a），历史强震几乎都发生在这些断裂带上，而宾川地震空区、小金河断裂地震空区和中甸断裂地震空区分别位于SN向的程海断裂中南段、NE向的小金河断裂带南段和NW向中甸断裂带上。因此将断裂带附近的中小地震密集区以003°的间隔进行网格化，以每个网格节点为圆心，在半径20 km的范围内进行地震统计并估算b值。计算条件要求每个网格的地震数量不少于30个，参与拟合的有效震级分档数不低于5档，不满足条件的网格在图中用空白显示，即没有b值计算结果。

分别以22，13为起始震级扫描计算研究区Ⅰ，Ⅱ时段的b值，空间分布如图2b，c所示。由图2b明显可见，Ⅰ时段的b值在中甸断裂地震空区、小金河断裂地震空区和程海断裂中段普遍低于10，表明该时段这些区域应力水平相对于宁蒗—盐源（b>10）和弥渡—祥云地区（b>10）要高，而在中甸断裂地震空区东南端和永胜—鹤庆一带的低b值异常尤为显著（b<07）;图2c所示，Ⅱ时段b值在中甸断裂地震空区、程海断裂中北段和剑川—洱源—漾濞地区低于10，而中甸断裂地震空区南端、宾川地震空区北端与永胜地区、程海断裂北段与小金河断裂地震空区东南角交汇区域b值低于07，几个局部地区b值甚至低于06。

对比Ⅱ时段与Ⅰ时段b值空间分布可看出，中甸地区低b值（b<07）区域有所扩展，永胜—鹤庆地区低b值（b<07）区域明显向程海断裂方向收缩，在程海断裂北段新增一个低b值（b=05～07）区域。而丽江地区b值从Ⅰ时段低于10显著升高为11～15，宾川地震空区南部的b   值也从Ⅰ时段的08～10显著升高为13以上。因此，统一以22（Ⅰ时段的最小完整性震级）为起始震级计算Ⅱ时段相对于Ⅰ时段的b值变化量Δb（图2d）。图2d显示，Ⅱ时段宁蒗南部—永胜—宾川一带和鹤庆—剑川—洱源—漾濞一带b值降低（Δb<0）。由此可见，2000年以来，b值显著降低的区域几乎沿着程海断裂、龙蟠—乔后断裂和维西—乔后断裂南段展布，这与2000年以来滇西地区4级以上地震空间分布格局相吻合。

32 震源机制一致性参数空间变化特征

本文收集的1999年1月—2018年10月滇西北地区的1 158次ML≥25地震震源机制空间分布如图3a所示，为了便于与b值计算结果进行对比，震源机制一致性参数扫描区域与b值扫描区域一致。将拟扫描区域按照005°×005°进行网格化，选取每个网格节点及其周围至少15个地震的震源机制解反演每个节点的应力张量方差。研究区应力张量方差空间分布特征如图3b所示，从图中可见，永胜—鹤庆—剑川—洱源—宾川一带、小金河断裂

4 结论与讨论

本文利用1965年1月—2018年11月云南地震台网提供的月报目录和1999年1月—2018年10月中小地震震源机制解，分析了滇西北主要断裂带和3个历史地震空区的强震危险性，主要得出以下结论：

（1）2000年1月以来，研究区的低b值（<07）主体区域为中甸断裂地震空区南端、宾川地震空区北端与永胜地区、程海断裂中北段与小金河断裂地震空区交汇地区，上述区域应力水平较高，具备发生强震的条件，而丽江地区b值较高，发生强震的可能性不大。

（2）对比Ⅱ时段与Ⅰ时段b值空间分布可知：中甸地区低b值区域有所扩展，这可能与Ⅱ时段MC降低，该地区参与计算的地震相对增多有关;永胜—鹤庆地区低b值区域往程海断裂方向大幅收缩，且在程海断裂北段新增一个低b值区域，而在其南段b值升高且空间上明显出现扩展，这表明程海断裂中北段的强震危险性在不断增加;丽江中东部和鹤庆北部地区b值显著升高，这可能与1996年丽江7级地震能量强释放有一定关系。Δb计算结果显示，Ⅱ时段永胜—宾川一带和鹤庆—剑川—洱源—漾濞一带b值相对于Ⅰ时段降低。由此可见2000年以来，b值显著降低的区域几乎沿着程海断裂、龙蟠—乔后断裂和维西—乔后断裂南段分布，这与2000年以来滇西地区4级以上地震空间分布格局基本一致。

（3）应力张量方差低值主体区域为永胜仁和地区、宾川—鹤庆—剑川—洱源地区，其值为010～015。这表明上述地区應力是相对均匀的，震源机制一致性较好，应力水平较高，发生强震的可能性较大。

综上，滇西北地区b值和应力张量方差的空间扫描结果显示程海断裂、龙蟠—乔后断裂和维西—乔后断裂南段及邻区低值变化特征相对显著，这可能是区域应力水平增强的表现。因此分析认为，未来滇西北地区的强震发生在程海断裂、龙蟠—乔后断裂和维西—乔后断裂南段及其邻区的危险性较大。

根据b值和震源机制一致性参数结果可判别出强震的相对危险地点，但在时间尺度上属于中-长期类，如果要在时间判断方面给出更明确的结论，需要结合其它地震学参数和地球物理观测资料开展进一步的研究后再进行综合判定。

参考文献：

陈阳，吕悦军，谢卓娟，等2013地震活动性参数b值的研究[C]北京：地壳构造与地壳应力所，38-47

陈颙1978用震源机制一致性作为描述地震活动性的新参数[J].地球物理学报，21（2）：142-159

程万正，阮祥，张永久2006川滇次级地块震源机制类型与一致性参数[J].地震学报，28（6）：561-573

刁桂苓，赵英萍，啜永清，等2004大同晚期强余震前震源机制解的一致性特征[J].内陆地震，18（4）：202-206

段云歌，陈天长，苏金蓉，等20182013年芦山MS70地震余震序列震源机制一致性的时空变化[J].地震研究，41（2）：252-257

韩晓明，张文韬，王树波2016河套地震带的b值时空变化特征分析[J].中国地震，32（3）：522-532

皇甫岗，李忠华201020世纪云南地区地震记录完整性评价[J].地震研究，33（1）：1-6

黄小巾，吴中海，李家存，等2014滇西北裂陷带的构造地貌特征与第四纪构造活动性[J].地质通报，33（4）：578-593

李永莉，蔡静观，曹刻2002云南地区强震活动过程中的调制比b值[J].地震研究，25（1）：25-30

罗睿洁，吴中海，黄小龙，等2015滇西北宾川地区主要活动断裂及其活动构造体系[J].地质通报，34（1）：155-170

毛燕，刘自凤，叶建庆2016小江断裂带强震危险性分析[J].地震研究，39（2）：213-217

钱晓东，秦嘉政2008小江断裂带及周边地区强震危险性分析[J]地震研究，31（4）：354-361

邵延秀，袁道阳，梁明剑2015滇西南地区龙陵—澜沧断裂带地震危险性评价[J].地震学报，37（6）：1011-1023

王辉，曹建玲，荆燕，等2012川滇地区强震活动前b值的时空分布特征[J].地震地质，34（3）：531-543

王俊国，刁桂苓2005千岛岛弧大震前哈佛大学矩心矩张量（CMT）解一致性的预测意义[J].地震学报，27（2）：178-183

吴萍萍，李振，叶庆东2015郯庐断裂带南段及邻区地震b值的空间分布特征[J].中国地震，31（2）：372-381

吴中海，张永双，胡道功，等2008滇西北哈巴玉龙雪山东麓断裂的晚第四纪正断层作用及其动力学机制探讨[J].中国科学：地球科学，38（11）：1361-1375

徐彦2013云南地区ML30以上中小地震震源机制解汇编[M].昆明：云南科技出版社，168-186

易桂喜，闻学泽，苏有锦2008川滇活动地块东边界强震危险性研究[J].地球物理学报，51（6）：1719-1725

易桂喜，闻学泽，辛华等2013龙门山断裂带南段应力状态与强震危险性研究[J].地球物理学报，56（4）：1112-1120

易桂喜，闻学泽，张致伟2010川南马边地区强震危险性分析[J].地震地质，32（2）：282-293

泽仁志玛，刁桂苓，李志雄，等2010大震前显示的地震震源机制趋于一致的变化[J].地震，30（1）：108-114

张广伟2016云南地区地震的重新定位及b值研究[J].中国地震，32（1）：54-62

赵英萍，刁桂苓，高景春，等2004张北强余震前震源机制解的一致性特征[J].华北地震学报，22（1）：1-4

朱航，刘杰，陈天长2006采用体波谱振幅相关系数方法研究地震序列的震源机制变化过程[J].地震，26（2）：1-11

M7专项工作组2012中国大陆大地震中-长期危险性研究[M].北京：地震出版社，234-235