Porosity and permeability alterations in processes of biomineralization in porous media - microfluidic investigations and their interpretation

Abstract

Motivation: Biomineralization refers to microbially induced processes resulting in mineral formations. In addition to complex biomineral structures frequently formed by marine organisms, like corals or mussels, microbial activities may also indirectly induce mineralization. A famous example is the formation of stromatolites, which result from biofilm activities that locally alter the chemical and physical properties of the environment in favor of carbonate precipitation. Recently, biomineralization gained attention as an engineering application. Especially with the background of global warming and the objective to reduce CO2 emissions, biomineralization offers an innovative and sustainable alternative to the usage of conventional Portland cement, whose production currently contributes significantly to global CO2 emissions. The most widely used method of biomineralization in engineering applications, is ureolytic calcium carbonate precipitation, which relies on the hydrolysis of urea and the subsequent precipitation of calcium carbonate. The hydrolysis of urea at moderate temperatures is relatively slow and therefore needs to be catalyzed by the enzyme urease to be practical for applications. Urease can be extracted from plants, for example from ground jack beans, and the process is consequently referred to as enzyme-induced calcium carbonate precipitation (ECIP). Another method is microbially induced calcium carbonate precipitation (MICP), which uses ureolytic bacteria that produce the enzyme in situ. EICP and MICP applications allow for producing various construction materials, stabilizing soils, or creating hydraulic barriers in the subsurface. The latter can be used, for example, to remediate leakages at the top layer of gas storage reservoirs, or to contain contaminant plumes in aquifers. Especially when remediating leakages in the subsurface, the most crucial parameter to be controlled is its intrinsic permeability. A valuable tool for predicting and planning field applications is the use of numerical simulation at the scale of representative elementary volumes (REV). For that, the considered domain is subdivided into several REV’s, which do not resolve the pore space in detail, but represent it by averaged parameters, such as the porosity and permeability. The porosity describes the ratio of the pore space to the considered bulk volume, and the permeability quantifies the ease of fluid flow through a porous medium. A change in porosity generally also affects permeability. Therefore, for REV-scale simulations, constitutive relationships are utilized to describe permeability as a function of porosity. There are several porosity-permeability relationships in the literature, such as the Kozeny-Carman relationship, Verma-Pruess, or simple power-law relationships. These constitutive relationships can describe individual states but usually do not include the underlying processes. Different boundary conditions during biomineralization may influence the course of porosity-permeability relationships. However, these relationships have not yet been adequately addressed. Pore-scale simulations are, in principle, very well suited to investigate pore space changes and their effects on permeability systematically. However, these simulations also rely on simplifications and assumptions. Therefore, it is essential to conduct experimental studies to investigate the complex processes during calcium carbonate precipitation in detail at the pore scale. Recent studies have shown that microfluidic methods are particularly suitable for this purpose. However, previous microfluidic studies have not explicitly addressed the impact of biomineralization on hydraulic effects. Therefore, this work aims to identify relevant phenomena at the pore scale to conclude on the REV-scale parameters, porosity and permeability, and their relationship. Contributions: This work comprises three publications. First, a suitable microfluidic setup and workflow were developed in Weinhardt et al. [2021a] to study pore space changes and the associated hydraulic effects reliably. This paper illustrated the benefits and insights of combining optical microscopy and micro X-ray computed tomography (micro XRCT) with hydraulic measurements in microfluidic chips. The elaborated workflow allowed for quantitative analysis of the evolution of calcium carbonate precipitates in terms of their size, shape, and spatial distribution. At the same time, their influence on differential pressure could be observed as a measure of flow resistance. Consequently, porosity and permeability changes could be determined. Along with this paper, we published two data sets [Weinhardt et al., 2021b, Vahid Dastjerdi et al., 2021] and set the basis for two other publications. In the second publication [von Wolff et al., 2021], the simulation results of a pore-scale numerical model, developed by Lars von Wolff, were compared to the experimental data of the first paper [Weinhardt et al., 2021b]. We observed a good agreement between the experimental data and the model results. The numerical studies complemented the experimental observations in allowing for accurate analysis of crystal growth as a function of local velocity profiles. In particular, we observed that crystal aggregates tend to grow toward the upstream side, where the supply of reaction products is higher than on the downstream side. Crystal growth during biomineralization under continuous inflow is thus strongly dependent on the locally varying velocities in a porous medium. In the third publication [Weinhardt et al., 2022a], we conducted further microfluidic experiments based on the experimental setup and workflow of the first contribution and published another data set [Weinhardt et al., 2022b]. We used microfluidic cells with a different, more realistic pore structure and investigated the influence of different injection strategies. We found that the development of preferential flow paths during EICP application may depend on the given boundary conditions. Constant inflow rates can lead to the development of preferential flow paths and keep them open. Gradually reduced inflow rates can mitigate this effect. In addition, we concluded that the coexistence of multiple calcium carbonate polymorphs and their transformations could influence the temporal evolution of porosity-permeability relationships.