Pauw Pieter S, van der Zee Sjoerd E A T M, Leijnse Anton, Oude Essink Gualbert H P
Department of Soil and Groundwater, Deltares, P.O. Box 85467, 3508 AL, Utrecht, The Netherlands.
Department of Soil Physics and Land Management, Wageningen University, P.O. Box 47, Wageningen, 6700 AA, The Netherlands.
Ground Water. 2016 Jul;54(4):521-31. doi: 10.1111/gwat.12382. Epub 2015 Oct 22.
This article deals with the quantification of saltwater upconing below horizontal wells in freshwater lenses using analytical solutions as a computationally fast alternative to numerical simulations. Comparisons between analytical calculations and numerical simulations are presented regarding three aspects: (1) cyclic pumping; (2) dispersion; and (3) finite horizontal wells in a finite domain (a freshwater lens). Various hydrogeological conditions and pumping regimes within a dry half year are considered. The results show that the influence of elastic and phreatic storage (which are not taken into account in the analytical solutions) on the upconing of the interface is minimal. Furthermore, the analytical calculations based on the interface approach compare well with numerical simulations as long as the dimensionless interface upconing is below 1/3, which is in line with previous studies on steady pumping. Superimposing an analytical solution for mixing by dispersion below the well over an analytical solution based on the interface approach is appropriate in case the vertical flow velocity around the interface is nearly constant but should not be used for estimating the salinity of the pumped groundwater. The analytical calculations of interface upconing below a finite horizontal well compare well with the numerical simulations in case the distance between the horizontal well and the initial interface does not vary significantly along the well and in case the natural fluctuation of the freshwater lens is small. In order to maintain a low level of salinity in the well during a dry half year, the dimensionless analytically calculated interface upconing should stay below 0.25.
本文利用解析解对淡水透镜体中水平井下方的盐水上溯进行量化,作为数值模拟在计算上的快速替代方法。在三个方面对解析计算和数值模拟进行了比较:(1)循环抽水;(2)弥散;(3)有限域(淡水透镜体)中的有限水平井。考虑了旱季半年内的各种水文地质条件和抽水工况。结果表明,弹性储存和潜水储存(解析解中未考虑)对界面上溯的影响极小。此外,只要无量纲界面上溯低于1/3,基于界面方法的解析计算与数值模拟结果吻合良好,这与之前关于稳定抽水的研究一致。在界面周围垂直流速近似恒定的情况下,将井下方弥散混合的解析解叠加在基于界面方法的解析解之上是合适的,但不应将其用于估算抽水井中地下水的盐度。当水平井与初始界面之间的距离沿井的方向变化不显著且淡水透镜体的自然波动较小时,有限水平井下方界面上溯的解析计算与数值模拟结果吻合良好。为了在旱季半年内保持井中较低的盐度水平,解析计算得到的无量纲界面上溯应保持在0.25以下。