VS2DRTI: Simulating Heat and Reactive Solute Transport in Variably Saturated Porous Media          下载免费PDF全文

Richard W. Healy  Sosina S. Haile  David L. Parkhurst  Scott R. Charlton 《Ground water》2018,56(5):810-815

Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid‐rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult‐to‐use models. To address the need for a simple and easy‐to‐use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two‐dimensional, constant‐density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature‐dependent cation exchange. VS2DRTI is freely available public domain software.  相似文献   

Surface and subsurface water balance estimation by the groundwater recharge model and a 3-D two-phase flow model/Estimation de bilan hydrologique de surface et de subsurface à l’aide de modèles de recharge de nappe et d’écoulement diphasique 3-D     

《水文科学杂志》2013,58(2)

Abstract

Abstract Land development often results in adverse environmental impact for surface and subsurface water systems. For areas close to the coast, land changes may also result in seawater intrusion into coastal aquifers. Due to this, it is important to evaluate potential adverse effects in advance of any land development. For evaluation purposes a combined groundwater recharge model is proposed with a quasi three-dimensional unconfined groundwater flow equation. The catchment water balance for a planned new campus area of Kyushu University in southern Japan, was selected as a case study to test the model approach. Since most of the study area is covered with forest, the proposed groundwater recharge model considers rainfall interception by forest canopy. The results show that simulated groundwater and surface runoff agree well with observations. It is also shown that actual evapotranspiration, including rainfall interception by forest canopy, is well represented in the proposed simulation model. Several hydrological components such as direct surface runoff rate, groundwater spring flow rate to a ground depression, trans-basin groundwater flow etc., were also investigated.  相似文献   

1DTempPro: Analyzing Temperature Profiles for Groundwater/Surface‐water Exchange     

Emily B. Voytek  Anja Drenkelfuss  Frederick D. Day‐Lewis  Richard Healy  John W. Lane Jr.  Dale Werkema 《Ground water》2014,52(2):298-302

A new computer program, 1DTempPro, is presented for the analysis of vertical one‐dimensional (1D) temperature profiles under saturated flow conditions. 1DTempPro is a graphical user interface to the U.S. Geological Survey code Variably Saturated 2‐Dimensional Heat Transport (VS2DH), which numerically solves the flow and heat‐transport equations. Pre‐ and postprocessor features allow the user to calibrate VS2DH models to estimate vertical groundwater/surface‐water exchange and also hydraulic conductivity for cases where hydraulic head is known.  相似文献   

Recharge process and aquifer models of a small watershed     

A. SRINIVAS  B. VENKATESWARA RAO  V. V. S. GURUNADHA RAO 《水文科学杂志》2013,58(5):681-692

Abstract

Kanchanapally watershed covering an area of about 11 km² in Nalgonda district, Andhra Pradesh, India is located in granitic terrain. Groundwater recharge has been estimated from a water balance model using hydrometeorological data from 1978–1994. The monthly recharge estimates obtained from the water balance model formed input for the groundwater flow model during transient model testing. The groundwater flow model has been prepared to simulate steady state groundwater conditions of 1977 using the nested squares finite difference method. The transient groundwater flow model has been tested during 1977–1994 using the estimated recharge values. The present study helped verify the usefulness of monthly recharge estimates for accounting dynamic variations in recharge as reflected in water level fluctuations in hydrographs.  相似文献   

Performance Analysis of the χMD Matrix Solver Package for MODFLOW-USG     

Motomu Ibaraki  Yiding Zhang  Richard G. Niswonger  Sorab Panday 《Ground water》2021,59(6):905-912

The χMD matrix solver package is incorporated into USGS groundwater modeling software, such as MODFLOW-NWT, MODFLOW-USG, and MT3D. The solver is used to solve matrices assembled through numerical discretization of the groundwater flow equation, and solute transport equations. χMD has demonstrated its higher robustness, faster execution speed, and more efficient memory usage compared to the existing solvers for many types of groundwater flow problems. χMD uses preconditioned iterative Krylov-subspace methods and consists of preconditioning and acceleration modules. Because the solver package uses a variety of preconditioning features including level-based incomplete lower-upper (ILU) factorization method with a drop tolerance scheme, users must choose optimal preconditioning parameters to improve execution speed and robustness. In order to examine how the preconditioning parameters, ILU factorization level, and drop tolerance values affect the overall performance of the matrix solver, we evaluated five different groundwater model applications using MODFLOW-USG that include different numerical complexities. For those five cases, the number of discretization nodes varied from 10,000 cells to 730,300 cells. From the analysis, we found that the preconditioning parameters greatly affect execution times and memory usage of the preconditioning and acceleration procedures. In addition, a combination of the ILU level between five to seven and the drop tolerance value between 10⁻² and 10⁻³ usually resulted in shorter overall execution time. Our study suggests that the users can elicit higher performance and robustness of the χMD matrix solver using this combination of the parameters and enhance computational efficiency of solving groundwater and solute transport problems.  相似文献   

