[發(fā)明專利]一種基于小波變換的高頻延拓方法

在審

• 申請?zhí)枺?/td>: 201410475359.4
• 申請日: 2014-09-17
• 公開（公告）號: CN105487115A
• 公開（公告）日: 2016-04-13
• 發(fā)明（設(shè)計(jì)）人: 謝金娥;段心標(biāo);孫晶梅
• 申請（專利權(quán)）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術(shù)研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/28;G06F19/00
• 代理公司: 北京思創(chuàng)畢升專利事務(wù)所 11218
• 代理人: 郭韞

1.一種基于小波變換的高頻延拓方法，其特征在于：所述方法包括：

(1)逐道輸入地震數(shù)據(jù)，對地震數(shù)據(jù)進(jìn)行連續(xù)小波變換，從而將地震數(shù)據(jù) 分解在不同小波尺度上，并獲得主軸頻率；

(2)按倍頻程逐個頻段延拓：在各個小波尺度上，對主軸頻率按倍頻程在 小波域進(jìn)行頻譜延拓，得到主諧波和一系列延拓諧波；

(3)重構(gòu)地震數(shù)據(jù)：將主諧波和延拓諧波利用小波變換反變換進(jìn)行重構(gòu)， 得到高頻延拓后的地震數(shù)據(jù)。

2.根據(jù)權(quán)利要求1所述的基于小波變換的高頻延拓方法，其特征在于：所述步 驟(1)是這樣實(shí)現(xiàn)的：

首先，利用Morlet小波對地震數(shù)據(jù)進(jìn)行連續(xù)小波變換：對每一道地震數(shù)據(jù) 逐道進(jìn)行小波變換，將地震數(shù)據(jù)分解在不同小波尺度上；

然后，選取優(yōu)勢區(qū)域數(shù)據(jù)進(jìn)行分析，獲得小波變換的頻率范圍和期望延拓的 頻譜范圍，從而確定主軸頻率值。

3.根據(jù)權(quán)利要求2所述的基于小波變換的高頻延拓方法，其特征在于：所述 Morlet小波的時間域和頻率率顯示如下：

W(τ,s)=∫-∞∞f(t)1|s|ψ*(t-τs)dt---(4)]]>

其中，W是小波變換后的地震信號，f(t)是輸入的地震信號，t是時間，τ＞0 是尺度因子，s反映位移，其值為正或負(fù)，ψ(t)是基本小波。

4.根據(jù)權(quán)利要求2所述的基于小波變換的高頻延拓方法，其特征在于：所述步 驟(2)是這樣實(shí)現(xiàn)的：

外推以主軸頻率的小波為基準(zhǔn)，按倍頻程向高頻端延拓，得到延拓諧波，具 體如下：

振幅用每個諧波振幅替換超諧波振幅，相位用來更新，這樣 便得到主諧波和一系列延拓諧波；

其中是主軸頻率相位，是噪聲信號相位。

5.根據(jù)權(quán)利要求4所述的基于小波變換的高頻延拓方法，其特征在于：所述步 驟(3)中的小波變換反變換如下：

f(t)=1CΨ∫0∞∫-∞∞1sW(τ,s)ψ(t-τs)dsdτs2---(6)]]>

其中，f(t)是重構(gòu)的地震信號，W是小波變換后得到主諧波和延拓諧波信號， t是時間，τ＞0是尺度因子，s反映位移，其值為正或負(fù)，ψ(t)是基本小波。