Pore-scale modeling of two-phase flow: A comparison of the generalized network model to direct numerical simulation.

Despite recent advances in pore-scale modeling of two-phase flow through porous media, the relative strengths and limitations of various modeling approaches have been largely unexplored. In this work, two-phase flow simulations from the generalized network model (GNM) [Phys. Rev. E 96, 013312 (2017)2470-004510.1103/PhysRevE.96.013312; Phys. Rev. E 97, 023308 (2018)2470-004510.1103/PhysRevE.97.023308] are compared with a recently developed lattice-Boltzmann model (LBM) [Adv. Water Resour. 116, 56 (2018)0309-170810.1016/j.advwatres.2018.03.014; J. Colloid Interface Sci. 576, 486 (2020)0021-979710.1016/j.jcis.2020.03.074] for drainage and waterflooding in two samples-a synthetic beadpack and a micro-CT imaged Bentheimer sandstone-under water-wet, mixed-wet, and oil-wet conditions. Macroscopic capillary pressure analysis reveals good agreement between the two models, and with experiments, at intermediate saturations but shows large discrepancy at the end-points. At a resolution of 10 grid blocks per average throat, the LBM is unable to capture the effect of layer flow which manifests as abnormally large initial water and residual oil saturations. Critically, pore-by-pore analysis shows that the absence of layer flow limits displacement to invasion-percolation in mixed-wet systems. The GNM is able to capture the effect of layers, and exhibits predictions closer to experimental observations in water and mixed-wet Bentheimer sandstones. Overall, a workflow for the comparison of pore-network models with direct numerical simulation of multiphase flow is presented. The GNM is shown to be an attractive option for cost and time-effective predictions of two-phase flow, and the importance of small-scale flow features in the accurate representation of pore-scale physics is highlighted.

New type of pore-snap-off and displacement correlations in imbibition.

HYPOTHESIS: Imbibition of a fluid into a porous material involves the invasion of a wetting fluid in the pore space through piston-like displacement, film and corner flow, snap-off and pore bypassing. These processes have been studied extensively in two-dimensional (2D) porous systems; however, their relevance to three-dimensional (3D) natural porous media is poorly understood. Here, we investigate these pore-scale processes in a natural rock sample using time-resolved 3D (i.e., four-dimensional or 4D) X-ray imaging. EXPERIMENTS: We performed a capillary-controlled drainage-imbibition experiment on an initially brine-saturated carbonate rock sample. The sample was imaged continuously during imbibition using 4D X-ray imaging to visualize and analyze fluid displacement and snap-off processes at the pore-scale. FINDINGS: We discover a new type of snap-off that occurs in pores, resulting in the entrapment of a small portion of the non-wetting phase in pore corners. This contrasts with previously-observed snap-off in throats which traps the non-wetting phase in pore centers. We relate the new type of pore-snap-off to the pinning of fluid-fluid interfaces at rough surfaces, creating contact angles close to 90°. Subsequently, we provide correlations for displacement events as a function of pore-throat geometry. Our findings indicate that having a small throat does not necessarily favor snap-off: the key criterion is the throat radius in relation to the pore radius involved in a displacement event, captured by the aspect ratio.

Dynamic fluid configurations in steady-state two-phase flow in Bentheimer sandstone.

Fast synchrotron tomography is used to study the impact of capillary number, Ca, on fluid configurations in steady-state two-phase flow in porous media. Brine and n-decane were co-injected at fixed fractional flow, f_{w}=0.5, in a cylindrical Bentheimer sandstone sample for a range of capillary numbers 2.1×10^{-7}≤Ca≤4.2×10^{-5}, while monitoring the pressure differential. As we have demonstrated in Gao et al. [Phys. Rev. Fluids 5, 013801 (2020)2469-990X10.1103/PhysRevFluids.5.013801], dependent on Ca, different flow regimes have been identified: at low Ca only fixed flow pathways exist, while after a certain threshold dynamic effects are observed resulting in intermittent fluctuations in fluid distribution which alter fluid connectivity. Additionally, the flow paths, for each capillary number, were imaged multiple times to quantify the less frequent changes in fluid occupancy, happening over timescales longer than the duration of our scans (40 s). In this paper we demonstrate how dynamic connectivity results from the interaction between oil ganglia populations. At low Ca connected pathways of ganglia are fixed with time-independent small, medium, and large ganglia populations. However, with an increase in Ca we see fluctuations in the size and numbers of the larger ganglia. With the onset of intermittency, fluctuations occur mainly in pores and throats of intermediate size. When Ca is further increased, we see rapid changes in occupancy in pores of all size. By combining observations on pressure fluctuations and flow regimes at various capillary numbers, we summarize a phase diagram over a range of capillary numbers for the wetting and nonwetting phases, Ca_{w} and Ca_{nw}, respectively, to quantify the degree of intermittent flow. These different regimes are controlled by a competition between viscous forces on the flowing fluids and the capillary forces acting in the complex pore space. Furthermore, we plot the phase diagrams of the transition from Darcy flow to intermittent flow over a range of Reynolds and Weber numbers for the wetting and nonwetting phases to evaluate the balance among capillary, viscous, and inertial forces, incorporating data from the literature. We demonstrate that pore geometry has a significant control on flow regime.

