宏观各向异性多孔介质中流体变形及溶质运移

为深入探究含水层中多孔介质的不规则形态和分布对于地下水流场及溶质运移的作用,根据沉积环境中已被发现的倾斜交叉分层地质单元结构,人为构造多孔介质双层交叉鱼骨结构,使其形成宏观各向异性水力传导特征,通过数值模拟研究该结构及其空间位置对于流场和溶质运移的影响。研究发现:宏观各向异性多孔介质引发了螺旋状流动,导致流线的拉伸和折叠,使得溶质羽发生不规则变形,显著提升了非反应性溶质的稀释程度;宏观各向异性结构和溶质羽在三维空间中的垂向相对位置对于溶质运移具有显著影响。     

Large-diameter,non-pumped wells filled with reactive media for groundwater remediation

A groundwater flow and contaminant transport model was used to simulate arrays of non-pumped wells with reactive media for remediating contaminated groundwater. Each array featured a minimum number of wells, with identical diameter, capable of removing a contaminant plume within a hypothetical site. Simulated well diameters ranged from 0.25 m (similar to typical remediation wells) to 1.25 m (similar to large-diameter, bucket-augered wells). Both arrays occupied a linear transect located approximately 5 m downgradient of the front of a polluted enclave and oriented 90° to the hydraulic gradient. The minimum smallest diameter array contained 23 wells, whereas the minimum largest diameter array contained only four wells. Results of this study suggest that bucket-augering technology, adapted to install non-pumped wells with reactive media, may be an effective alternative for remediating contaminated groundwater in some environments.     

岩溶区地下水数值模拟研究进展

岩溶含水介质的不均一性导致岩溶地下水流动、溶质运移和热量迁移的数学模拟研究成为地下水模拟的难点。本文综述了岩溶区地下水流模拟的几种方法,重点阐述了等效多孔介质法、双重连续介质法和三重介质法的定义、发展过程和适用范围,并回顾了这几种方法的研究成果。从等效多孔介质法到三重介质法,模拟精度不断提高,适用范围也逐渐由大区域实际问题向小区域理论研究过渡。介绍了溶质运移模拟和热迁移模拟的研究方法及实例。溶质运移模拟以对流弥散方程为基础,其中尺度效应是溶质运移模拟的重点研究问题;热量迁移模拟应考虑地下热水密度变化对地下热水运动的影响。溶质运移模拟和热量迁移模拟往往是将迁移模型和已经调试成功的地下水流动模型相耦合,从而达到模拟溶质及热量迁移的目的。由于溶质运移和热量迁移的复杂性,现阶段水流模型多数处于等效多孔介质模型阶段。综合理论及实际应用,指出精确刻画裂隙及管道和注重基础数学算法是岩溶水数值模拟进步的关键。      

Simulation on forecast and control for groundwater contamination of hazardous waste landfill

The groundwater flow and solute transport models were established by Visual Modflow, which was used to forecast the transport process of Cr⁶⁺ in groundwater and simulate the effects of three control measures of contaminants transport after percolation solution leakage happened in the impermeable layer of the landfill. The results showed that the contamination plume of Cr⁶⁺ would reach the pool’s boundary in 10 years, and the distance of contamination transport was 1,450 m. However, the contamination plume will not be obviously expanded between 10 and 20 years. While the ground was covered by hardened concrete, the contamination plume would not reach the pool’s boundary in 20 years. When the leakage-proof barrier was set in the bottom of an unconfined aquifer, the concentration of Cr⁶⁺ was higher than that of the leakage-proof barrier unset, but the result was opposite when setting the leakage-proof barrier in the bottom of confined aquifer. The range of the contamination plume was effectively controlled by setting drainage ditches in which water discharge was 2,298.05 m³ d^?1, which produced monitoring wells which are not contaminated in 20 years. In sum, combining ground hardening with drainage ditches could produce the best effect in controlling contaminants diffusion, and meanwhile, the drainage ditches daily discharge was reduced to 1,710.19 m³ d^?1.     

Numerical study of boundary conditions for solute transport through a porous medium

