本發(fā)明實施例涉及石油天然氣地震勘探技術(shù)，尤其涉及一種微地震事件的檢測和定位方法與裝置。

背景技術(shù)：

石油天然氣地震勘探過程中，需要比較快速地得到壓裂結(jié)果以對壓裂效果進(jìn)行評價，進(jìn)而對開采方案進(jìn)行及時調(diào)整，因此，需要對水力壓裂微地震監(jiān)測資料進(jìn)行實時的處理解釋，準(zhǔn)確檢測和定位出微地震事件的位置。

現(xiàn)有井中微地震事件的檢測是基于STA/LTA(Short Term Average/Long Term Average，長短時窗能量比法)的方法，其定位主要是通過拾取有效事件的初至波，利用初至?xí)r間進(jìn)行反演來確定微地震事件的位置，對于單井監(jiān)測系統(tǒng)，還要利用P波或S波方位角的信息反演確定微地震事件的位置。

然而，現(xiàn)有技術(shù)其檢測微地震事件的數(shù)量和定位微地震事件的精度強(qiáng)烈依賴于微地震事件的信噪比，對于低信噪比的微地震事件，現(xiàn)有技術(shù)不能準(zhǔn)確檢測到，更不能準(zhǔn)確的拾取事件的初至波，進(jìn)而無法實現(xiàn)對低信噪比的微地震事件的準(zhǔn)確定位。

技術(shù)實現(xiàn)要素：

本發(fā)明實施例提供一種微地震事件的檢測和定位方法與裝置，用以解決現(xiàn)有技術(shù)無法準(zhǔn)確檢測和定位低信噪比的微地震事件的位置的技術(shù)問題。

第一方面，本發(fā)明實施例提供一種微地震事件的檢測和定位方法，包括：

獲取模板事件的網(wǎng)格結(jié)構(gòu)中的每個網(wǎng)格與不同觀測站之間的P波校正時差、所述網(wǎng)格結(jié)構(gòu)中的每個網(wǎng)格與不同觀測站之間的S波校正時差、所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)和所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)；其中，所述網(wǎng)格結(jié)構(gòu)包括多個網(wǎng)格，所述模板事件為觀測站的原始數(shù)據(jù)中信噪比大于預(yù)設(shè)閾值的事件；

根據(jù)所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)、所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)、每個所述網(wǎng)格與不同觀測站之間的P波校正時差和每個所述網(wǎng)格與不同觀測站之間的S波校正時差，確定每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)，所述第一疊加互相系數(shù)用于表征所述網(wǎng)格為震源的概率；

根據(jù)每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)，確定所述模板事件對應(yīng)的微地震事件的位置和時刻。

第二方面，本發(fā)明實施例提供一種微地震事件的檢測和定位裝置，包括：

獲取模塊，用于獲取模板事件的網(wǎng)格結(jié)構(gòu)、所述模板事件的P波和S波、所述網(wǎng)格結(jié)構(gòu)中的每個網(wǎng)格與不同觀測站之間的P波校正時差、所述網(wǎng)格結(jié)構(gòu)中的每個網(wǎng)格與不同觀測站之間的S波校正時差、所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)和所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)；其中，所述網(wǎng)格結(jié)構(gòu)包括多個網(wǎng)格，所述模板事件為原始數(shù)據(jù)中信噪比大于預(yù)設(shè)閾值的事件；

處理模塊，用于根據(jù)所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)、所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)、每個所述網(wǎng)格與不同觀測站之間的P波校正時差和每個所述網(wǎng)格與不同觀測站之間的S波校正時差，確定每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)，所述第一疊加互相系數(shù)用于表征所述網(wǎng)格的地震強(qiáng)度；

確定模塊，用于根據(jù)每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)和第一預(yù)設(shè)閾值，確定所述模板事件對應(yīng)的微地震事件的位置和時刻。。

本發(fā)明實施例提高的微地震事件的檢測和定位方法與裝置，通過使用每個網(wǎng)格的與各觀測站之間P波校正時差對模板事件的各P波互相關(guān)系數(shù)進(jìn)行校正，使用每個網(wǎng)格與各觀測站之間的S波校正時差對模板事件的各S波互相關(guān)系數(shù)進(jìn)行校正，并將校正后的P波互相關(guān)系數(shù)和S波互相關(guān)系數(shù)進(jìn)行疊加，獲得每個網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)，并獲得滿足第一預(yù)設(shè)閾值的第一疊加互相關(guān)系數(shù)所對應(yīng)的網(wǎng)格的位置和時刻，進(jìn)而實現(xiàn)對低信噪比的微地震事件的準(zhǔn)確檢測和定位。

附圖說明

為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作一簡單地介紹，顯而易見地，下面描述中的附圖是本發(fā)明的一些實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。

圖1為本發(fā)明提供的微地震事件的檢測和定位方法實施例一的流程示意圖；

圖1a為模板事件與原始數(shù)據(jù)在E分量上的P波和S波互相關(guān)系數(shù)；

圖1b為模板事件與原始數(shù)據(jù)在N分量上的P波和S波互相關(guān)系數(shù)；

圖1c為模板事件與原始數(shù)據(jù)在Z分量上的P波和S波互相關(guān)系數(shù)；

圖2為本發(fā)明提供的微地震事件的檢測和定位方法實施例二的流程示意圖；

圖3為本發(fā)明提供的微地震事件的檢測和定位方法實施例三的流程示意圖；

圖4為本發(fā)明提供的微地震事件的檢測和定位方法實施例四的流程示意圖；

圖5為本發(fā)明提供的微地震事件的檢測和定位裝置實施例一的結(jié)構(gòu)示意圖；

圖6為本發(fā)明提供的微地震事件的檢測和定位裝置實施例二的結(jié)構(gòu)示意圖；

圖7為本發(fā)明提供的微地震事件的檢測和定位裝置實施例三的結(jié)構(gòu)示意圖。

具體實施方式

