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

Evaluation of three-point correlation functions from structural images on CPU and GPU architectures: Accounting for anisotropy effects

Structures, or spatial arrangements of matter and energy, including some fields (e.g., velocity or pressure) are ubiquitous in research applications and frequently require description for subsequent analysis, or stochastic reconstruction from limited data. The classical descriptors are two-point correlation functions (CFs), but the computation of three-point statistics is known to be advantageous in some cases as they can probe non-Gaussian signatures, not captured by their two-point counterparts. Moreover, n-point CFs with n≥3 are believed to possess larger information content and provide more information about studied structures. In this paper, we have developed algorithms and code to compute S_{3},C_{3},F_{sss},F_{ssv}, and F_{svv} with a right-angle and arbitrary triangle pattern. The former was believed to be faster to compute, but with the help of precomputed regular positions we achieved the same speed for arbitrary pattern. In this work we also implement and demonstrate computations of directional three-point CFs-for this purpose right-triangular pattern seems to be superior due to explicit orientation and high coverage. Moreover, we assess the errors in CFs' evaluation due to image or pattern rotations and show that they have minor effect on accuracy of computations. The execution times of our algorithms for the same number of samples are orders of magnitude lower than in existing published counterparts. We show that the volume of data produced gets unwieldy very easily, especially if computations are performed in frequency domain. For these reasons until information content of different sets of correlation functions with different "n-pointness" is known, advantages of CFs with n>3 are not clear. Nonetheless, developed algorithms and code are universal enough to be easily extendable to any n with increasing computational and random access memory (RAM) burden. All results are available as part of open-source package correlationfunctions.jl [V. Postnicov et al., Comput. Phys. Commun. 299, 109134 (2024)10.1016/j.cpc.2024.109134.]. As described in this paper, three-point CFs computations can be immediately applied in a great number of research applications, for example: (1) flow and transport velocity fields analysis or any data with non-Gaussian signatures, (2) deep learning for structural and physical properties, and (3) structure taxonomy and categorization. In all these and numerous other potential cases the ability to compute directional three-point functions may be crucial. Notably, the organization of the code functions allows computation of cross correlation, i.e., one can compute three-point CFs for multiphase images (while binary structures were used in this paper for simplicity of explanations).

Measuring structural nonstationarity: The use of imaging information to quantify homogeneity and inhomogeneity

Heterogeneity is the concept we encounter in numerous research areas and everyday life. While "not mixing apples and oranges" is easy to grasp, a more quantitative approach to such segregation is not always readily available. Consider the problem from a different angle: To what extent does one have to make apples more orange and oranges more "apple-shaped" to put them into the same basket (according to their appearance alone)? This question highlights the central problem of the blurred interface between heterogeneous and homogeneous, which also depends on the metrics used for its identification. This work uncovers the physics of structural stationarity quantification, based on correlation functions (CFs) and clustering based on CFs different between image subregions. By applying the methodology to a wide variety of synthetic and real images of binary porous media, we confirmed computationally that only periodically unit-celled structures and images produced by stationary processes with resolutions close to infinity are strictly stationary. Natural structures without recurring unit cells are only weakly stationary. We established a physically meaningful definition for these stationarity types and their distinction from nonstationarity. In addition, the importance of information content of the chosen metrics is highlighted and discussed. We believe the methodology as proposed in this contribution will find its way into numerous research areas dealing with materials, structures, and measurements and modeling based on structural imaging information.

Robust surface-correlation-function evaluation from experimental discrete digital images

Correlation functions (CFs) are universal structural descriptors; surface-surface F_{ss} and surface-void F_{sv} CFs are a subset containing additional information about the interface between the phases. The description of the interface between pores and solids in porous media is of particular importance and recently Ma and Torquato [Phys. Rev. E 98, 013307 (2018)2470-004510.1103/PhysRevE.98.013307] proposed an elegant way to compute these functions for a wide variety of cases. However, their "continuous" approach is not always applicable to digital experimental 2D and 3D images of porous media as obtained using x-ray tomography or scanning electron microscopy due to nonsingularities in chemical composition or local solid material's density and partial volume effects. In this paper we propose to use edge-detecting filters to compute surface CFs in the "digital" fashion directly in the images. Computed this way, surface correlation functions are the same as analytically known for Poisson disks in case the resolution of the image is adequate. Based on the multiscale image analysis we developed a C_{0.5} criterion that can predict if the imaging resolution is enough to make an accurate evaluation of the surface CFs. We also showed that in cases when the input image contains all major features, but do not pass the C_{0.5} criterion, it is possible with the help of image magnification to sample CFs almost similar to those obtained for high-resolution image of the same structure with high C_{0.5}. The computational framework as developed here is open source and available within the CorrelationFunctions.jl package developed by our group. Our "digital" approach was applied to a wide variety of real porous media images of different quality. We discuss critical aspects of surface correlation functions computations as related to different applications. The developed methodology allows applying surface CFs to describe the structure of porous materials based on their experimental images and enhance stochastic reconstructions or super-resolution procedures, or serve as an efficient metrics in machine learning applications due to computationally effective GPU implementation.

