[發明專利]一種高效的地震數據隨機噪聲衰減方法

有效

• 申請號: 201310356433.6
• 申請日: 2013-08-15
• 公開（公告）號: CN104375184B
• 公開（公告）日: 2017-04-05
• 發明（設計）人: 賈繼偉;周興元;羅國安
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京康信知識產權代理有限責任公司11240
• 代理人: 吳貴明,張永明

1.一種高效的地震數據隨機噪聲衰減方法，特點是采用以下步驟實現：

1）采集三維疊前地震數據；

2）根據地震數據空間各方向反射波復雜程度的不同，分別給定各方向的預測算子長度，根據預測算子長度給定各空間方向計算窗長度，并給定時間窗大小；

3）根據各方向計算窗大小將數據分成小計算窗；

4）在當前計算窗內沿時間方向做Fourier變換得到頻率空間域原始數據{s_i,j,k;f}，其中i,j,k分別表示X、Y、O三個空間方向，f表示頻率；

5）用下式計算三維預測算子序列{x_l,m,n}：

Σl,m,n(Σi,j,ksi-l,j-m,k-n;f·s‾i-p,j-q,k-r;f)xl,m,n=Σi,j,ksi,j,k;f·s‾i-p,j-q,k-r;f]]>

其中，l,m,n;p,q,r為自然數,變化范圍分別為1,…,LX;1,…,LY;1,…,LO；

LX、LY、LO分別是X、Y、O方向預測算子長度；

6）利用預測算子序列{x_l,m,n}與頻率空間域原始數據{s_i,j,k;f}進行褶積，得到該計算窗隨機噪聲衰減結果{s^O_i,j,k;f}；

7）利用下式計算得到權系數q_f：

θf=Σi,j,k|si,j,k;f|2/Σi,j,k|sOi,j,k;f|2]]>

8）將計算窗隨機噪聲衰減結果{s^O_i,j,k;f}進行加權得到該計算窗隨機噪聲衰減加權結果

9）對當前窗以及前一計算窗進行均值處理，或者線性插值；

10）在該窗內做反Fourier變換，得到當前計算窗時間－空間域隨機噪聲衰減結果；

11）當前計算窗如果是最末計算窗則完成隨機噪聲衰減；如不是最末計算窗則跳至下一計算窗，并返回步驟4）循環，直至最末計算窗結束。