[发明专利]一种叠前地震贝叶斯反演方法和装置

说明书

本发明公开了一种叠前地震贝叶斯反演方法和装置。该方法利用地震数据中的时间域叠前地震数据和所述频率域叠前地震数据进行正演算推倒，得到第一正演方程和第二正演方程。进一步，以贝叶斯反演模型为基础，利用第一正演方程和第二正演方程，建立时频域联合叠前地震反演目标函数，预先设定的参数约束条件的约束下，确定时频域联合叠前地震反演目标函数的最优解，得到地层弹性参数。本发明充分利用叠前地震数据中的时间域及频率域数据，使得反演结果的分辨率比现有方法有较大提高。

技术领域

本发明涉及地震监测领域，更具体的说是涉及一种叠前地震贝叶斯反演方法和装置。

背景技术

地震反演是获取地下介质弹性参数的有效途径。根据所采用地震资料的不同，地震反演可分为叠后反演和叠前反演。

在叠前地震反演中，贝叶斯反演在考虑模型及数据参数不完备性方面更有优势。贝叶斯反演理论在地球物理反演中应用广泛，该方法分别假设模型参数和似然函数服从某种先验分布，通过求解后验概率密度函数实现反演估计。

现有的叠前地震贝叶斯反演的实现主要在时间域进行，反演得到地层介弹性参数，其反演结果的分辨率低。

发明内容

有鉴于此，本发明提供一种叠前地震贝叶斯反演方法和装置，采用时频域数据联合反演的方式，得到地层弹性参数，提高了反演结果的分辨率。

为实现上述目的，本发明提供如下技术方案：

一种叠前地震贝叶斯反演方法，包括：

获取时间域叠前地震数据和频率域叠前地震数据；

分别利用所述时间域叠前地震数据和所述频率域叠前地震数据进行正演算推倒，得到第一正演方程和第二正演方程；

以贝叶斯反演模型为基础，利用所述第一正演方程和所述第二正演方程，建立时频域联合叠前地震反演目标函数；

基于预先设定的参数约束条件，确定所述时频域联合叠前地震反演目标函数的最优解，得到地层弹性参数。

优选的，所述分别利用所述时间域叠前地震数据和所述频率域叠前地震数据进行正演算推倒，得到第一正演方程和第二正演方程，之后还包括：

利用预设算法，对所述第一正演方程和所述第二正演方程中的参数进行去相关处理。

优选的，所述获取时间域叠前地震数据和频率域叠前地震数据，包括：

获取时间域叠前地震数据；

对所述时间域叠前地震数据进行傅里叶变换，得到所述频率域叠前地震数据。

优选的，所述时间域叠前地震数据为至少三个角叠加地震数据。

一种叠前地震贝叶斯反演装置，包括：

数据采集单元，用于获取时间域叠前地震数据和频率域叠前地震数据；

正演算单元，用于分别利用所述时间域叠前地震数据和所述频率域叠前地震数据进行正演算推倒，得到第一正演方程和第二正演方程；

反演目标函数建立单元，用于以贝叶斯反演模型为基础，利用所述第一正演方程和所述第二正演方程，建立时频域联合叠前地震反演目标函数；

地层弹性参数计算单元，用于基于预先设定的参数约束条件，确定所述时频域联合叠前地震反演目标函数的最优解，得到地层弹性参数。

优选的，还包括：连接与所述正演单元和所述反演目标函数建立单元之间的参数去相关处理单元，用于利用预设算法，对所述第一正演方程和所述第二正演方程中的参数进行去相关处理。

优选的，所述数据采集单元包括：