首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
Physical properties of alluvial environments typically feature a high degree of anisotropy and are characterized by dynamic interactions between the surface and the subsurface. Hydrogeological models are often calibrated under the assumptions of isotropic hydraulic conductivity fields and steady-state conditions. We aim at understanding how these simplifications affect predictions of the water table using physically based models and advanced calibration and uncertainty analysis approaches based on singular value decomposition and Bayesian analysis. Specifically, we present an analysis of the information content provided by steady-state hydraulic data compared to transient data with respect to the estimation of aquifer and riverbed hydraulic properties. We show that assuming isotropy or fixed anisotropy may generate biases both in the estimation of aquifer and riverbed parameters as well as in the predictive uncertainty of the water table. We further demonstrate that the information content provided by steady-state hydraulic heads is insufficient to jointly estimate the aquifer anisotropy together with the aquifer and riverbed hydraulic conductivities and that transient data can help to reduce the predictive uncertainty to a greater extent. The outcomes of the synthetic analysis are applied to the calibration of a dynamic and anisotropic alluvial aquifer in Switzerland (The Rhône River). The results of the synthetic and real world modeling and calibration exercises documented herein provide insight on future data acquisition as well as modeling and calibration strategies for these environments. They also provide an incentive for evaluation and estimation of commonly made simplifying assumptions in order to prevent underestimation of the predictive uncertainty.  相似文献   

2.
Numerical experiments with steady-state ground water flow models show that spiraling flow lines occur in layered aquifers that have different anisotropic horizontal hydraulic conductivities in adjacent layers. Bundles of such flow lines turning in the same direction can be referred to as ground water whirls. An anisotropic layered block in a field of uniform horizontal flow results in one or more whirls with their axes in the uniform flow direction. The number of whirls depends on the number of interfaces between layers with different anisotropic properties. For flow to a well in an aquifer consisting of two anisotropic layers, with perpendicular major principal directions, whirls are found to occur in quadrants that are bounded by the principal directions of the hydraulic conductivity. The combined effect of flow to a well and a layered anisotropy implies that a single well in a system with a single anisotropic layer within an otherwise isotropic aquifer causes eight whirls. All adjacent whirls rotate in opposite directions.  相似文献   

3.
A steady/quasi-steady model is developed for predicting flow into a partially penetrating well with skin zone in a confined aquifer overlying an impervious layer. The model takes into account flow through the bottom of the wellbore, finite skin thickness and finite horizontal and vertical extent of the aquifer. Moreover, the solution can be easily extended to include the mixed-type boundary condition at the well face, where a Dirichlet in the form of a specified hydraulic head and a Neumann in the form of zero flux coexist at the same time at different portions of the well face. The validity of the proposed solution is tested by comparing a few results obtained from the developed model with corresponding results obtained by analytical and numerical means. The study shows that, among other factors remaining constant, both the horizontal and vertical extent of an artesian aquifer, thickness of the skin zone, bottom flow and conductivity contrast of the skin and formation zones, play an important part in deciding flow to a well dug in the aquifer, and hence these factors must be considered while analyzing the problem. The model proposed here can be used to estimate skin thickness as well as hydraulic conductivities of the skin and formation zones of a well with skin zone in an artesian aquifer underlain by an impervious layer by utilizing pumping test data falling in the steady or quasi-steady state of a typical pumping test. As the proposed solution is of a general nature in the sense that it can handle, apart from partial penetration and bottom flow, the finite size skin zone and finite horizontal and vertical extent of an artesian aquifer together with the mixed-type boundary condition at the well face, it is hoped that the predictions coming out of the model will be more realistic than those obtained using solutions developed with more stringent assumptions.  相似文献   