為使本發(fā)明實施例的目的、技術(shù)方案和優(yōu)點更加清楚，下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實施例是本發(fā)明一部分實施例，而不是全部的實施例?；诒景l(fā)明中的實施例，本領(lǐng)域普通技術(shù)人員在沒有作出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都屬于本發(fā)明保護(hù)的范圍。

本發(fā)明實施例提供的微地震事件的檢測和定位方法與裝置，適用于石油天然氣地震勘探中對微地震震源位置的檢測和定位，用于解決現(xiàn)有技術(shù)無法準(zhǔn)確檢測和定位低信噪比的微地震事件的位置的技術(shù)問題。

本實施例的技術(shù)方案，通過原始模板數(shù)據(jù)獲取模板事件，并獲得模板事件的網(wǎng)格結(jié)構(gòu)，使用每個網(wǎng)格的P波校正時差對模板事件的P波互相關(guān)系數(shù)進(jìn)行校正，使用每個網(wǎng)格的S波校正時差對模板事件的S波互相關(guān)系數(shù)進(jìn)行校正，并將校正后的P波互相關(guān)系數(shù)和S波互相關(guān)系數(shù)進(jìn)行疊加，獲得每個網(wǎng)格的第一疊加互相關(guān)系數(shù)，并跟進(jìn)每個網(wǎng)格的第一疊加互相系數(shù)和第一預(yù)設(shè)閾值，準(zhǔn)確獲得模板事件對應(yīng)的微地震事件的位置，進(jìn)而實現(xiàn)對低信噪比的微地震事件的準(zhǔn)確檢測和定位。

需要說明的是，本文中的術(shù)語“第一”、“第二”僅用于描述目的，而不能理解為指示或暗示相對重要性。

下面以具體地實施例對本發(fā)明的技術(shù)方案進(jìn)行詳細(xì)說明。下面這幾個具體的實施例可以相互結(jié)合，對于相同或相似的概念或過程可能在某些實施例不再贅述。

圖1為本發(fā)明提供的微地震事件的檢測和定位方法實施例一的流程示意圖。本實施例的方法具體可以應(yīng)用于微地震的檢測和定位過程中。本實施例提供的微地震的檢測和定位方法可以通過具有檢測和定位微地震事件的裝置(以下簡稱檢測和定位裝置)來執(zhí)行，該檢測和定位裝置可以集成在具有微地震定位功能的設(shè)備中，例如微地震檢測設(shè)備中，也可以是單獨設(shè)備。本實施例涉及的是檢測和定位裝置根據(jù)原始數(shù)據(jù)對模板事件對應(yīng)的微地震的位置進(jìn)行定位的具體過程。如圖1所示，本實施例可以包括：

S101、獲取模板事件的網(wǎng)格結(jié)構(gòu)中的每個網(wǎng)格與不同觀測站之間的P波校正時差、所述網(wǎng)格結(jié)構(gòu)中的每個網(wǎng)格與不同觀測站之間的S波校正時差、所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)和所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)。

其中，所述網(wǎng)格結(jié)構(gòu)包括多個網(wǎng)格，所述模板事件為原始數(shù)據(jù)中信噪比大于預(yù)設(shè)閾值的事件。

在實際情況中，高信噪比的微地震事件周圍會發(fā)生低信噪比的微地震事件，本實施例從各觀測站觀測的原始數(shù)據(jù)中選出高信噪比的模板事件，以該高信噪比的模板事件為基礎(chǔ)獲得低信噪比的微地震事件的發(fā)生位置。

具體的，本實施例的技術(shù)方案，從各觀測站的檢波器中獲得原始數(shù)據(jù)，從原始數(shù)據(jù)中獲取信噪比大于預(yù)設(shè)閾值的事件，將這些事件記為模板事件?？蛇x的，本實施例的檢測和定位裝置可以使用STA/LTA方法掃描原始數(shù)據(jù)，將檢測到的信噪比大于預(yù)設(shè)閾值的事件作為模板事件。其中預(yù)設(shè)閾值的具體設(shè)定根據(jù)實際需要進(jìn)行確定，只要保證根據(jù)該預(yù)設(shè)閾值可以從原始數(shù)據(jù)中獲取至少一個模板事件即可，本實施例對模板事件的數(shù)目不做限制。

需要說明的是，為了方便闡述本實施例的技術(shù)方案，本實施例以一個模板事件為例進(jìn)行說明，當(dāng)包括多個模板事件時，其定位過程只是一個模板事件的定位過程的重復(fù)，具體參考一個模板事件的定位過程即可。

進(jìn)一步的，上述檢測和定位裝置獲取模板事件的同時，還獲取該模板事件的P波和S波，該P波的到時和走時以及該S波的到時和走時，并根據(jù)P波的到時和走時以及S波的到時和走時，獲得模板事件的位置。

接著，檢測和定位裝置獲取上述模板事件的網(wǎng)格結(jié)構(gòu)，例如可以以該模板事件的位置為中心劃分網(wǎng)格，獲得該模板事件的網(wǎng)格結(jié)構(gòu)，可選的，該網(wǎng)絡(luò)結(jié)構(gòu)還可以是任意包括該模板事件的位置的任意網(wǎng)格結(jié)構(gòu)。其中，網(wǎng)格結(jié)構(gòu)的大小以及網(wǎng)格結(jié)構(gòu)中各網(wǎng)格的大小根據(jù)實際需要設(shè)定，其中網(wǎng)格越小則計算結(jié)果越準(zhǔn)確，網(wǎng)格結(jié)構(gòu)中各網(wǎng)格的大小可以相同也可以不同，例如越靠近模板事件的位置網(wǎng)格的尺寸越小。

其中，各網(wǎng)格的P波和S波分別與模板事件的P波和S波的波形相同。

檢測和定位裝置獲得模板事件的網(wǎng)格結(jié)構(gòu)之后，計算每個網(wǎng)格與各觀測站之間的P波校正時差和S波校正時差。具體的，檢測和定位裝置根據(jù)模板事件的P波相對于各觀測站的走時和到時，以及網(wǎng)格X的P波相對于各觀測站的走時，獲得該網(wǎng)格X與各觀測站之間的P波校正時差。同理，根據(jù)模板事件的S波相對于各觀測站的走時和到時，以及網(wǎng)格X的S波相對于各觀測站的走時，獲得該網(wǎng)格X與各觀測站之間的S波校正時差。

