共查询到20条相似文献，搜索用时 31 毫秒
1.
Closed‐form solutions are proposed for natural seepage in semiconfined (leaky) aquifers such as those existing below the massive Champlain Sea clay layers in the Saint‐Lawrence River Valley. The solutions are for an ideal horizontal leaky aquifer below an ideal aquitard that may have either a constant thickness and a constant hydraulic head at its surface, or a variable thickness and a variable hydraulic head at its surface. A few simplifying assumptions were needed to obtain the closed‐form solutions. These have been verified using a finite element method, which did not make any of the assumptions but gave an excellent agreement for hydraulic heads and groundwater velocities. For example, the difference between the two solutions was smaller than 1 mm for variations in the 5 to 8 m range for the hydraulic head in the semiconfined aquifer. Note that fitting the hydraulic head data of monitoring wells to the theoretical solutions gives only the ratio of the aquifer and aquitard hydraulic conductivities, a clear case of multiple solutions for an inverse problem. Consequently, field permeability tests in the aquitard and the aquifer, and pumping tests in the aquifer, are still needed to determine the hydraulic conductivity values. 相似文献
2.
Abstract This paper develops a new analytical solution for the aquifer system, which comprises an unconfined aquifer on the top, a semi-confined aquifer at the bottom and an aquitard between them. This new solution is derived from the Boussinesq equation for the unconfined aquifer and one-dimensional leaky confined flow equation for the lower aquifer using the perturbation method, considering the water table over-height at the remote boundary. The head fluctuation predicted from this solution is generally greater than the one solved from the linearized Boussinesq equation when the ratio of the tidal amplitude to the thickness of unconfined aquifer is large. It is found that both submarine groundwater discharges from upper and lower aquifers increase with tidal amplitude–aquifer thickness ratio and may be underestimated if the discharge is calculated based on the average head fluctuation. The effects of the aquifer parameters and linearization of the Boussinesq equation on the normalized head fluctuation are also investigated. Editor D. Koutsoyiannis; Associate editor J. Simunek Citation Chuang, M.-H., Mahdi, A.-A. and Yeh, H.-D., 2012. A perturbation solution for head fluctuations in a coastal leaky aquifer system considering water table over-height. Hydrological Sciences Journal, 57 (1), 162–172. 相似文献
3.
In situ Mixing of Organic Matter Decreases Hydraulic Conductivity of Denitrification Walls in Sand Aquifers 总被引：2，自引：0，他引：2
Gregory F. Barkle Louis A. Schipper Craig P. Burgess Brett D.M. Painter 《Ground Water Monitoring & Remediation》2008,28(1):57-64
In a previous study, a denitrification wall was constructed in a sand aquifer using sawdust as the carbon substrate. Ground water bypassed around this sawdust wall due to reduced hydraulic conductivity. We investigated potential reasons for this by testing two new walls and conducting laboratory studies. The first wall was constructed by mixing aquifer material in situ without substrate addition to investigate the effects of the construction technique (mixed wall). A second, biochip wall, was constructed using coarse wood chips to determine the effect of size of the particles in the amendment on hydraulic conductivity. The aquifer hydraulic conductivity was 35.4 m/d, while in the mixed wall it was 2.8 m/d and in the biochip wall 3.4 m/d. This indicated that the mixing of the aquifer sands below the water table allowed the particles to re-sort themselves into a matrix with a significantly lower hydraulic conductivity than the process that originally formed the aquifer. The addition of a coarser substrate in the biochip wall significantly increased total porosity and decreased bulk density, but hydraulic conductivity remained low compared to the aquifer. Laboratory cores of aquifer sand mixed under dry and wet conditions mimicked the reduction in hydraulic conductivity observed in the field within the mixed wall. The addition of sawdust to the laboratory cores resulted in a significantly higher hydraulic conductivity when mixed dry compared to cores mixed wet. This reduction in the hydraulic conductivity of the sand/sawdust cores mixed under saturated conditions repeated what occurred in the field in the original sawdust wall. This indicated that laboratory investigations can be a useful tool to highlight potential reductions in field hydraulic conductivities that may occur when differing materials are mixed under field conditions. 相似文献
4.
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. 相似文献
5.
We present an analytical solution of groundwater head response to tidal fluctuation in a coastal multilayered aquifer system consisting of an unconfined aquifer, a leaky confined aquifer and a semi‐permeable layer between them. The submarine outcrop of the confined aquifer is covered by a thin silt layer. A mathematical model and the analytical solution of this model are given. The silt layer reduces the amplitude of the hydraulic head fluctuation by a constant factor, and shifts the phase by a positive constant (time lag), both of which depend on the leakances of the silt layer and the semi‐permeable layer. The time lag is less than 1·5 h and 3·0 h for semi‐diurnal and diurnal sea tides respectively. When the leakance of the semi‐permeable layer or the silt layer assumes certain special values, the solution becomes the existing solutions derived by previous researchers. The amplitude of the hydraulic head fluctuation in the confined aquifer increases with the leakance of the silt layer and decreases with the leakance of the semi‐permeable layer, whereas the phase shift of the fluctuation decreases with both of them. A hypothetical example shows that neglecting the silt layer may result in significant parameter estimation discrepancy between the amplitude attenuation and the time‐lag fittings. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
6.
Zhixue Zhao Xiaoguang Wang Yonghong Hao Tongke Wang Abderrahim Jardani Herve Jourde Tian‐Chyi Jim Yeh Ming Zhang 《水文研究》2019,33(19):2551-2560
The study on the hydraulic properties of coastal aquifers has significant implications both in hydrological sciences and environmental engineering. Although many analytical solutions are available, most of them are based on the same basic assumption that assumes aquifers extend landward semi‐infinitely, which does not necessarily reflect the reality. In this study, the general solutions for a leaky confined coastal aquifer have been developed that consider both finitely landward constant‐head and no‐flow boundaries. The newly developed solutions were then used to examine theoretically the joint effects of leakage and aquifer length on hydraulic head fluctuations within the leaky confined aquifer, and the validity of using the simplified solution, which assumes the aquifer is semi‐infinite. The results illustrated that the use of the simplified solution may cause significant errors, depending on joint effects of leakage and aquifer length. A dimensionless characteristic parameter was then proposed as an index for judging the applicability of the simplified solution. In addition, practical application of the general solution for the constant‐head inland boundary was used to characterize the hydraulic properties of a leaky confined aquifer using the data collected from a field site at the Seine River estuary, France, and the versatility of the general solution was further justified. 相似文献
7.
Analyses of independent laboratory- and field-scale measurements from two sites on Sapelo Island, Georgia reveal heterogeneity in hydraulic parameters across the upland–estuary interface. Regardless of the method used (short-duration pumping tests, amplitude attenuation of tidal pumping data, sediment grain size distributions, and falling head permeameter tests), we obtain hydraulic conductivity of 10−4 m s−1 for the fine-grained, well-sorted, clean sands that make up the upland areas. Proximal to the upland–estuary boundary, the tidal pumping analyses and permeameter tests suggest that hydraulic conductivities decrease by more than two orders of magnitude, a result consistent with the presence of a clogging layer. Such a clogging layer may arise due to a variety of physical, chemical, or biological processes. The extent and orientation of the layers of reduced hydraulic conductivity near the upland–estuary boundary influence the nature of the aquifer's response to tidal forcing. Where the lower conductivity layer forms a relatively flat creek bank, tidal pumping produces a primarily mechanical response in the adjacent aquifer. Where the creek bank is nearly vertical, there is a more direct hydraulic connection between the tidal creek and the adjacent aquifer. The clogging layer likely contributes to the development of complicated flow pathways across the upland–estuary boundary. Effective flow paths calculated from tidal pumping data terminate within the marsh, beyond the boundary of the upland aquifer, suggesting a diffuse regime of groundwater discharge in the marsh. We postulate that, in many settings, submarsh flow may be as important as seepage faces for groundwater discharge into the marsh–estuary complex. 相似文献
8.
Darcian flow law in aquifers assumes that the aquifer hydraulic conductivity is constant and the groundwater movement is due only to the piezometric level changes through hydraulic gradient. In practice, after the well development the aquifer just around the well has comparatively larger hydraulic conductivity and gradient. Patchy aquifer solutions in the literature consider sudden hydraulic conductivity changes with distance for the steady state flow. The change of transmissivity is demonstrated by the application of slope‐matching procedure to actual field data. It is the main purpose of this paper to derive simple analytical expressions for aquifer parameter evaluations with steadily decreasing hydraulic conductivity around the well. Spatial nonlinear hydraulic conductivity changes around a large‐diameter well within the depression cone of a confined aquifer are considered as exponentially decreasing functions of the radial distance. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
9.
《水文科学杂志》2013,58(2):409-420
Abstract This work extends the algebraic expression of influence coefficients developed for one-dimensional aquifer models to a two-dimensional (2-D) case. First, the partial differential equation governing the flow in a 2-D semi-confined aquifer is discretized using a finite difference scheme. This results in a system of discrete equations presented in the form of water balance equations associated with a network of interconnected compartments centred on the grid nodes. The foregoing system is transformed into a series of uncoupled 1-D equations stated in terms of some generalized hydraulic head for which they are also solved. Second, the original hydraulic head is recovered from the generalized one via an appropriate linear transformation. Whence, the algebraic expression making the hydraulic head explicit versus sources and boundary conditions is derived. This discrete expression, mapped onto its continuous counterpart, helps to deduce an algebraic form of the inter-compartment influence coefficients. Finally, a comparison with the analytical Green function is carried out. 相似文献
10.
《水文科学杂志》2013,58(3):481-486
Abstract Explicit algebraic expressions of influence coefficients of an aquifer system are addressed. The methodology is straightforwardly exemplified using a one-dimensional (1D) semi-confined homogeneous aquifer model. It consists of the integration of the governing equation from two different viewpoints. On the one hand, the analytical solution to the governing equation is presented under an integral form involving the appropriate Green function. On the other hand, equations resulting from a finite difference discretization are processed to deduce a fully explicit algebraic expression of the nodal hydraulic head. This latter, arranged after the previously developed continuum solution, helps with deriving an explicit algebraic expression of influence coefficients corresponding to the studied aquifer model. Finally, the study is extended to an inhomogeneous aquifer model. 相似文献
11.
AbstractAn analytical solution is developed to delineate the capture zone of a pumping well in an aquifer with a regional flow perpendicular to a stream, assuming a leaky layer between the stream and the aquifer. Three different scenarios are considered for different pumping rates. At low pumping rates, the capture zone boundary will be completely contained in the aquifer. At medium pumping rates, the tip of the capture zone boundary will intrude into the leaky layer. Under these two scenarios, all the pumped water is supplied from the regional groundwater flow in the aquifer. At high pumping rates, however, the capture zone boundary intersects the stream and pumped water is supplied from both the aquifer and the stream. The two critical pumping rates which separate these three scenarios, as well as the proportion of pumped water from the stream and the aquifer, are determined for different hydraulic settings.Editor D. Koutsoyiannis; Associate editor A. KoussisCitation Asadi-Aghbolaghi, M., Rakhshandehroo, G.R., and Kompani-Zare, M., 2013. An analytical approach to capture zone delineation for a well near a stream with a leaky layer. Hydrological Sciences Journal, 58 (8), 1813–1823. 相似文献
12.
This paper addresses the issue of flow in heterogeneous leaky confined aquifers subject to leakage. The leakage into the confined aquifer is driven by spatial and periodic fluctuations of water table in an overlying phreatic aquifer. The introduction of leakage leads to non-uniformity in the mean head gradient and results in nonstationarity in hydraulic head and velocity fields. Therefore, a nonstationary spectral approach based on Fourier–Stieltjes representations for the perturbed quantities is adopted to account for the spatial variability of nonstationary head fields. Closed-form expressions for the variances of hydraulic head and specific discharge are developed in terms of statistical properties of hydraulic parameters. The results indicate that the spatiotemporal variations in leakage leads to enhanced variability of the hydraulic head and of the specific discharge, which increase with distance from any arbitrary reference point. The coefficient of leakage and the spatial structure of log transmissivity field and of the amplitude of water table fluctuation are critical in quantifying the variability of the hydraulic head and of the specific discharge. 相似文献
13.
Salihu Mustafa 《Advances in water resources》1984,7(2):81-84
Flow situation into a ditch resulting from an aquifer that is semi-confined and leaky is investigated. Theoretical solution employing Fourier Series is given and then compared with that obtained by the Laplace Transformation. 相似文献
14.
The HySuf‐FEM code (Hydrodynamic of Subsurface Flow by Finite Element Method) is proposed in this article in order to estimate the spatial variability of the transmissivity values of the Berrechid aquifer (Morocco). The calibration of the model is based on the hydraulic head, hydraulic conductivity and recharge. Three numerical tests are used to validate the model and verify its convergence. The first test case consists in using the steady analytical solution of the Poisson equation. In the second, the model has been compared with the hydrogeological system which is characterized by an unconfined monolayer (isotropic layer) and computed by using PMWIN‐MODFLOW software. The third test case is based on the comparison between the results of HySuf‐FEM and the multiple cell balance method in the aquifer system with natural boundaries case. Good agreement between the Hydrodynamic of Subsurface Flow, the numerical tests and the spatial distribution of the thickening of the hydrogeological system is deduced from the analysis and the interpretations of hydrogeological wells. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
15.
《水文科学杂志》2013,58(6)
Abstract Abstract Characterization of heterogeneity at the field scale generally requires detailed aquifer properties such as transmissivity and hydraulic head. An accurate delineation of these properties is expensive and time consuming, and for many if not most groundwater systems, is not practical. As an alternative approach, stochastic representation of random fields is used and presented in this paper. Specifically, an iterative stochastic conditional simulation approach was applied to a hypothetical and highly heterogeneous pre-designed aquifer system. The approach is similar to the classical co-kriging technique; it uses a linear estimator that depends on the covariance functions of transmissivity (T), and hydraulic head (h), as well as their cross-covariances. A linearized flow equation along with a conditional random field generator constitutes the iterative process of the conditional simulation. One hundred equally likely realizations of transmissivity fields with pre-specified geostatistical parameters were generated, and conditioned to both limited transmissivity and head data. The successful implementation of the approach resulted in conditioned flow paths and travel-time distribution under different degrees of aquifer heterogeneity. This approach worked well for fields exhibiting small variances. However, for random fields exhibiting large variances (greater than 1.0), an iterative procedure was used. The results show that, as the variance of the ln[T] increases, the flow paths tend to diverge, resulting in a wide spectrum of flow conditions, with no direct discernable relationship between the degree of heterogeneity and travel time. The applied approach indicates that high errors may result when estimation of particle travel times in a heterogeneous medium is approximated by an equivalent homogeneous medium. 相似文献
16.
The hydraulic gradient between aquifers and rivers is one of the most variable properties in a river/aquifer system. Detailed process understanding of bank storage under hydraulic gradients is obtained from a two‐dimensional numerical model of a variably saturated aquifer slice perpendicular to a river. Exchange between the river and the aquifer occurs first at the interface with the unsaturated zone. The proportion of total water exchanged through the river bank compared to the river bed is a function of aquifer hydraulic conductivity, partial penetration, and hydraulic gradient. Total exchange may be estimated to within 50% using existing analytical solutions provided that unsaturated zone processes do not strongly influence exchange. Model‐calculated bank storage is at a maximum when no hydraulic gradient is present and increases as the hydraulic conductivity increases. However, in the presence of a hydraulic gradient, the largest exchange flux or distance of penetration does not necessarily correspond to the highest hydraulic conductivity, as high hydraulic conductivity increases the components of exchange both into and out of an aquifer. Flood wave characteristics do not influence ambient groundwater discharge, and so in large floods, hydraulic gradients must be high to reduce the volume of bank storage. Practical measurement of bank storage metrics is problematic due to the limitations of available measurement technologies and the nested processes of exchange that occur at the river‐aquifer interface. Proxies, such as time series concentration data in rivers and groundwater, require further development to be representative and quantitative. 相似文献
17.
In this second of two papers, analytical step-response functions, developed in the companion paper for several cases of transient hydraulic interaction between a fully penetrating stream and a confined, leaky, or water-table aquifer, are used in the convolution integral to calculate aquifer heads, streambank seepage rates, and bank storage that occur in response to stream-stage fluctuations and basinwide recharge or evapotranspiration. Two computer programs developed on the basis of these step-response functions and the convolution integral are applied to the analysis of hydraulic interaction of two alluvial stream–aquifer systems in the northeastern and central United States. These applications demonstrate the utility of the analytical functions and computer programs for estimating aquifer and streambank hydraulic properties, recharge rates, streambank seepage rates, and bank storage. Analysis of the water-table aquifer adjacent to the Blackstone River in Massachusetts suggests that the very shallow depth of water table and associated thin unsaturated zone at the site cause the aquifer to behave like a confined aquifer (negligible specific yield). This finding is consistent with previous studies that have shown that the effective specific yield of an unconfined aquifer approaches zero when the capillary fringe, where sediment pores are saturated by tension, extends to land surface. Under this condition, the aquifer's response is determined by elastic storage only. Estimates of horizontal and vertical hydraulic conductivity, specific yield, specific storage, and recharge for a water-table aquifer adjacent to the Cedar River in eastern Iowa, determined by the use of analytical methods, are in close agreement with those estimated by use of a more complex, multilayer numerical model of the aquifer. Streambank leakance of the semipervious streambank materials also was estimated for the site. The streambank-leakance parameter may be considered to be a general (or lumped) parameter that accounts not only for the resistance of flow at the river–aquifer boundary, but also for the effects of partial penetration of the river and other near-stream flow phenomena not included in the theoretical development of the step-response functions. 相似文献
18.
Abstract This work makes explicit an algebraic expression giving the matrix of transient influence coefficients associated with a one-dimensional semi-confined aquifer model. The domain studied is divided into a series of connected and completely mixed compartments over which the governing equation is discretized. The discrete equations obtained are solved for the compartmental hydraulic head and used to derive the algebraic expression in question. The basic properties of the so-called algebraic influence coefficients are investigated. In particular, their consistency with the exact Green function is highlighted. Finally, the newly derived influence coefficients are applied to a simplified aquifer system in order to formulate and solve the problem of identifying illegal groundwater pumping. 相似文献
19.
Hydraulic diffusivity and macrodispersivity calculations embedded in a geographic information system
F.A.L. Pacheco 《水文科学杂志》2013,58(4):930-944
Abstract A numerical technique is presented whereby aquifer hydraulic diffusivities (D) and macrodispersivities (α) are calculated by linear equations rewritten from flow and solute transport differential equations. The approach requires a GIS to calculate spatial and temporal hydraulic head (h) and solute concentration gradients. The model is tested in Portugal, in a semi-confined aquifer periodically monitored for h and chloride/sulphate concentrations. Average D (0.46 m2/s) and α (1975 m) compare favourably with literature results. The relationship between α and scale (L) is also investigated. In this context, two aquifer groups could be identified: the first group is heterogeneous at the “macroscopic” scale (solute travelled distances <1 km), but homogeneous at the “megascopic” scale. The overall scale dependency in this case is given by an equation of logarithmic type. The second group is heterogeneous at the macroscopic and megascopic scales, with a scale dependency of linear type. Citation Pacheco, F.A.L., 2013. Hydraulic diffusivity and macrodispersivity calculations embedded in a geographic information system. Hydrological Sciences Journal, 58 (4), 930–944. 相似文献
20.
Analytic solutions are presented for steady interface flow in aquifers consisting of a confined and a semi-confined section. The total discharge is constant in the confined section and is directed towards the semi-confined section, which is bounded on top by a leaky layer that separates the aquifer from the sea. All solutions are based on the Dupuit approximation. The interface position is fully determined by two dimensionless parameters. The first parameter is the product of the uniform gradient towards the sea and the dimensionless leakage factor divided by the dimensionless density difference. The second dimensionless parameter is the length of the semi-confined section divided by the leakage factor. Four types of flow are distinguished. For cases I and II the interface does not reach the end of the semi-confined section, while it does for cases III and IV. For cases I and III the interface extends into the confined section, while for cases II and IV it remains entirely in the semi-confined section. Analytic solutions are presented for the position of the interface for all four cases. Diagrams are presented to determine the type of flow and the position of the interface toe based on the two dimensionless parameters. The pre-development position of the interface along the Georgia coast is computed as an illustration of the practical application of the presented formulas. Another practical application is the estimation of how far a numerical model of seawater intrusion should be extended into the sea for accurate simulations. 相似文献