李瓊 付虹 毛慧玲 朱榮歡 何德強(qiáng)

摘要：選取2009—2015年云南地區(qū)觀測(cè)連續(xù)、穩(wěn)定的水溫觀測(cè)資料共34項(xiàng)，采用異常過(guò)程和從屬函數(shù)分析方法聯(lián)合提取異常，分析水溫異常與MS≥50地震的關(guān)系。結(jié)果表明：該方法提取的單項(xiàng)異常，具有短期指示意義，與地震有較好的對(duì)應(yīng)關(guān)系，但地震的漏報(bào)率較高；采用異常頻次和權(quán)重集成方法提取的全省群體異常與全省的MS≥50地震對(duì)應(yīng)較好，減少了漏報(bào)，并可以通過(guò)置信度為975%的R值評(píng)分檢驗(yàn)。按構(gòu)造和地震活動(dòng)進(jìn)行分區(qū)后，提取小區(qū)域的水溫群體異常，并只考慮對(duì)應(yīng)本區(qū)的地震，則各區(qū)域內(nèi)出現(xiàn)異常的虛報(bào)和漏報(bào)率均顯著增加。分析其原因，認(rèn)為水溫異常和地震都是區(qū)域構(gòu)造活動(dòng)增強(qiáng)的產(chǎn)物，群體異常能更好地刻畫(huà)區(qū)域構(gòu)造活動(dòng)增強(qiáng)，在構(gòu)造活動(dòng)增強(qiáng)的背景下是有利于地震發(fā)生的，因此區(qū)域的群體異常有更好的對(duì)應(yīng)效能；異常過(guò)程和從屬函數(shù)方法提取的水溫異常能較好地表達(dá)震前水溫變化加劇的特征，區(qū)域的群體異?？梢宰鳛樵颇系貐^(qū)有MS≥50地震的短期時(shí)間判據(jù)。

關(guān)鍵詞：云南水溫觀測(cè)；從屬函數(shù)方法；MS≥50地震；區(qū)域構(gòu)造活動(dòng)增強(qiáng)；群體異常

中圖分類(lèi)號(hào)：P315723文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-0666（2017）02-0233-08

0引言

[KG（0.15mm]地震預(yù)測(cè)目前仍處在一個(gè)經(jīng)驗(yàn)積累的探索階段，大量的震例總結(jié)、震前前兆觀測(cè)異常特征研究可以為認(rèn)識(shí)地震孕育過(guò)程提供更多的資料。在云南中等地震前兆異常特征研究方面，眾多學(xué)者作了的工作（付虹等，2007，2008；張立，2013；李瓊等，2014，2016），得到了一些比較有意義的結(jié)果，但這些研究大部分是基于原始數(shù)據(jù)異常識(shí)別和震例研究進(jìn)行的。從屬函數(shù)法主要以觀測(cè)資料的斜率和相關(guān)系數(shù)作為參量，能較好地表達(dá)地震孕育進(jìn)入短臨階段的漲落變化特征（林輝等，2004），因此適用于前兆觀測(cè)資料的短臨異常信息提取。部分地震學(xué)者（葉秀薇，2004；李杰等，2003；高小其等，2012；范雪芳等，2010）采用該方法提取震前短臨異常，取得了較好的結(jié)果。[KG）][JP]

水溫觀測(cè)的是不同含水層由于裂隙發(fā)育導(dǎo)致串層溫度混合的結(jié)果，高精度水溫觀測(cè)在地震前兆監(jiān)測(cè)中發(fā)揮著重要作用（付子忠，1988）。大量的觀測(cè)結(jié)果表明，水溫觀測(cè)具有較強(qiáng)的映震能力（車(chē)用太等，1999；陳沅俊等，1994a，b）。眾多地震學(xué)者對(duì)水溫與地震的關(guān)系進(jìn)行了大量研究（車(chē)用太等，2003；金眾范等，2001；李加穩(wěn)，孫文波，2004；田華等，1994；尤傳俠，1990；劉耀煒等，2008；鮮述東等，2011；楊從杰等，2008；張惠等，2013），得出水溫異常主要在地震前的短期階段出現(xiàn)，時(shí)間進(jìn)程上具有顯著的短臨異常特征。

