Effect of injection velocity and particle concentration on transport of nanoscale zero-valent iron and hydraulic conductivity in saturated porous media.

Successful groundwater remediation by injecting nanoscale zero-valent iron (NZVI) particles requires efficient particle transportation and distribution in the subsurface. This study focused on the influence of injection velocity and particle concentration on the spatial NZVI particle distribution, the deposition processes and on quantifying the induced decrease in hydraulic conductivity (K) as a result of particle retention by lab tests and numerical simulations. Horizontal column tests of 2m length were performed with initial Darcy injection velocities (q0) of 0.5, 1.5, and 4.1m/h and elemental iron input concentrations (Fe(0)in) of 0.6, 10, and 17g/L. Concentrations of Fe(0) in the sand were determined by magnetic susceptibility scans, which provide detailed Fe(0) distribution profiles along the column. NZVI particles were transported farther at higher injection velocity and higher input concentrations. K decreased by one order of magnitude during injection in all experiments, with a stronger decrease after reaching Fe(0) concentrations of about 14-18g/kg(sand). To simulate the observed nanoparticle transport behavior the existing finite-element code OGS has been successfully extended and parameterized for the investigated experiments using blocking, ripening, and straining as governing deposition processes. Considering parameter relationships deduced from single simulations for each experiment (e.g. deposition rate constants as a function of flow velocity) one mean parameter set has been generated reproducing the observations in an adequate way for most cases of the investigated realistic injection conditions. An assessment of the deposition processes related to clogging effects showed that the percentage of retention due to straining and ripening increased during experimental run time resulting in an ongoing reduction of K. Clogging is mainly evoked by straining which dominates particle deposition at higher flow velocities, while blocking and ripening play a significant role for attachment, mainly at lower injection velocities. Since the injection of fluids at real sites leads to descending flow velocities with increasing radial distance from the injection point, the simulation of particle transport requires accounting for all deposition processes mentioned above. Thus, the derived mean parameter set can be used as a basis for quantitative and predictive simulations of particle distributions and clogging effects at both lab and field scale. Since decreases in K can change the flow system, which may have positive as well as negative implications for the in situ remediation technology at a contaminated site, a reliable simulation is thus of great importance for NZVI injection and prediction.