1.一种油气储层参数标签数据生成方法,其特征在于,包括:
根据研究区的测井数据,利用高斯混合分布函数,得到研究区的岩相的先验分布概率;
根据研究区的岩相的先验分布概率,进行马尔科夫链蒙特卡洛随机模拟生成多条随深度变化的岩相曲线;
根据生成的多条随深度变化的岩相曲线,进行序贯高斯模拟,生成多条随横向变化的岩相曲线;
根据多条随深度变化的岩相曲线、多条随横向变化的岩相曲线和研究区的油气储层参数测井曲线取值的先验特征,随机填充各个岩相,得到研究区的物性参数曲线;
根据研究区的物性参数曲线,进行岩相约束统计岩石物理建模,确定研究区的弹性参数曲线;
根据研究区的弹性参数曲线,利用zoeppritz反射方程与地震子波进行褶积,生成油气储层参数标签数据。
2.如权利要求1所述的方法,其特征在于,根据研究区的测井数据,利用高斯混合分布函数,得到研究区的岩相的先验分布概率,包括:
采用期望最大估算法,对研究区的物性参数测井曲线和弹性参数测井曲线,进行多元高斯混合分布函数估算,得到物性参数对应的高斯混合分布函数中的比例系数和弹性参数对应的高斯混合分布函数中的比例系数;
根据物性参数对应的高斯混合分布函数中的比例系数,确定研究区的基于物性划分的岩相的种类与先验分布概率;
根据弹性参数对应的高斯混合分布函数中的比例系数,确定研究区的基于弹性划分的岩相的种类与先验分布概率;
根据研究区的基于物性划分的岩相的种类与先验分布概率和研究区的基于弹性划分的岩相的种类与先验分布概率,得到研究区的岩相的先验分布概率。
3.如权利要求1所述的方法,其特征在于,根据生成的多条随深度变化的岩相曲线,进行序贯高斯模拟,生成多条随横向变化的岩相曲线,包括:
将每条随深度变化的岩相曲线,作为种子样本,利用变差函数进行序贯高斯模拟,生成多条随横向变化的岩相曲线;
其中,所述变差函数是利用研究区的测井处的岩相样本进行椭圆函数拟合得到的。
4.如权利要求1所述的方法,其特征在于,根据多条随深度变化的岩相曲线、多条随横向变化的岩相曲线和研究区的油气储层参数测井曲线取值的先验特征,随机填充各个岩相,得到研究区的物性参数曲线,包括:
根据研究区的油气储层参数测井曲线取值的先验特征,构建高斯混合分布函数,随机产生符合各个岩相层段的油气储层参数先验分布特征的多个储层参数;
根据所述多个储层参数、多条随深度变化的岩相曲线和多条随横向变化的岩相曲线,得到研究区的物性参数曲线。
5.如权利要求1所述的方法,其特征在于,按照如下公式,根据研究区的物性参数曲线,进行岩相约束统计岩石物理建模,确定研究区的弹性参数曲线:
e=f(r,f)+χ
其中,e表示弹性参数,r表示物性参数,f表示岩相;
f(·)表示岩石物理模型,由研究区的岩石物理实验或经验关系得到;
χ表示岩石物理模型与实际值间的统计误差,服从高斯截断分布,由实测的测井曲线与模拟的测井曲线间的误差分布特征统计得到。
6.如权利要求1所述的方法,其特征在于,根据研究区的弹性参数曲线,利用zoeppritz反射方程与地震子波进行褶积,生成油气储层参数标签数据,包括:
根据研究区的弹性参数曲线,利用zoeppritz反射方程与地震子波进行褶积,引入不同强度的随机噪声,合成叠前地震角度道集;其中,所述叠前地震角度道集包含多个不同噪声强度下振幅随入射角变化的叠前道集的样本;
提取所述样本的岩相、物性参数和弹性参数,作为油气储层参数标签数据。
7.一种油气储层参数标签数据生成装置,其特征在于,包括:
先验分布计算模块,用于根据研究区的测井数据,利用高斯混合分布函数,得到研究区的岩相的先验分布概率;
蒙特卡洛模拟模块,用于根据研究区的岩相的先验分布概率,进行马尔科夫链蒙特卡洛随机模拟生成多条随深度变化的岩相曲线;
序贯高斯模拟模块,用于根据生成的多条随深度变化的岩相曲线,进行序贯高斯模拟,生成多条随横向变化的岩相曲线;
物性填充模块,用于根据多条随深度变化的岩相曲线、多条随横向变化的岩相曲线和研究区的油气储层参数测井曲线取值的先验特征,随机填充各个岩相,得到研究区的物性参数曲线;
统计岩石物理建模模块,用于根据研究区的物性参数曲线,进行岩相约束统计岩石物理建模,确定研究区的弹性参数曲线;
油气储层参数标签数据生成模块,用于根据研究区的弹性参数曲线,利用zoeppritz反射方程与地震子波进行褶积,生成油气储层参数标签数据。
8.如权利要求7所述的装置,其特征在于,所述先验分布计算模块,包括:
比例系数确定单元,用于采用期望最大估算法,对研究区的物性参数测井曲线和弹性参数测井曲线,进行多元高斯混合分布函数估算,得到物性参数对应的高斯混合分布函数中的比例系数和弹性参数对应的高斯混合分布函数中的比例系数;
物性划分单元,用于根据物性参数对应的高斯混合分布函数中的比例系数,确定研究区的基于物性划分的岩相的种类与先验分布概率;
弹性划分单元,用于根据弹性参数对应的高斯混合分布函数中的比例系数,确定研究区的基于弹性划分的岩相的种类与先验分布概率;
先验分布确定单元,用于根据研究区的基于物性划分的岩相的种类与先验分布概率和研究区的基于弹性划分的岩相的种类与先验分布概率,得到研究区的岩相的先验分布概率。
9.如权利要求7所述的装置,其特征在于,所述序贯高斯模拟模块具体用于:
将每条随深度变化的岩相曲线,作为种子样本,利用变差函数进行序贯高斯模拟,生成多条随横向变化的岩相曲线;
其中,所述变差函数是利用研究区的测井处的岩相样本进行椭圆函数拟合得到的。
10.如权利要求7所述的装置,其特征在于,所述物性填充模块包括:
储层参数确定单元,用于根据研究区的油气储层参数测井曲线取值的先验特征,构建高斯混合分布函数,随机产生符合各个岩相层段的油气储层参数先验分布特征的多个储层参数;
物性参数曲线确定单元,用于根据所述多个储层参数、多条随深度变化的岩相曲线和多条随横向变化的岩相曲线,得到研究区的物性参数曲线。
11.如权利要求7所述的装置,其特征在于,所述统计岩石物理建模模块具体用于:
按照如下公式,根据研究区的物性参数曲线,进行岩相约束统计岩石物理建模,确定研究区的弹性参数曲线:
e=f(r,f)+χ
其中,e表示弹性参数,r表示物性参数,f表示岩相;
f(·)表示岩石物理模型,由研究区的岩石物理实验或经验关系得到;
χ表示岩石物理模型与实际值间的统计误差,服从高斯截断分布,由实测的测井曲线与模拟的测井曲线间的误差分布特征统计得到。
12.如权利要求7所述的装置,其特征在于,所述油气储层参数标签数据生成模块具体用于:
根据研究区的弹性参数曲线,利用zoeppritz反射方程与地震子波进行褶积,引入不同强度的随机噪声,合成叠前地震角度道集;其中,所述叠前地震角度道集包含多个不同噪声强度下振幅随入射角变化的叠前道集的样本;
提取所述样本的岩相、物性参数和弹性参数,作为油气储层参数标签数据。
13.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至6任一所述方法。
14.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至6任一所述方法的计算机程序。