Enhancing in situ biodegradation in groundwater using pump and treat remediation: a proof of concept and modelling analysis of controlling variables.

A remediation approach which uses pump and treatment (PAT) to enhance the biodegradation of organic contaminants by increasing dispersive mixing between plumes and groundwater was evaluated for a phenol-contaminated aquifer, using a reactive transport model which simulates kinetic reactions between an electron donor (ED) in the plume and electron acceptor (EA) in the groundwater. The influence of system design and operation was examined in six modelling scenarios. Injection or extraction of groundwater increases biodegradation above no action and the location, pumping rate, and distance between well(s) are important variables which influence biodegradation. An increase in pumping rate, distance of the wells from the plume centreline, and changing the flow direction increase dispersive mixing between the plume and groundwater. This increases plume spreading and the plume fringe interface, providing a greater flux of dissolved EAs for biodegradation. In general, injection of groundwater containing natural EAs enhances biodegradation more than extraction. The enhancement of biodegradation is sensitive to the relative fluxes of ED and EA, as controlled by the arrangement of the wells. In the best performing scenario, biodegradation was enhanced by 128%, compared with no action.

Optimal design of pump-and-treat systems under uncertain hydraulic conductivity and plume distribution.

In this work, we present a stochastic optimal control framework for assisting the management of the cleanup by pump-and-treat of polluted shallow aquifers. In the problem being investigated, hydraulic conductivity distribution and dissolved contaminant plume location are considered as the uncertain variables. The framework considers the subdivision of the cleanup horizon in a number of stress periods over which the pumping policy implemented until that stage is dynamically adjusted based upon new information that has become available in the previous stages. In particular, by following a geostatistical approach, we study the idea of monitoring the cumulative contaminant mass extracted from the installed recovery wells, and using these measurements to generate conditional realizations of the hydraulic conductivity field. These realizations are thus used to obtain a more accurate evaluation of the initial plume distribution, and modify accordingly the design of the pump-and-treat system for the remainder of the remedial process. The study indicates that measurements of contaminant mass extracted from pumping wells retain valuable information about the plume location and the spatial heterogeneity characterizing the hydraulic conductivity field. However, such an information may prove quite soft, particularly in the instances where recovery wells are installed in regions where contaminant concentration is low or zero. On the other hand, integrated solute mass measurements may effectively allow for reducing parameter uncertainty and identifying the plume distribution if more recovery wells are available, in particular in the early stages of the cleanup process.

Surface flow boundary conditions in modeling land subsidence due to fluid withdrawal.

Land subsidence due to subsurface fluid (water, gas, oil) withdrawal is often predicted by either finite element or finite difference numerical models based on coupled poroelastic theory, where the soil is represented as a semi-infinite medium bounded by the traction-free (ground) surface. One of the variables playing a most important role on the final outcome is the flow condition used on the traction-free boundary, which may be assumed as either permeable or impermeable. Although occasionally justified, the assumption of no-flow surface seems to be in general rather unrealistic. A permeable boundary where the fluid pressure is fixed to the external atmospheric pressure appears to be more appropriate. This paper addresses the response, in terms of land subsidence, obtained with a coupled poroelastic finite element model that simulates a distributed pumping from a horizontal aquifer confined between two relatively impervious layers, and takes either a permeable boundary surface, i.e., constant hydraulic potential, or an impermeable boundary, i.e., a zero Neumann flow condition. The analysis reveals that land subsidence is rather sensitive to the flow condition implemented on the traction-free boundary. In general, the no-flow condition leads to an overestimate of the predicted ground surface settlement, which could even be 1 order of magnitude larger than that obtained with the permeable boundary.