[发明专利]储层孔隙度、含水饱和度和泥质含量参数同时反演新方法

权利要求书

1.储层孔隙度、含水饱和度和泥质含量参数同时反演新方法，其特征在于，该储层孔隙度、含水饱和度和泥质含量参数同时反演新方法包括：

步骤1，采集岩石弹性参数、测井数据和岩心数据信息；

步骤2，将岩石弹性参数与储层物性参数相结合，建立岩石物理模型；

步骤3，建立储层孔隙度、含水饱和度和泥质含量参数三参数反演的目标函数，并求解目标函数；

步骤4，通过求解目标函数，输出储层孔隙度、含水饱和度和泥质含量参数反演结果；

在步骤1中，通过叠前反演获得岩石弹性参数，叠前反演是通过分析叠前道集来研究地下岩性的地震分析方法，借助于Zoeppritz方程或其近似式，最终反演出岩石弹性参数，反演出的岩石弹性参数包括岩石的体积模量、剪切模量和密度信息；

在步骤2中，首先，将岩石弹性参数与储层物性参数相结合，建立确定性的岩石物理模型；其次，针对不同储集类型的储层，以确定性岩石物理模型为基础，通过统计信息，建立统计的岩石物理模型，统计岩石物理模型用来描述测井上没有体现的所有可能的储层条件实现；所建立的岩石物理模型如公式(1)所示：

[K,μ,ρ]＝f_RPM(φ,s_w,C)+ε (1)

其中，f_RPM表示确定性岩石物理模型，K,μ,ρ表示岩石弹性参数，包括体积模量、剪切模量、密度，φ,s_w,C表示储层物性参数，包括储层孔隙度、含水饱和度、泥质含量，ε表示实际观测数据的误差；

在步骤3中，根据贝叶斯算法建立物性参数提取的目标函数，目标函数为储层物性参数在已知弹性参数条件下的最大后验概率分布，表示为：

[φ,sw,C]＝argMax{P([φ,sw,C]_j|[K,μ,ρ])} (3)

其中，φ,s_w,C分别为储层孔隙度、含水饱和度、泥质含量，P()表示概率密度，argMax表示求最大值,K表示体积模量，μ表示剪切模量，ρ表示密度，j表示各个储层物性参数的类别数；

P([φ,sw,C]_j|[K,μ,ρ])为储层物性参数的后验分布，该公式的意义是，已知弹性参数，分别计算储层孔隙度、含水饱和度和泥质含量分属于各个储层物性参数类的后验概率分布；当后验概率取最大值时，该储层物性参数即为最终反演结果；

在步骤3中，通过遗传算法和蒙特卡罗方法相结合求解目标函数，首先生成初始化群体，并计算每个个体的适应度函数，如果结果满意，则遗传算法结束，如果结果不满意，更新初始化的群体，通过蒙特卡罗进行仿真，并再次计算每个个体的适应度函数，更新群体，如此循环，直到计算结果满意为止。

2.根据权利要求1所述的储层孔隙度、含水饱和度和泥质含量参数同时反演新方法，其特征在于，在步骤1中，测井数据包括储层孔隙度、含水饱和度和泥质含量信息；岩心数据提供测井数据的辅助信息，包括岩石类型、岩石孔隙度。