劉 月 邵志剛

(中國地震局地震預測研究所，中國地震局地震預測重點實驗室，北京 100036)

2014年云南景谷MS6.6地震前地震活動變化分析

劉 月 邵志剛

(中國地震局地震預測研究所，中國地震局地震預測重點實驗室，北京 100036)

基于Region-Time-Length(RTL)算法，使用震前約6a和15a的地震目錄定量檢驗了2014年景谷MS6.6地震前的地震活動變化。研究時段近6a的分析結果顯示： 自2013年初在震中附近檢測到地震活動增強，空間異常分布在22.5°～24.5°N，99°～102°E內(nèi)，異常范圍和異常程度隨時間由小變大再變小； 綜合反映區(qū)域地震活動水平的參數(shù)RTL面積分——IRTL從2013年8月開始上升，達到峰值點后回落，峰值點至地震發(fā)生的時間間隔為9個月。研究時段近15a時，研究區(qū)內(nèi)發(fā)生了2007年寧洱MS6.4地震，結果表明： 寧洱地震前檢測到地震活動平靜； 景谷地震前，先檢測到地震活動平靜，后于震前1a檢測到地震活動增強，異常分布在22.5°～25°N，99.5°～101.5°E范圍內(nèi)，IRTL峰值與地震發(fā)生的時間間隔為7個月。這說明即使發(fā)震位置相距較近且震級相當?shù)牡卣穑鹎暗牡卣鸹顒幼兓矔黠@不同。對于景谷MS6.6，不同研究時段檢測到的結果雖有不同，但震前地震活動增強異常區(qū)的分布和參數(shù)IRTL出現(xiàn)峰值的時間較一致，所以探索RTL空間分布和參數(shù)IRTL與發(fā)震位置和時間的關系有重要的研究意義。

Region-Time-Length(RTL)算法RTL面積分 景谷地震 地震危險性分析

0 引言

2014年10月7日云南省普洱市景谷傣族彝族自治縣發(fā)生了MS6.6地震(簡稱景谷地震)，震中23.4°N，100.5°E，震源深度5km(據(jù)中國地震臺網(wǎng)中心)。景谷地震發(fā)生在川滇菱形塊體與密支那-西盟塊體或滇西南塊體交界附近的紅河斷裂西側的無量山區(qū)，位于瀾滄江斷裂以東的思茅-普洱地震帶和以西的耿馬-瀾滄地震帶之間(徐錫偉等，2014; 李永華等，2014)。自21世紀以來，該地區(qū)發(fā)生了多次5級以上地震，包括2007年5月3日的寧洱MS6.4地震。以2014年云南景谷MS6.6為例研究震前震中及周邊地區(qū)的地震活動變化，對認識地震發(fā)生規(guī)律及分析地震危險性有重要意義( Gutenbergetal.，1944; Aki，1965; Mogi，1979; 張國民等，2005; 聞學澤等，2008; Gardonioetal.，2015)。

Region-Time-Length(RTL)算法(Sobolevetal.，1997)以地震目錄為研究資料，根據(jù)RTL(x，y，t)值定量反映震中及周邊地區(qū)t時刻前的地震活動水平，并應用到震例回顧性檢驗及地震發(fā)生危險性研究中。在1995年1月17日阪神MS7.2(Huangetal.，2001)、 1999年8月17 日土耳其MW7.4(Huangetal.，2002)、 2008年5月12日汶川MS8.0(Huang，2008)等強震前檢測到地震活動平靜； 在意大利(Gentili，2010)、 中國華北(蔣海昆等，2004a，b)、 甘肅(梅秀蘋等，2013)的一些中等強度地震前，以及川滇地區(qū)部分MS≥6.0強震前(劉月等，2016a)檢測到地震活動平靜或增強。在地震危險性研究中，綜合衡量區(qū)域地震活動水平的物理參數(shù)RTL面積分(IRTL)，對2008年于田MS7.3、 2008年汶川MS8.0、 2013年蘆山MS7.0及2014年魯?shù)镸S6.5等強震的發(fā)震時間可能存在 “指示”意義(劉月等，2016c，d)，參數(shù)IRTL可為分析未來地震的發(fā)生時間提供一定的依據(jù)。但RTL方法對研究資料的選取時段是否敏感，即對同一研究震例，使用震前不同時段的地震目錄得到的結果是趨于一致，還是存在明顯差異，這決定了對未來研究點(將要發(fā)生的地震)分析結果的準確性。 為了探討RTL方法對分析資料起始時間的依賴性，本文以2014年云南景谷MS6.6地震為例，分別根據(jù)震前約6a和15a的地震資料分析震前RTL空間異常演化與發(fā)震位置及IRTL與發(fā)震時間的關系。