Pore-scale imaging and analysis of low salinity waterflooding in a heterogeneous carbonate rock at reservoir conditions.

X-ray micro-tomography combined with a high-pressure high-temperature flow apparatus and advanced image analysis techniques were used to image and study fluid distribution, wetting states and oil recovery during low salinity waterflooding (LSW) in a complex carbonate rock at subsurface conditions. The sample, aged with crude oil, was flooded with low salinity brine with a series of increasing flow rates, eventually recovering 85% of the oil initially in place in the resolved porosity. The pore and throat occupancy analysis revealed a change in fluid distribution in the pore space for different injection rates. Low salinity brine initially invaded large pores, consistent with displacement in an oil-wet rock. However, as more brine was injected, a redistribution of fluids was observed; smaller pores and throats were invaded by brine and the displaced oil moved into larger pore elements. Furthermore, in situ contact angles and curvatures of oil-brine interfaces were measured to characterize wettability changes within the pore space and calculate capillary pressure. Contact angles, mean curvatures and capillary pressures all showed a shift from weakly oil-wet towards a mixed-wet state as more pore volumes of low salinity brine were injected into the sample. Overall, this study establishes a methodology to characterize and quantify wettability changes at the pore scale which appears to be the dominant mechanism for oil recovery by LSW.

Pore-by-pore modeling, analysis, and prediction of two-phase flow in mixed-wet rocks.

A pore-network model is an upscaled representation of the pore space and fluid displacement, which is used to simulate two-phase flow through porous media. We use the results of pore-scale imaging experiments to calibrate and validate our simulations, and specifically to find the pore-scale distribution of wettability. We employ energy balance to estimate an average, thermodynamic, contact angle in the model, which is used as the initial estimate of contact angle. We then adjust the contact angle of each pore to match the observed fluid configurations in the experiment as a nonlinear inverse problem. The proposed algorithm is implemented on two sets of steady state micro-computed-tomography experiments for water-wet and mixed-wet Bentheimer sandstone. As a result of the optimization, the pore-by-pore error between the model and experiment is decreased to less than that observed between repeat experiments on the same rock sample. After calibration and matching, the model predictions for capillary pressure and relative permeability are in good agreement with the experiments. The proposed algorithm leads to a distribution of contact angle around the thermodynamic contact angle. We show that the contact angle is spatially correlated over around 4 pore lengths, while larger pores tend to be more oil-wet. Using randomly assigned distributions of contact angle in the model results in poor predictions of relative permeability and capillary pressure, particularly for the mixed-wet case.

Pore-scale mechanisms of CO2 storage in oilfields.

Rapid implementation of global scale carbon capture and storage is required to limit temperature rises to 1.5 °C this century. Depleted oilfields provide an immediate option for storage, since injection infrastructure is in place and there is an economic benefit from enhanced oil recovery. To design secure storage, we need to understand how the fluids are configured in the microscopic pore spaces of the reservoir rock. We use high-resolution X-ray imaging to study the flow of oil, water and CO2 in an oil-wet rock at subsurface conditions of high temperature and pressure. We show that contrary to conventional understanding, CO2 does not reside in the largest pores, which would facilitate its escape, but instead occupies smaller pores or is present in layers in the corners of the pore space. The CO2 flow is restricted by a factor of ten, compared to if it occupied the larger pores. This shows that CO2 injection in oilfields provides secure storage with limited recycling of gas; the injection of large amounts of water to capillary trap the CO2 is unnecessary.

The architectural design of smart ventilation and drainage systems in termite nests.

Termite nests have been widely studied as effective examples for ventilation and thermoregulation. However, the mechanisms by which these properties are controlled by the microstructure of the outer walls remain unclear. Here, we combine multiscale X-ray imaging with three-dimensional flow field simulations to investigate the impact of the architectural design of nest walls on CO2 exchange, heat transport and water drainage. We show that termites build outer walls that contain both small and percolating large pores at the microscale. The network of larger microscale pores enhances permeability by one to two orders of magnitude compared to the smaller pores alone, and it increases CO2 diffusivity up to eight times. In addition, the pore network offers enhanced thermal insulation and allows quick drainage of rainwater, thereby restoring the ventilation and providing structural stability to the wet nest.