4.
The transient flowmeter test (TFMT) provides more information about the well–aquifer system than the traditional quasi-steady-state flowmeter test (QFMT). The TFMT duration may be much shorter than that of a QFMT, which is desirable at highly contaminated sites where the extracted water has to be treated as hazardous waste. Here we present the TFMT model that accounts for inter-layer crossflow, a thick skin surrounding the well, and wellbore storage. The model is derived under the simplifying assumptions of the pseudo-steady-state inter-layer crossflow and the uniform wellface flux within each layer. The semi-analytic solution is inverted numerically from the Laplace domain to the time domain. Layer and skin parameters are estimated from the TFMT data via the modified Levenberg–Marquardt algorithm. The estimation is robust when the initial parameter guesses are close to their true values. Otherwise, a computationally expensive search among the local minima of the objective function is necessary to find the parameter estimates. The modeling errors and the associated parameter estimation errors are evaluated in a number of synthetic TFMTs and compared to the corresponding results obtained with a general numerical model that relaxes the two simplifying assumptions. The TFMT provides reasonably accurate estimates of hydraulic conductivities for the aquifer layers and the damaged skins and order-of-magnitude estimates of layer specific storativities and hydraulic conductivities for the normal skin. The skin specific storativities should not be estimated from a TFMT. Multi-rate TFMTs with a step-variable pumping rate yield significantly more accurate parameters than constant-pumping-rate TFMTs. The calculated modeling errors may be useful in estimating the magnitude of parameter estimation errors from the TFMT. Our field tests in a coastal aquifer at the Lizzie Site in North Carolina (USA) demonstrate the feasibility of a TFMT for aquifer characterization. The downhole hydraulic conductivity profiles from our field and synthetic TFMTs are consistent with the corresponding profiles from QFMTs.  相似文献   

5.
Wells in aquifers of loose collapsible sediment are cased so that they have a blind wall and gain water only from the bottom. The hydraulic gradient established at the bottom of these wells during pumping brings the aquifer materials in a quicksand state, which may cause abrasion of pipes and pumps and even the destruction of well structure. To examine the quicksand occurrence, an analytical solution for the steady flow to a partially penetrating blind‐wall well in a confined aquifer is developed. The validity of the proposed solution is evaluated numerically. The sensitivity of maximum vertical gradient along the well bottom in response to aquifer and well parameters is examined. The solution is presented in the form of dimensionless‐type curves and equations that can be easily used to design the safe pumping rate and optimum well geometry to protect the well against sand production. The solution incorporates the anisotropy of aquifer materials and can also be used to determine the hydraulic conductivity of the aquifer. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

6.
This study investigates the behavior of flux and head in a strongly heterogeneous three-dimensional aquifer system. The analyses relied on data from 520 slug tests together with 38,000 one-foot core intervals lithological data from the site of the General Separations Area in central Savannah River Site, South Carolina, USA. The skewness in the hydraulic conductivity histograms supported the geologic information for the top two aquifers, but revealed stronger clay content, than was reported for the bottom aquifer. The log-normal distribution model described adequately the hydraulic conductivity measurements for all three aquifers although, other distributions described equally well the bottom aquifer measurements. No apparent anisotropy on the horizontal plane was found for the three aquifers, but ratios of horizontal to vertical correlation lengths between 33 and 75 indicated a strong stratification at the site. Three-dimensional Monte Carlo stochastic simulations utilized a grid with larger elements than the support volume of measurements, but of sub-REV (representative elementary volume) dimensions. This necessitated, on one hand, the use of upscaled hydraulic conductivity expressions, but on the other hand did not allow for the use of anisotropic effective hydraulic conductivity expressions (Sarris and Paleologos in J Stoch Environ Res Risk Assess 18: 188–197, 2004). Flux mean and standard deviations components were evaluated on three vertical cross-sections. The mean and variance of the horizontal flux component normal to a no-flow boundary tended to zero at approximately two to three integral scales from that boundary. Close to a prescribed head boundary both the mean and variance of the horizontal flux component normal to the boundary increased from a stable value attained at a distance of about five integral scales from that boundary. The velocity field 〈qx〉 was found to be mildly anisotropic in the top two aquifers, becoming highly anisotropic in the bottom aquifer; 〈qy〉 was anisotropic in all three aquifers with directions of high continuity normal to those of the 〈qx〉 field; finally, 〈qz〉 was highly anisotropic in all three aquifers, with higher continuity along the east–west direction. The mean head field was found to be continuous, despite the high heterogeneity of the underlying hydraulic conductivity field. Directions of high continuity were in alignment with field boundaries and mean flow direction. Conditioning did not influence significantly the expected value of the flux terms, with more pronounced being the effect on the standard deviation of the flux vector components. Conditioning reduced the standard deviations of the horizontal flux components by as much as 50% in the bottom aquifer. Variability in the head cross-sections was affected only marginally, with an average 10% reduction in the respective standard deviation. Finally, the location of the conditioning data did not appear to have a significant effect on the surrounding area, with uniform reduction in standard deviations.  相似文献   