本文采用從屬函數(shù)分析方法提取云南水溫前兆異常，以2009—2015年云南省內(nèi)的17組MS≥50地震作為研究樣本，在不知道未來(lái)地震震中位置情況下，從全時(shí)空的角度對(duì)水溫的短期異常及效能進(jìn)行探討。

1水溫資料及地震樣本選取

云南共有52項(xiàng)水溫觀測(cè)資料用于日常地震跟蹤監(jiān)視，[HJ2.4mm]資料情況大致可分為3類(lèi)：（1）觀測(cè)不穩(wěn)定，正常和異常動(dòng)態(tài)不清晰，沒(méi)有對(duì)應(yīng)地震的實(shí)例；（2）觀測(cè)穩(wěn)定，呈趨勢(shì)性上升或下降，但映震效果不好，在地震監(jiān)測(cè)預(yù)報(bào)中沒(méi)有出現(xiàn)過(guò)異常對(duì)應(yīng)地震的實(shí)例；（3）觀測(cè)穩(wěn)定，有異常對(duì)應(yīng)地震實(shí)例，且有部分異常多次對(duì)應(yīng)地震的事實(shí)。本研究選取了觀測(cè)資料相對(duì)連續(xù)、穩(wěn)定、有地震對(duì)應(yīng)實(shí)例的34項(xiàng)水溫資料，采用從屬函數(shù)計(jì)算方法分別提取每項(xiàng)觀測(cè)資料的異常，34項(xiàng)水溫測(cè)點(diǎn)分布如圖1所示。

2008年汶川大地震后，云南省地震局依托“十五”項(xiàng)目，投入了一批新的水溫觀測(cè)力量，這些觀測(cè)資料從2009年開(kāi)始相對(duì)穩(wěn)定?？紤]到資料的豐富性，選用2009年以來(lái)的水溫資料和2009—2015年間發(fā)生在云南省內(nèi)的21次17組（發(fā)生在同一個(gè)月內(nèi)的2次地震視為一組）MS≥50地震作為研究樣本，地震參數(shù)和分布見(jiàn)表1和圖1。

2異常提取方法

前兆觀測(cè)漲落變化幅度增大和出現(xiàn)加速，是地震孕育進(jìn)入短臨階段的主要特征（林輝等，2004），因此選取了能表征這些特征的數(shù)學(xué)方法來(lái)提取短臨異常。

21AP（Abnormal process）方法

[JP2]“AP”值是一階差分絕對(duì)值的滑動(dòng)累加值，通過(guò)過(guò)程累加可以同時(shí)達(dá)到有效抑制較長(zhǎng)周期、消除長(zhǎng)趨勢(shì)變化、使異常變化更突出的效果。該方法可以直接用來(lái)提取異常，也可以用來(lái)對(duì)資料進(jìn)行預(yù)處理（石紹先，曹刻，1999）。計(jì)算公式如下：[KH*1][JP]

[JP2] APi+T[KG-*4]=[KG-*2]∑[DD（]i+Ti[DD）][KG-*2][JB<1|]xi+1-xi[KG-*2][JB>1|]（i=1，2，3，…，n）[JY]（1）[KH*1D][JP]

其中：xi為某一測(cè)項(xiàng)的觀測(cè)值，T為累加過(guò)程的時(shí)間窗長(zhǎng)。本研究選取10天觀測(cè)數(shù)據(jù)作為時(shí)間窗長(zhǎng)，意在凸顯短期變化，為下一步采用從屬函數(shù)法提取異常做準(zhǔn)備。圖2為保山觀測(cè)值與AP計(jì)算值曲線(xiàn)，從圖中可以看出，原始曲線(xiàn)呈緩升變化，異常不突出，通過(guò)計(jì)算后，一些在原始觀測(cè)曲線(xiàn)上無(wú)法識(shí)別的加速異常被提取出來(lái)。

22從屬函數(shù)方法

前兆觀測(cè)量是時(shí)間的函數(shù)，前兆異常表現(xiàn)為觀測(cè)量隨時(shí)間的突出變化，雖然異常形態(tài)多種多樣，但其本質(zhì)都是觀測(cè)值曲線(xiàn)隨時(shí)間的斜率變化。模糊數(shù)學(xué)方法中的模糊從屬函數(shù)μi就是反映各種地震前兆觀測(cè)曲線(xiàn)隨時(shí)間的斜率變化的量（賀仲雄，1983）。相應(yīng)的解析式為：[KH*1]

