[发明专利]三维观测系统噪声压制估算方法

权利要求书

1.一种三维地震采集观测系统压制噪音估算方法，包括以下步骤：

(1)收集用地震勘探设备获取的地震剖面、井资料、表层的地震波速度和频率、地下目的层的平均速度、地面次生干扰、地理信息；

(2)利用获取的地震剖面、井资料、表层的地震波速度和频率、地下目的层的平均速度、地面次生干扰、地理信息，结合地质任务设计多种三维地震勘探观测系统方案；

(3)计算CMP属性信息，其包括CMP面元位置与干扰源的距离、CMP道集内每个炮检对的偏移距、方位角以及炮检连线与CMP面元至次生干扰源连线的夹角，所述的夹角以下简称为“干扰源的相对方位角”；

(4)通过下面公式(1)计算的有效波与噪音之间的剩余时差：

δtij=12Vd(xij2+4li2+4xijlicosαij+xij2+4li2-4xijlicosαij)-t01+xij2V2t02+Δtij---(1)]]>

其中：i为下标，表示CMP的序号

j为下标，表示CMP中地震道的序号

δt_ij为第i个CMP道集中第j道的噪音的剩余时差

V_d为噪音的传播速度

x_ij为第i个CMP道集中第j道炮检距

l_i为第i个CMP道集与次生干扰源距离，通过次生干扰源位置坐标和CMP点位置坐标计算得到

α_ij为第i个CMP道集中第j道的干扰源相对方位角

t₀为双程旅行时表示的目的层埋深

V为反射波的传播速度

Δt_ij为低降速带对第i个CMP道集中第j道的时间影响；

(5)根据剩余时差计算噪音的叠加振幅特性；

(6)通过计算各个共中心点噪音的综合振幅压制特性，

以及通过公式计算各个共中心点噪音的分频振幅压制特性；

其中：k为下标，表示噪音频率采样的序号；为第i个CMP点的噪音的综合振幅特性；M为噪音频率的采样个数；K_i(ω_k)为第i个CMP点第k个频率采样的噪音叠加振幅特性，ω_k为噪音频率范围内的第k个频率值；为第i个CMP点噪音的分频平均振幅特性，m₁和m₂是计算噪音压制特性时限定的频率范围，且以频率采样的序号表示，即1≤m₁＜m₂≤M；

(7)根据综合振幅压制特性频数和分频振幅压制特性频数，得到相应的频数图；从得到的频数图导出频幅特性的平面变化和统计特征；

(8)用振幅压制特性最小的三维观测系统实施地震数据采集，然后将采集的资料去做资料处理，获到较高信噪比的地震剖面。

2.如权利要求1所述的三维地震采集观测系统压制噪音估算方法，其特征在于步骤(5)所述的剩余时差计算噪音的叠加振幅特性是计算三维地震数据采集时伴生的折射波、直达波和次生干扰的剩余时差，由此得到它们的叠加振幅特性。