1 相關理論及基礎資料

1.1 RTL算法簡介

Region-Time-Length(RTL)方法(Sobolevetal.，1997)以地震目錄為研究資料，根據(jù)RTL(x，y，t)值衡量研究點(x，y)在t時刻前的地震活動異常程度(方法的具體介紹請參照文獻： 劉月等2016a，b)。RTL(x，y，t)是震中距R(x，y，t)、 時間T(x，y，t)以及破裂尺度函數(shù)L(x，y，t)分別除以各自標準差后的乘積：

(1)

分析過程中涉及到的主要參數(shù)有： 掃描半徑Rmax，特征半徑r0=0.5Rmax，掃描時間窗Tmax，特征時間t0=0.5Tmax，使用地震目錄的最小震級Mmin。實際分析時為了排除正常擾動，把RTL(x，y，t)<-1看作地震活動平靜，RTL(x，y，t)>1看作地震活動增強。

1.2 震級完整性分析

通過擴散鏈法(李閩峰，2002)對地震目錄刪除余震，然后分析2000-01-01—2014-10-06期間，21.5°～25°N，98.5°～102°E的地震目錄的震級完整性。采用了震級-序號圖像法(Ogataetal.，1991)、 最大曲率法(MAXC)及擬合度分別為90%和95%的擬合度檢測(GFT-90%和GFT-95%)法(Wiemeretal.，2000)，分析區(qū)域地震最小完整性震級MC隨時間的變化。對結果選擇的優(yōu)先級為GFT-95%>GFT-90%>MAXC，擇優(yōu)后的MC以黑色曲線給出，標注為MC-Best(蔣長勝等，2011)。結果如圖1 所示，不同方法得到的區(qū)域地震最小完整性震級MC的時序變化結果較為一致，2008年之前MC值最大為ML2.1； 2008年之后為ML1.6。

2 使用2009-01-01—2014-10-06地震目錄分析震前地震活動變化

本節(jié)使用2009-01-01—2014-10-06期間的地震目錄，以每組數(shù)據(jù)的截止時間標記分析結果。分析過程中，令特征時間t0=12個月，則每組分析資料的起止時間相差Tmax=2t0，即Tmax=24個月。例如，第1組資料為2009-01-01—2010-12-31期間的地震目錄，分析結果以數(shù)據(jù)的截止時間來標記，記為2010-12-31(后文相同)。

2.1 震中的地震活動分析

以景谷地震震中為研究點，分析震前震中附近的地震活動變化。令震級下限等于最小完整性震級，即Mmin=ML1.6，特征半徑r0=50km，即掃描半徑Rmax=100km，特征時間t0=12個月，滑移步長為10d。結果如圖2 所示，2012年上半年之前，RTL值在0附近變化，即沒有檢測到地震活動平靜或增強，地震活動在背景水平波動； 自2013年初，RTL值開始上升，于2014年2月初達到峰值，之后回落。景谷地震發(fā)生時，RTL值已降至較低水平。

圖1 震級完整性分析Fig. 1 Magnitude of completeness analysis.a 震級-序號圖，色塊的顏色表示地震數(shù)目； b 不同方法給出的完整震級

圖2 景谷地震前震中的RTL值隨時間的變化Fig. 2 The temporal variation of the RTL of the epicenter prior to Jinggu earthquake.使用2000-01-01—2014-10-06期間的地震目錄，同圖3—8

為了比較震級下限對分析結果的影響，在滿足震級下限大于最小完整性震級的前提下，隨機選取1組數(shù)據(jù)，令Mmin=ML1.6、 1.8、 2.0、 2.3、 2.5、 2.7，其他參數(shù)不變。結果如圖3 所示，根據(jù)不同的Mmin得到的RTL值雖略有差別，但隨時間的變化趨勢一致，均于2013年初至發(fā)震前檢測到地震活動增強。

為了分析空間掃描范圍對結果的影響，分別令r0=25km、 50km、 75km、 100km，Mmin=ML2.0，特征時間t0=12個月，滑移步長為10d。結果如圖4 所示，根據(jù)不同特征半徑得到的RTL值隨時間的變化趨勢一致； 但r0=25km時，波動較大； 當r0=100km時，RTL峰值為13，明顯小于其他結果的峰值，而且自2014年5月檢測到地震活動平靜，在其他特征半徑下沒有檢測到這種特征。這說明r0太小或太大都會影響結果的穩(wěn)定性，甚至會影響對異常趨勢的判斷，所以后文分析時令r0=50km或75km。

