CN112233166A - A method for evaluating pore size distribution based on three-dimensional pore space images of porous media - Google Patents

A method for evaluating pore size distribution based on three-dimensional pore space images of porous media Download PDF

Info

Publication number
CN112233166A
CN112233166A CN202010954831.8A CN202010954831A CN112233166A CN 112233166 A CN112233166 A CN 112233166A CN 202010954831 A CN202010954831 A CN 202010954831A CN 112233166 A CN112233166 A CN 112233166A
Authority
CN
China
Prior art keywords
pore
size distribution
dimensional
image
pore size
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010954831.8A
Other languages
Chinese (zh)
Inventor
宋帅兵
王磊
刘和武
徐楠
涂庆毅
焦振华
王飞
刘瑜
张通
张贵生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Anhui University of Science and Technology
Original Assignee
Anhui University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Anhui University of Science and Technology filed Critical Anhui University of Science and Technology
Priority to CN202010954831.8A priority Critical patent/CN112233166A/en
Publication of CN112233166A publication Critical patent/CN112233166A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Generation (AREA)

Abstract

本发明公开一种基于多孔介质三维孔隙空间图像的孔径分布评估方法,所述孔径分布评估方法类比压汞法的测试原理,利用离散区间内部预定义的不同半径尺寸的孔隙球体,对真实孔隙空间进行填充,直至将孔隙完全填充,然后,对孔隙空间内部已填充的球体体积进行统计计算分析,进而获得孔径分布信息,计算步骤包括:定义离散三维图像中孔隙几何形状;原始三维图像矩阵的填充;最大孔隙半径值的计算;孔隙区域的标记;孔隙占比分数统计计算。本发明孔径分布评估方法计算结果与理论值的相符度高,精准的对多孔介质的孔径分布信息进行提取,通过逻辑算法,能够快速标记出孔隙区域,执行速度快,标记效率高。

Figure 202010954831

The invention discloses a pore size distribution evaluation method based on a three-dimensional pore space image of a porous medium. The pore size distribution evaluation method is analogous to the testing principle of mercury porosimetry. Filling is performed until the pores are completely filled, and then, statistical calculation and analysis are performed on the volume of the sphere filled in the pore space, and then the pore size distribution information is obtained. The calculation steps include: defining the geometric shape of the pores in the discrete 3D image; filling the original 3D image matrix ; Calculation of the maximum pore radius value; Marking of the pore area; Statistical calculation of the proportion of pores. The calculation result of the pore size distribution evaluation method of the present invention has a high degree of agreement with the theoretical value, accurately extracts the pore size distribution information of the porous medium, and can quickly mark the pore area through a logical algorithm, with fast execution speed and high marking efficiency.

Figure 202010954831

Description

