Nested multiresolution hierarchical simulated annealing algorithm for porous media reconstruction

Research output: Contribution to journalArticlepeer-review

Standard

Nested multiresolution hierarchical simulated annealing algorithm for porous media reconstruction. / Lemmens, Laurent; Rogiers, Bart; Jacques, Diederik; Huysmans, Marijke; Swennen, Rudy; Urai, Janos L.; Desbois, Guillaume; Laloy, Eric.

In: Physical Review E, Vol. 100, No. 5, 053316, 01.11.2019.

Research output: Contribution to journalArticlepeer-review

Author

Lemmens, Laurent ; Rogiers, Bart ; Jacques, Diederik ; Huysmans, Marijke ; Swennen, Rudy ; Urai, Janos L. ; Desbois, Guillaume ; Laloy, Eric. / Nested multiresolution hierarchical simulated annealing algorithm for porous media reconstruction. In: Physical Review E. 2019 ; Vol. 100, No. 5.

Bibtex - Download

@article{df58a4b17ada420599b89a6fb6701ed7,
title = "Nested multiresolution hierarchical simulated annealing algorithm for porous media reconstruction",
abstract = "Microstructure strongly influences flow and transport properties of porous media. Flow and transport simulations within porous media, therefore, requires accurate three-dimensional (3D) models of the pore and solid phase structure. To date, no imaging method can resolve all relevant heterogeneities from the nano- to the centimeter scale within complex heterogeneous materials such as clay, reservoir rocks (e.g., travertine, chalk, ...), hardened cement paste, and concrete. To reconstruct these porous materials it is thus necessary to merge information from different 2D and potentially 3D imaging methods. One porous media reconstruction methodology that has been around for at least two decades is simulated annealing (SA). However, realizations with SA typically suffer an artificially reduced long-range connectivity, while multiphase reconstructions are not feasible in most cases because of a prohibitive computational burden. To solve these problems we propose a hierarchical multiresolution and multiphase simulated annealing algorithm. To decrease the computational cost of multiphase simulation, our algorithm sequentially simulates one phase after another, in a hierarchical way, which enables handling multimodal distributions and topological relations. Building upon recent work, our algorithm improves long-range connectivity and CPU efficiency by simulating larger particles using a coarser resolution that is subsequently refined compared to standard SA; our proposed extension not only offers the possibility to perform multiphase reconstruction but also allows us (i) to improve binary reconstruction quality, as quantified, e.g., by multiple-point histograms by up to one order of magnitude and (ii) to achieve an overall speed-up. The proposed algorithm is also shown to outperform the direct sampling multiple-point statistics method for the generation of cement paste microstructure with respect to both generation time and quality.",
keywords = "porous media",
author = "Laurent Lemmens and Bart Rogiers and Diederik Jacques and Marijke Huysmans and Rudy Swennen and Urai, {Janos L.} and Guillaume Desbois and Eric Laloy",
note = "Score=10",
year = "2019",
month = nov,
day = "1",
doi = "10.1103/PhysRevE.100.053316",
language = "English",
volume = "100",
journal = "Physical Review E",
issn = "2470-0045",
publisher = "APS - American Physical Society",
number = "5",

}

RIS - Download

TY - JOUR

T1 - Nested multiresolution hierarchical simulated annealing algorithm for porous media reconstruction

AU - Lemmens, Laurent

AU - Rogiers, Bart

AU - Jacques, Diederik

AU - Huysmans, Marijke

AU - Swennen, Rudy

AU - Urai, Janos L.

AU - Desbois, Guillaume

AU - Laloy, Eric

N1 - Score=10

PY - 2019/11/1

Y1 - 2019/11/1

N2 - Microstructure strongly influences flow and transport properties of porous media. Flow and transport simulations within porous media, therefore, requires accurate three-dimensional (3D) models of the pore and solid phase structure. To date, no imaging method can resolve all relevant heterogeneities from the nano- to the centimeter scale within complex heterogeneous materials such as clay, reservoir rocks (e.g., travertine, chalk, ...), hardened cement paste, and concrete. To reconstruct these porous materials it is thus necessary to merge information from different 2D and potentially 3D imaging methods. One porous media reconstruction methodology that has been around for at least two decades is simulated annealing (SA). However, realizations with SA typically suffer an artificially reduced long-range connectivity, while multiphase reconstructions are not feasible in most cases because of a prohibitive computational burden. To solve these problems we propose a hierarchical multiresolution and multiphase simulated annealing algorithm. To decrease the computational cost of multiphase simulation, our algorithm sequentially simulates one phase after another, in a hierarchical way, which enables handling multimodal distributions and topological relations. Building upon recent work, our algorithm improves long-range connectivity and CPU efficiency by simulating larger particles using a coarser resolution that is subsequently refined compared to standard SA; our proposed extension not only offers the possibility to perform multiphase reconstruction but also allows us (i) to improve binary reconstruction quality, as quantified, e.g., by multiple-point histograms by up to one order of magnitude and (ii) to achieve an overall speed-up. The proposed algorithm is also shown to outperform the direct sampling multiple-point statistics method for the generation of cement paste microstructure with respect to both generation time and quality.

AB - Microstructure strongly influences flow and transport properties of porous media. Flow and transport simulations within porous media, therefore, requires accurate three-dimensional (3D) models of the pore and solid phase structure. To date, no imaging method can resolve all relevant heterogeneities from the nano- to the centimeter scale within complex heterogeneous materials such as clay, reservoir rocks (e.g., travertine, chalk, ...), hardened cement paste, and concrete. To reconstruct these porous materials it is thus necessary to merge information from different 2D and potentially 3D imaging methods. One porous media reconstruction methodology that has been around for at least two decades is simulated annealing (SA). However, realizations with SA typically suffer an artificially reduced long-range connectivity, while multiphase reconstructions are not feasible in most cases because of a prohibitive computational burden. To solve these problems we propose a hierarchical multiresolution and multiphase simulated annealing algorithm. To decrease the computational cost of multiphase simulation, our algorithm sequentially simulates one phase after another, in a hierarchical way, which enables handling multimodal distributions and topological relations. Building upon recent work, our algorithm improves long-range connectivity and CPU efficiency by simulating larger particles using a coarser resolution that is subsequently refined compared to standard SA; our proposed extension not only offers the possibility to perform multiphase reconstruction but also allows us (i) to improve binary reconstruction quality, as quantified, e.g., by multiple-point histograms by up to one order of magnitude and (ii) to achieve an overall speed-up. The proposed algorithm is also shown to outperform the direct sampling multiple-point statistics method for the generation of cement paste microstructure with respect to both generation time and quality.

KW - porous media

UR - http://ecm.sckcen.be/OTCS/llisapi.dll/open/36530147

U2 - 10.1103/PhysRevE.100.053316

DO - 10.1103/PhysRevE.100.053316

M3 - Article

VL - 100

JO - Physical Review E

JF - Physical Review E

SN - 2470-0045

IS - 5

M1 - 053316

ER -

ID: 5843651