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

权利要求书

1.一种低频补偿方法，其特征在于，包括如下步骤：

S10：设计低频补偿因子，所述低频补偿因子为将频率域中各频率点所对应的能量补偿倍数；

S20：将地震数据由第一时间域转换到频率域，与所述低频补偿因子相乘后，由频率域转换至第二时间域，以得到低频补偿后的第二时间域；

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。

2.根据权利要求1所述的低频补偿方法，其特征在于，在S20步骤中：

在地震数据内取n个采样点，第一时间域的地震数据为实数数组，变换到频率域后为n个采样点的复数数组，所述低频补偿因子为与地震数据频率域点数相同的实数数组，低频补偿因子各点频率大小与频率域地震数据相同；将频率域地震数据复数数组与低频补偿因子实数数组相乘后，得到补偿后的频率域数据；

将补偿后的频率域数据通过快速傅里叶逆变换转换回第二时间域，以得到低频补偿后的第二时间域的地震数据。

3.根据权利要求1所述的低频补偿方法，其特征在于，所述起始时间值的算法为：

给定矩形域上的n’*m个结点在两个方向上的坐标分别为：

x₀＜x₁＜u＜x_n′-1与y₀＜y₁＜v＜y_m-1

相应的函数值为：

z_kj＝z(x_k,y_j),k＝0,1,…,n′-1；j＝0,1,…,m-1

在X方向与Y方向上，以插值点(u，v)为中心，前后各取4个坐标，分别为：

x_p＜x_p+1＜x_p+2＜x_p+3＜u＜x_p+4＜x_p+5＜x_p+6＜x_p+7和

y_q＜y_q+1＜y_q+2＜y_q+3＜v＜y_q+4＜y_q+5＜y_q+6＜y_q+7；

然后，利用二元插值公式：

计算插值点(u，v)处的函数近似值，所述插值点(u，v)处的函数近似值为中深层层位发生变化点的起始时间值。

4.根据权利要求1所述的低频补偿方法，其特征在于，所述S30步骤中，将地层特征开始发生变化的多个时间起始点信息，通过二元全区间插值方法插出空间上相对应的时间起始点，以将插出后的空间上的每个点的时间起始点后的地震数据与低频补偿后的地震数据进行替换。

5.一种低频补偿装置，适用于权利要求1至4中任一项所述的低频补偿方法，其特征在于，所述低频补偿装置包括：

低频补偿因子设计模块，以用于设计低频补偿因子；

低频补偿模块，与所述低频补偿因子设计模块连接，以根据低频补偿因子设计模块中设计出的低频补偿因子对地震数据进行低频补偿，以得到低频补偿后的第二时间域数据；

中深层地震数据合并模块，与所述低频补偿模块连接，以将时间起始值与低频补偿模块中得到的补偿后的第二时间域数据进行合并，以得到补偿后的地震数据。