Multiphase mass transport with partitioning and inter-phase transport in porous media

Effect of Pore Space Stagnant Zones on Interphase Mass Transfer in Porous Media, for Two-Phase Flow Conditions

Interphase mass transfer is an important solute transport process in two-phase flow in porous media. During two-phase flow, hydrodynamically stagnant and flowing zones are formed, with the stagnant ones being adjacent to the interfaces through which the interphase mass transfer happens. Due to the existence of these stagnant zones in the vicinity of the interface, the mass transfer coefficient decreases to a certain extent. There seems to be a phenomenological correlation between the mass transfer coefficient and the extent of the stagnant zone which, however, is not yet fully understood. In this study, the phase-field method-based continuous species transfer model is applied to simulate the interphase mass transfer of a dissolved species from the immobile, residual, non-aqueous phase liquid (NAPL) to the flowing aqueous phase. Both scenarios, this of a simple cavity and this of a porous medium, are investigated. The effects of flow rates on the mass transfer coefficient are significantly reduced when the stagnant zone and the diffusion length are larger. It is found that the stagnant zone saturation can be a proxy of the overall diffusion length of the terminal menisci in the porous medium system. The early-stage mass transfer coefficient continuously decreases due to the depletion of the solute in the small NAPL clusters that are in direct contact with the flowing water. The long-term mass transfer mainly happens on the interfaces associated with large NAPL clusters with larger diffusion lengths, and the mass transfer coefficient is mainly determined by the stagnant zone saturation.

GeoChemFoam: Direct Modelling of Multiphase Reactive Transport in Real Pore Geometries with Equilibrium Reactions

GeoChemFoam is an open-source OpenFOAM-based toolbox that includes a range of additional packages that solve various flow processes from multiphase transport with interface transfer, to single-phase flow in multiscale porous media, to reactive transport with mineral dissolution. In this paper, we present a novel multiphase reactive transport solver for simulations on complex pore geometries, including microfluidic devices and micro-CT images, and its implementation in GeoChemFoam. The geochemical model includes bulk and surface equilibrium reactions. Multiphase flow is solved using the Volume-Of-Fluid method, and the transport of species is solved using the continuous species transfer method. The reactive transport equations are solved using a sequential operator splitting method, with the transport step solved using GeoChemFoam, and the reaction step solved using Phreeqc, the US geological survey’s geochemical software. The model and its implementation are validated by comparison with analytical solutions in 1D and 2D geometries. We then simulate multiphase reactive transport in two test pore geometries: a 3D pore cavity and a 3D micro-CT image of Bentheimer sandstone. In each case, we show the pore-scale simulation results can be used to develop upscaled models that are significantly more accurate than standard macro-scale equilibrium models.

Effects of reservoir heterogeneity on scaling of effective mass transfer coefficient for solute transport

Modeling transport process at large scale requires proper scale-up of subsurface heterogeneity and an understanding of its interaction with the underlying transport mechanisms. A technique based on volume averaging is applied to quantitatively assess the scaling characteristics of effective mass transfer coefficient in heterogeneous reservoir models. The effective mass transfer coefficient represents the combined contribution from diffusion and dispersion to the transport of non-reactive solute particles within a fluid phase. Although treatment of transport problems with the volume averaging technique has been published in the past, application to geological systems exhibiting realistic spatial variability remains a challenge. Previously, the authors developed a new procedure where results from a fine-scale numerical flow simulation reflecting the full physics of the transport process albeit over a sub-volume of the reservoir are integrated with the volume averaging technique to provide effective description of transport properties. The procedure is extended such that spatial averaging is performed at the local-heterogeneity scale. In this paper, the transport of a passive (non-reactive) solute is simulated on multiple reservoir models exhibiting different patterns of heterogeneities, and the scaling behavior of effective mass transfer coefficient (Keff) is examined and compared. One such set of models exhibit power-law (fractal) characteristics, and the variability of dispersion and Keff with scale is in good agreement with analytical expressions described in the literature. This work offers an insight into the impacts of heterogeneity on the scaling of effective transport parameters. A key finding is that spatial heterogeneity models with similar univariate and bivariate statistics may exhibit different scaling characteristics because of the influence of higher order statistics. More mixing is observed in the channelized models with higher-order continuity. It reinforces the notion that the flow response is influenced by the higher-order statistical description of heterogeneity. An important implication is that when scaling-up transport response from lab-scale results to the field scale, it is necessary to account for the scale-up of heterogeneity. Since the characteristics of higher-order multivariate distributions and large-scale heterogeneity are typically not captured in small-scale experiments, a reservoir modeling framework that captures the uncertainty in heterogeneity description should be adopted.