檢測和定位裝置根據(jù)模板事件的P波與每個觀測站的原始數(shù)據(jù)，獲得模板事件的P波與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)。如圖1a-圖1c所示，其中，ZA01-ZA20表示不同觀測站，波形數(shù)據(jù)為各觀測站中的檢波器在同一時間段(0.00-0.30)內(nèi)檢測到的原始數(shù)據(jù)，PCC表示模板事件的P波與各觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)，SCC表示模板事件的S波與各觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)。如圖1a-圖1c所示，由于本實施例采用的檢波器為三分量檢波器，因此檢波器獲得的原始數(shù)據(jù)包括三分量(即分量E、分量N和分量Z)上的波形數(shù)據(jù)，對應(yīng)的檢測和定位裝置獲得模板事件的P波與原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)也包括三個分量，將這三個分量上的互相關(guān)系數(shù)進(jìn)行相加，獲得模板事件在一個檢波器上的P波互相關(guān)系數(shù)。同理檢測和定位裝置將模板事件的S波與原始數(shù)據(jù)在三個分量上的互相關(guān)系數(shù)進(jìn)行疊加，獲得模板事件在一個檢波器上的S波互相關(guān)系數(shù)。

S102、根據(jù)所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)、所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)、每個所述網(wǎng)格與不同觀測站之間的P波校正時差和每個所述網(wǎng)格與不同觀測站之間的S波校正時差，確定每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)。

其中，該所述第一疊加互相系數(shù)用于表征所述網(wǎng)格為震源的概率，第一疊加互相系數(shù)越大，該第一疊加互相關(guān)系數(shù)對應(yīng)的網(wǎng)格對應(yīng)的位置越可能是地震源的位置，因此，可以通過第一疊加互相關(guān)系數(shù)的大小來確定震源的位置。

具體的，檢測和定位裝置獲取模板事件的各P波互相關(guān)系數(shù)和各S波互相關(guān)系數(shù)之后，使用每個網(wǎng)格的P波校正時差對模板事件的各P波互相關(guān)系數(shù)進(jìn)行校正，使用每個網(wǎng)格的S波校正時差對模板事件的各S波互相關(guān)系數(shù)進(jìn)行校正。接著，將校正后的各P波互相關(guān)系數(shù)和各S波互相關(guān)系數(shù)進(jìn)行平均或者加權(quán)平均，獲得每個網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)。

需要說明的是，本實施例的方法均是針對各觀測站進(jìn)行的，需要將每一個觀測站對應(yīng)的數(shù)據(jù)進(jìn)行疊加，獲得最終的結(jié)果。例如，使用每個網(wǎng)格與觀測站A的P波校正時差去校正模板事件與觀測站A在時刻t的P波互相關(guān)系數(shù)，使用每個網(wǎng)格與觀測站A的S波校正時差去校正模板事件與觀測站A在時刻t的S波互相關(guān)系數(shù)，接著，將每個觀測站對應(yīng)的校正后的P波互相關(guān)系數(shù)和校正后的S波互相關(guān)系數(shù)進(jìn)行疊加求平均值，獲得每個網(wǎng)格在不同在時刻t的第一疊加互相關(guān)系數(shù)。

需要說明的是，對于同一微地震事件各檢波器觀測到的地震波基本相同，但由于檢波器的位置不同使得各檢波器觀測到地震波的到時不同，因此可以根據(jù)各網(wǎng)格的校正時差對模板事件的互相關(guān)系數(shù)進(jìn)行校正，并將校正后的互相關(guān)系數(shù)進(jìn)行疊加，獲得每個網(wǎng)格的第一疊加互相關(guān)系數(shù)，每個網(wǎng)格的第一疊加互相關(guān)系數(shù)的大小可以反映每個網(wǎng)格為震源的概率大小，例如網(wǎng)格的第一疊加互相關(guān)系數(shù)越大，則該網(wǎng)格的位置時震源位置的概率也越大。

S103、根據(jù)每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)和第一預(yù)設(shè)閾值，獲得所述模板事件對應(yīng)的微地震事件的位置和時刻。

具體的，根據(jù)上述步驟獲得每個網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)，將各第一疊加系數(shù)與第一預(yù)設(shè)閾值進(jìn)行比較，獲得滿足第一預(yù)設(shè)閾值的第一疊加互相關(guān)系數(shù)。獲取這些滿足第一預(yù)設(shè)閾值的第一疊加互相關(guān)系數(shù)對應(yīng)的時間和網(wǎng)格的位置，將這些網(wǎng)格的位置作為該模板事件對應(yīng)的微地震事件的位置，對應(yīng)的將上述時間作為該微地震事件的發(fā)震時間，進(jìn)而實現(xiàn)對微地震事件的準(zhǔn)確檢測和定位，從而提高了油氣開采的可靠性。

本發(fā)明實施例提高的微地震事件的檢測和定位方法，通過從原始數(shù)據(jù)中獲取模板事件，并獲得模板事件的網(wǎng)格結(jié)構(gòu)、模板事件的P波和S波、以及每個網(wǎng)格對于不同觀測站的P波校正時差和S波校正時差，接著，獲取所述模板事件的各P波互相關(guān)系數(shù)和各S波互相關(guān)系數(shù)，并根據(jù)所述模板事件的各P波互相關(guān)系數(shù)和各S波的互相關(guān)系數(shù)、每個所述網(wǎng)格的P波校正時差和每個所述網(wǎng)格的S波校正時差，確定每個所述網(wǎng)格的在不同時刻的第一疊加互相關(guān)系數(shù)，最后根據(jù)每個第一疊加互相關(guān)系數(shù)和第一預(yù)設(shè)閾值，獲得所述模板事件對應(yīng)的微地震事件的位置。本實施例的方法使用每個網(wǎng)格的P波校正時差對模板事件的各P波互相關(guān)系數(shù)進(jìn)行校正，使用每個網(wǎng)格的S波校正時差對模板事件的各S波互相關(guān)系數(shù)進(jìn)行校正，并將校正后的P波互相關(guān)系數(shù)和S波互相關(guān)系數(shù)進(jìn)行疊加，獲得每個網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)，并獲得滿足第一預(yù)設(shè)閾值的第一疊加互相關(guān)系數(shù)所對應(yīng)的網(wǎng)格的位置和時刻，進(jìn)而實現(xiàn)對低信噪比的微地震事件的準(zhǔn)確檢測和定位。

