[发明专利]一种地震记录自动增益方法

说明书

本发明公开了一种地震记录自动增益方法，包括以下步骤：确定时窗，包括滑动时窗长度；获取地震道数据，计算地震道数据中各样点的绝对值；确定地震道数据始端前半时窗长度中各样点的增益值；确定地震道数据中部各样点的增益值；确定地震道数据末端后半时窗长度中各样点的增益值；获取目标地震道数据。本发明能够解决地震资料处理中由于能量衰减，导致振幅失真的问题。本发明方法能够进行振幅增益，实现振幅的自动恢复，使生成的地震道数据能量更加均衡，提高地震资料处理目标的精度，从而能够有效恢复地震信号，提高地震资料的分辨率和信噪比。

技术领域

本发明涉及一种地震记录自动增益方法，属于石油地震勘探技术领域。

背景技术

随着技术的进步，地震资料处理目标向精细化方向发展，精度要求越来越高，针对岩性，流体识别，孔隙度，饱和度的计算，特别是页岩气勘探的发展，地震资料处理对于地震资料的振幅保真要求越来越高。然而，在实际地震资料采集中由于大地吸收作用，使得单炮记录远道能量较近道能量弱，深层反射能量较浅层弱，或者由于上覆盖火成岩影响，能量难以透射，造成火成岩下方反射波极弱，因此，在地震资料处理阶段需要正确恢复振幅值，并在后续处理流程中加以保护。虽然，目前常规的恢复振幅的方法有很多，包括基于球面扩散的振幅补偿，地表一致性振幅补偿，道内振幅均衡，Q补偿等，但是，常规的恢复振幅的方法依然不能满足地震资料处理对于地震资料的振幅保真需求。

基于此，亟需提出一种地震记录自动增益方法，以解决远道能量较近道能量弱，深层反射能量较浅层弱，从而造成的数据能量较弱的问题。

发明内容

针对上述问题，本发明的目的是提供一种地震记录自动增益方法，其能够解决地震资料处理中由于能量衰减，导致振幅失真的问题。

为实现上述目的，本发明采取以下技术方案：一种地震记录自动增益方法，包括以下步骤：确定时窗，包括滑动时窗长度；获取地震道数据，计算地震道数据中各样点的绝对值；确定地震道数据始端前半时窗长度中各样点的增益值；确定地震道数据中部各样点的增益值；确定地震道数据末端后半时窗长度中各样点的增益值；获取目标地震道数据。

在一个具体实施例中，根据地震道数据始端前半时窗长度中各样点的增益值，地震道数据中部各样点的增益值，以及地震道数据末端后半时窗长度中各样点的增益值，获取振幅增益的目标地震道数据。

在一个具体实施例中，确定时窗的长度为W，前半时窗和后半时窗的长度分别为L＝W/2。

在一个具体实施例中，确定地震道数据的总样点数为N，根据地震道数据中各样点的绝对值，确定地震道数据中对应样点的绝对振幅平均值。

在一个具体实施例中，确定地震道数据始端前半时窗长度中各样点的样点序号为(0，1，2，......，L)，根据对应的各个样点序号，计算地震道数据始端前半时窗长度中L个样点的绝对振幅平均值。

在一个具体实施例中，确定地震道数据始端前半时窗长度中各样点的增益值所依据的公式为：

agcArray[i]＝ampArray[i]/QArray[i] (1)

式中，agcArray[i]为地震道数据中第i个样点的增益值，i依次取0，1，2，......，L，ampArray[i]为地震道数据中第i个样点的原始值，QArray[i]为地震道数据中前i+L个样点的绝对振幅平均值。

在一个具体实施例中，确定地震道数据中部各样点的样点序号为(L+1， L+2，......，N-L)，根据对应的各个样点序号，通过时窗沿地震道数据由前至后滑动，计算时窗位置确定时中心点的前L个样点和后L个样点的绝对振幅平均值。

在一个具体实施例中，确定地震道数据中部各样点的增益值所依据的公式为：

agcArray[i′]＝ampArray[i′]/QArray[i′] (2)