[发明专利]一种采用频变的复数相关系数估算瞬时慢度的方法

说明书

技术领域

本发明涉及地球物理勘探技术，属于地震资料处理中属性提取范畴，是一种采用频变的复数相关系数估算瞬时慢度的方法。

背景技术

在地震资料处理中，瞬时慢度是地震数据一个非常重要的属性。瞬时慢度可广泛应用于地震数据插值或外推、噪声压制、VSP数据上下行波分离、波束偏移、速度层析反演和地震解释中断层预测等。

估算瞬时慢度的方法一般分为三大类。第一类是基于相关性分析的慢度扫描方法；第二类基于局部Radon变换的慢度方法；第三类是基于平面波微分方程的慢度估算方法。

基于相关性分析的慢度扫描方法以当前地震道的当前采样点为中心点，按扫描慢度确定信号在其相邻地震道上的轨迹，然后以信号在地震道上的轨迹点为中心取时窗，计算所有时窗内数据之间的互相关系数并求和，这样就得到了以扫描慢度为变量的互相关函数，最后求互相关函数的最大值对应的扫描慢度就可以得到当前地震道的当前采样点处的瞬时慢度。这种方法的优点是算法简单，缺点是得到的慢度精度低。

发明内容

本发明目的在于提供一种可更准确估算瞬时慢度的采用频变的复数相关系数估算瞬时慢度的方法。

本发明的具体实施步骤如下：

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

2）对地震剖面进行频谱分析并确定地震剖面的主频f_c；

对地震剖面进行F-K谱分析并确定地震剖面的最小和最大视速度v_min和v_max；

3）在时间方向，把地震道投影到以地震剖面的主频f_c为频率的富里叶变换的基函数上，得到投影后的地震剖面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）根据步骤2)确定的地震剖面的主频f_c，用下式计算以采样点为中心的沿时间方向的处理时窗的长度：

wt=1fcΔt]]>

式中，f_c为地震剖面的主频，Δt为地震剖面的时间采样率，w_t为处理时窗的长度；

6）根据步骤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）根据慢度扫描个数，计算慢度扫描的增量：