μi=[KG-*4][JB<2[]1+[SX（]∝[JB<1|]ki [KG-*4][JB>1|][KG-*2][JB<1|]ri [KG-*4][JB>1|][SX）][JB>2]]-1（i=1，2，3，…，n）[JY]（2）[KH*1D]

式中：ki為觀測(cè)值變化斜率，表征的是數(shù)據(jù)變化的幅度；ri為滑動(dòng)平均值M（t）與時(shí)間t之間的相關(guān)系數(shù)，反映了計(jì)算序列內(nèi)在質(zhì)量的好壞；∝為經(jīng)驗(yàn)常數(shù)，不同觀測(cè)的∝值需要根據(jù)從屬函數(shù)μi與歷史震例的最佳對(duì)應(yīng)關(guān)系來(lái)確定，對(duì)同一個(gè)測(cè)項(xiàng)來(lái)說(shuō)，∝值越小，說(shuō)明測(cè)值越穩(wěn)定，通過(guò)計(jì)算得到的33個(gè)∝值均小于1，其中小于001的有22個(gè)，說(shuō)明采用從屬函數(shù)提取異常前先用“AP”方法進(jìn)行數(shù)據(jù)預(yù)處理是比較理想的；當(dāng)μi≥05時(shí)，可視為異常。即當(dāng)斜率增大到一定程度，可視為觀測(cè)值的漲落變化增大，發(fā)生地震的危險(xiǎn)性增大。圖3為使用從屬函數(shù)方法對(duì)AP值進(jìn)行計(jì)算后得到的μi值，從計(jì)算結(jié)果看，進(jìn)一步計(jì)算后，短期異常更為突出、容易識(shí)別。[FL）]

[SD30*1，1][JZ][XC李瓊3.TIF][KH*2/3D][HT5K][JZ]

圖3從屬函數(shù)計(jì)算值曲線(xiàn)[JZ]

Fig3Graph of calculated values by using the slave function[KH*1]

[BT（12+*4]3初步結(jié)果分析

31單項(xiàng)異常及預(yù)測(cè)效能分析[BT）]

對(duì)34項(xiàng)實(shí)測(cè)數(shù)據(jù)分別采用“AP”方法進(jìn)行計(jì)算，再以計(jì)算結(jié)果為基礎(chǔ)數(shù)據(jù)，進(jìn)行從屬函數(shù)計(jì)算后，將μ≥05的數(shù)據(jù)作為異常，對(duì)每項(xiàng)資料出現(xiàn)異常的起止時(shí)間進(jìn)行標(biāo)注，出現(xiàn)異常，規(guī)定異常結(jié)束3個(gè)月內(nèi)云南發(fā)生地震為對(duì)應(yīng)，否則視為虛報(bào)。研究時(shí)間段內(nèi)，用對(duì)應(yīng)地震的異常次數(shù)除以異?？偞螖?shù)表示對(duì)應(yīng)率，虛報(bào)地震次數(shù)除以地震總次數(shù)表示漏報(bào)率，按此原則得到了每個(gè)測(cè)項(xiàng)與地震預(yù)測(cè)的對(duì)應(yīng)效能（表2）。

從單項(xiàng)指標(biāo)分析，雖然地震對(duì)應(yīng)率參差不齊，但50%以上的異常都能對(duì)應(yīng)云南省內(nèi)的地震，表明提取的水溫短期異常有較好的對(duì)應(yīng)地震的能力，但同時(shí)對(duì)地震都有漏報(bào)。

32異常頻次與地震關(guān)系分析

對(duì)34項(xiàng)資料的從屬函數(shù)值進(jìn)行0-1化處理，μ≥05用“1”表示，μ<05用“0”表示，以旬為時(shí)間單位繪制頻次曲線(xiàn)圖（圖4）。規(guī)定頻次N≥6，并持續(xù)30 d以上視為異常，2009年以來(lái)共出現(xiàn)9組異常。

