[發明專利]一種多道約束的地震反射識別方法

有效

• 申請號: 201610984280.3
• 申請日: 2016-11-09
• 公開（公告）號: CN108061916B
• 公開（公告）日: 2020-08-25
• 發明（設計）人: 陳科;王鵬燕;徐雷鳴;龐銳;張楊
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建

本發明公開了一種多道約束的地震反射識別方法，包括：地震反射時間確定步驟，確定該道層位對應的地震反射的時間；初始參數確定步驟，確定初始相位和初始頻率，搜索地震反射子波的初始尺度；局部優化步驟，在局部范圍內搜索最優的相位、頻率和尺度；子波形態確定步驟，由地震反射時間和局部最優的相位、頻率和尺度獲得子波振幅，進而確定地震反射子波的時間域形態。本申請對地震反射時間采取小面元局部最優策略，對子波形態控制參數采取實部和虛部分別匹配策略；二者相結合能夠確保最優地震反射子波形態的準確性和能量變化的連續性。

技術領域

本發明涉及地球物理勘探領域，用于對地震數據進行預處理，具體地說，涉及一種多道約束的地震反射識別方法。

背景技術

中國國內油氣藏成因復雜，煤系地層、多級不整合發育等復雜地質背景，導致油氣層段有效地震反射受強背景反射干擾，給強反射背景下地震儲層預測及含烴檢測帶來了很大的挑戰。

目前相關強反射分離的文獻或專利多是基于多子波分解方法，將原始地震信號分解為多個子波，然后去除能量最強的一個或幾個子波，便得到去除強反射之后的地震信號。但是現有的技術存在一些問題：(1)去除能量較強的一個或幾個子波，相當于只能去除能量比較大的強反射，不能去除一些能量相對較小的同相軸，實際應用受到限制；(2)針對煤層或不整合面等引起的強反射，同相軸比較連續，通過時窗控制去除強反射之后可以得到同相軸比較連續的地震數據，但是針對碎屑巖地區，同相軸能量橫向變化比較快，去除強反射之后的地震數據會出現掛面條現象，影響去除強反射之后地震數據的使用；(3)多子波分解過程都是在時間域搜索局部最佳子波，將子波作為一個整體尋找最優的控制參數，而決定子波的各參數對子波頻譜的實部和虛部影響程度是不同的，在時間域搜索時相當于采用實部、虛部相同的權重，導致搜索的局部最優子波不夠準確；每一次迭代搜索的子波都不準確，隨著迭代的進行，導致后續搜索的子波誤差較大，不具可用性，這也是目前多子波分解方法存在的問題。

強反射分離的核心為多子波分解方法，子波的好壞直接影響到強反射分離后地震數據的質量，那么如何能夠根據地震響應的局部特征獲得更精確、更連續的地震反射子波是目前亟需解決的問題。

發明內容

本發明的目的在于解決現有技術中存在搜索子波的準確性、子波能量連續性和實際應用靈活性的問題，提供一種多道約束的地震反射識別方法，該方法能夠針對任意地震反射同相軸、搜索精確的地震反射子波、能量隨著原始地震數據的能量變化。

本發明提供一種多道約束的地震反射識別方法，包括以下步驟：

地震反射時間確定；

初始參數確定；

局部優化；

子波形態確定。

如上所述，所述地震反射時間確定是由面元內各道的層位與波形的關系確定該道層位處最佳所述地震反射時間μ，具體地，

指定一道地震信號x(t)為中心道及一個局部小面元，統計所述面元內所述各道對應的層位時間與臨近的極值間的關系，離波峰較近的道數為M，離波谷較近的道數為N；比較M和N的大小；若所述M的值較大，則記錄離所述中心道的所述層位最近的所述波峰對應的時間為μ；若所述N的值較大，則記錄離所述中心道的所述層位最近的所述波谷對應的時間為μ。

如上所述，所述初始參數確定是由所述地震信號x(t)確定所述地震反射時間μ處的初始相位φ₀和初始頻率f₀，搜索地震反射子波的初始尺度σ₀，其中，所述初始相位φ₀為所述地震反射時間μ對應的瞬時相位，所述初始頻率f₀為所述地震反射時間μ對應的瞬時頻率。

如上所述，在搜索所述地震信號的所述初始尺度σ₀的步驟中，