圖3 不同震級下限得到的RTL值隨時間的變化Fig. 3 The temporal variation of RTL by different Mmin.

圖4 不同特征半徑得到的RTL值隨時間的變化Fig. 4 The temporal variation of RTL by different r0.

RTL值由無量綱量R、T及L三者共同決定，為了分析它們與RTL的關系，圖5 給出了R、T及L隨時間的變化過程。結果顯示，自2012年三者總體呈上升趨勢，分別于2014年1月底或2月初達到峰值，之后回落。這與RTL隨時間的變化趨勢較為一致。

前文根據(jù)擴散鏈法刪除余震，為了探討分析結果對剔除余震方法的依賴性，分別采用K-K法(Keilis-Boroketal.，1980)和Reasenberg法(Reasenberg，1985)對地震目錄進行處理。令r0=50km，震級下限Mmin=ML2.0，特征時間t0=12個月，滑移步長為10d。3種刪除余震方法得到的景谷地震前的地震活動變化如圖6 所示，2013年3—9月，Reasenberg法檢測到的RTL波動幅度明顯高于其他2種方法。但總的來講，不同刪除余震方法得到的RTL值隨時間的變化相近，所以它們對地震活動變化的判斷結果是一致的。

圖5 無量綱量R、 T、 L隨時間的變化Fig. 5 The temporal variation of the dimensionless qualities R，T and L.

圖6 不同刪除余震方法得到的RTL值隨時間的變化Fig. 6 The temporal variation of RTL by different aftershock deletion methods.

2.2 區(qū)域地震活動分析

上一節(jié)研究了震中的地震活動隨時間的變化，下面將分析震前震中及周邊地區(qū)的地震活動情況，使用2009-01-01—2014-10-06期間的地震資料。空間掃描以網(wǎng)格形式，令研究點以0.25°N×0.25°E滑移，特征時間t0=12個月，滑移步長10d，特征半徑r0=50km，震級下限Mmin=ML2.0。

北緯21.5°～25°，東經(jīng)98.5°～102°的RTL空間分布隨時間的變化如圖7 所示，限于篇幅，僅列出了可代表地震活動變化過程的結果，圖中標識時間為分析資料的截止日期。冷色表示地震活動平靜，暖色表示地震活動增強。在震中以西100km的瀾滄-耿馬地震帶最早于2012年4月初檢測到小范圍的地震活動增強(圖7a)； 隨著時間的推移，異常逐漸向震中擴散，至2013年2月初，在震中及其西北150km范圍內(nèi)檢測到地震活動增強(圖7b，c)； 2013年8月，耿馬-瀾滄地震帶和思茅-普洱地震帶的異常連在一起，RTL值最高的區(qū)域也由之前的震中西北方，轉(zhuǎn)移至其西南方(圖7d)； 之后，異常范圍以震中為中心，集中在22.5°～24.5°N，99°～102°E，且RTL值不斷升高，即地震活動不斷增強(圖7e)； 隨著發(fā)震時間的臨近，異常范圍和異常程度逐漸減弱(圖7f)。

圖7 景谷地震前RTL空間分布隨時間的變化Fig. 7 The space-time evolution of RTL befor Jinggu earthquake.

根據(jù)研究區(qū)的RTL分布隨時間的演化圖像，只能定性地觀察地震活動變化，為了定量檢驗區(qū)域地震活動水平，計算了圖7e紅色矩形區(qū)的RTL面積分——IRTL隨時間的變化。由于異常范圍和異常程度在2014年初幾乎達到最大，故把此時的異常區(qū)看作地震活動增強的主體區(qū)域。結果如圖8 所示： 在2012年之前(含2012年)，IRTL在較低水平波動； 自2013年不斷上升，于2014年1月中旬達到峰值，這意味著區(qū)域地震活動水平不斷增強，并于2014年1月中旬達到頂峰，IRTL=8.85×105km2； 之后，IRTL回落，區(qū)域地震活動逐漸向背景水平恢復。IRTL峰值至地震發(fā)生的時間間隔為9個月。這說明IRTL峰值對發(fā)震時間可能有指示意義。

3 使用2000-01-01—2014-10-06地震目錄分析震前地震活動變化

