Reconstruction of three-dimensional porous media from a single two-dimensional image using three-step sampling

Three-dimensional reconstruction of granular porous media based on deep generative models

Reconstruction of microstructure in granular porous media, which can be viewed as granular assemblies, is crucial for studying their characteristics and physical properties in various fields concerned with the behavior of such media, including petroleum geology and computational materials science. In spite of the fact that many existing studies have investigated grain reconstruction, most of them treat grains as simplified individuals for discrete reconstruction, which cannot replicate the complex geometrical shapes and natural interactions between grains. In this work, a hybrid generative model based on a deep-learning algorithm is proposed for high-quality three-dimensional (3D) microstructure reconstruction of granular porous media from a single two-dimensional (2D) slice image. The method extracts 2D prior information from the given image and generates the grain set as a whole. Both a self-attention module and effective pattern loss are introduced in a bid to enhance the reconstruction ability of the model. Samples with grains of varied geometrical shapes are utilized for the validation of our method, and experimental results demonstrate that our proposed approach can accurately reproduce the complex morphology and spatial distribution of grains without any artificiality. Furthermore, once the model training is complete, rapid end-to-end generation of diverse 3D realizations from a single 2D image can be achieved.

Improved recurrent generative model for reconstructing large-size porous media from two-dimensional images

Modeling the three-dimensional (3D) structure from a given 2D image is of great importance for analyzing and studying the physical properties of porous media. As an intractable inverse problem, many methods have been developed to address this fundamental problems over the past decades. Among many methods, the deep learning-(DL) based methods show great advantages in terms of accuracy, diversity, and efficiency. Usually, the 3D reconstruction from the 2D slice with a larger field-of-view is more conducive to simulate and analyze the physical properties of porous media accurately. However, due to the limitation of reconstruction ability, the reconstruction size of most widely used generative adversarial network-based model is constrained to 64^{3} or 128^{3}. Recently, a 3D porous media recurrent neural network based method (namely, 3D-PMRNN) (namely 3D-PMRNN) has been proposed to improve the reconstruction ability, and thus the reconstruction size is expanded to 256^{3}. Nevertheless, in order to train these models, the existed DL-based methods need to down-sample the original computed tomography (CT) image first so that the convolutional kernel can capture the morphological features of training images. Thus, the detailed information of the original CT image will be lost. Besides, the 3D reconstruction from a optical thin section is not available because of the large size of the cutting slice. In this paper, we proposed an improved recurrent generative model to further enhance the reconstruction ability (512^{3}). Benefiting from the RNN-based architecture, the proposed model requires only one 3D training sample at least and generates the 3D structures layer by layer. There are three more improvements: First, a hybrid receptive field for the kernel of convolutional neural network is adopted. Second, an attention-based module is merged into the proposed model. Finally, a useful section loss is proposed to enhance the continuity along the Z direction. Three experiments are carried out to verify the effectiveness of the proposed model. Experimental results indicate the good reconstruction ability of proposed model in terms of accuracy, diversity, and generalization. And the effectiveness of section loss is also proved from the perspective of visual inspection and statistical comparison.

Three-dimensional reconstruction of porous media using super-dimension-based adjacent block-matching algorithm

As porous media play an essential role in a variety of industrial applications, it is essential to understand their physical properties. Nowadays, the super-dimensional (SD) reconstruction algorithm is used to stochastically reconstruct a three-dimensional (3D) structure of porous media from a given two-dimensional image. This algorithm exhibits superiority in accuracy compared with classical algorithms because it learns information from the real 3D structure. However, owing to the short development time of the SD algorithm, it also has some limitations, such as inexact porosity characterization, long run time, blocking artifacts, and suboptimal accuracy that may be improved. To mitigate these limitations, this study presents the design of a special template that contains two parts of data (i.e., adjacent blocks and a central block); the proposed method matches adjacent blocks during reconstruction and assigns the matched central block to the area to be reconstructed. Furthermore, we design two important mechanisms during reconstruction: one for block matching and the other for porosity control. To verify the effectiveness of the proposed method compared with an existing SD method, both methods were tested on silica particle material and three homogeneous sandstones with different porosities; meanwhile, we compared the proposed method with a multipoint statistics method and a simulated annealing method. The reconstructed results were then compared with the target both visually and quantitatively. The experimental results indicate that the proposed method can overcome the aforementioned limitations and further improve the accuracy of existing methods. This method achieved 4-6 speedup factor compared with the traditional SD method.

