[发明专利]基于多域混合优化三维曲波变换的非规则地震数据重构方法

权利要求书

1.基于多域混合优化三维曲波变换的非规则地震数据重构方法，其特征在于，所述方法包括：

步骤S1、采集非规则地震数据；

步骤S2、获取曲波系数门槛值；

步骤S3、根据改进型数据驱动阈值模型、自适应最小阈值计算模型分别确定每次迭代计算阈值、最小阈值；

步骤S4、再根据曲波系数门槛值、每次迭代计算阈值、数据特征混合使用时间域重构方法和频率域重构方法进行计算，求解迭代中的曲波系数，并对曲波系数做曲波反变换处理即可获得重构地震数据；

步骤S5、重复步骤S3-S4，至每次迭代计算阈值小于最小阈值结束迭代计算，即可获得最终的高精度非规则重构数据。

2.根据权利要求1所述的基于多域混合优化三维曲波变换的非规则地震数据重构方法，其特征在于，所述步骤S2中通过下式获得曲波系数门槛值；

式中：λ_mk为曲波系数门槛值；λ_max为最大阈值；ε为最小阈值。

3.根据权利要求1所述的基于多域混合优化三维曲波变换的非规则地震数据重构方法，其特征在于，所述改进型数据驱动阈值模型为：

式中：λ_k为每次迭代计算阈值；Idx_max为最大阈值在ρ中的位置；Idx_ε为最小阈值在ρ中的位置；ρ是稀疏系数幅值按降序排列所得序列；k为迭代次数；I_ter为总迭代次数；

所述自适应最小阈值计算模型为：

式中：M_size为参与计算的矩阵大小；j为多尺度变换的分解尺度；j’为多尺度变换的总分解尺度；median为中值函数；ε为最小阈值。

4.根据权利要求1所述的基于多域混合优化三维曲波变换的非规则地震数据重构方法，其特征在于，所述步骤S4中若每次迭代计算阈值大于曲波系数门槛值则采用频率域重构方法；若每次迭代计算阈值小于或等于曲波系数门槛值则采用时间域重构方法；若每次迭代计算阈值小于最小阈值则结束迭代计算。

5.根据权利要求4所述的基于多域混合优化三维曲波变换的非规则地震数据重构方法，其特征在于，所述时间域重构方法包括：根据每次迭代计算阈值在时间域对所述非规则地震数据通过快速迭代阈值算法进行重构迭代计算，求解迭代中的曲波系数，并对曲波系数做曲波反变换处理获得重构数据。

6.根据权利要求4所述的基于多域混合优化三维曲波变换的非规则地震数据重构方法，其特征在于，所述频率域重构方法包括以下步骤：

步骤S411、先对非规则地震数据沿时间方向做一维傅里叶变换，将非规则地震数据从时间域变换到频率域；

步骤S412、再利用傅里叶系数确定有效信号频宽范围并切割出待重构频率域数据；

步骤S413、最后根据每次迭代计算阈值在频率域对待重构频率域数据通过快速迭代阈值算法进行重构迭代计算，求解迭代中的曲波系数，并对曲波系数做曲波反变换处理获得重构数据。

7.根据权利要求6所述的基于多域混合优化三维曲波变换的非规则地震数据重构方法，其特征在于，所述步骤S412中通过下式确定有效信号频宽范围：

式中：n_t为时间采样点数；n_f为有效频率样点数；为表示向上取整；f_s为采样频率；f_c为地震信号的截止频率。

8.根据权利要求6所述的基于多域混合优化三维曲波变换的非规则地震数据重构方法，其特征在于，所述步骤S412中通过下式切割出待重构频率域数据：

f'_obs＝f_obs(n_f×n_z)

式中：n_z为道数；n_f为有效频率样点数；f_obs为傅里叶系数；f'_obs为待重构频率域数据。