Pore size distribution evaluation method based on porous medium three-dimensional pore space image
Technical Field
The invention relates to a pore size distribution evaluation method, in particular to a pore size distribution evaluation method based on a three-dimensional pore space image of a porous medium.
Background
Thanks to the rapidly developing digital microscopic imaging techniques (CT and FIB/SEM), the microscopic three-dimensional pore structure inside the porous media material is visualized. Based on the three-dimensional pore structure image, corresponding pore size distribution information can be directly extracted. The existing algorithm is mainly similar to the mercury intrusion testing principle, and calculates the pore diameter distribution by executing different types of distance transformation (Euclidean distance or chessboard distance transformation), however, due to the discrete nature of the digital image and the lack of clear and definite definition of pore geometry, the accuracy of the calculation result is poor in some cases, and besides, the execution efficiency of the existing algorithm is low for three-dimensional large-scale images.
Disclosure of Invention
The invention aims to provide a pore size distribution evaluation method based on a three-dimensional pore space image of a porous medium
The purpose of the invention can be realized by the following technical scheme:
a pore size distribution evaluation method based on a porous medium three-dimensional pore space image is similar to a mercury intrusion method test principle, and comprises the following specific calculation steps of filling real pore spaces by utilizing pore spheres with different predefined radius sizes in a discrete interval until the pores are completely filled, and then carrying out statistical calculation and analysis on the filled sphere volumes in the pore spaces to further obtain pore size distribution information:
s1, defining the geometrical shape of the pore in the discrete three-dimensional image;
s2, filling an original three-dimensional image matrix;
s3, calculating a maximum pore radius value;
s4, marking of pore regions;
and S5, statistically calculating the pore ratio fraction.
Further, the step S1 predefines spheres of different radius sizes according to a sphere formula by assuming a sphere of pore geometry in a three-dimensional discrete space.
Further, in order to prevent the three-dimensional image matrix from being read out beyond the boundary in the subsequent calculation in step S2, the size of the original three-dimensional image matrix is first expanded and filled.
Further, the expansion size is equal to the size of the three-dimensional image matrix in the same direction, and all filling values of the filling operation are selected as pore pixels.
Further, the step S3 is to traverse the pore pixel points in the three-dimensional image, and at each pore position, fill the discrete pore spheres with different radii that have been predefined previously, to find the sphere with the largest radius size that can be placed at the position, and record and store the radius value.
Further, in step S4, after obtaining the sphere with the maximum radius size that can be placed at each pore pixel position in the three-dimensional image, traversing the pore pixels in the image again, and marking the pixel points that can be covered by the maximum sphere at each position, where the marking value is the corresponding maximum radius size value.
Furthermore, the marking value follows a principle of priority of a larger marking value in the marking process, namely if the aperture radius value to be marked subsequently is larger, the previously marked aperture radius value is replaced and updated, otherwise, the original marking value is maintained unchanged.
Further, in step S5, the marked three-dimensional image is used as a research object, the number of pixel points having different pore radius mark values inside the three-dimensional image is respectively counted to obtain the corresponding ratio, and finally the pore size distribution information of the real pore structure image is obtained.
The invention has the beneficial effects that:
1. the method for evaluating the pore size distribution has high conformity of the calculation result and the theoretical value, and accurately extracts the pore size distribution information of the porous medium;
2. the aperture distribution evaluation method can rapidly mark the aperture area through a logic algorithm, and has high execution speed and high marking efficiency.
Drawings
The invention will be further described with reference to the accompanying drawings.
FIG. 1 is a schematic diagram of the discrete sphere construction of the present invention;
FIG. 2 is a schematic of an aperture pixel of the present invention;
FIG. 3 is a schematic diagram of the logic algorithm for void region marking of the present invention;
FIG. 4 is a schematic illustration of the pore size distribution of the present invention;
FIG. 5 is a schematic diagram of the aperture frequency of a pixel of the present invention;
FIG. 6 is a schematic diagram of the aperture frequency of a pixel of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
A pore size distribution evaluation method based on a porous medium three-dimensional pore space image is based on a testing principle of a mercury intrusion method, and comprises the following specific calculation steps of filling real pore space by utilizing pore spheres with different predefined radius sizes in a discrete interval until pores are completely filled, and then carrying out statistical calculation analysis on the volume of the filled spheres in the pore space to further obtain pore size distribution information:
s1, defining the pore geometry in the discrete three-dimensional image: assuming that the geometric shape of the pore space in the three-dimensional discrete space is a sphere, predefining the spheres with different radius sizes according to a sphere formula, wherein the calculation formula is as follows:
(x-a)2+(y-b)2+(z-c)2≤R2
wherein the X, Y, Z variables are each integers;
a series of discrete spheres with different radii are established according to the formula, as shown in fig. 1, the square values of the radius R in fig. 1 are 0, 1, 2:
s2, filling of an original three-dimensional image matrix: in order to prevent the three-dimensional image matrix from being read out from the boundary, the size of the original three-dimensional image matrix is expanded and filled, wherein the expanded size is equal to the size of the three-dimensional image matrix in the same direction, and the filling values are all selected as pore pixels, as shown in fig. 2, wherein 0 represents a pore, and 1 represents a matrix:
s3, calculation of maximum pore radius value: traversing pore pixel points in the three-dimensional image, filling discrete pore spheres with different predefined radiuses at each pore position successively to find a sphere with the maximum radius size which can be placed at the position, and recording and storing the radius value;
s4, labeling of pore region: after obtaining a sphere with the maximum radius size which can be placed at each pore pixel position in the three-dimensional image, traversing the pore pixels in the image again, and marking pixel points which can be covered by the maximum sphere at each position, wherein the marking value is the corresponding maximum radius size value;
in the marking process, a principle of priority of a larger marking value is followed, namely if the subsequent pore radius value to be marked is larger, the previously marked pore radius value is replaced and updated, otherwise, the original marking value is maintained unchanged, and the marking process is as shown in fig. 3;
s5, statistical calculation of pore ratio fraction: taking the marked three-dimensional image as a research object, respectively counting the number of pixel points with different pore radius mark values in the three-dimensional image to obtain corresponding proportion, and finally obtaining the aperture distribution information of the real pore structure image.
Example 1
Selecting an artificially synthesized porous medium material with known pore size distribution information, and verifying the accuracy of the algorithm, wherein the method specifically comprises the following steps:
p1, artificially synthesized porous medium material and pore size distribution thereof, as shown in fig. 4 and 5;
p2, comparing the calculation results of the algorithm and the existing algorithm (the Munch algorithm and the Yang algorithm) with theoretical values respectively, as shown in FIG. 6, and verifying the results to show that the calculation result of the aperture distribution of the algorithm is better consistent with the theoretical values, while the existing algorithm and the theoretical values have larger deviation in the two aspects of the aperture interval range and the corresponding volume ratio, and the algorithm can accurately extract the aperture distribution information of the porous medium.
In the description herein, references to the description of "one embodiment," "an example," "a specific example" or the like are intended to mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed.