Reconstruction of 3D greyscale image for reservoir rock from a single image based on pattern dictionary

Most methods that model 3D porous media from 2D images are based on binary images. In this paper, we propose a method for reconstructing 3D greyscale isotropic porous media images from a single image. Our proposed method incorporates a fast-sampling procedure to control the continuity and variability between adjoining reconstruction layers, a new similarity calculation method to obtain the most similar patterns from a pattern dictionary, and a central area simulation procedure to solve the block effect problem. The reconstruction results from application of our proposed method to a real reservoir 3D model obtained via computed tomography (CT) and a comparison with the original CT structure demonstrate that our proposed method can reproduce properties such as autocorrelation function, linear function, shape distribution, average shape factor, average pore radius size, average throat radius size, average pore volume, permeability and grey histogram. Further, the comparison results indicate that the statistical characteristics of the reconstructions match the training image and the CT model perfectly.

Reconstructing 3D digital model without distortion for poorly conductive porous rock by nanoprobe-assisted FIB-SEM tomography

Oil and natural gas prospecting requires precise pore characterisation of insulating rock samples, which involves severe charging problems in the state-of-art FIB-SEM tomography, such as overexposure, drift and distortion. For weak cemented samples with very poor conductivity, the conventional ways such as decreasing accelerating voltage or current as well as coating a thin layer of carbon or gold fail to eliminate all the detrimental effect, leading to image distortion in the form of lateral shift and longitudinal stretching. A new nanoprobe-assisted method is explored in FIB-SEM tomography to address this problem and improve image quality. To be specific, a metallic nanoprobe is induced and attached on the sample surface to create an express path for the export of excess electrons near the region of interest, which effectively removes distortion and drift when imaging. Two adjacent areas were characterised and reconstructed into 3D digital models by FIB-SEM tomography with nanoprobe-assisted method applied to one region only. The lateral shift creates zigzag feature for distorted region and the longitudinal stretching of undistorted object can reach 14%. Average pore size of distorted region is larger than that of the undistorted region, however considering the longitudinal stretching, the average pore size of distorted region can be corrected to the same level as the undistorted region. The systematic error caused by distortion for poorly conductive porous rock is hazardous for digital rock physics analysis. Therefore, the nanoprobe-assisted FIB-SEM tomography should be regarded as a one of the optional and feasible procedures in case decreasing accelerating voltage or current as well as coating a thin layer of conductive material does not work.

A Two-Stage Reconstruction of Microstructures with Arbitrarily Shaped Inclusions

The main goal of our research is to develop an effective method with a wide range of applications for the statistical reconstruction of heterogeneous microstructures with compact inclusions of any shape, such as highly irregular grains. The devised approach uses multi-scale extended entropic descriptors (ED) that quantify the degree of spatial non-uniformity of configurations of finite-sized objects. This technique is an innovative development of previously elaborated entropy methods for statistical reconstruction. Here, we discuss the two-dimensional case, but this method can be generalized into three dimensions. At the first stage, the developed procedure creates a set of black synthetic clusters that serve as surrogate inclusions. The clusters have the same individual areas and interfaces as their target counterparts, but random shapes. Then, from a given number of easy-to-generate synthetic cluster configurations, we choose the one with the lowest value of the cost function defined by us using extended ED. At the second stage, we make a significant change in the standard technique of simulated annealing (SA). Instead of swapping pixels of different phases, we randomly move each of the selected synthetic clusters. To demonstrate the accuracy of the method, we reconstruct and analyze two-phase microstructures with irregular inclusions of silica in rubber matrix as well as stones in cement paste. The results show that the two-stage reconstruction (TSR) method provides convincing realizations for these complex microstructures. The advantages of TSR include the ease of obtaining synthetic microstructures, very low computational costs, and satisfactory mapping in the statistical context of inclusion shapes. Finally, its simplicity should greatly facilitate independent applications.

Three-dimensional multiscale fusion for porous media on microtomography images of different resolutions

