本發(fā)明涉及地震監(jiān)測(cè)和地震預(yù)測(cè)領(lǐng)域，特別涉及一種快速獲取輸入地震事件的地震參數(shù)估計(jì)的方法和系統(tǒng)以及地震實(shí)時(shí)搜索引擎。

背景技術(shù)：
地震事件是地球內(nèi)部在某個(gè)局部區(qū)域產(chǎn)生的介質(zhì)構(gòu)造突變而發(fā)生的，包括工業(yè)活動(dòng)引發(fā)的微震事件，為增加石油與天然氣的產(chǎn)量而采用注水破裂引起的微震事件，強(qiáng)度較小或者較大的自然地震事件，核爆炸事件或者是炸藥爆炸事件。當(dāng)一個(gè)地震事件產(chǎn)生時(shí)，該地震事件對(duì)應(yīng)的地震波形信號(hào)數(shù)據(jù)會(huì)被地震監(jiān)測(cè)臺(tái)站記錄下來(lái)(地震波接收儀器)，我們想要盡可能快的推斷出該地震事件發(fā)生的位置，震級(jí)以及震源機(jī)制等信息。目前該領(lǐng)域最為先進(jìn)的技術(shù)能夠利用縱波和/或橫波的傳播時(shí)間信息在地震發(fā)生后馬上推算出地震發(fā)生的位置，但是其震源機(jī)制信息卻不得不花費(fèi)數(shù)小時(shí)或者是數(shù)天的時(shí)間才能夠通過(guò)分析得到。對(duì)于當(dāng)前地震定位技術(shù)和地震震級(jí)判斷技術(shù)來(lái)說(shuō)，仍然需要一定水平的專家來(lái)進(jìn)行推算才能保證估計(jì)結(jié)果的準(zhǔn)確程度。以發(fā)生在2004年12月26日蘇門答臘北部西海岸的9.1級(jí)的地震事件為例，因?yàn)樵摰卣鹗录l(fā)生時(shí)是圣誕節(jié)假期的緣故，所以該地震事件沒(méi)有被及時(shí)的報(bào)告，海嘯的警報(bào)沒(méi)有發(fā)出。實(shí)時(shí)地報(bào)告發(fā)生的地震事件的震源機(jī)制是很重要的。如Bouchon，M.，Karabulut，H.，Aktar，M.，Ozalaybey，S.，Schmittbuhl，J.，andBouin，M-P.(2011).Extendednucleationofthe1999Mw7.6IzmitEarthquake，Sciences331，877-88分析了1999年土耳其伊茲米特發(fā)生的7.6級(jí)地震的震前征兆，發(fā)現(xiàn)一系列前震的震源機(jī)制都顯示相似的斷層錯(cuò)動(dòng)，直至加速到動(dòng)態(tài)斷裂。實(shí)時(shí)地獲得前震的震源機(jī)制能夠幫助我們開(kāi)發(fā)出 前期預(yù)警系統(tǒng)或者更進(jìn)一步，甚至是預(yù)測(cè)地震的發(fā)生。然而，獲得地震震源機(jī)制是需要一些諸如波形模擬和波形反演等過(guò)程的處理，這一處理過(guò)程通常在地震發(fā)生之后的若干小時(shí)或者數(shù)天之內(nèi)進(jìn)行，并且這一過(guò)程通常要求使用高端的模擬與計(jì)算技術(shù)。需要一種快速獲取發(fā)生的地震事件的地震參數(shù)信息的方法。

