Felix Weinhardt

Abstract

Alterations in the pore morphology of porous materials cause changes to the characteristic hydraulic properties, which are mostly non-linear and inherently difficult to predetermine. Assuming the alterations are known with sufficient accuracy, the relation between the altered pore structure, measured in terms of porosity, and intrinsic permeability may be determined by simulations with enormous computational effort. We focus on microfluidic experiments during the course of which the pore space becomes increasingly occupied with solid precipitate over elapsed process time. To analyze these domains, we present a novel geometry-informed drag formulation which allows for solving pseudo-3D Stokes equations for image-based input data of clogging porous media with accuracy and efficiency. In a pre-processing step, local pore space properties are analyzed and employed to spatially vary the magnitude of the drag term, which reflects the influence of neglected 3D effects. Calibration and validation is achieved through fully 3D Finite Difference Stokes simulations of different benchmark cases. With the proposed formulation we achieve the high accuracy of the pseudo-3D methods as far as permeability is concerned (<30% deviation), but also with respect to local velocities, for a microfluidic domain throughout the clogging process. Noteworthy, the computational cost is being reduced to less than 1%. Combining the efficiency of a Stokes 2D simulation and accuracy of a 3D model the presented approach is rendered an interesting option to investigate remaining open questions, for example on anisotropy of effective hydraulic parameters during the clogging process.

Abstract

Many subsurface engineering technologies or natural processes cause porous medium properties, such as porosity or permeability, to evolve in time. Studying and understanding such processes on the pore scale is strongly aided by visualizing the details of geometric and morphological changes in the pores. For realistic 3D porous media, X-Ray Computed Tomography (XRCT) is the method of choice for visualization. However, the necessary high spatial resolution requires either access to limited high-energy synchrotron facilities or data acquisition times which are considerably longer (e.g. hours) than the time scales of the processes causing the pore geometry change (e.g. minutes). Thus, so far, conventional benchtop XRCT technologies are often too slow to allow for studying dynamic processes. Interrupting experiments for performing XRCT scans is also in many instances no viable approach. We propose a novel workflow for investigating dynamic precipitation processes in porous media systems in 3D using a convent.

Abstract

Enzymatically induced calcium carbonate precipitation is a promising geotechnique with the potential, for example, to seal leakage pathways in the subsurface or to stabilize soils. Precipitation of calcium carbonate in a porous medium reduces the porosity and, consequently, the permeability. With pseudo-2D microfluidic experiments, including pressure monitoring and, for visualization, optical microscopy and X-ray computed tomography, pore-space alterations were reliably related to corresponding hydraulic responses. The study comprises six experiments with two different pore structures, a simple, quasi-1D structure, and a 2D structure. Using a continuous injection strategy with either constant or step-wise reduced flow rates, we identified key mechanisms that significantly influence the relationship between porosity and permeability. In the quasi-1D structure, the location of precipitates is more relevant to the hydraulic response (pressure gradients) than the overall porosity change. In the quasi-2D structure, this is different, because flow can bypass locally clogged regions, thus leading to steadier porosity–permeability relationships. Moreover, in quasi-2D systems, during continuous injection, preferential flow paths can evolve and remain open. Classical porosity–permeability power-law relationships with constant exponents cannot adequately describe this phenomenon. We furthermore observed coexistence and transformation of different polymorphs of calcium carbonate, namely amorphous calcium carbonate, vaterite, and calcite and discuss their influence on the observed development of preferential flow paths. This has so far not been accounted for in the state-of-the-art approaches for porosity–permeability relationships during calcium carbonate precipitation in porous media.

Abstract

