共查询到20条相似文献，搜索用时 812 毫秒
1.
Uncertainty analysis of radioactive nuclide transport for one-dimensional single fracture has been studied. First order differential analysis is applied to introduce analytical form of output expectation and variance for contaminant transport equation, by regarding uncertainty of dispersion coefficient and retardation factor. Breakthrough curve of dimensionless concentration is demonstrated by taking I-129 as radioactive nuclide in fracture transport. It is possible to pick up critical ranges in spatial and temporal domain from the output variance. From the viewpoint of preliminary performance assessment for nuclear waste disposal the parameter importance in such system can be substantially measured in the site characterization in future. 相似文献
2.
The problem of one-dimensional transport of passive solute by a random steady velocity field is investigated. This problem is representative of solute movement in porous media, for example, in vertical flow through a horizontally stratified formation of variable porosity with a constant flux at the soil surface. Relating moments of particle travel time and displacement, exact expressions for the advection and dispersion coefficients in the Focker-Planck equation are compared with the perturbation results for large distances. The first- and second-order approximations for the dispersion coefficient are robust for a lognormal velocity field. The mean Lagrangian velocity is the harmonic mean of the Eulerian velocity for large distances. This is an artifact of one-dimensional flow where the continuity equation provides for a divergence free fluid flux, rather than a divergence free fluid velocity. 相似文献
3.
V. D. Cvetkovic G. Dagan A. M. Shapiro 《Stochastic Environmental Research and Risk Assessment (SERRA)》1991,5(1):45-54
The problem of one-dimensional transport of passive solute by a random steady velocity field is investigated. This problem is representative of solute movement in porous media, for example, in vertical flow through a horizontally stratified formation of variable porosity with a constant flux at the soil surface. Relating moments of particle travel time and displacement, exact expressions for the advection and dispersion coefficients in the Focker-Planck equation are compared with the perturbation results for large distances. The first- and second-order approximations for the dispersion coefficient are robust for a lognormal velocity field. The mean Lagrangian velocity is the harmonic mean of the Eulerian velocity for large distances. This is an artifact of one-dimensional flow where the continuity equation provides for a divergence free fluid flux, rather than a divergence free fluid velocity. 相似文献
4.
《Advances in water resources》2001,24(6):607-619
Solute discharge moments (mean and variance) are computed using numerical modeling of flow and advective transport in two-dimensional heterogeneous aquifers and are compared to theoretical results. The solute discharge quantifies the temporal evolution of the total contaminant mass crossing a certain compliance boundary. In addition to analyzing the solute discharge moments within a classical absolute dispersion framework, we also analyze relative dispersion formulation, whereby plume meandering (deviation from mean flow path caused by velocity variations at scales larger than plume size) is removed. This study addresses some important issues related to the computation of solute discharge moments from random walk particle tracking experiments, and highlights some of the important differences between absolute and relative dispersion frameworks. Relative dispersion formulation produces maximum uncertainty that coincides with the peak mean discharge. Absolute dispersion, however, results in earlier arrival of the uncertainty peak as compared to the first moment peak. Simulations show that the standard deviation of solute discharge in a relative dispersion framework requires increasingly large temporal sampling windows to smooth out some of the large fluctuations in breakthrough curves associated with advective transport. Using smoothing techniques in particle tracking to distribute the particle mass over a volume rather than at a point significantly reduces the noise in the numerical simulations and removes the need to use large temporal windows. Same effect can be obtained by adding a local dispersion process to the particle tracking experiments used to model advective transport. The effect of the temporal sampling window bears some relevance and important consequences for evaluating risk-related parameters. The expected value of peak solute discharge and its standard deviation are very sensitive to this sampling window and so will be the risk distribution relying on such numerical models. 相似文献
5.
Atrazine, an herbicide widely used for selective control of grassy weeds in the fields where maize is grown, is a potential ground water contaminant in China and, consequently, there is interest in predicting its mobility in agricultural soils. In this study, we determined the nonequilibrium transport characteristics of atrazine in Shahe soil (Beijing sandy loam) using the advection-dispersion equation, and using a sensitivity analysis, we evaluated the contribution of the uncertainty in a given input parameter to the overall uncertainty in model results. The asymmetrical shape and tailing of the atrazine breakthrough curve (BTC) showed that atrazine was subject to nonequilibrium transport. The observed atrazine BTC was best fitted by the chemical nonequilibrium model with a nonlinear least-squares optimization approach. Results from the sensitivity analysis indicated that the retardation factor was the most sensitive parameter. Considering the reliability of the estimated parameters, the best fit to the atrazine BTC was obtained by fixing the retardation factor based on the linear distribution coefficient, and by calculating the dispersion coefficient from the bromide BTC and the average pore water velocity from the measured data; nonequilibrium parameters were the only unknown parameters that were optimized. Model verification procedures were based on best-fit parameters optimized from one soil column experiment and applied to simulate the transport of atrazine in the duplicate experiment. The results showed there was good agreement between measured and simulated concentrations for atrazine leaching in the soil column. 相似文献
6.
《Advances in water resources》1998,21(6):487-498
Transport of a sorbing solute in a two-dimensional steady and uniform flow field is modeled using a particle tracking random walk method. The solute is initially introduced from an instantaneous point source. Cases of linear and nonlinear sorption isotherms are considered. Local pore velocity and mechanical dispersion are used to describe the solute transport mechanisms at the local scale. The numerical simulation of solute particle transport yields the large scale behavior of the solute plume. Behavior of the plume is quantified in terms of the center-of-mass displacement distance, relative velocity of the center-of-mass, mass breakthrough curves, spread variance, and longitudinal skewness. The nonlinear sorption isotherm affects the plume behavior in the following way relative to the linear isotherm: (1) the plume velocity decreases exponentially with time; (2) the longitudinal variance increases nonlinearly with time; (3) the solute front is steepened and tailing is enhanced 相似文献
7.
《Advances in water resources》2002,25(5):565-576
In the dispersion theory, a linear relationship has been verified between the coefficient of hydrodynamic dispersion and water velocity, both in saturated and in unsaturated porous media. But for unsaturated soils the variability of flow directions and microscopic velocities can be larger than in saturated soils because of the lower degree of water saturation. This leads to an increased dispersion. Therefore, relationships between water content and relative water velocity fluctuations and water content together with the coefficient of dispersivity in unsaturated porous media respectively have been investigated systematically by displacement experiments in glass beads and coarse-textured sandy soil columns. The breakthrough curves (BTCs) of chloride showed that an increase of solute mixing with a decrease of water content was caused by an increase of flow velocity fluctuations for different pathways. In order to explain the observed tailing effect in unsaturated flow, two mathematical models were used to fit theoretically derived nonlinear functions of water content dependent dispersivities for both porous media. The close agreement between the observed and computed results suggests that the theoretical model of hydrodynamic dispersion can be extended to transport in unsaturated porous media, providing that BTCs of the effluent water are used to estimate representative dispersivity parameters of soils. 相似文献
8.
An empirical hyperbolic scale-dependent dispersion model, which predicts a linear growth of dispersivity close to the origin and the attainment of an asymptotic dispersivity at large distances, is presented for deterministic modelling of field-scale solute transport and the analysis of solute transport experiments. A simple relationship is derived between local dispersivity, which is used in numerical simulations of solute transport, and effective dispersivity, which is estimated from the analysis of tracer breakthrough curves. The scale-dependent dispersion model is used to interpret a field tracer experiment by nonlinear least-squares inversion of a numerical solution for unsaturated transport. Simultaneous inversion of concentration-time data from several sampling locations indicates a linear growth of the dispersion process over the scale of the experiment. These findings are consistent with the results of an earlier analysis based on the use of a constant dispersion coefficient model at each of the sampling depths. 相似文献
9.
HUANG Guanhua HUANG Quanzhong ZHAN Hongbin CHEN Jing XIONG Yunwu FENG Shaoyuan 《中国科学D辑(英文版)》2005,48(Z2)
The newly developed Fractional Advection-Dispersion Equation (FADE), which is FADE was extended and used in this paper for modelling adsorbing contaminant transport by adding an adsorbing term. A parameter estimation method and its corresponding FORTRAN based program named FADEMain were developed on the basis of Nonlinear Least Square Algorithm and the analytical solution for one-dimensional FADE under the conditions of step input and steady state flow. Data sets of adsorbing contaminants Cd and NH4+-N transport in short homogeneous soil columns and conservative solute NaCI transport in a long homogeneous soil column, respectively were used to estimate the transport parameters both by FADEMain and the advection-dispersion equation (ADE) based program CXTFIT2.1. Results indicated that the concentration simulated by FADE agreed well with the measured data. Compared to the ADE model, FADE can provide better simulation for the concentration in the initial lower concentration part and the late higher concentration part of the breakthrough curves for both adsorbing contaminants. The dispersion coefficients for ADE were from 0.13 to 7.06 cm2/min, while the dispersion coefficients for FADE ranged from 0.119 to 3.05 cm1.856/min for NaCI transport in the long homogeneous soil column. We found that the dispersion coefficient of FADE increased with the transport distance, and the relationship between them can be quantified with an exponential function. Less scale-dependent was also found for the dispersion coefficient of FADE with respect to ADE. 相似文献
10.
To study contaminant transport in groundwater, an essential requirement is robust and accurate estimation of the transport parameters such as dispersion coefficient. The commonly used inverse error function method (IEFM) may cause unacceptable errors in dispersion coefficient estimation using the breakthrough curves (BTCs) data. We prove that the random error in the measured concentrations, which might be described by a normal distribution, would no longer follow the normal distribution after the IEFM transformation. In this study, we proposed a new method using the weighted least squares method (WLSM) to estimate the dispersion coefficient and velocity of groundwater. The weights were calculated based on the slope of the observed BTCs. We tested the new method against other methods such as genetic algorithm and CXTFIT program and found great agreement. This new method acknowledged different characteristics of solute transport at early, intermediate, and late time stages and divided BTCs into three sections for analysis. The developed method was applied to interpret three column tracer experiments by introducing continuous, constant‐concentration of sodium chloride (NaCl) into columns filled with sand, gravel, and sand‐gravel media. This study showed that IEFM performed well only when the observed data points were located in the linear (intermediate time) section of BTCs; it performed poorly when data points were in the early and late time stages. The new WLSM method, however, performed well for data points scattering over the entire BTCs and appeared promising in parameter estimation for solute transport in a column. 相似文献
11.
The groundwater remediation field has been changing constantly since it first emerged in the 1970s. The remediation field has evolved from a dissolved‐phase centric conceptual model to a DNAPL‐dominated one, which is now being questioned due to a renewed appreciation of matrix diffusion effects on remediation. Detailed observations about contaminant transport have emerged from the remediation field, and challenge the validity of one of the mainstays of the groundwater solute transport modeling world: the concept of mechanical dispersion (Payne et al. 2008). We review and discuss how a new conceptual model of contaminant transport based on diffusion (the usurper) may topple the well‐established position of mechanical dispersion (the status quo) that is commonly used in almost every groundwater contaminant transport model, and evaluate the status of existing models and modeling studies that were conducted using advection‐dispersion models. 相似文献
12.
Recognizing the spatial heterogeneity of hydraulic parameters, many researchers have studied the solute transport by both groundwater and channel flow in a stochastic framework. One of the methodologies used to up‐scale the stochastic solute transport equation, from a point‐location scale to a grid scale, is the cumulant expansion method combined with the calculus for the time‐ordered exponential and the calculus for the Lie operator. When the point‐location scale transport equation is scaled up to the grid scale, using the cumulant expansion method, a new dispersion coefficient emerges in the dispersive term of the solute transport equation in addition to the molecular dispersion coefficient. This velocity driven dispersion is called ‘macrodispersion’. The macrodispersion coefficient is the integral function of the time‐ordered covariance of the random velocity field. The integral is calculated over a Lagrangian trajectory of the flow. The Lagrangian trajectory depends on the following: (i) the spatial origin of the particle; (ii) the time when the macrodispersion is calculated; and (iii) the mean velocity field along the trajectory itself. The Lagrangian trajectory is a recursive function of time because the location of the particle along the trajectory at a particular time depends on the location of the particle at the previous time. This recursive functional form of the Lagrangian trajectory makes the calculation of the macrodispersion coefficient difficult. Especially for the unsteady, spatially non‐stationary, non‐uniform flow field, the macrodispersion coefficient is a highly complex expression and, so far, calculated using numerical methods in the discrete domains. Here, an analytical method was introduced to calculate the macrodispersion coefficient in the discrete domain for the unsteady and steady, spatially non‐stationary flow cases accurately and efficiently. This study can fill the gap between the theory of the ensemble averaged solute transport model and its numerical implementations. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
13.
Multi-offset phase analysis of surface wave data (MOPA) 总被引：1，自引：0，他引：1
Inaccuracy in the shear wave velocity profile inverted from surface wave data manifests from both modelling error and data uncertainty. An alternative method for dispersion curve evaluation by weighted linear regression of phase-offset data can be applied to both equispaced and non-equispaced data for objective identification of these often overlooked error sources.From field data, near-field effects are noted to at most half a wavelength and lateral discontinuities identified by marked changes in wavenumber with offset. Transition frequencies to dominant higher modes appear lower than when identified from standard plane-wave transform methods. Effects can be discriminated by their frequency, position or offset dependence.When a non-corrupt dispersion curve is extracted, the errors are up to 5% at low frequency. Through theoretical Gaussian error propagation analysis, the resulting shear wave velocity profile shows up to 18% uncertainty at depth. 相似文献
14.
We analyze the impact of conditioning to measurements of hydraulic transmissivity on the transport of a conservative solute. The effects of conditioning on solute transport are widely discussed in the literature, but most of the published works focuses on the reduction of the uncertainty in the prediction of the plume dispersion. In this study both ensemble and effective plume moments are considered for an instantaneous release of a solute through a linear source normal to the mean flow direction, by taking into account different sizes of the source. The analysis, involving a steady and spatially inhomogeneous velocity field, is developed by using the stochastic finite element method. Results show that conditioning reduces the ensemble moment in comparison with the unconditioned case, whereas the effective dispersion may increase because of the contribution of the spatial moments related to the lack of stationarity in the flow field. As the number of conditioning points increases, this effect increases and it is significant in both the longitudinal and transverse directions. Furthermore, we conclude that the moment derived from data collected in the field can be assessed by the conditioned second-order spatial moment only with a dense grid of measured data, and it is manifest for larger initial lengths of the plume. Nevertheless, it seems very likely that the actual dispersion of the plume may be underestimated in practical applications. 相似文献
15.
A macroscopic transport model is developed, following the Taylor shear dispersion analysis procedure, for a 2D laminar shear flow between parallel plates possessing a constant specified concentration. This idealized geometry models flow with contaminant dissolution at pore-scale in a contaminant source zone and flow in a rock fracture with dissolving walls. We upscale a macroscopic transient transport model with effective transport coefficients of mean velocity, macroscopic dispersion, and first-order mass transfer rate. To validate the macroscopic model the mean concentration, covariance, and wall concentration gradient are compared to the results of numerical simulations of the advection–diffusion equation and the Graetz solution. Results indicate that in the presence of local-scale variations and constant concentration boundaries, the upscaled mean velocity and macrodispersion coefficient differ from those of the Taylor–Aris dispersion, and the mass transfer flux described by the first-order mass transfer model is larger than the diffusive mass flux from the constant wall. In addition, the upscaled first-order mass transfer coefficient in the macroscopic model depends only on the plate gap and diffusion coefficient. Therefore, the upscaled first-order mass transfer coefficient is independent of the mean velocity and travel distance, leading to a constant pore-scale Sherwood number of 12. By contrast, the effective Sherwood number determined by the diffusive mass flux is a function of the Peclet number for small Peclet number, and approaches a constant of 10.3 for large Peclet number. 相似文献
16.
To more accurately predict the migration behavior of pollutants in porous media, we conduct laboratory scale experiments and model simulation. Aniline (AN) is used in one-dimensional soil column experiments designed under various media and hydrodynamic conditions. The advection-dispersion equation (ADE) and the continuous-time random walk (CTRW) were used to simulate the breakthrough curves (BTCs) of the solute transport. The results show that the media and hydrodynamic conditions are two important factors affecting solute transport and are related to the degree of non-Fickian transport. The simulation results show that CTRW can more effectively describe the non-Fickian phenomenon in the solute transport process than ADE. The sensitive parameter in the CTRW simulation process is , which can reflect the degree of non-Fickian diffusion in the solute transport. Understanding the relationship of with velocity and media particle size is conducive to improving the reactive solute transport model. The results of this study provide a theoretical basis for better prediction of pollutant transport in groundwater. 相似文献
17.
We present an analytical expression for the shear dispersion during solute transport in a coupled fracture–matrix system. The dispersion coefficient is obtained in a fracture with porous walls by taking into account an accurate boundary condition at the interface between the matrix and fracture, and the results were compared with those in a non-coupled system. The analysis presented identifies three regimes: diffusion-dominated, transition, and advection-dominated. The results showed that it is important to consider the exchange of solute between the fracture and matrix in development of the shear dispersion coefficient for the transition and advection-dominated regimes. The new dispersion coefficient is obtained by imposing the continuity of concentrations and mass fluxes along the porous walls. The resulting equivalent transport equation revealed that the effective velocity in a fracture increases while the dispersion coefficient decreases due to mass transfer between the matrix and fracture. A larger effective advection term leads to greater storage of mass in the matrix as compared with the classical double-porosity model with a non-coupled dispersion coefficient. The findings of this study can be used for modeling of tracer tests as well as fate, transport, and remediation of groundwater contaminants in fractured rocks. 相似文献
18.
《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. 相似文献
19.
Yong Zhang Dongbao Zhou Maosheng Yin HongGuang Sun Wei Wei Shiyin Li Chunmiao Zheng 《水文研究》2020,34(25):5104-5122
Modelling pollutant transport in water is one of the core tasks of computational hydrology, and various physical models including especially the widely used nonlocal transport models have been developed and applied in the last three decades. No studies, however, have been conducted to systematically assess the applicability, limitations and improvement of these nonlocal transport models. To fill this knowledge gap, this study reviewed, tested and improved the state-of-the-art nonlocal transport models, including their physical background, mathematical formula and especially the capability to quantify conservative tracers moving in one-dimensional sand columns, which represents perhaps the simplest real-world application. Applications showed that, surprisingly, neither the popular time-nonlocal transport models (including the multi-rate mass transfer model, the continuous time random walk framework and the time fractional advection-dispersion equation), nor the spatiotemporally nonlocal transport model (ST-fADE) can accurately fit passive tracers moving through a 15-m-long heterogeneous sand column documented in literature, if a constant dispersion coefficient or dispersivity is used. This is because pollutant transport in heterogeneous media can be scale-dependent (represented by a dispersion coefficient or dispersivity increasing with spatiotemporal scales), non-Fickian (where plume variance increases nonlinearly in time) and/or pre-asymptotic (with transition between non-Fickian and Fickian transport). These different properties cannot be simultaneously and accurately modelled by any of the transport models reviewed by this study. To bypass this limitation, five possible corrections were proposed, and two of them were tested successfully, including a time fractional and space Hausdorff fractal model which minimizes the scale-dependency of the dispersion coefficient in the non-Euclidean space, and a two-region time fractional advection-dispersion equation which accounts for the spatial mixing of solute particles from different mobile domains. Therefore, more efforts are still needed to accurately model transport in non-ideal porous media, and the five model corrections proposed by this study may shed light on these indispensable modelling efforts. 相似文献