CN115146215A - Multi-scale splicing method and system for micro-aperture data based on digital core - Google Patents
Multi-scale splicing method and system for micro-aperture data based on digital core Download PDFInfo
- Publication number
- CN115146215A CN115146215A CN202210953403.2A CN202210953403A CN115146215A CN 115146215 A CN115146215 A CN 115146215A CN 202210953403 A CN202210953403 A CN 202210953403A CN 115146215 A CN115146215 A CN 115146215A
- Authority
- CN
- China
- Prior art keywords
- pore
- data
- core
- scale
- volume
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 239000011148 porous material Substances 0.000 claims abstract description 140
- 238000002591 computed tomography Methods 0.000 claims abstract description 50
- 238000002474 experimental method Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 abstract description 7
- 238000011161 development Methods 0.000 abstract description 7
- 230000009977 dual effect Effects 0.000 abstract description 4
- 238000011156 evaluation Methods 0.000 abstract description 3
- 239000011435 rock Substances 0.000 description 22
- 239000003245 coal Substances 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 208000010392 Bone Fractures Diseases 0.000 description 2
- 206010017076 Fracture Diseases 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N2015/0846—Investigating permeability, pore-volume, or surface area of porous materials by use of radiation, e.g. transmitted or reflected light
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Dispersion Chemistry (AREA)
- Mathematical Optimization (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention discloses a method and a system for multi-scale splicing of micro-aperture data based on a digital core. The method comprises the following steps: respectively carrying out a double-resolution multi-scale CT scanning experiment on the core sample to obtain a CT core scanning data volume under double resolutions; selecting typical areas of CT core scanning data volumes under double resolution ratios, and analyzing pore data; calculating volume proportion coefficients of different pore data structures according to the selected typical region; and splicing the pore data acquired under different scans according to the volume proportion coefficient, thereby acquiring full-scale pore size distribution data. The invention aims to solve the problem of splicing multi-scale CT scanning data, so as to obtain comprehensive characterization of a compact core pore structure; the method provides powerful theoretical and technical support for fine and comprehensive evaluation of a compact core reservoir, particularly an unconventional reservoir with dual medium of pore gaps, and provides powerful power for efficient and accurate exploration and development of unconventional oil and gas.
Description
Technical Field
The invention relates to the field of scanning analysis of a micro-pore structure of a digital core, in particular to a micro-pore diameter data multi-scale splicing method based on the digital core.
Background
The digital rock core is an image of the rock core after X-ray scanning, and a simulated rock core is obtained through computer inversion, modeling and reconstruction, and the core modeling method can be divided into two categories: physical experimental methods and numerical reconstruction methods. In the physical experiment method, high-precision instruments such as a high-power optical microscope, a scanning electron microscope or a CT (computed tomography) imager are used for acquiring a plane image of the core, and then the plane image is subjected to three-dimensional reconstruction to obtain the digital core.
The CT is called Computed tomograph and is a Computed Tomography technology, the CT digital core plays more and more important roles in the aspect of rock core analysis in petroleum exploration and development of geological reservoirs, the physical principle is based on the interaction of rays and substances, the different densities in the core are different from the X-ray response degree, the rock framework is high in density and low in X-ray transmittance, the core pore density is small, the X-ray transmittance is high, a series of two-dimensional projection images are obtained by performing 360-degree omnibearing scanning on a sample, three-dimensional reconstruction is performed, a rock sample three-dimensional model is obtained, and the digital core is constructed. In the process of reconstructing the digital core, the physical property characterization information of the topological structure of the core pore and pore throat network is obtained, and the physical property characterization information generally comprises parameters such as core pore radius, throat radius, pore throat connectivity, porosity, surface porosity, fracture porosity/width, open porosity, closed porosity and the like, and the characteristics of the core micro-pore structure are evaluated through the parameters, so that the evaluation and characterization of the reservoir quality are more finely realized.
The existing digital core CT scanning analysis technology only carries out single data characterization, for example, 1 micron resolution CT scanning can only characterize the distribution of micropores (which mean 1-10 μm), and 13 micron resolution CT scanning can only characterize the distribution of macropores (which mean larger than 10 μm), so that the existing single-scale resolution data pore structure characterization has incomplete defects.
However, compact rock cores such as coal rock, shale, mudstone, volcanic rock and the like usually have multi-scale hole-crack dual medium characteristics, and have multi-scale holes and crack development in rock core reservoirs of the type, wherein the radius of the holes is mostly distributed in the range of 1-10 micrometers, and the width of the microcracks can reach dozens of even hundreds of micrometers. However, the resolution of CT scanning is in inverse proportion to the view of the slice, the view of the slice is only about 1mm under the high-precision 1-micron scanning resolution, and the structural information of macropores and microcracks cannot be effectively acquired; the low-precision 13-micron resolution scanning slice view can reach about 13mm, but the pore data with the pore diameter smaller than 13 mu m cannot be obtained. If the dense rock core is scanned twice with different resolutions, for example, a rock sample is scanned twice with a resolution of 1 micrometer and a resolution of 13 micrometers respectively, and then two sets of quantitative data of the pore structure under different scales are calculated respectively. Because the two groups of data respectively and independently represent pore crack sizes of different scales, the same sample can obtain two pore crack data conclusions with large differences, and the results can only represent the distribution characteristics of the pore diameters under the respective scales independently and cannot obtain the distribution rule of the pore diameters under the full scale. Therefore, unification of conclusions cannot be achieved during data interpretation, application and comparison among different samples, great difficulty is brought to popularization and application of the digital core technology in unconventional reservoirs, and the problem of splicing multi-scale data needs to be solved urgently.
Disclosure of Invention
The invention provides a digital core-based micro-aperture data multi-scale splicing method, which comprises the following steps of:
respectively carrying out a double-resolution multi-scale CT scanning experiment on the core sample to obtain a CT core scanning data volume under double resolutions;
selecting typical areas of CT core scanning data volumes under double resolution ratios, and analyzing pore data;
calculating volume proportion coefficients of different pore data structures according to the selected typical region;
and splicing the pore data acquired under different scans according to the volume ratio coefficient so as to acquire full-scale pore size distribution data.
According to the digital core-based micro-aperture data multi-scale splicing method, when a certain sample multi-scale pore structure is evaluated, a professional CT scanning instrument is used for carrying out 1-micron and 13-micron resolution multi-scale CT scanning experiments on the sample, and a 1-micron CT core scanning data body and a 13-micron CT core scanning data body are obtained.
The digital core-based micro-aperture data multi-scale splicing method comprises the steps that a 1-micrometer CT core scanning data body is used for identifying pore structures with the radius of 1-10 micrometers, and a 13-micrometer CT core scanning data body is used for identifying pore structures with the radius of more than 10 micrometers.
The method for splicing the data of the micro-pore diameters based on the digital core in the multi-scale mode comprises the steps of performing pore calculation by using a maximum sphere method, wherein the maximum sphere method is to fill a series of spheres with different sizes into pore spaces of the core, connection relations exist among the filled spheres with different sizes from large to small according to radiuses, the internal pore structure of the whole core is characterized by mutually overlapped and contained sphere strings, and the establishment of pores and throats in a pore network structure is completed by searching for the minimum sphere between a local maximum sphere and two maximum spheres in the sphere strings so as to form a pore-throat-pore pairing relation; finally, the whole ball string structure is simplified into a pore network structure model taking pores and throats as units, and corresponding pore throat quantitative data is obtained through the ball-stick model.
The method for multi-scale splicing of the data of the microporosity based on the digital core is characterized in that the volume ratio coefficient is obtained by dividing the volumes of the typical areas with two resolutions.
The digital core-based micro-pore diameter data multi-scale splicing method is characterized in that the pore frequency, the volume and the surface area ratio of the pore diameter in the range of (i to i + 1) mu m are calculated according to the volume proportionality coefficient, the number, the volume and the surface area value of the pore diameter in the range of (i to i + 1) mu m, the total number, the total volume and the total surface area value of the pore diameter in the range of 1 to 10 mu m, and the total number, the total volume and the total surface area value of the pore diameter larger than 10 mu m, wherein i =0,1,2,3,4,5,6,7,8,9.
The method for multi-scale splicing of the data of the micro pore diameter based on the digital core is characterized in that the frequency, the volume and the surface area ratio of pores with the pore diameter in the range of (m-10-m) mum are calculated according to a volume proportionality coefficient, the number, the volume and the surface area value of pores with the pore diameter of (m-10-m) mum, the total number, the total volume and the total surface area value of pores with the pore diameter of 1-10μm, and the total number, the total volume and the total surface area value of pores with the pore diameter of more than 10μm, wherein m =20,30,40,50,60,70,80,90,100.
The invention also provides a digital core-based micro-aperture data multi-scale splicing system, which comprises a CT scanning device and a data processing device, wherein the CT scanning device respectively performs double-resolution multi-scale CT scanning experiments on a core sample, and the data processing device executes any one of the above-mentioned digital core-based micro-aperture data multi-scale splicing methods.
The present invention also provides a computer storage medium, comprising: at least one memory and at least one processor;
the memory is to store one or more program instructions;
a processor configured to execute one or more program instructions to perform any one of the above methods for multi-scale stitching of digital core-based micro-pore size data.
The invention has the following beneficial effects: the invention aims to solve the problem of splicing multi-scale CT scanning data, so that the comprehensive characterization of a compact core pore structure is obtained. The method provides powerful theoretical and technical support for fine and comprehensive evaluation of a compact core reservoir, particularly an unconventional reservoir with dual medium of pore gaps, and provides powerful power for efficient and accurate exploration and development of unconventional oil and gas.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments described in the present invention, and other drawings can be obtained by those skilled in the art according to the drawings.
Fig. 1 is a flowchart of a multi-scale splicing method for micro-pore diameter data based on a digital core according to an embodiment of the present invention;
FIG. 2 is a 1-10 micron aperture distribution plot obtained from a 1 micron resolution CT scan;
FIG. 3 is a >10 micron aperture distribution map obtained from a 13 micron resolution CT scan;
FIG. 4 is a multi-scale tiling feature map of pore size distribution;
FIG. 5 is a gray scale diagram of the microcrack development characteristics of different coal samples obtained by CT scanning of samples A and B with different resolutions;
FIGS. 6 and 7 are typical profiles of the distribution before and after splicing of samples A and B.
Detailed Description
The technical solutions in the embodiments of the present invention are 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 some, not all, embodiments of the present invention. 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.
Example one
As shown in fig. 1, an embodiment of the present invention provides a digital core-based multi-scale pore diameter data splicing method, including:
when a multi-scale pore structure of a certain sample is evaluated, firstly, a professional CT scanning instrument is utilized to carry out a multi-scale CT scanning experiment with the resolution of 1 micrometer and 13 micrometers on the sample (a standard core plunger with the diameter of one inch is recommended to be adopted), 1 micrometer resolution CT scanning is carried out to obtain a core three-dimensional structure with the vision field size of 1mm x 1mm, the core three-dimensional structure can mainly identify the pore structure with the radius of 1-10 micrometers, and a CT core scanning data volume under the resolution of 1 micrometer is obtained; obtaining a rock core three-dimensional structure with a view field size of 13mm, 13mm by CT scanning with a resolution of 13 micrometers, wherein the rock core three-dimensional structure can mainly identify a pore structure with a radius larger than 10 micrometers, and obtaining a CT rock core scanning data volume under the resolution of 13 micrometers, for example, the CT rock core scanning data volume is a rock core projection image; the two sets of structural data can realize the dimensional complementation, which is the scanning resolution preferably adopted by the method.
120, selecting typical areas of the CT core scanning data volume under the double-resolution ratio, and analyzing pore data;
specifically, the areas selected for cores of different heterogeneity characteristics are of different sizes. To be of representative nature: selecting a larger area for the core slice with strong heterogeneity to ensure representativeness; for core slices with more uniform pore distribution, a smaller area can be properly selected to improve the analysis efficiency.
In the embodiment of the application, typical region selection is performed on a CT core scanning data volume under a double-resolution ratio, and the typical region selection method specifically comprises the following substeps:
step 121, slice reconstruction is carried out on CT core scanning projection under double resolution to obtain a two-dimensional gray slice;
preferably, professional projection data volume reconstruction software is used, a Feldkamp algorithm (Feldkamp algorithm) is adopted to perform three-dimensional reconstruction of scanning projection, the algorithm is formed by stacking cone beam rays as fan beam rays with different inclination angles along the Z-axis direction, data reconstruction on a central plane belongs to fan beam scanning accurate reconstruction, and for reconstruction of a non-central plane, fan beam reconstruction is corrected to complete accurate reconstruction.
Step 122, carrying out noise reduction on the two-dimensional gray slice by using a median filtering method;
the median filtering is a nonlinear smooth filtering technology, and can overcome the fuzzy problem caused by linear filtering under certain conditions. A two-dimensional gray scale slice may be considered as a set of two-dimensional gray value matrices, such as the following two-dimensional gray value matrices:
the set of gray values is sorted from small to large as 12, 14, 15, 16, 17, 18, 19, 20, where 16 is the median value, then the gray values in the region are all replaced by 16. By the method, each pixel value in each gray scale slice is very close to the surrounding (neighborhood) pixel value, so that the isolated abnormal gray scale noise points achieve a good denoising effect.
Step 123, selecting a typical area of the slice after the noise reduction treatment;
setting the typical area volume size of 1 micron resolution scanning gray scale sliceIs a V 1 Pixel is X 1 *Y 1 *Z 1 Typical area volume size of 13 micron resolution scanning gray scale slice is V 13 Pixel is X 13 *Y 13 *Z 13 In which V is 1 =1 3 X (number of pixels in X direction X 1 * Number of pixels Y in Y direction 1 * Number of pixels Z in Z direction 1 ),V 13 =13 3 X (number of pixels in X direction X 13 * Number of pixels Y in Y direction 13 * Number of pixels Z in Z direction 13 )。
Specifically, a maximum sphere method is used for analyzing pore data, wherein the maximum sphere method is used for filling a series of spheres with different sizes into a pore space of a three-dimensional core image, and connection relations exist among filling spheres with different sizes according to the radius from large to small. The internal pore structure of the whole core is characterized by mutually overlapped and contained ball strings, and the establishment of the pores and the throats in the pore network structure is completed by finding the minimum ball between the local maximum ball and the two maximum balls in the ball strings so as to form a pore-throat-pore pairing relation. Finally, simplifying the whole ball string structure into a pore network structure model taking pores and throats as units, and solving corresponding pore throat quantitative data through the ball-stick model;
for example, two sets of calculated pore size distribution data results are shown in tables 1 and 2 below:
TABLE 1 micrometer resolution scanning Aperture distribution data sheet
TABLE 2 data sheet of pore size distribution for 13 micron resolution scan
Pore size range mum | Number/number | Volume/. Mu.m3 | Surface area/. |
10~20 | 23547 | 394897518 | 52725153.25 |
20~30 | 17897 | 1039839956 | 115026790.7 |
30~40 | 4709 | 756541241.4 | 69016332.53 |
40~50 | 763 | 269211554.7 | 21268330.53 |
50~60 | 253 | 174186922.9 | 11995063.24 |
60~70 | 183 | 207502227.9 | 12946409.47 |
70~80 | 172 | 305400533 | 17590958.81 |
80~90 | 144 | 372962666 | 20040843.95 |
90~100 | 111 | 391979896.5 | 20047704.88 |
100+ | 172 | 1251446171 | 55666337.05 |
Total of | 47951 | 5163968687 | 396323924.4 |
FIG. 2 is a 1-10 micron aperture distribution plot obtained from a 1 micron resolution CT scan; figure 3 is a >10 micron aperture distribution map obtained from a 13 micron resolution CT scan.
Referring back to FIG. 1, step 130, calculating volume fraction factors for different pore data structures based on the selected representative region;
specifically, the volume proportion coefficient K is obtained by dividing the volumes of the two resolution typical regionsI.e. K = V 13 /V 1 In which V is 13 =13 3 X (number of pixels in X direction X 1 * Number of pixels Y in Y direction 1 * Number of pixels Z in Z direction 1 ),V 1 =1 3 X (number of pixels in X direction X 1 * Number of pixels Y in Y direction 1 * Number of pixels Z in Z direction 1 ). For example, the selected 13 micron resolution scan gray scale slice representative region has 658 × 684 × 671 and the 1 micron resolution scan gray scale slice representative region has 481 × 519 × 460, i.e., the pixels
the method comprises the steps of obtaining core pore distribution data information of more than 10 micrometers through 13-micrometer resolution CT scanning, obtaining 1-10-micrometer pore distribution through 1-micrometer resolution CT scanning, connecting data results of the two data, and representing distribution characteristics of a sample in a micrometer-scale full-pore-size range. Compared with conventional means such as mercury injection, nuclear magnetic resonance and the like, the multi-scale digital core technology reduces artificial damage of the porous pores in the experimental process, the obtained pore data are more direct and objective, and the distribution of the coal rock pores in the original state is reflected more truly.
In this application embodiment, carry out the concatenation operation with the pore data that obtains under the two resolution ratio scans, the data concatenation process needs the view of field normalization to same view of field size under two kinds of different resolution ratios, and concrete calculation method is as follows:
wherein f is i 、v i And s i Respectively, i =0,1,2,3,4,5,6,7,8,9, n i~i+1 、v i~i+1 、s i~i+1 Respectively representing the number, volume and surface area of pores with the pore diameter ranging from (i to i + 1) mu m, K representing the volume proportion coefficient of the two scale data volumes, N 1 Total 、V 1 total 、S 1 Total Respectively representing the total number of pores with the pore diameter of 1-10 mu m, the total volume and the total surface area value, N 13 total 、V 13 total 、S 13 total Respectively, the total number of pores with a pore diameter of more than 10 μm, the total volume and the total surface area value.
Wherein, F m 、V m And S m Representing the frequency, volume and surface area ratio of pores with a pore diameter in the range of (m-10 to m) μm, m =20,30,40,50,60,70,80,90,100, n m-10~m 、V m-10~m 、S m-10~m Respectively represent the number, volume and surface area of pores with a pore diameter (m-10 to m) mum.
If the 1 micron scanning data contains data larger than 10 microns, adding the part of data to N m-10~m 、V m-10~m And S m-10~m In the formula, the concrete formula is as follows:
the number of pores, the volume and the surface area ratio of each section can be calculated according to the formulas (1) to (9), namely:
……
……
the data splicing results of the volume V and the surface area S can be obtained in the same way, and are specifically shown in the following table 3:
TABLE 3 Multi-Scale stitching data results Table
FIG. 4 is a multi-scale stitching feature map of aperture distribution drawn from a multi-scale stitching data result.
The invention is applied to two different types of coal rock samples, and CT scanning experiments show that the sample A has obvious hole-seam dual media, and the other sample B has no micro-crack structure. As shown in fig. 5, the white and grey parts are the rock skeleton, and the black parts are the rock pores and cracks: it can be seen from the figure that the sample A develops not only micro cracks, but also micro cracks in both scales, the cracks have the characteristic of multiple scales, namely, micro cracks exist in both the micrometer cracks and the hundred-micrometer cracks, while the sample B does not develop micro cracks in both scanning scales, and the micro cracks and the sample B have obvious difference in the hole-slit structure.
FIGS. 6 and 7 are typical distribution maps of a sample A and a sample B before and after splicing according to the technical scheme of the invention. By utilizing the method for splicing two groups of data of two samples, the characteristic of pore radius distribution of micron full scale is obtained after splicing, and the pore size distribution has the typical characteristic of multiple peaks. Analysis shows that the first peak of the sample A after splicing mainly represents pores with the radius of 1-10 microns; the second peak represents microcracks with a radius of 10 to 30 microns; the third peak represents a hundred micron crack with a radius greater than 100 microns. The first peak of sample B after stitching represents predominantly pores with a radius of 1-10 microns, the second peak represents a seamless peak, and the third peak represents a seamless peak. The multi-scale combined map obtained by the method can well reflect the characteristics of multi-scale pore = fracture development of the coal rock sample, and the pore-size map with obvious characteristics is more favorable for comprehensively evaluating and classifying unconventional oil and gas reservoirs with complex structures and multi-scale, so that the application effect of the digital core technology in the field of geological exploration, particularly unconventional oil and gas exploration and development is greatly improved, and powerful theoretical support is provided for efficient exploration and exploitation of clean energy such as national coal bed gas, shale gas and the like.
In accordance with the embodiments described above, embodiments of the present invention provide a computer-readable storage medium having one or more program instructions embodied therein for execution by a processor to perform a method for multi-scale stitching of digital core-based microaperture data.
The disclosed embodiments provide a computer-readable storage medium having stored therein computer program instructions that, when executed on a computer, cause the computer to perform a method for multi-scale stitching of digital core-based micro-pore size data.
In an embodiment of the present invention, the processor may be an integrated circuit chip having signal processing capability. The Processor may be a general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other programmable logic device, discrete Gate or transistor logic device, discrete hardware component.
The various methods, steps, and logic blocks disclosed in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The processor reads the information in the storage medium and completes the steps of the method in combination with the hardware.
The storage medium may be a memory, for example, which may be volatile memory or nonvolatile memory, or which may include both volatile and nonvolatile memory.
The nonvolatile Memory may be a Read-Only Memory (ROM), a Programmable ROM (PROM), an Erasable PROM (EPROM), an Electrically Erasable PROM (EEPROM), or a flash Memory.
The volatile Memory may be a Random Access Memory (RAM) which serves as an external cache. By way of example, and not limitation, many forms of RAM are available, such as Static Random Access Memory (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), SLDRAM (SLDRAM), and Direct Rambus RAM (DRRAM).
The storage media described in connection with the embodiments of the invention are intended to comprise, without being limited to, these and any other suitable types of memory.
Those skilled in the art will appreciate that the functionality described in the present invention may be implemented in a combination of hardware and software in one or more of the examples described above. When software is applied, the corresponding functionality may be stored on or transmitted over as one or more instructions or code on a computer-readable medium. Computer-readable media includes both computer storage media and communication media including any medium that facilitates transfer of a computer program from one place to another. A storage media may be any available media that can be accessed by a general purpose or special purpose computer.
The above-mentioned embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above-mentioned embodiments are only examples of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made on the basis of the technical solutions of the present invention should be included in the scope of the present invention.
Claims (9)
1. A multi-scale splicing method of micro-aperture data based on a digital core is characterized by comprising the following steps:
respectively carrying out a double-resolution multi-scale CT scanning experiment on the core sample to obtain a CT core scanning data volume under double resolutions;
selecting typical areas of CT core scanning data volumes under double resolution ratios, and analyzing pore data;
calculating volume proportion coefficients of different pore data structures according to the selected typical region;
and splicing the pore data acquired under different scans according to the volume proportion coefficient, thereby acquiring full-scale pore size distribution data.
2. The digital core-based multi-scale pore diameter data splicing method as claimed in claim 1, wherein when a multi-scale pore structure of a certain sample is evaluated, a professional CT scanning instrument is used for carrying out 1 micron and 13 micron resolution multi-scale CT scanning experiments on the sample to obtain a 1 micron CT core scanning data body and a 13 micron CT core scanning data body.
3. The digital core-based multi-scale pore diameter data splicing method as claimed in claim 2, wherein a 1 micron CT core scan data body is used for identifying pore structures with a radius of 1-10 microns, and a 13 micron CT core scan data body is used for identifying pore structures with a radius of >10 microns.
4. The method for multi-scale splicing of the data of the microporosity based on the digital core according to claim 1, wherein the pore calculation is performed by using a maximum sphere method, wherein the maximum sphere method is to fill a series of spheres with different sizes into the pore space of the core, the filling spheres with different sizes have a connection relationship from large to small according to the radius, the internal pore structure of the whole core is characterized by mutually overlapped and contained sphere strings, and the establishment of the pores and the throat in the pore network structure is completed by searching the local maximum spheres and the minimum spheres between the two maximum spheres in the sphere strings so as to form a pore-throat-pore pairing relationship; finally, the whole ball string structure is simplified into a pore network structure model taking pores and throats as units, and corresponding pore throat quantitative data is obtained through the ball stick model.
5. The digital core-based multi-scale pore size data stitching method as claimed in claim 1, wherein the volume proportionality coefficient is obtained by dividing the typical area volumes of the two resolutions.
6. The method for multi-scale stitching of data of microporosity based on digital cores according to claim 1, wherein the pore frequency volume and the surface area ratio of the pore diameter in the range of (i to i + 1) μm are calculated according to the volume proportionality coefficient, the number, volume and surface area value of pores with the pore diameter in the range of (i to i + 1) μm, the total number, volume and total surface area value of pores with the pore diameter in the range of 1 to 10 μm, and the total number, volume and total surface area value of pores with the pore diameter of more than 10 μm, wherein i =0,1,2,3,4,5,6,7,8,9.
7. The method for multi-scale stitching of data of microporosity based on digital cores according to claim 1, wherein the frequency, volume and surface area ratio of pores with pore diameters in the range of (m-10 to m) μm is calculated according to a volume scaling factor, the number, volume and surface area values of pores with pore diameters in the range of (m-10 to m), and the total number, total volume and total surface area values of pores with pore diameters in the range of (1 to 10) μm, and the total number, total volume and total surface area values of pores with pore diameters greater than 10 μm, wherein m =20,30,40,50,60,70,80,90,100.
8. A multi-scale splicing system of micropore diameter data based on a digital core is characterized by comprising a CT scanning device and a data processing device, wherein the CT scanning device respectively carries out a double-resolution multi-scale CT scanning experiment on a core sample, and the data processing device executes the multi-scale splicing method of micropore diameter data based on the digital core according to any one of claims 1 to 7.
9. A computer storage medium, comprising: at least one memory and at least one processor;
the memory is used for storing one or more program instructions;
a processor for executing one or more program instructions to perform a digital core-based micro-pore diameter data multi-scale stitching method according to any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210953403.2A CN115146215A (en) | 2022-08-10 | 2022-08-10 | Multi-scale splicing method and system for micro-aperture data based on digital core |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210953403.2A CN115146215A (en) | 2022-08-10 | 2022-08-10 | Multi-scale splicing method and system for micro-aperture data based on digital core |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115146215A true CN115146215A (en) | 2022-10-04 |
Family
ID=83416187
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210953403.2A Pending CN115146215A (en) | 2022-08-10 | 2022-08-10 | Multi-scale splicing method and system for micro-aperture data based on digital core |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115146215A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116183875A (en) * | 2023-04-21 | 2023-05-30 | 煤炭科学研究总院有限公司 | Rock mass wave velocity determination method and device and electronic equipment |
CN117593299A (en) * | 2024-01-18 | 2024-02-23 | 北京大学 | Method, device, equipment and medium for evaluating space effectiveness of lamellar shale reservoir |
-
2022
- 2022-08-10 CN CN202210953403.2A patent/CN115146215A/en active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116183875A (en) * | 2023-04-21 | 2023-05-30 | 煤炭科学研究总院有限公司 | Rock mass wave velocity determination method and device and electronic equipment |
CN117593299A (en) * | 2024-01-18 | 2024-02-23 | 北京大学 | Method, device, equipment and medium for evaluating space effectiveness of lamellar shale reservoir |
CN117593299B (en) * | 2024-01-18 | 2024-04-19 | 北京大学 | Method, device, equipment and medium for evaluating space effectiveness of lamellar shale reservoir |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115146215A (en) | Multi-scale splicing method and system for micro-aperture data based on digital core | |
CN105115874B (en) | The multicomponent 3-dimensional digital rock core construction method of Multi-source Information Fusion | |
WO2020215524A1 (en) | Method for predicting coal porosity and permeability parameters based on fractal theory and ct scanning | |
US8155377B2 (en) | Method for determining rock physics relationships using computer tomographic images thereof | |
Lindquist et al. | Investigating 3D geometry of porous media from high resolution images | |
US8170799B2 (en) | Method for determining in-situ relationships between physical properties of a porous medium from a sample thereof | |
Alyafei et al. | Characterization of petrophysical properties using pore-network and lattice-Boltzmann modelling: Choice of method and image sub-volume size | |
AU2018202368A1 (en) | Digital rock physics-based trend determination and usage for upscaling | |
CN113609696B (en) | Multi-scale multi-component digital core construction method and system based on image fusion | |
CN109242985B (en) | Method for determining key parameters of pore structure from three-dimensional image | |
CN109345625B (en) | Rock core image self-adaptive partition three-dimensional reconstruction method | |
CN110320137A (en) | A kind of Multiscale Fusion method based on digital cores | |
Zhang et al. | GPU-accelerated 3D reconstruction of porous media using multiple-point statistics | |
CN109448104B (en) | Method and device for expanding and reconstructing based on core image | |
CN114565658A (en) | Pore size calculation method and device based on CT technology | |
CN109556542B (en) | CT size measurement method for complex lattice hollow structure | |
Kuo | Quantifying the fabric of granular materials: an image analysis approach | |
CN116012545B (en) | Multi-scale digital core modeling method, system, storage medium and application | |
CN112634428A (en) | Porous medium three-dimensional image reconstruction method based on bidirectional cycle generation network | |
CN114428040B (en) | Quantitative characterization and parameter acquisition method for shale oil reservoir storage and seepage space | |
Dong et al. | Verification of the carbonate double-porosity conductivity model based on digital cores | |
CN114705606A (en) | Blocking method of key seepage nodes in rock based on networked analysis | |
Zhang et al. | Fine characterization of pore structures in coral reef limestones based on three-dimensional geometrical reconstruction | |
CN113504172A (en) | Permeability evaluation method and device based on digital core | |
CN115147340A (en) | Rock crack deformation information acquisition method, device, equipment and storage medium |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20231115 Address after: 414, building 2, No. 1 jiaochangkou street, Xicheng District, Beijing 100032 Applicant after: Beijing Runze Innovation Technology Co.,Ltd. Applicant after: CHINA UNITED COALBED METHANE NATIONAL ENGINEERING RESEARCH CENTER Co.,Ltd. Address before: Room 414, Building 2, Yard 1, Jiaochangkou Street, Xicheng District, Beijing 100120 Applicant before: Beijing Runze Innovation Technology Co.,Ltd. |
|
TA01 | Transfer of patent application right |