7.
The effects of wellbore‐wall compression from rough excavation on monitored groundwater levels and qualities under natural hydraulic gradient conditions were investigated in a shallow clayey Andisol aquifer. Nine wellbores reaching the underlying aquitard at about 2.6‐m depth were constructed by dynamic cone penetrometry to mimic rough wellbore construction. Five of these were constructed under wet aquifer soil conditions and the remaining four under dry conditions. A 15‐month period monitoring showed that the groundwater levels in the wellbores constructed under wet conditions responded significantly in retard of, and in narrower ranges than, those constructed under dry conditions. The wellbore‐wall hydraulic conductivities at the former wellbores were calculated to be more than one to two orders of magnitude lower than those at the latter ones. Furthermore, remarkable nitrate removal attributable to the occurrence of a heterotrophic denitrification was observed in one of the former wellbores. In contrast, the groundwater levels and qualities in the latter wellbores appeared to be generally similar to those monitored in the conventional soil coring and augering‐derived wellbores. Our results suggest that the wellbore‐wall compression induced by rough excavation under wet and soft aquifer soil conditions leads to a substantial decrease in the wellbore‐wall hydraulic conductivity, which in turn can lead to unreliable groundwater levels and qualities. This problem can occur in clayey Andisols whenever the aquifer soil is wet; however, the problem can be largely avoided by constructing the wellbore under dry and hard aquifer soil conditions.  相似文献   

8.
Two analytical solution methods are presented for regional steady-state groundwater flow in a two-dimensional stratified aquifer cross section where the water table is approximated by the topographic surface. For the first solution, the surficial aquifer is represented as a set of dipping parallel layers with different, but piecewise constant, anisotropic hydraulic conductivities, where the anisotropy is aligned with the dip of the layered formation. The model may be viewed as a generalization of the solutions developed by [Tóth JA. A theoretical analysis of groundwater flows in small drainage basins. J Geophys Res 1963;68(16):4795–812; Freeze R, Witherspoon P. Theoretical analysis of regional groundwater flow 1) analytical and numerical solution to the mathematical model, water resources research. Water Resour Res 1966;2(4):641–56; Selim HM. Water flow through multilayered stratified hillside. Water Resour Res 1975;11:949–57] to an multi-layer aquifer with general anisotropy, layer orientation, and a topographic surface that may intersect multiple layers. The second solution presumes curved (syncline) layer stratification with layer-dependent anisotropy aligned with the polar coordinate system. Both solutions are exact everywhere in the domain except at the topographic surface, where a Dirichlet condition is met in a least-squared sense at a set of control points; the governing equation and no-flow/continuity conditions are met exactly. The solutions are derived and demonstrated on multiple test cases. The error incurred at the location where the layer boundaries intersect the surface is assessed.  相似文献   

