Search

Hypothesis The change of wettability toward more water-wet by the injection of low salinity water can improve oil recovery from porous rocks, which is known as low salinity water flooding. To simulate this process at the pore-scale, we propose that the alteration in surface wettability mediated by thin water films which are below the resolution of simulation grid blocks has to be considered, as observed in experiments. This is modeled by a wettability alteration model based on rate-limited adsorption of ions onto the rock surface. Simulations The wettability alteration model is developed and incorporated into a lattice Boltzmann simulator which solves both the Navier-Stokes equation for oil/water two-phase flow and the advection-diffusion equation for ion transport. The model is validated against two experiments in the literature, then applied to 3D micro-CT images of a rock. Findings Our model correctly simulated the experimental observations caused by the slow wettability alteration driven by the development of water films. In the simulations on the 3D rock pore structure, a distinct difference in the mixing of high and low salinity water is observed between secondary and tertiary low salinity flooding, resulting in different oil recoveries.

Hypothesis: The wetting behaviour is a key property of a porous medium that controls hydraulic conductivity in multiphase flow. While many porous materials, such as hydrocarbon reservoir rocks, are initially wetted by the aqueous phase, surface active components within the non-wetting phase can alter the wetting state of the solid. Close to the saturation endpoints wetting phase fluid films of nanometre thickness impact the wetting alteration process. The properties of these films depend on the chemical characteristics of the system. Here we demonstrate that surface texture can be equally important and introduce a novel workflow to characterize the wetting state of a porous medium. Experiments: We investigated the formation of fluid films along a rock surface imaged with atomic force microscopy using zeta-potential measurements and a computational model for drainage. The results were compared to spontaneous imbibition test to link sub-pore-scale and core-scale wetting characteristics of the rock. Findings: The results show a dependency between surface coverage by oil, which controls the wetting alteration, and the macroscopic wetting response. The surface-area coverage is dependent on the capillary pressure applied during primary drainage. Close to the saturation endpoint, where the change in saturation was minor, the oil-solid contact changed more than 80%. (C) 2019 The Authors. Published by Elsevier Inc.

Within this work, we upscale the equations that describe the pore-scale behaviour of nonlinear porous elastic composites, using the asymptotic homogenization technique in order to derive the macroscale effective governing equations. A porous hyperelastic composite can be thought of as being comprised of a matrix interacting with a number of subphases and percolated by a fluid flowing in the pores (which is chosen to be Newtonian and incompressible here). A general nonlinear macroscale model is derived and is then specified for a particular choice of strain energy function, namely the de Saint-Venant function. This leads to a macroscale system of PDEs, which is of poroelastic type with additional terms and transformations to account for the nonlinear behaviour of the material. Our new porohyperelastic-type model describes the effective behaviour of nonlinear porous composites by prescribing the stress balance equations, the conservation of mass and Darcy’s law. The coefficients of these macroscale equations encode the detailed microstructure of the material and are to be found by solving pore-scale differential problems. The model reduces to the following limit cases of (a) linear poroelastic composites when the deformation gradient approaches the identity, (b) nonlinear composites when there are no pores and (c) nonlinear poroelasticity when only the matrix–fluid interaction is considered. This model is applicable when the interactions between various hyperelastic solid phases occur at the pore-scale, as in biological tissues such as artery walls, the myocardium, lungs and liver.

Core flooding experiments to determine multiphase flow in properties of rock such as relative permeability can show significant fluctuations in terms of pressure, saturation, and electrical conductivity. That is typically not considered in the Darcy scale interpretation but treated as noise. However, in recent years, flow regimes that exhibit spatio-temporal variations in pore scale occupancy related to fluid phase pressure changes have been identified. They are associated with topological changes in the fluid configurations caused by pore-scale instabilities such as snap-off. The common understanding of Darcy-scale flow regimes is that pore-scale phenomena and their signature should have averaged out at the scale of representative elementary volumes (REV) and above. In this work, it is demonstrated that pressure fluctuations observed in centimeter-scale experiments commonly considered Darcy-scale at fractional flow conditions, where wetting and non-wetting phases are co-injected into porous rock at small (−6) capillary numbers are ultimately caused by pore-scale processes, but there is also a Darcy-scale fractional flow theory aspect. We compare fluctuations in fractional flow experiments conducted on samples of few centimeters size with respective experiments andin-situmicro-CT imaging at pore-scale resolution using synchrotron-based X-ray computed micro-tomography. On that basis we can establish a systematic causality from pore to Darcy scale. At the pore scale, dynamic imaging allows to directly observe the associated breakup and coalescence processes of non-wetting phase clusters, which follow “trajectories” in a “phase diagram” defined by fractional flow and capillary number and can be used to categorize flow regimes. Connected pathway flow would be represented by a fixed point, whereas processes such as ganglion dynamics follow trajectories but are still overall capillary-dominated. That suggests that the origin of the pressure fluctuations observed in centimeter-sized fractional flow experiments are capillary effects. The energy scale of the pressure fluctuations corresponds to 105-106times the thermal energy scale. This means the fluctuations are non-thermal. At the centimeter scale, there are non-monotonic and even oscillatory solutions permissible by the fractional flow theory, which allow the fluctuations to be visible and—depending on exact conditions—significant at centimeter scale, within the viscous limit of classical (Darcy scale) fractional flow theory. That also means that the phenomenon involves both capillary aspects from the pore or cluster scale and viscous aspects of fractional flow and occurs right at the transition, where the physical description concept changes from pore to Darcy scale.

