[发明专利]一种地震数据分解方法、系统及存储介质与终端

权利要求书

1.一种地震数据分解方法，其特征是：包括以下步骤：

步骤一、确定地震波形中绝对值最大的振幅值，该振幅值对应的反射时间即为零相位雷克子波的中心时间；

步骤二、以步骤一的最大绝对值振幅对应的反射时间为中心点，确定雷克子波的主周期；

步骤三、根据步骤二的主周期计算雷克子波的主频与幅度值；

步骤四、根据中心时间、主频和幅度值计算零相位雷克子波波形；

步骤五、从地震波形中减去步骤四所得的子波分量，得到残差信号；将残差信号作为输入信号，转至步骤一继续循环迭代，直到残差信号的能量小于给定的门槛值为止。

2.如权利要求1所述的一种地震数据分解方法，其特征是：所述步骤一中，对地震波形逐点求绝对值，确定出其中绝对值最大的振幅值，将其对应的反射时间设为零相位雷克子波的中心时间。

3.如权利要求1所述的一种地震数据分解方法，其特征是：所述步骤二中，以步骤一的最大绝对值振幅对应的反射时间为中心点，向其两侧寻找地震波形中的相邻极值点，确定雷克子波的主周期。

4.如权利要求1所述的一种地震数据分解方法，其特征是：所述步骤二中具体过程包括：

(1)若最大绝对值振幅为波峰，则寻找两侧相邻的波谷，此两个波谷之间的时间跨度为雷克子波的主周期；

(2)若最大绝对值振幅为波谷，则寻找两侧相邻的波峰，此两个波峰之间的时间跨度为雷克子波的主周期。

5.如权利要求1所述的一种地震数据分解方法，其特征是：所述步骤三中具体过程为对主周期取倒数，对主周期倒数值乘以得到雷克子波的主频f_m。

6.如权利要求1所述的一种地震数据分解方法，其特征是：所述步骤三中根据地震波形与雷克子波之间的相关系数计算雷克子波的幅度，计算雷克子波的幅度的具体过程包括：

(1)计算地震波形与雷克子波的互相关系数；

(2)计算雷克子波本身的自相关系数；

(3)计算互相关系数与自相关系数的比值，该比值为雷克子波的幅度。

7.如权利要求1所述的一种地震数据分解方法，其特征是：所述步骤四中，由中心时间、主频和幅度值得到零相位雷克子波波形；具体步骤为：将中心时间、主频、幅度代入雷克子波方程中，计算雷克子波的波形：

其中，A为雷克子波的幅度，f_m为雷克子波的主频，t₀为中心时间。

8.一种地震数据分解系统，其特征是：运行于处理器或存储器上，被配置为执行以下指令：

步骤一、确定地震波形中绝对值最大的振幅值，该振幅值对应的反射时间即为零相位雷克子波的中心时间；

步骤二、以步骤一的最大绝对值振幅对应的反射时间为中心点，确定雷克子波的主周期；

步骤三、根据步骤二的主周期计算雷克子波的主频与幅度值；

步骤四、根据中心时间、主频和幅度值计算零相位雷克子波波形；

步骤五、从地震波形中减去步骤四所得的子波分量，得到残差信号；将残差信号作为输入信号，转至步骤一继续循环迭代，直到残差信号的能量小于给定的门槛值为止。

9.一种计算机可读存储介质，其中存储有多条指令，其特征是：所述指令适于由终端设备的处理器加载并执行权利要求1-7中任一项所述的地震数据分解方法。

10.一种终端设备，其特征是：包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令适于由处理器加载并执行权利要求1-7中任一项所述的地震数据分解方法。