恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的方法和裝置制造方法
【專利摘要】本發(fā)明提供了一種恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的方法和裝置，該方法包括：獲取低頻檢波器輸出的第一地震數(shù)據(jù)和模擬檢波器輸出的第二地震數(shù)據(jù)，其中，低頻檢波器和模擬檢波器之間的距離小于預(yù)設(shè)的閾值，低頻檢波器和模擬檢波器接收到的是同源地震信號；根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率；根據(jù)所述最低頻率，對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。本發(fā)明解決了現(xiàn)有技術(shù)中地震數(shù)據(jù)中的低頻信號難以有效恢復(fù)的技術(shù)問題，達(dá)到了有效恢復(fù)地震數(shù)據(jù)中的低頻信號的技術(shù)效果。
【專利說明】恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的方法和裝置
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地震勘測【技術(shù)領(lǐng)域】，特別涉及一種恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的方法和裝置。
【背景技術(shù)】
[0002]低頻信號由于其特有的穿透能力和穩(wěn)定性，在復(fù)雜構(gòu)造成像和地震反演方面都起著非常重要的作用。然而，常規(guī)地震勘探應(yīng)用的模擬檢波器的自然頻率一般為10Hz，對于IOHz以下的地震數(shù)據(jù)信號衰減較為嚴(yán)重，因此，如何實現(xiàn)對地震數(shù)據(jù)中的低頻信號的有效恢復(fù)是應(yīng)用模擬檢波器實現(xiàn)地震勘探的關(guān)鍵。[0003]然而，目前針對如何有效實現(xiàn)對地震數(shù)據(jù)中的低頻信號的恢復(fù)尚未提出有效的解決方案。

【發(fā)明內(nèi)容】