A transition region may be defined as a region of rapid change in medium properties about the interface between two porous media or at the interface between a porous medium and a reservoir. Modelling the transition region between different porous media can assist in the selection of the most appropriate boundary conditions for the standard advection–dispersion equation (ADE). An advantage of modelling the transition region is that it removes the need for explicitly defining boundary conditions, though boundary conditions may be recovered as limiting cases. As the width of a transition region is reduced, the solution of the transition region model (TR model) becomes equivalent to the solution of the standard ADE model with correct boundary conditions. In this paper numerical simulations using the TR model are employed to select the most appropriate boundary conditions for the standard ADE under a variety of configurations and conditions. It is shown that at the inlet boundary between a reservoir and porous medium, continuity of solute mass flux should be used as the boundary condition. At the boundary interface between two porous media both continuity of solute concentration and solute mass flux should be used. Finally, in a finite porous medium where the solute is allowed to advect freely from the exit point, both continuity of solute concentration and solute mass flux should be used as the outlet boundary condition. The findings made here are discussed with reference to a detailed review of previous relevant theoretical and experimental observations. Copyright © 2001 John Wiley & Sons, Ltd.     

Modeling of transient groundwater flow,pollutant transport,and biodegradation in an aquifer with large hydraulic head variations

Industrially sourced dense non-aqueous phase liquids (DNAPLs) contaminated an alluvial aquifer in France decades ago. The location(s) and nature of the pollution source zone(s) were unknown, and the dissolved concentrations of volatile organic compounds in the monitoring wells varied greatly with time. The aquifer was in hydraulic equilibrium with an artificial canal whose water level was highly variable (up to 5 m). These variations propagated into the aquifer, causing changes in the groundwater flow direction; a transient numerical model of flow and solute transport showed that they correlate with the concentration variations because the changes in the flow direction resulted in the contaminant plume shifting. The transient hydrogeological numerical model was built, taking into account solvent biodegradation with first-order chain, since biodegradation has a significant influence on the pollutant concentration evolution. The model parameterization confirms the position of the source zones among the potential troughs in the bedrock where DNAPLs could have accumulated. The groundwater model was successfully calibrated to reproduce the observed concentration variations over several years and allowed a rapid validation of the hypotheses on the functioning of the polluted system.     

Evaluation of heterogeneity and field-scale groundwater flow regime in a leaky till aquitard

Recent work in southern Ontario, Canada, demonstrates anomalously high vertical groundwater flow velocities (>1 m/year) through a thick (as much as 60 m), sandy silt till aquitard (Northern till), previously assumed to be of very low permeability (hydraulic conductivity <10^–10 m/s). Rapid recharge is attributed to the presence of fractures and sedimentary heterogeneities within the till, but the field-scale flow regime is poorly understood. This study identifies the nature of physical groundwater pathways through the till and provides estimates of the associated groundwater fluxes. The aquitard groundwater flow system is characterized by integrating details of the outcrop and subsurface sedimentary characteristics of the till with field-based hydrogeologic investigation and numerical modeling. Outcrop and subsurface data identify a composite internal aquitard stratigraphy consisting of tabular till beds (till elements) separated by laterally continuous sheet-like sands and gravels (interbeds) and boulder pavements. Individual till elements contain sedimentary heterogeneities, including discontinuous sand and gravel lenses, vertical sand dikes, and zones of horizontal and vertical fractures. Hydrogeologic field investigations indicate a three-layer aquitard flow system, consisting of upper and lower zones of more hydraulically active and heterogeneous till separated by a middle unit of relatively lower hydraulic conductivity. Groundwater pathways and fluxes in the till were evaluated using a two-dimensional aquitard/aquifer flow model which indicates a step-wise flow mechanism whereby groundwater moves alternately downward along vertical pathways (fractures, sedimentary dikes) and laterally along horizontal sand interbeds within the till. This model is consistent with observed hydraulic-head and isotope profiles, and the presence of tritiated pore waters at various depths throughout the till. Simulations suggest that a bulk aquitard vertical hydraulic conductivity on the order of 1×10^–9 m/s is required to reproduce observed hydraulic-head and tritium profiles. Electronic Publication     

Hydrogeophysical approach for contamination assessment in NamSon landfill,Hanoi, Vietnam

