[发明专利]基于反射率法的多波AVO储层弹性参数反演方法及系统

权利要求书

1.一种基于反射率法的多波AVO储层弹性参数反演方法，其特征是，所述的方法包括：

采集工区内的地震叠前道集、测井数据以及实际井旁角度域地震道集；

基于所述的地震叠前道集、测井数据以及实际井旁角度域地震道集确定地震子波以及振幅缩放因子；

采集工区内测井数据的统计模型参数；

根据所述的统计模型参数确定符合该工区的模型参数先验分布函数；

采集地震构造解释资料；

根据所述的地震构造解释资料以及测井数据建立深度域的初始弹性参数模型；

根据所述深度域的初始弹性参数模型确定PP波与PS波反演残差，该步骤包括：将所述深度域的初始弹性参数模型作为输入；在频率域通过传递矩阵算法确定每一个角度的PP波与PS波总反射率矩阵，所述PP波与PS波总反射率矩阵包含反射系数与传播时间、层间多次波与透射损失；对所述PP波与PS波总反射率矩阵进行傅里叶反变换；将傅里叶反变换后的PP波与PS波总反射率矩阵与地震子波进行褶积，得到角度域的PP波与PS波角道集；采集实际角度域地震道集；将所述角度域的PP波与PS波角道集与实际角度域地震道集作差，得到PP波与PS波反演残差；

构建最大后验概率意义下的反演目标函数；

根据所述的反演目标函数以及PP波与PS波反演残差确定最优弹性参数模型。

2.根据权利要求1所述的方法，其特征是，基于所述的地震叠前道集、测井数据以及实际井旁角度域地震道集确定地震子波以及振幅缩放因子包括：

根据所述的地震叠前道集以及测井数据、采用统计方法提取依赖于入射角度的地震子波；

将所述的测井数据作为输入模型；

采用反射法正演模拟角度域的PP波与PS波道集；

将所述的PP波与PS波道集与所述的实际井旁角度域地震道集进行对比，得到振幅缩放因子；

将所述的振幅缩放因子应用于所述的地震子波。

3.根据权利要求1或2所述的方法，其特征是，根据所述的统计模型参数确定符合该工区的模型参数先验分布函数包括：

所述的统计模型参数满足三变量高斯分布；

分析所述的测井数据，得到所述统计模型参数的均值；

确定所述统计模型参数的自相关系数以及互相关系数；

根据所述的自相关系数以及互相关系数构建三参数相关的协方差矩阵；

根据所述的三参数相关的协方差矩阵形成符合该工区的模型参数先验分布函数。

4.根据权利要求3所述的方法，其特征是，根据所述的地震构造解释资料以及测井数据建立深度域的初始弹性参数模型包括：

根据所述的地震构造解释资料、基于沉积模式建立地质模型；

将所述的测井资料按照构造模式进行插值和外推，得到每条测线的时间域的初始弹性参数模型，所述的初始弹性参数模型包括纵波速度模型、横波速度模型以及密度模型；

通过时深转换关系将所述时间域的初始弹性参数模型转换为深度域的初始弹性参数模型。

5.根据权利要求4所述的方法，其特征是，构建最大后验概率意义下的反演目标函数包括：

利用所述模型参数先验分布函数对反演结果进行约束；

根据Bayesian原理、综合反演似然函数以及先验分布函数得到后验概率分布函数；

根据所述的后验概率分布函数构建反演目标函数。