[发明专利]一种利用自适应分析时窗提取地震属性的方法、存储介质

权利要求书

1.一种利用自适应分析时窗提取地震属性的方法，其特征在于，包括：

步骤S100，获取目的层的地震信号的平均瞬时频率；

步骤S200，根据平均瞬时频率求取分析时窗的长度以及确定分析时窗的顶端和底面的原始层位面；

步骤S300，对分析时窗的顶端和底面的原始层位面分别进行平滑处理，获得顶端和底面的优化层位面；其中，所述步骤S300包括以下步骤：

步骤S310，对于分析时窗的顶端/底面的原始层位面中的任意一点(x₀,y₀,t₀)，按照下式计算表征该点的曲面能量的六个分量：

基于上述六个分量形成对称的二阶张量，记为：

令：以建立该点的特征方程组，

其中，λ为求解所述特征方程组获得的特征值，λ为多解；

取上述特征值中最大的值λ_max作为该点的特征值；

步骤S320，当获得分析时窗的顶端/底面的原始层位面中所有点的特征值后，建立如下趋势拟合目标函数F：

F＝∑λ_max+∑(H-H′)

其中，H为分析时窗的顶端/底面的原始层位面，H′为顶端/底面的优化层位面，

步骤S330，令趋势拟合目标函数F的值达到全局最小，由此确定顶端和底面的优化层位面H′；

步骤S400，保持分析时窗的长度不变，结合顶端和底面的优化层位面构建新的分析时窗；

步骤S500，利用新的分析时窗提取目的层的地震属性，以进行地震属性分析。

2.根据权利要求1所述的利用自适应分析时窗提取地震属性的方法，其特征在于，所述步骤S100包括以下步骤：

步骤S110，将目的层的地震信号转换为离散地震信号；

步骤S120，利用小波变换对目的层的离散地震信号进行变换；

步骤S130，基于下式求取变换后的目的层的离散地震信号的平均瞬时频率：

其中，FT(f,t)为变换后的目的层的离散地震信号，f为瞬时频率，t为瞬时时间点，f_m为平均瞬时频率。

3.根据权利要求2所述的利用自适应分析时窗提取地震属性的方法，其特征在于，所述步骤S200中，根据下式确定分析时窗的长度：

其中，n为根据分析时窗的长度需求而选取的常数值，T为周期，L为分析时窗的长度。

4.根据权利要求3所述的利用自适应分析时窗提取地震属性的方法，其特征在于，n取值范围为0.1至2。

5.根据权利要求3所述的利用自适应分析时窗提取地震属性的方法，其特征在于，所述步骤S200中，分析时窗的顶端和底面的原始层位面分别为：

其中，H_top和H_bot分别为分析时窗的顶端和底面的原始层位面，H₁为目的层的层位面。

6.一种计算机可读存储介质，其中存储有系统程序，所述系统程序在被处理器执行时实现如权利要求1至5中任意一项所述的利用自适应分析时窗提取地震属性的方法。