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

摘要

本发明涉及油田的勘探、开发、开采技术，是消除地震数据单频干扰的方法。由快速频率扫描算法确定不等间隔频率采样频率，计算不等间隔频率采样线性调频谱频率样本、线性调频谱样本，将线性调频谱频率样本中最大值所对应的频率确定为单频干扰的频率，根据单频干扰余弦函数和正弦函数的振幅、地震记录的时间采样间隔，计算消除单频干扰后的地震有效信号，本发明速度快，不损害有效波，提高了频率分量的信噪比，本发明也可以消除地震数据中由周期性震动(如发电机等)产生的单频干扰。

主权项

一种消除地震勘探单频干扰的方法，其特征在于采用如下技术方案，包括以下步骤：1)激发地震震源，采集地震数据并做预处理；2)确定单频干扰的初始频率；3)由快速频率扫描算法确定不等间隔频率采样频率fk；4)计算不等间隔频率采样线性调频谱频率样本Xk；采用如下计算方法：Xk＝(gTc(k)‑hTs(k))‑j(hTc(k)+gTs(k))       (1)其中gT＝(g0，g1，K，gN‑1)(2)hT＝(h0，h1，K，hN‑1)(3)cT(k)＝(c0(k)，c1(k)，K，cN‑1(k))(4)sT(k)＝(s0(k)，s1(k)，K，sN‑1(k))(5)和gn＝xncosω0n(6)h0＝xnsinω0n(7)cn(k)＝cosωkn(8)sn(k)＝sinωkn    (9)其中xn‑‑‑‑‑表示地震数据样本，由输入地震数据确定；ω0‑‑‑‑‑表示单频干扰初始圆频率，由步骤2)确定；ωk ‑‑‑‑‑表示第k个频率样本的圆频率，由频率算法确定；sn(k)‑‑‑‑‑表示圆频率为ωk的正弦函数第n个样本，由公式(9)计算确定；cn(k)‑‑‑‑‑表示圆频率为ωk的余弦函数第n个样本，由公式(8)计算确定；hn‑‑‑‑‑表示地震数据与初始圆频率正弦函数乘积的第n个样本，由公式(7)计算确定；gn‑‑‑‑‑表示地震数据与初始圆频率余弦函数乘积的第n个样本，由公式(6)计算确定；s(k)‑‑‑‑‑表示圆频率为ωk的正弦函数向量，由公式(5)构成；c(k)‑‑‑‑‑表示圆频率为ωk的余弦函数向量，由公式(4)构成；h   ‑‑‑‑‑表示地震数据与初始圆频率正弦函数乘积向量，由公式(3)构成；g   ‑‑‑‑‑表示地震数据与初始圆频率余弦函数乘积向量，由公式(2)构成；N  ‑‑‑‑‑表示输入地震数据样本个数，由输入地震数据确定；M  ‑‑‑‑‑表示线性调频谱频率样本个数，由频率算法确定；5)计算线性调频谱样本Yk，公式为：Yk＝|Xk|，k＝1，2，3，K，M    (10)6)将线性调频谱频率样本中最大值所对应的频率确定为单频干扰的频率f；7)计算单频干扰余弦函数和正弦函数的振幅A和B；8)根据单频干扰频率f、单频干扰余弦函数和正弦函数的振幅A和B、地震记录的时间采样间隔Δt，按照以下公式计算单频干扰yi；9)计算消除单频干扰后的地震有效信号：10)根据得出消除单频干扰的地震数据绘制消除单频干扰后的地震数据剖面和存储消除单频干扰后的地震数据。