專利名稱：三維地震數(shù)據(jù)的連續(xù)自適應面波分析的制作方法
技術領域：
本發(fā)明總體上涉及用于三維地震數(shù)據(jù)的連續(xù)面波分析。
背景技術：
地震勘探涉及為尋找碳氫化合物儲藏而勘測地下的地質構造?？睖y典型地涉及在預定的位置布置地震源和地震傳感器。所述源產(chǎn)生地震波，該地震波傳播到地質構造中，沿其路線產(chǎn)生壓力變化和振動。地質構造的彈性特性的變化使地震波反射，折射或分散，改變了它們的傳播方向和其他特性。由所述源發(fā)射出的部分能量到達地震傳感器。一些地震傳感器對壓力變化敏感(水聽器)，其它地震傳感器對粒子運動敏感(例如，地震檢波器)，并且工業(yè)勘測可以只布置一種類型的傳感器或者兩種類型的傳感器都布置。響應于檢測到的地震事件，傳感器產(chǎn)生電信號，以產(chǎn)生地震數(shù)據(jù)。然后，對地震數(shù)據(jù)的分析能夠指示存在或者不存在可能的碳氫化合物儲藏位置。對于基于陸地的地震勘測，地震波包括面波和體波。面波以不輻射到地球內部的方式傳播，平行于地球表面并且與體波相比具有減小的幾何擴散。面波攜帶由地震源輻射于地球表面的大部分能量。面波構成了地震數(shù)據(jù)中的大部分的相干噪聲。以這種方式，面波是以相對低的速率和相對高的幅值為特征的源生成事件，并且面波疊加到有用的信號上。這種相干噪聲 (在陸地地震應用中通常稱為地滾波)可以采用許多不同的波類型的形式，例如具有多種傳播模式的瑞利波(Rayleigh wave)、蘭姆波(Lamb wave)、P導波(P-guided wave)、拉夫波(Love wave)和斯科爾特波(Scholte wave)。面波的傳播特性取決于“近地表(near-surface) ” (地球的淺層部分)的彈性特性，其是造成所獲取的地震數(shù)據(jù)的大部分擾動和劣化的原因。為了設計消減面波噪聲的濾波器，通常重要的是精確地識別面波的特性。另外，了解面波特性可以有利于其他目的，例如確定近地表的局部彈性特性以及估計靜態(tài)校正。

發(fā)明內容
在本發(fā)明的實施方式中，一種技術包括接收在區(qū)域的勘測中所獲取的地震數(shù)據(jù)， 所述區(qū)域與位置關聯(lián)；處理地震數(shù)據(jù)以估計該區(qū)域上的至少一個頻率相關的面波特性；至少部分地基于所估計的頻率相關的面波特性，為每個位置確定頻率相關的數(shù)據(jù)處理幾何配置(data processing geometry)；以及至少部分地基于所確定的數(shù)據(jù)處理幾何配置來處理地震數(shù)據(jù)以得出跨該區(qū)域的面波特性的空間連續(xù)表示。根據(jù)以下附圖、描述和權利要求，本發(fā)明的優(yōu)點和其它特征將變得顯而易見。


