共查询到20条相似文献，搜索用时 296 毫秒
1.
《Advances in water resources》1999,23(3):261-269
A study of the effects of grid discretization on the migration of DNAPL within a discrete-fracture network embedded in a porous rock matrix is presented. It is shown that an insufficiently fine discretization of the fracture elements can lead to an overprediction of the volume of DNAPL that continues to migrate vertically at the intersection of a vertical and horizontal fracture. Uniform discretization of elements at the scale of one centimetre (or less) accurately resolved the density and capillary pressure components of the head gradient in the DNAPL. An alternative, non-uniform method of discretization of elements within the discrete-fracture network is presented whereby only fracture elements immediately adjacent to fracture intersections are refined. To further limit the number of elements employed, the porous matrix elements adjacent to the fracture elements are not similarly refined. Results show this alternative method of discretization reduces the numerical error to an acceptable level, while allowing the simulation of field-scale DNAPL contamination problems. The results from two field-scale simulations of a DNAPL-contaminated carbonate bedrock site in Ontario, Canada are presented. These simulations compare different methods of grid discretization, and highlight the importance of grid refinement when simulating DNAPL migration problems in fractured porous media. 相似文献
2.
A Eulerian analytical method is developed for nonreactive solute transport in heterogeneous, dual-permeability media where the hydraulic conductivities in fracture and matrix domains are both assumed to be stochastic processes. The analytical solution for the mean concentration is given explicitly in Fourier and Laplace transforms. Instead of using the fast fourier transform method to numerically invert the solution to real space (Hu et al., 2002), we apply the general relationship between spatial moments and concentration (Naff, 1990; Hu et al., 1997) to obtain the analytical solutions for the spatial moments up to the second for a pulse input of the solute. Owing to its accuracy and efficiency, the analytical method can be used to check the semi-analytical and Monte Carlo numerical methods before they are applied to more complicated studies. The analytical method can be also used during screening studies to identify the most significant transport parameters for further analysis. In this study, the analytical results have been compared with those obtained from the semi-analytical method (Hu et al., 2002) and the comparison shows that the semi-analytical method is robust. It is clearly shown from the analytical solution that the three factors, local dispersion, conductivity variation in each domain and velocity convection flow difference in the two domains, play different roles on the solute plume spreading in longitudinal and transverse directions. The calculation results also indicate that when the log-conductivity variance in matrix is 10 times less than its counterpart in fractures, it will hardly influence the solute transport, whether the conductivity field is matrix is treated as a homogeneous or random field. 相似文献
3.
Analytical solutions for the water flow and solute transport equations in the unsaturated zone are presented. We use the Broadbridge and White nonlinear model to solve the Richards’ equation for vertical flow under a constant infiltration rate. Then we extend the water flow solution and develop an exact parametric solution for the advection-dispersion equation. The method of characteristics is adopted to determine the location of a solute front in the unsaturated zone. The dispersion component is incorporated into the final solution using a singular perturbation method. The formulation of the analytical solutions is simple, and a complete solution is generated without resorting to computationally demanding numerical schemes. Indeed, the simple analytical solutions can be used as tools to verify the accuracy of numerical models of water flow and solute transport. Comparison with a finite-element numerical solution indicates that a good match for the predicted water content is achieved when the mesh grid is one-fourth the capillary length scale of the porous medium. However, when numerically solving the solute transport equation at this level of discretization, numerical dispersion and spatial oscillations were significant. 相似文献
4.
Discrete-fracture and dual-porosity models are infrequently used to simulate solute transport through fractured unconsolidated deposits, despite their more common application in fractured rock where distinct flow regimes are hypothesized. In this study, we apply four fracture transport models--the mobile-immobile model (MIM), parallel-plate discrete-fracture model (PDFM), and stochastic and deterministic discrete-fracture models (DFMs)--to demonstrate their utility for simulating solute transport through fractured till. Model results were compared to breakthrough curves (BTCs) for the conservative tracers potassium bromide (KBr), pentafluorobenzoic acid (PFBA), and 1,4-piperazinediethanesulfonic acid (PIPES) in a large-diameter column of fractured till. Input parameters were determined from independent field and laboratory methods. Predictions of Br BTCs were not significantly different among models; however, the stochastic and deterministic DFMs were more accurate than the MIM or PDFM when predicting PFBA and PIPES BTCs. DFMs may be more applicable than the MIM for tracers with small effective diffusion coefficients (De) or for short timescales due to differences in how these models simulate diffusion or incorporate heterogeneities by their fracture networks. At large scales of investigation, the more computationally efficient MIM and PDFM may be more practical to implement than the three-dimensional DFMs, or a combination of model approaches could be employed. Regardless of the modeling approach used, fractures should be incorporated routinely into solute transport models in glaciated terrain. 相似文献
5.
Vertical Discretization Impact in Numerical Modeling of Matrix Diffusion in Contaminated Groundwater
Shahla K. Farhat David T. Adamson Arun R. Gavaskar Sophia A. Lee Ronald W. Falta Charles J. Newell 《Ground Water Monitoring & Remediation》2020,40(2):52-64
Understanding the effects of contaminants that can diffuse into low-permeability (“low-k”) zones is crucial for effective groundwater remedial decision-making. Because low-k zones can serve as low-level sources of contamination to more transmissive zones over time, an accurate evaluation of the impacts of matrix diffusion at contaminated sites is vital. This study compared numerical groundwater flow and transport simulations using MODFLOW/RT3D at a hypothetical site using three cases, each with increasing discretization of the vertical 10-m thick domain: (1) a coarse multilayer heterogeneous grid based on one layer for each of four different hydrogeological units, (2) a “low-resolution” discretization approach where the low-k units were divided into several sublayers giving the model 10 layers, and (3) a “high-resolution” numerical model with 199 layers that are a few centimeters thick. When comparing the results of each case, significant differences were observed between the discretizations used, even though all other model input data were identical. The conventional grid models (Cases 1 and 2) appeared to underestimate groundwater plume concentrations by a factor ranging from 1.1 to 36 when compared to the high-resolution grid model (Case 3), and underestimated predicted cleanup times by more than a factor of 10 for some of the hypothetical sampling points in the modeling domain. These results validate the implication of Chapman et al. (2012), that conventional vertical discretization of numerical groundwater flow and transport models at contaminated sites (with layers that are greater than 1 m thick) can lead to significant errors when compared to more accurate high-resolution vertical discretization schemes (layers that are centimeters thick). 相似文献
6.
《Advances in water resources》2004,27(11):1045-1059
Transient and steady-state analytical solutions are derived to investigate solute transport in a fractured porous medium consisting of evenly spaced, parallel discrete fractures. The solutions incorporate a finite width strip source, longitudinal and transverse dispersion in the fractures, source decay, aqueous phase decay, one-dimensional diffusion into the matrix, sorption to fracture walls, and sorption within the matrix. The solutions are derived using Laplace and Fourier transforms, and inverted by interchanging the order of integration and utilizing a numerical Laplace inversion algorithm. The solutions are verified for simplified cases by comparison to solutions derived by Batu [Batu V. A generalized two-dimensional analytical solution for hydrodynamic dispersion in bounded media with the first-type condition at the source. Wat Resour Res 1989;25(6):1125] and Sudicky and Frind [Sudicky EA, Frind EO. Contaminant transport in fractured porous media: analytical solutions for a system of parallel fractures. Wat Resour Res 1982;18(6):1634]. The application of the solutions to a fractured sandstone demonstrates that narrower source widths and larger values of transverse dispersivity both lead to lower downstream concentrations in the fractures and shorter steady-state plumes. The incorporation of aqueous phase decay and source concentration decay both lead to lower concentrations and shorter plumes, with even moderate amounts of decay significantly shortening the persistence of contamination. 相似文献
7.
Boundary conditions are required to close the mathematical formulation of unstable density‐dependent flow systems. Proper implementation of boundary conditions, for both flow and transport equations, in numerical simulation are critical. In this paper, numerical simulations using the FEFLOW model are employed to study the influence of the different boundary conditions for unstable density‐dependent flow systems. A similar set up to the Elder problem is studied. It is well known that the numerical simulation results of the standard Elder problem are strongly dependent on spatial discretization. This work shows that for the cases where a solute mass flux boundary condition is employed instead of a specified concentration boundary condition at the solute source, the numerical simulation results do not vary between different convective solution modes (i.e., plume configurations) due to the spatial discretization. Also, the influence of various boundary condition types for nonsource boundaries was studied. It is shown that in addition to other factors such as spatial and temporal discretization, the forms of the solute transport equation such as divergent and convective forms as well as the type of boundary condition employed in the nonsource boundary conditions influence the convective solution mode in coarser meshes. On basis of the numerical experiments performed here, higher sensitivities regarding the numerical solution stability are observed for the Adams‐Bashford/Backward Trapezoidal time integration approach in comparison to the Euler‐Backward/Euler‐Forward time marching approach. The results of this study emphasize the significant consequences of boundary condition choice in the numerical modeling of unstable density‐dependent flow. 相似文献
8.
The numerical simulation of long‐term large‐scale (field to regional) variably saturated subsurface flow and transport remains a computational challenge, even with today's computing power. Therefore, it is appropriate to develop and use simplified models that focus on the main processes operating at the pertinent time and space scales, as long as the error introduced by the simpler model is small relative to the uncertainties associated with the spatial and temporal variation of boundary conditions and parameter values. This study investigates the effects of various model simplifications on the prediction of long‐term soil salinity and salt transport in irrigated soils. Average root‐zone salinity and cumulative annual drainage salt load were predicted for a 10‐year period using a one‐dimensional numerical flow and transport model (i.e. UNSATCHEM) that accounts for solute advection, dispersion and diffusion, and complex salt chemistry. The model uses daily values for rainfall, irrigation, and potential evapotranspiration rates. Model simulations consist of benchmark scenarios for different hypothetical cases that include shallow and deep water tables, different leaching fractions and soil gypsum content, and shallow groundwater salinity, with and without soil chemical reactions. These hypothetical benchmark simulations are compared with the results of various model simplifications that considered (i) annual average boundary conditions, (ii) coarser spatial discretization, and (iii) reducing the complexity of the salt‐soil reaction system. Based on the 10‐year simulation results, we conclude that salt transport modelling does not require daily boundary conditions, a fine spatial resolution, or complex salt chemistry. Instead, if the focus is on long‐term salinity, then a simplified modelling approach can be used, using annually averaged boundary conditions, a coarse spatial discretization, and inclusion of soil chemistry that only accounts for cation exchange and gypsum dissolution–precipitation. We also demonstrate that prediction errors due to these model simplifications may be small, when compared with effects of parameter uncertainty on model predictions. The proposed model simplifications lead to larger time steps and reduced computer simulation times by a factor of 1000. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
9.
The objective of this work is to develop a new numerical approach for the three-dimensional modelling of flow and transient solute transport in fractured porous media which would provide an accurate and efficient treatment of 3D complex geometries and inhomogeneities. For this reason, and in order to eliminate as much as possible the number of degrees of freedom, the fracture network, fractures and their intersections, are solved with a coupled 2D–1D model while the porous matrix is solved independently with a 3D model. The interaction between both models is accounted for by a coupling iterative technique. In this way it is possible to improve efficiency and reduce CPU usage by avoiding 3D mesh refinements of the fractures. The approach is based on the discrete-fracture model in which the exact geometry and location of each fracture in the network must be provided as an input. The formulation is based on a multidimensional coupling of the boundary element method-multidomain (BEM-MD) scheme for the flow and boundary element dual reciprocity method-multidomain (BE-DRM-MD) scheme for the transport. Accurate results and high efficiency have been obtained and are reported in this paper. 相似文献
10.
Three-dimensional grids representing a heterogeneous, ground water system are generated at 10 different resolutions in support of a site-scale flow and transport modeling effort. These grids represent hydrostratigraphy near Yucca Mountain, Nevada, consisting of 18 stratigraphic units with contrasting fluid flow and transport properties. The grid generation method allows the stratigraphy to be modeled by numerical grids of different resolution so that comparison studies can be performed to test for grid quality and determine the resolution required to resolve geologic structure and physical processes such as fluid flow and solute transport. The process of generating numerical grids with appropriate property distributions from geologic conceptual models is automated, thus making the entire process easy to implement with fewer user-induced errors. The series of grids of various resolutions are used to assess the level at which increasing resolution no longer influences the flow and solute transport results. Grid resolution is found to be a critical issue for ground water flow and solute transport. The resolution required in a particular instance is a function of the feature size of the model, the intrinsic properties of materials, the specific physics of the problem, and boundary conditions. The asymptotic nature of results related to flow and transport indicate that for a hydrologic model of the heterogeneous hydrostratigraphy under Yucca Mountain, a horizontal grid spacing of 600 m and vertical grid spacing of 40 m resolve the hydrostratigraphic model with sufficient precision to accurately model the hypothetical flow and solute transport to within 5% of the value that would be obtained with much higher resolution. 相似文献
11.
We present an analytical expression for the shear dispersion during solute transport in a coupled fracture–matrix system. The dispersion coefficient is obtained in a fracture with porous walls by taking into account an accurate boundary condition at the interface between the matrix and fracture, and the results were compared with those in a non-coupled system. The analysis presented identifies three regimes: diffusion-dominated, transition, and advection-dominated. The results showed that it is important to consider the exchange of solute between the fracture and matrix in development of the shear dispersion coefficient for the transition and advection-dominated regimes. The new dispersion coefficient is obtained by imposing the continuity of concentrations and mass fluxes along the porous walls. The resulting equivalent transport equation revealed that the effective velocity in a fracture increases while the dispersion coefficient decreases due to mass transfer between the matrix and fracture. A larger effective advection term leads to greater storage of mass in the matrix as compared with the classical double-porosity model with a non-coupled dispersion coefficient. The findings of this study can be used for modeling of tracer tests as well as fate, transport, and remediation of groundwater contaminants in fractured rocks. 相似文献
12.
The geochemical computer model PHREEQC can simulate solute transport in fractured bedrock aquifers that can be conceptualized as dual-porosity flow systems subject to one-dimensional advective-dispersive transport in the bedrock fractures and diffusive transport in the bedrock matrix. This article demonstrates how the physical characteristics of such flow systems can be parameterized for use in PHREEQC, it provides a method for minimizing numerical dispersion in PHREEQC simulations, and it compares PHREEQC simulations with results of an analytical solution. The simulations assumed a dual-porosity conceptual model involving advective-reactive-dispersive transport in the mobile zone (bedrock fracture) and diffusive-reactive transport in the immobile zone (bedrock matrix). The results from the PHREEQC dual-porosity transport model that uses a finite-difference approach showed excellent agreement compared with an analytical solution. 相似文献
13.
Various numerical methods have been used in the literature to simulate single and multiphase flow in fractured media. A promising approach is the use of the discrete-fracture model where the fracture entities in the permeable media are described explicitly in the computational grid. In this work, we present a critical review of the main conventional methods for multiphase flow in fractured media including the finite difference (FD), finite volume (FV), and finite element (FE) methods, that are coupled with the discrete-fracture model. All the conventional methods have inherent limitations in accuracy and applications. The FD method, for example, is restricted to horizontal and vertical fractures. The accuracy of the vertex-centered FV method depends on the size of the matrix gridcells next to the fractures; for an acceptable accuracy the matrix gridcells next to the fractures should be small. The FE method cannot describe properly the saturation discontinuity at the matrix–fracture interface. In this work, we introduce a new approach that is free from the limitations of the conventional methods. Our proposed approach is applicable in 2D and 3D unstructured griddings with low mesh orientation effect; it captures the saturation discontinuity from the contrast in capillary pressure between the rock matrix and fractures. The matrix–fracture and fracture–fracture fluxes are calculated based on powerful features of the mixed finite element (MFE) method which provides, in addition to the gridcell pressures, the pressures at the gridcell interfaces and can readily model the pressure discontinuities at impermeable faults in a simple way. To reduce the numerical dispersion, we use the discontinuous Galerkin (DG) method to approximate the saturation equation. We take advantage of a hybrid time scheme to alleviate the restrictions on the size of the time step in the fracture network. Several numerical examples in 2D and 3D demonstrate the robustness of the proposed model. Results show the significance of capillary pressure and orders of magnitude increase in computational speed compared to previous works. 相似文献
14.
Matrix diffusion can attenuate the rate of plume migration in fractured bedrock relative to the rate of ground water flow for both conservative and nonconservative solutes of interest. In a system of parallel, equally spaced constant aperture fractures subject to steady-state ground water flow and an infinite source width, the degree of plume attenuation increases with time and travel distance, eventually reaching an asymptotic level. The asymptotic degree of plume attenuation in the absence of degradation can be predicted by a plume attenuation factor, beta, which is readily estimated as R' (phi(m)/phi(f)), where R' is the retardation factor in the matrix, phi(m) is the matrix porosity, and phi(f) is the fracture porosity. This dual-porosity relationship can also be thought of as the ratio of primary to secondary porosity. Beta represents the rate of ground water flow in fractures relative to the rate of plume advance. For the conditions examined in this study, beta increases with greater matrix porosity, greater matrix fraction organic carbon, larger fracture spacing, and smaller fracture aperture. These concepts are illustrated using a case study where dense nonaqueous phase liquid in fractured sandstone produced a dissolved-phase trichloroethylene (TCE) plume approximately 300 m in length. Transport parameters such as matrix porosity, fracture porosity, hydraulic gradient, and the matrix retardation factor were characterized at the site through field investigations. In the fractured sandstone bedrock examined in this study, the asymptotic plume attenuation factors (beta values) for conservative and nonconservative solutes (i.e., chloride and TCE) were predicted to be approximately 800 and 12,210, respectively. Quantitative analyses demonstrate that a porous media (single-porosity) solute transport model is not appropriate for simulating contaminant transport in fractured sandstone where matrix diffusion occurs. Rather, simulations need to be conducted with either a discrete fracture model that explicitly incorporates matrix diffusion, or a dual-continuum model that accounts for mass transfer between mobile and immobile zones. Simulations also demonstrate that back diffusion from the matrix to fractures will likely be the time-limiting factor in reaching ground water cleanup goals in some fractured bedrock environments. 相似文献
15.
16.
Numerical simulations of variable-density flow and solute transport have been conducted to investigate dense plume migration for various configurations of 2D fracture networks. For orthogonal fractures, simulations demonstrate that dispersive mixing in fractures with small aperture does not stabilize vertical plume migration in fractures with large aperture. Simulations in non-orthogonal 2D fracture networks indicate that convection cells form and that they overlap both the porous matrix and fractures. Thus, transport rates in convection cells depend on matrix and fracture flow properties. A series of simulations in statistically equivalent networks of fractures with irregular orientation show that the migration of a dense plume is highly sensitive to the geometry of the network. If fractures in a random network are connected equidistantly to the solute source, few equidistantly distributed fractures favor density-driven transport. On the other hand, numerous fractures have a stabilizing effect, especially if diffusive transport rates are high. A sensitivity analysis for a network with few equidistantly distributed fractures shows that low fracture aperture, low matrix permeability and high matrix porosity impede density-driven transport because these parameters reduce groundwater flow velocities in both the matrix and the fractures. Enhanced molecular diffusion slows down density-driven transport because it favors solute diffusion from the fractures into the low-permeability porous matrix where groundwater velocities are smaller. For the configurations tested, variable-density flow and solute transport are most sensitive to the permeability and porosity of the matrix, which are properties that can be determined more accurately than the geometry and hydraulic properties of the fracture network, which have a smaller impact on density-driven transport. 相似文献
17.
Recognizing the spatial heterogeneity of hydraulic parameters, many researchers have studied the solute transport by both groundwater and channel flow in a stochastic framework. One of the methodologies used to up‐scale the stochastic solute transport equation, from a point‐location scale to a grid scale, is the cumulant expansion method combined with the calculus for the time‐ordered exponential and the calculus for the Lie operator. When the point‐location scale transport equation is scaled up to the grid scale, using the cumulant expansion method, a new dispersion coefficient emerges in the dispersive term of the solute transport equation in addition to the molecular dispersion coefficient. This velocity driven dispersion is called ‘macrodispersion’. The macrodispersion coefficient is the integral function of the time‐ordered covariance of the random velocity field. The integral is calculated over a Lagrangian trajectory of the flow. The Lagrangian trajectory depends on the following: (i) the spatial origin of the particle; (ii) the time when the macrodispersion is calculated; and (iii) the mean velocity field along the trajectory itself. The Lagrangian trajectory is a recursive function of time because the location of the particle along the trajectory at a particular time depends on the location of the particle at the previous time. This recursive functional form of the Lagrangian trajectory makes the calculation of the macrodispersion coefficient difficult. Especially for the unsteady, spatially non‐stationary, non‐uniform flow field, the macrodispersion coefficient is a highly complex expression and, so far, calculated using numerical methods in the discrete domains. Here, an analytical method was introduced to calculate the macrodispersion coefficient in the discrete domain for the unsteady and steady, spatially non‐stationary flow cases accurately and efficiently. This study can fill the gap between the theory of the ensemble averaged solute transport model and its numerical implementations. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
18.
Simulation of heat transport has its applications in geothermal exploitation of aquifers and the analysis of temperature dependent chemical reactions. Under homogeneous conditions and in the absence of a regional hydraulic gradient, groundwater flow and heat transport from or to a well exhibit radial symmetry, and governing equations are reduced by one dimension (1D) which increases computational efficiency importantly. Solute transport codes can simulate heat transport and input parameters may be modified such that the Cartesian geometry can handle radial flow. In this article, SEAWAT is evaluated as simulator for heat transport under radial flow conditions. The 1971, 1D analytical solution of Gelhar and Collins is used to compare axisymmetric transport with retardation (i.e., as a result of thermal equilibrium between fluid and solid) and a large diffusion (conduction). It is shown that an axisymmetric simulation compares well with a fully three dimensional (3D) simulation of an aquifer thermal energy storage systems. The influence of grid discretization, solver parameters, and advection solution is illustrated. Because of the high diffusion to simulate conduction, convergence criterion for heat transport must be set much smaller (10?10) than for solute transport (10?6). Grid discretization should be considered carefully, in particular the subdivision of the screen interval. On the other hand, different methods to calculate the pumping or injection rate distribution over different nodes of a multilayer well lead to small differences only. 相似文献
19.
F.A. Radu N. Suciu J. HoffmannA. Vogel O. Kolditz C.-H. ParkS. Attinger 《Advances in water resources》2011,34(1):47-61
This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater. 相似文献
20.
In general, the accuracy of numerical simulations is determined by spatial and temporal discretization levels. In fractured porous media, the time step size is a key factor in controlling the solution accuracy for a given spatial discretization. If the time step size is restricted by the relatively rapid responses in the fracture domain to maintain an acceptable level of accuracy in the entire simulation domain, the matrix tends to be temporally over-discretized. Implicit sub-time stepping applies smaller sub-time steps only to the sub-domain where the accuracy requirements are less tolerant and is most suitable for problems where the response is high in only a small portion of the domain, such as within and near the fractures in fractured porous media. It is demonstrated with illustrative examples that implicit sub-time stepping can significantly improve the simulation efficiency with minimal loss in accuracy when simulating flow and transport in fractured porous media. The methodology is successfully applied to density-dependent flow and transport simulations in a Canadian Shield environment, where the flow and transport is dominated by discrete, highly conductive fracture zones. 相似文献