[发明专利]一种不依赖测井资料的含油气饱和度预测方法

权利要求书

1.一种不依赖测井资料的含油气饱和度预测方法，其特征在于，包括以下步骤：

步骤1：对目标区域进行地质背景调查，以获得目标区域的岩石物理参数，并对岩石物理参数中各个参数赋予概率，得到岩石物理参数的先验概率分布；

步骤2：从具有先验概率分布的所有岩石物理参数中对每一个岩石物理参数进行多次随机抽样若干数值，得到岩石物理参数抽样值；

步骤3：预设R-H-G岩石物理经验模型和岩石密度方程，R-H-G岩石物理经验模型如公式①，岩石密度方程如公式②：

V_p＝(1-φ)²V_ps+φV_pf------①

ρ＝(1-φ)ρ_m+φρ_f------②

式中，V_p为纵波速度，φ为孔隙度，V_ps为岩石骨架速度，V_pf为孔隙流体的速度，ρ表示岩石密度，ρ_m为岩石矿物骨架密度，ρ_f为孔隙流体密度；

步骤4：根据岩石物理参数抽样值计算出纵波速度V_p和岩石密度ρ，纵波速度V_p和岩石密度ρ相乘结果得到纵波阻抗，从而得到纵波阻抗与含油气饱和度之间的后验概率分布数据；

步骤5：在后验概率分布数据中，将纵波阻抗作为自变量、含油气饱和度作为因变量进行拟合优化后，得到通过纵波阻抗计算含油气饱和度的计算公式，该计算公式作为用于预测含油气饱和度的预测公式；

获得目标区域的纵波阻抗，将纵波阻抗代入预测公式，输出结果即是预测结果。

2.根据权利要求1所述的不依赖测井资料的含油气饱和度预测方法，其特征在于，地质背景调查的内容包括沉积盆地的类型，目的层的沉积年代和演化过程，目的层的现今埋藏深度和温压场特征，目的层的沉积、构造和储层特征以及目的层的含流体类型，得到相应的地质背景调查结果，根据地质背景调查结果综合分析，得到目标区域的岩石物理参数及对应的概率。

3.根据权利要求1所述的不依赖测井资料的含油气饱和度预测方法，其特征在于，岩石物理参数的先验概率分布包括孔隙度先验概率分布、油气饱和度先验概率分布、地层温度先验概率分布、地层压力先验概率分布、石英矿物含量先验概率分布、粘土矿物含量先验概率分布。

4.根据权利要求1所述的不依赖测井资料的含油气饱和度预测方法，其特征在于，随机抽样采用马尔可夫链蒙特卡洛随机抽样。

5.根据权利要求1所述的不依赖测井资料的含油气饱和度预测方法，其特征在于，在步骤4中，将多组岩石物理参数抽样值代入公式①和②，每一次代入采用一组岩石物理参数抽样值，一组岩石物理参数抽样值是指同一次抽样得到的岩石物理参数抽样值。

6.根据权利要求1所述的不依赖测井资料的含油气饱和度预测方法，其特征在于，在步骤5之后，还包括以下步骤：

步骤6：获得地震CMP道集数据，对地震CMP道集数据进行反演，得到层析速度数据体，将层析速度数据体代入Gardner方程，计算出密度数据体，再根据密度数据体计算出纵波阻抗数据体，并将该纵波阻抗数据体作为反演所需的初始模型，该初始模型记为纵波阻抗初始模型；

步骤7：获得地震叠后振幅数据体并代入到纵波阻抗初始模型中，采用基于模型的反演方法进行反演，得到纵波阻抗反演数据体；

步骤8：将纵波阻抗反演数据体代入所述预测公式，从而计算出含油气饱和度数据体，完成预测。

7.根据权利要求6所述的不依赖测井资料的含油气饱和度预测方法，其特征在于，步骤6中，对地震CMP道集数据进行层析网格反演，得到层析速度数据体。

8.根据权利要求1所述的不依赖测井资料的含油气饱和度预测方法，其特征在于，步骤4中，所述根据岩石物理参数抽样值计算出纵波速度V_p和岩石密度ρ的具体实现包括以下步骤：

将步骤岩石物理参数中地层温度、地层压力和油气饱和度的抽样值代入Batzle-Wang公式，计算出V_pf和ρ_f；

将石英矿物含量和粘土矿物含量的抽样值代入岩石矿物骨架密度计算公式，计算出ρ_m，岩石矿物骨架密度计算公式为：

ρ_m＝V_qρ_q+V_shρ_sh

式中，V_q代表石英含量，V_sh代表粘土矿物含量，ρ_q代表石英矿物密度，ρ_sh代表粘土矿物密度；

将石英矿物含量和粘土矿物含量的抽样值代入Voigt-Reuss-Hill公式计算出岩石骨架的体积模量K；

将上述计算得到的的ρ_m和K，通过公式③计算出V_ps：

根据计算得到的V_ps、ρ_m代入公式①和②，从而计算得到纵波速度V_p和岩石密度ρ。