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

在審

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

1.一種疊后地震數據高頻層序自動追蹤方法，其特征在于，所述方法包括以下步驟：

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

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

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

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

2.如權利要求1所述的疊后地震數據高頻層序自動追蹤方法，其特征在于，所述步驟二具體為：

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

3.如權利要求1所述的疊后地震數據高頻層序自動追蹤方法，其特征在于，所述步驟三具體為：

定義目標函數

式(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)的邊界條件，整理得到：

式中，和分別為φ_i,j和ρ_i,j的離散余弦變換，最后對進行離散余弦逆變換即可求得φ_i,j，即最終展開的相位。

4.如權利要求3所述的疊后地震數據高頻層序自動追蹤方法，其特征在于，所述步驟三中，采用分段方式計算相位展開體，是用步驟二中追蹤到的三級層序界面或體系域界面約束相位展開，使相位展開體具有精確等時性。

5.如權利要求1所述的疊后地震數據高頻層序自動追蹤方法，其特征在于，所述步驟四中，相位展開體等值面對應于地震同相軸等時面，因此可以作為等時的層序。