Accurately acquiring the three-dimensional (3D) image of a porous medium is an imperative issue for the prediction of multiple physical properties. Considering the inherent nature of the multiscale pores contained in porous media such as tight sandstones, to completely characterize the pore structure, one needs to scan the microstructure at different resolutions. Specifically, low-resolution (LR) images cover a larger field of view (FOV) of the sample, but are lacking small-scale features, whereas high-resolution (HR) images contain ample information, but sometimes only cover a limited FOV. To address this issue, we propose a method for fusing the spatial information from a two-dimensional (2D) HR image into a 3D LR image, and finally reconstructing an integrated 3D structure with added fine-scale features. In the fusion process, the large-scale structure depicted by the 3D LR image is fixed as background and the 2D image is utilized as training image to reconstruct a small-scale structure based on the background. To assess the performance of our method, we test it on a sandstone scanned with low and high resolutions. Statistical properties between the reconstructed image and the target are quantitatively compared. The comparison indicates that the proposed method enables an accurate fusion of the LR and HR images because the small-scale information is precisely reproduced within the large one.

Reconstruction of porous media from extremely limited information using conditional generative adversarial networks

Porous media are ubiquitous in both nature and engineering applications. Therefore, their modeling and understanding is of vital importance. In contrast to direct acquisition of three-dimensional (3D) images of this type of medium, obtaining its subregion (s) such as 2D images or several small areas can be feasible. Therefore, reconstructing whole images from limited information is a primary technique in these types of cases. Given data in practice cannot generally be determined by users and may be incomplete or only partially informed, thus making existing reconstruction methods inaccurate or even ineffective. To overcome this shortcoming, in this study we propose a deep-learning-based framework for reconstructing full images from their much smaller subareas. In particular, conditional generative adversarial network is utilized to learn the mapping between the input (a partial image) and output (a full image). To ensure the reconstruction accuracy, two simple but effective objective functions are proposed and then coupled with the other two functions to jointly constrain the training procedure. Because of the inherent essence of this ill-posed problem, a Gaussian noise is introduced for producing reconstruction diversity, thus enabling the network to provide multiple candidate outputs. Our method is extensively tested on a variety of porous materials and validated by both visual inspection and quantitative comparison. It is shown to be accurate, stable, and even fast (∼0.08 s for a 128×128 image reconstruction). The proposed approach can be readily extended by, for example, incorporating user-defined conditional data and an arbitrary number of object functions into reconstruction, while being coupled with other reconstruction methods.

Correlation-driven direct sampling method for geostatistical simulation and training image evaluation

Multiple-point geostatistics (MPS) is a competitive algorithm that produces a set of geologically realistic models. Viewing a training image (TI) as a prior model, MPS extracts patterns from the TI and reproduces patterns which are compatible with the hard data (HD). However, two challenges within the MPS framework are the geologically complex simulation and the TI evaluation. With the objective to achieve a high-quality simulation, we explore a way to address these two issues. First, correlation-driven direct sampling (CDS) is proposed to realize geostatistical simulation. We introduce the correlation-driven distance as a measure of similarity between two patterns. The weights in our distance measurement are derived by the concepts of the ellipse, correlation coefficient, Gaussian function, and affine transformation. Second, we fulfill TI evaluation on the basis of the consistency between TI and HD. Inspired by CDS, the minimum correlation-driven distance (MCD) is proposed to improve the evaluation accuracy. We suggest a conditioning pattern extraction history strategy to speed up the evaluation program. Third, the local consistency is presented to address nonstationarity. The program automatically divides the simulation domain into several subareas. A two-dimensional (2D) channelized reservoir image and a three-dimensional (3D) rock image are used to validate our proposed method. In comparison with previous methods, CDS yields better simulation quality. The further applications include a set of 2D TI evaluations and a 3D simulation domain segmentation. MCD exhibits sensible evaluation accuracy, excellent computational efficiency, and the ability to deal with nonstationarity.

A core-scale reconstructing method for shale

