黎明曉 王欽軍 盧顯 解孟雨 姚琪 薛艷 閆偉
1)中國地震臺網(wǎng)中心,北京 100045 2)中國科學(xué)院空天信息創(chuàng)新研究院,數(shù)字地球重點實驗室,北京 100094 3)中國地震局地震預(yù)測研究所,北京 100036
自20世紀80年代開始,人工神經(jīng)網(wǎng)絡(luò)進入發(fā)展期,衍生出不同的模型,在圖像識別、評估預(yù)測、醫(yī)藥衛(wèi)生、地質(zhì)等領(lǐng)域得到廣泛研究與應(yīng)用(Rumelhart et al,1986;曾坤等,2021;吳婷婷等,2021;王麗紅,2021;王良玉等,2021)。人工神經(jīng)網(wǎng)絡(luò)不需要假設(shè)條件或提前建立預(yù)測函數(shù)關(guān)系或模型,通過自學(xué)習(xí)過程總結(jié)出非線性因果關(guān)系,因此被引入到地震預(yù)測中,并得到一定的應(yīng)用(馮德益等,1992;李東升等,1995;王煒等,2000;項月文等,2016)。其中,BP(Back Propagation)神經(jīng)網(wǎng)絡(luò)是使用較廣泛的神經(jīng)網(wǎng)絡(luò)模型之一,又稱反向傳播前饋性神經(jīng)網(wǎng)絡(luò)。其算法理論與應(yīng)用技術(shù)成熟,具有高適應(yīng)性的非線性系統(tǒng),可表達各種復(fù)雜映射,包括地震多因子的非線性變化規(guī)律。BP神經(jīng)網(wǎng)絡(luò)采用梯度搜索技術(shù),誤差向后傳播的同時修正權(quán)系數(shù),使網(wǎng)絡(luò)的實際輸出與期望輸出均方差最小(馮德益等,1992;李東升等,1995)。
多年來,地震學(xué)家開展了強震發(fā)生前多種參數(shù)的特征研究,探索強震前的中短期(1年內(nèi))異常指標。常用的地震活動參數(shù)有:古登堡公式系數(shù)b值(Aki,1965;Gutenberg et al,1965;易桂喜等,2013;史海霞等,2018)、震級-頻度關(guān)系偏離度η(馮德益等,1994)、地震活動性指數(shù)A(b)(吳佳翼等,1983;張瑞芳等,2004)、地震強度因子Mf(王煒等,1994)、地震空間集中度C(王煒等,1989)、地震危險度D(王煒等,1987)、算法復(fù)雜性Cn(呂悅軍等,1993)、廣義時間分維Dt值(孫家林等,1995)等。這些參數(shù)的震前異常在時間、空間、強度上所呈現(xiàn)的特征不同,有的參數(shù)震前異常顯著,有的無異常,有的在空間上有異常但不一定發(fā)生強震等。這些參數(shù)的異常變化特征與未來強震之間存在較強的非線性關(guān)系,難以通過某種函數(shù)關(guān)系表達。因此,引入BP神經(jīng)網(wǎng)絡(luò)方法,將地震預(yù)報中常用的具有中期或中短期預(yù)測時效的地震學(xué)參數(shù)b、η、A(b)、Mf、C、D、Cn、Dt及其差值組合,作為BP神經(jīng)網(wǎng)絡(luò)的輸入端,將神經(jīng)網(wǎng)絡(luò)的輸出值稱為地震活動平靜特征參數(shù)Wq(王煒等,2000)。若Wq值接近或等于1,表明已處于地震孕育的中短期階段;若Wq值接近或等于0,表明未處于發(fā)震時段。該方法被應(yīng)用于華北、青海和新疆等地區(qū),結(jié)果表明Wq值在中強震前有較好的中短期預(yù)測效果(王煒等,2000;張軍等,2002;王培玲等,2008;孫麗等,2016)。
Wq值方法在中國地震臺網(wǎng)中心年度、年中地震趨勢預(yù)測中得以廣泛應(yīng)用,積累了大量震例。如2007年5月5日新疆于田6.1級、2008年8月30日四川會理6.1級、10月5日新疆烏恰6.8級、2016年12月8日新疆呼圖壁6.2級和2020年1月19日新疆烏恰6.4級等地震均發(fā)生在Wq值方法異常區(qū)域內(nèi)或附近(中國地震臺網(wǎng)中心,2006、2007、2016;孟令媛等,2020)。2019年年底的計算結(jié)果顯示,青海、新疆和西藏交界地區(qū)存在較大面積的Wq值異常;2020年年中該區(qū)域異常進一步增強。在異常區(qū)域先后發(fā)生了2020年6月26日新疆于田6.4級和7月23日西藏尼瑪6.6級地震,表明Wq值對該區(qū)的中期、中短期強震預(yù)測具有一定意義。但以往震例顯示W(wǎng)q異常區(qū)域往往較大,依據(jù)Wq值方法判定未來發(fā)生地震的強度和危險區(qū)域一直是難點。本研究借助豐富的震例建立震級與Wq值異常區(qū)面積之間的統(tǒng)計關(guān)系,以此估計地震強度,通過Wq值差值變化縮小地震危險區(qū)域。
為了更深入地研究2020年新疆于田6.4級和西藏尼瑪6.6級地震前Wq值的中短期演化特征,將計算Wq值的滑動步長由年度會商時的1年、年中會商時的半年縮短至本研究的1個月,計算時間窗長為1年。嘗試建立震級與Wq值異常區(qū)面積的統(tǒng)計關(guān)系,尋求強度判定的定量依據(jù);通過不同期Wq值異常變化空間分布特征,縮小圈定未來發(fā)生強震的危險地區(qū),為在地球物理觀測較少而強震頻發(fā)的西藏及周邊地區(qū)開展強震中短期預(yù)測提供借鑒。
青海、新疆和西藏地區(qū)構(gòu)造運動強烈,主體位于青藏高原及周邊,受印度板塊的推擠作用,自北向南形成了東昆侖—柴達木地塊、巴顏喀拉地塊、羌塘地塊和拉薩地塊(鄧起東等,2014),新疆處于中國大陸西側(cè),長期處于擠壓環(huán)境下,形成新生代再生造山帶和壓陷盆地(圖1)。該研究區(qū)中強及以上地震頻發(fā),歷史上曾發(fā)生多次7級、8級地震,2010—2020年發(fā)生5級以上地震122次,占大陸發(fā)生的5級以上地震的58%,年平均發(fā)生5級以上地震11次。其中,發(fā)生6.0~6.9級地震16次,7.0級以上地震2次,分別為2010年4月14日青海玉樹7.1級和2014年2月12日新疆于田7.3級地震。該研究區(qū)2020年1個月內(nèi)先后發(fā)生6月26日新疆于田6.4級和7月23日西藏尼瑪6.6級地震,分別位于巴顏喀拉地塊和羌塘地塊內(nèi),2次地震相距約500km。上述2020年2次6級地震為本次Wq值方法預(yù)測應(yīng)用實例(圖1)。
注:震源機制解采用中國地震臺網(wǎng)中心計算結(jié)果;①阿爾金斷裂帶;②向陽湖—布若錯斷裂;③拉雄錯—噶爾孔茶卡;④瑪依崗日山;⑤依布茶卡—日干配錯斷裂。
2020年6月26日新疆于田6.4級地震發(fā)生在阿爾金斷裂帶西段,為張性地震(圖1),位于2008年3月21日和2014年2月12日2次于田7.3級地震之間。阿爾金斷裂帶是青藏高原北緣一條左旋走滑為主的主控邊界斷裂帶(國家地震局《阿爾金活動斷裂帶》課題組,1992;李玉江等,2015),斷裂帶中、西段構(gòu)造運動較強,左旋走滑運動速率可達(17.5±2)mm/a(徐錫偉等,2003)。
2020年7月23日西藏尼瑪6.6級地震發(fā)生在西藏羌塘地塊中部,為依布茶卡—日干配錯斷裂帶次級構(gòu)造上的一次張性地震(圖1)。羌塘地塊是一橫貫青藏高原中部的近EW向的活動地塊,整體被向東擠出,通過其內(nèi)部的NE和NW向2組共軛剪切破裂進行EW向伸展(馬宗晉等,1998;陳祖安等,2008)。尼瑪6.6級地震周邊主要發(fā)育了2組構(gòu)造,一組為NEE-NE向依布茶卡—日干配錯、瑪依崗日山斷裂,另一組為NW向向陽湖—布若錯斷裂、拉雄錯—噶爾孔茶卡斷裂等,其中NW向斷裂的地震活動性不強。尼瑪6.6級地震發(fā)生在依布茶卡—日干配錯斷裂帶的分支斷裂上,該斷裂為左旋走滑斷裂,并兼具一定的張性活動特征,基本呈NE走向,與依布茶卡—日干配錯斷裂呈較小角度相交,這組斷裂具有很強的地震活動,曾發(fā)生過多次中強地震。2008年1月依布茶卡—日干配錯斷裂西段曾發(fā)生西藏改則6.9級、6.0級地震。
該研究區(qū)地球物理觀測臺站分布不均,天山地震帶中北部地球物理觀測密度相對較高,其他區(qū)域觀測密度較低,尤其是西藏中西部至新藏交界及附近地區(qū)(30°~38°N,78°~90°E)地球物理觀測臺站密度較低,96×104km2內(nèi)僅有13項定點觀測測項(圖1(b))。采用地震活動平靜特征參數(shù)Wq方法在研究區(qū)進行系統(tǒng)掃描具有一定的優(yōu)勢,研究區(qū)中強及以上地震頻繁,有利于開展地震活動參數(shù)的中短期特征及震例研究。通過總結(jié)青海、新疆和西藏地區(qū)2010—2020年5級以上地震前Wq值異常特征,研究震級與Wq值異常區(qū)面積的正相關(guān)統(tǒng)計關(guān)系、危險區(qū)域和發(fā)震時間特征,以期為震情跟蹤與預(yù)測提供參考。
圖2 2000—2020年各時段地震累積頻度-震級分布
根據(jù)以上地震完備性分析,采用Wq值方法對研究區(qū)不同時間段內(nèi)的地震進行空間掃描時,地震目錄數(shù)據(jù)采用中國地震臺網(wǎng)中心月報目錄,起始震級為ML3.0。
地震預(yù)報中常用的b、η、A(b)、Mf、C、D、Cn、Dt等地震學(xué)參數(shù)異常變化與未來地震之間表現(xiàn)為很強的非線性關(guān)系,并且A(b)、Mf、C、D值中短期特征表現(xiàn)為高值異常的恢復(fù),b、η、Cn、Dt值表現(xiàn)為低值異常的恢復(fù)。為綜合表征當前各參數(shù)的大小以及與其前、中期異常階段的差異,將地震預(yù)報中常用的b、η、A(b)、Mf、C、D、Cn、Dt值8個地震學(xué)參數(shù)與其前2個相鄰時段的差值(合計16個參數(shù))作為神經(jīng)網(wǎng)絡(luò)的輸入值。其差值ΔXi按如下方式確定
(1)
神經(jīng)網(wǎng)絡(luò)的輸出值反映該區(qū)域地震活動在中短期階段平靜的綜合特征,稱之為地震活動平靜特征參數(shù)Wq(王煒等,2000)。當Wq值接近1時,表明已處于地震孕育的中短期階段,為異常;若Wq值接近0,表明未處于發(fā)震時段,無異常。據(jù)王煒等(2000)研究結(jié)果,中強地震前1年左右時間,未來震中周圍一般開始出現(xiàn)Wq值的中短期異常,因此式(1)中Xi為當前值,Xi-1為當前值之前半年的參數(shù)值,Xi-2為當前值之前1年的參數(shù)值。
采用Wq值方法對研究區(qū)不同時間段內(nèi)的地震進行空間掃描時,除了地震目錄數(shù)據(jù)采用中國地震臺網(wǎng)中心月報目錄、起始震級為ML3.0外,其他計算參數(shù)選擇參考西部地區(qū)Wq研究成果(王培玲等,2008;孫麗等,2016),時空掃描計算的時間窗長為12個月,滑動步長為6個月,空間掃描單元窗長為300km,滑動步長為50km。在對2020年新疆于田6.4級地震和西藏尼瑪6.6級地震進行Wq值時空掃描時,為細致地分析震前Wq值異常的時間演化特征,將時間滑動步長定為1個月。
在2010—2020年年度及年中Wq值跟蹤研究積累震例的基礎(chǔ)上,按上述計算方法于每年5月底和11月底對青海、新疆和西藏地區(qū)進行Wq值空間掃描,并開展5級以上地震預(yù)測檢驗,獲得近十余年Wq值震例數(shù)據(jù)。在震例研究中做如下約定:假設(shè)采用計算窗長為2018年6月1日至2019年5月31日的數(shù)據(jù)計算的Wq值,其異常開始時間記為2019年5月31日;震例的Wq值異常面積選取1年內(nèi)震中附近Wq值≥0.5的最大異常面積;同一Wq值異常對應(yīng)多次地震時,對應(yīng)地震震例僅記錄為最大震級的地震,其他對應(yīng)的地震未列入統(tǒng)計;對于在Wq值主異常區(qū)內(nèi)的5級以上地震,或者在Wq值主要異常區(qū)外100km范圍內(nèi)的5級以上地震以及在Wq值主要異常區(qū)外200km范圍內(nèi)的6級及以上地震,均視為對應(yīng)震例。經(jīng)統(tǒng)計,共有36次5級以上地震前出現(xiàn)顯著的Wq值異常(圖3),占同時段122次5級以上地震的29.5%。在此基礎(chǔ)上,開展Wq震例的時空強三要素特征研究。
圖3 2010—2020年青海、新疆和西藏地區(qū)Wq異常對應(yīng)的5級以上地震震中分布
對2010—2020年青海、新疆和西藏地區(qū)5級以上地震前有顯著Wq值異常的36次震例進行震級-異常面積的統(tǒng)計分析,結(jié)果見圖4。地震震級(MS)與Wq值異常面積(SWq)統(tǒng)計正相關(guān),相關(guān)系數(shù)達0.85,關(guān)系表達式為
圖4 地震震級與Wq值異常面積關(guān)系圖
SWq=16.774×MS-76.836
(2)
將震例的異常面積代入式(2),推算未來發(fā)生地震的強度,如圖5所示。預(yù)測震級走勢與實際發(fā)生地震基本一致,震級差在±0.2級的震例約占47%,震級差在±0.3級的震例約占61%,震級差在±0.5級的震例約占86%。具體地震震例參數(shù)及預(yù)測震級等信息見表1,表中地震序號與圖5一致。
圖5 Wq值異常面積推算的震級與實際發(fā)生地震的震級
表1 Wq值預(yù)測的震例信息統(tǒng)計
確定Wq值異常最大面積是判別強度大小的關(guān)鍵。目前震情研判多采用強度區(qū)間預(yù)測,如6~7級、7級,按照這種方式估算震級的所在區(qū)間與實際發(fā)生地震的震級較為一致。
2010—2020年青海、新疆和西藏地區(qū)共發(fā)生5級以上地震107次(去除余震),其中5.0~5.9級地震87次,6.0級以上地震20次。各震級分檔在震前出現(xiàn)Wq異常的數(shù)量及比例情況見圖6,由圖可見,隨地震震級的增加,地震數(shù)量越來越少,其對應(yīng)Wq值異常預(yù)測地震的比例越來越高。5.0~5.4級和5.5~5.9級地震震前出現(xiàn)Wq年度異常預(yù)測地震的比例分別為26.1%、27.8%,6.0~6.4級和6.5~6.9級地震前出現(xiàn)Wq值異常預(yù)測地震的比例相對較高,分別為60.0%和62.5%。由此可見,在青海、新疆和西藏地區(qū),Wq值方法對6級以上地震的預(yù)測效果相對較好。
圖6 2010—2020年青海、新疆和西藏地區(qū)Wq異常對應(yīng)5級以上地震情況
統(tǒng)計5級以上地震與Wq值出現(xiàn)異常的時間間隔,各發(fā)震時間統(tǒng)計區(qū)間內(nèi)的地震占震例總數(shù)的百分比見圖7(a),盡管發(fā)震時間受統(tǒng)計對應(yīng)1年內(nèi)地震的影響,但發(fā)生在9個月內(nèi)的震例占總數(shù)的百分比約為81%,顯示W(wǎng)q值異常具有中短期特征。
圖7 2010—2020年5級以上地震與震前Wq值異常的統(tǒng)計關(guān)系
統(tǒng)計5級以上地震震中與Wq值異常區(qū)的關(guān)系,各空間統(tǒng)計區(qū)間內(nèi)的地震占震例總數(shù)的百分比見圖7(b)。發(fā)生在Wq值異常區(qū)內(nèi)的震例接近70%,22.3%的地震發(fā)生在異常區(qū)外100km范圍內(nèi),發(fā)生在異常區(qū)外100~200km的地震為8.3%,表明地震多數(shù)發(fā)生在Wq值異常區(qū)內(nèi),但通常異常面積較大,對危險區(qū)域判定帶來較大的不缺定性。為進一步縮小危險區(qū)域,將Wq值變化較大的2期計算結(jié)果做差值,分析應(yīng)用實例(于田6.4級和尼瑪6.6級地震)與Wq值異常區(qū)域的空間關(guān)系。
利用上述青海、新疆和西藏地區(qū)Wq值異常與未來地震之間關(guān)系,對2020年新疆于田6.4級和西藏尼瑪6.6級地震進行回溯式時空強三要素研判。
采用上述介紹的地震活動平靜特征參數(shù)Wq值方法,獲得2019—2020年研究區(qū)Wq值空間掃描、月滑動計算結(jié)果。在2020年6月26日新疆于田6.4級地震和2020年7月23日西藏尼瑪6.6級地震前、后部分時段的Wq值空間分布見圖8。
圖8 西藏中西部至新疆、西藏交界地區(qū)Wq值掃描結(jié)果
自2019年11月30日(于田6.4級地震前約7個月)開始至2020年1月31日,青海、新疆和西藏交界地區(qū)的Wq值顯著異常區(qū)出現(xiàn)沿新藏交界地區(qū)逐漸向西擴展(圖8(a)、8(b)),異常面積最大約33.4×104km2(圖8(b))。2020年2月(于田6.4級地震前約4個月)擴展區(qū)域的Wq值異常面積顯著減小并恢復(fù),2020年6月26日新疆于田6.4級地震就發(fā)生在Wq值異常面積擴展—減小恢復(fù)的地區(qū)(圖8(b)~8(d)),且在于田6.4級地震發(fā)生后震中附近Wq值異常消失,于田6.4級地震東側(cè)的Wq值異常仍持續(xù)存在(圖8(e)),該區(qū)域后續(xù)發(fā)生了尼瑪6.6級地震。
2019年10月31日(西藏尼瑪6.6級地震前近9個月)青海、新疆和西藏交界地區(qū)Wq值異常(2018年11月1日—2019年10月31日,圖8(a))除向西的擴展區(qū)發(fā)生了于田6.4級地震外,還向南擴展至西藏南部,Wq值最大異常面積約41.9×104km2(圖8(c)),震前約2個月異常面積略有減小(圖8(d)),西藏尼瑪6.6級地震發(fā)生在Wq值異常面積增大的區(qū)域。在尼瑪6.6級地震前的演化過程中,Wq值異常區(qū)最南端還發(fā)生了2020年3月20日西藏定日5.9級地震。西藏定日5.9級地震震前約3個月(2019年12月)開始,震中附近的西藏南部定日地區(qū)出現(xiàn)Wq值異常,異常持續(xù)約3個月后發(fā)震,定日5.9級地震發(fā)生在Wq值擴展區(qū)。至2020年10月底,尼瑪6.6級地震震中區(qū)域的Wq值異常完全消失,定日5.9級地震震中周邊區(qū)域的Wq值異常始終未消失(圖8(f))。
對于田6.4級和尼瑪6.6級地震前后Wq值月滑動結(jié)果分析表明,在2次6級地震前Wq值表現(xiàn)為中短期特征:地震發(fā)生在Wq值異常的擴展或擴展—減小時段,地點位于Wq值異常面積擴展區(qū)或擴展—減小恢復(fù)區(qū)。該結(jié)論與在青海、新疆地區(qū)中強震前Wq值的中短期研究成果一致(孫麗等,2016;王培玲等,2008)。
2020年7月23日,西藏尼瑪6.6級地震前的Wq值異??臻g范圍較大,且在演化過程中還疊加了2020年6月26日新疆于田6.4級地震和3月20日西藏定日5.9級地震的異常信息,對判定發(fā)生強震的危險地點和強度造成了一定干擾。以2019年5月31日(計算窗長為2018年6月1日—2019年5月31日)計算的Wq值為基準值,分別將震前Wq異常面積擴展較大時段2019年11月30日、12月31日、2020年3月31日、4月30日計算的Wq值與基準值做差值,得到Wq值變化值空間分布圖(圖9),從圖中可以看出,2020年6月26日于田6.4級地震發(fā)生在Wq值異常擴展區(qū)內(nèi)(圖9(a)),地震后異常持續(xù)向南擴展,7月23日西藏尼瑪6.6級地震位于Wq值異常擴展區(qū)中部,擴展區(qū)南側(cè)發(fā)生了2020年3月20日西藏定日5.9級地震(圖9(b)~9(d))。該結(jié)論與王培玲等(2008)、張軍等(2002)得到的地震一般發(fā)生在Wq值顯著擴展變化區(qū)域的認識一致。因此,Wq值異常顯著擴展區(qū)更靠近未來發(fā)生中強震的震中。
圖9 于田6.4級地震和尼瑪6.6級地震前Wq差值變化分布
將于田6.4級地震前Wq最大異常面積33.4×104km2、尼瑪6.6級地震前Wq最大異常面積41.9×104km2分別代入公式(2),計算震級分別為6.6級和7.1級。根據(jù)Wq值異常面積估算的于田地震震級比實際發(fā)生的震級高0.2,比較接近實際值;尼瑪?shù)卣鸸浪阏鸺壉葘嶋H發(fā)生的震級高0.5,估算震級略偏高,其主要原因是尼瑪6.6級地震前的Wq值異常中疊加了2020年3月20日西藏定日5.9級地震前的異常信息,可能導(dǎo)致尼瑪6.6級地震前Wq異常面積偏大,進而估算震級略高。由此可見,可以通過Wq異常面積推算可能發(fā)生地震的震級,并適用于震情研判中通常采用的強度區(qū)間(±0.5級)預(yù)測。
(1)在青海、新疆和西藏地區(qū),中強震前存在Wq值中短期異常,Wq值可作為該區(qū)強震前的中短期地震學(xué)綜合參數(shù)。
(2)提供了青海、新疆和西藏地區(qū)Wq值異常區(qū)面積與未來地震震級之間的定量統(tǒng)計關(guān)系。
本文建立的青海、新疆和西藏地區(qū)Wq值異常面積與震級的擬合關(guān)系式結(jié)果表明,地震震級與Wq值異常面積呈顯著正相關(guān),二者的判定系數(shù)為0.85,異常面積越大,可能未來發(fā)生地震的強度越高。根據(jù)式(2)統(tǒng)計關(guān)系,可以利用Wq異常區(qū)面積來定量預(yù)測將來發(fā)生地震的強度,震級差在±0.2級的震例約占47%,震級差在±0.3級的震例約占61%,震級差在±0.5級的震例約占86%。目前,一般在震情跟蹤中對強度采用“*級左右”(上下浮動0.5級)的方式,該方法可以運用到現(xiàn)有震情跟蹤研判中,正確率相對較高。
(3)Wq值差值方法可進一步縮小未來發(fā)生地震的危險區(qū)域。
未來地震震中一般位于Wq值異常區(qū),但實際上該異常區(qū)范圍較大,進一步縮小危險地區(qū)是關(guān)鍵。于田6.4級和尼瑪6.6級地震前,Wq值時空差值變化研究結(jié)果顯示,將Wq值變化較大的兩期計算結(jié)果做差值處理,可以更好地確定未來發(fā)生地震的危險地區(qū)。
(4)年度Wq值方法對青海、新疆和西藏地區(qū)6級以上地震的預(yù)測效果高于5級地震。
2010—2020年青海、新疆和西藏地區(qū)5級以上地震年度Wq值震例總結(jié)顯示,6.0~6.9級地震前出現(xiàn)Wq年度異常占61%,5.0~5.9級地震震前出現(xiàn)Wq年度異常占25.0%,年度Wq值異常對6級以上地震的對應(yīng)率較高,預(yù)測能力較好。這與青藏高原絕大多數(shù)6級地震前存在ML4.0以上地震大范圍平靜現(xiàn)象(陳玉華等,2014)較一致,從側(cè)面反映出6級地震前大多出現(xiàn)區(qū)域地震活動增強轉(zhuǎn)平靜的中短期特征。除此之外,可能與研究區(qū)的地震監(jiān)測能力相對較弱、掃描參數(shù)震級的選定有一定關(guān)系。因此,目前的研究成果反映地震前平靜參數(shù)的Wq異常對6級地震預(yù)測效果更好。
致謝:審稿專家對論文修改提出了中肯的意見與建議,本文撰寫得到中國地震臺網(wǎng)中心劉杰研究員的指導(dǎo),中國地震臺網(wǎng)中心國家地震科學(xué)數(shù)據(jù)中心(http://data.earthquake.cn)提供數(shù)據(jù)支撐,在此一并表示衷心感謝。