圖2為本發(fā)明提供的微地震事件的檢測和定位方法實施例二的流程示意圖。在上述實施例的基礎(chǔ)上，本實施例涉及的是檢測和定位裝置根據(jù)所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)、所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)、每個所述網(wǎng)格與不同觀測站之間的P波校正時差和每個所述網(wǎng)格與不同觀測站之間的S波校正時差，確定每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)的具體過程。如圖2所示，上述S102具體可以包括：

S201、根據(jù)每個所述網(wǎng)格與不同觀測站之間的P波校正時差和所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一P波互相關(guān)系數(shù)。

S202、根據(jù)每個所述網(wǎng)格與不同觀測站之間的S波校正時差和所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)。

具體的，檢測和定位裝置使用每個網(wǎng)格與不同觀測站之間的P波校正時差校正模板事件在不同時刻與不同觀測站之間的P波互相關(guān)系數(shù)，獲得每個網(wǎng)格在不同時刻與不同觀測站之間的第一P波互相關(guān)系數(shù)，同理使用每個網(wǎng)格與不同觀測站之間的S波校正時差校正模板事件在不同時刻與不同觀測站之間的S波互相關(guān)系數(shù)，獲得每個網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)。

在本實施例的一種可能的實現(xiàn)方式中，上述獲取每個網(wǎng)格的P波校正時差和每個網(wǎng)格的S波校正時差，可以根據(jù)公式：

獲得。

其中，表示網(wǎng)格X與第j個觀測站之間的P波(或S波)校正時差，表示第j個觀測站觀測的模板事件的P波(或S波)的實際到時，表示第j個觀測站計算的模板事件的P波(或S波)的理論走時，表示第j個觀測站計算的網(wǎng)格X的P波(或S波)的理論走時，X^M表示模板事件。

可選的，檢測和定位裝置可以根據(jù)公式：

獲得模板事件在不同時刻與不同觀測站之間的P波互相關(guān)系數(shù)和模板事件在不同時刻與不同觀測站之間的S波互相關(guān)系數(shù)。下標(biāo)i表示三分量(分量E、分量N和分量Z)上的第i個分量。表示時間t時模板事件的P波(或S波)與第j個觀測站的原始波形之間的互相關(guān)系數(shù)，u_ij為第j個觀測站檢測到的i分量上的原始數(shù)據(jù)，為第j個觀測站檢測到的i分量上的模板事件的P波(或S波)，T為模板事件的P波(或S波)的長度。

根據(jù)上述兩個公式可以獲得公式和根據(jù)公式獲得網(wǎng)格X在時間與第j個觀測站之間的第一P波互相關(guān)系數(shù)，根據(jù)公式獲得網(wǎng)格X在時間與第j個觀測站之間的第一S波互相關(guān)系數(shù)。

S203、根據(jù)每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一P波互相關(guān)系數(shù)和每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)。

具體的，檢測和定位裝置將上述每個網(wǎng)格在不同時刻與不同觀測站之間的第一P波互相關(guān)系數(shù)和每個網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)進(jìn)行疊加，求其平均值或者加權(quán)平均值，將該平均值或者加權(quán)平均值作為每個網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)。

在本實施例的一種可能的實現(xiàn)方式中，上述S203具體可以是：

檢測和定位裝置根據(jù)公式獲得每個網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)。

其中，所述S₁(t,X)為所述網(wǎng)格X在時間t的第一疊加互相關(guān)系數(shù)，所述N為檢波器的個數(shù)，所述為所述網(wǎng)格X與第j個觀測站之間的P波校正時差，所述為所述網(wǎng)格X與第j個觀測站之間的S波校正時差，所述為所述網(wǎng)格X在時間與第j個觀測站的原始數(shù)據(jù)之間的第一P波互相關(guān)系數(shù)，所述為所述網(wǎng)格X在時間與第j個觀測站的原始數(shù)據(jù)之間的第一S波互相關(guān)系數(shù)。

本實施例的技術(shù)方案，根據(jù)上述各公式獲得每個網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)。

本發(fā)明實施例提供的微地震事件的檢測和定位方法，通過每個所述網(wǎng)格與不同觀測站之間的P波校正時差校正模板事件在不同時刻與各觀測站之間的P波互相關(guān)系數(shù)，獲得每個網(wǎng)格在不同時刻與不同觀測站之間的第一P波互相關(guān)系數(shù)，并通過每個所述網(wǎng)格的S波校正時差校正模板事件在不同時刻與各觀測站之間的S波互相關(guān)系數(shù)，獲得每個網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)。將每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一P波互相關(guān)系數(shù)和每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)進(jìn)行疊加求平均，獲得每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)，進(jìn)而實現(xiàn)對第一疊加互相關(guān)系數(shù)的準(zhǔn)確獲取，從而使得基于第一疊加互相關(guān)系數(shù)的微地震事件的檢測和定位更加準(zhǔn)確。

圖3為本發(fā)明提供的微地震事件的檢測和定位方法實施例三的流程示意圖。在上述實施例的基礎(chǔ)上，本實施例涉及的是檢測和定位裝置根據(jù)每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)和第一預(yù)設(shè)閾值，獲得所述模板事件對應(yīng)的微地震事件的位置和時刻的具體過程。如圖3所示，上述S103具體可以包括：

S301、根據(jù)每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)、第一預(yù)設(shè)閾值和所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)。