圖1和圖3是示出根據(jù)本發(fā)明的實施方式的確定在勘測區(qū)域上的面波特性的空間連續(xù)三維表示的技術的流程圖。圖2是示出根據(jù)本發(fā)明的實施方式的對于給定頻率的示例性局部數(shù)據(jù)處理幾何配置的示意圖。圖4是示出根據(jù)本發(fā)明的實施方式的針對示例性局部數(shù)據(jù)處理幾何配置的偏移及方位角約束的地震勘測陣列的示意圖。圖5A和5B是示出根據(jù)本發(fā)明的實施方式的基于頻率相關的波傳播特性的示例性自適應偏移選擇及變換的示意圖。圖6是根據(jù)本發(fā)明的實施方式的三維條帶(swath)地震獲取幾何配置的示意圖。圖7是根據(jù)本發(fā)明的實施方式的數(shù)據(jù)處理系統(tǒng)的示意圖。
具體實施例方式在此描述了使用從任意地震獲取幾何配置獲取的地震數(shù)據(jù)來構造面波特性(例如速度或衰減)的空間連續(xù)三維(3D)表示的技術。使用在此所公開的技術來確定的面波特性可以用于許多不同的應用中，例如近地表建模、靜態(tài)校正估計、擾動校正、噪聲濾波器設計、淺鉆風險測繪等。如下面所描述的，首先在所勘測區(qū)域的給定位置(例如，在接收器位置)獨立地確定面波特性以產(chǎn)生面波特性的相應的空間連續(xù)的偽截面。然后合并偽截面以產(chǎn)生跨該勘測區(qū)域的面波特性的3D表示。如在此所描述的，為了確定每個位置的面波特性，首先基于考慮了頻率相關的面波特性估計的自適應處理對該位置確定頻率相關的數(shù)據(jù)處理幾何配置。 數(shù)據(jù)處理幾何配置管理著對用于確定面波特性而處理的軌跡的選擇。參考圖1，根據(jù)更具體的實例，技術100可以用來處理由具有任意幾何配置的地震獲取陣列所獲取的地震數(shù)據(jù)以生成給定的面波特性的空間連續(xù)的3D表示。技術100包括首先處理地震數(shù)據(jù)以估計在勘測區(qū)域上的面波特性，如在塊104中所描述的。該步驟可以包括識別一般的面波特性，包括識別波的類型(即，瑞利波、拉夫波、斯科爾特波、P導波等)， 識別波的模式以及識別面波的頻率及速度范圍。應當注意，面波可以是由多于一種類型的基礎面波形成的復合波。然后，技術100包括為勘測區(qū)域的每個位置(例如，每個接收器位置)確定(塊 108)頻率相關的數(shù)據(jù)處理幾何配置。所述數(shù)據(jù)處理幾何配置限定了總的勘測獲取陣列的子陣列。處理來自子陣列的軌跡以確定相關位置的面波特性。數(shù)據(jù)處理幾何配置可以由例如方位角范圍和偏移范圍來限定；如下面所進一步描述的，這些范圍隨頻率變化。同樣如下面所進一步描述的，對于每個位置，基于諸如源-接收器偏移的空間取樣、地震數(shù)據(jù)的幅值譜、在位置附近的面波特性的橫向變化、所希望的空間分辨率、所希望的信噪比(SNR)、幾何配置的尺寸以及在相鄰的幾何配置之間的交疊程度的因素來選擇偏移范圍、方位角范圍、 最大波數(shù)及頻率。如果面波特性在包含給定的子陣列的區(qū)域中是局部均勻的，那么局部面波特性可以表示為源-檢測器偏移的函數(shù)。影響波傳播的區(qū)域限定于該位置的小鄰域內。如下面所進一步描述的，數(shù)據(jù)選擇是自適應的，使得對于每個位置，為每個頻率或頻段選擇不同的軌跡以最小化近場效應并排除遠的偏移(far offset)。作為更具體的實例，圖2示出了對于給定頻率的給定位置的示例性數(shù)據(jù)處理幾何配置10。數(shù)據(jù)處理幾何配置10的邊界由方位角θ所定義，該方位角θ為該位置描繪了被分析的波前。偏移范圍同樣被選擇作為數(shù)據(jù)處理幾何配置10的參數(shù)?？紤]到近場和衰減， 偏移選擇包括了頻率相關的最小及最大偏移。關于這方面，在圖2中示出了示例性的偏移 12。例如，可以將針對該位置的數(shù)據(jù)處理幾何配置10約束在偏移1 到12b之間并且約束在方位角線15和16之間。因為根據(jù)本發(fā)明的某些實施方式數(shù)據(jù)處理幾何配置針對每個具體的接收器/源位置而設計，可以將其他方位角θ和偏移12用于其他數(shù)據(jù)處理幾何配置。仍然參考圖1中的技術100，技術100還包括針對關聯(lián)的數(shù)據(jù)處理幾何配置在每個位置確定(塊11 傳播參數(shù)的局部估計。更具體地，對于與數(shù)據(jù)處理幾何配置關聯(lián)的每組接收器/源，選擇與關聯(lián)的數(shù)據(jù)處理幾何配置相符的射擊集合。對于每個射擊或接收器集合，應用波場變換以將來自時間-源偏移(t-χ)域的數(shù)據(jù)映射到頻率-波數(shù)(f-k)域之內。 于是可以應用高分辨率的變換，而在原始域中不需要不變的間距。所估計頻譜的特性(包括最大波數(shù)和頻譜分辨率)可能隨頻率變化。頻譜被疊合(Stack)，然后提取傳播參數(shù)?？梢詫⑾嗤姆椒☉糜诠餐慕邮掌骷?。然后，對于衰減特性，可以執(zhí)行頻譜幅值的局部回歸。為了估計衰減，工作流的自適應性質還允許數(shù)據(jù)的最優(yōu)選擇。為了得出面波特性的連續(xù)3D空間表示，根據(jù)塊116，技術100還包括執(zhí)行聯(lián)合插值以及在空間上平滑所確定的局部特性。該操作考慮了由于地點橫向變化的規(guī)模(scale)所引起的局部特性的互相關以及它們的空間相關，并且還考慮在相鄰的數(shù)據(jù)處理幾何配置10 之間的交疊。應當注意，可以對于每種波模式和波類型重復塊108、112和116。根據(jù)本發(fā)明的特定實施方式，可以從面波數(shù)據(jù)中提取出若干不同類型的傳播參數(shù)。對于在此所描述的實例，頻散曲線(dispersion curve)(稱為“V(f) ”)和衰減曲線(稱為“D(f)”)是所考慮的傳播參數(shù)。因為面波是頻散的，所以不能夠使用單一傳播速度來描述波速。因此，對于每一個體模式，存在相速度和群速度；以及對波速的完整描述需要指定在整個頻段中的頻散曲線 (相位或群組)。相速度和群速度是相關的，并且因此，一種速度可以由另一種速度所確定。頻散曲線V(f)可以用來以不同的方法(例如頻率-波數(shù)(f-k)過濾、相位匹配等)過濾數(shù)據(jù)，并且還可以用來估計在所考慮的位置的近地表的速度分布。根據(jù)本發(fā)明的某些實施方式，速度(稱為“V”)可以表示如下V = Vph (f)公式 1其中“Vph(f)，，代表相速度。類似地，固有衰減D (f)取決于波數(shù)，并且衰減系數(shù)是頻率相關的。對于每個單一
模式，固有衰減項可以如下寫出e-a 咖=e-K_x 公式 2其中e_a(f)x代表衰減系數(shù)。可以將阻尼比曲線(稱為“D(f)”)看作是局部傳播參數(shù)。由于不同頻率的穿透不同，D(f)阻尼比(以及品質因子Q = O. 5D—1)是頻率相關的。
7
關于數(shù)據(jù)處理幾何配置的設計，設計了用于特定的數(shù)據(jù)處理幾何配置的接收器的每個個體群組的大小以針對所確定的橫向變化以及所希望的橫向分辨率提取傳播參數(shù)。根據(jù)本發(fā)明的某些實施方式，可以設計接收器的最小尺寸及數(shù)量并且將它們用作選擇標準。 例如，在3D幾何配置中，對于接收器的線(receiver line)，可以以在縱測線和橫測線方向上的最大尺寸來提取出接收器的精確群組。作為更具體的實例，圖3示出了流程圖120，其示出了根據(jù)本發(fā)明的實施方式用于得出面波特性的空間連續(xù)的3D表示的技術。根據(jù)技術120，在幾個所選的射擊集合中執(zhí)行模式識別(塊124)。關于這方面，可以選擇跨越相對大的偏移范圍的所選的射擊集合以便分析近地表模式并將它們分類。更具體地，分析所選的射擊集合以檢測多種模式的存在、它們的相速度及衰減。射擊集合可以使用例如頻率-波數(shù)(f_k)變換的波場變換來分析，隨后跟蹤能量的最大值。作為非限制性的實例，可以反演瑞利波的頻散曲線以獲得作為模擬完整的近地表波場的輸入的初步的近地表模型?？梢詫⑺５哪Ｊ脚c所測得的模式進行比較以識別事件。在模式識別之后，技術120包括執(zhí)行(塊128)傳播特性的初步評估以及在所述幾個所選的射擊集合中識別橫向變化。關于這方面，對于每個不同的模式和/或波類型，識別頻率及速度的范圍。確定相速度和群速度，并且估計衰減。然后估計其中在足夠大的SNR下存在待分析的事件的偏移范圍。在所選集合內的橫向變化可以使用相位分析或短間距f_k 變換(作為非限制性的實例)來分析。執(zhí)行對橫向變化的規(guī)模的初步評估，并且將結果與諸如非地震方法及遙感數(shù)據(jù)之類的其他信息進行比較。技術120包括選擇(塊132)要分析的模式和/或波類型以及對于要分析的每種模式/波類型，根據(jù)塊132為每個位置確定數(shù)據(jù)處理幾何配置。幾何配置是在獲得提高的頻譜分辨率與橫向變化所需要的空間分辨率之間的折中。例如，結合圖3來參考圖5A，如在示例性的數(shù)據(jù)處理幾何配置180中所示出的，在特定位置的所選偏移在空間上且隨頻率變化。圖5B示出了在頻率-波數(shù)(f-k)空間內的對應的數(shù)據(jù)處理幾何配置180。如在圖5A和5B所能夠看見的，自適應的選擇及計算允許在其中記錄較短的波長(較大波數(shù))的高頻范圍內的更高的橫向分辨率以及在其中記錄較長的波長(較小波數(shù))的較低的頻率下的更高精度。更具體地，如圖5A所示，近場效應隨頻率的降低而降低，由此導致在較低的頻率下排除了較小的偏移。同樣如圖5A所示，由于在較高頻率下信號劣化較大，因而較遠的偏移在較高的頻率下被排除得更多。根據(jù)本發(fā)明的某些實施方式，在設計數(shù)據(jù)處理幾何配置時可以具體考慮幾何配置的陣列的大小、在相鄰的獲取幾何配置之間的交疊以及最小及最大源偏移。而且，考慮到要處理的數(shù)據(jù)的幾何配置，還確定了最小及最大的方位角。方位角范圍的加寬增加了為每個位置所考慮的射擊或接收器的數(shù)量，由此提高了估計的穩(wěn)健性。通過考慮整個方位角范圍以及比較針對不同方位角集合(azimuths bins)的結果能夠評估方位角的各向異性的存在。但是，同時，沒有準確位于數(shù)據(jù)處理幾何配置的接收器的下方的計算時間和異常的沖擊增加了。因而，根據(jù)塊136，對于每個位置，選擇一組射擊集合或接收器集合；以及對于給定的位置，選擇滿足根據(jù)數(shù)據(jù)處理幾何配置的不同標準的數(shù)據(jù)。例如，在十字排列 (cross-spread)的3D幾何配置中，對于沿著接收器的線或者沿著由多個子線(subline)構成的宏線(macro-line)的每個位置，所考慮的軌跡是在所設計的偏移及方位角范圍之內的所有射擊點。類似地，對于與沿著源的線(source line)的源陣列的共同接收器集合的選擇，使用了偏移和方位角范圍。對于由幾個子線組成的宏線，在給定的位置周圍提取出一組檢測器，并且選擇從一定范圍的偏移和方位角中選擇一組源，如圖4所示。關于這方面，同樣參考圖4，地震勘測陣列150包括接收器子的線160和源的線 154。通過要求將偏移約束在半徑172和半徑176之間來約束對于給定的位置171的源和接收器的位置，并且它們的位置還由定義在線164和168之間的方位角范圍所約束。方位角各向異性能夠通過評估跨不同的方位角集合的特性來評價。因此，塊136的最終結果是得出每個位置的頻率相關的數(shù)據(jù)處理幾何配置。仍然參考圖3，接下來，根據(jù)塊140，對于每個位置，變換在時間-源偏移域內的每個射擊/接收器集合以產(chǎn)生高分辨率的波場變換。變換數(shù)據(jù)被疊合，并且最大的跟蹤 (tracking)被用來得出頻散曲線。作為更具體的實例，根據(jù)本發(fā)明的某些實施方式，對于每個位置，確定間隔不均勻的高分辨率的波場變換，例如在 Schmidt R. 0.的Multiple Emitter Location and Signal Parameter Estimation, IEEE TRANSACTIONS ON ANTENNAS PROPAGATION, Vol.AP-34, 276-280(1986)中所闡述的多重信號分類(MUSIC)算法。對于每個頻率或頻帶，將該變換應用于包括滿足最小及最大偏移的不同標準的軌跡的不同的數(shù)據(jù)集。輸出參數(shù)根據(jù)所預期的最大波數(shù)和所需要的頻譜分辨率(波數(shù)離散化)來設計并且是頻率相關的。波場變換將所獲取的數(shù)據(jù)從該空間變換至波數(shù)域，并且標準的快速傅里葉變換 (FFT)可以用來將數(shù)據(jù)從時間域變換為頻率域，如在塊142中所示出的。在變換之前可以對集合應用減弱(muting)以從面波信號錐的外部去除能量。考慮到數(shù)據(jù)的質量，然后可以通過平均所有個體的貢獻來計算疊合頻譜，并且之后可以進行加權求和。例如通過在個體頻譜中使用相同的拾取算法經(jīng)由拾取分布的統(tǒng)計分析來估計不確定性，在疊合頻譜中跟蹤極大值。然后，根據(jù)技術120，對于每個位置，可以根據(jù)塊144來確定頻散及衰減曲線。以這種方式，為每個位置以及每個集合確定頻譜幅值。在每個頻率下，可以執(zhí)行回歸以計算衰減系數(shù)。個體模式的衰減可以表示如下 ，/) =+⑴公式 3其中“A(x，f) ”代表衰減，“&⑴”代表對于零偏移的衰減。公式3可以如下寫出l0g(^」(x,/))=l0g(4(/))-/}(/).足(/)1 公式 4以偏移的平方根來校正的幅值的對數(shù)是偏移的線性函數(shù)。斜率α (f)可以如下式那樣來確定α (f) = D (f) · K (f)公式 5衰減曲線A(f)被確定，然后與頻散曲線合并以確定阻尼比D(f)。幅值的衰減系數(shù)和殘差被存儲以確定平均衰減曲線，并且殘差的統(tǒng)計量指示擾動。最后，根據(jù)技術120，根據(jù)塊148對頻散曲線和衰減曲線進行插值以建立偽截面， 并且整合偽截面以形成面波特性的空間連續(xù)的3D表示。關于這方面，在每個被考慮的位置，確定了一組頻散及衰減曲線，并且所確定的曲線與數(shù)據(jù)處理幾何配置的中心關聯(lián)。然后合并為不同的位置所獲得的結果以獲得面波特性的連續(xù)平滑的空間表示。由于空間結構和
9相關性并且由于在計算中所使用的相鄰的幾何配置之間的交疊，所獲得的圖像具有作為頻率的函數(shù)的橫向平滑性，并且該特性被用來去除異常值和噪聲部分。作為自適應的工作流的結果，交疊是頻率相關的，并且反映了不同波長的物理橫向分辨率。在速度和衰減之間的互相關可以用來提高規(guī)律性。橫向平滑約束在預期劇烈的橫向變化預期出現(xiàn)的位置是不嚴格的。參考圖6，例如，在條帶獲取陣列200中，垂直的源220和接收器230的線的集合被用來形成偽截面的對應集合。如圖6所示，可以將源220和接收器230的線以如圖6所示的網(wǎng)格210布置。在可以使用在此所描述的系統(tǒng)和技術實現(xiàn)的優(yōu)點當中，面波特性的空間分布可以使用任何幾何配置來獲得。該技術和系統(tǒng)允許以正常的3D幾何配置來處理在每個接收器線和每個源線下方的面波特性的偽截面。所獲得的特性可以用于諸如近地表建模、靜態(tài)校正、用于校準過濾掉相干噪聲的局部濾波器以及產(chǎn)生用于過濾過程的合成噪聲的目的。參考圖7，根據(jù)本發(fā)明的某些實施方式，數(shù)據(jù)處理系統(tǒng)520可以包括處理器550，該處理器550被構造為執(zhí)行至少一個程序544 (存儲于存儲器540內)，以便處理地震數(shù)據(jù)，從而執(zhí)行在此所描述的一種或多種技術。處理器550可以與通信接口 560耦接以便接收在地震勘測中所獲取的地震數(shù)據(jù)。除了存儲用于程序554的指令外，存儲器540還可以存儲涉及在此所公開的技術的初步的、中間的及最終的數(shù)據(jù)集。其中，數(shù)據(jù)處理系統(tǒng)520可以包括顯示接口 570和顯示器574，以便顯示如在此所描述的那樣生成的各種曲線和模型。雖然本發(fā)明已經(jīng)針對數(shù)量有限的實施方式進行了描述，但是受益于本公開內容的本領域技術人員將認識到根據(jù)其進行的多種修改和變化。希望所附權利要求覆蓋落在本發(fā)明真正主旨和范圍內的所有這類修改和變化。
權利要求
1.一種方法，包括接收在區(qū)域的勘測中所獲取的地震數(shù)據(jù)，所述區(qū)域與位置關聯(lián)； 處理所述地震數(shù)據(jù)以估計在所述區(qū)域上的至少一個頻率相關的面波特性； 至少部分地基于所述至少一個所估計的頻率相關的面波特性，為每個位置確定頻率相關的數(shù)據(jù)處理幾何配置；以及至少部分地基于所確定的數(shù)據(jù)處理幾何配置來處理所述地震數(shù)據(jù)以得出跨所述區(qū)域的面波特性的空間連續(xù)表示。
2.根據(jù)權利要求1所述的方法，其中所述至少一個所估計的頻率相關的面波特性包括下列各項中的至少一個近場效應、遠偏移衰減、面波特性的頻譜估計以及面波特性的橫向變化。
3.根據(jù)權利要求1所述的方法，進一步包括 處理所述地震數(shù)據(jù)以識別面波模式；以及執(zhí)行以下過程處理所述地震數(shù)據(jù)以估計所述至少一個頻率相關的面波特性，為每個位置確定所述頻率相關的數(shù)據(jù)處理幾何配置以及基于每個所識別的模式的所確定的數(shù)據(jù)處理幾何配置來處理所述地震數(shù)據(jù)。
4.根據(jù)權利要求1所述的方法，進一步包括處理所述地震數(shù)據(jù)以確定與足夠的信噪比關聯(lián)的且沒有近場效應的頻率相關的偏移范圍；以及進一步使為每個位置確定所述頻率相關的數(shù)據(jù)處理幾何配置的過程基于所確定的頻率范圍。
5.根據(jù)權利要求1所述的方法，進一步包括進一步使為每個位置確定所述頻率相關的數(shù)據(jù)處理幾何配置的過程基于與所述數(shù)據(jù)處理幾何配置關聯(lián)的接收器陣列的大小。
6.根據(jù)權利要求1所述的方法，進一步包括進一步使為每個位置確定所述頻率相關的數(shù)據(jù)處理幾何配置的過程基于同所述數(shù)據(jù)處理幾何配置關聯(lián)的接收器陣列與同至少一個其他數(shù)據(jù)處理幾何配置關聯(lián)的至少一個其他接收器陣列的交疊。
7.根據(jù)權利要求1所述的方法，其中確定所述頻率相關的數(shù)據(jù)處理幾何配置的過程包括確定針對所述幾何配置的方位角范圍。
8.根據(jù)權利要求1所述的方法，其中確定所述頻率相關的數(shù)據(jù)處理幾何配置的過程包括評估方位角各向異性的存在；以及確定針對所述幾何配置的偏移范圍。
9.根據(jù)權利要求1所述的方法，其中至少部分地基于所確定的數(shù)據(jù)處理幾何配置來處理所述地震數(shù)據(jù)的過程包括處理所述地震數(shù)據(jù)以確定相速度、群速度及衰減中的至少一個；以及至少部分地基于波數(shù)的頻率相關性來優(yōu)化所述處理。
10.根據(jù)權利要求1所述的方法，其中至少部分地基于所確定的數(shù)據(jù)處理幾何配置來處理所述地震數(shù)據(jù)的過程包括為每個位置，基于關聯(lián)的數(shù)據(jù)處理幾何配置來選擇集合以及處理所選擇的集合。
11.根據(jù)權利要求1所述的方法，其中至少部分地基于所確定的數(shù)據(jù)處理幾何配置來處理所述地震數(shù)據(jù)的過程包括為每個位置，確定在與所述位置關聯(lián)的子區(qū)域中的面波特性的空間連續(xù)表示；以及合并在所述子區(qū)域中的面波特性的空間連續(xù)表示，以生成跨所述區(qū)域的所述面波特性的空間連續(xù)表示。
12.—種系統(tǒng)，包括接收在區(qū)域的勘測中所獲取的地震數(shù)據(jù)的接口，所述區(qū)域與位置關聯(lián)；以及處理器，用于處理所述地震數(shù)據(jù)以估計在所述區(qū)域上的至少一個頻率相關的面波特性； 至少部分地基于所述至少一個所估計的頻率相關的面波特性，為每個位置確定頻率相關的數(shù)據(jù)處理幾何配置；以及至少部分地基于所確定的數(shù)據(jù)處理幾何配置來處理所述地震數(shù)據(jù)以得出跨所述區(qū)域的所述面波特性的空間連續(xù)表示。
13.根據(jù)權利要求12所述的系統(tǒng)，其中所述至少一個所估計的頻率相關的面波特性包括下列各項中的至少一個近場效應、遠偏移衰減、面波特性的頻譜估計以及面波特性的橫向變化。
14.根據(jù)權利要求12所述的系統(tǒng)，其中所述處理器適合于 處理所述地震數(shù)據(jù)以識別面波模式；以及基于所識別的模式來估計所述至少一個頻率相關的面波特性。
15.根據(jù)權利要求12所述的系統(tǒng)，其中所述處理器適合于進一步使為每個位置確定所述頻率相關的數(shù)據(jù)處理幾何配置的過程基于與所述數(shù)據(jù)處理幾何配置關聯(lián)的接收器陣列的大小。
16.根據(jù)權利要求12所述的系統(tǒng)，其中所述處理器適合于進一步使為每個位置確定所述頻率相關的數(shù)據(jù)處理幾何配置的過程基于同所述數(shù)據(jù)處理幾何配置關聯(lián)的接收器陣列與同至少一個其他數(shù)據(jù)處理幾何配置關聯(lián)的至少一個其他接收器陣列的交疊。
17.根據(jù)權利要求12所述的系統(tǒng)，其中所述處理器適合于 確定針對所述幾何配置的方位角范圍。
18.根據(jù)權利要求12所述的系統(tǒng)，其中所述處理器適合于 確定針對所述幾何配置的偏移范圍。
19.根據(jù)權利要求12所述的系統(tǒng)，其中所述處理器適合于處理所述地震數(shù)據(jù)以確定相速度、群速度及衰減中的至少一個。
20.根據(jù)權利要求12所述的系統(tǒng)，其中所述處理器適合于為每個位置，基于關聯(lián)的頻率相關的數(shù)據(jù)處理幾何配置來選擇集合以及處理所選擇的集合。
21.根據(jù)權利要求12所述的系統(tǒng)，其中所述處理器適合于為每個位置，確定在與所述位置關聯(lián)的子區(qū)域中的面波特性的空間連續(xù)表示；以及合并在所述子區(qū)域中的面波特性的空間連續(xù)表示以生成跨所述區(qū)域的面波特性的空間連續(xù)表示。
22. 一種包括存儲指令的計算機可讀的存儲介質的物品，所述指令在被計算機執(zhí)行時促使所述計算機執(zhí)行如權利要求1-11所述的方法。
全文摘要
一種技術，包括接收在區(qū)域的勘測中所獲取的地震數(shù)據(jù)，所述區(qū)域與位置關聯(lián)；處理地震數(shù)據(jù)以估計在該區(qū)域上的至少一個頻率相關的面波特性；至少部分地基于所估計的頻率相關的面波特性，為每個位置確定頻率相關的數(shù)據(jù)處理幾何配置；以及至少部分地基于所確定的數(shù)據(jù)處理幾何配置來處理地震數(shù)據(jù)以得出跨該區(qū)域的面波特性的空間連續(xù)表示。
文檔編號G01V1/28GK102227650SQ200980147176
公開日2011年10月26日 申請日期2009年11月20日 優(yōu)先權日2008年11月26日
發(fā)明者A·格魯施申科, C·L·斯特羅比亞 申請人:格庫技術有限公司