[發明專利]一種疊后地震數據高頻層序自動追蹤方法

在審

• 申請號: 201910270711.3
• 申請日: 2019-04-04
• 公開（公告）號: CN109991659A
• 公開（公告）日: 2019-07-09
• 發明（設計）人: 劉力輝
• 申請（專利權）人: 北京諾克斯達石油科技有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京匯信合知識產權代理有限公司 11335
• 代理人: 孫民興

本發明公開了一種疊后地震數據高頻層序自動追蹤方法，其包括以下步驟：步驟一、對疊后地震數據體去噪并計算瞬時相位體；步驟二、追蹤三級層序界面或體系域界面；步驟三、以步驟二中追蹤到的三級層序界面或體系域界面為格架，分段計算步驟一中所計算出的瞬時相位體的相位展開體，然后合并；步驟四、在合并后的相位展開體上提取等值面，作為高頻層序。采用本發明的方法，可自動追蹤高頻層序，效率高，且具有精確等時性。

技術領域

本發明涉及石油物探技術領域，尤其涉及一種疊后地震數據高頻層序自動追蹤方法。

背景技術

在地震數據解釋中，高頻層序是指四、五級層序。精確等時的高頻層序對于利用等時切片進行沉積分析和刻畫沉積體具有重要意義。

目前高頻層序解釋主要是人工解釋。根據工區內測井資料建立工區的高頻層序，或者通過定義簡單的平行關系確定高頻層序。這樣就完成了高頻層序解釋。

傳統高頻層序解釋的方法，若是通過測井資料，則依賴于人的解釋經驗，具有主觀性，并且耗時耗力；若通過定義平行關系確定，其與地震同相軸不一定一致，因此效率高但精確性低。

發明內容

針對上述現有技術中存在的不足之處，本發明提供一種疊后地震數據高頻層序自動追蹤方法。該疊后地震數據高頻層序自動追蹤方法包括以下步驟：

步驟一、對疊后地震數據體去噪并計算瞬時相位體；

步驟二、追蹤三級層序界面或體系域界面；

步驟三、以步驟二中追蹤到的三級層序界面或體系域界面為格架，分段計算步驟一中所計算出的瞬時相位體的相位展開體，然后合并；

步驟四、在合并后的相位展開體上提取等值面，作為高頻層序。

進一步地，所述步驟二具體為：

通過拾取的種子點，用線性插值方法或者自動追蹤方法，追蹤出三級層序界面或體系域界面。

優選地，所述步驟三具體為：

定義目標函數

式(1)中，φ_ij為真實相位Φ第i行第j列的元素，為未知量，和分別表示瞬時相位的差分，即為：

其中，瞬時相位ψ，大小為M×N，上標x和y分別表示水平和豎直方向上的差分，W為包裹算子，為一次2π的取模運算；

所述目標函數中φ是未知量，通過最小化目標函數可以得到；

φ是全局范圍內的最優解，滿足下面的關系：

(φ_i+1,j-2φ_i,j+φ_i-1,j)+(φ_i,j-1-2φ_i,j+φ_i,j-1)＝ρ_i,j (4)

其中：

式(5)即為Poisson方程

的離散形式；求解形如式(6)的離散Poisson方程，定義Neumann邊界條件：

在上式所規定的邊界條件下，按照式(7)對瞬時相位進行計算得到ρ_i,j，之后對等式(4)兩邊同時進行離散余弦變換(DCT)，并結合式(7)的邊界條件，整理得到：