[發明專利]一種疊后地震數據濾波方法

有效

• 申請號: 201010176516.3
• 申請日: 2010-05-13
• 公開（公告）號: CN101852866A
• 公開（公告）日: 2010-10-06
• 發明（設計）人: 王彥倉;葉秋焱;張樹森
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 任默聞

技術領域

本發明是關于石油勘探技術，特別是關于地震數據處理技術，具體的講是關于一種疊后地震數據濾波方法。

背景技術

在地震數據的采集過程中所使用的觀測系統會對采集到的地震數據產生干擾，可以形象地將這種干擾被喻為在地震數據中留下了腳印，即采集腳印干擾。這種腳印表現為在疊后地震數據體的時間切片上出現有規律的振幅變化，是調制在地層反射信號上的一種系統性噪聲，它嚴重的影響了地震屬性反演的精確度。

現有壓制采集腳印干擾的技術一般是根據觀測系統，設計出壓制采集腳印干擾的濾波器，在地震數據疊后壓制采集腳印干擾。但是在實際的運用過程中，由于斜交等不規則觀測系統的使用，以及受觀測系統改變，地表起伏(如地面遇到建筑物)等因素的影響，采集腳印干擾的周期性分布會受到擾亂，導致根據觀測系統設計濾波器變得非常困難。

發明內容

本發明提供一種疊后地震數據濾波方法，以有效去除有數據采集、處理過程中引入的腳印干擾。

為了實現上述目的，本發明提供了一種疊后地震數據濾波方法，所述的方法包括：

獲取疊后地震數據；

將獲取的疊后地震數據由時間域轉換到時間切片域；

將所述時間切片域中每一時間切片的地震數據從空間域轉換到波數域；在所述波數域上獲取滿足下述條件的地震數據采集點：

1)在以原點為中心，半徑為R0的圓形區域外；

2)所述地震數據采集點的振幅為局部振幅極大值；

3)所述地震數據采集點的振幅大于全局振幅最大值的d倍；

4)所述地震數據采集點的振幅大于以滿足條件1)、2)、3)的地震數據采集點為中心，半徑為R1的區域內的振幅平均值的f倍；

根據獲取的地震數據采集點構造濾波器；

利用所述的濾波器對所述的地震數據采集點進行濾波。

本發明的有益技術效果：

本發明可以有效的去除由采集、處理過程中引入的采集腳印，提高了信噪比，消除了采集腳印對地震數據屬性的影響，有利于開展地震資料的綜合研究，提高了儲層預測精度，能夠更加準確地刻畫斷裂系統，精確預測砂體的分布范圍和展布特征。

附圖說明

為了更清楚地說明本發明實施例或現有技術中的技術方案，下面將對實施例或現有技術描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發明的一些實施例，對于本領域普通技術人員來講，在不付出創造性勞動性的前提下，還可以根據這些附圖獲得其他的附圖。在附圖中：

圖1為本發明實施例一正演模型的疊后地震數據示意圖；

圖2為本發明實施例一切除第700ms處的時間切片的疊后地震數據示意圖；

圖3為本發明實施例一變換得到的二維波數譜示意圖；

圖4為本發明實施例一二維波數域濾波器示意圖；

圖5為本發明實施例一由二維波數域(Kx，Ky)上被濾波了的數據變換回時間切片上的二維空間域的結果示意圖；

圖6為本發明實施例的方法流程圖。

具體實施方式

為使本發明實施例的目的、技術方案和優點更加清楚明白，下面結合附圖對本發明實施例做進一步詳細說明。在此，本發明的示意性實施例及其說明用于解釋本發明，但并不作為對本發明的限定。

實施例一

輸入正演模型的疊后地震數據，總道數260×300，共中心點(Common?MidPoint?CMP)線數260，每條CMP線有CMP點數300，正演模型的疊后地震數據如圖1所示，圖1顯示了距離地面不同深度的不同疊后地震數據。

將上述正演模型的疊后地震數據由時間域轉換到時間切片域。根據上述疊后地震數據，可以獲得時間點的時間切片，如圖2所示，為第700ms處切除的時間切片。

利用二維傅立葉變換把該時間切片的地震數據由空間域(x，y)轉換到波數域(Kx，Ky)，圖3為由空間域(x，y)轉換到波數域(Kx，Ky)變換得到的二維波數譜示意圖。

在二維波數譜中，按如下條件查找被采集腳印噪音干擾的點(地震數據采集點)，用于構造二維波數域濾波器。

①在以原點(0，0)點為中心半徑為R0(單位：周期數/道)的圓形保護區域外。其中，半徑R0的取值范圍為0.03～0.08(周期數/道)，該取值范圍僅為本發明的實施例，并非用于限定本發明。

②地震數據采集點的振幅為局部振幅極大值的點。

③地震數據采集點的振幅大于全局振幅最大值的d倍；其中，d的取值范圍為0.0007～0.005，該取值范圍僅為本發明的實施例，并非用于限定本發明。