Physically Consistent 2D to 3D Pore Space Synthesis: 1. Dynamic Grain Packing Under Strict Morphology Constraints
Abstract
The 3D digital reconstruction of multi-scale, heterogeneous porous media is traditionally limited by the resolution boundaries and artifacts of X-ray computed tomography (XCT). This paper introduces a physically consistent 2D to 3D pore space synthesis framework that supersedes conventional geometric shuffling with dynamic Newtonian gravitational deposition using a high-performance rigid-body mechanics engine. To honor strict physical contact mechanics, statistical 2D targets, and morphological constraints, we develop a two-stage dynamic packing pipeline that combines realistic 3D grain shapes with a virtual sub-particle insertion engine. The framework is validated against synthetic and natural core benchmarks. First, we show that our stereological integral equation systematically converts size distributions; bypassing this inversion step induces a systematic 50\% deviation in permeability predictions despite identical total porosities, independently altering flow channel topology. Second, temporary virtual spheres serve as architectural placeholders to stabilize large non-equilibrium dissolution voids, ensuring the reconstructed model satisfies target flow metrics within the prescribed tolerance. Finally, when applied to a fine-grained, weakly consolidated sandstone, the workflow exposes and suppresses systematic XCT artifacts, including phantom internal porosity. Ensembles of two-point correlation functions demonstrate that the hydrodynamically and structurally synthesized digital twins quantitatively exhibit closer structural proximity to the uncorrupted rock geometry than the reference XCT dataset itself. This methodology establishes a scalable foundation for high-fidelity multi-phase rock reconstruction, with potential extensions toward complex amorphous cement phases.
Related
- Physically constrained Mamba forecasting supports multiobjective reservoir water quality decisions
- WP-MIP: An Artificial Intelligence, Hybrid, and Physically Based Model Intercomparison Project for Weather Prediction
- Propagating meteorological uncertainty in physically based mountain snow modeling
- A physically integrated GNN surrogate for microbially-mediated kinetic reactive transport: PHYGNET
- Physically crosslinked hydrogel coatings for tunable phosphorus release and enhanced wheat growth
BibTeX
@misc{kulygin2026physically,
title = {{Physically Consistent 2D to 3D Pore Space Synthesis: 1. Dynamic Grain Packing Under Strict Morphology Constraints}},
author = {Dmitry Kulygin and Andrey Ananev and Aleksey Khlyupin},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.13020},
eprint = {2609.13020},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.13020},
}
RIS
TY - UNPB
TI - Physically Consistent 2D to 3D Pore Space Synthesis: 1. Dynamic Grain Packing Under Strict Morphology Constraints
AU - Dmitry Kulygin
AU - Andrey Ananev
AU - Aleksey Khlyupin
PY - 2026
DA - 2026/09/11
PB - arXiv
DO - 10.48550/arxiv.2609.13020
UR - https://arxiv.org/abs/2609.13020
AB - The 3D digital reconstruction of multi-scale, heterogeneous porous media is traditionally limited by the resolution boundaries and artifacts of X-ray computed tomography (XCT). This paper introduces a physically consistent 2D to 3D pore space synthesis framework that supersedes conventional geometric shuffling with dynamic Newtonian gravitational deposition using a high-performance rigid-body mechanics engine. To honor strict physical contact mechanics, statistical 2D targets, and morphological constraints, we develop a two-stage dynamic packing pipeline that combines realistic 3D grain shapes with a virtual sub-particle insertion engine. The framework is validated against synthetic and natural core benchmarks. First, we show that our stereological integral equation systematically converts size distributions; bypassing this inversion step induces a systematic 50\% deviation in permeability predictions despite identical total porosities, independently altering flow channel topology. Second, temporary virtual spheres serve as architectural placeholders to stabilize large non-equilibrium dissolution voids, ensuring the reconstructed model satisfies target flow metrics within the prescribed tolerance. Finally, when applied to a fine-grained, weakly consolidated sandstone, the workflow exposes and suppresses systematic XCT artifacts, including phantom internal porosity. Ensembles of two-point correlation functions demonstrate that the hydrodynamically and structurally synthesized digital twins quantitatively exhibit closer structural proximity to the uncorrupted rock geometry than the reference XCT dataset itself. This methodology establishes a scalable foundation for high-fidelity multi-phase rock reconstruction, with potential extensions toward complex amorphous cement phases.
ER -