李璐 王寶善 侯金欣

中國(guó)地震局地球物理研究所（地震觀測(cè)與地球物理成像重點(diǎn)實(shí)驗(yàn)室），北京市民族大學(xué)南路5號(hào) 100081

0 引言

由于受背景噪聲水平、記錄儀器分布等多種因素的影響，許多微弱的地震信號(hào)被淹沒(méi)在噪聲中，無(wú)法直接觀測(cè)。另一方面，大地震發(fā)生后的短時(shí)間內(nèi)，由于受強(qiáng)震尾波的影響以及大量余震事件的發(fā)生，波形相互重疊，難以區(qū)分。如果能夠自動(dòng)識(shí)別出這些微弱的地震信號(hào)和遺漏的早期余震，則對(duì)于完善地震目錄、全面分析地震活動(dòng)性、監(jiān)測(cè)地震破裂區(qū)域的震后形變等具有重要意義。

目前，在大量地震自動(dòng)識(shí)別方法中，廣泛應(yīng)用的主要有2種：一種是工業(yè)界普遍使用的長(zhǎng)短時(shí)窗比（Short-Term-Average／Long-Term-Average，STA／LTA）方法（Stevenson，1976；Earle et al，1994）。該方法利用包絡(luò)線上短時(shí)窗內(nèi)信號(hào)的平均值與長(zhǎng)時(shí)窗內(nèi)信號(hào)平均值的比值來(lái)表征信號(hào)振幅的變化，其中，短時(shí)窗用于截取微弱的有用信號(hào)，長(zhǎng)時(shí)窗內(nèi)信號(hào)的平均值則反映背景噪聲水平。雖然算法簡(jiǎn)單，計(jì)算速度快，但對(duì)于信噪比過(guò)低的信號(hào)，這種方法并不適用。另一種方法是基于互相關(guān)技術(shù)的模板匹配濾波技術(shù)（Matched Filter Technique，MFT）（Gibbons et al，2006；Shelly et al，2007；Peng et al，2009）。該方法在滑動(dòng)窗互相關(guān)（Sliding-Window Cross-Correlation，SCC）檢測(cè)技術(shù)的基礎(chǔ)上發(fā)展起來(lái)（Yang et al，2009），是在信噪比較低的情況下提取微弱信號(hào)的一種有效方法。其中，Gibbons等（2006）分析了1997年8月16日Kara Sea地震主震后4hr發(fā)生的1個(gè)微地震，發(fā)現(xiàn)傳統(tǒng)的STA／LTA方法無(wú)法檢測(cè)出此次事件，而模板匹配濾波技術(shù)在單通道記錄上即可將其分辨出來(lái)。由此可見(jiàn)，模板匹配濾波技術(shù)在微地震的自動(dòng)檢測(cè)方面具有明顯優(yōu)勢(shì)。

除了地震學(xué)，模板匹配濾波技術(shù)還在其他領(lǐng)域得到廣泛應(yīng)用。2015年曾在科學(xué)界引起轟動(dòng)的重大發(fā)現(xiàn)——引力波的探測(cè)（Abbott et al，2016）中，科學(xué)家們利用干涉疊加提高探測(cè)器精度，于2015年9月14日09：50：45（UTC）探測(cè)到引力波信號(hào)，該信號(hào)極其微弱，應(yīng)變約為1.0×10－21。為了進(jìn)一步驗(yàn)證該微弱信號(hào)的來(lái)源，科學(xué)家們以黑洞模型模擬出理論引力波信號(hào)作為模板，在實(shí)際觀測(cè)的連續(xù)波形上掃描，發(fā)現(xiàn)了與模板波形一致的信號(hào)，且時(shí)間上與實(shí)際觀測(cè)到的引力波信號(hào)一致，從而驗(yàn)證了引力波信號(hào)的可靠性。Zhao等（2016）以爆炸事件波形為模板掃描連續(xù)波形，并未在2015年8月12日天津危險(xiǎn)品倉(cāng)庫(kù)爆炸前發(fā)現(xiàn)小爆炸事件。此外，在醫(yī)學(xué)、天體物理、電子通信等領(lǐng)域，模板匹配濾波技術(shù)還是用于造影、圖像識(shí)別、語(yǔ)音識(shí)別等的重要手段（Hoover et al，2000；Dong et al，2008；Avadhanulu et al，2013）。