由圖4及表3可見(jiàn)，頻次異常一般持續(xù)時(shí)間為40～250 d，異常開(kāi)始時(shí)間超前地震20～230 d，如果每組異常過(guò)程只考慮首發(fā)地震，依據(jù)震例總結(jié)規(guī)范（張肇誠(chéng)等，2013），異常超前地震20～140 d全部為短期異常。8組異常對(duì)應(yīng)地震，1組虛報(bào)；異常超前地震20～140 d，17組地震報(bào)準(zhǔn)13組，漏報(bào)4組。如果異常過(guò)程中發(fā)震，則后續(xù)50～70 d內(nèi)多有中等地震繼續(xù)發(fā)生；如果異常結(jié)束后才發(fā)生地震，則之后再次發(fā)生中等地震的可能性小。

33權(quán)重集成異常與地震關(guān)系分析

各單項(xiàng)前兆，其過(guò)去的預(yù)報(bào)效能是不同的，它們有優(yōu)劣之分。用權(quán)重（w）大小對(duì)其優(yōu)劣或信度加以區(qū)分，對(duì)第i種前兆的預(yù)報(bào)結(jié)論yi配一個(gè)權(quán)重系數(shù)wi，好的前兆權(quán)重大，差的權(quán)重小。用下式求取權(quán)重集成預(yù)報(bào)值（陳立德，1993）：[KH*4]

P=[KG-*3]∑[DD（]ni=1[DD）]w1y1+w2y2+……+wnyn[JY]（3）[KH*2D]

各單項(xiàng)前兆的權(quán)重可以用下式求得：[KH*4]

wi=[HT10.][SX（][SX（]1第i種前兆虛報(bào)數(shù)+1[SX）]∑[DD（]ni=1[DD）][SX（]1第i種前兆虛報(bào)數(shù)+1[SX）][SX）][JY]（4）[KH*2D]

權(quán)重wi符合條件：[KH*2]

∑[DD（]ni=1[DD）]wi=w1+w2+……+wn=1[JY]（5）[KH*2]

根據(jù)各項(xiàng)異常與地震的對(duì)應(yīng)關(guān)系，由式（4）計(jì)算得到各單項(xiàng)水溫的權(quán)重（表2），從屬函數(shù)值μ≥05，則y為“1”，μ<05則為“0”，同樣以旬為時(shí)間單位繪制權(quán)重集成信度值曲線(xiàn)（圖5）。當(dāng)權(quán)重集成信度值大于02，并持續(xù)20 d以上為異常，2009年以來(lái)共出現(xiàn)9組異常。

由圖5及表4可見(jiàn)，9組異常持續(xù)時(shí)間20～200 d，異常開(kāi)始時(shí)間超前地震10～180 d，為短期異常。8組異常對(duì)應(yīng)地震發(fā)生，1組虛報(bào)；17組地震報(bào)準(zhǔn)12組，漏報(bào)5組。

34R值評(píng)分結(jié)果

采用許紹燮（1989）的R值評(píng)分法，對(duì)頻次疊加值和權(quán)重集成值出現(xiàn)異常預(yù)報(bào)地震的效能進(jìn)行評(píng)價(jià)。計(jì)算公式如下：[KH*5/4]

R=[HT10.][SX（]報(bào)準(zhǔn)次數(shù)地震總次數(shù)[SX）]-[SX（]預(yù)報(bào)占用時(shí)間研究占用時(shí)間[SX）] [JY]（6）[KH*5/4]

根據(jù)表3和表4結(jié)果，R頻次=[SX（]1317[SX）]-[SX（]8902 280[SX）]=03744≈0374，R權(quán)重=[SX（]1217[SX）]-[SX（]6802 280[SX）]=0407 6≈0408，置信水平為975%的隨機(jī)預(yù)報(bào)R0值分別為0318和0299。R頻次>0318，R權(quán)重>0299，異常頻度和權(quán)重集成信度值預(yù)測(cè)方法，均能通過(guò)R值評(píng)分檢驗(yàn)，說(shuō)明云南地區(qū)水溫用AP方法提取異常過(guò)程后，再用從屬函數(shù)提取異常的研究結(jié)果對(duì)省內(nèi)5級(jí)以上中等地震具有較好的時(shí)間預(yù)測(cè)意義。