2007年6月3日發(fā)生了寧洱MS6.4(簡稱寧洱地震)，震中(23.0°N，101.1°E)，該地震的發(fā)震時間與景谷地震相隔約7a，震中距約80km。在一定程度上，寧洱地震或許影響了區(qū)域地震活動變化。為了進一步探索基于RTL算法的地震活動異?？赡艽嬖诘那罢滋卣?，以及寧洱地震對景谷地震的影響，后文將使用更長時間尺度(2000-01-01—2014-10-06)，近16a的地震目錄進行分析。根據(jù)1.3節(jié)可知，區(qū)域最小完整性震級MC=ML2.1。為了與第2節(jié)的分析結果對比，其他參數(shù)選取參照第2節(jié)。

圖8 景谷地震前地震活動異常區(qū)(圖7e紅色矩形包圍區(qū)域)的IRTL隨時間的變化Fig. 8 The temporal variation of IRTL of the anomalous region(highlighted by red rectangle in Figure 7e)prior to Jinggu earthquake.

圖9 不同震級下限得到的RTL值隨時間的變化Fig. 9 The temporal variation of RTL by different Mmin.使用2000-01-01—2014-10-06期間的地震目錄，同圖10—12

3.1 震中的地震活動分析

本節(jié)將分析景谷地震前震中的地震活動變化，特征半徑r0=50km，特征時間t0=12個月，滑移步長為10d。為了比較震級下限對分析結果的影響，分別令Mmin=ML2.1、 2.3、 2.5。結果如圖9 所示，不同Mmin得到的RTL值隨時間的變化趨勢一致，共檢測到3次異常： 1)2003年9月初至2008年4月底檢測到地震活動平靜； 2)2011年初至2012年7月底再次檢測到地震活動平靜； 3)2013年11月至景谷地震前檢測到地震活動增強。下面將分別展開分析。

2003年9月初至2008年4月底的地震活動平靜，可與2007年寧洱地震相對應。在之前的研究中(劉月等，2016a)，檢驗了川滇地區(qū)1976年以來的MS≥6.0強震前震中的RTL值隨時間的變化，對2007年寧洱地震使用了震前約6a的資料，檢測到地震活動平靜，這與本文的分析結果相符。

2008年5月至2010年底，RTL值在0附近波動，即無地震活動平靜，也無地震活動增強。于2011年初至2012年7月底檢測到地震活動平靜，但無強震與之對應。在2.1節(jié)，于相同時段僅檢測到RTL值在0附近波動，未檢測到地震活動平靜。這說明研究資料起始時間的不同影響了RTL值隨時間的變化，進而影響對地震活動趨勢的判斷。

2013年11月至景谷地震前檢測到地震活動增強，持續(xù)時間近1a。與2.1節(jié)相同的是，均于震前檢測到地震活動增強； 不同的是，本節(jié)結果的持續(xù)時間僅1a，<2.1節(jié)中的2a。

圖10 不同特征半徑得到的RTL值隨時間的變化Fig. 10 The temporal variation of RTL by different r0.

為了比較空間掃描范圍對結果的影響，分別令r0=25km、 50km、 75km、 100km，Mmin=ML2.1，特征時間t0=12個月，滑移步長為10d。結果如圖10 所示，r0=50km、 75km、 100km時，RTL值隨時間的變化趨勢一致。r0=25km時，即Rmax=50km，計算結果與其他參數(shù)的明顯不同，不僅RTL值波動較大，且于景谷地震前未檢測到地震活動增強。造成以上差異的原因，可能是掃描半徑過小導致包含的資料信息不足。2.1節(jié)中，r0=25km時，RTL值雖波動幅度較大，但隨時間的變化趨勢與其他r0的計算結果相同。這說明分析資料起始時間不同可能導致分析結果差異較大。

圖11 RTL空間分布隨時間的變化(2003-09-04—2014-07-28)Fig. 11 The time-space evolution of RTL(2003-09-04—2014-07-28).

3.2 區(qū)域地震活動分析

空間掃描時掃描點以0.25°×0.25°滑移，特征時間t0=12個月，滑移步長為10d，特征半徑r0=50km，震級下限Mmin=ML2.1。限于篇幅，文中僅列出了能夠反映寧洱地震前地震活動平靜，以及景谷地震前地震活動平靜和增強的代表性結果，如圖11 所示。五角星分別代表寧洱MS6.4和景谷MS6.6的震中。