9.
Fitts CR 《Ground water》2006,44(1):99-101
Although most current applications of the analytic element method are formulated for isotropic hydraulic conductivity, anisotropic domains can be modeled with analytic elements using the well-known coordinate transformation where one coordinate axis is scaled by the square root of the anisotropy ratio. If the standard analytic solution for steady radial flow to a well is used with this coordinate transformation, the resulting solution correctly models the far field but it does not meet the constant head boundary condition at the well radius. This could be a significant shortcoming if you are interested in the flow field close to the well or want to estimate the head at the pumping well. A new solution for two-dimensional steady flow to a well in an anisotropic domain is presented. This solution satisfies the governing equations exactly and meets the constant head boundary condition at the well radius exactly. It was derived using a conformal mapping.  相似文献   

10.
Sepúlveda N 《Ground water》2008,46(1):144-155
An analytical solution for three-dimensional (3D) flow in the storative semiconfining layers of a leaky aquifer fully penetrated by a production well is developed in this article to provide a method from which accurate hydraulic parameters in the semiconfining layers can be derived from aquifer test data. The analysis of synthetic aquifer test data with the 3D analytical solution in the semiconfining layers provided more accurate optimal hydraulic parameters than those derived using the available quasi-two-dimensional (2D) solution. Differences between the 3D and 2D flow solutions in the semiconfining layers become larger when a no flow boundary condition is imposed at either at the top of the upper semiconfining layer or at the bottom of the lower semiconfining layer or when the hydraulic conductivity ratio of the semiconfining layer to the aquifer is larger than 0.001. In addition, differences between the 3D and 2D flow solutions in the semiconfining layers are illustrated when the thickness ratio of the semiconfining layer to the aquifer is changed. Analysis of water level data from two hypothetical and one real aquifer test showed that the 3D solution in the semiconfining layers provides lower correlation coefficients among hydraulic parameters than the 2D solution.  相似文献   

11.
Inflow to a tunnel is a great public concern and is closely related to groundwater hydrology, geotechnical engineering, and mining engineering, among other disciplines. Rapid computation of inflow to a tunnel provides a timely means for quickly assessing the inflow discharge, thus is critical for safe operation of tunnels. Dewatering of tunnels is another engineering practice that should be planned. In this study, an analytical solution of the inflow to a tunnel in a fractured unconfined aquifer is obtained. The solution takes into account either the spherical or slab-shaped matrix block and the unsteady state interporosity flow. The instantaneous drainage water table and anisotropic hydraulic conductivities of the fractures network are also considered. Both uniform flux and uniform head boundary condition are considered to simulate the constant head boundary condition in the tunnel. The effects of the hydraulic parameters of the fractured aquifer on the inflow variation of the tunnel are explored. The application of the presented solution to obtain the optimum location and discharge of the well to minimize the inflow to a tunnel is illustrated.  相似文献   

12.
An analytical model is presented for the analysis of constant flux tests conducted in a phreatic aquifer having a partially penetrating well with a finite thickness skin. The solution is derived in the Laplace transform domain for the drawdown in the pumping well, skin and formation regions. The time-domain solution in terms of the aquifer drawdown is then obtained from the numerical inversion of the Laplace transform and presented as dimensionless drawdown–time curves. The derived solution is used to investigate the effects of the hydraulic conductivity contrast between the skin and formation, in addition to wellbore storage, skin thickness, delayed yield, partial penetration and distance to the observation well. The results of the developed solution were compared with those from an existing solution for the case of an infinitesimally thin skin. The latter solution can never approximate that for the developed finite skin. Dimensionless drawdown–time curves were compared with the other published results for a confined aquifer. Positive skin effects are reflected in the early time and disappear in the intermediate and late time aquifer responses. But in the case of negative skin this is reversed and the negative skin also tends to disguise the wellbore storage effect. A thick negative skin lowers the overall drawdown in the aquifer and leads to more persistent delayed drainage. Partial penetration increases the drawdown in the case of a positive skin; however its effect is masked by the negative skin. The influence of a negative skin is pronounced over a broad range of radial distances. At distant observation points the influence of a positive skin is too small to be reflected in early and intermediate time pumping test data and consequently the type curve takes its asymptotic form.  相似文献   

