[發(fā)明專利]一種確定地震數(shù)據(jù)反射界面位置的方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201410206701.0
• 申請(qǐng)日: 2014-05-15
• 公開(kāi)（公告）號(hào): CN105093305A
• 公開(kāi)（公告）日: 2015-11-25
• 發(fā)明（設(shè)計(jì)）人: 楊麗;畢進(jìn)娜;周中彪;李呈呈;朱童;胡華鋒
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油化工股份有限公司;中國(guó)石油化工股份有限公司石油物探技術(shù)研究院
• 主分類號(hào): G01V1/30
• 分類號(hào): G01V1/30;G01V1/28;G06F19/00
• 代理公司: 北京思創(chuàng)畢升專利事務(wù)所 11218
• 代理人: 郭韞

1.一種確定地震數(shù)據(jù)反射界面位置的方法，其特征在于：所述方法包括：

(1)設(shè)地震波為f(t)，對(duì)其進(jìn)行Hilbert變換，求取實(shí)地震信號(hào)的解析信號(hào)，得到地震波f(t)的瞬時(shí)振幅、瞬時(shí)相位剖面θ(t)和瞬時(shí)頻率；

(2)基于地震波瞬時(shí)相位剖面θ(t)進(jìn)行一階差分運(yùn)算，得到瞬時(shí)相位一階差分G(t)，則

G(t)＝θ(t+1)-θ(t)；

(3)對(duì)地震數(shù)據(jù)f(t)進(jìn)行極性反轉(zhuǎn)處理得到F(t)＝-f(t)，對(duì)F(t)進(jìn)行Hilbert變換，得到瞬時(shí)相位剖面θ₁(t)，再對(duì)瞬時(shí)相位剖面θ₁(t)進(jìn)行一階差分運(yùn)算，得到極性反轉(zhuǎn)后的地震波瞬時(shí)相位一階差分j(t)，對(duì)其求相反數(shù)得到J(t)＝-j(t)；

(4)將步驟(2)和步驟(3)的結(jié)果作為輸入，利用下式對(duì)兩個(gè)結(jié)果進(jìn)行加權(quán)相加運(yùn)算：

M(t)＝a₁*G(t)+a₂*J(t)

其中a₁+a₂＝1，a₁代表地震數(shù)據(jù)中波谷占的權(quán)重，a₂代表地震數(shù)據(jù)中波峰占的權(quán)重；

(5)將步驟(4)的結(jié)果M(t)與原始地震波f(t)相乘得到新的屬性Q(t)：Q(t)＝M(t)*f(t)，并選取閾值濾去新的屬性Q(t)中的較低值，獲得準(zhǔn)確的反射界面位置。

2.根據(jù)權(quán)利要求1所述的確定地震數(shù)據(jù)反射界面位置的方法，其特征在于：所述新的屬性Q(t)是由不同振幅的脈沖組成，其中小振幅脈沖為干擾信號(hào)，而大振幅脈沖為真實(shí)反射界面的位置。

3.根據(jù)權(quán)利要求1所述的確定地震數(shù)據(jù)反射界面位置的方法，其特征在于：所述選取閾值濾去新的屬性Q(t)中的較低值，獲得準(zhǔn)確的反射界面位置是這樣實(shí)現(xiàn)的：

通過(guò)實(shí)驗(yàn)選取閾值；

將新的屬性Q(t)中的振幅小于或等于閾值的脈沖濾去，得到的結(jié)果即為準(zhǔn)確的反射界面位置。