其中，所述P波約束互相關(guān)系數(shù)為對所述模板事件的P波互相關(guān)系數(shù)進(jìn)行方位角(或者出射角)約束后獲得的。

具體的，檢測和定位裝置獲取滿足第一預(yù)設(shè)閾值的第一疊加互相關(guān)系數(shù)對應(yīng)的時刻，獲取原始數(shù)據(jù)中該時刻對應(yīng)的事件與各觀測站之間的實際出射角。接著，獲取每個網(wǎng)格與各觀測站之間的理論出射角，其中出射角的計算方法為本領(lǐng)域技術(shù)人員的公知常識，在此不再贅述。根據(jù)上述各實際出射角和每個網(wǎng)格的理論出射角對模板事件在不同時刻與不同觀測站之間的P波互相關(guān)系數(shù)進(jìn)行約束，獲得每個網(wǎng)格的在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)。

可選的，檢測和定位裝置獲取滿足第一預(yù)設(shè)閾值的第一疊加互相系數(shù)對應(yīng)的時刻，獲取原始數(shù)據(jù)中該時刻對應(yīng)的事件與各觀測站之間的實際方位角，用該實際方位角約束模板事件的P波互相關(guān)系數(shù)，獲得每個網(wǎng)格的P波約束互相關(guān)系數(shù)?？蛇x的，檢測和定位裝置還可以根據(jù)其他的方法獲得每個網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)，本實施例對獲取每個網(wǎng)格的P波約束互相關(guān)系數(shù)的方法不做限制，只要是根據(jù)每個第一疊加互相關(guān)系數(shù)、第一預(yù)設(shè)閾值和模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)獲得的即可。

S302、根據(jù)每個所述網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)、每個所述網(wǎng)格與不同觀測站之間的P波校正時差和每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻的第二疊加互相關(guān)系數(shù)。

具體的，檢測和定位裝置使用每個網(wǎng)格與不同觀測站之間的P波校正時差校正每個網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)，例如用網(wǎng)格X與觀測站A之間的P波校正時差，校正網(wǎng)格X在不同時刻與觀測站A之間的P波約束互相關(guān)系數(shù)。將校正過的每個網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)和上述每個網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)進(jìn)行疊加求平均等處理，獲得每個網(wǎng)格在不同時刻的第二疊加互相關(guān)系數(shù)。

S303、根據(jù)每個所述網(wǎng)格在不同時刻的所述第二疊加互相關(guān)系數(shù)和第二預(yù)設(shè)閾值，獲得所述模板事件對應(yīng)的微地震事件的位置和時刻。

具體的，檢測和定位裝置從各第二疊加互相關(guān)系數(shù)中獲取滿足第二預(yù)設(shè)閾值的至少一個目標(biāo)第二疊加互相關(guān)系數(shù)，獲取每個目標(biāo)第二疊加互相關(guān)系數(shù)對應(yīng)的網(wǎng)格的位置和時刻，將該些網(wǎng)格的位置作為模板事件對應(yīng)的微地震事件的位置，將對應(yīng)的時間作為微地震事件發(fā)生的時間，例如獲取的一個目標(biāo)第二疊加互相關(guān)系數(shù)對應(yīng)的時間為t₁，其對應(yīng)的網(wǎng)格的位置為Y，則可獲得模板事件對應(yīng)的一個微地震事件的發(fā)生的時間為t₁，發(fā)生地震的位置為Y。

上述第二預(yù)設(shè)閾值的大小根據(jù)實際需要進(jìn)行設(shè)定，本實施例對此不做限制。

可選的，本實施例還可以是根據(jù)每個所述網(wǎng)格在不同時刻與不同觀測站之間的S波約束互相關(guān)系數(shù)、每個所述網(wǎng)格與不同觀測站之間的S波校正時差和每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一P波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻的第二疊加互相關(guān)系數(shù)，其具體過程同上述步驟相同。

本發(fā)明實施例提供的微地震事件的檢測和定位方法，通過獲取每個網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)，并根據(jù)每個網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)、每個網(wǎng)格與不同觀測站之間的P波校正時差和每個網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)，獲得每個網(wǎng)格在不同時刻的第二疊加互相關(guān)系數(shù)，并根據(jù)每個網(wǎng)格在不同時刻的第二疊加互相關(guān)系數(shù)和第二預(yù)設(shè)閾值，獲得模板事件對應(yīng)的微地震事件的位置和時刻，進(jìn)一步提高了低信噪比的微地震事件的定位準(zhǔn)確性。

圖4為本發(fā)明提供的微地震事件的檢測和定位方法實施例四的流程示意圖。在上述實施例的基礎(chǔ)上，本實施例涉及的是檢測和定位裝置根據(jù)每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)、第一預(yù)設(shè)閾值和所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)的具體過程。即如圖4所示，上述S302具體可以包括：

S401、從每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)中獲取滿足所述第一預(yù)設(shè)閾值的至少一個參考第一疊加互相關(guān)系數(shù)。

S402、獲取每個參考第一疊加互相關(guān)系數(shù)所對應(yīng)的參考時刻，并從各觀測站對應(yīng)的原始數(shù)據(jù)中獲取每個參考時刻對應(yīng)的參考事件。

S403、確定每個參考事件與每個觀測站之間的實際方位角。

具體的，檢測和定位裝置從上述獲得每個網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)中選出滿足第一預(yù)設(shè)閾值的至少一個參考第一疊加互相關(guān)系數(shù)，獲取每個參考第一疊加互相關(guān)系數(shù)所對應(yīng)的時間，將該時間記為參考時刻。從原始數(shù)據(jù)中獲取各參考時刻對應(yīng)的事件，將該事件記為參考事件，進(jìn)而確定各參考事件與每個觀測站之間的實際方位角?？蛇x的，檢測和定位裝置可以根據(jù)三分量波形法獲取各參考事件與每個觀測站之間的實際方位角。

