[发明专利]一种基于岩石微CT图像的孔隙迂曲度计算方法

权利要求书

1.一种基于岩石微CT图像的孔隙迂曲度计算方法，其特征在于，该方法包括：

(1)以岩心微CT图像为基础，在图像REV分析的基础上，通过图像像素与网格单元一一对应的关系构建用于扩散模拟的网格单元模型；

(2)利用构建的岩心REV网格单元模型，基于随机游走理论采用MATLAB工具箱开展孔隙扩散模拟，计算有效扩散系数；

(3)结合岩石样品本征扩散系数与模拟得到的有效扩散系数分析，可以计算得到岩心的三维迂曲度。

2.根据权利要求1所述的一种基于岩石微CT图像的孔隙迂曲度计算方法，其特征在于，所述步骤(1)中，以岩心微CT图像为基础，通过降噪、滤波、分割等图像处理技术，提取包含岩心骨架和孔隙信息的二值化图像。

3.根据权利要求1所述的一种基于岩石微CT图像的孔隙迂曲度计算方法，其特征在于，所述步骤(1)中，基于获得的岩心二值化图像，以孔隙度和渗透率等参数为参考依据，计算样品的REV尺寸。

4.根据权利要求1所述的一种基于岩石微CT图像的孔隙迂曲度计算方法，其特征在于，所述步骤(1)中，基于图像像素与网格单元在空间一一对应的关系，通过像素替换方法重建样品的REV三维计算网格模型。

5.根据权利要求1所述的一种基于岩石微CT图像的孔隙迂曲度计算方法，其特征在于，所述步骤(2)中，将重建得到的岩心REV网格模型的骨架和孔隙相网格数据矩阵导入到MATLAB中，基于随机游走理论开展孔隙扩散模拟。

6.根据权利要求1所述的一种基于岩石微CT图像的孔隙迂曲度计算方法，其特征在于，所述步骤(2)中，扩散模拟假想粒子在孔隙空间中不同时间步的重复随机运动行为，通过统计粒子群的均方位移和时间的导数，进而计算得到模型的有效扩散系数D_e。

7.根据权利要求1所述的一种基于岩石微CT图像的孔隙迂曲度计算方法，其特征在于，所述步骤(3)中，通过分析岩石的本征扩散系数D和有效扩散系数的关系，结合上述步骤中获得的粒子群均方位移和孔隙度等参数，可以计算求得岩石的三维迂曲度。