共查询到20条相似文献，搜索用时 572 毫秒
1.
Chuiyu Lu Bo Zhang Xin He Guoliang Cao Qingyan Sun Lingjia Yan Tao Qin Tianchen Li Zepeng Li 《Ground water》2021,59(1):90-99
MODFLOW is one of the most popular groundwater simulation tools available; however, the development of lake modules that can be coupled with MODFLOW is lacking apart from the LAK3 package. This study proposes a new approach for simulating lake - groundwater interaction under steady-state flow, referred to as the sloping lakebed method (SLM). In this new approach, discretization of the lakebed in the vertical direction is independent of the spatial discretization of the aquifer system, which can potentially solve the problem that the lake and groundwater are usually simulated at different scales. The lakebed is generalized by a slant at the bottom of each lake grid cell, which can be classified as fully submerged, dry, and partly submerged. The SLM method accounts for all lake sources and sinks, establishing a governing equation that can be solved using Newton's method. A benchmarking case study was conducted using a modified model setup in the LAK3 user manual. It was found that when there is a sufficient number of layers at the top of the groundwater model, SLM simulates an almost identical groundwater head as the LAK3-based model; when the number of layers decreases, SLM is unaffected while LAK3 may be at a risk of giving unrealistic results. Additionally, the SLM can reflect the relationship between the simulated lake surface area and lake water depth more accurately. Therefore, the SLM method is a promising alternative to the LAK3 package when simulating lake - groundwater interaction. 相似文献
2.
Regional finite‐difference models often have cell sizes that are too large to sufficiently model well‐stream interactions. Here, a steady‐state hybrid model is applied whereby the upper layer or layers of a coarse MODFLOW model are replaced by the analytic element model GFLOW, which represents surface waters and wells as line and point sinks. The two models are coupled by transferring cell‐by‐cell leakage obtained from the original MODFLOW model to the bottom of the GFLOW model. A real‐world test of the hybrid model approach is applied on a subdomain of an existing model of the Lake Michigan Basin. The original (coarse) MODFLOW model consists of six layers, the top four of which are aggregated into GFLOW as a single layer, while the bottom two layers remain part of MODFLOW in the hybrid model. The hybrid model and a refined “benchmark” MODFLOW model simulate similar baseflows. The hybrid and benchmark models also simulate similar baseflow reductions due to nearby pumping when the well is located within the layers represented by GFLOW. However, the benchmark model requires refinement of the model grid in the local area of interest, while the hybrid approach uses a gridless top layer and is thus unaffected by grid discretization errors. The hybrid approach is well suited to facilitate cost‐effective retrofitting of existing coarse grid MODFLOW models commonly used for regional studies because it leverages the strengths of both finite‐difference and analytic element methods for predictions in mildly heterogeneous systems that can be simulated with steady‐state conditions. 相似文献
3.
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. 相似文献
4.
Regional finite-difference models tend to have large cell sizes, often on the order of 1–2 km on a side. Although the regional flow patterns in deeper formations may be adequately represented by such a model, the intricate surface water and groundwater interactions in the shallower layers are not. Several stream reaches and nearby wells may occur in a single cell, precluding any meaningful modeling of the surface water and groundwater interactions between the individual features. We propose to replace the upper MODFLOW layer or layers, in which the surface water and groundwater interactions occur, by an analytic element model (GFLOW) that does not employ a model grid; instead, it represents wells and surface waters directly by the use of point-sinks and line-sinks. For many practical cases it suffices to provide GFLOW with the vertical leakage rates calculated in the original coarse MODFLOW model in order to obtain a good representation of surface water and groundwater interactions. However, when the combined transmissivities in the deeper (MODFLOW) layers dominate, the accuracy of the GFLOW solution diminishes. For those cases, an iterative coupling procedure, whereby the leakages between the GFLOW and MODFLOW model are updated, appreciably improves the overall solution, albeit at considerable computational cost. The coupled GFLOW–MODFLOW model is applicable to relatively large areas, in many cases to the entire model domain, thus forming an attractive alternative to local grid refinement or inset models. 相似文献
5.
The nonhorizontal‐model‐layer (NHML) grid system is more accurate than the horizontal‐model‐layer grid system to describe groundwater flow in an unconfined sloping aquifer on the basis of MODFLOW‐2000. However, the finite‐difference scheme of NHML was based on the Dupuit‐Forchheimer assumption that the streamlines were horizontal, which was acceptable for slope less than 0.10. In this study, we presented a new finite‐difference scheme of NHML based on the Boussinesq assumption and developed a new package SLOPE which was incorporated into MODFLOW‐2000 to become the MODFLOW‐SP model. The accuracy of MODFLOW‐SP was tested against solution of Mac Cormack (1969). The differences between the solutions of MODFLOW‐2000 and MODFLOW‐SP were nearly negligible when the slope was less than 0.27, and they were noticeable during the transient flow stage and vanished in steady state when the slope increased above 0.27. We established a model considering the vertical flow using COMSOL Multiphysics to test the robustness of constrains used in MODFLOW‐SP. The results showed that streamlines quickly became parallel with the aquifer base except in the narrow regions near the boundaries when the initial flow was not parallel to the aquifer base. MODFLOW‐SP can be used to predict the hydraulic head of an unconfined aquifer along the profile perpendicular to the aquifer base when the slope was smaller than 0.50. The errors associated with constrains used in MODFLOW‐SP were small but noticeable when the slope increased to 0.75, and became significant for the slope of 1.0. 相似文献
6.
Sahila Beegum Jiří Šimůnek Adam Szymkiewicz K. P. Sudheer Indumathi M. Nambi 《Ground water》2019,57(3):392-408
The “HYDRUS package for MODFLOW” is an existing MODFLOW package that allows MODFLOW to simultaneously evaluate transient water flow in both unsaturated and saturated zones. The package is based on incorporating parts of the HYDRUS-1D model (to simulate unsaturated water flow in the vadose zone) into MODFLOW (to simulate saturated groundwater flow). The coupled model is effective in addressing spatially variable saturated-unsaturated hydrological processes at the regional scale. However, one of the major limitations of this coupled model is that it does not have the capability to simulate solute transport along with water flow and therefore, the model cannot be employed for evaluating groundwater contamination. In this work, a modified unsaturated flow and transport package (modified HYDRUS package for MODFLOW and MT3DMS) has been developed and linked to the three-dimensional (3D) groundwater flow model MODFLOW and the 3D groundwater solute transport model MT3DMS. The new package can simulate, in addition to water flow in the vadose zone, also solute transport involving many biogeochemical processes and reactions, including first-order degradation, volatilization, linear or nonlinear sorption, one-site kinetic sorption, two-site sorption, and two-kinetic sites sorption. Due to complex interactions at the groundwater table, certain modifications of the pressure head (compared to the original coupling) and solute concentration profiles were incorporated into the modified HYDRUS package. The performance of the newly developed model is evaluated using HYDRUS (2D/3D), and the results indicate that the new model is effective in simulating the movement of water and contaminants in the saturated-unsaturated flow domains. 相似文献
7.
We present a contaminant treatment system (CTS) package for MODFLOW 6 that facilitates the simulation of pump-and-treat systems for groundwater remediation. Using the “nonintrusive” MODFLOW 6 application programming interface (API) capability, the CTS package can balance flows between extraction and injection wells within the outer flow solution loop and applies blended concentration/mass treatment efficiency within the outer transport solution loop. The former can be important when the requested extraction rate cannot be satisfied by the current simulated groundwater system conditions, while the latter can be important for simulating incomplete/imperfect treatment schemes. Furthermore, the CTS package allows users to temporally vary all aspects of a simulated CTS system, including the configuration and location of injection and extraction wells, and the CTS efficiency. This flexibility combined with the API-based implementation provide a generic and general CTS package that can be applied across the wide range of MODFLOW 6 simulation options and that evolves in step with MODFLOW 6 code modifications and advancements without needing to update the CTS package itself. 相似文献
8.
The ground water flow model MODFLOW inherently implements a nongeneralized integrated finite-difference (IFD) numerical scheme. The IFD numerical scheme allows for construction of finite-difference model grids with curvilinear (piecewise linear) rows. The resulting grid comprises model cells in the shape of trapezoids and is distorted in comparison to a traditional MODFLOW finite-difference grid. A version of MODFLOW-88 (herein referred to as MODFLOW IFD) with the code adapted to make the one-dimensional DELR and DELC arrays two dimensional, so that equivalent conductance between distorted grid cells can be calculated, is described. MODFLOW IFD is used to inspect the sensitivity of the numerical head and velocity solutions to the level of distortion in trapezoidal grid cells within a converging radial flow domain. A test problem designed for the analysis implements a grid oriented such that flow is parallel to columns with converging widths. The sensitivity analysis demonstrates MODFLOW IFD's capacity to numerically derive a head solution and resulting intercell volumetric flow when the internal calculation of equivalent conductance accounts for the distortion of the grid cells. The sensitivity of the velocity solution to grid cell distortion indicates criteria for distorted grid design. In the radial flow test problem described, the numerical head solution is not sensitive to grid cell distortion. The accuracy of the velocity solution is sensitive to cell distortion with error <1% if the angle between the nonparallel sides of trapezoidal cells is <12.5 degrees. The error of the velocity solution is related to the degree to which the spatial discretization of a curve is approximated with piecewise linear segments. Curvilinear finite-difference grid construction adds versatility to spatial discretization of the flow domain. MODFLOW-88's inherent IFD numerical scheme and the test problem results imply that more recent versions of MODFLOW 2000, with minor modifications, have the potential to make use of a curvilinear grid. 相似文献
9.
Complex stratigraphy can be difficult to simulate in MODFLOW models. MODFLOW uses a structured grid that requires that each grid layer be continuous throughout the model domain. This makes it difficult to explicitly represent common features such as pinchouts and embedded seams in a MODFLOW model. In this paper, we describe a method for automatically generating MODFLOW-compatible grids from boundary-representation solid models. Solid models are data structures developed originally for computer-aided design applications that define the geometry of three-dimensional objects. Solid models can be used to represent arbitrarily complex stratigraphy. The elevations defined by the solids are then extracted from the solids in a manner that preserves the continuous-layer requirement imposed by MODFLOW. Two basic approaches are described: The first method adjusts the MODFLOW grid dimensions (layer elevations) to fit the solid model boundaries, and the second method creates a regular MODFLOW grid and adjusts the material properties to match the changes in stratigraphy. One of the main benefits of using solid models to define stratigraphy for MODFLOW models is that it provides a grid-independent definition of the layer elevations that can be used to immediately re-create the MODFLOW grid geometry after any change to the grid resolution. 相似文献
10.
A new MODFLOW package (Nonlinear Flow Process; NLFP) simulating nonlinear flow following the Forchheimer equation was developed and implemented in MODLFOW‐2005. The method is based on an iterative modification of the conductance calculated and used by MODFLOW to obtain an effective Forchheimer conductance. The package is compatible with the different layer types, boundary conditions, and solvers as well as the wetting capability of MODFLOW. The correct implementation is demonstrated using four different benchmark scenarios for which analytical solutions are available. A scenario considering transient flow in a more realistic setting and a larger model domain with a higher number of cells demonstrates that NLFP performs well under more complex conditions, although it converges moderately slower than the standard MODFLOW depending on the nonlinearity of flow. Thus, this new tool opens a field of opportunities to groundwater flow simulation with MODFLOW, especially for core sample simulation or vuggy karstified aquifers as well as for nonlinear flow in the vicinity of pumping wells. 相似文献
11.
In a typical ground water flow model, lakes are represented by specified head nodes requiring that lake levels be known a priori. To remove this limitation, previous researchers assigned high hydraulic conductivity (K) values to nodes that represent a lake, under the assumption that the simulated head at the nodes in the high-K zone accurately reflects lake level. The solution should also produce a constant water level across the lake. We developed a model of a simple hypothetical ground water/lake system to test whether solutions using high-K lake nodes are sensitive to the value of K selected to represent the lake. Results show that the larger the contrast between the K of the aquifer and the K of the lake nodes, the smaller the error tolerance required for the solution to converge. For our test problem, a contrast of three orders of magnitude produced a head difference across the lake of 0.005 m under a regional gradient of the order of 10(-3) m/m, while a contrast of four orders of magnitude produced a head difference of 0.001 m. The high-K method was then used to simulate lake levels in Pretty Lake, Wisconsin. Results for both the hypothetical system and the application to Pretty Lake compared favorably with results using a lake package developed for MODFLOW (Merritt and Konikow 2000). While our results demonstrate that the high-K method accurately simulates lake levels, this method has more cumbersome postprocessing and longer run times than the same problem simulated using the lake package. 相似文献
12.
《水文科学杂志》2013,58(6):953-970
Abstract The 5000 km2 topographically closed Estancia basin in central New Mexico has been the focus of several palaeoclimatic studies based on changes in the level of late Pleistocene Lake Estancia. A large, unknown volume of surface runoff and groundwater from adjacent mountains contributed to the hydrological balance during highstands and lowstands. The US Department of Agriculture hydrological model SWAT (Soil and Water Assessment Tool) and the US Geological Survey groundwater flow model MODFLOW, with the LAK2 package, were used in this study to estimate runoff and water balance under present climate. A Geographic Information Systems (GIS) interface was used for SWAT, digitized data were applied for soils and vegetation, and limited streamflow data were used to obtain an approximate calibration for the model. Simulated streamflow is generally within 30% of observed values, and simulated runoff for the entire basin is about 8% of the annual inflow volume needed to support lowstands of the former Lake Estancia. Results from the combined models suggest application to other palaeoclimate investigations in semiarid lake basins. 相似文献
13.
SFRmaker and Linesink-Maker: Rapid Construction of Streamflow Routing Networks from Hydrography Data
Groundwater models have evolved to encompass more aspects of the water cycle, but the incorporation of realistic boundary conditions representing surface water remains time-consuming and error-prone. We present two Python packages that robustly automate this process using readily available hydrography data as the primary input. SFRmaker creates input for the MODFLOW SFR package, while Linesink-maker creates linesink string input for the GFLOW analytic element program. These programs can reduce weeks or even months of manual effort to a few minutes of execution time, and carry the added advantages of reduced potential for error, improved reproducibility and facilitation of step-wise modeling through reduced dependency on a particular conceptual model or discretization. Two real-world examples at the county to multi-state scales are presented. 相似文献
14.
Joshua D. Larsen Christian D. Langevin Joseph D. Hughes Richard G. Niswonger 《Ground water》2024,62(1):157-166
An agricultural water use package has been developed for MODFLOW 6 using the MODFLOW Application Programming Interface (API). The MODFLOW API Agricultural Water Use Package (API-Ag) was based on the approach to simulate irrigation demand in the MODFLOW-NWT and GSFLOW Agricultural Water Use (AG) Package. The API-Ag Package differs from the previous approach by implementing new features and support for additional irrigation providers. New features include representation of deficit and over-irrigation, Multi-Aquifer Well and Lake Package irrigation providers, and support for structured, vertex, and unstructured grid models. Three example problems are presented that demonstrate how the API-Ag Package improves representation of highly managed systems and are further used to validate the irrigation demand and delivery formulations. Irrigation volumes simulated in the three example problems show excellent agreement with the MODFLOW-NWT AG Package. 相似文献
15.
Fast transport simulation with an adaptive grid refinement 总被引：2，自引：0，他引：2
One of the main difficulties in transport modeling and calibration is the extraordinarily long computing times necessary for simulation runs. Improved execution time is a prerequisite for calibration in transport modeling. In this paper we investigate the problem of code acceleration using an adaptive grid refinement, neglecting subdomains, and devising a method by which the Courant condition can be ignored while maintaining accurate solutions. Grid refinement is based on dividing selected cells into regular subcells and including the balance equations of subcells in the equation system. The connection of coarse and refined cells satisfies the mass balance with an interpolation scheme that is implicitly included in the equation system. The refined subdomain can move with the average transport velocity of the subdomain. Very small time steps are required on a fine or a refined grid, because of the combined effect of the Courant and Peclet conditions. Therefore, we have developed a special upwind technique in small grid cells with high velocities (velocity suppression). We have neglected grid subdomains with very small concentration gradients (zero suppression). The resulting software, MODCALIF, is a three-dimensional, modularly constructed FORTRAN code. For convenience, the package names used by the well-known MODFLOW and MT3D computer programs are adopted, and the same input file structure and format is used, but the program presented here is separate and independent. Also, MODCALIF includes algorithms for variable density modeling and model calibration. The method is tested by comparison with an analytical solution, and illustrated by means of a two-dimensional theoretical example and three-dimensional simulations of the variable-density Cape Cod and SALTPOOL experiments. Crossing from fine to coarse grid produces numerical dispersion when the whole subdomain of interest is refined; however, we show that accurate solutions can be obtained using a fraction of the execution time required by uniformly fine-grid solutions. 相似文献
16.
Groundwater exchanges with lakes resulting from cyclical wet and dry climate extremes maintain lake levels in the environment in ways that are not well understood, in part because they remain difficult to simulate. To better understand the atypical groundwater interactions with lakes caused by climatic extremes, an original conceptual approach is introduced using MODFLOW‐2005 and a kinematic‐wave approximation to variably saturated flow that allows lake size and position in the basin to change while accurately representing the daily lake volume and three‐dimensional variably saturated groundwater flow responses in the basin. Daily groundwater interactions are simulated for a calibrated lake basin in Florida over a decade that included historic wet and dry departures from the average rainfall. The divergent climate extremes subjected nearly 70% of the maximum lakebed area and 75% of the maximum shoreline perimeter to both groundwater inflow and lake leakage. About half of the lakebed area subject to flow reversals also went dry. A flow‐through pattern present for 73% of the decade caused net leakage from the lake 80% of the time. Runoff from the saturated lake margin offset the groundwater deficit only about half of that time. A centripetal flow pattern present for 6% of the decade was important for maintaining the lake stage and generated 30% of all net groundwater inflow. Pumping effects superimposed on dry climate extremes induced the least frequent but most cautionary flow pattern with leakage from over 90% of the actual lakebed area. 相似文献
17.
byWillem Jan Zaadnoordijk 《Ground water》2009,47(5):723-726
The standard MODFLOW packages offer limited capabilities to model piecewise-linear boundary conditions to describe ground water–surface water interaction. Specifically, MODFLOW is incapable of representing a Cauchy-type boundary with different resistances for discharge or recharge conditions. Such a more sophisticated Cauchy boundary condition is needed to properly represent surface waters alternatively losing water through the bottom (high resistance) or gaining water mostly near the water surface (low resistance). One solution would be to create a new package for MODFLOW to accomplish this. However, it is also possible to combine multiple instances of standard packages in a single cell to the same effect. In this specific example, the general head boundary package is combined with the drain package to arrive at the desired piecewise-linear behavior. In doing so, the standard USGS MODFLOW version can be used without any modifications at the expense of a minor increase in preprocessing and postprocessing and computational effort. The extra preprocessing for creating the input and extra postprocessing to determine the water balance in terms of the physical entities from the MODFLOW cell fluxes per package can be taken care of by a user interface. 相似文献
18.
Christian D. Langevin Joseph D. Hughes Alden M. Provost Martijn J. Russcher Sorab Panday 《Ground water》2024,62(1):111-123
MODFLOW 6 is the latest in a line of six “core” versions of MODFLOW released by the U.S. Geological Survey. The MODFLOW 6 architecture supports incorporation of additional hydrologic processes, in addition to groundwater flow, and allows interaction between processes. The architecture supports multiple model instances and multiple types of models within a single simulation, a flexible approach to formulating and solving the equations that represent hydrologic processes, and recent advances in interoperability, which allow MODFLOW to be accessed and controlled by external programs. The present version of MODFLOW 6 consolidates popular capabilities available in MODFLOW variants, such as the unstructured grid support in MODFLOW-USG, the Newton-Raphson formulation in MODFLOW-NWT, and the support for partitioned stress boundaries in MODFLOW-CDSS. The flexible multi-model capability allows users to configure MODFLOW 6 simulations to represent the local-grid refinement (LGR) capabilities available in MODFLOW-LGR, the multi-species transport capabilities in MT3DMS, and the coupled variable-density capabilities available in SEAWAT. This paper provides a new, holistic and integrated overview of simulation capabilities made possible by the MODFLOW 6 architecture, and describes how ongoing and future development can take advantage of the program architecture to integrate new capabilities in a way that is minimally invasive and automatically compatible with the existing MODFLOW 6 code. 相似文献
19.
Charles R. Fitts Joshua Godwin Kathleen Feiner Charles McLane Seth Mullendore 《Ground water》2015,53(3):432-439
This paper presents the analytic element modeling approach implemented in the software AnAqSim for simulating steady groundwater flow with a sharp fresh‐salt interface in multilayer (three‐dimensional) aquifer systems. Compared with numerical methods for variable‐density interface modeling, this approach allows quick model construction and can yield useful guidance about the three‐dimensional configuration of an interface even at a large scale. The approach employs subdomains and multiple layers as outlined by Fitts (2010) with the addition of discharge potentials for shallow interface flow (Strack 1989). The following simplifying assumptions are made: steady flow, a sharp interface between fresh‐ and salt water, static salt water, and no resistance to vertical flow and hydrostatic heads within each fresh water layer. A key component of this approach is a transition to a thin fixed minimum fresh water thickness mode when the fresh water thickness approaches zero. This allows the solution to converge and determine the steady interface position without a long transient simulation. The approach is checked against the widely used numerical codes SEAWAT and SWI/MODFLOW and a hypothetical application of the method to a coastal wellfield is presented. 相似文献
20.
Theodore A. Endreny 《水文研究》2014,28(11):3824-3832
During the evolution of meander bends, the intra‐meander groundwater head gradients steepen and generate zones of accelerated water and nutrient intra‐meander fluxes important for ecosystem processes. This paper compares and contrasts three MODFLOW groundwater model packages based on their simulation of intra‐meander flux for two stages of meander evolution observed in a sandbox river table and one level of river bed clogging, where the hydraulic conductivity in the river bed is lower than in the adjacent aquifer. These packages are the Time‐Variant Specified Head package [constant head (CHD)], River package (RIV), and Streamflow‐Routing package (SFR2), each controlling the groundwater or river head bounding the intra‐meander region. The RIV and SFR2 packages fix river stage and allow for variation in groundwater head below the river, which is suggested for simulating intra‐meander flux for all sinuosities with and without river bed clogging whenever river bed parameters are available. The CHD package fixes below river groundwater head and fails to simulate intra‐meander head loss and flux in meanders with high sinuosity or river bed clogging. In low sinuosity meanders and in cases without river bed clogging, there were no significant differences between MODFLOW packages for simulating river intra‐meander head loss and flux. This research demonstrates why MODFLOW users need to consider the limitations of each package when simulating intra‐meander flux in reaches with river bed clogging, high sinuosity, or similarly steep hydraulic gradients. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献