首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 484 毫秒
1.
  相似文献   

2.
Within the framework of stochastic theory and the spectral perturbation techniques, three-dimensional dispersion in partially saturated soils with fractal log hydraulic conductivity distribution is analyzed. Our analysis is focused on the impact of fractal dimension of log hydraulic conductivity distribution, local dispersivity, and unsaturated flow parameters, such as the soil poresize distribution parameter and the moisture distribution parameter, on the spreading behavior of solute plume and the concentration variance. Approximate analytical solutions to the stochastic partial differential equations are derived for the variance of asymptotic solute concentration and asymptotic macrodispersivities.  相似文献   

3.
In this study, we derive analytical solutions of the first two moments (mean and variance) of pressure head for one-dimensional steady state unsaturated flow in a randomly heterogeneous layered soil column under random boundary conditions. We first linearize the steady state unsaturated flow equations by Kirchhoff transformation and solve the moments of the transformed variable up to second order in terms of σY and σβ, the standard deviations of log hydraulic conductivity Y=ln(Ks) and of the log pore size distribution parameter β=ln(α). In addition, we also give solutions for the mean and variance of the unsaturated hydraulic conductivity. The analytical solutions of moment equations are validated via Monte Carlo simulations.  相似文献   

4.
An excellent tool for checking numerical models of unsaturated flow in groundwater is analytical solutions. However, because of the highly nonlinear nature of the governing partial differential equation, only a limited number of analytical solutions are available. This paper first gives some simple 1-D solutions. Next, by use of a transformation, the nonlinear partial differential equation is converted to a linear one for a specific form of the moisture content vs. pressure head and relative hydraulic conductivity vs. pressure head curves. This allows both 2-D and 3-D solutions to be derived, which is done in this paper. Finally, computations from a finite element computer program are compared with results from one of the analytical solutions to illustrate the use of the derived equations.  相似文献   

5.
Spatiotemporal variations of groundwater level due to a white noise recharge time series and a random transmissivity field in a bounded unconfined aquifer was studied. The analytical solutions for the variance and covariance of groundwater level were derived with non-stationary spectral analyses and superposition principle. It was found that the fluctuations of groundwater level are spatially non-stationary due to a fixed head boundary condition and temporal non-stationary at early time but gradually became stationary as time progresses due to effect of the initial condition. The variation in groundwater level is mainly caused by the random source/sink in the case of temporally random recharge and spatially random transmissivity. The effect of heterogeneity is to increase the variation of groundwater level and the maximum effect occurs close to the constant head boundary because of the linear mean hydraulic gradient. The heterogeneity also enhances the correlation of groundwater level, especially at large time intervals and small spatial distances.  相似文献   

6.
The unconditional stochastic studies on groundwater flow and solute transport in a nonstationary conductivity field show that the standard deviations of the hydraulic head and solute flux are very large in comparison with their mean values (Zhang et al. in Water Resour Res 36:2107–2120, 2000; Wu et al. in J Hydrol 275:208–228, 2003; Hu et al. in Adv Water Resour 26:513–531, 2003). In this study, we develop a numerical method of moments conditioning on measurements of hydraulic conductivity and head to reduce the variances of the head and the solute flux. A Lagrangian perturbation method is applied to develop the framework for solute transport in a nonstationary flow field. Since analytically derived moments equations are too complicated to solve analytically, a numerical finite difference method is implemented to obtain the solutions. Instead of using an unconditional conductivity field as an input to calculate groundwater velocity, we combine a geostatistical method and a method of moment for flow to conditionally simulate the distributions of head and velocity based on the measurements of hydraulic conductivity and head at some points. The developed theory is applied in several case studies to investigate the influences of the measurements of hydraulic conductivity and/or the hydraulic head on the variances of the predictive head and the solute flux in nonstationary flow fields. The study results show that the conditional calculation will significantly reduce the head variance. Since the hydraulic head measurement points are treated as the interior boundary (Dirichlet boundary) conditions, conditioning on both the hydraulic conductivity and the head measurements is much better than conditioning only on conductivity measurements for reduction of head variance. However, for solute flux, variance reduction by the conditional study is not so significant.  相似文献   

