专利名称: |
岩石孔隙孔径的测定方法及设备 |
摘要: |
本发明提供一种岩石孔隙孔径的测定方法及设备,该方法包括:获取待测定岩石对应的三维岩心图像,并提取三维岩心图像中的所有孔隙,获取预设的几何体的尺寸,分别根据几何体的尺寸计算各孔隙最多可填充的几何体的数目,得到各孔隙对应的最多填充数目,分别根据各孔隙对应的最多填充数目和几何体的尺寸计算各孔隙对应的孔隙体积,分别对各孔隙对应的孔隙体积进行等体积计算,得到各孔隙的孔径值。本发明提供的岩石孔隙孔径的测定方法可以提高孔径计算的准确度。 |
专利类型: |
发明专利 |
国家地区组织代码: |
北京;11 |
申请人: |
中国科学院地质与地球物理研究所 |
发明人: |
李国梁;杨继进 |
专利状态: |
有效 |
申请日期: |
2019-05-14T00:00:00+0800 |
发布日期: |
2019-07-16T00:00:00+0800 |
申请号: |
CN201910397472.8 |
公开号: |
CN110018108A |
代理机构: |
北京同立钧成知识产权代理有限公司 |
代理人: |
朱颖;刘芳 |
分类号: |
G01N15/08(2006.01);G;G01;G01N;G01N15 |
申请人地址: |
100029 北京市朝阳区北土城西路19号中科院地质与地球物理研究所 |
主权项: |
1.一种岩石孔隙孔径的测定方法,其特征在于,包括: 获取待测定岩石对应的三维岩心图像,并提取所述三维岩心图像中的所有孔隙; 获取预设的几何体的尺寸; 分别根据所述几何体的尺寸计算各孔隙最多可填充的几何体的数目,得到各孔隙对应的最多填充数目; 分别根据各孔隙对应的最多填充数目和所述几何体的尺寸计算各孔隙对应的孔隙体积; 分别对各孔隙对应的孔隙体积进行等体积计算,得到各孔隙的孔径值。 2.根据权利要求1所述的方法,其特征在于,在所述提取所述三维岩心图像中的所有孔隙之前,还包括: 对所述三维岩心图像进行图像阈值分割得到所述三维岩心图像中的所有孔隙区域和矿物区域; 相应地,所述提取所述三维岩心图像中的所有孔隙,包括: 分别根据预设的分水岭算法对各孔隙区域进行分割,得到各孔隙区域中的所有孔隙。 3.根据权利要求1所述的方法,其特征在于,在所述分别对各孔隙对应的孔隙体积进行等体积计算,得到各孔隙的孔径值之后,还包括: 根据各孔隙的孔径值得到所述待测定岩石中的孔隙的孔径分布。 4.根据权利要求3所述的方法,其特征在于,所述根据各孔隙的孔径值得到所述待测定岩石中的孔隙的孔径分布,包括: 获取多个指定孔径值; 针对各个指定孔径值,统计孔径值为该指定孔径值的孔隙的个数; 在预设孔径数目坐标系中,对各个指定孔径值和各个指定孔径值对应的个数进行曲线拟合,生成孔径数目与孔径的关系曲线。 5.根据权利要求1所述的方法,其特征在于,所述分别根据所述几何体的尺寸计算各孔隙最多可填充的几何体的数目,得到各孔隙对应的最多填充数目,包括: 分别根据各孔隙包含的体素的坐标信息确定各孔隙的尺寸; 分别根据各孔隙的尺寸和所述几何体的尺寸计算各孔隙最多可填充的几何体的数目,得到各孔隙对应的最多填充数目。 6.根据权利要求2所述的方法,其特征在于,所述分别根据预设的分水岭算法对各孔隙区域进行分割,得到各孔隙区域中的所有孔隙之前,还包括: 获取各孔隙区域包含的体素数目; 筛选出体素数目大于预设体素数目的孔隙区域。 7.根据权利要求1-6任一项所述的方法,其特征在于,所述分别对各孔隙对应的孔隙体积进行等体积计算,得到各孔隙的孔径值,包括: 分别根据各孔隙对应的孔隙体积和预设的球体体积公式计算各孔隙对应的球体直径; 将各孔隙对应的球体直径作为对应的孔隙的孔径值。 8.一种岩石孔隙孔径的测定设备,其特征在于,包括: 孔隙提取模块,用于获取待测定岩石对应的三维岩心图像,并提取所述三维岩心图像中的所有孔隙; 尺寸获取模块,用于获取预设的几何体的尺寸; 填充数目计算模块,用于分别根据所述几何体的尺寸计算各孔隙最多可填充的几何体的数目,得到各孔隙对应的最多填充数目; 孔隙体积计算模块,用于分别根据各孔隙对应的最多填充数目和所述几何体的尺寸计算各孔隙对应的孔隙体积; 孔径值计算模块,用于分别对各孔隙对应的孔隙体积进行等体积计算,得到各孔隙的孔径值。 9.一种岩石孔隙孔径的测定设备,其特征在于,包括:至少一个处理器和存储器; 所述存储器存储计算机执行指令; 所述至少一个处理器执行所述存储器存储的计算机执行指令,使得所述至少一个处理器执行如权利要求1至7任一项所述的岩石孔隙孔径的测定方法。 10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质中存储有计算机执行指令,当处理器执行所述计算机执行指令时,实现如权利要求1至7任一项所述的岩石孔隙孔径的测定方法。 |
所属类别: |
发明专利 |