CN116342541B - Rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction - Google Patents
Rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction Download PDFInfo
- Publication number
- CN116342541B CN116342541B CN202310317640.4A CN202310317640A CN116342541B CN 116342541 B CN116342541 B CN 116342541B CN 202310317640 A CN202310317640 A CN 202310317640A CN 116342541 B CN116342541 B CN 116342541B
- Authority
- CN
- China
- Prior art keywords
- pore
- permeability
- image
- rock
- adjacent
- 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.)
- Active
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 103
- 230000035699 permeability Effects 0.000 title claims abstract description 54
- 239000002689 soil Substances 0.000 title claims abstract description 38
- 230000004927 fusion Effects 0.000 title claims abstract description 22
- 238000004364 calculation method Methods 0.000 title claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 27
- 238000004891 communication Methods 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 26
- 238000011156 evaluation Methods 0.000 claims description 9
- 239000002131 composite material Substances 0.000 claims description 6
- 238000000638 solvent extraction Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000010884 ion-beam technique Methods 0.000 description 2
- 230000003204 osmotic effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 102100032219 Cathepsin D Human genes 0.000 description 1
- 101000869010 Homo sapiens Cathepsin D Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004626 scanning electron microscopy Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20036—Morphological image processing
Abstract
The invention discloses a rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction, which comprises the following steps: acquiring a three-dimensional microscopic image of a rock-soil body sample to be detected, and processing the three-dimensional microscopic image to acquire a pore-matrix binary image of the rock-soil body sample to be detected; based on the pore-matrix binary images, acquiring pore communication components of adjacent pore-matrix binary images, and carrying out fusion reconstruction processing on pores in the adjacent pore-matrix binary images to acquire a synthesized binary image; acquiring the permeability of the synthesized binary image based on the synthesized binary image; and acquiring the permeability of the rock-soil body sample to be detected based on the permeability of the synthesized binary image. According to the invention, by improving the fusion reconstruction link of the adjacent binary images, the obtained pores are more reasonable, and the phenomenon that the permeability is smaller due to smaller pores is avoided.
Description
Technical Field
The invention belongs to the technical field of rock-soil body permeability measurement and calculation, and particularly relates to a rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction.
Background
The permeability is a key technical index of the rock-soil material, and has important significance for scientific research and engineering application of various related industries. The scientific and convenient acquisition of the permeability is one of important preconditions for subsequent research and analysis of the physical properties of the rock-soil body. Conventional methods use experimental equipment to test for permeability values. Because of the large difference between different kinds of materials, the applicable test methods and equipment are different. At the same time, the test may cause irreversible damage or fouling of the sample, which can be inconvenient for subsequent further research analysis of the sample.
The research finds that the internal pores of the rock-soil material are key determinants of the seepage capability of the rock-soil body. Many scholars extract and calculate the size and distribution of pore structures to predict permeability. Microscopic three-dimensional imaging techniques are a common means for evaluating three-dimensional pore structures described in recent years, which observe pore structures using microscopic three-dimensional imaging equipment, thereby evaluating permeability. Such as computed tomography, focused ion beam microscopy, and wide ion beam scanning electron microscopy. However, these methods take a long time to calculate, and some of them can take several hours, which is disadvantageous for the task of predicting a large number of samples and batches. In the prior invention patent, a rock-soil mass rapid prediction method based on three-dimensional images exists. According to the method, according to the relative change of the pores of each adjacent image, each adjacent binary image is synthesized, and a binary image which can reflect the change of the pore channel of the sample is obtained. However, the method only considers the intersection of the adjacent image pores, ignores the rest parts of the pores, and ensures that the volume of the pores participating in calculation is smaller and the predicted value of the permeability is lower.
Disclosure of Invention
Aiming at the defects that the existing three-dimensional image-based permeability calculation method is unreasonable in fusion and reconstruction of pore channels, long in calculation time, inapplicable to a plurality of batches of mass measurement tasks and the like, the invention provides a rock-soil mass permeability calculation method based on adjacent image pore fusion and reconstruction, which optimizes and improves pore fusion and reconstruction links of adjacent images, so that the calculated pores can better reflect the morphological positions and the changes of the pores in the region, and unreasonable pore channels are avoided.
In order to achieve the above purpose, the invention provides a rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction, which comprises the following steps:
acquiring a three-dimensional microscopic image of a rock-soil body sample to be detected, and processing the three-dimensional microscopic image to acquire a pore-matrix binary image of the rock-soil body sample to be detected;
based on the pore-matrix binary images, acquiring pore communication components of adjacent pore-matrix binary images, and carrying out fusion reconstruction processing on pores in the adjacent pore-matrix binary images to acquire a synthesized binary image;
acquiring the permeability of the synthesized binary image based on the synthesized binary image;
and acquiring the permeability of the rock-soil body sample to be detected based on the permeability of the synthesized binary image.
Optionally, the method for acquiring the synthesized binary image includes:
based on the adjacent pore-matrix binary images, acquiring pore communication components of the adjacent pore-matrix binary images according to neighborhood 8 pixel communication;
and dividing the pore pixels in the adjacent pore-matrix binary images based on the pore communication components, and acquiring the synthesized binary image.
Optionally, the method for dividing the pore pixels in the adjacent pore-matrix binary image based on the pore communication component comprises:
acquiring intersection, intersection edge pixels, union and union edge pixels of the pore connected components;
acquiring the shortest distance between each pixel in the union and the intersection edge pixel, and recording the shortest distance as a first distance;
acquiring the shortest distance between each pixel in the union and the edge pixel of the union, and recording the shortest distance as a second distance;
acquiring an evaluation index according to the first distance and the second distance;
based on the evaluation index, the pore pixels in the adjacent pore-matrix binary images are divided.
Optionally, the specific method for obtaining the first distance, the second distance and the evaluation index includes:
tp u ＝d u2 /(d u1 +d u2 )
wherein d u1 At a first distance d u2 For a second distance tp u For evaluation index, x u And y u For the u-th pixel image array coordinates, x u1 And y u1 Is the closest intersection edge pixel coordinate with the u-th pixel, x u2 And y u2 Is the closest union edge pixel coordinate to the u-th pixel.
Optionally, the method for acquiring the permeability of the synthesized binary image based on the synthesized binary image includes:
according to the synthesized binary image, pore size distribution of pores in the synthesized binary image is obtained through a continuous pore size algorithm;
and obtaining the permeability of the synthesized binary image according to the Hagen-Poiseuille law and the pore diameter distribution of the pores.
Optionally, the permeability of the synthesized binary image specifically includes:
wherein k is i The permeability of the corresponding region of the binary image is synthesized for the ith sheet, i=1, 2, … … and N-1; a is the cross-sectional area of the sample, R ik For the radius of each micro-pore of the region, S ik Is of radius R ik The pores occupy area.
Optionally, the obtaining the permeability of the rock-soil body sample to be measured specifically includes:
wherein N is the number of pore-matrix binary images, k i The permeability of the corresponding region of the binary image is synthesized for the ith sheet, i=1, 2, … … and N-1; a is the cross-sectional area of the sample; r is R ik For the radius of each micro-pore of the region, S ik Is of radius R ik The pores occupy area.
The invention has the technical effects that: (1) Compared with the prior art, the method has the advantages that the fusion reconstruction link of the adjacent binary images is improved, the form and the position of the pores in the original two adjacent binary images are considered, the obtained pores are more reasonable, and the permeability is prevented from being smaller due to smaller pores; (2) Compared with other permeability calculation methods based on three-dimensional microscopic images, the method has the advantages that the calculation time is remarkably reduced, and the method is better applied to large-size and multi-batch calculation tasks.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application, illustrate and explain the application and are not to be construed as limiting the application. In the drawings:
FIG. 1 is a schematic flow chart of a rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction in an embodiment of the invention;
FIG. 2 is a schematic diagram of the main links of a binary image synthesis in an embodiment of the present invention;
fig. 3 is a diagram showing an example of the result of pore size distribution measurement according to an embodiment of the present invention.
Detailed Description
It should be noted that, in the case of no conflict, the embodiments and features in the embodiments may be combined with each other. The present application will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.
It should be noted that the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer executable instructions, and that although a logical order is illustrated in the flowcharts, in some cases the steps illustrated or described may be performed in an order other than that illustrated herein.
As shown in fig. 1, in this embodiment, a method for calculating permeability of a rock-soil body based on fusion reconstruction of adjacent image pores is provided, including:
s1, acquiring a three-dimensional microscopic image of a rock-soil body sample to be calculated, and carrying out image noise reduction and image binarization to acquire a pore-matrix binary image.
S2, obtaining a binary image showing pore network change in order to realize scientific and rapid permeability, namely providing a synthetic binary image showing the volume and the position of a part of sample pore seepage channel according to the relative change of each image pore. Specifically, a composite binary image is obtained by pore fusion reconstruction of every two adjacent two-dimensional binary images. Namely, if the three-dimensional image contains N binary images, N-1 combined binary images are obtained. In addition, aiming at the defects of the prior invention patent, the link is improved and innovated, so that the pore volume in the synthesized binary image is more reasonable.
The method comprises the following steps:
s2.1 obtains the ith and the (i+1) th two-dimensional binary maps (i=1, 2, … …, N-1), and obtains binary communication components according to the neighborhood 8 pixel communication condition.
S2.2, fusing and reconstructing pores of adjacent binary images. Intersection and intersection edge pixels and union edge pixels of each aperture-connected component of the ith and (i+1) th two-dimensional binary maps, and their respective edge pixels, are first determined. Assume that the intersection is U pixels total. Then, the distance d between the u-th pixel and the nearest intersection edge pixel is obtained u1 Distance d from nearest union edge pixel u2 Dividing the evaluation index tp u (u=1, 2, … …, U), distance d of the U-th pixel from the nearest intersection edge pixel u1 I.e., first distance, distance d of the u-th pixel from nearest union edge pixel u2 I.e. the second distance.
tp u ＝d u2 /(d u1 +d u2 ) (3)
Wherein x is u And y u For the u-th pixel image array coordinates, x u1 And y u1 Is the closest intersection edge pixel coordinate with the u-th pixel, x u2 And y u2 Is the closest union edge pixel coordinate to the u-th pixel. Next, it is determined whether the pixel at the corresponding position of the ith pixel in the ith synthesized binary image belongs to a pore: if the u-th pixel belongs to the intersection or tp u More than or equal to 0.5, the pixel at the corresponding position in the ith synthesized binary image is a pore pixel; otherwise, the corresponding position is a matrix pixel. Finally, traversing u, and determining the ith synthesized binary image element classification by the method: which of the ith synthesized binary patterns is the aperture pixel and which is the matrix pixel.
As shown in fig. 2, taking two holes at the same position of two adjacent binary images as an example, the communication components corresponding to the two holes and the intersection and union of the communication components are obtained. And then according to the shortest distance d with the intersection edge u1 And the shortest distance d of the union edges u2 And determining whether to divide into aperture pixels in the composite binary image, and finally obtaining the composite binary image according to the principle.
S2.3 starting from i=2, S2.1 and S2.2 are repeated to obtain the remaining total N-2 synthetic binary maps.
S3, calculating permeability of synthesized binary image
In the subsequent links, the ith synthesized binary image represents the pore morphology between the ith Zhang Erzhi image and the (i+1) binary image, and the pore morphology of the ith layer of the sample and the change thereof are reflected.
S3.1, acquiring the pore size distribution of the ith synthesized binary image by using a continuous pore size algorithm (continue pore size distribution, CPSD). Assuming that the pore communication components are a continuum, the pore size variation is continuous. The area or volume filled is calculated by scanning the pore portion using a circle of specified radius to obtain the ratio of the size pores in the overall pore space, thereby obtaining the area and number of each size pore in principle. Specifically, starting from i=1. Taking the 1 st pore binary diagram as an example, the measurement result of the pore size distribution result is shown in fig. 3.
S3.2, carrying out gas permeability calculation by combining with the Hagen-Poiseuille law. The method comprises the following steps: for the calculated porous media, consider the flow rate Q of fluid through the sample's ith layer i The sum of the pore flows of different sizes through the layer is equal, namely:
wherein: ΔP i The seepage pressure loss of the region corresponding to the ith synthesized binary image is l is the seepage path length of the region, mu is the dynamic viscosity and Q i To flow through the pores of the region, Q ik R is the flow rate through the different size pores of the region in the sample ik Is the pore percolation path radius of the region in the sample. Meanwhile, the flow of fluid also follows darcy's law, and the expression is:
wherein: k (k) i The permeability (i=1, 2, … …, N-1) of the region corresponding to the i-th synthesized binary image, and a is the cross-sectional area of the sample. Is obtained by the formulas (3) and (4):
wherein: r is R ik For the radius of each micro-pore of the region, S ik Is of radius R ik The pores occupy area. In practice, however, macropores are the main channels for fluid seepage, i.e. the influence of different sizes of pores on seepage is different. Thus leading to infiltrationThe concept of flow contribution values can be obtained:
wherein: s is S i Is of radius R i Is occupied by the pores of the substrate.
S3.3 starting from i=2, S3.1 and S3.2 are repeated to obtain the permeability of the remaining N-2 synthetic binary image.
And S4, calculating the permeability of the rock-soil body.
Obtained by the formula (5):
the overall rule of the rock-soil mass is as follows:
wherein: Δp is the total osmotic pressure loss of the rock-soil mass, l is the osmotic path length, μ is the dynamic viscosity, and Q is the flow through the rock-soil mass. Since the total pressure loss of the rock-soil body is equal to the sum of the pressure loss of all layers, the following steps are:
because the total flow Q of the rock and soil mass is equal to the flow Q of each layer i And, from formula (10):
the permeability k of the rock-soil mass is obtained as follows:
the invention has the beneficial effects that: (1) Compared with the prior art, the method has the advantages that the fusion reconstruction link of the adjacent binary images is improved, the form and the position of the pores in the original two adjacent binary images are considered, the obtained pores are more reasonable, and the permeability is prevented from being smaller due to smaller pores; (2) Compared with other permeability calculation methods based on three-dimensional microscopic images, the method has the advantages that the calculation time is remarkably reduced, and the method is better applied to large-size and multi-batch calculation tasks.
The foregoing is merely a preferred embodiment of the present application, but the scope of the present application is not limited thereto, and any changes or substitutions easily conceivable by those skilled in the art within the technical scope of the present application should be covered in the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.
Claims (5)
1. A rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction is characterized by comprising the following steps:
acquiring a three-dimensional microscopic image of a rock-soil body sample to be detected, and processing the three-dimensional microscopic image to acquire a pore-matrix binary image of the rock-soil body sample to be detected;
based on the pore-matrix binary images, acquiring pore communication components of adjacent pore-matrix binary images, and carrying out fusion reconstruction processing on pores in the adjacent pore-matrix binary images to acquire a synthesized binary image;
the method for acquiring the synthesized binary image comprises the following steps: based on the adjacent pore-matrix binary images, acquiring pore communication components of the adjacent pore-matrix binary images according to neighborhood 8 pixel communication; dividing pore pixels in adjacent pore-matrix binary images based on the pore communication components to obtain the synthesized binary image;
based on the pore communication component, a method of partitioning pore pixels in adjacent pore-matrix binary images includes: acquiring intersection, intersection edge pixels, union and union edge pixels of the pore connected components; acquiring the shortest distance between each pixel in the union and the intersection edge pixel, and recording the shortest distance as a first distance; acquiring the shortest distance between each pixel in the union and the edge pixel of the union, and recording the shortest distance as a second distance; acquiring an evaluation index according to the first distance and the second distance; dividing the pore pixels in adjacent pore-matrix binary images based on the evaluation index;
acquiring the permeability of the synthesized binary image based on the synthesized binary image;
and acquiring the permeability of the rock-soil body sample to be detected based on the permeability of the synthesized binary image.
2. The method for calculating the permeability of the rock-soil body based on the fusion reconstruction of the pores of the adjacent images according to claim 1, wherein the specific method for obtaining the first distance, the second distance and the evaluation index comprises the following steps:
tp u ＝d u2 /(d u1 +d u2 )
wherein d u1 At a first distance d u2 For a second distance tp u For evaluation index, x u And y u For the u-th pixel image array coordinates, x u1 And y u1 Is the closest intersection edge pixel coordinate with the u-th pixel, x u2 And y u2 Is the closest union edge pixel coordinate to the u-th pixel.
3. A method of rock-soil mass permeability calculation based on adjacent image pore fusion reconstruction as claimed in claim 1, wherein the method of obtaining the permeability of the composite binary image based on the composite binary image comprises:
according to the synthesized binary image, pore size distribution of pores in the synthesized binary image is obtained through a continuous pore size algorithm;
and obtaining the permeability of the synthesized binary image according to the Hagen-Poiseuille law and the pore diameter distribution of the pores.
4. A method for computing permeability of a rock-soil body based on adjacent image pore fusion reconstruction as claimed in claim 3, wherein the permeability of the composite binary image specifically comprises:
wherein k is i The permeability of the corresponding region of the binary image is synthesized for the ith sheet, i=1, 2, … … and N-1; a is the cross-sectional area of the sample, R ik For the radius of each micro-pore of the region, S ik Is of radius R ik The pores occupy area.
5. The method for calculating the permeability of the rock-soil body based on the fusion reconstruction of adjacent image pores according to claim 4, wherein the step of obtaining the permeability of the rock-soil body sample to be measured comprises the following steps:
wherein N is the number of pore-matrix binary images, k i The permeability of the corresponding region of the binary image is synthesized for the ith sheet, i=1, 2, … … and N-1; a is the cross-sectional area of the sample; r is R ik For the radius of each micro-pore of the region, S ik Is of radius R ik The pores occupy area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310317640.4A CN116342541B (en) | 2023-03-29 | 2023-03-29 | Rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310317640.4A CN116342541B (en) | 2023-03-29 | 2023-03-29 | Rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116342541A CN116342541A (en) | 2023-06-27 |
| CN116342541B true CN116342541B (en) | 2024-03-22 |
Family
ID=86883652
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310317640.4A Active CN116342541B (en) | 2023-03-29 | 2023-03-29 | Rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116342541B (en) |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101344474A (en) * | 2008-08-29 | 2009-01-14 | 中国农业大学 | Method for measuring porosity of soil |
| EP2103964A1 (en) * | 2008-03-21 | 2009-09-23 | Ifp | Method for calculating the density of the fractures in a rocky medium |
| CN103226089A (en) * | 2013-03-26 | 2013-07-31 | 中国石油天然气股份有限公司 | Method for testing shale gas permeability and shale gas permeability tester |
| KR20130096313A (en) * | 2010-12-25 | 2013-08-29 | 인텔 코오퍼레이션 | Systems, apparatuses, and methods for a hardware and software system to automatically decompose a program to multiple parallel threads |
| CN105136639A (en) * | 2015-06-24 | 2015-12-09 | 中国石油天然气股份有限公司 | Representation system of reservoir micro-pore structure |
| WO2017092409A1 (en) * | 2015-12-02 | 2017-06-08 | 华为技术有限公司 | Method and device for identifying high-usage intermediate code in language virtual machine |
| CN108133474A (en) * | 2017-12-22 | 2018-06-08 | 西安石油大学 | Permeability Prediction method based on core sample two dimension aperture image |
| CN108829950A (en) * | 2018-05-31 | 2018-11-16 | 中国科学院力学研究所 | A kind of unconventional reservoir permeability evaluation method based on core image |
| CN109191423A (en) * | 2018-07-18 | 2019-01-11 | 中国矿业大学 | A kind of porous media Permeability Prediction method based on machine image intelligence learning |
| CN110322427A (en) * | 2019-04-02 | 2019-10-11 | 中国矿业大学 | A kind of Rock And Soil calculation of permeability based on 3 dimension imaging technology |
| CA3062555A1 (en) * | 2018-07-18 | 2020-01-18 | China University Of Mining And Technology | Porous medium permeability prediction method based on machine image intelligent learning |
| CN114998246A (en) * | 2022-05-27 | 2022-09-02 | 安徽理工大学 | Method for calculating absolute permeability of porous medium based on image layering processing technology |
| CN115083543A (en) * | 2022-06-10 | 2022-09-20 | 深圳市国艺园林建设有限公司 | Method for predicting capillary water absorption performance of cement stone as super-large-volume mineral admixture |
| CN115239600A (en) * | 2022-08-25 | 2022-10-25 | 西南石油大学 | Shale scanning electron microscope image multi-component division method |
| WO2022227822A1 (en) * | 2021-04-29 | 2022-11-03 | 中国华能集团清洁能源技术研究院有限公司 | Simulation method and system for porous permeable rock |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2940441B1 (en) * | 2008-12-18 | 2010-12-24 | Inst Francais Du Petrole | METHOD FOR DETERMINING THE EVOLUTION OF PETROPHYSICAL PROPERTIES OF A ROCK DURING THE DIAGENESIS |
| CA2871781C (en) * | 2012-05-18 | 2017-02-07 | Ingrain, Inc. | Method and system for estimating rock properties from rock samples using digital rock physics imaging |
| US10416064B2 (en) * | 2015-12-14 | 2019-09-17 | Saudi Arabian Oil Company | Methods and systems for determining gas permeability of a subsurface formation |
| US11670073B2 (en) * | 2020-08-25 | 2023-06-06 | Saudi Arabian Oil Company | System and method for detection of carbonate core features from core images |
-
2023
- 2023-03-29 CN CN202310317640.4A patent/CN116342541B/en active Active
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2103964A1 (en) * | 2008-03-21 | 2009-09-23 | Ifp | Method for calculating the density of the fractures in a rocky medium |
| CN101344474A (en) * | 2008-08-29 | 2009-01-14 | 中国农业大学 | Method for measuring porosity of soil |
| KR20130096313A (en) * | 2010-12-25 | 2013-08-29 | 인텔 코오퍼레이션 | Systems, apparatuses, and methods for a hardware and software system to automatically decompose a program to multiple parallel threads |
| CN103226089A (en) * | 2013-03-26 | 2013-07-31 | 中国石油天然气股份有限公司 | Method for testing shale gas permeability and shale gas permeability tester |
| CN105136639A (en) * | 2015-06-24 | 2015-12-09 | 中国石油天然气股份有限公司 | Representation system of reservoir micro-pore structure |
| WO2017092409A1 (en) * | 2015-12-02 | 2017-06-08 | 华为技术有限公司 | Method and device for identifying high-usage intermediate code in language virtual machine |
| CN108133474A (en) * | 2017-12-22 | 2018-06-08 | 西安石油大学 | Permeability Prediction method based on core sample two dimension aperture image |
| CN108829950A (en) * | 2018-05-31 | 2018-11-16 | 中国科学院力学研究所 | A kind of unconventional reservoir permeability evaluation method based on core image |
| CN109191423A (en) * | 2018-07-18 | 2019-01-11 | 中国矿业大学 | A kind of porous media Permeability Prediction method based on machine image intelligence learning |
| CA3062555A1 (en) * | 2018-07-18 | 2020-01-18 | China University Of Mining And Technology | Porous medium permeability prediction method based on machine image intelligent learning |
| CN110322427A (en) * | 2019-04-02 | 2019-10-11 | 中国矿业大学 | A kind of Rock And Soil calculation of permeability based on 3 dimension imaging technology |
| WO2022227822A1 (en) * | 2021-04-29 | 2022-11-03 | 中国华能集团清洁能源技术研究院有限公司 | Simulation method and system for porous permeable rock |
| CN114998246A (en) * | 2022-05-27 | 2022-09-02 | 安徽理工大学 | Method for calculating absolute permeability of porous medium based on image layering processing technology |
| CN115083543A (en) * | 2022-06-10 | 2022-09-20 | 深圳市国艺园林建设有限公司 | Method for predicting capillary water absorption performance of cement stone as super-large-volume mineral admixture |
| CN115239600A (en) * | 2022-08-25 | 2022-10-25 | 西南石油大学 | Shale scanning electron microscope image multi-component division method |
Non-Patent Citations (1)
| Title |
|---|
| Image feature recognition and gas permeability prediction of Gaomiaozi bentonite based on digital images and machine learning;Jiangfeng Liu et al.;《Advances in Geo-Energy Research》;第6卷(第4期);314-323 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116342541A (en) | 2023-06-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Schlueter et al. | The fractal dimension of pores in sedimentary rocks and its influence on permeability | |
| CN104619952B (en) | Have reliable multiphase permeability and determines digital rock analysis system and method | |
| AU2018202368A1 (en) | Digital rock physics-based trend determination and usage for upscaling | |
| Lux et al. | 3D characterization of wood based fibrous materials: an application | |
| CN112446327A (en) | Remote sensing image target detection method based on non-anchor frame | |
| Parkinson et al. | Nanoimaging cells using soft X-ray tomography | |
| Dong et al. | Pore network modeling: analysis of pore size distribution of Arabian core samples | |
| CN104641259A (en) | Digital rock analysis with multiphase flow REV determination and method | |
| CN113609696B (en) | Multi-scale multi-component digital core construction method and system based on image fusion | |
| CN104885124A (en) | Method for producing a three-dimensional characteristic model of a porous material sample for analysis of permeability characteristics | |
| CN110322427B (en) | Rock-soil body permeability calculation method based on three-dimensional imaging technology | |
| CN111476159A (en) | Method and device for training and detecting detection model based on double-angle regression | |
| CN104685382A (en) | Digital rock analysis systems and methods that reliably predict a porosity-permeability trend | |
| US20140052420A1 (en) | Digital Rock Analysis Systems and Methods that Estimate a Maturity Level | |
| Salzer et al. | Quantitative comparison of segmentation algorithms for FIB‐SEM images of porous media | |
| CN108062789A (en) | Core sample choosing method and device | |
| Leue et al. | Separation of soil macropore types in three‐dimensional X‐ray computed tomography images based on pore geometry characteristics | |
| Solorzano et al. | Whole slide image registration for the study of tumor heterogeneity | |
| CN101980304A (en) | Three-dimensional digital volume image distortion measuring method | |
| Krutko et al. | A new approach to clastic rocks pore-scale topology reconstruction based on automatic thin-section images and CT scans analysis | |
| Gutiérrez‐Castorena et al. | Thematic micro‐maps of soil components using high‐resolution spatially referenced mosaics from whole soil thin sections and image analysis | |
| CN116342541B (en) | Rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction | |
| CN104634710B (en) | A kind of rock mass hole exceedes the analysis method and device oozed with connecting feature | |
| CN104809336A (en) | Method for sampling region factor by considering spatial correlation | |
| Bernard et al. | On the use of autocorrelation functions, permeability tensors, and computed tomography to characterise the anisotropy of diesel particulate filter materials |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |