[發(fā)明專利]一種采用頻變的復數(shù)相關系數(shù)估算瞬時慢度的方法

有效

• 申請?zhí)枺?/td>: 201310304036.4
• 申請日: 2013-07-19
• 公開（公告）號: CN104297782B
• 公開（公告）日: 2017-03-15
• 發(fā)明（設計）人: 柯本喜
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京康信知識產(chǎn)權代理有限責任公司11240
• 代理人: 吳貴明,張永明

技術領域

本發(fā)明涉及地球物理勘探技術，屬于地震資料處理中屬性提取范疇，是一種采用頻變的復數(shù)相關系數(shù)估算瞬時慢度的方法。

背景技術

在地震資料處理中，瞬時慢度是地震數(shù)據(jù)一個非常重要的屬性。瞬時慢度可廣泛應用于地震數(shù)據(jù)插值或外推、噪聲壓制、VSP數(shù)據(jù)上下行波分離、波束偏移、速度層析反演和地震解釋中斷層預測等。

估算瞬時慢度的方法一般分為三大類。第一類是基于相關性分析的慢度掃描方法；第二類基于局部Radon變換的慢度方法；第三類是基于平面波微分方程的慢度估算方法。

基于相關性分析的慢度掃描方法以當前地震道的當前采樣點為中心點，按掃描慢度確定信號在其相鄰地震道上的軌跡，然后以信號在地震道上的軌跡點為中心取時窗，計算所有時窗內(nèi)數(shù)據(jù)之間的互相關系數(shù)并求和，這樣就得到了以掃描慢度為變量的互相關函數(shù)，最后求互相關函數(shù)的最大值對應的掃描慢度就可以得到當前地震道的當前采樣點處的瞬時慢度。這種方法的優(yōu)點是算法簡單，缺點是得到的慢度精度低。

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

本發(fā)明目的在于提供一種可更準確估算瞬時慢度的采用頻變的復數(shù)相關系數(shù)估算瞬時慢度的方法。

本發(fā)明的具體實施步驟如下：

1）采集得到二維地震剖面；

2）對地震剖面進行頻譜分析并確定地震剖面的主頻f_c；

對地震剖面進行F-K譜分析并確定地震剖面的最小和最大視速度v_min和v_max；

3）在時間方向，把地震道投影到以地震剖面的主頻f_c為頻率的富里葉變換的基函數(shù)上，得到投影后的地震剖面A（x，t，f_c）：

A(x,t,fc)=a(x,t)e-2πifct]]>

式中，a（x，t）是二維地震剖面，x和t分別為地震道的空間和時間坐標，f_c為地震剖面的主頻；A（x，t，f_c）為投影后的地震剖面；

4）對投影后的地震剖面A（x，t，f_c）先按空間道序再按時間采樣點順序逐點按步驟5）-18）估算瞬時慢度；

5）根據(jù)步驟2)確定的地震剖面的主頻f_c，用下式計算以采樣點為中心的沿時間方向的處理時窗的長度：

wt=1fcΔt]]>

式中，f_c為地震剖面的主頻，Δt為地震剖面的時間采樣率，w_t為處理時窗的長度；

6）根據(jù)步驟2)確定的地震剖面的最小和最大視速度v_min和v_max，計算出慢度掃描的最少值和最大值；

p_min=1/v_max

p_max=1/v_min

式中，p_min和p_max表示慢度掃描的最少值和最大值，v_min和v_max是由步驟2）確定的地震剖面的最小和最大視速度；

7）根據(jù)慢度掃描個數(shù)，計算慢度掃描的增量：