Pore-scale morphology effects on colloid deposition by trajectory tracking simulations

The migration and deposition of colloids in natural porous media is a phenomenon of high importance in hydrocarbon recovery, wastewater treatment and contaminant hydrogeology. Commonly described by colloid filtration theory (CFT), and using a single representative collector element, the theoretical framework is not accounting for morphological effects and heterogeneity of realistic porous media, including the local variations of balance of forces or presence of hydraulically stagnant regions in the vicinity of grain-to-grain contacts. So far developed numerical models of colloidal transport in porous media accounting for all major interactions and pore structure were applied to 2D systems. Known 3D simulations are limited to particulate flow with a reduced set of interactions.We propose a coupled colloid transport model combining computational fluid dynamics (CFD) and discrete element model (DEM) approaches, where the former provides the velocity field at a given time step, while the second calculates the updated colloid positions accounting for a set of forces acting on each colloid. Relationships between colloid retention and morphological characteristics, such as the pore coordination number and surface to volume ratio of a random sphere pack are investigated by pore partitioning of a digitized representation. The role of stagnation zones is quantified by expressing the capture probability of colloids as a function of distance to the nearest grain-to-grain contact measured along the curved surface.The analysis of pore space morphology of the collector and the simulated filtration process reveals a strong relation between colloid retention and pore connectivity. The simulations successfully mimic the expected effects of ionic strength of the carrier fluid, such as the increasing role of grain-to-grain contacts in colloid retention at low ionic strength.