姜博午,劉金朋,陸文凱
(1.北京化工大學機電工程學院,北京100029;2.智能技術與系統國家重點實驗室,北京 100084;3.北京信息科學與技術國家研究中心,北京 100084;4.清華大學自動化系,北京 100084;5.中海油田服務股份有限公司物探事業部特普公司,廣東湛江524057)
工業界常用的基于伴隨算子的地震偏移成像算法假設地震數據中只存在一次反射波(一次波),而多次反射波(多次波)的存在不符合傳統成像算法的假設條件,造成有效波振幅、頻率、相位畸變,影響地震解釋[1]。多次波廣泛存在于陸地和海洋地震數據中,特別是在海洋地震數據中,由于存在海面和海底兩個強反射界面,因而導致海洋地震數據的表面多次波干擾問題更為嚴重。對于陸上地震數據而言,當地下存在強反射界面,諸如基巖面、火成巖侵入體、不整合面等,地震波在強反射界面間或其內部發生多次反射,形成層間多次波[2]。由于多次波傳播路徑更加復雜,因而壓制層間多次波更有挑戰性。多次波壓制方法主要分為兩類:濾波法和預測相減法[3]。濾波類方法利用一次波和多次波在變換域內的差異壓制多次波,包括F-K濾波[4]、Radon變換濾波方法[5-8]等。通常當一次波和多次波速度差異較大時,濾波類方法效果較理想;但當地下介質復雜時,濾波類方法很難在保護有效信號的同時壓制多次波。預測相減法首先基于波動理論得到預測多次波模型,然后通過自適應相減壓制多次波。這類算法更加適用于復雜介質情況。……