Search NASASearch

DOE OSTI · 3010747

A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials

Abstract

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Liu, Mengwei [National Energy Technology Laboratory (NETL), Morgantown, WV (United States)] (ORCID:0000000311965222), Seol, Yongkoo [National Energy Technology Laboratory (NETL), Morgantown, WV (United States)] (ORCID:0000000295163405). 2025-12-15. A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials. https://doi.org/10.1016/j.advwatres.2025.105197

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

Isolated vesicles in submarine pumice: insights from the 2019 Volcano F eruption, Kingdom of Tonga

No reliable method can distinguish between subaerial and submarine pumice, but recent work shows that highly vesicular (total porosity, φ > 65%) submarine pumice from the 2012 Havre and 2021 Fukutoku Oka-no-Ba eruptions contains abundant isolated porosity (connectivity, C < 0.6). This differs from textural measurements of subaerial pumice with similarly high vesicularities, where high connectivities (C > 0.8) are often measured. To investigate the implications of this under-studied texture for bubble nucleation, growth, and coalescence dynamics at submarine eruptions, we analyze the textural characteristics of rafted pumice clasts from the 2019 submarine eruption of Volcano F, located along the Tonga-Kermadec Arc. We examined clasts collected while floating near the vent (50–200 km) and distally on the shores of Fiji (900 km). Using helium pycnometry and X-ray tomography (XRT), we quantified porosity and connectivity in 45 pumice lapilli, classifying them by vesicle macrotexture and edge morphology. Microvesicular clasts, the majority of which have breadcrust and cauliflower textures, exhibit the lowest connectivities (C = 0.61 ± 0.2), while macrovesicular clasts have high connectivities with C ~ 1. We propose that isolated porosity abundance varies according to clast location in the eruption column post-fragmentation. Microvesicular clasts erupted on the edges of the column and experienced high levels of seawater ingestion and subsequent quenching of the entire clast. Macrovesicular clasts, in contrast, erupted in the center of the eruption column and were thus thermally insulated, allowing magmatic bubbles in the molten pumice to nucleate, grow, and coalesce for longer. Our results and interpretation imply that low connectivity, microvesicular clasts preserve pre and syn-eruptive bubble processes at the 2019 Volcano F eruption.

Permeability