Changes in pore geometry and relative permeability caused by carbonate precipitation in porous media

Impact of hydraulic tortuosity on microporous and nanoporous media flow

Using two-dimensional porous structures made up of homogeneously arranged solid obstacles, we examine the effects of rarefaction on the hydraulic tortuosity in the slip and early transition flow regimes via extended lattice Boltzmann method. We observed that modification in either the obstacle's arrangement or the porosity led to a power-law relation between the porosity-tortuosity. Along with this, we also found that in the slip-flow regime, the exponent of this relation contains the effect of finite Knudsen number (Kn). In addition, we observed that on properly scaling Kn with porosity and hydraulic tortuosity, a generalized correlation can be obtained for apparent permeability.

Pore Scale Numerical Modelling of Geological Carbon Storage Through Mineral Trapping Using True Pore Geometries

Abstract

Mineral trapping (MT)is the most secure method of sequestering carbon for geologically significant periods of time. The processes behind MT fundamentally occur at the pore scale, therefore understanding which factors control MT at this scale is crucial. We present a finite elements advection–diffusion–reaction numerical model which uses true pore geometry model domains generated from $$\upmu$$ μ CT imaging. Using this model, we investigate the impact of pore geometry features such as branching, tortuosity and throat radii on the distribution and occurrence of carbonate precipitation in different pore networks over 2000 year simulated periods. We find evidence that a greater tortuosity, greater degree of branching of a pore network and narrower pore throats are detrimental to MT and contribute to the risk of clogging and reduction of connected porosity. We suggest that a tortuosity of less than 2 is critical in promoting greater precipitation per unit volume and should be considered alongside porosity and permeability when assessing reservoirs for geological carbon storage (GCS). We also show that the dominant influence on precipitated mass is the Damköhler number, or reaction rate, rather than the availability of reactive minerals, suggesting that this should be the focus when engineering effective subsurface carbon storage reservoirs for long term security.

Article Highlights

Three-Dimensional Membrane Imaging with X-ray Ptychography: Determination of Membrane Transport Properties for Membrane Distillation

Abstract

Membrane distillation (MD) is a desalination technique that uses a membrane to thermally separate potable water from sea or brackish water. The mass transport processes through the membrane are commonly described by the dusty gas model. These processes are modeled assuming uniform, ideally cylindrical capillaries and are adjusted for the membrane geometry by including porosity and tortuosity. The tortuosity is usually set to 2 or is used as an adjusting parameter to fit theoretical models to experimentally measured data. In this work, ptychographic X-ray computed tomography is employed to map the three-dimensional (3D) structure of three commercial state-of-the-art PTFE membranes in MD. The porosity, tortuosity and permeability (viscous flow coefficient) of the samples are computed using the lattice Boltzmann method. The intrinsic permeability is compared to the dusty gas model and an apparent permeability is proposed which is corrected for Knudsen slip effects at the membrane structure.

Article Highlights

A numerical study on oxygen adsorption in porous media of coal rock based on fractal geometry

In order to explore the factors affecting coal spontaneous combustion, the fractal characteristics of coal samples are tested, and a pore-scale model for oxygen adsorption in coal porous media is developed based on self-similar fractal model. The liquid nitrogen adsorption experiments show that the coal samples indicate evident fractal scaling laws at both low-pressure and high-pressure sections, and the fractal dimensions, respectively, represent surface morphology and pore structure of coal rock. The pore-scale model has been validated by comparing with available experimental data and numerical simulation. The present numerical results indicate that the oxygen adsorption depends on both the pore structures and temperature of coal rock. The oxygen adsorption increases with increased porosity, fractal dimension and ratio of minimum to maximum pore sizes. The edge effect can be clearly seen near the cavity/pore, where the oxygen concentration is low. The correlation between the oxygen adsorption and temperature is found to obey Langmuir adsorption theory, and a new formula for oxygen adsorption and porosity is proposed. This study may help understanding the mechanisms of oxygen adsorption and accordingly provide guidelines to lower the risk of spontaneous combustion of coal.

A simple construction strategy for fabrication of sulfur-doped silicate materials from attapulgite

In this work, a simple, economic and effective one step process in which sodium sulfide is used as a modifier and a source of sulfate doping has been proposed to fabricate sulfur-doped silicate materials (SATP).

The Impact of Mineral Dissolution on Drainage Relative Permeability and Residual Trapping in Two Carbonate Rocks