The main means of secondary oil recovery is water flooding, which has been widely used in various oilfields. Different flow rates have a great impact on the recovery ratio and the occurrence of remaining oil. Scholars have carried out extensive research on it, but mostly on the macro scale, and research on the three-dimensional micro scale is also limited by accuracy and a lack of accurate understanding. In this paper, micro-CT and core displacement experiments are used to intuitively show the occurrence state of remaining oil under different flow rates. Through a series of quantitative image processing methods and remaining oil classification methods, the occurrence characteristics of remaining oil under different flow rates are systematically evaluated and studied. The results show that: (1) As the displacement rate increases, the remaining oil saturation decreases (61%, 35%, 23%), but the remaining oil is more evenly distributed along the slice, (2) Two lower displacement speeds (0.003 mL/min, 0.03 mL/min) can reduce the volume of huge oil clusters under oil-saturated conditions, and the highest displacement speed (0.3 mL/min) can completely break up large oil clusters into small oil droplets. At the same time, the shape factor of the oil clusters also gradually increases, (3) The proportion of continuous remaining oil volume decreases, and the proportion of discontinuous remaining oil increases. Discontinuous remaining oil is the main production target of EOR, (4) After water flooding, the microscopic remaining oil is more inclined to the middle and corner parts of the larger pores.

Three-phase flow in porous media is encountered in many applications including subsurface carbon dioxide storage, enhanced oil recovery, groundwater remediation and the design of microfluidic devices. However, the pore-scale physics that controls three-phase flow under capillary dominated conditions is still not fully understood. Recent advances in three-dimensional pore-scale imaging have provided new insights into three-phase flow. Based on these findings, this paper describes the key pore-scale processes that control flow and trapping in a three-phase system, namely wettability order, spreading and wetting layers, and double/multiple displacement events. We show that in a porous medium containing water, oil and gas, the behaviour is controlled by wettability, which can either be water-wet, weakly oil-wet or strongly oil-wet, and by gas–oil miscibility. We provide evidence that, for the same wettability state, the three-phase pore-scale events are different under near-miscible conditions—where the gas–oil interfacial tension is ≤ 1 mN/m—compared to immiscible conditions. In a water-wet system, at immiscible conditions, water is the most-wetting phase residing in the corners of the pore space, gas is the most non-wetting phase occupying the centres, while oil is the intermediate-wet phase spreading in layers sandwiched between water and gas. This fluid configuration allows for double capillary trapping, which can result in more gas trapping than for two-phase flow. At near-miscible conditions, oil and gas appear to become neutrally wetting to each other, preventing oil from spreading in layers; instead, gas and oil compete to occupy the centre of the larger pores, while water remains connected in wetting layers in the corners. This allows for the rapid production of oil since it is no longer confined to movement in thin layers. In a weakly oil-wet system, at immiscible conditions, the wettability order is oil–water–gas, from most to least wetting, promoting capillary trapping of gas in the pore centres by oil and water during water-alternating-gas injection. This wettability order is altered under near-miscible conditions as gas becomes the intermediate-wet phase, spreading in layers between water in the centres and oil in the corners. This fluid configuration allows for a high oil recovery factor while restricting gas flow in the reservoir. Moreover, we show evidence of the predicted, but hitherto not reported, wettability order in strongly oil-wet systems at immiscible conditions, oil–gas–water, from most to least wetting. At these conditions, gas progresses through the pore space in disconnected clusters by double and multiple displacements; therefore, the injection of large amounts of water to disconnect the gas phase is unnecessary. We place the analysis in a practical context by discussing implications for carbon dioxide storage combined with enhanced oil recovery before suggesting topics for future work.

