地球物理学报  2016, Vol. 59 Issue (9): 3269-3279   PDF    
鲁甸MS6.5地震前地震活动性分析及川滇东边界中部未来地震危险性研究
刘月 , 田勤俭 , 吕晓健 , 邵志刚 , 李文巧     
中国地震局地震预测研究所, 中国地震局地震预测重点实验室, 北京 100036
摘要: 本文努力探索地震活动性定量化前兆的检测方法,针对Region-Time-Length(RTL)算法在空间异常分析中的不足,根据新提出的综合衡量区域地震活动水平的物理参数——RTL面积分(IRTL),定量分析了2014年鲁甸MS6.5地震前地震活动时空变化,并与2012年彝良MS 5.7、5.6地震作了对比,研究结果显示这些地震前都检测到地震活动增强.彝良地震前地震活动增强异常主要分布在则木河断裂和昭通-莲峰断裂带西段的交汇区,鲁甸地震前异常主要分布在昭通-莲峰断裂带、马边断裂带及周边地区.以上地震活动增强异常区的IRTL随时间呈现由小变大,至峰值后回落的特征.彝良地震和鲁甸地震均发生于IRTL峰值后数月,这表明IRTL峰值对地震发生可能有一定的指示意义.鲁甸地震后,对周边地区地震活动跟踪研究检测到地震活动增强现象,且自2015年IRTL呈上升趋势,目前IRTL值已超过鲁甸地震前的峰值水平,也许更强的地震正在孕育,所以川滇交界东边界中部仍存在发生强震的可能.以上研究为我们认识鲁甸地震的孕震过程和地震前兆研究提供了新的认识.
关键词: RTL算法      RTL面积分      鲁甸地震      地震危险性分析     
Analysis of the seismicity changes prior to the MS6.5 Ludian earthquake and study on the potential strong-earthquake risk for the middle part of eastern boundary of Sichuan-Yunnan
LIU Yue, TIAN Qin-Jian, LÜ Xiao-Jian, SHAO Zhi-Gang, LI Wen-Qiao     
Institute of Earthquake Science, China Earthquake Administration, Beijing 100036, China
Abstract: For lacking physical parameters to analyze space anomalies of seismic activity which detected by Region-Time-Length (RTL) algorithm, an effort to explore the precursory parameter for quantifying seismic activity in space was conducted. The RTL algorithm and the new physical parameter-surface integral of RTL (IRTL), which could synthetically evaluate the space-time changes of seismicity in a region, were applied to analyze the seismicity changes prior to the 2014 Ludian MS 6.5 earthquake and the 2012 Yiliang MS5.7, MS5.6 earthquakes. An enhancement of seismic activity was detected prior to the above earthquakes. Space anomalies mainly distributed in Zemuhe fault and west part of Zhaotong-Lianfeng fault prior to Yiliang earthquakes. Before Ludian earthquake, anomalies were detected in Zhaotong-Lianfeng fault, Mabian fault and surrounding areas. The parameter IRTL of anomalous region changed from small to large, after reaching its peak value, then reduced. Several months after the peak points, Yiliang and Ludian earthquakes occurred. IRTL may be an indicator for earthquake occurrence. After Ludian earthquake, an enhancement of seismic activity was detected in the middle part of eastern boundary of Sichuan-Yunnan, and IRTL was increasing since January 2015. Now, the value of IRTL was larger than the peak value prior to Ludian earthquake. Perhaps a stronger earthquake may be preparing in the middle part of the eastern boundary of Sichuan-Yunnan in the future. The above studies may provide some new evidence for better understanding the seismogenic process of Ludian earthquake and earthquake precursors..
Key words: Region-Time-Length (RTL)      Surface integral of RTL      Ludian Earthquake      Seismic risk analysis     
1 引言

2014年8月3日云南省鲁甸县发生了MS6.5地震(简称鲁甸地震),震中位于27.1°N,103.3°E,震源深度12 km(据中国地震台网中心),震源机制为走滑型(张广伟等,2014; 王未来等,2014).此次地震虽为中等强度,却造成重大的人员伤亡和财产损失.研究者在发震构造、孕震机制等方面展开了广泛探讨(张广伟等,2014; 王未来等,2014; 徐锡伟等,2014; 刘成利等,2014; 何骁慧等,2015; 许力生等,2014; 张勇等,2015; 徐涛等,2014; 陈石等,2014).地震发生在青藏高原东缘与华南地块相互作用的边界变形带上,位于小江断裂带东侧北东向的昭通—莲峰断裂带内,破裂过程具有一定的复杂性,发震构造是单向破裂还是共轭破裂存在争议(徐锡伟等,2014; 张勇等,2015).徐涛等(2014)通过上地壳速度结构研究显示震中及周边地区强度较低,小江断裂带两侧应变速率很高,因此小江断裂带和鲁甸—昭通断裂带具有发生强震的力学机制.三维重力反演研究显示震源区南北两侧块体运动受到昭通断裂不同程度阻隔,于震前5个月在震中附近出现重力异常,持续了4个月(陈石等,2014),这对识别与局部重力变化有关的地震前兆有重要意义.