Carbon dioxide injection into deep saline aquifers is governed by a number of physico-chemical processes including mineral dissolution and precipitation, multiphase fluid flow, and capillary trapping. These processes can be coupled; however, the impact of fluid–rock reaction on the multiphase flow properties is difficult to study and is not simply correlated with variations in porosity. We observed the impact of rock mineral dissolution on multiphase flow properties in two carbonate rocks with distinct pore structures. Observations of steady-state $$\hbox {N}_2$$N2–water relative permeability and residual trapping were obtained, along with mercury injection capillary pressure characteristics. These tests alternated with eight stages in which 0.5% of the mineral volume was uniformly dissolved into solution from the rock cores using an aqueous solution with a temperature-controlled acid. Variations in the multiphase flow properties did not relate simply to changes in porosity, but corresponded to the changes in the underlying pore structure. In the Ketton carbonate, dissolution resulted in an increase in the fraction of pore volume made up by the smallest pores and a decrease in the fraction made up by the largest pores. This resulted in an increase in the relative permeability to the nonwetting phase, a decrease in the relative permeability to the wetting phase, and a modest, but systematic decrease in residual trapping. In the Estaillades carbonate, dissolution resulted in an increase in the fraction of pore volume made up by pores in the central range of the initial pore size distribution, and a corresponding decrease in the fraction made up by both the smallest and largest pores. This resulted in a decrease in the relative permeability to both the wetting and nonwetting fluid phases and no discernible impact on the residual trapping. In summary, the impact of rock matrix dissolution will be strongly dependent on the impact of that dissolution on the underlying pore structure of the rock. However, if the variation in pore structure can be observed or estimated with modelling, then it should be possible to estimate the impacts on multiphase flow properties.

Lattice Boltzmann simulations of liquid CO2 displacing water in a 2D heterogeneous micromodel at reservoir pressure conditions

We employed the color-fluid lattice Boltzmann multiphase model to simulate liquid CO₂ displacing water documented in experiments in a 2D heterogeneous micromodel at reservoir pressure conditions. The main purpose is to investigate whether lattice Boltzmann simulation can reproduce the CO₂ invasion patterns observed in these experiments for a range of capillary numbers. Although the viscosity ratio used in the simulation matches the experimental conditions, the viscosity of the fluids in the simulation is higher than that of the actual fluids used in the experiments. Doing so is required to enhance numerical stability, and is a common strategy employed in the literature when using the lattice Boltzmann method to simulate CO₂ displacing water. The simulations reproduce qualitatively similar trends of changes in invasion patterns as the capillary number is increased. However, the development of secondary CO₂ pathways, a key feature of the invasion patterns in the simulations and experiments, is found to occur at a much higher capillary number in the simulations compared with the experiments. Additional numerical simulations were conducted to investigate the effect of the absolute value of viscosity on the invasion patterns while maintaining the viscosity ratio and capillary number fixed. These results indicate that the use of a high viscosity (which significantly reduces the inertial effect in the simulations) suppresses the development of secondary CO₂ pathways, leading to a different fluid distribution compared with corresponding experiments at the same capillary number. Therefore, inertial effects are not negligible in drainage process with liquid CO₂ and water despite the low Reynolds number based on the average velocity, as the local velocity can be much higher due to Haines jump events. These higher velocities, coupled with the low viscosity of CO₂, further amplifies the inertial effect. Therefore, we conclude that caution should be taken when using proxy fluids that only rely on the capillary number and viscosity ratio in both experiment and simulation.

Impact of the kinetic boundary condition on porous media flow in the lattice Boltzmann formulation

To emphasize the importance of the kinetic boundary condition for micro- to nanoscale flow, we present an ad hoc kinetic boundary condition suitable for torturous geological porous media. We found that the kinetic boundary condition is one of the essential features which should be supplemented to the standard lattice Boltzmann scheme in order to obtain accurate continuum observables. The claim is validated using a channel flow setup by showing the agreement of mass flux with analytical value. Further, using a homogeneous porous structure, the importance of the kinetic boundary condition is shown by comparing the permeability correction factor with the analytical value. Finally, the proposed alternate to the kinetic boundary condition is validated by showing its capability to capture the basic feature of the kinetic boundary condition.

Characterization of reactive flow-induced evolution of carbonate rocks using digital core analysis - part 2: Calculation of the evolution of percolation and transport properties

Percolation of reactive fluids in carbonate rocks affects the rock microstructure and hence changes the rock macroscopic properties. In Part 1 paper, we examined the voxel-wise evolution of microstructure of the rock in terms of mineral dissolution/detachment, mineral deposition, and unchanged regions. In the present work, we investigate the relationships between changes in two characteristic transport properties, i.e. permeability and electrical conductivity and two critical parameters of the pore phase, i.e. the fraction of the pore space connecting the inlet and outlet faces of the core sample and the critical pore-throat diameter. We calculate the aforementioned properties on the images of the sample, wherein a homogeneous modification of pore structure occurred in order to ensure the representativeness of the calculated transport properties at the core scale. From images, the evolution of pore connectivity and the potential role of micropores on the connectivity are quantified. It is found that the changing permeability and electrical conductivity distributions along the core length are generally in good agreement with the longitudinal evolution of macro-connected macroporosity and the critical pore-throat diameter. We incorporate microporosity into critical length and permeability calculations and show how microporosity locally plays a role in permeability. It is shown that the Katz-Thompson model reasonably predicts the post-alteration permeability in terms of pre-alteration simulated parameters. This suggests that the evolution of permeability and electrical conductivity of the studied complex carbonate core are controlled by the changes in the macro-connected macroporosity as well as the smallest pore-throats between the connected macropores.

Simulation of dissolution in porous media in three dimensions with lattice Boltzmann, finite-volume, and surface-rescaling methods