Effects of multi-annual climate variability on the hydrodynamic evolution (1833 to present) in a shallow aquifer system in northern Belgium     

Marc Van Camp  Marleen Coetsiers  Kristine Martens  Kristine Walraevens 《水文科学杂志》2013,58(5):763-779

Abstract

Using a groundwater flow model and long historical meteorological time series data, the evolution of the groundwater flow regime in a multi-layered groundwater flow basin in northern Belgium during the last one and a half centuries (since 1833) is reconstructed. Model output parameters such as piezometric levels, depth to water table, seepage fluxes in the valleys and calculated baseflow to the river system are presented and inter-annual and decadal variations are evaluated against seasonal fluctuations. The main time-varying boundary condition in the model is the aquifer recharge which was estimated using the method of Thornthwaite and Mather based on precipitation and temperature data. The model does not take into account changes in boundary conditions due to changes in land use (deforestation, drainage of cultivated land) or groundwater exploitation. Variations in model output parameters are therefore only due to climatological forcing. Only the natural non-exploited state of the aquifer is considered. Although few historical piezometric measurements are available to verify model output, the results give an indication of the natural hydrodynamic variations on a time scale of decades.

Citation Van Camp, M., Coetsiers, M., Martens, K. & Walraevens, K. (2010) Effects of multi-annual climate variability on the hydrodynamic evolution (1833 to present) in a shallow aquifer system in northern Belgium. Hydrol. Sci. J. 55(5), 763–779.  相似文献   

Incorporating groundwater recharge and discharge functions into an existing monthly rainfall–runoff model/Incorporation de fonctions de recharge et de vidange superficielle de nappes au sein d’un modèle pluie-débit mensuel existant     

《水文科学杂志》2013,58(2)

Abstract

Abstract There is an urgent need for an integrated surface water and groundwater modelling tool that is suitable for southern African conditions and can be applied at various basin scales for broad strategic water resource planning purposes. The paper describes two new components (recharge and groundwater discharge) that have been added to an existing monthly time-step rainfall–runoff model that is widely used in the southern African subcontinent. The new components are relatively simple, consistent with the existing model formulation, but based on accepted groundwater flow principles and well understood groundwater parameters. The application of the revised model on two basins in southern Africa with quite different baseflow characteristics has demonstrated that the new components have a great deal of potential, even if the improvement is only to be able to simulate the groundwater baseflow component of total runoff more explicitly. More comprehensive testing and comparison of the results with existing groundwater and geohydrological data is required, while some extensions to the new components need to be considered to ensure that the model can be considered applicable to a wide range of basin and climate types.  相似文献   

OPTIMAL MANAGEMENT OF THE SUBSURFACE ENVIRONMENT / L'aménagement optimal de l'environnement souterrain     

ROBERT WILLIS 《水文科学杂志》2013,58(2):333-344

Abstract

A mathematical model for the optimal, conjunctive management of a groundwater basin's quality and quantity resources is presented. The quality resource, the assimilative waste capacity, is the ability of the aquifer to degrade certain non-conservative constituents through the actions of molecular diffusion, hydrodynamic dispersion, adsorption, biochemical reactions and convective mass transport. The model is applicable to saturated, isothermal, porous media. The Galerkin method is used to express the spatial and temporal variations in hydraulic head and constituent concentrations as a function of possible recharge, pumping, and treatment decisions. The Galerkin procedure ensures that the constraint equations actually model the response characteristics (flow and mass transport) of the groundwater basin. In each planning period, the groundwater basin is required to meet an exogenous water demand while maintaining adequate water quality levels throughout the aquifer. The model, structured as a mathematical program, minimizes the sum of waste water treatment and pumping and recharge costs. The results indicate (1) the feasibility of conjunctive management of an aquifer's quality and quantity resources, and (2) that the assimilative waste capacity can function as a form of secondary or advanced waste water treatment.  相似文献   

Groundwater flow through a Miliolite limestone aquifer     