Geological storage of CO2 in deep saline aquifers mitigates atmospheric emissions. In situ storage is facilitated by several trapping mechanisms including residual trapping, which plays a major role in the containment of CO2. Understanding the underlying mechanisms of residual trapping is crucial for planning storage projects. Of special interest is the relationship between the initial and residual CO2 saturations-the so-called IR curve, needed for predictive macroscopic-scale simulations. This study aims to improve the understanding of residual trapping in sandstone from the Heletz site, where extensive field experiments have been performed, by using 3D-image analysis on core sample CT-data. This was done to gain knowledge on physical properties (such as radius, coordination number, aspect ratio, shape factor of pores, and pore connectivity) of importance to residual CO2 trapping. Pore-network flow modeling on a network representation, with the extracted pore-space topology, was employed to estimate the IR curve. The core sample exhibited pores with a large range of coordination numbers, a mean aspect ratio of 1.4, and shape factors mostly corresponding to triangular cross-sections. The estimated IR curve was monotonic, fitting an Aissaoui-type trapping model, displaying a lower sensitivity to the advancing contact angle than previously thought, and indicating a good ability to residually trap CO2. This study provides the first report of values for the three above mentioned properties for Heletz sandstone, and the first estimate of an IR curve for CO2/brine in Heletz sandstone based on pore-network modeling on a network with a topology retrieved from a core-sample CT-scan. (c) 2021 Society of Chemical Industry and John Wiley & Sons, Ltd.

We use synchrotron X-ray micro-tomography to investigate the displacement dynamics during three-phase—oil, water and gas—flow in a hydrophobic porous medium. We observe a distinct gas invasion pattern, where gas progresses through the pore space in the form of disconnected clusters mediated by double and multiple displacement events. Gas advances in a process we name three-phase Haines jumps, during which gas re-arranges its configuration in the pore space, retracting from some regions to enable the rapid filling of multiple pores. The gas retraction leads to a permanent disconnection of gas ganglia, which do not reconnect as gas injection proceeds. We observe, in situ , the direct displacement of oil and water by gas as well as gas–oil–water double displacement. The use of local in situ measurements and an energy balance approach to determine fluid–fluid contact angles alongside the quantification of capillary pressures and pore occupancy indicate that the wettability order is oil–gas–water from most to least wetting. Furthermore, quantifying the evolution of Minkowski functionals implied well-connected oil and water, while the gas connectivity decreased as gas was broken up into discrete clusters during injection. This work can be used to design CO 2 storage, improved oil recovery and microfluidic devices.

The diffusion–adsorption behavior of methane in coal is an important factor that both affecting the decay rate of gas production and the total gas production capacity. In this paper, we established a pore-scale Lattice Boltzmann (LB) model coupled with fluid flow, gas diffusion, and gas adsorption–desorption in the bi-dispersed porous media of coalbed methane. The Knudsen diffusion and dynamic adsorption–desorption of gas in clusters of coal particles were considered. Firstly, the model was verified by two classical cases. Then, three dimensionless numbers, Re, Pe, and Da, were adopted to discuss the impact of fluid velocity, gas diffusivity, and adsorption/desorption rate on the gas flow–diffusion–adsorption process. The effect of the gas adsorption layer in micropores on the diffusion–adsorption–desorption process was considered, and a Langmuir isotherm adsorption theory-based method was developed to obtain the dynamic diffusion coefficient, which can capture the intermediate process during adsorption/desorption reaches equilibrium. The pore-scale bi-disperse porous media of coal matrix was generated based on the RCP algorithm, and the characteristics of gas diffusion and adsorption in the coal matrix with different Pe, Da, and pore size distribution were discussed. The conclusions were as follows: (1) the influence of fluid velocity on the diffusion–adsorption process of coalbed methane at the pore-scale is very small and can be ignored; the magnitude of the gas diffusivity in macropores affects the spread range of the global gas diffusion and the process of adsorption and determines the position where adsorption takes place preferentially. (2) A larger Fickian diffusion coefficient or greater adsorption constant can effectively enhance the adsorption rate, and the trend of gas concentration- adsorption is closer to the Langmuir isotherm adsorption curve. (3) The gas diffusion–adsorption–desorption process is affected by the adsorption properties of coal: the greater the pL or Vm, the slower the global gas diffusivity decay. (4) The effect of the gas molecular adsorption layer has a great impact on the kinetic process of gas diffusion–adsorption–desorption. Coal is usually tight and has low permeability, so it is difficult to ensure that the gas diffusion and adsorption are sufficient, the direct use of a static isotherm adsorption equation may be incorrect.

