New Model of Reactive Transport in Single-Well Injection-Withdrawal Test with Aquitard Effect

Abstract

Abstract. The model of single-well injection-withdrawal (SWIW) test has been widely used to investigate reactive radial dispersion in remediation or parameter estimation of the in situ aquifers. Previous analytical solutions only focused on a completely isolated aquifer for the SWIW test, excluding any influence of aquitards bounding the tested aquifer. This simplification might be questionable in field applications when test durations are relatively long, because solute transport in or out of the bounding aquitards is inevitable due to molecular diffusion and cross-formational advective transport. Here, a new SWIW model is developed in an aquifer-aquitard system, and the analytical solution in the Laplace domain is derived. Four phases of the test are included: the injection phase, the chaser phase, the rest phase and the extraction phase. The Green's function method is employed for the solution in the extraction phase. As the permeability of aquitard is much smaller than the permeability of the aquifer, the flow is assumed to be perpendicular to the aquitard, thus only vertical dispersive and advective transports are considered for aquitard. The validity of this treatment is tested by a numerical solution. The sensitivity analysis demonstrates that the influence of vertical flow velocity and porosity in the aquitards, and radial dispersion of the aquifer is more sensitive to the SWIW test than other parameters. In the injection phase, the larger radial dispersivity of the aquifer could result in the smaller values of breakthrough curves (BTCs), while greater values of BTCs of the chaser and rest phases. In the extraction phase, it could lead to the smaller peak values of BTCs. The new model of this study performs better than previous studies excluding the aquitard effect for interpreting data of the field SWIW test.