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

在審

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

本發(fā)明公開基于二維Teager?Huang變換的疊前地震資料薄層識別方法，涉及油氣物探工程領(lǐng)域。包括如下步驟：S1、對單道集數(shù)據(jù)進行二維經(jīng)驗?zāi)B(tài)算子分解，獲得不同尺度分量；S2、對不同尺度分量算子進行遴選，選擇與原始信號相關(guān)度最大的分量；S3、對分量進行2D?TKEO算子計算；S4、對計算后的數(shù)據(jù)進行Hilbert變換，獲得Hilbert譜；S5、對Hilbert譜進行疊加，獲得時頻分析疊加譜；S6、對疊加譜進行高斯平滑，提取結(jié)果峰值點，對應(yīng)儲層分界面。采用本方案解決了大多數(shù)薄儲層定位識別技術(shù)效用低的問題，提高了薄層識別的準確性。

技術(shù)領(lǐng)域

發(fā)明公開基于二維Teager-Huang變換的疊前地震資料薄層識別方法，涉及油氣物探工程領(lǐng)域。

背景技術(shù)

隨著油氣勘探開發(fā)程度加深，從常規(guī)油氣到非常規(guī)油氣突破，已經(jīng)從淺層轉(zhuǎn)向深層油氣勘探領(lǐng)域中。深層油氣儲層往往儲層厚度比較薄，再加上地震信號經(jīng)過反射、折射，接收到的有效信號相當弱，導(dǎo)致現(xiàn)有的大多數(shù)薄儲層定位識別技術(shù)效用低。并且如果不知道儲層頂?shù)孜恢茫苯佑矛F(xiàn)有的方法進行檢測，很難發(fā)揮方法的有效性，因此薄層定位顯得格外重要。

發(fā)明內(nèi)容

本發(fā)明的基于二維Teager-Huang變換的疊前地震資料薄層識別方法，解決了大多數(shù)薄儲層定位識別技術(shù)效用低的問題。

為達到上述目的，采用的技術(shù)方案如下：

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

S1、對單道集數(shù)據(jù)進行二維經(jīng)驗?zāi)B(tài)算子分解，獲得不同尺度分量；

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

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

S4、對計算后的數(shù)據(jù)進行Hilbert變換，獲得Hilbert譜；

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

S6、對疊加譜進行高斯平滑，提取結(jié)果峰值點，對應(yīng)儲層分界面。

優(yōu)選的，步驟S1的二維經(jīng)驗?zāi)B(tài)算子分解方法如下：首先將疊前道集單炮記錄比作是一個M×N的矩陣，即f(x,y),x＝1,…,M；y＝1,…,N；

算法流程如下：

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

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

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

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

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

優(yōu)選的，步驟S1.2的方法如下：

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

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

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

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

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