13.
An aquifer containing a skin zone is considered as a two-zone system. A mathematical model describing the head distribution is presented for a slug test performed in a two-zone confined aquifer system. A closed-form solution for the model is derived by Laplace transforms and Bromwich integral. This new solution is used to investigate the effects of skin type, skin thickness, and the contrast of skin transmissivity to formation transmissivity on the distributions of dimensionless hydraulic head. The results indicate that the effect of skin type is marked if the slug-test data is obtained from a radial two-zone aquifer system. The dimensionless well water level increases with the dimensionless positive skin thickness and decreases as the dimensionless negative skin thickness increases. In addition, the distribution of dimensionless well water level due to the slug test depends on the hydraulic properties of both the wellbore skin and formation zones.  相似文献   

14.
A physically based inverse method is developed using hybrid formulation and coordinate transform to simultaneously estimate hydraulic conductivity tensors, steady‐state flow field, and boundary conditions for a confined aquifer under ambient flow or pumping condition. Unlike existing indirect inversion techniques, the physically based method does not require forward simulations to assess model‐data misfits. It imposes continuity of hydraulic head and Darcy fluxes in the model domain while incorporating observations (hydraulic heads, Darcy fluxes, or well rates) at measurement locations. Given sufficient measurements, it yields a well‐posed inverse system of equations that can be solved efficiently with coarse grids and nonlinear optimization. When pumping and injection are active, well rates are used as measurements and flux sampling is not needed. The method is successfully tested on synthetic aquifer problems with regular and irregular geometries, different hydrofacies and flow patterns, and increasing conductivity anisotropy ratios. All problems yield stable inverse solutions under increasing head measurement errors. For a given set of observations, inversion accuracy is strongly affected by the conductivity anisotropy ratio. Conductivity estimation is also affected by flow pattern: within a hydrofacies, when Darcy flux component is very small, the corresponding directional conductivity perpendicular to streamlines becomes less identifiable. Finally, inversion is successful even if the location of aquifer boundaries is unknown. In this case, the inversion domain is defined by the location of the measurements.  相似文献   

15.
Wang X  Xu J  Cai Z 《Ground water》2012,50(3):348-354
A new finite layer method (FLM) is presented in this paper for transient flow analysis in layered radial two-zone aquifer systems. A radial two-zone system is an aquifer configuration in which a circular aquifer with finite radius is surrounded by a matrix possessing different permeability and storage properties. The aquifers can be pumped from fully or partially penetrating wells of infinitesimal radius. The trial function for drawdown is obtained through the use of piecewise linear correction functions in the present method. The trial function can satisfy the continuity conditions of flow and possess an appropriate continuity of C(0) at the two-zone interface. On the basis of Galerkin's method and the continuity condition of flow, the finite layer formulation is derived. The proposed method can cope with the anisotropy and layered heterogeneity in radial two-zone aquifer systems. Several numerical examples are presented to verify the validity of the present method through comparison with the analytical solution and the numerical results based on the finite difference method, in which a test of three-dimensional (3D) flow to a partially penetrating well in anisotropic two-zone aquifers is included. Furthermore, an additional application in simulating the two-zone flow in aquitard-aquifer systems is presented to demonstrate the applicability of FLM in modeling flow in more complex aquifer systems.  相似文献   

