[发明专利]低频补偿方法及低频补偿装置

说明书

本发明提供了低频补偿方法及低频补偿装置，低频补偿方法包括如下步骤：S10：设计低频补偿因子，低频补偿因子为将频率域中各频率点所对应的能量补偿倍数；S20：将地震数据由第一时间域转换到频率域，与低频补偿因子相乘后，由频率域转换至第二时间域，以得到低频补偿后的第二时间域；S30：利用二元全区间插值方法计算出地层特征开始发生变化的起始时间值，将起始时间值后的地震数据与S20中的地震数据进行替换，以得到补偿后的地震数据，以解决现有技术中的地震数据中深层信噪比技术会使高频信号损失的问题。

技术领域

本发明涉及石油地震勘探技术领域，具体而言，涉及一种低频补偿方法及低频补偿装置。

背景技术

地震数据在采集过程中接收有效信号的同时，也会接收到大量噪声，有效信号能量会随着地震波传播距离的增强而变弱，而噪声能量尤其是高频噪声能量不会发生变化。这样就会造成地震数据中深层信噪比偏低。而大部分地震数据中深层地层结构复杂，反射波场杂乱，信噪比低，难以建立精细的速度场，偏移成像困难。

目前，解决上述技术问题主要通过低通滤波将地震数据高频噪声滤除，提高地震数据中深层信噪比，通过设计一个滤波因子，该因子在低频端系数为A＝1，高频端系数为0，低频端到高频端过渡带系数从1逐渐过渡到0；滤波因子通过快速傅里叶正变换从频率域变换到时间域，然后与地震数据进行褶积，得到低通滤波后的地震数据。

但是，现有的方法虽然滤掉了高频噪声，但同时也滤除了高频有效信号，不利于地震数据的提高分辨率处理。

发明内容

本发明的主要目的在于提供一种低频补偿方法及低频补偿装置，以解决现有技术中的地震数据中深层信噪比技术会使高频信号损失的问题。

为了实现上述目的，根据本发明的一个方面，提供了一种低频补偿方法，包括如下步骤：S10：设计低频补偿因子，低频补偿因子为将频率域中各频率点所对应的能量补偿倍数；S20：将地震数据由第一时间域转换到频率域，与低频补偿因子相乘后，由频率域转换至第二时间域，以得到低频补偿后的第二时间域；S30：利用二元全区间插值方法计算出地层特征开始发生变化的起始时间值，将起始时间值后的地震数据与S20中的地震数据进行替换，以得到补偿后的地震数据。

进一步地，低频补偿方法还包括如下步骤：S40：将S30中的起始时间值与S20中的第二时间域数据进行合并，以得到补偿后的地震数据。

进一步地，将起始时间值与第二时间域数据进行合并的公式为：

其中，

D(t)为合并后的数据；

B(t)为补偿前的时间域数据；

A(t)为补偿后的时间域数据；

T为时间起始点，L为常数。

进一步地，在频率域中，将频率f划分为三个范围，分别为0＜f≤f1，f1＜f≤f2，f＞f2；其中，f1为低频补偿起始频率，f2为低频补偿截止频率，f2大于f1。

进一步地，当0＜f≤f1时，低频补偿因子为A，其中，A为常数，A大于1.0。

进一步地，当f1＜f≤f2时，低频补偿因子为

进一步地，当f＞f2时，低频补偿因子为1.0。

进一步地，在S20步骤中：在地震数据内取n个采样点，第一时间域的地震数据为实数数组，变换到频率域后为n个采样点的复数数组，低频补偿因子为与地震数据频率域点数相同的实数数组，低频补偿因子各点频率大小与频率域地震数据相同；将频率域地震数据复数数组与低频补偿因子实数数组相乘后，得到补偿后的频率域数据；将补偿后的频率域数据通过快速傅里叶逆变换转换回第二时间域，以得到低频补偿后的第二时间域的地震数据。