Stochastic reconstruction using multiple correlation functions with different-phase-neighbor-based pixel selection

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.

Adaptive phase-retrieval stochastic reconstruction with correlation functions: Three-dimensional images from two-dimensional cuts

Precise characterization of three-dimensional (3D) heterogeneous media is indispensable in finding the relationships between structure and macroscopic physical properties (permeability, conductivity, and others). The most widely used experimental methods (electronic and optical microscopy) provide high-resolution bidimensional images of the samples of interest. However, 3D material inner microstructure registration is needed to apply numerous modeling tools. Numerous research areas search for cheap and robust methods to obtain the full 3D information about the structure of the studied sample from its 2D cuts. In this work, we develop an adaptive phase-retrieval stochastic reconstruction algorithm that can create 3D replicas from 2D original images, APR. The APR is free of artifacts characteristic of previously proposed phase-retrieval techniques. While based on a two-point S_{2} correlation function, any correlation function or other morphological metrics can be accounted for during the reconstruction, thus, paving the way to the hybridization of different reconstruction techniques. In this work, we use two-point probability and surface-surface functions for optimization. To test APR, we performed reconstructions for three binary porous media samples of different genesis: sandstone, carbonate, and ceramic. Based on computed permeability and connectivity (C_{2} and L_{2} correlation functions), we have shown that the proposed technique in terms of accuracy is comparable to the classic simulated annealing-based reconstruction method but is computationally very effective. Our findings open the possibility of utilizing APR to produce fast or crude replicas further polished by other reconstruction techniques such as simulated annealing or process-based methods. Improving the quality of reconstructions based on phase retrieval by adding additional metrics into the reconstruction procedure is possible for future work.

Descriptor-based method combined with partition to reconstruct three-dimensional complex microstructures

A descriptor-based method combined with a partition approach is proposed to reconstruct three-dimensional (3D) microstructures based on a set of two-dimensional (2D) scanning electron microscopy (SEM) images. The features in the SEM images are identified and partitioned into small features using the watershed algorithm. The watershed algorithm first finds the local gray-level maxima, and partitions the features through the gray-level local minima. The 3D size distribution and radial distribution of the small spherical elements are inferred, respectively, based on the 2D size distribution and radial distribution using stereological analysis. The 3D microstructures are reconstructed by matching the inferred size distribution and radial distribution through a simulated annealing-based procedure. Combining with the proposed partition approach, the descriptor-based method can be applied to complex microstructures and the computational efficiency of the reconstruction can be largely improved. A case study is presented using a set of 2D SEM images with nanoscale pore structure from the low-density CSH (calcium silicate hydrate) phase of a hardened cement paste. Cross sections were randomly selected from the reconstructed 3D microstructure and compared with the original SEM images using the pore descriptors and the two-point correlation function with satisfactory agreement. Using the 3D reconstructed model, the properties of the sample material can be investigated on such a small scale as demonstrated in this paper on quantifying the absolute permeability.

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.

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.

Effect of the image resolution on the statistical descriptors of heterogeneous media

The characterization and reconstruction of heterogeneous materials, such as porous media and electrode materials, involve the application of image processing methods to data acquired by scanning electron microscopy or other microscopy techniques. Among them, binarization and decimation are critical in order to compute the correlation functions that characterize the microstructure of the above-mentioned materials. In this study, we present a theoretical analysis of the effects of the image-size reduction, due to the progressive and sequential decimation of the original image. Three different decimation procedures (random, bilinear, and bicubic) were implemented and their consequences on the discrete correlation functions (two-point, line-path, and pore-size distribution) and the coarseness (derived from the local volume fraction) are reported and analyzed. The chosen statistical descriptors (correlation functions and coarseness) are typically employed to characterize and reconstruct heterogeneous materials. A normalization for each of the correlation functions has been performed. When the loss of statistical information has not been significant for a decimated image, its normalized correlation function is forecast by the trend of the original image (reference function). In contrast, when the decimated image does not hold statistical evidence of the original one, the normalized correlation function diverts from the reference function. Moreover, the equally weighted sum of the average of the squared difference, between the discrete correlation functions of the decimated images and the reference functions, leads to a definition of an overall error. During the first stages of the gradual decimation, the error remains relatively small and independent of the decimation procedure. Above a threshold defined by the correlation length of the reference function, the error becomes a function of the number of decimation steps. At this stage, some statistical information is lost and the error becomes dependent on the decimation procedure. These results may help us to restrict the amount of information that one can afford to lose during a decimation process, in order to reduce the computational and memory cost, when one aims to diminish the time consumed by a characterization or reconstruction technique, yet maintaining the statistical quality of the digitized sample.

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.

Microstructure and mechanical properties of hyperuniform heterogeneous materials