技術(shù)實(shí)現(xiàn)要素：
根據(jù)本發(fā)明實(shí)施例，提出了一種獲取地震事件的地震參數(shù)估計(jì)的方法，包括：輸入地震事件的地震圖數(shù)據(jù)，地震圖數(shù)據(jù)包括地震波形；在預(yù)先建立的歷史地震圖數(shù)據(jù)庫(kù)和理論地震圖數(shù)據(jù)庫(kù)中通過(guò)波形匹配，采用近似最近鄰搜索方法進(jìn)行并行搜索，按照預(yù)定的相似度條件找出與該輸入地震事件的地震圖數(shù)據(jù)相似的多個(gè)地震圖數(shù)據(jù)，作為搜索結(jié)果集合，其中所述歷史地震圖數(shù)據(jù)庫(kù)存儲(chǔ)實(shí)際發(fā)生的地震事件的地震圖數(shù)據(jù)，所述理論地震圖數(shù)據(jù)庫(kù)存儲(chǔ)人工模擬的地震圖數(shù)據(jù)；以及從搜索結(jié)果集合中確定匹配地震圖數(shù)據(jù)，并從中獲取輸入地震事件的地震參數(shù)估計(jì)。根據(jù)本發(fā)明實(shí)施例，提出了一種獲取地震事件的地震參數(shù)估計(jì)的系統(tǒng)，包括：歷史地震圖數(shù)據(jù)庫(kù)，存儲(chǔ)實(shí)際發(fā)生的地震事件的地震圖數(shù)據(jù)；理論地震圖數(shù)據(jù)庫(kù)，存儲(chǔ)人工模擬的地震圖數(shù)據(jù)；輸入裝置，輸入地震事件的地震圖數(shù)據(jù)，地震圖數(shù)據(jù)包括地震波形；搜索引擎，在歷史地震圖數(shù)據(jù)庫(kù)和理論地震圖數(shù)據(jù)庫(kù)中通過(guò)波形匹配，采用近似最近鄰方法進(jìn)行并行搜索，按照預(yù)定的相似度條件找出與該輸入地震事件的地震圖數(shù)據(jù)相似的多個(gè)地震圖數(shù)據(jù)，作為搜索結(jié)果集合，以及從搜索結(jié)果集合中確定匹配地震圖數(shù)據(jù)，并從中獲取輸入地震事件的地震參數(shù)估計(jì)。根據(jù)本發(fā)明實(shí)施例，提出了一種地震圖搜索引擎，包括：搜索裝置，針對(duì)輸入的地震事件的包括地震波形的地震圖數(shù)據(jù)，在預(yù)先建立的歷史地震圖數(shù)據(jù)庫(kù)和理論地震圖數(shù)據(jù)庫(kù)中通過(guò)波形匹配，采用近似最近鄰方法進(jìn)行并行搜索，按照預(yù)定的相似度條件找出與該輸入地震事件的地震圖數(shù)據(jù)相似的多個(gè)地震圖數(shù)據(jù)，作為搜索結(jié)果集合，其中所述歷史地震圖數(shù)據(jù)庫(kù)存儲(chǔ)實(shí)際發(fā)生的地震事件的地震圖數(shù)據(jù)，所述理論地震圖數(shù)據(jù)庫(kù)存儲(chǔ)人工模擬的地震圖數(shù)據(jù)；以及確定裝置，在搜索結(jié)果集合中，確定匹配地震圖數(shù)據(jù)，并從中獲取輸入地震事件的地震參數(shù)估計(jì)。附圖說(shuō)明通過(guò)下面結(jié)合附圖說(shuō)明本發(fā)明的示例實(shí)施例，將使本發(fā)明的上述及其它目的、特征和優(yōu)點(diǎn)更加清楚，其中：圖1示出了三分量地震圖數(shù)據(jù)的示例；圖2示出了根據(jù)本發(fā)明實(shí)施例的獲取地震事件的地震參數(shù)估計(jì)的系統(tǒng)的示意框圖；圖3示出了建立歷史地震圖數(shù)據(jù)庫(kù)的過(guò)程的示意流程圖；圖4示出了建立理論地震圖數(shù)據(jù)庫(kù)的過(guò)程的示意流程圖；圖5示出了具有多個(gè)時(shí)間維度的地震圖數(shù)據(jù)的示例；圖6示出了建立多棵隨機(jī)KD樹(shù)(MRKD-Tree)結(jié)構(gòu)的過(guò)程的示意圖；圖7示出了利用多棵隨機(jī)KD樹(shù)第一輪搜索過(guò)程；圖8示出了利用多棵隨機(jī)KD樹(shù)第二輪搜索過(guò)程；圖9示出了根據(jù)本發(fā)明實(shí)施例的地震事件搜索引擎的示意框圖；圖10示出了根據(jù)本發(fā)明實(shí)施例的獲取地震事件的地震參數(shù)估計(jì)的方法的流程圖；圖11示出了根據(jù)本發(fā)明實(shí)施例的搜索結(jié)果的示例；圖12示出了三種方法的地震圖搜索時(shí)間的比較：精確線性搜索(LinearSearch)，局部敏感哈希(LSH)方法，以及MRKD-Tree方法；以及圖13示出了輸入地震事件和使用線性搜索(LinearSearch)和MRKD-Tree方法分別得到的最佳匹配地震事件之間的最大互相關(guān)系數(shù)的比較。具體實(shí)施方式以下參照附圖，對(duì)本發(fā)明的示例實(shí)施例進(jìn)行詳細(xì)描述，本發(fā)明不限于下述示例實(shí)施例。為了清楚描述本發(fā)明的基本思想，附圖中僅示出了與本發(fā)明的技術(shù)方案密切相關(guān)的部件、功能或步驟，并且以下描述中省略了對(duì)已知技術(shù)、功能、部件或步驟的具體描述。地震事件是地球內(nèi)部在某個(gè)局部區(qū)域產(chǎn)生的介質(zhì)構(gòu)造突變而發(fā)生的。地震事件是由地下的巖石構(gòu)造瞬間破裂所造成的。地震破裂機(jī)制可以用幾個(gè)參數(shù)描述，例如破裂面的走向角度，傾角角度，滑動(dòng)角度，震級(jí)和震源深度?？紤]這些參數(shù)的結(jié)合，以及地震波記錄的震中距離，其相應(yīng)的地震圖有獨(dú)特的波形記錄，如橫波，縱波，面波都有不同的特征。地震圖是一個(gè)時(shí)間信號(hào)序列。對(duì)于給定的地球地質(zhì)構(gòu)造，每一個(gè)地震圖對(duì)應(yīng)著一個(gè)特定震源機(jī)制參數(shù)和震中距的組合。所以從地震圖可以反推震源機(jī)制參數(shù)和震中距。從地震臺(tái)或地震儀記錄到的是地震波傳播到記錄地點(diǎn)引起的介質(zhì)位移，或位移的速度(一階時(shí)間導(dǎo)數(shù))，或位移的加速度(二階時(shí)間導(dǎo)數(shù))。一般地震臺(tái)可采用三分量(垂直，東西，南北方向震動(dòng)記錄)地震圖數(shù)據(jù)，或一分量(垂直分量)地震圖數(shù)據(jù)。圖1顯示的是三分量地震圖數(shù)據(jù)的示意性示例，其中上面的圖表示南北振動(dòng)分量，中間的圖表示東西振動(dòng)分量，下面的圖表示垂直振動(dòng)分量，橫軸為相對(duì)時(shí)間(秒)，縱軸為地震儀所在處地球介質(zhì)振動(dòng)的速度。在地震學(xué)領(lǐng)域，地震模擬中存在一個(gè)基本假設(shè)：如果同一個(gè)地震監(jiān)測(cè)臺(tái)站記錄的兩個(gè)地震圖是相同的，那么這兩個(gè)地震圖對(duì)應(yīng)的震源解以及它們傳播過(guò)程的速度結(jié)構(gòu)也應(yīng)該是一樣的。因此，在地震學(xué)領(lǐng)域，采用正演和反演方法來(lái)查找波形匹配的地震圖是估計(jì)地震參數(shù)和地球結(jié)構(gòu)的標(biāo)準(zhǔn)方法，例如可以參見(jiàn)Xu，Y.，Herrmann，R.B.，andKeith，D.K.(2010).Sourceparametersofregionalsnall-to-moderate earthquakesintheYunnan-SichuanregionofChina，Bull.Seismol.Soc.Am.100，No.5B，pp.2518-2531，以及Bonner，J.，Herrmann.R.，andHarley，B.(2010).Variable-periodsurface-wavemagnitudes：arapidandrobustestimatorofseismicmoments.Bull.Seismol.Soc.Am.100，No.5A，pp.2301-2309。本發(fā)明實(shí)施例對(duì)于發(fā)生的地震事件，通過(guò)快速搜索歷史記錄的地震圖數(shù)據(jù)和人工模擬的地震圖數(shù)據(jù)，通過(guò)波形匹配找出與該地震事件的地震圖數(shù)據(jù)相似的地震圖數(shù)據(jù)，并從中獲取該地震事件的地震參數(shù)估計(jì)。根據(jù)本發(fā)明實(shí)施例，歷史記錄和人工模擬的地震圖數(shù)據(jù)是預(yù)先準(zhǔn)備好的。此外，該搜索方法不僅僅是找出一個(gè)地震波形的解，還可以找出在一定相似度條件下不同的地震事件的波形的解集。根據(jù)本發(fā)明實(shí)施例，對(duì)于每個(gè)地震事件的波形的解就像是網(wǎng)絡(luò)搜索引擎的結(jié)果中的每個(gè)網(wǎng)頁(yè)鏈接一樣，可以按相似度排序。這樣使得有機(jī)會(huì)檢測(cè)地震模擬中非唯一的地震事件的波形的解集。根據(jù)本發(fā)明實(shí)施例，采用近似最近鄰搜索方法，能夠在地震發(fā)生后幾秒鐘內(nèi)，快速獲得這些非唯一的地震事件的波形的解集，從而在極短的時(shí)間內(nèi)對(duì)發(fā)生的地震事件有全面的認(rèn)識(shí)。考慮到地震信號(hào)噪聲的影響，模型參數(shù)的相互作用，以及地球結(jié)構(gòu)的復(fù)雜性，需要獲取多個(gè)非唯一的解，并且可以從這些解的特征評(píng)估搜索和解集的可靠性和解析度。圖2示出了根據(jù)本發(fā)明實(shí)施例的獲取地震事件的地震參數(shù)估計(jì)的系統(tǒng)的示意框圖。該系統(tǒng)10包括：歷史地震圖數(shù)據(jù)庫(kù)110，存儲(chǔ)實(shí)際發(fā)生的地震事件的地震圖數(shù)據(jù)；理論地震圖數(shù)據(jù)庫(kù)120，存儲(chǔ)人工模擬的地震圖數(shù)據(jù)；輸入裝置130，輸入地震事件的地震圖數(shù)據(jù)，地震圖數(shù)據(jù)包括地震波形；搜索引擎140，在歷史地震圖數(shù)據(jù)庫(kù)110和理論地震圖數(shù)據(jù)庫(kù)120中通過(guò)波形匹配，采用近似最近鄰方法進(jìn)行并行搜索，按照預(yù)定的相似度條件找出與該輸入地震事件的地震圖數(shù)據(jù)相似的多個(gè)地震圖數(shù)據(jù)，作為搜索結(jié)果集合，以及從搜索結(jié)果集合中確定匹配地震圖數(shù)據(jù)，并從中獲取輸入地震事件的地震參數(shù)估計(jì)。歷史和理論地震圖數(shù)據(jù)庫(kù)可以是本地或遠(yuǎn)程存儲(chǔ)的數(shù)據(jù)庫(kù)，也可以是集中式或分布式的數(shù)據(jù)。根據(jù)本發(fā)明實(shí)施例，系統(tǒng)10還可以包括評(píng)估裝置(未示出)，通過(guò)比較所有匹配地震圖數(shù)據(jù)的互相關(guān)系數(shù)分布，來(lái)評(píng)估搜索結(jié)果的解析度和可靠性。地震圖數(shù)據(jù)可以是如圖1所示的三分量波形信號(hào)，也可以采用任何其他適合的形式，例如一分量垂直波形信號(hào)。地震圖數(shù)據(jù)可以是數(shù)字化記錄的，采用如標(biāo)準(zhǔn)的SAC，SEG2，SEGD或SEGY格式等存儲(chǔ)，在數(shù)據(jù)文件頭存儲(chǔ)國(guó)際地震組織認(rèn)定的地震參數(shù)，地震參數(shù)例如包括震中距、震源深度、震源機(jī)制和震級(jí)，其中震源機(jī)制可以包括走向角度、傾角角度和滑動(dòng)角度。歷史地震圖數(shù)據(jù)庫(kù)110歷史地震圖數(shù)據(jù)庫(kù)110是預(yù)先建立的。為了建立歷史地震圖數(shù)據(jù)庫(kù)110，可以利用單個(gè)或多個(gè)不同地震監(jiān)測(cè)臺(tái)站記錄的地震事件的數(shù)字化記錄。例如，可以針對(duì)某個(gè)特定的地區(qū)以及地震監(jiān)測(cè)臺(tái)網(wǎng)建立數(shù)據(jù)庫(kù)，或?qū)τ诖蟪叨鹊膮^(qū)域性以及超遠(yuǎn)地震距離的地震問(wèn)題，根據(jù)震中距、震源深度和震源機(jī)制建立數(shù)據(jù)庫(kù)。根據(jù)本發(fā)明實(shí)施例，可以在數(shù)據(jù)庫(kù)中分開(kāi)存儲(chǔ)每個(gè)地震監(jiān)測(cè)臺(tái)的數(shù)字化記錄。可以針對(duì)每個(gè)地震臺(tái)獲取地震參數(shù)估計(jì)，并將多個(gè)地震臺(tái)地震參數(shù)估計(jì)的交集確定為輸入地震事件的地震參數(shù)估計(jì)。根據(jù)本發(fā)明實(shí)施例，歷史地震圖數(shù)據(jù)庫(kù)中的地震圖數(shù)據(jù)是經(jīng)過(guò)索引處理的，可以采用局部敏感哈希(LSH)方法進(jìn)行搜索或者采用MRKD-Tree方法進(jìn)行搜索。LSH是計(jì)算機(jī)科學(xué)領(lǐng)域的一種應(yīng)用于文檔檢索和圖片檢索的常用的快速的搜索方法，在此省略對(duì)其的詳細(xì)描述。詳細(xì)方法描述參見(jiàn)Slaney，M.，andCasey，M.(2008)，Locality-sensitivehashingforfindingnearestneighbors，IEEESingalProcessingMagzine，March，pp128-131，圖3示出了建立歷史地震圖數(shù)據(jù)庫(kù)的一個(gè)示例的示意流程圖300，在該示例中，應(yīng)用MRKD-Tree方法，對(duì)地震圖數(shù)據(jù)進(jìn)行索引。如圖3所示，在從單個(gè)或多個(gè)不同地震監(jiān)測(cè)臺(tái)站獲得記錄的地震圖數(shù)據(jù)后，在步驟302，進(jìn)行信號(hào)處理、去噪，拾取P波初至?xí)r間。在步驟304，對(duì)地震圖數(shù)據(jù)進(jìn)行分段頻率濾波以分成多個(gè)頻率分段，例如如選擇頻率分段[0.01Hz，0.04Hz]，[0.03Hz，0.06Hz]，并拾取濾波后的面波的最大振幅到時(shí)。這 可以例如去掉不同臺(tái)站數(shù)據(jù)中的儀器響應(yīng)(參見(jiàn)Herrmann，R.(1973).Someaspectsofband-passfilteringofsurfacewaves.Bull.Seismol.Soc.Am.63，No.2，pp.663-671)。在步驟306，針對(duì)每個(gè)頻率分段，分別根據(jù)地震圖數(shù)據(jù)的初至到時(shí)和最大振幅到時(shí)來(lái)對(duì)齊波形，截取得到兩組地震圖數(shù)據(jù)。在步驟308，分別對(duì)得到的兩組地震圖數(shù)據(jù)應(yīng)用MRKD-Tree方法，以針對(duì)每一組地震圖數(shù)據(jù)建立一棵或多棵KD樹(shù)。因此，創(chuàng)建了按照MRKD-Tree結(jié)構(gòu)進(jìn)行索引的歷史地震圖數(shù)據(jù)庫(kù)。根據(jù)本發(fā)明實(shí)施例，可以對(duì)具有這種MRKD-Tree結(jié)構(gòu)的數(shù)據(jù)庫(kù)應(yīng)用MRKD-Tree方法與近似最近鄰搜索方法的結(jié)合進(jìn)行搜索。根據(jù)本發(fā)明實(shí)施例，地震圖數(shù)據(jù)可以是多維度數(shù)據(jù)，這里“維度”可以是時(shí)間維度，或者也稱為時(shí)間采樣。對(duì)于每一組地震圖數(shù)據(jù)建立KD樹(shù)可以包括：計(jì)算組內(nèi)所有地震圖數(shù)據(jù)在每個(gè)維度上的均值和方差，確定多個(gè)維度中具有相對(duì)于均值的最大方差的一個(gè)或多個(gè)維度，并且對(duì)于所確定的一個(gè)或多個(gè)維度中的每一個(gè)維度，以該維度上的均值作為中位數(shù)，將地震圖數(shù)據(jù)劃分成兩部分，對(duì)每一個(gè)部分遞歸地進(jìn)行計(jì)算、確定和劃分步驟，直到每一部分只剩下單個(gè)地震圖數(shù)據(jù)，從而建立一棵或多棵KD樹(shù)。下面參照?qǐng)D5和6進(jìn)一步描述如何對(duì)地震圖數(shù)據(jù)應(yīng)用MRKD-Tree方法的示例。圖5顯示的是一個(gè)數(shù)據(jù)庫(kù)的10個(gè)地震圖在3個(gè)最大方差分布的時(shí)間維度以及相應(yīng)的均方差值(分別為1.67e-004，2.14e-004，1.85e-004)，其中地震圖數(shù)據(jù)按初至到時(shí)對(duì)齊。本發(fā)明不限于此，可以具有多個(gè)時(shí)間維度。此外，可以在每一個(gè)時(shí)間維度上計(jì)算均方差，選出均方差最大的幾個(gè)維度，作為地震圖數(shù)據(jù)的時(shí)間維度。圖6示出了建立MRKD-Tree結(jié)構(gòu)的過(guò)程的示意圖。如圖所示，用三個(gè)地震圖作數(shù)據(jù)庫(kù)，每個(gè)地震圖有四個(gè)維度。表格顯示了地震圖的幅值(地球介質(zhì)振動(dòng)速度)在每個(gè)維度的均值，以及其均方差。均方差的計(jì)算公式如下：其中，Ai(t)是數(shù)據(jù)庫(kù)里第i個(gè)地震圖在維度t的振幅幅值，α(t)是所有地震圖在維度t的振幅，N是數(shù)據(jù)庫(kù)里地震圖的個(gè)數(shù)。在圖6 的表格中，最大的均方差是在維度2和4，同為0.04667。那么可以同時(shí)在維度2和4開(kāi)始建樹(shù)。對(duì)于維度2，任何地震圖的幅值在維度2小于均值0.3的放在左邊，大于這個(gè)均值的放在右邊，從而進(jìn)行了劃分。然后繼續(xù)在左邊或右邊剩余地震圖中應(yīng)用該規(guī)則重復(fù)這個(gè)分樹(shù)過(guò)程，直至最后只剩一個(gè)地震圖。同時(shí)也在維度4開(kāi)始建另外一棵樹(shù)，依據(jù)同樣的劃分方法。由此，得到如圖6下部所示的MRKD-Tree結(jié)構(gòu)。綜上，在針對(duì)每一組地震圖數(shù)據(jù)建立MRKD-Tree結(jié)構(gòu)時(shí)，根據(jù)方差取值最大的維度(時(shí)間)來(lái)將所有地震圖數(shù)據(jù)分?jǐn)偝蓛刹糠?，即地震圖數(shù)據(jù)空間被一個(gè)正交超平面切分成了兩個(gè)空間。這兩部分的數(shù)據(jù)接著分別遞歸的進(jìn)行分?jǐn)傔^(guò)程，即空間切分過(guò)程，最后創(chuàng)建了平衡二叉樹(shù)。理論地震圖數(shù)據(jù)庫(kù)120實(shí)際的地震圖數(shù)據(jù)在地球內(nèi)部的覆蓋范圍是有限制的。一個(gè)新發(fā)生的地震可能或者不可能與以往歷史中的任何一個(gè)地震事件相似，特別是由于實(shí)際地震圖的數(shù)據(jù)庫(kù)不很完整，常有可能搜索不到非常相似的實(shí)際地震圖。因此，本發(fā)明提出了采用理論模擬地震圖數(shù)據(jù)作為補(bǔ)充，建立更大的地球內(nèi)部覆蓋范圍。于是，根據(jù)本發(fā)明實(shí)施例，預(yù)先建立理論地震圖數(shù)據(jù)庫(kù)120。這樣，對(duì)于新產(chǎn)生的地震事件，搜索過(guò)程是搜索歷史數(shù)據(jù)庫(kù)與理論數(shù)據(jù)庫(kù)的并行過(guò)程。這兩個(gè)并行搜索過(guò)程的返回結(jié)果都能夠提供重要的獨(dú)立信息，從而用于幫助估計(jì)當(dāng)前地震事件的地震參數(shù)。根據(jù)本發(fā)明實(shí)施例，由于實(shí)際地震圖數(shù)據(jù)是真實(shí)的地震記錄，沒(méi)有任何地球簡(jiǎn)化的假設(shè)，而理論地震圖數(shù)據(jù)都是在各種簡(jiǎn)化的假設(shè)條件下計(jì)算出來(lái)的，因此優(yōu)選地不將兩種數(shù)據(jù)庫(kù)混合在一起做搜索。當(dāng)從歷史地震圖數(shù)據(jù)庫(kù)和理論地震圖數(shù)據(jù)庫(kù)分別得到的匹配地震圖數(shù)據(jù)與輸入地震時(shí)間的地震圖數(shù)據(jù)的互相關(guān)系數(shù)相同或相近時(shí)，搜索引擎優(yōu)先選擇歷史地震圖數(shù)據(jù)庫(kù)的匹配地震圖數(shù)據(jù)，來(lái)估計(jì)地震參數(shù)。為了準(zhǔn)備一個(gè)充足的理論計(jì)算的地震圖數(shù)據(jù)庫(kù)，需要計(jì)算數(shù)百萬(wàn)的到數(shù)千萬(wàn)的人工模擬的地震圖。地震圖的特征，如橫波和縱波的時(shí) 間差，面波的頻譜特征，地震記錄的時(shí)間系列長(zhǎng)短等，隨監(jiān)測(cè)記錄距離的不同而不同。如果目標(biāo)是為了監(jiān)測(cè)一個(gè)局部的區(qū)域，可以采用三維有限元或者有限差分方法，用三維地球速度模型，網(wǎng)格化的震源參數(shù)(如震中距和震源深度等)來(lái)計(jì)算短波長(zhǎng)(高頻)的理論地震圖。方法如Zhang，W.，andX.Chen.2006.Tractionimagemethodforirregularfreesurfaceboundariesinfinitedifferenceseismicwavesimulation.GeophysicalJournalInternational167(1)：337-353.對(duì)于區(qū)域性的和超遠(yuǎn)距離的地震，可以采用彈性波模擬方法，使用全球平均的一維地球模型來(lái)計(jì)算長(zhǎng)波長(zhǎng)(低頻)的地震圖。一維模型如：Dziewonski，A.M.andAnderson，D.L.，(1981)，Preliminaryreferenceearthmodel.Phys.EarthPlanet.Inter.，25：297-356.無(wú)論距離遠(yuǎn)近，以上所有的地震圖的波形數(shù)據(jù)都會(huì)用最初信號(hào)到時(shí)或者最大振幅到時(shí)進(jìn)行對(duì)齊。此外，可以對(duì)特定的地震參數(shù)范圍進(jìn)行格點(diǎn)化，例如，對(duì)于震中距10°至35°(1°≈111公里)，深度0到50公里，按照距離間隔5公里，深度間隔5公里，走向角度間隔10°，傾角角度間隔10°，滑動(dòng)角度間隔10°，并對(duì)每一個(gè)格點(diǎn)，人工模擬或合成理論地震圖。同樣，可以采用SAC格式存儲(chǔ)理論地震圖。圖4示出了建立理論地震圖數(shù)據(jù)庫(kù)的一個(gè)示例的示意流程圖400，在該示例中，也應(yīng)用MRKD-Tree方法，對(duì)地震圖數(shù)據(jù)進(jìn)行索引。在步驟402，對(duì)特定的地震參數(shù)范圍，進(jìn)行格點(diǎn)化，并對(duì)每一個(gè)格點(diǎn)，計(jì)算理論地震圖。在步驟404，將特定的地震臺(tái)儀器頻譜響應(yīng)與理論地震圖做卷積，使得理論數(shù)據(jù)與實(shí)際數(shù)據(jù)可比，因?yàn)閷?shí)際接收的數(shù)據(jù)有地震臺(tái)頻譜響應(yīng)的卷積效應(yīng)。在步驟406，拾取P波初至到時(shí)。在步驟408，做分段頻率濾波并拾取濾波后的最大振幅到時(shí)。在步驟410，針對(duì)每個(gè)頻率分段，分別根據(jù)地震圖數(shù)據(jù)的初至到時(shí)和最大振幅到時(shí)來(lái)對(duì)齊波形，截取得到兩組地震圖數(shù)據(jù)。在步驟412，分別對(duì)得到的兩組地震圖數(shù)據(jù)應(yīng)用MRKD-Tree方法，以針對(duì)每一組地震圖數(shù)據(jù)建立一棵或多棵KD樹(shù)。因此，創(chuàng)建了按照MRKD-Tree結(jié)構(gòu)進(jìn)行索引的理論地震圖數(shù)據(jù)庫(kù)。在理論地震圖數(shù)據(jù)庫(kù)中，針對(duì)每一組地震圖數(shù) 據(jù)建立MRKD-Tree結(jié)構(gòu)的過(guò)程與在歷史地震圖數(shù)據(jù)庫(kù)中執(zhí)行的過(guò)程相同。輸入裝置130輸入裝置130可以是用于輸入地震事件的地震圖數(shù)據(jù)的接口，例如可以與單個(gè)或多個(gè)地震監(jiān)測(cè)臺(tái)的輸出裝置或通信接口相連，以從單個(gè)或多個(gè)地震監(jiān)測(cè)臺(tái)接收產(chǎn)生的地震事件的地震圖記錄。輸入裝置130也可以包括例如鍵盤、觸摸屏等用戶接口，用戶或操作者可以通過(guò)輸入裝置130輸入系統(tǒng)設(shè)置、預(yù)定條件等，來(lái)控制和管理系統(tǒng)的運(yùn)行?？梢圆捎萌魏芜m當(dāng)結(jié)構(gòu)或形式來(lái)實(shí)現(xiàn)輸入裝置130。搜索引擎140搜索引擎140可以是分立實(shí)現(xiàn)的設(shè)備，也可以實(shí)現(xiàn)為與歷史地震圖數(shù)據(jù)庫(kù)110、理論地震圖數(shù)據(jù)庫(kù)120和輸入裝置130中的一個(gè)或多個(gè)形成在同一位置。搜索引擎140可以利用微處理器、現(xiàn)場(chǎng)可編程門陣列或者專用計(jì)算設(shè)備來(lái)實(shí)現(xiàn)。一個(gè)充足的地震事件的數(shù)據(jù)庫(kù)包含數(shù)以億計(jì)到數(shù)以十億計(jì)的地震圖。即使是使用最先進(jìn)的計(jì)算機(jī)設(shè)備，如果采用原始的線性搜索方法，也沒(méi)有足夠的處理能力來(lái)查詢?nèi)绱舜罅康臄?shù)據(jù)。諸如近似最近鄰搜索(例如，參見(jiàn)Muja，M.，Lowe，D.G.(2009).Fastapproximatenearestneighborswithautomaticalgorithmconfiguration，VISAPP(1)2009：331-340)的計(jì)算機(jī)搜索技術(shù)使得我們可以從大規(guī)模數(shù)據(jù)庫(kù)中在一定準(zhǔn)確率程度下進(jìn)行快速的查找，從而得到相似的地震圖，并且這一技術(shù)在互聯(lián)網(wǎng)的數(shù)據(jù)量級(jí)下，已經(jīng)成功地應(yīng)用于圖片檢索，音樂(lè)檢索以及視頻檢索領(lǐng)域中。搜索引擎140采用近似最近鄰方法在歷史地震圖數(shù)據(jù)庫(kù)110和理論地震圖數(shù)據(jù)庫(kù)120中進(jìn)行并行搜索，按照預(yù)定的相似度條件找出與該輸入地震事件的地震圖數(shù)據(jù)相似的多個(gè)地震圖數(shù)據(jù)，作為搜索結(jié)果集合。例如，預(yù)定的相似度條件可以是預(yù)定的搜索準(zhǔn)確率，例如90％，凡是滿足這一條件的地震圖數(shù)據(jù)都可以作為搜索結(jié)果。根據(jù)本發(fā)明實(shí)施例，搜索引擎140在搜索過(guò)程中還進(jìn)一步結(jié)合了局部敏感哈希(LSH)方法或者M(jìn)RKD-Tree方法，從而進(jìn)一步加快了搜索過(guò)程。圖 7示出了根據(jù)本發(fā)明實(shí)施例的搜索引擎140的示意框圖。搜索引擎140包括搜索裝置710和確定裝置720。針對(duì)通過(guò)輸入裝置130輸入的地震事件的包括地震波形的地震圖數(shù)據(jù)，搜索裝置710在歷史和理論地震圖數(shù)據(jù)庫(kù)中通過(guò)波形匹配，采用近似最近鄰方法進(jìn)行并行搜索，按照預(yù)定的相似度條件找出與該輸入地震事件的地震圖數(shù)據(jù)相似的多個(gè)地震圖數(shù)據(jù)，作為搜索結(jié)果集合。確定裝置720在搜索結(jié)果集合中，確定匹配地震圖數(shù)據(jù)，并從中獲取輸入地震事件的地震參數(shù)估計(jì)。根據(jù)本發(fā)明實(shí)施例，系統(tǒng)10還可以包括輸出裝置(未示出)，連接至搜索引擎140，輸出搜索引擎140獲取的地震參數(shù)估計(jì)。輸出裝置例如可以是顯示屏、揚(yáng)聲器等。根據(jù)本發(fā)明實(shí)施例，搜索裝置710可以在歷史和理論地震圖數(shù)據(jù)庫(kù)中，針對(duì)多組地震圖數(shù)據(jù)，并行地進(jìn)行近似最近鄰搜索，按照預(yù)定的相似度條件(例如，預(yù)定的搜索準(zhǔn)確率，例如90％)，得到每組獨(dú)立的多個(gè)搜索結(jié)果，并按照相似度從高到低得出一個(gè)整體的優(yōu)先隊(duì)列池，取出前面若干個(gè)(例如，20個(gè))最近鄰，作為最終搜索結(jié)果集合。這些搜索結(jié)果可以看成是該大規(guī)模數(shù)據(jù)庫(kù)所有數(shù)據(jù)的一個(gè)子集。這里，相似度可以基于輸入地震圖的地震圖數(shù)據(jù)與搜索結(jié)果集合中每個(gè)搜索結(jié)果之間的歐式距離，即，兩者之間的差值。因?yàn)樗阉鞯降牡卣饒D在預(yù)處理對(duì)齊時(shí)可能存在時(shí)間選擇上的偏差，所以搜索結(jié)果根據(jù)歐式距離大小順序進(jìn)行排列就可能不是很準(zhǔn)確。根據(jù)本發(fā)明實(shí)施例，相似度可以基于輸入地震圖的地震圖數(shù)據(jù)與搜索結(jié)果集合中每個(gè)搜索結(jié)果之間的互相關(guān)系數(shù)。(歸一化的)互相關(guān)系數(shù)可以如下計(jì)算：其中，x(i)和y(i)是兩個(gè)時(shí)間序列地震圖，每個(gè)都有N個(gè)分立值(i＝1，2，3，...N)，mx是x(i)的平均值，my是y(i)的平均值。相關(guān)產(chǎn)生的時(shí)間序列r(j)(j＝1，2，...N)值應(yīng)該在[-1，+1]之間。這個(gè)值越大，表明兩個(gè)地震圖越相似。雖然互相關(guān)系數(shù)的計(jì)算開(kāi)銷遠(yuǎn)大于歐式距離的 計(jì)算開(kāi)銷，但是從相似度度量來(lái)看，其準(zhǔn)確程度要?dú)W式距離更為準(zhǔn)確。由于只在近似最近鄰的搜索結(jié)果里計(jì)算互相關(guān)系數(shù)，其數(shù)據(jù)規(guī)模要遠(yuǎn)小于整個(gè)數(shù)據(jù)庫(kù)，因此計(jì)算互相關(guān)系數(shù)不會(huì)花費(fèi)很多的時(shí)間。根據(jù)本發(fā)明實(shí)施例，確定裝置720計(jì)算輸入地震圖的地震圖數(shù)據(jù)與搜索結(jié)果集合中每個(gè)搜索結(jié)果之間的互相關(guān)系數(shù)，并采用互相關(guān)系數(shù)最大的搜索結(jié)果作為匹配地震圖數(shù)據(jù)，從而確定匹配地震圖數(shù)據(jù)。如果搜索裝置710已按照互相關(guān)系數(shù)的大小形成了整體的優(yōu)先隊(duì)列池，則確定裝置720可以直接從優(yōu)先隊(duì)列池中選擇互相關(guān)系數(shù)最大的搜索結(jié)果，作為匹配地震圖數(shù)據(jù)。當(dāng)從歷史地震圖數(shù)據(jù)庫(kù)和理論地震圖數(shù)據(jù)庫(kù)分別得到的匹配地震圖數(shù)據(jù)與輸入地震時(shí)間的地震圖數(shù)據(jù)的互相關(guān)系數(shù)相同或相近(例如，相似度滿足預(yù)定閾值)時(shí)，確定裝置720選擇歷史地震圖數(shù)據(jù)庫(kù)的匹配地震圖數(shù)據(jù)。圖9示出了根據(jù)本發(fā)明實(shí)施例的獲取地震事件的地震參數(shù)估計(jì)的方法的流程圖。方法800包括在步驟802，輸入地震事件的地震圖數(shù)據(jù)。地震圖數(shù)據(jù)可以是在地震監(jiān)測(cè)臺(tái)記錄的，并且是經(jīng)過(guò)預(yù)處理的。預(yù)處理可以包括：進(jìn)行信號(hào)處理、去噪，拾取P波初至?xí)r間；對(duì)地震圖數(shù)據(jù)進(jìn)行分段頻率濾波以分成多個(gè)頻率分段，并拾取濾波后的面波的最大振幅到時(shí)；針對(duì)每個(gè)頻率分段，分別根據(jù)地震圖數(shù)據(jù)的初至到時(shí)和最大振幅到時(shí)來(lái)對(duì)齊波形，截取得到兩組地震圖數(shù)據(jù)。方法800還包括：在步驟804，在歷史和理論地震圖數(shù)據(jù)庫(kù)中通過(guò)波形匹配，采用近似最近鄰搜索方法進(jìn)行并行搜索，按照預(yù)定的相似度條件找出與該輸入地震事件的地震圖數(shù)據(jù)相似的多個(gè)地震圖數(shù)據(jù)，作為搜索結(jié)果集合。這里結(jié)合圖6所示示例進(jìn)行具體描述。如上所述，圖6示意性示出了針對(duì)數(shù)據(jù)庫(kù)中的每一組地震圖數(shù)據(jù)建立MRKD-Tree結(jié)構(gòu)。該圖顯示的是同時(shí)建立兩棵樹(shù)的例子。具體而言，在輸入數(shù)據(jù)庫(kù)地震圖建樹(shù)時(shí)，并行地分別判斷每個(gè)地震圖數(shù)據(jù)在其維度2和4的值是否大于或小于樹(shù)頂?shù)木?.3和0.6。例如對(duì)于建立從維度2開(kāi)始的樹(shù)，任何地震圖的幅值在維度2小于均值0.3的放在左邊(ID＝1，2)，大于這個(gè)均值的放在右邊(ID＝3)。放在左邊的所有地震圖(ID＝1，2)的最大均值是維度4的0.65。同樣，在維度4，均值 小于0.65的地震圖放在節(jié)點(diǎn)的左邊(ID＝1)，均值大于0.65的地震圖放在節(jié)點(diǎn)的右邊(ID＝2)，直到到達(dá)了底層的最后一個(gè)地震圖，這樣便建立了一棵MRKD-Tree，至此樹(shù)1建樹(shù)結(jié)束。用同樣的方式我們建立了右側(cè)(樹(shù)2)從維度4開(kāi)始的樹(shù)。圖7和圖8顯示兩輪搜索的過(guò)程。圖7顯示的是第一輪搜索過(guò)程，即在每棵MRKD-Tree中從根結(jié)點(diǎn)開(kāi)始查詢，具體在每棵MRKD-Tree中是按照如下方式進(jìn)行查詢的：首先，從根結(jié)點(diǎn)開(kāi)始，比較根結(jié)點(diǎn)的“均值”與待查詢地震圖數(shù)據(jù)的“選擇維度”的值，如果待查詢地震圖數(shù)據(jù)的“選擇維度”的值小于等于根結(jié)點(diǎn)的“均值”，則進(jìn)入根結(jié)點(diǎn)的左兒子結(jié)點(diǎn)，同時(shí)計(jì)算右兒子與待查詢地震圖之間的啟發(fā)式距離信息(hdist)并將右兒子結(jié)點(diǎn)插入到優(yōu)先隊(duì)列池中，反之，如果大于根結(jié)點(diǎn)的“均值”，則進(jìn)入根結(jié)點(diǎn)的右兒子結(jié)點(diǎn)，同時(shí)計(jì)算左兒子相應(yīng)的啟發(fā)式距離信息(hdist)并將左兒子結(jié)點(diǎn)插入到優(yōu)先隊(duì)列池中，如圖7中的a)所示，進(jìn)入樹(shù)1的“選擇維度＝4，均值＝0.65”的結(jié)點(diǎn)，并且計(jì)算樹(shù)1的ID＝3結(jié)點(diǎn)的啟發(fā)式距離信息并該結(jié)點(diǎn)插入到優(yōu)先隊(duì)列池，按照同樣的處理方式，遞歸的比較內(nèi)部結(jié)點(diǎn)的“均值”與待查詢地震圖在該結(jié)點(diǎn)對(duì)應(yīng)的“選擇維度”的值，進(jìn)入左兒子或者右兒子結(jié)點(diǎn)，并且計(jì)算“未進(jìn)入”的兒子結(jié)點(diǎn)的啟發(fā)式距離信息并將其插入優(yōu)先隊(duì)列池中，直到到達(dá)葉子結(jié)點(diǎn)為止，此時(shí)將該葉子結(jié)點(diǎn)對(duì)應(yīng)ID加入到結(jié)果隊(duì)列中去，如圖7中的b)所示；如果結(jié)果隊(duì)列中已經(jīng)包含了該葉子結(jié)點(diǎn)對(duì)應(yīng)的ID，那么將跳過(guò)此次操作，不會(huì)重復(fù)插入相同的ID到結(jié)果隊(duì)列中；至此完成了一棵樹(shù)的第一輪搜索，按照同樣的處理方式，搜索每棵MRKD-Tree，如圖7中的c)所示，演示了樹(shù)2的第一輪搜索結(jié)果。搜索完每棵MRKD-Tree標(biāo)志著第一輪搜索的結(jié)束。圖8顯示的是從第二輪開(kāi)始的搜索流程。首先，從優(yōu)先隊(duì)列池中取出隊(duì)列首部的結(jié)點(diǎn)，即啟發(fā)式距離信息(hdist)最小的結(jié)點(diǎn)；然后，根據(jù)此結(jié)點(diǎn)的信息，查找到對(duì)應(yīng)的樹(shù)標(biāo)號(hào)，以及結(jié)點(diǎn)位置信息，從而，由此結(jié)點(diǎn)開(kāi)始，按照第一輪的自上到下的方式，比較該結(jié)點(diǎn)的“均值”與待查詢地震圖在該結(jié)點(diǎn)對(duì)應(yīng)“選擇維度”的值，如果待查詢地震圖 的“選擇維度”的值小于等于該結(jié)點(diǎn)“均值”，則進(jìn)入左兒子結(jié)點(diǎn)，并且計(jì)算其右兒子結(jié)點(diǎn)的啟發(fā)式距離信息，將其加入優(yōu)先隊(duì)列池，若是大于該結(jié)點(diǎn)“均值”，則左右兒子互換，按照同樣方式進(jìn)行；如此反復(fù)，按照這種方式遞歸的進(jìn)入兒子結(jié)點(diǎn)，直到到達(dá)葉子結(jié)點(diǎn)為止，然后將葉子結(jié)點(diǎn)加入結(jié)果隊(duì)列，若隊(duì)列中已經(jīng)含有該結(jié)點(diǎn)，則不重復(fù)加入。這些操作都與第一輪中的相同。需要注意的是，如果我們指定了結(jié)果隊(duì)列大小的上限，那么結(jié)果隊(duì)列中結(jié)點(diǎn)個(gè)數(shù)一旦達(dá)到這個(gè)上限，則整個(gè)搜索過(guò)程結(jié)束，否則繼續(xù)按照如上方式進(jìn)行查詢，直到所有結(jié)點(diǎn)都加入結(jié)果隊(duì)列為止。圖8中的a)到d)每個(gè)圖對(duì)應(yīng)的是一次從優(yōu)先隊(duì)列池中取出結(jié)點(diǎn)并查詢的過(guò)程。對(duì)于這種多棵KD樹(shù)搜索，原始的基礎(chǔ)的計(jì)算機(jī)理論與方法在Silpa-Anan，C.，andHartley，R.(2008).OptimisedKD-treesforfastimagedescriptormatching，InCVPR論文中有詳盡描述。方法800還包括：在步驟806，從步驟804中找出的搜索結(jié)果集合中確定匹配地震圖數(shù)據(jù)，并從中獲取輸入地震事件的地震參數(shù)估計(jì)。地震參數(shù)可以包括如下中的至少一個(gè)：震中距、震源深度、震源機(jī)制和震級(jí)。根據(jù)本發(fā)明實(shí)施例，方法800還可以包括：通過(guò)比較所有匹配地震圖數(shù)據(jù)的互相關(guān)系數(shù)分布，來(lái)評(píng)估搜索結(jié)果的解析度和可靠性。下面描述利用本發(fā)明實(shí)施例的模擬示例，其中估計(jì)一個(gè)人工合成地震的地震參數(shù)。假定輸入人工合成地震圖的震中距300公里，震源深度25公里，震源機(jī)制：走向角度＝170°，傾角角度＝230°，滑動(dòng)角度＝170°。對(duì)于這個(gè)示例，理論數(shù)據(jù)庫(kù)可以采用如下參數(shù)來(lái)進(jìn)行格點(diǎn)化：震中距100至1000km，以5km為間隔，震源深度5-30公里，以5km為間隔，震源機(jī)制：走向角度＝0至360°，傾角角度＝0至90°，滑動(dòng)角度＝-180°至180°，間隔為10°。將人工合成地震圖的波形輸入根據(jù)本發(fā)明實(shí)施例的系統(tǒng)10，圖11顯示了一個(gè)搜索結(jié)果的示例，包括MRKD-Tree方法搜索的二十個(gè)最佳匹配結(jié)果。最上方的粗線的地震圖是搜索出的與輸入地震圖完全相同的地震圖，以及其它19個(gè)人工模擬的相似地震事件的地震圖。圖12示出了三種方法的地震圖搜索時(shí)間的比較：精確的線性搜索(LinearSearch)，局部敏感哈希(LSH)方法，以及MRKD-Tree方法。在搜索相同數(shù)量的地震圖的數(shù)據(jù)庫(kù)時(shí)，線性精確搜索返回100個(gè)最佳匹配結(jié)果，并且得到了100％的準(zhǔn)確率和最長(zhǎng)的搜索時(shí)間。LSH方法和MRKD-Tree方法得到的結(jié)果有90％的準(zhǔn)確率，但是MRKD-Tree方法所用時(shí)間遠(yuǎn)遠(yuǎn)小于LSH方法和線性精確搜索方法。對(duì)于實(shí)際應(yīng)用來(lái)說(shuō)，90％的準(zhǔn)確率已經(jīng)足夠?？梢?jiàn)，MRKD-Tree方法的速度和效率非常高，根據(jù)實(shí)驗(yàn)結(jié)果，在數(shù)以百萬(wàn)計(jì)的數(shù)據(jù)量的數(shù)據(jù)庫(kù)上，MRKD-Tree方法僅僅花費(fèi)數(shù)秒時(shí)間就可以得到搜索結(jié)果，實(shí)現(xiàn)了快速甚至實(shí)時(shí)的地震參數(shù)估計(jì)。圖13示出了輸入地震事件和使用線性搜索(LinearSearch)和MRKD-Tree方法分別得到的最佳匹配地震事件之間的最大互相關(guān)系數(shù)的比較。從圖中可以看出，盡管MRKD-Tree搜索得到的結(jié)果中缺失了50個(gè)匹配地震圖數(shù)據(jù)中的若干個(gè)，但是對(duì)于實(shí)際應(yīng)用，MRKD-Tree方法已經(jīng)能夠搜索得到了足夠多統(tǒng)計(jì)上需要的匹配地震圖數(shù)據(jù)?？紤]到MRKD-Tree方法相比于線性精確搜索的非常高的速度和效率，根據(jù)本發(fā)明實(shí)施例的應(yīng)用MRKD-Tree方法的方法和系統(tǒng)具有幾乎實(shí)時(shí)獲取地震參數(shù)估計(jì)的優(yōu)勢(shì)。以上描述了根據(jù)本發(fā)明實(shí)施例的快速獲取輸入地震事件的地震參數(shù)估計(jì)的方法和系統(tǒng)以及地震時(shí)間搜索引擎，能夠在得到輸入地震數(shù)據(jù)后數(shù)秒鐘內(nèi)估計(jì)出該輸入地震事件的地震參數(shù)。該方法和系統(tǒng)使用地震事件搜索引擎，通過(guò)建立大規(guī)模的歷史實(shí)際地震資料數(shù)據(jù)庫(kù)和建立能夠完整覆蓋所有震源參數(shù)的理論地震圖的大規(guī)模數(shù)據(jù)庫(kù)，快速?gòu)臄?shù)據(jù)庫(kù)中搜索出相似于輸入地震圖的所有地震事件。一旦搜索得到了輸入地震事件對(duì)應(yīng)的匹配結(jié)果，該系統(tǒng)立即顯示出這些匹配地震事件的地震參數(shù)，以便獲取例如震源機(jī)制和其他震源參數(shù)的估計(jì)。本發(fā)明實(shí)施例方法和系統(tǒng)能夠利用計(jì)算機(jī)搜索技術(shù)，應(yīng)用于天然地震的實(shí)時(shí)監(jiān)測(cè)，實(shí)時(shí)監(jiān)測(cè)石油與天然氣開(kāi)采生產(chǎn)所激發(fā)的微震事件，并且也可以用于實(shí)時(shí)監(jiān)測(cè)核爆或者炸藥爆炸等。以上所述是本發(fā)明的示例實(shí)施方式，應(yīng)當(dāng)指出，對(duì)于本技術(shù)領(lǐng)域的普通技術(shù)人員來(lái)說(shuō)，在不脫離本發(fā)明所述原理的前提下，還可以作出若干改進(jìn)和潤(rùn)飾，這些改進(jìn)和潤(rùn)飾也應(yīng)視為本發(fā)明的保護(hù)范圍。