共查询到20条相似文献，搜索用时 31 毫秒
1.
Writing Analytic Element Programs in Python 总被引：1，自引：0，他引：1
The analytic element method is a mesh-free approach for modeling ground water flow at both the local and the regional scale. With the advent of the Python object-oriented programming language, it has become relatively easy to write analytic element programs. In this article, an introduction is given of the basic principles of the analytic element method and of the Python programming language. A simple, yet flexible, object-oriented design is presented for analytic element codes using multiple inheritance. New types of analytic elements may be added without the need for any changes in the existing part of the code. The presented code may be used to model flow to wells (with either a specified discharge or drawdown) and streams (with a specified head). The code may be extended by any hydrogeologist with a healthy appetite for writing computer code to solve more complicated ground water flow problems. 相似文献
2.
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. 相似文献
3.
Canonical balanced dynamic equations involving vertically sheared horizontal flow with heat or mass sources have emerged recently in systematic multi-scale modeling of the equatorial wave guide on a wide range of spatio-temporal scales. Here, a new self-contained derivation of these equations is developed briefly in a context for potential applications to the hurricane embryo. These canonical balanced equations are studied through a combination of exact solutions and simple numerics. The results below include elementary exact solutions given by velocity fields that are linear in the spatial coordinates combined with an exact nonlinear stability analysis for vertical vorticity amplification in such a preconditioned environment. Other elementary solutions studied here include the evolution of radial eddies, which represent “hot towers” in the hurricane embryo in a suitable radial preconditioned background environment. 相似文献
4.
The methods behind the predefined impulse response function in continuous time (PIRFICT) time series model are extended to cover more complex situations where multiple stresses influence ground water head fluctuations simultaneously. In comparison to autoregressive moving average (ARMA) time series models, the PIRFICT model is optimized for use on hydrologic problems. The objective of the paper is twofold. First, an approach is presented for handling multiple stresses in the model. Each stress has a specific parametric impulse response function. Appropriate impulse response functions for other stresses than precipitation are derived from analytical solutions of elementary hydrogeological problems. Furthermore, different stresses do not need to be connected in parallel in the model, as is the standard procedure in ARMA models. Second, general procedures are presented for modeling and interpretation of the results. The multiple-input PIRFICT model is applied to two real cases. In the first one, it is shown that this model can effectively decompose series of ground water head fluctuations into partial series, each representing the influence of an individual stress. The second application handles multiple observation wells. It is shown that elementary physical knowledge and the spatial coherence in the results of multiple wells in an area may be used to interpret and check the plausibility of the results. The methods presented can be used regardless of the hydrogeological setting. They are implemented in a computer package named Menyanthes (www.menyanthes.nl). 相似文献
5.
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. 相似文献
6.
Discharges of polluted water from abandoned mines are a major cause of degradation of water resources worldwide. Pollution arises after abandoned workings flood up to surface level, by the process termed ground water rebound. As flow in large, open mine voids is often turbulent, standard techniques for modeling ground water flow (which assume laminar flow) are inappropriate for predicting ground water rebound. More physically realistic models are therefore desirable, yet these are often expensive to apply to all but the smallest of systems. An overall strategy for ground water rebound modeling is proposed, with models of decreasing complexity applied as the temporal and spatial scales of the systems under analysis increase. For relatively modest systems (area < 200 km2), a physically based modeling approach has been developed, in which 3-D pipe networks (representing major mine roadways, etc.) are routed through a variably saturated, 3-D porous medium (representing the country rock). For systems extending more than 100 to 3000 km2, a semidistributed model (GRAM) has been developed, which conceptualizes extensively interconnected volumes of workings as ponds, which are connected to other ponds only at discrete overflow points, such as major inter-mine roadways, through which flow can be efficiently modeled using the Prandtl-Nikuradse pipe-flow formulation. At the very largest scales, simple water-balance calculations are probably as useful as any other approach, and a variety of proprietary codes may be used for the purpose. 相似文献
7.
Robert G. THOMAS 《水文科学杂志》2013,58(4):50-69
Abstract Graphical solutions to linear and radial ground water flow problems are described. These methods are based on finite difference approximations which have been widely used in heat flow problems. The case of the leaky aquiclude and vertical recharge to unconfined aquifers is discussed as well as constant head at a control and constant discharge or recharge. These solutions are useful only when the physical factors are known and an estimate of water level changes under certain conditions is desired. The methods are not intended to replace more rigorous mathematical solutions. 相似文献
8.
Multilayer analytic element modeling of radial collector wells 总被引：1，自引：0，他引：1
A new multilayer approach is presented for the modeling of ground water flow to radial collector wells. The approach allows for the inclusion of all aspects of the unique boundary condition along the lateral arms of a collector well, including skin effect and internal friction losses due to flow in the arms. The hydraulic conductivity may differ between horizontal layers within the aquifer, and vertical anisotropy can be taken into account. The approach is based on the multilayer analytic element method, such that regional flow and local three-dimensional detail may be simulated simultaneously and accurately within one regional model. Horizontal flow inside a layer is computed analytically, while vertical flow is approximated with a standard finite-difference scheme. Results obtained with the proposed approach compare well to results obtained with three-dimensional analytic element solutions for flow in unconfined aquifers. The presented approach may be applied to predict the yield of a collector well in a regional setting and to compute the origin and residence time, and thus the quality, of water pumped by the collector well. As an example, the addition of three lateral arms to a collector well that already has three laterals is investigated. The new arms are added at an elevation of 2 m above the existing laterals. The yield increase of the collector well is computed as a function of the lengths of the three new arms. 相似文献
9.
Calculation of ground water ages--a comparative analysis 总被引：1，自引：0，他引：1
Ground water age is a fundamental, yet complex, concept in ground water hydrology. Discrepancies between results obtained through different modeling approaches for ground water age calculation have been reported, in particular, between ground water ages modeled by advection and direct simulation of ground water ages (e.g., age-mass approach), which includes effects of advection and dispersion. Here, through a series of two-dimensional (2D) simulations, the impact of water mixing through advection and dispersion on modeled 14C and directly simulated ground water ages is assessed. Impact of dispersion on modeled ages is systematically stronger in areas where water velocities are smaller and far more pronounced on 14C ages. This effect is also observed in one-dimensional models. 2D simulations show that longitudinal dispersion generally acts as a "source" of 14C, while vertical dispersion acts as a "sink," leading to apparent younger or older modeled 14C ages as compared to advective and directly simulated ground water ages. The presence of permeable and impermeable faults provides an equally important source for discrepancies, leading to major differences in modeled ages among the three methods considered. Overall, our results show that a 14C modeling approach using a solute transport model for calculating ground water age appears to be more reliable in ground water systems without faults and where water velocities are relatively high than in systems that are relatively more heterogeneous and those where faults are present. Among the three modeling approaches considered here, direct simulation of ground water age seems to yield the most consistent results in complex, heterogeneous ground water flow systems, giving a vertical age structure consistent with ages expected from consideration of the flow system. 相似文献
10.
Comparative study of methods for WHPA delineation 总被引：3，自引：0，他引：3
Paradis D Martel R Karanta G Lefebvre R Michaud Y Therrien R Nastev M 《Ground water》2007,45(2):158-167
Human activities, whether agricultural, industrial, commercial, or domestic, can contribute to ground water quality deterioration. In order to protect the ground water exploited by a production well, it is essential to develop a good knowledge of the flow system and to adequately delineate the area surrounding the well within which potential contamination sources should be managed. Many methods have been developed to delineate such a wellhead protection area (WHPA). The integration of more information on the geologic and hydrogeologic characteristics of the study area increases the precision of any given WHPA delineation method. From a practical point of view, the WHPA delineation methods allowing the simplest and least expensive integration of the available information should be favored. This paper presents a comparative study in which nine different WHPA delineation methods were applied to a well and a spring in an unconfined granular aquifer and to a well in a confined highly fractured rock aquifer. These methods range from simple approaches to complex computer models. Hydrogeological mapping and numerical modeling with MODFLOW-MODPATH were used as reference methods to respectively compare the delineation of the zone of contribution and the zone of travel obtained from the various WHPA methods. Although applied to simple ground water flow systems, these methods provided a relatively wide range of results. To allow a realistic delineation of the WHPA in aquifers of variable geometry, a WHPA delineation method should ensure a water balance and include observed or calculated regional flow characteristics. 相似文献
11.
Ground water recharge is often estimated through the calibration of ground water flow models. We examine the nature of calibration errors by considering some simple mathematical and numerical calculations. From these calculations, we conclude that calibrating a steady-state ground water flow model to water level extremes yields estimates of recharge that have the same value as the time-varying recharge at the time the water levels are measured. These recharge values, however, are a subdued version of the actual transient recharge signal. In addition, calibrating a steady-state ground water flow model to data collected during periods of rising water levels will produce recharge values that underestimate the actual transient recharge. Similarly, calibrating during periods of falling water levels will overestimate the actual transient recharge. We also demonstrate that average water levels can be used to estimate the actual average recharge rate provided that water level data have been collected for a sufficient amount of time. 相似文献
12.
Olsthoorn TN 《Ground water》2008,46(1):13-22
Convolution is a form of superposition that efficiently deals with input varying arbitrarily in time or space. It works whenever superposition is applicable, that is, for linear systems. Even though convolution is well-known since the 19th century, this valuable method is still missing in most textbooks on ground water hydrology. This limits widespread application in this field. Perhaps most papers are too complex mathematically as they tend to focus on the derivation of analytical expressions rather than solving practical problems. However, convolution is straightforward with standard mathematical software or even a spreadsheet, as is demonstrated in the paper. The necessary system responses are not limited to analytic solutions; they may also be obtained by running an already existing ground water model for a single stress period until equilibrium is reached. With these responses, high-resolution time series of head or discharge may then be computed by convolution for arbitrary points and arbitrarily varying input, without further use of the model. There are probably thousands of applications in the field of ground water hydrology that may benefit from convolution. Therefore, its inclusion in ground water textbooks and courses is strongly needed. 相似文献
13.
As part of its Wellhead Protection Program, the U.S. EPA mandates the delineation of "time-of-travel capture zones" as the basis for the definition of wellhead protection zones surrounding drinking water production wells. Depending on circumstances the capture zones may be determined using methods that range from simply drawing a circle around the well to sophisticated ground water flow and transport modeling. The simpler methods are attractive when faced with the delineation of hundreds or thousands of capture zones for small public drinking water supply wells. On the other hand, a circular capture zone may not be adequate in the presence of an ambient ground water flow regime. A dimensionless time-of-travel parameter T is used to determine when calculated fixed-radius capture zones can be used for drinking water production wells. The parameter incorporates aquifer properties, the magnitude of the ambient ground water flow field, and the travel time criterion for the time-of-travel capture zone. In the absence of interfering flow features, three different simple capture zones can be used depending on the value of T . A modified calculated fixed-radius capture zone proves protective when T < 0.1, while a more elongated capture zone must be used when T > 1. For values of T between 0.1 and 1, a circular capture zone can be used that is eccentric with respect to the well. Finally, calculating T allows for a quick assessment of the validity of circular capture zones without redoing the delineation with a computer model. 相似文献
14.
Laboratory and numerical investigation of flow and transport near a seepage-face boundary 总被引：3，自引：0，他引：3
Laboratory and numerical modeling investigations were completed to study the unconfined ground water flow and transport processes near a seepage-face boundary. The laboratory observations were made in a radial sand tank and included measurements of the height of the seepage face, flow velocity near the seepage face, travel time distribution of multiple tracer slugs, and streamlines. All the observations were reliably reproduced with a three-dimensional, axi-symmetric, variably saturated ground water flow model. Physical data presented in this work demonstrate and quantify the importance of three-dimensional transport patterns within a seepage-face zone. The results imply that vertically averaged flow models that employ Dupuit approximations might introduce error in the analysis of localized solute transport near a seepage-face boundary. The experimental dataset reported in this work will also be of interest for those who are attempting to validate a numerical algorithm for solving ground water and contaminant discharge patterns near a surface-water boundary. 相似文献
15.
The availability of powerful desktop computers and graphical user interfaces for ground water flow models makes possible the construction of ever more complex models. A proposed copper-zinc sulfide mine in northern Wisconsin offers a unique case in which the same hydrologic system has been modeled using a variety of techniques covering a wide range of sophistication and complexity. Early in the permitting process, simple numerical models were used to evaluate the necessary amount of water to be pumped from the mine, reductions in streamflow, and the drawdowns in the regional aquifer. More complex models have subsequently been used in an attempt to refine the predictions. Even after so much modeling effort, questions regarding the accuracy and reliability of the predictions remain. We have performed a new analysis of the proposed mine using the two-dimensional analytic element code GFLOW coupled with the nonlinear parameter estimation code UCODE. The new model is parsimonious, containing fewer than 10 parameters, and covers a region several times larger in areal extent than any of the previous models. The model demonstrates the suitability of analytic element codes for use with parameter estimation codes. The simplified model results are similar to the more complex models; predicted mine inflows and UCODE-derived 95% confidence intervals are consistent with the previous predictions. More important, the large areal extent of the model allowed us to examine hydrological features not included in the previous models, resulting in new insights about the effects that far-field boundary conditions can have on near-field model calibration and parameterization. In this case, the addition of surface water runoff into a lake in the headwaters of a stream while holding recharge constant moved a regional ground watershed divide and resulted in some of the added water being captured by the adjoining basin. Finally, a simple analytical solution was used to clarify the GFLOW model's prediction that, for a model that is properly calibrated for heads, regional drawdowns are relatively unaffected by the choice of aquifer properties, but that mine inflows are strongly affected. Paradoxically, by reducing model complexity, we have increased the understanding gained from the modeling effort. 相似文献
16.
17.
Mehl S 《Ground water》2006,44(4):583-594
This study examines the use of Picard and Newton iteration to solve the nonlinear, saturated ground water flow equation. Here, a simple three-node problem is used to demonstrate the convergence difficulties that can arise when solving the nonlinear, saturated ground water flow equation in both homogeneous and heterogeneous systems with and without nonlinear boundary conditions. For these cases, the characteristic types of convergence patterns are examined. Viewing these convergence patterns as orbits of an attractor in a dynamical system provides further insight. It is shown that the nonlinearity that arises from nonlinear head-dependent boundary conditions can cause more convergence difficulties than the nonlinearity that arises from flow in an unconfined aquifer. Furthermore, the effects of damping on both convergence and convergence rate are investigated. It is shown that no single strategy is effective for all problems and how understanding pitfalls and merits of several methods can be helpful in overcoming convergence difficulties. Results show that Picard iterations can be a simple and effective method for the solution of nonlinear, saturated ground water flow problems. 相似文献
18.
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. 相似文献
19.
20.
Hunt RJ 《Ground water》2006,44(1):5-15
Though powerful and easy to use, applications of the analytic element method are not as widespread as finite-difference or finite-element models due in part to their relative youth. Although reviews that focus primarily on the mathematical development of the method have appeared in the literature, a systematic review of applications of the method is not available. An overview of the general types of applications of analytic elements in ground water modeling is provided in this paper. While not fully encompassing, the applications described here cover areas where the method has been historically applied (regional, two-dimensional steady-state models, analyses of ground water-surface water interaction, quick analyses and screening models, wellhead protection studies) as well as more recent applications (grid sensitivity analyses, estimating effective conductivity and dispersion in highly heterogeneous systems). The review of applications also illustrates areas where more method development is needed (three-dimensional and transient simulations). 相似文献