共查询到20条相似文献，搜索用时 31 毫秒
1.
Modern ground water characterization and remediation projects routinely require calibration and inverse analysis of large three-dimensional numerical models of complex hydrogeological systems. Hydrogeologic complexity can be prompted by various aquifer characteristics including complicated spatial hydrostratigraphy and aquifer recharge from infiltration through an unsaturated zone. To keep the numerical models computationally efficient, compromises are frequently made in the model development, particularly, about resolution of the computational grid and numerical representation of the governing flow equation. The compromise is required so that the model can be used in calibration, parameter estimation, performance assessment, and analysis of sensitivity and uncertainty in model predictions. However, grid properties and resolution as well as applied computational schemes can have large effects on forward-model predictions and on inverse parameter estimates. We investigate these effects for a series of one- and two-dimensional synthetic cases representing saturated and variably saturated flow problems. We show that "conformable" grids, despite neglecting terms in the numerical formulation, can lead to accurate solutions of problems with complex hydrostratigraphy. Our analysis also demonstrates that, despite slower computer run times and higher memory requirements for a given problem size, the control volume finite-element method showed an advantage over finite-difference techniques in accuracy of parameter estimation for a given grid resolution for most of the test problems. 相似文献
2.
Flow and transport models in heterogeneous geological formations are usually large-scale with excessive computational complexity and uncertain characteristics. Uncertainty quantification for predicting subsurface flow and transport often entails utilizing a numerical Monte Carlo framework, which repeatedly simulates the model according to a random field parameter representing hydrogeological characteristics of the aquifer. The physical resolution (e.g. spatial grid resolution) for the simulation is customarily chosen based on recommendations in the literature, independent of the number of Monte Carlo realizations. This practice may lead to either excessive computational burden or inaccurate solutions. We develop an optimization-based methodology that considers the trade-off between the following conflicting objectives: time associated with computational costs, statistical convergence of the model prediction and physical errors corresponding to numerical grid resolution. Computational resources are allocated by considering the overall error based on a joint statistical–numerical analysis and optimizing the error model subject to a given computational constraint. The derived expression for the overall error explicitly takes into account the joint dependence between the discretization error of the physical space and the statistical error associated with Monte Carlo realizations. The performance of the framework is tested against computationally extensive simulations of flow and transport in spatially heterogeneous aquifers. Results show that modelers can achieve optimum physical and statistical resolutions while keeping a minimum error for a given computational time. The physical and statistical resolutions obtained through our analysis yield lower computational costs when compared to the results obtained with prevalent recommendations in the literature. Lastly, we highlight the significance of the geometrical characteristics of the contaminant source zone on the optimum physical and statistical resolutions. 相似文献
3.
Widely used numerical models of solute transport processes in subsurface aquifers are limited to nonlocally refined rectangular, or logically rectangular, structured grids. This presents an unsuitable option to efficient numerical simulations maintaining an acceptable level of accuracy. Optimal selection of locally refined cells for efficient solute transport models is challenging to the current generation of numerical models. We present a novel and relatively simple to implement algorithm addressing these shortcomings. This method operates in four steps involving travel times simulations, a grid coarsening stage followed by a selective local grid refinement based on a cell-wise indicator, and a final postprocessing step. The refinement index is the sum of weighted logarithmic distributions of scaled forward and backward travel times. We calculate representative flow and transport properties at the two scales of the composite grid with a flow-based upscaling technique. We present two test problems to demonstrate the performances of this new gridding algorithm. We obtain the most important speedups for composite grids generated with the highest indicator thresholds. When hydrodynamic dispersion effects increase, we obtain less important speedups. An important outcome of this work is that grid design depends on nature and strength of the underlying flow and solute transport processes. Therefore, we suggest developing solute transport workflows integrating this grid generation algorithm as an integral component to build comprehensive and efficient groundwater models. 相似文献
4.
A 3D ERT study of solute transport in a large experimental tank 总被引：2，自引：0，他引：2
L. Slater A. Binley R. Versteeg G. Cassiani R. Birken S. Sandberg 《Journal of Applied Geophysics》2002,49(4)
A high resolution, cross-borehole, 3D electrical resistivity tomography (ERT) study of solute transport was conducted in a large experimental tank. ERT voxels comprising the time sequence of electrical images were converted into a 3D array of ERT estimated fluid conductivity breakthrough curves and compared with direct measurements of fluid conductivity breakthrough made in wells. The 3D ERT images of solute transport behaviour were also compared with predictions based on a 3D finite-element, coupled flow and transport model, accounting for gravity induced flow caused by concentration differences.The tank (dimensions 185×245×186 cm) was filled with medium sand, with a gravel channel and a fine sand layer installed. This heterogeneous system was designed to complicate solute transport behaviour relative to a homogeneous sand tank, and to thus provide a challenging but insightful analysis of the ability of 3D ERT to resolve transport phenomena. Four ERT arrays and 20 piezometers were installed during filling. A NaCl tracer (conductivity 1.34 S/m) was injected and intensively monitored with 3D ERT and direct sampling of fluid chemistry in piezometers.We converted the bulk conductivity estimate for 250 voxels in the ERT imaged volume into ERT estimated voxel fluid conductivity by assuming that matrix conduction in the tank is negligible. In general, the ERT voxel response is in reasonable agreement with the shape of fluid conductivity breakthrough observed in six wells in which direct measurements of fluid conductivity were made. However, discrepancies occur, particularly at early times, which we attribute to differences between the scale of the image voxels and the fluid conductivity measurement, measurement errors mapped into the electrical inversion and artificial image roughness resulting from the inversion.ERT images revealed the 3D tracer distribution at 15 times after tracer injection. The general pattern and timing of solute breakthrough observed with ERT agreed with that predicted from the flow/transport modelling. However, the ERT images indicate a vertical component of tracer transport and preferential flow paths in the medium sand. We attribute this to transient vertical gradients established during tracer injection, and heterogeneity caused by sorting of the sand resulting from the filling procedure. In this study, ERT provided a unique dataset of 250 voxel breakthrough curves in 1.04 m3. The use of 3D ERT to generate an array of densely sampled estimated fluid conductivity breakthrough curves is a potentially powerful tool for quantifying solute transport processes. 相似文献
5.
Modeling effects of multinode wells on solute transport 总被引：1，自引：0，他引：1
Long-screen wells or long open boreholes with intraborehole flow potentially provide pathways for contaminants to move from one location to another in a ground water flow system. Such wells also can perturb a flow field so that the well will not provide water samples that are representative of ground water quality a short distance away from the well. A methodology is presented to accurately and efficiently simulate solute transport in ground water systems that include wells longer than the grid spacing used in a simulation model of the system and hence are connected to multiple nodes of the grid. The methods are implemented in a MODFLOW-compatible solute-transport model and use MODFLOW's Multi-Node Well Package but are generic and can be readily implemented in other solute-transport models. For nonpumping multinode wells (used to simulate open boreholes or observation wells, for example) and for low-rate pumping wells (in which the flow between the well and the ground water system is not unidirectional), a simple routing and local mixing model was developed to calculate nodal concentrations within the borehole. For high-rate pumping multinode wells (either withdrawal or injection, in which flow between the well and the ground water system is in the same direction at all well nodes), complete and instantaneous mixing in the wellbore of all inflows is assumed. 相似文献
6.
《Advances in water resources》1988,11(2):67-73
A Eulerian-Langrangian scheme is used to reformulate the equation of solute transport with ground water in saturated soils. The governing equation is decomposed into advection along characteristic path lines and propagation of the residue at a fixed grid.The method was employed to simulate transport of a conservative pollutant in a hypothetical aquifer, subject to the equivalence of real conditions. Implementation was based on data involving parameters of a heterogeneous aquifer, heavy flux stresses of densed pumpage/recharge wells, precipitation and seasonally changing flow regimes. Simulation, with coarse grid and high Peclet numbers yielded minute mass balance errors. 相似文献
7.
8.
《Advances in water resources》1998,21(1):1-9
The unsteady free surface flow caused by sudden collapse of a dam produces discontinuities in the flow variables. As the flow surges downstream, it forms a moving bore front with steep gradients of water height and velocity. In the numerical simulation of this flow, proper grid distribution can play a crucial part in the prediction and resolution of the solutions. The use of presently available numerical schemes to solve this problem on a uniform course grid system fails to resolve the characteristic flow features and hence do a poor job in simulating this flow. In this paper, an adaptive grid which adjusts itself as the solution evolves is used for a better resolution of the flow properties. Rai and Anderson's12 method is used to determine the grid speed; however, a different partial differential equation based on the conservative principle of grid arc lengths for clustering grids in one-dimensional flow is used along with the St. Venant equations to numerically simulate the flow. Both the subcritical and the supercritical flows under extreme boundary conditions are solved using this technique. With a specified number of grid points, this provides better quality solutions as compared to those obtained with uniformly distributed grids. 相似文献
9.
Cellular‐based approaches for flood inundation modelling have been extensively calibrated and evaluated for the prediction of flood flows on rural river reaches. However, there has only been limited application of these approaches to urban environments, where the need for flood management is greatest. Practical application of two‐dimensional (2D) flood inundation models is often limited by computation time and processing power on standard desktop PCs when attempting to resolve flows on the high‐resolution grids necessary to replicate urban features. Consequently, it is necessary to evaluate the effectiveness of coarse grids to represent flood flows through urban environments. To examine these effects, LISFLOOD‐FP, a 2D storage cell model, is applied to hypothetical flooding scenarios in Greenfields, Glasgow. Grid resampling techniques in GIS software packages are evaluated and a bilinear gridding technique appears to provide the most accurate and physically intuitive results. A gridding method maintaining sharp elevation changes at building interfaces and neighbouring land is presented and estimates of the discretization noise associated with the coarse resolution grids suggest little improvement over current gridding methods. The variation in model results from the friction sensitivity analysis suggests a non‐stationary response to Manning's n with changing model resolution. Model results suggests that a coarse resolution model for urban applications is limited by the representation of urban media in coarse model grids. Furthermore, critical length scales related to building dimensions and building separation distances exist in urban areas that determine maximum possible grid resolutions for hydraulic models of urban flooding. Copyright ©, 2008 John Wiley & Sons, Ltd. 相似文献
10.
Modeling axisymmetric flow and transport 总被引：3，自引：0，他引：3
Langevin CD 《Ground water》2008,46(4):579-590
Unmodified versions of common computer programs such as MODFLOW, MT3DMS, and SEAWAT that use Cartesian geometry can accurately simulate axially symmetric ground water flow and solute transport. Axisymmetric flow and transport are simulated by adjusting several input parameters to account for the increase in flow area with radial distance from the injection or extraction well. Logarithmic weighting of interblock transmissivity, a standard option in MODFLOW, can be used for axisymmetric models to represent the linear change in hydraulic conductance within a single finite-difference cell. Results from three test problems (ground water extraction, an aquifer push-pull test, and upconing of saline water into an extraction well) show good agreement with analytical solutions or with results from other numerical models designed specifically to simulate the axisymmetric geometry. Axisymmetric models are not commonly used but can offer an efficient alternative to full three-dimensional models, provided the assumption of axial symmetry can be justified. For the upconing problem, the axisymmetric model was more than 1000 times faster than an equivalent three-dimensional model. Computational gains with the axisymmetric models may be useful for quickly determining appropriate levels of grid resolution for three-dimensional models and for estimating aquifer parameters from field tests. 相似文献
11.
The U.S. Department of Energy is currently studying Yucca Mountain, Nevada, as a potential site for a geological high-level waste repository. In the current conceptual models of radionuclide transport at Yucca Mountain, part of the transport path to pumping locations would be through an alluvial aquifer. Interactions with minerals in the alluvium are expected to retard the downstream migration of radionuclides, thereby delaying arrival times and reducing ground water concentrations. We evaluate the effectiveness of the alluvial aquifer as a transport barrier using the stochastic Lagrangian framework. A transport model is developed to account for physical and chemical heterogeneities and rate-limited mass transfer between mobile and immobile zones. The latter process is caused by small-scale heterogeneity and is thought to control the macroscopic-scale retardation in some field experiments. A geostatistical model for the spatially varying sorption parameters is developed from a site-specific database created from hydrochemical measurements and a calibrated modeling approach (Turner and Pabalan 1999). Transport of neptunium is considered as an example. The results are sensitive to the rate of transfer between mobile and immobile zones, and to spatial variability in the hydraulic conductivity. Chemical heterogeneity has only a small effect, as does correlation between hydraulic conductivity and the neptunium distribution coefficient. These results illustrate how general sensitivities can be explored with modest effort within the Lagrangian framework. Such studies complement and guide the application of more detailed numerical simulations. 相似文献
12.
When managing large-scale ground water contamination problems, it is often necessary to model flow and transport using finely discretized domains--for instance (1) to simulate flow and transport near a contamination source area or in the area where a remediation technology is being implemented; (2) to account for small-scale heterogeneities; (3) to represent ground water-surface water interactions; or (4) some combination of these scenarios. A model with a large domain and fine-grid resolution will need extensive computing resources. In this work, a domain decomposition-based assembly model implemented in a parallel computing environment is developed, which will allow efficient simulation of large-scale ground water flow and transport problems using domain-wide grid refinement. The method employs common ground water flow (MODFLOW) and transport (RT3D) simulators, enabling the solution of almost all commonly encountered ground water flow and transport problems. The basic approach partitions a large model domain into any number of subdomains. Parallel processors are used to solve the model equations within each subdomain. Schwarz iteration is applied to match the flow solution at the subdomain boundaries. For the transport model, an extended numerical array is implemented to permit the exchange of dispersive and advective flux information across subdomain boundaries. The model is verified using a conventional single-domain model. Model simulations demonstrate that the proposed model operated in a parallel computing environment can result in considerable savings in computer run times (between 50% and 80%) compared with conventional modeling approaches and may be used to simulate grid discretizations that were formerly intractable. 相似文献
13.
Incorporating Super‐Diffusion due to Sub‐Grid Heterogeneity to Capture Non‐Fickian Transport 
下载免费PDF全文
下载免费PDF全文 Numerical transport models based on the advection‐dispersion equation (ADE) are built on the assumption that sub‐grid cell transport is Fickian such that dispersive spreading around the average velocity is symmetric and without significant tailing on the front edge of a solute plume. However, anomalous diffusion in the form of super‐diffusion due to preferential pathways in an aquifer has been observed in field data, challenging the assumption of Fickian dispersion at the local scale. This study develops a fully Lagrangian method to simulate sub‐grid super‐diffusion in a multidimensional regional‐scale transport model by using a recent mathematical model allowing super‐diffusion along the flow direction given by the regional model. Here, the time randomizing procedure known as subordination is applied to flow field output from MODFLOW simulations. Numerical tests check the applicability of the novel method in mapping regional‐scale super‐diffusive transport conditioned on local properties of multidimensional heterogeneous media. 相似文献
14.
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. 相似文献
15.
This article examines the required spatial discretization perpendicular to the fracture-matrix interface (FMI) for numerical simulation of solute transport in discretely fractured porous media. The discrete-fracture, finite-element model HydroGeoSphere ( Therrien et al. 2005 ) and a discrete-fracture implementation of MT3DMS ( Zheng 1990 ) were used to model solute transport in a single fracture, and the results were compared to the analytical solution of Tang et al. (1981) . To match analytical results on the relatively short timescales simulated in this study, very fine grid spacing perpendicular to the FMI of the scale of the fracture aperture is necessary if advection and/or dispersion in the fracture is high compared to diffusion in the matrix. The requirement of such extremely fine spatial discretization has not been previously reported in the literature. In cases of high matrix diffusion, matching the analytical results is achieved with larger grid spacing at the FMI. Cases where matrix diffusion is lower can employ a larger grid multiplier moving away from the FMI. The very fine spatial discretization identified in this study for cases of low matrix diffusion may limit the applicability of numerical discrete-fracture models in such cases. 相似文献
16.
Surface elevations represented in MODFLOW head-dependent packages are usually derived from digital elevation models (DEMs) that are available at much high resolution. Conventional grid refinement techniques to simulate the model at DEM resolution increases computational time, input file size, and in many cases are not feasible for regional applications. This research aims at utilizing the increasingly available high resolution DEMs for effective simulation of evapotranspiration (ET) in MODFLOW as an alternative to grid refinement techniques. The source code of the evapotranspiration package is modified by considering for a fixed MODFLOW grid resolution and for different DEM resolutions, the effect of variability in elevation data on ET estimates. Piezometric head at each DEM cell location is corrected by considering the gradient along row and column directions. Applicability of the research is tested for the lower Rio Grande (LRG) Basin in southern New Mexico. The DEM at 10 m resolution is aggregated to resampled DEM grid resolutions which are integer multiples of MODFLOW grid resolution. Cumulative outflows and ET rates are compared at different coarse resolution grids. Results of the analysis conclude that variability in depth-to-groundwater within the MODFLOW cell is a major contributing parameter to ET outflows in shallow groundwater regions. DEM aggregation methods for the LRG Basin have resulted in decreased volumetric outflow due to the formation of a smoothing error, which lowered the position of water table to a level below the extinction depth. 相似文献
17.
Hu BX 《Ground water》2006,44(2):222-233
A Lagrangian stochastic approach is applied to develop a method of moment for solute transport in a physically and chemically nonstationary medium. Stochastic governing equations for mean solute flux and solute covariance are analytically obtained in the first-order accuracy of log conductivity and/or chemical sorption variances and solved numerically using the finite-difference method. The developed method, the numerical method of moments (NMM), is used to predict radionuclide solute transport processes in the saturated zone below the Yucca Mountain project area. The mean, variance, and upper bound of the radionuclide mass flux through a control plane 5 km downstream of the footprint of the repository are calculated. According to their chemical sorption capacities, the various radionuclear chemicals are grouped as nonreactive, weakly sorbing, and strongly sorbing chemicals. The NMM method is used to study their transport processes and influence factors. To verify the method of moments, a Monte Carlo simulation is conducted for nonreactive chemical transport. Results indicate the results from the two methods are consistent, but the NMM method is computationally more efficient than the Monte Carlo method. This study adds to the ongoing debate in the literature on the effect of heterogeneity on solute transport prediction, especially on prediction uncertainty, by showing that the standard derivation of solute flux is larger than the mean solute flux even when the hydraulic conductivity within each geological layer is mild. This study provides a method that may become an efficient calculation tool for many environmental projects. 相似文献
18.
Incorporation of conceptual and parametric uncertainty into radionuclide flux estimates from a fractured granite rock mass 总被引：1，自引：1，他引：0
Donald M. Reeves Karl F. Pohlmann Greg M. Pohll Ming Ye Jenny B. Chapman 《Stochastic Environmental Research and Risk Assessment (SERRA)》2010,24(6):899-915
Detailed numerical flow and radionuclide simulations are used to predict the flux of radionuclides from three underground nuclear tests located in the Climax granite stock on the Nevada Test Site. The numerical modeling approach consists of both a regional-scale and local-scale flow model. The regional-scale model incorporates conceptual model uncertainty through the inclusion of five models of hydrostratigraphy and five models describing recharge processes for a total of 25 hydrostratigraphic–recharge combinations. Uncertainty from each of the 25 models is propagated to the local-scale model through constant head boundary conditions that transfer hydraulic gradients and flow patterns from each of the model alternatives in the vicinity of the Climax stock, a fluid flux calibration target, and model weights that describe the plausibility of each conceptual model. The local-scale model utilizes an upscaled discrete fracture network methodology where fluid flow and radionuclides are restricted to an interconnected network of fracture zones mapped onto a continuum grid. Standard Monte Carlo techniques are used to generate 200 random fracture zone networks for each of the 25 conceptual models for a total of 5,000 local-scale flow and transport realizations. Parameters of the fracture zone networks are based on statistical analysis of site-specific fracture data, with the exclusion of fracture density, which was calibrated to match the amount of fluid flux simulated through the Climax stock by the regional-scale models. Radionuclide transport is simulated according to a random walk particle method that tracks particle trajectories through the fracture continuum flow fields according to advection, dispersion and diffusional mass exchange between fractures and matrix. The breakthrough of a conservative radionuclide with a long half-life is used to evaluate the influence of conceptual and parametric uncertainty on radionuclide mass flux estimates. The fluid flux calibration target was found to correlate with fracture density, and particle breakthroughs were generally found to increase with increases in fracture density. Boundary conditions extrapolated from the regional-scale model exerted a secondary influence on radionuclide breakthrough for models with equal fracture density. The incorporation of weights into radionuclide flux estimates resulted in both noise about the original (unweighted) mass flux curves and decreases in the variance and expected value of radionuclide mass flux. 相似文献
19.
Characterizing heterogeneous soil water flow and solute transport using information measures 总被引：1，自引：0，他引：1
Heterogeneous water flow and solute transport in soils are an important phenomenon and difficult to be characterized. The objectives of this study were to investigate the heterogeneity of solute transport related to heterogeneous soil water flow using dye infiltration experiments, and to characterize heterogeneous water flow and solute transport in soils using the information theory. Field experiments of dye infiltration were performed in four plots. Various information measures were applied to characterize information content and complexity of water flow and solute transport in soils. Information contents and complexities of the maximum and apparent infiltration depths, and the mean and standard deviation of concentrations in the vertical direction of the plots were calculated. More heterogeneous processes of soil water flow and transport result in higher information/complexity values. The probability distributions of mean concentration were similar to those of the corresponding apparent infiltration depths for the plots, indicating that heterogeneity of dye concentrations was closely related to that of soil water flow. However, the range of information entropy and complexity of the water flow sequences was much narrower than that of the sequences of the concentrations. The results suggested that the transport processes were more heterogeneous than the water flow processes. Compared with the probability distributions of flow parameters, the information measures appeared to be a more versatile tool to describe flow and transport heterogeneities in soils. 相似文献
20.
Flow and transport simulation in karst aquifers remains a significant challenge for the ground water modeling community. Darcy's law–based models cannot simulate the inertial flows characteristic of many karst aquifers. Eddies in these flows can strongly affect solute transport. The simple two-region conduit/matrix paradigm is inadequate for many purposes because it considers only a capacitance rather than a physical domain. Relatively new lattice Boltzmann methods (LBMs) are capable of solving inertial flows and associated solute transport in geometrically complex domains involving karst conduits and heterogeneous matrix rock. LBMs for flow and transport in heterogeneous porous media, which are needed to make the models applicable to large-scale problems, are still under development. Here we explore aspects of these future LBMs, present simple examples illustrating some of the processes that can be simulated, and compare the results with available analytical solutions. Simulations are contrived to mimic simple capacitance-based two-region models involving conduit (mobile) and matrix (immobile) regions and are compared against the analytical solution. There is a high correlation between LBM simulations and the analytical solution for two different mobile region fractions. In more realistic conduit/matrix simulation, the breakthrough curve showed classic features and the two-region model fit slightly better than the advection-dispersion equation (ADE). An LBM-based anisotropic dispersion solver is applied to simulate breakthrough curves from a heterogeneous porous medium, which fit the ADE solution. Finally, breakthrough from a karst-like system consisting of a conduit with inertial regime flow in a heterogeneous aquifer is compared with the advection-dispersion and two-region analytical solutions. 相似文献