[发明专利]叠后信噪比估算方法及装置

说明书

技术领域

本发明本方法涉及地球物理勘探数据处理技术，特别是关于是一种叠后信噪比估算方法及装置。

背景技术

在地震数据处理技术中，噪音常常影响叠加成像效果，噪声衰减在很多地震数据处理中都变成了常规处理流程，去噪后信噪比的提升一直是一大难题。如何对野外数据采集和室内数据处理进行地震数据信噪比宏观定性分析和数字显示，准确地分析资料处理中的各个环节质量，一直是地震勘探研究人员想要实现的目标。

目前信噪比主要是在指定时间内，用信号与噪声振幅的比值来衡量信号的强弱的方法。现有技术中，主要信噪比估算方法有：1)时间域叠加法；2)时间域奇异值分减法；3)模型道相关估算法。这三种信噪比估算方法，在地下构造倾角变化较大时，信噪比估算精度不高。时间域叠加和模型道相关主要是因为地下构造倾角变化时同相轴不能准确地同相叠加。而时间域奇异值分减法本身也不适应地下构造倾角较大变化的地震数据噪声估计。

发明内容

本发明提供了一种叠后信噪比估算方法及装置，以精确地震估算地震数据的信噪比。

为了实现上述目的，本发明提供一种叠后信噪比估算方法，所述的叠后信噪比估算方法包括：

步骤1：采集地震数据；

步骤2：对所述地震数据进行预处理，形成预处理后的地震数据；

步骤3：对所述预处理后的地震数据进行叠加，形成叠后地震数据；

步骤4：在所述叠后地震数据上选择一个时间空间窗，确定所述时间空间窗的中间地震道；

步骤5：以所述中间地震道为中心，确定所述时间空间窗中第一地震道的第一样点，以所述第一样点为中心选定N个样点长的时间窗，并计算所述时间空间窗中剩余每一地震道的第一样点；

步骤6：将所述每一地震道第一个样点作为所述时间窗的中心样点向上向下分别移动N/2个样点，确定所述时间窗的位置；

步骤7：根据所述时间窗的位置从所述时间窗中选择所述中间地震道的N个样点长的数据，以所述中间地震道的样点数据与剩余地震道的样点数据做互相关，生成所述第一样点的相关比值之和；

步骤8：重复执行步骤5至步骤7，生成所述时间空间窗内中间地震道的所有样点的最大相关值之和及所述中间地震道的最大相关值总和；

步骤9：重复执行步骤4至步骤8，生成所述叠后的每个地震道在选择的时间空间窗内的最大相关值总和。

在一实施例中，对所述地震数据进行预处理，包括：对所述地震数据进行地震数据解编、静校正、噪声衰减、振幅补偿、反褶积、速度分析及叠前偏移处理。

在一实施例中，所述叠后信噪比估算方法还包括：

判断是否对所述预处理后的地震数据进行了叠前偏移处理；

如果否，对所述预处理后的地震数据进行动校正。

在一实施例中，在所述步骤7中，以所述中间地震道的样点数据与剩余地震道的样点数据做互相关，生成所述第一样点的相关比值之和，包括：以所述中间地震道的样点数据与剩余地震道的样点数据做互相关，分别获得所述剩余地震道的最大相关值，将所述剩余地震道的最大相关值按比例相加，获得所述相关比值之和。

在一实施例中，在所述步骤8中，生成所述时间空间窗内中间地震道的所有样点的最大相关值之和及所述中间地震道的最大相关值总和，包括：

生成所述时间空间窗内中间地震道中除了第一个样点外的所有样点的最大互相关值之和；

将所述时间空间窗内中间地震道中包括第一个样点的所有样点的最大互相关值之和相加，形成所述时间空间窗内中间地震道最大互相关值的总和。

在一实施例中，所述叠后信噪比估算方法还包括：对每个地震道的所述最大相关值总和计算常用对数。

在一实施例中，在所述步骤4中，选择所述时间空间窗时，时间窗中包括一个完整的同相轴，空间窗的地震道的个数为奇数，所述空间窗按所述叠后地震数据地震道顺序号依次滚动。

为了实现上述目的，本发明提供一种叠后信噪比估算装置，所述的叠后信噪比估算装置包括：

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

数据预处理单元，用于对所述地震数据进行预处理，形成预处理后的地震数据；

数据叠加单元，用于对所述预处理后的地震数据进行叠加，形成叠后地震数据；

时间窗选择单元，用于在所述叠后地震数据上选择一个时间空间窗，确定所述时间空间窗的中间地震道；