定模块， 与轮廓图处理模块连接，用于根据轮廓图的分布方式确定孔洞和孔隙在碳酸盐岩样本上的 分布方式。
[0047] 在一个实施例中，还包括：图块划分模块，用于对于每个轮廓图，将轮廓图划分为 多个图块，在每个图块中建立与图块边缘相切的椭圆形，多个图块中椭圆的短轴方向相同； 短轴确定模块，与图块划分模块连接，将每个图块短轴的长度与该图块面积相乘，得到每个 图块的乘积，将多个图块的乘积相加，得到的和为与该轮廓图边缘相切的椭圆的短轴长度， 确定为孔洞或孔隙的短轴长度。
[0048] 在一个实施例中，还包括：长轴确定模块，用于将每个图块长轴的长度与该图块面 积相乘，得到每个图块的乘积，将多个图块的乘积相加，得到的和为与该轮廓图边缘相切的 椭圆的长轴长度；图形形态确定模块，用于根据轮廓图的长轴长度与短轴长度的大小关系， 确定孔洞或孔隙的图形形态。
[0049] 在一个实施例中，还包括：骨架图处理模块，用于根据孔洞和孔隙的轮廓图，获取 孔洞和孔隙的骨架图，其中，骨架图是由线段组成的表示孔洞或孔隙的拓扑结构的图形；位 置确定模块，与骨架图处理模块连接，用于根据骨架图在碳酸盐岩样本的CT图像上的坐 标，确定孔洞或孔隙在碳酸盐岩样本上的位置。
[0050] 在本发明实施例中，通过对碳酸盐岩样本进行CT扫描，获取碳酸盐岩样本的CT图 像，并对CT图像进行二值化处理，基于经二值化处理的图像中孔隙的像素数与颗粒组成的 基质的像素数不同，来计算孔隙的总面积和碳酸盐岩样本的CT图像面积，进而可以基于孔 隙的总面积与碳酸盐岩样本的CT图像面积确定孔隙度。由于采用碳酸盐岩样本的CT图 像来确定孔隙度，与现有技术中铸体薄片、扫描电镜和压汞等手段相比，可以避免对原来的 孔隙及其内部的粘土矿物等造成破坏，以反映真实的碳酸盐岩孔隙结构特征；避免对操作 人员技术水平的依赖，避免操作过程中使用汞等有害物质；也可以得到碳酸盐岩样本整体 的CT图像，提高CT图像的代表性；同时，基于碳酸盐岩样本的CT图像，可以定量地计算孔 隙的面积和碳酸盐岩样本的CT图像面积，进而可以准确、可靠地计算碳酸盐岩样本的孔隙 度，有利于对碳酸盐岩储层的微观结构进行可靠的分析。
[0051] 显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以 用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置 所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它 们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执 行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个 模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬 件和软件结合。
[0052] 以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技 术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的 任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种碳酸盐岩微观结构的确定方法，其特征在于，包括： 获取碳酸盐岩样本的电子计算机X射线断层扫描CT图像，并对CT图像进行二值化处 理，经二值化处理的图像中碳酸盐岩样本的孔隙的像素数与碳酸盐岩样本的颗粒的像素数 不同； 根据孔隙的像素数计算所有孔隙的总面积，根据孔隙的像素数和颗粒的像素数计算所 述碳酸盐岩样本的CT图像面积； 将孔隙的总面积与碳酸盐岩样本的CT图像面积的比值，确定为孔隙度。
2. 如权利要求1所述的方法，其特征在于，还包括： 将经二值化处理的图像中孔洞和孔隙的图形处理为轮廓图，其中，轮廓图是孔洞或孔 隙的图形的轮廓线形成的封闭的图形，经二值化处理的图像中孔洞的像素数与孔隙的像素 数相同； 根据轮廓图的分布方式确定孔洞和孔隙在碳酸盐岩样本上的分布方式。
3. 如权利要求2所述的方法，其特征在于，还包括： 对于每个轮廓图，将轮廓图划分为多个图块，在每个图块中建立与图块边缘相切的楠 圆，多个图块中楠圆的短轴方向相同； 将每个图块短轴的长度与该图块面积相乘，得到每个图块的乘积，将多个图块的乘积 相加，得到的和为与该轮廓图边缘相切的楠圆的短轴长度，确定为孔洞或孔隙的短轴长度。
4. 如权利要求3所述的方法，其特征在于，还包括： 将每个图块长轴的长度与该图块面积相乘，得到每个图块的乘积，将多个图块的乘积 相加，得到的和为与该轮廓图边缘相切的楠圆的长轴长度； 根据轮廓图的长轴长度与短轴长度的大小关系，确定孔洞或孔隙的图形形态。
5. 如权利要求2至4中任一项所述的方法，其特征在于，还包括： 根据孔洞和孔隙的轮廓图，获取孔洞和孔隙的骨架图，骨架图是由线段组成的表示孔 洞或孔隙的拓扑结构的图形； 根据骨架图在碳酸盐岩样本的CT图像上的坐标，确定孔洞或孔隙在碳酸盐岩样本上 的位置。
6. -种碳酸盐岩微观结构的确定装置，其特征在于，包括： 图像处理模块，用于获取碳酸盐岩样本的电子计算机X射线断层扫描CT图像，并对CT 图像进行二值化处理，经二值化处理的图像中碳酸盐岩样本的孔隙的像素数与碳酸盐岩样 本的颗粒的像素数不同； 面积计算模块，用于根据孔隙的像素数计算所有孔隙的总面积，根据孔隙的像素数和 颗粒的像素数计算所述碳酸盐岩样本的CT图像面积； 孔隙度确定模块，用于将孔隙的总面积与碳酸盐岩样本的CT图像面积的比值，确定为 孔隙度。
7. 如权利要求6所述的装置，其特征在于，所述装置还包括： 轮廓图处理模块，用于将经二值化处理的图像中孔洞和孔隙的图形处理为轮廓图，其 中，轮廓图是孔洞或孔隙的图形的轮廓线形成的封闭的图形，经二值化处理的图像中孔洞 的像素数与孔隙的像素数相同； 分布方式确定模块，用于根据轮廓图的分布方式确定孔洞和孔隙在碳酸盐岩样本上的 分布方式。
8. 如权利要求7所述的装置，其特征在于，还包括： 图块划分模块，用于对于每个轮廓图，将轮廓图划分为多个图块，在每个图块中建立与 图块边缘相切的楠圆形，多个图块中楠圆的短轴方向相同； 短轴确定模块，将每个图块短轴的长度与该图块面积相乘，得到每个图块的乘积，将多 个图块的乘积相加，得到的和为与该轮廓图边缘相切的楠圆的短轴长度，确定为孔洞或孔 隙的短轴长度。
9. 如权利要求8所述的装置，其特征在于，还包括： 长轴确定模块，用于将每个图块长轴的长度与该图块面积相乘，得到每个图块的乘积， 将多个图块的乘积相加，得到的和为与该轮廓图边缘相切的楠圆的长轴长度. 图形形态确定模块，用于根据轮廓图的长轴长度与短轴长度的大小关系，确定孔洞或 孔隙的图形形态。
10. 如权利要求7至9中任一项所述的装置，其特征在于，还包括： 骨架图处理模块，用于根据孔洞和孔隙的轮廓图，获取孔洞和孔隙的骨架图，其中，骨 架图是由线段组成的表示孔洞或孔隙的拓扑结构的图形； 位置确定模块，用于根据骨架图在碳酸盐岩样本的CT图像上的坐标，确定孔洞或孔隙 在碳酸盐岩样本上的位置。
【专利摘要】本发明实施例提供了一种碳酸盐岩微观结构的确定方法及装置，其中，该方法包括：获取碳酸盐岩样本的电子计算机X射线断层扫描CT图像，并对CT图像进行二值化处理，经二值化处理的图像中碳酸盐岩样本的孔隙的像素数与碳酸盐岩样本的颗粒的像素数不同；根据孔隙的像素数计算所有孔隙的总面积，根据孔隙的像素数和颗粒的像素数计算所述碳酸盐岩样本的CT图像面积；将孔隙的总面积与碳酸盐岩样本的CT图像面积的比值，确定为孔隙度。该方案基于碳酸盐岩样本的CT图像，可以定量地计算孔隙的面积和碳酸盐岩样本的CT图像面积，进而可以准确、可靠地计算碳酸盐岩样本的孔隙度，有利于对碳酸盐岩储层的微观结构进行定量、可靠的分析。
【IPC分类】G01N15-08
【公开号】CN104535475
【申请号】CN201510009254
【发明人】蔡振忠, 董平川, 邓兴梁, 张正红, 于红枫, 董睿涛, 杨书, 吴子森
【申请人】中国石油大学(北京)
【公开日】2015年4月22日
【申请日】2015年1月8日