[0004]本發(fā)明提供了一種恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的方法和裝置，以至少解決現(xiàn)有技術(shù)中難以有效恢復(fù)地震數(shù)據(jù)中的低頻信號的技術(shù)問題。
[0005]本發(fā)明實施例提供了一種恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的方法，包括:獲取低頻檢波器輸出的第一地震數(shù)據(jù)和模擬檢波器輸出的第二地震數(shù)據(jù)，其中，低頻檢波器和模擬檢波器之間的距離小于預(yù)設(shè)的閾值，低頻檢波器和模擬檢波器接收到的是同源地震信號；根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率；根據(jù)所述最低頻率，對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。
[0006]在一個實施例中，根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率，包括:確定第一地震數(shù)據(jù)和第二地震數(shù)據(jù)是否具有相同的數(shù)據(jù)格式；如果數(shù)據(jù)格式不相同，則將第一地震數(shù)據(jù)和第二地震數(shù)據(jù)轉(zhuǎn)換為相同的數(shù)據(jù)格式；對相同數(shù)據(jù)格式的第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差進(jìn)行校正；確定進(jìn)行時差校正后的第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比；根據(jù)功率譜密度比得到模擬檢波器可恢復(fù)信號的最低頻率。
[0007]在一個實施例中，對相同數(shù)據(jù)格式的第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差進(jìn)行校正，包括:按照以下公式計算第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差:
[0008]Z(t) = X(t)? Y(t)
[0009]Z (ω) =FFT [Ζ (t)]
[0010]δ (ω) =arctan (imag (Ζ (ω), real (Ζ (ω))
[0011]δ?(ω)=δ (ω)/ω
[0012]其中，Z(t)表示第一地震數(shù)據(jù)和第二地震數(shù)據(jù)的互相關(guān)，X(t)表示第一地震數(shù)據(jù)，Y(t)表示第二地震數(shù)據(jù)，?表示相關(guān)運算，ζ(ω)表示對Z(t)進(jìn)行傅里葉變換后得到的結(jié)果，F(xiàn)FT[]表示傅里葉變換，ω表示角頻率，δ (ω)表示X(t)和Y(t)在角頻率為ω時的相位差，δ?(ω)表示X(t)和Y(t)在角頻率為ω時第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差；根據(jù)計算得到的時差，對第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差進(jìn)行校正。
[0013]在一個實施例中，根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率，包括:確定第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比；繪制第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比曲線，將繪制的功率譜密度比曲線與理論的功率譜密度比曲線進(jìn)行對比，確定模擬檢波器可恢復(fù)信號的最低頻率。
[0014]在一個實施例中，將繪制的功率譜密度比曲線與理論的功率譜密度比曲線進(jìn)行對t匕，確定模擬檢波器可恢復(fù)信號的最低頻率，包括:將繪制的功率譜密度比曲線與理論的功率譜密度比曲線相重合的曲線段對應(yīng)的頻率最小值，作為模擬檢波器可恢復(fù)信號的最低頻率。
[0015]在一個實施例中，在確定第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比之后，所述方法還包括:根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，按照以下公式計算得到用于對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)的反濾波函數(shù):
[0016]hi (O) 二 11 ^PSDRi(O) = I / ^ PSDt{(0) I PSD2(CJ)
[0017]其中，h(co)表示反濾波函數(shù)，PSDR(CO)表示第一地震數(shù)據(jù)和第二地震數(shù)據(jù)的功率譜密度比，PSD1(GJ)表示第一地震數(shù)據(jù)的功率譜密度，PSD2(GJ)表示第二地震數(shù)據(jù)的功率譜密度。
[0018]在一個實施例中，根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率，包括:計算所述模擬檢波器所在的觀測區(qū)中的多個模擬檢波器可恢復(fù)信號的最低頻率；對多個模擬檢波器對應(yīng)的可恢復(fù)信號的最低頻率進(jìn)行擬合，得到應(yīng)用于所述觀測區(qū)中所有模擬檢波器的可恢復(fù)信號的最低頻率；根據(jù)所述最低頻率，對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)，包括:根據(jù)得到的應(yīng)用于所述觀測區(qū)中所有模擬檢波器的可恢復(fù)信號的最低頻率對該觀測區(qū)中的模擬檢波器輸出的第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。`
[0019]在一個實施例中，低頻檢波器和模擬檢波器之間的距離小于預(yù)設(shè)的閾值，包括:在所述模擬檢波器為一個模擬檢波器的情況下，所述低頻檢波器布置在以該模擬檢波器為圓心，半徑小于預(yù)定的閾值的區(qū)域內(nèi)；或者，在所述模擬檢波器為多個模擬檢波器的情況下，所述低頻檢波器布置在所述多個模擬檢波器所在位置所構(gòu)成的幾何圖形的中心位置。
[0020]本發(fā)明實施例還提供了一種恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的裝置，包括:獲取模塊，用于獲取低頻檢波器輸出的第一地震數(shù)據(jù)和模擬檢波器輸出的第二地震數(shù)據(jù)，其中，低頻檢波器和模擬檢波器之間的距離小于預(yù)設(shè)的閾值，低頻檢波器和模擬檢波器接收到的是同源地震信號；確定模塊，用于根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率；恢復(fù)模塊，用于根據(jù)所述最低頻率，對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。
[0021]在一個實施例中，確定模塊包括:第一確定單兀，用于確定第一地震數(shù)據(jù)和第二地震數(shù)據(jù)是否具有相同的數(shù)據(jù)格式；轉(zhuǎn)換單元，用于在確定出數(shù)據(jù)格式不相同的情況下，將第一地震數(shù)據(jù)和第二地震數(shù)據(jù)轉(zhuǎn)換為相同的數(shù)據(jù)格式；校正單元，用于對相同數(shù)據(jù)格式的第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差進(jìn)行校正；第二確定單元，用于確定進(jìn)行時差校正后的第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比；獲取單元，用于根據(jù)功率譜密度比得到模擬檢波器可恢復(fù)信號的最低頻率。
[0022]在一個實施例中，所述校正單元包括:計算子單元，用于按照以下公式計算第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差:
[0023]Z(t) = X(t) ? Y(t)
[0024]Z ( ω ) =FFT [Ζ (t)]
[0025]δ (ω) =arctan (imag (Ζ (ω), real (Ζ (ω))
[0026]δ t ( ω ) = δ ( ω ) / ω
[0027]其中，Z(t)表示第一地震數(shù)據(jù)和第二地震數(shù)據(jù)的互相關(guān)，X(t)表示第一地震數(shù)據(jù)，Y(t)表示第二地震數(shù)據(jù)，?表示相關(guān)運算，Ζ(ω)表示對Z(t)進(jìn)行傅里葉變換后得到的結(jié)果，F(xiàn)FT[]表示傅里葉變換，ω表示角頻率，δ (ω)表示X(t)和Y(t)在角頻率為ω時的相位差，δ?(ω)表示X(t)和Y(t)在角頻率為ω時第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差；校正子單元，用于根據(jù)計算得到的時差，對第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差進(jìn)行校正。
[0028]在一個實施例中，所述確定模塊包括:第三確定單元，用于確定第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比；第四確定單元，用于繪制第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比曲線，將繪制的功率譜密度比曲線與理論的功率譜密度比曲線進(jìn)行對比，確定模擬檢波器可恢復(fù)信號的最低頻率。
[0029]在一個實施例中，所述第四確定單元具體用于將繪制的功率譜密度比曲線與理論的功率譜密度比曲線相重合的曲線段對應(yīng)的頻率最小值，作為模擬檢波器可恢復(fù)信號的最低頻率。
`[0030]在一個實施例中，所述確定模塊還包括:計算單元，用于根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，按照以下公式計算得到用于對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)的反濾波函數(shù):
[0031]h{(0) = I /^PSDR(O)) = I / ^jPSDs(ω) / PSD2(Ca)
[0032]其中，h(co)表示反濾波函數(shù)，PSDR(GJ)表示第一地震數(shù)據(jù)和第二地震數(shù)據(jù)的功率譜密度比，PSD1(GJ)表示第一地震數(shù)據(jù)的功率譜密度，PSD2(GJ)表示第二地震數(shù)據(jù)的功率譜密度。
[0033]在一個實施例中，所述確定模塊具體用于計算所述模擬檢波器所在的觀測區(qū)中的多個模擬檢波器可恢復(fù)信號的最低頻率；對多個模擬檢波器對應(yīng)的可恢復(fù)信號的最低頻率進(jìn)行擬合，得到應(yīng)用于所述觀測區(qū)中所有模擬檢波器的可恢復(fù)信號的最低頻率；所述恢復(fù)模塊具體用于根據(jù)得到的應(yīng)用于所述觀測區(qū)中所有模擬檢波器的可恢復(fù)信號的最低頻率對該觀測區(qū)中的模擬檢波器輸出的第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。
[0034]在一個實施例中，低頻檢波器和模擬檢波器之間的距離小于預(yù)設(shè)的閾值包括:在所述模擬檢波器為一個模擬檢波器的情況下，所述低頻檢波器布置在以該模擬檢波器為圓心，半徑小于預(yù)定的閾值的區(qū)域內(nèi)；或者，在所述模擬檢波器為多個模擬檢波器的情況下，所述低頻檢波器布置在所述多個模擬檢波器所在位置所構(gòu)成的幾何圖形的中心位置。
[0035]在本發(fā)明實施例中，通過布設(shè)的與模擬檢波器之間的距離小于預(yù)設(shè)的閾值的低頻檢波器來接收同源地震信號，然后根據(jù)模擬檢波器與低頻檢波器輸出的地震數(shù)據(jù)的功率譜密度比確定模擬檢波器可恢復(fù)信號的最低頻率，在確定出模擬檢波器可恢復(fù)信號的最低頻率后再對模擬檢波器輸出的地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。因為預(yù)先確定了可恢復(fù)信號的最低頻率，從而使得恢復(fù)出的低頻信號更為準(zhǔn)確。通過上述方式解決了現(xiàn)有技術(shù)中地震數(shù)據(jù)中的低頻信號難以有效恢復(fù)的技術(shù)問題，達(dá)到了有效恢復(fù)低頻地震信號的技術(shù)效果。
【專利附圖】