K. R. RUSHTON  S. V. RAGHAVA RAO 《水文科学杂志》2013,58(5):449-464

Abstract

This paper describes a study of groundwater flow in a coastal Miliolite limestone aquifer in western India. An examination of field information suggested that the transmissivity of the aquifer varies significantly between high and low groundwater heads. Pumping tests indicate that this is due to the development of major fissures in the upper part of the aquifer. A regional groundwater model with varying transmissivities is used to represent the field behaviour. The model is also used to examine the effect of artificial recharge on the alleviation of saline intrusion problems in the coastal area.  相似文献   

Tracing groundwater flow using chemical data     

I. VARSÁNYI 《水文科学杂志》2013,58(3):265-275

Abstract

In the southern part of the Great Hungarian Plain there are two different groundwater flow systems located in the study area. Between these two systems there is a narrow groundwater divide where groundwater flow cannot be detected. In the western part, groundwater of greater hardness moves from the northwest, west and southwest. In the eastern part, groundwater flow is from southeast to northwest. The directions of groundwater flow have been established on the basis of dissolved mono- and di-valent cation concentrations. The major direction of groundwater flow detected by a statistical evaluation of water chemical data agrees with previous geological investigations.  相似文献   

Modelling solute leaching during fingered flow by integrating and expanding various theoretical and empirical concepts     

GERRIT H. DE ROOIJ  HIROYUKI CHO 《水文科学杂志》2013,58(3):447-465

Abstract

Wetting front instability (fingered flow) accelerates solute transport through the unsaturated zone to the groundwater table. Whether fingers widen or dissipate close to the groundwater is unclear. Water flow in a two-dimensional artificial capillary fringe below a dry layer exhibiting fingered flow was investigated. The flow diverged strongly in the wet soil, suggesting that fingers dissipate. Expressions for the finger radius in dry and wet soil were combined and adapted to a soil hydraulic property parameterization popular in numerical modelling. The modified equation provided finger radii for soils in humid and arid climates. The fingers in the arid soil were excessively wide. The finger radii were used to model solute transport, assuming fingers dissipated in the subsoil. Modelling was cumbersome for the arid climate. One shower may often be insufficient to trigger fingering in arid regions with short, heavy showers. In soils with shallow groundwater, the diverging subsoil flow determines solute leaching.  相似文献   

Flow–duration curve integration into digital filtering algorithms for simulating climate variability based on river baseflow     

R. Mohammed 《水文科学杂志》2013,58(10):1558-1573

ABSTRACT

A baseflow separation methodology combining the outcomes of the flow–duration curve and the digital filtering algorithms to cope with the restrictions of traditional procedures has been assessed. Using this methodology as well as the monitored and simulated hydro-climatological data, the baseflow annual variations due to climate change and human-induced activities were determined. The outcomes show that the long-term baseflow index at the upstream sub-basin is nearly half of that at the downstream from October to April, whereas they are close to each other for the remaining months. Some of the groundwater reacts to precipitation and an evident rise in the groundwater contribution was detected for the hydrological years 1998–2001 and 2006–2008. The contrary was recorded for1987. The water released from the reservoir in the dry periods led to distinctions in the detected baseflow index between the pre-damming and post-damming periods of the river.  相似文献   

Application of geochemical and stable isotopic tracers to investigate groundwater salinity in the Ochi-Narkwa Basin,Ghana   总被引：1，自引：1，他引：0  

Samuel Y. Ganyaglo  Shiloh Osae  Tetteh Akiti  Thomas Armah  Laurence Gourcy  Tomas Vitvar 《水文科学杂志》2017,62(8):1301-1316

Rainwater, groundwater and soil-water samples were analysed to assess groundwater geochemistry and the origin of salinity in the Ochi-Narkwa basin of the Central Region of Ghana. The samples were measured for major ions and stable isotopes (δ¹⁸O, δ²H and δ¹³C). The Cl^? content in rainwater decreased with distance from the coast. The major hydrochemical facies were Na-Cl for the shallow groundwaters and Ca-Mg-HCO₃, Na-Cl and Ca-Mg-Cl-SO₄ for the deep groundwaters. Groundwater salinization is caused largely by halite dissolution and to a minor extent by silicate weathering and seawater intrusion. Stable isotope composition of the groundwaters followed a slope of 3.44, suggesting a mixing line. Chloride profiles in the soil zone revealed the existence of salt crusts, which support halite dissolution in the study area. A conceptual flow model developed to explain the mechanism of salinization showed principal groundwater flow in the NW–SE direction.