Characterization of shale cores with low and anisotropic permeability is complicated, due to the presence of multiscale pore structure and thin layers, and defies conventional methods. To accurately reproduce the morphology of multiscale pore structure of the shale core, a novel core-scale reconstructing method is proposed to reconstruct 3D digital-experimental models by means of the combination of SEM, EDS images, nitrogen adsorption and pressure pulse decay experiment result. In this method, the multiscale and multicomponent reconstructing algorithm is introduced to build the representative multiscale model for each layer, which can describe the complex 3D structures of nano organic pores, micro-nano inorganic pores, micro slits and several typical minerals. Especially, to reproduce the realistic morphology for shale, the optimization algorithm based on simulated annealing algorithm uses the experimental data as constrain conditions to adjust and optimize the model for each layer. To describe the bedding characteristics of the shale core, bedding fractures are constructed by analysis of the mineral distribution in the interface of two layers, and then the representative models for different layers are integrated together to obtain the final core-scale digital-experimental model. Finally, the model is validated by computing its morphological and flow properties and comparing them with those of the actual 3D shale sample. This method provide a way for systematically and continuously describe the multiscale and anisotropic pore structure (from nm-cm) of the shale core, and will be helpful for understanding the quality of the shale reservoir.

Improved multipoint statistics method for reconstructing three-dimensional porous media from a two-dimensional image via porosity matching

Reconstructing a three-dimensional (3D) structure from a single two-dimensional training image (TI) is a challenging issue. Multiple-point statistics (MPS) is an effective method to solve this problem. However, in the traditional MPS method, errors occur while statistical features of reconstruction, such as porosity, connectivity, and structural properties, deviate from those of TI. Due to the MPS reconstruction mechanism that the voxel being reconstructed is dependent on the reconstructed voxel, it may cause error accumulation during simulations, which can easily lead to a significant difference between the real 3D structure and the reconstructed result. To reduce error accumulation and improve morphological similarity, an improved MPS method based on porosity matching is proposed. In the reconstruction, we search the matching pattern in the TI directly. Meanwhile, a multigrid approach is also applied to capture the large-scale structures of the TI. To demonstrate its superiority over the traditional MPS method, our method is tested on different sandstone samples from many aspects, including accuracy, stability, generalization, and flow characteristics. Experimental results show that the reconstruction results by the improved MPS method effectively match the CT sandstone samples in correlation functions, local porosity distribution, morphological parameters, and permeability.

Markov prior-based block-matching algorithm for superdimension reconstruction of porous media

A superdimension reconstruction algorithm is used for the reconstruction of three-dimensional (3D) structures of a porous medium based on a single two-dimensional image. The algorithm borrows the concepts of "blocks," "learning," and "dictionary" from learning-based superresolution reconstruction and applies them to the 3D reconstruction of a porous medium. In the neighborhood-matching process of the conventional superdimension reconstruction algorithm, the Euclidean distance is used as a criterion, although it may not really reflect the structural correlation between adjacent blocks in an actual situation. Hence, in this study, regular items are adopted as prior knowledge in the reconstruction process, and a Markov prior-based block-matching algorithm for superdimension reconstruction is developed for more accurate reconstruction. The algorithm simultaneously takes into consideration the probabilistic relationship between the already reconstructed blocks in three different perpendicular directions (x, y, and z) and the block to be reconstructed, and the maximum value of the probability product of the blocks to be reconstructed (as found in the dictionary for the three directions) is adopted as the basis for the final block selection. Using this approach, the problem of an imprecise spatial structure caused by a point simulation can be overcome. The problem of artifacts in the reconstructed structure is also addressed through the addition of hard data and by neighborhood matching. To verify the improved reconstruction accuracy of the proposed method, the statistical and morphological features of the results from the proposed method and traditional superdimension reconstruction method are compared with those of the target system. The proposed superdimension reconstruction algorithm is confirmed to enable a more accurate reconstruction of the target system while also eliminating artifacts.

Accelerating simulation for the multiple-point statistics algorithm using vector quantization

Multiple-point statistics (MPS) is a prominent algorithm to simulate categorical variables based on a sequential simulation procedure. Assuming training images (TIs) as prior conceptual models, MPS extracts patterns from TIs using a template and records their occurrences in a database. However, complex patterns increase the size of the database and require considerable time to retrieve the desired elements. In order to speed up simulation and improve simulation quality over state-of-the-art MPS methods, we propose an accelerating simulation for MPS using vector quantization (VQ), called VQ-MPS. First, a variable representation is presented to make categorical variables applicable for vector quantization. Second, we adopt a tree-structured VQ to compress the database so that stationary simulations are realized. Finally, a transformed template and classified VQ are used to address nonstationarity. A two-dimensional (2D) stationary channelized reservoir image is used to validate the proposed VQ-MPS. In comparison with several existing MPS programs, our method exhibits significantly better performance in terms of computational time, pattern reproductions, and spatial uncertainty. Further demonstrations consist of a 2D four facies simulation, two 2D nonstationary channel simulations, and a three-dimensional (3D) rock simulation. The results reveal that our proposed method is also capable of solving multifacies, nonstationarity, and 3D simulations based on 2D TIs.