除了重力异常(陈石等,2014),在众多前兆中,地震活动性研究在探索地震发生规律、判断地震发生危险性中扮演着重要角色,根据地震活动图像演化预测地震三要素取得了一些经验成果(Aki,1965; Gutenberg and Richter,1944; Mogi,1979; 张国民等,2005; 刘蒲雄和吕晓健,2012; Gardonio et al.,2015; Ishibe et al.,2015; 闻学泽等,2013).但在地震预测实践中,对空区和条带的确定及变化特征分析又存在一些不确定性和复杂的问题.Region-Time-Length(RTL)算法(Sobolev and Tyupkin,1997)实现了对地震活动水平的定量分析,在一定程度上降低了地震活动图像识别中的不确定性.该方法在勘察加地区MS>7.0、阪神MS7.2、Nemuro半岛MW6.8、汶川MS8.0、芦山MS7.0等多个强震,以及意大利、我国华北和甘肃等地区的中强震前地震活动变化回顾性检验中取得较好效果(Sobolev and Tyupkin,1999; Huang et al.,20012002; Huang,2008; Chen and Wu,2006; Gentili,2010; 蒋海昆等,2004a2004b; 梅秀苹和刘小凤,2013; 刘月等,2016),但在地震发生危险性实践分析中却应用不多.究其原因是在未知震中位置和发震时间的前提下,尚不能通过较大空间范围RTL分布随时间的演化正向判断地震危险进程及可能的发震地点等.一方面是对震前RTL空间分布和异常演化特征与地震发生位置和时间关系的认识不够,另一方面是缺少定量衡量区域地震活动综合水平的物理参数,不能定量建立空间地震活动异常与地震发生时间之间的关系.鉴于此,本文根据RTL方法检验了鲁甸地震前地震活动时空演化特征,并与2012年9月7日彝良MS5.7、MS5.6地震对比,根据新提出的综合衡量区域地震活动水平的参数——RTL面积分(IRTL)(刘月等,2016),分析震前空间区域的地震活动演化与地震发生时间之间的关系,并进一步探讨川滇交界东边界中部地区未来强震危险性.

2 相关理论及基础资料 2.1 RTL算法介绍

Region-Time-Length(RTL)方法由Sobolev和Tyupkin(1997)提出,该方法以地震目录为研究资料,定量检测地震活动偏离背景水平的异常程度.对于给定点(x,y)t时刻的RTL值可定义为震中距、时间及破裂尺度减去背景值后分别除以各自标准差的乘积.在t时刻之前发生的地震与点(x,y)的震中距函数R(x,y,t)、时间函数T(x,y,t)、破裂尺度函数L(x,y,t)以及RTL(x,y,t)分别表示为:

(1)

地震的破裂尺度l可由(2)式得到(郭增建和秦保燕,1991)

(2)

ri表示第i个地震与研究点(x,y)的距离,r0为特征距离;tili分别表示第i个地震的发震时间和破裂尺度,t0为特征时间.Rbk(x,y,t)Tbk(x,y,t)Lbk(x,y,t)分别代表震中距、时间和破裂尺度的背景值.对震级、距离及时间,引入了如下约束:

(3)

式中Mi为第i个地震的震级,震级下限Mmin通常取最小完备性震级,上限小于研究震例的震级.满足以上约束条件的地震数目记为n.从理论上讲,距离点(x,y)越远,与t的时间间隔越长,且破裂尺度越小的地震对(x,y,t)的影响越小;距离点(x,y)越近,与t的时间间隔越短,且破裂尺度越大的地震对(x,y,t)的影响也就越大.

RTL(x,y,t)<0代表t时刻的地震活动水平低于t之前的背景水平,可看作地震活动平静;RTL(x,y,t)>0意味着t时刻的地震活动水平高于t之前的背景水平,即地震活动增强.在实际分析中为了排除正常扰动,通常把RTL(x,y,t)<-1看作地震活动平静,把RTL(x,y,t)>1看作地震活动增强.

2.2 RTL面积分——IRTL

RTL方法对研究点(震中)在t时刻前的地震活动变化的研究取得了较好的效果,但对空间异常分析存在明显不足,为了综合定量分析区域地震活动水平以及与地震发生危险性的关系,我们提出了RTL面积分——IRTL(刘月等,2016):

(4)

研究区域地震活动平静时,取计算阈值RTL<-1;地震活动增强时,取计算阈值RTL>1.IRTL为RTL值与异常区面积的乘积,可综合反映研究区的地震活动平静或增强的异常程度.IRTL数值大小及随时间的变化趋势与计算所取的空间范围直接相关.震例分析时,空间范围的选取要根据震前研究区的RTL空间分布随时间的演化结果来确定.若震前主要表现为地震活动平静,则选取震中附近地震活动平静最大的区域;反之,则选取震中附近地震活动增强最大的区域.

