[發明專利]基于二維Teager-Huang變換的疊前地震資料薄層識別方法

在審

• 申請號: 202310520320.9
• 申請日: 2023-05-10
• 公開（公告）號: CN116482761A
• 公開（公告）日: 2023-07-25
• 發明（設計）人: 蔣旭東;有楚鵬;趙泠森
• 申請（專利權）人: 成都理工大學
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/32;G06F17/14
• 代理公司: 重慶市信立達專利代理事務所(普通合伙) 50230
• 代理人: 于躍

1.基于二維Teager-Huang變換的疊前地震資料薄層識別方法，其特征在于：包括如下步驟：

S1、對單道集數據進行二維經驗模態算子分解，獲得不同尺度分量；

S2、對不同尺度分量算子進行遴選，選擇與原始信號相關度最大的分量；

S3、對分量進行2D-TKEO算子計算；

S4、對計算后的數據進行Hilbert變換，獲得Hilbert譜；

S5、對Hilbert譜進行疊加，獲得時頻分析疊加譜；

S6、對疊加譜進行高斯平滑，提取結果峰值點，對應儲層分界面。

2.如權利要求1所述基于二維Teager-Huang變換的疊前地震資料薄層識別方法，其特征在于：步驟S1的二維經驗模態算子分解方法如下：首先將疊前道集單炮記錄比作是一個M×N的矩陣，即f(x,y),x＝1,…,M；y＝1,…,N；

算法流程如下：

S1.1、初始化：將處理的疊前道集矩陣設置為：

r₀(x,y)＝f(x,y)；

S1.2、進行BEMD分解獲取Bimf分量；

S1.3、疊前單道集f(x,y)經分解為N層分量后得：

其中Bimf_i(x,y)為分解出的第i個Bimf分量，r_N(x,y)為余量。

3.如權利要求2所述基于二維Teager-Huang變換的疊前地震資料薄層識別方法，其特征在于：步驟S1.2的方法如下：

S1.2.1、初始化：設置變量j為分解的Bimf層數j＝1，定義初始化變量h：

h(x,y)＝r_j-1(x,y),

S1.2.2、計算獲取h(x,y)的局部極大值點與局部極小值點；

S1.2.3、對所求取的極大值點和極小值點分布進行包絡擬合，獲得二維極大值包絡曲面u_max(x,y)和二維極小值包絡曲面u_min(x,y)；

S1.2.4、求取極大值和極小值包絡曲面的均值，得到中間過度量：

S1.2.5、利用原始矩陣減去中間過度量，則得到對應的Bimf分量：

Bimf_i＝h(x,y)-m_j-1(x,y)；

S1.2.6、求取余量并作為下一次輸入量：

r_j(x,y)＝r_j-1(x,y)-Bimf_i(x,y)；

S1.2.7、令j＝j+1，重復步驟S1-S6，依次得到Bimf分量，直至滿足終止條件：

跳出循環得到最后的余量。

4.如權利要求1-3中任一所述基于二維Teager-Huang變換的疊前地震資料薄層識別方法，其特征在于：步驟S3的方法如下：定義如下：

其中I＝Bimf_i(x,y)屬于BEMD算子計算后的分量，表示梯度算子：

根據上述定義，可以直接得出以下結果：

其中或表示對一維信號f₀(t)的TKEO算子：

因此，二維TKEO算子是對二維信號沿x和y方向進行一維TKEO計算的結果，定義為：

其中：I(x_i,y_j)，0≤i≤(N-1),0≤j≤(M-1)表示二維信號，表示沿x軸方向的TKEO算子，表示沿y軸方向的TKEO算子，所以I在(x_i,y_j)點處的2D-TKEO算子定義為：