InterPore Journal
type: channel
description: InterPore Journal is a new international, peer-reviewed, open-access journal, fully owned and published by the International Society for Porous Media (InterPore). Like the Society, our mission is to advance knowledge and research on porous media systems and their technical and industrial applications. We provide a trusted medium for sharing insights into the processes and phenomena involved, bringing together contributions from across disciplines to strengthen the global porous media community. The journal publishes original research articles, technical notes, review papers, open datasets and source codes, as well as short communications, commentaries, and opinion pieces. Topics span the full breadth of porous media science and technology, including imaging and characterization, physical and numerical modeling, data-driven methods, material design, and applications in geology, industry, environment, and biology. Our overarching goal is to serve as the focal point for the worldwide porous media community by publishing and freely sharing high-quality, multidisciplinary, and multiscale research. Articles are published under a Creative Commons license and appear four times per year with the online release of each issue, following final publication approval by the authors. Information on all published articles is shared through our bi-weekly e-newsletters (reaching 3000+ recipients) and across our social media channels. Stay Connected Keep up with the latest publications, news, and updates from InterPore Journal by subscribing to our newsletter and following us on social media. Join our vibrant online community and engage with fellow researchers worldwide: LinkedIn X YouTube Submit your manuscript today and join us in advancing the field of porous media science and technology. We look forward to receiving your contributions and continuing to support the growth of knowledge in 2025 and beyond.
publisher: InterPore: International Society for Porous Media
language: en
publication name: InterPore Journal
issn: 3105-6369
copyright: Unless otherwise stated above, this is an open access article published by InterPore under either the terms of the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0) (https://creativecommons.org/licenses/by-nc-nd/4.0/). Article metadata are available under the CCo license.
items: Selected Porous Media Contributions from the 1st Iran InterPore Conference
items: A Novel Surfactant-Enhanced DNAPL Remediation Technology with Horizontal Wells
items: Unraveling salt precipitation mechanisms in geological CO2 storage: insights into dominant driving forces
items: Pareto-based Multi-Objective Optimization of Underground Hydrogen Storage Operational Condition Utilizing Machine Learning-Based Surrogate Models
items: Highlights from the new issue of InterPore Journal
items: A Comparative SPH Analysis of Interaction Between Viscous Debris Flow and Two Types of Barriers with Basal Clearance
items: Design of a Laboratory Setup to Determine the Dispersion Parameter in a Single-Well Chemical Tracer Test
items: Multiphase Image Segmentation of Naturally Fractured Media: Benchmarking Deep Learning and Conventional Approaches
items: Modelling copper mobilisation from low-grade ore controlled by mineral surface reactions
items: Exploring the Inter-Relationship between Climatic Conditions and Heat-Water Transfer in Deformable Cracked Surface Strata
items: MAGNET: Medial Axis Guided Network Extraction Tool
items: Heat-activated Epoxy Foam for Permeability Alteration in Fractured Geothermal Fields: Proof of Concept
items: Optimization of the Energy Extraction Process from Enhanced Geothermal Systems
items: Effect of fissure drainage on dispersive mixing in porous media gravity currents: an experimental investigation
items: Petrophysical Characterization of Indiana Limestone Using Medical Dual-Energy Computed Tomography Technique: Insights into Porosity, Bulk Density, and Effective Atomic Number
items: Suitability Mapping for Subsurface Floodwater Storage Schemes
items: Static Elastic Properties of Rocks Obtained by X-ray Microtomography, Inverse Modeling and Surrogate Model
items: A Hydro-Geomechanical Porous-Media Model to Study Effects of Engineered Carbonate Precipitation in Faults
items: On the Interpretation of Unsteady State Experiments in Heterogeneous Rock by Stochastic Methods
items: Experimental Investigation of Solubility Trapping in 3D Printed Micromodels
items: 3D Printing Reactive Porous Media: Calcite Precipitation Kinetics on Surface Functionalized Polymer Films and 3D-printed Cores
items: Experimental Study of Microbial Hydrogen Consumption Rates by Oleidesulfovibrio Alaskensis in Porous Media
items: Regional-Scale Subsurface Modeling
items: Highlights From the New Issue of InterPore Journal
items: Understanding Fracture Aperture and Permeability Evolution due to Carbonate Mineralization Utilizing 3D Printing
items: A novel robust optimization framework based on surrogate modeling for underground hydrogen storage in depleted natural gas reservoirs
items: Rubble Characteristics Associated with Room Collapse at the Waste Isolation Pilot Plant: Impact of Salt Clast Shapes and Size Distributions on the Depositional Pore System
items: Anchored Physics-Informed Neural Network for Two-Phase Flow Simulation in Heterogeneous Porous Media
items: Highlights From the Fourth issue of InterPore Journal

Selected Porous Media Contributions from the 1st Iran InterPore Conference
type: item
subject: Editorial
  topic: section
subject: Porous Media
  topic: keywords
subject: Iran InterPore Conference
  topic: keywords
subject: Editorial
  topic: keywords
subject: Editorial
  topic: subjects
creator: Hassan Mahani
rights: Copyright (c) 2026 Hassan Mahani https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-01
publication date: 2026-03-01
starting page: IPJ010326
ending page: 1
doi: 10.69631/fhmwhz36

A Novel Surfactant-Enhanced DNAPL Remediation Technology with Horizontal Wells
type: item
description: This study investigates surfactant-enhanced in-situ DNAPL (Dense Non-Aqueous Phase Liquids) remediation methods by using both vertical and horizontal wells. This novel technology is examined in two steps: (i) matching the pattern of DNAPL plume leakage and spreading and (ii) predicting the remediation efficiency in various design and optimization scenarios when the surfactant solution is continuously imbibed from the upper vertical well and the formation fluids (containing DNAPLs) are extracted from the lower horizontal wells. A dimensionless number comparing two capillary pressure values is introduced to understand pore-scale governing displacement mechanisms, which is combined with the concept of large-scale sweep efficiency. By positioning these extraction wells beneath the contaminant plume with the surfactant solution introduced from the top, the simulation results show that the system configuration and remediation process successfully create a sufficient pressure gradient to overcome capillary trapping and mobilize oil blobs, generating enough contact areas to sweep the contaminants. The results also demonstrate that the surfactant-enhanced remediation technique, if properly designed and implemented with horizontal wells, can achieve the cumulative oil recovery of up to 36 - 37% (resulting from gravity drainage and low interfacial tension) and up to 52 - 53% (with additional well-to-well spacing adjustment) more than simple water injection methods, underscoring the potential of this innovative technology in the real-world applications. Discussions on the future research directions are also included.
subject: Original Research Papers
  topic: section
subject: surfactant imbibition
  topic: keywords
subject: horizontal wells
  topic: keywords
subject: DNAPL
  topic: keywords
subject: in-situ remediaiton
  topic: keywords
subject: simulation
  topic: keywords
subject: surfactant imbibition
  topic: subjects
subject: horizontal wells
  topic: subjects
subject: DNAPL
  topic: subjects
subject: in-situ remediation
  topic: subjects
subject: flow in porous media
  topic: disciplines
subject: enhanced oil recovery
  topic: disciplines
subject: in-situ remedaiiton
  topic: disciplines
creator: Betty Cepeda-Salgado
creator: Gyu Lee
creator: Frank Tsai
creator: Seung Kam
rights: Copyright (c) 2026 Betty Cepeda-Salgado, Gyu Lee, Frank Tsai, Seung Kam https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-18
publication date: 2026-03-18
doi: 10.69631/241eb066

Unraveling salt precipitation mechanisms in geological CO2 storage: insights into dominant driving forces
type: item
description: Salt precipitation during geological CO2 storage (GCS) in saline aquifers can impair injectivity of large-scale projects. Laboratory studies indicate that this process is sensitive to capillary force, but scaling-up to field conditions is complex due to radial flow, gravity, and geological heterogeneity. This study utilizes 3D simulations to characterize salt precipitation mechanisms in realistic storage formations. Non-dimensional groupings, e.g. capillary number (Ca) and gravity number (Gr), are used to characterize the relative impact of capillary, viscous, and gravity forces on salt formation. Three distinct regimes emerge for homogeneous systems: (1) capillary-dominated regime for Ca < 1 with severe salt precipitation at the inlet; (2) viscous-dominated regine for Ca > 1 and Gr < 0.2 with low-to-moderate salt precipitation; and (3) gravity-dominated regime for Ca > 1 and Gr > 0.2 where gravity override leads to localization of salt deposits inside the reservoir. The Ca-Gr regimes generally hold for heterogeneous systems. Additional complexity includes increased salt in high-permeability zones due to lower Ca enhanced by capillary-induced imbibition of saline water from low permeability zones. Tortuous flow paths can also emerge that further exacerbate pressure build-up. These novel characterization can assist practitioners in designing mitigation strategies for injectivity impairment in large-scale operations.
subject: Original Research Papers
  topic: section
subject: Salt Precipitation
  topic: keywords
subject: Injectivity Impairment
  topic: keywords
subject: Capillary number
  topic: keywords
subject: Reservoir Heterogeneity
  topic: keywords
subject: Geological CO2 Storage
  topic: keywords
creator: Nematollah Zamani
creator: David Landa Marban
creator: Tor Harald Sandve
creator: Sarah Eileen Gasda
rights: Copyright (c) 2026 Nematollah Zamani, David Landa Marban, Tor Harald Sandve, Sarah Eileen Gasda https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-18
publication date: 2026-03-18
doi: 10.69631/zhna7r09

Pareto-based Multi-Objective Optimization of Underground Hydrogen Storage Operational Condition Utilizing Machine Learning-Based Surrogate Models
type: item
description: Underground hydrogen storage (UHS) in depleted hydrocarbon reservoirs, saline aquifers, and salt caverns has emerged as a promising solution for large-scale, long-term energy storage to support the energy transition. However, evaluating the feasibility of UHS requires addressing substantial uncertainties in geological, operational, and economic factors. Physics-based reservoir simulations, though accurate, are computationally expensive and impractical for exploring a wide range of operational scenarios and optimization problems. To overcome these challenges, this study employs an intelligent surrogate model constructed using data collected from Latin hypercube sampling and validated against detailed reservoir simulations to perform sensitivity analysis and formulate a multi-objective optimization framework. The primary objectives are to maximize hydrogen recovery factor and hydrogen purity in the produced gas, while also incorporating net present value (NPV) as an economic criterion. Multi-objective optimization is performed using the Non-dominated Sorting Genetic Algorithm II (NSGA-II), which enables efficient exploration of the trade-offs between conflicting objectives and provides a Pareto front of non-dominated solutions. This Pareto-based decision-making framework offers valuable insights into selecting optimal storage strategies under uncertainty. The results highlight two key findings. First, pure nitrogen (N2) as a cushion gas yields higher hydrogen mole fraction in production stream compared to CO2 and CH4, because of favorable segregation and flow behavior in porous media. Second, the optimal cushion gas mixture was determined to be 45% N2, 40% CO2, and 15% CH4, balancing recovery and purity objectives. Nevertheless, from both economic and operational efficiency perspectives, a CO₂ fraction of 100% emerges as the optimal solution. Additionally, injecting from bottom perforations and producing from top perforations was identified as the most effective confi…
subject: Special Issue Submission - Iran2024
  topic: section
subject: Underground Hydrogen Storage
  topic: keywords
subject: Depleted Gas Reservoirs
  topic: keywords
subject: Machine Learning
  topic: keywords
subject: Surrogate models
  topic: keywords
subject: Genetic algorithms
  topic: keywords
subject: Energy Transition
  topic: keywords
subject: Underground Hydrogen Storage
  topic: subjects
subject: Machine Learning
  topic: subjects
subject: Depleted Gas Reservoirs
  topic: subjects
subject: Multi-Objective Optimization
  topic: subjects
subject: Cushion Gas
  topic: subjects
subject: Geoenergy Engineering
  topic: disciplines
creator: Hossein Kheirollahi
creator: Shahab Ayatollahi
creator: Hassan Mahani
rights: Copyright (c) 2026 Hossein Kheirollahi, Shahab Ayatollahi, Hassan Mahani https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-01
publication date: 2026-03-01
starting page: IPJ010326
ending page: 6
doi: 10.69631/m0g16y60

Highlights from the new issue of InterPore Journal
type: item
description: As we open this new issue of the InterPore Journal, I am encouraged by how rapidly our field is developing and by the role the journal plays in this progress. Each contribution in this edition reflects not only technical advancement but also the expanding relevance of porous media research to a wide range of challenges. The work presented here demonstrates a shared commitment to curiosity, rigor, and collaboration, which are necessary to sustain our community and push the boundaries of what we can achieve together.
subject: Editorial
  topic: section
subject: Editorial
  topic: keywords
creator: Nima Shokri
rights: Copyright (c) 2025 Nima Shokri https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-12-01
publication date: 2025-12-01
starting page: IPJ011225
ending page: 1
doi: 10.69631/6nwzqr77

A Comparative SPH Analysis of Interaction Between Viscous Debris Flow and Two Types of Barriers with Basal Clearance
type: item
description: Debris flows are a type of geophysical flow and are known as one of the most destructive geohazards. To control and mitigate damage to residential areas and infrastructure downstream, rigid barriers are installed along the flow path to dissipate the momentum and energy of the flow. While basal clearance and curved barriers have shown potential for granular flows, their combined application to viscous debris flows has remained unstudied. Therefore, this study aims to incorporate the advantages of barriers with basal clearance and curved geometries for application to viscous debris flows. To this end, a series of seven numerical simulations were performed using the smoothed particle hydrodynamics (SPH) method in LS-DYNA[1] to compare the performance of vertical and curved barriers with a basal clearance in the second row. The results indicate that curved barriers enhance the interaction process through rebounding, reducing the maximum kinetic energy by 59% in the most effective scenario (CNB2) compared to free flow conditions. Consequently, a system featuring a curved barrier in the first row and a narrow basal clearance in the second row was found to be the most practical. Furthermore, a single curved barrier reduced momentum by 16% more than a single vertical wall, performing similarly to dual vertical systems. Finally, it is noted that where site constraints or costs preclude two-barrier systems, a well-designed single curved barrier offers a viable alternative. [1] https://lsdyna.ansys.com/
subject: Special Issue Submission - Iran2024
  topic: section
subject: Debris flow
  topic: keywords
subject: Geohazard mitigation
  topic: keywords
subject: Curved barrier
  topic: keywords
subject: Basal clearance
  topic: keywords
subject: Smoothed Particle Hydrodynamics (SPH)
  topic: keywords
subject: Debris flow
  topic: subjects
subject: Geohazard mitigation
  topic: subjects
subject: Geotechnical Engineering
  topic: disciplines
creator: Anahita Aghagoli
creator: Hamed Sadeghi
rights: Copyright (c) 2026 Anahita Aghagoli, Hamed Sadeghi https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-01
publication date: 2026-03-01
ending page: 010326-7
doi: 10.69631/qkep8004

Design of a Laboratory Setup to Determine the Dispersion Parameter in a Single-Well Chemical Tracer Test
type: item
description: The Single Well Chemical Tracer Test (SWCTT) is a precise in-situ method for estimating residual oil saturation (Sor) based on the chromatographic delay between a partitioning tracer (ester) and its hydrolyzed product (alcohol). In this study, laboratory flooding experiments were designed and conducted in a glass-bead packed system at a constant temperature of 149°F to determine the dispersion and kinetic parameters of Ethyl Acetate (EtAc). The experiments were performed under two salinity conditions (5,000 ppm and 50,000 ppm) and simulated using a 1D model in the UTCHEM simulator with a 1 mm cell size. The results demonstrated that salinity significantly influences reaction kinetics and transport properties. Increasing salinity from 5,000 to 50,000 ppm resulted in an increase in the partition coefficient (KD) from 3 to 5.1 and the hydrolysis rate (KH) from 0.1 to 0.5/day. Furthermore, the hydrodynamic dispersion coefficient was determined through matching to be 0.003 ft2/day for the test conducted at a salinity of 5000 ppm and 0.01 ft2/day for the test conducted at a salinity of 50,000 ppm, resulting in accurate agreement between the simulated and experimental data. The Sor values of 0.2414 and 0.2303 obtained from the combined analysis of the model and experimental data showed excellent agreement with the reference value determined by material balance. This consistency serves as robust validation of the laboratory setup and confirms that the developed methodology reliably quantifies both Sor and hydrodynamic dispersion. Crucially, these findings establish a solid foundation for upscaling the methodology to pilot tests and reservoir-scale applications.
subject: Special Issue Submission - Iran2024
  topic: section
subject: Single well chemical tracer test (SWCTT)
  topic: keywords
subject: Oil Saturation
  topic: keywords
subject: Hydrodynamic Dispersion Coefficient
  topic: keywords
subject: Partitioning Coefficient
  topic: keywords
subject: Chromatographic
  topic: keywords
subject: UTCHEM
  topic: keywords
creator: Mostafa Ramezani
creator: Rohaldin Miri
rights: Copyright (c) 2026 Mostafa Ramezani, Rohaldin Miri https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-01
publication date: 2026-03-01
starting page: IPJ010326
ending page: 4
doi: 10.69631/dev9c518

Multiphase Image Segmentation of Naturally Fractured Media: Benchmarking Deep Learning and Conventional Approaches
type: item
description: Characterizing subsurface reservoirs, more specifically naturally fractured subsurface reservoir rocks, is essential for the study of subsurface reservoir properties. Image segmentation is an important aspect of digital rock physics (DRP) workflows. Traditional image segmentation techniques are less accurate than deep learning-based segmentation algorithms. In this paper, we investigate the segmentation accuracy of a convolutional neural network U-net and compare it with traditional methods of Watershed and multi-Otsu thresholding for multiphase segmentation of grayscale images from a naturally fractured coal sample. The segmentation target involved multiphase classification of the matrix, fully-filled fractures with minerals, and open fractures. The results reveal that U-net outperformed the others with an Intersection over Union metric of 94.9%, a Dice metric of 97%, and a Recall metric of 97.5%. The results support the importance of multiphase, deep learning-based segmentation techniques to support DRP studies of naturally fractured rocks.
subject: Special Issue Submission - Iran2024
  topic: section
subject: Naturally Fractured Rocks
  topic: keywords
subject: Digital Rock Physics
  topic: keywords
subject: Segmentation
  topic: keywords
subject: U-Net
  topic: keywords
subject: Watershed
  topic: keywords
subject: Multi-Otsu
  topic: keywords
creator: Behrad Tabrizipour
creator: Saeid Sadeghnejad
creator: Mastaneh Hajipour
creator: Thorsten Schäfer
rights: Copyright (c) 2026 Behrad Tabrizipour, Saeid Sadeghnejad, Mastaneh Hajipour, Thorsten Schäfer https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-01
publication date: 2026-03-01
starting page: IPJ010326
ending page: 5
doi: 10.69631/j6zv1891

Modelling copper mobilisation from low-grade ore controlled by mineral surface reactions
type: item
description: Many metals of economic value are commonly recovered from low grade ores through heap leaching at mine sites. The recovery of metals in heaps is often limited for multiple reasons including mineral overgrowth of targeted minerals resulting in suboptimal mineral dissolution rates. Chalcopyrite is a major copper iron sulphide mineral found in copper sulfide ores and its dissolution is controlled by multiple dissolution reaction mechanisms and potential mineral overgrowth. Traditional continuum scale reactive-transport models including those applied to simulate heap leaching do not account for local, mineral surface-controlled processes which may have implications for large-systems such as copper recovery in heaps. A new finite element, multicomponent reactive-transport model was developed and enabled the simulation of concurrent reactions at the discrete chalcopyrite surface and its surrounding gangue mineral and the fluid phase. The model predicts the rate of chalcopyrite dissolution controlled by different reaction mechanisms (proton- and ferric-iron promoted) as well as the rate of secondary mineral formation (jarosite) at the mineral surface. The latter may increasingly inhibit copper mobilisation over time through surface passivation. The onset and growth of the surface passivation layer is largely controlled by the surrounding fluid composition and ion transport through the gangue mineral. The reduction in porosity of the surface passivation layer leads to the breakdown of ionic transport to the chalcopyrite surface and associated copper mobilisation. Maximum copper recovery is achieved by high ferric iron concentrations and low pH and by preventing the formation of a passivation layer by keeping the sulphate concentration low.
subject: Original Research Papers
  topic: section
subject: Reactive Transport Modeling
  topic: keywords
subject: pore-scale
  topic: keywords
subject: reactive transport modelling
  topic: subjects
subject: mineral surface reactions
  topic: subjects
subject: Cu leaching
  topic: subjects
subject: geochemistry
  topic: disciplines
subject: hydrometallurgy
  topic: disciplines
creator: Ralf R. Haese
creator: Apoorv Jyoti
rights: Copyright (c) 2026 Prof. Ralf Haese, Dr. Apoorv Jyoti https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-18
publication date: 2026-03-18
doi: 10.69631/q237yd92

Highlights from the new issue of InterPore Journal
type: item
subject: Editorial
  topic: section
subject: Editorial
  topic: keywords
creator: Nima Shokri
rights: Copyright (c) 2025 Nima Shokri https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-08-25
publication date: 2025-08-25
starting page: IPJ250825
ending page: 1
doi: 10.69631/ipj.v2i3nr99

Exploring the Inter-Relationship between Climatic Conditions and Heat-Water Transfer in Deformable Cracked Surface Strata
type: item
description: Climate-induced soil desiccation and cracking present escalating challenges to infrastructure stability, necessitating comprehensive investigations of soil-climate interactions. Hence, this study presents a thermo-hydro-mechanical (THM) modeling framework to investigate the response of desiccation-cracked soil to climatic variations. The model was constructed using field observations and statistical analysis and incorporated three years of meteorological data from a susceptible area in Iran. The analysis explored heat conduction, water advection, mechanical deformation, and crack evolution under prolonged environmental conditions. Results revealed that heat flux was highly dependent on temperature gradients and crack geometry. Heat migrates upward from the subsurface when air temperature is lower than soil temperature or during rainfall, conditions typically associated with wet or cold seasons, whereas during warmer periods heat is transferred into the soil. Water flux exhibited strong seasonal variability, with rainfall-driven infiltration through the surface and enhanced evaporation near cracks during dry periods. During rainfall, water flow velocity in soil pores was at least 50 times higher than during warm weather. Moreover, a direct relationship between volumetric strain, temperature, and relative humidity was identified. Crack deformation intensified, and long-term climatic exposure resulted in cumulative shrinkage, increasing on average from 12.5% in the first year to 17.5% in the third year.
subject: Special Issue Submission - Iran2024
  topic: section
subject: Arid climate
  topic: keywords
subject: Cracked soil
  topic: keywords
subject: Heat conduction
  topic: keywords
subject: Water advection
  topic: keywords
subject: Deformation
  topic: keywords
creator: Milad Jabbarzadeh
creator: Hamed Sadeghi
rights: Copyright (c) 2026 Milad Jabbarzadeh, Hamed Sadeghi https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-01
publication date: 2026-03-01
starting page: IPJ010326
ending page: 2
doi: 10.69631/rk5nyd57

MAGNET: Medial Axis Guided Network Extraction Tool
type: item
description: Pore network models are useful for studying transport in porous materials in a computationally efficient way. Extraction of networks from volumetric images has evolved over the years, starting with medial axis-based approaches to more recent watershed segmentation. This paper reconsiders the classic medial axis method, which offers several advantages such as speed and topological correctness, and develops a modernized, updated, and improved version. The new method is named Medial Axis Guided Network Extraction Tool (MAGNET). It works by analyzing the skeleton of a porous material to identify pore centers at junctions and endpoints. Additional pore bodies are found on long throats using two different approaches. This work includes an efficient tool for calculating the cross-sectional area of throats with irregular shape by using walkers with an infinite mean-free path to probe the geometry orthogonal to the medial axis at the point of the throat constriction. This extra step was critical for obtaining an equivalent diameter needed to calculate the permeability. Lastly, MAGNET was written with computational efficiency in mind. The skeletonization approach was itself 4.2X faster than the SNOW watershed segmentation for a 10003 image. Additionally, a parallelized skeletonization was applied by processing the image in blocks with sufficient overlap which resulted in a 5.5X speed-up compared to the serial approach. To validate the output, MAGNET was tested on a 4003 voxel image of a Berea sandstone, and the flow and capillary properties of the extracted network were compared to the results from SNOW and the lattice-Boltzmann method. Structural information such as pore and throat size distribution and mercury intrusion curves was compared, and noticeable similarity was achieved. Crucially, the permeability predicted by MAGNET was within 5% of the lattice-Boltzmann prediction on the same image.
subject: Original Research Papers
  topic: section
subject: Network extraction
  topic: keywords
subject: Pore network modelling
  topic: keywords
subject: Image processing
  topic: keywords
subject: Permeability
  topic: keywords
subject: Porous media
  topic: keywords
creator: Michael McKague
creator: Hamed Fathiannasab
creator: Mohammad Amin Sadeghi
creator: Jeff Gostick
rights: Copyright (c) 2025 Michael McKague, Hamed Fathiannasab, Mohammad Amin Sadeghi, Jeff Gostick https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-12-01
publication date: 2025-12-01
starting page: IPJ011225
ending page: 7
doi: 10.69631/g47x8w91

Heat-activated Epoxy Foam for Permeability Alteration in Fractured Geothermal Fields: Proof of Concept
type: item
description: Geothermal energy plays a growing role in the transition to renewable and carbon free energy sources. A challenge for many geothermal operations is how to enhance water-rock heat exchange either by creation of new fractures in a tight rock, or by blocking short-circuiting large conduits. Here we report a novel approach for blocking large aperture (cm scale) fractures using heat-activated epoxy resin foam. The foam is injected as discrete inert resin droplets that are transported to regions of the geothermal field and activated upon reaching sufficiently high temperatures, where they undergo foaming and curing, thereby locally reducing permeability. In contrast to alternative methods for reducing rock permeability such as silicate gels or heat responsive polymer microbeads that target fractures of small apertures (< 0.1 mm), the epoxy foam can reduce permeability of fractures with large apertures (about 10 mm) while attaining mechanical strength and thermal stability. Results from laboratory aluminum-glass fracture models provide insights by visualizing the transport phase of resin droplets and their subsequent temperature-induced foaming and curing transformations with associated flow pathway blocking. Modeling results for transport and foaming in a simple fracture considering rheological properties and foaming (volume expansion) behavior are compared with measurements of permeability changes. Challenges associated with upscaling to fracture networks and large transport distances of resin droplets are discussed.
subject: Original Research Papers
  topic: section
subject: geothermal energy
  topic: keywords
subject: epoxy foam
  topic: keywords
subject: permeability reduction
  topic: keywords
creator: Dani Or
creator: Rishi Parashar
creator: Ying Yang
creator: Manish Bishwokarma
creator: Satish Karra
creator: Yutong Cui
rights: Copyright (c) 2026 Dani Or, Rishi Parashar, Ying Yang, Manish Bishwokarma, Satish Karra, Yutong Cui https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-10
publication date: 2026-03-10
doi: 10.69631/3mjx5y19

Optimization of the Energy Extraction Process from Enhanced Geothermal Systems
type: item
description: Enhanced Geothermal Systems (EGS) have emerged as a viable solution for sustainable energy extraction from hot dry rock reservoirs. In this study, numerical simulations are employed to investigate the optimization of energy extraction from enhanced geothermal systems, focusing on the effects of fracture geometry, the number of fractures, and well spacing on overall system performance and profitability. By validating the model against previously established semi-analytical and analytical models, this research explores the importance of the mentioned parameters in determining cumulative energy production and economic outcomes across different project lifespans. The results demonstrate that increased well spacing enhances energy output by delaying thermal breakthrough and increasing the optimal injection rates. Although small fractures have limited thermal capacity, the primary influence of fracture geometry is to limit the applicable distance between the wells. Deeper wells require more pumping power, which reduces net energy gain and lowers the optimal injection rate. Multiple fractures along the well enhance energy production, particularly at higher injection rates and for smaller fractures. Economic analyses underscore the importance of operational design (including the well placement and injection rate) in maximizing profitability.
subject: Special Issue Submission - Iran2024
  topic: section
subject: Enhanced geothermal systems (EGS)
  topic: keywords
subject: Hydrothermal simulation
  topic: keywords
subject: Finite fracture
  topic: keywords
subject: Economic optimization
  topic: keywords
subject: Hot dry rock (HDR)
  topic: keywords
subject: Heat transfer in porous media
  topic: keywords
creator: Mohammad Norouzi Delaviz
creator: Mozhdeh Sajjadi
creator: Mohammad Emami Niri
rights: Copyright (c) 2026 Mohammad Norouzi Delaviz, Mozhdeh Sajjadi, Mohammad Emami Niri https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-01
publication date: 2026-03-01
starting page: IPJ010326
ending page: 3
doi: 10.69631/thzmfy53

Effect of fissure drainage on dispersive mixing in porous media gravity currents: an experimental investigation
type: item
description: The prospect of storing hydrogen in depleted natural gas reservoirs requires that the dispersive mixing between injected hydrogen and the cushion gas resident within the porous medium be well understood. Motivated by this observation, we report upon a series of similitude laboratory experiments employing salt and fresh water; experiments use optical techniques to characterize the dispersive mixing that results from porous media gravity current flow. Importantly, the gravity current is supposed to propagate over a `"`leaky" boundary, i.e. one containing an isolated fissure that allows drainage into an underlying (and dynamically-isolated) layer. Thus do we find that drainage and dispersion are connected: as more (relatively undiluted) fluid drains through the fissure, the remaining gravity current fluid is disproportionately comprised of fluid that has been significantly diluted through dispersive mixing. A novel metric for comparing the relative amount of dispersed fluid present is proposed. We then characterize the variation of this so-called dispersed buoyancy fraction with the gravity current density and source height, the dip angle and the fissure geometry. Thus do we find, for example, that dispersion increases sharply for down-dip flow whereas dispersion may plateau as the fissure area is increased. Select experimental data are contrasted against the predictions of a previously-published theoretical model. The implications of said comparison are then discussed e.g. in the context of model refinement.
subject: Original Research Papers
  topic: section
subject: Dispersion
  topic: keywords
subject: Gravity current in porous media
  topic: keywords
subject: Underground hydrogen storage
  topic: keywords
creator: Mansimran Singh
creator: Saeed Sheikhi
creator: K.S. Bharath
creator: Morris R. Flynn
rights: Copyright (c) 2026 Mansimran Singh, Saeed Sheikhi, K.S. Bharath, Morris R. Flynn https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-18
publication date: 2026-03-18
doi: 10.69631/b5x6xn52

Petrophysical Characterization of Indiana Limestone Using Medical Dual-Energy Computed Tomography Technique: Insights into Porosity, Bulk Density, and Effective Atomic Number
type: item
description: Computed Tomography (CT) enables non-destructive 3D reconstruction of pore structures and rock properties mapping. Typically, such images are obtained from micro-CT (μm-scale) or synchrotron imaging (nm-scale). Despite their accuracy, these high-resolution imaging methods are expensive, time-consuming, and limited in sample size, affecting representative volume analysis. This work investigates an alternative approach using medical-CT (100 μm resolution) applying dual-energy CT technique (DECT) to characterize petrophysical properties (total porosity, bulk density, and effective atomic number) of six Indiana limestone samples (3.81 cm diameter, 4.88 cm length). Samples with porosities between 17.5% and 19.1% were scanned using paired high-energy (130 kV) and low-energy (80 kV) protocols, allowing DECT analysis to generate detailed 2D and 3D property maps. Additionally, the single-energy CT (SECT) technique, enhanced by subtracting images of the rock sample in saturated and dry conditions, improved the estimation of effective porosity. DECT and SECT results, processed with Python scripts and Avizo 3D, demonstrated average differences of 3.34% for bulk density, 5.30% for effective porosity, and 4.65% for effective atomic number compared to basic petrophysics measurements. Although artifacts from low-energy scans presented limitations, their impact can be reduced by optimizing acquisition parameters, improving the experimental setup, and applying reconstruction filtering techniques. Overall, this study highlights medical-CT as a fast, cost-effective method for analyzing larger samples, providing a practical alternative to traditional high-resolution imaging to estimate key petrophysical properties for the identification and visualization of heterogeneity on carbonates rocks.
subject: Original Research Papers
  topic: section
subject: Digital Rock Physics
  topic: keywords
subject: Petrophysics
  topic: keywords
subject: Dual-Energy CT
  topic: keywords
subject: Medical-CT
  topic: keywords
creator: Walter Leonardo Flores Antelo
creator: Janeth Alina Vidal Vargas
creator: Rosangela Barros Zanoni Lopes Moreno
rights: Copyright (c) 2026 Walter Leonardo Flores Antelo, Janeth Alina Vidal Vargas, Rosangela Barros Zanoni Lopes Moreno https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2026-03-10
publication date: 2026-03-10
doi: 10.69631/s3tb8r43

Suitability Mapping for Subsurface Floodwater Storage Schemes
type: item
description: The hydrological conditions a changing climate imposes require updated measures to address extreme water events sustainably. These must evolve from isolated solutions, such as water retention basins, to holistic management strategies that consider at least two situations simultaneously, such as floods and droughts. One strategy for sustainable groundwater management is managed aquifer recharge (MAR), which copes with decreasing groundwater levels by a targeted recharge of aquifers. Using high water as a recharge water source (Flood-MAR) can additionally cover the aspect of flood protection, resulting in a multi-beneficial solution for the region. This study evaluated the site suitability for subsurface floodwater storage schemes, which must focus more on assessing the existing aquifer characteristics and the specific flood dynamics in nearby rivers compared to known MAR schemes. Potential sites for underground flood storage are characterized by decreasing groundwater levels and, thus, frequent water scarcity, combined with eminent risks of flooding by a nearby river. Suitability, on the other hand, must take the aquifer, surface, and water source characteristics into account. In this study, we present a workflow for generating suitability maps for implementing subsurface floodwater storage systems with a geographic information system-based multi-criteria decision analysis (MCDA). The workflow was intentionally and exclusively based on publicly available data, was implemented in Python, and provided as open-source software. The resulting suitability maps spatially depict the feasibility of underground flood storage, and thus form the basis for the implementation planning of such projects. The approach was demonstrated for the administrative district of Swabia, Germany, where approximately 35% of the area was identified as suitable at varying levels. A sensitivity analysis of the assigned weights was applied to show the high robustness of the underlying data. The resu…
subject: Original Research Papers
  topic: section
subject: Flood and drought management
  topic: keywords
subject: Managed aquifer recharge (MAR)
  topic: keywords
subject: Geographic information systems (GIS)
  topic: keywords
subject: GIS-based multi-criteria decision analysis (GIS-MCDA)
  topic: keywords
subject: Floodwater storage
  topic: keywords
subject: Landscape water balance
  topic: keywords
creator: Lea Augustin
creator: Thomas Baumann
rights: Copyright (c) 2024 Lea Augustin, Thomas Baumann https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2024-08-24
publication date: 2024-08-24
starting page: ipj240824
ending page: 7
doi: 10.69631/ipj.v1i2nr20

Static Elastic Properties of Rocks Obtained by X-ray Microtomography, Inverse Modeling and Surrogate Model
type: item
description: Digital rock analysis uses imaging techniques to obtain information concerning the internal structure of the rock. In the context of obtaining mechanical properties from digital images, much attention is given to determining the elastic parameters of rocks based on their mineralogical composition. While many applications for this type of simulation are available in the literature, the quantification of mineral property variations resulting from external influences has largely been overlooked. This work adopts an inverse modeling approach to estimate such elastic property variations in rocks. This methodology requires a predictive tool and an optimization algorithm to iteratively update the model’s parameters while reducing the discrepancies between modeled and measured data. In this study, a representative sandstone sample from the Botucatu Formation in Brazil was used to evaluate the methodology. Two case studies were considered: a synthetic example was designed for validation of the methodology while the second case applies real-world laboratory test data to estimate the properties that are most supported by the experimental evidence. A MATLAB code was built to integrate a finite element program and a genetic algorithm in a single framework. Additionally, an artificial neural network was used as a surrogate model to reduce the computational time of the numerical forward run. Overall, the results support the robustness of the approach and present a new alternative to obtain the mechanical properties of rock constituents at the pore scale.
subject: Original Research Papers
  topic: section
subject: Elastic properties
  topic: keywords
subject: X-ray microtomography
  topic: keywords
subject: Inverse modeling
  topic: keywords
subject: Surrogate model
  topic: keywords
subject: Pore-scale
  topic: keywords
creator: Ruan Gomes
creator: Sergio Fontoura
creator: Guilherme Righetto
creator: Luiza Fernandes
creator: Rafael Lopez
creator: Rafaella Sampaio
creator: Claudio Lima
creator: Marcel Naumann
creator: William Silva
rights: Copyright (c) 2025 Ruan Gomes, Sergio Fontoura, Guilherme Righetto, Luiza Fernandes, Rafael Lopez, Rafaella Sampaio, Claudio Lima, Marcel Naumann, William Silva https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-12-01
publication date: 2025-12-01
starting page: IPJ011225
ending page: 5
doi: 10.69631/v65nzv74

A Hydro-Geomechanical Porous-Media Model to Study Effects of Engineered Carbonate Precipitation in Faults
type: item
description: Hydro-geomechanical models are required to predict or understand the impact of subsurface engineering applications such as gas storage in geological formations. This study focuses on engineered carbonate precipitation through biomineralization in a fault zone of a cap-rock, aiming to reduce gas leakage from a reservoir. In addition to altering hydraulic properties such as porosity and permeability, precipitated carbonates also change the mechanical properties of the rock. We present a conceptual modeling approach implemented in the open-source simulator DuMux. After model verification, we applied the model to a CO2-storage scenario to investigate how biomineralization affects stress distribution within the rock, as well as how it may alter the risks of fault reactivation and induced seismic events. The generic study shows a tendency towards increased stiffness due to precipitated carbonate, which may cause shear failure events to occur earlier than in an untreated setup, while the magnitude of the seismicity is smaller.
subject: Original Research Papers
  topic: section
subject: Hydro-geomechanical coupling
  topic: keywords
subject: Biomineralization
  topic: keywords
subject: Fault reactivation
  topic: keywords
subject: Model verification
  topic: keywords
creator: Yue Wang
creator: Holger Class
rights: Copyright (c) 2025 Yue Wang, Holger Class https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-06-04
publication date: 2025-06-04
starting page: IPJ040625
ending page: 6
doi: 10.69631/ipj.v2i2nr41

On the Interpretation of Unsteady State Experiments in Heterogeneous Rock by Stochastic Methods
type: item
description: Rock heterogeneity has a significant effect on immiscible displacement. This is especially true when the mobility ratio of the two fluids is unfavorable, favoring unstable displacement. However, this is not taken into account in the numerical analysis of classical core flooding experiments to quantify two-phase flow properties using Special Core Analysis (SCAL). Our approach combines the modern interpretation of SCAL data with experimental data measured on rock samples for which the homogeneity assumption - a prerequisite for SCAL experiments - can no longer apply due to their size and heterogeneity. In contrast to other studies that take heterogeneities into account, we focus on simple-to-perform unsteady-state experiments. We analyze these experiments by numerical interpretation using homogeneous and heterogeneous simulation domains and by introducing porosity-based heterogeneity and permeability as well as capillary scaling. In the current study, we first question the applicability of standard relative permeability measurements to heterogeneous rocks and fluids with an excessively high mobility ratio, such as for CO2-brine displacement in heterogeneous rocks. However, we show that they describe two-phase flow very well when porosity-based heterogeneity is taken into account, which is equivalent to downscaling. The study thus shows a way to fall back on established standard measurements if it should be possible to account for subgrid heterogeneities in SCAL workflows. To this end, we propose an approach based on steady-state experiments and appropriate sample selection.
subject: Original Research Papers
  topic: section
subject: Special Core Analysis
  topic: keywords
subject: CO2 sequestration
  topic: keywords
subject: unsteady-state experiments
  topic: keywords
subject: relative pereability
  topic: keywords
subject: capillary heterogeneity
  topic: keywords
subject: upscaling
  topic: keywords
creator: Omidreza Amrollahinasab
creator: Boris Jammernegg
creator: Siroos Azizmohammadi
creator: Holger Ott
rights: Copyright (c) 2025 Omidreza Amrollahinasab, Boris Jammernegg, Siroos Azizmohammadi, Holger Ott https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-06-04
publication date: 2025-06-04
starting page: IPJ040625
ending page: 5
doi: 10.69631/ipj.v2i2nr44

Experimental Investigation of Solubility Trapping in 3D Printed Micromodels
type: item
description: Understanding interfacial mass transfer during dissolution of gas in a liquid is vital for optimizing large-scale carbon capture and storage operations. While the dissolution of CO2 bubbles in reservoir brine is a crucial mechanism towards safe CO2 storage, it is a process that occurs at the pore-scale and is not yet fully understood. Direct numerical simulation (DNS) models describing this type of dissolution exist and have been validated with semi-analytical models on simple cases like a rising bubble in a liquid column. However, DNS models have not been experimentally validated for more complicated scenarios such as dissolution of trapped CO2 bubbles in pore geometries where there are few experimental datasets. In this work, we present an experimental and numerical study of trapping and dissolution of CO2 bubbles in 3D printed micromodel geometries. We used 3D printing technology to generate three different geometries, a single cavity geometry, a triple cavity geometry, and a multiple channel geometry. To investigate the repeatability of the trapping and dissolution experimental results, each geometry was printed three times, and three identical experiments were performed for each geometry. The experiments were performed at a low capillary number (Ca = 3.33 x 10-6), representative of flow during CO2 storage applications. The DNS simulations were then performed and compared with the experimental results. Our results show experimental reproducibility and consistency in terms of CO2 trapping and the CO2 dissolution process. At such a low capillary number, our numerical simulator cannot model the process accurately due to parasitic currents and the strong time-step constraints associated with capillary waves. However, we show that, for the single and triple cavity geometry, the interfacial transfer and resulting bubble dissolution can be reproduced by a numerical strategy where the interfacial tension is divided by 100 to relax the capillary time-step constraints. Th…
subject: Original Research Papers
  topic: section
subject: Mass Transfer
  topic: keywords
subject: CO2 Dissolution
  topic: keywords
subject: CO2 Trapping
  topic: keywords
subject: 3D printing
  topic: keywords
subject: Direct Numerical Simulation
  topic: keywords
creator: Alexandros Patsoukis Dimou
creator: Mahdi Mansouri Boroujeni
creator: Sophie Roman
creator: Hannah P. Menke
creator: Julien Maes
rights: Copyright (c) 2025 Alexandros Patsoukis Dimou, Mahdi Mansouri Boroujeni, Sophie Roman, Hannah Menke, Julien Maes https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-06-04
publication date: 2025-06-04
starting page: IPJ040625
ending page: 4
doi: 10.69631/ipj.v2i2nr49

3D Printing Reactive Porous Media: Calcite Precipitation Kinetics on Surface Functionalized Polymer Films and 3D-printed Cores
type: item
description: Mineral precipitation reactions in porous media can change the porosity and permeability of the rock formations. Predicting the rate of reaction and impacts on formation properties is challenging due to a lack of understanding of mineral precipitation reaction kinetics and mechanisms in porous media. This is furthermore challenging due to the highly heterogeneous nature of natural porous media. Here, we aim to develop a novel experimental platform leveraging 3D printing to facilitate replicable mineral precipitation experiments in controlled, heterogenous porous media systems. This requires fundamental understanding of the kinetics of mineral precipitation on the polymer materials used to fabricate the 3D printed porous media. In this work, we manipulate (via sulfonation) material surfaces (high impact polystyrene, HIPS) to promote calcite precipitation from supersaturated solutions to inform the design of synthetic subsurface systems. Calcite precipitation on HIPS films of varied surface sulfonation is confirmed using X-ray diffraction (XRD) analysis and weight-based precipitation experiments where increased precipitation with increased surface functionalization and solution saturation index are observed. This approach is then applied to 3D-printed porous media to enhance understanding of geochemical reactions, specifically calcite precipitation. Three dimensional images of Bentheimer Sandstone are used as the basis for 3D-printed porous media samples. Two 3D-printed samples were functionalized with acid to activate the surface and promote mineral precipitation. Functionalized and unfunctionalized samples underwent calcite precipitation core flooding experiments with oversaturated calcite solutions for 96 hours. Three dimensional X-ray micro-CT imaging revealed calcite growth in functionalized samples, with a calcite volume fraction of approximately 2.6% and a substantial reduction in porosity. Unfunctionalized samples exhibited diminished calcite precipitation and…
subject: Original Research Papers
  topic: section
subject: Mineral precipitation
  topic: keywords
subject: Surface functionalization
  topic: keywords
subject: 3D printing
  topic: keywords
subject: Porosity-permeability evolution
  topic: keywords
subject: X-ray CT imaging
  topic: keywords
creator: Harrish Kumar Senthil Kumar
creator: Abdullah Al Nahian
creator: Nailah Braziel
creator: Zhoufan Shi
creator: Madelyn Torrance
creator: Vinita V. Shinde
creator: Lauren E. Beckingham
creator: Bryan S. Beckingham
rights: Copyright (c) 2025 Harrish Kumar Senthil Kumar, Abdullah Al Nahian, Nailah Braziel, Zhoufan Shi, Madelyn Torrance, Vinita V. Shinde, Lauren E. Beckingham, Bryan Beckingham https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-06-04
publication date: 2025-06-04
starting page: IPJ040625
ending page: 7
doi: 10.69631/ipj.v2i2nr71

Experimental Study of Microbial Hydrogen Consumption Rates by Oleidesulfovibrio Alaskensis in Porous Media
type: item
description: The recovery efficiency of short- and long-term cyclic operations of porous media underground hydrogen storage (UHS) is a key parameter for successful implementation, but anaerobic microbes autochthonous in the storage formation can consume hydrogen and adversely influence hydrogen recoverability and storage efficiency. Here we have experimentally measured hydrogen consumption rates by a model sulphate-reducing bacterium (Oleidesulfovibrio alaskensis G20) in drainage-storage cycles that mimic porous media UHS. Laboratory tests were performed in cylindrical sand pack columns as storage site analogues (inner diameter: 14.7 mm, length: 51.4 mm) with an average porosity of 28% at conditions of 37°C and 1.15 bara. The storage capacity (initial hydrogen saturation in place) of each sand pack was also analyzed and compared against sterilized benchmarks. We observed an exponential decay in microbial hydrogen consumption between storage cycles: 28 ± 12% hydrogen was lost during the first cycle (with a peak average rate of 1.26 ± 0.12µmol/hr/cm3), compared with 10 ± 5% in the second cycle and 7 ± 3% in the third cycle. The cumulative loss across the three cycles amounted to 15 ± 6%, even though nutrient and carbon source concentrations were adequate for full hydrogen consumption in each cycle. The reduced microbial activity after the first storage cycle was explained by the observed increase in brine pH from an initial 7.5 to 8.4 ± 0.2 at the end of the last storage cycle. We observed improvement in the average hydrogen in place saturations after the first non-sterile storage cycles. Our experimental data enhances the understanding of microbial hydrogen loss during UHS and its impact on recovery and storage efficiency.
subject: Original Research Papers
  topic: section
subject: Underground hydrogen storage
  topic: keywords
subject: pH
  topic: keywords
subject: Microbial hydrogen consumption
  topic: keywords
subject: Anaerobic sulphate reduction
  topic: keywords
subject: Storage capacity
  topic: keywords
subject: Cyclic drainage and storage
  topic: keywords
creator: Raymond Mushabe
creator: Na Liu
creator: Nicole Dopffel
creator: Geir Ersland
creator: Martin A. Fernø
rights: Copyright (c) 2025 Raymond Mushabe, Na Liu, Nicole Dopffel, Geir Ersland, Martin Fernø https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-06-04
publication date: 2025-06-04
starting page: IPJ040625
ending page: 3
doi: 10.69631/ipj.v2i2nr72

Regional-Scale Subsurface Modeling
type: item
description: The subsurface is a critical foundation for society, providing essential resources such as water, minerals, and energy, while also offering storage capacity for resources and waste. Current subsurface use is expanding rapidly, leading to aggregated effects at regional scales ranging from tens to hundreds of kilometers. Despite their growing importance, current subsurface models remain limited to field-scale analysis and are insufficient for understanding regional scale impacts. This commentary advocates for a shift in subsurface modeling research, highlighting the importance of developing regional scale models that can assess the cumulative and interconnected effects of multiple subsurface operations. We outline a roadmap for developing the mathematical and numerical foundations necessary for regional scale subsurface analysis, highlighting the need for multiscale modeling, hybrid physics- and data-driven modeling approaches, and efficient computational tools.
subject: Invited Commentaries
  topic: section
subject: Geological
  topic: keywords
subject: Subsurface
  topic: keywords
subject: Resources
  topic: keywords
creator: Inga Berre
creator: Paola F. Antonietti
creator: Rainer Helmig
creator: Marc Hesse
creator: Jan Martin Nordbotten
rights: Copyright (c) 2025 Inga Berre, Paola Antonietti, Rainer Helmig, Marc Hesse, Jan Martin Nordbotten https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-06-04
publication date: 2025-06-04
starting page: IPJ040625
ending page: 2
doi: 10.69631/ipj.v2i2nr79

Highlights From the New Issue of InterPore Journal
type: item
subject: Editorial
  topic: section
subject: Editorial
  topic: keywords
creator: Nima Shokri
rights: Copyright (c) 2025 Nima Shokri https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-06-04
publication date: 2025-06-04
starting page: IPJ040625
ending page: 1
doi: 10.69631/ipj.v2i2nr88

Understanding Fracture Aperture and Permeability Evolution due to Carbonate Mineralization Utilizing 3D Printing
type: item
description: Caprock formations are a crucial part of subsurface-engineered systems. Composed largely of shale, caprocks act as natural barriers that prevent the upward migration of fluids, thereby ensuring the containment of stored substances in subsurface formations. Fractures in these formations are potential leakage pathways for stored fluids. Mineral precipitation reactions in these fractures, particularly calcite, can significantly restrict the fluid permeability, reducing leakage potential. However, predictive capabilities of mineral precipitation in fractures and associated permeability evolution are limited due to a lack of fundamental understanding of such reactions in natural samples, complicated by mineral heterogeneity and the complexity of the fracture structure. In this study, 3D-printed fracture samples are used to understand the impact of carbonate mineralization on fracture aperture and permeability evolution. Samples were printed using a digital light processing (DLP) 3D printer and commercial liquid resin. Calcite precipitation was first tested on printed 2D films before conducting plug flow column experiments aimed to understand fracture permeability changes due to mineral precipitation. Contact angle measurement and Fourier transform infrared (FTIR) spectroscopy on printed 2D films show evidence of a substantial amount of surface energy for calcite nucleation and precipitation. Surface topography analysis of printed fractured surfaces reveals comparable values, highlighting the high replicability of the printed samples. During the column experiments, the permeability reduces exponentially due to a decrease in fracture aperture. Reductions in fracture aperture estimated from effluent concentration and 3D X-ray computed tomography (CT) show comparable results. Moreover, 3D X-ray CT images suggest the impact of local flow velocities on precipitation. The insights gained from this research contribute to a deeper understanding of the permeability evolution due t…
subject: Original Research Papers
  topic: section
subject: 3D printing
  topic: keywords
subject: Rock fracture
  topic: keywords
subject: X-ray computed tomography
  topic: keywords
subject: Carbonate mineralization
  topic: keywords
subject: Permeability evolution
  topic: keywords
creator: Md. Fahim Salek
creator: Harrish Kumar Senthil Kumar
creator: Abdullah Al Nahian
creator: Samsul Mahmood
creator: Bryan S. Beckingham
creator: Kyle Schulze
creator: Lauren E. Beckingham
rights: Copyright (c) 2025 Md. Fahim Salek, Harrish Kumar Senthil Kumar, Abdullah Al Nahian, Samsul Mahmood, Bryan S. Beckingham, Kyle Schulze, Lauren E. Beckingham https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-12-01
publication date: 2025-12-01
starting page: IPJ011225
ending page: 6
doi: 10.69631/dpd8ha46

A novel robust optimization framework based on surrogate modeling for underground hydrogen storage in depleted natural gas reservoirs
type: item
description: Underground hydrogen storage (UHS) plays a vital role in global net-zero energy systems, enabling the storage of excess renewable energy for future use. However, physical reservoir model-based optimization for UHS system design and operation is computationally expensive due to complex geological properties and well-operational controls. This study developed a novel, efficient framework for UHS stochastic optimization to address this challenge, integrating advanced compositional reservoir simulation, accurate surrogate modeling, and stochastic optimization techniques. First, a base reservoir simulation model was developed to capture compositional fluid flow, hydrogen methanation reactions, gravity segregation, hysteresis, and capillary effects. To rapidly evaluate various well controls and reservoir configurations, convolutional neural network (CNN)-bi-directional long short-term memory (BiLSTM)-Attention models were trained as surrogate models using a comprehensive dataset generated from reservoir simulations. The CNN transforms three-dimensional (3D) geological fields into one-dimensional (1D) vectors, effectively capturing spatial features. The BiLSTM network learns the temporal evolution of the input features over time by processing them in both forward and backward directions. Subsequently, the attention mechanism enhances prediction accuracy by identifying and emphasizing the most significant features at critical time steps. The well-trained surrogate models were seamlessly integrated into the stochastic optimization framework based on the genetic algorithm, aiming to maximize the net present value (NPV) from UHS projects. The results demonstrate that the surrogate model exhibits satisfactory performance in the context of prediction accuracy, computational efficiency, and scalability. Notably, the newly developed framework based on surrogate models achieves an approximate 4878-fold speedup compared to an approach relying solely on reservoir simulation, while ma…
subject: Invited Student Papers
  topic: section
subject: Underground hydrogen storage
  topic: keywords
subject: Robust optimization
  topic: keywords
subject: Surrogate modeling
  topic: keywords
subject: Deep learning
  topic: keywords
creator: Zhilei Han
creator: Zeeshan Tariq
creator: Bicheng Yan
rights: Copyright (c) 2025 Zhilei Han, Zeeshan Tariq, Bicheng Yan https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-08-25
publication date: 2025-08-25
starting page: IPJ250825
ending page: 6
doi: 10.69631/ipj.v2i3nr69

Rubble Characteristics Associated with Room Collapse at the Waste Isolation Pilot Plant: Impact of Salt Clast Shapes and Size Distributions on the Depositional Pore System
type: item
description: The fluid transport properties of rubble associated with disposal room collapse within layered salt deposits at the Waste Isolation Pilot plant (WIPP) in southeastern New Mexico, USA are of concern when evaluating repository performance. Although properties of crushed salt have been studied extensively, salt rubble characteristics are less well known, while also being more difficult to characterize and model given the broader associated range in clast sizes. We describe a methodology for assessing rubble porosity and three-dimensional internal structure through deposition simulations that use high-fidelity renditions of clast shapes and that consider a broad range in clast sizes. We demonstrate the application of this approach using clast size and shape data that we collected on run of mine (ROM) material from WIPP, which represents the best available proxy for rubble material. Clast sizes from sieve analysis span ~1–100 mm with a mass weighted mean value of 13 mm. Simulated rubble deposits have interclast porosity values ranging from 34.1–38.6 vol% (n = 6). The largest clasts in these simulations tend to be underlain by the largest pores. Thus, although these clasts act as barriers for flow, they also are associated with highly permeable but localized flow paths. Scenarios involving alternative input configurations reveal that porosity values are highly sensitive to clast size variability (~10 vol% greater for nearly uniform compared to highly variable clast size distributions), clast shape (~10–13 vol% greater when using ROM shapes compared to spheres), and the extent of induced rearrangement (~3 vol% lower with limited rearrangement and ~10 vol% lower for extensive rearrangement). Porosity also is strongly affected by frictional coefficient values (~10 vol% lower for μ = 0.05 compared to μ = 1.0), but less so for values within the range of uncertainty for salt clasts (~2 vol% lower for μ = 0.62 compared to μ = 1.0).
subject: Original Research Papers
  topic: section
subject: Nuclear waste storage
  topic: keywords
subject: Nuclear waste disposal
  topic: keywords
subject: Waste isolation pilot plant
  topic: keywords
subject: Granular media
  topic: keywords
subject: Porosity
  topic: keywords
subject: Grain shape
  topic: keywords
subject: Polydispersity
  topic: keywords
creator: Robert Lander
creator: Edward Matteo
creator: Linda Bonnell
creator: Thomas Dewers
creator: Melisa Mills
creator: James Guilkey
creator: Chven Mitchell
creator: John Stormont
rights: Copyright (c) 2025 Robert Lander, Edward Matteo, Linda Bonnell, Thomas Dewers, Melisa Mills, James Guilkey, Chven Mitchell, John Stormont https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-08-25
publication date: 2025-08-25
starting page: IPJ250825
ending page: 3
doi: 10.69631/ipj.v2i3nr45

Anchored Physics-Informed Neural Network for Two-Phase Flow Simulation in Heterogeneous Porous Media
type: item
description: In this study, we propose a tensorization-anchoring strategy based on adaptive architectures to accelerate the computation of Enriched Physics-Informed Neural Networks (EPINN) for two-phase flow simulations. Specifically, we design an adaptive tensorization mechanism for the adjacency matrix embedding, the activation function, and the skip-gated connection in EPINN, which collectively expand the neural network's (NN) parameter space for learning more generalized patterns. Moreover, we developed an anchoring strategy by establishing Anchors-EPINN (An-EPINN). By detaching tensorization parameters from the computational graph and anchoring weighted nodes to fixed positions, the NN can benefit from tensorized fusion effects while reducing high-dimensional matrix calculations during forward and backward propagation, thereby enhancing simulation efficiency. This approach reduces execution time by 31.47% in homogeneous cases and 27.91% in heterogeneous cases, while maintaining higher computational accuracy.
subject: Invited Student Papers
  topic: section
subject: Artificial intelligence for partial differential equations
  topic: keywords
subject: AI4PDE
  topic: keywords
subject: Simulation
  topic: keywords
subject: Two-phase flow
  topic: keywords
subject: Heterogeneous
  topic: keywords
subject: Tensorizing
  topic: keywords
subject: Adaptive architecture
  topic: keywords
subject: Neural networks
  topic: keywords
creator: Jingqi Lin
creator: Xia Yan
creator: Kai Zhang
creator: Zhao Zhang
creator: Jun Yao
rights: Copyright (c) 2025 Jingqi Lin, Xia Yan, Kai Zhang, Zhao Zhang, Jun Yao https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-08-25
publication date: 2025-08-25
starting page: IPJ250825
ending page: 5
doi: 10.69631/ipj.v2i3nr67

Highlights From the Fourth issue of InterPore Journal
type: item
subject: Editorial
  topic: section
subject: Editorial
  topic: keywords
creator: Nima Shokri
rights: Copyright (c) 2025 Nima Shokri https://creativecommons.org/licenses/by-nc-nd/4.0
date: 2025-02-26
publication date: 2025-02-26
starting page: IPJ260225
ending page: 1
doi: 10.69631/ipj.v2i1nr77