寧洱地震前，于2003年9月在紅河斷裂和鎮(zhèn)遠-普洱斷裂之間檢測到地震活動平靜(圖11a)； 至2004年底，異常范圍不斷縮小且逐漸向震中靠攏(圖11b)； 于2005年年中至2006年年中，在震中西北的耿馬-瀾滄斷裂帶之間檢測到大范圍的地震活動平靜(圖11c，d)； 隨著發(fā)震時間的臨近，異常程度減弱，異常范圍縮小，僅在距震中約150km的南汀河西支斷裂和龍陵-瑞麗斷裂之間檢測到地震活動平靜(圖11e)； 臨震前在紅河斷裂和鎮(zhèn)遠-普洱斷裂之間檢測到地震活動增強(圖11e)，持續(xù)數(shù)月后減弱； 之后，地震活動平靜和增強異常均不明顯(圖11f)。通過以上分析可知，寧洱地震前以地震活動平靜為主，異常程度呈現(xiàn)從小變大再變小的過程。

圖12 22°～25°N，89.75°～102°E范圍內(nèi)IRTL隨時間的變化Fig. 12 The temporal variation of IRTL of the region of 22°～25°N，89.75°～102°E.

景谷地震前，先檢測到地震活動平靜，后檢測到地震活動增強。2011年初于震中以西的耿馬-瀾滄斷裂帶檢測到地震活動平靜(圖11g)； 隨著時間的推移，異常的空間分布呈西弱東強的趨勢，于2011年11月集中到紅河斷裂和瀾滄江斷裂之間(圖11h)； 之后異常程度不斷減弱，震中周邊的地震活動逐漸恢復至背景水平，于2013年2月初，在震中以北約150km檢測到小范圍的地震活動增強(圖11i)； 至2013年8月，于震中以南約150km范圍內(nèi)檢測到地震活動增強(圖11j)； 2014年初，震中以南和以北的異常連在一起，空間異常以震中為中心，主要分布在22.5°～25°N，99.5°～101.5°E內(nèi)(圖11k)，這與2.2節(jié)的空間異常分布范圍較接近； 隨著發(fā)震時間的臨近，異常程度逐漸減弱(圖11l)。

以上研究區(qū)的綜合地震活動水平隨時間的變化如圖12 所示： 寧洱地震前以地震活動平靜為前兆異常，參數(shù)IRTL的時程曲線于2006年2月下旬降至低谷點，即IRTL出現(xiàn)低谷后16個月寧洱地震發(fā)震； 景谷地震前先檢測到地震活動平靜，IRTL于2010年至2012年初不斷降低(RTL<-1)，之后回升，自2013年10月IRTL值(RTL>1)逐漸升高，于2014年3月初達到峰值，之后回落，IRTL峰值至地震發(fā)生的時間間隔為7個月。2.2節(jié)的研究結果顯示，IRTL峰值至景谷地震發(fā)生的時間間隔為9個月，這與本節(jié)分析結果僅相差2個月，二者比較一致。這說明對于景谷地震，IRTL峰值與地震發(fā)生時間的關系受研究資料起始時間的影響不大。

4 總結與討論

本文通過RTL算法定量分析2014年云南景谷MS6.6地震前震中及周邊地區(qū)的地震活動變化，為了探索地震活動演化與地震發(fā)生位置和時間的關系，對比了不同研究時段和掃描參數(shù)的結果。分析表明，使用震前約6a的地震目錄時，僅檢測到地震活動增強，異常范圍以震中為中心，分布在22.5°～24.5°N，99°～102°E范圍內(nèi)，綜合反映區(qū)域地震活動水平的物理參數(shù)IRTL達到峰值點后9個月發(fā)震。使用震前約15a的地震目錄時： 景谷地震前先檢測到地震活動平靜，后檢測到地震活動增強，參數(shù)IRTL達到峰值后7個月發(fā)震； 研究區(qū)內(nèi)包括了2007年6月3日寧洱MS6.4地震，檢測到寧洱地震前主要表現(xiàn)為地震活動平靜，參數(shù)IRTL降至低谷點后16個月發(fā)震。景谷地震和寧洱地震的震級相差0.2，且震中距僅80km，但震前地震活動變化明顯不同，這說明基于RTL方法分析地震危險性時，不能僅根據(jù)地震活動平靜或增強判斷發(fā)震震級>7還是<7。