本文主要介紹模板匹配濾波技術(shù)的基本原理和數(shù)據(jù)處理流程。從尋找大地震的前震與余震、遠(yuǎn)程動(dòng)態(tài)觸發(fā)地震、誘發(fā)地震和非火山震顫等方面，詳述該方法在地震數(shù)據(jù)處理中的應(yīng)用。此外，對(duì)于模板匹配濾波技術(shù)的改進(jìn)與發(fā)展也進(jìn)行了簡(jiǎn)要討論。

1 MFT基本原理

模板匹配濾波技術(shù)（MFT）利用現(xiàn)有地震目錄中記錄到的事件作為模板，掃描連續(xù)波形，識(shí)別出與模板事件波形具有一定相似度的地震事件。匹配濾波檢測(cè)微地震的處理流程如圖1所示。

圖1 模板匹配濾波技術(shù)數(shù)據(jù)處理流程

第1步，選取合適的頻段對(duì)模板事件和連續(xù)波形作同樣的濾波處理。微地震檢測(cè)一般為區(qū)域地震，因此，選取的頻段應(yīng)相對(duì)高頻。

第2步，標(biāo)定每個(gè)模板事件三分量波形上的P波、S波到時(shí)。標(biāo)定方法可以根據(jù)已知速度模型計(jì)算理論到時(shí)，也可人為手動(dòng)拾取。再對(duì)每一個(gè)模板事件的三分量分別計(jì)算P波、S波的信噪比。其中，信號(hào)項(xiàng)選取P波或S波到時(shí)后某一時(shí)間窗內(nèi)的波形計(jì)算均方根振幅；噪聲項(xiàng)選取在P波到時(shí)前同樣長(zhǎng)度的時(shí)間窗內(nèi)的波形計(jì)算均方根振幅。分別用P波或S波的信號(hào)項(xiàng)除以噪聲項(xiàng)，得到對(duì)應(yīng)的P波或S波信噪比。在后面的波形掃描中，只選取信噪比高于某一閾值的模板事件波形。

第3步，對(duì)于某一模板事件，利用下式計(jì)算模板事件的波形和對(duì)應(yīng)分量的連續(xù)波形的互相關(guān)系數(shù)（Cross-Correlation coefficient，CC系數(shù)））

其中，t0、t1分別為用于計(jì)算互相關(guān)系數(shù)的窗口的開(kāi)始與結(jié)束時(shí)刻；X（t）、Y（t）分別為 t0到 t1時(shí)間段內(nèi)的模板事件波形和連續(xù)波形。計(jì)算時(shí)每次在連續(xù)波形上移動(dòng)1個(gè)采樣點(diǎn)，最終得到該分量上的一道互相關(guān)波形。

第4步，根據(jù)模板事件上P波或S波到時(shí)與該事件發(fā)震時(shí)刻之間的到時(shí)差，將每道互相關(guān)波形平移至發(fā)震時(shí)刻，再對(duì)所有臺(tái)站、所有分量的互相關(guān)波形進(jìn)行疊加，取平均，得到該模板事件的平均互相關(guān)波形。

第5步，設(shè)定互相關(guān)系數(shù)的閾值，當(dāng)平均互相關(guān)波形上某一時(shí)刻的互相關(guān)系數(shù)高于該閾值時(shí)，就認(rèn)為這是一個(gè)新檢測(cè)事件。選取閾值時(shí)，一般選擇絕對(duì)中位差（Median Absolute Deviation，MAD）的 9倍（Peng et al，2009）或 12倍（Yao et al，2015）。但是，當(dāng)臺(tái)站分布稀疏、波形信噪比較低時(shí)，可以根據(jù)需要適當(dāng)提高閾值（如15倍MAD），以減少模板事件與噪聲信號(hào)間的錯(cuò)誤檢測(cè)。圖2為利用MFT以2015年4月26日20：08：33發(fā)生于中國(guó)藏南地區(qū)的1次2.40級(jí)地震為模板，檢測(cè)到發(fā)生于4月25日19：43：39的1次1.81級(jí)地震。如圖2所示，通過(guò)掃描在連續(xù)波形上可以識(shí)別出與模板事件波形相似、沒(méi)有在地震目錄上出現(xiàn)的地震事件?；谥貜?fù)地震的特性，認(rèn)為新檢測(cè)事件位置與對(duì)應(yīng)的模板事件位置一致。根據(jù)模板事件和新檢測(cè)事件在連續(xù)波形上的到時(shí)差，可以計(jì)算出新檢測(cè)事件的發(fā)震時(shí)刻。根據(jù)某一時(shí)間窗內(nèi)新檢測(cè)事件與模板事件的振幅比，可確定新檢測(cè)事件的震級(jí)。