Self potential (SP) and electrical resistivity tomography (ERT) methods are used together with the results of groundwater samples hydrogeochemical analysis to assess the impact of the water leak from the landfill garbage site at NamSon located in Northern Hanoi on causing pollution to the surrounding environment and affecting geological structure. Selected survey area covers an area of 180 × 300 m lying in the low land of the NamSon site with a slope ranging about 8 m in direction NW–SE. There are three geophysical measurements lines denoted as T1, T2 and T3. Processing 180 SP data points has allowed to draw maps of equipotential epoch in the two periods in 2015 and 2016. The maps show four zones of SP positive anomalies with maximum amplitudes of about +20 mV where the groundwater flow direction is downward and five zones of SP negative anomalies with minimum values in a range from ?180 to ?260 mV where the groundwater flow direction is upward. Resistivity values of the subsurface layers of soils and rocks have been aquired from 2D inverse model for measuring ERT in March 2015 and March 2016. The results of the ERT allowed to define the low resistivity in the range 15–20 Ωm related to leachate plume from NamSon landfill site. Results of the physico-chemical analysis of groundwater samples from the existing six boreholes show increases in concentration of the measured pollutant parameters indicating contamination of the groundwater as a result of solid waste leachate accumulation. This result is affirmative evidence for the survey results by geophysical technique. The rapid decrease in quality of groundwater over the last year is probably due to the influence of the leachate from the NamSon landfill site.     

The fate and transport of nitrate in shallow groundwater in northwestern Mississippi, USA

Agricultural contamination of groundwater in northwestern Mississippi, USA, has not been studied extensively, and subsurface fluxes of agricultural chemicals have been presumed minimal. To determine the factors controlling transport of nitrate-N into the Mississippi River Valley alluvial aquifer, a study was conducted from 2006 to 2008 to estimate fluxes of water and solutes for a site in the Bogue Phalia basin (1,250 km²). Water-quality data were collected from a shallow water-table well, a vertical profile of temporary sampling points, and a nearby irrigation well. Nitrate was detected within 4.4 m of the water table but was absent in deeper waters with evidence of reducing conditions and denitrification. Recharge estimates from 6.2 to 10.9 cm/year were quantified using water-table fluctuations, a Cl^– tracer method, and atmospheric age-tracers. A mathematical advection-reaction model predicted similar recharge to the aquifer, and also predicted that 15% of applied nitrogen is leached into the saturated zone. With current denitrification and application rates, the nitrate-N front is expected to remain in shallow groundwater, less than 6–9 m deep. Increasing application rates resulting from intensifying agricultural demands may advance the nitrate-N front to 16–23 m, within the zone of groundwater pumping.     

考虑固结效应的溶质传输研究

传统的孔隙介质水动力学采用对流-扩散方程,研究溶质在流体中的迁移。在这个过程中,孔隙介质被认为是不变形的,因而是一个稳态问题。针对二维情况下孔隙介质变形对溶质传输的影响,给出了考虑孔隙介质固结效应的溶质传输方程,并且探讨了该类问题的求解方法。     

Impacts of climate, lake size, and supra- and sub-permafrost groundwater flow on lake-talik evolution, Yukon Flats, Alaska (USA)

In cold regions, hydrologic systems possess seasonal and perennial ice-free zones (taliks) within areas of permafrost that control and are enhanced by groundwater flow. Simulation of talik development that follows lake formation in watersheds modeled after those in the Yukon Flats of interior Alaska (USA) provides insight on the coupled interaction between groundwater flow and ice distribution. The SUTRA groundwater simulator with freeze–thaw physics is used to examine the effect of climate, lake size, and lake–groundwater relations on talik formation. Considering a range of these factors, simulated times for a through-going sub-lake talik to form through 90 m of permafrost range from ～200 to >?1,000  years (vertical thaw rates <?0.1–0.5  m?yr^?1). Seasonal temperature cycles along lake margins impact supra-permafrost flow and late-stage cryologic processes. Warmer climate accelerates complete permafrost thaw and enhances seasonal flow within the supra-permafrost layer. Prior to open talik formation, sub-lake permafrost thaw is dominated by heat conduction. When hydraulic conditions induce upward or downward flow between the lake and sub-permafrost aquifer, thaw rates are greatly increased. The complexity of ground-ice and water-flow interplay, together with anticipated warming in the arctic, underscores the utility of coupled groundwater-energy transport models in evaluating hydrologic systems impacted by permafrost.     

孔隙尺度多孔介质流体流动与溶质运移高性能模拟