【附圖說明】
[0036]此處所說明的附圖用來提供對本發(fā)明的進(jìn)一步理解，構(gòu)成本申請的一部分，并不構(gòu)成對本發(fā)明的限定。在附圖中:
[0037]圖1是本發(fā)明實施例的恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的方法的處理流程圖；
[0038]圖2是本發(fā)明實施例的恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的方法的處理流程圖一個具體實例的流程圖；
[0039]圖3是本發(fā)明實施例的低頻檢波器布設(shè)的示意圖；
[0040]圖4是本發(fā)明實施例的恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的裝置的結(jié)構(gòu)示意圖；
[0041]圖5是本發(fā)明實施例的恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的裝置的具體的結(jié)構(gòu)示意圖；
[0042]圖6是本發(fā)明實施例的恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的裝置的具體的結(jié)構(gòu)示意圖。
【具體實施方式】
[0043]為使本發(fā)明的目的、技術(shù)方案和優(yōu)點更加清楚明白，下面結(jié)合實施方式和附圖，對本發(fā)明做進(jìn)一步詳細(xì)說明。在此，本發(fā)明的示意性實施方式及其說明用于解釋本發(fā)明，但并不作為對本發(fā)明的限定。
[0044]本發(fā)明實施例提供了一種恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的方法，如圖1所示，包括以下步驟:
[0045]步驟101:獲取低頻檢波器輸出的第一地震數(shù)據(jù)和模擬檢波器輸出的第二地震數(shù)據(jù)，其中，低頻檢波器和模擬檢波器之間的距離小于預(yù)設(shè)的閾值，低頻檢波器和模擬檢波器接收到的是同源地震信號；
[0046]步驟102:根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率；
[0047]步驟103:根據(jù)所述最低頻率，對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。
[0048]在上述實施例中，通過布設(shè)的與模擬檢波器之間的距離小于預(yù)設(shè)的閾值的低頻檢波器來接收同源地震信號，然后根據(jù)模擬檢波器與低頻檢波器輸出的地震數(shù)據(jù)的功率譜密度比確定模擬檢波器可恢復(fù)信號的最低頻率，在確定出模擬檢波器可恢復(fù)信號的最低頻率后再對模擬檢波器輸出的地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。因為預(yù)先確定了可恢復(fù)信號的最低頻率，從而使得恢復(fù)出的低頻信號更為準(zhǔn)確。通過上述方式解決了現(xiàn)有技術(shù)中地震數(shù)據(jù)中的低頻信號難以有效恢復(fù)的技術(shù)問題，達(dá)到了有效恢復(fù)低頻地震信號的技術(shù)效果。
[0049]考慮到在實施的過程中，模擬檢波器與布設(shè)的低頻檢波器的數(shù)據(jù)格式有可能是不同的，例如布設(shè)的低頻檢波器采用的時序的數(shù)據(jù)格式記錄數(shù)據(jù)，而模擬檢波器采用的是道序的數(shù)據(jù)格式記錄數(shù)據(jù)。在這種情況下，為了實現(xiàn)對數(shù)據(jù)的有效處理，就需要對低頻檢波器和模擬檢波器的數(shù)據(jù)進(jìn)行格式轉(zhuǎn)換，以實現(xiàn)數(shù)據(jù)格式的統(tǒng)一。在一個實施例中，上述步驟102中，根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率可以包括:先確定第一地震數(shù)據(jù)和第二地震數(shù)據(jù)是否具有相同的數(shù)據(jù)格式；如果數(shù)據(jù)格式不相同，則將第一地震數(shù)據(jù)和所述第二地震數(shù)據(jù)轉(zhuǎn)換為相同的數(shù)據(jù)格式，然后，再根據(jù)相同格式的第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率。
[0050]同時，因為不同的儀器在獲取相同數(shù)據(jù)的時候會產(chǎn)生一定的時差，為了使得測量的結(jié)果更為準(zhǔn)確，在計算之前可以先對得到的數(shù)據(jù)進(jìn)行時差校正(又可以稱為TO時校正)，即，先對相同數(shù)據(jù)格式的第一地震數(shù)據(jù)和第二地震之間的時差進(jìn)行校正；然后再確定進(jìn)行時差校正后的第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，根據(jù)該功率譜密度比得到模擬檢波器可恢復(fù)信號的最低頻率。
[0051]在本例中，具體的對第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差進(jìn)行校正可以包括:
[0052]I)按照以下公式計算第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差:
【權(quán)利要求】
1.一種恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的方法，其特征在于，包括: 獲取低頻檢波器輸出的第一地震數(shù)據(jù)和模擬檢波器輸出的第二地震數(shù)據(jù)，其中，低頻檢波器和模擬檢波器之間的距離小于預(yù)設(shè)的閾值，低頻檢波器和模擬檢波器接收到的是同源地震信號； 根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率； 根據(jù)所述最低頻率，對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。
2.如權(quán)利要求1所述的方法，其特征在于，根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率，包括: 確定第一地震數(shù)據(jù)和第二地震數(shù)據(jù)是否具有相同的數(shù)據(jù)格式； 如果數(shù)據(jù)格式不相同，則將第一地震數(shù)據(jù)和第二地震數(shù)據(jù)轉(zhuǎn)換為相同的數(shù)據(jù)格式； 對相同數(shù)據(jù)格式的第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差進(jìn)行校正； 確定進(jìn)行時差校正后的第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比； 根據(jù)功率譜密度比得到模擬檢波器可恢復(fù)信號的最低頻率。
3.如權(quán)利要求2所述的方法，其特征在于，對相同數(shù)據(jù)格式的第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差進(jìn)行校正，包括: 按照以下公式計算第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差:`
Z(I) = X(t) ? Y(t)
Z(ω)=FFT[Ζ (t)]
δ (ω)=arctan(imag(Ζ(ω), real(Ζ(ω))
δ t ( ω ) = δ ( ω ) / ω 其中，Z(t)表示第一地震數(shù)據(jù)和第二地震數(shù)據(jù)的互相關(guān)，X(t)表示第一地震數(shù)據(jù)，Y(t)表示第二地震數(shù)據(jù)，?表示相關(guān)運算，ζ(ω)表示對Z(t)進(jìn)行傅里葉變換后得到的結(jié)果，F(xiàn)FT[]表示傅里葉變換，ω表示角頻率，δ (ω)表示X(t)和Y(t)在角頻率為ω時的相位差，δ?(ω)表示X(t)和Y(t)在角頻率為ω時第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差； 根據(jù)計算得到的時差，對第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差進(jìn)行校正。
4.如權(quán)利要求1所述的方法，其特征在于，根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率，包括: 確定第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比； 繪制第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比曲線，將繪制的功率譜密度比曲線與理論的功率譜密度比曲線進(jìn)行對比，確定模擬檢波器可恢復(fù)信號的最低頻率。
5.如權(quán)利要求4所述的方法，其特征在于，將繪制的功率譜密度比曲線與理論的功率譜密度比曲線進(jìn)行對比，確定模擬檢波器可恢復(fù)信號的最低頻率，包括: 將繪制的功率譜密度比曲線與理論的功率譜密度比曲線相重合的曲線段對應(yīng)的頻率最小值，作為模擬檢波器可恢復(fù)信號的最低頻率。
6.如權(quán)利要求4所述的方法，其特征在于，在確定第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比之后，所述方法還包括: 根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，按照以下公式計算得到用于對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)的反濾波函數(shù):
7.如權(quán)利要求1至6中任一項所述的方法，其特征在于: 根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率，包括:計算所述模擬檢波器所在的觀測區(qū)中的多個模擬檢波器可恢復(fù)信號的最低頻率；對多個模擬檢波器對應(yīng)的可恢復(fù)信號的最低頻率進(jìn)行擬合，得到應(yīng)用于所述觀測區(qū)中所有模擬檢波器的可恢復(fù)信號的最低頻率； 根據(jù)所述最低頻率，對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)，包括:根據(jù)得到的應(yīng)用于所述觀測區(qū)中所有模擬檢波器的可恢復(fù)信號的最低頻率對該觀測區(qū)中的模擬檢波器輸出的第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。
8.如權(quán)利要求1至6中任一項所述的方法，其特征在于，低頻檢波器和模擬檢波器之間的距離小于預(yù)設(shè)的閾值，包括: 在所述模擬檢波器為一個模擬檢波器的情況下，所述低頻檢波器布置在以該模擬檢波器為圓心，半徑小于預(yù)定的閾值的區(qū)域內(nèi)； 或者，在所述模擬檢波器為多個模擬檢波器的情況下，所述低頻檢波器布置在所述多個模擬檢波器所在位置所構(gòu)成的幾何圖形的中心位置。
9.一種恢復(fù)模擬檢波器地震數(shù)據(jù)低頻信號的裝置，其特征在于，包括: 獲取模塊，用于獲取低頻檢波器輸出的第一地震數(shù)據(jù)和模擬檢波器輸出的第二地震數(shù)據(jù)，其中，低頻檢波器和模擬檢波器之間的距離小于預(yù)設(shè)的閾值，低頻檢波器和模擬檢波器接收到的是同源地震信號； 確定模塊，用于根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，確定模擬檢波器可恢復(fù)信號的最低頻率； 恢復(fù)模塊，用于根據(jù)所述最低頻率，對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。
10.如權(quán)利要求9所述的裝置，其特征在于，確定模塊包括: 第一確定單元，用于確定第一地震數(shù)據(jù)和第二地震數(shù)據(jù)是否具有相同的數(shù)據(jù)格式； 轉(zhuǎn)換單元，用于在確定出數(shù)據(jù)格式不相同的情況下，將第一地震數(shù)據(jù)和第二地震數(shù)據(jù)轉(zhuǎn)換為相同的數(shù)據(jù)格式； 校正單元，用于對相同數(shù)據(jù)格式的第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差進(jìn)行校正； 第二確定單元，用于確定進(jìn)行時差校正后的第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比； 獲取單元，用于根據(jù)功率譜密度比得到模擬檢波器可恢復(fù)信號的最低頻率。
11.如權(quán)利要求10所述的裝置，其特征在于，所述校正單元包括: 計算子單元，用于按照以下公式計算第一地震數(shù)據(jù)和第二地震數(shù)據(jù)之間的時差:
12.如權(quán)利要求9所述的裝置，其特征在于，所述確定模塊包括: 第三確定單元，用于確定第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比； 第四確定單元，用于繪制第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比曲線，將繪制的功率譜密度比曲線與理論的功率譜密度比曲線進(jìn)行對比，確定模擬檢波器可恢復(fù)信號的最低頻率。
13.如權(quán)利要求12所述的裝置，其特征在于，所述第四確定單元具體用于將繪制的功率譜密度比曲線與理論的功率譜密度比曲線相重合的曲線段對應(yīng)的頻率最小值，作為模擬檢波器可恢復(fù)信號的最低頻 率。
14.如權(quán)利要求12所述的 裝置，其特征在于，所述確定模塊還包括: 計算單元，用于根據(jù)第一地震數(shù)據(jù)與第二地震數(shù)據(jù)的功率譜密度比，按照以下公式計算得到用于對所述第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)的反濾波函數(shù):h{m) = I / ^PSDR(m) = I / ^PSDi (m) i PSD2 (ω) 其中，h(?)表示反濾波函數(shù)，PSDR(O)表示第一地震數(shù)據(jù)和第二地震數(shù)據(jù)的功率譜密度比，PSD1(GJ)表示第一地震數(shù)據(jù)的功率譜密度，PSD2(GJ)表示第二地震數(shù)據(jù)的功率譜山/又ο
15.如權(quán)利要求9至14中任一項所述的裝置，其特征在于: 所述確定模塊具體用于計算所述模擬檢波器所在的觀測區(qū)中的多個模擬檢波器可恢復(fù)信號的最低頻率；對多個模擬檢波器對應(yīng)的可恢復(fù)信號的最低頻率進(jìn)行擬合，得到應(yīng)用于所述觀測區(qū)中所有模擬檢波器的可恢復(fù)信號的最低頻率； 所述恢復(fù)模塊具體用于根據(jù)得到的應(yīng)用于所述觀測區(qū)中所有模擬檢波器的可恢復(fù)信號的最低頻率對該觀測區(qū)中的模擬檢波器輸出的第二地震數(shù)據(jù)中的低頻信號進(jìn)行恢復(fù)。
16.如權(quán)利要求9至14中任一項所述的裝置，其特征在于，低頻檢波器和模擬檢波器之間的距離小于預(yù)設(shè)的閾值包括: 在所述模擬檢波器為一個模擬檢波器的情況下，所述低頻檢波器布置在以該模擬檢波器為圓心，半徑小于預(yù)定的閾值的區(qū)域內(nèi)； 或者，在所述模擬檢波器為多個模擬檢波器的情況下，所述低頻檢波器布置在所述多個模擬檢波器所在位置所構(gòu)成的幾何圖形的中心位置。
【文檔編號】G01V1/28GK103499834SQ201310424949
【公開日】2014年1月8日 申請日期:2013年9月17日 優(yōu)先權(quán)日:2013年9月17日
【發(fā)明者】張翊孟, 倪宇東, 蔡錫偉, 李培明, 張慕剛, 王彥娟 申請人:中國石油天然氣集團公司, 中國石油集團東方地球物理勘探有限責(zé)任公司