[发明专利]一种地震随机反演方法、设备及系统

权利要求书

1.一种地震随机反演方法，其特征在于，所述方法包括：

获取地震数据和低频模型；

构建模型参数的先验信息；

利用所述地震数据和低频模型构建融入低频信息的似然函数；

根据所述先验信息以及似然函数实现后验概率的采样，以得到反演结果。

2.根据权利要求1所述的方法，其特征在于，所述构建模型参数的先验信息包括：

利用直接序贯模拟算法构建模型参数的先验信息；

再现先验信息的直方图。

3.根据权利要求2所述的方法，其特征在于，构建的似然函数为：

其中，L(m)为似然函数，m为模型参数，d为观测数据，g为正演算子，σ_e为高斯噪声方差，α为可调因子，m_LFM为低频模型。

4.根据权利要求3所述的方法，其特征在于，根据所述先验信息以及似然函数实现后验概率的采样，以得到反演结果包括：

将吉布斯采样与序贯模拟相结合形成序贯吉布斯重采样算子；

使用所述序贯吉布斯重采样算子扰动所述先验信息；

构建改进的Metropolis算法；

结合所述序贯吉布斯重采样算子和改进的Metropolis算法进行后验概率的采样，以得到反演结果。

5.根据权利要求4所述的方法，其特征在于，改进的Metropolis算子的接受准则为：

其中，m_current为初始模型，m_propose为建议模型，L(m_current)为初始模型的似然函数，L(m_propose)为建议模型的似然函数。

6.根据权利要求4所述的方法，其特征在于，根据所述先验信息以及似然函数实现后验概率的采样，以得到反演结果还包括：

选择序贯模拟的邻域范围内的数据作为条件数据进行序贯模拟以及序贯吉布斯重采样。

7.一种地震随机反演系统，其特征在于，所述系统包括：

地震数据获取模块，用于获取地震数据和低频模型；

先验信息构建模块，用于构建模型参数的先验信息；

似然函数构建模块，用于利用所述地震数据和低频模型构建融入低频信息的似然函数；

反演结果确定模块，用于根据所述先验信息以及似然函数实现后验概率的采样，以得到反演结果。

8.根据权利要求7所述的系统，其特征在于，所述先验信息构建模块包括：

先验信息确定模块，用于利用直接序贯模拟算法构建模型参数的先验信息；

直方图再现模块，用于再现先验信息的直方图。

9.根据权利要求7所述的系统，其特征在于，构建的似然函数为：

其中，L(m)为似然函数，m为模型参数，d为观测数据，g为正演算子，σ_e是高斯噪声方差，α是一个可调因子，m_LFM是低频模型。

10.根据权利要求7所述的系统，其特征在于，所述反演结果确定模块包括：

采样算子形成模块，用于将吉布斯采样与序贯模拟相结合形成序贯吉布斯重采样算子；

先验信息扰动模块，用于使用所述序贯吉布斯重采样算子扰动所述先验信息；

改进算法构建模块，用于构建改进的Metropolis算法；

后验概率采集模块，用于结合所述序贯吉布斯重采样算子和改进的Metropolis算法进行后验概率的采样，以得到反演结果。