7.
As is well known, a complete stochastic solution of the stochastic differential equation governing saturated groundwater flow leads to an infinite hierarchy of equations in terms of higher-order moments. Perturbation techniques are commonly used to close this hierarchy, using power-series expansions. These methods are applied by truncating the series after a finite number of terms, and products of random gradients of conductivity and head potential are neglected. Uncertainty regarding the number or terms required to yield a sufficiently accurate result is a significant drawback with the application of power series-based perturbation methods for such problems. Low-order series truncation may be incapable of representing fundamental characteristics of flow and can lead to physically unreasonable and inaccurate solutions of the stochastic flow equation. To support this argument, one-dimensional, steady-state, saturated groundwater flow is examined, for the case of a spatially distributed hydraulic conductivity field. An ordinary power-series perturbation method is used to approximate the mean head, using second-order statistics to characterize the conductivity field. Then an interactive perturbation approach is introduced, which yields improved results compared to low-order, power-series perturbation methods for situations where strong interactions exist between terms in such approximations. The interactive perturbation concept is further developed using Feynman-type diagrams and graph theory, which reduce the original stochastic flow problem to a closed set of equations for the mean and the covariance functions. Both theoretical and practical advantages of diagrammatic solutions are discussed; these include the study of bounded domains and large fluctuations.  相似文献   

8.
We present an efficient numerical method for solving stochastic porous media flow problems. Single-phase flow with a random conductivity field is considered in a standard first-order perturbation expansion framework. The numerical scheme, based on finite element techniques, is computationally more efficient than traditional approaches because one can work with a much coarser finite element mesh. This is achieved by avoiding the common finite element representation of the conductivity field. Computations with the random conductivity field only arise in integrals of the log conductivity covariance function. The method is demonstrated in several two- and three-dimensional flow situations and compared to analytical solutions and Monte Carlo simulations. Provided that the integrals involving the covariance of the log conductivity are computed by higher-order Gaussian quadrature rules, excellent results can be obtained with characteristic element sizes equal to about five correlation lengths of the log conductivity field. Investigations of the validity of the proposed first-order method are performed by comparing nonlinear Monte Carlo results with linear solutions. In box-shaped domains the log conductivity standard deviation σY may be as large as 1.5, while the head variance is considerably influenced by nonlinear effects as σY approaches unity in more general domains.  相似文献   

9.
This paper describes a stochastic analysis of steady state flow in a bounded, partially saturated heterogeneous porous medium subject to distributed infiltration. The presence of boundary conditions leads to non-uniformity in the mean unsaturated flow, which in turn causes non-stationarity in the statistics of velocity fields. Motivated by this, our aim is to investigate the impact of boundary conditions on the behavior of field-scale unsaturated flow. Within the framework of spectral theory based on Fourier–Stieltjes representations for the perturbed quantities, the general expressions for the pressure head variance, variance of log unsaturated hydraulic conductivity and variance of the specific discharge are presented in the wave number domain. Closed-form expressions are developed for the simplified case of statistical isotropy of the log hydraulic conductivity field with a constant soil pore-size distribution parameter. These expressions allow us to investigate the impact of the boundary conditions, namely the vertical infiltration from the soil surface and a prescribed pressure head at a certain depth below the soil surface. It is found that the boundary conditions are critical in predicting uncertainty in bounded unsaturated flow. Our analytical expression for the pressure head variance in a one-dimensional, heterogeneous flow domain, developed using a nonstationary spectral representation approach [Li S-G, McLaughlin D. A nonstationary spectral method for solving stochastic groundwater problems: unconditional analysis. Water Resour Res 1991;27(7):1589–605; Li S-G, McLaughlin D. Using the nonstationary spectral method to analyze flow through heterogeneous trending media. Water Resour Res 1995; 31(3):541–51], is precisely equivalent to the published result of Lu et al. [Lu Z, Zhang D. Analytical solutions to steady state unsaturated flow in layered, randomly heterogeneous soils via Kirchhoff transformation. Adv Water Resour 2004;27:775–84].  相似文献   

10.
In this article, we are concerned with the statistics of steady unsaturated flow in soils with a fractal hydraulic conductivity distribution. It is assumed that the spatial distribution of log hydraulic conductivity can be described as an isotropic stochastic fractal process. The impact of the fractal dimension of this process, the soil pore-size distribution parameter, and the characteristic length scale on the variances of tension head and the effective conductivity is investigated. Results are obtained for one-dimensional and three-dimensional flows. Our results indicate that the tension head variance is scale-dependent for fractal distribution of hydraulic conductivity. Both tension head variance and effective hydraulic conductivity depend strongly on the fractal dimension. The soil pore-size distribution parameter is important in reducing the variability of the unsaturated hydraulic conductivity and of the fluxes.  相似文献   