進(jìn)一步的，在本實施例中，為了進(jìn)一步提高了定位的準(zhǔn)確性，則本實施例在上述A403之后，對模板事件的網(wǎng)格結(jié)構(gòu)進(jìn)行了重新劃分，獲得模板事件的新的網(wǎng)格結(jié)構(gòu)，該新的網(wǎng)格結(jié)構(gòu)與之前的網(wǎng)格結(jié)構(gòu)的整體大小相同。但該新的網(wǎng)格結(jié)構(gòu)中的網(wǎng)格更加細(xì)，包括的網(wǎng)格數(shù)目較多。進(jìn)而根據(jù)上述方法可以獲得的新的網(wǎng)格的第二疊加互相關(guān)系數(shù)的數(shù)目增多，從而使得基于該較多數(shù)目的第二疊加互相關(guān)系數(shù)獲得的微地震事件的位置更加精確。

在本實施例中，獲取模板事件的新的網(wǎng)絡(luò)結(jié)構(gòu)的方法可以是，首先獲取模板事件的位置，以該模板事件的位置為中心，獲取該模板事件的新的網(wǎng)格結(jié)構(gòu)。

S404、獲取每個所述網(wǎng)格與每個觀測站之間的理論方位角。

其中，檢測和定位裝置可以根據(jù)現(xiàn)有技術(shù)獲取每個網(wǎng)格與每個觀測站之間的理論方位角。

S405、根據(jù)每個觀測站與不同參考時刻的參考事件的實際方位角、每個網(wǎng)格與每個觀測站之間的理論方位角和所述模板事件在不同的第一時間段內(nèi)與每個觀測站之間的P波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同的第一時間段內(nèi)與不同觀測站之間的P波約束互相關(guān)系數(shù)。

其中，所述第一時間段的長度等于預(yù)設(shè)的時間窗口長度，所述不同的第一時間段包括不同的參考時刻，如圖1a-圖1c所示，預(yù)設(shè)的時間窗口長度為0.03s。

在本實施例中，檢測和定位裝置使用不同參考時刻對應(yīng)的參考事件與各觀測站之間的實際方位角和每個網(wǎng)格與各觀測站之間的理論方位角，約束模板事件在不同時刻與每個觀測站之間的P波互相關(guān)系數(shù)，獲得每個網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)。進(jìn)一步的，為了提高計算速度，本實施例的檢測和定位裝置只獲取每個網(wǎng)格在不同的第一時間段內(nèi)與不同觀測站之間的P波約束互相關(guān)系數(shù)，每個第一時間段的長度與預(yù)設(shè)的時間窗口長度相同，且每個第一時間段包括對應(yīng)的參考時刻。例如，參考時刻為t₀，預(yù)設(shè)的時間窗口長度為t_w，則參考時刻t₀對應(yīng)的第一時間為(t₀-t_w/2，t₀+t_w/2)。

可選的，檢測和定位裝置可以使用每個參考事件與不同觀測站之間的實際方位角與每個網(wǎng)格與每個觀測站之間的理論方位角之間的差值來約束模板事件與不同觀測站之間的P波互相關(guān)系數(shù)，進(jìn)而獲得每個所述網(wǎng)格在不同的第一時間段內(nèi)與不同觀測站之間的P波約束互相關(guān)系數(shù)。

在本實施例的一種可能的實現(xiàn)方式中，上述S406具體可以是：

根據(jù)公式獲得每個所述網(wǎng)格在不同的第一時間段內(nèi)與不同觀測站之間的P波約束互相關(guān)系數(shù)；

其中，所述為所述網(wǎng)格X在時間t與第j個觀測站之間的P波約束互相關(guān)系數(shù)，所述為所述模板事件在時間t與第j個觀測站之間的P波互相關(guān)系數(shù)，所述為所述網(wǎng)格點X與第j個觀測站之間的理論方位角，所述為所述參考時刻t₀對應(yīng)的參考事件與第j個觀測站之間的實際方位角，所述t_w為預(yù)設(shè)的時間窗口長度。

在本實施例中，當(dāng)根據(jù)上述方法獲得每個參考時刻對應(yīng)的第一時間段內(nèi)每個所述網(wǎng)格與每個觀測站之間的P波約束互相關(guān)系數(shù)時，則上述S302可以用下面S406替換：

S406、根據(jù)每個所述網(wǎng)格在不同的第一時間段內(nèi)與不同觀測站之間的P波約束互相關(guān)系數(shù)、每個所述網(wǎng)格與不同觀測站之間的P波校正時差和每個所述網(wǎng)格在不同的第一時間內(nèi)與不同觀測站之間的第一S波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同的第一時間段內(nèi)的第二疊加互相關(guān)系數(shù)。

具體的，檢測和定位裝置使用每個網(wǎng)格與各觀測站之間的P波校正時差校正每個所述網(wǎng)格在不同的第一時間段內(nèi)與不同觀測站之間的P波約束互相關(guān)系數(shù)，接著，將校正后的每個所述網(wǎng)格在不同的第一時間段內(nèi)與不同觀測站之間的P波約束互相關(guān)系數(shù)與上述每個網(wǎng)格與各觀測站之間的第一S波互相關(guān)系數(shù)進(jìn)行疊加求平均，即可獲得每個網(wǎng)格在不同的第一時間段內(nèi)的第二疊加互相關(guān)系數(shù)。

在本實施例的一種可能的實現(xiàn)方式中，上述S407具體可以是：

檢測和定位裝置根據(jù)公式

每個網(wǎng)格在不同的第一時間段內(nèi)的第二疊加互相關(guān)系數(shù)。

其中，所述S₂(t,X)為網(wǎng)格X在時間t與第j個觀測站之間第二疊加互相關(guān)系數(shù)，

在本實施例的一種可能的實現(xiàn)方式中，上述S303還具體可以包括：

S407、從每個所述網(wǎng)格在不同的第一時間段內(nèi)的各第二疊加互相關(guān)系數(shù)中，獲取每個第一時間段中滿足所述第二預(yù)設(shè)閾值的至少一個參考第二疊加互相關(guān)系數(shù)。

