Dispersion in doublet-type flows through highly anisotropic porous formations

Abstract

Steady doublet-type flow takes place in a porous formation, where the log-transform $Y = \ln K$ of the spatially variable hydraulic conductivity $K$ is regarded as a stationary random field of two-point autocorrelation $\rho _Y$ . A passive solute is injected at the source in the porous formation and we aim to quantify the resulting dispersion process between the two lines by means of spatial moments. The latter depend on the distance $\ell$ between the lines, the variance $\sigma ^2_Y$ of $Y$ and the (anisotropy) ratio $\lambda$ between the vertical and the horizontal integral scales of $Y$ . A simple (analytical) solution to this difficult problem is obtained by adopting a few simplifying assumptions: (i) a perturbative solution, which regards $\sigma ^2_Y$ as a small parameter, of the velocity field is sought; (ii) pore-scale dispersion is neglected; and (iii) we deal with a highly anisotropic formation ( $\lambda \lesssim 0.1$ ). We focus on the longitudinal spatial moment, as it is of most importance for the dispersion mechanism. A general expression is derived in terms of a single quadrature, which can be straightforwardly carried out once the shape of $\rho _Y$ is specified. Results permit one to grasp the main features of the dispersion processes as well as to assess the difference with similar mechanisms observed in other non-uniform flows. In particular, the dispersion in a doublet-type flow is observed to be larger than that generated by a single line. This effect is explained by noting that the advective velocity in a doublet, unlike that in source/line flows, is rapidly increasing in the far field owing to the presence there of the singularity. From the standpoint of the applications, it is shown that the solution pertaining to $\lambda \to 0$ (stratified formation) provides an upper bound for the dispersion mechanism. Such a bound can be used as a conservative limit when, in a remediation procedure, one has to select the strength as well as the distance $\ell$ of the doublet. Finally, the present study lends itself as a valuable tool for aquifer tests and to validate more involved numerical codes accounting for complex boundary conditions.