11.
In this study, the KLME approach, a moment-equation approach based on the Karhunen–Loeve decomposition developed by Zhang and Lu (Comput Phys 194(2):773–794, 2004), is applied to unconfined flow with multiple random inputs. The log-transformed hydraulic conductivity F, the recharge R, the Dirichlet boundary condition H, and the Neumann boundary condition Q are assumed to be Gaussian random fields with known means and covariance functions. The F, R, H and Q are first decomposed into finite series in terms of Gaussian standard random variables by the Karhunen–Loeve expansion. The hydraulic head h is then represented by a perturbation expansion, and each term in the perturbation expansion is written as the products of unknown coefficients and Gaussian standard random variables obtained from the Karhunen–Loeve expansions. A series of deterministic partial differential equations are derived from the stochastic partial differential equations. The resulting equations for uncorrelated and perfectly correlated cases are developed. The equations can be solved sequentially from low to high order by the finite element method. We examine the accuracy of the KLME approach for the groundwater flow subject to uncorrelated or perfectly correlated random inputs and study the capability of the KLME method for predicting the head variance in the presence of various spatially variable parameters. It is shown that the proposed numerical model gives accurate results at a much smaller computational cost than the Monte Carlo simulation.  相似文献   

12.
Steady interface flow in heterogeneous aquifer systems is simulated with single‐density groundwater codes by using transformed values for the hydraulic conductivity and thickness of the aquifers and aquitards. For example, unconfined interface flow may be simulated with a transformed model by setting the base of the aquifer to sea level and by multiplying the hydraulic conductivity with 41 (for sea water density of 1025 kg/m3). Similar transformations are derived for unconfined interface flow with a finite aquifer base and for confined multi‐aquifer interface flow. The head and flow distribution are identical in the transformed and original model domains. The location of the interface is obtained through application of the Ghyben‐Herzberg formula. The transformed problem may be solved with a single‐density code that is able to simulate unconfined flow where the saturated thickness is a linear function of the head and, depending on the boundary conditions, the code needs to be able to simulate dry cells where the saturated thickness is zero. For multi‐aquifer interface flow, an additional requirement is that the code must be able to handle vertical leakage in situations where flow in an aquifer is unconfined while there is also flow in the aquifer directly above it. Specific examples and limitations are discussed for the application of the approach with MODFLOW. Comparisons between exact interface flow solutions and MODFLOW solutions of the transformed model domain show good agreement. The presented approach is an efficient alternative to running transient sea water intrusion models until steady state is reached.  相似文献   

13.
The impact of groundwater withdrawal on surface water is a concern of water users and water managers, particularly in the arid western United States. Capture maps are useful tools to spatially assess the impact of groundwater pumping on water sources (e.g., streamflow depletion) and are being used more frequently for conjunctive management of surface water and groundwater. Capture maps have been derived using linear groundwater flow models and rely on the principle of superposition to demonstrate the effects of pumping in various locations on resources of interest. However, nonlinear models are often necessary to simulate head‐dependent boundary conditions and unconfined aquifers. Capture maps developed using nonlinear models with the principle of superposition may over‐ or underestimate capture magnitude and spatial extent. This paper presents new methods for generating capture difference maps, which assess spatial effects of model nonlinearity on capture fraction sensitivity to pumping rate, and for calculating the bias associated with capture maps. The sensitivity of capture map bias to selected parameters related to model design and conceptualization for the arid western United States is explored. This study finds that the simulation of stream continuity, pumping rates, stream incision, well proximity to capture sources, aquifer hydraulic conductivity, and groundwater evapotranspiration extinction depth substantially affect capture map bias. Capture difference maps demonstrate that regions with large capture fraction differences are indicative of greater potential capture map bias. Understanding both spatial and temporal bias in capture maps derived from nonlinear groundwater flow models improves their utility and defensibility as conjunctive‐use management tools.  相似文献   