Understanding transport phenomena and governing mechanisms of different physical and chemical processes in porous media has been a critical research area for decades. Correlating fluid flow behaviour at the micro-scale with macro-scale parameters, such as relative permeability and capillary pressure, is key to understanding the processes governing subsurface systems, and this in turn allows us to improve the accuracy of modelling and simulations of transport phenomena at a large scale. Over the last two decades, there have been significant developments in our understanding of pore-scale processes and modelling of complex underground systems. Microfluidic devices (micromodels) and imaging techniques, as facilitators to link experimental observations to simulation, have greatly contributed to these achievements. Although several reviews exist covering separately advances in one of these two areas, we present here a detailed review integrating recent advances and applications in both micromodels and imaging techniques. This includes a comprehensive analysis of critical aspects of fabrication techniques of micromodels, and the most recent advances such as embedding fibre optic sensors in micromodels for research applications. To complete the analysis of visualization techniques, we have thoroughly reviewed the most applicable imaging techniques in the area of geoscience and geo-energy. Moreover, the integration of microfluidic devices and imaging techniques was highlighted as appropriate. In this review, we focus particularly on four prominent yet very wide application areas, namely “fluid flow in porous media”, “flow in heterogeneous rocks and fractures”, “reactive transport, solute and colloid transport”, and finally “porous media characterization”. In summary, this review provides an in-depth analysis of micromodels and imaging techniques that can help to guide future research in the in-situ visualization of fluid flow in porous media.

We identify a distinct two-phase flow invasion pattern in a mixed-wet porous medium. Time-resolved high-resolution synchrotron X-ray imaging is used to study the invasion of water through a small rock sample filled with oil, characterized by a wide non-uniform distribution of local contact angles both above and below 90 ° . The water advances in a connected front, but throats are not invaded in decreasing order of size, as predicted by invasion percolation theory for uniformly hydrophobic systems. Instead, we observe pinning of the three-phase contact between the fluids and the solid, manifested as contact angle hysteresis, which prevents snap-off and interface retraction. In the absence of viscous dissipation, we use an energy balance to find an effective, thermodynamic, contact angle for displacement and show that this angle increases during the displacement. Displacement occurs when the local contact angles overcome the advancing contact angles at a pinned interface: it is wettability which controls the filling sequence. The product of the principal interfacial curvatures, the Gaussian curvature, is negative, implying well-connected phases which is consistent with pinning at the contact line while providing a topological explanation for the high displacement efficiencies in mixed-wet media.

Many natural building stones are altered when subjected to freeze-thaw (FT) cycles. On the one hand, the sensitivity of a material to FT damage has been quantified in the past by the existence of a material-specific critical water saturation. On the other hand, it was noticed that natural stones with a high volume of ink-bottle pores, normally hold a relatively large FT resistance. The relationship between the pore structure, the saturation and the FT resistance is however not well understood. In this paper, the influence of water content and the porosity on FT behaviour is investigated macroscopically with temperature and strain measurements, established techniques in FT related research, and X-ray computed microtomography (mu CT). Strain measurements performed on Savonnieres limestone samples clearly show an increasing expansion with increasing water content, with the existence of a critical water saturation level between 70 and 80%. Differential X-ray imaging on differently saturated limestone samples in an unfrozen and frozen state then aided in explaining the origin of this critical saturation degree. By draining water from surrounding micro-pores through cryo-suction, ink-bottle ooid voids served as expansion reservoirs. When the majority of the ooid voids is water saturated prior to freezing, these voids lose their expansion reservoir ability and damaging pressures arise when water freezes inside undrained micropores. These findings not only help to understand the FT resistance of this limestone, but also give insight in the general FT behaviour of materials with bi- or multimodal pore-size distributions. (C) 2020 Elsevier Ltd. All rights reserved.

Hypothesis The development of high-resolution in situ imaging has allowed contact angles to be measured directly inside porous materials. We evaluate the use of concepts in integral geometry to determine contact angle. Specifically, we test the hypothesis that it is possible to determine an average contact angle from measurements of the Gaussian curvature of the fluid/fluid meniscus using the Gauss-Bonnet theorem. Theory and simulation We show that it is not possible to unambiguously determine an average contact angle from the Gauss-Bonnet theorem. We instead present an approximate relationship: 2 π n ( 1 - cos θ ) = 4 π - ∫ κ G 12 dS 12 , where n is the number of closed loops of the three-phase contact line where phases 1 and 2 contact the surface, θ is the average contact angle, while κ G 12 is the Gaussian curvature of the fluid meniscus which is integrated over its surface S 12 . We then use the results of pore-scale lattice Boltzmann simulations to assess the accuracy of this approach to determine a representative contact angle for two-phase flow in porous media. Findings We show that in simple cases with a flat solid surface, the approximate expression works well. When applied to simulations on pore space images, the equation provides a robust estimate of contact angle, accurate to within 3 ° , when averaged over many fluid clusters, although individual values can have significant errors because of the approximations used in the calculation.

