Chen Dongdong, Teng Qizhi, He Xiaohai, Xu Zhi, Li Zhengji
College of Electronics and Information Engineering, Sichuan University, Chengdu 610065, China.
Phys Rev E Stat Nonlin Soft Matter Phys. 2014 Jan;89(1):013305. doi: 10.1103/PhysRevE.89.013305. Epub 2014 Jan 17.
In this paper, we introduce a stable-phase approach for hierarchical annealing which addresses the very large computational costs associated with simulated annealing for the reconstruction of large-scale binary porous media images. Our presented method, which uses the two-point correlation function as the morphological descriptor, involves the reconstruction of three-phase and two-phase structures. We consider reconstructing the three-phase structures based on standard annealing and the two-phase structures based on standard and hierarchical annealings. From the result of the two-dimensional (2D) reconstruction, we find that the 2D generation does not fully capture the morphological information of the original image, even though the two-point correlation function of the reconstruction is in excellent agreement with that of the reference image. For the reconstructed three-dimensional (3D) microstructure, we calculate its permeability and compare it to that of the reference 3D microstructure. The result indicates that the reconstructed structure has a lower degree of connectedness than that of the actual sandstone. We also compare the computation time of our presented method to that of the standard annealing, which shows that our presented method of orders of magnitude improves the convergence rate. That is because only a small part of the pixels in the overall hierarchy need to be considered for sampling by the annealer.
在本文中，我们介绍了一种用于分层退火的稳定相方法，该方法解决了与模拟退火相关的巨大计算成本问题，用于重建大规模二元多孔介质图像。我们提出的方法使用两点相关函数作为形态描述符，涉及三相和两相结构的重建。我们考虑基于标准退火重建三相结构，基于标准退火和分层退火重建两相结构。从二维（2D）重建结果来看，我们发现即使重建的两点相关函数与参考图像的两点相关函数非常吻合，二维生成也不能完全捕捉原始图像的形态信息。对于重建的三维（3D）微观结构，我们计算其渗透率并将其与参考3D微观结构的渗透率进行比较。结果表明，重建结构的连通性程度低于实际砂岩的连通性程度。我们还将我们提出的方法的计算时间与标准退火的计算时间进行了比较，结果表明我们提出的方法在数量级上提高了收敛速度。这是因为在整个层次结构中，退火器只需考虑一小部分像素进行采样。
