[发明专利]拓宽地震数据频带的方法及装置

权利要求书

1.一种拓宽地震数据频带的方法，包括以下步骤：

步骤1：对时间域单道地震数据进行傅里叶变换，得到频率域单道地震数据；

步骤2：计算所述频率域单道地震数据的振幅；

步骤3：对所述振幅进行拟合与平滑处理，获得拟合曲线，计算所述拟合曲线最大值；

步骤4：基于所述拟合曲线和拟合曲线最大值计算单道拓频滤波算子；

步骤5：在频率域，基于所述单道拓频滤波算子计算经过拓频的频率域单道地震数据；

步骤6：对经过拓频的频率域单道地震数据进行反傅里叶变换；

其中基于以下公式计算单道拓频滤波算子：

其中，f(w)表示所述单道拓频滤波算子，Amp'(w)表示所述拟合曲线，Amp_max'(w)表示所述拟合曲线最大值，ε表示拓频控制因子。

2.根据权利要求1所述的拓宽地震数据频带的方法，其中ε的取值范围在0.01至1之间。

3.根据权利要求1所述的拓宽地震数据频带的方法，其中基于以下公式计算经过拓频的频率域单道地震数据：

F₁(w)＝Fr₁(w)+iFi₁(w)

其中，

Fr₁(w)＝Fr(w)×f(w)

Fi₁(w)＝Fi(w)×f(w)

其中，F₁(w)表示经过拓频的频率域单道地震数据，f(w)表示所述单道拓频滤波算子，Fr(w)表示所述频率域单道地震数据的实部，Fi(w)表示所述频率域单道地震数据的虚部，Fr₁(w)表示经过拓频的频率域单道地震数据的实部，Fi₁(w)表示经过拓频的频率域单道地震数据的虚部。

4.一种拓宽地震数据频带的装置，包括：

傅里叶变换模块，用于对时间域单道地震数据进行傅里叶变换，得到频率域单道地震数据；

频率域计算模块，用于计算所述频率域单道地震数据的振幅；

拟合模块，用于对所述振幅进行拟合与平滑处理，获得拟合曲线，计算所述拟合曲线最大值；

滤波算子确定模块，用于基于所述拟合曲线和拟合曲线最大值计算单道拓频滤波算子；

拓频计算模块，用于在频率域，基于所述单道拓频滤波算子计算经过拓频的频率域单道地震数据；

傅里叶反变换模块，用于对经过拓频的频率域单道地震数据进行反傅里叶变换；

其中所述滤波算子确定模块基于以下公式计算单道拓频滤波算子：

其中，f(w)表示所述单道拓频滤波算子，Amp'(w)表示所述拟合曲线，Amp_max'(w)表示所述拟合曲线最大值，ε表示拓频控制因子。

5.根据权利要求4所述的拓宽地震数据频带的装置，其中ε的取值范围在0.01至1之间。

6.根据权利要求4所述的拓宽地震数据频带的装置，其中所述拓频计算模块基于以下公式计算经过拓频的频率域单道地震数据：

F₁(w)＝Fr₁(w)+iFi₁(w)

其中，

Fr₁(w)＝Fr(w)×f(w)

Fi₁(w)＝Fi(w)×f(w)

其中，F₁(w)表示经过拓频的频率域单道地震数据，f(w)表示所述单道拓频滤波算子，Fr(w)表示所述频率域单道地震数据的实部，Fi(w)表示所述频率域单道地震数据的虚部，Fr₁(w)表示经过拓频的频率域单道地震数据的实部，Fi₁(w)表示经过拓频的频率域单道地震数据的虚部。