一種三維地震觀測系統多頻率偏移噪聲獲取方法及裝置製造方法
2023-06-04 12:54:01
一種三維地震觀測系統多頻率偏移噪聲獲取方法及裝置製造方法
【專利摘要】本發明實施例提供一種三維地震觀測系統多頻率偏移噪聲獲取方法及裝置，所述方法包括：獲取三維地震觀測系統數據；根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲；對所述三維地震觀測系統的單頻率偏移噪聲進行插值，獲取所述三維地震觀測系統的單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲；對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲。本發明實施例僅僅進行一次單頻率偏移噪聲運算，在保證計算精度的同時不增加算法的整體計算量，從而實現三維地震觀測系統多頻率偏移噪聲的快速獲取。
【專利說明】一種三維地震觀測系統多頻率偏移噪聲獲取方法及裝置
【技術領域】
[0001]本發明涉及三維地震觀測系統分析評價技術，尤其涉及一種三維地震觀測系統多頻率偏移噪聲獲取方法及裝置。
【背景技術】
[0002]三維地震勘探是石油與天然氣勘探的主要工具。三維地震勘探的主要過程包括:
(I)數據資料採集。在陸上地震勘探數據採集現場一般進行下述三個工作:地震觀測系統設計、根據設計方案在野外布設震源和檢波器，地震波的激發和接收。首先在室內進行地震觀測系統設計，以確定震源點和檢波點的最佳擺放位置。然後根據設計方案，在野外布設震源和檢波器，陸上地震勘探的震源點一般採用炸藥震源，並沿地震測線等間距布置多個檢波器來接收地震波信號，現代地震勘探中檢波器的數量多在1000或10000個。震源點在爆炸後產生出地震波，地震波遇巖層界面反射回來被檢波器接收並傳到儀器車，儀器車將檢波器傳來的信號記錄下來，這就獲得了用以研究地下油氣埋藏情況的地震記錄。(2)地震數據處理。地震數據處理是把第一步驟採集到的地震數據資料輸入專用電子計算機，按不同要求用一系列功能不同的程序進行處理運算，把數據進行歸類編排，突出有效的，除去無效和幹擾的，最後把經過各種處理的數據進行疊加和偏移，最終得到二維或三維地震數據體文件。(3)資料解釋。資料解釋是把經過處理的地震信息變成地質成果的過程，包括運用波動理論和地質知識，綜合地質、鑽井、測井等各項資料，作出構造解釋、地層解釋、巖性和烴類檢測解釋及綜合解釋，繪出有關成果圖件，對勘探區域作出含油氣評價，提出鑽探井位置坐寸ο
[0003]在整個地震勘探的處理流程中，地震疊前偏移方法是最重要的地震資料處理流程之一，其效果直接決定了地震資料處理的成敗甚至整個地震勘探的成敗。前面的介紹曾提至IJ，三維地震觀測系統設計是整個地震勘探的第一個環節，因此，必須從三維地震觀測系統設計階段考慮地震疊前偏移的效果(即偏移噪聲)，以選擇合適的地震觀測系統設計方案來完成地震勘探任務。目前直接採用地震疊前偏移理論計算三維地震觀測系統偏移噪聲的方法主要有兩種:(1)基於Beylkin理論，結合檢波點與震源點的分析方法；(2)Berkhout共聚焦分析方法。目前這兩種方法均採用地震波射線追蹤或地震波KirchhofT算子實現，因此僅適用於均勻介質或層狀介質情況，無法適用於複雜的地下介質條件，且計算精度較低。
【發明內容】
[0004]本發明實施例提供一種三維地震觀測系統多頻率偏移噪聲獲取方法及裝置，通過僅僅進行一次單頻率偏移噪聲運算，利用快速偏移噪聲插值算法將其推廣至多頻率條件下，在保證計算精度的同時不增加算法的整體計算量，從而實現三維地震觀測系統多頻率偏移噪聲的快速獲取。
[0005]一方面，本發明實施例提供了一種三維地震觀測系統多頻率偏移噪聲獲取方法，所述三維地震觀測系統多頻率偏移噪聲獲取方法包括:獲取三維地震觀測系統數據；根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲；對所述三維地震觀測系統的單頻率偏移噪聲進行插值，獲取所述三維地震觀測系統的單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲；對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲。
[0006]另一方面，本發明實施例提供了一種三維地震觀測系統多頻率偏移噪聲獲取裝置，所述三維地震觀測系統多頻率偏移噪聲獲取裝置包括:數據獲取單元，用於獲取三維地震觀測系統數據；單頻率偏移噪聲單元，用於根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲；插值單元，用於對所述三維地震觀測系統的單頻率偏移噪聲進行插值，獲取所述三維地震觀測系統的單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲；多頻率偏移噪聲單元，用於對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲。
[0007]上述技術方案具有如下有益效果:因為採用所述三維地震觀測系統多頻率偏移噪聲獲取方法包括:獲取三維地震觀測系統數據；根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲；對所述三維地震觀測系統的單頻率偏移噪聲進行插值，獲取所述三維地震觀測系統的單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲；對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲的技術手段，所以達到了如下的技術效果:克服了現有三維地震觀測系統多頻率偏移噪聲獲取方法運算速度較慢的缺陷，提出了一種可以適用於簡單或複雜模型的三維地震觀測系統多頻率偏移噪聲的快速獲取方法。可將僅僅進行一次單頻率偏移噪聲運算，利用快速偏移噪聲插值算法將其推廣至多頻率條件下，在保證計算精度的同時不增加算法的整體計算量，從而實現三維地震觀測系統多頻率偏移噪聲的快速獲取。本發明實施例大幅度提升多頻率偏移噪聲計算的速度，有力地促進了偏移噪聲分析技術在實際地震勘探中的應用。
【專利附圖】
【附圖說明】
[0008]為了更清楚地說明本發明實施例或現有技術中的技術方案，下面將對實施例或現有技術描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發明的一些實施例，對於本領域普通技術人員來講，在不付出創造性勞動的前提下，還可以根據這些附圖獲得其他的附圖。
[0009]圖1為本發明實施例一種三維地震觀測系統多頻率偏移噪聲獲取方法流程圖；
[0010]圖2為本發明實施例一種三維地震觀測系統多頻率偏移噪聲獲取裝置結構示意圖；
[0011]圖3為本發明應用實例多頻率偏移噪聲計算過程的示意圖；
[0012]圖4為本發明應用實例某三維地震觀測系統多頻率偏移噪聲分析結果示意圖。
【具體實施方式】
[0013]下面將結合本發明實施例中的附圖，對本發明實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本發明一部分實施例，而不是全部的實施例。基於本發明中的實施例，本領域普通技術人員在沒有做出創造性勞動前提下所獲得的所有其他實施例，都屬於本發明保護的範圍。
[0014]如圖1所示，為本發明實施例一種三維地震觀測系統多頻率偏移噪聲獲取方法流程圖，所述三維地震觀測系統多頻率偏移噪聲獲取方法包括:
[0015]101、獲取三維地震觀測系統數據；
[0016]優選的，所述三維地震觀測系統數據包括:震源點的三維空間坐標、檢波點三維空間坐標、震源點與檢波點之間的對應關係；其中，所述震源點的三維空間坐標包括:高程；所述檢波點三維空間坐標包括:高程。
[0017]102、根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲；
[0018]優選的，所述根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲，包括:
[0019]對於均勻介質，利用如下公式計算單頻率偏移噪聲:
[0020](/): 士 If cxP [z (/f'Δ?+人'Av)] D'dk' dK
[0021]其中，i為虛數單位，kjPky分別為水平X方向和水平y方向波數，AdP Ay分別為檢波點到地下目標反射點X方向和y方向的水平距離分量，D'為檢波點矩陣D的波數域變換，Bd(f)為單頻率偏移噪聲。
[0022]優選的，所述根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲，包括:
[0023]對於非均勻介質，利用如下公式計算單頻率偏移噪聲:
[0024]Bd (f) ^ f2 I I exp [2 n if (ρχ Δ x+py Δ y) ] D;/ dpxdpy,
[0025]其中，i為虛數單位，f為頻率，Px和Py分別為水平X方向和水平y方向慢度，Δχ和Ay分別為檢波點到地下目標反射點X方向和y方向的水平距離分量，D"為檢波點矩陣D的Radon域變換，Bd(f)為單頻率偏移噪聲。
[0026]103、對所述三維地震觀測系統的單頻率偏移噪聲進行插值，獲取所述三維地震觀測系統的單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲；
[0027]104、對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲。
[0028]優選的，所述對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲，包括:
[0029]利用如下公式計算所述三維地震觀測系統的多頻率偏移噪聲:
[0030]J"㈧賴5〔令，令』胃)，
[0031]其中，B為多頻率偏移噪聲，N是頻率採樣點數，Af為頻率採樣間隔，w(f)為雷克
, 、2 (/Δ/)" 「-(/Δ/『)2
子波頻率域表達式W(.M/) = "-JT-exP —γ_~ fp為地震波主頻。
[0032]對應於上述方法實施例，如圖2所示，為本發明實施例一種三維地震觀測系統多頻率偏移噪聲獲取裝置結構示意圖，所述三維地震觀測系統多頻率偏移噪聲獲取裝置包括:
[0033]數據獲取單元21，用於獲取三維地震觀測系統數據；
[0034]單頻率偏移噪聲單元22，用於根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲；
[0035]插值單元23，用於對所述三維地震觀測系統的單頻率偏移噪聲進行插值，獲取所述三維地震觀測系統的單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲；
[0036]多頻率偏移噪聲單元24，用於對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲。
[0037]優選的，所述三維地震觀測系統數據包括:震源點的三維空間坐標、檢波點三維空間坐標、震源點與檢波點之間的對應關係；其中，所述震源點的三維空間坐標包括:高程；所述檢波點三維空間坐標包括:高程。
[0038]優選的，所述單頻率偏移噪聲單元22對於均勻介質，利用如下公式計算單頻率偏移噪聲:
[0039]^/,(/)? ^rjjcxp / (/γλ Δν + /f, Δν j] D'dk, dkv
[0040]其中，i為虛數單位，kx和ky分別為水平X方向和水平I方向波數，Λ X和Λ y分別為檢波點到地下目標反射點X方向和y方向的水平距離分量，D'為檢波點矩陣D的波數域變換，Bd(f)為單頻率偏移噪聲。
[0041]優選的，所述單頻率偏移噪聲單元22對於非均勻介質，利用如下公式計算單頻率偏移噪聲:
[0042]Bd (f) ^ f2 I I exp [2 n if (ρχ Δ x+py Δ y) ] D;/ dpxdpy,
[0043]其中，i為虛數單位，f為頻率，Px和Py分別為水平X方向和水平y方向慢度，Δχ和Ay分別為檢波點到地下目標反射點X方向和y方向的水平距離分量，D"為檢波點矩陣D的Radon域變換，Bd(f)為單頻率偏移噪聲。
[0044]優選的，所述多頻率偏移噪聲單元24利用如下公式計算所述三維地震觀測系統的多頻率偏移噪聲:
[0045]J"秦印技坪爾、，
[0046]其中，B為多頻率偏移噪聲，N是頻率採樣點數，Af為頻率採樣間隔，w (f)為雷克
2 (/'Δ/ ν 「-(/Δ/.):
子波頻率域表達式叫./#) = ;^^^… f2 &為地震波主頻。
PL ^ J3
[0047]本發明實施例上述技術方案具有如下有益效果:因為採用所述三維地震觀測系統多頻率偏移噪聲獲取方法包括:獲取三維地震觀測系統數據；根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲；對所述三維地震觀測系統的單頻率偏移噪聲進行插值，獲取所述三維地震觀測系統的單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲；對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲的技術手段，所以達到了如下的技術效果:克服了現有三維地震觀測系統多頻率偏移噪聲獲取方法運算速度較慢的缺陷，提出了一種可以適用於簡單或複雜模型的三維地震觀測系統多頻率偏移噪聲的快速獲取方法。可將僅僅進行一次單頻率偏移噪聲運算，利用快速偏移噪聲插值算法將其推廣至多頻率條件下，在保證計算精度的同時不增加算法的整體計算量，從而實現三維地震觀測系統多頻率偏移噪聲的快速獲取。本發明實施例大幅度提升多頻率偏移噪聲計算的速度，有力地促進了偏移噪聲分析技術在實際地震勘探中的應用。
[0048]以下結合應用實例對本發明實施例上述技術方案進行詳細說明:
[0049]通常偏移噪聲分析均在單頻率下進行，然而實際的地震波頻率總是會有一定的帶寬，即總是多頻率的。為了解決上述矛盾，必須將單頻率偏移噪聲分析擴展至多頻率偏移噪聲分析。如圖3所示，為本發明應用實例多頻率偏移噪聲計算過程的示意圖。多頻率偏移噪聲結果可用下式表示:
[0050]
N 『廠'⑴
[0051]方程(I)表示所有頻率波場的能量和，其中，B為多頻率偏移噪聲，N是頻卒採樣點數，Af為頻率採樣間隔，B (j Af)為單頻率偏移噪聲，w (f)為雷克子波頻率域表達式
, 、2 (/AfY 「-(/Δ/.)2Η』(./Δ/ ) = ^ CXP '/■- fP 為地震波主頻。
[0052]對於均勻介質，偏移噪聲公式可表示為:
[0053]
^/；(/)?^rjjcxp z(/v.Δν + /c,Δν)]DWZfic//f,
兀(2)
[0054]其中，i為虛數單位，kx和ky分別為水平X方向和水平I方向波數，Λ X和Λ y分別為檢波點到地下目標反射點X方向和y方向的水平距離分量，D'為檢波點矩陣D的波數域變換。
[0055]對於非均勻介質，偏移噪聲公式可表示為:
[0056]Bd (f) ^ f2 I I exp [2 n if (ρχ Δ x+py Δ y) ] D;/ dpxdpy (3)
[0057]其中，i為虛數單位，f為頻率，Px和Py分別為水平X方向和水平y方向慢度(和速度相對，是速度的倒數)，Λ X和Ay分別為檢波點到地下目標反射點X方向和y方向的水平距離分量，D"為檢波點矩陣D的Radon域變換。
[0058]綜合方程⑵和(3)，可得:
[0059]β (紅 Αν,./'Δ/'./Αν,./Δν, JAf)
■)(4)
[0060]代入方程(I)，即可得到多頻率偏移噪聲快速獲取算法:
[0061]B' Affj J2 w(jAf)B(NAf⑶
\ MV TV TVJ
[0062]方程(5)表明，進行一次單頻率偏移噪聲分析與插值運算，即可近似得到三維地震觀測系統多頻率偏移噪聲結果，如圖4所示，為本發明應用實例某三維地震觀測系統多頻率偏移噪聲分析結果示意圖。
[0063]參照圖1，本發明應用實例包括以下步驟:
[0064]I)導入複雜地表下某地震觀測系統設計方案
[0065]觀測系統設計方案包括檢波點的三維空間坐標(含高程)、震源點三維空間坐標(含高程)，和震源點與檢波點之間的對應關係。一個地震觀測系統設計方案包含大量的震源點和檢波點的三維坐標信息。這裡的震源點指地震勘探普遍採用的人工炸藥震源。檢波點指地震勘探中在地表布設的檢波器，一般至少有上千個，用於接收震源發出經地下反射後回到地表的地震波。
[0066]2)計算三維地震觀測系統單頻率偏移噪聲
[0067]利用方程(2)或(3)計算三維地震觀測系統單頻率偏移噪聲。
[0068]3)插值得到其它頻率偏移噪聲結果
[0069]利用方程(4)插值得到其它頻率下偏移噪聲結果。
[0070]4)對所有頻率成分求和得到多頻率偏移噪聲結果
[0071]利用方程(5)得到多頻率偏移噪聲，實現僅利用一次單頻率偏移噪聲分析與插值運算，即可近似得到多頻率下偏移噪聲分析結果。
[0072]參照圖4所示本發明應用實例某三維地震觀測系統多頻率偏移噪聲分析結果，與傳統多頻率偏移噪聲分析不同，該技術僅需要進行一次單頻率偏移噪聲運算，進而利用快速插值算法得到多頻率偏移噪聲結果，在保證計算精度的同時不增加算法的整體計算量，從而實現三維地震觀測系統多頻率偏移噪聲的快速獲取方法。該發明大幅度提升多頻率偏移噪聲計算的速度，有力的促進了偏移噪聲分析技術在實際地震勘探中的應用。
[0073]相對於現有技術，本發明應用實例技術方案不需要對每個頻率單獨進行一次偏移噪聲分析運算，再對所有的頻率成分進行求和得到多頻率下結果，效率較高，適宜在石油勘探領域的推廣和應用。
[0074]本領域技術人員還可以了解到本發明實施例列出的各種說明性邏輯塊(illustrative logical block),單元,和步驟可以通過電子硬體、電腦軟體,或兩者的結合進行實現。為清楚展示硬體和軟體的可替換性(interchangeability),上述的各種說明性部件(illustrative components),單元和步驟已經通用地描述了它們的功能。這樣的功能是通過硬體還是軟體來實現取決於特定的應用和整個系統的設計要求。本領域技術人員可以對於每種特定的應用，可以使用各種方法實現所述的功能，但這種實現不應被理解為超出本發明實施例保護的範圍。
[0075]本發明實施例中所描述的各種說明性的邏輯塊，或單元都可以通過通用處理器，數位訊號處理器，專用集成電路(ASIC)，現場可編程門陣列或其它可編程邏輯裝置，離散門或電晶體邏輯，離散硬體部件，或上述任何組合的設計來實現或操作所描述的功能。通用處理器可以為微處理器，可選地，該通用處理器也可以為任何傳統的處理器、控制器、微控制器或狀態機。處理器也可以通過計算裝置的組合來實現，例如數位訊號處理器和微處理器，多個微處理器，一個或多個微處理器聯合一個數位訊號處理器核，或任何其它類似的配置來實現。
[0076]本發明實施例中所描述的方法或算法的步驟可以直接嵌入硬體、處理器執行的軟體模塊、或者這兩者的結合。軟體模塊可以存儲於RAM存儲器、快閃記憶體、ROM存儲器、EPROM存儲器、EEPROM存儲器、寄存器、硬碟、可移動磁碟、⑶-ROM或本領域中其它任意形式的存儲媒介中。示例性地，存儲媒介可以與處理器連接，以使得處理器可以從存儲媒介中讀取信息，並可以向存儲媒介存寫信息。可選地，存儲媒介還可以集成到處理器中。處理器和存儲媒介可以設置於ASIC中，ASIC可以設置於用戶終端中。可選地，處理器和存儲媒介也可以設置於用戶終端中的不同的部件中。
[0077]在一個或多個示例性的設計中，本發明實施例所描述的上述功能可以在硬體、軟體、固件或這三者的任意組合來實現。如果在軟體中實現，這些功能可以存儲與電腦可讀的媒介上，或以一個或多個指令或代碼形式傳輸於電腦可讀的媒介上。電腦可讀媒介包括電腦存儲媒介和便於使得讓電腦程式從一個地方轉移到其它地方的通信媒介。存儲媒介可以是任何通用或特殊電腦可以接入訪問的可用媒體。例如，這樣的電腦可讀媒體可以包括但不限於RAM、ROM、EEPROM、CD-ROM或其它光碟存儲、磁碟存儲或其它磁性存儲裝置，或其它任何可以用於承載或存儲以指令或數據結構和其它可被通用或特殊電腦、或通用或特殊處理器讀取形式的程序代碼的媒介。此外，任何連接都可以被適當地定義為電腦可讀媒介，例如，如果軟體是從一個網站站點、伺服器或其它遠程資源通過一個同軸電纜、光纖電纜、雙絞線、數字用戶線(DSL)或以例如紅外、無線和微波等無線方式傳輸的也被包含在所定義的電腦可讀媒介中。所述的碟片(disk)和磁碟(disc)包括壓縮磁碟、鐳射盤、光碟、DVD、軟盤和藍光光碟，磁碟通常以磁性複製數據，而碟片通常以雷射進行光學複製數據。上述的組合也可以包含在電腦可讀媒介中。
[0078]以上所述的【具體實施方式】，對本發明的目的、技術方案和有益效果進行了進一步詳細說明，所應理解的是，以上所述僅為本發明的【具體實施方式】而已，並不用於限定本發明的保護範圍，凡在本發明的精神和原則之內，所做的任何修改、等同替換、改進等，均應包含在本發明的保護範圍之內。
【權利要求】
1.一種三維地震觀測系統多頻率偏移噪聲獲取方法，其特徵在於，所述三維地震觀測系統多頻率偏移噪聲獲取方法包括: 獲取三維地震觀測系統數據； 根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲； 對所述三維地震觀測系統的單頻率偏移噪聲進行插值，獲取所述三維地震觀測系統的單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲； 對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲。
2.如權利要求1所述三維地震觀測系統多頻率偏移噪聲獲取方法，其特徵在於，所述三維地震觀測系統數據包括:震源點的三維空間坐標、檢波點三維空間坐標、震源點與檢波點之間的對應關係；其中，所述震源點的三維空間坐標包括:高程；所述檢波點三維空間坐標包括:聞程。
3.如權利要求1所述三維地震觀測系統多頻率偏移噪聲獲取方法，其特徵在於，所述根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲，包括: 對於均勻介質，利用如下公式計算單頻率偏移噪聲:
Bl>CXP [Ζ.(/f'Δν + k' Δ>.)_| D'dk, dky 其中，i為虛數單位，kjPky分別為水平X方向和水平y方向波數，AdP Ay分別為檢波點到地下目標反射點X方向和y方向的水平距離分量，D'為檢波點矩陣D的波數域變換，BD(f)為單頻率偏移噪聲。
4.如權利要求1所述三維地震觀測系統多頻率偏移噪聲獲取方法，其特徵在於，所述根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲，包括: 對於非均勻介質，利用如下公式計算單頻率偏移噪聲:
Bd (f) ^ f2 f f exp [2 if (ρχ Δ x+py Δ y) ] D;/ dpxdpy, 其中，i為虛數單位，f為頻率，Pj^P Py分別為水平X方向和水平y方向慢度，Λχ和Ay分別為檢波點到地下目標反射點X方向和y方向的水平距離分量，D"為檢波點矩陣D的Radon域變換，Bd(f)為單頻率偏移噪聲。
5.如權利要求3或4所述三維地震觀測系統多頻率偏移噪聲獲取方法，其特徵在於，所述對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲，包括: 利用如下公式計算所述三維地震觀測系統的多頻率偏移噪聲: 廠I — 5? Δ/Χ J2 w(jA/)Bi^,^,m/)，
\ j=\ |_\ ^ ^ J 其中，B為多頻率偏移噪聲，N是頻率釆樣點數，為頻率釆樣間隔，w(f)為雷克子波
, 、2.(/Aff 「-?/Δ,)]頻率域表達式w(./_A/)=lh^exp ~ &為地震波主頻。
6.一種三維地震觀測系統多頻率偏移噪聲獲取裝置，其特徵在於，所述三維地震觀測系統多頻率偏移噪聲獲取裝置包括: 數據獲取單元，用於獲取三維地震觀測系統數據； 單頻率偏移噪聲單元，用於根據所述三維地震觀測系統數據，計算所述三維地震觀測系統的單頻率偏移噪聲； 插值單元，用於對所述三維地震觀測系統的單頻率偏移噪聲進行插值，獲取所述三維地震觀測系統的單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲； 多頻率偏移噪聲單元，用於對所述三維地震觀測系統的單頻率偏移噪聲和所述單頻率偏移噪聲之外的一種或多種頻率下偏移噪聲進行求和，獲取所述三維地震觀測系統的多頻率偏移噪聲。
7.如權利要求6所述三維地震觀測系統多頻率偏移噪聲獲取裝置，其特徵在於，所述三維地震觀測系統數據包括:震源點的三維空間坐標、檢波點三維空間坐標、震源點與檢波點之間的對應關係；其中，所述震源點的三維空間坐標包括:高程；所述檢波點三維空間坐標包括:聞程。
8.如權利要求6所述三維地震觀測系統多頻率偏移噪聲獲取裝置，其特徵在於，所述單頻率偏移噪聲單元對於均勻介質，利用如下公式計算單頻率偏移噪聲:
Bn U.) ?-y-v --exp「i(k' Ax + k' Δν)1 D'dk' dk'
4π~L v.-， 其中，i為虛數單位，kjPky分別為水平X方向和水平y方向波數，AdP Ay分別為檢波點到地下目標反射點X方向和y方向的水平距離分量，D'為檢波點矩陣D的波數域變換，BD(f)為單頻率偏移噪聲。
9.如權利要求6所述三維地震觀測系統多頻率偏移噪聲獲取裝置，其特徵在於，所述單頻率偏移噪聲單元對於非均勻介質，利用如下公式計算單頻率偏移噪聲:
Bd (f) ^ f2 f f exp [2 if (ρχ Δ x+py Δ y) ] D;/ dpxdpy, 其中，i為虛數單位，f為頻率，Pj^P Py分別為水平X方向和水平y方向慢度，Λχ和Ay分別為檢波點到地下目標反射點X方向和y方向的水平距離分量，D"為檢波點矩陣D的Radon域變換，Bd(f)為單頻率偏移噪聲。
10.如權利要求8或9所述三維地震觀測系統多頻率偏移噪聲獲取裝置，其特徵在於，所述多頻率偏移噪聲單元利用如下公式計算所述三維地震觀測系統的多頻率偏移噪聲:
警，譬，規/)[， 其中，B為多頻率偏移噪聲，N是頻率採樣點數，Af為頻率採樣間隔，w (f)為雷克子波
, 、2 (jAff 「-?/Δ/)2頻率域表達式Η』(./_Δ/產fi fp為地震波主頻。
7
【文檔編號】G01V1/28GK104237939SQ201410466541
【公開日】2014年12月24日 申請日期:2014年9月12日 優先權日:2014年9月12日
【發明者】符力耘, 魏偉, 管西竹 申請人:中國科學院地質與地球物理研究所