35分區(qū)異常與異常區(qū)內(nèi)地震相關(guān)關(guān)系分析

為了追蹤未來(lái)地震的發(fā)震地點(diǎn)，采用更小區(qū)域的異常與該區(qū)地震對(duì)應(yīng)，根據(jù)地震活動(dòng)特點(diǎn)，將云南分為多個(gè)不同的區(qū)域，以旬為單位給出分區(qū)的水溫異常頻次（圖6），選擇發(fā)生在區(qū)內(nèi)的MS≥50地震作為研究樣本，采用最佳對(duì)應(yīng)原則，確定各區(qū)的異常判定閾值，其分區(qū)和地震對(duì)應(yīng)情況如下：

（1）滇東北地區(qū)：26°N以北，102°E以東的云南境內(nèi)地區(qū)，共有東川、會(huì)澤、宣威、昭通、漁洞5個(gè)測(cè)點(diǎn)，將至少有2個(gè)測(cè)項(xiàng)出現(xiàn)異常（N≥2）視為異常，共出現(xiàn)6組異常，對(duì)應(yīng)地震3次、虛報(bào)3次。

（2）滇中至滇西北地區(qū)：紅河斷裂以東、小江斷裂帶以西的云南境內(nèi)，共有麗江、永勝、姚安、元謀、華寧、江川等15個(gè)測(cè)點(diǎn)，異常頻次N≥4，且持續(xù)時(shí)間在20 d以上視為異常，共出現(xiàn)8組異常、2組對(duì)應(yīng)地震、6組虛報(bào)、漏報(bào)地震4組。

（3）小滇西至滇西南地區(qū)：瀾滄江以西的云南境內(nèi)，共有保山、騰沖、孟連等8個(gè)測(cè)點(diǎn)，將異常頻次N≥3，且持續(xù)時(shí)間在20 d以上視為異常，共出現(xiàn)6組異常、2組對(duì)應(yīng)地震、4組虛報(bào)；區(qū)內(nèi)共5組地震，對(duì)應(yīng)2組、漏報(bào)3組。[JP2]

（4）滇南至滇西南地區(qū)：25°N以南，995°E以東的云南境內(nèi)，共有華寧、峨山、鎮(zhèn)沅、大寨等11個(gè)測(cè)點(diǎn)，將異常頻次N≥3且持續(xù)時(shí)間在20 d以上視為異常，共出現(xiàn)6組異常、1組對(duì)應(yīng)地震、5組虛報(bào)；區(qū)內(nèi)共2組地震，對(duì)應(yīng)1組、漏報(bào)1組。[JP]

分區(qū)的結(jié)果顯示，如果把云南劃分為更小的構(gòu)造單元，用更小構(gòu)造單元的異常對(duì)應(yīng)本單元的地震，虛報(bào)、漏報(bào)率會(huì)大大增加，這一結(jié)果表明：地震不一定都發(fā)生在異常集中的地區(qū)，異常集中的小區(qū)域也不一定發(fā)生地震，這一觀測(cè)事實(shí)表明如果跟著異常走，地點(diǎn)錯(cuò)報(bào)的可能性會(huì)加大。分析認(rèn)為當(dāng)云南地區(qū)構(gòu)造活動(dòng)增強(qiáng)時(shí)，表現(xiàn)為大范圍多點(diǎn)的前兆觀測(cè)出現(xiàn)異常，以云南地區(qū)大區(qū)域?yàn)閱卧?，?dāng)異常增加的時(shí)候是有利于地震發(fā)生的。

4結(jié)論及討論

利用從屬函數(shù)方法提取單項(xiàng)水溫異常，分析地震預(yù)測(cè)效能，對(duì)單項(xiàng)異常進(jìn)行頻次疊加，提取綜合異常，分析與地震的關(guān)系。根據(jù)本文研究結(jié)果，得出如下幾點(diǎn)認(rèn)識(shí)：[JP2]