Claims (8)

1.一种基于多孔介质三维孔隙空间图像的孔径分布评估方法,其特征在于,所述孔径分布评估方法类比压汞法的测试原理,利用离散区间内部预定义的不同半径尺寸的孔隙球体,对真实孔隙空间进行填充,直至将孔隙完全填充,然后,对孔隙空间内部已填充的球体体积进行统计计算分析,进而获得孔径分布信息,具体计算步骤如下:1. a pore size distribution evaluation method based on a three-dimensional pore space image of a porous medium, it is characterized in that, described pore size distribution evaluation method is analogous to the test principle of mercury porosimetry, utilizes the pore spheres of different radius sizes predefined inside the discrete interval, to The real pore space is filled until the pores are completely filled. Then, the volume of the filled spheres in the pore space is statistically calculated and analyzed to obtain the pore size distribution information. The specific calculation steps are as follows: S1、定义离散三维图像中孔隙几何形状;S1. Define the pore geometry in the discrete three-dimensional image; S2、原始三维图像矩阵的填充;S2, the filling of the original three-dimensional image matrix; S3、最大孔隙半径值的计算;S3. Calculation of the maximum pore radius value; S4、孔隙区域的标记;S4, the marking of the pore area; S5、孔隙占比分数统计计算。S5. Statistical calculation of the proportion of pores. 2.根据权利要求1所述的一种基于多孔介质三维孔隙空间图像的孔径分布评估方法,其特征在于,所述步骤S1通过假设三维离散空间中孔隙几何形状球体,根据球体公式对不同半径尺寸的球体进行预定义。2. A method for evaluating pore size distribution based on a three-dimensional pore space image of a porous medium according to claim 1, characterized in that, in step S1, by assuming a sphere of pore geometry in a three-dimensional discrete space, according to the sphere formula, different radius sizes are evaluated. sphere is predefined. 3.根据权利要求1所述的一种基于多孔介质三维孔隙空间图像的孔径分布评估方法,其特征在于,所述步骤S2中为防止后续计算中三维图像矩阵的越界读取,首先对原始三维图像矩阵的尺寸进行扩张和填充操作。3. A method for evaluating the pore size distribution based on a three-dimensional pore space image of a porous medium according to claim 1, wherein in the step S2, in order to prevent the out-of-bounds reading of the three-dimensional image matrix in the subsequent calculation, the original three-dimensional The dimensions of the image matrix are expanded and filled. 4.根据权利要求3所述的一种基于多孔介质三维孔隙空间图像的孔径分布评估方法,其特征在于,所述扩张尺寸与同方向三维图像矩阵的尺寸相等,填充操作的填充值全部选为孔隙像素。4. A method for evaluating pore size distribution based on three-dimensional pore space images of porous media according to claim 3, wherein the expansion size is equal to the size of the three-dimensional image matrix in the same direction, and the filling values of the filling operation are all selected as Pore pixels. 5.根据权利要求1所述的一种基于多孔介质三维孔隙空间图像的孔径分布评估方法,其特征在于,所述步骤S3是对三维图像中的孔隙像素点进行遍历,在每一个孔隙位置处,将先前已经预定义好的不同半径离散孔隙球体逐次对其进行填充,以寻找该位置处所能放置的具有最大半径尺寸的球体,并对该半径值进行记录和存储。5. A method for evaluating pore size distribution based on a three-dimensional pore space image of a porous medium according to claim 1, wherein the step S3 is to traverse the pore pixel points in the three-dimensional image, and at each pore position , fill the previously predefined discrete pore spheres with different radii successively to find the sphere with the largest radius size that can be placed at the position, and record and store the radius value. 6.根据权利要求1所述的一种基于多孔介质三维孔隙空间图像的孔径分布评估方法,其特征在于,所述步骤S4是获得三维图像中每一个孔隙像素位置处所能放置的最大半径尺寸的球体后,再次对图像中的孔隙像素进行遍历,并对每一位置处最大球体所能覆盖的像素点进行标记,标记值为其相应的最大半径尺寸值。6. A method for evaluating pore size distribution based on a three-dimensional pore space image of a porous medium according to claim 1, wherein the step S4 is to obtain the maximum radius size that can be placed at each pore pixel position in the three-dimensional image. After the sphere is formed, the pore pixels in the image are traversed again, and the pixels that can be covered by the largest sphere at each position are marked, and the marked value is its corresponding maximum radius size value. 7.根据权利要求6所述的一种基于多孔介质三维孔隙空间图像的孔径分布评估方法,其特征在于,所述标记值在标记过程中,遵循较大标记值优先原则,即若后续待标记孔隙半径值较大,则对先前已标记的孔隙半径值进行替换更新,否则维持原有标记值不变。7 . The method for evaluating pore size distribution based on a three-dimensional pore space image of a porous medium according to claim 6 , wherein, in the marking process, the marking value follows the principle of prioritizing the larger marking value, that is, if the marking value is to be marked subsequently. 8 . If the pore radius value is larger, the previously marked pore radius value will be replaced and updated, otherwise the original marked value will remain unchanged. 8.根据权利要求1所述的一种基于多孔介质三维孔隙空间图像的孔径分布评估方法,其特征在于,所述步骤S5是把已完成标记的三维图像作为研究对象,分别统计其内部具有不同孔隙半径标记值的像素点的数目,获得其相应的占比,最终得到真实孔隙结构图像的孔径分布信息。8. A method for evaluating the pore size distribution based on a three-dimensional pore space image of a porous medium according to claim 1, wherein the step S5 is to take the marked three-dimensional image as a research object, and count the internal The number of pixel points of the pore radius marked value is obtained, and its corresponding proportion is obtained, and finally the pore size distribution information of the real pore structure image is obtained.
CN202010954831.8A 2020-09-11 2020-09-11 A method for evaluating pore size distribution based on three-dimensional pore space images of porous media Pending CN112233166A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010954831.8A CN112233166A (en) 2020-09-11 2020-09-11 A method for evaluating pore size distribution based on three-dimensional pore space images of porous media

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010954831.8A CN112233166A (en) 2020-09-11 2020-09-11 A method for evaluating pore size distribution based on three-dimensional pore space images of porous media

