Numerical Simulation of the Enrichment of Chemotactic Bacteria in Oil-Water Two-Phase Transfer Fields of Heterogeneous Porous Media

Oil pollution in soil-groundwater systems is difficult to remove, and a large amount of residual oil is trapped in the low permeable layer of the heterogeneous aquifer. Aromatic hydrocarbons in oil have high toxicity and low solubility in water, which are harmful to the ecological environment. Chemotactic degrading bacteria can perceive the concentration gradient of non-aqueous phase liquid (NAPL) pollutants in the groundwater environment, and enrich and proliferate around the pollutants, thus achieving a more efficient and thorough remediation effect. However, the existing theoretical models are relatively simple. The physical fields of oil-water two-phase flow and oil-phase solute convection and diffusion in water are not coupled, which further restricts the accuracy of studies on bacterial chemotaxis to NAPL. In this study, geometric models based on the actual microfluidic experimental study were constructed. Based on the phase field model, diffusion convection equation and chemotaxis velocity equation, the effects of heterogeneity of porous media, wall wettability and groundwater flow rate on the residual oil and the concentration distribution of chemotaxis bacteria were studied. Under all of the simulation conditions, the residual oil in the high permeable area was significantly lower than that in the low permeable area, and the wall hydrophilicity enhanced the water flooding effect. Chemotactic bacteria could react to the concentration gradient of pollutants dissolved into water in the oil phase, and enrich near the oil-water interface with high concentration of NAPL, and the density of chemotactic bacteria at the oil-water interface can be up to 1.8-2 times higher than that in the water phase at flow rates from 1.13 to 6.78 m/d.