2.3 区域地震活动

鲁甸地震发生在地震活动较活跃地区,近十年来周边发生了一系列中强地震,如:2003年云南鲁甸MS5.0和MS5.1地震、2004年鲁甸MS5.6地震、2006年云南盐津两次MS5.1地震、2012年云南彝良MS5.6和MS5.7地震(吕坚等,2013)以及2014年云南永善MS5.3地震等(徐涛等,2014).图 1给出了2000年1月1日至2014年8月3日26°N—29°N,102°E—105°E范围内ML1.0以上地震分布(地震目录由中国地震台网中心提供).在地震活动性研究中,对地震目录删除余震以及震级完整性分析是非常重要的.扩散链法(李闽峰,2002)通过链式扩散剔除余震,根据我国震例主震震级、主震与最大余震间距离拟合得到余震的空间范围,时间范围选取K-K法(Keilis-Borok et al.,1980)的时间窗.李闽峰(2002)针对我国多个震例比较了多种删除余震方法,结果显示扩散链法效果较好.因此,本研究主要根据扩散链法剔除余震.

图 1 鲁甸及周边地震震中分布图(2000-01-01—2014-08-03,ML≥1.0) Fig. 1 Distribution of the earthquake epicenters in and around Ludian(2000-01-01—2014-08-03,ML≥1.0)

对地震目录删除余震后,首先进行震级完整性分析.采用了震级-序号图像法(Ogata et al.,1991),即按地震发生时间的先后顺序排序,由地震对应的序号和震级的分布关系讨论地震完整性变化.图 2a给出了地震序列的震级-序号图像,图中不同颜色的色块对应地震数目.同时选取最大曲率法(MAXC)和拟合度分别为90%和95%的拟合度检测(GFT)法(Wiemer and Wyss,2000)定量计算区域地震最小完整性震级MC随时间的变化.其中MAXC法将震级-频度分布曲线一阶导数的最大值对应的震级记为MC.GFT法通过搜索对应给定的实际和理论震级-频度分布下的拟合度来确定MC.本研究中同时分析了拟合度为90%和95%的MC,分别记为GFT-90%和GFT-95%.依据上述不同方法计算得到的MC变化如图 2b,三种方法选择的优先级为GFT-95%>GFT-90%> MAXC,择优后的MC以黑色曲线给出,标注为MC-Best(蒋长胜和吴忠良,2011).

图 2 震级完整性分析 (a)震级-序号图,图中色块的颜色表示地震数目;(b)不同方法给出的完整震级. Fig. 2 Magnitude of completeness analysis (a)Magnitude-serial number: color represents the number of earthquake;(b)Magnitude of completeness from different methods.

图 2可见,不同方法得到的区域地震最小完整性震级MC的时序变化波动较大,2010年之前MC值最大为2.0,之后略有降低.为了满足研究资料选取时段为震前6年的条件,令MC=ML2.0.

3 震前地震活动异常分析

研究资料从2008年1月1日至地震前一天,以每组数据的截止时间标记分析结果,如[2008-01-01—2009-12-31]期间资料的结果记为2009-12-31.

3.1 时间异常

以鲁甸地震震中为研究点,分析震前震中附近地震活动性变化.令特征半径r0=50 km,特征时间t0=12个月,滑移步长10天.为了比较不同震级下限对结果的影响,分别令Mmin=ML2.0、2.1、2.2,其他参数相同,震前RTL值随时间的变化如图 3所示.

图 3 鲁甸地震前震中的RTL值随时间的变化 Fig. 3 The temporal variation of the RTL parameter of the epicenter prior to Ludian earthquake

从2008年1月1日至鲁甸地震发震期间,还发生了2012年9月7日的彝良MS5.7、MS5.6地震(简称彝良地震).可以看出(图 3),取不同震级下限得到的结果数值虽略有差别,但RTL值随时间的变化趋势一致.2010年8月至2011年9月期间的地震活动增强异常可与彝良地震对应.鲁甸地震前,Mmin=ML2.1和ML2.2时,于震前两年检测到地震活动增强,且持续至发震;而Mmin=ML2.0时,于2013年底至2014年初检测到RTL出现负值,但仍接近于0.鉴于这种差别,后文分析中以Mmin= ML2.1和ML2.2的结果为主.

以上地震都发生在RTL值由高值回落到低点再回升的过程中,这反映了震前地震活动经历了不断增强,增强现象减弱,然后继续增强的过程.不同的是两次地震前RTL峰值差别较大:彝良地震前,RTL最大值为5,鲁甸地震前为15.

为了比较空间扫描范围对结果的影响,令r0=60 km和70 km,Mmin=ML2.1,其他参数不变.结果如图 4所示,不同特征半径下的RTL值随时间的变化趋势一致,即检测到的震前地震活动变化一致,但r0=50 km时,稳定性较好,优于其他半径参数的结果.

