[發明專利]地震數據的恢復方法

有效

• 申請號: 201410100768.6
• 申請日: 2014-03-18
• 公開（公告）號: CN103901469A
• 公開（公告）日: 2014-07-02
• 發明（設計）人: 劉奇琳;杜金虎;何光明;陳愛萍;金德剛;陳三平;張恩嘉;唐虎;周強
• 申請（專利權）人: 中國石油集團川慶鉆探工程有限公司地球物理勘探公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京銘碩知識產權代理有限公司 11286
• 代理人: 張川緒;王兆賡

1.一種地震數據的恢復方法，所述方法包括：對時間空間域的地震數據進行掃描，得到缺失道地震數據和有效道地震數據，生成一個數據窗，所述數據窗內的缺失道地震數據與數據窗內的地震數據的比值在設定閾值范圍內，對數據窗內的時間空間域的地震數據執行以下步驟：

（a）對時間空間域的地震數據進行傅里葉變換，得到頻率空間域的地震數據；

（b）對頻率空間域的地震數據進行恢復，得到恢復后的頻率空間域的地震數據；

（c）對恢復后的頻率空間域的地震數據進行傅里葉反變換，得到時間空間域的恢復后的地震數據。

2.根據權利要求1所述的方法，其中，對頻率空間域的地震數據進行恢復的步驟包括：

（b1）對頻率空間域的地震數據進行傅里葉變換，得到頻率波數譜；

（b2）求取頻率波數譜的最大半徑，根據頻率波數譜的最大半徑和預定變量k得到權函數的范圍界限；

（b3）利用權函數的范圍界限對頻率波數譜進行處理，得到權函數處理后的頻率波數譜；

（b4）求取權函數處理后的頻率波數譜的模值，從所述模值中找到權函數處理后的頻率波數譜的最大模值和最小模值，根據權函數處理后的頻率波數譜的最大模值和最小模值和預定變量k，得到閾值函數；

（b5）利用閾值函數對權函數處理后的頻率波數譜進行處理，得到閾值函數處理后的頻率波數譜；

（b6）對閾值函數處理后的頻率波數譜進行傅里葉反變換，得到頻率空間域的閾值函數處理后的地震數據；

（b7）利用步驟（a）得到的頻率空間域的地震數據中的有效道地震數據取代頻率空間域的閾值函數處理后的地震數據中對應道編號的地震數據；

（b8）使得k=k+1，檢測k是否等于m，m為大于零的自然數；

（b9）當k≠m時，返回執行步驟（b1）～步驟（b8）；

（b10）當k=m時，執行步驟（c），

其中，預定變量k的初始值為1。

3.根據權利要求1所述的方法，還包括：當對一個數據窗內的時間空間域的地震數據執行步驟（a）、（b）、（c）之后，將所述數據窗沿空間域的一個方向移動一個窗口的距離，然后對移動后的數據窗內的時間空間域的地震數據執行步驟（a）、（b）、（c）。

4.根據權利要求1所述的方法，還包括：為所述數據窗增加邊框，其中，對增加了邊框的數據窗內的時間空間域的地震數據執行步驟（a）、（b）、（c）。

5.根據權利要求4所述的方法，其中，所述邊框被填充零值。

6.根據權利要求4所述的方法，其中，增加的邊框的寬度為10道。

7.根據權利要求2所述的方法，其中，步驟（b2）中計算權函數的范圍界限的公式如下：

Mk=k·Rmaxm]]>

其中，M_k為權函數的范圍界線，R_max為頻率波數譜的最大半徑。

8.根據權利要求2所述的方法，其中，步驟（b3）包括：對在所述權函數的范圍界限內的頻率波數譜保留，對在所述權函數的范圍界限外的頻率波數譜去除，將在所述權函數的范圍界限內的頻率波數譜作為權函數處理后的頻率波數譜。