[發明專利]基于波動方程延拓的近地表噪音壓制方法

有效

• 申請號: 201811292104.9
• 申請日: 2018-10-31
• 公開（公告）號: CN109212609B
• 公開（公告）日: 2021-06-22
• 發明（設計）人: 李凌云;曲志鵬;王延光;王靜軒;張輝;武天祥;劉瑞合;揭景榮
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 濟南日新專利代理事務所(普通合伙) 37224
• 代理人: 崔曉艷

本發明提供一種基于波動方程延拓的近地表噪音壓制方法，該基于波動方程延拓的近地表噪音壓制方法包括：步驟１，輸入近地表速度厚度模型；步驟２，將數據進行選排，分選為共炮，共檢道集；步驟３，共炮道集對檢波點進行延拓，延拓到低速帶底界面；步驟４，共炮道集對檢波點進行延拓，向上延拓到地表；步驟５，將完成了炮點和檢波點延拓的數據，進行重新選排，重新抽取回共炮道集，完成輸出。該基于波動方程延拓的近地表噪音壓制方法幾乎不損失有效信號就可以將與近地表相關噪音進行壓制，從原理上就可以保護有效波，去除噪音信號。

技術領域

本發明涉及油氣勘探地震資料處理技術領域，特別是涉及到一種基于波動方程延拓的近地表噪音壓制方法。

背景技術

地震資料中的噪音主要與近地表相關，如強面波干擾、強折射波、聲波以及近地表繞射源產生的噪音，這些噪音的存在嚴重影響了地震資料的信噪比。針對不同的噪音特征，目前已經提出了許多種噪音壓制的方法，如低F-K域，Radon變換，曲波變換，小波變換等不同域進行去噪的處理。這些方法都具有各自的特點并取得了明顯的效果，但往往只是側重考慮某一方面的特性或依據某種假設條件，都有其局限性。F-K濾波根據有效信號和噪音的視速度差異進行濾波，在很多實際資料中面波和有效信號視速度差異不夠明顯，造成去噪效果差；Radon變換和F-K濾波器相似，也是依據經過變換后在Randon域噪音和有效信號會形成不同的聚焦能量，但是實際噪音往往無法和有效信號徹底分離，造成去噪效果并不理想；曲波變換和小波變換方法的問題和前述方法類似，均是通過有效信息和噪音的某種差異進行信噪分離，在不能精確分離時，必然傷及有效信號。

由于地震資料信噪比對處理解釋的影響巨大，而目前常用的去噪方法存在諸多問題。隨著油氣勘探工作的不斷深化，地震資料處理中對任何有效信號的傷害必然對后續的處理解釋帶來無法彌補的損失。

地震記錄中噪音種類繁多，強折射、面波、聲波等干擾各具特點，由于這些干擾的存在嚴重影響了地震的信噪比、分辨率，如果疊前不對其進行有效地壓制，最終必然會影響到整個資料處理的質量。在壓制噪音的同時，如何避免對有效信號的損害，一直是地震資料處理的重要內容。為此我們發明了一種新的基于波動方程延拓的近地表噪音壓制方法，解決了以上技術問題。

發明內容

本發明的目的是提供一種最大限度保護有效信息，同時壓制多種噪音干擾的基于波動方程延拓的近地表噪音壓制方法。

本發明的目的可通過如下技術措施來實現：基于波動方程延拓的近地表噪音壓制方法，該基于波動方程延拓的近地表噪音壓制方法包括：步驟1，輸入近地表速度厚度模型；步驟2，將數據進行選排，分選為共炮，共檢道集；步驟3，共炮道集對檢波點進行延拓，延拓到低速帶底界面；步驟4，共炮道集對檢波點進行延拓，向上延拓到地表；步驟5，將完成了炮點和檢波點延拓的數據，進行重新選排，重新抽取回共炮道集，完成輸出。

本發明的目的還可通過如下技術措施來實現：

在步驟1中，輸入的模型為離散網格數據，每個網格點的數值為當前位置縱波速度。

在步驟2中，將數據分別選排為共炮道集以進行檢波點延拓和共檢波點道集以進行共炮點延拓，進行延拓時采用的基本運算公式為：

式中，u表示在某一深度位置記錄波場的傅里葉變換，v表示縱波在介質中的傳播速度，ω表示圓頻率，k_x表示水平方向上的圓波數，ΔZ表示深度域延拓步長，Z表示當前深度位置，公式中指數位置取正為向上延拓，取負值為向下延拓。

在步驟3中，選擇有限差分方法在炮集和檢波點道集記錄中進行波場記錄的向下延拓，將檢波點和炮點延拓的低速帶底界面，實現波場向下延拓過程，在延拓過程，僅保留正時間的地震數據，即延拓過程中對當前延拓面之上的信息進行清除。