图 4 不同特征半径得到的RTL值随时间的变化 Fig. 4 The temporal variation of the RTL parameter by different r0
3.2 空间异常

3.1节研究了震中的RTL值随时间的变化,下面将分析震前震中及周边地区的地震活动情况.由于彝良地震的发震时间早于鲁甸地震近两年,因此先介绍彝良地震前地震活动空间异常,再对鲁甸地震进行分析.空间扫描以网格形式,研究点以0.25°N×0.25°E滑移,特征时间t0=12个月,滑移步长10天,特征半径r0=50 km,震级下限Mmin= ML2.1,震级上限小于研究震例的震级.

3.2.1 彝良地震前空间异常

彝良地震发生在西渔河—昭通断裂附近,震中(27.5°N,104.0°E),震源机制与鲁甸地震一致均为走滑型(吕坚等,2013).本文分析了26°N—29°N,102°E—104.5°E的RTL空间分布随时间的变化,计算时使用了[2008-01-01—2012-09-06]期间的地震资料,限于篇幅,文中仅给出了可代表震前地震活动变化过程的结果.冷色表示地震活动平静,暖色表示地震活动增强.结果如图 5所示,彝良地震前检测到地震活动增强,异常主要分布在震中西南侧的则木河断裂与昭通—莲峰断裂的交汇区.2011年5月,在震中西南约100 km检测到小范围异常(图 5a);随着时间的推移,异常范围和异常程度不断增大,于2011年12月达到顶峰(图 5b);之后,又逐渐减弱(图 5c).在整个地震活动增强阶段,异常范围和RTL值总体呈现由小变大再变小的特征.

图 5 彝良地震前RTL空间分布随时间的变化 Fig. 5 The space-time evolution of RTL befor Yiliang earthquakes

根据研究区的RTL分布随时间的演化图像,只能定性地观察地震活动变化,为了定量检验区域地震活动水平,我们计算了图 5b红色不规则多边形区域的RTL面积分,即IRTL随时间的变化.由于异常范围和异常程度在图 5b中达到最大,故根据此时的结果选择异常区,并把这个区域看作彝良地震前地震活动异常的主体区.结果如图 6所示,在研究时段内,IRTL经历了三次由低值升至高值再回落的过程,前两个峰值对应的数值较小,峰值之后分别有4.0≤MS<5.0地震发生;第三个峰值的IRTL=5.51×104 km2,明显高于前两个,峰值后IRTL回落,9个月后彝良地震发生.

图 6 彝良地震前地震活动异常区(图 5b红色曲线区域)的IRTLM-T Fig. 6 The temporal variation of IRTL of the anomalous region(highlighted by red curve in Fig. 5b)prior to Yiliang earthquakes and earthquake sequence
3.2.2 鲁甸地震前空间异常

彝良地震与鲁甸地震相隔近两年,最大震级为MS5.7,明显小于鲁甸地震的震级,计算时应考虑彝良地震的影响.本文使用了[2008-01-01—2014-08-02]期间的地震资料,令Mmin =ML2.1,r0=50 km,分析鲁甸地震前RTL空间分布随时间的变化.图 7给出了可以反映震前地震活动变化过程的结果.震前1年,在震中附近检测到小范围异常,RTL值为3左右;随着时间的推移,在马边构造带与昭通—莲峰构造带围成的区域内检测到地震活动增强,异常程度不断增强,于2014年4月达到顶峰,2014年4月5日永善MS5.3地震也发生在这个地震活动增强区内;至 2014年7月,莲峰断裂东段异常减弱,其他地区的异常较之前变化不明显.总体来讲,地震活动增强主要分布在大凉山次级块体(闻学泽等,2013),2013年8月至2014年4月,异常程度不断增加,之后略有减弱.

图 7 鲁甸地震前RTL空间分布随时间的演化过程 Fig. 7 The space-time evolution of RTL prior to Ludian earthquake

为了定量分析震前区域地震活动水平,计算了图 7e所示的地震活动异常区,即红色不规则多边形区域的IRTL随时间的变化.结果如图 8所示,2011年初至2013年7月,IRTL在背景水平波动,从2013年8月开始不断上升,在上升过程中发生了2014年4月5日永善MS5.3地震;IRTL于2014年5月上旬达到峰值,IRTL=3.45×105 km2,远大于彝良地震前的水平,峰值后3个月鲁甸地震发生.发震时IRTL下降,但数值仍较大,表明该地区的地震活动仍处于增强阶段.

图 8 鲁甸地震前地震活动异常区(图 7e红色曲线区域)的IRTLM-T Fig. 8 The temporal variation of IRTL of the anomalous region(highlighted by red curve in Fig. 7e) prior to Ludian earthquake and earthquake sequence
4 未来地震发生危险性分析

