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

说明书

本发明公开了一种地震分辨率提高方法及电子设备，涉及地震勘探技术领域。地震分辨率提高方法包括：在频率域中，首先，利用有效频段内的地震数据频谱特征和预设伪反射系数，基于自动相关判别先验信息和最大期望算法，得到谱约束伪反射系数频谱，并反推地震数据中主频段外的信息，重点恢复地震数据中的低频段信息；然后对谱约束伪反射系数频谱进行差异性滤波，获得准确的可信频谱范围；最后对可信频谱范围进行反傅里叶变换，以得到拓频处理后的地震数据，从而提高了地震分辨率。本发明提出的方法，解决了现有技术中对地震数据低频拓展能力低，拓频中无法兼顾信噪比和分辨率的同时提升，拓频地震数据不准确且噪音高的问题。

技术领域

本发明涉及地震勘探技术领域，具体而言，涉及一种地震分辨率提高方法及电子设备。

背景技术

地震数据的分辨率直接影响后续地震解释、地震反演的质量和效果，对油气地震勘探具有重要的意义。而地震数据的频段宽度和质量是决定地震分辨率的重要因素。一般假设地震数据的带限特征是由于子波、吸收衰减等因素引起的，通过拓频处理可以消除子波影响，有效地恢复宽频反射系数的频段。

经发明人研究发现，现有的拓频处理方法可以分为：反褶积类拓频处理方法、谱白化类拓频处理方法和基于稀疏反演的拓频处理方法。反褶积类拓频处理方法重点在于拓展地震数据的高频成分，存在低频拓展能力低的问题。谱白化类拓频处理方法在应用的时候，主频段中的有效信号质量会获得提升，主频段外的高频和低频段信号能量会获得补偿，但低频和高频段地震数据信噪比较低，导致在有效信号补偿的同时，频段中的噪音也将同样被放大，存在难以兼顾信噪比和分辨率的同时提升的问题。基于稀疏约束反演的拓频处理方法，当地震数据中存在薄层时，稀疏约束不成立，难以获得准确的反射系数反演结果，导致后续褶积获得的地震数据也不准确；并且此方法在时间域中应用，运用了地震数据的全频段信息，然而地震数据主频范围外的地震数据品质较差，此类方法会在补偿的频段中引入放大的噪音，存在地震数据不准确且噪音高的问题。因此，提供一种能够对地震数据进行低频拓展，同时去除噪音的方法是亟待解决的技术问题。

发明内容

针对现有技术中存在的低频拓展能力低，拓频中无法兼顾信噪比和分辨率的同时提升，拓展地震数据不准确且噪音高的问题，本发明提出了一种地震分辨率提高方法，以解决上述问题，达到准确拓展地震数据，同时提升信噪比的目的。

本发明提出了一种地震分辨率提高方法，该方法包括：

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

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

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

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

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

在本发明实施例较佳的选择中，在上述地震分辨率提高方法中，所述地震数据为叠后地震数据。