[发明专利]地震资料信噪比提高方法、计算机存储介质及系统

权利要求书

1.一种地震资料信噪比提高方法，其特征在于，包括以下步骤：

获得源数据和目标数据的多组对应道数据点；

针对每一组对应道数据点，确定该组对应道数据点与预设的最小二乘法拟合表达式的最小值对应的滤波器目标函数；

根据所述滤波器目标函数对该组对应道数据点的源数据点进行子波匹配，获得与该源数据点子波匹配后的源数据点；

将所述多组对应道数据点获得的多个子波匹配后的源数据点进行排序，获得子波匹配后的源数据，以消除子波不一致性所带来的差异，提高地震资料的信噪比。

2.根据权利要求1所述的地震资料信噪比提高方法，其特征在于：还包括，构建所述最小二乘法拟合表达式，其包括：

对所述源数据进行滤波处理得到滤波后的源数据；

对所述目标数据进行反滤波处理得到反滤波后的目标数据；

将所述目标数据与滤波后的源数据两者之差的平方作为第一表达式；

将所述源数据与反滤波后的目标数据两者之差的平方作为第二表达式；

将所述第一表达式与所述第二表达式之和构建为所述最小二乘法拟合表达式。

3.根据权利要求2所述的地震资料信噪比提高方法，其特征在于：所述最小二乘法拟合表达式满足：

E＝(F·X-Y)²+(X-F^-1·Y)²；

其中，F为频域中的滤波器目标函数，X为源数据，Y为目标数据。

4.根据权利要求3所述的地震资料信噪比提高方法，其特征在于：所述滤波器目标函数对应的滤波器为常相位滤波器。

5.根据权利要求3或4所述的地震资料信噪比提高方法，其特征在于：所述滤波器目标函数满足：

F＝A₀·exp(iP₀)；

其中，P₀为相位系数，A₀为振幅系数，i为虚数单位。

6.根据权利要求5所述的地震资料信噪比提高方法，其特征在于：通过第一预设表达式计算所述相位系数，其中，所述第一预设表达式为：

P₀＝arctan(-hilbert(xcor(Y,X))/xcor(Y,X))；

xcor(Y,X)是零延迟时刻所估计的时移后的源数据和目标数据的互相关函数，hilbert(xcor(Y,X))是零延迟时刻所估计的时移后的源数据和目标数据的互相关的希尔伯特变换函数。

7.根据权利要求5所述的地震资料信噪比提高方法，其特征在于：通过第二预设表达式计算所述振幅系数，其中，所述第二预设表达式为：

acor(X)是零延迟时刻的源数据自相关函数，acor(Y)是零延迟时刻的目标数据自相关函数。

8.根据权利要求1所述的地震资料信噪比提高方法，其特征在于：通过将所述滤波器目标函数与源数据点作褶积实现子波匹配。

9.一种计算机存储介质，其特征在于：所述计算机存储介质中存储有计算机程序，当所述计算机程序被一个或多个处理器执行时实现如权利要求1-8任一项所述的地震资料信噪比提高方法。

10.一种计算机系统，其特征在于：包括处理器和存储器，所述存储器中存储有计算机程序，所述计算机程序被所述处理器执行时实现如权利要求1-8任一项所述的地震资料信噪比提高方法。