16.
We have developed a new method to analyze the power law based non-Darcian flow toward a well in a confined aquifer with and without wellbore storage. This method is based on a combination of the linearization approximation of the non-Darcian flow equation and the Laplace transform. Analytical solutions of steady-state and late time drawdowns are obtained. Semi-analytical solutions of the drawdowns at any distance and time are computed by using the Stehfest numerical inverse Laplace transform. The results of this study agree perfectly with previous Theis solution for an infinitesimal well and with the Papadopulos and Cooper’s solution for a finite-diameter well under the special case of Darcian flow. The Boltzmann transform, which is commonly employed for solving non-Darcian flow problems before, is problematic for studying radial non-Darcian flow. Comparison of drawdowns obtained by our proposed method and the Boltzmann transform method suggests that the Boltzmann transform method differs from the linearization method at early and moderate times, and it yields similar results as the linearization method at late times. If the power index n and the quasi hydraulic conductivity k get larger, drawdowns at late times will become less, regardless of the wellbore storage. When n is larger, flow approaches steady state earlier. The drawdown at steady state is approximately proportional to r1−n, where r is the radial distance from the pumping well. The late time drawdown is a superposition of the steady-state solution and a negative time-dependent term that is proportional to t(1−n)/(3−n), where t is the time.  相似文献   

17.
A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies.  相似文献   

18.
《Advances in water resources》2007,30(4):1016-1026
We have proved that the Hantush’s model [Hantush MS. Wells near streams with semipervious beds. J Geophys Res 1965;70:2829–38] in a half-domain can be extended to a whole-domain and becomes identical to that of Hunt [Hunt B. Unsteady stream depletion from ground water pumping. Ground Water 1999;37(1):98–102] for a shallow and infinitely narrow stream, provided that the Dupuit assumption is adopted. This proof helps correct a false concept that regards the Hantush’s model as less useful because of its fully penetrating stream assumption. This study deals with interaction of an aquifer with two parallel streams based on the Hantush’s model. Semi-analytical solutions are obtained based on rigorous mass conservation requirement by maintaining continuity of flux and head along the aquifer–streambed boundaries. This study shows that the hydraulic conductivity ratio of the two streambeds appears to be the most important factor controlling the stream–aquifer interaction, followed by a less important role played by the thickness ratio of the two streambeds. When the low-permeability streambeds do not exist, the steady-state stream depletion from one stream is linearly proportional to the ratio of the shortest distance from the pumping well to the other stream over the shortest distance between the two streams. When the low-permeability streambeds are presented, similar conclusion can be drawn except that the stream depletion now also strongly depends on the hydraulic conductivity ratio of the two streambeds. When the values of the hydraulic conductivity of the two streambeds are different by an order of magnitude, the location of the pumping well that receives equal flux from two streams can be off the middle-line between the two streams by nearly 90%.  相似文献   

19.
A dome-shaped layer can be selected as a storage site for fluid injection. In this study, we develop a mathematical model for simulating transient head distribution in a heterogeneous and anisotropic dome-shaped layer due to a constant-head injection in a fully penetrating well. In the model, a form of step change is adopted to approximate the upper and lower boundaries of the dome and then the layer is split into two regions. The Laplace-domain solution of the model is developed using the Laplace transform and method of separation of variables. The transient injection rate at wellbore can then be obtained based on Darcy’s law and Bromwich integral method. The predicted head contours from the head solution show significant vertical flow components near the location of step change in the dome reservoir. The results of sensitivity analysis indicate that the hydraulic conductivity is the most sensitive parameter and the specific storage is the least sensitive one to the injection rate after a short period of injection time. In addition, the injection rate for a dome reservoir is also very sensitive to the change of the height for the reservoir near the injection well (first region) at a very early injection time. In contrast, the injection rate is more sensitive to the change of the height of the second region than that of the first region at late time. This analytical solution may be used as a primary tool to assess the capacity of fluid injection to various dome reservoirs.  相似文献   

20.
This study develops a mathematical model for describing the steady-state head response to fluid injection into a fully penetrating well in a heterogeneous and anisotropic anticline reservoir. In the model, the upper boundary of the anticline reservoir is approximated by a form of step change in reservoir thickness and the domain of the reservoir is divided into two regions with different hydraulic conductivities. By virtue of the properties of Fourier series, the method of separation of variables is employed to develop the analytical solution of the model.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号