-
Simulating carbon mineralization at pore scale in capillary networks of digital rock
Authors:
David A. Lazo Vasquez,
Jaione Tirapu Azpiroz,
Rodrigo Neumann Barros Ferreira,
Ronaldo Giro,
Manuela Fernandes Blanco Rodriguez,
Matheus Esteves Ferreira,
Mathias B. Steiner
Abstract:
Predicting the geometrical evolution of the pore space in geological formations due to fluid-solid interactions has applications in reservoir engineering, oil recovery, and geological storage of carbon dioxide. However, modeling frameworks that combine fluid flow with physical and chemical processes at a rock's pore scale are scarce. Here, we report a method for modeling a rock's pore space as a n…
▽ More
Predicting the geometrical evolution of the pore space in geological formations due to fluid-solid interactions has applications in reservoir engineering, oil recovery, and geological storage of carbon dioxide. However, modeling frameworks that combine fluid flow with physical and chemical processes at a rock's pore scale are scarce. Here, we report a method for modeling a rock's pore space as a network of connected capillaries and to simulate the capillary diameter modifications caused by reactive flow processes. Specifically, we model mineral erosion, deposition, dissolution, and precipitation processes by solving the transport equations iteratively, computing diameter changes within each capillary of the network simultaneously. Our automated modeling framework enables simulations on digital rock samples as large as (1.125mm)$^3$ with 125$\times 10^6$ voxels within seconds of CPU time per iteration. As an application of the computational method, we have simulated brine injection and calcium carbonate precipitation in sandstone. For quantitatively comparing simulation results obtained with models predicting either a constant or a flow-rate dependent precipitation, we track the time-dependent capillary diameter distribution as well as the permeability of the connected pore space. For validation and reuse, we have made the automated simulation workflow, the reactive flow model library, and the digital rock samples available in public repositories.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
pyMSER -- An open-source library for automatic equilibration detection in molecular simulations
Authors:
Felipe Lopes Oliveira,
Binquan Luan,
Pierre Mothé Esteves,
Mathias Steiner,
Rodrigo Neumann Barros Ferreira
Abstract:
Automated molecular simulations are used extensively for predicting material properties. Typically, these simulations exhibit two regimes: a dynamic equilibration part, followed by a steady state. For extracting observable properties, the simulations must first reach a steady state so that thermodynamic averages can be taken. However, as equilibration depends on simulation conditions, predicting t…
▽ More
Automated molecular simulations are used extensively for predicting material properties. Typically, these simulations exhibit two regimes: a dynamic equilibration part, followed by a steady state. For extracting observable properties, the simulations must first reach a steady state so that thermodynamic averages can be taken. However, as equilibration depends on simulation conditions, predicting the optimal number of simulation steps a priori is impossible. Here, we demonstrate the application of the Marginal Standard Error Rule (MSER) for automatically identifying the optimal truncation point in Grand Canonical Monte Carlo (GCMC) simulations. This novel automatic procedure determines the point in which steady state is reached, ensuring that figures-of-merits are extracted in an objective, accurate, and reproducible fashion. In the case of GCMC simulations of gas adsorption in metal-organic frameworks, we find that this methodology reduces the computational cost by up to 90%. As MSER statistics are independent of the simulation method that creates the data, this library is, in principle, applicable to any time series analysis in which equilibration truncation is required. The open-source Python implementation of our method, pyMSER, is publicly available for reuse and validation at https://github.com/IBM/pymser.
△ Less
Submitted 28 March, 2024;
originally announced March 2024.
-
Optimizing carbon dioxide trapping for geological storage
Authors:
Jaione Tirapu Azpiroz,
Ronaldo Giro,
Rodrigo Neumann Barros Ferreira,
Marcio Nogueira Pereira da Silva,
Manuela Fernandes Blanco Rodriguez,
Adolfo E. Correa Lopez,
David A. Lazo Vasquez,
Matheus Esteves Ferreira,
Mariana Del Grande,
Ademir Ferreira Da Silva,
Mathias B. Steiner
Abstract:
Carbon dioxide (CO2) trapping in capillary networks of reservoir rocks is a pathway to long-term geological storage. At pore scale, the CO2 trapping potential depends on injection pressure, temperature, and the rock's interaction with the surrounding fluids. Modeling this interaction requires adequate representations of both capillary volume and surface. For the lack of scalable representations, h…
▽ More
Carbon dioxide (CO2) trapping in capillary networks of reservoir rocks is a pathway to long-term geological storage. At pore scale, the CO2 trapping potential depends on injection pressure, temperature, and the rock's interaction with the surrounding fluids. Modeling this interaction requires adequate representations of both capillary volume and surface. For the lack of scalable representations, however, the prediction of a rock's CO2 storage potential has been challenging. Here, we report how to represent a rock's pore space by statistically sampled capillary networks (ssCN) that preserve morphological rock characteristics. We have used the ssCN method to simulate CO2 drainage within a representative sandstone sample at reservoir pressures and temperatures, exploring intermediate- and CO2-wet conditions. This wetting regime is often neglected, despite evidence of plausibility. By raising pressure and temperature we observe increasing CO2 penetration within the capillary network. For contact angles approaching 90 degrees, the CO2 saturation exhibits a pronounced maximum reaching 80 percent of the accessible pore volume. This is about twice as high as the saturation values reported previously. For enabling validation of our results and a broader application of our methodology, we have made available the rock tomography data, the digital rock computational workflows, and the ssCN models used in this study.
△ Less
Submitted 20 December, 2023;
originally announced December 2023.
-
Full scale, microscopically resolved tomographies of sandstone and carbonate rocks augmented by experimental porosity and permeability values
Authors:
Matheus Esteves Ferreira,
Mariana Del Grande,
Rodrigo Neumann Barros Ferreira,
Ademir Ferreira da Silva,
Márcio Nogueira Pereira da Silva,
Jaione Tirapu-Azpiroz,
Everton Lucas-Oliveira,
Arthur Gustavo de Araújo Ferreira,
Renato Soares,
Christian B. Eckardt,
Tito J Bonagamba,
Mathias Steiner
Abstract:
We report a dataset containing full-scale, 3D images of rock plugs augmented by petrophysical lab characterization data for application in digital rock and capillary network analysis. Specifically, we have acquired microscopically resolved tomography datasets of 18 cylindrical sandstone and carbonate rock samples having lengths of 25.4 mm and diameters of 9.5 mm, respectively. Based on the micro-t…
▽ More
We report a dataset containing full-scale, 3D images of rock plugs augmented by petrophysical lab characterization data for application in digital rock and capillary network analysis. Specifically, we have acquired microscopically resolved tomography datasets of 18 cylindrical sandstone and carbonate rock samples having lengths of 25.4 mm and diameters of 9.5 mm, respectively. Based on the micro-tomography data, we have computed porosity-values for each imaged rock sample. For validating the computed porosity values with a complementary lab method, we have measured porosity for each rock sample by using standard petrophysical characterization techniques. Overall, the tomography-based porosity values agree with the measurement results obtained from the lab, with values ranging from 8% to 30%. In addition, we provide for each rock sample the experimental permeabilities, with values ranging from 0.4 mD to above 5D. This dataset will be essential for establishing, benchmarking, and referencing the relation between porosity and permeability of reservoir rock at pore scale.
△ Less
Submitted 21 December, 2022; v1 submitted 19 December, 2022;
originally announced December 2022.
-
CRAFTED -- An exploratory database of simulated adsorption isotherms of metal-organic frameworks
Authors:
Felipe Lopes Oliveira,
Conor Cleeton,
Rodrigo Neumann Barros Ferreira,
Binquan Luan,
Amir H. Farmahini,
Lev Sarkisov,
Mathias Steiner
Abstract:
Grand Canonical Monte Carlo is an important method for performing molecular-level simulations and assisting the study and development of nanoporous materials for gas capture application. These simulations are based on the use of force fields and partial charges to model the interaction between the adsorbent molecules and the solid framework. The choice of the force field parameters and partial cha…
▽ More
Grand Canonical Monte Carlo is an important method for performing molecular-level simulations and assisting the study and development of nanoporous materials for gas capture application. These simulations are based on the use of force fields and partial charges to model the interaction between the adsorbent molecules and the solid framework. The choice of the force field parameters and partial charges can significantly impact the results obtained, however, there are very few databases available to support a comprehensive impact evaluation. Here, we present a database of simulations of CO2 and N2 adsorption isotherms on 726 metal-organic frameworks taken from the CoRE MOF 2014 database. We performed simulations with two force fields (UFF and DREIDING), four partial charge schemes (no charge, Qeq, EQeq, and DDEC), and three temperatures (273, 298, 323 K). The resulting isotherms compose the Charge-dependent, Reproducible, Accessible, Forcefield-dependent, and Temperature-dependent Exploratory Database (CRAFTED) of adsorption isotherms.
△ Less
Submitted 17 October, 2022;
originally announced October 2022.