Generalized network modeling of capillary-dominated two-phase flow.

We present a generalized network model for simulating capillary-dominated two-phase flow through porous media at the pore scale. Three-dimensional images of the pore space are discretized using a generalized network-described in a companion paper [A. Q. Raeini, B. Bijeljic, and M. J. Blunt, Phys. Rev. E 96, 013312 (2017)2470-004510.1103/PhysRevE.96.013312]-which comprises pores that are divided into smaller elements called half-throats and subsequently into corners. Half-throats define the connectivity of the network at the coarsest level, connecting each pore to half-throats of its neighboring pores from their narrower ends, while corners define the connectivity of pore crevices. The corners are discretized at different levels for accurate calculation of entry pressures, fluid volumes, and flow conductivities that are obtained using direct simulation of flow on the underlying image. This paper discusses the two-phase flow model that is used to compute the averaged flow properties of the generalized network, including relative permeability and capillary pressure. We validate the model using direct finite-volume two-phase flow simulations on synthetic geometries, and then present a comparison of the model predictions with a conventional pore-network model and experimental measurements of relative permeability in the literature.

A study to investigate viscous coupling effects on the hydraulic conductance of fluid layers in two-phase flow at the pore level.

This paper examines the role of momentum transfer across fluid-fluid interfaces in two-phase flow. A volume-of-fluid finite-volume numerical method is used to solve the Navier-Stokes equations for two-phase flow at the micro-scale. The model is applied to investigate viscous coupling effects as a function of the viscosity ratio, the wetting phase saturation and the wettability, for different fluid configurations in simple pore geometries. It is shown that viscous coupling effects can be significant for certain pore geometries such as oil layers sandwiched between water in the corner of mixed wettability capillaries. A simple parametric model is then presented to estimate general mobility terms as a function of geometric properties and viscosity ratio. Finally, the model is validated by comparison with the mobilities computed using direct numerical simulation.

Validation of model predictions of pore-scale fluid distributions during two-phase flow.

Pore-scale two-phase flow modeling is an important technology to study a rock's relative permeability behavior. To investigate if these models are predictive, the calculated pore-scale fluid distributions which determine the relative permeability need to be validated. In this work, we introduce a methodology to quantitatively compare models to experimental fluid distributions in flow experiments visualized with microcomputed tomography. First, we analyzed five repeated drainage-imbibition experiments on a single sample. In these experiments, the exact fluid distributions were not fully repeatable on a pore-by-pore basis, while the global properties of the fluid distribution were. Then two fractional flow experiments were used to validate a quasistatic pore network model. The model correctly predicted the fluid present in more than 75% of pores and throats in drainage and imbibition. To quantify what this means for the relevant global properties of the fluid distribution, we compare the main flow paths and the connectivity across the different pore sizes in the modeled and experimental fluid distributions. These essential topology characteristics matched well for drainage simulations, but not for imbibition. This suggests that the pore-filling rules in the network model we used need to be improved to make reliable predictions of imbibition. The presented analysis illustrates the potential of our methodology to systematically and robustly test two-phase flow models to aid in model development and calibration.

Generalized network modeling: Network extraction as a coarse-scale discretization of the void space of porous media.

A generalized network extraction workflow is developed for parameterizing three-dimensional (3D) images of porous media. The aim of this workflow is to reduce the uncertainties in conventional network modeling predictions introduced due to the oversimplification of complex pore geometries encountered in natural porous media. The generalized network serves as a coarse discretization of the surface generated from a medial-axis transformation of the 3D image. This discretization divides the void space into individual pores and then subdivides each pore into sub-elements called half-throat connections. Each half-throat connection is further segmented into corners by analyzing the medial axis curves of its axial plane. The parameters approximating each corner-corner angle, volume, and conductivity-are extracted at different discretization levels, corresponding to different wetting layer thickness and local capillary pressures during multiphase flow simulations. Conductivities are calculated using direct single-phase flow simulation so that the network can reproduce the single-phase flow permeability of the underlying image exactly. We first validate the algorithm by using it to discretize synthetic angular pore geometries and show that the network model reproduces the corner angles accurately. We then extract network models from micro-CT images of porous rocks and show that the network extraction preserves macroscopic properties, the permeability and formation factor, and the statistics of the micro-CT images.