CN112233166A - Pore size distribution evaluation method based on porous medium three-dimensional pore space image - Google Patents
Pore size distribution evaluation method based on porous medium three-dimensional pore space image Download PDFInfo
- 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
- radius
- dimensional image
- 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
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 109
- 238000011156 evaluation Methods 0.000 title claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 14
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims abstract description 5
- 229910052753 mercury Inorganic materials 0.000 claims abstract description 5
- 238000012360 testing method Methods 0.000 claims abstract description 5
- 238000011160 research Methods 0.000 claims description 3
- 238000002372 labelling Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000000284 extract Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis 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
The invention discloses a pore size distribution evaluation method based on a porous medium three-dimensional pore space image, which is similar to the testing principle of a mercury intrusion method, fills real pore space by utilizing pore spheres with different predefined radius sizes in a discrete interval until the pores are completely filled, then performs statistical calculation and analysis on the volume of the filled spheres in the pore space to further obtain pore size distribution information, wherein the calculation step comprises the following steps: defining a pore geometry in the discrete three-dimensional image; filling an original three-dimensional image matrix; calculating a maximum pore radius value; indicia of a pore region; and (4) carrying out statistical calculation on the pore fraction. The method for evaluating the pore size distribution has the advantages that the calculated result is high in conformity with the theoretical value, the pore size distribution information of the porous medium is accurately extracted, the pore area can be quickly marked through a logic algorithm, the execution speed is high, and the marking efficiency is high.
Description
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. A pore size distribution evaluation method based on a porous medium three-dimensional pore space image is characterized in that the pore size distribution evaluation method is similar to a mercury intrusion method test principle, pore spheres with different predefined radius sizes in a discrete interval are utilized to fill a real pore space until pores are completely filled, then, the volume of the filled spheres in the pore space is subjected to statistical calculation and analysis, and further pore size distribution information is obtained, and the specific calculation steps are as follows:
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.
2. The method for evaluating pore size distribution based on three-dimensional pore space image of porous medium as claimed in claim 1, wherein said step S1 is performed by assuming sphere of pore geometry in three-dimensional discrete space, and predefining spheres with different radius sizes according to sphere formula.
3. The method for evaluating the pore size distribution based on the three-dimensional pore space image of the porous medium as claimed in claim 1, wherein in step S2, to prevent the three-dimensional image matrix from being read out beyond the boundary in the subsequent calculation, the original three-dimensional image matrix is first expanded and filled in size.
4. The method for evaluating the pore size distribution based on the three-dimensional pore space image of the porous medium 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. The method for evaluating pore size distribution based on three-dimensional pore space image of porous medium as claimed in claim 1, wherein said step S3 is to traverse the pore pixel points in the three-dimensional image, at each pore position, fill the discrete pore spheres with different radius 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.
6. The method as claimed in claim 1, wherein the step S4 is to obtain a sphere with the maximum radius size that can be placed at each pore pixel position in the three-dimensional image, then traverse the pore pixels in the image again, and mark the pixel points that can be covered by the maximum sphere at each position, where the mark value is the corresponding maximum radius size value.
7. The method as claimed in claim 6, wherein the labeled value is based on a principle of priority of a larger labeled value during the labeling process, that is, if the pore radius value to be labeled subsequently is larger, the previously labeled pore radius value is replaced and updated, otherwise, the original labeled value is maintained.
8. The method according to claim 1, wherein the step S5 is to take the marked three-dimensional image as a research object, count the number of pixel points with different mark values of pore radius in the three-dimensional image respectively, obtain the corresponding ratio, and finally obtain the information of pore size distribution of the real pore structure image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010954831.8A CN112233166A (en) | 2020-09-11 | 2020-09-11 | Pore size distribution evaluation method based on porous medium three-dimensional pore space image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010954831.8A CN112233166A (en) | 2020-09-11 | 2020-09-11 | Pore size distribution evaluation method based on porous medium three-dimensional pore space image |
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 | Pore size distribution evaluation method based on porous medium three-dimensional pore space image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112233166A (en) |
Cited By (1)
| 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)
| 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 |
-
2020
- 2020-09-11 CN CN202010954831.8A patent/CN112233166A/en active Pending
Patent Citations (5)
| 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)
| Title |
|---|
| SHUAIBING SONG ET AL: "Improved algorithm for estimating pore size distribution from pore space images of porous media", 《PHYSICAL REVIEW》 * |
Cited By (1)
| 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 |
|---|---|---|
| CN104809456B (en) | A kind of Three-dimensional target recognition method based on two valued description symbol | |
| CN109242985B (en) | Method for determining key parameters of pore structure from three-dimensional image | |
| CN112233166A (en) | Pore size distribution evaluation method based on porous medium three-dimensional pore space image | |
| CN110688958A (en) | Intersection identification method based on GoogleLeNet neural network | |
| CN110084355B (en) | Grid scale optimization method of large-amount interaction particle motion simulation system | |
| CN111475597A (en) | Non-rigid grid coding, spatial object unique identification and query method and device | |
| CN110765963A (en) | Vehicle brake detection method, device, equipment and computer readable storage medium | |
| CN112308040A (en) | River sewage outlet detection method and system based on high-definition images | |
| CN114998108A (en) | Vector data optimization processing method and system | |
| CN112435288B (en) | Image-based pore feature calculation method | |
| CN110806371A (en) | Method for constructing single-weight coal rock pore structure model | |
| CN110008837B (en) | Quick realization method for calculating POI according to closed polygon of road network | |
| CN113673186A (en) | Cartesian grid rapid generation method based on STL file | |
| CN113987082A (en) | Three-dimensional monomer model data model expression method based on geographic entity coding | |
| CN111446968B (en) | Vector space data multistage compression method | |
| CN112164087A (en) | Super-pixel segmentation method and device based on edge constraint and segmentation boundary search | |
| CN109118565B (en) | Electric power corridor three-dimensional model texture mapping method considering shielding of pole tower power line | |
| CN114322865B (en) | Method and device for measuring semiconductor device and storage medium | |
| CN115717887A (en) | Star point fast extraction method based on gray distribution histogram | |
| CN116341875A (en) | Engineering detection system and method applied to building construction site | |
| CN115995024A (en) | Image classification method based on class diagram neural network | |
| CN115984687A (en) | River work moving bed model test water boundary measuring method, device, equipment and medium | |
| CN112017273B (en) | Multi-resolution core pore fusion method based on mode density function simulation algorithm | |
| CN106658009A (en) | Improved double-ring sampling method and device for panoramic video | |
| CN112199984A (en) | Target rapid detection method of large-scale remote sensing image |
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 |
Application publication date: 20210115 |
|
| RJ01 | Rejection of invention patent application after publication |