[发明专利]三维叠前炮检域随机噪声压制方法

说明书

技术领域

本发明涉及地球物理勘探地震数据处理方法，属于压制噪声的技术范畴，具体一种三维叠前炮检域随机噪声压制方法。

背景技术

在陆上地震数据采集中，由于地表复杂及地下地质结构等原因，采集的地震数据信噪比极低，特别是沙漠地区，随机噪声比较严重，常规做法是在CMP(沿测线方向)域中处理，这种方法很难满足需求。

在CMP域中的实现方法是首先将数据抽道CMP域中，一个CMP内的道数比较少，受实际覆盖次数的影响，预测方法使用的空间道数就少，预测效果就受影响，另外在CMP域预测有一个缺陷，使信噪比高的地方改进很大，但信噪比低的地方改进不大，而且后续再进行其它域中去噪效果不好。

发明内容

本发明目的是提供一种通过对观测系统的变形重排技术，去噪效果更加显著的三维叠前炮检域随机噪声压制方法。

本发明通过以下步骤实现：

1)采集地震数据，输入观测系统定义以后的数据，经过处理得到动校正之后的CMP道集数据，并保留；

步骤1)所述的处理是做常规叠加处理；所述的保留是将对动校正之后的数据记录保留，准备做三维叠前炮检域随机噪声压制。

2)将动校正之后的数据在线方向进行线号重排，在沿测线方向(CMP)按炮序重排，使重排之后每条CMP线相当于一条二维测线；

3)修改道头数据，顺序重排；

步骤3)所述的数据重排是在三维施工时，一炮一般对应多个接收排列，每个接收排列与每个平行炮线相当于一条二维线，把每个二维线按地表连续位置重排，每个排列对应一炮。

4)进行三维叠前炮检域随机噪声处理，得到全部模型道；

步骤4)所述的三维叠前炮检域随机噪声压制是首先对四维数据体按照空间上的参数为线数，炮数，道数，时间上给定一个时窗，这样就形成一个小的三维数据体，把这个三维数据体变换到相应的F-XYZ域，在F-XYZ域求取预测算子，然后对这个小数据体进行预测滤波，得到预测结果；

然后在空间上重复道数，在时间上重复样点数，得到下一个小的三维数据体，然后求取下一个三维数据体的预测结果，反复在空间和时间上滑动求取，完成对整个三维叠前数据的处理。

5)将全部模型道与步骤1)数据混波，完成噪声压制。

步骤5)所述的混波是首先将步骤4)全部模型道做叠加，与步骤1)步保留的常规叠加做对比，根据两者剖面噪比情况，确定混波比例。

所述的确定混波比例是信噪比低的部位，混入模型道比例加大。

本方法通过对观测系统的变形重排技术，可以在不同域中实现对随机噪声压制，去噪效果更加显著。

附图说明

图1是数据重排前的炮检位置图；

图2是数据重排后的炮检位置图；

图3是随机噪声压制的时窗滑动示意图；

图4是随机噪声压制前道集；

图5是随机噪声压制后道集；

图6是随机噪声压制前叠加；

图7是随机噪声压制后叠加。

具体实施方式

本发明是随机噪声衰减方法，用于三维地震采集的数据。本发明基于三维频率空间的F-XYZ域预测去噪技术。它假设地震记录中的有效波在F-XYZ域具有可预测性，而随机噪声无此特性。利用多道复数最小平方原理求取三维预测算子，并用该预测算子对该频率成分的四维地震数据体进行预测滤波，达到衰减随机噪声的目的。

本发明具体步骤如下：

第一步：采集地震数据，输入观测系统定义之后的数据，经过一系列处理，得到动校正之后的CMP道集数据。一系列处理包括：静校正，振幅保持处理，规则噪声压制处理，反褶积(可选)，动校正，剩余静校正等常规处理。最后得到动校正之后的CMP道集数据。

第二步：将动校正之后的数据在线方向进行线号重排，在沿测线方向(CMP)按炮序重排，使重排之后每条CMP线相当于一条二维测线；见附图1和2，重排前后的炮检位置图。从图中看到，数据重排之后，位置图更加规则。如果按三维显示的化，数据的最后一维是一个排列的炮集，所以叫炮检域RNA技术。这一步也叫数据准备。

第三步：进行三维叠前炮检域随机噪声压制处理，得到全部模型道，这一步也叫三维叠前炮检域随机噪声压制。

三维叠前炮检域随机噪声压制是首先对四维数据体按照空间上的参数为线数，炮数，道数，时间上给定一个时窗，这样就形成一个小的三维数据体，把这个三维数据体变换到相应的F-XYZ域，在F-XYZ域求取预测算子，然后对这个小数据体进行预测滤波，得到预测结果；