A hyperuniform random heterogeneous material is one in which the local volume fraction fluctuations in an observation window decay faster than the reciprocal window volume as the window size increases. Recent studies show that this class of materials are endowed with superior physical properties such as large isotropic photonic band gaps and optimal transport properties. Here we employ a stochastic optimization procedure to systematically generate realizations of hyperuniform heterogeneous materials with controllable short-range order, which is partially quantified using the two-point correlation function S_{2}(r) associated with the phase of interest. Specifically, our procedure generalizes the widely used Yeong-Torquato reconstruction procedure by including an additional constraint for hyperuniformity, i.e., the volume integral of the autocovariance function χ(r)=S_{2}(r)-ϕ^{2} over the whole space is zero. In addition, we only require the reconstructed S_{2} to match the target function up to a certain cutoff distance γ, in order to give the system sufficient degrees of freedom to satisfy the hyperuniform condition. By systematically increasing the γ value for a given S_{2}, one can produce a spectrum of hyperuniform heterogeneous materials with varying degrees of partial short-range order compatible with the specified S_{2}. The mechanical performance including both elastic and brittle fracture behaviors of the generated hyperuniform materials is analyzed using a volume-compensated lattice-particle method. For the purpose of comparison, the corresponding nonhyperuniform materials with the same short-range order (i.e., with S_{2} constrained up to the same γ value) are also constructed and their mechanical performance is analyzed. Here we consider two specific S_{2} including the positive exponential decay function and the correlation function associated with an equilibrium hard-sphere system. For the constructed systems associated with these two specific functions, we find that although the hyperuniform materials are softer than their nonhyperuniform counterparts, the former generally possess a significantly higher brittle fracture strength than the latter. This superior mechanical behavior is attributed to the lower degree of stress concentration in the material resulting from the hyperuniform microstructure, which is crucial to crack initiation and propagation.

Multigrid hierarchical simulated annealing method for reconstructing heterogeneous media

A reconstruction methodology based on different-phase-neighbor (DPN) pixel swapping and multigrid hierarchical annealing is presented. The method performs reconstructions by starting at a coarse image and successively refining it. The DPN information is used at each refinement stage to freeze interior pixels of preformed structures. This preserves the large-scale structures in refined images and also reduces the number of pixels to be swapped, thereby resulting in a decrease in the necessary computational time to reach a solution. Compared to conventional single-grid simulated annealing, this method was found to reduce the required computation time to achieve a reconstruction by around a factor of 70-90, with the potential of even higher speedups for larger reconstructions. The method is able to perform medium sized (up to 300(3) voxels) three-dimensional reconstructions with multiple correlation functions in 36-47 h.

Universal Stochastic Multiscale Image Fusion: An Example Application for Shale Rock

Spatial data captured with sensors of different resolution would provide a maximum degree of information if the data were to be merged into a single image representing all scales. We develop a general solution for merging multiscale categorical spatial data into a single dataset using stochastic reconstructions with rescaled correlation functions. The versatility of the method is demonstrated by merging three images of shale rock representing macro, micro and nanoscale spatial information on mineral, organic matter and porosity distribution. Merging multiscale images of shale rock is pivotal to quantify more reliably petrophysical properties needed for production optimization and environmental impacts minimization. Images obtained by X-ray microtomography and scanning electron microscopy were fused into a single image with predefined resolution. The methodology is sufficiently generic for implementation of other stochastic reconstruction techniques, any number of scales, any number of material phases, and any number of images for a given scale. The methodology can be further used to assess effective properties of fused porous media images or to compress voluminous spatial datasets for efficient data storage. Practical applications are not limited to petroleum engineering or more broadly geosciences, but will also find their way in material sciences, climatology, and remote sensing.

Dynamic reconstruction of heterogeneous materials and microstructure evolution

Reconstructing heterogeneous materials from limited structural information has been a topic that attracts extensive research efforts and still poses many challenges. The Yeong-Torquato procedure is one of the most popular reconstruction techniques, in which the material reconstruction problem based on a set of spatial correlation functions is formulated as a constrained energy minimization (optimization) problem and solved using simulated annealing [Yeong and Torquato, Phys. Rev. E 57, 495 (1998)]. The standard two-point correlation function S2 has been widely used in reconstructions, but can also lead to large structural degeneracy for certain nearly percolating systems. To improve reconstruction accuracy and reduce structural degeneracy, one can successively incorporate additional morphological information (e.g., nonconventional or higher-order correlation functions), which amounts to reshaping the energy landscape to create a deep (local) energy minimum. In this paper, we present a dynamic reconstruction procedure that allows one to use a series of auxiliary S2 to achieve the same level of accuracy as those incorporating additional nonconventional correlation functions. In particular, instead of randomly sampling the microstructure space as in the simulated annealing scheme, our procedure utilizes a series of auxiliary microstructures that mimic a physical structural evolution process (e.g., grain growth). This amounts to constructing a series auxiliary energy landscapes that bias the convergence of the reconstruction to a favored (local) energy minimum. Moreover, our dynamic procedure can be naturally applied to reconstruct an actual microstructure evolution process. In contrast to commonly used evolution reconstruction approaches that separately generate individual static configurations, our procedure continuously evolves a single microstructure according to a time-dependent correlation function. The utility of our procedure is illustrated by successfully reconstructing nearly percolating hard-sphere packings and particle-reinforced composites as well as the coarsening process in a binary metallic alloy.