2 M FT在地震學(xué)中的應(yīng)用

在算法逐漸完善的基礎(chǔ)上，模板匹配濾波技術(shù)已經(jīng)在多種地震事件的檢測(cè)中得到應(yīng)用。

2.1 運(yùn)用M FT尋找前震和余震

許多大地震發(fā)生前會(huì)有1個(gè)或多個(gè)前震事件發(fā)生，但前震與主震成核過(guò)程之間的聯(lián)系還存在爭(zhēng)議。建立更加完整的前震活動(dòng)序列，有助于加深對(duì)地震發(fā)生物理過(guò)程的理解，為地震預(yù)測(cè)提供指導(dǎo)（Jones et al，1979；Ben-zion，2008）。Kato等（2012）通過(guò)波形互相關(guān)的方法，在2011年日本東北9.0級(jí)地震之前識(shí)別出2個(gè)相互分離的前震地震群，并發(fā)現(xiàn)它們?cè)跁r(shí)間上有向主震初始破裂點(diǎn)緩慢移動(dòng)的趨勢(shì)。通過(guò)對(duì)滑動(dòng)速率的分析發(fā)現(xiàn)，第2個(gè)前震序列產(chǎn)生主要的應(yīng)力加載，促進(jìn)主震的破裂過(guò)程。Kato等（2014a）對(duì)2014年智利8.1級(jí)地震發(fā)生前后的連續(xù)波形進(jìn)行模板掃描，新檢測(cè)到數(shù)量10倍于原始地震目錄的地震事件，并發(fā)現(xiàn)前震序列中慢滑移事件同樣有向著主震成核點(diǎn)遷移的趨勢(shì)。另一方面，Wu等（2014）對(duì)1999年土耳其7.1級(jí)地震前65hr內(nèi)的連續(xù)波形進(jìn)行了匹配濾波分析，并未發(fā)現(xiàn)前震地震活動(dòng)對(duì)主震破裂的促進(jìn)作用，這說(shuō)明前震對(duì)主震成核的促進(jìn)作用以及對(duì)主震發(fā)震位置的影響并不是普遍現(xiàn)象。

大地震發(fā)生后的數(shù)小時(shí)內(nèi)，在距主震震中幾個(gè)破裂長(zhǎng)度的范圍內(nèi)，地震活動(dòng)性明顯升高，這些被觸發(fā)的地震稱(chēng)為余震。清楚認(rèn)識(shí)早期余震的活動(dòng)特性，是分析地震觸發(fā)機(jī)制的重要基礎(chǔ)（Helmstetter et al，2009）。模板匹配濾波技術(shù)有助于識(shí)別出那些淹沒(méi)在主震尾波中的早期余震事件，以完善余震序列。Peng等（2009）將該方法運(yùn)用到2004年帕克菲爾德6.0級(jí)地震發(fā)生后的3天內(nèi)，共檢測(cè)到數(shù)量11倍于地震目錄的余震事件（圖3）。檢測(cè)結(jié)果有效補(bǔ)充了早期余震序列，有助于清晰地刻畫(huà)出余震的遷移特性。

圖2 模板匹配濾波檢測(cè)地震示意圖