深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体(分辨率3.7μm)为研究对象,在空间节点数多达64 000 000的情况下,基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数,并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现,结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD(即10^-13m^2)量级,远低于规则堆叠的圆球结构(>20 000 mD);土壤的迂曲度为1.40~1.60,明显高于规则堆叠的圆球结构。研究结果表明,渗透率大的样品具有较小的迂曲度,这与结构的空间异质性有较强的关系;土壤的渗透率和迂曲度呈现各向异性;在水力梯度一定的前提下,渗透率较大的样品,纵向弥散系数也较大;同时,结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟,可用于土壤多孔介质在孔隙尺度下的水力学特性研究。     

Identifying sources of chlorinated aliphatic hydrocarbons in a residential area in Italy using the integral pumping test method

The results of integral pumping tests (IPTs) performed in the city of Fabriano, Italy, are presented. The IPT methodology was developed by the European Union project INCORE, as a tool for groundwater investigation and source localization in contaminated areas. This methodology consists of a multiple-well pumping test in which the wells are positioned along a control plane downstream of suspected contaminant source zones and perpendicular to the mean groundwater flow direction. During the pumping, concentration time series of target contaminants are measured. In Fabriano, two control planes were realized to identify a chlorinated aliphatic hydrocarbon plume, to estimate the mass fluxes and draw up a ranked list of the main contamination sources. A numerical flow model was implemented to support the IPT design and to interpret the results. This study revealed low-level trichloroethylene contamination (concentration below 8 μg/l), tetrachloroethylene contamination (mean concentration up to 500 μg/l) and a mass flow rate of about 300 g/day. Through the application of the IPT method, the mean contaminant concentrations, the spatial distribution of concentration values along the control planes, and the total contaminant mass flow rates were evaluated, and the investigation area was reduced for further and deeper investigation activities.     

孔隙尺度多孔介质流体流动与溶质运移高性能模拟

深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体（分辨率3.7 μm）为研究对象，在空间节点数多达64 000 000的情况下，基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数，并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现，结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD（即10^-13m²）量级，远低于规则堆叠的圆球结构（>20 000 mD）；土壤的迂曲度为1.40~1.60，明显高于规则堆叠的圆球结构。研究结果表明，渗透率大的样品具有较小的迂曲度，这与结构的空间异质性有较强的关系；土壤的渗透率和迂曲度呈现各向异性；在水力梯度一定的前提下，渗透率较大的样品，纵向弥散系数也较大；同时，结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟，可用于土壤多孔介质在孔隙尺度下的水力学特性研究。     

Numerical modeling of fracking fluid migration through fault zones and fractures in the North German Basin

Gas production from shale formations by hydraulic fracturing has raised concerns about the effects on the quality of fresh groundwater. The migration of injected fracking fluids towards the surface was investigated in the North German Basin, based on the known standard lithology. This included cases with natural preferential pathways such as permeable fault zones and fracture networks. Conservative assumptions were applied in the simulation of flow and mass transport triggered by a high pressure boundary of up to 50 MPa excess pressure. The results show no significant fluid migration for a case with undisturbed cap rocks and a maximum of 41 m vertical transport within a permeable fault zone during the pressurization. Open fractures, if present, strongly control the flow field and migration; here vertical transport of fracking fluids reaches up to 200 m during hydraulic fracturing simulation. Long-term transport of the injected water was simulated for 300 years. The fracking fluid rises vertically within the fault zone up to 485 m due to buoyancy. Progressively, it is transported horizontally into sandstone layers, following the natural groundwater flow direction. In the long-term, the injected fluids are diluted to minor concentrations. Despite the presence of permeable pathways, the injected fracking fluids in the reported model did not reach near-surface aquifers, either during the hydraulic fracturing or in the long term. Therefore, the probability of impacts on shallow groundwater by the rise of fracking fluids from a deep shale-gas formation through the geological underground to the surface is small.     

基于流网单元的污染物优势运移数值模型

垃圾填埋场滤出液、入侵海水、核废物及生产生活废水等污染物随地下水的迁移威胁人类生存。地下水渗流的随机性导致溶质运移问题更加复杂。根据流网特点,利用流线与水头等势线对求解域进行离散,基于质量守恒原理建立流网单元内溶质浓度求解的隐式有限体积差分格式。基于多孔介质孔隙流速分布规律,利用蒙特卡洛法建立流管单元随机流速场进行溶质运移过程的数值模拟,最后根据数值模拟和模型试验的结果对比,验证了数值模拟方法的准确性。基于流网单元的数值模型中沿流线方向的物质交换由对流和扩散共同作用,而流管间的物质交换只有扩散作用,因此,可在不使用弥散系数下进行污染物运移的模拟。引入随机方法确定流管内流速为研究非均匀流场中染污物的优势迁移提供了新的思路。     