A statistical technique for the pore-scale analyses of heterogeneity and representative elemental volume (REV) in unconventional shale rocks is hereby presented. First, core samples were obtained from shale formations. The images were scanned using microcomputed tomography (micro-CT) machine at 6.7 μm resolution with voxels of 990 × 990 × 1000. These were then processed, digitised, thresholded, segmented and features captured using numerical algorithms. This allows the segmentation of each sample into four distinct morphological entities consisting of pores, organic matter, shale grains and minerals. In order to analyse the degree of heterogeneity, Eagle Ford parallel sample was further cropped into 96 subsamples. Descriptive statistical approach was then used to evaluate the existence of heterogeneity within the subsamples. Furthermore, the Eagle Ford parallel and perpendicular samples were analysed for volumetric entities representative of the petrophysical variable, porosity, using corner point cropping technique. The results of porosity REV for Eagle ford parallel and perpendicular indicated sample representation at 300 μm voxel edge. Both pore volume distribution and descriptive statistical analyses suggested that a wide variation of heterogeneity exists at this scale of investigation. Furthermore, this experiment allows for adequate extraction of necessary information and structural parameters for pore-scale modelling and simulation. Additional studies focusing on re-evaluation at higher resolution are recommended.

We study the impact of brine acidity and initial pore structure on the dynamics of fluid/solid reaction at high Peclet numbers and low Damkohler numbers. A laboratory μ -CT scanner was used to image the dissolution of Ketton, Estaillades, and Portland limestones in the presence of CO 2 -acidified brine at reservoir conditions (10 MPa and 50 °C) at two injected acid strengths for a period of 4 h. Each sample was scanned between 6 and 10 times at ∼4 μm resolution and multiple effluent samples were extracted. The images were used as inputs into flow simulations, and analysed for dynamic changes in porosity, permeability, and reaction rate. Additionally, the effluent samples were used to verify the image-measured porosity changes. We find that initial brine acidity and pore structure determine the type of dissolution. Dissolution is either uniform where the porosity increases evenly both spatially and temporally, or occurs as channelling where the porosity increase is concentrated in preferential flow paths. Ketton, which has a relatively homogeneous pore structure, dissolved uniformly at pH = 3.6 but showed more channelized flow at pH = 3.1. In Estaillades and Portland, increasingly complex carbonates, channelized flow was observed at both acidities with the channel forming faster at lower pH. It was found that the effluent pH, which is higher than that injected, is a reasonably good indicator of effective reaction rate during uniform dissolution, but a poor indicator during channelling. The overall effective reaction rate was up to 18 times lower than the batch reaction rate measured on a flat surface at the effluent pH, with the lowest reaction rates in the samples with the most channelized flow, confirming that transport limitations are the dominant mechanism in determining reaction dynamics at the fluid/solid boundary.

Water flooding is an economic method commonly used in secondary recovery, but a large quantity of crude oil is still trapped in reservoirs after water flooding. A deep understanding of the distribution of residual oil is essential for the subsequent development of water flooding. In this study, a pore-scale model is developed to study the formation process and distribution characteristics of residual oil. The Navier–Stokes equation coupled with a phase field method is employed to describe the flooding process and track the interface of fluids. The results show a significant difference in residual oil distribution at different wetting conditions. The difference is also reflected in the oil recovery and water cut curves. Much more oil is displaced in water-wet porous media than oil-wet porous media after water breakthrough. Furthermore, enhanced oil recovery (EOR) mechanisms of both surfactant and polymer flooding are studied, and the effect of operation times for different EOR methods are analyzed. The surfactant flooding not only improves oil displacement efficiency, but also increases microscale sweep efficiency by reducing the entry pressure of micropores. Polymer weakens the effect of capillary force by increasing the viscous force, which leads to an improvement in sweep efficiency. The injection time of the surfactant has an important impact on the field development due to the formation of predominant pathway, but the EOR effect of polymer flooding does not have a similar correlation with the operation times. Results from this study can provide theoretical guidance for the appropriate design of EOR methods such as the application of surfactant and polymer flooding., Energies, 12 (19), ISSN:1996-1073

Fractures are conduits that can enable fast advective transfer of (fluid, solute, reactant, particle, etc.) mass and energy. Such fast transfer can significantly affect pore-scale physico-chemical processes, which can in turn affect macroscopic mass and energy transport characteristics. Here, flooding experiments are conducted in a well-characterized fractured porous medium, manufactured by 3D printing. Given steady-state flow conditions, the micro-structure of the two-dimensional pore fluid flow field is delineated to resolve fluid velocities on the order of a sub-millimeter per second. We demonstrate the capabilities of a new temporo-ensemble particle image velocimetry method by maximizing its spatial resolution, employing in-line illumination. This method is advantageous as it is capable of minimizing the number of pixels, required for velocity determinations, down to one pixel, thereby enabling resolving high spatial resolutions of velocity vectors in a large field of view. While the main goal of this study is to introduce a novel experimental and velocimetry framework, this new method is then applied to specifically improve the understanding of fluid flow through fractured porous media. Histograms of measured velocities indicate log-normal and Gaussian-type distributions of longitudinal and lateral velocities in fractures, respectively. The magnitudes of fluid velocities in fractures and the flow interactions between fractures and matrices are shown to be influenced by the permeability of the background matrix and the orientation of the fractures.

