[发明专利]基于双参数目标寻优的自适应谱模拟反褶积方法

权利要求书

1.基于双参数目标寻优的自适应谱模拟反褶积方法，其特征在于，该基于双参数目标寻优的自适应谱模拟反褶积方法包括：

步骤1，采集地震数据，并依据数据的时频特性分不同时窗计算地震记录振幅谱；

步骤2，依次选取不同谱模拟参数对该地震记录振幅谱进行拟合，拟合出一簇地震子波振幅谱；

步骤3，由该簇地震子波振幅谱，计算地震子波自相关；

步骤4，由自相关法从该地震记录振幅谱中估算出地震子波自相关；

步骤5，分别对谱模拟法估算出的地震子波自相关与自相关法估算出的地震子波自相关做归一化处理； 

步骤6，分别计算归一化后的谱模拟法估算出的地震子波自相关与自相关法估算出的地震子波自相关间的误差，误差能量最小者所对应的子波振幅谱即为估算出的最优子波振幅谱；以及

步骤7，依据该最优子波振幅谱，设计反褶积算子，对地震数据进行反褶积。

2.根据权利要求1所述的基于双参数目标寻优的自适应谱模拟反褶积方法，其特征在于，在步骤2中，假设子波振幅谱是光滑的，给定谱模拟参数范围，依次选取不同谱模拟参数对该地震记录振幅谱进行最小二乘拟合，拟合出该簇地震子波振幅谱。

3.根据权利要求2所述的基于双参数目标寻优的自适应谱模拟反褶积方法，其特征在于，在步骤2中，进行最小二乘拟合时选取的数学模型为：

式中，为频率，是频率的阶多项式。

4.根据权利要求3所述的基于双参数目标寻优的自适应谱模拟反褶积方法，其特征在于，在步骤2中，参数的遍历范围是,参数的遍历范围是。

5.根据权利要求1所述的基于双参数目标寻优的自适应谱模拟反褶积方法，其特征在于，在步骤2中，当该地震数据为叠加后数据时，对拟合残差视为一非平稳信号做经验模态分解，提取其中包含的趋势项，作为残留的子波谱参与到子波振幅谱的估计中。

6.根据权利要求1所述的基于双参数目标寻优的自适应谱模拟反褶积方法，其特征在于，在步骤3中，根据自相关能量的收敛情况，确定子波自相关长度。

7.根据权利要求6所述的基于双参数目标寻优的自适应谱模拟反褶积方法，其特征在于，在步骤4中，依据该子波自相关长度，由自相关法从该地震记录振幅谱中估算出地震子波自相关，并对结果加斜坡或加窗。

8.根据权利要求1所述的基于双参数目标寻优的自适应谱模拟反褶积方法，其特征在于，在步骤5中，对该误差能量采用递减型的权函数加权求和，以选取该误差能量最小者。