（1）云南目前上報(bào)數(shù)據(jù)用于會(huì)商的水溫觀測(cè)資料共有52項(xiàng)，其中挑選出34項(xiàng)相對(duì)較好的資料，表明云南水溫資料的利用率大致為65%，有35%的資料沒(méi)有應(yīng)用效能。34個(gè)測(cè)項(xiàng)資料中，通過(guò)采用AP和從屬函數(shù)計(jì)算方法提取異常，出現(xiàn)異常后，14個(gè)測(cè)項(xiàng)對(duì)應(yīng)率在70%～89%間，4個(gè)測(cè)項(xiàng)對(duì)應(yīng)率在60%～69%間；8個(gè)測(cè)項(xiàng)對(duì)應(yīng)率在50%～59%間；若全省有8個(gè)測(cè)項(xiàng)同時(shí)出現(xiàn)異常，其后云南省內(nèi)發(fā)生5級(jí)以上地震的對(duì)應(yīng)率為100%，表明水溫在震前短臨階段的漲落變化是存在的，井水溫對(duì)地震在時(shí)間上有相對(duì)較好的預(yù)測(cè)效能。[JP]

（2）全省水溫的異常頻次和權(quán)重集成方法對(duì)省內(nèi)MS≥5地震都有較好的預(yù)測(cè)能力，預(yù)測(cè)效能都能通過(guò)置信水平為975%的R值評(píng)分檢驗(yàn)。從2種方法的結(jié)果看，權(quán)重集成比直接用頻次預(yù)測(cè)，漏報(bào)多1次。說(shuō)明在云南的水溫觀測(cè)中，出現(xiàn)異常的數(shù)量更為重要，異常數(shù)量足矣刻畫(huà)云南的構(gòu)造活動(dòng)特征，并可用該特征進(jìn)行發(fā)震時(shí)間預(yù)測(cè)。

（3）如果把水溫測(cè)點(diǎn)按照構(gòu)造和地震活動(dòng)進(jìn)行分區(qū)，用區(qū)內(nèi)的異常對(duì)應(yīng)區(qū)內(nèi)的地震，虛報(bào)率和漏報(bào)率都大大增加，說(shuō)明在云南地區(qū)，采用水溫異常找未來(lái)的地震地點(diǎn)是困難的，水溫的異常變化更適合作短期尺度的時(shí)間預(yù)測(cè)。

（4）震前大范圍水溫異常增多的事實(shí)，與前兆觀測(cè)到的大量異常是區(qū)域構(gòu)造活動(dòng)增強(qiáng)的結(jié)果，異常與地震都是區(qū)域應(yīng)力水平增加的認(rèn)識(shí)是一致的（陳立德，付虹，2003），當(dāng)區(qū)域構(gòu)造活動(dòng)增強(qiáng)時(shí)，我們能觀測(cè)到整個(gè)區(qū)域的異常增多，這個(gè)結(jié)果為地震短期預(yù)測(cè)提供了可能。從2009年以來(lái)的水溫群體異常特征還發(fā)現(xiàn)，如果地震發(fā)生在異常結(jié)束以后，其后繼續(xù)發(fā)生中等地震的可能性小，反之則有繼續(xù)發(fā)生地震的危險(xiǎn)。表明前兆異常數(shù)量居高不下，可能反映的是區(qū)域應(yīng)力水平繼續(xù)維持高值，在這種狀態(tài)下有利于地震的發(fā)生，這與數(shù)值模擬得到的結(jié)果也是一致的（耿魯明等，1993），表明異常的持續(xù)發(fā)展可以作為判斷后續(xù)是否有地震的指標(biāo)。

（5）微觀前兆異常是監(jiān)視區(qū)域地殼應(yīng)力增強(qiáng)作用下地殼介質(zhì)產(chǎn)生應(yīng)變和應(yīng)力集中并達(dá)到一定量的結(jié)果和傳遞出的信息。水溫測(cè)值異?？杀灰暈楫?dāng)構(gòu)造活動(dòng)增強(qiáng)時(shí)，觀測(cè)站及附近地殼中含流體的介質(zhì)受力后其裂隙、孔隙、裂縫等發(fā)生變化，引起流體移動(dòng)速度及物質(zhì)交換發(fā)生變化的結(jié)果（陳立德等，2008）。