由前文可知,鲁甸地震发生在地震活跃阶段,震后又相继发生了2014年8月17日永善MS5.0、2014年10月1日四川越西MS5.0和2015年1月14日四川乐山MS5.0,这表明川滇东边界中部地区地震活动仍然活跃.为了分析鲁甸地震后,周边地区是否存在发生震级相当或更大的地震的危险,我们跟踪研究了该地区的地震活动变化.

根据他人对RTL方法的研究,当分析一个地震之后,未来更强地震之前的地震活动变化时,要考虑先前发生地震的影响,所以计算资料应包含鲁甸地震并删除余震.特征半径、特征时间、滑移步长、震级下限等参数选取与3.2.2 节相同,计算资料截至2015年10月.本文列出了2014年9月至2015年10月26°N—29.5°N,102°E—104.5°E的RTL分布图.结果如图 9所示,2014年9月至2014年11月,地震活动增强的异常程度逐渐减弱,异常范围缩小至莲峰断裂带中段;但2015年初的地震活动水平再次增强,至2015年10月,异常范围已扩散到昭通—莲峰断裂带、马边断裂带、大凉山断裂带、安宁河—小江断裂带东支围成的活动构造区内.异常区(图 9c红色曲线包围区域)的参数IRTL显示(如图 10):鲁甸地震后该地区的地震活动水平出现短幅上升,从2014年9月开始下降,直至2014年底降至较低水平,即鲁甸地震后区域地震活动水平下降;2015年初IRTL再次回升,于2015年5月底达到较高水平,IRTL=8.25×105 km2,超过了鲁甸地震前的峰值,之后出现小幅度下降,2015年6月至8月上旬,IRTL在一稳定水平波动;2015年8月中旬至10月底,IRTL又持续上升,目前未达到峰值.

图 9 鲁甸地震后川滇交界东边界中段的RTL空间分布随时间的变化 Fig. 9 After Ludian earthquake,the space-time evolution of RTL of middle part of eastern boundary of Sichuan-Yunnan
图 10 川滇东边界中部(图 9c红色曲线区域)2014-01-01—2015-10-26的IRTL随时间的变化 Fig. 10 The temporal variation of IRTL of the middle part of eastern boundary of Sichuan-Yunnan(highlighted by red curve in Fig. 9c)from 2014-01-01 to 2015-10-26

鉴于彝良地震和鲁甸地震发震前,异常范围和异常程度随时间由小变大再变小,区域地震活动水平参数IRTL呈现由低值上升至峰值后回落的过程,且地震都在峰值后数月发震,表明参数IRTL可以反映区域地震活动水平的变化,而且IRTL出现峰值后回落与地震发生可能有一定的相关性.目前大凉山次级块体(闻学泽等,2013)及小江断裂带北段IRTL不断上升,即区域地震活动水平不断增强,可能意味着地震能量的不断积累.再结合他人分析,如:小江断裂带两侧应变速率较高(徐涛等,2014),安宁河南段和则木河北段存在地震空区(闻学泽等,2008)而且分布有多个高应力积累的凹凸体等(易桂喜等,2008),表明安宁河南段、则木河—小江断裂北段、昭通—莲峰断裂带仍存在发生强震的危险,所以川滇交界东边界中段未来的地震形势仍值得深入研究和不断关注.

5 总结与讨论

本文通过RTL方法定量分析了鲁甸地震前地震活动变化,并简要介绍了2012年彝良地震前地震活动异常.根据鲁甸地震和彝良地震震中的RTL值随时间的变化检测到地震活动增强,且观察到地震均发生在RTL值由低值升至高值回落再回升的过程中.这反映了震前地震活动经历了不断增强,增强减弱,再继续增强的过程.此外,着重分析了震前地震活动空间异常演化与地震发生位置的关系:彝良地震前空间异常主要分布在则木河断裂与昭通—莲峰断裂带西段的交汇区,鲁甸地震前空间异常主要分布在大凉山次级块体.总的来讲,异常范围主要分布在震中及周围的潜在发震断裂带及周边,且异常范围和RTL值(异常程度)随时间呈现由小变大再变小的过程.探索了地震活动定量化前兆参数的检测方法,根据新提出的区域地震活动水平参数RTL面积分(IRTL),研究空间范围地震活动的综合异常水平与地震发生时间的关系.在彝良地震和鲁甸地震前,检测到地震活动增强异常区的IRTL随时间由低值升至峰值之后回落,地震发生.彝良地震前,IRTL峰值后9个月发震;鲁甸地震前,IRTL峰值后3个月发震.这表明IRTL峰值对地震发生可能有一定的指示意义.

鲁甸地震发生时,震中和周边地区的地震活动变化率仍高于背景水平,即处于地震活动增强阶段.对该地区跟踪研究显示在大凉山次级块体及小江断裂带北段存在地震活动增强异常,自2015年初该地区的IRTL呈上升趋势,目前IRTL值已明显高于鲁甸地震前的峰值水平,所以川滇交界东边界中段未来的地震形势仍值得深入研究和继续关注.

