[发明专利]一种消除地震数据单频干扰的方法

权利要求书

1.一种消除地震勘探单频干扰的方法，其特征在于采用如下技术方案，包括以下步骤：

1)激发地震震源，采集地震数据并做预处理；

2)确定单频干扰的初始频率；

3)由快速频率扫描算法确定不等间隔频率采样频率f_k；

4)计算不等间隔频率采样线性调频谱频率样本X_k；采用如下计算方法：

X_k＝(g^Tc(k)-h^Ts(k))-j(h^Tc(k)+g^Ts(k))       (1)

其中

g^T＝(g₀，g₁，K，g_N-1)(2)

h^T＝(h₀，h₁，K，h_N-1)(3)

c^T(k)＝(c₀(k)，c₁(k)，K，c_N-1(k))(4)

s^T(k)＝(s₀(k)，s₁(k)，K，s_N-1(k))(5)

和

g_n＝x_ncosω₀n(6)

h₀＝x_nsinω₀n(7)

c_n(k)＝cosω^kn(8)

s_n(k)＝sinω^kn    (9)

其中

x_n-----表示地震数据样本，由输入地震数据确定；

ω₀-----表示单频干扰初始圆频率，由步骤2)确定；

ω^k -----表示第k个频率样本的圆频率，由频率算法确定；

s_n(k)-----表示圆频率为ω^k的正弦函数第n个样本，由公式(9)计算确定；

c_n(k)-----表示圆频率为ω^k的余弦函数第n个样本，由公式(8)计算确定；

h_n-----表示地震数据与初始圆频率正弦函数乘积的第n个样本，由公式(7)计算确定；

g_n-----表示地震数据与初始圆频率余弦函数乘积的第n个样本，由公式(6)计算确定；

s(k)-----表示圆频率为ω^k的正弦函数向量，由公式(5)构成；

c(k)-----表示圆频率为ω^k的余弦函数向量，由公式(4)构成；

h   -----表示地震数据与初始圆频率正弦函数乘积向量，由公式(3)构成；

g   -----表示地震数据与初始圆频率余弦函数乘积向量，由公式(2)构成；

N  -----表示输入地震数据样本个数，由输入地震数据确定；

M  -----表示线性调频谱频率样本个数，由频率算法确定；

5)计算线性调频谱样本Y_k，公式为：

Y_k＝|X_k|，k＝1，2，3，K，M    (10)

6)将线性调频谱频率样本中最大值所对应的频率确定为单频干扰的频率f；

7)计算单频干扰余弦函数和正弦函数的振幅A和B；

8)根据单频干扰频率f、单频干扰余弦函数和正弦函数的振幅A和B、地震记录的时间采样间隔Δt，按照以下公式计算单频干扰y_i；

9)计算消除单频干扰后的地震有效信号：

10)根据得出消除单频干扰的地震数据绘制消除单频干扰后的地震数据剖面和存储消除单频干扰后的地震数据。

2.根据权利要求1所述的一种消除地震勘探单频干扰的方法，其特征在于步骤1)所述的预处理是对地震数据进行置标签、定义观测系统的处理。