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

摘要

本发明涉及石油地球物理勘探技术，是在地震数据采集施三维观测系统噪声压制估算方法，步骤为：设计多种三维地震勘探观测系统方案，计算CMP属性信息，计算的有效波与噪音之间的剩余时差，根据剩余时差计算噪音的叠加振幅特性；计算各个共中心点噪音的综合振幅压制特性和分频振幅压制特性，得到频数图统计特征；用振幅压制特性最小的三维观测系统实施地震数据采集，获到较高信噪比的地震剖面。本发明实现了有效地压制了噪声，提高信噪比，改善三维观测系统的数据采集质量，获得较高信噪比的地震剖面，从而提高实施钻探的准确度。

主权项

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道炮检距，参考图1      l_i为第i个CMP道集与次生干扰源距离，通过次生干扰源位置坐标和CMP点位置坐标计算得到，参考图1      α_ij为第i个CMP道集中第j道的干扰源相对方位角      t₀为双程旅行时表示的目的层埋深V为反射波的传播速度Δt_ij为低降速带对第i个CMP道集中第j道的时间影响；(5)根据剩余时差计算噪音的叠加振幅特性；(6)通过Ai‾=1MΣk=1MKi(ωk)]]>计算各个共中心点噪音的综合振幅压制特性，以及通过公式Awi‾=1(m2-m1+1)Σk=m1m2Ki(ωk)]]>计算各个共中心点噪音的分频振幅压制特性；其中：k为下标，表示噪音频率采样的序号；A_i为第i个CMP点的噪音的综合振幅特性；M为噪音频率的采样个数；K_i(ω_k)为第i个CMP点第k个频率采样的噪音叠加振幅特性，ω_k为噪音频率范围内的第k个频率值；Aw_i为第i个CMP点噪音的分频平均振幅特性，m₁和m₂是计算噪音压制特性时限定的频率范围，且以频率采样的序号表示，即1≤m₁＜m₂≤M；(7)根据综合振幅压制特性频数和分频振幅压制特性频数，得到相应的频数图；从得到的频数图导出频幅特性的平面变化和统计特征；(8)用振幅压制特性最小的三维观测系统实施地震数据采集，然后将采集的资料去做资料处理，获到较高信噪比的地震剖面。