Effect of the Agglomerate Geometry on the Effective Electrical Conductivity of a Porous Electrode

The study of the microstructure of random heterogeneous materials, related to an electrochemical device, is relevant because their effective macroscopic properties, e.g., electrical or proton conductivity, are a function of their effective transport coefficients (ETC). The magnitude of ETC depends on the distribution and properties of the material phase. In this work, an algorithm is developed to generate stochastic two-phase (binary) image configurations with multiple geometries and polydispersed particle sizes. The recognizable geometry in the images is represented by the white phase dispersed and characterized by statistical descriptors (two-point and line-path correlation functions). Percolation is obtained for the geometries by identifying an infinite cluster to guarantee the connection between the edges of the microstructures. Finally, the finite volume method is used to determine the ETC. Agglomerate phase results show that the geometry with the highest local current distribution is the triangular geometry. In the matrix phase, the most significant results are obtained by circular geometry, while the lowest is obtained by the 3-sided polygon. The proposed methodology allows to establish criteria based on percolation and surface fraction to assure effective electrical conduction according to their geometric distribution; results provide an insight for the microstructure development with high projection to be used to improve the electrode of a Membrane Electrode Assembly (MEA).

Does Rheology of Bingham Fluid Influence Upscaling of Flow through Tight Porous Media?

Non-Newtonian fluids may cause nonlinear seepage even for a single-phase flow. Through digital rock technologies, the upscaling of this non-Darcy flow can be studied; however, the requirements for scanning resolution and sample size need to be clarified very carefully. This work focuses on Bingham fluid flow in tight porous media by a pore-scale simulation on CT-scanned microstructures of tight sandstones. A bi-viscous model is used to depict the Bingham fluid. The results show that when the Bingham fluid flows through a rock sample, the flowrate increases at a parabolic rate when the pressure gradient is small and then increases linearly with the pressure gradient. As a result, an effective permeability and a start-up pressure gradient can be used to characterize this flow behavior. By conducting flow simulations at varying sample sizes, we obtain the representative element volume (REV) for effective permeability and start-up pressure gradient. It is found that the REV size for the effective permeability is almost the same as that for the absolute permeability of Newtonian fluid. The interesting result is that the REV size for the start-up pressure gradient is much smaller than that for the effective permeability. The results imply that the sample size, which is large enough to reach the REV size for Newtonian fluids, can be used to investigate the Bingham fluids flow through porous media as well.

Computed Tomography 3D Super-Resolution with Generative Adversarial Neural Networks: Implications on Unsaturated and Two-Phase Fluid Flow

Fluid flow characteristics are important to assess reservoir performance. Unfortunately, laboratory techniques are inadequate to know these characteristics, which is why numerical methods were developed. Such methods often use computed tomography (CT) scans as input but this technique is plagued by a resolution versus sample size trade-off. Therefore, a super-resolution method using generative adversarial neural networks (GANs) was used to artificially improve the resolution. Firstly, the influence of resolution on pore network properties and single-phase, unsaturated, and two-phase flow was analysed to verify that pores and pore throats become larger on average and surface area decreases with worsening resolution. These observations are reflected in increasingly overestimated single-phase permeability, less moisture uptake at lower capillary pressures, and high residual oil fraction after waterflooding. Therefore, the super-resolution GANs were developed which take low (12 µm) resolution input and increase the resolution to 4 µm, which is compared to the expected high-resolution output. These results better predicted pore network properties and fluid flow properties despite the overestimation of porosity. Relevant small pores and pore surfaces are better resolved thus providing better estimates of unsaturated and two-phase flow which can be heavily influenced by flow along pore boundaries and through smaller pores. This study presents the second case in which GANs were applied to a super-resolution problem on geological materials, but it is the first one to apply it directly on raw CT images and to determine the actual impact of a super-resolution method on fluid predictions.

Effect of An Image Resolution Change on the Effective Transport Coefficient of Heterogeneous Materials

Electrochemical electrodes comprise multiple phenomena at different scales. Several works have tried to model such phenomena using statistical techniques. This paper proposes a novel process to work with reduced size images to reconstruct microstructures with the Simulated Annealing method. Later, using the Finite Volume Method, it is verified the effect of the image resolution on the effective transport coefficient (ETC). The method can be applied to synthetic images or images from the Scanning Electron Microscope. The first stage consists of obtaining the image of minimum size, which contains at least 98% of the statistical information of the original image, allowing an equivalent statistical study. The image size reduction was made by applying an iterative decimation over the image using the normalized coarseness to compare the amount of information contained at each step. Representative improvements, especially in processing time, are achieved by reducing the size of the reconstructed microstructures without affecting their statistical behavior. The process ends computing the conduction efficiency from the microstructures. The simulation results, obtained from two kinds of images from different materials, demonstrate the effectivity of the proposed approach. It is important to remark that the controlled decimation allows a reduction of the processor and memory use during the reconstruction and ETC computation of electrodes.