（6）當(dāng)區(qū)域構(gòu)造活動(dòng)增強(qiáng)，應(yīng)力水平發(fā)生變化時(shí)，水溫短期群體異常增多，與地震有較好的相關(guān)性，因此可以作為云南有震的時(shí)間判別指標(biāo)，但地點(diǎn)需要用其他方法進(jìn)行確定，如云南4級(jí)地震活動(dòng)增強(qiáng)的區(qū)域多為未來(lái)6級(jí)地震發(fā)生的地區(qū)（付虹等，2015）。如果每個(gè)觀測(cè)手段均發(fā)揮自己的優(yōu)勢(shì)，統(tǒng)籌考慮其作用，可能對(duì)提高地震短期預(yù)測(cè)水平更行之有效。

參考文獻(xiàn)：

車(chē)用太，黃輔瓊，王吉易，等1999張北地震地下流體異常變化[M]//中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司一九九八年張北地震北京：地震出版社，111-114

車(chē)用太，劉喜蘭，姚寶樹(shù)，等2003首都圈地區(qū)井水溫度的動(dòng)態(tài)類(lèi)型及其成因分析[J].地震地質(zhì)，25（3）：403-420

陳立德，付虹，鄔成棟2008強(qiáng)震短臨前兆異常共性特征的物理基礎(chǔ)[J].地震研究，31（2）：99-102

陳立德，付虹2003地震預(yù)報(bào)基礎(chǔ)與實(shí)踐[M].北京：地震出版社，182

陳立德1993地震綜合預(yù)報(bào)的權(quán)重集成法[J].地震研究，16（3）：267-271

陳沅俊，楊修信，劉永銘，等1994b首都圈精密地溫動(dòng)態(tài)觀測(cè)研究[M]//國(guó)家地震局科技監(jiān)測(cè)司首都圈地震短臨預(yù)報(bào)新方法觀測(cè)與研究北京：地震出版社，134-140

陳沅俊，楊修信，趙京海1994a一個(gè)地震前兆敏感點(diǎn)[M]//國(guó)家地震局科技監(jiān)測(cè)司首都圈地震短臨預(yù)報(bào)新方法觀測(cè)與研究北京：地震出版社，153-155

范雪芳，劉耀煒，吳桂娥，等2010華北地區(qū)水位與水氡中期、中短期前兆異常特征研究[J].地震研究，33（12）：147-158

付 虹，錢(qián)曉東，毛玉平，等20152014 年云南魯?shù)镸S65地震異常及預(yù)測(cè)[J].地震研究，38（2）：181-188

付虹，李永莉，趙小艷，等2008云南M≥5地震震前異常的統(tǒng)計(jì)特征[J].地震研究，31（4）：335-339

付虹，王世芹，秦嘉政，等20072007年5～6月滇西南地區(qū)M≥5成組地震活動(dòng)中、短臨異常特征及預(yù)測(cè)[J].地震研究，30（4）：34-310

付子忠1988地?zé)釀?dòng)態(tài)觀測(cè)與地?zé)崆罢譡J].地殼構(gòu)造與地殼應(yīng)力，1（1）：1-7

高小其，王瓊，朱成英，等20127級(jí)地震前地下流體前兆異常特征及判定指標(biāo)的研究[J].內(nèi)陸地震，26（2）：123-137

耿魯明，石耀霖，張國(guó)民1993地震孕育發(fā)生的場(chǎng)源關(guān)系初步研究[J].中國(guó)地震，9（4）：310-319

賀仲雄1983模糊數(shù)學(xué)及其應(yīng)用[M].天津：天津科學(xué)技術(shù)出版社

金眾范，王勇，盛國(guó)松，等2001遼寧省地?zé)嵊^測(cè)現(xiàn)狀及映震分析[J].東北地震研究，17（3）：25-32

李加穩(wěn)，孫文波2004邦臘掌溫泉水氡水溫地震異常特征研究[J].防災(zāi)技術(shù)高等專(zhuān)科學(xué)校學(xué)報(bào)，6（2）：25-29

李杰，韓海華，馬玉香，等2003數(shù)字化形變觀測(cè)資料異常識(shí)別方法的應(yīng)用及評(píng)價(jià)[J].東北地震研究，19（1）：25-33

李瓊，付虹，錢(qián)曉東2016云南盈江MS≥61地震前兆異常特征分析[J].地震研究，35（2）：155-160

李瓊，付虹，朱榮歡，等2014滇東北MS≥50地震前地下流體異常特征與魯?shù)?5級(jí)地震研究[J].地震研究，37（4）：495-502

