共查询到20条相似文献，搜索用时 281 毫秒
1.
A transient flow modeling analysis for potential public-supply wells on western Cape Cod, Massachusetts, demonstrates the difference between transient and steady-state recharge areas can have important implications for wellhead protection. An example of a single pumping well illustrates that commonly used steady-state time-related capture areas do not represent the recharge area and travel times of water being pumped from the well until sufficient time has elapsed for steady-state flow conditions to be established. Until that time, transient recharge areas are needed to account for the portion of water discharging from the well that entered the aquifer before pumping started. An example of two pumping wells demonstrates the same area at the water table cannot supply water to more than one well under steady-state flow conditions. Transient recharge areas to multiple wells can overlap, however, until steady-state flow conditions are established. The same area can, therefore, be a source of water to more than one well during early pumping times, and the water pumped from a given well may derive from source areas, including contaminated areas, that do not lie within the well's steady-state recharge area. 相似文献
2.
This paper presents analytical solutions for steady-state, compressible two-phase flow through a wellbore under isothermal conditions using the drift flux conceptual model. Although only applicable to highly idealized systems, the analytical solutions are useful for verifying numerical simulation capabilities that can handle much more complicated systems, and can be used in their own right for gaining insight about two-phase flow processes in wells. The analytical solutions are obtained by solving the mixture momentum equation of steady-state, two-phase flow with an assumption that the two phases are immiscible. These analytical solutions describe the steady-state behavior of two-phase flow in the wellbore, including profiles of phase saturation, phase velocities, and pressure gradients, as affected by the total mass flow rate, phase mass fraction, and drift velocity (i.e., the slip between two phases). Close matching between the analytical solutions and numerical solutions for a hypothetical CO2 leakage problem as well as to field data from a CO2 production well indicates that the analytical solution is capable of capturing the major features of steady-state two-phase flow through an open wellbore, and that the related assumptions and simplifications are justified for many actual systems. In addition, we demonstrate the utility of the analytical solution to evaluate how the bottomhole pressure in a well in which CO2 is leaking upward responds to the mass flow rate of CO2-water mixture. 相似文献
3.
Timothy T. Eaton 《Ground water》2020,58(4):511-523
Water table fluctuation (WTF) methods are a primary and well-established way to determine groundwater recharge based on the direct response of the water table to precipitation input. An emerging complexity of recharge is whether it occurs as an episodic and transient process, or a continuous steady-state process, however, most studies have not focused on these short-term vs. long-term timescales, in part because of a lack of data resolution. Here, high-resolution (subhourly) precipitation and water-level data are analyzed for wells in the suburbs of New York City using two contrasting WTF approaches, with a common mathematical basis, that are suited to episodic and continuous processes. The resulting hourly recharge results, like the individual water-level records from comparable wells, are sensitive indicators of subtle differences in aquifer conditions such as thickness of the unsaturated zone, position in the flow system and localized preferential flow. While the episodic, transient approach excludes diffuse recharge by design, the continuous, steady-state approach is influenced by short-term precipitation events, and therefore integrates transient processes to some extent. However, the continuous, steady-state approach is subject to its own limitations relating to position in the aquifer system, and may overestimate recharge if aquifer conditions are not well understood. More widespread use of higher resolution data as well as understanding aquifer conditions and refining aquifer parameters would improve WTF recharge estimation. 相似文献
4.
The hyporheic zone is an ecologically important ecotone that describes the extent to which nutrient-rich surface waters penetrate the shallow subsurface adjacent to a flowing surface water body. Although steady-state models satisfactorily explain the incursion of surface water into the subsurface as a function of head gradients developed across streambed riffles, they fail to account for the depth that surface water is observed to penetrate the subsurface or for the extent to which the hyporheic zone develops adjacent to the stream channel. To investigate these issues, transient flow modeling has been conducted at the riffle scale and supported by data for an instrumented site in northern Ontario where stream-stage fluctuations are strictly regulated. Model results show that daily stream-stage fluctuations between 0.6 and 4 m produce oscillating solute flow paths that typically reduce residence times of water and solutes in the hyporheic zone from 60 days or more under steady-state conditions to less than 1 day. Furthermore, similar stream-stage fluctuations increase the depth that solutes pervade the subsurface and banks lateral to the stream from around 1 m under steady-state conditions to as much as 2 and 10 m, respectively. The results demonstrate that the transient flow conditions triggered in the subsurface by variable stream stage can exert a strong influence on hyporheic zone development and have important implications for the hyporheos. The results are especially important for hyporheic communities that may survive gradual changes to their living conditions by migrating to more hospitable aquatic habitats, but are unable to respond to rapid changes provoked by more extreme hydrological events. 相似文献
5.
6.
A new computer program, FLASH (Flow-Log Analysis of Single Holes), is presented for the analysis of borehole vertical flow logs. The code is based on an analytical solution for steady-state multilayer radial flow to a borehole. The code includes options for (1) discrete fractures and (2) multilayer aquifers. Given vertical flow profiles collected under both ambient and stressed (pumping or injection) conditions, the user can estimate fracture (or layer) transmissivities and far-field hydraulic heads. FLASH is coded in Microsoft Excel(5) with Visual Basic for Applications routines. The code supports manual and automated model calibration. 相似文献
7.
We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated. 相似文献
8.
《Earth and Planetary Science Letters》2006,241(1-2):202-210
Plumes rising from the core–mantle boundary (CMB) are often assumed to transport most, or all, of the heat conducted across the CMB. Here this assumption is explored using numerical convection models in idealized geometries that lead to a single plume under steady-state or near steady state conditions. Plume heat transport is calculated for different internal heating rates using two methods and compared to the CMB heat flux. For these conditions, it is found that the heat flux transported by plumes in the upper mantle is only a fraction of the core heat flux and, thus, core heat flow estimates derived from observed hotspots could be multiplied by a factor of several. 相似文献
9.
Sara M. Klemm Eric A. Seagren Jennifer G. Becker Alex S. Mayer 《Ground Water Monitoring & Remediation》2021,41(4):33-48
Porous aquifer materials are often characterized by layered heterogeneities that influence groundwater flow and present complexities in contaminant transport modeling. Such flow variations also have the potential to impact the dissolution flux from dense nonaqueous phase liquid (DNAPL) pools. This study examined how these heterogeneous flow conditions affected the dissolution of a tetrachloroethene (PCE) pool in a two-dimensional intermediate-scale flow cell containing coarse sand. A steady-state mass-balance approach was used to calculate the PCE dissolution rate at three different flow rates. As expected, aqueous PCE concentrations increased along the length of the PCE pool and higher flow rates decreased the aqueous PCE concentration in the effluent. Nonreactive tracer studies at two flow rates confirmed the presence of a vertical flow gradient, with the most rapid velocity located at the bottom of the tank. These results suggest that flow focusing occurred near the DNAPL pool. Effluent PCE concentrations and pool dissolution flux rates were compared to model predictions assuming local equilibrium (LE) conditions at the DNAPL pool/aqueous phase interface and a uniform distribution of flow. The LE model did not describe the data well, even over a wide range of PCE solubility and macroscopic transverse dispersivity values. Model predictions assuming nonequilibrium mass-transfer-limited conditions and accounting for vertical flow gradients, however, resulted in a better fit to the data. These results have important implications for evaluating DNAPL pool dissolution in the field where subsurface heterogeneities are likely to be present. 相似文献
10.
We present a methodology for identifying highly-localized flow channels embedded in a significantly less permeable medium using steady-state head and geometrical data. This situation is typical of fractured media where flows are often strongly channeled at the scales of interest (10 m–1 km). The objective is to identify both geometrical and hydraulic characteristics of the conducting structures. Channels are identified in decreasing order of importance by successive optimizations of an objective function. The identification strategy takes advantage of the hierarchical flow organization to restrict the dimension of the solution space of each individual optimization step. The characteristics of the secondary channels are strongly determined by the main flow channels. The latter are slightly modified by the secondary channels through the addition of a regularization term to the main channel characteristics in the objective function. As the objective function is strongly non-convex with numerous local minima, inversion is performed using a stochastic algorithm (simulated annealing). We assess the possibilities of the hierarchical identification strategy on simple synthetic steady-state flow configurations where hydraulic data are made up of 25 regularly spaced heads and of the boundary conditions. Those flow structures that are dominated by at most two simple channels can be identified with these head data only. Configurations comprising up to three complex and interconnected channels can still be identified with additional geometrical information including the distances of piezometers to their closest channel. The capabilities of the hierarchical identification strategy are limited to flow structures dominated by at most three equivalent flow channels. We finally discuss the perspectives of application of the method to transient-state data obtained on a more restricted number of piezometers. 相似文献
11.
Hydraulic Gradient Comparison Method to Estimate Aquifer Hydraulic Parameters Under Steady-State Conditions 总被引：1，自引：0，他引：1
The hydraulic gradient comparison method is an inverse method for estimation of aquifer hydraulic conductivity (or trans-missivity) and boundary conductance for a ground water flow model under steady-state conditions. This method, following formal optimization techniques, defines its objective function to minimize differences between interpreted (observed) and simulated hydraulic gradients, which results in minimization of differences between observed and simulated hydraulic heads. The key features of this method are that (1) the derived optimality conditions have an explicit form with a clear hydrology concept that is con-sistent with Darcy's law, and (2) the derived optimality conditions are spatially independent as they are a function of only local hydraulic conductivity and local hydraulic gradient. This second feature allows a multidimensional optimization problem to be solved by many one-dimensional optimization procedures simultaneously, which results in a substantial reduction in computation time. The results of the numerical performance testing on a heterogeneous hypothetical case confirm that minimizing gradient residuals in the entire model domain leads to minimizing head residuals. Application of the method in real-world projects requires rigorous conceptual model development, use of a global calibration target, and an iterative calibration proess. The conceptual model development includes interpretation of a potentiometric surface and estimation of other hydrologic parameters. This method has been applied to a wide range of real-world modeling projects, including the Rocky Mountain Arsenal and Rocky Flats sites in Colorado, which demonstrates that the method is efficient and practical. 相似文献
12.
by Tom Clemo 《Ground water》2010,48(1):68-78
A model coupling fluid hydraulics in a borehole with fluid flow in an aquifer is developed in this paper. Conservation of momentum is used to create a one-dimensional steady-state model of vertical flow in an open borehole combined with radially symmetric flow in an aquifer and with inflow to the well through the wellbore screen. Both laminar and turbulent wellbore conditions are treated. The influence of inflow through the wellbore screen on vertical flow in the wellbore is included, using a relation developed by Siwoń (1987) . The influence of inflow reduces the predicted vertical variation in head up to 15% compared to a calculation of head losses due to fluid acceleration and the conventional Colebrook-White formulation of friction losses in a circular pipe. The wellbore flow model is embedded into the MODFLOW-2000 ground water flow code. The nonlinear conservation of momentum equations are iteratively linearized to calculate the conductance terms for vertical flow in the wellbore. The resulting simulations agree favorably with previously published results when the model is adjusted to meet the assumptions of the previous coupled models. 相似文献
13.
Dongwei Sun Ning Luo Aaron Vandenhoff Wesley McCall Zhanfeng Zhao Chenxi Wang David L. Rudolph Walter A. Illman 《Ground water》2024,62(3):384-404
Significant efforts have been expended for improved characterization of hydraulic conductivity (K) and specific storage (Ss) to better understand groundwater flow and contaminant transport processes. Conventional methods including grain size analyses (GSA), permeameter, slug, and pumping tests have been utilized extensively, while Direct Push-based Hydraulic Profiling Tool (HPT) surveys have been developed to obtain high-resolution K estimates. Moreover, inverse modeling approaches based on geology-based zonations, and highly parameterized Hydraulic Tomography (HT) have also been advanced to map spatial variations of K and Ss between and beyond boreholes. While different methods are available, it is unclear which one yields K estimates that are most useful for high resolution predictions of groundwater flow. Therefore, the main objective of this study is to evaluate various K estimates at a highly heterogeneous field site obtained with three categories of characterization techniques including: (1) conventional methods (GSA, permeameter, and slug tests); (2) HPT surveys; and (3) inverse modeling based on geology-based zonations and highly parameterized approaches. The performance of each approach is first qualitatively analyzed by comparing K estimates to site geology. Then, steady-state and transient groundwater flow models are employed to quantitatively assess various K estimates by simulating pumping tests not used for parameter estimation. Results reveal that inverse modeling approaches yield the best drawdown predictions under both steady and transient conditions. In contrast, conventional methods and HPT surveys yield biased predictions. Based on our research, it appears that inverse modeling and data fusion are necessary steps in predicting accurate groundwater flow behavior. 相似文献
14.
This article outlines analytical solutions to quantify the length scale associated with “upstream dispersion,” the artificial movement of solutes in the opposite direction to groundwater flow, in solute transport models. Upstream dispersion is an unwanted artifact in common applications of the advection-dispersion equation (ADE) in problems involving groundwater flow in the direction of increasing solute concentrations. Simple formulae for estimating the one-dimensional distance of upstream dispersion are provided. These show that under idealized conditions (i.e., steady-state flow and transport, and a homogeneous aquifer), upstream dispersion may be a function of only longitudinal dispersivity. The scale of upstream dispersion in a selection of previously presented situations is approximated to highlight the utility of the presented formulae and the relevance of this ADE anomaly in common transport problems. Additionally, the analytical solution is applied in a hypothetical scenario to guide the modification of dispersion parameters to minimize upstream dispersion. 相似文献
15.
J. R. L. Allen 《地球表面变化过程与地形》1976,1(4):361-374
In natural environments, where unsteady flows prevail, the common one-way bed forms (ripples, dunes, antidunes, parting lineations on plane beds) show delayed responses to changes of flow. A dual classification may then be necessary for the features: 1. Geometric, with no implication of flow conditions. 2. Geometric-hydraulic, in which the shape of the observed forms is combined with the character of the concurrent flow. The delayed responses arise because, under the simplest steady-state equilibrium conditions, the bed forms behave deterministically as populations, which take average dimensions controlled by the bed-material and/or flow conditions. The response of the bed to a change of flow is a rearrangement of the particles forming the bed in an attempt to meet the new conditions. The rate of rearrangement, however, is controlled by the sediment transport rate, in turn determined by the flow conditions, and so is finite. Two specific mechanisms of change of opposing tendency are quantifiable: 1. Creation-destruction of forms, where the new forms are better adjusted than predecessors. 2. The imperfect modification of existing forms during their life-spans. Both rates may be controlled by a combination of flow variables, the attributes of the bed forms, and ‘constants’ specific to each kind of form. A better knowledge of the unsteady responses of these bed forms should increase our understanding of and control over river and tidal systems and may provide models for a better appreciation of still larger features in landscapes and waterscapes. 相似文献
16.
A three-dimensional prognostic hydrodynamic model in cross sectional form is used to examine the influence of bottom friction,
mixing and topography upon the spin-down and steady-state circulation in a cold water bottom-dome. Parameters characteristic
of the Irish Sea or Yellow Sea cold water domes are used. In all calculations, motion is induced by specifying an initial
temperature distribution characteristic of the dome, and an associated along frontal flow. The spin-down of the dome is found
to be influenced by the coefficient of bottom friction, with a typical time scale of order 10 days, and in general to be independent
of the chosen initial vertical profile of along frontal flow. However, in the case in which the along frontal flow is such
that the near bed velocity is zero, then bottom stress is also zero, and there is no appreciable spin-down. Calculations showed
that the formulation of viscosity and diffusivity had a greater effect upon the steady-state circulation than topography,
suggesting that background mixing of tidal origin is important. The lack of topographic influence was due mainly to the formulation
of the initial conditions which were taken to be independent of topography. The steady-state circulation was characterized
by a cyclonic flow in the surface region, with an anti-cyclonic current near the bed, where frictional effects produced a
bottom Ekman layer and an across frontal flow. This gave rise to vertical circulation cells in the frontal region of the dome
with prevailing downwelling motion inside the dome. A detailed analysis of the dynamic balance of the various terms in the
hydrodynamic equations yielded insight into the processes controlling the steady-state circulation in cold water domes.
Responsible Editor: Phil Dyke 相似文献
17.
This paper describes the use of single-well tracer dilution techniques to resolve the rate of light nonaqueous phase liquid (LNAPL) flow through wells and the adjacent geologic formation. Laboratory studies are presented in which a fluorescing tracer is added to LNAPL in wells. An in-well mixer keeps the tracer well mixed in the LNAPL. Tracer concentrations in LNAPL are measured through time using a fiber optic cable and a spectrometer. Results indicate that the rate of tracer depletion is proportional to the rate of LNAPL flow through the well and the adjacent formation. Tracer dilution methods are demonstrated for vertically averaged LNAPL Darcy velocities of 0.00048 to 0.11 m/d and LNAPL thicknesses of 9 to 24 cm. Over the range of conditions studied, results agree closely with steady-state LNAPL flow rates imposed by pumping. A key parameter for estimating LNAPL flow rates in the formation is the flow convergence factor alpha. Measured convergence factors for 0.030-inch wire wrap, 0.030-inch-slotted polyvinyl chloride (PVC), and 0.010-inch-slotted PVC are 1.7, 0.91, and 0.79, respectively. In addition, methods for using tracer dilution data to determine formation transmissivity to LNAPL are presented. Results suggest that single-well tracer dilution techniques are a viable approach for measuring in situ LNAPL flow and formation transmissivity to LNAPL. 相似文献
18.
This investigation was undertaken to develop an integrated method of downhole fracture characterization using a tracer. The method presented can be used to locate water-bearing fractures that intersect the well, to determine the ambient fracture flow rate and hydraulic head, and to calculate fracture transmissivity. The method was tested in two fractured crystalline bedrock wells located at the University of Connecticut in Storrs. The method entails injecting a tracer (uranine dye) into the well, while at the same time water is pumped out of the well. After steady-state conditions are reached, a borehole tracer concentration profile is developed. The dilution of the tracer is used to locate the inflowing fractures and to determine their flow rate. The fracture flow rate, plus the drawdown in the well, is then used to determine the fracture hydraulic head, transmissivity, and ambient flow rate. 相似文献
19.
Existing analytical solutions to 2D and 3D contaminant transport problems are limited by the mathematically convenient assumption of uniform flow. An approximate method is developed herein for coordinate mapping of 2D (vertically-averaged) transport solutions to non-uniform steady-state irrotational and divergence-free flow fields in single-layer aquifers. The method enables existing analytical transport solutions to be applied to aquifer systems with wells, non-uniform saturated thickness, surface water features, and (to a limited degree) heterogeneous hydraulic conductivity and recharge. This mass-conservative coordinate mapping approach is inexact in its approximation of the dispersion process but is still sufficiently accurate for many simple flow systems. The degree of model error is directly proportional to the variation of velocity magnitude within the domain. These mapped analytical solutions are compared to numerical simulation results and the coordinate mapping errors are investigated. The methods described herein may be used in the traditional capacity of analytical transport models, i.e., screening and preliminary site assessment, without sacrificing accuracy by assuming locally uniform flow conditions or applying an ad-hoc coordinate transformation. The solutions benefit from the traditional advantages of analytical methods, particularly the removal of artifacts due to spatial and temporal discretization: no time-stepping or numerical discretization is required. 相似文献
20.
A solution is given for the travel time to a well pumping an ideal, horizontal unconfined aquifer, under steady-state conditions, when recharge from infiltration is negligible. Three forms of the solution are provided: a closed-form solution, an integral to be calculated in a worksheet, and a simple equation. The three forms of the solution give travel times nearly identical to those obtained using a finite-element code for saturated and unsaturated flow and particle tracking. 相似文献