[发明专利]一种确定地层流体饱和度的方法、装置及设备

权利要求书

1.一种确定地层流体饱和度的方法，其特征在于，包括：

获取不同深度点的核磁共振谱组成的第一混合矩阵和第二混合矩阵；其中，所述不同深度点的核磁共振谱是基于目标地层不同深度点的回波数据反演获得；所述第二混合矩阵中的数据未进行线性转换；

对所述第一混合矩阵和所述第二混合矩阵进行分析，基于分析结果确定所述目标地层中流体类型的数量；其中，所述对所述第一混合矩阵和所述第二混合矩阵进行分析，基于分析结果确定所述目标地层中流体类型的数量，包括：利用小波PCA算法对所述第一混合矩阵进行分析，获得第一分析结果；其中，所述第一分析结果表示流体类型个数与贡献量的第一关系；利用2DPCA算法对所述第二混合矩阵进行分析，获得第二分析结果；其中，所述第二分析结果表示流体类型个数与贡献量的第二关系；根据所述第一分析结果和所述第二分析结果，确定所述目标地层中流体类型的数量；

基于所述流体类型的数量，利用几何学解混方法对所述第一混合矩阵进行处理，获得不同流体类型在核磁共振谱上的伪特征，并组合为第三混合矩阵；

利用非负矩阵分解算法对所述第三混合矩阵进行分解，获得不同流体类型在核磁共振谱上的特征；

基于所述第一混合矩阵和利用非负矩阵分解算法得到的不同流体类型在核磁共振谱上的特征，利用预设方式确定不同深度点中不同流体类型的饱和度；其中，所述预设方式满足指定约束条件。

2.根据权利要求1所述的方法，其特征在于，所述利用小波PCA算法对所述第一混合矩阵进行分析，获得第一分析结果，包括：

对所述第一混合矩阵进行小波分解，获得小波近似系数和小波细节系数；

对所述小波细节系数进行软阈值处理，获得处理后的小波细节系数；

基于所述小波近似系数和所述处理后的小波细节系数，重构所述第一混合矩阵对应的小波处理数据；

对所述小波处理数据进行主成分分析，获得第一分析结果。

3.根据权利要求1所述的方法，其特征在于，所述根据所述第一分析结果和所述第二分析结果，确定所述目标地层中流体类型的数量，包括：

根据所述第一分析结果和所述第二分析结果中相同流体类型个数对应的贡献量，确定流体类型个数与贡献量差值的关系；

基于流体类型个数与贡献量差值的关系，确定贡献量差值满足预设条件时所对应的流体类型个数；

将所述贡献量差值满足预设条件时所对应的流体类型个数作为所述目标地层中流体类型的数量。

4.根据权利要求1所述的方法，其特征在于，所述基于所述流体类型的数量，利用几何学解混方法对所述第一混合矩阵进行处理，获得不同流体类型在核磁共振谱上的伪特征，并组合为第三混合矩阵，包括：

基于所述流体类型的数量，利用不同几何学解混方法对所述第一混合矩阵进行处理，获取不同流体在核磁共振谱上的伪特征；其中，不同流体在核磁共振谱上的伪特征为几何学解混方法得到的端元矩阵；

对不同几何学解混方法对应的端元矩阵中数据进行硬阈值处理；

将硬阈值处理后的端元矩阵进行组合，获得第三混合矩阵。

5.根据权利要求4所述的方法，其特征在于，所述几何学解混方法至少包括下述之一：像元纯度指数法、顶点成分分析法、正交子空间投影法、N-FINDR法、单形体生长算法、正交基算法。

6.根据权利要求1所述的方法，其特征在于，所述预设方式为反演方法，所述反演方法包括最小二乘法；所述指定约束条件包括饱和度非负约束以及饱和度总和为1的约束。