據(jù)自2000年的資料，于2011年初和2011年7月檢測到地震活動平靜，據(jù)2009年以來的地震資料卻未于以上時段檢測到異常，說明研究資料的起始時間不同，可能會導致在同一時段檢測到不同的地震活動變化特征，體現(xiàn)了RTL方法對分析資料時間范圍選取具有一定的依賴性。但是2種分析結果均于景谷地震前1a檢測到地震活動增強，異常集中在耿馬-瀾滄地震帶和思茅-普洱地震帶之間，且IRTL出現(xiàn)峰值的時間僅相差2個月。這表明，基于RTL方法的地震活動增強現(xiàn)象 “預測”景谷地震發(fā)生的危險性時，對發(fā)震位置和時間的 “判斷”是一致的。限于RTL方法對研究資料起始時間存在一定的依賴性，即便檢測到地震活動平靜或增強現(xiàn)象，也不能說明一定會發(fā)生地震。但在分析危險區(qū)與發(fā)震時間時，由于RTL空間異常集中在震中或發(fā)震斷層周邊，這或許可為判斷地震發(fā)生的危險區(qū)域提供依據(jù)，同時IRTL曲線的峰值點或低谷點也可為發(fā)震時間提供參考。

致謝 田勤儉、 江在森、 聞學澤、 蔣海昆、 李閩峰等研究員和呂曉健副研究員對本文提供了指導，蔣長勝研究員提供了區(qū)域震級完整性分析計算程序，文中使用的地震目錄來自中國地震臺網(wǎng)中心，部分計算是在中國科學院超級計算中心完成的，審稿專家提出了建設性的意見，在此一并表示感謝。

郭增建，秦保燕. 1991. 地震成因和地震預報 [M]. 北京: 地震出版社. 177—178.

GUO Zeng-jian，QIN Bao-yan. 1991. Causes of Earthquakes and Earthquake Prediction [M]. Seismological Press，Beijing. 177—178(in Chinese).

蔣長勝，吳忠良. 2011. 2010 年玉樹MS7.1地震前的中長期加速矩釋放(AMR)問題 [J]. 地球物理學報，54(6): 1501—1510.

JIANG Chang-sheng，WU Zhong-liang. 2011. Intermediate-term medium-range accelerating moment release(AMR)prior to the 2010 YushuMS7.1 earthquake [J]. Chinese Journal of Geophysics，54(6): 1501—1510(in Chinese).

蔣海昆，侯海峰，周煥鵬，等. 2004a. “區(qū)域-時間-長度算法”及其在華北中強地震中短期前兆特征研究中的應用 [J]. 地震學報，26(2): 151—161.

JIANG Hai-kun，HOU Hai-feng，ZHOU Huan-peng，etal. 2004a. Region-time-length algorithm and its application to the study of intermediate-short term earthquake precursor in North China [J]. Acta Seismologica Sinica，26(2): 151—161(in Chinese).

蔣海昆，侯海峰，周煥鵬，等. 2004b. “區(qū)域-時間-長度算法”在華北地區(qū)中強地震地點判定中的應用研究 [J]. 地震，24(4): 17—26.

JIANG Hai-kun，HOU Hai-feng，ZHOU Huan-peng，etal. 2004b. Application study on Region-Time-Length algorithm for the site judgment of coming moderately strong earthquakes on North China area [J]. Earthquake，24(4): 17—26(in Chinese).

李閩峰. 2002. 震害預測快速服務平臺的模型與方法及一些相關理論的研究[D]. 北京: 中國地震局地球物理研究所.

LI Min-feng. 2002. Study on models & approaches and some associated theories of rapid internet service of earthquake loss estimation [D]. Institute of Geophysics，China Earthquake Administration，Beijing(in Chinese).

李永華，徐小明，張恩會，等. 2014. 青藏高原東南緣地殼結構及云南魯?shù)椤?景谷地震深部孕震環(huán)境 [J]. 地震地質(zhì)，36(4): 1204—1216. doi: 10.3969/j.issn.0253-4967.2014.04.021.

LI Yong-hua，XU Xiao-ming，ZHANG En-hui，etal. 2014. Three-dimensional crust structure beneath SE Tibetan plateau and its seismotectonic implications for the Ludian and Jinggu earthquakes [J]. Seismology and Geology，36(4): 1204—1216(in Chinese).

劉月，呂曉健. 2016b. 基于區(qū)域-時間-長度算法的區(qū)域地震活動與強震關系的研究 [J]. 地震學報，38(2): 208—219.