We present a pore-scale dissolution model for the simulation of reactive transport in complex porous media such as those encountered in carbon-storage injection processes. We couple a lattice Boltzmann model for flow calculation with a finite-volume method for solving chemical transport equations, and allow the computational grid to change as mineral surfaces are dissolved according to first-order reaction kinetics. We appraise this scheme for use with high Péclet number flows in three-dimensional geometries and show how the popular first-order convection scheme is affected by severe numerical diffusion when grid Péclet numbers exceed unity, and confirm that this can be overcome relatively easily by using a second-order method in conjunction with a flux-limiter function. We then propose a surface rescaling method which uses parabolic elements to counteract errors in surface area exposed by the Cartesian grid and avoid the use of more complex embedded surface methods when surface reaction kinetics are incorporated. Finally, we compute dissolution in an image of a real porous limestone rock sample injected with HCl for different Péclet numbers and obtain dissolution patterns in concordance with theory and experimental observation. A low injection flow rate was shown to lead to erosion of the pore space concentrated at the face of the rock, whereas a high flow rate leads to wormhole formation.

Machine learning framework for analysis of transport through complex networks in porous, granular media: A focus on permeability

We present a data-driven framework to study the relationship between fluid flow at the macroscale and the internal pore structure, across the micro- and mesoscales, in porous, granular media. Sphere packings with varying particle size distribution and confining pressure are generated using the discrete element method. For each sample, a finite element analysis of the fluid flow is performed to compute the permeability. We construct a pore network and a particle contact network to quantify the connectivity of the pores and particles across the mesoscopic spatial scales. Machine learning techniques for feature selection are employed to identify sets of microstructural properties and multiscale complex network features that optimally characterize permeability. We find a linear correlation (in log-log scale) between permeability and the average closeness centrality of the weighted pore network. With the pore network links weighted by the local conductance, the average closeness centrality represents a multiscale measure of efficiency of flow through the pore network in terms of the mean geodesic distance (or shortest path) between all pore bodies in the pore network. Specifically, this study objectively quantifies a hypothesized link between high permeability and efficient shortest paths that thread through relatively large pore bodies connected to each other by high conductance pore throats, embodying connectivity and pore structure.

Lattice Boltzmann simulation of the gas-solid adsorption process in reconstructed random porous media

The gas-solid adsorption process in reconstructed random porous media is numerically studied with the lattice Boltzmann (LB) method at the pore scale with consideration of interparticle, interfacial, and intraparticle mass transfer performances. Adsorbent structures are reconstructed in two dimensions by employing the quartet structure generation set approach. To implement boundary conditions accurately, all the porous interfacial nodes are recognized and classified into 14 types using a proposed universal program called the boundary recognition and classification program. The multiple-relaxation-time LB model and single-relaxation-time LB model are adopted to simulate flow and mass transport, respectively. The interparticle, interfacial, and intraparticle mass transfer capacities are evaluated with the permeability factor and interparticle transfer coefficient, Langmuir adsorption kinetics, and the solid diffusion model, respectively. Adsorption processes are performed in two groups of adsorbent media with different porosities and particle sizes. External and internal mass transfer resistances govern the adsorption system. A large porosity leads to an early time for adsorption equilibrium because of the controlling factor of external resistance. External and internal resistances are dominant at small and large particle sizes, respectively. Particle size, under which the total resistance is minimum, ranges from 3 to 7 μm with the preset parameters. Pore-scale simulation clearly explains the effect of both external and internal mass transfer resistances. The present paper provides both theoretical and practical guidance for the design and optimization of adsorption systems.

Anomalous reactive transport in porous media: Experiments and modeling

We analyze dynamic behavior of chemically reactive species in a porous medium, subject to anomalous transport. In this context, we present transport experiments in a refraction-index-matched, three-dimensional, water-saturated porous medium. A pH indicator (Congo red) was used as either a conservative or a reactive tracer, depending on the tracer solution pH relative to that of the background solution. The porous medium consisted of an acrylic polymer material formed as spherical beads that have pH-buffering capacity. The magnitude of reaction during transport through the porous medium was related to the color change of the Congo red, via image analysis. Here, we focused on point injection of the tracer into a macroscopically uniform flow field containing water at a pH different from that of the injected tracer. The setup yielded measurements of the temporally evolving spatial (local-in-space) concentration field. Parallel experiments with the same tracer, but without reactions (no changes in pH), enabled identification of the transport itself to be anomalous (non-Fickian); this was quantified by a continuous time random walk (CTRW) formulation. A CTRW particle tracking model was then used to quantify the spatial and temporal migration of both the conservative and reactive tracer plumes. Model parameters related to the anomalous transport were determined from the conservative tracer experiments. An additional term accounting for chemical reaction was established solely from analysis of the reactant concentrations, and significantly, no other fitting parameters were required. The measurements and analysis emphasized the localized nature of reaction, caused by small-scale concentration fluctuations and preferential pathways. In addition, a threshold radius for pH-controlled reactive transport processes was defined under buffering conditions, which delineated the region in which reactions occurred rapidly.