專利名稱：一種檢測地震信號的控制方法及裝置的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及地震信號采集領(lǐng)域，尤其涉及一種檢測地震信號的控制方法及裝置。
背景技術(shù)：
地震勘探主要分采集、處理和解釋三大部分，地震數(shù)據(jù)采集時(shí)需要拾取地震信號，地震波初至的拾取廣泛應(yīng)用在表層調(diào)查方法和地質(zhì)構(gòu)造的探測中，初至的拾取對礦物資源和石油資源的開發(fā)，都起著重要的作用。在地震觀測中，地震事件的檢測識別及初至到時(shí)拾取是一項(xiàng)重要而又繁瑣的工作，當(dāng)前出現(xiàn)了一種基于地震波形峭度和偏斜度的檢測方法。對于一個(gè)有限長度的離散實(shí)數(shù)序列{xj ,其偏斜度S (skewness)和峭度K(kurtosis,又譯峰度)定義為
權(quán)利要求
1.一種檢測地震信號的控制方法，其特征在于，包括 獲取地震波的N個(gè)樣點(diǎn)的離散序列IxJ，當(dāng)前時(shí)窗長度為M ； 計(jì)算當(dāng)前時(shí)窗長度為M的時(shí)窗內(nèi)樣點(diǎn)的算術(shù)和、平方和、立方和、4次方和；將時(shí)窗向前移動一個(gè)樣點(diǎn)，計(jì)算新時(shí)窗內(nèi)樣點(diǎn)的各值，是在上一個(gè)時(shí)窗的值基礎(chǔ)上，減去上一個(gè)時(shí)窗第I個(gè)樣點(diǎn)的值，并加上當(dāng)前時(shí)窗最后一個(gè)樣點(diǎn)的值，依次類推，得到所有時(shí)窗N個(gè)樣點(diǎn)的峭度和偏斜度的數(shù)值；并通過得到的所有時(shí)窗N個(gè)樣點(diǎn)的峭度和偏斜度的數(shù)值得到峭度曲線和偏斜度曲線； 通過得到的地震波的峭度曲線和偏斜度曲線，判斷地震波到時(shí)的位置。
2.如權(quán)利要求I所述的方法，其特征在于， 所述得到所有時(shí)窗N個(gè)樣點(diǎn)的峭度K的數(shù)值的步驟，包括 計(jì)算當(dāng)前時(shí)窗長度為M的時(shí)窗內(nèi)樣點(diǎn)的算術(shù)和、平方和、立方和、4次方和，其中，將當(dāng)前時(shí)窗長度為M的時(shí)窗內(nèi)離散序列的算術(shù)和、平方和、立方和、4次方和分別記為Smj, j=l, 2，3，4,公式如下所示
3.如權(quán)利要求I所述的方法，其特征在于， 所述得到所有時(shí)窗N個(gè)樣點(diǎn)的偏斜度S的數(shù)值的步驟，包括 計(jì)算當(dāng)前時(shí)窗長度為M的時(shí)窗內(nèi)樣點(diǎn)的算術(shù)和、平方和、立方和，其中，將當(dāng)前時(shí)窗長度為M的時(shí)窗內(nèi)離散序列的算術(shù)和、平方和、立方和分別記為smy j=l, 2，3，公式如下所示
4.如權(quán)利要求I所述的方法，其特征在于， 所述通過得到的地震波的峭度曲線和偏斜度曲線，判斷地震波到時(shí)的位置的步驟，包括 通過地震波形峭度和偏斜度檢測地震波到時(shí)來判斷峭度曲線和偏斜度曲線的極大值，根據(jù)極大值點(diǎn)前曲線斜率最大的位置判斷地震波到時(shí)的位置。
5.一種檢測地震信號的控制裝置，其特征在于，包括 獲取地震波樣點(diǎn)模塊，用于獲取地震波的N個(gè)樣點(diǎn)的離散序列IxJ，當(dāng)前時(shí)窗長度為M ； 計(jì)算模塊，用于計(jì)算當(dāng)前時(shí)窗長度為M的時(shí)窗內(nèi)樣點(diǎn)的算術(shù)和、平方和、立方和、4次方和；將時(shí)窗向前移動一個(gè)樣點(diǎn)，計(jì)算新時(shí)窗內(nèi)樣點(diǎn)的各值，是在上一個(gè)時(shí)窗的值基礎(chǔ)上，減去上一個(gè)時(shí)窗第I個(gè)樣點(diǎn)的值，并加上當(dāng)前時(shí)窗最后一個(gè)樣點(diǎn)的值，依次類推，得到所有時(shí)窗N個(gè)樣點(diǎn)的峭度和偏斜度的數(shù)值； 存儲模塊，用于存儲數(shù)據(jù) '及判斷模塊，用于通過得到的所有時(shí)窗N個(gè)樣點(diǎn)的峭度和偏斜度的數(shù)值得到峭度曲線和偏斜度曲線，并通過得到的地震波的峭度曲線和偏斜度曲線，判斷地震波到時(shí)的位置。
6.如權(quán)利要求5所述的裝置，其特征在于， 所述計(jì)算模塊得到所有時(shí)窗N個(gè)樣點(diǎn)的峭度K的數(shù)值，是指 所述計(jì)算模塊計(jì)算當(dāng)前時(shí)窗長度為M的時(shí)窗內(nèi)樣點(diǎn)的算術(shù)和、平方和、立方和、4次方和，其中，將當(dāng)前時(shí)窗長度為M的時(shí)窗內(nèi)離散序列的算術(shù)和、平方和、立方和、4次方和分別記為sm,-，j=l, 2，3，4,公式如下所示
7.如權(quán)利要求5所述的裝置，其特征在于， 所述計(jì)算模塊得到所有時(shí)窗N個(gè)樣點(diǎn)的偏斜度S的數(shù)值，是指 所述計(jì)算模塊計(jì)算當(dāng)前時(shí)窗長度為M的時(shí)窗內(nèi)樣點(diǎn)的算術(shù)和、平方和、立方和，其中，將當(dāng)前時(shí)窗長度為M的時(shí)窗內(nèi)離散序列的算術(shù)和、平方和、立方和分別記為Sn^，j=l, 2，3，公式如下所示
8.如權(quán)利要求5所述的裝置，其特征在于， 所述判斷模塊通過得到的地震波的峭度曲線和偏斜度曲線，判斷地震波到時(shí)的位置，是指所述判斷模塊通過地震波形峭度和偏斜度檢測地震波到時(shí)來判斷峭度曲線和偏斜度曲線的極大值，根據(jù)極大值點(diǎn)前曲線斜率最大的位置判斷地震波到時(shí)的位置。
全文摘要
一種檢測地震信號的控制方法及系統(tǒng)，包括獲取地震波的N個(gè)樣點(diǎn)的離散序列；計(jì)算當(dāng)前時(shí)窗內(nèi)樣點(diǎn)的算術(shù)和、平方和、立方和、4次方和；將時(shí)窗向前移動一個(gè)樣點(diǎn)，計(jì)算新時(shí)窗內(nèi)樣點(diǎn)的各值，是在上一個(gè)時(shí)窗的值基礎(chǔ)上，減去上一個(gè)時(shí)窗第1個(gè)樣點(diǎn)的值，并加上當(dāng)前時(shí)窗最后一個(gè)樣點(diǎn)的值，依次類推，得到所有時(shí)窗N個(gè)樣點(diǎn)的峭度和偏斜度的數(shù)值得到峭度曲線和偏斜度曲線，判斷地震波到時(shí)的位置。本發(fā)明解決了當(dāng)前傳統(tǒng)的基于地震波形峭度和偏斜度的檢測方法存在運(yùn)算量過大導(dǎo)致效率低的問題。
文檔編號G01V1/22GK102707312SQ20121020475
公開日2012年10月3日 申請日期2012年6月15日 優(yōu)先權(quán)日2012年6月15日
發(fā)明者劉勁松 申請人:中國科學(xué)院地質(zhì)與地球物理研究所