CIESC Journal ›› 2018, Vol. 69 ›› Issue (S1): 26-33.DOI: 10.11949/j.issn.0438-1157.20171382

Previous Articles     Next Articles

Medial axis extraction method of pore structure based on fully parallel thinning algorithm

HU Peiyu1, WANG Shugang1, SONG Shuanglin1,3, JIANG Shuang2, LIANG Yuntao3   

  1. 1 Faculty of Infrastructure Engineering, Dalian University of Technology, Dalian 116024, Liaoning, Liaoning, China;
    2 College of Civil Engineering, Dalian Minzu University, Dalian 116600, Liaoning, China;
    3 State Key Laboratory of Coal Safety, CCTEG Shenyang Research Institute, Shenyang 110016, Liaoning, China
  • Received:2017-10-17 Revised:2018-02-27 Online:2018-09-30 Published:2018-09-30
  • Supported by:

    supported by the National Natural Science Foundation of China (51774181) and the Opening Research Fund of State Key Laboratory of Coal Mine Safety Technology (SKLCMST104).

基于完全并行细化算法的孔隙网络中轴提取方法

胡沛裕1, 王树刚1, 宋双林1,3, 蒋爽2, 梁运涛3   

  1. 1 大连理工大学建设工程学部, 辽宁 大连 116024;
    2 大连民族大学土木工程学院, 辽宁 大连 116600;
    3 煤科集团沈阳研究院有限公司煤矿安全技术国家重点实验室, 辽宁 沈阳 110016
  • 通讯作者: 蒋爽,E-mail:shjiang@foxmail.com
  • 基金资助:

    国家自然科学基金项目(51774181);煤矿安全技术国家重点实验室开放课题基金资助项目(SKLCMST104)。

Abstract:

Building a precise pore network by using 3D reconstruction lays the foundation for studying the fluid flow as well as heat and mass transfer in porous media, in which the medial axis extraction of pore structure is the key to reconstruction. Among the thinning algorithms to extract medial axis, simple points method has its advantages over accuracy and parallel operation, whereas its application to pore structure analysis is relatively few. Based on the fully parallel thinning algorithm of the simple points method, a medial axis extraction method specializing in porous media is developed, consisting of the following four major steps:image processing, reconstruction matrix generating, layer-by-layer marking and parallel thinning. Applying the proposed method to 230 objects with typically ordered pore structures and a reconstruction model of an activated carbon accumulation, the medial axes topologically equivalent to each object are obtained and the pore radiuses are calculated fast and precisely. By comparing the pore volume calculated from the media axis and its corresponding radius with actual pore volume, it is found that the relative error of the pore volume is lower than 6% eventually. The results show that the proposed method is capable of obtaining the medial axis preserving the geometry and connectivity commendably and resisting the influence of surface irregular variation effectively, hence can be used for nondestructive quantitative analysis and pore network extraction.

Key words: porous media, model, pore network, medial axis, algorithm

摘要:

应用三维重构技术建立介尺度下的多孔介质孔隙网络模型是分析多孔介质内流体流动与传热传质过程的基础,其中孔隙网络中轴的提取是三维重构的关键环节。在提取中轴的细化算法中,简单点集合细化方法具有准确、可并行的特点,但较少应用于多孔介质的结构分析。基于简单点集合细化方法中的完全并行细化算法,建立了一种面向多孔介质的孔隙网络中轴提取方法,包含图像处理、重构矩阵生成、逐层标记、并行细化四个主要步骤。以计算机生成的230个典型的有序多孔介质模型及1个活性炭堆积试样重构模型为研究对象,应用本方法获得了与研究对象的孔隙空间具有拓扑等价性的孔隙网络中轴,快速准确地计算了孔隙半径,利用中轴形态与孔隙半径计算孔隙占据的空间大小与实际体积比较,发现研究对象的孔隙体积计算值相对误差最终稳定在6%以下。结果表明,由本方法获得的中轴反映了孔隙空间的几何形状、减小了孔隙与固体界面处不规则变化的干扰,可用于多孔介质孔隙结构的无损定量分析及孔隙网络建模。

关键词: 多孔介质, 模型, 孔隙网络, 中轴, 算法

CLC Number: