[发明专利]一种压制陆上地震资料多次波干扰的方法

说明书

技术领域

本发明涉及石油地震勘探的数据处理技术，是一种压制陆上地震资料多次波干扰的方法。

技术背景

地震勘探是寻找和勘探石油天然气的主要方法。主要工作包括地震数据采集、处理和解释三个步骤。地震数据处理主要为地震资料解释提供资料成果，地震数据处理的内容很多，主要有去噪、反褶积、动静校正、速度分析、叠加、偏移和反演等。为了提高信噪比，去噪成为地震资料处理必不可少的工作。地震资料中包含的噪声繁多，多次波干扰就是其中一种，压制多次波也就成为资料处理中重要的处理技术。陆上资料处理压制多次波主要就是基于多次波与一次波速度差异进行压制，常用的方法有radon变换、τp变换等，通过变换使多次波与一次波在速度时间域分离，然后在速度时间域保留有效波区域，在把多次波区域的数值置为零，最后反变换到时间空间域获得压制多次波后的地震资料。

速度时间域的压制方法对速度准确度依赖性大，而且要求多次波与一次波动校正具有足够大的时差。但是对于成因复杂的多次波，特别是陆上地震资料，产生多次波的强反射层受上覆地层物性变化（速度、表层条件、地表一致性、岩性、波动学特征等）影响，使其速度在空间上随之变化，这种速度的不确定性，以及同时伴生的各种层间多次波、混合多次波等因素影响，导致其产生的多次波特征与一次波混叠在一起，在速度时间域依然很难找到明确的分界，从而影响信噪分离效果，往往在压制多次波的同时，使有效波损失严重，特别对于信噪比非常低的弱反射信号损害严重，这种损害将直接影响到压制多次波的成败，这也是目前陆上地震资料压制多次波方法很难取得理想效果的一个重要原因。

发明内容

本发明目的在于提供一种提高地震资料成果的信噪比的压制陆上地震资料多次波干扰的方法。

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

1）采集地震数据，预处理的标准记录格式的二维、三维地震数据，以及叠加或偏移成像速度数据；

所述的预处理包括：定义观测系统、建立炮检关系及设置成像面元网格的规格；计算、应用野外静校正；叠前道集去噪、振幅处理、子波处理；

所述的预处理的地震数据是按照炮点文件、偏移距或者CMP、偏移距索引顺序排序的道集数据。

所述的叠加成像速度是以经过预处理的地震数据为基础，经过CMP道集叠加速度分析和剩余静校正迭代处理，拾取的一个以目标目的层成像为目的的速度，该速度具有一定CMP空间间隔的离散的文本速度数据，基于最终基准面或基于CMP基准面。

所述的偏移成像速度是以叠加成像速度乘以70%-100%的加权比例作为初始速度进行叠前偏移速度分析迭代，获得的以目标目的层偏移成像为目的的速度，该速度为CRP控制点离散速度和偏移速度体形式。

所述的偏移速度体通过CRP控制点离散速度经过内插、平滑处理方法获得，该速度体具有与地震数据采样间隔对应的地震数据格式的均方根速度体，基于最终基准面或基于CMP基准面。

2）将预处理的地震数据道集按照CMP、偏移距索引顺序排序，将道集按照偏移距顺序从小到大划分成一个以上的区带，对各区带地震数据按照步骤1）的成像速度分别进行叠加和偏移成像处理获得成像剖面，以成像剖面的波组出现多次波特征变化作为分辨多次波与一次波优势区带分界点的依据；

迭代以上过程，进一步细分区带，直到可量化确定一次波与多次波的优势区带分界点，从零偏移距到该分界点之间的区带即为多次波主要作用的优势区带范围；

所述区带是指地震数据的偏移距范围，最小区带宽度不小于500米。

所述的优势区带是指区带成像剖面中目标波组典型特征表现为主要特征的区带。

所述的区带的分界点为不等间隔。

所述的多次波特征是：多次波在成像剖面上表现为与其上覆产生多次波的强反射层具有相同的形态和波组特征，而且在出现的时间上表现为倍数或者近似倍数关系，特征能量强则表示该区带为多次波发育区，反之则反。

3）在多次波优势区带以外的区域确定一次波优势区带的范围；

所述的多次波优势区带和一次波优势区带接触的分界是200-400米过渡范围；

4）在产生多次波的强反射层时间段和多次波出现时间段选取相同长度的时窗进行互相关，然后计算每道零相位振幅值，根据振幅值大小进行加权，振幅大的地震道将其样点振幅乘以小的权值，震幅小的地震道将其样点振幅乘以大的权值，实现压制陆上地震资料多次波干扰；

所述的选取相同时窗长度为200ms-500ms。

所述的计算是指在给定相同长度时窗内对各道样点振幅进行的均方根求和。

所述的权值范围0-1之间。