[发明专利]孔隙度预测方法及装置

权利要求书

1.一种孔隙度预测方法，其特征在于，包括：

根据待分析区的叠前地震数据，确定多角度叠加的地震数据；

根据测井点的岩石岩性的剪切模量参数、体积模量参数和密度参数，确定测井点的岩石岩性的剪切模量低频参数、体积模量低频参数和密度低频参数；

根据所述多角度叠加的地震数据，以及所述测井点的岩石岩性的剪切模量低频参数、体积模量低频参数和密度低频参数，基于Gray近似的贝叶斯线性反演方法，确定待分析区的岩石岩性的剪切模量反演参数；

根据测井点的粘土含量参数，确定测井点的岩石基质的剪切模量参数；

根据所述测井点的岩石基质的剪切模量参数，确定待分析区的岩石基质的剪切模量参数；

根据所述待分析区的岩石岩性的剪切模量反演参数，以及所述待分析区的岩石基质的剪切模量参数，基于临界孔隙度模型，预测待分析区的孔隙度。

2.如权利要求1所述的方法，其特征在于，根据测井点的岩石岩性的剪切模量参数、体积模量参数和密度参数，确定测井点的岩石岩性的剪切模量低频参数、体积模量低频参数和密度低频参数，包括；

对测井点的岩石岩性的剪切模量曲线、体积模量曲线和密度曲线进行插值处理；

对插值处理后的测井点的岩石岩性的剪切模量曲线、体积模量曲线和密度曲线进行低通滤波，确定测井点的岩石岩性的剪切模量低频参数、体积模量低频参数和密度低频参数。

3.如权利要求1所述的方法，其特征在于，根据所述多角度叠加的地震数据，以及所述测井点的岩石岩性的剪切模量低频参数、体积模量低频参数和密度低频参数，基于Gray近似的贝叶斯线性反演方法，确定待分析区的岩石岩性的剪切模量反演参数，包括：

建立待分析区的岩石岩性的剪切模量反演参数的多变量高斯分布函数；

根据所述多角度叠加的地震数据，确定多角度叠加的地震数据的多变量高斯分布函数；

根据所述多角度叠加的地震数据的多变量高斯分布函数，以及所述测井点的岩石岩性的剪切模量低频参数、体积模量低频参数和密度低频参数，基于Gray近似的贝叶斯线性反演方法，求解所述岩石岩性的剪切模量反演参数的多变量高斯分布函数，确定待分析区的岩石岩性的剪切模量的反演参数。

4.如权利要求3所述的方法，其特征在于，按照如下方式确定待分析区的岩石岩性的剪切模量的反演参数，包括：

其中，为待分析区的岩石岩性的反演参数矩阵，其中，K为体积模量的反演参数，μ为剪切模量的反演参数，ρ为密度的反演参数；为测井点的岩石岩性的低频参数矩阵，其中，为体积模量的低频参数，为剪切模量的低频参数，为密度的低频参数；Σ_m为待分析区的岩石岩性的反演参数的协方差矩阵；为多角度叠加的地震数据，其中，d(θ₁),d(θ₂),...,d(θ_M)为不同角度的地震数据，M为角度数目；Σ_d为多角度叠加的地震数据的协方差，其中，为噪声的能量，I为单位矩阵；G为基于Gray近似的AVO正演矩阵。

5.如权利要求1所述的方法，其特征在于，根据测井点的粘土含量参数，确定测井点的岩石基质的剪切模量参数，包括：

根据测井点的粘土含量曲线，基于V-R-H平均模型，确定测井点的岩石基质的剪切模量曲线。

6.如权利要求5所述的方法，其特征在于，按照如下方式确定所述测井点的岩石基质的剪切模量曲线，包括：

其中，μ_m为测井点的岩石基质的剪切模量；V_clay为粘土含量；μ_clay、μ_qua分别为粘土和石英矿物的剪切模量。

7.如权利要求1所述的方法，其特征在于，根据所述测井点的岩石基质的剪切模量参数，确定待分析区的岩石基质的剪切模量参数，包括：

对测井点的岩石基质的剪切模量曲线进行滤波处理和深时转换；

对滤波处理和深时转换之后的测井点的岩石基质的剪切模量曲线，沿层位进行横向插值，确定待分析区的岩石基质的剪切模量参数。