[发明专利]高精度时变子波提取方法和系统

权利要求书

1.一种高精度时变子波提取方法，其特征在于，包括：

获取目标区域的目标地震数据；所述目标地震数据携带所述目标区域的地下空间地质信息；

对所述目标地震数据进行时频变换，得到时频地震数据；

利用对数变换对所述时频地震数据进行分解，得到时变子波时频谱和反射系数时频谱；

对所述时变子波时频谱进行傅里叶级数拟合，得到拟合之后的时变子波时频谱；

基于所述时频地震数据和所述拟合之后的时变子波时频谱，得到反射系数；

所述时频变换为预设广义S变换；对所述目标地震数据进行时频变换，得到时频地震数据，包括：

通过算式：对所述目标地震数据进行预设广义S变换，得到时频地震数据；

其中，τ为时间，f为频率，x(t)为所述目标地震数据，λ、b和p为三个预设常数，V_sx(τ,f)为所述时频地震数据；

利用对数变换对所述时频地震数据进行分解，得到时变子波时频谱和反射系数时频谱，包括：

将所述时频地震数据分解为时变子波时频谱与反射系数时频谱之积；

利用对数变换，将对数形式的时频地震数据分解为对数形式的时变子波时频谱与对数形式的反射系数时频谱之和；

基于所述时频地震数据和所述拟合之后的时变子波时频谱，得到反射系数，包括：

基于所述时频地震数据，所述拟合之后的时变子波时频谱，所述时频地震数据和所述时变子波时频谱之间的目标关系式，得到目标反射系数时频谱；所述目标关系式为：ln|V_sx(τ,f)|≈ln|σ(τ,f)|+ln|V_sr(τ,f)|，V_sx(τ,f)为所述时频地震数据，σ(τ,f)为所述时变子波时频谱，V_sr(τ,f)为所述反射系数时频谱；

对所述目标反射系数时频谱进行时频变换反变换，得到反射系数。

2.根据权利要求1所述的方法，其特征在于，对所述时变子波时频谱进行傅里叶级数拟合，得到拟合之后的时变子波时频谱，包括：

建立通过最小二乘法对所述时变子波时频谱进行傅里叶级数拟合时的目标函数；其中，所述目标函数的数学表达式为：f(τ,f_m)为傅里叶级数展开，为L2范数的平方，为拟合误差，m为频率点，f_m为第m个频率点的频率；

通过信赖域算法，对所述目标函数进行迭代求解，得到拟合之后的时变子波时频谱。

3.一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至2任一项所述的方法的步骤。

4.一种具有处理器可执行的非易失的程序代码的计算机可读介质，其特征在于，所述程序代码使所述处理器执行所述权利要求1-2任一项所述方法。