[发明专利]一种地震分辨率提高方法及电子设备

权利要求书

1.一种地震分辨率提高方法，包括：

步骤1，根据获取到的目标地震褶积模型、目标地震褶积模型的似然函数、目标地震褶积模型的边缘似然函数以及自动相关判别先验信息，得到目标后验概率表达式，并对目标地震褶积模型的边缘似然函数进行最大化处理，得到第一待反演伪反射系数估计值；

步骤2，将所述第一待反演伪反射系数估计值代入所述目标后验概率表达式进行迭代计算，得到第二待反演伪反射系数估计值和当前迭代变量，并获得当前迭代次数，判断该当前迭代次数是否大于预设迭代阈值，以及判断当前迭代变量的值是否小于预设值；

步骤3，当所述当前迭代次数大于预设迭代阈值且所述当前迭代变量的值小于预设值时，将所述第二待反演伪反射系数估计值作为新的第一待反演伪反射系数估计值，并返回执行步骤2，直至所述当前迭代次数等于预设迭代阈值或所述当前迭代变量的值小于预设值时，才将获得的第二待反演伪反射系数估计值作为谱约束伪反射系数估计值；

步骤4，根据所述谱约束伪反射系数估计值对获取到的地震数据进行处理，得到地震数据谱约束伪反射系数频谱，利用测井数据对所述地震数据谱约束伪反射系数频谱进行差异性滤波，以确定地震数据谱约束伪反射系数可信频谱范围；

步骤5，对所述地震数据谱约束伪反射系数可信频谱进行拓频处理，以获得拓频地震数据，其中，所述拓频地震数据的频带宽度大于拓频处理前的地震数据的频带宽度，从而提高地震分辨率。

2.根据权利要求1所述的地震分辨率提高方法，其特征在于，所述地震数据为叠后地震数据。

3.根据权利要求1所述的地震分辨率提高方法，其特征在于，在执行所述步骤1之前，所述方法还包括：

接收用户输入的目标时间段，并根据该目标时间段和地震数据中的子波构建时变地震物理子波库；

根据所述时变地震物理子波库和所述目标时间段，得到地震数据时间域表达式；

对所述地震数据时间域表达式进行傅里叶变换，以得到地震数据频率域表达式；

对所述地震数据频率域表达式在地震数据的指定频段内进行离散化处理，以得到频率域地震褶积模型；

根据所述频率域地震褶积模型和预设的谱约束伪反射系数，得到所述目标地震褶积模型；

当噪音服从零均值高斯分布时，根据贝叶斯原理和所述目标地震褶积模型得到该目标地震褶积模型的似然函数和目标地震褶积模型的边缘似然函数。

4.根据权利要求3所述的地震分辨率提高方法，其特征在于，根据该目标时间段和地震数据中的子波构建时变地震物理子波库的步骤包括：

在所述目标时间段内将随时间变化的非稳态地震数据作为稳态地震数据；

对所述稳态地震数据进行开窗处理，并在选定的时窗内采用井震标定的方法对子波进行提取；

对提取后的子波进行汇总，以得到所述时变地震物理子波库。