Characterization of Capillary Pressure–Saturation Relationships for Double-Porosity Medium Using Light Transmission Visualization Technique

Pore-Scale Study of Fluid Displacement in Parallel-Layered Porous Media and Implications of Associated Distinct Flow Dynamics

Fluid displacement within layered porous media is more complex than in nonlayered ones. Most of the previous studies placed a focus on the porous media with layerings perpendicular to the flow direction, and the effects of pore topology were often ignored. Therefore, this study aims to reveal the flow physics in porous media with layering parallel to the flow direction by accounting for the specific pore topology. Based on the phase-field method, a series of pore-scale numerical simulations were performed to investigate the dynamic displacement processes within layered porous media under various capillary numbers and wettability. The results showed that oil recovery was strongly affected by the heterogeneity of porous media in the capillary fingering regime compared to the viscous one. Besides, the alternative advancing phenomenon of the water fingers was suppressed in layered porous media. Compared with water wetting conditions, the viscous fingering regime of the invading water is more pronounced under oil wetting conditions. Last but not least, the center of mass of water flow paths can be used for fingering regime identification. These findings contribute to our better understanding of the flow dynamics of fluid displacement within layered porous media, which are of great significance to the industry related to oil recovery, underground storage of carbon dioxide and hydrogen gas, etc.

Factors influencing residual air saturation during consecutive imbibition processes in an air-water two-phase fine sandy medium - A laboratory-scale experimental study

The residual air saturation plays a crucial role in modeling hydrological processes of groundwater and the migration and distribution of contaminants in subsurface environments. However, the influence of factors such as media properties, displacement history, and hydrodynamic conditions on the residual air saturation is not consistent across different displacement scenarios. We conducted consecutive drainage-imbibition cycles in sand-packed columns under hydraulic conditions resembling natural subsurface environments, to investigate the impact of wetting flow rate, initial fluid state, and number of imbibition rounds (NIR) on residual air saturation. The results indicate that residual air saturation changes throughout the imbibition process, with variations separated into three distinct stages, namely, unstable residual air saturation (Sgr-u), momentary residual air saturation (Sgr-m), and stable residual air saturation (Sgr). The results also suggest that the transition from Sgr-u to Sgr is driven by changes in hydraulic pressure and gradient; the calculated values followed the following trend: Sgr > Sgr-u > Sgr-m. An increase in capillary number, which ranged from 1.46 × 10-7 to 3.07 × 10-6, increased Sgr-u and Sgr-m in some columns. The increase in Sgr ranged from 0.034 to 0.117 across all the experimental columns; this consistent increase can be explained by water film expansion at the primary wetting front along with a strengthening of the hydraulic gradient during water injection. Both the pre-covered water film on the sand grain surface and a pore-to-throat aspect ratio of up to 4.42 were identified as important factors for the increased residual air saturation observed during the imbibition process. Initial air saturation (Sai) positively influenced all three types of residual air saturation, while initial capillary pressure (Pci) exhibited a more pronounced inhibitory effect on residual air saturation, as it can partly characterized the initial connectivity of the air phase generated under different drying flow rates. Under identical wetting flow rate conditions, Sgr was higher during the second imbibition than during the first imbibition due to variations in initial fluid state, involving both fluid distribution and the concentration of dissolved air in the pore water. In contrast, NIR did not have an obvious effect on the three types of residual air saturation. This work aims to provide empirical evidences and offer further insights into the capture of non-wetting phases in groundwater environments, as well as to put forward some potential suggestion for future investigations on the retention and migration of contaminants that involves multiphase interface interactions in subsurface environments.

Towards predicting DNAPL source zone formation to improve plume assessment: Using robust laboratory and numerical experiments to evaluate the relevance of retention curve characteristics

We conducted multiple laboratory trials in a robust and repeatable experimental layout to study dense non-aqueous phase liquid (DNAPL) source zone formation. We extended an image processing and analysis framework to derive DNAPL saturation distributions from reflective optical imaging data, with volume balance deviations < 5.07%. We used a multiphase flow model to simulate source zone formation in a Monte Carlo approach, where the parameter space was defined by the variation of retention curve parameters. Integral and geometric measures were used to characterize the source zones and implemented into a multi-criteria objective function. The latter showed good agreement between observation data and simulation results for effective DNAPL saturation values > 0.04, especially for early stages of DNAPL migration. The common hypothesis that parameters defining the DNAPL-water retention curves are constant over time was not confirmed. Once DNAPL pooling started, the optimal fit in the parameter space was significantly different compared to the earlier DNAPL migration stages. We suspect more complex processes (e.g., capillary hysteresis, adsorption) to become relevant during pool formation. Our results reveal deficits in the grayscale-DNAPL saturation relationship definition and laboratory estimation of DNAPL-water retention curve parameters to overcome current limitations to describe DNAPL source zone formation.