Reconstruction of a digital core containing clay minerals based on a clustering algorithm

It is difficult to obtain a core sample and information for digital core reconstruction of mature sandstone reservoirs around the world, especially for an unconsolidated sandstone reservoir. Meanwhile, reconstruction and division of clay minerals play a vital role in the reconstruction of the digital cores, although the two-dimensional data-based reconstruction methods are specifically applicable as the microstructure reservoir simulation methods for the sandstone reservoir. However, reconstruction of clay minerals is still challenging from a research viewpoint for the better reconstruction of various clay minerals in the digital cores. In the present work, the content of clay minerals was considered on the basis of two-dimensional information about the reservoir. After application of the hybrid method, and compared with the model reconstructed by the process-based method, the digital core containing clay clusters without the labels of the clusters' number, size, and texture were the output. The statistics and geometry of the reconstruction model were similar to the reference model. In addition, the Hoshen-Kopelman algorithm was used to label various connected unclassified clay clusters in the initial model and then the number and size of clay clusters were recorded. At the same time, the K-means clustering algorithm was applied to divide the labeled, large connecting clusters into smaller clusters on the basis of difference in the clusters' characteristics. According to the clay minerals' characteristics, such as types, textures, and distributions, the digital core containing clay minerals was reconstructed by means of the clustering algorithm and the clay clusters' structure judgment. The distributions and textures of the clay minerals of the digital core were reasonable. The clustering algorithm improved the digital core reconstruction and provided an alternative method for the simulation of different clay minerals in the digital cores.

Evaluating the morphological completeness of a training image

Understanding the three-dimensional (3D) stochastic structure of a porous medium is helpful for studying its physical properties. A 3D stochastic structure can be reconstructed from a two-dimensional (2D) training image (TI) using mathematical modeling. In order to predict what specific morphology belonging to a TI can be reconstructed at the 3D orthogonal slices by the method of 3D reconstruction, this paper begins by introducing the concept of orthogonal chords. After analyzing the relationship among TI morphology, orthogonal chords, and the 3D morphology of orthogonal slices, a theory for evaluating the morphological completeness of a TI is proposed for the cases of three orthogonal slices and of two orthogonal slices. The proposed theory is evaluated using four TIs of porous media that represent typical but distinct morphological types. The significance of this theoretical evaluation lies in two aspects: It allows special morphologies, for which the attributes of a TI can be reconstructed at a special orthogonal slice of a 3D structure, to be located and quantified, and it can guide the selection of an appropriate reconstruction method for a special TI.

Pattern density function for reconstruction of three-dimensional porous media from a single two-dimensional image

Three-dimensional (3D) structures are useful for studying the spatial structures and physical properties of porous media. A 3D structure can be reconstructed from a single two-dimensional (2D) training image (TI) by using mathematical modeling methods. Among many reconstruction algorithms, an optimal-based algorithm was developed and has strong stability. However, this type of algorithm generally uses an autocorrelation function (which is unable to accurately describe the morphological features of porous media) as its objective function. This has negatively affected further research on porous media. To accurately reconstruct 3D porous media, a pattern density function is proposed in this paper, which is based on a random variable employed to characterize image patterns. In addition, the paper proposes an original optimal-based algorithm called the pattern density function simulation; this algorithm uses a pattern density function as its objective function, and adopts a multiple-grid system. Meanwhile, to address the key point of algorithm reconstruction speed, we propose the use of neighborhood statistics, the adjacent grid and reversed phase method, and a simplified temperature-controlled mechanism. The pattern density function is a high-order statistical function; thus, when all grids in the reconstruction results converge in the objective functions, the morphological features and statistical properties of the reconstruction results will be consistent with those of the TI. The experiments include 2D reconstruction using one artificial structure, and 3D reconstruction using battery materials and cores. Hierarchical simulated annealing and single normal equation simulation are employed as the comparison algorithms. The autocorrelation function, linear path function, and pore network model are used as the quantitative measures. Comprehensive tests show that 3D porous media can be reconstructed accurately from a single 2D training image by using the method proposed in this paper.