Publications (1)

Publication Number Publication Date
CN112233166A true CN112233166A (en) 2021-01-15

Family

ID=74117257

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010954831.8A Pending CN112233166A (en) 2020-09-11 2020-09-11 A method for evaluating pore size distribution based on three-dimensional pore space images of porous media

Country Status (1)

Country Link
CN (1) CN112233166A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112967378A (en) * 2021-03-03 2021-06-15 安徽理工大学 Rapid random reconstruction method for large-size pore crack structure in porous medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105806765A (en) * 2016-04-13 2016-07-27 南京大学(苏州)高新技术研究院 Refined characterization method of micro CT scanning soil body space pore structure
CN107991710A (en) * 2017-10-23 2018-05-04 中国石油天然气股份有限公司 Reservoir pore size distribution obtaining method and device
CN108491677A (en) * 2018-07-04 2018-09-04 河海大学 Pore character statistical method based on the micro pore model for improving maximum ball
CN110018108A (en) * 2019-05-14 2019-07-16 中国科学院地质与地球物理研究所 The measuring method and equipment in blowhole aperture
CN110853138A (en) * 2019-11-21 2020-02-28 科吉思石油技术咨询(北京)有限公司 Dual-medium carbonate rock pore-fracture dual-network model construction method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105806765A (en) * 2016-04-13 2016-07-27 南京大学(苏州)高新技术研究院 Refined characterization method of micro CT scanning soil body space pore structure
CN107991710A (en) * 2017-10-23 2018-05-04 中国石油天然气股份有限公司 Reservoir pore size distribution obtaining method and device
CN108491677A (en) * 2018-07-04 2018-09-04 河海大学 Pore character statistical method based on the micro pore model for improving maximum ball
CN110018108A (en) * 2019-05-14 2019-07-16 中国科学院地质与地球物理研究所 The measuring method and equipment in blowhole aperture
CN110853138A (en) * 2019-11-21 2020-02-28 科吉思石油技术咨询(北京)有限公司 Dual-medium carbonate rock pore-fracture dual-network model construction method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
SHUAIBING SONG ET AL: "Improved algorithm for estimating pore size distribution from pore space images of porous media", 《PHYSICAL REVIEW》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112967378A (en) * 2021-03-03 2021-06-15 安徽理工大学 Rapid random reconstruction method for large-size pore crack structure in porous medium

