專利名稱：用于地震道分析的系統和方法
技術領域：
本發明一般涉及地質數據的處理，尤其涉及用于減小采樣數據中的混疊影響的系 統。
背景技術：
地震勘測往往產生稀疏和/或非均勻采樣的數據。這可能源自非均勻的震源和接 收器位置、系統和導航不完善、地物障礙以及大型儀器陣列的過高費用。另一方面，許多數 據處理和可視化系統和方法需要密集和規則采樣的數據作為輸入。因此，常常需要變換數 據，使得數據具有適當的采樣基礎。
發明內容
本發明的實施例的一些方面提供了一種處理代表與地下區域有關的信息的采樣 地震道的方法，包括使用離散傅立葉變換來估計地震道的傅立葉系數；通過將相干準則 應用于估計的傅立葉系數來選擇估計的傅立葉系數，以減小所述傅立葉系數所代表的混疊 能量的量；在k_空間中重新正交化所述傅立葉系數；以及使用重新正交化的傅立葉系數來 產生地下區域的地震圖像。本發明的實施例的一些方面包括一種成像地下特征的方法，包括接收包括與地 下特征有關的信息的采樣地震道；使用離散傅立葉變換來估計地震道的傅立葉系數；通過 將相干準則應用于估計的傅立葉系數來選擇估計的傅立葉系數，以減小所述傅立葉系數所 代表的混疊能量的量；在k_空間中重新正交化所述傅立葉系數；以及應用逆傅立葉變換來 產生規則化地震數據。本發明的實施例的一些方面可以包括編有用于執行上述方法或控制上述系統的 計算機可執行指令的計算機可讀介質。本發明的實施例的一些方面可以包括一種合并有上述系統并被配置和安排成依 照上述方法來提供對所述系統的控制的系統。這樣的系統可以合并有例如被編程成允許用 戶依照所述方法或其它方法來控制設備的計算機。通過參照形成本說明書的一部分、相同標號在各個圖形中表示相應部件的附圖對 如下描述和所附權利要求書加以研究，本發明的這些和其它目的、特征、和特性，以及操作 方法、結構的相關元件的功能、部件的組合、和制造成本將變得更加顯而易見。但是，不言而 喻，這些附圖只用于例示和描述的目的，而無意作為限制本發明的定義。正如用在說明書和 權利要求書中的那樣，除非上下文另有明確指明，單數形式“一個”、“一種”、和“該”也包括 復數指示物。
圖1例示了采樣道的集合的例子；圖2例示了圖1的集合的f_k譜；
4[oo10] 圖3例示了在應用依照本發明一個實施例的方法之后圖2的[—k譜；[oo11] 圖4例示了圖3的[—k譜的逆傅立葉變換；以及[oo12] 圖5是執行依照本發明實施例的方法的系統的實施例的示意性例示圖。
具體實施方式
[oo13] 地震數據可被編譯成二維斷面或集合。在圖l中例示了一個例子。一般說來，可以將一組這樣的數據表示成P(t，X)，其中，t和X分別代表時間和空間。在圖l的例子中，X軸代表以米為單位的距離(空間)，而y軸代表以秒為單位的時間。對于這個例子，樣本的間隔是不規則的，并且是稀疏采樣的。為了清楚起見，這些道大部分是無噪聲的，只示出了許多記錄的事件lo，而在每個大振幅事件之外，道基本上是平的，具有零振幅。而且，已經將數據限制在窗口內，使得地震事件相當線性。[oo14] 對于被良好采樣(即，足夠的采樣頻率和在時間上足夠規則)的數據，可以依照方程l，使用快速傅立葉變換([[T)來計算傅立葉變換[oo]5] P(‘o，t)一[T+[P(t，X)](1)[oo16] 其中，[T+是時間傅立葉變換算子。[oo17] 對于每個頻率0，把數據看作X的函數，可以寫出方程2[oo]8] h(X)一P(0，X)(2)[oo19] 如上所述，可以作出數據被均勻采樣的假設，使得可以將這個函數的均勻采樣序列表達成方程3
h。一h(n厶X)，n—o，l，．．．，N—l，(3)
可以將正向離散傅立葉變換(D[T)表達成
權利要求
1.一種處理代表與地下區域有關的信息的采樣地震道的方法，包含 使用離散傅立葉變換來估計地震道的傅立葉系數；通過將相干準則應用于估計的傅立葉系數來選擇估計的傅立葉系數，以減小所述傅立 葉系數所代表的混疊能量的量；在k_空間中重新正交化所述傅立葉系數；以及 使用重新正交化的傅立葉系數來產生地下區域的地震圖像。
2.按照權利要求1所述的方法，進一步包含對于每個頻率，選擇混疊抑制函數值最大和離散傅立葉系數振幅最大的傅立葉系數。
3.按照權利要求1所述的方法，其中，所述相干準則包含抑制函數
4.按照權利要求1所述的方法，其中，所述相干準則包含計算沿著穿過k空間中的原點 的線段的相干性的函數。
5.按照權利要求1所述的方法，其中，所述相干準則包含振幅譜的疊加
6.按照權利要求1所述的方法，進一步包含在重新正交化傅立葉系數之后，應用逆傅 立葉變換以獲得規則化地震數據。
7.一種成像地下特征的方法，包含接收包括與地下特征有關的信息的采樣地震道； 使用離散傅立葉變換來估計地震道的傅立葉系數；通過將相干準則應用于估計的傅立葉系數來選擇估計的傅立葉系數，以減小所述傅立 葉系數所代表的混疊能量的量；在k-空間中重新正交化所述傅立葉系數；以及 應用逆傅立葉變換來產生規則化地震數據。
8.按照權利要求7所述的方法，進一步包含根據所述規則化地震數據來顯示地震道圖像。
9.按照權利要求7所述的方法，進一步包含對于每個頻率，選擇混疊抑制函數值最大和快速傅立葉系數振幅最大的傅立葉系數。
10.按照權利要求7所述的方法，其中，所述相干準則包含抑制函數
11.按照權利要求7所述的方法，其中，所述相干準則包含計算沿著穿過k空間中的原 點的線段的相干性的函數。
12.按照權利要求7所述的方法，進一步包含在確定在所述地下特征附近存在烴資源 的概率的過程中使用所述規則化地震數據。
13. 一種成像與地下特征有關的地震數據的系統，包含具有計算機可讀數據的數據存儲設備，所述計算機可讀數據包括與采樣地震道有關的 數據，所述采樣地震道包括與地下特征有關的信息；以及處理器，被配置和安排成執行存儲在處理器可訪問的存儲器中的機器可執行指令，以 便執行包含如下步驟的方法使用離散傅立葉變換來估計地震道的傅立葉系數；通過將相干準則應用于估計的傅立葉系數來選擇估計的傅立葉系數，以減小所述傅立 葉系數所代表的混疊能量的量；在k-空間中重新正交化所述傅立葉系數；以及 應用逆傅立葉變換來產生規則化地震數據。
全文摘要
處理地震數據以減小或消除由例如稀疏或不規則采樣引起的混疊。迭代方法包括與估計傅立葉系數的振幅的函數結合使用的抑制函數，這兩個函數一起作用以減小混疊能量的影響并優先選擇真正能量。計算步驟主要在k空間中執行，而不返回到x空間，從而減小計算成本。
文檔編號G01V1/32GK102105814SQ200980128648
公開日2011年6月22日 申請日期2009年6月15日 優先權日2008年6月19日
發明者L·L·張 申請人:雪佛龍美國公司