[发明专利]基于贝叶斯分类的储层参数预测方法及系统

权利要求书

1.一种基于贝叶斯分类的储层参数预测方法，其特征在于，该方法包括：

1)基于已知的储层弹性参数，建立储层物性参数第一反演方程；

2)基于贝叶斯公式，将所述第一反演方程转换为关于所述储层弹性参数的似然函数的第二反演方程；

3)基于贝叶斯-序贯高斯方法，求解所述储层弹性参数的似然函数；

4)基于所得到的所述储层弹性参数的似然函数值求解第二反演方程获得所述储层物性参数值，

其中，所述储层参数包括所述储层弹性参数和所述储层物性参数。

2.根据权利要求1所述的基于贝叶斯分类的储层参数预测方法，其中，

所述储层物性参数包括孔隙度φ，含水饱和度sw和泥质含量C；

所述储层弹性参数包括纵波速度V_p，横波速度V_s和密度ρ。

3.根据权利要求2所述的基于贝叶斯分类的储层参数预测方法，其中，反演所述储层物性参数是孔隙度φ、含水饱和度sw和泥质含量C中的任意一种或任意组合。

4.根据权利要求2所述的基于贝叶斯分类的储层参数预测方法，其中，所述第一反演方程为：

其中，[φ,sw,C]为储层物性参数，φ为孔隙度，sw为含水饱和度，为泥质含量；为储层弹性参数，V_p为纵波速度，V_s为横波速度，ρ为密度；j为储层物性参数的类别数。

5.根据权利要求4所述的基于贝叶斯分类的储层参数预测方法，其中，所述第二反演方程为：

其中，[φ,sw,C]为储层物性参数，φ为孔隙度，sw为含水饱和度，C为泥质含量；为储层弹性参数，V_p为纵波速度，V_s为横波速度，ρ为密度；j为储层物性参数的类别数；P([[φ,sw,C]_j)为储层物性参数先验分布；P(V_p|[φ,sw,C]_j)为纵波速度的似然函数，P(V_s|[φ,sw,C]_j)为横波速度的似然函数，P(ρ|[φ,sw,C]_j)为密度的似然函数。

6.根据权利要求1所述的基于贝叶斯分类的储层参数预测方法，其中，基于贝叶斯-序贯高斯方法，求解所述储层弹性参数的似然函数包括：

3.1)创建三维储层模型，将原始的井数据和包含待模拟网格单元的柱体内的储层参数均值变换成正态分布的数据；

3.2)指定一个访问待模拟网格单元的随机路径；

3.3)从局部分布函数p(x_o|x_s,z₁,z₂)中随机抽样得到x_o的模拟结果，其中，x_o表示待模拟的值，x_s表示邻近的单元中已经模拟的值，z₁和z₂是在包含x_o所在的柱体内的第一种储层参数和第二种储层参数；

3.4)将x_o作为已经模拟的值，根据步骤3.2)指定的随机路径继续模拟其他网格单元，直到访问完所有的网格单元；

3.5)对步骤3.4)模拟结果进行正态反变换，获得所述储层弹性参数的似然函数值。