3D characterization of porosity and minerals of low permeability uranium-bearing sandstone based on multi-resolution image fusion

3D characterization of porosity and minerals of low permeability uranium-bearing sandstone based on multi-resolution image fusion 1 first-author 1 2 corresp Corresponding author, usczengs@126.com Corresponding author, usczengs@126.com usczengs@126.com https://orcid.org/0000-0002-9283-6598 2 2 3 School of Civil Engineering, University of South China, Hengyang 421001, China School of Civil Engineering, University of South China, Hengyang 421001, China School of Resource Environment and Safety Engineering, University of South China, Hengyang 421001, China School of Resource Environment and Safety Engineering, University of South China, Hengyang 421001, China School of Resources and Safety Engineering, Central South University, Changsha 410083, China School of Resources and Safety Engineering, Central South University, Changsha 410083, China In the process of in situ leaching of uranium, the microstructure controls and influences the flow distribution, percolation characteristics, and reaction mechanism of lixivium in the pores of reservoir rocks, and directly affects the leaching of useful components. In this study, the pore throat, pore size distribution, and mineral composition of low-permeability uranium-bearing sandstone were quantitatively analyzed by high pressure mercury injection (HPMI), nuclear magnetic resonance (NMR), X-ray diffraction (XRD), and wavelength dispersive X-ray fluorescence (WD-XRF). The distribution characteristics of pores and minerals in the samples were qualitatively analyzed using energy dispersive scanning electron microscopy (EDS-SEM) and multi-resolution CT images. Image registration with the landmarks algorithm provided by FEI Avizo was used to accurately match the CT images with different resolutions. The multi-scale and multi-mineral digital core model of low-permeability uranium-bearing sandstone is reconstructed through pore segmentation and mineral segmentation of fusion core scanning images. The results show that the pore structure of low-permeability uranium-bearing sandstone is complex and has multi-scale and multi-crossing characteristics. The intergranular pores determine the main seepage channel in the pore space, and the secondary pores have poor connectivity with other pores. Pyrite and coffinite are isolated from the connected pores and surrounded by a large number of clay minerals and ankerite cements, which increases the difficulty of uranium leaching. Clays and a large amount of ankerite cement filled in the primary and secondary pores and pore-throats of the low-permeability uranium-bearing sandstone, which significantly reduces the porosity of the movable fluid and results in low overall permeability of the cores. The multi-scale and multi-mineral digital core proposed in this study provides a basis for characterizing macroscopic and microscopic pore-throat structures and mineral distributions of low-permeability uranium-bearing sandstone, and can better understand the seepage characteristics. Key words: Low permeability uranium-bearing sandstone Digital core Micro-CT SEM-EDS Image fusion 1. Introduction 1 The in-situ leaching process of sandstone uranium ore can be regarded as the flow process of lixivium in the connected pores of reservoir rocks. The uranium stored in the rock minerals is dissolved into the lixivium under the action of a chemical solution, and then migrates out from the surface of the rock pores with the leachate [ 1 1 , 2 2 ]. The mineral components in uranium bearing sandstones provide the material basis for pore formation and evolution, and directly control the quality of reservoirs [ 3 3 , 4 4 ]. The microstructure of uranium ore rocks and the connectivity of pores have an important impact on lixivium transport. Microscopic study of the mineral and pore distribution characteristics of uranium bearing sandstone is of great significance. A large number of successful studies of microstructure of reservoir rocks have been performed worldwide. The research methods can be summarized into three categories: experimental data analysis technology, image analysis technology, and digital core technology [ 5 5 ]. The experimental data analysis techniques for the microscopic pore structure of reservoir rock mainly include high pressure mercury injection (HPMI), constant speed mercury injection (CRMI), nuclear magnetic resonance (NMR), gas adsorption, and neutron small angle scattering (SANS) experiments. The throat size of a certain level of rock and the distribution of the pore throat volume controlled by the corresponding throat can be obtained through the HPMI experiment, which is the main method to measure the capillary pressure of rocks. However, due to the limitation of mercury injection pressure, it cannot fully reflect the reservoir pore space [ 6 6 , 7 7 ]. NMR transforms the T 2 spectrum of nuclear magnetic resonance into the distribution curve of pore radius, which can be used to evaluate the pore structure of cores and characterize the pore size distribution at multi-scales, but it cannot obtain the seepage capacity of reservoir rocks [ 8 8 , 9 9 ]. X-ray diffraction (XRD) and X-ray fluorescence (XRF) tests are conventional quantitative analysis methods for minerals in geological research [ 10 10 , 11 11 ]. However, these methods cannot analyze the morphology characteristics of minerals and only quantitatively analyze their content. Image analysis techniques collect images of core sections through high-resolution instruments, and facilitate qualitative analysis of the mineral morphology, and pore and throat size distribution and connectivity in rocks. The image analysis techniques of reservoir rock microstructure include cast thin sections, scanning electron microscopy (SEM), micro/nano-CT scanning, and focused ion beam SEM (FIB-SEM). As a low cost two-dimensional (2D) imaging technology with high resolution that can detect nano pores, SEM is widely used in the analysis of microscopic pore structure and mineral characteristics of reservoir rocks [ 12 12 , 13 13 ]. The SEM images combined with energy spectrometer, elemental scanning analyzer, and cathode luminescence measurements, can identify mineral types and generate mineral distribution diagrams. By means of energy dispersive scanning electron microscopy (EDS-SEM), the morphological characteristics, pore types, genesis, output, and occurrence of various minerals on the microscopic scale can be studied [ 14 14 - 16 16 14-16 ]. Compared with the cast thin sections and SEM, micro/nano-CT scanning has the advantages of non-destructive measurement, intuitive measurement results, and accurate description of spatial characteristics of microstructure, which is widely used in the construction of various rock digital cores and pore network models [ 17 17 - 19 19 17-19 ]. FIB-SEM can observe the distribution and connectivity of nanopores of unconventional tight rocks such as natural gas shale and fine sandstone more clearly. However, it is not widely used owing to the high cost of sample preparation and limited scope of observation [ 20 20 , 21 21 ]. Digital core technologies obtain the microstructure characteristics of rocks without destroying the real cores by using various micro-imaging techniques. On this basis, various algorithms are used to construct three-dimensional (3D) digital core models of rocks [ 22 22 , 23 23 ]. The combination of the reconstructed digital core and multi-field coupling finite element technologies is used widely for fluid seepage simulation of rock–fluid systems, which can realize data visualization at the pore level and deepen the understanding of various flow mechanisms [ 24 24 - 29 29 24-29 ]. Digital petrophysics experiments can not only be verified by laboratory experiments, but can also explain petrophysics and chemistry experiments from a microscopic perspective, which can replace laboratory experiments when appropriate samples cannot be obtained [ 30 30 - 32 32 30-32 ]. Numerical simulation of petrophysics based on digital core technologies has become an important means of petrophysics research. However, the core size is small and the representative volume cannot be obtained when the resolution of the micro-CT scan is high. In addition, micro-CT scanning has the disadvantage of not characterizing the mineral composition information, which is very important for physicochemical seepage simulation of reservoir rocks [ 33 33 ]. Therefore, it is necessary to describe the core structure using several methods. Image fusion refers to the process in which images of the same scene obtained by different instruments or the same scene obtained by the same instrument at different times are de-noised, time registered, spatially registered, or resampled by some mathematical algorithm to obtain a higher quality synthetic image [ 34 34 , 35 35 ]. It can overcome the existence of single sensor image in geometry, spectrum, and spatial resolution limitations and differences. Multi-resolution image fusion for reservoir microstructure analysis is currently a hot research topic [ 18 18 ]. For example, Markussen et al. [ 36 36 ] analyzed the pore distribution and the connectivity between pore and throat by using dry/wet difference mapping CT scan method, and discussed the influence of carbonate cementation on porosity and pore-throat characteristics. According to the registration of 2D to 3D images, the nano-scale data in 2D slices can be directly related to the micro-scale features in 3D images [ 37 37 , 38 38 ]. Combining micro-CT with SEM images can improve information on microporosity and microstructure of analyzed samples, and the mineral information in CT image can be recognized to allow the visualization of spatial distribution of different mineral phases based on automated mineralogy mapping by SEM-EDS [ 14 14 , 39 39 ]. Liu et al. [ 40 40 ] established multi-mineral digital core of tight sandstone by using SEM-EDS combined with CT images, and analyzed different pore types and their connectivity in tight sandstone by using expansion algorithm. Based on the superposition method, Lin et al. [ 41 41 ] reconstructed the double pore 3D digital cores of carbonate rocks by using micro-CT combined with SEM images and analyzed its seepage characteristics by Lattice Boltzmann method. The low permeability of uranium ores is one of the key factors restricting the successful production of in-situ leaching of sandstone-type uranium deposits. Low-permeability uranium-bearing sandstone forms different pore types under geological and sedimentary processes, and the pore structure is complex [ 42 42 , 43 43 ]. At the same time, its internal mineral composition changes greatly, while different mineral compositions and structures have different physicochemical properties, which have an important influence on the seepage of lixivium [ 44 44 ]. At present, studies on in-situ leaching solute transport in uranium ore are rarely carried out from a micro perspective [ 45 45 ]. In this study, the microstructure characteristics of low-permeability uranium-bearing sandstone were studied using comprehensive testing methods. The low permeability uranium-bearing sandstone cores were scanned by medium- and high resolution micro-CT, and the core surface was tested by EDS-SEM. The feature-based image registration method is used to accurately match CT images with different resolutions in the Avizo software. The low-permeability uranium-bearing digital core model was established after pore segmentation and mineral segmentation on fusion CT images by K-means clustering segmentation and the watershed segmentation method combined with EDS-SEM images. It is expected to provide a deep understanding of the microstructural characteristics of low-permeability uranium-bearing sandstone and provide a model basis for the subsequent study of the influence of pore structure and mineral distribution on the seepage characteristics of lixivium. 2. Samples and experimental methods 1 2.1 Sample source and collection 2 The Monchgur uranium deposit is from the interlayer oxidation zone of the Middle-lower Jurassic Shuixigou group zone sandstone type rolled uranium deposit located in the southern margin of the Yili Basin, Xinjiang, China. The hydraulic conductivity of the entire ore-bearing aquifer is less than 0.5 m/d, which belongs to the low permeability ore body type. The uranium deposit has a special structure and hydrodynamic conditions. At the same time, the sand bodies of the Sangonghe Formation are the best developed in the southern margin of the Yili Basin with good lithofacies and lithologic conditions, which is conducive to the infiltration of oxidation fluid and the formation of uranium ore bodies [ 2 2 ]. Therefore, two representative sandstone samples from the upper member of the Sangonghe Formation are selected in this study. The selected samples belong to gray medium sandstone, which is mainly composed of medium-grained sand chips with loose structure. After oxidation and hydrolysis, it is in the shape of "bean curd residue"( Figs. 1 Figs. 1 , 2 2 , 3 3 ; Table 1 Table 1 ). *Porosity was measured by the MIP method 1 1 Sample number 1 1 Well number 1 1 Sampling depth(m) 1 1 Porosity*(%) 1 1 Lithology 1 1 S1 1 1 J5 1 1 429.3 1 1 16.39 1 1 Grayish white medium sandstone 1 1 S7 1 1 J9 1 1 425.9 1 1 12.38 1 1 Grayish white medium sandstone 2.2 Experimental methods 2 2.2.1 HMPI testing 3 HMIP is one of the most commonly used experimental techniques to measure pore throat size distribution and effective porosity. Based on the mercury injection method, the lower limit of fluid flow, throat, and its connected pore volume can be quantitatively evaluated, which is helpful for understanding the law of fluid seepage [ 47 47 ]. An Autopore IV 9500 pore analyzer was used in this HPMI experiment. The test pressure ranges from 0.1 Psia to 33000 Psia, and the measured pore throat diameter ranges from 360 μm to 5 nm. The pore throats are assumed to be connected capillary tubes of different diameters, and the permeability of the pore throat of radius r j is [ 48 48 , 49 49 ] <math display="block" id="M1"><mrow><msub><mi>k</mi><mtext>j</mtext></msub><mo>=</mo><mfrac><mrow><msub><mi>ϕ</mi><mtext>j</mtext></msub></mrow><mn>8</mn></mfrac><msubsup><mi>r</mi><mtext>j</mtext><mtext>2</mtext></msubsup><mo>,</mo></mrow></math> where k j represents the permeability of the pore throat with radius r j , and ∅ j is the ratio of the volume of the pore throat with radius r j to the total pore volume. The core permeability contribution rate corresponding to a certain radius of the pore throat is characterized by <math display="block" id="M2"><mrow><mfrac><mrow><msub><mi>k</mi><mi>j</mi></msub></mrow><mi>k</mi></mfrac><mo>=</mo><msub><mi>ϕ</mi><mi>j</mi></msub><msubsup><mi>r</mi><mi>j</mi><mn>2</mn></msubsup><mo>/</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mi>i</mi><mi>N</mi></munderover><mrow><msub><mi>ϕ</mi><mi>i</mi></msub></mrow></mstyle><msubsup><mi>r</mi><mi>i</mi><mn>2</mn></msubsup><mo>.</mo></mrow></math> 2.2.2 NMR testing 3 NMR experiments transform the T 2 spectrum of nuclear magnetic resonance into the distribution curve of the pore radius, which can be used to evaluate the pore structure of cores and characterize the pore size distribution at multiple scales [ 50 50 ]. NMR experiments were carried out on the MacroMR12-150H-I core nondestructive and microcosmic testing imaging and analysis system. 2.2.3 XRD testing 3 XRD can obtain information on the mineral composition, crystal structure, or morphology, through the diffraction patterns of samples, so as to quickly and accurately identify the minerals of ore rocks. A representative portion of the samples (approximately 10 g) was crushed and ground into powder (particle size < 0.2 mm), and the main minerals in the sample were analyzed using a Bruker D8 Advance X-ray diffractometer. The Jade6.0 software was used to analyze the phases of the XRD patterns. After peak division, background value deduction, and smoothing treatment, combined with the geological background of the study area, the phase compositions were retrieved using standard PDF reference cards. 2.2.4 WD-XRF testing 3 Wavelength dispersive X-ray fluorescence spectrometer is a type of spectrometer that uses crystal spectrum. It is widely used in the determination of various elements and oxides in geological studies. The WD-XRF spectrometer Axios MAX was used for the quantitative analysis of mineral elements in solid disk samples with a diameter of 38 mm and a thickness of 10 mm. 2.2.5 SEM-EDS experiment 3 SEM-EDS can automatically identify the mineralogy of each image pixel according to the detected element spectrum with the known chemical compositions and generate a mineral distribution map [ 31 31 , 40 40 ]. The surfaces of the samples were gently polished with 1000, 2000, and 4000 grade polishing paper, and then polished to 0.5 μm in a silicon solution. The samples were ultrasonically cleaned in water to remove any silica-like material during polishing and drying. In order to improve the conductivity of the samples under an electron beam, several nanometers of amorphous carbon was sputtered on the surface of the samples. SEM-EDS experiments were carried out on an FEI Scios™ system, which is equipped with secondary and backscatter electron detectors and silicon drift detectors of the Oxford instrument for EDS. By adjusting the scanning voltage, it was found that the contrast of the SEM images was good at 10 kV, so 10 kV was selected to scan the end face of the sandstone sample in this study. The spot size of the scanning electron microscope image was 10 nm and the working distance was 6.8 mm. EDS mapping images of the same surface were also obtained using AZtec software. 2.2.6 Multi-resolution CT scans 3 Micro-CT scans are based on the different X-ray absorption coefficients of different density components in the rock, and the computer image reconstruction is used to make the material information of different densities on the designated layer displayed in high-resolution digital images. The low permeability uranium-bearing sandstone used in this study is easily damaged when sampling with a drill bit because of its low mechanical strength, so the rock sample is embedded in epoxy resin to ensure the strength. Then, samples of 15 mm × 15 mm × 20 mm and 1.5 mm - 1 mm microplugs were drilled by bits of different sizes. Micro-CT scans were performed on the nanoVoxel-2000 X-ray CT scanner system produced by Sanying Precision Instrument Co., Ltd. (Tianjin, China), which has a maximum resolution of 0.5 μm. The basic sample size used for the entire study were rectangular samples of 15 mm × 15 mm × 20 mm. First, the basic samples were scanned by micro-CT with medium resolution (12.58 μm) to observe the overall structure and internal porosity. According to the medium-resolution CT scan results, the volume of interest (VOI) was then selected for further imaging. The locations of the subsamples (1.5 mm - 1 mm) were marked and extracted from the basic cores and scanned at a higher resolution (1 μm) to study the finer pore structure. Table 2 Table 2 shows the settings of the scanning parameters of the micro-CT equipment at different resolutions. 1 1 Sample 1 1 Dimension(mm) 1 1 Resolution(μm) 1 1 Voltage(kV) 1 1 Electric current(uA) 1 1 Exposure time(s) 1 1 Model 1 1 15 mm×15 mm×20 mm 1 1 12.58 1 1 150 1 1 150 1 1 0.65 1 1 Reference 1 1 1.5 mm - 1 mm 1 1 1 1 1 80 1 1 60 1 1 0.7 3. Results 1 3.1 Mineral characteristics of low-permeability uranium-bearing sandstone 2 3.1.1 Mineral composition and chemical constituents 3 XRD combined with XRF analysis can improve the accuracy of mineral composition identification. XRD analysis of the samples ( Fig. 4 Fig. 4 ) shows that the main minerals in the samples are quartz, ankerite, feldspar, and calcite, followed by clay minerals such as kaolin, illite, and chlorite, and a small amount of metal minerals (owing to the small content, it is not shown in the figure), such as pyrite, rutile, and sphalerite. Clay minerals mainly include kaolinite, illite, and chlorite. Kaolinite is the most developed, followed by chlorite and illite, which reflects the dissolution of feldspar during diagenesis and mineralization. The results of petrochemistry analysis of low-permeability uranium-bearing sandstone are shown in Table 3 Table 3 . The K 2 O+Na 2 O content of the ore bearing rock samples ranges from 2.88 % to 3.30 %, and the K 2 O/Na 2 O value is greater than 1, which indicates that the feldspar clasts in the sandstone are mainly K-feldspar. 1 1 Compound 2 1 Content (%) 1 1 Compound 2 1 Content (%) 1 1 1 1 S1 1 1 S2 1 1 1 1 S1 1 1 S2 1 1 SiO 2 1 1 45.047 1 1 53.865 1 1 Cl 1 1 0.041 1 1 0.039 1 1 CaO 1 1 25.178 1 1 19.387 1 1 U 1 1 0.036 1 1 0.031 1 1 Al 2 O 3 1 1 18.623 1 1 16.495 1 1 Mo 1 1 0.019 1 1 0.008 1 1 Fe 2 O 3 1 1 6.875 1 1 4.110 1 1 Tb 4 O 7 1 1 0.025 1 1 0.030 1 1 K 2 O 1 1 3.290 1 1 2.759 1 1 Cr 2 O 3 1 1 0.022 1 1 0.032 1 1 MnO 1 1 0.552 1 1 0.714 1 1 Na 2 O 1 1 0.0141 1 1 0.019 1 1 MgO 1 1 0.531 1 1 0.881 1 1 SrO 1 1 0.013 1 1 0.016 1 1 SO 3 1 1 0.248 1 1 0.619 1 1 Rb 2 O 1 1 0.013 1 1 0.014 1 1 TiO 2 1 1 0.239 1 1 0.398 1 1 ZrO 2 1 1 0.013 1 1 0.017 1 1 P 2 O 5 1 1 0.214 1 1 0.491 1 1 ZnO 1 1 0.011 1 1 0.013 1 1 CeO 2 1 1 0.046 1 1 0.057 1 1 Y 2 O 3 1 1 0.003 1 1 0.004 3.1.2 Mineral distribution characteristics 3 The mineral and pore distribution images of the polished faces were obtained using EDS-SEM analysis. EDS mapping analysis showed the distribution of Si, Al, Ca, Mg, O, C, K, and S ( Fig. 5 Fig. 5 ). It can be seen that the main intergranular pores in the rock sample are mainly filled with clays (Al+Si+Mg), while quartz (Si+O), dolomite (Ca+Mg+C), K-feldspar (Si+Al+K), and other particles that constitute the rock sample skeleton are covered with many clays. According to the mineralogical composition and mass fraction in Fig. 5 Fig. 5 , it can be seen that the relative content of ankerite and dolomite is high, while the relative content of calcite is low, which indicates that there is significant dissolution and metasomatism in the diagenesis of low permeability uranium-bearing sandstone reservoir [ 42 42 ]. 3.1.3 Occurrence state of uranium 3 Under the SEM, uranium minerals are mostly bright white, and the particles are small and scattered in the rock minerals. In addition, coffinite often coexists with pyrite and disperses in sandstone cements in clumps ( Fig. 6a Fig. 6a ). The BSE image shows that the irregular pyrite grains are mainly distributed around clastic quartz and rock debris, while star point pyrite is mainly distributed in cement and clay minerals ( Fig. 6a Fig. 6a ). Due to the strong adsorption of clay minerals, coffinite is absorbed by clay minerals and then precipitated around the clay membrane or in clay minerals ( Fig. 6b Fig. 6b ). Coffinite also exists on the surface of pyrite and in intergranular pores ( Fig. 6c Fig. 6c ). Moreover, some intact coffinite occurs in the middle of clastic grains such as feldspar, which should be carried from the provenance area ( Fig. 6d Fig. 6d ). 3.2 Pore structure characteristics of low-permeability uranium-bearing sandstone 2 3.2.1 Pore structure types 3 The pore types and genesis of low-permeability uranium-bearing sandstone were analyzed by micro-CT scanning and SEM. It can be seen that there are four main pore types in the cores of low permeability uranium-bearing sandstone, which include intergranular pores, intragranular dissolution pores, intercrystalline micropores, and micro fractures, as shown in Fig. 7 Fig. 7 . Among them, intergranular pores mainly include intergranular pores of detrital minerals and intergranular dissolution pores. The number of intergranular pores of clastic minerals is small, but the spatial scale is relatively large, and most of them are in the shape of narrow triangles or multilateral shapes with smooth straight edges ( Fig. 7a Fig. 7a ). The intergranular pores of clastic minerals were altered by diagenesis, such as mechanical compaction or cemented filling of clay minerals, resulting in a large number of micro-nano scale intergranular dissolution pores ( Fig. 7b Fig. 7b ). The dissolution pores are mainly secondary pores produced by carbonate particles, feldspar particles, and interstitials (kaolinite, chlorite) formed by dissolution. Dissolution pores within the grains are widely developed and have irregular pore surfaces with diameters ranging from a few nanometers to a few microns ( Fig. 7c, d, and e Fig. 7c, d, and e ). Intergranular micropores are mostly formed by authigenic minerals and clay minerals in rocks ( Fig. 7f, g Fig. 7f, g ) with small pore diameters. Micro fractures mainly exist in carbonate minerals and are formed by dissolution as well as brittle mineral fracture and clay mineral shrinkage ( Fig. 7h, i Fig. 7h, i ) with a width of less than ten nm and up to several hundred microns. The existence of micro fractures enhances the connectivity of reservoir pores, increases the permeability, and plays an important role in the drainage of the solution. In summary, it can be seen that in the samples, even the same mineral may contain multiple pore types. The pore structure of low-permeability uranium-bearing sandstone is complex and has multi-scale characteristics. 3.2.2 Pore throat size distribution analysis by MIP 3 The MICP data show that the pore throat sizes of low-permeability uranium-bearing sandstone vary from millimeter to micrometer. According to the diameter, the term ''micro pore throat’’ is applied to pore throat sizes less than 1 μm, ''meso pore throat’’ is the term applied to pore throat sizes less than 10 μm, and ''macro pore throat’’ is the term applied to pore throat sizes greater than 10 μm [ 51 51 , 52 52 ]. Based on the statistics of the influence of different pore-throat sizes on the permeability of the samples, it was found that the contribution of permeability of the samples mainly comes from the large throat (99.51 %), and it can be inferred that the large pore throat has good connectivity. 3.2.3 Pore size distribution analysis by NMR 3 The NMR test results of the samples were converted into pore diameter distributions, as shown in Fig. 9 Fig. 9 . The pore diameter of the samples is in the range of 50–300 μm. Although the seepage ability of micropores is small, the pore volumes of the nanopores and micropores are larger. The micropores are not easily wetted and permeated by the lixivium, which limits the free flow and permeability of the lixivium in its interior. Moreover, the influence of clay minerals on porosity and pore flow capacity cannot be ignored when clay minerals disperse into fine particles along with the flow of lixivium and block micropores. 3.3 Construction and evaluation of low-permeability uranium-bearing sandstone digital core 2 In this study, the CT data were analyzed using the 3D visualization software Avizo 9.0.3. First, the contrast of the CT images is adjusted and the signal-to-noise ratio is boosted by a non-local means filter, so as to optimize the images into a mode that is more suitable for identification and analysis. According to the X-ray attenuation characteristics of different minerals, the mineral phases whose gray values differ greatly can be directly identified from the tomographic images [ 14 14 , 53 53 ]. For example, the high density pyrite is bright white in the CT image. The clay minerals, quartz, and feldspar can be identified because of the large difference in gray values. It can be seen that the medium-resolution CT scan is helpful for obtaining a basic understanding of the internal structure of rocks over a large range ( Fig. 10 a Fig. 10 a ). For example, there is a large amount of carbonate cements and clay minerals that may have a great influence on the total porosity and permeability in low-permeability uranium-bearing sandstone. High-resolution CT images can identify micropores and obtain better image quality ( Fig. 10 b Fig. 10 b ). At the same time, the correlation between quartz particles, feldspar, and clay minerals is significant in the high-resolution images, which is conducive to separating minerals from the skeleton. In order to establish a digital core model that can accurately reflect the pore structure characteristics of reservoir rocks, two types of resolution CT images can be registered in space to find the corresponding position of the subsample image in the plunger sample image and form multi-scale core scanning images. 3.3.1 Image registration 3 Image registration is a necessary preliminary procedure for multi-image fusion. The feature-based image registration method uses the feature points of the image to achieve image registration. Firstly, the interested image features are detected by the feature extraction operator. Secondly, the feature matching algorithm is used to select the feature pairs with matching relationships in the corresponding feature set of two images. Then the transformation model between images is calculated by using the coordinate relationship between pairs feature parameters to achieve image registration [ 54 54 ]. The main goal of spatial registration of different resolution 3D CT images is to align and overlap the two-image data into the same coordinate system. In this study, registration with the landmarks algorithm provided by FEI Avizo was used to register CT images of low-permeability uranium-bearing sandstone with different resolutions. In registration, typically a low-resolution CT image dataset is taken as the reference, and the model of high-resolution CT image data sets is transformed, moved, and possibly rescaled, until both data sets match. According to the location of the pre-marked subsample, the converter of the Transform Editor was used to translate, rotate, and scale the reference image, so that the model image was roughly aligned with the reference image. Landmark sets were then artificially selected and created to calculate the transformation equation matching the landmark sets. Using the Tcl command in Avizo console, the numerical values of the registration transformation can be retrieved. Owing to the distortion of the image in the process of sampling, the simple rigid body transform algorithm cannot satisfy the registration of the images. Therefore, the Bookstein image distortion algorithm was used to register the images fine. The points in the transformed group (model) were forced to fit to the reference points, and all image points in between were shifted by using the correlation of the surrounding points to achieve the best image registration [ 14 14 , 55 55 ]. After verification and confirmation of overlapping known features, the registration is successful, and features (texture, crack, void) with the same coordinates are found on the final registered image ( Fig. 11 c, d Fig. 11 c, d ). 3.3.2 Pore Segmentation 3 After fusion, CT images have both high- and low-resolution image data. For medium-resolution CT images, owing to the serious marginal fuzzy effect of pores and skeletons, they cannot be simply divided by the binary method [ 56 56 ]. The K-means method is a purely stochastic technique that estimates G min and G max based only on the gray-level distribution of the 3D microtomography images [ 24 24 , 57 57 ]. Using the two threshold levels G min and G max estimated by the K-means algorithm, the gray-level of the CT image can be divided into three segments: pores, skeleton, and unresolved pores. The linear function of the porosity and probability distribution function of the grayscale of CT images with different resolutions are shown in Fig. 12 Fig. 12 . By calculating the G min and G max of the high-resolution voxel corresponding to the low-resolution voxel, and changing the gray value of G min of the low-resolution CT image, although the pores and skeletons cannot be completely accurately divided, the overlapped gray levels in the low-resolution CT image can be better segmented and the pore characteristics of diameter greater than 1 μm in the high-resolution CT image can be fused [ 16 16 ]. Indeed, the contribution of the macro throat ( > 10 μm) to the permeability of the sample in this test is dominant, and that of the micro throat ( < 1 μm) can almost be ignored, so the digital core pore model established after integrating the pores of the high-resolution CT image meets the reliability of the next permeability simulation analysis. Finally, the three-dimensional digital core pore model of the plunger sample is obtained, as shown in Fig. 13 Fig. 13 . It can be seen from Fig. 8b Fig. 8b that the pore connectivity in the medium-resolution CT image has been significantly improved after the fusion of the pores identified by the high resolution. 3.3.3. Skeleton Segmentation 3 Owing to the blurring of edges caused by noise and volume effects in the process of image acquisition and reconstruction, the commonly used simple threshold segmentation method often produces inaccurate or even wrong segmentation results, especially in the phase transition region where more than two phases are segmented. Watershed segmentation can effectively solve these problems in many cases. The principle and operation of watershed segmentation are introduced in detail in [ 36 36 , 46 46 ], which will not be repeated here. However, watershed segmentation is prone to undersegmentation or oversegmentation, and the segmentation results usually have some errors with the measured values. BSE images have the advantages of high resolution, clear microscopic component morphology, and can be combined with energy spectrum analysis (EDS) to identify minerals. At the same time, the gray value of the SEM image pixels is related to the density of rock samples, so it can be used to verify the distribution of gray values of different mineral components in CT images [ 38 38 , 58 58 , 59 59 ]. In this study, SEM images, EDS spectra, and maps were used to calibrate the threshold segmentation of different components of the CT image. The SEM image obtained in this study is an 8-bit grayscale image. The histogram tool in ImageJ software was used to quantitatively count the gray range of each component in the sample, and the average value was obtained. The segmentation thresholds of pyrite and coffinite, feldspar, dolomite, quartz, dolomite, clay minerals, and pores are 242-255, 173-241, 172-146, 105-145, 104-85, 54-85, and 0-53, respectively. Because the gray value in CT images depends on the density of each mineral phase, it is difficult to distinguish between two minerals with similar densities, such as calcite and dolomite. In this study, the CT images were divided into 7 phases, namely pore space, clay minerals, calcite/dolomite, quartz, ankerite, feldspar, pyrite, and coffinite for further analysis and volume performance simulation. Comparing the sample EDS spectrum and map with the segmentation threshold of watershed segmentation, reasonable segmentation thresholds of minerals are obtained. Multi-mineral SEM images after watershed segmentation are shown in Fig. 14b Fig. 14b . Compared with the EDS image ( Fig. 14a Fig. 14a ), it can be seen that good segmentation results can be obtained by using this method. The mineral skeleton of high-resolution CT images is also divided by watershed segmentation. The threshold value of each mineral in the CT image was determined by comparing the gray value and segmentation threshold of the BSE image. The multi-mineral images of CT after watershed segmentation are shown in Fig. 14d Fig. 14d . 3.3.4 Digital core evaluation 3 A variety of morphological algorithms in Avizo Fire are applied to quantitatively characterize the micropore characteristics and mineral content proportion of low-permeability uranium-bearing sandstone. The fast watershed algorithm was used to mark the pores, and each pore was quantitatively analyzed. According to the volume of marked pores, the porosity of the digital core could be obtained [ 25 25 , 60 60 - 62 62 60-62 ]: <math display="block" id="M3"><mrow><mi>ϕ</mi><mtext>=</mtext><mfrac><mrow><mstyle displaystyle="true"><mo>∑</mo><mrow><msub><mi>V</mi><mtext>p</mtext></msub></mrow></mstyle></mrow><mrow><msub><mi>V</mi><mtext>v</mtext></msub></mrow></mfrac><mo>,</mo></mrow></math> where <math id="M4"><mrow><msub><mi>V</mi><mtext>p</mtext></msub></mrow></math> is the single pore volume, pix 3 ; <math id="M5"><mrow><msub><mi>V</mi><mtext>v</mtext></msub></mrow></math> is the total volume, pix 3 . On the basis of the 3D digital core pore model, isolated pores can be removed by using four connected region labeling algorithms, and connected pores can be obtained [ 63 63 ]. The porosity of open pores and close pores can be calculated according to formula (3). It can be seen that the CT images combined with the two resolutions can identify more pores, which is helpful to build an accurate digital core ( Table 4 Table 4 ). 1 1 Sample 10 1 Porosity (%) 1 1 2 1 Low resolution CT 1 1 Total 2 1 High resolution CT 1 1 Total 2 1 Fusion CT 1 1 Total 1 1 HMPI 1 1 1 1 Open pore 1 1 Close pore 1 1 1 1 Open pore 1 1 Close pore 1 1 1 1 Open pore 1 1 Close pore 1 1 1 1 1 1 S1 1 1 6.82 1 1 1.21 1 1 8.03 1 1 8.26 1 1 4.23 1 1 12.49 1 1 9.12 1 1 5.06 1 1 14.72 1 1 16.39 1 1 S2 1 1 5.45 1 1 1.03 1 1 6.48 1 1 5.82 1 1 3.51 1 1 9.33 1 1 6.48 1 1 4.59 1 1 11.07 1 1 12.38 According to the principle of equivalence, the equivalent pore diameter of each marked pore can be obtained [ 25 25 ]: <math display="block" id="M6"><mrow><msub><mi>d</mi><mrow><mtext>eq</mtext></mrow></msub><mo>=</mo><mroot><mrow><mfrac><mrow><mn>3</mn><msub><mi>V</mi><mtext>p</mtext></msub></mrow><mrow><mn>2</mn><mtext>π</mtext></mrow></mfrac></mrow><mn>3</mn></mroot><mtext>   ,</mtext></mrow></math> where d eq is the equivalent pore diameter, pix. The histogram of the equivalent pore diameter distribution is shown in Fig. 15 Fig. 15 . Due to the limitation of CT resolution, pore sizes with a diameter of less than 1μm were not counted, but the distribution of pore size with a diameter of more than 1 μm is in good agreement with the results of the NMR test ( Fig. 9 Fig. 9 ). Similar to the porosity, according to the volume of the marked minerals, the proportion of mineral components of the digital core can be obtained. The proportions of the six mineral phases are listed in Table 5 Table 5 . It can be seen that the digital core mineral volume fraction is consistent with the results of the XRD measurement, and the structural information of different components of the core is reconstructed. 1 1 Sample 1 1 Methods 1 1 Quartz 1 1 Ankerite 1 1 Feldspars 1 1 Calcite 1 1 Dolomite 3 1 Clays 1 1 Metallic minerals 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 Kaolinite 1 1 Illite 1 1 Chlorite 1 1 1 1 S1 1 1 XRD 1 1 46.5 1 1 38.9 1 1 3.9 1 1 0.8 1 1 1.2 1 1 5.4 1 1 0.1 1 1 0.5 1 1 0.7 1 1 1 1 Digital Core 1 1 51.4 1 1 30.17 1 1 3.6 1 1 7.1 1 1 1 1 6.8 1 1 1 1 1 1 1.1 1 1 S2 1 1 XRD 1 1 48.4 1 1 37.5 1 1 5.2 1 1 1.2 1 1 1.3 1 1 3.6 1 1 1.5 1 1 0.1 1 1 1.1 1 1 1 1 Digital Core 1 1 45.9 1 1 36.7 1 1 4.6 1 1 5.4 1 1 1 1 6.8 1 1 1 1 1 1 0.6 4. Discussion 1 The multi-scale and multi-component digital core overcomes the difficulty of balancing the resolution and core size of the single-scale digital core model, and can not only represent the microscopic pore structure of the core, but also describe the macroscopic heterogeneity of the rock sample [ 56 56 ]. The 3D spatial distribution of pores, clay minerals, ankerite cements, pyrite, and coffinite are extracted from the established digital core of low-permeability uranium-bearing sandstone, as shown in Fig. 16 c, d, e, and f Fig. 16 c, d, e, and f . It can be seen from Fig. 16 c Fig. 16 c that the pore distribution in the samples has microscopic heterogeneity characteristics. The primary pores in the area with flaky or banded distribution and relatively rich are mainly related to the residual intergranular pores or large intergranular dissolution pores, and the secondary pores in the area with isolated distribution that are relatively dispersed are mainly related to the intergranular dissolution pores. Combined with the 3D connected pore distribution in Fig. 16 g Fig. 16 g , it can be seen that the intergranular pore forms the main seepage channel in the pore space of the low-permeability uranium-bearing sandstone, and the secondary pore has poor connectivity with other pores [ 64 64 , 65 65 ]. That is to say, intergranular pores mainly provide flow pores and seepage capacity, which is consistent with the conclusion that the contribution of permeability in the samples is mainly from the large throat obtained by the mercury injection experiment. By comparing Fig. 16 j Fig. 16 j , Fig. 16 k Fig. 16 k , and Fig. 16 l Fig. 16 l of the 3D distribution images of the connected pores and pyrite and coffinite, clay minerals, pyrite, coffinite and pyrite, and coffinite and iron dolomite, it can be seen that pyrite and coffinite are isolated from the connected pores and surrounded by a large number of clay minerals and ankerite cements, which increases the difficulty of uranium leaching. Based on the 3D connected digital core model, the pore network model extracted by the maximum sphere algorithm can provide the 3D topological and geometric characteristics of the pore space [ 66 66 , 67 67 ]. The connectivity of the model can be determined by counting the coordination number, which refers to the number of contacts a pore has with its neighbors [ 36 36 ]. The pore network model and its coordination number distribution of low permeability uranium-bearing sandstone under different conditions are shown in Fig. 17 Fig. 17 and Fig. 18 Fig. 18 . The content, diagenetic variation, and occurrence of clay minerals are important factors affecting the porosity and permeability of sandstone. From the 3D distribution of pores and clays in Fig. 16h Fig. 16h , it can be seen that the clays are filled in the primary and secondary pores and pore throats of the samples. When the clay minerals were removed, the coordination number of pores in the pore network model increased, and the connectivity of pores improved ( Fig. 18 Fig. 18 ). Combined with Fig. 16i Fig. 16i , Fig. 17c Fig. 17c , and Fig. 18 Fig. 18 , it is found that the connectivity of the internal pores of the low-permeability uranium-bearing sandstone is strongly affected by a large amount of ankerite cement, and the pore connectivity was significantly improved after the removal of ankerite cements. It can be seen that the low overall permeability of the core is related to the dense ankerite cements filled with a large number of internal pores. In conclusion, it is known that the high relative content of clay minerals and ankerite cement in low-permeability uranium-bearing sandstone is an important factor affecting the permeability and leaching effect of ore and rock. The corresponding lixivium can be designed and prepared to accelerate the dissolution of clay minerals and dolomite, and more specific leaching tests can be carried out to deepen the leaching mechanism. After being verified by a condition test or an extended test, the research results can be applied to in-situ leaching of low-permeability uranium mines. 5. Conclusion 1 To better understand the effect of pore structure and mineral distribution on the seepage characteristics of low-permeability uranium-bearing sandstone, multi-scale and multi-component low permeability uranium-bearing sandstone digital core is reconstructed by fusing different scale CT image information and EDS-SEM images in this paper. The pore distribution of the reconstructed digital core is in good agreement with the NMR results, and the mineral volume fraction of the digital core is in high agreement with the XRD results, which can truly reduce the structural information of different components of the core. The porosity of low-permeability uranium-bearing sandstone was quantitatively and qualitatively analyzed by MICP, NMR, and SEM. The results show that there are four main pore types in the cores of low permeability uranium-bearing sandstone, which include intergranular pores, intragranular dissolution pores, intercrystalline micropores and micro fractures, and even the same mineral in the sample may contain multiple pore types. The pore structure of low-permeability uranium-bearing sandstone is complex and has multi-scale characteristics. Connectivity analysis of the 3D digital core shows that the intergranular pores determine the main seepage channel in the pore space of low-permeability uranium-bearing sandstone, and the secondary pores have poor connectivity with other pores. This is consistent with the conclusion that the contribution of permeability in the low-permeability uranium sandstone sample is mainly from the macro throat obtained by mercury injection experiment. Pyrite and coffinite in low-permeability uranium-bearing sandstone are isolated from the connected pores and surrounded by a large number of clay minerals and ankerite cements, which increases the difficulty of uranium leaching. Moreover, the high content of clays and ankerite cements filled in the primary and secondary pores and pore throats of the samples, which significantly reduces the porosity of the movable fluid in tight sandstone, resulting in the low overall permeability of the core. The multi-scale and multi-mineral digital core proposed in this study provides a basis for characterizing the macroscopic and microscopic pore and pore-throat structure and mineral distribution of low permeability uranium-bearing sandstone. In the numerical simulation, the pore structure can be extracted for seepage simulation, and one or several mineral components can be selectively extracted for physical property simulation. It is an important research direction to simulate the scaled-up physical properties of reservoir rocks using a multi-scale digital core model. Reference S. Zeng , N. Zhang , S.W. Zhang et al ., Fractal characteristics of uranium-bearing sandstone structure and their effects on acid leaching . Energy. Sci. Eng . 7 , 1852 - 1866 , ( 2019 ). https://doi.org/10.1002/ese3.396 https://doi.org/10.1002/ese3.396 Fractal characteristics of uranium-bearing sandstone structure and their effects on acid leaching F.X. Chen , F.J Nie , C.Y. Zhang et al ., Sedimentary facies of the Lower Jurassic Sangonghe Formation and its relationship with uranium mineralization in Mengqiguer area of Yili Basin, Xinjiang . J. Palaeog . 18 , 833 - 842 , ( 2016 ). https://doi.org/10.7605/gdlxb.2016.05.063 https://doi.org/10.7605/gdlxb.2016.05.063 (in Chinese) Sedimentary facies of the Lower Jurassic Sangonghe Formation and its relationship with uranium mineralization in Mengqiguer area of Yili Basin, Xinjiang B.R. Song , H. Sun , S.L. Yang et al ., Characteristics and uranium mineralization of ore bearing rock seriesin Qianjiadian sandstone -type uranium deposits, Songliao Basin . J. Palaeogeog.-Chinese . 22 , 309 - 320 , ( 2020 ). https://doi.org/10.2113/econgeo.106.5.835  https://doi.org/10.2113/econgeo.106.5.835  (in Chinese) Characteristics and uranium mineralization of ore bearing rock seriesin Qianjiadian sandstone -type uranium deposits, Songliao Basin H. Song , S.H. Ni , M.C. Hou et al ., The characteristic of clay minerals in sandstone-type uranium deposit in the Yili Basin, NW China and its relationship with uranium mineralization . ACTA. GEOL. SIN . 90 , 3352 - 3366 , ( 2016 ). https://doi.org/10.3969/j.issn.0001-5717.2016.12.005 https://doi.org/10.3969/j.issn.0001-5717.2016.12.005 (in Chinese) The characteristic of clay minerals in sandstone-type uranium deposit in the Yili Basin, NW China and its relationship with uranium mineralization Z.M. Yang , X.L. Zhao , C.S. Xiong et al ., Research progress on microstructure characterization of pore throat for tight oil reservoirs . Sci. Tech. Rev . 37 , 89 - 98 , ( 2019 ) https://doi.org/10.3981/j.issn.1000-7857.2019.05.012 https://doi.org/10.3981/j.issn.1000-7857.2019.05.012 (in Chinese) Research progress on microstructure characterization of pore throat for tight oil reservoirs X.H. Shao , X.Q. Pang , F.J. Jiang et al ., Reservoir characterization of tight sandstones using nuclear magnetic resonance and incremental pressure mercury injection experiments: implication for tight sand gas reservoir quality . Energ. Fuel . 31 , 10420 - 10431 , ( 2017 ), https://doi.org/10.1021/acs.energyfuels.7b01184 https://doi.org/10.1021/acs.energyfuels.7b01184 Reservoir characterization of tight sandstones using nuclear magnetic resonance and incremental pressure mercury injection experiments: implication for tight sand gas reservoir quality H.W. Zhou , J.C. Zhong , W.G. Ren et al ., Characterization of pore- fracture networks and their evolution at various measurement scales in coal samples using X-ray μCT and a fractal method . Int. J. Coal. Geol . 189 , 35 - 49 , ( 2018 ), https://doi.org/10.1016/j.coal.2018.02.007 https://doi.org/10.1016/j.coal.2018.02.007 Characterization of pore- fracture networks and their evolution at various measurement scales in coal samples using X-ray μCT and a fractal method H.X. Huang , W. Sun , W.M. Ji et al ., Effects of pore-throat structure on gas permeability in the tight sandstone reservoirs of the Upper Triassic Yanchang formation in the Western Ordos Basin, China . J Petrol. Sci. Eng . 162 , 602 - 616 , ( 2017 ). https://doi.org/10.1016/j.petrol.2017.10.076 https://doi.org/10.1016/j.petrol.2017.10.076 Effects of pore-throat structure on gas permeability in the tight sandstone reservoirs of the Upper Triassic Yanchang formation in the Western Ordos Basin, China J. Lai , G.W. Wang , Z.Y. Fa et al ., Insight into the pore structure of tight sandstones using NMR and HPMI measurements . Energ. Fuel . 30 , 10200 - 10214 , ( 2016 ) https://doi.org/10.1021/acs.energyfuels.6b01982 https://doi.org/10.1021/acs.energyfuels.6b01982 Insight into the pore structure of tight sandstones using NMR and HPMI measurements M. Schmitt , C.P. Fernandes , F.G. Wolf et al ., Characterization of Brazilian tight gas sandstones relating permeability and Angstrom-to micron-scale pore structures . J. Nat. Gas. Sci. Eng . 27 , 785 - 807 , ( 2015 ), https://doi.org/10.1016/j.jngse.2015.09.027 https://doi.org/10.1016/j.jngse.2015.09.027 Characterization of Brazilian tight gas sandstones relating permeability and Angstrom-to micron-scale pore structures J. Klus , P. Mikysek , D. Prochazka et al ., Multivariate approach to the chemical mapping of uranium in sandstone-hosted uranium ores analyzed using double pulse Laser-Induced Breakdown Spectroscopy . Spectrochim. Acta. B . 123 , 143 - 149 , ( 2016 ), https://doi.org/10.1016/j.sab.2016.08.014 https://doi.org/10.1016/j.sab.2016.08.014 Multivariate approach to the chemical mapping of uranium in sandstone-hosted uranium ores analyzed using double pulse Laser-Induced Breakdown Spectroscopy Z.L. Cheng , Z.F. Ning , H.W. Zhao et al ., A comprehensive characterization of North China tight sandstone using micro-CT, SEM imaging, and mercury intrusion . Arab. J. Geosci . 12 , 407 , ( 2019 ). https://doi.org/10.1007/s12517-019-4568-9 https://doi.org/10.1007/s12517-019-4568-9 A comprehensive characterization of North China tight sandstone using micro-CT, SEM imaging, and mercury intrusion A. Kazak , S. Chugunov , A. Chashkov . Integration of large-area scanning-electron-microscopy imaging and automated mineralogy/petrography data for selection of nanoscale pore-space characterization sites . Spe. Reserv. Eval. Eng . 21 , 821 - 836 , ( 2018 ). https://doi.org/10.2118/191369-Pa https://doi.org/10.2118/191369-Pa Integration of large-area scanning-electron-microscopy imaging and automated mineralogy/petrography data for selection of nanoscale pore-space characterization sites W. De. Boever , H. Derluyn , D. Van Loo et al ., Data-fusion of high resolution X-ray CT, SEM and EDS for 3D and pseudo-3D chemical and structural characterization of sandstone . Micron . 74 , 15 - 21 , ( 2015 ). https://doi.org/10.1016/j.micron.2015.04.003 https://doi.org/10.1016/j.micron.2015.04.003 Data-fusion of high resolution X-ray CT, SEM and EDS for 3D and pseudo-3D chemical and structural characterization of sandstone M. Shabaninejad , J. Middleton , S. Latham et al ., Pore-scale analysis of residual oil in a reservoir sandstone and its dependence on water flood salinity, oil composition and local mineralogy . Energ. Fuel . 31 , 13221 - 13232 , ( 2017 ), https://doi.org/10.1021/acs.energyfuels.7b01978 https://doi.org/10.1021/acs.energyfuels.7b01978 Pore-scale analysis of residual oil in a reservoir sandstone and its dependence on water flood salinity, oil composition and local mineralogy K.L. Xi , Y.C. Cao , K.Y. Liu et al ., Authigenic minerals related to wettability and their impacts on oil accumulation in tight sandstone reservoirs: An example from the Lower Cretaceous Quantou Formation in the southern Songliao Basin, China . J. Asian. Earth. Sci . 178 , 173 - 192 , ( 2018 ),   https://doi.org/10.1016/j.jas.2010.03.005 https://doi.org/10.1016/j.jas.2010.03.005 Authigenic minerals related to wettability and their impacts on oil accumulation in tight sandstone reservoirs: An example from the Lower Cretaceous Quantou Formation in the southern Songliao Basin, China L. Ganzer , S. Schatzmann , C.D. Sattler et al ., Evaluation of digital rock methodology to complement rock laboratory experiments . Oil. Gas-Euro. Mag . 39 , 43 - 47 , ( 2013 ). https://www.researchgate.net/publication/287946845 https://www.researchgate.net/publication/287946845 Evaluation of digital rock methodology to complement rock laboratory experiments L. Wang , Y.M. He , X. Peng et al ., Pore structure characteristics of an ultradeep carbonate gas reservoir and their effffects on gas storage and percolation capacities in the Deng IV member, Gaoshiti-Moxi Area, Sichuan Basin, SW China . Mar. Petrol. Geol . 111 , 44 - 65 , ( 2020 ). https://doi.org/10.1016/j.marpetgeo.2019.08.012 https://doi.org/10.1016/j.marpetgeo.2019.08.012 Pore structure characteristics of an ultradeep carbonate gas reservoir and their effffects on gas storage and percolation capacities in the Deng IV member, Gaoshiti-Moxi Area, Sichuan Basin, SW China X.J. Liu , H.L. Zhu , L.X. Liang . Digital rock physics of sandstone based on micro-CT technology . Chinese. Journal. Of. Geophysics . 57 , 1133 - 1140 , ( 2014 ). https://doi.org/10.6038/cjg20140411 https://doi.org/10.6038/cjg20140411 Digital rock physics of sandstone based on micro-CT technology B. Bai , R. Zhu , S. Wu et al ., Multi-scale method of Nano (Micro)-CT study on microscopic pore structure of tight sandstone of Yanchang Formation, Ordos Basin . Petrol. Explor. Dev+ . 40 , 354 - 358 , ( 2013 ). https://doi.org/10.1016/s1876-3804(13)60042-7 https://doi.org/10.1016/s1876-3804(13)60042-7 Multi-scale method of Nano (Micro)-CT study on microscopic pore structure of tight sandstone of Yanchang Formation, Ordos Basin S.W. Zhou , G. Yan , H.Q. Xue et al ., 2D and 3D nanopore characterization of gas shale in Longmaxi formation based on FIB-SEM . Mar. Petrol. Geol . 73 , 174 - 180 , ( 2016 ). https://doi.org/10.1016/j.marpetgeo.2016.02.033 https://doi.org/10.1016/j.marpetgeo.2016.02.033 2D and 3D nanopore characterization of gas shale in Longmaxi formation based on FIB-SEM Z. Zhao , X.P. Zhou , An integrated method for 3D reconstruction model of porous geomaterials through 2D CT images . Comput. Geosci . ( 2018 ). https://doi.org/10.1016/j.cageo.2018.11.012 https://doi.org/10.1016/j.cageo.2018.11.012 . An integrated method for 3D reconstruction model of porous geomaterials through 2D CT images C.Y. Lin , Q.Y. Wu , L.H. Ren et al ., Review of digital core modeling methods . Prog. Geophysics . 33 , 679 - 689 , ( 2018 ). https://doi.org/10.6038/pg2018BB0335 https://doi.org/10.6038/pg2018BB0335 (in Chinese) Review of digital core modeling methods T.F. Faisal , A. Islam , M.S. Jouini et al ., Numerical prediction of carbonate elastic properties based on multi-scale imaging . Geomech. Energy. Envir . 20 , 1 - 15 , ( 2019 ). https://doi.org/10.1016/j.gete.2019.100125 https://doi.org/10.1016/j.gete.2019.100125 Numerical prediction of carbonate elastic properties based on multi-scale imaging Z. Zhao , X.P. Zhou . Digital microstructure insights to phase evolution and thermal flow properties of hydrates by X-ray computed tomography . Sci. China. Technol. Sc . 1 - 16 . ( 2020 ). https://doi.org/10.1007/s11431-020-1579-x https://doi.org/10.1007/s11431-020-1579-x Digital microstructure insights to phase evolution and thermal flow properties of hydrates by X-ray computed tomography S.Y. An , J. Yao , Y.F. Yang et al ., Influence of pore structure parameters on flow characteristics based on a digital rock and the pore network model . J Nat. Gas. Sci. Eng . 31 , 156 - 163 , ( 2016 ). https://doi.org/10.1016/j.jngse.2016.03.009 https://doi.org/10.1016/j.jngse.2016.03.009 Influence of pore structure parameters on flow characteristics based on a digital rock and the pore network model L. Zhang , W.L. Jing , Y.F. Yang et al ., The Investigation of permeability calculation using digital core simulation technology . Energies . 12 , 3273 , ( 2019 ). https://doi.org/10.3390/en12173273 https://doi.org/10.3390/en12173273 The Investigation of permeability calculation using digital core simulation technology J.J. Li , Y. Liu , Y.J. Gao et al ., Pore-scale study of the pressure-sensitive effect of sandstone and its influence on multiphase flows . Petrol. Sci . 16 , 382 - 395 , ( 2019 ). https://doi.org/10.1007/s12182-018-0266-6 https://doi.org/10.1007/s12182-018-0266-6 Pore-scale study of the pressure-sensitive effect of sandstone and its influence on multiphase flows J. Lei , B.Z. Pan , L.H. Zhang . Advance of microscopic flow simulation based on digital cores and pore network . Prog. Geophysics . 33 , 653 - 660 , ( 2018 ). https://doi.org/10.6038/pg2018BB0108 https://doi.org/10.6038/pg2018BB0108 (in Chinese) Advance of microscopic flow simulation based on digital cores and pore network M.S. Jouini , S. Vega & A. Al-Ratrout . Numerical estimation of carbonate rock properties using multiscale images . Geophys. Prospect . 63 , 405 - 421 , ( 2015 ). https://doi.org/10.1111/1365-2478.12156 https://doi.org/10.1111/1365-2478.12156 Numerical estimation of carbonate rock properties using multiscale images M.J. Blunt , B. Bijeljic , H. Dong et al ., Pore-scale imaging and modelling . Adv. Water. Resour . 51 , 197 - 216 , ( 2013 ). https://doi.org/10.1016/j.advwatres.2012.03.003 https://doi.org/10.1016/j.advwatres.2012.03.003 Pore-scale imaging and modelling Z. Zhao , X.P. Zhou . Establishment of numerical cracking constitutive models using 3D reconstruction and X-ray CT images of geomaterials , Int. J. Mech. Sci . 183 , 105814 , ( 2020 ). https://doi.org/10.1016/j.ijmecsci.2020.105814 https://doi.org/10.1016/j.ijmecsci.2020.105814 Establishment of numerical cracking constitutive models using 3D reconstruction and X-ray CT images of geomaterials J.E. Elkhoury , R. Shankar , T.S. Ramakrishnan , Resolution and limitations of x-ray micro-ct with applications to sandstones and limestones . Transport. Porous. Med . 129 , 413 - 425 , ( 2019 ). https://doi.org/10.1007/s11242-019-01275-1 https://doi.org/10.1007/s11242-019-01275-1 Resolution and limitations of x-ray micro-ct with applications to sandstones and limestones R. S. Devarapalli , A. Islam , T. F. Faisal et al ., Micro-CT and FIB-SEM imaging and pore structure characterization of dolomite rock at multiple scales . Arab. J. Geosci . 10 , 360 - 372 , ( 2017 ). https://doi.org/10.1007/s12517-017-3120-z https://doi.org/10.1007/s12517-017-3120-z Micro-CT and FIB-SEM imaging and pore structure characterization of dolomite rock at multiple scales S. Latham , T. Varslot , A. Sheppard , Image registration: Enhancing and calibrating X-ray micro-CT imaging . Imaging Society of Core Analysts Paper . 2008 , 1 - 12 , ( 2008 ). https://www.researchgate.net/publication/254935632 https://www.researchgate.net/publication/254935632 Image registration: Enhancing and calibrating X-ray micro-CT imaging Ø. Markussen , H. Dypvik , E. Hammer et al ., 3D characterization of porosity and authigenic cementation in Triassic conglomerates/arenites in the Edvard Grieg field using 3D micro-CT imaging . Mar. Petro. Geol . 99 , 265 - 281 , ( 2019 ). https://doi.org/10.1016/j.marpetgeo.2018.10.015 https://doi.org/10.1016/j.marpetgeo.2018.10.015 3D characterization of porosity and authigenic cementation in Triassic conglomerates/arenites in the Edvard Grieg field using 3D micro-CT imaging G. Desbois , J.L. Urai , P.A. Kukla et al ., High-resolution 3D fabric and porosity model in a tight gas sandstone reservoir: A new approach to investigate microstructures from mm- to nm-scale combining argon beam cross-sectioning and SEM imaging . J. Petrol. Sci. Eng . 78 , 243 - 257 , ( 2011 ). https://doi.org/10.1016/j.petrol.2011.06.004 https://doi.org/10.1016/j.petrol.2011.06.004 High-resolution 3D fabric and porosity model in a tight gas sandstone reservoir: A new approach to investigate microstructures from mm- to nm-scale combining argon beam cross-sectioning and SEM imaging H. Izadi , M. Baniassadi , F. Hormozzade et al ., Effect of 2d image resolution on 3d stochastic reconstruction and developing petrophysical trend . Transport. Porous. Med . 125 , 41 - 58 , ( 2018 ). https://doi.org/10.1007/s11242-018-0997-2 https://doi.org/10.1007/s11242-018-0997-2 Effect of 2d image resolution on 3d stochastic reconstruction and developing petrophysical trend W. De. Boever , H. Derluyn , D. Van Loo et al ., Data-fusion of high resolution X-ray CT, SEM and EDS for 3D and pseudo-3D chemical and structural characterization of sandstone . Micron . 74 , 15 - 21 , ( 2015 ). https://doi.org/10.1016/j.micron.2015.04.003 https://doi.org/10.1016/j.micron.2015.04.003 Data-fusion of high resolution X-ray CT, SEM and EDS for 3D and pseudo-3D chemical and structural characterization of sandstone X.F. Liu , J.F. Wang , L. Ge et al ., Pore-scale characterization of tight sandstone in Yanchang Formation Ordos Basin China using micro-CT and SEM imaging from nm- to cm-scale . Fuel . 209 , 254 - 264 , ( 2017 ). https://doi.org/10.1016/j.fuel.2017.07.068 https://doi.org/10.1016/j.fuel.2017.07.068 Pore-scale characterization of tight sandstone in Yanchang Formation Ordos Basin China using micro-CT and SEM imaging from nm- to cm-scale W. Lin , X.Z. Li , Z.M. Yang et al ., Construction of dual pore 3-d digital cores with a hybrid method combined with physical experiment method and numerical reconstruction method . Transport. Porous. Med . 120 , 227 - 238 , ( 2017 ). https://doi.org/10.1007/s11242-017-0917-x https://doi.org/10.1007/s11242-017-0917-x Construction of dual pore 3-d digital cores with a hybrid method combined with physical experiment method and numerical reconstruction method C. Tian , L.C. Jia , S. Li et al ., Petrology characteristics of uranium-bearing sandstone in the southern Ordos basin . U. Geol . 2 , 71 - 76 , ( 2007 ). https://doi.org/CNKI:SUN:YKDZ.0.2007-02-001 https://doi.org/CNKI:SUN:YKDZ.0.2007-02-001 (in Chinese) Petrology characteristics of uranium-bearing sandstone in the southern Ordos basin S. Zeng , K.X. Tan , X.M. Du et al ., The effects of structural characteristics of low seepage sandstone uranium ore rock on permeability . Adv. Mater . 524 - 527 , 371 - 374 , ( 2012 ). https://doi.org/10.4028/www.scientific.net/AMR.524-527.371 https://doi.org/10.4028/www.scientific.net/AMR.524-527.371 The effects of structural characteristics of low seepage sandstone uranium ore rock on permeability H. Li , H. Zhou . Progress and prospect of research on solute transport during in-situ leaching of uranium . Chinese. J. Rare. Metals . 43 , 319 - 330 . https://doi.org/10.13373/j.cnki.cjrm.XY18120015 https://doi.org/10.13373/j.cnki.cjrm.XY18120015 Progress and prospect of research on solute transport during in-situ leaching of uranium C. Tang , X.Z. Sima , Q. Zhu et al ., Carbon and oxygen isotopic composition and uranium mineralization significance of calcite of Zhiluo Formation uranium-bearing sandstones in Dongsheng area . Geol. J. China. Universities . 22 , 698 - 706 , ( 2016 ). https://doi.org/10.16108/j.issn1006-7493.2016055 https://doi.org/10.16108/j.issn1006-7493.2016055 (in Chinese) Carbon and oxygen isotopic composition and uranium mineralization significance of calcite of Zhiluo Formation uranium-bearing sandstones in Dongsheng area G.W. Huang , F.C. Yu , J.Y. Pan et al ., Geochemical characteristics of trace elements and sulfur isotopes of pyrite from Mengqiguer uranium deposit, Xinjiang: Implication for uranium mineralization . Geo. China . 1 - 16 , ( 2019 ), https://kns.cnki.net/kcms/detail/11.1167.P.20191230.1721.009.html https://kns.cnki.net/kcms/detail/11.1167.P.20191230.1721.009.html (in Chinese) Geochemical characteristics of trace elements and sulfur isotopes of pyrite from Mengqiguer uranium deposit, Xinjiang: Implication for uranium mineralization X.X. Kong , D.S. Xiao , S. Jiang et al ., Application of the combination of high-pressure mercury injection and nuclear magnetic resonance to the classification and evaluation of tight sandstone reservoirs: A case study of the Linxing Block in the Ordos Basin . Nat. Gas. Ind . 40 , 38 - 47 , ( 2020 ) https://doi.org/10.3787/j.issn.1000-0976.2020.03.005 https://doi.org/10.3787/j.issn.1000-0976.2020.03.005 (in Chinese) Application of the combination of high-pressure mercury injection and nuclear magnetic resonance to the classification and evaluation of tight sandstone reservoirs: A case study of the Linxing Block in the Ordos Basin G. Wang , X.J. Qin , J.N. Shen et al ., Quantitative analysis of microscopic structure and gas seepage characteristics of low-rank coal based on CT three-dimensional reconstruction of CT images and fractal theory . Fuel . 256 , 1 - 11 , ( 2019 ). https://doi.org/10.1016/.fuel.2019.115900 https://doi.org/10.1016/.fuel.2019.115900 Quantitative analysis of microscopic structure and gas seepage characteristics of low-rank coal based on CT three-dimensional reconstruction of CT images and fractal theory H.W. Zhou , J.C. Zhong , W.G. Ren et al ., Characterization of pore fracture networks and their evolution at various measurement scales in coal samples using X-ray μCT and a fractal method . Int. J. Coal. Geol . 189 , 35 - 49 , ( 2018 ). https://doi.org/10.1016/j.coal.2018.02.007 https://doi.org/10.1016/j.coal.2018.02.007 Characterization of pore fracture networks and their evolution at various measurement scales in coal samples using X-ray μCT and a fractal method X.X. Kong , D.S. Xiao , S. Jiang et al ., Application of the combination of high-pressure mercury injection and nuclear magnetic resonance to the classification and evaluation of tight sandstone reservoirs: A case study of the Linxing Block in the Ordos Basin . Nat. Gas. Ind . 40 , 38 - 47 , ( 2020 ) https://doi.org/10.3787/j.issn.1000-0976.2020.03.005 https://doi.org/10.3787/j.issn.1000-0976.2020.03.005 (in Chinese) Application of the combination of high-pressure mercury injection and nuclear magnetic resonance to the classification and evaluation of tight sandstone reservoirs: A case study of the Linxing Block in the Ordos Basin Lai J , Wang G , Wang Z , et al . A review on pore structure characterization in tight sandstones . Earth-Science Reviews , 2018 , 177 : 436 - 457 . https://doi.org/10.1016/j.earscirev.2017.12.003 https://doi.org/10.1016/j.earscirev.2017.12.003 A review on pore structure characterization in tight sandstones P.H. Nelson . Pore-throat sizes in sandstones, tight sandstones, and shales . AAPG Bulletin , 2009 , 93 ( 3 ): 329 - 340 . https://doi.org/10.1306/10240808059 https://doi.org/10.1306/10240808059 . Pore-throat sizes in sandstones, tight sandstones, and shales Z.L. Cheng , Z.F. Ning , H.W. Zhao et al ., A comprehensive characterization of North China tight sandstone using micro-CT, SEM imaging, and mercury intrusion . Arab. J. Geosci . 12 , 407 , ( 2019 ). https://doi.org/10.1007/s12517-019-4568-9 https://doi.org/10.1007/s12517-019-4568-9 A comprehensive characterization of North China tight sandstone using micro-CT, SEM imaging, and mercury intrusion Avizo , 2019 . Avizo Software Users Guide . Available at: http://www.fei.com/software/avizo-user-guide http://www.fei.com/software/avizo-user-guide . Avizo Software Users Guide F.L. Bookstein . Size and shape spaces for landmark data in two dimensions . Stat. Sci . 1 , 181 - 222 , ( 1986 ). https://doi.org/10.1214/ss/1177013696 https://doi.org/10.1214/ss/1177013696 Size and shape spaces for landmark data in two dimensions L.C. Cui , J.M. Sun , W.C. Yan et al ., Construction of multi-scale and-component digital cores based on fusion of different resolution core images . Journal of Jilin University: Earth Science Edition 47 , 1904 - 1912 , ( 2017 ). https://doi.org/CNKI:SUN:CCDZ.0.2017-06-027 https://doi.org/CNKI:SUN:CCDZ.0.2017-06-027 (in Chinese) Construction of multi-scale and-component digital cores based on fusion of different resolution core images S.J. Redmond and C. Heneghan . A method for initialising the K-means clustering algorithm using kd-trees . Pattern. Recogn. Lett . 28 , 965 - 973 , ( 2007 ) https://doi.org/10.1016/j.patrec.2007.01.001 https://doi.org/10.1016/j.patrec.2007.01.001 A method for initialising the K-means clustering algorithm using kd-trees Y.W. Gao , W. Wu , Q. Zhang et al ., CT image segmentation using multi-Otsu thresholding algorithm and scanning electron microscopy . J. Yangtze. University. (Natural Science Edition) . 16 , 17-24+4-5 , ( 2019 ). https://doi.org/10.16772/j.cnki.1673-1409.2019.04.003 https://doi.org/10.16772/j.cnki.1673-1409.2019.04.003 (in Chinese) CT image segmentation using multi-Otsu thresholding algorithm and scanning electron microscopy C. Chen , S.H. Lu , J.Q. Li et al ., Digital core modeling construction of different lithofacies shale: a case study of Dongying depression . Geosci . 31 , 1069 - 1078 , ( 2017 ). https://doi.org/10.3969/j.issn.1000-8527.2017.05.016 https://doi.org/10.3969/j.issn.1000-8527.2017.05.016 (in Chinese) Digital core modeling construction of different lithofacies shale: a case study of Dongying depression X.P. Zhou , Y.X. Zhang and Q.L. Ha . Real-time computerized tomography (CT) experiments on limestone damage evolution during unloading‎ . Theor. Appl. Fract. Mech . 50 , 49 - 56 , ( 2008 ). https://doi.org/10.1016/j.tafmec.2008.04.005 https://doi.org/10.1016/j.tafmec.2008.04.005 Real-time computerized tomography (CT) experiments on limestone damage evolution during unloading‎ Z. Zhao , X.P. Zhou . Pore-scale effect on the hydrate variation and flow behaviors in microstructures using X-ray CT imaging . J. Hydrol . 584 , 124678 , ( 2020 ). https://doi.org/10.1016/j.jhydrol.2020.124678 https://doi.org/10.1016/j.jhydrol.2020.124678 Pore-scale effect on the hydrate variation and flow behaviors in microstructures using X-ray CT imaging Z. Zhao , X.P. Zhou . Digital measurement of 2D and 3D cracks in sandstones through improved pseudo color image enhancement and 3D reconstruction method . Int. J. Numl. Met . 43 , 2565 - 2584 , ( 2019 ). https://doi.org/10.1002/nag.2993 https://doi.org/10.1002/nag.2993 Digital measurement of 2D and 3D cracks in sandstones through improved pseudo color image enhancement and 3D reconstruction method J.J. Wang , Q.J. Kang , Y.Z. Wang et al ., Simulation of gas flow in micro- porous media with the regularized Lattice Boltzmann method . Fuel . 205 , 232 - 246 , ( 2017 ). https://doi.org/10.1016/j.fuel.2017.05.080 https://doi.org/10.1016/j.fuel.2017.05.080 Simulation of gas flow in micro- porous media with the regularized Lattice Boltzmann method ‏Z. Zhao , X.P. Zhou and Q.H. Qian . Fracture characterization and permeability prediction by pore scale variables extracted from X-ray CT images of porous geomaterials , Sci. China. Technol. Sc . 63 , 755 - 767 , ( 2020 ). https://doi.org/10.1007/s11431-019-1449-4 https://doi.org/10.1007/s11431-019-1449-4 Fracture characterization and permeability prediction by pore scale variables extracted from X-ray CT images of porous geomaterials X.P. Zhou , Z. Zhao and Z. Li , Cracking behaviors and hydraulic properties evaluation based on fractural microstructure models in geomaterials . Int. J. Rock. Mech. Min. Sci . 130 , 104304 . ( 2020 ). https://doi.org/10.1016/j.ijrmms.2020.104304 https://doi.org/10.1016/j.ijrmms.2020.104304 Cracking behaviors and hydraulic properties evaluation based on fractural microstructure models in geomaterials H.M. Dong , J.M. Sun , Z.Z. Lin et al ., Quantitative characterization and characteristics analysis of microscopic pore structure in natural gas hydrate based on CT scanning . J. China. Un. Petrol. Nat. Sci . 42 , 40 - 49 , ( 2018 ). https://doi.org/10.3969/j.issn.1673-5005.2018.06.005 https://doi.org/10.3969/j.issn.1673-5005.2018.06.005 (in Chinese) Quantitative characterization and characteristics analysis of microscopic pore structure in natural gas hydrate based on CT scanning D. Zheng , Z. Reza . Pore-network extraction algorithm for shale accounting for geometry-effect . J. Petrol. Sci. Eng . 176 , 74 - 84 , ( 2019 ). https://doi.org/10.1016/j.petrol.2019.01.046 https://doi.org/10.1016/j.petrol.2019.01.046 Pore-network extraction algorithm for shale accounting for geometry-effect 10.1007/s41365-020-00810-w This work was supported by the National Natural Science Foundation of China (No.11775107) and the Key Projects of Education Department of Hunan Province of China (No. 16A184). 2020-05-13 2020-09-03 2020-09-09 2020-10-15 2020-10-01 2021-04-10 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2020 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-020-00810-w http://dx.doi.org/10.1007/s41365-020-00810-w . 2020 Nuclear Science and Techniques 10 31