同樣運(yùn)用模板匹配濾波技術(shù)，Lengliné等（2012）在2011年日本東北9.0級(jí)地震之后的12hr內(nèi)檢測(cè)到數(shù)量12倍于原始地震目錄的事件，通過(guò)這些新檢測(cè)事件在時(shí)間和空間上的分布，得到余震衰減速率（p＝1.0）和余震區(qū)沿?cái)鄬幼呦虻臄U(kuò)展趨勢(shì)；Meng等（2013）通過(guò)匹配濾波后更加完善的地震目錄發(fā)現(xiàn)，2003年美國(guó)圣西蒙6.5級(jí)地震后，近場(chǎng)的觸發(fā)地震與由主震導(dǎo)致的靜態(tài)剪切應(yīng)力變化間呈正相關(guān)，說(shuō)明在圣安德烈亞斯斷層帕克菲爾德地區(qū)存在可能的低摩擦作用；Meng等（2014）通過(guò)對(duì)比2個(gè)距主震不同震中距的區(qū)域內(nèi)的地震活動(dòng)性認(rèn)為，靜態(tài)、動(dòng)態(tài)應(yīng)力的改變均對(duì)近場(chǎng)地震的觸發(fā)有影響，其中，靜態(tài)庫(kù)侖應(yīng)力可以作用于長(zhǎng)期地震活動(dòng)性的改變，動(dòng)態(tài)應(yīng)力則主要影響更大區(qū)域尺度內(nèi)的瞬時(shí)地震活動(dòng)性的改變；Kato等（2014b）在2007年日本能登半島6.7級(jí)地震發(fā)生后的32天內(nèi)，同樣觀測(cè)到新檢測(cè)的余震事件沿破裂斷層的走向擴(kuò)展，同時(shí)發(fā)現(xiàn)該擴(kuò)展趨勢(shì)呈階梯狀，認(rèn)為可能與主震發(fā)生后的抗震余滑（a seismic after slip）有關(guān)。

圖3 對(duì)2004年帕克菲爾德6.0級(jí)地震后余震序列的模板匹配濾波檢測(cè)結(jié)果（據(jù)Peng等（2009））

2.2 運(yùn)用M FT尋找遠(yuǎn)程動(dòng)態(tài)觸發(fā)地震

除了在近場(chǎng)帶來(lái)應(yīng)力改變并產(chǎn)生觸發(fā)地震之外，大地震還可以傳播到幾百至上千千米處，在遠(yuǎn)場(chǎng)觸發(fā)一些地方震或慢地震（Peng et al，2010；Hill et al，2015）。這種觸發(fā)多與大振幅面波引起的本地動(dòng)態(tài)應(yīng)力擾動(dòng)有關(guān)，因此被稱(chēng)為動(dòng)態(tài)觸發(fā)。選取合適頻段濾波后，人為手動(dòng)拾取地方震到時(shí)的方法工作量較大，不利于系統(tǒng)觀測(cè)。已有研究表明，火山活動(dòng)區(qū)由于巖漿入侵等流體運(yùn)動(dòng)，容易達(dá)到應(yīng)力的臨界狀態(tài)，是尋找動(dòng)態(tài)觸發(fā)現(xiàn)象的重點(diǎn)區(qū)域（Moran et al，2004；Prejean et al，2004）。運(yùn)用模板匹配濾波技術(shù)，在2011年日本東北9.0級(jí)地震的面波持續(xù)期間，Yukutake等（2013）在距主震約450km的箱根火山區(qū)域檢測(cè)到由主震觸發(fā)的地方震。另一方面，板塊交界區(qū)域由于需要密切監(jiān)測(cè)地震活動(dòng)，臺(tái)站分布密集，有利于動(dòng)態(tài)觸發(fā)現(xiàn)象的觀測(cè)，而對(duì)于內(nèi)陸板塊的動(dòng)態(tài)觸發(fā)研究則較少。Yao等（2015）通過(guò)模板匹配濾波發(fā)現(xiàn)，在西藏中南部由2004年蘇門(mén)答臘9.1級(jí)地震和2005年尼亞斯8.6級(jí)地震動(dòng)態(tài)觸發(fā)的地方震，并對(duì)2次大地震動(dòng)態(tài)觸發(fā)的持續(xù)時(shí)間和觸發(fā)區(qū)域進(jìn)行了對(duì)比。

2.3 運(yùn)用M FT尋找誘發(fā)地震