林輝，張平，唐采，等2004云南地區(qū)強(qiáng)震前兆短期綜合異常特征研究[J].地震研究，24（4）：113-118

劉耀煒，孫小龍，王世芹，等2008井孔水溫異常與2007年寧洱64級(jí)地震關(guān)系分析[J].地震研究，31（4），347-353

石紹先，曹刻，1999前兆群體非均勻度異常與地震短期預(yù)報(bào)[J].地震，19（1）：59-64

田華，仲陽(yáng)，張之立1994白家疃地溫異常變化與地震的對(duì)應(yīng)關(guān)系[J].地震地磁觀測(cè)與研究，15（3）：9-14

鮮述東，劉子哲，董艷濤，等2011高七井水溫映震規(guī)律淺析[J].內(nèi)陸地震，25（4）：349-353

許紹燮1989地震預(yù)報(bào)能力的評(píng)分[C]//國(guó)家地震局監(jiān)測(cè)司地震預(yù)報(bào)方法實(shí)用化研究文集北京：學(xué)術(shù)書(shū)刊出版社，586-590

楊從杰，鄭江蓉，肖飛，等2008金湖06井水溫震前變化特征研究[J].地震地磁觀測(cè)與研究，29（5）：26-31

葉秀薇2004地下流體從屬函數(shù)異常與地震關(guān)系的初步研究[J].防災(zāi)減災(zāi)工程學(xué)報(bào)，24（2）：195-201

尤傳俠1990瀾滄—耿馬地震前后的地溫異常變化[J].地震研究，3（2）：196-202

張慧，顧申宜，李志雄，等2013單井多層位水溫微動(dòng)態(tài)初步研究[J].地震，33（1）：101-110

張立20132011年緬甸72級(jí)地震前云南地下流體異常特征分析[J].地震，33（1）：119-126

張肇誠(chéng)，陳棋福，鄭大林2013震例總結(jié)研究探討[M].北京：地震出版社，35-36[ZK）][HJ][FL）]

[STHZ][WT4HZ][JZ]Study on the Relationship Between Abnormal Water[JZ]Temperature and Earthquakes in Yunnan Area

[WT5B1][STBZ][JZ]LI Qiong1，F(xiàn)U Hong1，MAO Huiling1，ZHU Ronghuan2，HE Deqiang3

[WT5"B1X][JZ]（1 Earthquake Administration of Yunnan Province，Kunming 650224，Yunnan，China）[JZ]

（2 Earthquake Administration of Kunming Municipality，Kunming 650500，Yunnan，China）[JZ]

（3 Earthquake Administration of Zhaotong Municipality，Zhaotong 657000，Yunnan，China）

[WT5HZ][JZ]Abstract[WTB1]

This paper selects the observations data of a total of 34 continuous and stable water temperature，using the aomalous process and subordinate function analysis method to extract abnormal jointly and analyze the relationship between water temperature and the MS≥50 earthquakes The results show that the single anomaly extracted by this method has short-term implications and good relationship with earthquake，but increased the omission rate of earthquake；The groups anomalies in Yunnan Province extracted by anomaly frequency and weight integration method have obvious correspondence with MS≥50 earthquakes，which reduce the omission，and can pass the R-value inspection of a confidence level of 975%；After Yunnan region has been divided into different areas by tectonic and seismic activity，we extracted small area groups anomaliesOnly considering earthquake at the same area，then the omission and false rate of earthquake significantly increased in each area；Analyzes its reason，we found that the water temperature anomaly and earthquake are the products of increasing regional tectonic activities，and the water temperature shows a better characteristics of increasing regional tectonic activities，and the earthquake will occur easily based on increasing regional tectonic activities，therefore regional groups anomalies have more correspondent potency；The temperature abnormity extracted by using the anomalous process and subordinate function analysis method shows the characteristic of accelerated water temperature before earthquake，and regional groups anomalies can be used as short-term time indicators of the earthquakes（ MS≥50）in Yunnan region

[WTHZ]Keywords：[WTB1]water temperature observation in Yunnan；subordinate function analysis method； MS≥50 earthquakes；increasing of regional tectonic activities；groups anomalies[LM]