Conservation of energy is used to derive a thermodynamically-consistent contact angle, θ t , when fluid phase 1 displaces phase 2 in a porous medium. Assuming no change in Helmholtz free energy between two local equilibrium states we find that Δ a 1 s cos θ t = κ ϕ Δ S 1 + Δ a 12 , where a is the interfacial area per unit volume, ϕ is the porosity, S is the saturation and κ the curvature of the fluid-fluid interface. The subscript s denotes the solid, and we consider changes, Δ , in saturation and area. With the advent of high-resolution time-resolved three-dimensional X-ray imaging, all the terms in this expression can be measured directly. We analyse imaging datasets for displacement of oil by water in a water-wet and a mixed-wet sandstone. For the water-wet sample, the curvature is positive and oil bulges into the brine with almost spherical interfaces. In the mixed-wet case, larger interfacial areas are found, as the oil resides in layers. The mean curvature is close to zero, but the interface tends to bulge into brine in one direction, while brine bulges into oil in the other. We compare θ t with the values measured geometrically in situ on the pore-scale images, θ g . The thermodynamic angle θ t provides a robust and consistent characterization of wettability. For the water-wet case the calculated value of θ t gives an accurate prediction of multiphase flow properties using pore-scale modelling.

Recent advances in high resolution X-ray tomography (μCT) technology have enabled in-situ dynamic μCT imaging (4D-μCT) of time-dependent processes inside 3D structures, non-destructively and non-invasively. This paper illustrates the application of 4D-μCT for visualizing the removal of fatty liquids from kitchen sponges made of polyurethane after rinsing (absorption), squeezing (desorption) and cleaning (adding detergents). For the first time, time-dependent imaging of this type of system was established with sufficiently large contrast gradient between water (with/without detergent) and olive oil (model fat) by the application of suitable fat-sensitive X-ray contrast agents. Thus, contrasted olive oil filled sponges were rinsed and squeezed in a unique laboratory loading device with a fluid flow channel designed to fit inside a rotating gantry-based X-ray μCT system. Results suggest the use of brominated vegetable oil as a preferred contrast agent over magnetite powder for enhancing the attenuation coefficient of olive oil in a multi fluid filled kitchen sponge. The contrast agent (brominated vegetable oil) and olive oil were mixed and subsequently added on to the sponge. There was no disintegration seen in the mixture of contrast agent and olive oil during the cleaning process by detergents. The application of contrast agents also helped in accurately tracking the movement and volume changes of soils in compressed open cell structures. With the in house-built cleaning device, it was quantified that almost 99% of cleaning was possible for contrasted olive oil (brominated vegetable oil with olive oil) dispersed in the sponge. This novel approach allowed for realistic mimicking of the cleaning process and provided closer evaluation of the effectiveness of cleaning by detergents to minimize bacterial growth.

We perform direct numerical simulation using a pore-scale fluid–fluid reactive transport model (Alhashmi et al. in J Contam Hydrol 179:171–181, 2015. doi: 10.1016/j.jconhyd.2015.06.004 ) to investigate the impact of pore structure heterogeneity on the effective reaction rate in different porous media. We simulate flow, transport, and reaction in three pore-scale images: a beadpack, Bentheimer sandstone, and Doddington sandstone for a range of transport and reaction conditions. We compute the reaction rate, velocity distributions, and dispersion coefficient comparing them with the results for non-reactive transport. The rate of reaction is a subtle combination of the amount of mixing and spreading that cannot be predicted from the non-reactive dispersion coefficient. We demonstrate how the flow field heterogeneity affects the effective reaction rate. Dependent on the intrinsic flow field heterogeneity characteristic and the flow rate, reaction may: (a) occur throughout the zones where both resident and injected particles exist (for low Peclet number and a homogeneous flow field), (b) preferentially occur at the trailing edge of the plume (for high Peclet number and a homogeneous flow field), or (c) be disfavored in slow-moving zones (for high Peclet number and a heterogeneous flow field).

We investigate the impact of initial pore structure and velocity field heterogeneity on the dynamics of fluid/solid reaction at high Peclet numbers (fast flow) and low Damkohler number (relatively slow reaction rates). The Diamond Lightsource Pink Beam was used to image dissolution of Estaillades and Portland limestones in the presence of CO 2 -saturated brine at reservoir conditions (10 MPa and 50 °C representing ~ 1 km aquifer depth) at two flow rates for a period of 2 h. Each sample was scanned between 51 and 94 times at 4.76-μm resolution and the dynamic changes in porosity, permeability, and reaction rate were examined using image analysis and flow modelling. We find that the porosity can either increase uniformly through time along the length of the samples, or may exhibit a spatially and temporally varying increase that is attributed to channel formation, a process that is distinct from wormholing, depending on initial pore structure and flow conditions. The dissolution regime was structure-dependent: Estaillades with a higher porosity showed more uniform dissolution, while the lower porosity Portland experienced channel formation. The effective reaction rates were up to two orders of magnitude lower than those measured on a flat substrate with no transport limitations, indicating that the transport of reactant and product is severely hampered away from fast flow channels.

