1.一种基于饱和介质精确zeoppritz方程的地震振幅反演方法,其特征在于,包括:
从地震数据提取依赖于入射角度的子波;基于测井数据和各向同性饱和介质精确zeoppritz方程正演模拟地震角道集,并结合实际井旁地震数据确定振幅缩放因子;
基于工区内所有测井数据统计模型参数的先验信息;
对测井数据进行平滑处理并沿解释层位外推,获得初始参数模型,并构建基于最大后验概率的反演目标函数;
以最小化反演目标函数为目的优化所述初始参数模型,获得最优参数模型,进而利用所述参数模型实现地震振幅的反演。
2.如权利要求1所述的一种基于饱和介质精确zeoppritz方程的地震振幅反演方法,其特征在于,所述各向同性饱和介质精确zeoppritz方程将弹性介质精确zeoppritz方程中的纵横波速度替换成饱和纵横波速度。
3.如权利要求1所述的一种基于饱和介质精确zeoppritz方程的地震振幅反演方法,其特征在于,所述基于最大后验概率的反演目标函数的构建,具体包括:基于bayesian原理、先验信息和正演算子构建最大后验概率的反演目标函数,通过求目标函数关于模型参数的导数,并求解目标函数的极小值等价于令其导数等于零,结合迭代最小二乘算法获得更新迭代的模型参数公式。
4.如权利要求1所述的一种基于饱和介质精确zeoppritz方程的地震振幅反演方法,其特征在于,所述基于工区内所有测井数据统计模型参数的先验信息,具体包括:采用柯西分布函数作为先验分布函数基于工区内所有测井数据统计各模型参数的均值,求取各参数的自相关系数和互相关系数,构建参数相关的协方差矩阵,形成符合该工区的模型参数先验分布函数。
5.如权利要求1所述的一种基于饱和介质精确zeoppritz方程的地震振幅反演方法,其特征在于,所述对测井数据进行平滑处理并沿解释层位外推,获得初始参数模型,具体包括:首先利用散点插值的方法对各个层位的数据进行插值,完成地质层位建模,然后根据地质层位进行模型参数横向插值,计算得到地下每个点上的模型参数值,完成初始新属性参数建模的任务。
6.一种基于饱和介质精确zeoppritz方程的地震振幅反演系统,其特征在于,包括:
正演算子获取单元,其用于从地震数据提取依赖于入射角度的子波;基于测井数据和各向同性饱和介质精确zeoppritz方程正演模拟地震角道集,并结合实际井旁地震数据确定振幅缩放因子;
先验信息获取单元,其用于基于工区内所有测井数据统计模型参数的先验信息;
模型构建单元,其用于对测井数据进行平滑处理并沿解释层位外推,获得初始参数模型,
目标函数构建单元,其用于构建基于最大后验概率的反演目标函数;
模型参数求解单元,其用于以最小化反演目标函数为目的优化所述初始参数模型,获得最优参数模型,进而利用所述参数模型实现地震振幅的反演。
7.如权利要求6所述的一种基于饱和介质精确zeoppritz方程的地震振幅反演系统,其特征在于,所述各向同性饱和介质精确zeoppritz方程将弹性介质精确zeoppritz方程中的纵横波速度替换成饱和纵横波速度。
8.如权利要求6所述的一种基于饱和介质精确zeoppritz方程的地震振幅反演系统,其特征在于,所述基于最大后验概率的反演目标函数的构建,具体包括:基于bayesian原理、先验信息和正演算子构建最大后验概率的反演目标函数,通过求目标函数关于模型参数的导数,并求解目标函数的极小值等价于令其导数等于零,结合迭代最小二乘算法获得更新迭代的模型参数公式。
9.一种电子设备,包括存储器、处理器及存储在存储器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1-7任一项所述的基于饱和介质精确zeoppritz方程的地震振幅反演方法。
10.一种非暂态计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如权利要求1-7任一项所述的基于饱和介质精确zeoppritz方程的地震振幅反演方法。