[發(fā)明專利]地震剖面的經(jīng)驗(yàn)?zāi)B(tài)分解方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201310203002.6
• 申請(qǐng)日: 2013-05-27
• 公開(kāi)（公告）號(hào): CN103308940A
• 公開(kāi)（公告）日: 2013-09-18
• 發(fā)明（設(shè)計(jì)）人: 馮偉;馮佩華
• 申請(qǐng)（專利權(quán)）人: 中山大學(xué)
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28
• 代理公司: 廣州粵高專利商標(biāo)代理有限公司 44102
• 代理人: 林麗明

1.一種地震剖面的經(jīng)驗(yàn)?zāi)B(tài)分解方法，其特征在于，包括以下步驟：

S1.初始化，r₀=Iori，i=1；給定迭代終止條件中ε的值；Iori表示原始剖面數(shù)據(jù)，r代表當(dāng)前殘差剖面數(shù)據(jù)，i為第i個(gè)BIMF；

S2.迭代篩分計(jì)算第i個(gè)BIMF；

S21.初始化：h₀=r_i-1，j=1；

S22.檢測(cè)h_j-1的二維局部極小值點(diǎn)和極大值點(diǎn)，構(gòu)造上下包絡(luò)曲面E_up和E_down；

S23.計(jì)算均值曲面m=(E_up+E_down)/2；

S24.h_j=h_j-1-m；

S25.判斷當(dāng)前h是否滿足篩分終止條件，若滿足，則得到第i個(gè)BIMFb_i=h_j，跳轉(zhuǎn)到步驟S3；否則，j=j+1，跳轉(zhuǎn)到步驟S22，繼續(xù)篩分迭代；

S3.計(jì)算當(dāng)前殘差，r_i=r_i-1-b_i；

S4.判斷r_i是否滿足BIMF終止條件，若滿足，則分解結(jié)束，記I=i，r_i是最終殘余分量；否則i=i+1，跳轉(zhuǎn)到步驟S2；

BEMD計(jì)算結(jié)果為：

上述采用時(shí)間方向兩點(diǎn)鄰域的極值檢測(cè)方法檢測(cè)h_j-1的極值；

上述采用等式約束的緊框架域稀疏化曲面插值方法獲取上下包絡(luò)曲面。

2.根據(jù)權(quán)利要求1所述的地震剖面的經(jīng)驗(yàn)?zāi)B(tài)分解方法，其特征在于，所述采用時(shí)間方向兩點(diǎn)鄰域的極值檢測(cè)方法具體為：若某點(diǎn)(x₀,y₀)處的波場(chǎng)值大于或等于其上、下兩個(gè)相鄰位置的波場(chǎng)值，則點(diǎn)(x₀,y₀)為剖面的極大值；相應(yīng)地，若某點(diǎn)(x₀,y₀)處的波場(chǎng)值小于或等于其上、下兩個(gè)相鄰位置的波場(chǎng)值，則點(diǎn)(x₀,y₀)處波場(chǎng)值為剖面的極小值。