Hypothesis: Capillary-dominated multiphase flow in porous materials is strongly affected by the pore walls' wettability. Recent micro-computed tomography (mCT) studies found unexpectedly wide contact angle distributions measured on static fluid distributions inside the pores. We hypothesize that analysis on time-resolved mCT data of fluid invasion events may be more directly relevant to the fluid dynamics. Experiment: We approximated receding contact angles locally in time and space on time-resolved mCT datasets of drainage in a glass bead pack and a limestone. Whenever a meniscus suddenly entered one or more pores, geometric and thermodynamically consistent contact angles in the surrounding pores were measured in the time step just prior to the displacement event. We introduced a new force-based contact angle, defined to recover the measured capillary pressure in the invaded pore throat prior to interface movement. Findings: Unlike the classical method, the new geometric and force-based contact angles followed plausible, narrower distributions and were mutually consistent. We were unable to obtain credible results with the thermodynamically consistent method, likely because of sensitivity to common imaging artifacts and neglecting dissipation. Time-resolved mCT analysis can yield a more appropriate wettability characterization for pore scale models, despite the need to further reduce image analysis uncertainties. (C) 2020 The Authors. Published by Elsevier Inc.

During CO2 storage in depleted oil fields, under immiscible conditions, CO2 can be trapped in the pore space by capillary forces, providing safe storage over geological times - a phenomenon named capillary trapping. Synchrotron X-ray imaging was used to obtain dynamic three-dimensional images of the flow of the three phases involved in this process - brine, oil and gas (nitrogen) - at high pressure and temperature, inside the pore space of Ketton limestone. First, using continuous imaging of the porous medium during gas injection, performed after waterflooding, we observed chains of multiple displacements between the three phases, caused by the connectivity of the pore space. Then, brine was re-injected and double capillary trapping - gas trapping by oil and oil trapping by brine - was the dominant double displacement event. We computed pore occupancy, saturations, interfacial area, mean curvature and Euler characteristic to elucidate these double capillary trapping phenomena, which lead to a high residual gas saturation. Pore occupancy and saturation results show an enhancement of gas trapping in the presence of both oil and brine, which potentially makes CO2 storage in depleted oil reservoirs attractive, combining safe storage with enhanced oil recovery through immiscible gas injection. Mean curvature measurements were used to assess the capillary pressures between fluid pairs during double displacements and these confirmed the stability of the spreading oil layers observed, which facilitated double capillary trapping. Interfacial area and Euler characteristic increased, indicating lower oil and gas connectivity, due to the capillary trapping events.

The vanadium redox flow battery (VRFB) has emerged as a promising technology for large-scale storage of intermittent power generated from renewable energy sources due to its advantages such as scalability, high energy efficiency and low cost. In the current study, a three-dimensional(3D) Lattice Boltzmann model is developed to simulate the transport mechanisms of electrolyte flow, species and charge in the vanadium redox flow battery at the micro pore scale. An electrochemical model using the Butler-Volmer equation is used to provide species and charge coupling at the surface of active electrode. The detailed structure of the carbon paper electrode is obtained using X-ray Computed Tomography(CT). The new model developed in the paper is able to predict the local concentration of different species, over-potential and current density profiles under charge/discharge conditions. The simulated capillary pressure as a function of electrolyte volume fraction for electrolyte wetting process in carbon paper electrode is presented. Different wet surface area of carbon paper electrode correspond to different electrolyte volume fraction in pore space of electrode. The model is then used to investigate the effect of wetting area in carbon paper electrode on the performance of vanadium redox flow battery. It is found that the electrochemical performance of positive half cell is reduced with air bubbles trapped inside the electrode.

A model capable of providing a reliable estimation of two-phase flow dynamics and mass-transfer coefficients, is lacking for the design of micropacked-bed reactors via correlations, especially when the particle size of the bed is around 100 μm. In this work, we present a validation of the use of the phase field method for reproducing two-phase flow experiments found in the literature. This numerical simulation strategy sheds light on the impact of the micropacked-bed geometry and wettability on the formation of preferential gas channels. Counterintuitively, to homogenize the two-phase flow hydrodynamics and reduce radial mass-transfer limitations, solvent wettability of the support needs to be restricted, showing best performance when the contact angle ranges to 60° and capillary forces are still dominant. The tuning of gas–liquid–solid interactions by surface wettability modification opens a new window of opportunity for the design and scale-up of micropacked-bed reactors. This research was partially funded by the EU project MAPSYN: Microwave, Acoustic and Plasma SYNtheses, under Grant Agreement No. CP-IP 309376 of the European Union Seventh Framework Program.

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.

