[发明专利]贝叶斯时移地震差异反演方法及装置

摘要

本发明提供了一种贝叶斯时移地震差异反演方法及装置。方法包括：对采集的地震数据进行匹配处理；利用匹配良好的地震数据提取角度依赖的地震子波，确定振幅缩放因子；提取得到多个弹性参数及对应的均值，计算包含多个弹性参数之间统计相关性的协方差矩阵；提取得到注采引起的弹性参数变化量曲线及对应的均值，计算包含多个弹性参数变化量统计相关性的协方差矩阵；根据沉积模式建立时间域的基础数据对应的初始弹性参数模型和差异地震数据对应的初始弹性参数变化量模型；估算正演算子矩阵；构建反演目标函数，对其求导并令导数为零，得到弹性参数变化量的求解表达式；对求解表达式进行迭代求解，获得最优弹性参数变化量反演结果。

主权项

1.一种贝叶斯时移地震差异反演方法，其特征在于，所述方法包括：对同一工区不同时间采集的地震数据进行匹配处理，将匹配良好的基础数据和监测数据相减，得到时移地震差异数据；利用匹配良好的地震数据提取角度依赖的地震子波，结合测井数据、精确Zoeppritz方程以及实际井旁道观测地震数据，确定振幅缩放因子；根据工区内基础数据对应的测井数据，提取得到多个弹性参数及对应的均值，并计算包含多个弹性参数之间统计相关性的协方差矩阵；利用测井手段或地球物理方法获取监测数据对应井位的实际测井曲线，根据工区内的测井数据提取得到注采引起的弹性参数变化量曲线及对应的均值，并计算包含多个弹性参数变化量统计相关性的协方差矩阵；利用地震构造解释资料和测井数据，根据沉积模式建立时间域的基础数据对应的初始弹性参数模型和差异地震数据对应的初始弹性参数变化量模型；推导基于精确Zoeppritz方程的时移地震差异数据正演方程，并利用基于Zoeppritz方程的反演方法估算所述时移地震差异数据正演方程的正演算子矩阵；根据Bayesian理论，同时引入包含垂向块约束项的先验模型，构建反演目标函数，对所述反演目标函数关于弹性参数变化量求导并令导数为零，得到弹性参数变化量的求解表达式；采用迭代重加权最小二乘算法对弹性参数变化量的求解表达式进行迭代求解，并通过反演残差控制最大迭代次数，获得最优弹性参数变化量反演结果。