共查询到20条相似文献，搜索用时 62 毫秒
1.
《Advances in water resources》2002,25(7):733-746
Backward location and travel time probabilities, which provide information about the former location of contamination in an aquifer, can be used to identify unknown contamination sources. Backward location probability describes the possible upgradient positions of contamination at a known time in the past, and backward travel time probability describes the time required for contamination to travel from a known upgradient location to an observation point. These probabilities are related to adjoint states of resident concentration, and their governing equation is the adjoint of a forward contaminant transport model. Using adjoint theory to obtain the appropriate governing equation, we extend the backward probability model for conservative solutes to more general non-uniform and transient flow fields. In particular, we address three important extensions, spatially-varying porosity, transient flow and temporally-varying porosity, and internal distributed sources and sinks of solute and water. For the first time we learn that forward and backward location and travel time probabilities are not necessarily equivalent to adjoint states, but are related to them. The extensions are illustrated using a vertically-integrated groundwater model, creating transient flow by a step change in pumping and using areal recharge as an internal distributed source. Both the movement and spread of probabilities are affected. With internal sources of water, there are two interpretations of backward probability, depending on whether or not the source of water is also a source of solute. The results demonstrate how the backward probability model can be applied to other, perhaps more important, non-uniform and transient flow conditions, with time- and space-varying water storage, such as time-varying pumping or unsaturated (or saturated–unsaturated) flow and transport with spatially- and temporally-varying moisture content. 相似文献
2.
Fausto Cupola Maria Giovanna Tanda Andrea Zanini 《Stochastic Environmental Research and Risk Assessment (SERRA)》2015,29(1):169-182
Inverse methods can be used to recover the pollutant source location from concentration data. In this paper, the relative effectiveness of two proposed methods, simultaneous release function and source location identification (SRSI) and backward probability model based on adjoint state method (BPM-ASM) are evaluated using real data collected by using experimental equipment. The device is a sandbox that reproduces an unconfined aquifer in which all the variables are controlled. A numerical model was calibrated using experimental observations. The SRSI is a stochastic procedure which finds the source location and the release history by means of a Bayesian geostatistical approach (GA). The BPM-ASM provides the backward probability location of the pollutant detected at a monitoring point by means of a reverse transport simulation. The results show that both methods perform well. While the simultaneous release function and SRSI method requires a preliminary delineation of a probable source area and some weak hypotheses about the statistical structure of the unknown release function, the backward probability model requires some hypothesis about the contaminant release time. A case study was performed using two observation points only, and despite the scarcity of data, both methodologies were able to accurately reconstruct the true source location. The GA has the advantage to recover the release history function too, whilst the backward probability model works well with fewer data. If there are many observations, both methodologies may be computationally heavy. A transfer function approach has been adopted for the numerical definition of the sensitivity matrix in the SRSI method. The reliability of the experimental equipment was tested in previous laboratory works, conducted under several different conditions. 相似文献
3.
Identification of contaminant point source in surface waters based on backward location probability density function method 总被引：2，自引：0，他引：2
A backward location probability density function (BL-PDF) method capable of identifying location of point sources in surface waters is presented in this paper. The relation of forward location probability density function (FL-PDF) and backward location probability density, based on adjoint analysis, is validated using depth-averaged free-surface flow and mass transport models and several surface water test cases. The solutions of the backward location PDF transport equation agreed well to the forward location PDF computed using the pollutant concentration at the monitoring points. Using this relation and the distribution of the concentration detected at the monitoring points, an effective point source identification method is established. The numerical error of the backward location PDF simulation is found to be sensitive to the irregularity of the computational meshes, diffusivity, and velocity gradients. The performance of identification method is evaluated regarding the random error and number of observed values. In addition to hypothetical cases, a real case was studied to identify the source location where a dye tracer was instantaneously injected into a stream. The study indicated the proposed source identification method is effective, robust, and quite efficient in surface waters; the number of advection–diffusion equations needed to solve is equal to the number of observations. 相似文献
4.
Contaminant transport models under random sources 总被引：1，自引：0，他引：1
5.
《Advances in water resources》1996,19(5):297-315
The contamination of groundwater by various hazardous materials has emerged as a primary environmental issue. The pollution of oil reservoirs is a closely related problem in that microorganisms are involved in the contaminant process. The mathematical models that describe these phenomena involve a set of nonlinear advective-diffusive-reactive transport equations, which may involve reactions with all the species and are themselves coupled to growth equations for the subsurface bacterial population. In this article, we discuss and compare different mathematical models, present Eulerian-Lagrangian localized adjoint methods (ELLAM) and combine them with specific linearization techniques to solve these nonlinear transport systems. The derived numerical schemes systematically adapt to the changing features of governing equations. The relative importance of advection, diffusion and reaction is directly incorporated into the schemes by judicious choice of the test functions in the variational formulations. Numerical experiments are presented to show the potential of these methods. 相似文献
6.
ABSTRACT Forward–backward solute dispersion with an intermediate point source in one-dimensional semi-infinite homogeneous porous media is studied in this paper. Solute transport under sorption conditions, first-order decay and zero-order production terms are included. The first type of boundary condition is taken as a constant point source at an intermediate point from where forward and backward solute dispersion is examined. The Laplace transform method is adopted to solve the governing equation analytically. All the analytical results are obtained in graphical form to investigate the forward–backward solute transport in porous media for various hydrological input data. The graphical nature of the analytical solution is compared with numerical data taken from existing literature and similar results are obtained. Also, numerical solution of the governing equation is obtained by the Crank-Nicolson finite difference scheme and validated with the analytical solution, which demonstrates good agreement between them. Accuracy of the solution is also observed by using RMSE. 相似文献
7.
Daniel J. Goode 《Ground water》1992,30(2):257-261
Abstract. During unsteady or transient ground-water flow, the fluid mass per unit volume of aquifer changes as the potentiometric head changes, and solute transport is affected by this change in fluid storage. Three widely applied numerical models of two-dimensional transport partially account for the effects of transient flow by removing terms corresponding to the fluid continuity equation from the transport equation, resulting in a simpler governing equation. However, fluid-storage terms remaining in the transport equation that change during transient flow are, in certain cases, held constant in time in these models. For the case of increasing heads, this approximation, which is unacknowledged in these models'documentation, leads to transport velocities that are too high, and increased concentration at fluid and solute sources. If heads are dropping in time, computed transport velocities are too low. Using parameters that somewhat exaggerate the effects of this approximation, an example numerical simulation indicates solute travel time error of about 14 percent but only minor errors due to incorrect dilution volume. For horizontal flow and transport models that assume fluid density is constant, the product of porosity and aquifer thickness changes in time: initial porosity times initial thickness plus the change in head times the storage coefficient. This formula reduces to the saturated thickness in unconfined aquifers if porosity is assumed to be constant and equal to specific yield. The computational cost of this more accurate representation is insignificant and is easily incorporated in numerical models of solute transport. 相似文献
8.
Laboratory and numerical investigation of flow and transport near a seepage-face boundary 总被引：3，自引：0，他引：3
Laboratory and numerical modeling investigations were completed to study the unconfined ground water flow and transport processes near a seepage-face boundary. The laboratory observations were made in a radial sand tank and included measurements of the height of the seepage face, flow velocity near the seepage face, travel time distribution of multiple tracer slugs, and streamlines. All the observations were reliably reproduced with a three-dimensional, axi-symmetric, variably saturated ground water flow model. Physical data presented in this work demonstrate and quantify the importance of three-dimensional transport patterns within a seepage-face zone. The results imply that vertically averaged flow models that employ Dupuit approximations might introduce error in the analysis of localized solute transport near a seepage-face boundary. The experimental dataset reported in this work will also be of interest for those who are attempting to validate a numerical algorithm for solving ground water and contaminant discharge patterns near a surface-water boundary. 相似文献
9.
The concept of vulnerability of drinking water sources is reviewed, and a quantitative approach for assessing well vulnerability for complex three-dimensional ground water systems is developed. The approach focuses on the relative expected impact of potential contaminant sources at unknown locations within a well capture zone, providing relative measures of intrinsic well vulnerability, including the expected times of arrival of a contaminant, the dispersion-related reduction in concentration, the time taken to breach a certain quality objective, and the corresponding exposure times. Thus, the result of the analysis includes the usual advective travel time information used in conventional wellhead protection analysis, plus a set of selected quantitative measures expressing the expected impact. The technique is based on adjoint theory and combines forward- and backward-in-time transport modeling using a standard numerical flow and transport code. The methodology is demonstrated using the case study of a complex glacial multiaquifer system in Ontario. The new approach will be useful in helping water managers develop more physically based and quantitative wellhead protection strategies. 相似文献
10.
Rubin Y. Cushey M. A. Bellin A. 《Stochastic Environmental Research and Risk Assessment (SERRA)》1994,8(1):57-77
This paper presents the principles underlying a recently developed numerical technique for modeling transport in heterogeneous porous media. The method is then applied to derive the concentration mean and variance, the concentration CDF, exceedance probabilities and exposure time CDF, which are required by various regulatory agencies for risk and performance assessment calculations. The dependence of the various statistics on elapsed travel time, location in space, the dimension of the detection volume, natural variability and pore-scale dispersion is investigated and discussed. 相似文献
11.
Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation 
下载免费PDF全文
下载免费PDF全文 In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field. 相似文献
12.
Solute transport is usually modeled by the advection-dispersion-reaction equation. In the standard approach, mechanical dispersion is a tensor with principal directions parallel and perpendicular to the flow vector. Since realistic scenarios include nonuniform and unsteady flow fields, the governing equation has full tensor mechanical dispersion. When conventional grid-based numerical methods are used, approximation of the cross terms arising from the off-diagonal terms cause nonphysical solution with oscillations. As an example, for the common scenario of contaminant input into a domain with zero initial concentration, the cross-dispersion terms can result in negative concentrations that can wreak havoc in reactive transport applications. To address this issue, we use the well-known flux-corrected-transport (FCT) technique for a standard finite volume method. Although FCT has most often been used to eliminate oscillations resulting from discretization of the advection term for explicit time stepping, we show that it can be adapted for full-tensor dispersion and implicit time stepping. Unlike other approaches based on new discretization techniques (e.g., mimetic finite difference, nonlinear finite volume), FCT has the advantage of being flexible and widely applicable. Implementation of FCT requires solving an additional system of equations at each time step, using a modified “low order” matrix and a modified right-hand-side vector. To demonstrate the flexibility of FCT, we have modified the well-known and widely used groundwater solute transport simulator, MT3DMS. We apply the new simulator, MT3DMS-FCT, to several benchmark problems that suffer from negative concentrations when using MT3DMS. The new results are mass conservative and strictly nonnegative. 相似文献
13.
F.A. Radu N. Suciu J. HoffmannA. Vogel O. Kolditz C.-H. ParkS. Attinger 《Advances in water resources》2011,34(1):47-61
This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater. 相似文献
14.
《Advances in water resources》1998,22(2):103-116
We consider a one-dimensional model biodegradation system consisting of two reaction–advection equations for nutrient and pollutant concentrations and a rate equation for biomass. The hydrodynamic dispersion is ignored. Under an explicit condition on the decay and growth rates of biomass, the system can be approximated by two component models by setting biomass kinetics to equilibrium. We derive closed form solutions for constant speed traveling fronts for the reduced two component models and compare their profiles in homogeneous media. For a spatially random velocity field, we introduce travel time and study statistics of degradation fronts via representations in terms of the travel time probability density function (pdf) and the traveling front profiles. The travel time pdf does not vary with the nutrient and pollutant concentrations and only depends on the random water velocity. The traveling front profiles are expressed analytically or semi-analytically as functions of the travel time. The problem of nonlinear transport by a random velocity reduces to two subproblems: one being nonlinear transport by a known (unit) velocity, and the other being linear (advective) transport by a random velocity. The approach is illustrated through some examples where the randomness in velocity stems from the spatial variability of porosity. 相似文献
15.
Alternative fractional models of contaminant transport lead to a new travel time formula for arbitrary concentration levels. For an evolving contaminant plume in a highly heterogeneous aquifer, the new formula predicts much earlier arrival at low concentrations. Travel times of contaminant fronts and plumes are often obtained from Darcy's law calculations using estimates of average pore velocities. These estimates only provide information about the travel time of the average concentration (or peak, for contaminant pulses). Recently, it has been shown that finding the travel times of arbitrary concentration levels is a straightforward process, and equations were developed for other portions of the breakthrough curve for a nonreactive contaminant. In this paper, we generalize those equations to include alternative fractional models of contaminant transport. 相似文献
16.
The spread of a passive contaminant in an open-channel reach is considered with use of a two-dimensional advection-diffusion
equation with the included off-diagonal dispersion coefficients. This paper presents the calculation of truncation errors,
namely numerical diffusion and numerical dispersion for various finite difference schemes. The accuracy of the considered
finite-difference approximations is analysed by deriving and studying the relevant modified partial differential equation. 相似文献
17.
Time and Magnitude of Peak Concentration of Reactive Groundwater Contaminants Discharged to a River 下载免费PDF全文
下载免费PDF全文 An analytical solution for calculating the contaminant discharge rate in an aquifer following an instantaneous release of reactive contaminant mass to groundwater is used to derive relationships for the time and magnitude of peak concentration in a river receiving the transported material. Relationships are developed for the time of peak concentration relative to the time of travel for the contaminant, and the magnitude of peak concentration relative to the concentration calculated at the time of travel. Both quantities are found to be a function of two dimensionless parameters characterizing advective‐dispersive‐reactive transport—the Peclet number and the Damkohler number. It is shown that the time to peak concentration may occur before the time of travel, considering advection and retardation only, depending on the magnitudes of the Peclet and Damkohler numbers. Similarly, the magnitude of peak concentration may exceed the concentration calculated assuming that the time of peak concentration coincides with time of travel for the contaminant. For large Damkohler numbers, equating the time of peak concentration with the time of travel for the contaminant can significantly underestimate peak concentrations. 相似文献
18.
A function for the bed-load sediment transport rate is derived. This function is obtained by using the entrainment probabilities of the rolling and lifted sediment grains, and by introducing two travel lengths, respectively. The predictions from the new bed-load function agree well with experimental results over the entire experimental range and show significant improvement over the commonly used formula for the bed-load transport rate. The new function shows that, in terms of contributing to the bed-load transport rate, the total entrainment probability of the sediment grains is a weighted summation of those for the lifted and rolling grains, rather than a simple addition of the two. The function is also used to predict the total entrainment probability, saltation length, and the bed layer thickness at a high bed-load transport rate. These predictions all agree well with the experimental results. It is found that, on average, the travel length for the rolling sand grains is about an order of magnitude less than that of the lifted grains. 相似文献
19.
George A. Keramidas 《Advances in water resources》1981,4(4):165-173
The development of a displacement finite element formulation and its application to convective transport problems is presented. The formulation is based on the introduction of a generalized quantity defined as transport displacement. The governing equation is expressed in terms of this quantity and by using generalized coordinates a variational form of the governing equation is obtained. This equation may be solved by any numerical method, though it is of particular interest for application of the finite element method. Two finite element models are derived for the solution of convection-diffusion boundary value problems. The performance of the two element models is discussed and numerical results are given for different cases of convection and diffusion with two types of boundary conditions. The numerical results obtained show not only the efficiency of the numerical models in handling pure convection, pure diffusion and mixed convection-diffusion problems, but also good stability and accuracy. The applications of the developed numerical models are not limited to diffusion-convection problems but can also be applied to other types of problems such as mass transfer, hydrodynamics and wave propagation. 相似文献
20.
Simultaneous identification of the pollutant release history and the source location in groundwater by means of a geostatistical approach 总被引：4，自引：1，他引：3
Ilaria Butera Maria Giovanna Tanda Andrea Zanini 《Stochastic Environmental Research and Risk Assessment (SERRA)》2013,27(5):1269-1280
This paper describes an innovative procedure that is able to simultaneously identify the release history and the source location of a pollutant injection in a groundwater aquifer (simultaneous release function and source location identification, SRSI). The methodology follows a geostatistical approach: it develops starting from a data set and a reliable numerical flow and transport model of the aquifer. Observations can be concentration data detected at a given time in multiple locations or a time series of concentration measurements collected at multiple locations. The methodology requires a preliminary delineation of a probably source area and results in the identification of both the sub-area where the pollutant injection has most likely originated, and in the contaminant release history. Some weak hypotheses have to be defined about the statistical structure of the unknown release function such as the probability density function and correlation structure. Three case studies are discussed concerning two-dimensional, confined aquifers with strongly non-uniform flow fields. A transfer function approach has been adopted for the numerical definition of the sensitivity matrix and the recent step input function procedure has been successfully applied. 相似文献