In the last decade, the study of fluid flow in porous media has developed considerably due to the combination of X-ray Micro Computed Tomography (micro-CT) and advances in computational methods for solving complex fluid flow equations directly or indirectly on reconstructed three-dimensional pore space images. In this study, we calculate porosity and single phase permeability using micro-CT imaging and Lattice Boltzmann (LB) simulations for 8 different porous media: beadpacks (with bead sizes 50 µm and 350 µm), sandpacks (LV60 and HST95), sandstones (Berea, Clashach and Doddington) and a carbonate (Ketton). Combining the observed porosity and calculated single phase permeability, we shed new light on the existence and size of the Representative Element of Volume (REV) capturing the different scales of heterogeneity from the pore-scale imaging. Our study applies the concept of the ‘Convex Hull’ to calculate the REV by considering the two main macroscopic petrophysical parameters, porosity and single phase permeability, simultaneously. The shape of the hull can be used to identify strong correlation between the parameters or greatly differing convergence rates. To further enhance computational efficiency we note that the area of the convex hull (for well-chosen parameters such as the log of the permeability and the porosity) decays exponentially with sub-sample size so that only a few small simulations are needed to determine the system size needed to calculate the parameters to high accuracy (small convex hull area). Finally we propose using a characteristic length such as the pore size to choose an efficient absolute voxel size for the numerical rock.

A fundamental understanding of flow in porous media at the pore-scale is necessary to be able to upscale average displacement processes from core to reservoir scale. The study of fluid flow in porous media at the pore-scale consists of two key procedures: Imaging - reconstruction of three-dimensional (3D) pore space images; and modelling such as with single and two-phase flow simulations with Lattice-Boltzmann (LB) or Pore-Network (PN) Modelling. Here we analyse pore-scale results to predict petrophysical properties such as porosity, single-phase permeability and multi-phase properties at different length scales. The fundamental issue is to understand the image resolution dependency of transport properties, in order to up-scale the flow physics from pore to core scale. In this work, we use a high resolution micro-computed tomography (micro-CT) scanner to image and reconstruct three dimensional pore-scale images of five sandstones (Bentheimer, Berea, Clashach, Doddington and Stainton) and five complex carbonates (Ketton, Estaillades, Middle Eastern sample 3, Middle Eastern sample 5 and Indiana Limestone 1) at four different voxel resolutions (4.4 µm, 6.2 µm, 8.3 µm and 10.2 µm), scanning the same physical field of view. Implementing three phase segmentation (macro-pore phase, intermediate phase and grain phase) on pore-scale images helps to understand the importance of connected macro-porosity in the fluid flow for the samples studied. We then compute the petrophysical properties for all the samples using PN and LB simulations in order to study the influence of voxel resolution on petrophysical properties. We then introduce a numerical coarsening scheme which is used to coarsen a high voxel resolution image (4.4 µm) to lower resolutions (6.2 µm, 8.3 µm and 10.2 µm) and study the impact of coarsening data on macroscopic and multi-phase properties. Numerical coarsening of high resolution data is found to be superior to using a lower resolution scan because it avoids the problem of partial volume effects and reduces the scaling effect by preserving the pore-space properties influencing the transport properties. This is evidently compared in this study by predicting several pore network properties such as number of pores and throats, average pore and throat radius and coordination number for both scan based analysis and numerical coarsened data.

We present a pore scale model capable of simulating fluid/fluid reactive transport on images of porous media from first principles. We use a streamline-based particle tracking method for simulating flow and transport, while for reaction to occur, both reactants must be within a diffusive distance of each other during a time-step. We assign a probability of reaction (Pr), as a function of the reaction rate constant (kr) and the diffusion length. Firstly, we validate our model for reaction against analytical solutions for the bimolecular reaction (A+B→C) in a free fluid. Then, we simulate transport and reaction in a beadpack to validate the model through predicting the fluid/fluid reaction experimental results provided by Gramling et al. (2002). Our model accurately predicts the experimental data, as it takes into account the degree of incomplete mixing present at the sub-pore (image voxel) level, in contrast to advection–dispersion–reaction equation (ADRE) model that over-predicts pore scale mixing. Finally, we show how our model can predict dynamic changes in the reaction rate accurately accounting for the local geometry, topology and flow field at the pore scale. We demonstrate the substantial difference between the predicted early-time reaction rate in comparison to the ADRE model.