Experimental study of thermally enhanced recovery of high-viscosity DNAPL in saturated porous media under non-isothermal conditions

Fractal geometry of fingering front morphologies of water and NAPL flow in homogeneous porous media: A comparison study

Fingering front morphologies for water and non-aqueous phase liquids (NAPLs) infiltrating into homogeneous unsaturated porous media were quantitatively described based on the fractal assumption. Correlations of fractal dimensions with physical properties of the fluids were studied. The implications of fractal dimensions for environmental systems are further discussed. Fingering front morphologies had fractal properties, and diesels which with high capillary numbers experienced high fractal dimensions than water. Fractal dimension was suggested as an indicator for dye coverage and infiltration depth, which represent pollution area and depth, respectively. The pollution areas showed positive correlations with fractal dimensions while the infiltration depths showed negative correlations. This information is useful for contaminated soil risk management and important in the effective design of recovery and remediation schemes.

Effects of fluid composition in fluid injection experiments in porous media

Experiments on fluid flow in porous media, using fluids loaded with solids of various grain sizes, have been conducted in a modified Hele-Shaw setup. This setup utilised weakly cemented porous media with specific hydraulic and mechanical properties. Fluid injection in coarse granular media with clean or low-concentration fine particles, results in infiltration only, with pressure close to the material tensile strength, while injection in finer granular material causes damage alongside infiltration, with the fluid pressure still close to the material tensile strength. When larger particle sizes or higher particle concentrations are used in the mixture, the fluid travels further within the porous medium, primarily influenced by the grain size of the granular medium. In the latter case, the Darcy flow equation with an effective permeability term can be employed to determine the pressure differential. For the largest particle sizes included in the fluid, the equation is still applicable, but the effective permeability requires adjustment for particle size within the fluid rather than the granular medium. This is crucial when the injection point is locally clogged. The experiments show that fracturing conditions are controlled by different mechanisms. Dimensional and statistical analysis was used to classify the injection pressures to regimes predicted by fracturing theory or by Darcy law with modified effective permeabilities. The findings show that both the material properties and fluid composition are important designing parameters.

Pore-scale investigation of surfactant-enhanced DNAPL mobilization and solubilization

Surfactant-enhanced aquifer remediation has been proved successful to remove dense non-aqueous phase liquids (DNAPLs) from contaminated sites. However, the underlying mechanisms of the DNAPL mobilization and solubilization at the pore scale remains to be addressed for efficient application to the field remediation system. In this work, the emerging microfluidic and imaging technologies are applied to investigate the dynamics of DNAPL remediation. Visualized experiments of the evolution of DNAPL remediation are performed to study the role of surfactant type, concentration and injection rate. The DNAPL remediation is dominated by mobilization followed by solubilization for most surfactants. Mobilization occurs as soon as surfactants and DNAPL are in contact until forming a new stable phase structure, and the solubilization continues until the end of injection. We observe the breakup behavior of long droplets and ganglia during the mobilization, which is attributed to the surfactant-reduced interfacial tension and thus expedites DNAPL mobilization and redistribution. During the solubilization, the formation of micelles incorporating DNAPL fractions increases the DNAPL concentration gradient and thus enhances the mass transfer, but the rate-limited diffusion of micelles reduces the mass transfer rate coefficient. Increasing the surfactant content and decreasing the injection rate can promote mobilization and solubilization. The DNAPL mobilization ability of the surfactants SDS and SDBS is stronger than SAOS and Tween 80 regardless of the injection rates. Tween 80 may be considered an ideal surfactant of only solubilization but not mobilization is desired. This work elucidates the pore-scale mechanisms during surfactant-enhanced DNAPL remediation, which are beneficial for upscaling studies, predictive modeling, and operation optimization of DNAPL remediation in the field.

The effects of water table fluctuation on LNAPL deposit in highly permeable porous media: A coupled numerical and experimental study

Light Non-Aqueous Phase Liquid (LNAPL) flow on the water table is highly mobile and is sensitive to the fluctuation of groundwater. This process is highly complex and involves the migration of three immiscible phases (i.e. water, LNAPL and air) which need the explicit definition of multiple parameters. A coupled experimental and numerical simulation methodology is performed by using Time Domain Reflectrometer (TDR) and multiphase simulation of a controlled environment to mimic the water table fluctuation and its effect on the LNAPL residual saturation. TDR probes are installed in different locations of a 2D tank (i.e. a cuboid box with relatively low off-plane thickness) and the bulk permittivity of the phases are measured through artificially imposed boundary conditions. The bulk permittivity is then translated into saturation of the three different phases. The translated residual saturations along with the previously measured porous media properties (e.g. porosity and saturated permeability) are then inserted into the numerical simulator (i.e. COMSOL Multiphysics®) and the migration of the three phase in porous media is simulated. The numerical exponents and entry pressures needed for the simulation of the multiphase flow are estimated using the temporal experimental values. The exponents of water LNAPL relative permeability were estimated to be around 2 while the exponents gas LNAPL relative permeability were estimated to be closer to 3. The results, simulated with the optimized parameters, are then evaluated with pictures taken from the transparent face of the 2D tank different stages of the experiment. The temporal evolution of different phase saturation has been compared and validated between the experimental results obtained and interpreted by the TDR probe measurements and the simulations. The relative error stays in the 5 % confidence level for most reported points and only in the highly dynamic flow time steps the error reaches around 12% which are discussed in the text and is accepted due to the highly nonlinear nature of the problem.

DNAPL flow and complex electrical resistivity evolution in saturated porous media: A coupled numerical simulation

Induced Polarization (IP) is a non-intrusive geophysical method to monitor Dense Non-Aqueous Phase Liquid (DNAPL) contamination and remediation processes underground. In this study, an advanced numerical code simulating DNAPL flow and complex electrical resistivity is presented. The model was validated against existing IP results and image measurements that were carried out previously in a series of 2D tank experiment. Multiphase flow modeling in porous media is coupled with electrical current modeling to simulate the process of DNAPL migration and the associated IP response. This brings a broader view of the contamination in space and time compared to surface and borehole measurements, especially when the results are supported by field measurements or laboratory experiments. The simulations are developed in 3D and are performed in COMSOL Multiphysics®. The simulations using petrophysical relationships for in-phase and quadrature resistivity and the results of the experiments are in complete accordance with each other in the parts of the tank where the saturation of DNAPL is relatively low (i.e. especially in the cone of depression in the pumping scenario). However, the parts associated with high saturation of DNAPL show high errors between the in-phase resistivity simulations and the results from experiments. The present work can be regarded as a preliminary study toward further applications of coupled IP-multiphase flow for more accurate detection and monitoring of DNAPLs. It is suggested that the choice of tool/approach in this study be extended to larger-scale studies for further investigation.