例如，當(dāng)S303中的第二預(yù)設(shè)閾值為求不同的第一時間段內(nèi)的第二疊加互相關(guān)系數(shù)中的最大值時，則檢測和定位裝置從每個第一時間段對應(yīng)的各第二疊加互相關(guān)系數(shù)中，獲取一個最大的第二疊加互相關(guān)系數(shù)，將該最大的第二疊加互相關(guān)系數(shù)記為該第一時間段對應(yīng)的參考第二疊加互相關(guān)系數(shù)。

S408、獲得每個參考第二疊加互相關(guān)系數(shù)所對應(yīng)的第一時刻和所述每個參考第二疊加互相關(guān)系數(shù)所對應(yīng)的網(wǎng)格的第一位置。

S409、根據(jù)所述第一位置和所述第一時刻，確定所述模板事件對應(yīng)的微地震事件的位置和時刻。

具體的，檢測和定位裝置獲取每個參考第二疊加互相關(guān)系數(shù)所對應(yīng)的網(wǎng)格的位置第一位置和每個參考第二疊加互相關(guān)系數(shù)所對應(yīng)的第一時刻。例如從參考第二疊加互相關(guān)系數(shù)S₂(t₂,X)中獲取網(wǎng)格X，將網(wǎng)格X的位置記為該模板事件對應(yīng)的微地震事件的一個位置，將時間t₂記為微地震事件X發(fā)生的時刻。根據(jù)上述方法，可以獲得該模板事件的多個微地震事件的位置和時刻，進(jìn)而實現(xiàn)對該模板事件對應(yīng)的低信噪比的微地震事件的準(zhǔn)確檢測和定位。

本發(fā)明實施例提供的微地震事件的檢測和定位方法，使用方位角來約束模板事件與各觀測站之間的P波互相關(guān)系數(shù)，獲得每個網(wǎng)格在不同的第一時間段內(nèi)與每個觀測站之間的P波約束互相關(guān)系數(shù)，根據(jù)每個網(wǎng)格在不同的第一時間段內(nèi)與每個觀測站之間的P波約束互相關(guān)系數(shù)和每個網(wǎng)格在不同的第一時間段內(nèi)與每個觀測站之間的第一S互相關(guān)系數(shù)，獲得每個網(wǎng)格在不同的第一時間段內(nèi)的第二疊加互相關(guān)系數(shù)，再根據(jù)每個第二疊加互相關(guān)系數(shù)和第二預(yù)設(shè)閾值獲得模板事件對應(yīng)的微地震事件的位置和時刻，進(jìn)一步提高了低信噪比的微地震事件的定位準(zhǔn)確性。

本領(lǐng)域普通技術(shù)人員可以理解：實現(xiàn)上述各方法實施例的全部或部分步驟可以通過程序指令相關(guān)的硬件來完成。前述的程序可以存儲于一計算機(jī)可讀取存儲介質(zhì)中。該程序在執(zhí)行時，執(zhí)行包括上述各方法實施例的步驟；而前述的存儲介質(zhì)包括：ROM、RAM、磁碟或者光盤等各種可以存儲程序代碼的介質(zhì)。

圖5為本發(fā)明提供的微地震事件的檢測和定位裝置實施例一的結(jié)構(gòu)示意圖。本實施例的檢測和定位裝置可以繼承在現(xiàn)有的地震監(jiān)測裝置中，可選的還可以單獨的裝置，本實施例的檢測和定位裝置可以軟件和/或硬件實現(xiàn)。如圖5所示，本實施例的檢測和定位裝置可以包括：

獲取模塊10，用于獲取模板事件的網(wǎng)格結(jié)構(gòu)中的每個網(wǎng)格與不同觀測站之間的P波校正時差、所述網(wǎng)格結(jié)構(gòu)中的每個網(wǎng)格與不同觀測站之間的S波校正時差、所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)和所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)；其中，所述網(wǎng)格結(jié)構(gòu)包括多個網(wǎng)格，所述模板事件為觀測站的原始數(shù)據(jù)中信噪比大于預(yù)設(shè)閾值的事件。

處理模塊20，用于根據(jù)所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)、所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)、每個所述網(wǎng)格與不同觀測站之間的P波校正時差和每個所述網(wǎng)格與不同觀測站之間的S波校正時差，確定每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)，所述第一疊加互相系數(shù)用于表征所述網(wǎng)格為震源的概率。

確定模塊30，用于根據(jù)每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)和第一預(yù)設(shè)閾值，確定所述模板事件對應(yīng)的微地震事件的位置和時刻。

本實施例的裝置，可以用于執(zhí)行上述所示方法實施例的技術(shù)方案，其實現(xiàn)原理和技術(shù)效果類似，此處不再贅述。

圖6為本發(fā)明提供的微地震事件的檢測和定位裝置實施例二的結(jié)構(gòu)示意圖。在上述實施例的基礎(chǔ)上，本實施例的處理模塊20具體可以包括：第一獲取單元201、第二獲取單元202和第一疊加單元203；

上述第一獲取單元201，用于根據(jù)每個所述網(wǎng)格與不同觀測站之間的P波校正時差和所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一P波互相關(guān)系數(shù)。

上述第二獲取單元202，用于根據(jù)每個所述網(wǎng)格與不同觀測站之間的S波校正時差和所述模板事件的S波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的S波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)；

上述第一疊加單元203，用于根據(jù)每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一P波互相關(guān)系數(shù)和每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)。

本實施例的裝置，可以用于執(zhí)行上述所示方法實施例的技術(shù)方案，其實現(xiàn)原理和技術(shù)效果類似，此處不再贅述。

圖7為本發(fā)明提供的微地震事件的檢測和定位裝置實施例三的結(jié)構(gòu)示意圖。在上述實施例的基礎(chǔ)上，本實施例的確定模塊30具體可以包括：第三獲取單元301、第二疊加單元302和確定單元303；

上述第三獲取單元301，用于根據(jù)每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)、第一預(yù)設(shè)閾值和所述模板事件的P波在不同時刻與每個觀測站的原始數(shù)據(jù)之間的P波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)；所述P波約束互相關(guān)系數(shù)為對所述模板事件的P波互相關(guān)系數(shù)進(jìn)行方位角約束后獲得的。