在遠(yuǎn)場(chǎng)觸發(fā)地震活動(dòng)中，由人為活動(dòng)導(dǎo)致的誘發(fā)地震是一個(gè)重要組成部分。其中，非常規(guī)天然氣能源——頁(yè)巖氣的勘探開(kāi)發(fā)工作最早起源于美國(guó)，并隨著水壓致裂等多項(xiàng)技術(shù)的成熟在多個(gè)國(guó)家和地區(qū)推廣使用。而在該技術(shù)迅速發(fā)展的過(guò)程中，與生產(chǎn)活動(dòng)和提高產(chǎn)量有關(guān)的廢水注入地下的處理方式是否直接導(dǎo)致了當(dāng)?shù)氐卣鸹顒?dòng)性的增加，一直是科學(xué)界和工業(yè)界爭(zhēng)論的話題（Ellsworth，2013）。在實(shí)際觀測(cè)中，注水帶來(lái)的應(yīng)力擾動(dòng)十分微弱（約1MPa）（Zoback et al，1997），可觸發(fā)的地震震級(jí)較小，不易觀測(cè)。運(yùn)用模板匹配濾波技術(shù)可以實(shí)現(xiàn)對(duì)當(dāng)?shù)匚⑿〉卣鸬母尤娴奶綔y(cè)，從而有助于判斷地震活動(dòng)性變化與注水活動(dòng)之間的關(guān)聯(lián)。van der Hilst等（2013）利用模板匹配濾波技術(shù)尋找美國(guó)中西部地區(qū)的觸發(fā)地震，發(fā)現(xiàn)有廢水注入的地區(qū)，巖石流體孔隙壓增加，斷層系統(tǒng)更接近臨界狀態(tài)，當(dāng)有遠(yuǎn)場(chǎng)的大地震發(fā)生時(shí)，更容易在應(yīng)力擾動(dòng)下產(chǎn)生觸發(fā)地震。Walter等（2016）同樣利用模板匹配濾波技術(shù)進(jìn)一步研究了美國(guó)得克薩斯州東部和路易斯安那州西北部的地震活動(dòng)性，新檢測(cè)到的40個(gè)地震事件在空間上可分為幾個(gè)地震群，雖然沒(méi)有強(qiáng)烈的證據(jù)表明當(dāng)?shù)氐卣鸹顒?dòng)增加與廢水注入間有直接關(guān)系，但也并不能完全否定兩者間的相關(guān)性，作者認(rèn)為該地區(qū)需要更多的地震學(xué)監(jiān)測(cè)。在我國(guó)，也有一些研究工作將模板匹配濾波技術(shù)應(yīng)用于微小誘發(fā)地震事件的檢測(cè)中。Wang等（2015）在北京房山地區(qū)于2011年日本東北9.0級(jí)地震和2012年印度洋8.6級(jí)地震發(fā)生后的短時(shí)間內(nèi)，分別新檢測(cè)到分布集中的大量淺部地震，并認(rèn)為可能與當(dāng)?shù)夭傻V活動(dòng)帶來(lái)的地下應(yīng)力變化有關(guān)。

2.4 運(yùn)用M FT尋找非火山震顫

非火山震顫（Non-Volcanic Tremor，NVT）主要發(fā)生在俯沖帶內(nèi)，由許多相互重疊的低頻事件（Low Frequency Event，LFEs）組成（Shelly et al，2007、2011），具有持續(xù)時(shí)間長(zhǎng)、頻率較低、體波震相不清晰等特點(diǎn)（Beroza et al，2011），其信噪比較低，不易觀測(cè)。Shelly等（2009）在美國(guó)帕克菲爾德地區(qū)第1次利用波形互相關(guān)觀測(cè)到LFEs，并對(duì)這一組LFEs組成的震顫進(jìn)行定位，發(fā)現(xiàn)圣安德烈亞斯斷層可能延伸至地殼底部，比震源深度最大的普通地震還要深約10km。利用模板匹配濾波在震顫中發(fā)現(xiàn)更多重復(fù)信號(hào)（LFEs），有助于更好地刻畫(huà)俯沖帶內(nèi)斷層的破裂方向與可能的觸發(fā)機(jī)制（Aguiar et al，2009；Brown et al，2010）。

3 M FT方法的改進(jìn)與發(fā)展