LIU Yue，Lü Xiao-jian. 2016b. Relationship between regional seismic activity and strong earthquake based on the region-time-length algorithm [J]. Acta Seismologica Sinica，38(2): 208—219.

劉月，呂曉健，田勤儉. 2016a. 基于 “區(qū)域-時間-長度算法”分析川滇地區(qū)強震前地震活動性變化 [J]. 地震，36(2): 94—104.

LIU Yue，Lü Xiao-jian，TIAN Qin-jian. 2016a. Analysis of seismicity changes prior to strong earthquakes in Sichuan-Yunnan region by the Region-Time-Length algorithm [J]. Earthquake，36(2): 94—104(in Chinese).

劉月，呂曉健，田勤儉. 2016c. 蘆山MS7.0地震前地震活動性分析及區(qū)域地震活動水平參數(shù)IRTL的應用 [J]. 地學前緣(待發(fā)表). http://epub.cnki.net/KNS/brief/result.aspx?dbprefix=CJFQ.

LIU Yue，Lü Xiao-jian，TIAN Qin-jian. 2016c. Seismicity changes prior to theMS7.0 Lushan earthquake and the application of regional seismic activity parameterIRTL[J]. Earth Science Frontiers(To be published)(in Chinese). http://epub.cnki.net/KNS/brief/result.aspx?dbprefix=CJFQ.

劉月，田勤儉，呂曉健，等. 2016d. 魯?shù)镸S6.5地震前地震活動性分析及川滇東邊界中部未來地震危險性研究 [J]. 地球物理學報，59(9): 3269—3279.

LIU Yue，TIAN Qin-jian，Lü Xiao-jian，etal. 2016d. Analysis of the seismicity changes prior to theMS6.5 Ludian earthquake and study on the potential strong-earthquake risk for the middle part of eastern boundary of Sichuan-Yunnan [J]. Chinese Journal of Geophysics，59(9): 3269—3279(in Chinese).

梅秀蘋，劉小鳳. 2013. “區(qū)域-時間-長度算法”用于甘肅及鄰近地區(qū)中強地震預測的回溯性檢驗 [J]. 西北地震學報，35(s1): 61— 66.

MEI Xiu-ping，LIU Xiao-feng. 2013. Retrospective test on application of region-time-length algorithm to monitoring moderate earthquakes in Gansu Province and its surrounding areas [J]. Northwestern Seismological Journal，35(s1): 61— 66(in Chinese).

聞學澤，范軍，易桂喜，等. 2008. 川西安寧河斷裂上的地震空區(qū) [J]. 中國科學(D輯)，38(7): 797—807.

WEN Xue-ze，F(xiàn)AN Jun，YI Gui-xi，etal. 2008. The seismicity gap on the Anninghe fault zone of Sichuan [J]. Science in China(Ser D)，38(7): 797—807(in Chinese).

徐錫偉，程佳，許沖，等. 2014. 青藏高原塊體運動模型與地震活動主體地區(qū)討論： 魯?shù)楹途肮鹊卣鸬膯⑹?[J]. 地震地質(zhì)，36(4): 1117—1134. doi: 10.3969/j.issn.0253-4967.2014.04.015.

XU Xi-wei，CHENG Jia，XU Chong，etal. 2014. Discussion on block kinematic model and future themed areas for earthquake occurrence in the Tibetan plateau: Inspiration from the Ludian and Jinggu earthquakes [J]. Seismology and Geology，36(4): 1117—1134(in Chinese).

張國民，馬宏生，王輝，等. 2005. 中國大陸活動地塊邊界帶與強震活動 [J]. 地球物理學報，48(3): 141—149.

ZHANG Guo-min，MA Hong-sheng，WANG Hui，etal. 2005. Boundaries between active-tectonic blocks and strong earthquakes in the China mainland [J]. Chinese Journal of Geophysics，48(3):141— 149(in Chinese).

Aki K. 1965. Maximum likelihood estimate of b in the formula logN=a-bMand its confidence limits [J]. Bulletin of Earthquake Research Institute，Tokyo University，43: 237—239.

Gardonio B，Marsan D，LenglinéO，etal. 2015. Changes in seismicity and stress loading on subduction faults in the Kanto region，Japan，2011-2014 [J]. Journal of Geophysical Research: Solid Earth，120: 2616—2626. doi: 10.1002/2014JB011798.

Gentili S. 2010. Distribution of seismicity before the larger earthquakes in Italy in the time interval 1994-2004 [J]. Pure and Applied Geophysics，167: 933—958.

