[發明專利]一種基于自適應最優核時頻分布的地震衰減定性估計方法

無效

• 申請號: 201310272627.8
• 申請日: 2013-07-01
• 公開（公告）號: CN103364832A
• 公開（公告）日: 2013-10-23
• 發明（設計）人: 高靜懷;王曉凱;陳文超
• 申請（專利權）人: 西安交通大學
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 西安通大專利代理有限責任公司 61200
• 代理人: 蔡和平

1.一種基于自適應最優核時頻分布的地震衰減定性估計方法，其特征在于，包括以下步驟：

1）采集原始地震數據

采集原始地震數據的道號、線號、以及采樣點號時間，對采集到的原始地震數據進行預處理，得到疊后地震數據體，記為S(x,y,t)，其中x表示道號，y表示線號，t表示采樣點號時間；

2）拾取地震數據目標層的層位

拾取目標層位及目標層上方的層位，記目標層位為H_T(x,y)，目標層位上面的層位為H_A(x,y)；

3）確定高頻頻率f_H和低頻頻率f_L

首先，對目的層鄰近的地震數據進行頻譜分析，然后選擇高頻頻率f_H和低頻頻率f_L；

4）對疊后地震數據體中每道信號計算自適應最優核時頻分布并得到高頻分量AOKTFR(x,y,t,f_H)和低頻分量AOKTFR(x,y,t,f_L)；

5）確定修正因子α（x，y）

α(x,y)=AOKTFR(x,y,HA(x,y),fL)AOKTFR(x,y,HA(x,y),fH)---(1);]]>

6）估計地震衰減

利用得到的修正因子α（x，y）、高頻分量AOKTFR(x,y,t,f_H)以及低頻分量AOKTFR(x,y,t,f_L)估計目標層位H_T(x,y)附近的衰減SA（x,y,t）：

SA(x,y,t)=AOKTFR(x,y,t,f_L)-α(x,y)AOKTFR(x,y,t,f_H)????（2）。

2.根據權利要求1所述的基于自適應最優核時頻分布的地震衰減定性估計方法，其特征在于，所述確定高頻頻率f_H和低頻頻率f_L具體包括：

首先，對目的層鄰近地震數據進行頻譜分析，得到振幅譜，然后，選取高頻的幅度與低頻的幅度接近、高頻與低頻之間的頻率間隔大、并且高頻與低頻都在地震子波的頻帶范圍之內的高頻頻率和低頻頻率為所需高頻頻率f_H和低頻頻率f_L。

3.根據權利要求1所述的基于自適應最優核時頻分布的地震衰減定性估計方法，其特征在于，所述對疊后地震數據體中每道信號計算自適應最優核時頻分布并得到高頻分量AOKTFR(x,y,t,f_H)和低頻分量AOKTFR(x,y,t,f_L)具體為：

定義信號s(t)的對稱短時模糊函數SSTAF為：

SSTAF(t:θ,τ)=12π∫-∞+∞s*(u-τ/2)w*(u-t-τ/2)s(u+τ/2)w(u-t+τ/2)e-jθudu]]>

=12π∫t-T+|τ/2|t+T-|τ/2|s*(u-τ/2)s(u+τ/2)e-jθudu---(3)]]>

其中，w是寬度為2T的對稱矩形窗函數,τ為時延，θ為多普勒頻率，利用SSTAF和最優化過程可以得到時變最優核函數其中r=τ2+θ2,]]>

在時變最優核函數得到以后，最優核時頻分布如下：

AOKTFR(t,ω)=12π∫-∞+∞∫-∞+∞SSTAF(t;θ,τ)Φopt(t;θ,τ)ejθt-jτωdθdτ---(4)]]>

記S(x,y,t)的自適應最優核時頻分布的結果為AOKTFR(x,y,t,f)，得到的高頻分量為AOKTFR(x,y,t,f_H)，低頻分量為AOKTFR(x,y,t,f_L)。