Leakage mitigation methods are an important part of reservoir engineering and subsurface fluid storage, in particular. In the context of multi-phase systems of subsurface storage, e.g., subsurface CO2 storage, a reduction in the intrinsic permeability is not the only parameter to influence the potential flow or leakage; multi-phase flow parameters, such as relative permeability and capillary pressure, are key parameters that are likely to be influenced by pore-space reduction due to leakage mitigation methods, such as induced precipitation. In this study, we investigate the effects of enzymatically induced carbonate precipitation on capillary pressure&ndash;saturation relations as the first step in accounting for the effects of induced precipitation on multi-phase flow parameters. This is, to our knowledge, the first exploration of the effect of enzymatically induced carbonate precipitation on capillary pressure&ndash;saturation relations thus far. First, pore-scale resolved microfluidic experiments in 2D glass cells and 3D sintered glass-bead columns were conducted, and the change in the pore geometry was observed by light microscopy and micro X-ray computed tomography, respectively. Second, the effects of the geometric change on the capillary pressure&ndash;saturation curves were evaluated by numerical drainage experiments using pore-network modeling on the pore networks extracted from the observed geometries. Finally, parameters of both the Brooks&ndash;Corey and Van Genuchten relations were fitted to the capillary pressure&ndash;saturation curves determined by pore-network modeling and compared with the reduction in porosity as an average measure of the pore geometry&rsquo;s change due to induced precipitation. The capillary pressures increased with increasing precipitation and reduced porosity. For the 2D setups, the change in the parameters of the capillary pressure&ndash;saturation relation was parameterized. However, for more realistic initial geometries of the 3D samples, while the general patterns of increasing capillary pressure may be observed, such a parameterization was not possible using only porosity or porosity reduction, likely due to the much higher variability in the pore-scale distribution of the precipitates between the experiments. Likely, additional parameters other than porosity will need to be considered to accurately describe the effects of induced carbonate precipitation on the capillary pressure&ndash;saturation relation of porous media.

Abstract

We present version 3 of the open-source simulator for flow and transport processes in porous media DuMux. DuMux is based on the modular C++ framework Dune (Distributed and Unified Numerics Environment) and is developed as a research code with a focus on modularity and reusability. We describe recent efforts in improving the transparency and efficiency of the development process and community-building, as well as efforts towards quality assurance and reproducible research. In addition to a major redesign of many simulation components in order to facilitate setting up complex simulations in DuMux, version 3 introduces a more consistent abstraction of finite volume schemes. Finally, the new framework for multi-domain simulations is described, and three numerical examples demonstrate its flexibility.

Abstract

Enzymatically induced calcite precipitation (EICP) is an engineering technology that allows for targeted reduction of porosity in a porous medium by precipitation of calcium carbonates. This might be employed for reducing permeability in order to seal flow paths or for soil stabilization. This study investigates the growth of calcium-carbonate crystals in a micro-fluidic EICP setup and relies on experimental results of precipitation observed over time and under flow-through conditions in a setup of four pore bodies connected by pore throats. A phase-field approach to model the growth of crystal aggregates is presented, and the corresponding simulation results are compared to the available experimental observations. We discuss the model's capability to reproduce the direction and volume of crystal growth. The mechanisms that dominate crystal growth are complex depending on the local flow field as well as on concentrations of solutes. We have good agreement between experimental data and model results. In particular, we observe that crystal aggregates prefer to grow in upstream flow direction and toward the center of the flow channels, where the volume growth rate is also higher due to better supply.

Abstract

Abstract Enzymatically induced calcite precipitation (EICP) in porous media can be used as an engineering option to achieve precipitation in the pore space, for example, aiming at a targeted sealing of existing flow paths. This is accomplished through a porosity and consequent permeability alteration. A major source of uncertainty in modeling EICP is in the quantitative description of permeability alteration due to precipitation. This report presents methods for investigating experimentally the time-resolved effects of growing precipitates on porosity and permeability on the pore scale, in a poly-di-methyl-siloxane microfluidic flow cell. These methods include the design and production of the microfluidic cells, the preparation and usage of the chemical solutions, the injection strategy, and the monitoring of pressure drops for given fluxes for the determination of permeability. EICP imaging methods are explained, including optical microscopy and X-ray microcomputed tomography (XRCT), and the corresponding image processing and analysis. We present and discuss a new experimental procedure using a microfluidic cell, as well as the general perspectives for further experimental and numerical simulation studies on induced calcite precipitation. The results of this study show the enormous benefits and insights achieved by combining both light microscopy and XRCT with hydraulic measurements in microfluidic chips. This allows for a quantitative analysis of the evolution of precipitates with respect to their size and shape, while monitoring their influence on permeability. We consider this to be an improvement of the existing methods in the literature regarding the interpretation of recorded data (pressure, flux, and visualization) during pore morphology alteration.