在模板匹配濾波技術(shù)的基本原理之上，一些研究工作對(duì)其算法進(jìn)行了優(yōu)化。首先，模板匹配濾波是在連續(xù)波形上逐個(gè)采樣點(diǎn)掃描，運(yùn)算量大，耗時(shí)長(zhǎng)。因此，可以通過(guò)降采樣的方式，適當(dāng)縮短掃描的時(shí)間。此外，Meng等（2012）運(yùn)用圖形加速卡（Graphic Processing Unit，GPU）的并行計(jì)算方法，將連續(xù)波形分成幾個(gè)區(qū)塊，每個(gè)區(qū)塊在GPU的不同線程下同時(shí)進(jìn)行匹配濾波掃描。若同時(shí)啟用N個(gè)線程，則計(jì)算速度可以提高到原來(lái)的N倍。其次，前人研究工作中，對(duì)模板事件波形的垂直分量和水平分量分別選取P波、S波到時(shí)附近的時(shí)間窗（Peng et al，2009），當(dāng)研究區(qū)域內(nèi)臺(tái)站分布稀疏時(shí)，可以如圖2所示，對(duì)模板事件波形的每一分量同時(shí)選取P波和S波附近的時(shí)間窗，則模板事件波形的數(shù)量翻倍，更有利于互相關(guān)波形疊加時(shí)對(duì)不相關(guān)信號(hào)的壓制。最后，與重復(fù)地震的特征類(lèi)似，模板匹配濾波的前提假設(shè)是，新檢測(cè)到的事件與模板事件的發(fā)震位置一致。然而，在實(shí)際觀測(cè)中，有些遺漏事件可能具有與模板事件高度相似的波形，但位置上略有偏差，這導(dǎo)致它們傳至臺(tái)站的用時(shí)略有不同，此時(shí)，如果將各臺(tái)站的互相關(guān)波形按模板事件的到時(shí)平移至發(fā)震時(shí)刻，則會(huì)導(dǎo)致該遺漏事件的互相關(guān)信號(hào)沒(méi)有完全對(duì)齊，此時(shí)的疊加波形上互相關(guān)系數(shù)減小，降低了該事件被檢測(cè)到的可能性。Zhang等（2015a）發(fā)展了一種Match and Locate（M＆L）改進(jìn)算法，在模板事件附近的1個(gè)小三維區(qū)域內(nèi)進(jìn)行網(wǎng)格搜索，將新檢測(cè)事件定位在每一個(gè)網(wǎng)格點(diǎn)上，計(jì)算此時(shí)該事件和對(duì)應(yīng)的模板事件之間的到時(shí)差，將各道互相關(guān)波形按對(duì)應(yīng)的到時(shí)差平移后再疊加，最后比較所有網(wǎng)格點(diǎn)的疊加互相關(guān)系數(shù)和信噪比，以確定事件的位置。由于是在小范圍內(nèi)進(jìn)行空間搜索，因此，對(duì)速度模型的依賴(lài)性較小。該方法不僅實(shí)現(xiàn)了對(duì)新檢測(cè)事件的重定位，而且提高了被檢測(cè)事件的互相關(guān)系數(shù)，從而增加了更小震級(jí)事件被探測(cè)到的可能性。Zhang等（2015c）利用此方法監(jiān)測(cè)到朝鮮地區(qū)的多次核爆實(shí)驗(yàn)，因而將全球的探測(cè)核實(shí)驗(yàn)?zāi)芰μ岣叩絿嵉募?jí)別。運(yùn)用同樣的方法，Zhang等（2015b）在日本Ontake火山2007、2014年2次噴發(fā)前均探測(cè)到密集的微震活動(dòng)，為火山活動(dòng)的預(yù)測(cè)提供了依據(jù)。

4 結(jié)語(yǔ)

通過(guò)模板匹配濾波技術(shù)可以在背景噪聲水平較高時(shí)利用互相關(guān)掃描的方法尋找到肉眼難以觀測(cè)到的信號(hào)，該方法適用于各類(lèi)微弱信號(hào)的檢測(cè)。在地震學(xué)中，運(yùn)用該方法監(jiān)測(cè)地震活動(dòng)性對(duì)于斷層走向的刻畫(huà)、地震破裂過(guò)程模擬、地下速度結(jié)構(gòu)成像等方面的研究具有重要意義?；谀０迤ヅ錇V波技術(shù)的自動(dòng)掃描、識(shí)別信號(hào)的特征，今后工作中，可將該技術(shù)應(yīng)用于超密集臺(tái)陣的自動(dòng)化處理，以提高地震監(jiān)測(cè)工作的效率。于2015年全面啟動(dòng)的國(guó)家地震烈度速報(bào)與預(yù)警工程，將在全國(guó)范圍內(nèi)建成由5000余個(gè)強(qiáng)震臺(tái)站組成的密集觀測(cè)臺(tái)網(wǎng)用于災(zāi)情的快速分析，并指導(dǎo)抗震救災(zāi)工作。利用模板匹配濾波技術(shù)，將地震事件的自動(dòng)識(shí)別與定位相結(jié)合將是處理這些海量地震數(shù)據(jù)的有效手段。