IRTL衡量的是地震活动异常主体区的综合地震活动水平,其变化趋势受到统计空间范围的影响.选定空间统计区之后,随着时间的演化,若异常区范围较之前变大,则需要修改统计区的范围,应选择新的更大的异常区;如果较之前变小或变化不大,则无需修改.以本文第4节为例,文中选择了图 9c红色曲线包围的区域作为统计区,至2015年10月,地震活动增强区的范围较2015年6月底结果变化不大,故无需修改.此外,根据IRTL变化判断未来地震趋势还受到统计时间的影响,若第4节计算资料截至2015年6月,结果显示峰值出现时间为2015年5月底;而本文计算资料截至2015年10月,IRTL仍处于上升趋势,所以得出了不同的结论.日后我们将开展更多的研究,针对不同震级范围、地区以及地震类型,探索RTL空间分布、区域地震活动水平参数IRTL与地震发生位置、震级及时间的关系.

致谢

在本文研究过程中,江在森研究员、闻学泽研究员、蒋海昆研究员、李闽峰研究员提供了指导,蒋长胜研究员提供了区域震级完整性分析计算程序,文中使用的地震目录来自于中国地震台网中心,部分计算是在中国科学院超级计算中心完成的,审稿专家提出了建设性的意见,在此一并表示感谢.

