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

权利要求书

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

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

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

3)计算谐波干扰余弦函数的振幅A；

4)根据地震数据向量x、谐波干扰余弦函数的振幅A、谐波干扰余弦函数向量c，计算目标函数目Q，标函数Q的计算公式为：

Q＝(x-Ac)^T(x-Ac)＝x^Tx-2Ax^Tc+A²c^Tc    (4)

5)根据地震数据向量x、谐波干扰余弦函数的振幅A、谐波干扰余弦函数向量c、谐波干扰正弦函数向量s、时延系数矩阵W、地震记录时间采样间隔Δt，计算目标函数目一阶导数和二阶导数和

 其中：Q是目标函数，和是目标函数对频率的一阶导数和目标函数对频率的二阶导数；和是目标函数、目标函数对时延的一阶导数和目标函数对时延的二阶导数；

目标函数Q对频率f的一阶导数为

▿fQ=4πAΔtxTWs-4πA2ΔtcTWs---(5)]]>

目标函数Q对频率f的二阶导数为

▿f2Q=8π2AΔt2xTWTWc-8π2A2Δt2(c+s)TWTW(c-s)---(6)]]>

目标函数Q对时延τ的一阶导数为

▿τQk=4πAfΔtxTs-4πA2fΔtcTs---(7)]]>

目标函数Q对时延τ的二阶导数为

▿τ2Q=8π2A2f2Δt2xTc-8π2A2f2Δt2(c+s)T(c-s)---(8)]]>

这里

s^T＝(sin 2πf(1+τ)Δt，sin 2πf(2+τ)Δt，Λ，sin 2πf(N+τ)Δt)(9)

W=1+τ0Λ002+τΛ0MMΛM00ΛN+τ---(10)]]>

式中：

A  ------表示谐波干扰余弦函数的振幅；

x  ------表示地震记录的时间序列向量，野外地震数据采集得到；

c   -----表示谐波干扰余弦函数向量，根据频率、时间采样率计算得到；

s   -----表示谐波干扰正弦函数向量，根据频率、时间采样率计算得到；

f   -----表示谐波干扰频率，由自适应频率计算算法确定；

τ  -----表示谐波干扰时延，由自适应时延计算算法确定；

Δt------表示地震记录时间采样间隔，野外地震数据采集得到；

T  -----表示向量转置或者矩阵转置；

W  -----表示系数矩阵，由谐波干扰时延确定；

N  ------表示地震记录时间采样长度，野外地震数据采集得到；

6)根据目标函数一阶导数和二阶导数和计算谐波干扰频率修正量Δf和时延修正量Δτ，其计算公式为：

采用以下公式计算确定谐波干扰频率修正量Δf，

Δf=1▿f2Q▿fQ---(11)]]>

采用以下公式计算确定谐波干扰时延修正量Δτ，

Δτ=1▿τ2Q▿τQ---(12)]]>

7)根据谐波干扰初始频率f₀、初始时延τ₀、频率修正量Δf和时延修正量Δτ，计算确定谐波干扰频率f和时延τ，采用以下频率修正公式确定谐波干扰频率f，

f＝f₀-γΔf    (13)

采用以下时延修正公式确定谐波干扰时延τ，

τ＝τ₀-μΔτ(14)

其中f₀是初始频率，γ是迭代频率自适应步长，通过实验确定频率自适应步长；τ₀是初始时延，在计算时初始时延值可以设置为0；μ是迭代时延自适应步长，通过实验确定频率自适应步长；初始频率f₀由步骤2)确定；计算第一步中，f为f₀，τ为τ₀；

8)根据谐波干扰频率f、谐波干扰时延τ、谐波干扰余弦函数的振幅A、地震记录的时间采样间隔Δt；

9)计算地震有效信号；

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