上述第二疊加單元302，用于根據(jù)每個所述網(wǎng)格在不同時刻與不同觀測站之間的P波約束互相關(guān)系數(shù)、每個所述網(wǎng)格與不同觀測站之間的P波校正時差和每個所述網(wǎng)格在不同時刻與不同觀測站之間的第一S波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同時刻的第二疊加互相關(guān)系數(shù)。

上述確定單元303，用于根據(jù)每個所述網(wǎng)格在不同時刻的所述第二疊加互相關(guān)系數(shù)和第二預(yù)設(shè)閾值，獲得所述模板事件對應(yīng)的微地震事件的位置和時刻。

本實施例的裝置，可以用于執(zhí)行上述所示方法實施例的技術(shù)方案，其實現(xiàn)原理和技術(shù)效果類似，此處不再贅述。

在本發(fā)明的一種可能的實現(xiàn)方式中，上述第三獲取單元301，具體用于：

從每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)中獲取滿足所述第一預(yù)設(shè)閾值的至少一個參考第一疊加互相關(guān)系數(shù)；

獲取每個參考第一疊加互相關(guān)系數(shù)所對應(yīng)的參考時刻，并從各觀測站對應(yīng)的原始數(shù)據(jù)中獲取每個參考時刻對應(yīng)的參考事件；

確定每個參考事件與每個觀測站之間的實際方位角，以及每個所述網(wǎng)格與每個觀測站之間的理論方位角；

根據(jù)每個觀測站與不同參考時刻的參考事件的實際方位角、每個網(wǎng)格與每個觀測站之間的理論方位角和所述模板事件在不同的第一時間段內(nèi)與每個觀測站之間的P波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同的第一時間段內(nèi)與不同觀測站之間的P波約束互相關(guān)系數(shù)，其中，所述不同的第一時間段包括不同的參考時刻。

上述第二疊加單元302，具體用于根據(jù)每個所述網(wǎng)格在不同的第一時間段內(nèi)與不同觀測站之間的P波約束互相關(guān)系數(shù)、每個所述網(wǎng)格與不同觀測站之間的P波校正時差和每個所述網(wǎng)格在不同的第一時間內(nèi)與不同觀測站之間的第一S波互相關(guān)系數(shù)，獲得每個所述網(wǎng)格在不同的第一時間段內(nèi)的第二疊加互相關(guān)系數(shù)。

上述確定單元303，具體用于：

從每個所述網(wǎng)格在不同的第一時間段內(nèi)的各第二疊加互相關(guān)系數(shù)中，獲取每個第一時間段中滿足所述第二預(yù)設(shè)閾值的至少一個參考第二疊加互相關(guān)系數(shù)；

獲得每個參考第二疊加互相關(guān)系數(shù)所對應(yīng)的第一時刻和所述每個參考第二疊加互相關(guān)系數(shù)所對應(yīng)的網(wǎng)格的第一位置；

根據(jù)所述第一位置和所述第一時刻，確定所述模板事件對應(yīng)的微地震事件的位置和時刻。

在本發(fā)明的另一種可能的實現(xiàn)方式中，上述獲取模塊，還用于在獲取每個所述網(wǎng)格與每個觀測站之間的理論方位角之前，獲取所述模板事件的新的網(wǎng)格結(jié)構(gòu)，所述新的網(wǎng)格結(jié)構(gòu)包括多個新的網(wǎng)格，所述新的網(wǎng)格的尺寸小于所述網(wǎng)格的尺寸；對應(yīng)的，上述第三獲取單元，還具體用于獲取每個所述新的網(wǎng)格與每個所述觀測站之間的理論方位角。

進(jìn)一步的，上述第一疊加單元203，具體用于：

根據(jù)公式獲得每個所述網(wǎng)格在不同時刻的第一疊加互相關(guān)系數(shù)；

其中，所述S₁(t,X)為所述網(wǎng)格X在時間t的第一疊加互相關(guān)系數(shù)，所述N為觀測站的個數(shù)，所述為所述網(wǎng)格X與第j個觀測站之間的P波校正時差，所述為所述網(wǎng)格X與所述第j個觀測站之間的S波校正時差，所述為所述網(wǎng)格X在時間與所述第j個觀測站之間的第一P波互相關(guān)系數(shù)，所述為所述網(wǎng)格X在時間與所述第j個觀測站之間的第一S波互相關(guān)系數(shù)。

進(jìn)一步的，上述第三獲取單元301，具體用于：

根據(jù)公式獲得每個所述網(wǎng)格在不同的第一時間段內(nèi)與不同觀測站之間的P波約束互相關(guān)系數(shù)；

其中，所述為所述網(wǎng)格X在時間t與所述第j個觀測站之間的P波約束互相關(guān)系數(shù)，所述為所述模板事件在時間t與所述第j個觀測站之間的P波互相關(guān)系數(shù)，所述為所述網(wǎng)格點X與所述第j個觀測站之間的理論方位角，所述為所述t₀時刻對應(yīng)的參考事件與所述第j個觀測站之間的實際方位角，所述t₀為參考時刻，所述t_w為所述預(yù)設(shè)的時間窗口長度。

本實施例的裝置，可以用于執(zhí)行上述所示方法實施例的技術(shù)方案，其實現(xiàn)原理和技術(shù)效果類似，此處不再贅述。

最后應(yīng)說明的是：以上各實施例僅用以說明本發(fā)明的技術(shù)方案，而非對其限制；盡管參照前述各實施例對本發(fā)明進(jìn)行了詳細(xì)的說明，本領(lǐng)域的普通技術(shù)人員應(yīng)當(dāng)理解：其依然可以對前述各實施例所記載的技術(shù)方案進(jìn)行修改，或者對其中部分或者全部技術(shù)特征進(jìn)行等同替換；而這些修改或者替換，并不使相應(yīng)技術(shù)方案的本質(zhì)脫離本發(fā)明各實施例技術(shù)方案的范圍。