Gutenberg B，Richter C F. 1944. Frequency of earthquakes in California [J]. Bulletin of the Seismological Society of America，34: 185—188.

Huang Q H. 2008. Seismicity changes prior to theMS8. 0 Wenchuan earthquake in Sichuan，China [J]. Geophysical Research Letters，35(23). doi: 10.1029/2008GL036270.

Huang Q H，?ncel A O，Sobolev G A. 2002. Precursory seismicity changes associated with theMW=7.4,1999 August 17 Izmit(Turkey)earthquake [J]. Geophysical Journal International，151(1): 235—242.

Huang Q H，Sobolev G A，Nagao T. 2001. Characteristics of the seismic quiescence and activation patterns before theM=7.2 Kobe earthquake，January 17，1995 [J]. Tectonophysics，337: 99—116.

Keilis-Borok V I，Knopoff L，Rotyain I M. 1980. Bursts of aftershocks，long-term precursors of strong earthquakes [J]. Nature，283: 259—263.

Mogi K. 1979. Two kinds of seismic gap [J]. Pure and Applied Geophysics，117(6): 1176—1186.

Ogata Y，Imoto M，Katsura K. 1991. 3-D spatial variation of b values of magnitude frequency distribution beneath the Kanto district，Japan [J]. Geophysical Journal International，14: 135—146.

Reasenberg P. 1985. Second-order moment of central California seismicity，1969-1982 [J]. Journal of Geophysical Research，90(B7): 5479—5495.

Sobolev G A，Tyupkin Y S. 1997. Low-seismicity precursors of large earthquakes in Kamchatka [J]. Volcanology and Seismology，18: 433— 446.

Wiemer S，Wyss M. 2000. Minimum magnitude of complete reporting in earthquake catalogues: Examples from Alaska，the Western United States，and Japan [J]. Bulletin of the Seismological Society of America，90: 859—869.

ANALYSIS OF SEISMICITY CHANGES PRIOR TO THE 2014 YUNNAN JINGGUMS6.6 EARTHQUAKE

LIU Yue SHAO Zhi-gang

(InstituteofEarthquakeScience，ChinaEarthquakeAdministration，Beijing100036，China)

According to the Region-Time-Length(RTL)algorithm，the analysis of seismicity changes prior to the 2014 Yunnan JingguMS6.6 earthquake was conducted by using the earthquake catalogues about 6 and 15 years before this earthquake，respectively. When the studied period was nearly 6 years，an enhancement of seismic activity was detected around the epicenter since the beginning of 2013. The anomalies mainly distributed in the region of 22.5°～24.5°N and 99°～102°E. The range and degree of anomalies changed from small to large，and then to small chronologically. As the surface integral in respect toRTL，the physical parameterIRTL，which could reflect the regional seismicity level，began to increase since August 2013，and then reduced after reaching the peak point. The time length from the peak point ofIRTLcurve to the earthquake occurrence was 9 months. When the analyzed catalogue was nearly 15 years，the 2007 NingerMS6.4 occurred in the studied region. Seismicity quiescence was detected prior to the NingerMS6.4. Before the JingguMS6.6，seismicity quiescence was detected firstly，and then enhanced activity was observed 1 year prior to the earthquake occurrence. The anomalies mainly distributed in the region of 22.5°～24.5°N and 99°～102°E. The time length from the peak point ofIRTLcurve to the earthquake occurrence was 7 months. The above study showed that even the earthquakes location was near and the magnitude was close to each other，a big difference in seismic activity before the earthquakes may exist. Before the JingguMS6.6，there was some difference in seismicity changes according to different beginning time of catalogues，but the distribution of anomalies and the time length from the peak point ofIRTLto the earthquake occurrence were uniform. So there was an important significance for exploring the relationship between the distribution of anomalies and the earthquake location，and the relationship between the time of the peak point ofIRTLand the earthquake occurrence time.

Region-Time-Length(RTL)algorithm，surface integral ofRTL，Jinggu earthquake，seismic hazard evaluation

10.3969/j.issn.0253- 4967.2016.04.021

2016-02-19收稿，2016-10-27改回。

中國地震局地震預測研究所基本科研業(yè)務專項(2015 IES010204)資助。

P315.2

A

0253-4967(2016)04-1070-12

劉月，女，1987年生，2014年于中國科學院力學研究所獲得固體力學專業(yè)博士學位，助理研究員，主要從事地球動力學和地震活動性研究，電話： 010-88015344，E-mail： liuyue@cea-ies.ac.cn。