Application of Visual MODFLOW to simulation of migration in Cr6+ contaminated site

On the basis of site investigation and data collection of a certain electro plating factory, the groundwater flow and solute transport coupled models were established by applying Visual MODFLOW 4.1 software, which was used to conduct a numerical simulation that forecast the transport process of Cr6+ in groundwater. The results show that contamination plume of Cr6+ transports with groundwater flow direction. Without control measures, in 3 650 days, 19 wells for drinking would be contaminated, and the range of transport would be 52 172 m2, the maximum contamination would be 35.8 mg/L     

Transport and deposition of suspended particles in saturated porous media: hydrodynamic effect

A laboratory study was undertaken to determine the transport and deposition rate of suspended particles in columns of saturated porous media (gravel and glass beads), where the porous media were subjected to steady-state flow. Silt particles with a mode of 14 μm diameter (used as the suspended particles) and fluorescein (as the conservative tracer) were injected into the columns in short pulses. The breakthrough curves were competently described with the analytical solution of a convection–dispersion equation with a first-order deposition rate. The experiments were performed using different flow rates. The suspended particle size distribution, the porous media, and the flow rates themselves were the main factors retained in this study to investigate the mechanisms governing the transport and deposition kinetics in detail. The results showed the existence of a flow rate, beyond which suspended particles travel faster than the conservative tracer. A decrease of the deposition rate of suspended particles beyond a critical flow velocity was also observed. Such behaviour led to consideration of the couple hydrodynamic-gravity forces at high flow rates. As the hydrodynamic force increases, particle deposition rates are reduced due to the effect of hydrodynamic forces inhibiting the deposition.     

The effect of heterogeneity of diffusion parameters on chloride transport in low-permeability argillites

Understanding flow and transport in low-permeability media is very important in the context of nuclear waste disposal, oil and gas reservoirs and long term evolution of groundwater systems. In low-permeability media, transport by diffusion is often the most important mass transport process. This study investigates the effect of the heterogeneity of diffusion parameters on mass transport in low-permeability media. A geostatistical approach for integrating heterogeneity of diffusion parameters in groundwater flow and transport models is proposed and applied to the Toarcian argillites in France which are studied in the framework of feasibility of storing radioactive waste in deep clayey massifs. Stochastic fields of the diffusion parameters of the Toarcian argillites (France) are generated based on 64 measured values of diffusion coefficient and diffusion accessible porosity and used as input for a 3D local-scale groundwater flow and transport model. The chloride concentrations computed by these heterogeneous models are compared to the measured chloride concentrations and to concentrations calculated with a model in which the Toarcian argillites are subdivided into several homogeneous zones. The heterogeneous simulations result in a slightly better correspondence between measured and calculated values and have the additional advantage that the measured diffusion coefficient values in the Toarcian are perfectly honored in the model. This study shows that small-scale variability of diffusion parameters has a significant effect on solute concentrations and omitting this heterogeneity may be a problem in transport calculations in low-permeability media, depending on the specific setting and objectives of the study.     

Exchange and pathways of deep and shallow groundwater in different climate and permafrost conditions using the Forsmark site, Sweden, as an example catchment

This study simulates and quantifies the exchange and the pathways of deep and shallow groundwater flow and solute transport under different climate and permafrost conditions, considering the example field case of the coastal Forsmark catchment in Sweden. A number of simulation scenarios for different climate and permafrost condition combinations have been performed with the three-dimensional groundwater flow and transport model MIKE SHE. Results show generally decreasing vertical groundwater flow with depth, and smaller vertical flow under permafrost conditions than under unfrozen conditions. Also the overall pattern of both the vertical and the horizontal groundwater flow, and the water exchange between the deep and shallow groundwater systems, change dramatically in the presence of permafrost relative to unfrozen conditions. However, although the vertical groundwater flow decreases significantly in the presence of permafrost, there is still an exchange of water between the unfrozen groundwater system below the permafrost and the shallow groundwater in the active layer, via taliks. ‘Through taliks’ tend to prevail in areas that constitute groundwater discharge zones under unfrozen conditions, which then mostly shift to net recharge zones (through taliks with net downward flow) under permafrost conditions.