14.
 3D groundwater flow at the fractured site of Asp? (Sweden) is simulated. The aim was to characterise the site as adequately as possible and to provide measures on the uncertainty of the estimates. A stochastic continuum model is used to simulate both groundwater flow in the major fracture planes and in the background. However, the positions of the major fracture planes are deterministically incorporated in the model and the statistical distribution of the hydraulic conductivity is modelled by the concept of multiple statistical populations; each fracture plane is an independent statistical population. Multiple equally likely realisations are built that are conditioned to geological information on the positions of the major fracture planes, hydraulic conductivity data, steady state head data and head responses to six different interference tests. The experimental information could be reproduced closely. The results of the conditioning are analysed in terms of ensemble averaged average fracture plane conductivities, the ensemble variance of average fracture plane conductivities and the statistical distribution of the hydraulic conductivity in the fracture planes. These results are evaluated after each conditioning stage. It is found that conditioning to hydraulic head data results in an increase of the hydraulic conductivity variance while the statistical distribution of log hydraulic conductivity, initially Gaussian, becomes more skewed for many of the fracture planes in most of the realisations.  相似文献   

15.
Reply     
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.  相似文献   

16.
Non-local stochastic moment equations are used successfully to analyze groundwater flow in randomly heterogeneous media. Here we present a moment equations-based approach to quantify the uncertainty associated with the estimation of well catchments. Our approach is based on the development of a complete second order formalism which allows obtaining the first statistical moments of the trajectories of conservative solute particles advected in a generally non-uniform groundwater flow. Approximate equations of moments of particles’ trajectories are then derived on the basis of a second order expansion in terms of the standard deviation of the aquifer log hydraulic conductivity. Analytical expressions are then obtained for the predictors of locations of mean stagnation points, together with their associated uncertainties. We implement our approach on heterogeneous media in bounded two-dimensional domains, with and without including the effect of conditioning on hydraulic conductivity information. The impact of domain size, boundary conditions, heterogeneity and non-stationarity of hydraulic conductivity on the prediction of a well catchment is explored. The results are compared against Monte Carlo simulations and semi-analytical solutions available in the literature. The methodology is applicable to both infinite and bounded domains and is free of distributional assumptions (and so applies to both Gaussian and non-Gaussian log hydraulic conductivity fields) and formally includes the effect of conditioning on available information.  相似文献   

17.
In this work, a stochastic methodology is applied to analyze the variability of the poroelastic response of the heterogeneous medium at the field scale. To solve the problem analytically, we restrict our attention to the one-dimensional models, where fluid flow as well as deformation occurs in one direction only under a constant applied stress. Assuming statistic homogeneity, the closed-form solutions that describe the variability of fluid pressure head, and a solid's strain and displacement are developed using a spectral approach based on Fourier–Stieltjes representations for the perturbed quantities. The influence of the correlation length of the log hydraulic conductivity on these results is investigated. It is found that the variances of the solid's strain and displacement increase with the correlation length of the log hydraulic conductivity, while the correlation length of the log hydraulic conductivity plays the role in reducing the variability of the specific discharge.  相似文献   

18.
《水文科学杂志》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.  相似文献   

19.
The groundwater variations in unconfined aquifers are governed by the nonlinear Boussinesq's equation. Analytical solution for groundwater fluctuations in coastal aquifers under tidal forcing can be solved using perturbation methods. However, the perturbation parameters should be properly selected and predefined for traditional perturbation methods. In this study, a new dimensional, higher‐order analytical solution for groundwater fluctuations is proposed by using the homotopy perturbation method with a virtual perturbation parameter. Parameter‐expansion method is used to remove the secular terms generated during the solution process. The solution does not require any predefined perturbation parameter and valid for higher values of amplitude parameter A/D, where A is the amplitude of the tide and D is the aquifer thickness.  相似文献   

20.
Assuming that the ln hydraulic conductivity in an aquifer is mathematically approximated by a spatial deterministic surface, or trend, plus a stationary random noise, we treat the problem of finding what the effective hydraulic conductivity of that aquifer is. This problem is tackled by spectral methods applied to a type of diffusion equation of groundwater flow, together with suitable coordinate transformations. Analytical (exact) solutions in terms of elementary functions are presented for one- and three-dimensional finite and infinite domains. Stability criteria are obtained for the solutions, in terms of a critical parameter, that turns out to involve the product of correlation scale and trend gradient. For the case of finite and symmetrical domains, additional provisions to insure the stability of numerical calculations of effective hydraulic conductivity are provided. Effective hydraulic conductivity is an important property, with potential applications in the calibrations of groundwater and transport numerical models.  相似文献   

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

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