参考文献
Ak iK. 1965. Maximum likelihood estimate of b in the formula logN=a-bM and its confidence limits. Bull. Earth. Res. Inst. Univ. Tokyo , 43 (2) : 237-239.
Chen C C, Wu Y X. 2006. An improved region-time-length algorithm applied to the 1999 Chi-Chi, Taiwan earthquake. Geophys. J. Int. , 166 (3) : 1144-1148. DOI:10.1111/gji.2006.166.issue-3
Chen S, Wang Q H, Wang Q S, et al. 2014. The 3D density structure and gravity change of Ludian MS6.5 Yunnanepicenter and surrounding regions. Chinese J. Geophys. (in Chinese) , 57 (9) : 3080-390. DOI:10.6038/cjg20140933
Gardonio B, Marsan D, Lengliné O, et al. 2015. Changesin seismicity and stress loading on subduction faults in the Kanto region, Japan, 2011-2014. J. Geophys. Res.:Solid Earth , 120 (4) : 2616-2626. DOI:10.1002/2014JB011798
Gentili S. 2010. Distribution of seismicity before the larger earthquakes in Italy in the time interval 1994-2004. Pure Appl. Geophys. , 167 (8-9) : 933-958. DOI:10.1007/s00024-010-0089-x
Guo Z J, Qin B Y. Causes of Earthquakes and Earthquake Prediction. (in Chinese) Beijing: Seismological Press, 1991 : 177 -178.
Gutenberg B, Richter C F. 1944. Frequency of earthquakes in California. Bull. Seism. Soc. Amer. , 34 (4) : 185-188.
He X H, Ni S D, Liu J. 2015. Rupture directivity of the August 3rd, 2014 Ludian earthquake (Yunnan, China). Science China:Earth Sciences , 58 (5) : 795-804. DOI:10.1007/s11430-015-5053-2
Huang Q H, Sobolev G A, Nagao T. 2001. Characteristics of the seismic quiescence and activation patterns before the M=7.2 Kobe earthquake, January 17, 1995. Tectonophysics , 1995337 (1-2) : 99-116.
Huang Q H, Öncel A O, Sobolev G A. 2002. Precursory seismicity changes associated with the MW=7.4 1999 August 17 Izmit (Turkey) earthquake. Geophys. J. Int. , 151 (1) : 235-242. DOI:10.1046/j.1365-246X.2002.01762.x
Huang Q H. 2008. Seismicity changes prior to the MS8.0Wenchuan earthquake in Sichuan, China. Geophys. Res. Lett , 35 (23) . DOI:10.1029/2008GL036270
Ishibe T, Satake K, Sakai S, et al. 2015. Correlation between Coulomb stress imparted by the 2011 Tohoku-Oki earthquake and seismicity rate change in Kanto, Japan. Geophys. J. Int. , 201 (1) : 112-134. DOI:10.1093/gji/ggv001
Jiang C S, Wu Z L. 2011. Intermediate-term medium-rangeaccelerating moment release (AMR) priori to the 2010 Yushu MS7.1 earthquake. Chinese J. Geophys. (in Chinese) , 54 (6) : 1501-1510. DOI:10.3969/j.issn.0001-5733.2011.06.009
Jiang H K, Hou H F, Zhou H P, et al. 2004a. Region-time-length algorithm and its application to the study of intermediate-short term earthquake precursor in North China. Acta Seismologica Sinica (in Chinese) , 26 (2) : 151-161.
Jiang H K, Hou H F, Zhou H P, et al. 2004b. Application study on Region-Time-Length algorithm for the site judgment of coming moderately strong earthquakes in North China area. Earthquake (in Chinese) , 24 (4) : 17-26.
Keilis-Borok V I, Knopoff L, Rotyain I M. 1980. Bursts of aftershocks, long-term precursors of strong earthquakes. Nature , 283 (5744) : 259-263. DOI:10.1038/283259a0
Li M F. 2002. Study on models & approaches and some associated theories of rapid internet service of earthquake loss estimation[Ph. D. thesis] (in Chinese). Beijing:Institute of Geophysics, China Earthquake Administration.
Liu C L, Zheng Y, Xiong X, et al. 2014. Rupture process of MS6.5 Ludian earthquake constrained by regional broadband seismograms. Chinese J. Geophys. (in Chinese) , 57 (9) : 3028-3037. DOI:10.6038/cjg20140927
Liu P X, Lü X J. 2012. On foreshock identification and short-term predictability of the 2010 Yushu MS7.1 earthquake. Earthquake (in Chinese) , 32 (3) : 47-51.
Liu Y, Lü X J, Tian Q J. 2016. Seismicity changes prior to the MS7.0 Lushan earthquake and the application of regional seismic activity parameter IRTL. Earth Science Frontiers (in Chinese),http://epub.cnki.net/KNS/brief/result.aspx?dbprefix=CJFQ.
Lü J, Zheng X F, Xiao J, et al. 2013. Rupture characteristics and seismogenic structures of the MS5.6 Yiliang earthquakes of Sep.7, 2012. Chinese J. Geophys. (in Chinese) , 56 (8) : 2645-2654. DOI:10.6038/cjg20130814
Mei X P, Liu X F. 2013. Retrospective test on application of region-time-length algorithm to monitoring moderate earthquakes inGansu province and its surrounding areas. Northwestern Seismological Journal (in Chinese) , 35 (S1) : 61-66.
Mogi K. 1979. Two kinds of seismic gaps. Pure Appl.. Geophys. , 117 (6) : 1172-1186.
Ogata Y, Imoto M, Katsura K. 1991. 3-D spatial variation of b-values of magnitude-frequency distribution beneath the Kanto district, Japan. Geophys. J. Int. , 104 (1) : 135-146. DOI:10.1111/gji.1991.104.issue-1
Sobolev G A, Tyupkin Y S. 1997. Low-seismicity precursors of large earthquakes in Kamchatka. Volc. Seismol. , 18 (4) : 433-446.
Sobolev G A, Tyupkin Y S. 1999. Precursory phases, seismicity precursors, and earthquake prediction in Kamchatka. Volc. Seismol. , 20 (6) : 615-627.
Wang W L, Wu J P, Fang L H, et al. 2014. Double difference location of the Ludian MS6.5 earthquake sequences in Yunnan province in 2014. Chinese J. Geophys. (in Chinese) , 57 (9) : 3042-3051. DOI:10.6038/cjg20140929
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
Wen X Z, Du F, Yi G X, et al. 2013. Earthquake potential of the Zhaotong and Lianfeng fault zones of the eastern Sichuan-Yunnan border region. Chinese J. Geophys. (in Chinese) , 56 (10) : 3361-3372. DOI:10.6038/cjg20131012
Wiemer S, Wyss M. 2000. Minimum magnitude of completeness in earthquake catalogs:examples from Alaska, the Western United States, and Japan. Bull. Seism. Soc. Amer. , 90 (4) : 859-869. DOI:10.1785/0119990114
Xu L S, Zhang X, Yan C, et al. 2014. Analysis of the Love waves for the source complexity of the Ludian MS6.5 earthquake. Chinese J. Geophys. (in Chinese) , 57 (9) : 3006-3017. DOI:10.6038/cjg20140925
Xu T, Zhang M H, Tian X B, et al. 2014. Upper crustal velocity of Lijiang-Qingzhen profile and its relationship with the seismogenicenvironment of the MS6.5 Ludian earthquake. Chinese J. Geophys. (in Chinese) , 57 (9) : 3069-3079. DOI:10.6038/cjg20140932
Xu X W, Jiang G Y, Yu G H, et al. 2014. Discussion on seismogenic fault of the Ludian MS6.5 earthquake and its tectonic attribution. Chinese J. Geophys. (in Chinese) , 57 (9) : 3060-3068. DOI:10.6038/cjg20140931
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. (in Chinese) , 51 (6) : 1719-1725.
Zhang G M, Ma H S, Wang H, et al. 2005. Boundaries between active-tectonic blocks and strong earthquakes in the China mainland. Chinese J. Geophys. (in Chinese) , 48 (3) : 602-610. DOI:10.1002/cjg2.693
Zhang G W, Lei J S, Liang S S, et al. 2014. Relocations and focal mechanism solutions of the 3 August 2014 Ludian, Yunnan MS6.5 earthquake sequence. Chinese J. Geophys. (in Chinese) , 57 (9) : 3018-3027. DOI:10.6038/cjg20140926
Zhang Y, Chen Y T, Xu L S, et al. 2015. The 2014 Mw6.1 Ludian, Yunnan, earthquake:A complex conjugated ruptured earthquake. Chinese J. Geophys. (in Chinese) , 58 (1) : 153-162. DOI:10.6038/cjg20150113
陈石, 王青华, 王谦身, 等. 2014. 云南鲁甸MS6.5地震震源区和周边三维密度结构及重力场变化. 地球物理学报 , 57 (9) : 3080–390. DOI:10.6038/cjg20140933
郭增建, 秦保燕. 地震成因和地震预报. 北京: 地震出版社, 1991 : 177 -178.
何骁慧, 倪四道, 刘杰. 2015. 2014年8月3日云南鲁甸M6.5地震破裂方向性研究. 中国科学:地球科学 , 45 (3) : 253–263.
蒋长胜, 吴忠良. 2011. 2010 年玉树MS7.1地震前的中长期加速矩释放(AMR)问题. 地球物理学报 , 54 (6) : 1501–1510. DOI:10.3969/j.issn.0001-5733.2011.06.009
蒋海昆, 侯海峰, 周焕鹏, 等. 2004a. "区域-时间-长度算法"及其在华北中强地震中短期前兆特征研究中的应用. 地震学报 , 26 (2) : 151–161.
蒋海昆, 侯海峰, 周焕鹏, 等. 2004b. "区域-时间-长度算法"在华北地区中强地震地点判定中的应用研究. 地震 , 24 (4) : 17–26.
李闽峰. 2002. 震害预测快速服务平台的模型与方法及一些相关理论的研究[博士论文]. 北京:中国地震局地球物理研究所.
刘成利, 郑勇, 熊熊, 等. 2014. 利用区域宽频带数据反演鲁甸MS6.5级地震震源破裂过程. 地球物理学报 , 57 (9) : 3028–3037. DOI:10.6038/cjg20140927
刘蒲雄, 吕晓健. 2012. 2010年青海玉树7.1级地震的前震识别及短临预报可能性研究. 地震 , 32 (3) : 47–51.
刘月, 吕晓健, 田勤俭. 2016. 芦山MS7.0地震前地震活动性分析及区域地震活动水平参数IRTL的应用. 地学前缘, http://epub.cnki.net/KNS/brief/result.aspx?dbprefix=CJFQ.
吕坚, 郑秀芬, 肖健, 等. 2013. 2012年9月7日云南彝良MS5.7、MS5.6地震震源破裂特征与发震构造研究. 地球物理学报 , 56 (8) : 2645–2654. DOI:10.6038/cjg20130814
梅秀苹, 刘小凤. 2013. "区域-时间-长度算法"用于甘肃及邻近地区中强地震预测的回溯性检验. 西北地震学报 , 35 (S1) : 61–66.
王未来, 吴建平, 房立华, 等. 2014. 2014年云南鲁甸MS6.5地震序列的双差定位. 地球物理学报 , 57 (9) : 3042–3051. DOI:10.6038/cjg20140929
闻学泽, 范军, 易桂喜, 等. 2008. 川西安宁河断裂上的地震空区. 中国科学 , 38 (7) : 797–807.
闻学泽, 杜方, 易桂喜, 等. 2013. 川滇交界东段昭通、莲峰断裂带的地震危险背景. 地球物理学报 , 56 (10) : 3361–3372. DOI:10.6038/cjg20131012
许力生, 张旭, 严川, 等. 2014. 基于勒夫波的鲁甸MS6.5地震震源复杂性分析. 地球物理学报 , 57 (9) : 3006–3017. DOI:10.6038/cjg20140925
徐涛, 张明辉, 田小波, 等. 2014. 丽江-清镇剖面上地壳速度结构及其与鲁甸MS6.5级地震孕震环境的关系. 地球物理学报 , 57 (9) : 3069–3079. DOI:10.6038/cjg20140932
徐锡伟, 江国焰, 于贵华, 等. 2014. 鲁甸6.5级地震发震断层判定及其构造属性讨论. 地球物理学报 , 57 (9) : 3060–3068. DOI:10.6038/cjg20140931
易桂喜, 闻学泽, 苏有锦. 2008. 川滇活动地块东边界强震危险性研究. 地球物理学报 , 51 (6) : 1719–1725.
张国民, 马宏生, 王辉, 等. 2005. 中国大陆活动地块边界带与强震活动. 地球物理学报 , 48 (3) : 602–610.
张广伟, 雷建设, 梁姗姗, 等. 2014. 2014年8月3日云南鲁甸MS6.5级地震序列重定位与震源机制研究. 地球物理学报 , 57 (9) : 3018–3027. DOI:10.6038/cjg20140926
张勇, 陈运泰, 许力生, 等. 2015. 2014年云南鲁甸MW6.1地震:一次共轭破裂地震. 地球物理学报 , 58 (1) : 153–162. DOI:10.6038/cjg20150113