[发明专利]三维地震数据面波衰减方法及装置

说明书

技术领域

本发明是关于地球物理勘探数据处理技术，特别是关于一种三维地震数据面波衰减方法及装置。

背景技术

在现有的物探地震数据处理技术中，目前使用的面波衰减方法主要有两种：1)采用高通滤波的方法；2)根据面波的局部线性特征进行面波衰减的方法。

对于采用高通滤波进行面波衰减的方法，主要由于面波与有效信号之间有较大的频率差异，面波的频率偏低，有效信号频率偏高，衰减面波的同时有效信号与面波频率重叠的部分也会被当成面波衰减，不适合低频勘探。

对于根据面波的局部线性特征进行面波衰减，可以通过倾角能量扫描方式确定面波的倾角，这种方法可以在面波与有效信号重叠的频率部分更好地保护有效信号，更适合低频勘探数据处理。

但是第2)种方法，只能处理局线性分布的面波，地震勘探野外数据采集技术的发展主要经历了一维地震勘探(勘探早期)、二维地震勘探(勘探中期)、三维地震勘探(目前主要勘探方式)三个阶段。不同的地震勘探，面波的分布也有变化。一维和二维地震勘探，面波在近炮检和远炮检距都是局部线性分布，而三维地震勘探远排列面波分布就成抛物线分布。抛物线分布的面波如果采用线性倾角扫描的方法去识别就比较困难，面波衰减效果也在很大程度上受到影响。

发明内容

本发明提供一种三维地震数据面波衰减方法及装置，以提高三维地震勘探面波衰减的性能。

为了实现上述目的，本发明提供一种三维地震数据面波衰减方法，所述的三维地震数据面波衰减方法包括：

步骤1：采集三维地震数据；

步骤2：选择所述三维地震数据中的一个炮记录，根据所述炮记录的炮点的坐标与检波点的坐标计算所述炮点到所述检波点的距离；

步骤3：根据所述距离从所述炮记录中选择一个与所述炮点最近的排列，并在所述排列上选择一段信噪比最高的初至波同相轴；

步骤4：根据该段初至波的地震道参数生成初至波的视速度V；

步骤5：将所述炮点到所述检波点的连线投影到一直线上，该直线在坐标系中经过所述炮点并平行于检波线；

步骤6：根据所述视速度V、所述炮点、所述检波点及所述检波点在所述直线上的投影点的坐标生成投影时间差T；

步骤7：根据三维地震数据采样率及所述投影时间差计算所述炮记录中每个排列的地震道的样点个数Lj，其中，j为地震道序号，j＝1,2,3,4,5……N；

步骤8：将所述炮记录中每个地震道的所有样点向时间减少的方向移动Lj个样点；

步骤9：对所述三维地震数据中所有炮记录执行所述步骤2至步骤8，生成第一地震数据；

步骤10：根据面波的局部线性分布特征对所述第一地震数据进行面波衰减，生成第二地震数据；

步骤11：将所述衰减后的三维地震数据中每个地震道的所有样点都向时间增加的方向移动Lj个样点，生成第三地震数据。

在一实施例中，在所述步骤4中，根据该段初至波的地震道参数生成初至波的视速度V，包括：

计算该段初至波最后一个地震道所在的炮检距与第一个地震道所在的炮检距的差，作为初至波的传播距离；

计算该段初至波最后一个地震道所在的时间与第一个地震道所在的时间的差，作为初至波的传播时间；

利用所述传播距离除以所述传播时间，生成所述初至波的视速度V。

在一实施例中，在所述步骤6中，根据所述视速度V、所述炮点、所述检波点及所述检波点在所述直线上的投影点的坐标生成投影时间差T，包括：

计算炮检距与投影偏移距的差值S；

根据所述视速度V及所述差值S计算所述投影时间差T，T＝S/V；

其中，所述炮检距为所述炮点与检波点之间的距离，所述投影偏移距为所述炮点与所述投影点之间的距离。

在一实施例中，所述的方法还包括：把每个地震道的样点个数Lj记录到地震道的道头中。

在一实施例中，在所述步骤10中，根据面波的局部线性分布特征对所述第一地震数据进行面波衰减，生成第二地震数据包括：根据面波的局部线性分布特征，采用倾角能量扫描方式确定面波的位置，对所述第一地震数据进行面波衰减，生成所述第二地震数据。

为了实现上述目的，本发明提供一种三维地震数据面波衰减装置，所述的三维地震数据面波衰减装置包括：

数据采集单元，用于采集三维地震数据；

距离计算单元，用于选择所述三维地震数据中的一个炮记录，根据所述炮记录的炮点的坐标与检波点的坐标计算所述炮点到所述检波点的距离；