Similar Documents

Publication Publication Date Title
CN109978901B (en) Quick and accurate circular detection and circle center positioning method
CN105844602B (en) A kind of airborne LIDAR point cloud three-dimensional filtering method based on volume elements
CN109345625B (en) Rock core image self-adaptive partition three-dimensional reconstruction method
CN114119917A (en) Visual hidden danger ranging method based on high-precision three-dimensional model of power transmission line
CN112308040A (en) River sewage outlet detection method and system based on high-definition images
CN112233166A (en) A method for evaluating pore size distribution based on three-dimensional pore space images of porous media
CN116452583A (en) Point cloud defect detection method, device, system and storage medium
CN116012545B (en) Multi-scale digital core modeling method, system, storage medium and application
CN112001453B (en) Accuracy calculation method and device for video event detection algorithm
CN114742838B (en) Voxelization plane segmentation method based on clustering algorithm
CN102722910B (en) Volume mesh scene continuous collision detection method based on separation axis removal
CN115861566A (en) A LiDAR Point Cloud Building Integrity Enhancement Method Based on Intensity Voxel Model
CN119180786B (en) Prefabricated box girder template detection method and system
CN112017273B (en) Multi-resolution Core Pore Fusion Method Based on Model Density Function Simulation Algorithm
CN113836484B (en) Self-adaptive point cloud thinning method based on path point adjacent domain and ground filtering
TW200408034A (en) Pattern measurement method, manufacturing method of semiconductor device using the same, and pattern measurement device
CN112435288B (en) Image-based pore feature calculation method
CN111223787B (en) Groove structure measuring method, device, equipment and medium of three-dimensional memory
CN112284287B (en) A Stereo Vision 3D Displacement Measurement Method Based on Grayscale Feature of Structure Surface
JP2009266198A (en) Method for calculating thickness direction and thickness of solid element, and apparatus for calculating thickness direction and thickness of solid element
CN113204901B (en) Unit cell region division method of variable cross-section three-dimensional five-direction braided composite member
CN114529675B (en) A 3D model optimization method for fast local Delaunay transformation
CN114627466B (en) A method and system for automatically positioning and detecting floating objects in a tilted real-life model
CN115937295B (en) Bubble shape type discriminating method based on aggregate ideas
CN116821398B (en) Data set acquisition method for road defect recognition model training

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210115